Effect of methanol on the phase-transition properties of glycerol-monopalmitate lipid bilayers investigated using molecular dynamics simulations: In quest of the biphasic effect

Free Energy Cost of Interdigitation of Lamellar Bilayers of Fatty Alcohols with Cationic Surfactants from Molecular Dynamics Simulations

Cationic surfactant mixed with fatty alcohol as cosurfactant in excess water can form stable emulsions, known as "lamellar gel networks," that contain extended and interconnected networks of swollen bilayers, including ones with in-plane liquidlike disorder (L_α phase) and solidlike order (L_β phase). To study their structure and thermodynamics, molecular dynamics (MD) simulations with lateral pressure and temperature scans along reversible pathways were used to drive reversible phase changes, including formation at negative lateral pressure of the L_βI phase with interdigitated tails of opposing leaflets. Thermodynamic integration, with extrapolations to infinitely slow scans, yielded a free energy difference between the interdigitated L_βI and non-interdigitated L_β phases of 2.4 ± 0.5 kJ/mol, which is consistent with the spontaneous formation of the L_β phase under atmospheric pressure in simulation. Thermodynamic cycles involving temperature and lateral pressure for which the free energy difference is identically zero were constructed as negative controls to verify the method. Using lateral pressure, including negative lateral pressure, helps avoid kinetic bottlenecks that occur when temperature alone is used as the control variable. The method, using negative lateral pressure, should be widely applicable to other bilayers to identify molecular properties that control interdigitation and other bilayer properties.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

Long-timescale motions in glycerol-monopalmitate lipid bilayers investigated using molecular dynamics simulation

The occurrence of long-timescale motions in glycerol-1-monopalmitate (GMP) lipid bilayers is investigated based on previously reported 600 ns molecular dynamics simulations of a 2×8×8 GMP bilayer patch in the temperature range 302-338 K, performed at three different hydration levels, or in the presence of the cosolutes methanol or trehalose at three different concentrations. The types of long-timescale motions considered are: (i) the possible phase transitions; (ii) the precession of the relative collective tilt-angle of the two leaflets in the gel phase; (iii) the trans-gauche isomerization of the dihedral angles within the lipid aliphatic tails; and (iv) the flipping of single lipids across the two leaflets. The results provide a picture of GMP bilayers involving a rich spectrum of events occurring on a wide range of timescales, from the 100-ps range isomerization of single dihedral angles, via the 100-ns range of tilt precession motions, to the multi-μs range of phase transitions and lipid-flipping events.