Time-resolved NMR studies of RNA folding

Effects of mRNA conformational switching on translational noise in gene circuits

Intragenic translational heterogeneity describes the variation in translation at the level of transcripts for an individual gene. A factor that contributes to this source of variation is the mRNA structure. Both the composition of the thermodynamic ensemble, i.e., the stationary distribution of mRNA structures, and the switching dynamics between those play a role. The effect of the switching dynamics on intragenic translational heterogeneity remains poorly understood. We present a stochastic translation model that accounts for mRNA structure switching and is derived from a Markov model via approximate stochastic filtering. We assess the approximation on various timescales and provide a method to quantify how mRNA structure dynamics contributes to translational heterogeneity. With our approach, we allow quantitative information on mRNA switching from biophysical experiments or coarse-grain molecular dynamics simulations of mRNA structures to be included in gene regulatory chemical reaction network models without an increase in the number of species. Thereby, our model bridges a gap between mRNA structure kinetics and gene expression models, which we hope will further improve our understanding of gene regulatory networks and facilitate genetic circuit design.

Recognition-Enabled Automated Analyte Identification via 19F NMR

Nuclear magnetic resonance (NMR) is an indispensable tool for structural elucidation and noninvasive analysis. Automated identification of analytes with NMR is highly pursued in metabolism research and disease diagnosis; however, this process is often complicated by the signal overlap and the sample matrix. We herein report a detection scheme based on ¹⁹F NMR spectroscopy and dynamic recognition, which effectively simplifies the detection signal and mitigates the influence of the matrix on the detection. It is demonstrated that this approach can not only detect and differentiate capsaicin and dihydrocapsaicin in complex real-world samples but also quantify the ibuprofen content in sustained-release capsules. Based on the ¹⁹F signals obtained in the detection using a set of three ¹⁹F probes, automated analyte identification is achieved, effectively reducing the odds of misrecognition caused by structural similarity.

The RNA chaperone StpA enables fast RNA refolding by destabilization of mutually exclusive base pairs within competing secondary structure elements

In bacteria RNA gene regulatory elements refold dependent on environmental clues between two or more long-lived conformational states each associated with a distinct regulatory state. The refolding kinetics are strongly temperature-dependent and especially at lower temperatures they reach timescales that are biologically not accessible. To overcome this problem, RNA chaperones have evolved. However, the precise molecular mechanism of how these proteins accelerate RNA refolding reactions remains enigmatic. Here we show how the RNA chaperone StpA of Escherichia coli leads to an acceleration of a bistable RNA's refolding kinetics through the selective destabilization of key base pairing interactions. We find in laser assisted real-time NMR experiments on photocaged bistable RNAs that the RNA chaperone leads to a two-fold increase in refolding rates at low temperatures due to reduced stability of ground state conformations. Further, we can show that upon interaction with StpA, base pairing interactions in the bistable RNA are modulated to favor refolding through the dominant pseudoknotted transition pathway. Our results shed light on the molecular mechanism of the interaction between RNA chaperones and bistable RNAs and are the first step into a functional classification of chaperones dependent on their biophysical mode of operation.

Folding dynamics of polymorphic G-quadruplex structures

G-quadruplexes (G4), found in numerous places within the human genome, are involved in essential processes of cell regulation. Chromosomal DNA G4s are involved for example, in replication and transcription as first steps of gene expression. Hence, they influence a plethora of downstream processes. G4s possess an intricate structure that differs from canonical B-form DNA. Identical DNA G4 sequences can adopt multiple long-lived conformations, a phenomenon known as G4 polymorphism. A detailed understanding of the molecular mechanisms that drive G4 folding is essential to understand their ambivalent regulatory roles. Disentangling the inherent dynamic and polymorphic nature of G4 structures thus is key to unravel their biological functions and make them amenable as molecular targets in novel therapeutic approaches. We here review recent experimental approaches to monitor G4 folding and discuss structural aspects for possible folding pathways. Substantial progress in the understanding of G4 folding within the recent years now allows drawing comprehensive models of the complex folding energy landscape of G4s that we herein evaluate based on computational and experimental evidence.

Real-time nuclear magnetic resonance spectroscopy in the study of biomolecular kinetics and dynamics

The review describes the application of nuclear magnetic resonance (NMR) spectroscopy to study kinetics of folding, refolding and aggregation of proteins, RNA and DNA. Time-resolved NMR experiments can be conducted in a reversible or an irreversible manner. In particular, irreversible folding experiments pose large requirements for (i) signal-to-noise due to the time limitations and (ii) synchronising of the refolding steps. Thus, this contribution discusses the application of methods for signal-to-noise increases, including dynamic nuclear polarisation, hyperpolarisation and photo-CIDNP for the study of time-resolved NMR studies. Further, methods are reviewed ranging from pressure and temperature jump, light induction to rapid mixing to induce rapidly non-equilibrium conditions required to initiate folding.

Spatiotemporal Resolution of Conformational Changes in Biomolecules by Combining Pulsed Electron-Electron Double Resonance Spectroscopy with Microsecond Freeze-Hyperquenching

The function of proteins is linked to their conformations that can be resolved with several high-resolution methods. However, only a few methods can provide the temporal order of intermediates and conformational changes, with each having its limitations. Here, we combine pulsed electron-electron double resonance spectroscopy with a microsecond freeze-hyperquenching setup to achieve spatiotemporal resolution in the angstrom range and lower microsecond time scale. We show that the conformational change of the C_α-helix in the cyclic nucleotide-binding domain of the Mesorhizobium loti potassium channel occurs within about 150 μs and can be resolved with angstrom precision. Thus, this approach holds great promise for obtaining 4D landscapes of conformational changes in biomolecules.

Review and prospect: NMR spectroscopy denoising and reconstruction with low-rank Hankel matrices and tensors

Nuclear magnetic resonance (NMR) spectroscopy is an important analytical tool in chemistry, biology, and life science, but it suffers from relatively low sensitivity and long acquisition time. Thus, improving the apparent signal-to-noise ratio and accelerating data acquisition became indispensable. In this review, we summarize the recent progress on low-rank Hankel matrix and tensor methods, which exploit the exponential property of free-induction decay signals, to enable effective denoising and spectra reconstruction. We also outline future developments that are likely to make NMR spectroscopy a far more powerful technique.

Refolding through a Linear Transition State Enables Fast Temperature Adaptation of a Translational Riboswitch

The adenine-sensing riboswitch from the Gram-negative bacterium Vibrio vulnificus is an RNA-based gene regulatory element that acts in response to both its cognate low-molecular weight ligand and temperature. The combined sensitivity to environmental temperature and ligand concentration is maintained by an equilibrium of three distinct conformations involving two ligand-free states and one ligand-bound state. The key structural element that undergoes refolding in the ligand-free states comprises a 35-nucleotide temperature response module. Here, we present the structural characterization of this temperature response module. We employ high-resolution NMR spectroscopy and photocaged RNAs as molecular probes to decipher the kinetic and thermodynamic framework of the secondary structure transition in the apo state of the riboswitch. We propose a model for the transition state adopted during the thermal refolding of the temperature response module that connects two mutually exclusive long-lived and stable conformational states. This transition state is characterized by a comparatively low free activation enthalpy. A pseudoknot conformation in the transition state, as commonly seen in RNA refolding, is therefore unlikely. More likely, the transition state of the adenine-sensing riboswitch temperature response module features a linear conformation.

Site-Specific Spin Labeling of RNA for NMR and EPR Structural Studies

Many RNA architectures were discovered to be involved in essential biological pathways acting as catalysts and/or regulators of gene expression, transcription, translation, splicing, or viral infection. The key to understand their diverse biological functions is to investigate their structure and dynamic. Nuclear Magnetic Resonance (NMR) is a powerful method to gain insight into these properties. However, the study of high-molecular-weight RNAs by NMR remains challenging. Advances in biochemical and NMR methods over the recent years allow to overcome the limitation of NMR. In particular, the incorporation of paramagnetic probes, coupled to the measurement of the induced effects on nuclear spins, has become an efficient tool providing long-range distance restraints and information on dynamic in solution. At the same time, the use of spin label enabled the application of Electron Paramagnetic Resonance (EPR) to study biological macromolecules. Combining NMR and EPR is emerging as a new approach to investigate the architecture of biological systems.Here, we describe an efficient protocol to introduce a paramagnetic probe into a RNA at a specific position. This method enables various combinations of isotopic labeling for NMR and is also of interest for EPR studies.

A 300-fold enhancement of imino nucleic acid resonances by hyperpolarized water provides a new window for probing RNA refolding by 1D and 2D NMR

NMR sensitivity-enhancement methods involving hyperpolarized water could be of importance for solution-state biophysical investigations. Hyperpolarized water (HyperW) can enhance the ¹H NMR signals of exchangeable sites by orders of magnitude over their thermal counterparts, while providing insight into chemical exchange and solvent accessibility at a site-resolved level. As HyperW's enhancements are achieved by exploiting fast solvent exchanges associated with minimal interscan delays, possibilities for the rapid monitoring of chemical reactions and biomolecular (re)folding are opened. HyperW NMR can also accommodate heteronuclear transfers, facilitating the rapid acquisition of 2-dimensional (2D) ¹⁵N-¹H NMR correlations, and thereby combining an enhanced spectral resolution with speed and sensitivity. This work demonstrates how these qualities can come together for the study of nucleic acids. HyperW injections were used to target the guanine-sensing riboswitch aptamer domain (GSR^apt) of the xpt-pbuX operon in Bacillus subtilis Unlike what had been observed in proteins, where residues benefited of HyperW NMR only if/when sufficiently exposed to water, these enhancements applied to every imino resonance throughout the RNA. The >300-fold enhancements observed in the resulting ¹H NMR spectra allowed us to monitor in real time the changes that GSR^apt undergoes upon binding hypoxanthine, a high-affinity interaction leading to conformational refolding on a ∼1-s timescale at 36 °C. Structural responses could be identified for several nucleotides by 1-dimensional (1D) imino ¹H NMR as well as by 2D HyperW NMR spectra acquired upon simultaneous injection of hyperpolarized water and hypoxanthine. The folding landscape revealed by this HyperW strategy for GSR^apt, is briefly discussed.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

Helix-Based RNA Landscape Partition and Alternative Secondary Structure Determination

RNA is a versatile macromolecule with the ability to fold into and interconvert between multiple functional conformations. The elucidation of the RNA folding landscape, especially the knowledge of alternative structures, is critical to uncover the physical mechanism of RNA functions. Here, we introduce a helix-based strategy for RNA folding landscape partition and alternative secondary structure determination. The benchmark test of 27 RNAs involving alternative stable structures shows that the model has the ability to divide the whole landscape into distinct partitions at the secondary structure level and predict the representative structures for each partition. Furthermore, the predicted structures and equilibrium populations of metastable conformations for the 2'dG-sensing riboswitch reveal the allosteric conformational switch on transcript length, which is consistent with the experimental study, indicating the importance of metastable states for RNA-based gene regulation. The model delivers a starting point for the landscape-based strategy toward the RNA folding mechanism and functions.

Kinetic Mechanism of RNA Helix-Terminal Basepairing-A Kinetic Minima Network Analysis

RNA functions are often kinetically controlled. The folding kinetics of RNAs involves global structural changes and local nucleotide movement, such as base flipping. The most elementary step in RNA folding is the closing and opening of a basepair. By integrating molecular dynamics simulation, master equation, and kinetic Monte Carlo simulation, we investigate the kinetics mechanism of RNA helix-terminal basepairing. The study reveals a six-state folding scheme with three dominant folding pathways of tens, hundreds, and thousands of nanoseconds of folding timescales, respectively. The overall kinetics is rate limited by the detrapping of a misfolded state with the overall folding time of 10^-5 s. Moreover, the analysis examines the different roles of the various driving forces, such as the basepairing and stacking interactions and the ion binding/dissociation effects on structural changes. The results may provide useful insights for developing a basepair opening/closing rate model and further kinetics models of large RNAs.

NMR characterization of RNA small molecule interactions

Exciting discoveries of naturally occurring ligand-sensing and disease-linked noncoding RNAs have promoted significant interests in understanding RNA-small molecule interactions. NMR spectroscopy is a powerful tool for characterizing intermolecular interactions. In this review, we describe protocols and approaches for applying NMR spectroscopy to investigate interactions between RNA and small molecules. We review protocols for RNA sample preparation, methods for identifying RNA-binding small molecules, approaches for mapping RNA-small molecule interactions, determining complex structures, and characterizing binding kinetics. We hope this review will provide a guideline to streamline NMR applications in studying RNA-small molecule interactions, facilitating both basic mechanistic understandings of RNA functions and translational efforts in developing RNA-targeted therapeutics.

Chemical Exchange

The phenomenon of chemical or conformational exchange in NMR spectroscopy has enabled detailed characterization of time-dependent aspects of biomolecular function, including folding, molecular recognition, allostery, and catalysis, on timescales from microsecond to second. Importantly, NMR methods based on a variety of spin relaxation parameters have been developed that provide quantitative information on interconversion kinetics, thermodynamic properties, and structural features of molecular states populated to a fraction of a percent at equilibrium and otherwise unobservable by other NMR approaches. The ongoing development of more sophisticated experimental techniques and the necessity to apply these methods to larger and more complex molecular systems engenders a corresponding need for theoretical advances describing such techniques and facilitating data analysis in applications. This review surveys current aspects of the theory of chemical exchange, as utilized in ZZ-exchange; Hahn and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo; and R_1ρ, chemical exchange saturation transfer (CEST), and dark state saturation transfer (DEST) spin-locking experiments. The review emphasizes theoretical results for kinetic topologies with more than two interconverting states, both to obtain compact analytical forms suitable for data analysis and to establish conditions for distinguishability between alternative kinetic schemes.

Identification and Rationalization of Kinetic Folding Intermediates for a Low-Density Lipoprotein Receptor Ligand-Binding Module

Many mutations that cause familial hypercholesterolemia localize to ligand-binding domain 5 (LA5) of the low-density lipoprotein receptor, motivating investigation of the folding and misfolding of this small, disulfide-rich, calcium-binding domain. LA5 folding is known to involve non-native disulfide isomers, yet these folding intermediates have not been structurally characterized. To provide insight into these intermediates, we used nuclear magnetic resonance (NMR) to follow LA5 folding in real time. We demonstrate that misfolded or partially folded disulfide intermediates are indistinguishable from the unfolded state when focusing on the backbone NMR signals, which provide information on the formation of only the final, native state. However, ¹³C labeling of cysteine side chains differentiated transient intermediates from the unfolded and native states and reported on disulfide bond formation in real time. The cysteine pairings in a dominant intermediate were identified using ¹³C-edited three-dimensional NMR, and coarse-grained molecular dynamics simulations were used to investigate the preference of this disulfide set over other non-native arrangements. The transient population of LA5 species with particular non-native cysteine connectitivies during folding supports the conclusion that cysteine pairing is not random and that there is a bias toward certain structural ensembles during the folding process, even prior to the binding of calcium.

Efficient Detection of Structure and Dynamics in Unlabeled RNAs: The SELOPE Approach

The knowledge of structure and dynamics is crucial to explain the function of RNAs. While nuclear magnetic resonance (NMR) is well suited to probe these for complex biomolecules, it requires expensive, isotopically labeled samples, and long measurement times. Here we present SELOPE, a new robust, proton-only NMR method that allows us to obtain site-specific overview of structure and dynamics in an entire RNA molecule using an unlabeled sample. SELOPE simplifies assignment and allows for cost-effective screening of the response of nucleic acids to physiological changes (e.g. ion concentration) or screening of drugs in a high throughput fashion. This single technique allows us to probe an unprecedented range of exchange time scales (the whole μs to ms motion range) with increased sensitivity, surpassing all current experiments to detect chemical exchange. For the first time we could describe an RNA excited state using an unlabeled RNA.

Life times of metastable states guide regulatory signaling in transcriptional riboswitches

Transcriptional riboswitches modulate downstream gene expression by a tight coupling of ligand-dependent RNA folding kinetics with the rate of transcription. RNA folding pathways leading to functional ON and OFF regulation involve the formation of metastable states within well-defined sequence intervals during transcription. The kinetic requirements for the formation and preservation of these metastable states in the context of transcription remain unresolved. Here, we reversibly trap the previously defined regulatory relevant metastable intermediate of the Mesoplasma florum 2′-deoxyguanosine (2′dG)-sensing riboswitch using a photocaging-ligation approach, and monitor folding to its native state by real-time NMR in both presence and absence of ligand. We further determine transcription rates for two different bacterial RNA polymerases. Our results reveal that the riboswitch functions only at transcription rates typical for bacterial polymerases (10–50 nt s⁻¹) and that gene expression is modulated by 40–50% only, while subtle differences in folding rates guide population ratios within the structural ensemble to a specific regulatory outcome.

Riboswitches are RNA-based regulatory elements, which regulate downstream gene expression by binding of small molecular weight ligands. Here the authors demonstrate the molecular mechanism of a transcriptional riboswitch that integrates changes in transcription rates, metabolite concentration, and kinetic on- and off-rates of ligand binding.

Rotation of Guanine Amino Groups in G-Quadruplexes: A Probe for Local Structure and Ligand Binding

Nucleic acids are dynamic molecules whose functions may depend on their conformational fluctuations and local motions. In particular, amino groups are dynamic components of nucleic acids that participate in the formation of various secondary structures such as G-quadruplexes. Here, we present a cost-efficient NMR method to quantify the rotational dynamics of guanine amino groups in G-quadruplex nucleic acids. An isolated spectrum of amino protons from a specific tetrad-bound guanine can be extracted from the nuclear Overhauser effect spectroscopy spectrum based on the close proximity between the intra-residue imino and amino protons. We apply the method in different structural contexts of G-quadruplexes and their complexes. Our results highlight the role of stacking and hydrogen-bond interactions in restraining amino-group rotation. The measurement of the rotation rate of individual amino groups could give insight into the dynamic processes occurring at specific locations within G-quadruplex nucleic acids, providing valuable probes for local structure, dynamics, and ligand binding.

Resolving sugar puckers in RNA excited states exposes slow modes of repuckering dynamics

Recent studies have shown that RNAs exist in dynamic equilibrium with short-lived low-abundance 'excited states' that form by reshuffling base pairs in and around non-canonical motifs. These conformational states are proposed to be rich in non-canonical motifs and to play roles in the folding and regulatory functions of non-coding RNAs but their structure proves difficult to characterize given their transient nature. Here, we describe an approach for determining sugar pucker conformation in RNA excited states through nuclear magnetic resonance measurements of C1΄ and C4΄ rotating frame spin relaxation (R1ρ) in uniformly 13C/15N labeled RNA samples. Application to HIV-1 TAR exposed slow modes of sugar repuckering dynamics at the μs and ms timescale accompanying transitions between non-helical (C2΄-endo) to helical (C3΄-endo) conformations during formation of two distinct excited states. In contrast, we did not obtain any evidence for slow sugar repuckering dynamics for nucleotides in a variety of structural contexts that do not undergo non-helical to helical transitions. Our results outline a route for significantly improving the conformational characterization of RNA excited states and suggest that slow modes of repuckering dynamics gated by transient changes in secondary structure are quite common in RNA.

Optimized fast mixing device for real-time NMR applications

We present an improved fast mixing device based on the rapid mixing of two solutions inside the NMR probe, as originally proposed by Hore and coworkers (J. Am. Chem. Soc. 125 (2003) 12484-12492). Such a device is important for off-equilibrium studies of molecular kinetics by multidimensional real-time NMR spectrsocopy. The novelty of this device is that it allows removing the injector from the NMR detection volume after mixing, and thus provides good magnetic field homogeneity independently of the initial sample volume placed in the NMR probe. The apparatus is simple to build, inexpensive, and can be used without any hardware modification on any type of liquid-state NMR spectrometer. We demonstrate the performance of our fast mixing device in terms of improved magnetic field homogeneity, and show an application to the study of protein folding and the structural characterization of transiently populated folding intermediates.

Studying Dynamics by Magic-Angle Spinning Solid-State NMR Spectroscopy: Principles and Applications to Biomolecules

Magic-angle spinning solid-state NMR spectroscopy is an important technique to study molecular structure, dynamics and interactions, and is rapidly gaining importance in biomolecular sciences. Here we provide an overview of experimental approaches to study molecular dynamics by MAS solid-state NMR, with an emphasis on the underlying theoretical concepts and differences of MAS solid-state NMR compared to solution-state NMR. The theoretical foundations of nuclear spin relaxation are revisited, focusing on the particularities of spin relaxation in solid samples under magic-angle spinning. We discuss the range of validity of Redfield theory, as well as the inherent multi-exponential behavior of relaxation in solids. Experimental challenges for measuring relaxation parameters in MAS solid-state NMR and a few recently proposed relaxation approaches are discussed, which provide information about time scales and amplitudes of motions ranging from picoseconds to milliseconds. We also discuss the theoretical basis and experimental measurements of anisotropic interactions (chemical-shift anisotropies, dipolar and quadrupolar couplings), which give direct information about the amplitude of motions. The potential of combining relaxation data with such measurements of dynamically-averaged anisotropic interactions is discussed. Although the focus of this review is on the theoretical foundations of dynamics studies rather than their application, we close by discussing a small number of recent dynamics studies, where the dynamic properties of proteins in crystals are compared to those in solution.

Applications of high dimensionality experiments to biomolecular NMR

High dimensionality NMR experiments facilitate resonance assignment and precise determination of spectral parameters such as coupling constants. Sparse non-uniform sampling enables acquisition of experiments of high dimensionality with high resolution in acceptable time. In this review we present and compare some significant applications of NMR experiments of dimensionality higher than three in the field of biomolecular studies in solution.

Rapid NMR screening of RNA secondary structure and binding

Determination of RNA secondary structures by NMR spectroscopy is a useful tool e.g. to elucidate RNA folding space or functional aspects of regulatory RNA elements. However, current approaches of RNA synthesis and preparation are usually time-consuming and do not provide analysis with single nucleotide precision when applied for a large number of different RNA sequences. Here, we significantly improve the yield and 3' end homogeneity of RNA preparation by in vitro transcription. Further, by establishing a native purification procedure with increased throughput, we provide a shortcut to study several RNA constructs simultaneously. We show that this approach yields μmol quantities of RNA with purities comparable to PAGE purification, while avoiding denaturation of the RNA.

Unveiling Inherent Degeneracies in Determining Population-Weighted Ensembles of Interdomain Orientational Distributions Using NMR Residual Dipolar Couplings: Application to RNA Helix Junction Helix Motifs

A growing number of studies employ time-averaged experimental data to determine dynamic ensembles of biomolecules. While it is well-known that different ensembles can satisfy experimental data to within error, the extent and nature of these degeneracies, and their impact on the accuracy of the ensemble determination remains poorly understood. Here, we use simulations and a recently introduced metric for assessing ensemble similarity to explore degeneracies in determining ensembles using NMR residual dipolar couplings (RDCs) with specific application to A-form helices in RNA. Various target ensembles were constructed representing different domain-domain orientational distributions that are confined to a topologically restricted (<10%) conformational space. Five independent sets of ensemble averaged RDCs were then computed for each target ensemble and a "sample and select" scheme used to identify degenerate ensembles that satisfy RDCs to within experimental uncertainty. We find that ensembles with different ensemble sizes and that can differ significantly from the target ensemble (by as much as ∑Ω ∼ 0.4 where ∑Ω varies between 0 and 1 for maximum and minimum ensemble similarity, respectively) can satisfy the ensemble averaged RDCs. These deviations increase with the number of unique conformers and breadth of the target distribution, and result in significant uncertainty in determining conformational entropy (as large as 5 kcal/mol at T = 298 K). Nevertheless, the RDC-degenerate ensembles are biased toward populated regions of the target ensemble, and capture other essential features of the distribution, including the shape. Our results identify ensemble size as a major source of uncertainty in determining ensembles and suggest that NMR interactions such as RDCs and spin relaxation, on their own, do not carry the necessary information needed to determine conformational entropy at a useful level of precision. The framework introduced here provides a general approach for exploring degeneracies in ensemble determination for different types of experimental data.

Structural fidelity and NMR relaxation analysis in a prototype RNA hairpin

RNA hairpins are widespread and very stable motifs that contribute decisively to RNA folding and biological function. The GTP1G2C3A4C5U6U7C8G9G10U11G12C13C14 construct (with a central UUCG tetraloop) has been extensively studied by solution NMR, and offers and excellent opportunity to evaluate the structure and dynamical description afforded by molecular dynamics (MD) simulations. Here, we compare average structural parameters and NMR relaxation rates estimated from a series of multiple independent explicit solvent MD simulations using the two most recent RNA AMBER force fields (ff99 and ff10). Predicted overall tumbling times are ∼20% faster than those inferred from analysis of NMR data and follow the same trend when temperature and ionic strength is varied. The Watson-Crick stem and the "canonical" UUCG loop structure are maintained in most simulations including the characteristic syn conformation along the glycosidic bond of G9, although some key hydrogen bonds in the loop are partially disrupted. Our analysis pinpoints G9-G10 backbone conformations as a locus of discrepancies between experiment and simulation. In general the results for the more recent force-field parameters (ff10) are closer to experiment than those for the older ones (ff99). This work provides a comprehensive and detailed comparison of state of the art MD simulations against a wide variety of solution NMR measurements.

RNA folding: structure prediction, folding kinetics and ion electrostatics

Beyond the "traditional" functions such as gene storage, transport and protein synthesis, recent discoveries reveal that RNAs have important "new" biological functions including the RNA silence and gene regulation of riboswitch. Such functions of noncoding RNAs are strongly coupled to the RNA structures and proper structure change, which naturally leads to the RNA folding problem including structure prediction and folding kinetics. Due to the polyanionic nature of RNAs, RNA folding structure, stability and kinetics are strongly coupled to the ion condition of solution. The main focus of this chapter is to review the recent progress in the three major aspects in RNA folding problem: structure prediction, folding kinetics and ion electrostatics. This chapter will introduce both the recent experimental and theoretical progress, while emphasize the theoretical modelling on the three aspects in RNA folding.

Characterization of conformational dynamics of bistable RNA by equilibrium and non-equilibrium NMR

Unlike proteins, a given RNA sequence can adopt more than a single conformation. The two (or more) conformations are long-lived and have similar stabilities, but interconvert only slowly. Such bi- or multistability is often linked to the biological functions of the RNA. This unit describes how nuclear magnetic resonance (NMR) spectroscopy can be used to characterize the conformational dynamics of bistable RNAs.

Real-time protein NMR spectroscopy and investigation of assisted protein folding

BACKGROUND

During protein-folding reactions toward the native structure, short-lived intermediate states can be populated. Such intermediates expose hydrophobic patches and can self-associate leading to non-productive protein misfolding. A major focus of current research is the characterization of short-lived intermediates and how molecular chaperones enable productive folding. Real-time NMR spectroscopy, together with the development of advanced methods, is reviewed here and the potential these methods have to characterize intermediate states as well as interactions with molecular chaperone proteins at single-residue resolution is highlighted.

SCOPE OF REVIEW

Various chaperone interactions can guide the protein-folding reaction and thus are important for protein structure formation, stability, and activity of their substrates. Chaperone-assisted protein folding, characterization of intermediates, and their molecular interactions using real-time NMR spectroscopy will be discussed. Additionally, recent advances in NMR methods employed for characterization of high-energy intermediates will be discussed.

MAJOR CONCLUSIONS

Real-time NMR combines high resolution with kinetic information of protein reactions, which can be employed not only for protein-folding studies and the characterization of folding intermediates but also to investigate the molecular mechanisms of assisted protein folding.

GENERAL SIGNIFICANCE

Real-time NMR spectroscopy remains an effective tool to reveal structural details about the interaction between chaperones and transient intermediates. Methodologically, it provides in-depth understanding of how kinetic intermediates and their thermodynamics contribute to the protein-folding reaction. This review summarizes the most recent advances in this field. This article is part of a Special Issue titled Proline-directed Foldases: Cell Signaling Catalysts and Drug Targets.

Memory-efficient RNA energy landscape exploration

MOTIVATION

Energy landscapes provide a valuable means for studying the folding dynamics of short RNA molecules in detail by modeling all possible structures and their transitions. Higher abstraction levels based on a macro-state decomposition of the landscape enable the study of larger systems; however, they are still restricted by huge memory requirements of exact approaches.

RESULTS

We present a highly parallelizable local enumeration scheme that enables the computation of exact macro-state transition models with highly reduced memory requirements. The approach is evaluated on RNA secondary structure landscapes using a gradient basin definition for macro-states. Furthermore, we demonstrate the need for exact transition models by comparing two barrier-based approaches, and perform a detailed investigation of gradient basins in RNA energy landscapes.

AVAILABILITY AND IMPLEMENTATION

Source code is part of the C++ Energy Landscape Library available at http://www.bioinf.uni-freiburg.de/Software/.

A fully enzymatic method for site-directed spin labeling of long RNA

Site-directed spin labeling is emerging as an essential tool to investigate the structural and dynamical features of RNA. We propose here an enzymatic method, which allows the insertion of a paramagnetic center at a specific position in an RNA molecule. The technique is based on a segmental approach using a ligation protocol with T4 RNA ligase 2. One transcribed acceptor RNA is ligated to a donor RNA in which a thio-modified nucleotide is introduced at its 5′-end by in vitro transcription with T7 RNA polymerase. The paramagnetic thiol-specific reagent is subsequently attached to the RNA ligation product. This novel strategy is demonstrated by introducing a paramagnetic probe into the 55 nucleotides long RNA corresponding to K-turn and Specifier Loop domains from the Bacillus subtilis tyrS T-Box leader RNA. The efficiency of the coupling reaction and the quality of the resulting spin-labeled RNA were assessed by Mass Spectrometry, Electron Paramagnetic Resonance (EPR) and Nuclear Magnetic Resonance (NMR). This method enables various combinations of isotopic segmental labeling and spin labeling schemes, a strategy that will be of particular interest to investigate the structural and dynamical properties of large RNA complexes by NMR and EPR spectroscopies.

The importance of helix P1 stability for structural pre-organization and ligand binding affinity of the adenine riboswitch aptamer domain

We report here an in-depth characterization of the aptamer domain of the transcriptional adenine-sensing riboswitch (pbuE) by NMR and fluorescence spectroscopy. By NMR studies, the structure of two aptamer sequences with different lengths of the helix P1, the central element involved in riboswitch conformational switching, was characterized. Hydrogen-bond interactions could be mapped at nucleotide resolution providing information about secondary and tertiary structure, structure homogeneity and dynamics. Our study reveals that the elongation of helix P1 has pronounced effects not only on the local but on the global structure of the apo aptamer domain. The structural differences induced by stabilizing helix P1 were found to be linked to changes of the ligand binding affinity as revealed from analysis of kinetic and thermodynamic data obtained from stopped-flow fluorescence studies. The results provide new insight into the sequence-dependent fine tuning of the structure and function of purine-sensing riboswitches.

Time-resolved multidimensional NMR with non-uniform sampling

Time-resolved experiments demand high resolution both in spectral dimensions and in time of the studied kinetic process. The latter requirement traditionally prohibits applications of the multidimensional experiments, which, although capable of providing invaluable information about structure and dynamics and almost unlimited spectral resolution, require too lengthy data collection. Our work shows that the problem has a solution in using modern methods of NMR data collection and signal processing. A continuous fast pulsing three-dimensional experiment is acquired using non-uniform sampling during full time of the studied reaction. High sensitivity and time-resolution of a few minutes is achieved by simultaneous processing of the full data set with the multi-dimensional decomposition. The method is verified and illustrated in realistic simulations and by measuring deuterium exchange rates of amide protons in ubiquitin. We applied the method for characterizing kinetics of in vitro phosphorylation of two tyrosine residues in an intrinsically disordered cytosolic domain of the B cell receptor protein CD79b. Signals of many residues including tyrosines in both phosphorylated and unmodified forms of CD79b are found in a heavily crowded region of 2D ¹H-¹⁵N correlation spectrum and the significantly enhanced spectral resolution provided by the 3D time-resolved approach was essential for the quantitative site-specific analysis.

A method of determining RNA conformational ensembles using structure-based calculations of residual dipolar couplings

We describe a method of determining the conformational fluctuations of RNA based on the incorporation of nuclear magnetic resonance (NMR) residual dipolar couplings (RDCs) as replica-averaged structural restraints in molecular dynamics simulations. In this approach, the alignment tensor required to calculate the RDCs corresponding to a given conformation is estimated from its shape, and multiple replicas of the RNA molecule are simulated simultaneously to reproduce in silico the ensemble-averaging procedure performed in the NMR measurements. We provide initial evidence that with this approach it is possible to determine accurately structural ensembles representing the conformational fluctuations of RNA by applying the reference ensemble test to the trans-activation response element of the human immunodeficiency virus type 1.

SMMRNA: a database of small molecule modulators of RNA

We have developed SMMRNA, an interactive database, available at http://www.smmrna.org, with special focus on small molecule ligands targeting RNA. Currently, SMMRNA consists of ∼770 unique ligands along with structural images of RNA molecules. Each ligand in the SMMRNA contains information such as Kd, Ki, IC50, ΔTm, molecular weight (MW), hydrogen donor and acceptor count, XlogP, number of rotatable bonds, number of aromatic rings and 2D and 3D structures. These parameters can be explored using text search, advanced search, substructure and similarity-based analysis tools that are embedded in SMMRNA. A structure editor is provided for 3D visualization of ligands. Advance analysis can be performed using substructure and OpenBabel-based chemical similarity fingerprints. Upload facility for both RNA and ligands is also provided. The physicochemical properties of the ligands were further examined using OpenBabel descriptors, hierarchical clustering, binning partition and multidimensional scaling. We have also generated a 3D conformation database of ligands to support the structure and ligand-based screening. SMMRNA provides comprehensive resource for further design, development and refinement of small molecule modulators for selective targeting of RNA molecules.

Characterizing slow chemical exchange in nucleic acids by carbon CEST and low spin-lock field R(1ρ) NMR spectroscopy

Quantitative characterization of dynamic exchange between various conformational states provides essential insights into the molecular basis of many regulatory RNA functions. Here, we present an application of nucleic-acid-optimized carbon chemical exchange saturation transfer (CEST) and low spin-lock field R(1ρ) relaxation dispersion (RD) NMR experiments in characterizing slow chemical exchange in nucleic acids that is otherwise difficult if not impossible to be quantified by the ZZ-exchange NMR experiment. We demonstrated the application on a 47-nucleotide fluoride riboswitch in the ligand-free state, for which CEST and R(1ρ) RD profiles of base and sugar carbons revealed slow exchange dynamics involving a sparsely populated (p ~ 10%) and shortly lived (τ ~ 10 ms) NMR "invisible" state. The utility of CEST and low spin-lock field R(1ρ) RD experiments in studying slow exchange was further validated in characterizing an exchange as slow as ~60 s(-1).

Probing high-affinity 11-mer DNA aptamer against Lup an 1 (β-conglutin)

Aptamers are synthetic nucleic acids with great potential as analytical tools. However, the length of selected aptamers (typically 60-100 bases) can affect affinity, due to the presence of bases not required for interaction with the target, and therefore, the truncation of these selected sequences and identification of binding domains is a critical step to produce potent aptamers with higher affinities and specificities and lowered production costs. In this paper we report the truncation of an aptamer that specifically binds to β-conglutin (Lup an 1), an anaphylactic allergen. Through comparing the predicted secondary structures of the aptamers, a hairpin structure with a G-rich loop was determined to be the binding motif. The highest affinity was observed with a truncation resulting in an 11-mer sequence that had an apparent equilibrium dissociation constant (K D) of 1.7 × 10(-9) M. This 11-mer sequence was demonstrated to have high specificity for β-conglutin and showed no cross-reactivity to other lupin conglutins (α-, δ-, γ-conglutins) and closely related proteins such as gliadin. Finally, the structure of the truncated 11-mer aptamer was preliminarily elucidated, and the GQRS Mapper strongly predicted the presence of a G-quadruplex, which was subsequently corroborated using one-dimensional NMR, thus highlighting the stability of the truncated structure.

Photo-CIDNP NMR spectroscopy of amino acids and proteins

Photo-chemically induced dynamic nuclear polarization (CIDNP) is a nuclear magnetic resonance (NMR) phenomenon which, among other things, is exploited to extract information on biomolecular structure via probing solvent-accessibilities of tryptophan (Trp), tyrosine (Tyr), and histidine (His) amino acid side chains both in polypeptides and proteins in solution. The effect, normally triggered by a (laser) light-induced photochemical reaction in situ, yields both positive and/or negative signal enhancements in the resulting NMR spectra which reflect the solvent exposure of these residues both in equilibrium and during structural transformations in "real time". As such, the method can offer - qualitatively and, to a certain extent, quantitatively - residue-specific structural and kinetic information on both the native and, in particular, the non-native states of proteins which, often, is not readily available from more routine NMR techniques. In this review, basic experimental procedures of the photo-CIDNP technique as applied to amino acids and proteins are discussed, recent improvements to the method highlighted, and future perspectives presented. First, the basic principles of the phenomenon based on the theory of the radical pair mechanism (RPM) are outlined. Second, a description of standard photo-CIDNP applications is given and it is shown how the effect can be exploited to extract residue-specific structural information on the conformational space sampled by unfolded or partially folded proteins on their "path" to the natively folded form. Last, recent methodological advances in the field are highlighted, modern applications of photo-CIDNP in the context of biological NMR evaluated, and an outlook into future perspectives of the method is given.

Kinetic mechanism of conformational switch between bistable RNA hairpins

Transitions between the different conformational states play a critical role in many RNA catalytic and regulatory functions. In this study, we use the Kinetic Monte Carlo method to investigate the kinetic mechanism for the conformational switches between bistable RNA hairpins. We find three types of conformational switch pathways for RNA hairpins: refolding after complete unfolding, folding through basepair-exchange pathways and through pseudoknot-assisted pathways, respectively. The result of the competition between the three types of pathways depends mainly on the location of the rate-limiting base stacks (such as the GC base stacks) in the structures. Depending on the structural relationships between the two bistable hairpins, the conformational switch can follow single or multiple dominant pathways. The predicted folding pathways are supported by the activation energy results derived from the Arrhenius plot as well as the NMR spectroscopy data.

Light-controlled tools

Spatial and temporal control over chemical and biological processes plays a key role in life, where the whole is often much more than the sum of its parts. Quite trivially, the molecules of a cell do not form a living system if they are only arranged in a random fashion. If we want to understand these relationships and especially the problems arising from malfunction, tools are necessary that allow us to design sophisticated experiments that address these questions. Highly valuable in this respect are external triggers that enable us to precisely determine where, when, and to what extent a process is started or stopped. Light is an ideal external trigger: It is highly selective and if applied correctly also harmless. It can be generated and manipulated with well-established techniques, and many ways exist to apply light to living systems--from cells to higher organisms. This Review will focus on developments over the last six years and includes discussions on the underlying technologies as well as their applications.

RNA plasticity and selectivity applicable to therapeutics and novel biosensor development

Aptamers are short, single-stranded nucleic acid sequences that are selected in vitro from large oligonucleotide libraries based on their high affinity to a target molecule. Hence, aptamers can be thought of as a nucleic acid analog to antibodies. However, several viewpoints hold that the potential of aptamers arises from interesting characteristics that are distinct from, or in some cases, superior to those of antibodies. This review summarizes the recent achievements in aptamer programs developed in our laboratory against basic and therapeutic protein targets. Through these studies, we became aware of the remarkable conformational plasticity and selectivity of RNA, on which the published report has not shed much light even though this is evidently a crucial feature for the strong specificity and affinity of RNA aptamers.

Mapping the landscape of RNA dynamics with NMR spectroscopy

Among the three major classes of biomacromolecules (DNA, RNA, and proteins) RNA's pronounced dynamics are the most explicitly linked to its wide variety of functions, which include catalysis and the regulation of transcription, translation, and splicing. These functions are mediated by a range of RNA biomachinery, including such varied examples as macromolecular noncoding RNAs, microRNAs, small interfering RNAs, riboswitch RNAs, and RNA thermometers. In each case, the functional dynamics of an interconversion is characterized by an associated rate constant. In this Account, we provide an introduction to NMR spectroscopic characterization of the landscape of RNA dynamics. We introduce strategies for measuring NMR parameters at various time scales as well as the underlying models for describing the corresponding rate constants. RNA exhibits significant dynamic motion, which can be modulated by (i) intermolecular interactions, including specific and nonspecific binding of ions (such as Mg(2+) and tertiary amines), (ii) metabolites in riboswitches or RNA aptamers, and (iii) macromolecular interactions within ribonucleic protein particles, including the ribosome and the spliceosome. Our understanding of the nature of these dynamic changes in RNA targets is now being incorporated into RNA-specific approaches in the design of RNA inhibitors. Interactions of RNA with proteins, other RNAs, or small molecules often occur through binding mechanisms that follow an induced fit mechanism or a conformational selection mechanism, in which one of several populated RNA conformations is selected through ligand binding. The extent of functional dynamics, including the kinetic formation of a specific RNA tertiary fold, is dependent on the messenger RNA (mRNA) chain length. Thus, during de novo synthesis of mRNA, both in prokaryotes and eukaryotes, nascent mRNA of various lengths will adopt different secondary and tertiary structures. The speed of transcription has a critical influence on the functional dynamics of the RNA being synthesized. In addition to modulating the local dynamics of a conformational RNA ensemble, a given RNA sequence may adopt more than one global, three-dimensional structure. RNA modification is one way to select among these alternative structures, which are often characterized by nearly equal stability, but with high energy barriers for conformational interconversion. The refolding of different secondary and tertiary structures has been found to be a major regulatory mechanism for transcription and translation. These conformational transitions can be characterized with NMR spectroscopy, for any given RNA sequence, in response to external stimuli.

Characterizing RNA dynamics at atomic resolution using solution-state NMR spectroscopy

Many recently discovered noncoding RNAs do not fold into a single native conformation but sample many different conformations along their free-energy landscape to carry out their biological function. Here we review solution-state NMR techniques that measure the structural, kinetic and thermodynamic characteristics of RNA motions spanning picosecond to second timescales at atomic resolution, allowing unprecedented insights into the RNA dynamic structure landscape. From these studies a basic description of the RNA dynamic structure landscape is emerging, bringing new insights into how RNA structures change to carry out their function as well as applications in RNA-targeted drug discovery and RNA bioengineering.

Transient RNA-protein interactions in RNA folding

The RNA folding trajectory features numerous off-pathway folding traps, which represent conformations that are often equally as stable as the native functional ones. Therefore, the conversion between these off-pathway structures and the native correctly folded ones is the critical step in RNA folding. This process, referred to as RNA refolding, is slow, and is represented by a transition state that has a characteristic high free energy. Because this kinetically limiting process occurs in vivo, proteins (called RNA chaperones) have evolved that facilitate the (re)folding of RNA molecules. Here, we present an overview of how proteins interact with RNA molecules in order to achieve properly folded states. In this respect, the discrimination between static and transient interactions is crucial, as different proteins have evolved a multitude of mechanisms for RNA remodeling. For RNA chaperones that act in a sequence-unspecific manner and without the use of external sources of energy, such as ATP, transient RNA–protein interactions represent the basis of the mode of action. By presenting stretches of positively charged amino acids that are positioned in defined spatial configurations, RNA chaperones enable the RNA backbone, via transient electrostatic interactions, to sample a wider conformational space that opens the route for efficient refolding reactions.

Probing RNA dynamics via longitudinal exchange and CPMG relaxation dispersion NMR spectroscopy using a sensitive 13C-methyl label

The refolding kinetics of bistable RNA sequences were studied in unperturbed equilibrium via ¹³C exchange NMR spectroscopy. For this purpose a straightforward labeling technique was elaborated using a 2′-¹³C-methoxy uridine modification, which was prepared by a two-step synthesis and introduced into RNA using standard protocols. Using ¹³C longitudinal exchange NMR spectroscopy the refolding kinetics of a 20 nt bistable RNA were characterized at temperatures between 298 and 310 K, yielding the enthalpy and entropy differences between the conformers at equilibrium and the activation energy of the refolding process. The kinetics of a more stable 32 nt bistable RNA could be analyzed by the same approach at elevated temperatures, i.e. at 314 and 316 K. Finally, the dynamics of a multi-stable RNA able to fold into two hairpin- and a pseudo-knotted conformation was studied by ¹³C relaxation dispersion NMR spectroscopy.