Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model

Molecular mechanism of ATP and RNA binding to Zika virus NS3 helicase and identification of repurposed drugs using molecular dynamics simulations

Congenital Zika virus syndrome has caused a public health emergency of international concern. So far, there are no drugs available to prevent or treat the infection caused by Zika virus. The Zika virus NS3 helicase is a potential protein target for drug discovery due to its vital role in viral genome replication. NS3 helicase unwinds the viral RNA to enable the reproduction of the viral genome by the NS5 protein. NS3 helicase has two crucial binding sites; the ATP binding site and the RNA binding site. We used molecular docking and molecular dynamics (MD) simulations to study the structural behavior of Zika virus NS3 helicase in its apo form and in the presence of ATP, single-stranded RNA, and both ATP-RNA to understand their potential implications in NS3 helicase activity. Further, we have carried out virtual screening of FDA approved drugs, followed by molecular docking to identify the ATP-competitive hit molecules as probable Zika virus NS3 helicase inhibitors. The MD simulations trajectories were analyzed using normal mode analysis and principal component analysis that reveals fluctuations in the R-loop. These findings aid in understanding the molecular mechanisms of the simultaneous binding of ATP and RNA, and guide the design and discovery of new inhibitors of the Zika virus NS3 helicase as a promising drug target to treat the Zika virus infection. Communicated by Ramaswamy H. Sarma.

RNA Binding by the KTS Splice Variants of Wilms' Tumor Suppressor Protein WT1

Wilms' tumor suppressor protein WT1 regulates the expression of multiple genes through binding of the Cys₂-His₂ zinc finger domain to promoter sites. WT1 has also been proposed to be involved in post-transcriptional regulation, by binding to RNA using the same set of zinc fingers. WT1 has two major splice variants, where the Lys-Thr-Ser (KTS) tripeptide is inserted into the linker between the third and fourth zinc fingers. To obtain insights into the mechanism by which the different WT1 splice variants recognize both DNA and RNA, we have determined the solution structure of the WT1 (-KTS) zinc finger domain in complex with a 29mer stem-loop RNA. Zinc fingers 1-3 bind in a widened major groove favored by the presence of a bulge nucleotide in the double-stranded helical stem. Fingers 2 and 3 make specific contacts with the nucleobases in a conserved AUGG sequence in the helical stem. Nuclear magnetic resonance chemical shift mapping and relaxation analysis show that fingers 1-3 of the two splice variants (-KTS and +KTS) of WT1 form similar complexes with RNA. Finger 4 of the -KTS isoform interacts weakly with the RNA loop, an interaction that is abrogated in the +KTS isoform, and both isoforms bind with similar affinity to the RNA. In contrast, finger 4 is required for high-affinity binding to DNA and insertion of KTS into the linker of fingers 3 and 4 abrogates DNA binding. While finger 1 is required for RNA binding, it is dispensable for binding to consensus DNA sites.

Single-Molecule Titration in a Protein Nanoreactor Reveals the Protonation/Deprotonation Mechanism of a C:C Mismatch in DNA

Measurement of single-molecule reactions can elucidate microscopic mechanisms that are often hidden from ensemble analysis. Herein, we report the acid-base titration of a single DNA duplex confined within the wild-type α-hemolysin (α-HL) nanopore for up to 3 h, while monitoring the ionic current through the nanopore. Modulation between two states in the current-time trace for duplexes containing the C:C mismatch in proximity to the latch constriction of α-HL is attributed to the base flipping of the C:C mismatch. As the pH is lowered, the rate for the C:C mismatch to flip from the intra-helical state to the extra-helical state ( k_intra-extra) decreases, while the rate for base flipping from the extra-helical state to the intra-helical state ( k_extra-intra) remains unchanged. Both k_intra-extra and k_extra-intra are on the order of 1 × 10^-2 s^-1 to 1 × 10^-1 s^-1 and remain stable over the time scale of the measurement (several hours). Analysis of the pH-dependent kinetics of base flipping using a hidden Markov kinetic model demonstrates that protonation/deprotonation occurs while the base pair is in the intra-helical state. We also demonstrate that the rate of protonation is limited by transport of H⁺ into the α-HL nanopore. Single-molecule kinetic isotope experiments exhibit a large kinetic isotope effect (KIE) for k_intra-extra ( k_H/ k_D ≈ 5) but a limited KIE for k_extra-intra ( k_H/ k_D ≈ 1.3), supporting our model. Our experiments correspond to the longest single-molecule measurements performed using a nanopore, and demonstrate its application in interrogating mechanisms of single-molecule reactions in confined geometries.

The human mitochondrial tRNAMet: structure/function relationship of a unique modification in the decoding of unconventional codons

Human mitochondrial mRNAs utilize the universal AUG and the unconventional isoleucine AUA codons for methionine. In contrast to translation in the cytoplasm, human mitochondria use one tRNA, hmtRNA(Met)(CAU), to read AUG and AUA codons at both the peptidyl- (P-), and aminoacyl- (A-) sites of the ribosome. The hmtRNA(Met)(CAU) has a unique post-transcriptional modification, 5-formylcytidine, at the wobble position 34 (f(5)C(34)), and a cytidine substituting for the invariant uridine at position 33 of the canonical U-turn in tRNAs. The structure of the tRNA anticodon stem and loop domain (hmtASL(Met)(CAU)), determined by NMR restrained molecular modeling, revealed how the f(5)C(34) modification facilitates the decoding of AUA at the P- and the A-sites. The f(5)C(34) defined a reduced conformational space for the nucleoside, in what appears to have restricted the conformational dynamics of the anticodon bases of the modified hmtASL(Met)(CAU). The hmtASL(Met)(CAU) exhibited a C-turn conformation that has some characteristics of the U-turn motif. Codon binding studies with both Escherichia coli and bovine mitochondrial ribosomes revealed that the f(5)C(34) facilitates AUA binding in the A-site and suggested that the modification favorably alters the ASL binding kinetics. Mitochondrial translation by many organisms, including humans, sometimes initiates with the universal isoleucine codons AUU and AUC. The f(5)C(34) enabled P-site codon binding to these normally isoleucine codons. Thus, the physicochemical properties of this one modification, f(5)C(34), expand codon recognition from the traditional AUG to the non-traditional, synonymous codons AUU and AUC as well as AUA, in the reassignment of universal codons in the mitochondria.

TbMP42 is a structure-sensitive ribonuclease that likely follows a metal ion catalysis mechanism

RNA editing in African trypanosomes is characterized by a uridylate-specific insertion and/or deletion reaction that generates functional mitochondrial transcripts. The process is catalyzed by a multi-enzyme complex, the editosome, which consists of approximately 20 proteins. While for some of the polypeptides a contribution to the editing reaction can be deduced from their domain structure, the involvement of other proteins remains elusive. TbMP42, is a component of the editosome that is characterized by two C(2)H(2)-type zinc-finger domains and a putative oligosaccharide/oligonucleotide-binding fold. Recombinant TbMP42 has been shown to possess endo/exoribonuclease activity in vitro; however, the protein lacks canonical nuclease motifs. Using a set of synthetic gRNA/pre-mRNA substrate RNAs, we demonstrate that TbMP42 acts as a topology-dependent ribonuclease that is sensitive to base stacking. We further show that the chelation of Zn(2+) cations is inhibitory to the enzyme activity and that the chemical modification of amino acids known to coordinate Zn(2+) inactivates rTbMP42. Together, the data are suggestive of a Zn(2+)-dependent metal ion catalysis mechanism for the ribonucleolytic activity of rTbMP42.

Molecular dynamics simulations of human tRNA Lys,3 UUU: the role of modified bases in mRNA recognition

Accuracy in translation of the genetic code into proteins depends upon correct tRNA–mRNA recognition in the context of the ribosome. In human tRNAUUULys,3 three modified bases are present in the anticodon stem–loop—2-methylthio-N6-threonylcarbamoyladenosine at position 37 (ms²t⁶A37), 5-methoxycarbonylmethyl-2-thiouridine at position 34 (mcm⁵s²U34) and pseudouridine (ψ) at position 39—two of which, ms²t⁶A37 and mcm⁵s²U34, are required to achieve wild-type binding activity of wild-type human tRNAUUULys,3 [C. Yarian, M. Marszalek, E. Sochacka, A. Malkiewicz, R. Guenther, A. Miskiewicz and P. F. Agris (2000) Biochemistry, 39, 13390–13395]. Molecular dynamics simulations of nine tRNA anticodon stem–loops with different combinations of nonstandard bases were performed. The wild-type simulation exhibited a canonical anticodon stair-stepped conformation. The ms²t⁶ modification at position 37 is required for maintenance of this structure and reduces solvent accessibility of U36. Ms²t⁶A37 generally hydrogen bonds across the loop and may prevent U36 from rotating into solution. A water molecule does coordinate to ψ39 most of the simulation time but weakly, as most of the residence lifetimes are <40 ps.

The naturally occurring N6-threonyl adenine in anticodon loop of Schizosaccharomyces pombe tRNAi causes formation of a unique U-turn motif

Modified nucleosides play an important role in structure and function of tRNA. We have determined the solution structure of the anticodon stem-loop (ASL) of initiator tRNA of Schizosaccharomyces pombe. The incorporation of N6-threonylcarbamoyladenosine at the position 3' to the anticodon triplet (t6A37) results in the formation of a U-turn motif and enhances stacking interactions within the loop and stem regions (i.e. between A35 and t6A37) by bulging out U36. This conformation was not observed in a crystal structure of tRNAi including the same modification in its anticodon loop, nor in the solution structure of the unmodified ASL. A t6A modification also occurs in the well studied anti-stem-loop of lys-tRNA(UUU). A comparison of this stem-loop with our structure demonstrates different effects of the modification depending on the loop sequence.

Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study

Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.

Structure and dynamics of an RNA tetraloop: a joint molecular dynamics and NMR study

Molecular dynamics simulations of the RNA tetraloop 5'-CGCUUUUGCG-3' with high melting temperature and significant conformational heterogeneity in explicit water solvent are presented and compared to NMR studies. The NMR data allow for a detailed test of the theoretical model, including the quality of the force field and the conformational sampling. Due to the conformational heterogeneity of the tetraloop, high temperature (350 K) and locally enhanced sampling simulations need to be invoked. The Amber98 force field leads to a good overall agreement with experimental data. Based on NMR data and a principal component analysis of the 350 K trajectory, the dynamic structure of the tetraloop is revealed. The principal component free energy surface exhibits four minima, which correspond to well-defined conformational structures that differ mainly by their base stacking in the loop region. No correlation between the motion of the sugar rings and the stacking dynamics of the loop bases is found.

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.

Structural effects of hypermodified nucleosides in the Escherichia coli and human tRNALys anticodon loop: the effect of nucleosides s2U, mcm5U, mcm5s2U, mnm5s2U, t6A, and ms2t6A

Previous nuclear magnetic resonance (NMR) studies of unmodified and pseudouridine39-modified tRNA(Lys) anticodon stem loops (ASLs) show that significant structural rearrangements must occur to attain a canonical anticodon loop conformation. The Escherichia coli tRNA(Lys) modifications mnm(5)s(2)U34 and t(6)A37 have indeed been shown to remodel the anticodon loop, although significant dynamic flexibility remains within the weakly stacked U35 and U36 anticodon residues. The present study examines the individual effects of mnm(5)s(2)U34, s(2)U34, t(6)A37, and Mg(2+) on tRNA(Lys) ASLs to decipher how the E. coli modifications accomplish the noncanonical to canonical structural transition. We also investigated the effects of the corresponding human tRNA(Lys,3) versions of the E. coli modifications, using NMR to analyze tRNA ASLs containing the nucleosides mcm(5)U34, mcm(5)s(2)U34, and ms(2)t(6)A37. The human wobble modification has a less dramatic loop remodeling effect, presumably because of the absence of a positive charge on the mcm(5) side chain. Nonspecific magnesium effects appear to play an important role in promoting anticodon stacking. Paradoxically, both t(6)A37 and ms(2)t(6)A37 actually decrease anticodon stacking compared to A37 by promoting U36 bulging. Rather than stack with U36, the t(6)A37 nucleotide in the free tRNAs is prepositioned to form a cross-strand stack with the first codon nucleotide as seen in the recent crystal structures of tRNA(Lys) ASLs bound to the 30S ribosomal subunit. Wobble modifications, t(6)A37, and magnesium each make unique contributions toward promoting canonical tRNA structure in the fundamentally dynamic tRNA(Lys)(UUU) anticodon.

Efficient search on energy minima for structure prediction of nucleic acid motifs

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during conformational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.

Closing loop base pairs in RNA loop-loop complexes: structural behavior, interaction energy and solvation analysis through molecular dynamics simulations

Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop-loop RNA duplexes neutralized by Mg2+ counterions in aqueous phases. Mismatches GA, CU and Watson-Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop-loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.

Methylene Blue Binding to DNA with Alternating AT Base Sequence: Minor Groove Binding is Favored over Intercalation

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5'-YpR-3' intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.

Sequence context dependence of tandem guanine:adenine mismatch conformations in RNA: a continuum solvent analysis

Guanine:adenine (G:A) mismatches and in particular tandem G:A (tG:A) mismatches are frequently observed in biological RNA molecules and can serve as sites for tertiary interaction, metal binding and protein recognition. Depending on the surrounding sequence tG:A mismatches can adopt different basepairing topologies. In the sequence context (5'-) GGAC (tandem G:A in bold) a face-to-face (imino or Watson-Crick-like) pairing is preferred whereas in the CGAG context, G and A adopt a sheared arrangement. Systematic conformational searches with a generalized Born continuum model and molecular dynamics simulations including explicit water molecules and ions have been used to generate face-to-face and sheared tG:A mismatches in both CGAG and GGAC sequence contexts. Conformations from both approaches were evaluated using the same force field and a Poisson-Boltzmann continuum solvent model. Although the substate analysis predicted the sheared arrangement to be energetically preferred in both sequence contexts, a significantly greater preference of the sheared form was found for the CGAG context. In agreement with the experimental observation, the analysis of molecular dynamics trajectories indicated a preference of the sheared form in the case of the CGAG-context and a favorization of the face-to-face form in the case of the GGAC context. The computational studies allowed to identify energetic contributions that stabilize or destabilize the face-to-face and sheared tandem mismatch topologies. The calculated nonpolar solvation and Lennard-Jones packing interaction were found to stabilize the sheared topology independent of the sequence context. Electrostatic contributions are predicted to make the most significant contribution to the sequence context dependence on the structural preference of tG:A mismatches.

The crystal structure of an alternating RNA heptamer r(GUAUACA) forming a six base-paired duplex with 3'-end adenine overhangs

The crystal structure of an alternating RNA heptamer r(GUAUACA) has been determined to 2.0 A resolution and refined to an R(work) of 17.1% and R(free) of 18.5% using 2797 reflections. The heptamer crystallized in the space group C222 with a unit cell of a = 25.74, b = 106.58, c = 30.26 A and two independent strands in the asymmetric unit. Each heptamer forms a duplex with its symmetry-related strand and each duplex contains six Watson-Crick base pairs and 3'-end adenosine overhangs. Therefore, two kinds of duplex (duplex 1 and duplex 2) are formed. Duplexes 1 stack on each other forming a pseudo-continuous column, which is typical of the RNA packing mode, while duplex 2 is typical of A-DNA packing with its termini in abutting interactions. Overhang adenine residues stack within the duplexes with C3'-endo sugar pucker and C2'-endo sugar pucker in duplexes 1 and 2, respectively. A Na+ ion in the crystal lattice is water bridged to two N1 atoms of symmetry-related A7 bases.

Site-selective reactions of imperfectly matched DNA with small chemical molecules: applications in mutation detection

The last decade has witnessed many exciting scientific publications associated with site-selective reactions of small chemical molecules with imperfectly matched DNA. Typical examples are carbodiimide, hydroxylamine, potassium permanganate, osmium tetroxide, chemical tagging probes, biotinylated, chemiluminescent and fluorescent probes, and all of them selectively react with imperfectly matched DNA. More recently, some therapeutic agents including DNA intercalating drugs and groove binders have been found to promote the in vivo repair system to recognize and repair the mismatch more effectively. The results have established a novel method for detection of mismatches. Development of new chemical reactions for detection of imperfectly matched DNA and mutations is a rapidly growing field and has attracted significant interest of scientists from both chemical and biological fields and it is the main focus of this review.

Modelling ion binding to AA platform motifs in RNA: a continuum solvent study including conformational adaptation

Binding of monovalent and divalent cations to two adenine-adenine platform structures from the Tetrahymena group I intron ribozyme has been studied using continuum solvent models based on the generalised Born and the finite-difference Poisson-Boltzmann approaches. The adenine-adenine platform RNA motif forms an experimentally characterised monovalent ion binding site important for ribozyme folding and function. Qualitative agreement between calculated and experimental ion placements and binding selectivity was obtained. The inclusion of solvation effects turned out to be important to obtain low energy structures and ion binding placements in agreement with the experiment. The calculations indicate that differences in solvation of the isolated ions contribute to the calculated ion binding preference. However, Coulomb attraction and van der Waals interactions due to ion size differences and RNA conformational adaptation also influence the calculated ion binding affinity. The calculated alkali ion binding selectivity for both platforms followed the order K(+) > Na(+) > Rb(+) > Cs(+) > Li(+) (Eisenman series VI) in the case of allowing RNA conformational relaxation during docking. With rigid RNA an Eisenman series V was obtained (K(+) > Rb(+) > Na(+) > Cs(+) > Li(+)). Systematic energy minimisation docking simulations starting from several hundred initial placements of potassium ions on the surface of platform containing RNA fragments identified a coordination geometry in agreement with the experiment as the lowest energy binding site. The approach could be helpful to identify putative ion binding sites in nucleic acid structures determined at low resolution or with experimental methods that do not allow identification of ion binding sites.

Computational methods for RNA structure determination

During the past year, major improvements have been made in methods used to solve RNA structures from crystals, find RNA patterns in sequence data and determine RNA secondary structure. Computational methods for assisting an interactive computer graphics human modeler, searching the conformational space of RNA tertiary structure, studying the dynamics of complexes involving RNA and simulating RNA catalytic activities have also been advanced.

Conformational analysis of DNA-trinucleotide-hairpin-loop structures using a continuum solvent model

A number of trinucleotide sequences in DNA can form compact and stable hairpin loops that may have significance for DNA replication and transcription. The conformational analysis of these motifs is important for an understanding of the function and design of nucleic acid structures. Extensive conformational searches have been performed on three experimentally known trinucleotide hairpin loops (AGC, AAA, and GCA) closed by a four-base-pair stem. An implicit solvation model based on the generalized Born method has been employed during energy minimization and conformational search. In addition, energy-minimized conformers were evaluated using a finite-difference Poisson-Boltzmann approach. For all three loop sequences, conformations close to experiment were found as lowest-energy structures among several thousand alternative energy minima. The inclusion of reaction-field contributions was found to be important for a realistic conformer ranking. Most generated hairpin loop structures within approximately 5 kcal x mol(-1) of the lowest-energy structure have a similar topology. Structures within approximately 10 kcal x mol(-1) could be classified into about five structural families representing distinct arrangements of loop nucleotides. Although a large number of backbone torsion angle combinations were compatible with each structural class, some specific patterns could be identified. Harmonic mode analysis was used to account for differences in conformational flexibility of low-energy sub-states. Class-specific differences in the pattern of atomic fluctuations along the sequence were observed; however, inclusion of conformational entropy contributions did not change ranking of structural classes. For an additional loop sequence (AAG) with no available experimental structure, the approach suggests a lowest-energy loop topology overall similar to the other three loop sequences but closed by a different non-canonical base-pairing scheme.

Conformational dynamics of a 5S rRNA hairpin domain containing loop D and a single nucleotide bulge

Molecular modeling and molecular dynamics have been employed to study the conformation and flexibility of a 15-nucleotide fragment of the plant 5S rRNA containing loop D and a single uridine bulge. Two different model built initial structures were used: one with the bulge localized inside the helical stem and another with the bulge pointing out from the helix. Several independent 700-ps-long trajectories in aqueous solution with Na(+) conterions were produced for each starting structure. The bulge nucleotide inside the helix stayed in two main conformations, both of which affected the geometry of the stem part opposite the bulge. When the bulge nucleotide was located outside the helix, we found high base mobility and local backbone flexibility. The dynamics of the hydrogen bond network and conformational changes from a direct to a water mediated hydrogen bond in the sheared G-A basepair in the tetraloop was described. Our results correlate with lead ion induced cleavage patterns in 5S rRNA. Sites resistant to nonspecific lead cleavage appeared in all our simulations as the most rigid fragments independent of the localization of the bulge nucleotide.

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.

Simulation of the structure and dynamics of nonhelical RNA motifs

Computer simulation methods are increasingly being used to study possible conformations and dynamics of structural motifs in RNA. Recent results of molecular dynamics simulations and continum solvent studies of RNA structures and RNA-ligand complexes show promising agreement with experimental data. Combined with the ongoing progress in the experimental characterization of RNA structure and thermodynamics, these computational approaches can help to better understand the mechanism of RNA structure formation and the binding of ligands.

Crystal structure of an adenine bulge in the RNA chain of a DNA.RNA hybrid, d(CTCCTCTTC).r(gaagagagag)

Crystal structure of a DNA.RNA hybrid, d(CTCCTCTTC).r(gaagagagag), with an adenine bulge in the polypurine RNA strand was determined at 2.3 A resolution. The structure was solved by the molecular replacement method and refined to a final R-factor of 19.9% (Rfree 22.2%). The hybrid duplex crystallized in the space group I222 with unit cell dimensions, a = 46.66 A, b = 47.61 A and c = 54.05 A, and adopts the A-form conformation. All RNA and DNA sugars are in the C3'-endo conformation, the glycosyl angles in anti conformation and the majority of the C4'-C5' torsion angles in g+ except two trans angles, in conformity with the C3'-endo rigid nucleotide hypothesis. The adenine bulge is looped out and it is also in the anti C3'-endo conformation. The bulge is involved in a base-triple (C.g)*a interaction with the end base-pair (C9.g10) in the minor groove of a symmetry-related molecule. The 2' hydroxyl group of g15 is hydrogen bonded to O2P and O5' of g17, skipping the bulged adenine a16 and stabilizing the sugar-phosphate backbone of the hybrid. The hydrogen bonding and the backbone conformation at the bulged adenine site is very similar to that found in the crystal structure of a protein-RNA complex.

Conformational deformability of RNA: a harmonic mode analysis

The harmonic mode analysis method was used to characterize the conformational deformability of regular Watson-Crick paired, mismatch- and bulge-containing RNA. Good agreement between atomic Debye-Waller factors derived from x-ray crystallography of a regular RNA oligonucleotide and calculated atomic fluctuations was obtained. Calculated helical coordinate fluctuations showed a small sequence dependence of up to approximately 30-50%. A negative correlation between motions at a given base pair step and neighboring steps was found for most helical coordinates. Only very few calculated modes contribute significantly to global motions such as bending, twisting, and stretching of the RNA molecules. With respect to a local helical description of the RNA helix our calculations suggest that RNA bending is mostly due to a periodic change in the base pair step descriptors slide and roll. The presence of single guanine:uridine or guanine:adenine mismatches had little influence on the calculated RNA flexibility. In contrast, for tandem guanine:adenine base pairs the harmonic mode approach predicts a significantly reduced conformational flexibility in the case of a sheared arrangement and slightly enhanced flexibility for a face-to-face (imino proton) pairing relative to regular RNA. The presence of a single extra adenine bulge nucleotide stacked between flanking sequences resulted in an increased local atomic mobility around the bulge site (approximately 40%) and a slightly enhanced global bending flexibility. For an adenine bulge nucleotide in a looped-out conformation a strongly enhanced bulge nucleotide mobility but no increased bending flexibility compared to regular RNA was found.

Abasic sites in duplex DNA: molecular modeling of sequence-dependent effects on conformation

Molecular modeling calculations using JUnction Minimization of Nucleic Acids (JUMNA) have been used to study sequence effects on the conformation of abasic sites within duplex DNA. We have considered lesions leading to all possible unpaired bases (X), adenine, guanine, cytosine, or thymine contained within two distinct sequence contexts, CXC and GXG. Calculations were carried out on DNA 11-mers using extensive conformational search techniques to locate the most stable abasic conformations and using Poisson-Boltzmann corrected electrostatics to account for solvation effects. The results, which are in very good agreement with available experimental data, point to strong sequence effects on both the position of the unpaired base (intra or extrahelical) and on the overall curvature induced by the abasic lesion. For CXC, unpaired purines are found to lie within the helix, while unpaired pyrimidines are either extrahelical or in equilibrium between the intra and extrahelical forms. For GXG, all unpaired bases lead to intrahelical forms, but with marked, sequence-dependent differences in induced curvature.

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.