Strong doublet preferences in nucleotide sequences and DNA geometry

Summary Analysis of the sequence data available today, comprising more than 500,000 bases, confirms the previously observed phenomenon that there are distinct dinucleotide preferences in DNA sequences. Consistent behaviour is observed in the major sequence groups analysed here in prokaryotes, eukaryotes and mitochondria. Some doublet preferences are common to all groups and are found in most sequences of the Los Alamos Library. The patterns seen in such large data sets are very significant statistically and biologically. Since they are present in numerous and diverse nucleotide sequences, one may conclude that they confer evolutionary advantages on the organism.In eukaryotes RR and YY dinucleotides are preferred over YR and RY (where R is a purine and Y a pyrimidine). Since opposite-chain nearest-neighbour purine clashes are major determinants of DNA structure, it appears that the tight packaging of DNA in nucleosomes disfavors, in general, such (YR and RY) steric repulsion.     

Strong patterns in homooligomer tracts occurrences in non-coding and in potential regulatory sites in eukaryotic genomes

Previous studies of the dinucleotides flanking both the 5' and 3' ends of homooligomer tracts have shown that some flanks are consistently preferred over others (1,2). In the first preferred group, the homooligomer tracts are flanked by the same nucleotide and/or the complementary nucleotides, e.g.,ATAn,TTAn,CCGn, where n = 2-5. Runs flanked by nucleotides with which they cannot base pair are distinctly disfavored. (In this group An/Tn are flanked by C and/or G; Gn/Cn are flanked by A/T, e.g.,CGAn,TnGG,GnAT). The frequencies of runs flanked by A or T, and G or C ("mixed"group) are as expected. Here we seek the origin of this effect and its relevance to protein-DNA interactions. Surprisingly, within the first group, runs flanked by their complements with a pyrimidine-purine junction (e.g.,TTAn,CnGG) are greatly preferred. The frequencies of their purine-pyrimidine junction mirror-images is just as expected. This effect, as well as additional ones enumerated below, is seen universally in eukaryotes and in prokaryotes, although it is stronger in the former. Detailed analysis of regulatory regions shows these strong trends, particularly in GC sequences. The potential relationship to DNA conformation and DNA-protein interaction is discussed.     

Strong Patterns in Homooligomer Tracts Occurrences in Non-Coding and in Potential Regulatory Sites in Eukaryotic Genomes

Abstract

Previous studies of the dinucleotides flanking both the 5′ and 3′ ends of homooligomer tracts have shown that some flanks are consistently preferred over others (1,2). In the first preferred group, the homooligomer tracts are flanked by the same nucleotide and/or the complementary nucleotides, e.g., ATA_n, TTA_n, CCG_n, where n=2–5. Runs flanked by nucleotides with which they cannot base pair are distinctly disfavored. (In this group A/T_n are flanked by C and/or G; G_n/C_n are flanked by A/T, e.g., CGA_n, T_nGG, G., AT). The frequencies of runs flanked by AorT, and G or C (“mixed” group) are as expected. Here we seek the origin of this effect and its relevance to protein-DNA interactions. Surprisingly, within the first group, runs flanked by their complements with a pyrimidine-purine junction (e.g., TTA_n, C_nGG) are greatly preferred. The frequencies of their purine-pyrimidine junction mirror-images is just as expected. This effect, as well as additional ones enumerated below, is seen universally in eukaryotes and in prokaryotes, although it is stronger in the former. Detailed analysis of regulatory regions shows these strong trends, particularly in GC sequences. The potential relationship to DNA conformation and DNA-protein interaction is discussed.     

Sequence context of oligomer tracts in eukaryotic DNA: biological and conformational implications

Recent studies of homooligomer tracts suggest different characteristics from random sequence DNA (dA).(dT) and (dG).(dC) tracts are frequent in upstream regions and in some cases have been shown to be essential for regulation. Here we examine homooligomer occurrences in non-coding and coding eukaryotic sequences, focusing on the context in which the homooligomers occur. This analysis of sequences in the junction areas yields distinct and consistent characteristics. In particular, the nucleotide interrupting a run is most frequently complementary to the run. The base next to it is most frequently identical to the one constituting the run. For A or T runs the least frequent nearest and next to nearest neighbors are G or C. For G or C tracts the least frequent are A or T. Complementary oligomers behave similarly. These and additional trends are strongest for run lengths greater than or equal to 3. The computations are carried out on the whole eukaryotic database of greater than 4 x 10(6) nucleotides, separately for coding and non-coding regions. These same trends are evident for both groups, but are somewhat stronger for the non-coding regions. The context in which the homooligomers occur may yield some clues to DNA conformation and its biological implications.     

DNA basepair step deformability inferred from molecular dynamics simulations

The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.     

Trends in the 5' vs. 3' flanks of oligonucleotides in eukaryotic and prokaryotic genomes: the asymmetric roles played by cytosine and guanine.

Studies of sequence context preferences of oligonucleotides composed of (G/C)n and (A/T)m blocks (n + m = 3,4,5) unravel strong patterns. Comparisons of the 5' and 3' nearest neighbor doublets flanking these oligomers reveal the preference of (G/C)2 to be positioned immediately next to the (A/T)m block, enclosing it by (G/C) nucleotides rather than extending the (G/C)n block. That is, for a (G/C)n(A/T)m oligomer and a (G/C)2 doublet, (G/C)n(A/T)m(G/C)2 greater than (G/C)n + 2 (A/T)m. Similarly for an (A/T)m(G/C)n oligomer, (G/C)2(A/T)m(G/C)n greater than (A/T)m(G/C)n + 2. In an analogous manner, (A/T)2 flanking doublets prefer enclosing the (G/C)n blocks, although these patterns are weaker. Here we show a strong, direct relationship between the magnitude of the trends and the presence of Cs in the (G/C)n block in the (G/C)n(A/T)m oligomer, and the presence of Gs in the complementary (A/T)m(G/C)n oligomers. The trends are stronger in eukaryotic than in prokaryotic sequences. They are stronger for longer (G/C)n and shorter (A/T)m blocks. We suggest that the preference for (A/T)m to be enclosed by (G/C) rather than be flanked by them on only one side is related to DNA structure and DNA-protein interaction. Sequences of the (G/C)(A/T)(G/C) type may have more homogeneous minor groove geometry. In particular, the strong G vs. C asymmetry in the trends may be related to pyrimidine-purine junctions, possibly to CG sequences.     

Zipper-like Watson-Crick base-pairs

A series of DNA heptadecamers containing the DNA analogues of RNA E-like 5'-d(GXA)/(AYG)-5' motifs (X/Y is complementary T/A, A/T, C/G, or G/C pair) were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers reveal excellent resolution in NMR spectra and exhibit many unusual nuclear Overhauser effects (NOEs) that allow for good characterization of an unusual zipper-like conformation with zipper-like Watson-Crick base-pairs; the potential canonical X.Y H-bonding is not present, and the central X/Y pairs are transformed instead into inter-strand stacks that are bracketed by sheared G.A base-pairs. Such phenomenal structural change is brought about mainly through two backbone torsional angle adjustments, i.e. delta from C2'-endo to C3'-endo for the sugar puckers of unpaired residues and gamma from gauche(+) to trans for the following 3'-adenosine residues. Such motifs are analogous to the previously studied (GGA)(2) motif presumably present in the human centromeric (TGGAA)(n) tandem repeat sequence. The novel zipper-like motifs are only 4-7 deg. C less stable than the (GGA)(2) motif, suggesting that inter-strand base stacking plays an important role in stabilizing unusual nucleic acid structures. The discovery that canonical Watson-Crick G.C or A.T hydrogen-bonded pairs can be transformed into stacking pairs greatly increases the repertoire for unusual nucleic acid structural motifs.     

Distinct patterns in the dinucleotide nearest neighbors to G/C and A/T oligomers in eukaryotic sequences

Summary The eukaryotic and prokaryotic databases are scanned for potential nearest-neighbor doublet preferences at the 5 and 3 flanks of some oligomers. Here we focus on oligomers containing alternating nucleotides, i.e., UV, UVUV, and UUVV where UV. Strong, consistent trends are observed in eukaryotic sequences. A/T alternation oligomers are preferentially flanked by A/T. G/C flanks are disfavored. G/C alternation oligomers are preferentially flanked by G/C. A/T flanks are disfavored. These trends are consistent with those observed previously for homooligomer tracts (Nussinov et al. 1989a,b). G/C tracts are preferentially flanked by G/C. A/T nearest neighbors are disfavored. The reverse holds for A/T tracts. Additional patterns are described here as well. The possible origin of these DNA composition and sequence trends is discussed. These trends are suggested to stem from protein-DNA interaction constraints.     

Eukaryotic topoisomerase II preferentially cleaves alternating purine-pyrimidine repeats.

Alternating purine-pyrimidine sequences (RY repeats) demonstrate considerable homology to the consensus sequence for vertebrate topoisomerase II (Spitzner and Muller (1988) Nucleic Acids Res. 16: 1533-1556). This is shown below and positions that can match are underscored. RYRYRYRYRYRYRYRYRY = alternating purine-pyrimidine 18 bp RNYNNCNNGYNGKTNYNY = topoisomerase II consensus sequence (R is purine, Y is pyrimidine, K is G or T.) Topoisomerase II cleavage reactions were performed (in the absence of inhibitors) on a plasmid containing a 54 base RY repeat and the single strong cleavage site mapped to the RY repeat. Analysis of this DNA on sequencing gels showed that the enzyme cleaved a number of sites, all within the 54 base pair RY repeat. Topoisomerase II also made clustered cleavages within other RY repeats that were examined. Quantitative analysis of homology to the consensus sequence, as measured by the match of a site to a matrix of base proportions from the consensus data base (the matrix mean), showed that both the locations and the frequencies of cleavage sites within RY repeats were proportional to homology scores. However, topoisomerase II cleaved RY repeats preferentially in comparison to non-RY sites with similar homology scores. The activity of the enzyme at RY repeats appears to be proportional to the length of the repeat; additionally, GT, AC and AT repeats were better substrates for cleavage than GC repeats.     

The 6-thioguanine/5-methyl-2-pyrimidinone base pair.

As part of a program to determine the physical possibility of expanding the number of types of base pairs in DNA, the pairing stabilities of the analog bases 6-thioguanine (GS) and 5-methyl-2-pyrimidinone (TH) in oligodeoxynucleotides were measured. Procedures were developed to synthesize oligodeoxynucleotides with the analog bases. The sequences of the synthesized oligomers were T-C-G-A-C-G-G-X-Y-C-C-G. An enzymatic procedure was developed to measure relative association constants of oligomer pairs with the self complementary reference oligomer, X = A and Y = T, K(T/A) = K. The results were K(C/G) = (5 +/- .5)K, K(TH/GS) = K/(1 +/- .5), K(T/G) = K/(9 +/- 3), K(TH/G) = K/(25 +/- 5), K(C/GS) less than K/30, K(TH/A) less than K/40, K(T/GS) less than K/40, K(C/A) less than K/40. The results with the standard bases are consistent with other methods of measurement. The stability of the base pair GS/TH is approximately the same as the standard base pair A/T.     

NMR structure of a parallel-stranded DNA duplex at atomic resolution

DNA dodecamers have been designed with two cytosines on each end and intervening A and T stretches, such that the oligomers have fully complementary A:T base pairs when aligned in the parallel orientation. Spectroscopic (UV, CD and IR), NMR and molecular dynamics studies have shown that oligomers having the sequences d(CCATAATTTACC) and d(CCTATTAAATCC) form a parallel-stranded duplex when dissolved at 1:1 stoichiometry in aqueous solution. This is due to the C:C⁺ clamps on either end and extensive mismatches in the antiparallel orientation. The structure is stable at neutral and acidic pH. At higher temperatures, the duplex melts into single strands in a highly cooperative fashion. All adenine, cytosine and thymine nucleotides adopt the anti conformation with respect to the glycosidic bond. The A:T base pairs form reverse Watson–Crick base pairs. The duplex shows base stacking and NOEs between the base protons T(H6)/A(H8) and the sugar protons (H1′/H2′/H2″) of the preceding nucleotide, as has been observed in antiparallel duplexes. However, no NOEs are observed between base protons H2/H6/H8 of sequential nucleotides, though such NOEs are observed between T(CH₃) and A(H8). A three-dimensional structure of the parallel-stranded duplex at atomic resolution has been obtained using molecular dynamics simulations under NMR constraints. The simulated structures have torsional angles very similar to those found in B-DNA duplexes, but the base stacking and helicoid parameters are significantly different.     

A duplex of the oligonucleotides d(GGGGGTTTTT) and d(AAAAACCCCC) forms an A to B conformational junction in concentrated salt solutions

The coexistence of both A form and B form tracts and formation of an A-B junction in the oligomer d(GGGGGTTTTT).d(AAAAACCCCC) in saturated sodium chloride solution have been detected by Raman spectroscopy. The entire duplex adopts the familiar B-form conformation in aqueous solution at low salt concentrations (0.1M NaCl). In 6M NaCl the adoption of an A form is observed within the G,C tract while a B-form is maintained in the A.T tract. The experimental results indicate that two different helical forms can co-exist in a rather short oligonucleotide and that formation of an A-B junction can occur over a fairly small span of bases. This is in agreement with recent rules governing the relation between base sequence and secondary structure of DNA published from this laboratory. The conformational preferences of each of the individual oligomers d(AAAAACCCCC) and d(GGGGGTTTTT) have also been investigated. The oligomer d(AAAAACCCCC) is single stranded but some evidence for base stacking is observed at 2 degrees C. In contrast, a double stranded B-form structure characterized by wobble G-T base pairing is observed for d(GGGGGTTTTT) in 0.1M and 6M NaCl.     

Sequence-dependent bending of DNA and phasing of nucleosomes

Conformational analysis has revealed anisotropic flexibility of the B-DNA double helix: it bends most easily into the grooves, being the most rigid when bent in a perpendicular direction. This result implies that DNA in a nucleosome is curved by means of relatively sharp bends ("mini-kinks") which are directed into the major and minor grooves alternatively and separated by 5-6 base pairs. The "mini-kink" model proved to be in keeping with the x-ray structure of the B-DNA dodecamer resolved later, which exhibits two "annealed kinks", also directed into the grooves. The anisotropy of B DNA is sequence-dependent: the pyrimidine-purine dimers (YR) favor bending into the minor groove, and the purine-pyrimidine dinucleotides (RY), into the minor one. The RR and YY dimers appear to be the most rigid dinucleotides. Thus, a DNA fragment consisting of the interchanging oligopurine and oligopyrimidine blocks 5-6 base pairs long should manifest a spectacular curvature in solution. Similarly, a nucleotide sequence containing the RY and YR dimers separated by a half-pitch of the double helix is the most suitable for wrapping around the nucleosomal core. Analysis of the numerous examples demonstrating the specific alignment of nucleosomes on DNA confirms this concept. So, the sequence-dependent "mechanical" properties of the double helix influence the spatial arrangement of DNA in chromatin.     

Influence of sequence context and length on the structure and stability of triplet repeat DNA oligomers

Genetic expansion diseases have been linked to the properties of triplet repeat DNA sequences during replication. The most common triplet repeats associated with such diseases are CAG, CCG, CGG, and CTG. It has been suggested that gene expansion occurs as a result of hairpin formation of long stretches of these sequences on the leading daughter strand synthesized during DNA replication [Gellibolian, R., Bacolla, A., and Wells, R. D. (1997) J. Biol. Chem. 272, 16793-7]. To test the biophysical basis for this model, oligonucleotides of general sequence (CNG)(n), where N = A, C, G, or T and n = 4, 5, 10, 15, or 25, were synthesized and characterized by circular dichroism (CD) spectropolarimetry, optical melting studies, and differential scanning calorimetry (DSC). The goal of these studies was to evaluate the influence of sequence context and oligomer length on their secondary structures and stabilities. The results indicate that all single oligomers, even those as short as 12 nucleotides, form stable hairpin structures at 25 degrees C. Such hairpins are characterized by the presence of N:N mismatched base pairs sandwiched between G:C base pairs in the stems and loops of three to four unpaired bases. Thermodynamic analysis of these structures reveals that their stabilities are influenced by both the sequence of the particular oligomer and its length. Specifically, the stability order of CGG > CTG > CAG > CCG was observed. In addition, longer oligomers were found to be more stable than shorter oligomers of the same sequence. However, a stability plateau above 45 nucleotides suggests that the length dependence reaches a maximum value where the stability of the G:C base pairs can no longer compensate the instability of the N:N mismatches in the stems of the hairpins. The results are discussed in terms of the above model proposed for gene expansion.     

A Duplex of the Oligonucleotides d(GGGGGTTTTT) and d(AAAAACCCCC) Forms an A to B Conformational Junction in Concentrated Salt Solutions

Abstract

The coexistence of both A form and B form tracts and formation of an A-B junction in the oligomer d(GGGGGTTTTT)· d(AAAAACCCCC) in saturated sodium chloride solution have been detected by Raman spectroscopy. The entire duplex adopts the familiar B-form conformation in aqueous solution at low salt concentrations (0.1M NaCl). In 6M NaCl the adoption of an A form is observed within the G,C tract while a B-form is maintained in the A,T tract. The experimental results indicate that two different helical forms can co-exist in a rather short oligonucleotide and that formation of an A-B junction can occur over a fairly small span of bases. This is in agreement with recent rules governing the relation between base sequence and secondary structure of DNA published from this laboratory.

The conformational preferences of each of the individual oligomers d(AAAAACCCCC) and d(GGGGGTTTTT) have also been investigated. The oligomer d(AAAAACCCCC) is single stranded but some evidence for base stacking is observed at 2°C. In contrast, a double stranded B-form structure characterized by wobble G-T base pairing is observed for d(GGGGGTTTTT) in 0.1M and 6M NaCl.     

Characterization of the biological activities of uridine diphosphate in human dendritic cells: Influence on chemotaxis and CXCL8 release

Uridine nucleotides are endogenous nucleotides which are released into the extracellular space from mechanical stressed endothelial and epithelial cells as well as lipopolysaccharide (lps)-stimulated monocytes. Here, we studied the biological activity of the selective purinoreceptor P2Y6 (P2YR6) agonist Uridine 5'diphosphate (UDP) as well as the P2YR2- and P2YR4-activating uridine 5'triphosphate (UTP) on human dendritic cells (DC). These cells in their immature state have the ability to migrate from blood to peripheral target sites where they sense dangerous signals and capture potential antigens. Moreover, mature DC induce innate immune responses and migrate from peripheral tissues to secondary lymphoid organs in order to activate naive T cells and initiate adaptive immunity. Here, we were able to show that uridine nucleotides stimulated Ca(2+) transients, actin polymerization, and chemotaxis in immature DC. Experiments with pertussis toxin, the stable pyrimidine agonist uridine 5'-O-(2-thiodiphosphate) (UDPgammaS) and receptor antagonists, as well as desensitization studies suggested that these uridine nucleotides activities were mediated by different G(i) protein-coupled receptors. During lps-induced maturation, DC lost their ability to respond towards uridine nucleotides with these activities. Instead, UDP, but not UTP, stimulated the release of the CXC-chemokine 8 (CXCL8) from mature DC in a reactive blue sensitive manner. Moreover, our study indicates that UDP stimulates different signaling pathways in immature and mature DC in order to favor the accumulation of immature DC and to augment the capacity to secrete CXCL8 in mature DC.     

Solution structure of an LNA hybridized to DNA: NMR study of the d(CT(L)GCT(L)T(L)CT(L)GC):d(GCAGAAGCAG) duplex containing four locked nucleotides

We have used two-dimensional (1)H NMR spectroscopy at 750 MHz to determine a high-resolution solution structure of an oligonucleotide containing restricted nucleotides with a 2'-O, 4'-C-methylene bridge (LNA) hybridized to the complementary DNA strand. The LNA:DNA duplex examined contained four thymidine LNA modifications (T(L), d(C1T(L)2G3C4T(L)5T(L)6C7T(L)8G9C10):d( G11C12A13G14A15A16G17C 18A19G20). A total relaxation matrix approach was used to obtain interproton distance bounds from NOESY cross-peak intensities. These distance bounds were used as restraints in molecular dynamics (rMD) calculations. Forty final structures were generated for the duplex from A-form and B-form DNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the 40 structures of the complex was 0.6 A. The sugar puckerings are averaged values of a dynamic interchange between N- and S-type conformation except in case of the locked nucleotides that were found to be fixed in the C3'-endo conformation. Among the other nucleotides in the modified strand, the furanose ring of C7 and G9 is predominantly in the N-type conformation whereas that of G3 is in a mixed conformation. The furanose rings of the nucleotides in the unmodified complementary strand are almost exclusively in the S-type conformation. Due to these different conformations of the sugars in the two strands, there is a structural strain between the A-type modified strand and the B-type unmodified complementary strand. This strain is relaxed by decreasing the value of rise and compensating with tip, buckle, and propeller twist. The values of twist vary along the strand but for a majority of the base pairs a value even lower than that of A-DNA is observed. The average twist over the sequence is 32+/-1 degrees. On the basis of the structure, we conclude that the high stability of LNA:DNA duplexes is caused by a local change of the phosphate backbone geometry that favors a higher degree of stacking.     

Human DNA polymerase iota utilizes different nucleotide incorporation mechanisms dependent upon the template base

Human DNA polymerase ι (Polι) is a member of the Y family of DNA polymerases involved in translesion DNA synthesis. Polι is highly unusual in that it possesses a high fidelity on template A, but has an unprecedented low fidelity on template T, preferring to misincorporate a G instead of an A. To understand the mechanisms of nucleotide incorporation opposite different template bases by Polι, we have carried out pre-steady-state kinetic analyses of nucleotide incorporation opposite templates A and T. These analyses have revealed that opposite template A, the correct nucleotide is preferred because it is bound tighter and is incorporated faster than the incorrect nucleotides. Opposite template T, however, the correct and incorrect nucleotides are incorporated at very similar rates, and interestingly, the greater efficiency of G misincorporation relative to A incorporation opposite T arises predominantly from the tighter binding of G. Based on these results, we propose that the incipient base pair is accommodated differently in the active site of Polι dependent upon the template base and that when T is the templating base, Polι accommodates the wobble base pair better than the Watson-Crick base pair.     

Hypermutation at A-T base pairs: the A nucleotide replacement spectrum is affected by adjacent nucleotides and there is no reverse complementarity of sequences flanking mutated A and T nucleotides

Hypermutation is thought to be a two-phase process. The first phase is via the action of activation-induced cytidine deaminase (AID), which deaminates C nucleotides in WRC motifs. This results in the RGYW/WRCY hot spot motifs for mutation from G and C observed in vivo. The resemblance between the hot spot for C mutations and the reverse complement of that for G mutations implies a process acting equally on both strands of DNA. The second phase of hypermutation generates mutations from A and T and exhibits strand bias, with more mutations from A than T. Although this does not concur with the idea of one mechanism acting equally on both strands, it has been suggested that the AT mutator also has a reversible motif; WA/TW. We show here that the motifs surrounding the different substitutions from A vary significantly; there is no single targeting motif for all A mutations. Sequence preferences associated with mutations from A more likely reflect an influence of adjacent nucleotides over what the A mutates "to." This influence tends toward "like" replacements: Purines (A or G) in the 5' position bias toward replacement by another purine (G), whereas replacement with pyrimidines (C or T) is more likely if the preceding base is also a pyrimidine. There is no reverse complementarity in these observations, in that similar influences of nucleotides adjacent to T are not seen. Hence, WA and TW should not be considered as reverse complement hot spot motifs for A and T mutations.