Dear modeller team,
I am using a variant of the align2d_mult.py script on your homepage.
I want to use a multi-sequence alignment, align several template structures and the query sequence to it without changing the initial multi-sequence alignment. On the webpage it is mentioned that one can use align_block to fix a specific part of the alignment. Just to make sure - by using the script:
#!/usr/bin/python
from modeller import *
log.verbose()#output=1, notes=1, warnings=1, errors=1, memory=0)
env = environ()
env.libs.topology.read(file='$(LIB)/top_heav.lib')
# Read aligned sequences:
aln = alignment(env)
aln.append(file='../alignments/test.pir', align_codes='all')
aln_block = len(aln)
# Read aligned structures:
aln.append(file='../templates/pdb_fit/alignment_test_structures.ali', align_codes='all')
# Structure sensitive variable gap penalty sequence-sequence alignment:
aln.salign(output='', max_gap_length=20,
gap_function=False, # to use structure-dependent gap penalty
alignment_type='PAIRWISE', align_block=aln_block,
feature_weights=(1., 0., 0., 0., 0., 0.), overhang=0,
gap_penalties_1d=(-450, 0),
gap_penalties_2d=(0.35, 1.2, 0.9, 1.2, 0.6, 8.6, 1.2, 0., 0.),
similarity_flag=True)
aln.write(file='../alignments/modeller_multi.ali', alignment_format='PIR')
aln.check()
---------------------------
where "alignment_test_structures.ali" contains an alignment of potential template structures; is each of the sequence in "alignment_test_structures.ali" independently aligned to the alignment given in "test.pir" without changing the "test.pir" alignment?
Thanks in advance!
Kind regards,
Simone