For each simulations, the last frame was used as input for 300 ps NPT MDs (1 atm, 300 K, 1 fs timestep) with the AMOEBA polarizable forcefield as implemented in TINKER 8. We used the AOMEBA parameters reported by Chu et al 108 for the POPE lipids. We used the Nose-Hoover thermostat and barostat, as well as the Beeman integration method.