Local-to-global signal transduction at the core of a Mn2+ sensing riboswitch

Foundation model for efficient biological discovery in single-molecule data

Modern data-intensive techniques offer ever deeper insights into biology, but render the process of discovery increasingly complex. For example, exploiting the unique ability of single-molecule fluorescence microscopy (SMFM)1-5. to uncover rare but critical intermediates often demands manual inspection of time traces and iterative ad hoc approaches that are difficult to systematize. To facilitate systematic and efficient discovery from SMFM data, we introduce META-SiM, a transformer-based foundation model pre-trained on diverse SMFM analysis tasks. META-SiM achieves high performance-rivaling best-in-class algorithms-on a broad range of analysis tasks including trace selection, classification, segmentation, idealization, and stepwise photobleaching analysis. Additionally, the model produces high-dimensional embedding vectors that encapsulate detailed information about each trace, which the web-based META-SiM Projector (https://www.simol-projector.org) casts into lower-dimensional space for efficient whole-dataset visualization, labeling, comparison, and sharing. Combining this Projector with the objective metric of Local Shannon Entropy enables rapid identification of condition-specific behaviors, even if rare or subtle. As a result, by applying META-SiM to an existing single-molecule Förster resonance energy transfer (smFRET) dataset6, we discover a previously unobserved intermediate state in pre-mRNA splicing. META-SiM thus removes bottlenecks, improves objectivity, and both systematizes and accelerates biological discovery in complex single-molecule data.

FRET-guided modeling of nucleic acids

The functional diversity of RNAs is encoded in their innate conformational heterogeneity. The combination of single-molecule spectroscopy and computational modeling offers new attractive opportunities to map structural transitions within nucleic acid ensembles. Here, we describe a framework to harmonize single-molecule Förster resonance energy transfer (FRET) measurements with molecular dynamics simulations and de novo structure prediction. Using either all-atom or implicit fluorophore modeling, we recreate FRET experiments in silico, visualize the underlying structural dynamics and quantify the reaction coordinates. Using multiple accessible-contact volumes as a post hoc scoring method for fragment assembly in Rosetta, we demonstrate that FRET can be used to filter a de novo RNA structure prediction ensemble by refuting models that are not compatible with in vitro FRET measurement. We benchmark our FRET-assisted modeling approach on double-labeled DNA strands and validate it against an intrinsically dynamic manganese(II)-binding riboswitch. We show that a FRET coordinate describing the assembly of a four-way junction allows our pipeline to recapitulate the global fold of the riboswitch displayed by the crystal structure. We conclude that computational fluorescence spectroscopy facilitates the interpretability of dynamic structural ensembles and improves the mechanistic understanding of nucleic acid interactions.

A nascent riboswitch helix orchestrates robust transcriptional regulation through signal integration

Widespread manganese-sensing transcriptional riboswitches effect the dependable gene regulation needed for bacterial manganese homeostasis in changing environments. Riboswitches - like most structured RNAs - are believed to fold co-transcriptionally, subject to both ligand binding and transcription events; yet how these processes are orchestrated for robust regulation is poorly understood. Through a combination of single-molecule and bulk approaches, we discover how a single Mn2+ ion and the transcribing RNA polymerase (RNAP), paused immediately downstream by a DNA template sequence, are coordinated by the bridging switch helix P1.1 in the representative Lactococcus lactis riboswitch. This coordination achieves a heretofore-overlooked semi-docked global conformation of the nascent RNA, P1.1 base pair stabilization, transcription factor NusA ejection, and RNAP pause extension, thereby enforcing transcription readthrough. Our work demonstrates how a central, adaptable RNA helix functions analogous to a molecular fulcrum of a first-class lever system to integrate disparate signals for finely balanced gene expression control.

A nascent riboswitch helix orchestrates robust transcriptional regulation through signal integration

Widespread manganese-sensing transcriptional riboswitches effect the dependable gene regulation needed for bacterial manganese homeostasis in changing environments. Riboswitches - like most structured RNAs - are believed to fold co-transcriptionally, subject to both ligand binding and transcription events; yet how these processes are orchestrated for robust regulation is poorly understood. Through a combination of single molecule and bulk approaches, we discovered how a single Mn 2+ ion and the transcribing RNA polymerase (RNAP), paused immediately downstream by a DNA template sequence, are coordinated by the bridging switch helix P1.1 in the paradigmatic Lactococcus lactis riboswitch. This coordination achieves a heretofore-overlooked semi-docked global conformation of the nascent RNA, P1.1 base pair stabilization, transcription factor NusA ejection, and RNAP pause extension, thereby enforcing transcription readthrough. Our work demonstrates how a central, adaptable RNA helix functions analogous to a molecular fulcrum of a first-class lever system to integrate disparate signals for finely balanced gene expression control.

Mechanism of Fluoride Ion Encapsulation by Magnesium Ions in a Bacterial Riboswitch

Riboswitches sense various ions in bacteria and activate gene expression to synthesize proteins that help maintain ion homeostasis. The crystal structure of the aptamer domain (AD) of the fluoride riboswitch shows that the F- ion is encapsulated by three Mg2+ ions bound to the ligand-binding domain (LBD) located at the core of the AD. The assembly mechanism of this intricate structure is unknown. To this end, we performed computer simulations using coarse-grained and all-atom RNA models to bridge multiple time scales involved in riboswitch folding and ion binding. We show that F- encapsulation by the Mg2+ ions bound to the riboswitch involves multiple sequential steps. Broadly, two Mg2+ ions initially interact with the phosphate groups of the LBD using water-mediated outer-shell coordination and transition to a direct inner-shell interaction through dehydration to strengthen their interaction with the LBD. We propose that the efficient binding mode of the third Mg2+ and F- is that they form a water-mediated ion pair and bind to the LBD simultaneously to minimize the electrostatic repulsion between three Mg2+ bound to the LBD. The tertiary stacking interactions among the LBD nucleobases alone are insufficient to stabilize the alignment of the phosphate groups to facilitate Mg2+ binding. We show that the stability of the whole assembly is an intricate balance of the interactions among the five phosphate groups, three Mg2+, and the encapsulated F- ion aided by the Mg2+ solvated water. These insights are helpful in the rational design of RNA-based ion sensors and fast-switching logic gates.

Direct observation of tRNA-chaperoned folding of a dynamic mRNA ensemble

T-box riboswitches are multi-domain noncoding RNAs that surveil individual amino acid availabilities in most Gram-positive bacteria. T-boxes directly bind specific tRNAs, query their aminoacylation status to detect starvation, and feedback control the transcription or translation of downstream amino-acid metabolic genes. Most T-boxes rapidly recruit their cognate tRNA ligands through an intricate three-way stem I-stem II-tRNA interaction, whose establishment is not understood. Using single-molecule FRET, SAXS, and time-resolved fluorescence, we find that the free T-box RNA assumes a broad distribution of open, semi-open, and closed conformations that only slowly interconvert. tRNA directly binds all three conformers with distinct kinetics, triggers nearly instantaneous collapses of the open conformations, and returns the T-box RNA to their pre-binding conformations upon dissociation. This scissors-like dynamic behavior is enabled by a hinge-like pseudoknot domain which poises the T-box for rapid tRNA-induced domain closure. This study reveals tRNA-chaperoned folding of flexible, multi-domain mRNAs through a Venus flytrap-like mechanism.

Discovering riboswitches: the past and the future

Riboswitches are structured noncoding RNA domains used by many bacteria to monitor the concentrations of target ligands and regulate gene expression accordingly. In the past 20 years over 55 distinct classes of natural riboswitches have been discovered that selectively sense small molecules or elemental ions, and thousands more are predicted to exist. Evidence suggests that some riboswitches might be direct descendants of the RNA-based sensors and switches that were likely present in ancient organisms before the evolutionary emergence of proteins. We provide an overview of the current state of riboswitch research, focusing primarily on the discovery of riboswitches, and speculate on the major challenges facing researchers in the field.

Probing Transient Riboswitch Structures via Single Molecule Accessibility Analysis

Riboswitches are a class of RNA motifs in the untranslated regions of bacterial messenger RNAs (mRNAs) that can adopt different conformations to regulate gene expression. The binding of specific small molecule or ion ligands, or other RNAs, influences the conformation the riboswitch adopts. Single Molecule Kinetic Analysis of RNA Transient Structure (SiM-KARTS) offers an approach for probing this structural isomerization, or conformational switching, at the level of single mRNA molecules. SiM-KARTS utilizes fluorescently labeled, short, sequence-complementary DNA or RNA oligonucleotide probes that transiently access a specific RNA conformation over another. Binding and dissociation to a surface-immobilized target RNA of arbitrary length are monitored by Total Internal Reflection Fluorescence Microscopy (TIRFM) and quantitatively analyzed, via spike train and burst detection, to elucidate the rate constants of isomerization, revealing mechanistic insights into riboswitching.

Ligand-Dependent Volumetric Characterization of Manganese Riboswitch Folding: A High-Pressure Single-Molecule Kinetic Study

Nanoscopic differences in free volume result in pressure-dependent changes in free energies which can therefore impact folding/unfolding stability of biomolecules. Although such effects are typically insignificant under ambient pressure conditions, they are crucially important for deep ocean marine life, where the hydraulic pressure can be on the kilobar scale. In this work, single molecule FRET spectroscopy is used to study the effects of pressure on both the kinetics and overall thermodynamics for folding/unfolding of the manganese riboswitch. Detailed pressure-dependent analysis of the conformational kinetics allows one to extract precision changes (σ ≲ 4-8 ų) in free volumes not only between the fully folded/unfolded conformations but also with respect to the folding transition state of the manganese riboswitch. This permits first extraction of a novel "reversible work" free energy (PΔV) landscape, which reveals a monotonic increase in manganese riboswitch volume along the folding coordinate. Furthermore, such a tool permits exploration of pressure-dependent effects on both Mn²⁺ binding and riboswitch folding, which demonstrate that ligand attachment stabilizes the riboswitch under pressure by decreasing the volume increase upon folding (ΔΔV < 0). Such competition between ligand binding and pressure-induced denaturation dynamics could be of significant evolutionary advantage, compensating for a weakening in riboswitch tertiary structure with pressure-mediated ligand binding and promotion of folding response.

Potential effects of metal ion induced two-state allostery on the regulatory mechanism of add adenine riboswitch

Riboswitches normally regulate gene expression through structural changes in response to the specific binding of cellular metabolites or metal ions. Taking add adenine riboswitch as an example, we explore the influences of metal ions (especially for K⁺ and Mg²⁺ ions) on the structure and dynamics of riboswitch aptamer (with and without ligand) by using molecular dynamic (MD) simulations. Our results show that a two-state transition marked by the structural deformation at the connection of J12 and P1 (C_J12-P1) is not only related to the binding of cognate ligands, but also strongly coupled with the change of metal ion environments. Moreover, the deformation of the structure at C_J12-P1 can be transmitted to P1 directly connected to the expression platform in multiple ways, which will affect the structure and stability of P1 to varying degrees, and finally change the regulation state of this riboswitch.

Molecular dynamic simulations are employed to assess the influence of metal ions on riboswitch structure and dynamics, suggesting a conformational control of riboswitch aptamers by metal ions before ligand binding.

Precise tuning of bacterial translation initiation by non-equilibrium 5'-UTR unfolding observed in single mRNAs

Noncoding, structured 5′-untranslated regions (5′-UTRs) of bacterial messenger RNAs (mRNAs) can control translation efficiency by forming structures that either recruit or repel the ribosome. Here we exploit a 5′-UTR embedded preQ₁-sensing, pseudoknotted translational riboswitch to probe how binding of a small ligand controls recruitment of the bacterial ribosome to the partially overlapping Shine-Dalgarno (SD) sequence. Combining single-molecule fluorescence microscopy with mutational analyses, we find that the stability of 30S ribosomal subunit binding is inversely correlated with the free energy needed to unfold the 5′-UTR during mRNA accommodation into the mRNA binding cleft. Ligand binding to the riboswitch stabilizes the structure to both antagonize 30S recruitment and accelerate 30S dissociation. Proximity of the 5′-UTR and stability of the SD:anti-SD interaction both play important roles in modulating the initial 30S-mRNA interaction. Finally, depletion of small ribosomal subunit protein S1, known to help resolve structured 5′-UTRs, further increases the energetic penalty for mRNA accommodation. The resulting model of rapid standby site exploration followed by gated non-equilibrium unfolding of the 5′-UTR during accommodation provides a mechanistic understanding of how translation efficiency is governed by riboswitches and other dynamic structure motifs embedded upstream of the translation initiation site of bacterial mRNAs.

Dynamics of metal binding and mutation in yybP-ykoY riboswitch of Lactococcus lactis

Riboswitch is a regulatory segment of messenger RNA (mRNA), which by binding to various cellular metabolites regulates the activity of mRNA via modulating transcription, translation, alternative splicing, and stability of the mRNA. yybP–ykoY riboswitch of Lactococcus lactis, which is present upstream of the yoaB gene, functions as a Mn²⁺-specific genetic ON-switch, and modulates expression of proteins which are significant for Mn²⁺ homeostasis. The P1.1 switch helix of the aptamer domain of the riboswitch contains an intrinsic transcription terminator structure, which gets stabilized with Mn²⁺ binding and causes disruption of terminator structure and allows the continuation of transcription. The current research work involved the evaluation of structural and dynamical properties of the yybP-ykoY riboswitch of L. lactis in its Mn²⁺-free, Mn²⁺-bound (wild-type), and Mn²⁺-bound mutant (A41U) states by applying molecular dynamics simulations. Based on the simulations, the effects of Mn²⁺ absence and A41U mutation were evaluated on the structure and dynamics of the riboswitches followed by the computation of the free energy of metal binding in the wild-type and the mutant riboswitches. The simulation results provided insights into the properties of the riboswitch with the focus on the dynamics of the P1.1 switch helix, and the manganese binding site designated as M_B site, as well as the relative stability of the wild-type and the mutant riboswitches, which helped to understand the structural and dynamical role of the metal ion involved in the function of Mn²⁺-sensing riboswitch.

The current research work involved the evaluation of structural and dynamical properties of yybP–ykoY riboswitch of L. lactis in Mn²⁺-free, Mn²⁺-bound (wild-type), and Mn²⁺-bound mutant (A41U) states by applying molecular dynamics simulations.

Refining the RNA Force Field with Small-Angle X-ray Scattering of Helix-Junction-Helix RNA

The growing recognition of the functional and therapeutic roles played by RNA and the difficulties in gaining atomic-level insights by experiments are paving the way for all-atom simulations of RNA. One of the main impediments to the use of all-atom simulations is the imbalance between the energy terms of the RNA force fields. Through exhaustive sampling of an RNA helix-junction-helix (HJH) model using enhanced sampling, we critically assessed the select Amber force fields against small-angle X-ray scattering (SAXS) experiments. The tested AMBER99SB, DES-AMBER, and CUFIX force fields show deviations from measured profiles. First, we identified parameters leading to inconsistencies. Then, as a way to balance the forces governing RNA folding, we adopted strategies to refine hydrogen bonding, backbone, and base-stacking parameters. We validated the modified force field (HB-CUFIX) against SAXS data of the HJH model in different ionic strengths. Moreover, we tested a set of independent RNA systems to cross-validate the force field. Overall, HB-CUFIX demonstrates improved performance in studying thermodynamics and structural properties of realistic RNA motifs.

An anionic ligand snap-locks a long-range interaction in a magnesium-folded riboswitch

The archetypical transcriptional crcB fluoride riboswitch from Bacillus cereus is an intricately structured non-coding RNA element enhancing gene expression in response to toxic levels of fluoride. Here, we used single molecule FRET to uncover three dynamically interconverting conformations appearing along the transcription process: two distinct undocked states and one pseudoknotted docked state. We find that the fluoride anion specifically snap-locks the magnesium-induced, dynamically docked state. The long-range, nesting, single base pair A40-U48 acts as the main linchpin, rather than the multiple base pairs comprising the pseudoknot. We observe that the proximally paused RNA polymerase further fine-tunes the free energy to promote riboswitch docking. Finally, we show that fluoride binding at short transcript lengths is an early step toward partitioning folding into the docked conformation. These results reveal how the anionic fluoride ion cooperates with the magnesium-associated RNA to govern regulation of downstream genes needed for fluoride detoxification of the cell.

Chemical Dual End-Labeling of Large Ribozymes

Fast and efficient site-specific labeling of long RNAs is one of the main bottlenecks limiting distance measurements by means of Förster resonance energy transfer (FRET) or electron paramagnetic resonance (EPR) spectroscopy. Here, we present an optimized protocol for dual end-labeling with different fluorophores at the same time meeting the restrictions of highly labile and degradation-sensitive RNAs. We describe in detail the dual-labeling of a catalytically active wild-type group II intron as a typical representative of long functional RNAs. The modular procedure chemically activates the 5'-phosphate and the 3'-ribose for bioconjugation with a pair of fluorophores, as shown herein, or with spin labels. The mild reaction conditions preserve the structural and functional integrity of the biomacromolecule and results in covalent, dual-labeled RNA in its pre-catalytic state in yields suitable for both ensemble and single-molecule FRET experiments.

The International Society of RNA Nanotechnology and Nanomedicine (ISRNN): The Present and Future of the Burgeoning Field

The International Society of RNA Nanotechnology and Nanomedicine (ISRNN) hosts an annual meeting series focused on presenting the latest research achievements involving RNA-based therapeutics and strategies, aiming to expand their current biomedical applications while overcoming the remaining challenges of the burgeoning field of RNA nanotechnology. The most recent online meeting hosted a series of engaging talks and discussions from an international cohort of leading nanotechnologists that focused on RNA modifications and modulation, dynamic RNA structures, overcoming delivery limitations using a variety of innovative platforms and approaches, and addressing the newly explored potential for immunomodulation with programmable nucleic acid nanoparticles. In this Nano Focus, we summarize the main discussion points, conclusions, and future directions identified during this two-day webinar as well as more recent advances to highlight and to accelerate this exciting field.

Fluorogenic aptamers resolve the flexibility of RNA junctions using orientation-dependent FRET

To further understand the transcriptome, new tools capable of measuring folding, interactions, and localization of RNA are needed. Although Förster resonance energy transfer (FRET) is an angle- and distance-dependent phenomenon, the majority of FRET measurements have been used to report distances, by assuming rotationally averaged donor-acceptor pairs. Angle-dependent FRET measurements have proven challenging for nucleic acids due to the difficulties in incorporating fluorophores rigidly into local substructures in a biocompatible manner. Fluorescence turn-on RNA aptamers are genetically encodable tags that appear to rigidly confine their cognate fluorophores, and thus have the potential to report angular-resolved FRET. Here, we use the fluorescent aptamers Broccoli and Mango-III as donor and acceptor, respectively, to measure the angular dependence of FRET. Joining the two fluorescent aptamers by a helix of variable length allowed systematic rotation of the acceptor fluorophore relative to the donor. FRET oscillated in a sinusoidal manner as a function of helix length, consistent with simulated data generated from models of oriented fluorophores separated by an inflexible helix. Analysis of the orientation dependence of FRET allowed us to demonstrate structural rigidification of the NiCo riboswitch upon transition metal-ion binding. This application of fluorescence turn-on aptamers opens the way to improved structural interpretation of ensemble and single-molecule FRET measurements of RNA.

Transcriptional Riboswitches Integrate Timescales for Bacterial Gene Expression Control

Transcriptional riboswitches involve RNA aptamers that are typically found in the 5' untranslated regions (UTRs) of bacterial mRNAs and form alternative secondary structures upon binding to cognate ligands. Alteration of the riboswitch's secondary structure results in perturbations of an adjacent expression platform that controls transcription elongation and termination, thus turning downstream gene expression "on" or "off." Riboswitch ligands are typically small metabolites, divalent cations, anions, signaling molecules, or other RNAs, and can be part of larger signaling cascades. The interconnectedness of ligand binding, RNA folding, RNA transcription, and gene expression empowers riboswitches to integrate cellular processes and environmental conditions across multiple timescales. For a successful response to an environmental cue that may determine a bacterium's chance of survival, a coordinated coupling of timescales from microseconds to minutes must be achieved. This review focuses on recent advances in our understanding of how riboswitches affect such critical gene expression control across time.

Automatic classification and segmentation of single-molecule fluorescence time traces with deep learning

Traces from single-molecule fluorescence microscopy (SMFM) experiments exhibit photophysical artifacts that typically necessitate human expert screening, which is time-consuming and introduces potential for user-dependent expectation bias. Here, we use deep learning to develop a rapid, automatic SMFM trace selector, termed AutoSiM, that improves the sensitivity and specificity of an assay for a DNA point mutation based on single-molecule recognition through equilibrium Poisson sampling (SiMREPS). The improved performance of AutoSiM is based on accepting both more true positives and fewer false positives than the conventional approach of hidden Markov modeling (HMM) followed by hard thresholding. As a second application, the selector is used for automated screening of single-molecule Förster resonance energy transfer (smFRET) data to identify high-quality traces for further analysis, and achieves ~90% concordance with manual selection while requiring less processing time. Finally, we show that AutoSiM can be adapted readily to novel datasets, requiring only modest Transfer Learning.

Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics

Nucleic acids play critical roles in carrying genetic information, participating in catalysis, and preserving chromosomal structure. Despite over a century of study, efforts to understand the dynamics and structure-function relationships of DNA and RNA at the atomic level are still ongoing. Molecular dynamics (MD) simulations augment experiments by providing atomistic resolution and quantitative relationships between structure and conformational energy. Steady advancements in computer hardware, software, and atomistic force fields (FFs) over 40 years have facilitated new discoveries. Here, we review nucleic acid FF development with emphasis on recent refinements that have improved descriptions of important nucleic acid properties. We then discuss several key examples of successes and challenges in modeling nucleic acid structure and dynamics using the latest FFs.