Understanding the Role of Three-Dimensional Topology in Determining the?Folding Intermediates of Group I Introns

Many RNA molecules exert their biological function only after folding to unique three-dimensional structures. For long, noncoding RNA molecules, the complexity of finding the native topology can be a major impediment to correct folding to the biologically active structure. An RNA molecule may fold to a near-native structure but not be able to continue to the correct structure due to a topological barrier such as crossed strands or incorrectly stacked helices. Achieving the native conformation thus requires unfolding and refolding, resulting in a long-lived intermediate. We investigate the role of topology in the folding of two phylogenetically related catalytic group I introns, the Twort and Azoarcus group I ribozymes. The kinetic models describing the Mg²⁺-mediated folding of these ribozymes were previously determined by time-resolved hydroxyl (⋅OH) radical footprinting. Two intermediates formed by parallel intermediates were resolved for each RNA. These data and analytical ultracentrifugation compaction analyses are used herein to constrain coarse-grained models of these folding intermediates as we investigate the role of nonnative topology in dictating the lifetime of the intermediates. Starting from an ensemble of unfolded conformations, we folded the RNA molecules by progressively adding native constraints to subdomains of the RNA defined by the ⋅OH time-progress curves to simulate folding through the different kinetic pathways. We find that nonnative topologies (arrangement of helices) occur frequently in the folding simulations despite using only native constraints to drive the reaction, and that the initial conformation, rather than the folding pathway, is the major determinant of whether the RNA adopts nonnative topology during folding. From these analyses we conclude that biases in the initial conformation likely determine the relative flux through parallel RNA folding pathways.     

Structural Rearrangements Linked to Global Folding Pathways of the Azoarcus Group I Ribozyme

Stable RNAs must fold into specific three-dimensional structures to be biologically active, yet many RNAs form metastable structures that compete with the native state. Our previous time-resolved footprinting experiments showed that Azoarcus group I ribozyme forms its tertiary structure rapidly (τ < 30 ms) without becoming significantly trapped in kinetic intermediates. Here, we use stopped-flow fluorescence spectroscopy to probe the global folding kinetics of a ribozyme containing 2-aminopurine in the loop of P9. The modified ribozyme was catalytically active and exhibited two equilibrium folding transitions centered at 0.3 and 1.6 mM Mg²⁺, consistent with previous results. Stopped-flow fluorescence revealed four kinetic folding transitions with observed rate constants of 100, 34, 1, and 0.1 s^− 1 at 37 °C. From comparison with time-resolved Fe(II)-ethylenediaminetetraacetic acid footprinting of the modified ribozyme under the same conditions, these folding transitions were assigned to formation of the I_C intermediate, tertiary folding and docking of the nicked P9 tetraloop, reorganization of the P3 pseudoknot, and refolding of nonnative conformers, respectively. The footprinting results show that 50-60% of the modified ribozyme folds in less than 30 ms, while the rest of the RNA population undergoes slow structural rearrangements that control the global folding rate. The results show how small perturbations to the structure of the RNA, such as a nick in P9, populate kinetic folding intermediates that are not observed in the natural ribozyme.     

Structure of a folding intermediate reveals the interplay between core and peripheral elements in RNA folding

Though the molecular architecture of many native RNA structures has been characterized, the structures of folding intermediates are poorly defined. Here, we present a nucleotide-level model of a highly structured equilibrium folding intermediate of the specificity domain of the Bacillus subtilis RNase P RNA, obtained using chemical and nuclease mapping, circular dichroism spectroscopy, small-angle X-ray scattering and molecular modeling. The crystal structure indicates that the 154 nucleotide specificity domain is composed of several secondary and tertiary structural modules. The structure of the intermediate contains modules composed of secondary structures and short-range tertiary interactions, implying a sequential order of tertiary structure formation during folding. The intermediate lacks the native core and several long-range interactions among peripheral regions, such as a GAAA tetraloop and its receptor. Folding to the native structure requires the local rearrangement of a T-loop in the core in concert with the formation of the GAAA tetraloop-receptor interaction. The interplay of core and peripheral structure formation rationalizes the high degree of cooperativity observed in the folding transition leading to the native structure.     

Bioactivity of folding intermediates studied by the recovery of enzymatic activity during refolding

The peptide bond preceding proline residues realizes a cis/trans conformational switch with high switching resistance in native proteins and folding intermediates. Therefore, individual isomers have the potential to differ in bioactivity. However, information about isomer-specific bioactivities is difficult to obtain because of the risk of affecting isomeric distribution by bioactivity assay components.Here we present an approach that allows for the measurement of the recovery of enzymatic activities of wild-type RNase T₁ and RNase T₁ variants during refolding under conditions where the population of enzyme-substrate or enzyme-product complexes is negligible. Recovery of enzymatic activity was continuously monitored within the visible range of the spectrum by addition of a fluorescence-labeled nucleotide substrate to the refolding sample. We found that a nonnative trans conformation at Pro39 renders the RNase T₁ almost completely inactive. A folding intermediate having a nonnative trans conformation at Pro55 shows about 46% of the enzymatic activity referred to the native state. Pro55, in contrast to the active site located Pro39, is situated in a solvent-exposed loop region remote from active-site residues. In both cases, peptidyl prolyl cis/trans isomerases accelerate the regain of nucleolytic activity. Our findings show that even if there is a considerable distance between the site of isomerization and the active site, conformational control of the bioactivity of proteins is likely to occur, and that the surface location of prolyl bonds suffices for the control of buried active sites mediated by peptidyl prolyl cis/trans isomerases.     

Confirmation of the hierarchical folding of RNase H: a protein engineering study.

The kinetic intermediate of RNase H is structured in a core region of the protein. To probe the role of this intermediate in the folding of RNase H, the folding kinetics of mutant proteins with altered native state stabilities were investigated. Mutations within the folding core destabilize the kinetic intermediate and slow refolding in a manner consistent with an obligatory intermediate model. Mutations outside of the folding core, however, do not affect the stability of the kinetic intermediate but do perturb the native state and transition state. These results indicate that interactions formed in the intermediate persist in the transition and native states and that RNase H folds through a hierarchical mechanism.     

Protein Roles in Group I Intron RNA Folding: The Tyrosyl-tRNA Synthetase CYT-18 Stabilizes the Native State Relative to a Long-Lived Misfolded Structure without Compromising Folding Kinetics

The Neurospora crassa CYT-18 protein is a mitochondrial tyrosyl-tRNA synthetase that also promotes self-splicing of group I intron RNAs by stabilizing the functional structure in the conserved core. CYT-18 binds the core along the same surface as a common peripheral element, P5abc, suggesting that CYT-18 can replace P5abc functionally. In addition to stabilizing structure generally, P5abc stabilizes the native conformation of the Tetrahymena group I intron relative to a globally similar misfolded conformation that has only local differences within the core and is populated significantly at equilibrium by a ribozyme variant lacking P5abc (E^ΔP5abc). Here, we show that CYT-18 specifically promotes formation of the native group I intron core from this misfolded conformation. Catalytic activity assays demonstrate that CYT-18 shifts the equilibrium of E^ΔP5abc toward the native state by at least 35-fold, and binding assays suggest an even larger effect. Thus, similar to P5abc, CYT-18 preferentially recognizes the native core, despite the global similarity of the misfolded core and despite forming crudely similar complexes, as revealed by dimethyl sulfate footprinting. Interestingly, the effects of CYT-18 and P5abc on folding kinetics differ. Whereas P5abc inhibits refolding of the misfolded conformation by forming peripheral contacts that must break during refolding, CYT-18 does not display analogous inhibition, most likely because it relies to a greater extent on direct interactions with the core. Although CYT-18 does not encounter this RNA in vivo, our results suggest that it stabilizes its cognate group I introns relative to analogous misfolded intermediates. By specifically recognizing native structural features, CYT-18 may also interact with earlier folding intermediates to avoid RNA misfolding or to trap native contacts as they form. More generally, our results highlight the ability of a protein cofactor to stabilize a functional RNA structure specifically without incurring associated costs in RNA folding kinetics.     

Autonomously Folding Protein Fragments Reveal Differences in the Energy Landscapes of Homologous RNases H

An important approach to understanding how a protein sequence encodes its energy landscape is to compare proteins with different sequences that fold to the same general native structure. In this work, we compare E. coli and T. thermophilus homologs of the protein RNase H. Using protein fragments, we create equilibrium mimics of two different potential partially-folded intermediates (I_core and I_core+1) hypothesized to be present on the energy landscapes of these two proteins. We observe that both T. thermophilus RNase H (ttRNH) fragments are folded and have distinct stabilities, indicating that both regions are capable of autonomous folding and that both intermediates are present as local minima on the ttRNH energy landscape. In contrast, the two E. coli RNase H (ecRNH) fragments have very similar stabilities, suggesting that the presence of additional residues in the I_core+1 fragment does not affect the folding or structure as compared to I_core. NMR experiments provide additional evidence that only the I_core intermediate is populated by ecRNH. This is one of the biggest differences that has been observed between the energy landscapes of these two proteins. Additionally, we used a FRET experiment in the background of full-length ttRNH to specifically monitor the formation of the I_core+1 intermediate. We determine that the ttRNH I_core+1 intermediate is likely the intermediate populated prior to the rate-limiting barrier to global folding, in contrast to E. coli RNase H for which I_core is the folding intermediate. This result provides new insight into the nature of the rate-limiting barrier for the folding of RNase H.     

The Folding Trajectory of RNase H Is Dominated by Its Topology and Not Local Stability: A Protein Engineering Study of Variants that Fold via Two-State and Three-State Mechanisms

Proteins can sample a variety of partially folded conformations during the transition between the unfolded and native states. Some proteins never significantly populate these high-energy states and fold by an apparently two-state process. However, many proteins populate detectable, partially folded forms during the folding process. The role of such intermediates is a matter of considerable debate. A single amino acid change can convert Escherichia coli ribonuclease H from a three-state folder that populates a kinetic intermediate to one that folds in an apparent two-state fashion. We have compared the folding trajectories of the three-state RNase H and the two-state RNase H, proteins with the same native-state topology but altered regional stability, using a protein engineering approach. Our data suggest that both versions of RNase H fold through a similar trajectory with similar high-energy conformations. Mutations in the core and the periphery of the protein affect similar aspects of folding for both variants, suggesting a common trajectory with folding of the core region followed by the folding of the periphery. Our results suggest that formation of specific partially folded conformations may be a general feature of protein folding that can promote, rather than hinder, efficient folding.     

Fast-Folding Pathways of the Thrombin-Binding Aptamer G-Quadruplex Revealed by a Markov State Model

G-quadruplex structures participate in many important cellular processes. For a better understanding of their functions, knowledge of the mechanism by which they fold into the functional native structures is necessary. In this work, we studied the folding process of the thrombin-binding aptamer G-quadruplex. Enabled by a computational paradigm that couples an advanced sampling method and a Markov state model, four folding intermediates were identified, including an antiparallel G-hairpin, two G-triplex structures, and a double-hairpin conformation. Likewise, a misfolded structure with a nonnative distribution of syn/anti guanines was also observed. Based on these states, a transition path analysis revealed three fast-folding pathways, along which the thrombin-binding aptamer would fold to the native state directly, with no evidence of potential nonnative competing conformations. The results also showed that the TGT-loop plays an important role in the folding process. The findings of this research may provide general insight about the folding of other G-quadruplex structures.     

Folding Pathways of the Tetrahymena Ribozyme

Like many structured RNAs, the Tetrahymena group I intron ribozyme folds through multiple pathways and intermediates. Under standard conditions in vitro, a small fraction reaches the native state (N) with k_obs ≈ 0.6 min^− 1, while the remainder forms a long-lived misfolded conformation (M) thought to differ in topology. These alternative outcomes reflect a pathway that branches late in folding, after disruption of a trapped intermediate (I_trap). Here we use catalytic activity to probe the folding transitions from I_trap to the native and misfolded states. We show that mutations predicted to weaken the core helix P3 do not increase the rate of folding from I_trap but they increase the fraction that reaches the native state rather than forming the misfolded state. Thus, P3 is disrupted during folding to the native state but not to the misfolded state, and P3 disruption occurs after the rate-limiting step. Interestingly, P3-strengthening mutants also increase native folding. Additional experiments show that these mutants are rapidly committed to folding to the native state, although they reach the native state with approximately the same rate constant as the wild-type ribozyme (~ 1 min^− 1). Thus, the P3-strengthening mutants populate a distinct pathway that includes at least one intermediate but avoids the M state, most likely because P3 and the correct topology are formed early. Our results highlight multiple pathways in RNA folding and illustrate how kinetic competitions between rapid events can have long-lasting effects because the “choice” is enforced by energy barriers that grow larger as folding progresses.     

Mapping metal-binding sites in the catalytic domain of bacterial RNase P RNA

Ribonuclease P (RNase P) is a ribonucleoprotein enzyme that contains a universally conserved, catalytically active RNA component. RNase P RNA requires divalent metal ions for folding, substrate binding, and catalysis. Despite recent advances in understanding the structure of RNase P RNA, no comprehensive analysis of metal-binding sites has been reported, in part due to the poor crystallization properties of this large RNA. We have developed an abbreviated yet still catalytic construct, Bst P7Δ RNA, which contains the catalytic domain of the bacterial RNase P RNA and has improved crystallization properties. We use this mutant RNA as well as the native RNA to map metal-binding sites in the catalytic core of the bacterial RNase P RNA, by anomalous scattering in diffraction analysis. The results provide insight into the interplay between RNA structure and focalization of metal ions, and a structural basis for some previous biochemical observations with RNase P. We use electrostatic calculations to extract the potential functional significance of these metal-binding sites with respect to binding Mg²⁺. The results suggest that with at least one important exception of specific binding, these sites mainly map areas of diffuse association of magnesium ions.     

The “Pre-Molten Globule,” a New Intermediate in Protein Folding

In vitro folding studies of several proteins revealed the formation, within 2–4 msec, of transient intermediates with a large far-UV ellipticity but no amide proton protection. To solve the contradiction between the secondary structure contents estimated by these two methods, we characterized the isolated C-terminal fragment F2 of the tryptophan synthase β₂ subunit. In β₂, F2 forms its tertiary interactions with the F1 N-terminal region. Hence, in the absence of F1, isolated F2 should remain at an early folding stage with no long-range interactions. We shall show that isolated F2 folds into, and remains in, a “state” called the pre-molten globule, that indeed corresponds to a 2- to 4-msec intermediate. This condensed, but not compact, “state” corresponds to an array of conformations in rapid equilibrium comprising native as well as nonnative secondary structures. It fits the “new view” on the folding process.     

A Single Mutation at Residue 25 Populates the Folding Intermediate of E. coli RNase H and Reveals a Highly Dynamic Partially Folded Ensemble

Understanding the nature of partially folded intermediates transiently populated during protein folding is important for understanding both protein folding and misfolding. These ephemeral species, however, often elude direct experimental characterization. The well-characterized protein ribonuclease H (RNase H) from Escherichia coli populates an on-pathway intermediate identified in both bulk studies and single-molecule mechanical unfolding experiments. Here, we set out to trap the transient intermediate of RNase H at equilibrium by selectively destabilizing the region of the protein known to be unfolded in this species. Surprisingly, a single change at Ile25 (I25A) resulted in the equilibrium population of the intermediate under near-native conditions. The intermediate was undetectable in a series of heteronuclear single quantum coherences, revealing the dynamic nature of this partially unfolded form on the timescale of NMR detection. This result is in contrast to studies in which the structures of trapped intermediates are solved by NMR, indicating that they are well packed and native-like. The dynamic nature of the RNase H intermediate may be important for its role as an on-pathway, productive species that promotes efficient folding.     

Nonnative Electrostatic Interactions Can Modulate Protein Folding: Molecular Dynamics with a Grain of Salt

In recent years, a growing number of protein folding studies have focused on the unfolded state, which is now recognized as playing a major role in the folding process. Some of these studies show that interactions occurring in the unfolded state can significantly affect the stability and kinetics of the protein folding reaction. In this study, we modeled the effect of electrostatic interactions, both native and nonnative, on the folding of three protein systems that underwent selective charge neutralization or reversal or complete charge suppression. In the case of the N-terminal L9 protein domain, our results directly attribute the increase in thermodynamic stability to destabilization of the unfolded ensemble, reaffirming the experimental observations. These results provide a deeper structural insight into the ensemble of the unfolded state and predict a new mutation site for increased protein stability. In the second case, charge reversal mutations of RNase Sa affected protein stability, with the destabilizing mutations being less destabilizing at higher salt concentrations, indicating the formation of charge-charge interactions in the unfolded state. In the N-terminal L9 and RNase Sa systems, changes in electrostatic interactions in the unfolded state that cause an increase in free energy had an overall compaction effect that suggests a decrease in entropy. In the third case, in which we compared the β-lactalbumin and hen egg-white lysozyme protein homologues, we successfully eliminated differences between the folding kinetics of the two systems by suppressing electrostatic interactions, supporting previously reported findings. Our coarse-grained molecular dynamics study not only reproduces experimentally reported findings but also provides a detailed molecular understanding of the elusive unfolded-state ensemble and how charge-charge interactions can modulate the biophysical characteristics of folding.     

Multiple Unfolding Events during Native Folding of the Tetrahymena Group I Ribozyme

Despite the ubiquitous nature of misfolded intermediates in RNA folding, little is known about their physical properties or the folding transitions that allow them to continue folding productively. Folding of the Tetrahymena group I ribozyme includes sequential accumulation of two intermediates, termed I_trap and misfolded (M). Here, we probe the structure and folding transition of I_trap and compare them to those of M. Hydroxyl radical and dimethyl sulfate footprinting show that both I_trap and M are extensively structured and crudely resemble the native RNA. However, regions of the core P3-P8 domain are more exposed to solvent in I_trap than in M. I_trap rearranges to continue folding nearly 1000-fold faster than M, and urea accelerates folding of I_trap much less than M. Thus, the rate-limiting transition from I_trap requires a smaller increase in exposed surface. Mutations that disrupt peripheral tertiary contacts give large and nearly uniform increases in re-folding of M, whereas the same mutations give at most modest increases in folding from I_trap. Intriguingly, mutations within the peripheral element P5abc give 5- to 10-fold accelerations in escape from I_trap, whereas ablation of P13, which lies on the opposite surface in the native structure, near the P3-P8 domain, has no effect. Thus, the unfolding required from I_trap appears to be local, whereas the unfolding of M appears to be global. Further, the modest effects from several mutations suggest that there are multiple pathways for escape from I_trap and that escape is aided by loosening nearby native structural constraints, presumably to facilitate local movements of nucleotides or segments that have not formed native contacts. Overall, these and prior results suggest a model in which the global architecture and peripheral interactions of the RNA are achieved relatively early in folding. Multiple folding and re-folding events occur on the predominant pathway to the native state, with increasing native core interactions and cooperativity as folding progresses.     

Uncovering Specific Electrostatic Interactions in the Denatured States of Proteins

The stability and folding of proteins are modulated by energetically significant interactions in the denatured state that is in equilibrium with the native state. These interactions remain largely invisible to current experimental techniques, however, due to the sparse population and conformational heterogeneity of the denatured-state ensemble under folding conditions. Molecular dynamics simulations using physics-based force fields can in principle offer atomistic details of the denatured state. However, practical applications are plagued with the lack of rigorous means to validate microscopic information and deficiencies in force fields and solvent models. This study presents a method based on coupled titration and molecular dynamics sampling of the denatured state starting from the extended sequence under native conditions. The resulting denatured-state pK_as allow for the prediction of experimental observables such as pH- and mutation-induced stability changes. I show the capability and use of the method by investigating the electrostatic interactions in the denatured states of wild-type and K12M mutant of NTL9 protein. This study shows that the major errors in electrostatics can be identified by validating the titration properties of the fragment peptides derived from the sequence of the intact protein. Consistent with experimental evidence, our simulations show a significantly depressed pK_a for Asp⁸ in the denatured state of wild-type, which is due to a nonnative interaction between Asp⁸ and Lys¹². Interestingly, the simulation also shows a nonnative interaction between Asp⁸ and Glu⁴⁸ in the denatured state of the mutant. I believe the presented method is general and can be applied to extract and validate microscopic electrostatics of the entire folding energy landscape.     

Conformation of P22 tailspike folding and aggregation intermediates probed by monoclonal antibodies.

The partitioning of partially folded polypeptide chains between correctly folded native states and off-pathway inclusion bodies is a critical reaction in biotechnology. Multimeric partially folded intermediates, representing early stages of the aggregation pathway for the P22 tailspike protein, have been trapped in the cold and isolated by nondenaturing polyacrylamide gel electrophoresis (PAGE) (speed MA, Wang DIC, King J. 1995. Protein Sci 4:900-908). Monoclonal antibodies against tailspike chains discriminate between folding intermediates and native states (Friguet B, Djavadi-Ohaniance L, King J, Goldberg ME. 1994. J Biol Chem 269:15945-15949). Here we describe a nondenaturing Western blot procedure to probe the conformation of productive folding intermediates and off-pathway aggregation intermediates. The aggregation intermediates displayed epitopes in common with productive folding intermediates but were not recognized by antibodies against native epitopes. The nonnative epitope on the folding and aggregation intermediates was located on the partially folded N-terminus, indicating that the N-terminus remained accessible and nonnative in the aggregated state. Antibodies against native epitopes blocked folding, but the monoclonal directed against the N-terminal epitope did not, indicating that the conformation of the N-terminus is not a key determinant of the productive folding and chain association pathway.     

Structure of an archaeal homolog of the human protein complex Rpp21-Rpp29 that is a key core component for the assembly of active ribonuclease P

Ribonuclease P (RNase P) is a ribonucleoprotein complex involved in the processing of the 5′-leader sequence of precursor tRNA. Human RNase P protein subunits Rpp21 and Rpp29, which bind to each other, with catalytic RNA (H1 RNA) are sufficient for activating endonucleolytic cleavage of precursor tRNA. Here we have determined the crystal structure of the complex between the Pyrococcus horikoshii RNase P proteins PhoRpp21 and PhoRpp29, the archaeal homologs of Rpp21 and Rpp29, respectively. PhoRpp21 and PhoRpp29 form a heterodimeric structure where the two N-terminal helices (α1 and α2) in PhoRpp21 predominantly interact with the N-terminal extended structure, the β-strand (β2), and the C-terminal helix (α3) in PhoRpp29. The interface is dominated by hydrogen bonds and several salt bridges, rather than hydrophobic interactions. The electrostatic potential on the surface of the heterodimer shows a positively charged cluster on one face, suggesting a possible RNA-binding surface of the PhoRpp21-PhoRpp29 complex. The present structure, along with the result of a mutational analysis, suggests that heterodimerization between PhoRpp21 and PhoRpp29 plays an important role in the function of P. horikoshii RNase P.     

The Folding Pathway of the Antibody VL Domain

Antibodies are modular proteins consisting of domains that exhibit a β-sandwich structure, the so-called immunoglobulin fold. Despite structural similarity, differences in folding and stability exist between different domains. In particular, the variable domain of the light chain V_L is unusual as it is associated with misfolding diseases, including the pathologic assembly of the protein into fibrillar structures. Here, we have analysed the folding pathway of a V_L domain with a view to determine features that may influence the relationship between productive folding and fibril formation. The V_L domain from MAK33 (murine monoclonal antibody of the subtype κ/IgG1) has not previously been associated with fibrillisation but is shown here to be capable of forming fibrils. The folding pathway of this V_L domain is complex, involving two intermediates in different pathways. An obligatory early molten globule-like intermediate with secondary structure but only loose tertiary interactions is inferred. The native state can then be formed directly from this intermediate in a phase that can be accelerated by the addition of prolyl isomerases. However, an alternative pathway involving a second, more native-like intermediate is also significantly populated. Thus, the protein can reach the native state via two distinct folding pathways. Comparisons to the folding pathways of other antibody domains reveal similarities in the folding pathways; however, in detail, the folding of the V_L domain is striking, with two intermediates populated on different branches of the folding pathway, one of which could provide an entry point for molecules diverted into the amyloid pathway.