Assignment of 15N chemical shifts and 15N relaxation measurements for oxidized and reduced iso-1-cytochrome c

Anionic Lipids Confine Cytochrome c2 to the Surface of Bioenergetic Membranes without Compromising Its Interaction with Redox Partners

Cytochrome c₂ (cyt. c₂) is a major element in electron transfer between redox proteins in bioenergetic membranes. While the interaction between cyt. c₂ and anionic lipids abundant in bioenergetic membranes has been reported, their effect on the shuttling activity of cyt. c₂ remains elusive. Here, the effect of anionic lipids on the interaction and binding of cyt. c₂ to the cytochrome bc₁ complex (bc₁) is investigated using a combination of molecular dynamics (MD) and Brownian dynamics (BD) simulations. MD is used to generate thermally accessible conformations of cyt. c₂ and membrane-embedded bc₁, which were subsequently used in multireplica BD simulations of diffusion of cyt. c₂ from solution to bc₁, in the presence of various lipids. We show that, counterintuitively, anionic lipids facilitate association of cyt. c₂ with bc₁ by localizing its diffusion to the membrane surface. The observed lipid-mediated bc₁ association is further enhanced by the oxidized state of cyt. c₂, in line with its physiological function. This lipid-mediated enhancement is salinity-dependent, and anionic lipids can disrupt cyt. c₂-bc₁ interaction at nonphysiological salt levels. Our data highlight the importance of the redox state of cyt. c₂, the lipid composition of the chromatophore membrane, and the salinity of the chromatophore in regulating the efficiency of the electron shuttling process mediated by cyt. c₂. The conclusions can be extrapolated to mitochondrial systems and processes, or any bioenergetic membrane, given the structural similarity between cyt. c₂ and bc₁ and their mitochondrial counterparts.

Enhancing the population of the encounter complex affects protein complex formation efficiency

Optimal charge distribution is considered to be important for efficient formation of protein complexes. Electrostatic interactions guide encounter complex formation that precedes the formation of an active protein complex. However, disturbing the optimized distribution by introduction of extra charged patches on cytochrome c peroxidase does not lead to a reduction in productive encounters with its partner cytochrome c. To test whether a complex with a high population of encounter complex is more easily affected by suboptimal charge distribution, the interactions of cytochrome c mutant R13A with wild‐type cytochrome c peroxidase and a variant with an additional negative patch were studied. The complex of the peroxidase and cytochrome c R13A was reported to have an encounter state population of 80%, compared to 30% for the wild‐type cytochrome c. NMR analysis confirms the dynamic nature of the interaction and demonstrates that the mutant cytochrome c samples the introduced negative patch. Kinetic experiments show that productive complex formation is fivefold to sevenfold slower at moderate and high ionic strength values for cytochrome c R13A but the association rate is not affected by the additional negative patch on cytochrome c peroxidase, showing that the total charge on the protein surface can compensate for less optimal charge distribution. At low ionic strength (44 mm), the association with the mutant cytochrome c reaches the same high rates as found for wild‐type cytochrome c, approaching the diffusion limit.

Distal Unfolding of Ferricytochrome c Induced by the F82K Mutation

It is well known that axial coordination of heme iron in mitochondrial cytochrome c has redox-dependent stability. The Met80 heme iron axial ligand in the ferric form of the protein is relatively labile and can be easily replaced by alternative amino acid side chains under non-native conditions induced by alkaline pH, high temperature, or denaturing agents. Here, we showed a redox-dependent destabilization induced in human cytochrome c by substituting Phe82-conserved amino acid and a key actor in cytochrome c intermolecular interactions-with a Lys residue. Introducing a positive charge at position 82 did not significantly affect the structure of ferrous cytochrome c but caused localized unfolding of the distal site in the ferric state. As revealed by ¹H NMR fingerprint, the ferric form of the F82K variant had axial coordination resembling the renowned alkaline species, where the detachment of the native Met80 ligand favored the formation of multiple conformations involving distal Lys residues binding to iron, but with more limited overall structural destabilization.

A Heme Propionate Staples the Structure of Cytochrome c for Methionine Ligation to the Heme Iron

Ligand-switch reactions at the heme iron are common in biological systems, but their mechanisms and the features of the polypeptide fold that support dual ligation are not well understood. In cytochrome c (cyt c), two low-stability loops (Ω-loop C and Ω-loop D) are connected by the heme propionate HP6. At alkaline pH, the native Met80 ligand from Ω-loop D switches to a Lys residue from the same loop. Deprotonation of an as yet unknown group triggers the alkaline transition. We have created the two cyt c variants T49V/K79G and T78V/K79G with altered connections of these two loops to HP6. Electronic absorption, NMR, and EPR studies demonstrate that at pH 7.4 ferric forms of these variants are Lys-ligated, whereas ferrous forms maintain the native Met80 ligation. Measurements of protein stability, cyclic voltammetry, pH-jump and gated electron-transfer kinetics have revealed that these Thr to Val substitutions greatly affect the alkaline transition in both ferric and ferrous proteins. The substitutions modify the stability of the Met-ligated species and reduction potentials of the heme iron. The kinetics of ligand-switch processes are also altered, and analyses of these effects implicate redox-dependent differences in metal-ligand interactions and the role of the protein dynamics, including cross-talk between the two Ω-loops. With the two destabilized variants, it is possible to map energy levels for the Met- and Lys-ligated species in both ferric and ferrous proteins and assess the role of the protein scaffold in redox-dependent preferences for these two ligands. The estimated shift in the heme iron reduction potential upon deprotonation of the "trigger" group is consistent with those associated with deprotonation of an HP, suggesting that HP6, on its own or as a part of a hydrogen-bonded cluster, is a likely "trigger" for the Met to Lys ligand switch.

About the use of 13C-13C NOESY in bioinorganic chemistry

Herein we present examples of the application of the ¹³C-¹³C Nuclear Overhauser Effect Spectroscopy (NOESY) experiment to the study of metalloproteins and we critically discuss the advantages and drawbacks of the method as a function of the molecular size of the investigated systems. The contribution is focused on a few case studies among the systems analyzed in the group of the corresponding author. The ¹³C-¹³C NOESY experiment represents the gold standard for the observation of NMR signals in the 480 kDa ferritin nanocage and for monitoring its interaction with iron. By decreasing the protein size, the experiment progressively loses its importance as a tool for the detection of the complete spin pattern of the amino acid side chains, as exemplified by nickel-dependent regulatory protein, NikR (molecular mass of the homo-tetramer ~80 kDa). In very small proteins, such as mitochondrial cytochrome c (12.3 kDa), we are only able to detect cross peaks between adjacent ¹³C nuclei; this feature turned out to be useful for the assignment of the ¹³C core resonances of the porphyrin in a uniformly enriched heme.

Histidine-Lysine Axial Ligand Switching in a Hemoglobin: A Role for Heme Propionates

The hemoglobin of Synechococcus sp. PCC 7002, GlbN, is a monomeric group I truncated protein (TrHb1) that coordinates the heme iron with two histidine ligands at neutral pH. One of these is the distal histidine (His46), a residue that can be displaced by dioxygen and other small molecules. Here, we show with mutagenesis, electronic absorption spectroscopy, and nuclear magnetic resonance (NMR) spectroscopy that at high pH and exclusively in the ferrous state, Lys42 competes with His46 for the iron coordination site. When b heme is originally present, the population of the lysine-bound species remains too small for detailed characterization; however, the population can be increased significantly by using dimethyl-esterified heme. Electronic absorption and NMR spectroscopies showed that the reversible ligand switching process occurs with an apparent pK_a of 9.3 and a Lys-ligated population of ∼60% at the basic pH limit in the modified holoprotein. The switching rate, which is slow on the chemical shift time scale, was estimated to be 20-30 s^-1 by NMR exchange spectroscopy. Lys42-His46 competition and attendant conformational rearrangement appeared to be related to weakened bis-histidine ligation and enhanced backbone dynamics in the ferrous protein. The pH- and redox-dependent ligand exchange process observed in GlbN illustrates the structural plasticity allowed by the TrHb1 fold and demonstrates the importance of electrostatic interactions at the heme periphery for achieving axial ligand selection. An analogy is drawn to the alkaline transition of cytochrome c, in which Lys-Met competition is detected at alkaline pH, but, in contrast to GlbN, in the ferric state only.

Increased dynamics in the 40-57 Ω-loop of the G41S variant of human cytochrome c promote its pro-apoptotic conformation

Thrombocytopenia 4 is an inherited autosomal dominant thrombocytopenia, which occurs due to mutations in the human gene for cytochrome c that results in enhanced mitochondrial apoptotic activity. The Gly41Ser mutation was the first to be reported. Here we report stopped-flow kinetic studies of azide binding to human ferricytochrome c and its Gly41Ser variant, together with backbone amide H/D exchange and ¹⁵N-relaxation dynamics using NMR spectroscopy, to show that alternative conformations are kinetically and thermodynamically more readily accessible for the Gly41Ser variant than for the wild-type protein. Our work reveals a direct conformational link between the 40–57 Ω-loop in which residue 41 resides and the dynamical properties of the axial ligand to the heme iron, Met80, such that the replacement of glycine by serine promotes the dissociation of the Met80 ligand, thereby increasing the population of a peroxidase active state, which is a key non-native conformational state in apoptosis.

Protein docking using an ensemble of spin labels optimized by intra-molecular paramagnetic relaxation enhancement

The effect of spin label mobility on the accuracy of protein–protein docking calculations was investigated using inter- and intra-molecular PRE data.

Comparison of the backbone dynamics of wild-type Hydrogenobacter thermophilus cytochrome c(552) and its b-type variant

Cytochrome c552 from the thermophilic bacterium Hydrogenobacter thermophilus is a typical c-type cytochrome which binds heme covalently via two thioether bonds between the two heme vinyl groups and two cysteine thiol groups in a CXXCH sequence motif. This protein was converted to a b-type cytochrome by substitution of the two cysteine residues by alanines (Tomlinson and Ferguson in Proc Natl Acad Sci USA 97:5156-5160, 2000a). To probe the significance of the covalent attachment of the heme in the c-type protein, (15)N relaxation and hydrogen exchange studies have been performed for the wild-type and b-type proteins. The two variants share very similar backbone dynamic properties, both proteins showing high (15)N order parameters in the four main helices, with reduced values in an exposed loop region (residues 18-21), and at the C-terminal residue Lys80. Some subtle changes in chemical shift and hydrogen exchange protection are seen between the wild-type and b-type variant proteins, not only for residues at and neighbouring the mutation sites, but also for some residues in the heme binding pocket. Overall, the results suggest that the main role of the covalent linkages between the heme group and the protein chain must be to increase the stability of the protein.

The IR-¹⁵N-HSQC-AP experiment: a new tool for NMR spectroscopy of paramagnetic molecules

A crucial factor for the understanding of structure-function relationships in metalloproteins is the identification of NMR signals from residues surrounding the metal cofactor. When the latter is paramagnetic, the NMR information in the proximity of the metal center may be scarce, because fast nuclear relaxation quenches signal intensity and coherence transfer efficiency. To identify residues at a short distance from a paramagnetic center, we developed a modified version of the ¹⁵N-HSQC experiment where (1) an inversion recovery filter is added prior to HSQC, (2) the INEPT period has been optimized according to fast relaxation of interested spins, (3) the inverse INEPT has been eliminated and signals acquired as antiphase doublets. The experiment has been successfully tested on a human [Fe₂S₂] protein which is involved in the biogenesis of iron-sulfur proteins. Thirteen HN resonances, unobserved with conventional HSQC experiments, could be identified. The structural arrangement of the protein scaffold in the proximity of the Fe/S cluster is fundamental to comprehend the molecular processes responsible for the transfer of Fe/S groups in the iron-sulfur protein assembly machineries.

Redox state dependence of axial ligand dynamics in Nitrosomonas europaea cytochrome c552

Analysis of NMR spectra reveals that the heme axial Met ligand orientation and dynamics in Nitrosomonas europaea cytochrome c552 (Ne cyt c) are dependent on the heme redox state. In the oxidized state, the heme axial Met is fluxional, interconverting between two conformers related to each other by inversion through the Met δS atom. In the reduced state, there is no evidence of fluxionality, with the Met occupying one conformation similar to that seen in the homologous Pseudomonas aeruginosa cytochrome c551. Comparison of the observed and calculated pseudocontact shifts for oxidized Ne cyt c using the reduced protein structure as a reference structure reveals a redox-dependent change in the structure of the loop bearing the axial Met (loop 3). Analysis of nuclear Overhauser effects (NOEs) and existing structural data provides further support for the redox state dependence of the loop 3 structure. Implications for electron transfer function are discussed.

Mechanistic insights into the superoxide-cytochrome c reaction by lysine surface scanning

This study summarizes results which have been obtained by a mutational study of human cytochrome c. The protein can be used as a recognition element in analytical assays and biosensors for superoxide radicals since ferricytochrome c reacts with superoxide to form ferrocytochrome c and oxygen. Here lysine mutagenesis of the distal surface (i.e., of exposed residues around the Met80 axial ligand) of human cytochrome c was pursued to evaluate the effect of the surface charges on the reaction rate with the superoxide anion radical and on the redox properties of the heme protein. The latter feature is particularly relevant when the protein is immobilized on a negatively charged self-assembled monolayer on an electrode to be used as a biosensor. The observed effects of the mutations are rationalized through structural investigations based on solution NMR spectroscopy and computational analysis of the surface electrostatics. The results suggest the presence of a specific path that guides superoxide toward an efficient reaction site. Localized positive charges at the rim of the entry channel are effective in increasing the reaction rate, whereas diffused positive charges or charges far from this area are not effective or are even detrimental, resulting in a misguided approach of the anion to the protein surface.

Electron self-exchange of cytochrome c measured via13C detected protonless NMR

The use of protonless 13C′–13C′ EXSY (COCO-EXSY) is proposed here to measure electron self-exchange rates. The experiment is compared to the commonly employed 1H and 15N EXSY experiments using as a reference system human cytochrome c. In COCO-EXSY, the exchange peaks are stronger than in the other experiments with respect to the self peaks and their intensity is less dependent on the choice of the EXSY mixing time. The use of 13C directed detection may be essential for all those cases where T2 relaxation is detrimental, as in the case of proteins containing highly paramagnetic metal centers, or rotating slowly in solution, or where the amide signals are difficult to detect due to chemical or conformational exchange. The proposed experiment has a general applicability and can be used to monitor exchange phenomena different from electron self-exchange.

Conformational properties of polyglutamine sequences in guanidine hydrochloride solutions

Two sets of iso-1-cytochrome c variants have been prepared with N-terminal insertions of pure polyglutamine, i.e., PolyQ variants, or polyglutamine interrupted with lysine every sixth residue, i.e., Gln-rich variants. The polymer properties of these pure polyGln or Gln-rich sequences have been evaluated using equilibrium and kinetic His-heme loop formation methods for loop sizes ranging from 22 to 46 in 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). In 6.0 M GdnHCl, the scaling exponent, ν(3), for the pure polyGln sequences, is ~1.7--significantly less than ν(3) ≈ 2.15 for the Gln-rich sequences. The stability of the His-heme loops becomes progressively greater for the pure polyGln sequences relative to the Gln-rich sequences as GdnHCl concentration decreases from 6.0 to 1.5 M. Thus, the context of the sequence effects the polymer properties of Gln repeats even in denaturing concentrations of GdnHCl. Comparison of data for the Gln-rich variants with previous results for Gly-rich and Ala-rich variants shows that ν(3) ~ 2.2 for the Gln-rich, Gly-rich, and Ala-rich sequences in 6.0 M GdnHCl, whereas ν(3) remains unchanged at 3.0 M GdnHCl concentration for the Gln-rich and Ala-rich sequences but decreases to ~1.7 for the Gly-rich sequences. Thus, the polymer properties of Gln-rich and Ala-rich sequences are less sensitive to solvent quality in denaturing solutions of GdnHCl than Gly-rich sequences. Evaluation of Flory's characteristic ratio, C(n), for the Gln-rich and Ala-rich sequences relative to the Gly-rich sequences shows that Gln-rich sequences are stiffer than Ala-rich sequences at both 3.0 and 6.0 M GdnHCl.

Flexibility of the metal-binding region in apo-cupredoxins

Protein-mediated electron transfer is an essential event in many biochemical processes. Efficient electron transfer requires the reorganization energy of the redox event to be minimized, which is ensured by the presence of rigid donor and acceptor sites. Electron transfer copper sites are present in the ubiquitous cupredoxin fold, able to bind one or two copper ions. The low reorganization energy in these metal centers has been accounted for by assuming that the protein scaffold creates an entatic/rack-induced state, which gives rise to a rigid environment by means of a preformed metal chelating site. However, this notion is incompatible with the need for an exposed metal-binding site and protein-protein interactions enabling metallochaperone-mediated assembly of the copper site. Here we report an NMR study that reveals a high degree of structural heterogeneity in the metal-binding region of the nonmetallated Cu(A)-binding cupredoxin domain, arising from microsecond to second dynamics that are quenched upon metal binding. We also report similar dynamic features in apo-azurin, a paradigmatic blue copper protein, suggesting a general behavior. These findings reveal that the entatic/rack-induced state, governing the features of the metal center in the copper-loaded protein, does not require a preformed metal-binding site. Instead, metal binding is a major contributor to the rigidity of electron transfer copper centers. These results reconcile the seemingly contradictory requirements of a rigid, occluded center for electron transfer, and an accessible, dynamic site required for in vivo copper uptake.

Redox-dependent conformational changes in eukaryotic cytochromes revealed by paramagnetic NMR spectroscopy

Cytochrome c (Cc) is a soluble electron carrier protein, transferring reducing equivalents between Cc reductase and Cc oxidase in eukaryotes. In this work, we assessed the structural differences between reduced and oxidized Cc in solution by paramagnetic NMR spectroscopy. First, we have obtained nearly-complete backbone NMR resonance assignments for iso-1-yeast Cc and horse Cc in both oxidation states. These were further used to derive pseudocontact shifts (PCSs) arising from the paramagnetic haem group. Then, an extensive dataset comprising over 450 measured PCSs and high-resolution X-ray and solution NMR structures of both proteins were used to define the anisotropic magnetic susceptibility tensor, Δχ. For most nuclei, the PCSs back-calculated from the Δχ tensor are in excellent agreement with the experimental PCS values. However, several contiguous stretches-clustered around G41, N52, and A81-exhibit large deviations both in yeast and horse Cc. This behaviour is indicative of redox-dependent structural changes, the extent of which is likely conserved in the protein family. We propose that the observed discrepancies arise from the changes in protein dynamics and discuss possible functional implications.

Cytochrome c folding / unfolding: a unifying picture

Recent NMR structural and dynamical data on partially folded forms of mono-heme cytochrome c provide a unifying picture of the behavior of the protein far from the native conditions and suggest useful hints to explain the redox dependent stability of the protein. A fragile hinge in the structure of mitochondrial cytochrome c is identified, which may not have correspondents in smaller type-1 cytochromes. Former spectroscopic and kinetic data are here discussed in terms of this new view.

Validation of a lanthanide tag for the analysis of protein dynamics by paramagnetic NMR spectroscopy

Paramagnetic lanthanide tags potentially can enhance the effects of microsecond to millisecond dynamics in proteins on NMR signals and provide structural information on lowly populated states encoded in the pseudocontact shifts. We have investigated the microsecond to millisecond mobility of a two-point attached lanthanide tag, CLaNP-5, using paramagnetic (1)H CPMG relaxation dispersion methods. CLaNP-5 loaded with Lu(3+), Yb(3+), or Tm(3+) was attached to three sites on the surface of two proteins, pseudoazurin and cytochrome c. The paramagnetic center causes large relaxation dispersion effects for two attachment sites, suggesting that local dynamics of the protein at the attachment site causes mobility of the paramagnetic center. At one site the relaxation dispersions are small and limited to the immediate environment of the tag. It is concluded that paramagnetic relaxation dispersion could represent a sensitive method to probe protein dynamics. However, the selection of a rigid attachment site is of critical importance.

Denatured states of low-complexity polypeptide sequences differ dramatically from those of foldable sequences

How the primary sequence of a protein encodes conformational preferences that operate early in folding to promote efficient formation of the correct native topology is still poorly understood. To address this issue, we have prepared a set of yeast iso-1-cytochrome c variants that contain polyalanine inserts ranging from 6 to 30 residues in length near the N terminus of the protein. We study the thermodynamics and kinetics of His-heme loop formation in the denatured state at 3 and 6 M guanidine-HCl concentration. We find that polyalanine closely approximates a random coil with excluded volume giving scaling exponents, nu(3), for equilibrium loop formation of 2.26 +/- 0.13 and 1.97 +/- 0.04 in 3 and 6 M guanidine-HCl, respectively. The rate of loop breakage initially decreases and then becomes independent of loop size as would be expected for a random coil. Comparison with previously reported data for denatured state His-heme loop formation for iso-1-cytochrome c and Rhodopseudomonas palustris cytochrome c', shows that foldable sequences deviate significantly from random coil behavior and that the deviation is fold-dependent.

Fast structural dynamics in reduced and oxidized cytochrome c

The sub-nanosecond structural dynamics of reduced and oxidized cytochrome c were characterized. Dynamic properties of the protein backbone measured by amide (15)N relaxation and side chains measured by the deuterium relaxation of methyl groups change little upon change in the redox state. These results imply that the solvent reorganization energy associated with electron transfer is small, consistent with previous theoretical analyses. The relative rigidity of both redox states also implies that dynamic relief of destructive electron transfer pathway interference is not operational in free cytochrome c.

Compressing the free energy range of substructure stabilities in iso-1-cytochrome c

Evolutionary conservation of substructure architecture between yeast iso-1-cytochrome c and the well-characterized horse cytochrome c is studied with limited proteolysis, the alkaline conformational transition and global unfolding with guanidine-HCl. Mass spectral analysis of limited proteolysis cleavage products for iso-1-cytochrome c show that its least stable substructure is the same as horse cytochrome c. The limited proteolysis data yield a free energy of 3.8 +/- 0.4 kcal mol(-1) to unfold the least stable substructure compared with 5.05 +/- 0.30 kcal mol(-1) for global unfolding of iso-1-cytochrome c. Thus, substructure stabilities of iso-1-cytochrome c span only approximately 1.2 kcal mol(-1) compared with approximately 8 kcal mol(-1) for horse cytochrome c. Consistent with the less cooperative folding thus expected for the horse protein, the guanidine-HCl m-values are approximately 3 kcal mol(-1)M(-1) versus approximately 4.5 kcal mol(-1)M(-1) for horse versus yeast cytochrome c. The tight free energy spacing of the yeast cytochrome c substructures suggests that its folding has more branch points than for horse cytochrome c. Studies on a variant of iso-1-cytochrome c with an H26N mutation indicate that the least and most stable substructures unfold sequentially and the two least stable substructures unfold independently as for horse cytochrome c. Thus, important aspects of the substructure architecture of horse cytochrome c, albeit compressed energetically, are preserved evolutionally in yeast iso-1-cytochrome c.

Submolecular unfolding units of Pseudomonas aeruginosa cytochrome c-551

Hydrogen exchange rates for backbone amide protons of oxidized Pseudomonas aeruginosa cytochrome c-551 (P. aeruginosa cytochrome c) have been measured in the presence of low concentrations of the denaturant guanidine hydrochloride. Analysis of the data has allowed identification of submolecular unfolding units known as foldons. The highest-energy foldon bears similarity to the proposed folding intermediate for P. aeruginosa cytochrome c. Parallels are seen to the foldons of the structurally homologous horse cytochrome c, although the heme axial methionine-bearing loop has greater local stability in P. aeruginosa cytochrome c, in accord with previous folding studies. Regions of low local stability are observed to correspond with regions that interact with redox partners, providing a link between foldon properties and function.

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.

Conformational stability and dynamics of cytochrome c affect its alkaline isomerization

The alkaline isomerization of horse heart ferricytochrome c (cyt c) has been studied by electronic absorption spectroscopy in the presence of the Hofmeister series of anions: chloride, bromide, rhodanide and perchlorate. The anions significantly affect the apparent pK (a) value of the transition in a concentration-dependent manner according to their position in the Hofmeister series. The Soret region of the absorption spectra is not affected by the presence of the salts and shows no significant structural perturbation of the heme crevice. In the presence of perchlorate and rhodanide anions, the cyanide exchange rate between the bulk solvent and the binding site is increased. These results imply higher flexibility of the protein structure in the presence of chaotropic salts. The thermal and isothermal denaturations monitored by differential scanning calorimetry and circular dichroism, respectively, showed a decrease in the conformational stability of cyt c in the presence of the chaotropic salts. A positive correlation between the stability, DeltaG, of cyt c and the apparent pK (a) values that characterize the alkaline transition indicates the presence of a thermodynamic linkage between these conformational transitions. In addition, the rate constant of the cyanide binding and the partial molar entropies of anions negatively correlate with the pK (a) values. This indicates the important role of anion-induced solvent reorganization on the structural flexibility of cyt c in the alkaline transitions.

Molecular statistics of cytochrome c: structural plasticity and molecular environment

Nuclear magnetic resonance experiments performed on yeast mitochondrial cytochrome c (Cytc), a paradigmatic electron transfer protein, reveal that the two oxidation states have similar structures, but different mobility: despite the few structural differences compared with the reduced form, the oxidized form displays a larger unfolding propensity. Molecular dynamics simulations performed on both NMR reduced and NMR oxidized forms show that the reduced form has a larger solvent-accessible surface area (SASA). Starting from this observation, a molecular statistical approach was then applied in order to correlate the molecular surface to molecular mobility. Simulations started from biased initial conditions corresponding to different molecular sizes were combined with the maximal constrained entropy method. The NMR structure of oxidized Cytc is more suited to expose a smaller SASA than the NMR structure of the reduced form, but the accessible conformational landscape at 300 K around the NMR oxidized structure is flatter than for the NMR reduced structure. Protein configurations of smaller SASA and size display larger plasticity when they resemble the NMR oxidized structure, whereas they are more rigid when they resemble the NMR reduced structure. Implications of the results for the protein properties during its functional process are discussed.

Redox-Coupled Dynamics and Folding in Cytochrome c

Cytochrome c functions as an electron carrier in the mitochondrial electron-transport chain using the Fe(II)-Fe(III) redox couple of a covalently attached heme prosthetic group, and it has served as a paradigm for both biological redox activity and protein folding. On the basis of a wide variety of biophysical techniques, it has been suggested that the protein is more flexible in the oxidized state than in the reduced state, which has led to speculation that it is the dynamics of the protein that has been evolved to control the cofactor's redox properties. To test this hypothesis, we incorporated carbon-deuterium bonds throughout cytochrome c and characterized their absorption frequencies and line widths using IR spectroscopy. The absorption frequencies of several residues on the proximal side of the heme show redox-dependent changes, but none show changes in line width, implying that the flexibility of the oxidized and reduced proteins is not different. However, the spectra demonstrate that folded protein is in equilibrium with a surprisingly large amount of locally unfolded protein, which increases with oxidation for residues localized to the proximal side of the heme. The data suggest that while the oxidized protein is not more flexible than the reduced protein, it is more locally unfolded. Local unfolding of cytochrome c might be one mechanism whereby the protein evolved to control electron transfer.

The Effects of Ligand Exchange and Mobility on the Peroxidase Activity of a Bacterial Cytochrome c upon Unfolding

The effect on the heme environment upon unfolding Paracoccus versutus ferricytochrome c-550 and two site-directed variants, K99E and H118Q, has been assessed through a combination of peroxidase activity increase and one-dimensional NMR spectroscopy. At pH 4.5, the data are consistent with a low- to high-spin heme transition, with the K99E mutation resulting in a protein with increased peroxidase activity in the absence of or at low concentrations of denaturant. Furthermore, the mobility of the polypeptide chain at pH 4.5 for the wild-type protein has been monitored in the absence and presence of denaturant through heteronuclear NMR experiments. The results are discussed in terms of local stability differences between bacterial and mitochondrial cytochromes c that are inferred from peroxidase activity assays. At pH 7.0, a mixture of misligated heme states arising from protein-based ligands assigned to lysine and histidine is detected. At low denaturant concentrations, these partially unfolded misligated heme forms inhibit the peroxidase activity. Data from the K99E mutation at pH 7.0 indicate that K99 is not involved in heme misligation, whereas histidine coordination is proven by the data from the H118Q variant.

Folding, Conformational Changes, and Dynamics of Cytochromes c Probed by NMR Spectroscopy

NMR spectroscopy has become a vital tool for studies of protein conformational changes and dynamics. Oxidized Fe(III)cytochromes c are a particularly attractive target for NMR analysis because their paramagnetism (S = (1)/(2)) leads to high (1)H chemical shift dispersion, even for unfolded or otherwise disordered states. In addition, analysis of shifts induced by the hyperfine interaction reveals details of the structure of the heme and its ligands for native and nonnative protein conformational states. The use of NMR spectroscopy to investigate the folding and dynamics of paramagnetic cytochromes c is reviewed here. Studies of nonnative conformations formed by denaturation and by anomalous in vivo maturation (heme attachment) are facilitated by the paramagnetic, low-spin nature of native and nonnative forms of cytochromes c. Investigation of the dynamics of folded cytochromes c also are aided by their paramagnetism. As an example of this analysis, the expression in Escherichia coli of cytochrome c(552) from Nitrosomonas europaea is reported here, along with analysis of its unusual heme hyperfine shifts. The results are suggestive of heme axial methionine fluxion in N. europaea ferricytochrome c(552). The application of NMR spectroscopy to investigate paramagnetic cytochrome c folding and dynamics has advanced our understanding of the structure and dynamics of both native and nonnative states of heme proteins.

Heme axial methionine fluxionality in Hydrogenobacter thermophilus cytochrome c552

The heme group in paramagnetic (S = 1/2) ferricytochromes c typically displays a markedly asymmetric distribution of unpaired electron spin density among the heme pyrrole beta substituents. This asymmetry is determined by the orientations of the heme axial ligands, histidine and methionine. One exception to this is ferricytochrome c(552) from Hydrogenobacter thermophilus, which has similar amounts of unpaired electron spin density at the beta substituents on all four heme pyrroles. Here, determination of the orientation of the magnetic axes and analysis of NMR line shapes for H. thermophilus ferricytochrome c(552) is performed. These data reveal that the unusual electronic structure for this protein is a result of fluxionality of the heme axial methionine. It is proposed that the ligand undergoes inversion at the pyramidal sulfur, and the rapid interconversion between two diastereomeric forms results in the unusual heme electronic structure. Thus a fluxional process for a metal-bound amino acid side chain has now been identified.

Cytochrome c and SDS: a molten globule protein with altered axial ligation

Saccharomices cerevisiae (yeast iso-1) cytochrome c has been investigated in the presence of 100 mM SDS in order to simulate the interaction of cytochrome c with membrane. Under these circumstances, a high spin species with detached methionine axial ligand is observed through NMR, in analogy to findings on the horse heart protein. However, at variance with the latter system, for the yeast protein also a low spin species is detected, which appears to be present with a concentration of about 40% with respect to that of the high spin species. The R(1), R(2), [1H]-15N NOE of backbone amides which are not affected by paramagnetism are homogeneous and allow a simultaneous analysis of the data for the two species. The result is that the rotational correlation time is larger than in water and larger than expected on the basis of viscosity of the SDS-containing solution. This finding suggests interactions of cytochrome c with SDS. Furthermore, it appears that there is subnanosecond backbone mobility, which also accounts for the decreased intensity of NOE cross-peaks and may be associated with equilibria between helical and random coil structure. The dynamic behavior appears to be a common feature of the high spin and low spin species and is consistent with the presence of a molten globule state. The molten globule nature of the protein could account for the presence of the different axial coordination of the heme iron. Such findings are meaningful with respect to the physiology of cytochrome c as electron transfer protein and as promoter of apoptosis.

Cooperative omega loops in cytochrome c: role in folding and function

Hydrogen exchange experiments under slow exchange conditions show that an omega loop in cytochrome c (residues 40-57) acts as a cooperative unfolding/refolding unit under native conditions. This unit behavior accounts for an initial step on the unfolding pathway, a final step in refolding, and a number of other structural, functional and evolutionary properties.

15N-1H Residual dipolar coupling analysis of native and alkaline-K79A Saccharomyces cerevisiae cytochrome c

Residual dipolar couplings (RDCs) and pseudocontact shifts are experimentally accessible properties in nuclear magnetic resonance that are related to structural parameters and to the magnetic susceptibility anisotropy. We have determined RDCs due to field-induced orientation of oxidized-K79A and reduced cytochrome c at pH 7.0 and oxidized-K79A cytochrome c at pH 11.1 through measurements of amide (15)N-(1)H (1)J couplings at 800 and 500 MHz. The pH 7.0 RDCs for Fe(III)- and Fe(II)-cytochrome c together with available nuclear Overhauser effects were used to recalculate solution structures that were consistent with both sets of constraints. Molecular magnetic susceptibility anisotropy values were calculated for both redox states of the protein. By subtracting the residual dipolar couplings (RDCs) of the reduced form from those of the oxidized form measured at the same magnetic field (800 MHz), we found the RDC contribution of the paramagnetic metal ion in the oxidized protein. The magnetic susceptibility anisotropy, which was calculated from the structure, was found to be the same as that of the paramagnetic metal ion obtained independently from pseudocontact shifts, thereby indicating that the elements of secondary structure either are rigid or display the same mobility in both oxidation states. The residual dipolar coupling values of the alkaline-K79A form are small with respect to those of oxidized native cytochrome, whereas the pseudocontact shifts are essentially of the same magnitude, indicating local mobility. Importantly, this is the first time that mobility has been found through comparison of RDCs with pseudocontact shifts.

Structural model for an alkaline form of ferricytochrome C

An (15)N-enriched sample of the yeast iso-1-ferricytochrome c triple variant (Lys72Ala/Lys79Ala/Cys102Thr) in an alkaline conformation was examined by NMR spectroscopy. The mutations were planned to produce a cytochrome c with a single conformer. Despite suboptimal conditions for the collection of spectra (i.e., pH approximately equal to 11), NMR remains a suitable investigation technique capable of taking advantage of paramagnetism. 76% of amino acids and 49% of protons were assigned successfully. The assignment was in part achieved through standard methods, in part through the identification of groups maintaining the same conformation as in the native protein at pH 7 and, for a few other residues, through a tentative analysis of internuclear distance predictions. Lys73 was assigned as the axial ligand together with His18. In this manner, 838 meaningful NOEs for 108 amino acids, 50 backbone angle constraints, and 203 pseudocontact shifts permitted the convergence of randomly generated structures to a family of conformers with a backbone RMSD of 1.5 +/- 0.2 A. Most of the native cytochrome c conformation is maintained at high pH. The NOE pattern that involves His18 clearly indicates that the proximal side of the protein, including the 20s and 40s loops, remains essentially intact. Structural differences are concentrated in the 70-80 loop, because of the replacement of Met80 by Lys73 as an axial ligand, and in the 50s helix facing that loop; as a consequence, there is increased exposure of the heme group to solvent. Based on several spectral features, we conclude that the folded polypeptide is highly fluxional.

NMR study of the conformational transition of cytochrome c upon the displacement of Met80 by exogenous ligand: structural and magnetic characterization of azidoferricytochrome c

As the exogenous ligand-cytochrome c complexes were purported to represent models for the unfolding intermediate of cytochrome c, NMR spectroscopy has been utilized to study the azide adduct of horse heart cytochrome c. The structure of azidoferricytochrome c was modeled by restrained energy minimization using paramagenetic pseudocontact shifts as constraints. The bound azide moiety was found to be tilted approximately 15 degrees from the heme normal. The displacement of Met80 by the exogenous azide molecule causes large structural rearrangement in the distal cavity. Furthermore, the conformation transition associated with the swing out of the loop containing Met80 and the shift of the 50s-helix increases the solvent accessibility of the heme group. To elucidate the heme electronic structure of the complex, the paramagnetic 13C shifts were analyzed in terms of a model based on the pi molecular orbitals of the heme under perturbed D(4) symmetry. It turned out that the His-Fe bonding provides the protein constraint that orients the in-plane anisotropy in the complex. The electronic properties are in accordance with the calculated magnetic susceptibility anisotropy and the structural information.

The unfolding of oxidized c-type cytochromes: the instructive case of Bacillus pasteurii

The reversible unfolding of oxidized Bacillus pasteurii cytochrome c(553) by guanidinium chloride under equilibrium conditions has been monitored by NMR and optical spectroscopy. The results obtained indicate that unfolding takes place through a mechanism involving the detachment from heme iron coordination of the sulfur of the Met71 axial ligand and yielding either a high spin (HS) or a low spin (LS(1)) species, depending on the pH value. In the LS(1) form the Met71 is replaced by another protein ligand, possibly Lys. The ligand exchange reaction does not reach completion until the protein backbone reaches a largely unfolded state, as monitored through 1H-15N NMR experiments, thus demonstrating that there is a significant correlation between formation of the Fe-S bond and native structure stability. 1H/2H exchange data, however, show that helix alpha(3), the C-terminal region of helix alpha(4), and helix alpha(5) maintain low exchangeability of the amide protons in the LS(1) form. This finding most likely implies that these regions maintain some ordered non-covalent structure, in which the amide moieties are involved in H-bonds. Finally, a folding mechanism is proposed and discussed in terms of analogies and differences with the larger mitochondrial cytochrome c proteins. It is concluded that the thermodynamic stability of the region around the metal cofactor is determined by the chemical nature of the residues around the axial methionine residue.

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

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

NMR solution structure, backbone mobility, and homology modeling of c-type cytochromes from gram-positive bacteria

The solution structure of oxidized cytochrome c(553) (71 amino acid residues) from the Gram-positive bacterium Bacillus pasteurii is here reported and compared with the available crystal structure. The solution structure is obtained from 1609 meaningful NOE data (22.7 per residue), 76 dihedral angles, and 59 pseudocontact shifts. The root mean square deviations from the average structure are 0.25+/-0.07 and 0.59+/-0.13 A for the backbone and all heavy atoms, respectively, and the quality assessment of the structure is satisfactory. The solution structure closely reproduces the fold observed in the crystal structure. The backbone mobility was then investigated through amide (15)N relaxation rate and (15)N-(1)H NOE measurements. The protein is rigid in both the sub-nanosecond and millisecond time scales, probably due to the relatively large heme:number of amino acids ratio. Modeling of eight c-type cytochromes from other Gram-positive bacteria with a high sequence identity (>30 %) to the present cytochrome c(553) was performed. Analysis of consensus features accounts for the relatively low reduction potential as being due to extensive heme hydration and indicates residues 34-35, 44-46, 69-72, and 75 as a conserved hydrophobic patch for the interaction with a protein partner. At variance with mitochondrial c-type cytochrome, this protein does not experience pH-dependent coordination equilibria. The reasons for this difference are analyzed.

Backbone dynamics of plastocyanin in both oxidation states. Solution structure of the reduced form and comparison with the oxidized state

A model-free analysis based on (15)N R(1), (15)N R(2), and (15)N-(1)H nuclear Overhauser effects was performed on reduced (diamagnetic) and oxidized (paramagnetic) forms of plastocyanin from Synechocystis sp. PCC6803. The protein backbone is rigid, displaying a small degree of mobility in the sub-nanosecond time scale. The loops surrounding the copper ion, involved in physiological electron transfer, feature a higher extent of flexibility in the longer time scale in both redox states, as measured from D(2)O exchange of amide protons and from NH-H(2)O saturation transfer experiments. In contrast to the situation for other electron transfer proteins, no significant difference in the dynamic properties is found between the two redox forms. A solution structure was also determined for the reduced plastocyanin and compared with the solution structure of the oxidized form in order to assess possible structural changes related to the copper ion redox state. Within the attained resolution, the structure of the reduced plastocyanin is indistinguishable from that of the oxidized form, even though small chemical shift differences are observed. The present characterization provides information on both the structural and dynamic behavior of blue copper proteins in solution that is useful to understand further the role(s) of protein dynamics in electron transfer processes.

A further clue to understanding the mobility of mitochondrial yeast cytochrome c: a (15)N T1rho investigation of the oxidized and reduced species

A new approach was developed to overproduce 15N-enriched yeast iso-1-cytochrome c in the periplasm of Escherichia coli in order to perform a study of the motions in the ms-micros time scale on the oxidized and reduced forms through rotating frame 15N relaxation rates and proton/deuterium exchange studies. It is confirmed that the reduced protein is rather rigid whereas the oxidized species is more flexible. The regions of the protein that display increased internal mobility upon oxidation are easily identified by the number of residues experiencing conformational equilibria and by their exchange rates. These data complement the information already available in the literature and provide a comprehensive picture of the mobility in the protein. In particular, oxidation mobilizes the loop containing Met80 and, through specific contacts, affects the mobility of helix 3 and possibly of helix 5, and of a section of protein connecting the heme propionates to helix 2. The relevance of internal motions to molecular recognition and to the early steps of the unfolding process of the oxidized species is also discussed. In agreement with the reported data, subnanosecond mobility is found to be less informative than the ms-micros with respect to redox dependent properties.

Solution structure of Methylophilus methylotrophus cytochrome c": insights into the structural basis of haem-ligand detachment

Cytochrome c" from Methylophilus methylotrophus is a monohaem protein with 124 amino acid residues. The iron has two histidine ligands in the oxidised form, one of which detaches and picks up a proton when the protein is reduced. Thus, both forms are paramagnetic. The structure of the oxidised form in solution, determined from NMR data is presented. The family of structures has an average backbone rmsd value of 0.53 A, and a heavy atom rmsd value of 0.95 A, within a target function range of 32 %. This structure is related to class I cytochromes with an additional helix at the N terminus. The haem-binding site occurs in a domain essentially lacking secondary structure motifs and the axial histidinyl residues were found in an unusual near perpendicular orientation. Moreover, a disulfide bridge is present, an uncommon structural feature among c-type cytochromes. The disulfide bridge, linking cysteine residues 96 and 104, forms a loop that confers rigidity and is essential to the detachment of the axial histidine (His95) as demonstrated by chemical disruption of the S-S bond. A route for protonation of the distal histidine involving haem propionate 17 is proposed and discussed in the light of available models for complex membrane proton pumps.

Heme ligation and conformational plasticity in the isolated c domain of cytochrome cd1 nitrite reductase

The heme ligation in the isolated c domain of Paracoccus pantotrophus cytochrome cd(1) nitrite reductase has been characterized in both oxidation states in solution by NMR spectroscopy. In the reduced form, the heme ligands are His69-Met106, and the tertiary structure around the c heme is similar to that found in reduced crystals of intact cytochrome cd1 nitrite reductase. In the oxidized state, however, the structure of the isolated c domain is different from the structure seen in oxidized crystals of intact cytochrome cd1, where the c heme ligands are His69-His17. An equilibrium mixture of heme ligands is present in isolated oxidized c domain. Two-dimensional exchange NMR spectroscopy shows that the dominant species has His69-Met106 ligation, similar to reduced c domains. This form is in equilibrium with a high-spin form in which Met106 has left the heme iron. Melting studies show that the midpoint of unfolding of the isolated c domain is 320.9 +/- 1.2 K in the oxidized and 357.7 +/- 0.6 K in the reduced form. The thermally denatured forms are high-spin in both oxidation states. The results reveal how redox changes modulate conformational plasticity around the c heme and show the first key steps in the mechanism that lead to ligand switching in the holoenzyme. This process is not solely a function of the properties of the c domain. The role of the d1 heme in guiding His17 to the c heme in the oxidized holoenzyme is discussed.