NMR characterization and ligand binding site of the stem-loop 2 motif from the Delta variant of SARS-CoV-2.

The stem-loop 2 motif (s2m) in SARS-CoV-2 (SCoV-2) is located in the 3'-UTR. Although s2m has been reported to display characteristics of a mobile genomic element that might lead to an evolutionary advantage, its function has remained unknown. The secondary structure of the original SCoV-2 RNA sequence (Wuhan-Hu-1) was determined by NMR in late 2020, delineating the base-pairing pattern and revealing substantial differences in secondary structure compared to SARS-CoV-1 (SCoV-1). The existence of a single G29742-A29756 mismatch in the upper stem of s2m leads to its destabilization and impedes a complete NMR analysis. With Delta, a variant of concern has evolved with one mutation compared to the original sequence that replaces G29742 by U29742. We show here that this mutation results in a more defined structure at ambient temperature accompanied by a rise in melting temperature. Consequently, we were able to identify >90% of the relevant NMR resonances using a combination of selective RNA labeling and filtered 2D NOESY as well as 4D NMR experiments. We present a comprehensive NMR analysis of the secondary structure, (sub)nanosecond dynamics, and ribose conformation of s2m Delta based on heteronuclear 13C NOE and T 1 measurements and ribose carbon chemical shift-derived canonical coordinates. We further show that the G29742U mutation in Delta has no influence on the druggability of s2m compared to the Wuhan-Hu-1 sequence. With the assignment at hand, we identify the flexible regions of s2m as the primary site for small molecule binding.

The 5'-terminal stem-loop RNA element of SARS-CoV-2 features highly dynamic structural elements that are sensitive to differences in cellular pH.

We present the nuclear magnetic resonance spectroscopy (NMR) solution structure of the 5'-terminal stem loop 5_SL1 (SL1) of the SARS-CoV-2 genome. SL1 contains two A-form helical elements and two regions with non-canonical structure, namely an apical pyrimidine-rich loop and an asymmetric internal loop with one and two nucleotides at the 5'- and 3'-terminal part of the sequence, respectively. The conformational ensemble representing the averaged solution structure of SL1 was validated using NMR residual dipolar coupling (RDC) and small-angle X-ray scattering (SAXS) data. We show that the internal loop is the major binding site for fragments of low molecular weight. This internal loop of SL1 can be stabilized by an A12-C28 interaction that promotes the transient formation of an A+•C base pair. As a consequence, the pKa of the internal loop adenosine A12 is shifted to 5.8, compared to a pKa of 3.63 of free adenosine. Furthermore, applying a recently developed pH-differential mutational profiling (PD-MaP) approach, we not only recapitulated our NMR findings of SL1 but also unveiled multiple sites potentially sensitive to pH across the 5'-UTR of SARS-CoV-2.

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.

Cell-free transcription-translation system: a dual read-out assay to characterize riboswitch function.

Cell-free protein synthesis assays have become a valuable tool to understand transcriptional and translational processes. Here, we established a fluorescence-based coupled in vitro transcription-translation assay as a read-out system to simultaneously quantify mRNA and protein levels. We utilized the well-established quantification of the expression of shifted green fluorescent protein (sGFP) as a read-out of protein levels. In addition, we determined mRNA quantities using a fluorogenic Mango-(IV) RNA aptamer that becomes fluorescent upon binding to the fluorophore thiazole orange (TO). We utilized a Mango-(IV) RNA aptamer system comprising four subsequent Mango-(IV) RNA aptamer elements with improved sensitivity by building Mango arrays. The design of this reporter assay resulted in a sensitive read-out with a high signal-to-noise ratio, allowing us to monitor transcription and translation time courses in cell-free assays with continuous monitoring of fluorescence changes as well as snapshots of the reaction. Furthermore, we applied this dual read-out assay to investigate the function of thiamine-sensing riboswitches thiM and thiC from Escherichia coli and the adenine-sensing riboswitch ASW from Vibrio vulnificus and pbuE from Bacillus subtilis, which represent transcriptional and translational on- and off-riboswitches, respectively. This approach enabled a microplate-based application, a valuable addition to the toolbox for high-throughput screening of riboswitch function.

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.

m6A Methylation of Transcription Leader Sequence of SARS-CoV-2 Impacts Discontinuous Transcription of Subgenomic mRNAs.

The SARS-CoV-2 genome has been shown to be m6A methylated at several positions in vivo. Strikingly, a DRACH motif, the recognition motif for adenosine methylation, resides in the core of the transcriptional regulatory leader sequence (TRS-L) at position A74, which is highly conserved and essential for viral discontinuous transcription. Methylation at position A74 correlates with viral pathogenicity. Discontinuous transcription produces a set of subgenomic mRNAs that function as templates for translation of all structural and accessory proteins. A74 is base-paired in the short stem-loop structure 5'SL3 that opens during discontinuous transcription to form long-range RNA-RNA interactions with nascent (-)-strand transcripts at complementary TRS-body sequences. A74 can be methylated by the human METTL3/METTL14 complex in vitro. Here, we investigate its impact on the structural stability of 5'SL3 and the long-range TRS-leader:TRS-body duplex formation necessary for synthesis of subgenomic mRNAs of all four viral structural proteins. Methylation uniformly destabilizes 5'SL3 and long-range duplexes and alters their relative equilibrium populations, suggesting that the m6A74 modification acts as a regulator for the abundance of viral structural proteins due to this destabilization.

The COVID19-NMR Consortium: A Public Report on the Impact of this New Global Collaboration.

The outbreak of COVID-19 in December 2019 required the formation of international consortia for a coordinated scientific effort to understand and combat the virus. In this Viewpoint Article, we discuss how the NMR community has gathered to investigate the genome and proteome of SARS-CoV-2 and tested them for binding to low-molecular-weight binders. External factors including extended lockdowns due to the global pandemic character of the viral infection triggered the transition from locally focused collaborative research conducted within individual research groups to digital exchange formats for immediate discussion of unpublished results and data analysis, sample sharing, and coordinated research between more than 50 groups from 18 countries simultaneously. We discuss key lessons that might pertain after the end of the pandemic and challenges that we need to address.

Characterization of Structure and Dynamics of the Guanidine-II Riboswitch from Escherichia coli by NMR Spectroscopy and Small-Angle X-ray Scattering (SAXS).

Riboswitches are regulatory RNA elements that undergo functionally important allosteric conformational switching upon binding of specific ligands. The here investigated guanidine-II riboswitch binds the small cation, guanidinium, and forms a kissing loop-loop interaction between its P1 and P2 hairpins. We investigated the structural changes to support previous studies regarding the binding mechanism. Using NMR spectroscopy, we confirmed the structure as observed in crystal structures and we characterized the kissing loop interaction upon addition of Mg2+ and ligand for the riboswitch aptamer from Escherichia coli. We further investigated closely related mutant constructs providing further insight into functional differences between the two (different) hairpins P1 and P2. Formation of intermolecular interactions were probed by small-angle X-ray scattering (SAXS) and NMR DOSY data. All data are consistent and show the formation of oligomeric states of the riboswitch induced by Mg2+ and ligand binding.

Secondary structure determination of conserved SARS-CoV-2 RNA elements by NMR spectroscopy.

The current pandemic situation caused by the Betacoronavirus SARS-CoV-2 (SCoV2) highlights the need for coordinated research to combat COVID-19. A particularly important aspect is the development of medication. In addition to viral proteins, structured RNA elements represent a potent alternative as drug targets. The search for drugs that target RNA requires their high-resolution structural characterization. Using nuclear magnetic resonance (NMR) spectroscopy, a worldwide consortium of NMR researchers aims to characterize potential RNA drug targets of SCoV2. Here, we report the characterization of 15 conserved RNA elements located at the 5' end, the ribosomal frameshift segment and the 3'-untranslated region (3'-UTR) of the SCoV2 genome, their large-scale production and NMR-based secondary structure determination. The NMR data are corroborated with secondary structure probing by DMS footprinting experiments. The close agreement of NMR secondary structure determination of isolated RNA elements with DMS footprinting and NMR performed on larger RNA regions shows that the secondary structure elements fold independently. The NMR data reported here provide the basis for NMR investigations of RNA function, RNA interactions with viral and host proteins and screening campaigns to identify potential RNA binders for pharmaceutical intervention.

19 F NMR-Based Fragment Screening for 14 Different Biologically Active RNAs and 10 DNA and Protein Counter-Screens.

We report here the nuclear magnetic resonance 19 F screening of 14 RNA targets with different secondary and tertiary structure to systematically assess the druggability of RNAs. Our RNA targets include representative bacterial riboswitches that naturally bind with nanomolar affinity and high specificity to cellular metabolites of low molecular weight. Based on counter-screens against five DNAs and five proteins, we can show that RNA can be specifically targeted. To demonstrate the quality of the initial fragment library that has been designed for easy follow-up chemistry, we further show how to increase binding affinity from an initial fragment hit by chemistry that links the identified fragment to the intercalator acridine. Thus, we achieve low-micromolar binding affinity without losing binding specificity between two different terminator structures.

No Abuse Potential of Silexan in Healthy Recreational Drug Users: A Randomized Controlled Trial.

BACKGROUND: Silexan is a lavender essential oil with established anxiolytic and calming efficacy. Here we asked whether there is a potential for abuse in human patients. METHODS: We carried out a phase I abuse liability single-center, double-blind, 5-way crossover study in healthy users of recreational central nervous system depressants. They received single oral doses of 80 mg (therapeutic dose) and 640 mg Silexan, 2 mg and 4 mg lorazepam (active control) and placebo in randomized order, with 4- to 14-day washout periods between treatments. Pharmacodynamic measures included validated visual analogue scales assessing positive, negative, and sedative drug effects and balance of effects; a short form of the Addiction Research Center Inventory; and a drug similarity assessment. The primary outcome measure was the individual maximum value on the drug liking visual analogue scale during 24 hours post-dose. RESULTS: Forty participants were randomized and 34 were evaluable for pharmacodynamic outcomes. In intraindividual head-to-head comparisons of the drug liking visual analogue scale maximum value, both doses of Silexan were rated similar to placebo whereas differences were observed between Silexan and lorazepam and between placebo and lorazepam (P < .001). These data were supported by all secondary measures of positive drug effects and of balance of effects. Differences between placebo and both doses of Silexan were always negligible in magnitude. Moreover, Silexan showed no sedative effects and was not perceived to be similar to commonly used drugs that participants had used in the past. CONCLUSIONS: Silexan did not exhibit any abuse potential in a standard abuse potential detection screen study and is unlikely to be recreationally abused.

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.

Ligand binding to 2΄-deoxyguanosine sensing riboswitch in metabolic context.

The mfl-riboswitch is a transcriptional off-switch, which down-regulates expression of subunit ß of ribonucleotide reductase in Mesoplasma florum upon 2΄-deoxyguanosine binding. We characterized binding of 2΄-deoxyguanosine to the mfl-aptamer domain (WT aptamer) and a sequence-stabilized aptamer (MT aptamer) under in vitro and 'in-cell-like' conditions by isothermal titration calorimetry (ITC) and nuclear magnetic resonance (NMR) spectroscopy. 'In-cell-like' environment was simulated by Bacillus subtilis cell extract, in which both aptamers remained sufficiently stable to detect the resonances of structural elements and ligand binding in 2D NMR experiments. Under 'in-cell-like'-environment, (i) the WT aptamer bound the endogenous metabolite guanosine and (ii) 2΄-deoxyguanosine efficiently displaced guanosine from the WT aptamer. In contrast, MT aptamer exhibited moderate binding to 2΄-deoxyguanosine and weak binding to guanosine. NMR experiments indicated that binding of guanosine was not limited to the aptamer domain of the riboswitch but also the full-length mfl-riboswitch bound guanosine, impacting on the regulation efficiency of the riboswitch and hinting that, in addition to 2΄-deoxyguanosine, guanosine plays a role in riboswitch function in vivo. Reporter gene assays in B. subtilis demonstrated the regulation capacity of the WT aptamer, whereas the MT aptamer with lower affinity to 2΄-deoxyguanosine was not able to regulate gene expression.

NMR Structural Profiling of Transcriptional Intermediates Reveals Riboswitch Regulation by Metastable RNA Conformations.

Gene repression induced by the formation of transcriptional terminators represents a prime example for the coupling of RNA synthesis, folding, and regulation. In this context, mapping the changes in available conformational space of transcription intermediates during RNA synthesis is important to understand riboswitch function. A majority of riboswitches, an important class of small metabolite-sensing regulatory RNAs, act as transcriptional regulators, but the dependence of ligand binding and the subsequent allosteric conformational switch on mRNA transcript length has not yet been investigated. We show a strict fine-tuning of binding and sequence-dependent alterations of conformational space by structural analysis of all relevant transcription intermediates at single-nucleotide resolution for the I-A type 2'dG-sensing riboswitch from Mesoplasma florum by NMR spectroscopy. Our results provide a general framework to dissect the coupling of synthesis and folding essential for riboswitch function, revealing the importance of metastable states for RNA-based gene regulation.

NMR 1H,19F-based screening of the four stem-looped structure 5_SL1-SL4 located in the 5'-untranslated region of SARS-CoV 2 RNA.

Development of new antiviral medication against the beta-coronavirus SARS-CoV-2 (SCoV2) is actively being pursued. Both NMR spectroscopy and crystallography as structural screening technologies have been utilised to screen the viral proteome for binding to fragment libraries. Here, we report on NMR screening of elements of the viral RNA genome with two different ligand libraries using 1H-NMR-screening experiments and 1H and 19F NMR-screening experiments for fluorinated compounds. We screened against the 5'-terminal 119 nucleotides located in the 5'-untranslated region of the RNA genome of SCoV2 and further dissected the four stem-loops into its constituent RNA elements to test specificity of binding of ligands to shorter and longer viral RNA stretches. The first library (DRTL-F library) is enriched in ligands binding to RNA motifs, while the second library (DSI-poised library) represents a fragment library originally designed for protein screening. Conducting screens with two different libraries allows us to compare different NMR screening methodologies, describe NMR screening workflows, validate the two different fragment libraries, and derive initial leads for further downstream medicinal chemistry optimisation.

Structure and dynamics of the deoxyguanosine-sensing riboswitch studied by NMR-spectroscopy.

The mfl-riboswitch regulates expression of ribonucleotide reductase subunit in Mesoplasma florum by binding to 2'-deoxyguanosine and thereby promoting transcription termination. We characterized the structure of the ligand-bound aptamer domain by NMR spectroscopy and compared the mfl-aptamer to the aptamer domain of the closely related purine-sensing riboswitches. We show that the mfl-aptamer accommodates the extra 2'-deoxyribose unit of the ligand by forming a more relaxed binding pocket than these found in the purine-sensing riboswitches. Tertiary structures of the xpt-aptamer bound to guanine and of the mfl-aptamer bound to 2'-deoxyguanosine exhibit very similar features, although the sequence of the mfl-aptamer contains several alterations compared to the purine-aptamer consensus sequence. These alterations include the truncation of a hairpin loop which is crucial for complex formation in all purine-sensing riboswitches characterized to date. We further defined structural features and ligand binding requirements of the free mfl-aptamer and found that the presence of Mg(2+) is not essential for complex formation, but facilitates ligand binding by promoting pre-organization of key structural motifs in the free aptamer.

Influence of ground-state structure and Mg2+ binding on folding kinetics of the guanine-sensing riboswitch aptamer domain.

Riboswitch RNAs fold into complex tertiary structures upon binding to their cognate ligand. Ligand recognition is accomplished by key residues in the binding pocket. In addition, it often crucially depends on the stability of peripheral structural elements. The ligand-bound complex of the guanine-sensing riboswitch from Bacillus subtilis, for example, is stabilized by extensive interactions between apical loop regions of the aptamer domain. Previously, we have shown that destabilization of this tertiary loop-loop interaction abrogates ligand binding of the G37A/C61U-mutant aptamer domain (Gsw(loop)) in the absence of Mg(2+). However, if Mg(2+) is available, ligand-binding capability is restored by a population shift of the ground-state RNA ensemble toward RNA conformations with pre-formed loop-loop interactions. Here, we characterize the striking influence of long-range tertiary structure on RNA folding kinetics and on ligand-bound complex structure, both by X-ray crystallography and time-resolved NMR. The X-ray structure of the ligand-bound complex reveals that the global architecture is almost identical to the wild-type aptamer domain. The population of ligand-binding competent conformations in the ground-state ensemble of Gsw(loop) is tunable through variation of the Mg(2+) concentration. We quantitatively describe the influence of distinct Mg(2+) concentrations on ligand-induced folding trajectories both by equilibrium and time-resolved NMR spectroscopy at single-residue resolution.

Mapping the landscape of RNA dynamics with NMR spectroscopy.

Among the three major classes of biomacromolecules (DNA, RNA, and proteins) RNA's pronounced dynamics are the most explicitly linked to its wide variety of functions, which include catalysis and the regulation of transcription, translation, and splicing. These functions are mediated by a range of RNA biomachinery, including such varied examples as macromolecular noncoding RNAs, microRNAs, small interfering RNAs, riboswitch RNAs, and RNA thermometers. In each case, the functional dynamics of an interconversion is characterized by an associated rate constant. In this Account, we provide an introduction to NMR spectroscopic characterization of the landscape of RNA dynamics. We introduce strategies for measuring NMR parameters at various time scales as well as the underlying models for describing the corresponding rate constants. RNA exhibits significant dynamic motion, which can be modulated by (i) intermolecular interactions, including specific and nonspecific binding of ions (such as Mg(2+) and tertiary amines), (ii) metabolites in riboswitches or RNA aptamers, and (iii) macromolecular interactions within ribonucleic protein particles, including the ribosome and the spliceosome. Our understanding of the nature of these dynamic changes in RNA targets is now being incorporated into RNA-specific approaches in the design of RNA inhibitors. Interactions of RNA with proteins, other RNAs, or small molecules often occur through binding mechanisms that follow an induced fit mechanism or a conformational selection mechanism, in which one of several populated RNA conformations is selected through ligand binding. The extent of functional dynamics, including the kinetic formation of a specific RNA tertiary fold, is dependent on the messenger RNA (mRNA) chain length. Thus, during de novo synthesis of mRNA, both in prokaryotes and eukaryotes, nascent mRNA of various lengths will adopt different secondary and tertiary structures. The speed of transcription has a critical influence on the functional dynamics of the RNA being synthesized. In addition to modulating the local dynamics of a conformational RNA ensemble, a given RNA sequence may adopt more than one global, three-dimensional structure. RNA modification is one way to select among these alternative structures, which are often characterized by nearly equal stability, but with high energy barriers for conformational interconversion. The refolding of different secondary and tertiary structures has been found to be a major regulatory mechanism for transcription and translation. These conformational transitions can be characterized with NMR spectroscopy, for any given RNA sequence, in response to external stimuli.

Mechanisms for differentiation between cognate and near-cognate ligands by purine riboswitches.

Riboswitches are elements in the 5'-untranslated region of mRNAs that regulate gene expression by directly interacting with metabolites related to their own gene products. A remarkable feature of this gene regulation mechanism is the high specificity of riboswitches for their cognate ligands. In this study, we used a combination of static and time-resolved NMR-spectroscopic methods to investigate the mechanisms for ligand specificity in purine riboswitches. We investigate the xpt-aptamer domain from a guanine-responsive riboswitch and the mfl-aptamer domain from a 2'-deoxyguanosine-responsive riboswitch. The xpt-aptamer binds the purine nucleobases guanine/hypoxanthine with high affinity, but, unexpectedly, also the nucleoside 2'-deoxyguanosine. On the other hand, the mfl-aptamer is highly specific for its cognate ligand 2'-deoxyguanosine, and does not bind purine ligands. We addressed the question of aptamer`s ligand specificity by real-time NMR spectroscopy. Our studies of ligand binding and subsequently induced aptamer folding revealed that the xpt-aptamer discriminates against non-cognate ligands by enhanced life-times of the cognate complex compared with non-cognate complexes, whereas the mfl-aptamer rejects non-cognate ligands at the level of ligand association, employing a kinetic proofreading mechanism.