Playing fast and loose with electrostatic anchors on RNA

Two weeks ago we discussed how to find ligand-binding sites in RNA. Last week we wrote about how difficult it is to find good ligands even for good binding sites in RNA. A recent open-access paper in J. Med. Chem. by Christian Kersten and colleagues at Johannes Gutenberg-University explores why targeting RNA is so tough.

 

The researchers were interested in two well-characterized riboswitches, naturally occurring RNA elements that bind to small molecules such as metabolites. Specifically, they chose to study a riboswitch that binds to S-adenosyl methionine (SAM, structure here) and a riboswitch that binds to prequeuosine-1 (PreQ₁) and prequeuosine-0 (PreQ₀). 

Due to the phosphate backbone, RNA is highly negatively charged. The researchers asked whether positively charged moieties on ligands can serve as “electrostatic anchors” to generally improve affinity, and if so whether this can lead to any design principles. Multiple biophysical techniques were used to study the interactions of the two riboswitches with various natural and synthetic ligands: surface plasmon resonance (SPR), isothermal titration calorimetry (ITC), and microscale thermophoresis (MST).

 

In the case of the SAM-VI riboswitch, the researchers compared the binding of SAM with closely related molecules having either one fewer positive charge (S-adenosyl homocysteine, or SAH) or synthetic ligands with the same or one more positive charge than SAM. Not surprisingly, SAM has the highest affinity, binding 20-50 fold more tightly than SAH. Further analysis suggested this is largely driven by an increased association rate, in which the positive charge accelerates the kinetics of binding. The driving energy for binding the ligands is enthalpic, but the favorable electrostatic interactions for more positively charged ligands are largely countered by an entropic penalty.

 

Similarly, the affinity of positively charged PreQ₁ for the PreQ₁ riboswitch is higher than the affinity of neutral PreQ₀, though not dramatically. As in the case of the SAM-VI riboswitch, the association rate of the positively charged ligand is more rapid than that of the neutral ligand. Binding for both ligands is highly enthalpic, with unfavorable entropy.

 

Previous reports had described other synthetic ligands for the PreQ₁ riboswitch, each with between one and three cationic centers. However these ligands showed no binding by ITC, questionable binding by MTC, and non-saturable, non-specific “loose binding” by SPR. Positive charges alone are not sufficient for high affinity, specific binding.

 

So what does it all mean? While adding positive charges can improve affinity of ligands for RNA, the increased affinity is usually not dramatic due to enthalpy-entropy compensation. The researchers note that, even for good ligands, the “thermodynamic binding profiles differ from typical protein-ligand interactions, where enthalpic and entropic contributions are usually more balanced.” Moreover, as we’ve noted, protein ligands often gain significant affinity with entropic gains by displacing "high energy water" molecules, but such opportunities are likely less common on the polar surface of RNA.

 

The affinity and ligand efficiency of PreQ₁ for its riboswitch are impressive, so clearly it is possible for small drug-like ligands to bind tightly to RNA. But this interaction is the product of countless eons of evolution. This careful paper suggests why building similarly effective synthetic ligands for most RNA will be difficult.

Targeting SARS-CoV-2 RNA – but not specifically

Last week we highlighted work suggesting that small molecule binding sites in RNA are most likely to be found in complex structures. A new open-access paper in Angew. Chem. Int. Ed. by Harald Schwalbe and collaborators at Goethe University Frankfurt and elsewhere provides both a case in point and an illustration of how difficult it is to target RNA.

 

The researchers had previously screened 15 RNAs from the SARS-CoV-2 virus, an effort we highlighted in 2021. In the new paper, the researchers focus on a portion of the frameshift element, which is important for directing viral replication from either of two partially overlapping open reading frames. The core of this RNA element is a roughly 69-nucleotide-long structure called a pseudoknot. Like most RNA sequences, this one can form multiple structures, including dimers, and the researchers used NMR, small-angle X-ray scattering (SAXS), and native gel electrophoresis to confirm that the construct was behaving as a homogenous monomer, consistent with three previously determined structures.

 

Based on some of the initial fragment hits, the researchers selected 50 similar molecules, of which only 14 were sufficiently soluble for screening. One of the more potent compounds, D05, initially showed promising activity in a ligand-detected NMR assay but turned out to be completely inactive when retested from a fresh stock. It turns out that D05 decomposes to compound 2, which was confirmed as active. Further modification led to compound 4, the most potent compound described. (Dissociation constants were determined by NMR, fluorescence, or both, and the two methods were in good agreement.)

Two-dimensional NMR with isotopically labeled RNA was used to try to determine the location of the binding site(s). Even with access to a 1.2 GHz magnet, the NMR peaks were severely overlapped, so the researchers used segmental isotopic labeling, in which just half of the RNA was labeled at a time. This exercise revealed potentially three different binding sites for compound 2.

 

The researchers also used two different computational approaches, Vina and RLdock, to predict binding sites, each of which could find one or two of the binding sites identified by NMR.

 

Several compounds were tested to see if they could block frameshifting in cell-lysates, and compound 2 showed 40% inhibition at 145 µM.

 

So far so good. But consistent with best practices, the researchers tested compounds 2 and 4 against phenylalanine tRNA. Unfortunately, the two ligands exhibited similar affinities to this control RNA as they did to the SARS-CoV-2 pseudoknot, despite the lack of sequence similarity. This suggests that these ligands bind to RNA nonspecifically. Perhaps this is not surprising given the three binding sites observed in a single 69-mer.

 

In the end, this is a thorough but sobering paper. Despite an impressive screening campaign with multiple biophysical methods, the best ligands seem to have modest affinity and low specificity. Drugging RNA still appears much more difficult than drugging proteins. But for either sort of target, this sort of careful work will be essential to find promising leads.

Identifying ligand-binding pockets in RNA, computationally and experimentally

Most drugs bind to proteins, but RNA provides many interesting targets. Unfortunately, finding drug-like small molecules that bind to RNA is difficult. A new paper in Proc. Nat. Acad. Sci. USA from Kevin Weeks and colleagues at University of North Carolina Chapel Hill provides tools to do so.

 

RNA presents several challenges for drug discovery. First, there are far fewer high-resolution structures than there are for proteins. This is in part due to the second challenge: RNA strands are often wriggly, able to form multiple conformations. And finally, RNA is highly charged and more polar than most proteins, so there are fewer opportunities for the hydrophobic interactions that often provide significant affinity in protein-ligand complexes.

 

These challenges have not deterred intrepid investigators: Practical Fragments first wrote about targeting RNA with fragments way back in 2009. However, examples of high-affinity ligands remain elusive, and in 2023 I wondered whether “most RNA is truly undruggable.”

 

The latest paper leaves me more optimistic. It describes a computational approach to find small-molecule binding sites in RNA. The researchers started with an open-source tool called fpocket, which was built for proteins. The fpocket program places virtual spheres all around a biomolecule, where each sphere contacts the center of four atoms. The size of each sphere depends on local curvature, and clusters of spheres define pockets.

 

To benchmark fpocket on RNA, the researchers first constructed a curated database of drug-like ligands bound to RNA. Of 538 RNA-ligand structures solved at the fairly low bar of ˂ 3.5 Å resolution, only 48 ligands were deemed drug-like by the quantitative estimate of drug-likeness (QED) score. (Although the QED score may be overly restrictive, and many approved drugs have low QED scores, setting a strict threshold means that any pockets identified are likely to be particularly attractive.)

 

Using default (protein-appropriate) parameters, fpocket identified just 63% of known ligand-binding sites in RNA, vs 83% for proteins. Worse, many predicted RNA pockets probably aren’t actually ligandable because they are too exposed to solvent. By tweaking parameters, the researchers improved performance of the program for RNA to 92%, and they also identified several attractive pockets that had previously been missed.

 

When the researchers applied the reparametrized program, redubbed fpocketR, to two bacterial ribosomes, they found several dozen pockets in each, including known antibiotic-binding sites. To assess whether the new pockets could bind fragments, they used an experimental approach called Frag-MaP, which uses fully functionalized fragment (FFF) probes containing a variable fragment, a photoreactive diazirine, and an alkyne. Treating bacterial cells with these FFF probes in the presence of UV light crosslinks them to nearby RNA. Crosslinked probes can then be isolated using click chemistry with the alkyne, and RNA sequencing reveals the sites of modification. Impressively, 89% of ligand binding sites found in the Frag-MaP experiments were predicted by fpocketR.

 

In another validation experiment, fpocketR identified pockets where 7 out of 17 antibiotics bind to bacterial ribosomes. Notably, all but one of the undetected pockets bind antibiotics such as aminoglycosides that don’t appear conventionally drug-like and indeed are not orally bioavailable.

 

Continuing to apply fpocketR to more RNAs led to the identification of dozens of new pockets. Interestingly, most of these pockets occur in complex RNA structures, such as multi-helix junctions or pseudoknots, rather than simpler structures such as bulges and consecutive loops. This could explain the paucity of fragment hits in a study we highlighted in 2023, which focused on simple loops.

 

Now that we know where to find attractive ligand-binding pockets in RNA, hopefully we will be more successful finding high-affinity ligands.

Small and simple, but novel and potent

Back in 2012 we wrote about GDB-17, a database of possible small molecules having up to 17 carbon, oxygen, nitrogen, sulfur, and halogen atoms, most of which have never been synthesized. Although novelty isn’t strictly necessary for fragments, as evidenced by the fact that 7-azaindole has given rise to three approved drugs, it’s certainly nice to have. In a new (open-access) J. Med. Chem. paper, Jürg Gertsch, Jean-Louis Reymond, and colleagues at the University of Bern synthesize fragments that had not been previously made and show that they are biologically active.

 

When you start drawing all possible small molecules you get lots of weird stuff, including an explosion of compounds containing multiple three- and four-membered rings, which may be difficult to make. The researchers wisely focused on “mono- and bicyclic ring systems containing only five-, six-, or seven-membered rings.” They further limited their search to molecules containing just carbon and one or two nitrogen atoms (as well as hydrogen, of course). Systematic enumeration led to 1139 scaffolds, ignoring stereochemistry, of which 680 had not been previously reported in PubChem. Out of these, three related scaffolds were chosen for investigation.

 

Computational retrosynthesis was used to devise routes to the three bicyclic scaffolds, and these were successfully synthesized, along with mono-benzylated versions, for a total of 14 molecules (including stereoisomers), all rule-of-three compliant. The online Polypharmacology Browser 2 (PPB2) was used to predict targets, and several monoamine transporters came up as potential hits. The molecules were tested against norepinephrine transporter (NET), dopamine transporter (DAT), serotonin transporter (SERT), and the σ-R1 receptor in radioligand displacement assays. None of the free diamines were active, but several of the benzylated compounds were, in particular compound 1a.

 

Compound 1a was initially made as a racemic mixture, and when the two enantiomers were resolved (R,R)-1a was found to be a mid-nanomolar inhibitor of NET while (S,S)-1a was 26-fold weaker. Compound (R,R)-1a was also a mid- to high nanomolar inhibitor of σ-R1, DAT, and SERT. Pharmacokinetic experiments in mice revealed that the molecule had poor oral bioavailability but remarkably high brain penetration and caused sedation. The researchers conducted additional mechanistic studies beyond the scope of this blog post and conclude that (R,R)-1a could be a lead for “neuropsychiatric disorders associated with monoamine dysregulation.”

 

There are several nice lessons in this paper. First, as we noted more than a decade ago, there is plenty of novelty at the bottom of chemical space. Moreover, and in contrast to our post last week, even small fragments can have high affinities. But novelty comes at a cost: synthesis of compound 2a required eight steps from an inexpensive starting material with an overall yield of just 9%, though this could certainly be optimized. Nonetheless, particularly for CNS-targeting drugs which usually need to be small in order to cross the blood brain barrier, the price might be worth paying.

 

Of course, even within this paper there are hundreds more scaffolds to look at than the three tested, and perhaps the researchers were lucky that their choices were biologically active. As computational methods continue to advance, it will be worthwhile turning them loose on GDB-17.