Challenges in the Accurate Modelling of Lipid Dynamics in Monolayers and Bilayers

Recent advances in hydrodynamic theory have revealed the severe effect of periodic boundary conditions (PBCs) on the diffusive dynamics of lipid membranes in molecular dynamics simulations. Even when accounting for PBC effects, the corrected lipid diffusion coefficients often severely overshoot the experimental estimates. Here, we investigate the underlying reasons for the exaggerated dynamics, and suggest potential ways for improvement. To this end, we examine the diffusion of four lipid types in both bilayers and monolayers using the CHARMM36 force field. We account for PBC effects using the full hydrodynamic treatment: for bilayers we use non-equilibrium simulations to extract the interleaflet friction parameter used in the correction; whereas monolayer hydrodynamics are treated by setting this parameter to zero. Our results suggest that the dynamics of bilayers are too fast, even if interleaflet friction is accounted for. However, the change of the water model to OPC leads to an excellent agreement with experiments. For monolayers, the dynamics with the TIP3P water model agree well with experiments, whereas they are undershot with OPC. As OPC and TIP3P differ in both shear viscosity and surface tension, we develop two new mass-scaled water models to clarify the roles of the thermodynamic and kinetic properties of the water model on lipid dynamics. Our results indicate that both of these quantities play a major role in lipid dynamics. Moreover, it seems that the accurate description of diffusion in both lipid bilayers and monolayers cannot be accounted for by changes in the water model alone, but likely also requires modifications in the lipid model.