Lamellar and hexagonal lipid structures are of particular importance in the biological processes such as membrane fusion and budding. Atomistic simulations of formation of these phases and ...transitions between them are computationally prohibitive, hence development of coarse-grained models is an important part of the methodological development in this area. Here we apply systematic bottom-up coarse-graining to model different phase structures formed by 1,2-dioleoylphosphatidylethanolamine (DOPE) lipid molecules. We started from atomistic simulations of DOPE lipids in water carried out at two different water/lipid molar ratio corresponding to the lamellar Lα and inverted hexagonal HII structures at low and high lipid concentrations respectively. The atomistic trajectories were mapped to coarse-grained trajectories, in which each lipid was represented by 14 coarse-grained sites. Then the inverse Monte Carlo method was used to compute the effective coarse-grained potentials which for the coarse-grain model reproduce the same structural properties as the atomistic simulations. The potentials derived from the low concentration atomistic simulation were only able to form a bilayer structure, while both Lα and HII lipid phases were formed in simulations with potentials obtained at high concentration. The typical atomistic configurations of lipids at high concentration combine fragments of both lamellar and non-lamellar structures, that is reflected in the extracted coarse-grained potentials which become transferable and can form a wide range of structures including the inverted hexagonal, bilayer, tubule, vesicle and micellar structures.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
It is well established that the presence of the trivalent cobalt(III)-hexammine cation (CoHex3+) at submillimolar concentrations leads to bundling (condensation) of double-stranded DNA molecules, ...which is caused by DNA–DNA attraction induced by the multivalent counterions. However, the detailed mechanism of this process is still not fully understood. Furthermore, in all-atom molecular dynamics (MD) simulations, spontaneous aggregation of several DNA oligonucleotides in the presence of CoHex3+ has previously not been demonstrated. In order to obtain a rigorous description of CoHex3+–nucleic acid interactions and CoHex3+-induced DNA condensation to be used in MD simulations, we have derived optimized force field parameters of the CoHex3+ ion. They were obtained from Car–Parrinello molecular dynamics simulation of a single CoHex3+ ion in the presence of 125 water molecules. The new set of force field parameters reproduces the experimentally known transition of DNA from B- to A-form, and qualitatively describes changes of DNA and RNA persistence lengths. We then carried out a 2 μs long atomistic simulation of four DNA oligomers each consisting of 36 base pairs in the presence of CoHex3+. We demonstrate that, in this system, DNA molecules display attractive interactions and aggregate into bundle-like structures. This behavior depends critically on the details of the CoHex3+ interaction with DNA. A control simulation with a similar setup but in the presence of Mg2+ does not induce DNA–DNA attraction, which is also in agreement with experiment.
Computer modeling of very large biomolecular systems, such as long DNA polyelectrolytes or protein-DNA complex-like chromatin cannot reach all-atom resolution in a foreseeable future and this ...necessitates the development of coarse-grained (CG) approximations. DNA is both highly charged and mechanically rigid semi-flexible polymer and adequate DNA modeling requires a correct description of both its structural stiffness and salt-dependent electrostatic forces. Here, we present a novel CG model of DNA that approximates the DNA polymer as a chain of 5-bead units. Each unit represents two DNA base pairs with one central bead for bases and pentose moieties and four others for phosphate groups. Charges, intra- and inter-molecular force field potentials for the CG DNA model were calculated using the inverse Monte Carlo method from all atom molecular dynamic (MD) simulations of 22 bp DNA oligonucleotides. The CG model was tested by performing dielectric continuum Langevin MD simulations of a 200 bp double helix DNA in solutions of monovalent salt with explicit ions. Excellent agreement with experimental data was obtained for the dependence of the DNA persistent length on salt concentration in the range 0.1-100 mM. The new CG DNA model is suitable for modeling various biomolecular systems with adequate description of electrostatic and mechanical properties.
The influence of the local anesthetic lidocaine on electrostatic properties of a lipid membrane bilayer was studied by molecular dynamics simulations. The electrostatic dipole potential, charge ...densities, and orientations of the headgroup angle have been examined in the presence of different amounts of charged or uncharged forms of lidocaine. Important changes in the membrane properties caused by the presence of both forms of lidocaine are presented and discussed. Our simulations have shown that both charged and uncharged lidocaine cause almost the same increase in the electrostatic potential in the middle of the membrane, although for different reasons. The increase, ∼90
mV for 9
mol % of lidocaine and 220
mV for 28
mol % of lidocaine, is of a size that may affect the functioning of voltage-gated ion channels.
In order to investigate structural and dynamical properties of local anesthetic articaine in a model lipid bilayer, a series of molecular dynamics simulations have been performed. Simulations were ...carried out for neutral and charged (protonated) forms of articaine inserted in fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer. For comparison purpose, a fully hydrated DMPC bilayer without articaine was also simulated. The length of each simulation was 200
ns. Various properties of the lipid bilayer systems in the presence of both charged and uncharged forms of articaine taken at two different concentrations have been examined: membrane area per lipid, mass density distributions, order parameters, radial distribution functions, head group tilt, diffusion coefficients, electrostatic potential, etc, and compared with results of previous simulations of DMPC bilayer in the presence of lidocaine. It was shown that addition of both charged and neutral forms of articaine causes increase of the dipole electrostatic potential in the membrane interior.
Display omitted
► Articaine molecules are localized in the head group region of the bilayer. ► Articaine molecules are oriented predominantly parallel to the bilayer surface. ► Neutral form of articaine can penetrate across membrane. ► Presence of articaine increases dipole electrostatic potential inside bilayer.
The physical mechanism of the folding and unfolding of chromatin is fundamentally related to transcription but is incompletely characterized and not fully understood. We experimentally and ...theoretically studied chromatin compaction by investigating the salt-mediated folding of an array made of 12 positioning nucleosomes with 177 bp repeat length. Sedimentation velocity measurements were performed to monitor the folding provoked by addition of cations Na+, K+, Mg2+, Ca2+, spermidine3+, Co(NH3)63+, and spermine4+. We found typical polyelectrolyte behavior, with the critical concentration of cation needed to bring about maximal folding covering a range of almost five orders of magnitude (from 2 μM for spermine4+ to 100 mM for Na+). A coarse-grained model of the nucleosome array based on a continuum dielectric description and including the explicit presence of mobile ions and charged flexible histone tails was used in computer simulations to investigate the cation-mediated compaction. The results of the simulations with explicit ions are in general agreement with the experimental data, whereas simple Debye-Hückel models are intrinsically incapable of describing chromatin array folding by multivalent cations. We conclude that the theoretical description of the salt-induced chromatin folding must incorporate explicit mobile ions that include ion correlation and ion competition effects.
Systematic construction of coarse-grained molecular models from detailed atomistic simulations, and even from ab initio simulations is discussed. Atomistic simulations are first performed to extract ...structural information about the system, which is then used to determine effective potentials for a coarse-grained model of the same system. The statistical-mechanical equations expressing the canonical properties in terms of potential parameters can be inverted and solved numerically according to the iterative Newton scheme. In our previous applications, known as the Inverse Monte Carlo, radial distribution functions were inverted to reconstruct pair potential, while in a more general approach the targets can be other canonical averages. We have considered several examples of coarse-graining; for the united atom water model we suggest an easy way to overcome the known problem of high pressure. Further, we have developed coarse-grained models for L- and D-prolines, dissolved here in an organic solvent (dimethylsulfoxide), keeping their enantiomeric properties from the corresponding all-atom proline model. Finally, we have revisited the previously developed coarse-grained lipid model based on an updated all-atomic force field. We use this model in large-scale meso-scale simulations demonstrating spontaneous formation of different structures, such as vesicles, micelles, and multi-lamellar structures, depending on thermodynamical conditions.
To study the effects of water on conformational dynamics of polyalcohols, Molecular Dynamics simulations of glycerol–water liquid mixtures have been carried out at different concentrations: 42.9 and ...60.0 wt % of glycerol, respectively. On the basis of the analysis of backbone conformer distributions, it is found that the surrounding water molecules have a large impact on the populations of the glycerol conformers. While the local structure of water in the liquid mixture is surprisingly close to that in pure liquid water, the behavior of glycerols can be divided into three different categories where roughly 25% of them occur in a structure similar to that in pure liquid of glycerol, ca. 25% of them exist as monomers, solvated by water, and the remaining 50% of glycerols in the mixture form H-bonded strings as remains of the glycerol H-bond network. The typical glycerol H-bond network still exists even at the lower concentration of 40 wt % of glycerol. The microheterogeneity of water–glycerol mixtures is analyzed using time-averaged distributions of the sizes of the water aggregates. At 40 wt % of glycerol, the cluster sizes from 3 to 10 water molecules are observed. The increase of glycerol content causes a depletion of clusters leading to smaller 3–5 molecule clusters domination. Translational diffusion coefficients have been calculated to study the dynamical behavior of both glycerol and water molecules. Rotational-reorientational motion is studied both in overall and in selected substructures on the basis of time correlation functions. Characteristic time scales for different motional modes are deduced on the basis of the calculated correlation times. The general conclusion is that the presence of water increases the overall mobility of glycerol, while glycerol slows the mobility of water.
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.
In this work we propose a new force field for modelling of adsorption of CO2, N2, O2 and Ar in all silica and Na+ exchanged Si-Al zeolites. The force field has a standard molecular-mechanical ...functional form with electrostatic and Lennard-Jones interactions satisfying Lorentz-Berthelot mixing rules and thus has a potential for further extension in terms of new molecular types. The parameters for the zeolite framework atom types are optimized by an iterative procedure minimizing the difference with experimental adsorption data for a number of different zeolite structures and Si:Al ratios. The new force field shows a good agreement with available experimental data including those not used in the optimization procedure, and which also shows a reasonable transferability within different zeolite topologies. We suggest a potential usage in screening of different zeolite structures for carbon capture and storage process, and more generally, for separation of other gases.