Re: [AMBER] tleap splitting the lipid

From: David A Case <>
Date: Mon, 8 Jun 2020 08:22:36 -0400

On Mon, Jun 08, 2020, zhiyi wu wrote:
>>>source leaprc.protein.ff14SB
>>>source leaprc.lipid17
>>>lipid = sequence {MY PC MY}
>>>savepdb lipid lipid.pdb
>However, I notice that instead of having a ter at the end of the
>molecule, there is a ter at the end of the first acryl chain (MY) another
>ter at the end of the head group (PC) and a final ter at the end of the

The simple workaround is to not use tleap to create the PDB file.
Instead, use ambpdb with the parm7 and rst7 files that tleap creates.
You can also use the writeCoordinates action in parmed. Both of these
examine the bond list to find where molecules begin and end.


AMBER mailing list
Received on Mon Jun 08 2020 - 05:30:02 PDT
Custom Search