Re: [modeller_usage] Filling in Missing atoms for a docking complex
On Wed, Apr 1, 2009 at 10:55 AM, Jack Shultz jshultz@hydrogenathome.org wrote: > Hello, > > I am new to modeller and I have not gone through the tutorial yet. I > have submitted a couple of FASTA format files to mod web. I am would > like to fix a problem with my docking complexes I get from autodock. > Some of them are missing atoms in a given residue. It was suggested I > either use Deep-View or Modeller to insert missing atoms. Can modeller > indeed achieve this? Are there any existing scripts to perform this? > > -- > Jack > > http://www.facebook.com/home.php#/profile.php?id=832713248 > http://hydrogenathome.org >
Jack Shultz wrote: > I am new to modeller and I have not gone through the tutorial yet. I > have submitted a couple of FASTA format files to mod web. I am would > like to fix a problem with my docking complexes I get from autodock. > Some of them are missing atoms in a given residue. It was suggested I > either use Deep-View or Modeller to insert missing atoms. Can modeller > indeed achieve this? Are there any existing scripts to perform this?
Certainly - that is straightforward, and in fact Modeller will do it automatically. Just build a homology model using your existing model as a template (the alignment will be simply 1:1 since the sequences are the same). Use the regular modeling script at http://salilab.org/modeller/9v6/manual/node16.html
If you want to keep the original conformation of part of the system, use the script at http://salilab.org/modeller/9v6/manual/node23.html instead.
You can also use the same technique to fill in entire missing *residues*; see http://salilab.org/modeller/wiki/Missing%20residues for an example.
Ben Webb, Modeller Caretaker
Good o know about those options. I was just skimming through section 2 in the manual over lunch and I came across this example script. Would this also work? If so, are there any advantages to using model-segment.py?
2.2.9 Using your own initial model Normally, automodel generates an initial model by transferring coordinates from the templates. However, if you have a prepared PDB file containing an initial model, you can have automodel use this instead with the inifile keyword, as in the example below. (This automatically sets automodel.generate method to generate.read xyz for you, which is necessary for this to work.) Example: examples/automodel/model-myini.py # Modeling using a provided initial structure file (inifile) 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, alnfile = ’alignment.ali’, # alignment filename knowns = ’5fd1’, # codes of the templates sequence = ’1fdx’, # code of the target inifile = ’my-initial.pdb’) # use ’my’ initial structure a.starting_model= 1 # index of the first model a.ending_model = 1 # index of the last model # (determines how many models to calculate) a.make() # do homology modeling
On Wed, Apr 1, 2009 at 12:15 PM, Modeller Caretaker modeller-care@salilab.org wrote: > Jack Shultz wrote: >> I am new to modeller and I have not gone through the tutorial yet. I >> have submitted a couple of FASTA format files to mod web. I am would >> like to fix a problem with my docking complexes I get from autodock. >> Some of them are missing atoms in a given residue. It was suggested I >> either use Deep-View or Modeller to insert missing atoms. Can modeller >> indeed achieve this? Are there any existing scripts to perform this? > > Certainly - that is straightforward, and in fact Modeller will do it > automatically. Just build a homology model using your existing model as > a template (the alignment will be simply 1:1 since the sequences are the > same). Use the regular modeling script at > http://salilab.org/modeller/9v6/manual/node16.html > > If you want to keep the original conformation of part of the system, use > the script at http://salilab.org/modeller/9v6/manual/node23.html instead. > > You can also use the same technique to fill in entire missing > *residues*; see http://salilab.org/modeller/wiki/Missing%20residues for > an example. > > Ben Webb, Modeller Caretaker > -- > modeller-care@salilab.org http://www.salilab.org/modeller/ > Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage >
Jack Shultz wrote: > Good o know about those options. I was just skimming through section 2 > in the manual over lunch and I came across this example script. Would > this also work? If so, are there any advantages to using > model-segment.py? > > > [2.2.9 Using your own initial model]
Sure, that would also work, but in that case you'd have to first generate a PDB file as the initial model. In this case, since your target and template are 100% identical, the default initial model will match the template exactly anyway, so there seems little point in using your own initial model. The initial model is simply the starting guess for the optimizer - it doesn't otherwise affect the way in which models are constructed.
In homology modeling, the output models will not be exactly the same as the input (sidechain dihedrals will vary, for example). If this bothers you, you can fix the positions of some or all of the atoms you already know. That's the advantage of using model-segment.py.
Ben Webb, Modeller Caretaker
participants (2)
-
Jack Shultz
-
Modeller Caretaker