Simulations of CYP51A from Aspergillus fumigatus in a model bilayer provide insights into triazole drug resistance

Azole antifungal drugs target CYP51A in Aspergillus fumigatus by binding with the active site of the protein, blocking ergosterol biosynthesis. Resistance to azole antifungal drugs is now common, with a leucine to histidine amino acid substitution at position 98 the most frequent, predominantly conferring resistance to itraconazole, although cross-resistance has been reported in conjunction with other mutations. In this study, we create a homology model of CYP51A using a recently published crystal structure of the paralog protein CYP51B. The derived structures, wild type, and L98H mutant are positioned within a lipid membrane bilayer and subjected to molecular dynamics simulations in order improve the accuracy of both models. The structural analysis from our simulations suggests a decrease in active site surface from the formation of hydrogen bonds between the histidine substitution and neighboring polar side chains, potentially preventing the binding of azole drugs. This study yields a biologically relevant structure and set of dynamics of the A. fumigatus Lanosterol 14 alpha-demethylase enzyme and provides further insight into azole antifungal drug resistance.