Kinetic mechanism of cytochrome c folding: involvement of the heme and its ligands

Probing the versatility of cytochrome c by spectroscopic means: A Laudatio on resonance Raman spectroscopy

Over the last 50 years resonance Raman spectroscopy has become an invaluable tool for the exploration of chromophores in biological macromolecules. Among them, heme proteins and metal complexes have attracted considerable attention. This interest results from the fact that resonance Raman spectroscopy probes the vibrational dynamics of these chromophores without direct interference from the surrounding. However, the indirect influence via through-bond and through-space chromophore-protein interactions can be conveniently probed and analyzed. This review article illustrates this point by focusing on class 1 cytochrome c, a comparatively simple heme protein generally known as electron carrier in mitochondria. The article demonstrates how through selective excitation of resonance Raman active modes information about the ligation, the redox state and the spin state of the heme iron can be obtained from band positions in the Raman spectra. The investigation of intensities and depolarization ratios emerged as tools for the analysis of in-plane and out-of-plane deformations of the heme macrocycle. The article further shows how resonance Raman spectroscopy was used to characterize partially unfolded states of oxidized cytochrome c. Finally, it describes its use for exploring structural changes due to the protein's binding to anionic surfaces like cardiolipin containing membranes.

A surprisingly simple three-state generic process for reversible protein denaturation by trifluoroethanol

Despite the rich knowledge of the influence of 2,2,2-trifluoroethanol (TFE) on the structure and conformation of peptides and proteins, the mode(s) of TFE-protein interactions and the mechanism by which TFE reversibly denatures a globular protein remain elusive. This study systematically examines TFE-induced equilibrium transition curves for six paradigmatic globular proteins by using basic fluorescence and circular dichroism measurements under neutral pH conditions. The results are remarkably simple. Low TFE invariably unfolds the tertiary structure of all proteins to produce the obligate intermediate (I) which retains nearly all of native-state secondary structure, but enables the formation of extra α-helices as the level of TFE is raised higher. Inspection of the transitions at once reveals that the tertiary structure unfolding is always a distinct process, necessitating the inclusion of at least one obligate intermediate in the TFE-induced protein denaturation. It appears that the intermediate in the minimal unfolding mechanism N⇌I⇌D somehow acquires higher α-helical propensity to generate α-helices in excess of that in the native state to produce the denatured state (D), also called the TFE state. The low TFE-populated intermediate I may be called a universal intermediate by virtue of its α-helical propensity. Contrary to many earlier suggestions, this study dismisses molten globule (MG)-like attribute of I or D.

Exploring the folding energy landscapes of heme proteins using a hybrid AWSEM-heme model

Heme is an active center in many proteins. Here we explore computationally the role of heme in protein folding and protein structure. We model heme proteins using a hybrid model employing the AWSEM Hamiltonian, a coarse-grained forcefield for the protein chain along with AMBER, an all-atom forcefield for the heme. We carefully designed transferable force fields that model the interactions between the protein and the heme. The types of protein-ligand interactions in the hybrid model include thioester covalent bonds, coordinated covalent bonds, hydrogen bonds, and electrostatics. We explore the influence of different types of hemes (heme b and heme c) on folding and structure prediction. Including both types of heme improves the quality of protein structure predictions. The free energy landscape shows that both types of heme can act as nucleation sites for protein folding and stabilize the protein folded state. In binding the heme, coordinated covalent bonds and thioester covalent bonds for heme c drive the heme toward the native pocket. The electrostatics also facilitates the search for the binding site.

Open-Bundle Structure as the Unfolding Intermediate of Cytochrome c' Revealed by Small Angle Neutron Scattering

The dynamic structure changes, including the unfolding, dimerization, and transition from the compact to the open-bundle unfolding intermediate structure of Cyt c', were detected by a small-angle neutron scattering experiment (SANS). The structure of Cyt c' was changed into an unstructured random coil at pD = 1.7 (Rg = 25 Å for the Cyt c' monomer). The four-α-helix bundle structure of Cyt c' at neutral pH was transitioned to an open-bundle structure (at pD ~13), which is given by a numerical partial scattering function analysis as a joint-clubs model consisting of four clubs (α-helices) connected by short loops. The compactly folded structure of Cyt c' (radius of gyration, Rg = 18 Å for the Cyt c' dimer) at neutral or mildly alkaline pD transited to a remarkably larger open-bundle structure at pD ~13 (Rg = 25 Å for the Cyt c' monomer). The open-bundle structure was also supported by ab initio modeling.

New insights into the folding–unfolding mechanism and conformations of cytochrome C

Optical trapping experiments offer new insights into the folding and unfolding of cytochrome C.

Use of 3D domain swapping in constructing supramolecular metalloproteins

Supramolecules, which are formed by assembling multiple molecules by noncovalent intermolecular interactions instead of covalent bonds, often show additional properties that cannot be exhibited by a single molecule. Supramolecules have evolved into molecular machines in the field of chemistry, and various supramolecular proteins are responsible for life activities in the field of biology. The design and creation of supramolecular proteins will lead to development of new enzymes, functional biomaterials, drug delivery systems, etc.; thus, the number of studies on the regulation of supramolecular proteins is increasing year by year. Several methods, including disulfide bond, metal coordination, and surface-surface interaction, have been utilized to construct supramolecular proteins. In nature, proteins have been shown to form oligomers by 3D domain swapping (3D-DS), a phenomenon in which a structural region is exchanged between molecules of the same protein. We have been studying the mechanism of 3D-DS and utilizing 3D-DS to construct supramolecular metalloproteins. Cytochrome c forms cyclic oligomers and polymers by 3D-DS, whereas other metalloproteins, such as various c-type cytochromes and azurin form small oligomers and myoglobin forms a compact dimer. We have also utilized 3D-DS to construct heterodimers with different active sites, a protein nanocage encapsulating a Zn-SO₄ cluster in the internal cavity, and a tetrahedron with a designed building block protein. Protein oligomer formation was controlled for the 3D-DS dimer of a dimer-monomer transition protein. This article reviews our research on supramolecular metalloproteins.

Dissecting Multiple Pathways in the Relaxation Dynamics of Helix <==> Coil Transitions with Optimum Dimensionality Reduction

Fast kinetic experiments with dramatically improved time resolution have contributed significantly to understanding the fundamental processes in protein folding pathways involving the formation of a-helices and b-hairpin, contact formation, and overall collapse of the peptide chain. Interpretation of experimental results through application of a simple statistical mechanical model was key to this understanding. Atomistic description of all events observed in the experimental findings was challenging. Recent advancements in theory, more sophisticated algorithms, and a true long-term trajectory made way for an atomically detailed description of kinetics, examining folding pathways, validating experimental results, and reporting new findings for a wide range of molecular processes in biophysical chemistry. This review describes how optimum dimensionality reduction theory can construct a simplified coarse-grained model with low dimensionality involving a kinetic matrix that captures novel insights into folding pathways. A set of metastable states derived from molecular dynamics analysis generate an optimally reduced dimensionality rate matrix following transition pathway analysis. Analysis of the actual long-term simulation trajectory extracts a relaxation time directly comparable to the experimental results and confirms the validity of the combined approach. The application of the theory is discussed and illustrated using several examples of helix <==> coil transition pathways. This paper focuses primarily on a combined approach of time-resolved experiments and long-term molecular dynamics simulation from our ongoing work.

Real-time nuclear magnetic resonance spectroscopy in the study of biomolecular kinetics and dynamics

The review describes the application of nuclear magnetic resonance (NMR) spectroscopy to study kinetics of folding, refolding and aggregation of proteins, RNA and DNA. Time-resolved NMR experiments can be conducted in a reversible or an irreversible manner. In particular, irreversible folding experiments pose large requirements for (i) signal-to-noise due to the time limitations and (ii) synchronising of the refolding steps. Thus, this contribution discusses the application of methods for signal-to-noise increases, including dynamic nuclear polarisation, hyperpolarisation and photo-CIDNP for the study of time-resolved NMR studies. Further, methods are reviewed ranging from pressure and temperature jump, light induction to rapid mixing to induce rapidly non-equilibrium conditions required to initiate folding.

Electrical unfolding of cytochrome c during translocation through a nanopore constriction

Many small proteins move across cellular compartments through narrow pores. In order to thread a protein through a constriction, free energy must be overcome to either deform or completely unfold the protein. In principle, the diameter of the pore, along with the effective driving force for unfolding the protein, as well as its barrier to translocation, should be critical factors that govern whether the process proceeds via squeezing, unfolding/threading, or both. To probe this for a well-established protein system, we studied the electric-field-driven translocation behavior of cytochrome c (cyt c) through ultrathin silicon nitride (SiN_x) solid-state nanopores of diameters ranging from 1.5 to 5.5 nm. For a 2.5-nm-diameter pore, we find that, in a threshold electric-field regime of ∼30 to 100 MV/m, cyt c is able to squeeze through the pore. As electric fields inside the pore are increased, the unfolded state of cyt c is thermodynamically stabilized, facilitating its translocation. In contrast, for 1.5- and 2.0-nm-diameter pores, translocation occurs only by threading of the fully unfolded protein after it transitions through a higher energy unfolding intermediate state at the mouth of the pore. The relative energies between the metastable, intermediate, and unfolded protein states are extracted using a simple thermodynamic model that is dictated by the relatively slow (∼ms) protein translocation times for passing through the nanopore. These experiments map the various modes of protein translocation through a constriction, which opens avenues for exploring protein folding structures, internal contacts, and electric-field-induced deformability.

The effect of cardiolipin side chain composition on cytochrome c protein conformation and peroxidase activity

Skeletal muscle, a highly active tissue, makes up 40% of the total body weight. This tissue relies on mitochondria for ATP production, calcium homeostasis, and programed cell death. Mitochondrial phospholipid composition, namely, cardiolipin (CL), influences the functional efficiency of mitochondrial proteins, specifically cytochrome c. The interaction of CL with cytochrome c in the presence of free radicals induces structural and functional changes promoting peroxidase activity and cytochrome c release, a key event in the initiation of apoptosis. The CL acyl chain degree of saturation has been implicated in the cytochrome c to cytochrome c peroxidase transition in liposomal models. However, mitochondrial membranes are composed of differing CL acyl chain composition. Currently, it is unclear how differing CL acyl chain composition utilizing liposomes will influence the cytochrome c form and function as a peroxidase. Thus, this study examined the role of CL acyl chain saturation within liposomes broadly reflecting the relative CL composition of mitochondrial membranes from healthy and dystrophic mouse muscle on cytochrome c conformation and function. Despite no differences in protein conformation or function between healthy and dystrophic liposomes, cytochrome c's affinity to CL increased with greater unsaturation. These findings suggest that increasing CL acyl chain saturation, as implicated in muscle wasting diseases, may not influence cytochrome c transformation and function as a peroxidase but may alter its interaction with CL, potentially impacting further downstream effects.

pH-Induced Switch between Different Modes of Cytochrome c Binding to Cardiolipin-Containing Liposomes

Fluorescence, visible circular dichroism (CD), absorption, and resonance Raman spectroscopy techniques were combined to explore structural changes of ferricytochrome c upon its binding to cardiolipin-containing liposomes (20% 1,1',1,2'-tetraoleyolcardiolipin and 1,2-deoleyol-sn-glycero-3-phosphocholine) at acidic pH (6.5). According to the earlier work of Kawai [J. Biol. Chem.2005, 280, 34709-347171],cytochrome c binding at this pH is governed by interactions between the phosphate head groups of cardiolipin and amino acid side chains of the so-called L-site, which contains the charged residues K22, K25, K27, and potentially H26 and H33. We found that L-site binding causes a conformational transition that involves a change of the protein's ligation and spin state. In this paper, we report spectroscopic responses to an increasing number of cardiolipin-containing liposomes at pH 6.5 in the absence and presence of NaCl. The latter was found to mostly inhibit protein binding already with 50 mM concentration. The inhibition effect can be quantitatively reproduced by applying the electrostatic theory of Heimburg [Biophys. J.1995, 68, 536-546]. A comparison with corresponding spectroscopic response data obtained at pH 7.4 reveals major differences in that the latter indicates hydrophobic binding, followed by an electrostatically driven conformational change. Visible CD data suggest that structural changes in the heme pocket of liposome-bound ferricytochrome c resemble to some extent those in the denatured protein in urea at neutral and acidic pH. The measured noncoincidence between absorption and CD Soret band of cytochrome c in the presence of a large access of cardiolipin is caused by the electric field at the membrane surface. The very fact that its contribution to the internal electric field in the heme pocket is detectable by spectroscopic means suggests some penetration of the protein into membrane surface.

Uncovering dehydration in cytochrome c refolding from urea- and guanidine hydrochloride-denatured unfolded state by high pressure spectroscopy

To investigate the dehydration associated with protein folding, the partial molar volume changes for protein unfolding (ΔV_u) in cytochrome c (Cyt c) were determined using high pressure absorption spectroscopy. ΔV_u values for the unfolding to urea- and guanidine hydrochloride (GdnHCl)-denatured Cyt c were estimated to be 56±5 and 29±1 mL mol⁻¹, respectively. Considering that the volume change for hydration of hydrophobic groups is positive and that Cyt c has a covalently bonded heme, a positive ΔV_u reflects the primary contribution of the hydration of heme. Because of the marked tendency of guanidium ions to interact with hydrophobic groups, a smaller number of water molecules were hydrated with hydrophobic groups in GdnHCl-denatured Cyt c than in urea-denatured Cyt c, resulting in the smaller positive ΔV_u. On the other hand, urea is a relatively weak denaturant and urea-denatured Cyt c is not completely hydrated, which retains the partially folded structures. To unfold such partial structures, we introduced a mutation near the heme binding site, His26, to Gln, resulting in a negatively shifted ΔV_u (4±2 mL mol⁻¹) in urea-denatured Cyt c. The formation of the more solvated and less structured state in the urea-denatured mutant enhanced hydration to the hydrophilic groups in the unfolding process. Therefore, we confirmed the hydration of amino acid residues in the protein unfolding of Cyt c by estimating ΔV_u, which allows us to discuss the hydrated structures in the denatured states of proteins.

Relating the multi-functionality of cytochrome c to membrane binding and structural conversion

Cytochrome c is known as an electron-carrying protein in the respiratory chain of mitochondria. Over the last 20 years, however, alternative functions of this very versatile protein have become the focus of research interests. Upon binding to anionic lipids such as cardiolipin, the protein acquires peroxidase activity. Multiple lines of evidence suggest that this requires a conformational change of the protein which involves partial unfolding of its tertiary structure. This review summarizes the current state of knowledge of how cytochrome c interacts with cardiolipin-containing surfaces and how this affects its structure and function. In this context, we delineate partially conflicting results regarding the affinity of cytochrome c binding to cardiolipin-containing liposomes of different size and its influence on the structure of the protein and the morphology of the membrane.

Probing Cytochrome c Folding Transitions upon Phototriggered Environmental Perturbations Using Time-Resolved X-ray Scattering

Direct tracking of protein structural dynamics during folding-unfolding processes is important for understanding the roles of hierarchic structural factors in the formation of functional proteins. Using cytochrome c (cyt c) as a platform, we investigated its structural dynamics during folding processes triggered by local environmental changes (i.e., pH or heme iron center oxidation/spin/ligation states) with time-resolved X-ray solution scattering measurements. Starting from partially unfolded cyt c, a sudden pH drop initiated by light excitation of a photoacid caused a structural contraction in microseconds, followed by active site restructuring and unfolding in milliseconds. In contrast, the reduction of iron in the heme via photoinduced electron transfer did not affect conformational stability at short timescales (<1 ms), despite active site coordination geometry changes. These results demonstrate how different environmental perturbations can change the nature of interaction between the active site and protein conformation, even within the same metalloprotein, which will subsequently affect the folding structural dynamics.

Remote Perturbations in Tertiary Contacts Trigger Ligation of Lysine to the Heme Iron in Cytochrome c

Perturbations in protein structure define the mechanism of allosteric regulation and biological information transfer. In cytochrome c (cyt c), ligation of Met80 to the heme iron is critical for the protein's electron-transfer (ET) function in oxidative phosphorylation and for suppressing its peroxidase activity in apoptosis. The hard base Lys is a better match for the hard ferric iron than the soft base Met is, suggesting the key role of the protein scaffold in favoring Met ligation. To probe the role of the protein structure in the maintenance of Met ligation, mutations T49V and Y67R/M80A were designed to disrupt hydrogen bonding and packing of the heme coordination loop, respectively. Electronic absorption, nuclear magnetic resonance, and electron paramagnetic resonance spectra reveal that ferric forms of both variants are Lys-ligated at neutral pH. A minor change in the tertiary contacts in T49V, away from the heme coordination loop, appears to be sufficient to execute a change in ligation, suggesting a cross-talk between the different regions of the protein structure and a possibility of built-in conformational switches in cyt c. Analyses of thermodynamic stability, kinetics of Lys binding and dissociation, and the pH-dependent changes in ligation provide a detailed characterization of the Lys coordination in these variants and relate these properties to the extent of structural perturbations. The findings emphasize the importance of the hydrogen-bonding network in controlling ligation of the native Met80 to the heme iron.

pH Dependence of Ferricytochrome c Conformational Transitions during Binding to Cardiolipin Membranes: Evidence for Histidine as the Distal Ligand at Neutral pH

The conformational changes of ferricytochrome c upon binding to cardiolipin-containing small unilamellar vesicles were studied at slightly acidic pH using fluorescence, visible circular dichroism, UV-visible absorption, and resonance Raman spectroscopy. The obtained spectroscopic response data suggest a mode of interaction, which is clearly distinct from the binding process observed at neutral pH. Evidence of a reversible and electrostatic binding mechanism under these conditions is provided through binding inhibition in the presence of 150 mM NaCl. Moreover, UV-visible absorption and resonance Raman spectra reveal that the conformational ensemble of membrane bound cytochrome c is dominated by a mixture of conformers with pentacoordinated and hexacoordinated high-spin heme irons, which contrast with the dominance of low-spin species at neutral pH. While our results confirm the L-site binding proposed by Kawai et al., they point to the protonation of a histidine ligand (H33) as conformational trigger.

Time resolved calorimetry of photo-induced folding in horse heart cytochrome c at high pH

Here the molar volume and enthalpy changes associated with the early events in the folding of ferrocytochrome c (Cc) at high pH have been examined using time resolved photoacoustic calorimetry (PAC). The data reveal an overall volume change of 1.3 ± 0.3 mL mol^-1 and an enthalpy change of 13 ± 7 kcal mol ^-1 occurring subsequent to photodissociation of the unfolded CO bound Cc species in <∼20 ns. Two additional kinetic phases are observed that are associated with non-native His binding (ΔH and ΔV of 2 ± 4 kcal mol^-1 and -0.5 mL mol^-1, τ ∼ 2.5 μs ) and Met binding (ΔH and ΔV -0.4 ± 2 kcal mol^-1 and -0.1 ± 0.1 mL mol^-1, τ∼ 660 ns). Considering only protein conformational changes (excluding volume and enthalpies associated with heme ligation events) the initial conformational event exhibits a ΔH and ΔV of 6 ± 3 kcal mol^-1 and -3±0.1 mL mol^-1, respectively, that are attributed to a small contraction of the unfolded protein. The corresponding enthalpy associated with both native and non-native ligand binding are found to be -5±4 kcal mol^-1 (Fe-Met) and +20 ± 4 kcal mol^-1 (Fe-His) with the change in volume for both phases being essential negligible. This would indicate that non-native ligand binding likely occurs from an already collapsed conformation.

Autoxidation of Reduced Horse Heart Cytochrome c Catalyzed by Cardiolipin-Containing Membranes

Visible circular dichroism, absorption, and fluorescence spectroscopy were used to probe the binding of horse heart ferrocytochrome c to anionic cardiolipin (CL) head groups on the surface of 1,1',2,2'-tetraoleoyl cardiolipin (TOCL)/1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) (20%:80%) liposomes in an aerobic environment. We found that ferrocytochrome c undergoes a conformational transition upon binding that leads to complete oxidation of the protein at intermediate and high CL concentrations. At low lipid concentrations, the protein maintains a structure that is only slightly different from its native one, whereas an ensemble of misligated predominantly hexacoordinated low-spin states become increasingly populated at high lipid concentrations. A minor fraction of conformations with either high- or quantum-mixed-spin states were detected at a CL to protein ratio of 200 (the largest one investigated). The population of the non-native state is less pronounced than that found for cytochrome c-CL interactions initiated with oxidized cytochrome c. Under anaerobic conditions, the protein maintains its reduced state but still undergoes some conformational change upon binding to CL head groups on the liposome surface. Our data suggest that CL-containing liposomes function as catalysts by reducing the activation barrier for a Fe²⁺ → O₂ electron transfer. Adding NaCl to the existing cytochrome-liposome mixtures under aerobic conditions inhibits protein autoxidation of ferrocytochrome c and stabilizes the reduced state of the membrane-bound protein.

The role of key residues in structure, function, and stability of cytochrome-c

Cytochrome-c (cyt-c), a multi-functional protein, plays a significant role in the electron transport chain, and thus is indispensable in the energy-production process. Besides being an important component in apoptosis, it detoxifies reactive oxygen species. Two hundred and eighty-five complete amino acid sequences of cyt-c from different species are known. Sequence analysis suggests that the number of amino acid residues in most mitochondrial cyts-c is in the range 104 ± 10, and amino acid residues at only few positions are highly conserved throughout evolution. These highly conserved residues are Cys14, Cys17, His18, Gly29, Pro30, Gly41, Asn52, Trp59, Tyr67, Leu68, Pro71, Pro76, Thr78, Met80, and Phe82. These are also known as "key residues", which contribute significantly to the structure, function, folding, and stability of cyt-c. The three-dimensional structure of cyt-c from ten eukaryotic species have been determined using X-ray diffraction studies. Structure analysis suggests that the tertiary structure of cyt-c is almost preserved along the evolutionary scale. Furthermore, residues of N/C-terminal helices Gly6, Phe10, Leu94, and Tyr97 interact with each other in a specific manner, forming an evolutionary conserved interface. To understand the role of evolutionary conserved residues on structure, stability, and function, numerous studies have been performed in which these residues were substituted with different amino acids. In these studies, structure deals with the effect of mutation on secondary and tertiary structure measured by spectroscopic techniques; stability deals with the effect of mutation on T m (midpoint of heat denaturation), ∆G D (Gibbs free energy change on denaturation) and folding; and function deals with the effect of mutation on electron transport, apoptosis, cell growth, and protein expression. In this review, we have compiled all these studies at one place. This compilation will be useful to biochemists and biophysicists interested in understanding the importance of conservation of certain residues throughout the evolution in preserving the structure, function, and stability in proteins.

Constructing a man-made c-type cytochrome maquette in vivo: electron transfer, oxygen transport and conversion to a photoactive light harvesting maquette

The successful use of man-made proteins to advance synthetic biology requires both the fabrication of functional artificial proteins in a living environment, and the ability of these proteins to interact productively with other proteins and substrates in that environment. Proteins made by the maquette method integrate sophisticated oxidoreductase function into evolutionarily naive, non-computationally designed protein constructs with sequences that are entirely unrelated to any natural protein. Nevertheless, we show here that we can efficiently interface with the natural cellular machinery that covalently incorporates heme into natural cytochromes c to produce in vivo an artificial c-type cytochrome maquette. Furthermore, this c-type cytochrome maquette is designed with a displaceable histidine heme ligand that opens to allow functional oxygen binding, the primary event in more sophisticated functions ranging from oxygen storage and transport to catalytic hydroxylation. To exploit the range of functions that comes from the freedom to bind a variety of redox cofactors within a single maquette framework, this c-type cytochrome maquette is designed with a second, non-heme C, tetrapyrrole binding site, enabling the construction of an elementary electron transport chain, and when the heme C iron is replaced with zinc to create a Zn porphyrin, a light-activatable artificial redox protein. The work we describe here represents a major advance in de novo protein design, offering a robust platform for new c-type heme based oxidoreductase designs and an equally important proof-of-principle that cofactor-equipped man-made proteins can be expressed in living cells, paving the way for constructing functionally useful man-made proteins in vivo.

Investigations of the low-frequency spectral density of cytochrome c upon equilibrium unfolding

The equilibrium unfolding process of ferric horse heart cytochrome c (cyt c), induced by guanidinium hydrochloride (GdHCl), was studied using UV-vis absorption spectroscopy, resonance Raman spectroscopy, and vibrational coherence spectroscopy (VCS). The unfolding process was successfully fit using a three-state model which included the fully folded (N) and unfolded (U) states, along with an intermediate (I) assigned to a Lys bound heme. The VCS spectra revealed for the first time several low-frequency heme modes that are sensitive to cyt c unfolding: γ(a) (~50 cm(-1)), γ(b) (~80 cm(-1)), γ(c) (~100 cm(-1)), and ν(s)(His-Fe-His) at 205 cm(-1). These out-of-plane modes have potential functional relevance and are activated by protein-induced heme distortions. The free energies for the N-I and the I-U transitions at pH 7.0 and 20 °C were found to be 4.6 kcal/M and 11.6 kcal/M, respectively. Imidazole was also introduced to replace the methionine ligand so the unfolding can be modeled as a two-state system. The intensity of the mode γ(b)~80 cm(-1) remains nearly constant during the unfolding process, while the amplitudes of the other low frequency modes track with spectral changes observed at higher frequency. This confirms that the heme deformation changes are coupled to the protein tertiary structural changes that take place upon unfolding. These studies also reveal that damping of the coherent oscillations depends sensitively on the coupling between heme and the surrounding water solvent.

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.

Conformational toggling of yeast iso-1-cytochrome C in the oxidized and reduced states

To convert cyt c into a peroxidase-like metalloenzyme, the P71H mutant was designed to introduce a distal histidine. Unexpectedly, its peroxidase activity was found even lower than that of the native, and that the axial ligation of heme iron was changed to His71/His18 in the oxidized state, while to Met80/His18 in the reduced state, characterized by UV-visible, circular dichroism, and resonance Raman spectroscopy. To further probe the functional importance of Pro71 in oxidation state dependent conformational changes occurred in cyt c, the solution structures of P71H mutant in both oxidation states were determined. The structures indicate that the half molecule of cyt c (aa 50-102) presents a kind of "zigzag riveting ruler" structure, residues at certain positions of this region such as Pro71, Lys73 can move a big distance by altering the tertiary structure while maintaining the secondary structures. This finding provides a molecular insight into conformational toggling in different oxidation states of cyt c that is principle significance to its biological functions in electron transfer and apoptosis. Structural analysis also reveals that Pro71 functions as a key hydrophobic patch in the folding of the polypeptide of the region (aa 50-102), to prevent heme pocket from the solvent.

Temperature dependent equilibrium native to unfolded protein dynamics and properties observed with IR absorption and 2D IR vibrational echo experiments

Dynamic and structural properties of carbonmonoxy (CO)-coordinated cytochrome c(552) from Hydrogenobacter thermophilus (Ht-M61A) at different temperatures under thermal equilibrium conditions were studied with infrared absorption spectroscopy and ultrafast two-dimensional infrared (2D IR) vibrational echo experiments using the heme-bound CO as the vibrational probe. Depending on the temperature, the stretching mode of CO shows two distinct bands corresponding to the native and unfolded proteins. As the temperature is increased from low temperature, a new absorption band for the unfolded protein grows in and the native band decreases in amplitude. Both the temperature-dependent circular dichroism and the IR absorption area ratio R(A)(T), defined as the ratio of the area under the unfolded band to the sum of the areas of the native and unfolded bands, suggest a two-state transition from the native to the unfolded protein. However, it is found that the absorption spectrum of the unfolded protein increases its inhomogeneous line width and the center frequency shifts as the temperature is increased. The changes in line width and center frequency demonstrate that the unfolding does not follow simple two-state behavior. The temperature-dependent 2D IR vibrational echo experiments show that the fast dynamics of the native protein are virtually temperature independent. In contrast, the fast dynamics of the unfolded protein are slower than those of the native protein, and the unfolded protein fast dynamics and at least a portion of the slower dynamics of the unfolded protein change significantly, becoming faster as the temperature is raised. The temperature dependence of the absorption spectrum and the changes in dynamics measured with the 2D IR experiments confirm that the unfolded ensemble of conformers continuously changes its nature as unfolding proceeds, in contrast to the native state, which displays a temperature-independent distribution of structures.

Role of protein stabilizers on the conformation of the unfolded state of cytochrome c and its early folding kinetics: investigation at single molecular resolution

An insight into the conformation and dynamics of unfolded and early intermediate states of a protein is essential to understand the mechanism of its aggregation and to design potent inhibitor molecules. Fluorescence correlation spectroscopy has been used to study the effects of several model protein stabilizers on the conformation of the unfolded state and early folding dynamics of tetramethyl rhodamine-labeled cytochrome c from Saccharomyces cerevisiae at single molecular resolution. Special attention has been given to arginine, which is a widely used stabilizer for improving refolding yield of different proteins. The value of the hydrodynamic radius (r(H)) obtained by analyzing the intensity fluctuations of the diffusing molecules has been found to increase in a two-state manner as the protein is unfolded by urea. The results further show that the presence of arginine and other protein stabilizers favors a relatively structured conformation of the unfolded states (r(H) of 29 A) over an extended one (r(H) of 40 A), which forms in their absence. Also, the time constant of a kinetic component (tau(R)) of about 30 micros has been observed by analyzing the correlation functions, which represents formation of a collapsed state. This time constant varies with urea concentration representing an inverted Chevron plot that shows a roll-over and behavior in the absence of arginine. To the best of our knowledge, this is one of the first applications of fluorescence correlation spectroscopy to study direct folding kinetics of a protein.

Time-dependent changes in side-chain solvent accessibility during cytochrome c folding probed by pulsed oxidative labeling and mass spectrometry

The current work employs a novel approach for characterizing structural changes during the refolding of acid-denatured cytochrome c (cyt c). At various time points (ranging from 10 ms to 5 min) after a pH jump from 2 to 7, the protein is exposed to a microsecond hydroxyl radical (.OH) pulse that induces oxidative labeling of solvent-exposed side chains. Most of the covalent modifications appear as +16-Da adducts that are readily detectable by mass spectrometry. The overall extent of labeling decreases as folding proceeds, reflecting dramatic changes in the accessibility of numerous residues. Peptide mapping and tandem mass spectrometry reveal that the side chains of C14, C17, H33, F46, Y48, W59, M65, Y67, Y74, M80, I81, and Y97 are among the dominant sites of oxidation. Temporal changes in the accessibility of these residues are consistent with docking of the N- and C-terminal helices as early as 10 ms. However, structural reorganization at the helix interface takes place up to at least 1 s. Initial misligation of the heme iron by H33 leads to distal crowding, giving rise to low solvent accessibility of the displaced (native) M80 ligand and the adjacent I81. W59 retains a surprisingly high level of accessibility long into the folding process, indicating the presence of packing defects in the hydrophobically collapsed core. Overall, the results of this work are consistent with previous hydrogen/deuterium exchange studies that proposed a foldon-mediated mechanism. The structural data obtained by .OH labeling monitor the packing and burial of side chains, whereas hydrogen/deuterium exchange primarily monitors the formation of secondary structure elements. Hence, the two approaches yield complementary information. Considering the very short time scale of pulsed oxidative labeling, an extension of the approach used here to sub-millisecond folding studies should be feasible.

Cytochrome C on silica nanoparticles: influence of nanoparticle size on protein structure, stability, and activity

The structure, thermodynamic and kinetic stability, and activity of cytochrome c (cyt c) on silica nanoparticles (SNPs) of different sizes have been studied. Adsorption of cyt c onto larger SNPs results in both greater disruption of the cyt c global structure and more significant changes of the local heme microenvironment than upon adsorption onto smaller SNPs. The disruption of the heme microenvironment leads to a more solvent-accessible protein active site, as suggested by Soret circular dichroism spectroscopy and through an increase in peroxidase activity as a function of increased SNP size. Similarly, the stability of cyt c decreases more dramatically upon adsorption onto larger SNPs. These results are consistent with changes in protein-nanoparticle interactions that depend on the size or surface curvature of the supporting nanostructure. This study provides further fundamental insights into the effects of nanoscale surfaces on adsorbed protein structure and function.

Stability of the heme Fe-N-terminal amino group coordination bond in denatured cytochrome c

In the denatured states of Hydrogenobacter thermophilus cytochrome c(552) (HT) and Pseudomonas aeruginosa cytochrome c(551) (PA), and their mutants, the N-terminal amino group of the polypeptide chain is coordinated to heme Fe in place of the axial Met, the His-N(term) form being formed. The coordination of the N-terminal amino group to heme Fe leads to loop formation by the N-terminal stretch preceding the first Cys residue bound to the heme, and the N-terminal stretches of HT and PA are different from each other in terms of both the sequence and the number of constituent amino acid residues. The His-N(term) form was shown to be rather stable, and hence it can influence the stability of the denatured state. We have investigated the heme Fe coordination structures and stabilities of the His-N(term) forms emerging upon guanidine hydrochloric acid-induced unfolding of the oxidized forms of the proteins. The Fe-N(term) coordination bond in the His-N(term) form with a 9-residue N-terminal stretch of HT proteins was found to be tilted to some extent away from the heme normal, as reflected by the great heme methyl proton shift spread. On the other hand, the small heme methyl proton shift spread of the His-N(term) form with an 11-residue stretch of PA proteins indicated that its Fe-N(term) bond is nearly parallel with the heme normal. The stability of the His-N(term) form was found to be affected by the structural properties of the N-terminal stretch, such as its length and the N-terminal residue. With a given N-terminal residue, the stability of the His-N(term) form is higher for a 9-residue N-terminal stretch than an 11-residue one. In addition, with a given length of the N-terminal stretch, the His-N(term) form with an N-terminal Glu is stabilized by a few kJ mol(-1) relative to that with an N-terminal Asn. These results provide a novel insight into the stabilizing interactions in the denatured cyts c that will facilitate elucidation of the folding/unfolding mechanisms of the proteins.

Folding mechanism of reduced Cytochrome c: equilibrium and kinetic properties in the presence of carbon monoxide

Despite close structural similarity, the ferric and ferrous forms of cytochrome c differ greatly in terms of their ligand binding properties, stability, folding, and dynamics. The reduced heme iron binds diatomic ligands such as CO only under destabilizing conditions that promote weakening or disruption of native methionine-iron linkage. This makes CO a useful conformational probe for detecting partially structured states that cannot be observed in the absence of endogenous ligands. Heme absorbance, circular dichroism, and NMR were used to characterize the denaturant-induced unfolding equilibrium of ferrocytochrome c in the presence and in the absence of CO. In addition to the native state (N), which does not bind CO, and the unfolded CO complex (U-CO), a structurally distinct CO-bound form (M-CO) accumulates to high levels (approximately 75% of the population) at intermediate guanidine HCl concentrations. Comparison of the unfolding transitions for different conformational probes reveals that M-CO is a compact state containing a native-like helical core and regions of local disorder in the segment containing the native Met80 ligand and adjacent loops. Kinetic measurements of CO binding and dissociation under native, partially denaturing, and fully unfolded conditions indicate that a state M that is structurally analogous to M-CO is populated even in the absence of CO. The binding energy of the CO ligand lowers the free energy of this high-energy state to such an extent that it accumulates even under mildly denaturing equilibrium conditions. The thermodynamic and kinetic parameters obtained in this study provide a fully self-consistent description of the linked unfolding/CO binding equilibria of reduced cytochrome c.

Native and unfolded cytochrome c--comparison of dynamics using 2D-IR vibrational echo spectroscopy

Unfolded vs native CO-coordinated horse heart cytochrome c (h-cyt c) and a heme axial methionine mutant cyt c552 from Hydrogenobacter thermophilus ( Ht-M61A) are studied by IR absorption spectroscopy and ultrafast 2D-IR vibrational echo spectroscopy of the CO stretching mode. The unfolding is induced by guanidinium hydrochloride (GuHCl). The CO IR absorption spectra for both h-cyt c and Ht-M61A shift to the red as the GuHCl concentration is increased through the concentration region over which unfolding occurs. The spectra for the unfolded state are substantially broader than the spectra for the native proteins. A plot of the CO peak position vs GuHCl concentration produces a sigmoidal curve that overlays the concentration-dependent circular dichroism (CD) data of the CO-coordinated forms of both Ht-M61A and h-cyt c within experimental error. The coincidence of the CO peak shift curve with the CD curves demonstrates that the CO vibrational frequency is sensitive to the structural changes induced by the denaturant. 2D-IR vibrational echo experiments are performed on native Ht-M61A and on the protein in low- and high-concentration GuHCl solutions. The 2D-IR vibrational echo is sensitive to the global protein structural dynamics on time scales from subpicosecond to greater than 100 ps through the change in the shape of the 2D spectrum with time (spectral diffusion). At the high GuHCl concentration (5.1 M), at which Ht-M61A is essentially fully denatured as judged by CD, a very large reduction in dynamics is observed compared to the native protein within the approximately 100 ps time window of the experiment. The results suggest the denatured protein may be in a glassy-like state involving hydrophobic collapse around the heme.

Multiple routes and structural heterogeneity in protein folding

Experimental studies show that many proteins fold along sequential pathways defined by folding intermediates. An intermediate may not always be a single population of molecules but may consist of subpopulations that differ in their average structure. These subpopulations are likely to fold via independent pathways. Parallel folding and unfolding pathways appear to arise because of structural heterogeneity. For some proteins, the folding pathways can effectively switch either because different subpopulations of an intermediate get populated under different folding conditions, or because intermediates on otherwise hidden pathways get stabilized, leading to their utilization becoming discernible, or because mutations stabilize different substructures. Therefore, the same protein may fold via different pathways in different folding conditions. Multiple folding pathways make folding robust, and evolution is likely to have selected for this robustness to ensure that a protein will fold under the varying conditions prevalent in different cellular contexts.