HSQCTROSY RDC Measurement: Difference between revisions

From NESG Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(18 intermediate revisions by the same user not shown)
Line 1: Line 1:
'''NMR Data Processing >  2D HSQC-TROSY for RDCs Measurement'''
'''NMR Data Processing &gt;&nbsp; 2D HSQC-TROSY for RDCs Measurement''' <br>  
<br>
 
===== '''Theory/Background''' =====
 
<br>


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;
===== '''Theory/Background'''  =====


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


===== '''Recent Reviews''' =====
<br>


<br>
===== '''Recent Reviews'''  =====


Prestegard, Bougault&amp; Kishore, ''Chemical Reviews'', '''104''', 3519-3540 (2004)&nbsp;
Prestegard, Bougault&amp; Kishore, ''Chemical Reviews'', '''104''', 3519-3540 (2004)&nbsp;  


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


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


<span style="text-decoration: underline">[http://tesla.ccrc.uga.edu/courses/bionmr/lectures http://tesla.ccrc.uga.edu/courses/bionmr/lectures]</span>
<span style="text-decoration: underline">[http://tesla.ccrc.uga.edu/courses/bionmr/lectures http://tesla.ccrc.uga.edu/courses/bionmr/lectures]</span>
<br>
===== '''Brief Experimental Description''' =====
<br>
<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). &nbsp;


<br>
===== '''Brief Experimental Description'''  =====


===== '''NMR Data Processing''' =====
<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).


<br>
===== '''NMR Data Processing'''  =====


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.
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.
<br>
===== '''Software Information''' =====
<br>
NMRPipe (download information and user manual)
<span style="text-decoration: underline">[http://%22 http://spin.niddk.nih.gov/NMRPipe/]</span>
<span class="Apple-tab-span" style="white-space:pre"> </span>
Brief descriptions of specific functions are accessible via nmrPipe –
<span style="font: 12.0px Times New Roman"><span class="Apple-tab-span" style="white-space:pre"> </span></span>‘nmrPipe –help’ will list most functions.
<span class="Apple-tab-span" style="white-space:pre"> </span>‘nmrPipe –fn GM –help’ will give description of the GM function.
<br>
Supported Platforms
Linux (RedHat Linux/Fedora)&nbsp;
Mac OS X (10.3.4 and up)
SGI Irix&nbsp; (6.2 and up)
&nbsp;<span class="Apple-tab-span" style="white-space:pre"> </span>Sparc Solaris (2 and up)
Windows XP Pro with Microsoft Services for UNIX (SFU 3.5).
<br>
<br>


===== '''Varian to NMRPipe Data Conversion''' =====
===== Software Information  =====


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


#!/bin/csh
:‘nmrPipe –help’ will list most functions.
:‘nmrPipe –fn GM –help’ will give description of the GM function.


# conversion of varian 2d fids into nmrPipe format for a TROSY/hsqc data set
<br> Supported Platforms


# with array= ‘(TROSY, dm2),phase’
:Linux (RedHat Linux/Fedora)
:Mac OS X (10.3.4 and up)
:SGI Irix &nbsp;(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'''  =====
<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 \
var2pipe -in ./fid \
 
  -xN           2048   -yN     1024   \
&nbsp;&nbsp; -xN &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 2048&nbsp; &nbsp; -yN&nbsp; &nbsp; &nbsp; 1024 &nbsp; \
  -xT           1024   -yT     256   \
 
  -xMODE     Complex   -yMODE   Rance-Kay \
&nbsp;&nbsp; -xT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 1024&nbsp; &nbsp; -yT&nbsp; &nbsp; &nbsp; 256 &nbsp; \
  -xSW       10000.00   -ySW     2800.00 \
 
  -xOBS     599.520   -yOBS     60.756   \
&nbsp;&nbsp; -xMODE &nbsp; &nbsp; Complex&nbsp; &nbsp; -yMODE &nbsp; Rance-Kay&nbsp; \
  -xCAR         4.773   -yCAR   118.262   \
 
  -xLAB           HN   -yLAB         N   \
&nbsp;&nbsp; -xSW &nbsp; &nbsp; &nbsp; 10000.00 &nbsp; -ySW&nbsp; &nbsp; &nbsp; 2800.00&nbsp; \
  -ndim           2   -aq2D   States   \
 
  | nmrPipe -fn QMIX -ic 4 -oc 2 -cList 1 0 0 1 0 0 0 0   \
&nbsp;&nbsp; -xOBS&nbsp; &nbsp; &nbsp; 599.520&nbsp; &nbsp; -yOBS&nbsp; &nbsp; &nbsp; 60.756 &nbsp; \
  -out trosy.fid -verb -ov
 
&nbsp;&nbsp; -xCAR &nbsp; &nbsp; &nbsp; &nbsp; 4.773&nbsp; &nbsp; -yCAR&nbsp; &nbsp; 118.262 &nbsp; \
 
&nbsp;&nbsp; -xLAB &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; HN&nbsp; &nbsp; -yLAB &nbsp; &nbsp; &nbsp; &nbsp; N &nbsp; \
 
&nbsp;&nbsp; -ndim&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; -aq2D&nbsp; &nbsp; States &nbsp; \
 
&nbsp;&nbsp; | nmrPipe -fn QMIX -ic 4 -oc 2 -cList 1 0 0 1 0 0 0 0 &nbsp; \
 
&nbsp;&nbsp; -out trosy.fid -verb -ov
 
<br>


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


<br>
===== '''Bruker to NMRPipe Data Conversion''' =====
 
===== '''Bruker to NMRPipe Data Conversion''' =====
 
<br>


<span style="font: 12.0px Times New Roman"><span class="Apple-tab-span" style="white-space:pre"> </span></span>The conversion of data collected on most Bruker spectrometers can be accomplished with the following script.&nbsp; The same precautions pertain?
The conversion of data collected on most Bruker spectrometers can be accomplished with the following script.  


(insert script)
(insert script)


<br>
===== '''Fourier Transformation of Converted HSQC and TROSY spectra''' =====
 
===== '''Fourier Transformation of Converted HSQC and TROSY spectra''' =====
<br>
 
<span style="font: 8.0px Courier">'''<span class="Apple-tab-span" style="white-space:pre"> </span>'''</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.
 
<br>
 
'''#!/bin/csh'''
 
# This script is to process one of the TROSY/
 
| nmrPipe&nbsp; -in hsqc.fid&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; \
 
| nmrPipe -fn SOL &nbsp; \
 
| nmrPipe&nbsp; -fn SP -off 0.5 -end 0.99 -pow 2 -c 0.5&nbsp; \
 
| nmrPipe&nbsp; -fn ZF -auto&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; \
 
| nmrPipe&nbsp; -fn FT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; \
 
| nmrPipe -fn PS -p0 -61.6 -p1 -4.6&nbsp; -di&nbsp; &nbsp; \
 
| nmrPipe&nbsp; -fn EXT -x1 6.00ppm -xn 11.50ppm -verb 2 -sw&nbsp; \
 
| nmrPipe -fn TP&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; \
 
| nmrPipe -fn LP -x1 1 -xn 256 -pred 256 -after -verb &nbsp; \
 
| nmrPipe -fn SP -off 0.5 -end 0.99 -pow 2 -c 1.0&nbsp; \
 
| nmrPipe -fn ZF -auto &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; \
 
| nmrPipe -fn FT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; \
 
| nmrPipe -fn PS -p0 -91.8 -p1 7.0 &nbsp; -di &nbsp; &nbsp; \
 
| nmrPipe -fn TP&nbsp; &nbsp; &nbsp; &nbsp; \
 
| nmrPipe -fn EXT -y1 100.00ppm -yn&nbsp; &nbsp; 138.00ppm -sw \


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
| nmrPipe -out hsqc.ft2 -verb 2 -ov
</pre>


<br>
===== '''RDC Assignment and Calculation'''  =====


<br>
'''General Concepts'''&nbsp;


===== '''RDC Assignment and Calculation''' =====
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.
<br>


'''General Concepts'''&nbsp;
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.


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


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.
===== '''RDC Assignment and Calculation Scripts'''  =====


<br>
'''Step 1: step1_hsqc_ass.tcl''' <br>  


===== '''RDC Assignment and Calculation Scripts''' =====
'''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.
<br>


'''Step 1: step1_hsqc_ass.tcl'''
<br>  
<br>


'''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.
'''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>


<br>
<br>  


'''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>
'''Command:&nbsp;'''  
<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.  


<br>
<br>  


'''Command:&nbsp;'''
Download: [[Media:Step1_hsqc_ass.tcl|Step1_hsqc_ass.tcl]] &nbsp;  


./step1_hsqc_ass.tcl hsqc_ass.tab hsqc.tab
<br>


Where “hsqc_ass.tab” is the output file and “hsqc.tab” is the input file to this script.
'''Step 2: i-hsqctrosy_step2.tcl'''


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


Download: step1_hsqc_ass.tcl &nbsp;
[[Image:Step2-spectrum.png|Step2-spectrum.png]] [[Image:Step2-navi.png|Step2-navi.png]] <br>


<br>
'''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;
 
'''Step 2: i-hsqctrosy_step2.tcl'''


<br>
&nbsp;


'''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.
'''Command:&nbsp;'''  
<pre>./i-hsqctrosy_step2.tcl</pre>
<br>


<br>
Download: [[Media:I-hsqctrosy_step2.tcl|I-hsqctrosy_step2.tcl]]


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


&nbsp;
'''Step 3: i-hsqctrosy_step3.tcl'''


'''Command:&nbsp;'''
'''Function:''' This script calculates the coupling values for the protein.


./i-hsqctrosy_step2.tcl
<br>


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


Download: i-hsqctrosy_step2.tcl
<br>


<br>
'''Command:'''
<pre>./i-hsqctrosy_step3.tcl</pre>
<br>  


'''Step 3: i-hsqctrosy_step3.tcl'''
Download: [[Media:I-hsqctrosy_step3.tcl|I-hsqctrosy_step3.tcl]]


<br>
<br>  


'''Function:''' This script calculates the coupling values for the protein.
'''Step 4: i-hsqctrosy_step4.tcl'''  


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


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


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


'''Command:'''
<br>


./i-hsqctrosy_step3.tcl
'''Command:'''
<pre>./i-hsqctrosy_step4.tcl</pre>
<br>


<br>
Download: [[Media:I-hsqctrosy_step4.tcl|I-hsqctrosy_step4.tcl]]


Download: i-hsqctrosy_step3.tcl
<br>


<br>
'''Step 5: i-hsqctrosy_step5.tcl'''


'''Step 4: i-hsqctrosy_step4.tcl'''
'''Function:''' This script calculates the rdc values for the protein.


<br>
<br>  


'''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;
'''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).  


<br>
'''&nbsp;'''


'''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;
'''Command:'''  
<pre>./i-hsqctrosy_step5.tcl</pre>
<br>


Download: [[Media:I-hsqctrosy_step5.tcl|I-hsqctrosy_step5.tcl]]
<br>
<br>
'''Command:'''
./i-hsqctrosy_step4.tcl
<br>
<br>
Download: i-hsqctrosy_step4.tcl
<br>
'''Step 5: i-hsqctrosy_step5.tcl'''
<br>
'''Function:''' This script calculates the rdc values for the protein.
<br>
<br>
 
Updated by Hsiau-Wei Lee, 2011
'''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).
 
'''&nbsp;'''
 
'''Command:'''
 
./i-hsqctrosy_step5.tcl
 
<br>
 
Download: i-hsqctrosy_step5.tcl
<div><font class="Apple-style-span" face="Cambria" size="3"><span class="Apple-style-span" style="font-size: 12px; line-height: normal;">
</span></font></div>

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