gf712 / abpytools Goto Github PK

File "/home/misha/anaconda3/lib/python3.6/site-packages/abpytools/core/chain.py", line 71, in load_from_string
    new_chain.load()
  File "/home/misha/anaconda3/lib/python3.6/site-packages/abpytools/core/chain.py", line 92, in load
    self.load()
  File "/home/misha/anaconda3/lib/python3.6/site-packages/abpytools/core/chain.py", line 101, in load
    self.hydrophobicity_matrix = self.ab_hydrophobicity_matrix()
  File "/home/misha/anaconda3/lib/python3.6/site-packages/abpytools/core/chain.py", line 217, in ab_hydrophobicity_matrix
    sequence=self._sequence)
  File "/home/misha/anaconda3/lib/python3.6/site-packages/abpytools/core/chain.py", line 434, in calculate_hydrophobicity_matrix
    else 0 for x in whole_sequence])
  File "/home/misha/anaconda3/lib/python3.6/site-packages/abpytools/core/chain.py", line 434, in <listcomp>
    else 0 for x in whole_sequence])
KeyError: 'X'
We've got an error while stopping in post-mortem: <class 'KeyboardInterrupt'>

def calculate_hydrophobicity_matrix(whole_sequence, numbering, aa_hydrophobicity_scores, sequence):
    # instantiate numpy array (whole sequence includes all the amino acid positions of the VH/VL, 
    # even the ones
    # that aren't occupied -> these will be filled with zeros
    # hydrophobicity_matrix = np.zeros(len(whole_sequence))
    # #
    # # # iterate through each position
    # for i, position in enumerate(whole_sequence):
    # 
    #     if position in numbering:
    #         position_in_data = numbering.index(position)
    #         hydrophobicity_matrix[i] = aa_hydrophobicity_scores[sequence[position_in_data]]

    #  return hydrophobicity_matrix
    # same thing as above but in a comprehension list
    return np.array([aa_hydrophobicity_scores[sequence[numbering.index(x)]] if x in numbering
                     else 0 for x in whole_sequence])

 for sequence_alignment in sequences:
    chain = Chain.load_from_string(sequence=sequence_alignment)
    cdrs = dict()
    try:
      for cdr_name, cdr_indexes in chain.cdr[0].items():
        cdrs[cdr_name] = sequence_alignment[cdr_indexes[0]:(cdr_indexes[-1] + 1)]
    except Exception as error:
         time.sleep(0.1)
         if chain.cdr == None:
            print("ERROR::", error)