Interrogating RNA-Small Molecule Interactions with Structure Probing and Artificial Intelligence-Augmented Molecular Simulations

Mechanistic Studies of Small Molecule Ligands Selective to RNA Single G Bulges

Small-molecule RNA binders have emerged as an important pharmacological modality. A profound understanding of the ligand selectivity, binding mode, and influential factors governing ligand engagement with RNA targets is the foundation for rational ligand design. Here, we report a novel class of coumarin derivatives exhibiting selective binding affinity towards single G RNA bulges. Harnessing the computational power of all-atom Gaussian accelerated Molecular Dynamics (GaMD) simulations, we unveiled a rare minor groove binding mode of the ligand with a key interaction between the coumarin moiety and the G bulge. This predicted binding mode is consistent with results obtained from structure-activity-relationship (SAR) studies and transverse relaxation measurements by NMR spectroscopy. We further generated 444 molecular descriptors from 69 coumarin derivatives and identified key contributors to the binding events, such as charge state and planarity, by lasso (least absolute shrinkage and selection operator) regression. Strikingly, small structure perturbations on these key contributors, such as the addition of a methyl group that disrupts the planarity of the ligand resulted in > 100-fold reduction in the binding affinity. Our work deepened the understanding of RNA-small molecule interactions and integrated a new generalizable platform for the rational design of selective small-molecule RNA binders.

Mechanistic analysis of Riboswitch Ligand interactions provides insights into pharmacological control over gene expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Machine learning-augmented molecular dynamics simulations (MD) reveal insights into the disconnect between affinity and activation of ZTP riboswitch ligands

The challenge of targeting RNA with small molecules necessitates a better understanding of RNA-ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand-induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure-informed design, crystallography, and machine learning-augmented all-atom molecular dynamics simulations (MD) we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the ZTP riboswitch, a ligand-binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on-rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.

Embracing exascale computing in nucleic acid simulations

This mini-review reports the recent advances in biomolecular simulations, particularly for nucleic acids, and provides the potential effects of the emerging exascale computing on nucleic acid simulations, emphasizing the need for advanced computational strategies to fully exploit this technological frontier. Specifically, we introduce recent breakthroughs in computer architectures for large-scale biomolecular simulations and review the simulation protocols for nucleic acids regarding force fields, enhanced sampling methods, coarse-grained models, and interactions with ligands. We also explore the integration of machine learning methods into simulations, which promises to significantly enhance the predictive modeling of biomolecules and the analysis of complex data generated by the exascale simulations. Finally, we discuss the challenges and perspectives for biomolecular simulations as we enter the dawning exascale computing era.

Mechanism of phosphate release from actin filaments

After ATP-actin monomers assemble filaments, the ATP's [Formula: see text]-phosphate is hydrolyzedwithin seconds and dissociates over minutes. We used all-atom molecular dynamics simulations to sample the release of phosphate from filaments and study residues that gate release. Dissociation of phosphate from Mg2+ is rate limiting and associated with an energy barrier of 20 kcal/mol, consistent with experimental rates of phosphate release. Phosphate then diffuses within an internal cavity toward a gate formed by R177, as suggested in prior computational studies and cryo-EM structures. The gate is closed when R177 hydrogen bonds with N111 and is open when R177 forms a salt bridge with D179. Most of the time, interactions of R177 with other residues occlude the phosphate release pathway. Machine learning analysis reveals that the occluding interactions fluctuate rapidly, underscoring the secondary role of backdoor gate opening in Pi release, in contrast with the previous hypothesis that gate opening is the primary event.

Characterization of the dual regulation by a c-di-GMP riboswitch Bc1 with a long expression platform from Bacillus thuringiensis

A riboswitch generally regulates the expression of its downstream genes through conformational change in its expression platform (EP) upon ligand binding. The cyclic diguanosine monophosphate (c-di-GMP) class I riboswitch Bc1 is widespread and conserved among Bacillus cereus group species. In this study, we revealed that Bc1 has a long EP with two typical ρ-independent terminator sequences 28 bp apart. The upstream terminator T1 is dominant in vitro, while downstream terminator T2 is more efficient in vivo. Through mutation analysis, we elucidated that Bc1 exerts a rare and incoherent "transcription-translation" dual regulation with T2 playing a crucial role. However, we found that Bc1 did not respond to c-di-GMP under in vitro transcription conditions, and the expressions of downstream genes did not change with fluctuation in intracellular c-di-GMP concentration. To explore this puzzle, we conducted SHAPE-MaP and confirmed the interaction of Bc1 with c-di-GMP. This shows that as c-di-GMP concentration increases, T1 unfolds but T2 remains almost intact and functional. The presence of T2 masks the effect of T1 unwinding, resulting in no response of Bc1 to c-di-GMP. The high Shannon entropy values of EP region imply the potential alternative structures of Bc1. We also found that zinc uptake regulator can specifically bind to the dual terminator coding sequence and slightly trigger the response of Bc1 to c-di-GMP. This work will shed light on the dual-regulation riboswitch and enrich our understanding of the RNA world.IMPORTANCEIn nature, riboswitches are involved in a variety of metabolic regulation, most of which preferentially regulate transcription termination or translation initiation of downstream genes in specific ways. Alternatively, the same or different riboswitches can exist in tandem to enhance regulatory effects or respond to multiple ligands. However, many putative conserved riboswitches have not yet been experimentally validated. Here, we found that the c-di-GMP riboswitch Bc1 with a long EP could form a dual terminator and exhibit non-canonical and incoherent "transcription-translation" dual regulation. Besides, zinc uptake regulator specifically bound to the coding sequence of the Bc1 EP and slightly mediated the action of Bc1. The application of SHAPE-MaP to the dual regulation mechanism of Bc1 may establish the foundation for future studies of such complex untranslated regions in other bacterial genomes.

Mechanism of Phosphate Release from Actin Filaments

After ATP-actin monomers assemble filaments, the ATP's γ-phosphate is hydrolyzed within seconds and dissociates over minutes. We used all-atom molecular dynamics simulations to sample the release of phosphate from filaments and study residues that gate release. Dissociation of phosphate from Mg2+ is rate limiting and associated with an energy barrier of 20 kcal/mol, consistent with experimental rates of phosphate release. Phosphate then diffuses in an internal cavity toward a gate formed by R177 suggested in prior computational studies and cryo-EM structures. The gate is closed when R177 hydrogen bonds with N111 and is open when R177 forms a salt bridge with D179. Most of the time interactions of R177 with other residues occludes the phosphate release pathway. Machine learning analysis reveals that the occluding interactions fluctuate rapidly, underscoring the secondary role of backdoor gate opening in Pi release, in contrast with the previous hypothesis that gate opening is the primary event.

Mechanistic Analysis of Riboswitch Ligand Interactions Provides Insights into Pharmacological Control over Gene Expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Mechanism of Ligand Binding to Theophylline RNA Aptamer

Studying RNA-ligand interactions and quantifying their binding thermodynamics and kinetics are of particular relevance in the field of drug discovery. Here, we combined biochemical binding assays and accelerated molecular simulations to investigate ligand binding and dissociation in RNA using the theophylline-binding RNA as a model system. All-atom simulations using a Ligand Gaussian accelerated Molecular Dynamics method (LiGaMD) have captured repetitive binding and dissociation of theophylline and caffeine to RNA. Theophylline's binding free energy and kinetic rate constants align with our experimental data, while caffeine's binding affinity is over 10,000 times weaker, and its kinetics could not be determined. LiGaMD simulations allowed us to identify distinct low-energy conformations and multiple ligand binding pathways to RNA. Simulations revealed a "conformational selection" mechanism for ligand binding to the flexible RNA aptamer, which provides important mechanistic insights into ligand binding to the theophylline-binding model. Our findings suggest that compound docking using a structural ensemble of representative RNA conformations would be necessary for structure-based drug design of flexible RNA.

Small molecule approaches to targeting RNA

The development of innovative methodologies to identify RNA binders has attracted enormous attention in chemical biology and drug discovery. Although antibiotics targeting bacterial ribosomal RNA have been on the market for decades, the renewed interest in RNA targeting reflects the need to better understand complex intracellular processes involving RNA. In this context, small molecules are privileged tools used to explore the biological functions of RNA and to validate RNAs as therapeutic targets, and they eventually are to become new drugs. Despite recent progress, the rational design of specific RNA binders requires a better understanding of the interactions which occur with the RNA target to reach the desired biological response. In this Review, we discuss the challenges to approaching this underexplored chemical space, together with recent strategies to bind, interact and affect biologically relevant RNAs.

LncRNAs: the good, the bad, and the unknown

Long non-coding RNAs (lncRNAs) are significant contributors in maintaining genomic integrity through epigenetic regulation. LncRNAs can interact with chromatin-modifying complexes in both cis and trans pathways, drawing them to specific genomic loci and influencing gene expression via DNA methylation, histone modifications, and chromatin remodeling. They can also operate as building blocks to assemble different chromatin-modifying components, facilitating their interactions and gene regulatory functions. Deregulation of these molecules has been associated with various human diseases, including cancer, cardiovascular disease, and neurological disorders. Thus, lncRNAs are implicated as potential diagnostic indicators and therapeutic targets. This review discusses the current understanding of how lncRNAs mediate epigenetic control, genomic integrity, and their putative functions in disease pathogenesis.

RNA-Selective Small-Molecule Ligands: Recent Advances in Live-Cell Imaging and Drug Discovery

RNA structures, including those formed from coding and noncoding RNAs, alternative to protein-based drug targets, could be a promising target of small molecules for drug discovery against various human diseases, particularly in anticancer, antibacterial and antivirus development. The normal cellular activity of cells is critically dependent on the function of various RNA molecules generated from DNA transcription. Moreover, many studies support that mRNA-targeting small molecules may regulate the synthesis of disease-related proteins via the non-covalent mRNA-ligand interactions that do not involve gene modification. RNA-ligand interaction is thus an attractive approach to address the challenge of "undruggable" proteins in drug discovery because the intracellular activity of these proteins is hard to be suppressed with small molecule ligands. We selectively surveyed a specific area of RNA structure-selective small molecule ligands in fluorescence live cell imaging and drug discovery because the area was currently underexplored. This state-of-the-art review thus mainly focuses on the research published within the past three years and aims to provide the most recent information on this research area; hopefully, it could be complementary to the previously reported reviews and give new insights into the future development on RNA-specific small molecule ligands for live cell imaging and drug discovery.

Targeting Riboswitches with Beta-Axial-Substituted Cobalamins

RNA-targeting small-molecule therapeutics is an emerging field hindered by an incomplete understanding of the basic principles governing RNA-ligand interactions. One way to advance our knowledge in this area is to study model systems where these interactions are better understood, such as riboswitches. Riboswitches bind a wide array of small molecules with high affinity and selectivity, providing a wealth of information on how RNA recognizes ligands through diverse structures. The cobalamin-sensing riboswitch is a particularly useful model system, as similar sequences show highly specialized binding preferences for different biological forms of cobalamin. This riboswitch is also widely dispersed across bacteria and therefore holds strong potential as an antibiotic target. Many synthetic cobalamin forms have been developed for various purposes including therapeutics, but their interaction with cobalamin riboswitches is yet to be explored. In this study, we characterize the interactions of 11 cobalamin derivatives with three representative cobalamin riboswitches using in vitro binding experiments (both chemical footprinting and a fluorescence-based assay) and a cell-based reporter assay. The derivatives show productive interactions with two of the three riboswitches, demonstrating simultaneous plasticity and selectivity within these RNAs. The observed plasticity is partially achieved through a novel structural rearrangement within the ligand binding pocket, providing insight into how similar RNA structures can be targeted. As the derivatives also show in vivo functionality, they serve as several potential lead compounds for further drug development.

Metadynamics simulations for the investigation of drug loading on functionalized inorganic nanoparticles

Inorganic nanoparticles show promising properties that allow them to be efficiently used as drug carriers. The main limitation in this type of application is currently the drug loading capacity, which can be overcome with a proper functionalization of the nanoparticle surface. In this study, we present, for the first time, a computational approach based on metadynamics to estimate the binding free energy of the doxorubicin drug (DOX) to a functionalized TiO₂ nanoparticle under different pH conditions. On a thermodynamic basis, we demonstrate the robustness of our approach to capture the overall mechanism behind the pH-triggered release of DOX due to environmental pH changes. Notably, binding free energy estimations align well with what is expected for a pH-sensitive drug delivery system. Based on our results, we envision the use of metadynamics as a promising computational tool for the rational design and in silico optimization of organic ligands with improved drug carrier properties.

Internal Normal Mode Analysis Applied to RNA Flexibility and Conformational Changes

We investigated the capability of internal normal modes to reproduce RNA flexibility and predict observed RNA conformational changes and, notably, those induced by the formation of RNA-protein and RNA-ligand complexes. Here, we extended our iNMA approach developed for proteins to study RNA molecules using a simplified representation of the RNA structure and its potential energy. Three data sets were also created to investigate different aspects. Despite all the approximations, our study shows that iNMA is a suitable method to take into account RNA flexibility and describe its conformational changes opening the route to its applicability in any integrative approach where these properties are crucial.

A riboswitch separated from its ribosome-binding site still regulates translation

Riboswitches regulate downstream gene expression by binding cellular metabolites. Regulation of translation initiation by riboswitches is posited to occur by metabolite-mediated sequestration of the Shine-Dalgarno sequence (SDS), causing bypass by the ribosome. Recently, we solved a co-crystal structure of a prequeuosine1-sensing riboswitch from Carnobacterium antarcticum that binds two metabolites in a single pocket. The structure revealed that the second nucleotide within the gene-regulatory SDS, G34, engages in a crystal contact, obscuring the molecular basis of gene regulation. Here, we report a co-crystal structure wherein C10 pairs with G34. However, molecular dynamics simulations reveal quick dissolution of the pair, which fails to reform. Functional and chemical probing assays inside live bacterial cells corroborate the dispensability of the C10-G34 pair in gene regulation, leading to the hypothesis that the compact pseudoknot fold is sufficient for translation attenuation. Remarkably, the C. antarcticum aptamer retained significant gene-regulatory activity when uncoupled from the SDS using unstructured spacers up to 10 nucleotides away from the riboswitch-akin to steric-blocking employed by sRNAs. Accordingly, our work reveals that the RNA fold regulates translation without SDS sequestration, expanding known riboswitch-mediated gene-regulatory mechanisms. The results infer that riboswitches exist wherein the SDS is not embedded inside a stable fold.