EP using Phenix

From Phaserwiki
Revision as of 16:38, 20 August 2010 by Airlie (talk | contribs) (Experimental phasing with Phaser in Phenix)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

Experimental phasing with Phaser in Phenix

Data for this tutorial are found here.

We will assume that you are using version 1.6 of Phenix or later, as the graphical interface does not work as smoothly in earlier versions. Further documentation on Phenix can be found on the Phenix webpage.

Tutorial 1: SAD

SAD phasing and model-building in Phenix using HySS, Phaser and Resolve.

Reflection data: iod_scala-unique.mtz 
Sequence file: hewl.seq 

Lysozyme will readily crystallise from a wide range of conditions and yield several crystal forms. One of those, tetragonal lysozyme, is particularly well suited for halide soaking, since it grows from a high concentration of sodium chloride. A dataset has been collected from a lysozyme crystal soaked in 0.5 M potassium iodide. We will use the Phenix AutoSol wizard to find the anomalous substructure (HySS), complete the substructure and compute SAD phases (Phaser) and use a combination of map improvement and model-building to construct an initial model (Resolve).

  1. Download the data file, un-tar it and cd into the directory.
  2. Type the command "phenix" to open the Phenix GUI. Create a new project, which you could name "SAD".
  3. Open the AutoSol GUI by choosing the AutoSol button under Phasing on the right hand side.
  4. Use the "Browse" button to open the MTZ file containing the SAD data (iod_scala-unique.mtz) and the sequence file (hewl.seq). Enter the wavelength (1.5418 for Cu-Ka radiation), the atom type (I) and a guess for the number of sites from the soaking experiment (10). The guess will be used by HySS when it determines the substructure, but the subsequent refinement and completion in Phaser will add as many sites as there is good evidence for in log-likelihood-gradient maps. If you like, you can click the button "Guess missing f'/f" values" and the values will be filled in by table lookup, but Phaser will carry out a table lookup from the wavelength and atom type in any case.
  5. Click the "Run" button at the top of the window.
  6. Once phenix.xtriage has been run to check the data quality, the logfile and graphs will be available from the Summary tab. Click on "Results and graphs" to see that the data set appears to be of good quality. In particular, note that the "Anomalous measurability" graph indicates useful anomalous signal to at least 2.2Å.
  7. Once phenix.hyss has been run to determine the substructure and Phaser has been run to complete the substructure and carry out phasing, the number of sites and some quality indicators will be available under the "Heavy-atom search and phasing" tab. The "Graphs" tab will show the results of phasing at an initial low resolution then at the full resolution.
  8. After Resolve has been run to carry out a quick model-building job, the "Model-building" tab will show the statistics for the results. At this point it is probably a good idea to first go into the "Summary" tab and click "Open in Coot" to view the initial model in the density-modified map.
  9. If you have time, go back to the "Model-building" tab and click "Run AutoBuild" to launch the AutoBuild GUI for a more extensive round of model-building.

Tutorial 2: MR+SAD

This tutorial illustrates a common molecular replacement/experimental phasing scenario, when refinement is hindered by very strong model bias, but there is some experimental phasing signal available.

Reflection data: lyso2001_scala1.mtz 
Lactalbumin model: 1fkq_prot.pdb 
Sequence file: hewl.seq 

Goat α-lactalbumin is 40% identical to hen egg-white lysozyme. Although it is possible to solve lysozyme using α-lactalbumin as a model, it is very difficult to refine the structure, partly because of model bias. Unfortunately, the low solvent content of this crystal form limits the ability of density modification to remove the bias. However, one can use anomalous scattering from intrinsic sulfur atoms to improve phases dramatically. It is noteworthy that the anomalous signal from the sulfur atoms in this data set is not sufficient for ab initio phasing (it is not possible to locate the anomalous scatterers from the data alone).

If you have run the pure SAD tutorial above, you shouldn't need to set up a new project.

  1. Solve the structure with the α-lactalbumin model by running Phaser with the AutoMR wizard in Phenix. Use the "Browse" button to choose the reflection file (lyso2001_scala1.mtz). Under the "Search models" tab, define the MR model by specifying the PDB file (1fkq_prot.pdb) and its sequence identity to the target (0.4). Under the "Asymmetric unit contents" tab, specify the content by providing the sequence file (hewl.seq). Note that you are only expecting one copy in the asymmetric unit. Uncheck the box "Run AutoBuild after MR" because we will want to run AutoSol instead, so that the phases can be improved by addition of the S-SAD phase information. (If you forget, you will get a chance to cancel the AutoBuild run!) Now click "Run".
  2. This job results in a single solution with a convincing translation-function Z-score, but it would be relatively difficult (though not impossible) to rebuild and refine. We can add the S-SAD phase information easily by going to the "Summary" tab and clicking the "Run MR-SAD" button to launch the AutoSol GUI. Note that all of the required input files have already been selected. All you need to do is to specify the wavelength (1.5418Å) and the atom type (S) before clicking the Run button.
  3. You could look through the output as for the iodide test case. In particular, note that phenix.xtriage gives a much more cautious assessment of the quality of the anomalous differences.
  4. Once the "Density modification" section appears under the "Heavy-atom search and phasing" tab, you can view the S atom sites within the density-modified map by selecting solution 1 then clicking "View sites and density-modified map".
  5. Once the initial model-building has finished, as before, you can go to the "Summary" tab and click "Open in Coot" to view the initial model.
  6. If you have time, you can open the "Model-building" tab and click "Run AutoBuild" to carry out a more extensive model-building.