Structure Calculation With RDC's Using CYANA

Introduction

IMPORTANT DISCLAIMER: A number of NESG NMR groups are currently beta-testing sites for CYANA-3.0 and the information in this page is intended for use by the licensed members of the NESG consortium, other beta testers, and is to be used in accordance to the program licensing agreement.

The following page describes the setup and analysis using RDC constraints in the framework of CYANA-3.0.

CYANA version 3.0 incorporates many new features including:

a) inclusion of residual dipolar coupling in structure calculation (RDC)

b) inclusion of pseudocontact shifts from paramagnetic centers in structure calculation (PCS)

c) improved function for symmetric dimers and annotation of intermolecular NOE contacts in the peaklist using the xeasy color code notation.  

Good agreement between global orientational constraints from RDC and distance information from NOEs is very important to achieve better structures by NMR. The program is an excellent step forward.  Results should always be accompanied by energy refinement in CNS or NIH-XPLOR.  Further information about the program and publication references can be found in the CYANA WIKI page.

Automated NOE and RDC and Structure Calculation Setup

The structure calculation with automated noesy assignments and RDC restraints follows the canonical CYANA recipe.  Simple annealing calcualtions starting from a set of constraints that include RDC are easily derived by simplifying the scripts below and following the demo scripts. The program requires a sequence file (name.seq), a proton assignment list (name.prot), a noesy peaklist set (name.peaks), an RDC list (name.rdc), a CALC.cya script and an init.cya script.  

The sequence file now includes the RDC tensor origin separated by dummy linker residues:

MET      1
THR      2
SER      3
THR      4
PHE      5
ASP      6
ARG      7
VAL      8
ALA      9
THR     10

PL     350
LL5    351
LL5    352
LL5    353
LL5    354
LL5    355
ORI    360

The RDC list supports multiple interatomic vectors in multiple media. RDC with distinct scaling factors and distinct ORI residue numbers are listed in a single file. The program supports the Da (magnitude) and R (Rhombicity) notation typical of programs such as PALES, REDCAT etc. Below is a sample RDC file that includes N-H, N-CA (intra), and N-C' (sequential) vectors in one medium with appropriate errors and scaling factors:

# Orientation  Magnitude  Rhombicity  ORI residue number
       1     5.39535        0.63125       360
#  First atom      Second atom                    RDC      Error  Weight  Orientation 
     7    ARG      H       7    ARG     N        5.936     2.000   1.000  1
     8    VAL      H       8    VAL     N        3.827     2.000   1.000  1
     9    ALA      H       9    ALA     N       -2.822     2.000   1.000  1
    10    THR      H      10    THR     N       -0.674     2.000   1.000  1
    11    ILE      H      11    ILE     N        4.945     2.000   1.000  1
    12    ILE      H      12    ILE     N        1.709     2.000   1.000  1
    13    ALA      H      13    ALA     N       -1.336     2.000   1.000  1
#
    4     THR      N       3    SER     C        1.267     0.095   8.330  1
    5     PHE      N       4    THR     C       -0.246     0.207   8.330  1   
    7     ARG      N       6    ASP     C        0.161     0.052   8.330  1   
    8     VAL      N       7    ARG     C        0.439     0.034   8.330  1   
    9     ALA      N       8    VAL     C       -0.076     0.040   8.330  1   
   10     THR      N       9    ALA     C       -0.957     0.048   8.330  1   
   11     ILE      N      10    THR     C        1.123     0.022   8.330  1   
   12     ILE      N      11    ILE     C       -0.440     0.037   8.330  1   
   13     ALA      N      12    ILE     C        0.065     0.026   8.330  1   
#
    3     SER      N      3     SER    CA        0.251     0.199   8.330  1
    4     THR      N      4     THR    CA       -0.265     0.258   8.330  1
    5     PHE      N      5     PHE    CA       -0.499     0.281   8.330  1
    7     ARG      N      7     ARG    CA       -0.457     0.200   8.330  1
    8     VAL      N      8     VAL    CA       -0.481     0.154   8.330  1
    9     ALA      N      9     ALA    CA        0.349     0.083   8.330  1
   10     THR      N     10     THR    CA        0.548     0.121   8.330  1
   11     ILE      N     11     ILE    CA       -0.091     0.111   8.330  1
   12     ILE      N     12     ILE    CA       -0.678     0.078   8.330  1
   13     ALA      N     13     ALA    CA        1.014     0.108   8.330  1

Multiple media (e.g. orientations) should be listed as follows:

# Orientation  Magnitude  Rhombicity  ORI residue number
       1     5.39535        0.63125       360
       2     7.55656        0.58200       370

and the sequence should be modified to include further links and ORI:

ALA      9
THR     10

PL     350
LL5    351
LL5    352
LL5    353
LL5    354
LL5    355
ORI    360
LL5    361
LL5    362
LL5    363
LL5    364
LL5    365
ORI    370

If the values of Da and R are not known they can be determined using any desired software OR using the FindTensor.cya script below. The program yields results that are equivalent to PALES assiuming the same fitting method is employed.

## 8DEMOS: FindTensor - Determine alignment tensor
##
## Determine magnitude and rhombicity of the alignment tensor
## from input RDCs


# determine tensor from histogram, no structure needed

read rdc phage_all_mono.rdc
print "    Input alignment tensor:"
do i 1 orientations
  print "    Orientation $i: magnitude = $magnitude(i) Hz, rhombicity = $rhombicity(i)."
end do

rdc fittensor method=simplex       # (can take several minutes)
#rdc fittensor method=gridsearch   # systematic search (very slow)


# alternatively, determine tensor from given structure by SVD

read rdc phage_all_mono.rdc
read pdb final.pdb
overview

The dummy values in the rdc list are read initially and they can be updated after running the FindTensor.cya routine. If no models are present (e.g. final.pdb) the program will terminate with a warning.

The working directory that contains all the files necessary to start the calcualtion is ready and the CALC.cya modified for the presence of RDC constraints is used:

peaks       := ali5.peaks,aro5.peaks,n3.peaks  # names of NOESY peak lists
prot        := RpR324.prot               # names of chemical shift lists
restraints  := ssa.cya,phage_all_mono.rdc # additional (non-NOE) constraints
tolerance   := 0.04,0.025,0.3            # chemical shift tolerances: HX2-HX1-X1
calibration :=                           # NOE calibration parameters
structures  := 100,20                    # number of initial, final structures
steps       := 10000                     # number of torsion angle dynamics steps
rmsdrange   := 10..80                    # residue range for RMSD calculation
randomseed  := 56231       # random number generator seed

weight_rdc   = 0.02               # weight for RDC restraints
cut_rdc      = 0.2                # cutoff for RDC violation output

ssa
noeassign peaks=$peaks prot=$prot autoaco

Notice the restraint contain the forced stereospecifically assigned methyls and sidechain NH's and the RDC set.  Also, the NOE vs. RDC weight is set by the weight_rdc and cut_rdc functions. The remaining instructions are identical to the CYANA-2.1 file.  Further parameters are specified in the init.cya file below:

name:=RpR324
rmsdrange:=10-80
cyanalib
nproc:=8
read seq $name.seq
rdcdistances

The above script is intended to run off a single dual-quad core machine (nproc=8).  Please note the rdcdistance.cya macro is being called by the init.cya setup file.  This file, located in the cyana-3.0/macro directory contains the supported RDC vectors, more vectors could potentially be added such as Trp NE1-HE1 that maybe useful in deuterated sample to direct the large hydrophobic sidechain. 

# Copyright (c) 2002-08 Peter Guntert. All rights reserved.
## 7MACROS: rdcdistances - CYANA macro
##
## Parameters: (none)
##
# dipole definition format: atom1_name atom2_name atom1_index atom2_index
# if indexes are missing, zeros are assumed

var info echo

syntax

info:=none; echo:=off
rdc distance "N  H"    distance=1.041
rdc distance "CA HA"   distance=1.117
rdc distance "C  CA"   distance=1.525
rdc distance "C  N"    distance=2.461
rdc distance "C  N -1" distance=1.329
rdc distance "CA N"    distance=1.458
rdc distance "CA N -1" distance=2.425
rdc distance "CA H"    distance=2.117
rdc distance "CA H -1" distance=2.533
rdc distance "C  H -1" distance=2.000
rdc distance "C  HA"   distance=2.144
rdc distance "CB HB"   distance=1.080
rdc distance "CA CB"   distance=1.532
unset info
print "    Standard RDC distances defined."

The command line execution form single machine or cluster using MPI implememntation (highly recommended) is carried on as usual:

/cyana-3.0/cyana CALC > & log &

the MPI is launched using the script called, for example, submit_cyana:

#!/bin/bash
#PBS -S /bin/bash
#PBS -N cyana
#PBS -lnodes=6:ppn=8
lamboot ~/bhost.def
cd /farm/users/prossi/RpR324_structure/cyana_new_mono2
/opt/openmpi/tcp-gnu/bin/mpirun /farm/software/cyana-3.0-mpi/cyana CALC.cya
lamhalt

with the command:

qsub -q @master3 submit_cyana

The starting scripts are highly system specific they are almost guaranteed NOT to work on your sytem and are given here for general information only.

Output analysis

The output analysis is carried out in the usual manner, it should be noted that, during the calculation the specified values Da and R are kept fixed.  Following the final cycle a new model-based estimate of Da and R is calculated and used to compute the RDC violations and their contribution to the target function and the quality factor (Q). The resulting target function will be increased by the number and extent of RDC violations in addition to other violations from dihedral, vdw, and NOEs restraints.

A partial output file is given below (final.ovw):

 
    Structural statistics:
 
    str   target     upper limits    van der Waals             RDCs
        function   #    rms   max   #    sum   max   #    rms   max
      1     8.98   5 0.0089  0.40  12   18.5  0.37  13 0.2508  2.27
      2     8.96   7 0.0083  0.31  15   17.4  0.38  19 0.2656  2.27
      3     9.69  11 0.0122  0.54  15   18.6  0.38  18 0.2595  2.28
      4     9.63   6 0.0068  0.20  16   18.9  0.38  18 0.2604  2.27
      5     9.41  14 0.0089  0.22  17   20.4  0.37  14 0.2521  2.26
      6     9.80  10 0.0085  0.24  17   19.7  0.38  15 0.2563  2.27
      7    10.52  12 0.0158  0.78  15   19.1  0.38  13 0.2554  2.29
      8    10.05  10 0.0084  0.18  18   20.3  0.42  18 0.2587  2.28
      9    10.57   7 0.0079  0.29  19   18.3  0.53  16 0.3169  2.43
     10    10.38  13 0.0088  0.21  22   20.5  0.37  21 0.2636  2.28
     11    10.33   7 0.0067  0.16  17   19.6  0.64  15 0.2663  2.28
     12     9.93   9 0.0090  0.29  19   20.6  0.38  20 0.2541  2.27
     13    10.12  12 0.0098  0.23  21   19.9  0.47  13 0.2681  2.30
     14    10.53   7 0.0077  0.19  23   19.3  0.38  11 0.3145  2.40
     15    10.96   8 0.0089  0.31  21   22.6  0.38  12 0.2639  2.28
     16    10.50  12 0.0100  0.29  23   20.1  0.37  15 0.2696  2.27
     17    10.56  14 0.0119  0.30  23   21.1  0.37  16 0.2736  2.33
     18    10.75  17 0.0146  0.60  21   20.0  0.37  12 0.2547  2.27
     19    10.88  18 0.0125  0.38  19   23.3  0.38  10 0.2506  2.30
     20    10.88   8 0.0083  0.26  23   21.4  0.42  15 0.2644  2.28
 
    Ave    10.17  10 0.0097  0.32  19   20.0  0.41  15 0.2660  2.29
    +/-     0.59   4 0.0024  0.15   3    1.4  0.07   3 0.0177  0.04
    Min     8.96   5 0.0067  0.16  12   17.4  0.37  10 0.2506  2.26
    Max    10.96  18 0.0158  0.78  23   23.3  0.64  21 0.3169  2.43
    Cut                      0.10             0.20             0.20
 
    Constraints violated in 6 or more structures:
                                                   #   mean   max.  1   5   10   15   20
    Upper HA    PRO   19 - HB3   ARG   20   5.50  15   0.11   0.21  +++++++  ++++   ++*+  peak 970
    Upper HA    ILE   23 - QB    SER   27   5.34   6   0.09   0.16         + ++   +  +*   peak 276
    Upper HA    ILE   30 - HB2   LEU   33   4.95  10   0.10   0.15      ++++  ++  ++ +*   peak 313
    VdW   CB    ALA   69 - H     THR   70   2.55  18   0.24   0.28  ++++++++ ++++ *+++++
    VdW   O     THR   71 - N     PHE   75   2.75  12   0.20   0.32   ++* +++  +   ++ +++
    VdW   O     PHE   75 - C     VAL   76   2.80  13   0.21   0.31      ++++ +  *+++++++
    VdW   CG1   VAL   76 - HG2   LYS   78   2.60   9   0.17   0.25   ++*   +++++   +
    VdW   HG3   LYS   78 - C     LYS   78   2.50  14   0.22   0.36  +++*++ +++++ + +  +
    Ori 1 N     ALA   69 - CA    ALA   69  -2.84  20   1.95   2.04  +++++++++++++*++++++
    Ori 1 H     ALA   69 - N     ALA   69  -3.32  20   2.01   2.33  ++++++++++++++++*+++
    Ori 1 N     THR   70 - C     ALA   69   0.25  16   0.31   1.06  ++++++++*+++ + +  ++
    Ori 1 H     ASN   79 - N     ASN   79  -4.71  10   0.17   0.32  + *+ +  ++++   +  +
    Ori 1 H     GLY   92 - N     GLY   92  -5.54   8   0.15   0.51     ++++*       +  ++
    Ori 1 H     LEU   94 - N     LEU   94   0.42   6   0.14   0.54   +*+     +++
    3 violated distance restraints.
    5 violated van der Waals restraints.
    6 violated residual dipolar coupling restraints.
 
 
    RDC statistics:
    Correlation coefficient      :    0.906 +/-  0.003    (0.899..0.909, best in conformer 4)
    Q = rms(Dcalc-Dobs)/rms(Dobs):   42.709 +/-  0.583 %  (42.073..44.194)
    Q normalized by tensor       :   32.943 +/-  0.638 %  (32.231..34.776)
    Alignment tensor magnitude   :    5.881 +/-  0.045 Hz (5.760..5.950, best 5.898; input 5.898)
    Alignment tensor rhombicity  :    0.537 +/-  0.007    (0.524..0.556, best 0.539; input 0.539)
 
 

-- PaoloRossi - 14 Dec 2009