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.

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.

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.

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.

Targeting a noncanonical, hairpin-containing G-quadruplex structure from the MYCN gene.

The MYCN gene encodes the transcription factor N-Myc, a driver of neuroblastoma (NB). Targeting G-quadruplexes (G4s) with small molecules is attractive strategy to control the expression of undruggable proteins such as N-Myc. However, selective binders to G4s are challenging to identify due to the structural similarity of many G4s. Here, we report the discovery of a small molecule ligand (4) that targets the noncanonical, hairpin containing G4 structure found in the MYCN gene using small molecule microarrays (SMMs). Unlike many G4 binders, the compound was found to bind to a pocket at the base of the hairpin region of the MYCN G4. This compound stabilizes the G4 and has affinity of 3.5 ± 1.6 µM. Moreover, an improved analog, MY-8, suppressed levels of both MYCN and MYCNOS (a lncRNA embedded within the MYCN gene) in NBEB neuroblastoma cells. This work indicates that the approach of targeting complex, hybrid G4 structures that exist throughout the human genome may be an applicable strategy to achieve selectivity for targeting disease-relevant genes including protein coding (MYCN) as well as non-coding (MYCNOS) gene products.

Small molecule microarray identifies inhibitors of tyrosyl-DNA phosphodiesterase 1 that simultaneously access the catalytic pocket and two substrate binding sites.

Tyrosyl-DNA phosphodiesterase 1 (TDP1) is a member of the phospholipase D family of enzymes, which catalyzes the removal of both 3'- and 5'-DNA phosphodiester adducts. Importantly, it is capable of reducing the anticancer effects of type I topoisomerase (TOP1) inhibitors by repairing the stalled covalent complexes of TOP1 with DNA. It achieves this by promoting the hydrolysis of the phosphodiester bond between the Y723 residue of human TOP1 and the 3'-phosphate of its DNA substrate. Blocking TDP1 function is an attractive means of enhancing the efficacy of TOP1 inhibitors and overcoming drug resistance. Previously, we reported the use of an X-ray crystallographic screen of more than 600 fragments to identify small molecule variations on phthalic acid and hydroxyquinoline motifs that bind within the TDP1 catalytic pocket. Yet, the majority of these compounds showed limited (millimolar) TDP1 inhibitory potencies. We now report examining a 21 000-member library of drug-like Small Molecules in Microarray (SMM) format for their ability to bind Alexa Fluor 647 (AF647)-labeled TDP1. The screen identified structurally similar N,2-diphenylimidazo[1,2-a]pyrazin-3-amines as TDP1 binders and catalytic inhibitors. We then explored the core heterocycle skeleton using one-pot Groebke-Blackburn-Bienayme multicomponent reactions and arrived at analogs having higher inhibitory potencies. Solving TDP1 co-crystal structures of a subset of compounds showed their binding at the TDP1 catalytic site, while mimicking substrate interactions. Although our original fragment screen differed significantly from the current microarray protocol, both methods identified ligand-protein interactions containing highly similar elements. Importantly inhibitors identified through the SMM approach show competitive inhibition against TDP1 and access the catalytic phosphate-binding pocket, while simultaneously providing extensions into both the substrate DNA and peptide-binding channels. As such, they represent a platform for further elaboration of trivalent ligands, that could serve as a new genre of potent TDP1 inhibitors.

Three-Color Imaging Enables Simultaneous Screening of Multiple RNA Targets on Small Molecule Microarrays.

Small molecule microarray (SMM) technology has become a powerful tool used in high-throughput screening for target-based drug discovery. One area in which SMMs have found use is the identification of small molecule ligands for RNA. RNAs with unique secondary or tertiary three-dimensional structures are considered to be attractive targets for small molecules. Complex RNA structures can form hydrophobic pockets suitable for small molecule binding, representing an opportunity for developing novel therapeutics. Our lab has previously taken a target-based approach, screening a single target against many small molecules on an SMM platform. Here, we report a screening protocol for SMMs to investigate multiple RNAs simultaneously using multi-color imaging. By introducing a mixture containing different fluorophore-labeled RNAs, the fluorescence signal of each binding event can be observed simultaneously. Thus, the specificity of a hit compound binding to one RNA target over other highly abundant RNAs (such as tRNA or rRNA) can be easily evaluated. © 2020 Wiley Periodicals LLC. Basic Protocol: RNA screening on SMMs by multi-color imaging Support Protocol 1: Preparation of SMM slides Support Protocol 2: Fluorophore labeling of RNA through maleimide chemistry.

A Small Molecule Stabilizer of the MYC G4-Quadruplex Induces Endoplasmic Reticulum Stress, Senescence and Pyroptosis in Multiple Myeloma.

New approaches to target MYC include the stabilization of a guanine-rich, G-quadruplex (G4) tertiary DNA structure in the NHE III region of its promoter. Recent screening of a small molecule microarray platform identified a benzofuran, D089, that can stabilize the MYC G4 and inhibit its transcription. D089 induced both dose- and time-dependent multiple myeloma cell death mediated by endoplasmic reticulum induced stress. Unexpectedly, we uncovered two mechanisms of cell death: cellular senescence, as evidenced by increased levels of p16, p21 and γ-H2AX proteins and a caspase 3-independent mechanism consistent with pyroptosis. Cells treated with D089 exhibited high levels of the cleaved form of initiator caspase 8; but failed to show cleavage of executioner caspase 3, a classical apoptotic marker. Cotreatment with the a pan-caspase inhibitor Q-VD-OPh did not affect the cytotoxic effect of D089. In contrast, cleaved caspase 1, an inflammatory caspase downstream of caspases 8/9, was increased by D089 treatment. Cells treated with D089 in addition to either a caspase 1 inhibitor or siRNA-caspase 1 showed increased IC50 values, indicating a contribution of cleaved caspase 1 to cell death. Downstream effects of caspase 1 activation after drug treatment included increases in IL1B, gasdermin D cleavage, and HMGB1 translocation from the nucleus to the cytoplasm. Drug treated cells underwent a 'ballooning' morphology characteristic of pyroptosis, rather than 'blebbing' typically associated with apoptosis. ASC specks colocalized with NLRP3 in proximity ligation assays after drug treatment, indicating inflammasome activation and further confirming pyroptosis as a contributor to cell death. Thus, the small molecule MYC G4 stabilizer, D089, provides a new tool compound for studying pyroptosis. These studies suggest that inducing both tumor senescence and pyroptosis may have therapeutic potential for cancer treatment.