Electrostatic interaction between NADH-cytochrome b5 reductase and cytochrome b5 studied by site-directed mutagenesis

Cytochrome b5 reductases: Redox regulators of cell homeostasis

The cytochrome-b5 reductase (CYB5R) family of flavoproteins is known to regulate reduction-oxidation (redox) balance in cells. The five enzyme members are highly compartmentalized at the subcellular level and function as "redox switches" enabling the reduction of several substrates, such as heme and coenzyme Q. Critical insight into the physiological and pathophysiological significance of CYB5R enzymes has been gleaned from several human genetic variants that cause congenital disease and a broad spectrum of chronic human diseases. Among the CYB5R genetic variants, CYB5R3 is well-characterized and deficiency in expression and activity is associated with type II methemoglobinemia, cancer, neurodegenerative disorders, diabetes, and cardiovascular disease. Importantly, pharmacological and genetic-based strategies are underway to target CYB5R3 to circumvent disease onset and mitigate severity. Despite our knowledge of CYB5R3 in human health and disease, the other reductases in the CYB5R family have been understudied, providing an opportunity to unravel critical function(s) for these enzymes in physiology and disease. In this review, we aim to provide the broad scientific community an up-to-date overview of the molecular, cellular, physiological, and pathophysiological roles of CYB5R proteins.

Transmembrane helices mediate the formation of a stable ternary complex of b5R, cyt b5, and SCD1

Mammalian cytochrome b₅ (cyt b₅) and cytochrome b₅ reductase (b₅R) are electron carrier proteins for membrane-embedded oxidoreductases. Both b₅R and cyt b₅ have a cytosolic domain and a single transmembrane (TM) helix. The cytosolic domains of b₅R and cyt b₅ contain cofactors required for electron transfer, but it is not clear if the TM helix has function beyond being an anchor to the membrane. Here we show that b₅R and cyt b₅ form a stable binary complex, and so do cyt b₅ and stearoyl-CoA desaturase-1 (SCD1). We also show that b₅R, cyt b₅ and SCD1 form a stable ternary complex. We demonstrate that the TM helices are required for the assembly of stable binary and ternary complexes where electron transfer rates are greatly enhanced. These results reveal a role of the TM helix in cyt b₅ and b₅R, and suggest that an electron transport chain composed of a stable ternary complex may be a general feature in membrane-embedded oxidoreductases that require cyt b₅ and b₅R.

The transmembrane domains of mammalian cytochrome b₅ (cyt b₅), cyt b₅ reductase (b₅R), and stearoyl-CoA desaturase-1 (SCD1) form stable binary complexes between cyt b₅/b₅R or cyt b₅/SCD1 and a ternary complex, which enhance electron transfer rates.

Mutation update: Variants of the CYB5R3 gene in recessive congenital methemoglobinemia

NADH-cytochrome b5 reductase 3 deficiency is an important genetic cause of recessive congenital methemoglobinemia (RCM) and occurs worldwide in autosomal recessive inheritance. In this Mutation Update, we provide a comprehensive review of all the pathogenic mutations and their molecular pathology in RCM along with the molecular basis of RCM in 21 new patients from the Indian population, including four novel variants: c.103A>C (p.Thr35Pro), c.190C>G (p.Leu64Val), c.310G>T (p.Gly104Cys), and c.352C>T (p.His118Tyr). In this update, over 78 different variants have been described for RCM globally. Molecular modeling of all the variants reported in CYB5R3 justifies association with the varying severity of the disease. The majority of the mutations associated with the severe form with a neurological disorder (RCM Type 2) were associated with the FAD-binding domain of the protein while the rest were located in another domain of the protein (RCM Type 1).

Crystal structures of the naturally fused CS and cytochrome b5 reductase (b5R) domains of Ncb5or reveal an expanded CS fold, extensive CS-b5R interactions and productive binding of the NAD(P)+ nicotinamide ring

Ncb5or (NADH-cytochrome b5 oxidoreductase), a cytosolic ferric reductase implicated in diabetes and neurological diseases, comprises three distinct domains, cytochrome b5 (b5) and cytochrome b5 reductase (b5R) domains separated by a CHORD-Sgt1 (CS) domain, and a novel 50-residue N-terminal region. Understanding how interdomain interactions in Ncb5or facilitate the shuttling of electrons from NAD(P)H to heme, and how the process compares with the microsomal b5 (Cyb5A) and b5R (Cyb5R3) system, is of interest. A high-resolution structure of the b5 domain (PDB entry 3lf5) has previously been reported, which exhibits substantial differences in comparison to Cyb5A. The structural characterization of a construct comprising the naturally fused CS and b5R domains with bound FAD and NAD+ (PDB entry 6mv1) or NADP+ (PDB entry 6mv2) is now reported. The structures reveal that the linker between the CS and b5R cores is more ordered than predicted, with much of it extending the β-sandwich motif of the CS domain. This limits the flexibility between the two domains, which recognize one another via a short β-sheet motif and a network of conserved side-chain hydrogen bonds, salt bridges and cation-π interactions. Notable differences in FAD-protein interactions in Ncb5or and Cyb5R3 provide insight into the selectivity for docking of their respective b5 redox partners. The structures also afford a structural explanation for the unusual ability of Ncb5or to utilize both NADH and NADPH, and represent the first examples of native, fully oxidized b5R family members in which the nicotinamide ring of NAD(P)+ resides in the active site. Finally, the structures, together with sequence alignments, show that the b5R domain is more closely related to single-domain Cyb5R proteins from plants, fungi and some protists than to Cyb5R3 from animals.

Structural and enzymatic analysis of the cytochrome b5 reductase domain of Ulva prolifera nitrate reductase

Rapid accumulations of unattached green macroalgae, referred to as blooms, constitute ecological disasters and occur in many coastal regions. Ulva are a major cause of blooms, owing to their high nitrogen utilization capacity, which requires nitrate reductase (NR) activity; however, molecular characterization of Ulva NR remains lacking. Herein we determined the crystal structure and performed an enzymatic analysis of the cytochrome b₅ reductase domain of Ulva prolifera NR (UpCbRNR). The structural analysis revealed an N-terminal FAD-binding domain primarily consisting of six antiparallel β strands, a C-terminal NADH-binding domain forming a Rossmann fold, and a three β-stranded linker region connecting these two domains. The FAD cofactor was located in the cleft between the two domains and interacted primarily with the FAD-binding domain. UpCbRNR shares similarities in overall structure and cofactor interactions with homologs, and its catalytic ability is comparable to that of higher plant CbRNRs. Structure and sequence comparisons of homologs revealed two regions of sequence length variation potentially useful for phylogenetic analysis: one in the FAD-binding domain, specific to U. prolifera, and another in the linker region that may be used to differentiate between plant, fungi, and animal homologs. Our data will facilitate molecular-level understanding of nitrate assimilation in Ulva.

Topography of human cytochrome b5/cytochrome b5 reductase interacting domain and redox alterations upon complex formation

Cytochrome b₅ is the main electron acceptor of cytochrome b₅ reductase. The interacting domain between both human proteins has been unidentified up to date and very little is known about its redox properties modulation upon complex formation. In this article, we characterized the protein/protein interacting interface by solution NMR and molecular docking. In addition, upon complex formation, we measured an increase of cytochrome b₅ reductase flavin autofluorescence that was dependent upon the presence of cytochrome b_5. Data analysis of these results allowed us to calculate a dissociation constant value between proteins of 0.5±0.1μM and a 1:1 stoichiometry for the complex formation. In addition, a 30mV negative shift of cytochrome b₅ reductase redox potential in presence of cytochrome b₅ was also measured. These experiments suggest that the FAD group of cytochrome b₅ reductase increase its solvent exposition upon complex formation promoting an efficient electron transfer between the proteins.

Role of cytochrome b5 in the modulation of the enzymatic activities of cytochrome P450 17α-hydroxylase/17,20-lyase (P450 17A1)

Cytochrome b5 (cyt b5) is a small hemoprotein that plays a significant role in the modulation of activities of an important steroidogenic enzyme, cytochrome P450 17α-hydroxylase/17,20-lyase (P450 17A1, CYP17A1). Located in the zona fasciculata and zona reticularis of the adrenal cortex and in the gonads, P450 17A1 catalyzes two different reactions in the steroidogenic pathway; the 17α-hydroxylation and 17,20-lyase, in the endoplasmic reticulum of these respective tissues. The activities of P450 17A1 are regulated by cyt b5 that enhances the 17,20-lyase reaction by promoting the coupling of P450 17A1 and cytochrome P450 reductase (CPR), allosterically. Cyt b5 can also act as an electron donor to enhance the 16-ene-synthase activity of human P450 17A1. In this review, we discuss the many roles of cyt b5 and focus on the modulation of CYP17A1 activities by cyt b5 and the mechanisms involved.

Contribution of Electrostatics to the Kinetics of Electron Transfer from NADH-Cytochrome b5 Reductase to Fe(III)-Cytochrome b5

Brownian dynamics (BD) simulations provide here a theoretical atomic-level treatment of the reduction of human ferric cytochrome b5 (cyt b5) by NADH-cytochrome b5 reductaste (cyt b5r) and several of its mutants. BD is used to calculate the second-order rate constant of electron transfer (ET) between the proteins for direct correlation with experiments. Interestingly, the inclusion of electrostatic forces dramatically increases the reaction rate of the native proteins despite the overall negative charge of both proteins. The role played by electrostatic charge distribution in stabilizing the ET complexes and the role of mutations of several amino acid residues in stabilizing or destabilizing the complexes are analyzed. The complex with the shortest ET reaction distance (d = 6.58 Å) from rigid body BD is further subjected to 1 ns of molecular dynamics (MD) in a periodic box of TIP3P water to produce a more stable complex allowed by flexibility and with a shorter average reaction distance d = 6.02 Å. We predict a docking model in which the following ion-ion interactions are dominant (cyt b5r/cyt b5): Lys162-Heme O1D/Lys163-Asp64/Arg91-Heme O1A/Lys125-Asp70.

Human cytochrome b5 reductase: structure, function, and potential applications

Cytochrome b5 reductase is a flavoprotein that is produced as two different isoforms that have different localizations. The amphipathic microsomal isoform, found in all cell types with the exception of erythrocytes, consists of one hydrophobic membrane-anchoring domain and a larger hydrophilic flavin catalytic domain. The soluble cytochrome b5 reductase isoform, found in human erythrocytes, is a truncated protein that is encoded by an alternative transcript and consists of the larger domain only. Cytochrome b5 reductase is involved in the transfer of reducing equivalents from the physiological electron donor, NADH, via an FAD domain to the small molecules of cytochrome b5. This protein has received much attention from researchers due to its involvement in many oxidation and reduction reactions, such as the reduction of methemoglobin to hemoglobin. Autosomal cytochrome b5 reductase gene deficiency manifests with the accumulation of oxidized Fe+3 and recessive congenital methemoglobinemia in humans. In this article, we provide a comprehensive overview of the structure and function of cytochrome b5 reductase from different eukaryotic sources and its potential use in the food industry, biosensor, and diagnostic areas.

Stimulation and clustering of cytochrome b5 reductase in caveolin-rich lipid microdomains is an early event in oxidative stress-mediated apoptosis of cerebellar granule neurons

The apoptosis of cerebellar granule neurons (CGN) induced by low potassium in the extracellular medium is a model of neuronal apoptosis where an overshot of reactive oxygen species (ROS) triggers the neuronal death. In this work, using dihydroethidium and L-012 as specific dyes for superoxide anion detection we show that this ROS overshot can be accounted by an increased release of superoxide anion to the extracellular medium. The amplitude and time course of the increase of superoxide anion observed early during apoptosis correlated with the increase of the content of soluble cytochrome b(5), a substrate of the NADH-dependent oxidase activity of the cytochrome b(5) reductase associated with lipid rafts in CGN. Western blotting and immunofluorescence microscopy approaches, including fluorescence energy transfer, pointed out an enhanced clustering of cytochrome b(5) reductase within caveolins-rich lipid rafts microdomains. Protein/protein docking analysis suggests that cytochrome b(5) reductase can form complexes with caveolins 1α, 1β and 2, playing electrostatic interactions a major role in this association. In conclusion, our results indicate that overstimulation of cytochrome b(5) reductase associated with lipid rafts can account for the overshot of plasma membrane-focalized superoxide anion production that triggers the entry of CGN in the irreversible phase of apoptosis. This article is part of a Special Issue entitled: Proteomics: The clinical link.

Study of the individual cytochrome b5 and cytochrome b5 reductase domains of Ncb5or reveals a unique heme pocket and a possible role of the CS domain

NADH cytochrome b(5) oxidoreductase (Ncb5or) is found in animals and contains three domains similar to cytochrome b(5) (b(5)), CHORD-SGT1 (CS), and cytochrome b(5) reductase (b(5)R). Ncb5or has an important function, as suggested by the diabetes and lipoatrophy phenotypes in Ncb5or null mice. To elucidate the structural and functional properties of human Ncb5or, we generated its individual b(5) and b(5)R domains (Ncb5or-b(5) and Ncb5or-b(5)R, respectively) and compared them with human microsomal b(5) (Cyb5A) and b(5)R (Cyb5R3). A 1.25 Å x-ray crystal structure of Ncb5or-b(5) reveals nearly orthogonal planes of the imidazolyl rings of heme-ligating residues His(89) and His(112), consistent with a highly anisotropic low spin EPR spectrum. Ncb5or is the first member of the cytochrome b(5) family shown to have such a heme environment. Like other b(5) family members, Ncb5or-b(5) has two helix-loop-helix motifs surrounding heme. However, Ncb5or-b(5) differs from Cyb5A with respect to location of the second heme ligand (His(112)) and of polypeptide conformation in its vicinity. Electron transfer from Ncb5or-b(5)R to Ncb5or-b(5) is much less efficient than from Cyb5R3 to Cyb5A, possibly as a consequence of weaker electrostatic interactions. The CS linkage probably obviates the need for strong interactions between b(5) and b(5)R domains in Ncb5or. Studies with a construct combining the Ncb5or CS and b(5)R domains suggest that the CS domain facilitates docking of the b(5) and b(5)R domains. Trp(114) is an invariant surface residue in all known Ncb5or orthologs but appears not to contribute to electron transfer from the b(5)R domain to the b(5) domain.

Mechanism of electron transfer in fusion protein cytochrome b5-NADH-cytochrome b5 reductase

In the present work we summarize results on construction of expression plasmid, heterologous expression in Escherichia coli, isolation and purification, as well as physicochemical characterization of chimeric protein consisting of hydrophilic domain of cytochrome b(5) and truncated from the N-terminal sequence (Delta(23)) form of NADH-cytochrome b(5) reductase. The kinetics and mechanism of electron transfer between NADH-cytochrome b(5) reductase and cytochrome b(5) in the frames of fusion protein consisting of cytochrome b(5) (94 amino acids) and truncated form of NADH-cytochrome b(5) reductase (277 amino acids) have been studied. It is shown that electron transfer takes place between redox partners belonging to two different molecules of the chimeric protein. Using computer modeling, we built the model of the tertiary structure of the fusion protein, which is in agreement with experimental data. By using Marcus theory of electron transfer in polar media, we demonstrate the inability of the hypothesis of electrostatic repulsions to explain the increase of electron transfer rate on increase of ion concentration in media due to elimination of the repulsion of similar charges. The real reason for the increase of the first order rate constant in some oxidation-reduction reactions between proteins, as shown in the present work, is a decrease of the media reorganization energy resulting in decrease of activation energy for oxidation-reduction reactions.

Molecular and enzymatic characterizations of novel bifunctional 3beta-hydroxysteroid dehydrogenases/C-4 decarboxylases from Arabidopsis thaliana

We have isolated two cDNAs from Arabidopsis thaliana encoding bifunctional 3beta-hydroxysteroid dehydrogenase/C-4 decarboxylases (3betaHSD/D) involved in sterol synthesis, termed At3betaHSD/D1 and At3betaHSD/D2. Transformation of the yeast ergosterol auxotroph erg26 mutant, which lacks 3betaHSD/D activity, with the At3betaHSD/D1 isoform or with At3betaHSD/D2 isoform containing a C-terminal At3betaHSD/D1 endoplasmic reticulum-retrieval sequence restored growth and ergosterol synthesis in erg26. An in vitro enzymatic assay revealed high 3betaHSD/D activity for both isoenzymes in the corresponding microsomal extracts. The two At3betaHSD/D isoenzymes showed similar substrate specificities that required free 3beta-hydroxyl and C-4-carboxyl groups but were quite tolerant in terms of variations of the sterol nucleus and side chain structures. Data obtained with 4alpha-carboxy-cholest-7-en-3beta-ol and its 3alpha-deuterated analog revealed that 3alpha-hydrogen-carbon bond cleavage is not the rate-limiting step of the reaction. In planta reduction on the expression of the 3betaHSD/D gene as a consequence of VIGS-mediated gene silencing in Nicotiana benthamiana led to a substantial accumulation of 3beta-hydroxy-4beta,14-dimethyl-5alpha-ergosta-9beta,19-cyclo-24(24(1))-en-4alpha-carboxylic acid, consistent with a decrease in 3betaHSD/D activity. These two novel oxidative decarboxylases constitute the first molecularly and functionally characterized HSDs from a short chain dehydrogenase/reductase family in plants.

Human cytochrome b5 requires residues E48 and E49 to stimulate the 17,20-lyase activity of cytochrome P450c17

Cytochrome P450c17 (CYP17) catalyzes both the 17alpha-hydroxylase and 17,20-lyase reactions in human steroid biosynthesis. Cytochrome b5 (b5) stimulates the rate of the 17,20-lyase reaction 10-fold with little influence on 17alpha-hydroxylase activity. Studies with apo-b5 suggest that stimulation of 17,20-lyase activity results from an allosteric action on the hCYP17 x POR complex, rather than electron transfer by b5. We hypothesized that specific residues on b5 interact with the hCYP17 x POR complex and that targeted mutation of surface-exposed residues might identify b5 residues critical for stimulating 17,20-lyase activity. We constructed, expressed, and purified 14 single plus 3 double b5 mutations and assayed their ability to stimulate 17,20-lyase activity. Most mutations did not alter the capacity of b5 to stimulate 17,20-lyase activity or appeared to modestly alter the affinity of b5 for the hCYP17 x POR complex. In contrast, mutation of E48, E49, or R52 reduced the maximal stimulation of 17,20-lyase activity. In particular, b5 mutation E48G + E49G lost over 95% of the capacity to stimulate 17,20-lyase activity, yet this mutation retained normal electron transfer properties. In addition, mutation E48G + E49G did not impair stimulation of 17,20-lyase activity by wild-type b5, suggesting that the mutation binds poorly to the site of the hCYP17 x POR complex occupied by b5. These data suggest that a specific allosteric binding site on b5, which includes residues E48, E49, and possibly R52, mediates the stimulation of 17,20-lyase activity.

The structure and biochemistry of NADH-dependent cytochrome b5 reductase are now consistent

Cytochrome b5 reductase (cb5r) (EC 1.6.6.2) catalyzes the reduction of two molecules of cytochrome b5 using NADH as the physiological electron donor. The structure of pig cb5r at 2.4 A resolution was previously reported in the literature, but it was inconsistent with the biochemistry; for example, K83 and C245 were both implicated in the mechanism, but were not located at the active site. To address this problem, we have determined the structures of cb5r from rat at 2.0 A resolution and in a complex with NAD+ at 2.3 A resolution. We found significant differences throughout the rat structure compared to that of pig, including the locations of the lysine and cysteine residues mentioned above. To test the structural models, we made single amino acid substitutions of this lysine and showed that all substitutions produced correctly folded proteins and exhibited normal flavin behavior. However, the apparent kcat(NADH) decreased, and the apparent K(m) for NADH increased; the K(m)'s for cytochrome b5 were unchanged relative to that of the wild type. The largest effect was for the glutamate-substituted protein, which was further characterized using a charge transfer assay and found to be less efficient at NADH utilization than the wild type. These results are consistent with a role for this lysine in stabilizing the NADH-bound form of cb5r. We have concluded that the pig structure was mistraced in several regions and have reinterpreted mutants in these regions that give rise to the hereditary disease methemoglobinemia.

Seven new mutations in the nicotinamide adenine dinucleotide reduced-cytochrome b(5) reductase gene leading to methemoglobinemia type I

Cytochrome b(5) reductase (b5R) deficiency manifests itself in 2 distinct ways. In methemoglobinemia type I, the patients only suffer from cyanosis, whereas in type II, the patients suffer in addition from severe mental retardation and neurologic impairment. Biochemical data indicate that this may be due to a difference in mutations, causing enzyme instability in type I and complete enzyme deficiency or enzyme inactivation in type II. We have investigated 7 families with methemoglobulinemia type I and found 7 novel mutations in the b5R gene. Six of these mutations predicted amino acid substitutions at sites not involved in reduced nicotinamide adenine dinucleotide (NADH) or flavin adenine dinucleotide (FAD) binding, as deduced from a 3-dimensional model of human b5R. This model was constructed from comparison with the known 3-dimensional structure of pig b5R. The seventh mutation was a splice site mutation leading to skipping of exon 5 in messenger RNA, present in heterozygous form in a patient together with a missense mutation on the other allele. Eight other amino acid substitutions, previously described to cause methemoglobinemia type I, were also situated in nonessential regions of the enzyme. In contrast, 2 other substitutions, known to cause the type II form of the disease, were found to directly affect the consensus FAD-binding site or indirectly influence NADH binding. Thus, these data support the idea that enzyme inactivation is a cause of the type II disease, whereas enzyme instability may lead to the type I form.

Systematic mutations of highly conserved His49 and carboxyl-terminal of recombinant porcine liver NADH-cytochrome b5 reductase solubilized domain

The cDNA encoding solubilized porcine liver NADH-cytochrome b5 reductase catalytic domain (Pb5R) was cloned and overexpressed in Escherichia coli. A highly conserved His49 and a C-terminal Phe272 of Pb5R, which are located near the isoalloxazine moiety of the FAD, were systematically modulated by site-directed mutagenesis. Large structural change was not detected on the absorption and circular dichroism spectra of mutant proteins. Drastic changes in enzymatic properties were not observed, but the apparent Km value for soluble form of porcine liver cytochrome b5 (Pb5) was affected by the substitutions of His49 with glutamic acid and with lysine, deletion of C-terminal Phe272, and addition of Gly273. The values of the catalytic constant (kcat) were obviously decreased by the substitution of His49 with glutamic acid or the addition of Gly273. In these two mutants, the rate for reduction of FAD was decreased, and the rate for autoxidation of reduced FAD was increased. These results showed that His49 and C-terminal carboxyl group in Pb5R are not critical for the electron transfer to Pb5, but the electrostatic environmental changes at these positions could affect the recognition of Pb5 and modulate the catalytic function of the enzyme by changing the stability of reduced FAD.

Characterization of the second form of NADH-Cytochrome b5 reductase gene from arachidonic acid-producing fungus Mortierella alpina 1S-4

The second type of cytochrome b5 reductase (Cb5R-II) gene was characterized in the arachidonic acid-producing fungus Mortierella alpina 1S-4. Its cDNA (897 bp) and predicted amino acid (298 aa) sequences show more than 70% similarity to the previously isolated first type of Cb5R. Highly conserved exon-intron organization suggests that the two genes evolved from the duplication of a common ancestral gene. Cb5R-II has a flavin-binding domain at its highly hydrophobic N-terminal and an NADH-binding domain at the C-terminal. In comparison with Cb5R genes from other sources, high homology (46-54%) was found for yeast and plant genes. Phylogenetic analysis revealed that microbial and plant Cb5R genes represent a gene family evolved from one prototype and are different from mammalian Cb5R genes.