Resonance Assignment/Abacus/FMCGUI commands: Difference between revisions
Jump to navigation
Jump to search
Data Menu
No edit summary |
|||
| (23 intermediate revisions by one other user not shown) | |||
| Line 13: | Line 13: | ||
<div></div><div> :: Save the current state of the project</div><div> </div><div>What is currently in the computer memory is saved in the file PROJECTNAME.prj located in the root directory of the project. </div><div> </div> | <div></div><div> :: Save the current state of the project</div><div> </div><div>What is currently in the computer memory is saved in the file PROJECTNAME.prj located in the root directory of the project. </div><div> </div> | ||
=== '''<span>Poject>Quit</span>''' === | === '''<span>Poject>Quit</span>''' === | ||
<div>:: Save the current state of the project and quit. | <div> </div><div>:: Save the current state of the project and quit.</div><div> </div> | ||
== Data Menu<br> == | == Data Menu<br> == | ||
<div>The DATA section serves to load&save the data (such as protein sequence and peak lists). Since there are different formats of data-files that could be loaded in memory or saved on disk, one can use this section as format converter as well.</div><div> </div> | <div>The DATA section serves to load&save the data (such as protein sequence and peak lists). Since there are different formats of data-files that could be loaded in memory or saved on disk, one can use this section as format converter as well.</div><div> </div> | ||
| Line 239: | Line 240: | ||
=== '''Structure>Constraints>H-bonds>Specify''' === | === '''Structure>Constraints>H-bonds>Specify''' === | ||
<div>''':''': Specify H-bond distance constraints manually.</div><div> </div><div>''Prerequisites'': </div><div> ♦ protein sequence is loaded in memory</div><div> </div><div>[[Image:Fmcgui fig2.16.jpg|thumb|left|100px]]A new window “O/HN Pairs” is opened (see Figure | <div>''':''': Specify H-bond distance constraints manually.</div><div> </div><div>''Prerequisites'': </div><div> ♦ protein sequence is loaded in memory</div><div> </div><div>[[Image:Fmcgui fig2.16.jpg|thumb|left|100px]]A new window “O/HN Pairs” is opened (see Figure ). </div><div> </div><div> </div><div></div><div> </div><div> <br> </div><div> </div><div> </div><div></div><div> </div><div></div><div>For each H-bond constraint, user have to specify O-HN pair by typing in the ID of residues corresponding to O and HN atoms. Pressing “OK” button will save H-Bond constraints on disk in two formats: files 'hbond.upl' and 'prot_hbond.tbl' for CYANA and CNS calculations, respectively. Both files are saved in the project root project directory.</div><div> <br></div><div> <br> </div><div></div> | ||
=== '''Structure>Constraints>H-Bonds>load''' === | === '''Structure>Constraints>H-Bonds>load''' === | ||
<div>''':''' | <div>'''::''' Load HB-ond distance constraints in CYANA format (upl-file) </div><div> </div><div> </div> | ||
=== '''Structure>Calculate>Cyana''' === | === '''Structure>Calculate>Cyana''' === | ||
<div>''':''' | <div>''':''': Set up new structure calculation run with CAYANA.</div><div> </div><div>'' Prerequisites'': </div><div> ♦ protein sequence is loaded in memory</div><div> ♦ assigned chemical shifts are loaded in memory</div><div> ♦ N15-NOESY, C13-NOESY, and Arom-NOESY peak lists are loaded in memory</div><div> ♦ specified tolerances</div><div> ♦ dihedral angle constraints are created ( i.e. file “belok.aco” is present in the project root directory)</div><div> ♦'' Optional:''</div><div> ♦ HBond distance constraints (file “hbond.upl” is present in the project root directory)</div><div> ♦ file that contains ZN ion ligands (file “zn_ligands” is present in the project root directory)</div><div> </div><div>A new directory under the name that is specified by user (normally “crun#”) will be created inside project root directory PROJECTNAME. This directory contains all files required to start automatic structure calculations with CYANA.</div><div> </div> | ||
=== '''Structure>Calculate>ABCUS''' === | === '''Structure>Calculate>ABCUS''' === | ||
<div>''':''' | <div>''':''': Sset up ABACUS structure calculations</div><div> </div><div> </div> | ||
=== '''Structure>RPF>RP''' === | === '''Structure>RPF>RP''' === | ||
<div>''':''' | <div>''':''': Perform Recall/Precision analysis of structural ensemble.</div><div> </div><div> ''Prerequisites'': </div><div> ♦ protein sequence is loaded in memory</div><div> ♦ assigned shemical shifts are loaded in memory</div><div> ♦ N15-NOESY, C13-NOESY, and Arom-NOESY peak lists are loaded in memory</div><div> ♦ specified tolerances</div><div> ♦ coordinates of structural ensemble in CYANA format (final.pdb) or in CNS format <br> </div><div> </div> | ||
<br> | ''Results:''<br> | ||
<div><span> </span></div><div> </div><div> </div><div> | |||
♦ summary report<br> | |||
<div> ♦ flase negative and false positive peak lists<span> </span></div><div> </div><div> </div><div> The parameters to set up (see Figure ):</div><div>[[Image:Fmcgui_fig2.17.jpg|thumb|left|400px]] <br></div> | |||
*<span><span> </span></span>“RP directory name”. Normally the name is rp#. A new directory under this name will be created within PROJECTNAME/assign directory. The results of the RPFanalysis will be stored in this directory. | |||
*<span><span> </span></span>“sequence gap”. Residue pairs separated by less than the value of sequence gap will be excluded from generating expected peak lists. | |||
*<span><span> </span></span>“cutting distance for recall”. Distance threshold for evaluating matching of an experimental peak to a structural ensemble (Recall score) | |||
*<span><span> </span></span>“cutting distance for precition”. Distance threshold for generating expected peak from structural ensemble (Precision score) | |||
<div> </div><div>The results of the RPF analysis will be saved in a new directory (for example, directory "rp1" on Figure) that is located in the project root directory. Results include peak lists in the SPARKY format of both false negative and false positive peaks for different NOESY spectra in a separate files.</div><div> </div> | |||
=== '''Structure>RPF>DP''' === | === '''Structure>RPF>DP''' === | ||
:''':''' Set up calculations of DP score with AutoStructure | |||
<div> </div><div>''Prerequisites'': </div><div> | <div> </div><div>''Prerequisites'': </div><div> ♦ protein sequence is loaded in memory</div><div> ♦ assigned shemical shifts are loaded in memory</div><div> ♦ N15-NOESY, C13-NOESY, and Arom-NOESY peak lists are loaded in memory</div><div> ♦ specified tolerances</div><div> ♦ coordinates of structural ensemble in CYANA format (final.pdb) or in CNS format </div><div> </div><div> </div><div> </div> | ||
=== '''Structure>Water refinement>calculate''' === | === '''Structure>Water refinement>calculate''' === | ||
<div> ''':''' | <div> ''':''': Set up water refinement calculations</div><div> </div><div>''Prerequisites'': </div><div> ♦ structure calculation with CYANA should be done</div><div> ♦ ''Optional:''</div><div> ♦ RDC data in PALES format.</div><div> ♦ file with ZN ion ligands (file [[FMCGUI_commands#Structure.3EAdd_ZN_ligands|“zn_ligands”]] )</div><div> <br></div><div> [[Image:Fmcgui fig2.18.jpg|thumb|right|400px]]<br> </div> <div>User will be asked to :</div> | ||
< | * specify the name of directory for water refinement calculations, WATDIR; | ||
<div> | * select a number of files with coordinates and constraints; | ||
* indicate cisProline residues (if there are any) | |||
* specify protonation state of HIS residues (which is double protonated by default). <br> | |||
<div> </div><div> </div><div></div><div>In the result, a new directory WATDIR will be created inside project root directory. I will contain all files and scripts required to carry out water refinement calculations with CNS. </div><div>The calculations should be started outside FMCGUI. (linux cluster is is recommended)</div><div> <br> </div><div> </div><div> <br> </div><div> </div><div> </div><div></div><div></div><div><br></div><div> </div> | |||
=== '''Structure>Water refinement>summary''' === | === '''Structure>Water refinement>summary''' === | ||
<div>''':''' | <div>'''::''' Create a summary for water refined structural ensemble.</div><div> </div><div>''Prerequisites'': </div><div> ♦ water refinement with CNS should be done.</div><div> </div><div>User will be asked to specify residues used for structure superposition and to select the water refinement directory, WATDIR (see Figure).</div><div> [[Image:Fmcgui_fig2.19.jpg|thumb|right|400px]]</div><div> <span> </span></div><div> </div><div>The refined structural models will be superimpose and combined in one file. The created summary reports values of different energy components and constraint violation statistic for each structural model. A new directory WATDIR_results will be created. The directory contains final superimposed coordinates, distance and dihedral angle constrains in a format suitable for PDB deposition, summary and constraint violations report.</div><div> </div> | ||
=== '''Structure>Add ZN ligands''' === | === '''Structure>Add ZN ligands''' === | ||
<div>''':''' | <div>''':: '''Create "zn_ligands" file.</div><div> </div><div>[[Image:Fmcgui fig2.20.jpg|thumb|left|100px]]“ZN ligands” window pops up. User have to type in IDs of residues that ligate zinc ion. If there are a few zinc ions, information for each ion should be provided in a separate row (see Figure).</div><div><span> </span></div><div> </div><div> </div><div> </div><div>The file “zn_ligands” will be created inside the project root directory by pressing “OK” button.</div><div> </div> | ||
=== '''Structure>RCI''' === | === '''Structure>RCI''' === | ||
& | ''':''': Calculate Random Coil Index | ||
<div> </div><div> </div><div>Random Coil Index is calculated using rci_v_1c.py script. New directory | <div> </div><div> </div><div>Random Coil Index is calculated using rci_v_1c.py script. New directory "rci", that contains results of the calculations, is created </div><div>inside project root directory.</div><div> </div><div> </div> | ||
== View menu == | == View menu == | ||
<div>This section provides means to visualize data and calculation results.</div><div> </div> | <div>This section provides means to visualize data and calculation results.</div><div> </div> | ||
=== '''View>Fragment''' === | === '''View>Fragment''' === | ||
<div> ''':''' | <div> ''':''': Display current properties of selected fragment<br></div><div> </div><div> This command opens “Fragment Graph” (FG) window were all properties of a selected fragment that currently are loaded in memory will be displayed (see Figure). </div><div> [[Image:Fmcgui 2.21.jpg|thumb|left|550px]]</div><div>The chemical shifts making up the fragment are shown in the middle part of the window. The heading line (for example, line “Fragment #21 rem_run1 L108” shown on Figure 2.21) shows fragment ID (#21), the name of directory that contains assignment probabilities used to assign the fragment (rem_run1) and sequence position to which the fragment is assigned (L108). All other fragment properties are shown in the form of graphs. Typing probabilities the top section of the window one can see typing probabilities <span> and both assignment probabilities and are shown on three graph on the top part of the FG window. It is also indicated the name of the directories from which the displayed assignment probabilities were loaded (for example, on the Figure, the P<sub>SA</sub> and P<sub>REM</sub> assignmet probabilities were loaded from directories "sa_run2" and "rem_run2", respectively). Two graphs on the bottom part of the FG window shows a column ( C<sup>f</sup><sub>NOE_B</sub>(''U_id'') or C<sup>f</sup><sub>NOE_F</sub>(''U_id'') ) and a row (C<sup>''U_id''</sup><sub>NOE_B</sub>(f) or </span><span>C<sup>''U_id''</sup><sub>NOE_F</sub>(f) </span><span> ) of contact map calculated from NOESY data, respectively. Here ''U_id'' is selected fragment ID, while ''f'' is any fragment ID for which the corresponding score is > 0. What contact map, C<sub>NOE_B</sub> or C<sub>NOE_F</sub>, is shown on the graphs is indicated on the bar at the bottom of FG window. Clicking on this bar using mouse will switch from one contact map to another.</span></div><div>The contact map <span>C<sub>HNCA</sub> is shown on this graph implicitly as well by a color of graph bars.</span></div><div> For example, on the Figure the elements of contact map <span>C<sub>NOE_F</sub></span> <span>related to the fragment 21 are shown. Scores for the fragment 21 to be before the fragments 37, 20 and 17 in the protein sequence are shown in magenta, which means that these connections being derived from NOESY data are supported by HNCA spectra as well. At the same time, scores for the fragment 21 to be before the fragments 22 and 4 are shown in red meaning that theses connections are not supported by HNCA spectra.</span></div><div> <br></div><div> </div><div> <br></div> | ||
=== '''View>Assignment''' === | === '''View>Assignment''' === | ||
''':''': Display and analize assignment results | |||
<div> </div><div>The user is asked to select a directory were assignment probabilities were calculated | <div> </div><div>The user is asked to select a directory were assignment probabilities were calculated (sa_run# or rem_run#). </div><div> New “Assignments Graphs” (AG) window pops up (see Figure 1, left). This window provides user with graphical means to manipulate PB fragment’s assignments without making changes of assignment status of PB fragments in memory. The current fragment’s assignment, taken from the memory, altogether with different scores associated with this assignment is shown when AG window is opened. The fragments that currently are not assigned are also shown in the bottom-right part of the AG window. User can modify the current assignment and to observe the resulting changes in the scores.</div><div>[[Image:Fmcgui fig2.22.jpg|thumb|left|445px]][[Image:Fmcgui 2.24a.jpg|thumb|right|445px]]</div><div></div><div></div><div> <br></div> | ||
<br> | <br> | ||
<div>Once a new assignment for the sequence segment is selected the current fragments assignment shown on the left part of AG window can be modified by pressing button “Update” (see Figure | <div><br> </div><div> A particular fragment’s assignment is displayed on the left part of the AG window. The graph at the bottom shows ID of fragments assigned to each sequence position, while four graphs on the top show different scores that correspond to this assignment, namely, HNCA score, NOE scores, typing and assignment probabilities that currently loaded in memory. The scroll bar at the bottom allows user to move along protein sequence. <br></div><div> <br></div><div>In order to modify the assignment user have first to select protein sequence segment of interest by clicking on IDs of two residues that correspond to the edges of the segment. The colour of these residues, for example K56 and D59, will change to red (see Figure 1,right). Then pressing on “Select Segment” bar at the bottom of AG window will result in possible assignment for the selected window to be shown on top-right corner of AG window. <div>[[Image:Fmcgui fig2.24.jpg|thumb|right|250px]]</div></div><div>The possible assignments for the selected sequence segment are taken form the selected directory shown on the title bar and correspond to sub-optimal fragment assignments sampled during SA or REM calculations. The new assignment for the sequence segment should be selected from the list of possible assignments using mouse. The color of the selected line will turn red (see Figure 2.23). In the case user want to consider an assignment that is not present in the list, he should select any assignment from the list ant modify it by pressing the button “Edit” at the bottom of the AG window. A new window “Edit Assignment” pops up, and user can modify the assignment by typing changes in this window see Figure 2). </div><div> <br></div><div><br> </div> <div>Once a new assignment for the sequence segment is selected the current fragments assignment shown on the left part of AG window can be modified by pressing button “Update” (see Figure 3)</div><div>[[Image:Fmcgui_fig2.25.jpg|thumb|center|500px]]</div><div> </div> | ||
Latest revision as of 23:33, 5 January 2010
Project Menu
Project>New
:: Start a new project
User have to provide a name of the project PROJECTNAME,and to select a directory that will host the project. The project root directory with the same name PROJECTNAME is created.
Project>Load
:: Continue to work with previously saved project.
User have to select file PROJECTNAME.prj in the directory PROJECTNAME, where PROJECTNAME is the name of the root directory of the project.
Project>Save
:: Save the current state of the project
What is currently in the computer memory is saved in the file PROJECTNAME.prj located in the root directory of the project.
Poject>Quit
:: Save the current state of the project and quit.
Data Menu
The DATA section serves to load&save the data (such as protein sequence and peak lists). Since there are different formats of data-files that could be loaded in memory or saved on disk, one can use this section as format converter as well.
Data>Protein Sequence>Load
:: Load protein sequence into memory.
The input formats:
1-letter code (fasta format);
3-letter code (standard format).
User have to select the file with sequence and to specify the first residue ID in case it is not specified in the input file. If there is His-tag in the sequence file, it is recommended to set the first residue ID to a negative number such that the first residue of a protein has ID of 1.
Data>Protein Sequence>Save as
:: Save protein sequence in file on disk.
The output formats:
1-letter code (fasta format);
3-letter code (standard format, for cyana)
3-letter code (for AutoStructure);
3-letter code (for RCI);
Data>N15 NOESY>
:: Load or save peak list.
Referencing is not required for NOESY peak list
Input and output formats:
Data>C13 NOESY>
:: Load or save peak list.
Referencing is not required for NOESY peak list
Input and output formats:
Data>Arom NOESY>
:: Load or save peak list.
Referencing is not required for NOESY peak list
Input and output formats:
Data>N15 HSQC>
:: Load or save peak list.
This peak list should be referenced.
Input and output formats:
Data>C13 HSQC>
Data>HNCA>
:: Load or save peak list.
It is recommended this peak list to be referenced.
Input and output formats:
Data>HNCO>
:: Load or save peak list.
It is not required this peak list to be referenced.
Input and output formats:
Data>CBCACONHN>
:: Load or save peak list.
Following these warnings user can check and modify generated PB fragments in the left section of “Create Fragment’ window.
Then the graph will show typing probabilities Tt1(f) <span />that correspond to the selected residue type t1 for all available fragments IDs f .
Results:
Results:
It is not reauired this peak list to be referenced.
Input and output formats:
Data>HBHACONH>
Data>Tolerances
:: Set tolerances for chemical shift matching in different spectral dimensions.
There are six tolerances to set up:
NX - tolerance for matching resonances in N15 dimension
CX - tolerance for matching resonances in C13 dimension
HN - tolerance for matching resonances in HN direct dimension
HC - tolerance for matching resonances in HC direct dimension
Hi - tolerance for matching resonances in H indirect dimension
Ci - tolerance for matching resonances in C13 indirect dimension
Fragment>Load>assigned
:: Load assigned chemical shifts (spin-systems) in memory.
Prerequisites:
♦ Loaded sequence
Input formats:
Fragment>Load>PB fragments
:: Load unassigned spin-systems (PB fragments) in memory.
Input format :
Fragment>Save>PB fragments
:: Save PB-fragments in a file on disk.
All PB fragments are saved in standard format.
The name of the saved file and it’s location are specified by user.
There are 3 options to save PB-fragments in the file:
- In the order of fragments index, that is in the order by which fragments are stored in memory;
- In the order of fragments user ID
- In the order of fragments assignment ID, A_id. In this case two files are saved. One file, with user specified name 'user_name', contains only fragments assigned to protein sequence positions, that is to positions with residue ID of >= 1. The second file, with the name 'user_name_na', contains all not assigned fragments (that is fragments with A_id = -99).
Fragment>Save>cyana
:: Save assigned chemical shifts (that is fragments with assignment ID > 0 ) in CYANA format.
Fragment>Save>bmrb
:: Save assigned chemical shifts (that is fragments with assignment ID > 0 ) in the format suitable for BMRB deposition (star2.1)
Fragment>Save>talos
:: Save assigned chemical shifts (that is fragments with assignment ID > 0 ) in the format suitable for TALOS/CS-Rosetta;
Fragment>Save>abacus
:: Save assigned chemical shifts (that is fragments with assignment ID > 0 ) in the format suitable for BACUS;
Fragment>Create>fawn
:: Create and evaluate bPB fragments
Prerequisites:
♦ loaded in memory referenced peak lists of CBCA(CO)HN, HBHA(CO)HN, N15HSQC, and HNCA spectra;
♦ Specified tolerances.
There are two steps in executing this command.
- On the first step, a fake C13HSQC peak list is created and shown in the popped up window “fake C13HSQC” .
User can use the information shown in the main FMCGUI window and check /edit the list in the entry section of “fake C13HSQC” window accordingly.
Pressing OK will result in loading the peak list from the entry window into memory as C13HSQC peak list.
2. On the second step, a number of bPB fragments corresponding to 20 different amino acid types are generated from user-identified spin-systems.
Each generated bPB fragment is evaluated by a score S(T) that measure how good the spin-system chemical shifts match corresponding statistical chemical shifts derived from BMRB database (see Figure 1.2). The bPB fragment with highest score is selected to form a list of bPB-fragments.
The warning messages of the command are shown in the main FMCGUI window (see for example, Figure 2.4).
In addition to that, a new window ‘Create Fragment’ pops up:
.
"Create Fragment" window consists of three sections. The left sections contains suggested bPB fragments, while the other sections contains two reports of fragments scoring with both C and H resonances and with only C resonances, respectively. Following the warning messages shown in the main FMCGUI window (see Figure 2.4), user can accept/modify generated bPB fragments. Alternatively, when ‘poor’ bPB fragments are present, user can go back to spectra, fix the pick lists accordingly, and repeat the fragment generation again.
User-approved bPB fragments will be loaded in the memory by pressing OK button.
Results:
♦ C13HSQC peak list loaded in memory
♦ bPB fragments are loaded in memory
Fragment>Create>abacus
:: Create and avaluate PB fragments.
Prerequisites:
♦ Loaded in memory referenced C13HSQC, N15HSQC, and HNCA peak lists and not referenced CBCA(CO)HN peak list.
( As an option, HNCA peak list could be not referenced as well)
♦ Specified tolerances.
Results:
♦ PB fragments are loaded in memory
A number of PB fragments corresponding to 20 different amino acid types are generated from user-identified spin-systems. Each generated PB fragment is evaluated by a score S (see Figure 1.2) that measure how good the spin-system chemical shifts match corresponding statistical chemical shifts derived from BMRB database. The PB fragment with highest score is selected to form a list of PB fragments.
Spin-system which have Smax less than 10-4 are reported in the warning lines shown in the main FMCGUI window:
Figure 2.4
Following these warnings user can check and modify generated PB fragments in the left section of “Create Fragment’ window.
.
"Create Fragment" window consists of three sections. The left sections contains suggested PB fragments, while the other two sections contain report of fragment's scoring with both C and H atom resonances and with only C atom resonances, respectively. Following the warning messages shown in the main FMCGUI window (see Figure 2.4), user can accept/modify generated bPB fragments.
Pressing OK button, user-approved PB fragments shown in the left section of "Create Fragment" window will be loaded in the memory.
Alternatively, when ‘poor’ PB fragments are present, user can go back to spectra, fix the pick lists accordingly, and repeat the fragment generation again.
Fragment>Type>Calculate>
:: Probabilistic typing of bPB (fawn) or PB (abacus) fragments
Prerequisites:
♦ loaded in memory protein sequence
♦ loaded in memory PB-fragments
♦ specified tolerances.
Results:
♦ Fragment typing probabilities are calculated and loaded in memory.
The main FMCGUI window displays (see Figure ):
- the summary table that shows how many fragments of each AA-residue type are expected and how many fragments were actually recognized by the typing script
- warning messages that suggest user to check and possibly modify typing manually of some fragments manually (see command {Fragment>Type>fix})
Fragment>Type>fix
:: Modify typing probabilities manually.
Prerequisites:
♦ loaded in memory protein sequence
♦ loaded in memory PB fragments
The command opens "Fragment Property Modification" (FPM) window (see Figure).
This window has 3 sections.
In the top section of FPM window user can select fragment user ID, U_id.
Then typing probabilities Tt(U_id) for all 20 amino acid types t will be shown on the graph. The chemical shifts of the fragments and its assignment status (A_id) are shown as well. User can modify typing probabilities of the fragment U_id . Selecting amino acid types on the graph using right mouse button and pressing ‘Update’ button on the bottom of FPM window will set probabilities Tt(U_id) corresponding to the selected amino acid types to the same non-zero value (such that the normalization condition is satisfied).
In the middle section of FPM window user can select amino acid residue type t1 ( for example, on the following Figure, t1 = P ).
Then the graph will show typing probabilities Tt1(f) <span />that correspond to the selected residue type t1 for all available fragments IDs f .
Selecting a particular fragment U_id on the graph using right mouse button (the color of U_id will be changed to red) and then pressing “Update” button will set the probability Tt1(U_id ) to 1.
Selecting a particular fragment U_id on the graph using left mouse button (the color of U_id wil be changed to blue) and then pressing “Update” button will set the probability Tt1(U_id ) to 0.
Fragment>Expected Peaks>
:: Generate lists of peaks expected from covalent structure of PB fragments.
Prerequisites:
♦ protein sequence loaded in memory
♦ PB-fragments loaded in memory
Expected peak lists of the following spectra could be generated: N15-NOESY, C13-NOESY, H(C)CH-TOCSY, (H)CCH-TOCSY,
N15-HSQC, and C13-HSQC. Generated peak list is saved to the file on disk in SPARKY format.
Fragment>Modify assigned
:: Correct assigned fragments.
Prerequisites:
♦ loaded in memory protein sequence
♦ loaded in memory PB-fragments
♦ loaded HNCO, CBCACONH, and HNCA peak lists.
♦ specified tolerances.
Resuls:
♦ CO chemical shifts are assigned
♦ chemical shift names are corrected for assigned PB-fragments (the fragments with A_id > -99 )
Assignment>Contacts>HNCA
:: Calculate fragments contact map CHNCA
Prerequisites:
♦ PB-fragments are loaded in memory
♦ typing probabilities are calculated
♦ HNCA peak list (recommended to be referenced) is loaded in memory
♦ tolerances are specified
Results:
♦ Contact map CHNCA is calculated and loaded in memory
During the execution of this command the loaded HNCA peak list is validated. The validation results are shown in the project main window.
User have to check suspicious peaks going back to the spectra and correct the input peak list if it is needed.
Assignment>Contacts>NOE>fawn
:: Calculate fragment contact map CNOE_F.
Prerequisites:
♦ PB-fragments are loaded in memory;
♦ typing probabilities are calculated;
♦ N15-NOESY peak list is loaded in memory.
♦ tolerances are specified.
Results:
♦ contact map CNOE_F is calculated using N15-NOESY peak list and loaded in memory
Assignment>Contacts>NOE>abacus
:: Calculate both CNOE_B and CNOE_F contact maps.
Prerequisites:
♦ PB-fragments are loaded in memory
♦ typing probabilities are calculated
♦ N15-NOESY and C13-NOESY peak lists are loaded in memory
♦ tolerances are specified
Results:
♦ CNOE_F contact map is calculated using only N15_NOESY peak list
♦ CNOE_B contact map is calculated using both N15_NOESY and C13_NOESY peak lists that are interpreted by BACUS procedure
♦ Both calculated maps are loaded in memory
Assignment>Calculate Probabilities>SA
:: Calculate assignment probabilities
Prerequisites:
♦ protein sequence is loaded in memory
♦ PB-fragments are loaded in memory
♦ typing probabilities are calculated
♦ fragments contact map CHNCA is calculated
♦ at least one contact map CNOE_F or CNOE_B is calculated
Results:
♦ Assignment probability map PSA is calculated and loaded in memory.
Probabilistic mapping of PB fragments onto protein sequence is performed using Simulated Annealing (SA) Monte Carlo simulations.
A new window “Calculate SA” is open were user can specify different parameters in the control file of the SA simulations (see Figure 2.9).
The main parameters to consider are:
<span />
- “Name of the SA run”. Normally the name is sa_run#. A new directory under this name will be created within PROJECTNAME/assign directory. SA calculations will be curried out and the results will be stored in this directory.
- “Size of the pool for unassigned fragments”. The number of positions that are appended to the protein sequence and discarded (unassigned) fragments, if there are any, will be located there. It is safe to over-estimate this number. (If this number is under-estimated, this will force the mapping of spurious spin-systems onto protein sequence);
- “Number of SA trajectories”. The time needed for calculations is proportional to this number. On the other hand, having more SA trajectories the assignment probabilities could be calculated more accurately. In the case of good data, when all SA trajectories converge to assignments with the same energy, 10-15 trajectories should be enough. In the case of poor data, it is better to calculate 40-50 SA trajectories.
- “NOE bbcmap type”. User should specify which one NOE contact map, (abacus) or (fawn) should be used in the calculations;
- “Fixing position flag”. If the flag is set to 1, sequence position of all fragments which has assignment ID > -99 will be fixed during the simulation.
- “Final temperature”. Setting the optimal final temperature will provide all SA trajectories converge to optimal or sub-optimal assignments (the assignments that are in the vicinity of the global energy minimum). The optimal final temperature could be find by running one or a few SA runs with 3-4 trajectories and by analysing convergence of the trajectories from the report shown in the main FMCGUI window after each run (see Figure 2.10).
In the result of the SA calculations assignment probability map PSA is calculated and loaded in memory. The map is also saved in the file 'sa.probmap' located into PROJECTNAME/assignment/sa_run# directory.
Assignment>Calculate Probabilities>REM
:: Calculate assignment probabilities.
Prerequisites:
♦ protein sequence is loaded in memory
♦ PB fragments are loaded in memory
♦ typing probabilities are calculated
♦ fragments contact map CHNCA is calculated
♦ at least one contact map CNOE_F or CNOE_B is calculated
Results:
♦ assignment probbilities PREM are calculated and loaded in memory
Probabilistic mapping of PB fragments onto protein sequence is performed using Replica Exchange Method Monte Carlo simulations.
A new window “Calculate REM” is open were user can specify different parameters in the control file of the REM simulations (see Figure 2.11)
The main parameters to consider are:
- “Name of the REM run”. Normally the name is rem_run#. A new directory under this name will be created within PROJECTNAME/assign directory. REM calculations will be curried out and the results will be stored in this directory.
- “Size of the pool for unassigned fragments”. The number of positions that are appended to the protein sequence and discarded (unassigned) fragments, if there are any, will be located there. It is safe to over-estimate this number. (If this number is under-estimated, this will force the mapping of spurious spin-systems onto protein sequence);
- “Number of REM steps”. The time needed for calculations is proportional to this number. On the other hand, with more REM steps more extensive sampling of assignment space wil be achieved, which in turn results in more accurate assignment probabilities. This number should be increased for large proteins.
- “NOE bbcmap type”. User should specify which one NOE contact map, (abacus) or (fawn) should be used in the calculations;
- “Fixing position flag”. If the flag is set to 1, sequence position of all fragments which has assignment ID > -99 will be fixed during the simulation.
- “Low temperature”. The optimal low temperature will provide extensive sampling of should sub-optimal assignments during REM simulation.
User can check a correct setting of the low temperature as well as the number of REM steps by analysing a report shown in the main FMCGUI window after the calculations are done.
In the result, the 50 lowest energy assignments are used to calculate assignment probabilities . The probabilities PREM are loaded in memory and saved in the file ‘rem.prbmap’ located into PROJECTNAME/assign/rem_run# directory.
Assignment>Fix Assignment>Using Probability map
:: Perform sequence specific assignment of PB fragments using results of SA or REM calculations.
Prerequisites:
♦ protein sequence is loaded in memory;
♦ PB fragments loaded in memory;
♦ at least one SA or REM calculation of assignment probabilities is done
Results:
♦ Assignment ID for part (or all) of PB -fragments is specified
Calculation of assignment probabilities with FMCGUI could be repeated a few times using different methods and parameters. Results of each calculation are stored in a separate directory with the user-specified name. Therefore, there could be a few different directories (for example, sa_run1, sa_run2, rem_run0, rem_run1, rem_run2) located within PROJECTNAME/assign/ directory that contain different assignment probability maps.
User will be asked to select the calculation directory (sa_run# or rem_run#) and to specify the value of acceptance probability Pmin. Normally, Pmin =0.9 is appropriate. A fragment is considered to be assigned to a sequence position if the corresponding assignment probability (taken from the selected directory) is >= Pmin.
In the result, part of the PB-fragments will be assigned, namely, their assignment IDs A_id will be specified. The assignment report will be shown in the main FMCGUI window (see Figure 2.13) and saved in the corresponding simulation directory (‘sa.fix ‘or ‘rema.fix’ files) as well.
For each sequence position, the IDs of both unambiguously and ambiguously assigned fragments are shown in the report. The list of discarded fragments is also presented (see Figure 2.13).
Assignment>Fix Assignment>Manually
:: Fix sequence position of PB fragments manually
Prerequisites:
♦- protein sequence is loaded in memory;
♦- PB-fragments are loaded in memory;
Results:
♦ Assignment ID for selectedf PB fragments is specified
♦ Assignment ID for selectedf PB fragments is specified
“Fragment Property Modification” (FPM) window will be open (see Figure 2.14). User can select a fragment in the bottom section of the window and the ingormation regarding the fragment assignemts will be displayed in this section. Namely, the graph shows assignment probabilities that are currently loaded in memory and the text part shows the chemical shifts making up the fragment and it’s assignment ID (A_id). Which assignment probabilities, PSA or PREM, are shown on the graph can be selectes by pressing left button at the bottom of FPM window.
To modify the current fragment assignment user have to select sequence position on the graph (using mouse right button) and then to press "Update" button. In the result, Assignment ID of the selected fragment will be set to the selected sequence position.
There are two special positions “U” and “B” shown on the graph at the end of the protein sequence (see Figure 2.15). Selecting “U” and pressing “Update” button results in changing assignment status of the fragments to Unassigned, that is A_id is set to -99. Selecting “B” and pressing “Update” button results in fixing fragment position in the pool of discarded fragments.
Assignment>Fix Assignment>Reset all
:: Change assignment status of all PB fragments to “Unassigned”.
Prerequisites:
♦ protein sequence is loaded in memory
♦ PB fragments are loaded in memory
Results:
♦ Assignment ID for all PB fragments is set to -99
Assignment>Load Probabilities
:: Load assignment probabilities
Prerequisites:
♦ protein sequence is loaded in memory;
♦ PB fragments are loaded in memory;
♦ SA or REM calculations of assignment probabilities are done
User have to select a directory where SA or REM calculations were done (sa_run# or rem_run#).
The assignment probabilities from the selected directory will be loaded in memory.
Structure>Constraints>Talos>calculate
:: Generate dihedral angle constraints.
Prerequisites:
♦ protein sequence is loaded in memory
♦ assigned shemical shifts are loaded in memory
Backbone dihedral angles are predicted using TALOS and then transformed to dihedral angle contraints.
In the result, the constraints are saved in two formats: files 'belok.aco' and 'prot_dihe.tbl' for CYANA and CNS calculations, respectively. Both files are saved in the project root directory.
Structure>Constraints>Talos>load
:: Load dihedral angle constraints in CYANA format (aco-file).
Structure>Constraints>H-bonds>Specify
:: Specify H-bond distance constraints manually.
Prerequisites:
♦ protein sequence is loaded in memory
For each H-bond constraint, user have to specify O-HN pair by typing in the ID of residues corresponding to O and HN atoms. Pressing “OK” button will save H-Bond constraints on disk in two formats: files 'hbond.upl' and 'prot_hbond.tbl' for CYANA and CNS calculations, respectively. Both files are saved in the project root project directory.
Structure>Constraints>H-Bonds>load
:: Load HB-ond distance constraints in CYANA format (upl-file)
Structure>Calculate>Cyana
:: Set up new structure calculation run with CAYANA.
Prerequisites:
♦ protein sequence is loaded in memory
♦ assigned chemical shifts are loaded in memory
♦ N15-NOESY, C13-NOESY, and Arom-NOESY peak lists are loaded in memory
♦ specified tolerances
♦ dihedral angle constraints are created ( i.e. file “belok.aco” is present in the project root directory)
♦ Optional:
♦ HBond distance constraints (file “hbond.upl” is present in the project root directory)
♦ file that contains ZN ion ligands (file “zn_ligands” is present in the project root directory)
A new directory under the name that is specified by user (normally “crun#”) will be created inside project root directory PROJECTNAME. This directory contains all files required to start automatic structure calculations with CYANA.
Structure>Calculate>ABCUS
:: Sset up ABACUS structure calculations
Structure>RPF>RP
:: Perform Recall/Precision analysis of structural ensemble.
Prerequisites:
♦ protein sequence is loaded in memory
♦ assigned shemical shifts are loaded in memory
♦ N15-NOESY, C13-NOESY, and Arom-NOESY peak lists are loaded in memory
♦ specified tolerances
♦ coordinates of structural ensemble in CYANA format (final.pdb) or in CNS format
Results:
♦ summary report
♦ flase negative and false positive peak lists
The parameters to set up (see Figure ):
- “RP directory name”. Normally the name is rp#. A new directory under this name will be created within PROJECTNAME/assign directory. The results of the RPFanalysis will be stored in this directory.
- “sequence gap”. Residue pairs separated by less than the value of sequence gap will be excluded from generating expected peak lists.
- “cutting distance for recall”. Distance threshold for evaluating matching of an experimental peak to a structural ensemble (Recall score)
- “cutting distance for precition”. Distance threshold for generating expected peak from structural ensemble (Precision score)
The results of the RPF analysis will be saved in a new directory (for example, directory "rp1" on Figure) that is located in the project root directory. Results include peak lists in the SPARKY format of both false negative and false positive peaks for different NOESY spectra in a separate files.
Structure>RPF>DP
- : Set up calculations of DP score with AutoStructure
Prerequisites:
♦ protein sequence is loaded in memory
♦ assigned shemical shifts are loaded in memory
♦ N15-NOESY, C13-NOESY, and Arom-NOESY peak lists are loaded in memory
♦ specified tolerances
♦ coordinates of structural ensemble in CYANA format (final.pdb) or in CNS format
Structure>Water refinement>calculate
:: Set up water refinement calculations
Prerequisites:
♦ structure calculation with CYANA should be done
♦ Optional:
♦ RDC data in PALES format.
♦ file with ZN ion ligands (file “zn_ligands” )
User will be asked to :
- specify the name of directory for water refinement calculations, WATDIR;
- select a number of files with coordinates and constraints;
- indicate cisProline residues (if there are any)
- specify protonation state of HIS residues (which is double protonated by default).
In the result, a new directory WATDIR will be created inside project root directory. I will contain all files and scripts required to carry out water refinement calculations with CNS.
The calculations should be started outside FMCGUI. (linux cluster is is recommended)
Structure>Water refinement>summary
:: Create a summary for water refined structural ensemble.
Prerequisites:
♦ water refinement with CNS should be done.
User will be asked to specify residues used for structure superposition and to select the water refinement directory, WATDIR (see Figure).
The refined structural models will be superimpose and combined in one file. The created summary reports values of different energy components and constraint violation statistic for each structural model. A new directory WATDIR_results will be created. The directory contains final superimposed coordinates, distance and dihedral angle constrains in a format suitable for PDB deposition, summary and constraint violations report.
Structure>Add ZN ligands
:: Create "zn_ligands" file.
“ZN ligands” window pops up. User have to type in IDs of residues that ligate zinc ion. If there are a few zinc ions, information for each ion should be provided in a separate row (see Figure).
The file “zn_ligands” will be created inside the project root directory by pressing “OK” button.
Structure>RCI
:: Calculate Random Coil Index
Random Coil Index is calculated using rci_v_1c.py script. New directory "rci", that contains results of the calculations, is created
inside project root directory.
This section provides means to visualize data and calculation results.
View>Fragment
:: Display current properties of selected fragment
This command opens “Fragment Graph” (FG) window were all properties of a selected fragment that currently are loaded in memory will be displayed (see Figure).
The chemical shifts making up the fragment are shown in the middle part of the window. The heading line (for example, line “Fragment #21 rem_run1 L108” shown on Figure 2.21) shows fragment ID (#21), the name of directory that contains assignment probabilities used to assign the fragment (rem_run1) and sequence position to which the fragment is assigned (L108). All other fragment properties are shown in the form of graphs. Typing probabilities the top section of the window one can see typing probabilities and both assignment probabilities and are shown on three graph on the top part of the FG window. It is also indicated the name of the directories from which the displayed assignment probabilities were loaded (for example, on the Figure, the PSA and PREM assignmet probabilities were loaded from directories "sa_run2" and "rem_run2", respectively). Two graphs on the bottom part of the FG window shows a column ( CfNOE_B(U_id) or CfNOE_F(U_id) ) and a row (CU_idNOE_B(f) or CU_idNOE_F(f) ) of contact map calculated from NOESY data, respectively. Here U_id is selected fragment ID, while f is any fragment ID for which the corresponding score is > 0. What contact map, CNOE_B or CNOE_F, is shown on the graphs is indicated on the bar at the bottom of FG window. Clicking on this bar using mouse will switch from one contact map to another.
The contact map CHNCA is shown on this graph implicitly as well by a color of graph bars.
For example, on the Figure the elements of contact map CNOE_F related to the fragment 21 are shown. Scores for the fragment 21 to be before the fragments 37, 20 and 17 in the protein sequence are shown in magenta, which means that these connections being derived from NOESY data are supported by HNCA spectra as well. At the same time, scores for the fragment 21 to be before the fragments 22 and 4 are shown in red meaning that theses connections are not supported by HNCA spectra.
View>Assignment
:: Display and analize assignment results
The user is asked to select a directory were assignment probabilities were calculated (sa_run# or rem_run#).
New “Assignments Graphs” (AG) window pops up (see Figure 1, left). This window provides user with graphical means to manipulate PB fragment’s assignments without making changes of assignment status of PB fragments in memory. The current fragment’s assignment, taken from the memory, altogether with different scores associated with this assignment is shown when AG window is opened. The fragments that currently are not assigned are also shown in the bottom-right part of the AG window. User can modify the current assignment and to observe the resulting changes in the scores.
A particular fragment’s assignment is displayed on the left part of the AG window. The graph at the bottom shows ID of fragments assigned to each sequence position, while four graphs on the top show different scores that correspond to this assignment, namely, HNCA score, NOE scores, typing and assignment probabilities that currently loaded in memory. The scroll bar at the bottom allows user to move along protein sequence.
In order to modify the assignment user have first to select protein sequence segment of interest by clicking on IDs of two residues that correspond to the edges of the segment. The colour of these residues, for example K56 and D59, will change to red (see Figure 1,right). Then pressing on “Select Segment” bar at the bottom of AG window will result in possible assignment for the selected window to be shown on top-right corner of AG window.
The possible assignments for the selected sequence segment are taken form the selected directory shown on the title bar and correspond to sub-optimal fragment assignments sampled during SA or REM calculations. The new assignment for the sequence segment should be selected from the list of possible assignments using mouse. The color of the selected line will turn red (see Figure 2.23). In the case user want to consider an assignment that is not present in the list, he should select any assignment from the list ant modify it by pressing the button “Edit” at the bottom of the AG window. A new window “Edit Assignment” pops up, and user can modify the assignment by typing changes in this window see Figure 2).
Once a new assignment for the sequence segment is selected the current fragments assignment shown on the left part of AG window can be modified by pressing button “Update” (see Figure 3)
