Hi!<br>&nbsp;I have a question . I am working with a protein that they will count two<br>centers FES co-ordinated for cisteins, these cisteins were not on to iron<br>atoms, thus I modified the archive itp only with respect the bonds, no
<br>angles or diedrals, to bind to these two residuos. After the minimização, I<br>made one md with restricted position, when I make the total dynamic with I<br>tie occurs error in the algorithm shake and it stops in the first steps
<br>informing that I did not obtain to set molecules of water. I am using the<br>field of force of gromacs (option 4) and the water is spc216.<br>&nbsp;ESS<br>