Structural analysis of the FMN binding domain of NADPH-cytochrome P-450 oxidoreductase by site-directed mutagenesis

A Novel, Highly Potent NADPH-Dependent Cytochrome P450 Reductase from Waste Liza klunzingeri Liver

The use of marine enzymes as catalysts for biotechnological applications is a topical subject. Marine enzymes usually display better operational properties than their animal, plant or bacterial counterparts, enlarging the range of possible biotechnological applications. Due to the fact that cytochrome P450 enzymes can degrade many different toxic environmental compounds, these enzymes have emerged as valuable tools in bioremediation processes. The present work describes the isolation, purification and biochemical characterization of a liver NADPH-dependent cytochrome P450 reductase (CPR) from the marine fish Liza klunzingeri (LkCPR). Experimental results revealed that LkCPR is a monomer of approximately 75 kDa that is active in a wide range of pH values (6-9) and temperatures (40-60 °C), showing the highest catalytic activity at pH 8 and 50 °C. The activation energy of the enzyme reaction was 16.3 kcal mol^-1 K^-1. The K_M values for cytochrome C and NADPH were 8.83 μM and 7.26 μM, and the k_cat values were 206.79 s^-1 and 202.93 s^-1, respectively. LkCPR displayed a specific activity versus cytochrome C of 402.07 µmol min^-1 mg¹, the highest activity value described for a CPR up to date (3.2-4.7 times higher than the most active reported CPRs) and showed the highest thermostability described for a CPR. Taking into account all these remarkable catalytic features, LkCPR offers great potential to be used as a suitable biocatalyst.

CRABPs Alter all-trans-Retinoic Acid Metabolism by CYP26A1 via Protein-Protein Interactions

Cellular retinoic acid binding proteins (CRABP1 and CRABP2) bind all-trans-retinoic acid (atRA), the active metabolite of vitamin A, with high affinity. CRABP1 and CRABP2 have been shown to interact with the atRA-clearing cytochrome P450 enzymes CYP26B1 and CYP26C1 and with nuclear retinoic acid receptors (RARs). We hypothesized that CRABP1 and CRABP2 also alter atRA metabolism and clearance by CYP26A1, the third key atRA-metabolizing enzyme in the CYP26 family. Based on stopped-flow experiments, atRA bound CRABP1 and CRABP2 with K_d values of 4.7 nM and 7.6 nM, respectively. The unbound atRA K_m values for 4-OH-atRA formation by CYP26A1 were 4.7 ± 0.8 nM with atRA, 6.8 ± 1.7 nM with holo-CRABP1 and 6.1 ± 2.7 nM with holo-CRABP2 as a substrate. In comparison, the apparent k_cat value was about 30% lower (0.71 ± 0.07 min⁻¹ for holo-CRABP1 and 0.75 ± 0.09 min⁻¹ for holo-CRABP2) in the presence of CRABPs than with free atRA (1.07 ± 0.08 min⁻¹). In addition, increasing concentrations in apo-CRABPs decreased the 4-OH-atRA formation rates by CYP26A1. Kinetic analyses suggest that apo-CRABP1 and apo-CRABP2 inhibit CYP26A1 (K_i = 0.39 nM and 0.53 nM, respectively) and holo-CRABPs channel atRA for metabolism by CYP26A1. These data suggest that CRABPs play a critical role in modulating atRA metabolism and cellular atRA concentrations.

Prediction of FMN Binding Sites in Electron Transport Chains Based on 2-D CNN and PSSM Profiles

Flavin mono-nucleotides (FMNs) are cofactors that hold responsibility for carrying and transferring electrons in the electron transport chain stage of cellular respiration. Without being facilitated by FMNs, energy production is stagnant due to the interruption in most of the cellular processes. Investigation on FMN's functions, therefore, can gain holistic understanding about human diseases and molecular information on drug targets. We proposed a deep learning model using a two-dimensional convolutional neural network and position specific scoring matrices that could identify FMN interacting residues with the sensitivity of 83.7 percent, specificity of 99.2 percent, accuracy of 98.2 percent, and Matthews correlation coefficients of 0.85 for an independent dataset containing 141 FMN binding sites and 1,920 non-FMN binding sites. The proposed method outperformed other previous studies using similar evaluation metrics. Our positive outcome can also promote the utilization of deep learning in dealing with various problems in bioinformatics and computational biology.

Roles of cytochrome P450 2A6 in the oxidation of flavone, 4'-hydroxyflavone, and 4'-, 3'-, and 2'-methoxyflavones by human liver microsomes

Nine forms of recombinant cytochrome P450 (P450 or CYP) enzymes were used to study roles of individual P450 enzymes in the oxidation of flavone and some other flavonoids, 4'-hydroxyflavone and 4'-, 3'-, and 2'-methoxyflavones, by human liver microsomes using LC-MS/MS analysis.As has been reported previously , 4'-, 3'-, and 2'-methoxyflavones were preferentially O-demethylated by human liver P450 enzymes to form 4'-, 3'-, and 2'-hydroxylated flavones and also 3',4'-dihydroxyflavone from the former two substrates.In comparisons of product formation by oxidation of these methoxylated flavones, CYP2A6 was found to be a major enzyme catalysing flavone 4'- and 3'-hydroxylations by human liver microsomes but did not play significant roles in 2'-hydroxylation of flavone, O-demethylations of three methoxylated flavones, and the oxidation of 4'-hydroxyflavone to 3',4'-dihydroxyflavone.The effects of anti-CYP2A6 IgG and chemical P450 inhibitors suggested that different P450 enzymes, as well as CYP2A6, catalysed oxidation of these flavonoids at different positions by liver microsomes.These studies suggest that CYP2A6 catalyses flavone 4'- and 3'-hydroxylations in human liver microsomes and that other P450 enzymes have different roles in oxidizing these flavonoids.

Structural and kinetic investigations of the carboxy terminus of NADPH-cytochrome P450 oxidoreductase

The influence of the side chains and positioning of the carboxy-terminal residues of NADPH-cytochrome P450 oxidoreductase (CYPOR) on catalytic activity, structure of the carboxy terminus, and interaction with cofactors has been investigated. A tandem deletion of residues Asp675 and Val676, that was expected to shift the position of the functionally important Trp677, resulted in higher cytochrome c reductase activity than that expected from previous studies on the importance of Asp675 and Trp677 in catalysis. Crystallographic determination of the structure of this variant revealed two conformations of the carboxy terminus. In one conformation (Mol A), the last α-helix is partially unwound, resulting in repositioning of all subsequent residues in β-strand 21, from Arg671 to Leu674 (corresponding to Ser673 and Val676 in the wild type structure). This results in the two C-terminal residues, Trp677 and Ser678, being maintained in their wild type positions, with the indole ring of Trp677 stacked against the isoalloxazine ring of FAD as seen in the wild type structure, and Ser673 occupying a similar position to the catalytic residue, Asp675. The other, more disordered conformation is a mixture of the Mol A conformation and one in which the last α-helix is not unwound and the nicotinamide ring is in one of two conformations, out towards the protein surface as observed in the wild type structure (1AMO), or stacked against the flavin ring, similar to that seen in the W677X structure that lacks Trp677 and Ser678 (1JA0). Further kinetic analysis on additional variants showed deletion or substitution of alanine or glycine for Trp677 in conjunction with deletion of Ser678 produced alterations in interactions of CYPOR with NADP+, 2'5'-ADP, and 2'-AMP, as well as the pH dependence of cytochrome c reductase activity. We postulate that deletion of bulky residues at the carboxy terminus permits increased mobility leading to decreased affinity for the 2'5'-ADP and 2'-AMP moieties of NADP+ and subsequent domain movement.

Amphipol-facilitated elucidation of the functional tetrameric complex of full-length cytochrome P450 CYP2B4 and NADPH-cytochrome P450 oxidoreductase

Interactions of membrane-bound mammalian cytochromes P450 (CYPs) with NADPH-cytochrome P450 oxidoreductase (POR), which are required for metabolism of xenobiotics, are facilitated by membrane lipids. A variety of membrane mimetics, such as phospholipid liposomes and nanodiscs, have been used to simulate the membrane to form catalytically active CYP:POR complexes. However, the exact mechanism(s) of these interactions are unclear because of the absence of structural information of full-length mammalian CYP:POR complexes in membranes. Herein, we report the use of amphipols (APols) to form a fully functional, soluble, homogeneous preparation of full-length CYP:POR complexes amenable to biochemical and structural study. Incorporation of CYP2B4 and POR into APols resulted in a CYP2B4:POR complex with a stoichiometry of 1:1, which was fully functional in demethylating benzphetamine at a turnover rate of 37.7 ± 2.2 min^-1, with a coupling efficiency of 40%. Interestingly, the stable complex had a molecular weight (M_w) of 338 ± 22 kDa determined by multiangle light scattering, suggestive of a tetrameric complex of 2CYP2B4:2POR embedded in one APol nanoparticle. Moreover, negative stain electron microscopy (EM) validated the homogeneity of the complex and allowed us to generate a three-dimensional EM map and model consistent with the tetramer observed in solution. This first report of the full-length mammalian CYP:POR complex by transmission EM not only reveals the architecture that facilitates electron transfer but also highlights a potential use of APols in biochemical and structural studies of functional CYP complexes with redox partners.

Heteromeric complex formation between human cytochrome P450 CYP1A1 and heme oxygenase-1

P450 and heme oxygenase-1 (HO-1) receive their necessary electrons by interaction with the NADPH-cytochrome P450 reductase (POR). As the POR concentration is limiting when compared with P450 and HO-1, they must effectively compete for POR to function. In addition to these functionally required protein-protein interactions, HO-1 forms homomeric complexes, and several P450s have been shown to form complexes with themselves and with other P450s, raising the question, 'How are the HO-1 and P450 systems organized in the endoplasmic reticulum?' Recently, CYP1A2 was shown to associate with HO-1 affecting the function of both proteins. The goal of this study was to determine if CYP1A1 formed complexes with HO-1 in a similar manner. Complex formation among POR, HO-1, and CYP1A1 was measured using bioluminescence resonance energy transfer, with results showing HO-1 and CYP1A1 form a stable complex that was further stabilized in the presence of POR. The POR•CYP1A1 complex was readily disrupted by the addition of HO-1. CYP1A1 also was able to affect the POR•HO-1 complex, although the effect was smaller. This interaction between CYP1A1 and HO-1 also affected function, where the presence of CYP1A1 inhibited HO-1-mediated bilirubin formation by increasing the KmPOR•HO-1 without affecting the Vmaxapp. In like manner, HO-1 inhibited CYP1A1-mediated 7-ethoxyresorufin dealkylation by increasing the KmPOR•CYP1A1. Based on the mathematical simulation, the results could not be explained by a model where CYP1A1 and HO-1 simply compete for POR, and are consistent with the formation of a stable CYP1A1•HO-1 complex that affected the functional characteristics of both moieties.

Liquid chromatography-tandem mass spectrometry analysis of oxidation of 2'-, 3'-, 4'- and 6-hydroxyflavanones by human cytochrome P450 enzymes

2'-Hydroxyflavanone (2'OHFva), 3'OHFva, 4'OHFva, and 6OHFva, the major oxidative products of flavanone by human cytochrome P450 (P450, CYP) enzymes, were studied in regard to further oxidation by human CYP1A1, 1A2, 1B1.1, 1B1.3, and 2A6. The products formed were analyzed with LC-MS/MS and characterized by their positive ion fragmentations on mass spectrometry. Several di-hydroxylated flavanone (diOHFva) and di-hydroxylated flavone (diOHFvo) products, detected by analyzing parent ions at m/z 257 and 255, respectively, were found following incubation of these four hydroxylated flavanones with P450s. The m/z 257 products were produced at higher levels than the latter with four substrates examined. The structures of the m/z 257 products were characterized by LC-MS/MS product ion spectra, and the results suggest that 3'OHFva and 4'OHFva are further oxidized mainly at B-ring by P450s while 6OHFva oxidation was at A-ring. Different diOHFvo products (m/z 255) were also characterized by LC-MS/MS, and the results suggested that most of these diOHFvo products were formed through oxidation or desaturation of the diOHFva products (m/z 257) by P450s. Only when 4'OHFva (m/z 241) was used as a substrate, formation of 4'OHFvo (m/z 239) was detected, indicating that diOHFvo might also be formed through oxidation of 4'OHFvo by P450s. Finally, our results indicated that CYP1 family enzymes were more active than CYP2A6 in catalyzing the oxidation of these four hydroxylated flavanones, and these findings were supported by molecular docking studies of these chemicals with active sites of P450 enzymes.

Mechanisms of Herb-Drug Interactions Involving Cinnamon and CYP2A6: Focus on Time-Dependent Inhibition by Cinnamaldehyde and 2-Methoxycinnamaldehyde

Information is scarce regarding pharmacokinetic-based herb-drug interactions (HDI) with trans-cinnamaldehyde (CA) and 2-methoxycinnamaldehyde (MCA), components of cinnamon. Given the presence of cinnamon in food and herbal treatments for various diseases, HDIs involving the CYP2A6 substrates nicotine and letrozole with MCA (K_S = 1.58 µM; Hill slope = 1.16) and CA were investigated. The time-dependent inhibition (TDI) by MCA and CA of CYP2A6-mediated nicotine metabolism is a complex process involving multiple mechanisms. Molecular dynamic simulations showed that CYP2A6's active site accommodates two dynamic ligands. The preferred binding orientations for MCA and CA were consistent with the observed metabolism: epoxidation, O-demethylation, and aromatic hydroxylation of MCA and cinnamic acid formation from CA. The percent remaining activity plots for TDI by MCA and CA were curved, and they were analyzed with a numerical method using models of varying complexity. The best-fit models support multiple inactivator binding, inhibitor depletion, and partial inactivation. Deconvoluted mass spectra indicated that MCA and CA modified CYP2A6 apoprotein with mass additions of 156.79 (142.54-171.04) and 132.67 (123.37-141.98), respectively, and it was unaffected by glutathione. Heme degradation was observed in the presence of MCA (48.5% ± 13.4% loss; detected by liquid chromatography-tandem mass spectrometry). In the absence of clinical data, HDI predictions were made for nicotine and letrozole using inhibition parameters from the best-fit TDI models and parameters scaled from rats. Predicted area under the concentration-time curve fold changes were 4.29 (CA-nicotine), 4.92 (CA-letrozole), 4.35 (MCA-nicotine), and 5.00 (MCA-letrozole). These findings suggest that extensive exposure to cinnamon (corresponding to ≈ 275 mg CA) would lead to noteworthy interactions. SIGNIFICANCE STATEMENT: Human exposure to cinnamon is common because of its presence in food and cinnamon-based herbal treatments. Little is known about the risk for cinnamaldehyde and methoxycinnamaldehyde, two components of cinnamon, to interact with drugs that are eliminated by CYP2A6-mediated metabolism. The interactions with CYP2A6 are complex, involving multiple-ligand binding, time-dependent inhibition of nicotine metabolism, heme degradation, and apoprotein modification. An herb-drug interaction prediction suggests that extensive exposure to cinnamon would lead to noteworthy interactions with nicotine.

Haplotype Diversity of NADPH-Cytochrome P450 Reductase Gene of Ophiocordyceps sinensis and the Effect on Fungal Infection in Host Insects

Ophiocordyceps sinensis Berk. is a fungal parasite that parasitizes the larvae of Hepialidae and is used as a traditional Chinese medicine. However, it is not clear how O. sinensis infects its host. The encoding gene haplotype diversity and predicted function of the nicotinamide adenine dinucleotide phosphate (NADPH)-cytochrome P450 reductase (CPR) related to the fungal pathogenicity was analyzed for 219 individuals from 47 O. sinensis populations. Two NADPH CPR genes of O. sinensis were detected and their dominant haplotypes were widely distributed throughout the entire distribution range in Western China. Only 5.43% of all O. sinensis individuals possessed the specific private haplotypes of NADPH CPR-1 and CPR-2 genes. Bioinformatic analyses predicted that the phosphorylation sites, motifs, and domains of NADPH CPR of O. sinensis were different between those encoding by the dominant and private gene haplotypes. The one-to-one match fungus–host correspondence of the same individual suggested that the widely distributed O. sinensis with the dominant NADPH CPR gene haplotypes may strongly infect almost all host insects through a random infection by oral or respiratory pores. Conversely, O. sinensis with the specific private NADPH CPR gene haplotypes is likely to infect only a few corresponding host insects by breaching the cuticle, due to the changed NADPH CPR structure and function.

Toward a systems approach to cytochrome P450 ensemble: interactions of CYP2E1 with other P450 species and their impact on CYP1A2

In this study, we investigate the ability of ethanol-inducible CYP2E1 to interact with other cytochrome P450 species and affect the metabolism of their substrates. As a model system, we used CYP2E1-enriched human liver microsomes (HLM) obtained by the incorporation of purified CYP2E1. Using a technique based on homo-FRET in oligomers of CYP2E1 labeled with BODIPY 577/618 maleimide we demonstrated that the interactions of CYP2E1 with HLM result in the formation of its mixed oligomers with other P450 species present in the microsomal membrane. Incorporation of CYP2E1 results in a multifold increase in the rate of metabolism of CYP2E1-specific substrates p-Nitrophenol and Chlorzaxozone. The rate of their oxidation remains proportional to the amount of incorporated CYP2E1 up to the content of 0.3-0.4 nmol/mg protein (or ∼50% CYP2E1 in the P450 pool). The incorporated CYP2E1 becomes a fully functional member of the P450 ensemble and do not exhibit any detectable functional differences with the endogenous CYP2E1. Enrichment of HLM with CYP2E1 results in pronounced changes in the metabolism of 7-ethoxy-4-cyanocoumarin (CEC), the substrate of CYP2C19 and CYP1A2 suggesting an increase in the involvement of the latter in its metabolism. This effect goes together with an augmentation of the rate of dealkylation of CYP1A2-specific substrate 7-ethoxyresorufin. Furthermore, probing the interactions of CYP2E1 with model microsomes containing individual P450 enzymes we found that CYP2E1 efficiently interacts with CYP1A2, but lacks any ability to form complexes with CYP2C19. This finding goes inline with CYP2E1-induced redirection of the main route of CEC metabolism from CYP2C19 to CYP1A2.

Preference for O-demethylation reactions in the oxidation of 2'-, 3'-, and 4'-methoxyflavones by human cytochrome P450 enzymes

2'-, 3'-, and 4'-Methoxyflavones (MeFs) were incubated with nine forms of recombinant human cytochrome P450 (P450 or CYP) enzymes in the presence of an NADPH-generating system and the products formed were analyzed with LC-MS/MS methods.CYP1B1.1 and 1B1.3 were highly active in demethylating 4'MeF to form 4'-hydroxyflavone (rate of 5.0 nmol/min/nmol P450) and further to 3',4'-dihydroxyflavone (rates of 2.1 and 0.66 nmol/min/nmol P450, respectively). 3'MeF was found to be oxidized by P450s to m/z 239 (M-14) products (presumably 3'-hydroxyflavone) and then to 3',4'-dihydroxyflavone. P450s also catalyzed oxidation of 2'MeF to m/z 239 (M-14) and m/z 255 (M-14, M-14 + 16) products, presumably mono- and di-hydroxylated products, respectively.At least two types of ring oxidation products having m/z 269 fragments were formed, although at slower rates than the formation of mono- and di-hydroxylated products, on incubation of these MeFs with P450s; one type was products oxidized at the C-ring, having m/z 121 fragments, and the other one was the products oxidized at the A-ring (having m/z 137 fragments).Molecular docking analysis indicated the preference of interaction of O-methoxy moiety of methoxyflavones in the active site of CYP1A2.These results suggest that 2'-, 3'-, and 4'-methoxyflavones are principally demethylated by human P450s to form mono- and di-hydroxyflavones and that direct oxidation occurs in these MeFs to form mono-hydroxylated products, oxidized at the A- or B-ring of MeF.

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₅ 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.

Dynamics and Location of the Allosteric Midazolam Site in Cytochrome P4503A4 in Lipid Nanodiscs

Promiscuous and allosteric drug interactions with cytochrome P450 3A4 (CYP3A4) are ubiquitous but incompletely understood at the molecular level. A classic allosteric CYP3A4 drug interaction includes the benzodiazepine midazolam (MDZ). MDZ exhibits homotropic and heterotropic allostery when metabolized to 1'-hydroxy and 4-hydroxy metabolites in varying ratios. The combination of hydrogen-deuterium exchange mass spectrometry (HDX-MS) and Gaussian accelerated molecular dynamics (GaMD) simulations of CYP3A4 in lipid nanodiscs and in a lipid bilayer, respectively, reveals MDZ-dependent changes in dynamics in a membrane environment. The F-, G-, and intervening helices, as well as the loop preceding the β1-sheets, display the largest observed changes in HDX. The GaMD suggests a potential allosteric binding site for MDZ in the F'- and G'-regions, which undergo significant increases in HDX at near-saturating MDZ concentrations. The HDX-MS and GaMD results confirm that changes in dynamics are most significant near the developing consensus allosteric site, and these changes are distinct from those observed previously with the nonallosteric inhibitor ketoconazole. The results suggest that the allosteric MDZ remains mobile in its binding site at the Phe-cluster. The results further suggest that this binding site remains dynamic or changes the depth of insertion in the membrane.

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.

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.

Lipid composition and macromolecular crowding effects on CYP2J2-mediated drug metabolism in nanodiscs

Lipid composition and macromolecular crowding are key external effectors of protein activity and stability whose role varies between different proteins. Therefore, it is imperative to study their effects on individual protein function. CYP2J2 is a membrane-bound cytochrome P450 in the heart involved in the metabolism of fatty acids and xenobiotics. In order to facilitate this metabolism, cytochrome P450 reductase (CPR), transfers electrons to CYP2J2 from NADPH. Herein, we use nanodiscs to show that lipid composition of the membrane bilayer affects substrate metabolism of the CYP2J2-CPR nanodisc (ND) system. Differential effects on both NADPH oxidation and substrate metabolism by CYP2J2-CPR are dependent on the lipid composition. For instance, sphingomyelin containing nanodiscs produced more secondary substrate metabolites than discs of other lipid compositions, implying a possible conformational change leading to processive metabolism. Furthermore, we demonstrate that macromolecular crowding plays a role in the lipid-solubilized CYP2J2-CPR system by increasing the K_m and decreasing the V_max , and effect that is size-dependent. Crowding also affects the CYP2J2-CPR-ND system by decreasing both the K_m and V_max for Dextran-based macromolecular crowding agents, implying an increase in substrate affinity but a lack of metabolism. Finally, protein denaturation studies show that crowding agents destabilize CYP2J2, while the multidomain protein CPR is stabilized. Overall, these studies are the first report on the role of the surrounding lipid environment and macromolecular crowding in modulating enzymatic function of CYP2J2-CPR membrane protein system.

Molecular mechanism of metabolic NAD(P)H-dependent electron-transfer systems: The role of redox cofactors

NAD(P)H-dependent electron-transfer (ET) systems require three functional components: a flavin-containing NAD(P)H-dehydrogenase, one-electron carrier and metal-containing redox center. In principle, these ET systems consist of one-, two- and three-components, and the electron flux from pyridine nucleotide cofactors, NADPH or NADH to final electron acceptor follows a linear pathway: NAD(P)H → flavin → one-electron carrier → metal containing redox center. In each step ET is primarily controlled by one- and two-electron midpoint reduction potentials of protein-bound redox cofactors in which the redox-linked conformational changes during the catalytic cycle are required for the domain-domain interactions. These interactions play an effective ET reactions in the multi-component ET systems. The microsomal and mitochondrial cytochrome P450 (cyt P450) ET systems, nitric oxide synthase (NOS) isozymes, cytochrome b₅ (cyt b₅) ET systems and methionine synthase (MS) ET system include a combination of multi-domain, and their organizations display similarities as well as differences in their components. However, these ET systems are sharing of a similar mechanism. More recent structural information obtained by X-ray and cryo-electron microscopy (cryo-EM) analysis provides more detail for the mechanisms associated with multi-domain ET systems. Therefore, this review summarizes the roles of redox cofactors in the metabolic ET systems on the basis of one-electron redox potentials. In final Section, evolutionary aspects of NAD(P)H-dependent multi-domain ET systems will be discussed.

Kinetic mechanism of time-dependent inhibition of CYP2D6 by 3,4-methylenedioxymethamphetamine (MDMA): Functional heterogeneity of the enzyme and the reversibility of its inactivation

We investigate the mechanism of time-dependent inhibition (TDI) of human cytochrome P450 2D6 (CYP2D6) by 3,4-methylenedioxymethamphetamine (MDMA, ecstasy), one of the most widespread recreational drugs of abuse. In an effort to unravel the kinetic mechanism of the formation of metabolic inhibitory complex (MIC) of CYP2D6 with MDMA-derived carbene we carried out a series of spectrophotometric studies paralleled with registration of the kinetics of time-dependent inhibition (TDI) in CYP2D6-incorporated proteoliposomes. The high amplitude of spectral signal in this system allowed us to characterize the spectral properties of the formed MIC in details and obtain an accurate spectral signature of MIC formation. This information was then used in the studies with CYP2D6-containing microsomes of insect cells (CYP2D6 Supersomes™). Our results demonstrate that in both systems the formation of the ferrous carbene-derived MIC is relatively slow, reversible and is not associated with the accumulation of the ferric carbene intermediate, as takes place in the case of CYP3A4 and podophylotoxin. Furthermore, the limited amplitude of MIC formation suggests that only a fraction (∼50%) of spectrally detectable CYP2D6 in both proteoliposomes and Supersomes participates in the formation of MIC and is therefore involved in the MDMA metabolism. This observation reveals yet another example of a cytochrome P450 that exhibits persistent functional heterogeneity of its population in microsomal membranes. Our study provides a solid methodological background for further mechanistic studies of MIC formation in human liver microsomes and demonstrates that the potency and physiological relevance of MDMA-dependent TDI of CYP2D6 may be overestimated.

Comparison of Antifungal Azole Interactions with Adult Cytochrome P450 3A4 versus Neonatal Cytochrome P450 3A7

Adult drug metabolism is dominated by cytochrome P450 3A4 (CYP3A4), which is often inhibited by antifungal azole drugs, resulting in potential alterations in drug metabolism and adverse drug/drug interactions. In the fetal and neonatal stages of life, the 87% identical cytochrome P450 3A7 (CYP3A7) is expressed but not CYP3A4. Azole antifungals developed for adults are also used in neonates, assuming they interact similarly with both enzymes, but systematic information is lacking. Herein a method was developed for generating recombinant purified CYP3A7. Thirteen different azoles were then evaluated for binding and inhibition of purified human CYP3A4 versus CYP3A7. All imidazole-containing azoles bound both enzymes via coordination to the heme iron and inhibited both with IC₅₀ values ranging from 180 nM for clotrimazole to the millimolar range for imidazole itself. Across this wide range of potencies, CYP3A4 was consistently inhibited more strongly than CYP3A7, with clotrimazole being the least selective (1.5-fold) inhibitor and econazole the most selective (12-fold). Observations for 1,2,4-triazole-containing azoles were more varied. Most bound to CYP3A4 via coordination to the heme iron, but several also demonstrated evidence of a distinct binding mode at low concentrations. However, only posaconazole inhibited CYP3A4. Of the triazoles, only posaconazole inhibited CYP3A7, again less potently than CYP3A4. Spectral evidence for binding was weak or nonexistent for all triazoles. Overall, although the details of binding interactions do vary, the same azole compounds inhibit both enzymes, albeit with weaker interactions with CYP3A7 compared with CYP3A4.

CYP2D6 Allelic Variants *34, *17-2, *17-3, and *53 and a Thr309Ala Mutant Display Altered Kinetics and NADPH Coupling in Metabolism of Bufuralol and Dextromethorphan and Altered Susceptibility to Inactivation by SCH 66712

Metabolic phenotype can be affected by multiple factors, including allelic variation and interactions with inhibitors. Human CYP2D6 is responsible for approximately 20% of cytochrome P450-mediated drug metabolism but consists of more than 100 known variants; several variants are commonly found in the population, whereas others are quite rare. Four CYP2D6 allelic variants-three with a series of mutations distal to the active site (*34, *17-2, *17-3) and one ultra-metabolizer with mutations near the active site (*53), along with reference *1 and an active site mutant of *1 (Thr309Ala)-were expressed, purified, and studied for interactions with the typical substrates dextromethorphan and bufuralol and the inactivator SCH 66712. We found that *34, *17-2, and *17-3 displayed reduced enzyme activity and NADPH coupling while producing the same metabolites as *1, suggesting a possible role for Arg296 in NADPH coupling. A higher-activity variant, *53, displayed similar NADPH coupling to *1 but was less susceptible to inactivation by SCH 66712. The Thr309Ala mutant showed similar activity to that of *1 but with greatly reduced NADPH coupling. Overall, these results suggest that kinetic and metabolic analysis of individual CYP2D6 variants is required to understand their possible contributions to variable drug response and the complexity of personalized medicine.

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.

Use of Phenoxyaniline Analogues To Generate Biochemical Insights into the Interactio n of Polybrominated Diphenyl Ether with CYP2B Enzymes

Human hepatic cytochromes P450 (CYP) are integral to xenobiotic metabolism. CYP2B6 is a major catalyst of biotransformation of environmental toxicants, including polybrominated diphenyl ethers (PBDEs). CYP2B substrates tend to contain halogen atoms, but the biochemical basis for this selectivity and for species specific determinants of metabolism has not been identified. Spectral binding titrations and inhibition studies were performed to investigate interactions of rat CYP2B1, rabbit CYP2B4, and CYP2B6 with a series of phenoxyaniline (POA) congeners that are analogues of PBDEs. For most congeners, there was a <3-fold difference between the spectral binding constants (K_S) and IC₅₀ values. In contrast, large discrepancies between these values were observed for POA and 3-chloro-4-phenoxyaniline. CYP2B1 was the enzyme most sensitive to POA congeners, so the Val-363 residue from that enzyme was introduced into CYP2B4 or CYP2B6. This substitution partially altered the protein-ligand interaction profiles to make them more similar to that of CYP2B1. Addition of cytochrome P450 oxidoreductase (POR) to titrations of CYP2B6 with POA or 2'4'5'TCPOA decreased the affinity of both ligands for the enzyme. Addition of cytochrome b₅ to a recombinant enzyme system containing POR and CYP2B6 increased the POA IC₅₀ value and decreased the 2'4'5'TCPOA IC₅₀ value. Overall, the inconsistency between K_S and IC₅₀ values for POA versus 2'4'5'TCPOA is largely due to the effects of redox partner binding. These results provide insight into the biochemical basis of binding of diphenyl ethers to human CYP2B6 and changes in CYP2B6-mediated metabolism that are dependent on POA congener and redox partner identity.

Human P450 CYP17A1: Control of Substrate Preference by Asparagine 202

CYP17A1 is a key steroidogenic enzyme known to conduct several distinct chemical transformations on multiple substrates. In its hydroxylase activity, this enzyme adds a hydroxyl group at the 17α position of both pregnenolone and progesterone at approximately equal rates. However, the subsequent 17,20 carbon-carbon scission reaction displays variable substrate specificity in the numerous CYP17A1 isozymes operating in vertebrates, manifesting as different K_d and k_cat values when presented with 17α-hydroxypregnenlone (OHPREG) versus 17α-hydroxyprogesterone (OHPROG). Here we show that the identity of the residue at position 202 in human CYP17A1, thought to form a hydrogen bond with the A-ring alcohol substituent on the pregnene- nucleus, is a key driver of this enzyme's native preference for OHPREG. Replacement of asparagine 202 with serine completely reverses the preference of CYP17A1, more than doubling the rate of turnover of the OHPROG to androstenedione reaction and substantially decreasing the rate of formation of dehydroepiandrosterone from OHPREG. In a series of resonance Raman experiments, it was observed that, in contrast with the case for the wild-type protein, in the mutant the 17α alcohol of OHPROG tends to form a H-bond with the proximal rather than terminal oxygen of the oxy-ferrous complex. When OHPREG was a substrate, the mutant enzyme was found to have a H-bonding interaction with the proximal oxygen that is substantially weaker than that of the wild type. These results demonstrate that a single-point mutation in the active site pocket of CYP17A1, even when far from the heme, has profound effects on steroidogenic selectivity in androgen biosynthesis.

Toward a systems approach to the human cytochrome P450 ensemble: interactions between CYP2D6 and CYP2E1 and their functional consequences

Functional cross-talk among human drug-metabolizing cytochrome P450 through their association is a topic of emerging importance. Here, we studied the interactions of human CYP2D6, a major metabolizer of psychoactive drugs, with one of the most prevalent human P450 enzymes, ethanol-inducible CYP2E1. Detection of P450-P450 interactions was accomplished through luminescence resonance energy transfer between labeled proteins incorporated into human liver microsomes and the microsomes of insect cells containing NADPH-cytochrome P450 reductase. The potential of CYP2D6 to form oligomers in the microsomal membrane is among the highest observed with human cytochrome P450 studied up to date. We also observed the formation of heteromeric complexes of CYP2D6 with CYP2E1 and CYP3A4, and found a significant modulation of these interactions by 3,4-methylenedioxymethylamphetamine, a widespread drug of abuse metabolized by CYP2D6. Our results demonstrate an ample alteration of the catalytic properties of CYP2D6 and CYP2E1 caused by their association. In particular, we demonstrated that preincubation of microsomes containing co-incorporated CYP2D6 and CYP2E1 with CYP2D6-specific substrates resulted in considerable time-dependent activation of CYP2D6, which presumably occurs via a slow substrate-induced reorganization of CYP2E1-CYP2D6 hetero-oligomers. Furthermore, we demonstrated that the formation of heteromeric complexes between CYP2E1 and CYP2D6 affects the stoichiometry of futile cycling and substrate oxidation by CYP2D6 by means of decreasing the electron leakage through the peroxide-generating pathways. Our results further emphasize the role of P450-P450 interactions in regulatory cross-talk in human drug-metabolizing ensemble and suggest a role of interactions of CYP2E1 with CYP2D6 in pharmacologically important instances of alcohol-drug interactions.

Altered CYP19A1 and CYP3A4 Activities Due to Mutations A115V, T142A, Q153R and P284L in the Human P450 Oxidoreductase

All cytochromes P450s in the endoplasmic reticulum rely on P450 oxidoreductase (POR) for their catalytic activities. Mutations in POR cause metabolic disorders of steroid hormone biosynthesis and affect certain drug metabolizing P450 activities. We studied mutations A115V, T142A, Q153R identified in the flavin mononucleotide (FMN) binding domain of POR that interacts with partner proteins and P284L located in the hinge region that is required for flexibility and domain movements in POR. Human wild-type (WT) and mutant POR as well as CYP3A4 and CYP19A1 proteins in recombinant form were expressed in bacteria, and purified proteins were reconstituted in liposomes for enzyme kinetic assays. Quality of POR protein was checked by cytochrome c reduction assay as well as flavin content measurements. We found that proteins carrying mutations A115V, T142A located close to the FMN binding site had reduced flavin content compared to WT POR and lost almost all activity to metabolize androstenedione via CYP19A1 and showed reduced CYP3A4 activity. The variant P284L identified from apparently normal subjects also had severe loss of both CYP19A1 and CYP3A4 activities, indicating this to be a potentially disease causing mutation. The mutation Q153R initially identified in a patient with disordered steroidogenesis showed remarkably increased activities of both CYP19A1 and CYP3A4 without any significant change in flavin content, indicating improved protein–protein interactions between POR Q153R and some P450 proteins. These results indicate that effects of mutations on activities of individual cytochromes P450 can be variable and a detailed analysis of each variant with different partner proteins is necessary to accurately determine the genotype-phenotype correlations of POR variants.

Molecular cloning and functional characterization of NADPH-dependent cytochrome P450 reductase from the green microalga Botryococcus braunii, B race

The green microalga Botryococcus braunii of the B race accumulates various lipophilic compounds containing a 10,11-oxidosqualene epoxide moiety in addition to large amounts of triterpene hydrocarbons. While 2,3-squalene epoxidases have already been isolated and characterized from the alga, the enzyme that catalyzes the 10,11-epoxidation of squalene has remained elusive. In order to obtain a molecular tool to explore a 10,11-squalene epoxidase, cDNA cloning of an NADPH-dependent cytochrome P450 reductase (CPR) that is required by both squalene epoxidases and cytochrome P450 enzymes was carried out. The isolated cDNA contained an open reading frame (1998 bp) that encoded for a protein with 665 amino acid residues with a predicted molecular weight of 71.46 kDa and a theoretical pI of 5.49. Analysis of the deduced amino acid sequence revealed the presence of conserved motifs, including FMN, FAD, and NADPH binding domains, which are typical of other CPRs and necessary for enzyme activity. By truncation of the N-terminal transmembrane anchor and addition of a 6× His-tag, BbCPR was heterologously produced in Escherichia coli and purified by Ni-NTA affinity chromatography. The purified recombinant enzyme showed optimal reducing activity of cytochrome c at around a neutral pH at a temperature range of 30-37°C. For steady state kinetic parameters, the recombinant enzyme had a k_m for cytochrome c and NADPH of 11.7±1.6 and 9.4±1.4 μM, and a k_cat for cytochrome c and NADPH of 2.78±0.09 and 3.66±0.11 μmol/min/mg protein, respectively. This is the first study to perform the functional characterization of a CPR from eukaryotic microalgae.

Modulation of the interaction between human P450 3A4 and B. megaterium reductase via engineered loops

Chimerogenesis involving cytochromes P450 is a successful approach to generate catalytically self-sufficient enzymes. However, the connection between the different functional modules should allow a certain degree of flexibility in order to obtain functional and catalytically efficient proteins. We previously applied the molecular Lego approach to develop a chimeric P450 3A4 enzyme linked to the reductase domain of P450 BM3 (BMR). Three constructs were designed with the connecting loop containing no glycine, 3 glycine or 5 glycine residues and showed a different catalytic activity and coupling efficiency. Here we investigate how the linker affects the ability of P450 3A4 to bind substrates and inhibitors. We measure the electron transfer rates and the catalytic properties of the enzyme also in the presence of ketoconazole as inhibitor. The data show that the construct 3A4-5GLY-BMR with the longest loop better retains the binding ability and cooperativity for testosterone, compared to P450 3A4. In both 3A4-3GLY-BMR and 3A4-5GLY-BMR, the substrate induces an increase in the first electron transfer rate and a shorter lag phase related to a domain rearrangements, when compared to the construct without Gly. These data are consistent with docking results and secondary structure predictions showing a propensity to form helical structures in the loop of the 3A4-BMR and 3A4-3GLY-BMR. All three chimeras retain the ability to bind the inhibitor ketoconazole and show an IC₅₀ comparable with those reported for the wild type protein. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Functional characterization of a soluble NADPH-cytochrome P450 reductase from Fusarium graminearum

Fusarium head blight is a devastating disease in wheat caused by some fungal pathogens of the Fusarium genus mainly F. graminearum, due to accumulation of toxic trichothecenes. Most of the trichothecene biosynthetic pathway has been mapped, although some proteins of the pathway remain uncharacterized, including an NADPH-cytochrome P450 reductase. We subcloned a F. graminearum cytochrome P450 reductase that might be involved in the trichothecene biosynthesis. It was expressed heterologously in E. coli as N-terminal truncated form with an octahistidine tag for purification. The construct yielded a soluble apoprotein and its incubation with flavins yielded the corresponding monomeric holoprotein. It was characterized for activity in the pH range 5.5-9.5, using thiazolyl blue tetrazolium bromide (MTT) or cytochrome c as substrates. Binding of the small molecule MTT was weaker than for cytochrome c, however, the rate of MTT reduction was faster. Contrary to other studies of cytochrome reductase proteins, MTT reduction proceeded in a cooperative manner in our studies. Optimum kinetic activity was found at pH 7.5-8.5 for bothMTT and cytochrome c. This is the first paper presenting characterization of a cytochrome P450 reductase from F. graminearum which most likely is involved in mycotoxin biosynthesis or some primary metabolic pathway such as sterol biosynthesis in F. graminearum.

Liver microsomal lipid enhances the activity and redox coupling of colocalized cytochrome P450 reductase-cytochrome P450 3A4 in nanodiscs

The haem-containing mono-oxygenase cytochrome P450 3A4 (CYP3A4) and its redox partner NADPH-dependent cytochrome P450 oxidoreductase (CPR) are among the most important enzymes in human liver for metabolizing drugs and xenobiotic compounds. They are membrane-bound in the endoplasmic reticulum (ER). How ER colocalization and the complex ER phospholipid composition influence enzyme activity are not well understood. CPR and CYP3A4 were incorporated into phospholipid bilayer nanodiscs, both singly, and together in a 1 : 1 ratio, to investigate the significance of membrane insertion and the influence of varying membrane composition on steady-state reaction kinetics. Reaction kinetics were analysed using a fluorimetric assay with 7-benzyloxyquinoline as substrate for CYP3A4. Full activity of the mono-oxygenase system, with electron transfer from NADPH via CPR, could only be reconstituted when CPR and CYP3A4 were colocalized within the same nanodiscs. No activity was observed when CPR and CYP3A4 were each incorporated separately into nanodiscs then mixed together, or when soluble forms of CPR were mixed with preassembled CYP3A4-nanodiscs. Membrane integration and colocalization are therefore essential for electron transfer. Liver microsomal lipid had an enhancing effect compared with phosphatidylcholine on the activity of CPR alone in nanodiscs, and a greater enhancing effect on the activity of CPR-CYP3A4 nanodisc complexes, which was not matched by a phospholipid mixture designed to mimic the ER composition. Furthermore, liver lipid enhanced redox coupling within the system. Thus, natural ER lipids possess properties or include components important for enhanced catalysis by CPR-CYP3A4 nanodisc complexes. Our findings demonstrate the importance of using natural lipid preparations for the detailed analysis of membrane protein activity.

Role of active site loop in coenzyme binding and flavin reduction in cytochrome P450 reductase

Cytochrome P450 reductase (CPR) contains a loop within the active site (comprising Asp(634), Ala(635), Arg(636) and Asn(637); human CPR numbering) that relocates upon NADPH binding. Repositioning of the loop triggers the reorientation of an FAD-shielding tryptophan (Trp(679)) to a partially stacked conformer, reducing the energy barrier for displacement of the residue by the NADPH nicotinamide ring: an essential step for hydride transfer. We used site-directed mutagenesis and kinetic analysis to investigate if the amino acid composition of the loop influences the catalytic properties of CPR. The D634A and D634N variants elicited a modest increase in coenzyme binding affinity coupled with a 36- and 10-fold reduction in cytochrome c(3+) turnover and a 17- and 3-fold decrease in the pre-steady state rate of flavin reduction. These results, in combination with a reduction in the kinetic isotope effect for hydride transfer, suggest that diminished activity is due to destabilization of the partially stacked conformer of Trp(677) and slower release of NADP(+). In contrast, R636A, R636S and an A635G/R636S double mutant led to a modest increase in cytochrome c(3+) reduction, which is linked to weaker coenzyme binding and faster interflavin electron transfer. A potential mechanism by which Arg(636) influences catalysis is discussed.

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.

Selective Targeting of Heme Protein in Cytochrome P450 and Nitric Oxide Synthase by Diphenyleneiodonium

Cytochrome P450 (CYP) enzymes mediate mixed-function oxidation reactions important in drug metabolism. The aromatic heterocyclic cation, diphenyleneiodonium (DPI), binds flavin in cytochrome P450 reductase and inhibits CYP-mediated activity. DPI also inhibits CYP by directly interacting with heme. Herein, we report that DPI effectively inhibits a number of CYP-related monooxygenase reactions including NADPH oxidase, a microsomal enzyme activity that generates hydrogen peroxide in the absence of metabolizing substrates. Inhibition of monooxygenase by DPI was time and concentration dependent with IC50's ranging from 0.06 to 1.9 μM. Higher (4.6-23.9 μM), but not lower (0.06-1.9 μM), concentrations of DPI inhibited electron flow via cytochrome P450 reductase, as measured by its ability to reduce cytochrome c and mediate quinone redox cycling. Similar results were observed with inducible nitric oxide synthase (iNOS), an enzyme containing a C-terminal reductase domain homologous to cytochrome P450 reductase that mediates reduction of cytochrome c, and an N-terminal heme-thiolate oxygenase domain mediating nitric oxide production. Significantly greater concentrations of DPI were required to inhibit cytochrome c reduction by iNOS (IC50 = 3.5 µM) than nitric oxide production (IC50 = 0.16 µM). Difference spectra of liver microsomes, recombinant CYPs, and iNOS demonstrated that DPI altered heme-carbon monoxide interactions. In the presence of NADPH, DPI treatment of microsomes and iNOS yielded a type II spectral shift. These data indicate that DPI interacts with both flavin and heme in CYPs and iNOS. Increased sensitivity for inhibition of CYP-mediated metabolism and nitric oxide production by iNOS indicates that DPI targets heme moieties within the enzymes.

Inactivation of CYP2A6 by the Dietary Phenylpropanoid trans-Cinnamic Aldehyde (Cinnamaldehyde) and Estimation of Interactions with Nicotine and Letrozole

Human exposure to trans-cinnamic aldehyde [t-CA; cinnamaldehyde; cinnamal; (E)-3-phenylprop-2-enal] is common through diet and through the use of cinnamon powder for diabetes and to provide flavor and scent in commercial products. We evaluated the likelihood of t-CA to influence metabolism by inhibition of P450 enzymes. IC50 values from recombinant enzymes indicated that an interaction is most probable for CYP2A6 (IC50 = 6.1 µM). t-CA was 10.5-fold more selective for human CYP2A6 than for CYP2E1; IC50 values for P450s 1A2, 2B6, 2C9, 2C19, 2D6, and 3A4 were 15.8-fold higher or more. t-CA is a type I ligand for CYP2A6 (KS = 14.9 µM). Inhibition of CYP2A6 by t-CA was metabolism-dependent; inhibition required NADPH and increased with time. Glutathione lessened the extent of inhibition modestly and statistically significantly. The carbon monoxide binding spectrum was dramatically diminished after exposure to NADPH and t-CA, suggesting degradation of the heme or CYP2A6 apoprotein. Using a static model and mechanism-based inhibition parameters (K(I) = 18.0 µM; k(inact) = 0.056 minute(-1)), changes in the area under the concentration-time curve (AUC) for nicotine and letrozole were predicted in the presence of t-CA (0.1 and 1 µM). The AUC fold-change ranged from 1.1 to 3.6. In summary, t-CA is a potential source of pharmacokinetic variability for CYP2A6 substrates due to metabolism-dependent inhibition, especially in scenarios when exposure to t-CA is elevated due to high dietary exposure, or when cinnamon is used as a treatment of specific disease states (e.g., diabetes).

Active Site Mutations as a Suitable Tool Contributing to Explain a Mechanism of Aristolochic Acid I Nitroreduction by Cytochromes P450 1A1, 1A2 and 1B1

Aristolochic acid I (AAI) is a plant drug found in Aristolochia species that causes aristolochic acid nephropathy, Balkan endemic nephropathy and their associated urothelial malignancies. AAI is activated via nitroreduction producing genotoxic N-hydroxyaristolactam, which forms DNA adducts. The major enzymes responsible for the reductive bioactivation of AAI are

NAD(P)H

quinone oxidoreductase and cytochromes P450 (CYP) 1A1 and 1A2. Using site-directed mutagenesis we investigated the possible mechanisms of CYP1A1/1A2/1B1-catalyzed AAI nitroreduction. Molecular modelling predicted that the hydroxyl groups of serine122/threonine124 (Ser122/Thr124) amino acids in the CYP1A1/1A2-AAI binary complexes located near to the nitro group of AAI, are mechanistically important as they provide the proton required for the stepwise reduction reaction. In contrast, the closely related CYP1B1 with no hydroxyl group containing residues in its active site is ineffective in catalyzing AAI nitroreduction. In order to construct an experimental model, mutant forms of CYP1A1 and 1A2 were prepared, where Ser122 and Thr124 were replaced by Ala (CYP1A1-S122A) and Val (CYP1A2-T124V), respectively. Similarly, a CYP1B1 mutant was prepared in which Ala133 was replaced by Ser (CYP1B1-A133S). Site-directed mutagenesis was performed using a quickchange approach. Wild and mutated forms of these enzymes were heterologously expressed in Escherichia coli and isolated enzymes characterized using UV-vis spectroscopy to verify correct protein folding. Their catalytic activity was confirmed with CYP1A1, 1A2 and 1B1 marker substrates. Using (32)P-postlabelling we determined the efficiency of wild-type and mutant forms of CYP1A1, 1A2, and 1B1 reconstituted with NADPH:CYP oxidoreductase to bioactivate AAI to reactive intermediates forming covalent DNA adducts. The S122A and T124V mutations in CYP1A1 and 1A2, respectively, abolished the efficiency of CYP1A1 and 1A2 enzymes to generate AAI-DNA adducts. In contrast, the formation of AAI-DNA adducts was catalyzed by CYP1B1 with the A133S mutation. Our experimental model confirms the importance of the hydroxyl group possessing amino acids in the active center of CYP1A1 and 1A2 for AAI nitroreduction.

Kinetic analysis of electron flux in cytochrome P450 reductases reveals differences in rate-determining steps in plant and mammalian enzymes

Herein, we compare the kinetic properties of CPR from Arabidopsis thaliana (ATR2), with CPR from Artemisia annua (aaCPR) and human CPR (hCPR). While all three CPR forms elicit comparable rates for cytochrome c(3+) turnover, NADPH reduction of the FAD cofactor is ∼50-fold faster in aaCPR and ATR2 compared to hCPR, with a kobs of ∼500 s(-1) (6 °C). Stopped-flow analysis of the isolated FAD-domains reveals that NADP(+)-FADH2 charge-transfer complex formation is also significantly faster in the plant enzymes, but the rate of its decay is comparable for all three proteins. In hCPR, transfer of a hydride ion from NADPH to FAD is tightly coupled to subsequent FAD to FMN electron transfer, indicating that the former catalytic event is slow relative to the latter. In contrast, interflavin electron transfer is slower than NADPH hydride transfer in aaCPR and ATR2, occurring with an observed rate constant of ∼50 s(-1). Finally, the transfer of electrons from FMN to cytochrome c(3+) is rapid (>10(3) s(-1)) in all three enzymes and does not limit catalytic turnover. In combination, the data reveal differences in rate-determining steps between plant CPR and their mammalian equivalent in mediating the flux of reducing equivalents from NADPH to external electron acceptors.

Anion-Dependent Stimulation of CYP3A4 Monooxygenase

We co-crystallized human cytochrome P450 3A4 (CYP3A4) with progesterone (PRG) under two different conditions, but the resulting complexes contained only one PRG molecule bound to the previously identified peripheral site. A novel feature in one of our structures is a citrate ion, originating from the crystallization solution. The citrate-binding site is located in an area where the N-terminus splits from the protein core and, thus, is suitable for the interaction with the anionic phospholipids of the microsomal membrane. We investigated how citrate affects the function of a soluble CYP3A4 monooxygenase system consisting of equimolar amounts of CYP3A4 and cytochrome P450 reductase (CPR). Citrate was found to affect the properties of both redox partners and stimulated their catalytic activities in a concentration-dependent manner via a complex mechanism. CYP3A4-substrate binding, reduction of CPR with NADPH, and interflavin and interprotein electron transfer were identified as citrate-sensitive steps. Comparative analysis of various negatively charged organic compounds indicated that, in addition to alterations caused by changes in ionic strength, anions modulate the properties of CYP3A4 and CPR through specific anion-protein interactions. Our results help to better understand previous observations and provide new mechanistic insights into CYP3A4 function.

The role of cytochrome b5 structural domains in interaction with cytochromes P450

To understand the role of the structural elements of cytochrome b5 in its interaction with cytochrome P450 and the catalysis performed by this heme protein, we carried out comparative structural and functional analysis of the two major mammalian forms of membrane-bound cytochrome b5 - microsomal and mitochondrial, designed chimeric forms of the heme proteins in which the hydrophilic domain of one heme protein is replaced by the hydrophilic domain of another one, and investigated the effect of the highly purified native and chimeric heme proteins on the enzymatic activity of recombinant cytochromes P4503A4 and P45017A1 (CYP3A4 and CYP17A1). We show that the presence of a hydrophobic domain in the structure of cytochrome b5 is necessary for its effective interaction with its redox partners, while the nature of the hydrophobic domain has no significant effect on the ability of cytochrome b5 to stimulate the activity of cytochrome P450-catalyzed reactions. Thus, the functional properties of cytochrome b5 are mainly determined by the structure of the heme-binding domain.

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.

Mechanism-based inactivation of human cytochrome P450 3A4 by two piperazine-containing compounds

Human cytochrome P450 3A4 (CYP3A4) is responsible for the metabolism of more than half of pharmaceutic drugs, and inactivation of CYP3A4 can lead to adverse drug-drug interactions. The substituted imidazole compounds 5-fluoro-2-[4-[(2-phenyl-1H-imidazol-5-yl)methyl]-1-piperazinyl]pyrimidine (SCH 66712) and 1-[(2-ethyl-4-methyl-1H-imidazol-5-yl)methyl]-4-[4-(trifluoromethyl)-2-pyridinyl]piperazine (EMTPP) have been previously identified as mechanism-based inactivators (MBI) of CYP2D6. The present study shows that both SCH 66712 and EMTPP are also MBIs of CYP3A4. Inhibition of CYP3A4 by SCH 66712 and EMTPP was determined to be concentration, time, and NADPH dependent. In addition, inactivation of CYP3A4 by SCH 66712 was shown to be unaffected by the presence of electrophile scavengers. SCH 66712 displays type I binding to CYP3A4 with a spectral binding constant (Ks) of 42.9 ± 2.9 µM. The partition ratios for SCH 66712 and EMTPP were 11 and 94, respectively. Whole protein mass spectrum analysis revealed 1:1 binding stoichiometry of SCH 66712 and EMTPP to CYP3A4 and a mass increase consistent with adduction by the inactivators without addition of oxygen. Heme adduction was not apparent. Multiple mono-oxygenation products with each inactivator were observed; no other products were apparent. These are the first MBIs to be shown to be potent inactivators of both CYP2D6 and CYP3A4.

Functional characterization of steroid hydroxylase CYP106A1 derived from Bacillus megaterium

In this study, we examined the catalytic activity of CYP106A1 from the Bacillus megaterium American Type Culture Collection 14581 strain. The CYP106A1 gene was cloned from B. megaterium, heterologously expressed in Escherichia coli, and purified. Potential electron partners and possible bacterial CYP106A1 substrates were identified by examining the oxidative activity toward a set of steroids in the presence of several reductase systems. The activities of CYP106A1 in a reconstituted system could not be achieved using rat NADPH-P450 reductase or a putidaredoxin reductase-putidaredoxin pair. However, the spinach redox proteins, a ferredoxin reductase-ferredoxin pair, were found to be efficient redox partners for CYP106A1. CYP106A1 catalyzes the hydroxylation of a set of steroids including testosterone, progesterone, 17α-hydroxyprogesterone, 11-deoxycorticosterone, corticosterone, and 11-deoxycortisol to produce monohydroxylated products as the major metabolites. These results suggest that CYP106A1 would be useful for the bioconversion of steroid hormones to hydroxylated products that can be used for industrial applications.

Crystal structure of cindoxin, the P450cin redox partner

The crystal structure of the flavin mononucleotide (FMN)-containing redox partner to P450cin, cindoxin (Cdx), has been determined to 1.3 Å resolution. The overall structure is similar to that of the FMN domain of human cytochrome P450 reductase. A Brownian dynamics–molecular dynamics docking method was used to produce a model of Cdx with its redox partner, P450cin. This Cdx–P450cin model highlights the potential importance of Cdx Tyr96 in bridging the FMN and heme cofactors as well P450cin Arg102 and Arg346. Each of the single-site Ala mutants exhibits ∼10% of the wild-type activity, thus demonstrating the importance of these residues for binding and/or electron transfer. In the well-studied P450cam system, redox partner binding stabilizes the open low-spin conformation of P450cam and greatly decreases the stability of the oxy complex. In sharp contrast, Cdx does not shift P450cin to a low-spin state, although the stability of oxy-P450cin is decreased 10-fold in the presence of Cdx. This indicates that Cdx may have a modest effect on the open–closed equilibrium in P450cin compared to that in P450cam. It has been postulated that part of the effector role of Pdx on P450cam is to promote a significant structural change that makes available a proton relay network involving Asp251 required for O₂ activation. The structure around the corresponding Asp in P450cin, Asp241, provides a possible structural reason for why P450cin is less dependent on its redox partner for functionally important structural changes.

Genomic and bioinformatic analysis of NADPH-cytochrome P450 reductase in Anopheles stephensi (Diptera: Culicidae)

The cytochrome P450 monooxygenase (P450) enzyme system is a major mechanism of xenobiotic biotransformation. The nicotinamide adenine dinucleotide phosphate (NADPH)-cytochrome P450 reductase (CPR) is required for transfer of electrons from NADPH to P450. One CPR gene was identified in the genome of the malaria-transmitting mosquito Anopheles stephensi Liston (Diptera: Culicidae). The gene encodes a polypeptide containing highly conserved flavin mononucleotide-, flavin adenine dinucleotide-, and NADPH-binding domains, a unique characteristic of the reductase. Phylogenetic analysis revealed that the A. stephensi and other known mosquito CPRs belong to a monophyletic group distinctly separated from other insects in the same order, Diptera. Amino acid residues of CPRs involved in binding of P450 and cytochrome c are conserved between A. stephensi and the Norway rat Rattus norvegicus Berkenhout (Rodentia: Muridae). However, gene structure particularly within the coding region is evidently different between the two organisms. Such difference might arise during the evolution process as also seen in the difference of P450 families and isoforms found in these organisms. CPR in the mosquito A. stephensi is expected to be active and serve as an essential component of the P450 system.

Cytochrome P450-dependent catabolism of vitamin K: ω-hydroxylation catalyzed by human CYP4F2 and CYP4F11

Vitamin K plays an essential role in many biological processes including blood clotting, maintenance of bone health, and inhibition of arterial calcification. A menaquinone form of vitamin K, MK4, is increasingly recognized for its key roles in mitochondrial electron transport, as a ligand for the nuclear receptor SXR, which controls the expression of genes involved in transport and metabolism of endo- and xenobiotics, and as a pharmacotherapeutic in the treatment of osteoporosis. Although cytochrome P450 (CYP) 4F2 activity is recognized as an important determinant of phylloquinone (K1) metabolism, the enzymes involved in menaquinone catabolism have not been studied previously. CYP4F2 and CYP4F11 were expressed and purified and found to be equally efficient as in vitro catalysts of MK4 ω-hydroxylation. CYP4F2, but not CYP4F11, catalyzed sequential metabolism of MK4 to the ω-acid without apparent release of the intermediate aldehyde. The ω-alcohol could also be metabolized to the acid by microsomal NAD(+)-dependent alcohol and aldehyde dehydrogenases. LC-MS/MS analysis of trypsinized human liver microsomes (using a surrogate peptide approach) revealed the mean concentrations of CYP4F2 and CYP4F11 to be 14.3 and 8.4 pmol/mg protein, respectively. Microsomal MK4 ω-hydroxylation activities correlated with the CYP4F2 V433M genotype but not the CYP4F11 D446N genotype. Collectively, these data expand the lexicon of vitamin K ω-hydroxylases to include the 'orphan' P450 CYP4F11 and identify a common variant, CYP4F2 (rs2108622), as a major pharmacogenetic variable influencing MK4 catabolism.