Cationic DMPC/DMTAP lipid bilayers: molecular dynamics study (original) (raw)
Related papers
Biophysical Journal, 2007
Molecular dynamics simulations of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using the CHARMM27 force field in the tensionless isothermal-isobaric (NPT) ensemble give highly ordered, gel-like bilayers with an area per lipid of ;48 Å 2. To obtain fluid (L a) phase properties of DPPC bilayers represented by the CHARMM energy function in this ensemble, we reparameterized the atomic partial charges in the lipid headgroup and upper parts of the acyl chains. The new charges were determined from the electron structure using both the Mulliken method and the restricted electrostatic potential fitting method. We tested the derived charges in molecular dynamics simulations of a fully hydrated DPPC bilayer. Only the simulation with the new restricted electrostatic potential charges shows significant improvements compared with simulations using the original CHARMM27 force field resulting in an area per lipid of 60.4 6 0.1 Å 2. Compared to the 48 Å 2 , the new value of 60.4 Å 2 is in fair agreement with the experimental value of 64 Å 2. In addition, the simulated order parameter profile and electron density profile are in satisfactory agreement with experimental data. Thus, the biologically more interesting fluid phase of DPPC bilayers can now be simulated in all-atom simulations in the NPT ensemble by employing our modified CHARMM27 force field.
Effect of monovalent salt on cationic lipid membranes as revealed by molecular dynamics simulations
2005
An atomic-scale understanding of cationic lipid membranes is required for development of gene delivery agents based on cationic liposomes. To address this problem, we recently performed molecular dynamics (MD) simulations of mixed lipid membranes comprised of cationic dimyristoyltrimethylammonium propane (DMTAP) and zwitterionic dimyristoylphosphatidylcholine (DMPC) (Biophys. J. 2004, 86, 3461-3472). Given that salt ions are always present under physiological conditions, here we focus on the effects of monovalent salt (NaCl) on cationic (DMPC/DMTAP) membranes. Using atomistic MD simulations, we found that saltinduced changes in membranes depend strongly on their composition. When the DMTAP mole fraction is small (around 6%), the addition of monovalent salt leads to a considerable compression of the membrane and to a concurrent enhancement of the ordering of lipid acyl chains. That is accompanied by reorientation of phosphatidylcholine headgroups in the outward normal direction and slight changes in electrostatic properties. We attribute these changes to complexation of DMPC lipids with Na + ions which penetrate deep into the membrane and bind to the carbonyl region of the DMPC lipids. In contrast, at medium and high molar fractions of cationic DMTAP (50 and 75%) a substantial positive surface charge density of the membranes prevents the binding of Na + ions, making such membranes almost insensitive to monovalent salt. Finally, we compare our results to the Poisson-Boltzmann theory. With the exception of the immediate vicinity of the bilayer plane, we found excellent agreement with the theory. This is as expected since unlike in the theoretical description the surface is now structured due to its atomic scale nature.
Molecular dynamics simulation of the fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer
2003
The structural properties of lipid bilayers in biological membranes are of great interest in biochemistry, biophysics, and medicine. The main goal of the present study was to use molecular dynamic (MD) techniques to investigate physical properties of the hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer. -- The bilayer model consists of 25 DPPC molecules per each monolayer and 44.8% water by total weight. A modified version of AMBER MD suit of programs with CHARMM22 force field for phospholipids was used in simulation. The isothermal-isobaric or NPT ensemble with a fully flexible simulation box in ROAR program was used in this study. Simulations were performed under different pressure and temperature conditions. -- According to experimental results, a liquid crystal phase (Lα) is expected with the DPPC bilayer simulated under 1 atm pressure and 323 K temperature conditions. However, area per lipid, bilayer thickness, chain tilt, and the order parameters resulting from the prese...
Molecular dynamics simulation studies of lipid bilayer systems
Acta biochimica Polonica, 2000
The main structural element of biological membranes is a liquid-crystalline lipid bilayer. Other constituents, i.e. proteins, sterols and peptides, either intercalate into or loosely attach to the bilayer. We applied a molecular dynamics simulation method to study membrane systems at various levels of compositional complexity. The studies were started from simple lipid bilayers containing a single type phosphatidylcholine (PC) and water molecules (PC bilayers). As a next step, cholesterol (Chol) molecules were introduced to the PC bilayers (PC-Chol bilayers). These studies provided detailed information about the structure and dynamics of the membrane/water interface and the hydrocarbon chain region in bilayers built of various types of PCs and Chol. This enabled studies of membrane systems of higher complexity. They included the investigation of an integral membrane protein in its natural environment of a PC bilayer, and the antibacterial activity of magainin-2. The latter study req...
Molecular Dynamics Simulations of Phospholipid Bilayers
Journal of Biomolecular Structure and Dynamics, 1994
To investigate the microscopic interactions between cholesterol and lipids in biological membranes, we have performed a series of molecular dynamics simulations of large membranes with different levels of cholesterol content. The simulations extend to 10 ns, and were performed with hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers. The bilayers contain 1024 lipids of which 0-40% were cholesterol and the rest DPPC. The effects of cholesterol on the structure and mesoscopic dynamics of the bilayer were monitored as a function of cholesterol concentration. The main effects observed are a significant ordering of the DPPC chains (as monitored by NMR type order parameters), a reduced fraction of gauche bonds, a reduced surface area per lipid, less undulations-corresponding to an increased bending modulus for the membrane, smaller area fluctuations, and a reduced lateral diffusion of DPPC-lipids as well as cholesterols.
Colloids and Surfaces B: Biointerfaces, 2009
A molecular dynamics simulation study of the steady and dynamic properties of an asymmetric phospholipid bilayer was carried out in the presence of sodium or calcium ions. The asymmetric lipid bilayer was seen to resemble a cellular membrane of an eukaryotic cell, which was modeled by dipalmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylserine (DPPS), placing the DPPS in one of the two leaflets of the lipid bilayer. From a numerical analysis of the simulated trajectories, information was obtained with atomic resolution for both membrane leaflet concerning the effect of bilayer asymmetry on different properties of the lipid/water interface, such as the translational diffusion coefficient and rotational relaxation time of the water molecules, lipid hydration, and residence time of water around different lipid atoms. In addition, information related to lipid conformation, and lipid-lipid interactions was also analyzed.
Biophysical Journal, 2003
We study the influence of truncating the electrostatic interactions in a fully hydrated pure dipalmitoylphosphatidylcholine (DPPC) bilayer through 20ns molecular dynamics simulations. The computations in which the electrostatic interactions were truncated are compared to similar simulations using the particle-mesh Ewald (PME) technique. All examined truncation distances (1.8–2.5nm) lead to major effects on the bilayer properties, such as enhanced order of
Examining the Contributions of Lipid Shape and Headgroup Charge on Bilayer Behavior
Biophysical Journal, 2008
To better understand bilayer property dependency on lipid electrostatics and headgroup size, we use atomistic molecular dynamics simulations to study negatively charged and neutral lipid membranes. We compare the negatively charged phosphatidic acid (PA), which at physiological pH and salt concentration has a negative spontaneous curvature, with the negatively charged phosphatidylglycerol (PG) and neutrally charged phosphatidylcholine (PC), both of which have zero spontaneous curvature. The PA lipids are simulated using two different sets of partial charges for the headgroup and the varied charge distribution between the two PA systems results in significantly different locations for the Na 1 ions relative to the water/membrane interface. For one PA system, the Na 1 ions are localized around the phosphate group. In the second PA system, the Na 1 ions are located near the ester carbonyl atoms, which coincides with the preferred location site for the PG Na 1 ions. We find that the Na 1 ion location has a larger effect on bilayer fluidity properties than lipid headgroup size, where the A lipid and acyl chain order parameter values are more similar between the PA and PG bilayers that have Na 1 ions located near the ester groups than between the two PA bilayers.
Journal of Chemical Theory and Computation, 2010
Molecular dynamics simulations of fully hydrated pure bilayers of four widely studied phospholipids, 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and 2-oleoyl-1palmitoyl-sn-glycero-3-phosphocholine (POPC) using a recent revision of the GROMOS96 force field are reported. It is shown that the force field reproduces the structure and the hydration of bilayers formed by each of the four lipids with high accuracy. Specifically, the solvation and the orientation of the dipole of the phosphocholine headgroup and of the ester carbonyls show that the structure of the primary hydration shell in the simulations closely matches experimental findings. This work highlights the need to reproduce a broad range of properties beyond the area per lipid, which is poorly defined experimentally, and to consider the effect of system size and sampling times well beyond those commonly used.
The Journal of Chemical Physics, 2009
The electrostatic properties of lipid membranes are of profound importance as they are directly associated with membrane potential and, consequently, with numerous membrane-mediated biological phenomena. Here we address a number of methodological issues related to the computation of the electrostatic potential from atomic-scale molecular dynamics simulations of lipid bilayers. We discuss two slightly different forms of Poisson equation that are normally used to calculate the membrane potential: ͑i͒ a classical form when the potential and the electric field are chosen to be zero on one of the sides of a simulation box and ͑ii͒ an alternative form, when the potential is set to be the same on the opposite sides of a simulation box. Both forms differ by a position-dependent correction term, which has been shown to be proportional to the overall dipole moment of a bilayer system ͑for neutral systems͒. For symmetric bilayers we demonstrate that both approaches give essentially the same potential profiles, provided that simulations are long enough ͑a production run of at least 100 ns is required͒ and that fluctuations of the center of mass of a bilayer are properly accounted for. In contrast, for asymmetric lipid bilayers, the second approach is no longer appropriate due to a nonzero net dipole moment across a simulation box with a single asymmetric bilayer. We demonstrate that in this case the electrostatic potential can adequately be described by the classical form of Poisson equation, provided that it is employed in conjunction with tin-foil boundary conditions, which exactly balance a nonzero surface charge of a periodically replicated multibilayer system. Furthermore, we show that vacuum boundary conditions give qualitatively similar potential profiles for asymmetric lipid bilayers as compared to the conventional periodic boundaries, but accurate determination of the transmembrane potential difference is then hindered due to detachment of some water dipoles from bulk aqueous solution to vacuum.