Identification of receptors for UNCG and GNRA Z-turns and their occurrence in rRNA

Biophysical characterization and design of a minimal version of the Hoechst RNA aptamer

RNA aptamers are oligonucleotides, selected through Systematic Evolution of Ligands by EXponential Enrichment (SELEX), that can bind to specific target molecules with high affinity. One such molecule is the RNA aptamer that binds to a blue-fluorescent Hoechst dye that was modified with bulky t-Bu groups to prevent non-specific binding to DNA. This aptamer has potential for biosensor applications; however, limited information is available regarding its conformation, molecular interactions with the ligand, and binding mechanism. The study presented here aims to biophysically characterize the Hoechst RNA aptamer when complexed with the t-Bu Hoechst dye and to further optimize the RNA sequence by designing and synthesizing new sequence variants. Each variant aptamer-t-Bu Hoechst complex was evaluated through a combination of fluorescence emission, native polyacrylamide gel electrophoresis, fluorescence titration, and isothermal titration calorimetry experiments. The results were used to design a minimal version of the aptamer consisting of only 21 nucleotides. The performed study also describes a more efficient method for synthesizing the t-Bu Hoechst dye derivative. Understanding the biophysical properties of the t-Bu Hoechst dye-RNA complex lays the foundation for nuclear magnetic resonance spectroscopy studies and its potential development as a building block for an aptamer-based biosensor that can be used in medical, environmental or laboratory settings.

Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions

Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r^-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.

Occurrence and stability of anion-π interactions between phosphate and nucleobases in functional RNA molecules

We present a systematic structural and energetic characterization of phosphate(OP)-nucleobase anion…π stacking interactions in RNAs. We observed OP-nucleobase stacking contacts in a variety of structural motifs other than regular helices and spanning broadly diverse sequence distances. Apart from the stacking between a phosphate and a guanine or a uracil two-residue upstream in specific U-turns, such interactions in RNA have been scarcely characterized to date. Our QM calculations showed an energy minimum at a distance between the OP atom and the nucleobase plane centroid slightly below 3 Å for all the nucleobases. By sliding the OP atom over the nucleobase plane we localized the optimal mutual positioning of the stacked moieties, corresponding to an energy minimum below -6 kcal•mol-1, for all the nucleobases, consistently with the projections of the OP atoms over the different π-rings we observed in experimental occurrences. We also found that the strength of the interaction clearly correlates with its electrostatic component, pointing to it as the most relevant contribution. Finally, as OP-uracil and OP-guanine interactions represent together 86% of the instances we detected, we also proved their stability under dynamic conditions in model systems simulated by state-of-the art DFT-MD calculations.

Lone Pair…π Contacts and Structure Signatures of r(UNCG) Tetraloops, Z-Turns, and Z-Steps: A WebFR3D Survey

Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.

Engineering Crystal Packing in RNA Structures I: Past and Future Strategies for Engineering RNA Packing in Crystals

X-ray crystallography remains a powerful method to gain atomistic insights into the catalytic and regulatory functions of RNA molecules. However, the technique requires the preparation of diffraction-quality crystals. This is often a resource- and time-consuming venture because RNA crystallization is hindered by the conformational heterogeneity of RNA, as well as the limited opportunities for stereospecific intermolecular interactions between RNA molecules. The limited success at crystallization explains in part the smaller number of RNA-only structures in the Protein Data Bank. Several approaches have been developed to aid the formation of well-ordered RNA crystals. The majority of these are construct-engineering techniques that aim to introduce crystal contacts to favor the formation of well-diffracting crystals. A typical example is the insertion of tetraloop-tetraloop receptor pairs into non-essential RNA segments to promote intermolecular association. Other methods of promoting crystallization involve chaperones and crystallization-friendly molecules that increase RNA stability and improve crystal packing. In this review, we discuss the various techniques that have been successfully used to facilitate crystal packing of RNA molecules, recent advances in construct engineering, and directions for future research in this vital aspect of RNA crystallography.

Special Issue: A, B and Z: The Structure, Function and Genetics of Z-DNA and Z-RNA

It is now difficult to believe that a biological function for the left-handed Z-DNA and Z-RNA conformations was once controversial. The papers in this Special Issue, "Z-DNA and Z-RNA: from Physical Structure to Biological Function", are based on presentations at the ABZ2021 meeting that was held virtually on 19 May 2021 and provide evidence for several biological functions of these structures. The first of its kind, this international conference gathered over 200 scientists from many disciplines to specifically address progress in research involving Z-DNA and Z-RNA. These high-energy left-handed conformers of B-DNA and A-RNA are associated with biological functions and disease outcomes, as evidenced from both mouse and human genetic studies. These alternative structures, referred to as "flipons", form under physiological conditions, regulate type I interferon responses and induce necroptosis during viral infection. They can also stimulate genetic instability, resulting in adaptive evolution and diseases such as cancer. The meeting featured cutting-edge science that was, for the most part, unpublished. We plan for the ABZ meeting to reconvene in 2022.

Recognition of non-CpG repeats in Alu and ribosomal RNAs by the Z-RNA binding domain of ADAR1 induces A-Z junctions

Adenosine-to-inosine (A-to-I) editing of eukaryotic cellular RNAs is essential for protection against auto-immune disorders. Editing is carried out by ADAR1, whose innate immune response-specific cytoplasmic isoform possesses a Z-DNA binding domain (Zα) of unknown function. Zα also binds to CpG repeats in RNA, which are a hallmark of Z-RNA formation. Unexpectedly, Zα has been predicted — and in some cases even shown — to bind to specific regions within mRNA and rRNA devoid of such repeats. Here, we use NMR, circular dichroism, and other biophysical approaches to demonstrate and characterize the binding of Zα to mRNA and rRNA fragments. Our results reveal a broad range of RNA sequences that bind to Zα and adopt Z-RNA conformations. Binding is accompanied by destabilization of neighboring A-form regions which is similar in character to what has been observed for B-Z-DNA junctions. The binding of Zα to non-CpG sequences is specific, cooperative and occurs with an affinity in the low micromolar range. This work allows us to propose a model for how Zα could influence the RNA binding specificity of ADAR1.

ADAR1 is an interferon-induced enzyme that catalyzes editing of adenine to inosine across the transcriptome as part of the immune response. Here the authors establish how ADAR1 recognizes non-CpG RNA sequences to facilitate the formation of A-Z junctions.

RNA stabilization by a poly(A) tail 3'-end binding pocket and other modes of poly(A)-RNA interaction

Polyadenylate [poly(A)] tail addition to the 3' end of a wide range of RNAs is a highly conserved modification that plays a central role in cellular RNA function. Elements for nuclear expression (ENEs) are cis-acting RNA elements that stabilize poly(A) tails by sequestering them in RNA triplex structures. A crystal structure of a double ENE from a rice hAT transposon messenger RNA complexed with poly(A)₂₈ at a resolution of 2.89 angstroms reveals multiple modes of interaction with poly(A), including major-groove triple helices, extended minor-groove interactions with RNA double helices, a quintuple-base motif that transitions poly(A) from minor-groove associations to major-groove triple helices, and a poly(A) 3'-end binding pocket. Our findings both expand the repertoire of motifs involved in long-range RNA interactions and provide insights into how polyadenylation can protect an RNA's extreme 3' end.

Structural distinctions between NAD+ riboswitch domains 1 and 2 determine differential folding and ligand binding

Riboswitches are important gene regulatory elements frequently encountered in bacterial mRNAs. The recently discovered nadA riboswitch contains two similar, tandemly arrayed aptamer domains, with the first domain possessing high affinity for nicotinamide adenine dinucleotide (NAD⁺). The second domain which comprises the ribosomal binding site in a putative regulatory helix, however, has withdrawn from detection of ligand-induced structural modulation thus far, and therefore, the identity of the cognate ligand and the regulation mechanism have remained unclear. Here, we report crystal structures of both riboswitch domains, each bound to NAD⁺. Furthermore, we demonstrate that ligand binding to domain 2 requires significantly higher concentrations of NAD⁺ (or ADP retaining analogs) compared to domain 1. Using a fluorescence spectroscopic approach, we further shed light on the structural features which are responsible for the different ligand affinities, and describe the Mg²⁺-dependent, distinct folding and pre-organization of their binding pockets. Finally, we speculate about possible scenarios for nadA RNA gene regulation as a putative two-concentration sensor module for a time-controlled signal that is primed and stalled by the gene regulation machinery at low ligand concentrations (domain 1), and finally triggers repression of translation as soon as high ligand concentrations are reached in the cell (domain 2).

Conformational adaptation of UNCG loops upon crowding

If the A-form helix is the major structural motif found in RNA, the loops that cap them constitute the second most important family of motifs. Among those, two are overrepresented, GNRA and UNCG tetraloops. Recent surveys of RNA structures deposited in the PDB show that GNRA and UNCG tetraloops can adopt tertiary folds that are very different from their canonical conformations, characterized by the presence of a U-turn of a Z-turn, respectively. Crystallographic data from both a lariat-capping (LC) ribozyme and a group II intron ribozyme reveal that a given UUCG tetraloop can adopt a distinct fold depending on its structural environment. Specifically, when the crystal packing applies relaxed constraints on the loop, the canonical Z-turn conformation is observed. In contrast, a highly packed environment induces "squashing" of the tetraloop by distorting its sugar-phosphate backbone in a specific way that expels the first and fourth nucleobases out of the loop, and falls in van der Waals distance of the last base pair of the helix, taking the place of the pair formed between the first and fourth residues in Z-turn loops. The biological relevance of our observations is supported by the presence of similarly deformed loops in the highly packed environment of the ribosome and in a complex between a dsRNA and a RNase III. The finding that Z-turn loops change conformation under higher molecular packing suggests that, in addition to their demonstrated role in stabilizing RNA folding, they may contribute to the three-dimensional structure of RNA by mediating tertiary interactions with distal residues.

Deducing putative ancestral forms of GNRA/receptor interactions from the ribosome

Stable RNAs rely on a vast repertoire of long-range interactions to assist in the folding of complex cellular machineries such as the ribosome. The universally conserved L39/H89 interaction is a long-range GNRA-like/receptor interaction localized in proximity to the peptidyl transferase center of the large subunit of the ribosome. Because of its central location, L39/H89 likely originated at an early evolutionary stage of the ribosome and played a significant role in its early function. However, L39/H89 self-assembly is impaired outside the ribosomal context. Herein, we demonstrate that structural modularity principles can be used to re-engineer L39/H89 to self-assemble in vitro. The new versions of L39/H89 improve affinity and loop selectivity by several orders of magnitude and retain the structural and functional features of their natural counterparts. These versions of L39/H89 are proposed to be ancestral forms of L39/H89 that were capable of assembling and folding independently from proteins and post-transcriptional modifications. This work demonstrates that novel RNA modules can be rationally designed by taking advantage of the modular syntax of RNA. It offers the prospect of creating new biochemical models of the ancestral ribosome and increases the tool kit for RNA nanotechnology and synthetic biology.

Z-DNA and Z-RNA in human disease

Left-handed Z-DNA/Z-RNA is bound with high affinity by the Zα domain protein family that includes ADAR (a double-stranded RNA editing enzyme), ZBP1 and viral orthologs regulating innate immunity. Loss-of-function mutations in ADAR p150 allow persistent activation of the interferon system by Alu dsRNAs and are causal for Aicardi-Goutières Syndrome. Heterodimers of ADAR and DICER1 regulate the switch from RNA- to protein-centric immunity. Loss of DICER1 function produces age-related macular degeneration, a different type of Alu-mediated disease. The overlap of Z-forming sites with those for the signal recognition particle likely limits invasion of primate genomes by Alu retrotransposons.