MR using CCP4i:TOXD
Data for this tutorial are found here
Reflection data: toxd.mtz Structure files: 1BIK.pdb, 1D0D_B.pdb Sequence file: toxd.seq
This tutorial demonstrates the ensembling procedure in Phaser.
α-Dendrotoxin (TOXD, 7139Da) is a small neurotoxin from green mamba venom. You have two models for the structure. One is in the file 1BIK.pdb, which contains the protein chain from PDB entry 1BIK, and the other is in the file 1D0D_B.pdb, which contains chain B from PDB entry 1D0D. 1BIK is the structure of Bikunin, a serine protease inhibitor from the human inter-α-inhibitor complex, with sequence identity 37.7% to TOXD. 1DOD is the complex between tick anticoagulant protein (chain A) and bovine pancreatic trypsin inhibitor (BPTI, chain B). BPTI has a sequence identity of 36.4% to TOXD. Note that models making up an ensemble must be superimposed on each other, which has not yet been done with these two structures. At relatively low sequence identity, it is also a good idea to trim unconserved loops and the ends of non-identical side chains, so that the model comprises a conserved core.
- Start the ccp4 GUI by typing ccp4i at the command line.
- Make a new project called "phaser_tute" using the Directories&ProjectDir button on the RHS of the GUI. Set the "Project" to phaser_tute and "uses directory" to the directory where the files for this tutorial are located, and make this the "Project for this session of the CCP4Interface". You will then be able to go directly to this directory in the GUI using the pull-down menu that appears before every file selection.
- Go to the Molecular Replacement module, in the yellow pull-down on the LHS of the GUI
- Under Model Generation within the Molecular Replacement module, bring up the GUI for Sculptor
- Run sculptor on each of the MR models, using the alignment information in toxd_sequences.aln. Just use protocol 12 (the best overall in a series of trials) of the Predefined protocols to modify the starting models.
- It is a good idea to fill in the TITLE.
- Under Model Generation, now bring up the GUI for Ensembler
- Run ensembler to combine the two sculpted models.
- Check the resulting merged ensemble PDB file using a graphics program such as coot.
- Bring up the GUI for Phaser
- All the yellow boxes need to be filled in.
- It is a good idea to change the Ensemble id from the default.
- By default, Phaser will attempt to use intensities rather than amplitudes. Because there are only amplitudes available for this data set, uncheck the box labelled "Input data are merged intensities".
- The merged ensemble file from ensembler contains REMARK records containing the sequence identities, so you can choose "from pdb REMARK" in the pulldown for "Similarity of PDB..."
- Don't forget to enter a TITLE.
- When you have entered all the information, run Phaser.
- Has Phaser solved the structure? What was the LLG of the best solution? What was the Z‑score of the best translation function solution?
- The meaning of the Z‑score is given in the documentation
- Look though the log file and identify the anisotropy correction, rotation function, translation function, packing, and refinement modes. Draw a flow diagram of the search strategy.
- How many potential solutions did Phaser find or reject at each stage? What were the selection criteria for carrying potential solutions forward to the next step in the rotation and translation functions? How many other selection criteria could have been used, and what are they?
- Use the documentation
- Run Phaser again without using ensembling i.e. run two jobs, one using 1BIK only and the other using 1D0D only as models. What are the LLGs of the final solutions? What are the Z‑scores of the translation functions? Was ensembling a good idea? You could also run Phaser on the models before sculpting, to see if this has improved their quality.