Long-range Tertiary Interactions Research Articles

A comprehensive survey of long-range tertiary interactions and motifs in non-coding RNA structures.

Understanding the 3D structure of RNA is key to understanding RNA function. RNA 3D structure is modular and can be seen as a composition of building blocks of various sizes called tertiary motifs. Currently, long-range motifs formed between distant loops and helical regions are largely less studied than the local motifs determined by the RNA secondary structure. We surveyed long-range tertiary interactions and motifs in a non-redundant set of non-coding RNA 3D structures. A new dataset of annotated LOng-RAnge RNA 3D modules (LORA) was built using an approach that does not rely on the automatic annotations of non-canonical interactions. An original algorithm, ARTEM, was developed for annotation-, sequence- and topology-independent superposition of two arbitrary RNA 3D modules. The proposed methods allowed us to identify and describe the most common long-range RNA tertiary motifs. Along with the prevalent canonical A-minor interactions, a large numberof previously undescribed staple interactions were observed. The most frequent long-range motifs were found to belong to three main motif families: planar staples, tilted staples, and helical packing motifs.

RNAcanvas: interactive drawing and exploration of nucleic acid structures.

Two-dimensional drawing of nucleic acid structures, particularly RNA structures, is fundamental to the communication of nucleic acids research. However, manually drawing structures is laborious and infeasible for structures thousands of nucleotides long. RNAcanvas automatically arranges residues into strictly shaped stems and loops while providing robust interactive editing features, including click-and-drag layout adjustment. Drawn elements are highly customizable in a point-and-click manner, including colours, fonts, size and shading, flexible numbering, and outlining of bases. Tertiary interactions can be drawn as draggable, curved lines. Leontis-Westhof notation for depicting non-canonical base-pairs is fully supported, as well as text labels for structural features (e.g. hairpins). RNAcanvas also has many unique features and performance optimizations for large structures that cannot be correctly predicted and require manual refinement based on the researcher's own analyses and expertise. To this end, RNAcanvas has point-and-click structure editing with real-time highlighting of complementary sequences and motif search functionality, novel features that greatly aid in the identification of putative long-range tertiary interactions, de novo analysis of local structures, and phylogenetic comparisons. For ease in producing publication quality figures, drawings can be exported in both SVG and PowerPoint formats. URL: https://rnacanvas.app.

Context-dependence of T-loop Mediated Long-range RNA Tertiary Interactions

The architecture and folding of complex RNAs is governed by a limited set of highly recurrent structural motifs that form long-range tertiary interactions. One of these motifs is the T-loop, which was first identified in tRNA but is broadly distributed across biological RNAs. While the T-loop has been examined in detail in different biological contexts, the various receptors that it interacts with are not as well defined. In this study, we use a cell-based genetic screen in concert with bioinformatic analysis to examine three different, but related, T-loop receptor motifs found in the flavin mononucleotide (FMN) and cobalamin (Cbl) riboswitches. As a host for different T-loop receptors, we employed the env8 class-II Cbl riboswitch, an RNA that uses two T-loop motifs for both folding and supporting the ligand binding pocket. A set of libraries was created in which select nucleotides that participate in the T-loop/T-loop receptor (TL/TLR) interaction were fully randomized. Library members were screened for their ability to support Cbl-dependent expression of a reporter gene. While T-loops appear to be variable in sequence, we find that the functional sequence space is more restricted in the Cbl riboswitch, suggesting that TL/TLR interactions are context dependent. Our data reveal clear sequence signatures for the different types of receptor motifs that align with phylogenic analysis of these motifs in the FMN and Cbl riboswitches. Finally, our data suggest the functional contribution of various nucleobase-mediated long-range interactions within the riboswitch subclass of TL/TLR interactions that are distinct from those found in other RNAs.

Orchestrating ribosomal RNA folding during ribosome assembly.

Construction of the eukaryotic ribosome is a complex process in which a nascent ribosomal RNA (rRNA) emerging from RNA Polymerase I hierarchically folds into a native three-dimensional structure. Modular assembly of individual RNA domains through interactions with ribosomal proteins and a myriad of assembly factors permit efficient disentanglement of the error-prone RNA folding process. Following these dynamic events, long-range tertiary interactions are orchestrated to compact rRNA. A combination of genetic, biochemical, and structural studies is now providing clues into how a nascent rRNA is transformed into a functional ribosome with high precision. With this essay, we aim to draw attention to the poorly understood process of establishing correct RNA tertiary contacts during ribosome formation.

IsRNA1: De Novo Prediction and Blind Screening of RNA 3D Structures.

Modeling structures and functions of large ribonucleic acid (RNAs) especially with complicated topologies is highly challenging due to the inefficiency of large conformational sampling and the presence of complicated tertiary interactions. To address this problem, one highly promising approach is coarse-grained modeling. Here, following an iterative simulated reference state approach to decipher the correlations between different structural parameters, we developed a potent coarse-grained RNA model named as IsRNA1 for RNA studies. Molecular dynamics simulations in the IsRNA1 can predict the native structures of small RNAs from a sequence and fold medium-sized RNAs into near-native tertiary structures with the assistance of secondary structure constraints. A large-scale benchmark test on RNA 3D structure prediction shows that IsRNA1 exhibits improved performance for relatively large RNAs of complicated topologies, such as large stem-loop structures and structures containing long-range tertiary interactions. The advantages of IsRNA1 include the consideration of the correlations between the different structural variables, the appropriate characterization of canonical base-pairing and base-stacking interactions, and the better sampling for the backbone conformations. Moreover, a blind screening protocol was developed based on IsRNA1 to identify good structural models from a pool of candidates without prior knowledge of the native structures.

Time-Resolved, Single-Molecule, Correlated Chemical Probing of RNA.

Capturing the folding dynamics of large, functionally important RNAs has relied primarily on global measurements of structure or on per-nucleotide chemical probing. These approaches infer, but do not directly measure, through-space structural interactions. Here we introduce trimethyloxonium (TMO) as a chemical probe for RNA. TMO alkylates RNA at high levels in seconds, and thereby enables time-resolved, single-molecule, through-space probing of RNA folding using the RING-MaP correlated chemical probing framework. Time-resolved correlations in the RNase P RNA-a functional RNA with a complex structure stabilized by multiple noncanonical interactions-revealed that a long-range tertiary interaction guides native RNA folding for both secondary and tertiary structure. This unanticipated nonhierarchical folding mechanism was directly validated by examining the consequences of concise disruption of the through-space interaction. Single-molecule, time-resolved RNA structure probing with TMO is poised to reveal a wide range of dynamic RNA folding processes and principles of RNA folding.

Using an L7Ae-Tethered, Hydroxyl Radical-Mediated Footprinting Strategy to Identify and Validate Kink-Turns in RNAs.

Kink-turns are important RNA structural modules that facilitate long-range tertiary interactions and form binding sites for members of the L7Ae family of proteins. Present in a wide variety of functional RNAs, kink-turns play key organizational roles in many RNA-based cellular processes, including translation, modification, and tRNA biogenesis. It is important to determine the contribution of kink-turns to the overall architecture of resident RNAs, as these modules dictate ribonucleoprotein (RNP) assembly and function. This chapter describes a site-directed, hydroxyl radical-mediated footprinting strategy that utilizes L7Ae-tethered chemical nucleases to experimentally validate computationally identified kink-turns in any RNA and under a wide variety of conditions. The work plan described here uses the catalytic RNase P RNA as an example to provide a blueprint for using this footprinting method to map RNA-protein interactions in other RNP complexes.

Metal ions and sugar puckering balance single-molecule kinetic heterogeneity in RNA and DNA tertiary contacts

The fidelity of group II intron self-splicing and retrohoming relies on long-range tertiary interactions between the intron and its flanking exons. By single-molecule FRET, we explore the binding kinetics of the most important, structurally conserved contact, the exon and intron binding site 1 (EBS1/IBS1). A comparison of RNA-RNA and RNA-DNA hybrid contacts identifies transient metal ion binding as a major source of kinetic heterogeneity which typically appears in the form of degenerate FRET states. Molecular dynamics simulations suggest a structural link between heterogeneity and the sugar conformation at the exon-intron binding interface. While Mg2+ ions lock the exon in place and give rise to long dwell times in the exon bound FRET state, sugar puckering alleviates this structural rigidity and likely promotes exon release. The interplay of sugar puckering and metal ion coordination may be an important mechanism to balance binding affinities of RNA and DNA interactions in general.

Mining for recurrent long-range interactions in RNA structures reveals embedded hierarchies in network families.

The wealth of the combinatorics of nucleotide base pairs enables RNA molecules to assemble into sophisticated interaction networks, which are used to create complex 3D substructures. These interaction networks are essential to shape the 3D architecture of the molecule, and also to provide the key elements to carry molecular functions such as protein or ligand binding. They are made of organised sets of long-range tertiary interactions which connect distinct secondary structure elements in 3D structures. Here, we present a de novo data-driven approach to extract automatically from large data sets of full RNA 3D structures the recurrent interaction networks (RINs). Our methodology enables us for the first time to detect the interaction networks connecting distinct components of the RNA structure, highlighting their diversity and conservation through non-related functional RNAs. We use a graphical model to perform pairwise comparisons of all RNA structures available and to extract RINs and modules. Our analysis yields a complete catalog of RNA 3D structures available in the Protein Data Bank and reveals the intricate hierarchical organization of the RNA interaction networks and modules. We assembled our results in an online database (http://carnaval.lri.fr) which will be regularly updated. Within the site, a tool allows users with a novel RNA structure to detect automatically whether the novel structure contains previously observed RINs.

Is the Conformational Ensemble of Alzheimer's Aβ10-40 Peptide Force Field Dependent?

By applying REMD simulations we have performed comparative analysis of the conformational ensembles of amino-truncated Aβ10-40 peptide produced with five force fields, which combine four protein parameterizations (CHARMM36, CHARMM22*, CHARMM22/cmap, and OPLS-AA) and two water models (standard and modified TIP3P). Aβ10-40 conformations were analyzed by computing secondary structure, backbone fluctuations, tertiary interactions, and radius of gyration. We have also calculated Aβ10-40 3JHNHα-coupling and RDC constants and compared them with their experimental counterparts obtained for the full-length Aβ1-40 peptide. Our study led us to several conclusions. First, all force fields predict that Aβ adopts unfolded structure dominated by turn and random coil conformations. Second, specific TIP3P water model does not dramatically affect secondary or tertiary Aβ10-40 structure, albeit standard TIP3P model favors slightly more compact states. Third, although the secondary structures observed in CHARMM36 and CHARMM22/cmap simulations are qualitatively similar, their tertiary interactions show little consistency. Fourth, two force fields, OPLS-AA and CHARMM22* have unique features setting them apart from CHARMM36 or CHARMM22/cmap. OPLS-AA reveals moderate β-structure propensity coupled with extensive, but weak long-range tertiary interactions leading to Aβ collapsed conformations. CHARMM22* exhibits moderate helix propensity and generates multiple exceptionally stable long- and short-range interactions. Our investigation suggests that among all force fields CHARMM22* differs the most from CHARMM36. Fifth, the analysis of 3JHNHα-coupling and RDC constants based on CHARMM36 force field with standard TIP3P model led us to an unexpected finding that in silico Aβ10-40 and experimental Aβ1-40 constants are generally in better agreement than these quantities computed and measured for identical peptides, such as Aβ1-40 or Aβ1-42. This observation suggests that the differences in the conformational ensembles of Aβ10-40 and Aβ1-40 are small and the former can be used as proxy of the full-length peptide. Based on this argument, we concluded that CHARMM36 force field with standard TIP3P model produces the most accurate representation of Aβ10-40 conformational ensemble.

HIV-1 Pr55Gag binds genomic and spliced RNAs with different affinity and stoichiometry

ABSTRACTThe HIV-1 Pr55Gag precursor specifically selects genomic RNA (gRNA) from a large variety of cellular and spliced viral RNAs (svRNAs), however the molecular mechanisms of this selective recognition remains poorly understood. To gain better understanding of this process, we analyzed the interactions between Pr55Gag and a large panel of viral RNA (vRNA) fragments encompassing the main packaging signal (Psi) and its flanking regions by fluorescence spectroscopy. We showed that the gRNA harbors a high affinity binding site which is absent from svRNA species, suggesting that this site might be crucial for selecting the HIV-1 genome. Our stoichiometry analysis of protein/RNA complexes revealed that few copies of Pr55Gag specifically associate with the 5′ region of the gRNA. Besides, we found that gRNA dimerization significantly impacts Pr55Gag binding, and we confirmed that the internal loop of stem-loop 1 (SL1) in Psi is crucial for specific interaction with Pr55Gag. Our analysis of gRNA fragments of different length supports the existence of a long-range tertiary interaction involving sequences upstream and downstream of the Psi region. This long-range interaction might promote optimal exposure of SL1 for efficient Pr55Gag recognition. Altogether, our results shed light on the molecular mechanisms allowing the specific selection of gRNA by Pr55Gag among a variety of svRNAs, all harboring SL1 in their first common exon.

Single-molecule analysis reveals multi-state folding of a guanine riboswitch.

Guanine-responsive riboswitches undergo ligand-dependent structural rearrangements to control gene expression by transcription termination. While the molecular basis for ligand recognition is well established, the associated structural rearrangements and the kinetics involved in the formation of the aptamer domain are less well understood. Using high-resolution optical tweezers, we followed the folding trajectories of a single molecule of the xpt-pbuX guanine aptamer from Bacillus subtilis. We report a rapid six-state conformational rearrangement, in which three of the states are guanine dependent, during the transition from the linear to the native receptor conformation. The folding completes in <1 s. The force-dependent equilibrium kinetics and the mutational data indicated that the flexible J2-J3 junction undergoes a ligand-dependent conformational switching, which triggers the formation of the long-range tertiary interactions and the P1 helix. In the absence of the right ligand, the junction failed to initiate the series of conformational rearrangements required for the riboswitch activities.

Differential copper binding to alpha-synuclein and its disease-associated mutants affect the aggregation and amyloid formation

BackgroundCopper is an essential trace element required for the proper functioning of various enzymes present in the central nervous system. An imbalance in the copper homeostasis results in the pathology of various neurodegenerative disorders including Parkinson’s Disease. Hence, residue specific interaction of Cu2+ to α-Syn along with the familial mutants H50Q and G51D needs to be studied in detail. MethodsWe investigated the residue specific mapping of Cu2+ binding sites and binding strength using solution-state NMR and ITC respectively. The aggregation kinetics, secondary structural changes, and morphology of the formed fibrils in the presence and absence of Cu2+ were studied using fluorescence, CD, and AFM respectively. ResultsCopper binding to α-Syn takes place at three different sites with a higher affinity for the region 48-53. While one of the sites got abolished in the case of H50Q, the mutant G51D showed a binding pattern similar to WT. The aggregation kinetics of these proteins in the presence of Cu2+ showed an enhanced rate of fibril formation with a pronounced effect for G51D. ConclusionCu2+ binding results in the destabilization of long-range tertiary interactions in α-Syn leading to the exposure of highly amyloidogenic NAC region which results in the increased rate of fibril formation. Although the residues 48-53 have a stronger affinity for Cu2+ in case of WT and G51D, the binding is not responsible for enhancing the rate of fibril formation in case of H50Q. General SignificanceThese findings will help in the better understanding of Cu2+ catalyzed aggregation of synucleins.

Observation of long-range tertiary interactions during ligand binding by the TPP riboswitch aptamer.

The thiamine pyrophosphate (TPP) riboswitch is a cis-regulatory element in mRNA that modifies gene expression in response to TPP concentration. Its specificity is dependent upon conformational changes that take place within its aptamer domain. Here, the role of tertiary interactions in ligand binding was studied at the single-molecule level by combined force spectroscopy and Förster resonance energy transfer (smFRET), using an optical trap equipped for simultaneous smFRET. The 'Force-FRET' approach directly probes secondary and tertiary structural changes during folding, including events associated with binding. Concurrent transitions observed in smFRET signals and RNA extension revealed differences in helix-arm orientation between two previously-identified ligand-binding states that had been undetectable by spectroscopy alone. Our results show that the weaker binding state is able to bind to TPP, but is unable to form a tertiary docking interaction that completes the binding process. Long-range tertiary interactions stabilize global riboswitch structure and confer increased ligand specificity.

Challenge of mimicking the influences of the cellular environment on RNA structure by PEG-induced macromolecular crowding.

There are large differences between the cellular environment and the conditions widely used to study RNA in vitro. SHAPE RNA structure probing in Escherichia coli cells has shown that the cellular environment stabilizes both long-range and local tertiary interactions in the adenine riboswitch aptamer domain. Synthetic crowding agents are widely used to understand the forces that stabilize RNA structure and in efforts to recapitulate the cellular environment under simplified experimental conditions. Here, we studied the structure and ligand binding ability of the adenine riboswitch in the presence of the macromolecular crowding agent, polyethylene glycol (PEG). Ethylene glycol and low-molecular mass PEGs destabilized RNA structure and caused the riboswitch to sample secondary structures different from those observed in simple buffered solutions or in cells. In the presence of larger PEGs, longer-range loop-loop interactions were more similar to those in cells than in buffer alone, consistent with prior work showing that larger PEGs stabilize compact RNA states. Ligand affinity was weakened by low-molecular mass PEGs but increased with high-molecular mass PEGs, indicating that PEG cosolvents exert complex chemical and steric effects on RNA structure. Regardless of polymer size, however, nucleotide-resolution structural characteristics observed in cells were not recapitulated in PEG solutions. Our results reveal that the cellular environment is difficult to recapitulate in vitro; mimicking the cellular state will likely require a combination of crowding agents and other chemical species.

RNA folding retrospective: lessons from ribozymes big and small.

RNA folding retrospective: lessons from ribozymes big and small.

The "complex" RNA post-transcriptional element of a "simple" retrovirus.

Replication of retroviruses and transposition of endogenous retroelements exploits a unique mechanism of post-transcriptional regulation as a means of exporting their incompletely-spliced mRNAs (which serve as both the genomic RNA and the template for protein synthesis). Following discovery of the Rev response element (RRE) that mediates nucleocytoplasmic export of the full-length and singly-spliced human immunodeficiency virus type 1 (HIV-1) genome, equivalent cis-acting regulatory elements have been characterized for both complex and simple retroviruses and retroelements, together with the obligate viral and host proteins with which they interact. The exception to this is the gammaretrovirus family of simple retroviruses, exemplified by reticuloendotheliosis virus (REV), murine leukemia virus (MLV) and xenotropic MLV-related retrovirus (XMRV). In this commentary, we discuss our recent data that reported structural and functional data on the MLV/XMRV post-transcriptional regulatory element (designated the PTE). The PTE was characterized by a highly-structured region of multiple stem-loops (SL1 – SL7) overlapping the pro and 5′ portion of the pol open reading frames, comprising a bipartite export signal whose structures are separated by ∼1400 nt. In addition, structural probing suggested that SL3 nucleotides were involved in pseudoknot formation. These data, when compared with RNA transport elements of complex retroviruses (HIV) and simple murine retrotransposons (musD), collectively present an emerging picture that long-range tertiary interactions are critical mediators of their biological function.

Nuclear magnetic resonance structure of the III-IV-V three-way junction from the Varkud satellite ribozyme and identification of magnesium-binding sites using paramagnetic relaxation enhancement.

The VS ribozyme is a catalytic RNA found within some natural isolates of Neurospora that is being used as a model system to improve our understanding of RNA structure, catalysis, and engineering. The catalytic domain contains five helical domains (SLII-SLVI) that are organized by two three-way junctions. The III-IV-V junction is required for high-affinity binding of the substrate domain (SLI) through formation of a kissing loop interaction with SLV. Here, we determine the high-resolution nuclear magnetic resonance (NMR) structure of a 47-nucleotide RNA containing the III-IV-V junction (J345). The J345 RNA adopts a Y-shaped fold typical of the family C three-way junctions, with coaxial stacking between stems III and IV and an acute angle between stems III and V. The NMR structure reveals that the core of the III-IV-V junction contains four stacked base triples, a U-turn motif, a cross-strand stacking interaction, an A-minor interaction, and a ribose zipper. In addition, the NMR structure shows that the cCUUGg tetraloop used to stabilize stem IV adopts a novel RNA tetraloop fold, different from the known gCUUGc tetraloop structure. Using Mn(2+)-induced paramagnetic relaxation enhancement, we identify six Mg(2+)-binding sites within J345, including one associated with the cCUUGg tetraloop and two with the junction core. The NMR structure of J345 likely represents the conformation of the III-IV-V junction in the context of the active VS ribozyme and suggests that this junction functions as a dynamic hinge that contributes to substrate recognition and catalysis. Moreover, this study highlights a new role for family C three-way junctions in long-range tertiary interactions.

Roles of long-range tertiary interactions in limiting dynamics of the Tetrahymena group I ribozyme.

We determined the effects of mutating the long-range tertiary contacts of the Tetrahymena group I ribozyme on the dynamics of its substrate helix (referred to as P1) and on catalytic activity. Dynamics were assayed by fluorescence anisotropy of the fluorescent base analogue, 6-methyl isoxanthopterin, incorporated into the P1 helix, and fluorescence anisotropy and catalytic activity were measured for wild type and mutant ribozymes over a range of conditions. Remarkably, catalytic activity correlated with P1 anisotropy over 5 orders of magnitude of activity, with a correlation coefficient of 0.94. The functional and dynamic effects from simultaneous mutation of the two long-range contacts that weaken P1 docking are cumulative and, based on this RNA’s topology, suggest distinct underlying origins for the mutant effects. Tests of mechanistic predictions via single molecule FRET measurements of rate constants for P1 docking and undocking suggest that ablation of the P14 tertiary interaction frees P2 and thereby enhances the conformational space explored by the undocked attached P1 helix. In contrast, mutation of the metal core tertiary interaction disrupts the conserved core into which the P1 helix docks. Thus, despite following a single correlation, the two long-range tertiary contacts facilitate P1 helix docking by distinct mechanisms. These results also demonstrate that a fluorescence anisotropy probe incorporated into a specific helix within a larger RNA can report on changes in local helical motions as well as differences in more global dynamics. This ability will help uncover the physical properties and behaviors that underlie the function of RNAs and RNA/protein complexes.

Loop-Sequence Features and Stability Determinants in Antibody Variable Domains by High-Throughput Experiments

Protein loops are frequently considered as critical determinants in protein structure and function. Recent advances in high-throughput methods for DNA sequencing and thermal stability measurement have enabled effective exploration of sequence-structure-function relationships in local protein regions. Using these data-intensive technologies, we investigated the sequence-structure-function relationships of six complementarity-determining regions (CDRs) and ten non-CDR loops in the variable domains of a model vascular endothelial growth factor (VEGF)-binding single-chain antibody variable fragment (scFv) whose sequence had been optimized via a consensus-sequence approach. The results show that only a handful of residues involving long-range tertiary interactions distant from the antigen-binding site are strongly coupled with antigen binding. This implies that the loops are passive regions in protein folding; the essential sequences of these regions are dictated by conserved tertiary interactions and the consensus local loop-sequence features contribute little to protein stability and function.