PYMOL: Output Polar Contacts Between Chains to Text File
I am new to PyMOL but have a very specific task that I need to do.
I have a PDB structure file of protein homo-oligomer, and I want to use PyMol to determine polar contacts between chain A and all other chains.
However, I not only want to visualize these polar contacts (which is easy to do via PyMOL GUI), but I also want to be able to store/save pairs of interacting residues in a dictionary or whatever may be the most convenient data structure to do this. The end goal to create a sort of “interaction map” in a csv or text file where each interacting pair of residues has its own line.
I would ideally like to do this via a Python script where I call the PyMOL API or through a PYMOL script – whichever is easier.
I have reviewed a similar question on Stack Overflow (stackoverflow.com/questions/49633809/pymol-select-residues-for-which-i-have-found-hydrogen-bonds) but found it unhelpful because I not only want to select and visualize these polar contacts but also want to save/record them in a list or dictionary (and then to a text or csv file), which the above question does not address.
• 17 views