GROMACS version: 2020.1
GROMACS modification: Yes/No
Dear Gromacs users,
I downloaded the last version of the CHARMM36 ff from the MacKerell lab website and I intended to use this ff with pdb2gmx to obtain the topology of a pdb file that I have for a RNA sequence.
When running pdb2gmx it shows me the following error:
"atom N not found in building block 1ADE while combining tdb and rtp"
Now when checking the rtp for the FF I realized that there are not definition for 3’- and 5’-terminal nucleotides of a RNA sequence, or at least I can not find them. So it seems to me that the FF does not recognize the 1ADE (5’-adenosine) and the same will happen for the 3’-teminal CYT140. So the program can not find P or O1P or O2P, but this is my guess.
How could I correct this problem? Is there any way to add the 5- and 3’-terminal nucleotides to the .rtp file of the FF, or if they are there how to located them? What would be the name of these nucleotides in the rtp for CHARMM36?
As alternative solution I also generated the Gromacs compatible CHARMM36 topologies for this pdb using CHARMM-GUI. Is this a good enough solution?
Would appreciate any help with this.
Kindly,
Lazaro