We irst explained how to estimate the spontaneous curvature caused by asymmetric Ca 2+ concentration from the MD simulations. When two lipid monolayers were combined to form a bilayer, we used the mid-plane between the two monolayers to represent the shape of the bilayer, which was parameterized as X. This position was denoted as z = 0 along the normal direction of the bilayer.