Structural basis for human NADPH-cytochrome P450 oxidoreductase deficiency

Mitochondria-endoplasmic reticulum crosstalk in apoptosis: The interactions of cytochrome c with monooxygenase and its reductase

The crosstalk between endoplasmic reticulum and mitochondria is of significance in apoptosis, in which cytochrome b5 (Cyt b5) is thought to be a major target for cytochrome c (Cyt c) upon its release from the mitochondria. In the absence of Cyt b5, the role of interactions of Cyt c with CYP-dependent monooxygenase system in apoptotic regulation was explored in this study. NADPH-dependent and Cyt c-induced formation of reactive oxygen species (ROS) and NADPH-independent Cyt c unfolding were revealed. With the aid of a CPR inhibitor and CYP antibodies, the interactions among Cyt c, cytochrome P450 reductase (CPR) and cytochrome P450 (CYP) are evidenced, which are found crucial for monooxygenase-derived ROS formation. The underlying structural basis of Cyt c-CYP complex was unveiled by molecular dynamics simulations. This study provides novel insights into how Cyt c regulates ROS formation through the interactions with CPR and CYP, and is implicated for a deeper understanding of the regulation mechanism in the mitochondria-endoplasmic reticulum apoptotic pathway.

Sterol 14-alpha demethylase (CYP51) activity in Leishmania donovani is likely dependent upon cytochrome P450 reductase 1

Liposomal amphotericin B is an important frontline drug for the treatment of visceral leishmaniasis, a neglected disease of poverty. The mechanism of action of amphotericin B (AmB) is thought to involve interaction with ergosterol and other ergostane sterols, resulting in disruption of the integrity and key functions of the plasma membrane. Emergence of clinically refractory isolates of Leishmania donovani and L. infantum is an ongoing issue and knowledge of potential resistance mechanisms can help to alleviate this problem. Here we report the characterisation of four independently selected L. donovani clones that are resistant to AmB. Whole genome sequencing revealed that in three of the moderately resistant clones, resistance was due solely to the deletion of a gene encoding C24-sterol methyltransferase (SMT1). The fourth, hyper-resistant resistant clone (>60-fold) was found to have a 24 bp deletion in both alleles of a gene encoding a putative cytochrome P450 reductase (P450R1). Metabolic profiling indicated these parasites were virtually devoid of ergosterol (0.2% versus 18% of total sterols in wild-type) and had a marked accumulation of 14-methylfecosterol (75% versus 0.1% of total sterols in wild-type) and other 14-alpha methylcholestanes. These are substrates for sterol 14-alpha demethylase (CYP51) suggesting that this enzyme may be a bona fide P450R specifically involved in electron transfer from NADPH to CYP51 during catalysis. Deletion of P450R1 in wild-type cells phenocopied the metabolic changes observed in our AmB hyper-resistant clone as well as in CYP51 nulls. Likewise, addition of a wild type P450R1 gene restored sterol profiles to wild type. Our studies indicate that P450R1 is essential for L. donovani amastigote viability, thus loss of this gene is unlikely to be a driver of clinical resistance. Nevertheless, investigating the mechanisms underpinning AmB resistance in these cells provided insights that refine our understanding of the L. donovani sterol biosynthetic pathway.

Computational identification and analysis of deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene: a structural and functional impact

Cytochrome P450 oxidoreductase (POR) protein is essential for steroidogenesis, and POR gene mutations are frequently associated with P450 Oxidoreductase Deficiency (PORD), a disorder of hormone production. To our knowledge, no previous attempt has been made to identify and analyze the deleterious/pathogenic non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene through an extensive computational approach. Computational algorithms and tools were employed to identify, characterize, and validate the pathogenic SNPs associated with certain diseases. To begin with, all the high-confidence SNPs were collected, and their structural and functional impacts on the protein structures were explored. The results of various in silico analyses affirm that the A287P and R457H variants of POR could destabilize the interactions between the amino acids and the hydrogen bond networks, resulting in functional deviations of POR. The literature study further confirms that the pathogenic mutations (A287P and R457H) are associated with the onset of PORD. Molecular dynamics simulations (MDS) and essential dynamics (ED) studies characterized the structural consequences of prioritized deleterious mutations, representing the structural destabilization that might disrupt POR biological function. The identified deleterious mutations at the cofactor's binding domains might interfere with the essential interactions between the protein and cofactors, thus inhibiting POR catalytic activity. The consolidated insights from the computational analyses can be used to predict potential deleterious mutants and understand the disease's pathological basis and the molecular mechanism of drug metabolism for the application of personalized medication. HIGHLIGHTSNADPH cytochrome P450 oxidoreductase (POR) mutations are associated with a broad spectrum of human diseasesIdentified and analyzed the most deleterious nsSNPs of POR through the sequence and structure-based prediction toolsInvestigated the structural and functional impacts of the most significant mutations (A287P and R457H) associated with PORDMolecular dynamics and PCA-based FEL analysis were utilized to probe the mutation-induced structural alterations in PORCommunicated by Ramaswamy H. Sarma.

Exploring Novel Variants of the Cytochrome P450 Reductase Gene (POR) from the Genome Aggregation Database by Integrating Bioinformatic Tools and Functional Assays

Cytochrome P450 oxidoreductase (POR) is an essential redox partner for steroid and drug-metabolizing cytochromes P450 located in the endoplasmic reticulum. Mutations in POR lead to metabolic disorders, including congenital adrenal hyperplasia, and affect the metabolism of steroids, drugs, and xenobiotics. In this study, we examined approximately 450 missense variants of the POR gene listed in the Genome Aggregation Database (gnomAD) using eleven different in silico prediction tools. We found that 64 novel variants were consistently predicted to be disease-causing by most tools. To validate our findings, we conducted a population analysis and selected two variations in POR for further investigation. The human POR wild type and the R268W and L577P variants were expressed in bacteria and subjected to enzyme kinetic assays using a model substrate. We also examined the activities of several cytochrome P450 proteins in the presence of POR (WT or variants) by combining P450 and reductase proteins in liposomes. We observed a decrease in enzymatic activities (ranging from 35% to 85%) of key drug-metabolizing enzymes, supported by POR variants R288W and L577P compared to WT-POR. These results validate our approach of curating a vast amount of data from genome projects and provide an updated and reliable reference for diagnosing POR deficiency.

Biological Role of Pazopanib and Sunitinib Aldehyde Derivatives in Drug-Induced Liver Injury

Tyrosine kinase inhibitors pazopanib and sunitinib are both used to treat advanced renal cell carcinoma but expose patients to an increased risk of hepatotoxicity. We have previously identified two aldehyde derivatives for pazopanib and sunitinib (P-CHO and S-CHO, respectively) in liver microsomes. In this study, we aimed to decipher their role in hepatotoxicity by treating HepG2 and HepaRG hepatic cell lines with these derivatives and evaluating cell viability, mitochondrial dysfunction, and oxidative stress accumulation. Additionally, plasma concentrations of P-CHO were assessed in a cohort of patients treated with pazopanib. Results showed that S-CHO slightly decreased the viability of HepG2, but to a lesser extent than sunitinib, and affected the maximal respiratory capacity of the mitochondrial chain. P-CHO decreased viability and ATP production in HepG2. Traces of P-CHO were detected in the plasma of patients treated with pazopanib. Overall, these results showed that P-CHO and S-CHO affect hepatocyte integrity and could be involved in the pazopanib and sunitinib hepatotoxicity.

Characterization of POR haplotype distribution in African populations and comparison with other global populations

Background & aim: POR is an enzyme that mediates electron transfer to enable the drug-metabolizing activity of CYP450 proteins. However, POR has been understudied in pharmacogenomics despite this vital role. This study aimed to characterize the genetic variation in POR across African populations and to compare the star allele (haplotype) distribution with that in other global populations. Materials & methods: POR star alleles were called from whole-genome sequencing data using the StellarPGx pipeline. Results: In addition to the common POR*1 and *28 (defined by rs1057868), five novel rare haplotypes were computationally inferred. No significant frequency differences were observed among the majority of African populations. However, POR*28 was observed at a higher frequency in individuals of non-African ancestry. Conclusion: This study highlights the distribution of POR alleles in Africa and across global populations with a view toward informing future precision medicine implementation.

NADPH–Cytochrome P450 Reductase Mediates the Fatty Acid Desaturation of ω3 and ω6 Desaturases from Mortierella alpina

Fatty acid desaturases play an important role in maintaining the appropriate structure and function of biological membranes. The biochemical characterization of integral membrane desaturases, particularly ω3 and ω6 desaturases, has been limited by technical difficulties relating to the acquisition of large quantities of purified proteins, and by the fact that functional activities of these proteins were only tested in an NADH-initiated reaction system. The main aim of this study was to reconstitute an NADPH-dependent reaction system in vitro and investigate the kinetic properties of Mortierella alpina ω3 and ω6 desaturases in this system. After expression and purification of the soluble catalytic domain of NADPH–cytochrome P450 reductase, the NADPH-dependent fatty acid desaturation was reconstituted for the first time in a system containing NADPH, NADPH–cytochrome P450 reductase, cytochrome b5, M. alpina ω3 and ω6 desaturase and detergent. In this system, the maximum activity of ω3 and ω6 desaturase was 213.4 ± 9.0 nmol min⁻¹ mg⁻¹ and 10.0 ± 0.5 nmol min⁻¹ mg⁻¹, respectively. The highest k_cat/K_m value of ω3 and ω6 desaturase was 0.41 µM⁻¹ min⁻¹ and 0.09 µM⁻¹ min⁻¹ when using linoleoyl CoA (18:2 ω6) and oleoyl CoA (18:1 ω9) as substrates, respectively. M. alpina ω3 and ω6 desaturases were capable of using NADPH as reductant when mediated by NADPH–cytochrome P450 reductase; although, their efficiency is distinguishable from NADH-dependent desaturation. These results provide insights into the mechanisms underlying ω3 and ω6 fatty acid desaturation and may facilitate the production of important fatty acids in M. alpina.

Neutron scattering maps the higher-order assembly of NADPH-dependent assimilatory sulfite reductase

Precursor molecules for biomass incorporation must be imported into cells and made available to the molecular machines that build the cell. Sulfur-containing macromolecules require that sulfur be in its S^2- oxidation state before assimilation into amino acids, cofactors, and vitamins that are essential to organisms throughout the biosphere. In α-proteobacteria, NADPH-dependent assimilatory sulfite reductase (SiR) performs the final six-electron reduction of sulfur. SiR is a dodecameric oxidoreductase composed of an octameric flavoprotein reductase (SiRFP) and four hemoprotein metalloenzyme oxidases (SiRHP). SiR performs the electron transfer reduction reaction to produce sulfide from sulfite through coordinated domain movements and subunit interactions without release of partially reduced intermediates. Efforts to understand the electron transfer mechanism responsible for SiR's efficiency are confounded by structural heterogeneity arising from intrinsically disordered regions throughout its complex, including the flexible linker joining SiRFP's flavin-binding domains. As a result, high-resolution structures of SiR dodecamer and its subcomplexes are unknown, leaving a gap in the fundamental understanding of how SiR performs this uniquely large-volume electron transfer reaction. Here, we use deuterium labeling, in vitro reconstitution, analytical ultracentrifugation (AUC), small-angle neutron scattering (SANS), and neutron contrast variation (NCV) to observe the relative subunit positions within SiR's higher-order assembly. AUC and SANS reveal SiR to be a flexible dodecamer and confirm the mismatched SiRFP and SiRHP subunit stoichiometry. NCV shows that the complex is asymmetric, with SiRHP on the periphery of the complex and the centers of mass between SiRFP and SiRHP components over 100 Å apart. SiRFP undergoes compaction upon assembly into SiR's dodecamer and SiRHP adopts multiple positions in the complex. The resulting map of SiR's higher-order structure supports a cis/trans mechanism for electron transfer between domains of reductase subunits as well as between tightly-bound or transiently-interacting reductase and oxidase subunits.

Conformational Rearrangements in the Redox Cycling of NADPH-Cytochrome P450 Reductase from Sorghum bicolor Explored with FRET and Pressure-Perturbation Spectroscopy

Simple Summary

NADPH-cytochrome P450 reductase (CPR) enzymes are known to undergo an ample conformational transition between the closed and open states in the process of their redox cycling. To explore the conformational landscape of CPR from the potential biofuel crop Sorghum bicolor (SbCPR), we incorporated a FRET donor/acceptor pair into the enzyme and employed rapid scanning stop-flow and pressure perturbation spectroscopy to characterize the equilibrium between its open and closed states at different stages of the redox cycle. Our results suggest the presence of several open conformational sub-states differing in the system volume change associated with the opening transition (ΔV⁰). Although the closed conformation always predominates in the conformational landscape, the population of the open conformations increases by order of magnitude upon the two-electron reduction and the formation of the disemiquinone state of the enzyme. In addition to elucidating the functional choreography of plant CPRs, our study demonstrates the high exploratory potential of a combination of the pressure-perturbation approach with the FRET-based monitoring of protein conformational transitions.

Abstract

NADPH-cytochrome P450 reductase (CPR) from Sorghum bicolor (SbCPR) serves as an electron donor for cytochrome P450 essential for monolignol and lignin production in this biofuel crop. The CPR enzymes undergo an ample conformational transition between the closed and open states in their functioning. This transition is triggered by electron transfer between the FAD and FMN and provides access of the partner protein to the electron-donating FMN domain. To characterize the electron transfer mechanisms in the monolignol biosynthetic pathway better, we explore the conformational transitions in SbCPR with rapid scanning stop-flow and pressure-perturbation spectroscopy. We used FRET between a pair of donor and acceptor probes incorporated into the FAD and FMN domains of SbCPR, respectively, to characterize the equilibrium between the open and closed states and explore its modulation in connection with the redox state of the enzyme. We demonstrate that, although the closed conformation always predominates in the conformational landscape, the population of open state increases by order of magnitude upon the formation of the disemiquinone state. Our results are consistent with several open conformation sub-states differing in the volume change (ΔV⁰) of the opening transition. While the ΔV⁰ characteristic of the oxidized enzyme is as large as −88 mL/mol, the interaction of the enzyme with the nucleotide cofactor and the formation of the double-semiquinone state of CPR decrease this value to −34 and −18 mL/mol, respectively. This observation suggests that the interdomain electron transfer in CPR increases protein hydration, while promoting more open conformation. In addition to elucidating the functional choreography of plant CPRs, our study demonstrates the high exploratory potential of a combination of the pressure-perturbation approach with the FRET-based monitoring of protein conformational transitions.

Functional and structural insight into the flexibility of cytochrome P450 reductases from Sorghum bicolor and its implications for lignin composition

Plant NADPH-dependent cytochrome P450 reductase (CPR) is a multidomain enzyme that donates electrons for hydroxylation reactions catalyzed by class II cytochrome P450 monooxygenases involved in the synthesis of many primary and secondary metabolites. These P450 enzymes include trans-cinnamate-4-hydroxylase, p-coumarate-3′-hydroxylase, and ferulate-5-hydroxylase involved in monolignol biosynthesis. Because of its role in monolignol biosynthesis, alterations in CPR activity could change the composition and overall output of lignin. Therefore, to understand the structure and function of three CPR subunits from sorghum, recombinant subunits SbCPR2a, SbCPR2b, and SbCPR2c were subjected to X-ray crystallography and kinetic assays. Steady-state kinetic analyses demonstrated that all three CPR subunits supported the oxidation reactions catalyzed by SbC4H1 (CYP73A33) and SbC3′H (CYP98A1). Furthermore, comparing the SbCPR2b structure with the well-investigated CPRs from mammals enabled us to identify critical residues of functional importance and suggested that the plant flavin mononucleotide–binding domain might be more flexible than mammalian homologs. In addition, the elucidated structure of SbCPR2b included the first observation of NADP⁺ in a native CPR. Overall, we conclude that the connecting domain of SbCPR2, especially its hinge region, could serve as a target to alter biomass composition in bioenergy and forage sorghums through protein engineering.

Molecular mechanism of interactions between ACAD9 and binding partners in mitochondrial respiratory complex I assembly

The dual function protein ACAD9 catalyzes α,β-dehydrogenation of fatty acyl-CoA thioesters in fatty acid β-oxidation and is an essential chaperone for mitochondrial respiratory complex I (CI) assembly. ACAD9, ECSIT, and NDUFAF1 interact to form the core mitochondrial CI assembly complex. Current studies examine the molecular mechanism of ACAD9/ECSIT/NDUFAF1interactions. ACAD9 binds to the carboxy-terminal half and NDUFAF1 to the amino-terminal half of ECSIT. Binary complexes are unstable and aggregate easily, while the ACAD9/ECSIT/NDUFAF1 ternary complex is soluble and highly stable. Molecular modeling and small-angle X-ray scattering studies identified intra-complex interaction sites and binding sites for other assembly factors. Binding of ECSIT at the ETF binding site in the amino-terminal domain of ACAD9 is consistent with observed loss of FAD and enzymatic activity and demonstrates that the two functions of ACAD9 are mutually exclusive. Mapping of 42 known pathogenic mutations onto the homology-modeled ACAD9 structure provides structural insights into pathomechanisms of CI deficiency.

In Silico Structural, Functional, and Phylogenetic Analysis of Cytochrome (CYPD) Protein Family

Cytochrome (CYP) enzymes catalyze the metabolic reactions of endogenous and exogenous compounds. The superfamily of enzymes is found across many organisms, regardless of type, except for plants. Information was gathered about CYP2D enzymes through protein sequences of humans and other organisms. The secondary structure was predicted using the SOPMA. The structural and functional study of human CYP2D was conducted using ProtParam, SOPMA, Predotar 1.03, SignalP, TMHMM 2.0, and ExPASy. Most animals shared five central motifs according to motif analysis results. The tertiary structure of human CYP2D, as well as other animal species, was predicted by Phyre2. Human CYP2D proteins are heavily conserved across organisms, according to the findings. This indicates that they are descended from a single ancestor. They calculate the ratio of alpha-helices to extended strands to beta sheets to random coils. Most of the enzymes are alpha-helix, but small amounts of the random coil were also found. The data were obtained to provide us with a better understanding of mammalian proteins' functions and evolutionary relationships.

Non-classic cytochrome P450 oxidoreductase deficiency strongly linked with menstrual cycle disorders and female infertility as primary manifestations

STUDY QUESTION

Can cytochrome P450 oxidoreductase deficiency (PORD) be revealed in adult women with menstrual disorders and/or infertility?

SUMMARY ANSWER

PORD was biologically and genetically confirmed in five adult women with chronically elevated serum progesterone (P) who were referred for oligo-/amenorrhea and/or infertility.

WHAT IS KNOWN ALREADY

PORD is an autosomal recessive disease typically diagnosed in neonates and children with ambiguous genitalia and/or skeletal abnormalities. It is responsible for the decreased activity of several P450 enzymes, including CYP21A2, CYP17A1 and CYP19A1, that are involved in adrenal and/or gonadal steroidogenesis. Little is known about the optimal way to investigate and treat patients with adult-onset PORD.

STUDY DESIGN, SIZE, DURATION

In this series, we report five adult females who were evaluated in three tertiary endocrine reproductive departments between March 2015 and September 2018.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Five women aged 19-38 years were referred for unexplained oligo-/amenorrhea and/or infertility. Genetic testing excluded 21-hydroxylase deficiency (21OH-D), initially suspected due to the increased 17-hydroxyprogesterone (17-OHP) levels. Extensive phenotyping, steroid profiling by mass spectrometry, pelvic imaging and next-generation sequencing of 84 genes involved in gonadal and adrenal disorders were performed in all patients. IVF followed by frozen embryo transfer (ET) under glucocorticoid suppression therapy was performed for two patients.

MAIN RESULTS AND THE ROLE OF CHANCE

All patients had oligomenorrhea or amenorrhea. None had hyperandrogenism. Low-normal serum estradiol (E2) and testosterone levels contrasted with chronically increased serum P and 17-OHP levels, which further increased after adrenocorticotrophic hormone (ACTH) administration. Despite excessive P, 17OH-P and 21-deoxycortisol rise after ACTH stimulation suggesting non-classic 21OH-D, CYP21A2 sequencing did not support this hypothesis. Basal serum cortisol levels were low to normal, with inadequate response to ACTH in some women, suggesting partial adrenal insufficiency. All patients harbored rare biallelic POR mutations classified as pathogenic or likely pathogenic according to the American College of Medical Genetics and Genomics standards. Pelvic imaging revealed bilateral ovarian macrocysts in all women. IVF was performed for two women after retrieval of a normal oocyte number despite very low E2 levels during ovarian stimulation. Frozen ET under glucocorticoid suppression therapy led to successful pregnancies.

LIMITATIONS, REASONS FOR CAUTION

The number of patients described here is limited and these data need to be confirmed on a larger number of women with non-classic PORD.

WIDER IMPLICATIONS OF THE FINDINGS

The diagnosis of PORD must be considered in infertile women with chronically elevated P and 17OH-P levels and ovarian macrocysts. Differentiation of this entity from non-classic 21OH-D is important, as the multiple enzyme deficiency requires a specific management. Successful fertility induction is possible by IVF, providing that P levels be sufficiently suppressed by glucocorticoid therapy prior to implantation.

STUDY FUNDING/COMPETING INTEREST(S)

No specific funding was used for this study. There are no potential conflicts of interest.

TRIAL REGISTRATION NUMBER

N/A.

Biased cytochrome P450-mediated metabolism via small-molecule ligands binding P450 oxidoreductase

Metabolic control is mediated by the dynamic assemblies and function of multiple redox enzymes. A key element in these assemblies, the P450 oxidoreductase (POR), donates electrons and selectively activates numerous (>50 in humans and >300 in plants) cytochromes P450 (CYPs) controlling metabolism of drugs, steroids and xenobiotics in humans and natural product biosynthesis in plants. The mechanisms underlying POR-mediated CYP metabolism remain poorly understood and to date no ligand binding has been described to regulate the specificity of POR. Here, using a combination of computational modeling and functional assays, we identify ligands that dock on POR and bias its specificity towards CYP redox partners, across mammal and plant kingdom. Single molecule FRET studies reveal ligand binding to alter POR conformational sampling, which results in biased activation of metabolic cascades in whole cell assays. We propose the model of biased metabolism, a mechanism akin to biased signaling of GPCRs, where ligand binding on POR stabilizes different conformational states that are linked to distinct metabolic outcomes. Biased metabolism may allow designing pathway-specific therapeutics or personalized food suppressing undesired, disease-related, metabolic pathways.

Small-angle neutron scattering solution structures of NADPH-dependent sulfite reductase

Sulfite reductase (SiR), a dodecameric complex of flavoprotein reductase subunits (SiRFP) and hemoprotein oxidase subunits (SiRHP), reduces sulfur for biomass incorporation. Electron transfer within SiR requires intra- and inter-subunit interactions that are mediated by the relative position of each protein, governed by flexible domain movements. Using small-angle neutron scattering, we report the first solution structures of SiR heterodimers containing a single copy of each subunit. These structures show how the subunits bind and how both subunit binding and oxidation state impact SiRFP's conformation. Neutron contrast matching experiments on selectively deuterated heterodimers allow us to define the contribution of each subunit to the solution scattering. SiRHP binding induces a change in the position of SiRFP's flavodoxin-like domain relative to its ferredoxin-NADP⁺ reductase domain while compacting SiRHP's N-terminus. Reduction of SiRFP leads to a more open structure relative to its oxidized state, re-positioning SiRFP's N-terminal flavodoxin-like domain towards the SiRHP binding position. These structures show, for the first time, how both SiRHP binding to, and reduction of, SiRFP positions SiRFP for electron transfer between the subunits.

An electron transfer competent structural ensemble of membrane-bound cytochrome P450 1A1 and cytochrome P450 oxidoreductase

Cytochrome P450 (CYP) heme monooxygenases require two electrons for their catalytic cycle. For mammalian microsomal CYPs, key enzymes for xenobiotic metabolism and steroidogenesis and important drug targets and biocatalysts, the electrons are transferred by NADPH-cytochrome P450 oxidoreductase (CPR). No structure of a mammalian CYP-CPR complex has been solved experimentally, hindering understanding of the determinants of electron transfer (ET), which is often rate-limiting for CYP reactions. Here, we investigated the interactions between membrane-bound CYP 1A1, an antitumor drug target, and CPR by a multiresolution computational approach. We find that upon binding to CPR, the CYP 1A1 catalytic domain becomes less embedded in the membrane and reorients, indicating that CPR may affect ligand passage to the CYP active site. Despite the constraints imposed by membrane binding, we identify several arrangements of CPR around CYP 1A1 that are compatible with ET. In the complexes, the interactions of the CPR FMN domain with the proximal side of CYP 1A1 are supplemented by more transient interactions of the CPR NADP domain with the distal side of CYP 1A1. Computed ET rates and pathways agree well with available experimental data and suggest why the CYP-CPR ET rates are low compared to those of soluble bacterial CYPs.

Ferroptosis-Related Flavoproteins: Their Function and Stability

Ferroptosis has been described recently as an iron-dependent cell death driven by peroxidation of membrane lipids. It is involved in the pathogenesis of a number of diverse diseases. From the other side, the induction of ferroptosis can be used to kill tumor cells as a novel therapeutic approach. Because of the broad clinical relevance, a comprehensive understanding of the ferroptosis-controlling protein network is necessary. Noteworthy, several proteins from this network are flavoenzymes. This review is an attempt to present the ferroptosis-related flavoproteins in light of their involvement in anti-ferroptotic and pro-ferroptotic roles. When available, the data on the structural stability of mutants and cofactor-free apoenzymes are discussed. The stability of the flavoproteins could be an important component of the cellular death processes.

Next-Generation Sequencing Revealed Disease-Causing Variants in Two Genes in a Patient With Combined Features of Spherocytosis and Antley-Bixler Syndrome With Genital Anomalies and Disordered Steroidogenesis

Conventionally, patients with combined rare diseases are often difficult to diagnose. This is because some clinicians tend to consider the multiple disease symptoms as the presentation of a complicated “syndrome.” This pattern of thinking also confines their way of filtering pathogenic mutations. Some real pathogenic mutations might be ignored due to not covering all disease presentations. Here we report the case of a girl who was suffering from spherocytosis and Antley-Bixler syndrome with genital anomalies and disordered steroidogenesis. She remained undiagnosed even after targeted gene detection before. However, after performing next-generation sequencing and analyzing the sequencing data, we identified two mutations: c.2978T > A in ANK1 and c.1370G > A in POR. Our findings and experiences in diagnosing these mutations could contribute to the existing knowledge on the clinical and genetic diagnosis of patients with disease presentations in multiple systems.

Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding

The activity of microsomal cytochromes P450 (CYP) is strictly dependent on the supply of electrons provided by NADPH cytochrome P450 oxidoreductase (CPR). The variant nature of the isoform-specific proximal interface of microsomal CYPs implies that the interacting interface between the two proteins is degenerated. Recently, we demonstrated that specific CPR mutations in the FMN-domain (FD) may induce a gain in activity for a specific CYP isoform. In the current report, we confirm the CYP isoform dependence of CPR's degenerated binding by demonstrating that the effect of four of the formerly studied FD mutants are indeed exclusive of a specific CYP isoform, as verified by cytochrome c inhibition studies. Moreover, the nature of CYP's substrate seems to have a modulating role in the CPR:CYP interaction. In silico molecular dynamics simulations of the FD evidence that mutations induces very subtle structural alterations, influencing the characteristics of residues formerly implicated in the CPR:CYP interaction or in positioning of the FMN moiety. CPR seems therefore to be able to form effective interaction complexes with its structural diverse partners via a combination of specific structural features of the FD, which are functional in a CYP isoform dependent manner, and dependent on the substrate bound.

Insight into the structural and functional analysis of the impact of missense mutation on cytochrome P450 oxidoreductase

Cytochrome P450 oxidoreductase (POR) is a steroidogenic and drug-metabolizing enzyme. It helps in the NADPH dependent transfer of electrons to cytochrome P450 (CYP) enzymes for their biological activity. In this study, we employed integrative computational approaches to decipher the impact of proline to leucine missense mutation at position 384 (P384L) in the connecting/hinge domain region which is essential for the catalytic activity of POR. Analysis of protein stability using DUET, MUpro, CUPSAT, I-Mutant2.0, iStable and SAAFEC servers predicted that mutation might alter the structural stability of POR. The significant conformational changes induced by the mutation to the POR structure were analyzed by long-range molecular dynamics simulation. The results revealed that missense mutation decreased the conformational stability of POR as compared to wild type (WT). The PCA based FEL analysis described the mutant-specific conformational alterations and dominant motions essential for the biological activity of POR. The LIGPLOT interaction analysis showed the different binding architecture of FMN, FAD, and NADPH as a result of mutation. The increased number of hydrogen bonds in the FEL conformation of WT proved the strong binding of cofactors in the binding pocket as compared to the mutant. The porcupine plot analysis associated with cross-correlation analysis depicted the high-intensity flexible motion exhibited by functionally important FAD and NADPH binding domain regions. The computational findings unravel the impact of mutation at the structural level, which could be helpful in understanding the molecular mechanism of drug metabolism.

Exploring the catalytic cascade of cembranoid biosynthesis by combination of genetic engineering and molecular simulations

While chemical steps involved in bioactive cembranoid biosynthesis have been examined, the corresponding enzymatic mechanisms leading to their formation remain elusive. In the tobacco plant, Nicotiana tabacum, a putative cembratriene-ol synthase (CBTS) initiates the catalytic cascade that lead to the biosynthesis of cembratriene-4,6-diols, which displays antibacterial- and anti-proliferative activities. We report here on structural homology models, functional studies, and mechanistic explorations of this enzyme using a combination of biosynthetic and computational methods. This approach guided us to develop an efficient de novo production of five bioactive non- and monohydroxylated cembranoids. Our homology models in combination with quantum and classical simulations suggested putative principles of the CBTS catalytic cycle, and provided a possible rationale for the formation of premature olefinic side products. Moreover, the functional reconstruction of a N. tabacum-derived class II P450 with a cognate CPR, obtained by transcriptome mining provided for production of bioactive cembratriene-4,6-diols. Our combined findings provide mechanistic insights into cembranoid biosynthesis, and a basis for the sustainable industrial production of highly valuable bioactive cembranoids.

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.

Biochemical and structural insights into the cytochrome P450 reductase from Candida tropicalis

Cytochrome P450 reductases (CPRs) are diflavin oxidoreductases that supply electrons to type II cytochrome P450 monooxygenases (CYPs). In addition, it can also reduce other proteins and molecules, including cytochrome c, ferricyanide, and different drugs. Although various CPRs have been functionally and structurally characterized, the overall mechanism and its interaction with different redox acceptors remain elusive. One of the main problems regarding electron transfer between CPRs and CYPs is the so-called “uncoupling”, whereby NAD(P)H derived electrons are lost due to the reduced intermediates’ (FAD and FMN of CPR) interaction with molecular oxygen. Additionally, the decay of the iron-oxygen complex of the CYP can also contribute to loss of reducing equivalents during an unproductive reaction cycle. This phenomenon generates reactive oxygen species (ROS), leading to an inefficient reaction. Here, we present the study of the CPR from Candida tropicalis (CtCPR) lacking the hydrophobic N-terminal part (Δ2–22). The enzyme supports the reduction of cytochrome c and ferricyanide, with an estimated 30% uncoupling during the reactions with cytochrome c. The ROS produced was not influenced by different physicochemical conditions (ionic strength, pH, temperature). The X-ray structures of the enzyme were solved with and without its cofactor, NADPH. Both CtCPR structures exhibited the closed conformation. Comparison with the different solved structures revealed an intricate ionic network responsible for the regulation of the open/closed movement of CtCPR.

Biochemical and structural insights into the cytochrome P450 reductase from Candida tropicalis

Cytochrome P450 reductases (CPRs) are diflavin oxidoreductases that supply electrons to type II cytochrome P450 monooxygenases (CYPs). In addition, it can also reduce other proteins and molecules, including cytochrome c, ferricyanide, and different drugs. Although various CPRs have been functionally and structurally characterized, the overall mechanism and its interaction with different redox acceptors remain elusive. One of the main problems regarding electron transfer between CPRs and CYPs is the so-called "uncoupling", whereby NAD(P)H derived electrons are lost due to the reduced intermediates' (FAD and FMN of CPR) interaction with molecular oxygen. Additionally, the decay of the iron-oxygen complex of the CYP can also contribute to loss of reducing equivalents during an unproductive reaction cycle. This phenomenon generates reactive oxygen species (ROS), leading to an inefficient reaction. Here, we present the study of the CPR from Candida tropicalis (CtCPR) lacking the hydrophobic N-terminal part (Δ2-22). The enzyme supports the reduction of cytochrome c and ferricyanide, with an estimated 30% uncoupling during the reactions with cytochrome c. The ROS produced was not influenced by different physicochemical conditions (ionic strength, pH, temperature). The X-ray structures of the enzyme were solved with and without its cofactor, NADPH. Both CtCPR structures exhibited the closed conformation. Comparison with the different solved structures revealed an intricate ionic network responsible for the regulation of the open/closed movement of CtCPR.

Variability in Loss of Multiple Enzyme Activities Due to the Human Genetic Variation P284T Located in the Flexible Hinge Region of NADPH Cytochrome P450 Oxidoreductase

Cytochromes P450 located in the endoplasmic reticulum require NADPH cytochrome P450 oxidoreductase (POR) for their catalytic activities. Mutations in POR cause multiple disorders in humans related to the biosynthesis of steroid hormones and also affect drug-metabolizing cytochrome P450 activities. Electron transfer in POR occurs from NADH to FAD to FMN, and the flexible hinge region in POR is essential for domain movements to bring the FAD and FMN close together for electron transfer. We tested the effect of variations in the hinge region of POR to check if the effects would be similar across all redox partners or there will be differences in activities. Here we are reporting the effects of a POR genetic variant P284T located in the hinge region of POR that is necessary for the domain movements and internal electron transfer between co-factors. Human wild-type and P284T mutant of POR and cytochrome P450 proteins were expressed in bacteria, purified, and reconstituted for enzyme assays. We found that for the P284T variant of POR, the cytochrome c reduction activity was reduced to 47% of the WT and MTT reduction was reduced to only 15% of the WT. No impact on ferricyanide reduction activity was observed, indicating intact direct electron transfer from FAD to ferricyanide, but a severe loss of CYP19A1 (aromatase) activity was observed (9% of WT). In the assays of drug-metabolizing cytochrome P450 enzymes, the P284T variant of POR showed 26% activity for CYP2C9, 44% activity for CYP2C19, 23% activity for CYP3A4, and 44% activity in CYP3A5 assays compared to the WT POR. These results indicate a severe effect on several cytochrome P450 activities due to the P284T variation in POR, which suggests a negative impact on both the steroid as well as drug metabolism in the individuals carrying this variation. The negative impact of P284T mutation in the hinge region of POR seems to be due to disruption of FAD to FMN electron transfer. These results further emphasize the importance of hinge region in POR for protein flexibility and electron transfer within POR as well as the interaction of POR with different redox partners.

Rational evolution of the cofactor-binding site of cytochrome P450 reductase yields variants with increased activity towards specific cytochrome P450 enzymes

NADPH-cytochrome P450 reductase (CPR) is the natural redox partner of microsomal cytochrome P450 enzymes. CPR shows a stringent preference for NADPH over the less expensive cofactor, NADH, economically limiting its use as a biocatalyst. The complexity of cofactor-linked CPR protein dynamics and the incomplete understanding of the interaction of CPR with both cofactors and electron acceptors present challenges for the successful rational engineering of a CPR with enhanced activity with NADH. Here, we report a rational evolution approach to enhance the activity of CPR with NADH, in which mutations were introduced into the NADPH-binding flavin adenine dinucleotide (FAD) domain. Multiple CPR mutants that used NADH more effectively than the wild-type CPR in the reduction of the surrogate electron acceptor, cytochrome c were found. However, most were inactive in supporting P450 activity, arguing against the use of cytochrome c as a surrogate electron acceptor. Unexpectedly, several mutants showed significantly improved activity towards CYP2C9 (mutant 1-014) and/or CYP2A6 (mutants 1-014, 1-015, 1-053 and 1-077) using NADPH, even though the mutations were introduced at locations remote from the putative CPR-P450 interaction face. Therefore, mutations at sites in the FAD domain of CPR may be promising future engineering targets to enhance P450-mediated substrate turnover. ENZYMES: NADPH-cytochrome P450 reductase - EC 1.6.2.4; cytochrome P450 - EC 1.14.14.1.

Activation of the p38/MAPK pathway regulates autophagy in response to the CYPOR-dependent oxidative stress induced by zearalenone in porcine intestinal epithelial cells

Zearalenone (ZEA) can widely contaminate crops and agricultural products. The ingestion of ZEA-contaminated food or feed affects the integrity and functions of the intestines. In this study, we aimed to find the potential protective mechanism against ZEA ingestion. We found that ZEA induced cell death in IPEC-J2 cells. Meanwhile, the cytoprotective autophagy was activated in ZEA-treated cells. Further studies demonstrated that a p38/MAPK inhibitor down-regulated autophagy and increased cell death compared to those of the controls. Furthermore, ZEA could induce the accumulation of ROS, and eliminating ROS with NAC resulted in a decline in cell death, p38/MAPK phosphorylation, and the expression of LC3-II compared to those of ZEA-group. In addition, cytochrome P450 reductase (CYPOR) was significantly increased in ZEA-treated cells compared to that in the controls, and an inhibitor of CYPOR decreased ROS levels and mitigated cell death compared to those of the ZEA-group. More importantly, we found that blocking both p38/MAPK signalling and autophagy could enhance CYPOR expression and elevate ROS levels. Overall, our study indicated that the p38/MAPK pathway could activate protective autophagy in response to the CYPOR-dependent oxidative stress that was induced by ZEA in IPEC-J2 cells.

Structural insight into the high reduction potentials observed for Fusobacterium nucleatum flavodoxin

Flavodoxins are small flavin mononucleotide (FMN)-containing proteins that mediate a variety of electron transfer processes. The primary sequence of flavodoxin from Fusobacterium nucleatum, a pathogenic oral bacterium, is marked with a number of distinct features including a glycine to lysine (K13) substitution in the highly conserved phosphate-binding loop (T/S-X-T-G-X-T), variation in the aromatic residues that sandwich the FMN cofactor, and a more even distribution of acidic and basic residues. The E_ox/sq (oxidized/semiquinone; -43 mV) and E_sq/hq (semiquinone/hydroquinone; -256 mV) are the highest recorded reduction potentials of known long-chain flavodoxins. These more electropositive values are a consequence of the apoprotein binding to the FMN hydroquinone anion with ~70-fold greater affinity compared to the oxidized form of the cofactor. Inspection of the FnFld crystal structure revealed the absence of a hydrogen bond between the protein and the oxidized FMN N5 atom, which likely accounts for the more electropositive E_ox/sq . The more electropositive E_sq/hq is likely attributed to only one negatively charged group positioned within 12 Å of the FMN N1. We show that natural substitutions of highly conserved residues partially account for these more electropositive reduction potentials.

Structural and Functional Studies of the Membrane-Binding Domain of NADPH-Cytochrome P450 Oxidoreductase

NADPH-cytochrome P450 oxidoreductase (CYPOR), the essential flavoprotein of the microsomal cytochrome P450 monooxygenase system, is anchored in the phospholipid bilayer by its amino-terminal membrane-binding domain (MBD), which is necessary for efficient electron transfer to cytochrome P450. Although crystallographic and kinetic studies have established the structure of the soluble catalytic domain and the role of conformational motions in the control of electron transfer, the role of the MBD is largely unknown. We examined the role of the MBD in P450 catalysis through studies of amino-terminal deletion mutants and site-directed spin labeling. We show that the MBD spans the membrane and present a model for the orientation of CYPOR on the membrane capable of forming a complex with cytochrome P450. EPR power saturation measurements of CYPOR mutants in liposomes containing a lipid/Ni(II) chelate identified a region of the soluble domain interacting with the membrane. The deletion of more than 29 residues from the N-terminus of CYPOR decreases cytochrome P450 activity concomitant with alterations in electrophoretic mobility and an increased resistance to protease digestion. The altered kinetic properties of these mutants are consistent with electron transfer through random collisions rather than via formation of a stable CYPOR-P450 complex. Purified MBD binds weakly to cytochrome P450, suggesting that other interactions are also required for CYPOR-P450 complex formation. We propose that the MBD and flexible tether region of CYPOR, residues 51-63, play an important role in facilitating the movement of the soluble domain relative to the membrane and in promoting multiple orientations that permit specific interactions of CYPOR with its varied partners.

Post-Translational Regulation of CYP450s Metabolism As Revealed by All-Atoms Simulations of the Aromatase Enzyme

Phosphorylation by kinases enzymes is a widespread regulatory mechanism able of rapidly altering the function of target proteins. Among these are cytochrome P450s (CYP450), a superfamily of enzymes performing the oxidation of endogenous and exogenous substrates thanks to the electron supply of a redox partner. In spite of its pivotal role, the molecular mechanism by which phosphorylation modulates CYP450s metabolism remains elusive. Here by performing microsecond-long all-atom molecular dynamics simulations, we disclose how phosphorylation regulates estrogen biosynthesis, catalyzed by the Human Aromatase (HA) enzyme. Namely, we unprecedentedly propose that HA phosphorylation at Y361 markedly stabilizes its adduct with the flavin mononucleotide domain of CYP450s reductase (CPR), the redox partner of microsomal CYP450s, and a variety of other proteins. With CPR present at physiological conditions in a limiting ratio with respect to its multiple oxidative partners, the enhanced stability of the CPR/HA adduct may favor HA in the competition with the other proteins requiring CPR's electron supply, ultimately facilitating the electron transfer and estrogen biosynthesis. As a result, our work elucidates at atomic-level the post-translational regulation of CYP450s catalysis. Given the potential for rational clinical management of diseases associated with steroid metabolism disorders, unraveling this mechanism is of utmost importance, and raises the intriguing perspective of exploiting this knowledge to devise novel therapies.

Solution structure of the cytochrome P450 reductase-cytochrome c complex determined by neutron scattering

Electron transfer in all living organisms critically relies on formation of complexes between the proteins involved. The function of these complexes requires specificity of the interaction to allow for selective electron transfer but also a fast turnover of the complex, and they are therefore often transient in nature, making them challenging to study. Here, using small-angle neutron scattering with contrast matching with deuterated protein, we report the solution structure of the electron transfer complex between cytochrome P450 reductase (CPR) and its electron transfer partner cytochrome c This is the first reported solution structure of a complex between CPR and an electron transfer partner. The structure shows that the interprotein interface includes residues from both the FMN- and FAD-binding domains of CPR. In addition, the FMN is close to the heme of cytochrome c but distant from the FAD, indicating that domain movement is required between the electron transfer steps in the catalytic cycle of CPR. In summary, our results reveal key details of the CPR catalytic mechanism, including interactions of two domains of the reductase with cytochrome c and motions of these domains relative to one another. These findings shed light on interprotein electron transfer in this system and illustrate a powerful approach for studying solution structures of protein-protein complexes.

Structural and Kinetic Studies of Asp632 Mutants and Fully Reduced NADPH-Cytochrome P450 Oxidoreductase Define the Role of Asp632 Loop Dynamics in the Control of NADPH Binding and Hydride Transfer

Conformational changes in NADPH-cytochrome P450 oxidoreductase (CYPOR) associated with electron transfer from NADPH to electron acceptors via FAD and FMN have been investigated via structural studies of the four-electron-reduced NADP⁺-bound enzyme and kinetic and structural studies of mutants that affect the conformation of the mobile Gly631-Asn635 loop (Asp632 loop). The structure of four-electron-reduced, NADP⁺-bound wild type CYPOR shows the plane of the nicotinamide ring positioned perpendicular to the FAD isoalloxazine with its carboxamide group forming H-bonds with N1 of the flavin ring and the Thr535 hydroxyl group. In the reduced enzyme, the C8-C8 atoms of the two flavin rings are ∼1 Å closer than in the fully oxidized and one-electron-reduced structures, which suggests that flavin reduction facilitates interflavin electron transfer. Structural and kinetic studies of mutants Asp632Ala, Asp632Phe, Asp632Asn, and Asp632Glu demonstrate that the carboxyl group of Asp632 is important for stabilizing the Asp632 loop in a retracted position that is required for the binding of the NADPH ribityl-nicotinamide in a hydride-transfer-competent conformation. Structures of the mutants and reduced wild type CYPOR permit us to identify a possible pathway for NADP(H) binding to and release from CYPOR. Asp632 mutants unable to form stable H-bonds with the backbone amides of Arg634, Asn635, and Met636 exhibit decreased catalytic activity and severely impaired hydride transfer from NADPH to FAD, but leave interflavin electron transfer intact. Intriguingly, the Arg634Ala mutation slightly increases the cytochrome P450 2B4 activity. We propose that Asp632 loop movement, in addition to facilitating NADP(H) binding and release, participates in domain movements modulating interflavin electron transfer.

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.

Orchestrated Domain Movement in Catalysis by Cytochrome P450 Reductase

NADPH-cytochrome P450 reductase is a multi-domain redox enzyme which is a key component of the P450 mono-oxygenase drug-metabolizing system. We report studies of the conformational equilibrium of this enzyme using small-angle neutron scattering, under conditions where we are able to control the redox state of the enzyme precisely. Different redox states have a profound effect on domain orientation in the enzyme and we analyse the data in terms of a two-state equilibrium between compact and extended conformations. The effects of ionic strength show that the presence of a greater proportion of the extended form leads to an enhanced ability to transfer electrons to cytochrome c. Domain motion is intrinsically linked to the functionality of the enzyme, and we can define the position of the conformational equilibrium for individual steps in the catalytic cycle.

Single-Protein Tracking Reveals That NADPH Mediates the Insertion of Cytochrome P450 Reductase into a Biomimetic of the Endoplasmic Reticulum

Cytochrome P450 reductase (CPR) is the redox partner for most human cytochrome P450 enzymes. It is also believed that CPR is an integral membrane protein exclusively. Herein, we report that, contrary to this belief, CPR can exist as a peripheral membrane protein in the absence of NADPH and will transition to an integral membrane protein in the presence of stoichiometric amounts of NADPH or greater. All experiments were performed in a solid-supported cushioned lipid bilayer that closely matched the chemical composition of the human endoplasmic reticulum and served as an ER biomimetic. The phase characteristics and fluidity of the ER biomimetic was characterized with fluorescence micrographs and temperature-dependent fluorescence recovery after photobleaching. The interactions of CPR with the ER biomimetic were directly observed by tracking single CPR molecules using time-lapse single-molecule fluorescence imaging and subsequent analysis of tracks. These studies revealed dramatic changes in diffusion coefficient and the degree of partitioning of CPR as a function of NADPH concentration.

Impact on CYP19A1 activity by mutations in NADPH cytochrome P450 oxidoreductase

Cytochrome P450 aromatase (CYP19A1), in human placenta metabolizes androgens to estrogens and uses reduced nicotinamide adenine dinucleotide phosphate through cytochrome P450 oxidoreductase (POR) for the energy requirements of its metabolic activities. Mutations in the human POR lead to congenital adrenal hyperplasia due to loss of activities of several steroid metabolizing enzymatic reactions conducted by the cytochrome P450 proteins located in the endoplasmic reticulum. Effect of POR mutations on different P450 activities depend on individual partner proteins. In this report we have studied the impact of mutations found in the POR on the enzymatic activity of CYP19A1. We expressed wild type as well mutant human POR proteins in bacteria and purified the recombinant proteins, which were then used in an in vitro reconstitution system in combination with CYP19A1 and lipids for enzymatic analysis. We found that several mutations as well as polymorphisms in human POR can cause reduction of CYP19A1 activity. This would affect metabolism of estrogens in people with variations of POR allele. The POR mutants Y181D and R616X were found to have no activity in supporting CYP19A1 reactions. The POR mutations Y607C and delF646 showed a loss of 60-90% activity and two polymorphic forms of POR, R316W and G413S showed similar to WT activity. One POR variant, Q153R had almost double the activity of WT. Loss of CYP19A1 activity may contribute to disordered steroidogenesis in female patients with POR mutations as well as in mothers with POR variants carrying a male child.

P450 Oxidoreductase deficiency: Analysis of mutations and polymorphisms

Cytochrome P450 oxidoreductase (POR) is required for metabolic reactions of steroid and drug metabolizing cytochrome P450 proteins located in endoplasmic reticulum. Mutations in POR cause a complex set of disorders resembling combined deficiencies of multiple steroid metabolizing enzymes. The P450 oxidoreductase deficiency (PORD) was first reported in patients with symptoms of defects in steroidogenic cytochrome P450 enzymes and ambiguous genitalia, and bone malformation features resembling Antley-Bixler syndrome. POR is now classified as a separate and rare form of congenital adrenal hyperplasia (CAH), which may cause disorder of sexual development (DSD). Since the initial description of PORD in 2004, a large number of POR mutations and polymorphisms have been described. In this report we have performed computational analysis of mutations and polymorphisms in POR linked to metabolism of steroids and xenobiotics and pathology of PORD from the reported cases. The mutations in POR that were identified in patients with disruption of steroidogenesis also have severe effects on cytochrome P450 proteins involved in metabolism of drugs. Different variations in POR show a range of diverse effects on different partner proteins that are often linked to the location of the particular variants. The variations in POR that cause defective binding of co-factors always have damaging effects on all partner proteins, while the mutations causing subtle structural changes may lead to altered interaction with partner proteins and the overall effect may be different for each individual partner. Computational analysis of available sequencing data and mutation analysis shows that Japanese (R457H), Caucasian (A287P) and Turkish (399-401) populations can be linked to unique founder mutations. Other mutations identified so far were identified as rare alleles or in single isolated reports. The common polymorphism of POR is the variant A503V which can be found in about 27% of alleles in general population but there are remarkable differences among different sub populations.

Instability of the Human Cytochrome P450 Reductase A287P Variant Is the Major Contributor to Its Antley-Bixler Syndrome-like Phenotype

Human NADPH-cytochrome P450 oxidoreductase (POR) gene mutations are associated with severe skeletal deformities and disordered steroidogenesis. The human POR mutation A287P presents with disordered sexual development and skeletal malformations. Difficult recombinant expression and purification of this POR mutant suggested that the protein was less stable than WT. The activities of cytochrome P450 17A1, 19A1, and 21A2, critical in steroidogenesis, were similar using our purified, full-length, unmodified A287P or WT POR, as were those of several xenobiotic-metabolizing cytochromes P450, indicating that the A287P protein is functionally competent in vitro, despite its functionally deficient phenotypic behavior in vivo Differential scanning calorimetry and limited trypsinolysis studies revealed a relatively unstable A287P compared with WT protein, leading to the hypothesis that the syndrome observed in vivo results from altered POR protein stability. The crystal structures of the soluble domains of WT and A287P reveal only subtle differences between them, but these differences are consistent with the differential scanning calorimetry results as well as the differential susceptibility of A287P and WT observed with trypsinolysis. The relative in vivo stabilities of WT and A287P proteins were also examined in an osteoblast cell line by treatment with cycloheximide, a protein synthesis inhibitor, showing that the level of A287P protein post-inhibition is lower than WT and suggesting that A287P may be degraded at a higher rate. Current studies demonstrate that, unlike previously described mutations, A287P causes POR deficiency disorder due to conformational instability leading to proteolytic susceptibility in vivo, rather than through an inherent flavin-binding defect.

Functioning of Microsomal Cytochrome P450s: Murburn Concept Explains the Metabolism of Xenobiotics in Hepatocytes

Using oxygen and NADPH, the redox enzymes cytochrome P450 (CYP) and its reductase (CPR) work in tandem to carry out the phase I metabolism of a vast majority of drugs and xenobiotics. As per the erstwhile understanding of the catalytic cycle, binding of the substrate to CYP's heme distal pocket allows CPR to pump electrons through a CPR-CYP complex. In turn, this trigger (a thermodynamic push of electrons) leads to the activation of oxygen at CYP's heme-center, to give Compound I, a two-electron deficient enzyme reactive intermediate. The formation of diffusible radicals and reactive oxygen species (DROS, hitherto considered an undesired facet of the system) was attributed to the heme-center. Recently, we had challenged these perceptions and proposed the murburn ("mured burning" or "mild unrestricted burning") concept to explain heme enzymes' catalytic mechanism, electron-transfer phenomena and the regulation of redox equivalents' consumption. Murburn concept incorporates a one-electron paradigm, advocating obligatory roles for DROS. The new understanding does not call for high-affinity substrate-binding at the heme distal pocket of the CYP (the first and the most crucial step of the erstwhile paradigm) or CYP-CPR protein-protein complexations (the operational backbone of the erstwhile cycle). Herein, the dynamics of reduced nicotinamide nucleotides' consumption, peroxide formation and depletion, product(s) formation, etc. was investigated with various controls, by altering reaction variables, environments and through the incorporation of diverse molecular probes. In several CYP systems, control reactions lacking the specific substrate showed comparable or higher peroxide in milieu, thereby discrediting the foundations of the erstwhile hypothesis. The profiles obtained by altering CYP:CPR ratios and the profound inhibitions observed upon the incorporation of catalytic amounts of horseradish peroxidase confirm the obligatory roles of DROS in milieu, ratifying murburn as the operative concept. The mechanism of uncoupling (peroxide/water formation) was found to be dependent on multiple one and two electron equilibriums amongst the reaction components. The investigation explains the evolutionary implications of xenobiotic metabolism, confirms the obligatory role of diffusible reactive species in routine redox metabolism within liver microsomes and establishes that a redox enzyme like CYP enhances reaction rates (achieves catalysis) via a novel (hitherto unknown) modality.

Mutants of Cytochrome P450 Reductase Lacking Either Gly-141 or Gly-143 Destabilize Its FMN Semiquinone

NADPH-cytochrome P450 oxidoreductase transfers electrons from NADPH to cytochromes P450 via its FAD and FMN. To understand the biochemical and structural basis of electron transfer from FMN-hydroquinone to its partners, three deletion mutants in a conserved loop near the FMN were characterized. Comparison of oxidized and reduced wild type and mutant structures reveals that the basis for the air stability of the neutral blue semiquinone is protonation of the flavin N5 and strong H-bond formation with the Gly-141 carbonyl. The ΔGly-143 protein had moderately decreased activity with cytochrome P450 and cytochrome c It formed a flexible loop, which transiently interacts with the flavin N5, resulting in the generation of both an unstable neutral blue semiquinone and hydroquinone. The ΔGly-141 and ΔG141/E142N mutants were inactive with cytochrome P450 but fully active in reducing cytochrome c In the ΔGly-141 mutants, the backbone amide of Glu/Asn-142 forms an H-bond to the N5 of the oxidized flavin, which leads to formation of an unstable red anionic semiquinone with a more negative potential than the hydroquinone. The semiquinone of ΔG141/E142N was slightly more stable than that of ΔGly-141, consistent with its crystallographically demonstrated more rigid loop. Nonetheless, both ΔGly-141 red semiquinones were less stable than those of the corresponding loop in cytochrome P450 BM3 and the neuronal NOS mutant (ΔGly-810). Our results indicate that the catalytic activity of cytochrome P450 oxidoreductase is a function of the length, sequence, and flexibility of the 140s loop and illustrate the sophisticated variety of biochemical mechanisms employed in fine-tuning its redox properties and function.

A well-balanced preexisting equilibrium governs electron flux efficiency of a multidomain diflavin reductase

Diflavin reductases are bidomain electron transfer proteins in which structural reorientation is necessary to account for the various intramolecular and intermolecular electron transfer steps. Using small-angle x-ray scattering and nuclear magnetic resonance data, we describe the conformational free-energy landscape of the NADPH-cytochrome P450 reductase (CPR), a typical bidomain redox enzyme composed of two covalently-bound flavin domains, under various experimental conditions. The CPR enzyme exists in a salt- and pH-dependent rapid equilibrium between a previously described rigid, locked state and a newly characterized, highly flexible, unlocked state. We further establish that maximal electron flux through CPR is conditioned by adjustable stability of the locked-state domain interface under resting conditions. This is rationalized by a kinetic scheme coupling rapid conformational sampling and slow chemical reaction rates. Regulated domain interface stability associated with fast stochastic domain contacts during the catalytic cycle thus provides, to our knowledge, a new paradigm for improving our understanding of multidomain enzyme function.

Molecular Cloning, Expression Pattern and Polymorphisms of NADPH-Cytochrome P450 Reductase in the Bird Cherry-Oat Aphid Rhopalosiphum padi (L.)

NADPH–cytochrome P450 reductase (CPR) plays an important role in the cytochrome P450 (CYP)-mediated metabolism of endogenous and exogenous substrates. CPR has been found to be associated with insecticide metabolism and resistance in many insects. However, information regarding CPR in the bird cherry-oat aphid, Rhopalosiphum padi, is unavailable. In the current study, a full-length cDNA (2,476 bp) of CPR (RpCPR) encoding 681 amino acids was cloned from R. padi. Nucleotide sequence and deduced amino acid sequence analysis showed that RpCPR exhibits characteristics of classical CPRs and shares high identities with those of other insects, especially with the pea aphid, Acyrthosiphon pisum. The mRNA of RpCPR was expressed at all developmental stages, with the highest expression level found in the second instar and the lowest in adult. Expression levels of RpCPR in isoprocarb-resistant and imidacloprid-resistant strains were 3.74- and 3.53-fold higher, respectively, than that of a susceptible strain. RpCPR expression could also be induced by low concentrations (LC₃₀) of isoprocarb and imidacloprid. Moreover, we sequenced the open reading frame (ORF) of RpCPR from 167 field samples collected in 11 geographical populations. Three hundred and thirty-four SNPs were detected, of which, 65 were found in more than two individuals. One hundred and ninety-four missense mutations were present in the amino acid sequence, of which, the P484S mutant had an allele frequency of 35.1%. The present results suggest that RpCPR may play an important role in the P450-mediated insecticide resistance of R. padi to isoprocarb and imidacloprid and possibly other insecticides. Meanwhile, RpCPRmaintains high genetic diversity in natural individuals, which provides the possibility of studying potential correlations between variants and certain special physiological characters.

Mapping of the Allosteric Site in Cholesterol Hydroxylase CYP46A1 for Efavirenz, a Drug That Stimulates Enzyme Activity

Cytochrome P450 46A1 (CYP46A1) is a microsomal enzyme and cholesterol 24-hydroxylase that controls cholesterol elimination from the brain. This P450 is also a potential target for Alzheimer disease because it can be activated pharmacologically by some marketed drugs, as exemplified by efavirenz, the anti-HIV medication. Previously, we suggested that pharmaceuticals activate CYP46A1 allosterically through binding to a site on the cytosolic protein surface, which is different from the enzyme active site facing the membrane. Here we identified this allosteric site for efavirenz on CYP46A1 by using a combination of hydrogen-deuterium exchange coupled to MS, computational modeling, site-directed mutagenesis, and analysis of the CYP46A1 crystal structure. We also mapped the binding region for the CYP46A1 redox partner oxidoreductase and found that the allosteric and redox partner binding sites share a common border. On the basis of the data obtained, we propose the mechanism of CYP46A1 allostery and the pathway for the signal transmission from the P450 allosteric site to the active site.