The scripts and input files that accompany this demo can be found in the
demos/public directory of the Rosetta weekly releases.
KEYWORDS: LIGANDS DOCKING INTERFACES
We are predicting the conformation of the complex of FKBP12, FRAP, and rapamycin. Rapamycin is a dimerizer that allows FK506-binding-protein (FKBP12) to form an interface with FKBP-rapamycin-associated protein (FRAP). The first section of this tutorial demonstrates how to prepare input files for the proteins and small molecule. The second section describes docking rapamycin to FKBP12. In the third section we will dock the result from section 2 with FRAP.
To complete this quest you must have at least a 3 member party, including a cleric (level 43, must have blessing spell), a warrior proficient with hammer, and a thief (unlock skill level 10).
This demo was written by Gordon Lemmon, Sergey Lyskov, and Loren Looger.
Unzip the PDB file 1FAP.pdb.gz. In order for the ligand docking to work correctly the 1-letter chain identifier for the ligand must be different from the protein chain ids. Look inside the file 1FAP.pdb. On line 307 and 308 we find that chain A is FKBP12 and chain B is FRAP. Toward the bottom of the file Rapamycin is specified by the residue id RAP (2375-2442). Make a new file with just the RAP lines.
grep HETATM 1FAP.pdb | grep RAP > rap.pdb
Using your favorite text editor change the chain id found in rap.pdb from A to X. Use clean_pdb.py (where to find?) to prepare 1FAP.pdb for Rosetta.
clean_pdb.py 1FAP.pdb A # this should output a file with only atom records from chain A, 1FAP_A.pdb clean_pdb.py 1FAP.pdb B # this should output a file with only atom records from chain B, 1FAP_B.pdb
We now have a separate PDB file for both proteins and an additional PDB file
for the ligand. We now must add hydrogens to our rapamycin molecule and save it
in MOL format. Simply open the file rap.pdb in Pymol and add hydrogens using
the action menu
all->A->hydrogens->add (or type
h_add on the command line).
Then and save it as a MOL file by using the file menu (
molecule->ok, select type as MOL file, and change the extension to
.mol). You should now find the file rap.mol in your directory.
Rosetta requires a PARAMS file for each ligand. These files describe the
atoms, bonds and bond angles within the ligand. To make a params file for
rapamycin use the script
molfile_to_params.py found here:
/rosetta_source/src/python/apps/public/molfile_to_params.py -c -n RAP rap.mol
We use the -c option to produce centroid mode params used in Part 3 of this demo.
Notice the warnings that are produced by the script. These are informing us that the ligand we are using is large and flexible, which means we will struggle to sample all of its flexibility during docking. Since we are starting with the correct conformation of Rapamycin we can ignore these warnings.
mol_to_params.py should have created a file called RAP_0001.pdb which has the
same coordinates as rap.pdb but has been prepared for use with Rosetta.
Combine 1FAP_A.pdb and RAP_0001.pdb into a new file:
cat 1FAP_A.pdb RAP_0001.pdb > FKBP+RAP.pdb
ligand_dock.xml files from
rosetta_source/src/test/integration/tests/ligand_dock_scripts to your
directory. We will use these as a starting point for our docking script.
We have modified the flags file to be specific for our study.
We now modify the options in this file to be specific for our study. First we
comment out the start_from lines, since our ligand is already in the correct
starting position. Other important options to consider optimizing include the
angstroms option of
Translate should represent
the size of your binding pocket (your ligand will move within a sphere with a
radius of this size).
Now we are ready to run our ligand docking protocol:
This should produce a file with a model of rapamycin docked to FKBP:
FKBP+RAP_0001.pdb. This file serves as an input to protein docking.
Combine 1FAP_B.pdb with FKBP+RAP_0001.pdb. Put ATOM lines from 1FAP_B first, followed by ATOM lines from FKBP+RAP.pdb, and then HETATM lines from FKBP+RAP.pdb.
egrep 'ATOM|HETATM' 1FAP_B.pdb FKBP+RAP_0001.pdb > combined.pdb
Prepare a flag file that specifies the centroid and full-atom PARAMS files for rapamycin. Also specify combined.pdb as the input file. Run the docking protocol:
This should produce an output file,
combined_0001.pdb. Using pymol you can
see that the FKBP/RAP complex has moved relative to FRAP.
For a production run you will want to run this protocol 10,000 or more times. Then find your best scoring models. An alternative strategy would be to produce thousands of models with Part 1 of this tutorial, then filter for the top few models of FKBP with RAP. Use each of the top models as inputs for part 2, producing several thousand models for each of these inputs.