Simulations of Asymmetric Membranes Illustrate Cooperative Leaflet Coupling and Lipid Adaptability

Abstract

Biological membranes are composed of lipid bilayers that are often asymmetric with regards to the lipid composition and/or aqueous solvent they separate. Studying lipid asymmetry both experimentally and computationally is challenging. Molecular dynamics simulations of lipid bilayers with asymmetry are difficult due to finite system sizes and time scales accessible to simulations. Due to the very slow flip-flop rate for phospholipids, one must first choose how many lipids are on each side of the bilayer, but the resulting bilayer may be unstable (or metastable) due to differing tensile and compressive forces between leaflets. Here we use molecular dynamics simulations to investigate a number of different asymmetric membrane systems, both with atomistic and coarse-grained models. Asymmetries studied include differences in number of lipids, lipid composition (unsaturated and saturated tails and different headgroups), and chemical gradients between the aqueous phases. Extensive analysis of the bilayers' properties such as area per lipid, density, and lateral pressure profiles are used to characterize bilayer asymmetry. We also address how cholesterol (which flip-flops relatively quickly) influences membrane asymmetries. Our results show how each leaflet is influenced by the other and can mitigate the structural changes to the bilayer overall structure. Cholesterol can respond to changes in bilayer asymmetry to alleviate some of the effect on the bilayer structure, but that will alter its leaflet distribution, which in turn affects its chemical potential. Ionic imbalances are shown to have a modest change in bilayer structure, despite large changes in the electrostatic potential. Bilayer asymmetry can also induce a modest electrostatic potential across the membrane. Our results highlight the importance of membrane asymmetry on bilayer properties, the influence of lipid headgroups, tails and cholesterol on asymmetry, and the ability of lipids to adapt to different environments.

Two Liberibacter Proteins Combine to Suppress Critical Innate Immune Defenses in Citrus

“…The Coulomb interactions were calculated using the smooth particle–mesh Ewald method (Cerutti et al, 2009 ; Kratz et al, 2016 ; Boateng, 2020 ) with a Fourier grid spacing of 0.12 nm (Fischer et al, 2015 ). The simulation in the NpT ensemble was achieved by semi-isotropic coupling at 1 bar with coupling constants of 4 ps (Aoki and Yonezawa, 1992 ; Blumer et al, 2020 ) and temperature-coupling the simulation system using velocity Langevin dynamics with a coupling constant of 1 ps (Washio et al, 2018 ). The integration time step was 2 fs.…”

How stereochemistry of lipid components can affect lipid organization and the route of liposome internalization into cells

“…The number of lipid molecules (128 total) was chosen based on recent reports. [65][66][67] Starting coordinates of the lipid bilayer, composed of 128 molecules of DMPC, were taken from the website of Biocomputing laboratory at the University of Calgary (http:// wcm.ucalgary.ca/tieleman/downloads) and solvated with 5227 water molecules in a rectangular box of 6.5 × 6.5 × 7.0 nm 3 . The mixed bilayers of DMPC and gemini components were obtained by a random selection of 52 molecules of DMPC replaced by the same number of gemini molecules.…”

Model Membrane Systems Used to Study Plasma Membrane Lipid Asymmetry

“…Published results from MD simulations of asymmetric bilayers include investigations of the lateral organization of lipids within the two leaflets of asymmetric bilayers [79,80,84,85], the effects of asymmetry on dipole potential [5,21,86], cholesterol distribution [87,88], permeation [89], peptide-induced membrane deformation [39], and interleaflet coupling [79,[90][91][92][93], to name a few. The properties of any simulated bilayer are subject to the underlying set of force field parameters, which determine the forces exerted on the individual atoms and thus the dynamics of the system.…”