Dear All,
I hope you all know that there is script available in modeller "mutate_model.py". To generate a mutant models and then optimize the mutated residue following refinement in sidechain alone but I want to make this process (optimization/refinement) including the backbone/near residues instead of restricting only with sidechain.
So, Does anyone know how to modify the default script
I hope the modification needs on the following lines (139 to 153)
#only optimize the selected residue (in first pass, just atoms in selected #residue, in second pass, include nonbonded neighboring atoms) #set up the mutate residue selection segment #s = selection(mdl1.chains[chain].residues[respos]) s = selection(mdl1.chains[chain].residues[respos])
mdl1.restraints.unpick_all() mdl1.restraints.pick(s)
s.energy()
s.randomize_xyz(deviation=4.0)
mdl1.env.edat.nonbonded_sel_atoms=2 optimize(s,sched)
Thanking you in advance