On 10/16/23 8:15 AM, Isabel Elliott via modeller_usage wrote: > I am trying to model missing residues into my crystal structure > following the method on the website: > https://salilab.org/modeller/wiki/Missing_residues. I am also trying > to restrain the crystal coordinates so that only missing residues move > (and the rest of the original crystal coordinates remain the same) > during refinement using the select_atoms method. However when I compare > my final models produced by my script to my original crystal structure, > I find that the side chains of several amino acids three of four amino > acids away from the missing residues have been flipped. How can I ensure > that only the missing residues are refined and the crystal coordinates > remain the same? I am wondering if I’ve missed something important out > of my script. See my script below.
In your script you have defined the select_loop_atoms method, *not* select_atoms. The script is essentially doing three things:
1. Building an initial model, using your template where the sequence is the same, and internal coordinates where you have gaps or sequence differences. This file is generated in PDB format with a '.ini' extension.
2. Refining a subset of the model atoms, defined by select_atoms (by default all atoms will be refined).
3. Further refinement using Modeller's loop modeling protocol, potentially of a different subset (defined by select_loop_atoms).
If your template is missing some sidechains then they will be filled in in step (1) above - and you'll see that in the .ini file. Otherwise, defining select_atoms, as described in the wiki page at https://salilab.org/modeller/wiki/Missing_residues, is what you need.
Ben Webb, Modeller Caretaker