Re: [AMBER] cpptraj - rmsd/rmsf & reference structure

From: colvin chin <colvin4367.gmail.com>
Date: Thu, 26 May 2016 21:56:03 +0800

Hi.

Thank you for the reply.

I copied the wrong output for case 1. For case 1, I used rms reference, and the output didn't state the reference mask as in the output for case 2 when I use rms first. I did specified reference structure before the rms reference command. Am I missing anything here?

Thanks!

> On 26 May 2016, at 9:22 PM, Daniel Roe <daniel.r.roe.gmail.com> wrote:
>
> Hi,
>
>> On Thu, May 26, 2016 at 3:19 AM, colvin <colvin4367.gmail.com> wrote:
>> Is this strange bug fixed?
>>
>> http://archive.ambermd.org/201509/0399.html
>>
>> I am still getting the same problem as described. Currently using CPPTRAJ ver. 15.
>
> That bug is indeed fixed in version 16 (and I recommend upgrading).
> However, your problem does not appear to be the same one.
>
>>
>> Does this means that in case 1, the rmsd calculation is not rmsfit to the reference structure? and in case 2, it is rmsfit into the first frame of the input trajectory instead of my reference structure - min2_worst.rst?
>
> In both cases you're specifying 'first', which means use the first
> frame read in as the reference structure. You need to specify
> 'reference' or 'ref <refname>' etc. if you want to use a specific
> reference structure. See the manual for full details.
>
>>
>> For the atomicfluct, if rms reference command put before atomicfluct, the calculation will rmsfit to the reference structure, right?
>
> Yes, that is what the 'rms' command does by default unless 'nofit' or
> 'nomod' is specified.
>
> Hope this helps,
>
> -Dan
>
>>
>> My input:
>> trajin ../combine_traj/combine-every-10ps-production.nc
>> center :1-400
>> image familiar
>> reference ../../min2_worst.rst
>> rms reference
>> atomicfluct out backbonermsfluct_all :1-400.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_all.arg :1-400.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_all.gnu :1-400.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_rec :1-394.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_rec.arg :1-394.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_rec.gnu :1-394.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_pep :395-400.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_pep.arg :395-400.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_pep.gnu :395-400.C,CA,N,O byres
>> atomicfluct out backbonermsfluct_activesite :395-400<:4.0 byres
>> go
>>
>> Pls advise.
>>
>> Thank you very much!!
>> _______________________________________________
>> AMBER mailing list
>> AMBER.ambermd.org
>> http://lists.ambermd.org/mailman/listinfo/amber
>
>
>
> --
> -------------------------
> Daniel R. Roe, PhD
> Department of Medicinal Chemistry
> University of Utah
> 30 South 2000 East, Room 307
> Salt Lake City, UT 84112-5820
> http://home.chpc.utah.edu/~cheatham/
> (801) 587-9652
> (801) 585-6208 (Fax)
>
> _______________________________________________
> 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 May 26 2016 - 07:00:03 PDT
Custom Search