You are here

Sequence design with mixed fixbb-relax protocol

5 posts / 0 new
Last post
Sequence design with mixed fixbb-relax protocol
#1

Hello,

I would like to design a protein sequence using a mixed fixbb-relax protocol.

Basically, I want to do something like in this publication: http://www.sciencedirect.com/science/article/pii/S0969212612001517

(The idea is to iterate between fixbb and relax protocols as long as the Rosetta energy difference between models is bigger than 1.0 REU)

Does anybody have an idea how to use RosettaScripts to implement this?
I understand that my script(below) runs fixbb, takes the best trajectory and inputs it into relax and I get output model that has gone through fixbb and relax only once. But I would like more cycles to occur and each time energy difference between cycle i and i-1 should be calculated to determine if next cycle should be made.

I attach what I have done so far, XML file (as txt attachment, no other way), command and flags:

~/rosetta3.5/rosetta_scripts.
default.linuxgccrelease @flags_designrelax -out:file:silent results.out -parser:protocol design_relax.xml

-s INPUT.pdb
-nstruct 2
-database /mnt/db/rosetta_database3.5/
-ex1
-ex2
-ndruns 3

best wishes,
Krzysiek

AttachmentSize
design_relax.txt741 bytes
Post Situation: 
Thu, 2013-10-03 08:11
kszczepaniak

So, I checked out the supplemental info, expecting a RosettaScript, and I saw a command line that involves an unreleased PilotApp.

However, if you take a look at Rosetta/Demos/public/ideal_proteins, there is a mover called in the example XML flxbbdesign that iterates between relax and design. Use the example from the ideal proteins script in your XML. The reference is here:

### Principles for designing ideal protein structures ###
### Nobuyasu Koga, Rie Tatsumi-Koga, Gaohua Liu, Rong Xiao, Thomas B. Acton, ###
### Gaetano T. Montelione, and David Baker ###
### Nature 2012

You may want to look at the packstat filter as well...
You could use the LoopOver mover with a score filter, but I'm not exactly sure how to get exactly what you want. Maybe someone else can chime in who has more experience with RS than I do..

Hope this helps at least a bit...

-J

Thu, 2013-10-03 09:20
jadolfbr

If you do want to automate this protocol, you may want to use the IteratedConvergence mover. ( https://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/Mov...(RosettaScripts)#IteratedConvergence ) This will repeat a sub-mover (which can be a multi-step ParsedProtocol https://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/Mov...(RosettaScripts)#ParsedProtocol_.28formerly_DockDesign.29 mover) until a given filter (e.g a ScoreType filter for total_score https://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/Fil...(RosettaScripts)#ScoreType ) doesn't change by more than a given amount. Your ParsedProtocol mover can contain both the PackRotamersMover ( Which does basically the same thing as the fixbb design https://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/Mov...(RosettaScripts)#PackRotamersMover ) and the FastRelax mover (which does basically the same thing as the relax application https://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/Mov...(RosettaScripts)#FastRelax )

Tue, 2013-10-08 08:28
rmoretti

The flxbbdesign application is indeed handy, but has no documentation yet therefore it is hard to use it, but thanks for the idea.
In the end I applied procedure proposed by rmoretti and it works great!

Thank you both for help :)

Krzysiek

Fri, 2013-10-11 05:23
kszczepaniak

Hi,

could you share the .xml file you used to apply this procedure?

 

 
 
 
Wed, 2020-09-23 11:24
fede