Dear Amber developers,
We've just found a possible "harmless" bug in Lipid11-17 libraries causing tleap to split lipid's head and tail fragments by TER statements when saving a PDB file.
The issue was previously reported in these threads:
+ [
http://archive.ambermd.org/201905/0138.html] ... [
http://archive.ambermd.org/201905/0197.html]
+ [
http://archive.ambermd.org/201610/0143.html]
A simple way to show the issue is executing the following lines in tleap (the same is obtained by loading a PDB):
source leaprc.lipid17
popc = sequence {PA PC OL}
savepdb popc popc.pdb
saveamberparm popc popc.prmtop popc.rst7
quit
Head and tail fragments of POPC (PA, PC and OL) are split by TER statements in the output "popc.pdb" file despite belonging to the same chain. The connectivity is correctly set in popc.prmtop.
I've found that the field "!entry.*.unit.residueconnect" (where * = any fragment) in lipid17.lib was set to " 0 0 0 0 0 ", the same happens in lipid11.lib and lipid14.lib.
Setting this entry to the same atom index in "!entry.*.unit.connect" fix the problem:
!entry.OL.unit.residueconnect table int c1x int c2x int c3x int c4x int c5x int c6x
1 1 0 0 0 0
!entry.PA.unit.residueconnect table int c1x int c2x int c3x int c4x int c5x int c6x
44 44 0 0 0 0
!entry.PC.unit.residueconnect table int c1x int c2x int c3x int c4x int c5x int c6x
1 37 0 0 0 0
It's not very clear to me the practical function of that entry, but the fix modifies the TREE_CHAIN_CLASSIFICATION section of the prmtop from "BLA BLA BLA BLA..." to something that makes more sense "3 E E E 3 E E ... M ... B ..."
Although this isn't a dramatic issue, neither the simulation nor the analysis with cpptraj is affected, it may constitute a pitfall in the AMBER's pipeline, as the output PDB file from leap can't be properly used again as a new input for leap.
Best regards,
Matias Machado
------------------------------------
PhD.
Researcher at Biomolecular Simulations Lab.
Institut Pasteur de Montevideo | Uruguay
[
http://pasteur.uy/en/labs/biomolecular-simulations-laboratory]
[
http://www.sirahff.com]
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jan 15 2020 - 09:00:02 PST
