We have written a complete suite of Varian protein experiments in an easy-to-understand style.
The first sequences and basic architecture were developed in
a collaboration between Erik Zuiderweg and Ananya Majumdar.
The sequences all use the same symbols, and can be easily translated to Bruker
and vice-versa.
Main philosophy: one sequence per experiment. No flags or if statements.
An example, Watergate fast-HSQC is shown below. If you are interested in this approach,
you should contact Erik Zuiderweg (zuiderwe@umich.edu).
It is free for academia, but you need to sign some documents with respect
to liability.
We have available
3D Triple resonance experiments
3D NOESYs in all kinds of flavors
3D 13C correlations in all kinds of flavors
15N, 1H and 13C Relaxation experiments
RDC experiments
Water NOE/ROE experiments
/* hsqc-up.c */
#include <standard.h>
static double d2_init = 0.0;
static int phi1[2] = {0,2},
psirec[2]={0,2};
pulsesequence()
{
#include "erpz_include.h"
/* Load Phase Table */
settable(t11,2,phi1);
settable(t31,2,psirec);
/* create phase correction for watergate */
if(phcor12 < 360) phcor12 = phcor12+360;
initval(phcor12,v12);
/* Phase incrementation for hypercomplex data */
if ( phase1 == 2 ) /* Hypercomplex in t1
*/
tsadd(t11, 1, 4);
/* calculate modification to phases based on current t1 values
to achieve States-TPPI acquisition */
if(ix == 1)
d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5);
if(t1_counter %2) {
tsadd(t11,2,4);
tsadd(t31,2,4);
}
tau1 = d2/2.;
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
rcvroff();
delay(5e-6); /* recycel delay time */
delay(del18);
obspower(hpwr0); xmtrphase(zero) ; decpower(cpwr11);
dec2power(npwr0);
delay(10.e-6);
dec2rgpulse(p90n,zero,2.0e-6,0.);
delay(10.e-6);
zgradpulse(grpwr0,grdel0);
delay(del10);
rgpulse(p90h,zero,2.0e-6,0.); /* INEPT */
delay(10.e-6);
zgradpulse(grpwr10,grdel10);
delay(del3-grdel10-10.e-6);
rgpulse(p90h*2.0,one,2.0e-6,0.);
dec2rgpulse(p90n*2.0,zero,2.0e-6,0.);
delay(10.e-6);
zgradpulse(grpwr10,grdel10);
delay(del3-grdel10-10.e-6-p90n*2);
rgpulse(p90h,one,2.0e-6,0.);
delay(10.e-6);
zgradpulse(grpwr1,grdel1);
delay(del10);
dec2rgpulse(p90n,t11,2.0e-6,0.); /* evolution with 13C dec without 1st order phase corr */
delay(tau1);
decshaped_pulse(shname6,shpul2,zero,2.0e-6,0.); /* 180 CO from Ca */
rgpulse(p90h*2.0,one,2.0e-6,0.);
decshaped_pulse(shname2,shpul2,zero,2.0e-6,0.); /* 180 CA */
delay(tau1);
dec2rgpulse(p90n*2.0,zero,2.0e-6,0.);
delay(d2_init+shpul4+shpul6+p90h*2.);
dec2rgpulse(p90n,zero,2.0e-6,0.);
delay(10.e-6);
zgradpulse(grpwr2,grdel2);
delay(del10);
rgpulse(p90h,two,2.0e-6,0.); /*
reverse INEPT with water gate */
delay(del4-grdel3-del10-pul8-20.e-6-p90n*2);
delay(10.e-6);
zgradpulse(grpwr3,grdel3);
delay(del10);
obspower(hpwr1); obsstepsize(1.0); xmtrphase(v12) ; txphase(two);
delay(5.e-6);
rgpulse(pul8,two,2.0e-6,0.);
obspower(hpwr0); xmtrphase(zero) ;
delay(5.e-6);
delay(p90n*2.0);
rgpulse(p90h*2.0,zero,2.0e-6,0.);
dec2rgpulse(p90n*2.0,zero,2.0e-6,0.);
obspower(hpwr1); obsstepsize(1.0); xmtrphase(v12) ; txphase(two);
delay(5.e-6);
rgpulse(pul8,two,2.0e-6,0.);
delay(10.e-6);
zgradpulse(grpwr3,grdel3);
delay(del10);
dec2power(cpdpwr2);
delay(del4-grdel3-del10-pul8-20.e-6-p90n*2-15.e-6);
status(B);
setreceiver(t31);
}