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

Structure and Internal Dynamics of Short RNA Duplexes Determined by a Combination of Pulsed EPR Methods and MD Simulations

We used EPR spectroscopy to characterize the structure of RNA duplexes and their internal twist, stretch and bending motions. We prepared eight 20-base-pair-long RNA duplexes containing the rigid spin-label Çm, a cytidine analogue, at two positions and acquired orientation-selective PELDOR/DEER data. By using different frequency bands (X-, Q-, G-band), detailed information about the distance and orientation of the labels was obtained and provided insights into the global conformational dynamics of the RNA duplex. We used 19F Mims ENDOR experiments on three singly Çm- and singly fluorine-labeled RNA duplexes to determine the exact position of the Çm spin label in the helix. In a quantitative comparison to MD simulations of RNA with and without Çm spin labels, we found that state-of-the-art force fields with explicit parameterization of the spin label were able to describe the conformational ensemble present in our experiments. The MD simulations further confirmed that the Çm spin labels are excellent mimics of cytidine inducing only small local changes in the RNA structure. Çm spin labels are thus ideally suited for high-precision EPR experiments to probe the structure and, in conjunction with MD simulations, motions of RNA.

Cellular Applications of DNP Solid-State NMR - State of the Art and a Look to the Future

Sensitivity enhanced dynamic nuclear polarization solid-state NMR is emerging as a powerful technique for probing the structural properties of conformationally homogenous and heterogenous biomolecular species irrespective of size at atomic resolution within their native environments. Herein we detail advancements that have made acquiring such data, specifically within the confines of intact bacterial and eukaryotic cell a reality and further discuss the type of structural information that can presently be garnered by the technique's exploitation. Subsequently, we discuss bottlenecks that have thus far curbed cellular DNP-ssNMR's broader adoption namely due a lack of sensitivity and spectral resolution. We also explore possible solutions ranging from utilization of new pulse sequences, design of better performing polarizing agents, and application of additional biochemical/ cell biological methodologies.

Probing RNA structure and dynamics using nanopore and next generation sequencing

It has become increasingly evident that the structures RNAs adopt are conformationally dynamic; the various structured states that RNAs sample govern their interactions with other nucleic acids, proteins, and ligands to regulate a myriad of biological processes. Although several biophysical approaches have been developed and used to study the dynamic landscape of structured RNAs, technical limitations have limited their application to all classes of RNA due to variable size and flexibility. Recent advances combining chemical probing experiments with next-generation- and direct sequencing have emerged as an alternative approach to exploring the conformational dynamics of RNA. In this review, we provide a methodological overview of the sequencing-based techniques used to study RNA conformational dynamics. We discuss how different techniques have enabled us to better understand the propensity of RNAs from a variety of different classes to sample multiple conformational states. Finally, we present examples of the ways these techniques have reshaped how we think about RNA structure.

RNA-Small-Molecule Interaction: Challenging the "Undruggable" Tag

RNA targeting, specifically with small molecules, is a relatively new and rapidly emerging avenue with the promise to expand the target space in the drug discovery field. From being "disregarded" as an "undruggable" messenger molecule to FDA approval of an RNA-targeting small-molecule drug Risdiplam, a radical change in perspective toward RNA has been observed in the past decade. RNAs serve important regulatory functions beyond canonical protein synthesis, and their dysregulation has been reported in many diseases. A deeper understanding of RNA biology reveals that RNA molecules can adopt a variety of structures, carrying defined binding pockets that can accommodate small-molecule drugs. Due to its functional diversity and structural complexity, RNA can be perceived as a prospective target for therapeutic intervention. This perspective highlights the proof of concept of RNA-small-molecule interactions, exemplified by targeting of various transcripts with functional modulators. The advent of RNA-oriented knowledge would help expedite drug discovery.

Visualizing a two-state conformational ensemble in stem-loop 3 of the transcriptional regulator 7SK RNA

Structural plasticity is integral to RNA function; however, there are currently few methods to quantitatively resolve RNAs that have multiple structural states. NMR spectroscopy is a powerful approach for resolving conformational ensembles but is size-limited. Chemical probing is well-suited for large RNAs but provides limited structural and kinetics information. Here, we integrate the two approaches to visualize a two-state conformational ensemble for the central stem-loop 3 (SL3) of 7SK RNA, a critical element for 7SK RNA function in transcription regulation. We find that the SL3 distal end exchanges between two equally populated yet structurally distinct states in both isolated SL3 constructs and full-length 7SK RNA. We rationally designed constructs that lock SL3 into a single state and demonstrate that both chemical probing and NMR data fit to a linear combination of the two states. Comparison of vertebrate 7SK RNA sequences shows either or both states are highly conserved. These results provide new insights into 7SK RNA structural dynamics and demonstrate the utility of integrating chemical probing with NMR spectroscopy to gain quantitative insights into RNA conformational ensembles.

Z-Form Adoption of Nucleic Acid is a Multi-Step Process Which Proceeds through a Melted Intermediate

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of innate immune response proteins. Zα domains stabilize this higher-energy conformation by making specific interactions with the unique geometry of Z-DNA/Z-RNA. However, the mechanism by which a right-handed helix contorts to become left-handed in the presence of proteins, including the intermediate steps involved, is poorly understood. Through a combination of nuclear magnetic resonance (NMR) and other biophysical measurements, we have determined that in the absence of Zα, under low salt conditions at room temperature, d(CpG) and r(CpG) constructs show no observable evidence of transient Z-conformations greater than 0.5% on either the intermediate or slow NMR time scales. At higher temperatures, we observed a transient unfolded intermediate. The ease of melting a nucleic acid duplex correlates with Z-form adoption rates in the presence of Zα. The largest contributing factor to the activation energies of Z-form adoption as calculated by Arrhenius plots is the ease of flipping the sugar pucker, as required for Z-DNA and Z-RNA. Together, these data validate the previously proposed "zipper model" for Z-form adoption in the presence of Zα. Overall, Z-conformations are more likely to be adopted by double-stranded DNA and RNA regions flanked by less stable regions and by RNAs experiencing torsional/mechanical stress.

Dynamic geometry design of cyclic peptide architectures for RNA structure

Designing inhibitors for RNA is still challenging due to the bottleneck of maintaining the binding interaction of inhibitor-RNA accompanied by subtle RNA flexibility. Thus, the current approach usually needs to screen thousands of candidate inhibitors for binding. Here, we propose a dynamic geometry design approach to enrich the hits with only a tiny pool of designed geometrically compatible scaffold candidates. First, our method uses graph-based tree decomposition to explore the complementarity rigid binding cyclic peptide and design the amino acid side chain length and charge to fit the RNA pocket. Then, we perform an energy-based dynamical network algorithm to optimize the inhibitor-RNA hydrogen bonds. Dynamic geometry-guided design yields successful inhibitors with low micromolar binding affinity scaffolds and experimentally competes with the natural RNA chaperone. The results indicate that the dynamic geometry method yields higher efficiency and accuracy than traditional methods. The strategy could be further optimized to design the length and chirality by adopting nonstandard amino acids and facilitating RNA engineering for biological or medical applications.

Chemical Analysis of an Isotopically Labeled Molecule Using Two-Dimensional NMR Spectroscopy at 34 μT

Low-field nuclear magnetic resonance (NMR) spectroscopy, conducted at or below a few millitesla, provides only limited spectral information due to its inability to resolve chemical shifts. Thus, chemical analysis based on this technique remains challenging. One potential solution to overcome this limitation is the use of isotopically labeled molecules. However, such compounds, particularly their use in two-dimensional (2D) NMR techniques, have rarely been studied. This study presents the results of both experimental and simulated correlation spectroscopy (COSY) on 1-13C-ethanol at 34.38 μT. The strong heteronuclear coupling in this molecule breaks the magnetic equivalence, causing all J-couplings, including homonuclear coupling, to split the 1H spectrum. The obtained COSY spectrum clearly shows the spectral details. Furthermore, we observed that homonuclear coupling between 1H spins generated cross-peaks only when the associated 1H spins were coupled to identical 13C spin states. Our findings demonstrate that a low-field 2D spectrum, even with a moderate spectral line width, can reveal the J-coupling networks of isotopically labeled molecules.

Visualizing a two-state conformational ensemble in stem-loop 3 of the transcriptional regulator 7SK RNA

Structural plasticity is integral to RNA function; however, there are currently few methods to quantitatively resolve RNAs that have multiple structural states. NMR spectroscopy is a powerful approach for resolving conformational ensembles but is size-limited. Chemical probing is well-suited for large RNAs but provides limited structural and no kinetics information. Here, we integrate the two approaches to visualize a two-state conformational ensemble for the central stem-loop 3 (SL3) of 7SK RNA, a critical element for 7SK RNA function in transcription regulation. We find that the SL3 distal end exchanges between two equally populated yet structurally distinct states in both isolated SL3 constructs and full-length 7SK RNA. We rationally designed constructs that lock SL3 into a single state and demonstrate that both chemical probing and NMR data fit to a linear combination of the two states. Comparison of vertebrate 7SK RNA sequences shows conservation of both states, suggesting functional importance. These results provide new insights into 7SK RNA structural dynamics and demonstrate the utility of integrating chemical probing with NMR spectroscopy to gain quantitative insights into RNA conformational ensembles.

RPflex: A Coarse-Grained Network Model for RNA Pocket Flexibility Study

RNA regulates various biological processes, such as gene regulation, RNA splicing, and intracellular signal transduction. RNA’s conformational dynamics play crucial roles in performing its diverse functions. Thus, it is essential to explore the flexibility characteristics of RNA, especially pocket flexibility. Here, we propose a computational approach, RPflex, to analyze pocket flexibility using the coarse-grained network model. We first clustered 3154 pockets into 297 groups by similarity calculation based on the coarse-grained lattice model. Then, we introduced the flexibility score to quantify the flexibility by global pocket features. The results show strong correlations between the flexibility scores and root-mean-square fluctuation (RMSF) values, with Pearson correlation coefficients of 0.60, 0.76, and 0.53 in Testing Sets I–III. Considering both flexibility score and network calculations, the Pearson correlation coefficient was increased to 0.71 in flexible pockets on Testing Set IV. The network calculations reveal that the long-range interaction changes contributed most to flexibility. In addition, the hydrogen bonds in the base–base interactions greatly stabilize the RNA structure, while backbone interactions determine RNA folding. The computational analysis of pocket flexibility could facilitate RNA engineering for biological or medical applications.

Combining NMR, SAXS and SANS to characterize the structure and dynamics of protein complexes

Understanding the structure and dynamics of biological macromolecules is essential to decipher the molecular mechanisms that underlie cellular functions. The description of structure and conformational dynamics often requires the integration of complementary techniques. In this review, we highlight the utility of combining nuclear magnetic resonance (NMR) spectroscopy with small angle scattering (SAS) to characterize these challenging biomolecular systems. NMR can assess the structure and conformational dynamics of multidomain proteins, RNAs and biomolecular complexes. It can efficiently provide information on interaction surfaces, long-distance restraints and relative domain orientations at residue-level resolution. Such information can be readily combined with high-resolution structural data available on subcomponents of biomolecular assemblies. Moreover, NMR is a powerful tool to characterize the dynamics of biomolecules on a wide range of timescales, from nanoseconds to seconds. On the other hand, SAS approaches provide global information on the size and shape of biomolecules and on the ensemble of all conformations present in solution. Therefore, NMR and SAS provide complementary data that are uniquely suited to investigate dynamic biomolecular assemblies. Here, we briefly review the type of data that can be obtained by both techniques and describe different approaches that can be used to combine them to characterize biomolecular assemblies. We then provide guidelines on which experiments are best suited depending on the type of system studied, ranging from fully rigid complexes, dynamic structures that interconvert between defined conformations and systems with very high structural heterogeneity.

NMR spectroscopy, excited states and relevance to problems in cell biology - transient pre-nucleation tetramerization of huntingtin and insights into Huntington's disease

Solution nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for analyzing three-dimensional structure and dynamics of macromolecules at atomic resolution. Recent advances have exploited the unique properties of NMR in exchanging systems to detect, characterize and visualize excited sparsely populated states of biological macromolecules and their complexes, which are only transient. These states are invisible to conventional biophysical techniques, and play a key role in many processes, including molecular recognition, protein folding, enzyme catalysis, assembly and fibril formation. All the NMR techniques make use of exchange between sparsely populated NMR-invisible and highly populated NMR-visible states to transfer a magnetization property from the invisible state to the visible one where it can be easily detected and quantified. There are three classes of NMR experiments that rely on differences in distance, chemical shift or transverse relaxation (molecular mass) between the NMR-visible and -invisible species. Here, I illustrate the application of these methods to unravel the complex mechanism of sub-millisecond pre-nucleation oligomerization of the N-terminal region of huntingtin, encoded by exon-1 of the huntingtin gene, where CAG expansion leads to Huntington's disease, a fatal autosomal-dominant neurodegenerative condition. I also discuss how inhibition of tetramerization blocks the much slower (by many orders of magnitude) process of fibril formation.

Accelerating the Ensemble Convergence of RNA Hairpin Simulations with a Replica Exchange Structure Reservoir

RNA is a key participant in many biological processes, but studies of RNA using computer simulations lag behind those of proteins, largely due to less-developed force fields and the slow dynamics of RNA. Generating converged RNA ensembles for force field development and other studies remains a challenge. In this study, we explore the ability of replica exchange molecular dynamics to obtain well-converged conformational ensembles for two RNA hairpin systems in an implicit solvent. Even for these small model systems, standard REMD remains computationally costly, but coupling to a pre-generated structure library using the reservoir REMD approach provides a dramatic acceleration of ensemble convergence for both model systems. Such precise ensembles could facilitate RNA force field development and validation and applications of simulation to more complex RNA systems. The advantages and remaining challenges of applying R-REMD to RNA are investigated in detail.

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Hidden intermediates in Mango III RNA aptamer folding revealed by pressure perturbation

Fluorescent RNA aptamers have the potential to enable routine quantitation and localization of RNA molecules and serve as models for understanding biologically active aptamers. In recent years, several fluorescent aptamers have been selected and modified to improve their properties, revealing that small changes to the RNA or the ligands can modify significantly their fluorescent properties. Although structural biology approaches have revealed the bound, ground state of several fluorescent aptamers, characterization of low-abundance, excited states in these systems is crucial to understanding their folding pathways. Here we use pressure as an alternative variable to probe the suboptimal states of the Mango III aptamer with both fluorescence and NMR spectroscopy approaches. At moderate KCl concentrations, increasing pressure disrupted the G-quadruplex structure of the Mango III RNA and led to an intermediate with lower fluorescence. These observations indicate the existence of suboptimal RNA structural states that still bind the TO1-biotin fluorophore and moderately enhance fluorescence. At higher KCl concentration as well, the intermediate fluorescence state was populated at high pressure, but the G-quadruplex remained stable at high pressure, supporting the notion of parallel folding and/or binding pathways. These results demonstrate the usefulness of pressure for characterizing RNA folding intermediates.

Exciton Delocalization and Scaffold Stability in Bridged Nucleotide-Substituted, DNA Duplex-Templated Cyanine Aggregates

Molecular excitons play a foundational role in chromophore aggregates found in light-harvesting systems and offer potential applications in engineered excitonic systems. Controlled aggregation of chromophores to promote exciton delocalization has been achieved by covalently tethering chromophores to deoxyribonucleic acid (DNA) scaffolds. Although many studies have documented changes in the optical properties of chromophores upon aggregation using DNA scaffolds, more limited work has investigated how structural modifications of DNA via bridged nucleotides and chromophore covalent attachment impact scaffold stability as well as the configuration and optical behavior of attached aggregates. Here we investigated the impact of two types of bridged nucleotides, LNA and BNA, as a structural modification of duplex DNA-templated cyanine (Cy5) aggregates. The bridged nucleotides were incorporated in the domain of one to four Cy5 chromophores attached between adjacent bases of a DNA duplex. We found that bridged nucleotides increase the stability of DNA scaffolds carrying Cy5 aggregates in comparison with natural nucleotides in analogous constructs. Exciton coupling strength and delocalization in Cy5 aggregates were evaluated via steady-state absorption, circular dichroism, and theoretical modeling. Replacing natural nucleotides with bridged nucleotides resulted in a noticeable increase in the coupling strength (≥10 meV) between chromophores and increased H-like stacking behavior (i.e., more face-to-face stacking). Our results suggest that bridged nucleotides may be useful for increasing scaffold stability and coupling between DNA templated chromophores.

A Carbazole-Derived Nitroxide That Is an Analogue of Cytidine: A Rigid Spin Label for DNA and RNA

A variety of semirigid and rigid spin labels comprise a valuable arsenal for measurements of biomolecular structures and dynamics by electron paramagnetic resonance (EPR) spectroscopy. Here, we report the synthesis and characterization of rigid spin labels Ċ and Ċm for DNA and RNA, respectively, that are carbazole-derived nitroxides and analogues of cytidine. Ċ and Ċm were converted to their phosphoramidites and used for their incorporation into oligonucleotides by solid-phase synthesis. Analysis of Ċ and Ċm by single-crystal X-ray crystallography verified their identity and showed little deviation from planarity of the nucleobase. Analysis of the continuous-wave (CW) EPR spectra of the spin-labeled DNA and RNA duplexes confirmed their incorporation into the nucleic acids and the line-shape was characteristic of rigid spin labels. Circular dichroism (CD) and thermal denaturation studies of the Ċ-labeled DNAs and Ċm-labeled RNAs indicated that the labels are nonperturbing of duplex structure.

Ensemble Switching Unveils a Kinetic Rheostat Mechanism of the Eukaryotic Thiamine Pyrophosphate Riboswitch

Thiamine pyrophosphate (TPP) riboswitches regulate thiamine metabolism by inhibiting the translation of enzymes essential to thiamine synthesis pathways upon binding to thiamine pyrophosphate in cells across all domains of life. Recent work on the Arabidopsis thaliana TPP riboswitch suggests a multi-step TPP binding process involving multiple riboswitch configurational ensembles and that Mg2+ dependence underlies the mechanism of TPP recognition and subsequent transition to the expression-inhibiting state of the aptamer domain followed by changes in the expression platform. However, details of the relationship between TPP riboswitch conformational changes and interactions with TPP and Mg2+ ¬¬in the aptamer domain constituting this mechanism are unknown. Therefore, we integrated single-molecule multiparameter fluorescence and force spectroscopy with atomistic molecular dynamics simulations and found that conformational transitions within the aptamer domain's sensor helices associated with TPP and Mg2+ ligand binding occurred between at least five different ensembles on timescales ranging from µs to ms. These dynamics are orders of magnitude faster than the 10 second-timescale folding kinetics associated with expression-state switching in the switch sequence. Together, our results show that a TPP and Mg2+ dependent mechanism determines dynamic configurational state ensemble switching of the aptamer domain's sensor helices that regulates the stability of the switch helix, which ultimately may lead to the expression-inhibiting state of the riboswitch. Additionally, we propose that two pathways exist for ligand recognition and that this mechanism underlies a kinetic rheostat-like behavior of the Arabidopsis thaliana TPP riboswitch.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA's diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

Dynamics-Function Analysis in Catalytic RNA Using NMR Spin Relaxation and Conformationally Restricted Nucleotides

A full understanding of biomolecular function requires an analysis of both the dynamic properties of the system of interest and the identification of those dynamics that are required for function. We describe NMR methods based on metabolically directed specific isotope labeling for the identification of molecular disorder and/or conformational transitions on the RNA backbone ribose groups. These analyses are complemented by the use of synthetic covalently modified nucleotides constrained to a single sugar pucker, which allow functional assessment of dynamics by selectively removing a minor conformer identified by NMR from the structural ensemble.

A perspective on the molecular simulation of DNA from structural and functional aspects

As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.

Divalent Cation Dependence Enhances Dopamine Aptamer Biosensing

Oligonucleotide receptors (aptamers), which change conformation upon target recognition, enable electronic biosensing under high ionic-strength conditions when coupled to field-effect transistors (FETs). Because highly negatively charged aptamer backbones are influenced by ion content and concentration, biosensor performance and target sensitivities were evaluated under application conditions. For a recently identified dopamine aptamer, physiological concentrations of Mg²⁺ and Ca²⁺ in artificial cerebrospinal fluid produced marked potentiation of dopamine FET-sensor responses. By comparison, divalent cation-associated signal amplification was not observed for FET sensors functionalized with a recently identified serotonin aptamer or a previously reported dopamine aptamer. Circular dichroism spectroscopy revealed Mg²⁺- and Ca²⁺-induced changes in target-associated secondary structure for the new dopamine aptamer, but not the serotonin aptamer nor the old dopamine aptamer. Thioflavin T displacement corroborated the Mg²⁺ dependence of the new dopamine aptamer for target detection. These findings imply allosteric binding interactions between divalent cations and dopamine for the new dopamine aptamer. Developing and testing sensors in ionic environments that reflect intended applications are best practices for identifying aptamer candidates with favorable attributes and elucidating sensing mechanisms.

NMR spectroscopy captures the essential role of dynamics in regulating biomolecular function

Biomolecules are in constant motion. To understand how they function, and why malfunctions can cause disease, it is necessary to describe their three-dimensional structures in terms of dynamic conformational ensembles. Here, we demonstrate how nuclear magnetic resonance (NMR) spectroscopy provides an essential, dynamic view of structural biology that captures biomolecular motions at atomic resolution. We focus on examples that emphasize the diversity of biomolecules and biochemical applications that are amenable to NMR, such as elucidating functional dynamics in large molecular machines, characterizing transient conformations implicated in the onset of disease, and obtaining atomic-level descriptions of intrinsically disordered regions that make weak interactions involved in liquid-liquid phase separation. Finally, we discuss the pivotal role that NMR has played in driving forward our understanding of the biomolecular dynamics-function paradigm.

Cooperativity and Allostery in RNA Systems

Allostery is among the most basic biological principles employed by biological macromolecules to achieve a biologically active state in response to chemical cues. Although initially used to describe the impact of small molecules on the conformation and activity of protein enzymes, the definition of this term has been significantly broadened to describe long-range conformational change of macromolecules in response to small or large effectors. Such a broad definition could be applied to RNA molecules, which do not typically serve as protein-free cellular enzymes but fold and form macromolecular assemblies with the help of various ligand molecules, including ions and proteins. Ligand-induced allosteric changes in RNA molecules are often accompanied by cooperative interactions between RNA and its ligand, thus streamlining the folding and assembly pathways. This chapter provides an overview of the interplay between cooperativity and allostery in RNA systems and outlines methods to study these two biological principles.

Slow motions in A·T rich DNA sequence

In free B-DNA, slow (microsecond-to-millisecond) motions that involve equilibrium between Watson-Crick (WC) and Hoogsteen (HG) base-pairing expand the DNA dynamic repertoire that could mediate DNA-protein assemblies. R_1ρ relaxation dispersion NMR methods are powerful tools to capture such slow conformational exchanges in solution using ¹³C/¹⁵ N labelled DNA. Here, these approaches were applied to a dodecamer containing a TTAAA element that was assumed to facilitate nucleosome formation. NMR data and inferred exchange parameters assign HG base pairs as the minor, transient conformers specifically observed in three successive A·T base pairs forming the TAA·TTA segment. The abundance of these HG A·T base pairs can be up to 1.2% which is high compared to what has previously been observed. Data analyses support a scenario in which the three adenines undergo non-simultaneous motions despite their spatial proximity, thus optimising the probability of having one HG base pair in the TAA·TTA segment. Finally, revisiting previous NMR data on H2 resonance linewidths on the basis of our results promotes the idea of there being a special propensity of A·T base pairs in TAA·TTA tracts to adopt HG pairing. In summary, this study provides an example of a DNA functional element submitted to slow conformational exchange. More generally, it strengthens the importance of the role of the DNA sequence in modulating its dynamics, over a nano- to milli-second time scale.

Advances in RNA 3D Structure Modeling Using Experimental Data

RNA is a unique bio-macromolecule that can both record genetic information and perform biological functions in a variety of molecular processes, including transcription, splicing, translation, and even regulating protein function. RNAs adopt specific three-dimensional conformations to enable their functions. Experimental determination of high-resolution RNA structures using x-ray crystallography is both laborious and demands expertise, thus, hindering our comprehension of RNA structural biology. The computational modeling of RNA structure was a milestone in the birth of bioinformatics. Although computational modeling has been greatly improved over the last decade showing many successful cases, the accuracy of such computational modeling is not only length-dependent but also varies according to the complexity of the structure. To increase credibility, various experimental data were integrated into computational modeling. In this review, we summarize the experiments that can be integrated into RNA structure modeling as well as the computational methods based on these experimental data. We also demonstrate how computational modeling can help the experimental determination of RNA structure. We highlight the recent advances in computational modeling which can offer reliable structure models using high-throughput experimental data.

On an enhancement of RNA probing data using information theory

Identifying the secondary structure of an RNA is crucial for understanding its diverse regulatory functions. This paper focuses on how to enhance target identification in a Boltzmann ensemble of structures via chemical probing data. We employ an information-theoretic approach to solve the problem, via considering a variant of the Rényi-Ulam game. Our framework is centered around the ensemble tree, a hierarchical bi-partition of the input ensemble, that is constructed by recursively querying about whether or not a base pair of maximum information entropy is contained in the target. These queries are answered via relating local with global probing data, employing the modularity in RNA secondary structures. We present that leaves of the tree are comprised of sub-samples exhibiting a distinguished structure with high probability. In particular, for a Boltzmann ensemble incorporating probing data, which is well established in the literature, the probability of our framework correctly identifying the target in the leaf is greater than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$90\%$$\end{document}90%.

Structural transitions in the RNA 7SK 5' hairpin and their effect on HEXIM binding

7SK RNA, as part of the 7SK ribonucleoprotein complex, is crucial to the regulation of transcription by RNA-polymerase II, via its interaction with the positive transcription elongation factor P-TEFb. The interaction is induced by binding of the protein HEXIM to the 5′ hairpin (HP1) of 7SK RNA. Four distinct structural models have been obtained experimentally for HP1. Here, we employ computational methods to investigate the relative stability of these structures, transitions between them, and the effects of mutations on the observed structural ensembles. We further analyse the results with respect to mutational binding assays, and hypothesize a mechanism for HEXIM binding. Our results indicate that the dominant structure in the wild type exhibits a triplet involving the unpaired nucleotide U40 and the base pair A43-U66 in the GAUC/GAUC repeat. This conformation leads to an open major groove with enough potential binding sites for peptide recognition. Sequence mutations of the RNA change the relative stability of the different structural ensembles. Binding affinity is consequently lost if these changes alter the dominant structure.

Exposing Hidden High-Affinity RNA Conformational States

RNA recognition frequently results in conformational changes that optimize intermolecular binding. As a consequence, the overall binding affinity of RNA to its binding partners depends not only on the intermolecular interactions formed in the bound state but also on the energy cost associated with changing the RNA conformational distribution. Measuring these "conformational penalties" is, however, challenging because bound RNA conformations tend to have equilibrium populations in the absence of the binding partner that fall outside detection by conventional biophysical methods. In this study we employ as a model system HIV-1 TAR RNA and its interaction with the ligand argininamide (ARG), a mimic of TAR's cognate protein binding partner, the transactivator Tat. We use NMR chemical shift perturbations and relaxation dispersion in combination with Bayesian inference to develop a detailed thermodynamic model of coupled conformational change and ligand binding. Starting from a comprehensive 12-state model of the equilibrium, we estimate the energies of six distinct detectable thermodynamic states that are not accessible by currently available methods. Our approach identifies a minimum of four RNA intermediates that differ in terms of the TAR conformation and ARG occupancy. The dominant bound TAR conformation features two bound ARG ligands and has an equilibrium population in the absence of ARG that is below detection limit. Consequently, even though ARG binds to TAR with an apparent overall weak affinity (Kdapp ≈ 0.2 mM), it binds the prefolded conformation with a Kd in the nM range. Our results show that conformational penalties can be major determinants of RNA-ligand binding affinity as well as a source of binding cooperativity, with important implications for a predictive understanding of how RNA is recognized and for RNA-targeted drug discovery.

Probing excited conformational states of nucleic acids by nitrogen CEST NMR spectroscopy

Characterizing low-populated and short-lived excited conformational states has become increasingly important for understanding mechanisms of RNA function. Interconversion between RNA ground and excited conformational states often involves base pairing rearrangements that lead to changes in the hydrogen-bond network. Here, we present two 15N chemical exchange saturation transfer (CEST) NMR experiments that utilize protonated and non-protonated nitrogens, which are key hydrogen-bond donors and acceptors, for characterizing excited conformational states in RNA. We demonstrated these approaches on the B. Cereus fluoride riboswitch, where 15N CEST profiles complement 13C CEST profiles in depicting a potential pathway for ligand-dependent allosteric regulation of the excited conformational state of the fluoride riboswitch.

Asymmetric base-pair opening drives helicase unwinding dynamics

The opening of a Watson-Crick double helix is required for crucial cellular processes, including replication, repair, and transcription. It has long been assumed that RNA or DNA base pairs are broken by the concerted symmetric movement of complementary nucleobases. By analyzing thousands of base-pair opening and closing events from molecular simulations, here, we uncover a systematic stepwise process driven by the asymmetric flipping-out probability of paired nucleobases. We demonstrate experimentally that such asymmetry strongly biases the unwinding efficiency of DNA helicases toward substrates that bear highly dynamic nucleobases, such as pyrimidines, on the displaced strand. Duplex substrates with identical thermodynamic stability are thus shown to be more easily unwound from one side than the other, in a quantifiable and predictable manner. Our results indicate a possible layer of gene regulation coded in the direction-dependent unwindability of the double helix.

Systems NMR: single-sample quantification of RNA, proteins and metabolites for biomolecular network analysis

Cellular behavior is controlled by the interplay of diverse biomolecules. Most experimental methods, however, can monitor only a single molecule class or reaction type at a time. We developed an in vitro Nuclear Magnetic Resonance spectroscopy (NMR) approach, which permitted dynamic quantification of an entire “heterotypic” network – simultaneously monitoring three distinct molecule classes (metabolites, proteins, RNA) and all elementary reaction types (bimolecular interactions, catalysis, unimolecular changes). Focusing on an 8-reaction co-transcriptional RNA folding network, in a single sample we recorded over 35 time-points with over 170 observables each, and accurately determined 5 core reaction constants in multiplex. This reconstruction revealed unexpected cross-talk between the different reactions. We further observed dynamic phase-separation in a system of five distinct RNA binding domains in the course of the RNA transcription reaction. Our Systems NMR approach provides a deeper understanding of biological network dynamics by combining the dynamic resolution of biochemical assays and the multiplexing ability of “omics”.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

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.

An NMR View of Protein Dynamics in Health and Disease

Biological molecules are often highly dynamic, and this flexibility can be critical for function. The large range of sampled timescales and the fact that many of the conformers that are continually explored are only transiently formed and sparsely populated challenge current biophysical approaches. Solution nuclear magnetic resonance (NMR) spectroscopy has emerged as a powerful method for characterizing biomolecular dynamics in detail, even in cases where excursions involve short-lived states. Here, we briefly review a number of NMR experiments for studies of biomolecular dynamics on the microsecond-to-second timescale and focus on applications to protein and nucleic acid systems that clearly illustrate the functional relevance of motion in both health and disease.

Advances that facilitate the study of large RNA structure and dynamics by nuclear magnetic resonance spectroscopy

The characterization of functional yet nonprotein coding (nc) RNAs has expanded the role of RNA in the cell from a passive player in the central dogma of molecular biology to an active regulator of gene expression. The misregulation of ncRNA function has been linked with a variety of diseases and disorders ranging from cancers to neurodegeneration. However, a detailed molecular understanding of how ncRNAs function has been limited; due, in part, to the difficulties associated with obtaining high-resolution structures of large RNAs. Tertiary structure determination of RNA as a whole is hampered by various technical challenges, all of which are exacerbated as the size of the RNA increases. Namely, RNAs tend to be highly flexible and dynamic molecules, which are difficult to crystallize. Biomolecular nuclear magnetic resonance (NMR) spectroscopy offers a viable alternative to determining the structure of large RNA molecules that do not readily crystallize, but is itself hindered by some technical limitations. Recently, a series of advancements have allowed the biomolecular NMR field to overcome, at least in part, some of these limitations. These advances include improvements in sample preparation strategies as well as methodological improvements. Together, these innovations pave the way for the study of ever larger RNA molecules that have important biological function. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Towards an understanding of RNA structural modalities: a riboswitch case study

A riboswitch is a type of RNA molecule that regulates important biological functions by changing structure, typically under ligand-binding. We assess the extent that these ligand-bound structural alternatives are present in the Boltzmann sample, a standard RNA secondary structure prediction method, for three riboswitch test cases. We use the cluster analysis tool RNAStructProfiling to characterize the different modalities present among the suboptimal structures sampled. We compare these modalities to the putative base pairing models obtained from independent experiments using NMR or fluorescence spectroscopy. We find, somewhat unexpectedly, that profiling the Boltzmann sample captures evidence of ligand-bound conformations for two of three riboswitches studied. Moreover, this agreement between predicted modalities and experimental models is consistent with the classification of riboswitches into thermodynamic versus kinetic regulatory mechanisms. Our results support cluster analysis of Boltzmann samples by RNAStructProfiling as a possible basis for de novo identification of thermodynamic riboswitches, while highlighting the challenges for kinetic ones.

Overhauser Dynamic Nuclear Polarization for the Study of Hydration Dynamics, Explained

We outline the physical properties of hydration water that are captured by Overhauser Dynamic Nuclear Polarization (ODNP) relaxometry and explore the insights that ODNP yields about the water and the surface that this water is coupled to. As ODNP relies on the pairwise cross-relaxation between the electron spin of a spin probe and a proton nuclear spin of water, it captures the dynamics of single-particle diffusion of an ensemble of water molecules moving near the spin probe. ODNP principally utilizes the same physics as other nuclear magnetic resonance (NMR) relaxometry (i.e., relaxation measurement) techniques. However, in ODNP, electron paramagnetic resonance (EPR) excites the electron spins probes and their high net polarization acts as a signal amplifier. Furthermore, it renders ODNP parameters highly sensitive to water moving at rates commensurate with the EPR frequency of the spin probe (typically 10GHz). Also, ODNP selectively enhances the NMR signal contributions of water moving within close proximity to the spin label. As a result, ODNP can capture ps-ns movements of hydration waters with high sensitivity and locality, even in samples with protein concentrations as dilute as 10 µM. To date, the utility of the ODNP technique has been demonstrated for two major applications: the characterization of the spatial variation in the properties of the hydration layer of proteins or other surfaces displaying topological diversity, and the identification of structural properties emerging from highly disordered proteins and protein domains. The former has been shown to correlate well with the properties of hydration water predicted by MD simulations and has been shown capable of evaluating the hydrophilicity or hydrophobicity of a surface. The latter has been demonstrated for studies of an interhelical loop of proteorhodopsin, the partial structure of α-synuclein embedded at the lipid membrane surface, incipient structures adopted by tau proteins en route to fibrils, and the structure and hydration profile of a transmembrane peptide. This chapter focuses on offering a mechanistic understanding of the ODNP measurement and the molecular dynamics encoded in the ODNP parameters. In particular, it clarifies how the electron-nuclear dipolar coupling encodes information about the molecular dynamics in the nuclear spin self-relaxation and, more importantly, the electron-nuclear spin cross-relaxation rates. The clarification of the molecular dynamics underlying ODNP should assist in establishing a connection to theory and computer simulation that will offer far richer interpretations of ODNP results in future studies.

The RNA Base-Pairing Problem and Base-Pairing Solutions

RNA molecules are folded into structures and complexes to perform a wide variety of functions. Determination of RNA structures and their interactions is a fundamental problem in RNA biology. Most RNA molecules in living cells are large and dynamic, posing unique challenges to structure analysis. Here we review progress in RNA structure analysis, focusing on methods that use the "cross-link, proximally ligate, and sequence" principle for high-throughput detection of base-pairing interactions in living cells. Beginning with a comparison of commonly used methods in structure determination and a brief historical account of psoralen cross-linking studies, we highlight the important features of cross-linking methods and new biological insights into RNA structures and interactions from recent studies. Further improvement of these cross-linking methods and application to previously intractable problems will shed new light on the mechanisms of the "modern RNA world."

Coupling between conformational dynamics and catalytic function at the active site of the lead-dependent ribozyme

In common with other self-cleaving RNAs, the lead-dependent ribozyme (leadzyme) undergoes dynamic fluctuations to a chemically activated conformation. We explored the connection between conformational dynamics and self-cleavage function in the leadzyme using a combination of NMR spin-relaxation analysis of ribose groups and conformational restriction via chemical modification. The functional studies were performed with a North-methanocarbacytidine modification that prevents fluctuations to C2'-endo conformations while maintaining an intact 2'-hydroxyl nucleophile. Spin-relaxation data demonstrate that the active-site Cyt-6 undergoes conformational exchange attributed to sampling of a minor C2'-endo state with an exchange lifetime on the order of microseconds to tens of microseconds. A conformationally restricted species in which the fluctuations to the minor species are interrupted shows a drastic decrease in self-cleavage activity. Taken together, these data indicate that dynamic sampling of a minor species at the active site of this ribozyme, and likely of related naturally occurring motifs, is strongly coupled to catalytic function. The combination of NMR dynamics analysis with functional probing via conformational restriction is a general methodology for dissecting dynamics-function relationships in RNA.

Extending the Applicability of Exact Nuclear Overhauser Enhancements to Large Proteins and RNA

Distance-dependent nuclear Overhauser enhancements (NOEs) are one of the most popular and important experimental restraints for calculating NMR structures. Despite this, they are mostly employed as semiquantitative upper distance bounds, and this discards the wealth of information that is encoded in the cross-relaxation rate constant. Information that is lost includes exact distances between protons and dynamics that occur on the sub-millisecond timescale. Our recently introduced exact measurement of the NOE (eNOE) requires little additional experimental effort relative to other NMR observables. So far, we have used eNOEs to calculate multistate ensembles of proteins up to approximately 150 residues. Here, we briefly revisit eNOE methodology and present two new directions for the use of eNOEs: applications to large proteins and RNA.

High-resolution small RNA structures from exact nuclear Overhauser enhancement measurements without additional restraints

RNA not only translates the genetic code into proteins, but also carries out important cellular functions. Understanding such functions requires knowledge of the structure and dynamics at atomic resolution. Almost half of the published RNA structures have been solved by nuclear magnetic resonance (NMR). However, as a result of severe resonance overlap and low proton density, high-resolution RNA structures are rarely obtained from nuclear Overhauser enhancement (NOE) data alone. Instead, additional semi-empirical restraints and labor-intensive techniques are required for structural averages, while there are only a few experimentally derived ensembles representing dynamics. Here we show that our exact NOE (eNOE) based structure determination protocol is able to define a 14-mer UUCG tetraloop structure at high resolution without other restraints. Additionally, we use eNOEs to calculate a two-state structure, which samples its conformational space. The protocol may open an avenue to obtain high-resolution structures of small RNA of unprecedented accuracy with moderate experimental efforts.

Studying sparsely populated conformational states in RNA combining chemical synthesis and solution NMR spectroscopy

Using chemical synthesis and solution NMR spectroscopy, RNA structural ensembles including a major ground state and minor populated excited states can be studied at atomic resolution. In this work, atom-specific ¹³C labeled RNA building blocks - a 5-¹³C-uridine and a 2,8-¹³C₂-adenosine building block - are used to introduce isolated ¹³C-¹H-spin topologies into a target RNA to probe such structural ensembles via NMR spectroscopy. First, the 5-¹³C-uridine 2'-O-TBDMS-phosphoramidite building block was introduced into a 21 nucleotide (nt) tP5c stem construct of the tP5abc subdomain of the Tetrahymena group I ribozyme. Then, the 2,8-¹³C₂-adenosine 2'-O-TBDMS-phosphoramidite building block was incorporated into a 9 kDa and a 15 kD construct derived from the epsilon (ε) RNA element of the duck Hepatitis B virus. The 2,8-¹³C₂-adenosine resonances of the 9 kDa 28 nt sequence could be mapped to the full-length 53 nt construct. The isolated NMR active nuclei pairs were used to probe for low populated excited states (<10%) via ¹³C-Carr-Purcell-Meiboom-Gill (CPMG)-relaxation dispersion NMR spectroscopy. The ¹³C-CPMG relaxation dispersion experiment recapitulated a secondary structure switching event in the P5c hairpin of the group I intron construct previously revealed by ¹⁵N relaxation dispersion experiments. In the ε-HBV RNA an unfolding event occurring on the millisecond time scale was found in the upper stem in-line with earlier observations. This unpaired conformational state is presumed to be important for the binding of the epsilon reverse transcriptase (RT) enzyme. Thus, a full description of an RNA's folding landscape helps to obtain a deeper understanding of its function, as these high energy conformational states often represent functionally important intermediates involved in (un)folding or ribozyme catalysis.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.