I am working with a protein having phosphorylated serine residue. During topology generation in Gromacs I faced this problem with the above mentioned residue. Please suggest a forcefield, which have parameters for phosphoserine residue.