Molecular recognition and function of riboswitches

RNA cis-regulators are important for Streptococcus pneumoniae in vivo success

Bacteria have evolved complex transcriptional regulatory networks, as well as many diverse regulatory strategies at the RNA level, to enable more efficient use of metabolic resources and a rapid response to changing conditions. However, most RNA-based regulatory mechanisms are not well conserved across different bacterial species despite controlling genes important for virulence or essential biosynthetic processes. Here, we characterize the activity of, and assess the fitness benefit conferred by, twelve cis-acting regulatory RNAs (including several riboswitches and a T-box), in the opportunistic pathogen Streptococcus pneumoniae TIGR4. By evaluating native locus mutants of each regulator that result in constitutively active or repressed expression, we establish that growth defects in planktonic culture are associated with constitutive repression of gene expression, while constitutive activation of gene expression is rarely deleterious. In contrast, in mouse nasal carriage and pneumonia models, strains with either constitutively active and repressed gene expression are significantly less fit than matched control strains. Furthermore, two RNA-regulated pathways, FMN synthesis/transport and pyrimidine synthesis/transport display exceptional sensitivity to mis-regulation or constitutive gene repression in both planktonic culture and in vivo environments. Thus, despite lack of obvious phenotypes associated with constitutive gene expression in vitro, the fitness benefit conferred on bacteria via fine-tuned metabolic regulation through cis-acting regulatory RNAs is substantial in vivo, and therefore easily sufficient to drive the evolution and maintenance of diverse RNA regulatory mechanisms.

Intricate 3D architecture of a DNA mimic of GFP

Numerous studies have shown how RNA molecules can adopt elaborate three-dimensional (3D) architectures1-3. By contrast, whether DNA can self-assemble into complex 3D folds capable of sophisticated biochemistry, independent of protein or RNA partners, has remained mysterious. Lettuce is an in vitro-evolved DNA molecule that binds and activates4 conditional fluorophores derived from GFP. To extend previous structural studies5,6 of fluorogenic RNAs, GFP and other fluorescent proteins7 to DNA, we characterize Lettuce-fluorophore complexes by X-ray crystallography and cryogenic electron microscopy. The results reveal that the 53-nucleotide DNA adopts a four-way junction (4WJ) fold. Instead of the canonical L-shaped or H-shaped structures commonly seen8 in 4WJ RNAs, the four stems of Lettuce form two coaxial stacks that pack co-linearly to form a central G-quadruplex in which the fluorophore binds. This fold is stabilized by stacking, extensive nucleobase hydrogen bonding-including through unusual diagonally stacked bases that bridge successive tiers of the main coaxial stacks of the DNA-and coordination of monovalent and divalent cations. Overall, the structure is more compact than many RNAs of comparable size. Lettuce demonstrates how DNA can form elaborate 3D structures without using RNA-like tertiary interactions and suggests that new principles of nucleic acid organization will be forthcoming from the analysis of complex DNAs.

Ionic Cooperativity between Lysine and Potassium in the Lysine Riboswitch: Single-Molecule Kinetic and Thermodynamic Studies

Functionality in many biological systems, including proteins and nucleic acid structures, including protein and nucleic acid riboswitch structures, can depend on cooperative kinetic behavior between multiple small molecule ligands. In this work, single-molecule FRET data on the Bacillus subtilis lysine riboswitch reveals that affinity for the cognate lysine ligand increases significantly with K⁺, providing evidence for synergism between lysine/K⁺ binding to the aptamer and successful folding of the riboswitch. To describe/interpret this more complex kinetic scenario, we explore the conventional 4-state ("square") model for aptamer binding as a function of K⁺. Extension into this additional dimension generates a novel "cube" model for riboswitch folding dynamics with respect to lysine/K⁺ binding, revealing that riboswitch folding (k_fold) and unfolding (k_unfold) rate constants increase and decrease dramatically with K⁺, respectively. Furthermore, temperature-dependent single-molecule kinetic studies indicate that the presence of K⁺ entropically enhances the transition state barrier to folding but partially compensates for this by increasing the overall exothermicity for lysine binding. We rationalize this behavior as evidence that K⁺ facilitates hydrogen bonding between the negatively charged carboxyl group of lysine and the RNA, increasing structural rigidity and lowering entropy in the binding pocket. Finally, we explore the effects of cation size with Na⁺ and Cs⁺ studies to demonstrate that K⁺ is optimally suited for bridging interactions between lysine and the riboswitch aptamer domain. Regulation of lysine production and transport, dictated by the riboswitch's ability to recognize and bind lysine, is therefore intimately tied to the presence of K⁺ in the binding pocket and is strongly modulated by local cation conditions. The results suggest an increase in lysine riboswitch functionality by sensitivity to additional species in the cellular riboswitch environment.

A variant of guanidine-IV riboswitches exhibits evidence of a distinct ligand specificity

Riboswitches are regulatory RNAs that specifically bind a small molecule or ion. Like metabolite-binding proteins, riboswitches can evolve new ligand specificities, and some examples of this phenomenon have been validated. As part of work based on comparative genomics to discover novel riboswitches, we encountered a candidate riboswitch with striking similarities to the recently identified guanidine-IV riboswitch. This candidate riboswitch, the Gd4v motif, is predicted in four distinct bacterial phyla, thus almost as widespread as the guanidine-IV riboswitch. Bioinformatic and experimental analysis suggest that the Gd4v motif is a riboswitch that binds a ligand other than guanidine. It is found associated with gene classes that differ from genes regulated by confirmed guanidine riboswitches. In inline-probing assays, we showed that free guanidine binds only weakly to one of the tested sequences of the variant. Further tested compounds did not show binding, attenuation of transcription termination, or activation of a genetic reporter construct. We characterized an N-acetyltransferase frequently associated with the Gd4v motif and compared its substrate preference to an N-acetyltransferase that occurs under control of guanidine-IV riboswitches. The substrates of this Gd4v-motif-associated enzyme did not show activity for Gd4v RNA binding or transcription termination. Hence, the ligand of the candidate riboswitch motif remains unidentified. The variant RNA motif is predominantly found in gut metagenome sequences, hinting at a ligand that is highly relevant in this environment. This finding is a first step to determining the identity of this unknown ligand, and understanding how guanidine-IV-riboswitch-like structures can evolve to bind different ligands.

The promise of cryo-EM to explore RNA structural dynamics

Conformational dynamics are essential to macromolecular function. This is certainly true of RNA, whose ability to undergo programmed conformational dynamics is essential to create and regulate complex biological processes. However, methods to easily and simultaneously interrogate both the structure and conformational dynamics of fully functional RNAs in isolation and in complex with proteins have not historically been available. Due to its ability to image and classify single particles, cryogenic electron microscopy (cryo-EM) has the potential to address this gap and may be particularly amenable to exploring structural dynamics within the three-dimensional folds of biologically active RNAs. We discuss the possibilities and current limitations of applying cryo-EM to simultaneously study RNA structure and conformational dynamics, and present one example that illustrates this (as of yet) not fully realized potential.

Incorporation and Utility of a Responsive Ribonucleoside Analogue in Probing the Conformation of a Viral RNA Motif by Fluorescence and 19 F NMR Spectroscopy

Development of versatile probes that can enable the study of different conformations and recognition properties of therapeutic nucleic acid motifs by complementing biophysical techniques can greatly aid nucleic acid analysis and therapy. Here, we report the design, synthesis and incorporation of an environment-sensitive ribonucleoside analogue, which serves as a two-channel biophysical platform to investigate RNA structure and recognition by fluorescence and ¹⁹ F NMR spectroscopy techniques. The nucleoside analogue is based on a 5-fluorobenzofuran-uracil core and its fluorescence and ¹⁹ F NMR chemical shifts are highly sensitive to changes in solvent polarity and viscosity. Notably, the modified ribonucleotide and phosphoramidite substrates can be efficiently incorporated into RNA oligonucleotides (ONs) by in vitro transcription and standard solid-phase ON synthesis protocol, respectively. Fluorescence and ¹⁹ F readouts of the nucleoside incorporated into model RNA ONs are sensitive to the neighbouring base environment. The responsiveness of the probe was aptly utilized in detecting and quantifying the metal ion-induced conformational change in an internal ribosome entry site RNA motif of hepatitis C virus, which is an important therapeutic target. Taken together, our probe is a good addition to the nucleic acid analysis toolbox with the advantage that it can be used to study nucleic acid conformation and recognition simultaneously by two biophysical techniques.

Modeling Noncanonical RNA Base Pairs by a Coarse-Grained IsRNA2 Model

Noncanonical base pairs contribute crucially to the three-dimensional architecture of large RNA molecules; however, how to accurately model them remains an open challenge in RNA 3D structure prediction. Here, we report a promising coarse-grained (CG) IsRNA2 model to predict noncanonical base pairs in large RNAs through molecular dynamics simulations. By introducing a five-bead per nucleotide CG representation to reserve the three interacting edges of nucleobases, IsRNA2 accurately models various base-pairing interactions, including both canonical and noncanonical base pairs. A benchmark test indicated that IsRNA2 achieves a comparable performance to the atomic model in de novo modeling of noncanonical RNA structures. In addition, IsRNA2 was able to refine the 3D structure predictions for large RNAs in RNA-puzzle challenges. Finally, the graphics processing unit acceleration was introduced to speed up the sampling efficiency in IsRNA2 for very large RNA molecules. Therefore, the CG IsRNA2 model reported here offers a reliable approach to predict the structures and dynamics of large RNAs.

A natural riboswitch scaffold with self-methylation activity

Methylation is a prevalent post-transcriptional modification encountered in coding and non-coding RNA. For RNA methylation, cells use methyltransferases and small organic substances as methyl-group donors, such as S-adenosylmethionine (SAM). SAM and other nucleotide-derived cofactors are viewed as evolutionary leftovers from an RNA world, in which riboswitches have regulated, and ribozymes have catalyzed essential metabolic reactions. Here, we disclose the thus far unrecognized direct link between a present-day riboswitch and its inherent reactivity for site-specific methylation. The key is O⁶-methyl pre-queuosine (m⁶preQ₁), a potentially prebiotic nucleobase which is recognized by the native aptamer of a preQ₁ class I riboswitch. Upon binding, the transfer of the ligand’s methyl group to a specific cytidine occurs, installing 3-methylcytidine (m³C) in the RNA pocket under release of pre-queuosine (preQ₁). Our finding suggests that nucleic acid-mediated methylation is an ancient mechanism that has offered an early path for RNA epigenetics prior to the evolution of protein methyltransferases. Furthermore, our findings may pave the way for the development of riboswitch-descending methylation tools based on rational design as a powerful alternative to in vitro selection approaches.

In humans, protein methyltransferase is responsible for RNA methylation using S-adenosylmethionine as a methyl group donor. Here the authors report a self-methylation activity of a bacterial riboswitch.

Synthesis of DNA and RNA Oligonucleotides Containing a Dual-Purpose Selenium-Modified Fluorescent Nucleoside Probe

Development of efficient tools that would enable direct correlation of nucleic acid structure and recognition in solution and in solid state at atomic resolution is highly desired. In this context, we recently developed dual-purpose nucleoside probes made of a 5-selenophene-modified uracil core, which serves both as a conformation-sensitive fluorophore and as an X-ray crystallography phasing agent. In this article, we provide a detailed synthetic procedure to synthesize the phosphoramidites of 5-selenophene-modified 2'-deoxyuridine and 5-selenophene-modified uridine analogs. We also describe their site-specific incorporation into therapeutically relevant DNA and RNA oligonucleotide motifs by an automated solid support synthesis protocol. The dual-purpose and minimally invasive nature of the probes enables efficient analysis of the conformation and ligand binding abilities of bacterial decoding site RNA (A-site) and G-quadruplex structures of the human telomeric overhang in real time by fluorescence and in 3D by X-ray crystallography. © 2020 by John Wiley & Sons, Inc. Basic Protocol 1: Synthesis of 5-selenophene-2'-deoxyuridine 2 and its phosphoramidite 5 Support Protocol 1: Synthesis of 2-(tri-n-butylstannyl) selenophene Support Protocol 2: Synthesis of 5'-O-DMT-protected 5-iodo-2'-deoxyuridine 3 Basic Protocol 2: Synthesis of 5-selenophene-modified uridine 7 and its phosphoramidite 11 Basic Protocol 3: Synthesis of DNA oligonucleotides containing 5-selenophene-modified 2'-deoxyuridine 2 Basic Protocol 4: Synthesis of an RNA oligonucleotide containing 5-selenophene-modified uridine 7.

The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding

The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.

IsRNA1: De Novo Prediction and Blind Screening of RNA 3D Structures

Modeling structures and functions of large ribonucleic acid (RNAs) especially with complicated topologies is highly challenging due to the inefficiency of large conformational sampling and the presence of complicated tertiary interactions. To address this problem, one highly promising approach is coarse-grained modeling. Here, following an iterative simulated reference state approach to decipher the correlations between different structural parameters, we developed a potent coarse-grained RNA model named as IsRNA1 for RNA studies. Molecular dynamics simulations in the IsRNA1 can predict the native structures of small RNAs from a sequence and fold medium-sized RNAs into near-native tertiary structures with the assistance of secondary structure constraints. A large-scale benchmark test on RNA 3D structure prediction shows that IsRNA1 exhibits improved performance for relatively large RNAs of complicated topologies, such as large stem-loop structures and structures containing long-range tertiary interactions. The advantages of IsRNA1 include the consideration of the correlations between the different structural variables, the appropriate characterization of canonical base-pairing and base-stacking interactions, and the better sampling for the backbone conformations. Moreover, a blind screening protocol was developed based on IsRNA1 to identify good structural models from a pool of candidates without prior knowledge of the native structures.

A Bayes-inspired theory for optimally building an efficient coarse-grained folding force field

Because of their potential utility in predicting conformational changes and assessing folding dynamics, coarse-grained (CG) RNA folding models are appealing for rapid characterization of RNA molecules. Previously, we reported the iterative simulated RNA reference state (IsRNA) method for parameterizing a CG force field for RNA folding, which consecutively updates the simulation force field to reflect marginal distributions of folding coordinates in the structure database and extract various energy terms. While the IsRNA model was validated by showing close agreement between the IsRNA-simulated and experimentally observed distributions, here, we expand our theoretical understanding of the model and, in doing so, improve the parameterization process to optimize the subset of included folding coordinates, which leads to accelerated simulations. Using statistical mechanical theory, we analyze the underlying, Bayesian concept that drives parameterization of the energy function, providing a general method for developing predictive, knowledge-based, polymer force fields on the basis of limited data. Furthermore, we propose an optimal parameterization procedure, based on the principal of maximum entropy.

Fundamental studies of functional nucleic acids: aptamers, riboswitches, ribozymes and DNAzymes

This review aims at juxtaposing common versus distinct structural and functional strategies that are applied by aptamers, riboswitches, and ribozymes/DNAzymes. Focusing on recently discovered systems, we begin our analysis with small-molecule binding aptamers, with emphasis on in vitro-selected fluorogenic RNA aptamers and their different modes of ligand binding and fluorescence activation. Fundamental insights are much needed to advance RNA imaging probes for detection of exo- and endogenous RNA and for RNA process tracking. Secondly, we discuss the latest gene expression-regulating mRNA riboswitches that respond to the alarmone ppGpp, to PRPP, to NAD⁺, to adenosine and cytidine diphosphates, and to precursors of thiamine biosynthesis (HMP-PP), and we outline new subclasses of SAM and tetrahydrofolate-binding RNA regulators. Many riboswitches bind protein enzyme cofactors that, in principle, can catalyse a chemical reaction. For RNA, however, only one system (glmS ribozyme) has been identified in Nature thus far that utilizes a small molecule - glucosamine-6-phosphate - to participate directly in reaction catalysis (phosphodiester cleavage). We wonder why that is the case and what is to be done to reveal such likely existing cellular activities that could be more diverse than currently imagined. Thirdly, this brings us to the four latest small nucleolytic ribozymes termed twister, twister-sister, pistol, and hatchet as well as to in vitro selected DNA and RNA enzymes that promote new chemistry, mainly by exploiting their ability for RNA labelling and nucleoside modification recognition. Enormous progress in understanding the strategies of nucleic acids catalysts has been made by providing thorough structural fundaments (e.g. first structure of a DNAzyme, structures of ribozyme transition state mimics) in combination with functional assays and atomic mutagenesis.

Crystal structure of an atypical cobalamin riboswitch reveals RNA structural adaptability as basis for promiscuous ligand binding

Cobalamin riboswitches encompass a structurally diverse group of cis-acting, gene regulatory elements found mostly in bacterial messenger RNA and are classified into subtypes based on secondary and tertiary characteristics. An unusual variant of the cobalamin riboswitch with predicted structural features was identified in Bacillus subtilis over a decade ago, but its structure and mechanisms of cobalamin selectivity and translational control have remained unsolved. We present the crystal structure of the aptamer domain of this atypical cobalamin riboswitch and a model for the complete riboswitch, including its expression platform domain. We demonstrate that this riboswitch binds to multiple cobalamin derivatives and correlate its promiscuous behavior to its structure and unique arrangement of peripheral elements. Comparative structural analyses between conventional cobalamin riboswitches and the B. subtilis cobalamin riboswitch reveal that the likely basis for this promiscuous ligand binding is intrinsic structural adaptability encoded in the RNA structure. It suggests that cobalamin selectivity might ultimately be viewed as existing on a spectrum of affinity for each derivative rather than as belonging to distinct types based on ligand specificities. Our work provides an interesting and notable example of functional coupling of ligand-sensing and adaptive folding by a structured RNA molecule.

Tumor Cell-Intrinsic Immunometabolism and Precision Nutrition in Cancer Immunotherapy

One of the greatest challenges in the cancer immunotherapy field is the need to biologically rationalize and broaden the clinical utility of immune checkpoint inhibitors (ICIs). The balance between metabolism and immune response has critical implications for overcoming the major weaknesses of ICIs, including their lack of universality and durability. The last decade has seen tremendous advances in understanding how the immune system's ability to kill tumor cells requires the conspicuous metabolic specialization of T-cells. We have learned that cancer cell-associated metabolic activities trigger shifts in the abundance of some metabolites with immunosuppressory roles in the tumor microenvironment. Yet very little is known about the tumor cell-intrinsic metabolic traits that control the immune checkpoint contexture in cancer cells. Likewise, we lack a comprehensive understanding of how systemic metabolic perturbations in response to dietary interventions can reprogram the immune checkpoint landscape of tumor cells. We here review state-of-the-art molecular- and functional-level interrogation approaches to uncover how cell-autonomous metabolic traits and diet-mediated changes in nutrient availability and utilization might delineate new cancer cell-intrinsic metabolic dependencies of tumor immunogenicity. We propose that clinical monitoring and in-depth molecular evaluation of the cancer cell-intrinsic metabolic traits involved in primary, adaptive, and acquired resistance to cancer immunotherapy can provide the basis for improvements in therapeutic responses to ICIs. Overall, these approaches might guide the use of metabolic therapeutics and dietary approaches as novel strategies to broaden the spectrum of cancer patients and indications that can be effectively treated with ICI-based cancer immunotherapy.

Probing G-quadruplex topologies and recognition concurrently in real time and 3D using a dual-app nucleoside probe

Comprehensive understanding of structure and recognition properties of regulatory nucleic acid elements in real time and atomic level is highly important to devise efficient therapeutic strategies. Here, we report the establishment of an innovative biophysical platform using a dual-app nucleoside analog, which serves as a common probe to detect and correlate different GQ structures and ligand binding under equilibrium conditions and in 3D by fluorescence and X-ray crystallography techniques. The probe (SedU) is composed of a microenvironment-sensitive fluorophore and an excellent anomalous X-ray scatterer (Se), which is assembled by attaching a selenophene ring at 5-position of 2'-deoxyuridine. SedU incorporated into the loop region of human telomeric DNA repeat fluorescently distinguished subtle differences in GQ topologies and enabled quantify ligand binding to different topologies. Importantly, anomalous X-ray dispersion signal from Se could be used to determine the structure of GQs. As the probe is minimally perturbing, a direct comparison of fluorescence data and crystal structures provided structural insights on how the probe senses different GQ conformations without affecting the native fold. Taken together, our dual-app probe represents a new class of tool that opens up new experimental strategies to concurrently investigate nucleic acid structure and recognition in real time and 3D.

Folding intermediate states of the parallel human telomeric G-quadruplex DNA explored using Well-Tempered Metadynamics

An increasingly comprehension of the folding intermediate states of DNA G-quadruplexes (G4s) is currently an important scientific challenge, especially for the human telomeric (h-tel) G4s-forming sequences, characterized by a highly polymorphic nature. Despite the G-triplex conformation was proposed as one of the possible folding intermediates for the antiparallel and hybrid h-tel G4s, for the parallel h-tel topology with an all-anti guanine orientation, a vertical strand-slippage involving the G-triplets was proposed in previous works through microseconds-long standard molecular dynamics simulations (MDs). Here, in order to get further insights into the vertical strand-slippage and the folding intermediate states of the parallel h-tel G4s, we have carried out a Well-Tempered Metadynamics simulation (WT-MetaD), which allowed us to retrieve an ensemble of six G4s having two/G-tetrad conformations derived by the G-triplets vertical slippage. The insights highlighted in this work are aimed at rationalizing the mechanistic characterisation of the parallel h-tel G4 folding process.

Small synthetic molecule-stabilized RNA pseudoknot as an activator for -1 ribosomal frameshifting

Programmed –1 ribosomal frameshifting (−1PRF) is a recoding mechanism to make alternative proteins from a single mRNA transcript. −1PRF is stimulated by cis-acting signals in mRNA, a seven-nucleotide slippery sequence and a downstream secondary structure element, which is often a pseudoknot. In this study we engineered the frameshifting pseudoknot from the mouse mammary tumor virus to respond to a rationally designed small molecule naphthyridine carbamate tetramer (NCTn). We demonstrate that NCTn can stabilize the pseudoknot structure in mRNA and activate –1PRF both in vitro and in human cells. The results illustrate how NCTn-inducible –1PRF may serve as an important component of the synthetic biology toolbox for the precise control of gene expression using small synthetic molecules.

Using Rosetta for RNA homology modeling

The three-dimensional structures of RNA molecules provide rich and often critical information for understanding their functions, including how they recognize small molecule and protein partners. Computational modeling of RNA 3D structure is becoming increasingly accurate, particularly with the availability of growing numbers of template structures already solved experimentally and the development of sequence alignment and 3D modeling tools to take advantage of this database. For several recent "RNA puzzle" blind modeling challenges, we have successfully identified useful template structures and achieved accurate structure predictions through homology modeling tools developed in the Rosetta software suite. We describe our semi-automated methodology here and walk through two illustrative examples: an adenine riboswitch aptamer, modeled from a template guanine riboswitch structure, and a SAM I/IV riboswitch aptamer, modeled from a template SAM I riboswitch structure.

Identification of bioactive metabolites using activity metabolomics

The metabolome, the collection of small-molecule chemical entities involved in metabolism, has traditionally been studied with the aim of identifying biomarkers in the diagnosis and prediction of disease. However, the value of metabolome analysis (metabolomics) has been redefined from a simple biomarker identification tool to a technology for the discovery of active drivers of biological processes. It is now clear that the metabolome affects cellular physiology through modulation of other 'omics' levels, including the genome, epigenome, transcriptome and proteome. In this Review, we focus on recent progress in using metabolomics to understand how the metabolome influences other omics and, by extension, to reveal the active role of metabolites in physiology and disease. This concept of utilizing metabolomics to perform activity screens to identify biologically active metabolites - which we term activity metabolomics - is already having a broad impact on biology.

Bacterial riboswitches and RNA thermometers: Nature and contributions to pathogenesis

Bacterial pathogens are always challenged by fluctuations of chemical and physical parameters that pose serious threats to cellular integrity and metabolic status. Sudden deprivation of nutrients or key metabolites, changes in surrounding pH, and temperature shifts are the most important examples of such parameters. To elicit a proper response to such fluctuations, bacterial cells coordinate the expression of parameter-relevant genes. Although protein-mediated control of gene expression is well appreciated since many decades, RNA-based regulation has been discovered in early 2000s as a parallel level of regulation. Small regulatory RNAs have emerged as one of the most widespread and important gene regulatory systems in bacteria with rare representatives found in Archaea and Eukarya. Riboswitches and thermosensors are cis-encoded RNA regulatory elements that employ different mechanisms to regulate the expression of related genes controlling key metabolic pathways and genes of temperature relevant proteins including virulence factors. The extent of RNA contributions to gene regulation is not completely known even in well-studied models such E. coli and B. subtilis. In depth understanding of riboswitches is promising for opportunity to discover a narrow spectrum antibacterial drugs that target riboswitches of essential metabolic pathways.

Efficient computation of co-transcriptional RNA-ligand interaction dynamics

Riboswitches form an abundant class of cis-regulatory RNA elements that mediate gene expression by binding a small metabolite. For synthetic biology applications, they are becoming cheap and accessible systems for selectively triggering transcription or translation of downstream genes. Many riboswitches are kinetically controlled, hence knowledge of their co-transcriptional mechanisms is essential. We present here an efficient implementation for analyzing co-transcriptional RNA-ligand interaction dynamics. This approach allows for the first time to model concentration-dependent metabolite binding/unbinding kinetics. We exemplify this novel approach by means of the recently studied I-A 2^'-deoxyguanosine (2^'dG)-sensing riboswitch from Mesoplasma florum.

IsRNA: An Iterative Simulated Reference State Approach to Modeling Correlated Interactions in RNA Folding

Coarse-grained RNA folding models promise great potential for RNA structure prediction. A key component in a coarse-grained folding model is the force field. One of the challenges in the coarse-grained force field calculation is how to treat the correlation between the different degrees of freedoms. Here, we describe a new approach (IsRNA) to extract the correlated energy functions from the known structures. Through iterative molecular dynamics simulations, we build the correlation effects into the reference states, from which we extract the energy functions. The validity of IsRNA is supported by the close agreement between the simulated Boltzmann-like probability distributions for all the structure parameters and those observed from the experimentally determined structures. The correlated energy functions derived here may provide a new tool for RNA 3D structure prediction.

ROC'n'Ribo: Characterizing a Riboswitching Expression System by Modeling Single-Cell Data

RNA-engineered systems offer simple and versatile control over gene expression in many organisms. In particular, the design and implementation of riboswitches presents a unique opportunity to manipulate any reporter device in cis, executing tight temporal and spatial control at low metabolic costs. Assembled to higher order genetic circuits, such riboswitch-regulated devices may efficiently process logical operations. Here, we propose a hierarchical stochastic modeling approach to characterize an in silico repressor gate based on neomycin- and tetracycline-sensitive riboswitches. The model was calibrated on rich, transient in vivo single-cell data to account for cell-to-cell variability. To capture the effect of this variability on gate performance we employed the well-known ROC-analysis and derived a novel performance indicator for logic gates. Introduction of such a performance measure is necessary, since we aimed to assess the correct functionality of the gate at the single-cell level-a prerequisite for its further adaption to a genetic circuitry. Our results may be applied to other genetic devices to analyze their efficiency and ensure their correct performance in the light of cell-to-cell variability.

Co-Transcriptional Folding and Regulation Mechanisms of Riboswitches

Riboswitches are genetic control elements within non-coding regions of mRNA. These self-regulatory elements have been found to sense a range of small metabolites, ions, and other physical signals to exert regulatory control of transcription, translation, and splicing. To date, more than a dozen riboswitch classes have been characterized that vary widely in size and secondary structure. Extensive experiments and theoretical studies have made great strides in understanding the general structures, genetic mechanisms, and regulatory activities of individual riboswitches. As the ligand-dependent co-transcriptional folding and unfolding dynamics of riboswitches are the key determinant of gene expression, it is important to investigate the thermodynamics and kinetics of riboswitches both in the presence and absence of metabolites under the transcription. This review will provide a brief summary of the studies about the regulation mechanisms of the pbuE, S_MK, yitJ, and metF riboswitches based on the ligand-dependent co-transcriptional folding of the riboswitches.

Site-Specific Self-Catalyzed DNA Depurination: A Biological Mechanism That Leads to Mutations and Creates Sequence Diversity

Self-catalyzed DNA depurination is a sequence-specific physiological mechanism mediated by spontaneous extrusion of a stem-loop catalytic intermediate. Hydrolysis of the 5'G residue of the 5'GA/TGG loop and of the first 5'A residue of the 5'GAGA loop, together with particular first stem base pairs, specifies their hydrolysis without involving protein, cofactor, or cation. As such, this mechanism is the only known DNA catalytic activity exploited by nature. The consensus sequences for self-depurination of such G- and A-loop residues occur in all genomes examined across the phyla, averaging one site every 2,000-4,000 base pairs. Because apurinic sites are subject to error-prone repair, leading to substitution and short frameshift mutations, they are both a source of genome damage and a means for creating sequence diversity. Their marked overrepresentation in genomes, and largely unchanging density from the lowest to the highest organisms, indicate their selection over the course of evolution. The mutagenicity at such sites in many human genes is associated with loss of function of key proteins responsible for diverse diseases.

Bioinformatic analysis of riboswitch structures uncovers variant classes with altered ligand specificity

Riboswitches are RNAs that form complex, folded structures that selectively bind small molecules or ions. As with certain groups of protein enzymes and receptors, some riboswitch classes have evolved to change their ligand specificity. We developed a procedure to systematically analyze known riboswitch classes to find additional variants that have altered their ligand specificity. This approach uses multiple-sequence alignments, atomic-resolution structural information, and riboswitch gene associations. Among the discoveries are unique variants of the guanine riboswitch class that most tightly bind the nucleoside 2'-deoxyguanosine. In addition, we identified variants of the glycine riboswitch class that no longer recognize this amino acid, additional members of a rare flavin mononucleotide (FMN) variant class, and also variants of c-di-GMP-I and -II riboswitches that might recognize different bacterial signaling molecules. These findings further reveal the diverse molecular sensing capabilities of RNA, which highlights the potential for discovering a large number of additional natural riboswitch classes.

Structure of the Ribosomal RNA Decoding Site Containing a Selenium-Modified Responsive Fluorescent Ribonucleoside Probe

Comprehensive understanding of the structure–function relationship of RNA both in real time and at atomic level will have a profound impact in advancing our understanding of RNA functions in biology. Here, we describe the first example of a multifunctional nucleoside probe, containing a conformation‐sensitive fluorophore and an anomalous X‐ray diffraction label (5‐selenophene uracil), which enables the correlation of RNA conformation and recognition under equilibrium and in 3D. The probe incorporated into the bacterial ribosomal RNA decoding site, fluorescently reports antibiotic binding and provides diffraction information in determining the structure without distorting native RNA fold. Further, by comparing solution binding data and crystal structure, we gained insight on how the probe senses ligand‐induced conformational change in RNA. Taken together, our nucleoside probe represents a new class of biophysical tool that would complement available tools for functional RNA investigations.

Single-molecule analysis reveals multi-state folding of a guanine riboswitch

Guanine-responsive riboswitches undergo ligand-dependent structural rearrangements to control gene expression by transcription termination. While the molecular basis for ligand recognition is well established, the associated structural rearrangements and the kinetics involved in the formation of the aptamer domain are less well understood. Using high-resolution optical tweezers, we followed the folding trajectories of a single molecule of the xpt-pbuX guanine aptamer from Bacillus subtilis. We report a rapid six-state conformational rearrangement, in which three of the states are guanine dependent, during the transition from the linear to the native receptor conformation. The folding completes in <1 s. The force-dependent equilibrium kinetics and the mutational data indicated that the flexible J2-J3 junction undergoes a ligand-dependent conformational switching, which triggers the formation of the long-range tertiary interactions and the P1 helix. In the absence of the right ligand, the junction failed to initiate the series of conformational rearrangements required for the riboswitch activities.

Multianalytical Study of the Binding between a Small Chiral Molecule and a DNA Aptamer: Evidence for Asymmetric Steric Effect upon 3'- versus 5'-End Sequence Modification

Nucleic acid aptamers are involved in a broad field of applications ranging from therapeutics to analytics. Deciphering the binding mechanisms between aptamers and small ligands is therefore crucial to improve and optimize existing applications and to develop new ones. Particularly interesting is the enantiospecific binding mechanism involving small molecules with nonprestructured aptamers. One archetypal example is the chiral binding between l-tyrosinamide and its 49-mer aptamer for which neither structural nor mechanistic information is available. In the present work, we have taken advantage of a multiple analytical characterization strategy (i.e., using electroanalytical techniques such as kinetic rotating droplet electrochemistry, fluorescence polarization, isothermal titration calorimetry, and quartz crystal microbalance) for interpreting the nature of binding process. Screening of the binding thermodynamics and kinetics with a wide range of aptamer sequences revealed the lack of symmetry between the two ends of the 23-mer minimal binding sequence, showing an unprecedented influence of the 5' aptamer modification on the bimolecular binding rate constant k_on and no significant effect on the dissociation rate constant k_off. The results we have obtained lead us to conclude that the enantiospecific binding reaction occurs through an induced-fit mechanism, wherein the ligand promotes a primary nucleation binding step near the 5'-end of the aptamer followed by a directional folding of the aptamer around its target from 5'-end to 3'-end. Functionalization of the 5'-end position by a chemical label, a polydA tail, a protein, or a surface influences the kinetic/thermodynamic constants up to 2 orders of magnitude in the extreme case of a surface immobilized aptamer, while significantly weaker effect is observed for a 3'-end modification. The reason is that steric hindrance must be overcome to nucleate the binding complex in the presence of a modification near the nucleation site.

Tuning RNA Flexibility with Helix Length and Junction Sequence

The increasing awareness of RNA's central role in biology calls for a new understanding of how RNAs, like proteins, recognize biological partners. Because RNA is inherently flexible, it assumes a variety of conformations. This conformational flexibility can be a critical aspect of how RNA attracts and binds molecular partners. Structurally, RNA consists of rigid basepaired duplexes, separated by flexible non-basepaired regions. Here, using an RNA system consisting of two short helices, connected by a single-stranded (non-basepaired) junction, we explore the role of helix length and junction sequence in determining the range of conformations available to a model RNA. Single-molecule Förster resonance energy transfer reports on the RNA conformation as a function of either mono- or divalent ion concentration. Electrostatic repulsion between helices dominates at low salt concentration, whereas junction sequence effects determine the conformations at high salt concentration. Near physiological salt concentrations, RNA conformation is sensitive to both helix length and junction sequence, suggesting a means for sensitively tuning RNA conformations.

Fluorescence-Based Strategies to Investigate the Structure and Dynamics of Aptamer-Ligand Complexes

In addition to the helical nature of double-stranded DNA and RNA, single-stranded oligonucleotides can arrange themselves into tridimensional structures containing loops, bulges, internal hairpins and many other motifs. This ability has been used for more than two decades to generate oligonucleotide sequences, so-called aptamers, that can recognize certain metabolites with high affinity and specificity. More recently, this library of artificially-generated nucleic acid aptamers has been expanded by the discovery that naturally occurring RNA sequences control bacterial gene expression in response to cellular concentration of a given metabolite. The application of fluorescence methods has been pivotal to characterize in detail the structure and dynamics of these aptamer-ligand complexes in solution. This is mostly due to the intrinsic high sensitivity of fluorescence methods and also to significant improvements in solid-phase synthesis, post-synthetic labeling strategies and optical instrumentation that took place during the last decade. In this work, we provide an overview of the most widely employed fluorescence methods to investigate aptamer structure and function by describing the use of aptamers labeled with a single dye in fluorescence quenching and anisotropy assays. The use of 2-aminopurine as a fluorescent analog of adenine to monitor local changes in structure and fluorescence resonance energy transfer (FRET) to follow long-range conformational changes is also covered in detail. The last part of the review is dedicated to the application of fluorescence techniques based on single-molecule microscopy, a technique that has revolutionized our understanding of nucleic acid structure and dynamics. We finally describe the advantages of monitoring ligand-binding and conformational changes, one molecule at a time, to decipher the complexity of regulatory aptamers and summarize the emerging folding and ligand-binding models arising from the application of these single-molecule FRET microscopy techniques.

Coarse-grained modeling of RNA 3D structure

Functional RNA molecules depend on three-dimensional (3D) structures to carry out their tasks within the cell. Understanding how these molecules interact to carry out their biological roles requires a detailed knowledge of RNA 3D structure and dynamics as well as thermodynamics, which strongly governs the folding of RNA and RNA-RNA interactions as well as a host of other interactions within the cellular environment. Experimental determination of these properties is difficult, and various computational methods have been developed to model the folding of RNA 3D structures and their interactions with other molecules. However, computational methods also have their limitations, especially when the biological effects demand computation of the dynamics beyond a few hundred nanoseconds. For the researcher confronted with such challenges, a more amenable approach is to resort to coarse-grained modeling to reduce the number of data points and computational demand to a more tractable size, while sacrificing as little critical information as possible. This review presents an introduction to the topic of coarse-grained modeling of RNA 3D structures and dynamics, covering both high- and low-resolution strategies. We discuss how physics-based approaches compare with knowledge based methods that rely on databases of information. In the course of this review, we discuss important aspects in the reasoning process behind building different models and the goals and pitfalls that can result.

3D RNA and Functional Interactions from Evolutionary Couplings

Non-coding RNAs are ubiquitous, but the discovery of new RNA gene sequences far outpaces the research on the structure and functional interactions of these RNA gene sequences. We mine the evolutionary sequence record to derive precise information about the function and structure of RNAs and RNA-protein complexes. As in protein structure prediction, we use maximum entropy global probability models of sequence co-variation to infer evolutionarily constrained nucleotide-nucleotide interactions within RNA molecules and nucleotide-amino acid interactions in RNA-protein complexes. The predicted contacts allow all-atom blinded 3D structure prediction at good accuracy for several known RNA structures and RNA-protein complexes. For unknown structures, we predict contacts in 160 non-coding RNA families. Beyond 3D structure prediction, evolutionary couplings help identify important functional interactions-e.g., at switch points in riboswitches and at a complex nucleation site in HIV. Aided by increasing sequence accumulation, evolutionary coupling analysis can accelerate the discovery of functional interactions and 3D structures involving RNA.

Conformational Rearrangements of Individual Nucleotides during RNA-Ligand Binding Are Rate-Differentiated

A pronounced rate differentiation has been found for conformational rearrangements of individual nucleobases that occur during ligand recognition of the preQ1 class-I riboswitch aptamer from Thermoanaerobacter tengcongensis. Rate measurements rely on the 2ApFold approach by analyzing the fluorescence response of riboswitch variants, each with a single, strategically positioned 2-aminopurine nucleobase substitution. Observed rate discrimination between the fastest and the slowest conformational adaption is 22-fold, with the largest rate observed for the rearrangement of a nucleoside directly at the binding site and the smallest rate observed for the 3'-unpaired nucleoside that stacks onto the pseudo-knot-closing Watson-Crick base pair. Our findings provide novel insights into how compact, prefolded RNAs that follow the induced-fit recognition mechanism adapt local structural elements in response to ligand binding on a rather broad time scale and how this process culminates in a structural signal that is responsible for efficient downregulation of ribosomal translation.

Quantitative and predictive model of kinetic regulation by E. coli TPP riboswitches

Riboswitches are non-coding elements upstream or downstream of mRNAs that, upon binding of a specific ligand, regulate transcription and/or translation initiation in bacteria, or alternative splicing in plants and fungi. We have studied thiamine pyrophosphate (TPP) riboswitches regulating translation of thiM operon and transcription and translation of thiC operon in E. coli, and that of THIC in the plant A. thaliana. For all, we ascertained an induced-fit mechanism involving initial binding of the TPP followed by a conformational change leading to a higher-affinity complex. The experimental values obtained for all kinetic and thermodynamic parameters of TPP binding imply that the regulation by A. thaliana riboswitch is governed by mass-action law, whereas it is of kinetic nature for the two bacterial riboswitches. Kinetic regulation requires that the RNA polymerase pauses after synthesis of each riboswitch aptamer to leave time for TPP binding, but only when its concentration is sufficient. A quantitative model of regulation highlighted how the pausing time has to be linked to the kinetic rates of initial TPP binding to obtain an ON/OFF switch in the correct concentration range of TPP. We verified the existence of these pauses and the model prediction on their duration. Our analysis also led to quantitative estimates of the respective efficiency of kinetic and thermodynamic regulations, which shows that kinetically regulated riboswitches react more sharply to concentration variation of their ligand than thermodynamically regulated riboswitches. This rationalizes the interest of kinetic regulation and confirms empirical observations that were obtained by numerical simulations.

SimRNA: a coarse-grained method for RNA folding simulations and 3D structure prediction

RNA molecules play fundamental roles in cellular processes. Their function and interactions with other biomolecules are dependent on the ability to form complex three-dimensional (3D) structures. However, experimental determination of RNA 3D structures is laborious and challenging, and therefore, the majority of known RNAs remain structurally uncharacterized. Here, we present SimRNA: a new method for computational RNA 3D structure prediction, which uses a coarse-grained representation, relies on the Monte Carlo method for sampling the conformational space, and employs a statistical potential to approximate the energy and identify conformations that correspond to biologically relevant structures. SimRNA can fold RNA molecules using only sequence information, and, on established test sequences, it recapitulates secondary structure with high accuracy, including correct prediction of pseudoknots. For modeling of complex 3D structures, it can use additional restraints, derived from experimental or computational analyses, including information about secondary structure and/or long-range contacts. SimRNA also can be used to analyze conformational landscapes and identify potential alternative structures.

G-quadruplexes with (4n - 1) guanines in the G-tetrad core: formation of a G-triad·water complex and implication for small-molecule binding

G-quadruplexes are non-canonical structures of nucleic acids, in which guanine bases form planar G-tetrads (G·G·G·G) that stack on each other in the core of the structure. G-quadruplexes generally contain multiple times of four (4n) guanines in the core. Here, we study the structure of G-quadruplexes with only (4n - 1) guanines in the core. The solution structure of a DNA sequence containing 11 guanines showed the formation of a parallel G-quadruplex involving two G-tetrads and one G-triad with a vacant site. Molecular dynamics simulation established the formation of a stable G-triad·water complex, where water molecules mimic the position of the missing guanine in the vacant site. The concept of forming G-quadruplexes with missing guanines in the core broadens the current definition of G-quadruplex-forming sequences. The potential ability of such structures to bind different metabolites, including guanine, guanosine and GTP, in the vacant site, could have biological implications in regulatory functions. Formation of this unique binding pocket in the G-triad could be used as a specific target in drug design.

Recognition of the bacterial alarmone ZMP through long-distance association of two RNA subdomains

The bacterial alarmone 5-aminoimidazole-4-carboxamide riboside 5'-triphosphate (ZTP), derived from the monophosphorylated purine precursor ZMP, accumulates during folate starvation. ZTP regulates genes involved in purine and folate metabolism through a cognate riboswitch. The linker connecting this riboswitch’s two sub-domains varies in length by over 100 nucleotides. We report the co-crystal structure of the Fusobacterium ulcerans riboswitch bound to ZMP, which spans the two sub-domains whose interface also comprises a pseudoknot and ribose zipper. The riboswitch recognizes the carboxamide oxygen of ZMP through an unprecedented inner-sphere coordination with a Mg²⁺ ion. We demonstrate that the affinity of the riboswitch for ZMP is modulated by the linker length. Notably, ZMP can bind to the two sub-domains together even when synthesized as separate RNAs. The ZTP riboswitch demonstrates how specific small-molecule binding can drive association of distant non-coding RNA domains to regulate gene expression.

Regulation of bacterial virulence by Csr (Rsm) systems

Most bacterial pathogens have the remarkable ability to flourish in the external environment and in specialized host niches. This ability requires their metabolism, physiology, and virulence factors to be responsive to changes in their surroundings. It is no surprise that the underlying genetic circuitry that supports this adaptability is multilayered and exceedingly complex. Studies over the past 2 decades have established that the CsrA/RsmA proteins, global regulators of posttranscriptional gene expression, play important roles in the expression of virulence factors of numerous proteobacterial pathogens. To accomplish these tasks, CsrA binds to the 5' untranslated and/or early coding regions of mRNAs and alters translation, mRNA turnover, and/or transcript elongation. CsrA activity is regulated by noncoding small RNAs (sRNAs) that contain multiple CsrA binding sites, which permit them to sequester multiple CsrA homodimers away from mRNA targets. Environmental cues sensed by two-component signal transduction systems and other regulatory factors govern the expression of the CsrA-binding sRNAs and, ultimately, the effects of CsrA on secretion systems, surface molecules and biofilm formation, quorum sensing, motility, pigmentation, siderophore production, and phagocytic avoidance. This review presents the workings of the Csr system, the paradigm shift that it generated for understanding posttranscriptional regulation, and its roles in virulence networks of animal and plant pathogens.

Probing RNA-protein networks: biochemistry meets genomics

RNA-protein interactions are pervasive. The specificity of these interactions dictates which RNAs are controlled by what protein. Here we describe a class of revolutionary new methods that enable global views of RNA-binding specificity in vitro, for both single proteins and multiprotein complexes. These methods provide insight into central issues in RNA regulation in living cells, including understanding the balance between free and bound components, the basis for exclusion of binding sites, detection of binding events in the absence of discernible regulatory elements, and new approaches to targeting endogenous transcripts by design. Comparisons of in vitro and in vivo binding provide a foundation for comprehensive understanding of the biochemistry of protein-mediated RNA regulatory networks.

Heavy atom containing fluorescent ribonucleoside analog probe for the fluorescence detection of RNA-ligand binding

Although numerous biophysical tools have provided effective systems to study nucleic acids, our current knowledge on how RNA structure complements its function is limited. Therefore, development of robust tools to study the structure–function relationship of RNA is highly desired. Toward this endeavor, we have developed a new ribonucleoside analog, based on a (selenophen-2-yl)pyrimidine core, which could serve as a fluorescence probe to study the function of RNA in real time and as an anomalous scattering label (selenium atom) for the phase determination in X-ray crystallography. The fluorescent selenophene-modified uridine analog is minimally perturbing and exhibits probe-like properties such as sensitivity to microenvironment and conformation changes. Utilizing these properties and amicability of the corresponding ribonucleotide analog to enzymatic incorporation, we have synthesized a fluorescent bacterial ribosomal decoding site (A-site) RNA construct and have developed a fluorescence binding assay to effectively monitor the binding of aminoglycoside antibiotics to the A-site. Our results demonstrate that this simple approach of building a dual probe could provide new avenues to study the structure–function relationship of not only nucleic acids, but also other biomacromolecules.

Site-specific fluorescence dynamics in an RNA 'thermometer' reveals the role of ribosome binding in its temperature-sensitive switch function

RNA thermometers control the translation of several heat shock and virulence genes by their temperature-sensitive structural transitions. Changes in the structure and dynamics of MiniROSE RNA, which regulates translation in the temperature range of 20–45°C, were studied by site specifically replacing seven adenine residues with the fluorescent analog, 2-aminopurine (2-AP), one at a time. Dynamic fluorescence observables of 2-AP-labeled RNAs were compared in their free versus ribosome-bound states for the first time. Noticeably, position dependence of fluorescence observables, which was prominent at 20°C, was persistent even at 45ºC, suggesting the persistence of structural integrity up to 45ºC. Interestingly, position-dependent dispersion of fluorescence lifetime and quenching constant at 45°C was ablated in ribosome-bound state, when compared to those at 20°C, underscoring loss of structural integrity at 45°C, in ribosome-bound RNA. Significant increase in the value of mean lifetime for 2-AP corresponding to Shine–Dalgarno sequences, when the temperature was raised from 20 to 45°C, to values seen in the presence of urea at 45°C was a strong indicator of melting of the 3D structure of MiniROSE RNA at 45°C, only when it was ribosome bound. Taken all together, we propose a model where we invoke that ribosome binding of the RNA thermometer critically regulates temperature sensing functions in MiniROSE RNA.

Molecular dynamics simulations identify time scale of conformational changes responsible for conformational selection in molecular recognition of HIV-1 transactivation responsive RNA

The HIV-1 Tat protein and several small molecules bind to HIV-1 transactivation responsive RNA (TAR) by selecting sparsely populated but pre-existing conformations. Thus, a complete characterization of TAR conformational ensemble and dynamics is crucial to understand this paradigmatic system and could facilitate the discovery of new antivirals targeting this essential regulatory element. We show here that molecular dynamics simulations can be effectively used toward this goal by bridging the gap between functionally relevant time scales that are inaccessible to current experimental techniques. Specifically, we have performed several independent microsecond long molecular simulations of TAR based on one of the most advanced force fields available for RNA, the parmbsc0 AMBER. Our simulations are first validated against available experimental data, yielding an excellent agreement with measured residual dipolar couplings and order parameter S(2). This contrast with previous molecular dynamics simulations (Salmon et al., J. Am. Chem. Soc. 2013 135, 5457-5466) based on the CHARMM36 force field, which could achieve only modest accord with the experimental RDC values. Next, we direct the computation toward characterizing the internal dynamics of TAR over the microsecond time scale. We show that the conformational fluctuations observed over this previously elusive time scale have a strong functionally oriented character in that they are primed to sustain and assist ligand binding.