Residue Name for C-terminated Amide in CHARMM Force Fields

User discussions
1

GROMACS version: 2021.4
GROMACS modification: No

I am trying to perform MD simulations using CHARMM27 force field via GROMACS 2021.4.
The input structure file is a .pdb built from sequence using AmberTools20+tleap. The structure is c-terminated with NH2 and, by default, the “amino acid” is named as NHE (see last 5 lines of .pdb below).
When I attempt to generate GROMACS-required files for simulation, “NHE” residue is not found in the CHARM22 force field residue database. This issue isn’t seen in AMBER force fileds not OPLS as the topology file is buil

I tried changing residue name to “NH2”, the obvious answer, and got no results. I checked the force field rtp, itp and tdb files, and found no solution.
There might be a simple answer I am missing or not understanding. Perhaps it is the difference in CHARMM’s termini definition and patching, but it seems that amide termination is not there for CHARMM force fields.
Please let me know if there is a solution.

“”""
ATOM 138 N NHE 11 24.340 34.120 25.340 1.00 0.00
ATOM 139 H1 NHE 11 24.450 34.930 25.930 1.00 0.00
ATOM 140 H2 NHE 11 23.460 33.920 24.870 1.00 0.00
TER
ENDMDL
“”""

2

You do not need a special terminal residue for the case of an amidated C-terminus, you simply need to choose CT2 as the terminal patch when prompted by pdb2gmx. I will recommend that you upgrade to CHARMM36, unless you have some specific reason for using the old CHARMM22/CMAP parameters. The C36 port is available at MacKerell Lab

1 Like
3

Thank you so much Dr. Lemkul.
This isn’t the first time your forum replies have saved me from a headache and it is highly appreciated.