Search

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. [ABSTRACT FROM AUTHOR]

In this work, the properties of saturated and unsaturated fatty acid acyl chains 16:0, 18:0, 18:1(n-9)cis, 18:2(n-6)cis, 18:3(n-3)cis, 18:4(n-3)cis, 18:5(n-3)cis, 20:4(n-6)cis, 20:5(n-3)cis and 22:6(n-3)cis in a bilayer liquid crystalline state and similar hydrocarbon chains (with CH3 terminal groups instead of C=O groups) in the unperturbed state characterised by a lack of long-range interaction were investigated. The unperturbed hydrocarbon chains were modelled by Monte Carlo simulations at temperature T=303 K; sixteen fully hydrated homogeneous liquid crystalline phosphatidylcholine bilayers containing these chains were studied by molecular dynamics simulations at the same temperature. To eliminate effects of the simulation parameters, the molecular dynamics and Monte Carlo simulations were carried out using the same structural data and force field coefficients. From these computer simulations, the average distances between terminal carbon atoms of the chains (end-to-end distances) were calculated and compared. The trends in the end-to-end distances obtained for the unperturbed chains were found to be qualitatively similar to those obtained for the same lipid chains in the bilayers. So, for understanding of a number of processes in biological membranes (e.g., changes in fatty acid composition caused by environmental changes such as temperature and pressure), it is possible to use, at least as a first approximation, the relationships between the structure and properties for unperturbed or isolated hydrocarbon chains. [ABSTRACT FROM AUTHOR]

To provide insight into the molecular mechanisms of local anesthetic action, we have carried out an extensive investigation of two amide type local anesthetics, lidocaine and articaine in both charged and uncharged forms, interacting with DMPC lipid membrane. We have applied both standard molecular dynamics simulations and metadynamics simulations to provide a detailed description of the free energy landscape of anesthetics embedded in the lipid bilayer. The global minimum of the free energy surface (equilibrium position of anesthetics in the lipid membrane) occurred around 1nm of the bilayer center. The uncharged anesthetics show more affinity to bind to this region compared to the charged drugs. The binding free energy of uncharged lidocaine in the membrane (−30.3kJ/mol) is higher than uncharged articaine (−24.0kJ/mol), which is in good agreement with higher lipid solubility of lidocaine relative to the articaine. The octanol/water partition coefficient of uncharged drugs was also investigated using expanded ensemble simulations. In addition, complementary standard MD simulations were carried out to study the partitioning behavior of multiple anesthetics inside the lipid bilayer. The results obtained here are in line with previously reported simulations and suggest that the different forms of anesthetics induce different structural modifications in the lipid bilayer, which can provide new insights into their complex membrane translocation phenomena. [ABSTRACT FROM AUTHOR]

The all-atomic force field Slipids (Stockholm Lipids) for lipid bilayers simulations has been extended to polyunsaturated lipids. Following the strategy adopted in the development of previous versions of the Slipids force field, the parametrization was essentially based on high-level ab initio calculations. Atomic charges and torsion angles related to polyunsaturated lipid tails were parametrized using structures of dienes molecules. The new parameters of the force field were validated in simulations of bilayers composed of seven polyunsaturated lipids. An overall good agreement was found with available experimental data on the areas per lipids, volumetric properties of bilayers, deuterium order parameters, and scattering form factors. Furthermore, simulations of bilayers consisting of highly polyunsaturated lipids and cholesterol molecules have been carried out. The majority of cholesterol molecules were found in a position parallel to bilayer normal with the hydroxyl group directed to the membrane surface, while a small fraction of cholesterol was found in the bilayer center parallel to the membrane plane. Furthermore, a tendency of cholesterol molecules to form chain-like clusters in polyunsaturated bilayers was qualitatively observed. [ABSTRACT FROM AUTHOR]

With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg. [ABSTRACT FROM AUTHOR]

We have used systematic structure-based coarse graining to derive effective site-site potentials for a 10-site coarse-grained dimyristoylphosphatidylcholine (DMPC) lipid model and investigated their state point dependence. The potentials provide for the coarse-grained model the same site-site radial distribution functions, bond and angle distributions as those computed in atomistic simulations carried out at four different lipid-water molar ratios. It was shown that there is a non-negligible dependence of the effective potentials on the concentration at which they were generated, which is also manifested in the properties of the lipid bilayers simulated using these potentials. Thus, effective potentials computed at low lipid concentration favor to more condensed and ordered structure of the bilayer with lower average area per lipid, while potentials obtained at higher lipid concentrations provide more fluid-like structure. The best agreement with the reference data and experiment was achieved using the set of potentials derived from atomistic simulations at 1:30 lipid:water molar ratio providing fully saturated hydration of DMPC lipids. Despite theoretical limitations of pairwise coarse-grained potentials expressed in their state point dependence, all the resulting potentials provide a stable bilayer structure with correct partitioning of different lipid groups across the bilayer as well as acceptable values of the average lipid area, compressibility and orientational ordering. In addition to bilayer simulations, the model has proven its robustness in modeling of self-aggregation of lipids from randomly dispersed solution to ordered bilayer structures, bicelles, and vesicles. © 2014 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

An all-atomistic force field (FF) has been developedfor fullysaturated phospholipids. The parametrization has been largely basedon high-level ab initio calculations in order to keep the empiricalinput to a minimum. Parameters for the lipid chains have been developedbased on knowledge about bulk alkane liquids, for which thermodynamicand dynamic data are excellently reproduced. The FFs ability to simulatelipid bilayers in the liquid crystalline phase in a tensionless ensemblewas tested in simulations of three lipids: 1,2-diauroyl-sn-glycero-3-phospocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1,2-dipalmitoyl-sn-glycero-3-phospcholine (DPPC). Computed areas and volumes per lipid,and three different kinds of bilayer thicknesses, have been investigated.Most importantly NMR order parameters and scattering form factorsagree in an excellent manner with experimental data under a rangeof temperatures. Further, the compatibility with the AMBER FF forbiomolecules as well as the ability to simulate bilayers in gel phasewas demonstrated. Overall, the FF presented here provides the importantbalance between the hydrophilic and hydrophobic forces present inlipid bilayers and therefore can be used for more complicated studiesof realistic biological membranes with protein insertions. [ABSTRACT FROM AUTHOR]

Abstract: 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 200ns. 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. [ABSTRACT FROM AUTHOR]

The CHARMM force field for DMPC lipids was modified in order to improve agreement with experiment for a number of important properties of hydrated lipid bilayer. The modification consists in introduction of a scaling factor 0.83 for 1–4 electrostatic interactions (between atoms separated by three covalent bonds), which provides correct transgauche ratio in the alkane tails, and recalculation of the headgroup charges on the basis of HF/6-311(d,p) ab-initio computations. Both rigid TIP3P and flexible SPC water models were used with the new lipid model, showing similar results. The new model in a 75 ns simulation has shown a correct value of the area per lipid at zero surface tension, as well as good agreement with the experiment for the electron density, structure factor, and order parameters, including those in the headgroup part of lipids. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2008 [ABSTRACT FROM AUTHOR]

Abstract: Interactions between anesthetics (lidocaine and short chain alcohols) and lipid membranes formed by dimyristoylphosphatidylcholine (DMPC) were studied using NMR spectroscopy. The orientational order of lidocaine was investigated using deuterium NMR on a selectively labelled compound whereas segmental ordering in the lipids was probed by two-dimensional 1H–13C separated local field experiments under magic-angle spinning conditions. In addition, trajectories generated in molecular dynamics (MD) computer simulations were used for interpretation of the experimental results. Separate simulations were carried out with charged and uncharged lidocaine molecules. Reasonable agreement between experimental dipolar interactions and the calculated counterparts was observed. Our results clearly show that charged lidocaine affects significantly the lipid headgroup. In particular the ordering of the lipids is increased accompanied by drastic changes in the orientation of the P–N vector in the choline group. [Copyright &y& Elsevier]