GROMACS version: 2023.1 GROMACS modification: Yes/No Dear users, I have to deal with a simulation of a protein containing two zinc ions, tetrahedrally coordinated but with two different sets of ligands. In the first case, one zinc is bound by four cysteines residues - and I guess all of them should be in anioni form - whereas in the second case, one zinc is coordinated by two cy…