Molecular structure of cytochrome c2 isolated from Rhodobacter capsulatus determined at 2.5 A resolution

Carotenoids and Photosynthesis

Carotenoids are ubiquitous and essential pigments in photosynthesis. They absorb in the blue-green region of the solar spectrum and transfer the absorbed energy to (bacterio-)chlorophylls, and so expand the wavelength range of light that is able to drive photosynthesis. This is an example of singlet-singlet energy transfer, and so carotenoids serve to enhance the overall efficiency of photosynthetic light reactions. Carotenoids also act to protect photosynthetic organisms from the harmful effects of excess exposure to light. Triplet-triplet energy transfer from chlorophylls to carotenoids plays a key role in this photoprotective reaction. In the light-harvesting pigment-protein complexes from purple photosynthetic bacteria and chlorophytes, carotenoids have an additional role of structural stabilization of those complexes. In this article we review what is currently known about how carotenoids discharge these functions. The molecular architecture of photosynthetic systems will be outlined first to provide a basis from which to describe carotenoid photochemistry, which underlies most of their important functions in photosynthesis.

Microscopic insights into the NMR relaxation-based protein conformational entropy meter

Conformational entropy is a potentially important thermodynamic parameter contributing to protein function. Quantitative measures of conformational entropy are necessary for an understanding of its role but have been difficult to obtain. An empirical method that utilizes changes in conformational dynamics as a proxy for changes in conformational entropy has recently been introduced. Here we probe the microscopic origins of the link between conformational dynamics and conformational entropy using molecular dynamics simulations. Simulation of seven proteins gave an excellent correlation with measures of side-chain motion derived from NMR relaxation. The simulations show that the motion of methyl-bearing side chains are sufficiently coupled to that of other side chains to serve as excellent reporters of the overall side-chain conformational entropy. These results tend to validate the use of experimentally accessible measures of methyl motion--the NMR-derived generalized order parameters--as a proxy from which to derive changes in protein conformational entropy.

Structural characterization of nitrosomonas europaea cytochrome c-552 variants with marked differences in electronic structure

Nitrosomonas europaea cytochrome c-552 (Ne c-552) variants with the same His/Met axial ligand set but with different EPR spectra have been characterized structurally, to aid understanding of how molecular structure determines heme electronic structure. Visible light absorption, Raman, and resonance Raman spectroscopy of the protein crystals was performed along with structure determination. The structures solved are those of Ne c-552, which displays a "HALS" (or highly anisotropic low-spin) EPR spectrum, and of the deletion mutant Ne N64Δ, which has a rhombic EPR spectrum. Two X-ray crystal structures of wild-type Ne c-552 are reported; one is of the protein isolated from N. europaea cells (Ne c-552n, 2.35 Å resolution), and the other is of recombinant protein expressed in Escherichia coli (Ne c-552r, 1.63 Å resolution). Ne N64Δ crystallized in two different space groups, and two structures are reported [monoclinic (2.1 Å resolution) and hexagonal (2.3 Å resolution)]. Comparison of the structures of the wild-type and mutant proteins reveals that heme ruffling is increased in the mutant; increased ruffling is predicted to yield a more rhombic EPR spectrum. The 2.35 Å Ne c-552n structure shows 18 molecules in the asymmetric unit; analysis of the structure is consistent with population of more than one axial Met configuration, as seen previously by NMR. Finally, the mutation was shown to yield a more hydrophobic heme pocket and to expel water molecules from near the axial Met. These structures reveal that heme pocket residue 64 plays multiple roles in regulating the axial ligand orientation and the interaction of water with the heme. These results support the hypothesis that more ruffled hemes lead to more rhombic EPR signals in cytochromes c with His/Met axial ligation.

Inversion of the stereochemistry around the sulfur atom of the axial methionine side chain through alteration of amino acid side chain packing in Hydrogenobacter thermophilus cytochrome C552 and its functional consequences

In cytochrome c, the coordination of the axial Met Sδ atom to the heme Fe atom occurs in one of two distinctly different stereochemical manners, i.e., R and S configurations, depending upon which of the two lone pairs of the Sδ atom is involved in the bond; hence, the Fe-coordinated Sδ atom becomes a chiral center. In this study, we demonstrated that an alteration of amino acid side chain packing induced by the mutation of a single amino acid residue, i.e., the A73V mutation, in Hydrogenobacter thermophilus cytochrome c552 (HT) forces the inversion of the stereochemistry around the Sδ atom from the R configuration [Travaglini-Allocatelli, C., et al. (2005) J. Biol. Chem. 280, 25729-25734] to the S configuration. Functional comparison between the wild-type HT and the A73V mutant possessing the R and S configurations as to the stereochemistry around the Sδ atom, respectively, demonstrated that the redox potential (Em) of the mutant at pH 6.00 and 25 °C exhibited a positive shift of ∼20 mV relative to that of the wild-type HT, i.e., 245 mV, in an entropic manner. Because these two proteins have similar enthalpically stabilizing interactions, the difference in the entropic contribution to the Em value between them is likely to be due to the effect of the conformational alteration of the axial Met side chain associated with the inversion of the stereochemistry around the Sδ atom due to the effect of mutation on the internal mobility of the loop bearing the axial Met. Thus, the present study demonstrated that the internal mobility of the loop bearing the axial Met, relevant to entropic control of the redox function of the protein, is affected quite sensitively by the contextual stereochemical packing of amino acid side chains in the proximity of the axial Met.

The heme chaperone ApoCcmE forms a ternary complex with CcmI and apocytochrome c

Cytochrome c maturation (Ccm) is a post-translational process that occurs after translocation of apocytochromes c to the positive (p) side of energy-transducing membranes. Ccm is responsible for the formation of covalent bonds between the thiol groups of two cysteines residues at the heme-binding sites of the apocytochromes and the vinyl groups of heme b (protoporphyrin IX-Fe). Among the proteins (CcmABCDEFGHI and CcdA) required for this process, CcmABCD are involved in loading heme b to apoCcmE. The holoCcmE thus formed provides heme b to the apocytochromes. Catalysis of the thioether bonds between the apocytochromes c and heme b is mediated by the heme ligation core complex, which in Rhodobacter capsulatus contains at least the CcmF, CcmH, and CcmI components. In this work we show that the heme chaperone apoCcmE binds to the apocytochrome c and the apocytochrome c chaperone CcmI to yield stable binary and ternary complexes in the absence of heme in vitro. We found that during these protein-protein interactions, apoCcmE favors the presence of a disulfide bond at the apocytochrome c heme-binding site. We also establish using detergent-dispersed membranes that apoCcmE interacts directly with CcmI and CcmH of the heme ligation core complex CcmFHI. Implications of these findings are discussed with respect to heme transfer from CcmE to the apocytochromes c during heme ligation assisted by the core complex CcmFHI.

CcmI subunit of CcmFHI heme ligation complex functions as an apocytochrome c chaperone during c-type cytochrome maturation

Cytochrome c maturation (Ccm) is a sophisticated post-translational process. It occurs after translocation of apocytochromes c to the p side of energy transducing membranes and forms stereo-specific thioether bonds between the vinyl groups of heme b (protoporphyrin IX-Fe) and the thiol groups of cysteines at their conserved heme binding sites. In many organisms this process involves up to 10 (CcmABCDEFGHI and CcdA) membrane proteins. One of these proteins is CcmI, which has an N-terminal membrane-embedded domain with two transmembrane helices and a large C-terminal periplasmic domain with protein-protein interaction motifs. Together with CcmF and CcmH, CcmI forms a multisubunit heme ligation complex. How the CcmFHI complex recognizes its apocytochrome c substrates remained unknown. In this study, using Rhodobacter capsulatus apocytochrome c(2) as a Ccm substrate, we demonstrate for the first time that CcmI binds apocytochrome c(2) but not holocytochrome c(2). Mainly the C-terminal portions of both CcmI and apocytochrome c(2) mediate this binding. Other physical interactions via the conserved structural elements in apocytochrome c(2), like the heme ligating cysteines or heme iron axial ligands, are less crucial. Furthermore, we show that the N-terminal domain of CcmI can also weakly bind apocytochrome c(2), but this interaction requires a free thiol group at apocytochrome c(2) heme binding site. We conclude that the CcmI subunit of the CcmFHI complex functions as an apocytochrome c chaperone during the Ccm process used by proteobacteria, archaea, mitochondria of plants and red algae.

Thiol redox requirements and substrate specificities of recombinant cytochrome c assembly systems II and III

The reconstitution of biosynthetic pathways from heterologous hosts can help define the minimal genetic requirements for pathway function and facilitate detailed mechanistic studies. Each of the three pathways for the assembly of cytochrome c in nature (called systems I, II, and III) has been shown to function recombinantly in Escherichia coli, covalently attaching heme to the cysteine residues of a CXXCH motif of a c-type cytochrome. However, recombinant systems I (CcmABCDEFGH) and II (CcsBA) function in the E. coli periplasm, while recombinant system III (CCHL) attaches heme to its cognate receptor in the cytoplasm of E. coli, which makes direct comparisons between the three systems difficult. Here we show that the human CCHL (with a secretion signal) attaches heme to the human cytochrome c (with a signal sequence) in the E. coli periplasm, which is bioenergetically (p-side) analogous to the mitochondrial intermembrane space. The human CCHL is specific for the human cytochrome c, whereas recombinant system II can attach heme to multiple non-cognate c-type cytochromes (possessing the CXXCH motif.) We also show that the recombinant periplasmic systems II and III use components of the natural E. coli periplasmic DsbC/DsbD thiol-reduction pathway. This article is part of a Special Issue entitled: Biogenesis/Assembly of Respiratory Enzyme Complexes.

Reorientation of cytochrome c2 upon interaction with oppositely charged macromolecules probed by SR EPR: implications for the role of dipole moment to facilitate collisions in proper configuration for electron transfer

The reaction of water-soluble cytochrome c (c(2)) with its physiological redox partners is facilitated by electrostatic attractions between the two protein surfaces. Using spin-labeled cytochrome c(2) from Rhodobacter capsulatus and pulse electron paramagnetic resonance (EPR) measurements we compared spatial orientation of cytochrome c(2) upon its binding to surfaces of opposite charge. We observed that cytochrome c(2) can use its negatively charged "back" side when exposed to interact with positively charged surfaces (DEAE resin) which is the opposite to the use of its positively charged "front" side in physiological interaction with negatively charged binding domain of cytochrome bc(1). The later orientation is also adopted upon non-physiological binding of cytochrome c(2) to negatively charged carboxymethyl cellulose resin. These results directly demonstrate how the electric dipolar nature of cytochrome c(2) influences its orientation in interactions with charged surfaces, which may facilitate collisions with other redox proteins in a proper orientation to support physiologically-competent electron transfer. Saturation recovery EPR provides an attractive tool for monitoring spatial orientation of proteins in their interaction with surfaces in liquid phase. It is particularly valuable for metalloproteins engaged in redox reactions as a means to monitor the geometry and dynamics of formation of protein complexes in measurements that are independent of electron transfer processes.

NMR and DFT investigation of heme ruffling: functional implications for cytochrome c

Out-of-plane (OOP) deformations of the heme cofactor are found in numerous heme-containing proteins and the type of deformation tends to be conserved within functionally related classes of heme proteins. We demonstrate correlations between the heme ruffling OOP deformation and the (13)C and (1)H nuclear magnetic resonance (NMR) hyperfine shifts of heme aided by density functional theory (DFT) calculations. The degree of ruffling in the heme cofactor of Hydrogenobacter thermophilus cytochrome c(552) has been modified by a single amino acid mutation in the second coordination sphere of the cofactor. The (13)C and (1)H resonances of the cofactor have been assigned using one- and two-dimensional NMR spectroscopy aided by selective (13)C-enrichment of the heme. DFT has been used to predict the NMR hyperfine shifts and electron paramagnetic resonance (EPR) g-tensor at several points along the ruffling deformation coordinate. The DFT-predicted NMR and EPR parameters agree with the experimental observations, confirming that an accurate theoretical model of the electronic structure and its response to ruffling has been established. As the degree of ruffling increases, the heme methyl (1)H resonances move upfield while the heme methyl and meso (13)C resonances move downfield. These changes are a consequence of altered overlap of the Fe 3d and porphyrin pi orbitals, which destabilizes all three occupied Fe 3d-based molecular orbitals and decreases the positive and negative spin density on the beta-pyrrole and meso carbons, respectively. Consequently, the heme ruffling deformation decreases the electronic coupling of the cofactor with external redox partners and lowers the reduction potential of heme.

Evidence from the structure and function of cytochromes c(2) that nonsulfur purple bacterial photosynthesis followed the evolution of oxygen respiration

Cytochromes c(2) are the nearest bacterial homologs of mitochondrial cytochrome c. The sequences of the known cytochromes c(2) can be placed in two subfamilies based upon insertions and deletions, one subfamily is most like mitochondrial cytochrome c (the small C2s, without significant insertions and deletions), and the other, designated large C2, shares 3- and 8-residue insertions as well as a single-residue deletion. C2s generally function between cytochrome bc(1) and cytochrome oxidase in respiration (ca 80 examples known to date) and between cytochrome bc(1) and the reaction center in nonsulfur purple bacterial photosynthesis (ca 21 examples). However, members of the large C2 subfamily are almost always involved in photosynthesis (12 of 14 examples). In addition, the gene for the large C2 (cycA) is associated with those for the photosynthetic reaction center (pufBALM). We hypothesize that the insertions in the large C2s, which were already functioning in photosynthesis, allowed them to replace the membrane-bound tetraheme cytochrome, PufC, that otherwise mediates between the small C2 or other redox proteins and photosynthetic reaction centers. Based upon our analysis, we propose that the involvement of C2 in nonsulfur purple bacterial photosynthesis was a metabolic feature subsequent to the evolution of oxygen respiration.

Demonstration of short-lived complexes of cytochrome c with cytochrome bc1 by EPR spectroscopy: implications for the mechanism of interprotein electron transfer

One of the steps of a common pathway for biological energy conversion involves electron transfer between cytochrome c and cytochrome bc1. To clarify the mechanism of this reaction, we examined the structural association of those two proteins using the electron transfer-independent electron paramagnetic resonance (EPR) techniques. Drawing on the differences in the continuous wave EPR spectra and saturation recoveries of spin-labeled bacterial and mitochondrial cytochromes c recorded in the absence and presence of bacterial cytochrome bc1, we have exposed a time scale of dynamic equilibrium between the bound and the free state of cytochrome c at various ionic strengths. Our data show a successive decrease of the bound cytochrome c fraction as the ionic strength increases, with a limit of approximately 120 mm NaCl above which essentially no bound cytochrome c can be detected by EPR. This limit does not apply to all of the interactions of cytochrome c with cytochrome bc1 because the cytochrome bc1 enzymatic activity remained high over a much wider range of ionic strengths. We concluded that EPR monitors just the tightly bound state of the association and that an averaged lifetime of this state decreases from over 100 micros at low ionic strength to less than 400 ns at an ionic strength above 120 mm. This suggests that at physiological ionic strength, the tightly bound complex on average lasts less than the time needed for a single electron exchange between hemes c and c1, indicating that productive electron transfer requires several collisions of the two molecules. This is consistent with an early idea of diffusion-coupled reactions that link the soluble electron carriers with the membranous complexes, which, we believe, provides a robust means of regulating electron flow through these complexes.

Insights into porphyrin chemistry provided by the microperoxidases, the haempeptides derived from cytochrome c

The water-soluble haem-containing peptides obtained by proteolytic digestion of cytochrome c, the microperoxidases, have been used to explore aspects of the chemistry of iron porphyrins, and as mimics for some reactions catalysed by the haemoproteins, including the reactions catalysed by the peroxidases and the cytochromes P450. The preparation of the microperoxidases, their physical and chemical properties including their electronic structure, the kinetics and thermodynamics of their reactions with ligands, electrochemical studies and examples of their uses as haemoproteins mimics, is reviewed.

Plasmon waveguide resonance spectroscopic evidence for differential binding of oxidized and reduced Rhodobacter capsulatus cytochrome c2 to the cytochrome bc1 complex mediated by the conformation of the Rieske iron-sulfur protein

The dissociation constants for the binding of Rhodobacter capsulatus cytochrome c2 and its K93P mutant to the cytochrome bc1 complex embedded in a phospholipid bilayer were measured by plasmon waveguide resonance spectroscopy in the presence and absence of the inhibitor stigmatellin. The reduced form of cytochrome c2 strongly binds to reduced cytochrome bc1 (Kd = 0.02 microM) but binds much more weakly to the oxidized form (Kd = 3.1 microM). In contrast, oxidized cytochrome c2 binds to oxidized cytochrome bc1 in a biphasic fashion with Kd values of 0.11 and 0.58 microM. Such a biphasic interaction is consistent with binding to two separate sites or conformations of oxidized cytochrome c2 and/or cytochrome bc1. However, in the presence of stigmatellin, we find that oxidized cytochrome c2 binds to oxidized cytochrome bc1 in a monophasic fashion with high affinity (Kd = 0.06 microM) and reduced cytochrome c2 binds less strongly (Kd = 0.11 microM) but approximately 30-fold more tightly than in the absence of stigmatellin. Structural studies with cytochrome bc1, with and without the inhibitor stigmatellin, have led to the proposal that the Rieske protein is mobile, moving between the cytochrome b and cytochrome c1 components during turnover. In one conformation, the Rieske protein binds near the heme of cytochrome c1, while the cytochrome c2 binding site is also near the cytochrome c1 heme but on the opposite side from the Rieske site, where cytochrome c2 cannot directly interact with Rieske. However, the inhibitor, stigmatellin, freezes the Rieske protein iron-sulfur cluster in a conformation proximal to cytochrome b and distal to cytochrome c1. We conclude from this that the dual conformation of the Rieske protein is primarily responsible for biphasic binding of oxidized cytochrome c2 to cytochrome c1. This optimizes turnover by maximizing binding of the substrate, oxidized cytochrome c2, when the iron-sulfur cluster is proximal to cytochrome b and minimizing binding of the product, reduced cytochrome c2, when it is proximal to cytochrome c1.

Local stability of Rhodobacter capsulatus cytochrome c2 probed by solution phase hydrogen/deuterium exchange and mass spectrometry

The hydrogen/deuterium exchange kinetics of Rhodobacter capsulatus cytochrome c2 have been determined using mass spectrometry. As expected, the relative domain stability was generally similar to that of the cytochrome c2 structural homolog, horse heart cytochrome c, but we were able to find evidence to support the presence of a second, small beta-sheet not found in the horse cytochrome, which stabilizes a structural region dominated by Omega loops. Importantly, we find that the so-called hinge region, comprised of 15 amino acids, which include the methionine sixth heme ligand (M96), is destabilized on oxidation, and this destabilization is propagated to a portion of the second Omega loop, most likely through perturbation of two hydrogen bonds that couple these two domains in the three dimensional structure. The mutation of a lysine at position 93 to proline amplifies the destabilization observed on oxidation of the wild-type cytochrome c2 and results in further destabilization observed in regions 52-60, 75-82, and 83-97. This suggests that hydrogen bond interactions involving two bound waters, the T94 hydroxyl, the front heme propionate and the Y75 hydroxyl, are significantly compromised upon mutation. In summary, these observations are consistent with the approximately 20-fold increase in the movement of the hinge away from the heme face in the oxidized cytochrome c2 as determined by ligand binding kinetics. Thus, H/D exchange kinetics can be used to identify relatively subtle structural features and at least in some cases facilitate the understanding of the structural basis of the dynamic properties of proteins.

The effect of replacing the axial methionine ligand with a lysine residue in cytochrome c-550 from Paracoccus versutus assessed by X-ray crystallography and unfolding

The structure of cytochrome c-550 from the nonphotosynthetic bacteria Paraccocus versutus has been solved by X-ray crystallography to 1.90 A resolution, and reveals a high structural homology to other bacterial cytochromes c(2). The effect of replacing the axial heme-iron methionine ligand with a lysine residue on protein structure and unfolding has been assessed using the M100K variant. From X-ray structures at 1.95 and 1.55 A resolution it became clear that the amino group of the lysine side chain coordinates to the heme-iron. Structural differences compared to the wild-type protein are confined to the lysine ligand loop connecting helices four and five. In the heme cavity an additional water molecule is found which participates in an H-bonding interaction with the lysine ligand. Under cryo-conditions extra electron density in the lysine ligand loop is revealed, leading to residues K97 to T101 being modeled with a double main-chain conformation. Upon unfolding, dissociation of the lysine ligand from the heme-iron is shown to be pH dependent, with NMR data consistent with the occurrence of a ligand exchange mechanism similar to that seen for the wild-type protein.

Simulation of Electron Transfer between Cytochrome c2 and the Bacterial Photosynthetic Reaction Center: Brownian Dynamics Analysis of the Native Proteins and Double Mutants

Electron transfer is essential for bacterial photosynthesis which converts light energy into chemical energy. This paper theoretically studies the interprotein electron transfer from cytochrome c(2) of Rhodobacter capsulatus to the photosynthetic reaction center of Rhodobacter sphaeroides in native and mutated systems. Brownian dynamics is used with an exponential distance-dependent electron-transfer rate model to compute bimolecular rate constants, which are consistent with experimental data when reasonable prefactors and decay constants are used. Interestingly, switching of the reaction mechanism from the diffusion-controlled limit in the native proteins to the activation-controlled limit in one of the mutants (DK(L261)/KE(C99)) was found. We also predict that the second-order rate for the native reaction center/cytochrome c(2) system will decrease with increasing ionic strength, a characteristic of electrostatically controlled docking.

Structural and functional studies on the tetraheme cytochrome subunit and its electron donor proteins: the possible docking mechanisms during the electron transfer reaction

The photosynthetic reaction centers (RCs) classified as the group II possess a peripheral cytochrome (Cyt) subunit, which serves as the electron mediator to the special-pair. In the cycle of the photosynthetic electron transfer reactions, the Cyt subunit accepts electrons from soluble electron carrier proteins, and re-reduces the photo-oxidized special-pair of the bacteriochlorophyll. Physiologically, high-potential cytochromes such as the cytochrome c2 and the high-potential iron-sulfur protein (HiPIP) function as the electron donors to the Cyt subunit. Most of the Cyt subunits possess four heme c groups, and it was unclear which heme group first accepts the electron from the electron donor. The most distal heme to the special-pair, the heme-1, has a lower redox potential than the electron donors, which makes it difficult to understand the electron transfer mechanism mediated by the Cyt subunit. Extensive mutagenesis combined with kinetic studies has made a great contribution to our understanding of the molecular interaction mechanisms, and has demonstrated the importance of the region close to the heme-1 in the electron transfer. Moreover, crystallographic studies have elucidated two high-resolution three-dimensional structures for the RCs containing the Cyt subunit, the Blastochloris viridis and Thermochromatium tepidum RCs, as well as the structures of their electron donors. An examination of the structural data also suggested that the binding sites for both the cytochrome c2 and the HiPIP are located adjacent to the solvent-accessible edge of the heme-1. In addition, it is also indicated by the structural and biochemical data that the cytochrome c2 and the HiPIP dock with the Cyt subunit by c2 is recognized through electrostatic interactions while hydrophobic interactions are important in the HiPIP docking.

Transition state and encounter complex for fast association of cytochrome c2 with bacterial reaction center

Electrostatic interactions strongly enhance the electron transfer reaction between cytochrome (Cyt) c(2) and reaction center (RC) from photosynthetic bacteria, yielding a second-order rate constant, k(2) approximately 10(9) s(-1).M(-1), close to the diffusion limit. The proposed mechanism involves an encounter complex (EC) stabilized by electrostatic interactions, followed by a transition state (TS), leading to the bound complex active in electron transfer. The effect of electrostatic interactions was previously studied by Tetreault et al. [Tetreault, M., Cusanovich, M., Meyer, T., Axelrod, H. & Okamura, M. Y. (2002) Biochemistry 41, 5807-5815] by measuring k(2) for RC and Cyt molecules with modified charged residues at the binding interface. The present work is a computational analysis of this kinetic study to determine the ensemble of configurations of the TS and EC. Changes in the TS energies due to different mutations were compared with differences in the calculated electrostatic energies for a wide range of Cyt/RC configurations. The TS ensemble, obtained from structures having the highest correlation coefficients in the comparison with experimental data, has the Cyt displaced by approximately 10 A from its position in x-ray crystal structure, close to the average position of the EC ensemble, with strong electrostatic interactions between Cyt on the M subunit side of the RC surface. The heme of the Cyt is oriented toward Tyr L162 on the RC, the tunneling contact in the bound final state on the RC. The similarity between the structures of the EC, TS, and bound state can account for the rapid rate of association responsible for fast diffusion-controlled electron transfer.

Characterization of the interaction of Rhodobacter capsulatus cytochrome c peroxidase with charge reversal mutants of cytochrome c(2)

Steady-state kinetics for the reaction of Rhodobacter capsulatus bacterial cytochrome c peroxidase (BCCP) with its substrate cytochrome c(2) were investigated. The Rb. capsulatus BCCP is dependent on calcium for activation as previously shown for the Pseudomonas aeruginosa BCCP and Paracoccus denitrificans enzymes. Furthermore, the activity shows a bell-shaped pH dependence with optimum at pH 7.0. Enzyme activity is greatest at low ionic strength and drops off steeply as ionic strength increases, resulting in an apparent interaction domain charge product of -13. All cytochromes c(2) show an asymmetric distribution of surface charge, with a concentration of 14 positive charges near the exposed heme edge of Rb. capsulatus c(2) which potentially may interact with approximately 6 negative charges, localized near the edge of the high-potential heme of the Rb. capsulatus BCCP. To test this proposal, we constructed charge reversal mutants of the 14 positively charged residues located on the front face of Rb. capsulatus cytochrome c(2) and examined their effect on steady-state kinetics with BCCP. Mutated residues in Rb. capsulatus cytochrome c(2) that showed the greatest effects on binding and enzyme activity are K12E, K14E, K54E, K84E, K93E, and K99E, which is consistent with the site of electron transfer being located at the heme edge. We conclude that a combination of long-range, nonspecific electrostatic interactions as well as localized salt bridges between, e.g., cytochrome c(2) K12, K14, K54, and K99 with BCCP D194, D241, and D6, account for the observed kinetics.

Conservation of the free energy change of the alkaline isomerization in mitochondrial and bacterial cytochromes c

The thermodynamic parameters of the alkaline transition for oxidized native yeast iso-1 cytochrome c and Rhodopseudomonas palustris cytochrome c(2) (cytc(2)) have been determined through direct electrochemistry experiments carried out at variable pH and temperature and compared to those for horse and beef heart cytochromes c. We have found that both transition enthalpy and entropy are remarkably species dependent, following the order R. palustris cytc(2) >> beef (horse) heart cytc>yeast iso-1 cytc. Considering the high homology at the heme-protein interface in the native species, this variability is likely to be mainly determined by differences in the structural and solvation properties and the relative abundance of the various alkaline conformers. Notably, changes in transition enthalpy and entropy among these cytochromes c are compensative and result in small variations in the free energy change of the process (which amounts approximately to +50 kJ mol(-1)) and consequently in the apparent pK(a) value. This compensation indicates that solvent reorganization effects play an important role in the thermodynamics of the transition. This mechanism is functional to ensure a relatively high pK(a) value for the alkaline transition, which is needed to preserve His,Met ligation to the heme iron in cytochrome c at physiological pH and temperature, hence the E(o) value required for the biological function.

Double mutant studies identify electrostatic interactions that are important for docking cytochrome c2 onto the bacterial reaction center

Cytochrome c2 (cyt) is the mobile electron donor to the reaction center (RC) in photosynthetic bacteria. The electrostatic interactions involved in the dynamics of docking of cyt onto the RC were examined by double mutant studies of the rates of electron transfer between six modified Rhodobacter sphaeroides RCs in which negatively charged acid residues were replaced with Lys and five modified Rhodobacter capsulatus Cyt c2 molecules in which positively charged Lys residues were replaced with Glu. We measured the second-order rate constant, k2, for electron transfer from the reduced cyt to the oxidized primary donor on the RC, which reflects the energy of the transition state for the formation of the active electron transfer complex. Strong interactions were found between Lys C99 and Asp M184/Glu M95, and between Lys C54 and Asp L261/Asp L257. The interacting residues were found to be located close to each other in the recently determined crystal structure of the cyt-RC complex [Axelrod, H., et al. (2002) J. Mol. Biol. (in press)]. The interaction energies were approximately inversely proportional to the distances between charges. These results support earlier suggestions [Tetreault, M., et al. (2001) Biochemistry 40, 8452-8462] that the structure of the transition state in solution resembles the structure of the cyt-RC complex in the cocrystal and indicate that specific electrostatic interactions facilitate docking of the cyt onto the RC in a configuration optimized for both binding and electron transfer. The specific interaction between Asp M184 and Lys C99 may help to nucleate short-range hydrophobic contacts.

Cleavage of the iron-methionine bond in c-type cytochromes: crystal structure of oxidized and reduced cytochrome c(2) from Rhodopseudomonas palustris and its ammonia complex

The three-dimensional structures of the native cytochrome c(2) from Rhodopseudomonas palustris and of its ammonia complex have been obtained at pH 4.4 and pH 8.5, respectively. The structure of the native form has been refined in the oxidized state at 1.70 A and in the reduced state at 1.95 A resolution. These are the first high-resolution crystal structures in both oxidation states of a cytochrome c(2) with relatively high redox potential (+350 mV). The differences between the two oxidation states of the native form, including the position of internal water molecules, are small. The unusual six-residue insertion Gly82-Ala87, which precedes the heme binding Met93, forms an isolated 3(10)-helix secondary structural element not previously observed in other c-type cytochromes. Furthermore, this cytochrome shows an external methionine residue involved in a strained folding near the exposed edge of the heme. The structural comparison of the present cytochrome c(2) with other c-type cytochromes has revealed that the presence of such a residue, with torsion angles phi and psi of approximately -140 and -130 degrees, respectively, is a typical feature of this family of proteins. The refined crystal structure of the ammonia complex, obtained at 1.15 A resolution, shows that the sulphur atom of the Met93 axial ligand does not coordinate the heme iron atom, but is replaced by an exogenous ammonia molecule. This is the only example so far reported of an X-ray structure with the heme iron coordinated by an ammonia molecule. The detachment of Met93 is accompanied by a very localized change in backbone conformation, involving mainly the residues Lys92, Met93, and Thr94. Previous studies under typical denaturing conditions, including high-pH values and the presence of exogenous ligands, have shown that the detachment of the Met axial ligand is a basic step in the folding/unfolding process of c-type cytochromes. The ammonia adduct represents a structural model for this important step of the unfolding pathway. Factors proposed to be important for the methionine dissociation are the strength of the H-bond between the Met93 and Tyr66 residues that stabilizes the native form, and the presence in this bacterial cytochrome c(2) of the rare six-residue insertion in the helix 3(10) conformation that increases Met loop flexibility.

An Alternative to the Accepted Phylogeny of Purple Bacteria Based on 16S rRNA: Analyses of the Amino Acid Sequences of Cytochromes C2 and C556 from Rhodobacter (Rhodovulum) sulfidophilus

It is becoming increasingly apparent from complete genome sequences that 16S rRNA data, as currently interpreted, does not provide an unambiguous picture of bacterial phylogeny. In contrast, we have found that analysis of insertions and deletions in the amino acid sequences of cytochrome c2 has some advantages in establishing relationships and that this approach may have broad utility in acquiring a better understanding of bacterial relationships. The amino acid sequences of cytochromes c2 and c556 have been determined in whole or in part from four strains of Rhodobacter sulfidophilus. The cytochrome c2 contains three- and eight-residue insertions as well as a single-residue deletion in common with the large cytochromes c2 but in contrast to the small cytochromes c2 and mitochondrial cytochromes. In addition, the Rb. sulfidophilus protein shares a rare six- to seven-residue insertion with other Rhodobacter cytochromes c2. The cytochrome c556 is a low-spin class II cytochrome c homologous to the greater family of cytochromes c', which are usually high-spin. The similarity of cytochrome c556 to other species of class II cytochromes is consistent with the relationships deduced from comparisons of cytochromes c2. Thus, our results do not support placement of Rb. sulfidophilus in a separate genus, Rhodovulum, which was proposed primarily on the basis of 16S rRNA sequences. Instead, the Rhodobacter cytochromes c2 are distinct from those of other genera and species of purple bacteria and show a different pattern of relationships among species than reported for 16S rRNA.

Crystal structure of the oxidized cytochrome c(2) from Blastochloris viridis

The crystal structure of the oxidized cytochrome c(2) from Blastochloris (formerly Rhodopseudomonas) viridis was determined at 1.9 A resolution. Structural comparison with the reduced form revealed significant structural changes according to the oxidation state of the heme iron. Slight perturbation of the polypeptide chain backbone was observed, and the secondary structure and the hydrogen patterns between main-chain atoms were retained. The oxidation state-dependent conformational shifts were localized in the vicinity of the methionine ligand side and the propionate group of the heme. The conserved segment of the polypeptide chain in cytochrome c and cytochrome c(2) exhibited some degree of mobility, interacting with the heme iron atom by the hydrogen bond network. These results indicate that the movement of the internal water molecule conserved in various c-type cytochromes drives the adjustments of side-chain atoms of nearby residue, and the segmental temperature factor changes along the polypeptide chain.

Redox-related conformational changes in Rhodobacter capsulatus cytochrome c2

WEFT-NOESY and transfer WEFT-NOESY NMR spectra were used to determine the heme proton assignments for Rhodobacter capsulatus ferricytochrome c2. The Fermi contact and pseudo-contact contributions to the paramagnetic effect of the unpaired electron in the oxidized state were evaluated for the heme and ligand protons. The chemical shift assignments for the 1H and 15N NMR spectra were obtained by a combination of 1H-1H and 1H-15N two-dimensional NMR spectroscopy. The short-range nuclear Overhauser effect (NOE) data are consistent with the view that the secondary structure for the oxidized state of this protein closely approximates that of the reduced form, but with redox-related conformational changes between the two redox states. To understand the decrease in stability of the oxidized state of this cytochrome c2 compared to the reduced form, the structural difference between the two redox states were analyzed by the differences in the NOE intensities, pseudo-contact shifts and the hydrogen-deuterium exchange rates of the amide protons. We find that the major difference between redox states, although subtle, involve heme protein interactions, orientation of the heme ligands, differences in hydrogen bond networks and, possible alterations in the position of some internal water molecules. Thus, it appears that the general destabilization of cytochrome c2, which occurs on oxidation, is consistent with the alteration of hydrogen bonds that result in changes in the internal dynamics of the protein.

Structure and stability effects of the mutation of glycine 34 to serine in Rhodobacter capsulatus cytochrome c(2)

Gly 34 and the adjacent Pro 35 of Rhodobacter capsulatus cytochrome c(2) (or Gly 29 and Pro 30 in vertebrate cytochrome c) are highly conserved side chains among the class I c-type cytochromes. The mutation of Gly 34 to Ser in Rb. capsulatus cytochrome c(2) has been characterized in terms of physicochemical properties and NMR in both redox states. A comparison of the wild-type cytochrome c(2), the G34S mutation, and the P35A mutation is presented in the context of differences in chemical shifts, the differences in NOE patterns, and structural changes resulting from oxidation of the reduced cytochrome. G34S is substantially destabilized relative to wild-type (2.2 kcal/mol in the oxidized state) but similarly destabilized relative to P35A. Nevertheless, differences in terms of the impact of the mutations on specific structural regions are found when comparing G34S and P35A. Although available data indicates that the overall secondary structure of G34S and wild-type cytochrome c(2) are similar, a number of both perturbations of hydrogen bond networks and interactions with internal waters are found. Thus, the impact of the mutation at position 35 is propagated throughout the cytochrome but with alterations at defined sites within the molecule. Interestingly, we find that the substitution of serine at position 34 results in a perturbation of the heme beta meso and the methyl-5 protons. This suggests that the hydroxyl and beta carbon are positioned away from the solvent and toward the heme. This has the consequence of preferentially stabilizing the oxidized state in G34S, thus, altering hydrogen bond networks which involve the heme propionate, internal waters, and key amino acid side chains. The results presented provide important new insights into the stability and solution structure of the cytochrome c(2).

Protein dynamics: imidazole binding to class I C-type cytochromes

The oxidized cytochrome c(2) from the purple phototrophic bacteria, Rhodobacter sphaeroides and Rhodobacter capsulatus, bind the neutral species of imidazole (K(a) = 1440 +/- 40 M(-1)) 50 times more strongly than does horse mitochondrial cytochrome c (K(a) = 30 +/- 1 M(-1)). The kinetics of imidazole binding are consistent with a change in rate-limiting step at high ligand concentrations for all three proteins. This is attributed to a conformational change leading to breakage of the iron-methionine bond which precedes imidazole binding. The three-dimensional structure of the Rb. sphaeroides cytochrome c(2) imidazole complex (Axelrod et al., Acta Crystalogr. D50, 596-602) supports the view that the conformational changes are essentially localized to approximately seven residues on either side of the ligated methionine and there is a hydrogen bond between the Phe 102 carbonyl, an internal water, and the bound imidazole. Insertions and deletions in this region of cytochrome c(2), the presence of a proline near the methionine, and the smaller size of the dynamic region of horse cytochrome c suggest that the stabilizing hydrogen bond is not present in horse cytochrome c, hence, the dramatic difference in affinity for imidazole. The kinetics of ligand binding do not correlate with either the strength of the iron-methionine bond as measured by the pK of the 695-nm absorption band or the overall stability of the cytochromes studied. However, the very similar imidazole binding properties of the two cytochromes c(2) indicate that the Rb. sphaeroides cytochrome c(2)-imidazole complex structure is an excellent model for the corresponding Rb. capsulatus cytochrome c(2) complex. It is notable that the movement of the peptide chain in the vicinity of the ligated methionine has been preserved throughout evolution and suggests a role in the function of c-type cytochromes.

Distinct structures and environments for the three hemes of the cytochrome bc1 complex from Rhodospirillum rubrum. A resonance Raman study using B-band excitations

The B-band excited resonance Raman (RR) spectra (100-1700 cm-1) of the bacterial cytochrome bc1 complex purified from Rhodospirillum rubrum are reported. Four redox states, i.e., the persulfate-oxidized, "as prepared", and ascorbate- and dithionite-reduced states of the complex, were investigated with the laser excitations at 406.7, 413.1, and 441.6 nm. Following the different absorption properties of the b- and c-type hemes and the different resonance enhancements of the vibrational modes of oxidized and reduced hemes, RR contributions from the b- and c-type hemes were characterized. For the nu2, nu10, and nu8 porphyrin vibrational modes, individual contributions of hemes c1, bH, and bL were determined. The data show that the macrocycle conformation of the three hemes of the cytochrome bc1 complex is different. In particular, the frequencies assigned to ferrous heme bL (1580, 1610, and 352 cm-1, respectively) reveal that its porphyrin is more strongly distorted than that of ferrous heme bH (1584, 1614, and 344 cm-1, respectively). The frequencies of the nu11 modes (1543, 1536, and 1526 cm-1 for ferrous heme c1, heme bH, and heme bL, respectively) confirm that the axial histidylimidazole ligands of heme bL have a marked anionic character. Strong differences in the peripheral interactions of the three hemes with the proteins were also detected through the frequency differences of the nu5, nu13, nu14, and nu42 modes. Considering that hemes bH and bL are inserted into a four-helice bundle, the RR data are interpreted in the frame of a strong protein constraint on heme bL.

Imidazole binding to Rhodobacter capsulatus cytochrome c2. Effect of site-directed mutants on ligand binding

Although ligand binding in c-type cytochromes is not directly related to their physiological function, it has the potential to provide valuable information on protein stability and dynamics, particularly in the region of the methionine sixth heme ligand and the nearby peptide chain that has been implicated in electron transfer. Thus, we have measured the equilibrium and kinetics of binding of imidazole to eight mutants of Rhodobacter capsulatus cytochrome c2 that differ in overall protein stability. We found that imidazole binding affinity varies 70-fold, but does not correlate with overall protein stability. Instead, each mutant exerts an effect at the local level, with the largest change due to mutant G95E (glycine substituted by glutamate), which shows 30-fold stronger binding as compared with the wild-type protein. The kinetics of imidazole binding are monophasic and reach saturation at high ligand concentrations for all the mutants and wild-type protein, which is attributed to a rate-limiting conformational change leading to breakage of the iron-methionine bond and providing a binding site for imidazole. The mutants show as much as an 18-fold variation in the first-order rate constant for the conformational change, with the largest effect found with mutant G95E. The kinetics also show a lack of correlation with overall protein stability, but are consistent with localized effects on the dynamics of hinge region 88-102 of the protein, which changes conformation to permit ligand binding. These results are consistent with R. capsulatus cytochrome c2 stabilizing the complex through hydrogen bonding to the imidazole. The larger effects of mutant G95E on equilibrium and kinetics are likely to be due to its location within the hinge region adjacent to heme ligand methionine 96, which is displaced by imidazole.

Solution structure, rotational diffusion anisotropy and local backbone dynamics of Rhodobacter capsulatus cytochrome c2

The solution structure, backbone dynamics and rotational diffusion of the Rhodobacter capsulatus cytochrome c2 have been determined using heteronuclear NMR spectroscopy. In all, 1204 NOE-derived distances were used in the structure calculation to give a final ensemble with 0.59(+/-0.08) A rms deviation for the backbone atoms (C, Calpha and N) with respect to the mean coordinates. There is no major difference between the solution structure and the previously solved X-ray crystal structure (1.07(+/-0.07) A rms difference for the backbone atoms), although certain significant local structural differences have been identified. This protein contains five helical regions and a histidine-heme binding domain, connected by a series of structured loops. The orientation of the helices provides an excellent sampling of angular space and thus allows a precise characterization of the anisotropic diffusion tensor. Analysis of the hydrodynamics of the protein has been performed by interpretation of the 15N relaxation data using isotropic, axially asymmetric and fully anisotropic diffusion tensors. The protein can be shown to exhibit significant anisotropic reorientation with a diffusion tensor with principal axes values of 1.405(+/-0.031)x10(7) s-1, 1.566(+/-0.051)x10(7) s-1 and 1.829(+/-0.054)x10(7) s-1. Hydrodynamic calculations performed on the solution structure predict values of 1.399x10(7) s-1, 1.500x10(7) s-1 and 1.863x10(7) s-1 when a solvent shell of 3.5 A is included in the calculation. The optimal orientation of the diffusion tensor has been incorporated into a hybrid Lipari-Szabo type local motion-anisotropic rotational diffusion model to characterize the local mobility in the molecule. The mobility parameters thus extracted show a quantitative improvement with respect to the model-free analysis assuming isotropic reorientation; helical regions exhibit similar dynamic properties and fewer residues require more complex models of internal motion. While the molecule is essentially rigid, a tripeptide loop region (residues 101 to 103) exhibits flexibility in the range of 20 to 30 ps, which appears to be correlated with the order in the NMR solution structure.

Molecular mechanisms of cytochrome c biogenesis: three distinct systems

The past 10 years have heralded remarkable progress in the understanding of the biogenesis of c-type cytochromes. The hallmark of c-type cytochrome synthesis is the covalent ligation of haem vinyl groups to two cysteinyl residues of the apocytochrome (at a Cys-Xxx-Yyy-Cys-His signature motif). From genetic, genomic and biochemical studies, it is clear that three distinct systems have evolved in nature to assemble this ancient protein. In this review, common principles of assembly for all systems and the molecular mechanisms predicted for each system are summarized. Prokaryotes, plant mitochondria and chloroplasts use either system I or II, which are each predicted to use dedicated mechanisms for haem delivery, apocytochrome ushering and thioreduction. Accessory proteins of systems I and II co-ordinate the positioning of these two substrates at the membrane surface for covalent ligation. The third system has evolved specifically in mitochondria of fungi, invertebrates and vertebrates. For system III, a pivotal role is played by an enzyme called cytochrome c haem lyase (CCHL) in the mitochondrial intermembrane space.

Purification and primary structure analysis of two cytochrome c2 isozymes from the purple phototrophic bacterium Rhodospirillum centenum

The isolation and amino acid sequences of two cytochromes c-552 from the thermotolerant bacterium Rhodospirillum (R.) centenum have been determined. They are very similar to one another with 85% identity. They are homologous to the cytochromes c2 from purple bacteria with approximately 67% identity to that from Rhodopseudomonas (Rps.) palustris compared to only 42% identity with others of the c2 subclass. In addition, they share an unusual six-residue insertion with Rps. palustris cytochrome c2 not found in any other cytochrome. The relationship with Rps. palustris is thus highly significant. The redox potentials of the R. centenum isozymes are 293 and 316 mV. Although the proteins have strongly different iso-electric points, both have three conserved lysine residues at the proposed site of electron transfer. These results suggest that they may be functionally interchangeable.

Protein folding and protein evolution: common folding nucleus in different subfamilies of c-type cytochromes?

Amino acid sequences of seven subfamilies of cytochromes c (mitochondrial cytochromes c, c1; chloroplast cytochromes c6, cf; bacterial cytochromes c2, c550, c551; in total 164 sequences) have been compared. Despite extensive homology within eukaryotic subfamilies, homology between different subfamilies is very weak. Other than the three heme-binding residues (Cys13, Cys14, His18, in numeration of horse cytochrome c) there are only four positions which are conserved in all subfamilies: Gly/Ala6, Phe/Tyr10, Leu/Val/Phe94 and Tyr/Trp/Phe97. In all 17 cytochromes c with known 3D-structures, these residues form a network of conserved contacts (6-94, 6-97, 10-94, 10-97 and 94-97). Especially strong is the contact between aromatic groups in positions 10 and 97, which corresponds to 13 interatomic contacts. As residues 6, 10 and residues 94, 97 are in (i, i+4) and (i, i+3) positions in the N and C-terminal helices, respectively, the above mentioned system of conserved contacts consists mainly of contacts between one turn of N-terminal helix and one turn of C-terminal helix. The importance of the contacts between interfaces of these helices has been confirmed by the existence of these contacts in both equilibrium and kinetic molten globule-like folding intermediates, as well as by mutational evidence that these contacts are involved in tight packing between the N and C-helices. Since these four residues are not involved in heme binding and have no other apparent functional role, their conservation in highly diverged cytochromes c suggests that they are of a critical importance for protein folding. The author assumes that they are involved in a common folding nucleus of all subfamilies of c-type cytochromes.

Conservation of the conformation of the porphyrin macrocycle in hemoproteins

The out-of-plane distortions of porphyrins in hemoproteins are characterized by displacements along the lowest-frequency out-of-plane normal coordinates of the D4h-symmetric macrocycle. X-ray crystal structures are analyzed using a computational procedure developed for determining these orthogonal displacements. The x-ray crystal structures of the heme groups are described within experimental error, using the set composed of only the lowest frequency normal coordinate of each out-of-plane symmetry type. That is, the distortion is accurately simulated by a linear combination of these orthonormal deformations, which include saddling (B2u), ruffling (B1u), doming (A2u), waving (Eg), and propellering (A1u). For example, orthonormal structural decomposition of the hemes in deoxymyoglobins reveals a predominantly dom heme deformation combined with a smaller wav(y) deformation. Generally, the heme conformation is remarkably similar for proteins from different species. For cytochromes c, the conformation is conserved as long as the amino acids between the cysteine linkages to the heme are homologous. Differences occur if this short segment varies in the number or type of residues, suggesting that this small segment causes the nonplanar distortion. Some noncovalently linked hemes like those in the peroxidases also have highly conserved characteristic distortions. Conservation occurs even for some proteins with a large natural variation in the amino acid sequence.

Redox thermodynamics of the native and alkaline forms of eukaryotic and bacterial class I cytochromes c

The reduction potentials of beef heart cytochrome c and cytochromes c2 from Rhodopseudomonas palustris, Rhodobacter sphaeroides, and Rhodobacter capsulatus were measured through direct electrochemistry at a surface-modified gold electrode as a function of temperature in nonisothermal experiments carried out at neutral and alkaline pH values. The thermodynamic parameters for protein reduction (DeltaS degrees rc and DeltaH degrees rc) were determined for the native and alkaline conformers. Enthalpy and entropy terms underlying species-dependent differences in E degrees and pH- and temperature-induced E degrees changes for a given cytochrome were analyzed. The difference of about +0.1 V in E degrees between cytochromes c2 and the eukaryotic species can be separated into an enthalpic term (-DeltaDeltaH degrees rc/F) of +0.130 V and an entropic term (TDeltaDeltaS degrees rc/F) of -0.040 V. Hence, the higher potential of the bacterial species appears to be determined entirely by a greater enthalpic stabilization of the reduced state. Analogously, the much lower potential of the alkaline conformer(s) as compared to the native species is by far enthalpic in origin for both protein families, and is largely determined by the substitution of Met for Lys in axial heme ligation. Instead, the biphasic E degrees /temperature profile for the native cytochromes is due to a difference in reduction entropy between the conformers at low and high temperatures. Temperature-dependent 1H NMR experiments suggest that the temperature-induced transition also involves a change in orientation of the axial methionine ligand with respect to the heme plane.

Influence of conserved amino acids on the structure and environment of the heme of cytochrome c2. A resonance Raman study

Resonance Raman spectra using Soret excitations of oxidized and reduced Rhodobacter capsulatus cytochrome c2 at pH 7.5 were studied. The spectra of oxidized cytochrome c2 show three components for the v10 mode at 1638, 1633, and 1629 cm(-1). The intensities of these components are sensitive to the excitation wavelength. This effect is explained in the context of a conformational equilibrium of the ferriheme between a nearly planar structure and two ruffled structures. In the case of reduced cytochrome c2, the absolute frequencies as well as the excitation-dependent frequency dispersion of the v10 mode (1618-1621 cm(-1)) indicate a displacement of the conformational equilibrium of heme toward the more planar structures. To measure the influence of some key amino acid residues on the heme-protein interaction of cytochrome c2, four site-directed mutants of Rb. capsulatus cytochrome c2 have been studied by resonance Raman spectroscopy and their spectra compared with the spectra obtained for the wild type cytochrome. The mutants studied are K14E/K32E, P35A, W67Y, and Y75F. The spectral changes induced by the mutations are interpreted in terms of alterations in the structure and/or environment of the cytochrome c2 heme in the framework of the expected role of the different amino acid residues in the stability and redox potential.

Anion binding to cytochrome c2: implications on protein-Ion interactions in class I cytochromes c

The binding of several inorganic and carboxylate anions to cytochrome c2 from Rhodopseudomonas palustris has been investigated by monitoring the salt-induced changes in the redox potential of the heme, using an interpretative model based on the extended Debye-Hückel equation. Most anions were found to interact specifically with the protein at one or multiple sites. Binding constants to the oxidized protein in the range 10(1)-10(2) m-1 were determined from the anion concentration dependence of the chemical shift of the isotropically shifted heme methyl resonances. For several anions the stoichiometry and strength of the binding to cytochrome c2 were found comparable with those determined for mitochondrial cytochromes c, in spite of the limited sequence similarity (less than 40%) and the lower positive charge of the bacterial protein. These analogies were interpreted as indicative of the existence of common binding sites which are proposed to be located in the conserved lysine-rich domain around the solvent-exposed heme edge, which is also the surface area likely involved in the interaction with redox partners. The changes in E degrees due to partial neutralization of the positive charge of cytochrome c2 due to specific anion binding were found comparable with those for the mitochondrial species.

Kinetic mechanism of folding and unfolding of Rhodobacter capsulatus cytochrome c2

In spite of marginal sequence homology, cytochrome c2 from photosynthetic bacteria and the mitochondrial cytochromes c exhibit some striking structural similarities, including the tertiary arrangement of the three main helices. To compare the folding mechanisms for these two distantly related groups of proteins, equilibrium and kinetic measurements of the folding/unfolding reaction of cytochrome c2 from Rhodobacter capsulatus were performed as a function of guanidine hydrochloride (GuHCl) concentration in the absence and presence of a stabilizing salt, sodium sulfate. Quenching of the fluorescence of Trp67 by the heme was used as a conformational probe. Kinetic complexities due to non-native histidine ligation are avoided, since cytochrome c2 contains only one histidine, His17, which forms the axial heme ligand under native and denaturing conditions. Quantitative kinetic modeling showed that both equilibrium and kinetic results are consistent with a minimal four-state mechanism with two sequential intermediates. The observation of a large decrease in fluorescence during the 2-ms dead-time of the stopped-flow measurement (burst phase) at low GuHCl concentration, followed by a sigmoidal recovery of the initial amplitude toward the unfolding transition region, is attributed to a well-populated compact folding intermediate in rapid exchange with unfolded molecules. A nearly denaturant-independent process at low GuHCl concentrations reflects the rate-limiting conversion of a compact intermediate to the native state. At high GuHCl concentrations, a process with little denaturant dependence is attributed to the rate-limiting Met96-iron deligation process during unfolding, which is supported by the kinetics of imidazole binding. The strong GuHCl-dependence of folding and unfolding rates near the midpoint of the equilibrium transition is attributed to destabilization of each intermediate and their transition states in folding and unfolding. Addition of sodium sulfate shifts the rate profile to higher denaturant concentration, which can be understood in terms of the relative stabilizing effect of the salt on partially and fully folded states.

An optimized g-tensor for Rhodobacter capsulatus cytochrome c2 in solution: a structural comparison of the reduced and oxidized states

The optimized g-tensor parameters for the oxidized form of Rhodobacter capsulatus cytochrome c2 in solution were obtained using a set (50) of backbone amide protons. Dipolar shifts for more than 500 individual protons of R. capsulatus cytochrome c2 have been calculated by using the optimized g-tensor and the X-ray crystallographic coordinates of the reduced form of R. capsulatus cytochrome c2. The calculated results for dipolar shifts are compared with the observed paramagnetic shifts. The calculated and the observed data are in good agreement throughout the entire protein, but there are significant differences between calculated and experimental results localized to the regions in the immediate vicinity of the heme ligand and the region of the front crevice of the protein (residues 44-50, 53-57, and 61-68). The results not only indicate that the overall solution structures are very similar in both the reduced and oxidized states, but that these structures in solution are similar to the crystal structure. However, there are small structural changes near the heme and the rearrangement of certain residues that result in changes in their hydrogen bonding concomitant with the change in the oxidation states; this was also evident in the data for the NH exchange rate measurements for R. capsulatus cytochrome c2.

Stability study of Rhodobacter capsulatus ferrocytochrome c2 wild-type and site-directed mutants using hydrogen/deuterium exchange monitored by electrospray ionization mass spectrometry

To estimate the stability of Rhodobacter capsulatus ferrocytochrome c2 wild-type and site-directed mutants, charge state distributions and hydrogen/deuterium exchange rates were monitored by electrospray ionization mass spectrometry. The relative stability of the mutants was observed with the order: V11 insert > Y75F > wild-type = K32E > K12D = K14E > or = K52E > K14E/K32E > W67Y > P35A > I57N > G34S. (A preliminary account has been presented for mutants G34S and P35A [Jaquinod et al. (1995) Rapid Commun. Mass Spectrom. 9, 1135-1140].) This approach is shown to be a useful tool for rapid characterization of mutational effects on protein conformation.

Comparison of low oxidoreduction potential cytochrome c553 from Desulfovibrio vulgaris with the class I cytochrome c family

The cytochrome c553 from Desulfovibrio vulgaris (DvH c553) is of importance in the understanding of the relationship of structure and function of cytochrome c due to its lack of sequence homology with other cytochromes, and its abnormally low oxido-reduction potential. In evolutionary terms, this protein also represents an important reference point for the understanding of both bacterial and mitochondrial cytochromes c. Using the recently determined nuclear magnetic resonance (NMR) structure of the reduced protein we compare the structural, dynamic, and functional characteristics of DvH c553 with members of both the mitochondrial and bacterial cytochromes c to characterize the protein in the context of the cytochrome c family, and to understand better the control of oxide-reduction potential in electron transfer proteins. Despite the low sequence homology, striking structural similarities between this protein and representatives of both eukaryotic [cytochrome c from tuna (tuna c)] and prokaryotic [Pseudomonas aeruginosa c551 (Psa c551)] cytochromes c have been recognized. The previously observed helical core is also found in the DvH c553. The structural framework and hydrogen bonding network of the DvH c553 is most similar to that of the tuna c, with the exception of an insertion loop of 24 residues closing the heme pocket and protecting the propionates, which is absent in the DvH c553. In contrast, the Psa c551 protects the propionates from the solvent principally by extending the methionine ligand arm. The electrostatic distribution at the recognized encounter surface around the heme in the mitochondrial cytochrome is reproduced in the DvH c553, and corresponding hydrogen bonding networks, particularly in the vicinity of the heme cleft, exist in both molecules. Thus, although the cytochrome DvH c553 exhibits higher primary sequence homology to other bacterial cytochromes c, the structural and physical homology is significantly greater with respect to the mitochondrial cytochrome c. The major structural and functional difference is the absence of solvent protection for the heme, differentiating this cytochrome from both reference cytochromes, which have evolved different mechanisms to cover the propionates. This suggests that the abnormal redox potential of the DvH c553 is linked to the raised accessibility of the heme and supports the theory that redox potential in cytochromes is controlled by heme propionate solvent accessibility.

Knowledge-based protein secondary structure assignment

We have developed an automatic algorithm STRIDE for protein secondary structure assignment from atomic coordinates based on the combined use of hydrogen bond energy and statistically derived backbone torsional angle information. Parameters of the pattern recognition procedure were optimized using designations provided by the crystallographers as a standard-of-truth. Comparison to the currently most widely used technique DSSP by Kabsch and Sander (Biopolymers 22:2577-2637, 1983) shows that STRIDE and DSSP assign secondary structural states in 58 and 31% of 226 protein chains in our data sample, respectively, in greater agreement with the specific residue-by-residue definitions provided by the discoverers of the structures while in 11% of the chains, the assignments are the same. STRIDE delineates every 11th helix and every 32nd strand more in accord with published assignments.

Cyclic voltammetry and 1H-NMR of Rhodopseudomonas palustris cytochrome c2 pH-dependent conformational states

The pH-induced protein conformational transitions and changes in the ligation state of the heme iron in cytochrome c2 from Rhodopseudomonas palustris were monitored by electrochemical and spectroscopic measurements. In the pH range 1.5-11, the E degree values (and/or the peak potentials) determined by cyclic voltammetry, the electronic spectra and the hyperfine-shifted 1H-NMR resonances of the protein are sensitive to a number of acid/base equilibria. In particular, four equilibria have been determined for the oxidized protein with pKa values of 2.5, 5.5, 6.6 and 9.0. The lowest pKa most probably involves disruption of both axial heme iron bonds and protein unfolding. The subsequent pKa is associated with a low-pH oxidation of the protein by dioxygen, which is accompanied by a conformational change. The equilibrium with an apparent pKa of 6.6 modulates the E degree values without determining any detectable spectral change and most likely involves the acid/base equilibrium of an histidine residue in close vicinity of the heme (possibly His53). Finally, the alkaline ionization is due to the replacement of the methionine axially bound to the heme iron with a stronger (most probably N-donor) ligand. The reduced alkaline form is unstable and spontaneously converts to the neutral reduced form with a kinetic constant of 0.98 s-1 at pH 9.2.

1H and 13C NMR assignments and structural aspects of a ferrocytochrome c-551 from the purple phototrophic bacterium Ectothiorhodospira halophila

Two-dimensional nuclear magnetic resonance was used to assign the 1H and 13C resonances of ferrocytochrome c-551 from Ectothiorhodospira halophila, a halophilic phototrophic purple bacterium. This 78-residue protein belongs to a small subgroup of class I cytochromes c together with the analogous cytochromes c-551 from E. halochloris and E. abdelmalekii. A nearly complete assignment of 13C resonances was obtained at natural abundance using a gradient-enhanced 1H-detected heteronuclear single quantum coherence experiment (HSQC). This was found to be extremely useful for the unambigous assignment of side chain protons. The secondary structure of the protein was determined from analyses of short- and medium-range nuclear Overhauser enhancements (NOE), amide proton exchange and 13C alpha chemical shifts. Three helices could be identified which are well conserved among the class I cytochromes c. There is some evidence for two other regions of less well defined helical structure. From a preliminary analysis of long-range NOE it is shown that in the E. halophila cytochrome c-551 the general cytochrome c fold is well conserved, including the three conserved helices (residues 2-8, 41-50, 63-76), the regions around the heme ligands (Cys14-Ser15-Ser16-Cys17-His18, Met55) and the omega loop (residues 18-28). In addition, three variable segments of the protein are discussed in detail, one of those including a cis-proline, a feature so far unique in the cytochrome c family. Structural alignments of the E. halophila cytochrome c-551 with two other Pseudomonas cytochrome c5 homologs (Azotobacter vinelandii cytochrome c5 and Chlorobium limicola cytochrome c-555) are provided which are based on sequence similarities and secondary structure alignments.