NOE Calibration Using CYANA: Difference between revisions

From NESG Wiki
Jump to navigation Jump to search
(Created page with '== '''Introduction''' == NOE calibration converts NOESY cross-peak intensities into upper distance limits. A certain functional relationship is assumed between peak intensities …')
 
No edit summary
Line 1: Line 1:
== '''Introduction''' ==
== '''Introduction''' ==


NOE calibration converts NOESY cross-peak intensities into upper distance limits. A certain functional relationship is assumed between peak intensities and UPLs for a given group of peaks. The most common is NOE peak volume proportional to interproton distance to the minus six: <br>  
NOE calibration converts NOESY cross-peak intensities into upper distance limits. A certain functional relationship is assumed between peak intensities and UPLs for a given group of peaks. The most common is NOE peak volume proportional to interproton distance to the minus six: <br>  
<pre> V=A/d^6</pre>
<pre> V=A/d^6</pre>  
where <code>V</code> is the peak volume (or intensity), <code>d</code> is the upper distance limit, and <code>A</code> is the calibration constant.
where <code>V</code> is the peak volume (or intensity), <code>d</code> is the upper distance limit, and <code>A</code> is the calibration constant.  


<br>


== '''NOE Calibration Using CYANA 2.1'''  ==


== '''NOE Calibration Using CYANA 2.1''' ==
This section describes how to determine NOE calibration constants using CYANA 2.1.


This section describes how to determine NOE calibration constants using CYANA 2.1.
Before running the calibration, the N chemical shifts in the peak list must be updated, and the peak lists separated into a list of N noes and C noes, using the UBNMR macro, '''macro1'''.  


Before running the calibration, the N chemical shifts in the peak list must be updated, and the peak lists separated into a list of N noes and C noes, using the UBNMR macro, '''macro1'''.
=== '''Command "calibrate"''' ===
 
=== '''Command "calibrate"''' ===


A CYANA command <code>calibrate</code> is the basic calibration tool, unchanged since DYANA.  
A CYANA command <code>calibrate</code> is the basic calibration tool, unchanged since DYANA.  


 
<br>
<pre>    Parameters: f(d)                      (required)
<pre>    Parameters: f(d)                      (required)
                 dmin [d2 ...] dmax        (default: 2.4 5.5)
                 dmin [d2 ...] dmax        (default: 2.4 5.5)
Line 32: Line 32:
Type <code>help calibrate</code> at CYANA prompt for more information or see the DYANA documentation.  
Type <code>help calibrate</code> at CYANA prompt for more information or see the DYANA documentation.  


<br>


 
=== '''Macro "caliba"''' ===
=== ''''''Macro "caliba"'''''' ===


CYANA macro <code>caliba.cya</code> is functionally the same as <code>caliba.dya</code> of DYANA. It is a more sophisticated calibration tool aimed at getting more nuanced calibration than a uniform <code>1/d</code><sup><code>6</code></sup> dependence. CALIBA is recommended for manual structure calculation with CYANA 2.1. The <code>caliba</code> macro in fact uses <code>calibrate</code>, therefore, it also uses "center" averaging.  
CYANA macro <code>caliba.cya</code> is functionally the same as <code>caliba.dya</code> of DYANA. It is a more sophisticated calibration tool aimed at getting more nuanced calibration than a uniform <code>1/d</code><sup><code>6</code></sup> dependence. CALIBA is recommended for manual structure calculation with CYANA 2.1. The <code>caliba</code> macro in fact uses <code>calibrate</code>, therefore, it also uses "center" averaging.  


 
<br>
<pre>Caliba parameters:
<pre>Caliba parameters:
                 dmin=dmin                (default: 2.4)
                 dmin=dmin                (default: 2.4)
Line 52: Line 52:
This macro divides each peaklist into three classes of peaks with different calibration functions and constants:  
This macro divides each peaklist into three classes of peaks with different calibration functions and constants:  


<br>


 
{| rules="rows" cellspacing="0" cellpadding="0" border="1" style="width: 638px; height: 76px;" id="table1" class="twikiTable"
{| rules="rows" cellspacing="0" cellpadding="0" border="1" class="twikiTable" id="table1" style="width: 638px; height: 76px;"
|- class="twikiTableOdd twikiTableRowdataBgSorted0 twikiTableRowdataBg0"
|- class="twikiTableOdd twikiTableRowdataBgSorted0 twikiTableRowdataBg0"
! valign="top" bgcolor="#687684" class="twikiTableCol0 twikiFirstCol" | [http://www.nsm.buffalo.edu/Chem/HTP_twiki43/bin/view/NESG/NOECalibrationCYANA?sortcol=0;table=1;up=0#sorted_table <font color="#ffffff">class</font>]
! valign="top" bgcolor="#687684" class="twikiTableCol0 twikiFirstCol" | <font color="#ffffff">class</font>  
! valign="top" bgcolor="#687684" class="twikiTableCol1" | [http://www.nsm.buffalo.edu/Chem/HTP_twiki43/bin/view/NESG/NOECalibrationCYANA?sortcol=1;table=1;up=0#sorted_table <font color="#ffffff">peaks/constraints</font>]
! valign="top" bgcolor="#687684" class="twikiTableCol1" | <font color="#ffffff">peaks/constraints</font>  
! valign="top" bgcolor="#687684" class="twikiTableCol2 twikiLastCol" | [http://www.nsm.buffalo.edu/Chem/HTP_twiki43/bin/view/NESG/NOECalibrationCYANA?sortcol=2;table=1;up=0#sorted_table <font color="#ffffff">function</font>]
! valign="top" bgcolor="#687684" class="twikiTableCol2 twikiLastCol" | <font color="#ffffff">function</font>
|- class="twikiTableEven twikiTableRowdataBgSorted0 twikiTableRowdataBg0"
|- class="twikiTableEven twikiTableRowdataBgSorted0 twikiTableRowdataBg0"
| valign="top" bgcolor="#ffffff" class="twikiTableCol0 twikiFirstCol" |  
| valign="top" bgcolor="#ffffff" class="twikiTableCol0 twikiFirstCol" |  
backbone&nbsp;
backbone&nbsp;  


| valign="top" bgcolor="#ffffff" class="twikiTableCol1" |  
| valign="top" bgcolor="#ffffff" class="twikiTableCol1" |  
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; all HN/HA— HN/HA, and HN(i)/HA(i) — HB(j) with abs(i - j) &lt; 5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; all HN/HA— HN/HA, and HN(i)/HA(i) — HB(j) with abs(i - j) &lt; 5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>  


| valign="top" bgcolor="#ffffff" class="twikiTableCol2 twikiLastCol" |  
| valign="top" bgcolor="#ffffff" class="twikiTableCol2 twikiLastCol" |  
<code></code><code>V = A/d</code><sup><code>6</code></sup>
<code></code><code>V = A/d</code><sup><code>6</code></sup>  


|- class="twikiTableOdd twikiTableRowdataBgSorted1 twikiTableRowdataBg1"
|- class="twikiTableOdd twikiTableRowdataBgSorted1 twikiTableRowdataBg1"
| valign="top" bgcolor="#edf4f9" class="twikiTableCol0 twikiFirstCol" |  
| valign="top" bgcolor="#edf4f9" class="twikiTableCol0 twikiFirstCol" |  
sidechain
sidechain  


| valign="top" bgcolor="#edf4f9" class="twikiTableCol1" |  
| valign="top" bgcolor="#edf4f9" class="twikiTableCol1" |  
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; neither "backbone" nor "methyl"
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; neither "backbone" nor "methyl"  


| valign="top" bgcolor="#edf4f9" class="twikiTableCol2 twikiLastCol" |  
| valign="top" bgcolor="#edf4f9" class="twikiTableCol2 twikiLastCol" |  
<code>V = B/d<sup>4</sup></code><code></code>
<code>V = B/d<sup>4</sup></code><code></code>  


|- class="twikiTableEven twikiTableRowdataBgSorted0 twikiTableRowdataBg0"
|- class="twikiTableEven twikiTableRowdataBgSorted0 twikiTableRowdataBg0"
| valign="top" bgcolor="#ffffff" class="twikiTableCol0 twikiFirstCol twikiLast" |  
| valign="top" bgcolor="#ffffff" class="twikiTableCol0 twikiFirstCol twikiLast" |  
methyl
methyl  


| valign="top" bgcolor="#ffffff" class="twikiTableCol1 twikiLast" |  
| valign="top" bgcolor="#ffffff" class="twikiTableCol1 twikiLast" |  
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; all involving methyl groups
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; all involving methyl groups  


| valign="top" bgcolor="#ffffff" class="twikiTableCol2 twikiLastCol twikiLast" |  
| valign="top" bgcolor="#ffffff" class="twikiTableCol2 twikiLastCol twikiLast" |  
<code>V = C/d<sup>4</sup></code><code></code>
<code>V = C/d<sup>4</sup></code><code></code>  


|}
|}


By default, constant <code>A</code> is calculated from <code>avedis</code> - the assumed average distance between backbone atoms. It is also possible to specify the backbone constant <code>A</code> explicitly. Constant <code>B</code> is then calculated from <code>A</code>, and <code>C</code> from <code>B</code>  
By default, constant <code>A</code> is calculated from <code>avedis</code> - the assumed average distance between backbone atoms. It is also possible to specify the backbone constant <code>A</code> explicitly. Constant <code>B</code> is then calculated from <code>A</code>, and <code>C</code> from <code>B</code>  
Line 108: Line 110:
Type <code>help caliba</code> at CYANA prompt or see the DYANA manual for more information. You can also look at the script itself in the <code>macro</code> subdirectory of a CYANA installation.  
Type <code>help caliba</code> at CYANA prompt or see the DYANA manual for more information. You can also look at the script itself in the <code>macro</code> subdirectory of a CYANA installation.  


 
<br>
====

Revision as of 23:00, 13 November 2009

Introduction

NOE calibration converts NOESY cross-peak intensities into upper distance limits. A certain functional relationship is assumed between peak intensities and UPLs for a given group of peaks. The most common is NOE peak volume proportional to interproton distance to the minus six:

 V=A/d^6

where V is the peak volume (or intensity), d is the upper distance limit, and A is the calibration constant.


NOE Calibration Using CYANA 2.1

This section describes how to determine NOE calibration constants using CYANA 2.1.

Before running the calibration, the N chemical shifts in the peak list must be updated, and the peak lists separated into a list of N noes and C noes, using the UBNMR macro, macro1.

Command "calibrate"

A CYANA command calibrate is the basic calibration tool, unchanged since DYANA.


    Parameters: f(d)                      (required)
                dmin [d2 ...] dmax        (default: 2.4 5.5)
                weight=w                  (default: 1.0)
                plot=file  log  minimal   (default: none)

calibrate accepts a calibration function f(d) as the required argument.

This command uses "center" averaging. Explicit pseudoatom corrections are added to UPLs:

  • Multiplicity correction is applied by dividing the peak volume by the numbers of 1H spins in pseudoatoms assigned to the peak. For instance, the volume of a cross peak between a Leu QQD pseudo atom and a Tyr QD pseudo atom is divided by a factor of 6 x 2 = 12 prior to applying the calibration function. The resulting UPL is the subject to the upper and lower cutoffs.
  • Distance correction is applied by adding a distance between the pseudoatom and its constituent spins. It is applied after the application of upper and lower cutoffs. For example, for a Tyr QD pseudoatom this correction is equal to half the distance between the HD1 and HD2 spins. Distance correction can be manually performed with the distance correct command.

Type help calibrate at CYANA prompt for more information or see the DYANA documentation.


Macro "caliba"

CYANA macro caliba.cya is functionally the same as caliba.dya of DYANA. It is a more sophisticated calibration tool aimed at getting more nuanced calibration than a uniform 1/d6 dependence. CALIBA is recommended for manual structure calculation with CYANA 2.1. The caliba macro in fact uses calibrate, therefore, it also uses "center" averaging.


Caliba parameters:
                dmin=dmin                 (default: 2.4)
                dmax=dmax                 (default: 5.5)
                vmin=Vmin                 (default: 0.0)
                bb=A                      (default: calculated automatically)
                sc=B                      (default: A/dmin2)
                methyl=C                  (default: B/3)
                weight=w                  (default: 1.0)
                avedis=d                  (default: 3.4)
                plot=file                 (optional)

This macro divides each peaklist into three classes of peaks with different calibration functions and constants:


class peaks/constraints function

backbone 

       all HN/HA— HN/HA, and HN(i)/HA(i) — HB(j) with abs(i - j) < 5      

V = A/d6

sidechain

       neither "backbone" nor "methyl"

V = B/d4

methyl

       all involving methyl groups

V = C/d4


By default, constant A is calculated from avedis - the assumed average distance between backbone atoms. It is also possible to specify the backbone constant A explicitly. Constant B is then calculated from A, and C from B

Parameters dmin and dmax are the lower and upper cutoff limits for UPLs, respectively.

Parameter vmin is useful to filter out peaks with intensities below a certain threshold. Without it even incorrectly picked peaks with intensities below S/N level will produce UPLs set at the upper cutoff.

plot is an optional parameter, used to create logarithmic plots of volumes versus corresponding minimal distances in the selected structures together with the calibration functions. It is useful for fine-tuning the calibration during structure refinement. Of course, a structure must be then loaded into CYANA before running caliba.

VERY IMPORTANT!!!. The default behavior of structure calculation (anneal, calc_all, noeassign) in CYANA 2.x is to use "sum of r-6" averaging for UPLs. Setting expand=1 declares "sum of r-6" averaging, and expand=0 declares "center" averaging. If your UPLs reflect "center" averaging (that is, they were created with caliba or calibrate), then

  • They can be used with DYANA, CYANA 1.x and AutoStructure.
  • They can be used with legacy macros in CYANA 2.1, such as habas.
  • They can be used for manual structure calculation in CYANA 2.1, if expand=0 is set before calling anneal.
  • Do not use them for automated structure calculation in CYANA 2.1 with noeassign. The noeassign macro employs peaks calibrate (see below) and sets exapand=1. In general, do not mix them with UPLs calibrated with peaks calibrate.

Type help caliba at CYANA prompt or see the DYANA manual for more information. You can also look at the script itself in the macro subdirectory of a CYANA installation.