<div dir="ltr"><div class="gmail_extra">Hi,</div><div class="gmail_extra"><br></div><div class="gmail_extra">Just one other question if I could trouble you? Is it possible to have two colors of H-bond in the per-frame script. That way I could highlight bonds of particular importance (red), while the ones that maintain a steady structure are yellow. That would help draw the users eye and save wordy annotations. I tried the following to no avail. I suspect that maybe this is not possible? If you could offer any guidance, I would be much obliged. :)</div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra">###########</div><div class="gmail_extra">select :1-158; hbonds line 5 color gold ; ~select<br></div><div class="gmail_extra">select :64,80,76,50,21; hbonds selrestrict both retainCurrent true line 5 color red ; ~select</div><div class="gmail_extra">############</div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, May 17, 2015 at 3:13 AM, George Green <span dir="ltr"><<a href="mailto:soyo.green@gmail.com" target="_blank">soyo.green@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Hi Elaine,<div><br></div><div>Thanks for your detailed reply. Only just saw your email. I will test and check. I think your tips should sort my problem out!</div><div><br></div><div>Many thanks!</div></div><div class=""><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 14, 2015 at 5:18 PM, Elaine Meng <span dir="ltr"><<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi George,<br>
I see at least one definite problem and suspect an additional problem.  It is impossible to tell for sure without having your data, which you should usually include with such a question to allow your problem to be reproduced so we don’t have to waste time guessing.<br>
<br>
(1) The definite problem is that your “hbonds” command isn't using the selection in any way.  It simply acts on all atoms, which is fine if that’s what you want it to do.   But based on your description, it is not actually working on all atoms, so there must be some misunderstanding or incomplete description of what’s happening.  To limit “hbonds” to the selection, you’d need the “selRestrict” option.  See “hbonds” documentation for the possible values of that option:<br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/hbonds.html" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/hbonds.html</a>><br>
<br>
It’s possible that hbonds is working on those waters but that they aren’t displayed due to suspected problem #2, below.  To force undisplayed items to be displayed if they are involved in H-bonds, you’d use the hbonds option “reveal true” … but that might accrete more and more displayed residues throughout your trajectory, so I’m not sure you would want it.  “intermodel” and “intramodel” are both default true so you don’t have to put them in your command.<br>
<br>
Also, you don’t need to display the same stuff at every frame since presumably it was displayed in the previous frame.  You would just display it once at the beginning before running your per-frame script.<br>
<br>
(2) The suspected secondary problem is that “:159-162” may not actually specify the waters, because if they don’t have a chain ID, Chimera automatically processes them into chain “water”.  Thus they may not be displayed by your first “disp” command that list residue numbers, but they will be displayed by the second one, “disp :WAT”.  Your select command may not be selecting them (try typing it in the command line to see) and in that case you’d need something like:<br>
<br>
select :1-158:159-162.water<br>
- OR  -<br>
select :1-158:159-162.*<br>
- OR -<br>
select :1-158:wat<br>
(where the latter would select all residues named WAT, so you’d only use that if there weren’t other WAT residues you didn’t want to show)<br>
<br>
 The automatic processing into chain “water” was put into Chimera because some data had protein residues and water residues with the same residue numbers, and the added chain ID provided a handle for specifying them independently.<br>
<br>
I hope this helps,<br>
Elaine<br>
----------<br>
Elaine C. Meng, Ph.D.<br>
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<div><div><br>
> On May 14, 2015, at 7:42 AM, George Green <<a href="mailto:soyo.green@gmail.com" target="_blank">soyo.green@gmail.com</a>> wrote:<br>
><br>
> Hi,<br>
> I'm trying to get a per-frame script going which shows the inter-molecular H-bonds within a protein binding site. It is working in that all H-bonds show apart from water-water h-bonds. Is there any fix for this at all? My code is below:<br>
><br>
> Water molecules are residues :159-162<br>
><br>
> I would be very grateful for any assistance. :)<br>
> sg<br>
><br>
> ################################################<br>
> ####################<br>
> # per-frame H-bond<br>
> repr bs :158<br>
> color hot pink :158<br>
> color byhet :158<br>
> #<br>
> disp :60,37,85,112,109,35,158,84,60,40,38,114,69,90,58,61-68,56,21,45,52,101,15,92,40,105,116,24,42,54,40,24,10,97,120,159,160,161,162<br>
> disp :WAT<br>
> #<br>
> select :1-162; hbonds intramodel true intermodel true line 5 color gold ; ~select<br>
> select :90<br>
><br>
> ################################################<br>
<br>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div>