Xsight



5       Tutorial


Overview of tutorial lessons

The purpose of these lessons is to introduce you to the commands and functionalities of the Xsight module of the Insight II molecular modeling program. You will perform several common procedures normally performed in macromolecular crystallographic analysis. This should acquaint you with the most commonly used features and calculation strategies of Xsight. Lessons illustrating MIR phasing (Isomorphous replacement), Phase determination from MAD data, Density modification, Density fitting, Molecular replacement, Refinement, Completing the atomic model and Visualization are included in this set of tutorials.

Throughout the tutorial sessions that involve XtalView modules, the following conventions are observed when discussing mouse actions:

Hardcopy tutorial lessons

All of the Xsight tutorial lessons are provided in this user guide (there are no pilot tutorials).

This user guide contains the following lessons:

Isomorphous replacement

Lesson 1:Merging native and derivative data
Lesson 2:Finding heavy atom sites for the first derivative
Lesson 3:Adding a second derivative

Phase determination from MAD data

Lesson 4:MAD phase determination

Density modification

Lesson 5:Solvent flattening and histogram matching
Lesson 6:Non-crystallographic symmetry averaging

Density fitting

Lesson 7:Making maps and basic fitting techniques
Lesson 8:De novo map fitting

Molecular replacement

Lesson 9:Calculating a self-rotation function
Lesson 10:Generating structure factors for the search model
Lesson 11:Ordinary cross-rotation function
Lesson 12:Translation search
Lesson 13:Applying the solutions

Refinement

Lesson 14:Rigid-body refinement
Lesson 15:Multi-rigid-body refinement
Lesson 16:Simulated annealing refinement
Lesson 17:Least-squares minimization

Completing the atomic model

Lesson 18:Generating ligand dictionaries
Lesson 19:Locating ordered water molecules

Visualization

Lesson 20:Richardson diagrams
Lesson 21:Display of atomic models with electron density maps
Lesson 22:Display and creation of symmetry related molecules


Setting up to run tutorials

Before starting to use Xsight, you must carry out the following setup steps. For a more detailed description of the directory structure referenced by the $CRYSTALDATA environment variable, see Chapter 3, Implementation.

1. Create a directory that will be referenced by the $CRYSTALDATA environment variable by typing the UNIX command:


>	mkdir $HOME/crystal_info
2. Copy the crystals and projects, and the cvccp files, from the $XSIGHTTUTORIAL directory to the $CRYSTALDATA directory:


>	cp $XSIGHTTUTORIAL/crystals $HOME/crystal_info

>	cp $XSIGHTTUTORIAL/projects $HOME/crystal_info

>	cp $XSIGHTTUTORIAL/cvccp $HOME/crystal_info
3. Make sure that your crystals and projects files are writeable.


>	chmod +w $HOME/crystal_info/crystals

>	chmod +w $HOME/crystal_info/projects
4. Edit your .cshrc file to include two environment variables:


setenv CRYSTALDATA $HOME/crystal_info
setenv CRYSTAL cvccp
5. Log out and log back in.

Each time you start a new project or create a new crystal, the information will be updated in the $CRYSTALDATA directory. The cvccp crystal file is used for many of the tutorials. You will want to change your $CRYSTAL environment variable, depending on which crystal you most frequently use.

Do not do your work in the $CRYSTALDATA directory. This directory is only for storing information about crystals and projects. Define your working directory through the Utilities/Project command in the Xsight interface.


Isomorphous replacement

Isomorphous replacement is one of the key techniques used in solving protein crystal structures. The following lessons lead you through the process of merging native and derivative data, calculating Patterson maps, finding heavy atom sites, and combining phases from multiple derivatives.

In Lesson 1:Merging native and derivative data, you learn to convert processed data into the .fin file format. You then merge and scale a native data set with an Au derivative. Finally, the lesson takes you through the different statistical plots that you can generate to evaluate the quality of the derivative.

In Lesson 2:Finding heavy atom sites for the first derivative, you calculate and contour a difference Patterson map for the Au derivative. You then locate the two major Au sites using an automated search method. You verify these positions against the difference Patterson map. Finally, you refine the Au sites and calculate SIR phases.

In Lesson 3:Adding a second derivative, you merge and scale an Ir derivative with the native data set. You then calculate a cross-Fourier using the SIR phases from the Au derivative and the isomorphous differences between the Ir data and the native data. From this map you locate an Ir site in the Ir derivative crystal. Finally, you refine the Ir site and calculate the MIR phases using the Au and Ir positions.


Lesson 1:   Merging native and derivative data

This lesson uses the following commands, accessed from the Xsight module:

Xsight module

Utilities

Data_Control

MIR

Project

Import_Data

Merge_Data

Statistics

The Scenario: You have just collected a putative derivative data set of C. vinosum Cytochrome c' (cvccp) and you want to merge it with an existing native data set. You then want to calculate statistics that will indicate whether you have been successful in producing a derivative.

1.   Preparing the tutorial directory

In order to run this tutorial (and the other associated MIR tutorials), you must set up a directory which contains the ccp_raxis_au1.dat, ccp_raxis_ir1.dat, and ccp_raxis_n1.dat files, that is, two derivative data sets and a native data set.

Enter the following commands, in the order shown, at the command prompt:

>	mkdir tutorial_mir

>	cd tutorial_mir

>	cp $XSIGHTTUTORIAL/*raxis* . 

The following steps are only necessary if you do not already have a cvccp crystal file. If this is the first Xsight tutorial that you have run, then you probably need to enter these commands. If you have already completed these steps once before, there is no need to repeat them.

Use your favorite text editor to add the following line to the crystals file:

>	cvccp
Alternatively, you can add the line with the cat command by typing the following lines at the UNIX command prompt:

>	cd $CRYSTALDATA

>	cp $XSIGHTTUTORIAL/cvccp .

>	cat >> crystals
Type cvccp in the textport, even though there is no UNIX prompt, then press <Ctrl-D> to exit the cat command mode.

The >> symbol allows you to append lines to an existing file. Once you type the cat command, all typed entries in the textport are appended to the crystals file. For this lesson, you simply want to add the cvccp crystal to the file and exit from the cat append mode. Be sure not to type anything else in the textport while you are in the cat command mode.

2.   Getting started

Invoke the Insight II software.

At the system prompt issue the command:

>	insightII
and press <Enter>.

Wait a few moments while Insight II loads.

3.   Creating a project

Go to the Xsight pulldown by clicking the MSI logo and selecting Xsight from the list that appears (the Xsight module commands are now displayed on the lower menu bar).

Select the Utilities/Project command and type in the name of the project (for example, mir).

Toggle on Edit.

Click the Crystal parameter field, and select cvccp from the Crystal_List, then enter the name of the Project_Directory you created in Step 1. (You must use the absolute path for this directory, since ~usr is not accepted in this field--for example: /net/machine/usr/people/pns/tutorial_mir.)

Select Execute.

The dialog box disappears automatically.

4.   Importing and converting a processed data file

This lesson assumes that you have your processed reflection files in a format other than the .fin reflection format. In particular, it assumes you are starting with three reflection files in the standard R-AXIS II format. In this lesson, you use the native data set and one of the derivative data sets. The second derivative data set is used in Lesson 3:Adding a second derivative.

First, you will convert the native data set to the .fin file format.

Select the Data_Control/Import_Data command

Select the hkl_F1_SIG1_DELANOM option for the Input_File_Format.

Enter ccp_raxis_n1.dat in the Input_File_Name parameter field.

Select Fin(.fin) in the Output_File_Format parameter field.

Enter ccp_n1.fin in the Output_File_Name parameter field.

Select Execute.

5.   Converting the Au derivative data set

Now you will convert the Au derivative data set to a .fin format.

Select the Data_Control/Import_Data command.

Select the hkl_F1_SIG1_DELANOM option for the Input_File_Format.

Enter ccp_raxis_au1.dat in the Input_File_Name parameter field.

Select Fin(.fin) in the Output_File_Format parameter field.

Enter ccp_au1.fin in the Output_File_Name parameter field.

Select Execute.

6.   Merging the native and Au derivative data files

You now want to scale and merge the native and Au derivative data files so you can evaluate the quality of the derivative.

Select the MIR/Merge_Data command.

Enter ccp_n1.fin as the Fin_File_1 parameter.

Enter ccp_au1.fin as the Fin_File_2 parameter.

Enter ccp_n1.ccp_au1.fin as the Merge_Output_File parameter.

Select Execute.

This displays the Xmerge window

In the Xmerge menu, set Scaling Type to Anisotropic.

Click Scale to begin the merging procedure.

After a few moments, two graphs are generated at the end of the merging procedure: merging R-factor versus resolution and the average |Fp-Fph| versus resolution. These two graphs help you evaluate your derivative data.

As a general rule of thumb, a merging R-factor above 10% indicates that you have a derivative crystal with significant phasing power. An R-factor above 7% may have significant phasing power. An R-factor below 7% indicates a crystal that cannot be used for phasing.

The |Fp-Fph| differences should decrease as a function of resolution. If they do not, it is a symptom of non-isomorphism. In addition, the centric differences should be larger than the acentric differences.

Select Quit, then click Yes to confirm.

7.   Generating more graphs

This is an optional step in which you may replot the comparisons between the native and derivative data.

The graphs that are generated from the Xmerge program are not optimal because reflections are not evenly distributed into resolution bins. A better analysis can be made using the Xstat program.

Select the MIR/Statistics command.

Set Input_File to ccp_n1.ccp_au1.fin.

Select Execute.

The Xfinstat window now appears.

In the Xfinstat menu toggle on R-factor vs. resolution and |F1-F2| vs. resolution.

Click Graph it.

After a few moments, two graphs are calculated and generated using an even reflection distribution between bins. These plots can be easily sent to a PostScript printer, if you desire.

At the top of either graph, there is a parameter field labeled Printer/File. Enter the command that you normally use to print a PostScript file, except do not include the filename.

Click Print.

The graphs can be deleted by double-clicking the top left button of each graph window.

Now exit the Xstat program.

Select Quit, then click Yes to confirm.

8.   Ending the tutorial

You can now continue to the next tutorial to locate and refine the heavy atoms of the first derivative, or you can quit.

To quit Insight II, type quit on at the command line, and press <Enter>.


Lesson 2:   Finding heavy atom sites for the first derivative

This lesson uses the following commands, accessed from the Xsight module:

Xsight module

Utility

MIR

Project

Calc_Fourier

Contour_Map_3D

Find_Heavy_Atoms

Calc_Phases

You need to have run lesson 1 before doing this tutorial lesson.

The Scenario: In Lesson 1:Merging native and derivative data, you merged an Au derivative data set of C. vinosum Cytochrome c' (cvccp) with a native data set. The intensity statistics indicated a derivative with good phasing power. In this lesson, you will locate and refine the positions of the two major heavy atom sites for this first derivative.

1.   Starting the tutorial

If you have just finished Lesson 1:Merging native and derivative data, and you are still running Insight II, you can skip to step 3.

If you have run Lesson 1:Merging native and derivative data, but need to restart Insight II carry out step 2.

2.   Getting started

At the system prompt issue the command:

>	insightII
and press <Enter>.

Insight II takes a few moments to load.

Go to the Xsight pulldown by clicking the MSI logo and selecting Xsight from the list that appears.

Select the Utilities/Project command and enter mir as the Project_Name parameter.

Select Execute.

3.   Calculating a difference Patterson map

In the last lesson you created a data file that contains scaled and merged reflections from a native data set and an Au derivative data set. You now want to calculate a Patterson map using coefficients based on the differences between the two data sets.

Select the MIR/Calc_Fourier command.

Enter ccp_n1.ccp_au1.fin for the Phase_File parameter.

Make sure that ccp_n1.ccp_au1.map is entered as the Map_File parameter.

Select Execute.

Since you are calculating a Patterson map, you need to generate the Fourier coefficients from the observed differences between the native and derivative data sets. Your input file to the Fourier program therefore comes from a .fin file (structure factor amplitude file) rather than a .phs file (phase angle file).

4.   Start the Patterson calculation

When the Xfft menu appears you will need to change a few of the default parameters. You should type <Enter> after altering the numerical values in order to ensure that your changes are properly entered.

Set the following parameter values in the Xfft menu:

Map Type: Fo*Fo (Patterson)

Resolution filter: 999.999

to: 5.0

Outlier Filter: differences> 100

Using MB3, click and hold the Phase File Type button. Drag the cursor to the Extra for Pattersons only selection and then drag to the right. (An additional pulldown menu appears.)

Select the fin(Fp s(Fp) Fph s(Fph)) option.

The label Patterson map from h,k,l,fp,s(fp),fm,s(fm) should now be loaded next to the Phase File Type button.

The other parameters should remain at their default values.

Click Calculate to create the difference Patterson map.

The calculation takes a few moments. The message FFT finished tells you when it is completed.

Click Quit when you are finished, then click Yes to confirm.

5.   Display the difference Patterson map

Now you want to contour the Patterson map and look at the Harker sections for significant peaks. These peaks correspond to self-vectors between crystallographically related copies of bound heavy atoms and confirm the bounding of the heavy atoms in the derivative crystal.

Select the MIR/Contour_Map_3D command

Check that Contour_Operation is set to Setup.

Set the Map_File_Type to Patterson_Map.

Set Map_File to ccp_n1.ccp_au1.map.

Check that the Contour_Sigma_Level is set to 3.0, the Two_Level_Contour option is toggled on, the High_Sigma_Level is set to 6.0 and the Input_Atom_Sites option is toggled off.

Select Execute.

After a few moments the three-dimensional representation of the Patterson map appears.

The space group for this structure is P212121, and therefore the three Harker sections occur at the x = 1/2 section, the y = 1/2 section, and the z = 1/2 section. In this display the Harker sections are represented by sets of green grid lines. The appearance of a large slightly non-spherical peak on each of the Harker sections in the Patterson map suggests that there are probably one or two major Au sites in the crystal asymmetric unit.

Now that you have confirmed the presence of bound Au atoms you may leave this command

Set Contour_Operation to Quit_Command.

Select Execute.

6.   Locate the heavy atoms

Select the MIR/Find_Heavy_Atoms command.

Set Input_Fin_File to ccp_n1.ccp_au1.fin.

Make sure that Input_Data_Type is set to Isomorphous_Diff, Output_Sol_File is set to ccp_n1.ccp_au1.sol, Scattering_Element is set to AU, Low_Res_Cutoff is set to 10.0 and High_Res_Cutoff is set to 5.0.

The Search Limits parameters should be set so that Start_X is 0.0, Stop_X is 0.5, Start_Y is 0.0, Stop_Y is 0.5, Start_Z is 0.0 and Stop_Z is 1.0.

Make sure that Set_Search_Type is set to Double_Site.

Set Min_Cor_On_I to 0.1.

Select Execute.

These input parameters set up a two-site heavy atom search using data to 5 Å resolution over a unique volume in the crystal cell. The Au sites that are located will be written to a solution file called ccp_n1.ccp_au1.sol.

You will now run the search program.

Make sure that Heavy_Operation is set to Run.

MIR_Run_Mode
should be set to Run_Now.

Select Execute.

The program run takes about 15 minutes on an SGI R4400. Diagnostic information in the log file called find_heavy.log will indicate that two sites were found.

When the Notifier message appears to inform you that the background job is complete, select Continue.

7.   Checking the heavy atom solution

You should now confirm that the predicted sites are consistent with the difference Patterson function.

Select the MIR/Contour_Map_3D command.

Make sure that Map_File_Type is set to Patterson_Map.

Make sure that the Map_File parameter contains ccp_n1.ccp_au1.map.

Set Input_Atom_Sites to on.

Set In_Solution_File to ccp_n1.ccp_au1.sol.

Select Execute.

The graphical display of the difference Patterson map will appear. A table reports the positions of the two Au atoms that were located by the heavy atom search program. A pickable table will report the interatomic vectors that should correspond to peaks in the difference Patterson map.

When you have checked that the self-vectors (yellow) and cross-vectors (purple) lie in or are very close to significant Patterson density you may leave this command.

Set Contour_Operation to Quit_Command.

Select Execute.

8.   Starting the heavy atom refinement and phase calculation program

Select the MIR/Calc_Phases command.

Select ccp_n1.ccp_au1.sol for the Solution_File parameter.

Enter ccp_n1.ccp_au1.phs for the Phase_Output parameter.

Select Execute.

The heavy atom search program used in the previous step wrote the derivative name `AU' and the name of the .fin file containing the native and derivative data into the solution file header. Therefore, this information will be passed automatically to the Xheavy program. You only need to use the derivative editor to select the resolution limits for the data that will be used in the heavy atom refinement

Use MB1 to select the derivative AU.

Click Edit.

The Xheavy Derivative Editor will appear.

In the Xheavy Derivative Editor menu change the high Resolution parameter from 1.0 to 5.0.

Select Apply.

Next you will refine the two sites. By default, Method in the Xheavy menu is set to Refine Selected Derivative Only. Since there is only one derivative, this option is the correct one.

In the main Xheavy menu select Apply.

The two Au sites are refined by the correlation method. This calculation takes about half a minute on an SGI R4400. When five refinement cycles are completed the correlation coefficient is reported as 0.8214.

Since the refinement did not completely converge you will want to run further refinement cycles.

In the main Xheavy menu select Apply again

After three refinement cycles the refinement converges with a correlation coefficient of 0.8386.

You will now want to save the refined parameters.

Click the Save Derivative button, then click Overwrite at the prompt.

9.   Calculating SIR phases

You can now calculate SIR phases based on the refined heavy atom sites.

Using MB3, click and hold the arrow button next to Method.

Select Calculate Protein Phases and release MB3.

Click Apply.

The calculated phases are now stored in the Output Phases file (ccp_n1.ccp_au1.phs).

Now exit the Xheavy program.

Select Quit, then click Yes to confirm that you want to exit.

10.   Ending the tutorial

You can now continue to the next tutorial to locate and refine the heavy atom of the second derivative, or you can quit.

To quit Insight II, type quit on the command line, and press <Enter>. Press <Enter> again at the next prompt.


Lesson 3:   Adding a second derivative

This lesson uses the following commands, accessed from the Xsight modules:

Xsight module

Utilities

Data_Control

MIR

Project

Import_Data

Merge_Data

Calc_Fourier

Contour_Map_3D

Calc_Phases

Merge_Phases

You need to have run lessons 1 and 2 before doing this tutorial.

The Scenario: You have just collected a second derivative data set of C. vinosum Cytochrome c' (cvccp) and you want to merge it with an existing native data set. You then want to calculate statistics that will indicate whether you have been successful in producing a derivative. Secondly, you want to calculate a cross-Fourier map to locate the heavy atom position. Finally, you want to combine the phasing information from both derivatives in order to calculate MIR phases.

1.   Starting the tutorial

If you have just finished Lesson 1:Merging native and derivative data, and you are still running Insight II, you can skip to step 3.

If you have previously run Lesson 2:Finding heavy atom sites for the first derivative but have left Insight II go to step 2.

2.   Getting started

At the system prompt issue the command:

>	insightII
and press <Enter>.

Wait a few moments while Insight II loads.

Go to the Xsight pulldown by clicking the MSI logo and selecting Xsight from the list that appears.

Select the Utilities/Project command and enter mir as the Project_Name parameter.

Select Execute.

3.   Importing a processed data file

This lesson assumes that you have your processed reflection files in a format other than the .fin reflection file format. In particular, it assumes you are starting with a reflection file that is in the standard R-AXIS II format.

Select the Data_Control/Import_Data command.

Select the hkl_F1_SIG1_DELANOM option for the Input_File_Format.

Enter ccp_raxis_ir1.dat in the Input_File_Name parameter field.

Select Fin(.fin) in the Output_File_Name parameter field.

Enter ccp_ir1.fin in the Output_File_Name parameter field.

Select Execute.

4.   Merging the native and Ir derivative data files

You now want to scale and merge the native and Ir derivative data file, so that you can evaluate the quality of the derivative.

Select the MIR/Merge_Data command.

Select ccp_n1.fin from the value-aid for the Fin_File_1 parameter.

Select ccp_ir1.fin from the value-aid for the Fin_File_2 parameter.

Enter ccp_n1.ccp_ir1.fin as the Merge_Output_File parameter.

Select Execute

This displays the Xmerge window.

In the Xmerge menu, set Scaling Type to Anisotropic.

Click Scale to begin the merging procedure.

Two graphs are generated at the end of the merging procedure: one merging R factor versus resolution, and the other showing average |Fp-Fph| versus resolution. These two graphs help you evaluate your derivative data.

As a general rule of thumb, a merging R factor above 10% indicates that you have a derivative crystal with significant phasing power. An R factor above 7% may have significant phasing power. An R factor below 7% indicates a crystal that cannot be used for phasing.

The |Fp-Fph| differences should decrease as a function of resolution. If they do not, this is a symptom of non-isomorphism. In addition, the centric differences should be larger than the acentric differences.

Now exit the Xmerge program.

Select Quit in the Xmerge menu, then click Yes to confirm.

5.   Merging the Au phases with the Ir and native .fin file

In order to use a cross-Fourier map to locate the Ir site, you must first merge the SIR phases calculated from the Au derivative with the new merged Ir and native data.

Select the MIR/Merge_Phases command.

Select ccp_n1.ccp_ir1.fin from the value-aid for the Fin_File_1 parameter.

Select ccp_n1.ccp_au1.phs from the value-aid for the Phase_File parameter.

Enter ccp_ir1.ccp_au1.phs in the Phase_Output parameter field.

Select Execute.

The ccp_n1.ccp_ir1.fin file that was created in an earlier step has the native data stored first. In order to locate the Ir position as a positive peak, you need to subtract the native structure factor amplitudes from the Ir structure factor amplitudes. This requires a reversal of F values when the new .phs file is generated.

Toggle on the Swap f1 and f2 (isomorphous Fourier) check box and click Merge.

The creation of the new .phs file will only take a few seconds.

Select Quit, then click Yes to confirm.

6.   Calculating a cross-Fourier map

You can now calculate a cross-Fourier map, which when contoured should reveal the position of the major Ir site with the origin properly set relative to the Au sites.

Select the MIR/Calc_Fourier command.

Select ccp_ir1.ccp_au1.phs from the value-aid for the Phase_File parameter.

Make sure that ccp_ir1.ccp_au1.map is entered in the Map_File parameter field.

Select Execute.

7.   Calculating the Fourier using the Xfft menu

When the Xfft menu appears you will need to change a few of the default parameters. You should type <Enter> after altering the numerical values in order to ensure that your changes are properly entered.

Set the following parameter values in the Xfft menu:

Map Type: mFo-Fc
Resolution filter: 999.999
to: 5.0
Outlier Filter: differences> 100

You switched the position of the native and Ir data when you created the .phs file in the last step. That is why an Fo-Fc map will generate a positive peak for an Ir site.

The other parameters should remain at their default values.

Click the Calculate button to create the cross-Fourier map.

Select Quit when you are finished, then click Yes to confirm.

8.   Display the cross-difference Fourier map

Now you want to use the cross-Fourier map to look for significant peaks that indicate Ir sites.

Select the MIR/Contour_Map_3D command.

Set Map_File_Type to Difference_Map.

In the Map_File parameter enter ccp_ir1.ccp_au1.map.

Make sure that Contour_Sigma_Level is set to 4.0, Two_Level_Contour is toggled on and Second_Sigma_Level is set to 8.0.

Set Input_Atom_Sites to off.

Select Execute.

You will see one large peak (and its symmetry related copies) that exceeds the 8 sigma level in the display of the cross-difference Fourier map. You will pick an atom at this site and write out a solution file containing the co-ordinates.

Make sure the Site_Edit_Operation is set to Peak_Search.

Use the value-aid to change the Scattering_Element to IR.

Change the Site_Peak_Threshold value to 8.0.

Select Execute.

The unique Ir site will be identified by a yellow sphere and the sites related by the crystallographic symmetry are marked as purple spheres. The site co-ordinates are also be reported in a table.

Now you will write a solution file containing this site and leave the Contour_Map_3D command.

Set Site_Edit_Operation to Write_Solution_File.

Enter ccp_n1.ccp_ir1.sol in the Out_Solution_File parameter block.

Select Execute.

Set Contour_Operation to Quit_Command.

Select Execute.

9.   Starting the heavy atom refinement and phase calculation program

Select the MIR/Calc_Phases command.

Use the value-aid to enter ccp_n1.ccp_ir1.sol as the Solution_File parameter.

Enter ccp_mir.phs as the Phase_File parameter.

Select Execute.

When you analyzed the cross-difference Fourier map for Ir sites the derivative name IR was written into the header of the solution file. You now need to use the Xheavy derivative editor to enter the name of the structure factor file containing the native and derivative data and to control the resolution range for the heavy atom refinement.

Use MB1 to select the derivative IR.

Click Edit.

The Xheavy Derivative Editor appears.

In the DataFile field enter ccp_n1.ccp_ir1.fin.

Change the high Resolution parameter from 1.00 to 5.00.

Select Apply.

Now you will refine the Ir site using data to 5Å resolution. By default, Method in the Xheavy menu is set to Refine Selected Derivative Only. Since there is only one derivative, this option is the correct one.

Select Apply in the main Xheavy menu.

The Ir site is refined by the correlation method.

Now refine the Ir site using data to 3 Å resolution.

In the Xheavy Derivative Editor menu change the high Resolution parameter from 5.00 to 3.00.

Select Apply.

In the main Xheavy menu select Apply.

You will now save the refined Ir parameters.

Click the Save Derivative button, then click Overwrite at the prompt.

10.   Loading the Au derivative file and calculating MIR phases

Now you load the Au derivative.

In the main Xheavy menu enter the filename, ccp_n1.ccp_au1.sol, in the Derivative File parameter field.

Click the Load Derivative button.

The Au derivative needs to be refined using 3.0 Å data (previously only 5.0 Å data was used to refine it).

Select the AU derivative from the derivative list in the main Xheavy menu using MB1, then click Edit.

Now you will use the Xheavy Derivative Editor menu.

In the Xheavy Derivative Editor menu change the high Resolution parameter from 5.00 to 3.00.

Click Apply.

Next you return to the main menu in order to refine the derivative. The Au derivative should still be selected in the scrolled derivative list. The Method should still be set to Refine Selected Derivative Only.

In the main Xheavy menu, click Apply.

This step takes about a minute. At the end of the refinement the correlation coefficient is reported as 0.7742.

You will now save the refined Au parameters.

Click the Save Derivative button, then click Overwrite at the prompt.

11.   Calculating MIR phases

You can now calculate MIR phases based on the refined heavy atom sites.

Use MB3 to click and hold the Method arrow button. Select the Calculate Protein Phases button and release MB3.

Click Apply.

After about one minute, the calculation of the protein phases will be completed. Because of the limited amount of heavy atom derivative information this phase set is relatively poor, with a figure of merit of 0.497. The calculated phases are stored in the Output Phases file, ccp_mir.phs.

Note: The MIR map produced in this tutorial is enantiomorphic to the correct MIR map because of an enantiomorph ambiguity inherent in heavy atom site determination. In an interpretable MIR map this ambiguity is most easily recognized by noting the wrong handedness of density corresponding to alpha helices. The correct MIR map is obtained by changing the enantiomorph of the heavy atoms sites and repeating the final phasing calculation. The Contour_Map_3D command contains an Invert_Sites option in the Sites_Edit_Operation list that is accessed for the Difference_Map display mode for carrying out the enantiomorph inversion.

12.   Ending the tutorial

Exit the Xheavy program.

Select Quit in the Xheavy menu, then click Yes to confirm that you want to exit.

If you now wish to, quit Insight II.

To quit Insight II, type quit on the command line, and press <Enter>. Press <Enter> again at the prompt.




Last updated November 16, 1998 at 11:59AM Pacific Standard Time.
Copyright © 1998, Molecular Simulations, Inc. All rights reserved.