Computational simulations uncover enantioselective metabolism of chiral triazole fungicides by human CYP450 enzymes: A case study of tebuconazole

Tebuconazole (TEB), a prominent chiral triazole fungicide, has been extensively utilized for plant pathogen control globally. Despite experimental evidence of TEB metabolism in mammals, the enantioselectivity in the biotransformation of R- and S-TEB enantiomers by specific CYP450s remains elusive. In this work, integrated in silico simulations were employed to unveil the binding interactions and enantioselective metabolic fate of TEB enantiomers within human CYP1A2, 2B6, 2E1, and 3A4. Molecular dynamics (MD) simulations clearly delineated the binding specificity of R- and S-TEB to the four CYP450s, crucially determining their differences in metabolic activity and enantioselectivity. The primary driving force for robust ligand binding was identified as van der Waals interactions with CYP450s, particularly involving the hydrophobic residues. Mechanistic insights derived from quantum mechanics/molecular mechanics (QM/MM) calculations established C₂-methyl hydroxylation as the predominant route of R-/S-TEB metabolism, while C₆-hydroxylation and triazol epoxidation were deemed kinetically infeasible pathways. Specifically, the resulting hydroxy-R-TEB metabolite primarily originates from R-TEB biotransformation by 1A2, 2E1 and 3A4, whereas hydroxy-S-TEB is preferentially produced by 2B6. These findings significantly contribute to our comprehension of the binding specificity and enantioselective metabolic fate of chiral TEB by CYP450s, potentially informing further research on human health risk assessment associated with TEB exposure.