The cotranscriptional folding landscape for two cyclic di-nucleotide-sensing riboswitches with highly homologous aptamer domains acting either as ON- or OFF-switches.

Riboswitches are gene regulatory elements located in untranslated mRNA regions. They bind inducer molecules with high affinity and specificity. Cyclic-di-nucleotide-sensing riboswitches are major regulators of genes for the environment, membranes and motility (GEMM) of bacteria. Up to now, structural probing assays or crystal structures have provided insight into the interaction between cyclic-di-nucleotides and their corresponding riboswitches. ITC analysis, NMR analysis and computational modeling allowed us to gain a detailed understanding of the gene regulation mechanisms for the Cd1 (Clostridium difficile) and for the pilM (Geobacter metallireducens) riboswitches and their respective di-nucleotides c-di-GMP and c-GAMP. Binding capability showed a 25 nucleotide (nt) long window for pilM and a 61 nt window for Cd1. Within this window, binding affinities ranged from 35 µM to 0.25 µM spanning two orders of magnitude for Cd1 and pilM showing a strong dependence on competing riboswitch folds. Experimental results were incorporated into a Markov simulation to further our understanding of the transcriptional folding pathways of riboswitches. Our model showed the ability to predict riboswitch gene regulation and its dependence on transcription speed, pausing and ligand concentration.

RNA modifications stabilize the tertiary structure of tRNAfMet by locally increasing conformational dynamics.

A plethora of modified nucleotides extends the chemical and conformational space for natural occurring RNAs. tRNAs constitute the class of RNAs with the highest modification rate. The extensive modification modulates their overall stability, the fidelity and efficiency of translation. However, the impact of nucleotide modifications on the local structural dynamics is not well characterized. Here we show that the incorporation of the modified nucleotides in tRNAfMet from Escherichia coli leads to an increase in the local conformational dynamics, ultimately resulting in the stabilization of the overall tertiary structure. Through analysis of the local dynamics by NMR spectroscopic methods we find that, although the overall thermal stability of the tRNA is higher for the modified molecule, the conformational fluctuations on the local level are increased in comparison to an unmodified tRNA. In consequence, the melting of individual base pairs in the unmodified tRNA is determined by high entropic penalties compared to the modified. Further, we find that the modifications lead to a stabilization of long-range interactions harmonizing the stability of the tRNA's secondary and tertiary structure. Our results demonstrate that the increase in chemical space through introduction of modifications enables the population of otherwise inaccessible conformational substates.

Integrated NMR/Molecular Dynamics Determination of the Ensemble Conformation of a Thermodynamically Stable CUUG RNA Tetraloop.

Both experimental and theoretical structure determinations of RNAs have remained challenging due to the intrinsic dynamics of RNAs. We report here an integrated nuclear magnetic resonance/molecular dynamics (NMR/MD) structure determination approach to describe the dynamic structure of the CUUG tetraloop. We show that the tetraloop undergoes substantial dynamics, leading to averaging of the experimental data. These dynamics are particularly linked to the temperature-dependent presence of a hydrogen bond within the tetraloop. Interpreting the NMR data by a single structure represents the low-temperature structure well but fails to capture all conformational states occurring at a higher temperature. We integrate MD simulations, starting from structures of CUUG tetraloops within the Protein Data Bank, with an extensive set of NMR data, and provide a structural ensemble that describes the dynamic nature of the tetraloop and the experimental NMR data well. We thus show that one of the most stable and frequently found RNA tetraloops displays substantial dynamics, warranting such an integrated structural approach.

A Protonated Cytidine Stabilizes the Ligand-Binding Pocket in the PreQ1 Riboswitch in Thermophilic Bacteria.

Riboswitches are bacterial mRNA structure elements regulating either transcription or translation of downstream genes in response to high-affinity binding of a low molecular weight ligand. Among this diverse group of RNA structures, the class-I preQ1 sensing riboswitches (QSW) stand out since they are the smallest known natural riboswitches. The preQ1 sensing riboswitches combine ligand sensing and functional control within a single structural domain that adopts a pseudoknot conformation encapsulating both the cognate ligand and the ribosome binding site. preQ1 sensing riboswitches also occur in thermophilic bacteria. In these cases, their tertiary structures have to be stable even at temperatures above 60 °C to be functional at the organism's optimal growth temperatures. Despite the available high-resolution structures of these riboswitches, it is not yet understood which tertiary interactions are primarily responsible for their exceptional temperature stability. Here, we show that an intricate three-dimensional network of non-canonical interactions involving various non-neighboring nucleobases is the origin of the riboswitch's thermostability. An essential part of this network is a so far undetected stably protonated cytidine. It is characterized by an exceptional high pKA value of >9.7 and could be unambiguously identified through the application of modern heteronuclear detected NMR experiments. Thus, the presence or absence of a single proton can modulate the formation of an RNA tertiary structure and ligand binding capacity under extreme environmental conditions.

Modulation of Aß42 Aggregation Kinetics and Pathway by Low-Molecular-Weight Inhibitors.

The aggregation of amyloid-ß 42 (Aß42) is directly related to the pathogenesis of Alzheimer's disease. Here, we have investigated the early stages of the aggregation process, during which most of the cytotoxic species are formed. Aß42 aggregation kinetics, characterized by the quantification of Aß42 monomer consumption, were tracked by real-time solution NMR spectroscopy (RT-NMR) allowing the impact that low-molecular-weight (LMW) inhibitors and modulators exert on the aggregation process to be analysed. Distinct differences in the Aß42 kinetic profiles were apparent and were further investigated kinetically and structurally by using thioflavin T (ThT) and transmission electron microscopy (TEM), respectively. LMW inhibitors were shown to have a differential impact on early-state aggregation. Insight provided here could direct future therapeutic design based on kinetic profiling of the process of fibril formation.

Site-Specific Labeling of RNAs with Modified and 19 F-Labeled Nucleotides by Chemo-Enzymatic Synthesis.

More than 170 post-transcriptional modifications of RNAs have currently been identified. Detailed biophysical investigations of these modifications have been limited since large RNAs containing these post-transcriptional modifications are difficult to produce. Further, adequate readout of spectroscopic fingerprints are important, necessitating additional labeling procedures beyond the naturally occurring RNA modifications. Here, we report the chemo-enzymatic synthesis of RNA modifications and several structurally similar fluorine-modified analogs further optimizing a recently developed methodology.[1] This chemo-enzymatic method allows synthesis of also large RNAs. We were able to incorporate 16 modified nucleotides and 6 19 F-labeled nucleotides. To showcase the applicability of such modified large RNAs, we incorporated a 19 F-labeled cytidine into the aptamer domain of the 2'dG sensing riboswitch (2'dG-sw) from Mesoplasma florum, enabling characterizing RNA fold, ligand binding and kinetics. Thanks to the large chemical shift dispersion of 19 F, we can detect conformational heterogeneity in the apo state of the riboswitch.

The RNA chaperone StpA enables fast RNA refolding by destabilization of mutually exclusive base pairs within competing secondary structure elements.

In bacteria RNA gene regulatory elements refold dependent on environmental clues between two or more long-lived conformational states each associated with a distinct regulatory state. The refolding kinetics are strongly temperature-dependent and especially at lower temperatures they reach timescales that are biologically not accessible. To overcome this problem, RNA chaperones have evolved. However, the precise molecular mechanism of how these proteins accelerate RNA refolding reactions remains enigmatic. Here we show how the RNA chaperone StpA of Escherichia coli leads to an acceleration of a bistable RNA's refolding kinetics through the selective destabilization of key base pairing interactions. We find in laser assisted real-time NMR experiments on photocaged bistable RNAs that the RNA chaperone leads to a two-fold increase in refolding rates at low temperatures due to reduced stability of ground state conformations. Further, we can show that upon interaction with StpA, base pairing interactions in the bistable RNA are modulated to favor refolding through the dominant pseudoknotted transition pathway. Our results shed light on the molecular mechanism of the interaction between RNA chaperones and bistable RNAs and are the first step into a functional classification of chaperones dependent on their biophysical mode of operation.

NMR structure of the Vibrio vulnificus ribosomal protein S1 domains D3 and D4 provides insights into molecular recognition of single-stranded RNAs.

The ribosomal S1 protein (rS1) is indispensable for translation initiation in Gram-negative bacteria. rS1 is a multidomain protein that acts as an RNA chaperone and ensures that mRNAs can bind the ribosome in a single-stranded conformation, which could be related to fast recognition. Although many ribosome structures were solved in recent years, a high-resolution structure of a two-domain mRNA-binding competent rS1 construct is not yet available. Here, we present the NMR solution structure of the minimal mRNA-binding fragment of Vibrio Vulnificus rS1 containing the domains D3 and D4. Both domains are homologues and adapt an oligonucleotide-binding fold (OB fold) motif. NMR titration experiments reveal that recognition of miscellaneous mRNAs occurs via a continuous interaction surface to one side of these structurally linked domains. Using a novel paramagnetic relaxation enhancement (PRE) approach and exploring different spin-labeling positions within RNA, we were able to track the location and determine the orientation of the RNA in the rS1-D34 bound form. Our investigations show that paramagnetically labeled RNAs, spiked into unmodified RNA, can be used as a molecular ruler to provide structural information on protein-RNA complexes. The dynamic interaction occurs on a defined binding groove spanning both domains with identical ß2-ß3-ß5 interfaces. Evidently, the 3'-ends of the cis-acting RNAs are positioned in the direction of the N-terminus of the rS1 protein, thus towards the 30S binding site and adopt a conformation required for translation initiation.

A 300-fold enhancement of imino nucleic acid resonances by hyperpolarized water provides a new window for probing RNA refolding by 1D and 2D NMR.

NMR sensitivity-enhancement methods involving hyperpolarized water could be of importance for solution-state biophysical investigations. Hyperpolarized water (HyperW) can enhance the 1H NMR signals of exchangeable sites by orders of magnitude over their thermal counterparts, while providing insight into chemical exchange and solvent accessibility at a site-resolved level. As HyperW's enhancements are achieved by exploiting fast solvent exchanges associated with minimal interscan delays, possibilities for the rapid monitoring of chemical reactions and biomolecular (re)folding are opened. HyperW NMR can also accommodate heteronuclear transfers, facilitating the rapid acquisition of 2-dimensional (2D) 15N-1H NMR correlations, and thereby combining an enhanced spectral resolution with speed and sensitivity. This work demonstrates how these qualities can come together for the study of nucleic acids. HyperW injections were used to target the guanine-sensing riboswitch aptamer domain (GSRapt) of the xpt-pbuX operon in Bacillus subtilis Unlike what had been observed in proteins, where residues benefited of HyperW NMR only if/when sufficiently exposed to water, these enhancements applied to every imino resonance throughout the RNA. The >300-fold enhancements observed in the resulting 1H NMR spectra allowed us to monitor in real time the changes that GSRapt undergoes upon binding hypoxanthine, a high-affinity interaction leading to conformational refolding on a ∼1-s timescale at 36 °C. Structural responses could be identified for several nucleotides by 1-dimensional (1D) imino 1H NMR as well as by 2D HyperW NMR spectra acquired upon simultaneous injection of hyperpolarized water and hypoxanthine. The folding landscape revealed by this HyperW strategy for GSRapt, is briefly discussed.

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.

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.

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.

Refolding through a Linear Transition State Enables Fast Temperature Adaptation of a Translational Riboswitch.

The adenine-sensing riboswitch from the Gram-negative bacterium Vibrio vulnificus is an RNA-based gene regulatory element that acts in response to both its cognate low-molecular weight ligand and temperature. The combined sensitivity to environmental temperature and ligand concentration is maintained by an equilibrium of three distinct conformations involving two ligand-free states and one ligand-bound state. The key structural element that undergoes refolding in the ligand-free states comprises a 35-nucleotide temperature response module. Here, we present the structural characterization of this temperature response module. We employ high-resolution NMR spectroscopy and photocaged RNAs as molecular probes to decipher the kinetic and thermodynamic framework of the secondary structure transition in the apo state of the riboswitch. We propose a model for the transition state adopted during the thermal refolding of the temperature response module that connects two mutually exclusive long-lived and stable conformational states. This transition state is characterized by a comparatively low free activation enthalpy. A pseudoknot conformation in the transition state, as commonly seen in RNA refolding, is therefore unlikely. More likely, the transition state of the adenine-sensing riboswitch temperature response module features a linear conformation.

More than Proton Detection-New Avenues for NMR Spectroscopy of RNA.

Ribonucleic acid oligonucleotides (RNAs) play pivotal roles in cellular function (riboswitches), chemical biology applications (SELEX-derived aptamers), cell biology and biomedical applications (transcriptomics). Furthermore, a growing number of RNA forms (long non-coding RNAs, circular RNAs) but also RNA modifications are identified, showing the ever increasing functional diversity of RNAs. To describe and understand this functional diversity, structural studies of RNA are increasingly important. However, they are often more challenging than protein structural studies as RNAs are substantially more dynamic and their function is often linked to their structural transitions between alternative conformations. NMR is a prime technique to characterize these structural dynamics with atomic resolution. To extend the NMR size limitation and to characterize large RNAs and their complexes above 200 nucleotides, new NMR techniques have been developed. This Minireview reports on the development of NMR methods that utilize detection on low-γ nuclei (heteronuclei like 13 C or 15 N with lower gyromagnetic ratio than 1 H) to obtain unique structural and dynamic information for large RNA molecules in solution. Experiments involve through-bond correlations of nucleobases and the phosphodiester backbone of RNA for chemical shift assignment and make information on hydrogen bonding uniquely accessible. Previously unobservable NMR resonances of amino groups in RNA nucleobases are now detected in experiments involving conformational exchange-resistant double-quantum 1 H coherences, detected by 13 C NMR spectroscopy. Furthermore, 13 C and 15 N chemical shifts provide valuable information on conformations. All the covered aspects point to the advantages of low-γ nuclei detection experiments in RNA.

Quantitative modeling of the function of kinetically driven transcriptional riboswitches.

Riboswitches are cis-acting regulatory mRNA elements in bacteria, that modulate the expression of their associated genes in response to a cognate metabolite, operating either on the level of translation or transcription. Transcriptional riboswitches have to fold into functional structures as they are being synthesized and, only if transcription rates and ligand binding kinetics match, structured transcription intermediates are enabled to undergo ligand-dependent conformational refolding as a prerequisite for ligand-mediated gene expression. Therefore, transcription rates are of essential importance for functional riboswitch-mediated gene regulation. Here, we propose a generalized modeling framework for the kinetic mechanisms of transcriptional riboswitches. The formalism accommodates time-dependent transcription rates and changes of metabolite concentration and permits incorporation of variations in transcription rate depending on transcript length. We derive explicit analytical expressions for the fraction of transcripts that determine repression or activation of gene expression as a function of pause site location and its slowing down of transcription for the case of the (2'dG)-sensing riboswitch from Mesoplasma florum. Our modeling challenges the current view on the exclusive importance of metabolite binding to transcripts containing only the aptamer domain. Numerical simulations of transcription proceeding in a continuous manner under time-dependent changes of metabolite concentration further suggest that rapid modulations in concentration result in a reduced dynamic range for riboswitch function regardless of transcription rate, while a combination of slow modulations and small transcription rates ensures a wide range of finely tuneable regulatory outcomes.

Conformational switch in the ribosomal protein S1 guides unfolding of structured RNAs for translation initiation.

Initiation of bacterial translation requires that the ribosome-binding site in mRNAs adopts single-stranded conformations. In Gram-negative bacteria the ribosomal protein S1 (rS1) is a key player in resolving of structured elements in mRNAs. However, the exact mechanism of how rS1 unfolds persistent secondary structures in the translation initiation region (TIR) is still unknown. Here, we show by NMR spectroscopy that Vibrio vulnificus rS1 displays a unique architecture of its mRNA-binding domains, where domains D3 and D4 provide the mRNA-binding platform and cover the nucleotide binding length of the full-length rS1. D5 significantly increases rS1's chaperone activity, although it displays structural heterogeneity both in isolation and in presence of the other domains, albeit to varying degrees. The heterogeneity is induced by the switch between the two equilibrium conformations and is triggered by an order-to-order transition of two mutually exclusive secondary structures (ß-strand-to-α-helix) of the 'AERERI' sequence. The conformational switching is exploited for melting of structured 5'-UTR's, as the conformational heterogeneity of D5 can compensate the entropic penalty of complex formation. Our data thus provides a detailed understanding of the intricate coupling of protein and RNA folding dynamics enabling translation initiation of structured mRNAs.

Dynamics of Bacteriorhodopsin in the Dark-Adapted State from Solution Nuclear Magnetic Resonance Spectroscopy.

To achieve efficient proton pumping in the light-driven proton pump bacteriorhodopsin (bR), the protein must be tightly coupled to the retinal to rapidly convert retinal isomerization into protein structural rearrangements. Methyl group dynamics of bR embedded in lipid nanodiscs were determined in the dark-adapted state, and were found to be mostly well ordered at the cytosolic side. Methyl groups in the M145A mutant of bR, which displays only 10 % residual proton pumping activity, are less well ordered, suggesting a link between side-chain dynamics on the cytosolic side of the bR cavity and proton pumping activity. In addition, slow conformational exchange, attributed to low frequency motions of aromatic rings, was indirectly observed for residues on the extracellular side of the bR cavity. This may be related to reorganization of the water network. These observations provide a detailed picture of previously undescribed equilibrium dynamics on different time scales for ground-state bR.

Three-state mechanism couples ligand and temperature sensing in riboswitches.

Riboswitches are cis-acting gene-regulatory RNA elements that can function at the level of transcription, translation and RNA cleavage. The commonly accepted molecular mechanism for riboswitch function proposes a ligand-dependent conformational switch between two mutually exclusive states. According to this mechanism, ligand binding to an aptamer domain induces an allosteric conformational switch of an expression platform, leading to activation or repression of ligand-related gene expression. However, many riboswitch properties cannot be explained by a pure two-state mechanism. Here we show that the regulation mechanism of the adenine-sensing riboswitch, encoded by the add gene on chromosome II of the human Gram-negative pathogenic bacterium Vibrio vulnificus, is notably different from a two-state switch mechanism in that it involves three distinct stable conformations. We characterized the temperature and Mg(2+) dependence of the population ratios of the three conformations and the kinetics of their interconversion at nucleotide resolution. The observed temperature dependence of a pre-equilibrium involving two structurally distinct ligand-free conformations of the add riboswitch conferred efficient regulation over a physiologically relevant temperature range. Such robust switching is a key requirement for gene regulation in bacteria that have to adapt to environments with varying temperatures. The translational adenine-sensing riboswitch represents the first example, to our knowledge, of a temperature-compensated regulatory RNA element.

Ligand-modulated folding of the full-length adenine riboswitch probed by NMR and single-molecule FRET spectroscopy.

The full-length translation-regulating add adenine riboswitch (Asw) from Vibrio vulnificus has a more complex conformational space than its isolated aptamer domain. In addition to the predicted apo (apoA) and holo conformation that feature the conserved three-way junctional purine riboswitch aptamer, it adopts a second apo (apoB) conformation with a fundamentally different secondary structure. Here, we characterized the ligand-dependent conformational dynamics of the full-length add Asw by NMR and by single-molecule FRET (smFRET) spectroscopy. Both methods revealed an adenine-induced secondary structure switch from the apoB-form to the apoA-form that involves no tertiary structural interactions between aptamer and expression platform. This strongly suggests that the add Asw triggers translation by capturing the apoA-form secondary structure in the holo state. Intriguingly, NMR indicated a homogenous, docked aptamer kissing loop fold for apoA and holo, while smFRET showed persistent aptamer kissing loop docking dynamics between comparably stable, undocked and docked substates of the apoA and the holo conformation. Unraveling the folding of large junctional riboswitches thus requires the integration of complementary solution structural techniques such as NMR and smFRET.