Author: Kevin Drew, (kdrew@nyu.edu)
Last updated December 7, 2012;
rosetta/main/source/src/apps/public/noncanonical_backbones/hbs_design.cc
rosetta/main/tests/integration/tests/hbs_design/ ) and demo folder ( rosetta/demos/public/hbs_design ).The main references for the hbs_design app are to be published in the PlosOne 2013 Rosetta Special Collection
Hydrogen bond surrogates (HBS) are stabilized alpha helices. The first hydrogen bond is replaced by a covalent linker. Many protein interaction interfaces involve helices and there is interest in using hbs for protein interaction inhibitors. The rational design of hbs using the hbs_design app attempts to find high affinity binders to target proteins.
hbs_design requires the following inputs:
| Flag | Description | Type | Default |
|---|---|---|---|
| -hbs_design_positions | positions on hbs chain to be designed | list of numbers [ex. 1 2 3 6] | None, only repack, no design |
| -pert_num | number of perturbations made during perturbation phase | integer | 10 |
| -design_loop_num | number of perturbation + design cycles | integer | 10 |
More information on common Rosetta flags can be found in the relevant rosetta manual pages. In particular, flags related to the job-distributor (jd2), scoring function, constraint files and packing resfiles are identical to those in any other Rosetta protocol).
| Flag | Description |
|---|---|
| -in::file::s Or -in:file:silent | Specify starting structure (in::file::s for PDB format, in:file:silent for silent file format). |
| -in::file::silent_struct_type AND -out::file::silent_struct_type | Format of silent file to be read in/out. For silent output, use the binary file type since other types may not support ideal form |
| -nstruct | Number of models to create in the simulation |
| -use_input_sc | Include rotamer conformations from the input structure during side-chain repacking. Unlike the -unboundrot flag, not all rotamers from the input structure are added each time to the rotamer library, only those conformations accepted at the end of each round are kept and the remaining conformations are lost. |
| -ex1/-ex1aro -ex2/-ex2aro -ex3 -ex4 | Adding extra side-chain rotamers |
| -database | The Rosetta database |
| -include_patches | Turn on patch files which are off by default. Turning on patches/hbs_pre.txt and patches/hbs_post.txt is a requirement for hbs applications. |
generate 1000 (or more) models
hbs_design.{ext} -database ${rosetta_db} -include_patches patches/hbs_pre.txt patches/hbs_post.txt -s start.pdb -nstruct 1000 -hbs_design_positions 1 2 3 4 -pert_num 100 -design_loop_num 10Sort models by total score, take top 5%, sort by REPACK_ENERGY_DIFF. Inspect top models.
Start with good binding mode: This application does not do large perturbations and therefore only minimally explores the space around the input structure. Use experimental evidience, a full docking app or estimate the binding mode by hotspot residues of native complex.
Target model: This application assumes the backbone of the target structure is in a fixed conformation and does not do flexible backbone moves.
Typical running time: With default settings and a typical input structure (~ 100 residues), a single model is estimated around 15 seconds. For a recommended 1000x models, the runtime should be around < 25 mins.
Multichain receptors: This app currently supports only the use of single-chain targets. If you want to use a multichain target, you may be able to get away with hacking your input PDB to re-label all target chains as one big chain (removing termini atoms and TER cards). Renumbering is necessary only to resolve numbering conflicts.
The output of a hbs_design run is a score file (score.sc by default) and k model structures (as specified by the -nstruct flag and the other common Rosetta input and output flags). The score of each model is the second column of the score file.
Model selection should be made based on sorting by total score, take top 5% and sort by REPACK_ENERGY_DIFF in the score file.