Free energy landscape and multiple folding pathways of an H-type RNA pseudoknot.
Ontology highlight
ABSTRACT:
SUBMITTER: Bian Y
PROVIDER: S-EPMC4451515 | biostudies-literature | 2015
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The functions of RNA pseudoknots (PKs), which are minimal tertiary structural motifs and an integral part of several ribozymes and ribonucleoprotein complexes, are determined by their structure, stability, and dynamics. Therefore, it is important to elucidate the general principles governing their thermodynamics/folding mechanisms. Here, we combine laser temperature-jump experiments and coarse-grained simulations to determine the folding/unfolding pathways of VPK, a variant of the mouse mammary tumor virus (MMTV) PK involved in ribosomal frameshifting. Fluorescent nucleotide analogs (2-aminopurine and pyrrolocytidine) placed at different stem/loop positions in the PK serve as local probes allowing us to monitor the order of assembly of VPK that has two constituent hairpins with different intrinsic stabilities. We show that at 50 mM KCl, the dominant folding pathway populates only the more stable hairpin intermediate; as the salt concentration is increased, a parallel folding pathway emerges involving the less stable hairpin as an alternate intermediate. Notably, the flux between the pathways is modulated by the ionic strength. Our findings support the principle that the order of PK structure formation is determined by the relative stabilities of the hairpins, which can be altered by sequence variations or salt concentrations. The experimental results of salt effects on the partitioning between the two folding pathways are in remarkable agreement with simulations that were performed with no adjustable parameters. Our study not only unambiguously demonstrates that VPK folds by parallel pathways but also showcases the power of combining experiments and simulations for a more enriched description of RNA self-assembly.
| S-EPMC6077692 | biostudies-other
Project description:Based on the experimentally determined atomic coordinates for RNA helices and the self-avoiding walks of the P (phosphate) and C4 (carbon) atoms in the diamond lattice for the polynucleotide loop conformations, we derive a set of conformational entropy parameters for RNA pseudoknots. Based on the entropy parameters, we develop a folding thermodynamics model that enables us to compute the sequence-specific RNA pseudoknot folding free energy landscape and thermodynamics. The model is validated through extensive experimental tests both for the native structures and for the folding thermodynamics. The model predicts strong sequence-dependent helix-loop competitions in the pseudoknot stability and the resultant conformational switches between different hairpin and pseudoknot structures. For instance, for the pseudoknot domain of human telomerase RNA, a native-like and a misfolded hairpin intermediates are found to coexist on the (equilibrium) folding pathways, and the interplay between the stabilities of these intermediates causes the conformational switch that may underlie a human telomerase disease.
| S-EPMC1463895 | biostudies-literature
Project description:There are potentially several ways Mg2+ might promote formation of an RNA tertiary structure: by causing a general "collapse" of the unfolded ensemble to more compact conformations, by favoring a reorganization of structure within a domain to a form with specific tertiary contacts, and by enhancing cooperative linkages between different sets of tertiary contacts. To distinguish these different modes of action, we have studied Mg2+ interactions with the adenine riboswitch, in which a set of tertiary interactions that forms around a purine-binding pocket is thermodynamically linked to the tertiary "docking" of two hairpin loops in another part of the molecule. Each of four RNA forms with different extents of tertiary structure were characterized by small-angle X-ray scattering. The free energy of interconversion between different conformations in the absence of Mg2+ and the free energy of Mg2+ interaction with each form have been estimated, yielding a complete picture of the folding energy landscape as a function of Mg2+ concentration. At 1 mM Mg2+ (50 mM K+), the overall free energy of stabilization by Mg2+ is large, -9.8 kcal/mol, and about equally divided between its effect on RNA collapse to a partially folded structure and on organization of the binding pocket. A strong cooperative linkage between the two sets of tertiary contacts is intrinsic to the RNA. This quantitation of the effects of Mg2+ on an RNA with two distinct sets of tertiary interactions suggests ways that Mg2+ may work to stabilize larger and more complex RNA structures.
| S-EPMC3507322 | biostudies-literature
Project description:In spite of the success of computational methods for predicting RNA secondary structure, the problem of predicting RNA tertiary structure folding remains. Low-resolution structural models show promise as they allow for rigorous statistical mechanical computation for the conformational entropies, free energies, and the coarse-grained structures of tertiary folds. Molecular dynamics refinement of coarse-grained structures leads to all-atom 3D structures. Modeling based on statistical mechanics principles also has the unique advantage of predicting the full free energy landscape, including local minima and the global free energy minimum. The energy landscapes combined with the 3D structures form the basis for quantitative predictions of RNA functions. In this chapter, we present an overview of statistical mechanical models for RNA folding and then focus on a recently developed RNA statistical mechanical model -- the Vfold model. The main emphasis is placed on the physics underpinning the models, the computational strategies, and the connections to RNA biology.
| S-EPMC4902161 | biostudies-literature
Project description:Trp-cage is a 20-residue miniprotein, which is believed to be the fastest folder known so far. In this study, the folding free energy landscape of Trp-cage has been explored in explicit solvent by using an OPLSAA force field with periodic boundary condition. A highly parallel replica exchange molecular dynamics method is used for the conformation space sampling, with the help of a recently developed efficient molecular dynamics algorithm P3ME/RESPA (particle-particle particle-mesh Ewald/reference system propagator algorithm). A two-step folding mechanism is proposed that involves an intermediate state where two correctly formed partial hydrophobic cores are separated by an essential salt-bridge between residues Asp-9 and Arg-16 near the center of the peptide. This metastable intermediate state provides an explanation for the superfast folding process. The free energy landscape is found to be rugged at low temperatures, and then becomes smooth and funnel-like above 340 K. The lowest free energy structure at 300 K is only 1.50 A Calpha-RMSD (Calpha-rms deviation) from the NMR structures. The simulated nuclear Overhauser effect pair distances are in excellent agreement with the raw NMR data. The temperature dependence of the Trp-cage population, however, is found to be significantly different from experiment, with a much higher melting transition temperature above 400 K (experimental 315 K), indicating that the current force fields, parameterized at room temperature, need to be improved to correctly predict the temperature dependence.
| S-EPMC263783 | biostudies-literature
Project description:The dynamic mechanisms by which RNAs acquire biologically functional structures are of increasing importance to the rapidly expanding fields of RNA therapeutics and biotechnology. Large energy barriers separating misfolded and functional states arising from alternate base pairing are a well-appreciated characteristic of RNA. In contrast, it is typically assumed that functionally folded RNA occupies a single native basin of attraction that is free of deeply dividing energy barriers (ergodic hypothesis). This assumption is widely used as an implicit basis to interpret experimental ensemble-averaged data. Here, we develop an experimental approach to isolate persistent sub-populations of a small RNA enzyme and show by single molecule fluorescence resonance energy transfer (smFRET), biochemical probing and high-resolution mass spectrometry that commitment to one of several catalytically active folds occurs unexpectedly high on the RNA folding energy landscape, resulting in partially irreversible folding. Our experiments reveal the retention of molecular heterogeneity following the complete loss of all native secondary and tertiary structure. Our results demonstrate a surprising longevity of molecular heterogeneity and advance our current understanding beyond that of non-functional misfolds of RNA kinetically trapped on a rugged folding-free energy landscape.
| S-EPMC2602785 | biostudies-literature
Project description:Folding funnel is the essential concept of the free energy landscape for ordered proteins. How does this concept apply to intrinsically disordered proteins (IDPs)? Here, we address this fundamental question through the explicit characterization of the free energy landscapes of the representative α-helical (HP-35) and β-sheet (WW domain) proteins and of an IDP (pKID) that folds upon binding to its partner (KIX). We demonstrate that HP-35 and WW domain indeed exhibit the steep folding funnel: the landscape slope for these proteins is ca. -50 kcal/mol, meaning that the free energy decreases by ~5 kcal/mol upon the formation of 10% native contacts. On the other hand, the landscape of pKID is funneled but considerably shallower (slope of -24 kcal/mol), which explains why pKID is disordered in free environments. Upon binding to KIX, the landscape of pKID now becomes significantly steep (slope of -54 kcal/mol), which enables otherwise disordered pKID to fold. We also show that it is the pKID-KIX intermolecular interactions originating from hydrophobic residues that mainly confer the steep folding funnel. The present work not only provides the quantitative characterization of the protein folding free energy landscape, but also establishes the usefulness of the folding funnel concept to IDPs.
| S-EPMC6797787 | biostudies-literature
Project description:Protein folding is a central problem in biological physics. Energetic roughness is an important aspect that controls protein-folding stability and kinetics. The roughness is associated with conflicting interactions in the protein and is also known as frustration. Recent studies indicate that an addition of a small amount of energetic frustration may enhance folding speed for certain proteins. In this study, we have investigated the conditions under which frustration increases the folding rate. We used a Cα structure-based model to simulate a group of proteins. We found that the free-energy barrier at the transition state (ΔF) correlates with nonnative-contact variation (ΔA), and the simulated proteins are clustered according to their fold motifs. These findings are corroborated by the Clementi-Plotkin analytical model. As a consequence, the optimum frustration regime for protein folding can be predicted analytically.
| S-EPMC4968343 | biostudies-literature
Project description:A complete description of the pathways and mechanisms of protein folding requires a detailed structural and energetic characterization of the conformational ensemble along the entire folding reaction coordinate. Simulations can provide this level of insight for small proteins. In contrast, with the exception of hydrogen exchange, which does not monitor folding directly, experimental studies of protein folding have not yielded such structural and energetic detail. NMR can provide residue specific atomic level structural information, but its implementation in protein folding studies using chemical or temperature perturbation is problematic. Here we present a highly detailed structural and energetic map of the entire folding landscape of the leucine-rich repeat protein, pp32 (Anp32), obtained by combining pressure-dependent site-specific 1H-15N HSQC data with coarse-grained molecular dynamics simulations. The results obtained using this equilibrium approach demonstrate that the main barrier to folding of pp32 is quite broad and lies near the unfolded state, with structure apparent only in the C-terminal region. Significant deviation from two-state unfolding under pressure reveals an intermediate on the folded side of the main barrier in which the N-terminal region is disordered. A nonlinear temperature dependence of the population of this intermediate suggests a large heat capacity change associated with its formation. The combination of pressure, which favors the population of folding intermediates relative to chemical denaturants; NMR, which allows their observation; and constrained structure-based simulations yield unparalleled insight into protein folding mechanisms.
| S-EPMC5153537 | biostudies-literature
Project description:The folding free energy landscape of the C-terminal beta hairpin of protein G has been explored in this study with explicit solvent under periodic boundary condition and OPLSAA force field. A highly parallel replica exchange method that combines molecular dynamics trajectories with a temperature exchange Monte Carlo process is used for sampling with the help of a new efficient algorithm P3ME/RESPA. The simulation results show that the hydrophobic core and the beta strand hydrogen bond form at roughly the same time. The free energy landscape with respect to various reaction coordinates is found to be rugged at low temperatures and becomes a smooth funnel-like landscape at about 360 K. In contrast to some very recent studies, no significant helical content has been found in our simulation at all temperatures studied. The beta hairpin population and hydrogen-bond probability are in reasonable agreement with the experiment at biological temperature, but both decay more slowly than the experiment with temperature.
| S-EPMC64961 | biostudies-literature