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.     

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.     

The crystal structure of the dimerization initiation site of genomic HIV-1 RNA reveals an extended duplex with two adenine bulges

BACKGROUND: An important step in retroviral replication is dimerization of the genomic RNA prior to encapsidation. Dimerization is initiated by the formation of a transient 'kissing-loop complex' that is thought to be subsequently matured into an extended duplex by the nucleocapsid protein (NCp). Although chemical probing and nuclear magnetic resonance spectroscopy have provided insight into the structure of the kissing-loop structure, no structural information concerning the extended-duplex state is available so far. RESULTS: The structure of a minimal HIV-1 RNA dimerization initiation site has been solved at 2.3 A resolution in two different space groups. It reveals a 22 base pair extended duplex with two noncanonical Watson-Crick-like G-A mismatches, each adjacent to a bulged-out adenine. The structure shows significant asymmetry in deep groove width and G-A base-pair conformations. A network of eight magnesium cations was clearly identified, one being unusually chelated by the 3' phosphate of each bulge across an extremely narrowed deep major groove. CONCLUSIONS: These crystal structures represent the putative matured form of the initial kissing-loop complex. They show the ability of this self-complementary RNA hairpin loop to acquire a more stable extended duplex structure. Both bulged adenines form a striking 'base grip' that could be a recognition signal, either in cis for another viral RNA sequence, or in trans for a protein, possibly the NCp. Magnesium binding might be important to promote and stabilize the observed extrahelical conformation of these bulges.     

Role of conserved nucleotides in building the 16S rRNA binding site of E. coli ribosomal protein S8.

Ribosomal protein S8 specifically recognizes a helical and irregular region of 16S rRNA that is highly evolutionary constrained. Despite its restricted size, the precise conformation of this region remains a question of debate. Here, we used chemical probing to analyze the structural consequences of mutations in this RNA region. These data, combined with computer modelling and previously published data on protein binding were used to investigate the conformation of the RNA binding site. The experimental data confirm the model in which adenines A595, A640 and A642 bulge out in the deep groove. In addition to the already proposed non canonical U598-U641 interaction, the structure is stabilized by stacking interactions (between A595 and A640) and an array of hydrogen bonds involving bases and the sugar phosphate backbone. Mutations that alter the ability to form these interdependent interactions result in a local destabilization or reorganization. The specificity of recognition by protein S8 is provided by the irregular and distorted backbone and the two bulged adenines 640 and 642 in the deep groove. The third adenine (A595) is not a direct recognition site but must adopt a bulged position. The U598-U641 pair should not be directly in contact with the protein.     

Dimer initiation sequence of HIV-1Lai genomic RNA: NMR solution structure of the extended duplex.

The genome of all retrovirus consists of two copies of genomic RNA which are noncovalently linked near their 5' end. A sequence localized immediately upstream from the splice donor site inside the HIV-1 psi-RNA region was identified as the domain responsible for the dimerization initiation. It was shown that a kissing complex and a stable dimer are both involved in the HIV-1Lai RNA dimerization process in vitro. The NCp7 protein activates the dimerization by converting a transient loop-loop complex into a more stable dimer. The structure of this transitory loop-loop complex was recently elucidated by Mujeeb et al. In work presented here, we use NMR spectroscopy to determine the stable extended dimer structure formed from a 23 nucleotides RNA fragment, part of the 35 nucleotides SL1 sequence. By heating to 90 degrees C, then slowly cooling this sequence, we were able to show that an extended dimer is formed. We present evidence for the three dimensional structure of this dimer. NMR data yields evidence for a zipper like motif A8A9.A16 existence. This motif enables the surrounding bases to be positioned more closely and permit the G7 and C17 bases to be paired. This is different to other related sequences where only the kissing complex is observed, we suggest that the zipper like motif AA.A could be an important stabilization factor of the extended duplex.     

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.     

NMR structure of the full-length linear dimer of stem-loop-1 RNA in the HIV-1 dimer initiation site

The packaging signal of HIV-1 RNA contains a stem-loop structure, SL1, which serves as the dimerization initiation site for two identical copies of the genome and is important for packaging of the RNA genome into the budding virion and for overall infectivity. SL1 spontaneously dimerizes via a palindromic hexanucleotide sequence in its apical loop, forming a metastable kissing dimer form. Incubation with nucleocapsid protein causes this form to refold to a thermodynamically stable mature linear dimer. Here, we present an NMR structure of the latter form of the full-length SL1 sequence of the Lai HIV-1 isolate. The structure was refined using nuclear Overhauser effect and residual dipolar coupling data. The structure presents a symmetric homodimer of two RNA strands of 35 nucleotides each; it includes five stems separated by four internal loops. The central palindromic stem is surrounded by two symmetric adenine-rich 1-2 internal loops, A-bulges. All three adenines in each A-bulge are stacked inside the helix, consistent with the solution structures of shorter SL1 constructs determined previously. The outer 4-base pair stems and, proximal to them, purine-rich 1-3 internal loops, or G-bulges, are the least stable parts of the molecule. The G-bulges display high conformational variability in the refined ensemble of structures, despite the availability of many structural restraints for this region. Nevertheless, most conformations share a similar structural motif: a guanine and an adenine from opposite strands form a GA mismatch stacked on the top of the neighboring stem. The two remaining guanines are exposed, one in the minor groove and another in the major groove side of the helix, consistent with secondary structure probing data for SL1. These guanines may be recognized by the nucleocapsid protein, which binds tightly to the G-bulge in vitro.     

Solution conformation of a bulged adenosine base in an RNA duplex by relaxation matrix refinement

Bulges are common structural motifs in RNA secondary structure and are thought to play important roles in RNA-protein and RNA-drug interactions. Adenosine bases are the most commonly occurring unpaired base in double helical RNA secondary structures. The solution conformation and dynamics of a 25-nucleotide RNA duplex containing an unpaired adenosine, r(GGCAGAGUGCCGC): r(GCGGCACCUGCC) have been studied by NMR spectroscopy and MORASS iterative relaxation matrix structural refinement. The results show that the bulged adenosine residue stacks into the RNA duplex with little perturbation around the bulged region. Most of the bases in the RNA duplex adopt C(3)'-endo conformation, exhibiting the N-type sugar pucker as found in the A form helices. The sugars of the bulged residue and the 5' flanking residue to it are found to exhibit C(2)'-endo conformation. None of the residues are in syn conformation.     

The 5′ Stem-Loop From Human U4 snRNA Adopts a Novel Conformation Stabilized by G-C and G-U Base Pairs

Abstract

¹H NMR and molecular modeling studies of the 5′ stem-loop from human U4 snRNA were undertaken to determine the conformation of this stem-loop that is essential for spliceosome formation and pre-mRNA splicing. Sixteen of the 35 nucleotides of this stem-loop are in the loop region and inspection of the loop sequence revealed no decomposition into elements of secondary structure commonly found in other RNA stem-loops. An analysis of possible base pairing interactions for this stem-loop using the methods of Zuker revealed the lowest energy secondary structure for the 16 nucleotide loop consisted of four base pairs at the base of a non-canonical tetraloop (UUUA). This shorter stem-loop was joined to the nine base pair stem by two A residues on the 5′ side and a single bulged A on the 3′ side. Both stems also had bulged A residues. ¹H NMR experiments performed on solutions of the 35mer stem-loop, the stem region, and the loop region confirmed the 35mer adopted this secondary structure in solution. A 3D molecular model of this structure consistent with the NMR data was generated to assist in visualization of this novel structure.     

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.     

Thermodynamic examination of trinucleotide bulged RNA in the context of HIV-1 TAR RNA

RNA structures contain many bulges and loops that are expected to be sites for inter- and intra-molecular interactions. Nucleotides in the bulge are expected to influence the structure and recognition of RNA. The same stability is assigned to all trinucleotide bulged RNA in the current secondary structure prediction models. In this study thermal denaturation experiments were performed on four trinucleotide bulged RNA, in the context of HIV-1 TAR RNA, to determine whether the bulge sequence affects RNA stability and its divalent ion interactions. Cytosine-rich bulged RNA were more stable than uracil-rich bulged RNA in 1 M KCl. Interactions of divalent ions were more favorable with uracil-rich bulged RNA by ~2 kcal/mol over cytosine-rich bulged RNA. The UCU-TAR RNA (wild type) is stabilized by 1.7 kcal/mol in 9.5 mM Ca²⁺ as compared with 1 M KCl, whereas no additional gain in stability is measured for CCC-TAR RNA. These results have implications for base substitution experiments traditionally employed to identify metal ion binding sites. To our knowledge, this is the first systematic study to quantify the effect of small sequence changes on RNA stability upon interactions with divalent ions.     

Change of RNase P RNA function by single base mutation correlates with perturbation of metal ion binding in P4 as determined by NMR spectroscopy

The solution structures of two 27 nt RNA hairpins and their complexes with cobalt(III)-hexammine [Co(NH₃)₆³⁺] were determined by NMR spectroscopy. The RNA hairpins are variants of the P4 region from Escherichia coli RNase P RNA: a U-to-A mutant changing the identity of the bulged nucleotide, and a U-to-C, C-to-U double mutant changing only the bulge position. Structures calculated from NMR constraints show that the RNA hairpins adopt different conformations. In the U-to-C, C-to-U double mutant, the conserved bulged uridine in the P4 wild-type stem is found to be shifted in the 3′-direction by one nucleotide when compared with the wild-type structure. Co(NH₃)₆³⁺ is used as a spectroscopic probe for Mg(H₂O)₆²⁺ binding sites because both complexes have octahedral symmetry and have similar radii. Intermolecular NOE crosspeaks between Co(NH₃)₆³⁺ and RNA protons were used to locate the site of Co(NH₃)₆³⁺ binding to both RNA hairpins. The metal ion binds in the major groove near a bulge loop in both mutants, but is shifted 3′ by about one base pair in the double mutant. The change of the metal ion binding site is compared with results obtained on corresponding mutant RNase P RNA molecules as reported by Harris and co-workers (RNA, 1, 210–218).     

Modeling the dynamics of a mutated stem-loop in the SL1 domain of HIV-1Lai genomic RNA by 1H-NOESY spectra

The cross-peaks of ¹H-NOESY spectra at different time delays are compared to a mode-coupling diffusion (MCD) calculation, including the evaluation of the full ¹H relaxation matrix, in the case of a 23 nucleotide fragment of the stem-loop SL1 domain of HIV-1_Lai genomic RNA mutated in a single position. The MCD theory gives significant agreement with ¹H relaxation experiments enabling a thorough understanding of the differential local dynamics along the sequence and particularly of the dynamics of nucleotides in the stem and in the loop. The differential dynamics of this hairpin structure is important in directing the dimerization of the retroviral genome, a fundamental step in the infectious process. The demonstration of a reliable use of time dependent NOE cross-peaks, largely available from NMR solution structure determination, coupled to MCD analysis, to probe the local dynamics of biological macromolecules, is a result of general interest of this paper.     

NMR structure of the mature dimer initiation complex of HIV-1 genomic RNA

The two identical genomic RNA strands inside each HIV-1 viral particle are linked through homodimerization of an RNA stem-loop, termed SL1, near their 5' ends. SL1 first dimerizes through a palindromic sequence in its loop, forming a transient kissing-loop complex which then refolds to a mature, linear duplex. We previously reported the NMR structure of a 23-base truncate of SLI in kissing-dimer form, and here report the high-resolution structure of its linear isoform. This structure comprises three short duplex regions--derived from the central palindrome and two stem regions of each strand, respectively--separated by two bulges that each encompass three unpaired adenines flanking the palindromes. The stacking pattern of these adenines differs from that seen in the kissing-loop complex, and leads to greater colinear base stacking overall. Moreover, the mechanical distortion of the palindrome helix is reduced, and base pairs ruptured during formation of the kissing-loop complex are re-established, so that all potential Watson-Crick pairs are intact. These features together likely account for the greater thermodynamic stability of the mature dimer as compared to its kissing-loop precursor.     

Binding of 9-aminoacridine to bulged-base DNA oligomers from a frame-shift hot spot

Complexes of 9-aminoacridine and two derivatives with oligomers based on the sequence of a hot spot for frame-shift mutations, 5'dGATGGGGCAG, are investigated by proton NMR and equilibrium dialysis. Competition dialysis experiments show that the drug binds bulge-containing oligomers more strongly than regular duplexes of similar sequence and length, with one apparent strong site. A duplex containing an extra cytidine in a run of C's has the highest affinity for 9-aminoacridine among the sequences tested. An oligomer containing five consecutive G.C pairs shows cooperative drug binding, indicating that G tracts of this length may have an altered helical structure. Complexes of a regular 8-mer and a 9-mer containing a bulged guanosine are examined in detail by two-dimensional NMR techniques. 9-Aminoacridine preferentially binds at TpG sites in the 8-mer but binds primarily at the bulged guanosine in the G-bulge 9-mer. Drug-DNA NOE's in the 8-mer complex are compared with the crystal structure of 9-aminoacridine and 5-iodo-CpG [Sakore et al. (1979) J. Mol. Biol. 135, 763-785]. The NMR data suggest that the drug intercalates across the base pairs of both strands with the amino group projecting into the minor groove.     

Bulged adenosine influence on the RNA duplex conformation in solution

The RNA single bulge motif is an unpaired residue within a strand of several complementary base pairs. To gain insight into structural changes induced by the presence of the adenosine bulge on RNA duplex, the solution structures of RNA duplex containing a single adenine bulge (5'-GCAGAAGAGCG-3'/5'-CGCUCUCUGC-3') and a reference duplex with all Watson-Crick base pairs (5'-GCAGAGAGCG-3'/5'-CGCUCUCUGC-3') have been determined by NMR spectroscopy. The reference duplex structure is a regular right-handed helix with all of the attributes of an A-type helix. In the bulged duplex, single adenine bulge stacks into the helix, and the bulge region forms a well-defined structure. Both structures were analyzed by the use of calculated helical parameters. Distortions induced by the accommodation of unpaired residue into the helical structure propagate over the entire structure and are manifested as the reduced base pairs inclination and x-displacement. Intrahelical position of bulged adenine A5 is stabilized by efficient stacking with 5'-neighboring residues G4.     

Solution structure of the pseudo-5' splice site of a retroviral splicing suppressor

Control of Rous sarcoma virus RNA splicing depends in part on the interaction of U1 and U11 snRNPs with an intronic RNA element called the negative regulator of splicing (NRS). A 23mer RNA hairpin (NRS23) of the NRS directly binds U1 and U11 snRNPs. Mutations that disrupt base-pairing between the loop of NRS23 and U1 snRNA abolish its negative control of splicing. We have determined the solution structure of NRS23 using NOEs, torsion angles, and residual dipolar couplings that were extracted from multidimensional heteronuclear NMR spectra. Our structure showed that the 6-bp stem of NRS23 adopts a nearly A-form duplex conformation. The loop, which consists of 11 residues according to secondary structure probing, was in a closed conformation. U913, the first residue in the loop, was bulged out or dynamic, and loop residues G914-C923, G915-U922, and U916-A921 were base-paired. The remaining UUGU tetraloop sequence did not adopt a stable structure and appears flexible in solution. This tetraloop differs from the well-known classes of tetraloops (GNRA, CUYG, UNCG) in terms of its stability, structure, and function. Deletion of the bulged U913, which is not complementary to U1 snRNA, increased the melting temperature of the RNA hairpin. This hyperstable hairpin exhibited a significant decrease in binding to U1 snRNP. Thus, the structure of the NRS RNA, as well as its sequence, is important for interaction with U1 snRNP and for splicing suppression.     

The NMR structure of an internal loop from 23S ribosomal RNA differs from its structure in crystals of 50s ribosomal subunits

Internal loops play an important role in structure and folding of RNA and in recognition of RNA by other molecules such as proteins and ligands. An understanding of internal loops with propensities to form a particular structure will help predict RNA structure, recognition, and function. The structures of internal loops 5' 1009CUAAG1013 3'/3' 1168GAAGC1164 5' and 5' 998CUAAG1002 3'/3' 1157GAAGC1153 5' from helix 40 of the large subunit rRNA in Deinococcus radiodurans and Escherichia coli, respectively, are phylogenetically conserved, suggesting functional relevance. The energetics and NMR solution structure of the loop were determined in the duplex 5' 1GGCUAAGAC9 3'/3' 18CCGAAGCUG10 5'. The internal loop forms a different structure in solution and in the crystal structures of the ribosomal subunits. In particular, the crystal structures have a bulged out adenine at the equivalent of position A15 and a reverse Hoogsteen UA pair (trans Watson-Crick/Hoogsteen UA) at the equivalent of U4 and A14, whereas the solution structure has a single hydrogen bond UA pair (cis Watson-Crick/sugar edge A15U4) between U4 and A15 and a sheared AA pair (trans Hoogsteen/sugar edge A14A5) between A5 and A14. There is cross-strand stacking between A6 and A14 (A6/A14/A15 stacking pattern) in the NMR structure. All three structures have a sheared GA pair (trans Hoogsteen/sugar edge A6G13) at the equivalent of A6 and G13. The internal loop has contacts with ribosomal protein L20 and other parts of the RNA in the crystal structures. These contacts presumably provide the free energy to rearrange the base pairing in the loop. Evidently, molecular recognition of this internal loop involves induced fit binding, which could confer several advantages. The predicted thermodynamic stability of the loop agrees with the experimental value, even though the thermodynamic model assumes a Watson-Crick UA pair.