Validation of Consensus Run: Difference between revisions

From NESG Wiki
Jump to navigation Jump to search
(Created page with '== '''Validation of Consensus Run''' == === '''Verification of Consensus Assignments''' === # Create a working subdirectory, for example, <tt>cyana21/calc2</tt>. # Copy consens…')
 
No edit summary
Line 1: Line 1:
== '''Validation of Consensus Run''' ==
== '''Validation of Consensus Run''' ==


=== '''Verification of Consensus Assignments''' ===
=== '''Verification of Consensus Assignments''' ===


# Create a working subdirectory, for example, <tt>cyana21/calc2</tt>.
#Create a working subdirectory, for example, <tt>cyana21/calc2</tt>.  
# Copy consensus peaklists (<tt>n_cons.peaks</tt>, <tt>ali_cons.peaks</tt> and <tt>aro_cons.peaks</tt>) into this subdirectory.
#Copy consensus peaklists (<tt>n_cons.peaks</tt>, <tt>ali_cons.peaks</tt> and <tt>aro_cons.peaks</tt>) into this subdirectory.
** If using CARA, re-integrate the peakslist and export them.
 
# Copy <tt>init.cya</tt>, <tt>XXXX.seq</tt>, <tt>XXXX.prot</tt> and <tt>stereofound.cya</tt> from the CYANA automated run (<tt>cyana21/calc2</tt>). Rename <tt>stereofound.cya</tt> into <tt>stereo.cya</tt>
*If using CARA, re-integrate the peakslist and export them.
# Copy <tt>GridSearch.cya</tt>, <tt>CALC.cya</tt> and <tt>talos.aco</tt> from the <tt>found</tt> subdirectory. Modify <tt>GridSearch.cya</tt> to read <tt>cons.upl</tt> instead of <tt>short.upl</tt>.
 
# Create UPL file <tt>cons.upl</tt> using run_caliba.cya (see below).
#Copy <tt>init.cya</tt>, <tt>XXXX.seq</tt>, <tt>XXXX.prot</tt> and <tt>stereofound.cya</tt> from the CYANA automated run (<tt>cyana21/calc2</tt>). Rename <tt>stereofound.cya</tt> into <tt>stereo.cya</tt>  
# In CYANA 2.1 run <tt>GridSearch.cya</tt>
#Copy <tt>GridSearch.cya</tt>, <tt>CALC.cya</tt> and <tt>talos.aco</tt> from the <tt>found</tt> subdirectory. Modify <tt>GridSearch.cya</tt> to read <tt>cons.upl</tt> instead of <tt>short.upl</tt>.  
# In CYANA 2.1 run CALC.cya for manual structure calculation with only consensus constraints.
#Create UPL file <tt>cons.upl</tt> using run_caliba.cya (see below).  
# Clean violations the same was is in [[NESG:UsingFOUNDakaHABAS|Verifying Local NOE Network with FOUND (HABAS)]] and try to keep TF < 1.0.
#In CYANA 2.1 run <tt>GridSearch.cya</tt>  
** For the violated NOE, re-integrate the peaks first if it is overintegrated (in XEASY).
#In CYANA 2.1 run CALC.cya for manual structure calculation with only consensus constraints.  
** Unassign peaks, which appear at spectral artifacts.
#Clean violations the same was is in [[NESG:UsingFOUNDakaHABAS|Verifying Local NOE Network with FOUND (HABAS)]] and try to keep TF &lt; 1.0.
** If violation persisit while peak integrals are correct, then modify the calibration parameters.
*For the violated NOE, re-integrate the peaks first if it is overintegrated (in XEASY).  
# Verify the additional stereospecific assignments made in this run. Swap the atom groups if needed.
*Unassign peaks, which appear at spectral artifacts.  
*If violation persisit while peak integrals are correct, then modify the calibration parameters.
 
#Verify the additional stereospecific assignments made in this run. Swap the atom groups if needed.


You may want to run an additional iteration of consensus assignment, though it is usually not necessary:  
You may want to run an additional iteration of consensus assignment, though it is usually not necessary:  


=== '''Acceptance Criteria''' ===
=== '''Acceptance Criteria''' ===


A successful consensus run shall satisfy:
A successful consensus run shall satisfy:  
# Both program CYANA and AutoStructure automatic structure calculation deliver converged structure.
# The structures from program CYANA and AutoStructure are similar, the RMSD value between the mean structures from both programs should be less than the sum of pairwise RMSD values of the two structure families.
# The manual run of CYANA structure calculation with distance constraints from consensus assigned NOEs delivers converged structure, the pairwise RMSD of the structure family shall be less than 3A.
# The yield of consensus assigned NOEs shall be large than 50%.


=== '''Example scripts''' ===
#Both program CYANA and AutoStructure automatic structure calculation deliver converged structure.
#The structures from program CYANA and AutoStructure are similar, the RMSD value between the mean structures from both programs should be less than the sum of pairwise RMSD values of the two structure families.
#The manual run of CYANA structure calculation with distance constraints from consensus assigned NOEs delivers converged structure, the pairwise RMSD of the structure family shall be less than 3A.
#The yield of consensus assigned NOEs shall be large than 50%.


==== '''run_caliba.cya''' ====
=== '''Example scripts'''  ===
 
==== '''run_caliba.cya''' ====


<nowiki>
<nowiki>
Line 52: Line 56:
write upl cons.upl
write upl cons.upl
protocol:=
protocol:=
</nowiki>
</nowiki>  
 
%COMMENT%


%COMMENT%


-- Main.GaohuaLiu - 26 Jan 2007
<br> -- Main.GaohuaLiu - 26 Jan 2007

Revision as of 20:18, 4 November 2009

Validation of Consensus Run

Verification of Consensus Assignments

  1. Create a working subdirectory, for example, cyana21/calc2.
  2. Copy consensus peaklists (n_cons.peaks, ali_cons.peaks and aro_cons.peaks) into this subdirectory.
  • If using CARA, re-integrate the peakslist and export them.
  1. Copy init.cya, XXXX.seq, XXXX.prot and stereofound.cya from the CYANA automated run (cyana21/calc2). Rename stereofound.cya into stereo.cya
  2. Copy GridSearch.cya, CALC.cya and talos.aco from the found subdirectory. Modify GridSearch.cya to read cons.upl instead of short.upl.
  3. Create UPL file cons.upl using run_caliba.cya (see below).
  4. In CYANA 2.1 run GridSearch.cya
  5. In CYANA 2.1 run CALC.cya for manual structure calculation with only consensus constraints.
  6. Clean violations the same was is in Verifying Local NOE Network with FOUND (HABAS) and try to keep TF < 1.0.
  • For the violated NOE, re-integrate the peaks first if it is overintegrated (in XEASY).
  • Unassign peaks, which appear at spectral artifacts.
  • If violation persisit while peak integrals are correct, then modify the calibration parameters.
  1. Verify the additional stereospecific assignments made in this run. Swap the atom groups if needed.

You may want to run an additional iteration of consensus assignment, though it is usually not necessary:

Acceptance Criteria

A successful consensus run shall satisfy:

  1. Both program CYANA and AutoStructure automatic structure calculation deliver converged structure.
  2. The structures from program CYANA and AutoStructure are similar, the RMSD value between the mean structures from both programs should be less than the sum of pairwise RMSD values of the two structure families.
  3. The manual run of CYANA structure calculation with distance constraints from consensus assigned NOEs delivers converged structure, the pairwise RMSD of the structure family shall be less than 3A.
  4. The yield of consensus assigned NOEs shall be large than 50%.

Example scripts

run_caliba.cya

protocol:=caliba.log read prot $name read peaks n_cons peaks select "**" peaks set volume=abs(volume) caliba bb=1.7E+06 read peaks ali_cons peaks select "**" peaks set volume=abs(volume) caliba bb=1.7E+06 read peaks aro_cons peaks select "**" peaks set volume=abs(volume) caliba bb=1.7E+06 write upl cons.upl protocol:=

%COMMENT%


-- Main.GaohuaLiu - 26 Jan 2007