-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathT1rhoN15_tr_calib.cw
332 lines (294 loc) · 7.02 KB
/
T1rhoN15_tr_calib.cw
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
; T1rhoN15_tr_calib.cw
;
; calibration of 15N spin-lock by nutation method
; adapted from Guenneugues (1999) J Magn Reson 136 118-126
;
; 15N-T1rho relaxation experiment with TROSY read-out
; for 15N, 15N13C and 2H15N13C labelled proteins
; written by NL 10/25/11
; see footnotes
#include <Avance.incl>
#include <Grad.incl>
#include <Delay.incl>
;#define LABEL_CN ; switch on for 13C labelled samples
#define TEMP_COMPENSATION
"in0=inf1*0.5"
# ifdef LABEL_CN
"d0=97u-p4*2+p7*0.66-p1*0.5"
#else
"d0=100u+p7*0.66-p1*0.5"
#endif /*LABEL_CN*/
"d11=30m"
"DELTA=2.65m"
"DELTA1=2.65m"
"DELTA2=2.65m-p22-p11-300u"
"DELTA3=2.65m-p23-p10-300u"
"DELTA4=260u-p24-p1*0.66"
"d27=p24+35u"
"l1=1"
"l2=1"
"cnst21=176"
"cnst22=56"
"cnst18=-800"
"spoff4=bf2*((cnst22-cnst21)/1000000)"
1 ze
1m
2 d11 ;do:f2
1m LOCKH_OFF
3m
3 1m
1m
4 3m
5 2m BLKGRAD
10u pl1:f1
10u ;pl4:f2
10u pl7:f3
(p7 ph0):f3
10u
;---------temperature compensation---------
# ifdef TEMP_COMPENSATION
"p17=p18+1m-vp"
"d17=d1-p17"
# endif
"d8=vp*0.25-p1"
"d9=vp*0.25-p1-6u"
"d29=p8"
;-------- Temperature compensation and d1 recovery delay----------------
#ifdef TEMP_COMPENSATION
10u fq=cnst18(bf ppm):f3
10u pl8:f3
(p17 ph0):f3 ; 15N pulse is applied far off-resonance
10u
10u fq=0:f3
d17
#else
d1
#endif
1m UNBLKGRAD
10u pl7:f3
;------- kill steady state 15N ------------
(p7 ph0):f3
5u
p20:gp6
200u
;------- first INEPT Hz-> 2HxNz -----------
(p1 ph0):f1
5u
DELTA gron0 ; soft gradient to prevent radiation damping
5u groff
(center(p1*2 ph0):f1 (p7*2 ph0):f3)
5u
DELTA gron0
5u groff
;------- rephase 2HxNz to Nz------ ---------
(p1 ph5):f1 (p7 ph0):f3
5u
DELTA1 gron1 ; soft gradient to prevent radiation damping
5u groff
(center (p1*2 ph0):f1 (p7*2 ph0):f3)
5u
DELTA1 gron1
5u groff
(p7 ph6):f3
5u
;---------------------------------------------
(p1 ph2):f1 ; purge pulse to kill any residual HzNz
;goto 999
5u
p21:gp7 ; cleaning gradient
100u
100u pl8:f3 ; change to pl0:120 dB for Avance I or DRX systems
;------ N15 T1rho relaxation period-----------
# ifdef LABEL_CN
;(d8 p1 ph0 3u p1*2.3 ph1 3u p1 ph0 3u d9*0.5 gron9 3u groff 6u d9*0.5 gron9*-1 3u groff 6u p1 ph0 3u p1*2.3 ph1 3u p1 ph0):f1 (p8:sp8 ph0 3u pl8 vp ph0 p8:sp9 ph0):f3 (d29 d9 p4*2 ph0 3u 3u pl2 p4*2:sp4 ph0):f2
(d8 p1*2 ph0 3u d9 gron9 6u groff d9 gron9*-1 3u groff p1*2 ph0):f1 (vp ph0):f3 ;NB no carbon implemented yet!
#else
;(d8 p1 ph0 3u p1*2.3 ph1 3u p1 ph0 3u d9*0.5 gron9 3u groff 6u d9*0.5 gron9*-1 3u groff 6u p1 ph0 3u p1*2.3 ph1 3u p1 ph0):f1
if "vp<1m"
{
(d8 p1*2 ph0 3u d9 6u d9 3u p1*2 ph0):f1 (vp ph0):f3
}
else
{
(d8 p1*2 ph0 3u d9 gron9 6u groff d9 gron9*-1 3u groff p1*2 ph0):f1 (vp ph0):f3
}
#endif
5u
p21:gp8
200u
;------Echo/ Anti-echo encoding for TROSY read-out------------
3u pl7:f3
3u ;pl4:f2
3u pl1:f1
if "l1==1"
{
(p7 ph7):f3
10u
p25:gp5
200u
(p7*2 ph7):f3
10u
p25:gp5*-1
}
else
{
(p7 ph17):f3
10u
p25:gp5*-1
200u
(p7*2 ph17):f3
10u
p25:gp5
}
;------ t1 (15N) evolution period ------------------------------
d0
# ifdef LABEL_CN
(p4*2 ph0 3u 3u pl2 p4*2:sp4 ph0):f2
# endif
d0
;------ start TROSY read-out------------------------------------
if "l1==1"
{
(p1 ph1):f1 ; Echo
3u
3u pl0:f1
(p11:sp11 ph11:r):f1
6u
}
else
{
(p1 ph3):f1 ; Anti-Echo
3u
3u pl0:f1
(p11:sp11 ph13:r):f1
6u
}
5u pl1:f1
;goto 999 ; optimization of water supression (sp11)
DELTA2
p22:gp2
300u
(center (p1*2 ph0):f1 (p7*2 ph0):f3)
7u
p22:gp2
DELTA2
300u pl0:f1
;-------------------------------------------------
(p11:sp12 ph12:r):f1
5u
3u pl1:f1
if "l1==1"
{
(p1 ph0):f1 (p7 ph1):f3 ; Echo
}
else
{
(p1 ph0):f1 (p7 ph3):f3 ; Anti-Echo
}
;goto 999 ; for optimization of water supression (sp12)
DELTA3
p23:gp3
200u
100u pl10:f1
(center(p10 ph10:r 5u pl1 p1*2 ph0 5u pl10 p10 ph10:r):f1 (p7*2 ph0 d27):f3)
5u
;goto 999 ; for optimization of water supression (pl10)
p23:gp3
DELTA3
DELTA4
(p7 ph0):f3
5u
p24:gp4 ; Echo/Anti-echo decoding gradient
999 5u
5u ;pl31:f2
20u BLKGRAMP
go=2 ph31 ;cpds2:f2
1m ;do:f2
d11 wr #0 if #0 zd
1m LOCKH_OFF
2m ivp
lo to 3 times l6
;3m iu1
;lo to 4 times 2
;1m id0
;1m ru1
;lo to 5 times l3
1m
1m BLKGRAD
exit
ph0=0
ph1=1 ;check right phase for Boltzmann !!!!!
ph2=2
ph3=3
ph5= 1
ph6= 1 1 1 1 3 3 3 3
ph9=2
ph10=2
ph11=3
ph12=0
ph13=1
ph7=1 0 3 2
ph17=1 2 3 0
ph31=1 2 3 0 3 0 1 2
;-------------NOTES----------------------
;o1p = 4.7 ppm
;o2p=176 ppm (CO)
;o3p=119 ppm
;NS=8*n
;in0=inf/2
;SW=1/(2*in0)
;echo-antiecho in N15 (process as Complex in NmrDraw before splitting the spectra)
; 1H pulses
;p1: 90 deg hard 1H pulse @pl1
;pl1: 1H 90 deg
;pl0: 120 dB
;p10: 1200u (@ 600 MHz) 90 deg soft rectangular water flip-back pulse
;p11: 1900u (@ 600 MHz) 90 deg Sinc1.1000 water flip-back pulse (sp11,sp12)
;sp11: 90 deg Sinc1.1000 water flip-back pulse
;sp12: 90 deg Sinc1.1000 water flip-back pulse
;spnam11 Sinc1.1000
;spnam12 Sinc1.1000
; 13C pulses
;p4: 13CO selective 180 deg (23.7*2us @ 600 MHz) @pl4
;pl4: 13C 90 deg
;sp4: 13CA selective 180 deg (23.7*2us @ 600 MHz)
;CPDPRG2: garp (aq C' decoupling)
;pcpd5: C' decoupling (140u or 280u @pl31)
;pl31: C' decoupling power
;15N pulses
;p7 : 90 deg hard 15N pulse @pl7
;pl7 :15N 90 deg
;p8: adiabatic half-passage pulse (sp8,sp9) to rotate magnetization along spin-lock axis (sp8) and back to z (sp9)
;p18: maximum duration of spin-lock; temperature compensation
;pl8: 15N spin-lock power
;sp8 @pl8 spin-lock power
;sp9 @pl8 spin-lock power
;spnam8 AHP TanhTan (1st half)
;spnam9 AHP TanhTan (2nd half)
;spoal8 1
;spoal9 0
; gradients
;p20: 1000u
;p21: 200u
;p22: 300u
;p23: 1000u
;p24: 60.8u Echo/Anti-echo decoding gradient
;p25: 300u Echo/Anti-echo half-encoding gradient
;for z-only gradients
;gpz0: 3%
;gpz1: 2%
;gpz2: 10%
;gpz3: 50%
;gpz4: 33%
;gpz5: -33%
;gpz6: 30%
;gpz7: -50%
;gpz8: 40%
;gpz9: 0.5%
;gpnam2 SINE.10
;gpnam3 SINE.50
;gpnam4 SINE.10
;gpnam5 SINE.10
;gpnam6 SINE.50
;gpnam7 SINE.10
;gpman8 SINE.10