We study the influence of truncating the electrostatic interactions in a fully hydrated pure dipalmitoylphosphatidylcholine (DPPC) bilayer through 20 ns 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.5 nm) lead to major effects on the bilayer properties, such as enhanced order of acyl chains together with decreased areas per lipid. The results obtained using PME, on the other hand, are consistent with experiments. These artifacts are interpreted in terms of radial distribution functions g(r) of molecules and molecular groups in the bilayer plane. Pronounced maxima or minima in g(r) appear exactly at the cutoff distance indicating that the truncation gives rise to artificial ordering between the polar phosphatidyl and choline groups of the DPPC molecules. In systems described using PME, such artificial ordering is not present.

Extensive microscopic molecular dynamics simulations have been performed to study the effects of short-chain alcohols, methanol and ethanol, on two different fully hydrated lipid bilayer systems (POPC and DPPC) in the fluid phase at 323 K. It is found that ethanol has a stronger effect on the structural properties of the membranes. In particular, the bilayers become more fluid and permeable: ethanol molecules are able to penetrate through the membrane in typical timescales of approximately 200 ns, whereas for methanol that timescale is considerably longer, at least of the order of microseconds. A closer examination exposes a number of effects due to ethanol. Hydrogen-bonding analysis reveals that a large fraction of ethanols is involved in hydrogen bonds with lipids. This in turn is intimately coupled to the ordering of hydrocarbon chains: we find that binding to an ethanol decreases the order of the chains. We have also determined the dependence of lipid-chain ordering on ethanol concentration and found that to be nonmonotonous. Overall, we find good agreement with NMR and micropipette studies.

We employ 100-ns molecular dynamics simulations to study the influence of cholesterol on structural and dynamic properties of dipalmitoylphosphatidylcholine bilayers in the fluid phase. The effects of the cholesterol content on the bilayer structure are considered by varying the cholesterol concentration between 0 and 50%. We concentrate on the free area in the membrane and investigate quantities that are likely to be affected by changes in the free area and free volume properties. It is found that cholesterol has a strong impact on the free area properties of the bilayer. The changes in the amount of free area are shown to be intimately related to alterations in molecular packing, ordering of phospholipid tails, and behavior of compressibility moduli. Also the behavior of the lateral diffusion of both dipalmitoylphosphatidylcholine and cholesterol molecules with an increasing amount of cholesterol can in part be understood in terms of free area. Summarizing, our results highlight the central role of free area in comprehending the structural and dynamic properties of membranes containing cholesterol.

We provide compelling evidence that different treatments of electrostatic interactions in molecular dynamics simulations may dramatically affect dynamic properties of lipid bilayers. To this end, we consider a fully hydrated pure dipalmitoylphosphatidylcholine bilayer through 50-ns molecular dynamics simulations and study various dynamic properties of individual lipids in a membrane, including the velocity autocorrelation function, the lateral and rotational diffusion coefficients, and the autocorrelation function for the area per molecule. We compare the results based on the Particle-Mesh Ewald (PME) and reaction field (RF) techniques with those obtained by an approach where the electrostatic interactions are truncated at r cut = 1.8, 2.0, and 2.5 nm. We find that the overall performance of PME is very good; its results are consistent with the expected behavior. The RF method performs rather well, too, despite certain inherent problems and the fact that its results differ from those obtained by PME. Nevertheless, the largest differences are found for the truncation methods, for which all examined truncation methods lead to results distinctly different from those obtained by PME. The lateral diffusion coefficients obtained by PME and truncation at 1.8 nm, for example, differ by a factor of 10, while the PME results are consistent with experimental values. The observed deviations can be interpreted in terms of artificial ordering due to truncation and highlight the important role of electrostatic interactions in the dynamics of systems composed of lipids and other biologically relevant molecules such as proteins and DNA.

Cationic lipid membranes are known to form compact complexes with DNA and to be effective as gene delivery agents both in vitro and in vivo. Here we employ molecular dynamics simulations for a detailed atomistic study of lipid bilayers consisting of a mixture of cationic dimyristoyltrimethylammonium propane (DMTAP) and zwitterionic dimyristoylphosphatidylcholine (DMPC). Our main objective is to examine how the composition of the DMPC/DMTAP bilayers affects their structural and electrostatic properties in the liquid-crystalline phase. By varying the mole fraction of DMTAP, we have found that the area per lipid has a pronounced nonmonotonic dependence on the DMTAP concentration, with a minimum around the point of equimolar DMPC/DMTAP mixture. We show that this behavior has an electrostatic origin and is driven by the interplay between positively charged TAP headgroups and the zwitterionic phosphatidylcholine (PC) heads. This interplay leads to considerable reorientation of PC headgroups for an increasing DMTAP concentration, and gives rise to major changes in the electrostatic properties of the lipid bilayer, including a significant increase of total dipole potential across the bilayer and prominent changes in the ordering of water in the vicinity of the membrane. Moreover, chloride counterions are bound mostly to PC nitrogens implying stronger screening of PC heads by Cl ions compared to TAP headgroups. The implications of these findings are briefly discussed.

In this paper we compare different force fields that are widely used (Gromacs, Charmm-22/xPlor, Charmm-27, Amber-1999, OPLS-AA) in biophysical simulations containing aqueous NaCl. We show that the uncertainties of the microscopic parameters of, in particular, sodium and, to a lesser extent, chloride translate into large differences in the computed radial-distribution functions. This uncertainty reflects the incomplete experimental knowledge of the structural properties of ionic aqueous solutions at finite molarity. We discuss possible implications on the computation of potential of mean force and effective potentials.

We construct a coarse-grained (CG) model for dipalmitoylphosphatidylcholine (DPPC) / cholesterol bilayers and apply it to large-scale simulation studies of lipid membranes. Our CG model is a two-dimensional representation of the membrane, where the individual lipid and sterol molecules are described by point-like particles. The effective intermolecular interactions used in the model are systematically derived from detailed atomic-scale molecular dynamics simulations using the Inverse Monte Carlo technique, which guarantees that the radial distribution properties of the CG model are consistent with those given by the corresponding atomistic system. We find that the coarse-grained model for the DPPC / cholesterol bilayer is substantially more efficient than atomistic models, providing a speed-up of approximately eight orders of magnitude. The results are in favor of formation of cholesterol-rich and cholesterol-poor domains at intermediate cholesterol concentrations, in agreement with the experimental phase diagram of the system. We also explore the limits of the novel coarse-grained model, and discuss the general validity and applicability of the present approach.

Free volume pockets or voids are important to many biological processes in cell membranes. Free volume fluctuations are a prerequisite for diffusion of lipids and other macromolecules in lipid bilayers. Permeation of small solutes across a membrane, as well as diffusion of solutes in the membrane interior are further examples of phenomena where voids and their properties play a central role. Cholesterol has been suggested to change the structure and function of membranes by altering their free volume properties. We study the effect of cholesterol on the properties of voids in dipalmitoylphosphatidylcholine (DPPC) bilayers by means of atomistic molecular dynamics simulations. We find that an increasing cholesterol concentration reduces the total amount of free volume in a bilayer. The effect of cholesterol on individual voids is most prominent in the region where the steroid ring structures of cholesterol molecules are located. Here a growing cholesterol content reduces the number of voids, completely removing voids of the size of a cholesterol molecule. The voids also become more elongated. The broad orientational distribution of voids observed in pure DPPC is, with a 30% molar concentration of cholesterol, replaced by a distribution where orientation along the bilayer normal is favored. Our results suggest that instead of being uniformly distributed to the whole bilayer, these effects are localized to the close vicinity of cholesterol molecules.