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.

The reversible condensation and expansion of the rotavirus genome

Understanding the structural organization of the genome is particularly relevant in segmented double-stranded RNA viruses, which exhibit endogenous transcription activity. These viruses are molecular machines capable of repeated cycles of transcription within the intact capsid. Rotavirus, a major cause of infantile gastroenteritis, is a prototypical segmented double-stranded RNA virus. From our three-dimensional structural analyses of rotavirus examined under various chemical conditions using electron cryomicroscopy, we show here that the viral genome exhibits a remarkable conformational flexibility by reversibly changing its packaging density. In the presence of ammonium ions at high pH, the genome condenses to a radius of ≈180 Å from ≈220 Å. Upon returning to physiological conditions, the genome re-expands and fully maintains its transcriptional properties. These studies provide further insights into the genome organization and suggest that the observed isometric and concentric nature of the condensation is due to strong interactions between the genome core and the transcription enzymes anchored to the capsid inner surface. The ability of the genome to condense beyond what is normally observed in the native virus indicates that the negative charges on the RNA in the native state may be only partially neutralized. Partial neutralization may be required to maintain appropriate interstrand spacing for templates to move around the enzyme complexes during transcription. Genome condensation was not observed either with increased cation concentrations at normal pH or at high pH without ammonium ions. This finding indicates that the observed genome condensation is a synergistic effect of hydroxyl and ammonium ions involving disruption of protein–RNA interactions that perhaps facilitate further charge neutralization and consequent reduction in the interstrand spacing.

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.

Destabilization of Ca2+-free gelsolin may not be responsible for proteolysis in Familial Amyloidosis of Finnish Type

Mutations at position 187 in secreted gelsolin enable aberrant proteolysis at the 172–173 and 243–244 amide bonds, affording the 71-residue amyloidogenic peptide deposited in Familial Amyloidosis of Finnish Type (FAF). Thermodynamic comparisons of two different domain 2 constructs were carried out to study possible effects of the mutations on proteolytic susceptibility. In the construct we consider to be most representative of domain 2 in the context of the full-length protein (134–266), the D187N FAF variant is slightly destabilized relative to wild type (WT) under the conditions of urea denaturation, but exhibits a Tm identical to WT. The D187Y variant is less stable to intermediate urea concentrations and exhibits a Tm that is estimated to be ≈5°C lower than WT (pH 7.4, Ca2+-free). Although the thermodynamic data indicate that the FAF mutations may slightly destabilize domain 2, these changes are probably not sufficient to shift the native to denatured state equilibrium enough to enable the proteolysis leading to FAF. Biophysical data indicate that these two FAF variants may have different native state structures and possibly different pathways of amyloidosis.

Kinetic stability as a mechanism for protease longevity

The folding of the extracellular serine protease, α-lytic protease (αLP; EC 3.4.21.12) reveals a novel mechanism for stability that appears to lead to a longer functional lifetime for the protease. For αLP, stability is based not on thermodynamics, but on kinetics. Whereas this has required the coevolution of a pro region to facilitate folding, the result has been the optimization of native-state properties independent of their consequences on thermodynamic stability. Structural and mutational data lead to a model for catalysis of folding in which the pro region binds to a conserved β-hairpin in the αLP C-terminal domain, stabilizing the folding transition state and the native state. The pro region is then proteolytically degraded, leaving the active αLP trapped in a metastable conformation. This metastability appears to be a consequence of pressure to evolve properties of the native state, including a large, highly cooperative barrier to unfolding, and extreme rigidity, that reduce susceptibility to proteolytic degradation. In a test of survival under highly proteolytic conditions, homologous mammalian proteases that have not evolved kinetic stability are much more rapidly degraded than αLP. Kinetic stability as a means to longevity is likely to be a mechanism conserved among the majority of extracellular bacterial pro-proteases and may emerge as a general strategy for intracellular eukaryotic proteases subject to harsh conditions as well.

Chaperone rings in protein folding and degradation

Chaperone rings play a vital role in the opposing ATP-mediated processes of folding and degradation of many cellular proteins, but the mechanisms by which they assist these life and death actions are only beginning to be understood. Ring structures present an advantage to both processes, providing for compartmentalization of the substrate protein inside a central cavity in which multivalent, potentially cooperative interactions can take place between the substrate and a high local concentration of binding sites, while access of other proteins to the cavity is restricted sterically. Such restriction prevents outside interference that could lead to nonproductive fates of the substrate protein while it is present in non-native form, such as aggregation. At the step of recognition, chaperone rings recognize different motifs in their substrates, exposed hydrophobicity in the case of protein-folding chaperonins, and specific “tag” sequences in at least some cases of the proteolytic chaperones. For both folding and proteolytic complexes, ATP directs conformational changes in the chaperone rings that govern release of the bound polypeptide. In the case of chaperonins, ATP enables a released protein to pursue the native state in a sequestered hydrophilic folding chamber, and, in the case of the proteases, the released polypeptide is translocated into a degradation chamber. These divergent fates are at least partly governed by very different cooperating components that associate with the chaperone rings: that is, cochaperonin rings on one hand and proteolytic ring assemblies on the other. Here we review the structures and mechanisms of the two types of chaperone ring system.

Folding funnels and frustration in off-lattice minimalist protein landscapes

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a β-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.

beta-Lactamase binds to GroEL in a conformation highly protected against hydrogen/deuterium exchange.

Escherichia coli RTEM beta-lactamase reversibly forms a stable complex with GroEL, devoid of any enzymatic activity, at 48 degrees C. When beta-lactamase is diluted from this complex into denaturant solution, its unfolding rate is identical to that from the native state, while the unfolding rate from the molten globule state is too fast to be measured. Electrospray mass spectrometry shows that the rate of proton exchange in beta-lactamase in the complex at 48 degrees C is slower than in the absence of GroEL at the same temperature, and resembles the exchange of the native state at 25 degrees C. Similarly, the final number of protected deuterons is higher in the presence of GroEL than in its absence. We conclude that, for beta-lactamase, a state with significant native structure is bound to GroEL. Thus, different proteins are recognized by GroEL in very different states, ranging from totally unfolded to native-like, and this recognition may depend on which state can provide sufficient accessible hydrophobic amino acids in a suitably clustered arrangement. Reversible binding of native-like states with hydrophobic patches may be an important property of GroEL to protect the cell from aggregating protein after heat-shock.

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.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

X-ray structure of a vanadium-containing enzyme: chloroperoxidase from the fungus Curvularia inaequalis.

The chloroperoxidase (EC 1.11.1.-) from the fungus Curvularia inaequalis belongs to a class of vanadium enzymes that oxidize halides in the presence of hydrogen peroxide to the corresponding hypohalous acids. The 2.1 A crystal structure (R = 20%) of an azide chloroperoxidase complex reveals the geometry of the catalytic vanadium center. Azide coordinates directly to the metal center, resulting in a structure with azide, three nonprotein oxygens, and a histidine as ligands. In the native state vanadium will be bound as hydrogen vanadate(V) in a trigonal bipyramidal coordination with the metal coordinated to three oxygens in the equatorial plane, to the OH group at one apical position, and to the epsilon 2 nitrogen of a histidine at the other apical position. The protein fold is mainly alpha-helical with two four-helix bundles as main structural motifs and an overall structure different from other structures. The helices pack together to a compact molecule, which explains the high stability of the protein. An amino acid sequence comparison with vanadium-containing bromoperoxidase from the seaweed Ascophyllum nodosum shows high similarities in the regions of the metal binding site, with all hydrogen vanadate(V) interacting residues conserved except for lysine-353, which is an asparagine.

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.

Kinetic traps in lysozyme folding.

Folding of lysozyme from hen egg white was investigated by using interrupted refolding experiments. This method makes use of a high energy barrier between the native state and transient folding intermediates, and, in contrast to conventional optical techniques, it enables one to specifically monitor the amount of native molecules during protein folding. The results show that under strongly native conditions lysozyme can refold on parallel pathways. The major part of the lysozyme molecules (86%) refold on a slow kinetic pathway with well-populated partially folded states. Additionally, 14% of the molecules fold faster. The rate constant of formation of native molecules on the fast pathway corresponds well to the rate constant expected for folding to occur by a two-state process without any detectable intermediates. The results suggest that formation of the native state for the major fraction of lysozyme molecules is retarded compared with the direct folding process. Partially structured intermediates that transiently populate seem to be kinetically trapped in a conformation that can only slowly reach the native structure.

Lattice model for rapidly folding protein-like heteropolymers.

Protein folding is a relatively fast process considering the astronomical number of conformations in which a protein could find itself. Within the framework of a lattice model, we show that one can design rapidly folding sequences by assigning the strongest attractive couplings to the contacts present in a target native state. Our protein design can be extended to situations with both attractive and repulsive contacts. Frustration is minimized by ensuring that all the native contacts are again strongly attractive. Strikingly, this ensures the inevitability of folding and accelerates the folding process by an order of magnitude. The evolutionary implications of our findings are discussed.

On constructing folding heteropolymers.

Simplified models of the protein-folding process have led to valuable insights into the generic properties of the folding of heteropolymers. On the basis of theoretical arguments, Shakhnovich and Gutin [(1993) Proc. Natl. Acad. Sci. USA 90, 7195-7199] have proposed a specific method to generate folding sequences for one of these. Here we present a model of folding in heteropolymers that is comparable in simplicity but different in spirit to the one studied by Shakhnovich and Gutin. In our model, the proposed recipe for constructing folding sequence fails. We find that, as a rule, the construction of folding sequences is impossible to achieve by looking at the native conformation only. Rather, competing conformations have to be taken into account too. An evolutionary algorithm that generates folding sequences by optimizing both stability of the native state and folding time is described. Remarkably, this algorithm produces, among others, sequences that fold reproducibly to metastable states.