Hi, I saw that already. in yasara software, they assign the protonation states of amino acids and ligands depending on the chosen pH.i know that that classical MD does not include dissociation and reassociation of protons, so residues will not change their protonation states during the MD. possibility in gromacs ?