parameter suffix suffix2 rseed

set dir      = /data/elgar/cesg/at3g51030/3garant/6
set dir2     = /data/elgar/cesg/at3g51030/3garant/6/{$suffix}
set dir3     = /data/elgar/cesg/at3g51030/3garant/6/{$suffix2}

# Input data
set atomlist = at3g51030
set sequence = at3g51030
#set cbpl         = hncacb 	# peak lists hncacb
set hsqcpl       = hs15n
set hncapl       = hnca
set hncocapl     = hncoca
#set hncopl       = hnco
#set cconhpl      = cconh
#set hcconhpl      = hcconh
#set hncacopl     = hncaco
#set hncocacbpl  = hncocacb
#set n15noesypl   = no15n_na_new.peaks~	# peak lists n15noesy
#set c13noesypl   = c13no2_na_new
#set hcchpl       = hcch_sol_unfold

set refhsqcpl       = hs15n_{$suffix2}
set refhncapl       = hnca_{$suffix2}
set refhncocapl     = hncoca_{$suffix2}
#set refhncopl       = hnco_{$suffix2}
#set refcconhpl      = cconh_{$suffix2}
#set refhcconhpl     = hcconh_{$suffix2}
#set refhncocacbpl   = hncocacb_{$suffix2}
#set refcbpl         = hncacb_{$suffix2}
#set refhncacopl     = hncaco_{$suffix2}

# macro to be used for writing out best assignment
set savemacro   = write_042204 

set name     = {$suffix}
set sizePop  = 100

seed random numbers $rseed

N15H1
   set N15H1_w1         = 2
   set N15H1_w2         = 1
#   set N15H1_acc_w1     = 0.02
#   set N15H1_acc_w2     = 0.2
   load peaks N15H1 $dir/$hsqcpl   

HNCA
#   set HNCA_acc_w1      = 0.04
#   set HNCA_acc_w2      = 0.15
#   set HNCA_acc_w3      = 0.05
   set HNCA_w1          = 1
   set HNCA_w2          = 3
   set HNCA_w3          = 2
   load peaks HNCA      $dir/$hncapl

HN_CO_CA
#    set HN_CO_CA_acc_w1  = 0.02
#    set HN_CO_CA_acc_w2  = 0.2
#    set HN_CO_CA_acc_w3  = 0.2
    set HN_CO_CA_w1      = 1
    set HN_CO_CA_w2      = 3
    set HN_CO_CA_w3      = 2
    load peaks HN_CO_CA  $dir/$hncocapl

#CBCACONH
#    set CBCACONH_acc_w1  = 0.02
#    set CBCACONH_acc_w2  = 0.2
#    set CBCACONH_acc_w3  = 0.2
#    set CBCACONH_w1      = 1
#    set CBCACONH_w2      = 3
#    set CBCACONH_w3      = 2
#    load peaks CBCACONH  $dir/$hncocacbpl

#HNCO
#    set HNCO_acc_w1      = 0.04
#    set HNCO_acc_w2      = 0.15
#    set HNCO_acc_w3      = 0.05
#    set HNCO_w1          = 1
#    set HNCO_w2          = 3
#    set HNCO_w3          = 2
#    load peaks HNCO      $dir/$hncopl

#CBCANH
#   set CBCANH_w1     = 1	# w1 of XEASY-list -> w1 (garant) == HN
#   set CBCANH_w2     = 3	# w3 of XEASY-list -> w2 (garant) == N
#   set CBCANH_w3     = 2	# w2 of XEASY-list -> w3 (garant) == CA/CB
#   set CBCANH_acc_w1 = 0.04	# 
#   set CBCANH_acc_w2 = 0.12
#   set CBCANH_acc_w3 = 0.20
#   load peaks CBCANH $dir/$cbpl

#C_CO_NH
#    set C_CO_NH_acc_w1   = 0.04
#    set C_CO_NH_acc_w2   = 0.12
#    set C_CO_NH_acc_w3   = 0.2
#    set C_CO_NH_w1       = 1
#    set C_CO_NH_w2       = 3
#    set C_CO_NH_w3       = 2
#    load peaks C_CO_NH   $dir/$cconhpl

#HC_CO_NH
#    set HC_CO_NH_acc_w1   = 0.025
#    set HC_CO_NH_acc_w2   = 0.3
#    set HC_CO_NH_acc_w3   = 0.025
#    set HC_CO_NH_w1       = 1
#    set HC_CO_NH_w2       = 3
#    set HC_CO_NH_w3       = 2
#    load peaks HC_CO_NH   $dir/$hcconhpl

#HN_CA_CO
#    set HN_CA_CO_acc_w1   = 0.025
#    set HN_CA_CO_acc_w2   = 0.3
#    set HN_CA_CO_acc_w3   = 0.3
#    set HN_CA_CO_w1       = 1
#    set HN_CA_CO_w2       = 3
#    set HN_CA_CO_w3       = 2
#    load peaks HN_CA_CO   $dir/$hncacopl

#N15NOESY 
#  set N15NOESY_w1     = 3    # N15 axis
#  set N15NOESY_w2     = 2
#  set N15NOESY_w3     = 1	
#  set N15NOESY_acc_w1 = 0.12
#  set N15NOESY_acc_w2 = 0.04
#  set N15NOESY_acc_w3 = 0.05
#  load peaks N15NOESY $dir/$n15noesypl
#  set N15NOESY_imp    = 10.0    # importance of NOESY peaks

#C13NOESY
#   set C13NOESY_acc_w1 = 0.3
#   set C13NOESY_acc_w2 = 0.025
#   set C13NOESY_acc_w3 = 0.025
#   set C13NOESY_w1     = 3
#   set C13NOESY_w2     = 1
#   set C13NOESY_w3     = 2
#   load peaks C13NOESY $dir/$c13noesypl

#HCCH7
#   set HCCH7_w1     = 3
#   set HCCH7_w2     = 2
#   set HCCH7_w3     = 1
#   load peaks HCCH7 $dir/$hcchpl

#define peak importance "N15NOESY: * * *" 0.1  
#define peak importance "N15NOESY: N(+0) HN(+0) H(-1)" 1.5
#define peak importance "N15NOESY: N(+0) HN(+0) HN(+1)" 3.0
#define peak importance "N15NOESY: N(+0) HN(+0) HN(-1)" 3.0
#define peak importance "N15NOESY: N(+0) HN(+0) HA(-1)" 3.0
#define peak importance "N15NOESY: N(+0) HN(+0) HA(+0)" 3.0

set pick_acc      = 1.0 #sets maximal difference between chemical shifts of the same atom in different spectra

load sequence $dir/$sequence
load comparison sequence  $dir/$sequence 

#pseudocorrection #corrections applied when generating expected NOEs from homologous protein

define atom probability "HE(GLN)" 1.0   # exclude NH2 side chains
define atom probability "HD(ASN)" 1.0   # exclude NH2 side chains
define atom probability "HE(ARG)" 1.0   # exclude NH side chains
define atom probability "NZ(LYS)" 1.0   # exclude NH side chains
define atom probability "CD(GLN)" 1.0   # exclude NH2 side chains
define atom probability "CG(ASN)" 1.0   # exclude NH2 side chains
define atom probability "CD(GLU)" 1.0   # exclude NH2 side chains
define atom probability "CG(ASP)" 1.0   # exclude NH2 side chains

make coherences			# generate all possible coherences
make couplings			# generate couplings
#make noe			# generate noe's if no homologous structure
make peaks			# generate expected peaks

# define shifts according to secondary structure
#define secondary shifts $dir/CalcSecondarytrunc.txt

load comparison assignments N15H1 $dir3/$refhsqcpl $dir/$atomlist 0.0
load comparison assignments HNCA $dir3/$refhncapl $dir/$atomlist 0.0
#load comparison assignments HNCO $dir3/$refhncopl $dir/$atomlist 0.0
load comparison assignments HN_CO_CA $dir3/$refhncocapl $dir/$atomlist 0.0
#load comparison assignments CBCANH $dir3/$refcbpl $dir/$atomlist 0.0
#load comparison assignments CBCACONH $dir3/$refhncocacbpl $dir/$atomlist 0.0
#load comparison assignments C_CO_NH $dir3/$refcconhpl $dir/$atomlist 0.0
#load comparison assignments HC_CO_NH $dir3/$refhcconhpl $dir/$atomlist 0.0
#load comparison assignments HN_CA_CO $dir3/$refhncacopl $dir/$atomlist 0.0

check inconsistencies

#reset user assignments HNCA
#reset user assignments HNCO
#reset user assignments HN_CO_CA
#reset user assignments CBCANH
#reset user assignments C_CO_NH

# insert here fixed assignments


criteria fragment intrainter intra 0 0
criteria fragment intrainter seq 1 1
criteria fragment correct ok
criteria comment on ok
criteria fragment correct okN N
criteria comment on okN
criteria fragment correct okHN HN
criteria comment on okHN
criteria fragment correct okHA HA
criteria comment on okHA
criteria fragment correct okCA CA
criteria comment on okCA
criteria fragment correct okHB HB
criteria comment on okHB
criteria fragment correct okCB CB
criteria comment on okCB
#criteria fragment correct okCD CD
#criteria comment on okCD
#criteria fragment correct okCG CG
#criteria comment on okCG
#criteria fragment SScorrect SSok
#criteria comment on SSok
criteria fragment mapping map
criteria comment on map
criteria comment on intra
criteria threshold intra 0.6 0.8
#criteria coherence mapping mapco
#criteria comment on mapco

opt $sizePop $savemacro

quit