Re: [AMBER] Lipid14 falling apart

From: Benjamin Madej <bdmadej.gmail.com>
Date: Thu, 20 Mar 2014 14:19:42 -0700

Hi Filip,

We have done a lot of simulation of the six lipid types included in Lipid14
without the constant surface tension term. In particular, we were
interested in running with the tensionless, anisotropic NPT ensemble for
lipid bilayers.

We have been also working on membrane-bound protein systems. My
collaborators may be able to shed more light on what ensemble they have
used for their lipid membrane and protein systems and whether that ensemble
was appropriate.

All the best,
Ben Madej
Walker Molecular Dynamics Lab


On Thu, Mar 20, 2014 at 11:38 AM, filip fratev <filipfratev.yahoo.com>wrote:

> Hi Mustafa,
>
> Recently I performed several simulations with lipid14 POPC and did not
> detect such problem thus probably this is something related to your input
> file.
>
>
> Benjamin,
> Did you try Lipid14 in protein environment and whether the surface tension
> is still required as in the case of Lipid11?
>
>
> Regards,
> Filip
>
>
> ________________________________
> From: Benjamin Madej <bdmadej.gmail.com>
> To: AMBER Mailing List <amber.ambermd.org>
> Sent: Thursday, March 20, 2014 8:13 PM
> Subject: Re: [AMBER] Lipid14 falling apart
>
>
> Hi Mustafa,
>
> Yes, it sounds like the prmtop/inpcrd did not correctly link the lipid
> "residues" together into a lipid molecule.
>
> Verify that you are sourcing your force field files correctly, and read the
> leap.log file carefully. It also could be an issue with the residue and
> atom naming and TER card format in your input file. Double-check that your
> PDB file follows the Lipid11/Lipid14 format as described in the AmberTools
> manual and the Lipid11 paper.
>
> The warnings about "name change in pdb file residue" are normal and to be
> expected if you're using the charmmlipid2amber.x script to format your
> input PDB.
>
> It also appears that the residue and molecule types are not being found in
> the name map. Can you confirm that you are correctly sourcing the
> leaprc.lipid14 and leaprc.ff12SB files?
>
> When you say "there is TER flag introduced by tleap", how did you generate
> a PDB? The only way that I know of that correctly inserts TER cards for the
> conversion of Lipid11/Lipid14 prmtop/inpcrd to PDBs is the program
> "ambpdb". See the AmberTools manual for more details about ambpdb. If in a
> PDB generated by ambpdb, you see TER cards after each residue, then the
> connectivity is wrong.
>
> If all of that doesn't work, we'll need more information about your input
> PDB file format as well as the leap.log file to diagnose your problem.
>
> Finally, I strongly recommend that you use Amber for the actual production
> of your system. Lipid14 is fully compatible with GPU accelerated Amber and
> we ran all of our validation calculations with GPU accelerated Amber13.
>
> All the best,
> Ben Madej
> Walker Molecular Dynamics Lab
>
>
> On Thu, Mar 20, 2014 at 8:50 AM, Ghulam Mustafa <ghulam.mustafa.h-its.org
> >wrote:
>
> > Dear All
> >
> > I am trying to use new amber lipid14 ff on my protein-membrane complex
> > system but i get strange results after running heating and equilibration
> > (1.5ns) step. The head group of POPC lipid molecules fall apart from the
> > tails and keep floating in water box above and below the z axes of
> > bilayer. I have run initial minimization using Amber and then switching
> > to NAMD using amber forcefields (ff12SB and LIPID14).
> >
> > Although i find some warning in tleap while loading lipid-protein
> > complex but i can successfully prepare coordinate and topology files.
> > These warnings are as below:
> >
> > Enter zPdbReadScan from call depth 0.
> > Warning: name change in pdb file residue 492 ;
> > this residue is split into OL and PC.
> > Warning: name change in pdb file residue 492 ;
> > this residue is split into PC and PA.
> > Warning: name change in pdb file residue 493 ;
> > this residue is split into OL and PC.
> > Warning: name change in pdb file residue 493 ;
> > this residue is split into PC and PA.
> > Warning: name change in pdb file residue 494 ;
> > so on ..........
> > and for proteins:
> >
> > (Residue 1: ASP, Nonterminal, was not found in name map.)
> > (Residue 2: SER, Nonterminal, was not found in name map.)
> > (Residue 3: LEU, Nonterminal, was not found in name map.)
> > (Residue 4: VAL, Nonterminal, was not found in name map.)
> > (Residue 5: VAL, Nonterminal, was not found in name map.)
> > (Residue 6: LEU, Nonterminal, was not found in name map.)
> > (Residue 7: VAL, Nonterminal, was not found in name map.)
> >
> > similarly for ions
> >
> > (Residue 881: Na+, Nonterminal, was not found in name map.)
> > (Residue 882: Cl-, Nonterminal, was not found in name map.)
> > (Residue 883: Cl-, Nonterminal, was not found in name map.)
> > (Residue 884: Cl-, Nonterminal, was not found in name map.)
> > (Residue 885: Na+, Nonterminal, was not found in name map.)
> > (Residue 886: Na+, Nonterminal, was not found in name map.)
> > (Residue 887: Na+, Nonterminal, was not found in name map.)
> > (Residue 888: Na+, Nonterminal, was not found in name map.)
> > (Residue 889: Na+, Nonterminal, was not found in name map.)
> >
> > for solvent
> >
> > (Residue 1461: WAT, Nonterminal, was not found in name map.)
> > (Residue 1462: WAT, Nonterminal, was not found in name map.)
> > (Residue 1463: WAT, Nonterminal, was not found in name map.)
> > (Residue 1464: WAT, Nonterminal, was not found in name map.)
> > (Residue 1465: WAT, Nonterminal, was not found in name map.)
> > (Residue 1466: WAT, Nonterminal, was not found in name map.)
> > (Residue 1467: WAT, Nonterminal, was not found in name map.)
> > (Residue 1468: WAT, Nonterminal, was not found in name map.)
> > (Residue 1469: WAT, Nonterminal, was not found in name map.)
> > (Residue 1470: WAT, Nonterminal, was not found in name map.)
> >
> > and for lipids :
> >
> > (Residue 79964: PC, Terminal/last, was not found in name map.)
> > (Residue 79965: PA, Terminal/last, was not found in name map.)
> > (Residue 79966: OL, Terminal/last, was not found in name map.)
> > (Residue 79967: PC, Terminal/last, was not found in name map.)
> > (Residue 79968: PA, Terminal/last, was not found in name map.)
> > (Residue 79969: OL, Terminal/last, was not found in name map.)
> > (Residue 79970: PC, Terminal/last, was not found in name map.)
> > (Residue 79971: PA, Terminal/last, was not found in name map.)
> > (Residue 79972: OL, Terminal/last, was not found in name map.)
> > (Residue 79973: PC, Terminal/last, was not found in name map.)
> > (Residue 79974: PA, Terminal/last, was not found in name map.)
> > (Residue 79975: OL, Terminal/last, was not found in name map.)
> > (Residue 79976: PC, Terminal/last, was not found in name map.)
> > (Residue 79977: PA, Terminal/last, was not found in name map.)
> >
> > Is there something wrong with my initial coordinate and topology files?
> > as it seems there is no connectivity between head and tail groups, and
> > also there is TER flag introduced by tleap between each tail, head and
> > tail molecules.
> >
> > Looking for your kind input in this regard
> >
> > Mustafa
> >
> >
> > _______________________________________________
> > 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
> _______________________________________________
> 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 Mar 20 2014 - 14:30:03 PDT
Custom Search