Bending Modulus of Lipid Membranes from Density Correlation Functions.

Area of Science:

• Biophysics
• Materials Science
• Computational Chemistry

Background:

• The bending modulus (κ) is crucial for understanding biological membrane elasticity and its response to deformation.
• Current methods using molecular dynamics (MD) simulations to calculate κ often rely on identifying mesoscopic surfaces, introducing uncertainties.
• Filtering non-mesoscopic fluctuation modes is a challenge in existing simulation-based approaches.

Purpose of the Study:

• To develop and validate a novel method for computing the bending modulus (κ) of lipid membranes using molecular simulations.
• To overcome limitations of existing methods that require explicit surface or orientation field definitions.
• To analyze the impact of system size, lipid models, and cholesterol on membrane bending rigidity.

Main Methods:

• Utilizing molecular dynamics (MD) simulations to analyze thermal fluctuations of lipid bilayers.
• Applying a novel method based on the Bedeaux and Weeks (BW) theory for the density correlation function (DCF).
• Leveraging the coupled undulatory (CU) mode analysis to extract elastic properties without defining mesoscopic surfaces.

Main Results:

• Successfully computed the bending modulus (κ) for lipid membranes of varying sizes (500–6000 lipids) using coarse-grained and all-atom models.
• Quantified the influence of cholesterol on the bending modulus of DPPC bilayers.
• Demonstrated that the BW-DCF method, combined with CU modes, effectively extracts elastic properties.

Conclusions:

• The proposed BW-DCF method provides a robust and direct approach to calculate membrane bending modulus (κ) from MD simulations.
• This method eliminates the need for complex surface identification, reducing uncertainties in κ calculations.
• The findings offer a more reliable computational tool for studying lipid membrane mechanics and the effects of composition.