Re: [phenixbb] rosetta fragment files -- adjusting residue numbers to combine them

...

[[email protected]] On Behalf Of Leonid Sazanov [[email protected]] Sent: 16 June 2011 13:11 To: [email protected] Subject: [phenixbb] rosetta fragment files

Hi, I am trying to use mr_rosetta with phenix-1.7.1-743 and rosetta 3.2.1. All rosetta tests passed OK. I have several large protein chains in the starting model (already

...

So I had to get several pairs of fragment files from robetta server and

Hi Leonid and others wanting to use phenix.mr_rosetta with large chains, After looking at the Rosetta code and conferring with Frank DiMaio from the Baker laboratory... You can paste several fragment files together to represent your entire chain, but you need to edit all but the first one to set the position in the sequence file that each set of fragments goes with. Here is a little script that will take a fragments file that you have created that started with residue 651 and edit it in a way that you can just add to the end of your fragments file for residues 1-650. Run it with: phenix.python adjust_robetta_resid.py \ fragments_in fragments_adjusted 651 -------------------------------- # adjust_robetta_resid.py # use: phenix.python adjust_robetta_resid.py fragments_in fragments_adjusted 651 # to adjust the first residue of fragments_in to be 651 import sys args=sys.argv[1:] if len(args) != 3: print "USE: phenix.python adjust_robetta_resid.py fragments_in fragments_adjusted " sys.exit() file_in=args[0] file_out=args[1] start_resno=int(args[2]) print "Copying %s to %s offsetting start residue number to %d" %( file_in,file_out,start_resno) f=open(file_out,'w') new_resno=0 for line in open(file_in).readlines(): if line.find('position:') > -1 and line.find('neighbors:')>-1: spl=line.split() current_resno=int(spl[1]) new_resno=current_resno+start_resno print "Residue %d changed to %d" %(current_resno,new_resno) spl[1]=" %d " %(new_resno) new_line=" ".join(spl) else: new_line=line print >>f, new_line.rstrip() f.close() -------------------------------- I will add this to phenix as "phenix.adjust_robetta_resid" as well. However...there still is a problem if you have multiple chains...as mr_rosetta has only one set of fragments files the same fragments files are being applied to all your chains. I will add a way to specify which fragments files go with which chains. All the best, Tom T placed). list all of them in the script.

...

Then mr_rosetta fails at the REBUILD_1 stage with the message: could not open file (fragment file), although file is definitely there. Is the reason that one can have only two input fragment files (3/9)? Can several fragment files be combined into two 3/9 ones? Thanks

-- Dr. Leonid A. Sazanov Research group leader Medical Research Council Mitochondrial Biology Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY

_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb _______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb