VMD-L Mailing List

From: Ashar Malik (asharjm_at_gmail.com)
Date: Fri Dec 15 2017 - 06:10:00 CST

For simplicity I am going to break this problem down.

Your starting frame for each of the two trajectories is the same. Lets call
this *"0"*.
Then your two trajectories can be represented by:

Trajectory - A: {0, 1A, 2A, 3A, ... }
Trajectory - B: {0, 1B, 2B, 3B, ... }

Lets say that trajectory *"A"* is with the constraints and *"B"* without
the constraints.

RMSD is calculated by comparing each frame in a trajectory to a reference
frame. If you didn't change this, the default should have been frame
*"0" *which
is the starting frame.
So the following comparisons will have been made.

Trajectory - 1(RMSD): {0-0, 1A-0, 2A-0, 3A-0, ... }
Trajectory - 2(RMSD): {0-0, 1B-0, 2B-0, 3B-0, ... }

If I understand correctly your RMSD plots start at non-zero values. Right?
Here are some situations.

   1. Do your plots start at frame 0? If they do it must be a value of 0
   for both plots. In which case the starting values should be identical for
   both plots because frame 0 is compared to itself in both plots.
   2. Do your plots start at frame 1? In which case why would want them to
   be the same? Frame 1A and 1B are not the same. So why should you expect the
   same RMSD? (This is what is likely happening).

There may be other reasons - e.g. your trajectory alignment method which
may also result in a value mismatch, but I think the previous two points
are important. If you provide more evidence and that doesn't explain the
difference you can write back and some one will help.

On Fri, Dec 15, 2017 at 11:55 PM, Monika Madhavi <monikamadhavi_at_gmail.com>
wrote:

> Hi all,
>
> I minimised and equilibrated a protein water system. Then I used the final
> configuration of this system to start two simulations with the same
> parameters except the constraints on backbone atoms. When I calculated the
> rmsd using Rmsd visualiser tool, the rmsd starts from two values for the
> two trajectories. I can't figure out what could be the reason for this. Any
> insight would be appreciated.
>
> Thank you
> Best,
> Monika
> 

-- 
Best,
/A