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 semiempirical 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: ...

Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types

http://www.bio21.bas.bg/ibf/PC_review.pdf

Phases and phase transitions of the phosphatidylcholines

Structure of the inverted hexagonal (HII) phase, and non-lamellar phase transitions of lipids.

The AMBER lipid force field has been updated to create Lipid14, allowing tensionless simulation of a number of lipid types with the AMBER MD package. The modular nature of this force field allows numerous combinations of head and tail groups to create different lipid types, enabling the easy insertion of new lipid species. The Lennard-Jones and torsion parameters of both the head and tail groups have been revised and updated partial charges calculated. The force field has been validated by simulating bilayers of six different lipid types for a total of 0.5 μs each without applying a surface tension; with favorable comparison to experiment for properties such as area per lipid, volume per lipid, bilayer thickness, NMR order parameters, scattering data, and lipid lateral diffusion. As the derivation of this force field is consistent with the AMBER development philosophy, Lipid14 is compatible with the AMBER protein, nucleic acid, carbohydrate, and small molecule force fields.

/pdf/lipid14-the-amber-lipid-force-field-1gj0t3c4bx.pdf

Lipid14: The Amber Lipid Force Field

Quantitative structures are obtained at 30 degrees C for the fully hydrated fluid phases of palmitoyloleoylphosphatidylcholine (POPC), with a double bond on the sn-2 hydrocarbon chain, and for dierucoylphosphatidylcholine (di22:1PC), with a double bond on each hydrocarbon chain. The form factors F(qz) for both lipids are obtained using a combination of three methods. (1) Volumetric measurements provide F(0). (2) X-ray scattering from extruded unilamellar vesicles provides /F(qz)/ for low q(z). (3) Diffuse X-ray scattering from oriented stacks of bilayers provides /F(qz)/ for high q(z). Also, data using method (2) are added to our recent data for dioleoylphosphatidylcholine (DOPC) using methods (1) and (3); the new DOPC data agree very well with the recent data and with (4) our older data obtained using a liquid crystallographic X-ray method. We used hybrid electron density models to obtain structural results from these form factors. The result for area per lipid (A) for DOPC 72.4 +/- 0.5 A(2) agrees well with our earlier publications, and we find A = 69.3 +/- 0.5 A2 for di22:1PC and A = 68.3 +/- 1.5 A2 for POPC. We obtain the values for five different average thicknesses: hydrophobic, steric, head-head, phosphate-phosphate and Luzzati. Comparison of the results for these three lipids and for our recent dimyristoylphosphatidylcholine (DMPC) determination provides quantitative measures of the effect of unsaturation on bilayer structure. Our results suggest that lipids with one monounsaturated chain have quantitative bilayer structures closer to lipids with two monounsaturated chains than to lipids with two completely saturated chains.

http://www.norbbi.com/public/publications/Kuc_JMB(2006).pdf

Structure of Fully Hydrated Fluid Phase Lipid Bilayers with Monounsaturated Chains

The thermotropic behavior and molecular properties of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) and 1-palmitoyl-2-dihydrosterculoyl-sn-glycero-3-phosphoethanolamine (PDSPE) have been investigated by 2H NMR spectroscopy using samples selectively labeled at the 5'-, 9'-, 10'-, and 16'-positions of the sn-2 chains. Comparison with the corresponding phosphocholine analogues (POPC and PDSPC), obtained as intermediate synthetic products, was used to monitor the role of the polar head group. Replacement of the choline moiety by ethanolamine increased the gel to liquid-crystal transition temperature by 10-32 degrees C and led to a significantly higher ordering of the fatty acyl chains in the liquid-crystalline bilayer state. The lateral compression effect, due to the smaller area per polar head group in PE, results in a bilayer to hexagonal phase transition at elevated temperatures. The effects on both PC and PE due to replacement of the olefinic group by a cyclopropane unit are similar. A decrease in the temperature of the gel to liquid-crystal phase transition, Tc, is observed upon introduction of a cyclopropane ring; it goes from 26 degrees C in POPE to approximately 10 degrees C in PDSPE. In addition, a very significant broadening of the transition profile is observed. These observations are consistent with the poor packing ability of mixed saturated and cyclopropane-containing chains due to the bulky substituent effect. The temperature of the bilayer-hexagonal phase transition of PE samples was decreased by 15-20 degrees C on replacement of oleoyl chains by dihydrosterculoyl chains at the sn-2 position.(ABSTRACT TRUNCATED AT 250 WORDS)

Effects of replacement of a double bond by a cyclopropane ring in phosphatidylethanolamines: a 2H NMR study of phase transitions and molecular organization.

The dynamical behavior of the acyl chains of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and 1-palmitoyl-2-dihydrosterculoyl-sn-glycero-3-phosphoethanolamine has been investigated by using 2H T1 and T2 relaxation times. Lipids were labeled at the 5-,9-,10-, and 16-positions of the sn-2 acyl chain. The profile of deuterium spin-lattice relaxation rate (T1(-1) vs. chain position is characterized in all systems by a marked discontinuity at the positions of the carbon-carbon double bond and the cyclopropane ring; the deuterons at these positions have relaxation rates which are greater than at any other labeled position of the sn-2 chain. For both types of sn-2 acyl chain, assuming a single-exponential correlation time and that the motion is within the rapid regime, the phosphatidylcholine lipid systems are less mobile than their phosphatidylethanolamine analogues. Systems containing an oleoyl chain are more dynamic than their analogues containing a dihydrosterculoyl chain. The rates of motion of the sn-2 acyl chains of phosphatidylethanolamine in a bilayer structure are slower than those of the lipid in an inverted hexagonal structure. In the hexagonal phase, the motional rates of a dihydrosterculoyl chain are slower than those of the corresponding positions of an oleoyl chain.

Acyl chain dynamics of phosphatidylethanolamines containing oleic acid and dihydrosterculic acid: 2H NMR relaxation studies.

Improved reliable methods for the partial chemical synthesis of phosphatidylcholines and phosphatidylethanolamines are described. Glycerophosphorylcholine (GPC) as its cadmium chloride complex was acylated rapidly (3 hr) with deuterated acyl anhydrides in dimethylformamide (DMF) to give phosphatidylcholines in > 90% yields. The products were free of sn-1 and sn-2 isomers. 2-Lysophosphatidylcholines were rapidly acylated with acyl anhydrides, used in nearly stoichiometric amounts, in chloroform affording mixed-chain phosphatidylcholines in > 85% yield. Pure phosphatidylethanolamines are prepared routinely in > 70% yield by an improved procedure from the corresponding phospatidylcholines using phospholipase-D.

Bruno Perly

Papers

Effects of replacement of a double bond by a cyclopropane ring in phosphatidylethanolamines: a 2H NMR study of phase transitions and molecular organization.

Acyl chain dynamics of phosphatidylethanolamines containing oleic acid and dihydrosterculic acid: 2H NMR relaxation studies.

Facile and high yielding syntheses of phosphatidylcholines and phosphatidylethanolamines containing 2H-labeled acyl chains