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.

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.

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.

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.

HNRNPH1-dependent splicing of a fusion oncogene reveals a targetable RNA G-quadruplex interaction.

The primary oncogenic event in ∼85% of Ewing sarcomas is a chromosomal translocation that generates a fusion oncogene encoding an aberrant transcription factor. The exact genomic breakpoints within the translocated genes, EWSR1 and FLI1, vary; however, in EWSR1, breakpoints typically occur within introns 7 or 8. We previously found that in Ewing sarcoma cells harboring EWSR1 intron 8 breakpoints, the RNA-binding protein HNRNPH1 facilitates a splicing event that excludes EWSR1 exon 8 from the EWS-FLI1 pre-mRNA to generate an in-frame mRNA. Here, we show that the processing of distinct EWS-FLI1 pre-mRNAs by HNRNPH1, but not other homologous family members, resembles alternative splicing of transcript variants of EWSR1 We demonstrate that HNRNPH1 recruitment is driven by guanine-rich sequences within EWSR1 exon 8 that have the potential to fold into RNA G-quadruplex structures. Critically, we demonstrate that an RNA mimetic of one of these G-quadruplexes modulates HNRNPH1 binding and induces a decrease in the growth of an EWSR1 exon 8 fusion-positive Ewing sarcoma cell line. Finally, we show that EWSR1 exon 8 fusion-positive cell lines are more sensitive to treatment with the pan-quadruplex binding molecule, pyridostatin (PDS), than EWSR1 exon 8 fusion-negative lines. Also, the treatment of EWSR1 exon 8 fusion-positive cells with PDS decreases EWS-FLI1 transcriptional activity, reversing the transcriptional deregulation driven by EWS-FLI1. Our findings illustrate that modulation of the alternative splicing of EWS-FLI1 pre-mRNA is a novel strategy for future therapeutics against the EWSR1 exon 8 containing fusion oncogenes present in a third of Ewing sarcoma.

Characterization of clinically used oral antiseptics as quadruplex-binding ligands.

Approaches to characterize the nucleic acid-binding properties of drugs and druglike small molecules are crucial to understanding the behavior of these compounds in cellular systems. Here, we use a Small Molecule Microarray (SMM) profiling approach to identify the preferential interaction between chlorhexidine, a widely used oral antiseptic, and the G-quadruplex (G4) structure in the KRAS oncogene promoter. The interaction of chlorhexidine and related drugs to the KRAS G4 is evaluated using multiple biophysical methods, including thermal melt, fluorescence titration and surface plasmon resonance (SPR) assays. Chlorhexidine has a specific low micromolar binding interaction with the G4, while related drugs have weaker and/or less specific interactions. Through NMR experiments and docking studies, we propose a plausible binding mode driven by both aromatic stacking and groove binding interactions. Additionally, cancer cell lines harbouring oncogenic mutations in the KRAS gene exhibit increased sensitivity to chlorhexidine. Treatment of breast cancer cells with chlorhexidine decreases KRAS protein levels, while a KRAS gene transiently expressed by a promoter lacking a G4 is not affected. This work confirms that known ligands bind broadly to G4 structures, while other drugs and druglike compounds can have more selective interactions that may be biologically relevant.

Custom G4 Microarrays Reveal Selective G-Quadruplex Recognition of Small Molecule BMVC: A Large-Scale Assessment of Ligand Binding Selectivity.

G-quadruplexes (G4) are considered new drug targets for human diseases such as cancer. More than 10,000 G4s have been discovered in human chromatin, posing challenges for assessing the selectivity of a G4-interactive ligand. 3,6-bis(1-Methyl-4-vinylpyridinium) carbazole diiodide (BMVC) is the first fluorescent small molecule for G4 detection in vivo. Our previous structural study shows that BMVC binds to the MYC promoter G4 (MycG4) with high specificity. Here, we utilize high-throughput, large-scale custom DNA G4 microarrays to analyze the G4-binding selectivity of BMVC. BMVC preferentially binds to the parallel MycG4 and selectively recognizes flanking sequences of parallel G4s, especially the 3'-flanking thymine. Importantly, the microarray results are confirmed by orthogonal NMR and fluorescence binding analyses. Our study demonstrates the potential of custom G4 microarrays as a platform to broadly and unbiasedly assess the binding selectivity of G4-interactive ligands, and to help understand the properties that govern molecular recognition.

Evidence for ligandable sites in structured RNA throughout the Protein Data Bank.

RNA has attracted considerable attention as a target for small molecules. However, methods to identify, study, and characterize suitable RNA targets have lagged behind strategies for protein targets. One approach that has received considerable attention for protein targets has been to utilize computational analysis to investigate ligandable "pockets" on proteins that are amenable to small molecule binding. These studies have shown that selected physical properties of pockets are important parameters that govern the ability of a structure to bind to small molecules. This work describes a similar analysis to study pockets on all RNAs in the Protein Data Bank (PDB). Using parameters such as buriedness, hydrophobicity, volume, and other properties, the set of all RNAs is analyzed and compared to all proteins. Considerable overlap is observed between the properties of pockets on RNAs and proteins. Thus, many RNAs are capable of populating conformations with pockets that are likely suitable for small molecule binding. Further, principal moment of inertia (PMI) calculations reveal that liganded RNAs exist in diverse structural space, much of which overlaps with protein structural space. Taken together, these results suggest that complex folded RNAs adopt unique structures with pockets that may represent viable opportunities for small molecule targeting.

Measuring RNA-Ligand Interactions with Microscale Thermophoresis.

In recent years, there has been dramatic growth in the study of RNA. RNA has gone from being known as an intermediate in the central dogma of molecular biology to a molecule with a large diversity of structure and function that is involved in all aspects of biology. As new functions are rapidly discovered, it has become clear that there is a need for RNA-targeting small molecule probes to investigate RNA biology and clarify the potential for therapeutics based on RNA-small molecule interactions. While a host of techniques exist to measure RNA-small molecule interactions, many of these have drawbacks that make them intractable for routine use and are often not broadly applicable. A newer technology called microscale thermophoresis (MST), which measures the directed migration of a molecule and/or molecule-ligand complex along a temperature gradient, can be used to measure binding affinities using very small amounts of sample. The high sensitivity of this technique enables measurement of affinity constants in the nanomolar and micromolar range. Here, we demonstrate how MST can be used to study a range of biologically relevant RNA interactions, including peptide-RNA interactions, RNA-small molecule interactions, and displacement of an RNA-bound peptide by a small molecule.

Macrophilones from the Marine Hydroid Macrorhynchia philippina Can Inhibit ERK Cascade Signaling.

Six new macrophilone-type pyrroloiminoquines were isolated and identified from an extract of the marine hydroid Macrorhynchia philippina. The proton-deficient and heteroatom-rich structures of macrophilones B-G (2-7) were elucidated by spectroscopic analysis and comparison of their data with those of the previously reported metabolite macrophilone A (1). Compounds 1-7 are the first pyrroloiminoquines to be reported from a hydroid. The macrophilones were shown to inhibit the enzymatic conjugation of SUMO to peptide substrates, and macrophilones A (1) and C (3) exhibit potent and selective cytotoxic properties in the NCI-60 anticancer screen. Bioinformatic analysis revealed a close association of the cytotoxicity profiles of 1 and 3 with two known B-Raf kinase inhibitory drugs. While compounds 1 and 3 showed no kinase inhibitory activity, they resulted in a dramatic decrease in cellular protein levels of selected components of the ERK signal cascade. As such, the chemical scaffold of the macrophilones could provide small-molecule therapeutic leads that target the ERK signal transduction pathway.

Microarray-based technologies for the discovery of selective, RNA-binding molecules.

The identification of small molecules that bind specifically to RNA is a challenge. However, the recent explosion in knowledge about the role RNA plays in a number of physiological processes apart from coding for protein sequences makes it a highly interesting target for chemical probes and therapeutics. One technology that has played an important role in the discovery of RNA-binding molecules is microarrays. Microarrays have been broadly employed to screen, profile, and quantify RNA interactions, and will likely play an important role in the discovery of new classes of ligands going forward. Here, we discuss the development of microarray technologies, including aminoglycoside, peptide, peptoid, and small molecule microarrays, and their use in studying RNA-interacting molecules.

Insights Into the Allosteric Inhibition of the SUMO E2 Enzyme Ubc9.

Conjugation of the small ubiquitin-like modifier (SUMO) to protein substrates is an important disease-associated posttranslational modification, although few inhibitors of this process are known. Herein, we report the discovery of an allosteric small-molecule binding site on Ubc9, the sole SUMO E2 enzyme. An X-ray crystallographic screen was used to identify two distinct small-molecule fragments that bind to Ubc9 at a site distal to its catalytic cysteine. These fragments and related compounds inhibit SUMO conjugation in biochemical assays with potencies of 1.9-5.8 mm. Mechanistic and biophysical analyses, coupled with molecular dynamics simulations, point toward ligand-induced rigidification of Ubc9 as a mechanism of inhibition.

Synthetic derivatives of the SUMO consensus sequence provide a basis for improved substrate recognition.

Protein sumoylation is a dynamic posttranslational modification that regulates a diverse subset of the proteome. The mechanism by which sumoylation enzymes recognize their cognate substrates is unclear, and the consequences of sumoylation remain difficult to predict. While small molecule probes of the sumoylation process could be valuable for understanding SUMO biology, few small molecules that modulate this process exist. Here, we report the synthesis and evaluation of over 600 oxime-containing peptide sumoylation substrates. Our work demonstrates that higher modification efficiency can be achieved with non-natural side chains that deviate substantially from the consensus site requirement of a hydrophobic substituent. Furthermore, docking studies suggest that these improved substrates mimic binding interactions that are used by other endogenous protein sequences through tertiary interactions. The development of these high efficiency substrates provides key mechanistic insights toward specific recognition of low molecular weight species in the sumoylation pathway.

Identification of biologically active, HIV TAR RNA-binding small molecules using small molecule microarrays.

Identifying small molecules that selectively bind to structured RNA motifs remains an important challenge in developing potent and specific therapeutics. Most strategies to find RNA-binding molecules have identified highly charged compounds or aminoglycosides that commonly have modest selectivity. Here we demonstrate a strategy to screen a large unbiased library of druglike small molecules in a microarray format against an RNA target. This approach has enabled the identification of a novel chemotype that selectively targets the HIV transactivation response (TAR) RNA hairpin in a manner not dependent on cationic charge. Thienopyridine 4 binds to and stabilizes the TAR hairpin with a Kd of 2.4 µM. Structure-activity relationships demonstrate that this compound achieves activity through hydrophobic and aromatic substituents on a heterocyclic core, rather than cationic groups typically required. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis was performed on a 365-nucleotide sequence derived from the 5' untranslated region (UTR) of the HIV-1 genome to determine global structural changes in the presence of the molecule. Importantly, the interaction of compound 4 can be mapped to the TAR hairpin without broadly disrupting any other structured elements of the 5' UTR. Cell-based anti-HIV assays indicated that 4 inhibits HIV-induced cytopathicity in T lymphocytes with an EC50 of 28 µM, while cytotoxicity was not observed at concentrations approaching 1 mM.

Synthesis of 2',3',4'-trihydroxyflavone (2-D08), an inhibitor of protein sumoylation.

Protein sumoylation is a dynamic posttranslational modification involved in diverse biological processes during cellular homeostasis and development. Recently sumoylation has been shown to play a critical role in cancer, although to date there are few small molecule probes available to inhibit enzymes involved in the SUMO conjugation process. As part of a program to identify and study inhibitors of sumoylation we recently reported the discovery that 2',3',4'-trihydroxyflavone (2-D08) is a cell permeable, mechanistically unique inhibitor of protein sumoylation. The work reported herein describes an efficient synthesis of 2-D08 as well as a structurally related but inactive isomer. We also report an unanticipated Wessely-Moser rearrangement that occurs under vigorous methyl ether deprotection conditions. This rearrangement likely gave rise to 2-D08 during a deprotection step, resulting in 2-D08 appearing as a contaminant in a screening well from a commercial supplier.

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.

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.

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.