Abstract
Cellular lipid membranes are spatially inhomogeneous soft materials. Materials properties such as pressure and surface tension thus show important microscopic-scale variation that is critical to many biological functions. We present a means to calculate pressure and surface tension in a 3D-resolved manner within molecular-dynamics simulations and show how such measurements can yield important insight. We also present the first corrections to local virial and pressure fields to account for the constraints typically used in lipid simulations that otherwise cause problems in highly oriented systems such as bilayers. Based on simulations of an asymmetric bacterial ion channel in a POPC bilayer, we demonstrate how 3D-resolved pressure can probe for both short-range and long-range effects from the protein on the membrane environment. We also show how surface tension is a sensitive metric for inter-leaflet equilibrium and can be used to detect even subtle imbalances between bilayer leaflets in a membrane-protein simulation. Since surface tension is known to modulate the function of many proteins, this effect is an important consideration for predictions of ion channel function. We outline a strategy by which our local pressure measurements, which we make available within a version of the GROMACS simulation package, may be used to design optimally equilibrated membrane-protein simulations.