Difference between revisions of "Custom peak list"

From relax wiki
Jump to navigation Jump to search
(Forced creation of a TOC - this will improve the formatting on the main page 'Did you know...' section.)
 
(One intermediate revision by the same user not shown)
Line 1: Line 1:
 +
__TOC__
 +
 
== Custom peak list ==
 
== Custom peak list ==
 
<source lang="text">
 
<source lang="text">
Line 53: Line 55:
 
= See also =
 
= See also =
  
[[Category:NOE]]
 
 
[[Category:Peak_lists]]
 
[[Category:Peak_lists]]
 
[[Category:Relaxation]]
 
[[Category:Relaxation]]
[[Category:Relaxation_dispersion]]
+
[[Category:Relaxation dispersion analysis]]
 
[[Category:Relaxation_rates]]
 
[[Category:Relaxation_rates]]
 +
[[Category:Steady-state NOE analysis]]

Latest revision as of 13:46, 15 October 2015

Custom peak list

protein 3 R 3 N 3.642647e+06 6.056554e+06 3.753433e+06 ...
protein 5 E 5 N 1.609356e+06 2.927111e+06 1.726433e+06 ...
protein 6 V 6 N 1.697771e+06 3.015788e+06 1.771777e+06 ...
protein 7 N 7 N 1.535896e+06 3.005234e+06 1.683477e+06 ...
protein 8 I 8 N 1.332059e+06 2.611723e+06 1.519182e+06 ...
protein 9 V 9 N 1.532644e+06 2.834876e+06 1.597882e+06 ...
protein 11 N 11 N 1.219399e+06 2.479556e+06 1.289455e+06 ...
protein 15 D 15 N 2.808704e+06 4.812759e+06 2.935094e+06 ...
protein 16 Q 16 N 3.410860e+06 5.823456e+06 3.427248e+06 ...
protein 22 E 22 N 3.681818e+06 6.378086e+06 3.663392e+06 ...
 ...

Sequence read

sequence.read(file='Custom_model.txt', dir=None, spin_id_col=None, mol_name_col=1, res_num_col=2, res_name_col=3, spin_num_col=4, spin_name_col=5, sep=None, spin_id=None)

Use awk to make a selected residues file

tcsh ;
set IN=table_ser_files_model.txt ;
set SELRESIS=sel_resi.txt ;
awk '{print $2}' $IN > $SELRESIS ;
gedit $SELRESIS

After this step, you should remove those lines with residues you dont want to have in the new file.

tcsh ;
set IN=table_ser_files_model.txt ;
set SELRESIS=sel_resi.txt ;
set RESIS=`cat $SELRESIS` ;
set DATE=`date '+%Y%m%d_%H%M'` ;
set TMP=${IN}_${DATE} ;
cp $IN $TMP ;
rm $IN ;
set T=$ ;
echo $RESIS ;
foreach RESI ($RESIS) ;
cat  $TMP | awk "{if (${T}2 == $RESI ) { print ${T}0; } }" >> $IN
end ;
cat $IN ;

Make an execution script

Start relax GUI with logfile

relax_disp -g -l LOGFILE.txt


See also