Member Site › Forums › Rosetta 3 › Rosetta 3 – Applications › mover “RepackMinimize” does not support chains > 2?
- This topic has 6 replies, 3 voices, and was last updated 11 years, 3 months ago by
Anonymous.
-
AuthorPosts
-
-
May 20, 2013 at 9:17 am #1587
Anonymous
In RosettaScript, the mover “RepackMinimize” has attributes “design_partner1” and “design_partner2” to define which partner will be designed, but it does not have attribute like “jumps” or “partners” (just like “DockingProtocol”) to let you set the jump in fold tree. It assumes to use jump=1, this is OK for two-chain complex. But it is not OK for more than two chains.
I applied this mover on a 4-chain complex with “minimize_rb=1”:
<RepackMinimize name=des1 scorefxn_repack=soft_rep scorefxn_minimize=soft_rep minimize_bb=0 minimize_rb=1/>But the program halted with errors:
protocols.rosetta_scripts.ParsedProtocol: =======================BEGIN MOVER RepackMinimize=======================
{
core.pack.task: Packer task: initialize from command line()
core.pack.interaction_graph.interaction_graph_factory: Instantiating PDInteractionGraph
core.pack.pack_rotamers: built 6972 rotamers at 79 positions.
core.pack.pack_rotamers: IG: 17487320 bytesERROR: min_rb.size() == pose.num_jump()
ERROR:: Exit from: src/protocols/protein_interface_design/design_utils.cc line: 365 -
May 20, 2013 at 6:38 pm #8790
Anonymous
First off, even the documentation (https://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/Movers_%28RosettaScripts%29#RepackMinimize) doesn’t really recommend using RepackMinimize – using a separate PackRotamersMover and MinMover (or TaskAwareMinMover) will allow for better control.
As for the error – it’s a horrible bug. It’s expecting min_rb to be vector of bools, one for each jump, but it’s assigning a single bool to the entry. And no, there’s no way to get around it. Either don’t set minimize_rb (it should default to all true anyway), or use use the separate PackRotamersMover/MinMover.
-
May 20, 2013 at 8:28 pm #8792
Anonymous
Rocco filed the bug here: https://carbon.structbio.vanderbilt.edu/mantisbt/view.php?id=274
-
May 21, 2013 at 1:19 am #8795
Anonymous
Thank you for pointing out the bug.
Actually, the XML script I used here was adopted from protocol_capture/RosettaScripts/pp-interfaceredesign/flexbb-interfacedesign.xml, which was originally published in the paper ‘RosettaScripts: A Scripting Language Interface to the Rosetta Macromolecular Modeling Suite’. This script has four ‘RepackMinimize’ mover with one ‘Backrub’ mover to do the design. -
May 21, 2013 at 1:56 am #8796
Anonymous
I also found mover “BackrubDD” has attribute “partner1” and “partner2” but no “jumps”. I have checked the source codes of “BackrubDD” in protocols/protein_interface_design/movers/BackrubDDMover.cc, it has conditional block:
if( pose.conformation().num_chains() == 2 ){
}
else if ( !residues_.size() ) { // pose does not have 2 chains, backrub all (protein)
}
which means for structures with more than 2 chains, the “partner1” and “partner2” are meaningless, all residues will be moved. -
May 21, 2013 at 6:43 pm #8804
Anonymous
Submitted to the tracker: https://carbon.structbio.vanderbilt.edu/mantisbt/view.php?id=275
By the way, I believe the Rosetta bug tracker should be open from submissions from non-developers. If you run across other bugs, feel free to make a bug tracker account and submit them yourself.
-
May 22, 2013 at 1:58 am #8806
Anonymous
many thanks!
Hope to get version 3.5 released soon.
-
-
AuthorPosts
- You must be logged in to reply to this topic.