NMR studies of nucleic acid dynamics

Molecular insights into the interaction between a disordered protein and a folded RNA

Intrinsically disordered protein regions (IDRs) are well established as contributors to intermolecular interactions and the formation of biomolecular condensates. In particular, RNA-binding proteins (RBPs) often harbor IDRs in addition to folded RNA-binding domains that contribute to RBP function. To understand the dynamic interactions of an IDR-RNA complex, we characterized the RNA-binding features of a small (68 residues), positively charged IDR-containing protein, Small ERDK-Rich Factor (SERF). At high concentrations, SERF and RNA undergo charge-driven associative phase separation to form a protein- and RNA-rich dense phase. A key advantage of this model system is that this threshold for demixing is sufficiently high that we could use solution-state biophysical methods to interrogate the stoichiometric complexes of SERF with RNA in the one-phase regime. Herein, we describe our comprehensive characterization of SERF alone and in complex with a small fragment of the HIV-1 Trans-Activation Response (TAR) RNA with complementary biophysical methods and molecular simulations. We find that this binding event is not accompanied by the acquisition of structure by either molecule; however, we see evidence for a modest global compaction of the SERF ensemble when bound to RNA. This behavior likely reflects attenuated charge repulsion within SERF via binding to the polyanionic RNA and provides a rationale for the higher-order assembly of SERF in the context of RNA. We envision that the SERF-RNA system will lower the barrier to accessing the details that support IDR-RNA interactions and likewise deepen our understanding of the role of IDR-RNA contacts in complex formation and liquid-liquid phase separation.

Molecular insights into the interaction between a disordered protein and a folded RNA

Intrinsically disordered protein regions (IDRs) are well-established as contributors to intermolecular interactions and the formation of biomolecular condensates. In particular, RNA-binding proteins (RBPs) often harbor IDRs in addition to folded RNA-binding domains that contribute to RBP function. To understand the dynamic interactions of an IDR-RNA complex, we characterized the RNA-binding features of a small (68 residues), positively charged IDR-containing protein, SERF. At high concentrations, SERF and RNA undergo charge-driven associative phase separation to form a protein- and RNA-rich dense phase. A key advantage of this model system is that this threshold for demixing is sufficiently high that we could use solution-state biophysical methods to interrogate the stoichiometric complexes of SERF with RNA in the one-phase regime. Herein, we describe our comprehensive characterization of SERF alone and in complex with a small fragment of the HIV-1 TAR RNA (TAR) with complementary biophysical methods and molecular simulations. We find that this binding event is not accompanied by the acquisition of structure by either molecule; however, we see evidence for a modest global compaction of the SERF ensemble when bound to RNA. This behavior likely reflects attenuated charge repulsion within SERF via binding to the polyanionic RNA and provides a rationale for the higher-order assembly of SERF in the context of RNA. We envision that the SERF-RNA system will lower the barrier to accessing the details that support IDR-RNA interactions and likewise deepen our understanding of the role of IDR-RNA contacts in complex formation and liquid-liquid phase separation.

Structural and computational studies of HIV-1 RNA

Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.

A beginner's handbook to identify and characterize i-motif DNA

Genomic DNA exhibits an innate ability to manifest diverse sequence-dependent secondary structures, serving crucial functions in gene regulation and cellular equilibrium. While extensive research has confirmed the formation of G-quadruplex structures by guanine-rich sequences in vitro and in cells, recent investigations have turned the quadruplex community's attention to the cytosine (C)-rich complementary strands that can adopt unique tetra-stranded conformation, termed as intercalated motif or i-motif. I-motifs are stabilized by hemi-protonated C:CH+ base pairs under acidic conditions. Initially, the in vivo occurrence of i-motifs was underestimated because their formation is favored at non-physiological pH. However, groundbreaking research utilizing the structure-specific iMab antibody and high-throughput sequencing have recently detected their conserved dispersion throughout the genome, challenging previous assumptions. Given the evolving nature of this research field, it becomes imperative to conduct independent in vitro experiments aimed at identifying potential i-motif formation in C-rich sequences and consolidating the findings to address the properties of i-motifs. This chapter serves as an introductory guide for the swift identification of novel i-motifs, where we present an experimental framework for investigating and characterizing i-motif sequences in vitro. In this chapter, we selected a synthetic oligonucleotide (C7T3) sequence and outlined appropriate methodologies for annealing the i-motif structure into suitable buffers. Then, we validated its formation by CD (Circular Dichroism) and NMR (Nuclear Magnetic Resonance) spectroscopy. Finally, we provided a thorough account of the step-by-step procedures to investigate the effect of i-motif formation on the stalling or retardation of DNA replication using high resolution primer extension assays.

Conformational Plasticity of Parallel G-Quadruplex─Implications on Duplex-Quadruplex Motifs

DNA G-quadruplexes are essential motifs in molecular biology performing a wide range of functions enabled by their unique and diverse structures. In this study, we focus on the conformational plasticity of the most abundant and biologically relevant parallel G-quadruplex topology. A multipronged approach of structure survey, solution-state NMR spectroscopy, and molecular dynamics simulations unravels subtle yet essential features of the parallel G-quadruplex topology. Stark differences in flexibility are observed for the nucleotides depending upon their positioning in the tetrad planes that are intricately correlated with the conformational sampling of the propeller loop. Importantly, the terminal nucleotides in the 5'-end versus the 3'-end of the parallel quadruplex display differential dynamics that manifests their ability to accommodate a duplex on either end of the G-quadruplex. The conformational plasticity characterized in this study provides essential cues toward biomolecular processes such as small molecular binding, intermolecular quadruplex stacking, and implications on how a duplex influences the structure of a neighboring quadruplex.

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.

Transient Hoogsteen Base Pairs Observed in Unbiased Molecular Dynamics Simulations of DNA

Duplex DNA is modeled as canonical B-DNA displaying the characteristic Watson-Crick base pairs. A less common and short-lived pairing of the nucleobases is the Hoogsteen (HG) conformation. The low population of the HG base pairs (<1%) necessitates extended sampling times in order to analyze through unbiased molecular dynamics (MD) simulations. We have discovered that with extended sampling times using multiple independent copies of an 18-mer sequence, the MD trajectories reproduce the expected and transient HG base pairing. Consistent with experiment, the percentage of the HG events are within the range of ∼0.1-1.0% over the combined aggregate sampling time of more than 3.6 ms. We present the reliability of the current AMBER set of nucleic acid force fields and tools to accurately simulate naturally occurring base-pairing and opening events without any bias or restraints. The mechanism consists of base pair fraying, flipping of the purine, and reformation with HG base pairs.

Use of anion-exchange HPLC to study DNA conformational polymorphism

Anion-exchange chromatography carried out under non-denaturing conditions is a versatile tool to differentiate DNA conformations. In this work, the utility of this form of HPLC was demonstrated in four examples. The hairpin and duplex forms of d(CG)₉ were readily resolved, which allowed for the studies of the influence of salt on the equilibrium of these two forms of secondary structures. Similarly, the minimum size of Tn in the loop region required for the sequence 5'-d(CCCAA-(T)n-TTGGG)-3' to form hairpin was established to be two nucleotides using anion-exchange HPLC and fluorescence resonance energy transfer. Furthermore, the efficiency of hybridization of partially self-complementary sequences d[(CG)₆N_x] was readily monitored by non-denaturing anion-exchange HPLC. Finally, different structures adopted by quadruplex-forming sequences were resolved in the same manner.

Role of conformational heterogeneity in ligand recognition by viral RNA molecules

Ribonucleic acid (RNA) molecules are known to undergo conformational changes in response to various environmental stimuli including temperature, pH, and ligands. In particular, viral RNA molecules are a key example of conformationally adapting molecules that have evolved to switch between many functional conformations. The transactivation response element (TAR) RNA from the type-1 human immunodeficiency virus (HIV-1) is a viral RNA molecule that is being increasingly explored as a potential therapeutic target due to its role in the viral replication process. In this work, we have studied the dynamics in TAR RNA in apo and liganded states by performing explicit-solvent molecular dynamics (MD) simulations initiated with 27 distinct structures. We determined that the TAR RNA structure is significantly stabilized on ligand binding with especially decreased fluctuations in its two helices. This rigidity is further coupled with the decreased flipping of bulge nucleotides, which were observed to flip more frequently in the absence of ligands. We found that initially-distinct structures of TAR RNA converged to similar conformations on removing ligands. We also report that conformational dynamics in unliganded TAR structures leads to the formation of binding pockets capable of accommodating ligands of various sizes.

Challenges and perspectives for structural biology of lncRNAs-the example of the Xist lncRNA A-repeats

Following the discovery of numerous long non-coding RNA (lncRNA) transcripts in the human genome, their important roles in biology and human disease are emerging. Recent progress in experimental methods has enabled the identification of structural features of lncRNAs. However, determining high-resolution structures is challenging as lncRNAs are expected to be dynamic and adopt multiple conformations, which may be modulated by interaction with protein binding partners. The X-inactive specific transcript (Xist) is necessary for X inactivation during dosage compensation in female placental mammals and one of the best-studied lncRNAs. Recent progress has provided new insights into the domain organization, molecular features, and RNA binding proteins that interact with distinct regions of Xist. The A-repeats located at the 5′ end of the transcript are of particular interest as they are essential for mediating silencing of the inactive X chromosome. Here, we discuss recent progress with elucidating structural features of the Xist lncRNA, focusing on the A-repeats. We discuss the experimental and computational approaches employed that have led to distinct structural models, likely reflecting the intrinsic dynamics of this RNA. The presence of multiple dynamic conformations may also play an important role in the formation of the associated RNPs, thus influencing the molecular mechanism underlying the biological function of the Xist A-repeats. We propose that integrative approaches that combine biochemical experiments and high-resolution structural biology in vitro with chemical probing and functional studies in vivo are required to unravel the molecular mechanisms of lncRNAs.

Ligand Recognition in Viral RNA Necessitates Rare Conformational Transitions

Ribonucleic acids (RNAs) are conformationally flexible molecules that fold into three-dimensional structures and play an important role in different cellular processes as well as in the development of many diseases. RNA has therefore become an important target for developing novel therapeutic approaches. The biophysical processes underlying RNA function are often associated with rare structural transitions that play a key role in ligand recognition. In this work, we probe these rarely occurring transitions using nonequilibrium simulations by characterizing the dissociation of a ligand molecule from an HIV-1 viral RNA element. Specifically, we observed base-flipping rare events that are coupled with ligand binding/unbinding and also provided mechanistic details underlying these transitions.

PyShifts: A PyMOL Plugin for Chemical Shift-Based Analysis of Biomolecular Ensembles

Here, we present PyShifts-a PyMOL plugin for chemical shift-based analysis of biomolecular ensembles. With PyShifts, users can compare and visualize differences between experimentally measured and computationally predicted chemical shifts. When analyzing multiple conformations of a biomolecule with PyShifts, users can also sort a set of conformations based on chemical shift differences and identify the conformers that exhibit the best agreement between measured and predicted chemical shifts. Although we have integrated PyShifts with the chemical shift predictors LARMOR^D and LARMOR^Cα, PyShifts can read in chemical shifts from any source, and so, users can employ PyShifts to analyze biomolecular structures using chemical shifts computed by any chemical shift predictor. We envision, therefore, that PyShifts (https://github.com/atfrank/PyShifts) will find utility as a general-purpose tool for exploring chemical shift-structure relationships in biomolecular ensembles.

Helix-Based RNA Landscape Partition and Alternative Secondary Structure Determination

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

Development of Vibrational Frequency Maps for Nucleobases

Vibrational spectroscopy provides a powerful tool to probe the structure and dynamics of nucleic acids because specific normal modes, particularly the base carbonyl stretch modes, are highly sensitive to the hydrogen bonding patterns and stacking configurations in these biomolecules. In this work, we develop vibrational frequency maps for the C═O and C═C stretches in nucleobases that allow the calculations of their site frequencies directly from molecular dynamics simulations. We assess the frequency maps by applying them to nucleobase derivatives in aqueous solutions and nucleosides in organic solvents and demonstrate that the predicted infrared spectra are in good agreement with experimental measurements. The frequency maps can be readily used to model the linear and nonlinear vibrational spectroscopy of nucleic acids and elucidate the molecular origin of the experimentally observed spectral features.

Relating Structure and Dynamics in RNA Biology

Recent advances in structural biology methods have enabled a surge in the number of RNA and RNA-protein assembly structures available at atomic or near-atomic resolution. These complexes are often trapped in discrete conformational states that exist along a mechanistic pathway. Single-molecule fluorescence methods provide temporal resolution to elucidate the dynamic mechanisms of processes involving complex RNA and RNA-protein assemblies, but interpretation of such data often requires previous structural knowledge. Here we highlight how single-molecule tools can directly complement structural approaches for two processes--translation and reverse transcription-to provide a dynamic view of molecular function.

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.

Energy Landscape and Pathways for Transitions between Watson-Crick and Hoogsteen Base Pairing in DNA

The recent discovery that Hoogsteen (HG) base pairs are widespread in DNA across diverse sequences and positional contexts could have important implications for understanding DNA replication and DNA-protein recognition. While evidence is emerging that the Hoogsteen conformation could be a thermodynamically accessible conformation of the DNA duplex and provide a means to expand its functionality, relatively little is known about the molecular mechanism underlying the Watson-Crick (WC) to HG transition. In this Perspective, we describe pathways and kinetics for this transition at an atomic level of detail, using the energy landscape perspective. We show that competition between the duplex conformations results in a double funnel landscape, which explains some recent experimental observations. The interconversion pathways feature a number of intermediates, with a variable number of WC and HG base pairs. The relatively slow kinetics, with possible deviations from two-state behavior, suggest that this conformational switch is likely to be a challenging target for both simulation and experiment.

A coarse-grained model for assisting the investigation of structure and dynamics of large nucleic acids by ion mobility spectrometry-mass spectrometry

Ion Mobility Spectrometry-Mass Spectrometry (IMS-MS) is a rapidly emerging tool for the investigation of nucleic acid structure and dynamics. IMS-MS determinations can provide valuable information regarding alternative topologies, folding intermediates, and conformational heterogeneities, which are not readily accessible to other analytical techniques. The leading strategies for data interpretation rely on computational and experimental approaches to correctly assign experimental observations to putative structures. A very effective strategy involves the application of molecular dynamics (MD) simulations to predict the structure of the analyte molecule, calculate its collision cross section (CCS), and then compare this computational value with the corresponding experimental data. While this approach works well for small nucleic acid species, analyzing larger nucleic acids of biological interest is hampered by the computational cost associated with capturing their extensive structure and dynamics in all-atom detail. In this report, we describe the implementation of a coarse graining (CG) approach to reduce the cost of the computational methods employed in the data interpretation workflow. Our framework employs a five-bead model to accurately represent each nucleotide in the nucleic acid structure. The beads are appropriately parameterized to enable the direct calculation of CCS values from CG models, thus affording the ability to pursue the analysis of larger, highly dynamic constructs. The validity of this approach was successfully confirmed by the excellent correlation between the CCS values obtained in parallel by all-atom and CG workflows.

Integrated structural biology to unravel molecular mechanisms of protein-RNA recognition

Recent advances in RNA sequencing technologies have greatly expanded our knowledge of the RNA landscape in cells, often with spatiotemporal resolution. These techniques identified many new (often non-coding) RNA molecules. Large-scale studies have also discovered novel RNA binding proteins (RBPs), which exhibit single or multiple RNA binding domains (RBDs) for recognition of specific sequence or structured motifs in RNA. Starting from these large-scale approaches it is crucial to unravel the molecular principles of protein-RNA recognition in ribonucleoprotein complexes (RNPs) to understand the underlying mechanisms of gene regulation. Structural biology and biophysical studies at highest possible resolution are key to elucidate molecular mechanisms of RNA recognition by RBPs and how conformational dynamics, weak interactions and cooperative binding contribute to the formation of specific, context-dependent RNPs. While large compact RNPs can be well studied by X-ray crystallography and cryo-EM, analysis of dynamics and weak interaction necessitates the use of solution methods to capture these properties. Here, we illustrate methods to study the structure and conformational dynamics of protein-RNA complexes in solution starting from the identification of interaction partners in a given RNP. Biophysical and biochemical techniques support the characterization of a protein-RNA complex and identify regions relevant in structural analysis. Nuclear magnetic resonance (NMR) is a powerful tool to gain information on folding, stability and dynamics of RNAs and characterize RNPs in solution. It provides crucial information that is complementary to the static pictures derived from other techniques. NMR can be readily combined with other solution techniques, such as small angle X-ray and/or neutron scattering (SAXS/SANS), electron paramagnetic resonance (EPR), and Förster resonance energy transfer (FRET), which provide information about overall shapes, internal domain arrangements and dynamics. Principles of protein-RNA recognition and current approaches are reviewed and illustrated with recent studies.

Simultaneous NMR characterisation of multiple minima in the free energy landscape of an RNA UUCG tetraloop

RNA molecules in solution tend to undergo structural fluctuations of relatively large amplitude and to populate a range of different conformations some of which with low populations. It is still very challenging, however, to characterise the structures of these low populated states and to understand their functional roles. In the present study, we address this problem by using NMR residual dipolar couplings (RDCs) as structural restraints in replica-averaged metadynamics (RAM) simulations. By applying this approach to a 14-mer RNA hairpin containing the prototypical UUCG tetraloop motif, we show that it is possible to construct the free energy landscape of this RNA molecule. This free energy landscapes reveals the surprisingly rich dynamics of the UUCG tetraloop and identifies the multiple substates that exist in equilibrium owing to thermal fluctuations. The approach that we present is general and can be applied to the study of the free energy landscapes of other RNA or RNA-protein systems.

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.

NMR solution structure determination of large RNA-protein complexes

Structure determination of RNA-protein complexes is essential for our understanding of the multiple layers of RNA-mediated posttranscriptional regulation of gene expression. Over the past 20years, NMR spectroscopy became a key tool for structural studies of RNA-protein interactions. Here, we review the progress being made in NMR structure determination of large ribonucleoprotein assemblies. We discuss approaches for the design of RNA-protein complexes for NMR structural studies, established and emerging isotope and segmental labeling schemes suitable for large RNPs and how to gain distance restraints from NOEs, PREs and EPR and orientational information from RDCs and SAXS/SANS in such systems. The new combination of NMR measurements with MD simulations and its potential will also be discussed. Application and combination of these various methods for structure determination of large RNPs will be illustrated with three large RNA-protein complexes (>40kDa) and other interesting complexes determined in the past six and a half years.

Structure of a low-population binding intermediate in protein-RNA recognition

The interaction of the HIV-1 protein transactivator of transcription (Tat) and its cognate transactivation response element (TAR) RNA transactivates viral transcription and represents a paradigm for the widespread occurrence of conformational rearrangements in protein-RNA recognition. Although the structures of free and bound forms of TAR are well characterized, the conformations of the intermediates in the binding process are still unknown. By determining the free energy landscape of the complex using NMR residual dipolar couplings in replica-averaged metadynamics simulations, we observe two low-population intermediates. We then rationally design two mutants, one in the protein and another in the RNA, that weaken specific nonnative interactions that stabilize one of the intermediates. By using surface plasmon resonance, we show that these mutations lower the release rate of Tat, as predicted. These results identify the structure of an intermediate for RNA-protein binding and illustrate a general strategy to achieve this goal with high resolution.

A dynamic look backward and forward

The 2015 Gunther Laukien Prize recognized solution NMR studies of protein dynamics and thermodynamics. This Perspective surveys aspects of the development and application of NMR spin relaxation for investigations of protein flexibility and function over multiple time scales in solution. Methods highlighted include analysis of overall rotational diffusion, theoretical descriptions of R1ρ relaxation, and molecular dynamics simulations to interpret NMR spin relaxation. Applications are illustrated for the zinc-finger domain Xfin-31, the calcium-binding proteins calbindin D9k and calmodulin, and the bZip transcription factor of GCN4.

The solution structural ensembles of RNA kink-turn motifs and their protein complexes

With the growing number of crystal structures of RNA and RNA-protein complexes, a critical next step is understanding the dynamic solution behavior of these entities in terms of conformational ensembles and energy landscapes. To this end, we have used X-ray scattering interferometry (XSI) to probe the ubiquitous RNA kink-turn motif and its complexes with the canonical kink-turn binding protein L7Ae. XSI revealed that the folded kink-turn is best described as a restricted conformational ensemble. The ions present in solution alter the nature of this ensemble, and protein binding can perturb the kink-turn ensemble without collapsing it to a unique state. This study demonstrates how XSI can reveal structural and ensemble properties of RNAs and RNA-protein complexes and uncovers the behavior of an important RNA-protein motif. This type of information will be necessary to understand, predict and engineer the behavior and function of RNAs and their protein complexes.

RNA recognition by Roquin in posttranscriptional gene regulation

Posttranscriptional regulation of gene expression plays a central role in the initiation of innate and adaptive immune responses. This is exemplified by the protein Roquin, which has attracted great interest during the past decade owing to its ability to prevent autoimmunity. Roquin controls T-cell activation and T helper cell differentiation by limiting the induced expression of costimulatory receptors on the surface of T cells. It does so by recognizing cis regulatory RNA-hairpin elements in the 3' UTR of target transcripts via its ROQ domain-a novel RNA-binding fold-and triggering their degradation through recruitment of factors that mediate deadenylation and decapping. Recent structural studies have revealed molecular details of the recognition of RNA hairpin structures by the ROQ domain. Surprisingly, it was found that Roquin mainly relies on shape-specific recognition of the RNA. This observation implies that a much broader range of RNA motifs could interact with the protein, but it also complicates systematic searches for novel mRNA targets of Roquin. Thus, large-scale approaches, such as crosslinking and immunoprecipitation or systematic evolution of ligands by exponential enrichment experiments coupled with next-generation sequencing, will be required to identify the complete spectrum of its target RNAs. Together with structural analyses of their binding modes, this will enable us to unravel the intricate complexity of 3' UTR regulation by Roquin and other trans-acting factors. Here, we review our current understanding of Roquin-RNA interactions and their role for Roquin function. WIREs RNA 2016, 7:455-469. doi: 10.1002/wrna.1333 For further resources related to this article, please visit the WIREs website.

Impact of nucleic acid self-alignment in a strong magnetic field on the interpretation of indirect spin-spin interactions

Heteronuclear and homonuclear direct (D) and indirect (J) spin-spin interactions are important sources of structural information about nucleic acids (NAs). The Hamiltonians for the D and J interactions have the same functional form; thus, the experimentally measured apparent spin-spin coupling constant corresponds to a sum of J and D. In biomolecular NMR studies, it is commonly presumed that the dipolar contributions to Js are effectively canceled due to random molecular tumbling. However, in strong magnetic fields, such as those employed for NMR analysis, the tumbling of NA fragments is anisotropic because the inherent magnetic susceptibility of NAs causes an interaction with the external magnetic field. This motional anisotropy is responsible for non-zero D contributions to Js. Here, we calculated the field-induced D contributions to 33 structurally relevant scalar coupling constants as a function of magnetic field strength, temperature and NA fragment size. We identified two classes of Js, namely (1)JCH and (3)JHH couplings, whose quantitative interpretation is notably biased by NA motional anisotropy. For these couplings, the magnetic field-induced dipolar contributions were found to exceed the typical experimental error in J-coupling determinations by a factor of two or more and to produce considerable over- or under-estimations of the J coupling-related torsion angles, especially at magnetic field strengths >12 T and for NA fragments longer than 12 bp. We show that if the non-zero D contributions to J are not properly accounted for, they might cause structural artifacts/bias in NA studies that use solution NMR spectroscopy.

Conformational Selection in a Protein-Protein Interaction Revealed by Dynamic Pathway Analysis

Molecular recognition plays a central role in biology, and protein dynamics has been acknowledged to be important in this process. However, it is highly debated whether conformational changes happen before ligand binding to produce a binding-competent state (conformational selection) or are caused in response to ligand binding (induced fit). Proposals for both mechanisms in protein/protein recognition have been primarily based on structural arguments. However, the distinction between them is a question of the probabilities of going via these two opposing pathways. Here, we present a direct demonstration of exclusive conformational selection in protein/protein recognition by measuring the flux for rhodopsin kinase binding to its regulator recoverin, an important molecular recognition in the vision system. Using nuclear magnetic resonance (NMR) spectroscopy, stopped-flow kinetics, and isothermal titration calorimetry, we show that recoverin populates a minor conformation in solution that exposes a hydrophobic binding pocket responsible for binding rhodopsin kinase. Protein dynamics in free recoverin limits the overall rate of binding.

Rich RNA Structure Landscapes Revealed by Mutate-and-Map Analysis

Landscapes exhibiting multiple secondary structures arise in natural RNA molecules that modulate gene expression, protein synthesis, and viral infection [corrected]. We report herein that high-throughput chemical experiments can isolate an RNA's multiple alternative secondary structures as they are stabilized by systematic mutagenesis (mutate-and-map, M2) and that a computational algorithm, REEFFIT, enables unbiased reconstruction of these states' structures and populations. In an in silico benchmark on non-coding RNAs with complex landscapes, M2-REEFFIT recovers 95% of RNA helices present with at least 25% population while maintaining a low false discovery rate (10%) and conservative error estimates. In experimental benchmarks, M2-REEFFIT recovers the structure landscapes of a 35-nt MedLoop hairpin, a 110-nt 16S rRNA four-way junction with an excited state, a 25-nt bistable hairpin, and a 112-nt three-state adenine riboswitch with its expression platform, molecules whose characterization previously required expert mutational analysis and specialized NMR or chemical mapping experiments. With this validation, M2-REEFFIT enabled tests of whether artificial RNA sequences might exhibit complex landscapes in the absence of explicit design. An artificial flavin mononucleotide riboswitch and a randomly generated RNA sequence are found to interconvert between three or more states, including structures for which there was no design, but that could be stabilized through mutations. These results highlight the likely pervasiveness of rich landscapes with multiple secondary structures in both natural and artificial RNAs and demonstrate an automated chemical/computational route for their empirical characterization.

Conformational dynamics of nucleic acid molecules studied by PELDOR spectroscopy with rigid spin labels

Nucleic acid molecules can adopt a variety of structures and exhibit a large degree of conformational flexibility to fulfill their various functions in cells. Here we describe the use of Pulsed Electron-Electron Double Resonance (PELDOR or DEER) to investigate nucleic acid molecules where two cytosine analogs have been incorporated as spin probes. Because these new types of spin labels are rigid and incorporated into double stranded DNA and RNA molecules, there is no additional flexibility of the spin label itself present. Therefore the magnetic dipole-dipole interaction between both spin labels encodes for the distance as well as for the mutual orientation between the spin labels. All of this information can be extracted by multi-frequency/multi-field PELDOR experiments, which gives very precise and valuable information about the structure and conformational flexibility of the nucleic acid molecules. We describe in detail our procedure to obtain the conformational ensembles and show the accuracy and limitations with test examples and application to double-stranded DNA.

TEMPO-derived spin labels linked to the nucleobases adenine and cytosine for probing local structural perturbations in DNA by EPR spectroscopy

Three 2´-deoxynucleosides containing semi-flexible spin labels, namely (T)A, (U)A and (U)C, were prepared and incorporated into deoxyoligonucleotides using the phosphoramidite method. All three nucleosides contain 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) connected to the exocyclic amino group; (T)A directly and (U)A as well as (U)C through a urea linkage. (T)A and (U)C showed a minor destabilization of a DNA duplex, as registered by a small decrease in the melting temperature, while (U)A destabilized the duplex by more than 10 °C. Circular dichroism (CD) measurements indicated that all three labels were accommodated in B-DNA duplex. The mobility of the spin label (T)A varied with different base-pairing partners in duplex DNA, with the (T)A•T pair being the least mobile. Furthermore, (T)A showed decreased mobility under acidic conditions for the sequences (T)A•C and (T)A•G, to the extent that the EPR spectrum of the latter became nearly superimposable to that of (T)A•T. The reduced mobility of the (T)A•C and (T)A•G mismatches at pH 5 is consistent with the formation of (T)AH(+)•C and (T)AH(+)•G, in which protonation of N1 of A allows the formation of an additional hydrogen bond to N3 of C and N7 of G, respectively, with G in a syn-conformation. The urea-based spin labels (U)A and (U)C were more mobile than (T)A, but still showed a minor variation in their EPR spectra when paired with A, G, C or T in a DNA duplex. (U)A and (U)C had similar mobility order for the different base pairs, with the lowest mobility when paired with C and the highest when paired with T.

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.

Evaluating the uncertainty in exchange parameters determined from off-resonance R1ρ relaxation dispersion for systems in fast exchange

Spin relaxation in the rotating frame (R1ρ) is a powerful NMR technique for characterizing fast microsecond timescale exchange processes directed toward short-lived excited states in biomolecules. At the limit of fast exchange, only k(ex)=k(1)+k(-1) and Φex=p(G)p(E)(Δω)(2) can be determined from R1ρ data limiting the ability to characterize the structure and energetics of the excited state conformation. Here, we use simulations to examine the uncertainty with which exchange parameters can be determined for two state systems in intermediate-to-fast exchange using off-resonance R1ρ relaxation dispersion. R1ρ data computed by solving the Bloch-McConnell equations reveals small but significant asymmetry with respect to offset (R1ρ (ΔΩ)≠R1ρ (-ΔΩ)), which is a hallmark of slow-to-intermediate exchange, even under conditions of fast exchange for free precession chemical exchange line broadening (k(ex)/Δω>10). A grid search analysis combined with bootstrap and Monte-Carlo based statistical approaches for estimating uncertainty in exchange parameters reveals that both the sign and magnitude of Δω can be determined at a useful level of uncertainty for systems in fast exchange (k(ex)/Δω<10) but that this depends on the uncertainty in the R1ρ data and requires a thorough examination of the multidimensional variation of χ(2) as a function of exchange parameters. Results from simulations are complemented by analysis of experimental R1ρ data measured in three nucleic acid systems with exchange processes occurring on the slow (k(ex)/Δω=0.2; pE=∼0.7%), fast (k(ex)/Δω=∼10-16; p(E)=∼13%) and very fast (k(ex)=39,000 s(-1)) chemical shift timescales.

A tensor-free method for the structural and dynamical refinement of proteins using residual dipolar couplings

Residual dipolar couplings (RDCs) are parameters measured in nuclear magnetic resonance spectroscopy that can provide exquisitely detailed information about the structure and dynamics of biological macromolecules. We describe here a method of using RDCs for the structural and dynamical refinement of proteins that is based on the observation that the RDC between two atomic nuclei depends directly on the angle ϑ between the internuclear vector and the external magnetic field. For every pair of nuclei for which an RDC is available experimentally, we introduce a structural restraint to minimize the deviation from the value of the angle ϑ derived from the measured RDC and that calculated in the refinement protocol. As each restraint involves only the calculation of the angle ϑ of the corresponding internuclear vector, the method does not require the definition of an overall alignment tensor to describe the preferred orientation of the protein with respect to the alignment medium. Application to the case of ubiquitin demonstrates that this method enables an accurate refinement of the structure and dynamics of this protein to be obtained.

Chemical exchange in biomacromolecules: past, present, and future

The perspective reviews quantitative investigations of chemical exchange phenomena in proteins and other biological macromolecules using NMR spectroscopy, particularly relaxation dispersion methods. The emphasis is on techniques and applications that quantify the populations, interconversion kinetics, and structural features of sparsely populated conformational states in equilibrium with a highly populated ground state. Applications to folding, molecular recognition, catalysis, and allostery by proteins and nucleic acids are highlighted.

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

Complex RNA structures are constructed from helical segments connected by flexible loops that move spontaneously and in response to binding of small molecule ligands and proteins. Understanding the conformational variability of RNA requires the characterization of the coupled time evolution of interconnected flexible domains. To elucidate the collective molecular motions and explore the conformational landscape of the HIV-1 TAR RNA, we describe a new methodology that utilizes energy-minimized structures generated by the program “Fragment Assembly of RNA with Full-Atom Refinement (FARFAR)”. We apply structural filters in the form of experimental residual dipolar couplings (RDCs) to select a subset of discrete energy-minimized conformers and carry out principal component analyses (PCA) to corroborate the choice of the filtered subset. We use this subset of structures to calculate solution T₁ and T_1ρ relaxation times for ¹³C spins in multiple residues in different domains of the molecule using two simulation protocols that we previously published. We match the experimental T₁ times to within 2% and the T_1ρ times to within less than 10% for helical residues. These results introduce a protocol to construct viable dynamic trajectories for RNA molecules that accord well with experimental NMR data and support the notion that the motions of the helical portions of this small RNA can be described by a relatively small number of discrete conformations exchanging over time scales longer than 1 μs.