Hi Ben,
Thank you for your comment:
The error message indicates that one of your residues in your PDB has an
incorrect number of atoms. This is one of
the checks that the script does on the molecules while processing it.
Basically, the script does a substitution on each
atom in the residue and if it finds the incorrect number of atoms then it
stops.
It well correlates with the error message from charmmlipid2amber.py:
Error: Number of atoms in residue does not match number of atoms in residue
in replacement data file
Unfortunately your comment and the error message are too general to work
the problem around.
I can't find any specific place of incorrect atom number in the
step5_assembly.pdb file that causes the error, because the atom numbers are
smoothly increasing from 1 to 91154 as generated by Charmm-Gui.
Regarding the recommended in the A16 tutorial TER insertion after protein.
In original step5_assembly.pdb generated by Charmm-Gui there is no TER
after protein:
ATOM 5461 OT1 ALA 348 -10.030 -0.328 -36.334 1.00 0.00 PROA
ATOM 5462 OT2 ALA 348 -9.337 -1.565 -38.055 1.00 0.00 PROA
ATOM 5463 C3 CHL1 1 -28.933 18.011 17.293 1.00 0.00 MEMB
ATOM 5464 H3 CHL1 1 -27.890 18.323 17.067 1.00 0.00 MEMB
Following the A16 tutorial direction I have added it manually by text
editor:
ATOM 5461 OT1 ALA 348 -10.030 -0.328 -36.334 1.00 0.00 PROA
ATOM 5462 OT2 ALA 348 -9.337 -1.565 -38.055 1.00 0.00 PROA
TER
ATOM 5463 C3 CHL1 1 -28.933 18.011 17.293 1.00 0.00 MEMB
ATOM 5464 H3 CHL1 1 -27.890 18.323 17.067 1.00 0.00 MEMB
You see that indeed there is an interruption in atom numbers in the TER row
which doesn't contain any atom number in contrast to the standard PDB
format where
it shoud be: TER 5463 ALA 348
Unfortunately, I can't manually renumber all remaining ~90000 atoms after
TER by text editor,
so I used Chimera software to correct the situation:
ATOM 5461 OT1 ALA 348 -10.030 -0.328 -36.334 1.00 0.00 O
ATOM 5462 OT2 ALA 348 -9.337 -1.565 -38.055 1.00 0.00 O
TER 5463 ALA 348
HETATM 5464 C3 CHL1 1 -28.933 18.011 17.293 1.00 0.00 C
HETATM 5465 H3 CHL1 1 -27.890 18.323 17.067 1.00 0.00 H
As you can see the Chimera changes the original format of Charmm-Gui,
that could be a potential error origin.
I have examined all three variants of PDB file. The answer from
charmmlipid2amber.py was absolutely identical:
Error:
Number of atoms in residue does not match number of atoms in residue
in replacement data file
I would appreciate any your help in the problem solution.
Thank you,
Michael
*****************************
Michael Shokhen, PhD
Associate Professor
Department of Chemistry
Bar Ilan University,
Ramat Gan, 52900
Israel
email: michael.shokhen.gmail.com
email: shokhen.mail.biu.ac.il
On Thu, Apr 9, 2015 at 6:15 AM, Benjamin D Madej <bmadej.ucsd.edu> wrote:
> Hi Michael,
>
> The error message indicates that one of your residues in your PDB has an
> incorrect number of atoms. This is one of the checks that the script does
> on the molecules while processing it. Basically, the script does a
> substitution on each atom in the residue and if it finds the incorrect
> number of atoms then it stops.
>
> Residues are split in the PDB by the residue number field of the PDB and
> identified by the residue sequence field and TER cards. The script only
> processes C36 lipids, water, and some ions. All other residues are ignored.
>
> So, first check if the residues to be converted (lipids, water, ions)
> match the template of substitution file (charmmlipid2amber.csv). If there
> is no problem there, then it is probably a problem with residue sequence or
> TER. The script assumes that the normal PDB residue sequence and TER cards
> are used.
>
> Best,
> Ben
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Apr 09 2015 - 03:30:03 PDT
