pdb2sql icon indicating copy to clipboard operation
pdb2sql copied to clipboard

compute_irmsd_fast() only works on two chains

Open LilySnow opened this issue 3 years ago • 6 comments

>>> ref = "ref.pdb"
>>> model = "model.pdb"
>>> sim = StructureSimilarity(model,ref)
>>> irmsd_fast = sim.compute_irmsd_fast()
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/home/lixue/tools/anaconda3/lib/python3.7/site-packages/pdb2sql/StructureSimilarity.py", line 277, in compute_irmsd_fast
    resData = self.compute_izone(cutoff, save_file=False)
  File "/home/lixue/tools/anaconda3/lib/python3.7/site-packages/pdb2sql/StructureSimilarity.py", line 328, in compute_izone
    'exactly two chains are needed for irmsd calculation but we found %d' % len(chains), chains)
ValueError: ('exactly two chains are needed for irmsd calculation but we found 5', ['A', 'B', 'C', 'D', 'E'])

Archive.zip

LilySnow avatar Jul 11 '22 17:07 LilySnow

This issue is stale because it has been open for 30 days with no activity. Remove stale label or comment or this will be closed in 7 days.

github-actions[bot] avatar Aug 12 '22 03:08 github-actions[bot]

This issue was closed because it has been inactive for 7 days since being marked as stale.

github-actions[bot] avatar Aug 19 '22 04:08 github-actions[bot]

To solve this issue, I am planning to add an option chain_pairs = ['D:E', 'AB:C'] to compute_izone(). This way, pdb2sql knows which interfaces are of interest when calculating i-RMSD for a multiple-chain complex.

LilySnow avatar Aug 23 '22 12:08 LilySnow

i've suggested to make a few change to Li's work see #80. This basically allows to use get_contact_atoms using multiple chains as :

db.get_contact_atoms(chain1='AB', chain2='C')

NicoRenaud avatar Aug 24 '22 14:08 NicoRenaud

This issue is stale because it has been open for 30 days with no activity. Remove stale label or comment or this will be closed in 7 days.

github-actions[bot] avatar Sep 26 '22 05:09 github-actions[bot]

Related to PR #80 and #79

NicoRenaud avatar Oct 03 '22 12:10 NicoRenaud