If you do not have an initial alignment between your templates and target sequence, MODELLER can derive one for you, fully automatically. All MODELLER requires is a a PIR file containing the target sequence and the template PDB codes (their sequences are not required — just use a single '*' character — as MODELLER will read these from the PDBs). Use the AutoModel class as per usual, but call the AutoModel.auto_align() method before AutoModel.make(); see the example below. (MODELLER has a variety of other alignment methods which you can use instead for this purpose; see Section 6.16 for more details.)
Please be aware that the single most important factor that determines the quality of a model is the quality of the alignment. If the alignment is incorrect, the model will also be incorrect. For this reason, automated alignment for comparative modeling should not be used unless the sequences are so similar that the calculated alignment is likely to be correct, which usually requires more than 50% sequence identity. Instead, the alignment should be carefully inspected, optimized by hand, and checked by the Alignment.check() command before used in modeling. Moreover, several iterations of alignment and modeling may be necessary in general.
# A sample script for fully automated comparative modeling from modeller import * from modeller.automodel import * # Load the AutoModel class log.verbose() env = Environ() # directories for input atom files env.io.atom_files_directory = ['.', '../atom_files'] a = AutoModel(env, # file with template codes and target sequence alnfile = 'alignment.seg', # PDB codes of the templates knowns = ('5fd1', '1fdn', '1fxd', '1iqz'), # code of the target sequence = '1fdx') a.auto_align() # get an automatic alignment a.make() # do comparative modeling
>P1;1fdx sequence::::::ferredoxin:Peptococcus aerogenes:-1.00:-1.00 AYVINDSCIACGACKPECPVNIIQGSIYAIDADSCIDCGSCASVCPVGAPNPED* >P1;1fdn structureX:1fdn:FIRST:@:55:@:ferredoxin:Clostrodium acidiurici: 1.84:-1.0 * >P1;5fd1 structureX:5fd1:FIRST:@:60:@:ferredoxin:Azotobacter vinelandii: 1.90:0.192 * >P1;1fxd structureX:1fxd:FIRST:@:58:@:ferredoxin:Desolfovibrio gigas: 1.70:-1.0 * >P1;1iqz structureX:1iqz:FIRST:@:60:@:ferredoxin:Bacillus thermoproteolyticus: 2.30:-1.0 *