Towards structural genomics of RNA: rapid NMR resonance assignment and simultaneous RNA tertiary structure determination using residual dipolar couplings

We report a new residual dipolar couplings (RDCs) based NMR procedure for rapidly determining RNA tertiary structure demonstrated on a uniformly (15)N/(13)C-labeled 27 nt variant of the trans-activation response element (TAR) RNA from HIV-I. In this procedure, the time-consuming nuclear Overhauser enhancement (NOE)-based sequential assignment step is replaced by a fully automated RDC-based assignment strategy. This approach involves examination of all allowed sequence-specific resonance assignment permutations for best-fit agreement between measured RDCs and coordinates for sub-structures in a target RNA. Using idealized A-form geometries to model Watson-Crick helices and coordinates from a previous X-ray structure to model a hairpin loop in TAR, the best-fit RDC assignment solutions are determined very rapidly (相似文献   

Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.     

Mg2+-induced variations in the conformation and dynamics of HIV-1 TAR RNA probed using NMR residual dipolar couplings

The effects of divalent Mg(2+) on the conformation and dynamics of the stem-loop transactivation response element (TAR) RNA from HIV-1 have been characterized using NMR residual dipolar couplings (RDCs). Order matrix analysis of one bond 13C-1H RDCs measured in TAR at [Mg(2+)]:[TAR] stoichiometric ratios of approximately 3:1 (TAR(3.0Mg)) and approximately 4.5:1 (TAR(4.5Mg)) revealed that Mg(2+) reduces the average inter-helical angle from 47(+/-5) degrees in TAR(free) to 5(+/-7) degrees in TAR(4.5Mg). In contrast to the TAR(free) state, the generalized degree of order for the two stems in TAR(4.5Mg) is found to be identical within experimental uncertainty, indicating that binding of Mg(2+) leads to an arrest of inter-helical motions in TAR(free). Results demonstrate that RDC-NMR methodology can provide new insight into the effects of Mg(2+) on both the conformation and dynamics of RNA.     

Argininamide binding arrests global motions in HIV-1 TAR RNA: comparison with Mg2+-induced conformational stabilization

The structure and dynamics of the stem-loop transactivation response element (TAR) RNA from the human immunodeficiency virus type-1 (HIV-1) bound to the ligand argininamide (ARG) has been characterized using a combination of a large number of residual dipolar couplings (RDCs) and trans-hydrogen bond NMR methodology. Binding of ARG to TAR changes the average inter-helical angle between the two stems from approximately 47 degrees in the free state to approximately 11 degrees in the bound state, and leads to the arrest of large amplitude (+/-46 degrees ) inter-helical motions observed previously in the free state. While the global structural dynamics of TAR-ARG is similar to that previously reported for TAR bound to Mg2+, there are substantial differences in the hydrogen bond alignment of bulge and neighboring residues. Based on a novel H5(C5)NN experiment for probing hydrogen-mediated 2hJ(N,N) scalar couplings as well as measured RDCs, the TAR-ARG complex is stabilized by a U38-A27.U23 base-triple involving an A27.U23 reverse Hoogsteen hydrogen bond alignment as well as by a A22-U40 Watson-Crick base-pair at the junction of stem I. These hydrogen bond alignments are not observed in either the free or Mg2+ bound forms of TAR. The combined conformational analysis of TAR under three states reveals that ligands and divalent ions can stabilize similar RNA global conformations through distinct interactions involving different hydrogen bond alignments in the RNA.     

Probing Na(+)-induced changes in the HIV-1 TAR conformational dynamics using NMR residual dipolar couplings: new insights into the role of counterions and electrostatic interactions in adaptive recognition

Many regulatory RNAs undergo large changes in structure upon recognition of proteins and ligands, but the mechanism by which this occurs remains poorly understood. Using NMR residual dipolar coupling (RDCs), we characterized Na+-induced changes in the structure and dynamics of the bulge-containing HIV-1 transactivation response element (TAR) RNA that mirrors changes induced by small molecules bearing a different number of cationic groups. Increasing the Na+ concentration from 25 to 320 mM led to a continuous reduction in the average inter-helical bend angle (from 46 degrees to 22 degrees ), inter-helical twist angle (from 66 degrees to -18 degrees ), and inter-helix flexibility (as measured by an increase in the internal generalized degree of order from 0.56 to 0.74). Similar conformational changes were observed with Mg2+, indicating that nonspecific electrostatic interactions drive the conformational transition, although results also suggest that Na+ and Mg2+ may associate with TAR in distinct modes. The transition can be rationalized on the basis of a population-weighted average of two ensembles comprising an electrostatically relaxed bent and flexible TAR conformation that is weakly associated with counterions and a globally rigid coaxial conformation that has stronger electrostatic potential and association with counterions. The TAR inter-helical orientations that are stabilized by small molecules fall around the metal-induced conformational pathway, indicating that counterions may help predispose the TAR conformation for target recognition. Our results underscore the intricate sensitivity of RNA conformational dynamics to environmental conditions and demonstrate the ability to detect subtle conformational changes using NMR RDCs.     

Discovery of selective bioactive small molecules by targeting an RNA dynamic ensemble

Current approaches used to identify protein-binding small molecules are not suited for identifying small molecules that can bind emerging RNA drug targets. By docking small molecules onto an RNA dynamic ensemble constructed by combining NMR spectroscopy and computational molecular dynamics, we virtually screened small molecules that target the entire structure landscape of the transactivation response element (TAR) from HIV type 1 (HIV-1). We quantitatively predict binding energies for small molecules that bind different RNA conformations and report the de novo discovery of six compounds that bind TAR with high affinity and inhibit its interaction with a Tat peptide in vitro (K(i) values of 710 nM-169 μM). One compound binds HIV-1 TAR with marked selectivity and inhibits Tat-mediated activation of the HIV-1 long terminal repeat by 81% in T-cell lines and HIV replication in an HIV-1 indicator cell line (IC(50) ～23.1 μM).     

Impact of static and dynamic A-form heterogeneity on the determination of RNA global structural dynamics using NMR residual dipolar couplings

We examined how static and dynamic deviations from the idealized A-form helix propagate into errors in the principal order tensor parameters determined using residual dipolar couplings (rdcs). A 20-ns molecular dynamics (MD) simulation of the HIV-1 transactivation response element (TAR) RNA together with a survey of spin relaxation studies of RNA dynamics reveals that pico-to-nanosecond local motions in non-terminal Watson-Crick base-pairs will uniformly attenuate base and sugar one bond rdcs by approximately 7%. Gaussian distributions were generated for base and sugar torsion angles through statistical comparison of 40 RNA X-ray structures solved to <3.0 A resolution. For a typical number (>or=11) of one bond C-H base and sugar rdcs, these structural deviations together with rdc uncertainty (1.5 Hz) lead to average errors in the magnitude and orientation of the principal axis of order that are <9% and <4 degrees, respectively. The errors decrease to <5% and <4 degrees for >or=17 rdcs. A protocol that allows for estimation of error in A-form order tensors due to both angular deviations and rdc uncertainty (Aform-RDC) is validated using theoretical simulations and used to analyze rdcs measured previously in TAR in the free state and bound to four distinct ligands. Results confirm earlier findings that the two TAR helices undergo large changes in both their mean relative orientation and dynamics upon binding to different targets.     

Dynamic ensemble view of the conformational landscape of HIV-1 TAR RNA and allosteric recognition

RNA conformational dynamics and the resulting structural heterogeneity play an important role in RNA functions, e.g., recognition. Recognition of HIV-1 TAR RNA has been proposed to occur via a conformational capture mechanism. Here, using ultrafast time-resolved fluorescence spectroscopy, we have probed the complexity of the conformational landscape of HIV-1 TAR RNA and monitored the position-dependent changes in the landscape upon binding of a Tat protein-derived peptide and neomycin B. In the ligand-free state, the TAR RNA samples multiple families of conformations with various degrees of base stacking around the three-nucleotide bulge region. Some subpopulations partially resemble those ligand-bound states, but the coaxially stacked state is below the detection limit. When Tat or neomycin B binds, the bulge region as an ensemble undergoes a conformational transition in a position-dependent manner. Tat and neomycin B induce mutually exclusive changes in the TAR RNA underlying the mechanism of allosteric inhibition at an ensemble level with residue-specific details. Time-resolved anisotropy decay measurements revealed picosecond motions of bases in both ligand-free and ligand-bound states. Mutation of a base pair at the bulge--stem junction has differential effects on the conformational distributions of the bulge bases. A dynamic model of the ensemble view of the conformational landscape for HIV-1 TAR RNA is proposed, and the implication of the general mechanism of RNA recognition and its impact on RNA-based therapeutics are discussed.     

Structural features in TAR RNA of human and simian immunodeficiency viruses: a phylogenetic analysis.

A comparative analysis of TAR RNA structures in human and simian immunodeficiency viruses reveals the conservation of certain structural features despite the divergence in sequence. Both the TAR elements of HIV-1 and SIV-chimpanzee can be folded into relatively simple one-stem hairpin structures. Chemical and RNAase probes were used to analyze the more complex structure of HIV-2 TAR RNA, which folds into a branched hairpin structure. A surprisingly similar RNA conformation can be proposed for SIV-mandrill, despite considerable divergence in nucleotide sequence. A third structural presentation of TAR sequences is seen for SIV-african green monkey. These results are generally consistent with the classification of HIV-SIV viruses in four subgroups based on sequence analyses (both nucleotide- and amino acid-sequences). However, some conserved TAR structures were detected for members of different virus subgroups. It is therefore proposed that RNA structure analysis might provide an additional tool for determining phylogenetic relationships among the HIV-SIV viruses.     

NMR studies of RNA dynamics and structural plasticity using NMR residual dipolar couplings

An increasing number of RNAs are being discovered that perform their functions by undergoing large changes in conformation in response to a variety of cellular signals, including recognition of proteins and small molecular targets, changes in temperature, and RNA synthesis itself. The measurement of NMR residual dipolar couplings (RDCs) in partially aligned systems is providing new insights into the structural plasticity of RNA through combined characterization of large-amplitude collective helix motions and local flexibility in noncanonical regions over a wide window of biologically relevant timescales (相似文献   

Coupled nucleotide covariations reveal dynamic RNA interaction patterns

Evolutionarily conserved structures in related RNA molecules contain coordinated variations (covariations) of paired nucleotides. Analysis of covariations is a very powerful approach to deduce phylogenetically conserved (i.e., functional) conformations, including tertiary interactions. Here we discuss conserved RNA folding pathways that are revealed by covariation patterns. In such pathways, structural requirements for alternative pairings cause some nucleotides to covary with two different partners. Such "coupled" covariations between three or more nucleotides were found in various types of RNAs. The analysis of coupled covariations can unravel important features of RNA folding dynamics and improve phylogeny reconstruction in some cases. Importantly, it is necessary to distinguish between multiple covariations determined by mutually exclusive structures and those determined by tertiary contacts.     

Improved prediction of RNA tertiary structure with insights into native state dynamics

The importance of RNA tertiary structure is evident from the growing number of published high resolution NMR and X-ray crystallographic structures of RNA molecules. These structures provide insights into function and create a knowledge base that is leveraged by programs such as Assemble, ModeRNA, RNABuilder, NAST, FARNA, Mc-Sym, RNA2D3D, and iFoldRNA for tertiary structure prediction and design. While these methods sample native-like RNA structures during simulations, all struggle to capture the native RNA conformation after scoring. We propose RSIM, an improved RNA fragment assembly method that preserves RNA global secondary structure while sampling conformations. This approach enhances the quality of predicted RNA tertiary structure, provides insights into the native state dynamics, and generates a powerful visualization of the RNA conformational space. RSIM is available for download from http://www.github.com/jpbida/rsim.     

Efficient in silico exploration of RNA interhelical conformations using Euler angles and WExplore

HIV-1 TAR RNA is a two-helix bulge motif that plays a critical role in HIV viral replication and is an important drug target. However, efforts at designing TAR inhibitors have been challenged by its high degree of structural flexibility, which includes slow large-amplitude reorientations of its helices with respect to one another. Here, we use the recently introduced algorithm WExplore in combination with Euler angles to achieve unprecedented sampling of the TAR conformational ensemble. Our ensemble achieves similar agreement with experimental NMR data when compared with previous TAR computational studies, and is generated at a fraction of the computational cost. It clearly emerges from configuration space network analysis that the intermittent formation of the A22-U40 base pair acts as a reversible switch that enables sampling of interhelical conformations that would otherwise be topologically disallowed. We find that most previously determined ligand-bound structures are found in similar location in the network, and we use a sample-and-select approach to guide the construction of a set of novel conformations which can serve as the basis for future drug development efforts. Collectively, our findings demonstrate the utility of WExplore in combination with suitable order parameters as a method for exploring RNA conformational space.     

A conditional random fields method for RNA sequence-structure relationship modeling and conformation sampling

Accurate tertiary structures are very important for the functional study of non-coding RNA molecules. However, predicting RNA tertiary structures is extremely challenging, because of a large conformation space to be explored and lack of an accurate scoring function differentiating the native structure from decoys. The fragment-based conformation sampling method (e.g. FARNA) bears shortcomings that the limited size of a fragment library makes it infeasible to represent all possible conformations well. A recent dynamic Bayesian network method, BARNACLE, overcomes the issue of fragment assembly. In addition, neither of these methods makes use of sequence information in sampling conformations. Here, we present a new probabilistic graphical model, conditional random fields (CRFs), to model RNA sequence-structure relationship, which enables us to accurately estimate the probability of an RNA conformation from sequence. Coupled with a novel tree-guided sampling scheme, our CRF model is then applied to RNA conformation sampling. Experimental results show that our CRF method can model RNA sequence-structure relationship well and sequence information is important for conformation sampling. Our method, named as TreeFolder, generates a much higher percentage of native-like decoys than FARNA and BARNACLE, although we use the same simple energy function as BARNACLE. CONTACT: zywang@ttic.edu; j3xu@ttic.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.