Structure-based discovery of ligands targeted to the RNA double helix

RNA-binding of Ru(II) complexes [Ru(phen)2(7-OCH3-dppz)]2+ and [Ru(phen)2(7-NO2-dppz)]2+: The former serves as a molecular "light switch" for poly(A)•poly(U)

To further determine the factors that affect the binding properties of ruthenium(II) polypyridine complexes with RNA duplex and to find excellent RNA-binding agents, the binding properties of ruthenium(II) complexes [Ru(phen)₂(7-OCH₃-dppz)]²⁺ (Ru1, phen = 1,10-phenan- throline, 7-OCH₃-dppz = 7-methoxy-dipyrido-[3,2-a,2',3'-c]-phenazine) and [Ru(phen)₂(7-NO₂- dppz)]²⁺ (Ru2, 7-NO₂-dppz = 7-nitro-dipyrido-[3,2-a,2',3'-c]-phenazine) with RNA poly(A)•poly(U) duplex have been investigated by spectroscopic methods and viscosity measurements in this work. The results show that complexes Ru1 and Ru2 bind to poly(A)•poly(U) through intercalation and the binding affinity between Ru2 and poly(A)•poly(U) is greater than that of Ru1. Thermal denaturation experiments suggest that both ruthenium(II) complexes exhibit a significant stabilizing effect on poly(A)•poly(U) duplex. Moreover, fluorescence emission spectra exhibit that, deviating from Ru2, Ru1 exhibits a "light switch" effect for poly(A)•poly(U). This effect can be observed by the naked eye under UV light and adjusted by pH, meaning that Ru1 may act as a reversible pH controlled molecular "light switch". The results obtained in this work will contribute to our understanding of the significant influence of the intercalative ligand substituent effect in the binding process of ruthenium(II) complexes with RNA duplex.

How 'Protein-Docking' Translates into the New Emerging Field of Docking Small Molecules to Nucleic Acids?

In this review, we retraced the '40-year evolution' of molecular docking algorithms. Over the course of the years, their development allowed to progress from the so-called 'rigid-docking' searching methods to the more sophisticated 'semi-flexible' and 'flexible docking' algorithms. Together with the advancement of computing architecture and power, molecular docking's applications also exponentially increased, from a single-ligand binding calculation to large screening and polypharmacology profiles. Recently targeting nucleic acids with small molecules has emerged as a valuable therapeutic strategy especially for cancer treatment, along with bacterial and viral infections. For example, therapeutic intervention at the mRNA level allows to overcome the problematic of undruggable proteins without modifying the genome. Despite the promising therapeutic potential of nucleic acids, molecular docking programs have been optimized mostly for proteins. Here, we have analyzed literature data on nucleic acid to benchmark some of the widely used docking programs. Finally, the comparison between proteins and nucleic acid targets docking highlighted similarity and differences, which are intrinsically related to their chemical and structural nature.

Molecular Docking of Biologically Active Substances to Double Helical Nucleic Acids

Molecular docking of ligands to DNA-targets is of great importance for the design of new anticancer drugs. Unfortunately, most docking programs were developed for protein-ligand docking which raises a question about their applicability for the DNA-ligand docking. In this study, the popular docking programs AutoDock Vina, AutoDock4 and AutoDock3 were compared for a test set of 50 DNA-ligand complexes taken from the Nucleic Acid Database. It was shown that the version 3.05 of the AutoDock program was the most successful in reproducing the structures of intercalation and minor-groove complexes. The program AutoDock4 was able to re-dock to within 2 Å RMSD most of the intercalation complexes of the test set, but showed poor performance for minor groove binders. While Vina, on the contrary, failed to construct six intercalation complexes of the test set, but showed satisfactory results for DNA-ligand minor-groove complexes when small search space was used.

¹⁹F-Site-Specific-Labeled Nucleotides for Nucleic Acid Structural Analysis by NMR

Naturally occurring RNA lacks fluorine-19 ((19)F), thus, their specifically fluorinated counterparts are particularly well suited to noninvasively monitoring the dynamic conformational properties and ligand-binding interactions of the RNA. For nuclear magnetic resonance (NMR) spectroscopy, (19)F-NMR of fluorine-substituted RNA provides an attractive, site-specific probe for structure determination in solution. Advantages of (19)F include high NMR sensitivity (83% of (1)H), high natural abundance (100%), and the extreme sensitivity of (19)F to the chemical environment leading to a large range of chemical shifts. The preparation of base-substituted 2-fluoropurine and 5-fluoropyrimidine 5'-triphosphates (2F-ATP/5F-CTP/5F-UTP) can be carried out using efficient enzymatic synthesis methods. Both pyrimidine analogs, 5-fluorouridine and 5-fluorocytidine, as well as, 2-fluoroadenosine are readily incorporated into RNA transcribed in vitro using T7 RNA polymerase.

Charge reduction and thermodynamic stabilization of substrate RNAs inhibit RNA editing

African trypanosomes cause a parasitic disease known as sleeping sickness. Mitochondrial transcript maturation in these organisms requires a RNA editing reaction that is characterized by the insertion and deletion of U-nucleotides into otherwise non-functional mRNAs. Editing represents an ideal target for a parasite-specific therapeutic intervention since the reaction cycle is absent in the infected host. In addition, editing relies on a macromolecular protein complex, the editosome, that only exists in the parasite. Therefore, all attempts to search for editing interfering compounds have been focused on molecules that bind to proteins of the editing machinery. However, in analogy to other RNA-driven biochemical pathways it should be possible to stall the reaction by targeting its substrate RNAs. Here we demonstrate inhibition of editing by specific aminoglycosides. The molecules bind into the major groove of the gRNA/pre-mRNA editing substrates thereby causing a stabilization of the RNA molecules through charge compensation and an increase in stacking. The data shed light on mechanistic details of the editing process and identify critical parameters for the development of new trypanocidal compounds.

Synthesis and properties of vitamin E analog-conjugated neomycin for delivery of RNAi drugs to liver cells

RNA interference (RNAi) is a promising tool to regulate gene expression by external double stranded RNAs (dsRNAs) such as siRNAs. As an efficient method to deliver siRNAs to liver cells, we propose a novel strategy using vitamin E (VE)-conjugated neomycin derivatives. With the aim of delivering RNAi-based drugs to liver cells, several tripod-type and prodrug-type neomycin derivatives were synthesized, all of which were thermodynamically stabilized RNA duplexes. The prodrug-type derivative 7 and the tripod-type derivative 10 were delivered to liver cancer cells and successfully induced RNAi activity. These results indicated the potential use of natural aminoglycosides as carriers of RNAi drugs.

Structure-based virtual screening for the identification of RNA-binding ligands

Structure-based virtual screening exploits the 3D structure of the target as a template for the discovery of new ligands. It is a key method for hit discovery and was originally developed for protein targets. Recently, this method has also been applied to RNA targets. This chapter gives an overview of this method and its application in the context of ligand discovery for RNA. In addition, it describes in detail how to conduct virtual screening for RNA targets, making use of software that is free for noncommercial use. Some advice on how to avoid common pitfalls in virtual screening is also given.

Novel insights of structure-based modeling for RNA-targeted drug discovery

Substantial progress in RNA biology highlights the importance of RNAs (e.g., microRNAs) in diseases and the potential of targeting RNAs for drug discovery. However, the lack of RNA-specific modeling techniques demands the development of new tools for RNA-targeted rational drug design. Herein, we implemented integrated approaches of accurate RNA modeling and virtual screening for RNA inhibitor discovery with the most comprehensive evaluation to date of five docking and 11 scoring methods. For the first time, statistical analysis was heavily employed to assess the significance of our predictions. We found that GOLD:GOLD Fitness and rDock:rDock_solv could accurately predict the RNA ligand poses, and ASP rescoring further improved the ranking of ligand binding poses. Due to the weak correlations (R(2) < 0.3) of existing scoring with experimental binding affinities, we implemented two new RNA-specific scoring functions, iMDLScore1 and iMDLScore2, and obtained better correlations with R(2) = 0.70 and 0.79, respectively. We also proposed a multistep virtual screening approach and demonstrated that rDock:rDock_solv together with iMDLScore2 rescoring obtained the best enrichment on the flexible RNA targets, whereas GOLD:GOLD Fitness combined with rDock_solv rescoring outperformed other methods for rigid RNAs. This study provided practical strategies for RNA modeling and offered new insights into RNA-small molecule interactions for drug discovery.

Virtual Screening for RNA-Interacting Small Molecules

Computational virtual screening is useful and powerful strategy for rapid discovery of small biologically active molecules in the early stage of drug discovery. The discovery of a broad range of important biological processes involved by RNA has increased the importance of RNA as a new drug target. To apply structure-based virtual screening methods to the discovery of RNA-binding ligands, many RNA 3D structure prediction programs and optimized docking algorithms have been developed. In this chapter, a number of successful cases of virtual screening targeting RNA will be introduced.

Target Flexibility in RNA-Ligand Docking Modeled by Elastic Potential Grids

The highly flexible nature of RNA provides a formidable challenge for structure-based drug design approaches that target RNA. We introduce an approach for modeling target conformational changes in RNA-ligand docking based on potential grids that are represented as elastic bodies using Navier's equation. This representation provides an accurate and efficient description of RNA-ligand interactions even in the case of a moving RNA structure. When applied to a data set of 17 RNA-ligand complexes, filtered out of the largest validation data set used for RNA-ligand docking so far, the approach is twice as successful as docking into an apo structure and still half as successful as redocking to the holo structure. The approach allows considering RNA movements of up to 6 Å rmsd and is based on a uniform and robust parametrization of the properties of the elastic potential grids, so that the approach is applicable to different RNA-ligand complex classes.

Synthesis of oligodiaminosaccharides having α-glycoside bonds and their interactions with oligonucleotide duplexes

Syntheses of the novel oligodiaminosaccharides, α-(1→4)-linked-2,6-diamino-2,6-dideoxy-D-glucopyranose oligomers, and their interactions with nucleic acid duplexes DNA-DNA, RNA-RNA, and DNA-RNA are described. Monomers to tetramers of oligodiaminoglucose derivatives having α-glycosyl bonds were successfully synthesized using a chain elongation cycle including glycosylation reactions of a 6-phthalimide glycosyl donor. UV melting experiments for a variety of nucleic acid duplexes in the absence and presence of the oligodiaminosaccharides were performed. The synthesized oligodiaminosaccharides exhibited notable thermodynamic stabilization effects on A-type RNA-RNA and DNA-RNA duplexes, whereas B-type DNA-DNA duplexes were not stabilized by the synthesized oligodiaminosaccharides. Among the oligodiaminosaccharides, the tetramer exhibited the highest ability to stabilize A-type duplexes, and the increase in T(m) values induced by the tetramer were higher than those induced by neomycin B and tobramycin, which are known aminoglycosides having ability to bind and stabilize a variety of RNA molecules. CD spectrometry experiments revealed that the oligodiaminosaccharides caused small structural changes in RNA-RNA duplexes, whereas no appreciable changes were observed in the structure of DNA-DNA duplexes. ITC (isothermal titration calorimetry) experiments demonstrated that the amount of heat generated by the interaction between RNA-RNA duplexes and the tetradiaminosaccharides was approximately double that generated by that between DNA-DNA duplexes and the tetradiaminosaccharides. These results strongly suggested the existence of an A-type nucleic acid specific-binding mode of the oligodiaminosaccharides, which bind to these duplexes and cause small structural changes.

Binding-interaction prediction of RNA-binding ligands

RNA molecules are involved in a wide range of biological processes and have been recognized as very important therapeutic targets. Mainly owing to the scarcity of information and experimental studies, the application of computational approaches and, in particular, of docking methodologies in the RNA field has developed slowly. However, in recent years the docking of RNA-binding ligands has experienced significant expansion. This article focuses attention on the docking of RNA-binding ligands, analyzing the development of RNA-docking approaches, the reliability of the docking methods and, finally, evaluating the results of docking-based virtual screening studies reported in the literature.

A Discovery Funnel for Nucleic Acid Binding Drug Candidates

Computational approaches are becoming increasingly popular for the discovery of drug candidates against a target of interest. Proteins have historically been the primary targets of many virtual screening efforts. While in silico screens targeting proteins has proven successful, other classes of targets, in particular DNA, remain largely unexplored using virtual screening methods. With the realization of the functional importance of many non-cannonical DNA structures such as G-quadruplexes, increased efforts are underway to discover new small molecules that can bind selectively to DNA structures. Here, we describe efforts to build an integrated in silico and in vitro platform for discovering compounds that may bind to a chosen DNA target. Millions of compounds are initially screened in silico for selective binding to a particular structure and ranked to identify several hundred best hits. An important element of our strategy is the inclusion of an array of possible competing structures in the in silico screen. The best hundred or so hits are validated experimentally for binding to the actual target structure by a high-throughput 96-well thermal denaturation assay to yield the top ten candidates. Finally, these most promising candidates are thoroughly characterized for binding to their DNA target by rigorous biophysical methods, including isothermal titration calorimetry, differential scanning calorimetry, spectroscopy and competition dialysis.This platform was validated using quadruplex DNA as a target and a newly discovered quadruplex binding compound with possible anti-cancer activity was discovered. Some considerations when embarking on virtual screening and in silico experiments are also discussed.

Molecular recognition of RNA: challenges for modelling interactions and plasticity

There is growing interest in molecular recognition processes of RNA because of RNA's widespread involvement in biological processes. Computational approaches are increasingly used for analysing and predicting binding to RNA, fuelled by encouraging progress in developing simulation, free energy and docking methods for nucleic acids. These developments take into account challenges regarding the energetics of RNA-ligand binding, RNA plasticity, and the presence of water molecules and ions in the binding interface. Accordingly, we will detail advances in force field and scoring function development for molecular dynamics (MD) simulations, free energy computations and docking calculations of nucleic acid complexes. Furthermore, we present methods that can detect moving parts within RNA structures based on graph-theoretical approaches or normal mode analysis (NMA). As an example of the successful use of these developments, we will discuss recent structure-based drug design approaches that focus on the bacterial ribosomal A-site RNA as a drug target.

Calorimetric and spectroscopic studies of aminoglycoside binding to AT-rich DNA triple helices

Calorimetric and fluorescence techniques were used to characterize the binding of aminoglycosides-neomycin, paromomycin, and ribostamycin, with 5'-dA(12)-x-dT(12)-x-dT(12)-3' intramolecular DNA triplex (x = hexaethylene glycol) and poly(dA).2poly(dT) triplex. Our results demonstrate the following features: (1) UV thermal analysis reveals that the T(m) for triplex decreases with increasing pH value in the presence of neomycin, while the T(m) for the duplex remains unchanged. (2) The binding affinity of neomycin decreases with increased pH, although there is an increase in observed binding enthalpy. (3) ITC studies conducted in two buffers (sodium cacodylate and MOPS) yield the number of protonated drug amino groups (Deltan) as 0.29 and 0.40 for neomycin and paromomycin interaction with 5'-dA(12)-x-dT(12)-x-dT(12)-3', respectively. (4) The specific heat capacity change (DeltaC(p)) determined by ITC studies is negative, with more negative values at lower salt concentrations. From 100 mM to 250 mM KCl, the DeltaC(p) ranges from -402 to -60 cal/(mol K) for neomycin. At pH 5.5, a more positive DeltaC(p) is observed, with a value of -98 cal/(mol K) at 100 mM KCl. DeltaC(p) is not significantly affected by ionic strength. (5) Salt dependence studies reveal that there are at least three amino groups of neomycin participating in the electrostatic interactions with the triplex. (6) FID studies using thiazole orange were used to derive the AC(50) (aminoglycoside concentration needed to displace 50% of the dye from the triplex) values. Neomycin shows a seven fold higher affinity than paromomycin and eleven fold higher affinity than ribostamycin at pH 6.8. (7) Modeling studies, consistent with UV and ITC results, show the importance of an additional positive charge in triplex recognition by neomycin. The modeling and thermodynamic studies indicate that neomycin binding to the DNA triplex depends upon significant contributions from charge as well as shape complementarity of the drug to the DNA triplex Watson-Hoogsteen groove.

Discovery of novel triple helical DNA intercalators by an integrated virtual and actual screening platform

Virtual Screening is an increasingly attractive way to discover new small molecules with potential medicinal value. We introduce a novel strategy that integrates use of the molecular docking software Surflex with experimental validation by the method of competition dialysis. This integrated approach was used to identify ligands that selectively bind to the triplex DNA poly(dA)-[poly(dT)]₂. A library containing ∼2 million ligands was virtually screened to identify compounds with chemical and structural similarity to a known triplex intercalator, the napthylquinoline MHQ-12. Further molecular docking studies using compounds with high structural similarity resulted in two compounds that were then demonstrated by competition dialysis to have a superior affinity and selectivity for the triplex nucleic acid than MHQ-12. One of the compounds has a different chemical backbone than MHQ-12, which demonstrates the ability of this strategy to ‘scaffold hop’ and to identify small molecules with novel binding properties. Biophysical characterization of these compounds by circular dichroism and thermal denaturation studies confirmed their binding mode and selectivity. These studies provide a proof-of-principle for our integrated screening strategy, and suggest that this platform may be extended to discover new compounds that target therapeutically relevant nucleic acid morphologies.

Molecular docking of intercalators and groove-binders to nucleic acids using Autodock and Surflex

The molecular docking tools Autodock and Surflex accurately reproduce the crystallographic structures of a collection of small molecule ligands that have been shown to bind nucleic acids. Docking studies were performed with the intercalators daunorubicin and ellipticine and the minor groove binders distamycin and pentamidine. Autodock and Surflex dock daunorubicin and distamycin to their nucleic acid targets within a resolution of approximately 2 A, which is similar to the limit of the crystal structure resolution. However, for the top ranked poses, Autodock and Surflex both dock ellipticine into the correct site but in a different orientation compared to the crystal structure. This appears not only to be partly related to the symmetry of the target nucleic acid, as ellipticine is able to dock from either side of the intercalation site, but also due to the shape of the ligand and docking accuracy. Surflex docks pentamidine in a symmetrically equivalent orientation relative to the crystal structure, while Autodock was able to dock this molecule in the original orientation. In the case of the Surflex docking of pentamidine, the initial rmsd is misleading, given the symmetrical structure of pentamidine. Importantly, the ranking functions of both of these programs are able to return a top pose within approximately 2 A rmsd for daunorubicin, distamycin, and pentamidine and approximately 3 A rmsd for ellipticine compared to their respective crystal structures. Some docking challenges and potential pitfalls are explored, such as the importance of hydrogen treatment on ligands as well as the scoring functions of Autodock and Surflex. Overall for this set of complexes, Surflex is preferred over Autodock for virtual screening, as although the results are comparable, Surflex has significantly faster performance and ease of use under the optimal software conditions tested. These experiments show that molecular docking techniques can be successfully extended to include nucleic acid targets, a finding which has important implications for virtual screening applications and in the design of new small molecules to target therapeutically relevant morphologies of nucleic acids.

Docking to RNA via root-mean-square-deviation-driven energy minimization with flexible ligands and flexible targets

Structure-based drug design is now well-established for proteins as a key first step in the lengthy process of developing new drugs. In many ways, RNA may be a better target to treat disease than a protein because it is upstream in the translation pathway, so inhibiting a single mRNA molecule could prevent the production of thousands of protein gene products. Virtual screening is often the starting point for structure-based drug design. However, computational docking of a small molecule to RNA seems to be more challenging than that to protein due to the higher intrinsic flexibility and highly charged structure of RNA. Previous attempts at docking to RNA showed the need for a new approach. We present here a novel algorithm using molecular simulation techniques to account for both nucleic acid and ligand flexibility. In this approach, with both the ligand and the receptor permitted some flexibility, they can bind one another via an induced fit, as the flexible ligand probes the surface of the receptor. A possible ligand can explore a low-energy path at the surface of the receptor by carrying out energy minimization with root-mean-square-distance constraints. Our procedure was tested on 57 RNA complexes (33 crystal and 24 NMR structures); this is the largest data set to date to reproduce experimental RNA binding poses. With our procedure, the lowest-energy conformations reproduced the experimental binding poses within an atomic root-mean-square deviation of 2.5 A for 74% of tested complexes.

Towards the development of universal, fast and highly accurate docking/scoring methods: a long way to go

Accelerating the drug discovery process requires predictive computational protocols capable of reducing or simplifying the synthetic and/or combinatorial challenge. Docking-based virtual screening methods have been developed and successfully applied to a number of pharmaceutical targets. In this review, we first present the current status of docking and scoring methods, with exhaustive lists of these. We next discuss reported comparative studies, outlining criteria for their interpretation. In the final section, we describe some of the remaining developments that would potentially lead to a universally applicable docking/scoring method.

Guanidinoneomycin B recognition of an HIV-1 RNA helix

Aminoglycoside antibiotics are small-molecule drugs that bind RNA. The affinity and specificity of aminoglycoside binding to RNA can be increased through chemical modification, such as guanidinylation. Here, we report the binding of guanidinoneomycin B (GNB) to an RNA helix from the HIV-1 frameshift site. The binding of GNB increases the melting temperature (T(m)) of the frameshift-site RNA by at least 10 degrees C, to a point at which a melting transition is not even observed in 2 M urea. A structure of the complex was obtained by using multidimensional heteronuclear NMR spectroscopic methods. We also used a novel paramagnetic-probe assay to identify the site of GNB binding to the surface of the RNA. GNB makes major-groove contacts to two sets of Watson-Crick bases and is in van der Waals contact with a highly structured ACAA tetraloop. Rings I and II of GNB fit into the major groove and form the binding interface with the RNA, whereas rings III and IV are exposed to the solvent and disordered. The binding of GNB causes a broadening of the major groove across the binding site.

Synthesis and properties of 4'-ThioDNA: unexpected RNA-like behavior of 4'-ThioDNA

The synthesis and properties of fully modified 4′-thioDNAs, oligonucleotides consisting of 2′-deoxy-4′-thionucleosides, were examined. In addition to the known literature properties (preferable hybridization with RNA and resistance to endonuclease hydrolysis), we also observed higher resistance of 4′-thioDNA to 3′-exonuclease cleavage. Furthermore, we found that fully modified 4′-thioDNAs behaved like RNA molecules in their hybridization properties and structural aspect, at least in the case of the 4′-thioDNA duplex. This observation was confirmed by experiments using groove binders, in which a 4′-thioDNA duplex interacts with an RNA major groove binder, lividomycin A, but not with DNA groove binders, to give an increase in its thermal stability. Since a 4′-thioDNA duplex competitively inhibited the hydrolysis of an RNA duplex by RNase V₁, it was not only the physical properties but also this biological data suggested that a 4′-thioDNA duplex has an RNA-like structure.

Aminoglycoside-Quinacridine Conjugates: Towards Recognition of the P6.1 Element of Telomerase RNA

A modular synthesis has been developed which allows easy and rapid attachment of one or two aminoglycoside units to a quinacridine intercalator, thereby leading to monomeric and dimeric conjugates. Melting temperature (Tm) experiments show that the tobramycin dimeric conjugate TD1 exhibits strong binding to the P6.1 element of human telomerase RNA. By contrast, tobramycin alone is much less efficient and the monomeric compound TM1 elicits a poor binding ability. Monitoring of the interaction by an electrophoretic mobility shift assay shows a 1:1 stoichiometry for the binding of the dimeric compound to the hairpin structure and confirms the lower affinity for a control duplex. Protection experiments with RNase T1 indicate interaction of the drug both in the stem and in the loop of the hairpin. Taken together, the data suggest a binding of TD1 inside the hairpin at the stem-loop junction. The same trends are observed with paromomycin and kanamycin analogues but with a lower affinity.

Discovery of Ligands by a Combination of Computational and NMR-Based Screening: RNA as an Example Target

NMR for screening of knowledge-based focused libraries of compounds provides an efficient, cost-effective method to develop promising drug leads that target functionally important RNA structures. A knowledge-based focused library may be constructed from virtual (i.e., computational) screening of commercial or proprietary databases of available compounds for binding to the three-dimensional structure of a selected RNA target. Alternatively, the library may be constructed from compounds with properties deemed desirable, e.g., molecular moiety commonly found in drugs or known to bind RNA. The library ideally should be composed of small water-soluble, nonpeptide, nonnucleotide organic compounds. Various simple, robust NMR experiments are described that enable experimental screening of such a library for binding to a selected RNA structure. Some of the NMR experiments enable rapid mapping of the interaction site on the RNA to verify that the targeted structure is hit rather than the double helical region or a commonly occurring tetraloop. Other experiments enable elucidation of the ligand's binding moiety. Of course, any compounds thus identified should represent promising scaffolds suitable for easy chemical modification to enhance their pharmaceutical properties for subsequent drug development.

A structure-based strategy to identify new molecular scaffolds targeting the bacterial ribosomal A-site

The need for novel antibiotics is widely recognized. A well validated target of antibiotics is the bacterial ribosome. Recent X-ray structures of the ribosome bound to antibiotics have shed new light on the binding sites of these antibiotics, providing fresh impetus for structure-based strategies aiming at identifying new ribosomal ligands. In that respect, the ribosomal decoding region of the aminoacyl-tRNA acceptor site (A-site) is of particular interest because oligonucleotide model systems of this site are available for crystallography, NMR and compound binding assays. This work presents how these different resources can be combined in a hierarchical screening strategy which has led to the identification of new A-site ligands. The approach exploits an X-ray structure of the A-site against which large and diverse libraries of compounds were computationally docked. The complementarity of the compounds to the A-site was assessed using a scoring function specifically calibrated for RNA targets. Starting from approximately 1 million compounds, the computational selection of candidate ligands allowed us to focus the experimental work on 129 compounds, 34 of which showed affinity for the A-site in a FRET-based binding assay. NMR experiments confirmed binding to the A-site for some compounds. For the most potent compound in the FRET assay, a tentative binding mode is suggested, which is compatible with the NMR data and the limited SAR in this series. Overall, the results validate the screening strategy.

Recognition of Planar and Nonplanar Ligands in the Malachite Green-RNA Aptamer Complex

Ribonucleic acids are an attractive drug target owing to their central role in many pathological processes. Notwithstanding this potential, RNA has only rarely been successfully targeted with novel drugs. The difficulty of targeting RNA is at least in part due to the unusual mode of binding found in most small-molecule-RNA complexes: the ligand binding pocket of the RNA is largely unstructured in the absence of ligand and forms a defined structure only with the ligand acting as scaffold for folding. Moreover, electrostatic interactions between RNA and ligand can also induce significant changes in the ligand structure due to the polyanionic nature of the RNA. Aptamers are ideal model systems to study these kinds of interactions owing to their small size and the ease with which they can be evolved to recognize a large variety of different ligands. Here we present the solution structure of an RNA aptamer that binds triphenyl dyes in complex with malachite green and compare it with a previously determined crystal structure of a complex formed with tetramethylrosamine. The structures illustrate how the same RNA binding pocket can adapt to accommodate both planar and nonplanar ligands. Binding studies with single- and double-substitution mutant aptamers are used to correlate three-dimensional structure with complex stability. The two RNA-ligand complex structures allow a discussion of structural changes that have been observed in the ligand in the context of the overall complex structure. Base pairing and stacking interactions within the RNA fold the phosphate backbone into a structure that results in an asymmetric charge distribution within the binding pocket that forces the ligand to adapt through a redistribution of the positive partial charge.

Fluorescence-based methods for evaluating the RNA affinity and specificity of HIV-1 Rev-RRE inhibitors

RNA plays a pivotal role in the replication of all organisms, including viral and bacterial pathogens. The development of small molecules that selectively interfere with undesired RNA activity is a promising new direction for drug design. Currently, there are no anti-HIV treatments that target nucleic acids. This article presents the HIV-1 Rev response element (RRE) as an important focus for the development of antiviral agents that target RNA. The Rev binding site on the RRE is highly conserved, even between different groups of HIV-1 isolates. Compounds that inhibit HIV replication by binding to the RRE and displacing Rev are therefore expected to retain activity across groups of genetically diverse HIV infections. Systematic evaluations of both the RRE affinity and specificity of numerous small molecule inhibitors are essential for deciphering the parameters that govern effective RRE recognition. This article discusses fluorescence-based techniques that are useful for probing a small molecule's RRE affinity and its ability to inhibit Rev-RRE binding. Rev displacement experiments can be conducted by observing the fluorescence anisotropy of a fluorescein-labeled Rev peptide, or by quantifying its displacement from a solid-phase immobilized RRE. Experiments conducted in the presence of competing nucleic acids are useful for evaluating the RRE specificity of Rev-RRE inhibitors. The discovery and characterization of new RRE ligands are described. Eilatin is a polycyclic aromatic heterocycle that has at least one binding site on the RRE (apparent Kd is approximately 0.13 microM), but it does not displace Rev upon binding the RRE (IC50 > 3 microM). In contrast, ethidium bromide and two eilatin-containing metal complexes show better consistency between their RRE affinity and their ability to displace a fluorescent Rev peptide from the RRE. These results highlight the importance of conducting orthogonal binding assays that establish both the RNA affinity of a small molecule and its ability to inhibit the function of the RNA target. Some Rev-RRE inhibitors, including ethidium bromide, Lambda-[Ru(bpy)(2)eilatin]2+, and Delta-[Ru(bpy)(2)eilatin]2+ also inhibit HIV-1 gene expression in cell cultures (IC50 = 0.2-3 microM). These (and similar) results should facilitate the future discovery and implementation of anti-HIV drugs that are targeted to viral RNA sites. In addition, a deeper general understanding of RNA-small molecule recognition will assist in the effective targeting of other therapeutically important RNA sites.

RNA as a target for developing antivirals

The base of knowledge concerning RNA structure and function has been expanding rapidly in recent years. Simultaneously, an increasing awareness of the pivotal role RNA plays in viral diseases has prompted many researchers to apply new technologies in high-throughput screening and molecular modelling to the design of antiviral drugs that target RNA. While the two RNA viruses with the greatest unmet medical need, HIV and HCV, have been most actively pursued, the approaches discussed in this review are relevant to all virus infections. Both traditional small-molecule and large-molecule therapeutics, such as antisense, ribozymes and interfering dsRNAs have been described, and several molecules are under development for commercialization. The purpose of this review is to summarize the current state of the art in this field and to postulate new directions in the future.

Aminoglycoside-derived cationic lipids as efficient vectors for gene transfection in vitro and in vivo

BACKGROUND

Cationic lipids are at present very actively investigated for gene transfer studies and gene therapy applications. Basically, they rely on the formation of DNA/lipid aggregates via electrostatic interactions between their cationic headgroup and the negatively charged DNA. Although their structure/activity relationships are not well understood, it is generally agreed that the nature of the positive headgroup impacts on their transfection activity. Thus, we have directed our efforts toward the development of cationic lipids with novel cationic moieties. In the present work, we have explored the transfection potential of the lipophilic derivatives of the aminoglycoside kanamycin A. Indeed, aminoglycosides, which are natural polyamines known to bind to nucleic acids, provide a favorable scaffold for the synthesis of a variety of cationic lipids because of their structural features and multifunctional nature.

METHODS AND RESULTS

We report here the synthesis of a cationic cholesterol derivative characterized by a kanamycin A headgroup and of its polyguanidinylated derivative. The amino-sugar-based cationic lipid is highly efficient for gene transfection into a variety of mammalian cell lines when used either alone or as a liposomal formulation with the neutral phospholipid dioleoylphosphatidylethanolamine (DOPE). Its polyguanidinylated derivative was also found to mediate in vitro gene transfection. In addition, colloidally stable kanamycin-cholesterol/DOPE lipoplexes were found to be efficient for gene transfection into the mouse airways in vivo.

CONCLUSIONS

These results reveal the usefulness of cationic lipids characterized by headgroups composed of an aminoglycoside or its guanidinylated derivative for gene transfection in vitro and in vivo.

Structure-based computational database screening, in vitro assay, and NMR assessment of compounds that target TAR RNA

There has been little prior effort to discover new drugs on the basis of a unique RNA structure. Binding of the viral transactivator Tat to the 5' bulge of the transactivation response (TAR) element is necessary for HIV-1 replication, so TAR RNA is a superb target. A computational approach was developed to screen a large chemical library for binding to a three-dimensional RNA structure. Scoring function development, flexible ligand docking, and limited target flexibility were essential. From the ranked list of compounds predicted to bind TAR, 43 were assayed for inhibition of the Tat-TAR interaction via electrophoretic mobility shift assays. Eleven compounds (between 0.1 and 1 microM) inhibited the Tat-TAR interaction, and some inhibited Tat transactivation in cells. NMR spectra verified specific binding to the 5' bulge and no interaction with other regions of TAR.

Identification of ligands for RNA targets via structure-based virtual screening: HIV-1 TAR

Binding of the Tat protein to TAR RNA is necessary for viral replication of HIV-1. We screened the Available Chemicals Directory (ACD) to identify ligands to bind to a TAR RNA structure using a four-step docking procedure: rigid docking first, followed by three steps of flexible docking using a pseudobrownian Monte Carlo minimization in torsion angle space with progressively more detailed conformational sampling on a progressively smaller list of top-ranking compounds. To validate the procedure, we successfully docked ligands for five RNA complexes of known structure. For ranking ligands according to binding avidity, an empirical binding free energy function was developed which accounts, in particular, for solvation, isomerization free energy, and changes in conformational entropy. System-specific parameters for the function were derived on a training set of RNA/ligand complexes with known structure and affinity. To validate the free energy function, we screened the entire ACD for ligands for an RNA aptamer which binds L-arginine tightly. The native ligand ranked 17 out of ca. 153,000 compounds screened, i.e., the procedure is able to filter out >99.98% of the database and still retain the native ligand. Screening of the ACD for TAR ligands yielded a high rank for all known TAR ligands contained in the ACD and suggested several other potential TAR ligands. Eight of the highest ranking compounds not previously known to be ligands were assayed for inhibition of the Tat-TAR interaction, and two exhibited a CD50 of ca. 1 microM.

Simulation of the structure and dynamics of nonhelical RNA motifs

Computer simulation methods are increasingly being used to study possible conformations and dynamics of structural motifs in RNA. Recent results of molecular dynamics simulations and continum solvent studies of RNA structures and RNA-ligand complexes show promising agreement with experimental data. Combined with the ongoing progress in the experimental characterization of RNA structure and thermodynamics, these computational approaches can help to better understand the mechanism of RNA structure formation and the binding of ligands.

Aminoglycoside binding in the major groove of duplex RNA: the thermodynamic and electrostatic forces that govern recognition

We use a combination of spectroscopic, calorimetric, viscometric and computer modeling techniques to characterize the binding of the aminoglycoside antibiotic, tobramycin, to the polymeric RNA duplex, poly(rI).poly(rC), which exhibits the characteristic A-type conformation that is conserved among natural and synthetic double-helical RNA sequences. Our results reveal the following significant features: (i) CD-detected binding of tobramycin to poly(rI).poly(rC) reveals an apparent site size of four base-pairs per bound drug molecule; (ii) tobramycin binding enhances the thermal stability of the host poly(rI).poly(rC) duplex, the extent of which decreases upon increasing in Na(+) concentration and/or pH conditions; (iii) the enthalpy of tobramycin- poly(rI).poly(rC) complexation increases with increasing pH conditions, an observation consistent with binding-induced protonation of one or more drug amino groups; (iv) the affinity of tobramycin for poly(rI).poly(rC) is sensitive to both pH and Na(+) concentration, with increases in pH and/or Na(+) concentration resulting in a concomitant reduction in binding affinity. The salt dependence of the tobramycin binding affinity reveals that the drug binds to the host RNA duplex as trication. (v) The thermodynamic driving force for tobramycin- poly(rI).poly(rC) complexation depends on pH conditions. Specifically, at pH< or =6.0, tobramycin binding is entropy driven, but is enthalpy driven at pH > 6.0. (vi) Viscometric data reveal non-intercalative binding properties when tobramycin complexes with poly(rI).poly(rC), consistent with a major groove-directed mode of binding. These data also are consistent with a binding-induced reduction in the apparent molecular length of the host RNA duplex. (vii) Computer modeling studies reveal a tobramycin-poly(rI). poly(rC) complex in which the drug fits snugly at the base of the RNA major groove and is stabilized, at least in part, by an array of hydrogen bonding interactions with both base and backbone atoms of the host RNA. These studies also demonstrate an inability of tobramycin to form a stable low-energy complex with the minor groove of the poly(rI).poly(rC) duplex. In the aggregate, our results suggest that tobramycin-RNA recognition is dictated and controlled by a broad range of factors that include electrostatic interactions, hydrogen bonding interactions, drug protonation reactions, and binding-induced alterations in the structure of the host RNA. These modulatory effects on tobramycin-RNA complexation are discussed in terms of their potential importance for the selective recognition of specific RNA structural motifs, such as asymmetric internal loops or hairpin loop-stem junctions, by aminoglycoside antibiotics and their derivatives.

Structural basis for recognition of the RNA major groove in the tau exon 10 splicing regulatory element by aminoglycoside antibiotics

Drug-like molecules that bind RNA with sequence selectivity would provide valuable tools to elucidate gene expression pathways and new avenues to the treatment of degenerative and chronic conditions. Efforts at discovering such agents have been hampered, until recently, by the limited knowledge of RNA recognition principles. Several recent structures of aminoglycoside-RNA complexes have begun to reveal the structural basis for RNA-drug recognition. However, the absence of suitable chemical scaffolds known to bind the RNA major groove, where specificity could be provided by the diversity of functional groups exposed on the RNA bases, has represented a major obstacle. Here we report an investigation of the structural basis for recognition of an RNA stem-loop by neomycin, a naturally occurring aminoglycoside antibiotic. We found that neomycin binds the RNA stem-loop that regulates alternative splicing of exon 10 within the gene coding for human tau protein. Mutations within this splicing regulatory element destabilise the RNA structure and cause frontotemporal dementia and Parkinsonism linked to chromosome 17 (FTDP-17), an autosomal dominant condition leading to neurodegeneration and death. The three-dimensional structure of the RNA-neomycin complex shows interaction of the drug in the major groove of the short RNA duplex, where familial mutations cluster. Analysis of the structure shows how aminoglycosides and related drugs bind to the RNA major groove, adding to our understanding of the principles of drug-RNA recognition.

The H(2)(18)O solvent-induced isotope shift in (19)F NMR

The H(2)(16)O/H(2)(18)O solvent-induced isotope shifts ((18)O SIIS) of the (19)F NMR signals of a number of fluorine compounds have been measured. These isotope shifts are observed to be upfield, downfield, or zero, depending on the specific compound and the precise solution conditions. At 25 degrees C and with an (18)O enrichment of 86%, the (18)O SIIS of several fluorinated amino acids were in the range of 0.0014-0.0018 ppm downfield. 5-Fluorouridine displays a significantly wider range of (18)O SIIS values. A 5-fluorouridine-labeled 16-mer RNA also displayed observable (18)O SIIS values, but the characteristics of these were significantly modified from those of free 5-fluorouridine. The experimental observations are consistent with the (18)O SIIS being composed of upfield and downfield components, with the relative contributions of these determining the size and direction of the overall isotope shift. This is discussed in terms of a combination of van der Waals interactions between the fluorine atom and the solvent, electrical and hydrogen bonding effects, and the perturbations to these due to (18)O substitution in the solvent water. This isotope effect promises to be a highly useful tool in a range of (19)F NMR studies.

tRNA(Phe) binds aminoglycoside antibiotics

Aminoglycoside antibiotics have recently been found to bind to a variety of unrelated RNA molecules, including sequences that are important for retroviral replication. We report the binding of neomycin B, kanamycin A, and Neo-Neo (a synthetic neomycin-neomycin dimer) to tRNA(Phe). Using thermal denaturation studies, fluorescence spectroscopy, Pb2+-mediated tRNA(Phe) cleavage, and gel mobility shift assays, we have established that aminoglycosides interact with yeast tRNA(Phe) and are likely to induce a conformational change. Thermal denaturation studies revealed that aminoglycosides have a substantial stabilizing effect on tRNA(Phe) secondary and tertiary structures, much greater than the stabilization effect of spermine, an unstructured polyamine. Aminoglycoside-induced inhibition of Pb2+-mediated tRNA(Phe) cleavage yielded IC50 values of: 5 microM for Neo-Neo, 100 microM for neomycin B, > 1 mM for kanamycin A, and > 10 mM for spermine. Enzymatic and chemical footprinting indicate that the anticodon stem as well as the junction of the TpsiC and D loops are preferred aminoglycoside binding sites.

Role of sugar re-puckering in the transition of A and B forms of DNA in solution. A molecular dynamics study

State of the art molecular dynamic simulations show that simple modification of the sugar puckering of 2'deoxyriboses leads to a reversible change between two stable forms of DNA which resemble very closely the canonical A and B duplex forms. Analysis of the A-type and B-type structures reveals interesting, and previously unknown features of these two families of conformations of the DNA.

Molecular basis of HIV-1 TAR RNA specific recognition by an acridine tat-antagonist

We investigated the interaction of a highly potent acridine-based tat-antagonist with the TAR RNA of HIV-1. The wild type TAR RNA and three mutants with U-->C23, G x C-->C x G26-39 or G x C-->A x U26-39 substitutions were used as substrates to study the molecular basis of drug-TAR RNA complex formation. Melting temperature and RNase protection experiments reveal that the G x C26-39 pair is a critical element for specific major groove recognition of TAR at the pyrimidine bulge. The results provide a rational basis for future design of optimized tat/TAR inhibitors.

Translation of hepatitis C virus RNA

The 341-nucleotide 5' non-translated region is the most conserved part of the hepatitis C virus (HCV) genome. It contains a highly structured internal ribosomal entry site (IRES) that mediates cap-independent initiation of translation of the viral polyprotein by a mechanism that is unprecedented in eukaryotes. The first step in translation initiation is assembly of eukaryotic initiation factor (eIF) 3, eIF2, GTP, initiator tRNA and a 40S ribosomal subunit into a 43S preinitiation complex. The HCV IRES recruits this complex and directs its precise attachment at the initiation codon to form a 48S complex in a process that does not involve eIFs 4A, 4B or 4F. The IRES contains sites that bind independently with the eIF3 and 40S subunit components of 43S complexes, and structural determinants that ensure the correct spatial orientation of these binding sites so that the 48S complex assembles precisely at the initiation codon.

Templates for design of inhibitors for serine proteases: application of the program DOCK to the discovery of novel inhibitors for thrombin

The program DOCK was used to search for novel inhibitors for alpha-thrombin. Four among the top twelve best scoring compounds from the Cambridge Structural Data Base inhibited this enzyme, and three of them inhibited alpha-thrombin in a competitive mode. These molecules are expected to serve as general templates for structural elaboration in targeting diverse serine proteases for selective inhibition.

RNA as a drug target: chemical, modelling, and evolutionary tools

Dramatic technical progress in RNA synthesis and structure determination has allowed several difficulties inherent to the preparation, handling and structural analysis of RNA to be overcome, and this has led to a wealth of information about RNA structure and its relationship with biological function. It is now fully recognized that RNA molecules intervene at all stages of cell life, not only because of key sequence motifs but also because of intricate three-dimensional folds. This realization has promoted RNA to a potential therapeutic target. As in protein motifs recognizing nucleic acids, groups of the molecule interacting with RNA contribute to specific binding through defined hydrogen bonds and van der Waals docking, while other parts contribute to the driving force of binding via less specific electrostatic interactions accompanied by water and ion displacement.