The CHARMM force field for DMPC lipids was modified in order to improve agreement with experiment for a number of important properties of hydrated lipid bilayer. The modification consists in introduction of a scaling factor 0.83 for 1-4 electrostatic interactions (between atoms separated by three covalent bonds), which provides correct transgauche ratio in the alkane tails, and recalculation of the headgroup charges on the basis of HF/6-311(d,p) ab-initio computations. Both rigid TIP3P and flexible SPC water models were used with the new lipid model, showing similar results. The new model in a 75 ns simulation has shown a correct value of the area per lipid at zero surface tension, as well as good agreement with the experiment for the electron density, structure factor, and order parameters, including those in the headgroup part of lipids.

A series of molecular dynamics simulations of dimyristoylphosphatidylcholine bilayers, with different levels of hydration and temperature, were performed to examine the influence of hydration on properties of lipid membranes. Structural and dynamical properties such as area per lipid, electron densities, order parameters for all CH bonds and water, diffusion, and reorientation autocorrelation functions were determined and were all found to be affected by changes in the hydration level. The simulations give an overall picture of the bilayer going to a more ordered state when the hydration level is reduced. Lipid headgroups were found to adopt an orientation more parallel to the membrane plane when the water content was decreased. Dynamical properties such as lipid diffusion and relaxation of reorientation time correlation functions were found to become slower with the removal of water. Our simulation results generally agree with experimental data in cases where such data are available. One important conclusion drawn is that while structural properties are affected only moderately dynamical properties are affected very strongly by a decrease of water content.

The influence of the local anesthetic lidocaine on electrostatic properties of a lipid membrane bilayer was studied by molecular dynamics simulations. The electrostatic dipole potential, charge densities, and orientations of the headgroup angle have been examined in the presence of different amounts of charged or uncharged forms of lidocaine. Important changes in the membrane properties caused by the presence of both forms of lidocaine are presented and discussed. Our simulations have shown that both charged and uncharged lidocaine cause almost the same increase in the electrostatic potential in the middle of the membrane, although for different reasons. The increase, approximately 90 mV for 9 mol % of lidocaine and 220 mV for 28 mol % of lidocaine, is of a size that may affect the functioning of voltage-gated ion channels.

No abstract

This paper describes an analysis of NMR dipolar couplings in a bilayer formed by dimyristoylphosphatidylcholine (DMPC). The couplings are calculated from a trajectory generated in a molecular dynamics (MD) simulation based on a realistic atom-atom interaction potential. The analysis is carried out employing a recently developed approach that focuses on the construction of the conformational distribution function. This approach is a combination of two models, the additive potential (AP) model and the maximum entropy (ME) method, and is therefore called APME. In contrast to the AP model, the APME procedure does not require an intuition-based choice of the functional form of the torsional potential and is, unlike the ME method, applicable to weakly ordered systems. The conformational distribution function for the glycerol moiety of the DMPC molecule derived from the APME analysis of the dipolar couplings is in reasonable agreement with the "true" distributions calculated from the trajectory. Analyses of dipolar couplings derived from MD trajectories can, in general, serve as guidelines for experimental investigations of bilayers and other complex biological systems.