pymol加载两个结构,align一下。

rms_residue = cmd.rms_cur("chainblii and resid 16", "chainbsi and resid 16")

print("RMSD for residue 28:", rms_residue, "Å")