On Mon, May 27, 2019, Hirdesh Kumar wrote:
>
>I have simulated a dimer (A: 1-370, B: 371-540)
>Now, I want to study the motion of 1 monomer with respect to the second
>monomer.
>
>parm mol.parm
>trajin traj.nc
>
>#align to first frame, backbone atoms
>reference :1:1-370.CA,N,C,O
>rms ref out rmsd.dat :1-370.CA,C,N,O
>rms ref out rmsd.dat :371-540.CA,C,N,O
You probably want to add the "nofit" keyword to the final rms command.
The first rms command will align monomer A to the reference. Then the
second command (adding in "nofit") will compute how far monomer B has
moved away from the reference state, without doing any fitting (since
alignment was already done with the first rms command.)
...good luck....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon May 27 2019 - 19:00:02 PDT
