Functional interaction of cytochrome P450 with its redox partners: a critical assessment and update of the topology of predicted contact regions

Systems engineering Escherichia coli for efficient production p-coumaric acid from glucose

P-coumaric acid (p-CA), a pant metabolite with antioxidant and anti-inflammatory activity, is extensively utilized in biomedicine, food, and cosmetics industry. In this study, a synthetic pathway (PAL) for p-CA was designed, integrating three enzymes (AtPAL2, AtC4H, AtATR2) into a higher l-phenylalanine-producing strain Escherichia coli PHE05. However, the lower soluble expression and activity of AtC4H in the PAL pathway was a bottleneck for increasing p-CA titers. To overcome this limitation, the soluble expression of AtC4H was enhanced through N-terminal modifications. And an optimal mutant, AtC4HL373T/G211H, which exhibited a 4.3-fold higher kcat/Km value compared to the wild type, was developed. In addition, metabolic engineering strategies were employed to increase the intracellular NADPH pool. Overexpression of ppnk in engineered E. coli PHCA20 led to a 13.9-folds, 1.3-folds, and 29.1% in NADPH content, the NADPH/NADP+ ratio and p-CA titer, respectively. These optimizations significantly enhance p-CA production, in a 5-L fermenter using fed-batch fermentation, the p-CA titer, yield and productivity of engineered strain E. coli PHCA20 were 3.09 g/L, 20.01 mg/g glucose, and 49.05 mg/L/h, respectively. The results presented here provide a novel way to efficiently produce the plant metabolites using an industrial strain.

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

Unravelling the role of transient redox partner complexes in P450 electron transfer mechanics

The molecular evolution of cytochromes P450 and associated redox-driven oxidative catalysis remains a mystery in biology. It is widely believed that sterol 14α-demethylase (CYP51), an essential enzyme of sterol biosynthesis, is the ancestor of the whole P450 superfamily given its conservation across species in different biological kingdoms. Herein we have utilized X-ray crystallography, molecular dynamics simulations, phylogenetics and electron transfer measurements to interrogate the nature of P450-redox partner binding using the naturally occurring fusion protein, CYP51-ferredoxin found in the sterol-producing bacterium Methylococcus capsulatus. Our data advocates that the electron transfer mechanics in the M. capsulatus CYP51-ferredoxin fusion protein involves an ensemble of ferredoxin molecules in various orientations and the interactions are transient. Close proximity of ferredoxin, however, is required to complete the substrate-induced large-scale structural switch in the P450 domain that enables proton-coupled electron transfer and subsequent oxygen scission and catalysis. These results have fundamental implications regarding the early evolution of electron transfer proteins and for the redox reactions in the early steps of sterol biosynthesis. They also shed new light on redox protein mechanics and the subsequent diversification of the P450 electron transfer machinery in nature.

Structure-based virtual screening of CYP1A1 inhibitors: towards rapid tier-one assessment of potential developmental toxicants

Cytochrome P450 1A1 (CYP1A1) metabolizes estrogens, melatonin, and other key endogenous signaling molecules critical for embryonic/fetal development. The enzyme has increasing expression during pregnancy, and its inhibition or knockout increases embryonic/fetal lethality and/or developmental problems. Here, we present a virtual screening model for CYP1A1 inhibitors based on the orthosteric and predicted allosteric sites of the enzyme. Using 1001 reference compounds with CYP1A1 activity data, we optimized the decision thresholds of our model and classified the training compounds with 68.3% balanced accuracy (91.0% sensitivity and 45.7% specificity). We applied our final model to 11 known CYP1A1 orthosteric binders and related compounds, and found that our ranking of the known orthosteric binders generally agrees with the relative activity of CYP1A1 in metabolizing these compounds. We also applied the model to 22 new test compounds with unknown/unclear CYP1A1 inhibitory activity, and predicted 16 of them are CYP1A1 inhibitors. The CYP1A1 potency and modes of inhibition of these 22 compounds were experimentally determined. We confirmed that most predicted inhibitors, including drugs contraindicated during pregnancy (amiodarone, bicalutamide, cyproterone acetate, ketoconazole, and tamoxifen) and environmental agents suspected to be endocrine disruptors (bisphenol A, diethyl and dibutyl phthalates, and zearalenone), are indeed potent inhibitors of CYP1A1. Our results suggest that virtual screening may be used as a rapid tier-one method to screen for potential CYP1A1 inhibitors, and flag them out for further experimental evaluations.

Enzymatic Characterization and Comparison of Two Steroid Hydroxylases CYP154C3-1 and CYP154C3-2 from Streptomyces Species

Bacterial cytochrome P450 (CYP) enzymes are responsible for the hydroxylation of diverse endogenous substances with a heme molecule used as a cofactor. This study characterized two CYP154C3 proteins from Streptomyces sp. W2061 (CYP154C3-1) and Streptomyces sp. KCCM40643 (CYP154C3-2). The enzymatic activity assays of both CYPs conducted using heterologous redox partners' putidaredoxin and putidaredoxin reductase showed substrate flexibility with different steroids and exhibited interesting product formation patterns. The enzymatic characterization revealed good activity over a pH range of 7.0 to 7.8 and the optimal temperature range for activity was 30 to 37°C. The major product was the C16-hydroxylated product and the kinetic profiles and patterns of the generated hydroxylated products differed between the two enzymes. Both enzymes showed a higher affinity toward progesterone, with CYP154C3-1 demonstrating slightly higher activity than CYP154C3-2 for most of the substrates. Oxidizing agents (diacetoxyiodo) benzene (PIDA) and hydrogen peroxide (H₂O₂) were also utilized to actively support the redox reactions, with optimum conversion achieved at concentrations of 3 mM and 65 mM, respectively. The oxidizing agents affected the product distribution, influencing the type and selectivity of the CYP-catalyzed reaction. Additionally, CYP154C3s also catalyzed the C-C bond cleavage of steroids. Therefore, CYP154C3s may be a good candidate for the production of modified steroids for various biological uses.

An overview of the factors playing a role in cytochrome P450 monooxygenase and ferredoxin interactions

Cytochrome P450 monooxygenases (CYPs/P450s) are heme-thiolate proteins that are ubiquitously present in organisms, including non-living entities such as viruses. With the exception of self-sufficient P450s, all other P450 enzymes need electrons to perform their enzymatic activity and these electrons are supplied by P450 redox proteins. Different types of P450 redox proteins can be found in organisms and are classified into different classes. Bacterial P450s (class I) receive electrons from ferredoxins which are iron-sulfur cluster proteins. The presence of more than one copy and different types of ferredoxins within a bacterial species poses fundamental questions about the selectivity of P450s and ferredoxins in relation to each other. Apart from transferring electrons, ferredoxins have also been found to modulate P450 functions. Achieving an understanding of the interaction between ferredoxins and P450s is required to harness their biotechnological potential for designing a universal electron transfer protein. A brief overview of factors playing a role in ferredoxin and P450 interactions is presented in this review article.

The Role of the FMN-Domain of Human Cytochrome P450 Oxidoreductase in Its Promiscuous Interactions With Structurally Diverse Redox Partners

NADPH cytochrome P450 oxidoreductase (CPR) is the obligatory electron supplier that sustains the activity of microsomal cytochrome P450 (CYP) enzymes. The variant nature of the isoform-specific proximal interface of microsomal CYPs indicates that CPR is capable of multiple degenerated interactions with CYPs for electron transfer, through different binding mechanisms, and which are still not well-understood. Recently, we showed that CPR dynamics allows formation of open conformations that can be sampled by its structurally diverse redox partners in a CYP-isoform dependent manner. To further investigate the role of the CPR FMN-domain in effective binding of CPR to its diverse acceptors and to clarify the underlying molecular mechanisms, five different CPR-FMN-domain random mutant libraries were created. These libraries were screened for mutants with increased activity when combined with specific CYP-isoforms. Seven CPR-FMN-domain mutants were identified, supporting a gain in activity for CYP1A2 (P117H, G144C, A229T), 2A6 (P117L/L125V, G175D, H183Y), or 3A4 (N151D). Effects were evaluated using extended enzyme kinetic analysis, cytochrome b 5 competition, ionic strength effect on CYP activity, and structural analysis. Mutated residues were located either in or adjacent to several acidic amino acid stretches - formerly indicated to be involved in CPR:CYP interactions - or close to two tyrosine residues suggested to be involved in FMN binding. Several of the identified positions co-localize with mutations found in naturally occurring CPR variants that were previously shown to cause CYP-isoform-dependent effects. The mutations do not seem to significantly alter the geometry of the FMN-domain but are likely to cause very subtle alterations leading to improved interaction with a specific CYP. Overall, these data suggest that CYPs interact with CPR using an isoform specific combination of several binding motifs of the FMN-domain.

Probing the Role of the Hinge Segment of Cytochrome P450 Oxidoreductase in the Interaction with Cytochrome P450

NADPH-cytochrome P450 reductase (CPR) is the unique redox partner of microsomal cytochrome P450s (CYPs). CPR exists in a conformational equilibrium between open and closed conformations throughout its electron transfer (ET) function. Previously, we have shown that electrostatic and flexibility properties of the hinge segment of CPR are critical for ET. Three mutants of human CPR were studied (S243P, I245P and R246A) and combined with representative human drug-metabolizing CYPs (isoforms 1A2, 2A6 and 3A4). To probe the effect of these hinge mutations different experimental approaches were employed: CYP bioactivation capacity of pre-carcinogens, enzyme kinetic analysis, and effect of the ionic strength and cytochrome b₅ (CYB5) on CYP activity. The hinge mutations influenced the bioactivation of pre-carcinogens, which seemed CYP isoform and substrate dependent. The deviations of Michaelis-Menten kinetic parameters uncovered tend to confirm this discrepancy, which was confirmed by CYP and hinge mutant specific salt/activity profiles. CPR/CYB5 competition experiments indicated a less important role of affinity in CPR/CYP interaction. Overall, our data suggest that the highly flexible hinge of CPR is responsible for the existence of a conformational aggregate of different open CPR conformers enabling ET-interaction with structural varied redox partners.

The Hinge Segment of Human NADPH-Cytochrome P450 Reductase in Conformational Switching: The Critical Role of Ionic Strength

NADPH-cytochrome P450 reductase (CPR) is a redox partner of microsomal cytochromes P450 and is a prototype of the diflavin reductase family. CPR contains 3 distinct functional domains: a FMN-binding domain (acceptor reduction), a linker (hinge), and a connecting/FAD domain (NADPH oxidation). It has been demonstrated that the mechanism of CPR exhibits an important step in which it switches from a compact, closed conformation (locked state) to an ensemble of open conformations (unlocked state), the latter enabling electron transfer to redox partners. The conformational equilibrium between the locked and unlocked states has been shown to be highly dependent on ionic strength, reinforcing the hypothesis of the presence of critical salt interactions at the interface between the FMN and connecting FAD domains. Here we show that specific residues of the hinge segment are important in the control of the conformational equilibrium of CPR. We constructed six single mutants and two double mutants of the human CPR, targeting residues G240, S243, I245 and R246 of the hinge segment, with the aim of modifying the flexibility or the potential ionic interactions of the hinge segment. We measured the reduction of cytochrome c at various salt concentrations of these 8 mutants, either in the soluble or membrane-bound form of human CPR. All mutants were found capable of reducing cytochrome c yet with different efficiency and their maximal rates of cytochrome c reduction were shifted to lower salt concentration. In particular, residue R246 seems to play a key role in a salt bridge network present at the interface of the hinge and the connecting domain. Interestingly, the effects of mutations, although similar, demonstrated specific differences when present in the soluble or membrane-bound context. Our results demonstrate that the electrostatic and flexibility properties of the hinge segment are critical for electron transfer from CPR to its redox partners.

Modulating proposed electron transfer pathways in P450BM3 led to improved activity and coupling efficiency

Electrochemical in vitro reduction of P450 enzymes is a promising alternative to in vivo applications. Previously we presented three engineered P450_BM3 variants for aniline hydroxylation, equipped with a carbon nanotube binding-peptide (CNT-tag) for self-assembly on CNT electrodes. Compared to wildtype P450_BM3 the NADPH-dependent activity was enhanced, but the coupling efficiency remained low. For P450_BM3 Verma, Schwaneberg and Roccatano (2014, Biopolymers 101, 197-209) calculated putative electron transfer pathways (eTPs) by MD simulations. We hypothesised that knockouts of these transfer pathways would alter the coupling efficiency of the system. The results revealed no improved system for the electrically-driven P450s. For the NADPH-driven P450s, however, the most active eTP-mutant showed a 13-fold increased activity and a 32-fold elevated coupling efficiency using NADPH as reducing equivalent. This suggests an alternative principle of electron transport for the reduction by NADPH and an electrode, respectively. The work presents moreover a tool to improve the coupling and activity of P450s with non-natural substrates.

NADH:Cytochrome b5 Reductase and Cytochrome b5 Can Act as Sole Electron Donors to Human Cytochrome P450 1A1-Mediated Oxidation and DNA Adduct Formation by Benzo[a]pyrene

Benzo[a]pyrene (BaP) is a human carcinogen that covalently binds to DNA after activation by cytochrome P450 (P450). Here, we investigated whether NADH:cytochrome b5 reductase (CBR) in the presence of cytochrome b5 can act as sole electron donor to human P450 1A1 during BaP oxidation and replace the canonical NADPH:cytochrome P450 reductase (POR) system. We also studied the efficiencies of the coenzymes of these reductases, NADPH as a coenzyme of POR, and NADH as a coenzyme of CBR, to mediate BaP oxidation. Two systems containing human P450 1A1 were utilized: human recombinant P450 1A1 expressed with POR, CBR, epoxide hydrolase, and cytochrome b5 in Supersomes and human recombinant P450 1A1 reconstituted with POR and/or with CBR and cytochrome b5 in liposomes. BaP-9,10-dihydrodiol, BaP-7,8-dihydrodiol, BaP-1,6-dione, BaP-3,6-dione, BaP-9-ol, BaP-3-ol, a metabolite of unknown structure, and two BaP-DNA adducts were generated by the P450 1A1-Supersomes system, both in the presence of NADPH and in the presence of NADH. The major BaP-DNA adduct detected by (32)P-postlabeling was characterized as 10-(deoxyguanosin-N(2)-yl)-7,8,9-trihydroxy-7,8,9,10-tetrahydro-BaP (assigned adduct 1), while the minor adduct is probably a guanine adduct derived from 9-hydroxy-BaP-4,5-epoxide (assigned adduct 2). BaP-3-ol as the major metabolite, BaP-9-ol, BaP-1,6-dione, BaP-3,6-dione, an unknown metabolite, and adduct 2 were observed in the system using P450 1A1 reconstituted with POR plus NADPH. When P450 1A1 was reconstituted with CBR and cytochrome b5 plus NADH, BaP-3-ol was the predominant metabolite too, and an adduct 2 was also generated. Our results demonstrate that the NADH/cytochrome b5/CBR system can act as the sole electron donor both for the first and second reduction of P450 1A1 during the oxidation of BaP in vitro. They suggest that NADH-dependent CBR can replace NADPH-dependent POR in the P450 1A1-catalyzed metabolism of BaP.

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.

Catalytically relevant electrostatic interactions of cytochrome P450c17 (CYP17A1) and cytochrome b5

Two acidic residues, Glu-48 and Glu-49, of cytochrome b5 (b5) are essential for stimulating the 17,20-lyase activity of cytochrome P450c17 (CYP17A1). Substitution of Ala, Gly, Cys, or Gln for these two glutamic acid residues abrogated all capacity to stimulate 17,20-lyase activity. Mutations E49D and E48D/E49D retained 23 and 38% of wild-type activity, respectively. Using the zero-length cross-linker ethyl-3-(3-dimethylaminopropyl)carbodiimide, we obtained cross-linked heterodimers of b5 and CYP17A1, wild-type, or mutations R347K and R358K. In sharp contrast, the b5 double mutation E48G/E49G did not form cross-linked complexes with wild-type CYP17A1. Mass spectrometric analysis of the CYP17A1-b5 complexes identified two cross-linked peptide pairs as follows: CYP17A1-WT: (84)EVLIKK(89)-b5: (53)EQAGGDATENFEDVGHSTDAR(73) and CYP17A1-R347K: (341)TPTISDKNR(349)-b5: (40)FLEEHPGGEEVLR(52). Using these two sites of interaction and Glu-48/Glu-49 in b5 as constraints, protein docking calculations based on the crystal structures of the two proteins yielded a structural model of the CYP17A1-b5 complex. The appositional surfaces include Lys-88, Arg-347, and Arg-358/Arg-449 of CYP17A1, which interact with Glu-61, Glu-42, and Glu-48/Glu-49 of b5, respectively. Our data reveal the structural basis of the electrostatic interactions between these two proteins, which is critical for 17,20-lyase activity and androgen biosynthesis.

Cross-linking mass spectrometry and mutagenesis confirm the functional importance of surface interactions between CYP3A4 and holo/apo cytochrome b(5)

Cytochrome b(5) (cyt b(5)) is one of the key components in the microsomal cytochrome P450 monooxygenase system. Consensus has not been reached about the underlying mechanism of cyt b(5) modulation of CYP catalysis. Both cyt b(5) and apo b(5) are reported to stimulate the activity of several P450 isoforms. In this study, the surface interactions of both holo and apo b(5) with CYP3A4 were investigated and compared for the first time. Chemical cross-linking coupled with mass spectrometric analysis was used to identify the potential electrostatic interactions between the protein surfaces. Subsequently, the models of interaction of holo/apo b(5) with CYP3A4 were built using the identified interacting sites as constraints. Both cyt b(5) and apo b(5) were predicted to bind to the same groove on CYP3A4 with close contacts to the B-B' loop of CYP3A4, a substrate recognition site. Mutagenesis studies further confirmed that the interacting sites on CYP3A4 (Lys96, Lys127, and Lys421) are functionally important. Mutation of these residues reduced or abolished cyt b(5) binding affinity. The critical role of Arg446 on CYP3A4 in binding to cyt b(5) and/or cytochrome P450 reductase was also discovered. The results indicated that electrostatic interactions on the interface of the two proteins are functionally important. The results indicate that apo b(5) can dock with CYP3A4 in a manner analogous to that of holo b(5), so electron transfer from cyt b(5) is not required for its effects.

Reaction of human cytochrome P450 3A4 with peroxynitrite: nitrotyrosine formation on the proximal side impairs its interaction with NADPH-cytochrome P450 reductase

The reaction of peroxynitrite (PN) with purified human cytochrome P450 3A4 (CYP3A4) resulted in the loss of the reduced-CO difference spectrum, but the absolute absorption spectrum of the heme was not significantly altered. The loss of 7-benzyloxy-4-(trifluoromethyl)coumarin (BFC) O-debenzylation activity of CYP3A4 was concentration-dependent with respect to PN, and the loss of BFC activity supported by NADPH-cytochrome P450 reductase (CPR) was much greater than that supported by tert-butyl hydroperoxide. Moreover, the PN-treated CYP3A4 exhibited a reduced-CO spectrum when reduced by CPR that was much smaller than when it was reduced by dithionite. These results suggest that modification of CYP3A4 by PN may impair its interaction with CPR, leading to the loss of catalytic activity. Tyrosine nitration, as measured by an increase in mass of 45 Da due to the addition of a nitro group, was used as a biomarker for protein modification by PN. PN-treated CYP3A4 was digested by trypsin and endoproteinase Glu C, and nitrotyrosine formation was then determined by using electrospray ionization-liquid chromatography-tandem mass spectrometry. Tyr residues 99, 307, 347, 430, and 432 were found to be nitrated. Using the GRAMM-X docking program, the structure for the CYP3A4-CPR complex shows that Tyr99, Tyr347, and Tyr430 are on the proximal side of CYP3A4 and are in close contact with three acidic residues in the FMN domain of CPR, suggesting that modification of one or more of these tyrosine residues by PN may influence CPR binding or the transfer of electrons to CYP3A4. Mutagenesis of Tyr430 to Phe or Val revealed that both the aromatic and the hydroxyl groups of Tyr are required for CPR-dependent catalytic activity and thus support the idea that the proximal side Tyr participates in the 3A4-CPR interaction. In conclusion, modification of tyrosine residues by PN and their subsequent identification can be used to enhance our knowledge of the structure/function relationships of the P450s with respect to the electron transfer steps, which are critical for P450 activity.

Three-dimensional structure of steroid 21-hydroxylase (cytochrome P450 21A2) with two substrates reveals locations of disease-associated variants

Steroid 21-hydroxylase (cytochrome P450 21A2, CYP21A2) deficiency accounts for ∼95% of individuals with congenital adrenal hyperplasia, a common autosomal recessive metabolic disorder of adrenal steroidogenesis. The effects of amino acid mutations on CYP21A2 activity lead to impairment of the synthesis of cortisol and aldosterone and the excessive production of androgens. In order to understand the structural and molecular basis of this group of diseases, the bovine CYP21A2 crystal structure complexed with the substrate 17-hydroxyprogesterone (17OHP) was determined to 3.0 Å resolution. An intriguing result from this structure is that there are two molecules of 17OHP bound to the enzyme, the distal one being located at the entrance of the substrate access channel and the proximal one bound in the active site. The substrate binding features locate the key substrate recognition residues not only around the heme but also along the substrate access channel. In addition, orientation of the skeleton of the proximal molecule is toward the interior of the enzyme away from the substrate access channel. The 17OHP complex of CYP21A2 provides a good relationship between the crystal structure, clinical data, and genetic mutants documented in the literature, thereby enhancing our understanding of congenital adrenal hyperplasia. In addition, the location of certain CYP21A2 mutations provides general understanding of structure/function relationships in P450s.

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.

Structure-based analysis of five novel disease-causing mutations in 21-hydroxylase-deficient patients

Congenital adrenal hyperplasia (CAH) due to 21-hydroxylase deficiency is the most frequent inborn error of metabolism, and accounts for 90–95% of CAH cases. The affected enzyme, P450C21, is encoded by the CYP21A2 gene, located together with a 98% nucleotide sequence identity CYP21A1P pseudogene, on chromosome 6p21.3. Even though most patients carry CYP21A1P-derived mutations, an increasing number of novel and rare mutations in disease causing alleles were found in the last years. In the present work, we describe five CYP21A2 novel mutations, p.R132C, p.149C, p.M283V, p.E431K and a frameshift g.2511_2512delGG, in four non-classical and one salt wasting patients from Argentina. All novel point mutations are located in CYP21 protein residues that are conserved throughout mammalian species, and none of them were found in control individuals. The putative pathogenic mechanisms of the novel variants were analyzed in silico. A three-dimensional CYP21 structure was generated by homology modeling and the protein design algorithm FoldX was used to calculate changes in stability of CYP21A2 protein. Our analysis revealed changes in protein stability or in the surface charge of the mutant enzymes, which could be related to the clinical manifestation found in patients.

Conformational changes of the NADPH-dependent cytochrome P450 reductase in the course of electron transfer to cytochromes P450

The NADPH-dependent cytochrome P450 reductase (CPR) is a key electron donor to eucaryotic cytochromes P450 (CYPs). CPR shuttles electrons from NADPH through the FAD and FMN-coenzymes into the iron of the prosthetic heme-group of the CYP. In the course of these electron transfer reactions, CPR undergoes large conformational changes. This mini-review discusses the new evidence provided for such conformational changes involving a combination of a "swinging" and "rotating" model and highlights the molecular mechanisms by which formation of these conformations are controlled and thereby enables CPR to serve as an effective electron transferring "nano-machine".

Application of the Amplex red/horseradish peroxidase assay to measure hydrogen peroxide generation by recombinant microsomal enzymes

The formation of reactive oxygen species by the cytochrome P450 monooxygenase system is thought to be due to autoxidation of NADPH-cytochrome P450 reductase and the nonproductive decay of oxygen-bound cytochrome P450 intermediates. To characterize this process in recombinant microsomal enzymes, we used a highly sensitive hydrogen peroxide assay based on Amplex red oxidation. This assay is 20 times more sensitive (LLD=5.0pmol/assay and LLQ=30pmol/assay) than the standard ferrous thiocyanate assay for detection of hydrogen peroxide. We found low, but detectable, spontaneous generation of hydrogen peroxide by recombinant human NADPH-cytochrome P450 reductase complexes (0.09nmol hydrogen peroxide/min/100Units of NADPH-cytochrome P450 reductase). Significantly higher rates of hydrogen peroxide production were observed when recombinant cytochrome P450 enzymes were coexpressed with NADPH-cytochrome P450 reductase (0.31nmol of hydrogen peroxide/min/100Units of NADPH-cytochrome P450 reductase). This was independent of the addition of any exogenous cytochrome P450 substrates. These data demonstrate that cytochrome P450s are a major source of hydrogen peroxide in the recombinant cytochrome P450 monooxygenase system. Moreover, substrate binding is not required for the cytochrome P450s to generate reactive oxygen species.

Human cytochrome P450 oxidoreductase deficiency caused by the Y181D mutation: molecular consequences and rescue of defect

Patients with congenital adrenal hyperplasia, exhibiting combined CYP17 and CYP21 deficiency, were shown by Arlt et al. (2004) to harbor a 541T-->G mutation in exon 5 of POR (encoding NADPH-cytochrome P450 reductase, CYPOR), which resulted in a Y181D substitution that obliterated electron transfer capacity. Using bacterial expression models, we examined catalytic and physical properties of the human CYPOR Y181D variant. As purified, Y181D lacked flavin mononucleotide (FMN) and NADPH-cytochrome c reductase (NCR) activity but retained normal flavin adenine dinucleotide binding and NADPH utilization. Titration of the purified protein with FMN restored 64 of wild-type (WT) NCR activity in Y181D with an activation constant of approximately 2 microM. As determined by FMN fluorescence quenching, Y181D had K(d)(FMN) = 7.3 microM. Biplasmid coexpression of CYPOR and CYP1A2, at the physiological ratio of approximately 1:10 in the engineered MK_1A2_POR Escherichia coli strain, showed the compromised capacity of Y181D to support CYP1A2-catalyzed metabolism of the procarcinogens 2-aminoanthracene, 2-amino-3-methylimidazo(4,5-f)quinoline, and 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone. Isolated MK1A2_POR membranes confirmed FMN stimulation of Y181D NCR activity with a 1.6 microM activation constant. CYP1A2 ethoxyresorufin-O-dealkylase activity of the MK1A2_POR(Y181D) membranes, undetectable in the absence of added FMN, increased to 37% of MK1A2_POR(WT) membranes with a 1.2 microM FMN activation constant. Therefore, we conclude that compromised FMN binding is the specific molecular defect causing POR deficiency in patients with Y181D mutation and that this defect, in large part, can be overcome in vitro by FMN addition.

Organization of cytochrome P450 enzymes involved in sex steroid synthesis: PROTEIN-PROTEIN INTERACTIONS IN LIPID MEMBRANES

Mounting evidence underscores the importance of protein-protein interactions in the functional regulation of drug-metabolizing P450s, but few studies have been conducted in membrane environments, and none have examined P450s catalyzing sex steroid synthesis. Here we report specific protein-protein interactions for full-length, human, wild type steroidogenic cytochrome P450 (P450, CYP) enzymes: 17alpha-hydroxylase/17,20-lyase (P450c17, CYP17) and aromatase (P450arom, CYP19), as well as their electron donor NADPH-cytochrome P450 oxidoreductase (CPR). Fluorescence resonance energy transfer (FRET)(3) in live cells, coupled with quartz crystal microbalance (QCM), and atomic force microscopy (AFM) studies on phosphatidyl choline +/- cholesterol (mammalian) biomimetic membranes were used to investigate steroidogenic P450 interactions. The FRET results in living cells demonstrated that both P450c17 and P450arom homodimerize but do not heterodimerize, although they each heterodimerize with CPR. The lack of heteroassociation between P450c17 and P450arom was confirmed by QCM, wherein neither enzyme bound a membrane saturated with the other. In contrast, the CPR bound readily to either P450c17- or P450arom-saturated surfaces. Interestingly, N-terminally modified P450arom was stably incorporated and gave similar results to the wild type, although saturation was achieved with much less protein, suggesting that the putative transmembrane domain is not required for membrane association but for orientation. In fact, all of the proteins were remarkably stable in the membrane, such that high resolution AFM images were obtained, further supporting the formation of P450c17, P450arom, and CPR homodimers and oligomers in lipid bilayers. This unique combination of in vivo and in vitro studies has provided strong evidence for homodimerization and perhaps some higher order interactions for both P450c17 and P450arom.

Steroid 17alpha-hydroxylase deficiency: functional characterization of four mutations (A174E, V178D, R440C, L465P) in the CYP17A1 gene

CONTEXT

Steroid 17alpha-hydroxylase (CYP17A1, alias P450c17) deficiency (17OHD) is a rare form of congenital adrenal hyperplasia. The CYP17A1 enzyme catalyzes two distinct reactions, 17alpha-hydroxylase and 17,20-lyase activities.

OBJECTIVE

The aim of the study was to analyze the structural and functional consequences of three novel (A174E, V178D, and L465P) and one previously reported (R440C) CYP17A1 mutation found in three patients clinically and biochemically presenting with 17OHD.

PATIENTS AND METHODS

Two patients suffering from 46,XY disordered sex development presented at ages 5.5 and 8.8 yr, respectively, with tall stature and hypertension. Mutation analysis revealed compound heterozygous CYP17A1 mutations (A174E/K388X; V178D/R440C). The third patient (46,XX) presented with primary amenorrhea and hypertension at age 15 yr. She was homozygous for the novel L465P mutation. Functional studies employing a yeast microsomal expression system compared wild-type and mutant CYP17A1 both with regard to 17alpha-hydroxylase and 17,20-lyase activity. Mutants were examined in a computational three-dimensional model of the CYP17A1 protein.

RESULTS

The activity assays showed that all three mutants retain only 0-7% of both 17alpha-hydroxylase and 17,20-lyase activity relative to CYP17A1 wild-type activity, corresponding to the in vivo situation. Enzyme kinetic studies proved the impairment of both reactions, respectively. Computer-based three-dimensional model analysis of CYP17A1 using CYP2B4 as template showed that three of the mutations had no direct effect on the active center, whereas one affects the heme coordination.

CONCLUSION

The functional studies revealed that the described missense mutations result in severe 17OHD. Our data are important to predict the phenotypic expressions and provide important information for patient management and genetic counseling.

Functional and structural consequences of a novel point mutation in the CYP21A2 gene causing congenital adrenal hyperplasia: potential relevance of helix C for P450 oxidoreductase-21-hydroxylase interaction

BACKGROUND

Congenital adrenal hyperplasia is caused by insufficient adrenal steroid biosynthesis due to impaired steroidogenic enzymes. The majority of patients suffer from deficiency of 21-hydroxylase (CYP21) coded by the CYP21A2 gene.

OBJECTIVE

Our objective was to study the functional and structural consequences of the novel CYP21A2 missense mutation c.364A > C (K121Q) detected in a female patient with nonclassical 21-hydroxylase deficiency. The patient was compound heterozygous for the novel K121Q mutation and the mild P453S mutation.

RESULTS

In vitro expression analysis of the mutant K121Q enzyme in transiently transfected COS-7 cells revealed reduced CYP21 activity of 14.0 +/- 5% for the conversion of 17-hydroxyprogesterone and 19.5 +/- 4% for the conversion of progesterone. K121 is located on helix C in the CYP21 protein, which is part of the heme coordinating system. In addition, helix C is involved in the interaction with the electron-providing enzyme P450 oxidoreductase. Protein modeling revealed that the substitution of glutamine for the basic amino acid lysine introduces an electrostatic change on the surface of CYP21 and may additionally change heme coordination. We hypothesize that the electron flux between P450 oxidoreductase and CYP21 is impaired and, moreover, that substrate affinity is altered due to heme dislocation with K121Q.

CONCLUSION

Both the interaction of P450 oxidoreductase and CYP21 as well as heme coordination are likely to be disturbed due to the K121Q mutation. Our data exemplify how the combination of in vitro expression and structural protein analysis provide novel insights into molecular mechanisms of reduced CYP21 activity, eventually explaining the patient's phenotype.

Impairment of human CYP1A2-mediated xenobiotic metabolism by Antley-Bixler syndrome variants of cytochrome P450 oxidoreductase

Y459H and V492E mutations of cytochrome P450 reductase (CYPOR) cause Antley-Bixler syndrome due to diminished binding of the FAD cofactor. To address whether these mutations impaired the interaction with drug-metabolizing CYPs, a bacterial model of human liver expression of CYP1A2 and CYPOR was implemented. Four models were generated: POR(null), POR(wt), POR(YH), and POR(VE), for which equivalent CYP1A2 and CYPOR levels were confirmed, except for POR(null), not containing any CYPOR. The mutant CYPORs were unable to catalyze cytochrome c and MTT reduction, and were unable to support EROD and MROD activities. Activity was restored by the addition of FAD, with V492E having a higher apparent FAD affinity than Y459H. The CYP1A2-activated procarcinogens, 2-aminoanthracene, 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone, and 2-amino-3-methylimidazo(4,5-f)quinoline, were significantly less mutagenic in POR(YH) and POR(VE) models than in POR(wt), indicating that CYP1A2, and likely other drug-metabolizing CYPs, are impaired by ABS-related POR mutations as observed in the steroidogenic CYPs.

Structural analysis of CYP2R1 in complex with vitamin D3

The activation of vitamin D to its hormonal form is mediated by cytochrome P450 enzymes. CYP2R1 catalyzes the initial step converting vitamin D into 25-hydroxyvitamin D. A CYP2R1 gene mutation causes an inherited form of rickets due to 25-hydroxylase deficiency. To understand the narrow substrate specificity of CYP2R1 we obtained the hemeprotein in a highly purified state, confirmed the enzyme as a vitamin D 25-hydroxylase, and solved the crystal structure of CYP2R1 in complex with vitamin D3. The CYP2R1 structure adopts a closed conformation with the substrate access channel being covered by the ordered B'-helix and slightly opened to the surface, which defines the substrate entrance point. The active site is lined by conserved, mostly hydrophobic residues. Vitamin D3 is bound in an elongated conformation with the aliphatic side-chain pointing toward the heme. The structure reveals the secosteroid binding mode in an extended active site and allows rationalization of the molecular basis of the inherited rickets associated with CYP2R1.

Global analysis of protein-protein interactions reveals multiple CYP2E1-reductase complexes

Although a single binary functional complex between cytochrome P450 (P450 or CYP for a specific isoform) and cytochrome P450 reductase (CPR) has been generally accepted in the literature, this simple model failed to explain the experimentally observed catalytic activity of recombinant CYP2E1 in dependence on the total concentration of the added CPR-K56Q mutant. Our rejection of the simplest 1:1 binding model was based on two independent lines of experimental evidence. First, under the assumption of the 1:1 binding model, separate analyses of titration curves obtained while varying either P450 or CPR concentrations individually produced contradictory results. Second, an asymmetric Job plot suggested the existence of higher order molecular complexes. To identify the most probable complexation mechanism, we generated a comprehensive data set where the concentrations of both P450 and P450 were varied simultaneously, rather than one at a time. The resulting two-dimensional data were globally fit to 32 candidate mechanistic models, involving the formation of binary, ternary, and quaternary P450.CPR complexes, in the absence or presence or P450 and CPR homodimers. Of the 32 candidate models (mechanisms), two models were approximately equally successful in explaining our experimental data. The first plausible model involves the binary complex P450.CPR, the quaternary complex (P450)2.(CPR)2, and the homodimer (P450)2. The second plausible model additionally involves a weakly bound ternary complex (P450)2.CPR. Importantly, only the binary complex P450.CPR seems catalytically active in either of the two most probable mechanisms.

Molecular mechanism of phase I and phase II drug-metabolizing enzymes: implications for detoxification

Enzymes that catalyze the biotransformation of drugs and xenobiotics are generally referred to as drug-metabolizing enzymes (DMEs). DMEs can be classified into two main groups: oxidative or conjugative. The NADPH-cytochrome P450 reductase (P450R)/cytochrome P450 (P450) electron transfer systems are oxidative enzymes that mediate phase I reactions, whereas the UDP-glucuronosyltransferases (UGTs) are conjugative enzymes that mediate phase II enzymes. Both enzyme systems are localized to the endoplasmic reticulum (ER) where a number of drugs are sequentially metabolized. DMEs, including P450s and UGTs, generally have a highly plastic active site that can accommodate a wide variety of substrates. The P450 and UGT genes constitute a supergene family, in which UGT proteins are encoded by distinct genes and a complex gene. Both the P450 and UGT genes have evolved to diversify their functions. This chapter reviews advances in understanding the structure and function of the P450R/P450 and UGT enzyme systems. In particular, the coordinate biotransformation of xenobiotics by phase I and II enzymes in the ER membrane is examined.

Allelic Variation in the Depressaria pastinacella CYP6AB3 Protein Enhances Metabolism of Plant Allelochemicals by Altering a Proximal Surface Residue and Potential Interactions with Cytochrome P450 Reductase

CYP6AB3v1, a cytochrome P450 monooxygenase in Depressaria pastinacella (parsnip webworm), is highly specialized for metabolizing imperatorin, a toxic furanocoumarin in the apiaceous host plants of this insect. Cloning and heterologous expression of CYP6AB3v2, an allelic variant identified in D. pastinacella, reveals that it metabolizes imperatorin at a rate (V(max) of 10.02 pmol/min/pmol of cytochrome P450 monooxygenase (P450)) significantly higher than CYP6AB3v1 (V(max) of 2.41 pmol/min/pmol) when supplemented with even low levels of cytochrome P450 reductase. Comparisons of the NADPH consumption rates for these variants indicate that CYP6AB3v2 utilizes this electron source at a faster rate than does CYP6AB3v1. Molecular modeling of the five amino acid differences between these variants and their potential interactions with P450 reductase suggests that replacement of Val(92) on the proximal face of CYP6AB3v1 with Ala(92) in CYP6AB3v2 affects interactions with P450 reductase so as to enhance its catalytic activity. Allelic variation at this locus potentially allows D. pastinacella to adapt to both intraspecific and interspecific variation in imperatorin concentrations in its host plants.

Adrenodoxin supports reactions catalyzed by microsomal steroidogenic cytochrome P450s

The interaction of adrenodoxin (Adx) and NADPH cytochrome P450 reductase (CPR) with human microsomal steroidogenic cytochrome P450s was studied. It is found that Adx, mitochondrial electron transfer protein, is able to support reactions catalyzed by human microsomal P450s: full length CYP17, truncated CYP17, and truncated CYP21. CPR, but not Adx, supports activity of truncated CYP19. Truncated and the full length CYP17s show distinct preference for electron donor proteins. Truncated CYP17 has higher activity with Adx compared to CPR. The alteration in preference to electron donor does not change product profile for truncated enzymes. The electrostatic contacts play a major role in the interaction of truncated CYP17 with either CPR or Adx. Similarly electrostatic contacts are predominant in the interaction of full length CYP17 with Adx. We speculate that Adx might serve as an alternative electron donor for CYP17 at the conditions of CPR deficiency in human.

Mapping electrostatic potential of a protein on its hydrophobic surface: Implications for crystallization of Cytochrome P450scc

Calculation and combined visualization of electrostatic and hydrophobic properties of Cytochrome P450scc based on two very different homology models allowed to identify extensive hydrophobic patches with neutral electrostatic potential and mutations removing such patches and thus expecting to facilitate crystallization of Cytochrome P450scc, especially for the nanotemplate crystallization method. Implications are discussed for optimizing crystallization and other aspects of protein surface properties and protein recognition.

Molecular model of human CYP21 based on mammalian CYP2C5: structural features correlate with clinical severity of mutations causing congenital adrenal hyperplasia

Enhanced understanding of structure-function relationships of human 21-hydroxylase, CYP21, is required to better understand the molecular causes of congenital adrenal hyperplasia. To this end, a structural model of human CYP21 was calculated based on the crystal structure of rabbit CYP2C5. All but two known allelic variants of missense type, a total of 60 disease-causing mutations and six normal variants, were analyzed using this model. A structural explanation for the corresponding phenotype was found for all but two mutants for which available clinical data are also discrepant with in vitro enzyme activity. Calculations of protein stability of modeled mutants were found to correlate inversely with the corresponding clinical severity. Putative structurally important residues were identified to be involved in heme and substrate binding, redox partner interaction, and enzyme catalysis using docking calculations and analysis of structurally determined homologous cytochrome P450s (CYPs). Functional and structural consequences of seven novel mutations, V139E, C147R, R233G, T295N, L308F, R366C, and M473I, detected in Scandinavian patients with suspected congenital adrenal hyperplasia of different severity, were predicted using molecular modeling. Structural features deduced from the models are in good correlation with clinical severity of CYP21 mutants, which shows the applicability of a modeling approach in assessment of new CYP21 mutations.

Identification of the interactions between cytochrome P450 2E1 and cytochrome b5 by mass spectrometry and site-directed mutagenesis

The reaction cycles of cytochrome P450s (P450) require input of two electrons. Electrostatic interactions are considered important driving forces in the association of P450s with their redox partners, which in turn facilitates the transfer of the two electrons. In this study, the cross-linking reagent, 1-ethyl-3-[3-dimethylaminopropyl]carbodiimide hydrochloride (EDC), was used to covalently link cytochrome P450 2E1 (CYP2E1) with cytochrome b(5) (b(5)) through the formation of specific amide bonds between complementary charged residue pairs. Cross-linked peptides in the resulting protein complex were distinguished from non-cross-linked peptides using an (18)O-labeling method on the basis that cross-linked peptides incorporate twice as many (18)O atoms as non-cross-linked peptides during proteolysis conducted in (18)O-water. Subsequent tandem mass spectrometric (MS/MS) analysis of the selected cross-linked peptide candidates led to the identification of two intermolecular cross-links, Lys(428)(CYP2E1)-Asp(53)(b(5)) and Lys(434)(CYP2E1)-Glu(56)(b(5)), which provides the first direct experimental evidence for the interacting orientations of a microsomal P450 and its redox partner. The biological importance of the two ion pairs for the CYP2E1-b(5) interaction, and the stimulatory effect of b(5), was confirmed by site-directed mutagenesis. Based on the characterized cross-links, a CYP2E1-b(5) complex model was constructed, leading to improved insights into the protein interaction. The described method is potentially useful for mapping the interactions of various P450 isoforms and their redox partners, because the method is relatively rapid and sensitive, and is capable of suggesting not only protein interacting regions, but also interacting orientations.

Functional studies of two novel and two rare mutations in the 21-hydroxylase gene

Congenital adrenal hyperplasia (CAH) is most commonly due to 21-hydroxylase deficiency and presents with a wide spectrum of clinical manifestations, from prenatal virilization and salt-wasting in the neonatal period to precocious pubarche and late-onset hyperandrogenic symptoms during adulthood. A limited number of mutations account for the majority of all mutated alleles, but a growing number of rare mutations are responsible for the disease in some patients. By sequence analysis of the CYP21A2 gene, we identified two novel (I171N and L446P) and two rare (R341P and R426H) mutations in seven Italian patients with CAH. One of the patients was diagnosed with mild non-classical CAH and was found to be a compound heterozygote (I171N/V281L), while all other patients showed severe phenotypes with latent or manifest salt-wasting. The residual activities measured after expression of the four mutant enzymes in COS-1 cells were all below 1% towards both natural substrates (17-OH-progesterone and progesterone) compared with the wild-type protein. All four mutations are, thus, associated with severe enzyme deficiency and are predicted to cause classic CAH if found in trans with other mutations causing severe enzyme deficiency.

Thermal inactivation of the reductase domain of cytochrome P450 BM3

Although the reductase domain of cytochrome P450 BM3 (BMR) catalyzes the reduction of cytochrome c and 2,6-dichlorophenolindophenol, we observed a catalytically independent loss of activity. By varying the incubation time for the enzyme prior to reaction initiation, we measured an inactivation rate of 0.22 min(-1). We hypothesized that either an active BMR dimer dissociates to an inactive monomer or BMR undergoes denaturation. We were not able to trap or destabilize a dimer, and BMR inactivation proved to be irreversible. Addition of excess FMN only slightly decreased the rate of inactivation from 0.22 to 0.13 min(-1), indicating inactivation likely does not reflect loss of flavin. When inactivation rates as a function of temperature were fit to the Arrhenius equation, the energy required to inactivate BMR was 9.9 kcal mol(-1)--equivalent to a few hydrogen bonds. The potential instability of BMR under certain conditions raises concerns for the use of BMR as a model or surrogate P450 reductase in other systems.