This is done by using the AllHModel class rather than AutoModel. Otherwise, operation is identical. Note that the AllHModel class automatically turns on env.io.hydrogen for you and selects the all-atom topology and radii files.
from modeller import * from modeller.automodel import * log.verbose() env = Environ() env.io.atom_files_directory = ['.', '../atom_files'] a = AllHModel(env, alnfile='alignment.ali', knowns='5fd1', sequence='1fdx') a.starting_model = a.ending_model = 4 a.make()