Y. C. Park, G. Fiorin, J. D. Faraldo-Gomez
A broad range of cellular functions involve transient or persistent changes in the morphology of lipid membranes, from the organellar to the molecular scale. By and large, the thermodynamics of protein-mediated remodeling, which often distort the membrane beyond its thermally accessible undulations, remain to be understood. Molecular Dynamics simulations enhanced by advanced sampling methods are uniquely suited to examine and quantitate these phenomena. Here, we leverage the Multi-Map simulation method to directly quantify how applied lateral tension impacts the energetics of both global and localized membrane perturbations induced extrinsically. We obtain free-energy profiles of curvature and thickness deformations with and without tension, which allow us to test the quadratic approximations of Helfrich-based continuum theories. We find that these continuum models are applicable up to ~ 20 kBT, beyond which they under-estimate the free-energy cost of bending. However, we find that the effect of tension on these deformations is primarily driven by changes in the projected area of the membrane, and therefore consistent with continuum theory, even beyond 30 kBT. These changes outpace the effects of tension on lipid packing and diffusion, and importantly we find that tension levels preceding rupture ([≤] 10 mN/m) do not impede complete lipid mixing. Finally, we investigate the free-energy landscapes of membrane thickness and lipid enrichment with a heterogeneous lipid composition and show that strong enrichment can be sustained even under steady exchange of the constituent lipid molecules.