RESUMEN

The relation between the accuracy of centroid molecular dynamics correlation functions, and the geometry of the centroid potential is investigated. It is shown that, depending on the temperature, there exist several regimes, and in each of them certain features of the exact Kubo correlation functions are reproduced. The change of regimes is related to the emergence of barriers in the centroid potential. In order to clarify how the above described picture of regimes is modified in real systems when dissipation is important, a methodology is developed to test the accuracy of centroid correlation functions for the model of a particle coupled to a harmonic heat bath. A modification of the centroid molecular dynamics method to include the influence of the heat bath is introduced. Preliminary results of comparison of centroid molecular dynamics with the numerically exact results of filtered propagator functional method are presented.

RESUMEN

The behaviour of a flexible anionic chain of 150 univalent and negatively charged beads connected by a harmonic-like potential with each other in the presence of an equal number of positive and free counterions, is studied in molecular dynamics simulations with Langevin thermostat in a wide range of temperatures. Simulations were carried out for several values of the bending parameter, corresponding to fully flexible polyion, moderately and strongly stiff polyion as well as for the case when bend conformation is preferable to the straight one. We have found that in all cases three regimes can be distinguished, which can be characterized as "random coil", observed at high temperatures; "extended conformation" observed at moderate temperatures (of the order of 1 in reduced units), and compact "globular conformation" attained at low temperatures. While the transition between high-temperature random and extended conformations is gradual, the transition from the extended coil to the globular state, taking place at a temperature of about 0.2 in reduced units, is of abrupt character resembling a phase transition.

RESUMEN

Small-angle X-ray scattering (SAXS) is used to demonstrate the presence of density fluctuations in ambient water on a physical length-scale of approximately 1 nm; this is retained with decreasing temperature while the magnitude is enhanced. In contrast, the magnitude of fluctuations in a normal liquid, such as CCl(4), exhibits no enhancement with decreasing temperature, as is also the case for water from molecular dynamics simulations under ambient conditions. Based on X-ray emission spectroscopy and X-ray Raman scattering data we propose that the density difference contrast in SAXS is due to fluctuations between tetrahedral-like and hydrogen-bond distorted structures related to, respectively, low and high density water. We combine our experimental observations to propose a model of water as a temperature-dependent, fluctuating equilibrium between the two types of local structures driven by incommensurate requirements for minimizing enthalpy (strong near-tetrahedral hydrogen-bonds) and maximizing entropy (nondirectional H-bonds and disorder). The present results provide experimental evidence that the extreme differences anticipated in the hydrogen-bonding environment in the deeply supercooled regime surprisingly remain in bulk water even at conditions ranging from ambient up to close to the boiling point.

Asunto(s)

RESUMEN

We analyze the validity of the commonly used electric-field (E-field) approximation to vibrational OH stretch Raman spectra of dilute HOD in D(2)O by computing the OH stretch frequency of all molecules in several different structure models, each containing around 2000 molecules. The calculations are done at the B3LYP level using clusters containing 32 molecules centered around the molecule for which the frequencies are calculated; the large cluster size is required due to significant nonlocal contributions influencing the computed frequencies. The vibrational frequencies are determined using a six-point potential optimized discrete variable representation. Raman and infrared intensities are furthermore computed to generate the spectra. We find that a quadratic fit of E-field versus frequency gives a reasonable representation of the calculated distribution of frequencies. However, the mapping depends significantly on the structural model and is thus not universal. Anharmonic couplings are calculated for several optimized clusters showing a general trend to compress the computed frequency distributions, which is in agreement with dynamical simulations (motional narrowing).

Asunto(s)

RESUMEN

Expanded-ensemble Monte Carlo method with Wang-Landau algorithm was used for calculations of the ratio of partition functions for classes of permutations in the problem of several interacting quantum particles (fermions) in an external field. Simulations for systems consisting of 2 up to 7 interacting particles in harmonic or Coulombic field were performed. The presented approach allows one to carry out calculations for low enough temperatures that makes it possible to extract data for the ground-state energy and low-temperature thermodynamics.

Asunto(s)

RESUMEN

We extend the Monte Carlo methods developed in our previous papers [J. Phys. A 37, 1573 (2004); Macromol. Theory-Simul. 14, 491 (2005)] and based on entropic sampling within the Wang-Landau algorithm to simulation of a lattice model of flexible polyelectrolytes. We consider a strongly charged polyelectrolyte chain accompanied by neutralizing counterions on a simple cubic lattice with periodic boundary conditions. The Coulomb potential and the excluded volume condition between different ions or beads are taken into account. The obtained energy distributions make possible the calculation of canonical properties such as conformational energy, heat capacity, entropy, free energy, and mean-square end-to-end distance over a wide temperature range in a single simulation.

RESUMEN

The association of spermine(4+) (Spm(4+)), Mg(2+) and monovalent (M(+)) ions with DNA in crystal form, have been studied using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) computer simulations. GCMC calculations were used to calculate the distribution of Spm(4+), Mg(2+), and M(+) between the equilibrating solvent and the DNA crystal under conditions mimicking the crystal-growing protocols reported in a number of recent X-ray diffraction studies of DNA oligomers. The GCMC simulations show that the composition of ions neutralizing the negative charge of DNA can vary in a broad range. The GCMC simulations were used to provide appropriate conditions for subsequent 6 ns constant pressure and temperature MD simulations of DNA in a typical crystalline environment consisting of three DNA double helix decamers in a periodic hexagonal cell, containing 1200 water molecules, eight Spm(4+), 32 Na(+) and four Cl(-) ions. Based on the simulation results, it seems possible to give an explanation why spermine molecules are usually not detected in X-ray studies in spite of their high concentration in the preparatory samples used as the crystallizing agent. It appears that this flexible polyamine molecule has several binding modes, interacting in fairly irregular manner with different sites on DNA and showing no regular ordering in the DNA crystals. Ions of Na(+) and Spm(4+) compete with each other and with water molecules in binding to bases in the minor groove and they influence the structure of the DNA hydration shell in different ways.

Asunto(s)

RESUMEN

Competition of the trivalent cation, Co(NH3)(3+)(6), with K+ and Na+ ions in binding to DNA was studied by equilibrating oriented DNA fibers with ethanol/water solutions (65 and 52% v/v EtOH), containing different combinations and concentrations of KCl and NaCl and constant concentration (0.8 mM) of Co(NH3)(6)Cl(3). The degree of Co(NH3)(3+)(6) binding to DNA does not depend significantly on the ethanol concentration or on the kind of univalent cation (Na+ or K+). The ion exchange selectivity coefficient of monovalent-trivalent ion competition, D(1)(c3), increases with the concentration of Me+, C(o)(+), and the monotonic dependence of log D(1)(c3) vs log C(o)(+) has an inflection between 100 and 300 mM that is caused by a structural transformation of DNA from A- to B-form. The ion exchange experimental data are compared with results of grand canonical Monte Carlo (GCMC) simulations of systems of parallel and hexagonally ordered, discretely charged polyions with density and spatial distribution of the charged groups modeling B- and A-forms of DNA. The GCMC method for discretely charged models of the DNA polyion produces a quantitative agreement with experimental data on trivalent-monovalent ion competition in dependence on DNA structural state and salt concentration. Based on this and previous studies it is concluded that the affinity of DNA for the cations decreases in the order Co(NH3)(3+)(6) >> Ca2+ > Mg2+ >> Na+ approximately K+ > Li+. DNA does not exhibit selectivity for Na+ or K+ in ethanol/water solutions either in the absence or in the presence of Co(NH3)(3+)(6), Ca2+, and Mg2+.

Asunto(s)

RESUMEN

Competitive binding of the most common cations of the cytoplasm (K(+), Na(+), Ca(2+), and Mg(2+)) with DNA was studied by equilibrating oriented DNA fibers with ethanol/water solutions (65 and 52% v/v EtOH) containing different combinations and concentrations of the counterions. The affinity of DNA for the cations decreases in the order Ca > Mg >> Na approximately K. The degree of Ca(2+) and/or Mg(2+) binding to DNA displays maximum changes just at physiological concentrations of salts (60-200 mM) and does not depend significantly on the ethanol concentration or on the kind of univalent cation (Na(+) or K(+)). Ca(2+) is more tightly bound to DNA and is replaced by the monovalent cations to a lesser extent than is Mg(2+). Similarly, Ca(2+) is a better competitor for binding to DNA than Mg(2+): the ion exchange equilibrium constant for a 1:1 mixture of Ca(2+) and Mg(2+) ions, K(c)(Ca)(Mg), changes from K(c)(Ca)(Mg) approximately 2 in 65% EtOH (in 3-30 mM NaCl and/or KCl) to K(c)(Ca)(Mg) approximately 1.2-1.4 in 52% EtOH (in 300 mM NaCl and/or KCl). DNA does not exhibit selectivity for Na(+) or K(+) in ethanol/water solutions either in the absence or in the presence of Ca(2+) and/or Mg(2+). The ion exchange experimental data are compared with results of grand canonical Monte Carlo (GCMC) simulations of systems of parallel and hexagonally ordered, uniformly and discretely charged polyions with the density and spatial distribution of the charged groups modeling B DNA. A quantitative agreement with experimental data on divalent-monovalent competition has been obtained for discretely charged models of the DNA polyion (for the uniformly charged cylinder model, coincidence with experiment is qualitative). The GCMC method gives also a qualitative description of experimental results for DNA binding competitions of counterions of the same charge (Ca(2+) with Mg(2+) or K(+) with Na(+)).

Asunto(s)

RESUMEN

The hydration of charged Lennard-Jones spheres by simple point charge water is considered. Molecular dynamics and expanded ensemble simulations were used to compare the hydration structures surrounding solutes with extreme solvation entropy. The variations in the solvation entropy were analyzed in terms of changes in the spatial and topological structure of the hydration shells. The solvation entropy was found to be maximal for solutes that can replace water molecules in the hydrogen-bond network. Further, using a Kirkwood-type factorization, the solvation entropy was expanded as a sum over the partial n-body distribution functions. The two-body solute-water contribution to the solvation entropy was found to exceed the full solvation entropy for solutes with low charge, whereas the converse is true for the other solutes. This is consistent with the idea that water-water correlations are enhanced by solvation of, for example, noble gases, whereas they are disrupted by solvation of ions. Further, the orientational and radial parts of the two-body solute-water entropy were calculated as functions of the charge of the solute. The orientational part has a single maximum, whereas the radial part maintains the bimodal form of the full solvation entropy.

RESUMEN

Numerical calculations, using Poisson-Boltzmann (PB) and counterion condensation (CC) polyelectrolyte theories, of the electrostatic free energy difference, DeltaGel, between single-stranded (coil) and double-helical DNA have been performed for solutions of NaDNA + NaCl with and without added MgCl2. Calculations have been made for conditions relevant to systems where experimental values of helix coil transition temperature (Tm) and other thermodynamic quantities have been measured. Comparison with experimental data has been possible by invoking values of Tm for solutions containing NaCl salt only. Resulting theoretical values of enthalpy, entropy, and heat capacity (for NaCl salt-containing solutions) and of Tm as a function of NaCl concentration in NaCl + MgCl2 solutions have thus been obtained. Qualitative and, to a large extent, quantitative reproduction of the experimental Tm, DeltaHm, DeltaSm, and DeltaCp values have been found from the results of polyelectrolyte theories. However, the quantitative resemblance of experimental data is considerably better for PB theory as compared to the CC model. Furthermore, some rather implausible qualitative conclusions are obtained within the CC results for DNA melting in NaCl + MgCl2 solutions. Our results argue in favor of the Poisson-Boltzmann theory, as compared to the counterion condensation theory.

Asunto(s)

RESUMEN

Molecular dynamics simulations of the [d(ATGCAGTCAG]2 fragment of DNA, in water and in the presence of three different counter-ions (Li+, Na+ and Cs+) are reported. Three-dimensional hydration structure and ion distribution have been calculated using spatial distribution functions for a detailed picture of local concentrations of ions and water molecules around DNA. According to the simulations, Cs+ ions bind directly to the bases in the minor groove, Na+ ions bind prevailing to the bases in the minor groove through one water molecule, whereas Li+ ions bind directly to the phosphate oxygens. The different behavior of the counter-ions is explained by specific hydration structures around the DNA and the ions. It is proposed how the observed differences in the ion binding to DNA may explain different conformational behavior of DNA. Calculated self-diffusion coefficients for the ions agree well with the available NMR data.

Asunto(s)

RESUMEN

A new time saving numerical method for calculation of equilibrium potential and density distribution of mobile ions around the polyion in a polyelectrolyte system is proposed: the region around the polyion is being divided into two zones-internal and external; in the internal zone all the ions are accounted explicity with the aid of Monte-Carlo procedure; in the external zone the combined Monte-Carlo-self consistent field method proposed earlier is applied, an exchange of ions between regions is being implied. For 1:1 electrolyte the optimal choice of the boundary between the zones has been demonstrated. As an example of a more complicated system calculation for 2:2:1:1 electrolyte was carried out.

Asunto(s)