Conservation in the CYP51 family. Role of the B' helix/BC loop and helices F and G in enzymatic function

CYP51 (sterol 14 alpha-demethylase) is an essential enzyme in sterol biosynthetic pathways and the only P450 gene family having catalytically identical orthologues in different biological kingdoms. The proteins have low sequence similarity across phyla, and the whole family contains about 40 completely conserved amino acid residues. Fifteen of these residues lie in the secondary structural elements predicted to form potential substrate recognition sites within the P450 structural fold. The role of 10 of these residues, in the B' helix/BC loop, helices F and G, has been studied by site-directed mutagenesis using as a template the soluble sterol 14 alpha-demethylase of known structure, CYP51 from Mycobacterium tuberculosis (MT) and the human orthologue. Single amino acid substitutions of seven residues (Y76, F83, G84, D90, L172, G175, and R194) result in loss of the ability of the mutant MTCYP51 to metabolize lanosterol. Residual activity of D195A is very low, V87A is not expressed as a P450, and A197G has almost 1 order of magnitude increased activity. After purification, all of the mutants show normal spectral properties, heme incorporation, and the ability to be reduced enzymatically and to interact with azole inhibitors. Profound influence on the catalytic activity correlates well with the spectral response to substrate binding, effect of substrate stabilization on the reduced state of the P450, and substrate-enhanced efficiency of enzymatic reduction. Mutagenesis of corresponding residues in human CYP51 implies that the conserved amino acids might be essential for the evolutionary conservation of sterol 14 alpha-demethylation from bacteria to mammals.     

Conformational dynamics in the F/G segment of CYP51 from Mycobacterium tuberculosis monitored by FRET

A cysteine was introduced into the FG-loop (P187C) of CYP51 from Mycobacterium tuberculosis (MT) for selective labeling with BODIPY and fluorescence energy transfer (FRET) analysis. F?rster radius for the BODIPY-heme pair was calculated assuming that the distance between the heme and Cys187 in solution corresponds to that in the crystal structure of ligand free MTCYP51. Interaction of MTCYP51 with azole inhibitors ketoconazole and fluconazole or the substrate analog estriol did not influence the fluorescence, but titration with the substrate lanosterol quenched BODIPY emission, the effect being proportional to the portion of substrate bound MTCYP51. The detected changes correspond to approximately 10A decrease in the calculated distance between BODIPY-Cys187 and the heme. The results confirm (1) functional importance of conformational motions in the MTCYP51 F/G segment and (2) applicability of FRET to monitor them in solution.     

Substrate recognition sites in cytochrome P450 family 2 (CYP2) proteins inferred from comparative analyses of amino acid and coding nucleotide sequences.

The substrate recognition regions in cytochrome P450 family 2 (CYP2) proteins were inferred by group-to-group alignment of CYP2 sequences and those of bacterial P450s, including Pseudomonas putida P450 101A (P450cam), whose substrate-binding residues have been definitely identified by x-ray crystallography of a substrate-bound form (Poulos T. L., Finzel, B. C., and Howard, A. J. (1987) J. Mol. Biol. 195, 687-700). The six putative substrate recognition sites, SRSs, thus identified are dispersively located along the primary structure and constitute about 16% of the total residues. All the reported point mutations and chimeric fragments that significantly affect the substrate specificities of the parental CYP2 enzymes fell within or overlapped some of the six SRSs. Analysis of nucleotide substitution patterns in closely related members in four subfamilies, CYP2A, 2B, 2C, and 2D, consistently indicated that the SRSs have accumulated more nonsynonymous (amino acid-changing) substitutions than the rest of the sequence. This observation supports the idea that diversification of duplicate genes of drug-metabolizing P450s occurs primarily in substrate recognition regions to cope with an increasing number of foreign compounds.     

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.     

Structure-function relationships of rat liver CYP3A9 to its human liver orthologs: site-directed active site mutagenesis to a progesterone dihydroxylase

CYP3A9 is an estrogen-inducible ortholog of human liver CYP3A4 with 76.5% sequence identity to CYP3A4. Unlike CYP3A4, it is a very poor testosterone 6beta- and 2beta-hydroxylase, but a relatively better catalyst of progesterone monohydroxylation largely at 6beta, 16alpha, and 21 positions with negligible 6beta, 21-dihydroxylation. We reasoned that such differences in substrate catalyses must be due to differences in the active site architecture of each CYP3A enzyme. Indeed, alignment of CYP3A4 substrate recognition sites (SRSs) with the corresponding regions of CYP3A9 sequence revealed that of the 22 fully divergent residues, 4 reside in SRS regions [P107N (SRS-1), M371G (SRS-5), and L479K and G480Q (SRS-6)]. Accordingly, we substituted these and other divergent CYP3A9 SRS residues with the corresponding residues of CYP3A4 and/or CYP3A5. Our findings of the influence of these site-directed mutations of the CYP3A9 active site on its catalysis of testosterone and three other established but structurally different CYP3A substrates (progesterone, imipramine, and carbamazepine) are described. These findings revealed that some mutations (N107P, N107S, V207T, G371M, and Q480G) not only improved the ability of CYP3A9 to hydroxylate testosterone at the 6beta and 2beta positions, but also converted it into a robust progesterone 6beta, 21-dihydroxylase. The latter in the case of CYP3A9N107P was accompanied by a shift from sigmoidal to hyperbolic enzyme-substrate kinetics. In contrast, the catalytic potential of CYP3A9 mutants K206N, K206S, M240V, and K479L/Q480G was either relatively unchanged or negligible to nonexistent. Together these findings attest to the unique substrate-active site fit of each CYP3A enzyme.     

Sterol 14-demethylase P450 (CYP51) provides a breakthrough for the discussion on the evolution of cytochrome P450 gene superfamily

Biodiversity is the most characteristic feature of cytochrome P450. Finding of CYP51 distributing widely in biological kingdoms provided breakthroughs for the discussion on the evolution and diversification of P450. Molecular phylogenetic analysis demonstrated that CYP51 appeared in the prokaryotic era and distributed into most kingdoms concomitant with phylogenetic divergence. This is the first evolutionary evidence indicating the prokaryotic origin of P450. Modification of substrate specificity of eukaryotic CYP51s occurred independently to adapt to the different sterol precursors existing in each kingdom. Formation of CYP51 variants through the mutation of active site and the selection of the advantageous ones from them were demonstrated by the emergence of azole-resistant CYP51s in Candida albicans under the environments rich in azole antifungal agents. These findings illustrate the most probable core process of P450 diversification consisting of modification of active site and selection of the resulting variants through interaction with endogenous and exogenous chemicals.     

Impact of incorporating the 2C5 crystal structure into comparative models of cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) metabolizes approximately one third of the drugs in current clinical use. To gain insight into its structure and function, we have produced four different sets of comparative models of 2D6: one based on the structures of P450s from four different microorganisms (P450 terp, P450 eryF, P450 cam, and P450 BM3), another on the only mammalian P450 (2C5) structure available, and the other two based on alternative amino acid sequence alignments of 2D6 with all five of these structures. Principal component analysis suggests that inclusion of the 2C5 crystal structure has a profound effect on the modeling process, altering the general topology of the active site, and that the models produced differ significantly from all of the templates. The four models of 2D6 were also used in conjunction with molecular docking to produce complexes with the substrates codeine and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP); this identified Glu 216 [in the F-helix; substrate recognition site (SRS) 2] as a key determinant in the binding of the basic moiety of the substrate. Our studies suggest that both Asp 301 and Glu 216 are required for metabolism of basic substrates. Furthermore, they suggest that Asp 301 (I-helix, SRS-4), a residue thought from mutagenesis studies to bind directly to the basic moiety of substrates, may play a key role in positioning the B'-C loop (SRS-1) and that the loss of activity on mutating Asp 301 may therefore be the result of an indirect effect (movement of the B'-C loop) on replacing this residue.     

Essential requirements for substrate binding affinity and selectivity toward human CYP2 family enzymes

A detailed analysis of substrate selectivity within the cytochrome P450 2 (CYP2) family is reported. From a consideration of specific interactions between drug substrates for human CYP2 family enzymes and the putative active sites of CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, and CYP2E1, it is likely that the number and disposition of hydrogen bond donor/acceptors and aromatic rings within the various P450 substrate molecules determines their enzyme selectivity and binding affinity, together with directing their preferred routes of metabolism by the CYP2 enzymes concerned. Although many aliphatic residues are present in most P450 active sites, it would appear that their main contribution centers around hydrophobic interactions and desolvation processes accompanying substrate binding. Molecular modeling studies based on the recent CYP2C5 crystal structure appear to show close agreement with site-directed mutagenesis experiments and with information on substrate metabolism and selectivity within the CYP2 family.     

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.     

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 evidence for a functionally relevant second camphor binding site in P450cam: model for substrate entry into a P450 active site

P450cam has long served as a prototype for the cytochrome P450 (CYP) gene family. But, little is known about how substrate enters its active site pocket, and how access is achieved in a way that minimizes exposure of the reactive heme. We hypothesize that P450cam may first bind substrate transiently near the mobile F-G helix that covers the active site pocket. Such a two-step binding process is kinetically required if P450cam rarely populates an open conformation-as suggested by previous literature and the inability to obtain a crystal structure of P450cam in an open conformation. Such a mechanism would minimize exposure of the heme by allowing P450cam to stay in a closed conformation as long as possible, since only brief flexing into an open conformation would be required to allow substrate entry. To test this model, we have attempted to dock a second camphor molecule into the crystal structure of camphor-bound P450cam. The docking identified only one potential entry site pocket, a well-defined cavity on the F-helix side of the F-G flap, 16 A from the heme iron. Location of this entry site pocket is consistent with our NMR T1 relaxation-based measurements of distances for a camphor that binds in fast exchange (active site camphor is known to bind in slow exchange). Presence of a second camphor binding site is also confirmed with [(1)H-(13)C] HSQC titrations of (13)CH3-threonine labeled P450cam. To confirm that camphor can bind outside of the active site pocket, (13)CH3-S-pyridine was bound to the heme iron to physically block the active site, and to serve as an NMR chemical shift probe. Titration of this P450cam-pyridine complex confirms that camphor can bind to a site outside the active site pocket, with an estimated Kd of 43 microM. The two-site binding model that is proposed based on these data is analogous to that recently proposed for CYP3A4, and is consistent with recent crystal structures of P450cam bound to tethered-substrates, which force a partially opened conformation.     

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.     

Evidence for communality in the primary determinants of CYP74 catalysis and of structural similarities between CYP74 and classical mammalian P450 enzymes

In silico structural analysis of CYP74C3, a membrane-associated P450 enzyme from the plant Medicago truncatula (barrel medic) with hydroperoxide lyase (HPL) specificity, showed that it had strong similarities to the structural folds of the classical microsomal P450 enzyme from rabbits (CYP2C5). It was not only the secondary structure predictions that supported the analysis but site directed mutagenesis of the substrate interacting residues was also consistent with it. This led us to develop a substrate-binding model of CYP74C3 which predicted three amino acid residues, N285, F287, and G288 located in the putative I-helix and distal haem pocket of CYP74C3 to be in close proximity to the preferred substrate 13-HPOTE. These residues were judged to be in equivalent positions to those identified in SRS-4 of CYP2C5. Significance of the residues and their relevance to the model were further assessed by site directed mutagenesis of the three residues followed by EPR spectroscopic and detailed kinetic investigations of the mutated proteins in the presence and absence of detergent. Although point mutation of the residues had no effect on the haem content of the mutated proteins, significant effects on the spin state equilibrium of the haem iron were noted. Kinetic effects of the mutations, which were investigated using three different substrates, were dramatic in nature. In the presence of detergent with the preferred substrate (13-HPOTE), the catalytic center activities and substrate binding affinities of the mutant proteins were reduced by a factor of 8-32 and 4-12, respectively, compared with wild-type--a two orders of magnitude reduction in catalytic efficiencies. We believe this is the first report where primary determinants of catalysis for any CYP74 enzyme, which are fully consistent with our model, have been identified. Our working model predicts that N285 is close enough to suggest that a hydrogen bond with the peroxy group of the enzyme substrate 13-HPOTE is warranted, whereas significance of F287 may arise from a strong hydrophobic interaction between the alkyl group(s) of the substrate and the phenyl ring of F287. We believe that G288 is crucial because of its size. Any other residue with a relatively bulky side chain will hinder the access of substrate to the active site. The effects of the mutations suggests that subtle protein conformational changes in the putative substrate-binding pocket regulate the formation of a fully active monomer-micelle complex with low-spin haem iron and that structural communication exists between the substrate- and micelle-binding sites of CYP74C3. Conservation in CYP74 sequence alignments suggests that N285, F287, and G288 in CYP74C3 and the equivalent residues at positions in other CYP74 enzymes are likely to be critical to catalysis. To support this we show that G324 in CYP74D4 (Arabidopsis AOS), equivalent to G288 in CYP74C3, is a primary determinant of positional specificity. We suggest that the overall structure of CYP74 enzymes is likely to be very similar to those described for classical P450 monooxygenase enzymes.     

Structure-function analysis of cytochromes P450 2B

In the last 4 years, breakthroughs were made in the field of P450 2B (CYP2B) structure-function through determination of one ligand-free and two inhibitor-bound X-ray crystal structures of CYP2B4, which revealed many of the structural features required for binding ligands of different size and shape. Large conformational changes of several plastic regions of CYP2B4 can dramatically reshape the active site of the enzyme to fit the size and shape of the bound ligand without perturbing the overall P450 fold. Solution biophysical studies using isothermal titration calorimetry (ITC) have revealed the large difference in the thermodynamic parameters of CYP2B4 in binding inhibitors of different ring chemistry and side chains. Other studies have revealed that the effects of site-specific mutations on steady-state kinetic parameters and mechanism-based inactivation are often substrate dependent. These findings agree with the structural data that the enzymes adopt different conformations to bind various ligands. Thus, the substrate specificity of an individual enzyme is determined not only by active site residues but also non-active site residues that modulate conformational changes that are important for substrate access and rearrangement of the active site to accommodate the bound substrate.     

Insect cytochromes P450: topology of structural elements predicted to govern catalytic versatility

Among eukaryotic P450s, the greatest expansion has been in insects, providing useful model systems for the study of enzyme evolution in response to natural and anthropogenic pressures such as the chemical warfare against plant toxins and synthetic insecticides. To better understand diversification of the catalytic properties in the various P450 clades, insight into the molecular principles governing biotransformation of the array of endogenous and exogenous compounds is of paramount importance. Based on a general, CYP102A1-related construct, the majority of prospective substrate-docking residues were found to cluster near the distal heme face within the six known substrate recognition sites (SRSs) made up by the α-helical B′, F, G and I tetrad as well as the B′-C turn and strands of certain β-sheets. Highest proportion of contact sites, having a mean frequency of conservation of not more than 10%, was detected to be housed in the variable SRS-1 and SRS-2 domains, affording conformational mobility to meet the structural diversity in the bulk of substrates. Some key determinants may have a function in governing substrate access, positioning ligands toward the catalytic center to allow regioselective attack, and promoting dioxygen activation. Moreover, non-SRS elements occupying regions on the bowl-like proximal heme face may serve in recognition of electron delivery partners. Physico-chemical factors interfering with substrate anchoring mainly relate to the variably expressed lipophilicity/hydrophilicity profile of the diverse binding sites and bulkiness of the side chains of critical amino acids. Detailed knowledge about structure-function relationships in insect P450s, thus, may offer an incentive to tailoring of novel insecticides and synergistic enzyme inhibitors to more efficiently control insect pests.     

Crystal structure of the Mycobacterium tuberculosis P450 CYP121-fluconazole complex reveals new azole drug-P450 binding mode

Azole and triazole drugs are cytochrome P450 inhibitors widely used as fungal antibiotics and possessing potent antimycobacterial activity. We present here the crystal structure of Mycobacterium tuberculosis cytochrome P450 CYP121 in complex with the triazole drug fluconazole, revealing a new azole heme ligation mode. In contrast to other structurally characterized cytochrome P450 azole complexes, where the azole nitrogen directly coordinates the heme iron, in CYP121 fluconazole does not displace the aqua sixth heme ligand but occupies a position that allows formation of a direct hydrogen bond to the aqua sixth heme ligand. Direct ligation of fluconazole to the heme iron is observed in a minority of CYP121 molecules, albeit with severe deviations from ideal geometry due to close contacts with active site residues. Analysis of both ligand-on and -off structures reveals the relative position of active site residues derived from the I-helix is a key determinant in the relative ratio of on and off states. Regardless, both ligand-bound states lead to P450 inactivation by active site occlusion. This previously unrecognized means of P450 inactivation is consistent with spectroscopic analyses in both solution and in the crystalline form and raises important questions relating to interaction of azoles with both pathogen and human P450s.     

A survey of active site access channels in cytochromes P450

In cytochrome P450s, the active site is situated deep inside the protein next to the heme cofactor, and is often completely isolated from the surrounding solvent. To identify routes by which substrates may enter into and products exit from the active site, random expulsion molecular dynamics simulations were performed for three cytochrome P450s: CYP101, CYP102A1 and CYP107A1 [J. Mol. Biol. 303 (2000) 797; Proc. Natl. Acad. Sci. USA 99 (2002) 5361]. Amongst the different pathways identified, one pathway was found to be common to all three cytochrome P450s although the mechanism of ligand passage along it was different in each case and apparently adapted to the substrate specificity of the enzyme. Recently, a number of new crystal structures of cytochrome P450s have been solved. Here, we analyse the open channels leading to the active site that these structures reveal. We find that in addition to showing the common pathway, they provide experimental evidence for the existence of three additional channels that were identified by simulation. We also discuss how the location of xenon binding sites in CYP101 suggests a role for one of the pathways identified by molecular dynamics simulations as a route for gaseous species, such as oxygen, to access the active site.     

Structure–function analysis of cytochromes P450 2B

In the last 4 years, breakthroughs were made in the field of P450 2B (CYP2B) structure–function through determination of one ligand-free and two inhibitor-bound X-ray crystal structures of CYP2B4, which revealed many of the structural features required for binding ligands of different size and shape. Large conformational changes of several plastic regions of CYP2B4 can dramatically reshape the active site of the enzyme to fit the size and shape of the bound ligand without perturbing the overall P450 fold. Solution biophysical studies using isothermal titration calorimetry (ITC) have revealed the large difference in the thermodynamic parameters of CYP2B4 in binding inhibitors of different ring chemistry and side chains. Other studies have revealed that the effects of site-specific mutations on steady-state kinetic parameters and mechanism-based inactivation are often substrate dependent. These findings agree with the structural data that the enzymes adopt different conformations to bind various ligands. Thus, the substrate specificity of an individual enzyme is determined not only by active site residues but also non-active site residues that modulate conformational changes that are important for substrate access and rearrangement of the active site to accommodate the bound substrate.     

In silico analysis of cytochrome P450 monooxygenases in chronic granulomatous infectious fungus Sporothrix schenckii: Special focus on CYP51

Sporotrichosis is an emerging chronic, granulomatous, subcutaneous, mycotic infection caused by Sporothrix species. Sporotrichosis is treated with the azole drug itraconazole as ketoconazole is ineffective. It is a well-known fact that azole drugs act by inhibiting cytochrome P450 monooxygenases (P450s), heme-thiolate proteins. To date, nothing is known about P450s in Sporothrix schenckii and the molecular basis of its resistance to ketoconazole. Here we present genome-wide identification, annotation, phylogenetic analysis and comprehensive P450 family-level comparative analysis of S. schenckii P450s with pathogenic fungi P450s, along with a rationale for ketoconazole resistance by S. schenckii based on in silico structural analysis of CYP51. Genome data-mining of S. schenckii revealed 40 P450s in its genome that can be grouped into 32 P450 families and 39 P450 subfamilies. Comprehensive comparative analysis of P450s revealed that S. schenckii shares 11 P450 families with plant pathogenic fungi and has three unique P450 families: CYP5077, CYP5386 and CYP5696 (novel family). Among P450s, CYP51, the main target of azole drugs was also found in S. schenckii. 3D modeling of S. schenckii CYP51 revealed the presence of characteristic P450 motifs with exceptionally large reductase interaction site 2. In silico analysis revealed number of mutations that can be associated with ketoconazole resistance, especially at the channel entrance to the active site. One of possible reason for better stabilization of itraconazole, compared to ketoconazole, is that the more extended molecule of itraconazole may form a hydrogen bond with ASN-230. This in turn may explain its effectiveness against S. schenckii vis-a-vis resistant to ketoconazole. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.     

Evaluation of structural features in fungal cytochromes P450 predicted to rule catalytic diversification

Fungi belong to the large kingdom of lower eukaryotic organisms encompassing yeasts along with filamentous and dimorphic members. Microbial P450 enzymes have contributed to exploration of and adaptation to diverse ecological niches such as conversion of lipophilic compounds to more hydrophilic derivatives or degradation of a vast array of environmental toxicants. To better understand diversification of the catalytic behavior of fungal P450s, detailed insight into the molecular machinery steering oxidative attack on the distinctly structured endogenous and xenobiotic substrates is of preeminent interest. Based on a general, CYP102A1-related template the bulk of predicted substrate/inhibitor-binding determinants were shown to cluster near the distal heme face within the six known substrate recognition sites (SRSs) made up by the α-helical B′/F/G/I tetrad, the B′–C interhelical loop and strands of the β6-sheet, population density being highest in the structurally flexible SRS-1 and SRS-4 domains, showing a low degree of conservation. Reactivity toward ligands favorably coincides with the lipophilicity/hydrophilicity profile and bulkiness of critical amino acids acting as selective filters. Some decisive elements may also serve in maintenance of catalytic competence via their action as gatekeepers directing substrate access/positioning or stabilizers of the heme environment enabling dioxygen activation. Non-SRS residues seem to control spin state equilibria and attract redox partners by electrostatic forces. Of note, the inhibitory potency of azole-type fungicides is likely to arise from perturbation of the complex interplay of the mechanistic principles addressed above. Knowledge-supported exploitation of the topological data will be helpful in the manufacture of commodity/specialty chemicals as well as therapeutic agents. Also, engineered fungal P450s may be used to improve pollutant-specific bioremediation of contaminated soils.