A significant modification to the additive all-atom CHARMM lipid force field (FF) is developed and applied to phospholipid bilayers with both choline and ethanolamine containing head groups and with both saturated and unsaturated aliphatic chains. Motivated by the current CHARMM lipid FF (C27 and C27r) systematically yielding values of the surface area per lipid that are smaller than experimental estimates and gel-like structures of bilayers well above the gel transition temperature, selected torsional, Lennard-Jones and partial atomic charge parameters were modified by targeting both quantum mechanical (QM) and experimental data. QM calculations ranging from high-level ab initio calculations on small molecules to semi-empirical QM studies on a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer in combination with experimental thermodynamic data were used as target data for parameter optimization. These changes were tested with simulations of pure bilayers at high hydration of the following six lipids: DPPC, 1,2-dimyristoyl-sn-phosphatidylcholine (DMPC), 1,2-dilauroyl-sn-phosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-phosphatidylcholine (POPC), 1,2-dioleoyl-sn-phosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-sn-phosphatidylethanolamine (POPE); simulations of a low hydration DOPC bilayer were also performed. Agreement with experimental surface area is on average within 2%, and the density profiles agree well with neutron and x-ray diffraction experiments. NMR deuterium order parameters (SCD) are well predicted with the new FF, including proper splitting of the SCD for the aliphatic carbon adjacent to the carbonyl for DPPC, POPE, and POPC bilayers. The area compressibility modulus and frequency dependence of 13C NMR relaxation rates of DPPC, and the water distribution of low hydration DOPC bilayers also agree well with experiment. Accordingly, the presented lipid FF, referred to as C36, allows for molecular dynamics simulations to be run in the tensionless ensemble (NPT), and is anticipated to be of utility for simulations of pure lipids systems as well as heterogeneous systems including membrane proteins.

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

CHARMM-GUI Membrane Builder, http://www.charmm-gui.org/input/membrane, is a web-based user interface designed to interactively build all-atom protein/membrane or membrane-only systems for molecular dynamics simulation through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builderthat allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types such as phosphoinositides, cardiolipin, sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance-based algorithm for faster initial ion displacement, (5) CHARMM inputs for P21 image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: cardiolipin lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM-GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments.

The CHARMM-GUI Membrane Builder (http://www.charmm-gui.org/input/membrane), an intuitive, straightforward, web-based graphical user interface, was expanded to automate the building process of heterogeneous lipid bilayers, with or without a protein and with support for up to 32 different lipid types. The efficacy of these new features was tested by building and simulating lipid bilayers that resemble yeast membranes, composed of cholesterol, dipalmitoylphosphatidylcholine, dioleoylphosphatidylcholine, palmitoyloleoylphosphatidylethanolamine, palmitoyloleoylphosphatidylamine, and palmitoyloleoylphosphatidylserine. Four membranes with varying concentrations of cholesterol and phospholipids were simulated, for a total of 170 ns at 303.15 K. Unsaturated phospholipid chain concentration had the largest influence on membrane properties, such as average lipid surface area, density profiles, deuterium order parameters, and cholesterol tilt angle. Simulations with a high concentration of unsaturated chains (73%, membrane(unsat)) resulted in a significant increase in lipid surface area and a decrease in deuterium order parameters, compared with membranes with a high concentration of saturated chains (60-63%, membrane(sat)). The average tilt angle of cholesterol with respect to bilayer normal was largest, and the distribution was significantly broader for membrane(unsat). Moreover, short-lived cholesterol orientations parallel to the membrane surface existed only for membrane(unsat). The membrane(sat) simulations were in a liquid-ordered state, and agree with similar experimental cholesterol-containing membranes.

Energies of 119 conformations of normal alkanes from butane to heptane were calculated at approximately the CCSD(T)/cc-pVQZ level. Energies of gauche (g) conformers relative to trans (t) decrease as chain length increases. In what is termed the "positive pentane effect", adjacent gauche conformers of the same sign are stabilized compared to nonadjacent conformers; e.g., for hexane the energies of tgt, tgg, and gtg are 0.600, 0.930, and 1.18 kcal/mol, respectively. Torsional terms in the CHARMM27 (C27) force field were fit to the calculated QM energies to yield a revised potential, C27r. Molecular dynamics simulations of normal alkanes (heptane, decane, tridecane, and pentadecane) with C27r yield higher populations of gauche states, increased transition rates, and improved agreement with experiment as compared to C27. In addition, C27r simulations of a hydrated DPPC lipid bilayer yield improved agreement with the experimental NMR deuterium order parameters for the aliphatic chain ends.

No abstract

No abstract

Glycolipids (such as glycoglycerolipids, glycosphingolipids, and glycosylphosphatidylinositol) and lipoglycans (such as lipopolysaccharides (LPS), lipooligosaccharides (LOS), mycobacterial lipoarabinomannan, and mycoplasma lipoglycans) are typically found on the surface of cell membranes and play crucial roles in various cellular functions. Characterizing their structure and dynamics at the molecular level is essential to understand their biological roles, but systematic generation of glycolipid and lipoglycan structures is challenging because of great variations in lipid structures and glycan sequences (i.e., carbohydrate types and their linkages). To facilitate the generation of all-atom glycolipid/LPS/LOS structures, we have developed Glycolipid Modeler and LPS Modeler in CHARMM-GUI (http://www.charmm-gui.org), a webbased interface that simplifies building of complex biological simulation systems. In addition, we have incorporated these modules into Membrane Builder, so that users can readily build a complex symmetric or asymmetric biological membrane system with various glycolipids and LPS/LOS. These tools are expected to be useful in innovative and novel glycolipid/LPS/LOS modeling and simulation research by easing tedious and intricate steps in modeling complex biological systems, and shall provide insight into structures, dynamics, and underlying mechanisms of complex glycolipid-/LPS-/LOS-containing biological membrane systems.