Structural Basis of Human CYP51 Inhibition by Antifungal Azoles

The obligatory step in sterol biosynthesis in eukaryotes is demethylation of sterol precursors at the C14-position, which is catalyzed by CYP51 (sterol 14-alpha demethylase) in three sequential reactions. In mammals, the final product of the pathway is cholesterol, while important intermediates, meiosis-activating sterols, are produced by CYP51. Three crystal structures of human CYP51, ligand-free and complexed with antifungal drugs ketoconazole and econazole, were determined, allowing analysis of the molecular basis for functional conservation within the CYP51 family. Azole binding occurs mostly through hydrophobic interactions with conservative residues of the active site. The substantial conformational changes in the B′ helix and F-G loop regions are induced upon ligand binding, consistent with the membrane nature of the protein and its substrate. The access channel is typical for mammalian sterol-metabolizing P450 enzymes, but is different from that observed in Mycobacterium tuberculosis CYP51. Comparison of the azole-bound structures provides insight into the relative binding affinities of human and bacterial P450 enzymes to ketoconazole and fluconazole, which can be useful for the rational design of antifungal compounds and specific modulators of human CYP51.     

Structural diversities of active site in clinical azole-bound forms between sterol 14alpha-demethylases (CYP51s) from human and Mycobacterium tuberculosis

To gain insights into the molecular basis of the design for the selective azole anti-fungals, we compared the binding properties of azole-based inhibitors for cytochrome P450 sterol 14alpha-demethylase (CYP51) from human (HuCYP51) and Mycobacterium tuberculosis (MtCYP51). Spectroscopic titration of azoles to the CYP51s revealed that HuCYP51 has higher affinity for ketoconazole (KET), an azole derivative that has long lipophilic groups, than MtCYP51, but the affinity for fluconazole (FLU), which is a member of the anti-fungal armamentarium, was lower in HuCYP51. The affinity for 4-phenylimidazole (4-PhIm) to MtCYP51 was quite low compared with that to HuCYP51. In the resonance Raman spectra for HuCYP51, the FLU binding induced only minor spectral changes, whereas the prominent high frequency shift of the bending mode of the heme vinyl group was detected in the KET- or 4-PhIm-bound forms. On the other hand, the bending mode of the heme propionate group for the FLU-bound form of MtCYP51 was shifted to high frequency as found for the KET-bound form, but that for 4-PhIm was shifted to low frequency. The EPR spectra for 4-PhIm-bound MtCYP51 and FLU-bound HuCYP51 gave multiple g values, showing heterogeneous binding of the azoles, whereas the single gx and gz values were observed for other azole-bound forms. Together with the alignment of the amino acid sequence, these spectroscopic differences suggest that the region between the B' and C helices, particularly the hydrophobicity of the C helix, in CYP51s plays primary roles in determining strength of interactions with azoles; this differentiates the binding specificity of azoles to CYP51s.     

Fluconazole binding and sterol demethylation in three CYP51 isoforms indicate differences in active site topology

14alpha-Demethylase (CYP51) is a key enzyme in all sterol biosynthetic pathways (animals, fungi, plants, protists, and some bacteria), catalyzing the removal of the C-14 methyl group following cyclization of squalene. Based on mutations found in CYP51 genes from Candida albicans azole-resistant isolates obtained after fluconazole treatment of fungal infections, and using site-directed mutagenesis, we have found that fluconazole binding and substrate metabolism vary among three different CYP51 isoforms: human, fungal, and mycobacterial. In C. albicans, the Y132H mutant from isolates shows no effect on fluconazole binding, whereas the F145L mutant results in a 5-fold increase in its IC(50) for fluconazole, suggesting that F145 (conserved only in fungal 14alpha-demethylases) interacts with this azole. In C. albicans, F145L accounts, in part, for the difference in fluconazole sensitivity reported between mammals and fungi, providing a basis for treatment of fungal infections. The C. albicans Y132H and human Y145H CYP51 mutants show essentially no effect on substrate metabolism, but the Mycobacterium tuberculosis F89H CYP51 mutant loses both its substrate binding and metabolism. Because these three residues align in the three isoforms, the results indicate that their active sites contain important structural differences, and further emphasize that fluconazole and substrate binding are uncoupled properties.     

CYP51 from Trypanosoma brucei is obtusifoliol-specific

New isoforms of CYP51 (sterol 14alpha-demethylase), an essential enzyme in sterol biosynthesis and primary target of azole antimycotic drugs, are found in pathogenic protists, Trypanosoma brucei(TB), T. vivax, T. cruzi, and Leishmania major. The sequences share approximately 80% amino acid identity and are approximately 25% identical to sterol 14alpha-demethylases from other biological kingdoms. Differences of residues conserved throughout the rest of the CYP51 family that align with the BC-loop and helices F and G of CYP51 from Mycobacterium tuberculosis (MT)) imply possible alterations in the topology of the active site cavity of the protozoan enzymes. CYP51 and cytochrome P450 reductase (CPR) from TB were cloned, expressed in Escherichia coli, and purified. The P450 has normal spectral features (including absolute absorbance, carbon monoxide, and ligand binding spectra), is efficiently reduced by TB and rat CPR but demonstrates altered specificity in comparison with human CYP51 toward three tested azole inhibitors, and contrary to the human, Candida albicans, and MT isoforms, reveals profound substrate preference toward obtusifoliol (turnover 5.6 min(-1)). It weakly interacts with the other known CYP51 substrates; slow lanosterol conversion predominantly produces the 14alpha-carboxyaldehyde intermediate. Although obtusifoliol specificity is typical for plant isoforms of CYP51, the set of sterol biosynthetic enzymes in the protozoan genomes together with available information about sterol composition of kinetoplastid cells suggest that the substrate preference of TBCYP51 may reflect a novel sterol biosynthetic pathway in Trypanosomatidae.     

Functional expression and characterization of CYP51 from dandruff-causing Malassezia globosa

Malassezia globosa is one of the most common yeasts to cause various human skin diseases including dandruff and seborrheic dermatitis. Genomic analysis of M. globosa revealed four putative cytochrome P450 (CYP) enzymes. Here, we report the purification and characterization of recombinant CYP51, a putative lanosterol 14α-demethylase, from M. globosa. The M. globosa CYP51 was expressed heterologously in Escherichia coli, followed by purification. Purified CYP51 showed a typical reduced CO-difference spectrum of P450, with a maximum absorption at 447?nm. Purified CYP51 exhibited tight binding to azole antifungal agents such as ketoconazole, econazole, fluconazole, or itraconazole, with K(d) values around 0.26-0.84?μM, which suggests that CYP51 is an orthologous target for antifungal agents in the M. globosa. In addition, three mutations (Y127F, A169S, and K176N) in the amino acid sequence of M. globosa CYP51 were identified in one of the azole-resistant strains. Homology modeling of M. globosa CYP51 suggested that the Y127F mutation may influence the resistance to azoles by blocking substrate access channels. Taken together, functional expression and characterization of the CYP51 enzyme can provide a fundamental basis for a specific antifungal drug design for dandruff caused by M. globosa.     

Molecular modeling of human lanosterol 14α-demethylase complexes with substrates and their derivatives

Lanosterol 14α-demethylase (CYP51A1) is a key enzyme in sterol biosynthesis. In humans, this enzyme is involved in the cholesterol biosynthesis pathway. The majority of antifungal drugs are aimed at the inhibition of CYP51 in fungi. To elucidate the molecular mechanisms of highly specific protein-ligand recognition, we have developed a full-atomic model of human CYP51A1 and performed docking of natural substrates and their derivatives to the active site of the enzyme. The parameters of the binding enthalpy of substrates, intermediates, and final products of the reaction of 14α-demethylation were estimated using the MMPB(GB)SA algorithm. Dynamic properties and conformational changes of the protein globule upon binding of the ligand near the active site have been investigated by the molecular dynamics method. Our studies reveal that hydroxylated intermediate reaction products have a greater affinity than the initial substrates, which facilitates the multistage reaction without accumulation of intermediate products. The contribution to the free energy of steroid ligand binding of 30 amino acids forming the substrate-binding region of CYP51A1, as well as the influence of their substitutions to alanine on the stability of the protein molecule, has been clarified using alanine scanning modeling. We demonstrate that the most serious weakening of the binding is observed in the case of substitutions Y137A, F145A, V149A, I383A, and R388A. The results of molecular modeling are in agreement with the data obtained by analysis of primary sequences of representatives of the CYP51 family.     

Molecular dynamic modeling of CYP51B in complex with azole inhibitors

Cytochrome P450 14α-sterol demethylase (CYP51), the key enzyme in sterol biosynthesis pathway, is an important target protein of cholesterol-lowering agents, antifungal drugs, and herbicides. CYP51B enzyme is one of the CYP51 family members. In the present study, we have chosen four representative inhibitors of CYP51B: diniconazole (Din), fluconazole (Flu), tebuconazole (Teb), and voriconazole (Vor), and launched to investigate the binding features of CYP51B-inhibitor and gating characteristics via molecular docking and molecular dynamics (MD) simulations. The results show that the ranking of binding affinities among these inhibitors is Din > Teb > Vor > Flu. Din shows the strongest binding affinity, whereas Flu shows the weakest binding affinity. More importantly, based on the calculated binding free energy results, we hypothesize that the nonpolar interactions are the most important contributors, and three key residues (Thr77, Ala258, and Lys454) play crucial role in protein-inhibitor binding. Besides, residue Phe180 may play a molecular switch role in the process of the CYP51B-Teb and CYP51B-Vor binding. Additionally, Tunnel analysis results show that the major tunnel of Din, Flu, and Teb is located between helix K, FG loop, and β4 hairpin (Tunnel II).The top ranked possible tunnel (Tunnel II) corresponds to Vor exits through helix K, F and helix J. This study further revealed the CYP51B relevant structural characteristics at the atomic level as well as provided a basis for rational design of new and more efficacious antifungal agents.     

CYP51 from Trypanosoma cruzi: a phyla-specific residue in the B' helix defines substrate preferences of sterol 14alpha-demethylase

A potential drug target for treatment of Chagas disease, sterol 14alpha-demethylase from Trypanosoma cruzi (TCCYP51), was found to be catalytically closely related to animal/fungi-like CYP51. Contrary to the ortholog from Trypanosoma brucei (TB), which like plant CYP51 requires C4-monomethylated sterol substrates, TCCYP51 prefers C4-dimethylsterols. Sixty-six CYP51 sequences are known from bacteria to human, their sequence homology ranging from approximately 25% between phyla to approximately 80% within a phylum. TC versus TB is the first example of two organisms from the same phylum, in which CYP51s (83% amino acid identity) have such profound differences in substrate specificity. Substitution of animal/fungi-like Ile105 in the B' helix to Phe, the residue found in this position in all plant and the other six CYP51 sequences from Trypanosomatidae, dramatically alters substrate preferences of TCCYP51, converting it into a more plant-like enzyme. The rates of 14alpha-demethylation of obtusifoliol and its 24-demethyl analog 4alpha-,4alpha-dimethylcholesta-8,24-dien-3beta-ol(norlanosterol) increase 60- and 150-fold, respectively. Turnover of the three 4,4-dimethylated sterol substrates is reduced approximately 3.5-fold. These catalytic properties correlate with the sterol binding parameters, suggesting that Phe in this position provides necessary interactions with C4-monomethylated substrates, which Ile cannot. The CYP51 substrate preferences imply differences in the post-squalene portion of sterol biosynthesis in TC and TB. The phyla-specific residue can be used to predict preferred substrates of new CYP51 sequences and subsequently for the development of new artificial substrate analogs, which might serve as highly specific inhibitors able to kill human parasites.     

Estriol bound and ligand-free structures of sterol 14alpha-demethylase

Sterol 14alpha-demethylases (CYP51) are essential enzymes in sterol biosynthesis in eukaryotes and drug targets in antifungal therapy. Here, we report CYP51 structures in ligand-free and estriol bound forms. Using estriol as a probe, we determined orientation of the substrate in the active site, elucidated protein contacts with the invariant 3beta-hydroxy group of a sterol, and identified F78 as a key discriminator between 4alpha-methylated and 4alpha,beta-dimethylated substrates. Analysis of CYP51 dynamics revealed that the C helix undergoes helix-coil transition upon binding and dissociation of a ligand. Loss of helical structure of the C helix in the ligand-free form results in an unprecedented opening of the substrate binding site. Upon binding of estriol, the BC loop loses contacts with molecular surface and tends to adopt a closed conformation. A mechanism for azole resistance in the yeast pathogen Candida albicans associated with mutations in the ERG11 gene encoding CYP51 is suggested based on CYP51 protein dynamics.     

Structural complex of sterol 14α-demethylase (CYP51) with 14α-methylenecyclopropyl-Delta7-24, 25-dihydrolanosterol

Sterol 14α-demethylase (CYP51) that catalyzes the removal of the 14α-methyl group from the sterol nucleus is an essential enzyme in sterol biosynthesis, a primary target for clinical and agricultural antifungal azoles and an emerging target for antitrypanosomal chemotherapy. Here, we present the crystal structure of Trypanosoma (T) brucei CYP51 in complex with the substrate analog 14α-methylenecyclopropyl-Δ7-24,25-dihydrolanosterol (MCP). This sterol binds tightly to all protozoan CYP51s and acts as a competitive inhibitor of F105-containing (plant-like) T. brucei and Leishmania (L) infantum orthologs, but it has a much stronger, mechanism-based inhibitory effect on I105-containing (animal/fungi-like) T. cruzi CYP51. Depicting substrate orientation in the conserved CYP51 binding cavity, the complex specifies the roles of the contact amino acid residues and sheds new light on CYP51 substrate specificity. It also provides an explanation for the effect of MCP on T. cruzi CYP51. Comparison with the ligand-free and azole-bound structures supports the notion of structural rigidity as the characteristic feature of the CYP51 substrate binding cavity, confirming the enzyme as an excellent candidate for structure-directed design of new drugs, including mechanism-based substrate analog inhibitors.     

X-ray structure of 4,4'-dihydroxybenzophenone mimicking sterol substrate in the active site of sterol 14alpha-demethylase (CYP51)

A universal step in the biosynthesis of membrane sterols and steroid hormones is the oxidative removal of the 14alpha-methyl group from sterol precursors by sterol 14alpha-demethylase (CYP51). This enzyme is a primary target in treatment of fungal infections in organisms ranging from humans to plants, and development of more potent and selective CYP51 inhibitors is an important biological objective. Our continuing interest in structural aspects of substrate and inhibitor recognition in CYP51 led us to determine (to a resolution of 1.95A) the structure of CYP51 from Mycobacterium tuberculosis (CYP51(Mt)) co-crystallized with 4,4'-dihydroxybenzophenone (DHBP), a small organic molecule previously identified among top type I binding hits in a library screened against CYP51(Mt). The newly determined CYP51(Mt)-DHBP structure is the most complete to date and is an improved template for three-dimensional modeling of CYP51 enzymes from fungal and prokaryotic pathogens. The structure demonstrates the induction of conformational fit of the flexible protein regions and the interactions of conserved Phe-89 essential for both fungal drug resistance and catalytic function, which were obscure in the previously characterized CYP51(Mt)-estriol complex. DHBP represents a benzophenone scaffold binding in the CYP51 active site via a type I mechanism, suggesting (i) a possible new class of CYP51 inhibitors targeting flexible regions, (ii) an alternative catalytic function for bacterial CYP51 enzymes, and (iii) a potential for hydroxybenzophenones, widely distributed in the environment, to interfere with sterol biosynthesis. Finally, we show the inhibition of M. tuberculosis growth by DHBP in a mouse macrophage model.     

Differential inhibition of human CYP3A4 and Candida albicans CYP51 with azole antifungal agents

The inhibition by azole antifungals of human cytochrome CYP3A4, the major form of drug metabolising enzyme within the liver, was compared with their inhibitory activity against their target enzyme, Candida albicans sterol 14alpha-demethylase (CYP51), following heterologous expression in Saccharomyces cerevisiae. IC(50) values for ketoconazole and itraconazole CYP3A4 inhibition were 0.25 and 0. 2 microM. These values compared with much lower doses required for the complete inhibition of C. albicans CYP51, where IC(50) values of 0.008 and 0.0076 microM were observed for ketoconazole and itraconazole, respectively. Additionally, stereoselective inhibition of CYP3A4 and CYP51 was observed with enantiomers of the azole antifungal compounds diclobutrazol and SCH39304. In both instances, the RR(+) configuration at their asymmetric carbon centres was most active. Interestingly, the SS(-) enantiomeric form of SCH39304 was inactive and failed to bind CYP3A4, as demonstrable by Type II binding spectra.     

Rapid P450 heme iron reduction by laser photoexcitation of Mycobacterium tuberculosis CYP121 and CYP51B1. Analysis of CO complexation reactions and reversibility of the P450/P420 equilibrium

We demonstrate that photoexcitation of NAD(P)H reduces heme iron of Mycobacterium tuberculosis P450s CYP121 and CYP51B1 on the microsecond time scale. Rates of formation for the ferrous-carbonmonoxy (Fe(II)-CO) complex were determined across a range of coenzyme/CO concentrations. CYP121 reaction transients were biphasic. A hyperbolic dependence on CO concentration was observed, consistent with the presence of a CO binding site in ferric CYP121. CYP51B1 absorption transients for Fe(II)-CO complex formation were monophasic. The reaction rate was second order with respect to [CO], suggesting the absence of a CO-binding site in ferric CYP51B1. In the absence of CO, heme iron reduction by photoexcited NAD(P)H is fast ( approximately 10,000-11,000 s(-1)) with both P450s. For CYP121, transients revealed initial production of the thiolate-coordinated (P450) complex (absorbance maximum at 448 nm), followed by a slower phase reporting partial conversion to the thiol-coordinated P420 species (at 420 nm). The slow phase amplitude increased at lower pH values, consistent with heme cysteinate protonation underlying the transition. Thus, CO binding occurs to the thiolate-coordinated ferrous form prior to cysteinate protonation. For CYP51B1, slow conversions of both the ferrous/Fe(II)-CO forms to species with spectral maxima at 423/421.5 nm occurred following photoexcitation in the absence/presence of CO. This reflected conversion from ferrous thiolate- to thiol-coordinated forms in both cases, indicating instability of the thiolate-coordinated ferrous CYP51B1. CYP121 Fe(II)-CO complex pH titrations revealed reversible spectral transitions between P450 and P420 forms. Our data provide strong evidence for P420 formation linked to reversible heme thiolate protonation, and demonstrate key differences in heme chemistry and CO binding for CYP121 and CYP51B1.     

Homology modeling,molecular docking and spectra assay studies of sterol 14α-demethylase from <Emphasis Type="Italic">Penicillium digitatum</Emphasis>

Sterol 14α-demethylase from Penicillium digitatum (PdCYP51) is a prime target of antifungal drugs for citrus disease in plants. To design novel antifungal compounds, a homology model of PdCYP51 was constructed using the recently reported crystal structure of human CYP51 as the template. Molecular docking was performed to investigate the interaction of four commercial fungicides with the modeled enzyme. The side chain of these compounds interplayed with PdCYP51 mainly through hydrophobic and van der Waals interactions. Biochemical spectra analysis of inhibitors combined with PdCYP51 are also compatible with the docking results. This is the first molecular modeling for PdCYP51 based on the eukaryotic crystal structure of CYP51. The structural information and binding site mapping of PdCYP51 for different inhibitors obtained from this study could aid in screening and designing new antifungal compounds targeting this enzyme.     

Biophysical characterization of the sterol demethylase P450 from Mycobacterium tuberculosis, its cognate ferredoxin, and their interactions

Mycobacterium tuberculosis encodes a P450 of the sterol demethylase family (CYP51) chromosomally located adjacent to a ferredoxin (Fdx). CYP51 and Fdx were purified to homogeneity and characterized. Spectroscopic analyses were consistent with cysteinate- and aqua-ligated heme iron in CYP51. An epsilon419 of 134 mM(-1) cm(-1) was determined for oxidized CYP51. Analysis of interactions of 1-, 2-, and 4-phenylimidazoles with CYP51 showed that the 1- and 4-forms were heme iron-coordinating inhibitors, while 2-phenylimidazole induced a substrate-like optical shift. The 2-phenyimidazole-bound CYP51 demonstrated unusual decreases in high-spin heme iron content at elevated temperatures and an almost complete absence of high-spin heme iron by low-temperature EPR. These data suggest thermally induced alterations in CYP51 active site structure and/or binding modes for the small ligand. Reduction of CYP51 in the presence of carbon monoxide leads to formation of an Fe(II)-CO complex with a Soret absorption maximum at 448.5 nm, which collapses (at 0.246 min(-1) at pH 7.0) forming a species with a Soret maximum at 421.5 nm (the inactive P420 form). The rate of P420 formation is accelerated at lower pH, consistent with protonation of the cysteinate (Cys 394) to a thiol underlying the P450-P420 transition. The P450 form is stabilized by estriol, which induces a type I spectral shift on binding CYP51 (Kd = 21.7 microM). Nonstandard spectral changes occur on CYP51 reduction (using either dithionite or natural redox partners), including a blue-shifted Soret band and development of a strong feature at approximately 558.5 nm, suggestive of cysteine thiol ligation. Thus, ligand-free ferrous CYP51 is prone to thiolate ligand protonation even in the absence of carbon monoxide. Analysis of reoxidized CYP51 demonstrates that the enzyme re-forms P450, indicating that Cys 394 thiol is readily deprotonated to thiolate in the ferric form. Spectroscopic analysis of Fdx by EPR (resonance at g = 2.03) and magnetic CD (intensity for oxidized and reduced forms and signal intensity dependence on field strength and temperature) demonstrated that Fdx binds a [3Fe-4S] iron-sulfur cluster. Potentiometric studies show that the midpoint potential for ligand-free CYP51 is -375 mV, increasing to -225 mV in the estriol-bound form. The Fdx potential is -31 mV. Fdx forms a productive electron transfer complex with CYP51 and reduces it at a rate of 3.0 min(-1) in the ligand-free form and 4.3 min(-1) in the estriol-bound form, despite a thermodynamic barrier. Steady-state analysis of a M. tuberculosis class I redox system comprising flavoprotein reductase A (FprA), Fdx, and estriol-bound CYP51 indicates heme iron reduction as a rate-limiting step.     

Dynamics of CYP51: implications for function and inhibitor design

Sterol 14α‐demethylase (cytochrome P450 family 51 (CYP51)) is an essential enzyme occurring in all biological kingdoms. In eukaryotes, it is located in the membrane of the endoplasmic reticulum. Selective inhibitors of trypanosomal CYP51s that do not affect the human CYP51 have been discovered in vitro and found to cure acute and chronic mouse Chagas disease without severe side effects in vivo. Crystal structures indicate that CYP51 may be more rigid than most CYPs, and it has been proposed that this property may facilitate antiparasitic drug design. Therefore, to investigate the dynamics of trypanosomal CYP51, we built a model of membrane‐bound Trypanosoma brucei CYP51 and then performed molecular dynamics simulations of T. brucei CYP51 in membrane‐bound and soluble forms. We compared the dynamics of T. brucei CYP51 with those of human CYP51, CYP2C9, and CYP2E1. In the simulations, the CYP51s display low mobility in the buried active site although overall mobility is similar in all the CYPs studied. The simulations suggest that in CYP51, pathway 2f serves as the major ligand access tunnel, and both pathways 2f (leading to membrane) and S (leading to solvent) can serve as ligand egress tunnels. Compared with the other CYPs, the residues at the entrance of the ligand access tunnels in CYP51 have higher mobility that may be necessary to facilitate the passage of its large sterol ligands. The water (W) tunnel is accessible to solvent during most of the simulations of CYP51, but its width is affected by the conformations of the heme's two propionate groups. These differ from those observed in the other CYPs studied because of differences in their hydrogen‐bonding network. Our simulations give insights into the dynamics of CYP51 that complement the available experimental data and have implications for drug design against CYP51 enzymes. Copyright © 2015 John Wiley & Sons, Ltd.