Biomedical subjects

A P Lyubartsev

Publications and source records attributed to A P Lyubartsev.

7 recordsLinked to original sources

PubMed2001

Spermine: an "invisible" component in the crystals of B-DNA. A grand canonical Monte Carlo and molecular dynamics simulation study.

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.

Algorithms↗

PubMed2001

Competitive substitution of hexammine cobalt(III) for Na+ and K+ ions in oriented DNA fibers.

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

Cobalt↗

PubMed1999

Competitive binding of Mg2+, Ca2+, Na+, and K+ ions to DNA in oriented DNA fibers: experimental and Monte Carlo simulation results.

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(+)).

Animals↗

PubMed1999

Topological and spatial aspects of the hydration of solutes of extreme solvation entropy.

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.

Journal Article↗

PubMed1998

Application of polyelectrolyte theories for analysis of DNA melting in the presence of Na+ and Mg2+ ions.

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.

Animals↗

PubMed1998

Molecular dynamics simulations of DNA in solutions with different counter-ions.

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.

Cesium↗

PubMed1989

Monte-Carlo-self consistent field method in the polyelectrolyte theory.

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.

Algorithms↗