Hi,
On Fri, Sep 20, 2019 at 3:25 PM Lenka <Vellryba.seznam.cz> wrote:
> I have stripped all the atoms from trajectories and pdb (and associated
> prmtops) so that there is an equal number of atoms.
>
> For a sanity check, I got a few frames from the trajectory and compared the
> rmsd using Chimera.
>
> My issues is, that amber calculates the rmsd to be around 7 and Chimera
> around 1. Since the structures are very similar, Chimera is right.
As mentioned previously, given the same input and reference
structures, cpptraj and chimera should return a best-fit RMSD within
numerical precision.
There are a few places your calculation could have gone awry - most
notably with exactly what atoms you stripped from each structure (you
don't say exactly what you did so no way to check if it's ok or not).
In order for the RMSD calculation to work properly, not only the
number of atoms needs to match but the atom ordering as well. Could
you send me those files as well? Preferably the originals too if you
have them. Also any other relevant information (like what atoms were
stripped). Thanks,
-Dan
>
>
> Any idea as to what am I doing wrong?
>
>
>
>
>
> Thank you very much
>
>
>
>
>
> L.
> _______________________________________________
> 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 Mon Sep 23 2019 - 06:30:05 PDT
