HA and HB Assignment with GFT in XEASY: Difference between revisions

From NESG Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(4 intermediate revisions by 2 users not shown)
Line 5: Line 5:
=== '''Analysis of the(4,3)D GFT HABCAB(CO)NHN'''  ===
=== '''Analysis of the(4,3)D GFT HABCAB(CO)NHN'''  ===


The sequence-specific 15N /1HN /13CAB resonance assignments obtained as described in backbone assignment are first used to also obtain 1HA and 1HB shift assignments by analyzing (4,3)D HABCAB(CO)NHN. Then, shifts of more peripheral aliphatic spins are obtained from (4,3)D HCCH.  
The sequence-specific 15N / 1HN / 13CAB resonance assignments obtained as described in backbone assignment are first used to also obtain 1HA and 1HB shift assignments by analyzing (4,3)D HABCAB(CO)NHN. Then, shifts of more peripheral aliphatic spins are obtained from (4,3)D HCCH.  


#Go to the directory /analisys/xeasy/habcab and copy the final-clean.seq and final-clean .prot files from the backbone directory. In UBNMR, run <tt>makeHabcabPeaks</tt> to generate an extended AtomList <tt>habcabconhI1.prot</tt>, which contains linear combinations of 13CAB and 1HAB shifts, and to generate a starting HABCAB peak list <tt>habcabconhI1.peaks</tt> to guide manual peak identification in (4,3)D HABCAB(CO)NHN. Peaks are colored according to type of sub-spectrum and type of proton involved in the GFT dimension.
#Go to the directory <tt>/analisys/xeasy/habcab</tt> and copy the <tt>final-clean.seq</tt> and <tt>final-clean.prot</tt> files from the backbone directory  
#Run <tt>makeHabcabPeaks</tt> in [[UBNMR|UBNMR]] to generate an extended [[XEASY Atom List|AtomList]] <tt>habcabconhI1.prot</tt> (containing linear combinations of 13CAB and 1HAB shifts) and the initial HABCAB [[XEASY Peak List|PeakList]] <tt>habcabconhI1.peaks</tt> to guide manual peak identification in (4,3)D HABCAB(CO)NHN. Peaks are colored according to type of sub-spectrum and type of proton involved in the GFT dimension.  
#In XEASY
##use <tt>ns</tt> to load '''subspectrum I''', <tt>HABCABcoNH1</tt>, from <tt>/protName/analysis/xeasy/habcab/</tt>
##<tt>cp</tt> to display the spectrum as a contour plot
##<tt>ls</tt> to load the [[XEASY Sequence List|SequencList]] (<tt>protein.seq</tt>)
##<tt>lc</tt> to load [[XEASY Atom List|AtomList]] (<tt>habcabconhI1.prot</tt>)
##<tt>lp</tt> to load the starting [[XEASY Peak List|PeakList]] (<tt>habcabconhI1.peaks</tt>)
##<tt>se</tt> to create the strips for '''subspectrum I'''
##<tt>ns</tt> to load '''subspectrum II''', <tt>HABCABcoNH2</tt>
##<tt>cp</tt> to display the spectrum as a contour plot
##<tt>se</tt> to create the strips for '''subspectrum II''' and <tt>gs</tt> to display the strips (see Figure 1, below)
##<tt>pw</tt> and click the box for '''assignment y''' to display the assignments in the GFT HabCab dimension
##<tt>mr</tt> to accurately position the peaks along w1(13C) and <tt>bs</tt> / <tt>fs</tt> to move through the strips
###when moving peaks, be sure to move the proper peak for each of the subspectra&nbsp;('''NOTE''':&nbsp;if the spectra were processed using Buffalo's PERL scripts, then '''subspectrum I'''&nbsp;should use the "plus" peaks -- CA+HA, CB+HB2, CB+HB3, etc.)
###if the HA1/2 or HB2/3 peaks are degenerate, move both of the peak markers to the same position and try to resolve them in the NOESY<br>
###for non-degenerate HB2/3 (or HA1/2) peaks, place the CB+HB3 and CB-HB3 peaks on the "outside" and the CB+HB2 and CB-HB2 peaks on the "inside" of the peak pairs
###if you are unable to correctly identify the peak pairs&nbsp;(e.g., CA+HA&nbsp;and CA-HA), use the 13C single-quantum chemical shifts from the backbone experiments to identify the center of the peak pair
#In UBNMR, run <tt>updateHabAtom</tt> to update single-quantum 1HAB and 13CAB shifts as <tt>habcabconhO2.prot</tt>. A least-squares fit provides the single-quantum HA/HB chemical shifts (consistency of peaks representing different linear combinations of shifts is checked by UBNMR in oder to identify inadvertently 'mis-picked' peaks; a warning is then provided).
#In XEASY, double-check the peak positions for residues that gaving large error. Repeat previous step and this step till no improvement can be made.
#Now it is ready to move to [[XEASY Side Chain Assignments|Side Chain Assignments]]
#'''Figure 1: Example of (4,3)D GFT HABCAB(CO)NHN analysis before (left) and after (right) peak position adjustment by <tt>mr</tt>.''' <br> [[Image:XEASY hab5.jpg]] [[Image:XEASY hab4.jpg]] <br>


<br>  
<br>  
Line 13: Line 34:
<br>  
<br>  


#In XEASY, use <tt>ns</tt>, <tt>ls</tt>, <tt>lc</tt> and <tt>lp</tt> to load the sub-spectra of (4,3)D HABCAB(CO)NHN and the starting HABCAB-peak list; use <tt>se</tt>, <tt>gs</tt> to display [w1(13C,1H),w3(1HN)]-strips; use <tt>mr</tt> to accurately position peaks; use <tt>aa</tt>, <tt>ac</tt>, <tt>wp</tt> and <tt>wc</tt> to update the peak list as <tt>habcabconhO1.peaks</tt> and AtomList as <tt>habcabconhO1.prot</tt>. See the analysis procedure below.
*[[Media:XEASY_makeHabcabPeaks.txt|makeHabcabPeaks]]: UBNMR macro
##In XEASY, use <tt>ns</tt> to load subspectrum I, <tt>HABCABcoNH1</tt>, from <tt>/protName/analysis/xeasy/data/</tt>; <tt>ls</tt> to load the sequence file, e.g. <tt>protein.seq</tt>; <tt>lc</tt> to load AtomList from <tt>habcabconhI1.prot</tt>; <tt>lp</tt> to load the starting peaklist <tt>habcabconhI1.peaks</tt>; use <tt>ns</tt> to load sub-spectrum II, <tt>HABCABcoNH2</tt>.
##In XEASY, use <tt>cp</tt> to display the spectrum as a countour plot; use <tt>fp</tt> or <tt>bp</tt> to step through the planes until you see representative peaks for setting the coutour levels; use <tt>pm</tt> to display [w1(13C/1H),w2(15N)]-planes; use <tt>mr</tt> to accurately position the peaks along w2(15N); use <tt>pm</tt> to display [w2(15N),w3(1HN)]-planes.
##In XEASY, use <tt>se</tt> and <tt>gs</tt> to create strips from entries of loaded PeakList and display, for example, ten at once; navigate through the strips by using <tt>fs</tt> and <tt>bs</tt>; use <tt>pw</tt> to set the '''assignment y''' in the peak window so that only the assignment in the GFT dimension is displayed. One should see strips from both sub-spectra assigned to the same residue side-by-side.
##In XEASY, use <tt>gs</tt>, <tt>sf</tt>, <tt>fs</tt> and <tt>bs</tt> to go through strips and use <tt>mr</tt> to move the peaks to the actual peak positions. Make sure that the peaks are moved in the right sub-spectrum. If two HB-proton shifts are apparently degenerate, move the corresponding peak to the same position (note that in the better resolved NOESY, the two shifts may turn out to be non-degenerate). If the assignment is ambiguous, load 'central peak spectrum' to identify peak pairs in the sub-spectra. <br> '''Figure 1: Example of (4,3)D GFT HABCAB(CO)NHN analysis before (left) and after (right) peak position adjustment by <tt>mr</tt>.''' <br> [[File:XEASY_hab5.jpg]] [[File:XEASY_hab4.jpg]] <br>
#In UBNMR, run <tt>updateHabAtom</tt> to update single-quantum 1HAB and 13CAB shifts as <tt>habcabconhO2.prot</tt>. A least-squares fit provides the single-quantum HA/HB chemical shifts (consistency of peaks representing different linear combinations of shifts is checked by UBNMR in oder to identify inadvertently 'mis-picked' peaks; a warning is then provided).
#In XEASY, double-check the peak positions for residues that gaving large error. Repeat previous step and this step till no improvement can be made.
#Now it is ready to move to [[XEASY Side Chain Assignments|Side Chain Assignments]]
 
=== '''[[NESG:Assignment validation|Assignment validation]]'''  ===
 
Please refer to the assignment validation in [[NESG:BackboneAssignment|BackboneAssignment]]
 
<br> %COMMENT%
 
<br> -- Main.GaohuaLiu - 16 Feb 2007
 
*[[NESG:%ATTACHURL%/makeHabcabPeaks|makeHabcabPeaks]]: UBNMR macro


*[[NESG:%ATTACHURL%/updateHabAtom|updateHabAtom]]: UBNMR macro
*[[Media:XEASY_updateHabAtom.txt|updateHabAtom]]: UBNMR macro

Latest revision as of 18:07, 1 December 2009

HA And HB Assignments with XEASY/UBNMR

HA and HB assignments provide the bridge from Backbone Assignment to the complete Side Chain Assignments.

Analysis of the(4,3)D GFT HABCAB(CO)NHN

The sequence-specific 15N / 1HN / 13CAB resonance assignments obtained as described in backbone assignment are first used to also obtain 1HA and 1HB shift assignments by analyzing (4,3)D HABCAB(CO)NHN. Then, shifts of more peripheral aliphatic spins are obtained from (4,3)D HCCH.

  1. Go to the directory /analisys/xeasy/habcab and copy the final-clean.seq and final-clean.prot files from the backbone directory
  2. Run makeHabcabPeaks in UBNMR to generate an extended AtomList habcabconhI1.prot (containing linear combinations of 13CAB and 1HAB shifts) and the initial HABCAB PeakList habcabconhI1.peaks to guide manual peak identification in (4,3)D HABCAB(CO)NHN. Peaks are colored according to type of sub-spectrum and type of proton involved in the GFT dimension.
  3. In XEASY
    1. use ns to load subspectrum I, HABCABcoNH1, from /protName/analysis/xeasy/habcab/
    2. cp to display the spectrum as a contour plot
    3. ls to load the SequencList (protein.seq)
    4. lc to load AtomList (habcabconhI1.prot)
    5. lp to load the starting PeakList (habcabconhI1.peaks)
    6. se to create the strips for subspectrum I
    7. ns to load subspectrum II, HABCABcoNH2
    8. cp to display the spectrum as a contour plot
    9. se to create the strips for subspectrum II and gs to display the strips (see Figure 1, below)
    10. pw and click the box for assignment y to display the assignments in the GFT HabCab dimension
    11. mr to accurately position the peaks along w1(13C) and bs / fs to move through the strips
      1. when moving peaks, be sure to move the proper peak for each of the subspectra (NOTE: if the spectra were processed using Buffalo's PERL scripts, then subspectrum I should use the "plus" peaks -- CA+HA, CB+HB2, CB+HB3, etc.)
      2. if the HA1/2 or HB2/3 peaks are degenerate, move both of the peak markers to the same position and try to resolve them in the NOESY
      3. for non-degenerate HB2/3 (or HA1/2) peaks, place the CB+HB3 and CB-HB3 peaks on the "outside" and the CB+HB2 and CB-HB2 peaks on the "inside" of the peak pairs
      4. if you are unable to correctly identify the peak pairs (e.g., CA+HA and CA-HA), use the 13C single-quantum chemical shifts from the backbone experiments to identify the center of the peak pair
  4. In UBNMR, run updateHabAtom to update single-quantum 1HAB and 13CAB shifts as habcabconhO2.prot. A least-squares fit provides the single-quantum HA/HB chemical shifts (consistency of peaks representing different linear combinations of shifts is checked by UBNMR in oder to identify inadvertently 'mis-picked' peaks; a warning is then provided).
  5. In XEASY, double-check the peak positions for residues that gaving large error. Repeat previous step and this step till no improvement can be made.
  6. Now it is ready to move to Side Chain Assignments
  7. Figure 1: Example of (4,3)D GFT HABCAB(CO)NHN analysis before (left) and after (right) peak position adjustment by mr.
    XEASY hab5.jpg XEASY hab4.jpg