The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

We present a new dual-resolution approach for coupling atomistic and coarse-grained models in molecular dynamics simulations of hydrated systems. In particular, a coarse-grained point dipolar water model is used to solvate molecules represented with standard all-atom force fields. A unique characteristic of our methodology is that the mixing of resolutions is direct, meaning that no additional or ad hoc scaling factors, intermediate regions, or extra sites are required. To validate the methodology, we compute the hydration free energy of 14 atomistic small molecules (analogs of amino acid side chains) solvated by the coarse-grained water. Remarkably, our predictions reproduce the experimental data as accurately as the predictions from state-of-the-art fully atomistic simulations. We also show that the hydration free energy of the coarse-grained water itself is in comparable or better agreement with the experimental value than the predictions from all but one of the most common multisite atomistic models. The coarse-grained water is then applied to solvate a typical atomistic protein containing both α-helix and β-strand elements. Moreover, parallel tempering simulations are performed to investigate the folding free energy landscape of a representative α helical and a β hairpin structure. For the simulations considered in this work, our dual-resolution method is found to be 3 to 6 times more computationally efficient than corresponding fully atomistic approaches.

Biological bilayer membranes typically contain varying amounts of lamellar and nonlamellar lipids. Lamellar lipids, such as dioleoylphosphatidylcholine (DOPC), are defined by their tendency to form the lamellar phase, ubiquitous in biology. Nonlamellar lipids, such as dioleoylphosphatidylethanolamine (DOPE), prefer instead to form nonlamellar phases, which are mostly nonbiological. However, nonlamellar lipids mix with lamellar lipids in biomembrane structures that remain overall lamellar. Importantly, changes in the lamellar vs nonlamellar lipid composition are believed to affect membrane function and modulate membrane proteins. In this work, we employ atomistic molecular dynamics simulations to quantify how a range of bilayer properties are altered by variations in the lamellar vs nonlamellar lipid composition. Specifically, we simulate five DOPC/DOPE bilayers at mixing ratios of 1/0, 3/1, 1/1, 1/3, and 0/1. We examine properties including lipid area and bilayer thickness, as well as the transmembrane profiles of electron density, lateral pressure, electric field, and dipole potential. While the bilayer structure is only marginally altered by lipid composition changes, dramatic effects are observed for the lateral pressure, electric field, and dipole potential profiles. Possible implications for membrane function are discussed.

Membrane permeation depends on fat content, and (permeant) size also matters.

The resistance of pathogens to traditional antibiotics is currently a global issue of enormous concern. As the discovery and development of new antibiotics become increasingly challenging, synthetic antimicrobial lipopeptides (AMLPs) are now receiving renewed attention as a new class of antimicrobial agents. In contrast to traditional antibiotics, AMLPs act by physically disrupting the cell membrane (rather than targeting specific proteins), thus reducing the risk of inducing bacterial resistance. In this study, we use microsecond-timescale atomistic molecular dynamics simulations to quantify the interaction of a short AMLP (C16-KKK) with model bacterial lipid bilayers. In particular, we investigate how fundamental transmembrane properties change in relation to a range of lipopeptide concentrations. A number of structural, mechanical, and dynamical features are found to be significantly altered in a non-linear fashion. At 10 mol% concentration, lipopeptides have a condensing effect on bacterial bilayers, characterized by a decrease in the area per lipid and an increase in the bilayer order. Higher AMLP concentrations of 25 and 40 mol% destabilize the membrane by disrupting the bilayer core structure, inducing membrane thinning and water leakage. Important transmembrane properties such as the lateral pressure and dipole potential profiles are also affected. Potential implications on membrane function and associated proteins are discussed.