Modeling the dynamics of the solvated SL1 domain of HIV-1 genomic RNA

A mode-coupling solution of the Smoluchowski diffusion equation (MCD theory), designed to describe the dynamics of wobbling macromolecules in water, is applied to a macromolecular bead model including water beads in the nearest layers. The necessary statistical averages are evaluated by time averaging along a molecular dynamics (MD) trajectory where both solute and water are introduced as atomistic models. The cross peaks in (1)H nuclear Overhauser effect spectroscopy (NOESY) NMR spectra that are routinely measured to determine biological structures are here calculated for the mutated 23 nucleotides stem-loop fragment of the SL1 domain in the HIV-1(Lai) genomic RNA. The calculations are in acceptable agreement with experiments without requiring any screening of the hydrodynamic interactions. The screening of hydrodynamics was necessary in previous MCD calculations obtained by using the same full atomistic MD trajectory, but a nonsolvated frictional model.     

Insight into the intrinsic flexibility of the SL1 stem-loop from genomic RNA of HIV-1 as probed by molecular dynamics simulation

The SL1 stem-loop is the dimerization initiation site for linking the two copies of the RNA forming the HIV-1 genome. The 26 nucleotides stem contains a defect consisting on a highly conserved G-rich 1-3 asymmetrical internal loop, which is a major site for nucleocapsid protein binding. Several NMR attempts were undertaken to determine the internal loop structure in the SL1 monomer. However, the RNA constructs used in the different studies were largely mutated, in particular with replacement of the nine nucleotides apical loop by a tetraloop, and divergent results were obtained ranging from a rigid structure to a particularly large flexibility. To investigate the reasons for such discrepancies, we used molecular dynamics simulation of the SL1 monomer to probe the effect of mutations on the conformational stability of the internal loop and of the whole stem. It is found that in the wild-type sequence, the internal loop displays conformational variability originating mainly from the nine nucleotides apical loop flexibility that causes large conformational fluctuations (without changing the average structure) in the 7 bp duplex linking the apical and internal loops. The large amplitude atomic motions in the duplex are transmitted to the internal loop in which they induce conformational changes characterized by a labile hydrogen bond network such as G5 successively H-bonded to A29 and G30. On the contrary, with a four nucleotides apical loop, conformational fluctuations in the duplex are reduced by a factor of 2 and are not sufficiently energizing for promoting changes in the internal loop structure at the time scale of the simulations.     

Molecular Dynamics Simulation for Probing the Flexibility of the 35 Nucleotide SL1 Sequence Kissing Complex from HIV-1Lai Genomic RNA

Abstract

The SL1 stem-loop located in the encapsidation domain is responsible for initiating the dimerisation of HIV-1 genomic RNA by means of a loop-loop interaction known as Kissing Complex (KC). The SL1 secondary structure has been predicted as a 35 nucleotides [K. G. Murti, M. Bondurant, and A. Tereba. J Virol 37, 411–419 (1981)] stem-loop composed of a 4 base pairs (bp) terminal duplex, a 4 nt asymmetrical internal loop, a 7 bp internal duplex, and a 9 nt apical loop. Several high resolution structures of the monomer and of KC of a 23 nt sequence containing only the internal duplex and the apical loop of SL1 are available in the literature. No experimental high resolution structure of the complete native SL1 sequence has been reported so far, either for the monomer or for KC. The asymmetrical internal loop has been described from NMR studies of different monomeric hairpin sequences, leading to divergent results, which suggests its high flexibility. In this work, we built a SL1₃₅ KC model which was submitted to a 31 ns molecular dynamics simulation (MD).

Our results allows to describe the internal dynamics of SL1₃₅ KC and the differences of behavior of the different parts of the dimer. Thus, we could show the stability of the interactions between the two apical loops and of the terminal duplexes, the destabilization of the internal duplexes and the high flexibility of the asymmetrical internal loops.     

On the stability of different experimental dimeric structures of the SL1 sequence from the genomic RNA of HIV-1 in solution: a molecular dynamics simulation and electrophoresis study

SL1 is a stem-loop RNA sequence from the genome of HIV-1 thought to be the initiation site for the dimerization of the retroviral genomic RNA. The aim of this study is to check the stability in solution of different experimental dimeric structures available in the literature. Two kinds of dimer have been evidenced: an extended duplex looking like a double helix with two internal bulges and a kissing complex in which the monomers with a stem/loop conformation are linked by intermolecular loop-loop interactions. Two divergent experimental structures of the kissing complex from the Lai isolate are reported in the literature, one obtained from NMR (Mujeeb et al., Nature Structural Biology, 1998, Vol. 5, pp. 432-436) and the other one from x-ray crystallography (Ennifar et al., Nature Structural Biology, 2001, Vol. 8, pp. 1064-1068). A crystallographic structure of the Mal isolate was also reported (Ennifar et al., Nature Structure Biology, 2001, Vol. 8, pp. 1064-1068). Concerning the extended duplex, a NMR structure is available for Lai (Girard et al., Journal of Biomolecular Structure and Dynamics, 1999, Vol. 16, pp. 1145-1157) and a crystallographic structure for Mal (Ennifar et al., Structure, 1999, Vol. 7, pp. 1439-1449). Using a molecular dynamics technique, all these experimental structures have been simulated in solution with explicit water and counterions. We show that both extended duplex structures are stable. On the contrary, the crystallographic structures of the Lai and Mal kissing complexes are rapidly destabilized in aqueous environment. Finally, the NMR structure of the Lai loop-loop kissing complex remains globally stable over a 20 ns MD simulation, although large rearrangements occur at the level of the stem/loop junctions that are flexible, as shown from free energy calculations. These results are compared to electrophoresis experiments on dimer formation.     

Conformational pathway for the kissing complex-->extended dimer transition of the SL1 stem-loop from genomic HIV-1 RNA as monitored by targeted molecular dynamics techniques

HIV-1 retroviral genomic RNA dimerization is initiated by loop-loop interactions between the SL1 stem-loops of two identical RNA molecules. The SL1-SL1 unstable resulting kissing complex (KC) then refolds irreversibly into a more stable complex called extended dimer (ED). Although the structures of both types of complex have been determined, very little is known about the conformational pathway corresponding to the transition, owing to the difficulty of observing experimentally intermediate conformations. In this study, we applied targeted molecular dynamics simulation techniques (TMD) to the phosphorus atoms for monitoring this pathway for the backbone, and a two-step strategy was adopted. In a first step, called TMD(-1), the dimer structure was constrained to progressively move away from KC without indicating the direction, until the RMSD from KC reaches 36A. A total of 20 TMD(-1) simulations were performed under different initial conditions and different simulation parameters. For RMSD ranging between 0 and 22A, the whole set of TMD(-1) simulations follows a similar pathway, then divergences are observed. None of the simulations leads to the ED structure. At RMSD=22A, the dimers look like two parallel Us, still linked by the initial loop-loop interaction, but the strands of the stems (the arms of the Us) are positioned in such a manner that they can form intramolecular as well as intermolecular Watson-Crick base-pairs. This family of structure is called UU. In a second step (TMD simulations), 18 structures were picked up along the pathways generated with TMD(-1) and were constrained to move toward ED by decreasing progressively their RMSD from ED. We found that only structures from the UU family are able to easily reach ED-like conformations of the backbones without exhibiting a large constraint energy.     

HIV-1(Lai) genomic RNA: combined used of NMR and molecular dynamics simulation for studying the structure and internal dynamics of a mutated SL1 hairpin

The genome of all retroviruses consists of two identical copies of an RNA sequence associated in a non-covalent dimer. A region upstream from the splice donor (SL1) comprising a self-complementary sequence is responsible for the initiation of the dimerization. This region is able to dimerize in two conformations: a loop-loop complex or an extended duplex. Here, we solve by 2D NMR techniques the solution structure of a 23-nucleotide sequence corresponding to HIV-1 SL1(Lai) in which the mutation G12-->A12 is included to prevent dimerization. It is shown that this monomer adopts a stem-loop conformation with a seven base pairs stem and a nine nucleotide loop containing the G10 C11 A12 C13 G14 C15 sequence. The stem is well structured in an A-form duplex, while the loop is more flexible even though elements of structure are evident. We show that the structure adopted by the stem can be appreciably different from its relaxed structure when the adenines A8, A9 and A16 in the loop are mechanically constrained. This point could be important for the efficiency of the dimerization. This experimental study is complemented with a 10 ns molecular dynamics simulation in the presence of counterions and explicit water molecules. This simulation brings about information on the flexibility of the loop, such as a hinge motion between the stem and the loop and a labile lattice of hydrogen bonds in the loop. The bases of the nucleotides G10 to C15 were found outside of the loop during a part of the trajectory, which is certainly necessary to initiate the dimerization process of the genuine SL1(Lai) sequence.     

Molecular dynamics simulation for probing the flexibility of the 35 nucleotide SL1 sequence kissing complex from HIV-1Lai genomic RNA

The SL1 stem-loop located in the encapsidation domain is responsible for initiating the dimerisation of HIV-1 genomic RNA by means of a loop-loop interaction known as Kissing Complex (KC). The SL1 secondary structure has been predicted as a 35 nucleotides [K. G. Murti, M. Bondurant, and A. Tereba. J Virol 37, 411-419 (1981)] stem-loop composed of a 4 base pairs (bp) terminal duplex, a 4 nt asymmetrical internal loop, a 7 bp internal duplex, and a 9 nt apical loop. Several high resolution structures of the monomer and of KC of a 23 nt sequence containing only the internal duplex and the apical loop of SL1 are available in the literature. No experimental high resolution structure of the complete native SL1 sequence has been reported so far, either for the monomer or for KC. The asymmetrical internal loop has been described from NMR studies of different monomeric hairpin sequences, leading to divergent results, which suggests its high flexibility. In this work, we built a SL1(35) KC model which was submitted to a 31 ns molecular dynamics simulation (MD). Our results allows to describe the internal dynamics of SL1(35) KC and the differences of behavior of the different parts of the dimer. Thus, we could show the stability of the interactions between the two apical loops and of the terminal duplexes, the destabilization of the internal duplexes and the high flexibility of the asymmetrical internal loops.     

Base pairing at the stem-loop junction in the SL1 kissing complex of HIV-1 RNA: a thermodynamic study probed by molecular dynamics simulation

The thermodynamics of the opening/closure process of a GC base pair located at the stem-loop junction of the SL1 sequence from HIV-1(Lai) genomic RNA was investigated in the context of a loop-loop homodimer (or kissing complex) using molecular dynamics simulation. The free energy, enthalpy and entropy changes for the closing reaction are 0 kcal x mol(-1), -11 kcal x mol(-1) and -0.037 kcal x mol(-1) x K(-1) at 300 degrees K respectively. Furthermore it is found that the free energy change is the same for the formation of a 11 nucleotide loop closed with UG and for the formation of a 9 nucleotide loop closed with GC. Our study evidences the high flexibility of the nucleotides at the stem-loop junction explaining the weak stability of this structure.     

The SL1-SL2 (stem-loop) domain is the primary determinant for stability of the gamma retroviral genomic RNA dimer

Retroviral genomes are assembled from two sense-strand RNAs by noncovalent interactions at their 5' ends, forming a dimer. The RNA dimerization domain is a potential target for antiretroviral therapy and represents a compelling RNA folding problem. The fundamental dimerization unit for the Moloney murine sarcoma gamma retrovirus spans a 170-nucleotide minimal dimerization active sequence. In the dimer, two self-complementary sequences, PAL1 and PAL2, form intermolecular duplexes, and an SL1-SL2 (stem-loop) domain forms loop-loop base pairs, mediated by GACG tetraloops, and extensive tertiary interactions. To develop a framework for assembly of the retroviral RNA dimer, we quantified the stability of and established nucleotide resolution secondary structure models for sequence variants in which each motif was compromised. Base pairing and tertiary interactions between SL1-SL2 domains contribute a large free energy increment of -10 kcal/mol. In contrast, even though the PAL1 and PAL2 intermolecular duplexes span 10 and 16 bp in the dimer, respectively, they contribute only -2.5 kcal/mol to stability, roughly equal to a single new base pair. First, these results emphasize that the energetic costs for disrupting interactions in the monomer state nearly balance the PAL1 and PAL2 base pairing interactions that form in the dimer. Second, intermolecular duplex formation plays a biological role distinct from simply stabilizing the structure of the retroviral genomic RNA dimer.     

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.     

Direct quantification of HIV-1 RNA by capillary electrophoresis with laser-induced fluorescence

HIV-1 RNA was quantitated directly by capillary electrophoresis with laser-induced fluorescence (CE-LIF). CE-LIF was used to analyze cellular RNA and various nucleotide complexes. A fluorescently labeled DNA probe (DNA/RNA complex) in conjunction with thiazole orange intercalator was determined to have optimal stability and sensitivity for RNA analysis. Based on this observation, a hybridization method using a HIV-specific fluorescently labeled probe with analysis by CE-LIF was developed. Plasma samples from a HIV-seropositive patient were lysed to obtain RNA, hybridized with the HIV-specific probe and analyzed by CE-LIF. As little as 19 fg (1710 copies per 1 ml of starting plasma) of HIV RNA can be reliably and quantitatively detected. CE-LIF appears to be an efficient and sensitive method to quantitatively analyze RNA from a variety of sources.     

NMR structure of stem-loop SL2 of the HIV-1 psi RNA packaging signal reveals a novel A-U-A base-triple platform

The genome of the human immunodeficiency virus type-1 (HIV-1) contains a stretch of approximately 120 nucleotides known as the psi-site that is essential for RNA packaging during virus assembly. These nucleotides have been proposed to form four stem-loops (SL1-SL4) that have both independent and overlapping functions. Stem-loop SL2 is important for efficient recognition and packaging of the full-length, unspliced viral genome, and also contains the major splice-donor site (SD) for mRNA splicing. We have determined the structure of the 19-residue SL2 oligoribonucleotide by heteronuclear NMR methods. The structure is generally consistent with the most recent of two earlier secondary structure predictions, with residues G1-G2-C3-G4 and C6-U7 forming standard Watson Crick base-pairs with self-complementary residues C16-G17-C18-C19 and A12-G13, respectively. However, residue A15, which is located near the center of the stem, does not form a predicted bulge, and residues A5 and U14 do not form an expected Watson-Crick base-pair. Instead, these residues form a novel A5-U14-A15 base-triple that appears to be stabilized by hydrogen bonds from A15-H61 and -H62 to A5-N1 and U14-O2, respectively; from A5-H61 to U14-O2, and from C16-H42 to U14-O2'. A kink in the backbone allows the aromatic rings of the sequential U14-A15 residues to be approximately co-planar, adopting a stable "platform motif" that is structurally similar to the A-A (adenosine) platforms observed in the P4-P6 ribozyme domain of the Tetrahymena group I intron. Platform motifs generally function in RNA by mediating long-range interactions, and it is therefore possible that the A-U-A base-triple platform mediates long-range interactions that either stabilize the psi-RNA or facilitate splicing and/or packaging. Residue G8 of the G8-G9-U10-G11 tetraloop is stacked above the U7-A12 base-pair, and the remaining tetraloop residues are disordered and available for potential interactions with either other RNA or protein components.     

YB-1 stabilizes HIV-1 genomic RNA and enhances viral production

HIV-1 utilizes cellular factors for efficient replication. The viral RNA is different from cellular mRNAs in many aspects, and is prone to attacks by cellular RNA quality control systems. To establish effective infection, the virus has evolved multiple mechanisms to protect its RNA. Here, we show that expression of the Y-box binding protein 1 (YB-1) enhanced the production of HIV-1. Downregulation of endogenous YB-1 in producer cells decreased viral production. YB-1 increased viral protein expression by stabilizing HIV-1 RNAs. The stem loop 2 in the HIV-1 RNA packaging signal was mapped to be the YB-1-responsive element. Taken together, these results indicate that YB-1 stabilizes HIV-1 genomic RNA and thereby enhances HIV-1 gene expression and viral production.     

A new NMR solution structure of the SL1 HIV-1Lai loop-loop dimer

Dimerization of genomic RNA is directly related with the event of encapsidation and maturation of the virion. The initiating sequence of the dimerization is a short autocomplementary region in the hairpin loop SL1. We describe here a new solution structure of the RNA dimerization initiation site (DIS) of HIV-1_Lai. NMR pulsed field-gradient spin-echo techniques and multidimensional heteronuclear NMR spectroscopy indicate that this structure is formed by two hairpins linked by six Watson–Crick GC base pairs. Hinges between the stems and the loops are stabilized by intra and intermolecular interactions involving the A8, A9 and A16 adenines. The coaxial alignment of the three A-type helices present in the structure is supported by previous crystallography analysis but the A8 and A9 adenines are found in a bulged in position. These data suggest the existence of an equilibrium between bulged in and bulged out conformations in solution.     

Targeting to the HIV-1 RRE of the influenza virus NS1 protein effector domain as a potent, specific anti-HIV agent

The Rev axis of HIV autoregulation is one of two critical viral regulatory pathways required for expression of viral genomic and mRNA and for replication. Consequently it is an attractive therapeutic target. Previous studies have investigated the anti-HIV efficacy of targeting to the RRE (the viral RNA target sequence of the Rev axis) a trans-dominant negative inhibitor mutant Rev, M10. In this study we have fused a portion of the influenza virus NS1 protein (which normally inhibits polyA(+) mRNA transport and splicing) to the Rev M10 gene while deleting the NS1 poly(A) binding region. The resulting chimera demonstrates specific and enhanced inhibition of viral-RRE-containing RNA expression.     

核酸定量检测试验应用于HIV-1感染实验室诊断的探讨

目的:探讨核酸定量检测在HIV-1感染实验室诊断中的应用。方法:选取145例第四代抗原/抗体联合诊断筛查试验为阳性反应的血浆样本,分别用Western印迹和HIV-1核酸定量方法进行检测,综合对比分析2种方法检测结果。结果:Western印迹检出阳性样本120例,不确定样本17例,阴性样本8例;HIV-1核酸定量试验检出结果大于检测限样本131例,其中包括12例Western印迹不确定样本、2例Western印迹阴性样本;有3例Western印迹阳性样本用HIV-1核酸定量检测试验未能检出。结论:核酸定量检测试验对于HIV-1感染阳性样本是一种有效的实验室诊断方法;对HIV-1核酸定量检测结果为"TND"的样本,建议加做Western印迹或结合其他补充试验结果进行综合诊断。     

cis-Acting determinants of 7SL RNA packaging by HIV-1

The host noncoding RNA 7SL is highly enriched in the virions of retroviruses. We examined the regions of 7SL that mediate packaging by HIV-1. Both the Alu domain and the S domain were sufficient to mediate specific packaging when expressed separately as truncations of 7SL. However, while the Alu domain competed with endogenous 7SL for packaging in proportion to Gag, the S domain was packaged additively, implying that the Alu and S domains are packaged via separate mechanisms and that the Alu domain is packaged by the same mechanism as endogenous 7SL. Further truncations of the Alu domain or mutation of the Alu domain helix 5c region significantly reduced packaging efficiency, implicating helix 5c as critical for packaging, reinforcing the finding that 7SL packaging is highly selective, and confirming that 7SL is not passively acquired. Surprisingly, when the Alu domain was mutated so that it no longer contained a binding site for the SRP protein heterodimer SRP9/14, it was no longer packaged in a competitive manner but instead was packaged additively with endogenous 7SL. These data support a model in which 7SL RNA is packaged via interactions between Gag and a 7SL RNA structure that exists transiently at a discrete stage of SRP biogenesis. Our data further indicate that a secondary "additive" pathway exists that can result in the packaging of certain 7SL derivatives in molar excess to endogenously packaged 7SL.