Re: [AMBER] cpptraj: 1) align using a certain region 2) calc rmsd for another region

From: Jason Swails <jason.swails.gmail.com>
Date: Tue, 30 Aug 2011 11:04:18 -0400

On Tue, Aug 30, 2011 at 10:41 AM, Jan-Philip Gehrcke <
jgehrcke.googlemail.com> wrote:

> Hey.
>
> Using cpptraj, I am interested in how to get the RMSD of a certain
> region within two structures after fitting these structures using
> another region. The same question was already answered for ptraj:
> http://archive.ambermd.org/201102/0128.html
>
> Quote by Dan:
> "
> reference refstruct.rst7
> rms reference :1-100
> rms reference :50-55 nofit out rmsd.dat
>
> The first rms command performs the best fit, then the second rms
> command calculates the rmsd you want without changing the fit from the
> first rms command.
> "
>
> This is what I propose for cpptraj (for my system):
>
> parm topology.top
> trajin traj1.mdcrd
> rmsd :1-150 first
>

This aligns all structures by minimizing the RMS of that structure against
the first structure using only residues 1-150 when calculating the RMS.


> rmsd :151-161 first nofit out rmsd.dat
>

This calculates the RMS between each structure using residues 151-161
*without* first doing a fit to that region. Thus, the structures will still
be aligned according to the RMS command above (since that one is executed
first).


> Would that for each frame
> 1) align the residues 1-150 of the current frame to the residues 1-150
> of the first frame and then
> 2) calculate the RMSD for residues 151 to 161 with the alignment from
> (1) (again current frame vs. first frame)?
>

Yes.

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 Tue Aug 30 2011 - 08:30:02 PDT
Custom Search