Member Site › Forums › PyRosetta › PyRosetta – General › DockingHighResLegacy repacks all sidechains on first apply(), interface sidechains on subsequent apply()
- This topic has 4 replies, 4 voices, and was last updated 5 years, 4 months ago by
Anonymous.
-
AuthorPosts
-
-
February 23, 2012 at 6:46 am #1171
Anonymous
[Using PyRosetta.ScientificLinux-r47452.64Bit]
docker = DockingHighResLegacy(jump, scoreFxn)
docker.apply(pose)
When I use the above in a python script, I get:
…
protocols.docking.DockingHighRes: in DockingHighRes.apply
protocols.docking.DockingHighRes: ::::::::::::::::::DOCK_MCM:::::::::::::::::::
core.pack.task: Packer task: initialize from command line()
core.pack.interaction_graph.interaction_graph_factory: Instantiating DensePDInteractionGraph
core.pack.pack_rotamers: built 4141 rotamers at 256 positions.
…
which repacks sidechains of my entire structure. However, the next time I use .apply(pose), I get:
…
protocols.docking.DockingHighRes: in DockingHighRes.apply
protocols.docking.DockingHighRes: Using user-defined TaskFactory.
protocols.docking.DockingHighRes: ::::::::::::::::::DOCK_MCM:::::::::::::::::::
core.pack.task: Packer task: initialize from command line()
core.pack.task: Packer task: initialize from command line()
core.pack.task: Packer task: initialize from command line()
core.pack.interaction_graph.interaction_graph_factory: Instantiating DensePDInteractionGraph
core.pack.pack_rotamers: built 30 rotamers at 2 positions.
…
which just repacks the interface.I see in Rosetta 3.3’s DockingHighResLegacy constructor that “init_task_factory_=NULL” is active, while in ver 3.2.1 “//init_task_factory_=NULL”.
Is there a way to turn off pre-packing in the DockingHighResLegacy class? Is there now a better class to dock with?
-
March 2, 2012 at 4:58 am #6734
Anonymous
Hello,
I would definitely recommend using the standard DockingHighRes mover for docking rather than the legacy version. I’m not quite sure how to prevent the initial full-repack of DockingHighResLegacy… You could create a packer task and set “or_include_current” to True (task=TaskFactor.create_packer_task(pose); task.or_include_current(True)). This will ensure that the original rotamers in your pose are considered during packing; so if you are starting from an xtal structure or something stable there is a good chance Rosetta won’t touch many of the side chains. The reason the full-prepack is necessary is to remove any high-energy clashes from the initial structure, if you skip this step you may get poor results as a result of rounding errors (during docking the total energy drop might be minuscule compared to the high energy clash). Hope this helps!
Mike
-
March 10, 2012 at 7:44 pm #6799
Anonymous
Appreciate the help.
I was using
DockingHighResLegacy(...)
because I am unable to instantiateDockingHighRes(...)
with PyRosetta. Have you been able to get this working?As for repacking sidechains, I’m doing this before docking. Since it is stochastic, I would rather not have the Docker repack all of them on its own before each run. This leads to docked structures that are not consistent, though only slightly.
I’m currently using
DockingHighResMCM(..)
and manually minimizing after. -
March 17, 2019 at 1:18 pm #14606
Anonymous
7 years later,i found the same problem,do u know how to instantiate
DockingHighRes with PyRosetta?
-
May 9, 2019 at 10:03 pm #14714
Anonymous
One option you might try is to initialize the mover through the XML interface, as that is often better tested/supported than the direct Python initialization.
An example of how to use the XML interface to initialize a mover can be found at https://github.com/RosettaCommons/pyrosettascripts_demo/blob/master/data_generation/ddG_pssm.ipynb (You would just need to adjust the tags for the mover you’re interested in initializing.)
-
-
AuthorPosts
- You must be logged in to reply to this topic.