pymol加载两个结构,align一下。
rms_residue = cmd.rms_cur("chainblii and resid 16", "chainbsi and resid 16")
print("RMSD for residue 28:", rms_residue, "Å")