| The following script optimize.py is
not working for me. It is giving an error in the line: is giving an error in this line: mdl.generate_topology<aln, sequence=code+'ini'> Please suggest. Script env = environ() env.edat.dynamic_sphere = True env.libs.topology.read(file='$(LIB)/top_heav.lib') env.libs.parameters.read(file='$(LIB)/par.lib') code = 'unbound-scrwl' aln = alignment(env) mdl = model(env, file=code) aln.append_model(mdl, align_codes=code, atom_files=code) aln.append_model(mdl, align_codes=code+'-ini', atom_files=code+'.ini') mdl.generate_topology(aln, sequence=code+'-ini') mdl.transfer_xyz(aln) mdl.build(initialize_xyz=False, build_method='INTERNAL_COORDINATES') mdl.write(file=code+'.ini') # Generate the restraints: mdl.restraints.make(aln, restraint_type='stereo', spline_on_site=False) mdl.restraints.write(file=code+'.rsr') mpdf = mdl.energy() mdl.switch_trace(file=code+'.trc') mdl.optimize(optimization_method=1, max_iterations=20, output='REPORT', trace_output=1) mdl.optimize(optimization_method=3, temperature=300, max_iterations=50, output='REPORT', trace_output=1) mdl.optimize(optimization_method=1, max_iterations=20, output='REPORT', trace_output=1) mpdf = mdl.energy() mdl.write(file=code+'.B') Dr. Jayashree Ramana Lecturer, Department of Biotechnology and BioinformaticsJaypee University of Information Technology Waknaghat, Solan Himachal Pradesh-173234 http://www.juit.ac.in/Department/bio/faculty_bio.php?faculty_id=118 |