The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding.

The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.

Statistical mechanical prediction of ligand perturbation to RNA secondary structure and application to riboswitches.

The realization that noncoding RNA is implicated in numerous cellular processes, makes it imperative to understand and predict RNA-folding. RNA secondary structure prediction is more tractable than tertiary structure or protein structure. Yet insights into RNA structure-function relationships are complicated by coupling between RNA-folding and ligand-binding. Here, perturbations to equilibrium secondary structure conformational distributions for two riboswitches are calculated in the presence of bound cognate ligands. This work incorporates a key factor coupling ligand binding to RNA conformation but not considered in most previous calculations: the differential affinity of the ligand for a range of RNA-folding intermediates. Significant shifts in the free energy landscape (FEL) due to the ligand occur for transcripts of lengths corresponding to the "decision window," following transcription of the so-called anti-terminator helix. The results suggest how ligand perturbation can stabilize the formation of an intermediate conformation, readily facilitating terminator hairpin formation in the full-length riboswitch.

Basis for ligand discrimination between ON and OFF state riboswitch conformations: the case of the SAM-I riboswitch.

Riboswitches are RNA elements that bind to effector ligands and control gene expression. Most consist of two domains. S-Adenosyl Methionine (SAM) binds the aptamer domain of the SAM-I riboswitch and induces conformational changes in the expression domain to form an intrinsic terminator (transcription OFF state). Without SAM the riboswitch forms the transcription ON state, allowing read-through transcription. The mechanistic link between the SAM/aptamer recognition event and subsequent secondary structure rearrangement by the riboswitch is unclear. We probed for those structural features of the Bacillus subtilis yitJ SAM-I riboswitch responsible for discrimination between the ON and OFF states by SAM. We designed SAM-I riboswitch RNA segments forming "hybrid" structures of the ON and OFF states. The choice of segment constrains the formation of a partial P1 helix, characteristic of the OFF state, together with a partial antiterminator (AT) helix, characteristic of the ON state. For most choices of P1 vs. AT helix lengths, SAM binds with micromolar affinity according to equilibrium dialysis. Mutational analysis and in-line probing confirm that the mode of SAM binding by hybrid structures is similar to that of the aptamer. Altogether, binding measurements and in-line probing are consistent with the hypothesis that when SAM is present, stacking interactions with the AT helix stabilize a partially formed P1 helix in the hybrids. Molecular modeling indicates that continuous stacking between the P1 and the AT helices is plausible with SAM bound. Our findings raise the possibility that conformational intermediates may play a role in ligand-induced aptamer folding.

The impact of a ligand binding on strand migration in the SAM-I riboswitch.

Riboswitches sense cellular concentrations of small molecules and use this information to adjust synthesis rates of related metabolites. Riboswitches include an aptamer domain to detect the ligand and an expression platform to control gene expression. Previous structural studies of riboswitches largely focused on aptamers, truncating the expression domain to suppress conformational switching. To link ligand/aptamer binding to conformational switching, we constructed models of an S-adenosyl methionine (SAM)-I riboswitch RNA segment incorporating elements of the expression platform, allowing formation of an antiterminator (AT) helix. Using Anton, a computer specially developed for long timescale Molecular Dynamics (MD), we simulated an extended (three microseconds) MD trajectory with SAM bound to a modeled riboswitch RNA segment. Remarkably, we observed a strand migration, converting three base pairs from an antiterminator (AT) helix, characteristic of the transcription ON state, to a P1 helix, characteristic of the OFF state. This conformational switching towards the OFF state is observed only in the presence of SAM. Among seven extended trajectories with three starting structures, the presence of SAM enhances the trend towards the OFF state for two out of three starting structures tested. Our simulation provides a visual demonstration of how a small molecule (<500 MW) binding to a limited surface can trigger a large scale conformational rearrangement in a 40 kDa RNA by perturbing the Free Energy Landscape. Such a mechanism can explain minimal requirements for SAM binding and transcription termination for SAM-I riboswitches previously reported experimentally.

A mechanism for S-adenosyl methionine assisted formation of a riboswitch conformation: a small molecule with a strong arm.

The S-adenosylmethionine-1 (SAM-I) riboswitch mediates expression of proteins involved in sulfur metabolism via formation of alternative conformations in response to binding by SAM. Models for kinetic trapping of the RNA in the bound conformation require annealing of nonadjacent mRNA segments during a transcriptional pause. The entropic cost required to bring nonadjacent segments together should slow the folding process. To address this paradox, we performed molecular dynamics simulations on the SAM-I riboswitch aptamer domain with and without SAM, starting with the X-ray coordinates of the SAM-bound RNA. Individual trajectories are 200 ns, among the longest reported for an RNA of this size. We applied principle component analysis (PCA) to explore the global dynamics differences between these two trajectories. We observed a conformational switch between a stacked and nonstacked state of a nonadjacent dinucleotide in the presence of SAM. In the absence of SAM the coordination between a bound magnesium ion and the phosphate of A9, one of the nucleotides involved in the dinucleotide stack, is destabilized. An electrostatic potential map reveals a 'hot spot' at the Mg binding site in the presence of SAM. These results suggest that SAM binding helps to position J1/2 in a manner that is favorable for P1 helix formation.

Protocol for efficient fluorescence 3' end-labeling of native noncoding RNA domains.

Noncoding RNAs (ncRNAs) comprise a class of versatile transcripts that are highly involved in the regulation of a wide range of biological processes. Functional long ncRNAs (> 200 nts in length) often adopt secondary structures that arise co-transcriptionally. To maintain the secondary structure elements as well as preparation homogeneity of such transcripts, native-like conditions should be maintained throughout the in vitro synthesis, purification and chemical tagging processes. In this optimized protocol, we describe a simple method for obtaining homogenous samples followed by chemically tagging the 3' termini of natively-purified structured ncRNA domains that are longer than 200 nts. This protocol replaces traditional hazardous radioactive labeling with fluorescence tagging, and eliminates laborious and time consuming RNA purification and concentration steps and replaces them with straightforward recovery of RNA through centrifugal filtration, preserving the homogeneity and mono-dispersion of the preparations. The protocol provides:•An integrative, simple and straightforward approach for synthesis, purification and labeling of structured ncRNAs whilst maintaining their secondary structure intact.•Replacing hazardous, laborious and time-consuming radioactive labeling of RNA with much simpler fluorescence tagging, thereby facilitating potential downstream applications such as electrophoretic mobility shift assay (EMSA).•A versatile protocol that could be applicable to a wide-range of chemical tags and in principle could be used to label DNA or RNA.

Crystal structures of complexes between aminoglycosides and decoding A site oligonucleotides: role of the number of rings and positive charges in the specific binding leading to miscoding.

The crystal structures of six complexes between aminoglycoside antibiotics (neamine, gentamicin C1A, kanamycin A, ribostamycin, lividomycin A and neomycin B) and oligonucleotides containing the decoding A site of bacterial ribosomes are reported at resolutions between 2.2 and 3.0 A. Although the number of contacts between the RNA and the aminoglycosides varies between 20 and 31, up to eight direct hydrogen bonds between rings I and II of the neamine moiety are conserved in the observed complexes. The puckered sugar ring I is inserted into the A site helix by stacking against G1491 and forms a pseudo base pair with two H-bonds to the Watson-Crick sites of the universally conserved A1408. This central interaction helps to maintain A1492 and A1493 in a bulged-out conformation. All these structures of the minimal A site RNA complexed to various aminoglycosides display crystal packings with intermolecular contacts between the bulging A1492 and A1493 and the shallow/minor groove of Watson-Crick pairs in a neighbouring helix. In one crystal, one empty A site is observed. In two crystals, two aminoglycosides are bound to the same A site with one bound specifically and the other bound in various ways in the deep/major groove at the edge of the A sites.

Towards the discovery of drug-like RNA ligands?

Targeting RNA with small molecule drugs is an area of great potential for therapeutic treatment of infections and possibly genetic and autoimmune diseases. However, a mature set of precedents and established methodology is lacking. The physicochemical properties of RNA raise specific issues and obstacles to development, and contribute to explain the distinct characteristics of natural RNA ligands, including antibiotics. Yet, RNA-targeting strategies are being implemented to reinvigorate antibacterial discovery by using the ribosomal X-ray structures to modify known antibiotics. To exploit further these structures, we suggest the use of existing protein kinase-directed libraries of drug-like compounds to target the A-site of the bacterial ribosome, on the basis of a specific structural hypothesis.

Rational design of inhibitors of HIV-1 TAR RNA through the stabilisation of electrostatic "hot spots".

The targeting of RNA for the design of novel anti-viral compounds has until now proceeded largely without incorporating direct input from structure-based design methodology, partly because of lack of structural data, and complications arising from substrate flexibility. We propose a paradigm to explain the physical mechanism for ligand-induced refolding of trans-activation response element (TAR RNA) from human immunodeficiency virus 1 (HIV-1). Based upon Poisson-Boltzmann analysis of the TAR structure, as bound by a peptide derived from the transcriptional activator protein, Tat, our hypothesis shows that two specific electrostatic interactions are necessary to stabilise the conformation. This result contradicts the belief that a single argininamide residue is responsible for stabilising the TAR fold, as well as the conventional wisdom that electrostatic interactions with RNA are non-specific or dominated by phosphates. We test this hypothesis by using NMR and computational methods to model the interaction of a series of novel inhibitors of the in vitro RNA-binding activities for a peptide derived from Tat. A subset of inhibitors, including the bis-guanidine compound rbt203 and its analogues, induce a conformation in TAR similar to that brought about by the protein. Comparison of the interactions of two of these ligands with the RNA and structure-activity relationships observed within the compound series, confirm the importance of the two specific electrostatic interactions in the stabilisation of the Tat-bound RNA conformation. This work illustrates how the use of medicinal chemistry and structural analysis can provide a rational basis for prediction of ligand-induced conformational change, a necessary step towards the application of structure-based methods in the design of novel RNA or protein-binding drugs.

Structure-based drug design targeting an inactive RNA conformation: exploiting the flexibility of HIV-1 TAR RNA.

The targeting of RNA for the design of novel anti-viral compounds represents an area of vast potential. We have used NMR and computational methods to model the interaction of a series of synthetic inhibitors of the in vitro RNA binding activities of a peptide derived from the transcriptional activator protein, Tat, from human immunodeficiency virus type 1. Inhibition has been measured through the monitering of fluorescence resonance energy transfer between fluorescently labeled peptide and RNA components. A series of compounds containing a bi-aryl heterocycle as one of the three substituents on a benzylic scaffold, induce a novel, inactive TAR conformation by stacking between base-pairs at the site of a three-base bulge within TAR. The development of this series resulted in an enhancement in potency (with Ki < 100 nM in an in vitro assay) and the removal of problematic guanidinium moieties. Ligands from this series can act as inhibitors of Tat-induced transcription in a cell-free system. This study validates the drug design strategy of using a ligand to target the RNA receptor in a non-functional conformation.

Structural basis for contrasting activities of ribosome binding thiazole antibiotics.

Thiostrepton and micrococcin inhibit protein synthesis by binding to the L11 binding domain (L11BD) of 23S ribosomal RNA. The two compounds are structurally related, yet they produce different effects on ribosomal RNA in footprinting experiments and on elongation factor-G (EF-G)-dependent GTP hydrolysis. Using NMR and an assay based on A1067 methylation by thiostrepton-resistance methyltransferase, we show that the related thiazoles, nosiheptide and siomycin, also bind to this region. The effect of all four antibiotics on EF-G-dependent GTP hydrolysis and EF-G-GDP-ribosome complex formation was studied. Our NMR and biochemical data demonstrate that thiostrepton, nosiheptide, and siomycin share a common profile, which differs from that of micrococcin. We have generated a three-dimensional (3D) model for the interaction of thiostrepton with L11BD RNA. The model rationalizes the differences between micrococcin and the thiostrepton-like antibiotics interacting with L11BD.

Targeting the A site RNA of the Escherichia coli ribosomal 30 S subunit by 2'-O-methyl oligoribonucleotides: a quantitative equilibrium dialysis binding assay and differential effects of aminoglycoside antibiotics.

The bacterial ribosome comprises 30 S and 50 S ribonucleoprotein subunits, contains a number of binding sites for known antibiotics and is an attractive target for selection of novel antibacterial agents. On the 30 S subunit, for example, the A site (aminoacyl site) close to the 3'-end of 16 S rRNA is highly important in the decoding process. Binding by some aminoglycoside antibiotics to the A site leads to erroneous protein synthesis and is lethal for bacteria. We targeted the A site on purified 30 S ribosomal subunits from Escherichia coli with a set of overlapping, complementary OMe (2'-O-methyl) 10-mer oligoribonucleotides. An equilibrium dialysis technique was applied to measure dissociation constants of these oligonucleotides. We show that there is a single high-affinity region, spanning from A1493 to C1510 (Kd, 29-130 nM), flanked by two lower-affinity regions, within a span from U1485 to G1516 (Kd, 310-4300 nM). Unexpectedly, addition of the aminoglycoside antibiotic paromomycin (but not hygromycin B) caused a dose-dependent increase of up to 7.5-fold in the binding of the highest affinity 10-mer 1493 to 30 S subunits. Oligonucleotides containing residues complementary to A1492 and/or A1493 showed particularly marked stimulation of binding by paromomycin. The results are consistent with high-resolution structures of antibiotic binding to the A site and with greater accessibility of residues of A1492 and A1493 upon paromomycin binding. 10-mer 1493 binding is thus a probe of the conformational switch to the 'closed' conformation triggered by paromomycin that is implicated in the discrimination by 30 S subunits of cognate from non-cognate tRNA and the translational misreading caused by paromomycin. Finally, we show that OMe oligonucleotides targeted to the A site are moderately good inhibitors of in vitro translation and that there is a limited correlation of inhibition activity with binding strength to the A site.

Linking aptamer-ligand binding and expression platform folding in riboswitches: prospects for mechanistic modeling and design.

The power of riboswitches in regulation of bacterial metabolism derives from coupling of two characteristics: recognition and folding. Riboswitches contain aptamers, which function as biosensors. Upon detection of the signaling molecule, the riboswitch transduces the signal into a genetic decision. The genetic decision is coupled to refolding of the expression platform, which is distinct from, although overlapping with, the aptamer. Early biophysical studies of riboswitches focused on recognition of the ligand by the aptamer-an important consideration for drug design. A mechanistic understanding of ligand-induced riboswitch RNA folding can further enhance riboswitch ligand design, and inform efforts to tune and engineer riboswitches with novel properties. X-ray structures of aptamer/ligand complexes point to mechanisms through which the ligand brings together distal strand segments to form a P1 helix. Transcriptional riboswitches must detect the ligand and form this P1 helix within the timescale of transcription. Depending on the cell's metabolic state and cellular environmental conditions, the folding and genetic outcome may therefore be affected by kinetics of ligand binding, RNA folding, and transcriptional pausing, among other factors. Although some studies of isolated riboswitch aptamers found homogeneous, prefolded conformations, experimental, and theoretical studies point to functional and structural heterogeneity for nascent transcripts. Recently it has been shown that some riboswitch segments, containing the aptamer and partial expression platforms, can form binding-competent conformers that incorporate an incomplete aptamer secondary structure. Consideration of the free energy landscape for riboswitch RNA folding suggests models for how these conformers may act as transition states-facilitating rapid, ligand-mediated aptamer folding.

RNA as a drug target.

In the antiviral and antibacterial area, increasing drug resistance means that there is an ever growing need for novel approaches towards structures and mechanisms which avoid the current problems. The huge increase in high resolution structural data is set to make a dramatic impact on targeting RNA as a drug target. The examples of the RNA binding antibiotics, particularly, the totally synthetic oxazolidinones, should help persuade the skceptics that clinically useful, selective drugs can be obtained from targeting RNA directly.

Conformational heterogeneity of the SAM-I riboswitch transcriptional ON state: a chaperone-like role for S-adenosyl methionine.

Riboswitches are promising targets for the design of novel antibiotics and engineering of portable genetic regulatory elements. There is evidence that variability in riboswitch properties allows tuning of expression for genes involved in different stages of biosynthetic pathways by mechanisms that are not currently understood. Here, we explore the mechanism for tuning of S-adenosyl methionine (SAM)-I riboswitch folding. Most SAM-I riboswitches function at the transcriptional level by sensing the cognate ligand SAM. SAM-I riboswitches orchestrate the biosynthetic pathways of cysteine, methionine, SAM, and so forth. We use base-pair probability predictions to examine the secondary-structure folding landscape of several SAM-I riboswitch sequences. We predict different folding behaviors for different SAM-I riboswitch sequences. We identify several "decoy" base-pairing interactions involving 5' riboswitch residues that can compete with the formation of a P1 helix, a component of the ligand-bound "transcription OFF" state, in the absence of SAM. We hypothesize that blockage of these interactions through SAM contacts contributes to stabilization of the OFF state in the presence of ligand. We also probe folding patterns for a SAM-I riboswitch RNA using constructs with different 3' truncation points experimentally. Folding was monitored through fluorescence, susceptibility to base-catalyzed cleavage, nuclear magnetic resonance, and indirectly through SAM binding. We identify key decision windows at which SAM can affect the folding pathway towards the OFF state. The presence of decoy conformations and differential sensitivities to SAM at different transcript lengths is crucial for SAM-I riboswitches to modulate gene expression in the context of global cellular metabolism.

Strategies for the design of RNA-binding small molecules.

Bacterial ribosomal RNA is the target of clinically important antibiotics, while biologically important RNAs in viral and eukaryotic genomes present a range of potential drug targets. The physicochemical properties of RNA present difficulties for medicinal chemistry, particularly when oral availability is needed. Peptidic ligands and analysis of their RNA-binding properties are providing insight into RNA recognition. RNA-binding ligands include far more chemical classes than just aminoglycosides. Chemical functionalities from known RNA-binding small molecules are being exploited in fragment- and ligand-based projects. While targeting of RNA for drug design is very challenging, continuing advances in our understanding of the principles of RNA-ligand interaction will be necessary to realize the full potential of this class of targets.

Design and implementation of an ribonucleic acid (RNA) directed fragment library.

The design of RNA binding ligands is complicated by issues of specificity, target flexibility, and the tractability of known RNA inhibitors toward chemical derivitization. To address these difficulties, an RNA-directed fragment compound library is presented. We began with an analysis of 120 small molecules with reported RNA-binding activity. Calculated physical and chemical properties for the RNA ligands are comparable to those of ligands for established protein drug targets. To ensure that our library contained RNA-binding functionalities that might not be detected by the above comparisons, 114 fragment compounds were purchased on the basis of similarity to substructures of RNA ligands. Five "hits" were identified for the decoding site from the bacterial ribosome by NMR. These included fragments derived from A-site binding ligands but also compounds not previously identified as A-site binders. Hits generated in this manner can be used to probe the interaction surface of RNA and its conformational plasticity, facilitating structure-based optimization.

Design and characterization of libraries of molecular fragments for use in NMR screening against protein targets.

We have designed four generations of a low molecular weight fragment library for use in NMR-based screening against protein targets. The library initially contained 723 fragments which were selected manually from the Available Chemicals Directory. A series of in silico filters and property calculations were developed to automate the selection process, allowing a larger database of 1.79 M available compounds to be searched for a further 357 compounds that were added to the library. A kinase binding pharmacophore was then derived to select 174 kinase-focused fragments. Finally, an additional 61 fragments were selected to increase the number of different pharmacophores represented within the library. All of the fragments added to the library passed quality checks to ensure they were suitable for the screening protocol, with appropriate solubility, purity, chemical stability, and unambiguous NMR spectrum. The successive generations of libraries have been characterized through analysis of structural properties (molecular weight, lipophilicity, polar surface area, number of rotatable bonds, and hydrogen-bonding potential) and by analyzing their pharmacophoric complexity. These calculations have been used to compare the fragment libraries with a drug-like reference set of compounds and a set of molecules that bind to protein active sites. In addition, an analysis of the overall results of screening the library against the ATP binding site of two protein targets (HSP90 and CDK2) reveals different patterns of fragment binding, demonstrating that the approach can find selective compounds that discriminate between related binding sites.

A conserved RNA structure within the HCV IRES eIF3-binding site.

The hepatitis C virus (HCV) internal ribosome entry site (IRES) is recognized specifically by the small ribosomal subunit and eukaryotic initiation factor 3 (eIF3) before viral translation initiation. Using extensive mutagenesis and structure probing analysis, we show that the eIF3-binding domain of the HCV IRES contains an internal loop structure (loop IIIb) and an adjacent mismatched helix that are important for IRES-dependent initiation of translation. NMR studies reveal a unique three-dimensional structure for this internal loop that is conserved between viral isolates of varying primary sequence in this region. These data indicate that internal loop IIIb may be an attractive target for structure-based design of new antiviral agents.