NOE Calibration Using CYANA

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.