A binding mode hypothesis for prothioconazole binding to CYP51 derived from first principles quantum chemistry

In order to assess safety and efficacy of small molecule drugs as well as agrochemicals, it is key to understanding the nature of protein-ligand interaction on an atomistic level. Prothioconazole (PTZ), although commonly considered to be an azole-like inhibitor of sterol 14-α demethylase (CYP51), differs from classical azoles with respect to how it binds its target. The available evidence is only indirect, as crystallographic elucidation of CYP51 complexed with PTZ have not yet been successful. We derive a binding mode hypothesis for PTZ binding its target, compare to DPZ, a triazole-type metabolite of PTZ, and set our findings into context of its biochemistry and spectroscopy. Quantum Theory of Atoms in Molecules (QTAIM) analysis of computed DFT electron densities is used to qualitatively understand the topology of binding, revealing significant differences of how R- and S-enantiomers are binding and, in particular, how the thiozolinthione head of PTZ binds to heme compared to DPZ's triazole head. The difference of binding enthalpy is calculated at coupled cluster (DLPNO-CCSD(T)) level of theory, and we find that DPZ binds stronger to CYP51 than PTZ by more than ΔH ~ 11 kcal/mol.