RNA Conformational Changes Research Articles

Deciphering structure-activity relationships in a series of Tat/TAR inhibitors

A series of pentameric “Polyamide Amino Acids” (PAAs) compounds derived from the same trimeric precursor have been synthesized and investigated as HIV TAR RNA ligands, in the absence and in the presence of a Tat fragment. All PAAs bind TAR with similar sub-micromolar affinities but their ability to compete efficiently with the Tat fragment strongly differs, IC50 ranging from 35 nM to >2 μM. While NMR and CD studies reveal that all PAA interact with TAR at the same site and induce globally the same RNA conformational change upon binding, a comparative thermodynamic study of PAA/TAR equilibria highlights distinct TAR binding modes for Tat competitor and non-competitor PAAs. This led us to suggest two distinct interaction modes that have been further validated by molecular modeling studies. While the binding of Tat competitor PAAs induces a contraction at the TAR bulge region, the binding of non-competitor ones widens it. This could account for the distinct PAA ability to compete with Tat fragment. Our work illustrates how comparative thermodynamic studies of a series of RNA ligands of same chemical family are of value for understanding their binding modes and for rationalizing structure-activity relationships.

Traditional Chemical Mapping of RNA Structure In Vitro and In Vivo.

Chemical probing is often used to gain knowledge on the secondary and tertiary structures of RNA molecules either free or engaged in complexes with ligands. The method monitors the reactivity of each nucleotide towards chemicals of various specificities reflecting the hydrogen bonding environment of each nucleotide within the RNA molecule. In addition, information can be obtained on the binding site of a ligand (noncoding RNAs, protein, metabolites), and on RNA conformational changes that accompanied ligand binding or perturbation of the environmental cues. The detection of the modifications can be obtained either by using end-labeled RNA molecules or by primer extension using reverse transcriptase. The goal of this chapter is to provide the reader with an experimental guide to probe the structure of RNA in vitro and in vivo with the most suitable chemical probes.

Mg(2+) shifts ligand-mediated folding of a riboswitch from induced-fit to conformational selection.

Bacterial riboswitches couple small-molecule ligand binding to RNA conformational changes that widely regulate gene expression, rendering them potential targets for antibiotic intervention. Despite structural insights, the ligand-mediated folding mechanisms of riboswitches are still poorly understood. Using single-molecule fluorescence resonance energy transfer (smFRET), we have investigated the folding mechanism of an H-type pseudoknotted preQ1 riboswitch in dependence of Mg(2+) and three ligands of distinct affinities. We show that, in the absence of Mg(2+), both weakly and strongly bound ligands promote pseudoknot docking through an induced-fit mechanism. By contrast, addition of as low as 10 μM Mg(2+) generally shifts docking toward conformational selection by stabilizing a folded-like conformation prior to ligand binding. Supporting evidence from transition-state analysis further highlights the particular importance of stacking interactions during induced-fit and of specific hydrogen bonds during conformational selection. Our mechanistic dissection provides unprecedented insights into the intricate synergy between ligand- and Mg(2+)-mediated RNA folding.

Site-Directed Spin Labeling Evidence of the Leader-Linker Interaction in the Glycine Riboswitch using Electron Paramagnetic Resonance Spectroscopy

Riboswitches are mRNA transcripts that function to modulate genetic expression through selective recognition and binding of cognate ligands which, subsequently, induces RNA conformational changes. Continuous wave electron paramagnetic resonance (CW EPR) spectroscopy, when employed with site-directed spin labeling (SDSL), is a useful technique for investigating changes in site-specific dynamics within biological systems. In this work, SDSL CW EPR was used to study the leader-linker interaction in the Vibrio cholerae glycine riboswitch. The glycine riboswitch binds two molecules of glycine to regulate the expression of genes associated with glycine metabolism. The recently described leader-linker interaction in the glycine riboswitch has been investigated using biochemical methods and was shown to play a functional role in the ligand binding process. To probe local RNA backbone dynamics of select sites within the leader-linker interaction SDSL, using the R5 spin label, was employed. Incorporation of spin labels was achieved through the use of optimized ligation methodologies that allow small, synthetically modified RNA to be joined to the larger riboswitch RNA sequence. Empirical analysis of X-band EPR line shapes was used to characterize dynamics of the interaction at varying temperatures for differing folded states of the riboswitch in the absence or presence of salts and glycine ligand. Spectral variation at the labeled sites is in agreement with postulated secondary structural elements and has provided spectroscopic evidence in support of the leader-linker interaction.

Binding Induced RNA Conformational Changes Control Substrate Recognition and Catalysis by the Thiostrepton Resistance Methyltransferase (Tsr)

Ribosomal RNA (rRNA) post-transcriptional modifications are essential for ribosome maturation, translational fidelity, and are one mechanism used by both antibiotic-producing and pathogenic bacteria to resist the effects of antibiotics that target the ribosome. The thiostrepton producer Streptomyces azureus prevents self-intoxication by expressing the thiostrepton-resistance methyltransferase (Tsr), which methylates the 2'-hydroxyl of 23 S rRNA nucleotide adenosine 1067 within the thiostrepton binding site. Tsr is a homodimer with each protomer containing an L30e-like amino-terminal domain (NTD) and a SPOUT methyltransferase family catalytic carboxyl-terminal domain (CTD). We show that both enzyme domains are required for high affinity RNA substrate binding. The Tsr-CTD has intrinsic, weak RNA affinity that is necessary to direct the specific high-affinity Tsr-RNA interaction via NTDs, which have no detectable RNA affinity in isolation. RNA structure probing experiments identify the Tsr footprint on the RNA and structural changes in the substrate, induced specifically upon NTD binding, which are necessary for catalysis by the CTD. Additionally, we identify a key amino acid in each domain responsible for CTD-RNA binding and the observed NTD-dependent RNA structural changes. These studies allow us to develop a model for Tsr-RNA interaction in which the coordinated substrate recognition of each Tsr structural domain is an obligatory pre-catalytic recognition event. Our findings underscore the complexity of substrate recognition by RNA modification enzymes and the potential for direct involvement of the RNA substrate in controlling the process of its modification.

Ribosome RNA assembly intermediates visualized in living cells.

In cells, RNAs likely adopt numerous intermediate conformations prior to formation of functional RNA-protein complexes. We used single-nucleotide resolution selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) to probe the structure of Escherichia coli 16S rRNA in healthy growing bacteria. SHAPE-directed modeling indicated that the predominant steady-state RNA conformational ensemble in dividing cells had a base-paired structure different from that expected on the basis of comparative sequence analysis and high-resolution studies of the 30S ribosomal subunit. We identified the major cause of these differences by stopping ongoing in-cell transcription (in essence, an in-cell RNA structure pulse-chase experiment) which caused the RNA to chase into a structure that closely resembled the expected one. Most helices that formed alternate RNA conformations under growth conditions interact directly with tertiary-binding ribosomal proteins and form a C-shape that surrounds the mRNA channel and decoding site. These in-cell experiments lead to a model in which ribosome assembly factors function as molecular struts to preorganize this intermediate and emphasize that the final stages of ribonucleoprotein assembly involve extensive protein-facilitated RNA conformational changes.

Fluorescence methods in the investigation of the DEAD-box helicase mechanism.

DEAD-box proteins catalyze the ATP-dependent unwinding of RNA duplexes and accompany RNA molecules throughout their cellular life. Conformational changes in the helicase core of DEAD-box proteins are intimately linked to duplex unwinding. In the absence of ligands, the two RecA domains of the helicase core are separated. ATP and RNA binding induces a closure of the cleft between the RecA domains that is coupled to the distortion of bound RNA, leading to duplex destabilization and dissociation of one RNA strand. Reopening of the helicase core occurs after ATP hydrolysis and is coupled to phosphate release and dissociation of the second RNA strand.Fluorescence spectroscopy provides an array of approaches to study intermolecular interactions, local structural rearrangements, or large conformational changes of biomolecules. The fluorescence intensity of a fluorophore reports on its environment, and fluorescence anisotropy reflects the size of the molecular entity the fluorophore is part of. Fluorescence intensity and anisotropy are therefore sensitive probes to report on binding and dissociation events. Fluorescence resonance energy transfer (FRET) reports on the distance between two fluorophores and thus on conformational changes. Single-molecule FRET experiments reveal the distribution of conformational states and the kinetics of their interconversion. This chapter summarizes fluorescence approaches for monitoring individual aspects of DEAD-box protein activity, from nucleotide and RNA binding and RNA unwinding to protein and RNA conformational changes in the catalytic cycle, and illustrates exemplarily how fluorescence-based methods have contributed to understanding the mechanism of DEAD-box helicase-catalyzed RNA unwinding.

Uncoating of common cold virus is preceded by RNA switching as determined by X-ray and cryo-EM analyses of the subviral A-particle

During infection, viruses undergo conformational changes that lead to delivery of their genome into host cytosol. In human rhinovirus A2, this conversion is triggered by exposure to acid pH in the endosome. The first subviral intermediate, the A-particle, is expanded and has lost the internal viral protein 4 (VP4), but retains its RNA genome. The nucleic acid is subsequently released, presumably through one of the large pores that open at the icosahedral twofold axes, and is transferred along a conduit in the endosomal membrane; the remaining empty capsids, termed B-particles, are shuttled to lysosomes for degradation. Previous structural analyses revealed important differences between the native protein shell and the empty capsid. Nonetheless, little is known of A-particle architecture or conformation of the RNA core. Using 3D cryo-electron microscopy and X-ray crystallography, we found notable changes in RNA-protein contacts during conversion of native virus into the A-particle uncoating intermediate. In the native virion, we confirmed interaction of nucleotide(s) with Trp(38) of VP2 and identified additional contacts with the VP1 N terminus. Study of A-particle structure showed that the VP2 contact is maintained, that VP1 interactions are lost after exit of the VP1 N-terminal extension, and that the RNA also interacts with residues of the VP3 N terminus at the fivefold axis. These associations lead to formation of a well-ordered RNA layer beneath the protein shell, suggesting that these interactions guide ordered RNA egress.

Deriving quantitative dynamics information for proteins and RNAs using ROTDIF with a graphical user interface

To facilitate rigorous analysis of molecular motions in proteins, DNA, and RNA, we present a new version of ROTDIF, a program for determining the overall rotational diffusion tensor from single- or multiple-field nuclear magnetic resonance relaxation data. We introduce four major features that expand the program's versatility and usability. The first feature is the ability to analyze, separately or together, (13)C and/or (15)N relaxation data collected at a single or multiple fields. A significant improvement in the accuracy compared to direct analysis of R2/R1 ratios, especially critical for analysis of (13)C relaxation data, is achieved by subtracting high-frequency contributions to relaxation rates. The second new feature is an improved method for computing the rotational diffusion tensor in the presence of biased errors, such as large conformational exchange contributions, that significantly enhances the accuracy of the computation. The third new feature is the integration of the domain alignment and docking module for relaxation-based structure determination of multi-domain systems. Finally, to improve accessibility to all the program features, we introduced a graphical user interface that simplifies and speeds up the analysis of the data. Written in Java, the new ROTDIF can run on virtually any computer platform. In addition, the new ROTDIF achieves an order of magnitude speedup over the previous version by implementing a more efficient deterministic minimization algorithm. We not only demonstrate the improvement in accuracy and speed of the new algorithm for synthetic and experimental (13)C and (15)N relaxation data for several proteins and nucleic acids, but also show that careful analysis required especially for characterizing RNA dynamics allowed us to uncover subtle conformational changes in RNA as a function of temperature that were opaque to previous analysis.

Identification of the binding site of Rlp7 on assembling 60S ribosomal subunits in Saccharomyces cerevisiae

Eukaryotic ribosome assembly requires over 200 assembly factors that facilitate rRNA folding, ribosomal protein binding, and pre-rRNA processing. One such factor is Rlp7, an essential RNA binding protein required for consecutive pre-rRNA processing steps for assembly of yeast 60S ribosomal subunits: exonucleolytic processing of 27SA3 pre-rRNA to generate the 5' end of 5.8S rRNA and endonucleolytic cleavage of the 27SB pre-rRNA to initiate removal of internal transcribed spacer 2 (ITS2). To better understand the functions of Rlp7 in 27S pre-rRNA processing steps, we identified where it crosslinks to pre-rRNA. We found that Rlp7 binds at the junction of ITS2 and the ITS2-proximal stem, between the 3' end of 5.8S rRNA and the 5' end of 25S rRNA. Consistent with Rlp7 binding to this neighborhood during assembly, two-hybrid and affinity copurification assays showed that Rlp7 interacts with other assembly factors that bind to or near ITS2 and the proximal stem. We used in vivo RNA structure probing to demonstrate that the proximal stem forms prior to Rlp7 binding and that Rlp7 binding induces RNA conformational changes in ITS2 that may chaperone rRNA folding and regulate 27S pre-rRNA processing. Our findings contradict the hypothesis that Rlp7 functions as a placeholder for ribosomal protein L7, from which Rlp7 is thought to have evolved in yeast. The binding site of Rlp7 is within eukaryotic-specific RNA elements, which are not found in bacteria. Thus, we propose that Rlp7 coevolved with these RNA elements to facilitate eukaryotic-specific functions in ribosome assembly and pre-rRNA processing.

Thermodynamic studies of a series of homologous HIV-1 TAR RNA ligands reveal that loose binders are stronger Tat competitors than tight ones

RNA is a major drug target, but the design of small molecules that modulate RNA function remains a great challenge. In this context, a series of structurally homologous ‘polyamide amino acids’ (PAA) was studied as HIV-1 trans-activating response (TAR) RNA ligands. An extensive thermodynamic study revealed the occurence of an enthalpy–entropy compensation phenomenon resulting in very close TAR affinities for all PAA. However, their binding modes and their ability to compete with the Tat fragment strongly differ according to their structure. Surprisingly, PAA that form loose complexes with TAR were shown to be stronger Tat competitors than those forming tight ones, and thermal denaturation studies demonstrated that loose complexes are more stable than tight ones. This could be correlated to the fact that loose and tight ligands induce distinct RNA conformational changes as revealed by circular dichroism experiments, although nuclear magnetic resonance (NMR) experiments showed that the TAR binding site is the same in all cases. Finally, some loose PAA also display promising inhibitory activities on HIV-infected cells. Altogether, these results lead to a better understanding of RNA interaction modes that could be very useful for devising new ligands of relevant RNA targets.

Specificity of the Double-Stranded RNA-Binding Domain from the RNA-Activated Protein Kinase PKR for Double-Stranded RNA: Insights from Thermodynamics and Small-Angle X-ray Scattering

The interferon-inducible, double-stranded (ds) RNA-activated protein kinase (PKR) contains a dsRNA-binding domain (dsRBD) and plays key roles in viral pathogenesis and innate immunity. Activation of PKR is typically mediated by long dsRNA, and regulation of PKR is disfavored by most RNA imperfections, including bulges and internal loops. Herein, we combine isothermal titration calorimetry (ITC), electrophoretic mobility shift assays, and small-angle X-ray scattering (SAXS) to dissect the thermodynamic basis for the specificity of the dsRBD termed "p20" for various RNAs and to detect any RNA conformational changes induced upon protein binding. We monitor binding of p20 to chimeric duplexes containing terminal RNA-DNA hybrid segments and a central dsRNA segment, which was either unbulged ("perfect") or bulged. The ITC data reveal strong binding of p20 to the perfect duplex (K(d) ~ 30 nM) and weaker binding to the bulged duplex (K(d) ~ 2-5 μM). SAXS reconstructions and p(r) distance distribution functions further uncover that p20 induces no significant conformational change in perfect dsRNA but largely straightens bulged dsRNA. Together, these observations support the dsRBD's ability to tightly bind to only A-form RNA and suggest that in a noninfected cell, PKR may be buffered via weak interactions with various bulged and looped RNAs, which it may straighten. This work suggests that PKR-regulating RNAs with complex secondary and tertiary structures likely mimic dsRNA and/or engage portions of PKR outside of the dsRBD.

Visualizing transient low-populated structures of RNA

The visualization of RNA conformational changes has provided fundamental insights into how regulatory RNAs carry out their biological functions. The RNA structural transitions that have been characterized to date involve long-lived species that can be captured by structure characterization techniques. Here, we report the Nuclear Magnetic Resonance visualization of RNA transitions towards invisible ‘excited states’ (ES), which exist in too little abundance (2–13%) and for too short periods of time (45–250 μs) to allow structural characterization by conventional techniques. Transitions towards ESs result in localized rearrangements in base-pairing that alter building block elements of RNA architecture, including helix-junction-helix motifs and apical loops. The ES can inhibit function by sequestering residues involved in recognition and signaling or promote ATP-independent strand exchange. Thus, RNAs do not adopt a single conformation, but rather exist in rapid equilibrium with alternative ESs, which can be stabilized by cellular cues to affect functional outcomes.

Rational design of modular allosteric aptamer sensor for label-free protein detection

An aptamer can be redesigned to new functional molecules by conjugating with other oligonucleotides. However, it requires experimental trials to optimize the conjugating module with the sensitivity and selectivity toward a target. To reduce these efforts, we report rationally-designed modular allosteric aptamer sensor (MAAS), which is composed of coupled two aptamers and the regulator. For label-free protein detection, the protein–aptamer was conjugated with the malachite green (MG) aptamer for signaling. The MAAS additionally has the regulator domain which is designed to hybridize to a protein binding domain. The regulator makes MAAS to be inactive by destructing the original structure of the two aptamers. However, its conformation becomes active by dissociating the hybridization from the protein recognition signal, thereby inducing the binding of MG emitting the enhanced fluorescence. The design of regulator is based on the thermodynamic energy difference by the RNA conformational change and protein–aptamer affinity. Here we first demonstrated the MAAS for hepatitis C helicase and replicase. The target proteins were detected up to 250nM with minimized blank signals and displayed high specificities 10-fold greater than in non-specific proteins. The MAAS provides valuable tools that can be adapted to a wide range of configurations in bioanalytical applications.

Sub-Millisecond RNA Collapse Observed in a Microfluidic Mixer

Measuring the fast collapse of RNA molecules can elucidate the first steps of RNA folding. The P5abc sub-domain of the Tetrahymena ribozyme is an example of an RNA molecule that folds rapidly following the addition of Mg2+ ions. We studied changes in truncated-P5abc (tP5abc) RNA as reported by Forster Resonance Energy Transfer (FRET) in a microfluidic mixing device using confocal microscopy. With sub-millisecond time resolution, we measured the ion species and valence dependent collapse timescales of tP5abc. In this presentation, we discuss differences between collapse times of tP5abc and those of a mutant which does not fully fold in the presence of magnesium ions. We also discuss experimental issues that arise when using an internally labeled fluorophore to report changes in RNA conformation. To account for changes in the FRET efficiency we consider results of fluorescence anisotropy measurements, as well as fluorescence correlation spectroscopy measurements.

MBNL1-RNA Interactions: Binding-Induced Rna Conformational Changes

The muscleblind-like protein 1 (MBNL1) is an RNA binding protein that regulates mRNA alternative splicing. It is implicated to be involved in myotonic dystrophy, an inherited multisystem disease, through interacting with poly(CUG) or poly(CCUG) RNA. In the disease state, MBNL1 is abnormally expressed and translocated, which leads to many mis-splicing events. MBNL1 can bind various RNA substrates, most of which share a stem-loop structure containing 5’YGCY3’(Y=U or C) motifs. Using FRET together with a combination of stopped-flow kinetic fluorescence measurements and gel shift assays, we have investigated the relationship between MBNL1-RNA interactions and RNA secondary structure. This work suggested that binding of MBNL1 destabilizes the RNA secondary structure and leads to a partial melting of the stem loop RNA. The RNA secondary structure stability is directly correlated with the kinetic and thermodynamic properties of the MBNL1-RNA complex. The data presented here provides insight into the role of RNA secondary structure in regulation of alternative splicing by MBNL1 and suggests new perspectives for drug development for myotonic dystrophy.

Toward a Digital Gene Response: RNA G-Quadruplexes with Fewer Quartets Fold with Higher Cooperativity

Changes in RNA conformation can alter gene expression. The guanine quadruplex sequence (GQS) is an RNA motif that folds in the presence of K(+) ions. Changes in the conformation of this motif could be especially important in regulating gene expression in plants because intracellular K(+) concentrations often increase during drought stress. Little is known about the folding thermodynamics of RNA GQS. We show here that RNA GQS with tracts containing three G's [e.g., (GGGxx)(4)] have a modest dependence on the K(+) concentration, folding with no or even negative cooperativity (Hill coefficients ≤1), and are associated with populated folding intermediates. In contrast, GQS with tracts containing just two G's [e.g., (GGxx)(4)] have a steep dependence on the K(+) concentration and fold with positive cooperativity (Hill coefficients of 1.7-2.7) without significantly populating intermediate states. We postulate that in plants, the more stable G3 sequences are largely folded even under unstressed conditions, while the less stable G2 sequences fold only at the higher K(+) concentrations associated with cellular stress, wherein they respond sharply to changing K(+) concentrations. Given the binary nature of their folding, G2 sequences may find application in computation with DNA and in engineering of genetic circuits.

A template trigger for translocation

The telomerase ribonucleoprotein enzyme maintains eukaryotic chromosomes through unique means. Its telomerase reverse transcriptase (TERT) enzymatic subunit synthesizes telomeric DNA repeats from the template RNA subunit (TER) by 'repeat addition processivity'. How telomerase recycles this internal template, and how the template progresses through the TERT active site, was unclear. Berman et al. addressed these questions in Tetrahymena thermophila. They used biochemical assays to assess which regions of TER are important for enzyme activity and assayed changes in RNA conformation during the enzymatic cycle using single molecule fluorescence resonance energy transfer. They propose that the RNA template actively drives translocation through the TERT active site by an 'accordion' movement, in which the RNA sequences that flank the template ensure its correct position and translocation.

Molecular Mechanics Investigation of an Adenine–Adenine Non-Canonical Pair Conformational Change

Conformational changes are important in RNA for binding and catalysis and understanding these changes is important for understanding how RNA functions. Computational techniques using all-atom molecular models can be used to characterize conformational changes in RNA. These techniques are applied to an RNA conformational change involving a single base pair within a nine base pair RNA duplex. The Adenine-Adenine (AA) non-canonical pair in the sequence 5'GGUGAAGGCU3' paired with 3'PCCGAAGCCG5', where P is Purine, undergoes conformational exchange between two conformations on the timescale of tens of microseconds, as demonstrated in a previous NMR solution structure [Chen, G., et al., Biochemistry, 2006. 45: 6889-903]. The more populated, major, conformation was estimated to be 0.5 to 1.3 kcal/mol more stable at 30 °C than the less populated, minor, conformation. Both conformations are trans-Hoogsteen/sugar edge pairs, where the interacting edges on the adenines change with the conformational change. Targeted Molecular Dynamics (TMD) and Nudged Elastic Band (NEB) were used to model the pathway between the major and minor conformations using the AMBER software package. The adenines were predicted to change conformation via intermediates in which they are stacked as opposed to hydrogen-bonded. The predicted pathways can be described by an improper dihedral angle reaction coordinate. Umbrella sampling along the reaction coordinate was performed to model the free energy profile for the conformational change using a total of 1800 ns of sampling. Although the barrier height between the major and minor conformations was reasonable, the free energy difference between the major and minor conformations was the opposite of that expected based on the NMR experiments. Variations in the force field applied did not improve the misrepresentation of the free energies of the major and minor conformations. As an alternative, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approximation was applied to predict free energy differences between the two conformations using a total of 800 ns of sampling. MM-PBSA also incorrectly predicted the major conformation to be higher in free energy than the minor conformation.

Assembly of Bacterial Ribosomes

The assembly of ribosomes from a discrete set of components is a key aspect of the highly coordinated process of ribosome biogenesis. In this review, we present a brief history of the early work on ribosome assembly in Escherichia coli, including a description of in vivo and in vitro intermediates. The assembly process is believed to progress through an alternating series of RNA conformational changes and protein-binding events; we explore the effects of ribosomal proteins in driving these events. Ribosome assembly in vivo proceeds much faster than in vitro, and we outline the contributions of several of the assembly cofactors involved, including Era, RbfA, RimJ, RimM, RimP, and RsgA, which associate with the 30S subunit, and CsdA, DbpA, Der, and SrmB, which associate with the 50S subunit.