The crystal structure of the Rev binding element of HIV-1 reveals novel base pairing and conformational variability

Structure of HIV-1 RRE stem-loop II identifies two conformational states of the high-affinity Rev binding site

During HIV infection, specific RNA-protein interaction between the Rev response element (RRE) and viral Rev protein is required for nuclear export of intron-containing viral mRNA transcripts. Rev initially binds the high-affinity site in stem-loop II, which promotes oligomerization of additional Rev proteins on RRE. Here, we present the crystal structure of RRE stem-loop II in distinct closed and open conformations. The high-affinity Rev-binding site is located within the three-way junction rather than the predicted stem IIB. The closed and open conformers differ in their non-canonical interactions within the three-way junction, and only the open conformation has the widened major groove conducive to initial Rev interaction. Rev binding assays show that RRE stem-loop II has high- and low-affinity binding sites, each of which binds a Rev dimer. We propose a binding model, wherein Rev-binding sites on RRE are sequentially created through structural rearrangements induced by Rev-RRE interactions.

Structural and computational studies of HIV-1 RNA

Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.

Dynamic geometry design of cyclic peptide architectures for RNA structure

Designing inhibitors for RNA is still challenging due to the bottleneck of maintaining the binding interaction of inhibitor-RNA accompanied by subtle RNA flexibility. Thus, the current approach usually needs to screen thousands of candidate inhibitors for binding. Here, we propose a dynamic geometry design approach to enrich the hits with only a tiny pool of designed geometrically compatible scaffold candidates. First, our method uses graph-based tree decomposition to explore the complementarity rigid binding cyclic peptide and design the amino acid side chain length and charge to fit the RNA pocket. Then, we perform an energy-based dynamical network algorithm to optimize the inhibitor-RNA hydrogen bonds. Dynamic geometry-guided design yields successful inhibitors with low micromolar binding affinity scaffolds and experimentally competes with the natural RNA chaperone. The results indicate that the dynamic geometry method yields higher efficiency and accuracy than traditional methods. The strategy could be further optimized to design the length and chirality by adopting nonstandard amino acids and facilitating RNA engineering for biological or medical applications.

ADAR activation by inducing a syn conformation at guanosine adjacent to an editing site

ADARs (adenosine deaminases acting on RNA) can be directed to sites in the transcriptome by complementary guide strands allowing for the correction of disease-causing mutations at the RNA level. However, ADARs show bias against editing adenosines with a guanosine 5′ nearest neighbor (5′-GA sites), limiting the scope of this approach. Earlier studies suggested this effect arises from a clash in the RNA minor groove involving the 2-amino group of the guanosine adjacent to an editing site. Here we show that nucleosides capable of pairing with guanosine in a syn conformation enhance editing for 5′-GA sites. We describe the crystal structure of a fragment of human ADAR2 bound to RNA bearing a G:G pair adjacent to an editing site. The two guanosines form a G_syn:G_anti pair solving the steric problem by flipping the 2-amino group of the guanosine adjacent to the editing site into the major groove. Also, duplexes with 2′-deoxyadenosine and 3-deaza-2′-deoxyadenosine displayed increased editing efficiency, suggesting the formation of a G_syn:AH⁺_anti pair. This was supported by X-ray crystallography of an ADAR complex with RNA bearing a G:3-deaza dA pair. This study shows how non-Watson–Crick pairing in duplex RNA can facilitate ADAR editing enabling the design of next generation guide strands for therapeutic RNA editing.

Could gene therapy cure HIV?

The Human Immunodeficiency Virus (HIV)/Acquired Immune Deficiency Syndrome (AIDS) continues to be a major global public health issue, having claimed almost 33 million lives so far. According to the recent report of the World Health Organization (WHO) in 2019, about 38 million people are living with AIDS. Hence, finding a solution to overcome this life-threatening virus can save millions of lives. Scientists and medical doctors have prescribed HIV patients with specific drugs for many years. Methods such antiretroviral therapy (ART) or latency-reversing agents (LRAs) have been used for a while to treat HIV patients, however they have some side effects and drawbacks causing their application to be not quite successful. Instead, the application of gene therapy which refers to the utilization of the therapeutic delivery of nucleic acids into a patient's cells as a drug to treat disease has shown promising results to control HIV infection. Therefore, in this review, we will summarize recent advances in gene therapy approach against HIV.

Native mass spectrometry reveals the initial binding events of HIV-1 rev to RRE stem II RNA

Nuclear export complexes composed of rev response element (RRE) ribonucleic acid (RNA) and multiple molecules of rev protein are promising targets for the development of therapeutic strategies against human immunodeficiency virus type 1 (HIV-1), but their assembly remains poorly understood. Using native mass spectrometry, we show here that rev initially binds to the upper stem of RRE IIB, from where it is relayed to binding sites that allow for rev dimerization. The newly discovered binding region implies initial rev recognition by nucleotides that are not part of the internal loop of RRE stem IIB RNA, which was previously identified as the preferred binding region. Our study highlights the unique capability of native mass spectrometry to separately study the binding interfaces of RNA/protein complexes of different stoichiometry, and provides a detailed understanding of the mechanism of RRE/rev association with implications for the rational design of potential drugs against HIV-1 infection.

The HIV-1 RNA-binding protein rev facilitates nuclear export of viral RNA. Here, the authors use native mass spectrometry to study the interactions between rev-derived peptides and rev response elements of HIV-1 RNA, providing mechanistic insights into rev recognition and recruitment.

m6A minimally impacts the structure, dynamics, and Rev ARM binding properties of HIV-1 RRE stem IIB

N6-methyladenosine (m6A) is a ubiquitous RNA post-transcriptional modification found in coding as well as non-coding RNAs. m6A has also been found in viral RNAs where it is proposed to modulate host-pathogen interactions. Two m6A sites have been reported in the HIV-1 Rev response element (RRE) stem IIB, one of which was shown to enhance binding to the viral protein Rev and viral RNA export. However, because these m6A sites have not been observed in other studies mapping m6A in HIV-1 RNA, their significance remains to be firmly established. Here, using optical melting experiments, NMR spectroscopy, and in vitro binding assays, we show that m6A minimally impacts the stability, structure, and dynamics of RRE stem IIB as well as its binding affinity to the Rev arginine-rich-motif (ARM) in vitro. Our results indicate that if present in stem IIB, m6A is unlikely to substantially alter the conformational properties of the RNA. Our results add to a growing view that the impact of m6A on RNA depends on sequence context and Mg2+.

Dynamic ensemble of HIV-1 RRE stem IIB reveals non-native conformations that disrupt the Rev-binding site

The HIV-1 Rev response element (RRE) RNA element mediates the nuclear export of intron containing viral RNAs by forming an oligomeric complex with the viral protein Rev. Stem IIB and nearby stem II three-way junction nucleate oligomerization through cooperative binding of two Rev molecules. Conformational flexibility at this RRE region has been shown to be important for Rev binding. However, the nature of the flexibility has remained elusive. Here, using NMR relaxation dispersion, including a new strategy for directly observing transient conformational states in large RNAs, we find that stem IIB alone or when part of the larger RREII three-way junction robustly exists in dynamic equilibrium with non-native excited state (ES) conformations that have a combined population of ∼20%. The ESs disrupt the Rev-binding site by changing local secondary structure, and their stabilization via point substitution mutations decreases the binding affinity to the Rev arginine-rich motif (ARM) by 15- to 80-fold. The ensemble clarifies the conformational flexibility observed in stem IIB, reveals long-range conformational coupling between stem IIB and the three-way junction that may play roles in cooperative Rev binding, and also identifies non-native RRE conformational states as new targets for the development of anti-HIV therapeutics.

HIV Rev Assembly on the Rev Response Element (RRE): A Structural Perspective

HIV-1 Rev is an ~13 kD accessory protein expressed during the early stage of virus replication. After translation, Rev enters the nucleus and binds the Rev response element (RRE), a ~350 nucleotide, highly structured element embedded in the env gene in unspliced and singly spliced viral RNA transcripts. Rev-RNA assemblies subsequently recruit Crm1 and other cellular proteins to form larger complexes that are exported from the nucleus. Once in the cytoplasm, the complexes dissociate and unspliced and singly-spliced viral RNAs are packaged into nascent virions or translated into viral structural proteins and enzymes, respectively. Rev binding to the RRE is a complex process, as multiple copies of the protein assemble on the RNA in a coordinated fashion via a series of Rev-Rev and Rev-RNA interactions. Our understanding of the nature of these interactions has been greatly advanced by recent studies using X-ray crystallography, small angle X-ray scattering (SAXS) and single particle electron microscopy as well as biochemical and genetic methodologies. These advances are discussed in detail in this review, along with perspectives on development of antiviral therapies targeting the HIV-1 RRE.

Structural features of a 3' splice site in influenza a

Influenza A is an RNA virus with a genome of eight negative sense segments. Segment 7 mRNA contains a 3′ splice site for alternative splicing to encode the essential M2 protein. On the basis of sequence alignment and chemical mapping experiments, the secondary structure surrounding the 3′ splice site has an internal loop, adenine bulge, and hairpin loop when it is in the hairpin conformation that exposes the 3′ splice site. We report structural features of a three-dimensional model of the hairpin derived from nuclear magnetic resonance spectra and simulated annealing with restrained molecular dynamics. Additional insight was provided by modeling based on ¹H chemical shifts. The internal loop containing the 3′ splice site has a dynamic guanosine and a stable imino (cis Watson–Crick/Watson–Crick) GA pair. The adenine bulge also appears to be dynamic with the A either stacked in the stem or forming a base triple with a Watson–Crick GC pair. The hairpin loop is a GAAA tetraloop closed by an AC pair.

How does the long G·G* Watson-Crick DNA base mispair comprising keto and enol tautomers of the guanine tautomerise? The results of a QM/QTAIM investigation

The double proton transfer (DPT) in the long G·G* Watson-Crick base mispair (|C6N1(G*)N1C6(G)| = 36.4°; C1 symmetry), involving keto and enol tautomers of the guanine (G) nucleobase, along two intermolecular neighboring O6H···O6 (8.39) and N1···HN1 (6.14 kcal mol(-1)) H-bonds that were established to be slightly anti-cooperative, leads to its transformation into the G*·G base mispair through a single transition state (|C6N1N1C6| = 37.1°; C1), namely to the interconversion into itself. It was shown that the G·G* ↔ G*·G tautomerisation via the DPT is assisted by the third specific contact, that sequentially switches along the intrinsic reaction coordinate (IRC) in an original way: (G)N2H···N2(G*) H-bond (-25.13 to -10.37) → N2···N2 van der Waals contact (-10.37 to -9.23) → (G)N2···HN2(G*) H-bond (-9.23 to 0.79) → (G*)N2···HN2(G) H-bond (0.79 to 7.35 Bohr). The DPT tautomerisation was found to proceed through the asynchronous concerted mechanism by employing the QM/QTAIM approach and the methodology of the scans of the geometric, electron-topological, energetic, polar and NBO properties along the IRC. Nine key points, that can be considered as part of the tautomerisation repertoire, have been established and analyzed in detail. Furthermore, it was shown that the G·G* or G*·G base mispair is a thermodynamically and dynamically stable structure with a lifetime of 8.22 × 10(-10) s and all 6 low-frequency intermolecular vibrations are able to develop during this time span. Lastly, our results highlight the importance of the G·G* ↔ G*·G DPT tautomerisation, which can have implications for biological and chemical sensing applications.

Dioxaphosphorinane-constrained nucleic Acid dinucleotides as tools for structural tuning of nucleic acids

We describe a rational approach devoted to modulate the sugar-phosphate backbone geometry of nucleic acids. Constraints were generated by connecting one oxygen of the phosphate group to a carbon of the sugar moiety. The so-called dioxaphosphorinane rings were introduced at key positions along the sugar-phosphate backbone allowing the control of the six-torsion angles α to ζ defining the polymer structure. The syntheses of all the members of the D-CNA family are described, and we emphasize the effect on secondary structure stabilization of a couple of diastereoisomers of α,β-D-CNA exhibiting wether B-type canonical values or not.

The influenza A segment 7 mRNA 3' splice site pseudoknot/hairpin family

The 3′ splice site of the influenza A segment 7 transcript is utilized to produce mRNA for the critical M2 ion-channel protein. In solution a 63 nt fragment that includes this region can adopt two conformations: a pseudoknot and a hairpin. In each conformation, the splice site, a binding site for the SF2/ASF exonic splicing enhancer and a polypyrimidine tract, each exists in a different structural context. The most dramatic difference occurs for the splice site. In the hairpin the splice site is between two residues that are involved in a 2 by 2 nucleotide internal loop. In the pseudoknot, however, these bases are canonically paired within one of the pseudoknotted helices. The conformational switching observed in this region has implications for the regulation of splicing of the segment 7 mRNA. A measure of stability of the structures also shows interesting trends with respect to host specificity: avian strains tend to be the most stable, followed by swine and then human.

The 3' splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin

The 3′ splice site of influenza A segment 7 is used to produce mRNA for the M2 ion-channel protein, which is critical to the formation of viable influenza virions. Native gel analysis, enzymatic/chemical structure probing, and oligonucleotide binding studies of a 63 nt fragment, containing the 3′ splice site, key residues of an SF2/ASF splicing factor binding site, and a polypyrimidine tract, provide evidence for an equilibrium between pseudoknot and hairpin structures. This equilibrium is sensitive to multivalent cations, and can be forced towards the pseudoknot by addition of 5 mM cobalt hexammine. In the two conformations, the splice site and other functional elements exist in very different structural environments. In particular, the splice site is sequestered in the middle of a double helix in the pseudoknot conformation, while in the hairpin it resides in a two-by-two nucleotide internal loop. The results suggest that segment 7 mRNA splicing can be controlled by a conformational switch that exposes or hides the splice site.

Role of wobble base pair geometry for codon degeneracy: purine-type bases at the anticodon wobble position

Codon degeneracy is a key feature of the genetic code, explained by Crick (J Mol Biol 19:548-555, 1966) in terms of imprecision of base pairing at the codon third position (the wobble position) of the codon-anticodon duplex. The Crick wobble rules define, but do not explain, which base pairs are allowed/disallowed at the wobble position of this duplex. This work examines whether the H-bonded configurations of solitary RNA base pairs can in themselves help decide which base pairs are allowed at the wobble position during codon-anticodon pairing. Taking the purine-type bases guanine, hypoxanthine, queuine and adenine as anticodon wobble bases, H-bonded pairing energies and optimized configurations of numerous RNA base pairs are calculated in gas and modeled aqueous phase at the B3LYP/6-31 G(d,p) level. Calculated descriptors of alignment of these solitary base pairs are able to screen between allowed and disallowed base pairs for all cases studied here, except two cases which invoke base-sugar interactions in the codon wobble nucleoside. The exclusion of adenine from the anticodon wobble position cannot be explained on the basis of pairing facility or base pair geometry. These DFT results thus account for the specificity and degeneracy of the genetic code for all cases involving guanine, hypoxanthine and queuine as anticodon wobble bases.

Extensive purifying selection acting on synonymous sites in HIV-1 Group M sequences

Background

Positive selection pressure acting on protein-coding sequences is usually inferred when the rate of nonsynonymous substitution is greater than the synonymous rate. However, purifying selection acting directly on the nucleotide sequence can lower the synonymous substitution rate. This could result in false inference of positive selection because when synonymous changes at some sites are under purifying selection, the average synonymous rate is an underestimate of the neutral rate of evolution. Even though HIV-1 coding sequences contain a number of regions that function at the nucleotide level, and are thus likely to be affected by purifying selection, studies of positive selection assume that synonymous substitutions can be used to estimate the neutral rate of evolution.

Results

We modelled site-to-site variation in the synonymous substitution rate across coding regions of the HIV-1 genome. Synonymous substitution rates were found to vary significantly within and between genes. Surprisingly, regions of the genome that encode proteins in more than one frame had significantly higher synonymous substitution rates than regions coding in a single frame. We found evidence of strong purifying selection pressure affecting synonymous mutations in fourteen regions with known functions. These included an exonic splicing enhancer, the rev-responsive element, the poly-purine tract and a transcription factor binding site. A further five highly conserved regions were located within known functional domains. We also found four conserved regions located in env and vpu which have not been characterized previously.

Conclusion

We provide the coordinates of genomic regions with markedly lower synonymous substitution rates, which are putatively under the influence of strong purifying selection pressure at the nucleotide level as well as regions encoding proteins in more than one frame. These regions should be excluded from studies of positive selection acting on HIV-1 coding regions.

Noncanonical G(syn)-G(anti) base pairs stabilized by sulphate anions in two X-ray structures of the (GUGGUCUGAUGAGGCC) RNA duplex

The structures of two crystal forms of the RNA 16-mer with the sequence GUGGUCUGAUGAGGCC, grown in the presence of a high concentration of sulphate ions, have been determined using synchrotron radiation at 1.4- and 2.0-A resolution. RNA with this sequence is known as one of the two strands of the noncleavable form of the hammerhead ribozyme. In both crystal structures, two G(syn)-G(anti) noncanonical base pairs are observed in the middle of a 14 base-pair (bp) duplex having 5'-dangling GU residues. Both structures contain sulphate anions interacting with the G-G bp stabilizing G in its syn conformation and bridging the two RNA strands. In both cases the interactions take place in the major groove, although the anions are accommodated within different helix geometries, most pronounced in the changing width of the major groove. In one structure, where a single sulphate spans both G-G pairs, the major groove is closed around the anion, while in the other structure, where each of the two G-G pairs is associated with a separate sulphate, the groove is open. This work provides the first examples of a G-G pair in syn-anti conformation, which minimizes the purine-purine clash in the center of the duplex, while utilizing its residual hydrogen bonding potential in specific interactions with sulphate anions.

The 1.3 A resolution structure of the RNA tridecamer r(GCGUUUGAAACGC): metal ion binding correlates with base unstacking and groove contraction

Metal ions play a key role in RNA folding and activity. Elucidating the rules that govern the binding of metal ions is therefore an essential step for better understanding the RNA functions. High-resolution data are a prerequisite for a detailed structural analysis of ion binding on RNA and, in particular, the observation of monovalent cations. Here, the high-resolution crystal structures of the tridecamer duplex r(GCGUUUGAAACGC) crystallized under different conditions provides new structural insights on ion binding on GAAA/UUU sequences that exhibit both unusual structural and functional properties in RNA. The present study extends the repertory of RNA ion binding sites in showing that the two first bases of UUU triplets constitute a specific site for sodium ions. A striking asymmetric pattern of metal ion binding in the two equivalent halves of the palindromic sequence demonstrates that sequence and its environment act together to bind metal ions. A highly ionophilic half that binds six metal ions allows, for the first time, the observation of a disodium cluster in RNA. The comparison of the equivalent halves of the duplex provides experimental evidences that ion binding correlates with structural alterations and groove contraction.

Base-tetrad swapping results in dimerization of RNA quadruplexes: implications for formation of the i-motif RNA octaplex

Nucleic acids adopt different multistranded helical architectures to perform various biological functions. Here, we report a crystal structure of an RNA quadruplex containing "base-tetrad swapping" and bulged nucleotide at 2.1-Angstroms resolution. The base-tetrad swapping results in a dimer of quadruplexes with an intercalated octaplex fragment at the 5' end junction. The intercalated base tetrads provide the basic repeat unit for constructing a model of intercalated RNA octaplex. The model we obtained shows fundamentally different characteristics from duplex, triplex, and quadruplex. We also observed two different orientations of bulged uridine residues that are related to the interaction with surroundings. This structural evidence reflects the conformational flexibility of bulged nucleotides in RNA quadruplexes and implies the potential roles of bulged nucleotides as recognition and interaction sites in RNA-protein and RNA-RNA interactions.

Crystal structure of a bulged RNA tetraplex at 1.1 a resolution: implications for a novel binding site in RNA tetraplex

Bulges are an important structural motif in RNA and can be used as recognition and interaction sites in RNA-protein interaction and RNA-RNA interaction. Here we report the first crystal structure of a bulged RNA tetraplex at 1.1 A resolution. The hexamer r(U)(BrdG)r(UGGU) forms a parallel tetraplex with the uridine sandwiched by guanines bulging out. The bulged uridine adopts the syn glycosidic conformation and its O2 and N3 atoms face outwards, serving as an effective recognition and interaction site. The bulge formation both widens the groove width and changes the groove hydrogen-bonding pattern on its 5' side. However, the bulge does not make any bends or kinks in the tetraplex structure. The present study demonstrates the dramatic difference between uridine and guanine in forming tetraplex structure. In addition, both G(syn) tetrad and G(anti) tetrad have been observed. They display the same base-pairing pattern and similar C1'-C1' distance but different hydrogen-bonding patterns in the groove.

Selenium derivatization and crystallization of DNA and RNA oligonucleotides for X-ray crystallography using multiple anomalous dispersion

We report here the solid phase synthesis of RNA and DNA oligonucleotides containing the 2'-selenium functionality for X-ray crystallography using multiwavelength anomalous dispersion. We have synthesized the novel 2'-methylseleno cytidine phosphoramidite and improved the accessibility of the 2'-methylseleno uridine phosphoramidite for the synthesis of many selenium-derivatized DNAs and RNAs in large scales. The yields of coupling these Se-nucleoside phosphoramidites into DNA or RNA oligonucleotides were over 99% when 5-(benzylmercapto)-1H-tetrazole was used as the coupling reagent. The UV melting study of A-form dsDNAs indicated that the 2'-selenium derivatization had no effect on the stability of the duplexes with the 3'-endo sugar pucker. Thus, the stems of functional RNA molecules with the same 3'-endo sugar pucker appear to be the ideal sites for the selenium derivatization with 2'-Se-C and 2'-Se-U. Crystallization of the selenium-derivatized oligonucleotides is also reported here. The results demonstrate that this 2'-selenium functionality is suitable for RNA and A-form DNA derivatization in X-ray crystallography.

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

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

Rev response elements (RRE) in lentiviruses: an RNAMotif algorithm-based strategy for RRE prediction

Lentiviruses (a sub-family of the retroviridae family) include primate and non-primate viruses associated with chronic diseases of the immune system and the central nervous system. All lentiviruses encode a regulatory protein Rev that is essential for post-transcriptional transport of the unspliced and incompletely spliced viral mRNAs from nuclei to cytoplasm. The Rev protein acts via binding to an RNA structural element known as the Rev responsive element (RRE). The RRE location and structure and the mechanism of the Rev-RRE interaction in primate and non-primate lentiviruses have been analyzed and compared. Based on structural data available for RRE of HIV-1, a two step computational strategy for prediction of putative RRE regions in lentivirus genomes has been developed. First, the RNAMotif algorithm was used to search genomic sequence for highly structured regions (HSR). Then the program RNAstructure, version 3.6 was used to calculate the structure and thermodynamic stability of the region of approximately 350 nucleotides encompassing the HSR. Our strategy correctly predicted the locations of all previously reported lentivirus RREs. We were able also to predict the locations and structures of potential RREs in four additional lentiviruses.

Structure and stability of wild-type and mutant RNA internal loops from the SL-1 domain of the HIV-1 packaging signal

The packaging signal (Psi) of the human immunodeficiency virus type 1 (HIV-1) enables encapsidation of the full-length genomic RNA against a background of a vast excess of cellular mRNAs. The core HIV-1 Psi is approximately 109 nucleotides and contains sequences critical for viral genomic dimerisation and splicing, in addition to the packaging signal. It consists of a series of stem-loops (termed SL-1 to SL-4), which can be arranged in a cloverleaf secondary structure. Using a combination of NMR spectroscopy, UV melting experiments, molecular modeling and phylogenetic analyses, we have explored the structure of two conserved internal loops proximal to the palindromic sequence of SL-1. Internal loop A, composed of six purines, forms a flexible structure that is strikingly similar to the Rev responsive element motif when bound to Rev protein. This result suggests that it may function as a protein-binding site. The absolutely conserved four-purine internal loop B is instead conformationally and thermodynamically unstable, and exhibits multiple conformations in solution. By introducing a double AGG to GGA mutation within this loop, its conformation is stabilised to form a new intra-molecular G:A:G base-triplet. The structure of the GGA mutant explains the relative instability of the wild-type loop. In a manner analogous to SL-3, we propose that conformational flexibility at this site may facilitate melting of the structure during Gag protein capture or genomic RNA dimerisation.

Protein-dependent ribozymes report molecular interactions in real time

Most approaches to monitoring interactions between biological macromolecules require large amounts of material, rely upon the covalent modification of an interaction partner, or are not amenable to real-time detection. We have developed a generalizable assay system based on interactions between proteins and reporter ribozymes. The assay can be configured in a modular fashion to monitor the presence and concentration of a protein or of molecules that modulate protein function. We report two applications of the assay: screening for a small molecule that disrupts protein binding to its nucleic acid target and screening for protein protein interactions. We screened a structurally diverse library of antibiotics for small molecules that modulate the activity of HIV-1 Rev-responsive ribozymes by binding to Rev. We identified an inhibitor that subsequently inhibited HIV-1 replication in cells. A simple format switch allowed reliable monitoring of domain-specific interactions between the blood-clotting factor thrombin and its protein partners. The rapid identification of interactions between proteins or of compounds that disrupt such interactions should have substantial utility for the drug-discovery process.

Solid-state NMR data support a helix-loop-helix structural model for the N-terminal half of HIV-1 Rev in fibrillar form

Rev is a 116 residue basic protein encoded by the genome of human immunodeficiency virus type 1 (HIV-1) that binds to multiple sites in the Rev response element (RRE) of viral mRNA transcripts in nuclei of host cells, leading to transport of incompletely spliced and unspliced viral mRNA to the cytoplasm of host cells in the latter phases of the HIV-1 life cycle. Rev is absolutely required for viral replication. Because Rev aggregates and fibrillizes in solution at concentrations required for crystal growth or liquid state NMR measurements, high-resolution structural characterization of full-length Rev has not been possible. Previously, circular dichroism studies have shown that approximately 50 % of the Rev sequence adopts helical secondary structure, predicted to correspond to a helix-loop-helix structural motif in the N-terminal half of the protein. We describe the application of solid-state NMR techniques to Rev fibrils as a means of obtaining site-specific, atomic-level structural constraints without requiring a high degree of solubility or crystallinity. Solid-state NMR measurements, using the double-quantum chemical shift anisotropy and constant-time double-quantum-filtered dipolar recoupling techniques, provide constraints on the phi and psi backbone dihedral angles at sites in which consecutive backbone carbonyl groups are labeled with (13)C. Quantitative analysis of the solid-state NMR data, by comparison with numerical simulations, indicates helical phi and psi angles at residues Leu13 and Val16 in the predicted helix 1 segment, and at residues Arg39, Arg 42, Arg43, and Arg44 in the predicted helix 2 segment. These data represent the first site-specific structural constraints from NMR spectroscopy on full-length Rev, and support the helix-loop-helix structural model for its N-terminal half.

Merged screening for human immunodeficiency virus Tat and Rev inhibitors

In addition to "conventional" drug discovery targets used in modern strategies, mainly focusing on proteins, recent insights into gene regulation as a novel drug concept have begun to invite the targeting of biomolecular interactions between proteins and RNA. Because two protein-RNA interactions (Tat and trans-activation-responsive element, Rev and Rev-responsive element) are essential for any productive replication of human immunodeficiency virus, this important human pathogen was used as a model system for our studies. The design of a fluorescence-based high throughput assay, in which both targets were presented in the same vessel, enabled us to simultaneously interrogate two characteristics of a potential inhibitor: potency of interference and selectivity toward each of the interactions. Although related systems have been reported for several DNA binders, an extension into interference with transcription events would open a new dimension of cellular regulation. Here we describe the setup of the screening assay for over 110,000 compounds as well as a primary characterization of identified hits. The assay's characteristics demonstrate that a microwell-based dual screening system for RNA binders may add a powerful tool to modern drug discovery.

X-ray crystallographic observation of "in-line" and "adjacent" conformations in a bulged self-cleaving RNA/DNA hybrid

The RNA strand in an RNA/DNA duplex with unpaired ribonucleotides can undergo self-cleavage at bulge sites in the presence of a variety of divalent metal ions (Hüsken et al., Biochemistry, 1996, 35:16591-16600). Transesterification proceeds via an in-line mechanism, with the 2'-OH of the bulged nucleotide attacking the 3'-adjacent phosphate group. The site-specificity of the reaction is most likely a consequence of the greater local conformational freedom of the RNA backbone in the bulge region. A standard A-form backbone geometry prohibits formation of an in-line arrangement between 2'-oxygen and phosphate. However, the backbone in the region of an unpaired nucleotide appears to be conducive to an in-line approach. Therefore, the bulge-mediated phosphoryl transfer reaction represents one of the simplest RNA self-cleavage systems. Here we focus on the conformational features of the RNA that underlie site-specific cleavage. The structures of an RNA/DNA duplex with single ribo-adenosyl bulges were analyzed in two crystal forms, permitting observation of 10 individual conformations of the RNA bulge moiety. The bulge geometries cover a range of relative arrangements between the 2'-oxygen of the bulged nucleotide and the P-O5' bond (including adjacent and near in-line) and give a detailed picture of the conformational changes necessary to line up the 2'-OH nucleophile and scissile bond. Although metal ions are of crucial importance in the catalysis of analogous cleavage reactions by ribozymes, it is clear that local strain or conformational flexibility in the RNA also affect cleavage selectivity and rate (Soukup & Breaker, RNA, 1999, 5:1308-1325). The geometries of the RNA bulges frozen out in the crystals provide snapshots along the reaction pathway prior to the transition state of the phosphoryl transfer reaction.

The chromosomal protein HMG-D binds to the TAR and RBE RNA of HIV-1

The high mobility group protein HMG-D is known to bind preferentially to DNA of irregular structures with little or no sequence specificity. Upon binding to DNA, this HMG-box protein widens the minor groove of the double helix and induces a significant bending of the helix. We show here that HMG-D can strongly bind to double-stranded RNA. Electrophoretic mobility shift assays show that HMG-D100 interacts with the transactivation response region (TAR) RNA from HIV-1. Strong interaction with a high affinity Rev protein binding element (RBE) RNA was also characterized. Gel shift experiments performed with several TAR RNA constructs lacking the lateral pyrimidine bulge or with modified apical loop regions indicate that the protein does not recognize the single-strand domains of the RNA but apparently interacts directly with the double-stranded stem regions. No protein-RNA complexes could be detected when using single-stranded oligoribonucleotides. HMG-D protein could bind to the wide minor groove of the A-form TAR RNA. The comparison of the amino acid sequence of HMG-D with that of known RNA binding proteins suggests that the interaction of the protein with a double-stranded RNA implicates the basic region of HMG-D as well as its HMG-box domain. From the in vitro data reported here, we propose a novel functional role for proteins of the HMG-1 family. The results suggest that architectural HMG proteins can be recruited by double-stranded RNA for the development of HIV-1 in the host cell.