Parametrization of pyroglutamate

GROMACS version: 2026
GROMACS modification: No

Hi everyone,

I am currently setting up simulations of Aβ(pE3-42) using the CHARMM36m force field in GROMACS.

The N-terminal residue of this peptide is pyroglutamate (pGlu/PCA). Because this is a cyclized residue lacking a free N-terminal charge, it is not included in the standard CHARMM36m force field.

Could anyone advise on the best approach to parametrize and integrate pyroglutamate into CHARMM36m for use with GROMACS?

Any guidance, relevant literature, or suggestions would be greatly appreciated.

Thanks in advance for your help!

The latest Tinker distribution has charmm36m parameters for pGlu, convert it over to gromacs format