The Intramolecular Impact to the Sequence Specificity of B→A Transition: Low Energy Conformational Variations in AA/TT and GG/CC Steps

Abstract It is well known, that local B→A transformation in DNA is involved in several biological processes. In vitro B?A transition is sequence-specific. The physical basis of this specificity is not known yet. Here we analyze the effect of intramolecular interactions on the structural behavior of the GG/CC and AA/TT steps. These steps exemplify sequence specific bias to the B- or A-form structure. Optimization of potential energy of the molecular systems composed of an octanucle-otide, neutralized by Na(+) and solvated with TIP3P water molecules in rectangular box with periodic boundary conditions gives the statistically representative sets of low energy structures for GG/CC and AA/TT steps in the middle of the diverse flanking sequences. Permissible 3D variations of GG/CC and AA/TT, and correlation of the relative motion of base pairs in these steps were analyzed. AA/TT step permits high variability for low energy conformers in the B-form DNA and small variability for low energy conformers in the A-form DNA. In contrast GG/CC step permits high variability for low energy conformers in the A-form DNA and small variability for low energy conformers in the B-form DNA. The relative motion of base pairs in GG/CC step is high correlated, while in AA/TT step this correlation is notably less. Atom-atom interactions inside-the-step always favors the B-form and their component - stacking interactions (atomatom interactions between nucleic bases) is crucial for the duplex stabilization. Formation of the A-form for both steps is a result of interactions with the flanking sequences and water-cation environment in the box. The average energy difference between conformations presenting B-form and A-form for the GG/CC step is high, while for the AA/TT step it is rather low. Thus, intramolecular interactions in GG/CC and AA/TT steps affect the possible conformational diversity ("conformational entropy") of the A- and B- type structures of DNA step. This determines the known bias of the A-form DNA depending on the enrichment of sequences with GG/CC. If structural tuning during the process of protein-DNA complex formation lead to the local B→A transformation of DNA, it is largely directed by high conformational diversity of GG/CC step in the A-form. In such a case the presence in the target site of both kinds of examined steps ensures the reversible character of ligand binding.     

Selective binding of the TATA box-binding protein to the TATA box-containing promoter: analysis of structural and energetic factors.

We report the results of an energy-based exploration of the components of selective recognition of the TATA box-binding protein (TBP) to a TATA box sequence that includes 1) the interaction between the hydrophobic Leu, Pro, and Phe residues of TBP with the TA, AT, AA, TT, and CG steps, by ab initio quantum mechanical calculations; and 2) the free energy penalty, calculated from molecular dynamics/potential of mean force simulations, for the conformational transition from A-DNA and B-DNA into the TA-DNA form of DNA observed in a complex with TBP. The GTAT, GATT, GAAT, and GTTT tetramers were explored. The results show that 1) the discrimination of TA, AT, AA, TT, or CG steps by TBP cannot rest on their interaction with the inserting Phe side chains; 2) the steric clash between the bulky and hydrophobic Pro and Leu residues and the protruding -NH2 group of guanine is responsible for the observed selectivity against any Gua-containing basepair; 3) the Pro and Leu residues cannot selectively discriminate among TA, AT, AA, or TT steps; and 4) the calculated energy required to achieve the TA-DNA conformation of DNA that is observed in the complex with TBP appears to be a key determinant for the observed selectivity against the AT, AA, and TT steps. The simulations also indicate that only the TA step can form a very efficient interbase hydrogen bond network in the TA-DNA conformation. Such an energetically stabilizing network is not achievable in the AA and TT steps. While it is viable in the AT step, structural constraints render the hydrogen bonding network energetically ineffective there.     

C-H.O hydrogen bonds in minor groove of A-tracts in DNA double helices

Analysis of available B-DNA type oligomeric crystal structures as well as protein-bound DNA fragments (solved using data with resolution <2.6 A) indicates that in both data sets, a majority of the (3'-Ade) H2..O2(3'-Thy/Cyt) distances in AA.TT and GA.TC dinucleotide steps, are considerably shorter than their values in a uniform fibre model, and are smaller than their optimum separation distance. Since the electropositive C2-H2 group of adenine is in close proximity of the electronegative keto oxygen atoms of both pyrimidine bases in the antiparallel strand of the double-helical DNA structures, it suggests the possibility of intra-base-pair as well as cross-strand C-H..O hydrogen bonds in the minor groove. The C2-H2..O2 hydrogen bonds within the A.T base-pairs could be a natural consequence of Watson-Crick pairing. However, the close cross-strand interactions between the bases at the 3'-ends of the AA.TT and GA.TC steps arise due to the local sequence-dependent geometry of these steps. While the base-pair propeller twist in these steps is comparable to the fibre model, some of the other local parameters such as base-pair opening angle and inter-base-pair slide show coordinated changes, leading to these shorter C2-H2..O2 distances. Hence, in addition to the well-known minor groove hydration, it appears that favourable C2-H2..O2 cross-strand interactions may play a role in imparting a characteristic geometry to AA.TT and GA.TC steps, as well as An.Tn and GAn.TnC tracts, which leads to a narrow minor groove in these regions.     

Characterization of inherent curvature in DNA lacking polyadenine runs.

Sequence-directed DNA curvature is most commonly associated with AA dinucleotides in the form of polyadenine runs. We demonstrate inherent curvature in DNA which lacks AA/TT dinucleotides using the criteria of polyacrylamide gel mobility and efficiency of DNA cyclization. These studies are based upon two 21-base pair synthetic DNA fragments designed to exhibit fixed curvature according to deflections made to the helical axis by non-AA dinucleotide stacks. Repeats of these sequences display anomalously slow migration in polyacrylamide gels. Moreover, both sequences describe helical conformations that are closed into circles by DNA ligase at much smaller sizes than is typical of nondeformed DNA. Chemical cleavage of these DNA molecules with hydroxyl radical is also consistent with local variation in helical conformation at specific dinucleotide steps.     

Sequence dependence of DNA conformational flexibility

By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.     

Three sequence rules for chromatin

Extensive DNA sequence analysis of three eukaryotes, S. cerevisiae, C. elegans, and D. melanogaster, reveals two different AA/TT periodical patterns associated with the nucleosome positioning. The first pattern is the counter-phase oscillation of AA and TT dinucleotides, which has been frequently considered as the nucleosome DNA pattern. This represents the sequence rule I for chromatin structure. The second pattern is the in-phase oscillation of the AA and TT dinucleotides with the same nucleosome DNA period, 10.4 bases. This pattern apparently corresponds to curved DNA, that also participates in the nucleosome formation, and represents the sequence rule II for chromatin. The positional correlations of AA and TT dinucleotides also indicate that the nucleosomes are separated by specific linker sizes (preferably 8, 18, ... bases), dictated by the steric exclusion rules. Thus, the sequence positions of the neighboring nucleosomes are correlated, and this represents the sequence rule III.     

Three Sequence Rules for Chromatin

Abstract

Extensive DNA sequence analysis of three eukaryotes, S. cerevisiae, C. elegans, and D. melanogaster, reveals two different AA/TT periodical patterns associated with the nucleosome positioning. The first pattern is the counter-phase oscillation of AA and TT dinucleotides, which has been frequently considered as the nucleosome DNA pattern. This represents the sequence rule I for chromatin structure. The second pattern is the in-phase oscillation of the AA and TT dinucleotides with the same nucleosome DNA period, 10.4 bases. This pattern apparently corresponds to curved DNA, that also participates in the nucleosome formation, and represents the sequence rule II for chromatin. The positional correlations of AA and TT dinucleotides also indicate that the nucleosomes are separated by specific linker sizes (preferably 8, 18,…bases), dictated by the steric exclusion rules. Thus, the sequence positions of the neighboring nucleosomes are correlated, and this represents the sequence rule III.     

Methylenetetrahydrofolate reductase genotypes and haplotypes associated with susceptibility to colorectal cancer in an eastern Chinese Han population

Methylenetetrahydrofolate reductase (MTHFR) plays an important role in folate metabolism and is involved in DNA synthesis, DNA repair and DNA methylation. The two common functional polymorphisms of MTHFR, C677T and A1298C have been associated with several diseases, including cancer. We made a case-control study to analyze a possible association of MTHFR gene polymorphisms C677T and A1298C with risk for colorectal cancer in an eastern Chinese Han population of 137 patients with a confirmed histopathological diagnosis of CRC and 145 age- and gender-matched controls with no history of cancer. DNA was isolated from peripheral blood samples and the genotypes were determined by PCR-RFLP. The concentrations of folate in plasma were measured by chemiluminescence immunoassay. The MTHFR 677TT genotype had a protective effect against colorectal cancer, with an odds ratio (OR) = 0.467 (95% confidence interval (CI) = 0.225-0.966). The 1298CC genotype was significantly correlated with a reduced risk of colorectal cancer (OR = 0.192; 95%CI = 0.040-0.916). Compared with the MTHFR 677CC and MTHFR 1298 AA genotypes, for individuals who carried both MTHFR 677CC and 1298CC genotypes, the OR of colorectal cancer was 0.103 (95%CI = 0.012-0.900); among individuals who carried both MTHFR 677TT and 1298AC genotypes, the OR for risk of colorectal cancer was 0.169 (95%CI = 0.044-0.654). MTHFR 677TT+CT genotypes had a significantly lower plasma folate concentration than those with the MTHFR 677CC genotype. MTHFR 1298AC+CC genotypes had a lower plasma folate concentration than those with the MTHFR 1298AA genotype (P < 0.05). In conclusion, subjects with the MTHFR 677TT and MTHFR 1298CC genotypes appeared to have a significantly lower risk for colorectal cancer. MTHFR haplotypes 677CC/1298CC and 677TT/1298AC were less common in cases than in controls. These haplotypes, when compared to the most common haplotype 677CC/1298AA, were associated with a decreased risk for colorectal cancer. We conclude that plasma folate level is influenced by MTHFR genotypes.     

DNase I susceptibility of bent DNA and its alteration by ditercalinium and distamycin

The bending of kinetoplast DNA from Crithidia fasciculata is thought to be related to the periodic distribution of AA or TT cluster sequences. The sensitivity to DNase I of the two strands of this DNA was analyzed at nucleotide resolution by sequencing gel electrophoresis. The effect on the DNase I cleavage pattern of two drugs, ditercalinium and distamycin, that are able to remove bending was analyzed. The same analysis was done on a pBR 322 DNA fragment of random sequence as a control. The periodic distribution of the AA or TT clusters in the bent DNA fragment was first analyzed by computing the autocorrelation function of the AA or TT clusters in the bent DNA fragment. It is shown that the AT tracts are on average 10.5 base pairs apart. This value is almost identical with that of the B-DNA helix pitch in solution [10.5 (Wang, 1979); 10.6 +/- 0.1 (Rhodes & Klug, 1980)]. To reveal the periodic pattern of DNase I cleavage on this bent DNA, alone or in presence of drugs, the cross correlation between the different bands obtained from DNAse I cleavage and the presence of AA or TT sequences was computed. This shows that GC and mixed sequences are the most sensitive regions. These data also suggest that there is a periodic fluctuation in the width of the minor groove in the bent fragment. Ditercalinium and distamycin alter the DNase I cutting pattern of the bent DNA fragment but in an inverse fashion.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structural characterization of intrinsically curved AT-rich DNA sequences.

AT-rich DNA sequences other than AnTm tracts (n + m > or = 4) are known to be intrinsically curved. The AATAT-element constitutes one known example of these sequences. In this paper, the elucidation of the structural basis of the curvature induced by this sequence element was addressed. As judged by the patterns of cleavage by the hydroxyl radical and DNase I, the AATAT sequence shows a narrow minor groove. Furthermore, the 5' adenine residue of the AA dinucleotide contained within the sequence is hyperreactive to diethylpyrocarbonate. Similar structural properties are shown by several sequences inducing intrinsic DNA curvature, such as an A5-tract or the closely related ATAAT, AATATA and TAATAT sequences, which are also shown here to induce curvature. On the other hand, other related sequences, such as TATAA and ATATA, that do not induce curvature, show different structural characteristics.     

Sequence structure of human nucleosome DNA

Positional distributions of various dinucleotides in experimentally derived human nucleosome DNA sequences are analyzed. Nucleosome positioning in this species is found to depend largely on GG and CC dinucleotides periodically distributed along the nucleosome DNA sequence, with the period of 10.4 bases. The GG and CC dinucleotides oscillate counterphase, i.e., their respective preferred positions are shifted about a half-period from one another, as it was observed earlier for AA and TT dinucleotides. Other purine-purine and pyrimidine-pyrimidine dinucleotides (RR and YY) display the same periodical and counterphase pattern. The dominance of oscillating GG and CC dinucleotides in human nucleosomes and the contribution of AG(CT), GA(TC), and AA(TT) suggest a general nucleosome DNA sequence pattern - counterphase oscillation of RR and YY dinucleotides. AA and TT dinucleotides, commonly accepted as major players, are only weak contributors in the case of human nucleosomes.     

Evidence for opposite groove-directed curvature of GGGCCC and AAAAA sequence elements.

The repetitive sequence (AGGGCCCTAGAGGGGCCC-TAG)n was previously shown to be curved by gel mobility assays. Here we show, using hydroxy radical/DNase I digestion and differential helical phasing experiments that the curvature is directed towards the major groove and is located in the GGGCCC, but not the CTAGAG segments. The effect of the GC step in the context of the GGGCCC motif is apparently about as large as that of AA/TT, i.e. enough to cancel the macroscopic curvature of helically phased A-tracts. These data are in agreement with positive roll-like curvature of the GCC/GGC motif, predicted from nucleosome packing data and the 3D structure of the GGGGCCCC octamer, but they are not in agreement with the dinucleotide-based roll angle values predicted for AG/CT, TA, GG/CC and GC steps. Our results thus indicate the importance of interactions beyond the dinucleotide steps in predictive models of DNA curvature.     

Preferred positions of AA and TT dinucleotides in aligned nucleosomal DNA sequences.

Multiple alignment of 118 nucleosomal DNA sequences by maximizing simultaneously match of AA dinucleotides and match of TT dinucleotides results in a pattern of the dinucleotide distributions which is characteristic of the nucleosomal DNA sequences. The AA dinucleotides are found to be distributed symmetrically relative to the TT dinucleotide distribution, around the middle point of the nucleosomal DNA sequence. The distances between major peaks of the distributions are multiples of about 10.4 bases. The peaks of the TT distribution are shifted by 6 bases downstream from the peaks of the AA distribution.     

Stability of DNA duplexes containing GG, CC, AA, and TT mismatches

We employed salt-dependent differential scanning calorimetric measurements to characterize the stability of six oligomeric DNA duplexes (5'-GCCGGAXTGCCGG-3'/5'-CCGGCAYTCCGGC-3') that contain in the central XY position the GC, AT, GG, CC, AA, or TT base pair. The heat-induced helix-to-coil transitions of all the duplexes are associated with positive changes in heat capacity, DeltaC(p), ranging from 0.43 to 0.53 kcal/mol. Positive values of DeltaC(p) result in strong temperature dependences of changes in enthalpy, DeltaH degrees, and entropy, DeltaS degrees , accompanying duplex melting and cause melting free energies, DeltaG degrees, to exhibit characteristically curved shapes. These observations suggest that DeltaC(p) needs to be carefully taken into account when the parameters of duplex stability are extrapolated to temperatures distant from the transition temperature, T(M). Comparison of the calorimetric and van't Hoff enthalpies revealed that none of the duplexes studied in this work exhibits two-state melting. Within the context of the central AXT/TYA triplet, the thermal and thermodynamic stabilities of the duplexes in question change in the following order: GC > AT > GG > AA approximately TT > CC. Our estimates revealed that the thermodynamic impact of the GG, AA, and TT mismatches is confined within the central triplet. In contrast, the thermodynamic impact of the CC mismatch propagates into the adjacent helix domains and may involve 7-9 bp. We discuss implications of our results for understanding the origins of initial recognition of mismatched DNA sites by enzymes of the DNA repair machinery.     

Propeller-Twisting of Base-pairs and the Conformational Mobility of Dinucleotide Steps in DNA

When DNA is bent around a protein, it must distort. The distortion occurs by changes in the conformation of successive dinucleotide steps. Bending does not necessarily occur uniformly: some steps might remain particularly rigid, i.e. they might deform relatively little, while others might take more than their proportional share of deformation. We investigate here the deformational capacity of specific dinucleotide steps by examining a database of crystallized oligomers. Dividing the steps into ten types by sequence (AA(=TT), AC(=GT), AG(=CT), AT, CA(=TG), CG, GA(=TC), GC, GG(=CC) and TA), we find that some step types are practically rigid, while others have considerable internal mobility or conformational flexibility. Now in general base-pairs are not planar, but have Propeller-Twist. We find a clear empirical correlation between the level of Propeller-Twist in the base-pairs and the flexibility of the dinucleotide step which they constitute. Propeller-Twist in the base-pairs makes stacking into a dinucleotide step more awkward than in plane base-pairs. In particular, it provides a stereochemical “locking” effect which can make steps with highly Propeller-Twisted base-pairs rigid. Although the origins of Propeller-Twist are not yet clearly understood, this result provides a key to understanding the flexibility of DNA in bending around proteins.     

Structural basis of stable bending in DNA containing An tracts. Different types of bending

Structural determinants of DNA bending of different types have been studied by theoretical conformational analysis of duplexes. Their terminal parts were fixed either in an ordinary low-energy B-like conformation or in "anomalous" conformations with a narrowed minor groove typical of An tracts. The anomalous conformations had different negative tilt angles (up to about zero), different propeller twists and minor groove widths. Calculations have been performed for DNA fragments AnTm, TnAm, AnGCTm, AnCGTm, TmGCAn, TmCGAn which are the models of the junction of two anomalous structures on An and Tm tracts. On the AT step of the AnTm fragment the minor groove can be easily narrowed so that a whole unbent fragment of anomalous structure is formed on AnTm. According to our energy estimates, there should not be any reliable bending on AnTm. In contrast, in all other cases there was a pronounced roll-like bending into the major groove in the chemical symmetry region. Calculations of the junction between the anomalous and ordinary B-like structure for GnTm and CnAm have shown that there is an equilibrium bending with a tilt component towards the chain having the anomalous structure at the 5'-end. From our calculations it is impossible to determine precisely the direction of bending, though it can be suggested that the roll component of bending might be directed towards the major groove. The anomalous structure is the main reason of bending; alternations of pyrimidines and purines can modulate the value and the direction of equilibrium bending (only the value in the case of self-complementary fragments).(ABSTRACT TRUNCATED AT 250 WORDS)     

Sequence motifs and free energies of selected natural and non-natural nucleosome positioning DNA sequences.

Our laboratories recently completed SELEX experiments to isolate DNA sequences that most-strongly favor or disfavor nucleosome formation and positioning, from the entire mouse genome or from even more diverse pools of chemically synthetic random sequence DNA. Here we directly compare these selected natural and non-natural sequences. We find that the strongest natural positioning sequences have affinities for histone binding and nucleosome formation that are sixfold or more lower than those possessed by many of the selected non-natural sequences. We conclude that even the highest-affinity sequence regions of eukaryotic genomes are not evolved for the highest affinity or nucleosome positioning power. Fourier transform calculations on the selected natural sequences reveal a special significance for nucleosome positioning of a motif consisting of approximately 10 bp periodic placement of TA dinucleotide steps. Contributions to histone binding and nucleosome formation from periodic TA steps are more significant than those from other periodic steps such as AA (=TT), CC (=GG) and more important than those from the other YR steps (CA (=TG) and CG), which are reported to have greater conformational flexibility in protein-DNA complexes even than TA. We report the development of improved procedures for measuring the free energies of even stronger positioning sequences that may be isolated in the future, and show that when the favorable free energy of histone-DNA interactions becomes sufficiently large, measurements based on the widely used exchange method become unreliable.