Alcohols dehydrate lipid membranes: an infrared study on hydrogen bonding.

The effects of alcohols (methanol, ethanol, and n-butanol) on the hydrogen bonding of dipalmitoylphosphatidylcholine (DPPC) were studied by Fourier-transform infrared spectroscopy (FTIR) in water-in-oil (carbon tetrachloride) reversed micelles. The bound O-H stretching mode of water, bonded to DPPC, appeared as a broad band at around 3400 cm-1. The O-H bending mode of this complex appeared as a weak broad band at 1644 cm-1. No free O-H signal was observed. When alcohols were added, a part of DPPC-bound water was replaced by the alcohols. The released 'free' water appeared at 3680 cm-1. This free O-H stretching band represents water-alcohol complex. A new broad band of O-H stretching appeared at 3235 cm-1, which represents the alcohol molecules bound to the phosphate moiety of DPPC. When the alcohol concentration was increased, the intensities of the free O-H stretching and bending bands increased. The P = O- antisymmetric stretching band at 1238 cm-1 became broader and shifted to lower frequencies. This means that alcohols interacted with the phosphate moiety and replaced the bound water. In the deconvoluted spectra of the C = O stretching mode, the ratio between the free sn-2 and the hydrogen-bonded sn-2 bands increased; a part of the bound water at the sn-2 carbon in the glycerol skeleton is also released and the free sn-2 signal increased. From the change in the intensity of the P = O- stretching band, the partition coefficients of alcohols between the phosphate region of DPPC and water were estimated: methanol 7.8, ethanol 16.7 at 22.0 degrees C in mole fraction bases. In molality, these values translates into methanol 0.21 and ethanol 0.45. These results indicate that short-chain alcohols interact with lipid membranes at the phosphate moiety at the hydrophilic head, weaken the membrane-water interaction, and destabilize membranes.     

Water structural changes in the bacteriorhodopsin photocycle: analysis by Fourier transform infrared spectroscopy.

The Fourier transform infrared difference spectra between light-adapted bacteriorhodopsin (BR) and its photointermediates, L and M, were analyzed for the 3750-3450-cm-1 region. The O-H stretching vibrational bands were identified from spectra upon substitution with 2H2O. Among them, the 3642-cm-1 band of BR was assigned to water by substitution with H2(18)O. By a comparison with the published infrared spectra of the water in model systems [Mohr, S.C., Wilk, W.D., & Barrow, G.M. (1965) J. Am. Chem. Soc. 87, 3048-3052], it is shown that the O-H bonds of the water in BR interact very weakly. Upon formation of L, the interaction becomes stronger. The O-H bonds of the protein side chain undergo similar changes. On the other hand, M formation further weakens the interaction of the same water molecules in BR. The appearance of a sharp band at 3486 cm-1, which was assigned tentatively to the N-H stretching vibration of the peptide bond, is unique to L. The results suggest that the water molecules are involved in the perturbation of Asp-96 in the L intermediate and that they are exerted from the protonated Schiff base which changes position upon the light-induced reaction.     

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)     

A Raman scattering study of the interactions of DNA with its water of hydration

Raman spectroscopy is used to probe the nature of the hydrogen bonds which hold the water of hydration to DNA. The ～ 3450?cm^?1 molecular O–H stretching mode shows that the first six water molecules per base pair of the primary hydration shell are very strongly bound to the DNA. The observed shift in the peak position of this mode permits a determination of the length of the hydrogen bonds for these water molecules. These hydrogen bonds appear to be about 0.3?Å shorter than the hydrogen bonds in bulk water. The linewidth of this mode shows no significant changes above water contents of about 15 water molecules per base pair. This technique of using a vibrational spectroscopy to obtain structural information about the hydration shells of DNA could be used to study the hydration shells of other biomolecules.     

Resonance Raman spectroscopy of ribonucleotide reductase. Evidence for a deprotonated tyrosyl radical and photochemistry of the binuclear iron center

Native ribonucleotide reductase from Escherichia coli exhibits a resonance-enhanced Raman mode at 1498 cm-1 that is characteristic of a tyrosyl radical. The Raman frequency as well as the absorption maximum at 410 nm identifies the radical as being in a deprotonated state. The B2 subunit of ribonucleotide reductase shows an additional resonance Raman mode at 493 cm-1 that has been assigned to the symmetric stretch of an Fe-O-Fe moiety. When samples of active B2 or metB2 are exposed to a tightly focused laser beam at 406.7 nm, there is a loss of intensity at 493 cm-1 and the appearance of a new peak at 595 cm-1. Although the 595-cm-1 feature was previously assigned to an Fe-OH vibration on the basis of its 23-cm-1 shift to lower energy in H2(18)O and the apparent dependence of its intensity on pH [Sj?berg, B. M., Loehr, T. M., & Sanders-Loehr, J. (1987) Biochemistry 26, 4242], the present studies indicate that the intensity of this mode is dependent primarily on input laser power. The peak at 595 cm-1 is more plausibly assigned to a new vs(Fe-O-Fe) mode in view of its lack of the deuterium isotope dependence expected for an Fe-OH mode and its resonant scattering cross section which is comparable to that of the 493-cm-1 mode. This new species has a calculated Fe-O-Fe angle of approximately 113 degrees compared to approximately 138 degrees calculated for the Fe-O-Fe unit in unmodified protein B2. One possible explanation for the photoinduced vibrational mode is that a bridging solvent molecule has been inserted in place of a bridging carboxylate.     

Raman spectroscopic study of left-handed Z-RNA

The solvent conditions that induce the formation of a left-handed Z form of poly[r(G-C)] have been extended to include 6.5 M NaBr at 35 degrees C and 3.8 M MgCl2 at room temperature. The analysis of the A----Z transition in RNA by circular dichroism (CD), 1H and 31P NMR, and Raman spectroscopy shows that two distinct forms of left-handed RNA exist. The ZR-RNA structure forms in high concentrations of NaBr and NaClO4 and exhibits a unique CD signature. ZD-RNA is found in concentrated MgCl2 and has a CD signature similar to the Z form of poly[d(G-C)]. The loss of Raman intensity of the 813-cm-1 A-form marker band in both the A----ZR-RNA and A----ZD-RNA transitions parallels the loss of intensity at 835 cm-1 in the B----Z transition of DNA. A guanine vibration that is sensitive to the glycosyl torsion angle shifts from 671 cm-1 in A-RNA to 641 cm-1 in both ZD- and ZR-RNA, similar to the B----Z transition in DNA in which this band shifts from 682 to 625 cm-1. Significant differences in the glycosyl angle and sugar pucker between Z-DNA and Z-RNA are suggested by the 16-cm-1 difference in the position of this band. The Raman evidence for structural difference between ZD- and ZR-RNA comes from two groups of bands: First, Raman intensities between 1180 and 1600 cm-1 of ZD-RNA differ from those for ZR-RNA, corroborating the CD evidence for differences in base-stacking geometry. Second, the phosphodiester stretching bands near 815 cm-1 provide evidence of differences in backbone geometry between ZD- and ZR-RNA.     

Ionic strength dependence of cytochrome c structure and Trp-59 H/D exchange from ultraviolet resonance Raman spectroscopy

Ultraviolet resonance Raman spectra are reported for cytochrome c (cyt c) in FeII and FeIII oxidation states at low (0.005 M) and high (0.9-1.5 M) ionic strength. With 200-nm excitation the amide band intensities are shown to remain constant, establishing that redox state and ionic strength have no influence on the alpha-helical content. The tyrosine 830/850-cm-1 doublet, however, shows a loss in 830-cm-1 intensity at I = 0.005 M for the FeIII protein, suggesting a weakening or a loss of H-bonding from an internal tyrosine, probably Tyr-48, which is H-bonded to a heme propionate group in cyt c crystals. Excitation profiles of tryptophan peak at approximately 229 nm for both FeII and FeIII forms of cyt c, but at approximately 218 nm for aqueous tryptophan. The approximately 2200-cm-1 red shift of the resonant electronic transition is attributed to the Trp-59 residue being buried and H-bonded. Consistent with this Trp environment, the H-bond-sensitive 877-cm-1 Trp band is strong and sharp, and the 1357/1341-cm-1 doublet has a large intensity ratio, approximately 1.5, for both FeII and FeIII cyt c. The 877-cm-1-band frequency shifts to 860 cm-1 when the Trp indole proton is replaced by a deuteron. This band was used to show that Trp H/D exchange in D2O is much faster for FeIII than FeII cyt c. The half-time for exchange at room temperature is estimated to be approximately 30 and approximately 5 h, respectively, for FeII and FeIII when examined at I = 0.005.(ABSTRACT TRUNCATED AT 250 WORDS)     

The A--B transition: temperature and base composition effects on hydration of DNA

Natural DNAs and some polynucleotides organised in fiber present the A--B form transition at a relative humidity (r.h.) which depends on the temperature. A shift of the midpoint of that helix--helix transition to higher r.h. values is observed when the temperature is risen. It is shown that the average number of water molecules associated to a nucleotide pair is the relevant parameter for the A-B transition and that this parameter can be given a precise value by a combination of different r.h. and temperature values. The minimum number of water molecules necessary to get the B form depends on the base composition of the DNA. It is observed that AT base pairs have a higher affinity toward water molecules than GC base pairs. In the B form there are 27 water molecules per GC nucleotide pair and 44 per AT pair. Moreover, we noted that the fraction of nucleotides in the B form as a function of the average number of water molecules associated per base pair does not depend on the temperature. The A helical form is obtained with about 11 water molecules per nucleotide pair and this number is not very sensitive to the base composition of DNA.     

Water structural changes in the L and M photocycle intermediates of bacteriorhodopsin as revealed by time-resolved step-scan Fourier transform infrared (FTIR) spectroscopy

In previous Fourier transform infrared (FTIR) studies of the photocycle intermediates of bacteriorhodopsin at cryogenic temperatures, water molecules were observed in the L intermediate, in the region surrounded by protein residues between the Schiff base and Asp96. In the M intermediate, the water molecules had moved away toward the Phe219-Thr46 region. To evaluate the relevance of this scheme at room temperature, time-resolved FTIR difference spectra of bacteriorhodopsin, including the water O-H stretching vibration frequency regions, were recorded in the micro- and millisecond time ranges. Vibrational changes of weakly hydrogen-bonded water molecules were observed in L, M, and N. In each of these intermediates, the depletion of a water O-H stretching vibration at 3645 cm-1, originating from the initial unphotolyzed bacteriorhodopsin, was observed as a trough in the difference spectrum. This vibration is due to the dangling O-H group of a water molecule, which interacts with Asp85, and its absence in each of these intermediates indicates that there is perturbation of this O-H group. The formation of M is accompanied by the appearance of water O-H stretching vibrations at 3670 and 3657 cm-1, the latter of which persists to N. The 3670 cm-1 band of M is due to water molecules present in the region surrounded by Thr46, Asp96, and Phe219. The formation of L at 298 K is accompanied by the perturbations of Asp96 and the Schiff base, although in different ways from what is observed at 170 K. Changes in a broad water vibrational feature, centered around 3610 cm-1, are kinetically correlated with the L-M transition. These results imply that, even at room temperature, water molecules interact with Asp96 and the Schiff base in L, although with a less rigid structure than at cryogenic temperatures.     

Identification of a hydroxide ligand at the iron center of ribonucleotide reductase by resonance Raman spectroscopy

The resonance Raman spectrum of protein B2 of ribonucleotide reductase from Escherichia coli shows several features to its oxo-bridged binuclear iron center. A peak at 492 cm-1 is assigned to the symmetric stretch of the Fe-O-Fe moiety on the basis of its 13-cm-1 shift to lower energy upon 18O substitution. The 18O species shows an additional peak at 731 cm-1, which is a good candidate for the asymmetric stretch of the Fe-O-Fe moiety. Its exact location in the 16O species is obscured by the presence of a protein tryptophan vibration at 758 cm-1. A third resonance-enhanced peak at 598 cm-1 is identified as an Fe-OH vibration on the basis of its 24-cm-1 shift to lower energy in H2 18O, its 2-cm-1 shift to lower energy in D2O, and its pH-dependent intensity. A hydrogen-bonded mu-oxo bridge similar to that in hemerythrin is suggested by the unusually low frequency for the Fe-O-Fe symmetric stretch and the 3-cm-1 shift to higher energy of vs(Fe-O-Fe) in D2O. From the oxygen isotope dependence of vs(Fe-O-Fe), an Fe-O-Fe angle of 138 degrees can be calculated. This small angle suggests that the iron center consists of a tribridged core as in hemerythrin. A model for the binuclear iron center of ribonucleotide reductase is presented in which the hydroxide ligand sites provide an explanation for the half-of-sites reactivity of the enzyme.     

Ion solvation and water structure in potassium halide aqueous solutions

The structure of water and the nature of ionic hydration is explored in aqueous solutions of potassium fluoride, chloride, bromide and iodide over a range of concentrations up to 4.8 ion pairs per 100 water molecules, using the combined techniques of neutron diffraction with hydrogen isotope substitution. The diffraction data are interpreted using the method of empirical potential structure refinement, which attempts to build a three-dimensional model of the scattering system consistent with the diffraction data. The water structure is strongly perturbed in the first hydration shells of both anion and cation, but is found to be only mildly perturbed outside of this region, with the largest effects occurring with the smallest anion and highest concentrations. For the potassium ion there are strong orientational correlations in the first hydration shell, with the water molecules lying with their dipole moments pointing almost directly away from the cation on average, but with an angular spread of approximately +/-60 degrees which is mildly dependent on the anion type present. For all the anions the water molecules in the first shell are strongly oriented with one O-H vector pointing directly towards the anion on average, with an angular spread of approximately +/-10 degrees for F(-), increasing to approximately +/-22 degrees for I(-). For both anions and cations the second hydration shell is much more disordered than the first, but there is a weak pattern of orientational correlation which becomes more pronounced with the larger anions. There is some evidence that the fluoride ion structures water significantly in its first hydration shell, but not beyond. The findings throw further light on recent findings that the orientational relaxation time for water outside the first shell of dissolved ions is the same as in the bulk liquid.     

Characteristics and variations of B-type DNA conformations in solution: a quantitative analysis of Raman band intensities of eight DNAs

Raman spectra were obtained from four bacterial DNAs varying in GC content and four periodic DNA polymers in 0.1 M NaCl at 25 degrees C. A curve fitting procedure was employed to quantify and compare Raman band characteristics (peak location, height, and width) from 400 to 1600 cm-1. This procedure enabled us to determine the minimum number of Raman bands in regions with overlapping peaks. Quantitative comparison of the Raman bands of the eight DNAs provided several new results. All of the DNAs examined required bands near 809 (+/- 7) and 835 (+/- 5) cm-1 to accurately reproduce the experimental spectra. Since bands at these frequencies are associated with A-family and B-family conformations, respectively, this result indicates that all DNAs in solution have a mixture of conformations on the time scale of the Raman scattering process. Band characteristics in the 800-850-cm-1 region exhibited some dependence on CG content and base pair sequence. As previously noted by Thomas and Peticolas [Thomas, G. A., & Peticolas, W. L. (1983) J. Am. Chem. Soc. 105, 993], the poly[d(A)].poly[d(T)] spectra were qualitatively distinct in this region. The A-family band is clearly observed at 816 cm-1. The intensity of this band and that of the B-family band at 841 cm-1 were similar, however, to intensities in the natural DNA spectra. Three bands at 811, 823, and 841 cm-1 were required to reproduce the 800-850-cm-1 region of the poly[d(A-T)].poly[d(A-T)] spectra. This may indicate the presence of three backbone conformations in this DNA polymer. Analysis of intensity vs. GC content for 42 Raman bands confirmed previous assignments of base and backbone vibrations and provided additional information on a number of bands.     

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.     

Hydration of Escherichia coli lipids. Deuterium T1 relaxation time studies of phosphatidylglycerol, phosphatidylethanolamine and phosphatidylcholine

The hydration properties of Escherichia coli lipids (phosphatidylglycerol, phosphatidylethanolamine) and synthetic 1,2-dioleoyl-sn-glycero-3-phosphocholine in H2O/2H2O mixtures (9:1, v/v) were investigated with 2H-NMR. Comparison of the 2H2O spin lattice relaxation time (T1) as a function of the water content revealed a remarkable quantitative similarity of all three lipid-H2O systems. Two distinct hydration regions could be discerned in the T1 relaxation time profile. (1) A minimum of 11-16 water molecules was needed to form a primary hydration shell, characterized by an average relaxation time of T1 approximately equal to 90 ms. (2) Additional water was found to be in exchange with the primary hydration shell. The exchange process could be described in terms of a two-site exchange model, assuming rapid exchange between bulk water with T1 = 500 ms and hydration water with T1 = 80-120 ms. Analysis of the linewidth and the residual quadrupole splitting (at low water content) confirmed the size of the primary hydration layer. However, each lipid-water system exhibited a somewhat different linewidth behavior, and a detailed molecular interpretation appeared to be preposterous.     

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 of single-stranded phosphodiester and phosphorothioate oligodeoxyribonucleotides.

Infrared spectroscopy was used to identify hydration-sensitive structural differences between single- stranded phosphorothioate (PS) and phosphodiester (PO) oligodeoxyribonucleotides. Spectra were recorded in the mid-infrared region, 500-1800 cm-1, at relative humidities between 0 and 98%; the PS and PO spectra are substantially different. The hydration effects on spectral bands in these single-stranded oligodeoxyribonucleotides is markedly different from such behavior in double- and triple-stranded oligodeoxyribonucleotides. A strong absorption occurs at 656 cm-1 in the phosphorothioate sample which is completely absent from the PO spectra. Gravimetric measurements were carried out on one PS and one PO sample to monitor and confirm hydration. The calculated BET adsorption constants [Brunauer, S., Emmett, RH. and Teller, E. (1938) J. Am. Chem. Soc., 60, 309-319] are 1.2 and 1.4 water molecules per nucleotide in the first hydration layer of PS and PO respectively. While the gravimetric data indicate that the single-stranded oligodeoxyribonucleotides hydrate very similarly to duplex DNA, the mid-infrared conformational marker bands are strikingly different from those observed for duplex DNA. In particular, the Vas of the phosphate group (PO2) at 1222 cm-1 in the single-stranded PO spectra is independent of relative humidity.     

Low-frequency Raman spectra of lysozyme crystals and oriented DNA films: dynamics of crystal water.

We observed low-frequency Raman spectra of tetragonal lysozyme crystals and DNA films, with varying water content of the samples. The spectra are fitted well by sums of relaxation modes and damped harmonic oscillators in the region from approximately 1 cm(-1) to 250 cm(-1). The relaxation modes are due to crystal water, and the distribution of relaxation times is determined. In wet samples, the relaxation time of a small part of the water molecules is a little longer than that of bulk water. The relaxation time of a considerable part of the crystal water, which belongs mainly to the secondary hydration shell, is an order of magnitude longer than that of bulk water. Furthermore, the relaxation time of some water molecules in the primary hydration shell of semidry samples is shorter than we expected. Thus we have shown that low-frequency Raman measurements combined with properly oriented samples can give specific information on the dynamics of hydration water in the ps range. On the other hand, we concluded, based on polarized Raman spectra of lysozyme crystals, that the damped oscillators correspond to essentially intramolecular vibrational modes.     

Two-dimensional/ATR infrared correlation spectroscopic study on water diffusion in a poly(epsilon-caprolactone) matrix

In the present contribution, two-dimensional ATR-FTIR spectroscopy was used to study the diffusion of water in poly(epsilon-caprolactone) (PCL). In the spectral region of the nu(OH) stretching vibration of water, four absorption bands (3635, 3560, 3411, and 3220 cm(-1)) can be identified. The higher wavenumber band pair at 3635 and 3560 cm(-1) is assigned to the antisymmetric and symmetric OH stretching vibrations, respectively, of water which is partially hydrogen-bonded to the carbonyl groups of PCL, whereas the lower frequency band pair at 3411 (antisymmetric) and 3220 cm(-1) (symmetric) is attributed to the OH stretching vibrations of bulk water which is fully hydrogen-bonded to other water molecules. From the asynchronous map of a 2D correlation analysis of spectra recorded during the diffusion of water into PCL, it was concluded that during the diffusion process the water molecules first penetrate into the free volume (microvoids) of the PCL matrix or are molecularly dispersed in the polymer matrix and then form hydrogen bonds with the C=O groups of the polymer.     

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.     

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

Monte Carlo simulations [(N, V, T)-ensemble] were performed for the hydration shell of poly(dA-dT).poly(dA-dT) in canonical B form and for the hydration shell of poly(dA).poly(dT) in canonical B conformation and in a conformation with narrow minor groove, highly inclined bases, but with a nearly zero-inclined base pair plane (B' conformation). We introduced helical periodic boundary conditions with a rather small unit cell and a limited number of water molecules to reduce the dimensionality of the configuration space. The coordinates of local maxima of water density and the properties of one- and two-membered water bridges between polar groups of the DNA were obtained. The AT-alternating duplex hydration mirrors the dyad symmetry of polar group distribution. At the dApdT step, a water bridge between the two carbonyl oxygens O2 of thymines is formed as in the central base-pair step of Dickerson's dodecamer. In the major groove, 5-membered water chains along the tetranucleotide pattern d(TATA).d(TATA) are observed. The hydration geometry of poly(dA).poly(dT) in canonical B conformation is distinguished by autonomous primary hydration of the base-pair edges in both grooves. When this polymer adopts a conformation with highly inclined bases and narrow minor groove, the water density distribution in the minor groove is in excellent agreement with Dickerson's spine model. One local maximum per base pair of the first layer is located near the dyad axis between adjacent base pairs, and one local maximum per base pair in the second shell lies near the dyad axis of the base pair itself. The water bridge between the two strands formed within the first layer was observed with high probability. But the water molecules of the second layer do not have a statistically favored orientation necessary for bridging first layer waters. In the major groove, the hydration geometry of the (A.T) base-pair edge resembles the main features of the AT-pair hydration derived from other sequences for the canonical B form. The preference of the B' conformation for oligo(dA).oligo(dT) tracts may express the tendency to common hydration of base-pair edges of successive base pairs in the grooves of B-type DNA. The mean potential energy of hydration of canonical B-DNA was estimated to be -60 to -80 kJ/mole nucleotides in dependence on the (G.C) contents. Because of the small system size, this estimation is preliminary.