[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[modeller_usage] problem of build_profile
- To: Modeller Caretaker <>,
- Subject: [modeller_usage] problem of build_profile
- From: Bishwajit Das <>
- Date: Tue, 8 Mar 2011 00:51:49 +0530
When i run this program showing an error
from modeller import *
log.verbose()
env = environ()
#-- Prepare the input files
#-- Read in the sequence database
sdb = sequence_db(env)
sdb.read(seq_database_file='pdb95.fsa', seq_database_format='FASTA',
chains_list='ALL', minmax_db_seq_len=(1, 40000), clean_sequences=True)
#-- Write the sequence database in binary form
sdb.write(seq_database_file='pdb95.bin', seq_database_format='BINARY',
chains_list='ALL')
#-- Now, read in the binary database
sdb.read(seq_database_file='pdb95.bin', seq_database_format='BINARY',
chains_list='ALL')
#-- Read in the target sequence/alignment
aln = alignment(env)
aln.append(file='toxin.ali', alignment_format='PIR', align_codes='ALL')
#-- Convert the input sequence/alignment into
# profile format
prf = aln.to_profile()
#-- Scan sequence database to pick up homologous sequences
prf.build(sdb, matrix_offset=-450, rr_file='${LIB}/blosum62.sim.mat',
gap_penalties_1d=(-500, -50), n_prof_iterations=5,
check_profile=False, max_aln_evalue=0.01, gaps_in_target=False)
#-- Write out the profile
prf.write(file='buildprofile.prf', profile_format='TEXT')
#-- Convert the profile back to alignment format
aln = prf.to_alignment()
#-- Write out the alignment file
aln.write(file='buildprofile.ali', alignment_format='PIR')
could not open the file pdb_95.fsa
where i get this file??
--
BISHWAJIT DAS
+91-9433237920
KOLKATA