Abstract

Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity ({kappa}C) and area compressibility (KA). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. While computational tools such as Molecular Dynamics (MD) simulations can provide alternative means to estimate KA values, current approaches suffer significant technical limitations. Here, we present a novel computational framework that allows for a direct estimation of KA values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives KA from real-space analysis of local thickness fluctuations sampled in MD simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid composition and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model (PBM) and confirm the PBMs predicted linear relationship with proportionality coefficient of 24 using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the new method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular, its dependence on acyl chain unsaturation and the presence of cholesterol.