Bending Modulus of Lipid Membranes from Density Correlation Functions

The bending modulus κ quantifies the elasticity of biological membranes in terms of the free energy cost of increasing the membrane corrugation. Molecular dynamics (MD) simulations provide a powerful approach to quantify κ by analyzing the thermal fluctuations of the lipid bilayer. However, existing methods require the identification and filtering of non-mesoscopic fluctuation modes. State of the art methods rely on identifying a smooth surface to describe the membrane shape. These methods introduce uncertainties in calculating κ since they rely on different criteria to select the relevant fluctuation modes. Here, we present a method to compute κ using molecular simulations. Our approach circumvents the need to define a mesoscopic surface or an orientation field for the lipid tails explicitly. The bending and tilt moduli can be extracted from the analysis of the density correlation function (DCF). The method introduced here builds on the Bedeaux and Weeks (BW) theory for the DCF of fluctuating interfaces and on the coupled undulatory (CU) mode introduced by us in previous work. We test the BW-DCF method by computing the elastic properties of lipid membranes with different system sizes (from 500 to 6000 lipid molecules) and using coarse-grained (for POPC and DPPC lipids) and fully atomistic models (for DPPC). Further, we quantify the impact of cholesterol on the bending modulus of DPPC bilayers. We compare our results with bending moduli obtained with X-ray diffraction data and different computer simulation methods.