Estimating a transfer function with a known subset of poles/zeros

Hello, I want to estimate a transfer function from measured time domain data. My special case is, that I know exactly the place of a Pole and a Zero of the System. (High Pass filter included in the measurement device)

So I'm looking for the System without the High Pass.

Right now I'm estimating the system with 3 Zeros and 1 Pole (with tfest) and afterwards I cancel the closest Pole and Zero from the HP.

But as the canceled estimated poles are not equal to the HP poles, there might be a better estimation when Matlab knows in advance which poles are included in the system. Is there a way to specify given poles/zeros in the estimation progress?