Hi, Pymol Users, Would anyone help me how to count the salt bridges in all models?
After my script, I have the salt bridge list for each model. I would like to know all the possible salt bridges in all models, and their frequency of each salt bridge in all the models. for example, /A/ASP/446/OD2 to /A/LYS/239/NZ: 3.776 exist in both models, while A/GLU/328/OE1 to/A/ARG/322/NE: 2.377 exist in only one model. I have total 50 models. I cannot image a manual way to isolate each salt bridge and count their frequency in 50 models. My wild guess is that text processing might solve this issue. Any suggestions? Thanks! The output of model00000 is as follows model00000/A/GLU/328/OE1 to model00000/A/ARG/322/NE: 2.377 model00000/A/GLU/430/OE2 to model00000/A/ARG/368/NE: 3.883 model00000/A/ASP/446/OD1 to model00000/A/LYS/239/NZ: 2.818 model00000/A/ASP/446/OD2 to model00000/A/LYS/239/NZ: 3.877 model00000/B/GLU/328/OE1 to model00000/B/ARG/322/NE: 2.376 model00000/B/GLU/430/OE2 to model00000/B/ARG/368/NE: 3.884 model00000/B/ASP/446/OD1 to model00000/B/LYS/239/NZ: 2.818 model00000/B/ASP/446/OD2 to model00000/B/LYS/239/NZ: 3.877 model00000/C/GLU/328/OE1 to model00000/C/ARG/322/NE: 2.377 model00000/C/GLU/430/OE2 to model00000/C/ARG/368/NE: 3.883 model00000/C/ASP/446/OD1 to model00000/C/LYS/239/NZ: 2.818 model00000/C/ASP/446/OD2 to model00000/C/LYS/239/NZ: 3.876 model00000/D/GLU/328/OE1 to model00000/D/ARG/322/NE: 2.377 model00000/D/GLU/430/OE2 to model00000/D/ARG/368/NE: 3.883 model00000/D/ASP/446/OD1 to model00000/D/LYS/239/NZ: 2.819 model00000/D/ASP/446/OD2 to model00000/D/LYS/239/NZ: 3.876 model00000/E/GLU/328/OE1 to model00000/E/ARG/322/NE: 2.378 model00000/E/GLU/430/OE2 to model00000/E/ARG/368/NE: 3.883 model00000/E/ASP/446/OD1 to model00000/E/LYS/239/NZ: 2.819 model00000/E/ASP/446/OD2 to model00000/E/LYS/239/NZ: 3.877 The output of model00001 is as follows model00001/A/ASP/429/OD1 to model00001/B/ARG/424/NE: 2.634 model00001/A/ASP/429/OD2 to model00001/B/ARG/424/NE: 3.672 model00001/A/GLU/437/OE2 to model00001/A/ARG/310/NE: 3.310 model00001/A/ASP/446/OD2 to model00001/A/LYS/239/NZ: 3.776 model00001/B/ASP/429/OD1 to model00001/C/ARG/424/NE: 2.600 model00001/B/ASP/429/OD2 to model00001/C/ARG/424/NE: 3.629 model00001/B/GLU/437/OE2 to model00001/B/ARG/310/NE: 3.311 model00001/B/ASP/446/OD2 to model00001/B/LYS/239/NZ: 3.777 model00001/C/ASP/429/OD1 to model00001/D/ARG/424/NE: 2.597 model00001/C/ASP/429/OD2 to model00001/D/ARG/424/NE: 3.621 model00001/C/GLU/437/OE2 to model00001/C/ARG/310/NE: 3.311 model00001/C/ASP/446/OD2 to model00001/C/LYS/239/NZ: 3.776 model00001/D/ASP/429/OD1 to model00001/E/ARG/424/NE: 2.574 model00001/D/ASP/429/OD2 to model00001/E/ARG/424/NE: 3.604 model00001/D/GLU/437/OE2 to model00001/D/ARG/310/NE: 3.310 model00001/D/ASP/446/OD2 to model00001/D/LYS/239/NZ: 3.776 model00001/E/ASP/429/OD1 to model00001/A/ARG/424/NE: 2.589 model00001/E/ASP/429/OD2 to model00001/A/ARG/424/NE: 3.604 model00001/E/GLU/437/OE2 to model00001/E/ARG/310/NE: 3.311 model00001/E/ASP/446/OD2 to model00001/E/LYS/239/NZ: 3.776 Message: 3 Date: Fri, 12 Feb 2021 14:07:13 +0000 From: "Chen, Qiang" <q...@pitt.edu> To: "pymol-users@lists.sourceforge.net" <pymol-users@lists.sourceforge.net> Subject: Re: [PyMOL] PyMOL-users Digest, Vol 177, Issue 5 Message-ID: <bl0pr04mb466024e36e795ae57cbbc6a7fa...@bl0pr04mb4660.namprd04.prod.outlook.com> Content-Type: text/plain; charset="us-ascii" In case anyone has a similar problem, here is my solution from pymol import cmd from glob import glob cmd.do("run list_hb.py") for file in glob("model*.pdb"): print(file) cmd.load(file) obj=cmd.get_object_list('all') print(obj) #hydrogen bonds cmd.do("list_hb all, all, write_distances_file=hb_"+str(obj[0])+".txt") #salt bridge cmd.do("select negative, (resn ASP+GLU and name OD*+OE*)") cmd.do("select positive, (resn Lys and name NZ) or (resn ARG and name NE + NH*)") cmd.do("pairwise_dist negative, positive, 4, show=N, output="+str(obj[0]))
_______________________________________________ PyMOL-users mailing list Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net Unsubscribe: https://sourceforge.net/projects/pymol/lists/pymol-users/unsubscribe