Calculate Rmsd Between Two Sets Of Residues In Pymol
4
0
Entering edit mode
7.5 years ago
GreenDiamond ▴ 70

I have two peptides in PyMol, and each has one alpha helix that spans several residues. Is it possible to calculate the RMSD just between the two alpha helices?

protein-structure pymol • 22k views
ADD COMMENT
2
Entering edit mode
7.5 years ago
B. Arman Aksoy ★ 1.2k

I think you just need to run the align command with the appropriate parameters (to set a cutoff) and the syntax (to define residues for the alignment base). This tutorial gives a lot of details on these subjects: http://pldserver1.biochem.queensu.ca/~rlc/work/teaching/BCHM823/pymol/alignment/ and although pretty limited in scope, this wiki page might also come handy: http://www.pymolwiki.org/index.php/Align

For example in your case, all you need is to set the cutoff=1; e.g.:

PyMOL> align 1kao & i. 100-200, 1ctq & i. 100-200, cutoff=1

and you will see the RMSD score in the output printed within the console.

ADD COMMENT
2
Entering edit mode
7.5 years ago
Pappu ★ 1.9k

align command is not so good. There are a few other commands of this purpose: cealign, super, tmalign etc. http://www.pymolwiki.org/index.php/Tmalign

ADD COMMENT
1
Entering edit mode
7.5 years ago
João Rodrigues ★ 2.5k

Align, Super, CEalign, etc, all these command take selections so you can use for example align protein_A and ss h, protein_B and ss h to align the helical regions of each proteins. Be aware that PyMOL will try to optimize the fit (by removing atoms) and reduce the RMSD value. You can tweak this using the cycles option (not sure if available in all alignment commands).

ADD COMMENT
0
Entering edit mode
4.3 years ago
crooveck ▴ 10

If you only want to calculate RMSD without making any structure alignment or fitting, you should use rms_cur function.

ADD COMMENT

Login before adding your answer.

Traffic: 2762 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6