I'm looking for the easiest way of calculating RMSD between two homologous proteins. (i.e. sequences do not match perfectly). Until now I was using Chimera software, where steps are as follow:
- Calculating pairwise alignment
- Then structures are superimposed
- RMSD is calculated between corresponding amino acids in alignment.
And this work quite well, but I'm looking for something simpler that I can use in python script. Preferably in with no need of running subprocess.