Resolving the motional modes that code for RNA adaptation

Sensitivity-enhanced NMR 15N R1 and R1ρ relaxation experiments for the investigation of intrinsically disordered proteins at high magnetic fields

NMR relaxation experiments provide residue-specific insights into the structural dynamics of proteins. Here, we present an optimized set of sensitivity-enhanced 15N R1 and R1ρ relaxation experiments applicable to fully protonated proteins. The NMR pulse sequences are conceptually similar to the set of TROSY-based sequences and their HSQC counterpart (Lakomek et al., J. Biomol. NMR 2012). Instead of the TROSY read-out scheme, a sensitivity-enhanced HSQC read-out scheme is used, with improved and easier optimized water suppression. The presented pulse sequences are applied on the cytoplasmic domain of the SNARE protein Synpatobrevin-2 (Syb-2), which is intrinsically disordered in its monomeric pre-fusion state. A two-fold increase in the obtained signal-to-noise ratio is observed for this intrinsically disordered protein, therefore offering a four-fold reduction of measurement time compared to the TROSY-detected version. The inter-scan recovery delay can be shortened to two seconds. Pulse sequences were tested at 600 MHz and 1200 MHz 1H Larmor frequency, thus applicable over a wide magnetic field range. A comparison between protonated and deuterated protein samples reveals high agreement, indicating that reliable 15N R1 and R1ρ rate constants can be extracted for fully protonated and deuterated samples. The presented pulse sequences will benefit not only for IDPs but also for an entire range of low and medium-sized proteins.

Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods

Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.

Kinetic Resolution of the Atomic 3D Structures Formed by Ground and Excited Conformational States in an RNA Dynamic Ensemble

Knowing the 3D structures formed by the various conformations populating the RNA free-energy landscape, their relative abundance, and kinetic interconversion rates is required to obtain a quantitative and predictive understanding of how RNAs fold and function at the atomic level. While methods integrating ensemble-averaged experimental data with computational modeling are helping define the most abundant conformations in RNA ensembles, elucidating their kinetic rates of interconversion and determining the 3D structures of sparsely populated short-lived RNA excited conformational states (ESs) remains challenging. Here, we developed an approach integrating Rosetta-FARFAR RNA structure prediction with NMR residual dipolar couplings and relaxation dispersion that simultaneously determines the 3D structures formed by the ground-state (GS) and ES subensembles, their relative abundance, and kinetic rates of interconversion. The approach is demonstrated on HIV-1 TAR, whose six-nucleotide apical loop was previously shown to form a sparsely populated (∼13%) short-lived (lifetime ∼ 45 μs) ES. In the GS, the apical loop forms a broad distribution of open conformations interconverting on the pico-to-nanosecond time scale. Most residues are unpaired and preorganized to bind the Tat-superelongation protein complex. The apical loop zips up in the ES, forming a narrow distribution of closed conformations, which sequester critical residues required for protein recognition. Our work introduces an approach for determining the 3D ensemble models formed by sparsely populated RNA conformational states, provides a rare atomic view of an RNA ES, and kinetically resolves the atomic 3D structures of RNA conformational substates, interchanging on time scales spanning 6 orders of magnitude, from picoseconds to microseconds.

Elongated Bacterial Pili as a Versatile Alignment Medium for NMR Spectroscopy

In NMR spectroscopy, residual dipolar couplings (RDCs) have emerged as one of the most exquisite probes of biological structure and dynamics. The measurement of RDCs relies on the partial alignment of the molecule of interest, for example by using a liquid crystal as a solvent. Here, we establish bacterial type 1 pili as an alternative liquid-crystalline alignment medium for the measurement of RDCs. To achieve alignment at pilus concentrations that allow for efficient NMR sample preparation, we elongated wild-type pili by recombinant overproduction of the main structural pilus subunit. Building on the extraordinary stability of type 1 pili against spontaneous dissociation and unfolding, we show that the medium is compatible with challenging experimental conditions such as high temperature, the presence of detergents, organic solvents or very acidic pH, setting it apart from most established alignment media. Using human ubiquitin, HIV-1 TAR RNA and camphor as spectroscopic probes, we demonstrate the applicability of the medium for the determination of RDCs of proteins, nucleic acids and small molecules. Our results show that type 1 pili represent a very useful alternative to existing alignment media and may readily assist the characterization of molecular structure and dynamics by NMR.

DNA base order parameter determination without influence of chemical exchange

Nuclear magnetic resonance (NMR) spectroscopy is a versatile tool used to investigate the dynamic properties of biological macromolecules and their complexes. NMR relaxation data can provide order parameters S2, which represent the mobility of bond vectors reorienting within a molecular frame. Determination of S2 parameters typically involves the use of transverse NMR relaxation rates. However, the accuracy in S2 determination can be diminished by elevation of the transverse relaxation rates through conformational or chemical exchange involving protonation/deprotonation or non-Watson-Crick base-pair states of nucleic acids. Here, we propose an approach for determination of S2 parameters without the influence of exchange processes. This approach utilizes transverse and longitudinal 13C chemical shift anisotropy (CSA) - dipole-dipole (DD) cross-correlation rates instead of 13C transverse relaxation rates. Anisotropy in rotational diffusion is taken into consideration. An application of this approach to nucleotide base CH groups of a uniformly 13C/15N-labeled DNA duplex is demonstrated.

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.

Biomotors, viral assembly, and RNA nanobiotechnology: Current achievements and future directions

The International Society of RNA Nanotechnology and Nanomedicine (ISRNN) serves to further the development of a wide variety of functional nucleic acids and other related nanotechnology platforms. To aid in the dissemination of the most recent advancements, a biennial discussion focused on biomotors, viral assembly, and RNA nanobiotechnology has been established where international experts in interdisciplinary fields such as structural biology, biophysical chemistry, nanotechnology, cell and cancer biology, and pharmacology share their latest accomplishments and future perspectives. The results summarized here highlight advancements in our understanding of viral biology and the structure-function relationship of frame-shifting elements in genomic viral RNA, improvements in the predictions of SHAPE analysis of 3D RNA structures, and the understanding of dynamic RNA structures through a variety of experimental and computational means. Additionally, recent advances in the drug delivery, vaccine design, nanopore technologies, biomotor and biomachine development, DNA packaging, RNA nanotechnology, and drug delivery are included in this critical review. We emphasize some of the novel accomplishments, major discussion topics, and present current challenges and perspectives of these emerging fields.

Targeting RNA structures with small molecules

RNA adopts 3D structures that confer varied functional roles in human biology and dysfunction in disease. Approaches to therapeutically target RNA structures with small molecules are being actively pursued, aided by key advances in the field including the development of computational tools that predict evolutionarily conserved RNA structures, as well as strategies that expand mode of action and facilitate interactions with cellular machinery. Existing RNA-targeted small molecules use a range of mechanisms including directing splicing - by acting as molecular glues with cellular proteins (such as branaplam and the FDA-approved risdiplam), inhibition of translation of undruggable proteins and deactivation of functional structures in noncoding RNAs. Here, we describe strategies to identify, validate and optimize small molecules that target the functional transcriptome, laying out a roadmap to advance these agents into the next decade.

The dynamic, motile and deformative properties of RNA nanoparticles facilitate the third milestone of drug development

Besides mRNA, rRNA, and tRNA, cells contain many other noncoding RNA that display critical roles in the regulation of cellular functions. Human genome sequencing revealed that the majority of non-protein-coding DNA actually codes for non-coding RNAs. The dynamic nature of RNA results in its motile and deformative behavior. These conformational transitions such as the change of base-pairing, breathing within complemented strands, and pseudoknot formation at the 2D level as well as the induced-fit and conformational capture at the 3D level are important for their biological functions including regulation, translation, and catalysis. The dynamic, motile and catalytic activity has led to a belief that RNA is the origin of life. We have recently reported that the deformative property of RNA nanoparticles enhances their penetration through the leaky blood vessel of cancers which leads to highly efficient tumor accumulation. This special deformative property also enables RNA nanoparticles to pass the glomerulus, overcoming the filtration size limit, resulting in fast renal excretion and rapid body clearance, thus low or no toxicity. The biodistribution of RNA nanoparticles can be further improved by the incorporation of ligands for cancer targeting. In addition to the favorable biodistribution profiles, RNA nanoparticles possess other properties including self-assembly, negative charge, programmability, and multivalency; making it a great material for pharmaceutical applications. The intrinsic negative charge of RNA nanoparticles decreases the toxicity of drugs by preventing nonspecific binding to the negative charged cell membrane and enhancing the solubility of hydrophobic drugs. The polyvalent property of RNA nanoparticles allows the multi-functionalization which can apply to overcome drug resistance. This review focuses on the summary of these unique properties of RNA nanoparticles, which describes the mechanism of RNA dynamic, motile and deformative properties, and elucidates and prepares to welcome the RNA therapeutics as the third milestone in pharmaceutical drug development.

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

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

2-Amino-1,3-benzothiazole-6-carboxamide Preferentially Binds the Tandem Mismatch Motif r(UY:GA)

RNA helices are often punctuated with non-Watson-Crick features that may be targeted by chemical compounds, but progress toward identifying such compounds has been slow. We embedded a tandem UU:GA mismatch motif (5'-UG-3':5'-AU-3') within an RNA hairpin stem to identify compounds that bind the motif specifically. The three-dimensional structure of the RNA hairpin and its interaction with a small molecule identified through virtual screening are presented. The G-A mismatch forms a sheared pair upon which the U-U base pair stacks. The hydrogen bond configuration of the U-U pair involves O2 of the U adjacent to the G and O4 of the U adjacent to the A. The G-A and U-U pairs are flanked by A-U and G-C base pairs, respectively, and the stability of the mismatch is greater than when the motif is within the context of other flanking base pairs or when the 5'-3' orientation of the G-A and U-U pairs is swapped. Residual dipolar coupling constants were used to generate an ensemble of structures against which a virtual screen of 64480 small molecules was performed. The tandem mismatch was found to be specific for one compound, 2-amino-1,3-benzothiazole-6-carboxamide, which binds with moderate affinity but extends the motif to include the flanking A-U and G-C base pairs. The finding that the affinity for the UU:GA mismatch is dependent on flanking sequence emphasizes the importance of the motif context and potentially increases the number of small noncanonical features within RNA that can be specifically targeted by small molecules.

Hierarchical Conformational Dynamics Confers Thermal Adaptability to preQ1 RNA Riboswitches

Single-stranded noncoding regulatory RNAs, as exemplified by bacterial riboswitches, are highly dynamic. The conformational dynamics allow the riboswitch to reach maximum switching efficiency under appropriate conditions. Here we characterize the conformational dynamics of preQ₁ riboswitches from mesophilic and thermophilic bacterial species at various temperatures. With the integrative use of small-angle X-ray scattering, NMR, and molecular dynamics simulations, we model the ensemble-structures of the preQ₁ riboswitch aptamers without or with a ligand bound. We show that the preQ₁ riboswitch is sufficiently dynamic and fluctuating among multiple folding intermediates only near the physiological temperature of the microorganism. The hierarchical folding dynamics of the RNA involves the docking of 3'-tail to form a second RNA helix and the helical stacking to form an H-type pseudoknot structure. Further, we show that RNA secondary and tertiary dynamics can be modulated by temperature and by the length of an internal loop. The coupled equilibria between RNA folding intermediates are essential for preQ₁ binding, and a four-state exchange model can account for the change of ligand-triggered switching efficiency with temperature. Together, we have established a relationship between the hierarchical dynamics and riboswitch function, and illustrated how the RNA adapts to high temperature.

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.

Exposing Hidden High-Affinity RNA Conformational States

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

Protein Dynamics from Accurate Low-Field Site-Specific Longitudinal and Transverse Nuclear Spin Relaxation

Nuclear magnetic relaxation provides invaluable quantitative site-specific information on the dynamics of complex systems. Determining dynamics on nanosecond time scales requires relaxation measurements at low magnetic fields incompatible with high-resolution NMR. Here, we use a two-field NMR spectrometer to measure carbon-13 transverse and longitudinal relaxation rates at a field as low as 0.33 T (proton Larmor frequency 14 MHz) in specifically labeled side chains of the protein ubiquitin. The use of radiofrequency pulses enhances the accuracy of measurements as compared to high-resolution relaxometry approaches, where the sample is moved in the stray field of the superconducting magnet. Importantly, we demonstrate that accurate measurements at a single low magnetic field provide enough information to characterize complex motions on low nanosecond time scales, which opens a new window for the determination of site-specific nanosecond motions in complex systems such as proteins.

Functional protein dynamics on uncharted time scales detected by nanoparticle-assisted NMR spin relaxation

Protein function depends critically on intrinsic internal dynamics, which is manifested in distinct ways, such as loop motions that regulate protein recognition and catalysis. Under physiological conditions, dynamic processes occur on a wide range of time scales from subpicoseconds to seconds. Commonly used NMR spin relaxation in solution provides valuable information on very fast and slow motions but is insensitive to the intermediate nanosecond to microsecond range that exceeds the protein tumbling correlation time. Presently, very little is known about the nature and functional role of these motions. It is demonstrated here how transverse spin relaxation becomes exquisitely sensitive to these motions at atomic resolution when studying proteins in the presence of nanoparticles. Application of this novel cross-disciplinary approach reveals large-scale dynamics of loops involved in functionally critical protein-protein interactions and protein-calcium ion recognition that were previously unobservable.

An NMR View of Protein Dynamics in Health and Disease

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

Long-Range RNA Structural Information via a Paramagnetically Tagged Reporter Protein

NMR has provided a wealth of structural and dynamical information for RNA molecules of up to ∼50 nucleotides, but its application to larger RNAs has been hampered in part by difficulties establishing global structural features. A potential solution involves measurement of NMR perturbations after site-specific paramagnetic labeling. Although the approach works well for proteins, the inability to place the label at specific sites has prevented its application to larger RNAs transcribed in vitro. Here, we present a strategy in which RNA loop residues are modified to promote binding to a paramagnetically tagged reporter protein. Lanthanide-induced pseudocontact shifts are demonstrated for a 232-nucleotide RNA bound to tagged derivatives of the spliceosomal U1A RNA-binding domain. Further, the method is validated with a 36-nucleotide RNA for which measured NMR values agreed with predictions based on the previously known protein and RNA structures. The ability to readily insert U1A binding sites into ubiquitous hairpin and/or loop structures should make this approach broadly applicable for the atomic-level study of large RNAs.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

Accurate Measurement of Residual Dipolar Couplings in Large RNAs by Variable Flip Angle NMR

NMR approaches using nucleotide-specific deuterium labeling schemes have enabled structural studies of biologically relevant RNAs of increasing size and complexity. Although local structure is well-determined using these methods, definition of global structural features, including relative orientations of independent helices, remains a challenge. Residual dipolar couplings, a potential source of orientation information, have not been obtainable for large RNAs due to poor sensitivity resulting from rapid heteronuclear signal decay. Here we report a novel multiple quantum NMR method for RDC determination that employs flip angle variation rather than a coupling evolution period. The accuracy of the method and its utility for establishing interhelical orientations are demonstrated for a 36-nucleotide RNA, for which comparative data could be obtained. Applied to a 78 kDa Rev response element from the HIV-1 virus, which has an effective rotational correlation time of ca. 160 ns, the method yields sensitivity gains of an order of magnitude or greater over existing approaches. Solution-state access to structural organization in RNAs of at least 230 nucleotides is now possible.

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

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

CCR5 RNA Pseudoknots: Residue and Site-Specific Labeling correlate Internal Motions with microRNA Binding

Conformational dynamics of RNA molecules play a critical role in governing their biological functions. Measurements of RNA dynamic behavior sheds important light on sites that interact with their binding partners or cellular stimulators. However, such measurements using solution-state NMR are difficult for large RNA molecules (>70 nt; nt=nucleotides) owing to severe spectral overlap, homonuclear 13 C scalar couplings, and line broadening. Herein, a strategic combination of solid-phase synthesis, site-specific isotopic labeled phosphoramidites, and enzymatic ligation is introduced. This approach allowed the position-specific insertion of isotopic probes into a 96 nt CCR5 RNA fragment. Accurate measurements of functional dynamics using the Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion (RD) experiments enabled extraction of the exchange rates and populations of this RNA. NMR chemical shift perturbation analysis of the RNA/microRNA-1224 complex indicated that A90-C1' of the pseudoknot exhibits similar changes in chemical shift observed in the excited state. This work demonstrates the general applicability of a NMR-labeling strategy to probe functional RNA structural dynamics.

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.

Insights into Watson-Crick/Hoogsteen breathing dynamics and damage repair from the solution structure and dynamic ensemble of DNA duplexes containing m1A

In the canonical DNA double helix, Watson-Crick (WC) base pairs (bps) exist in dynamic equilibrium with sparsely populated (∼0.02-0.4%) and short-lived (lifetimes ∼0.2-2.5 ms) Hoogsteen (HG) bps. To gain insights into transient HG bps, we used solution-state nuclear magnetic resonance spectroscopy, including measurements of residual dipolar couplings and molecular dynamics simulations, to examine how a single HG bp trapped using the N1-methylated adenine (m1A) lesion affects the structural and dynamic properties of two duplexes. The solution structure and dynamic ensembles of the duplexes reveals that in both cases, m1A forms a m1A•T HG bp, which is accompanied by local and global structural and dynamic perturbations in the double helix. These include a bias toward the BI backbone conformation; sugar repuckering, major-groove directed kinking (∼9°); and local melting of neighboring WC bps. These results provide atomic insights into WC/HG breathing dynamics in unmodified DNA duplexes as well as identify structural and dynamic signatures that could play roles in m1A recognition and repair.

Structure of RNA Stem Loop B from the Picornavirus Replication Platform

The presumptive RNA cloverleaf at the start of the 5'-untranslated region of the picornavirus genome is an essential element in replication. Stem loop B (SLB) of the cloverleaf is a recognition site for the host polyC-binding protein, which initiates a switch from translation to replication. Here we present the solution structure of human rhinovirus isotype 14 SLB using nuclear magnetic resonance spectroscopy. SLB adopts a predominantly A-form helical structure. The stem contains five Watson-Crick base pairs and one wobble base pair and is capped by an eight-nucleotide loop. The wobble base pair introduces perturbations into the helical parameters but does not appear to introduce flexibility. However, the helix major groove appears to be accessible. Flexibility is seen throughout the loop and in the terminal nucleotides. The pyrimidine-rich region of the loop, the apparent recognition site for the polyC-binding protein, is the most disordered region of the structure.

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.

Inter-helical conformational preferences of HIV-1 TAR-RNA from maximum occurrence analysis of NMR data and molecular dynamics simulations

Detecting conformational heterogeneity in biological macromolecules is a key for the understanding of their biological function. We here provide a comparison between two independent approaches to assess conformational heterogeneity: molecular dynamics simulations, performed without inclusion of any experimental data, and maximum occurrence (MaxOcc) distribution over the topologically available conformational space. The latter only reflects the extent of the averaging and identifies regions which are most compliant with the experimentally measured NMR Residual Dipolar Couplings (RDCs). The analysis was performed for the HIV-1 TAR RNA, consisting of two helical domains connected by a flexible bulge junction, for which four sets of RDCs were available as well as an 8.2 μs all-atom molecular dynamics simulation. A sample and select approach was previously applied to extract from the molecular dynamics trajectory conformational ensembles in agreement with the four sets of RDCs. The MaxOcc analysis performed here identifies the most likely sampled region in the conformational space of the system which, strikingly, overlaps well with the structures independently sampled in the molecular dynamics calculations and even better with the RDC selected ensemble.

Structural biology of supramolecular assemblies by magic-angle spinning NMR spectroscopy

In recent years, exciting developments in instrument technology and experimental methodology have advanced the field of magic-angle spinning (MAS) nuclear magnetic resonance (NMR) to new heights. Contemporary MAS NMR yields atomic-level insights into structure and dynamics of an astounding range of biological systems, many of which cannot be studied by other methods. With the advent of fast MAS, proton detection, and novel pulse sequences, large supramolecular assemblies, such as cytoskeletal proteins and intact viruses, are now accessible for detailed analysis. In this review, we will discuss the current MAS NMR methodologies that enable characterization of complex biomolecular systems and will present examples of applications to several classes of assemblies comprising bacterial and mammalian cytoskeleton as well as human immunodeficiency virus 1 and bacteriophage viruses. The body of work reviewed herein is representative of the recent advancements in the field, with respect to the complexity of the systems studied, the quality of the data, and the significance to the biology.

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.

Shortening the HIV-1 TAR RNA Bulge by a Single Nucleotide Preserves Motional Modes over a Broad Range of Time Scales

Helix-junction-helix (HJH) motifs are flexible building blocks of RNA architecture that help define the orientation and dynamics of helical domains. They are also frequently involved in adaptive recognition of proteins and small molecules and in the formation of tertiary contacts. Here, we use a battery of nuclear magnetic resonance techniques to examine how deleting a single bulge residue (C24) from the human immunodeficiency virus type 1 (HIV-1) transactivation response element (TAR) trinucleotide bulge (U23-C24-U25) affects dynamics over a broad range of time scales. Shortening the bulge has an effect on picosecond-to-nanosecond interhelical and local bulge dynamics similar to that casued by increasing the Mg(2+) and Na(+) concentration, whereby a preexisting two-state equilibrium in TAR is shifted away from a bent flexible conformation toward a coaxial conformation, in which all three bulge residues are flipped out and flexible. Surprisingly, the point deletion minimally affects microsecond-to-millisecond conformational exchange directed toward two low-populated and short-lived excited conformational states that form through reshuffling of bases pairs throughout TAR. The mutant does, however, adopt a slightly different excited conformational state on the millisecond time scale, in which U23 is intrahelical, mimicking the expected conformation of residue C24 in the excited conformational state of wild-type TAR. Thus, minor changes in HJH topology preserve motional modes in RNA occurring over the picosecond-to-millisecond time scales but alter the relative populations of the sampled states or cause subtle changes in their conformational features.

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.

Insights into the mechanisms of RNA secondary structure destabilization by the HIV-1 nucleocapsid protein

The mature HIV-1 nucleocapsid protein NCp7 (NC) plays a key role in reverse transcription facilitating the two obligatory strand transfers. Several properties contribute to its efficient chaperon activity: preferential binding to single-stranded regions, nucleic acid aggregation, helix destabilization, and rapid dissociation from nucleic acids. However, little is known about the relationships between these different properties, which are complicated by the ability of the protein to recognize particular HIV-1 stem-loops, such as SL1, SL2, and SL3, with high affinity and without destabilizing them. These latter properties are important in the context of genome packaging, during which NC is part of the Gag precursor. We used NMR to investigate destabilization of the full-length TAR (trans activating response element) RNA by NC, which is involved in the first strand transfer step of reverse transcription. NC was used at a low protein:nucleotide (nt) ratio of 1:59 in these experiments. NMR data for the imino protons of TAR identified most of the base pairs destabilized by NC. These base pairs were adjacent to the loops in the upper part of the TAR hairpin rather than randomly distributed. Gel retardation assays showed that conversion from the initial TAR-cTAR complex to the fully annealed form occurred much more slowly at the 1:59 ratio than at the higher ratios classically used. Nevertheless, NC significantly accelerated the formation of the initial complex at a ratio of 1:59.

Tat-dependent production of an HIV-1 TAR-encoded miRNA-like small RNA

Evidence is accumulating that retroviruses can produce microRNAs (miRNAs). To prevent cleavage of their RNA genome, retroviruses have to use an alternative RNA source as miRNA precursor. The transacting responsive (TAR) hairpin structure in HIV-1 RNA has been suggested as source for miRNAs, but how these small RNAs are produced without impeding virus replication remained unclear. We used deep sequencing analysis of AGO2-bound HIV-1 RNAs to demonstrate that the 3′ side of the TAR hairpin is processed into a miRNA-like small RNA. This ∼21 nt RNA product is able to repress the expression of mRNAs bearing a complementary target sequence. Analysis of the small RNAs produced by wild-type and mutant HIV-1 variants revealed that non-processive transcription from the HIV-1 LTR promoter results in the production of short TAR RNAs that serve as precursor. These TAR RNAs are cleaved by Dicer and processing is stimulated by the viral Tat protein. This biogenesis pathway differs from the canonical miRNA pathway and allows HIV-1 to produce the TAR-encoded miRNA-like molecule without cleavage of the RNA genome.

Rapid Exchange Between Free and Bound States in RNA-Dendrimer Polyplexes: Implications on the Mechanism of Delivery and Release

A combination of solution NMR, dynamic light scattering (DLS), and fluorescence quenching assays were employed to obtain insights into the dynamics and structural features of a polyplex system consisting of HIV-1 transactivation response element (TAR) and PEGylated generation 5 poly(amidoamine) dendrimer (G5-PEG). NMR chemical shift mapping and (13)C spin relaxation based dynamics measurements depict the polyplex system as a highly dynamic assembly where the RNA, with its local structure and dynamics preserved, rapidly exchanges ( Collapse

Key Words

• therapeutic rna delivery
• pamam
• pegylation
• tar
• dls
• nmr spin relaxation

Collapse

MESH Headings

• Dendrimers/chemistry
• HIV-1/chemistry
• Magnetic Resonance Spectroscopy/methods
• Polyamines/chemistry
• Polyethylene Glycols/chemistry
• RNA/chemistry
• Response Elements/genetics
• Transfection/methods

Collapse

Grants

• R01 AI066975 NIAID NIH HHS
• P50 GM103297 NIGMS NIH HHS
• R01 AI087463-01 NIAID NIH HHS
• R01 EB005028 NIBIB NIH HHS
• R01EB005028 NIBIB NIH HHS
• R01 GM089846 NIGMS NIH HHS

Collapse

Affiliation(s)

Anisha Shakya Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Casey A. Dougherty Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Yi Xue Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Hashim M. Al-Hashimi Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Mark M. Banaszak Holl Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA

Collapse

Elucidating the Relation between Internal Motions and Dihedral Angles in an RNA Hairpin Using Molecular Dynamics

Molecular dynamics simulations were performed to characterize the internal motions of the ribonucleic acid apical stem loop of human hepatitis B virus. The NMR relaxation rates calculated directly from the trajectory are in good agreement with the experiment. Calculated order parameters follow the experimental pattern. Order parameters lower than 0.8 are observed for nucleotides that are weakly hydrogen bonded to their base pair partner, unpaired, or part of the loop. These residues show slow decay of the internal correlation functions of their base and sugar C-H vectors. Concerted motions around backbone dihedral angles influence the amplitude of motion of the sugar and base C-H vectors. The order parameters for base C-H vectors are also affected by the fluctuation of the glycosidic dihedral angle.

Charge screening in RNA: an integral route for dynamical enhancements

Electrostatic interactions of RNA are at the center of determining the dynamical flexibility and structural stability. By analysing neutron scattering spectroscopy, we show that fast dynamics of hydrated tRNA on ps to ns timescales increases with stronger charge screening, while its structural stability either increases or remains largely unchanged. An unprecedented electrostatic threshold for the onset of additional flexibility is induced from the correlation between the charge-screening density of counterions and the promoted dynamical properties. The results demonstrate that the enhanced dynamical flexibility of tRNA originates from local conformational relaxation coupled with stabilized charge screening rather than governed by fluctuation of hydrated counterions. The present study casts light on the specificity of electrostatic interactions in the thermodynamic balance between the dynamical flexibility and structural stability of RNA.

What NMR Relaxation Can Tell Us about the Internal Motion of an RNA Hairpin: A Molecular Dynamics Simulation Study

Classical molecular dynamics simulations of a 14-mer UUCG RNA hairpin are performed to study its conformational dynamics and corresponding NMR relaxation parameters. The direct calculation of the relaxation rates from the trajectory yields good agreement with experiment, indicating the validity of the theoretical model. Various ways to provide a link between theory and experiment are considered, including the "model-free approach" of Lipari and Szabo and Gaussian axial fluctuation model of Brüschweilwer. It is studied if the underlying assumptions of these approaches are satisfied in the case of a flexible RNA hairpin. Being consistent with the analysis of the NMR experiments, Lipari-Szabo fits of the first 100 or 1000 ps of the internal correlation functions lead to a nice agreement between calculated and experimental order parameters and internal correlation times. Finally, the relation between NMR order parameters and the underlying internal motion of the RNA hairpin is discussed in detail. A principal component analysis reveals that the principal motions of the molecule account only partially for the measured NMR order parameters, because the latter are insensitive to internal dynamics occurring on a nanosecond time scale due to molecular tumbling.

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

Stable isotope-labeled RNA phosphoramidites to facilitate dynamics by NMR

Given that Ribonucleic acids (RNAs) are a central hub of various cellular processes, methods to synthesize these RNAs for biophysical studies are much needed. Here, we showcase the applicability of 6-(13)C-pyrimidine phosphoramidites to introduce isolated (13)C-(1)H spin pairs into RNAs up to 40 nucleotides long. The method allows the incorporation of 6-(13)C-uridine and -cytidine residues at any desired position within a target RNA. By site-specific positioning of the (13)C-label using RNA solid phase synthesis, these stable isotope-labeling patterns are especially well suited to resolve resonance assignment ambiguities. Of even greater importance, the labeling pattern affords accurate quantification of important functional transitions of biologically relevant RNAs (e.g., riboswitch aptamer domains, viral RNAs, or ribozymes) in the μs- to ms time regime and beyond without complications of one bond carbon scalar couplings. We outline the chemical synthesis of the 6-(13)C-pyrimidine building blocks and their use in RNA solid phase synthesis and demonstrate their utility in Carr Purcell Meiboom Gill relaxation dispersion, ZZ exchange, and chemical exchange saturation transfer NMR experiments.

Structural Insights Reveal the Dynamics of the Repeating r(CAG) Transcript Found in Huntington's Disease (HD) and Spinocerebellar Ataxias (SCAs)

In humans, neurodegenerative disorders such as Huntington's disease (HD) and many spinocerebellar ataxias (SCAs) have been found to be associated with CAG trinucleotide repeat expansion. An important RNA-mediated mechanism that causes these diseases involves the binding of the splicing regulator protein MBNL1 (Muscleblind-like 1 protein) to expanded r(CAG) repeats. Moreover, mutant huntingtin protein translated from expanded r(CAG) also yields toxic effects. To discern the role of mutant RNA in these diseases, it is essential to gather information about its structure. Detailed insight into the different structures and conformations adopted by these mutant transcripts is vital for developing therapeutics targeting them. Here, we report the crystal structure of an RNA model with a r(CAG) motif, which is complemented by an NMR-based solution structure obtained from restrained Molecular Dynamics (rMD) simulation studies. Crystal structure data of the RNA model resolved at 2.3 Å reveals non-canonical pairing of adenine in 5´-CAG/3´-GAC motif samples in different syn and anti conformations. The overall RNA structure has helical parameters intermediate to the A- and B-forms of nucleic acids due to the global widening of major grooves and base-pair preferences near internal AA loops. The comprehension of structural behaviour by studying the spectral features and the dynamics also supports the flexible nature of the r(CAG) motif.

Thermodynamic examination of 1- to 5-nt purine bulge loops in RNA and DNA constructs

Bulge loops are common features of RNA structures that are involved in the formation of RNA tertiary structures and are often sites for interactions with proteins and ions. Minimal thermodynamic data currently exist on the bulge size and sequence effects. Using thermal denaturation methods, thermodynamic properties of 1- to 5-nt adenine and guanine bulge loop constructs were examined in 10 mM MgCl(2) or 1 M KCl. The [Formula: see text] loop parameters for 1- to 5-nt purine bulge loops in RNA constructs were between 3.07 and 5.31 kcal/mol in 1 M KCl buffer. In 10 mM magnesium ions, the ΔΔG° values relative to 1 M KCl were 0.47-2.06 kcal/mol more favorable for the RNA bulge loops. The [Formula: see text] loop parameters for 1- to 5-nt purine bulge loops in DNA constructs were between 4.54 and 5.89 kcal/mol. Only 4- and 5-nt guanine constructs showed significant change in stability for the DNA constructs in magnesium ions. A linear correlation is seen between the size of the bulge loop and its stability. New prediction models are proposed for 1- to 5-nt purine bulge loops in RNA and DNA in 1 M KCl. We show that a significant stabilization is seen for small bulge loops in RNA in the presence of magnesium ions. A prediction model is also proposed for 1- to 5-nt purine bulge loop RNA constructs in 10 mM magnesium chloride.

Slowdown of Interhelical Motions Induces a Glass Transition in RNA

RNA function depends crucially on the details of its dynamics. The simplest RNA dynamical unit is a two-way interhelical junction. Here, for such a unit--the transactivation response RNA element--we present evidence from molecular dynamics simulations, supported by nuclear magnetic resonance relaxation experiments, for a dynamical transition near 230 K. This glass transition arises from the freezing out of collective interhelical motional modes. The motions, resolved with site-specificity, are dynamically heterogeneous and exhibit non-Arrhenius relaxation. The microscopic origin of the glass transition is a low-dimensional, slow manifold consisting largely of the Euler angles describing interhelical reorientation. Principal component analysis over a range of temperatures covering the glass transition shows that the abrupt slowdown of motion finds its explanation in a localization transition that traps probability density into several disconnected conformational pools over the low-dimensional energy landscape. Upon temperature increase, the probability density pools then flood a larger basin, akin to a lakes-to-sea transition. Simulations on transactivation response RNA are also used to backcalculate inelastic neutron scattering data that match previous inelastic neutron scattering measurements on larger and more complex RNA structures and which, upon normalization, give temperature-dependent fluctuation profiles that overlap onto a glass transition curve that is quasi-universal over a range of systems and techniques.

Structural features of a 3' splice site in influenza a

Influenza A is an RNA virus with a genome of eight negative sense segments. Segment 7 mRNA contains a 3′ splice site for alternative splicing to encode the essential M2 protein. On the basis of sequence alignment and chemical mapping experiments, the secondary structure surrounding the 3′ splice site has an internal loop, adenine bulge, and hairpin loop when it is in the hairpin conformation that exposes the 3′ splice site. We report structural features of a three-dimensional model of the hairpin derived from nuclear magnetic resonance spectra and simulated annealing with restrained molecular dynamics. Additional insight was provided by modeling based on ¹H chemical shifts. The internal loop containing the 3′ splice site has a dynamic guanosine and a stable imino (cis Watson–Crick/Watson–Crick) GA pair. The adenine bulge also appears to be dynamic with the A either stacked in the stem or forming a base triple with a Watson–Crick GC pair. The hairpin loop is a GAAA tetraloop closed by an AC pair.

Computational docking simulations of a DNA-aptamer for argininamide and related ligands

The binding properties of sequence-specific nucleic acids (aptamers) to low-molecular-weight ligands, macromolecules and even cells attract substantial scientific interest. These ligand-DNA complexes found different applications for sensing, nanomedicine, and DNA nanotechnology. Structural information on the aptamer-ligand complexes is, however, scarce, even though it would open-up the possibilities to design novel features in the complexes. In the present study we apply molecular docking simulations to probe the features of an experimentally documented L-argininamide aptamer complex. The docking simulations were performed using AutoDock 4.0 and YASARA Structure software, a well-suited program for following intermolecular interactions and structures of biomolecules, including DNA. We explored the binding features of a DNA aptamer to L-argininamide and to a series of arginine derivatives or arginine-like ligands. We find that the best docking results are obtained after an energy-minimization of the parent ligand-aptamer complexes. The calculated binding energies of all mono-substituted guanidine-containing ligands show a good correlation with the experimentally determined binding constants. The results provide valuable guidelines for the application of docking simulations for the prediction of aptamer-ligand structures, and for the design of novel features of ligand-aptamer complexes.