Re-reading what I wrote, here is clarification: Protein X reads residue indices from pdb, starts at 103, ends at 173,
Your results are for sequence 1-70, You want to see in the structure why in your results segment 46-53 is unique. select 46-53 results in no atoms selected, because you are off with the numbering. Instead you need to calulate 46+103=149 (and perhaps +1, never sure about that) and 53+103=156 and re-issue the command select 149-156. Then you want to know, if the loop in positions 125-132 in the structure, which you pick by cursor on screen is somehow identified in your calculation based on the 70aa sequence. Now 125-103=22 (?) and 132-103=129 (?) is needed to convert the screen structure indexing into your tables. Ad so forth. Re C How Rasmo'sl renumber x worked was that if there were gaps in the readed indices e.g. ATOM 11 NH2 ARG A 103 -10.684 12.480 58.277 1.00 25.62 RADA N ATOM 12 N PRO A 108 -14.360 5.552 63.412 1.00 22.41 RADA N it renumbered all sequentially so the virtual result was something like ATOM 11 NH2 ARG A 1 -10.684 12.480 58.277 1.00 25.62 RADA N ATOM 12 N PRO A 2 -14.360 5.552 63.412 1.00 22.41 RADA N which of course was a problem, because the you have no clue where you are (bug in the rasmol sort of) The desired outcome of renumber 1 (if residue 103 is the first in the molecule) would be something like ATOM 11 NH2 ARG A 1 -10.684 12.480 58.277 1.00 25.62 RADA N ATOM 12 N PRO A 6 -14.360 5.552 63.412 1.00 22.41 RADA N 103->1 (104)->2 (105)->3 (106)->4 (107)->5 108-> 6 Petr Dr. Petr Pancoska Department of Pathology SUNY Stony Brook, NY 11794 phone: (631)-444-3030 ****************************************************************************** This e- mail message, including any attachments, is for the sole use of the intended recipient(s) and may contain confidential and privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not the intended recipient, please contact the sender by e-mail and destroy all copies of the original. ****************************************************************************** "Miguel" <[EMAIL PROTECTED]> Sent by: To jmol-users-admin@ jmol-users@lists.sourceforge.net lists.sourceforge cc .net Subject Re: [Jmol-users] renumbering in 12/21/2004 12:14 Jmol PM Please respond to [EMAIL PROTECTED] sourceforge.net Petr wrote: > Trying "renumber x" in the Rasmol Scripts console results in Script > command not implemented: renumber. True, that is not implemented. A year ago when I was working on the scripting engine I did not understand how it worked, and noone has asked for it. > Cannot find any replacement functionality in the docs. This is relatively > important function for transferring sequence analysis into 3D picture. Please explain ... I do not understand what you mean. > PDB files unfortunately sometimes contain numbering that starts not where > one would like to so it is painful to re-calculate everything to the new > indexing. > renumber was an easy remedy for that - any suggestions how to do that? I now understand that the residue sequence numbers are often not sequential, and that insertion codes add to the complexity. I can probably implement 'renumber' ... but the idea of changing the raw data makes me nervous ... I fear it will have unforseen repercussions. So first ... I do not understand what the value of 'renumber' is. It seems to me that it would break all of the scripts associated with that .pdb file. Q: What is the benefit of using 'renumber'? What are you trying to accomplish? Miguel ------------------------------------------------------- SF email is sponsored by - The IT Product Guide Read honest & candid reviews on hundreds of IT Products from real users. Discover which products truly live up to the hype. Start reading now. http://productguide.itmanagersjournal.com/ _______________________________________________ Jmol-users mailing list Jmol-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/jmol-users ------------------------------------------------------- SF email is sponsored by - The IT Product Guide Read honest & candid reviews on hundreds of IT Products from real users. Discover which products truly live up to the hype. Start reading now. http://productguide.itmanagersjournal.com/ _______________________________________________ Jmol-users mailing list Jmol-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/jmol-users