[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [modeller_usage] build short 1000 olioopeptid from fasta file



Thank for the script. I have collected into a module to get a smaller file. All my seqeunces files is build except the one I am sending. I get this error message. The seqeunces is sent as attachment and as the following patter MLXP where the third aminoacid can be anything.

Best regards
Knut J


              Number of sequences in the alignment:        1
Traceback (most recent call last):
  File "/home/knut/workspace/ligand_petide/src/create_seqeunce.py", line 35, in <module>
    buildseqeunce('hNatE.txt','hNatE')
  File "/home/knut/workspace/ligand_petide/src/create_seqeunce.py", line 27, in buildseqeunce
    mdl.build_sequence(line.rstrip('\n'))
  File "/usr/lib/python2.6/dist-packages/modeller/model.py", line 133, in build_sequence
    self.generate_topology(a[0], patch_default=patch_default)
  File "/usr/lib/python2.6/dist-packages/modeller/model.py", line 238, in generate_topology
    self.env.libs.modpt, alnseq._num, patch_default)
_modeller.ModellerError: generat_477E> No templates in the alignment. Number of sequences in the alignment:        1




from modeller import *


def buildseqeunce(filename,out):
    env = environ()
    aln = alignment(env)
    mdl=model(env)

    env.libs.topology.read('${LIB}/top_heav.lib')
    env.libs.parameters.read('${LIB}/par.lib')
   
    file=open(filename,'r')
   
  
   

   
    for (n, line) in enumerate(file):
       
       
        mdl.build_sequence(line.rstrip('\n'))
        mdl.write(file=out + "%g.pdb" %n)
    file.close()   
 

buildseqeunce('hNatE.txt','hNatE')
















> Date: Tue, 6 Apr 2010 15:39:34 -0700
> From:
> To:
> CC:
> Subject: Re: [modeller_usage] build short 1000 olioopeptid from fasta file
>
> On 04/06/2010 06:24 AM, Knut J Bjuland wrote:
> > I have 1000 or more short peptid seqeunce which are 4 amino acid long.
> > Can I use modeller to build them just as a linear structure without any
> > 3d structure.
>
> If by "linear structure without any 3d structure" you mean an extended
> chain, then sure, you can use model.build_sequence():
> http://salilab.org/modeller/9v7/manual/node166.html
>
> As an example, if your 1000 sequences are in a file called seqs.txt, one
> per line, then the following script will make a set of PDB files, one
> for each sequence:
>
>
> import modeller
>
> env = modeller.environ()
> env.libs.topology.read('${LIB}/top_heav.lib')
> env.libs.parameters.read('${LIB}/par.lib')
>
> for (n, line) in enumerate(open('seqs.txt')):
> mdl = modeller.model(env)
> mdl.build_sequence(line.rstrip('\r\n'))
> mdl.write(file='seq-%d.pdb' % n)
>
>
>
>
> Ben Webb, Modeller Caretaker
> --
> http://www.salilab.org/modeller/
> Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage


Your E-mail and More On-the-Go. Get Windows Live Hotmail Free. Sign up now.
MLKp
MLRp
MLHp
MLIp
MLDp
MLYp
MLAp
MLAp
MLFp
MLNp
MLMp
MLEp
MLNp
MLEp
MLDp
MLGp
MLKp
MLKp
MLKp
MLIp
MLPp
MLDp
MLNp
MLKp
MLTp
MLWp
MLSp
MLAp
MLYp
MLNp
MLTp
MLSp
MLAp
MLCp
MLTp
MLAp
MLLp
MLTp
MLHp
MLAp
MLAp
MLCp
MLDp
MLQp
MLSp
MLDp
MLKp
MLKp
MLGp
MLHp
MLRp
MLMp
MLMp
MLHp
MLWp
MLYp
MLCp
MLGp
MLCp
MLFp
MLAp
MLYp
MLGp
MLGp
MLNp
MLKp
MLTp
MLMp
MLMp
MLHp
MLYp
MLCp
MLPp
MLVp
MLCp
MLIp
MLHp
MLCp
MLAp
MLSp
MLQp
MLIp
MLVp
MLQp
MLIp
MLIp
MLSp
MLSp
MLPp
MLHp
MLDp
MLWp
MLDp
MLEp
MLLp
MLYp
MLYp
MLPp
MLRp
MLTp
MLSp
MLFp
MLWp
MLFp
MLKp
MLAp
MLDp
MLHp
MLHp
MLKp
MLPp
MLEp
MLMp
MLEp
MLGp
MLAp
MLIp
MLIp
MLNp
MLEp
MLHp
MLLp
MLCp
MLDp
MLAp
MLWp
MLQp
MLTp
MLMp
MLAp
MLCp
MLSp
MLRp
MLRp
MLPp
MLCp
MLPp
MLKp
MLPp
MLNp
MLKp
MLCp
MLDp
MLIp
MLAp
MLCp
MLKp
MLIp
MLGp
MLDp
MLLp
MLNp
MLMp
MLMp
MLRp
MLWp
MLQp
MLPp
MLKp
MLIp
MLMp
MLHp
MLFp
MLEp
MLCp
MLQp
MLDp
MLVp
MLIp
MLEp
MLAp
MLSp
MLWp
MLDp
MLTp
MLIp
MLKp
MLLp
MLGp
MLQp
MLPp
MLLp
MLHp
MLVp
MLMp
MLHp
MLRp
MLNp
MLSp
MLIp
MLDp
MLHp
MLLp
MLTp
MLPp
MLRp
MLFp
MLVp
MLTp
MLGp
MLIp
MLSp
MLGp
MLCp
MLEp
MLHp
MLAp
MLPp
MLDp
MLDp
MLFp
MLPp
MLNp
MLYp
MLHp
MLDp
MLSp
MLCp
MLVp
MLCp
MLPp
MLLp
MLRp
MLTp
MLHp
MLPp
MLSp
MLKp
MLFp
MLEp
MLNp
MLTp
MLRp
MLDp
MLMp
MLVp
MLAp
MLSp
MLHp
MLGp
MLQp
MLIp
MLHp
MLDp
MLKp
MLIp
MLHp
MLCp
MLRp
MLEp
MLGp
MLNp
MLCp
MLWp
MLHp
MLIp
MLLp
MLNp
MLVp
MLTp
MLCp
MLLp
MLMp
MLIp
MLEp
MLYp
MLNp
MLVp
MLYp
MLLp
MLWp
MLEp
MLLp
MLWp
MLMp
MLTp
MLLp
MLIp
MLPp
MLPp
MLIp
MLVp
MLDp
MLMp
MLGp
MLPp
MLQp
MLRp
MLPp
MLLp
MLKp
MLYp
MLGp
MLPp
MLNp
MLKp
MLDp
MLYp
MLAp
MLMp
MLCp
MLLp
MLVp
MLDp
MLYp
MLVp
MLRp
MLIp
MLVp
MLRp
MLKp
MLFp
MLSp
MLKp
MLYp
MLVp
MLHp
MLQp
MLAp
MLQp
MLWp
MLKp
MLIp
MLFp
MLLp
MLSp
MLYp
MLIp
MLTp
MLYp
MLDp
MLHp
MLNp
MLWp
MLEp
MLNp
MLMp
MLRp
MLLp
MLMp
MLYp
MLHp
MLRp
MLQp
MLEp
MLLp
MLIp
MLIp
MLWp
MLIp
MLCp
MLKp
MLQp
MLSp
MLLp
MLDp
MLSp
MLDp
MLCp
MLYp
MLMp
MLCp
MLIp
MLHp
MLAp
MLRp
MLNp
MLMp
MLRp
MLIp
MLWp
MLGp
MLTp
MLVp
MLDp
MLHp
MLTp
MLQp
MLAp
MLDp
MLIp
MLPp
MLTp
MLIp
MLRp
MLGp
MLHp
MLWp
MLYp
MLIp
MLDp
MLHp
MLMp
MLKp
MLTp
MLKp
MLDp
MLEp
MLSp
MLEp
MLYp
MLVp
MLEp
MLMp
MLAp
MLDp
MLPp
MLFp
MLFp
MLKp
MLYp
MLAp
MLLp
MLWp
MLDp
MLVp
MLIp
MLEp
MLGp
MLDp
MLIp
MLEp
MLFp
MLSp
MLGp
MLYp
MLWp
MLFp
MLAp
MLFp
MLYp
MLCp
MLVp
MLMp
MLRp
MLPp
MLSp
MLCp
MLTp
MLKp
MLEp
MLTp
MLPp
MLCp
MLTp
MLMp
MLPp
MLMp
MLTp
MLHp
MLTp
MLYp
MLCp
MLDp
MLEp
MLDp
MLNp
MLQp
MLVp
MLMp
MLAp
MLCp
MLEp
MLGp
MLAp
MLCp
MLVp
MLYp
MLNp
MLNp
MLNp
MLHp
MLRp
MLKp
MLDp
MLFp
MLQp
MLNp
MLFp
MLPp
MLCp
MLSp
MLGp
MLLp
MLNp
MLMp
MLWp
MLTp
MLKp
MLRp
MLIp
MLVp
MLYp
MLPp
MLMp
MLAp
MLTp
MLNp
MLGp
MLWp
MLHp
MLDp
MLSp
MLSp
MLDp
MLKp
MLIp
MLGp
MLTp
MLCp
MLNp
MLVp
MLRp
MLRp
MLYp
MLDp
MLFp
MLSp
MLKp
MLMp
MLGp
MLVp
MLSp
MLAp
MLHp
MLDp
MLEp
MLHp
MLCp
MLYp
MLPp
MLQp
MLEp
MLHp
MLLp
MLDp
MLVp
MLTp
MLPp
MLNp
MLMp
MLFp
MLVp
MLRp
MLCp
MLDp
MLTp
MLRp
MLYp
MLSp
MLCp
MLTp
MLEp
MLVp
MLIp
MLGp
MLNp
MLGp
MLGp
MLSp
MLKp
MLIp
MLQp
MLSp
MLAp
MLEp
MLTp
MLWp
MLLp
MLRp
MLWp
MLWp
MLFp
MLFp
MLIp
MLYp
MLCp
MLPp
MLTp
MLNp
MLMp
MLWp
MLHp
MLSp
MLIp
MLTp
MLWp
MLMp
MLTp
MLIp
MLYp
MLKp
MLHp
MLQp
MLQp
MLLp
MLCp
MLKp
MLKp
MLFp
MLRp
MLCp
MLWp
MLRp
MLIp
MLGp
MLMp
MLPp
MLRp
MLYp
MLYp
MLYp
MLLp
MLEp
MLKp
MLDp
MLQp
MLKp
MLRp
MLRp
MLRp
MLKp
MLAp
MLHp
MLWp
MLQp
MLKp
MLLp
MLQp
MLAp
MLAp
MLPp
MLGp
MLEp
MLFp
MLIp
MLLp
MLEp
MLSp
MLEp
MLNp
MLQp
MLDp
MLTp
MLKp
MLDp
MLYp
MLAp
MLDp
MLSp
MLCp
MLHp
MLSp
MLVp
MLKp
MLDp
MLHp
MLHp
MLEp
MLAp
MLFp
MLCp
MLCp
MLDp
MLTp
MLAp
MLAp
MLKp
MLDp
MLEp
MLWp
MLMp
MLSp
MLDp
MLMp
MLDp
MLIp
MLRp
MLPp
MLTp
MLYp
MLTp
MLRp
MLGp
MLMp
MLGp
MLRp
MLFp
MLNp
MLNp
MLLp
MLCp
MLTp
MLIp
MLTp
MLYp
MLSp
MLCp
MLIp
MLEp
MLSp
MLMp
MLHp
MLHp
MLPp
MLQp
MLHp
MLKp
MLMp
MLYp
MLLp
MLTp
MLQp
MLCp
MLQp
MLKp
MLLp
MLQp
MLNp
MLDp
MLNp
MLLp
MLIp
MLYp
MLPp
MLEp
MLKp
MLGp
MLKp
MLKp
MLSp
MLAp
MLKp
MLCp
MLHp
MLCp
MLDp
MLGp
MLTp
MLKp
MLTp
MLAp
MLEp
MLGp
MLVp
MLQp
MLDp
MLHp
MLTp
MLRp
MLHp
MLGp
MLRp
MLGp
MLHp
MLIp
MLWp
MLAp
MLCp
MLGp
MLHp
MLFp
MLIp
MLRp
MLCp
MLVp
MLSp
MLSp
MLSp
MLGp
MLGp
MLGp
MLGp
MLWp
MLFp
MLFp
MLAp
MLCp
MLFp
MLVp
MLGp
MLHp
MLLp
MLYp
MLWp
MLVp
MLWp
MLVp
MLNp
MLLp
MLVp
MLCp
MLYp
MLGp
MLVp
MLVp
MLNp
MLKp
MLSp
MLGp
MLTp
MLLp
MLMp
MLRp
MLCp
MLWp
MLEp
MLMp
MLVp
MLTp
MLQp
MLKp
MLSp
MLLp
MLEp
MLMp
MLKp
MLCp
MLEp
MLVp
MLWp
MLYp
MLEp
MLRp
MLQp
MLGp
MLVp
MLMp
MLHp
MLLp
MLDp
MLWp
MLTp
MLCp
MLLp
MLFp
MLVp
MLDp
MLRp
MLKp
MLIp
MLNp
MLKp
MLTp
MLHp
MLKp
MLQp
MLGp
MLMp
MLNp
MLIp
MLCp
MLPp
MLEp
MLFp
MLGp
MLQp
MLCp
MLHp
MLDp
MLDp
MLDp
MLTp
MLKp
MLSp
MLKp
MLKp
MLTp
MLNp
MLRp
MLQp
MLNp
MLMp
MLTp
MLPp
MLGp
MLLp
MLDp
MLSp
MLEp
MLAp
MLDp
MLVp
MLNp
MLMp
MLRp
MLPp
MLCp
MLSp
MLDp
MLQp
MLGp
MLIp
MLKp
MLNp
MLDp
MLIp
MLCp
MLWp
MLLp
MLRp
MLHp
MLCp
MLGp
MLSp
MLEp
MLPp
MLWp
MLNp
MLNp
MLSp
MLYp
MLLp
MLCp
MLMp
MLEp
MLAp
MLGp
MLVp
MLMp
MLWp
MLVp
MLEp
MLVp
MLFp
MLYp
MLMp
MLNp
MLFp
MLRp
MLKp
MLLp
MLFp
MLMp
MLDp
MLEp
MLPp
MLAp
MLLp
MLIp
MLKp
MLVp
MLFp
MLEp
MLRp
MLWp
MLMp
MLCp
MLFp
MLKp
MLYp
MLLp
MLMp
MLFp
MLPp
MLWp
MLFp
MLLp
MLGp
MLHp
MLHp
MLVp
MLLp
MLPp
MLPp
MLYp
MLKp
MLDp
MLVp
MLTp
MLAp
MLIp
MLSp
MLTp
MLFp
MLIp
MLSp
MLNp
MLIp
MLHp
MLLp
MLGp
MLEp
MLEp
MLVp
MLHp
MLIp
MLVp
MLCp
MLAp
MLWp
MLKp
MLVp
MLYp
MLTp
MLYp
MLGp
MLIp
MLDp
MLTp