Re: [AMBER] tleap splitting the lipid

From: Matias Machado <mmachado.pasteur.edu.uy>
Date: Mon, 08 Jun 2020 10:12:34 -0300 (UYT)

Dear Zhiyi Wu,

The issue and the solution is documented in this post:

http://archive.ambermd.org/202001/0172.html

Best,

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]

----- Mensaje original -----
De: "David A Case" <david.case.rutgers.edu>
Para: "AMBER Mailing List" <amber.ambermd.org>
Enviados: Lunes, 8 de Junio 2020 9:22:36
Asunto: Re: [AMBER] tleap splitting the lipid

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
>molecular.

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.

....dac


_______________________________________________
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 Mon Jun 08 2020 - 06:30:03 PDT
Custom Search