Rainer A. Böckmann

Biomolecular simulations of membranes: Physical properties from different force fields

Phospholipid force fields are of ample importance for the simulation of artificial bilayers, membranes, and also for the simulation of integral membrane proteins. Here, we compare the two most applied atomic force fields for phospholipids, the all-atom CHARMM27 and the united atom Berger force field, with a newly developed all-atom generalized AMBER force field (GAFF) for dioleoylphosphatidylcholine molecules. Only the latter displays the experimentally observed difference in the order of the C2 atom between the two acyl chains. The interfacial water dynamics is smoothly increased between the lipid carbonyl region and the bulk water phase for all force fields; however, the water order and with it the electrostatic potential across the bilayer showed distinct differences between the force fields. Both Berger and GAFF underestimate the lipid self-diffusion. GAFF offers a consistent force field for the atomic scale simulation of biomembranes.

Predicting free energy changes using structural ensembles

editorial office: 75 Varick Street, Fl 9, New York, NY 10013-1917. Tel: (212) 726 9200, Fax: (212) 689 9702. annual subscription rates: USA/ Canada: US

Kinetics, Statistics, and Energetics of Lipid Membrane Electroporation Studied by Molecular Dynamics Simulations

150 (personal), US

Going Backward: A Flexible Geometric Approach to Reverse Transformation from Coarse Grained to Atomistic Models

2,185 (institution), Canada add 7% GST #104911595RT001; Euro-zone: €153 (personal), €1,736. (institution); UK and Europe £99 (personal), £1,120 (institution); Rest of world (excluding China, Japan, Korea): £99 (personal), £1,120 (institution); Japan: Contact NPG Nature Asia-Pacific, Chiyoda Building, 2-37 Ichigayatamachi, Shinjuku-ku, Tokyo 162-0843. Tel: 81 (03) 3267 8751, Fax: 81 (03) 3267 8746. Back issues: US

Computational Lipidomics with insane : A Versatile Tool for Generating Custom Membranes for Molecular Simulations

20, Canada add 7% for GST.

Optimization of the OPLS-AA Force Field for Long Hydrocarbons

Membrane electroporation is the method to directly transfer bioactive substances such as drugs and genes into living cells, as well as preceding electrofusion. Although much information on the microscopic mechanism has been obtained both from experiment and simulation, the existence and nature of possible intermediates is still unclear. To elucidate intermediates of electropore formation by direct comparison with measured prepore formation kinetics, we have carried out 49 atomistic electroporation simulations on a palmitoyl-oleoyl-phosphatidylcholine bilayer for electric field strengths between 0.04 and 0.7 V/nm. A statistical theory is developed to facilitate direct comparison of experimental (macroscopic) prepore formation kinetics with the (single event) preporation times derived from the simulations, which also allows us to extract an effective number of lipids involved in each pore formation event. A linear dependency of the activation energy for prepore formation on the applied field is seen, with quantitative agreement between experiment and simulation. The distribution of preporation times suggests a four-state pore formation model. The model involves a first intermediate characterized by a differential tilt of the polar lipid headgroups on both leaflets, and a second intermediate (prepore), where a polar chain across the bilayer is formed by 3-4 lipid headgroups and several water molecules, thereby providing a microscopic explanation for the polarizable volume derived previously from the measured kinetics. An average pore radius of 0.47 +/- 0.15 nm is seen, in favorable agreement with conductance measurements and electrooptical experiments of lipid vesicles.

The closure of Pak1‐dependent macropinosomes requires the phosphorylation of CtBP1/BARS

The conversion of coarse-grained to atomistic models is an important step in obtaining insight about atomistic scale processes from coarse-grained simulations. For this process, called backmapping or reverse transformation, several tools are available, but these commonly require libraries of molecule fragments or they are linked to a specific software package. In addition, the methods are usually restricted to specific molecules and to a specific force field. Here, we present an alternative method, consisting of geometric projection and subsequent force-field based relaxation. This method is designed to be simple and flexible, and offers a generic solution for resolution transformation. For simple systems, the conversion only requires a list of particle correspondences on the two levels of resolution. For special cases, such as nondefault protonation states of amino acids and virtual sites, a target particle list can be specified. The mapping uses simple building blocks, which list the particles on the different levels of resolution. For conversion to higher resolution, the initial model is relaxed with several short cycles of energy minimization and position-restrained MD. The reconstruction of an atomistic backbone from a coarse-grained model is done using a new dedicated algorithm. The method is generic and can be used to map between any two particle based representations, provided that a mapping can be written. The focus of this work is on the coarse-grained MARTINI force field, for which mapping definitions are written to allow conversion to and from the higher-resolution force fields GROMOS, CHARMM, and AMBER, and to and from a simplified three-bead lipid model. Together, these offer the possibility to simulate mesoscopic membrane structures, to be transformed to MARTINI and subsequently to an atomistic model for investigation of detailed interactions. The method was tested on a set of systems ranging from a simple, single-component bilayer to a large protein-membrane-solvent complex. The results demonstrate the efficiency and the efficacy of the new approach.

Effects of Alkali Cations and Halide Anions on the DOPC Lipid Membrane

For simulations of membranes and membrane proteins, the generation of the lipid bilayer is a critical step in the setup of the system. Membranes comprising multiple components pose a particular challenge, because the relative abundances need to be controlled and the equilibration of the system may take several microseconds. Here we present a comprehensive method for building membrane containing systems, characterized by simplicity and versatility. The program uses preset, coarse-grain lipid templates to build the membrane, and also allows on-the-fly generation of simple lipid types by specifying the headgroup, linker, and lipid tails on the command line. The resulting models can be equilibrated, after which a relaxed atomistic model can be obtained by reverse transformation. For multicomponent membranes, this provides an efficient means for generating equilibrated atomistic models. The method is called insane, an acronym for INSert membrANE. The program has been made available, together with the complementary method for reverse transformation, at http://cgmartini.nl/ . This work highlights the key features of insane and presents a survey of properties for a large range of lipids as a start of a computational lipidomics project.

Mechanism of interaction of monovalent ions with phosphatidylcholine lipid membranes.

The all-atom optimized potentials for liquid simulations (OPLS-AA) force field is a popular force field for simulating biomolecules. However, the current OPLS parameters for hydrocarbons developed using short alkanes cannot reproduce the liquid properties of long alkanes in molecular dynamics simulations. Therefore, the extension of OPLS-AA to (phospho)lipid molecules required for the study of biological membranes was hampered in the past. Here, we optimized the OPLS-AA force field for both short and long hydrocarbons. Following the framework of the OPLS-AA parametrization, we refined the torsional parameters for hydrocarbons by fitting to the gas-phase ab initio energy profiles calculated at the accurate MP2/aug-cc-pVTZ theory level. Additionally, the depth of the Lennard-Jones potential for methylene hydrogen atoms was adjusted to reproduce the densities and the heats of vaporization of alkanes and alkenes of different lengths. Optimization of partial charges finally allowed to reproduce the gel-to-liquid-phase transition temperature for pentadecane and solvation free energies. It is shown that the optimized parameter set (L-OPLS) yields improved hydrocarbon diffusion coefficients, viscosities, and gauche-trans ratios. Moreover, its applicability for lipid bilayer simulations is shown for a GMO bilayer in its liquid-crystalline phase.

HLA-B27 Subtypes Differentially Associated with Disease Exhibit Conformational Differences in Solution

Membrane fission is an essential process in membrane trafficking and other cellular functions. While many fissioning and trafficking steps are mediated by the large GTPase dynamin, some fission events are dynamin independent and involve C‐terminal‐binding protein‐1/brefeldinA‐ADP ribosylated substrate (CtBP1/BARS). To gain an insight into the molecular mechanisms of CtBP1/BARS in fission, we have studied the role of this protein in macropinocytosis, a dynamin‐independent endocytic pathway that can be synchronously activated by growth factors. Here, we show that upon activation of the epidermal growth factor receptor, CtBP1/BARS is (a) translocated to the macropinocytic cup and its surrounding membrane, (b) required for the fission of the macropinocytic cup and (c) phosphorylated on a specific serine that is a substrate for p21‐activated kinase, with this phosphorylation being essential for the fission of the macropinocytic cup. Importantly, we also show that CtBP1/BARS is required for macropinocytic internalization and infection of echovirus 1. These results provide an insight into the molecular mechanisms of CtBP1/BARS activation in membrane fissioning, and extend the relevance of CtBP1/BARS‐induced fission to human viral infection.