Identification of amino acids that promote specific and rigid TAR RNA-tat protein complex formation

The Tat protein and the transactivation responsive (TAR) RNA form an essential complex in the HIV lifecycle, and mutations in the basic region of the Tat protein alter this RNA-protein molecular recognition. Here, EPR spectroscopy was used to identify amino acids, flanking an essential arginine of the Tat protein, which contribute to specific and rigid TAR-Tat complex formation by monitoring changes in the mobility of nitroxide spin-labeled TAR RNA nucleotides upon binding. Arginine to lysine N-terminal mutations did not affect TAR RNA interfacial dynamics. In contrast, C-terminal point mutations, R56 in particular, affected the mobility of nucleotides U23 and U38, which are involved in a base-triple interaction in the complex. This report highlights the role of dynamics in specific molecular complex formation and demonstrates the ability of EPR spectroscopy to study interfacial dynamics of macromolecular complexes.     

Mapping the interaction forces between TAR RNA and TAT peptides on GaAs surfaces using chemical force microscopy

The complexation of the HIV transactivation response element (TAR) RNA with the viral regulatory protein TAT is of enormous interest for the design of new sensing and therapeutic strategies. In this work, we anchored TAT peptides on GaAs surfaces using microcontact printing. Atomic force microscopy was used to quantify the interaction between TAR RNA and model TAT peptide sequences. Different pH conditions were utilized in order to assess specific vs nonspecific interactions. AFM tips functionalized with TAR RNA molecules were used to collect adhesion maps that displayed stronger interaction with peptide sequences that contained a greater number of arginine residues. All of the studies consistently showed a pH dependence of the interaction between the surface bound peptides and the TAR RNA on the AFM tips. This work quantifies the TAR RNA/TAT peptide interaction after one of the molecules is anchored on a surface. The conclusions in this paper are consistent with previous work and demonstrate that cationic residues are responsible for the polyelectrolyte-like affinity of TAT peptides for TAR RNA.     

Structural mimicry of retroviral tat proteins by constrained beta-hairpin peptidomimetics: ligands with high affinity and selectivity for viral TAR RNA regulatory elements

An approach is described to the design of beta-hairpin peptidomimetic ligands for bovine immunodeficiency virus (BIV) Tat protein, which inhibit binding to its transactivator response element (TAR) RNA. A library of peptidomimetics was derived by grafting onto a hairpin-inducing d-Pro-l-Pro template sequences related to the RNA recognition element in Tat. One hairpin mimetic was identified that binds tightly (K(d) approximately 150 nM) to BIV TAR, and another that binds also to HIV-1 TAR RNA (K(d) approximately 1-2 microM). (In the same assay, the wild-type BIV Tat(65-81) peptide binds to BIV TAR with K(d) approximately 50 nM.) The high-affinity BIV-Tat mimetic was shown to adopt a stable beta-hairpin conformation in free solution by NMR methods. Amino acid substitutions in this mimetic were shown to impact on the hairpin structure and to disrupt binding to the RNA. This family of conformationally constrained peptidomimetics affords insights into the structural requirements for binding to TAR RNA and provides a basis for the design of new ligands with increased inhibitory activity and specificity to both BIV and HIV TAR RNAs.     

13C/15N-19F intermolecular REDOR NMR study of the interaction of TAR RNA with Tat peptides

The complex of the HIV TAR RNA with the viral regulatory protein Tat is of considerable interest, but the plasticity of this interaction has made it impossible so far to establish the structure of that complex. In order to explore a new approach to obtain structural information on protein-RNA complexes, we performed (13)C/(15)N-(19)F REDOR NMR experiments in the solid state on TAR bound to a peptide comprising the RNA-binding section of Tat. A critical arginine in the peptide was uniformly (13)C and (15)N labeled, and 5-fluorouridine was incorporated at the U23 position of TAR. REDOR irradiation resulted in dephasing of the (13)C and (15)N resonances, indicating the proximity of the U23(5F)-C and U23(5F)-N spin pairs. Best fits to the REDOR data show the U23(5F)-C distances and the U23(5F)-N distances are in good agreement with the distances obtained from solution NMR structures of partial complexes of Tat with TAR. These results demonstrate that it is possible to study protein-RNA complexes using solid-state REDOR NMR measurements, adding to a growing list of solid state techniques for studying protein-nucleic acid complexes.     

Strong correlation between SHAPE chemistry and the generalized NMR order parameter (S2) in RNA

The functions of most RNA molecules are critically dependent on the distinct local dynamics that characterize secondary structure and tertiary interactions and on structural changes that occur upon binding by proteins and small molecule ligands. Measurements of RNA dynamics at nucleotide resolution set the foundation for understanding the roles of individual residues in folding, catalysis, and ligand recognition. In favorable cases, local order in small RNAs can be quantitatively analyzed by NMR in terms of a generalized order parameter, S2. Alternatively, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemistry measures local nucleotide flexibility in RNAs of any size using structure-sensitive reagents that acylate the 2'-hydroxyl position. In this work, we compare per-residue RNA dynamics, analyzed by both S2 and SHAPE, for three RNAs: the HIV-1 TAR element, the U1A protein binding site, and the Tetrahymena telomerase stem loop 4. We find a very strong correlation between the two measurements: nucleotides with high SHAPE reactivities consistently have low S2 values. We conclude that SHAPE chemistry quantitatively reports local nucleotide dynamics and can be used with confidence to analyze dynamics in large RNAs, RNA-protein complexes, and RNAs in vivo.     

A distance ruler for RNA using EPR and site-directed spin labeling

As a basic model study for measuring distances in RNA molecules using continuous wave (CW) EPR spectroscopy, site-directed spin-labeled 10-mer RNA duplexes and HIV-1 TAR RNA motifs with various interspin distances were examined. The spin labels were attached to the 2'-NH2 positions of appropriately placed uridines in the duplexes, and interspin distances were measured from both molecular dynamics simulations (MD) and Fourier deconvolution methods (FD). The 10-mer duplexes have interspin distances ranging from 10 A to 30 A based on MD; however, dipolar line broadening of the CW EPR spectrum is only observed for the RNAs for predicted interspin distances of 10-21 A and not for distances over 25 A. The conformational changes in TAR (transactivating responsive region) RNA in the presence and in the absence of different divalent metal ions were monitored by measuring distances between two nucleotides in the bulge region. The predicted interspin distances obtained from the FD method and those from MD calculations match well for both the model RNA duplexes and the structural changes predicted for TAR RNA. These results demonstrate that distance measurement using EPR spectroscopy is a potentially powerful method to help predict the structures of RNA molecules.     

Studying aminoglycoside antibiotic binding to HIV-1 TAR RNA by electrospray ionization mass spectrometry

The recognition of the aminoglycosides neomycin and streptomycin by HIV-1 TAR RNA was studied by electrospray ionization mass spectrometry (ESI-MS). Members of the aminoglycoside family of antibiotics are known to target a wide variety of RNA molecules. Neomycin and streptomycin inhibit the formation of the Tat protein–TAR RNA complex, an assembly that is believed to be necessary for HIV replication. The noncovalent complexes formed by the binding of aminoglycosides to TAR RNA and the Tat–TAR complex were detected by ESI-MS. Neomycin has a maximum binding stoichiometry of three and two to TAR RNA and to the Tat–TAR complex, respectively. Data from the ESI-MS experiments suggest that a high affinity binding site of neomycin is located near the three-nucleotide bulge region of TAR RNA. This is consistent with previous solution phase footprinting measurements [H.-Y. Mei et al., Biochemistry 37 (1998) 14204]. Neomycin has a higher affinity toward TAR RNA than streptomycin, as measured by ESI-MS competition binding experiments. A noncovalent complex formed between a small molecule inhibitor of TAR RNA, which has a similar solution binding affinity as the aminoglycosides, and TAR RNA is much less stable than the RNA–aminoglycoside complexes to collisional dissociation in the gas phase. It is believed that the small molecule inhibitor interacts with TAR RNA via hydrophobic interactions, whereas the aminoglycosides bind to RNAs through electrostatic forces. This difference in gas phase stabilities may prove useful for discerning the types of noncovalent forces holding complexes together.     

Driving in vitro selection of anti-HIV-1 TAR aptamers by magnesium concentration and temperature

In vitro selection with either DNA or RNA libraries was performed against the TAR RNA element of HIV-1. The role of the selection conditions on the outcome of the selection was evaluated by varying the magnesium concentration and the temperature. The selection stringency was demonstrated to determine i) the affinity of the best identified aptamers for the TAR target, and ii) the type of interaction between the two partners. Selections performed with a DNA library under low (4 degrees C, 10 mM magnesium) and high stringency (23 degrees C, 3 mM magnesium) led to the emergence of "kissing aptamers"; but even if the motif interacting directly with the TAR loop were identical in the two kinds of aptamers, the consensus was extended from eight to thirteen nucleotides when the Mg(2+) concentration was decreased from 10 to 3 mM. Similar kissing aptamers were selected at 23 degrees C and 37 degrees C starting with two different RNA libraries under identical ionic conditions. In addition, selection performed at 37 degrees C yielded a significant proportion of antisense sequences. Only antisense RNAs complementary to the TAR loop competitively inhibited the association of a Tat peptide with TAR.     

Structure of TAR RNA complexed with a Tat-TAR interaction nanomolar inhibitor that was identified by computational screening

HIV-1 TAR RNA functions critically in viral replication by binding the transactivating regulatory protein Tat. We recently identified several compounds that experimentally inhibit the Tat-TAR interaction completely at a 100 nM concentration. We used computational screening of the 181,000-compound Available Chemicals Directory against the three-dimensional structure of TAR [1]. Here we report the NMR-derived structure of TAR complexed with acetylpromazine. This structure represents a new class of compounds with good bioavailability and low toxicity that bind with high affinity to TAR. NMR data unambiguously show that acetylpromazine binds only to the unique 5' bulge site to which the Tat protein binds. Specificity and affinity of binding are conferred primarily by a network of base stacking and hydrophobic interactions. Acetylpromazine alters the structure of free TAR less than Tat peptides and neomycin do.     

Measurement of the effect of monovalent cations on RNA hairpin stability

Using optical tweezers, we have measured the effect of monovalent cation concentration and species on the folding free energy of five large (49-124 nt) RNA hairpins, including HIV-1 TAR and molecules approximating A.U and G.C homopolymers. RNA secondary structure thermodynamics are accurately described by a model consisting of nearest-neighbor interactions and additive loop and bulge terms. Melting of small (<15 bp) duplexes and hairpins in 1 M NaCl has been used to determine the parameters of this model, which is now used extensively to predict structure and folding dynamics. Few systematic measurements have been made in other ionic conditions or for larger structures. By applying mechanical force, we measured the work required to fold and unfold single hairpins at room temperature over a range of cation concentrations from 50 to 1000 mM. Free energies were then determined using the Crooks fluctuation theorem. We observed the following: (1) In most cases, the nearest-neighbor model accurately predicted the free energy of folding at 1 M NaCl. (2) Free energy was proportional to the logarithm of salt concentration. (3) Substituting potassium ions for sodium slightly decreased hairpin stability. The TAR hairpin also misfolded nearly twice as often in KCl, indicating a differential kinetic response. (4) Monovalent cation concentration affects RNA stability in a sequence-dependent manner. G.C helices were unaffected by changing salt concentration, A.U helices were modestly affected, and the hairpin loop was very sensitive. Surprisingly, the U.C.U bulge of TAR was found to be equally stable in all conditions tested. We also report a new estimate for the elastic parameters of single-stranded RNA.     

New RNA labeling methods offer dramatic sensitivity enhancements in 2H NMR relaxation spectra

A new labeling strategy is presented that greatly facilitates the measurement of 2H spin relaxation rates in RNA molecules as a probe of pico- to nanosecond time scale dynamics. In this labeling scheme the sugar positions are uniformly 13C-labeled, with position 2' protonated and all other sites on the sugar deuterated. Pulse sequences are presented for measurement of 2H R1 and R2 relaxation rates at positions 1', 3', and 4' with sensitivity gains that are on the order of 5-fold relative to previous methods that employed random fractional deuteration. The improved sensitivity is transformative and facilitates the study of motion in moderately sized RNA molecules with good sensitivity. The utility of the approach is demonstrated with an application to HIV-2 TAR, where the site-specific measures of molecular dynamics at sugar positions obtained here complement previous studies of dynamics at aromatic sites in the molecule.     

Native Top‐Down Mass Spectrometry of TAR RNA in Complexes with a Wild‐Type tat Peptide for Binding Site Mapping

Ribonucleic acids (RNA) frequently associate with proteins in many biological processes to form more or less stable complex structures. The characterization of RNA–protein complex structures and binding interfaces by nuclear magnetic resonance (NMR) spectroscopy, X‐ray crystallography, or strategies based on chemical crosslinking, however, can be quite challenging. Herein, we have explored the use of an alternative method, native top‐down mass spectrometry (MS), for probing of complex stoichiometry and protein binding sites at the single‐residue level of RNA. Our data show that the electrostatic interactions between HIV‐1 TAR RNA and a peptide comprising the arginine‐rich binding region of tat protein are sufficiently strong in the gas phase to survive phosphodiester backbone cleavage of RNA by collisionally activated dissociation (CAD), thus allowing its use for probing tat binding sites in TAR RNA by top‐down MS. Moreover, the MS data reveal time‐dependent 1:2 and 1:1 stoichiometries of the TAR–tat complexes and suggest structural rearrangements of TAR RNA induced by binding of tat peptide.     

Dynamics of large elongated RNA by NMR carbon relaxation

We present an NMR strategy for characterizing picosecond-to-nanosecond internal motions in uniformly 13C/15N-labeled RNAs that combines measurements of R1, R1rho, and heteronuclear 13C{1H} NOEs for protonated base (C2, C5, C6, and C8) and sugar (C1') carbons with a domain elongation strategy for decoupling internal from overall motions and residual dipolar coupling (RDC) measurements for determining the average RNA global conformation and orientation of the principal axis of the axially symmetric rotational diffusion. TROSY-detected pulse sequences are presented for the accurate measurement of nucleobase carbon R1 and R1rho rates in large RNAs. The relaxation data is analyzed using a model free formalism which takes into account the very high anisotropy of overall rotational diffusion (Dratio approximately 4.7), asymmetry of the nucleobase CSAs and noncollinearity of C-C, C-H dipolar and CSA interactions under the assumption that all interaction tensors for a given carbon experience identical isotropic internal motions. The approach is demonstrated and validated on an elongated HIV-1 TAR RNA (taum approximately 18 ns) both in free form and bound to the ligand argininamide (ARG). Results show that, while ARG binding reduces the amplitude of collective helix motions and local mobility at the binding pocket, it leads to a drastic increase in the local mobility of "spacer" bulge residues linking the two helices which undergo virtually unrestricted internal motions (S2 approximately 0.2) in the ARG bound state. Our results establish the ability to quantitatively study the dynamics of RNAs which are significantly larger and more anisotropic than customarily studied by NMR carbon relaxation.     

Target-Directed Azide-Alkyne Cycloaddition for Assembling HIV-1 TAR RNA Binding Ligands

The highly conserved HIV-1 transactivation response element (TAR) binds to the trans-activator protein Tat and facilitates viral replication in its latent state. The inhibition of Tat–TAR interactions by selectively targeting TAR RNA has been used as a strategy to develop potent antiviral agents. Therefore, HIV-1 TAR RNA represents a paradigmatic system for therapeutic intervention. Herein, we have employed biotin-tagged TAR RNA to assemble its own ligands from a pool of reactive azide and alkyne building blocks. To identify the binding sites and selectivity of the ligands, the in situ cycloaddition has been further performed using control nucleotide (TAR DNA and TAR RNA without bulge) templates. The hit triazole-linked thiazole peptidomimetic products have been isolated from the biotin-tagged target templates using streptavidin beads. The major triazole lead generated by the TAR RNA presumably binds in the bulge region, shows specificity for TAR RNA over TAR DNA, and inhibits Tat–TAR interactions.     

Solid-phase synthesis and evaluation of TAR RNA targeted beta-carboline-nucleoside conjugates

Four types of beta-carboline-nucleoside conjugates were synthesized. The binding affinities of these beta-carboline-nucleoside conjugates , and to TAR RNA were evaluated by affinity capillary electrophoresis. The data of binding affinities to TAR RNA show that conjugates and are stronger binders than the parent compound . Computer modeling indicates that the beta-carboline-nucleoside conjugate can fit to the UCU three-nucleotide bulge region of TAR RNA.     

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 CD₅₀ of ca. 1 M.     

Target‐Directed Azide‐Alkyne Cycloaddition for Assembling HIV‐1 TAR RNA Binding Ligands

The highly conserved HIV‐1 transactivation response element (TAR) binds to the trans‐activator protein Tat and facilitates viral replication in its latent state. The inhibition of Tat–TAR interactions by selectively targeting TAR RNA has been used as a strategy to develop potent antiviral agents. Therefore, HIV‐1 TAR RNA represents a paradigmatic system for therapeutic intervention. Herein, we have employed biotin‐tagged TAR RNA to assemble its own ligands from a pool of reactive azide and alkyne building blocks. To identify the binding sites and selectivity of the ligands, the in situ cycloaddition has been further performed using control nucleotide (TAR DNA and TAR RNA without bulge) templates. The hit triazole‐linked thiazole peptidomimetic products have been isolated from the biotin‐tagged target templates using streptavidin beads. The major triazole lead generated by the TAR RNA presumably binds in the bulge region, shows specificity for TAR RNA over TAR DNA, and inhibits Tat–TAR interactions.     

Adding energy minimization strategy to peptide‐design algorithm enables better search for RNA‐binding peptides: Redesigned λ N peptide binds boxB RNA

Our previously developed peptide‐design algorithm was improved by adding an energy minimization strategy which allows the amino acid sidechains to move in a broad configuration space during sequence evolution. In this work, the new algorithm was used to generate a library of 21‐mer peptides which could substitute for λ N peptide in binding to boxB RNA. Six potential peptides were obtained from the algorithm, all of which exhibited good binding capability with boxB RNA. Atomistic molecular dynamics simulations were then conducted to examine the ability of the λ N peptide and three best evolved peptides, viz. Pept01, Pept26, and Pept28, to bind to boxB RNA. Simulation results demonstrated that our evolved peptides are better at binding to boxB RNA than the λ N peptide. Sequence searches using the old (without energy minimization strategy) and new (with energy minimization strategy) algorithms confirm that the new algorithm is more effective at finding good RNA‐binding peptides than the old algorithm. © 2016 Wiley Periodicals, Inc.