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

Structure analysis and characterization of the cytochrome c-554 from thermophilic green sulfur photosynthetic bacterium Chlorobaculum tepidum

The cytochrome (Cyt) c-554 in thermophilic green photosynthetic bacterium Chlorobaculum tepidum serves as an intermediate electron carrier, transferring electrons to the membrane-bound Cyt c z from various enzymes involved in the oxidations of sulfide, thiosulfate, and sulfite compounds. Spectroscopically, this protein exhibits an asymmetric α-absorption band for the reduced form and particularly large paramagnetic (1)H NMR shifts for the heme methyl groups with an unusual shift pattern in the oxidized form. The crystal structure of the Cyt c-554 has been determined at high resolution. The overall fold consists of four α-helices and is characterized by a remarkably long and flexible loop between the α3 and α4 helices. The axial ligand methionine has S-chirality at the sulfur atom with its C(ε)H3 group pointing toward the heme pyrrole ring I. This configuration corresponds to an orientation of the lone-pair orbital of the sulfur atom directed at the pyrrole ring II and explains the lowest-field (1)H NMR shift arising from the 18(1) heme methyl protons. Differing from most other class I Cyts c, no hydrogen bond was formed between the methionine sulfur atom and polypeptide chain. Lack of this hydrogen bond may account for the observed large paramagnetic (1)H NMR shifts of the heme methyl protons. The surface-exposed heme pyrrole ring II edge is in a relatively hydrophobic environment surrounded by several electronically neutral residues. This portion is considered as an electron transfer gateway. The structure of the Cyt c-554 is compared with those of other Cyts c, and possible interactions of this protein with its electron transport partners are discussed.

Characterization of Heme-Coordinating Histidyl Residues of Cytochrome b5 Based on the Reactivity with Diethylpyrocarbonate: A Mechanism for the Opening of Axial Imidazole Rings

We investigated the reactivity of heme-coordinating imidazole with diethylpyrocarbonate using a soluble domain of cytochrome b(5). Analyses with various spectroscopic methods including MALDI-TOF-MS indicated that two axial His residues (His44 and His68) of cytochrome b(5) were protected from the modification by several factors, i.e., limited steric exposure of the axial imidazole to the solvent, the Fe-N(epsilon2) coordination bond, and protonation of the N(delta1) position by forming a hydrogen bond with its immediate surroundings. However, once N-carbethoxylation at the N(epsilon2) position of the axial His residues occurred with a higher concentration of diethylpyrocarbonate, displacement of heme prosthetic group from the protein moiety continued. Simultaneously, it facilitated the second N-carbethoxylation to take place at the N(epsilon1) position of the same imidazole ring, leading to a bis-N-carbethoxylated derivative and further to a ring-opened derivative. A similar mechanism seemed in operation for one non-axial His residue (His85), in which the N(delta1) atom works as a hydrogen acceptor in a strong hydrogen-bond and the other N(epsilon2) atom is in a protonated form, resulting in a formation of the ring-opened derivative upon treatment with a higher concentration of diethylpyrocarbonate. These results suggested that the use of diethylpyrocarbonate for MALDI-TOF-MS analysis might provide a unique method to characterize the protonation state of His residues and the strength of their hydrogen-bondings at the active site of enzymes.

Binding of Ferric Heme by the Recombinant Globin from the Cyanobacterium Synechocystis sp. PCC 6803

The product of the cyanobacterium Synechocystis sp. PCC 6803 gene slr2097 is a 123 amino acid polypeptide chain belonging to the truncated hemoglobin family. Recombinant, ferric heme-reconstituted Synechocystis sp. PCC 6803 hemoglobin is a low-spin complex whose endogenous hexacoordination gives rise to optical and NMR characteristics reminiscent of cytochrome b(5) [Scott, N. L., and Lecomte, J. T. J. (2000) Protein Sci. 9, 587-597]. In this work, the sequential assignments using (15)N-(13)C-labeled protein, (1)H nuclear Overhauser effects, and longitudinal relaxation data identified His70 as the proximal histidine and His46 as the sixth ligand to the iron ion. It was also found that one of two possible heme orientations within the protein matrix is highly preferred (>90%) and that this orientation is the same as in vertebrate myoglobins. The rate constant for the 180 degrees rotation of the heme within a protein cage to produce the favored isomer was 0.5 h(-1) at 25 degrees C, approximately 35 times faster than in sperm whale myoglobin. Variable temperature studies revealed an activation energy of 132 +/- 4 kJ mol(-1), similar to the value in metaquomyoglobin at the same pH. The rate constant for heme loss from the major isomer was estimated to be 0.01 h(-1) by optical spectroscopy, close to the value for myoglobin and decades slower than in the related Nostoc commune cyanoglobin. The slow heme loss was attributed in part to the additional coordination bond to His46, whereas the relatively fast rate of heme reorientation suggested that this bond was weaker than the proximal His70-Fe bond. The standard reduction potential of the hexacoordinated protein was measured with and without poly-L-lysine as a mediator and found to be approximately -150 mV vs SHE, indicating a stabilization of the ferric state compared to most hemoglobins and b(5) cytochromes.

Structural and dynamic perturbations induced by heme binding in cytochrome b5

The water-soluble domain of rat hepatic cytochrome b(5) undergoes marked structural changes upon heme removal. The solution structure of apocytochrome b(5) shows that the protein is partially folded in the absence of the heme group, exhibiting a stable module and a disordered heme-binding loop. The quality of the apoprotein structure in solution was improved with the use of heteronuclear NMR data. Backbone amide hydrogen exchange was studied to characterize cooperative units in the protein. It was found that this criterion distinguished the folded module from the heme-binding loop in the apoprotein, in contrast to the holoprotein. The osmolyte trimethylamine N-oxide (TMAO) did not affect the structure of the apoprotein in the disordered region. TMAO imparted a small stabilization consistent with an unfolded state effect correlating with the extent of buried surface area in the folded region of the native apoprotein. The failure of the osmolyte to cause large conformational shifts in the disordered loop supported the view that the specificity of the local sequence for the holoprotein fold was best developed with the stabilization of the native state through heme binding. To dissect the role of the heme prosthetic group in forcing the disordered region into the holoprotein conformation, the axial histidine belonging to the flexible loop (His63) was replaced with an alanine, and the structural properties of the protein with carbon-monoxide-ligated reduced iron were studied. The His63Ala substitution resulted in a protein with lower heme affinity but nevertheless capable of complete refolding. This indicated that the coordination bond was not necessary to establish the structural features of the holoprotein. In addition, the weak binding of the heme in this protein resulted in conformational shifts at a location distant from the binding site. The data suggested an uneven distribution of cooperative elements in the structure of the cytochrome.

Engineering out motion: a surface disulfide bond alters the mobility of tryptophan 22 in cytochrome b5 as probed by time-resolved fluorescence and 1H NMR experiments

In the accompanying paper [Storch et al. (1999) Biochemistry 38, 5054-5064] equilibrium denaturation studies and molecular dynamics (MD) simulations were used to investigate localized dynamics on the surface of cytochrome b5 (cyt b5) that result in the formation of a cleft. In those studies, an S18C:R47C disulfide mutant was engineered to inhibit cleft mobility. Temperature- and urea-induced denaturation studies revealed significant differences in Trp 22 fluorescence between the wild-type and mutant proteins. On the basis of the results, it was proposed that wild type populates a conformational ensemble that is unavailable to the disulfide mutant and is mediated by cleft mobility. As a result, the solvent accessibility of Trp 22 is decreased in S18C:R47C, suggesting that the local environment of this residue is less mobile due to the constraining effects of the disulfide on cleft dynamics. To further probe the structural effects on the local environment of Trp 22 caused by inhibition of cleft formation, we report here the results of steady-state and time-resolved fluorescence quenching, differential phase/modulation fluorescence anisotropy, and 1H NMR studies. In Trp fluorescence experiments, the Stern-Volmer quenching constant increases in wild type versus the oxidized disulfide mutant with increasing temperature. At 50 degrees C, KSV is nearly 1.5-fold greater in wild type compared to the oxidized disulfide mutant. In the reduced disulfide mutant, KSV was the same as wild type. The bimolecular collisional quenching constant, kq, for acrylamide quenching of Trp 22 increases 2.7-fold for wild type and only 1.8-fold for S18C:R47C, upon increasing the temperature from 25 to 50 degrees C. The time-resolved anisotropy decay at 25 degrees C was fit to a double-exponential decay for both the wild type and S18C:R47C. Both proteins exhibited a minor contribution from a low-amplitude fast decay, consistent with local motion of Trp 22. This component was more prevalent in the wild type, and the fractional contribution increased significantly upon raising the temperature. The fast rotational component of the S18C:R47C mutant was less sensitive to increasing temperature. A comparison of the 1H NMR monitored temperature titration of the delta-methyl protons of Ile 76 for wild type and oxidized disulfide mutant, S18C:R47C, showed a significantly smaller downfield shift for the mutant protein, suggesting that Trp 22 in the mutant protein experiences comparatively decreased cleft dynamics in core 2 at higher temperatures. Furthermore, comparison of the delta'-methyl protons of Leu 25 in the two proteins revealed a difference in the ratio of the equilibrium heme conformers of 1.2:1 for S18C:R47C versus 1.5:1 for wild type at 40 degrees C. The difference in equilibrium heme orientations between wild type and S18C:R47C suggests that the disulfide bond affects heme binding within core 1, possibly through damped cleft fluctuations. Taken together, the NMR and fluorescence studies support the proposal that an engineered disulfide bond inhibits the formation of a dynamic cleft on the surface of cyt b5.

Structure, interaction and electron transfer between cytochrome b5, its E44A and/or E56A mutants and cytochrome c

Site-directed mutagenesis has been used to produce variants of a tryptic fragment of bovine liver cytochrome b5 in which Glu44 and Glu56 are mutated to alanine. The reduction potentials measured by spectroelectrochemical titration (in the presence of 1 mM (Ru(NH3)6)3+, pH 7.0 and I=0.1 M) are 4.5, 6.0, 6.0 and 7.5 mV versus the standard hydrogen electrode (SHE) for the wild-type and E44A, E56A and E44/56A mutants of cytochrome b5, respectively. A comparative two-dimensional NMR study of cytochrome b5 and its E44/56A mutant in water solution has been achieved. Resonance assignments of side-chains have been completed successfully. The NMR results suggest that the secondary structures and global folding of the E44/56A mutant remain unchanged, but the mutation of both Glu44 and Glu56 to hydrophobic alanine may lead to the two helices containing mutated residues contracting towards the heme center. The inner mobility of the Gly42 approximately Glu44 segment in cytochrome b5 may be responsible for the difference of the binding mode between Glu44 and Glu56 with cytochrome c. The binding between cytochrome c and cytochrome b5 was studied by optical difference spectra of cytochrome c and variants of cytochrome b5. The association constants (KA) for the wild-type, E44A, E56A, and E44/56A mutants of cytochrome b5 with cytochrome c, are 4.70(+/-0. 10)x10(6) M-1, 1.88(+/-0.03)x10(6) M-1, 2.70(+/-0.13)x10(6) M-1, and 1.14(+/-0.05)x10(6) M-1, respectively. This is indicative that both Glu44 and Glu56 are involved in the complex formation between cytochrome b5 and cytochrome c. The reduction of horse heart ferricytochrome c by recombinant ferrocytochrome b5 and its mutants has been studied. The rate constant of the electron transfer reaction between ferricytochrome c and wild-type ferrocytochrome b5 (1.074(+/-0.49)x10(7) M-1 s-1) is higher than those of the mutant protein E44A (8.98(+/-0.20)x10(6) M-1 s-1), E56A (8.76(+/-0. 39)x10(6) M-1 s-1), and E44/56A (8.02(+/-0.38)x10(6) M-1 s-1) at 15 degreesC, pH 7.0, I=0.35 M. The rate constants are strongly dependent on ionic strength and temperature. These studies, by means of a series of techniques, provide conclusive results that the interaction between cytochrome b5 and cytochrome c is electrostatically guided, and, more importantly, that both Glu44 and Glu56 participate in the electron transfer 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.

Characterization of a site-directed mutant of cytochrome b5 designed to alter axial imidazole ligand plane orientation

Mutants of cytochrome b5 were designed to achieve reorientation of individual axial imidazole ligands. The orientation of the axial ligand planes is thought to modulate the reduction potential of bis(imidazole) axially ligated heme proteins. The A67V mutation achieved this goal through the substitution of a bulkier, hydrophobic ligand for a residue, in the sterically hindered hydrophobic heme binding pocket. Solution structures of mutant and wild-type proteins in the region of the mutation were calculated using restraints obtained from 1H and 15N 2D homonuclear and heteronuclear NMR spectra and 1H-15N 3D heteronuclear NMR spectra. More than 10 restraints per residue were used in the refinement of both structures. Average local rmsd for 20 refined structures was 0.30 A for the wild-type structure and 0.38 A for the A67V mutant. The transfer of amide proton resonance assignments from wild-type to the mutant protein was achieved through overlays of 15N-1H heteronuclear correlation spectra of the reduced proteins. Side chain assignments and sequential assignments were established using conventional assignment strategies. Calculation of the orientation of the components of the anisotropic paramagnetic susceptibility tensor, using methods similar to procedures applied to the wild-type protein, shows that the orientation of the in-plane components are identical in the wild-type and mutant proteins. However, the orientation of the z-component of the susceptibility tensor calculated for the mutant protein differs by 17 degrees for the A-form and by 11 degrees for the B-form from the orientation calculated for the wild-type protein. The rotation of the z-component of the susceptibility tensor (toward the delta meso proton) is in the same direction and is of the same magnitude as the rotation of the H63 imidazole ring induced by mutation.

Molecular cloning, overexpression in Escherichia coli, structural and functional characterization of house fly cytochrome b5

A microsomal cytochrome b5 cDNA from the house fly, Musca domestica, was cloned and sequenced. The deduced amino acid sequence of the full-length house fly cytochrome b5 (134 residues) is 48% identical to that of rat microsomal cytochrome b5. The house fly cytochrome b5 protein was overexpressed in Escherichia coli, purified, and characterized. Absorption and EPR spectroscopy reveal properties very similar to cytochromes b5 from vertebrates. NMR spectra indicate that the orientation of the heme in the protein relative to its alpha,gamma meso axis is about 1:1. A redox potential of -26 mV versus standard hydrogen electrode was measured by cyclic voltammetry on a modified gold electrode in the presence of hexamminechromium(III) chloride. The cytochrome b5 is reduced by house fly cytochrome P450 reductase in a reconstituted system at a high rate (5.5 s-1), and it stimulates heptachlor epoxidation when reconstituted with house fly cytochrome P450 reductase, cytochrome P450 6A1, phospholipid, and detergent. Cytochrome b5 decreases the apparent Km for P450 reductase and increases the Vmax for heptachlor epoxidation at constant cytochrome P450 6A1 concentrations. The results indicate that cytochrome b5 stimulates a step following the first electron transfer during cytochrome P450 6A1 turnover.

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.

Assignment of hyperfine-shifted heme carbon resonances in ferricytochrome b5

The reverse detection heteronuclear multiple quantum coherence, HMQC, study of native bovine ferricytochrome b5 has provided the complete assignment of hyperfine shifted resonances of heme carbons attached proton(s). The dominant delocalized pi-spin density to vinyl groups gives rise to contact shifts which have opposite direction for a carbon and its attached proton(s). The most hyperfine shifted 13C heme signals are mainly generated from 3rd heme pyrrole ring substituents which identifies that the molecular orbital for facile electron transfer is oriented to exposed heme edge. Magnetic/electronic asymmetry of heme induced by two axial His makes spread the hyperfine shifted heme carbon resonances over the range of 280 ppm at 25 degrees C, which would be the more sensitive probe than those of proton resonances in characterizing the nature of heme electronic structure of ferricytochrome b5.

Hydrogen isotope effects on the proton nuclear magnetic resonance spectrum of bovine ferricytochrome b5: axial hydrogen bonding involving the axial His-39 imidazole ligand

The potential role of hydrogen bonding interactions in modulating the molecular and electronic structure of the active site of solubilized bovine ferricytochrome b5 has been investigated by monitoring solvent isotope effects on proton-NMR spectral parameters. It is observed that the hyperfine shifts of both the heme prosthetic group and one coordinated His are sensitive, while those for the other axial His and non-coordinated residues are insensitive, to 2H for 1H exchange. Two types of isotope influences are characterized; one whose chemical shift influence is time-resolved on the NMR time scale, and involves a single proton on one axial ligand, and a second effect which involves multiple protons, is not time resolved, and influences primarily the heme. A large isotope effect on the hyperfine shift is identified for the C beta H signals of His-39 but not His-63. The exchangeable ring NH of His-39 is assigned, and the pH influence on the exchange properties of heme pocket labile protons, when compared to the rate of base catalyzed averaging of the His-39 C beta H isotope effect, lead to the conclusion that the axial hydrogen bond which is responsible for this isotope effect is that between His-39 ring NH and Gly-42 carbonyl. The more rapid exchange of labile protons with solvent for His-63 than His-39 confirms a less solvent accessible and stronger hydrogen bonded His-39 than His-63. The stronger His-39-Gly-42 than His-63-Phe-58 hydrogen bond involving the ring NH leads to more extensive His-39 imidazolate character and hence a stronger iron-His-39 than iron-His-63 bond. The much larger hyperfine shifts for His-39 than His-63 imidazole ring non-labile protons support the stronger bonding of the former ligand, and account for the orientation of the rhombic magnetic axes by His-39 rather than His-63. The solvent isotope effect on the heme leads to rotation of the prosthetic group about the His-Fe-His bond by approximately 0.5 degrees so as to shorten the 7-propionate link to Ser-64. This suggests that the hydrogen bonds between the 7-propionate group and Ser-64 are responsible for the effect.

Proton NMR study of the heme complex of hemopexin

Proton nuclear magnetic resonance spectroscopy of the complex of heme with hemopexin, a plasma protein with an exceptionally high affinity for heme, is reported. Characteristic spectra are shown for heme.hemopexin of cow, human, rabbit, and rat. Each of these spectra demonstrate that the iron of heme bound by hemopexin is paramagnetic and low-spin. Rabbit heme.hemopexin, which exhibits the best signal-to-noise ratio, is studied in detail. Deuterium isotope labeling experiments indicate that the methyls in heme positions 1-, 3-, and 8- are resolved downfield from the protein envelope of resonances; the 5-methyl may lie in the -5 to +12 ppm region. Two-dimensional nuclear Overhauser effect spectroscopy locates other protons of the heme periphery, including from the 2-vinyl. Strongly relaxed upfield resonances are identified and assigned to protons on the axial ligands. Cyanide interaction with heme.hemopexin produces an additional low-spin adduct.

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.

Interpretation of hyperfine shift patterns in ferricytochromes b5 in terms of angular position of the heme: a sensitive probe for peripheral heme protein interactions

The 1H-NMR hyperfine shift pattern of the heme in a variety of low-spin ferricytochromes b5 has been analyzed in terms of the angular position of the prosthetic group within a structurally and magnetically-conserved protein matrix. A simple model is presented in which the changes in the spread of the predominantly contact shifted methyl and predominantly dipolar shifted meso-H signals of the heme, as well as shift trends for individual signals, provide sensitive indicators of the orientation of the heme relative to the orbital hole (singly-occupied d orbital), which in turn is related to the rhombic magnetic axes. The invariance of the axial His and non-coordinated residue hyperfine shifts show that it is the heme within a relatively rigid protein matrix, rather than the magnetic coordinate system, which is displaced angularly about the heme normal in order to accommodate variations in the polypeptide, orientation of the heme about the alpha,gamma-meso axis, and the length of heme carboxylate chains. Native heme shows increased counterclockwise rotation about the heme normal in the order rat-->beef-->chicken ferricytochrome b5, which is attributed largely to increased bulk of a variable sequence hydrophobic cluster consisting of residues 23, 25 and 32. The two alternate heme orientations about the alpha,gamma-meso axis are shown to also differ by rotation about the heme normal. A semiquantitative estimate of the degree of angular accommodation based on the spread of the meso-H rhombic dipolar shifts indicate rotations of 2-10 degrees. Possible functional consequences of such angular accommodation in relation to the role of these proteins in electron transfer are discussed.

Characterization of the haem environment in Methylophilus methylotrophus ferricytochrome c" by 1H-NMR

Two-dimensional NMR techniques have been used to assign proton resonances in the haem cavity of Methylophilus methylotrophus cytochrome c", a monohaem protein with bis-histidinyl ligation which has been shown to couple electron and proton transfer. All the assignments were made directly for the oxidized paramagnetic form of the cytochrome. Nearly all of the haem protons (90%) and the protons of both axial ligands have been assigned; the side-chain protons from four other residues in the haem pocket have also been identified. The data indicate a highly symmetric unpaired-electron distribution in the haem group, which agrees with a perpendicular orientation of the axial imidazole planes. The two haem propionate groups have contrasting degrees of exposure to the solvent, with the propionate group at position 13 being highly exposed. To obtain information on the dynamics of the haem environment, measurements of the 1H/2H-exchange rates of amide protons located in the haem cavity were performed. The two faces of the haem are found to differ markedly with respect to water accessibility. All of this information, together with additional protein sequencing data, indicates that His52 remains attached upon reduction and that the redox-linked protonation occurs via a channel running through the haem cleft on the opposite face.

Gene synthesis, bacterial expression, and 1H NMR spectroscopic studies of the rat outer mitochondrial membrane cytochrome b5

The gene coding for the water-soluble domain of the outer mitochondrial membrane cytochrome b5 (OM cytochrome b5) from rat liver has been synthetized and expressed in Escherichia coli. The DNA sequence was obtained by back-translating the known amino acid sequence [Lederer, F., Ghrir, R., Guiard, B., Cortial, S., & Ito, A. (1983) Eur. J. Biochem. 132, 95-102]. The recombinant OM cytochrome b5 was characterized by UV-visible, EPR, and 1H NMR spectroscopy. The UV-visible and EPR spectra of the OM cytochrome b5 are almost identical to the ones obtained from the overexpressed rat microsomal cytochrome b5 [Bodman, S. B. V., Schyler, M. A., Jollie, D. R., & Sligar, S. G. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 9443-9447]. The one-dimensional 1H NMR spectrum of the OM cytochrome b5 indicates that the rhombic perturbation of the ferric center is essentially identical to that in the microsomal beef, rabbit, chicken, and rat cytochromes b5. Two-dimensional 1H NMR spectroscopy (NOESY) and one-dimensional NOE difference spectroscopy were used to assign the contact-shifted resonances that correspond to each of the two isomers that result from the rotation of the heme around its alpha-gamma-meso axis. The assignment of the resonances allowed the determination of the heme orientation ratio in the OM cytochrome b5, which was found to be 1.0 +/- 0.1. It is noteworthy that the two cytochromes b5 that have similar populations of the two heme isomers (large heme disorder) originate from the rat liver.

Sequential resonance assignments of oxidized high-potential iron-sulfur protein from Chromatium vinosum

2D NMR spectra of the high-potential iron-sulfur protein (HiPIP) from Chromatium vinosum have been used to obtain partial resonance assignments for the oxidized paramagnetic redox state of the protein. Sequence-specific assignments were made using NOESY and COSY spectra in H2O and D2O of the following backbone segments: Asn-5-Arg-33, Glu-39-Asp-45, Gly-55-Cys-63, Gly-68-Ala-78, and Leu-82-Gly-85. NOESY spectra with a spectral width wide enough to include the hyperfine-shifted resonances revealed numerous NOE contacts between these signals and those in the main envelope of the proton spectrum. With the aid of the X-ray crystal structure [Carter, C.W., Kraut, J., Freer, S. T., Xuong, N. H., Alden, R. A., & Bartsch, R. G. (1974) J. Biol. Chem. 249, 4212], these NOEs permitted seven of the nine hyperfine-shifted signals to be assigned to three of the cysteine residues liganded to the metal cluster (Cys-43, Cys-46, and Cys-77). The other two hyperfine-shifted signals produced no detectable NOEs to other resonances in the spectrum and were tentatively assigned to the remaining cysteinyl ligand (Cys-63). These assignments, in conjunction with recent theoretical models of the electronic structure of the Fe4S4 cluster [Noodleman, L. (1988) Inorg. Chem. 27, 3677; Bertini, I., Briganti, F., Luchinat, C., Scozzafava, A., & Sola, M. (1991) J. Am. Chem. Soc. 113, 1237], indicate that the iron atoms coordinated to Cys-63 and Cys-77 are those of the mixed-valence Fe(3+)-Fe2+ pair whereas Cys-43 and Cys-46 are ligands to the Fe(3+)-Fe3+ metal pair.

Subunit interactions change the heme active-site geometry in p-cresol methylhydroxylase

The enzyme p-cresol methylhydroxylase [4-cresol: (acceptor) oxidoreductase (methyl-hydroxylating), EC 1.17.99.1] contains two subunits: a cytochrome c (electron transfer) subunit (cytochrome cpc) and a flavin (catalytic) subunit. When these subunits are separated by isoelectric focusing, a stable cytochrome subunit is obtained. Significant differences are observed between the one-dimensional NMR spectra of oxidized cytochrome cpc and of oxidized p-cresol methylhydroxylase. Analysis of the two-dimensional nuclear Overhauser enhancement and exchange spectroscopy (NOESY) spectrum of reduced cytochrome cpc suggests that the axial ligand, Met-50, of the stable subunit reorients by a rotation about the C gamma-S delta bond when cytochrome cpc binds to the flavin subunit. This reorientation must result in a change in bonding at the heme, which is reflected both in the para-magnetically shifted resonances and in the redox potential. p-Cresol methylhydroxylase thereby provides an interesting example of the coupling of subunit interactions to active-site structure and reactivity.

1H NMR study of the solution molecular and electronic structure of Escherichia coli ferricytochrome b562: evidence for S = 1/2 in equilibrium S = 5/2 spin equilibrium for intact His/Met ligation

The solution 500-MHz 1H NMR spectral parameters for ferricytochrome b562, a soluble 12-kDa electron carrier from Escherichia coli with axial His/Met coordination, are shown to be strongly influenced by protein concentration and ionic strength at low pH and 25 degrees C in a manner consistent with significant aggregation at low ionic strength. At high ionic strength a well-resolved 1H NMR spectrum reveals over 40 hyperfine-shifted resonances which arise from two isomeric species in the ratio 2:1. 2D COSY and NOESY maps at 25 degrees C for the hyperfine-shifted resonances allow the assignment of a number of axial His resonances and all heme peripheral substituent peaks. The resulting asymmetric heme contact shift patterns, together with the halving of the number of lines when reconstituting with 2-fold symmetric hemin, demonstrate the molecular basis of the solution heterogeneity to be heme orientational disorder. The strongly upfield-shifted axial Met-7 resonances, characteristic of low-spin ferricytochromes c with His/Met ligation, appear upfield only at very low temperatures. At elevated temperatures, all resonances, in particular those of the axial Met, move strongly downfield. Detailed analysis of the deviation from Curie behavior for different functional groups demonstrates the presence of a low spin in equilibrium high spin equilibrium with an intact His-Fe-Met coordination. The weaker axial field in ferricytochrome b562, relative to the purely low-spin ferricytochromes c, is attributed to a perturbed iron-Met bond. The contact shifts for a coordinated Met in the high-spin state are estimated. A link between equatorial hemin and axial ligand interactions is indicated by a differential population of the high-spin form for the two hemin orientations.

1H NMR study of the role of heme carboxylate side chains in modulating heme pocket structure and the mechanism of reconstitution of cytochrome b5

1H nuclear magnetic resonance spectroscopy was used to assign the hyperfine-shifted resonances and determine the position of a side chain in the heme cavity of wild-type rat apocytochrome b5 reconstituted with a series of synthetic hemins possessing systematically perturbed carboxylate side chains. The hemins included protohemin derivatives with individually removed or pairwise shortened and lengthened carboxylate side chains, as well as (propionate)n(methyl)8-nporphine-iron(III) isomers with n = 1-3 designed to force occupation of nonnative propionate sites. The resonance assignments were effected on the basis of available empirical heme contact shift correlations and steady-state nuclear Overhauser effect measurements in the low-spin oxidized proteins. The failure to detect holoproteins with certain hemins dictates that the stable holoproteins, unlike the case of myoglobin, demand the axial iron-His bonds and cannot accommodate carboxylate side chains at interior positions in the binding pocket. Hence, the heme pocket interior in cytochrome b5 is judged much less polar and less sterically accommodating than that of myoglobin. The propionate occupational preference was greatest as the native 7-propionate site, but also possible at the nonnative crystallographic 5-methyl or 8-methyl positions. Only for a propionate at the crystallographic 8-methyl position was a significant perturbation of the native molecular/electronic structure observed, and this was attributed to an alternative propionate-protein hydrogen bond at the crystallographic 8-methyl position. The structures of the transient protein complexes detected only shortly after reconstitution reveal that the initial encounter complexes during assembly of holoprotein from apoprotein and hemin involve one of the two alternate propionate-protein links at either the 7-propionate or native 8-methyl position. In a monopropionate hemin, this leads to the characterization of a new type of heme orientational disorder involving rotation about a N-Fe-N axis.

The formation of protein complexes between ferricytochrome b5 and ferricytochrome c studied using high-resolution 1H-NMR spectroscopy

The association of the tryptic fragment of bovine microsomal cytochrome b5 with cytochrome c has been studied by one- and two-dimensional 1H-NMR spectroscopy. The association of cytochromes to form protein complexes is apparent from the increase in linewidths for resonances of ferricytochrome b5 as well as small perturbations in their chemical shifts that occur upon increasing the cytochrome c/b5 molar ratio. The changes in the chemical shifts of hyperfine shifted resonances of ferricytochrome b5 with increasing ratios of ferricytochrome c indicate the formation of binary 1:1 complexes and ternary 1:2 complexes. Similarly, titrations of the linewidth of resolved resonances of ferricytochrome b5 are consistent with stoichiometries of 1:1 and 1:2 for complexes formed between cytochromes b5 and c. Surprisingly, in the 1:1 complex, mobility is shown to be a function of ionic strength. Two-dimensional correlated spectroscopy (COSY) and nuclear Overhauser enhancement spectroscopy (NOESY) of the binary complex formed between ferricytochrome b5 and c indicate that the positions of many resonances attributable to amino acids are unaltered by protein association, although distinctive chemical shift changes are detected in the alpha-CH of the haem C17 propionate. The protein complex detected by NMR is discussed with respect to the model for the binary complex proposed by Salemme and possible mechanisms of electron transfer.

An analysis of pseudocontact shifts and their relationship to structural features of the redox states of cytochrome b5

The assignment of proton resonances in both redox states of a heme protein is necessary for the evaluation of pseudocontact shift data. Many new assignments are presented here for cytochrome b5, particularly in the paramagnetic oxidised state, thereby allowing both the calculation of electronic g-tensor values with the magnetic axis orientation and a comparison of observed and calculated pseudocontact shifts utilising a computational procedure. The possible redox linked conformational changes are found to be minimal in contrast with cytochrome c although the procedure additionally highlights aspects of the mobility of certain residues in cytochrome b5. In this respect the residue Gly-42 appears mobile both by this method and by the observation from NMR spectra of a major and minor conformation in this region.

Structural properties of apocytochrome b5: presence of a stable native core

Upon removal of the heme group, the water-soluble fragment of cytochrome b5 adopts a conformation less stable and compact than that of the holoprotein [Huntley, T. E., & Strittmatter, P. (1972) J. Biol. Chem. 247, 4641-4647]. This conformation, imposed by the amino acid sequence alone, has not been described in detail. One- and two-dimensional proton nuclear magnetic resonance spectroscopy techniques were applied to the apoprotein of the soluble fragment of rat liver cytochrome b5 in an effort to characterize the structure of the apoprotein. Nuclear Overhauser spectroscopy revealed a number of short interresidue distances and demonstrated that, in spite of the increased flexibility, at least one cluster of side chains exists on a time scale long enough for study. Several residues participating in the cluster, in particular the only Trp (Trp 22), were identified. Similarities with the spectrum of the reduced holoprotein were observed that led to the inspection of the cytochrome b5 crystal structure for assigning resonances. It appeared that the environment of this residue maintains its integrity in the apoprotein. Since in the holoprotein Trp 22 belongs to a hydrophobic core formed in part by beta-strands, it is proposed that some of this beta-structure is stable in the absence of the heme-protein interactions. Implications for structure and folding are discussed.

Solution structural characteristics of cyanometmyoglobin: resonance assignment of heme cavity residues by two-dimensional NMR

Steady-state nuclear Overhauser effects (NOE), two-dimensional (2D) nuclear Overhauser effect spectroscopy (NOESY), and 2D spin correlation spectroscopy (COSY) have been applied to the fully paramagnetic low-spin, cyanide-ligated complex of sperm whale ferric myoglobin to assign the majority of the heme pocket side-chain proton signals and the remainder of the heme signals. It is shown that the 2D NOESY map reveals essentially all dipolar connectivities observed in ordinary 1D NOE experiments and expected on the basis of crystal coordinates, albeit often more weakly than in a diamagnetic analogue. For extremely broad (approximately 600-Hz) and rapidly relaxing (Tf1 approximately 3 ms) signals which show no NEOSY peaks, we demonstrate that conventional steady-state NOEs obtained under very rapid pulsing conditions still allow detection of the critical dipoar connectivities that allow unambiguous assignments. The COSY map was found to be generally less useful for the hyperfine-shifted residues, with cross peaks detected only for protons greater than 6 A from the iron. Nevertheless, numerous critical COSY cross peaks between strongly hyperfine-shifted peaks were resolved and assigned. In all, 95% (53 of 56 signals) of the total proton sets within approximately 7.5 A of the iron, the region experiencing the strongest hyperfine shifts and paramagnetic relaxation, are now unambiguously assigned. Hence it is clear that the 2D methods can be profitably applied to paramagnetic proteins. The scope and limitations of such application are discussed. The resulting hyperfine shift pattern for the heme confirmed expectations based on model compounds. In contrast, while exhibiting fortuitous 1H NMR spectral similarities, a major discrepancy was uncovered between the hyperfine shift pattern of the axially bound (F8 histidyl) imidazole in the protein and that of the imidazole in a relevant model compound [Chacko, V.P., & La Mar, G. N. (1982) J. Am. Chem. Soc. 104, 7002-7007], providing direct evidence for a protein-based deformation of axial bonding in the protein.

Structural studies of cytochrome b5: complete sequence-specific resonance assignments for the trypsin-solubilized microsomal ferrocytochrome b5 obtained from pig and calf

We report complete sequence-specific proton resonance assignments for the trypsin-solubilized microsomal ferrocytochrome b5 obtained from calf liver. In addition, sequence-specific resonance assignments for the main-chain amino acid protons (i.e., C alpha, C beta, and amide protons) are also reported for the porcine cytochrome b5. Assignment of the majority of the main-chain resonances was rapidly accomplished by automated procedures that used COSY and HOHAHA peak coordinates as input. Long side chain amino acid spin system identification was facilitated by long-range coherence-transfer experiments (HOHAHA). Problems with resonance overlap were resolved by examining differences between the two-dimensional 500-MHz NMR spectra of rabbit, pig, and calf proteins and by examining the temperature-dependent variation of amide proton resonances. Calculations of the aromatic ring-current shifts for protons that the X-ray crystal structure indicated were proximal to aromatic residues were found to be useful in corroborating assignments, especially those due to the large shifts induced by the heme. Assignment of NOESY cross peaks was greatly facilitated by a prediction of intensities using a complete relaxation matrix analysis based on the crystal structure. These results suggest that the single-crystal X-ray structure closely resembles that of the solution structure although there is evidence that the solution structure has a more dynamic character.

Determination of the functionally important heme peripheral vinyl group orientation in paramagnetic hemoprotein by 2D NMR

2D NMR spectroscopies have been successfully used to characterize the heme peripheral vinyl groups in paramagnetic hemoprotein in spite of the difficulties from the rapid paramagnetic relaxation and the low digital resolution of the 2D NMR map. The scalar coupling network system among the vinyl protons is clearly identified in the COSY spectra from its characteristic cross-peak pattern and the dipolar coupling connectivities of the vinyl proton resonances with other heme side-chain proton resonances not only provide the specific assignment of vinyl beta-proton resonances but also allow the determination of the vinyl group orientation with respect to the heme plane.