[PyMOL] How to automate RMSD calculation for large no of structures

2014-01-07 Thread Om Shiv
Hello All I am a new user to Pymol and Python. I have a PDB Structure with single chain (A) and I have modeled 500 single point mutants of this wild structure. Now for each such 500 modeled structures, I would like to calculate RMSD (ALL ATOMS) with the wild type PDB structure. Can anybody help

Re: [PyMOL] How to automate RMSD calculation for large no of structures

2014-01-07 Thread Jed Goldstone
Check out Robert Campbell's Pymol script repository. I think align_all_to_all.py should do what you want, including exporting RMSD values. http://pldserver1.biochem.queensu.ca/~rlc/work/pymol/ Jed On 1/7/2014 11:58 AM, Om Shiv wrote: > Hello All > > I am a new user to Pymol and Python. > > I h

Re: [PyMOL] How to automate RMSD calculation for large no of structures

2014-01-07 Thread Robert Campbell
On Tue, 2014-01-07 12:05 EST, Jed Goldstone wrote: > Check out Robert Campbell's Pymol script repository. > I think align_all_to_all.py should do what you want, including exporting > RMSD values. > > http://pldserver1.biochem.queensu.ca/~rlc/work/pymol/ Actually, I think Om only wants to ali

Re: [PyMOL] How to automate RMSD calculation for large no of structures

2014-01-07 Thread Tsjerk Wassenaar
If the Wild-Type is called wildtype: rmsds = [(i,cmd.align(i,"wildtype")[0]) for i in cmd.get_object_list("not wildtype")] print rmsds or, writing out the results directly open("rmsds.dat","w").writelines("%s %f\n"%(i,cmd.align(i,"wildtype")[0]) for i in cmd.get_object_list("not wildtype")) Che

Re: [PyMOL] How to automate RMSD calculation for large no of structures

2014-01-07 Thread Thomas Holder
I like to add that to get the ALL ATOMS RMSD you need to pass cycles=0 to cmd.align or to Rob's align_allfiles. Otherwise, the outlier rejection refinement might drop atoms from the alignment. Cheers, Thomas On 07 Jan 2014, at 13:58, Tsjerk Wassenaar wrote: > If the Wild-Type is called wild

Re: [PyMOL] How to automate RMSD calculation for large no of structures

2014-01-07 Thread Sampson, Jared
Hi Om - Here is another option that loads only one mutant at a time. You can save the following as a Python script (e.g., "mut-rmsd.py") and run it from the PyMOL command line by typing `run mut-rmsd.py` (with a full or relative path if it's not in PyMOL's working directory). For more informa