You will need an invertible realization. Try high order transfer function estimations with no delays. Or, if the estimated model with delays is a good/accurate model, compute its frequency response (freqresp), invert it, and fit a transfer function to the inverse frequency response.
For literally just removing the delay you can do model.ioDelay= 0