NCIR: a database of non-canonical interactions in known RNA structures

The secondary and tertiary structure of an RNA molecule typically includes a number of non-canonical base–base interactions. The known occurrences of these interactions are tabulated in the NCIR database, which can be accessed from http://prion.bchs.uh.edu/bp_type/. The number of examples is now over 1400, which is an increase of >700% since the database was first published. This dramatic increase reflects the addition of data from the recently published crystal structures of the 50S (2.4 Å) and 30S (3.0 Å) ribosomal subunits. In addition, non-canonical interactions observed in published crystal and NMR structures of tRNAs, group I introns, ribozymes, RNA aptamers and synthetic oligonucleotides are included. Properties associated with these interactions, such as sequence context, sugar pucker conformation, glycosidic angle conformation, melting temperature, chemical shift and free energy, are also reported when available. Out of the 29 anticipated pairs with at least two hydrogen bonds, 28 have been observed to date. In addition, several novel examples, not generally predicted, have also been encountered, bringing the total of such pairs to 36. Added to this list are a variety of single, bifurcated, triple and quadruple interactions. The most common non-canonical pairs are the sheared GA, GA imino, AU reverse Hoogsteen, and the GU and AC wobble pairs. The most frequent triple interaction connects N3 of an A with the amino of a G that is also involved in a standard Watson–Crick pair.     

Structural basis of DNA flexibility

It has been recently shown by us, on the basis of crystal structure database that the flexibility of B-DNA double helices depends significantly on their base sequence. Our model building studies further indicated that the existence of bifurcated cross-strand hydrogen bonds between successive base pairs is possibly the main factor behind the sequence directed DNA flexibility. These cross-strand hydrogen bonds are, of course, weaker than the usual Watson-Crick hydrogen bonds and their bond geometry is characterized by relatively larger bond lengths and smaller bond angles. We have tried to improve our model structures by incorporating non-planarity of the amino groups in DNA bases due to the presence of lone pair electrons at the nitrogen atoms. Energy minimization studies have been carried out by using different quantum chemical methods, whereby it is found that in all cases of N-H....O type cross-strand hydrogen bonds, the bond geometry improves significantly. In the cases of N-H....N type hydrogen bonds, however, no such consistent improvements can be noticed. Perhaps the true picture would emerge only if all the other interactions present in the DNA macromolecule could be appropriately taken into account.     

Binding of the antitumor drug nogalamycin and its derivatives to DNA: structural comparison

The three-dimensional molecular structures of the complexes between a novel antitumor drug nogalamycin and its derivative U-58872 with a modified DNA hexamer d[m5CGT(pS)Am5CG] have been determined at 1.7- and 1.8-A resolution, respectively, by X-ray diffraction analyses. Both structures (in space group P6(1)) have been refined with constrained refinement procedure to final R factors of 0.208 (3386 reflections) and 0.196 (2143 reflections). In both complexes, two nogalamycins bind to the DNA hexamer double helix in a 2:1 ratio with the elongated aglycon chromophore intercalated between the CpG steps at both ends of the helix. The aglycon chromophore spans across the GC Watson-Crick base pairs with its nogalose lying in the minor groove and the aminoglucose lying in the major groove of the distorted B-DNA double helix. Most of the sugars remain in the C2'-endo pucker family, except three deoxycytidine residues (terminal C1, C7, and internal C5). All nucleotides are in the anti conformation. Specific hydrogen bonds are found in the complex between the drug and guanine-cytosine bases in both grooves of the helix. One hydroxyl group of the aminoglucose donates a hydrogen bond to the N7 of guanine, while the other receives a hydrogen bond from the N4 amino group of cytosine. The orientation of these two hydrogen bonds suggests that nogalamycin prefers a GC base pair with its aglycon chromophore intercalating at the 5'-side of a guanine (between NpG), or at the 3'-side of a cytosine (between CpN) with the sugars pointing toward the GC base pair. The binding of nogalamycin to DNA requires that the base pairs in DNA open up transiently to allow the bulky sugars to go through, suggesting that nogalamycin prefers GC sequences embedded in a stretch of AT sequences.     

Simulation of interactions in the co-planar nucleic acid base pairs using atom-atom potential functions

Calculations of the energy of nucleic acid base interactions as a function of parameters determining mutual position of two bases in a plane have been performed. Atom-atom potential functions used include terms proportional to the first (electrostatic), sixth (or tenth for the atoms of hydrogen bond) and 12th power of interatomic distance. The calculations have shown the existence of 27 energy minima which correspond to the formation of co-planar pairs with two (or three for G : C pair) almost linear N--H...O and N--H...N hydrogen bonds. The positions of nitrogen bases bound by two hydrogen bonds in every crystal of nucleic acid components, in the complexes of polynucleotides and in tRNA are near to the positions in one of these minima. In addition for every pair there exist energy minima which correspond to the formation of one N--H...O or N--H...N and one C--H...O or C--H...N hydrogen bond. Energy behavior near minima have been investigated. The results of our calculations are in agreement with experimental data and with the calculations which employ quantum mechanical results.     

Non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis

Non-canonical base pairs, mostly present in the RNA, often play a prominent role towards maintaining their structural diversity. Higher order structures like base triples are also important in defining and stabilizing the tertiary folded structure of RNA. We have developed a new program BPFIND to analyze different types of canonical and non-canonical base pairs and base triples involving at least two direct hydrogen bonds formed between polar atoms of the bases or sugar O2' only. We considered 104 possible types of base pairs, out of which examples of 87 base pair types are found to occur in the available RNA crystal structures. Analysis indicates that approximately 32.7% base pairs in the functional RNA structures are non-canonical, which include different types of GA and GU Wobble base pairs apart from a wide range of base pair possibilities. We further noticed that more than 10.4% of these base pairs are involved in triplet formation, most of which play important role in maintaining long-range tertiary contacts in the three-dimensional folded structure of RNA. Apart from detection, the program also gives a quantitative estimate of the conformational deformation of detected base pairs in comparison to an ideal planar base pair. This helps us to gain insight into the extent of their structural variations and thus assists in understanding their specific role towards structural and functional diversity.     

Probing the role of hydrogen bonds in the stability of base pairs in double-helical DNA

Aromatic stacking and hydrogen bonding between nucleobases are two of the key interactions responsible for stabilization of DNA double-helical structures. The present work aims at defining the specific contributions of these interactions to the stability of individual base pairs in DNA. The two DNA double helices investigated are formed, respectively, by the palindromic base sequences 5'-dCCAACGTTGG-3' and 5'-dCGCAGATCTGCG-3'. The strength of the N==H...N inter-base hydrogen bond in each base pair is characterized from the measurement of the protium-deuterium fractionation factor of the corresponding imino proton using NMR spectroscopy. The structural stability of each base pair is evaluated from the exchange rate of the imino proton, measured by NMR. The results reveal that the fractionation factors of the imino protons in the two DNA double helices investigated fall within a narrow range of values, between 0.92 and 1.0. In contrast, the free energies of structural stabilization for individual base pairs span 3.5 kcal/mol, from 5.2 to 8.7 kcal/mol (at 15 degrees C). These findings indicate that, in the two DNA double helices investigated, the strength of N==H...N inter-base hydrogen bonds does not change significantly depending on the nature or the sequence context of the base pair. Hence, the variations in structural stability detected by proton exchange do not involve changes in the strength of inter-base hydrogen bonds. Instead, the results suggest that the energetic identity of a base pair is determined by the number of inter-base hydrogen bonds, and by the stacking interactions with neighboring base pairs.     

NMR studies of the stable mismatch purine-thymine in the self-complementary d(CGPuAATTTCG) duplex in solution

One- and two-dimensional nuclear Overhauser effect experiments demonstrate that a single hydrogen bond between a T imino proton and purine N3 is sufficient to hold the base pair dPu.dT in d(CGPuAATTTCG) by a Watson-Crick fashion rather than a Hoogsteen type. In addition, the dPu.dT base pair is well stacked with neighboring base pairs. The spin-lattice relaxation measurements at 30 and 35 degrees C of two decamers, d(CGPuAATTTCG) and d(CGAAATTTCG), reveal that the elimination of two single hydrogen bonds of dA.dT base pairs (due to the substitution of adenine for purine) in the sequence results in an increase in the overall imino proton exchange rate from 7 to 36 s-1 at the site of mismatch.     

Nuclear magnetic resonance spectroscopy and molecular modeling reveal that different hydrogen bonding patterns are possible for G.U pairs: one hydrogen bond for each G.U pair in r(GGCGUGCC)(2) and two for each G.U pair in r(GAGUGCUC)(2)

G.U pairs occur frequently and have many important biological functions. The stability of symmetric tandem G.U motifs depends both on the adjacent Watson-Crick base pairs, e.g., 5'G > 5'C, and the sequence of the G.U pairs, i.e., 5'-UG-3' > 5'-GU-3', where an underline represents a nucleotide in a G.U pair [Wu, M., McDowell, J. A., and Turner, D. H. (1995) Biochemistry 34, 3204-3211]. In particular, at 37 degrees C, the motif 5'-CGUG-3' is less stable by approximately 3 kcal/mol compared with other symmetric tandem G.U motifs with G-C as adjacent pairs: 5'-GGUC-3', 5'-GUGC-3', and 5'-CUGG-3'. The solution structures of r(GAGUGCUC)(2) and r(GGCGUGCC)(2) duplexes have been determined by NMR and restrained simulated annealing. The global geometry of both duplexes is close to A-form, with some distortions localized in the tandem G.U pair region. The striking discovery is that in r(GGCGUGCC)(2) each G.U pair apparently has only one hydrogen bond instead of the two expected for a canonical wobble pair. In the one-hydrogen-bond model, the distance between GO6 and UH3 is too far to form a hydrogen bond. In addition, the temperature dependence of the imino proton resonances is also consistent with the different number of hydrogen bonds in the G.U pair. To test the NMR models, U or G in various G.U pairs were individually replaced by N3-methyluridine or isoguanosine, respectively, thus eliminating the possibility of hydrogen bonding between GO6 and UH3. The results of thermal melting studies on duplexes with these substitutions support the NMR models.     

Structural features and hydration of d(C-G-C-G-A-A-T-T-A-G-C-G); a double helix containing two G.A mispairs

Single crystal X-ray diffraction techniques have been used to characterise the molecular structure of the title compound to 2.5A resolution. The structure consists of ten standard Watson-Crick base pairs and two G.A mismatched base pairs. The purine-purine mismatches have guanine in the usual anti orientation with respect to the sugar and adenine in syn orientation. There are two hydrogen bonds formed between the mismatch bases, N-1 and O-6 of guanine with N-7 and N-6 of adenine respectively. The bulky purine-purine mismatches are accommodated with minor perturbation of the sugar-phosphate backbone. There is a slight improvement in base pair overlap at the mismatch sites. Details of the backbone conformation, base stacking interactions and hydration are presented and compared with those of the parent compound d(C-G-C-G-A-A-T-T-C-G-C-G).     

Deoxyguanosine-deoxyadenosine pairing in the d(C-G-A-G-A-A-T-T-C-G-C-G) duplex: conformation and dynamics at and adjacent to the dG X dA mismatch site

Nuclear magnetic resonance (NMR) has been used to monitor the conformation and dynamics of the d-(C1-G2-A3-G4-A5-A6-T6-T5-C4-G3-C2-G1) self-complementary dodecanucleotide (henceforth called 12-mer GA) that contains a dG X dA purine-purine mismatch at position 3 in the sequence. These results are compared with the corresponding d(C-G-C-G-A-A-T-T-C-G-C-G) dodecamer duplex (henceforth called 12-mer) containing standard Watson-Crick base pairs at position 3 [Patel, D.J., Kozlowski, S.A., Marky, L.A., Broka, C., Rice, J.A., Itakura, K., & Breslauer, K.J. (1982) Biochemistry 21, 428-436]. The dG X dA interaction at position 3 was monitored at the guanosine exchangeable H-1 and nonexchangeable H-8 protons and the nonexchangeable adenosine H-2 proton. We demonstrate base-pair formation between anti orientations of the guanosine and adenosine rings on the basis of nuclear Overhauser effects (NOE) observed between the H-2 proton of adenosine 3 and the imino protons of guanosine 3 (intra base pair) and guanosines 2 and 4 (inter base pair). The dG(anti) X dA(anti) pairing should result in hydrogen-bond formation between the guanosine imino H-1 and carbonyl O-6 groups and the adenosine N-1 and NH2-6 groups, respectively. The base pairing on either side of the dG X dA pair remains intact at low temperature, but these dG X dC pairs at positions 2 and 4 are kinetically destabilized in the 12-mer GA compared to the 12-mer duplex. We have estimated the hydrogen exchange kinetics at positions 4-6 from saturation-recovery measurements on the imino protons of the 12-mer GA duplex between 5 and 40 degrees C. The measured activation energies for imino proton exchange in the 12-mer GA are larger by a factor of approximately 2 compared to the corresponding values in the 12-mer duplex. This implies that hydrogen exchange in the 12-mer GA duplex results from a cooperative transition involving exchange of several base pairs as was previously reported for the 12-mer containing a G X T wobble pair at position 3 [Pardi, A., Morden, K.M., Patel, D.J., & Tinoco, I., Jr. (1982) Biochemistry 21, 6567-6574]. We have assigned the nonexchangeable base protons by intra and inter base pair NOE experiments and monitored these assigned markers through the 12-mer GA duplex to strand transition.(ABSTRACT TRUNCATED AT 400 WORDS)     

Systematic Analysis of Possible Hydrogen Bonds between Amino Acid Side Chains and B-form DNA

Abstract

The ways in which amino acid side chains could make a pair of hydrogen bonds within the major groove of B DNA are systematically analyzed. Hydrogen bond donors within the major groove are characterized by determining the idealized position of the hydrogen bond acceptors that they might bond with. It appears that an amino acid side chain could, at most, contact a pair of base pairs. The ten possible pairs of base pairs are analyzed to determine how they could be recognized by the amino acid side chains.     

Structural diversity and isomorphism of hydrogen-bonded base interactions in nucleic acids

The wide structural diversity of RNA results in part from the diversity of non-Watson-Crick interactions between bases. To examine the repertoire of possible hydrogen bond interactions among bases, we computed databases of base-pairs and base-triples by systematically matching all possible hydrogen-bond donors and acceptors between bases and evaluating the geometries of each planar configuration. For base-pairs, we find 53 arrangements having at least two hydrogen bonds, including 23 pairs with protonated bases that have not previously been modeled. A comparison with experimentally observed base-pairs reveals an unexpected G:U pair recently observed in the ribosome. For base-triples, we find 840 arrangements in which the three bases are constrained by a total of at least three hydrogen bonds. Base-triples in particular exhibit a wide range of structural diversity, suggesting how compact or elongated nucleic acid structures may be constructed using different hydrogen-bonding patterns. Base-pair and base-triple conformations were systematically compared to identify structurally isomorphic combinations, and the experimentally observed arrangements within double and triple helices are among the most isomorphic. Unexpectedly, however, other combinations in the database are even more isomorphic, including several in which all-purine arrangements overlap with all-pyrimidine arrangements. These studies highlight some of the combinatoric and geometric versatility of base interactions and help provide a framework for analyzing and modeling isomorphic interactions and potentially for designing novel nucleic acid structures.     

Simulation of interactions between nucleic acid bases by refined atom-atom potential functions

Energy of interaction between nitrogen bases of nucleic acid has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H . . . O and (or) N-H . . . N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H . . . O or N-H . . . N and one C-H . . . O or C-H . . . N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum-mechanical approach.     

Structural Features and Hydration of d(C-G-C-G-A-A-T-T-A-G-C-G); a Double Helix Containing Two G.A Mispairs

Abstract

Single crystal X-ray diffraction techniques have been used to characterise the molecular structure of the title compound to 2.5Å resolution. The structure consists of ten standard Watson-Crick base pairs and two G.A mismatched base pairs. The purine-purine mismatches have guanine in the usual anti orientation with respect to the sugar and adenine in syn orientation. There are two hydrogen bonds formed between the mismatch bases, N-l and 0–6 of guanine with N-7 and N-6 of adenine respectively. The bulky purine-purine mismatches are accommodated with minor perturbation of the sugar-phosphate backbone. There is a slight improvement in base pair overlap at the mismatch sites. Details of the backbone conformation, base stacking interactions and hydration are presented and compared with those of the parent compound d(C-G-C-G-A-A-T-T-C-G-C-G).     

Simulation of Interactions between Nucleic Acid Bases by Refined Atom-Atom Potential Functions

Abstract

Energy of interaction between nitrogen bases of nucleic acids has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H … O and (or) N-H … N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component, in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H … O or N-H … N and one C-H … O or C-H … N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum- mechanical approach.     

Silver(I)-mediated Hoogsteen-type base pairs

Metal-mediated Hoogsteen-type base pairs are useful for the construction of DNA duplexes containing contiguous stretches of metal ions along the helical axis. To fine-tune the stability of such base pairs and the selectivity toward different metal ions, the availability of a selection of artificial nucleobases is highly desirable. In this study, we follow a theoretical approach utilizing dispersion-corrected density functional methods to evaluate a variety of artificial nucleobases as candidates for metal-mediated Hoogsteen-type base pairs. We focus on silver(I)-mediated Hoogsteen- and reverse Hoogsteen-type base pairs formed between 1-deaza- and 1,3-dideazapurine-derived nucleobases, respectively, and cytosine. Apart from two coordinative bonds, these base pairs are stabilized by a hydrogen bond. We elucidate the impact of different substituents at the C6 position and the presence or absence of an endocyclic N3 nitrogen atom on the overall stability of a base pair and concomitantly on the strength of the hydrogen and coordinative bonds. All artificial base pairs investigated in this study are less stable than the experimentally established benchmark base pair C–Ag⁺–G. The base pair formed from 1,3-dideaza-6-methoxypurine is isoenergetic to the experimentally observed C–Ag⁺–C base pair. This makes 1,3-dideaza-6-methoxypurine a promising candidate for the use as an artificial nucleobase in DNA.     

NMR structures of r(GCAGGCGUGC)2 and determinants of stability for single guanosine-guanosine base pairs

Nucleotides in RNA that are not Watson-Crick-paired form unique structures for recognition or catalysis, but determinants of these structures and their stabilities are poorly understood. A single noncanonical pair of two guanosines (G) is more stable than other noncanonical pairs and can potentially form pairing structures with two hydrogen bonds in four different ways. Here, the energetics and structure of single GG pairs are investigated in several sequence contexts by optical melting and NMR. The data for r(5'GCAGGCGUGC3')(2), in which G4 and G7 are paired, are consistent with a model in which G4 and G7 alternate syn glycosidic conformations in a two-hydrogen-bond pair. The two distinct structures are derived from nuclear Overhauser effect spectroscopic distance restraints coupled with simulated annealing using the AMBER 95 force field. In each structure, the imino and amino protons of the anti G are hydrogen bonded to the O6 and N7 acceptors of the syn G, respectively. An additional hydrogen-bond connects the syn G amino group to the 5' nonbridging pro-R(p) phosphate oxygen. The GG pair fits well into a Watson-Crick helix. In r(5'GCAGGCGUGC3')(2), the G4(anti), G7(syn) structure is preferred over G4(syn), G7(anti). For single GG pairs in other contexts, exchange processes make interpretation of spectra more difficult but the pairs are also G(syn), G(anti). Thermodynamic data for a variety of duplexes containing pairs of G, inosine, and 7-deazaguanosine flanked by GC pairs are consistent with the structural and energetic interpretations for r(5'GCAGGCGUGC3')(2), suggesting similar GG conformations.     

NMR structures of (rGCUGAGGCU)2 and (rGCGGAUGCU)2: probing the structural features that shape the thermodynamic stability of GA pairs

The NMR structures of [see text] and [see text] are reported. The internal loop, [see text], is about 2 kcal/mol more stable than [see text] at 37 degrees C. The duplexes assemble into similar global folds characterized by the formation of tandem sheared GA pairs. The different stabilities of the loops are accompanied by differences in the local structure of the closing GU pairs. In the [see text] internal loop, the GU pairs form canonical wobble configurations with two hydrogen bonds, whereas in [see text], the GU pairs form a single hydrogen bond involving the amino group, GH22, and the carbonyl group, UO4. This pairing is similar to the GU closing pair of the 690 hairpin loop found in E. coli 16S rRNA. The [see text] and [see text] structures reveal how the subtle interplay between stacking and hydrogen bonding determines sequence dependent conformation and thermodynamic stability. Thus, this work provides structural and thermodynamic benchmarks for theoreticians in the ongoing effort to understand the sequence dependence of RNA physicochemical properties.     

Formation of purine-purine mispairs by Sulfolobus solfataricus DNA polymerase IV

DNA damage that stalls replicative polymerases can be bypassed with the Y-family polymerases. These polymerases have more open active sites that can accommodate modified nucleotides. The lack of protein-DNA interactions that select for Watson-Crick base pairs correlate with the lowered fidelity of replication. Interstrand hydrogen bonds appear to play a larger role in dNTP selectivity. The mechanism by which purine-purine mispairs are formed and extended was examined with Solfolobus solfataricus DNA polymerase IV, a member of the RAD30A subfamily of the Y-family polymerases, as is pol eta. The structures of the purine-purine mispairs were examined by comparing the kinetics of mispair formation with adenine versus 1-deaza- and 7-deazaadenine and guanine versus 7-deazaguanine at four positions in the DNA, the incoming dNTP, the template base, and both positions of the terminal base pair. The time course of insertion of a single dNTP was examined with a polymerase concentration of 50 nM and a DNA concentration of 25 nM with various concentrations of dNTP. The time courses were fitted to a first-order equation, and the first-order rate constants were plotted against the dNTP concentration to produce k pol and K d (dNTP) values. A decrease in k pol/ K d (dNTP) associated with the deazapurine substitution would indicate that the position is involved in a crucial hydrogen bond. During correct base pair formation, the adenine to 1-deazaadenine substitution in both the incoming dNTP and template base resulted in a >1000-fold decrease in k pol/ K d (dNTP), indicating that interstrand hydrogen bonds are important in correcting base pair formation. During formation of purine-purine mispairs, the k pol/ K d (dNTP) values for the insertion of dATP and dGTP opposite 7-deazaadenine and 7-deazaguanine were decreased >10-fold with respect to those of the unmodified nucleotides. In addition, the rate of incorporation of 1-deaza-dATP opposite guanine was decreased 5-fold. These results suggest that during mispair formation the newly forming base pair is in a Hoogsteen geometry with the incoming dNTP in the anti conformation and the template base in the syn conformation. These results indicate that Dpo4 holds the incoming dNTP in the normal anti conformation while allowing the template nucleotide to change conformations to allow reaction to occur. This result may be functionally relevant in the replication of damaged DNA in that the polymerase may allow the template to adopt multiple configurations.