Dear User, I would like to constraint the Iron bound active site of the protein in octahedral shape in EM and MD steps ... Is there a direct routine to create such a constraint in gromacs by feeding information of ligating atoms of amino acids with Iron at the center of the octahedral geometry. Thanks, Raja