Anomalous structure and properties of poly (dA).poly(dT). Computer simulation of the polynucleotide structure with the spine of hydration in the minor groove.

The results of the search for low-energy conformations of poly(dA).poly(dT) and of the poly(dA).poly(dT) "complex" with the spine of hydration similar to that found by Dickerson and co-workers (Kopka, M.L., Fratini, A.V., Drew, H.R. and Dickerson, R.E. (1983) J. Mol. Biol. 163, 129-146) in the minor groove of the CGCGAATTCGCG crystals are described. It is shown that the existence of such a spine in the minor groove of poly(dA).poly(dT) is energetically favourable. Moreover, the spine of hydration makes the polynucleotide conformation similar to the poly(dA).poly(dT) structure in fibers and to the conformation of the central part of CGCGAATTCGCG in crystals; it also acquires features characteristic of the structure of poly(dA).poly(dT) and DNA oligo(dA)-tracts in solution. It is shown that the existence of the TpA step in conformations characteristic of the poly(dA).poly(dT) complex with the spine of hydration is energetically unfavourable (in contrast to the ApT step) and therefore this step should result in destabilization of the spine of hydration in the DNA minor groove. Thus, it appears that the spine of hydration as described by Dickerson and co-workers is unlikely to exist in the poly d(A-T).poly d(A-T) structure. The data obtained permit us to interpret a large body of experimental facts concerning the unusual structure and properties of poly(dA).poly(dT) and oligo(dA)-tracts in DNA both in fibers and in solution. The results provide evidence of the existence of the minor groove spine of hydration both in fibers and in solution on A/T tracts of DNA which do not contain the TpA step. The spine plays an active role in the formation of the anomalous conformation of these tracts.     

Experimental study of DNA hydratation

The present article as a review of experimental investigations of water-DNA interaction in solutions, films, fibres, and crystals. The discussion of the experimental data is directed to the analysis of the sizes of the hydration shells of the DNA atomic groups, the distances at which atomic groups of DNA affect the hydration of each other; the structural, thermodynamic and kinetic characteristics of the water in the DNA hydration shell and their relations with DNA structure. The modern views on the mechanisms by means of which the water affects the physico-chemical properties of DNA are discussed.     

Sequence Dependencies of DNA Deformability and Hydration in the Minor Groove

DNA deformability and hydration are both sequence-dependent and are essential in specific DNA sequence recognition by proteins. However, the relationship between the two is not well understood. Here, systematic molecular dynamics simulations of 136 DNA sequences that differ from each other in their central tetramer revealed that sequence dependence of hydration is clearly correlated with that of deformability. We show that this correlation can be illustrated by four typical cases. Most rigid basepair steps are highly likely to form an ordered hydration pattern composed of one water molecule forming a bridge between the bases of distinct strands, but a few exceptions favor another ordered hydration composed of two water molecules forming such a bridge. Steps with medium deformability can display both of these hydration patterns with frequent transition. Highly flexible steps do not have any stable hydration pattern. A detailed picture of this correlation demonstrates that motions of hydration water molecules and DNA bases are tightly coupled with each other at the atomic level. These results contribute to our understanding of the entropic contribution from water molecules in protein or drug binding and could be applied for the purpose of predicting binding sites.     

Dependence of the hydration shell structure in the minor groove of the DNA double helix on the groove width as revealed by Monte Carlo simulation.

The hydration shell of several conformations of the polynucleotides poly(dA).poly(dT), poly(dA).poly(dU), and poly(dA-dI).poly(dT-dC) has been simulated using the Monte Carlo method (Metropolis sampling). Calculations have shown that the structure of the hydration shell of the minor groove greatly depends on its width. In conformations with a narrowed minor groove, the first layer of the hydration shell of this groove has only one molecule per nucleotide pair that forms H bonds with purine N3 of one pair and pyrimidine O2 of the next pair. The second layer of the hydration shell of such conformations contains molecules that form H bonds between two adjacent molecules of the first layer. The probability of formation of hydration spine is about 20% while the bridges of the first layer are formed with a probability of about 70%. In the first layer of the minor groove of the B-DNA conformation with wide minor groove there are approximately two water molecules per base pair that form H bonds with purine N3 or pyrimidine O2 and with the sugar ring oxygen of the adjacent nucleotide. The probability of simultaneous H bonding of a water molecule with N3 (or O2) and O of sugar ring is about 30%. The results of simulation suggest that hydration spine proposed for the narrowed minor groove of oligonucleotide crystals [H. R. Drew, and R. E. Dickerson (1981) Journal of Molecular Biology, Vol. 151, pp. 535-556] can be formed in fibers of poly(dA).poly(dT), poly(dA).poly(dU), and poly(dA-dI).poly(dT-dC) as well as in DNA fragments of these sequences in solution.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular dynamics simulations of A. T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A. T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na(+) intrusions in the first layer of hydration of the minor groove, with long residence times, up to approximately 3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove.     

NMR observation of individual molecules of hydration water bound to DNA duplexes: direct evidence for a spine of hydration water present in aqueous solution.

The residence times of individual hydration water molecules in the major and minor grooves of DNA were measured by nuclear magnetic resonance (NMR) spectroscopy in aqueous solutions of d-(CGCGAATTCGCG)2 and d-(AAAAATTTTT)2. The experimental observations were nuclear Overhauser effects (NOE) between water protons and the protons of the DNA. The positive sign of NOEs with the thymine methyl groups shows that the residence times of the hydration water molecules near these protons in the major groove of the DNA must be shorter than about 500 ps, which coincides with the behavior of surface hydration water in peptides and proteins. Negative NOEs were observed with the hydrogen atoms in position 2 of adenine in both duplexes studied. This indicates that a 'spine of hydration' in the minor groove, as observed by X-ray diffraction in DNA crystals, is present also in solution, with residence times significantly longer than 1 ns. Such residence times are reminiscent of 'interior' hydration water molecules in globular proteins, which are an integral part of the molecular architecture both in solution and in crystals.     

Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.     

On the efficiency of hole and electron transfer from the hydration layer to DNA: An EPR study of crystalline DNA X-irradiated at 4 K

The aim of this project was to gain an improved understanding of how the efficiency of hole and electron transfer from the solvation layer to DNA decreases as a function of distance from DNA. The packing of DNA in crystals of known structure makes it possible to calculate the degree of DNA hydration with a precision that is significantly greater than that achievable for amorphous samples. Previous work on oligodeoxynucleotide crystals has demonstrated that the efficiency of free radical trapping by DNA exposed to ionizing radiation at 4 K is relatively insensitive to base sequence, conformation, counterion, or base stacking continuity. Having eliminated these confounding variables, it is now possible to ascertain the degree of radical transfer that occurs from ionized water as a function of DNA hydration (Gamma, in mol water/mol nucleotide). EPR is used to measure the hydroxyl radical concentration in crystals irradiated at 4 K. From a lack of hydroxyl radicals trapped in the inner hydration mantle, we determine that hole transfer to DNA is complete for water molecules located within 8 A. This corresponds to Gamma = 9-11 and indicates that hole transfer is 100% (as efficient as direct ionization of DNA) for water molecules adjacent to DNA. Beyond approximately 8 A (Gamma > 10), hydroxyl radicals are observed; thus deprotonation of the water radical cation is seen to compete with hole transfer to DNA as soon as one water intervenes between the ionized water and DNA. The boundary for 0% hole transfer is projected to occur somewhere between 15 and 20 waters per nucleotide. Electron transfer, on the other hand, is 100% efficient across the entire range studied, 4.2 相似文献   

The hydration of nucleic acid duplexes as assessed by a combination of volumetric and structural techniques

Using high precision densimetric and ultrasonic measurements, we have determined, at 25°C, the apparent molar volumes ΦV and the apparent molar compressibilities ΦK_S of four nucleic acid duplexes—namely, the DNA duplex, poly(dIdC)poly(dIdC); the RNA duplex, poly(rA)poly(rU); and the two DNA/RNA hybrid duplexes, poly(rA)poly(dT) and poly(dA)poly(rU). Using available fiber diffraction data on these duplexes, we have calculated the molecular volumes as well as the solvent‐accessible surface areas of the constituent charged, polar, and nonpolar atomic groups. We found that the hydration properties of these nucleic acid duplexes do not correlate with the extent and the chemical nature of the solvent‐exposed surfaces, thereby suggesting a more specific set of duplex–water interactions beyond general solvation effects. A comparative analysis of our volumetric data on the four duplexes, in conjunction with available structural information, suggests the following features of duplex hydration: (a) The four duplexes exhibit different degrees of hydration, in the order poly(dIdC)poly(dIdC) > poly(dGdC)poly(dGdC) > poly(dAdT)poly(dAdT) ≈ poly(dA)poly(dT). (b) Repetitive AT and IC sequences within a duplex are solvated beyond general effects by a spine of hydration in the minor groove, with this sequence‐specific water network involving about 8 additional water molecules from the second and, perhaps, even the third hydration layers. (c) Repetitive GC and IC sequences within a duplex are solvated beyond general effects by a “patch of hydration” in the major groove, with this water network involving about 13 additional water molecules from the second and, perhaps, even the third hydration layers. (d) Random sequence, polymeric DNA duplexes, which statistically lack extended regions of repetitive AT, GC, or IC sequences, do not experience such specific enhancements of hydration. Consequently, consistent with our previous observations (T. V. Chalikian, A. P. Sarvazyan, G. E. Plum, and K. J. Breslauer, Biochemistry, 1994, Vol. 33, pp. 2394–2401), duplexes with approximately 50% AT content exhibit the weakest hydration, while an increase or decrease from this AT content causes enhancement of hydration, either due to stronger hydration of the minor groove (an increase in AT content) or due to stronger hydration of the major groove (an increase in GC content). (e) In dilute aqueous solutions, a B‐DNA duplex is more hydrated than an A‐DNA duplex, a volumetric‐based conclusion that is in agreement with previous results obtained on crystals, fibers, and DNA solutions in organic solvent–water mixtures. (f) the A‐like, RNA duplex poly(rA)poly(rU) and the structurally similar A‐like, hybrid duplex poly(rA)poly(dT), exhibit similar hydration properties, while the structurally distinct A‐like, hybrid duplex poly(rA)poly(dT) and non‐A‐like, hybrid duplex poly(dA)poly(rU) exhibit differential hydration properties, consistent with structural features dictating hydration characteristics. We discuss how volumetric characterizations, in conjunction with structural studies, can be used to describe, define, and resolve the general and sequence/conformation‐specific hydration properties of nucleic acid duplexes. © 1999 John Wiley & Sons, Inc. Biopoly 50: 459–471, 1999     

DNA curvature does not require bifurcated hydrogen bonds or pyrimidine methyl groups.

Short tracts of the homopolymer dA.dT confer intrinsic curvature on the axis of the DNA double helix. This phenomenon is assumed to be a consequence of such tracts adopting a stable B'-DNA conformation that is distinct from B-form structure normally assumed by other DNA sequences. The more stable B' structure of dA.dT tracts has been attributed to several possible stabilizing factors: (1) optimal base stacking interactions consequent upon the high propeller twist, (2) bifurcated hydrogen bonds between adjacent dA.dT base-pairs, (3) stacking interactions involving the dT methyl groups, and finally (4) a putative spine of ordered water molecules in the minor groove. DNA oligodeoxynucleotides have been synthesized that enable these hypotheses to be tested; of particular interest is the combination of effects due to bifurcation (2) and methylation of the pyrimidines nucleotides (3). The data indicate that neither bifurcated hydrogen bonds nor pyrimidine methyl groups nor both are essential for DNA curvature. The data further suggest that the influence of the minor groove spine of hydration on the B'-formation is small. The experiments favor the hypothesis that base stacking interactions are the dominant force in stabilizing the B'-form structure.     

Molecular dynamics simulation study of oriented polyamine- and Na-DNA: sequence specific interactions and effects on DNA structure

Molecular dynamics (MD) computer simulations have been carried out on four systems that correspond to an infinite array of parallel ordered B-DNA, mimicking the state in oriented DNA fibers and also being relevant for crystals of B-DNA oligonucleotides. The systems were all comprised of a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charges. The sequence of the double helical DNA decamer was d(5'-ATGCAGTCAG)xd(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine(3+) (Spd(3+)) and putrescine(2+) (Put(2+)), the synthetic polyamine diaminopropane(2+) (DAP(2+)), and the simple monovalent cation Na(+). This work compares the specific structures of the polyamine- and Na-DNA systems and how they are affected by counterion interactions. It also describes sequence-specific hydration and interaction of the cations with DNA. The local DNA structure is dependent on the nature of the counterion. Even the very similar polyamines, Put(2+) and DAP(2+), show clear differences in binding to DNA and in effect on hydration and local structure. Generally, the polyamines disorder the hydration of the DNA around their binding sites whereas Na(+) being bound to DNA attracts and organizes water in its vicinity. Cation binding at the selected sites in the minor and in the major groove is compared for the different polyamines and Na(+). We conclude that the synthetic polyamine (DAP(2+)) binds specifically to several structural and sequence-specific motifs on B-DNA, unlike the natural polyamines, Spd(3+) and Put(2+). This specificity of DAP(2+) compared to the more dynamic behavior of Spd(3+) and Put(2+) may explain why the latter polyamines are naturally occurring in cells.     

Mechanism of DNA Recognition by the Restriction Enzyme EcoRV

EcoRV, a restriction enzyme in Escherichia coli, destroys invading foreign DNA by cleaving it at the center step of a GATATC sequence. In the EcoRV-cognate DNA crystallographic complex, a sharp kink of 50° has been found at the center base-pair step (TA). Here, we examine the interplay between the intrinsic propensity of the cognate sequence to kink and the induction by the enzyme by performing all-atom molecular dynamics simulations of EcoRV unbound and interacting with three DNA sequences: the cognate sequence, GATATC (TA); the non-cognate sequence, GAATTC (AT); and with the cognate sequence methylated on the first adenine GA_CH3TATC (TA-CH₃). In the unbound EcoRV, the cleft between the two C-terminal subdomains is found to be open. Binding to AT narrows the cleft and forms a partially bound state. However, the intrinsic bending propensity of AT is insufficient to allow tight binding. In contrast, the cognate TA sequence is easier to bend, allowing specific, high-occupancy hydrogen bonds to form in the complex. The absence of cleavage for this methylated sequence is found to arise from the loss of specific hydrogen bonds between the first adenine of the recognition sequence and Asn185. On the basis of the results, we suggest a three-step recognition mechanism. In the first step, EcoRV, in an open conformation, binds to the DNA at a random sequence and slides along it. In the second step, when the two outer base pairs, GAxxTC, are recognized, the R loops of the protein become more ordered, forming strong hydrogen-bonding interactions, resulting in a partially bound EcoRV-DNA complex. In the third step, the flexibility of the center base pair is probed, and in the case of the full cognate sequence the DNA bends, the complex strengthens and the protein and DNA interact more closely, allowing cleavage.     

Overview of the structure of all-AT oligonucleotides: organization in helices and packing interactions

We present the crystalline organization of 33 all-AT deoxyoligonucleotide duplexes, studied by x-ray diffraction. Most of them have very similar structures, with Watson-Crick basepairs and a standard average twist close to 36 degrees. The molecules are organized as parallel columns of stacked duplexes in a helical arrangement. Such organization of duplexes is very regular and repetitive: all sequences show the same pattern. It is mainly determined by the stacking of the terminal basepairs, so that the twist in the virtual TA base step between neighbor duplexes is always negative, approximately -22 degrees. The distance between the axes of parallel columns is practically identical in all cases, approximately 26 A. Interestingly, it coincides with that found in DNA viruses and fibers in their hexagonal phase. It appears to be a characteristic distance for ordered parallel DNA molecules. This feature is due to the absence of short range intermolecular forces, which are usually due to the presence of CG basepairs at the end of the oligonucleotide sequence. The duplexes apparently interact only through their diffuse ionic atmospheres. The results obtained can thus be considered as intermediate between liquid crystals, fibers, and standard crystal structures. They provide new information on medium range DNA-DNA interactions.     

Sequence Recognition of DNA by Protein-Induced Conformational Transitions

The binding of proteins to specific sequences of DNA is an important feature of virtually all DNA transactions. Proteins recognize specific DNA sequences using both direct readout (sensing types and positions of DNA functional groups) and indirect readout (sensing DNA conformation and deformability). Previously we showed that the P22 c2 repressor N-terminal domain (P22R NTD) forces the central non-contacted 5'-ATAT-3' sequence of the DNA operator into the B′ state, a state known to affect DNA hydration, rigidity and bending. Usually the B′ state, with a narrow minor groove and a spine of hydration, is reserved for A-tract DNA (TpA steps disrupt A-tracts). Here, we have co-crystallized P22R NTD with an operator containing a central 5′-ACGT-3′ sequence in the non-contacted region. C·G base pairs have not previously been observed in the B′ state and are thought to prevent it. However, P22R NTD induces a narrow minor groove and a spine of hydration to 5'-ACGT-3'. We observe that C·G base pairs have distinctive destabilizing and disordering effects on the spine of hydration. It appears that the reduced stability of the spine results in a higher energy cost for the B to B′ transition. The differential effect of DNA sequence on the barrier to this transition allows the protein to sense the non-contacted DNA sequence.     

Isopycnic sedimentation of DNA in metrizamide: the effect of low concentrations of ions on buoyant density and hydration

Metrizamide, an inert, non-ionic organic compound, dissolves in water to give a dense solution in which DNA bands isopycnically at a density corresponding to that of fully hydrated DNA. Density-gradient centrifugation in solutions of metrizamide has been used to determine the effects of very dilute solutions of salts on the buoyant density of native and denatured DNA. It has been shown that the buoyant density of DNA is dependent on both the counter-cation and the anion present. Interpretation of the data in terms of the degree of hydration of the macromolecule indicates that (i), NaDNA is more highly hydrated than CsDNA; and (ii), the hydration of NaDNA varies with anion in the order sulphate< fluoride< chloride< bromide< iodide.     

Hydration of B-DNA: comparison between the water network around poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC) on the basis of Monte Carlo computations

A computational method is elaborated for studying the water environment around regular polynucleotide duplexes; it allows rigorous structural information on the hydration shell of DNA to be obtained. The crucial aspect of this Monte Carlo simulation is the use of periodical boundary conditions. The output data consists of local maxima of water density in the space near the DNA molecule and the properties of one- and two-membered water bridges as function of pairs of polar groups of DNA. In the present paper the results for poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC) are presented. The differences in their hydration shells are of a purely structural nature and are caused by the symmetry of the polar groups of the polymers under study, the symmetry being reflected by the hydration shell. The homopolymer duplex hydration shell mirrors the mononucleotide repeat. The water molecules contacting the polynucleotide in the minor groove are located nearly in the plane midway between the planes of successive base pairs. One water molecule per base pair forms a water bridge facing two polar groups of bases from adjacent base pairs and on different strands making a "spine"-like structure. In contrast, the major groove hydration is stabilized exclusively by two-membered water bridges; the water molecules deepest in the groove are concentrated near the plane of the corresponding base pair. The alternating polymer is characterized by a marked dyad symmetry of the hydration shell corresponding to the axis between two successive base pairs. The minor groove hydration of the dCpdG step resembles the characteristic features of the homopolymer, but the bridge between the O2 oxygens of the other base-stacking type is formed by two water molecules. The major groove hydration is characterized by high probability of one-membered water bridges and by localization of a water molecule on the dyad axis of the dGpdC step. The found structural elements are discussed as reasonable invariants of a dynamic hydration shell.     

Theory of acoustic mode vibrations of DNA fibers

A previous model for acoustic mode vibrations of a DNA molecule in water is extended to the case of an array of many DNA molecules, as occurs in the fibers studied in most experimental work on DNA. The acoustic modes of this system are found to consist of coupled modes of water sound vibrations and DNA acoustic modes. This model is used to study the electrostatic coupling of acoustic vibrations to the relaxational modes of the orientational degrees of freedom of the water molecules. It is found that the long-range or macroscopic electric field generated by the acoustic mode vibrations of the water-DNA system gives too small a damping and frequency shift of the acoustic modes to account for the observations on DNA fibers. Therefore, the observed damping and frequency shifts are most likely due to either friction between the surrounding water and the vibrating DNA, or coupling to the water orientation degrees of freedom resulting from the short range (i.e., screened) Coulomb interaction. The latter explanation (which is most likely the correct one) implies that the relaxation time of the hydration shell water is longer than the observed relaxation time by a factor of the static dielectric constant of the hydration water.     

Complex between triple helix of collagen and double helix of DNA in aqueous solution

We demonstrate in this paper that one example of a biologically important and molecular self-assembling complex system is a collagen–DNA ordered aggregate which spontaneously forms in aqueous solutions. Interaction between the collagen and the DNA leads to destruction of the hydration shell of the triple helix and stabilization of the double helix structure. From a molecular biology point of view this nano-scale self-assembling superstructure could increase the stability of DNA against the nucleases during collagen diseases and the growth of collagen fibrills in the presence of DNA.     

Complicated water orientations in the minor groove of the B-DNA decamer d(CCATTAATGG)2 observed by neutron diffraction measurements

It has long been suspected that the structure and function of a DNA duplex can be strongly dependent on its degree of hydration. By neutron diffraction experiments, we have succeeded in determining most of the hydrogen (H) and deuterium (D) atomic positions in the decameric d(CCATTAATGG)₂ duplex. Moreover, the D positions in 27 D₂O molecules have been determined. In particular, the complex water network in the minor groove has been observed in detail. By a combined structural analysis using 2.0 Å resolution X-ray and 3.0 Å resolution neutron data, it is clear that the spine of hydration is built up, not only by a simple hexagonal hydration pattern (as reported in earlier X-ray studies), but also by many other water bridges hydrogen-bonded to the DNA strands. The complexity of the hydration pattern in the minor groove is derived from an extraordinary variety of orientations displayed by the water molecules.