Index of /ftp/pub/bmrb/entry_lists/timedomain_entries/bmr6464/timedomain/analysis/3garant
Procedure for running GARANT
1. Copy the initial.gar, comp.gar and write.gar default scripts from
/data/elgar/cesg/methods/garant/ into /data/elgar/target#/3garant dir,
along with [target].seq and [target].prot and peaklists from 2spscan.
2. Make directories 1 (initial.gar output) and 1a,1b,1c (comp.gar output)
3. Execute the first garant run by typing in the main garant dir:
garant initial 1 [seed] [target#] >! run1.out &
4. When the first run is complete, begin three comparison runs by typing:
garant comp 1a 1 [seed] [target#] >! run1a.out &
^
|
Substitute 1b, 1c and a different seed for other two runs.
5. After verifying the assignments from the first run by looking at strips in
hnca, cconh, hncoca, etc. the regions with correct assignments can be fixed
using the macro initial_fix.gar. There are three important sets of changes to
the Garant macro that must be made. (1) Peaklists with verified assignments are
named hnca_ass, cconh_ass, etc., and are loaded with the 'load assignments'
commands near the bottom of the macro. (2) Make sure that a 'reset user assignments'
command is executed for each peaklist that contains assignments to be held fixed.
(3) A section from the initial_fix.gar file below illustrates how to fix assignments
for a portion of the sequence, in this case residues 20 through 52.
# insert here fixed assignments
do i 20 52
fix assignments "HN($i)"
fix assignments "N($i)"
# fix assignments "C($i)"
fix assignments "CA($i)"
fix assignments "CB($i)"
end do
# insert here fixed assignments
do i 20 52
fix assignments "*($i)"
end do