This might help: use this as a shell script replace $i with model
name. This will find the best three models and move them into a
directory best3 and superimpose them. It will also move the best
model into directory best1. You will just need to replace the
impose section with loop modeling and add the 4th and 5th models.
>>>>
grep OBJECTIVE $i.B9999????|cut -d':'
-f1 > temp1
grep OBJECTIVE $i.B9999????|cut -d' ' -f10 > temp2
paste temp1 temp2 >temp3
sort -n +1 temp3 >temp4
sed -n '1,3w temp5' temp4
cut -f1 temp5 >temp6
sed -n '1w temp6a' temp6
sed -n '2w temp6b' temp6
sed -n '3w temp6c' temp6
read SEQ1 < temp6a
read SEQ2 < temp6b
read SEQ3 < temp6c
echo $SEQ1
echo $SEQ2
echo $SEQ3
cp $SEQ1 best1/$i.pdb
echo "Finished"
rm temp1 temp2 temp3 temp4 temp5 temp6 temp6a temp6b temp6c
echo "#Superimpose models
SET ATOM_FILES = '$SEQ1' '$SEQ2' '$SEQ3'
SET WRITE_FIT = ON
MALIGN3D" > impose$i.top
mod4 impose$i.top
cp $SEQ1 $SEQ2 $SEQ3 best3
cp $i.imposed impbest3
<<<<<<
At 09:19 AM 1/10/02 -0500, Andrej Sali wrote:
Yes,
but it would require running a unix script out of the top file to do the
sorting and file renaming.
Andrej
-----Original Message-----
From: owner-modeller_usage@salilab.org
[mailto:owner-modeller_usage@salilab.org]
On Behalf Of Dan Thomas Major
Sent: Thursday, January 10, 2002 6:53 AM
To: modeller_usage
Subject: stepwise top
Hi,
Is it possible in a modeller top file to calculate, let's say 100 models,
and then tell modeller to refine the best 10 and then to add loops for
the best 5 (all in one run).
Thanks,
Dan
--
Dan Thomas Major (at Dr. B. Fischer's lab)
Bar-Ilan University
Ramat-Gan, Israel
************************************************
Michael Buck
NCSU Genetics
mjbuck@unity.ncsu.edu
Phone (919)515-5759
Fax
(919)515-3355
*************************************************