Conformational transitions of flanking purines in HIV-1 RNA dimerization initiation site kissing complexes studied by CHARMM explicit solvent molecular dynamics

Dimerization of HIV-1 genomic RNA is initiated by kissing loop interactions at the Dimerization Initiation Site (DIS). Dynamics of purines that flank the 5' ends of the loop-loop helix in HIV-1 DIS kissing complex were explored using explicit solvent molecular dynamics (MD) simulations with the CHARMM force field. Multiple MD simulations (200 ns in total) of X-ray structures for HIV-1 DIS Subtypes A, B, and F revealed conformational variability of flanking purines. In particular, the flanking purines, which in the starting X-ray structures are bulged-out and stack in pairs, formed a consecutive stack of four bulged-out adenines at the beginning of several simulations. This conformation is seen in the crystal structure of DIS Subtype F with no interference from crystal packing, and was frequently reported in our preceding MD studies performed with the AMBER force field. However, as CHARMM simulations progressed, the four continuously stacked adenines showed conformational transitions from the bulged-out into the bulged-in geometries. Although such an arrangement has not been seen in any X-ray structure, it has been suggested by a recent NMR investigation. In CHARMM simulations, in the longer time scale, the flanking purines display the tendency to move to bulged-in conformations. This is in contrast with the AMBER simulations, which indicate a modest prevalence for bulged-out flanking base positions in line with the X-ray data. The simulations also suggest that the intermolecular stacking between purines from the opposite hairpins can additionally stabilize the kissing complex.     

Polymorphism of bulged-out residues in HIV-1 RNA DIS kissing complex and structure comparison with solution studies

All retroviruses encapsidate their genome as a dimer of homologous single-stranded RNAs. The dimerization initiation site (DIS) of human immunodeficiency virus type 1 (HIV-1) is located in the 5'-untranslated region of the viral genome and consists of a hairpin with a 6 nt self-complementary loop sequence. Genomic RNA dimerization, a crucial step for virion infectivity, is promoted by the formation of a loop-loop complex (or kissing complex) between two DIS hairpins. Crystal structures for the subtypes A, B and F of the HIV-1 DIS kissing complex have now been solved at 2.3 A, 1.9 A and 1.6 A, respectively. They revealed a polymorphism of bulged-out residues showing clearly that their conformation is not a mere consequence of crystal packing. They also provide more insights into ion binding, hydration, and RNA conformation and flexibility. In particular, we observed the binding of spermine to the loop-loop helix, which displaced a magnesium cation important for subtype A DIS dimerization. The excellent agreement between X-ray structures and the results of chemical probing and interference data on larger viral RNA fragments shows that the crystal structures are relevant for the DIS kissing complex present in solution and in viral particles. Accordingly, these structures will be helpful for designing new drugs derived from aminoglycoside antibiotics and targeted against the RNA dimerization step of the viral life-cycle.     

Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets

Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing–loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem–loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs. The core region of the kissing complex forms a cation-binding pocket with highly negative electrostatic potential. The pocket shows nanosecond-scale breathing motions coupled with oscillations of the whole molecule. Additional simulations were carried out for 6 bp kissing complexes of the DIS HIV-1 subtypes A and B. The simulated structures agree well with the X-ray data. The subtype B forms a novel four-base stack of bulged-out adenines. Both 6 bp kissing complexes have extended cation-binding pockets in their central parts. While the pocket of subtype A interacts with two hexacoordinated Mg²⁺ ions and one sodium ion, pocket of subtype B is filled with a string of three delocalized Na⁺ ions with residency times of individual cations 1–2 ns. The 6 bp complexes show breathing motions of the cation-binding pockets and loop major grooves.     

Stabilities of HIV-1 DIS type RNA loop-loop interactions in vitro and in vivo

RNA loop–loop interactions are a prevalent motif in the formation of tertiary structure and are well suited to trigger molecular recognition between RNA molecules. We determined the stabilities of several loop–loop interactions with a constant 6 bp core sequence and varying unpaired flanking nucleotides and found that the flanking bases have a strong influence on the stability and ion dependence of the kissing complex. In general, the stabilities determined in 1 M Na⁺ are equivalent to those in the presence of near physiological Mg²⁺ concentrations. Therefore we further tested whether the stabilities determined in vitro and within yeast cells correlate, using a recently developed yeast RNA-hybrid system. For the majority of the loop types analyzed here, the melting temperatures determined in vitro are in good agreement with the relative β-galactosidase activity in yeast cells, showing that data derived from in vitro measurements reflect in vivo properties. The most stable interactions are the naturally occurring HIV-1 DIS MAL and LAI derived loops with the motif (5′ A^A/_GN₆A 3′), emphasizing the crucial role of stable kissing complexes in HIV genome dimerization.     

MD studies of the DIS/DIS kissing complex solution and x-ray structures

As in all retroviruses, human immunodeficiency virus (HIV) genomic RNA is packaged into virions as a dimer. The two copies of the genome are noncovalently linked by their 5'-ends in the dimerization initiating site (DIS), which folds as a hairpin containing an apical autocomplementary sequence. In vitro, DIS is able to dimerize in two conformations: a kissing complex and an extended dimer. Both conformations have been resolved by NMR and x-ray diffraction. Here, we report molecular dynamics (MD) studies of the two available structures for the DIS/DIS kissing complex in aqueous solution and in the presence of sodium counterions. The two structures behave in two different manners. On one hand, the NMR structure displays a very stable behavior, and the simulated structure remains very close to the starting structure. On the other hand, the structure issued from crystallography displays a more dynamic behavior, in which residues A8 and A9 are seen in a new and surprising bulge-in conformation. The transition from the bulge-out to the bulge-in conformation is analyzed, and a new and simple dimerization process is proposed.     

Crystal structures of coaxially stacked kissing complexes of the HIV-1 RNA dimerization initiation site

We describe the crystal structures of the RNA dimerization initiation sites (DIS) of HIV-1 subtypes A and B. Both molecules adopt a hairpin conformation, with loop sequences consisting of 272-AGGUGCACA-280 and 272-AAGCGCGCA-280, respectively. This includes a six-base self-complementary stretch (underlined) that allows homodimerization through the formation of a loop-loop, or 'kissing-loop', complex. The DISs for the two sequences have identical conformations, and the two interacting hairpins show a perfect coaxial alignment. The conserved purines, A272 and R273, are stacked in a bulged-out conformation and symmetrically join the upward and downward strands of each hairpin by crossing the helix major groove. There is good agreement between these structures and previous results from chemical probing in solution, as well as with differences in magnesium dependence for dimerization. The overall shape of the kissing-loop complex is very similar to that of the previously determined subtype A DIS duplex form. Unexpectedly, the purine R273 is the only base seen at a different position and is responsible for the difference in topology between the two forms. We propose that the transition from kissing-loop duplex could occur by a recombination mechanism based on symmetrical chain cleavage at R273 of each hairpin and subsequent cross-religation, rather than by base-pair melting and subsequent reannealing.     

Solution studies of the dimerization initiation site of HIV-1 genomic RNA.

Dimerization of HIV-1 genomic RNA is an essential step of the viral cycle, initiated at a conserved stem-loop structure which forms a 'kissing complex' involving loop-loop interactions (dimerization initiation site, DIS). A 19mer RNA oligonucleotide (DIS-19) has been synthesized which forms a stable symmetrical dimer in solution at millimolar concentrations. Dimerization does not depend on addition of Mg2+. RNA ligation experiments unambiguously indicate that the formed dimer is a stable kissing complex under the NMR experimental conditions.1H NMR resonance assignments were obtained for DIS-19 at 290 K and pH 6.5. Analysis of the pattern of NOE connectivities reveals that the helix formed by loop-loop base pairing is stacked onto the two terminal stems. Non-canonical base pairs between two essential and conserved adenines are found at the junctions between the two intramolecular and the single intramolecular helices.     

A Short Sequence Motif in the 5′ Leader of the HIV-1 Genome Modulates Extended RNA Dimer Formation and Virus Replication

The 5′ leader of the HIV-1 RNA genome encodes signals that control various steps in the replication cycle, including the dimerization initiation signal (DIS) that triggers RNA dimerization. The DIS folds a hairpin structure with a palindromic sequence in the loop that allows RNA dimerization via intermolecular kissing loop (KL) base pairing. The KL dimer can be stabilized by including the DIS stem nucleotides in the intermolecular base pairing, forming an extended dimer (ED). The role of the ED RNA dimer in HIV-1 replication has hardly been addressed because of technical challenges. We analyzed a set of leader mutants with a stabilized DIS hairpin for in vitro RNA dimerization and virus replication in T cells. In agreement with previous observations, DIS hairpin stability modulated KL and ED dimerization. An unexpected previous finding was that mutation of three nucleotides immediately upstream of the DIS hairpin significantly reduced in vitro ED formation. In this study, we tested such mutants in vivo for the importance of the ED in HIV-1 biology. Mutants with a stabilized DIS hairpin replicated less efficiently than WT HIV-1. This defect was most severe when the upstream sequence motif was altered. Virus evolution experiments with the defective mutants yielded fast replicating HIV-1 variants with second site mutations that (partially) restored the WT hairpin stability. Characterization of the mutant and revertant RNA molecules and the corresponding viruses confirmed the correlation between in vitro ED RNA dimer formation and efficient virus replication, thus indicating that the ED structure is important for HIV-1 replication.     

Molecular modeling and dynamics studies of HIV-1 kissing loop structures

Recognition of an RNA loop by another RNA loop is involved in several biological functions. The dimerization of two copies of the HIV-1 genomic RNA is thought to be involved in several steps of the retroviral life cycle. It has been shown that the dimerization of the two HIV-1 RNA genomes is initiated by the so called kissing loop. The 9nt kissing loop consists of a palindromic 6nt sequence that forms Watson-Crick base-pairs at the kissing site in HIV-1. We report the results of our molecular modeling and dynamics studies on two major subtype isolates (MAL and LAI) of HIV-1 kissing loop structures. From our modeling studies, we conclude that the conformation of the loop in the monomer might be closer to the A-RNA-like conformation in order to form an initial kissing structure. This is achieved by the stacking interactions of the bases at the 3' end of the loop and by the intramolecular tertiary interactions of a single linker nucleotide. We discuss the effect of the loop size and the structural limitations on the formation of kissing loop structures. Also, we propose a possible mechanism to convert the kissing loop structure to a stable extended duplex structure without unwinding the stems.     

Analysis of the stability of looped-out and stacked-in conformations of an adenine bulge in DNA using a continuum model for solvent and ions

A combination of conformational search, energy minimization, and energetic evaluation using a continuum solvent treatment has been employed to study the stability of various conformations of the DNA fragment d(CGCAGAA)/d(TTCGCG) containing a single adenine bulge. The extra-helical (looped-out) bulge conformation derived from a published x-ray structure and intra-helical (stacked bulge base) model structures partially based on nuclear magnetic resonance (NMR) data were used as start structures for the conformational search. Solvent-dependent contributions to the stability of the conformations were calculated from the solvent exposed molecular surface area and by using the finite difference Poisson-Boltzmann approach. Three classes (I-III) of bulge conformations with calculated low energies can be distinguished. The lowest-energy conformations were found in class I, corresponding to structures with the bulge base stacked between flanking helices, and class II, composed of structures forming a triplet of the bulge base and a flanking base pair. All extra-helical bulge structures, forming class III, were found to be less stable compared with the lowest energy structures of class I and II. The results are consistent with NMR data on an adenine bulge in the same sequence context indicating an intra-helical or triplet bulge conformation in solution. Although the total energies and total electrostatic energies of the low-energy conformations show only relatively modest variations, the energetic contributions to the stability were found to vary significantly among the classes of bulge structures. All intra-helical bulge structures are stabilized by a more favorable Coulomb charge-charge interaction but destabilized by a larger electrostatic reaction field contribution compared with all extra-helical and most triplet bulge structures. Van der Waals packing interactions and nonpolar surface-area-dependent contributions appear to favor triplet class II structures and to a lesser degree also the intra-helical stacked bulge conformations. The large conformational variation found for class III conformers might add a favorable entropic contribution to the stability of the extra-helical bulge form.     

Mechanism of nucleocapsid protein catalyzed structural isomerization of the dimerization initiation site of HIV-1

Dimerization of two homologous strands of genomic RNA is an essential feature of retroviral replication. In the human immunodeficiency virus type 1 (HIV-1), a conserved stem-loop sequence, the dimerization initiation site (DIS), has been identified as the domain primarily responsible for initiation of this aspect of viral assembly. The DIS loop contains an autocomplementary hexanucleotide sequence flanked by highly conserved 5' and 3' purines and can form a homodimer through a loop-loop kissing interaction. In a structural rearrangement activated by the HIV-1 nucleocapsid protein (NCp7) and considered to be associated with viral particle maturation, the DIS dimer converts from an intermediate kissing to an extended duplex isoform. Using 2-aminopurine (2-AP) labeled sequences derived from the DIS(Mal) variant and fluorescence methods, the two DIS dimer isoforms have been unambiguously distinguished, allowing a detailed examination of the kinetics of this RNA structural isomerization and a characterization of the role of NCp7 in the reaction. In the presence of divalent cations, the DIS kissing dimer is found to be kinetically trapped and converts to the extended duplex isoform only upon addition of NCp7. NCp7 is demonstrated to act catalytically in inducing the structural isomerization by accelerating the rate of strand exchange between the two hairpin stem helices, without disruption of the loop-loop helix. Observation of an apparent maximum conversion rate for NCp7-activated DIS isomerization, however, requires protein concentrations in excess of the 2:1 stoichiometry estimated for high-affinity NCp7 binding to the DIS kissing dimer, indicating that transient interactions with additional NCp7(s) may be required for catalysis.     

HIV-1 DIS stem loop forms an obligatory bent kissing intermediate in the dimerization pathway

The HIV-1 dimerization initiation sequence (DIS) is a conserved palindrome in the apical loop of a conserved hairpin motif in the 5′-untranslated region of its RNA genome. DIS hairpin plays an important role in genome dimerization by forming a ‘kissing complex’ between two complementary hairpins. Understanding the kinetics of this interaction is key to exploiting DIS as a possible human immunodeficiency virus (HIV) drug target. Here, we present a single-molecule Förster resonance energy transfer (smFRET) study of the dimerization reaction kinetics. Our data show the real-time formation and dissociation dynamics of individual kissing complexes, as well as the formation of the mature extended duplex complex that is ultimately required for virion packaging. Interestingly, the single-molecule trajectories reveal the presence of a previously unobserved bent intermediate required for extended duplex formation. The universally conserved A272 is essential for the formation of this intermediate, which is stabilized by Mg²⁺, but not by K⁺ cations. We propose a 3D model of a possible bent intermediate and a minimal dimerization pathway consisting of three steps with two obligatory intermediates (kissing complex and bent intermediate) and driven by Mg²⁺ ions.     

The dimer initiation site hairpin mediates dimerization of the human immunodeficiency virus, type 2 RNA genome

The untranslated leader of retroviral RNA genomes encodes multiple structural signals that are critical for virus replication. In the human immunodeficiency virus, type 1 (HIV-1) leader, a hairpin structure with a palindrome-containing loop is termed the dimer initiation site (DIS), because it triggers in vitro RNA dimerization through base pairing of the loop-exposed palindromes (kissing loops). Controversy remains regarding the region responsible for HIV-2 RNA dimerization. Different studies have suggested the involvement of the transactivation region, the primer binding site, and a hairpin structure that is the equivalent of the HIV-1 DIS hairpin. We have performed a detailed mutational analysis of the HIV-2 leader RNA, and we also used antisense oligonucleotides to probe the regions involved in dimerization. Our results unequivocally demonstrate that the DIS hairpin is the main determinant for HIV-2 RNA dimerization. The 6-mer palindrome sequence in the DIS loop is essential for dimer formation. Although the sequence can be replaced by other 6-mer palindromes, motifs that form more than two A/U base pairs do not dimerize efficiently. The inability to form stable kissing-loop complexes precludes formation of dimers with more extended base pairing. Structure probing of the DIS hairpin in the context of the complete HIV-2 leader RNA suggests a 5-base pair elongation of the DIS stem as it is proposed in current RNA secondary structure models. This structure is supported by phylogenetic analysis of leader RNA sequences from different viral isolates, indicating that RNA genome dimerization occurs by a similar mechanism for all members of the human and simian immunodeficiency viruses.     

Solution structure of the SL1 RNA of the M1 double-stranded RNA virus of Saccharomyces cerevisiae

The 20-nucleotide SL1 VBS RNA, 5'-GGAGACGC[GAUUC]GCGCUCC (bulged A underlined and loop bases in brackets), plays a crucial role in viral particle binding to the plus strand and packaging of the RNA. Its structure was determined by NMR spectroscopy. Structure calculations gave a precisely defined structure, with an average pairwise root mean square deviation (RMSD) of 1.28 A for the entire molecule, 0.57 A for the loop region (C8-G14), and 0.46 A for the bulge region (G4-G7, C15-C17). Base stacking continues for three nucleotides on the 5' side of the loop. The final structure contains a single hydrogen bond involving the guanine imino proton and the carbonyl O(2) of the cytosine between the nucleotides on the 5' and 3' ends of the loop, although they do not form a Watson-Crick base pair. All three pyrimidine bases in the loop point toward the major groove, which implies that Cap-Pol protein may recognize the major groove of the SL1 loop region. The bulged A5 residue is stacked in the stem, but nuclear Overhauser enhancements (NOEs) suggest that A5 spends part of the time in the bulged-out conformation. The rigid conformation of the upper stem and loop regions may allow the SL1 VBS RNA to interact with Cap-Pol protein without drastically changing its own conformation.     

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.     

2-Aminopurine Fluorescence: Discrimination Between Specific and Unspecific Ligand Binding to the Kissing-Loop Dimer of the HIV-1 RNA

Abstract

The fluorescent 2-aminopurine probe (2-AP) incorporated into the loop of 23-mer RNA hairpin of HIV-1 genome dimerization initiation site (DIS) was used for discrimination of specific and unspecific binding of paromomycin and spermine to the kissing loop dimer (KD) formed in solution. While both ligands stabilized the KD RNA structure, only paromomycin binding resulted in significant increase of 2-AP fluorescence. These observations suggest that the 2-AP fluorescent RNA construct might be useful for selecting ligands specifically binding the HIV-1 kissing loop RNA dimer.     

Dissecting the protein-RNA and RNA-RNA interactions in the nucleocapsid-mediated dimerization and isomerization of HIV-1 stemloop 1

The specific binding of HIV-1 nucleocapsid protein (NC) to the different forms assumed in vitro by the stemloop 1 (Lai variant) of the genome's packaging signal has been investigated using electrospray ionization-Fourier transform mass spectrometry (ESI-FTMS). The simultaneous observation of protein-RNA and RNA-RNA interactions in solution has provided direct information about the role of NC in the two-step model of RNA dimerization and isomerization. In particular, two distinct binding sites have been identified on the monomeric stemloop structure, corresponding to the apical loop and stem-bulge motifs. These sites share similar binding affinities that are intermediate between those of stemloop 3 (SL3) and the putative stemloop 4 (SL4) of the packaging signal. Binding to the apical loop, which contains the dimerization initiation site (DIS), competes directly with the annealing of self-complementary sequences to form a metastable kissing-loop (KL) dimer. In contrast, binding to the stem-bulge affects indirectly the monomer-dimer equilibrium by promoting the rearrangement of KL into the more stable extended duplex (ED) conformer. This process is mediated by the duplex-melting activity of NC, which destabilizes the intramolecular base-pairs surrounding the KL stem-bulges and enables their exchange to form the inter-strand pairs that define the ED structure. In this conformer, high-affinity binding takes place at stem-bulge sites that are identical to those present in the monomeric and KL forms. In this case, however, the NC-induced "breathing" does not result in dissociation of the double-stranded structure because of the large number of intermolecular base-pairs. The different binding modes manifested by conformer-specific mutants have shown that NC can also provide low affinity interactions with the bulged-out adenine bases flanking the DIS region of the ED conformer, thus supporting the hypothesis that these exposed nucleotides may constitute "base-grips" for protein contacts during the late stages of the viral lifecycle.     

2-aminopurine fluorescence: discrimination between specific and unspecific ligand binding to the kissing-loop dimer of the HIV-1 RNA

The fluorescent 2-aminopurine probe (2-AP) incorporated into the loop of 23-mer RNA hairpin of HIV-1 genome dimerization initiation site (DIS) was used for discrimination of specific and unspecific binding of paromomycin and spermine to the kissing loop dimer (KD) formed in solution. While both ligands stabilized the KD RNA structure, only paromomycin binding resulted in significant increase of 2-AP fluorescence. These observations suggest that the 2-AP fluorescent RNA construct might be useful for selecting ligands specifically binding the HIV-1 kissing loop RNA dimer.