An all-atomistic force field (FF) has been developed for fully saturated phospholipids. The parametrization has been largely based on high-level ab initio calculations in order to keep the empirical input to a minimum. Parameters for the lipid chains have been developed based on knowledge about bulk alkane liquids, for which thermodynamic and dynamic data are excellently reproduced. The FFs ability to simulate lipid bilayers in the liquid crystalline phase in a tensionless ensemble was 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 factors agree in an excellent manner with experimental data under a range of temperatures. Further, the compatibility with the AMBER FF for biomolecules as well as the ability to simulate bilayers in gel phase was demonstrated. Overall, the FF presented here provides the important balance between the hydrophilic and hydrophobic forces present in lipid bilayers and therefore can be used for more complicated studies of realistic biological membranes with protein insertions.

Biological membranes are versatile in composition and host intriguing molecular processes. In order to be able to study these systems, an accurate model Hamiltonian or force field (FF) is a necessity. Here, we report the results of our extension of earlier developed all-atomistic FF parameters for fully saturated phospholipids that complements an earlier parameter set for saturated phosphatidylcholine lipids (J. Phys. Chem. B, 2012, 116, 3164-3179). The FF, coined Slipids (Stockholm lipids), now also includes parameters for unsaturated phosphatidylcholine and phosphatidylethanolamine lipids, e.g., POPC, DOPC, SOPC, POPE, and DOPE. As the extended set of parameters is derived with the same philosophy as previously applied, the resulting FF has been developed in a fully consistent manner. The capabilities of Slipids are demonstrated by performing long simulations without applying any surface tension and using the correct isothermal-isobaric (NPT) ensemble for a range of temperatures and carefully comparing a number of properties with experimental findings. Results show that several structural properties are very well reproduced, such as scattering form factors, NMR order parameters, thicknesses, and area per lipid. Thermal dependencies of different thicknesses and area per lipid are reproduced as well Lipid diffusion is systematically slightly underestimated, whereas the normalized lipid diffusion follows the experimental trends. This is believed to be due to the lack of collective movement in the relatively small bilayer patches used Furthermore, the compatibility with amino acid FFs from the AMBER family is tested in explicit transmembrane complexes of the WALP23 peptide with DLPC and DOPC bilayers, and this shows that Slipids can be used to study more complex and biologically relevant systems.

To be able to model complex biological membranes in a more realistic manner, the force field Slipids (Stockholm lipids) has been extended to include parameters for sphingomyelin (SM), phosphatidylglycerol (PG), phosphatidylserine (PS) lipids, and cholesterol. Since the parametrization scheme was faithful to the scheme used in previous editions of Slipids, all parameters are consistent and fully compatible. The results of careful validation of a number of key structural properties for one and two component lipid bilayers are in excellent agreement with experiments. Potentials of mean force for transferring water across binary mixtures of lipids and cholesterol were also computed in order to compare water permeability rates to experiments. In agreement with experimental and simulation studies, it was found that the permeability and partitioning of water is affected by cholesterol in lipid bilayers made of saturated lipids to the largest extent. With the extensions of Slipids presented here, it is now possible to study complex systems containing many different lipids and proteins in a fully atomistic resolution in the isothermic-isobaric (NPT) ensemble, which is the proper ensemble for membrane simulations.

We propose an effective and straightforward way of including atomic polarization in simulations of the partitioning of small molecules in inhomogenous media based on classical molecular dynamics with non-polarizable force fields. The approach presented here takes advantage of the relatively fast sampling of phase space obtained with additive force fields by adding the polarization effects afterwards. By using pre-polarized charges for the polar and non-polar phases together with a polarization correction term the effects of atomic polarization are effectively taken into account. The results show a clear improvement compared to using the more common setup with one set of charges obtained from gas phase ab initio calculations. It is shown that when proper measures are taken into account computer simulations with non-polarizable force fields are able to accurately determine water-membrane partitioning and preferential location of small molecules in the membrane interior. We believe that the approach presented here can be useful in rational drug design and in investigations of molecular mechanisms of anesthetic or toxic action.

Free energy calculations are vital for our understanding of biological processes on an atomistic scale and can offer insight to various mechanisms. However, in some cases, degrees of freedom (DOFs) orthogonal to the reaction coordinate have high energy barriers and/or long equilibration times, which prohibit proper sampling. Here we identify these orthogonal DOFs when studying the transfer of a solute from water to a model membrane. Important DOFs are identified in bulk liquids of different dielectric nature with metadynamics simulations and are used as reaction coordinates for the translocation process, resulting in two- and three-dimensional space of reaction coordinates. The results are in good agreement with experiments and elucidate the pitfalls of using one-dimensional reaction coordinates. The calculations performed here offer the most detailed free energy landscape of solutes embedded in lipid bilayers to date and show that free energy calculations can be used to study complex membrane translocation phenomena.