First, one of the compromises we made for the presented lipid parameter set, was a small systematic decrease in APL across most lipid types in favor of drastically improved phase transition behavior. One axis requiring further exploration is the properties of the phosphate bead, for which we currently use the Q5 bead, a bead type assignment which is shared with other chemical groups such as carboxylates. It is possible that phosphate could be independently optimized to represent the distinct behavior of phospholipids. Additionally, fine-tuning the interactions between the phosphate and metal ions, like calcium 174 , could also be explored to improve the model's accuracy. Such parameterization could focus on accurate reproduction of experimental APL values while preserving the improved phase transition properties achieved with the current models. Additional refinements of the Martini 3 bead interaction matrix, specifically for the lipid headgroup region, are expected to further improve consistency with experimental and atomistic benchmarks.