You are here

Issue with Rosetta LinkResidues Script

2 posts / 0 new
Last post
Issue with Rosetta LinkResidues Script
#1

Hi All,

I am trying to design a 2 fold-symmetrical protein where I want to design a selected number of residues in each half and not mutate the remaining residues in both the halves. For e.g., what I want Rosetta to do is in "Half-A, the residues at position 32, 49, 56, 65, 66 should be mutated to certain residues, and the exact same mutated residues should appear after designing position 32, 49, 56, 65, 66 in Half-B.

I found that the Rosetta LinkResidues can do it. So, I tried to write a Rosetta Script with all the necessary information from the LinkResidues page of Rosetta. Note: I took into consideration the suggestions for possible bugs to avoid Segmentation fault.

But when I run my script with my resfile and flags files, I am getting "Segmentation fault (core dumped)" error without any further information.

I am attaching all the files and scripts that I have used. If any expert can help me in resolving this probelm, that would be really helpful!

Script I am using

<ROSETTASCRIPTS>
    <TASKOPERATIONS>

        <ReadResfile name="resfile" filename="my_resfile.res"/>
        <LinkResidues name="linkres">
            <LinkGroup group="32,49,56,65,66"/>
        </LinkResidues>

    </TASKOPERATIONS>

    <MOVERS>

        <PackRotamersMover name="packer" scorefxn="talaris2014" task_operations="resfile,linkres"/>

    </MOVERS>

    <PROTOCOLS>
        Runs movers and filters in this order
        <Add mover_name="packer"/>
    </PROTOCOLS>

</ROSETTASCRIPTS>

 

Sincerely

Aditya

AttachmentSize
flags file276 bytes
resfile100 bytes
input PDB file50.77 KB
Category: 
Post Situation: 
Wed, 2019-11-20 22:00
adityapadhi

Hi,

I can't answer your question, but I share your frustration with the link residue system. I created a task operation that replaces the residue linking for simple systems (see link below) but it looks like it won't help your case. Are you comfortable running a debugger like GDB to provide a backtrace when the program segfaults?

 

https://www.rosettacommons.org/docs/latest/scripting_documentation/RosettaScripts/TaskOperations/taskoperations_pages/KeepSequenceSymmetry

 

Best,

-Jack Maguire

Fri, 2019-11-22 13:45
JackMaguire