Ali Algarrous wrote:
I'm sorry but it seems I didn't explain myself clearer. What I asked for last time is different than what I'm asking now. This time I want the script to transfer the coordinates for aligned residues only and nothing that is not aligned. So the MSFIEKMIG in my sequance is not going to have any coordinates because it was not aligned with anything from my templae.
That's straightforward - leave out the mdl.build(...) command from your script. transfer_xyz only transfers the known coordinates, so the coordinates for your loop atoms will end up 'undefined' (999).
(On the other hand, if you are calling mdl.build() because you have missing atoms in some of your non-loop residues which you want Modeller to fill in for you, you could leave mdl.build() in to build all coordinates, then select all of your loops using model.loops(), then use selection.unbuild() to reset the coordinates of the loop atoms to 'undefined'.)
Ben Webb, Modeller Caretaker
--
http://www.salilab.org/modeller/
Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage