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

Coarse grained models reveal essential contributions of topological constraints to the conformational free energy of RNA bulges

Recent studies have shown that simple stereochemical constraints encoded at the RNA secondary structure level significantly restrict the orientation of RNA helices across two-way junctions and yield physically reasonable distributions of RNA 3D conformations. Here we develop a new coarse-grain model, TOPRNA, that is optimized for exploring detailed aspects of these topological constraints in complex RNA systems. Unlike prior models, TOPRNA effectively treats RNAs as collections of semirigid helices linked by freely rotatable single strands, allowing us to isolate the effects of secondary structure connectivity and sterics on 3D structure. Simulations of bulge junctions show that TOPRNA captures new aspects of topological constraints, including variations arising from deviations in local A-form structure, translational displacements of the helices, and stereochemical constraints imposed by bulge-linker nucleotides. Notably, these aspects of topological constraints define free energy landscapes that coincide with the distribution of bulge conformations in the PDB. Our simulations also quantitatively reproduce NMR RDC measurements made on HIV-1 TAR at low salt concentrations, although not for different TAR mutants or at high salt concentrations. Our results confirm that topological constraints are an important determinant of bulge conformation and dynamics and demonstrate the utility of TOPRNA for studying the topological constraints of complex RNAs.

A method of determining RNA conformational ensembles using structure-based calculations of residual dipolar couplings

We describe a method of determining the conformational fluctuations of RNA based on the incorporation of nuclear magnetic resonance (NMR) residual dipolar couplings (RDCs) as replica-averaged structural restraints in molecular dynamics simulations. In this approach, the alignment tensor required to calculate the RDCs corresponding to a given conformation is estimated from its shape, and multiple replicas of the RNA molecule are simulated simultaneously to reproduce in silico the ensemble-averaging procedure performed in the NMR measurements. We provide initial evidence that with this approach it is possible to determine accurately structural ensembles representing the conformational fluctuations of RNA by applying the reference ensemble test to the trans-activation response element of the human immunodeficiency virus type 1.

NMR studies of nucleic acid dynamics

Nucleic acid structures have to satisfy two diametrically opposite requirements; on one hand they have to adopt well-defined 3D structures that can be specifically recognized by proteins; on the other hand, their structures must be sufficiently flexible to undergo very large conformational changes that are required during key biochemical processes, including replication, transcription, and translation. How do nucleic acids introduce flexibility into their 3D structure without losing biological specificity? Here, I describe the development and application of NMR spectroscopic techniques in my laboratory for characterizing the dynamic properties of nucleic acids that tightly integrate a broad set of NMR measurements, including residual dipolar couplings, spin relaxation, and relaxation dispersion with sample engineering and computational approaches. This approach allowed us to obtain fundamental new insights into directional flexibility in nucleic acids that enable their structures to change in a very specific functional manner.

Thermodynamic studies of a series of homologous HIV-1 TAR RNA ligands reveal that loose binders are stronger Tat competitors than tight ones

RNA is a major drug target, but the design of small molecules that modulate RNA function remains a great challenge. In this context, a series of structurally homologous 'polyamide amino acids' (PAA) was studied as HIV-1 trans-activating response (TAR) RNA ligands. An extensive thermodynamic study revealed the occurence of an enthalpy-entropy compensation phenomenon resulting in very close TAR affinities for all PAA. However, their binding modes and their ability to compete with the Tat fragment strongly differ according to their structure. Surprisingly, PAA that form loose complexes with TAR were shown to be stronger Tat competitors than those forming tight ones, and thermal denaturation studies demonstrated that loose complexes are more stable than tight ones. This could be correlated to the fact that loose and tight ligands induce distinct RNA conformational changes as revealed by circular dichroism experiments, although nuclear magnetic resonance (NMR) experiments showed that the TAR binding site is the same in all cases. Finally, some loose PAA also display promising inhibitory activities on HIV-infected cells. Altogether, these results lead to a better understanding of RNA interaction modes that could be very useful for devising new ligands of relevant RNA targets.

Riboswitch control of aminoglycoside antibiotic resistance

The majority of riboswitches are regulatory RNAs that regulate gene expression by binding small-molecule metabolites. Here we report the discovery of an aminoglycoside-binding riboswitch that is widely distributed among antibiotic-resistant bacterial pathogens. This riboswitch is present in the leader RNA of the resistance genes that encode the aminoglycoside acetyl transferase (AAC) and aminoglycoside adenyl transferase (AAD) enzymes that confer resistance to aminoglycoside antibiotics through modification of the drugs. We show that expression of the AAC and AAD resistance genes is regulated by aminoglycoside binding to a secondary structure in their 5' leader RNA. Reporter gene expression, direct measurements of drug RNA binding, chemical probing, and UV crosslinking combined with mutational analysis demonstrate that the leader RNA functions as an aminoglycoside-sensing riboswitch in which drug binding to the leader RNA leads to the induction of aminoglycosides antibiotic resistance.

Strategies to Block HIV Transcription: Focus on Small Molecule Tat Inhibitors

After entry into the target cell, the human immunodeficiency virus type I (HIV) integrates into the host genome and becomes a proviral eukaryotic transcriptional unit. Transcriptional regulation of provirus gene expression is critical for HIV replication. Basal transcription from the integrated HIV promoter is very low in the absence of the HIV transactivator of transcription (Tat) protein and is solely dependent on cellular transcription factors. The 5' terminal region (+1 to +59) of all HIV mRNAs forms an identical stem-bulge-loop structure called the Transactivation Responsive (TAR) element. Once Tat is made, it binds to TAR and drastically activates transcription from the HIV LTR promoter. Mutations in either the Tat protein or TAR sequence usually affect HIV replication, indicating a strong requirement for their conservation. The necessity of the Tat-mediated transactivation cascade for robust HIV replication renders Tat one of the most desirable targets for transcriptional therapy against HIV replication. Screening based on inhibition of the Tat-TAR interaction has identified a number of potential compounds, but none of them are currently used as therapeutics, partly because these agents are not easily delivered for an efficient therapy, emphasizing the need for small molecule compounds. Here we will give an overview of the different strategies used to inhibit HIV transcription and review the current repertoire of small molecular weight compounds that target HIV transcription.

The interactions and recognition of cyclic peptide mimetics of Tat with HIV-1 TAR RNA: a molecular dynamics simulation study

The interaction of HIV-1 trans-activator protein Tat with its cognate trans-activation response element (TAR) RNA is critical for viral transcription and replication. Therefore, it has long been considered as an attractive target for the development of antiviral compounds. Recently, the conformationally constrained cyclic peptide mimetics of Tat have been tested to be a promising family of lead peptides. Here, we focused on two representative cyclic peptides termed as L-22 and KP-Z-41, both of which exhibit excellent inhibitory potency against Tat and TAR interaction. By means of molecular dynamics simulations, we obtained a detailed picture of the interactions between them and HIV-1 TAR RNA. In results, it is found that the binding modes of the two cyclic peptides to TAR RNA are almost identical at or near the bulge regions, whereas the binding interfaces at the apical loop exhibit large conformational heterogeneity. In addition, it is revealed that electrostatic interaction energy contributes much more to KP-Z-41 complex formation than to L-22 complex, which is the main source of energy that results in a higher binding affinity of KP-Z-41 over-22 for TAR RNA. Furthermore, we identified a conserved motif RRK (Arg-Arg-Lys) that is shown to be essential for specific binding of this class of cyclic peptides to TAR RNA. This work can provide a useful insight into the design and modification of cyclic peptide inhibitors targeting the association of HIV-1 Tat and TAR RNA.

Conformation and dynamics of nucleotides in bulges and symmetric internal loops in duplex DNA studied by EPR and fluorescence spectroscopies

The dynamics and conformation of base bulges and internal loops in duplex DNA were studied using the bifunctional spectroscopic probe Ç, which becomes fluorescent (Ç(f)) upon reduction of the nitroxide functional group, along with EPR and fluorescence spectroscopies. A one-base bulge was in a conformational equilibrium between looped-out and stacked states, the former favored at higher temperature and the latter at lower temperature. Stacking of bulge bases was favored in two- and three-base bulges, independent of temperature, resulting in DNA bending as evidenced by increased fluorescence of Ç(f). EPR spectra of Ç-labeled three-, four- and five-base symmetrical interior DNA bulges at 20 °C showed low mobility, indicating that the spin-label was stacked within the loop. The spin-label mobility at 37 °C increased as the loops became larger. A considerable variation in fluorescence between different loops was observed, as well as a temperature-dependence within constructs. Fluorescence unexpectedly increased as the size of the loop decreased at 2 °C. Fluorescence of the smallest loops, where a single T·T mismatch was located between the stem region and the probe, was even larger than for the single strand, indicating a considerable local structural deformation of these loops from regular B-DNA. These results show the value of combining EPR and fluorescence spectroscopy to study non-helical regions of nucleic acids.

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.

Toward targeting RNA structure: branched peptides as cell-permeable ligands to TAR RNA

Rational design of RNA ligands continues to be a formidable challenge, but the potential powerful applications in biology and medicine catapults it to the forefront of chemical research. Indeed, small molecule and macromolecular intervention are attractive approaches, but selectivity and cell permeability can be a hurdle. An alternative strategy is to use molecules of intermediate molecular weight that possess large enough surface area to maximize interaction with the RNA structure but are small enough to be cell-permeable. Herein, we report the discovery of nontoxic and cell-permeable branched peptide (BP) ligands that bind to TAR RNA in the low micromolar range from on-bead high-throughput screening of 4,096 compounds. TAR is a short RNA motif in the 5'-UTR of HIV-1 that is responsible for efficient generation of full RNA transcripts. We demonstrate that BPs are selective for the native TAR RNA structure and that "branching" in peptides provides multivalent interaction, which increases binding affinity to RNA.

Interactions of protein side chains with RNA defined with REDOR solid state NMR

Formation of the complex between human immunodeficiency virus type-1 Tat protein and the transactivation response region (TAR) RNA is vital for transcriptional elongation, yet the structure of the Tat-TAR complex remains to be established. The NMR structures of free TAR, and TAR bound to Tat-derived peptides have been obtained by solution NMR, but only a small number of intermolecular NOEs could be identified unambiguously, preventing the determination of a complete structure. Here we show that a combination of multiple solid state NMR REDOR experiments can be used to obtain multiple distance constraints from (15)N to (13)C spins within the backbone and side chain guanidinium groups of arginine in a Tat-derived peptide, using (19)F spins incorporated into the base of U23 in TAR and (31)P spins in the P22 and P23 phosphate groups. Distances between the side chain of Arg52 and the base and phosphodiester backbone near U23 measured by REDOR NMR are comparable to distances observed in solution NMR-derived structural models, indicating that interactions of TAR RNA with key amino acid side chains in Tat are the same in the amorphous solid state as in solution. This method is generally applicable to other protein-RNA complexes where crystallization or solution NMR has failed to provide high resolution structural information.

A small-molecule probe induces a conformation in HIV TAR RNA capable of binding drug-like fragments

The HIV-1 transactivation response (TAR) element-Tat interaction is a potentially valuable target for treating HIV infection, but efforts to develop TAR-binding antiviral drugs have not yet yielded a successful candidate for clinical development. In this work, we describe a novel approach toward screening fragments against RNA that uses a chemical probe to target the Tat-binding region of TAR. This probe fulfills two critical roles in the screen: by locking the RNA into a conformation capable of binding other fragments, it simultaneously allows the identification of proximal binding fragments by ligand-based NMR. Using this approach, we have discovered six novel TAR-binding fragments, three of which were docked relative to the probe-RNA structure using experimental NMR restraints. The consistent orientations of functional groups in our data-driven docked structures and common electrostatic properties across all fragment leads reveal a surprising level of selectivity by our fragment-sized screening hits. These models further suggest linking strategies for the development of higher-affinity lead compounds for the inhibition of the TAR-Tat interaction.

6-desfluoroquinolones as HIV-1 Tat-mediated transcription inhibitors

The current anti-HIV treatments fail to completely eradicate the virus in HIV-infected individuals, mainly as a result of a small pool of latently infected cells. This issue, together with the emergence of multidrug-resistant viruses, clearly highlights the need to find additional strategies. An overview of the Tat-mediated transcription inhibitors 6-desfluoroquinolones (6-DFQs), identified by our group, is given in this review along with a critical appraisal of their advantages and drawbacks. Attempts are also made to place them within the context of new potential anti-HIV therapeutics. Due to their innovative mechanism of action, the 6-DFQs could be interesting candidates for use in association with the currently used cocktail of drugs. Their potential as antivirals deserves further investigation.

Inhibition of both HIV-1 reverse transcription and gene expression by a cyclic peptide that binds the Tat-transactivating response element (TAR) RNA

The RNA response element TAR plays a critical role in HIV replication by providing a binding site for the recruitment of the viral transactivator protein Tat. Using a structure-guided approach, we have developed a series of conformationally-constrained cyclic peptides that act as structural mimics of the Tat RNA binding region and block Tat-TAR interactions at nanomolar concentrations in vitro. Here we show that these compounds block Tat-dependent transcription in cell-free systems and in cell-based reporter assays. The compounds are also cell permeable, have low toxicity, and inhibit replication of diverse HIV-1 strains, including both CXCR4-tropic and CCR5-tropic primary HIV-1 isolates of the divergent subtypes A, B, C, D and CRF01_AE. In human peripheral blood mononuclear cells, the cyclic peptidomimetic L50 exhibited an IC₅₀ ∼250 nM. Surprisingly, inhibition of LTR-driven HIV-1 transcription could not account for the full antiviral activity. Timed drug-addition experiments revealed that L-50 has a bi-phasic inhibition curve with the first phase occurring after HIV-1 entry into the host cell and during the initiation of HIV-1 reverse transcription. The second phase coincides with inhibition of HIV-1 transcription. Reconstituted reverse transcription assays confirm that HIV-1 (−) strand strong stop DNA synthesis is blocked by L50-TAR RNA interactions in-vitro. These findings are consistent with genetic evidence that TAR plays critical roles both during reverse transcription and during HIV gene expression. Our results suggest that antiviral drugs targeting TAR RNA might be highly effective due to a dual inhibitory mechanism.

The HIV-1 transactivator protein (Tat), together with the elongation factor P-TEFb binds to an HIV-1 RNA secondary structure in the 5′-UTRs of nascent viral mRNAs (TAR) and promotes transcription elongation. This process has been an attractive target for drug development but previous inhibitors that bind either Tat or TAR have been plagued by poor inhibition of virus replication, limited cell penetration, and off-target effects. In this article, we describe a series of rationally designed cyclic peptides that block Tat-TAR interactions. L50, the most potent of these compounds, inhibits a wide range of HIV-1 strains from around the world. Remarkably, L50 inhibits two distinct steps in the HIV-1 lifecycle. As expected, L50 inhibits Tat-dependent HIV-1 transcription, but the majority of its anti-HIV activity is due to a block in reverse transcription, i.e. synthesis of the proviral DNA from the RNA genome. L50 inhibition of reverse transcription reveals an important role for TAR RNA during reverse transcription as well as providing one of first examples of a drug with a dual mechanism of action.

Dynamic ensemble view of the conformational landscape of HIV-1 TAR RNA and allosteric recognition

RNA conformational dynamics and the resulting structural heterogeneity play an important role in RNA functions, e.g., recognition. Recognition of HIV-1 TAR RNA has been proposed to occur via a conformational capture mechanism. Here, using ultrafast time-resolved fluorescence spectroscopy, we have probed the complexity of the conformational landscape of HIV-1 TAR RNA and monitored the position-dependent changes in the landscape upon binding of a Tat protein-derived peptide and neomycin B. In the ligand-free state, the TAR RNA samples multiple families of conformations with various degrees of base stacking around the three-nucleotide bulge region. Some subpopulations partially resemble those ligand-bound states, but the coaxially stacked state is below the detection limit. When Tat or neomycin B binds, the bulge region as an ensemble undergoes a conformational transition in a position-dependent manner. Tat and neomycin B induce mutually exclusive changes in the TAR RNA underlying the mechanism of allosteric inhibition at an ensemble level with residue-specific details. Time-resolved anisotropy decay measurements revealed picosecond motions of bases in both ligand-free and ligand-bound states. Mutation of a base pair at the bulge--stem junction has differential effects on the conformational distributions of the bulge bases. A dynamic model of the ensemble view of the conformational landscape for HIV-1 TAR RNA is proposed, and the implication of the general mechanism of RNA recognition and its impact on RNA-based therapeutics are discussed.

Therapeutic targeting of HCV internal ribosomal entry site RNA

HCV infection is a significant human disease, leading to liver cirrhosis and cancer, and killing >10,000 people in the US annually. Translation of the viral RNA genome is initiated by ribosomal binding to a highly structured RNA element, the internal ribosomal entry site (IRES), which presents a novel target for therapeutic intervention. We will first discuss studies of oligonucleotide therapeutics targeting various regions of the 340-nucleotide IRES, many of which have effectively blocked IRES function in vitro and are active against virus replication in cell culture. Although low nanomolar potencies have been obtained for DNA- and RNA-based molecules, stability and drug delivery challenges remain to be addressed for these particular HCV compounds. Several classes of small molecule inhibitors have been identified from screening protocols or designed from established RNA therapeutic scaffolds. In particular, small molecule IRES inhibitors based on a benzimidazole scaffold bind specifically to the IRES, and inhibit viral replication in cell culture at micromolar concentrations with low toxicity. The structure of the RNA target in complex with a representative member of these small molecule inhibitors demonstrates that a large RNA conformational change occurs upon inhibitor binding. The RNA complex shows how the inhibitor alters the global RNA structure and provides a framework for structure-based drug design of novel HCV therapeutics.

HIV-1 TAR RNA spontaneously undergoes relevant apo-to-holo conformational transitions in molecular dynamics and constrained geometrical simulations

We report all-atom molecular dynamics and replica exchange molecular dynamics simulations on the unbound human immunodeficiency virus type-1 (HIV-1) transactivation responsive region (TAR) RNA structure and three TAR RNA structures in bound conformations of, in total, approximately 250 ns length. We compare the extent of observed conformational sampling with that of the conceptually simpler and computationally much cheaper constrained geometrical simulation approach framework rigidity optimized dynamic algorithm (FRODA). Atomic fluctuations obtained by replica-exchange molecular dynamics (REMD) simulations agree quantitatively with those obtained by molecular dynamics (MD) and FRODA simulations for the unbound TAR structure. Regarding the stereochemical quality of the generated conformations, backbone torsion angles and puckering modes of the sugar-phosphate backbone were reproduced equally well by MD and REMD simulations, but further improvement is needed in the case of FRODA simulations. Essential dynamics analysis reveals that all three simulation approaches show a tendency to sample bound conformations when starting from the unbound TAR structure, with MD and REMD simulations being superior with respect to FRODA. These results are consistent with the experimental view that bound TAR RNA conformations are transiently sampled in the free ensemble, following a conformation selection model. The simulation-generated TAR RNA conformations have been successfully used as receptor structures for docking. This finding has important implications for RNA-ligand docking in that docking into an ensemble of simulation-generated RNA structures is shown to be a valuable means to cope with large apo-to-holo conformational transitions of the receptor structure.

Essential structural requirements for specific recognition of HIV TAR RNA by peptide mimetics of Tat protein

The pharmacological disruption of the interaction between the HIV Tat protein and its cognate transactivation response RNA (TAR) would generate novel anti-viral drugs with a low susceptibility to drug resistance, but efforts to discover ligands with sufficient potency to warrant pharmaceutical development have been unsuccessful. We have previously described a family of structurally constrained β-hairpin peptides that potently inhibits viral growth in HIV-infected cells. The nuclear magnetic resonance (NMR) structure of an inhibitory complex revealed that the peptide makes intimate contacts with the 3-nt bulge and the upper helix of the RNA hairpin, but that a single residue contacts the apical loop where recruitment of the essential cellular co-factor cyclin T₁ occurs. Attempting to extend the peptide to form more interactions with the RNA loop, we examined a library of longer peptides and achieved >6-fold improvement in affinity. The structure of TAR bound to one of the extended peptides reveals that the peptide slides down the major groove of the RNA, relative to our design, in order to maintain critical interactions with TAR. These conserved contacts involve three amino acid side chains and identify critical interaction points required for potent and specific binding to TAR RNA. They constitute a template of essential interactions required for inhibition of this RNA.

Base pairs and pseudo pairs observed in RNA-ligand complexes

Previously, a geometric nomenclature was proposed in which RNA base pairs were classified by their interaction edges (Watson-Crick, Hoogsteen or sugar-edge) and the glycosidic bond orientations relative to the hydrogen bonds formed (cis or trans). Here, base pairs and pseudo pairs observed in RNA-ligand complexes are classified in a similar manner. Twenty-one basic geometric families are geometrically possible (18 for base pairs formed between a nucleic acid base and a ligand containing heterocycle and 3 families for pseudo pairs). Of those, 16 of them have been observed in X-ray and/or NMR structures.

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.

Domain-elongation NMR spectroscopy yields new insights into RNA dynamics and adaptive recognition

By simplifying the interpretation of nuclear magnetic resonance spin relaxation and residual dipolar couplings data, recent developments involving the elongation of RNA helices are providing new atomic insights into the dynamical properties that allow RNA structures to change functionally and adaptively. Domain elongation, in concert with spin relaxation measurements, has allowed the detailed characterization of a hierarchical network of local and collective motional modes occurring at nanosecond timescale that mirror the structural rearrangements that take place following adaptive recognition. The combination of domain elongation with residual dipolar coupling measurements has allowed the experimental three-dimensional visualization of very large amplitude rigid-body helix motions in HIV-1 transactivation response element (TAR) that trace out a highly choreographed trajectory in which the helices twist and bend in a correlated manner. The dynamic trajectory allows unbound TAR to sample many of its ligand bound conformations, indicating that adaptive recognition occurs by "conformational selection" rather than "induced fit." These studies suggest that intrinsic flexibility plays essential roles directing RNA conformational changes along specific pathways.

Multivalent binding oligomers inhibit HIV Tat-TAR interaction critical for viral replication

We describe the development of a new type of scaffold to target RNA structures. Multivalent binding oligomers (MBOs) are molecules in which multiple sidechains extend from a polyamine backbone such that favorable RNA binding occurs. We have used this strategy to develop MBO-based inhibitors to prevent the association of a protein-RNA complex, Tat-TAR, that is essential for HIV replication. In vitro binding assays combined with model cell-based assays demonstrate that the optimal MBOs inhibit Tat-TAR binding at low micromolar concentrations. Antiviral studies are also consistent with the in vitro and cell-based assays. MBOs provide a framework for the development of future RNA-targeting molecules.

Structures of HIV TAR RNA-ligand complexes reveal higher binding stoichiometries

Target TAR by NMR: Tripeptides containing arginines as terminal residues and non-natural amino acids as central residues are good leads for drug design to target the HIV trans-activation response element (TAR). The structural characterization of the RNA-ligand complex by NMR spectroscopy reveals two specific binding sites that are located at bulge residue U23 and around the pyrimidine-stretch U40-C41-U42 directly adjacent to the bulge.

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.

Simultaneous recognition of HIV-1 TAR RNA bulge and loop sequences by cyclic peptide mimics of Tat protein

The interaction of the HIV-1 transactivator protein Tat with its transactivation response (TAR) RNA is an essential step in viral replication and therefore an attractive target for developing antivirals with new mechanisms of action. Numerous compounds that bind to the 3-nt bulge responsible for binding Tat have been identified in the past, but none of these molecules had sufficient potency to warrant pharmaceutical development. We have discovered conformationally-constrained cyclic peptide mimetics of Tat that are specific nM inhibitors of the Tat-TAR interaction by using a structure-based approach. The lead peptides are nearly as active as the antiviral drug nevirapine against a variety of clinical isolates in human lymphocytes. The NMR structure of a peptide-RNA complex reveals that these molecules interfere with the recruitment to TAR of both Tat and the essential cellular cofactor transcription elongation factor-b (P-TEFb) by binding simultaneously at the RNA bulge and apical loop, forming an unusually deep pocket. This structure illustrates additional principles in RNA recognition: RNA-binding molecules can achieve specificity by interacting simultaneously with multiple secondary structure elements and by inducing the formation of deep binding pockets in their targets. It also provides insight into the P-TEFb binding site and a rational basis for optimizing the promising antiviral activity observed for these cyclic peptides.

An approach to the construction of tailor-made amphiphilic peptides that strongly and selectively bind to hairpin RNA targets

The hairpin RNA motif is one of the most frequently observed secondary structures and is often targeted by therapeutic agents. An amphiphilic peptide with seven lysine and eight leucine residues and its derivatives were designed for use as ligands against RNA hairpin motifs. We hypothesized that variations in both the hydrophobic leucine-rich and hydrophilic lysine-rich spheres of these amphiphilic peptides would create extra attractive interactions with hairpin RNA targets. A series of alanine-scanned peptides were probed to identify the most influential lysine residues in the hydrophilic sphere. The binding affinities of these modified peptides with several hairpins, such as RRE, TAR from HIV, a short hairpin from IRES of HCV, and a hairpin from the 16S A-site stem from rRNA, were determined. Since the hairpin from IRES of HCV was the most susceptible to the initial series of alanine-scanned peptides, studies investigating how further variations in the peptides effect binding employed the IRES hairpin. Next, the important Lys residues were substituted by shorter chain amines, such as ornithine, to place the peptide deeper into the hairpin groove. In a few cases, a 70-fold improved binding was observed for peptides that contained the specifically located shorter amine side chains. To further explore changes in binding affinities brought about by alterations in the hydrophobic sphere, tryptophan residues were introduced in place of leucine. A few peptides with tryptophan in specific positions also displayed 70-fold improved binding affinities. Finally, double mutant peptides incorporating both specifically located shorter amine side chains in the hydrophilic region and tryptophan residues in the hydrophobic region were synthesized. The binding affinities of peptides containing the simple double modification were observed to be 80 times lower, and their binding specificities were increased 40-fold. The results of this effort provide important information about strategies that can be used to prepare peptides that both strongly and selectively target hairpin RNAs. Specifically, the findings indicate that tailor-made amphiphilic peptide ligands against certain hairpin RNAs can be obtained if the RNA target possesses a deep groove in which both the hydrophobic and hydrophilic spheres of the peptide interact.

Small molecule microarrays of RNA-focused peptoids help identify inhibitors of a pathogenic group I intron

Peptoids that inhibit the group I intron RNA from Candida albicans, an opportunistic pathogen that kills immunocompromised hosts, have been identified using microarrays. The arrayed peptoid library was constructed using submonomers with moieties similar to ones found in small molecules known to bind RNA. Library members that passed quality control analysis were spotted onto a microarray and screened for binding to the C. albicans group I intron ribozyme. Each ligand binder identified from microarray-based screening inhibited self-splicing in the presence of 1 mM nucleotide concentration of bulk yeast tRNA with IC(50)'s between 150 and 2200 microM. The binding signals and the corresponding IC(50)'s were used to identify features in the peptoids that predispose them for RNA binding. After statistical analysis of the peptoids' structures that bind, a second generation of inhibitors was constructed using these important features; all second generation inhibitors have improved potencies with IC(50)'s of <100 microM. The most potent inhibitor is composed of one phenylguanidine and three tryptamine submonomers and has an IC(50) of 31 microM. This compound is 6-fold more potent than pentamidine, a clinically used drug that inhibits self-splicing. These results show that (i) modulators of RNA function can be identified by designing RNA-focused chemical libraries and screening them via microarray; (ii) statistical analysis of ligand binders can identify features in leads that predispose them for binding to their targets; and (iii) features can then be programmed into second generation inhibitors to design ligands with improved potencies.

Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.

How binding of small molecule and peptide ligands to HIV-1 TAR alters the RNA motional landscape

The HIV-1 TAR RNA represents a well-known paradigm to study the role of dynamics and conformational change in RNA function. This regulatory RNA changes conformation in response to binding of Tat protein and of a variety of peptidic and small molecule ligands, indicating that its conformational flexibility and intrinsic dynamics play important roles in molecular recognition. We have used ¹³C NMR relaxation experiments to examine changes in the motional landscape of HIV-1 TAR in the presence of three ligands of different affinity and specificity. The ligands are argininamide, a linear peptide mimic of the Tat basic domain and a cyclic peptide that potently inhibits Tat-dependent activation of transcription. All three molecules induce the same motional characteristics within the three nucleotides bulge that represents the Tat-binding site. However, the cyclic peptide has a unique motional signature in the apical loop, which represents a binding site for the essential host co-factor cyclin T1. These results suggest that all peptidic mimics of Tat induce the same dynamics in TAR within this protein binding site. However, the new cyclic peptide mimic of Tat represents a new class of ligands with a unique effect on the dynamics and the structure of the apical loop.

Characterizing complex dynamics in the transactivation response element apical loop and motional correlations with the bulge by NMR, molecular dynamics, and mutagenesis

The HIV-1 transactivation response element (TAR) RNA binds a variety of proteins and is a target for developing anti-HIV therapies. TAR has two primary binding sites: a UCU bulge and a CUGGGA apical loop. We used NMR residual dipolar couplings, carbon spin relaxation (R(1) and R(2)), and relaxation dispersion (R(1rho)) in conjunction with molecular dynamics and mutagenesis to characterize the dynamics of the TAR apical loop and investigate previously proposed long-range interactions with the distant bulge. Replacement of the wild-type apical loop with a UUCG loop did not significantly affect the structural dynamics at the bulge, indicating that the apical loop and the bulge act largely as independent dynamical recognition centers. The apical loop undergoes complex dynamics at multiple timescales that are likely important for adaptive recognition: U31 and G33 undergo limited motions, G32 is highly flexible at picosecond-nanosecond timescales, and G34 and C30 form a dynamic Watson-Crick basepair in which G34 and A35 undergo a slow (approximately 30 mus) likely concerted looping in and out motion, with A35 also undergoing large amplitude motions at picosecond-nanosecond timescales. Our study highlights the power of combining NMR, molecular dynamics, and mutagenesis in characterizing RNA dynamics.

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.

Tripeptides from synthetic amino acids block the Tat-TAR association and slow down HIV spread in cell cultures

Non-natural amino acids with aromatic or heteroaromatic side chains were incorporated into tripeptides of the general structure Arg-X-Arg and tested as ligands of the HIV RNA element TAR. Some of these compounds could compete efficiently with the association of TAR and Tat and downregulated a TAR-controlled reporter gene in HeLa cells. Peptide 7, which contains a 2-pyrimidinyl-alkyl chain, also inhibited the spread of HIV-1 in cell cultures. NMR studies of 7 bound to HIV-2-TAR gave evidence for contacts in the bulge region.

Visualizing spatially correlated dynamics that directs RNA conformational transitions

RNAs fold into three-dimensional (3D) structures that subsequently undergo large, functionally important, conformational transitions in response to a variety of cellular signals. RNA structures are believed to encode spatially tuned flexibility that can direct transitions along specific conformational pathways. However, this hypothesis has proved difficult to examine directly because atomic movements in complex biomolecules cannot be visualized in 3D by using current experimental methods. Here we report the successful implementation of a strategy using NMR that has allowed us to visualize, with complete 3D rotational sensitivity, the dynamics between two RNA helices that are linked by a functionally important trinucleotide bulge over timescales extending up to milliseconds. The key to our approach is to anchor NMR frames of reference onto each helix and thereby directly measure their dynamics, one relative to the other, using 'relativistic' sets of residual dipolar couplings (RDCs). Using this approach, we uncovered super-large amplitude helix motions that trace out a surprisingly structured and spatially correlated 3D dynamic trajectory. The two helices twist around their individual axes by approximately 53 degrees and 110 degrees in a highly correlated manner (R = 0.97) while simultaneously (R = 0.81-0.92) bending by about 94 degrees. Remarkably, the 3D dynamic trajectory is dotted at various positions by seven distinct ligand-bound conformations of the RNA. Thus even partly unstructured RNAs can undergo structured dynamics that directs ligand-induced transitions along specific predefined conformational pathways.

NMR studies of RNA dynamics and structural plasticity using NMR residual dipolar couplings

An increasing number of RNAs are being discovered that perform their functions by undergoing large changes in conformation in response to a variety of cellular signals, including recognition of proteins and small molecular targets, changes in temperature, and RNA synthesis itself. The measurement of NMR residual dipolar couplings (RDCs) in partially aligned systems is providing new insights into the structural plasticity of RNA through combined characterization of large-amplitude collective helix motions and local flexibility in noncanonical regions over a wide window of biologically relevant timescales (<milliseconds). Here, we review RDC methodology for studying RNA structural dynamics and survey what has been learnt thus far from application of these methods. Future methodological challenges are also identified.

Characterizing the relative orientation and dynamics of RNA A-form helices using NMR residual dipolar couplings

We present a protocol for determining the relative orientation and dynamics of A-form helices in 13C/15N isotopically enriched RNA samples using NMR residual dipolar couplings (RDCs). Non-terminal Watson-Crick base pairs in helical stems are experimentally identified using NOE and trans-hydrogen bond connectivity and modeled using the idealized A-form helix geometry. RDCs measured in the partially aligned RNA are used to compute order tensors describing average alignment of each helix relative to the applied magnetic field. The order tensors are translated into Euler angles defining the average relative orientation of helices and order parameters describing the amplitude and asymmetry of interhelix motions. The protocol does not require complete resonance assignments and therefore can be implemented rapidly to RNAs much larger than those for which complete high-resolution NMR structure determination is feasible. The protocol is particularly valuable for exploring adaptive changes in RNA conformation that occur in response to biologically relevant signals. Following resonance assignments, the procedure is expected to take no more than 2 weeks of acquisition and data analysis time.

Cyclin box structure of the P-TEFb subunit cyclin T1 derived from a fusion complex with EIAV tat

The positive transcription elongation factor b (P-TEFb) is an essential regulator of viral gene expression during the life cycle of human immunodeficiency virus type 1 (HIV-1). Its cyclin T1 subunit forms a ternary complex with the viral transcriptional transactivator (Tat) protein and the transactivation response (TAR) RNA element thereby activating cyclin dependent kinase 9 (Cdk9), which stimulates transcription at the level of chain elongation. We report the structure of the cyclin box domain of human cyclin T1 at a resolution of 2.67 A. The structure was obtained by crystallographic analysis of a fusion protein composed of cyclin T1 linked to the transactivator protein Tat from equine infectious anemia virus (EIAV), which is functionally and structurally related to HIV-1 Tat. The conserved cyclin box domain of cyclin T1 exhibits structural features for interaction with physiological binding partners such as Cdk9. A recognition site for Cdk/Cyclin substrates is partly covered by a cyclin T-specific insert, suggesting specific interactions with regulatory factors. The previously identified Tat/TAR recognition motif (TRM) forms a C-terminal helix that is partly occluded in the cyclin box repeat interface, while cysteine 261 is accessible to form an intermolecular zinc finger with Tat. Residues of the TRM contribute to a positively charged groove that may directly attract RNA molecules. The EIAV Tat protein instead appeared undefined from the electron density map suggesting that it is highly disordered. Functional experiments confirmed the TAR binding properties of the fusion protein and suggested residues on the second cyclin box repeat to contribute to Tat stimulated transcription.

Probing Na(+)-induced changes in the HIV-1 TAR conformational dynamics using NMR residual dipolar couplings: new insights into the role of counterions and electrostatic interactions in adaptive recognition

Many regulatory RNAs undergo large changes in structure upon recognition of proteins and ligands, but the mechanism by which this occurs remains poorly understood. Using NMR residual dipolar coupling (RDCs), we characterized Na+-induced changes in the structure and dynamics of the bulge-containing HIV-1 transactivation response element (TAR) RNA that mirrors changes induced by small molecules bearing a different number of cationic groups. Increasing the Na+ concentration from 25 to 320 mM led to a continuous reduction in the average inter-helical bend angle (from 46 degrees to 22 degrees ), inter-helical twist angle (from 66 degrees to -18 degrees ), and inter-helix flexibility (as measured by an increase in the internal generalized degree of order from 0.56 to 0.74). Similar conformational changes were observed with Mg2+, indicating that nonspecific electrostatic interactions drive the conformational transition, although results also suggest that Na+ and Mg2+ may associate with TAR in distinct modes. The transition can be rationalized on the basis of a population-weighted average of two ensembles comprising an electrostatically relaxed bent and flexible TAR conformation that is weakly associated with counterions and a globally rigid coaxial conformation that has stronger electrostatic potential and association with counterions. The TAR inter-helical orientations that are stabilized by small molecules fall around the metal-induced conformational pathway, indicating that counterions may help predispose the TAR conformation for target recognition. Our results underscore the intricate sensitivity of RNA conformational dynamics to environmental conditions and demonstrate the ability to detect subtle conformational changes using NMR RDCs.

Influence of generation 2-5 of PAMAM dendrimer on the inhibition of Tat peptide/ TAR RNA binding in HIV-1 transcription

The special binding of Tat protein to TAR RNA leads to the transcription of HIV-1 virus. In this study, the influence of 2-5 generation of PAMAM dendrimers on the inhibition of Tat protein/TAR RNA binding has been investigated. The absorption of PAMAM dendrimers on TAR RNA, fixed on a gold substrate through an avidin-biotion connection, was carried out by using a quartz crystal microbalance. Experimental result shows a Langmuir-type isotherm could be used to describe this kind of binding, implying a specific and monolayer adsorption existed. The combination coefficient (K(D)(-1))s can be calculated according to Langmuir Equation, having the order of G3 > G4 > G5 >Tat > G2, indicating that PAMAM G3, G4 and G5 having the possibility to be the inhibitors of HIV-1 transcription. The migration time (T(migra)) of capillary electrophoretic technique has the same sequence as (K(D)(-1))s. These two parameters could be used as simple and quantitative criteria for the selection of possible drugs from numerous candidates for HIV therapy in vitro.

Novel in vivo model for the study of human immunodeficiency virus type 1 transcription inhibitors: evaluation of new 6-desfluoroquinolone derivatives

Two novel 6-desfluoroquinolone derivatives, HM-12 and HM-13, were evaluated for anti-human immunodeficiency virus (anti-HIV) activity in acutely, chronically, and latently HIV type 1 (HIV-1)-infected cell cultures and were found to behave as potent HIV-1 transcription inhibitors. In order to extend this result in vivo, we developed an artificial hu-SCID mouse model for HIV-1 latency based on SCID mice engrafted with latently HIV-1-infected promyelocytic OM-10.1 cells in which HIV-1 can be reactivated in vivo by the administration of human tumor necrosis factor alpha (hTNF-alpha). Treating these SCID mice with HM-12 or HM-13 prior to hTNF-alpha stimulation resulted in a pronounced suppressive effect on viral reactivation. Since both quinolone derivatives were able to inhibit the reactivation of HIV-1 from this artificial viral reservoir in vivo, we provide encouraging evidence for the use of quinolones in the control of HIV-1 infections.

Synthesis and testing of a focused phenothiazine library for binding to HIV-1 TAR RNA

We have synthesized a series of phenothiazine derivatives, which were used to test the structure-activity relationship of binding to HIV-1 TAR RNA. Variations from our initial compound, 2-acetylphenothiazine, focused on two moieties: ring substitutions and n-alkyl substitutions. Binding characteristics were ascertained via NMR, principally by saturation transfer difference spectra of the ligand and imino proton resonance shifts of the RNA. Both ring and alkyl substitutions manifested NMR changes upon binding. In general, the active site, while somewhat flexible, has regions that can be capitalized for increased binding through van der Waals interactions and others that can be optimized for solubility in subsequent stages of development. However, binding can be nontrivially enhanced several-fold through optimization of van der Waals and hydrophilic sites of the scaffold.

Cross-interaction between JC virus agnoprotein and human immunodeficiency virus type 1 (HIV-1) Tat modulates transcription of the HIV-1 long terminal repeat in glial cells

The human polyomavirus JC virus (JCV) is the causative agent of the fatal demyelinating disease progressive multifocal leukoencephalopathy (PML), which is commonly seen in AIDS patients. The bicistronic viral RNA, which is transcribed at the late phase of infection, is responsible for expressing the viral capsid proteins and a small regulatory protein, agnoprotein. Immunohistochemical analysis of brain tissue from subjects with AIDS/PML revealed colocalization of the human immunodeficiency virus type 1 (HIV-1) transactivator, Tat, and JCV agnoprotein in nucleus and cytoplasm of "bizarre" astrocytes. In accord with this observation, we detected the copresence of agnoprotein and Tat in human astrocytes upon infection with JCV and HIV-1 or in astrocytic cells expressing these proteins after transfection. Interestingly, results from infection of human astrocytes with HIV-1 and JCV showed a decrease in the level of HIV-1 replication in cells that are coinfected with JCV. Conversely, a slight increase in the level of JCV replication was observed in the presence of HIV-1. The copresence of JCV and HIV-1 in astrocytes prompted us to investigate the possible cross-interaction of agnoprotein with Tat and its impact on HIV-1 gene transcription. Our results demonstrate that agnoprotein through its N-terminal domain associates with Tat and the interaction causes the suppression of Tat-mediated enhancement of HIV-1 promoter activity in these cells. Results from RNA and protein binding assays showed that agnoprotein can inhibit the association of Tat with its target RNA sequence, TAR, and with cyclin T1. Furthermore, agnoprotein is able to interfere with cross-interaction of Tat with the p65 subunit of NF-kappaB and Sp1, whose functions are critical for Tat activation of the long terminal repeat. These observations unravel a new pathway for the molecular interaction of these two viruses in biologically relevant cells in the brains of AIDS/PML patients.