The apical loop of the HIV-1 TAR RNA hairpin is stabilized by a cross-loop base pair

Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods

Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.

Conformational dynamics and energetics of viral RNA recognition by lab-evolved proteins

The conserved and structured elements in viral RNA genomes interact with proteins to regulate various events in the viral life cycle and have become key targets for developing novel therapeutic approaches. We probe physical interactions between lab-evolved proteins and a viral RNA element from the HIV-1 genome. Specifically, we study the role of an arginine-rich loop in recognition of designed proteins by the viral RNA element. We report free energy calculations to quantitatively estimate the protein/RNA binding energetics, focusing on the mutations of arginine residues involved in recognition of the major groove of RNA by proteins.

Strategies for Covalent Labeling of Long RNAs

The introduction of chemical modifications into long RNA molecules at specific positions for visualization, biophysical investigations, diagnostic and therapeutic applications still remains challenging. In this review, we present recent approaches for covalent internal labeling of long RNAs. Topics included are the assembly of large modified RNAs via enzymatic ligation of short synthetic oligonucleotides and synthetic biology approaches preparing site-specifically modified RNAs via in vitro transcription using an expanded genetic alphabet. Moreover, recent approaches to employ deoxyribozymes (DNAzymes) and ribozymes for RNA labeling and RNA methyltransferase based labeling strategies are presented. We discuss the potentials and limits of the individual methods, their applicability for RNAs with several hundred to thousands of nucleotides in length and indicate future directions in the field.

Targeting structural features of viral genomes with a nano-sized supramolecular drug

RNA targeting is an exciting frontier for drug design. Intriguing targets include functional RNA structures in structurally-conserved untranslated regions (UTRs) of many lethal viruses. However, computational docking screens, valuable in protein structure targeting, fail for inherently flexible RNA. Herein we harness MD simulations with Markov state modeling to enable nanosize metallo-supramolecular cylinders to explore the dynamic RNA conformational landscape of HIV-1 TAR untranslated region RNA (representative for many viruses) replicating experimental observations. These cylinders are exciting as they have unprecedented nucleic acid binding and are the first supramolecular helicates shown to have anti-viral activity in cellulo: the approach developed in this study provides additional new insight about how such viral UTR structures might be targeted with the cylinder binding into the heart of an RNA-bulge cavity, how that reduces the conformational flexibility of the RNA and molecular details of the insertion mechanism. The approach and understanding developed represents a new roadmap for design of supramolecular drugs to target RNA structural motifs across biology and nucleic acid nanoscience.

RNA Helicase A Regulates the Replication of RNA Viruses

The RNA helicase A (RHA) is a member of DExH-box helicases and characterized by two double-stranded RNA binding domains at the N-terminus. RHA unwinds double-stranded RNA in vitro and is involved in RNA metabolisms in the cell. RHA is also hijacked by a variety of RNA viruses to facilitate virus replication. Herein, this review will provide an overview of the role of RHA in the replication of RNA viruses.

Regulation of Promoter Proximal Pausing of RNA Polymerase II in Metazoans

Regulation of transcription is a tightly choreographed process. The establishment of RNA polymerase II promoter proximal pausing soon after transcription initiation and the release of Pol II into productive elongation are key regulatory processes that occur in early elongation. We describe the techniques and tools that have become available for the study of promoter proximal pausing and their utility for future experiments. We then provide an overview of the factors and interactions that govern a multipartite pausing process and address emerging questions surrounding the mechanism of RNA polymerase II's subsequent advancement into the gene body. Finally, we address remaining controversies and future areas of study.

HIV-1 Latency and Viral Reservoirs: Existing Reversal Approaches and Potential Technologies, Targets, and Pathways Involved in HIV Latency Studies

Eradication of latent human immunodeficiency virus (HIV) infection is a global health challenge. Reactivation of HIV latency and killing of virus-infected cells, the so-called "kick and kill" or "shock and kill" approaches, are a popular strategy for HIV cure. While antiretroviral therapy (ART) halts HIV replication by targeting multiple steps in the HIV life cycle, including viral entry, integration, replication, and production, it cannot get rid of the occult provirus incorporated into the host-cell genome. These latent proviruses are replication-competent and can rebound in cases of ART interruption or cessation. In general, a very small population of cells harbor provirus, serve as reservoirs in ART-controlled HIV subjects, and are capable of expressing little to no HIV RNA or proteins. Beyond the canonical resting memory CD4⁺ T cells, HIV reservoirs also exist within tissue macrophages, myeloid cells, brain microglial cells, gut epithelial cells, and hematopoietic stem cells (HSCs). Despite a lack of active viral production, latently HIV-infected subjects continue to exhibit aberrant cellular signaling and metabolic dysfunction, leading to minor to major cellular and systemic complications or comorbidities. These include genomic DNA damage; telomere attrition; mitochondrial dysfunction; premature aging; and lymphocytic, cardiac, renal, hepatic, or pulmonary dysfunctions. Therefore, the arcane machineries involved in HIV latency and its reversal warrant further studies to identify the cryptic mechanisms of HIV reservoir formation and clearance. In this review, we discuss several molecules and signaling pathways, some of which have dual roles in maintaining or reversing HIV latency and reservoirs, and describe some evolving strategies and possible approaches to eliminate viral reservoirs and, ultimately, cure/eradicate HIV infection.

RNA Drugs and RNA Targets for Small Molecules: Principles, Progress, and Challenges

RNA-based therapies, including RNA molecules as drugs and RNA-targeted small molecules, offer unique opportunities to expand the range of therapeutic targets. Various forms of RNAs may be used to selectively act on proteins, transcripts, and genes that cannot be targeted by conventional small molecules or proteins. Although development of RNA drugs faces unparalleled challenges, many strategies have been developed to improve RNA metabolic stability and intracellular delivery. A number of RNA drugs have been approved for medical use, including aptamers (e.g., pegaptanib) that mechanistically act on protein target and small interfering RNAs (e.g., patisiran and givosiran) and antisense oligonucleotides (e.g., inotersen and golodirsen) that directly interfere with RNA targets. Furthermore, guide RNAs are essential components of novel gene editing modalities, and mRNA therapeutics are under development for protein replacement therapy or vaccination, including those against unprecedented severe acute respiratory syndrome coronavirus pandemic. Moreover, functional RNAs or RNA motifs are highly structured to form binding pockets or clefts that are accessible by small molecules. Many natural, semisynthetic, or synthetic antibiotics (e.g., aminoglycosides, tetracyclines, macrolides, oxazolidinones, and phenicols) can directly bind to ribosomal RNAs to achieve the inhibition of bacterial infections. Therefore, there is growing interest in developing RNA-targeted small-molecule drugs amenable to oral administration, and some (e.g., risdiplam and branaplam) have entered clinical trials. Here, we review the pharmacology of novel RNA drugs and RNA-targeted small-molecule medications, with a focus on recent progresses and strategies. Challenges in the development of novel druggable RNA entities and identification of viable RNA targets and selective small-molecule binders are discussed. SIGNIFICANCE STATEMENT: With the understanding of RNA functions and critical roles in diseases, as well as the development of RNA-related technologies, there is growing interest in developing novel RNA-based therapeutics. This comprehensive review presents pharmacology of both RNA drugs and RNA-targeted small-molecule medications, focusing on novel mechanisms of action, the most recent progress, and existing challenges.

Human Immunodeficiency Virus-Associated Exosomes Promote Kaposi's Sarcoma-Associated Herpesvirus Infection via the Epidermal Growth Factor Receptor

Kaposi’s sarcoma-associated herpesvirus (KSHV) is the causal agent for Kaposi’s sarcoma (KS), the most common malignancy in HIV/AIDS patients. Oral transmission through saliva is considered the most common route for spreading the virus among HIV/AIDS patients. However, the role of HIV-specific components in the cotransfection of KSHV is unclear. We demonstrate that exosomes purified from the saliva of HIV-positive patients and secreted by HIV-infected T-cell lines promote KSHV infectivity in immortalized and primary oral epithelial cells. HIV-associated exosomes promote KSHV infection, which depends on HIV trans-activation response element (TAR) RNA and EGFR of oral epithelial cells, which can be targeted for reducing KSHV infection. These results reveal that HIV-associated exosomes are a risk factor for KSHV infection in the HIV-infected population.

Kaposi’s sarcoma-associated herpesvirus (KSHV) is the causal agent for Kaposi’s sarcoma (KS), the most common malignancy in people living with human immunodeficiency virus (HIV)/AIDS. The oral cavity is a major route for KSHV infection and transmission. However, how KSHV breaches the oral epithelial barrier for spreading to the body is not clear. Here, we show that exosomes purified from either the saliva of HIV-positive individuals or the culture supernatants of HIV-1-infected T-cell lines promote KSHV infectivity in immortalized and primary human oral epithelial cells. HIV-associated saliva exosomes contain the HIV trans-activation response element (TAR), Tat, and Nef RNAs but do not express Tat and Nef proteins. The TAR RNA in HIV-associated exosomes contributes to enhancing KSHV infectivity through the epidermal growth factor receptor (EGFR). An inhibitory aptamer against TAR RNA reduces KSHV infection facilitated by the synthetic TAR RNA in oral epithelial cells. Cetuximab, a monoclonal neutralizing antibody against EGFR, blocks HIV-associated exosome-enhanced KSHV infection. Our findings reveal that saliva containing HIV-associated exosomes is a risk factor for the enhancement of KSHV infection and that the inhibition of EGFR serves as a novel strategy for preventing KSHV infection and transmission in the oral cavity.

IMPORTANCE Kaposi’s sarcoma-associated herpesvirus (KSHV) is the causal agent for Kaposi’s sarcoma (KS), the most common malignancy in HIV/AIDS patients. Oral transmission through saliva is considered the most common route for spreading the virus among HIV/AIDS patients. However, the role of HIV-specific components in the cotransfection of KSHV is unclear. We demonstrate that exosomes purified from the saliva of HIV-positive patients and secreted by HIV-infected T-cell lines promote KSHV infectivity in immortalized and primary oral epithelial cells. HIV-associated exosomes promote KSHV infection, which depends on HIV trans-activation response element (TAR) RNA and EGFR of oral epithelial cells, which can be targeted for reducing KSHV infection. These results reveal that HIV-associated exosomes are a risk factor for KSHV infection in the HIV-infected population.

Face-time with TAR: Portraits of an HIV-1 RNA with diverse modes of effector recognition relevant for drug discovery

Small molecules and short peptides that potently and selectively bind RNA are rare, making the molecular structures of these complexes highly exceptional. Accordingly, several recent investigations have provided unprecedented structural insights into how peptides and proteins recognize the HIV-1 transactivation response (TAR) element, a 59-nucleotide-long, noncoding RNA segment in the 5' long terminal repeat region of viral transcripts. Here, we offer an integrated perspective on these advances by describing earlier progress on TAR binding to small molecules, and by drawing parallels to recent successes in the identification of compounds that target the hepatitis C virus internal ribosome entry site (IRES) and the flavin-mononucleotide riboswitch. We relate this work to recent progress that pinpoints specific determinants of TAR recognition by: (i) viral Tat proteins, (ii) an innovative lab-evolved TAR-binding protein, and (iii) an ultrahigh-affinity cyclic peptide. New structural details are used to model the TAR-Tat-super-elongation complex (SEC) that is essential for efficient viral transcription and represents a focal point for antiviral drug design. A key prediction is that the Tat transactivation domain makes modest contacts with the TAR apical loop, whereas its arginine-rich motif spans the entire length of the TAR major groove. This expansive interface has significant implications for drug discovery and design, and it further suggests that future lab-evolved proteins could be deployed to discover steric restriction points that block Tat-mediated recruitment of the host SEC to HIV-1 TAR.

Structural mechanism for HIV-1 TAR loop recognition by Tat and the super elongation complex

Promoter-proximal pausing by RNA polymerase II (Pol II) is a key regulatory step in human immunodeficiency virus-1 (HIV-1) transcription and thus in the reversal of HIV latency. By binding to the nascent transactivating response region (TAR) RNA, HIV-1 Tat recruits the human super elongation complex (SEC) to the promoter and releases paused Pol II. Structural studies of TAR interactions have been largely focused on interactions between the TAR bulge and the arginine-rich motif (ARM) of Tat. Here, the crystal structure of the TAR loop in complex with Tat and the SEC core was determined at a 3.5-Å resolution. The bound TAR loop is stabilized by cross-loop hydrogen bonds. It makes structure-specific contacts with the side chains of the Cyclin T1 Tat-TAR recognition motif (TRM) and the zinc-coordinating loop of Tat. The TAR loop phosphate backbone forms electrostatic and VDW interactions with positively charged side chains of the CycT1 TRM. Mutational analysis showed that these interactions contribute importantly to binding affinity. The Tat ARM was present in the crystallized construct; however, it was not visualized in the electron density, and the TAR bulge was not formed in the RNA construct used in crystallization. Binding assays showed that TAR bulge-Tat ARM interactions contribute less to TAR binding affinity than TAR loop interactions with the CycT1 TRM and Tat core. Thus, the TAR loop evolved to make high-affinity interactions with the TRM while Tat has three roles: scaffolding and stabilizing the TRM, making specific interactions through its zinc-coordinating loop, and making electrostatic interactions through its ARM.

Exosomes derived from HIV-1-infected cells promote growth and progression of cancer via HIV TAR RNA

People living with HIV/AIDS on antiretroviral therapy have increased risk of non-AIDS-defining cancers (NADCs). However, the underlying mechanism for development and progression of certain NADCs remains obscure. Here we show that exosomes released from HIV-infected T cells and those purified from blood of HIV-positive patients stimulate proliferation, migration and invasion of oral/oropharyngeal and lung cancer cells. The HIV transactivation response (TAR) element RNA in HIV-infected T-cell exosomes is responsible for promoting cancer cell proliferation and inducing expression of proto-oncogenes and Toll-like receptor 3 (TLR3)-inducible genes. These effects depend on the loop/bulge region of the molecule. HIV-infected T-cell exosomes rapidly enter recipient cells through epidermal growth factor receptor (EGFR) and stimulate ERK1/2 phosphorylation via the EGFR/TLR3 axis. Thus, our findings indicate that TAR RNA-containing exosomes from HIV-infected T cells promote growth and progression of particular NADCs through activation of the ERK cascade in an EGFR/TLR3-dependent manner.

HIV patients have an increased risk of developing non-AIDS-defining cancers but the molecular mechanisms underlying this predisposition are unclear. Here the authors show that exosomes secreted by HIV-infected T cells or isolated from the blood of HIV-positive patients, stimulate oncogenic properties of cancer cells through the activation of ERK1/2 signaling pathway.

HIV-1 Tat interactions with cellular 7SK and viral TAR RNAs identifies dual structural mimicry

The HIV Tat protein competes with the 7SK:HEXIM interaction to hijack pTEFb from 7SK snRNP and recruit it to the TAR motif on stalled viral transcripts. Here we solve structures of 7SK stemloop-1 and TAR in complex with Tat’s RNA binding domain (RBD) to gain insights into this process. We find that 7SK is peppered with arginine sandwich motifs (ASM)—three classical and one with a pseudo configuration. Despite having similar RBDs, the presence of an additional arginine, R52, confers Tat the ability to remodel the pseudo configuration, required for HEXIM binding, into a classical sandwich, thus displacing HEXIM. Tat also uses R52 to remodel the TAR bulge into an ASM whose structure is identical to that of the remodeled ASM in 7SK. Together, our structures reveal a dual structural mimicry wherein viral Tat and TAR have co-opted structural motifs present in cellular HEXIM and 7SK for productive transcription of its genome.

The HIV Tat protein recruits a host elongation factor from the cellular 7SK complex to the viral TAR RNA to ensure transcriptional elongation. Here, Pham et al. solve the structures of both 7SK and TAR RNAs in complex with Tat’s RNA binding domain and gain mechanistic insights into the process.

Structure of HIV TAR in complex with a Lab-Evolved RRM provides insight into duplex RNA recognition and synthesis of a constrained peptide that impairs transcription

Natural and lab-evolved proteins often recognize their RNA partners with exquisite affinity. Structural analysis of such complexes can offer valuable insight into sequence-selective recognition that can be exploited to alter biological function. Here, we describe the structure of a lab-evolved RNA recognition motif (RRM) bound to the HIV-1 trans-activation response (TAR) RNA element at 1.80 Å-resolution. The complex reveals a trio of arginines in an evolved β2-β3 loop penetrating deeply into the major groove to read conserved guanines while simultaneously forming cation-π and salt-bridge contacts. The observation that the evolved RRM engages TAR within a double-stranded stem is atypical compared to most RRMs. Mutagenesis, thermodynamic analysis and molecular dynamics validate the atypical binding mode and quantify molecular contributions that support the exceptionally tight binding of the TAR-protein complex (KD,App of 2.5 ± 0.1 nM). These findings led to the hypothesis that the β2-β3 loop can function as a standalone TAR-recognition module. Indeed, short constrained peptides comprising the β2-β3 loop still bind TAR (KD,App of 1.8 ± 0.5 μM) and significantly weaken TAR-dependent transcription. Our results provide a detailed understanding of TAR molecular recognition and reveal that a lab-evolved protein can be reduced to a minimal RNA-binding peptide.

Integrative, dynamic structural biology at atomic resolution--it's about time

Biomolecules adopt a dynamic ensemble of conformations, each with the potential to interact with binding partners or perform the chemical reactions required for a multitude of cellular functions. Recent advances in X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy and other techniques are helping us realize the dream of seeing--in atomic detail--how different parts of biomolecules shift between functional substates using concerted motions. Integrative structural biology has advanced our understanding of the formation of large macromolecular complexes and how their components interact in assemblies by leveraging data from many low-resolution methods. Here, we review the growing opportunities for integrative, dynamic structural biology at the atomic scale, contending there is increasing synergistic potential between X-ray crystallography, NMR and computer simulations to reveal a structural basis for protein conformational dynamics at high resolution.

Strong epistatic selection on the RNA secondary structure of HIV

A key question in evolutionary genomics is how populations navigate the adaptive landscape in the presence of epistasis, or interactions among loci. This problem can be directly addressed by studying the evolution of RNA secondary structures, for which there is constraint to maintain pairing between Watson-Crick (WC) sites. Replacement of a nucleotide at one site of a WC pair reduces fitness by disrupting binding, which can be restored via a compensatory replacement at the interacting site. Here, I present the first genome-scale analysis of epistasis on the RNA secondary structure of human immunodeficiency virus type 1 (HIV-1). Comparison of polymorphism frequencies at ancestrally conserved sites reveals that selection against replacements is ∼2.7 times stronger at WC than at non-WC sites, such that nearly 50% of constraint can be attributed to epistasis. However, almost all epistatic constraint is due to selection against conversions of WC pairs to unpaired (UP) nucleotides, whereas conversions to GU wobbles are only slightly deleterious. This disparity is also evident in pairs with second-site compensatory replacements; conversions from UP nucleotides to WC pairs increase median fitness by ∼4.2%, whereas conversions from GU wobbles to WC pairs only increase median fitness by ∼0.3%. Moreover, second-site replacements that convert UP nucleotides to GU wobbles also increase median fitness by ∼4%, indicating that such replacements are nearly as compensatory as those that restore WC pairing. Thus, WC peaks of the HIV-1 epistatic adaptive landscape are connected by high GU ridges, enabling the viral population to rapidly explore distant peaks without traversing deep UP valleys.

Epistasis is an evolutionary process in which the effect of a nucleotide at one site in the genome is dependent on the presence or absence of particular nucleotides at other sites in the genome. One of the simplest types of epistasis occurs between Watson-Crick (WC) nucleotides in RNA secondary structures, which are under constraint to maintain base-pairing. In this study, I examine the effects of mutations at WC sites in the RNA secondary structure of HIV-1. I show that while epistasis plays a major role in the evolution of the HIV-1 secondary structure, different types of mutations have variable effects on fitness. Therefore, by favoring certain mutational trajectories, HIV-1 can evolve rapidly despite strong epistatic constraint on its RNA secondary structure.

Characterizing RNA ensembles from NMR data with kinematic models

Functional mechanisms of biomolecules often manifest themselves precisely in transient conformational substates. Researchers have long sought to structurally characterize dynamic processes in non-coding RNA, combining experimental data with computer algorithms. However, adequate exploration of conformational space for these highly dynamic molecules, starting from static crystal structures, remains challenging. Here, we report a new conformational sampling procedure, KGSrna, which can efficiently probe the native ensemble of RNA molecules in solution. We found that KGSrna ensembles accurately represent the conformational landscapes of 3D RNA encoded by NMR proton chemical shifts. KGSrna resolves motionally averaged NMR data into structural contributions; when coupled with residual dipolar coupling data, a KGSrna ensemble revealed a previously uncharacterized transient excited state of the HIV-1 trans-activation response element stem-loop. Ensemble-based interpretations of averaged data can aid in formulating and testing dynamic, motion-based hypotheses of functional mechanisms in RNAs with broad implications for RNA engineering and therapeutic intervention.

Structural and thermodynamic signatures that define pseudotriloop RNA hairpins

Pseudotriloop (PTL) structures in RNAs have been recognized as essential elements in RNA folding and recognition of proteins. PTL structures are derived from hexaloops by formation of a cross-loop base pair leaving a triloop and 3' bulged out residue. Despite their common presence and functional importance, insufficient structural and thermodynamic data are available that can be used to predict formation of PTLs from sequence alone. Using NMR spectroscopy and UV-melting data we established factors that contribute to the formation and stability of PTL structures derived from hepatitis B virus and human foamy virus. The NMR data show that, besides the cross-loop base pair, also a 3' pyrimidine bulge and a G-C loop-closing base pair are primary determinants of PTL formation. By changing the G-C closing base pair into C-G, the PTL switches into a hexaloop. Comparison of these rules with regular triloop hairpins and PTLs from other sources is discussed as well as the conservation of a PTL in human foamy virus and other spumaretroviruses.

Visualizing transient low-populated structures of RNA

The visualization of RNA conformational changes has provided fundamental insights into how regulatory RNAs carry out their biological functions. The RNA structural transitions that have been characterized to date involve long-lived species that can be captured by structure characterization techniques. Here, we report the Nuclear Magnetic Resonance visualization of RNA transitions towards invisible ‘excited states’ (ES), which exist in too little abundance (2–13%) and for too short periods of time (45–250 μs) to allow structural characterization by conventional techniques. Transitions towards ESs result in localized rearrangements in base-pairing that alter building block elements of RNA architecture, including helix-junction-helix motifs and apical loops. The ES can inhibit function by sequestering residues involved in recognition and signaling or promote ATP-independent strand exchange. Thus, RNAs do not adopt a single conformation, but rather exist in rapid equilibrium with alternative ESs, which can be stabilized by cellular cues to affect functional outcomes.

Dominant role of the 5' TAR bulge in dimerization of HIV-1 genomic RNA, but no evidence of TAR-TAR kissing during in vivo virus assembly

The 5' untranslated region of HIV-1 genomic RNA (gRNA) contains two stem-loop structures that appear to be equally important for gRNA dimerization: the 57-nucleotide 5' TAR, at the very 5' end, and the 35-nucleotide SL1 (nucleotides 243-277). SL1 is well-known for containing the dimerization initiation site (DIS) in its apical loop. The DIS is a six-nucleotide palindrome. Here, we investigated the mechanism of TAR-directed gRNA dimerization. We found that the trinucleotide bulge (UCU24) of the 5' TAR has dominant impacts on both formation of HIV-1 RNA dimers and maturation of the formed dimers. The ΔUCU trinucleotide deletion strongly inhibited the first process and blocked the other, thus impairing gRNA dimerization as severely as deletion of the entire 5' TAR, and more severely than deletion of the DIS, inactivation of the viral protease, or most severe mutations in the nucleocapsid protein. The apical loop of TAR contains a 10-nucleotide palindrome that has been postulated to stimulate gRNA dimerization by a TAR-TAR kissing mechanism analogous to the one used by SL1 to stimulate dimerization. Using mutations that strongly destabilize formation of the TAR palindrome duplex, as well as compensatory mutations that restore duplex formation to a wild-type-like level, we found no evidence of TAR-TAR kissing, even though mutations nullifying the kissing potential of the TAR palindrome could impair dimerization by a mechanism other than hindering of SL1. However, nullifying the kissing potential of TAR had much less severe effects than ΔUCU. By not uncovering a dimerization mechanism intrinsic to TAR, our data suggest that TAR mutations exert their effect 3' of TAR, yet not on SL1, because TAR and SL1 mutations have synergistic effects on gRNA dimerization.

Bioinformatics for RNomics

Rapid improvements in high-throughput experimental technologies make it nowadays possible to study the expression, as well as changes in expression, of whole transcriptomes under different environmental conditions in a detailed view. We describe current approaches to identify genome-wide functional RNA transcripts (experimentally as well as computationally), and focus on computational methods that may be utilized to disclose their function. While genome databases offer a wealth of information about known and putative functions for protein-coding genes, functional information for novel non-coding RNA genes is almost nonexistent. This is mainly explained by the lack of established software tools to efficiently reveal the function and evolutionary origin of non-coding RNA genes. Here, we describe in detail computational approaches one may follow to annotate and classify an RNA transcript.

Efficient algorithms for probing the RNA mutation landscape

The diversity and importance of the role played by RNAs in the regulation and development of the cell are now well-known and well-documented. This broad range of functions is achieved through specific structures that have been (presumably) optimized through evolution. State-of-the-art methods, such as McCaskill's algorithm, use a statistical mechanics framework based on the computation of the partition function over the canonical ensemble of all possible secondary structures on a given sequence. Although secondary structure predictions from thermodynamics-based algorithms are not as accurate as methods employing comparative genomics, the former methods are the only available tools to investigate novel RNAs, such as the many RNAs of unknown function recently reported by the ENCODE consortium. In this paper, we generalize the McCaskill partition function algorithm to sum over the grand canonical ensemble of all secondary structures of all mutants of the given sequence. Specifically, our new program, RNAmutants, simultaneously computes for each integer k the minimum free energy structure MFE(k) and the partition function Z(k) over all secondary structures of all k-point mutants, even allowing the user to specify certain positions required not to mutate and certain positions required to base-pair or remain unpaired. This technically important extension allows us to study the resilience of an RNA molecule to pointwise mutations. By computing the mutation profile of a sequence, a novel graphical representation of the mutational tendency of nucleotide positions, we analyze the deleterious nature of mutating specific nucleotide positions or groups of positions. We have successfully applied RNAmutants to investigate deleterious mutations (mutations that radically modify the secondary structure) in the Hepatitis C virus cis-acting replication element and to evaluate the evolutionary pressure applied on different regions of the HIV trans-activation response element. In particular, we show qualitative agreement between published Hepatitis C and HIV experimental mutagenesis studies and our analysis of deleterious mutations using RNAmutants. Our work also predicts other deleterious mutations, which could be verified experimentally. Finally, we provide evidence that the 3′ UTR of the GB RNA virus C has been optimized to preserve evolutionarily conserved stem regions from a deleterious effect of pointwise mutations. We hope that there will be long-term potential applications of RNAmutants in de novo RNA design and drug design against RNA viruses. This work also suggests potential applications for large-scale exploration of the RNA sequence-structure network. Binary distributions are available at http://RNAmutants.csail.mit.edu/.

Evolution is a central concept in biology. This phenomenon can be observed at all levels of the organization of life—from single molecules to multicellular organisms. Here, we focus our attention on the implication of evolution at the level of nucleic acid sequences. In this context, RNA sequences presumably have been optimized by evolution to achieve specific functions. These functions are supported by a structure that can be determined using thermodynamics-based models and energy minimization techniques. In this work, we develop efficient algorithms for predicting energetically favorable mutations and study their impact on the stability of the structure. We use these techniques to reveal sequences under evolutionary pressure and design new methods to predict lethal mutations. Applications of our tool lead to a better understanding of the mutational process of some key regulatory elements of two important pathogenic RNA viruses—human immunodeficiency virus and hepatitis C virus.

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.

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.

Synthesis of 5-fluoropyrimidine nucleotides as sensitive NMR probes of RNA structure

Enzymatic synthesis methods for the fluorinated 5'-triphosphate analogues 5F-UTP and 5F-CTP have been developed to facilitate 19F-labeling of RNAs for biophysical studies. HIV-2 TAR RNAs were synthesized using these analogues by in vitro transcription reactions using T7 RNA polymerase. The uniform incorporation of 5F-U or 5F-C analogues into HIV-2 TAR RNA transcripts does not significantly alter the RNA structure or thermodynamic stability. Fluorine observed homonuclear 19F-19F and heteronuclear 19F-1H NOE experiments providing selective distance information are presented and discussed. The availability of efficient synthesis of 5F-UTP, and for the first time, 5F-CTP, will facilitate the use of 5F-labeled RNAs in structural, ligand binding, and dynamic studies of RNAs using the advantages of 19F-labeling.

iRED analysis of TAR RNA reveals motional coupling, long-range correlations, and a dynamical hinge

The HIV-1 transactivation response RNA element (TAR), which is essential to the lifecycle of the virus, has been suggested, based on NMR and hydrodynamic measurements, to undergo substantial, collective, structural dynamics that are important for its function. To deal with the significant coupling between overall diffusional rotation and internal motion expected to exist in TAR, here we utilize an isotropic reorientational eigenmode dynamics analysis of simulated molecular trajectories to obtain a detailed description of TAR dynamics and an accurately quantified pattern of dynamical correlations. The analysis demonstrates the inseparability of internal and overall motional modes, confirms the existence and reveals the nature of collective domain dynamics, and additionally reveals that the hinge for these motions is centered on residues U23, C24, and C41. Results also indicate the existence of long-range communication between the loop and the core of the RNA, and between the loop and the bulge. Additionally, the isotropic reorientational eigenmode dynamics analysis explains, from a dynamical perspective, several existing biochemical mutational studies and suggests new mutations for future structural dynamics studies.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

TAR-RNA recognition by a novel cyclic aminoglycoside analogue

The formation of the Tat-protein/TAR-RNA complex is a crucial step in the regulation of human immunodeficiency virus (HIV)-gene expression. To obtain full-length viral transcripts the Tat/TAR complex has to recruit the positive transcription elongation factor complex (P-EFTb), which interacts with TAR through its cyclin T1 (CycT1) component. Mutational studies identified the TAR hexanucleotide loop as a crucial region for contacting CycT1. Interfering with the interaction between the Tat/CycT1 complex and the TAR-RNA is an attractive strategy for the design of anti-HIV drugs. Positively charged molecules, like aminoglycosides or peptidomimetics, bind the TAR-RNA, disrupting the Tat/TAR complex. Here, we investigate the complex between the HIV-2 TAR-RNA and a neooligoaminodeoxysaccharide by NMR spectroscopy. In contrast to other aminoglycosides, this novel aminoglycoside analogue contacts simultaneously the bulge residues required for Tat binding and the A35 residue of the hexanucleotide loop. Upon complex formation, the loop region undergoes profound conformational changes. The novel binding mode, together with the easy accessibility of derivatives for the neooligoaminodeoxysaccharide, could open the way to the design of a new class of TAR-RNA binders, which simultaneously inhibit the formation of both the Tat/TAR binary complex and the Tat/TAR/CycT1 ternary complex by obstructing both the bulge and loop regions of the RNA.

New, extended hairpin form of the TAR-2 RNA domain points to the structural polymorphism at the 5' end of the HIV-2 leader RNA

The HIV-2 TAR RNA domain (TAR-2) plays a key role in the trans-activation of HIV-2 transcription as it is the target for the Tat-2 protein and several cell factors. Here, we show that the TAR-2 domain exists in vitro in two global, alternative forms: a new, extended hairpin form with two conformers and the already proposed branched hairpins form. This points strongly to the structural polymorphism of the 5′ end of the HIV-2 leader RNA. The evidence comes from the non-denaturing PAGE mobility assay, 2D structure prediction, enzymatic and Pb²⁺- or Mg²⁺-induced RNA cleavages. Existence of the TAR-2 extended form was further proved by the examination of engineered TAR-2 mutants stabilized either in the branched or extended structure. The TAR-2 extended form predominates with an increasing magnesium concentration. Gel retardation assays reveal that both TAR-2 wt and its mutant, unable to form branched structure, bind Tat-2 protein with comparable, high affinity, while RNA hairpins I and II, derived from TAR-2 branched structure model, show much less protein binding. We propose that an internal loop region of the TAR-2 extended hairpin form is a potential Tat-2 binding site.

Specific Interactions Between HIV-1 Nucleocapsid Protein and the TAR Element

During retroviral reverse transcription, the minus-strand strong-stop DNA (ss-cDNA) is transferred to the 3' end of the genomic RNA and this requires the repeat (R) sequences present at both ends of the genome. In vitro, the human immunodeficiency virus type 1 (HIV-1) R sequence can promote DNA strand transfer when present in ectopic internal positions. Using HIV-1 model systems, the R sequences and nucleocapsid protein (NC) were found to be key determinants of ss-cDNA transfer. To gain insights into specific interactions between HIV-1 NC and RNA and the influence of NC on R folding, we investigated the secondary structures of R in two natural contexts, namely at the 5' or 3' end of RNAs representing the terminal regions of the genome, and in two ectopic internal positions that also support efficient minus-strand transfer. To investigate the roles of NC zinc fingers and flanking basic domains in the NC/RNA interactions, we used NC mutants. Analyses of the viral RNA/NC complexes by chemical and enzymatic probings, and gel retardation assays were performed under conditions allowing ss-cDNA transfer by reverse transcriptase. We report that NC binds the TAR apical loop specifically in the four genetic contexts without changing the folding of the TAR hairpin and R region significantly, and this requires the NC zinc fingers. In addition, we show that efficient annealing of cTAR DNA to the 3' R relies on sequence complementarities between TAR and cTAR terminal loops. These findings suggest that the TAR apical loop in the acceptor RNA is the initiation site for the annealing reaction that is chaperoned by NC during the minus-strand transfer.

Structural Features of MicroRNA (miRNA) Precursors and Their Relevance to miRNA Biogenesis and Small Interfering RNA/Short Hairpin RNA Design

We have established the structures of 10 human microRNA (miRNA) precursors using biochemical methods. Eight of these structures turned out to be different from those that were computer-predicted. The differences localized in the terminal loop region and at the opposite side of the precursor hairpin stem. We have analyzed the features of these structures from the perspectives of miRNA biogenesis and active strand selection. We demonstrated the different thermodynamic stability profiles for pre-miRNA hairpins harboring miRNAs at their 5'- and 3'-sides and discussed their functional implications. Our results showed that miRNA prediction based on predicted precursor structures may give ambiguous results, and the success rate is significantly higher for the experimentally determined structures. On the other hand, the differences between the predicted and experimentally determined structures did not affect the stability of termini produced through "conceptual dicing." This result confirms the value of thermodynamic analysis based on mfold as a predictor of strand section by RNAi-induced silencing complex (RISC).

Role of the structure of the top half of HIV-1 cTAR DNA on the nucleic acid destabilizing activity of the nucleocapsid protein NCp7

The viral nucleic acid chaperone protein NCp7 of HIV-1 assists the two obligatory strand transfers required for the conversion of the genomic RNA into double-stranded DNA by reverse transcriptase. The first strand transfer necessitates the annealing of the early product of cDNA synthesis, the minus strand strong stop DNA (ss-cDNA) to the 3' end of the genomic RNA. The hybridization reaction involves regions containing imperfect stem-loop (SL) structures, namely the TAR RNA at the 3' end of the genomic RNA and the complementary sequence cTAR at the 3' end of ss-cDNA. To pursue the characterization of the interaction between NCp7 and cTAR DNA, we investigated by absorbance, steady-state and time-resolved fluorescence spectroscopy, the interaction of NCp7 with wild-type and mutated DNAs representing the top half of cTAR. NCp7 was found to activate the transient melting of this cTAR DNA structure but less efficiently than that of cTAR lower half. The NCp7-induced destabilization of cTAR top half is dependent upon the three nucleotides bulging out of the stem, which thus represent a melting initiation site. In contrast, despite its ability to bind NCp7, the top loop does not play any significant role in NCp7-mediated melting. Thermodynamic data further suggest that NCp7-mediated destabilization of this cTAR structure correlates with the free energy changes afforded by destabilizing motifs like loops and bulges within the SL secondary structure. Interestingly, since NCp7 melts only short double-stranded sequences, destabilizing motifs need to be regularly positioned along the genomic sequence in order to promote strand transfer and thus genetic recombination during proviral DNA synthesis.

Evidence for a base triple in the free HIV-1 TAR RNA

We propose the existence of a novel base triple in the HIV-1 TAR hairpin. This triple is supported by covariation of loop residue 31 with residue 22, which is part of an unusual base pair with U40 below the 3-nucleotide bulge. A set of mutants was constructed to test the involvement of bases A22, U31, and U40 in a triple interaction. RNA structure probing, trans-activation assays, and structure modeling are consistent with the existence of this base triple in a bent conformation of the free TAR element. However, disruption of the base triple does not affect binding of a Tat-derived peptide. We therefore compared the structure of free and Tat-bound TAR RNA by footprinting and site-specific cross-linking analyses. These studies indicate that the Tat arginine-rich motif, in addition to its known binding site at the bulge, is in close contact with U31 in the TAR loop. Because binding of Tat to TAR is known to coincide with the formation of a base triple with residues U23, A27, and U38, we hypothesize that Tat binding and the associated straightening of TAR triggers the disruption of the (A22-U40)U31 triple.