Side chain assignment with aliphatic (4,3)D HCCH-COSY in XEASY

From NESG Wiki
Jump to navigation Jump to search


Aliphatic Side-chain Resonance Assignment Using the Aliphatic HCCHCOSY Spectrum

The sequence specific 15N /1HN /13CAB resonance assignments obtained as described in Backbone Assignment are first used to also obtain 1HA and 1HB Assignment by analyzing (4,3)D HABCAB(CO)NHN. Then, shifts of more peripheral aliphatic spins are obtained from (4,3)D HCCH or other spectra. Two approaches by using (4,3)D HCCH are described below and the later one requires more experience. NOTE: Here atom A, B, G, D, E stand for atom alpha α, beta β, gama γ, delta δ, and epsilon ε. Dimensions of HCCH are w1, 13C; _w2, (13C+-1H), the GFT dimension;_ and w3, 1H, direct detect dimension.

Approach I, using UBNMR to update SQ chemical shift

  1. Go to analysis/xeasy/hcchcosy.In UBNMR, run makeCosyPeak to generate an extended AtomList as bbscgftI1.prot which contains all required linear combinations to assign more peripheral aliphatic spins; A complete starting HC-peak list as hcchcosyI1.peaks is generated for analysis of (4,3)D HCCH, taking averaged HG/CG-, HD/CD-, HE/Main.CE- shifts from the BMRB for thus far unassigned side-chain chemical shifts. Peaks predicted based on the BMRB shifts are colored according to being detected on =HG/CG=-,=HD/CD=-,=HE/Main.CE=-moieties.
  2. In XEASY, use ns to load HCCH sub-spectrum I (comprising peaks at w2(13C+1H)), sub-spectrum II (comprising peaks at w2(13C-1H)), and the central peak sub-spectrum; use ls, lc and lp to load the corresponding SequenceList protein.seq AtomList bbscgftI1.prot and starting PeakList hcchcosyI1.peaks; use dw to set the displayed color for different types of peaks; use sp to select the peaks for display.
  3. In XEASY, use cp to display the sub-spectra as countour plots; use es, se, rs, se, rs, se and gs to create strips for all three spectra (used rs here to change the spectra, followed by se to search and sort strips) and display strips in the order of residue sequence (Figure 1).
    Figure 1. Example of GFT (4,3)D HCCH analysis. Side by side [w2(13C;1H),w3(1HA)]-strip view of (4,3)D HCCH spectra in the order of central spectum and two basic spectra (sub-spectrum I and II).
    XEASY cos1.jpg

  4. In XEASY, for each residue start with [w2(13C;1H),w3(1HA)]-strip containing signals detected on 1HA (refered to in the following as HA strip). Provided that backbone and CB/HB chemical shifts were accurately obtained in previous steps, each HA strip in the basic spectra of (4,3)D HCCH should have peaks accurately positioned at w(13CA +- 1HA) and w(13CB +- 1HB) along w2 (the GFT dimension) as shown in Figure 1. In case peaks are not accurately positioned, use mr along X-axis (w3(HA)) or Y-axis (w2(13C;1H)) to move peaks to the correct position. After adjusting all peaks of all HA strips, use ac, wc and wp to save updated AtomList and Peaklist. If the expected CA shift definig the 13C-plane of the peaks is different from the actually observed shift along w2(13C), peak positions need to be adjusted the correct plane; use fp and bp to search for correct plane; use pm and gs to display [w2(13C;1H),w1(13CA)]-strips (shown in Figure 2); use mr to move the peaks along X-axis( w1(13C)) into the correct 13C-shift. If two B-proton shifts are apparently degenerate, move both corresponding peaks to the same position (note that in better resolved NOESY the two shifts may appear to be non-degenerate).
    Figure 2. Example of GFT (4,3)D HCCH analysis. Side by side [w2(13C;1H),w1(13C)]-strip view of (4,3)D HCCH spectra in the order of central spectum and two basic spectra (sub-spectrum I and II).
    XEASY cos2.jpg

  5. In UBNMR, run updatecosyGFT to update the SQ chemical shifts as hcchcosyI2.prot and re-generate a new peak list hcchcosyI2.peaks to proceed further with side-chain assignments.
  6. In XEASY, use lp to re-load the PeakList; es, se; and gs; to re-sort the strips of the three sub-spectra as described above; use sf to display the [w1(13C;1H),w3(1HB)]-strip (CB/HB strips). The peaks are colored according to the proton on which signals are detected, that is, different colors are used for peaks detected on A-, B-, G-, D-, and E- protons; use mr to accurately position peaks at w(13CG) and w(13CG+-1HG) (Figure 3); use ac, wc and wp to save updated AtomList hcchcosyO2.prot and PeakList hcchcosyO2.peaks.
    Figure 3 Example of GFT (4,3)D HCCH analysis. Side by side [w2(13C;1HB),w3(1HB)]-strip view of (4,3)D HCCH spectra in the order of central spectum and two basic spectra (sub-spectrum I and II) before (left) and after (right) peak position adjustment of w(13CG) and w(13CG+-1HG).
    XEASY cos4.jpg XEASY cos3.jpg

  7. In UBNMR, run updatecosyGFT to update the SQ shifts and to generate a PeakList (e.g. hcchcosyI3.peaks) using updated shifts.
  8. In XEASY, repeat step 6 for [w1(13C;1H),w3(1HG)]-strips (CG/HG strips) to confirm peaks accurately positioned at w(13CG+- 1HG) as show in Figure 4, re-adjust the peak position if necessary.
    Figure 4 Example of GFT (4,3)D HCCH analysis. Side by side [w2(13C;1HB),w3(1HB)] (first three strips) and [w2(13C;1HG),w3(1HG)] (later three stips) strips view of (4,3)D HCCH spectra in the order of central spectum and two basic spectra (sub-spectrum I and II). Top: before peak position adjustment and before updating shifts of w(13CG) and w(13CG+-1HG); Middle: after peak position adjustment and before updating shifts of w(13CG) and w(13CG+-1HG); Bottom: after peak position adjustment and after updating shifts of w(13CG) and w(13CG+-1HG).

    XEASY cos5.jpg

    XEASY cos6.jpg

    XEASY cos7.jpg


  9. In UBNMR, run updatecosyGFT to update SQ shifts and to generate a PeakList (e.g. hcchcosyI4.peaks) using updated shifts.
  10. In XEASY, for peripheral nuclei beyond CD, repeat steps 6-8 for [w1(13C;1H),w3(1HD)]-strips (CD/HD strips) to obtain peaks accurately positioned at w(13CD+- 1HD).
  11. In UBNMR, run updatecosyGFT to update SQ shifts and to generate a PeakList (e.g. hcchcosyI5.peaks) using updated shifts.
  12. In XEASY, for Main.CE, repeat steps 6-8, confirm that peaks are accurately positioned in [w1(13C;1H),w3(1HE)]-strips (Main.CE/HE strips) at w(13Main.CE+-E 1HE).
  13. In XEASY, use ac, wc and wp to save final AtomList (e.g. hcchcosyO5.prot) and PeakList (e.g. hcchcosyO5.peaks).
  14. In UBNMR, run cleanCosyGft. The AtomList is loaded, SQ 1H and 13C are calculated and written into the AtomList while all linear combinations of shifts are removed. The thus generated final SQ-only AtomList as bbsc.prot contains all shifts and represents the starting point for analysis of NOESY.


Approach II, using mr to find the SQ chemical shift

The procotol describe above updates SQ shifts step by step by using UBNMR, there is another approach that update SQ shifts step by step by directly adjust peak position with two orthognal display of the spectra. The later approach is given below and may require addtional experience.

  1. In XEASY, use ns to load two copies of the three sub-spectra of (4,3) HCCHand display the two copies in two orthogonal view (X-axis: w3(1H), Y-axis: w2(13C;1H), Z-axis: w1(13C) & X-axis: w3(13C), Y-axis:w2(1H), z: w3(1H)); use ls, lc, and lp, respectively, to load corresponding SequenceList, AtomList and PeakList; use sp, se and gs to display [w3(1H),w2(13C;1H)]-strips and [w1(13C),w2(13C;1H)]-strips for the interested residue.
  2. In XEASY, start from both views (X-axis: w3(1H), Y-axis: w2(13C;1H), Z-axis: w1(13C) & X-axis: w1(13C), Y-axis, w2(13C;1H), Z-axis: w3(1H)) of both HA and HB strips, confirm the assignment for CA/HA and CB/HB, and re-adjust the peak position if necessarys by using mr. Adjusting peak position along X-axis in strips with view (X-axis: w3(1H), Y-axis: w2(13C;1H), Z-axis: w1(13C) for 1H chemical shift, along X-axis in strips with view (X-axis: w1(13C), Y-axis: w2(13C;1H), Z-axis: w3(1H))for 13C chemical shifts, along Y-axis for GFT chemical shift w(C13;1H).
  3. In XEASY, go to the CB/HB strips, use mr to accurately position peaks at w(13CG) and w(13CG+-1HG) along Y-axis (Figure 5). While adjusting the peak position, pay attention to wether the chemical shift difference value of (w(13CG+HG) - w(13CG-HG)) become larger or small. In the case of water frenquency as the GFT dimension 1H carrier, the real HG shift is smaller than initial value if the difference become small (more negative).
    Figure 5. Example of GFT (4,3)D HCCH analysis. Side by side orthogonal views of [w2(13C;1HB),w3(1HB)] and [w2(13C;1HG),w3(1HG)] strips view of (4,3)D HCCH spectra in the order of X-axis as HB central, CB central, HB sub-spectrum I, CB sub-spectrum I, HB sub-spectrum II, CB sub-spectrum II, HG central, CG central, HG sub-spectrum I, CG sub-spectrum I, HG sub-spectrum II and CG sub-spectrum II. Top: before peak position adjustment of w(13CG) and w(13CG+-1HG); Middle: after peak position adjustment of w(13CG) and w(13CG+-1HG) and before resorting strips; Bottom: after peak position adjustment of w(13CG), and w(13CG+-1HG) and after resorting strips.
    XEASY cos8.jpg
    XEASY cos9.jpg
    XEASY cos10.jpg

  4. In XEASY use es, se and gs to resort and redisplay the strips in both view. As shown in Figure 6 Top, since the chemical shift of CG is correct, identify the correct HG chemical shift in the HG strip with view (X-axis: w3(13C), Y-axis:w2(1H); z: w3(1H)) based on the expected peak pattern and expected HG peak position. Use zd to duplicate the HG central strip if the expect position of HG is out of this strip view (Figure 6, middle), use gl or gr to change the dulpilicate strip view to find the expected HG.
    Figure 6. Example of GFT (4,3)D HCCH analysis. Side by side orthogonal views of [w2(13C;1HB),w3(1HB)] and [w2(13C;1HG),w3(1HG)] strips view of (4,3)D HCCH spectra in the order of X-axis as HB central, CB central, HB sub-spectrum I, CB sub-spectrum I, HB sub-spectrum II, CB sub-spectrum II, HG central (has two HG central in the top and middle figures), CG central, HG sub-spectrum I, CG sub-spectrum I, HG sub-spectrum II and CG sub-spectrum II. Top: after zd, and before peak position adjustment of w(1HG); Middle: after zd, after peak position adjustment of w(1HG) and before resorting strips; Bottom: after peak position adjustment of w(1HG) and after resorting strips.
    XEASY cos11.jpg
    XEASY cos12.jpg
    XEASY cos13.jpg

  5. In XEASY use es, se and gs to resort and redisplay the strips in both view, and confirm the CG/HG shift if both view.
  6. In XEASY, repeat the above steps for assignments of peripheral nuclei beyond CG of this residue, such as CD/HD and Main.CE/HE.
  7. In XEASU, use sp to select peaks for the next residues and repeat the above steps for assignments.
  8. In XEASY, use ac, wc and wp to save final AtomList and PeakList.
  9. In UBNMR, run cleanCosyGft. The AtomList is loaded, SQ 1H and 13C are calculated and written into the AtomList while all linear combinations of shifts are removed. The thus generated final SQ-only AtomList as bbsc.prot contains all shifts and represents the starting point for analysis of NOESY.


-- Main.GaohuaLiu - 31 Jan 2007@