| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |

* Centro de Investigação em Química (Universidade do Porto), Department of Chemistry, Faculty of Sciences, University of Porto, P-4169-007 Porto, Portugal; and
Laboratoire Léon Brillouin (Commissariat à I'Énergie Atomique-Centre National de la Recherche Scientific), CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France
Correspondence: Address reprint requests to Margarida Bastos, Dept. of Chemistry, R. Campo Alegre, 687 P-4169-007 Porto, Portugal. Tel.: 351-22-6082811; E-mail: mbastos{at}fc.up.pt.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
1 ps ("bulk-like") and another
1012 ps (weakly bound type). The authors stress the importance of hydration water in structural integrity maintenance and recognition. It has also been pointed out in another work that the different hydration found for the different motifs for the tandem wobble G U mismatches may reflect that hydration plays a role in their stability; and further, that their results imply that hydration of nucleic acids must be important in the genetic code recognition (Sundaralingam and Pan, 2002
Despite the increasing knowledge and the considerable progress made in the x-ray analysis of short oligonucleotides (Egli et al., 1998
), the complex surface of the double helix and its hydration are still the most poorly defined parts of the structures obtained by crystallographic methods, because they are topologically disordered or dynamically averaged and have conformational fluctuations (especially for the high molecular, semiflexible, highly charged DNA chains in aqueous solution). Hydration water is supposed to interact directly with phosphate groups, sugar oxygen atoms, and the polar groups of the bases. It has been shown that water molecules are mainly concentrated in six hydration sites per phosphate and that the positions and occupancies of these sites depend on the DNA conformation and on its nucleotide composition (Schneider et al., 1998
). Hydration is also very important regarding protein-DNA interactions, because it has been suggested that hydration sites mark the binding sites at protein-DNA interfaces (Woda et al., 1998
).
In this respect neutron quasielastic scattering (NQES) is probably the most powerful technique that can be used for the study of the mechanism of biopolymerswater interactions in hydrated powders, proteins, and DNA fibers (see, e.g., Schreiner et al., 1988
; Grimm and Rupprecht, 1989
; Bellissent-Funel et al., 1992
; Zanotti et al., 1997
; Foucat et al., 2002
). DNA hydration has also been studied by neutron diffraction, and chains of ordered water in DNA grooves were thereby identified (Shotton et al., 1997
).
In this article results of high-resolution NQES experiments are reported for hydrated samples of DNA in the double-stranded (ds) and single-stranded (ss) states, at equal guanine-cytosine content. According to previous calorimetric data, the transformation "double-stranded helix"
"single-stranded chains" was found to be accompanied by "melting" of the structural part of the hydration shell of ds-DNA (Mrevlishvili et al., 2001
, 2002
). Moreover, formation of the duplex from mixing of their complementary single strands was found to be mostly accompanied by the uptake of structural water molecules (Mrevlishvili et al., 2002
; Kankia and Marky, 1999
). We wanted then to study ds- and ss-DNA, at low, limiting hydration levels by NQES, to see if we could characterize the type of water present in each case, as well as to find whether a difference could be found between ds- and ss-DNA as regarding water behavior and mobility.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Sample preparation
Samples of ds-DNA and ss-DNA were hydrated by exposing them to a high relative humidity atmosphere, from a saturated solution of potassium sulfate (relative humidity 98.5% at 5°C). The initial water content of the "dry" samples was determined by coulometry (Karl Fischer titration), using a Karl Fischer coulometer, model XX (Metrohm, Herisau, Switzerland). This initial water content was found to be 18.8% (w/w) for ds-DNA and 14.4% (w/w) for ss-DNA. The samples were thereafter maintained in an exicator under the high relative humidity atmosphere and were weighted daily, until they reached the desired water content, 0.559 g H2O/g ds-DNA (B-DNA) and 0.434 g H2O/g ss-DNA. This final water content of each sample was calculated taking into account the amount of water already present in the starting sample.
In both cases, after the desired water content was attained, each sample was placed in a rectangular aluminum plate, which was the container to be used in the scattering experiments, (45 x 75 mm and 0.5-mm thickness) and tightly closed with eight screws, by use of an indium wire.
Neutron-scattering experiments
Data collection
Neutron-scattering experiments were performed at Laboratoire Léon Brillouin (Saclay, France) using a time-of-flight (TOF) spectrometer (MIBEMOL). In a neutron-scattering experiment, one measures the exchange of momentum, Q, and the exchange of energy,
, due to the interaction of the incident neutron beam with the nuclei of the sample. Momenta and energy conservation are expressed by:
![]() | (1) |
![]() | (2) |
Thermalized neutrons were used, with two different incident wavelengths, so as to obtain different resolutions. When the wavelength was 5.2 Å, a triangular resolution function with halfwidth at half-maximum (HWHM) of 64 µeV was obtained, whereas a value of HWHM of 17 µeV was obtained for the wavelength of 8 Å. The resolution function was in both cases Q independent. The aluminum plate (sample holder) was placed at an angle of 135° to the incident neutron beam. The scattered neutrons were detected at angles from 23.5° to 141.8°. The corresponding range of accessible values of the elastic momentum exchange Q was 0.3461.254 Å1 (8 Å) and 0.5331.929 Å1 (5.2 Å). Sampling times of 8 h were used for ds-DNA. For ss-DNA, the small amount of sample required a larger sampling time to improve the statistics, and therefore 16-h sampling times were used. The measurements were performed at 298 K. The resolution at these two wavelengths was determined by running the spectra of vanadium. The spectrum of the empty container was also recorded, and used to correct all experiments.
Data analysis
The TOF spectra were analyzed with programs, available at Laboratoire Léon Brillouin, that allow calibration of the detectors with the vanadium spectra, correction for the empty cell, transformation of the TOF spectra into energy spectra, and finally data grouping (GROUPE software) to improve signal/noise ratio. The data so treated was thereafter analyzed with the FITMIB software, from which the scattering laws S(Q,
) could be obtained. The protons that appear as static within the spectrometer resolution give rise to an elastic contribution. In addition, mobile protons with similar dynamics appear in the spectra as Lorentzians, characterized by a HWHM. The scattering function, S(Q,
), at each Q and
can then be written as the sum of elastic and broadened terms:
![]() | (3) |
(
), weighted by the elastic incoherent structure factor (EISF) Ao(Qa), and the second term is in principle an infinite sum of Lorentzians. In practice, due to the finite resolution of the instrument and to the experimental statistics, it can be approximated by a single Lorentzian having a HWHM,
1/2
![]() | (4) |
The program FITMIB provides both the HWHM,
1/2 and the intensity of the Lorentzian, as well as the intensity of the elastic contribution.
From the experimentally obtained integrated intensities, the experimental EISF were calculated as
![]() | (5) |
This general formulation can be extended to extract more information from the obtained data. The goal is to characterize solvent dynamics, trying to discriminate between different types of solvent motion, and to quantify the respective fraction of protons. According to the model of Volino and Dianoux for the diffusion of a point particle inside a sphere of radius a (Volino and Dianoux, 1980
) the EISF can be calculated according to:
![]() | (6) |
![]() | (7) |
After resolution broadening by a resolution function R(
),
![]() | (8) |
means convolution in
space. If the water molecules strongly attached to the DNA surface have a Lorentzian line narrower than the resolution function R(
), then to a good approximation we can write
![]() | (9) |
![]() | (10) |
The obtained scattering results were fitted to this model. The parameters
1/2 as a function of Q were extracted by use of the FITMIB program, and the experimental EISF values calculated (Eq. 5) from the integrated intensities also obtained from FITMIB. The calculation of the EISF values, coupled with the observation of the profile of the HWHM of the Lorentzian, allows the analysis of the geometry of the movements that give rise to the quasielastic broadening. The experimental values of EISF were used along with the ones that can be calculated from the model of Volino and Dianoux (1980)
(Eq. 6) to obtain a and P by minimization (nonlinear minimum square fitting), as from Eq. 9 we can see that
![]() | (11) |
After obtaining the value of P, the parameter q (fraction of protons (from hydration water) strongly attached to the DNA surface), could also be calculated, because p (fraction of immobile protons in DNA) was known from the sample stoichiometry (see Eq. 10 above).
| RESULTS |
|---|
|
|
|---|
1/2 value was obtained, which was Q independent. Further, the
1/2 values for the narrow Lorentzian were identical, within error. It turned out, then, that for this analysis, either could be used. The results presented correspond to the fit obtained with two Lorentzians at 5.2 Å, because the EISF values obtained with this model fit perfectly with those obtained at 8 Å.
In Fig. 1 one example of obtained fits for ds-DNA is presented, at 8 Å, at a scattering angle
= 23.5° (for 8 Å). These results are representative of the general goodness of fit obtained at the other scattering angles and wavelengths.
|
1/2 for the narrow Lorentzian at 5.2 Å are represented as a function of Q. As can be seen in the figure,
1/2 presents, within error bars, a plateau value of 0.080 meV until 1.12 Å1. This feature can be accounted for by the model of Volino and Dianoux referred to above as being due to confinements, and represented by diffusion within a sphere (Volino and Dianoux, 1980
/a, allows the calculation of a, the radius of the sphere, as
![]() | (12) |
|
![]() | (13) |
As explained in the data-analysis section, the experimental values of EISF (Eq. 5) along with the ones that can be calculated from the model of Volino and Dianoux (1980)
(Eq. 6) were used to obtain a and P by minimization (nonlinear minimum square fitting) We did the fitting either by fixing a to the value of 2.8 Å obtained above and just fitting for the value of P, or by allowing both parameters to vary. The obtained results were in perfect agreement. The P-value obtained for ds-DNA was 0.51 and a = 2.8 Å.
The experimental values of EISF, at 8 and 5.2 Å, as well as the calculated EISF according to Eq. 6 (using the parameters a and P just obtained) are plotted in Fig. 3, as a function of Q.
|
The same initial data treatment with FITMIB was applied to the results obtained at 5.2 Å for ss-DNA, and the
1/2 values as a function of Q as well as the integrated intensities were obtained. In Fig. 4 we can see the plot of
1/2 as a function of Q and in Fig. 5 the EISF(exp) values are represented also as a function of Q. Indeed, the error bars are much larger for ss-DNA, as the amount of substance in the sample was reduced, and although we used double sampling time (16 h), we have rather poor statistics. Nevertheless, we can say that no momentum dependence is observed for the obtained
1/2 values. The same is apparent from the EISF values plotted on Fig. 5. Even so, we did the minimization of the EISF to obtain the P-value (Eq. 11), considering a = 2.8 Å, just for the sake of comparison with the ds-DNA results. The value obtained was P = 0.78, which equals, within uncertainty, the p-value calculated from the sample stoichiometry, p = 0.76.
|
|
| DISCUSSION |
|---|
|
|
|---|
30% (v/v) of solvent, which corresponds to
0.38 g H2O/g protein) at 300 K (Bon et al., 2002
Regarding the water diffusion, our value D = 2.2 x 105 cm2 s1 is similar to the value for bulk water at this temperature, D = 2.3 x 105 cm2 s1 (Teixeira et al., 1985
). It is also similar to the value D = 2.6 x 105 cm2 s1 obtained for fully hydrated C-phycocyanin at 300 K (Bellissent-Funel et al., 1992
) and to the range of values estimated for ds-DNA with water contents between 0.5 and 2.5 g H2O/g ds-DNA, 1.44.2 x 105 cm2 s1 (Beta et al., 2003
), and larger than the values D = 0.68 x 105 cm2 s1 obtained for parvalbumin with 0.65 g D2O/g protein, at 298 K (Zanotti et al., 1999
), D = 0.39 x 105 cm2 s1 for hen egg-white lyzozyme crystals with 0.38 g H2O/g protein (Bon et al., 2002
) and D = 0.5 x 105 cm2 s1 for DNA fibers. Regarding this last value, we should note that they have used a lower hydration regime (9.5 mol H2O/mol bp whereas we have 20 mol H2O/mol bp), which reflects in a lower limiting
1/2 value (20 µeV) and should correspond to a lower diffusion coefficient. If we except the value obtained for parvalbumin, which is very low taking into account the water content of the sample, the other results, including ours, are not surprising, as it is known that only with degrees of hydration of
0.6 (w/w) the properties of interfacial water become close to bulk solvent properties (Goldanskii and Krupanskii., 1995
).
The use of a refined model, introducing the parameters p and q, allowed us to distinguish between two types of water present in ds-DNA. Our results indicate that a fraction q = 0.28 of the water protons are strongly attached to the ds-DNA surface, whereas the remaining water has a motion that can be interpreted within the model diffusive motion confined within a sphere of radius 2.8 Å. Very interestingly, a recent study of water at DNA surfaces by femtosecond-resolved fluorescence (Pal et al., 2003b
) also showed the presence of two types of water, that the authors call a "bimodal hydration behavior", with two distinct time constants:
1 ps (that they classify as bulk type, labile water, that represents 60% of the total water) and 20 ps (classified by the authors as weakly bound, ordered water, representing 40% of total water). Not only have they observed the same behavior as we find here, namely of two types of water with different dynamics, but also our percentages are close to theirs; our value for the percentage of protons strongly attached to the ds-DNA surface (q = 0.28) is
30% and corresponds to their fraction of 40% of "weakly bound, ordered water."
In the samples of ss-DNA, on the contrary, we could not find any indication of confinement, from the momentum transfer dependences of either the EISF or
1/2. This shows that the water present, although contributing to the quasielastic peak, does not present the typical behavior of confined water.
It is very enlightening to compare these results with the ones obtained recently by one of us (G.M.) in a differential scanning calorimetry study of low hydration samples of ds-DNA and ss-DNA (Mrevlishvili at al., 2001
,2002
). In these studies, the authors use the concept of "unfreezable water", that is, the amount of water present in the DNA sample that does not undergo the liquid-to-solid phase transition when the sample is cooled down in a DSC experiment. This amount was found to be 0.55 g H2O/g DNA for ds-DNA and 0.40 g H2O/g DNA for ss-DNA. We chose our water contents for the scattering experiments accordingly. In the interpretation of the obtained calorimetric data, the transformation "double-stranded helix"
"single-stranded chains" is suggested to be accompanied by dehydration of ds-DNA, which amounts to 0.15 g H2O/g DNA. Further, this value is equal to the one obtained in the first study (Mrevlishvili et al., 2001
) for the hydration change upon thermal denaturation of ds-DNA. This was interpreted as suggesting that the "double helix"
"single-stranded" thermal transition is accompanied by disruption of the ordered part of the water present in the hydration shell of the double helix. Now, if we consider the result obtained in this study by neutron scattering, we did obtain a value for q, the fraction of strongly attached protons to the ds-DNA surface, of 0.28. This value corresponds to 0.16 g H2O/g DNA. As referred to above, in the case of ss-DNA no evidence was found of strongly attached protons, or better said, we could not distinguish different types of protons in the ss-DNA sample. This shows that we have obtained, by a completely independent measurement, a confirmation of the previous results obtained by calorimetry.
| CONCLUSION |
|---|
|
|
|---|
70%), within a sphere of radius 2.8 Å, with local diffusion coefficient D = 2.2 x 105 cm2 s1, that is, slightly smaller than bulk water (probably water inside the grooves), and another type of water, attached to the ds-DNA surface (
30%). This second type of water could be inside or outside the grooves. This water is lost upon ds-DNA thermal denaturation, as opposite to the first one, that is kept even after thermal denaturation. The fact that we found this second type of water to be absent in ss-DNA, and that the numerical values obtained here for the difference do agree with the ones obtained by calorimetry, strengthens the interpretation put forward before, i.e., that the thermal transition between ds-DNA and ss-DNA is accompanied by disruption of the ordered water fraction in the inner hydration of the double helix, and stresses the importance of hydration effects on the maintenance of the double-helix structure. | ACKNOWLEDGEMENTS |
|---|
|
|
|---|
M.B. and V.C. also thank, for the economical support at Laboratoire Léon Brillouin, the European Union through the Access to Research Infrastructures action of the Human Potential Program Access to Research Infrastructures. We acknowledge Fundação para a Ciência e Tecnologia for financial support to Centro de Investigação em Química (Universidad do Porto), Unidade de Investigação 81.
| FOOTNOTES |
|---|
Submitted on January 5, 2004; accepted for publication February 10, 2004.
| REFERENCES |
|---|
|
|
|---|
Bellissent-Funel, M.-C., J. Teixeira, K. F. Bradley, and S. H. Chen. 1992. Dynamics of hydration water in protein. J. Phys. I. 2:9951001.[CrossRef]
Beta, I. A., I. Michalarias, R. C. Ford, J. C. Li, and M.-C. Bellissent-Funel. 2003. Quasi-elastic neutron scattering study of hydrated DNA. Chem. Phys. 292:451454.[CrossRef]
Blackburn, G. M., and M. J. Gait. 1990. Nucleic Acids in Chemistry and Biology. Oxford University Press, Oxford, UK.
Bloomfield, V. A., D. M. Crothers, and I. Tinoco, Jr. 2000. Nucleic Acids. University Science Books, Sausalito, CA.
Bon, C., A. J. Dianoux, M. Ferrand, and M. S. Lehmann. 2002. A model for water motion in crystals of lyzozyme based on an incoherent quasielastic neutron scaettering study. Biophys. J. 83:15781588.
Egli, M., V. Tereshko, M. Teplova, G. Minasov, A. Joachimiak, R. Sanishvili, C. M. Weeks, R. Miller, M. Maier, H. An, P. Dan Cook, and M. Manoharan. 1998. X-ray crystallographic analysis of the hydration of A- and B-form DNA at atomic resolution. Biopolymers. 48:234252.[CrossRef][Medline]
Foucat, L., J.-P. Renou, C. Tengroth, S. Janssen, and H. D. Middendorf. 2002. Dynamics of proteins at low temperatures: fibrous vs. globular. Appl. Phys. A. 74:S1290S1292.[CrossRef]
Goldanskii, V. I., and Y. F. Krupanskii. 1995. Protein Dynamics in Protein-Solvent Interactions. R. B. Gregory, editor. Marcel Dekker, New York/Basel, Switzerland/Hong Kong.
Grimm, H., and A. Rupprecht. 1989. Hydration structure in natural DNA observed by thermal neutron scattering. Eur. Biophys. J. 17:173186.[Medline]
Kankia, B., and L. Marky. 1999. DNA, RNA, and DNA/RNA oligomer duplex: a comparative study of their stability, heat, hydration and Mg binding properties. J. Phys. Chem. B. 103:87598767.
Liepinsh, E., G. Otting, and K. Wuthrich. 1992. NMR observation of individual molecules of hydration water bound to DNA duplexes: direct evidence for a spine of hydration water present in aqueous solution. Nucleic Acids Res. 20:65496553.
Mrevlishvili, G. M., A. P. S. M. C. Carvalho, M. A. V. Ribeiro da Silva, T. D. Mdzinarashvili, G. Z. Razmadze, and T. O. Tarielashvili. 2001. The role of bound water on the energetics of DNA duplex melting. J. Therm. Anal. Cal. 66:133144.[CrossRef]
Mrevlishvili, G. M., A. P. S. M. C. Carvalho, and M. A. V. Ribeiro da Silva. 2002. Low-temperature DSC study of the hydration of ss-DNA and ds-DNA and the role of hydrogen-bonded network to the duplex transition thermodynamics. Thermochim. Acta. 394:7382.[CrossRef]
Pal, S. K., L. Zhao, T. Xia, and A. H. Zewail. 2003a. Site- and sequence selective ultrafast hydration of DNA. Proc. Natl. Acad. Sci. USA. 100:1374613751.
Pal, S. K., L. Zhao, and A. H. Zewail. 2003b. Water at DNA surfaces: ultrafast dynamics in minor groove recognition. Proc. Natl. Acad. Sci. USA. 100:81138118.
Saenger, W. 1984. Principles of Nucleic Acid Structure. Springer-Verlag, New York/Berlin/Hidelberg/Tokyo.
Schneider, B., K. Patel, and H. M. Berman. 1998. Hydration of the phosphate group in double-helical DNA. Biophys. J. 75:24222434.
Schreiner, I. J., M. M. Pintar, A. J. Dianoux, P. Volino, and A. Rupprecht. 1988. Hydration of NaDNA by neutron quasi-elastic scattering. Biophys. J. 53:119122.
Shotton, M. W., L. H. Pope, T. Forsyth, P. Langan, R. C. Denny, U. Giesen, M.-T. Dauvergne, and W. Fuller. 1997. A high-angle neutron fibre diffraction study of the hydration of deuterated A-DNA. Biophys. Chem. 69:8596.[CrossRef][Medline]
Stent, G. S. 1995. The aperiodic crystal of heredity. Ann. N. Y. Acad. Sci. 758:2531.[Medline]
Sundaralingam, M., and B. Pan. 2002. Hydrogen and hydration of DNA and RNA nucleotides. Biophys. Chem. 95:273282.[CrossRef][Medline]
Teixeira, J., M.-C. Bellissent-Funel, S. H. Chen, and A. J. Dianoux. 1985. Experimental determination of the nature of diffusive motions of water molecules at low temperatures. Phys. Rev. A. 31:19131916.[CrossRef][Medline]
Volino, F., and A. J. Dianoux. 1980. Neutron incoherent scattering law for diffusion in a potential of spherical symmetry: general formalism and application to diffusion within a sphere. Mol. Phys. 41:271279.[CrossRef]
Woda, J., B. Schneider, K. Patel, K. Mistry, and H. M. Berman. 1998. An analysis of the relationship between hydration and protein-DNA interactions. Biophys. J. 75:21702177.
Zanotti, J.-M., M.-C. Bellissent-Funel, and J. Parello. 1997. Dynamics of a globular protein as studied by neutron scattering and solid-state NMR. Phys. B. 234236:228230.
Zanotti, J.-M., M.-C. Bellissent-Funel, and J. Parello. 1999. Hydration coupled dynamics in proteins studied by neutron scattering and NMR: the case of the typical EF-hand calcium-binding parvalbumin. Biophys. J. 76:23902411.
This article has been cited by other articles:
![]() |
J. Zhou, S. Krueger, and S. K. Gregurick A coarse graining approach to determine nucleic acid structures from small angle neutron scattering profiles in solution Nucleic Acids Res., November 10, 2005; 33(19): 6361 - 6371. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |