Molecular dynamics simulation of cytochrome c3: studying the reduction processes using free energy calculations

Graphite/gold nanoparticles electrode for direct protein attachment: characterization and gas sensing application

In this work, graphite/gold nanoparticles (G/AuNPs) were synthesized through a facile chemical method, and its potential application for direct protein attachment for electrochemical detection of carbon monoxide (CO) was investigated. The preparation of G/AuNPs electrodes was optimized by synthesizing the nanoparticles in different concentration of HAuCl₄.3H₂O at various temperatures. The G/AuNPs electrode was subsequently modified by four types of mercaptopropionic acid, including 1-mercaptopropionic, 3-mercaptopropionic, 6-mercaptopropionic, and 11-mercaptopropionic acid, to achieve the best structure for protein attachment. Visible absorption and electrochemical studies showed that 3-mercaptopropionic acid possesses the best performance regarding the electrical conductivity between electrode and protein redox center. The cyclic voltammetry results revealed that the modified electrode has an appropriate performance for CO detection at very low concentrations while keeping a linear response. The limit of detection for the modified electrode was calculated to be about 0.2 ppb. Finally, the interactions of cytochrome C and carbon monoxides were simulated using molecular dynamics (MD), and the effect of protein conformation changes on the electrochemical signal was thoroughly examined. The simulation results suggested that the proposed electrochemical sensor has an acceptable performance for the detection of CO due to less fluctuation of amino acids near the protein chain in the presence of CO molecules.

A brief survey of the "cytochromome"

Multihaem cytochromes c are widespread in nature where they perform numerous roles in diverse anaerobic metabolic pathways. This is achieved in two ways: multihaem cytochromes c display a remarkable diversity of ways to organize multiple hemes within the protein frame; and the hemes possess an intrinsic reactive versatility derived from diverse spin, redox and coordination states. Here we provide a brief survey of multihaem cytochromes c that have been characterized in the context of their metabolic role. The contribution of multihaem cytochromes c to dissimilatory pathways handling metallic minerals, nitrogen compounds, sulfur compounds, organic compounds and phototrophism are described. This aims to set the stage for the further exploration of the vast unknown "cytochromome" that can be anticipated from genomic databases.

Possible dynamically gated conductance along heme wires in bacterial multiheme cytochromes

The staggered cross decaheme configuration of electron transfer cofactors in the outer-membrane cytochrome MtrF serves as a prototype for conformationally gated multiheme electron transport. Derived from the bacterium Shewanella oneidensis, the staggered cross configuration reveals intersecting c-type octaheme and tetraheme "wires" containing thermodynamic "hills" and "valleys" (Proc. Natl. Acad. Sci. U. S. A. 2014, 11, 611-616), suggesting that the protein structure may include a dynamical mechanism for conductance and pathway switching depending on enzymatic functional need. Here, we applied classical molecular and statistical mechanics calculations of large-amplitude protein dynamics in MtrF, to address its potential to modulate pathway conductance, including assessment of the effect of the total charge state. Explicit solvent molecular dynamics simulations of fully oxidized and fully reduced MtrF showed that the slowest mode of collective decaheme motion is 90% similar between the oxidized and reduced states and consists primarily of interheme separation with minor rotational contributions. The frequency of this motion is 1.7 × 10(7) s(-1), both for fully oxidized and fully reduced MtrF, slower than the downhill electron transfer rates between stacked heme pairs at the octaheme termini and faster than the electron transfer rates between parallel hemes in the tetraheme chain. This implies that MtrF uses slow conformational fluctuations to modulate electron flow along the octaheme pathway, apparently for the purpose of increasing the residence time of electrons on lowest potential hemes 4 and 9. This apparent gating mechanism should increase the success rate of electron transfer from MtrF to low potential environmental acceptors via these two solvent-exposed hemes.

Charge parametrization of the DvH-c3 heme group: validation using constant-(pH,E) molecular dynamics simulations

We studied the effect of using different heme group charge parametrization methods and schemes (Merz-Kollman, CHELPG, and single- and multiconformational RESP) on the quality of the results produced by the constant-(pH,E) MD method, applied to the redox titration of Desulfovibrio vulgaris Hildenborough cytochrome c(3). These new and more accurate charge sets enabled us to overcome the previously reported dependence of the method's performance on the dielectric constant, ε, assigned to the protein region. In particular, we found a systematic, clear shift of the E(mod) toward more negative values than those previously reported, in agreement with an electrostatics based reasoning. The simulations showed strong coupling between protonating/redox sites. We were also able to capture significant direct and, especially, indirect interactions between hemes, such as those mediated by histidine 67. Our results highlight the importance of having a good quantum description of the system before deriving atomic partial charges for classic force fields.

Free energies for biological electron transfer from QM/MM calculation: method, application and critical assessment

Computer simulations of biological electron transfer reactions are reviewed with a focus on the calculation of reaction free energy (driving force) and reorganization free energy. Then a mixed quantum mechanical/molecular mechanical (QM/MM) approach is described which is designed for computation of these quantities for pure electron transfer reactions with large donor-acceptor separation distances. The method is applied to intra-protein electron transfer in Ru(bpy)(2)(im)His33 cytochrome c and the results compared to experimental data. Several modeling aspects which are important for successful calculation of free energies with QM/MM are discussed in detail.

Cooperative interaction of high-potential hemes in the cytochrome subunit of the photosynthetic reaction center of bacterium Ectothiorhodospira shaposhnikovii

Cooperative interaction of the high-potential hemes (C(h)) in the cytochrome subunit of the photosynthesizing bacterium Ectothiorhodospira shaposhnikovii was studied by comparing redox titration curves of the hemes under the conditions of pulse photoactivation inducing single turnover of electron-transport chain and steady-state photoactivation, as well as by analysis of the kinetics of laser-induced oxidation of cytochromes by reaction center (RC). A mathematical model of the processes of electron transfer in cytochrome-containing RC was considered. Theoretical analysis revealed that the reduction of one heme C(h) facilitated the reduction of the other heme, which was equivalent to a 60 mV positive shift of the midpoint potential. In addition, reduction of the second heme C(h) caused a three- to four-fold acceleration of the electron transfer from the cytochrome subunit to RC.

Proton thrusters: overview of the structural and functional features of soluble tetrahaem cytochromes c 3

Tetrahaem cytochromes c (3) from sulfate-reducing bacteria have revealed exquisite complexity in their ligand binding properties and they couple the cooperative binding of two electrons with the binding of protons. In this review, the molecular mechanisms for these cooperative effects are described, and the functional consequences of these cooperativities are discussed in the context of the general mechanisms of biological energy transduction and the specific physiological metabolism of Desulfovibrio.

Electrostatic environment of hemes in proteins: pK(a)s of hydroxyl ligands

The pK(a)s of ferric aquo-heme and aquo-heme electrochemical midpoints (E(m)s) at pH 7 in sperm whale myoglobin, Aplysia myoblogin, hemoglobin I, heme oxygenase 1, horseradish peroxidase and cytochrome c oxidase were calculated with Multi-Conformation Continuum Electrostatics (MCCE). The pK(a)s span 3.3 pH units from 7.6 in heme oxygenase 1 to 10.9 in peroxidase, and the E(m)s range from -250 mV in peroxidase to 125 mV in Aplysia myoglobin. Proteins with higher in situ ferric aquo-heme pK(a)s tend to have lower E(m)s. Both changes arise from the protein stabilizing a positively charged heme. However, compared with values in solution, the protein shifts the aquo-heme E(m)s more than the pK(a)s. Thus, the protein has a larger effective dielectric constant for the protonation reaction, showing that electron and proton transfers are coupled to different conformational changes that are captured in the MCCE analysis. The calculations reveal a breakdown in the classical continuum electrostatic analysis of pairwise interactions. Comparisons with DFT calculations show that Coulomb's law overestimates the large unfavorable interactions between the ferric water-heme and positively charged groups facing the heme plane by as much as 60%. If interactions with Cu(B) in cytochrome c oxidase and Arg 38 in horseradish peroxidase are not corrected, the pK(a) calculations are in error by as much as 6 pH units. With DFT corrected interactions calculated pK(a)s and E(m)s differ from measured values by less than 1 pH unit or 35 mV, respectively. The in situ aquo-heme pK(a) is important for the function of cytochrome c oxidase since it helps to control the stoichiometry of proton uptake coupled to electron transfer [Song, Michonova-Alexova, and Gunner (2006) Biochemistry 45, 7959-7975].

Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations

A protein structure should provide the information needed to understand its observed properties. Significant progress has been made in developing accurate calculations of acid/base and oxidation/reduction reactions in proteins. Current methods and their strengths and weaknesses are discussed. The distribution and calculated ionization states in a survey of proteins is described, showing that a significant minority of acidic and basic residues are buried in the protein and that most of these remain ionized. The electrochemistry of heme and quinones are considered. Proton transfers in bacteriorhodopsin and coupled electron and proton transfers in photosynthetic reaction centers, 5-coordinate heme binding proteins and cytochrome c oxidase are highlighted as systems where calculations have provided insight into the reaction mechanism.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.

Reorganization and conformational changes in the reduction of tetraheme cytochromes

Molecular dynamics simulation (MD) constitutes an alternative to time-consuming experiments for studying conformational changes. We apply MD on a redox system where experimental information exists for the fully oxidized and fully reduced states: tetraheme cytochrome c3. Instead of doing one simulation for each state, we apply 10 4-ns replicas for both states, which provides robust statistics to characterize the redox changes. Besides these long simulations, we perform 120 short ones (50 ps), where an equilibrated oxidized state is perturbed to a reduced state. This allows the application of a nonequilibrium method, the subtraction technique, which makes it possible to characterize the different timescales of conformational changes. Reduction induces conformational changes in the N-terminus and on the loops spanning residues 36-42 and 88-93, which correlate very well with experiments, demonstrating the applicability of this methodology. We also analyze the effect of reduction on hydrogen bonds, solvent accessible surface and bound water, the changes being found to involve the hemes and propionate groups. Redox-induced protonation is also investigated, by protonating the propionates D from hemes I and IV. Although this change in the former does not have major conformational consequences, it induces in the latter conformational changes beyond the ones obtained with reduction.

Modeling electron transfer thermodynamics in protein complexes: interaction between two cytochromes c(3)

Redox protein complexes between type I and type II tetraheme cytochromes c(3) from Desulfovibrio vulgaris Hildenborough are here analyzed using theoretical methodologies. Various complexes were generated using rigid-body docking techniques, and the two lowest energy complexes (1 and 2) were relaxed using molecular dynamics simulations with explicit solvent and subjected to further characterization. Complex 1 corresponds to an interaction between hemes I from both cytochromes c(3). Complex 2 corresponds to an interaction between the heme IV from type I and the heme I from type II cytochrome c(3). Binding free energy calculations using molecular mechanics, Poisson-Boltzmann, and surface accessibility methods show that complex 2 is more stable than complex 1. Thermodynamic calculations on complex 2 show that complex formation induces changes in the reduction potential of both cytochromes c(3), but the changes are larger in the type I cytochrome c(3) (the largest one occurring on heme IV, of approximately 80 mV). These changes are sufficient to invert the global titration curves of both cytochromes, generating directionally in electron transfer from type I to type II cytochrome c(3), a phenomenon of obvious thermodynamic origin and consequences, but also with kinetic implications. The existence of processes like this occurring at complex formation may constitute a natural design of efficient redox chains.

Molecular basis for redox-Bohr and cooperative effects in cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774: Crystallographic and modeling studies of oxidized and reduced high-resolution structures at pH 7.6

The tetraheme cytochrome c3 is a small metalloprotein with ca. 13,000 Da found in sulfate-reducing bacteria, which is believed to act as a partner of hydrogenase. The three-dimensional structure of the oxidized and reduced forms of cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774 at pH 7.6 were determined using high-resolution X-ray crystallography and were compared with the previously determined oxidized form at pH 4.0. Theoretical calculations were performed with both structures, using continuum electrostatic calculations and Monte Carlo sampling of protonation and redox states, in order to understand the molecular basis of the redox-Bohr and cooperativity effects related to the coupled transfer of electrons and protons. We were able to identify groups that showed redox-linked conformational changes. In particular, Glu61, His76, and propionate D of heme II showed important contributions to the redox-cooperativity, whereas His76, propionate A of heme I, and propionate D of heme IV were the key residues for the redox-Bohr effect. Upon reduction, an important movement of the backbone region surrounding hemes I and II was also identified, that, together with a few redox-linked conformational changes in side-chain residues, results in a significant decrease in the solvent accessibility of hemes I and II.

Redox-Bohr and other cooperativity effects in the nine-heme cytochrome C from Desulfovibrio desulfuricans ATCC 27774: crystallographic and modeling studies

The nine-heme cytochrome c is a monomeric multiheme cytochrome found in Desulfovibrio desulfuricans ATCC 27774. The polypeptide chain comprises 296 residues and wraps around nine hemes of type c. It is believed to take part in the periplasmic assembly of proteins involved in the mechanism of hydrogen cycling, receiving electrons from the tetraheme cytochrome c3. With the purpose of understanding the molecular basis of electron transfer processes in this cytochrome, we have determined the crystal structures of its oxidized and reduced forms at pH 7.5 and performed theoretical calculations of the binding equilibrium of protons and electrons in these structures. This integrated study allowed us to observe that the reduction process induced relevant conformational changes in several residues, as well as protonation changes in some protonatable residues. In particular, the surroundings of hemes I and IV constitute two areas of special interest. In addition, we were able to ascertain the groups involved in the redox-Bohr effect present in this cytochrome and the conformational changes that may underlie the redox-cooperativity effects on different hemes. Furthermore, the thermodynamic simulations provide evidence that the N- and C-terminal domains function in an independent manner, with the hemes belonging to the N-terminal domain showing, in general, a lower redox potential than those found in the C-terminal domain. In this way, electrons captured by the N-terminal domain could easily flow to the C-terminal domain, allowing the former to capture more electrons. A notable exception is heme IX, which has low redox potential and could serve as the exit path for electrons toward other proteins in the electron transfer pathway.

How cytochromes with different folds control heme redox potentials

The electrochemical midpoint potentials (E(m)'s) of 13 cytochromes, in globin (c, c(2), c(551), c(553)), four-helix bundle (c', b(562)), alpha beta roll (b(5)), and beta sandwich (f) motifs, with E(m)'s spanning 450 mV were calculated with multiconformation continuum electrostatics (MCCE). MCCE calculates changes in oxidation free energy when a heme-axial ligand complex is moved from water into protein. Calculated and experimental E(m)'s are in good agreement for cytochromes with His-Met and bis-His ligated hemes, where microperoxidases provide reference E(m)'s. In all cytochromes, E(m)'s are raised by 130-260 mV relative to solvated hemes by the loss of reaction field (solvation) energy. However, there is no correlation between E(m) and heme surface exposure. Backbone amide dipoles in loops or helix termini near the axial ligands raise E(m)'s, but amides in helix bundles contribute little. Heme propionates lower E(m)'s. If the propionic acids are partially protonated in the reduced cytochrome, protons are released on heme oxidation, contributing to the pH dependence of the E(m). In all cytochromes studied except b(5)'s and low potential globins, buried side chains raise E(m)'s. MCCE samples ionizable group protonation states, heme redox states, and side chain rotamers simultaneously. Globins show the largest structural changes on heme oxidation and four-helix bundles the least. Given the calculated protein-induced E(m) shift and measured cytochrome E(m) the five-coordinate, His heme in c' is predicted to have a solution E(m) between that of isolated bis-His and His-Met hemes, while the reference E(m) for His-Ntr ligands in cytochrome f should be near that of His-Met hemes.

Molecular dynamics study of Desulfovibrio africanus cytochrome c3 in oxidized and reduced forms

A 5-ns molecular dynamics study of a tetraheme cytochrome in fully oxidized and reduced forms was performed using the CHARMM molecular modeling program, with explicit water molecules, Langevin dynamics thermalization, Particle Mesh Ewald long-range electrostatics, and quantum mechanical determination of heme partial charges. The simulations used, as starting points, crystallographic structures of the oxidized and reduced forms of the acidic cytochrome c(3) from Desulfovibrio africanus obtained at pH 5.6. In this paper we also report structures for the two forms obtained at pH 8. In contrast to previous cytochrome c(3) dynamics simulations, our model is stable. The simulation structures agree reasonably well with the crystallographic ones, but our models show higher flexibility and the water molecules are more labile. We have compared in detail the differences between the simulated and experimental structures of the two redox states and observe that the hydration structure is highly dependent on the redox state. We have also analyzed the interaction energy terms between the hemes, the protein residues, and water. The direct electrostatic interaction between hemes is weak and nearly insensitive to the redox state, but the remaining terms are large and contribute in a complex way to the overall potential energy differences that we see between the redox states.

Solution structure and dynamics of the functional domain of Paracoccus denitrificans cytochrome c(552) in both redox states

A soluble and fully functional 10.5 kDa fragment of the 18.2 kDa membrane-bound cytochrome c(552) from Paracoccus denitrificans has been heterologously expressed and (13)C/(15)N-labeled to study the structural features of this protein in both redox states. Well-resolved solution structures of both the reduced and oxidized states have been determined using high-resolution heteronuclear NMR. The overall protein topology consists of two long terminal helices and three shorter helices surrounding the heme moiety. No significant redox-induced structural differences have been observed. (15)N relaxation rates and heteronuclear NOE values were determined at 500 and 600 MHz. Several residues located around the heme moiety display increased backbone mobility in both oxidation states, while helices I, III, and V as well as the two concatenated beta-turns between Leu30 and Arg36 apparently form a less flexible domain within the protein structure. Major redox-state-dependent differences of the internal backbone mobility on the picosecond-nanosecond time scale were not evident. Hydrogen exchange experiments demonstrated that the slow-exchanging amide proton resonances mainly belong to the helices and beta-turns, corresponding to the regions with high order parameters in the dynamics data. Despite this correlation, the backbone amide protons of the oxidized cytochrome c(552) exchange considerably faster with the solvent compared to the reduced protein. Using both differential scanning calorimetry as well as temperature-dependent NMR spectroscopy, a significant difference in the thermostabilities of the two redox states has been observed, with transition temperatures of 349.9 K (76.8 degrees C) for reduced and 307.5 K (34.4 degrees C) for oxidized cytochrome c(552). These results suggest a clearly distinct backbone stability between the two oxidation states.

Free energy calculations and molecular dynamics simulations of wild-type and variants of the DNA-EcoRI complex

Molecular dynamics simulations and free energy calculations of the wild-type EcoRI-DNA complex and several variants have been performed in aqueous solvent. In general, he theoretical estimations of the free energy differences (DeltaDeltaA) qualitatively agree well with the corresponding experimental data. The modifications which were experimentally found unfavorable compared to the wild-type complex were also found to be so in theoretical estimates. The mutant where the amino group of the base Ade(6) was replaced by a hydrogen atom eliminating one H-bond between the DNA and the protein, was experimentally found to be more stable than the wild-type complex. It was speculated that the modification also caused a structural relaxation in the DNA making DeltaDeltaA favorable. Our theoretical estimate yields a positive DeltaDeltaA in this case, but the difference is small, and no significant local structural relaxation was observed. The major H-bonds between the DNA and the protein in the wild-type complex are found to be maintained in the different mutants although the specific and non-specific interaction energies between the interacting the DNA bases and the protein residues are different in different mutants. The interaction pattern of the other nearby nucleotides are significantly influenced by each modification. Thus, the alteration of the non-specific interactions may also play an indirect role in determining the specificity of the complex. The interaction of the Gua(4) of the DNA with the protein is found to be most sensitive to any alteration in the recognition site. Because Gua(4) is the nucleotide closest to the scissile bond, this extra sensitivity seems to play an important role in altering the functional activity of the complex.

Crystal structure of the oxidised and reduced acidic cytochrome c3from Desulfovibrio africanus

Unique among sulphate-reducing bacteria, Desulfovibrio africanus has two periplasmic tetraheme cytochromes c3, one with an acidic isoelectric point which exhibits an unusually low reactivity towards hydrogenase, and another with a basic isoelectric point which shows the usual cytochrome c3reactivity. The crystal structure of the oxidised acidic cytochrome c3of Desulfovibrio africanus (Dva.a) was solved by the multiple anomalous diffraction (MAD) method and refined to 1.6 A resolution. Its structure clearly belongs to the same family as the other known cytochromes c3, but with weak parentage with those of the Desulfovibrio genus and slightly closer to the cytochromes c3of Desulfomicrobium norvegicum. In Dva.a, one edge of heme I is completely exposed to the solvent and surrounded by a negatively charged protein surface. Heme I thus seems to play an important role in electron exchange, in addition to heme III or heme IV which are the electron exchange ports in the other cytochromes c3. The function of Dva.a and the nature of its redox partners in the cell are thus very likely different. By alignment of the seven known 3D structures including Dva.a, it is shown that the structure which is most conserved in all cytochromes c3is the four-heme cluster itself. There is no conserved continuous protein structure which could explain the remarkable invariance of the four-heme cluster. On the contrary, the proximity of the heme edges is such that they interact directly by hydrophobic and van der Waals contacts. This direct interaction, which always involves a pyrrole CA-CB side-chain and its bound protein cysteine Sgammaatom, is probably the main origin of the four-heme cluster stability. The same kind of interaction is found in the chaining of the hemes in other multihemic redox proteins.The crystal structure of reduced Dva. a was solved at 1.9 A resolution. The comparison of the oxidised and reduced structures reveals changes in the positions of water molecules and polar residues which probably result from changes in the protonation state of amino acids and heme propionates. Water molecules are found closer to the hemes and to the iron atoms in the reduced than in the oxidised state. A global movement of a chain fragment in the vicinity of hemes III and IV is observed which result very likely from the electrostatic reorganization of the polypeptide chain induced by reduction.

Simulation of electron-proton coupling with a Monte Carlo method: application to cytochrome c3 using continuum electrostatics

A new method is presented for simulating the simultaneous binding equilibrium of electrons and protons on protein molecules, which makes it possible to study the full equilibrium thermodynamics of redox and protonation processes, including electron-proton coupling. The simulations using this method reflect directly the pH and electrostatic potential of the environment, thus providing a much closer and realistic connection with experimental parameters than do usual methods. By ignoring the full binding equilibrium, calculations usually overlook the twofold effect that binding fluctuations have on the behavior of redox proteins: first, they affect the energy of the system by creating partially occupied sites; second, they affect its entropy by introducing an additional empty/occupied site disorder (here named occupational entropy). The proposed method is applied to cytochrome c3 of Desulfovibrio vulgaris Hildenborough to study its redox properties and electron-proton coupling (redox-Bohr effect), using a continuum electrostatic method based on the linear Poisson-Boltzmann equation. Unlike previous studies using other methods, the full reduction order of the four hemes at physiological pH is successfully predicted. The sites more strongly involved in the redox-Bohr effect are identified by analysis of their titration curves/surfaces and the shifts of their midpoint redox potentials and pKa values. Site-site couplings are analyzed using statistical correlations, a method much more realistic than the usual analysis based on direct interactions. The site found to be more strongly involved in the redox-Bohr effect is propionate D of heme I, in agreement with previous studies; other likely candidates are His67, the N-terminus, and propionate D of heme IV. Even though the present study is limited to equilibrium conditions, the possible role of binding fluctuations in the concerted transfer of protons and electrons under nonequilibrium conditions is also discussed. The occupational entropy contributions to midpoint redox potentials and pKa values are computed and shown to be significant.

Functional and mechanistic studies of cytochrome c3 from Desulfovibrio gigas: thermodynamics of a "proton thruster"

Nuclear magnetic resonance and visible spectroscopies were used to determine the thermodynamic parameters of the four hemes in cytochrome c3 from Desulfovibrio gigas at 298 and 277 K and to investigate the mechanism of electron/proton energy transduction. Data obtained in the pH range from 5 to 9 were analyzed according to a model in which the hemes interact with each other (redox cooperativities) and with an ionizable center (redox-Bohr cooperativities). The results obtained at the two temperatures allow the deconvolution of the entropic contribution to the free energy of the four hemes, to the acid-base equilibrium of the ionizable center, and to the network of cooperativities among the five centers. The redox potentials of the hemes are modulated by the enthalpic contribution to the free energy, and evidence for the participation of the propionates of heme I in the redox-Bohr effect is presented. The network of interactions between the centers in this protein facilitates the concerted transfer of electrons and protons, in agreement with the "proton thruster" mechanism proposed for electronic to protonic energy transduction by cytochromes c3.

Solution structure of Desulfovibrio vulgaris (Hildenborough) ferrocytochrome c3: structural basis for functional cooperativity

Desulfovibrio vulgaris cytochrome c3 is a 14 kDa tetrahaem cytochrome that plays a central role in energy transduction. The three-dimensional structure of the ferrocytochrome at pH 8.5 was solved through two-dimensional 1H-NMR. The structures were calculated using a large amount of experimental information, which includes upper and lower distance limits as well as dihedral angle restraints. The analysis allows for fast-flipping aromatic residues and flexibility in the haem plane. The structure was determined using 2289 upper and 2390 lower distance limits, 63 restricted ranges for the phi torsion angle, 88 stereospecific assignments out of the 118 stereopairs with non-degenerate chemical shifts (74.6%), and 115 out of the 184 nuclear Overhauser effects to fast-flipping aromatic residues (62.5%), which were pseudo-stereospecifically assigned to one or the other side of the ring. The calculated NMR structures are very well defined, with an average root-mean-square deviation value relative to the mean coordinates of 0.35 A for the backbone atoms and 0.70 A for all heavy-atoms. Comparison of the NMR structures of the ferrocytochrome at pH 8.5 with the available X-ray structure of the ferricytochrome at pH 5.5 reveals that the general fold of the molecule is very similar, but that there are some distinct differences. Calculation of ring current shifts for the residues with significantly different conformations confirms that the NMR structures represent better its solution structure in the reduced form. Some of the localised differences, such as a reorientation of Thr24, are thought to be state-dependent changes that involve alterations in hydrogen bond networks. An important rearrangement in the vicinity of the propionate groups of haem I and involving the covalent linkage of haem II suggests that this is the critical region for the functional cooperativities of this protein.