Regularities in formation of the spine of hydration in the DNA minor groove and its influence on the DNA structure

Computer calculations as well as an analysis of space-filling models and literature data allowed the following conclusions to be made: an ordered spine of water in the DNA minor groove, similar to that revealed in the CGCGAATTCGCG crystal, seems to exist in DNA crystals, fibers and solutions; it is shown that this spine may be formed on A/T runs containing no TA step while on the TA step the spine is disrupted; the existence of this spine changes the double helix structure stabilizing a definite DNA conformation; the spine of hydration makes the DNA more stable to conformational transitions. These conclusions permit us to interpret a large body of experimental data on DNA crystals, fibers and solutions. The role of water bridges constituting the first hydration shell of the ordered spine of water is discussed in connection with the B-to-A transition.     

A Monte Carlo simulation study of the aqueous hydration of r(GpC)2: comparison with crystallographic ordered water sites

Monte Carlo computer simulation is described for the dinucleotide duplex rGpC together with 562 water molecules at an environmental density of 1 g/cc in a cubic cell under periodic boundary conditions. Water-water interactions were treated using the TIP4P potential and the solute water interactions by TIP4P spliced with the nonbonded interactions from the AMBER 3.0 force field. The simulation was subjected to proximity analysis to obtain solute coordinate numbers and pair interaction energies for each solute atom. Hydration density distributions partitioned into contributions from the major groove side, the minor groove side and the sugar-phosphate backbone were examined, and the probabilities of occurrence for one- and two-water bridges in the simulation were enumerated. The results were compared with observations of crystallographic ordered water sites from x-ray diffraction studies on G and C containing small molecules, and in crystal structure determinations of the sodium, calcium, and ammonium salts of rGpC. The calculated results are generally consistent with the observed sites, except for cytosine N4, where a hydration site is predicted yet none observed in rGpC salts, and for guanine N3, which appears in this calculation to compete unfavorably with the adjacent donor site at guanine N2. There is, however, a significant probability of finding a one-water G-N3-W-G-N2 bridge indicated in the simulation. An explanation for the guanine N3 discrepancy in terms of electrostatic potentials is also offered. The calculated one- and two-water bridges in the rGpC hydration complex coincide in a number of cases to those observed in the ordered water structure of the sodium rGpC crystal hydrate.     

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.     

The B- to A-DNA transition and the reorganization of solvent at the DNA surface

DNA geometry depends on relative humidity. Using the CHARMM22 force field to push B-DNA to A-DNA, a molecular dynamics simulation of a mixed-sequence 24-basepair DNA double-stranded oligomer, starting from B-DNA, was carried out to explore both the mechanism of the transition and the evolution of hydration patterns on the surface of DNA. Over the 11-ns trajectory, the transition recapitulates the slide-first, roll-later mechanism, is opposed by DNA electrostatics, and is favored by an increasing amount of condensed sodium ions. Hydration was characterized by counting the hydrogen bonds between water and DNA, and by the number of water bridges linking two DNA atoms. The number of hydrogen bonds between water and DNA remains constant during the transition, but there is a 40% increase in the number of water bridges, in agreement with the principle of economy of hydration. Water bridges emerge as delicate sensors of both structure and dynamics of DNA. Both local flexibility and the frustration of the water network on the surface of DNA probably account for the low populations and short residence times of the bridges, and for the lubricant role of water in ligand-DNA interactions.     

Transcription of the his3 gene region in Saccharomyces cerevisiae

The dodecamer d(CpGpCpGpApApTpTpCpGpCpG) or C-G-C-G-A-A-T-T-C-G-C-G crystallizes as slightly more than one full turn of right-handed B-DNA. It is surrounded in the crystal by one bound spermine molecule and 72 ordered water molecules, most of which associate with polar N and O atoms at the exposed edges of base-pairs. Hydration within the major groove is principally confined to a monolayer of water molecules associated with exposed N and O groups on the bases, with most association being monodentate. Waters hydrating backbone phosphate oxygens tend not to be ordered, except where they are immobilized by 5-methyl groups from nearby thymines. In contrast, the minor groove is hydrated in an extensive and regular manner, with a zigzag “spine” of first- and second-shell hydration along the floor of the groove serving as a foundation for less-regular outer shells extending beyond the radius of the phosphate backbone. This spine network bridges purine N-3 and pyrimidine O-2 atoms in adjacent base-pairs. It is particularly regular in the A-A-T-T center, and is disrupted at the C-G-C-G ends, in part by the presence of the N-2 amino groups on guanine residues. The minor groove hydration spine may be responsible for the stability of the B form of polymers containing only A · T and I · C base-pairs, and its disruption may explain the ease of transition to the A form of polymers with G · C pairs.     

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.     

Hydration and structure of hemiprotonated polycytidylic acid in films

Structural transitions of poly(rC)-Ka+ in humid films with different water content were studied by infrared spectroscopy and piezogravimetry. From analysis of the hydration isotherms and the dependence of spectral parameters (frequencies and intensities of the main bands) on n the hydration sites of the polynucleotide were determined (C2O, O4', N4H2, N1, PO2-, C2'OH). It was found that the transition of the polynucleotide from the unordered state to a double-stranded complex poly(rC+).poly(rC) occurs in the interval of n from 2 to 8. The value n = 8 corresponds to the total hydration of poly(rC). A model of hydration of poly(rC+).poly(rC) based on the experimental results and known X-ray parameters of this double helix complex is proposed. The most important feature of the model is the presence of single water bridges between PO2(-)-groups in the first hydration shell of each chain and triple water bridges between O4', N4H2 and C2'OH- atomic groups of opposite chains. The experimental results obtained and the proposed structure of hydration environment of poly(rC+).poly(rC) suggest that the stabilization of this complex is stabilized by the intra- and inter-chain water bridges and hydrogen bonds between pairs of cytosine bases.     

The molecular stoichiometric hydration model (SHM) as applied to tendon/collagen, globular proteins and cells

This report describes and documents the presence of multiple water-of-hydration fractions on proteins and in cells. Initial studies of hydration fractions in g of water/g of DM (dry mass) for tendon/collagen led to the development of the molecular SHM (stoichiometric hydration model) and the development of methods for calculating the size of hydration fractions on a number of different proteins of known amino acid composition. The water fractions have differences in molecular motion and other physical properties due to electrostatic interactions of polar water molecules with electric fields generated by covalently bound pairs of opposite partial charge on the protein backbone. The methods allow calculation of the size of four hydration fractions: single water bridges, double water bridges, dielectric water clusters over polar-hydrophilic surfaces and water clusters over hydrophobic surfaces. These four fractions provide monolayer water coverage. The predicted SHM hydration fractions match closely measured hydration fraction values for collagen and for globular proteins. This report also presents water sorption findings that support the SHM. The SHM is applicable for cell systems where it has been studied. In seven cell systems studied, more than half of all of the cell water had properties unlike those of bulk water. The SHM predicts and explains the commonly cited and measured bound water fraction of 0.2-0.4 g of water/g of DM on proteins. The commonly accepted concept that water beyond this bound water fraction can be considered bulk-like water in its physical properties is unwarranted.     

Monte Carlo simulation of hydration of the nucleic acid fragments

Systems containing a base or a base pair and 25 water molecules, as well as a helical stack and 30 water molecules per base pair, have been simulated. Changes in the base hydration shell structure, after the bases have been included into the pair and then into the base pair stack, are discussed. Hydration shells of several configurations of the base pair stacks are discussed. Probabilities of formation of the hydrogen-bonded bridges of 1, 2 and 3 water molecules between hydrophilic centres have been estimated. The hydration shell structure was shown to depend on the nature of the base pair and on the stack configuration, while dependence of the global hydration shell characteristics on the stack configuration has been proved to be rather slight. The most typical structural elements of hydration shells, in the glycosidic (minor in B-like conformation) and non-glycosidic (major) grooves, for different configurations of AU and GC stacks, have been found and discussed. The number of hydrogen bonds between water molecules and bases per water molecule was shown to change upon transformation of the stack from A to B configuration. This result is discussed in connection with the reasons for B to A conformational transition and the concept of "water economy". Hydration shell patterns of NH2-groups of AU and GC helical stacks differ significantly.     

Hydration dynamics of the collagen triple helix by NMR

The hydration of the collagen-like Ac-(Gly-Pro-Hyp)(6)-NH(2) triple-helical peptide in solution was investigated using an integrated set of high-resolution NMR hydration experiments, including different recently developed exchange-network editing methods. This approach was designed to explore the hydration dynamics in the proximity of labile groups, such as the hydroxyproline hydroxyl group, and revealed that the first shell of hydration in collagen-like triple helices is kinetically labile with upper limits for water molecule residence times in the nanosecond to sub-nanosecond range. This result is consistent with a "hopping" hydration model in which solvent molecules are exchanged in and out of solvation sites at a rate that is not directly correlated to the degree of site localization. The hopping model thus reconciles the dynamic view of hydration revealed by NMR with the previously suggested partially ordered semi-clathrate-like cylinder of hydration. In addition, the nanosecond to sub-nanosecond upper limits for water molecule residence times imply that hydration-dehydration events are not likely to be the rate-limiting step for triple helix self-recognition, complementing previous investigations on water dynamics in collagen fibers. This study has also revealed labile proton features expected to facilitate the characterization of the structure and folding of triple helices in collagen peptides.     

Structure of the hydration envelope of the B-form of polydeoxyribonucleotides poly(dA-dC).poly(dG-dT) and poly(dA-dG).poly(dTs-dT) from the data of Monte-Carlo simulation

The results of a Monte Carlo simulation of the hydration shell of two polynucleotides poly (dA-dC).poly(dG-dT) and poly(dA-dG).poly(dC-dT) are reported. This study is a part of a series of Monte Carlo computations of the hydration of regular polydeoxyribonucleotides with dinucleotide repeat aimed at looking for dependences of hydration shell structure on base sequence. The coordinates of the main local maximal of water density near the polymers and the topology of the most probable one- and two-membered water bridges are published. For most of the sequences a common primary hydration of base edges of successive base pairs is characteristic. The AT-homopolymeric sequence represents an exception with autonomous primary hydration of a base pair in both grooves, which correlates with the sequence-dependent flexibility and the occurrence of bends of DNA.     

Hydrogen and hydration of DNA and RNA oligonucleotides

In this paper, hydrogen bonding interaction and hydration in crystal structures of both DNA and RNA oligonucleotides are discussed. Their roles in the formation and stabilization of oligonucleotides have been covered. Details of the Watson-Crick base pairs G.C and A.U in DNA and RNA are illustrated. The geometry of the wobble (mismatched) G.U base pairs and the cis and almost trans conformations of the mismatched U.U base pairs in RNA is described. The difference in hydration of the Watson-Crick base pairs G.C, A.U and the wobble G.U in different sequences of codon-anticodon interaction in double helical molecules are indicative of the effect of hydration. The hydration patterns of the phosphate, the 2'-hydroxyl groups, the water bridges linking the phosphate group, N7 (purine) and N4 of Cs or O4 of Us in the major groove, the water bridges between the 2'-hydroxyl group and N3 (purine) and O2 (pyrimidine) in the minor groove are discussed.     

Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.     

Monte Carlo Simulation of Hydration of the Nucleic Acid Fragments

Abstract

Systems containing a base or a base pair and 25 water molecules, as well as a helical stack and 30 water molecules per base pair, have been simulated. Changes in the base hydration shell structure, after the bases have been included into the pair and then into the base pair stack are discussed. Hydration shells of several configurations of the base pair stacks are discussed. Probabilities of formation of the hydrogen-bonded bridges of 1, 2 and 3 water molecules between hydrophilic centres have been estimated. The hydration shell structure was shown to depend on the nature of the base pair and on the stack configuration, while dependence of the global hydration shell characteristics on the stack configuration has been proved to be rather slight. The most typical structural elements of hydration shells, in the glycosidic (minor in B-like conformation) and non-glycosidic (major) grooves, for different configurations of AU and GC stacks, have been found and discussed. The number of hydrogen bonds between water molecules and bases per water molecule was shown to change upon transformation of the stack from A to B configuration. This result is discussed in connection with the reasons for B to A conformational transition and the concept of “water economy”. Hydration shell patterns of NH₂-groups of AU and GC helical stacks differ significantly.     

Aqueous hydration of nucleic acid constituents: Monte Carlo computer simulation studies

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 form 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.     

Phase behavior and structure of aqueous dispersions of sphingomyelin

The phase behavior of bovine brain sphingomyelin in water has been determined by polarizing light microscopy, differential scanning calorimetry, and X-ray diffraction. Lamellar phases, in which water is intercalated between sheets of lipid molecules arranged in the classical bilayer fashion, are present over much of the phase diagram. An order-disorder transition separates the high temperature, liquid crystalline, lamellar phase from a more ordered lamellar phase at low temperatures. The hydration characteristics of sphingomyelin are similar to the structurally related lecithin in that only limited amounts of water are incorporated above and below the transition. Above the transition at 47 degrees C, a maximum of 35% by weight of water can be incorporated between the lipid bilayers, the total thickness at maximum hydration being 60.2 A, the lipid thickness 38 A, and the surface area per lipid molecule at the interface 60 A(2). Water in excess of 35% by weight is present as a separate phase. Below the phase transition, at 25 degrees C a maximum of 42% by weight of water may be incorporated between the lipid bilayers. On increasing the hydration, the lamellar repeat distance increases from 63.5 A to a limiting value of 76 A. Within this hydration range the calculated lipid thickness decreases from 63.5 to 42.5 A, and the surface area per lipid molecule increases from 36.1 to 53.6 A(2). Although these changes may be accounted for by a structure in which the hexagonally packed ordered hydrocarbon chains tilt progressively with respect to the normal to the bilayer plane on increasing hydration, it is possible that changes in other more complex lamellar structures may be responsible for these variations in lipid thickness and surface area.     

Molecular dynamics simulations of trehalose as a 'dynamic reducer' for solvent water molecules in the hydration shell

Systematic computational work for a series of 13 disaccharides was performed to provide an atomic-level insight of unique biochemical role of the alpha,alpha-(1-->1)-linked glucopyranoside dimer over the other glycosidically linked sugars. Superior osmotic and cryoprotective abilities of trehalose were explained on the basis of conformational and hydration characteristics of the trehalose molecule. Analyses of the hydration number and radial distribution function of solvent water molecules showed that there was very little hydration adjacent to the glycosidic oxygen of trehalose and that the dynamic conformation of trehalose was less flexible than any of the other sugars due to this anisotropic hydration. The remarkable conformational rigidity that allowed trehalose to act as a sugar template was required for stable interactions with hydrogen-bonded water molecules. Trehalose made an average of 2.8 long-lived hydrogen bonds per each MD step, which was much larger than the average of 2.1 for the other sugars. The stable hydrogen-bond network is derived from the formation of long-lived water bridges at the expense of decreasing the dynamics of the water molecules. Evidence for this dynamic reduction of water by trehalose was also established based on each of the lowest translational diffusion coefficients and the lowest intermolecular coulombic energy of the water molecules around trehalose. Overall results indicate that trehalose functions as a 'dynamic reducer' for solvent water molecules based on its anisotropic hydration and conformational rigidity, suggesting that macroscopic solvent properties could be modulated by changes in the type of glycosidic linkages in sugar molecules.     

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.     

Orientational order and dynamics of hydration water in a single crystal of bovine pancreatic trypsin inhibitor.

The orientational order and dynamics of the water molecules in form II crystals of bovine pancreatic trypsin inhibitor (BPTI) are studied by (2)H NMR in the temperature range 6-50 degrees C. From the orientation dependence of the single crystal quadrupole splitting and linewidth, the principal components of the motionally averaged quadrupole interaction tensor and the irreducible linewidth components for the orthorhombic crystal are determined. With the aid of water orientations derived from neutron and x-ray diffraction, it is shown that the NMR data can be accounted for by a small number of highly ordered crystal waters, some of which have residence times in the microsecond range. Most of these specific hydration sites must be located at intermolecular contacts. The surface hydration layer that is also present in dilute solution is likely to be only weakly ordered and would then not contribute significantly to the splitting and linewidth from the protein crystal. To probe water dynamics on shorter time scales, the (2)H longitudinal relaxation dispersion is measured for a polycrystalline BPTI sample. The observed dispersion is dominated by rapidly exchanging deuterons in protein side chains, undergoing restricted rotational motions on a time scale of 10 ns.     

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.