HSQCTROSY RDC Measurement: Difference between revisions

From NESG Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(22 intermediate revisions by the same user not shown)
Line 1: Line 1:
&nbsp;'''<span style="font-size:14.0pt;mso-bidi-font-size:12.0pt">NMR Data Processing </span>''''''<span style="font-size:14.0pt">&gt;<span style="mso-spacerun: yes">&nbsp; </span>RDCs</span>'''<!--StartFragment-->  
'''NMR Data Processing &gt;&nbsp; 2D HSQC-TROSY for RDCs Measurement''' <br>  


<font class="Apple-style-span" size="6"><span class="Apple-style-span" style="font-size: 19px;">'''
===== '''Theory/Background''' =====
'''</span></font>


'''2D HSQC-TROSY'''
Residual Dipolar Couplings (RDCs) originate from the anisotropic component of the dipolar interaction, which is dependent on the angle between an internuclear vector and the magnetic field. When a molecule samples orientations uniformly, as it does in normal solution NMR, RDCs average to zero and are not observable. However, if a molecule is dissolved in a dilute liquid crystalline medium, or other orienting medium, it becomes partially aligned, and the dipolar couplings are not completely averaged to zero.&nbsp; This leads to a small contribution to the splitting of NMR signals. The angular dependence of these contributions can provide valuable structural information.&nbsp; RDCs can be used to validate protein structures, to refine structures to improve quality, and to provide constraints as a part of an initial structure determination<span style="font: 12.0px Times New Roman">.</span>&nbsp; RDCs can also be used to determine the rotational symmetry axis of oligomers and this information can be use to properly orient the positions of subunits within an oligomeric complex.&nbsp;


'''<br>'''
<br>  


===== '''1.1<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Theory/Background''' =====
===== '''Recent Reviews''' =====


<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">Residual Dipolar Couplings (RDCs) originate from the anisotropic
Prestegard, Bougault&amp; Kishore, ''Chemical Reviews'', '''104''', 3519-3540 (2004)&nbsp;  
component of the dipolar interaction, which is dependent on the angle between
an internuclear vector and the magnetic field. When a molecule samples
orientations uniformly, as it does in normal solution NMR, RDCs average to zero
and are not observable. However, if a molecule is dissolved in a dilute liquid
crystalline medium, or other orienting medium, it becomes partially aligned,
and the dipolar couplings are not completely averaged to zero.<span style="mso-spacerun: yes">&nbsp; </span>This leads to a small contribution to the splitting of NMR signals. The angular dependence of these contributions can provide valuable structural information. <span style="mso-spacerun: yes">&nbsp;</span>RDCs can be used to validate protein structures, to refine structures to improve quality, and to provide constraints as a part of an initial structure determination.<span style="mso-spacerun: yes">&nbsp; </span>RDCs can also be used to determine the rotational symmetry axis of oligomers and this information can be use to properly orient the positions of subunits within an oligomeric complex.</span>


<br>
Lipsitz&amp; Tjandra, ''Ann. Rev. Biophys. Biomol. Struct''., '''33''', 387-413 (2004)&nbsp;


===== '''1.2<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Recent Reviews'''<br> =====
Additional information can also be found at the following web site.  


<span style="mso-bidi-font-size: 28.0pt">Prestegard, Bougault&amp; Kishore, ''Chemical Reviews'', '''104''',
<span style="text-decoration: underline">[http://tesla.ccrc.uga.edu/courses/bionmr/lectures http://tesla.ccrc.uga.edu/courses/bionmr/lectures]</span>
3519-3540 (2004)</span>


<span style="mso-bidi-font-size: 28.0pt">Lipsitz&amp; Tjandra, ''Ann. Rev. Biophys. Biomol. Struct''., '''33''',
===== '''Brief Experimental Description''' =====
387-413 (2004)</span><span style="mso-bidi-font-family:Helvetica">&nbsp;</span>


<span style="mso-bidi-font-family:Helvetica" />Additional information can also be found at the following web site.
<span style="font: 12.0px Times New Roman">'''<span class="Apple-tab-span" style="white-space:pre"> </span>'''</span>HSQC-TROSY pairs provide one way of measuring RDCs.&nbsp; The procedure is applicable to proteins of small and moderate size (10-30 kDa).&nbsp; The acquisition time in the indirect t1 should be at least T2 of nitrogen.&nbsp; The measurements are easily visualized as displacements of peaks in 2D frequency domain spectra.&nbsp; The HSQC member of the pair gives a single peak at the center of the four peaks expected in a fully coupled HSQC experiment. The TROSY member gives only the sharpest component. The offset in either H or N direction is one half (J).&nbsp; Procedures for collecting this pair of experiments are described elsewhere in this Twiki (Experimental setup for HSQC-TROSY).


[http://tesla.ccrc.uga.edu/courses/bionmr/lectures http://tesla.ccrc.uga.edu/courses/bionmr/lectures]
===== '''NMR Data Processing'''  =====


'''<br>'''
The data processing described below is specific for data collected on Varian spectrometers having the BioPack (version: 2009<span style="font: 12.0px Times New Roman">-0</span>9<span style="font: 12.0px Times New Roman">-0</span>3) sequence (gNhsqc.c) installed.&nbsp; Or on Bruker spectrometers having the …….. It is applied to both isotropic and aligned HSQC-TROSY pairs.


===== '''1.3<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Brief Experimental Description'''<br> =====
===== Software Information  =====


'''<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>'''HSQC-TROSY pairs provide one way of measuring RDCs.<span style="mso-spacerun: yes">&nbsp;
NMRPipe (download information and user manual) <br>[http://spin.niddk.nih.gov/NMRPipe/ http://spin.niddk.nih.gov/NMRPipe/] <br> <br> Brief descriptions of specific functions are accessible via nmrPipe –
</span>The procedure is applicable to proteins of small and moderate size (10-30 kDa).<span style="mso-spacerun: yes">&nbsp; </span>The acquisition time in the indirect t1 should be at least T2 of nitrogen.<span style="mso-spacerun: yes">&nbsp; </span>The measurements are easily visualized as displacements of peaks in 2D frequency domain spectra. <span style="mso-spacerun: yes">&nbsp;</span><span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">The HSQC member of the pair gives a single peak at the center of the
four peaks expected in a fully coupled HSQC experiment. The TROSY member gives
only the sharpest component. The offset in either H or N direction is one half
(J). <span style="mso-spacerun: yes">&nbsp;</span>Procedures for collecting this pair of experiments are described elsewhere in this Twiki (Experimental setup for HSQC-TROSY).</span>


<br>
:‘nmrPipe –help’ will list most functions.
:‘nmrPipe –fn GM –help’ will give description of the GM function.


===== '''2.1<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>NMR Data Processing''' =====
<br> Supported Platforms


<span style="mso-bidi-font-size: 14.0pt;mso-bidi-font-family:Arial">The data processing described below is
:Linux (RedHat Linux/Fedora)
specific for data collected on Varian spectrometers having the BioPack
:Mac OS X (10.3.4 and up)  
(version: 2009-09-03) sequence (gNhsqc.c) installed.<span style="mso-spacerun: yes">&nbsp; </span>Or on Bruker spectrometers having the …….. It is applied to both isotropic and aligned HSQC-TROSY pairs.</span>
:SGI Irix &nbsp;(6.2 and up)
:Sparc Solaris (2 and up)
:Windows XP Pro with Microsoft Services for UNIX (SFU 3.5)


'''<br>'''
===== '''Varian to NMRPipe Data Conversion''' =====
<pre>#!/bin/csh
##conversion of varian 2d fids into nmrPipe format for a TROSY/hsqc data set
##with array= ‘(TROSY, dm2),phase’
var2pipe -in ./fid \
  -xN          2048    -yN      1024  \
  -xT          1024    -yT      256  \
  -xMODE    Complex    -yMODE  Rance-Kay  \
  -xSW      10000.00  -ySW      2800.00  \
  -xOBS      599.520    -yOBS      60.756  \
  -xCAR        4.773    -yCAR    118.262  \
  -xLAB          HN    -yLAB        N  \
  -ndim            2    -aq2D    States  \
  | nmrPipe -fn QMIX -ic 4 -oc 2 -cList 1 0 0 1 0 0 0 0  \
  -out trosy.fid -verb -ov


'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">2.2<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Required Software:&lt;o:p&gt;&lt;/o:p&gt;</span>'''  
var2pipe -in ./fid \
  -xN          2048    -yN      1024  \
  -xT          1024    -yT      256  \
  -xMODE    Complex    -yMODE  Rance-Kay  \
  -xSW      10000.00  -ySW    2800.00  \
  -xOBS      599.520    -yOBS      60.756  \
  -xCAR        4.773    -yCAR    118.262  \
  -xLAB          HN    -yLAB        N  \
  -ndim            2   -aq2D    States  \
  | nmrPipe -fn QMIX -ic 4 -oc 2 -cList 0 0 0 0 1 0 0 1   \
  -out hsqc.fid -verb –ov
</pre>  
'''Cautions:''' The usual Varian to NMRPipe macros will read in experimental parameters, but will not be able to perform data conversion since both the HSQC and TROSY spectra are collected in a interleaved manner and are stored in a single FID data file.&nbsp; The above script includes the QMIX function to separate fids.


<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
===== '''Bruker to NMRPipe Data Conversion'''  =====


<span style="mso-bidi-font-size: 14.0pt;mso-bidi-font-family:Arial">NMRPipe &lt;o:p&gt;&lt;/o:p&gt;</span>
The conversion of data collected on most Bruker spectrometers can be accomplished with the following script.  


<span style="mso-bidi-font-size: 14.0pt;mso-bidi-font-family:Arial">[http://spin.niddk.nih.gov/NMRPipe/ http://spin.niddk.nih.gov/NMRPipe/]&lt;o:p&gt;&lt;/o:p&gt;</span>
(insert script)


<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
===== '''Fourier Transformation of Converted HSQC and TROSY spectra'''  =====


'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">2.3<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Supported Platforms:&lt;o:p&gt;&lt;/o:p&gt;</span>'''
The following script contains some data specific parameters including phasing parameters and convolution functions that will need to be optimized for each specific set.
<pre>#!/bin/csh
#This script is to process one of the TROSY/
| nmrPipe  -in hsqc.fid              \
| nmrPipe -fn SOL  \
| nmrPipe  -fn SP -off 0.5 -end 0.99 -pow 2 -c 0.5  \
| nmrPipe  -fn ZF -auto                    \
| nmrPipe  -fn FT                    \
| nmrPipe -fn PS -p0 -61.6 -p1 -4.6  -di    \
| nmrPipe  -fn EXT -x1 6.00ppm -xn 11.50ppm -verb 2 -sw  \
| nmrPipe -fn TP                      \
| nmrPipe -fn LP -x1 1 -xn 256 -pred 256 -after -verb  \
| nmrPipe -fn SP -off 0.5 -end 0.99 -pow 2 -c 1.0  \
| nmrPipe -fn ZF -auto          \
| nmrPipe -fn FT                  \
| nmrPipe -fn PS -p0 -91.8 -p1 7.0  -di    \
| nmrPipe -fn TP        \
| nmrPipe -fn EXT -y1 100.00ppm -yn    138.00ppm -sw \
| nmrPipe -out hsqc.ft2 -verb 2 -ov
</pre>


&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
===== '''RDC Assignment and Calculation'''  =====


<span style="mso-bidi-font-size: 16.0pt;mso-bidi-font-family:Arial">Linux (RedHat Linux/Fedora) &lt;o:p&gt;&lt;/o:p&gt;</span>
'''General Concepts'''&nbsp;  


<span style="mso-bidi-font-size: 16.0pt;mso-bidi-font-family:Arial">Mac OS X (10.3.4 and up)&lt;o:p&gt;&lt;/o:p&gt;</span>
Transfer all assignments (or cross peak labels) from a reference isotropic HSQC spectrum to TROSY and HSQC spectra obtained under both isotropic and aligned conditions.&nbsp; The positions of&nbsp; cross-peaks in the aligned HSQC spectrum usually change little with addition of an alignment medium; if they do change, there is often a constant offset for all peaks due to changes in the lock resonance<span style="font: 12.0px Times New Roman">.</span>&nbsp; Hence<span style="font: 12.0px Times New Roman">,</span> transferring assignments is usually straightforward. Corresponding peaks for TROSY assignments are diagonally displaced (down and to the right) by 38-58 Hz.&nbsp; The isotropic coupling (J) is calculated from the difference in Hz&nbsp; between the <sup>15</sup>N dimension of isotropic HSQC and TROSY spectra (the offset is one half (J), so the final value needs to be multiply by 2).&nbsp; Similarly, the aligned coupling (J+D) is calculated from differences in peak positions for the aligned HSQC and TROSY spectra.&nbsp; The Hz values of the corresponding peaks can be subtracted using a spreadsheet or script to calculate the RDC values.


<span style="mso-bidi-font-size: 16.0pt;mso-bidi-font-family:Arial">SGI Irix <span style="mso-spacerun: yes">&nbsp;</span>(6.2 and up)&lt;o:p&gt;&lt;/o:p&gt;</span>
The peak picking and assignment process can be done using any graphical NMR assignment program of user’s choice.&nbsp; NMRViewJ has a special feature that would match and transfer assignment from one spectrum to another. &nbsp; For this tutorial, we have included a Tcl/TK script for obtaining RDCs using NMRDraw and NMRWish; the process is broken down into five steps.


<span style="mso-bidi-font-size: 16.0pt;mso-bidi-font-family:Arial">Sparc Solaris (2 and up)&lt;o:p&gt;&lt;/o:p&gt;</span>  
<br>  


<span style="mso-bidi-font-size: 16.0pt;mso-bidi-font-family:Arial">Windows XP Pro with Microsoft Services for
===== '''RDC Assignment and Calculation Scripts'''  =====
UNIX (SFU 3.5).</span>&lt;o:p&gt;&lt;/o:p&gt;


&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
'''Step 1: step1_hsqc_ass.tcl''' <br>


'''2.4<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><sup></sup>Varian to NMRPipe data conversion&lt;o:p&gt;&lt;/o:p&gt;'''
'''Function:''' Automatically pick peaks in NMRDraw for the isotropic HSQC spectrum and save the peak list as hsqc.tab to be use as reference peak list.&nbsp; This script cuts out vast amounts of&nbsp; non-essential data that are stored in the hsqc.tab file after initial peak picking, leaving only the labels and, HN and N ppm values.


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>  
<br>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">#!/bin/csh&lt;o:p&gt;&lt;/o:p&gt;</span>  
'''Execution:''' Executing this script does not require code modification. The only requirement is that you add the output and input script names in the execution line of a LINUX operating system<span style="font: 12.0px Times New Roman">.</span>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"># conversion of
<br>  
varian 2d fids into nmrPipe format for a TROSY/hsqc data set&lt;o:p&gt;&lt;/o:p&gt;</span>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"># with array= ‘(TROSY,
'''Command:&nbsp;'''
dm2),phase’&lt;o:p&gt;&lt;/o:p&gt;</span>  
<pre>./step1_hsqc_ass.tcl hsqc_ass.tab hsqc.tab</pre>  
Where “hsqc_ass.tab” is the output file and “hsqc.tab” is the input file to this script.


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">var2pipe -in
<br>  
./fid \&lt;o:p&gt;&lt;/o:p&gt;</span>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-xN<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Download: [[Media:Step1_hsqc_ass.tcl|Step1_hsqc_ass.tcl]] &nbsp;  
</span>2048<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yN<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>1024<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span></span><span lang="FR" style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700;mso-ansi-language:FR">-xT<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<br>  
</span>1024<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yT<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>256<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>  


<span lang="FR" style="font-size:8.0pt;mso-bidi-font-size: 11.0pt;font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700; mso-ansi-language:FR"><span style="mso-spacerun: yes">&nbsp;&nbsp;
'''Step 2: i-hsqctrosy_step2.tcl'''
</span>-xMODE<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;
</span>Complex<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;
</span>-yMODE<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>Rance-Kay<span style="mso-spacerun: yes">&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>


<span lang="FR" style="font-size:8.0pt;mso-bidi-font-size: 11.0pt;font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700; mso-ansi-language:FR"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span></span><span style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700">-xSW<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>10000.00<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-ySW<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>2800.00<span style="mso-spacerun: yes">&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
'''Function:''' This script takes the HSQC and TROSY peak tables that you created while peak picking and combines them into a single table. It also transfers the assignment from the isotropic HSQC to HSQC aligned, TROSY isotropic and TROSY aligned spectra. On execution, it opens “NMRWish” which can further help in transferring the assignment peak by peak from HSQC to HSQC-TROSY. All the blue peaks in the NMRWish display belong to the HSQC spectrum and the white peaks belong to the TROSYspectrum. The peak under consideration is on the crosswire of the NMRWish window.&nbsp; A navigation bar will also appear along with the NMRWish display to allow user to visually inspect the assignment of each peak and make corrections if necessary.  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-xOBS<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>599.520<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yOBS<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>60.756<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>  
[[Image:Step2-spectrum.png|Step2-spectrum.png]] [[Image:Step2-navi.png|Step2-navi.png]] <br>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-xCAR<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>4.773<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yCAR<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>118.262<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
'''Execution:''' Executing this script requires editing the script file itself and changing a few values. Make sure that '''assName '''(from step 1)''', tabName1'''(peak table for HSQC)''', tabName2''' (peak table for TROSY)''', and outName '''(output file name) have the correct input table names. You must also enter the height (contour level) at which you picked the peaks for hsqc in the '''hi1 '''field and the height at which you picked peaks for trosy in the '''hi2''' field. Also make sure that the correct names for the .ft2 spectrum are specified in the script.&nbsp;  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-xLAB<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;  
</span>HN<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yLAB<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
</span>N<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-ndim<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
'''Command:&nbsp;'''
</span>2<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-aq2D<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>States<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>  
<pre>./i-hsqctrosy_step2.tcl</pre>  
<br>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>| nmrPipe -fn QMIX -ic 4 -oc 2 -cList 1 0 0 1 0 0 0 0<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
Download: [[Media:I-hsqctrosy_step2.tcl|I-hsqctrosy_step2.tcl]]


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span></span><span lang="NO-BOK" style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700;mso-ansi-language:NO-BOK">-out
<br>  
trosy.fid -verb -ov&lt;o:p&gt;&lt;/o:p&gt;</span>  


<span lang="NO-BOK" style="font-size:8.0pt;mso-bidi-font-size: 11.0pt;font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700; mso-ansi-language:NO-BOK">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
'''Step 3: i-hsqctrosy_step3.tcl'''


<span lang="NO-BOK" style="font-size:8.0pt;mso-bidi-font-size: 11.0pt;font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700; mso-ansi-language:NO-BOK">var2pipe -in ./fid \&lt;o:p&gt;&lt;/o:p&gt;</span>
'''Function:''' This script calculates the coupling values for the protein.  


<span lang="NO-BOK" style="font-size:8.0pt;mso-bidi-font-size: 11.0pt;font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700; mso-ansi-language:NO-BOK"><span style="mso-spacerun: yes">&nbsp; </span></span><span style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700">-xN<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<br>  
</span>2048<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yN<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>1024<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span></span><span lang="FR" style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700;mso-ansi-language:FR">-xT<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
'''Execution:''' Execution of this script requires that you open up the script and change the '''inName''' to the output of step 2. It also requires that you once again input the '''hi1''' and '''hi2''' values just like you did in step 2.  
</span>1024<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yT<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>256<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>


<span lang="FR" style="font-size:8.0pt;mso-bidi-font-size: 11.0pt;font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700; mso-ansi-language:FR"><span style="mso-spacerun: yes">&nbsp;&nbsp;
<br>  
</span>-xMODE<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;
</span>Complex<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;
</span>-yMODE<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>Rance-Kay<span style="mso-spacerun: yes">&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>  


<span lang="FR" style="font-size:8.0pt;mso-bidi-font-size: 11.0pt;font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700; mso-ansi-language:FR"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span></span><span style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700">-xSW<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>10000.00<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-ySW<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp; </span>2800.00<span style="mso-spacerun: yes">&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>  
'''Command:'''
<pre>./i-hsqctrosy_step3.tcl</pre>  
<br>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-xOBS<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>599.520<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yOBS<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>60.756<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
Download: [[Media:I-hsqctrosy_step3.tcl|I-hsqctrosy_step3.tcl]]


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-xCAR<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>4.773<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yCAR<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>118.262<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>  
<br>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-xLAB<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
'''Step 4: i-hsqctrosy_step4.tcl'''
</span>HN<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-yLAB<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
</span>N<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>-ndim<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
'''Function:''' This script cleans up the output from the previous script (step 3) and outputs the necessary information (chemical shift of H and N in ppm and coupling values) into a new text file.&nbsp;  
</span>2<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>-aq2D<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>States<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp;&nbsp; </span>| nmrPipe -fn QMIX -ic 4 -oc 2 -cList 0 0 0 0 1 0 0 1<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>  
<br>  


<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-spacerun: yes">&nbsp; </span><span style="mso-spacerun: yes">&nbsp;</span>-out hsqc.fid -verb –ov&lt;o:p&gt;&lt;/o:p&gt;</span>
'''Execution:''' Execution of this script requires that you change a single variable in the script. Make sure that you set the '''outName''' variable to the proper media type (gel, peg, iso, etc…).&nbsp;  


'''<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family:Courier">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>'''
<br>  


'''<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family:Courier">Cautions:</span>'''<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family:Courier"> The usual
'''Command:'''  
Varian to NMRPipe macros will read in experimental parameters, but will not be
<pre>./i-hsqctrosy_step4.tcl</pre>  
able to perform data conversion since both the HSQC and TROSY spectra are
<br>  
collected in a interleaved manner and are stored in a single FID data
file.<span style="mso-spacerun: yes">&nbsp; </span>The above script includes the QMIX function to separate fids.&lt;o:p&gt;&lt;/o:p&gt;</span>  


<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family: Courier">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
Download: [[Media:I-hsqctrosy_step4.tcl|I-hsqctrosy_step4.tcl]]


'''2.5<sup><span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span></sup>Bruker to NMRPipe data conversion&lt;o:p&gt;&lt;/o:p&gt;'''
<br>  


<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family: Courier">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
'''Step 5: i-hsqctrosy_step5.tcl'''


<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family: Courier"><span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>The conversion of data collected on most Bruker spectrometers can be accomplished with the following script.<span style="mso-spacerun: yes">&nbsp; </span>The same precautions pertain?&lt;o:p&gt;&lt;/o:p&gt;</span>
'''Function:''' This script calculates the rdc values for the protein.  


<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family: Courier">(insert script)&lt;o:p&gt;&lt;/o:p&gt;</span>  
<br>  


<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family: Courier">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
'''Execution:''' We recommend that a new folder be created in the protein directory titled '''rdc_name?''' to organize and retain all of the files required to run this script. Copy or Move all of the required files into the '''rdc '''folder; these include: '''hsqc_ass.tab''' file and all of the '''*_j.tab''' files generated from isotropic data and each medium. Make sure that the '''assName, isoName, '''and '''othName '''We recommend the '''outName''' field to have the naming convention '''protein.name_rdc_media.tab''' (ex. rnr77_rdc_gel).


'''<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family:Courier">2.6<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Fourier transformation of converted HSQC and TROSY spectra&lt;o:p&gt;&lt;/o:p&gt;</span>'''  
'''&nbsp;'''  


'''<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>'''
'''Command:'''
<pre>./i-hsqctrosy_step5.tcl</pre>  
<br>  


'''<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700"><span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span></span>'''<span style="mso-bidi-font-family:Courier;color:#7D4700">The following script
Download: [[Media:I-hsqctrosy_step5.tcl|I-hsqctrosy_step5.tcl]]
contains some data specific parameters including phasing parameters and
<br>
convolution functions that will need to be optimized for each specific set.&lt;o:p&gt;&lt;/o:p&gt;</span>
<br>
 
<br>
<span style="mso-bidi-font-family:Courier;color:#7D4700">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
Updated by Hsiau-Wei Lee, 2011
 
'''<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt;font-family:Courier; mso-bidi-font-family:Courier;color:#7D4700">#!/bin/csh&lt;o:p&gt;&lt;/o:p&gt;</span>'''
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700"># This script
is to process one of the TROSY/&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe<span style="mso-spacerun: yes">&nbsp; </span>-in hsqc.fid<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
</span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
SOL<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe<span style="mso-spacerun: yes">&nbsp; </span>-fn SP -off 0.5 -end 0.99 -pow 2 -c 0.5<span style="mso-spacerun: yes">&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe<span style="mso-spacerun: yes">&nbsp; </span>-fn ZF -auto<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
</span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe<span style="mso-spacerun: yes">&nbsp; </span>-fn FT<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
</span><span style="mso-spacerun: yes">&nbsp;&nbsp;</span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
PS -p0 -61.6 -p1 -4.6<span style="mso-spacerun: yes">&nbsp; </span>-di<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe<span style="mso-spacerun: yes">&nbsp; </span>-fn EXT -x1 6.00ppm -xn 11.50ppm -verb 2 -sw<span style="mso-spacerun: yes">&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
TP<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
</span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
LP -x1 1 -xn 256 -pred 256 -after -verb<span style="mso-spacerun: yes">&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
SP -off 0.5 -end 0.99 -pow 2 -c 1.0<span style="mso-spacerun: yes">&nbsp;
</span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
ZF -auto<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
FT<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
</span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
PS -p0 -91.8 -p1 7.0<span style="mso-spacerun: yes">&nbsp;&nbsp;
</span>-di<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp; </span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
TP<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
</span>\&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -fn
EXT -y1 100.00ppm -yn<span style="mso-spacerun: yes">&nbsp;&nbsp;&nbsp;
</span>138.00ppm -sw \&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; font-family:Courier;mso-bidi-font-family:Courier;color:#7D4700">| nmrPipe -out
hsqc.ft2 -verb 2 -ov</span><span style="font-size:8.0pt;mso-bidi-font-size: 11.0pt;mso-bidi-font-family:Courier">&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="font-size:8.0pt;mso-bidi-font-size:11.0pt; mso-bidi-font-family:Courier">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
 
<span style="mso-bidi-font-size:11.0pt;mso-bidi-font-family: Courier">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
 
'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">3.0<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>RDC Assignment and Calculation&lt;o:p&gt;&lt;/o:p&gt;</span>'''
 
'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>'''
 
'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">3.1<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>General Concepts</span>'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial"> &lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">Transfer all assignments (or cross peak labels) from a reference isotropic
HSQC spectrum to TROSY and HSQC spectra obtained under both isotropic and
aligned conditions. <span style="mso-spacerun: yes">&nbsp;</span>The positions of<span style="mso-spacerun: yes">&nbsp; </span>cross-peaks in the aligned HSQC spectrum usually change little with addition of an alignment medium; if they do change, there is often a constant offset for all peaks due to changes in the lock resonance.<span style="mso-spacerun: yes">&nbsp; </span>Hence, transferring assignments is usually straightforward. Corresponding peaks for TROSY assignments are diagonally displaced (down and to the right) by 38-58 Hz.<span style="mso-spacerun: yes">&nbsp; </span>The isotropic coupling (J) is calculated from the difference in Hz<span style="mso-spacerun: yes">&nbsp;
</span>between the <sup>15</sup>N dimension of isotropic HSQC and TROSY spectra (the offset is one half (J), so the final value needs to be multiply by 2).<span style="mso-spacerun: yes">&nbsp; </span>Similarly, the aligned coupling (J+D) is calculated from differences in peak positions for the aligned HSQC and TROSY spectra.<span style="mso-spacerun: yes">&nbsp; </span>The Hz values of the corresponding peaks can be subtracted using a spreadsheet or script to calculate the RDC values.&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">The peak picking and assignment process can be done using any graphical
NMR assignment program of user’s choice.<span style="mso-spacerun: yes">&nbsp;
</span>NMRViewJ has a special feature that would match and transfer assignment from one spectrum to another.<span style="mso-spacerun: yes">&nbsp;&nbsp;
</span>For this tutorial, we have included a Tcl/TK script for obtaining RDCs using NMRDraw and NMRWish; the process is broken down into five steps.&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
 
'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">3.2<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>RDC Scripts&lt;o:p&gt;&lt;/o:p&gt;</span>'''
 
'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>'''
 
'''<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">3.2.1<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Step 1: step1_hsqc_ass.tcl&lt;o:p&gt;&lt;/o:p&gt;</span>'''
 
'''&lt;o:p&gt;&nbsp;&lt;/o:p&gt;'''
 
'''Function:''' <span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family:Arial">Automatically pick
peaks in NMRDraw for the isotropic HSQC spectrum and save the peak list as
hsqc.tab to be use as reference peak list.</span> <span style="mso-spacerun: yes">&nbsp;</span>This script cuts out vast amounts of<span style="mso-spacerun: yes">&nbsp; </span>non-essential data that are stored in the hsqc.tab file after initial peak picking, leaving only the labels and, HN and N ppm values.'''<span style="mso-bidi-font-size: 14.0pt;mso-bidi-font-family:Arial">&lt;o:p&gt;&lt;/o:p&gt;</span>'''
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Execution:''' Executing this script does not require code modification. The only requirement is that you add the output and input script names in the execution line of a LINUX operating system.&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Command: &lt;o:p&gt;&lt;/o:p&gt;'''
 
./step1_hsqc_ass.tcl hsqc_ass.tab hsqc.tab&lt;o:p&gt;&lt;/o:p&gt;
 
Where “hsqc_ass.tab” is the output file and “hsqc.tab” is the input file to this script.&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
Download: step1_hsqc_ass.tcl <span style="mso-spacerun: yes">&nbsp;</span><span style="font-size:11.0pt;mso-bidi-font-size:14.0pt; mso-bidi-font-family:Arial">&lt;o:p&gt;&lt;/o:p&gt;</span>  
 
<span style="font-size:11.0pt;mso-bidi-font-size:14.0pt; mso-bidi-font-family:Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>  
 
'''3.2.2<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Step 2: i-hsqctrosy_step2.tcl&lt;o:p&gt;&lt;/o:p&gt;'''
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Function:''' This script takes the HSQC and TROSY peak tables that you created while peak picking and combines them into a single table. It also transfers the assignment from the isotropic HSQC to HSQC aligned, TROSY isotropic and TROSY aligned spectra. On execution, it opens “NMRWish” which can further help in transferring the assignment peak by peak from HSQC to HSQC-TROSY. All the blue peaks in the NMRWish display belong to the HSQC spectrum and the white peaks belong to the TROSYspectrum. The peak under consideration is on the crosswire of the NMRWish window.<span style="mso-spacerun: yes">&nbsp; </span>A navigation bar will also appear along with the NMRWish display to allow user to visually inspect the assignment of each peak and make corrections if necessary.&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Execution:''' Executing this script requires editing the script file itself and changing a few values. Make sure that '''assName '''(from step 1)''', tabName1'''(peak table for HSQC)''', tabName2''' (peak table for TROSY)''', and outName '''(output file name) have the correct input table names. You must also enter the height (contour level) at which you picked the peaks for hsqc in the '''hi1 '''field and the height at which you picked peaks for trosy in the '''hi2''' field. Also make sure that the correct names for the .ft2 spectrum are specified in the script. &lt;o:p&gt;&lt;/o:p&gt;
 
<span style="mso-spacerun: yes">&nbsp;</span>&lt;o:p&gt;&lt;/o:p&gt;
 
'''Command: &lt;o:p&gt;&lt;/o:p&gt;'''
 
./i-hsqctrosy_step2.tcl&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
Download: i-hsqctrosy_step2.tcl&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''3.2.3<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Step 3: i-hsqctrosy_step3.tcl'''&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Function:''' This script calculates the coupling values for the protein.&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Execution:''' Execution of this script requires that you open up the script and change the '''inName''' to the output of step 2. It also requires that you once again input the '''hi1''' and '''hi2''' values just like you did in step 2.&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Command:&lt;o:p&gt;&lt;/o:p&gt;'''
 
./i-hsqctrosy_step3.tcl&lt;o:p&gt;&lt;/o:p&gt;
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">Download: i-hsqctrosy_step3.tcl&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
 
'''3.2.4<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Step 4: i-hsqctrosy_step4.tcl'''&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Function:''' This script cleans up the output from the previous script (step 3) and outputs the necessary information (chemical shift of H and N in ppm and coupling values) into a new text file. &lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Execution:''' Execution of this script requires that you change a single variable in the script. Make sure that you set the '''outName''' variable to the proper media type (gel, peg, iso, etc…). &lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Command:&lt;o:p&gt;&lt;/o:p&gt;'''
 
./i-hsqctrosy_step4.tcl&lt;o:p&gt;&lt;/o:p&gt;
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">Download: i-hsqctrosy_step4.tcl&lt;o:p&gt;&lt;/o:p&gt;</span>
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">&lt;o:p&gt;&nbsp;&lt;/o:p&gt;</span>
 
'''3.2.5<span style="mso-tab-count:1">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Step 5: i-hsqctrosy_step5.tcl'''&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Function:''' This script calculates the rdc values for the protein.&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
'''Execution:''' We recommend that a new folder be created in the protein directory titled '''rdc_name?''' to organize and retain all of the files required to run this script. Copy or Move all of the required files into the '''rdc '''folder; these include: '''hsqc_ass.tab''' file and all of the '''*_j.tab''' files generated from isotropic data and each medium. Make sure that the '''assName, isoName, '''and '''othName '''We recommend the '''outName''' field to have the naming convention '''protein.name_rdc_media.tab''' (ex. rnr77_rdc_gel).&lt;o:p&gt;&lt;/o:p&gt;
 
'''<span style="mso-spacerun: yes">&nbsp;</span>&lt;o:p&gt;&lt;/o:p&gt;'''
 
'''Command:&lt;o:p&gt;&lt;/o:p&gt;'''
 
./i-hsqctrosy_step5.tcl&lt;o:p&gt;&lt;/o:p&gt;
 
&lt;o:p&gt;&nbsp;&lt;/o:p&gt;
 
<span style="mso-bidi-font-size:14.0pt;mso-bidi-font-family: Arial">Download: i-hsqctrosy_step5.tcl&lt;o:p&gt;&lt;/o:p&gt;</span> <!--EndFragment-->

Latest revision as of 19:35, 17 August 2011

NMR Data Processing >  2D HSQC-TROSY for RDCs Measurement

Theory/Background

Residual Dipolar Couplings (RDCs) originate from the anisotropic component of the dipolar interaction, which is dependent on the angle between an internuclear vector and the magnetic field. When a molecule samples orientations uniformly, as it does in normal solution NMR, RDCs average to zero and are not observable. However, if a molecule is dissolved in a dilute liquid crystalline medium, or other orienting medium, it becomes partially aligned, and the dipolar couplings are not completely averaged to zero.  This leads to a small contribution to the splitting of NMR signals. The angular dependence of these contributions can provide valuable structural information.  RDCs can be used to validate protein structures, to refine structures to improve quality, and to provide constraints as a part of an initial structure determination.  RDCs can also be used to determine the rotational symmetry axis of oligomers and this information can be use to properly orient the positions of subunits within an oligomeric complex. 


Recent Reviews

Prestegard, Bougault& Kishore, Chemical Reviews, 104, 3519-3540 (2004) 

Lipsitz& Tjandra, Ann. Rev. Biophys. Biomol. Struct., 33, 387-413 (2004) 

Additional information can also be found at the following web site.

http://tesla.ccrc.uga.edu/courses/bionmr/lectures

Brief Experimental Description

HSQC-TROSY pairs provide one way of measuring RDCs.  The procedure is applicable to proteins of small and moderate size (10-30 kDa).  The acquisition time in the indirect t1 should be at least T2 of nitrogen.  The measurements are easily visualized as displacements of peaks in 2D frequency domain spectra.  The HSQC member of the pair gives a single peak at the center of the four peaks expected in a fully coupled HSQC experiment. The TROSY member gives only the sharpest component. The offset in either H or N direction is one half (J).  Procedures for collecting this pair of experiments are described elsewhere in this Twiki (Experimental setup for HSQC-TROSY).

NMR Data Processing

The data processing described below is specific for data collected on Varian spectrometers having the BioPack (version: 2009-09-03) sequence (gNhsqc.c) installed.  Or on Bruker spectrometers having the …….. It is applied to both isotropic and aligned HSQC-TROSY pairs.

Software Information

NMRPipe (download information and user manual)
http://spin.niddk.nih.gov/NMRPipe/

Brief descriptions of specific functions are accessible via nmrPipe –

‘nmrPipe –help’ will list most functions.
‘nmrPipe –fn GM –help’ will give description of the GM function.


Supported Platforms

Linux (RedHat Linux/Fedora)
Mac OS X (10.3.4 and up)
SGI Irix  (6.2 and up)
Sparc Solaris (2 and up)
Windows XP Pro with Microsoft Services for UNIX (SFU 3.5)
Varian to NMRPipe Data Conversion
#!/bin/csh
##conversion of varian 2d fids into nmrPipe format for a TROSY/hsqc data set
##with array= ‘(TROSY, dm2),phase’
var2pipe -in ./fid \
   -xN           2048    -yN      1024   \
   -xT           1024    -yT      256   \
   -xMODE     Complex    -yMODE   Rance-Kay  \
   -xSW       10000.00   -ySW      2800.00  \
   -xOBS      599.520    -yOBS      60.756   \
   -xCAR         4.773    -yCAR    118.262   \
   -xLAB           HN    -yLAB         N   \
   -ndim            2    -aq2D    States   \
   | nmrPipe -fn QMIX -ic 4 -oc 2 -cList 1 0 0 1 0 0 0 0   \
   -out trosy.fid -verb -ov

var2pipe -in ./fid \
  -xN           2048    -yN      1024   \
   -xT           1024    -yT      256   \
   -xMODE     Complex    -yMODE   Rance-Kay  \
   -xSW       10000.00   -ySW     2800.00  \
   -xOBS      599.520    -yOBS      60.756   \
   -xCAR         4.773    -yCAR    118.262   \
   -xLAB           HN    -yLAB         N   \
   -ndim            2    -aq2D    States   \
   | nmrPipe -fn QMIX -ic 4 -oc 2 -cList 0 0 0 0 1 0 0 1   \
   -out hsqc.fid -verb –ov 

Cautions: The usual Varian to NMRPipe macros will read in experimental parameters, but will not be able to perform data conversion since both the HSQC and TROSY spectra are collected in a interleaved manner and are stored in a single FID data file.  The above script includes the QMIX function to separate fids.

Bruker to NMRPipe Data Conversion

The conversion of data collected on most Bruker spectrometers can be accomplished with the following script.

(insert script)

Fourier Transformation of Converted HSQC and TROSY spectra

The following script contains some data specific parameters including phasing parameters and convolution functions that will need to be optimized for each specific set.

#!/bin/csh
#This script is to process one of the TROSY/
| nmrPipe  -in hsqc.fid              \
| nmrPipe -fn SOL   \
| nmrPipe  -fn SP -off 0.5 -end 0.99 -pow 2 -c 0.5  \
| nmrPipe  -fn ZF -auto                    \
| nmrPipe  -fn FT                     \
| nmrPipe -fn PS -p0 -61.6 -p1 -4.6  -di    \
| nmrPipe  -fn EXT -x1 6.00ppm -xn 11.50ppm -verb 2 -sw  \
| nmrPipe -fn TP                      \
| nmrPipe -fn LP -x1 1 -xn 256 -pred 256 -after -verb   \
| nmrPipe -fn SP -off 0.5 -end 0.99 -pow 2 -c 1.0  \
| nmrPipe -fn ZF -auto           \
| nmrPipe -fn FT                   \
| nmrPipe -fn PS -p0 -91.8 -p1 7.0   -di     \
| nmrPipe -fn TP        \
| nmrPipe -fn EXT -y1 100.00ppm -yn    138.00ppm -sw \
| nmrPipe -out hsqc.ft2 -verb 2 -ov
RDC Assignment and Calculation

General Concepts 

Transfer all assignments (or cross peak labels) from a reference isotropic HSQC spectrum to TROSY and HSQC spectra obtained under both isotropic and aligned conditions.  The positions of  cross-peaks in the aligned HSQC spectrum usually change little with addition of an alignment medium; if they do change, there is often a constant offset for all peaks due to changes in the lock resonance.  Hence, transferring assignments is usually straightforward. Corresponding peaks for TROSY assignments are diagonally displaced (down and to the right) by 38-58 Hz.  The isotropic coupling (J) is calculated from the difference in Hz  between the 15N dimension of isotropic HSQC and TROSY spectra (the offset is one half (J), so the final value needs to be multiply by 2).  Similarly, the aligned coupling (J+D) is calculated from differences in peak positions for the aligned HSQC and TROSY spectra.  The Hz values of the corresponding peaks can be subtracted using a spreadsheet or script to calculate the RDC values.

The peak picking and assignment process can be done using any graphical NMR assignment program of user’s choice.  NMRViewJ has a special feature that would match and transfer assignment from one spectrum to another.   For this tutorial, we have included a Tcl/TK script for obtaining RDCs using NMRDraw and NMRWish; the process is broken down into five steps.


RDC Assignment and Calculation Scripts

Step 1: step1_hsqc_ass.tcl

Function: Automatically pick peaks in NMRDraw for the isotropic HSQC spectrum and save the peak list as hsqc.tab to be use as reference peak list.  This script cuts out vast amounts of  non-essential data that are stored in the hsqc.tab file after initial peak picking, leaving only the labels and, HN and N ppm values.


Execution: Executing this script does not require code modification. The only requirement is that you add the output and input script names in the execution line of a LINUX operating system.


Command: 

./step1_hsqc_ass.tcl hsqc_ass.tab hsqc.tab

Where “hsqc_ass.tab” is the output file and “hsqc.tab” is the input file to this script.


Download: Step1_hsqc_ass.tcl  


Step 2: i-hsqctrosy_step2.tcl

Function: This script takes the HSQC and TROSY peak tables that you created while peak picking and combines them into a single table. It also transfers the assignment from the isotropic HSQC to HSQC aligned, TROSY isotropic and TROSY aligned spectra. On execution, it opens “NMRWish” which can further help in transferring the assignment peak by peak from HSQC to HSQC-TROSY. All the blue peaks in the NMRWish display belong to the HSQC spectrum and the white peaks belong to the TROSYspectrum. The peak under consideration is on the crosswire of the NMRWish window.  A navigation bar will also appear along with the NMRWish display to allow user to visually inspect the assignment of each peak and make corrections if necessary.

Step2-spectrum.png Step2-navi.png

Execution: Executing this script requires editing the script file itself and changing a few values. Make sure that assName (from step 1), tabName1(peak table for HSQC), tabName2 (peak table for TROSY), and outName (output file name) have the correct input table names. You must also enter the height (contour level) at which you picked the peaks for hsqc in the hi1 field and the height at which you picked peaks for trosy in the hi2 field. Also make sure that the correct names for the .ft2 spectrum are specified in the script. 

 

Command: 

./i-hsqctrosy_step2.tcl


Download: I-hsqctrosy_step2.tcl


Step 3: i-hsqctrosy_step3.tcl

Function: This script calculates the coupling values for the protein.


Execution: Execution of this script requires that you open up the script and change the inName to the output of step 2. It also requires that you once again input the hi1 and hi2 values just like you did in step 2.


Command:

./i-hsqctrosy_step3.tcl


Download: I-hsqctrosy_step3.tcl


Step 4: i-hsqctrosy_step4.tcl

Function: This script cleans up the output from the previous script (step 3) and outputs the necessary information (chemical shift of H and N in ppm and coupling values) into a new text file. 


Execution: Execution of this script requires that you change a single variable in the script. Make sure that you set the outName variable to the proper media type (gel, peg, iso, etc…). 


Command:

./i-hsqctrosy_step4.tcl


Download: I-hsqctrosy_step4.tcl


Step 5: i-hsqctrosy_step5.tcl

Function: This script calculates the rdc values for the protein.


Execution: We recommend that a new folder be created in the protein directory titled rdc_name? to organize and retain all of the files required to run this script. Copy or Move all of the required files into the rdc folder; these include: hsqc_ass.tab file and all of the *_j.tab files generated from isotropic data and each medium. Make sure that the assName, isoName, and othName We recommend the outName field to have the naming convention protein.name_rdc_media.tab (ex. rnr77_rdc_gel).

 

Command:

./i-hsqctrosy_step5.tcl


Download: I-hsqctrosy_step5.tcl


Updated by Hsiau-Wei Lee, 2011