Aqueous Hydration of Nucleic Acid Constituents: Monte Carlo Computer Simulation Studies

Abstract

Monte Carlo computer simulations were performed on dilute aqueous solutions of thymine, cytosine, uracil, adenine, guanine, the dimethyl phosphate anion in the gauche-gauche conformation and a ribose and deoxyribose derivative. The aqueous hydration of each molecule was analysed in terms of quasi-component distribution functions based on the Proximity Criterion, and partitioned into hydrophobic, hydrophilic and ionic contributions. Color stereo views of selected hydration complexes are also presented. A preliminary discussion of the transferability of functional group coordination numbers is given. The results enable to comment on two current problems related to the hydration of nucleic acids: a) the theory of Dickerson and coworkers on the role of water in the relative stability of the A and B forms of DNA and b) the idea of water bridges and filaments emerging from the computer simulation results on the hydration of DNA fragments by Clementi.     

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.     

Monte Carlo studies on water in the dCpG/proflavin crystal hydrate

The extensive water network identified in the crystallographic studies of the dCpG/Proflavin hydrate by Neidle, Berman and Shieh (Nature 288, 129, 1980) forms an ideal test case for a) assessing the accuracy of theoretical calculations on nucleic acid--water systems based on statistical thermodynamic computer simulation, and b) the possible use of computer simulation in predicting the water positions in crystal hydrates for use in the further refinement and interpretation of diffraction data. Monte Carlo studies have been carried out on water molecules in the unit cell of dCpG/proflavin, with the nucleic acid complex fixed and the condensed phase environment of the system treated by means of periodic boundary conditions. Intermolecular interactions are described by potential functions representative of quantum mechanical calculations developed by Clementi and coworkers, and widely used in recent studies of the aqueous hydration of various forms of DNA fragments. The results are analyzed in terms of hydrogen bond topology, hydrogen bond distances and energies, mean water positions, and water crystal probability density maps. Detailed comparison of calculated and experimentally observed results are given, and the sensitivity of results to choice of potential is determined by comparison with simulation results based on a set of empirical potentials.     

Monte Carlo Studies on Water in the dCpG/Proflavin Crystal Hydrate

Abstract

The extensive water network identified in the crystallographic studies of the dCpG/Proflavin hydrate by Neidle, Berman and Shieh (Nature 288, 129, 1980) forms an ideal test case for a) assessing the accuracy of theoretical calculations on nucleic acid—water systems based on statistical thermodynamic computer simulation, and b) the possible use of computer simulation in predicting the water positions in crystal hydrates for use in the further refinement and interpretation of diffraction data. Monte Carlo studies have been carried out on water molecules in the unit cell of dCpG/proflavin, with the nucleic acid complex fixed and the condensed phase environment of the system treated by means of periodic boundary conditions. Intermolecular interactions are described by potential functions representative of quantum mechanical calculations developed by Clementi and coworkers, and widely used in recent studies of the aqueous hydration of various forms of DNA fragments. The results are analyzed in terms of hydrogen bond topology, hydrogen bond distances and energies, mean water positions, and water crystal probability density maps. Detailed comparison of calculated and experimentally observed results are given, and the sensitivity of results to choice of potential is determined by comparison with simulation results based on a set of empirical potentials.     

Water and ion binding around RNA and DNA (C,G) oligomers

The dynamics, hydration, and ion-binding features of two duplexes, the A(r(CG)(12)) and the B(d(CG)(12)), in a neutralizing aqueous environment with 0.25 M added KCl have been investigated by molecular dynamics (MD) simulations. The regular repeats of the same C=G base-pair motif have been exploited as a statistical alternative to long MD simulations in order to extend the sampling of the conformational space. The trajectories demonstrate the larger flexibility of DNA compared to RNA helices. This flexibility results in less well defined hydration patterns around the DNA than around the RNA backbone atoms. Yet, 22 hydration sites are clearly characterized around both nucleic acid structures. With additional results from MD simulations, the following hydration scale for C=G pairs can be deduced: A-DNA相似文献   

Hydration properties of DNA-lysine gels by microwave dielectric measurements as a function of temperature

Dielectric measurements by a cavity perturbation method at 10 GHz in the temperature range from-20°C to +45°C are reported for aqueous gels of herring sperm DNA in the presence of 1 or 3 lysine molecules per nucleotide. Measurements for lysine-water and DNA-water systems are also reported. The experimental results can be accounted for by the presence of interfacial water, with dielectric properties different from those of bulk water, and are analyzed in terms of a three component equation (solute molecules, interfacial water and bulk water) to calculate hydration parameters of the systems. The lysine molecule is found to coordinate a particular number of water molecules, in agreement with the literature. The specific hydration of DNA is reduced by the presence of lysine, indicating a direct interaction between the polyion and the aminoacid: a decrease to about 50% was observed at a ratio of one molecule of lysine per nucleotide. A suggestion is made that the interaction is mainly electrostatic in nature.     

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.     

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     

Hydration effects accompanying the formation of DNA complexes with some ligands

In the range of millimeter wavelengths the dielectric properties of aqueous solutions of some biologically active ligands (potential anticarcinogen chlorophyllin; pharmacological drug caffeine; polyamine putrescine; mutagens proflavine and ethidium bromide; actinocin derivative, an analogue of antitumor antibiotic actinomycin D) and DNA complexes with these substances were studied. It was shown that complex formation is accompanied by the change in dielectric properties of the solution. These changes during interaction of DNA with the first three compounds correspond to a decrease in hydration (compared with the total hydration of free components), and in other cases they cause an increase in hydration. The number of water molecules bound with both the ligand and DNA nucleotide in the complex was estimated. The results were compared with existing models of DNA interaction with the studied substances.     

Hydration of Cytidine, 2′-Deoxycytidine and their Phosphate Salts in the Aqueous Solutions. A Molecular Dynamics Computer Simulation Study

Abstract

To compare the hydration pattern of the cytidine (Cyd) and 2′-deoxycytidine (dCyd) in the aqueous solutions at the level of microscopic interactions, Molecular Dynamics (MD) computer simulations have been undertaken. The results indicate that the hydration of the heterocyclic base moiety in cytidine and 2′-deoxycytidine has a hydrophobic character. None of the three potential Watson—;Crick base pair centres hydrogen bonds with the water molecules and the formation of something akin to a clathrate cage structure of water around base moieties of nucleosides in the aqueous solution is suggested. In contrast, the hydration of Cyd and dCyd sugar moieties shows a hydrophilic character and the three-dimensional networks of H-bonds involving all hydrophilic centres are formed differently around the ribose and 2′-deoxyribose. The sugar hydroxyl groups participate in the hydrogen bonding with water both as H-donor and as H-acceptor. Their donor-acceptor abilities have been evaluated and compared. The coordination numbers, the geometrical data of the first hydration shell, and the number of hydrogen bonds have been calculated. The changes in the pattern of hydration with the increased concentration of nucleosides and upon nucleoside protonation are discussed. The analysis of the pairwise interaction energies are also presented.     

The effect of a low-frequency electromagnetic field on DNA molecules in aqueous solutions

A chemiluminescence study showed that hepatitis B virus (HBV) and hepatitis C virus (HCV) DNA amplicons are capable of induced radiation when exposed to electromagnetic fields (EMFs) that range from 7.5 to 30 Hz in frequency and from 24 to 40 A/m in field strength. An EMF with a frequency of 9 Hz was shown to exert the greatest effect on aqueous solutions of the hepatitis virus DNA amplicons. The hydration shell of the DNA amplicons was observed to change. The change in the DNA hydration shell on exposure to a low-frequency EMF was presumed to restore hydrogen bonds, to induce crosslinks, and to facilitate DNA repair.     

Hydration of nucleic acid fragments: comparison of theory and experiment for high-resolution crystal structures of RNA, DNA, and DNA-drug complexes.

A computationally efficient method to describe the organization of water around solvated biomolecules is presented. It is based on a statistical mechanical expression for the water-density distribution in terms of particle correlation functions. The method is applied to analyze the hydration of small nucleic acid molecules in the crystal environment, for which high-resolution x-ray crystal structures have been reported. Results for RNA [r(ApU).r(ApU)] and DNA [d(CpG).d(CpG) in Z form and with parallel strand orientation] and for DNA-drug complexes [d(CpG).d(CpG) with the drug proflavine intercalated] are described. A detailed comparison of theoretical and experimental data shows positional agreement for the experimentally observed water sites. The presented method can be used for refinement of the water structure in x-ray crystallography, hydration analysis of nuclear magnetic resonance structures, and theoretical modeling of biological macromolecules such as molecular docking studies. The speed of the computations allows hydration analyses of molecules of almost arbitrary size (tRNA, protein-nucleic acid complexes, etc.) in the crystal environment and in aqueous solution.     

Hydration of enzyme in nonaqueous media is consistent with solvent dependence of its activity

Water plays an important role in enzyme structure and function in aqueous media. That role becomes even more important when one focuses on enzymes in low water media. Here we present results from molecular dynamics simulations of surfactant-solubilized subtilisin BPN' in three organic solvents (octane, tetrahydrofuran, and acetonitrile) and in pure water. Trajectories from simulations are analyzed with a focus on enzyme structure, flexibility, and the details of enzyme hydration. The overall enzyme and backbone structures, as well as individual residue flexibility, do not show significant differences between water and the three organic solvents over a timescale of several nanoseconds currently accessible to large-scale molecular dynamics simulations. The key factor that distinguishes molecular-level details in different media is the partitioning of hydration water between the enzyme and the bulk solvent. The enzyme surface and the active site region are well hydrated in aqueous medium, whereas with increasing polarity of the organic solvent (octane --> tetrahydrofuran --> acetonitrile) the hydration water is stripped from the enzyme surface. Water stripping is accompanied by the penetration of tetrahydrofuran and acetonitrile molecules into crevices on the enzyme surface and especially into the active site. More polar organic solvents (tetrahydrofuran and acetonitrile) replace mobile and weakly bound water molecules in the active site and leave primarily the tightly bound water in that region. In contrast, the lack of water stripping in octane allows efficient hydration of the active site uniformly by mobile and weakly bound water and some structural water similar to that in aqueous solution. These differences in the active site hydration are consistent with the inverse dependence of enzymatic activity on organic solvent polarity and indicate that the behavior of hydration water on the enzyme surface and in the active site is an important determinant of biological function especially in low water media.     

Structural transitions of deoxyribonucleic acid in aqueous electrolyte solutions. II. The role of hydration.

The data and approach reported in paper I (Hanlon et al., 1975, preceding paper) have been used to calculate the fractional changes in secondary structure of calf thymus deoxyribonucleic acid which occur in aqueous solutions as a function of the concentration of NaCl, KCl, LiCl, CsCl, and NH4Cl. There is a continuous loss in the "B" character of the nucleic acid with concomitant production of the C and, in some instances, an A form, as well, as the salt concentration increases. Sedimentation velocity studies suggest that there is an accompanying change in the hydrodynamic characteristics of the DNA molecules, as well. Utilizing the existing hydration data in the literature (Hearst and Vinograd, 1961a,b; Hearst, 1965; Tunis and Hearst, 1968a; Cohen and Eisenberg, 1968; Falk et al., 1962, 1963a,b), we have found that a gradual loss of "B" character and a decrease in the frictional coefficient of DNA occur as the net hydration of DNA is reduced from the fully hydrated from (60-80 mol of H2O/mol of nucleotide) to values of ca. 12-14 mol of H2O/mol of nucleotide. Below that value, a more precipitous decrease in these properties occurs. Extrapolation of the linear relationship observed between the fractional B content and the net hydration in the latter regions yield values of ca. 18 mol of H2O/mol of nucleotide at 100% B and ca. 4 mol of H2O/mol of nucleotide at 0% B (i.e., 100% C or C + A) for the alkali metal salts of DNA. The ammonium salt retains somewhat more H2O in the C and A forms (ca. 7). These results together with the hydration site assignments of Falk et al. (1962, 1963a,b) are interpreted in terms of a hydration model for DNA in aqueous solution in which an intact primary hydration shell of ca. 18 mol of H2O/mol of nucleotide is required for the maintenance of the "B" conformation. Removal of all but those water molecules solvating the phosphate groups results in the conversion to the C forms, predominantly, with a small amount of A structure formed as well in some salts. The accompanying changes in the sedimentation coefficients suggest that the DNA molecule assumes a more compact and/or flexible form under these conditions in which it is mainly in the C and A structures. The combination of these two events which ensue upon dehydration create a polymeric structure which can be more easily packaged in biological systems.     

Dynamic properties of water bound to matrices of natural DNA and model complexes

From experimental data on the hydration energetics of nucleic acids obtained by differential scanning calorimetry under isothermal conditions, dielectric relaxation time tau d and "free volume" Vf occupied by water molecules in hydration shells of natural DNA and model polyribonucleotides were calculated. In addition, systems consisting of dinucleotides ApA, TpT, UpU, TpU, UpT and water clusters of various sizes (from 20 to 400 water molecules) were studied by Monte Carlo computer simulation. It was shown that, as water content in systems increases, the dynamic characteristics of bound water obtained with both methods approached the values for bulk water.     

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.     

Spin isomeric selectivity of water molecules upon DNA hydration

Four-photon coherent scattering of laser radiation was used to study the influence of DNA on the content of quasi-free ortho and para isomers of water molecules in its aqueous solution. It was shown that the concentration of quasi-free molecules that form the rotational spectrum of spin isomers increases considerably in the hydration shell of the DNA molecule as compared with pure water. The increase in the concentration of spin isomers occurs disproportionally. In the presence of DNA, the intensity of the rotational spectrum of ortho isomers is on the average much greater than that of para isomers. It was also demonstrated that the character of hydration and the ortho/para ratio change noticeably upon DNA denaturation, which may be evidence of changes in preferable solvation of DNA during its denaturation. The data obtained allowed us to assume that the stability of different biologically important states of macromolecules can be changed by varying the relative concentration of water spin isomers in solution.     

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.     

Monte Carlo computer simulation of the aqueous hydration of the glycine zwitterion at 25 degree C

Monte Carlo computer simulation on a dilute aqueous solution of the glycine zwitterion are reported. The results are presented in terms of the Quasi-Component Distribution Functions (QCDF) of Ben Naim and partitioned into atomic and functional group contributions using the Proximity Criterion. The Proximity Criterion analysis has been extended to orientational properties and a new normalization procedure has been introduced for the radial distribution functions obtained by the Proximity Criterion. The solvation environment of the glycine zwitterion is found to contain, on the average, 14.4 water molecules out of which 3.2 belong to the ammonium group, 6.1 to the methylene group and 5.1 to the carboxyl group. The importance of the many-body statistical mechanical approach to hydration is emphasized by our finding that the configuration corresponding to the absolute minimum of the glycine zwitterion-water potential surface was found to have negligible statistical weight in the aqueous simulation.