Structure of an internal loop motif with three consecutive U•U mismatches from stem-loop 1 in the 3'-UTR of the SARS-CoV-2 genomic RNA.

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups. However, their conformational and dynamical properties cannot be directly inferred from chemical probing and are difficult to predict theoretically. A mismatch motif (SL1-motif) consisting of three consecutive U•U base pairs is located in stem-loop 1 of the 3'-UTR. We combined NMR-spectroscopy and MD-simulations to investigate its structure and dynamics. All three U•U base pairs feature two direct hydrogen bonds and are as stable as Watson-Crick A:U base pairs. Plasmodium falciparum 25S rRNA contains a triple U•U mismatch motif (Pf-motif) differing from SL1-motif only with respect to the orientation of the two closing base pairs. Interestingly, while the geometry of the outer two U•U mismatches was identical in both motifs the preferred orientation of the central U•U mismatch was different. MD simulations and potassium ion titrations revealed that the potassium ion-binding mode to the major groove is connected to the different preferred geometries of the central base pair in the two motifs.

Structural and dynamic effects of pseudouridine modifications on noncanonical interactions in RNA.

Pseudouridine is the most frequently naturally occurring RNA modification, found in all classes of biologically functional RNAs. Compared to uridine, pseudouridine contains an additional hydrogen bond donor group and is therefore widely regarded as a structure stabilizing modification. However, the effects of pseudouridine modifications on the structure and dynamics of RNAs have so far only been investigated in a limited number of different structural contexts. Here, we introduced pseudouridine modifications into the U-turn motif and the adjacent U:U closing base pair of the neomycin-sensing riboswitch (NSR)-an extensively characterized model system for RNA structure, ligand binding, and dynamics. We show that the effects of replacing specific uridines with pseudouridines on RNA dynamics crucially depend on the exact location of the replacement site and can range from destabilizing to locally or even globally stabilizing. By using a combination of NMR spectroscopy, MD simulations and QM calculations, we rationalize the observed effects on a structural and dynamical level. Our results will help to better understand and predict the consequences of pseudouridine modifications on the structure and function of biologically important RNAs.

Development of a novel tobramycin dependent riboswitch.

We herein report the selection and characterization of a new riboswitch dependent on the aminoglycoside tobramycin. Its dynamic range rivals even the tetracycline dependent riboswitch to be the current best performing, synthetic riboswitch that controls translation initiation. The riboswitch was selected with RNA Capture-SELEX, a method that not only selects for binding but also for structural changes in aptamers on binding. This study demonstrates how this method can fundamentally reduce the labour required for the de novo identification of synthetic riboswitches. The initially selected riboswitch candidate harbours two distinct tobramycin binding sites with KDs of 1.1 nM and 2.4 µM, respectively, and can distinguish between tobramycin and the closely related compounds kanamycin A and B. Using detailed genetic and biochemical analyses and 1H NMR spectroscopy, the proposed secondary structure of the riboswitch was verified and the tobramycin binding sites were characterized. The two binding sites were found to be essentially non-overlapping, allowing for a separate investigation of their contribution to the activity of the riboswitch. We thereby found that only the high-affinity binding site was responsible for regulatory activity, which allowed us to engineer a riboswitch from only this site with a minimal sequence size of 33 nt and outstanding performance.

High-resolution structure of stem-loop 4 from the 5'-UTR of SARS-CoV-2 solved by solution state NMR.

We present the high-resolution structure of stem-loop 4 of the 5'-untranslated region (5_SL4) of the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) genome solved by solution state nuclear magnetic resonance spectroscopy. 5_SL4 adopts an extended rod-like structure with a single flexible looped-out nucleotide and two mixed tandem mismatches, each composed of a G•U wobble base pair and a pyrimidine•pyrimidine mismatch, which are incorporated into the stem-loop structure. Both the tandem mismatches and the looped-out residue destabilize the stem-loop structure locally. Their distribution along the 5_SL4 stem-loop suggests a role of these non-canonical elements in retaining functionally important structural plasticity in particular with regard to the accessibility of the start codon of an upstream open reading frame located in the RNA's apical loop. The apical loop-although mostly flexible-harbors residual structural features suggesting an additional role in molecular recognition processes. 5_SL4 is highly conserved among the different variants of SARS-CoV-2 and can be targeted by small molecule ligands, which it binds with intermediate affinity in the vicinity of the non-canonical elements within the stem-loop structure.

Multi-Site Conformational Exchange in the Synthetic Neomycin-Sensing Riboswitch Studied by 19 F NMR.

The synthetic neomycin-sensing riboswitch interacts with its cognate ligand neomycin as well as with the related antibiotics ribostamycin and paromomycin. Binding of these aminoglycosides induces a very similar ground state structure in the RNA, however, only neomycin can efficiently repress translation initiation. The molecular origin of these differences has been traced back to differences in the dynamics of the ligand:riboswitch complexes. Here, we combine five complementary fluorine based NMR methods to accurately quantify seconds to microseconds dynamics in the three riboswitch complexes. Our data reveal complex exchange processes with up to four structurally different states. We interpret our findings in a model that shows an interplay between different chemical groups in the antibiotics and specific bases in the riboswitch. More generally, our data underscore the potential of 19 F NMR methods to characterize complex exchange processes with multiple excited states.

Solution structure and RNA-binding of a minimal ProQ-homolog from Legionella pneumophila (Lpp1663).

Small regulatory RNAs (sRNAs) play an important role for posttranscriptional gene regulation in bacteria. sRNAs recognize their target messenger RNAs (mRNAs) by base-pairing, which is often facilitated by interactions with the bacterial RNA-binding proteins Hfq or ProQ. The FinO/ProQ RNA-binding protein domain was first discovered in the bacterial repressor of conjugation, FinO. Since then, the functional role of FinO/ProQ-like proteins in posttranscriptional gene regulation was extensively studied in particular in the enterobacteria E. coli and Salmonella enterica and a wide range of sRNA-targets was identified for these proteins. In addition, enterobacterial ProQ homologs also recognize and protect the 3'-ends of a number of mRNAs from exonucleolytic degradation. However, the RNA-binding properties of FinO/ProQ proteins with regard to the recognition of different RNA targets are not yet fully understood. Here, we present the solution NMR structure of the so far functionally uncharacterized ProQ homolog Lpp1663 from Legionella pneumophila as a newly confirmed member and a minimal model system of the FinO/ProQ protein family. In addition, we characterize the RNA-binding preferences of Lpp1663 with high resolution NMR spectroscopy and isothermal titration calorimetry (ITC). Our results suggest a binding preference for single-stranded uridine-rich RNAs in the vicinity of stable stem-loop structures. According to chemical shift perturbation experiments, the single-stranded U-rich RNAs interact mainly with a conserved RNA-binding surface on the concave site of Lpp1663.

Structure of an RNA aptamer in complex with the fluorophore tetramethylrhodamine.

RNA aptamers-artificially created RNAs with high affinity and selectivity for their target ligand generated from random sequence pools-are versatile tools in the fields of biotechnology and medicine. On a more fundamental level, they also further our general understanding of RNA-ligand interactions e. g. in regard to the relationship between structural complexity and ligand affinity and specificity, RNA structure and RNA folding. Detailed structural knowledge on a wide range of aptamer-ligand complexes is required to further our understanding of RNA-ligand interactions. Here, we present the atomic resolution structure of an RNA-aptamer binding to the fluorescent xanthene dye tetramethylrhodamine. The high resolution structure, solved by NMR-spectroscopy in solution, reveals binding features both common and different from the binding mode of other aptamers with affinity for ligands carrying planar aromatic ring systems such as the malachite green aptamer which binds to the tetramethylrhodamine related dye malachite green or the flavin mononucleotide aptamer.

Identification of the 3-amino-3-carboxypropyl (acp) transferase enzyme responsible for acp3U formation at position 47 in Escherichia coli tRNAs.

tRNAs from all domains of life contain modified nucleotides. However, even for the experimentally most thoroughly characterized model organism Escherichia coli not all tRNA modification enzymes are known. In particular, no enzyme has been found yet for introducing the acp3U modification at position 47 in the variable loop of eight E. coli tRNAs. Here we identify the so far functionally uncharacterized YfiP protein as the SAM-dependent 3-amino-3-carboxypropyl transferase catalyzing this modification and thereby extend the list of known tRNA modification enzymes in E. coli. Similar to the Tsr3 enzymes that introduce acp modifications at U or m1Ψ nucleotides in rRNAs this protein contains a DTW domain suggesting that acp transfer reactions to RNA nucleotides are a general function of DTW domain containing proteins. The introduction of the acp3U-47 modification in E. coli tRNAs is promoted by the presence of the m7G-46 modification as well as by growth in rich medium. However, a deletion of the enzymes responsible for the modifications at position 46 and 47 in the variable loop of E. coli tRNAs did not lead to a clearly discernible phenotype suggesting that these two modifications play only a minor role in ensuring the proper function of tRNAs in E. coli.

The structure of the SAM/SAH-binding riboswitch.

S-adenosylmethionine (SAM) is a central metabolite since it is used as a methyl group donor in many different biochemical reactions. Many bacteria control intracellular SAM concentrations using riboswitch-based mechanisms. A number of structurally different riboswitch families specifically bind to SAM and mainly regulate the transcription or the translation of SAM-biosynthetic enzymes. In addition, a highly specific riboswitch class recognizes S-adenosylhomocysteine (SAH)-the product of SAM-dependent methyl group transfer reactions-and regulates enzymes responsible for SAH hydrolysis. High-resolution structures are available for many of these riboswitch classes and illustrate how they discriminate between the two structurally similar ligands SAM and SAH. The so-called SAM/SAH riboswitch class binds both ligands with similar affinities and is structurally not yet characterized. Here, we present a high-resolution nuclear magnetic resonance structure of a member of the SAM/SAH-riboswitch class in complex with SAH. Ligand binding induces pseudoknot formation and sequestration of the ribosome binding site. Thus, the SAM/SAH-riboswitches are translational 'OFF'-switches. Our results establish a structural basis for the unusual bispecificity of this riboswitch class. In conjunction with genomic data our structure suggests that the SAM/SAH-riboswitches might be an evolutionary late invention and not a remnant of a primordial RNA-world as suggested for other riboswitches.

Cooperation between a T†Domain and a Minimal C-Terminal Docking Domain to Enable Specific Assembly in a Multiprotein NRPS.

Non-ribosomal peptide synthetases (NRPS) produce natural products from amino acid building blocks. They often consist of multiple polypeptide chains which assemble in a specific linear order via specialized N- and C-terminal docking domains (N/C DDs). Typically, docking domains function independently from other domains in NRPS assembly. Thus, docking domain replacements enable the assembly of "designer" NRPS from proteins that normally do not interact. The multiprotein "peptide-antimicrobial-Xenorhabdus" (PAX) peptide-producing PaxS NRPS is assembled from the three proteins PaxA, PaxB and PaxC. Herein, we show that the small C DD of PaxA cooperates with its preceding thiolation (T1 ) domain to bind the N DD of PaxB with very high affinity, establishing a structural and thermodynamical basis for this unprecedented docking interaction, and we test its functional importance in vivo in a truncated PaxS assembly line. Similar docking interactions are apparently present in other NRPS systems.

Exploring the Druggability of Conserved RNA Regulatory Elements in the SARS-CoV-2 Genome.

SARS-CoV-2 contains a positive single-stranded RNA genome of approximately 30 000 nucleotides. Within this genome, 15 RNA elements were identified as conserved between SARS-CoV and SARS-CoV-2. By nuclear magnetic resonance (NMR) spectroscopy, we previously determined that these elements fold independently, in line with data from in vivo and ex-vivo structural probing experiments. These elements contain non-base-paired regions that potentially harbor ligand-binding pockets. Here, we performed an NMR-based screening of a poised fragment library of 768 compounds for binding to these RNAs, employing three different 1 H-based 1D NMR binding assays. The screening identified common as well as RNA-element specific hits. The results allow selection of the most promising of the 15 RNA elements as putative drug targets. Based on the identified hits, we derive key functional units and groups in ligands for effective targeting of the RNA of SARS-CoV-2.

Adenine protonation enables cyclic-di-GMP binding to cyclic-GAMP sensing riboswitches.

In certain structural or functional contexts, RNA structures can contain protonated nucleotides. However, a direct role for stably protonated nucleotides in ligand binding and ligand recognition has not yet been demonstrated unambiguously. Previous X-ray structures of c-GAMP binding riboswitch aptamer domains in complex with their near-cognate ligand c-di-GMP suggest that an adenine of the riboswitch either forms two hydrogen bonds to a G nucleotide of the ligand in the unusual enol tautomeric form or that the adenine in its N1 protonated form binds the G nucleotide of the ligand in its canonical keto tautomeric state. By using NMR spectroscopy we demonstrate that the c-GAMP riboswitches bind c-di-GMP using a stably protonated adenine in the ligand binding pocket. Thereby, we provide novel insights into the putative biological functions of protonated nucleotides in RNA, which in this case influence the ligand selectivity in a riboswitch.

An intricate balance of hydrogen bonding, ion atmosphere and dynamics facilitates a seamless uracil to cytosine substitution in the U-turn of the neomycin-sensing riboswitch.

The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.

Small but large enough: structural properties of armless mitochondrial tRNAs from the nematode Romanomermis culicivorax.

As adapter molecules to convert the nucleic acid information into the amino acid sequence, tRNAs play a central role in protein synthesis. To fulfill this function in a reliable way, tRNAs exhibit highly conserved structural features common in all organisms and in all cellular compartments active in translation. However, in mitochondria of metazoans, certain dramatic deviations from the consensus tRNA structure are described, where some tRNAs lack the D- or T-arm without losing their function. In Enoplea, this miniaturization comes to an extreme, and functional mitochondrial tRNAs can lack both arms, leading to a considerable size reduction. Here, we investigate the secondary and tertiary structure of two such armless tRNAs from Romanomermis culicivorax. Despite their high AU content, the transcripts fold into a single and surprisingly stable hairpin structure, deviating from standard tRNAs. The three-dimensional form is boomerang-like and diverges from the standard L-shape. These results indicate that such unconventional miniaturized tRNAs can still fold into a tRNA-like shape, although their length and secondary structure are very unusual. They highlight the remarkable flexibility of the protein synthesis apparatus and suggest that the translational machinery of Enoplea mitochondria may show compensatory adaptations to accommodate these armless tRNAs for efficient translation.

Novel 13 C-detected NMR Experiments for the Precise Detection of RNA Structure.

Up to now, NMR spectroscopic investigations of RNA have utilized imino proton resonances as reporters for base pairing and RNA structure. The nucleobase amino groups are often neglected, since most of their resonances are broadened beyond detection due to rotational motion around the C-NH2 bond. Here, we present 13 C-detected NMR experiments for the characterization of all RNA amino groups irrespective of their motional behavior. We have developed a C(N)H-HDQC experiment that enables the observation of a complete set of sharp amino resonances through the detection of proton-NH2 double quantum coherences. Further, we present an "amino"-NOESY experiment to detect NOEs to amino protons, which are undetectable by any other conventional NOESY experiment. Together, these experiments allow the exploration of additional chemical shift information and inter-residual proton distances important for high-resolution RNA secondary and tertiary structure determination.

An intermolecular G-quadruplex as the basis for GTP recognition in the class V-GTP aptamer.

Many naturally occurring or artificially created RNAs are capable of binding to guanine or guanine derivatives with high affinity and selectivity. They bind their ligands using very different recognition modes involving a diverse set of hydrogen bonding and stacking interactions. Apparently, the potential structural diversity for guanine, guanosine, and guanine nucleotide binding motifs is far from being fully explored. Szostak and coworkers have derived a large set of different GTP-binding aptamer families differing widely in sequence, secondary structure, and ligand specificity. The so-called class V-GTP aptamer from this set binds GTP with very high affinity and has a complex secondary structure. Here we use solution NMR spectroscopy to demonstrate that the class V aptamer binds GTP through the formation of an intermolecular two-layered G-quadruplex structure that directly incorporates the ligand and folds only upon ligand addition. Ligand binding and G-quadruplex formation depend strongly on the identity of monovalent cations present with a clear preference for potassium ions. GTP binding through direct insertion into an intermolecular G-quadruplex is a previously unobserved structural variation for ligand-binding RNA motifs and rationalizes the previously observed specificity pattern of the class V aptamer for GTP analogs.

Ribosome biogenesis factor Tsr3 is the aminocarboxypropyl transferase responsible for 18S rRNA hypermodification in yeast and humans.

The chemically most complex modification in eukaryotic rRNA is the conserved hypermodified nucleotide N1-methyl-N3-aminocarboxypropyl-pseudouridine (m(1)acp(3)Ψ) located next to the P-site tRNA on the small subunit 18S rRNA. While S-adenosylmethionine was identified as the source of the aminocarboxypropyl (acp) group more than 40 years ago the enzyme catalyzing the acp transfer remained elusive. Here we identify the cytoplasmic ribosome biogenesis protein Tsr3 as the responsible enzyme in yeast and human cells. In functionally impaired Tsr3-mutants, a reduced level of acp modification directly correlates with increased 20S pre-rRNA accumulation. The crystal structure of archaeal Tsr3 homologs revealed the same fold as in SPOUT-class RNA-methyltransferases but a distinct SAM binding mode. This unique SAM binding mode explains why Tsr3 transfers the acp and not the methyl group of SAM to its substrate. Structurally, Tsr3 therefore represents a novel class of acp transferase enzymes.

NMR experiments for the rapid identification of P=O···H-X type hydrogen bonds in nucleic acids.

Hydrogen bonds involving the backbone phosphate groups occur with high frequency in functional RNA molecules. They are often found in well-characterized tertiary structural motifs presenting powerful probes for the rapid identification of these motifs for structure elucidation purposes. We have shown recently that stable hydrogen bonds to the phosphate backbone can in principle be detected by relatively simple NMR-experiments, providing the identity of both the donor hydrogen and the acceptor phosphorous within the same experiment (Duchardt-Ferner et al., Angew Chem Int Ed Engl 50:7927-7930, 2011). However, for imino and hydroxyl hydrogen bond donor groups rapidly exchanging with the solvent as well as amino groups broadened by conformational exchange experimental sensitivity is severely hampered by extensive line broadening. Here, we present improved methods for the rapid identification of hydrogen bonds to phosphate groups in nucleic acids by NMR. The introduction of the SOFAST technique into 1H,31P-correlation experiments as well as a BEST-HNP experiment exploiting 3hJN,P rather than 2hJH,P coupling constants enables the rapid and sensitive identification of these hydrogen bonds in RNA. The experiments are applicable for larger RNAs (up to ~ 100-nt), for donor groups influenced by conformational exchange processes such as amino groups and for hydrogen bonds with rather labile hydrogens such as 2'-OH groups as well as for moderate sample concentrations. Interestingly, the size of the through-hydrogen bond scalar coupling constants depends not only on the type of the donor group but also on the structural context. The largest coupling constants were measured for hydrogen bonds involving the imino groups of protonated cytosine nucleotides as donors.

Evaluation of 15N-detected H-N correlation experiments on increasingly large RNAs.

Recently, 15N-detected multidimensional NMR experiments have been introduced for the investigation of proteins. Utilization of the slow transverse relaxation of nitrogen nuclei in a 15N-TROSY experiment allowed recording of high quality spectra for high molecular weight proteins, even in the absence of deuteration. Here, we demonstrate the applicability of three 15N-detected H-N correlation experiments (TROSY, BEST-TROSY and HSQC) to RNA. With the newly established 15N-detected BEST-TROSY experiment, which proves to be the most sensitive 15N-detected H-N correlation experiment, spectra for five RNA molecules ranging in size from 5 to 100 kDa were recorded. These spectra yielded high resolution in the 15N-dimension even for larger RNAs since the increase in line width with molecular weight is more pronounced in the 1H- than in the 15N-dimension. Further, we could experimentally validate the difference in relaxation behavior of imino groups in AU and GC base pairs. Additionally, we showed that 15N-detected experiments theoretically should benefit from sensitivity and resolution advantages at higher static fields but that the latter is obscured by exchange dynamics within the RNAs.

A Stably Protonated Adenine Nucleotide with a Highly Shifted pKa Value Stabilizes the Tertiary Structure of a GTP-Binding RNA Aptamer.

RNA tertiary structure motifs are stabilized by a wide variety of hydrogen-bonding interactions. Protonated A and C nucleotides are normally not considered to be suitable building blocks for such motifs since their pKa values are far from physiological pH. Here, we report the NMR solution structure of an in vitro selected GTP-binding RNA aptamer bound to GTP with an intricate tertiary structure. It contains a novel kind of base quartet stabilized by a protonated A residue. Owing to its unique structural environment in the base quartet, the pKa value for the protonation of this A residue in the complex is shifted by more than 5 pH units compared to the pKa for A nucleotides in single-stranded RNA. This is the largest pKa shift for an A residue in structured nucleic acids reported so far, and similar in size to the largest pKa shifts observed for amino acid side chains in proteins. Both RNA pre-folding and ligand binding contribute to the pKa shift.