Time-resolved RNA SHAPE chemistry

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry yields quantitative RNA secondary and tertiary structure information at single nucleotide resolution. SHAPE takes advantage of the discovery that the nucleophilic reactivity of the ribose 2'-hydroxyl group is modulated by local nucleotide flexibility in the RNA backbone. Flexible nucleotides are reactive toward hydroxyl-selective electrophiles, whereas constrained nucleotides are unreactive. Initial versions of SHAPE chemistry, which employ isatoic anhydride derivatives that react on the minute time scale, are emerging as the ideal technology for monitoring equilibrium structures of RNA in a wide variety of biological environments. Here, we extend SHAPE chemistry to a benzoyl cyanide scaffold to make possible facile time-resolved kinetic studies of RNA in approximately 1 s snapshots. We then use SHAPE chemistry to follow the time-dependent folding of an RNase P specificity domain RNA. Tertiary interactions form in two distinct steps with local tertiary contacts forming an order of magnitude faster than long-range interactions. Rate-determining tertiary folding requires minutes despite that no non-native interactions must be disrupted to form the native structure. Instead, overall folding is limited by simultaneous formation of interactions approximately 55 A distant in the RNA. Time-resolved SHAPE holds broad potential for understanding structural biogenesis and the conformational interconversions essential to the functions of complex RNA molecules at single nucleotide resolution.     

Topology-based models and NMR structures in protein folding simulations

Topology-based interaction potentials are simplified models that use the native contacts in the folded structure of a protein to define an energetically unfrustrated folding funnel. They have been widely used to analyze the folding transition and pathways of different proteins through computer simulations. Obviously, they need a reliable, experimentally determined folded structure to define the model interactions. In structures elucidated through NMR spectroscopy, a complex treatment of the raw experimental data usually provides a series of models, a set of different conformations compatible with the available experimental data. Here, we use an efficient coarse-grained simulation technique to independently consider the contact maps from every different NMR model in a protein whose structure has been resolved by the use of NMR spectroscopy. For lambda-Cro repressor, a homodimeric protein, we have analyzed its folding characteristics with a topology-based model. We have focused on the competition between the folding of the individual chains and their binding to form the final quaternary structure. From 20 different NMR models, we find a predominant three-state folding behavior, in agreement with experimental data on the folding pathway for this protein. Individual NMR models, however, show distinct characteristics, which are analyzed both at the level of the interplay between tertiary/quaternary structure formation and also regarding the thermal stability of the tertiary structure of every individual chain.     

Importance of native-state topology for determining the folding rate of two-state proteins

Understanding the relationship between amino acid sequences and folding rate of proteins is a challenging task similar to protein folding problem. In this work, we have analyzed the relative importance of protein sequence and structure for predicting the protein folding rates in terms of amino acid properties and contact distances, respectively. We found that the parameters derived with protein sequence (physical-chemical, energetic, and conformational properties of amino acid residues) show very weak correlation (|r| < 0.39) with folding rates of 28 two-state proteins, indicating that the sequence information alone is not sufficient to understand the folding rates of two-state proteins. However, the maximum positive correlation obtained for the properties, number of medium-range contacts, and alpha-helical tendency reveals the importance of local interactions to initiate protein folding. On the other hand, a remarkable correlation (r varies from -0.74 to -0.88) has been obtained between structural parameters (contact order, long-range order, and total contact distance) and protein folding rates. Further, we found that the secondary structure content and solvent accessibility play a marginal role in determining the folding rates of two-state proteins. Multiple regression analysis carried out with the combination of three properties, beta-strand tendency, enthalpy change, and total contact distance improved the correlation to 0.92 with protein folding rates. The relative importance of existing methods along with multiple-regression model proposed in this work will be discussed. Our results demonstrate that the native-state topology is the major determinant for the folding rates of two-state proteins.     

Slow Unfolded-State Structuring in Acyl-CoA Binding Protein Folding Revealed by Simulation and Experiment

Protein folding is a fundamental process in biology, key to understanding many human diseases. Experimentally, proteins often appear to fold via simple two- or three-state mechanisms involving mainly native-state interactions, yet recent network models built from atomistic simulations of small proteins suggest the existence of many possible metastable states and folding pathways. We reconcile these two pictures in a combined experimental and simulation study of acyl-coenzyme A binding protein (ACBP), a two-state folder (folding time ～10 ms) exhibiting residual unfolded-state structure, and a putative early folding intermediate. Using single-molecule FRET in conjunction with side-chain mutagenesis, we first demonstrate that the denatured state of ACBP at near-zero denaturant is unusually compact and enriched in long-range structure that can be perturbed by discrete hydrophobic core mutations. We then employ ultrafast laminar-flow mixing experiments to study the folding kinetics of ACBP on the microsecond time scale. These studies, along with Trp-Cys quenching measurements of unfolded-state dynamics, suggest that unfolded-state structure forms on a surprisingly slow (～100 μs) time scale, and that sequence mutations strikingly perturb both time-resolved and equilibrium smFRET measurements in a similar way. A Markov state model (MSM) of the ACBP folding reaction, constructed from over 30 ms of molecular dynamics trajectory data, predicts a complex network of metastable stables, residual unfolded-state structure, and kinetics consistent with experiment but no well-defined intermediate preceding the main folding barrier. Taken together, these experimental and simulation results suggest that the previously characterized fast kinetic phase is not due to formation of a barrier-limited intermediate but rather to a more heterogeneous and slow acquisition of unfolded-state structure.     

The Construction of New Proteins and Enzymes-a Prospect for the Future?

Only a vanishingly small proportion of the almost infinite number of possible proteins occur in nature. Can this remaining potential of structural and functional diversity be used in the construction of new proteins? Is a “second evolution” of proteins and enzymes about to occur? These questions have suddenly become of interest because the recombinant DNA technique allows the synthesis of any given amino acid sequence. Examples of enzyme models demonstrate clearly that the unusual catalytic properties of enzymes are associated with the presence of a specifically folded polypeptide chain which has a complex three-dimensional form. The critical hurdle in the path of artificial proteins is thus the design of amino acid sequences which are able to fold into tertiary structures. — Recent studies on the topology and the mechanism of folding have provided considerable insight into the occurrence of, and the rules governing the three-dimensional architecture of proteins. Secondary structures apparently play a key role in the folding process; helices and “β-structures” act as nucleation centers directing folding and account for the surprisingly small number of different folding topologies. The problem of secondary structure formation can be investigated directly by means of conformational studies on model peptides. Oligopeptides with tailormade physicochemical, structural and conformational properties can already be designed. The theoretical and experimental basis for the construction of polypeptides with stable tertiary structures is therefore established. The path to macromolecules with an immense variety of novel properties lays before us.     

Probing Tertiary Structures of Proteins by Limited Proteolysis and HPLC-MS

It is generally accepted that one or more distinct, populated intermediates are involved in the folding of most proteins. Although the understanding of the tertiary structure of the intermediate is important to elucidate the hierarchic folding pathway, it is difficult to characterize the structure because it exists only transiently during the folding process. Recently, it has been shown that several globular proteins, when placed under certain denaturing conditions, can exist in equilibrium molten globule states.     

A lattice protein with an amyloidogenic latent state: stability and folding kinetics

We have designed a model lattice protein that has two stable folded states, the lower free energy native state and a latent state of somewhat higher energy. The two states have a sizable part of their structures in common (two "alpha-helices") and differ in the content of "alpha-helices" and "beta-strands" in the rest of their structures; i.e. for the native state, this part is alpha-helical, and for the latent state it is composed of beta-strands. Thus, the lattice protein free energy surface mimics that of amyloidogenic proteins that form well organized fibrils under appropriate conditions. A Go-like potential was used and the folding process was simulated with a Monte Carlo method. To gain insight into the equilibrium free energy surface and the folding kinetics, we have combined standard approaches (reduced free energy surfaces, contact maps, time-dependent populations of the characteristic states, and folding time distributions) with a new approach. The latter is based on a principal coordinate analysis of the entire set of contacts, which makes possible the introduction of unbiased reaction coordinates and the construction of a kinetic network for the folding process. The system is found to have four characteristic basins, namely a semicompact globule, an on-pathway intermediate (the bifurcation basin), and the native and latent states. The bifurcation basin is shallow and consists of the structure common to the native and latent states, with the rest disorganized. On the basis of the simulation results, a simple kinetic model describing the transitions between the characteristic states was developed, and the rate constants for the essential transitions were estimated. During the folding process the system dwells in the bifurcation basin for a relatively short time before it proceeds to the native or latent state. We suggest that such a bifurcation may occur generally for proteins in which native and latent states have a sizable part of their structures in common. Moreover, there is the possibility of introducing changes in the system (e.g., mutations), which guide the system toward the native or misfolded state.     

Prediction of protein folding rates from primary sequences using hybrid sequence representation

The ability to predict protein folding rates constitutes an important step in understanding the overall folding mechanisms. Although many of the prediction methods are structure based, successful predictions can also be obtained from the sequence. We developed a novel method called prediction of protein folding rates (PPFR), for the prediction of protein folding rates from protein sequences. PPFR implements a linear regression model for each of the mainstream folding dynamics including two-, multi-, and mixed-state proteins. The proposed method provides predictions characterized by strong correlations with the experimental folding rates, which equal 0.87 for the two- and multistate proteins and 0.82 for the mixed-state proteins, when evaluated with out-of-sample jackknife test. Based on in-sample and out-of-sample tests, the PPFR's predictions are shown to be better than most of other sequence only and structure-based predictors and complementary to the predictions of the most recent sequence-based QRSM method. We show that simultaneous incorporation of several characteristics, including the sequence, physiochemical properties of residues, and predicted secondary structure provides improved quality. This hybridized prediction model was analyzed to reveal the complementary factors that can be used in tandem to predict folding rates. We show that bigger proteins require more time for folding, higher helical and coil content and the presence of Phe, Asn, and Gln may accelerate the folding process, the inclusion of Ile, Val, Thr, and Ser may slow down the folding process, and for the two-state proteins increased beta-strand content may decelerate the folding process. Finally, PPFR provides strong correlation when predicting sequences with low similarity.     

Using ESI-MS to probe protein structure by site-specific noncovalent attachment of 18-crown-6

A new method for probing the equilibrium structures and folding states of proteins utilizing electrospray ionization mass spectrometry is described. Protein structure is explored as a function of side-chain availability as determined by a specific interaction between lysine and 18-crown-6 ether (18C6). Various intramolecular interactions are competitive with the lysine/18C6 interaction and can prevent noncovalent attachment of 18C6. Changes to protein structure modify these inhibiting intramolecular interactions, which leads to a change in the number of 18C6s that attach to the protein. Experiments conducted with cytochrome c, ubiquitin, and melittin reveal that the method is sensitive to changes in both tertiary and secondary structure. In addition, the structure of each charge state can be examined independently. Experiments can be performed under conditions where the pH and amount of organic cosolvent are varied. Control experiments conducted with pentalysine, which lacks structural organization, are also presented.     

How does the coupling of secondary and tertiary interactions control the folding of helical macromolecules?

The authors study how the simultaneous presence of short-range secondary and long-range tertiary interactions controls the folding and collapse behavior of a helical macromolecule. The secondary interactions stabilize the helical conformation of the chain, while the tertiary interactions govern its overall three-dimensional shape. The authors have carried out Monte Carlo simulations to study the effect of chain length on the folding and collapse behavior of the chain. They have calculated state diagrams for four chain lengths and found that the physics is very rich with a plethora of stable conformational states. In addition to the helix-coil and coil-globule transitions, their model describes the coupling between them which takes place at low temperatures. Under these conditions, their model predicts a cascade of continuous, conformational transitions between states with an increase in the strength of the tertiary interactions. During each transition the chain shrinks, i.e., collapses, in a rapid and specific manner. In addition, the number of the transitions increases with increasing chain length. They have also found that the low-temperature regions of the state diagram between the transition lines cannot be associated with specific structures of the chain, but rather, with ensembles of various configurations of the chain with similar characteristics. Based on these results the authors propose a mechanism for the folding and collapse of helical macromolecules which is further supported by the analysis of configurational, configurational, and thermodynamic properties of the chain.     

Folding dynamics of Trp-cage in the presence of chemical interference and macromolecular crowding. I

Proteins fold and function in the crowded environment of the cell's interior. In the recent years it has been well established that the so-called "macromolecular crowding" effect enhances the folding stability of proteins by destabilizing their unfolded states for selected proteins. On the other hand, chemical and thermal denaturation is often used in experiments as a tool to destabilize a protein by populating the unfolded states when probing its folding landscape and thermodynamic properties. However, little is known about the complicated effects of these synergistic perturbations acting on the kinetic properties of proteins, particularly when large structural fluctuations, such as protein folding, have been involved. In this study, we have first investigated the folding mechanism of Trp-cage dependent on urea concentration by coarse-grained molecular simulations where the impact of urea is implemented into an energy function of the side chain and/or backbone interactions derived from the all-atomistic molecular dynamics simulations with urea through a Boltzmann inversion method. In urea solution, the folding rates of a model miniprotein Trp-cage decrease and the folded state slightly swells due to a lack of contact formation between side chains at the terminal regions. In addition, the equilibrium m-values of Trp-cage from the computer simulations are in agreement with experimental measurements. We have further investigated the combined effects of urea denaturation and macromolecular crowding on Trp-cage's folding mechanism where crowding agents are modeled as hard-spheres. The enhancement of folding rates of Trp-cage is most pronounced by macromolecular crowding effect when the extended conformations of Trp-cast dominate at high urea concentration. Our study makes quantitatively testable predictions on protein folding dynamics in a complex environment involving both chemical denaturation and macromolecular crowding effects.     

The Construction of New Proteins. Part III. Artificial folding units by assembly of amphiphilic secondary structures on a template

A new general strategy for the construction of artificial proteins with predetermined tertiary structure is presented. Amphiphilic α-helix and β-sheet-forming oligopeptides are assembled on a multifunctional ttemplate molecule which directs the peptide blocks to adopt characteristic folding topologies. The design, synthesis, and conformational properties of these template-assembled synthetic proteins (TASP) are exemplified for βαβ-, α-helix-bundle- and β-barrel-like tertiary structures using specially designed oligopeptides as template molecules. In contrast to linear polypeptide chains of comparable molecular weights, these conceptually novel marcromolecules are readily accessible to chemical synthesis and exhibit excellent solubility in a number of solvents. Experimental evidence is provided for a template-induced intramolecular folding to secondary and tertiary structures in aqueous solutions. This approach opens new prospects for the chemical construction of biomacromolecules with tailormade structural and functional properties.     

Evidence of a two-stage thermal denaturation process in lysozyme: a Raman scattering and differential scanning calorimetry investigation

Raman spectroscopy (in the low-frequency range and the amide I band region) and modulated differential scanning calorimetry investigations have been used to analyze temperature-induced structural changes in lysozyme dissolved in 1H2O and 2H2O in the thermal denaturation process. Low-frequency Raman data reveal a change in tertiary structure without concomitant unfolding of the secondary structure. Calorimetric data show that this structural change is responsible for the configurational entropy change associated with the strong-to-fragile liquid transition and correspond to about 1/3 of the native-denaturated transition enthalpy. This is the first stage of the thermal denaturation which is a precursor of the secondary structure change and is determined to be strongly dependent on the stability of the hydrogen-bond network in water. Low-frequency Raman spectroscopy provides information on the flexibility of the tertiary structure (in the native state and the transient folding state) in relation to the fragility of the mixture. The unfolding of the secondary structure appears as a consequence of the change in the tertiary structure and independent of the solvent. Protein conformational stability is directly dependent on the stability of the native tertiary structure. The structural transformation of tertiary structure can be detected through the enhanced 1H/2H exchange inhibited in native proteins. Taking into account similar features reported in the literature observed for different proteins it can be considered that the two-stage transformation observed in lysozyme dissolved in water is a general mechanism for the thermal denaturation of proteins.     

Folding of elongated proteins: conventional or anomalous?

The complexity of the mechanisms by which proteins fold has been shown by many studies to be governed by their native-state topologies. This was manifested in the ability of the native topology-based model to capture folding mechanisms and the success of folding rate predictions based on various topological measures, such as the contact order. However, while the finer details of topological complexity have been thoroughly examined and related to folding kinetics, simpler characteristics of the protein, such as its overall shape, have been largely disregarded. In this study, we investigated the folding of proteins with an unusual elongated geometry that differs substantially from the common globular structure. To study the effect of the elongation degree on the folding kinetics, we used repeat proteins, which become more elongated as they include more repeating units. Some of these have apparently anomalous experimental folding kinetics, with rates that are often less than expected on the basis of rates for globular proteins possessing similar topological complexity. Using experimental folding rates and a larger set of rates obtained from simulations, we have shown that as the protein becomes increasingly elongated, its folding kinetics becomes slower and deviates more from the rate expected on the basis of topology measures fitted for globular proteins. The observed slow kinetics is a result of a more complex pathway in which stable intermediates composed of several consecutive repeats can appear. We thus propose a novel measure, an elongation-sensitive contact order, that takes into account both the extent of elongation and the topological complexity of the protein. This new measure resolves the apparent discrimination between the folding of globular and elongated repeat proteins. Our study extends the current capabilities of folding-rate predictions by unifying the kinetics of repeat and globular proteins.     

A conventional an 2DCOS infrared approach to the kinetics of protein misfolding

Cell viability depends on the correct folding of the proteins involved in metabolism. Proteins are synthesized on the endoplasmic reticulum and must follow a pathway to a correct, metastable, tridimensional structure. Changes in structure or in environmental conditions can drive an instability of the folding conditions and produce non-active aggregates that in principle are proteolysed by the cellular mechanisms. However, these aggregates can be even more stable than the native proteins, escaping the cellular control. They can be classified as amorphous, if there is not a well-organized structural pattern, or ordered if a repetitive pattern is produced. These ordered structures, known as fibrils, are involved in many diseases. Infrared spectroscopy is a method of choice to study its formation because it is not affected by turbidity or the formation of high molecular weight aggregates. Moreover, in both cases, two bands characteristic of intermolecular β-sheets allow the monitoring of the aggregate formation. In both cases, the appearance of these bands involves a non-reversible path in protein folding. It has been suggested that a difference in the ordered structures involves an increasing in band intensity. This change can be the origin in variations on the 2DCOS maps. The synchronous map gives an overall idea of the process involved. The asynchronous is more informative because reflects the kinetic changes produced. The outcome of both processes, amorphous or ordered is that 2DCOS can provide a further insight to the knowledge of the kinetic processes giving rise to aggregated structures. This outcome could consist on the order in which the different secondary structures are prone to form the aggregates.     

Tertiary interactions determine the accuracy of RNA folding

RNAs must fold into unique three-dimensional structures to function in the cell, but how each polynucleotide finds its native structure is not understood. To investigate whether the stability of the tertiary structure determines the speed and accuracy of RNA folding, docking of a tetraloop with its receptor in a bacterial group I ribozyme was perturbed by site-directed mutagenesis. Disruption of the tetraloop or its receptor destabilizes tertiary interactions throughout the ribozyme by 2-3 kcal/mol, demonstrating that tertiary interactions form cooperatively in the transition from a native-like intermediate to the native state. Nondenaturing PAGE and RNase T1 digestion showed that base pairs form less homogeneously in the mutant RNAs during the transition from the unfolded state to the intermediate. Thus, tertiary interactions between helices bias the ensemble of secondary structures toward native-like conformations. Time-resolved hydroxyl radical footprinting showed that the wild-type ribozyme folds completely within 5-20 ms. By contrast, only 40-60% of a tetraloop mutant ribozyme folds in 30-40 ms, with the remainder folding in 30-200 s via nonnative intermediates. Therefore, destabilization of tetraloop-receptor docking introduces an alternate folding pathway in the otherwise smooth energy landscape of the wild-type ribozyme. Our results show that stable tertiary structure increases the flux through folding pathways that lead directly and rapidly to the native structure.     

Trigger factor slows co-translational folding through kinetic trapping while sterically protecting the nascent chain from aberrant cytosolic interactions

The E. coli chaperone trigger factor (TF) interacts directly with nascent polypeptide chains as they emerge from the ribosome exit tunnel. Small protein domains can fold under the cradle created by TF, but the co-translational folding of larger proteins is slowed down by its presence. Because of the great experimental challenges in achieving high spatial and time resolution, it is not yet known whether or not TF alters the folding properties of small proteins and if the reduced rate of folding of larger proteins is the result of kinetic or thermodynamic effects. We show, by molecular simulations employing a coarse-grained model of a series of ribosome nascent-chain complexes, that TF does not alter significantly the co-translational folding process of a small protein G domain but delays that of a large β-galactosidase domain as a result of kinetic trapping of its unfolded ensemble. We demonstrate that this trapping occurs through a combination of three distinct mechanisms: a decrease in the rate of structural rearrangements within the nascent chain, an increase in the effective exit tunnel length due to folding outside the cradle, and entanglement of the nascent chain with TF. We present evidence that this TF-induced trapping represents a trade-off between promoting co-translational folding and sterically shielding the nascent chain from aberrant cytosolic interactions that could lead to its aggregation or degradation.     

Multi-nucleation and vectorial folding pathways of large helix protein

At present we have already had the detailed knowledge of the folding of small model proteins, but a unified picture of how large proteins fold is still absent. We simulated the folding of a large eight-helix-bundle protein with a length of 145 amino acids by using a united-residue protein model. We observed a multiple nucleation folding pathway: the formation of secondary structures was followed by the nucleation of helices at the two terminal parts and also at the middle of the chain, and then the nuclei grew and combined with each other to form the tertiary structure. Surprisingly, we also found a vectorial folding pathway that was shown recently for co-translational folding in the ribosome exit tunnel. Furthermore, we found that all three-helix subunits in the chain can fold into native-like conformations independently, especially those at the two terminal parts and the middle of the chain, which may be responsible for the nucleation's. These results may be helpful to understand the folding mechanism of large repeat helical proteins.