Models for the complexes formed between cytochrome b5 and the subunits of methemoglobin

Redox-Regulation of α-Globin in Vascular Physiology

Interest in the structure, function, and evolutionary relations of circulating and intracellular globins dates back more than 60 years to the first determination of the three-dimensional structure of these proteins. Non-erythrocytic globins have been implicated in circulatory control through reactions that couple nitric oxide (NO) signaling with cellular oxygen availability and redox status. Small artery endothelial cells (ECs) express free α-globin, which causes vasoconstriction by degrading NO. This reaction converts reduced (Fe2+) α-globin to the oxidized (Fe3+) form, which is unstable, cytotoxic, and unable to degrade NO. Therefore, (Fe3+) α-globin must be stabilized and recycled to (Fe2+) α-globin to reinitiate the catalytic cycle. The molecular chaperone α-hemoglobin-stabilizing protein (AHSP) binds (Fe3+) α-globin to inhibit its degradation and facilitate its reduction. The mechanisms that reduce (Fe3+) α-globin in ECs are unknown, although endothelial nitric oxide synthase (eNOS) and cytochrome b5 reductase (CyB5R3) with cytochrome b5 type A (CyB5a) can reduce (Fe3+) α-globin in solution. Here, we examine the expression and cellular localization of eNOS, CyB5a, and CyB5R3 in mouse arterial ECs and show that α-globin can be reduced by either of two independent redox systems, CyB5R3/CyB5a and eNOS. Together, our findings provide new insights into the regulation of blood vessel contractility.

Methemoglobin formation in mutant hemoglobin α chains: electron transfer parameters and rates

Hemoglobin-mediated transport of dioxygen (O₂) critically depends on the stability of the reduced (Fe²⁺) form of the heme cofactors. Some protein mutations stabilize the oxidized (Fe³⁺) state (methemoglobin, Hb M), causing methemoglobinemia, and can be lethal above 30%. The majority of the analyses of factors influencing Hb oxidation are retrospective and give insights only for inner-sphere mutations of heme (His58, His87). Herein, we report the first all-atom molecular dynamics simulations on both redox states and calculations of the Marcus electron transfer (ET) parameters for the α chain Hb oxidation and reduction rates for Hb M. The Hb wild-type (WT) and most of the studied α chain variants maintain globin structure except the Hb M Iwate (H87Y). The mutants forming Hb M tend to have lower redox potentials and thus stabilize the oxidized (Fe³⁺) state (in particular, the Hb Miyagi variant with K61E mutation). Solvent reorganization (λ_solv 73-96%) makes major contributions to reorganization free energy, whereas protein reorganization (λ_prot) accounts for 27-30% except for the Miyagi and J-Buda variants (λ_prot ∼4%). Analysis of heme-solvent H-bonding interactions among variants provide insights into the role of Lys61 residue in stabilizing the Fe²⁺ state. Semiclassical Marcus ET theory-based calculations predict experimental k_ET for the Cyt b5-Hb complex and provide insights into relative reduction rates for Hb M in Hb variants. Thus, our methodology provides a rationale for the effect of mutations on the structure, stability, and Hb oxidation reduction rates and has potential for identification of mutations that result in methemoglobinemia.

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.

Cytochrome b₅ coexpression increases Tetrahymena thermophila Δ6 fatty acid desaturase activity in Saccharomyces cerevisiae

Very-long-chain polyunsaturated fatty acids such as arachidonic, eicosapentaenoic, and docosahexaenoic acids, are important to the physiology of many microorganisms and metazoans and are vital to human development and health. The production of these and related fatty acids depends on Δ6 desaturases, the final components of an electron transfer chain that introduces double bonds into 18-carbon fatty acid chains. When a Δ6 desaturase identified from the ciliated protist Tetrahymena thermophila was expressed in Saccharomyces cerevisiae cultures supplemented with the 18:2(Δ9,12) substrate, only 4% of the incorporated substrate was desaturated. Cytochrome b₅ protein sequences identified from the genome of T. thermophila included one sequence with two conserved cytochrome b₅ domains. Desaturation by the Δ6 enzyme increased as much as 10-fold when T. thermophila cytochrome b₅s were coexpressed with the desaturase. Coexpression of a cytochrome b₅ from Arabidopsis thaliana with the Δ6 enzyme also increased desaturation. A split ubiquitin growth assay indicated that the strength of interaction between cytochrome b₅ proteins and the desaturase plays a vital role in fatty acid desaturase activity, illustrating the importance of protein-protein interactions in this enzyme activity.

Surface electrochemistry and electrocatalytic activity of ion pairs formed from oppositely-charged phthalocyanine and porphyrin species

Self-assembled ion pairs are deposited onto a graphite electrode surface by dipping the electrode sequentially into phthalocyanine cations and anions, or a phthalocyanine cation and porphyrin anion. The electrochemical properties of these surface-bound ion pairs are compared with the properties of the individual surface-bound components. The order of deposition is shown to be crucial, implying that these species lie flat on the surface. A detailed study of the electrocatalytic oxidation of oxalate anion is used to explore their reactivity as a function of deposition order.

The nitrite anion binds to human hemoglobin via the uncommon O-nitrito mode

The nitrite anion is known to oxidize and degrade hemoglobin (Hb). Recent literature reports suggest a nitrite reductase activity for Hb, converting nitrite into nitric oxide. Surprisingly, no structural information about Hb-nitrite interactions has been reported. We have determined the crystal structure of the ferric Hb-nitrite complex at 1.80 A resolution. The nitrite ligand adopts the uncommon O-nitrito binding mode. In addition, the nitrito conformations in the alpha and beta subunits are different, reflecting subtle effects of the distal His in orienting the nitrite ligand in the O-nitrito binding mode.

Phospholipid membranes affect tertiary structure of the soluble cytochrome b5 heme-binding domain

The influence of charged phospholipid membranes on the conformational state of the water-soluble fragment of cytochrome b5 has been investigated by a variety of techniques at neutral pH. The results of this work provide the first evidence that aqueous solutions with high phospholipid/protein molar ratios (pH 7.2) induce the cytochrome to undergo a structural transition from the native conformation to an intermediate state with molten-globule like properties that occur in the presence of an artificial membrane surface and that leads to binding of the protein to the membrane. At other phospholipid/protein ratios, equilibrium was observed between cytochrome free in solution and cytochrome bound to the surface of vesicles. Inhibition of protein binding to the vesicles with increasing ionic strength indicated for the most part an electrostatic contribution to the stability of cytochrome b5-vesicle interactions at pH 7.2. The possible physiological role of membrane-induced conformational change in the structure of cytochrome b5 upon the interaction with its redox partners is discussed.

Unique structure of Ascaris suum b5-type cytochrome: an additional alpha-helix and positively charged residues on the surface domain interact with redox partners

Cytochrome b5 of the body wall of adult Ascaris suum, a porcine parasitic nematode, is a soluble protein that lacks a C-terminal membrane-anchoring domain, but possesses an N-terminal pre-sequence of 30 amino acids. During the maturation of cytochrome b5, the N-terminal pre-sequence is proteolytically cleaved to form the mature protein of 82 amino acid residues. A. suum cytochrome b5 is a basic protein containing more lysine residues and exhibiting a higher midpoint redox potential than its mammalian counterparts. We developed an expression system for the production of the recombinant nematode cytochrome b5, which is chemically and functionally identical with the native protein. Using this recombinant protein, we have determined the X-ray crystal structure of A. suum cytochrome b5 at 1.8 A (1 A=0.1 nm) resolution, and we have shown that this protein is involved in the reduction of nematode body-wall metmyoglobin. The crystal structure of A. suum cytochrome b5 consists of six alpha-helices and five beta-strands. It differs from its mammalian counterparts by having a head-to-tail disulphide bridge, as well as a four-residue insertion in the vicinity of the sixth ligating histidine, which forms an additional alpha-helix, alpha4A, between helices alpha4 and alpha5. A. suum cytochrome b5 exists predominantly as a haem-orientation B isomer. Furthermore, the haem plane is rotated approx. 80 degrees relative to the axis formed by haem-Fe and N atoms of the two histidine residues that are ligated to haem-Fe. The charge distribution around the haem crevice of A. suum cytochrome b5 is remarkably different from that of mammalian cytochrome b5 in that the nematode protein bears positively charged lysine residues surrounding the haem crevice. Using immunohistochemistry, we found that A. suum cytochrome b5 is present in the nematode hypodermis. Based on this histochemical and structural information, the physiological function of A. suum cytochrome b5 and its interaction with nematode metmyoglobin can be hypothesized.

Characterization and calculation of a cytochrome c-cytochrome b5 complex using NMR data

To identify the binding site for bovine cytochrome b(5) (cyt b(5)) on horse cytochrome c (cyt c), cross-saturation transfer NMR experiments were performed with (2)H- and (15)N-enriched cyt c and unlabeled cyt b(5). In addition, chemical shift changes of the cyt c backbone amide and side chain methyl resonances were monitored as a function of cyt b(5) concentration. The chemical shift changes indicate that the complex is in fast exchange, and are consistent with a 1:1 stoichiometry. A K(a) of (4 +/- 3) x 10(5) M(-)(1) was obtained with a lower limit of 855 s(-)(1) for the dissociation rate of the complex. Mapping of the chemical shift variations and intensity changes upon cross-saturation NMR experiments in the complex reveals a single, contiguous interaction interface on cyt c. Using NMR data as constraints, a protein docking program was used to calculate two low-energy model complex clusters. Independent calculations of the effect of the cyt b(5) heme ring current-induced magnetic dipole on cyt c were used to discriminate between the different models. The interaction surface of horse cyt c in the current experimentally constrained model of the cyt c-cyt b(5) complex is similar but not identical to the interface predicted in yeast cyt c by Brownian dynamics and docking calculations. The occurrence of different amino acids at the protein-protein interface and the dissimilar assumptions employed in the calculations can largely account for the nonidentical interfaces.

Mapping the binding interface of the cytochrome b5-cytochrome c complex by nuclear magnetic resonance

The interaction between bovine cytochrome b(5) (cyt b(5)) and horse heart cytochrome c (cyt c) is investigated by NMR spectroscopy. Chemical shifts of cyt b(5) backbone resonances and side chain methyl resonances were monitored as a function of cyt c concentration. The shifts are small but saturatable and indicate that the binding of cyt b(5) with cyt c is in fast exchange. An equilibrium association constant of (6 +/- 3) x 10(4) M(-1) was obtained with a lower limit of 180 s(-1) for the dissociation rate of the complex. To resolve considerable ambiguities in the interpretation of the chemical shift mapping, (15)N relaxation experiments and cross-saturation experiments were used as alternative methods to map the cyt b(5)-cyt c binding interface. Results from the three experiments combined demonstrate that the conserved negatively charged region of cyt b(5) surrounding the solvent-exposed heme edge is involved in the interaction with cyt c. These data support the models proposed by Salemme and Mauk [(1976) J. Mol. Biol. 102, 563-568; (1993) Biochemistry 32, 6613-6623].

Functional interaction of cytochrome P450 with its redox partners: a critical assessment and update of the topology of predicted contact regions

The problem of donor-acceptor recognition has been the most important and intriguing one in the area of P450 research. The present review outlines the topological background of electron-transfer complex formation, showing that the progress in collaborative investigations, combining physical techniques with chemical-modification and immunolocalization studies as well as site-directed mutagenesis experiments, has increasingly enabled the substantiation of hypothetical work resulting from homology modelling of P450s. Circumstantial analysis reveals the contact regions for redox proteins to cluster on the proximal face of P450s, constituting parts of the highly conserved, heme-binding core fold. However, more variable structural components located in the periphery of the hemoprotein molecules also participate in donor docking. The cross-reactivity of electron carriers, purified from pro- and eukaryotic sources, with a diversity of P450 species points at a possible evolutionary conservation of common anchoring domains. While electrostatic mechanisms appear to dominate orientation toward each other of the redox partners to generate pre-collisional encounter complexes, hydrophobic forces are likely to foster electron transfer events by through-bonding or pi-stacking interactions. Moreover, electron-tunneling pathways seem to be operative as well. The availability of new P450 crystal structures together with improved validation strategies will undoubtedly permit the production of increasingly satisfactory three-dimensional donor-acceptor models serving to better understand the molecular principles governing functional association of the redox proteins.

X-ray crystallography, CD and kinetic studies revealed the essence of the abnormal behaviors of the cytochrome b5 Phe35-->Tyr mutant

Conserved phenylalanine 35 is one of the hydrophobic patch residues on the surface of cytochrome b5 (cyt b5). This patch is partially exposed on the surface of cyt b5 while its buried face is in direct van der Waals' contact with heme b. Residues Phe35 and Phe/Tyr74 also form an aromatic channel with His39, which is one of the axial ligands of heme b. By site-directed mutagenesis we have produced three mutants of cyt b5: Phe35-->Tyr, Phe35-->Leu, and Phe35-->His. We found that of these three mutants, the Phe35-->Tyr mutant displays abnormal properties. The redox potential of the Phe35-->Tyr mutant is 66 mV more negative than that of the wild-type cyt b5 and the oxidized Phe35-->Tyr mutant is more stable towards thermal and chemical denaturation than wild-type cyt b5. In this study we studied the most interesting mutant, Phe35-->Tyr, by X-ray crystallography, thermal denaturation, CD and kinetic studies of heme dissociation to explore the origin of its unusual behaviors. Analysis of crystal structure of the Phe35-->Tyr mutant shows that the overall structure of the mutant is basically the same as that of the wild-type protein. However, the introduction of a hydroxyl group in the heme pocket, and the increased van der Waals' and electrostatic interactions between the side chain of Tyr35 and the heme probably result in enhancement of stability of the Phe35-->Tyr mutant. The kinetic difference of the heme trapped by the heme pocket also supports this conclusion. The detailed conformational changes of the proteins in response to heat have been studied by CD for the first time, revealing the existence of the folding intermediate.

Electron transfer between cytochrome b(5) and some oxidised haemoglobins: the role of ionic strength

We have compared experimental measurements and Brownian dynamic calculations for the reduction of oxidised adult human haemoglobin by reduced bovine cytochrome b(5) over a range of ionic strengths. Our calculations suggest that the presence of molecular electrostatic fields have a significant role to play in the formation of the electron transfer complexes. These results predict that electron transfer occurs within an ensemble of similarly weakly docked complexes, the formation of which is strongly ionic strength dependent. Application of electron tunneling analysis to the complexes allows us to predict the rates of electron transfer within each ensemble of complexes as a function of ionic strength. The outcome of this theoretical study is compared with experimental rate measurements. A comparison of the results obtained from adult and embryonic haemoglobins, at a fixed ionic strength, indicates a significant difference in the characteristics of complex formation. These data emphasise the role played by electrostatic interactions in this important physiological reaction.

Determination of the hemoglobin surface domains that react with cytochrome b5

We have compared the photoinitiated electron-transfer (ET) reaction between cytochrome b(5) (b(5)) and zinc mesoporphyrin-substituted hemoglobin [(ZnM)Hb] and Hb variants in order to determine whether b(5) binds to the subunit surface of either or both Hb chains, or to sites which span the dimer--dimer interface. Because the dimer--dimer interface would be disrupted for monomers or alpha beta dimers, we studied the reaction of b(5) with alpha ZnM chains and (ZnM)Hb beta W37E, which exists as alpha beta dimers in solution. Triplet quenching titrations of the ZnHb proteins with Fe(3+)b(5) show that the binding affinity and ET rate constants for the alpha-chains are the same when they are incorporated into a Hb tetramer or dimer, or exist as monomers. Likewise, the parameters for beta-chains in tetramers and dimers differ minimally. In parallel, we have modified the surface of the Hb chains by neutralizing the heme propionates through the preparation of zinc deuterioporphyrin dimethyl ester hemoglobin, (ZnD-DME)Hb. The charge neutralization increases the ET rate constants 100-fold for the alpha-chains and 40-fold for the beta-chains (but has has little effect on the affinity of either chain type for b(5), similar to earlier results for myoglobin). Together, these results indicate that b(5) binds to sites at the subunit surface of each chain rather than to sites which span the dimer-dimer interface. The charge-neutralization results further suggest that b(5) binds over a broad area of the subunit face, but reacts only in a minority population of binding geometries.

Three-dimensional structure of the nonaheme cytochrome c from Desulfovibrio desulfuricans Essex in the Fe(III) state at 1.89 A resolution

A nine heme group containing cytochrome c isolated from the soluble and membrane fractions of Desulfovibrio desulfuricans Essex, termed nonaheme cytochrome c, was crystallized, and the structure was solved using the multiple wavelength anomalous dispersion (MAD) phasing method. Refinement was carried out to a resolution of 1.89 A, and anisotropic temperature factors were addressed to the iron and sulfur atoms in the model. The structure revealed two cytochrome c(3) like domains with the typical arrangement of four heme centers. Both domains flanked an extra heme buried under the protein surface. This heme is held in position by loop extensions in each of the two domains. Although both the N- and C-terminal tetraheme domains exhibit a fold and heme arrangement very similar to that of cytochrome c(3), they differ considerably in their loop extensions and electrostatic surface. Analysis of the structure provides evidence for a different function of both domains, namely, anchoring the protein in a transmembranous complex with the N-terminal domain and formation of an electron-transfer complex with hydrogenase by the C-terminal domain.

The effect of mutation at valine-45 on the stability and redox potentials of trypsin-cleaved cytochrome b5

In an attempt to elucidate the determinants of redox potential and protein stability in cytochrome b5, three mutants at a highly conserved residue Val45, which is a member of heme hydrophobic pocket residues have been characterized. The V45Y mutant was designed to introduce a bulkier residue and a hydroxyl group to the heme pocket. The mutants V45H and V45E were constructed to test the effect of positive and negative charge on the stability and redox potential of proteins. The influence of these mutants on the protein stability towards thermal, urea, acid, ethanol and on the redox potential were studied. It is concluded that the decrease of hydrophobic free energy and the larger volume of the tyrosine make the phenylhydroxyl group of tyrosine still sitting inside the hydrophobic pocket, while the side chain of the mutant V45E and V45H shift away from the heme pocket. The redox potentials of mutants V45Y, V45H, V45E and wild-type of cytochrome b5 are -35 mV, 8 mV, -26 mV and -3 mV, respectively. The bigger change of the V45Y on redox potential is due to the close contact between the hydroxyl group and the heme, while the changes of the V45E and V45H result from the alteration of charge density and distribution around the heme. Different relative stability of these mutants towards heat have been observed with the order: WT > V45Y-V45H > V45E being both in the oxidized and reduced state. The relative stability induced by addition of urea decreases in the order: WT > V45Y > V45H > V45E. These results suggest that the difference in the hydrophobic free energy is a major factor contributing to the stability of the Val45 mutants. Also the loose of the helix III in the mutant V45E makes it more unstable. These results indicate that residue Val45 plays an important role in the stability and redox potential of the protein.

Effect of mutation at valine 61 on the three-dimensional structure, stability, and redox potential of cytochrome b5

To elucidate the role played by Val61 of cytochrome b(5), this residue of the tryptic fragment of bovine liver cytochrome b(5) was chosen for replacement with tyrosine (Val61Tyr), histidine (Val61His), glutamic acid (Val61Glu), and lysine (Val61Lys) by means of site-directed mutagenesis. The mutants Val61Tyr, Val61Glu, Val61His, and Val61Lys exhibit electronic spectra identical to that of the wild type, suggesting that mutation at Val61 did not affect the overall protein structure significantly. The redox potentials determined by differential pulse voltammetry were -10 (wild type), -25 (Val61Glu), -33 (Val61Tyr), 12 (Val61His), and 17 mV (Val61Lys) versus NHE. The thermal stabilities and urea-mediated denaturation of wild-type cytochrome b(5) and its mutants were in the following order: wild type > Val61Glu > Val61Tyr > Val61His > Val61Lys. The kinetics of denaturation of cytochrome b(5) by urea was also analyzed. The first-order rate constants of heme transfer between cytochrome b(5) and apomyoglobin at 20 +/- 0.2 degrees C were 0.25 +/- 0.01 (wild type), 0.42 +/- 0.02 (Val61Tyr), 0.93 +/- 0.04 (Val61Glu), 2.88 +/- 0.01 (Val61His), and 3.88 +/- 0.02 h(-)(1) (Val61Lys). The crystal structure of Val61His was determined using the molecular replacement method and refined at 2.1 A resolution, showing that the imidazole side chain of His61 points away from the heme-binding pocket and extends into the solvent, the coordination distances from Fe to NE2 atoms of two axial ligands are approximately 0.6 A longer than the reported value, and the hydrogen bond network involving Val61, the heme propionates, and three water molecules no longer exists. We conclude that the conserved residue Val61 is located at one of the key positions, the "electrostatic potential" around the heme-exposed area and the hydrophobicity of the heme pocket are determinant factors modulating the redox potential of cytochrome b(5), and the hydrogen bond network around the exposed heme edge is also an important factor affecting the heme stability.

Methanol-induced unfolding and refolding of cytochrome b5 and its P40V mutant monitored by UV-visible, CD, and fluorescence spectra

In order to illustrate the structural importance of proline-40 of cytochrome b5 (Cyt b5), the P40V mutant gene was constructed. Unfolding and refolding of Cyt b5 induced by methanol was investigated by means of the UV-visible spectrum, circular dichroism, and the fluorescence spectrum. Methanol denaturation of Cyt b5 is a cooperative process, that is, the heme group dissociates from the heme pocket accompanied by unfolding of the polypeptide chain both in the secondary and tertiary structures. Substitution of proline by valine reduces the stability of the mutant under methanol denaturation. The unfolding process is almost reversible by dilution. During refolding, the denatured polypeptide must be folded to a more ordered structure prior to the heme capture. Pro40 plays an important role in modulating the protein's stability. The role of tyrosine in the unfolding and refolding of Cyt b5 is evaluated for the first time. A mechanism of methanol denaturation is also proposed.

Engineering out motion: introduction of a de novo disulfide bond and a salt bridge designed to close a dynamic cleft on the surface of cytochrome b5

A previous molecular dynamics (MD) simulation of cytochrome b5 (cyt b5) at 25 degrees C displayed localized dynamics on the surface of the protein giving rise to the periodic formation of a cleft that provides access to the heme through a protected hydrophobic channel [Storch and Daggett (1995) Biochemistry 34, 9682]. Here we describe the production and testing of mutants designed to prevent the cleft from opening using a combination of experimental and theoretical techniques. Two mutants have been designed to close the surface cleft: S18D to introduce a salt bridge and S18C:R47C to incorporate a disulfide bond. The putative cleft forms between two separate cores of the protein: one is structural in nature and can be monitored through the fluorescence of Trp 22, and the other binds the heme prosthetic group and can be tracked via heme absorbance. An increase in motion localized to the cleft region was observed for each protein, except for the disulfide-containing variant, in MD simulations at 50 degrees C compared to simulations at 25 degrees C. For the disulfide-containing variant, the cleft remained closed. Both urea and temperature denaturation curves were nearly identical for wild-type and mutant proteins when heme absorbance was monitored. In contrast, fluorescence studies revealed oxidized S18C:R47C to be considerably more stable based on the midpoints of the denaturation transitions, Tm and U1/2. Moreover, the fluorescence changes for each protein were complete at approximately 50 degrees C and a urea concentration of approximately 3.9 M, significantly below the temperature and urea concentration (62 degrees C, 5 M urea) required to observe heme release. In addition, solvent accessibility based on acrylamide quenching of Trp 22 was lower in the S18C:R47C mutant, particularly at 50 degrees C, before heme release [presented in the accompanying paper (58)]. The results suggest that a constraining disulfide bond can be designed to inhibit dynamic cleft formation on the surface of cyt b5. Located near the heme, the native dynamics of the cleft may be functionally important for protein-protein recognition and/or complex stabilization.

The influence of mutation at Glu44 and Glu56 of cytochrome b5 on the protein's stabilization and interaction between cytochrome c and cytochrome b5

To characterize the roles played by Glu44 and Glu56 of cytochrome b5 in the formation of the electrostatic complex between cytochrome c and cytochrome b5, the Glu44, Glu56, or both sites were changed to alanine by site-directed mutagenesis. The influence of these two residues on the protein stability was probed by investigating the kinetic behaviors of protein denaturation in urea or upon heating and the heme-transfer reactions between apo-myoglobin and the variants of cytochrome b5. It has been found that when the Glu44 and/or Glu56 are mutated to alanine, the protein stability increases slightly due to the fact that the hydrophilic residue is changed to a hydrophobic one, resulting in the two pairs of helices surrounding the heme taking a more compact conformation. The difference in voltammetric behavior of cytochrome c, cytochrome b5, and its three mutants, Cyt b5 E44A, E56A, and E44/56A, alone and in 1:1 protein complexes demonstrates that both Glu44 and Glu56 of cytochrome b5 take part in the electrostatic interaction with cytochrome c. The entropy changes, DeltaS degreesrc and enthalpy changes, DeltaH degrees, derived from the temperature dependence of the formal reduction potentials of each protein in different protein systems suggest that, because of the mutual interaction with cytochrome c, cytochrome b5 mutants, especially the E44A-containing mutants, in the protein complexes suffer greater conformational changes upon reduction than that of the wild type. The variation of these thermodynamic parameters indicates that the strength of mutual interactions between cytochrome c and cytochrome b5 or its mutants has the following order: Cyt c/Cyt b5 > Cyt c/Cyt b5 E56A > Cyt c/Cyt b5 E44A > Cyt c/Cyt b5 E44/56A.

Probing the backbone dynamics of oxidized and reduced rat microsomal cytochrome b5 via 15N rotating frame NMR relaxation measurements: biological implications

Rotating frame 15N relaxation NMR experiments have been performed to study the local mobility of the oxidized and reduced forms of rat microsomal cytochrome b5, in the microsecond to millisecond time range. Measurements of rotating frame relaxation rates (R1rho) were performed as a function of the effective magnetic field amplitude by using off-resonance radio frequency irradiation. Detailed analysis of the two data sets resulted in the identification of slow motions along the backbone nitrogens for both oxidation states of the protein. The local mobility of reduced and oxidized cytochrome b5 turned out to be significantly different; 28 backbone nitrogens of the oxidized form were shown to participate in a conformational exchange process, while this number dropped to 12 in the reduced form. The correlation time, tauex, for the exchange processes could be determined for 21 and 9 backbone nitrogens for oxidized and reduced cytochrome b5, respectively, with their values ranging between 70 and 280 microseconds. The direct experimental evidence provided in this study for the larger mobility of the oxidized form of the protein is consistent with the different backbone NH solvent exchangeability recently documented for the two oxidation states [Arnesano, F., et al. (1998) Biochemistry 37, 173-184]. Our experimental observations may have significant biological implications. The differential local mobility between the two oxidation states is proposed to be an important factor controlling the molecular recognition processes in which cytochrome b5 is involved.

The origin of differences in the physical properties of the equilibrium forms of cytochrome b5 revealed through high-resolution NMR structures and backbone dynamic analyses

On the basis of a comparison of high-resolution solution structures calculated for both equilibrium forms of rat ferrocytochrome b5, differences in reduction potential and thermodyanmic stability have been characterized in terms of significant structural and dynamic differences between the two forms. The dominant difference between A and B conformations has long been known to be due to a 180 degrees rotation of the heme in the binding pocket about an axis defined by the alpha- and gamma-meso carbons, however, the B form has not been structurally characterized until now. The most significant differences observed between the two forms were the presence of a hydrogen bond between the 7-propionate and the S64 amide in the A form but not the B form and surprisingly a displacement of the heme out of the binding pocket by 0.9 A in the B form relative to the A form. The magnitude of other factors which could contribute to the known difference in reduction potentials in the bovine protein [Walker, F. A., Emrick, D., Rivera, J. E., Hanquet, B. J., and Buttlaire, D. H. (1988) J. Am. Chem. Soc. 110, 6234-6240], such as differences in the orientation of the axial imidazoles and differences in hydrogen bond strength to the imidazoles, have been evaluated. The dominant effector of the reduction potential would appear to be the lack of the hydrogen bond to the S64 amide in the B form which frees up the propionate to charge stabilize the iron in the oxidized state and thus lower the reduction potential of the B form. The structure we report for the A form, based on heteronuclear NMR restraints, involving a total of 1288 restraints strongly resembles both the X-ray crystal structure of the bovine protein and a recently reported structure for the A form of the rat protein based on homonuclear data alone [Banci, L., Bertini, I., Ferroni, F., and Rosato, A. (1997) Eur. J. Biochem. 249, 270-279]. The rmsd for the backbone atoms of the A form is 0.54 A (0.92 A for all non-hydrogens). The rmsd for the backbone of the B form is 0.51 A (0. 90 A for all non-hydrogen atoms). An analysis of backbone dynamics based on a model-free analysis of 15N relaxation data, which incorporated axially symmetric diffusion tensor modeling of the cytochrome, indicates that the protein is more rigid in the reduced state relative to the oxidized state, based on a comparison with order parameters reported for the bovine protein in the oxidized state [Kelly, G. P., Muskett, F. W., and Whitford, D. (1997) Eur. J. Biochem. 245, 349-354].

Role of carboxyl residues surrounding heme of human cytochrome b5 in the electrostatic interaction with NADH-cytochrome b5 reductase

To identify the cytochrome b5 residues responsible for the electrostatic interaction with NADH-cytochrome b5 reductase (b5R), we prepared and characterized the cytochrome b5 mutants in which Glu41, Glu42, Glu63, Asp70, and Glu73 were replaced by Ala, utilizing site-directed mutagenesis and the expression system for cytochrome b5 in Escherichia coli. Apparent Km values of the wild type b5R for Glu42Ala cytochrome b5 and Asp70Ala cytochrome b5 were approximately three-fold and six-fold higher than that for the wild type cytochrome b5, respectively, while the kcat values for those mutants were not remarkably affected. In contrast, Glu41Ala, Glu63Ala, and Glu73Ala cytochrome b5 showed almost the same kinetic properties as the wild type cytochrome b5. Furthermore, kinetic studies on combinations of the cytochrome b5 and b5R mutants suggested the interaction between Glu42 and Asp70 of cytochrome b5 and Lys125 and Lys41 of b5R, respectively, in the reaction.

Effect of axial ligand plane reorientation on electronic and electrochemical properties observed in the A67V mutant of rat cytochrome b5

Mutational studies directed at evaluating the effect of the axial ligand plane orientation on electrochemical properties of cytochrome b5 have been performed. As described in the previous paper, structural consequences of one of these mutations, the A67V mutation, have been evaluated using NMR solution methods. The lack of large shifts relative to the wild-type protein in both the imidazole Ndelta nitrogen and proton resonances of the H63 imidazole ring indicates that the hydrogen bond between the carbonyl of F58 and the imidazole ring of H63 remains intact in this mutant. Effects of the imidazole plane reorientation on the Fe d-orbitals were evaluated on the basis of interpretation of EPR spectra, near-infrared bands associated with ligand-to-metal charge transfer transitions, reorientation of the anisotropy of the paramagnetic center determined by calculation of pseudocontact shifts, and the temperature dependence of the contact-shifted resonances. The dominant effect of the imidazole reorientation appears to have been a destabilization of the d(xz) orbital energy and a reorientation of the d(pi) orbitals. This is surprising in light of the -20 mV shift in the reduction potential of the mutant relative to the wild-type protein and indicates that a destabilization of d(yz)-orbital energy level of the reduced state dictates the observed change in reduction potential. Measured values for the reorganizational energy and heterogeneous electron transfer rates were indistinguishable for wild-type and mutant proteins. This is perhaps surprising, given significant differences in the pattern of electron delocalization into the porphyrin ring observed as significantly altered contact shift patterns. Mutational studies perturbing the H39 imidazole were also performed but with more limited success.

Molecular dynamics simulation of cytochrome b5: implications for protein-protein recognition

Cytochrome b5 participates in electron-transfer reactions with a variety of different proteins. To explore how this protein might discern between structurally varied proteins, we have performed a molecular dynamics simulation focusing on its structural stability and dynamic behavior in solution. The protein was simulated in water at 298 K and pH 6.9 for 2.5 ns. The protein deviated significantly from the crystal structure midway through the simulation, but ultimately the crystalline conformation was regained. The simulation was at all times well behaved as judged by comparison to structural NMR data obtained in solution. One region of the protein backbone that deviated from the crystal conformation contains acidic residues implicated in electrostatic-based protein-protein recognition. The mobility in this region caused the protein to display different patterns of residues at the surface with time, as well as the formation of a large cleft partially exposing the hydrophobic core lining the heme pocket. Furthermore, the position and cyclical formation of this cleft suggest that hydrophobic interactions may be important in protein-protein recognition events and possibly even electron transfer, as the cleft allows for easy access to the heme group. These results indicate that thermal motion could provide a low-energy mechanism for controlling recognition events. Thus, the dynamical behavior observed through the varying solution conformations sampled may be important in influencing the diverse range of protein-protein interactions in which cytochrome b5 participates.

Experimental and theoretical analysis of the interaction between cytochrome c and cytochrome b5

Experimental and theoretical investigation of the interaction of cytochrome c and cytochrome b5 performed over nearly twenty years has produced considerable insight into the manner in which these proteins recognize and bind to each other. The results of these studies and the experimental and theoretical strategies that have been developed to achieve these results have significant implications for understanding the behavior of similar complexes formed by more complex and less-well characterized electron transfer proteins. The current review provides a comprehensive summary and critical evaluation of the literature on which the current status of our understanding of the interaction of cytochrome c and cytochrome b5 is based. The general issues related to the study of electron transfer complexes of this type are discussed and some new directions for future investigation of such systems are considered.

Cytochrome b5, its functions, structure and membrane topology

The first part of the present communication reviews recent advances in our understanding of the known physiological functions of cytochrome b5. In addition, one section is devoted to a description of a recently discovered function of cytochrome b5, namely its involvement in the synthesis of the oncofetal antigen N-glycolylneuraminic acid. The second part of the article summarizes site-directed mutagenesis studies, primarily conducted in the author's laboratory, in both the catalytic heme-binding and membrane-binding domain of cytochrome b5. These studies have shown that: 1) the membrane binding domain of cytochrome b5 spans the bilayer; 2) cytochrome b5 lacking 19 COOH-terminal amino acids does not bind to membrane bilayers; and 3) specific amino acids in the membrane binding domain have been mutated and shown not to be essential for the function of cytochrome b5 with its redox partners.

Tertiary structure of the heme-binding domain of rat cytochrome b5 based on homology modeling

The in vitro complexes formed between cytochrome b5 and other proteins (e.g. cytochrome c) have served as a useful means to probe electrostatic contributions to macromolecular recognition. Extensive experimentation has been carried out to test the specificity and stability of these complexes, including site-directed mutagenesis based on the heterologous expression of rat cytochrome b5 in E. coli. Despite this interest, there has not been a determination of the complete structure of cytochrome b5. Here we report coordinates for the complete tertiary structure of the heme-binding domain of rat cytochrome b5 based on homology modeling. Protein Data Bank (PDB) coordinates derived from the crystal structure of the highly homologous bovine cytochrome b5 were used for main chain scaffolding. Secondary structures for the termini missing in the bovine structure were generated using homologous sequences derived from an exhaustive search of the PDB database. The model structure was solvated and further refined using energy minimization techniques. The N-terminal residues of the model appear to be in a beta sheet conformation while the carboxy terminus is in a helical conformation. The rest of the rat model is folded virtually identically to the bovine x-ray crystal structure (r.m.s. deviation 1.28 A), despite six sequence differences between the two cores. This homology-based structure should be useful for structure-function analyses of molecular recognition involving cytochrome b5.

The cytochrome b5-fold: an adaptable module

The family of b5-like cytochromes encompasses, besides cytochrome b5 itself, hemoprotein domains covalently associated with other redox proteins, in flavocytochrome b2 (L-lactate dehydrogenase), sulfite oxidase and assimilatory nitrate reductase. A comparison of about 40 amino acid sequences deposited in data banks shows that eight residues are invariant and about 15 positions carry strongly conservative substitutions. Examination of the location of these invariant and conserved positions in the light of the three-dimensional structures of beef cytochrome b5 and S cerevisiae flavocytochrome b2 suggests a strongly conserved protein structure for the b5-like heme-binding domain throughout evolution. Numerous NMR studies have demonstrated the existence of a positional isomerism for the heme, which involves both a 180 degree-rotation around the heme alpha,gamma-meso carbon atoms and a rotation through an axis normal to the heme plane at the iron. NMR studies did not detect significant differences in protein structure between reduced and oxidized states, or between species. The role of a number of side chains was probed by site-directed mutagenesis. Studies of complex formation and of electron transfer rates between cytochrome b5 and redox partners have led to the idea that complexation is driven by electrostatic forces, that it is generally the exposed heme edge which makes contact with electron donors and acceptors, but that there are multiple overlapping sites within this general area. For the bi- and trifunctional members of the family, extrapolation of available data would suggest a mobile heme-binding domain within a complex structure. In these cases the existence of a single interaction area for both electron donor and acceptor, or of two different ones, remains open to discussion.

A hypothetical complex between crystalline flavocytochrome b2 and cytochrome c

Flavocytochrome b2 and cytochrome c are physiological electron transfer partners in yeast mitochondria. The formation of a stable complex between them has been demonstrated both in solution and in the crystalline state. On the basis of the three-dimensional structures, using molecular modeling and energy minimization, we have generated a hypothetical model for the interaction of these redox partners in the crystal lattice. General criteria such as good charge and surface complementarity, plausible orientation, and separation distance of the prosthetic groups, as well as more specific criteria such as the stoichiometry determined in the crystal, and the involvement of both domains and of more than one subunit of flavocytochrome b2 led us to discriminate between several possible interaction sites. In the hypothetical model we present, four cytochrome c molecules interact with a tetramer of flavocytochrome b2. The b2 and c hemes are coplanar, with an edge-to-edge distance of 14 A. The contact surface area is ca. 800 A2. Several electrostatic interactions involving the flavin and the heme domains of flavocytochrome b2 stabilize the binding of cytochrome c.

Cytochrome C (Fe2+) as a competitive inhibitor of NADPH-dependent reduction of cytochrome P450 LM2: locating protein-protein interaction sites in microsomal electron carriers

The kinetics of NADPH-dependent reduction of cytochrome P450 LM2 in the soluble monomeric reconstituted system in the absence of any substrate is shown to be monophasic. We show that ferrous cytochrome c acts as a competitive inhibitor of the reduction. In the presence of 1 mM benzphetamine an additional extremely fast phase was observed. Under these conditions ferrous cytochrome c was found to be a competitive inhibitor of the slow phase of the reduction process, which accounted for 80% of the total reduction amplitude. Inhibition experiments yield a dissociation constant for the LM2-reductase complex of 3.0 +/- 1.5 microM. This constant was the same both in the presence and in the absence of benzphetamine. Based on these data we conclude that cytochromes P450 and c bind to the same center on the NADPH-cytochrome P450 reductase molecule. Comparative analysis of the amino acid sequences reveals a detectable similarity between cytochrome c and cytochrome P450 LM2 at positions 68-87 and 121-145, respectively. In addition, a substantial similarity was shown for sequence fragments 204-224 of NADPH-cytochrome P450 reductase and 40-60 of cytochrome b5. Based on these findings a hypothesis for the location of the centers of intermolecular interactions on the molecules of cytochrome P450 LM2 and NADPH-cytochrome P450 reductase is proposed.

Electron-transfer self-exchange kinetics of trimethylphosphine horse-heart myoglobin

Electron self-exchange has been measured by an NMR technique for horse-heart myoglobin. The rate is 3.1 x 10(3) M-1 s-1 at 23 degrees in 0.1 M phosphate at pH 6.9. The rate was weakly dependent on ionic strength up to 0.7 M in added KCl (3.9 x 10(3) M-1 s-1). The enthalpy of activation was 12.1 +/- 0.5 kcal mol-1, and the entropy of activation was -1.2 +/- 0.5 cal mol-1 deg-1. Analysis of the data in terms of the Marcus theory gives a reorganization energy, lambda, for self-exchange of 1.6 eV mol-1.

Mapping electrostatic interactions in macromolecular associations

In the association of electron transfer proteins, electrostatics has been proposed to play a role in maintaining the stability and specificity of the biomolecular complexes formed. An excellent model system is the interaction between mammalian cytochrome b5 and cytochrome c, in which the X-ray structures of the individual components reveal a complementary asymmetry of charges surrounding their respective redox centers. Determining the exact extent of the electrostatic interactions and identifying the specific residues involved in the formation of the electron transfer complex has proved more elusive. We report herein the utilization of high-pressure techniques, together with site-directed mutagenesis, to provide a map of the interaction domains in biomolecular complex formation. The application of high pressure disrupts macromolecular associations since dissociation of the complex results in a decreased volume of the system due to the solvation of charges that had been previously sequestered in the interface region and force solvation of hydrophobic surfaces. Site-directed mutagenesis of a totally synthetic gene for rat liver cytochrome b5, which expresses this mammalian protein in Escherichia coli as a hemecontaining soluble component, was used to selectively alter negatively charged residues of cytochrome b5 to neutral amide side-chains. We have demonstrated that the interaction domain of cytochrome b5 with cytochrome c can be mapped from a comparison of dissociation volumes of these modified cytochrome b5-cytochrome c complexes with the native complex. Using these techniques we can specifically investigate the role of particular residues in the equilibrium association of these two electron transfer proteins. Single-point mutations in the interaction domain give nearly identical effects on the measured dissociation volumes, yet removal of acidic residues outside the recognition surface yield volumes similar to wild-type protein. Multiple mutations in the proposed protein-protein interaction site are found to allow greater solvent-accessibility of the interface as reflected in a diminution in the volume changes on subsequent charge removal. This is indicative that the interprotein salt-bridges in this complex provide a mechanism for a greater exclusion of solvent from the interfacial domain of the complex, resulting in a more stable association.

Characterisation of the electron self-exchange rates in hexametaphosphate-cytochrome-c aggregates measured using high-resolution 1H-NMR spectroscopy

1H-NMR spectroscopy has been used to measure the rate of unimolecular electron exchange between cytochrome c molecules in protein aggregates stabilised by the addition of sodium hexametaphosphate. The average intracomplex electron exchange rate is measured from line broadening of hyperfine-shifted resonances of ferricytochrome c in an equimolar mixture of reduced and oxidised protein. The line-broadening due to electron exchange is significantly greater than that due to protein aggregation and reaches a maximum value between 1-2 mol hexametaphosphate/mol protein. Significantly the exchange-induced broadening is a first-order process and is directly proportional to the size of the cytochrome c oligomer. From the temperature dependence of exchange broadening the activation enthalpy was estimated to be 75.8 kJ mol-1 whereas the activation entropy was 295 J mol-1 K-1 for a dimer of cytochrome c at a hexametaphosphate/protein molar ratio of 1. Both activation parameters decrease in magnitude as the order of the cytochrome c oligomer increases. The rates of intracomplex electron exchange in Saccharomyces cerevisiae iso-2 and Candida krusei cytochromes c are lower than that of the horse protein, implying that primary sequence plays a fundamental part in determining the rate of exchange. The relevance of these observations is discussed in terms of the function of cytochrome c.

Kinetic properties of purified sheep lung microsomal NADH-cytochrome b5 reductase

1. Lung NADH-cytochrome b5 reductase was saturated with its artificial substrate, potassium ferricyanide at approximately 0.1 mM ferricyanide concentration, and the activity of the lung enzyme was inhibited by the higher concentrations of potassium ferricyanide. Ferricyanide at 0.5 and 1.0 mM inhibited the activity of the enzyme by about 20 and 61% respectively. The apparent Km value was calculated as 13.7 microM potassium ferricyanide and 4.3 microM NADH. 2. The Michaelis constants for cytochrome b5 and NADH were determined to be 1.67 and 7.7 microM from the Lineweaver-Burk plots. These results demonstrate that affinity of the lung reductase for its natural substrate is almost 10 times higher than that for potassium ferricyanide. 3. Addition of non-ionic detergent stimulated the rate of reductase-catalyzed reduction of lung cytochrome b5 up to 8.2-fold. 4. Kinetic studies performed with lung reductase by varying NADH and cytochrome b5 concentrations at different fixed concentrations at cytochrome b5 or NADH showed a series of parallel lines indicating a "ping-pong" type of kinetic mechanism for interaction of NADH and cytochrome b5 with lung cytochrome b5 reductase.

Molecular structure of flavocytochrome b2 at 2.4 A resolution

The crystal structure of flavocytochrome b2 has been solved at 3.0 A resolution by the method of multiple isomorphous replacement with anomalous scattering. Area detector data from native and two heavy-atom derivative crystals were used. The phases were refined by the B.C. Wang phase-filtering procedure utilizing the 67% (v/v) solvent content of the crystals. A molecular model was built first on a minimap and then on computer graphics from a combination of maps both averaged and not averaged about the molecular symmetry axis. The structure was extended to 2.4 A resolution using film data recorded at a synchrotron and refined by the Hendrickson-Konnert procedure. The molecule, a tetramer of Mr 230,000, is located on a crystallographic 2-fold axis and possesses local 4-fold symmetry. Each subunit is composed of two domains, one binding a heme and the other an FMN prosthetic group. In subunit 1, both the cystochrome and the flavin-binding domain are visible in the electron density map. In subunit 2 the cytochrome domain is disordered. However, in the latter, a molecule of pyruvate, the product of the enzymatic reaction, is bound at the active site. The cytochrome domain consists of residues 1 to 99 and is folded in a fashion similar to the homologous soluble fragment of cytochrome b5. The flavin binding domain contains a parallel beta 8 alpha 8 barrel structure and is composed of residues 100 to 486. The remaining 25 residues form a tail that wraps around the molecular 4-fold axis and is in contact with each remaining subunit. The FMN moiety, which is located at the C-terminal end of the central beta-barrel, is mostly sequestered from solvent; it forms hydrogen bond interactions with main- and side-chain atoms from six of the eight beta-strands. The interaction of Lys349 with atoms N-1 and O-2 of the flavin ring is probably responsible for stabilization of the anionic form of the flavin semiquinone and hydroquinone and enhancing the reactivity of atom N-5 toward sulfite. The binding of pyruvate at the active site in subunit 2 is stabilized by interaction of its carboxylate group with the side-chain atoms of Arg376 and Tyr143. Residues His373 and Tyr254 interact with the keto-oxygen atom and are involved in catalysis. In contrast, four water molecules occupy the substrate-binding site in subunit 1 and Tyr143 forms a hydrogen bond to the ordered heme propionate group. Otherwise the two flavin-binding domains are identical within experimental error.(ABSTRACT TRUNCATED AT 400 WORDS)

One- and two-dimensional nuclear Overhauser effect studies of the electronic/molecular structure of the heme cavity of ferricytochrome b5

A nuclear Overhauser effect, NOE, study of solubilized native bovine ferricytochrome b5 has provided the complete assignment of the heme resonances as well as those of the majority of the amino acid side-chains making contact with the prosthetic group. The resonances which could not be identified are those from positions very close to the iron (less than 5 A) for which paramagnetic relaxation is sufficiently strong to significantly decrease the NOEs. The observed 1H-1H dipolar contacts generally confirm a solution structure unchanged from that described in single crystals, except for the detailed orientation of the heme side-chains. The 2-vinyl group is found in both the cis and trans in-plane orientation as opposed to exclusively cis in the crystal, and the 7-propionate group is rotated by 30 degrees in solution towards the 6-propionate group. Identification of resonances for the individual axial histidine residues indicates non-equivalent interaction with the heme iron, and the patterns of meso-H, pyrrole substituent and amino acid dipolar shifts allow the location of the principal magnetic axes in the protein coordinate system. This identifies His-39 as the dominant influence in determining the electronic ground state that orients the molecular orbital for facile electron transfer via the exposed heme edge. The complete two-dimensional NOESY map for ferricytochrome b5 is presented that yields all the cross peaks expected on the basis of the one-dimensional NOE studies, and indicates that such two-dimensional methods should have profitable extension to strongly hyperfine-shifted resonances in paramagnetic proteins.