Biomedical subjects

Alexander P Lyubartsev

Publications and source records attributed to Alexander P Lyubartsev.

16 recordsLinked to original sources

PubMed2006

Dynamical and structural properties of charged and uncharged lidocaine in a lipid bilayer.

Molecular dynamics computer simulations have been performed to investigate dynamical and structural properties of a lidocaine local anesthetic. Both charged and uncharged forms of the lidocaine molecule were investigated. Properties such as membrane area per lipid, diffusion, mass density, bilayer penetration and order parameters have been examined. An analysis of the lidocaine interaction with the lipid surrounding according to a simple mean field theory has also been performed. Almost all examined properties were found to depend on which of the two forms of lidocaine, charged or uncharged, is studied. The overall picture is a rather static behavior determined by the lipids for the charged molecules and more mobile situation of the uncharged form with higher diffusion and lower orientational and positional order.

Anesthetics, Local↗

PubMed2006

A molecular dynamics investigation of the influence of hydration and temperature on structural and dynamical properties of a dimyristoylphosphatidylcholine bilayer.

A series of molecular dynamics simulations of dimyristoylphosphatidylcholine bilayers, with different levels of hydration and temperature, were performed to examine the influence of hydration on properties of lipid membranes. Structural and dynamical properties such as area per lipid, electron densities, order parameters for all CH bonds and water, diffusion, and reorientation autocorrelation functions were determined and were all found to be affected by changes in the hydration level. The simulations give an overall picture of the bilayer going to a more ordered state when the hydration level is reduced. Lipid headgroups were found to adopt an orientation more parallel to the membrane plane when the water content was decreased. Dynamical properties such as lipid diffusion and relaxation of reorientation time correlation functions were found to become slower with the removal of water. Our simulation results generally agree with experimental data in cases where such data are available. One important conclusion drawn is that while structural properties are affected only moderately dynamical properties are affected very strongly by a decrease of water content.

Computer Simulation↗

PubMed2006

Computer modeling demonstrates that electrostatic attraction of nucleosomal DNA is mediated by histone tails.

We conducted molecular dynamics computer simulations of charged histone tail-DNA interactions in systems mimicking nucleosome core particles (NCP) . In a coarse-grained model, the NCP is modeled as a negatively charged spherical particle with flexible polycationic histone tails attached to it in a dielectric continuum with explicit mobile counterions and added salt. The size, charge, and distribution of the tails relative to the core were built to mimick real NCP. In this way, we incorporate attractive ion-ion correlation effects due to fluctuations in the ion cloud and the attractive entropic and energetic tail-bridging effects. In agreement with experimental data, increase of monovalent salt content from salt-free to physiological concentration leads to the formation of NCP aggregates; likewise, in the presence of MgCl2, the NCPs form condensed systems via histone-tail bridging and accumulation of counterions. More detailed mechanisms of the histone tail-DNA interactions and dynamics have been obtained from all-atom molecular dynamics simulations (including water), comprising three DNA 22-mers and 14 short fragments of the H4 histone tail (amino acids 5-12) carrying three positive charges on lysine+ interacting with DNA. We found correlation of the DNA-DNA distance with the presence and association of the histone tail between the DNA molecules.

Binding Sites↗

PubMed2005

Multiscale modeling of lipids and lipid bilayers.

A multiscale modeling approach is applied for simulations of lipids and lipid assemblies on mesoscale. First, molecular dynamics simulation of initially disordered system of lipid molecules in water within all-atomic model was carried out. On the next stage, structural data obtained from the molecular dynamics (MD) simulation were used to build a coarse-grained (ten sites) lipid model, with effective interaction potentials computed by the inverse Monte Carlo method. Finally, several simulations of the coarse-grained model on longer length- and time-scale were performed, both within Monte Carlo and molecular dynamics simulations: a periodical sample of lipid molecules ordered in bilayer, a free sheet of such bilayer without periodic boundary conditions, formation of vesicle from a plain membrane, process of self-assembly of lipids randomly dispersed in volume. It was shown that the coarse-grained model, developed exclusively from all-atomic simulation data, reproduces well all the basic features of lipids in water solution.

Computer Simulation↗

PubMed2005

Simulations of one- and two-electron systems by Bead-Fourier path integral molecular dynamics.

The Bead-Fourier path integral molecular dynamics technique introduced earlier [S. D. Ivanov, A. P. Lyubartsev, and A. Laaksonen, Phys. Rev. E 67 066710 (2003)] is applied for simulation of electrons in the simplest molecules: molecular hydrogen, helium atom, and their ions. Special attention is paid to the correct description of electrons in the core region of a nucleus. In an attempt to smooth the Coulomb potential at small distances, a recipe is suggested. The simulation results are in excellent agreement with the analytical solution for the "harmonic helium atom", as well as with the vibrational potential of the H(2) molecule and He ionization energies. It is demonstrated, that the Bead-Fourier path integral molecular dynamics technique is able to provide the accuracy required for the description of electron structure and chemical bonds in cases when electron exchange effects need not be taken into account.

Journal Article↗

PubMed2005

Solvation structure of hydroxyl radical by Car-Parrinello molecular dynamics.

Car-Parrinello molecular dynamics simulations of a hydroxyl radical in liquid water have been performed. Structural and dynamical properties of the solvated structure have been studied in details. The partial atom-atom radial distribution functions for the hydrated hydroxyl do not show drastic differences with the radial distribution functions for liquid water. The OH is found to be a more active hydrogen bond donor and acceptor than the water molecule, but the accepted hydrogen bonds are much weaker than for the hydroxide OH- ion. The first solvation shell of the OH is less structured than the water's one and contains a considerable fraction of water molecules that are not hydrogen bonded to the hydroxyl. Part of them are found to come closer to the solvated radical than the hydrogen bonded molecules do. The lifetime of the hydrogen bonds accepted by the hydroxyl is found to be shorter than the hydrogen bond lifetime in water. A hydrogen transfer between a water molecule and the OH radical has been observed, though it is a much rarer event than a proton transfer between water and an OH- ion. The velocity autocorrelation power spectrum of the hydroxyl hydrogen shows the properties both of the OH radical in clusters and of the OH- ion in liquid.

Computer Simulation↗

PubMed2004

A molecular dynamics simulation study of polyamine- and sodium-DNA. Interplay between polyamine binding and DNA structure.

Four different molecular dynamics (MD) simulations have been performed for infinitely long ordered DNA molecules with different counterions, namely the two natural polyamines spermidine(3+) (Spd3+) and putrescine(2+) (Put2+), the synthetic polyamine diaminopropane(2+) (DAP2+), and the simple monovalent cation Na+. All systems comprised a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charge. The simulation setup mimics the DNA state in oriented DNA fibers, previously studied using NMR and other experimental methods. In this paper the interplay between polyamine binding and local DNA structure is analyzed by investigating how and if the minor groove width of DNA depends on the presence and dynamics of the counterions. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural [spermine(4+), Spd3+, Put2+] and the synthetic (DAP2+) polyamines.

Base Sequence↗

PubMed2004

Molecular dynamics simulation study of oriented polyamine- and Na-DNA: sequence specific interactions and effects on DNA structure.

Molecular dynamics (MD) computer simulations have been carried out on four systems that correspond to an infinite array of parallel ordered B-DNA, mimicking the state in oriented DNA fibers and also being relevant for crystals of B-DNA oligonucleotides. The systems were all comprised of a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charges. The sequence of the double helical DNA decamer was d(5'-ATGCAGTCAG)xd(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine(3+) (Spd(3+)) and putrescine(2+) (Put(2+)), the synthetic polyamine diaminopropane(2+) (DAP(2+)), and the simple monovalent cation Na(+). This work compares the specific structures of the polyamine- and Na-DNA systems and how they are affected by counterion interactions. It also describes sequence-specific hydration and interaction of the cations with DNA. The local DNA structure is dependent on the nature of the counterion. Even the very similar polyamines, Put(2+) and DAP(2+), show clear differences in binding to DNA and in effect on hydration and local structure. Generally, the polyamines disorder the hydration of the DNA around their binding sites whereas Na(+) being bound to DNA attracts and organizes water in its vicinity. Cation binding at the selected sites in the minor and in the major groove is compared for the different polyamines and Na(+). We conclude that the synthetic polyamine (DAP(2+)) binds specifically to several structural and sequence-specific motifs on B-DNA, unlike the natural polyamines, Spd(3+) and Put(2+). This specificity of DAP(2+) compared to the more dynamic behavior of Spd(3+) and Put(2+) may explain why the latter polyamines are naturally occurring in cells.

Base Sequence↗

PubMed2004

Determination of solvation free energies by adaptive expanded ensemble molecular dynamics.

A new method of calculating absolute free energies is presented. It was developed as an extension to the expanded ensemble molecular dynamics scheme and uses probability density estimation to continuously optimize the expanded ensemble parameters. The new method is much faster as it removes the time-consuming and expertise-requiring step of determining balancing factors. Its efficiency and accuracy are demonstrated for the dissolution of three qualitatively very different chemical species in water: methane, ionic salts, and benzylamine. A recently suggested optimization scheme by Wang and Landau [Phys. Rev. Lett. 86, 2050 (2001)] was also implemented and found to be computationally less efficient than the proposed adaptive expanded ensemble method.

Journal Article↗

PubMed2004

Electrostatic background of chromatin fiber stretching.

We have carried out an investigation of the electrostatic forces involved in gradual removal of the DNA from the histone proteins in chromatin. Two simple models of DNA-histone core dissociation were considered. Calculations of the electrostatic free energy within the Poisson-Boltzmann theory gave similar results for the both models, which turned out to be in a qualitative agreement with recent optical tweezers stretching experiments measuring the force necessary to unwrap DNA from the histone core. Our analysis shows that the electrostatic interactions between the highly negatively charged polymeric DNA and the positively charged histones play a determining role in stabilizing the nucleosomes at physiological conditions.

Chromatin↗

PubMed2003

A molecular dynamics simulation study of oriented DNA with polyamine and sodium counterions: diffusion and averaged binding of water and cations.

Four different molecular dynamics (MD) simulations have been performed for ordered DNA decamers, d(5'-ATGCAGTCAG)*d(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine3+ (Spd3+) and putrescine2+ (Put2+), the synthetic polyamine diaminopropane2+ (DAP2+) and Na+. The simulation set-up corresponds to an infinite array of parallel DNA mimicking the state in oriented DNA fibers or crystals. This work describes general properties of polyamine and Na+ binding to DNA. Simulated diffusion coefficients show satisfactory agreement with experimental NMR diffusion data of comparable systems. The interaction of the polyamines with DNA is dynamic in character and the cations mostly form short-lived contacts with the electronegative binding sites of DNA. Polyamines, Na+ and water interact most frequently with the charged phosphate atoms with preference for association from the minor groove side with O1P over O2P. There is a strong anti-correlation in the cation binding to the electronegative groups of DNA, i.e. the presence of a cation near one of the DNA sites repels other cations from binding to this and to the other sites separated by <7.5 A from each other. In contrast to the other polyamines, DAP2+ is able to form 'bridges' connecting neighboring phosphate groups along the DNA strand. A small fraction of DAP2+ and Put2+ can be found in the major grooves, while Spd3+ is absent there. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural (Spd3+, Put2+ and spermine4+) and synthetic (DAP2+) polyamines.

Binding Sites↗

PubMed2003

Bead-Fourier path integral molecular dynamics.

Molecular dynamics formulation of Bead-Fourier path integral method for simulation of quantum systems at finite temperatures is presented. Within this scheme, both the bead coordinates and Fourier coefficients, defining the path representing the quantum particle, are treated as generalized coordinates with corresponding generalized momenta and masses. Introduction of the Fourier harmonics together with the center-of-mass thermostating scheme is shown to remove the ergodicity problem, known to pose serious difficulties in standard path integral molecular dynamics simulations. The method is tested for quantum harmonic oscillator and hydrogen atom (Coulombic potential). The simulation results are compared with the exact analytical solutions available for both these systems. Convergence of the results with respect to the number of beads and Fourier harmonics is analyzed. It was shown that addition of a few Fourier harmonics already improves the simulation results substantially, even for a relatively small number of beads. The proposed Bead-Fourier path integral molecular dynamics is a reliable and efficient alternative to simulations of quantum systems.

Journal Article↗

PubMed2003

Application of the Poisson Boltzmann polyelectrolyte model for analysis of thermal denaturation of DNA in the presence of Na+ and polyamine cations.

The Poisson Boltzmann (PB) cell model of polyelectrolyte solution has been used for calculation of the electrostatic free energy difference, Delta G(el), between double- and single-stranded DNA. The calculations have been performed for conditions relevant to describe the DNA helix-coil transition in NaCl solution in the presence of the natural polyamines putrescine(2+), spermidine(3+), spermine(4+) and their synthetic homologs with different spacing between the charged amino groups, for which experimental values of the DNA 'melting' transition temperature (T(m)) are available. Using the PB theory and the polyamine ion radius as an adjusting parameter provides quantitative agreement between experimental and theoretical T(m)--salt concentration dependencies only by using physically unreasonable radii for the polyamine. Thus, modeling the linear and flexible polyamines as charged spheres within the PB cell model is an implausible oversimplification. We propose another explanation for the experimental observations, still within the frame of the 'primitive' PB polyelectrolyte theory. This explanation is based on an analysis of the Delta G(el) dependence on the stoichiometry of polyamine-polyanion binding to double- and single-stranded DNA.

Cations↗

PubMed2002

Evaluation of effective ion-ion potentials in aqueous electrolytes.

We present high-accuracy simulation data for ion-ion radial distribution functions in NaCl aqueous electrolyte. From the data we evaluate ion-ion effective potentials via inverse Monte Carlo simulation procedure. Alternatively, we first evaluate effective direct correlation functions and then obtain effective potentials using the hypernetted-chain (HNC) approximation. Both methods are in excellent agreement, indicating that HNC approximation performs well with effective ion-ion potentials in aqueous solutions.

Journal Article↗

PubMed2002

On the competition between water, sodium ions, and spermine in binding to DNA: a molecular dynamics computer simulation study.

The interaction of DNA with the polyamine spermine(4+) (Spm(4+)), sodium ions, and water molecules has been studied using molecular dynamics computer simulations in a system modeling a DNA crystal. The simulation model consisted of three B-DNA decamers in a periodic hexagonal cell, containing 1200 water molecules, 8 Spm(4+), 32 Na(+), and 4 Cl(-) ions. The present paper gives a more detailed account of a recently published report of this system and compares results on this mixed Spm(4+)/Na(+)-cation system with an molecular dynamics simulation carried out for the same DNA decamer under similar conditions with only sodium counterions (Korolev et al., J. Mol. Biol. 308:907). The presence of Spm(4+) makes significant influence on the DNA hydration and on the interaction of the sodium ions with DNA. Spermine pushes water molecules out of the minor groove, whereas Na(+) attracts and organizes water around DNA. The major binding site of the Spm(4+) amino groups and the Na(+) ions is the phosphate group of DNA. The flexible polyamine spermine displays a high presence in the minor groove but does not form long-lived and structurally defined complexes. Sodium ions compete with Spm(4+) for binding to the DNA bases in the minor groove. Sodium ions also have several strong binding sites in the major groove. The ability of water molecules, Spm(4+), and Na(+) to modulate the local structure of the DNA double helix is discussed.

Base Sequence↗

PubMed2002

Application of the Poisson Boltzmann polyelectrolyte model for analysis of equilibria between single-, double-, and triple-stranded polynucleotides in the presence of K(+), Na(+), and Mg(2+) ions.

The Poisson Boltzmann (PB) cell model of polyelectrolyte solution has been used for numerical calculations of the change in electrostatic free energy, DeltaG(el), for transformations between different structural forms (tri-, double-, and single-stranded) of the polyribonucleotides poly(rA).poly(rU)2, poly(rA).poly(rU), poly(rA), and poly(rU). In particular, the dependence on monovalent salt concentration, MCl (M = Na or K) in the absence and in the presence of MgCl(2) has been calculated. The results were obtained for conditions relevant to available experimental values of structural transition ("melting") temperatures (T(m)) and other thermodynamic quantities. Using the experimental T(m) values and theoretical electrostatic DeltaG(el), DeltaH(el), and DeltaS(el) functions, non-electrostatic contributions to the corresponding thermodynamic parameters of the poly(rA)/poly(rU) melting transitions were determined in MCl solutions in the absence of Mg(2+). Qualitative and to a large extent quantitative reproduction of the experimental calorimetry enthalpy, entropy and heat capacity values was found from the results of the PB theory. Furthermore, dependencies of T(m) on MCl concentration in the presence of MgCl(2) for the poly(rA)/poly(rU) transitions were also calculated. Compared to a model considering Mg(2+) as a fully-hydrated ion, much better agreement between experimental and PB theory was achieved by assuming the ion size of Mg(2+) to be given by that of a bare non-hydrated ion smaller than that of hydrated Na(+) or K(+). In agreement with data of experimental studies reported in literature, this result indicates that magnesium(II) can bind to RNA as a bare ion in a way that is different from that of DNA. Generally, we can conclude that the PB polyelectrolyte theory can provide an adequate description of thermally induced structural transitions in polydeoxy- and polyribonucleotides in different salt solutions in spite of the rather simplified model treating the solvent as dielectric continuum, the polyion as a uniformly charged cylinder, and the mobile ions as hard spheres in the absence of excluded volume effects.

DNA↗