I relaxed the homodimeric biological assembly of PDB 1ISA, and am trying to use Interface Analyzer to calculate the interfacial energy between the subunits. I have cleaned the strucutre to removed all the water, metals, et.c and leave only the protein residues. Below are the command line flags I am using:
/media/Rosetta/rosetta_src_2016.13.58602_bundle/main/source/bin/InterfaceAnalyzer.linuxgccrelease -s /home/hannah/ISA_BIO_ignorechain_0001_0001.pdb -interface A_B -compute_packstat true -pack_input false -pack_separated true -no_optH false -optH_MCA true -flip_HNQ -use_input_sc -ex1 -ex2 -ex1aro_exposed -ex2aro_exposed -extrachi_cutoff 1 -atomic_burial_cutoff 0.01 -pose_metrics::interface_cutoff 8.0 -tracer_data_print false -out:file:score_only Interface_Analyzer_1ISA.sc -add_regular_scores_to_scorefile true
When I run Interface Analyzer on the file, it produces a segmentation fault. It fills the pose, calculates dSASA, calculates per-res data, builds rotamers at the identified positions, computes packstats, computes delta unsat polar residues, and calculates shape complementarity score. The error occurs in the next step when the function designates the upstream chain and the downstream chain.
I increased the interface cutoff from 8 Angstroms to as high as 20 Angstroms to test if the problem was interface detection. However, it produced the same segmentation fault each time. This protocol has worked on ~14 other protein structures. This is the only one to cause an error. Any insights would be helpful.