While functional roles of several long non-coding RNAs (lncRNAs) have been determined, the molecular mechanisms are not well understood. Here, we report the first experimentally derived secondary structure of a human lncRNA, the steroid receptor RNA activator (SRA), 0.87 kB in size. The SRA RNA is a non-coding RNA that coactivates several human sex hormone receptors and is strongly associated with breast cancer. Coding isoforms of SRA are also expressed to produce proteins, making the SRA gene a unique bifunctional system. Our experimental findings (SHAPE, in-line, DMS and RNase V1 probing) reveal that this lncRNA has a complex structural organization, consisting of four domains, with a variety of secondary structure elements. We examine the coevolution of the SRA gene at the RNA structure and protein structure levels using comparative sequence analysis across vertebrates. Rapid evolutionary stabilization of RNA structure, combined with frame-disrupting mutations in conserved regions, suggests that evolutionary pressure preserves the RNA structural core rather than its translational product. We perform similar experiments on alternatively spliced SRA isoforms to assess their structural features.

Decoding is the key step during protein synthesis that enables information transfer from RNA to protein, a process critical for the survival of all organisms. We have used large-scale (2.64 × 10 ⁶ atoms) all-atom simulations of the entire ribosome to understand a critical step of decoding. Although the decoding problem has been studied for more than four decades, the rate-limiting step of cognate tRNA selection has only recently been identified. This step, known as accommodation, involves the movement inside the ribosome of the aminoacyl-tRNA from the partially bound “A/T” state to the fully bound “A/A” state. Here, we show that a corridor of 20 universally conserved ribosomal RNA bases interacts with the tRNA during the accommodation movement. Surprisingly, the tRNA is impeded by the A-loop (23S helix 92), instead of enjoying a smooth transition to the A/A state. In particular, universally conserved 23S ribosomal RNA bases U2492, C2556, and C2573 act as a 3D gate, causing the acceptor stem to pause before allowing entrance into the peptidyl transferase center. Our simulations demonstrate that the flexibility of the acceptor stem of the tRNA, in addition to flexibility of the anticodon arm, is essential for tRNA selection. This study serves as a template for simulating conformational changes in large (>10 ⁶ atoms) biological and artificial molecular machines.

The ribosome is a massive ribonucleoprotein complex (∼2.4 MDa) that utilizes large-scale structural fluctuations to produce unidirectional protein synthesis. Accommodation is a key conformational change during transfer RNA (tRNA) selection that allows movement of tRNA into the ribosome. Here, we address the structure–function relationship that governs accommodation using all-atom molecular simulations and single-molecule fluorescence resonance energy transfer (smFRET). Simulations that employ an all-atom, structure-based (Gō-like) model illuminate the interplay between configurational entropy and effective enthalpy during the accommodation process. This delicate balance leads to spontaneous reversible accommodation attempts, which are corroborated by smFRET measurements. The dynamics about the endpoints of accommodation (the A/T and A/A conformations) obtained from structure-based simulations are validated by multiple 100–200 ns explicit-solvent simulations (3.2 million atoms for a cumulative 1.4 μs), and previous crystallographic analysis. We find that the configurational entropy of the 3′-CCA end of aminoacyl-tRNA resists accommodation, leading to a multistep accommodation process that encompasses a distribution of parallel pathways. The calculated mechanism is robust across simulation methods and protocols, suggesting that the structure of the accommodation corridor imposes stringent limitations on the accessible pathways. The identified mechanism and observed parallel pathways establish an atomistic framework for interpreting a large body of biochemical data and demonstrate that conformational changes during translation occur through a stochastic trial-and-error process, rather than in concerted lock-step motions.

There is considerable debate about the functionality of long non-coding RNAs (lncRNAs). Lack of sequence conservation has been used to argue against functional relevance. We investigated antisense lncRNAs, called COOLAIR, at the A. thaliana FLC locus and experimentally determined their secondary structure. The major COOLAIR variants are highly structured, organized by exon. The distally polyadenylated transcript has a complex multi-domain structure, altered by a single non-coding SNP defining a functionally distinct A. thaliana FLC haplotype. The A. thaliana COOLAIR secondary structure was used to predict COOLAIR exons in evolutionarily divergent Brassicaceae species. These predictions were validated through chemical probing and cloning. Despite the relatively low nucleotide sequence identity, the structures, including multi-helix junctions, show remarkable evolutionary conservation. In a number of places, the structure is conserved through covariation of a non-contiguous DNA sequence. This structural conservation supports a functional role for COOLAIR transcripts rather than, or in addition to, antisense transcription.