RNA folding kinetics regulates translation of phage MS2 maturation gene

The gene for the maturation protein of the single-stranded RNA coliphage MS2 is preceded by an untranslated leader of 130 nt, which folds into a cloverleaf, i.e., three stem–loop structures enclosed by a long distance interaction (LDI). This LDI prevents translation because its 3′ moiety contains the Shine–Dalgarno sequence of the maturation gene. Previously, several observations suggested that folding of the cloverleaf is kinetically delayed, providing a time window for ribosomes to access the RNA. Here we present direct evidence for this model. In vitro experiments show that ribosome binding to the maturation gene is faster than refolding of the denatured cloverleaf. This folding delay appears related to special properties of the leader sequence. We have replaced the three stem–loop structures by a single five nt loop. This change does not affect the equilibrium structure of the LDI. Nevertheless, in this construct, the folding delay has virtually disappeared, suggesting that now the RNA folds faster than ribosomes can bind. Perturbation of the cloverleaf by an insertion makes the maturation start permanently accessible. A pseudorevertant that evolved from an infectious clone carrying the insertion had overcome this defect. It showed a wild-type folding delay before closing down the maturation gene. This experiment reveals the biological significance of retarded cloverleaf formation.

Protein folding kinetics exhibit an Arrhenius temperature dependence when corrected for the temperature dependence of protein stability

The anomalous temperature dependence of protein folding has received considerable attention. Here we show that the temperature dependence of the folding of protein L becomes extremely simple when the effects of temperature on protein stability are corrected for; the logarithm of the folding rate is a linear function of 1/T on constant stability contours in the temperature–denaturant plane. This convincingly demonstrates that the anomalous temperature dependence of folding derives from the temperature dependence of the interactions that stabilize proteins, rather than from the super Arrhenius temperature dependence predicted for the configurational diffusion constant on a rough energy landscape. However, because of the limited temperature range accessible to experiment, the results do not rule out models with higher order temperature dependences. The significance of the slope of the stability-corrected Arrhenius plots is discussed.

Solvent effects on the energy landscapes and folding kinetics of polyalanine

The effect of a solvation on the thermodynamics and kinetics of polyalanine (Ala12) is explored on the basis of its energy landscapes in vacuum and in an aqueous solution. Both energy landscapes are characterized by two basins, one associated with α-helical structures and the other with coil and β-structures of the peptide. In both environments, the basin that corresponds to the α-helical structure is considerably narrower than the basin corresponding to the β-state, reflecting their different contributions to the entropy of the peptide. In vacuum, the α-helical state of Ala12 constitutes the native state, in agreement with common helical propensity scales, whereas in the aqueous medium, the α-helical state is destabilized, and the β-state becomes the native state. Thus solvation has a dramatic effect on the energy landscape of this peptide, resulting in an inverted stability of the two states. Different folding and unfolding time scales for Ala12 in hydrophilic and hydrophobic chemical environments are caused by the higher entropy of the native state in water relative to vacuum. The concept of a helical propensity has to be extended to incorporate environmental solvent effects.

Simple model of protein folding kinetics.

A simple model of the kinetics of protein folding is presented. The reaction coordinate is the "correctness" of a configuration compared with the native state. The model has a gap in the energy spectrum, a large configurational entropy, a free energy barrier between folded and partially folded states, and a good thermodynamic folding transition. Folding kinetics is described by a master equation. The folding time is estimated by means of a local thermodynamic equilibrium assumption and then is calculated both numerically and analytically by solving the master equation. The folding time has a maximum near the folding transition temperature and can have a minimum at a lower temperature.

Probing the kinetics of single molecule protein folding

We propose an approach to integrate the theory, simulations, and experiments in protein-folding kinetics. This is realized by measuring the mean and high-order moments of the first-passage time and its associated distribution. The full kinetics is revealed in the current theoretical framework through these measurements. In the experiments, information about the statistical properties of first-passage times can be obtained from the kinetic folding trajectories of single molecule experiments ( for example, fluorescence). Theoretical/simulation and experimental approaches can be directly related. We study in particular the temperature-varying kinetics to probe the underlying structure of the folding energy landscape. At high temperatures, exponential kinetics is observed; there are multiple parallel kinetic paths leading to the native state. At intermediate temperatures, nonexponential kinetics appears, revealing the nature of the distribution of local traps on the landscape and, as a result, discrete kinetic paths emerge. At very low temperatures, exponential kinetics is again observed; the dynamics on the underlying landscape is dominated by a single barrier. The ratio between first-passage-time moments is proposed to be a good variable to quantitatively probe these kinetic changes. The temperature-dependent kinetics is consistent with the strange kinetics found in folding dynamics experiments. The potential applications of the current results to single-molecule protein folding are discussed.

Observation of strange kinetics in protein folding

Highly nonexponential folding kinetics in aqueous solution have been observed during temperature jump-induced refolding of two proteins, yeast phosphoglycerate kinase and a ubiquitin mutant. The observations are most easily interpreted in terms of downhill folding, which posits a heterogeneous ensemble of structures en route to the folded state. The data are also reconciled with exponential kinetics measured under different experimental conditions and with titration experiments indicating cooperative folding.

Magnesium-dependent folding of self-splicing RNA: Exploring the link between cooperativity, thermodynamics, and kinetics

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

Coordinate and time-dependent diffusion dynamics in protein folding

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Configuration-Dependent Diffusion Dynamics of Downhill and Two-State Protein Folding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion. Including the configurational dependence will challenge the transition state theory of protein folding. The classical transition state theory will have to be modified to be consistent. The more detailed folding mechanistic studies involving phi value analysis based on the classical transition state theory also will have to be modified quantitatively.

STABILITY AND FOLDING OF MUTANT FORMS OF ISO-1 CYTOCHROME-C FROM SACCHAROMYCES CEREVISIAE

Partially functional forms of iso-1-cytochrome c from Saccharomyces cerevisiae were obtained by replacements of the evolutionarily conserved proline 71 with valine, isoleucine and threonine (Ernst et.al.,1985). Pro-71 lies at the juncture of two short helical regions and is believed to be important for proper local polypeptide chain folding within the iso-1-cytochrome c structure.^ To study folding in the absence of intermolecular disulfide dimer formation the free sulfhydryl group of Cys-102 was modified in both wild type and mutant proteins with an alkylating reagent, methyl methanethiosulfonate. Spectral analysis of the wild type and mutant proteins shows that the native-like functional (or partially functional) folded structure of cytochrome c is retained in the chemically modified derivatives. The replacement of Pro-71 with valine, isoleucine or threonine reduces the intensity of the 696 nm absorbance band which is an indicator of the Met-80 ligation to the heme. Thermal stability and guanidine hydrochloride unfolding studies of the mutant proteins shows a destabilization of the protein as a result of mutation. The degree of destabilization depends on the chemical nature of the substituent amino acid in the mutant protiens.^ Kinetics of folding/unfolding reactions of the proteins were monitored by fluorescence changes using stopped flow mixing to obtain guanidine hydrochloride concentration jumps ending below, within, and above the transition zone. The replacement of Pro-71 alters the rate on one of the fastest phases, $\tau\sb3$, while the two other phases, $\tau\sb1$ & $\tau\sb2$, remain the same.^ Slow refolding kinetic studies indicate that replacement of Pro-71 does not completely eliminate the absorbance or fluorescence detected slow phases leading to the conclusion that Pro-71 is not involved in the generation of the slow phases in the folding kinetics of iso-1-cytochrome c.^ The alkaline conformational change involving the disappearance of the 696 nm absorbance band occurs with increasing pH in the alkaline pH region (Davis et al., 1974). The apparent pK of this conformational change in mutant proteins is shifted as much as two pH units compared to wild type. The equilibrium and kinetic data of alkaline transition for the wild type follows a simple mechanism proposed by Davis et al., (1974) for horse heart cytochrome c. A more complex mechanism is proposed for the behavior of the mutant proteins. ^

An optimal Mg2+ concentration for kinetic folding of the Tetrahymena ribozyme

Divalent metal ions, such as Mg2+, are generally required for tertiary structure formation in RNA. Although the role of Mg2+ binding in RNA-folding equilibria has been studied extensively, little is known about the role of Mg2+ in RNA-folding kinetics. In this paper, we explore the effect of Mg2+ on the rate-limiting step in the kinetic folding pathway of the Tetrahymena ribozyme. Analysis of these data reveals the presence of a Mg2+-stabilized kinetic trap that slows folding at higher Mg2+ concentrations. Thus, the Tetrahymena ribozyme folds with an optimal rate at 2 mM Mg2+, just above the concentration required for stable structure formation. These results suggest that thermodynamic and kinetic folding of RNA are cooptimized at a Mg2+ concentration that is sufficient to stabilize the folded form but low enough to avoid kinetic traps and misfolding.

Cytochrome c′ folding triggered by electron transfer: Fast and slow formation of four-helix bundles

Reduced (FeII) Rhodopseudomonas palustris cytochrome c′ (Cyt c′) is more stable toward unfolding ([GuHCl]1/2 = 2.9(1) M) than the oxidized (FeIII) protein ([GuHCl]1/2 = 1.9(1) M). The difference in folding free energies (ΔΔGf° = 70 meV) is less than half of the difference in reduction potentials of the folded protein (100 mV vs. NHE) and a free heme in aqueous solution (≈−150 mV). The spectroscopic features of unfolded FeII–Cyt c′ indicate a low-spin heme that is axially coordinated to methionine sulfur (Met-15 or Met-25). Time-resolved absorption measurements after CO photodissociation from unfolded FeII(CO)–Cyt c′ confirm that methionine can bind to the ferroheme on the microsecond time scale [kobs = 5(2) × 104 s−1]. Protein folding was initiated by photoreduction (two-photon laser excitation of NADH) of unfolded FeIII–Cyt c′ ([GuHCl] = 2.02–2.54 M). Folding kinetics monitored by heme absorption span a wide time range and are highly heterogeneous; there are fast-folding (≈103 s−1), intermediate-folding (102–101 s−1), and slow-folding (10−1 s−1) populations, with the last two likely containing methionine-ligated (Met-15 or Met-25) ferrohemes. Kinetics after photoreduction of unfolded FeIII–Cyt c′ in the presence of CO are attributable to CO binding [1.4(6) × 103 s−1] and FeII(CO)–Cyt c′ folding [2.8(9) s−1] processes; stopped-flow triggered folding of FeIII–Cyt c′ (which does not contain a protein-derived sixth ligand) is adequately described by a single kinetics phase with an estimated folding time constant of ≈4 ms [ΔGf° = −33(3) kJ mol−1] at zero denaturant.

Chaperonin-facilitated protein folding: optimization of rate and yield by an iterative annealing mechanism.

We develop a heuristic model for chaperonin-facilitated protein folding, the iterative annealing mechanism, based on theoretical descriptions of "rugged" conformational free energy landscapes for protein folding, and on experimental evidence that (i) folding proceeds by a nucleation mechanism whereby correct and incorrect nucleation lead to fast and slow folding kinetics, respectively, and (ii) chaperonins optimize the rate and yield of protein folding by an active ATP-dependent process. The chaperonins GroEL and GroES catalyze the folding of ribulose bisphosphate carboxylase at a rate proportional to the GroEL concentration. Kinetically trapped folding-incompetent conformers of ribulose bisphosphate carboxylase are converted to the native state in a reaction involving multiple rounds of quantized ATP hydrolysis by GroEL. We propose that chaperonins optimize protein folding by an iterative annealing mechanism; they repeatedly bind kinetically trapped conformers, randomly disrupt their structure, and release them in less folded states, allowing substrate proteins multiple opportunities to find pathways leading to the most thermodynamically stable state. By this mechanism, chaperonins greatly expand the range of environmental conditions in which folding to the native state is possible. We suggest that the development of this device for optimizing protein folding was an early and significant evolutionary event.

Submillisecond folding of monomeric lambda repressor.

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.