Abstract Computational investigations into the structure and function of metalloenzymes with transition metal cofactors require proper preparation model, which requires obtaining reliable force field parameters for cofactor. Here, we present a test case where several methods were used to derive amber bonded model Fe(II) cofactor ectoine synthase. Moreover, spin ground state was probed by DFT po...