The modes of interaction of Ca 2+ and Mg 2+ with 11 pre-selected conformations of serine phosphate (SP) are investigated by using an additive procedure based on ab initio Self Consistent Field computations for the calculation of intermolecular interaction energies. Possible models for the arrangements, SP-Ca 2+-SP and SP-Mg 2+-SP, are investigated. The comparison between the binding energetics of Mg 2+ and Ca 2+ to one and two serine phosphates is discussed. It appears that some specific arrangements, SP-M 2+-SP ( M 2+ = Ca 2+ or Mg 2+ ), are able to account for the displayed marked selectivity of phosphatidylserine for Ca 2+, in keeping with the distinctive features of this complex in model membranes.