Structure of the hydration shells of oligo(dA-dT).oligo(dA-dT) and oligo(dA).oligo(dT) tracts in B-type conformation on the basis of Monte Carlo calculations

Attractive hydration forces in DNA-dendrimer interactions on the nanometer scale

The energetic contribution of attractive hydration forces arising from water ordering is an interesting but often neglected aspect of macromolecular interactions. Ordering effects of water can bring about cooperativity in many intermolecular transactions, in both the short and long range. Given its high charge density, this is of particular importance for DNA. For instance, in nanotechnology, highly charged dendrimers are used for DNA compaction and transfection. Hypothesizing that water ordering and hydration forces should be maximal for DNA complexes that show charge complementarity (positive-negative), we present here an analysis of water ordering from molecular dynamics simulations and free energy calculations of the interaction between DNA and a nanoparticle with a high positive charge density. Our results indicate not only that complexation of the dendrimer with DNA affects the local water structure but also that ordered water molecules facilitate long-range interactions between the molecules. This contributes significantly to the free energy of binding of dendrimers to DNA and extends the interaction well beyond the electrostatic range of the DNA. Such water effects are of potentially substantial importance in cases when molecules appear to recognize each other across sizable distances, or for which kinetic rates are too fast to be due to pure diffusion. Our results are in good agreement with experiments on the role of solvent in DNA condensation by multivalent cations and exemplify a microscopic realization of mean-field phenomenological theories for hydration forces between mesoscopic surfaces.

Adiabatic Electron Affinities of the Polyhydrated Adenine−Thymine Base Pair: A Density Functional Study

Adiabatic electron affinities (AEAs) of the adenine-thymine (AT) base pair surrounded by 5 and 13 water molecules have been studied by density functional theory (DFT). Geometries of neutral AT x nH2O and anionic (AT x nH2O)- complexes (n = 5 and 13) were fully optimized, and vibrational frequency analysis was performed at the B3LYP/6-31+G** level of theory. The optimized structures of the neutral (AT x nH2O) and (AT x nH2O)- complexes were found to be somewhat nonplanar. Some of the water molecules are displaced away from the AT ring plane and linked with one another by hydrogen bonds. The optimized structures of the complexes are found to be in a satisfactory agreement with the observed experimental and molecular dynamics simulation results. In the optimized anionic complexes, the thymine (T) moiety was found to be puckered, whereas the adenine (A) moiety remained almost planar. Natural population analysis (NPA) performed using the B3LYP/6-31+G** method shows that the thymine moiety in the anionic (AT x nH2O)- complexes (n = 5 and 13) has most of the excess electronic charge, i.e., approximately -0.87 and approximately -0.81 (in the unit of magnitude of the electronic charge), respectively. The zero-point energy corrected adiabatic electron affinities of the hydrated AT base pair were found to be positive both for n = 5 and 13 and have the values of 0.97 and 0.92 eV, respectively, which are almost three times the AEA of the AT base pair. The results show that the presence of water molecules appreciably enhances the EA of the base pair.

Comparative studies of high resolution Z-DNA crystal structures. Part 1: Common hydration patterns of alternating dC-dG

The water structure in three crystal forms of the left-handed Z-DNA hexamer [d(CGCGCG)]2 has been analyzed. Several common motifs have been found in the first hydration shells. On the convex surface, the major groove of the left-handed conformation, water molecules bridge the guanine O-6 keto groups at GpC steps. Cytosine N-4 nitrogens of opposite strands are hydrated by tandem water molecules. At the bottom of the minor groove, a string of water molecules connects the cytosine O-2 keto groups. Across the minor groove guanine N-2 nitrogens are bridged to phosphate oxygens of cytosine and guanine residues by one or two water molecules. In contrast to the very regular geometry of the water structure around the bases, the arrangement of water molecules between phosphate groups appears to be less ordered. However, there is a strong correlation between the interphosphate distances and the number of water molecules or ions which link the phosphate groups. In all three structures various ions, such as sodium and magnesium ions, as well as the protonated amino and imino groups of the polycation spermine displace and replace water molecules in the first hydration shell. Nevertheless, the analysis reveals that numerous first hydration shell water molecules in Z-DNA crystals can be regarded as part of the DNA structure. Their positions and thermal parameters are generally independent of changes in the local crystallographic environment.

Molecular dynamics simulations of poly(dA).poly(dT): comparisons between implicit and explicit solvent representations

The program AMBER 3.0 has been used to generate molecular dynamics trajectories of a poly(dA).poly(dT) decamer. The simulations were performed using different methods to treat solvent effects. Results of a simulation including 18 counterions NH4+ and 4109 water molecules under (N, P, T) conditions were compared to simulation runs with implicit solvent representation in which solvent screening effects were represented by the use of a sigmoidal distance-dependent dielectric function. In the latter case, the system was simulated under microcanonical (N, V, E) and canonical (N, V, T) conditions. For the fully hydrated system simulation, a preequilibration protocol was developed since it was observed that long and progressive periods of heating and equilibration on the overall system were necessary in order to avoid energetic collisions between the solute and the solvent molecules, leading to severe irreversible deformation of the solute. A detailed analysis of DNA conformations, sugar puckers, and stability of the hydrogen bonds, Watson-Crick and three-center H bonds, is reported. The results show that DNA remains essentially in the B conformer with a tendency in the hydrated model to adopt a slightly distorted, unwound, and stretched conformation in comparison to standard B-DNA. Concerning sugar puckers, the mean pseudorotation phases of the adenine residues are systematically higher than those of the thymine residues, except in the case of the hydrated model for which a articular behavior is observed for the adenine strand. In this case, the terminal bases oscillate between C2'-endo and O4'-endo and the central ones stay in the C3'-endo domain. The mean lifetimes of the internal Watson-Crick H-bond (A) HN6...O4(T) are also dependent on the base pairs included in the calculation, excepted for the implicit solvent simulation at constant temperature. The three-center H bonds have very small mean lifetimes in all three cases of MD simulation. In the minor groove of the hydrated model, a spine of hydration is found as observed by x-ray crystallography and other theoretical simulations. On the basis of the rms deviations, it appears that the fully hydrated simulation has not reached a plateau at the end of the run, while the implicit simulation at constant energy seems to have converged. At constant temperature, very large oscillations in rms deviations are observed.

A molecular dynamics study of conformational changes and hydration of left-handed d(CGCGCGCGCGCG)2 in a nonsalt solution

Twelve dinucleotides (one complete turn) of left-handed, flexible, double-helix poly(dG-dC) Z-DNA have been simulated in aqueous solution with K+ counterions for 70 ps. Most of the d(GpC) phosphates have rotated in accordance with a ZI----ZII transition. The ZII conformation was probably partly stabilized by counterions, which coordinate one of the anionic oxygens and the guanine-N7 of the next (5'----3' direction) base. The presence of base-coordinating ions close to the helical axis rotated and pulled about half of the d(CpG) phosphates further into the groove. These ions also gave rise to rather large deviations from the crystal structure (ZI) with their tendency of pulling the bases closer toward the helical axis. A flipping of the orientation about the glycosyl bond from the +sc to the -sc region was observed for one guanosine, also leading to deviations from the crystal structure. Many bridges containing one or two water molecules were found, with a dominance for the latter. They essentially formed a network of intra- and interstrand bridges between anionic and esterified phosphate oxygens. A "spine" of water molecules could be distinguished as a dark zig-zag pattern in the water density map. The lifetime of a bridge containing one water was about twice as long as that of a two-water bridge and it lasted 5-15 times longer than a hydrogen bond in water. The lifetimes were also calculated for a selection of bridge types, in order of decreasing stability: O1P/O2P ... W ... O'4 much greater than O1P/O2P ... W ... guanine-N2 greater than O1P/O2P ... W ... O1P/O2P. The reorientational motion of water molecules in the first hydration shell around selected groups was slowed down considerably compared to bulk water and the decreasing order of correlation times was guanine-N2 greater than O'4 greater than O'3/O'5 greater than O1P/O2P.

Hydration of DNA bases: analysis of crystallographic data

We present a systematic analysis of water structure around nucleic acid bases. We have examined 28 crystal structures of oligonucleotides, and have studied the patterns of water around the four bases, guanine, cytosine, adenine, and thymine. The geometries of water positions were calculated up to 4.00 A from base atoms. We have found conformation-dependent differences in both the geometry and extent of hydration of the bases.

Stabilities of double- and triple-strand helical nucleic acids

In this selected literature survey, we have seen that the stabilities of duplexes and triplexes are governed by the vertical base stacking, the horizontal specific base-paired H-bonding and the environmental parameters. The entropic contribution in the solvation/desolvation process is important in driving the aggregation of NA strands and duplex formation, but base stacking and specific H-bonding maintain the helical order. Triplex formation shares most of the physical environmental prerequisites with those of duplex NAs. However, some additional environmental conditions are often needed. Only in low pH solution is the polycytidylic strand protonated and, thus, it is possible for the strand to bind to a G.C duplex sequence to give the C+(G.C) triplex. High ionic strength is often necessary for the screening of inter-phosphate repulsion due to the high linear charge density in triplexes. The presence of specific counterions is important for complexation. In the absence of negative supercoiling, existence of an intramolecular triplex is rare except under very acidic conditions for the formation of C+(G.C)-type intramolecular triplex. As expected, the stabilities of both inter- and intramolecular triplexes increase with sequence length. The thermodynamic principles of helix-coil transition of oligo-duplex may be described by the van't Hoff relationship, which assumes a two-state cooperative melting profile. Thus, the enthalpy, entropy and free energy of transition can be evaluated from the experimental melting curves (e.g. OD, DSC). For polynucleotides, because of the non-two-state nature of transition, the simple van't Hoff relationship is no longer valid, and direct calorimetry is needed to obtain reliable thermodynamic parameters. The pH and salt concentration dependence of duplex stability can be formulated and derived from a van't Hoff equation. Base-stacking patterns are simple in duplexes but not so in triplexes due to the diversity in triplet schemes. The sequence dependence of base stacking for duplexes has been characterized and employed to predict the stability of an arbitrary sequence. In conclusion, the stability of duplex is relatively well-characterized by thermodynamic data in terms of both base stacking and specific H-bonding. Thermodynamic studies of triplexes have been far fewer in number. Oligonucleotides have found application in the detection and localization of a mRNA or its gene, the detection of bacterial or viral sequences, and the inhibition of the translation of mRNA and the transcription and replication of DNA (Englisch and Gauss, 1991). In a different approach, oligonucleotides have been targeted directly to a DNA duplex motif of a gene in order to inhibit the expression at the beginning of the transcriptional process.(ABSTRACT TRUNCATED AT 400 WORDS)