Molecular dynamics studies of the behaviour of water molecules and small ions in concentrated solutions of polymeric B-DNA

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.

Ion and solvent density distributions around canonical B-DNA from integral equations

We calculate the water and ion spatial distributions around charged oligonucleotides using a renormalized three-dimensional reference interaction site theory coupled with the HNC closure. Our goal is to understand the balance between inter-DNA strand forces and solvation forces as a function of oligonucleotide length in the short strand limit. The DNA is considered in aqueous electrolyte solutions of 1 M KCl, 0.1 M KCl, or 0.1 M NaCl. The current theoretical results are compared to molecular dynamics (MD) simulations and experiments. It is found that the integral equation (IE) theory replicates the MD and the experimental results for the base-specific hydration patterns in both the major and the minor grooves. We are also able to discern characteristic structural pattern differences between Na(+) and K(+) ions. When compared to Poisson-Boltzmann methods, the IE theory, like simulation, predicts a richly structured ion environment, which is better described as multilayer rather than double layer.

Dielectric properties of fully hydrated nucleotides in the terahertz frequency range

We use terahertz time domain spectroscopy (THz-TDS) to determine the complex frequency-dependent dielectric response of all four nucleotides at different dilute concentrations. In addition, the suspension model's ability to extract the dielectric response of just the nucleotide with the hydration shell epsilon(b) excluding the dielectric information relating to the bulk will be verified. The suspension model enables us to make the determination that the nucleotides have influences on the water molecules out to the fourth hydration shell. We use a two Debye relaxation fit model for water, all concentrations and all epsilon(b) values. We observed how the nucleotides affect the relaxation parameters in relation to that of pure bulk water. With this information, we notice a transition between purines and pyrimidines, where one is a hydrogen-bond network structure building type material with a low concentration increment and the other is a structure breaking type material with a low concentration decrement. Due to conductivity measurements, we determine that kinetic depolarization is a negligible affect compared to that of dielectric saturation, which we find to dominate where a decrement is found.

Monte Carlo simulations of flexible polyelectrolytes inside viral capsids with dodecahedral charge distribution

Structural properties of encapsidated flexible polyelectrolytes in viral capsids with dodecahedral charge distribution have been investigated by Monte Carlo simulations using a coarse-grained model. Several capsid charge distributions ranging from a homogeneous surface charge distribution (lambda=0) to a complete dodecahedral distribution (lambda=1) at constant total capsid charge and fixed radial location of the capsid charges have been considered. The radial and lateral organizations of the polyelectrolyte have been examined as a function of the polyelectrolyte length and capsid charge distribution. With short polyelectrolytes a single polyelectrolyte layer was formed at the inner capsid surface, whereas at increasing polyelectrolyte length also a uniform polyelectrolyte density inside the surface layer was established. At low lambda , the polyelectrolyte layer was laterally isotropic, but at lambda> or =0.05 a dodecahedral structure started to appear. At lambda=1 , the polyelectrolyte followed essentially a path along the edges of a dodecahedron. With sufficiently long chains, all edges were decorated with polyelectrolyte, facilitated by loop formation. For an undercharged capsid, the capsid counterions inside the capsid also adopted a dodecahedral distribution.

The effect of salt concentration on DNA conformation transition: a molecular-dynamics study

We performed three 3-ns molecular dynamics simulations of d(CGCGAATTCGCG)2 using the AMBER 8 package to determine the effect of salt concentration on DNA conformational transitions. All the simulations were started with A-DNA, with different salt concentrations, and converged with B-DNA with similar conformational parameters. However, the dynamic processes of the three MD simulations were very different. We found that the conformation transition was slow in the solution with higher salt concentration. To determine the cause of this retardation, we performed three additional 1.5-ns simulations starting with B-DNA and with the salt concentrations corresponding to the simulations mentioned above. However, astonishingly, there was no delayed conformation evolution found in any of the three simulations. Thus, our simulation conclusion is that higher salt concentrations slows the A --> B conformation transition, but have no effect on the final stable structure. [Figure: see text].

Brownian dynamics simulations of ion atmospheres around polyalanine and B-DNA: effects of biomolecular dielectric

We have extended an earlier Brownian dynamics simulation algorithm for simulating the structural dynamics of ions around biomolecules to accommodate dielectric inhomogeneity. The electrostatic environment of a biomolecule immersed in water was obtained by numerically solving the Poisson equation with the biomolecule treated as a low dielectric region and the solvent treated as a high dielectric region. Instead of using the mean-field type approximations of ion interactions as in the Poisson-Boltzmann model, the ions were treated explicitly by allowing them to evolve dynamically under the electrostatic field of the biomolecule. This model thus accounts for ion-ion correlations and the finite-size effects of the ions. For a 13-residue alpha-helical polyalanine and a 12-base-pair bp B-form DNA, we found that the choice of the dielectric constant of the biomolecule has much larger effects on the mean ionic structure around the biomolecule than on the fluctuational and dynamical properties of the ions surrounding the biomolecule.

Molecular dynamics simulation of nucleic acids

We review molecular dynamics simulations of nucleic acids, including those completed from 1995 to 2000, with a focus on the applications and results rather than the methods. After the introduction, which discusses recent advances in the simulation of nucleic acids in solution, we describe force fields for nucleic acids and then provide a detailed summary of the published literature. We emphasize simulations of small nucleic acids ( approximately 6 to 24 mer) in explicit solvent with counterions, using reliable force fields and modern simulation protocols that properly represent the long-range electrostatic interactions. We also provide some limited discussion of simulation in the absence of explicit solvent. Absent from this discussion are results from simulations of protein-nucleic acid complexes and modified DNA analogs. Highlights from the molecular dynamics simulation are the spontaneous observation of A B transitions in duplex DNA in response to the environment, specific ion binding and hydration, and reliable representation of protein-nucleic acid interactions. We close by examining major issues and the future promise for these methods.

Sodium and chlorine ions as part of the DNA solvation shell

The distribution of sodium and chlorine ions around DNA is presented from two molecular dynamics simulations of the DNA fragment d(C(5)T(5)). (A(5)G(5)) in explicit solvent with 0.8 M additional NaCl salt. One simulation was carried out for 10 ns with the CHARMM force field that keeps the DNA structure close to A-DNA, the other for 12 ns with the AMBER force field that preferentially stabilizes B-DNA conformations (, Biophys. J. 75:134-149). From radial distributions of sodium and chlorine ions a primary ion shell is defined. The ion counts and residence times of ions within this shell are compared between conformations and with experiment. Ordered sodium ion sites were found in minor and major grooves around both A and B-DNA conformations. Changes in the surrounding hydration structure are analyzed and implications for the stabilization of A-DNA and B-DNA conformations are discussed.

Modeling high-resolution hydration patterns in correlation with DNA sequence and conformation

Hydration around the DNA fragment d(C5T5).(A5G5) is presented from two molecular dynamics simulations of 10 and 12 ns total simulation time. The DNA has been simulated as a flexible molecule with both the CHARMM and AMBER force fields in explicit solvent including counterions and 0.8 M additional NaCl salt. From the previous analysis of the DNA structure B-DNA conformations were found with the AMBER force-field and A-DNA conformations with CHARMM parameters. High-resolution hydration patterns are compared between the two conformations and between C.G and T.A base-pairs from the homopolymeric parts of the simulated sequence. Crystallographic results from a statistical analysis of hydration sites around DNA crystal structures compare very well with the simulation results. Differences between the crystal sites and our data are explained by variations in conformation, sequence, and limitations in the resolution of water sites by crystal diffraction. Hydration layers are defined from radial distribution functions and compared with experimental results. Excellent agreement is found when the measured experimental quantities are compared with the equivalent distribution of water molecules in the first hydration shell. The number of water molecules bound to DNA was found smaller around T.A base-pairs and around A-DNA as compared to B-DNA. This is partially offset by a larger number of water molecules in hydrophobic contact with DNA around T.A base-pairs and around A-DNA. The numbers of water molecules in minor and major grooves have been correlated with helical roll, twist, and inclination angles. The data more fully explain the observed B-->A transition at low humidity.

RNA hydration: three nanoseconds of multiple molecular dynamics simulations of the solvated tRNA(Asp) anticodon hairpin

The hydration of the tRNA(Asp) anticodon hairpin was investigated through the analysis of six 500 ps multiple molecular dynamics (MMD) trajectories generated by using the particle mesh Ewald method for the treatment of the long-range electrostatic interactions. Although similar in their dynamical characteristics, these six trajectories display different local hydration patterns reflecting the landscape of the "theoretical" conformational space being explored. The statistical view gained through the MMD strategy allowed us to characterize the hydration patterns around important RNA structural motifs such as a G-U base-pair, the anticodon U-turn, and two modified bases: pseudouridine and 1-methylguanine. The binding of ammonium counterions to the hairpin has also been investigated. No long-lived hydrogen bond between water and a 2'-hydroxyl has been observed. Water molecules with long-residence times are found bridging adjacent pro-Rp phosphate atoms. The conformation of the pseudouridine is stiffened by a water-mediated base-backbone interaction and the 1-methylguanine is additionally stabilized by long-lived hydration patterns. Such long-lived hydration patterns are essential to ensure the structural integrity of this hairpin motif. Consequently, our simulations confirm the conclusion reached from an analysis of X-ray crystal structures according to which water molecules form an integral part of nucleic acid structure. The fact that the same conclusion is reached from a static and a dynamic point of view suggests that RNA and water together constitute the biologically relevant functional entity.

Theoretical studies of the intercalation of 9-hydroxyellipticine in DNA

Extensive molecular dynamics (MD) simulations have been used to investigate the intercalative binding of 9-hydroxyellipticine to the DNA oligonucleotide d(ATATATATATAT)2. Four independent simulations differing in the initial orientation of the drug at the intercalation site were carried out, and compared both with each other and a control simulation of the free DNA sequence. The structure of the latter was compared with structures obtained from x-ray crystallography and nmr spectroscopy, as well as the theoretically derived "alternating B-DNA" model [A. Klug et al. (1979), Journal of Molecular Biology, Vol. 131, p. 669]. The alternation of twist angles observed in experimental structures was reproduced in the simulation. All four independent simulations of the drug-DNA intercalation complex converged in placing the pyridine ring of the ellipticine chromophore in the major groove; in one case this involved a 180 degrees rotation of the drug at the intercalation site. At a more detailed level, the drug is seen to be capable of adopting several distinct orientations, each stable over a period of hundreds of pico-seconds. Despite the presence of several polar groups in the drug, however, no direct hydrogen bonding to the DNA occurs; instead, interactions between the methyl groups of the drug and the thymine bases at the intercalation site appear important in determining the orientational preferences of the drug. Comparison of the intercalation complexes with the free DNA sequence shows a degree of unwinding resulting from intercalation, in good agreement with experimental results, but spread over the three central base-pair steps, not confined to the intercalation site itself. Measurements of torsional rigidity indicate only a slight stiffening of the DNA restricted to the immediate site of intercalation. The structures obtained from the MD simulations were used to calculate theoretical CD spectra, with separate simulations giving very different results. This appears to indicate that given an accurate assignment of the main electronic transition dipole moment of the ellipticine chromophore, discrimination of the more realistic binding geometries may be possible. The relative merits of the various drug orientations observed in the simulations are discussed and a perpendicular orientation of the drug at the intercalation site is considered to be the most consistent with experimental data. While the simulations themselves represent a total of over 2 ns, however, the differences apparent between independent runs indicate that longer simulation times will be required before a complete, unequivocal view of DNA intercalation is obtained.

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.

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)