Mammalian microsomal cytochrome P450 monooxygenase: structural adaptations for membrane binding and functional diversity

Effects of membrane mimetics on cytochrome P450-cytochrome b5 interactions characterized by NMR spectroscopy

Mammalian cytochrome P450 (P450) is a membrane-bound monooxygenase whose catalytic activities require two electrons to be sequentially delivered from its redox partners: cytochrome b5 (cytb5) and cytochrome P450 reductase, both of which are membrane proteins. Although P450 functional activities are known to be affected by lipids, experimental evidence to reveal the effect of membrane on P450-cytb5 interactions is still lacking. Here, we present evidence for the influence of phospholipid bilayers on complex formation between rabbit P450 2B4 (CYP2B4) and rabbit cytb5 at the atomic level, utilizing NMR techniques. General line broadening and modest chemical shift perturbations of cytb5 resonances characterize CYP2B4-cytb5 interactions on the intermediate time scale. More significant intensity attenuation and a more specific protein-protein binding interface are observed in bicelles as compared with lipid-free solution, highlighting the importance of the lipid bilayer in stabilizing stronger and more specific interactions between CYP2B4 and cytb5, which may lead to a more efficient electron transfer. Similar results observed for the interactions between CYP2B4 lacking the transmembrane domain (tr-CYP2B4) and cytb5 imply interactions between tr-CYP2B4 and the membrane surface, which might assist in CYP2B4-cytb5 complex formation by orienting tr-CYP2B4 for efficient contact with cytb5. Furthermore, the observation of weak and nonspecific interactions between CYP2B4 and cytb5 in micelles suggests that lipid bilayer structures and low curvature membrane surface are preferable for CYP2B4-cytb5 complex formation. Results presented in this study provide structural insights into the mechanism behind the important role that the lipid bilayer plays in the interactions between P450s and their redox partners.

Interactions among cytochromes P450 in microsomal membranes: oligomerization of cytochromes P450 3A4, 3A5, and 2E1 and its functional consequences

The body of evidence of physiologically relevant P450-P450 interactions in microsomal membranes continues to grow. Here we probe oligomerization of human CYP3A4, CYP3A5, and CYP2E1 in microsomal membranes. Using a technique based on luminescence resonance energy transfer, we demonstrate that all three proteins are subject to a concentration-dependent equilibrium between the monomeric and oligomeric states. We also observed the formation of mixed oligomers in CYP3A4/CYP3A5, CYP3A4/CYP2E1, and CYP3A5/CYP2E1 pairs and demonstrated that the association of either CYP3A4 or CYP3A5 with CYP2E1 causes activation of the latter enzyme. Earlier we hypothesized that the intersubunit interface in CYP3A4 oligomers is similar to that observed in the crystallographic dimers of some microsomal drug-metabolizing cytochromes P450 (Davydov, D. R., Davydova, N. Y., Sineva, E. V., Kufareva, I., and Halpert, J. R. (2013) Pivotal role of P450-P450 interactions in CYP3A4 allostery: the case of α-naphthoflavone. Biochem. J. 453, 219-230). Here we report the results of intermolecular cross-linking of CYP3A4 oligomers with thiol-reactive bifunctional reagents as well as the luminescence resonance energy transfer measurements of interprobe distances in the oligomers of labeled CYP3A4 single-cysteine mutants. The results provide compelling support for the physiological relevance of the dimer-specific peripheral ligand-binding site observed in certain CYP3A4 structures. According to our interpretation, these results reveal an important general mechanism that regulates the activity and substrate specificity of the cytochrome P450 ensemble through interactions between multiple P450 species. As a result of P450-P450 cross-talk, the catalytic properties of the cytochrome P450 ensemble cannot be predicted by simple summation of the properties of the individual P450 species.

Mechanistic basis of electron transfer to cytochromes p450 by natural redox partners and artificial donor constructs

Cytochromes P450 (P450s) are hemoproteins catalyzing oxidative biotransformation of a vast array of natural and xenobiotic compounds. Reducing equivalents required for dioxygen cleavage and substrate hydroxylation originate from different redox partners including diflavin reductases, flavodoxins, ferredoxins and phthalate dioxygenase reductase (PDR)-type proteins. Accordingly, circumstantial analysis of structural and physicochemical features governing donor-acceptor recognition and electron transfer poses an intriguing challenge. Thus, conformational flexibility reflected by togging between closed and open states of solvent exposed patches on the redox components was shown to be instrumental to steered electron transmission. Here, the membrane-interactive tails of the P450 enzymes and donor proteins were recognized to be crucial to proper orientation toward each other of surface sites on the redox modules steering functional coupling. Also, mobile electron shuttling may come into play. While charge-pairing mechanisms are of primary importance in attraction and complexation of the redox partners, hydrophobic and van der Waals cohesion forces play a minor role in docking events. Due to catalytic plasticity of P450 enzymes, there is considerable promise in biotechnological applications. Here, deeper insight into the mechanistic basis of the redox machinery will permit optimization of redox processes via directed evolution and DNA shuffling. Thus, creation of hybrid systems by fusion of the modified heme domain of P450s with proteinaceous electron carriers helps obviate the tedious reconstitution procedure and induces novel activities. Also, P450-based amperometric biosensors may open new vistas in pharmaceutical and clinical implementation and environmental monitoring.

Computational Identification and Systematic Classification of Novel Cytochrome P450 Genes in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most economically important medicinal plants. Cytochrome P450 (CYP450) genes have been implicated in the biosynthesis of its active components. However, only a dozen full-length CYP450 genes have been described, and there is no systematic classification of CYP450 genes in S. miltiorrhiza. We obtained 77,549 unigenes from three tissue types of S. miltiorrhiza using RNA-Seq technology. Combining our data with previously identified CYP450 sequences and scanning with the CYP450 model from Pfam resulted in the identification of 116 full-length and 135 partial-length CYP450 genes. The 116 genes were classified into 9 clans and 38 families using standard criteria. The RNA-Seq results showed that 35 CYP450 genes were co-expressed with CYP76AH1, a marker gene for tanshinone biosynthesis, using r≥0.9 as a cutoff. The expression profiles for 16 of 19 randomly selected CYP450 obtained from RNA-Seq were validated by qRT-PCR. Comparing against the KEGG database, 10 CYP450 genes were found to be associated with diterpenoid biosynthesis. Considering all the evidence, 3 CYP450 genes were identified to be potentially involved in terpenoid biosynthesis. Moreover, we found that 15 CYP450 genes were possibly regulated by antisense transcripts (r≥0.9 or r≤–0.9). Lastly, a web resource (SMCYP450, http://www.herbalgenomics.org/samicyp450) was set up, which allows users to browse, search, retrieve and compare CYP450 genes and can serve as a centralized resource.

Targeting of splice variants of human cytochrome P450 2C8 (CYP2C8) to mitochondria and their role in arachidonic acid metabolism and respiratory dysfunction

In this study, we found that the full-length CYP2C8 (WT CYP2C8) and N-terminal truncated splice variant 3 (∼ 44-kDa mass) are localized in mitochondria in addition to the endoplasmic reticulum. Analysis of human livers showed that the mitochondrial levels of these two forms varied markedly. Molecular modeling based on the x-ray crystal structure coordinates of CYP2D6 and CYP2C8 showed that despite lacking the N-terminal 102 residues variant 3 possessed nearly complete substrate binding and heme binding pockets. Stable expression of cDNAs in HepG2 cells showed that the WT protein is mostly targeted to the endoplasmic reticulum and at low levels to mitochondria, whereas variant 3 is primarily targeted to mitochondria and at low levels to the endoplasmic reticulum. Enzyme reconstitution experiments showed that both microsomal and mitochondrial WT CYP2C8 efficiently catalyzed paclitaxel 6-hydroxylation. However, mitochondrial variant 3 was unable to catalyze this reaction possibly because of its inability to stabilize the large 854-Da substrate. Conversely, mitochondrial variant 3 catalyzed the metabolism of arachidonic acid into 8,9-, 11,12-, and 14,15-epoxyeicosatrienoic acids and 20-hydroxyeicosatetraenoic acid when reconstituted with adrenodoxin and adrenodoxin reductase. HepG2 cells stably expressing variant 3 generated higher levels of reactive oxygen species and showed a higher level of mitochondrial respiratory dysfunction. This study suggests that mitochondrially targeted variant 3 CYP2C8 may contribute to oxidative stress in various tissues.

Structural characterization of human cholesterol 7α-hydroxylase

Hepatic conversion to bile acids is a major elimination route for cholesterol in mammals. CYP7A1 catalyzes the first and rate-limiting step in classic bile acid biosynthesis, converting cholesterol to 7α-hydroxycholesterol. To identify the structural determinants that govern the stereospecific hydroxylation of cholesterol, we solved the crystal structure of CYP7A1 in the ligand-free state. The structure-based mutation T104L in the B' helix, corresponding to the nonpolar residue of CYP7B1, was used to obtain crystals of complexes with cholest-4-en-3-one and with cholesterol oxidation product 7-ketocholesterol (7KCh). The structures reveal a motif of residues that promote cholest-4-en-3-one binding parallel to the heme, thus positioning the C7 atom for hydroxylation. Additional regions of the binding cavity (most distant from the access channel) are involved to accommodate the elongated conformation of the aliphatic side chain. Structural complex with 7KCh shows an active site rigidity and provides an explanation for its inhibitory effect. Based on our previously published data, we proposed a model of cholesterol abstraction from the membrane by CYP7A1 for metabolism. CYP7A1 structural data provide a molecular basis for understanding of the diversity of 7α-hydroxylases, on the one hand, and cholesterol-metabolizing enzymes adapted for their specific activity, on the other hand.

Structural and energetic analysis to provide insight residues of CYP2C9, 2C11 and 2E1 involved in valproic acid dehydrogenation selectivity

Docking and molecular dynamics (MD) simulation have been two computational techniques used to gain insight about the substrate orientation within protein active sites, allowing to identify potential residues involved in the binding and catalytic mechanisms. In this study, both methods were combined to predict the regioselectivity in the binding mode of valproic acid (VPA) on three cytochrome P-450 (CYP) isoforms CYP2C9, CYP2C11, and CYP2E1, which are involved in the biotransformation of VPA yielding reactive hepatotoxic intermediate 2-n-propyl-4-pentenoic acid (4nVPA). There are experimental data about hydrogen atom abstraction of the C4-position of VPA to yield 4nVPA, however, there are not structural evidence about the binding mode of VPA and 4nVPA on CYPs. Therefore, the complexes between these CYP isoforms and VPA or 4nVPA were studied to explore their differences in binding and energetic stabilization. Docking results showed that VPA and 4nVPA are coupled into CYPs binding site in a similar conformation, but it does not explain the VPA hydrogen atom abstraction. On the other hand, MD simulations showed a set of energetic states that reorient VPA at the first ns, then making it susceptible to a dehydrogenation reaction. For 4nVPA, multiple binding modes were observed in which the different states could favor either undergo other reaction mechanism or ligand expulsion from the binding site. Otherwise, the energetic and entropic contribution point out a similar behavior for the three CYP complexes, showing as expected a more energetically favorable binding free energy for the complexes between CYPs and VPA than with 4nVPA.

Toward reduction in animal sacrifice for drugs: molecular modeling of Macaca fascicularis P450 2C20 for virtual screening of Homo sapiens P450 2C8 substrates

Macaca fascicularis P450 2C20 shares 92% identity with human cytochrome P450 2C8, which is involved in the metabolism of more than 8% of all prescribed drugs. To date, only paclitaxel and amodiaquine, two substrate markers of the human P450 2C8, have been experimentally confirmed as M. fascicularis P450 2C20 drugs. To bridge the lack of information on the ligands recognized by M. fascicularis P450 2C20, in this study, a three-dimensional homology model of this enzyme was generated on the basis of the available crystal structure of the human homologue P450 2C8 using YASARA. The results indicated that 90.0%, 9.0%, 0.5%, and 0.5% of the residues of the P450 2C20 model were located in the most favorable, allowed, generously allowed, and disallowed regions, respectively. The root-mean-square deviation of the C-alpha superposition of the M. fascicularis P450 2C20 model with the Homo sapiens P450 2C8 was 0.074 Å, indicating a very high similarity of the two structures. Subsequently, the 2C20 model was used for in silico screening of 58 known P450 2C8 substrates and 62 inhibitors. These were also docked in the active site of the crystal structure of the human P450 2C8. The affinity of each compound for the active site of both cytochromes proved to be very similar, meaning that the few key residues that are mutated in the active site of the M. fascicularis P450 do not prevent the P450 2C20 from recognizing the same substrates as the human P450 2C8.

Human cytochrome P450 17A1 conformational selection: modulation by ligand and cytochrome b5

Crystallographic studies of different membrane cytochrome P450 enzymes have provided examples of distinct structural conformations, suggesting protein flexibility. It has been speculated that conformational selection is an integral component of substrate recognition and access, but direct evidence of such substate interconversion has thus far remained elusive. In the current study, solution NMR revealed multiple and exchanging backbone conformations for certain structural features of the human steroidogenic cytochrome P450 17A1 (CYP17A1). This bifunctional enzyme is responsible for pregnenolone C17 hydroxylation, followed by a 17,20-lyase reaction to produce dehydroepiandrosterone, the key intermediate in human synthesis of androgen and estrogen sex steroids. The distribution of CYP17A1 conformational states was influenced by temperature, binding of these two substrates, and binding of the soluble domain of cytochrome b5 (b5). Notably, titration of b5 to CYP17A1·pregnenolone induced a set of conformational states closely resembling those of CYP17A1·17α-hydroxypregnenolone without b5, providing structural evidence consistent with the reported ability of b5 to selectively enhance 17,20-lyase activity. Solution NMR thus revealed a set of conformations likely to modulate human steroidogenesis by CYP17A1, demonstrating that this approach has the potential to make similar contributions to understanding the functions of other membrane P450 enzymes involved in drug metabolism and disease states.

Structure-based ligand design to overcome CYP inhibition in drug discovery projects

Cytochrome P450 (CYP) enzymes are key players in xenobiotic metabolism, and inhibition of CYPs can therefore result in unwanted drug-drug interactions. Within drug discovery, CYP inhibition can cause delays in the progression of candidate drugs, or even premature closure of projects. During the past decade, a massive effort in the pharmaceutical industry and academic research has produced a wealth of structural information in the CYP field. In this short review, we will describe how structure-based approaches can be used in the pharmaceutical industry to work away from CYP inhibition, with a focus on the opportunities and challenges. We will show two examples from our own work where structural information on CYP2C9 and CYP3A4 inhibitor complexes have been successfully exploited in ongoing drug discovery projects.

Genome-wide association study for feedlot average daily gain in Nellore cattle (Bos indicus)

The genome-wide association study (GWAS) results are presented for average daily gain (ADG) in Nellore cattle. Phenotype of 720 male Bos indicus animals with information of ADG in feedlots and 354,147 single-nucleotide polymorphisms (SNPs) obtained from a database added by information from Illumina Bovine HD (777,962 SNPs) and Illumina BovineSNP50 (54,609) by imputation were used. After quality control and imputation, 290,620 SNPs remained in the association analysis, using R package Genome-wide Rapid Association using Mixed Model and Regression method GRAMMAR-Gamma. A genomic region with six significant SNPs, at Bonferroni-corrected significance, was found on chromosome 3. The most significant SNP (rs42518459, BTA3: 85849977, p = 9.49 × 10(-8)) explained 5.62% of the phenotypic variance and had the allele substitution effect of -0.269 kg/day. Important genes such as PDE4B, LEPR, CYP2J2 and FGGY are located near this region, which is overlapped by 12 quantitative trait locus (QTLs) described for several production traits. Other regions with markers with suggestive effects were identified in BTA6 and BTA10. This study showed regions with major effects on ADG in Bos indicus in feedlots. This information may be useful to increase the efficiency of selecting this trait and to understand the physiological processes involved in its regulation.

Architecture of a single membrane spanning cytochrome P450 suggests constraints that orient the catalytic domain relative to a bilayer

Bitopic integral membrane proteins with a single transmembrane helix play diverse roles in catalysis, cell signaling, and morphogenesis. Complete monospanning protein structures are needed to show how interaction between the transmembrane helix and catalytic domain might influence association with the membrane and function. We report crystal structures of full-length Saccharomyces cerevisiae lanosterol 14α-demethylase, a membrane monospanning cytochrome P450 of the CYP51 family that catalyzes the first postcyclization step in ergosterol biosynthesis and is inhibited by triazole drugs. The structures reveal a well-ordered N-terminal amphipathic helix preceding a putative transmembrane helix that would constrain the catalytic domain orientation to lie partly in the lipid bilayer. The structures locate the substrate lanosterol, identify putative substrate and product channels, and reveal constrained interactions with triazole antifungal drugs that are important for drug design and understanding drug resistance.

Hepatic cytochrome P450 ubiquitination: conformational phosphodegrons for E2/E3 recognition?

Hepatic endoplasmic reticulum (ER) integral cytochromes P450 (P450s) are monooxygenases engaged in the biotransformation and elimination of endo- as well as xenobiotics. Of the human liver P450s, CYP3A4 is the major and most dominant catalyst responsible for the biotransformation of over 50% of clinically prescribed drugs. CYP2E1 metabolizes smaller molecular weight compounds (EtOH), carcinogens, environmental toxins, and endobiotics, and is justly implicated in various toxigenic/pathogenic mechanisms of human disease. Both P450s are notorious for their potential to generate pathogenic reactive oxygen species (ROS) during futile oxidative cycling and/or oxidative uncoupling. Such ROS not only oxidatively damage the P450 catalytic cage, but on their escape into the cytosol, also the P450 outer surface and any surrounding cell organelles. Given their ER-monotopic topology coupled with this high potential to acquire oxidative lesions in their cytosolic (C) domain, not surprisingly these P450 proteins exhibit shorter lifespans and are excellent prototype substrates of ER-associated degradation ("ERAD-C") pathway. Indeed, we have shown that both CYP3A4 and CYP2E1 incur ERAD-C, during which they are first phosphorylated by protein kinases A and C, which greatly enhance/accelerate their ubiquitination by UBC7/gp78 and UbcH5a/CHIP/Hsp70/Hsp40 E2/E3 ubiquitin ligase complexes. Such P450 phosphorylation occurs on Ser/Thr residues within linear sequences as well as spatially clustered acidic (Asp/Glu) residues. We propose that such S/T phosphorylation within these clusters creates negatively charged patches or conformational phosphodegrons for interaction with positively charged E2/E3 domains. Such P450 S/T phosphorylation we posit serves as a molecular switch to turn on its ubiquitination and ERAD-C.

Polymorphisms of CYP51A1 from cholesterol synthesis: associations with birth weight and maternal lipid levels and impact on CYP51 protein structure

We investigated the housekeeping cytochrome P450 CYP51A1 encoding lanosterol 14α-demethylase from cholesterol synthesis that was so far not directly linked to human disorders. By direct sequencing of CYP51A1 in 188 women with spontaneous preterm delivery and 188 unrelated preterm infants (gestational age <37 weeks) we identified 22 variants where 10 are novel and rare. In infants there were two novel CYP51A1 variants where damaging effects of p.Tyr145Asp from the substrate recognition region, but not p.Asn193Asp, were predicted by PolyPhen2 and SIFT. This was confirmed by molecular modeling showing that Tyr145Asp substitution results in changed electrostatic potential of the CYP51 protein surface and lengthened distance to the heme which prevents hydrogen bonding. The CYP51 Tyr145Asp mutation is rare and thus very interesting for further structure/function relationship studies. From the 12 identified known variants rs6465348 was chosen for family based association studies due to its high minor allele frequency. Interestingly, this CYP51A1 common variant associates with small for gestational age weight in newborns (p = 0.028) and lower blood total cholesterol and low density lipoprotein cholesterol levels in mothers in 2nd trimester of pregnancy (p = 0.042 and p = 0.046 respectively). Our results indicate a new link between a cholesterol synthesis gene CYP51A1 and pregnancy pathologies.

Identifying cytochrome p450 functional networks and their allosteric regulatory elements

Cytochrome P450 (CYP) enzymes play key roles in drug metabolism and adverse drug-drug interactions. Despite tremendous efforts in the past decades, essential questions regarding the function and activity of CYPs remain unanswered. Here, we used a combination of sequence-based co-evolutionary analysis and structure-based anisotropic thermal diffusion (ATD) molecular dynamics simulations to detect allosteric networks of amino acid residues and characterize their biological and molecular functions. We investigated four CYP subfamilies (CYP1A, CYP2D, CYP2C, and CYP3A) that are involved in 90% of all metabolic drug transformations and identified four amino acid interaction networks associated with specific CYP functionalities, i.e., membrane binding, heme binding, catalytic activity, and dimerization. Interestingly, we did not detect any co-evolved substrate-binding network, suggesting that substrate recognition is specific for each subfamily. Analysis of the membrane binding networks revealed that different CYP proteins adopt different membrane-bound orientations, consistent with the differing substrate preference for each isoform. The catalytic networks were associated with conservation of catalytic function among CYP isoforms, whereas the dimerization network was specific to different CYP isoforms. We further applied low-temperature ATD simulations to verify proposed allosteric sites associated with the heme-binding network and their role in regulating metabolic fate. Our approach allowed for a broad characterization of CYP properties, such as membrane interactions, catalytic mechanisms, dimerization, and linking these to groups of residues that can serve as allosteric regulators. The presented combined co-evolutionary analysis and ATD simulation approach is also generally applicable to other biological systems where allostery plays a role.

Correlating structure and function of drug-metabolizing enzymes: progress and ongoing challenges

This report summarizes a symposium sponsored by the American Society for Pharmacology and Experimental Therapeutics at Experimental Biology held April 20-24 in Boston, MA. Presentations discussed the status of cytochrome P450 (P450) knowledge, emphasizing advances and challenges in relating structure with function and in applying this information to drug design. First, at least one structure of most major human drug-metabolizing P450 enzymes is known. However, the flexibility of these active sites can limit the predictive value of one structure for other ligands. A second limitation is our coarse-grain understanding of P450 interactions with membranes, other P450 enzymes, NADPH-cytochrome P450 reductase, and cytochrome b5. Recent work has examined differential P450 interactions with reductase in mixed P450 systems and P450:P450 complexes in reconstituted systems and cells, suggesting another level of functional control. In addition, protein nuclear magnetic resonance is a new approach to probe these protein/protein interactions, identifying interacting b5 and P450 surfaces, showing that b5 and reductase binding are mutually exclusive, and demonstrating ligand modulation of CYP17A1/b5 interactions. One desired outcome is the application of such information to control drug metabolism and/or design selective P450 inhibitors. A final presentation highlighted development of a CYP3A4 inhibitor that slows clearance of human immunodeficiency virus drugs otherwise rapidly metabolized by CYP3A4. Although understanding P450 structure/function relationships is an ongoing challenge, translational advances will benefit from continued integration of existing and new biophysical approaches.

Characterizing the membrane-bound state of cytochrome P450 3A4: structure, depth of insertion, and orientation

Cytochrome P450 3A4 (CYP3A4) is the most abundant membrane-associated isoform of the P450 family in humans and is responsible for biotransformation of more than 50% of drugs metabolized in the body. Despite the large number of crystallographic structures available for CYP3A4, no structural information for its membrane-bound state at an atomic level is available. In order to characterize binding, depth of insertion, membrane orientation, and lipid interactions of CYP3A4, we have employed a combined experimental and simulation approach in this study. Taking advantage of a novel membrane representation, highly mobile membrane mimetic (HMMM), with enhanced lipid mobility and dynamics, we have been able to capture spontaneous binding and insertion of the globular domain of the enzyme into the membrane in multiple independent, unbiased simulations. Despite different initial orientations and positions of the protein in solution, all the simulations converged into the same membrane-bound configuration with regard to both the depth of membrane insertion and the orientation of the enzyme on the surface of the membrane. In tandem, linear dichroism measurements performed on CYP3A4 bound to Nanodisc membranes were used to characterize the orientation of the enzyme in its membrane-bound form experimentally. The heme tilt angles measured experimentally are in close agreement with those calculated for the membrane-bound structures resulted from the simulations, thereby verifying the validity of the developed model. Membrane binding of the globular domain in CYP3A4, which appears to be independent of the presence of the transmembrane helix of the full-length enzyme, significantly reshapes the protein at the membrane interface, causing conformational changes relevant to access tunnels leading to the active site of the enzyme.

[Metabolism of xenobiotics: beneficial and adverse effects]

The systems developed by living organisms for the metabolism of xenobiotics play a key role in the adaptation of living species to their chemical environment. Recent data about mammalian cytochrome P450 structures have led to a better understanding of the molecular basis for the adaptability of these enzymes to xenobiotics exhibiting highly variable structures. The action of these enzymes on xenobiotics leads to other beneficial effects such as the bioactivation of some drugs, but also to adverse effects with the formation of aggressive metabolites for the cell that are responsible for the appearance of many toxicities.

Structural diversity of eukaryotic membrane cytochrome p450s

X-ray crystal structures are available for 29 eukaryotic microsomal, chloroplast, or mitochondrial cytochrome P450s, including two non-monooxygenase P450s. These structures provide a basis for understanding structure-function relations that underlie their distinct catalytic activities. Moreover, structural plasticity has been characterized for individual P450s that aids in understanding substrate binding in P450s that mediate drug clearance.

Mechanism-based inhibition profiles of erythromycin and clarithromycin with cytochrome P450 3A4 genetic variants

Inhibition of cytochrome P450 (CYP) 3A4 is the major cause of drug-drug interactions (DDI). We have previously reported that the genetic variation of CYP3A4 significantly affected the inhibitory profiles of typical competitive inhibitors. In addition to competitive inhibition, some clinically significant DDI are attributable to mechanism-based inhibition (MBI). However, the differences in the MBI kinetics among CYP3A4 genetic variants remain to be characterized. In this study, we quantitatively investigated the inhibition kinetics of MBI inhibitors, erythromycin and clarithromycin, on the CYP3A4 variants CYP3A4.1, 4.2, 4.7, 4.16, and 4.18. The activity of CYP3A4 was assessed using testosterone 6β-hydroxylation with recombinant CYP3A4. Both erythromycin and clarithromycin decreased the activity of CYP3A4 in a time-dependent manner. The maximum inactivation rate constants, k(inact,max), of erythromycin for CYP3A4.2 and CYP3A4.7 were 0.5-fold that for CYP3A4.1, while that for CYP3A4.16 and CYP3A4.18 were similar to that for CYP3A4.1. The K(I) values of erythromycin for CYP3A4.2, 4.7, 4.16, and 4.18 were 1.2-, 0.4-, 2.2- and 0.72-fold those of CYP3A4.1, respectively. Similar results were obtained for clarithromycin. In conclusion, the inhibitory profiles of MBI inhibitors, as well as competitive inhibitors, may possibly differ among CYP3A4 variants. This difference may contribute to interindividual differences in the extent of DDI based on MBI.

π-π Stacking mediated drug-drug interactions in human CYP2E1

Because of having many low molecular mass substrates, CYP2E1 is of particular interests to the pharmaceutical industry. Many evidences showed that this enzyme can adopt multiple substrates to significantly reduce the oxidation rate of the substrates. The detailed mechanism for this observation is still unclear. In the current study, we employed GPU-accelerated molecular dynamics simulations to study the multiple-binding mode of human CYP2E1, with an aim of offering a mechanistic explanation for the unexplained multiple-substrate binding. Our results showed that Thr303 and Phe478 were key factors for the substrate recognition and multiple-substrate binding. The former can form a significant hydrogen bond to recognize and position the substrate in the productive binding orientation in the active site. The latter acted as a mediator for the substrate communications via π-π stacking interactions. In the multiple-binding mode, the aforementioned π-π stacking interactions formed by the aromatic rings of both substrates and Phe478 drove the first substrate far away from the catalytic center, orienting in an additional binding position and going against the substrate metabolism. All these findings could give atomic insights into the detailed mechanism for the multiple-substrate binding in human CYP2E1, providing useful information for the drug metabolism mechanism and personalized use of clinical drugs.

Structural insights into aldosterone synthase substrate specificity and targeted inhibition

Aldosterone is a major mineralocorticoid hormone that plays a key role in the regulation of electrolyte balance and blood pressure. Excess aldosterone levels can arise from dysregulation of the renin-angiotensin-aldosterone system and are implicated in the pathogenesis of hypertension and heart failure. Aldosterone synthase (cytochrome P450 11B2, CYP11B2) is the sole enzyme responsible for the production of aldosterone in humans. Blocking of aldosterone synthesis by mediating aldosterone synthase activity is thus a recently emerging pharmacological therapy for hypertension, yet a lack of structural information has limited this approach. Here, we present the crystal structures of human aldosterone synthase in complex with a substrate deoxycorticosterone and an inhibitor fadrozole. The structures reveal a hydrophobic cavity with specific features associated with corticosteroid recognition. The substrate binding mode, along with biochemical data, explains the high 11β-hydroxylase activity of aldosterone synthase toward both gluco- and mineralocorticoid formation. The low processivity of aldosterone synthase with a high extent of intermediates release might be one of the mechanisms of controlled aldosterone production from deoxycorticosterone. Although the active site pocket is lined by identical residues between CYP11B isoforms, most of the divergent residues that confer additional 18-oxidase activity of aldosterone synthase are located in the I-helix (vicinity of the O(2) activation path) and loops around the H-helix (affecting an egress channel closure required for retaining intermediates in the active site). This intrinsic flexibility is also reflected in isoform-selective inhibitor binding. Fadrozole binds to aldosterone synthase in the R-configuration, using part of the active site cavity pointing toward the egress channel. The structural organization of aldosterone synthase provides critical insights into the molecular mechanism of catalysis and enables rational design of more specific antihypertensive agents.

The synthesis, characterization, and application of ¹³C-methyl isocyanide as an NMR probe of heme protein active sites

The cytochromes P450 (CYPs) play a central role in a variety of important biological oxidations, such as steroid synthesis and the metabolism of xenobiotic compounds, including most drugs. Because CYPs are frequently assayed as drug targets or as anti-targets, tools that provide confirmation of active-site binding and information on binding orientation would be of great utility. Of greatest value are assays that are reasonably high throughput. Other heme proteins, too-such as the nitric oxide synthases (NOSs), with their importance in signaling, regulation of blood pressure, and involvement in the immune response-often display critical roles in the complex functions of many higher organisms, and also require improved assay methods. To this end, we have developed an analog of cyanide, with a (13)CH3-reporter group attached to make methyl isocyanide. We describe the synthesis and use of (13)C-methyl isocyanide as a probe of both bacterial (P450cam) and membrane-bound mammalian (CYP2B4) CYPs. The (13)C-methyl isocyanide probe can be used in a relatively high-throughput 1-D experiment to identify binders, but it can also be used to detect structural changes in the active site based on chemical shift changes, and potentially nuclear Overhauser effects between probe and inhibitor.

Directed evolution reveals requisite sequence elements in the functional expression of P450 2F1 in Escherichia coli

Cytochrome P450 2F1 (P450 2F1) is expressed exclusively in the human respiratory tract and is implicated in 3-methylindole (3MI)-induced pneumotoxicity via dehydrogenation of 3MI to a reactive electrophilic intermediate, 3-methyleneindolenine (3-MEI). Studies of P450 2F1 to date have been limited by the failure to express this enzyme in Escherichia coli. By contrast, P450 2F3, a caprine homologue that shares 84% sequence identity with P450 2F1 (86 amino acid differences), has been expressed in E. coli at yields greater than 250 nmol/L culture. We hypothesized that a limited number of sequence differences between P450s 2F1 and 2F3 could limit P450 2F1 expression in E. coli and that problematic P450 2F1 sequence elements could be identified by directed evolution. A library of P450 2F1/2F3 mutants was created by DNA family shuffling and screened for expression in E. coli. Three generations of DNA shuffling revealed a mutant (named JH_2F_F3_1_007) with 96.5% nucleotide sequence identity to P450 2F1 and which expressed 119 ± 40 pmol (n = 3, mean ± SD) hemoprotein in 1 mL microaerobic cultures. Across all three generations, two regions were observed where P450 2F3-derived sequence was consistently substituted for P450 2F1 sequence in expressing mutants, encoding nine amino acid differences between P450s 2F1 and 2F3: nucleotides 191-278 (amino acids 65-92) and 794-924 (amino acids 265-305). Chimeras constructed to specifically test the importance of these two regions confirmed that P450 2F3 sequence is essential in both regions for expression in E. coli but that other non-P450 2F1 sequence elements outside of these regions also improved the expression of mutant JH_2F_F3_1_007. Mutant JH_2F_F3_1_007 catalyzed the dehydrogenation of 3MI to 3-MEI as indicated by the observation of glutathione adducts after incubation in the presence of glutathione. The JH_2F_F3_1_007 protein differs from P450 2F1 at only 20 amino acids and should facilitate further studies of the structure-activity relationships of P450s of the 2F subfamily.

Cytochrome P450 1A1-mediated anticancer drug discovery: in silico findings

INTRODUCTION

Target-specific drugs may offer fewer side/adverse effects in comparison with other anticancer agents and thus save normal healthy cells to a greater extent. The selective overexpression of cytochrome P450 1A1 (CYP1A1) in tumor cells induces the metabolism of benzothiazole and aminoflavone compounds to their reactive species, which are responsible for DNA adduct formation and cell death. This review encompasses the novelty of CYP1A1 as an anticancer drug target and explores the possible in silico strategies that would be applicable in the discovery and development of future antitumor compounds.

AREAS COVERED

This review highlights the various ligand-based and target-based in silico methodologies that were efficiently used in exploration of CYP1A1 as a novel antitumor target. These methodologies include electronic structure analysis, CoMFA studies, homology modeling, molecular docking, molecular dynamics analysis, pharmacophore mapping and quantitative structure activity relationship (QSAR) studies. It also focuses on the various approaches used in the development of the lysyl amide prodrug of 5F-203 (NSC710305) and dimethanesulfonate salt of 5-aminoflavone (NSC710464) as clinical candidates from their less potent analogues.

EXPERT OPINION

Selective overexpression of CYP1A1 in cancer cells offers tumor-specific drug design to ameliorate the current adverse effects associated with existing antitumor agents. Medicinal chemistry and in vitro driven approaches, in combination with knowledge-based drug design and by using the currently available tools of in silico methodologies, would certainly make it possible to design and develop novel anticancer compounds targeting CYP1A1.

Plasticity of CYP2B enzymes: structural and solution biophysical methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

Classification and characterization of putative cytochrome P450 genes from Panax ginseng C. A. Meyer

In plants heme containing cytochrome P450 (P450) is a superfamily of monooxygenases that catalyze the addition of one oxygen atom from O2 into a substrate, with a substantial reduction of the other atom to water. The function of P450 families is attributed to chemical defense mechanism under terrestrial environmental conditions; several are involved in secondary and hormone metabolism. However, the evolutionary relationships of P450 genes in Panax ginseng remain largely unknown. In the present study, data mining methods were implemented and 116 novel putative P450 genes were identified from Expressed Sequence Tags (ESTs) of a ginseng database. These genes were classified into four clans and 22 families by sequence similarity conducted at amino acid level. The representative putative P450 sequences of P. ginseng and known P450 family from other plants were used to construct a phylogenetic tree. By comparing with other genomes, we found that most of the P450 genes from P. ginseng can be found in other dicot species. Depending on P450 family functions, seven P450 genes were selected, and for that organ specific expression, abiotic, and biotic studies were performed by quantitative reverse transcriptase-polymerase chain reaction. Different genes were found to be expressed differently in different organs. Biotic stress and abiotic stress transcript level was regulated diversely, and upregulation of P450 genes indicated the involvement of certain genes under stress conditions. The upregulation of the P450 genes under methyl jasmonate and fungal stress justifies the involvement of specific genes in secondary metabolite biosynthesis. Our results provide a foundation for further elucidating the actual function and role of P450 involved in various biochemical pathways in P. ginseng.

Impact of six fruits--banana, guava, mangosteen, pineapple, ripe mango and ripe papaya--on murine hepatic cytochrome P450 activities

The effects of six Thai fruits, namely banana, guava, mangosteen, pineapple, ripe mango and ripe papaya, on cytochrome P450 (P450) activities were investigated. The median inhibitory concentrations (IC(50) ) of each of the fruit juices on CYP1A1, CYP1A2, CYP2E1 and CYP3A11 activities were determined. Pineapple juice showed the strongest inhibitory effect against all the evaluated P450 isozyme activities in mouse hepatic microsomes, followed by mangosteen, guava, ripe mango, ripe papaya and banana. The study was further performed in male ICR mice given pineapple juice intragastrically at doses of 10, 20 and 40 mg kg(-1) per day for 7 or 28 days. In a concentration-dependent fashion, the pineapple juice raised ethoxyresorufin O-deethylase, aniline hydroxylase and erythromycin N-demethylase activities, which are marker enzymatic reactions responsible for CYP1A1, CYP2E1 and CYP3A11, respectively. The effect of pineapple juice on the expression of CYP1A1, CYP2E1 and CYP3A11 mRNAs corresponded to their enzymatic activities. However, the pineapple juice significantly decreased methoxyresorufin O-demethylase activity. These observations supported that the six Thai fruits were a feasible cause of food-drug interaction or adverse drug effects owing to their potential to modify several essential P450 activities. Individuals consuming large quantities of pineapple for long periods of time should be cautioned of these potential adverse effects.

Prediction of sites of metabolism in a substrate molecule, instanced by carbamazepine oxidation by CYP3A4

In drug discovery process, improvement of ADME/Tox properties of lead compounds including metabolic stability is critically important. Cytochrome P450 (CYP) is one of the major metabolizing enzymes and the prediction of sites of metabolism (SOM) on the given lead compounds is key information to modify the compounds to be more stable against metabolism. There are two factors essentially important in SOM prediction. First is accessibility of each substrate atom to the oxygenated Fe atom of heme in a CYP protein, and the other is the oxidative reactivity of each substrate atom. To predict accessibility of substrate atoms to the heme iron, conventional protein-rigid docking simulations have been applied. However, the docking simulations without consideration of protein flexibility often lead to incorrect answers in the case of very flexible proteins such as CYP3A4. In this study, we demonstrated an approach utilizing molecular dynamics (MD) simulation for SOM prediction in which multiple MD runs were executed using different initial structures. We applied this strategy to CYP3A4 and carbamazepine (CBZ) complex. Through 10 ns MD simulations started from five different CYP3A4-CBZ complex models, our approach correctly predicted SOM observed in experiments. The experimentally known epoxidized sites of CBZ by CYP3A4 were successfully predicted as the most accessible sites to the heme iron that was judged from a numerical analysis of calculated ΔG(binding) and the frequency of appearance. In contrast, the predictions using protein-rigid docking methods hardly provided the correct SOM due to protein flexibility or inaccuracy of the scoring functions. Our strategy using MD simulation with multiple initial structures will be one of the reliable methods for SOM prediction.

Homology modeling of the three membrane proteins of the dhurrin metabolon: catalytic sites, membrane surface association and protein-protein interactions

Formation of metabolons (macromolecular enzyme complexes) facilitates the channelling of substrates in biosynthetic pathways. Metabolon formation is a dynamic process in which transient structures mediated by weak protein-protein interactions are formed. In Sorghum, the cyanogenic glucoside dhurrin is derived from l-tyrosine in a pathway involving the two cytochromes P450 (CYPs) CYP79A1 and CYP71E1, a glucosyltransferase (UGT85B1), and the redox partner NADPH-dependent cytochrome P450 reductase (CPR). Experimental evidence suggests that the enzymes of this pathway form a metabolon. Homology modeling of the three membrane bound proteins was carried out using the Sybyl software and available relevant crystal structures. Residues involved in tight positioning of the substrates and intermediates in the active sites of CYP79A1 and CYP71E1 were identified. In both CYPs, hydrophobic surface domains close to the N-terminal trans-membrane anchor and between the F' and G helices were identified as involved in membrane anchoring. The proximal surface of both CYPs showed positively charged patches complementary to a negatively charged bulge on CPR carrying the FMN domain. A patch of surface exposed, positively charged amino acid residues positioned on the opposite face of the membrane anchor was identified in CYP71E1 and might be involved in binding UGT85B1 via a hypervariable negatively charged loop in this protein.

Cytochromes p450

There are 244 cytochrome P450 genes (and 28 pseudogenes) in the Arabidopsis genome. P450s thus form one of the largest gene families in plants. Contrary to what was initially thought, this family diversification results in very limited functional redundancy and seems to mirror the complexity of plant metabolism. P450s sometimes share less than 20% identity and catalyze extremely diverse reactions leading to the precursors of structural macromolecules such as lignin, cutin, suberin and sporopollenin, or are involved in biosynthesis or catabolism of all hormone and signaling molecules, of pigments, odorants, flavors, antioxidants, allelochemicals and defense compounds, and in the metabolism of xenobiotics. The mechanisms of gene duplication and diversification are getting better understood and together with co-expression data provide leads to functional characterization.

Discovery of orteronel (TAK-700), a naphthylmethylimidazole derivative, as a highly selective 17,20-lyase inhibitor with potential utility in the treatment of prostate cancer

A novel naphthylmethylimidazole derivative 1 and its related compounds were identified as 17,20-lyase inhibitors. Based on the structure-activity relationship around the naphthalene scaffold and the results of a docking study of 1a in the homology model of 17,20-lyase, the 6,7-dihydro-5H-pyrrolo[1,2-c]imidazole derivative (+)-3c was synthesized and identified as a potent and highly selective 17,20-lyase inhibitor. Biological evaluation of (+)-3c at a dose of 1mg/kg in a male monkey model revealed marked reductions in both serum testosterone and dehydroepiandrosterone concentrations. Therefore, (+)-3c (termed orteronel [TAK-700]) was selected as a candidate for clinical evaluation and is currently in phase III clinical trials for the treatment of castration-resistant prostate cancer.

Structures of cytochrome P450 2B6 bound to 4-benzylpyridine and 4-(4-nitrobenzyl)pyridine: insight into inhibitor binding and rearrangement of active site side chains

The biochemical, biophysical, and structural analysis of the cytochrome P450 2B subfamily of enzymes has provided a wealth of information regarding conformational plasticity and substrate recognition. The recent X-ray crystal structure of the drug-metabolizing P450 2B6 in complex with 4-(4-chlorophenyl)imidazole (4-CPI) yielded the first atomic view of this human enzyme. However, knowledge of the structural basis of P450 2B6 specificity and inhibition has remained limited. In this study, structures of P450 2B6 were determined in complex with the potent inhibitors 4-benzylpyridine (4-BP) and 4-(4-nitrobenzyl)pyridine (4-NBP). Comparison of the present structures with the previous P450 2B6-4-CPI complex showed that reorientation of side chains of the active site residue Phe206 on the F-helix and Phe297 on the I-helix was necessary to accommodate the inhibitors. However, P450 2B6 does not require any major side chain rearrangement to bind 4-NBP compared with 4-BP, and the enzyme provides no hydrogen-bonding partners for the polar nitro group of 4-NBP within the hydrophobic active site. In addition, on the basis of these new structures, substitution of residue 172 with histidine as observed in the single nucleotide polymorphism Q172H and in P450 2B4 may contribute to a hydrogen bonding network connecting the E- and I-helices, thereby stabilizing active site residues on the I-helix. These results provide insight into the role of active site side chains upon inhibitor binding and indicate that the recognition of the benzylpyridines in the closed conformation structure of P450 2B6 is based solely on hydrophobicity, size, and shape.

Structure and dynamics of the membrane-bound cytochrome P450 2C9

The microsomal, membrane-bound, human cytochrome P450 (CYP) 2C9 is a liver-specific monooxygenase essential for drug metabolism. CYPs require electron transfer from the membrane-bound CYP reductase (CPR) for catalysis. The structural details and functional relevance of the CYP-membrane interaction are not understood. From multiple coarse grained molecular simulations started with arbitrary configurations of protein-membrane complexes, we found two predominant orientations of CYP2C9 in the membrane, both consistent with experiments and conserved in atomic-resolution simulations. The dynamics of membrane-bound and soluble CYP2C9 revealed correlations between opening and closing of different tunnels from the enzyme's buried active site. The membrane facilitated the opening of a tunnel leading into it by stabilizing the open state of an internal aromatic gate. Other tunnels opened selectively in the simulations of product-bound CYP2C9. We propose that the membrane promotes binding of liposoluble substrates by stabilizing protein conformations with an open access tunnel and provide evidence for selective substrate access and product release routes in mammalian CYPs. The models derived here are suitable for extension to incorporate other CYPs for oligomerization studies or the CYP reductase for studies of the electron transfer mechanism, whereas the modeling procedure is generally applicable to study proteins anchored in the bilayer by a single transmembrane helix.

We describe the first atomic-detail models and simulations of a full-length, membrane-bound mammalian cytochrome P450. To date, all the structural studies of microsomal, drug-metabolizing cytochrome P450s have been performed using engineered, solubilized forms of the enzymes and it is not yet understood how the membrane influences their structure, dynamics, and ability to bind substrates. We focused on CYP2C9, the second most abundant cytochrome P450 in the human liver which oxidizes 20% of all marketed drugs. Here, we have derived models of CYP2C9-membrane complexes with a modeling procedure based on molecular dynamics simulations started with arbitrary configurations of the protein in the membrane and performed using both coarse grained and atomic-detail molecular representations. This procedure is expected to be generally applicable to proteins that are anchored in the membrane with a single transmembrane helix. The models and simulations provide evidence for selective substrate access and product release routes in membrane-bound CYPs. This observation may contribute to new strategies to manipulate the activity of CYPs and other enzymes with buried active sites. We anticipate that this study will bring about a paradigm shift towards studying microsomal CYPs as dynamic structures in their natural, lipid environment rather than in artificially solubilized forms.

Membrane position of ibuprofen agrees with suggested access path entrance to cytochrome P450 2C9 active site

Cytochrome P450 2C9 (CYP2C9) is a membrane-anchored human microsomal protein involved in the drug metabolism in liver. CYP2C9 consists of an N-terminal transmembrane anchor and a catalytic cytoplasmic domain. While the structure of the catalytic domain is well-known from X-ray experiments, the complete structure and its incorporation into the membrane remains unsolved. We constructed an atomistic model of complete CYP2C9 in a dioleoylphosphatidylcholine membrane and evolved it by molecular dynamics simulations in explicit water on a 100+ ns time-scale. The model agrees well with known experimental data about membrane positioning of cytochromes P450. The entry to the substrate access channel is proposed to be facing the membrane interior while the exit of the product egress channel is situated above the interface pointing toward the water phase. The positions of openings of the substrate access and product egress channels correspond to free energy minima of CYP2C9 substrate ibuprofen and its metabolite in the membrane, respectively.

Structural features of cytochromes P450 and ligands that affect drug metabolism as revealed by X-ray crystallography and NMR

Cytochromes P450 (P450s) play a major role in the clearance of drugs, toxins, and environmental pollutants. Additionally, metabolism by P450s can result in toxic or carcinogenic products. The metabolism of pharmaceuticals by P450s is a major concern during the design of new drug candidates. Determining the interactions between P450s and compounds of very diverse structures is complicated by the variability in P450-ligand interactions. Understanding the protein structural elements and the chemical attributes of ligands that dictate their orientation in the P450 active site will aid in the development of effective and safe therapeutic agents. The goal of this review is to describe P450-ligand interactions from two perspectives. The first is the various structural elements that microsomal P450s have at their disposal to assume the different conformations observed in X-ray crystal structures. The second is P450-ligand dynamics analyzed by NMR relaxation studies.

Crystal structure of H2O2-dependent cytochrome P450SPalpha with its bound fatty acid substrate: insight into the regioselective hydroxylation of fatty acids at the alpha position

Cytochrome P450(SPα) (CYP152B1) isolated from Sphingomonas paucimobilis is the first P450 to be classified as a H(2)O(2)-dependent P450. P450(SPα) hydroxylates fatty acids with high α-regioselectivity. Herein we report the crystal structure of P450(SPα) with palmitic acid as a substrate at a resolution of 1.65 Å. The structure revealed that the C(α) of the bound palmitic acid in one of the alternative conformations is 4.5 Å from the heme iron. This conformation explains the highly selective α-hydroxylation of fatty acid observed in P450(SPα). Mutations at the active site and the F-G loop of P450(SPα) did not impair its regioselectivity. The crystal structures of mutants (L78F and F288G) revealed that the location of the bound palmitic acid was essentially the same as that in the WT, although amino acids at the active site were replaced with the corresponding amino acids of cytochrome P450(BSβ) (CYP152A1), which shows β-regioselectivity. This implies that the high regioselectivity of P450(SPα) is caused by the orientation of the hydrophobic channel, which is more perpendicular to the heme plane than that of P450(BSβ).

Structural basis for pregnenolone biosynthesis by the mitochondrial monooxygenase system

In humans, the precursor to all steroid hormones, pregnenolone, is synthesized from cholesterol by an enzyme complex comprising adrenodoxin reductase (AdR), adrenodoxin (Adx), and a cytochrome P450 (P450scc or CYP11A1). This complex not only plays a key role in steroidogenesis, but also has long been a model to study electron transfer, multistep catalysis, and C-C bond cleavage performed by monooxygenases. Detailed mechanistic understanding of these processes has been hindered by a lack of structural information. Here we present the crystal structure of the complex of human Adx and CYP11A1--the first of a complex between a eukaryotic CYP and its redox partner. The structures with substrate and a series of reaction intermediates allow us to define the mechanism underlying sequential hydroxylations of the cholesterol and suggest the mechanism of C-C bond cleavage. In the complex the [2Fe-2S] cluster of Adx is positioned 17.4 Å away from the heme iron of CYP11A1. This structure suggests that after an initial protein-protein association driven by electrostatic forces, the complex adopts an optimized geometry between the redox centers. Conservation of the interaction interface suggests that this mechanism is common for all mitochondrial P450s.

Gene engineering, purification, crystallization and preliminary X-ray diffraction of cytochrome P450 p-coumarate-3-hydroxylase (C3H), the Arabidopsis membrane protein

Cytochrome P450s (P450s) are the most versatile biological catalysts in plants; however, because the structure of the P450s has not been fully established, their broad substrate specificity has been limitedly discussed. p-coumarate-3-hydroxylase (C3H) is an essential enzyme for the biosynthesis of phenolic natural products in plants, but all attempts to express and purify C3H, have failed. In this research, we developed a bacterial expression of Arabidopsis C3H by combinational mutagenesis and purified C3H as a catalytically active form. The modified C3H could be purified in the absence of detergent, and crystallized in two forms (orthorhombic and trigonal space group) under different conditions. X-ray diffraction was processed to a 4.0 Å resolution (first type crystal) and a 3.8 Å resolution (second type crystal). Although the diffraction results of C3H(mod) crystals are not enough to determine crystallographic structure due to low resolution, the simplicity and rapidity of this technology are competitive advantages in comparison with other methods, and may contribute to structural analyses of other membrane proteins including P450s family.

Why human cytochrome P450c21 is a progesterone 21-hydroxylase

Human cytochrome P450c21 (steroid 21-hydroxylase, CYP21A2) catalyzes the 21-hydroxylation of progesterone (P4) and its preferred substrate 17α-hydroxyprogestrone (17OHP4). CYP21A2 activities, which are required for cortisol and aldosterone biosynthesis, involve the formation of energetically disfavored primary carbon radicals. Therefore, we hypothesized that the binding of P4 and 17OHP4 to CYP21A2 restricts access of the reactive heme-oxygen complex to the C-21 hydrogen atoms, suppressing oxygenation at kinetically more favorable sites such as C-17 and C-16, which are both hydroxylated by cytochrome P450c17 (CYP17A1). We reasoned that expansion of the CYP21A2 substrate-binding pocket would increase substrate mobility and might yield additional hydroxylation activities. We built a computer model of CYP21A2 based principally on the crystal structure of CYP2C5, which also 21-hydroxylates P4. Molecular dynamics simulations indicate that binding of the steroid nucleus perpendicular to the plane of the CYP21A2 heme ring limits access of the heme oxygen to the C-21 hydrogen atoms. Residues L107, L109, V470, I471, and V359 were found to contribute to the CYP21A2 substate-binding pocket. Mutation of V470 and I471 to alanine or glycine preserved P4 21-hydroxylase activity, and mutations of L107 or L109 were inactive. Mutations V359A and V359G, in contrast, acquired 16α-hydroxylase activity, accounting for 40% and 90% of the P4 metabolites, respectively. We conclude that P4 binds to CYP21A2 in a fundamentally different orientation than to CYP17A1 and that expansion of the CYP21A2 substrate-binding pocket allows additional substrate trajectories and metabolic switching.

Structure and function of cytochromes P450 2B: from mechanism-based inactivators to X-ray crystal structures and back

This article reviews work from the author dating back to 1978 and focuses on the structural basis of cytochrome P450 (P450) function using available contemporary techniques. Early studies used mechanism-based inactivators that bound to the protein moiety of hepatic P450s to try to localize the active site. Subsequent studies used cDNA cloning, heterologous expression, site-directed mutagenesis, and homology modeling based on multiple bacterial P450 X-ray crystal structures to predict the active sites of CYP2B enzymes with considerable accuracy. Breakthroughs in engineering and expression of mammalian P450s enabled us to determine our first X-ray crystal structure of ligand-free rabbit CYP2B4. To date, we have solved 11 CYP2B4 and three human CYP2B6 structures, which represent four significantly different conformations. The plasticity of CYP2B4 has been confirmed by deuterium exchange mass spectrometry and is substantiated by molecular dynamics simulations. In addition to major movement of secondary structure elements, more subtle reorientation of active site side chains, especially Phe206, Phe297, and Glu301, contributes to the ability of CYP2B enzymes to bind various ligands. Isothermal titration calorimetry has proven to be a useful tool for studying the thermodynamics of ligand binding to CYP2B4 and CYP2B6, and NMR has enabled study of ligand binding orientation in solution as an adjunct to X-ray crystallography. A major challenge remains to harness the power of the various approaches to facilitate prediction of CYP2B specificity and inhibition.

Cytochrome P450: taming a wild type enzyme

Protein engineering of cytochrome P450 monooxygenases (P450s) has been very successful in generating valuable non-natural activities and properties, allowing these powerful catalysts to be used for the synthesis of drug metabolites and in biosynthetic pathways for the production of precursors of artemisinin and paclitaxel. Collected experience indicates that the P450s are highly 'evolvable' - they are particularly robust to mutation in their active sites and readily accept new substrates and exhibit new selectivities. Their ability to adapt to new challenges upon mutation may reflect the nonpolar nature of their active sites as well as their high degree of conformational variability.