Estimation of microscopic redox potentials of a tetraheme protein, cytochrome c3 of Desulfovibrio vulgaris, Miyazaki F, and partial assignments of heme groups

Exploring Coupled Redox and pH Processes with a Force-Field-Based Approach: Applications to Five Different Systems

Coupled redox and pH-driven processes are at the core of many important biological mechanisms. As the distribution of protonation and redox states in a system is associated with the pH and redox potential of the solution, having efficient computational tools that can simulate under these conditions becomes very important. Such tools have the potential to provide information that complement and drive experiments. In previous publications we have presented the implementation of the constant pH and redox potential molecular dynamics (C(pH,E)MD) method in AMBER and we have shown how multidimensional replica exchange can be used to significantly enhance the convergence efficiency of our simulations. In the current work, after an improvement in our C(pH,E)MD approach that allows a given residue to be simultaneously pH- and redox-active, we have employed our methodologies to study five different systems of interest in the literature. We present results for capped tyrosine dipeptide, two maquette systems containing one pH- and redox-active tyrosine (α₃Y and peptide A), and two proteins that contain multiple heme groups (diheme cytochrome c from Rhodobacter sphaeroides and Desulfovibrio vulgaris Hildenborough cytochrome c₃). We show that our results can provide new insights into previous theoretical and experimental findings by using a fully force-field-based and GPU-accelerated approach, which allows the simulations to be executed with high computational performance.

Design and fine-tuning redox potentials of metalloproteins involved in electron transfer in bioenergetics

Redox potentials are a major contributor in controlling the electron transfer (ET) rates and thus regulating the ET processes in the bioenergetics. To maximize the efficiency of the ET process, one needs to master the art of tuning the redox potential, especially in metalloproteins, as they represent major classes of ET proteins. In this review, we first describe the importance of tuning the redox potential of ET centers and its role in regulating the ET in bioenergetic processes including photosynthesis and respiration. The main focus of this review is to summarize recent work in designing the ET centers, namely cupredoxins, cytochromes, and iron-sulfur proteins, and examples in design of protein networks involved these ET centers. We then discuss the factors that affect redox potentials of these ET centers including metal ion, the ligands to metal center and interactions beyond the primary ligand, especially non-covalent secondary coordination sphere interactions. We provide examples of strategies to fine-tune the redox potential using both natural and unnatural amino acids and native and nonnative cofactors. Several case studies are used to illustrate recent successes in this area. Outlooks for future endeavors are also provided. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Analysis of the electrochemistry of hemes with E(m)s spanning 800 mV

The free energy of heme reduction in different proteins is found to vary over more than 18 kcal/mol. It is a challenge to determine how proteins manage to achieve this enormous range of E(m)s with a single type of redox cofactor. Proteins containing 141 unique hemes of a-, b-, and c-type, with bis-His, His-Met, and aquo-His ligation were calculated using Multi-Conformation Continuum Electrostatics (MCCE). The experimental E(m)s range over 800 mV from -350 mV in cytochrome c(3) to 450 mV in cytochrome c peroxidase (vs. SHE). The quantitative analysis of the factors that modulate heme electrochemistry includes the interactions of the heme with its ligands, the solvent, the protein backbone, and sidechains. MCCE calculated E(m)s are in good agreement with measured values. Using no free parameters the slope of the line comparing calculated and experimental E(m)s is 0.73 (R(2) = 0.90), showing the method accounts for 73% of the observed E(m) range. Adding a +160 mV correction to the His-Met c-type hemes yields a slope of 0.97 (R(2) = 0.93). With the correction 65% of the hemes have an absolute error smaller than 60 mV and 92% are within 120 mV. The overview of heme proteins with known structures and E(m)s shows both the lowest and highest potential hemes are c-type, whereas the b-type hemes are found in the middle E(m) range. In solution, bis-His ligation lowers the E(m) by approximately 205 mV relative to hemes with His-Met ligands. The bis-His, aquo-His, and His-Met ligated b-type hemes all cluster about E(m)s which are approximately 200 mV more positive in protein than in water. In contrast, the low potential bis-His c-type hemes are shifted little from in solution, whereas the high potential His-Met c-type hemes are raised by approximately 300 mV from solution. The analysis shows that no single type of interaction can be identified as the most important in setting heme electrochemistry in proteins. For example, the loss of solvation (reaction field) energy, which raises the E(m), has been suggested to be a major factor in tuning in situ E(m)s. However, the calculated solvation energy vs. experimental E(m) shows a slope of 0.2 and R(2) of 0.5 thus correlates weakly with E(m)s. All other individual interactions show even less correlation with E(m). However the sum of these terms does reproduce the range of observed E(m)s. Therefore, different proteins use different aspects of their structures to modulate the in situ heme electrochemistry. This study also shows that the calculated E(m)s are relatively insensitive to different heme partial charges and to the protein dielectric constant used in the simulation.

Strategic roles of axial histidines in structure formation and redox regulation of tetraheme cytochrome c3

Tetraheme cytochrome c 3 (cyt c 3) exhibits extremely low reduction potentials and unique properties. Since axial ligands should be the most important factors for this protein, every axial histidine of Desulfovibrio vulgaris Miyazaki F cyt c 3 was replaced with methionine, one by one. On mutation at the fifth ligand, the relevant heme could not be linked to the polypeptide, revealing the essential role of the fifth histidine in heme linking. The fifth histidine is the key residue in the structure formation and redox regulation of a c-type cytochrome. A crystal structure has been obtained for only H25M cyt c 3. The overall structure was not affected by the mutation except for the sixth methionine coordination at heme 3. NMR spectra revealed that each mutated methionine is coordinated to the sixth site of the relevant heme in the reduced state, while ligand conversion takes place at hemes 1 and 4 during oxidation at pH 7. The replacement of the sixth ligand with methionine caused an increase in the reduction potential of the mutated heme of 222-244 mV. The midpoint potential of a triheme H52M cyt c 3 is higher than that of the wild type by approximately 50 mV, suggesting a contribution of the tetraheme architecture to the lowering of the reduction potentials. The hydrogen bonding of Thr24 with an axial ligand induces a decrease in reduction potential of approximately 50 mV. In conclusion, the bis-histidine coordination is strategically essential for the structure formation and the extremely low reduction potential of cyt c 3.

Roles of charged residues in pH-dependent redox properties of cytochrome c3 from Desulfovibrio vulgaris Miyazaki F

Complicated pH-properties of the tetraheme cytochrome c3 (cyt c3) from Desulfovibrio vulgaris Miyazaki F (DvMF) were examined by the pH titrations of 1H-15N HSQC spectra in the ferric and ferrous states. The redox-linked pKa shift for the propionate group at C13 of heme 1 was observed as the changes of the NH signals around it. This pKa shift is consistent with the redox-linked conformational alteration responsible for the cooperative reduction between hemes 1 and 2. On the other hand, large chemical shift changes caused by the protonation/deprotonation of Glu41 and/or Asp42, and His67 were redox-independent. Nevertheless, these charged residues affect the redox properties of the four hemes. Furthermore, one of interesting charged residues, Glu41, was studied by site-directed mutagenesis. E41K mutation increased the microscopic redox potentials of heme 1 by 46 and 34 mV, and heme 2 by 35 and 30 mV at the first and last reduction steps, respectively. Although global folding in the crystal structure of E41K cyt c3 is similar to that of wild type, local change was observed in 1H NMR spectrum. Glu41 is important to keep the stable conformation in the region between hemes 1 and 2, controlling the redox properties of DvMF cyt c3. In contrast, the kinetic parameters for electron transfer from DvMF [NiFe] hydrogenase were not influenced by E41K mutation. This suggests that the region between hemes 1 and 2 is not involved in the interaction with [NiFe] hydrogenase, and it supports the idea that heme 4 is the exclusive entrance gate to accept the electron in the initial reduction stage.

Solution structures of tetrahaem ferricytochrome c3 from Desulfovibrio vulgaris (Hildenborough) and its K45Q mutant: The molecular basis of cooperativity

The NMR structure of the oxidised wild-type cytochrome c3 from Desulfovibrio vulgaris Hildenborough was determined in solution. Using a newly developed methodology, NMR data from the K45Q mutant was then grafted onto data from the wild-type protein to determine the structure in the region of the mutation. The structural origins of the redox-Bohr effect and haem-haem cooperativities are discussed with respect to the redox-related conformational changes observed in solution.

Phospholipid-dependent regulation of cytochrome c3-mediated electron transport across membranes

Cytochrome c3 (cyt c3) can mediate electron transport across phosphatidylcholine (PC)/cardiolipin (CL) and PC/phosphatidylglycerol (PG) membranes. A two-molecule process is involved in the electron transport across PC/CL membranes in the liquid-crystalline state. In contrast, a single-molecule process dominates the electron transport across PC/CL membranes in the gel state and PC/PG membranes in the liquid-crystalline and gel states. Namely, the electron transport mechanism differs with the phospholipid composition and membrane fluidity. The rate-limiting step of the two-molecule process was lateral diffusion of cyt c3 in membranes. The rate constants for the three single-molecule process cases were similar to each other. To elucidate these reaction processes, interactions between cyt c3 and phosphate groups and between cyt c3 and the glycerol backbones of phospholipid bilayers were investigated by means of 31P and 2H solid-state NMR, respectively, for CL and PC/CL membranes. The results showed that the polar headgroups of both phosphatidylcholine and CL are involved in the binding of cyt c3. Also, cyt c3 penetrates into membranes, which would induce distortion of the lipid bilayer. The molecular mechanisms underlying the single- and two-molecule processes are discussed in terms of membrane structure.

Role of the aromatic ring of Tyr43 in tetraheme cytochrome c(3) from Desulfovibrio vulgaris Miyazaki F

Tyrosine 43 is positioned parallel to the fifth heme axial ligand, His34, of heme 1 in the tetraheme cytochrome c(3). The replacement of tyrosine with leucine increased the redox potential of heme 1 by 44 and 35 mV at the first and last reduction steps, respectively; its effects on the other hemes are small. In contrast, the Y43F mutation hardly changed the potentials. It shows that the aromatic ring at this position contributes to lowering the redox potential of heme 1 locally, although this cannot be the major contribution to the extremely low redox potentials of cytochrome c(3). Furthermore, temperature-dependent line-width broadening in partially reduced samples established that the aromatic ring at position 43 participates in the control of the kinetics of intramolecular electron transfer. The rate of reduction of Y43L cytochrome c(3) by 5-deazariboflavin semiquinone under partially reduced conditions was significantly different from that of the wild type in the last stage of the reduction, supporting the involvement of Tyr43 in regulation of reduction kinetics. The mutation of Y43L, however, did not induce a significant change in the crystal structure.

Redox-coupled conformational alternations in cytochrome c(3) from D. vulgaris Miyazaki F on the basis of its reduced solution structure

Heteronuclear NMR spectroscopy was performed to determine the solution structure of (15)N-labeled ferrocytochrome c(3) from Desulfovibrio vulgaris Miyazaki F (DvMF). Although the folding of the reduced cytochrome c(3) in solution was similar to that of the oxidized one in the crystal structure, the region involving hemes 1 and 2 was different. The redox-coupled conformational change is consistent with the reported solution structure of D. vulgaris Hildenborough ferrocytochrome c(3), but is different from those of other cytochromes c(3). The former is homologous with DvMF cytochrome c(3) in amino acid sequence. Small displacements of hemes 1 and 2 relative to hemes 3 and 4 were observed. This observation is consistent with the unusual behavior of the 2(1)CH(3) signal of heme 3 reported previously. As shown by the (15)N relaxation parameters of the backbone, a region between hemes 1 and 2 has more flexibility than the other regions. The results of this work strongly suggest that the cooperative reduction of hemes 1 and 2 is based on the conformational changes of the C-13 propionate of heme 1 and the aromatic ring of Tyr43, and the interaction between His34 and His 35 through covalent and coordination bonds.

Structure-function relationships in heme-proteins

Biological systems rely on heme-proteins to carry out a number of basic functions essential for their survival. Hemes, or iron-porphyrin complexes, are the versatile and ubiquitous active centers of these proteins. In the past decade, discovery of new heme-proteins, together with functional and structural research, provided a wealth of information on these diverse and biologically important molecules. Structure determination work has shown that nature has used a variety of different scaffolds and architectures to bind heme and modulate functions such as redox properties. Structural data have also provided insights into the heme-linked protein conformational changes required in many regulatory heme-proteins. Remarkable efforts have been made towards the understanding of factors governing redox potentials. Site-directed mutagenesis studies and theoretical calculations on heme environments investigated the roles of hydrophobic and electrostatic residues, and analyzed the effect of heme solvent accessibility. This review focuses on the structure-function relationships underlying the association of heme in signaling and iron metabolism proteins. In addition, an account is given about molecular features affecting heme's redox properties; this briefly revisits previous conclusions in the light of some more recent reports.

Ionic strength-dependent physicochemical factors in cytochrome c3 regulating the electron transfer rate

The effect of ionic strength on the macroscopic and microscopic redox potentials and the heme environment of cytochrome c3 from Desulfovibrio vulgaris Miyazaki F have been investigated by NMR and electrochemical methods. The redox potentials of this tetraheme protein are found to be ionic strength-dependent. Especially, the microscopic redox potentials of hemes 2 and 3 at the fourth reduction step increase significantly with increasing ionic strength, which is in contraction to the theoretical expectation. The coordinated imidazole proton signals are unaffected by ionic strength. However, the methyl and propionate proton signals of hemes 1 and 4 showed significant ionic strength dependencies that are distinct from those for hemes 2 and 3. This heme classification is the same as that found in the ionic strength dependencies of the microscopic redox potentials at the fourth reduction step. Furthermore, the effect of ionic strength on the electrostatic potentials at the heme irons has been examined on the theoretical basis. The electrostatic potential at heme 4 changes up to 1 M ionic strength, which was not expected from the observations reported on cytochromes so far. These results are discussed in connection with the reported anomalous ionic strength dependency of the reduction rate of cytochrome c3.

Use of paramagnetic NMR probes for structural analysis in cytochrome c3 from Desulfovibrio vulgaris

The dipolar field generated by each of the four haems in the tetrahaem ferricytochrome c3 from Desulfovibrio vulgaris (Hildenborough) (c3DvH) is determined by means of a novel procedure. In this method the 13C chemical shifts of the nuclei directly bound to the haems are used to determine the in-plane orientations of the rhombic perturbation in each of the four haems with respect to a model of molecular orbitals of e(g) symmetry which are subject to a rhombic perturbation [Turner, D. L., Salgueiro, C. A., Schenkels, P., LeGall, J. & Xavier, A. V. (1995) Biochim. Biophys. Acta 1246, 24-281. These orientations, together with the components of the magnetic susceptibility tensors obtained from the EPR g values and the crystal structure of c3DvH, can be used to calculate the dipolar shifts induced by each haem throughout the protein. Thus the observed 13C paramagnetic shifts of the c3DvH haem substituents were fitted considering both the pseudocontact and contact shifts of each haem simultaneously. The dipolar shifts calculated by this method were tested against the observed dipolar shifts for some amino acid residues strategically placed in the protein and also for the haem propionate groups. The effect of considering the calculated dipolar extrinsic shifts on the behaviour of the chemical shifts of the haem methyl groups in the intermediate stages of oxidation at different pH values was also analysed. The several tests applied to the calculated dipolar shifts have shown that the method is extremely useful for predicting chemical shifts as an aid to complete proton assignment, and to add further constraints in the refinement of solution structures of paramagnetic proteins and hence to probe subtle structural rearrangements around the haem pocket.

NMR studies of cooperativity in the tetrahaem cytochrome c3 from Desulfovibrio vulgaris

The thermodynamic properties of the Desulfovibrio vulgaris (Hildenborough) tetrahaem cytochrome c3 (Dvc3) are rationalised by a model which involves both homotropic (e-/e-) and heterotropic (e-/H+) cooperativity. The paramagnetic shifts of a methyl group from each haem of the Dvc3 have been determined in each stage of oxidation at several pH values by means of two-dimensional exchange NMR. The thermodynamic parameters are obtained by fitting the model to the NMR data and to redox titrations followed by visible spectroscopy. They show significant positive cooperativity between two of the haems whereas the remaining interactions appear to be largely electrostatic in origin. These parameters imply that the protein undergoes a proton-assisted two-electron transfer which can be used for energy transduction. Comparison with the crystal structure together with measurement of the kinetics of proton exchange suggest that the pH dependence is mediated by a charged residue(s) readily acessible to the solvent and close to haem I.

Structural and functional characterization of cytochrome c3 from D. desulfuricans ATCC 27774 by 1H-NMR

Cooperativity between redox and protonation centres is known to be crucial for the function of complex proteins, but it is often difficult to describe in terms of thermodynamic parameters. Cytochrome c3 is a good model for these studies since, while retaining the overall complexity of larger systems, it is suitable for detailed crystallographic and spectroscopic studies. Assignment of the haem substituent NMR resonances, together with NMR redox titrations of cytochrome c3 from D. desulfuricans ATCC 27774, was used to correlate relative redox potentials to specific haems in the structure: haem II approximately equal to haem I < haem IV < haem III. This order is different from that determined for the homologous proteins studied and in disagreement with that previously reported for this cytochrome (Morais, J., Palma, N., Frazäo, C., Caldeira, J., LeGall, J., Moura, I., Moura, J.J.G. and Carrondo, M.A. (1995) Biochemistry 34, 12830-12841).

Regulation of the redox order of four hemes by pH in cytochrome c3 from D. vulgaris Miyazaki F

The assignment of 1H-NMR signals of the heme methyl and propionate groups of cytochrome c3 of D. vulgaris Miyazaki F was performed. The heme assignment was revised for hemes 2 and 3 (sequential heme numbering). Namely, heme 4 is mainly reduced at first with hemes 1, 2 and 3 following it in this order. The p2H titration of heme methyl signals in four macroscopic oxidation states was performed in the p2H range of 5.2 to 9.0. While the heme methyl resonances in the fully oxidized state showed just small changes with p2H, most resonances in the intermediate oxidation states revealed clear p2H dependence. In particular, the methyl resonances of heme 1 shifted significantly in the acidic region. Then, the chemical shifts of beta-CH2 (next to the carboxyl group) of all propionate groups in the fully oxidized state were observed at various p2H in the range of 4.5 to 9.0. Only the propionate group at C-13 (IUPAC-IUB nomenclature) of heme 1 showed a clear change in this p2H range, its titration curve being similar to those of the methyl resonances of heme 1 in the intermediate oxidation states. pKa of the propionate group was 5.95 +/- 0.05. Analysis of the microscopic formal redox potentials was carried out for the observations at p2H 5.2, 7.1 and 9.0. The redox potentials of heme 1 showed the most remarkable p2H dependence, resulting in the change of the order of the redox potentials of four hemes. A significant change was also found in the interacting potential between hemes 1 and 2. In the light of the p2H-titration experiments, the propionate at C-13 of heme 1 was identified as the most plausible ionizable group responsible for the p2H dependence of microscopic redox potentials of heme 1 in the acidic region.

Control of the redox potential in c-type cytochromes: importance of the entropic contribution

The enthalpic and entropic components of the redox free energy variation of cytochrome c553 from Desulfovibrio vulgaris Hildenborough and its mutant Y64V, flavocytochrome b2 from Saccharomyces cerevisiae, and the different hemes of cytochromes c3 from Desulfovibrio vulgaris Miyazaki and Desulfovibrio desulfuricans Norway have been determined in 0.1 M Tris-HCl pH 7.0 (7.6 for cytochromes c3) at 25 degrees C by using nonisothermal potentiometric titrations. The set of available experimental data demonstrates that the entropic component plays an important role in the control of the redox potential in c-type and b-type cytochromes. The variation of the entropic component within the class of cytochromes characterized by a positive value of E degrees ' is proposed to be mainly determined by the variation of the exposure of the heme propionates to the solvent. In the case of tetraheme cytochromes c3, the thermodynamic characteristics vary largely among the hemes belonging to the same molecule, which reflects the environmental peculiarities of each heme and also the heme-heme redox interactions. This study substantiates the existence of compensatory effects between large and opposite contributions to E degree ' predicted by all the current theoretical models which are based on electrostatic free energy calculations.

Crystal structure of cytochrome c3 from Desulfovibrio desulfuricans Norway at 1.7 A resolution

The crystal structure of cytochrome c3 (M(r) 13,000) from Desulfovibrio desulfuricans (118 residues, four heme groups) has been crystallographically refined to 1.7 A resolution using a simulated annealing method, based on the structure-model at 2.5 A resolution, already published. The final R-factor for 10,549 reflections was 0.198 covering the range from 5.5 to 1.7 A resolution. The individual temperature factors were refined for a total of 1059 protein atoms, together with 126 bound solvent molecules. The structure has been analyzed with respect to its detailed conformational properties, secondary structure features, temperature factor behaviour, bound solvent sites and heme geometry and ligation. The characteristic secondary structures of the polypeptide chain of this molecule are one extended alpha-helix, a short beta-strand and 13 reverse turns. The four heme groups are located in different structural environments, all highly exposed to solvent. The particular structural features of the heme environments are compared to the four hemes of the cytochrome c3 from Desulfovibrio vulgaris Miyazaki.

Homotropic and heterotropic cooperativity in the tetrahaem cytochrome c3 from Desulfovibrio vulgaris

The thermodynamic parameters which govern the homotropic (e-/e-) and heterotropic (e-/H+) cooperativity in the tetrahaem cytochrome c3 isolated from Desulfovibrio vulgaris (Hildenborough) were determined, using the paramagnetic shifts of haem methyl groups in the NMR spectra of intermediate oxidized states at different pH levels. A model is put forward to explain how the network of positive and negative cooperativities between the four haems and acid/base group(s) enables the protein to achieve a proton-assisted 2e- step.

Design and synthesis of multi-haem proteins

A water-soluble, 62-residue, di-alpha-helical peptide has been synthesized which accommodates two bis-histidyl haem groups. The peptide assembles into a four-helix dimer with 2-fold symmetry and four parallel haems that closely resemble native haems in their spectral and electrochemical properties, including haem-haem redox interaction. This protein is an essential intermediate in the synthesis of molecular 'maquettes', a novel class of simplified versions of the metalloproteins involved in redox catalysis and in energy conversion in respiratory and photosynthetic electron transfer.

Individual redox characteristics and kinetic properties of the hemes in cytochromes c3: new methods of investigation

The elucidation of the role of the four hemes in cytochromes c3 requires several complementary approaches. The measurements and the assignment of the redox potentials resort to magnetic spectroscopies, EPR and NMR, which are able to discriminate the hemes. The origin of the differences between the redox properties of the hemes can be studied by comparing their thermodynamic parameters delta S and delta H, as measured by the temperature dependence of their individual potentials. Lastly, the available data concerning the electron exchange between cytochromes c3 and their redox partners can be analysed through a detailed kinetic model which provides important information on the role of the different hemes.

Molecular and structural basis of electron transfer in tetra- and octa-heme cytochromes

The first three-dimensional structure of a dimeric, octa-heme cytochrome c3 (M(r) 26000) from Desulfovibrio desulfuricans Norway, established at 2.2 A resolution, is briefly presented and compared to the known 3-D-structures of different C3-type tetraheme cytochromes, in order to contribute to a better understanding of the function of multiheme clusters and of the role of conserved amino acids implicated in possible electron transfer pathways. The dimeric protein crystallizes in the space group P3(1)21 with a = 73.01 A, c = 61.81 A and the asymmetric unit contains one monomer subunit, the dimer being generated by the crystallographic two-fold axis. The 3-D-structure was solved using the molecular replacement method with a model based on the structure of the tetraheme cytochrome c3 (M(r) 13000) from D desulfuricans Norway, presently refined at 1.7 A resolution. The monomeric subunit has the same overall fold as all cytochromes c3 (M(r) 13000). Moreover, the heme core of all examined cytochromes c3 is highly conserved, but differences appear concerning the heme environments and the histidines, axial ligands of the heme-iron atoms.

Involvement of electrostatic interactions in cytochrome c complex formations

Structural studies on various electron transfer complexes involving the tetrahemic cytochrome c3 provided evidence that one of the hemes (heme 4) is the interacting site on the molecule. The reactivity of this particular heme is allocated to the positive charges found around the heme group which are strongly involved in the electrostatic interaction processes. Electrostatic and hydrophobic effects in complex formation are considered on the basis of two electron transfer complex examples: the soluble cytochrome c-cytochrome c peroxidase and the membrane bound photosynthetic reaction center.

A detailed comparison of the refined structures of cytochrome c3 molecules from two strains in Desulfovibrio vulgaris: the relationship between the heme structures and their redox properties

The refined structures of the cytochrome c3 molecules from Desulfovibrio vulgaris Miyazaki and Hildenborough have been compared in detail. Though there are no significant differences of the overall structure and spatial arrangement of four heme groups between two molecules, there are some unique features with regard to the local structures near the heme pockets. Two of the heme groups show significant differences including a hydrogen bonding scheme between the imidazole rings and water molecules in heme pocket. The water molecule near the heme 1 pocket in Miyazaki cytochrome c3 is replaced by the C delta 1 atom of Leu9 in the Hildenborough structure. The rotation around the C alpha-C beta bond of Leu9 results in the increase of the solvent accessibility of the portion of the heme 1 edge. On the other hand, in heme3 one water molecule is additionally included in the hydrogen bonding network with a N delta 1 atom of the imidazole ring and carbonyl oxygen atom of the main chain in Hildenborough cytochrome c3. This additional hydrogen bonding network will contribute to the possible passage between the imidazole ring and the molecular surface. These features have been related to the upfield shift of the C2 proton signals of NMR spectra and the low redox potentials of the cytochrome c3 molecules of Desulfovibrio vulgaris Hildenborough in comparison with those from Miyazaki.

Nuclear magnetic resonance studies of cytochromes c in solution

Cytochromes c are small soluble proteins, which have been extensively studied by nuclear magnetic resonance spectroscopy. The specific NMR features of paramagnetic proteins are discussed for the oxidized form (paramagnetic shift and line broadening). Early NMR studies have focused on the electronic structure of the heme and its direct environment. The conformations of cytochromes c are now investigated by two-dimensional 1H NMR spectroscopy combined with restrained molecular dynamics. 15N and 13C NMR, which greatly benefit from isotopic enrichment, may help in obtaining reliable 1H assignments and thus high quality solution structure. Finally, hydrogen exchange rates provide insight in the rigidity (and stability) of cytochromes c in both redox states at the atomic level.

Characterization of the structure and redox behaviour of cytochrome c3 from Desulfovibrio baculatus by 1H-nuclear-magnetic-resonance spectroscopy

Complete assignment of the aromatic and haem proton resonances in the cytochromes c3 isolated from Desulfovibrio baculatus strains (Norway 4, DSM 1741) and (DSM 1743) was achieved using one- and two-dimensional 1H n.m.r. Nuclear Overhauser enhancements observed between haem and aromatic resonances and between resonances due to different haems, together with the ring-current contributions to the chemical shifts of haem resonances, support the argument that the haem core architecture is conserved in the various cytochromes c3, and that the X-ray structure of the D. baculatus cytochrome c3 is erroneous. The relative orientation of the haems for both cytochromes was determined directly from n.m.r. data. The n.m.r. structures have a resolution of approximately 0.25 nm and are found to be in close agreement with the X-ray structure from D. vulgaris cytochrome c3. The proton assignments were used to relate the highest potential to a specific haem in the three-dimensional structure by monitoring the chemical-shift variation of several haem resonances throughout redox titrations followed by 1H n.m.r. The haem with highest redox potential is not the same as that in other cytochromes c3.

Structural studies on Desulfovibrio gigas cytochrome c3 by two-dimensional 1H-nuclear-magnetic-resonance spectroscopy

Several aromatic amino acid residues and haem resonances in the fully reduced form of Desulfovibrio gigas cytochrome c3 are assigned, using two-dimensional 1H n.m.r., on the basis of the interactions between the protons of the aromatic amino acids and the haem protons as well as the intrahaem distances known from the X-ray structure [Kissinger (1989) Ph.D. Thesis, Washington State University]. The interhaem interactions observed in the n.m.r. spectra are in full agreement with the D. gigas X-ray structure and also with the n.m.r. data from Desulfovibrio vulgaris (Hildenborough) [Turner, Salgueiro, LeGall and Xavier (1992) Eur. J. Biochem. 210, 931-936]. The good correlation between the calculated ring-current shifts and the observed chemical shifts strongly supports the present assignments. Observation of the two-dimensional nuclear-Overhauser-enhancement spectra of the protein in the reduced, intermediate and fully oxidized stages led to the ordering of the haems in terms of their midpoint redox potentials and their identification in the X-ray structure. The first haem to oxidize is haem I, followed by haems II, III and IV, numbered according to the Cys ligand positions in the amino acid sequences [Mathews (1985) Prog. Biophys. Mol. Biol. 54, 1-56]. Although the haem core architecture is the same for the different Desulfovibrio cytochromes c3, the order of redox potentials is different.

Structural studies of Desulfovibrio vulgaris ferrocytochrome c3 by two-dimensional NMR

Two-dimensional NMR has been used to make specific assignments for the four haems in Desulfovibrio vulgaris (Hildenborough) ferrocytochrome c3 and to determine their haem core architecture. The NMR signals from the haem protons were assigned according to type using two-dimensional NMR experiments which led to four sets of signals, one for each of the haems. Specific assignments were obtained by calculating the ring current shifts which arise from other haems and aromatic residues. Observation of interhaem NOEs confirmed the assignments and established that the relative orientation of the haems is identical to that found in the crystal structure of D. vulgaris (Miyazaki F.) ferricytochrome c3. Assignments were also made for all the aromatic residues except for the haem ligands and F20, which is shifted under the main envelope of signals. The NOEs observed between these aromatic protons and haem protons confirm the similarity between the structures in solution and in the crystal. The assignments reported here are the basis for the cross-assignments of the four microscopic haem redox potentials to specific haems in the protein structure [Salgueiro, C. A., Turner, D. L., Santos, H., LeGall, J. and Xavier, A. V. (1992) FEBS Lett., in the press]

Assignment of the redox potentials to the four haems in Desulfovibrio vulgaris cytochrome c3 by 2D-NMR

Using 2D-NMR the four haems of Desulfovibrio vulgaris (Hildenborough) cytochrome c3 within the X-ray structure were fully cross assigned according to their redox potential. The strategy used was based on a complete network of chemical exchange connectivities between the NMR signals obtained for all oxidation levels to the corresponding ones in the fully reduced spectrum [1992, Eur. J. Biochem., in press]. This unequivocal cross-assignment disagrees with earlier results obtained for the similar protein from Desulfovibrio vulgaris (Miyazaki F.).

Reduction kinetics of the four hemes of cytochrome c3 from Desulfovibrio vulgaris by flash photolysis

The reduction of the tetraheme cytochrome c3 (from Desulfovibrio vulgaris, strains Miyazaki F and Hildenbourough) by flavin semiquinone and reduced methyl viologen follows a monophasic kinetic profile, even though the four hemes do not have equivalent reduction potentials. Rate constants for reduction of the individual hemes are obtained subsequent to incrementally reducing the cytochrome by phototitration. The dependence of each rate constant on the reduction potential difference between the heme and the reductant can be described by outer sphere electron transfer theroy. Thus, the very low reduction potentials of the cytochrome c3 hemes compensate for the very large solvent accessibility of the hemes. The relative rate constants for electron transfer to the four hemes of cytochrome c3 are consistent with the assignments of reduction potential to hemes previously made by Park et al. (Park, J.-S., Kano, K., Niki, S. and Akutsu, H. (1991) FEBS Lett. 285, 149-151) using NMR techniques. The ionic strength dependence of the observed rate constant for reduction by the methyl viologen radical cation indicates that ionic strength substantially alters the structure and/or the heme reduction potentials of the cytochrome. This result is confirmed by reduction with a neutral flavin species (5-deazariboflavin semiquinone) in which the reactivity of the highest potential heme decreases and the reactivity of the lowest potential heme increases at high (500 mM) ionic strength, and by the sensitivity of heme methyl resonances to ionic strength as observed by 1H-NMR. These unusual ionic strength-dependent effects may be due to a combination of structural changes in the cytochrome and alterations of the electrostatic fields at elevated ionic strengths.

Revision of the haem-core architecture in the tetraheam cytochrome c3 from Desulfovibrio baculatus by two-dimensional 1H NMR

The haem-core architecture in cytochrome c3 isolated from Desulfovibrio baculatus (Norway 4) was probed using two-dimensional 1H NMR. Interhaem connectivities detected in NOE spectroscopy experiments performed at short mixing times are incompatible with the structure of the protein determined by X-ray crystallography, but agree instead with the haem arrangement found in cytochrome c3 from Desulfovibrio vulgaris (Miyazaki). These experiments show unequivocally that the relative orientation of the four haems in the two proteins is the same and does not involve the 180 degrees rotation of haems I and IV indicated in the X-ray structure determined for the cytochrome c3 from D. baculatus (Norway 4).

Non-equivalent natures of the coordinated imidazole rings of cytochrome c3 from D. vulgaris Miyazaki F as studied by 1H NMR

All of the C2 proton signals of the coordinated histidine residues in the 1H NMR spectrum of cytochrome c3 from D. vulgaris Miyazaki F were assigned by specific deuteration. They appeared at extremely high fields and scattered in a wide range from -4 to -22 ppm. This clearly shows that the chemical properties of the imidazole groups are quite different from one another. The extremely high-field shift of the C2 signal indicates that some of them must carry the imidazolate-like nature to some extent. This might be responsible for the extremely low redox potentials of the four hemes. On changing temperature, most of them showed Curie-type change. All of the C2 signals showed a small p2H dependence in the range of p2H 4.8-10.0.

Mössbauer characterization of the tetraheme cytochrome c3 from Desulfovibrio baculatus (DSM 1743). Spectral deconvolution of the heme components

Mössbauer spectroscopy was used to study the tetraheme cytochrome c3 from Desulfovibrio baculatus (DSM 1743). Samples with different degrees of reduction were prepared using a redoxtitration technique. In the reduced cytochrome c3, all four hemes are reduced and exhibit diamagnetic Mössbauer spectra typical for low-spin ferrous hemes (S = 0). In the oxidized protein, the hemes are low-spin ferric (S = 1/2) and exhibit overlapping magnetic Mössbauer spectra. A method of differential spectroscopy was applied to deconvolute the four overlapping heme spectra and a crystal-field model was used for data analysis. Characteristic Mössbauer spectral components for each heme group are obtained. Hyperfine and crystal-field parameters for all four hemes are determined from these deconvoluted spectra.

A thermodynamic model for the cooperative functional properties of the tetraheme cytochrome c3 from Desulfovibrio gigas

A thermodynamic model is presented to describe the redox behaviour of the tetraheme cytochrome c3 from Desulfovibrio gigas. This molecule displays different intrinsic redox potentials for the four hemes and during the redox titration process, interactions among different hemes occur, thus altering the values of redox potentials according to which of the hemes are oxidized [Santos, H., Moura, J.J.G., Moura, I., LeGall, J. & Xavier, A.V. (1984) Eur. J. Biochem. 141, 283-296]. This complex cooperative behaviour [Xavier, A.V. (1986) J. Inorg. Biochem. 28, 239-243] has been analyzed here using an I2H4-interaction network [Cornish-Bowden, A. & Koshland, D.E. Jr (1970) J. Biol. Chem. 245, 6241-6250] coupled to a proton-linked equilibrium between two tertiary structures. Such a formalism, which requires a reduced number of parameters, is able to fully account quantitatively for the pH dependence of the NMR redox-titration curves. The 'redox-Bohr' effect is discussed in terms of the available structure and thermodynamic data and a functional mechanism is proposed.