Ion dynamics and water percolation effects in DNA polymorphism

The dynamics of ions and water at the surface of DNA are studied by computer simulations in a wide range of hydrations involving the zone of low-hydration polymorphism in DNA. The long-range mobility of ions exhibits a stepwise increase at three distinct hydration levels. The first of them is close to the midpoint of the water percolation transition as well as the midpoint of the transition between A- and B-DNA forms. It coincides with the onset of the dissociation of ion pairs on the DNA surface probably caused by the increase in the water dielectric permittivity due to the appearance of the spanning hydrogen-bonding network. The other two steps are attributed to the formation of percolating water layers on the surface of DNA accompanied by the progressive escape of ions from the DNA surface. The results agree with earlier experimental data and further corroborate the suggested universal mechanism of the low hydration polymorphism in DNA including intraduplex electrostatic condensation close to the water percolation threshold.     

The Electrostatic Origin of Low‐Hydration Polymorphism in DNA

In recent years, significant progress has been made towards uncovering the physical mechanisms of low‐hydration polymorphism in double‐helical DNA. The effect appears to be mechanistically similar in different biological systems, and it is due to the ability of water to form spanning H‐bonded networks around biomacromolecules via a quasi‐two‐dimensional percolation transition. In the case of DNA, disintegration of the spanning H‐bonded network leads to electrostatic condensation of DNA strands because, below the percolation threshold, water loses its high dielectric permittivity, whereas the concentration of neutralizing counterions becomes high. In this Concept article arguments propose that this simple electrostatic mechanism represents the universal origin of low‐hydration polymorphism in DNA.     

Thermal breaking of spanning water networks in the hydration shell of proteins

The presence of a spanning hydrogen-bonded network of water at the surface of biomolecules is important for their conformational stability, dynamics, and function. We have studied by computer simulations the clustering and percolation of water in the hydration shell of a small elastinlike peptide (ELP) and the medium-size protein staphylococcal nuclease (SNase), in aqueous solution. We have found that in both systems a spanning network of hydration water exists at low temperatures and breaks up with increasing temperature via a quasi-two-dimensional percolation transition. The thermal breaking of the spanning water network occurs at biologically relevant temperatures, in the temperature range, which is close to the temperature of the "inverse temperature transition" of ELP and the unfolding temperature of SNase, respectively.     

Formation of spanning water networks on protein surfaces via 2D percolation transition

The formation of spanning hydrogen-bonded water networks on protein surfaces by a percolation transition is closely connected with the onset of their biological activity. To analyze the structure of the hydration water at this important threshold, we performed the first computer simulation study of the percolation transition of water in a model protein powder and on the surface of a single protein molecule. The formation of an infinite water network in the protein powder occurs as a 2D percolation transition at a critical hydration level, which is close to the values observed experimentally. The formation of a spanning 2D water network on a single rigid protein molecule can be described by adapting the cluster analysis of conventional percolation studies to the characterization of the connectivity of the hydration water on the surface of finite objects. Strong fluctuations of the surface water network are observed close to the percolation threshold. Our simulations also furnish a microscopic picture for understanding the specific values of the experimentally observed hydration levels, where different steps of increasing mobility in the hydrated powder are observed.     

Origin of the dynamic transition upon pressurization of crystalline proteins

We study the role of hydration water in the dynamic transition of low-hydrated proteins upon pressurization found recently (Meinhold, L.; Smith, J. C. Phys. Rev. E 2005, 72, 061908). Clustering and percolation of water in the hydration shells of protein molecules in crystalline Staphylococcal nuclease are analyzed at various pressures. The number of water molecules in the hydration shell increases and the hydrogen-bonded network of hydration water spans with increasing pressure. The dynamic transition of protein occurs when the spanning water network exists with the probability of about 50% and hydration water shows large density fluctuations. Formation of a spanning water network upon pressurization promotes protein dynamics as in the case of the dynamic transition with increasing hydration. Properties of hydration water in various thermodynamic states and their influence on biological function are discussed.     

Titration in silico of reversible B <=> A transitions in DNA

Reversible transitions between the A- and B-forms of DNA are obtained in free molecular dynamics simulations of a single double helix immersed in a water drop with Na(+) counterions. The dynamics of the transitions agrees with their supposed cooperative character. In silico titration of the transitions was carried out by smooth variation of the drop size. The estimated range of hydration numbers corresponding to the transition roughly agrees with experimental data. The chain length dependence was studied for double helices from 6 to 16 base pairs. It appeared that the B --> A transition is hindered for DNA shorter than one helical turn. With increased NaCl concentration in the drop, stabilization of the B-form is observed accompanied by the salt crystallization. The results strongly suggest that the B --> A transition at low hydration is caused by Na(+) ions sandwiched between phosphate strands in the major groove and is driven by direct medium range electrostatic interactions. The role of the reduced water shell apparently consists of increasing the counterion concentration in the opening of the major groove. Analysis of the available experimental data suggests that this mechanism is perhaps generally responsible for the A/B polymorphism in DNA.     

Which Properties of a Spanning Network of Hydration Water Enable Biological Functions?

The central role of water in biological functions is well‐recognized, but numerous questions concerning the physical mechanisms behind the importance of water for life remain unanswered. Water in biosystems exists mainly as hydration water. Analysis of the phase diagram of hydration water shows that biological functions are possible only when the surfaces of biomolecules are covered by spanning hydrogen‐bonded networks of hydration water. The comparative studies of the various properties of hydrated biosystems in the presence and in the absence of a spanning water network should clarify its specific physical properties, which are crucial for biological functions. Herein, we summarize the recent progress in these studies. The biological activity of the living organisms is maximal in a narrow temperature interval, where the spanning network of hydration water breaks up with heating via a percolation transition. The entropy of the hydration water related to the diversity of cluster size diverges at this percolation threshold. The possible role of this phenomenon in life processes is discussed.     

Iso-energy contour maps for the interaction of water with DNA double helix in the B conformation

The interaction energy between water with B-DNA double helix is computed for few cylindrical surfaces (enclosing the helix) using analytical pair potentials. The iso-energy contour maps indicate a strong attraction for water extending to three water layers surrounding DNA and very stable bridging structure of water molecules connecting two successive phosphate groups along a single helix in the innermost layer.     

A theoretical study on the hydration of B- and Z-DNA double helices

A theoretical study on the hydration of B- and Z-DNA double helices has been carried out using empirical potential energy functions. The interaction energy between water and the model compounds has been computed considering only the first hydration shell.The results show the number of binding water molecules to be thirty-six and twenty-five in B- and in Z-DNA, respectively. The water molecules in the first hydration shell of B-DNA are very well ordered along the phosphate groups of the backbone whereas those of Z-DNA are more disordered than in B-DNA and are more strongly bound. The water molecules near the first hydration shell of Z-DNA are thought to move more freely than those of B-DNA.     

Properties of spanning water networks at protein surfaces

The formation of a spanning two-dimensional hydrogen-bonded water network at the surface of proteins via a percolation transition enables their biological function. We show in detail how the spanning (percolating) water network appears at the surfaces of model hydrophilic spheres and at the surface of a single protein (lysozyme) molecule. We have found essential correlations of the linear extension, radius of gyration, and position of the center of mass of the largest water cluster with its size. The specific two-peak structure of the probability distribution of the largest cluster size allowed us to study various properties separately for spanning and nonspanning largest clusters. The radius of gyration of the spanning cluster always exceeds the radii of the spheres or the effective radius of the protein. Any spanning cluster envelops essentially more than half of the surface area. The temporal decay of the spanning networks shows a stretched exponential character. Their average lifetime at the percolation threshold is about the lifetime of a water-water hydrogen bond.     

The role of water H-bond imbalances in B-DNA substate transitions and peptide recognition revealed by time-resolved FTIR spectroscopy

The conformational substates B(I) and B(II) of the phosphodiester backbone in B-DNA are thought to contribute to DNA flexibility and protein recognition. We have studied by rapid scan FTIR spectroscopy the isothermal B(I)-B(II) transition on its intrinsic time scale. Correlation analysis of IR absorption changes occurring within seconds after a reversible incremental growth of the DNA hydration shell identifies water populations w(1) (PO(2)(-)-bound) and w(2) (non-PO(2)(-)-bound) exhibiting weaker and stronger H-bonds, respectively, than those dominating in bulk water. The B(II) substate is stabilized by w(2). The water H-bond imbalance of 3-4 kJ mol(-1) is equalized at little enthalpic cost upon formation of a contiguous water network (at 12-14 H(2)O molecules per DNA phosphate) of reduced ν(OH) bandwidth. In this state, hydration water cooperatively stabilizes the B(I) conformer via the entropically favored replacement of w(2)-DNA interactions by additional w(2)-water contacts, rather than binding to B(I)-specific hydration sites. Such water rearrangements contribute to the recognition of DNA by indolicidin, an antimicrobial 13-mer peptide from bovine neutrophils which, despite little intrinsic structure, preferentially binds to the B(I) conformer in a water-mediated induced fit. The FTIR spectra resolve sequential steps leading from PO(2)(-)-solvation to substate transition and eventually to base stacking changes in the complex. In combination with CD-spectral titrations, the data indicate that, in the absence of a bulk aqueous phase, as in molecular crowded environments, water relocation within the DNA hydration shell allows for entropic contributions similar to those assigned to water upon DNA ligand recognition in solution.     

Dynamics of surface water in ZrO2 studied by quasielastic neutron scattering

A quasielastic neutron scattering experiment has revealed the dynamics of surface water in a high surface area zirconium oxide in the temperature range of 300-360 K. The characteristic times of the rotational (picoseconds) and translational (tens of picoseconds) components of diffusion motion are well separated. The rotational correlation time shows an Arrhenius-type behavior with an activation energy of 4.48 kJ/mol, which is lower compared to bulk water. The rotational diffusion at room temperature is slower by about a factor of 2 compared to bulk water, whereas the translational diffusion slows down by a factor of 40. In contrast to bulk water, the translational correlation time exhibits an Arrhenius-type temperature dependence with an activation energy of 11.38 kJ/mol. Comparison of different models for jump diffusion processes suggests that water molecules perform two-dimensional jumps at a well-defined, almost temperature-independent distance of 4.21-4.32 A. Such a large jump distance indicates a low molecular density of the layer of diffusing molecules. We argue that undissociated water molecules on an average form two hydrations layers on top of the surface layer of hydroxyl groups, and all the layers have similar molecular density. Quasielastic neutron scattering experiment assesses the dynamics of the outermost hydration layer, whereas slower motion of the water molecules in the inner hydration layer contributes to the elastic signal.     

Hydrogen-bonded network of hydration water around model solutes

Clustering of water molecules in the hydration shells of spherical structureless solutes was studied in dependence on thermodynamic state, solute radius R(sp) and strength U(0) of water-solute interaction. Two qualitatively different clustering states of hydration water have been found: an "ordered" state with a hydrogen-bonded (H-bonded) network, which includes most of the hydration water, and a "disordered" state with small H-bonded clusters of hydration water. The transition from the ordered to disordered state occurs upon increasing temperature and decreasing pressure. This percolation transition is rounded due to the finite solute size and occurs in some temperature (pressure) interval. A finite-size scaling was applied to determine the transition temperature T(∞) in the limit R(sp)→∞. Strengthening of the water-solute interaction strongly enhances the stability of the ordered state: the transition temperature increases by about 35 °C, when U(0) decreases by 1 kcal mol(-1). At T > T(∞) and fixed U(0), the stability of the H-bonded water network increases upon decreasing solute size.     

Percolation transition in supercritical water: a Monte Carlo simulation study

Computer simulations of water have been performed on the canonical ensemble at 15 different molecular number densities, ranging from 0.006 to 0.018 A-3, along the supercritical isotherm of 700 K, in order to characterize the percolation transition in the system. It is found that the percolation transition occurs at a somewhat higher density than what is corresponding to the supercritical extension of the boiling line. We have shown that the fractal dimension of the largest cluster and the probability of finding a spanning cluster are the most appropriate properties for the location of the true percolation threshold. Thus, percolation transition occurs when the fractal dimension of the largest cluster reaches 2.53, and the probability of finding a cluster that spans the system in at least one dimension and in all the three dimensions reaches 0.97 and 0.65, respectively. On the other hand, the percolation threshold cannot be accurately located through the cluster size distribution, as it is distorted by appearance of clusters crossing the finite simulated system even far below the percolation threshold. The structure of the largest water cluster is dominated by a linear, chainlike arrangement, which does not change noticeably until the largest cluster becomes infinite.     

B-DNA Helix Stability in a Solvent-Free Environment

B-DNA is the most common DNA helix conformation under physiological conditions. However, when the amount of water in a DNA solution is decreased, B-to-A helix transitions have been observed. To understand what type of helix conformations exist in a solvent-free environment, a series of poly d(CG)(n) and mixed sequence DNA duplexes from 18 to 30 bp were examined with circular dichroism (CD), ESI-MS, ion mobility, and molecular dynamics. From the CD spectra, it was observed that all sequences had B-form helices in solution. However, the solvent-free results were more complex. For the poly d(CG)(n) series, the 18 bp duplex had an A-form helix conformation, both A- and B-helices were present for the 22 bp duplex, and only B-helices were observed for the 26 and 30 bp duplexes. Since these sequences were all present as B-DNA in solution, the observed solvent-free structures illustrate that smaller helices with fewer base pairs convert to A-DNA more easily than larger helices in the absence of solvent. A similar trend was observed for the mixed sequence duplexes where both an A- and B-helix were present for the 18 bp duplex, while only B-helices occur for the larger 22, 26, and 30 bp duplexes. Since the solvent-free B-helices appear at smaller sizes for the mixed sequences than for the pure d(CG)(n) duplexes, the pure d(CG)(n) duplexes have a greater A-philicity.     

Changes in hydration of lanthanide ions on binding to DNA in aqueous solution

The interaction of the trivalent lanthanides Ce(III), Eu(III), and Tb(III) with sodium deoxyribonucleic acid (DNA) in aqueous solution has been studied using their luminescence spectra and decays. Complexation with DNA is indicated by changes in luminescence intensity. In the system terbium(III)-DNA, changes in luminescence with pH are suggested to be due to the protonation of phosphate groups. The degree of hydration of Tb(III) on binding to DNA is followed by luminescence lifetime measurements in water and deuterium oxide solutions, and it is found that the lanthanide ion loses at least one hydration water on binding to long double stranded DNA at pH 4.7 and pH 7. Rather different behavior is observed on binding to long or short single stranded DNA, where six water molecules are lost, independent of pH. It is suggested that in this case the lanthanide probably binds to the bases of the DNA backbone. The DNA conformation seems to be an important factor in the binding. In addition, the isotopic effect on terbium luminescence lifetime may provide a useful method to distinguish between single and double stranded DNA. DSC results are consistent with cleavage of the double helix of DNA at pH 9 in the presence of terbium.     

Hydration and temperature interdependence of protein picosecond dynamics

We investigate the nature of the solvent motions giving rise to the rapid temperature dependence of protein picoseconds motions at 220 K, often referred to as the protein dynamical transition. The interdependence of picoseconds dynamics on hydration and temperature is examined using terahertz time domain spectroscopy to measure the complex permittivity in the 0.2-2.0 THz range for myoglobin. Both the real and imaginary parts of the permittivity over the frequency range measured have a strong temperature dependence at >0.27 h (g water per g protein), however the permittivity change is strongest for frequencies <1 THz. The temperature dependence of the real part of the permittivity is not consistent with the relaxational response of the bound water, and may reflect the low frequency protein structural vibrations slaved to the solvent excitations. The hydration necessary to observe the dynamical transition is found to be frequency dependent, with a critical hydration of 0.19 h for frequencies >1 THz, and 0.27 h for frequencies <1 THz. The data are consistent with the dynamical transition solvent fluctuations requiring only clusters of ~5 water molecules, whereas the enhancement of lowest frequency motions requires a fully spanning water network.     

Molecular dynamics simulation of the A-DNA to B-DNA transition in aqueous RbCl solution

Unrestrained molecular dynamics (MD) simulations have been carried out to characterize the stability of DNA conformations and the dynamics of A-DNA→B-DNA conformational transitions in aqueous RbCl solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbCl concentration on the dynamics of the A→B conformational transition in the DNA duplex d(CGCGAATTCGCG)2 . Canonical Aand B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was monitored for 3.0 ns and the distances between the C5′ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simulations. In addition, increasing the RbCl concentration in aqueous solution hindered the A→B conformational transition and the transition in aqueous RbCl solution was faster than that in aqueous NaCl solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A→B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The rubidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help explain the effect of RbCl concentration on the dynamics of the A→B conformational transition.     

MD simulation of the transitions between B-DNA and A-DNA in the framework of a coarse-grained model

Transitions between the B and A forms of a short DNA double helix (12 base pairs) at different salt concentrations in an aqueous solution have been studied by the molecular dynamics method in the framework of a coarse-grained model with explicit ions but without friction. It has been shown that the A-DNA, stable at high salt concentrations, is a dynamic conglomerate of the molecule and the ions coming from the solution into the deep major groove and then leaving it. In such a short helix, in the model without friction, even at low salt concentrations, transitions from B-DNA to A-DNA and back are frequent and fast. Stable ADNA (without transitions to B-DNA) forms at salt concentrations greater than 0.45 mol/L.     

Ultrafast Vibrational Dynamics and Local Interactions of Hydrated DNA

Biochemical processes occur mainly in aqueous environments, where interactions with water molecules play a key role for both the structure and function of biomolecules. Deoxyribonucleic acid (DNA), the basic carrier of genetic information, is characterized by an equilibrium double helix structure which is held together by intermolecular hydrogen bonds between base pairs and hydrated by an environment of water molecules with fluctuating hydrogen bonds. Basic vibrational motions of hydrated DNA and the fastest changes in the DNA–water interactions and hydration geometries occur in less than 1 ps. These processes can be accessed by mapping the vibrational dynamics of DNA and water in a time‐resolved way by nonlinear ultrafast vibrational spectroscopy. Recent studies provide a detailed understanding of DNA vibrations and their dynamics, and give insight into nonequilibrium properties and structures of hydrated DNA.