From 0a8b0bd1378fb78aec377991ef73338d4156732e Mon Sep 17 00:00:00 2001
From: Gabriel Studer <gabriel.studer@unibas.ch>
Date: Thu, 22 Sep 2022 09:55:46 +0200
Subject: [PATCH] chain mapping: handle case of no surviving residues in
ProcessStructure
---
modules/mol/alg/pymod/chain_mapping.py | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/modules/mol/alg/pymod/chain_mapping.py b/modules/mol/alg/pymod/chain_mapping.py
index 1d3661340..2307ee2b1 100644
--- a/modules/mol/alg/pymod/chain_mapping.py
+++ b/modules/mol/alg/pymod/chain_mapping.py
@@ -1339,6 +1339,10 @@ class ChainMapper:
polypep_seqs = seq.CreateSequenceList()
polynuc_seqs = seq.CreateSequenceList()
+ if len(view.residues) == 0:
+ # no residues survived => return
+ return (view, polypep_seqs, polynuc_seqs)
+
for ch in view.chains:
n_res = len(ch.residues)
n_pep = sum([r.IsPeptideLinking() for r in ch.residues])
--
GitLab