Difference between revisions of "Lab Journal of Task 3 (MSUD)"

From Bioinformatikpedia
m (Disordered protein)
(Disordered protein)
Line 26: Line 26:
 
done
 
done
 
done
 
done
  +
</source>
  +
  +
* '''MetaDisorder(MD)'''
  +
** As the man page of <tt>metadisorder</tt> describes, the prediction of disordered region is based on the results of other programs such as NORSnet, PROFbval etc. Rather than directly call <tt>metadisorder</tt> we have used the wrapper program <tt>predictprotein</tt> as is described on the exercise page. Following script was called for the task:
  +
  +
<source lang="bash">
  +
#!/bin/sh -e
  +
INPUT=$HOME/master-practical-2013/task03
  +
OUTPUT=$HOME/master-practical-2013/task03/metadisorder
  +
EXE=predictprotein
  +
#make output directory
  +
if [ ! -d $OUTPUT ]; then
  +
mkdir $OUTPUT
  +
fi
  +
#call metadisorder for all query sequences
  +
for seq in $INPUT/*.fasta
  +
do
  +
filename=`basename $seq`
  +
$EXE --seqfile $seq --target metadisorder -p metadisorder_$filename \
  +
-o $OUTPUT
  +
done
  +
echo Done!
 
</source>
 
</source>
   

Revision as of 19:41, 9 May 2013

Secondary structure

Disordered protein

  • IUPred
    • Predictions were performed through the web server of IUPred. Graphical profiles of the results were downloaded.
    • Output of IUPred are stored in the directory /mnt/home/student/weish/master-practical-2013/task03/02-disordered-protein/iupred
    • We have also performed the prediction from command-line, following is the bash script:

<source lang="bash">

  1. !/bin/sh -e

INPUT=$HOME/master-practical-2013/task03 OUTPUT=$HOME/master-practical-2013/task03/02-disordered-protein/iupred PARAMS="long short glob"

if [ ! -d $OUTPUT ]; then

       mkdir $OUTPUT

fi

for seq in $INPUT/*.fasta do

       filename=`basename $seq`
       for param in $PARAMS
       do
               iupred $seq $param > $OUTPUT/iupred_${filename}_$param.tsv
       done

done </source>

  • MetaDisorder(MD)
    • As the man page of metadisorder describes, the prediction of disordered region is based on the results of other programs such as NORSnet, PROFbval etc. Rather than directly call metadisorder we have used the wrapper program predictprotein as is described on the exercise page. Following script was called for the task:

<source lang="bash">

  1. !/bin/sh -e

INPUT=$HOME/master-practical-2013/task03 OUTPUT=$HOME/master-practical-2013/task03/metadisorder EXE=predictprotein

  1. make output directory

if [ ! -d $OUTPUT ]; then

       mkdir $OUTPUT

fi

  1. call metadisorder for all query sequences

for seq in $INPUT/*.fasta do

       filename=`basename $seq`
       $EXE --seqfile $seq --target metadisorder -p metadisorder_$filename \
       -o $OUTPUT

done echo Done! </source>

Transmembrane helices

Signal peptides

GO terms