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.

Structure-Informed Design of an Ultra Bright RNA-activated Fluorophore.

Fluorogenic RNAs such as the Mango aptamers are uniquely powerful tools for imaging RNA. A central challenge has been to develop brighter, more specific, and higher affinity aptamer-ligand systems for cellular imaging. Here, we report an ultra-bright fluorophore for the Mango II system discovered using a structure-informed, fragment-based small molecule microarray approach. The new dye, Structure informed, Array-enabled LigAnD 1 (SALAD1) exhibits 3.5-fold brighter fluorescence than TO1-Biotin and subnanomolar aptamer affinity. Improved performance comes solely from alteration of dye-RNA interactions, without alteration of the chromophore itself. Multiple high-resolution structures reveal a unique and specific binding mode for the new dye resulting from improved pocket occupancy, a more defined binding pose, and a novel bonding interaction with potassium. The dye notably improves in-cell confocal RNA imaging. This work provides both introduces a new RNA-activated fluorophore and also a powerful demonstration of how to leverage fragment-based ligand discovery against RNA targets.

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.

Competitive Microarray Screening Reveals Functional Ligands for the DHX15 RNA G-Quadruplex.

RNAs are increasingly considered valuable therapeutic targets, and the development of methods to identify and validate both RNA targets and ligands is more important than ever. Here, we utilized a bioinformatic approach to identify a hairpin-containing RNA G-quadruplex (rG4) in the 5' untranslated region (5' UTR) of DHX15 mRNA. By using a novel competitive small molecule microarray (SMM) approach, we identified a compound that specifically binds to the DHX15 rG4 (K D = 12.6 ± 1.0 µM). This rG4 directly impacts translation of a DHX15 reporter mRNA in vitro, and binding of our compound (F1) to the structure inhibits translation up to 57% (IC50 = 22.9 ± 3.8 µM). This methodology allowed us to identify and target the mRNA of a cancer-relevant helicase with no known inhibitors. Our target identification method and the novelty of our screening approach make our work informative for future development of novel small molecule cancer therapeutics for RNA targets.

KDM3B inhibitors disrupt the oncogenic activity of PAX3-FOXO1 in fusion-positive rhabdomyosarcoma.

Fusion-positive rhabdomyosarcoma (FP-RMS) is an aggressive pediatric sarcoma driven primarily by the PAX3-FOXO1 fusion oncogene, for which therapies targeting PAX3-FOXO1 are lacking. Here, we screen 62,643 compounds using an engineered cell line that monitors PAX3-FOXO1 transcriptional activity identifying a hitherto uncharacterized compound, P3FI-63. RNA-seq, ATAC-seq, and docking analyses implicate histone lysine demethylases (KDMs) as its targets. Enzymatic assays confirm the inhibition of multiple KDMs with the highest selectivity for KDM3B. Structural similarity search of P3FI-63 identifies P3FI-90 with improved solubility and potency. Biophysical binding of P3FI-90 to KDM3B is demonstrated using NMR and SPR. P3FI-90 suppresses the growth of FP-RMS in vitro and in vivo through downregulating PAX3-FOXO1 activity, and combined knockdown of KDM3B and KDM1A phenocopies P3FI-90 effects. Thus, we report KDM inhibitors P3FI-63 and P3FI-90 with the highest specificity for KDM3B. Their potent suppression of PAX3-FOXO1 activity indicates a possible therapeutic approach for FP-RMS and other transcriptionally addicted cancers.

Decoding complexity in biomolecular recognition of DNA i-motifs with microarrays.

DNA i-motifs (iMs) are non-canonical C-rich secondary structures implicated in numerous cellular processes. Though iMs exist throughout the genome, our understanding of iM recognition by proteins or small molecules is limited to a few examples. We designed a DNA microarray containing 10976 genomic iM sequences to examine the binding profiles of four iM-binding proteins, mitoxantrone and the iMab antibody. iMab microarray screens demonstrated that pH 6.5, 5% BSA buffer was optimal, and fluorescence was correlated with iM C-tract length. hnRNP K broadly recognizes diverse iM sequences, favoring 3-5 cytosine repeats flanked by thymine-rich loops of 1-3 nucleotides. Array binding mirrored public ChIP-Seq datasets, in which 35% of well-bound array iMs are enriched in hnRNP K peaks. In contrast, other reported iM-binding proteins had weaker binding or preferred G-quadruplex (G4) sequences instead. Mitoxantrone broadly binds both shorter iMs and G4s, consistent with an intercalation mechanism. These results suggest that hnRNP K may play a role in iM-mediated regulation of gene expression in vivo, whereas hnRNP A1 and ASF/SF2 are possibly more selective in their binding preferences. This powerful approach represents the most comprehensive investigation of how biomolecules selectively recognize genomic iMs to date.

Discovery of Small Molecules Targeting the Frameshifting Element RNA in SARS-CoV-2 Viral Genome.

Targeting structured RNA elements in the SARS-CoV-2 viral genome with small molecules is an attractive strategy for pharmacological control over viral replication. In this work, we report the discovery of small molecules that target the frameshifting element (FSE) in the SARS-CoV-2 RNA genome using high-throughput small-molecule microarray (SMM) screening. A new class of aminoquinazoline ligands for the SARS-CoV-2 FSE are synthesized and characterized using multiple orthogonal biophysical assays and structure-activity relationship (SAR) studies. This work reveals compounds with mid-micromolar binding affinity (KD = 60 ± 6 µM) to the FSE RNA and supports a binding mode distinct from previously reported FSE binders MTDB and merafloxacin. In addition, compounds are active in in vitro dual-luciferase and in-cell dual-fluorescent-reporter frameshifting assays, highlighting the promise of targeting structured elements of RNAs with druglike compounds to alter expression of viral proteins.

Proximity-Induced Nucleic Acid Degrader (PINAD) Approach to Targeted RNA Degradation Using Small Molecules.

Nature has evolved intricate machinery to target and degrade RNA, and some of these molecular mechanisms can be adapted for therapeutic use. Small interfering RNAs and RNase H-inducing oligonucleotides have yielded therapeutic agents against diseases that cannot be tackled using protein-centered approaches. Because these therapeutic agents are nucleic acid-based, they have several inherent drawbacks which include poor cellular uptake and stability. Here we report a new approach to target and degrade RNA using small molecules, proximity-induced nucleic acid degrader (PINAD). We have utilized this strategy to design two families of RNA degraders which target two different RNA structures within the genome of SARS-CoV-2: G-quadruplexes and the betacoronaviral pseudoknot. We demonstrate that these novel molecules degrade their targets using in vitro, in cellulo, and in vivo SARS-CoV-2 infection models. Our strategy allows any RNA binding small molecule to be converted into a degrader, empowering RNA binders that are not potent enough to exert a phenotypic effect on their own. PINAD raises the possibility of targeting and destroying any disease-related RNA species, which can greatly expand the space of druggable targets and diseases.

DECODING COMPLEXITY IN BIOMOLECULAR RECOGNITION OF DNA I-MOTIFS.

DNA i-motifs (iMs) are non-canonical C-rich secondary structures implicated in numerous cellular processes. Though iMs exist throughout the genome, our understanding of iM recognition by proteins or small molecules is limited to a few examples. We designed a DNA microarray containing 10,976 genomic iM sequences to examine the binding profiles of four iM-binding proteins, mitoxantrone, and the iMab antibody. iMab microarray screens demonstrated that pH 6.5, 5% BSA buffer was optimal, and fluorescence was correlated with iM C-tract length. hnRNP K broadly recognizes diverse iM sequences, favoring 3-5 cytosine repeats flanked by thymine-rich loops of 1-3 nucleotides. Array binding mirrored public ChIP-Seq datasets, in which 35% of well-bound array iMs are enriched in hnRNP K peaks. In contrast, other reported iM-binding proteins had weaker binding or preferred G-quadruplex (G4) sequences instead. Mitoxantrone broadly binds both shorter iMs and G4s, consistent with an intercalation mechanism. These results suggest that hnRNP K may play a role in iM-mediated regulation of gene expression in vivo, whereas hnRNP A1 and ASF/SF2 are possibly more selective in their binding preferences. This powerful approach represents the most comprehensive investigation of how biomolecules selectively recognize genomic iMs to date.

Investigating the NRAS 5' UTR as a target for small molecules.

Neuroblastoma RAS (NRAS) is an oncogene that is deregulated and highly mutated in cancers including melanomas and acute myeloid leukemias. The 5' untranslated region (UTR) (5' UTR) of the NRAS mRNA contains a G-quadruplex (G4) that regulates translation. Here we report a novel class of small molecule that binds to the G4 structure located in the 5' UTR of the NRAS mRNA. We used a small molecule microarray screen to identify molecules that selectively bind to the NRAS-G4 with submicromolar affinity. One compound inhibits the translation of NRAS in vitro but showed only moderate effects on the NRAS levels in cellulo. Rapid Amplification of cDNA Ends and RT-PCR analysis revealed that the predominant NRAS transcript does not possess the G4 structure. Thus, although NRAS transcripts lack a G4 in many cell lines the concept of targeting folded regions within 5' UTRs to control translation remains a highly attractive strategy.

Crystal Structure of an i-Motif from the HRAS Oncogene Promoter.

An i-motif is a non-canonical DNA structure implicated in gene regulation and linked to cancers. The C-rich strand of the HRAS oncogene, 5'-CGCCCGTGCCCTGCGCCCGCAACCCGA-3' (herein referred to as iHRAS), forms an i-motif in vitro but its exact structure was unknown. HRAS is a member of the RAS proto-oncogene family. About 19 % of US cancer patients carry mutations in RAS genes. We solved the structure of iHRAS at 1.77 Šresolution. The structure reveals that iHRAS folds into a double hairpin. The two double hairpins associate in an antiparallel fashion, forming an i-motif dimer capped by two loops on each end and linked by a connecting region. Six C-C+ base pairs form each i-motif core, and the core regions are extended by a G-G base pair and a cytosine stacking. Extensive canonical and non-canonical base pairing and stacking stabilizes the connecting region and loops. The iHRAS structure is the first atomic resolution structure of an i-motif from a human oncogene. This structure sheds light on i-motifs folding and function in the cell.

Machine Learning Informs RNA-Binding Chemical Space.

Small molecule targeting of RNA has emerged as a new frontier in medicinal chemistry, but compared to the protein targeting literature our understanding of chemical matter that binds to RNA is limited. In this study, we reported Repository Of BInders to Nucleic acids (ROBIN), a new library of nucleic acid binders identified by small molecule microarray (SMM) screening. The complete results of 36 individual nucleic acid SMM screens against a library of 24 572 small molecules were reported (including a total of 1 627 072 interactions assayed). A set of 2 003 RNA-binding small molecules was identified, representing the largest fully public, experimentally derived library of its kind to date. Machine learning was used to develop highly predictive and interpretable models to characterize RNA-binding molecules. This work demonstrates that machine learning algorithms applied to experimentally derived sets of RNA binders are a powerful method to inform RNA-targeted chemical space.

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

While there is increasing interest in the study of RNA as a therapeutic target, efforts to understand RNA-ligand recognition at the molecular level lag far behind our understanding of protein-ligand recognition. This problem is complicated due to the more than 10 orders of magnitude in time scales involved in RNA dynamics and ligand binding events, making it not straightforward to design experiments or simulations. Here, we make use of artificial intelligence (AI)-augmented molecular dynamics simulations to directly observe ligand dissociation for cognate and synthetic ligands from a riboswitch system. The site-specific flexibility profiles from our simulations are compared with in vitro measurements of flexibility using selective 2' hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP). Our simulations reproduce known relative binding affinity profiles for the cognate and synthetic ligands, and pinpoint how both ligands make use of different aspects of riboswitch flexibility. On the basis of our dissociation trajectories, we also make and validate predictions of pairs of mutations for both the ligand systems that would show differing binding affinities. These mutations are distal to the binding site and could not have been predicted solely on the basis of structure. The methodology demonstrated here shows how molecular dynamics simulations with all-atom force-fields have now come of age in making predictions that complement existing experimental techniques and illuminate aspects of systems otherwise not trivial to understand.

Dataset of high-throughput ligand screening against the RNA Packaging Signals regulating Hepatitis B Virus nucleocapsid formation.

Multiple ssRNA viruses which infect bacteria, plants or humans use RNA Packaging Signal (PS)-mediated regulation during assembly to package their genomes faithfully and efficiently. PSs typically comprise short nucleotide recognition motifs, most often presented in the unpaired region of RNA stem-loops, and often bind their cognate coat proteins (CPs) with nanomolar affinity. PSs identified to date are resilient in the face of the typical error prone replication of their virus-coded polymerases, making them potential drug targets. An immobilised array of small molecular weight, drug-like compounds was panned against a fluorescently-labelled oligonucleotide encompassing the most conserved Hepatitis B Virus (HBV) PS, PS1, known to be a major determinant in nucleocapsid formation. This identified > 70 compounds that bind PS1 uniquely in the array. The commercially available 66 of these were tested for their potential effect(s) on HBV nucleocapsid-like particle (NCP) assembly in vitro, which identified potent assembly inhibitors. Here, we describe a high-throughput screen for such effects using employing fluorescence anisotropy in a 96-well microplate format. HBV genomic RNAs (gRNA) and short oligonucleotides encompassing PS1 were 5' labelled with an Alexa Fluor 488 dye. Excess (with respect to stoichiometric T = 4 NCP formation) HBV core protein (Cp) dimers were titrated robotically into solutions containing each of these RNAs stepwise, using a Biomek 4000 liquid handling robot. The anisotropy values of these mixtures were monitored using a POLARstar microplate reader. NCP-like structures were challenged with RNase A to identify reactions that did not result in complete NCP formation. The results imply that ∼50% of the compounds prevent complete NCP formation, highlighting both PS-meditated assembly and the PS-binding compounds as potential directly-acting anti-virals with a novel molecular target. Importantly, this method allows high-throughput in vitro screening for assembly inhibitors in this major human pathogen.

Dysregulation of Hepatitis B Virus Nucleocapsid Assembly in vitro by RNA-binding Small Ligands.

RNA sequences/motifs dispersed across the genome of Hepatitis B Virus regulate formation of nucleocapsid-like particles (NCPs) by core protein (Cp) in vitro, in an epsilon/polymerase-independent fashion. These multiple RNA Packaging Signals (PSs) can each form stem-loops encompassing a Cp-recognition motif, -RGAG-, in their loops. Drug-like molecules that bind the most important of these PS sites for NCP assembly regulation with nanomolar affinities, were identified by screening an immobilized ligand library with a fluorescently-labelled, RNA oligonucleotide encompassing this sequence. Sixty-six of these "hits", with affinities ranging from low nanomolar to high micromolar, were purchased as non-immobilized versions. Their affinities for PSs and effects on NCP assembly were determined in vitro by Surface Plasmon Resonance. High-affinity ligand binding is dependent on the presence of an -RGAG- motif within the loop of the PS, consistent with ligand cross-binding between PS sites. Simple structure-activity relationships show that it is also dependent on the presence of specific functional groups in these ligands. Some compounds are potent inhibitors of in vitro NCP assembly at nanomolar concentrations. Despite appropriate logP values, these ligands do not inhibit HBV replication in cell culture. However, modelling confirms the potential of using PS-binding ligands to target NCP assembly as a novel anti-viral strategy. This also allows for computational exploration of potential synergic effects between anti-viral ligands directed at distinct molecular targets in vivo. HBV PS-regulated assembly can be dysregulated by novel small molecule RNA-binding ligands opening a novel target for developing directly-acting anti-virals against this major pathogen.

Heterobifunctional molecules tackle targeted protein dephosphorylation.

Recently, Chen and colleagues reported the development of phosphatase-targeting chimeric molecules (PhosTACs), heterobifunctional small molecules that promote targeted, proximity-induced protein dephosphorylation. This strategy represents an innovative approach to selectively manipulate phosphoprotein function and provides proof-of-concept for a new class of bifunctional small molecules in the chemical biologist's toolbox.

Structural insights of the conserved "priming loop" of hepatitis B virus pre-genomic RNA.

Initiation of protein-primed (-) strand DNA synthesis in hepatitis B virus (HBV) requires interaction of the viral polymerase with a cis-acting regulatory signal, designated epsilon (ε), located at the 5'-end of its pre-genomic RNA (pgRNA). Binding of polymerase to ε is also necessary for pgRNA encapsidation. While the mechanistic basis of this interaction remains elusive, mutagenesis studies suggest its internal 6-nt "priming loop" provides an important structural contribution. ε might therefore be considered a promising target for small molecule interventions to complement current nucleoside-analog based anti-HBV therapies. An ideal prerequisite to any RNA-directed small molecule strategy would be a detailed structural description of this important element. Herein, we present a solution NMR structure for HBV ε which, in combination with molecular dynamics and docking simulations, reports on a flexible ligand "pocket", reminiscent of those observed in proteins. We also demonstrate the binding of the selective estrogen receptor modulators (SERMs) Raloxifene, Bazedoxifene, and a de novo derivative to the priming loop.Communicated by Ramaswamy H. Sarma.

Dynamic bulge nucleotides in the KSHV PAN ENE triple helix provide a unique binding platform for small molecule ligands.

Cellular and virus-coded long non-coding (lnc) RNAs support multiple roles related to biological and pathological processes. Several lncRNAs sequester their 3' termini to evade cellular degradation machinery, thereby supporting disease progression. An intramolecular triplex involving the lncRNA 3' terminus, the element for nuclear expression (ENE), stabilizes RNA transcripts and promotes persistent function. Therefore, such ENE triplexes, as presented here in Kaposi's sarcoma-associated herpesvirus (KSHV) polyadenylated nuclear (PAN) lncRNA, represent targets for therapeutic development. Towards identifying novel ligands targeting the PAN ENE triplex, we screened a library of immobilized small molecules and identified several triplex-binding chemotypes, the tightest of which exhibits micromolar binding affinity. Combined biophysical, biochemical, and computational strategies localized ligand binding to a platform created near a dinucleotide bulge at the base of the triplex. Crystal structures of apo (3.3 Å) and ligand-soaked (2.5 Å) ENE triplexes, which include a stabilizing basal duplex, indicate significant local structural rearrangements within this dinucleotide bulge. MD simulations and a modified nucleoside analog interference technique corroborate the role of the bulge and the base of the triplex in ligand binding. Together with recently discovered small molecules that reduce nuclear MALAT1 lncRNA levels by engaging its ENE triplex, our data supports the potential of targeting RNA triplexes with small molecules.

A chemical probe based on the PreQ1 metabolite enables transcriptome-wide mapping of binding sites.

The role of metabolite-responsive riboswitches in regulating gene expression in bacteria is well known and makes them useful systems for the study of RNA-small molecule interactions. Here, we study the PreQ1 riboswitch system, assessing sixteen diverse PreQ1-derived probes for their ability to selectively modify the class-I PreQ1 riboswitch aptamer covalently. For the most active probe (11), a diazirine-based photocrosslinking analog of PreQ1, X-ray crystallography and gel-based competition assays demonstrated the mode of binding of the ligand to the aptamer, and functional assays demonstrated that the probe retains activity against the full riboswitch. Transcriptome-wide mapping using Chem-CLIP revealed a highly selective interaction between the bacterial aptamer and the probe. In addition, a small number of RNA targets in endogenous human transcripts were found to bind specifically to 11, providing evidence for candidate PreQ1 aptamers in human RNA. This work demonstrates a stark influence of linker chemistry and structure on the ability of molecules to crosslink RNA, reveals that the PreQ1 aptamer/ligand pair are broadly useful for chemical biology applications, and provides insights into how PreQ1, which is similar in structure to guanine, interacts with human RNAs.