Mass spectrometry identification of cytochrome P450 2B4 interaction sites for NADPH: Cytochrome P450 reductase

The interaction sites for protein partners, cytochrome P450 2B4 (d-2B4) and NADPH: cytochrome P450 reductase (d-Fp), have been identified. These proteins form complexes during their functioning. Nonspecific covalent cross-linking of the d-2B4 complexes with d-Fp in the Emulgen 913 monomerized system was achieved by 4,4′- dithiobis-phenyl azide. Covalently cross-linked peptides of this complex were identified by ESI-MS/MS. Several binding sites have been identified for these proteins. Based on these sites a model for intermolecular interaction between these proteins has been proposed. This model includes 5 contact sites on d-2B4 for d-Fp (stabilized by the cross-linker); these include the following pairs of corresponding peptides of d-2B4 and d-Fp: 1) d-2B4_324–336 and d-Fp_570–585; 2) d-2B4_423–433 and d-Fp_102–109; 3) d-2B4_327–336 and d Fp_452–464; 4) d-2B4_192–197 and d-Fp_456–464; 5) d-2B4_134–139 and d-Fp_406–425. In the two last cases d-Fp peptides are located in the interdomain cleft and stabilize the protein-protein complex via the cross-linker and so the d 2B4/d-Fp complex formation by these sites may involve amino acid residues of the peptides d-Fp_456–464 and d-Fp_406–425, which surround the interdomain cleft.     

Involvement of NADPH in the interaction between heme oxygenase-1 and cytochrome P450 reductase

Heme oxygenase-1 (HO-1) catalyzes the physiological degradation of heme at the expense of molecular oxygen using electrons donated by NADPH-cytochrome P450 reductase (CPR). In this study, we investigated the effect of NADP(H) on the interaction of HO-1 with CPR by surface plasmon resonance. We found that HO-1 associated with CPR more tightly in the presence of NADP(+) (K(D) = 0.5 microm) than in its absence (K(D) = 2.4 microm). The HO-1 mutants, K149A, K149A/K153A, and R185A, showed almost no heme degradation activity with NADPH-CPR, whereas they exhibited activity comparable to that of the wild type when sodium ascorbate was used. R185A showed a 100-fold decreased affinity for CPR compared with wild type, even in the presence of NADP(+) (K(D) = 36.3 microm). The affinities of K149A and K149A/K153A for CPR were decreased 7- and 9-fold (K(D) = 16.8 and 21.8 microm), respectively. In contrast to R185A, the affinities of K149A and K149A/K153A were improved by the addition of NADP(+) (K(D) = 5.2 and 9.6 microm, respectively), as was the case with wild type. Computer modeling of the HO-1/CPR complex showed that the guanidino group of Arg(185) is located within the hydrogen bonding distance of 2'-phosphate of NADPH, suggesting that Arg(185) contributes to the binding to CPR through an electrostatic interaction with the phosphate group. On the other hand, Lys(149) is close to a cluster of acidic amino acids near the FMN binding site of CPR. Thus, Lys(149) and Lys(153) appear to interact with CPR in such a way as to orient the redox partners for optimal electron transfer from FMN of CPR to heme of HO-1.     

Interflavin electron transfer in human cytochrome P450 reductase is enhanced by coenzyme binding. Relaxation kinetic studies with coenzyme analogues.

The role of coenzyme binding in regulating interflavin electron transfer in human cytochrome P450 reductase (CPR) has been studied using temperature-jump spectroscopy. Previous studies [Gutierrez, A., Paine, M., Wolf, C.R., Scrutton, N.S., & Roberts, G.C.K. Biochemistry (2002) 41, 4626-4637] have shown that the observed rate, 1/tau, of interflavin electron transfer (FADsq - FMNsq-->FADox - FMNhq) in CPR reduced at the two-electron level with NADPH is 55 +/- 2 s-1, whereas with dithionite-reduced enzyme the observed rate is 11 +/- 0.5 s-1, suggesting that NADPH (or NADP+) binding has an important role in controlling the rate of internal electron transfer. In relaxation experiments performed with CPR reduced at the two-electron level with NADH, the observed rate of internal electron transfer (1/tau = 18 +/- 0.7 s-1) is intermediate in value between those seen with dithionite-reduced and NADPH-reduced enzyme, indicating that the presence of the 2'-phosphate is important for enhancing internal electron transfer. To investigate this further, temperature jump experiments were performed with dithionite-reduced enzyme in the presence of 2',5'-ADP and 2'-AMP. These two ligands increase the observed rate of interflavin electron transfer in two-electron reduced CPR from 1/tau = 11 s-1 to 35 +/- 0.2 s-1 and 32 +/- 0.6 s-1, respectively. Reduction of CPR at the two-electron level by NADPH, NADH or dithionite generates the same spectral species, consistent with an electron distribution that is equivalent regardless of reductant at the initiation of the temperature jump. Spectroelectrochemical experiments establish that the redox potentials of the flavins of CPR are unchanged on binding 2',5'-ADP, supporting the view that enhanced rates of interdomain electron transfer have their origin in a conformational change produced by binding NADPH or its fragments. Addition of 2',5'-ADP either to the isolated FAD-domain or to full-length CPR (in their oxidized and reduced forms) leads to perturbation of the optical spectra of both the flavins, consistent with a conformational change that alters the environment of these redox cofactors. The binding of 2',5'-ADP eliminates the unusual dependence of the observed flavin reduction rate on NADPH concentration (i.e. enhanced at low coenzyme concentration) observed in stopped-flow studies. The data are discussed in the context of previous kinetic studies and of the crystallographic structure of rat CPR.     

Differential effect of phenobarbital and beta-naphthoflavone on the mRNAs coding for cytochrome P450 and NADPH cytochrome P450 reductase

The induction in rat liver of a specific variant(s) of cytochrome P450 (PB-P450) by phenobarbital and its repression by β-naphthoflavone occur through corresponding changes in the levels of mRNA coding for the protein(s). The level of translatable mRNA coding for NADPH-cytochrome P450 reductase in rat liver increases on treatment with phenobarbital but not β-naphthoflavone.     

Conformational dynamics and the energetics of protein--ligand interactions: role of interdomain loop in human cytochrome P450 reductase

A combination of mutagenesis, calorimetry, kinetics, and small-angle X-ray scattering (SAXS) has been used to study the mechanism of ligand binding energy propagation through human cytochrome P450 reductase (CPR). Remarkably, the energetics of 2',5'-ADP binding to R597 at the FAD-binding domain are affected by mutations taking place at an interdomain loop located 60 A away. Either deletion of a 7 amino acid long segment (T236-G237-E238-E239-S240-S241-I242) or its replacement by poly-proline repeats (5 and 10 residues) results in a significant increase in 2',5'-ADP enthalpy of binding (DeltaHB). This is accompanied by a decrease in the number of thermodynamic microstates available for the ligand-CPR complex. Moreover, the estimated heat capacity change (DeltaCp) for this interaction changes from -220 cal mol-1 K-1 in the wild-type enzyme to -580 cal mol-1 K-1 in the deletion mutant. Pre-steady-state kinetics measurements reveal a 50-fold decrease in the microscopic rate for interdomain (FAD --> FMN) electron transfer in the deletion mutant (kobs = 0.4 s-1). Multiple turnover cytochome c reduction assays indicate that these mutations impair the ability of the FMN-binding domain to shuttle electrons from the FAD-binding domain to the cytochrome partner. Binding of 2',5'-ADP to wild-type CPR triggers a large-scale structural rearrangement resulting in the complex having a more compact domain organization, and the maximum molecular dimension (Dmax) decreases from 110 A in ligand-free enzyme to 100 A in the ligand-bound CPR. The SAXS experiments also demonstrate that what is affected by the mutations is indeed the relative diffusional motion of the domains. Furthemore, ab initio shape reconstruction and homology modeling would suggest that-in the deletion mutant-hindering of domain motion occurs concomitantly with dimerization. The results presented here show that the energetics of this highly localized interaction (2',5'-ADP binding) have a global character, and are highly sensitive to functional structural dynamics involving distal domains. These findings support early theoretical studies which postulate a single protein molecule to be a real, independent thermodynamic ensemble.     

Structure of the open conformation of a functional chimeric NADPH cytochrome P450 reductase

Two catalytic domains, bearing FMN and FAD cofactors, joined by a connecting domain, compose the core of the NADPH cytochrome P450 reductase (CPR). The FMN domain of CPR mediates electron shuttling from the FAD domain to cytochromes P450. Together, both enzymes form the main mixed‐function oxidase system that participates in the metabolism of endo‐ and xenobiotic compounds in mammals. Available CPR structures show a closed conformation, with the two cofactors in tight proximity, which is consistent with FAD‐to‐FMN, but not FMN‐to‐P450, electron transfer. Here, we report the 2.5 Å resolution crystal structure of a functionally competent yeast–human chimeric CPR in an open conformation, compatible with FMN‐to‐P450 electron transfer. Comparison with closed structures shows a major conformational change separating the FMN and FAD cofactors from 86 Å.     

Conformational studies of NADPH cytochrome P-450 reductase by circular dichroism: interaction with phospholipids

Ultraviolet circular dichroism spectrum of purified NADPH cytochrome P-450 reductase was characterized by two negative bands centered at 208 and 222 nm. The approximation of the alpha-helical content from the value of the mean residue ellipticity at 222 nm indicated 28% of alpha-helical structures. Heat inactivation of the enzyme was associated to a drastic change in the secondary structure of the protein. Membrane reconstitution experiments by inclusion of the enzyme into liposomes revealed that the conformation of NADPH cytochrome P-450 reductase was sensitive to its phospholipid environment. Egg lecithin as well as synthetic phosphatidylcholines, at the optimal phospholipid-enzyme molar ratio 200, was able to increase up to 37% the mean residue ellipticity at 222 nm. Addition of phosphatidylserine or phosphatidylethanolamine produced no effect. Non-ionic detergent such as Emulgen 913 weakly enhanced the mean residue ellipticity.     

The role of cytochrome P450 lysine residues in the interaction between cytochrome P450IA1 and NADPH-cytochrome P450 reductase.

Cytochrome P450IA1 (purified from hepatic microsomes of beta-naphthoflavone-treated rats) has been covalently modified with the lysine-modifying reagent acetic anhydride. Different levels of lysine residue modification in cytochrome P450IA1 can be achieved by varying the concentration of acetic anhydride. Modification of lysine residues in P450IA1 greatly inhibits the interaction of P450IA1 with NADPH-cytochrome P450 reductase. Modification of 1.0 and 3.3 mol lysine residues per mole P450IA1 resulted in 30 and 95% decreases, respectively, in 7-ethoxycoumarin hydroxylation by a reconstituted P450IA1/reductase complex. However, modification of 3.3 mol lysine residues per mole P450IA1 decreased only cumene hydroperoxide-supported P450-dependent 7-ethoxycoumarin hydroxylation by 30%. Spectral and fluorescence studies showed no indication of global conformational change of P450IA1 even with up to 8.8 mol lysine residues modified per mole P450IA1. These data suggest that at least three lysine residues in P450IA1 may be involved in the interaction with reductase. Identification of lysine residues in P450IA1 possibly involved in this interaction was carried out by [14C]acetic anhydride modification, trypsin digestion, HPLC separation, and amino acid sequencing. The lysine residue candidates identified in this manner were K97, K271, K279, and K407.     

Cross-linking of human cytochrome P450 2B6 to NADPH-cytochrome P450 reductase: Identification of a potential site of interaction

The site(s) of interaction between human cytochrome P450 2B6 and NADPH-cytochrome P450 reductase (P450 reductase) have yet to be identified. To investigate this, the cross-linking agent 1-ethyl-3-[3-dimethylaminopropyl]carbodiimide hydrochloride (EDC) was used to covalently link P450 2B6-P450 reductase. Following digestion with trypsin, the cross-linked peptides were identified by reconstituting the peptides in ¹⁸O-water based on the principle that the cross-linked peptides would be expected to incorporate twice as many ¹⁸O atoms as the non-cross-linked peptides. Subsequent mass spectrometric analyses of the resulting peptides led to the identification of one cross-linked peptide candidate. De novo sequencing of the peptide indicated that it is a complex between residues in the C-helix of the P450 (based upon solved X-ray crystal structures of P450 2B4) and the connecting domain of the P450 reductase. To confirm this experimentally, the P450 2B6 peptide identified through the cross-linking studies was synthesized and peptide competition studies were performed. In the presence of the synthetic peptide, P450 catalytic activity was decreased by up to 60% when compared to competition studies performed using a nonsense peptide. Taken together, these studies indicate that residues in the C-helix of P450 2B6 play a major role in the interaction with the P450 reductase.     

Altered heme catabolism by heme oxygenase-1 caused by mutations in human NADPH cytochrome P450 reductase

Human heme oxygenase-1 (HO-1) carries out heme catabolism supported by electrons supplied from the NADPH through NADPH P450 reductase (POR, CPR). Previously we have shown that mutations in human POR cause a rare form of congenital adrenal hyperplasia. In this study, we have evaluated the effects of mutations in POR on HO-1 activity. We used purified preparations of wild type and mutant human POR and in vitro reconstitution with purified HO-1 to measure heme degradation in a coupled assay using biliverdin reductase. Here we show that mutations in POR found in patients may reduce HO-1 activity, potentially influencing heme catabolism in individuals carrying mutant POR alleles. POR mutants Y181D, A457H, Y459H, V492E and R616X had total loss of HO-1 activity, while POR mutations A287P, C569Y and V608F lost 50-70% activity. The POR variants P228L, R316W and G413S, A503V and G504R identified as polymorphs had close to WT activity. Loss of HO-1 activity may result in increased oxidative neurotoxicity, anemia, growth retardation and iron deposition. Further examination of patients affected with POR deficiency will be required to assess the metabolic effects of reduced HO-1 activity in affected individuals.     

Escherichia coli MTC, a NADPH cytochrome P450 reductase competent mutagenicity tester strain for the expression of human cytochrome P450: comparison of three types of expression systems

Currently three different methods have been taken to develop new mutagenicity tester strains containing human cytochrome P450s (CYPs). Each of these use a single expression vector. In this paper we describe a fourth approach, i.e., the coexpression of a CYP and its electron-transfer flavoprotein, NADPH CYP reductase (RED), encoded by two different expression vectors. The Escherichia coli mutagenicity tester strain BMX100 has been expanded to a strain, MTC which stably expresses human RED. This new tester strain permits the biplasmid coexpression of human CYP1A2 and RED (MTC1A2). This novel strain can be used for the determination of the mutagenicity of chemicals known to be procarcinogens and metabolized by CYP1A2. The mutagenicity tester strain MTC1A2 was compared with: (i) BMX100 using the post-mitochondrial rat liver fraction (S9); (ii) BMX100 with expressing CYP1A2 alone (iii) or with expressing CYP1A2 fused to rat RED or (iv) with expressing CYP1A2, bicistronically coexpressed with rat RED. The biplasmid RED/CYP coexpression system generated a strain with the highest methoxy- and ethoxy-resorufin dealkylase activities and the highest mutagenic activities for the procarcinogens 2-aminoanthracene (2AA), aflatoxin B1 (AFB1) and 2-amino-3-methylimidazo(4,5-f)quinoline (IQ). Furthermore, the metabolism of 2AA and IQ was detected more efficiently using the MTC1A2 strain than with the BMX100 strain plus the standard rodent liver S9 metabolic system.     

Examining the mechanism of stimulation of cytochrome P450 by cytochrome b5: the effect of cytochrome b5 on the interaction between cytochrome P450 2B4 and P450 reductase

Dissociation constants K(d) for cytochrome P450 reductase (reductase) and cytochrome P450 2B4 are measured in the presence of various substrates. Aminopyrine increases the dissociation constant for binding of the two proteins. Furthermore, cytochrome b(5) (b(5)) stimulates metabolism of this substrate and dramatically decreases the substrate-related K(d) values. Experiments are performed to test if the b(5)-mediated stimulation is effected through a conformational change of P450. The effects of a redox-inactive analogue of b(5) (Mn b(5)) on product formation and reaction stoichiometry are determined. Variations in the concentration of Mn b(5) stock solution that have been shown to effect the aggregation state of the protein alter the rate of P450-mediated NADPH oxidation but have no effect on the rate of product formation. Thus, the electron transfer capability of b(5) is necessary for stimulation of metabolism. Furthermore, stopped flow spectrometry measurements of the rate of first electron reduction of the P450 by reductase indicate that the coupling of P450 2B4-mediated metabolism improves, in the presence of Mn b(5), with slower delivery of the first electron of the catalytic cycle by the reductase. These results are consistent with a model involving the regulation of the P450 catalytic cycle by conformational changes of the P450 enzyme. We propose that the conformational change(s) necessary for progression of the catalytic cycle is inhibited when reduced, but not oxidized, reductase is bound to the P450.     

Electrostatic interaction between cytochrome P450 and NADPH-P450 reductase: comparison of mixed and fused systems consisting of rat cytochrome P450 1A1 and yeast NADPH-P450 reductase

The electrostatic interaction between rat cytochrome P450 1A1 and yeast NADPH-P450 reductase was analyzed by using recombinant yeast microsomes containing both native enzymes or their fused enzyme. The Vmax of the 7-ethoxycoumarin O-deethylation in the recombinant microsomes containing both rat cytochrome P4501A1 and yeast NADPH-P450 reductase (the mixed system) was maximal when the ionic strength of the reaction mixture was 0.1-0.15. However, on the fused enzyme between rat cytochrome P450 1A1 and yeast NADPH-P450 reductase (the fused system), the activity was uniformly reduced with increasing ionic strength. The pH profiles of Vmax were also different between the mixed and the fused systems. Based on these results, we propose a hypothesis that cytochrome P450 and NADPH-P450 reductase have more than one binding mode. The maximal activity of the mixed system at ionic strength of 0.1-0.15 is explained by change of the binding mode. On the other hand, the fused enzyme appears to have only one binding mode due to the limited topology of cytochrome P450 and NADPH-P450 reductase domains.     

Relaxation kinetics of cytochrome P450 reductase: internal electron transfer is limited by conformational change and regulated by coenzyme binding

The kinetics of internal electron transfer in human cytochrome P450 reductase have been studied using temperature-jump relaxation spectroscopy. Temperature perturbation of CPR reduced at the two-electron level with NADPH yields biphasic absorption transients at 450 and 600 nm. The observed rate, 1/tau, for the fast phase is 2200 +/- 300 s(-1). The absence of this phase in fluorescence transients and in absorption transients collected with dithionite-reduced enzyme indicates this phase does not report on electron/hydride transfer and is consistent with its origin in local conformational change in the vicinity of the FAD isoalloxazine ring. The slow phase (1/tau = 55 +/- 2 s(-1)) observed in the absorption transients obtained with CPR reduced at the two-electron level with NADPH reports on internal electron transfer: FAD(sq)-FMN(sq) --> FAD(ox)-FMN(hq). The observed rate of this transient is slower (1/tau = 11 +/- 0.5 s(-1)) in CPR reduced to the two-electron level by dithionite rather than NADPH, demonstrating that coenzyme binding has an important influence on the observed rate of internal electron transfer. Temperature perturbation experiments with CPR reduced with 10-fold molar excess of NADPH produce monophasic absorption transients (1/tau = 20 +/- 0.2 s(-1)) reporting on internal electron transfer: FAD(sq)-FMN(hq) --> FAD(hq)-FMN(sq). The observed rate constants for electron transfer are substantially less than those expected from analysis of CPR by electron-transfer theory (approximately 10(10) s(-1)). Potential gating mechanisms have been investigated using the temperature-jump method. Observed rates for electron transfer were unaffected in experiments performed in deuterated solvent, indicating that deprotonation does not gate the reaction. Introduction of glycerol into the sample significantly decreased the observed rate for internal electron transfer, suggesting conformational gating of the reaction. Replacement of Trp-676 with His-676 reduces approximately 2-fold the observed rate of internal electron transfer in two-electron-reduced enzyme, whereas the observed rate for FAD(sq)-FMN(hq) --> FAD(hq)-FMN(sq) transfer is increased approximately 13-fold in the W676H mutant reduced with a 10-fold molar excess of NADPH. The studies reveal altered redox properties of the FAD in W676H CPR. The data are discussed in the context of previous stopped-flow studies of human CPR and the X-ray crystallographic structure of rat CPR.     

High diversity and complex evolution of fungal cytochrome P450 reductase: cytochrome P450 systems

Cytochrome P450 reductase (CPR) is the redox partner of P450 monooxygenases, involved in primary and secondary metabolism of eukaryotes. Two novel CPR genes, sharing 34% amino acid identity, were found in the filamentous ascomycete Cochliobolus lunatus. Fungal genomes were searched for putative CPR enzymes. Phylogenetic analysis suggests that multiple independent CPR duplication events occurred in fungi, whereas P450-CPR fusion occurred before the diversification of Dikarya and Zygomycota. Additionally, losses of methionine synthase reductase were found in certain fungal taxa; a truncated form of this enzyme was conserved in Pezizomycotina. In fungi, high numbers of cytochrome P450 enzymes, multiple CPRs, and P450-CPR fusion proteins were associated with filamentous growth. Evolution of multiple CPR-like oxidoreductases in filamentous fungi might have been driven by the complexity of biochemical functions necessitated by their growth form, as opposed to yeast.     

Differential effect of phenobarbital and β-naphthoflavone on the mRNAs coding for cytochrome P450 and NADPH cytochrome P450 reductase

The induction in rat liver of a specific variant(s) of cytochrome P450 (PB-P450) by phenobarbital and its repression by β-naphthoflavone occur through corresponding changes in the levels of mRNA coding for the protein(s). The level of translatable mRNA coding for NADPH-cytochrome P450 reductase in rat liver increases on treatment with phenobarbital but not β-naphthoflavone.     

Expression and characterization of full-length human heme oxygenase-1: the presence of intact membrane-binding region leads to increased binding affinity for NADPH cytochrome P450 reductase

Heme oxygenase-1 (HO-1) is the chief regulatory enzyme in the oxidative degradation of heme to biliverdin. In the process of heme degradation, HO-1 receives the electrons necessary for catalysis from the flavoprotein NADPH cytochrome P450 reductase (CPR), releasing free iron and carbon monoxide. Much of the recent research involving heme oxygenase has been done using a 30 kDa soluble form of the enzyme, which lacks the membrane binding region (C-terminal 23 amino acids). The goal of this study was to express and purify a full-length human HO-1 (hHO-1) protein; however, due to the lability of the full-length form, a rapid purification procedure was required. This was accomplished by use of a glutathione-s-transferase (GST)-tagged hHO-1 construct. Although the procedure permitted the generation of a full-length HO-1, this form was contaminated with a 30 kDa degradation product that could not be eliminated. Therefore, attempts were made to remove a putative secondary thrombin cleavage site by a conservative mutation of amino acid 254, which replaces arginine with lysine. This mutation allowed the expression and purification of a full-length hHO-1 protein. Unlike wild type (WT) HO-1, the R254K mutant could be purified to a single 32 kDa protein capable of degrading heme at the same rate as the WT enzyme. The R254K full-length form had a specific activity of approximately 200-225 nmol of bilirubin h-1 nmol-1 HO-1 as compared to approximately 140-150 nmol of bilirubin h-1 nmol-1 for the WT form, which contains the 30 kDa contaminant. This is a 2-3-fold increase from the previously reported soluble 30 kDa HO-1, suggesting that the C-terminal 23 amino acids are essential for maximal catalytic activity. Because the membrane-spanning domain is present, the full-length hHO-1 has the potential to incorporate into phospholipid membranes, which can be reconstituted at known concentrations, in combination with other endoplasmic reticulum resident enzymes.     

Microsomal lipid peroxidation: the role of NADPH--cytochrome P450 reductase and cytochrome P450

The role of NADPH--cytochrome P450 reductase and cytochrome P450 in NADPH- and ADP--Fe3(+)-dependent lipid peroxidation was investigated by using the purified enzymes and liposomes prepared from either total rat-liver phospholipids or a mixture of bovine phosphatidyl choline and phosphatidyl ethanolamine (PC/PE liposomes). The results suggest that NADPH- and ADP--Fe3(+)-dependent lipid peroxidation involves both NADPH--cytochrome P450 reductase and cytochrome P450. Just as in the case of cytochrome P450-linked monooxygenations, the role of these enzymes in lipid peroxidation may be to provide two electrons for O2 reduction. The first electron is used for reduction of ADP--Fe3+ and subsequent addition of O2 to the perferryl radical (ADP--Fe3(+)-O2-), which then extracts an H atom from a polyunsaturated lipid (LH) giving rise to a free radical (LH.) that reacts with O2 yielding a peroxide free radical (LOO.). The second electron is then used to reduce LOO. to the lipid hydroperoxide (LOOH). In the latter capacity, reduced cytochrome P450 can be replaced by EDTA--Fe2+ or by the superoxide radical as generated through redox cycling of a quinone such as menadione.     

High diversity and complex evolution of fungal cytochrome P450 reductase: cytochrome P450 systems.

Cytochrome P450 reductase (CPR) is the redox partner of P450 monooxygenases, involved in primary and secondary metabolism of eukaryotes. Two novel CPR genes, sharing 34% amino acid identity, were found in the filamentous ascomycete Cochliobolus lunatus. Fungal genomes were searched for putative CPR enzymes. Phylogenetic analysis suggests that multiple independent CPR duplication events occurred in fungi, whereas P450-CPR fusion occurred before the diversification of Dikarya and Zygomycota. Additionally, losses of methionine synthase reductase were found in certain fungal taxa; a truncated form of this enzyme was conserved in Pezizomycotina. In fungi, high numbers of cytochrome P450 enzymes, multiple CPRs, and P450-CPR fusion proteins were associated with filamentous growth. Evolution of multiple CPR-like oxidoreductases in filamentous fungi might have been driven by the complexity of biochemical functions necessitated by their growth form, as opposed to yeast.