<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Nadir,<div><span class="Apple-tab-span" style="white-space:pre"> </span>I think it's fairly well known that hydrogen bonds are continuous phenomena and not discrete on/off entities as presented by Chimera (and most programs of its ilk). That's why Chimera allows you to loosen its H-bond criteria to explore interactions beyond the default cutoffs. Similarly, weak H-bond interactions can occur between chemical entities other than the ones examined by the Chimera H-bond code, <i>e.g.</i> carbon donors or pi-cloud acceptors. When Chimera's H-bond code was written, Chimera had no capability to assess the energetics of H-bond interactions and therefore simply uses the geometric criteria outlined in the Mills paper cited in Chimera's findhbond documentation to locate obvious H-bond interactions between well known donors and acceptors. C-H hydrogen bonds were excluded in the Mills paper due to their low occurrence and the fact that their geometry was frequently distorted to accommodate stronger H-bond interactions.</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>Chimera isn't really trying to reinvent the wheel here. I think if one needs sophisticated analysis of all types of hydrogen bonding interactions then a dedicated QM or MM program should be used, not a visualization program. Nonetheless, if geometric criteria similar to those in the Mills paper were determined for some of these weaker interactions, I would consider adding them to the code.</div><div><br></div><div>--Eric</div><div><br></div><div><div><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> Eric Pettersen</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> UCSF Computer Graphics Lab</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> <a href="http://www.cgl.ucsf.edu">http://www.cgl.ucsf.edu</a></font></div><br class="Apple-interchange-newline"></span></div><div><div>On Mar 26, 2013, at 7:25 AM, Nadir T. Mrabet wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Hi Elaine,<br><br>I gave a reference that should be read before initiating a fruitful discussion.<br>I give a second one, which is an (critical) essay and a synthetic summary :<br>Desiraju G.R. (2011) "A Bond by Any Other Name", Angew. Chem. Int. Ed. 50: 52-59.<br><br>The hydrogen bond is designated X-H°°°Y-Z (4 atoms; there have been so many articles where only the 3 atoms of a H-bond, C-H°°°O (3) are considered).<br>X must be more electronegative than H.<br>The H-bond acceptor, Y, is electron rich (either lone electron pair, or pi-bonded pair in Y-Z).<br>The evidence for H-bond may be experimental or theoretical, but better be both.<br><br>Now, the C atom can be recognized as an H-bond donor (if it fulfills the required criteria).<br>Aromatic rings become H-bond acceptors.<br><br>For the purpose of optimizing energy, the angle, X-H°°°Y, must be as close as possible to 180 ° (not new!), but there is no constraint imposed on angle H°°°Y-Z (a point on which I disagree).<br><br>Quantum-mechanical calculation is not adamant, but should help (provided one uses an appropriate force field).<br><br>Therefore, even if force fields attribute partial atomic charges (Note that Amber does not agree "qualitatively" with others on some atoms) they usually mention only NH and OH (and sometimes S) as H-bond donors and N-Z or O-Z, as acceptors. C-H is not considered as an H-bond donor; the pi-electron cloud is not recognized as acceptor.<br><br>The tough job resides in the choice of an empirical force field that best correlates with experimental data.<br>Redefining H-bond donors and acceptors and their partial charges is also to be done in the molecular modeling program(s).<br><br>I cannot be exhaustive in a matter of minutes.<br>I just wish to start the discussion and see whether (or not) if it will impact on Chimera.<br><br>Best regards,<br><br>Nadir<br><br>Pr. Nadir T. Mrabet<br>Structural & Molecular Biochemistry<br>N-gere - INSERM U-954<br>University of Lorraine, Nancy<br>School of Sciences and Technologies<br>& School of Medicine<br>9, Avenue de la Foret de Haye, BP 184<br>54505 Vandoeuvre-les-Nancy Cedex<br>France<br>Phone: +33 (0)3.83.68.32.73<br>Fax: +33 (0)3.83.68.32.79<br>E-mail: Nadir.Mrabet <at> <a href="http://univ-lorraine.fr">univ-lorraine.fr</a><br>Cell.: +33 (0)6.11.35.69.09<br><br>LEGAL NOTICE<br>Unless expressly stated otherwise, this message is confidential and may<br>be privileged. It is intended for the addressee(s) only.<br>Access to this E-mail by anyone else is unauthorized.<br>If you are not an addressee, any disclosure or copying of the contents<br>of this E-mail, or any action taken (or not taken) in reliance on it,<br>is unauthorized and may be unlawful.<br>If you are not an addressee, please inform the sender immediately.<br><br>On 25/03/2013 23:13, Elaine Meng wrote:<br><blockquote type="cite">Hi Nadir,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Maybe it would help if you explained the change in definition, why you think this definition is an improvement over currently used measures, and whether the criteria are things that can be measured or calculated by molecular graphics & analysis programs. For example, if it requires quantum-mechanical calculations, that would be beyond the capabilities of several programs.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Elaine<br></blockquote><blockquote type="cite">----------<br></blockquote><blockquote type="cite">Elaine C. Meng, Ph.D.<br></blockquote><blockquote type="cite">UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br></blockquote><blockquote type="cite">Department of Pharmaceutical Chemistry<br></blockquote><blockquote type="cite">University of California, San Francisco<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Mar 23, 2013, at 12:42 PM, Nadir Mrabet wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><blockquote type="cite">Hi Elaine,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">I addressed the question to CCP4BB but il did not seem to interest anyone.<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">I state it again :<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">I wonder how much of an impact should there be on structure refinement, validation and, last but not least, molecular modeling, if we take into account the recent (2011) re-definition of the hydrogen bond made by the IUPAC.<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Ref:<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">"Elangannan Arunan, Gautam R. Desiraju, Roger A. Klein,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Joanna Sadlej, Steve Scheiner, Ibon Alkorta, David C. Clary,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Robert H. Crabtree, Joseph J. Dannenberg, Pavel Hobza,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Henrik G. Kjaergaard, Anthony C. Legon, Benedetta Mennucci,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">and David J. Nesbitt<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Pure Appl. Chem., Vol. 83, No. 8, pp. 1619–1636, 2011.<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">doi:10.1351/PAC-REP-10-01-01<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">© 2011 IUPAC, Publication date (Web): 8 July 2011<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Defining the hydrogen bond: An account (IUPAC Technical Report)*"<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Tough job, isn't ?<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">My guess is that Chimera should be involved.<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Nadir Mrabet<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Dear Elaine,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"> Many thanks for the suggestion and reply. I will do as said<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">and will mail if any problem persists.<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">On Sat, Mar 23, 2013 at 11:41 PM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>> wrote:<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Dear Amit,<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Yes, there is a FindHBond tool (in menu under Tools... Structure Analysis) and command "findhbond" with the same capabilities. Click the Help button on the tool to see its full instructions, or enter command "help findhbond" to see the command instructions.<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">There is an example of using FindHBond in the "Structure Analysis and Comparison" tutorial (see Chimera menu, Help... Tutorials).<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">You can also use Help... Search Documentation in the Chimera menu, for example to search for "hydrogen bonds" ... this will give several links into the user manual including the things I mentioned above.<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">I hope this helps,<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Elaine<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">-----<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Elaine C. Meng, Ph.D.<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Department of Pharmaceutical Chemistry<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">University of California, San Francisco<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">On Mar 23, 2013, at 11:02 AM, Amit Jaiswal wrote:<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Dear Chimera Users,<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"> I am new to chimera and I have some<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">protein protein interaction results to be analysed. Most of<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">the results are from HADDOCK server which gives a<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">protein protein complex. So, is it possible to visualize<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">the hydrogen bonds formed between the proteins using<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">chimera? If yes, how? I am eagerly waiting for the answer.<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">-- <br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Yours Sincerely,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">~~~~~~~~~~~~~~~~~~~~~~~~<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Amit Jaiswal,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Department of Bioinformatics,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">School of Life Sciences,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Pondicherry Central University,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Kalapet, Pondicherry,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Puducherry - 605 014.<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">~~~~~~~~~~~~~~~~~~~~~~~~~<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">_______________________________________________<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Chimera-users mailing list<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">_______________________________________________<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Chimera-users mailing list<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br></blockquote></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">_______________________________________________<br></blockquote><blockquote type="cite">Chimera-users mailing list<br></blockquote><blockquote type="cite"><a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br></blockquote><blockquote type="cite"><a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br></blockquote><blockquote type="cite"><br></blockquote><br>_______________________________________________<br>Chimera-users mailing list<br><a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br>http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users<br><br></div></blockquote></div><br></div><br><br></body></html>