GROMACS version: :-) GROMACS - gmx pdb2gmx, 2022.4
GROMACS modification: Yes/No
Hello,
I am trying to convert a pdb file that is a nano-particle generated from ChimeraX to a gro file and generate the topology. It consists of many span 60. The problem is Span 60 is not on Charmm36 that I am aware of, so it suggests generating the hdb and rtp files. I have not been able to locate a readily available file of these needed files for Span 60 and the generators based on information given by PubChem have helped me develop an rtp and hdb, which I added it to the forcefield as UNL in aminoacid.rtp and aminoacid.hbd. But I am getting this error:
All occupancies are one
All occupancies are one
Opening force field file ./charmm36-jul2022.ff/atomtypes.atp
Reading residue database… (Charmm36-jul2022)
Opening force field file ./charmm36-jul2022.ff/aminoacids.rtp
Warning: file does not end with a newline, last line:
-C N CA C +N
Opening force field file ./charmm36-jul2022.ff/carb.rtp
Opening force field file ./charmm36-jul2022.ff/cgenff.rtp
Opening force field file ./charmm36-jul2022.ff/ethers.rtp
Opening force field file ./charmm36-jul2022.ff/lipid.rtp
Opening force field file ./charmm36-jul2022.ff/metals.rtp
Opening force field file ./charmm36-jul2022.ff/na.rtp
Opening force field file ./charmm36-jul2022.ff/silicates.rtp
Opening force field file ./charmm36-jul2022.ff/solvent.rtp
Opening force field file ./charmm36-jul2022.ff/aminoacids.hdb
Opening force field file ./charmm36-jul2022.ff/carb.hdb
Opening force field file ./charmm36-jul2022.ff/cgenff.hdb
Opening force field file ./charmm36-jul2022.ff/ethers.hdb
Opening force field file ./charmm36-jul2022.ff/lipid.hdb
Opening force field file ./charmm36-jul2022.ff/metals.hdb
Opening force field file ./charmm36-jul2022.ff/na.hdb
Opening force field file ./charmm36-jul2022.ff/silicates.hdb
Opening force field file ./charmm36-jul2022.ff/solvent.hdb
Opening force field file ./charmm36-jul2022.ff/aminoacids.n.tdb
Opening force field file ./charmm36-jul2022.ff/carb.n.tdb
Opening force field file ./charmm36-jul2022.ff/cgenff.n.tdb
Opening force field file ./charmm36-jul2022.ff/ethers.n.tdb
Opening force field file ./charmm36-jul2022.ff/lipid.n.tdb
Opening force field file ./charmm36-jul2022.ff/metals.n.tdb
Opening force field file ./charmm36-jul2022.ff/na.n.tdb
Opening force field file ./charmm36-jul2022.ff/silicates.n.tdb
Opening force field file ./charmm36-jul2022.ff/solvent.n.tdb
Opening force field file ./charmm36-jul2022.ff/aminoacids.c.tdb
Opening force field file ./charmm36-jul2022.ff/carb.c.tdb
Opening force field file ./charmm36-jul2022.ff/cgenff.c.tdb
Opening force field file ./charmm36-jul2022.ff/ethers.c.tdb
Opening force field file ./charmm36-jul2022.ff/lipid.c.tdb
Opening force field file ./charmm36-jul2022.ff/metals.c.tdb
Opening force field file ./charmm36-jul2022.ff/na.c.tdb
Opening force field file ./charmm36-jul2022.ff/silicates.c.tdb
Opening force field file ./charmm36-jul2022.ff/solvent.c.tdb
Processing chain 1 ‘B’ (124992 atoms, 63 residues)
Problem with chain definition, or missing terminal residues. This chain does not appear to contain a recognized chain molecule. If this is incorrect, you can edit residuetypes.dat to modify the behavior.
8 out of 8 lines of specbond.dat converted successfully
Program: gmx pdb2gmx, version 2022.4
Source file: src/gromacs/gmxpreprocess/resall.cpp (line 617)
Fatal error:
Residue ‘UNL’ not found in residue topology database
and I am unsure if I should
1.change the HETAM type to ATOM?
2. Change the rtp and hdb information to another type file like ligand or carb?
3. Or if I need to add something on top of the other two files?
4. Or if there may be someone who has experience. This is my first time trying to create rtp and hdb, plus running a fully Nano(Span60) model in Gromacs.
Thanks in advance and appreciate any help,
Heather