The scripts and input files that accompany this demo can be found in the
demos/public directory of the Rosetta weekly releases.
KEYWORDS: MEMBRANES DOCKING
This application assembles and docks symmetric protein complexes in the membrane bilayer. The symmetric native complex is first refined using the MP_Relax application. The lowest scoring native conformation (by total Rosetta score) is then used as input to the membrane symmetric docking application, which searches for possible conformations by reassembling and docking subunits together.
This application combines the membrane framework, symmetry machinery, and standard symmetric docking algorithm in Rosetta. Currently, docking of Cyclic (C) symmetries is supported.
Reference for this protocol capture:
Three initial input files are required for this protocol:
Steps for generating these inputs are provided below. These inputs can also be found in the example_inputs/ directory
PDB File: Generate a PDB file where the membrane protein structure is transformed into PDB coordinates (z-axis is membrane normal). This can be done either by downloading the transformed PDB directly from the PDBTM website (http://pdbtm.enzim.hu/) or by downloading a PDB file from the PDB and running it through the PPM server (http://opm.phar.umich.edu/server.php).
We use the potassium channel KcsA as an example here:
Clean the PDB using
Rosetta/tools/protein_tools/scripts/clean_pdb.py 1bl8_tr.pdb ignorechain
(the resulting PDB is referred as 1bl8_tr.pdb from here)
Full (symmetric) and asymmetric Span File: Generate a spanfile from the PDB structure using the spanfile_from_pdb application described in the MP_spanfile-from-pdb protocol capture in Rosetta/demos/protocol_captures/2015.
Rosetta/main/source/bin/spanfile_from_pdb.linuxgccrelease -database /path/to/Rosetta/main/database -in:file:s 1bl8_tr.pdb
For this example, this command will produce 5 output files:
1bl8_tr.span: Predicted trans-membrane spans for the full symmetric complex
1bl8_tr<A-D>.span: Predicted trans-membrane spans for each chain in the complex
Here, we describe the steps required to run the MP_SymDock protocol. As an example, all steps use a C4 Symmetric Potassium Channel (PDB ID: 1bl8)
$> $ROSETTA3/bin/rosetta_scripts.default.linuxgccrelease -parser:protocol membrane_relax.xml @relax_flags
The following output files will be generated:
1bl8_tr_<0001-0010>.pdb: 10 refined models of 1bl8
relax_scores_1bl8.sc : Rosetta scores for each resulting models
Examples of these outputs can be found in example_refined_models/
Note on timing: Depending on protein size, refinement is a time consuming step. Each decoy will take 0.5-1.0hrs depending on avaialble processing power.
Input model selection: Use the score file from the refinement step to select the lowest scoring refined model by total Rosetta score. In this example, the model 1bl8_tr_0009.pdb has the lowest total Rosetta score and will therefore be used as input to the next step. From this point forward, this model will be referred to as 1bl8_refined.pdb and is located in example_symmetry_files/.
Generate inputs for symmetry: To prepare the structure for assembly and docking in the protocol, a set of asymmetric inputs must be generated. These inputs describe the asymmetric unit, which will later be used to re-assemble the complex based on a generated symmetry definition. A version of these generated input files are provided in the example_symmetry_files/ directory
First, create the asymmetric input structure and symmetry definition file from the refined symmetric complex using the make_symmdef_file.pl script. An example commandline is provided below:
$> $ROSETTA3/src/apps/public/symmetry/make_symmdef_file.pl -p example_symmetry_inputs/1bl8_refined.pdb -a A -i B:4 > 1bl8.c4.symm
In this command-line:
-pspecifies the input PDB file
-aspecifies the chain or chains to use as the asymmetric unit,
-ispecifies how to organize the remaining chains. In this example, "B:4" means use chain B as the next subunit and arrange subunits as a C4 tetramer (4 subunits around the Z axis))
This command will generate various files, only two of which are key here:
1bl8_refined_INPUT.pdb: PDB file containing the asymmetric subunit
1bl8.c4.symm: Symmetry definition file describing the arrangement of subunits in the complex. This file specifically describes needed translations and rotations to regenerate and assemble this complex from the input file.
Note: To generate a correct symmetry, make_symmdef_file.pl requires all chains be of equal length. Subunits should also be close to <0.5Å rmsd to one another. Any asymmetry may result in an incorrect symmetry definiton. To check, you can visualize the example_symmetry_inputs/1bl8_refined_symm.pdb to ensure this initial setup is correct.
Next, you will need the 1bl8_trA.span file containing trans-membrane spans for only the asymmetric unit. If your asymmetric unit contains multiple chains, you may need to assemble this file yourself from the full set of spans.
Running the symmetric docking application: Using the asymmetric unit PDB, symmetry definition file, and asymmetric unit span file as inputs, you are ready to run the membrane symmetric docking application. Flags, recommended settings, and commandlines are described below:
Required Options: Options (flags) needed to run this application are described below. A file with these flags,
symdock_flags, is also provided for the 1bl8 example.
flags descriptions -------------------------------------------------------------------------------------------------- -in:file:s <pdbfile> Input PDB Structure: Asymmetric input structure -in:file:native <pdbfile> Structure of native symmetric complex for RMSD calculations -membrane_new:setup:spanfiles <spanfile> Spanfile describing spanning topology of asymmetric unit -membrane_new:scoring:hbond Turn on membrane depth-dependent hydrogen bonding weight -symmetry:symmetry_definition Symmetry definition file -symmetry:initialize_rigid_body_dofs Locally sample rigid body conformations during intial complex assembly (before docking algorithm) -nstruct Number of structures to generate -packing:pack_missing_sidechains 0 Wait to pack until the membrane mode is turned on -docking:dock_lowres_filter 5.0 10.0 Lower van der Waals scoring criteria during centroid stage to allow wider range of rigid body sampling; Both numbers are highest allowed scores for a low-resolution model to proceed to the high-resolution stage: 5.0 = highest allowed van der Waals score 10.0 = highest allowed interchain contact score
Recommended # of Decoys
Rosetta/main/source/bin/membrane_symdocking.linuxgccrelease -database /path/to/Rosetta/main/database @symdock_flags
The following outputs will be generated from the symmetric docking protocol. A version of these outputs are also provided in the example_outputs/ directory:
1bl8_refined_INPUT_0001.pdb: Symmetrically docked output model from the
score.sc: Scorefile output by Rosetta containing memrbane and symmetry
scores for this model
DiMaio F, Leaver-Fay A, Bradley P, Baker D, André I (2011) Modeling Symmetric Macromolecular Structures in Rosetta3. PLoS ONE 6: e20450.
Barth P, Schonbrun J, Baker D (2007) Toward high-resolution prediction and design of transmembrane helical protein structures. Proc Natl Acad Sci 104: 15682–15687.