You are here

Segfault when loading scorefxn() and load_from_pdb()

1 post / 0 new
Segfault when loading scorefxn() and load_from_pdb()
#1

Hi,

I am trying to load a structure from PDB and then calculate pairwise residue energies between residues in the structure. For this, I want to first load the scorefxn and then the structure (or viceversa), but whenever I try, either of them gives me a segfault, i.e. if I do scorefxn = get_fa_scorefxn(); pose_from_pdb() the latter crashes and if I do the opposite, get_fa_scorefxn() crashes. 

 

This is a code that used to work fine, so I am a bit confused as to what might be going wrong. I know the segfault could be caused by trying to load wrong files or pointing at objects that don't exist, but I can't see anything wrong. I am pasting my code below, I was wondering if there is anything obviously wrong?

Thank you!

 

from rosetta import *
rosetta.init()
from toolbox import cleanATOM
from rosetta import get_fa_scorefxn
from rosetta.core.scoring import *
from Bio.PDB import *

class TCRpep_select(Select): #function to create TCR + peptide PDB
    def accept_chain(self, chain, alpha_chain = "A", beta_chain = "B", pep_chain = "P"):
        if chain.full_id[2] == str(alpha_chain) or chain.full_id[2] == str(beta_chain) or chain.full_id[2] == str(pep_chain):
            return 1
        else:
            return 0

Parser = PDBParser()
structure = Parser.get_structure(structure_id, folder + filename)

scorefxn = get_fa_scorefxn()

io = PDBIO()
io.set_structure(structure)
io.save(folder + "TCR_peptide_only/" + structure_id +"_TCR_peptide.pdb", TCRpep_select())
cleanATOM(folder + "TCR_peptide_only/"structure_id +"_TCR_peptide.pdb")
pdb_TCR_pep = pose_from_pdb(str(folder + "TCR_peptide_only/" + structure_id + "_TCR_peptide.clean.pdb"))



 

 

Category: 
Post Situation: 
Wed, 2020-07-22 03:00
mm