|
Rosetta 3.1 Release Manual
|
- Author:
- Steven Lewis smlewi@gmail.com
Last edited 2/2/10. Code and documentation by Steven Lewis smlewi@gmail.com.
The code is at rosetta_source/src/apps/public/scenarios/FloppyTail/; there's an integration test which serves as a demo at rosetta_source/test/integration/tests/FloppyTail/
- Kleiger G, Saha A, Lewis S, Kuhlman B, Deshaies RJ. Rapid E2-E3 assembly and disassembly enable processive ubiquitylation of cullin-RING ubiquitin ligase substrates. Cell. 2009 Nov 25;139(5):957-68. PubMed PMID: 19945379.
This code was written for a relatively singular application. The system in question was a protein with a long (dozens of residues) flexible tail, which was not seen in crystal structures. Biochemical evidence suggested a particular binding site for the tail on a known binding partner (the two binding partners also had a known binding interface separate from this tail). The code was intended to model the reach of the long flexible tail and determine whether the hypothesized binding site was plausible.
The algorithm is fairly simple: small/shear/fragment moves in centroid mode to collapse the tail into some sort of folded conformation from an initially straight-out-into-space extended conformation, and small/shear moves with repacking to refine its position.
This code can also model internally flexible regions between domains, although that is not its intended purpose. It may not do a good job of detecting appropriate interfaces for repacking in this latter case.
The code is compatible with constraints. Early modeling proceeded using constraints and some small hacks to help guide models to the hypothesized tail-binding site. Ultimately this was not necessary for the original system, but the code retains the ability to use constraints, etc. Your mileage may vary.
This code is NOT intended to do "half-abinitio" where you know half a structure and want to fold the other half. Although it is modeled on abinitio, it is only tested on a truly floppy tail, and I have no idea if it is able to fold compact structures.
See test/integration/tests/FloppyTail/ for example usage.
This code was intended for a single purpose, but it may work if you have a similarly flexible tail. It can also model internal flexible regions between domains.
- You are likely to need a large number of trajectories if you have a long tail. The protocol gets trapped in bad conformations fairly easily even with the Metropolis criterion. Production runs for publication completed slightly less than 30000 trajectories in 1000 processor-days on 2.3 gigahertz processors.
- Constraints are only used in centroid mode in this protocol. Change that if you want it changed. Some constraint reporting is deactivated by comments, you should comment it back in if you want constraints.
- FloppyTail.cc includes extra_analysis.hh in the same folder. This code is riddled with system-based assumptions and magic numbers (in other words it directly prints data specific to the original system, like the distance between two residues of interest). You'll want to turn off the function create_extra_output at line 392 in FloppyTail, or else replace extra_analysis with your own code.
- This code is not compatible with silent-file input, but you could rewire the PDB-numbering based inputs into resid based inputs without much difficulty and it would work.
- When modeling a terminal flexible region (a tail), in the refinement phase, the protocol can be directed to model a shorter portion of the tail for part of refinement mode with the short_tail_xxx options. The reasoning is that the tail may be too close to a binding partner in the structure fresh from centroid mode, even with repacking (centroids tend to be a bit too small for this sort of docking). By only remodeling the tip of the tail in the first part of refinement, you can relax clashes without swinging the tail back out into space.
- The short_tail_xxx, constraint, and pair_off options were not used for production runs with this code. They were used for early experiments, controls, debugging, etc. They ought to still work.
FloppyTail supports three types of options: general Rosetta 3.x options (packing, etc.), generic protocol options like "how many cycles" borrowed from the AnchoredDesign application (not included in the release), and FloppyTail specific options.
FloppyTail options
- flexible_start_resnum - integer - this is the start of the flexible tail in PDB numbering.
- flexible_stop_resnum - integer - this is the end of the flexible region, in PDB numbering. Passing 0 or not using this option means the entire chain after flexible_start_resnum.
- flexible_chain - string - the first character of this string is interpreted as the PDB chain for the flexible region; any other characters are ignored.
- shear_on - real - In centroid mode, shear moves are completely nonproductive early on when the tail is still largely extended. This value gives the fraction of centroid cycles when shear moves will be allowed (introduced into the moveset of the RandomMover choosing perturbation moves). For example, passing 0.333 means that for the first third of centroid mode, shear moves will be disallowed.
- short_tail_fraction - real - Fraction of the tail used in the short tail fraction of refinement mode. 0.1 would mean the last tenth of the tail is flexible. Not compatible with non-terminal flexible regions.
- short_tail_off - real - Fraction of refinement cycles dedicated to refining only the short part of the tail. 0.33 means the first third of refinement cycles will be with the shorter flexible region.
- pair_off - boolean - If true, disable the electrostatic Epair (pair and fa_pair) terms. Used for a control experiment, not for general use.
AnchoredDesign options (borrowed for simplicity, not tied to AnchoredDesign in any other way); all are in the AnchoredDesign namespace
- AnchoredDesign::perturb_temp - real - Monte Carlo temperature for perturb phase (0.8 used for production)
- AnchoredDesign::perturb_cycles - unsigned integer - number of perturb phase cycles (5000 used for production)
- AnchoredDesign::perturb_show - boolean - if true, outputs centroid poses after perturbation
- AnchoredDesign::debug - debug - if true, outputs poses for each monte carlo cycle
- AnchoredDesign::refine_temp - real - Monte Carlo temperature for refine phase (0.8 used for production)
- AnchoredDesign::refine_cycles - unsigned integer - number of refine phase cycles (3000 used for production)
- AnchoredDesign::refine_repack_cycles - unsigned integer - Perform a repack/minimize every N cycles of refine mode (30 used for production)
General options: All packing namespace options loaded by the PackerTask are respected. jd2 namespace options are respected. Anything very low-level, like the database paths, is respected.
- packing::resfile - string - a resfile if you want one
- in::file::frag3 - string - fragments if you've got them
- run::min_type - string - minimizer type. dfpmin_armijo_nonmonotone used for production.
Generated on Tue Apr 20 07:50:06 2010 for Rosetta Projects by
1.5.2
© Copyright Rosetta Commons Member Institutions. For more information, see http://www.rosettacommons.org.