Coupling of Redox and Structural States in Cytochrome P450 Reductase Studied by Molecular Dynamics Simulation

The Dissociation Process of NADP+ from NADPH-Cytochrome P450 Reductase Studied by Molecular Dynamics Simulation

NADPH-cytochrome P450 reductase (CPR) plays a vital role as a redox partner for mammalian cytochrome P450 enzymes (P450s), facilitating the transfer of two electrons from NADPH to the P450 heme center in a sequential manner. Previous experimental studies revealed substantial domain movements of CPR, transitioning between closed and open states during the electron transfer (ET) cycle. These transitions are essential and are influenced by the binding of NADPH or the release of NADP+. However, the intricate molecular mechanisms governing the CPR-mediated ET cycle have largely remained elusive. This study employed molecular dynamics (MD) simulation techniques to explore the dissociation of NADP+ from CPR, a crucial step preceding the initial ET from CPR to a P450. Alongside the binding structure of NADP+ observed in a crystal structure (pose I), our MD simulations identified an alternative binding structure (pose II). Although pose II exhibits slightly lower stability than pose I, it can be formed through an approximate 210° counterclockwise rotation of the adenine group, with a free energy barrier of only 2.76 kcal/mol. The simulation results further suggest that NADP+ dissociation involves a tentative formation of pose II from pose I before complete dissociation, and that the binding of NADP+ to CPR is primarily governed by nonbonded interactions within the adenosine binding pocket.

Structural insights into the semiquinone form of human Cytochrome P450 reductase by DEER distance measurements between a native flavin and a spin labelled non-canonical amino acid

The flavoprotein Cytochrome P450 reductase (CPR) is the unique electron pathway from NADPH to Cytochrome P450 (CYPs). The conformational dynamics of human CPR in solution, which involves transitions from a "locked/closed" to an "unlocked/open" state, is crucial for electron transfer. To date, however, the factors guiding these changes remain unknown. By Site-Directed Spin Labelling coupled to Electron Paramagnetic Resonance spectroscopy, we have incorporated a non-canonical amino acid onto the flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) domains of soluble human CPR, and labelled it with a specific nitroxide spin probe. Taking advantage of the endogenous FMN cofactor, we successfully measured for the first time, the distance distribution by DEER between the semiquinone state FMNH• and the nitroxide. The DEER data revealed a salt concentration-dependent distance distribution, evidence of an "open" CPR conformation at high salt concentrations exceeding previous reports. We also conducted molecular dynamics simulations which unveiled a diverse ensemble of conformations for the "open" semiquinone state of the CPR at high salt concentration. This study unravels the conformational landscape of the one electron reduced state of CPR, which had never been studied before.

Non-covalent interactions that tune the reactivities of the flavins in bifurcating electron transferring flavoprotein

Bifurcating electron transferring flavoproteins (Bf-ETFs) tune chemically identical flavins to two contrasting roles. To understand how, we used hybrid quantum mechanical molecular mechanical calculations to characterize non-covalent interactions applied to each flavin by the protein. Our computations replicated the differences between the reactivities of the flavins: the electron transferring flavin (^ETflavin) was calculated to stabilize anionic semiquinone (ASQ) as needed to execute its single-electron transfers, whereas the Bf flavin (^Bfflavin) was found to disfavor the ASQ state more than does free flavin and to be less susceptible to reduction. The stability of ^ETflavin ASQ was attributed in part to H-bond donation to the flavin O2 from a nearby His side chain, via comparison of models employing different tautomers of His. This H-bond between O2 and the ET site was uniquely strong in the ASQ state, whereas reduction of ^ETflavin to the anionic hydroquinone (AHQ) was associated with side chain reorientation, backbone displacement and reorganization of its H-bond network including a Tyr from the other domain and subunit of the ETF. The Bf site was less responsive overall, but formation of the ^Bfflavin AHQ allowed a nearby Arg side chain to adopt an alternative rotamer that can H-bond to the ^Bfflavin O4. This would stabilize the anionic ^Bfflavin and rationalize effects of mutation at this position. Thus, our computations provide insights on states and conformations that have not been possible to characterize experimentally, offering explanations for observed residue conservation and raising possibilities that can now be tested.

Improved Bioproduction of the Nylon 12 Monomer by Combining the Directed Evolution of P450 and Enhancing Heme Synthesis

The nylon 12 (PA12) monomer ω-aminododecanoic acid (ω-AmDDA) could be synthesized from lauric acid (DDA) through multi-enzyme cascade transformation using engineered E. coli, with the P450 catalyzing terminal hydroxylation of DDA as a rate-limiting enzyme. Its activity is jointly determined by the heme domain and the reductase domain. To obtain a P450 mutant with higher activity, directed evolution was conducted using a colorimetric high-throughput screening (HTS) system with DDA as the real substrate. After two rounds of directed evolution, a positive double-site mutant (R14R/D629G) with 90.3% higher activity was obtained. Molecular docking analysis, kinetic parameter determination and protein electrophoresis suggested the improved soluble expression of P450 resulting from the synonymous mutation near the N-terminus and the shortened distance of the electron transfer between FMN and FAD caused by D629G mutation as the major reasons for activity improvement. The significantly increased k_cat and unchanged K_m provided further evidence for the increase in electron transfer efficiency. Considering the important role of heme in P450, its supply was strengthened by the metabolic engineering of the heme synthesis pathway. By combining P450-directed evolution and enhancing heme synthesis, 2.02 ± 0.03 g/L of ω-AmDDA was produced from 10 mM DDA, with a yield of 93.6%.

Modulating catalytic activity of a modified flavin analogue via judicially positioned metal ion toward aerobic sulphoxidation

This manuscript describes the synthesis, spectroscopic and crystallographic characterization of a cadmium complex of 10-propoylisoalloxazine-7-carboxylic acid (Flc-Cd). Catalytic activity of Flc-Cd towards aerobic sulphoxidation reaction was investigated in the presence of l-ascorbic acid as the reducing agent. Notably the neutral un-metalated flavin analogue did not show any significant catalytic activity. The design strategy for Flc provides a close proximity of the metal centre to the flavin core without compromising the catalytic site thereby assisting the product formation when compared to unmetallated Flc. Minor enantioselectivity is also observed in cases where unsymmetrical sulphides were used; indicative of the possible involvement of chiral l-ascorbic acid in the intermediate formation.

Design and synthesis of a catalytically efficient metal-flavin complex toward aerobic sulphoxidation was achieved via judicially positioning the metal ion near the catalytic site thereby assisting the intermediate formation.

Redox metabolism for improving whole-cell P450-catalysed terpenoid biosynthesis

The growing preference for producing cytochrome P450-mediated natural products in microbial systems stems from the challenging nature of the organic chemistry approaches. The P450 enzymes are redox-dependent proteins, through which they source electrons from reducing cofactors to drive their activities. Widely researched in biochemistry, most of the previous studies have extensively utilised expensive cell-free assays to reveal mechanistic insights into P450 functionalities in presence of commercial redox partners. However, in the context of microbial bioproduction, the synergic activity of P450- reductase proteins in microbial systems have not been largely investigated. This is mainly due to limited knowledge about their mutual interactions in the context of complex systems. Hence, manipulating the redox potential for natural product synthesis in microbial chassis has been limited. As the potential of redox state as crucial regulator of P450 biocatalysis has been greatly underestimated by the scientific community, in this review, we re-emphasize their pivotal role in modulating the in vivo P450 activity through affecting the product profile and yield. Particularly, we discuss the applications of widely used in vivo redox engineering methodologies for natural product synthesis to provide further suggestions for patterning on P450-based terpenoids production in microbial platforms.

Crystal structures of hydroxymethylbilane synthase complexed with a substrate analog: a single substrate-binding site for four consecutive condensation steps

Hydroxymethylbilane synthase (HMBS), which is involved in the heme biosynthesis pathway, has a dipyrromethane cofactor and combines four porphobilinogen (PBG) molecules to form a linear tetrapyrrole, hydroxymethylbilane. Enzyme kinetic study of human HMBS using a PBG-derivative, 2-iodoporphobilinogen (2-I-PBG), exhibited noncompetitive inhibition with the inhibition constant being 5.4 ± 0.3 µM. To elucidate the reaction mechanism of HMBS in detail, crystal structure analysis of 2-I-PBG-bound holo-HMBS and its reaction intermediate possessing two PBG molecules (ES2), and inhibitor-free ES2 was performed at 2.40, 2.31, and 1.79 Å resolution, respectively. Their overall structures are similar to that of inhibitor-free holo-HMBS, and the differences are limited near the active site. In both 2-I-PBG-bound structures, 2-I-PBG is located near the terminus of the cofactor or the tetrapyrrole chain. The propionate group of 2-I-PBG interacts with the side chain of Arg173, and its acetate group is associated with the side chains of Arg26 and Ser28. Furthermore, the aminomethyl group and pyrrole nitrogen of 2-I-PBG form hydrogen bonds with the side chains of Gln34 and Asp99, respectively. These amino acid residues form a single substrate-binding site, where each of the four PBG molecules covalently binds to the cofactor (or oligopyrrole chain) consecutively, ultimately forming a hexapyrrole chain. Molecular dynamics simulation of the ES2 intermediate suggested that the thermal fluctuation of the lid and cofactor-binding loops causes substrate recruitment and oligopyrrole chain shift needed for consecutive condensation. Finally, the hexapyrrole chain is hydrolyzed self-catalytically to produce hydroxymethylbilane.

Conformational Motion of Ferredoxin Enables Efficient Electron Transfer to Heme in the Full-Length P450TT

Cytochrome P450 monooxygenases (P450s) are versatile biocatalysts used in natural products biosynthesis, xenobiotic metabolisms, and biotechnologies. In P450s, the electrons required for O₂ activation are supplied by NAD(P)H through stepwise electron transfers (ETs) mediated by redox partners. While much is known about the machinery of the catalytic cycle of P450s, the mechanisms of long-range ET are largely unknown. Very recently, the first crystal structure of full-length P450_TT was solved. This enables us to decipher the interdomain ET mechanism between the [2Fe-2S]-containing ferredoxin and the heme, by use of molecular dynamics simulations. In contrast to the "distal" conformation characterized in the crystal structure where the [2Fe-2S] cluster is ∼28 Å away from heme-Fe, our simulations demonstrated a "proximal" conformation of [2Fe-2S] that is ∼17 Å [and 13.7 Å edge-to-edge] away from heme-Fe, which may enable the interdomain ET. Key residues involved in ET pathways and interdomain complexation were identified, some of which have already been verified by recent mutation studies. The conformational transit of ferredoxin between "distal" and "proximal" was found to be controlled mostly by the long-range electrostatic interactions between the ferredoxin domain and the other two domains. Furthermore, our simulations show that the full-length P450_TT utilizes a flexible ET pathway that resembles either P450_Scc or P450_cam. Thus, this study provides a uniform picture of the ET process between reductase domains and heme domain.

Conformational Landscape of Cytochrome P450 Reductase Interactions

Oxidative reactions catalyzed by Cytochrome P450 enzymes (CYPs), which constitute the most relevant group of drug-metabolizing enzymes, are enabled by their redox partner Cytochrome P450 reductase (CPR). Both proteins are anchored to the membrane of the endoplasmic reticulum and the CPR undergoes a conformational change in order to interact with the respective CYP and transfer electrons. Here, we conducted over 22 microseconds of molecular dynamics (MD) simulations in combination with protein-protein docking to investigate the conformational changes necessary for the formation of the CPR-CYP complex. While some structural features of the CPR and the CPR-CYP2D6 complex that we highlighted confirmed previous observations, our simulations revealed additional mechanisms for the conformational transition of the CPR. Unbiased simulations exposed a movement of the whole protein relative to the  membrane, potentially to facilitate interactions with its diverse set of redox partners. Further, we present a structural mechanism for the susceptibility of the CPR to different redox states based on the flip of a glycine residue disrupting the local interaction network that maintains inter-domain proximity. Simulations of the CPR-CYP2D6 complex pointed toward an additional interaction surface of the FAD domain and the proximal side of CYP2D6. Altogether, this study provides novel structural insight into the mechanism of CPR-CYP interactions and underlying conformational changes, improving our understanding of this complex machinery Cytochrome P450 reductase; CPR; conformational; dynamicsrelevant for drug metabolism.

Conformational Equilibrium of NADPH-Cytochrome P450 Oxidoreductase Is Essential for Heme Oxygenase Reaction

Heme oxygenase (HO) catalyzes heme degradation using electrons supplied by NADPH-cytochrome P450 oxidoreductase (CPR). Electrons from NADPH flow first to FAD, then to FMN, and finally to the heme in the redox partner. Previous biophysical analyses suggest the presence of a dynamic equilibrium between the open and the closed forms of CPR. We previously demonstrated that the open-form stabilized CPR (ΔTGEE) is tightly bound to heme-HO-1, whereas the reduction in heme-HO-1 coupled with ΔTGEE is considerably slow because the distance between FAD and FMN in ΔTGEE is inappropriate for electron transfer from FAD to FMN. Here, we characterized the enzymatic activity and the reduction kinetics of HO-1 using the closed-form stabilized CPR (147CC514). Additionally, we analyzed the interaction between 147CC514 and heme-HO-1 by analytical ultracentrifugation. The results indicate that the interaction between 147CC514 and heme-HO-1 is considerably weak, and the enzymatic activity of 147CC514 is markedly weaker than that of CPR. Further, using cryo-electron microscopy, we confirmed that the crystal structure of ΔTGEE in complex with heme-HO-1 is similar to the relatively low-resolution structure of CPR complexed with heme-HO-1 in solution. We conclude that the "open-close" transition of CPR is indispensable for electron transfer from CPR to heme-HO-1.

All-Atom Simulations Disclose How Cytochrome Reductase Reshapes the Substrate Access/Egress Routes of Its Partner CYP450s

Cytochromes P450 enzymes (CYP450s) promote the oxidative metabolism of a variety of substrates via the electrons supplied by the cytochrome P450 reductase (CPR) and upon formation of a CPR/CYP450 adduct. In spite of the pivotal regulatory importance of this process, the impact of CPR binding on the functional properties of its partner CYP450 remains elusive. By performing multiple microsecond-long all-atom molecular dynamics simulations of a 520 000-atom model of a CPR/CYP450 adduct embedded in a membrane mimic, we disclose the molecular terms for their interactions, considering the aromatase (HA) enzyme as a proxy of the CYP450 family. Our study strikingly unveils that CPR binding alters HA's functional motions, bolstering a change in the shape and type of the channels traveled by substrates/products during their access/egress to/from the enzyme's active site. Our outcomes unprecedentedly contribute to extricate the many entangled facets of the CYP450 metabolon, redrafting its intricate panorama from an atomic-level perspective.