Early events, kinetic intermediates and the mechanism of protein folding in cytochrome C

Genetically Encoded Fluorescent Probe for Detection of Heme-Induced Conformational Changes in Cytochrome c

Cytochrome c (Cytc) is a key redox protein for energy metabolism and apoptosis in cells. The activation of Cytc is composed of several steps, including its transfer to the mitochondrial membrane, binding to cytochrome c heme lyase (CCHL) and covalent attachment to heme. The spectroscopic methods are often applied to study the structural changes of Cytc. However, they require the isolation of Cytc from cells and have limited availability under physiological conditions. Despite recent studies to elucidate the tightly regulated folding mechanism of Cytc, the role of these events and their association with different conformational states remain elusive. Here, we provide a genetically encoded fluorescence method that allows monitoring of the conformational changes of Cytc upon binding to heme and CCHL. Cerulean and Venus fluorescent proteins attached at the N and C terminals of Cytc can be used to determine its unfolded, intermediate, and native states by measuring FRET amplitude. We found that the noncovalent interaction of heme in the absence of CCHL induced a shift in the FRET signal, indicating the formation of a partially folded state. The higher concentration of heme and coexpression of CCHL gave rise to the recovery of Cytc native structure. We also found that Cytc was weakly associated with CCHL in the absence of heme. As a result, a FRET-based fluorescence approach was demonstrated to elucidate the mechanism of heme-induced Cytc conformational changes with spatiotemporal resolution and can be applied to study its interaction with small molecules and other protein partners in living cells.

Dramatic Shape Changes Occur as Cytochrome c Folds

Extensive experimental studies on the folding of cytochrome c (Cyt c) make this small protein an ideal target for atomic detailed simulations for the purposes of quantitatively characterizing the structural transitions and the associated time scales for folding to the native state from an ensemble of unfolded states. We use previously generated atomically detailed folding trajectories by the stochastic difference equation in length to calculate the time-dependent changes in the small-angle X-ray scattering (SAXS) profiles. Excellent agreement is obtained between experiments and simulations for the time-dependent SAXS spectra, allowing us to identify the structures of the folding intermediates, which shows that Cyt c reaches the native state by a sequential folding mechanism. Using the ensembles of structures along the folding pathways, we show that compaction and the sphericity of Cyt c change dramatically from the prolate ellipsoid shape in the unfolded state to the spherical native state. Our data, which are in unprecedented quantitative agreement with all aspects of time-resolved SAXS experiments, show that hydrophobic collapse and amide group protection coincide on the 100 microseconds time scale, which is in accordance with ultrafast hydrogen/deuterium exchange studies. Based on these results, we propose that compaction of polypeptide chains, accompanied by dramatic shape changes, is a universal characteristic of globular proteins, regardless of the underlying folding mechanism.

The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins

From experimental studies of protein folding, it is now clear that there are two types of folding behavior, i.e., two-state folding and non-two-state folding, and understanding the relationships between these apparently different folding behaviors is essential for fully elucidating the molecular mechanisms of protein folding. This article describes how the presence of the two types of folding behavior has been confirmed experimentally, and discusses the relationships between the two-state and the non-two-state folding reactions, on the basis of available data on the correlations of the folding rate constant with various structure-based properties, which are determined primarily by the backbone topology of proteins. Finally, a two-stage hierarchical model is proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the molten globule state as a folding intermediate. The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

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.

On the preservation of non-covalent protein complexes during electrospray ionization

The application range of electrospray ionization mass spectrometry for the quantitative determination of stoichiometries and binding constants for non-covalent protein complexes is broadly discussed. The underlying fundamental question is whether or not the original molecular equilibrium can be preserved during the ionization process and be revealed by subsequent mass spectrometry analysis. Here, we take a new look at this question by discussing recent studies in droplet chemistry.This article is part of the themed issue 'Quantitative mass spectrometry'.

Microsecond protein folding events revealed by time-resolved fluorescence resonance energy transfer in a microfluidic mixer

We demonstrate the combination of the time-resolved fluorescence resonance energy transfer (tr-FRET) measurement and the ultrarapid hydrodynamic focusing microfluidic mixer. The combined technique is capable of probing the intermolecular distance change with temporal resolution at microsecond level and structural resolution at Angstrom level, and the use of two-photon excitation enables a broader exploration of FRET with spectrum from near-ultraviolet to visible wavelength. As a proof of principle, we used the coupled microfluidic laminar flow and time-resolved two-photon excitation microscopy to investigate the early folding states of Cytochrome c (cyt c) by monitoring the distance between the tryptophan (Trp-59)-heme donor-acceptor (D-A) pair. The transformation of folding states of cyt c in the early 500 μs of refolding was revealed on the microsecond time scale. For the first time, we clearly resolved the early transient state of cyt c, which is populated within the dead time of the mixer (<10 μs) and has a characteristic Trp-59-heme distance of ∼31 Å. We believe this tool can find more applications in studying the early stages of biological processes with FRET as the probe.

Heterogeneity of equilibrium molten globule state of cytochrome c induced by weak salt denaturants under physiological condition

While many proteins are recognized to undergo folding via intermediate(s), the heterogeneity of equilibrium folding intermediate(s) along the folding pathway is less understood. In our present study, FTIR spectroscopy, far- and near-UV circular dichroism (CD), ANS and tryptophan fluorescence, near IR absorbance spectroscopy and dynamic light scattering (DLS) were used to study the structural and thermodynamic characteristics of the native (N), denatured (D) and intermediate state (X) of goat cytochorme c (cyt-c) induced by weak salt denaturants (LiBr, LiCl and LiClO₄) at pH 6.0 and 25°C. The LiBr-induced denaturation of cyt-c measured by Soret absorption (Δε₄₀₀) and CD ([θ]₄₀₉), is a three-step process, N ↔ X ↔ D. It is observed that the X state obtained along the denaturation pathway of cyt-c possesses common structural and thermodynamic characteristics of the molten globule (MG) state. The MG state of cyt-c induced by LiBr is compared for its structural and thermodynamic parameters with those found in other solvent conditions such as LiCl, LiClO₄ and acidic pH. Our observations suggest: (1) that the LiBr-induced MG state of cyt-c retains the native Met80-Fe(III) axial bond and Trp59-propionate interactions; (2) that LiBr-induced MG state of cyt-c is more compact retaining the hydrophobic interactions in comparison to the MG states induced by LiCl, LiClO₄ and 0.5 M NaCl at pH 2.0; and (3) that there exists heterogeneity of equilibrium intermediates along the unfolding pathway of cyt-c as highly ordered (X1), classical (X2) and disordered (X3), i.e., D ↔ X3 ↔ X2 ↔ X1 ↔ N.

Microsecond barrier-limited chain collapse observed by time-resolved FRET and SAXS

It is generally held that random-coil polypeptide chains undergo a barrier-less continuous collapse when the solvent conditions are changed to favor the fully folded native conformation. We test this hypothesis by probing intramolecular distance distributions during folding in one of the paradigms of folding reactions, that of cytochrome c. The Trp59-to-heme distance was probed by time-resolved Förster resonance energy transfer in the microsecond time range of refolding. Contrary to expectation, a state with a Trp59-heme distance close to that of the guanidinium hydrochloride (GdnHCl) denatured state is present after ~27 μs of folding. A concomitant decrease in the population of this state and an increase in the population of a compact high-FRET (Förster resonance energy transfer) state (efficiency>90%) show that the collapse is barrier limited. Small-angle X-ray scattering (SAXS) measurements over a similar time range show that the radius of gyration under native favoring conditions is comparable to that of the GdnHCl denatured unfolded state. An independent comprehensive global thermodynamic analysis reveals that marginally stable partially folded structures are also present in the nominally unfolded GdnHCl denatured state. These observations suggest that specifically collapsed intermediate structures with low stability in rapid equilibrium with the unfolded state may contribute to the apparent chain contraction observed in previous fluorescence studies using steady-state detection. In the absence of significant dynamic averaging of marginally stable partially folded states and with the use of probes sensitive to distance distributions, barrier-limited chain contraction is observed upon transfer of the GdnHCl denatured state ensemble to native-like conditions.

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.

Ultrafast hydrogen exchange reveals specific structural events during the initial stages of folding of cytochrome c

Many proteins undergo a sharp decrease in chain dimensions during early stages of folding, prior to the rate-limiting step in folding. However, it remains unclear whether compact states are the result of specific folding events or a general hydrophobic collapse of the poly peptide chain driven by the change in solvent conditions. To address this fundamental question, we extended the temporal resolution of NMR-detected H/D exchange labeling experiments into the microsecond regime by adopting a microfluidics approach. By observing the competition between H/D exchange and folding as a function of labeling pH, coupled with direct measurement of exchange rates in the unfolded state, we were able to monitor hydrogen-bond formation for over 50 individual backbone NH groups within the initial 140 microseconds of folding of horse cytochrome c. Clusters of solvent-shielded amide protons were observed in two α-helical segments in the C-terminal half of the protein, while the N-terminal helix remained largely unstructured, suggesting that proximity in the primary structure is a major factor in promoting helix formation and association at early stages of folding, while the entropically more costly long-range contacts between the N- and C-terminal helices are established only during later stages. Our findings clearly indicate that the initial chain condensation in cytochrome c is driven by specific interactions among a subset of α-helical segments rather than a general hydrophobic collapse.

Proteins with highly similar native folds can show vastly dissimilar folding behavior when desolvated

Proteins can be exposed to vastly different environments such as the cytosol or membranes, but the delicate balance between external factors and intrinsic determinants of protein structure, stability, and folding is only poorly understood. Here we used electron capture dissociation to study horse and tuna heart Cytochromes c in the complete absence of solvent. The significantly different stability of their highly similar native folds after transfer into the gas phase, and their strikingly different folding behavior in the gas phase, can be rationalized on the basis of electrostatic interactions such as salt bridges. In the absence of hydrophobic bonding, protein folding is far slower and more complex than in solution.

Effects of macromolecular crowding on burst phase kinetics of cytochrome c folding

Excluded volume and viscosity effects of crowding agents that mimic crowded conditions in vivo on "classical" burst phase folding kinetics of cytochrome c are assessed in vitro. Upon electron transfer-triggered folding of reduced cytochrome c, far-UV time-resolved circular dichroism (TRCD) is used to monitor folding under different conditions. Earlier work has shown that folding of reduced cytochrome c from the guanidinium hydrochloride-induced unfolded ensemble in dilute phosphate buffer involves kinetic partitioning: one fraction of molecules folds rapidly, on a time scale identical to that of reduction, while the remaining population folds more slowly. In the presence of 220 mg/mL dextran 70, a synthetic macromolecular crowding agent that occupies space but does not interact with proteins, the population of the fast folding step for cytochrome c is greatly reduced. Increasing the viscosity with sucrose to the same microviscosity exhibited by the dextran solution showed no significant decrease in the amplitude of the fast-folding phase of cytochrome c. Experiments show that the unfolded-state heme ligation remains bis-His in the presence of dextran 70, but coarse-grained simulations suggest that the unfolded-state ensemble becomes more compact in the presence of crowders. We conclude that excluded volume effects alter unfolded cytochrome c such that access to fast-folding conformations is reduced.

Interconnection of salt-induced hydrophobic compaction and secondary structure formation depends on solution conditions: revisiting early events of protein folding at single molecule resolution

What happens in the early stage of protein folding remains an interesting unsolved problem. Rapid kinetics measurements with cytochrome c using submillisecond continuous flow mixing devices suggest simultaneous formation of a compact collapsed state and secondary structure. These data seem to indicate that collapse formation is guided by specific short and long range interactions (heteropolymer collapse). A contrasting interpretation also has been proposed, which suggests that the collapse formation is rapid, nonspecific, and a trivial solvent related compaction, which could as well be observed by a homopolymer (homopolymer collapse). We address this controversy using fluorescence correlation spectroscopy (FCS), which enables us to monitor the salt-induced compaction accompanying collapse formation and the associated time constant directly at single molecule resolution. In addition, we follow the formation of secondary structure using far UV CD. The data presented here suggest that both these models (homopolymer and heteropolymer) could be applicable depending on the solution conditions. For example, the formation of secondary structure and compact state is not simultaneous in aqueous buffer. In aqueous buffer, formation of the compact state occurs through a two-state co-operative transition following heteropolymer formalism, whereas secondary structure formation takes place gradually. In contrast, in the presence of urea, a compaction of the protein radius occurs gradually over an extended range of salt concentration following homopolymer formalism. The salt-induced compaction and the formation of secondary structure take place simultaneously in the presence of urea.

Nanosecond time-resolved polarization spectroscopies: tools for probing protein reaction mechanisms

Polarization methods, introduced in the 1800s, offered one of the earliest ways to examine protein structure. Since then, many other structure-sensitive probes have been developed, but circular dichroism (CD) remains a powerful technique because of its versatility and the specificity of protein structural information that can be explored. With improvements in time resolution, from millisecond to picosecond CD measurements, it has proven to be an important tool for studying the mechanism of folding and function in many biomolecules. For example, nanosecond time-resolved CD (TRCD) studies of the sub-microsecond events of reduced cytochrome c folding have provided direct experimental evidence of kinetic heterogeneity, which is an inherent property of the diffusional nature of early folding dynamics on the energy landscape. In addition, TRCD has been applied to the study of many biochemical processes, such as ligand rebinding in hemoglobin and myoglobin and signaling state formation in photoactive yellow protein and prototropin 1 LOV2. The basic approach to TRCD has also been extended to include a repertoire of nanosecond polarization spectroscopies: optical rotatory dispersion (ORD), magnetic CD and ORD, and linear dichroism. This article will discuss the details of the polarization methods used in this laboratory, as well as the coupling of time-resolved ORD with the temperature-jump trigger so that protein folding can be studied in a larger number of proteins.

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.

Magnetic circular dichroism spectroscopy as a probe of the structures of the metal sites in metalloproteins

Magnetic circular dichroism (MCD) is a powerful probe of the electronic and geometric structures of metal centres in metalloproteins. MCD has provided significant insight into the nature of the axial donors at haem centres and, more recently, sophisticated methods for the analysis of MCD spectra have had a major impact on the study of the electronic structures of the ground states of a range of Cu, non-haem iron and Mo-containing active sites. This detail, together with data from other complimentary spectroscopies, has played a major role in defining the chemistry underpinning the catalysis achieved by these metal centres.

The human alpha-lactalbumin molten globule: comparison of structural preferences at pH 2 and pH 7

Structural investigations of molten globules provide an important contribution towards understanding protein folding pathways. A close similarity between equilibrium molten globule states and kinetic species observed during refolding has been reported for several proteins. However, the experimental conditions, and in particular the pH, under which the equilibrium and kinetic species are studied often differ significantly. For human α-lactalbumin (α-LA), the equilibrium molten globule is most often studied at pH 2, the so-called A-state, while kinetic refolding experiments are performed at neutral pH. α-LA contains a large number of acidic amino acid residues that may influence the properties of the molten globule differently at low and neutral pH. In this study, we investigate the structural preferences of the α-LA molten globule at pH 7 at the level of individual residues using nuclear magnetic resonance spectroscopy and compare these data with previous results obtained at pH 2. We show that differences exist in the conformational ensemble that describes the α-LA molten globule at these two pH values. The molten globule at pH 7 is generally less stable than that at the low pH A-state. Most notable are differences in the stability of structure for the C-helix and the calcium-binding loop that precedes it and differences in the contribution of long-range hydrophobic contacts between the N-terminal and C-terminal regions of the α-domain to the stability of the molten globule. Our results are discussed in the context of previous studies of the α-LA molten globule and can be used to reconcile apparent discrepancies in published data relating to the C-helix. In the light of our results, the low pH A-state may not be the best model for the kinetic molten globule observed during refolding of α-LA. The pH-dependent effects reported here for α-LA may be of relevance in comparisons of equilibrium and kinetic molten globules of other proteins.