NOE Calibration Using CYANA: Difference between revisions

From NESG Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 2: Line 2:


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>  
<br>


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

Revision as of 23:05, 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.


Command "peaks calibrate"

Command peaks calibrate is a new feature of CYANA 2.x, which produced a "sum of r-6"-averaged calibration. The documentation for it is scarce. For example, the manual structure calculation script ~/demo/manual/CALC.cya contains these lines:


peaks      := c13,n15,aro            # names of peak lists
prot       := demo                   # names of proton lists
tolerance  := 0.040,0.030,0.45       # chemical shift tolerances
                                     # order: 1H(a), 1H(b), 13C/15N(b), 13C/15N(a)
calibration:= 6.7E5,8.2E5,8.0E4      # calibration constants (will be determined
                                     #   automatically, if commented out)
dref       := 4.2                    # average upper distance limit for
                                     #   automatic calibration
  
  calibration prot=$prot peaks=$peaks constant=$calibration dref=$dref
  peaks calibrate "**" simple
  write upl $name-in.upl

Line calibration prot=$prot peaks=$peaks constant=$calibration dref=$dref calls the macro calibration.cya, which iteratively sets calibration constants for each peaklist. If the calibration list variable is defined, its values are taken as calibration constants. If it is left blank, calibration constants are calibrated based on assumed average distance dref between all atoms of the molecule. Note that dref here is inherently different from avedis of caliba - the latter refers only to distances between backbone protons.

Line peaks calibrate "**" simple generates UPLs assuming simple C/d^6 functions with calibration constants C from the previous step. These calibration constants are thus closely related the "backbone" class constants A of caliba.

To modify upper and lower cutoffs for NOE calibration with peaks calibrate add a line set upl_values:=2.4,6.0. The defaults are 2.4 and 5.5.

The main difference from calibrate and caliba is that pseudoatom correction are not applied to resulting UPLs.


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 you UPLs use "sum of r-6" averaging (if created with peaks calibrate), then

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


The aim of this new calibration method is likely to facilitate automated NOE assignment and structure calculation with noeassign. Unlike caliba and calibrate, peaks calibrate effectively generates UPLs even for unassigned peaks assuming NOEs between single protons. If such a peak gets assigned to a pseudoatom, the required correction is applied implicitly during the calculation.