There is also intra_rms and instra_fit commands:
https://pymol.org/pymol-command-ref.html#intra_rms
https://pymol.org/pymol-command-ref.html#intra_fit
For example:
intra_fit name CA and resid 1326-1350, 0
Maybe they will be useful.
Best, Mat
On Mon, 19 Aug 2019 at 13:54, sunyeping <sunyep...@aliyun.com> wrote:
> Hello Mateusz,
>
> Thank you very much for your reply to my inquiry. Before I try the PyMOL
> extension you mentioned, I want to try scripting language in pymol. I
> wrote a attemptive one (align.py):
>
> cmd.load("A.gro")
> cmd.load("A.xtc")
> cmd.load("ref.pdb")
> for i in range(1, 5001):
> cmd.align(ref, A,target_state=i)
>
> I ran it in pymol command line, but it doesn't work. I am new to pymol
> scripting. Could you write a correct one for me? Thank you very much.
>
> Best regards,
>
> Yeping
>
>
> ------------------------------------------------------------------
> From:Mateusz Bieniek <bieniek...@gmail.com>
> Sent At:2019 Aug. 19 (Mon.) 20:32
> To:孙业平 <sunyep...@aliyun.com>
> Cc:pymol-users <pymol-users@lists.sourceforge.net>
> Subject:Re: [PyMOL] How to calculate rmsd between each state of a MD
> simulation trajectory and a reference structure?
>
> Hi Sunyeping,
>
> Besides using the scripting language in pymol, you might want to consider
> the experimental PyMOL extension created by the PyMOL Fellows (me and Paul
> Smith, https://pymol.org/fellowship/). The extension relies on the
> MDAnalysis package and so far includes RMSD as an example. Here is the link
> to the right github branch
> https://github.com/bieniekmateusz/pymol-open-source/tree/fellows_mp_2018
>
> If you compile it, you can use the RMSD this way:
> mda_load systemX.gro
> mda_load_traj systemX.xtc
> mda_rmsd
>
>
> Kind Regards,
> Mateusz
>
>
> On Mon, 19 Aug 2019 at 12:53, sunyeping via PyMOL-users <
> pymol-users@lists.sourceforge.net> wrote:
> Dear everyall,
>
> I loaded a molecular dynamics simulation trajectory (A.xtc) of 5000 frames
> and a reference structure (B.pdb) into pymol. I wish to get the rmsd value
> between each state of the MD simulation trajectory and the
> reference structure, but I don't know how. I can use the align command as:
> "align A, B", but this only gives one rmsd value. Could you tell me how to
> get the rmsd value between each state of the MD simulation trajectory and
> the reference structure?
> _______________________________________________
> PyMOL-users mailing list
> Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net
> Unsubscribe:
> https://sourceforge.net/projects/pymol/lists/pymol-users/unsubscribe
>
>
>
_______________________________________________
PyMOL-users mailing list
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net
Unsubscribe:
https://sourceforge.net/projects/pymol/lists/pymol-users/unsubscribe
