Difference between revisions of "Molecular replacement with MRage"

From Phaserwiki
(Input)
(Input)
Line 13: Line 13:
  
 
The following screenshot shows the this required settings to run the tutorial. Note, that almost everything is filled in automatically once the reflection file has been selected, or default. The only additional change one needs to make is to switch the '''space group exploration''' box to ''enantiomorph'' from its default value ''dataset'', as the space group is enantiomorphic, and it is not possible to know without solving the structure, which hand is the correct one. In this tutorial, only a single CPU is used, but feel free to use more if available! Note that the '''Overall count''' box has been left blank. This tells the program to determine the number of copies automatically.
 
The following screenshot shows the this required settings to run the tutorial. Note, that almost everything is filled in automatically once the reflection file has been selected, or default. The only additional change one needs to make is to switch the '''space group exploration''' box to ''enantiomorph'' from its default value ''dataset'', as the space group is enantiomorphic, and it is not possible to know without solving the structure, which hand is the correct one. In this tutorial, only a single CPU is used, but feel free to use more if available! Note that the '''Overall count''' box has been left blank. This tells the program to determine the number of copies automatically.
{| class="wikitable"
+
{| class="wikitable" width=10%
 
|[[File:Bb-basic-options_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/Basic-options.png]]
 
|[[File:Bb-basic-options_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/Basic-options.png]]
 
|-
 
|-
Line 22: Line 22:
 
First, the ''beta'' component is described with its sequence, and the only search model. Since there is no modification needed to this model (it is 100% identical), it is input through the '''Ensemble''' tab. This tells the program that there is one alternative model for the ''beta'' component.
 
First, the ''beta'' component is described with its sequence, and the only search model. Since there is no modification needed to this model (it is 100% identical), it is input through the '''Ensemble''' tab. This tells the program that there is one alternative model for the ''beta'' component.
  
{|class="wikitable"
+
{|class="wikitable" width=10%
 
|[[File:Bb-beta-composition_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/bb-beta-composition.png|centre]]
 
|[[File:Bb-beta-composition_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/bb-beta-composition.png|centre]]
 
||[[File:Bb-beta-ensembles_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/bb-beta-ensembles.png|centre]]
 
||[[File:Bb-beta-ensembles_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/bb-beta-ensembles.png|centre]]
Line 30: Line 30:
 
|}
 
|}
 
Then, a new component is created for ''blip'' (using the '''Add component''' button on the '''Basic info''' tab). This is filled out accordingly.
 
Then, a new component is created for ''blip'' (using the '''Add component''' button on the '''Basic info''' tab). This is filled out accordingly.
{|class="wikitable"
+
{|class="wikitable" width=10%
 
|[[File:Bb-blip-composition_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/bb-blip-composition.png|centre]]
 
|[[File:Bb-blip-composition_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/bb-blip-composition.png|centre]]
 
||[[File:Bb-blip-ensembles_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/bb-blip-ensembles.png|centre]]
 
||[[File:Bb-blip-ensembles_thumb.png|link=http://www-structmed.cimr.cam.ac.uk/phaser/tutorial/bb-blip-ensembles.png|centre]]
Line 80: Line 80:
 
'''Command line'''
 
'''Command line'''
 
  phaser.MRage --verbosity=VERBOSE betablip.eff
 
  phaser.MRage --verbosity=VERBOSE betablip.eff
 
  
 
=== Output ===
 
=== Output ===

Revision as of 11:48, 19 July 2013

This tutorial explains how to use to run phaser.MRage from the Phenix GUI (input files to run the same through the command line examples are also available).


Beta/blip

This is a complex between beta lactamase (beta) and beta lactamase inhibitor (blip). There are good models available, but the space group is ambiguous. This tutorial shows basic MRage functionality with a two-component system.

Experimental data for this tutorial is distributed with Phenix GUI and is also available from the Tutorials page.

Input

The MRage GUI is available from the main Phenix window under the Molecular replacement tab (select the Enable alpha features options under File/Preferences dialog for it to show up). On startup, the GUI presents the Basic options dialog.

The following screenshot shows the this required settings to run the tutorial. Note, that almost everything is filled in automatically once the reflection file has been selected, or default. The only additional change one needs to make is to switch the space group exploration box to enantiomorph from its default value dataset, as the space group is enantiomorphic, and it is not possible to know without solving the structure, which hand is the correct one. In this tutorial, only a single CPU is used, but feel free to use more if available! Note that the Overall count box has been left blank. This tells the program to determine the number of copies automatically.

Bb-basic-options thumb.png
The Basic options dialog

Next, we need to tell the program about the composition of the asymmetric unit and the available search models.

First, the beta component is described with its sequence, and the only search model. Since there is no modification needed to this model (it is 100% identical), it is input through the Ensemble tab. This tells the program that there is one alternative model for the beta component.

Bb-beta-composition thumb.png
Bb-beta-ensembles thumb.png
The Basic info dialog of Component1 (beta) The Ensembles dialog of Component1 (beta)

Then, a new component is created for blip (using the Add component button on the Basic info tab). This is filled out accordingly.

Bb-blip-composition thumb.png
Bb-blip-ensembles thumb.png
The Basic info dialog of Component2 (blip) The Ensembles dialog of Component2 (blip)

The job is ready to run!

Running from the command line

Input file (betablip.eff)

hklin = beta_blip.mtz
labin = Fobs,Sigma
symmetry_exploration = enantiomorph
composition
{
  component
  {
    sequence = beta.seq
    ensemble
    {
      coordinates
      {
        pdb = beta.pdb
        identity = 1.00
      }
    }
  }
  component
  {
    sequence = blip.seq
    ensemble
    {
      coordinates
      {
        pdb = blip.pdb
        identity = 1.00
      }
    }
  }
}

Command line

phaser.MRage --verbosity=VERBOSE betablip.eff

Output

The GUI displays the program output in the logfile window. First, the program starts searching in both possible space groups.


Bb-log-first-round-start thumb.png
Bb-log-space-group-evaluation thumb.png
Bb-log-second-round-start thumb.png
Initially, search starts in both P3₁21 and P3₂21 After finding the first component, the search in P3₁21 is dropped... ...and only P3₂21 is propagated with the partial solutions found.
Bb-log-final thumb.png
Bb-solutions thumb.png
The final solution: evaluation shows a high probability for it to be correct. Final solution shown in the GUI Top solution tab with possible further steps.

Lysozyme

As a frequently used test protein, there are good models for lysozyme. This tutorial demonstrates how to solve lysozyme starting from a homology search.

Input

Lyso-basic-options thumb.png
Lyso-composition thumb.png
Lyso-homology-searches thumb.png
Basic options dialog. Note that the space group is enantiomorphic, but we pretend we know the correct space group P4₃2₁2 to save time Setting up composition. It would be possible to tick the NCBI Blast option and run, but we want to avoid overloading the NCBI and provide a saved homology search file Specifying the homology search file. There are many hits, and the GUI by default selects the first three. Selecting any larger number, however, makes no difference!
Running from the command line

Input file (lysozyme.eff)

hklin = lyso.mtz
labin = "F_New,SIGF_New,DANO_New,SIGDANO_New,ISYM_New"
composition
{
  component
  {
    sequence = hewl.seq
    homology
    {
      file_name = hewl_ebi_wu_blast.xml
      max_hits = 3
    }
  }
}

Command line

phaser.MRage --verbosity=VERBOSE lysozyme.eff

Output

Lyso-fetch-1lsg thumb.png
Lyso-sculptor-1lsg thumb.png
Lyso-finish thumb.png
The program fetches the first homologue... ...then processes it with sculptor. Note that although all 13 protocols in sculptor are active, it only creates 7 models, because remaining protocols create models that are identical to the ones used. This model gives a clear solution, and processing stops. Although 3 hits were specified, the second and the third was never downloaded, because the first gave a clear solution.