Extending the Sobolev theory of quasiconformal and quasiregular maps to subdomains of the complex plane motivates our investigation of Sobolev regularity of singular integral operators on domains. We introduce new paraproducts which lead to higher order T1-type testing conditions. A special case provides weighted Sobolev estimates for the compressed Beurling transform which imply quantitative Sobolev estimates for the Beltrami resolvent. This is joint work with Francesco Di Plinio and Brett D. Wick.