Elucidating molecular motion through structural and dynamic filters of energy-minimized conformer ensembles

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

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

Ultraslow Domain Motions in HIV-1 TAR RNA Revealed by Solid-State Deuterium NMR

Intrinsic motions may allow HIV-1 transactivation response (TAR) RNA to change its conformation to form a functional complex with the Tat protein, which is essential for viral replication. Understanding the dynamic properties of TAR necessitates determining motion on the intermediate nanosecond-to-microsecond time scale. To this end, we performed solid-state deuterium NMR line-shape and T_1Z relaxation-time experiments to measure intermediate motions for two uridine residues, U40 and U42, within the lower helix of TAR. We infer global motions at rates of ∼10⁵ s^-1 in the lower helix, which are much slower than those in the upper helix (∼10⁶ s^-1), indicating that the two helical domains reorient independently of one another in the solid-state sample. These results contribute to the aim of fully describing the properties of functional motions in TAR RNA.

RNA folding pathways in stop motion

We introduce a method for predicting RNA folding pathways, with an application to the most important RNA tetraloops. The method is based on the idea that ensembles of three-dimensional fragments extracted from high-resolution crystal structures are heterogeneous enough to describe metastable as well as intermediate states. These ensembles are first validated by performing a quantitative comparison against available solution nuclear magnetic resonance (NMR) data of a set of RNA tetranucleotides. Notably, the agreement is better with respect to the one obtained by comparing NMR with extensive all-atom molecular dynamics simulations. We then propose a procedure based on diffusion maps and Markov models that makes it possible to obtain reaction pathways and their relative probabilities from fragment ensembles. This approach is applied to study the helix-to-loop folding pathway of all the tetraloops from the GNRA and UNCG families. The results give detailed insights into the folding mechanism that are compatible with available experimental data and clarify the role of intermediate states observed in previous simulation studies. The method is computationally inexpensive and can be used to study arbitrary conformational transitions.

Integrative, dynamic structural biology at atomic resolution--it's about time

Biomolecules adopt a dynamic ensemble of conformations, each with the potential to interact with binding partners or perform the chemical reactions required for a multitude of cellular functions. Recent advances in X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy and other techniques are helping us realize the dream of seeing--in atomic detail--how different parts of biomolecules shift between functional substates using concerted motions. Integrative structural biology has advanced our understanding of the formation of large macromolecular complexes and how their components interact in assemblies by leveraging data from many low-resolution methods. Here, we review the growing opportunities for integrative, dynamic structural biology at the atomic scale, contending there is increasing synergistic potential between X-ray crystallography, NMR and computer simulations to reveal a structural basis for protein conformational dynamics at high resolution.

NMR studies of dynamic biomolecular conformational ensembles

Multidimensional heteronuclear NMR approaches can provide nearly complete sequential signal assignments of isotopically enriched biomolecules. The availability of assignments together with measurements of spin relaxation rates, residual spin interactions, J-couplings and chemical shifts provides information at atomic resolution about internal dynamics on timescales ranging from ps to ms, both in solution and in the solid state. However, due to the complexity of biomolecules, it is not possible to extract a unique atomic-resolution description of biomolecular motions even from extensive NMR data when many conformations are sampled on multiple timescales. For this reason, powerful computational approaches are increasingly applied to large NMR data sets to elucidate conformational ensembles sampled by biomolecules. In the past decade, considerable attention has been directed at an important class of biomolecules that function by binding to a wide variety of target molecules. Questions of current interest are: "Does the free biomolecule sample a conformational ensemble that encompasses the conformations found when it binds to various targets; and if so, on what time scale is the ensemble sampled?" This article reviews recent efforts to answer these questions, with a focus on comparing ensembles obtained for the same biomolecules by different investigators. A detailed comparison of results obtained is provided for three biomolecules: ubiquitin, calmodulin and the HIV-1 trans-activation response RNA.

Characterizing RNA ensembles from NMR data with kinematic models

Functional mechanisms of biomolecules often manifest themselves precisely in transient conformational substates. Researchers have long sought to structurally characterize dynamic processes in non-coding RNA, combining experimental data with computer algorithms. However, adequate exploration of conformational space for these highly dynamic molecules, starting from static crystal structures, remains challenging. Here, we report a new conformational sampling procedure, KGSrna, which can efficiently probe the native ensemble of RNA molecules in solution. We found that KGSrna ensembles accurately represent the conformational landscapes of 3D RNA encoded by NMR proton chemical shifts. KGSrna resolves motionally averaged NMR data into structural contributions; when coupled with residual dipolar coupling data, a KGSrna ensemble revealed a previously uncharacterized transient excited state of the HIV-1 trans-activation response element stem-loop. Ensemble-based interpretations of averaged data can aid in formulating and testing dynamic, motion-based hypotheses of functional mechanisms in RNAs with broad implications for RNA engineering and therapeutic intervention.