Re: [AMBER] prmtop for stripped down coordinates

From: Jason Swails <jason.swails.gmail.com>
Date: Mon, 3 Jun 2013 09:18:16 -0400

On Mon, Jun 3, 2013 at 8:38 AM, Paul Driscoll <pdrisco.nimr.mrc.ac.uk>wrote:

> Hi,
>
> I plan to strip down my trajectory files to protein backbone amide NH
> atoms only. How do I prepare a prmtop for use in subsequent ptraj
> analysis? [In other words how do I prepare a prmtop for only a subset of
> protein atoms?]
>

As Marc suggested, parmstrip will do the trick (ParmEd also has a stripping
ability, but if you're already stripping a mask with cpptraj you might as
well use cpptraj to strip the topology, too).

I think it's even easier than that, though. In cpptraj, do your
stripping like this:

parm my.parm7
trajin mytraj.nc

strip !:N outprefix backbone_amide

trajout backbone_amide.nc

This will create a trajectory file 'backbone_amide.nc' and a matching
topology file backbone_amide.my.parm7.

However, I would be wary of stripping the coordinates down to _just_ a
single atom in each residue, since some actions in cpptraj make use of bond
information. If it was me, I would just strip out the solvent and leave
the entire solute intact and do any analyses you have using a specific atom
mask.

HTH,
Jason

-- 
Jason M. Swails
Quantum Theory Project,
University of Florida
Ph.D. Candidate
352-392-4032
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jun 03 2013 - 06:30:02 PDT
Custom Search