Article ID: 637 - Last Modified: August 4, 2011
I am trying to calculate the RMSD between the docked pose of a ligand, and the same co-crystallized ligand. Unfortunately the atom numbers don't match, resulting in an error. Do I have to renumber the atom numbers or is there any other way of doing it?
Yes, you can use the Script Center script, rmsd.py, which you can obtain from our web site (Resources & Downloads → Script Center, in the General Tools section). This script computes the RMSD after automatically renumbering the atoms. It is a command-line script; if you install it via Scripts → Update in Maestro, you can run it as follows:
$SCHRODINGER/run rmsd.py file1 file2
file1 and file2 are the Maestro files that contain the co-crystallized ligand and the pose file. If you have generated ionization or tautomeric states of your ligand, you should use the
-use_neutral_scaffold option to ensure that the multiple bonds patterns and ionization states are ignored when calculating the RMSD.
Keywords: calculate RMSD, renumber atoms, rmsd.py
Type the words or phrases on which you would like to search, or click here to view a list of all
Knowledge Base articles