From owner-chemistry \\at// ccl.net Mon Mar 5 18:32:00 2007 From: "=?ISO-8859-1?Q?Tom=E1s_Pe=F1a_Ruiz?= truiz++ujaen.es" <owner-chemistry|a|server.ccl.net> To: CCL Subject: CCL: Concerning NBO bond orders Message-Id: <-33735-070301110610-24044-z7sW5S1I8EsmP06e4T61SQ|a|server.ccl.net> X-Original-From: =?ISO-8859-1?Q?Tom=E1s_Pe=F1a_Ruiz?= <truiz_-_ujaen.es> Content-Type: multipart/alternative; boundary="------------060707080203000104070305" Date: Thu, 01 Mar 2007 16:09:06 +0100 MIME-Version: 1.0 Sent to CCL by: =?ISO-8859-1?Q?Tom=E1s_Pe=F1a_Ruiz?= [truiz(_)ujaen.es] This is a multi-part message in MIME format. --------------060707080203000104070305 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 8bit Hallo!! I'm studying a chemical reaction and trying to compare the modification of the bond orders among reactants, molecular complex and transition state. NBO provides a series of bond orders calculated under different assumptions. Which kind of bond orders could be more appropiate, Wiberg's indexes, atom-atom overlap NAO bond order, NLMO bond orders? Thanks Tom�s Konrad Hinsen hinsen!^!cnrs-orleans.fr wrote: >Sent to CCL by: "Konrad Hinsen" [hinsen]~[cnrs-orleans.fr] >On 28.02.2007, at 07:07, Pradipta Bandyopadhyay pradipta:iitg.ernet.in wrote: > > > >>Applying GNM to a protein (whose structure is a bit open), I found that in >>some regions B-factors obtained from calculation is different from >>experimental one (for the rest other regions these are close). This >>probably means using one gamma may not be enough. Does anyone know if >>there is any simple way to use more than one gamma in GNM (any other >>suggestion would be appreciated also)? >> >> > >It is certainly possible to use different gammas for the different interacting pairs, but I am not >aware of anyone having done so in a GNM. > >If you are reasonably sure that the experimental B factors are reliable, then my first suggestion >would be to try an elastic network model, i.e. a model taking into account the directions of atomic >displacements. This increases the number of coordinates in your system by a factor of three, but >that difference matters only for huge proteins. You gain in two points: > > 1) The description with three coordinates per C-alpha atom is in itself more precise and thus >more realistic. > > 2) The GNM potential is physically unrealistic in not being invariant under rotations. This >means that it penalizes global rotations (which you probably don't care about), but also rotational >domain motions which you may well have in your protein. > >There is also a wider choice of potential models that have been published for elastic network >models. The most widely used variety of elastic network models uses a force constant that depends >on the interatomic distance in the input configuration, nearby atom pairs interacting stronger than >more distant ones. In this family we find, in historical order: > > 1) A exp(-r^2) decrease with distance: > K Hinsen > Analysis of domain motions by approximate normal mode calculations > Proteins 33:417-429 (1998) > http://dirac.cnrs-orleans.fr/plone/publications/all-publications/H_1998 > > 2) A more complex function obtained from fitting to the Amber94 all-atom potential: > K Hinsen, A J Petrescu, S Dellerue, M C Bellissent-Funel, and G R Kneller > Harmonicity in slow protein dynamics > Chem. Phys. 261(1+2):25-37 (2000) > http://dirac.cnrs-orleans.fr/plone/publications/all-publications/Hinsen2000a > > 3) A step function: constant up to a cutoff distance, then zero (like for GNM): > A. R. Atilgan, S. R. Durell, R. L. Jernigan, M. C. Demirel, O. Keskin, and I. Bahar > Anisotropy of Fluctuation Dynamics of Proteins with an Elastic Network Model > Biophys. J. 80, 505515 (2001) > >The best model is probably the one described in 2). You can find a ready-to-use implementation in >the Molecular Modelling Toolkit (http://dirac.cnrs-orleans.fr/MMTK, it's called CalphaForceField >there). There is also a Web server that performs various normal-mode based calculations (including >B factors) using this model for any PDB file you submit; the address is > > http://www.bioinfo.no/tools/normalmodes > >I remember having seen a paper that modifies the elastic network model in choosing the force >constants in a more protein-specific way, but I don't have the reference at hand. > >Konrad. >-- >--------------------------------------------------------------------- >Konrad Hinsen >Centre de Biophysique Molculaire, CNRS Orlans >Synchrotron Soleil - Division Expriences >Saint Aubin - BP 48 >91192 Gif sur Yvette Cedex, France >Tel. +33-1 69 35 97 15 >E-Mail: hinsen||cnrs-orleans.fr >---------------------------------------------------------------------> > > > > > > --------------060707080203000104070305 Content-Type: text/html; charset=us-ascii Content-Transfer-Encoding: 7bit <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> <html> <head> <meta http-equiv="Content-Type" content="text/html;charset=ISO-8859-1"> <title></title> </head> <body text="#000000" bgcolor="#ffffff"> Hallo!!<br> <br> I'm studying a chemical reaction and trying to compare the modification of the bond orders among reactants, molecular complex and transition state. NBO provides a series of bond orders calculated under different assumptions. Which kind of bond orders could be more appropiate, Wiberg's indexes, atom-atom overlap NAO bond order, NLMO bond orders?<br> <br> Thanks<br> <br> Tomás<br> <br> Konrad Hinsen hinsen!^!cnrs-orleans.fr wrote:<br> <blockquote type="cite" cite="mid-id%233er-33712-070301064015-17279-hH5AkmdZw1b344wN%2FCpgZQ+/-server.ccl.net"> <pre wrap="">Sent to CCL by: "Konrad Hinsen" [hinsen]~[cnrs-orleans.fr] On 28.02.2007, at 07:07, Pradipta Bandyopadhyay pradipta:iitg.ernet.in wrote: </pre> <blockquote type="cite"> <pre wrap="">Applying GNM to a protein (whose structure is a bit open), I found that in some regions B-factors obtained from calculation is different from experimental one (for the rest other regions these are close). This probably means using one gamma may not be enough. Does anyone know if there is any simple way to use more than one gamma in GNM (any other suggestion would be appreciated also)? </pre> </blockquote> <pre wrap=""><!----> It is certainly possible to use different gammas for the different interacting pairs, but I am not aware of anyone having done so in a GNM. If you are reasonably sure that the experimental B factors are reliable, then my first suggestion would be to try an elastic network model, i.e. a model taking into account the directions of atomic displacements. This increases the number of coordinates in your system by a factor of three, but that difference matters only for huge proteins. You gain in two points: 1) The description with three coordinates per C-alpha atom is in itself more precise and thus more realistic. 2) The GNM potential is physically unrealistic in not being invariant under rotations. This means that it penalizes global rotations (which you probably don't care about), but also rotational domain motions which you may well have in your protein. There is also a wider choice of potential models that have been published for elastic network models. The most widely used variety of elastic network models uses a force constant that depends on the interatomic distance in the input configuration, nearby atom pairs interacting stronger than more distant ones. In this family we find, in historical order: 1) A exp(-r^2) decrease with distance: K Hinsen Analysis of domain motions by approximate normal mode calculations Proteins 33:417-429 (1998) <a class="moz-txt-link-freetext" href="http://dirac.cnrs-orleans.fr/plone/publications/all-publications/H_1998">http://dirac.cnrs-orleans.fr/plone/publications/all-publications/H_1998</a> 2) A more complex function obtained from fitting to the Amber94 all-atom potential: K Hinsen, A J Petrescu, S Dellerue, M C Bellissent-Funel, and G R Kneller Harmonicity in slow protein dynamics Chem. Phys. 261(1+2):25-37 (2000) <a class="moz-txt-link-freetext" href="http://dirac.cnrs-orleans.fr/plone/publications/all-publications/Hinsen2000a">http://dirac.cnrs-orleans.fr/plone/publications/all-publications/Hinsen2000a</a> 3) A step function: constant up to a cutoff distance, then zero (like for GNM): A. R. Atilgan, S. R. Durell, R. L. Jernigan, M. C. Demirel, O. Keskin, and I. Bahar Anisotropy of Fluctuation Dynamics of Proteins with an Elastic Network Model Biophys. J. 80, 505515 (2001) The best model is probably the one described in 2). You can find a ready-to-use implementation in the Molecular Modelling Toolkit (<a class="moz-txt-link-freetext" href="http://dirac.cnrs-orleans.fr/MMTK">http://dirac.cnrs-orleans.fr/MMTK</a>, it's called CalphaForceField there). There is also a Web server that performs various normal-mode based calculations (including B factors) using this model for any PDB file you submit; the address is <a class="moz-txt-link-freetext" href="http://www.bioinfo.no/tools/normalmodes">http://www.bioinfo.no/tools/normalmodes</a> I remember having seen a paper that modifies the elastic network model in choosing the force constants in a more protein-specific way, but I don't have the reference at hand. Konrad. -- --------------------------------------------------------------------- Konrad Hinsen Centre de Biophysique Molculaire, CNRS Orlans Synchrotron Soleil - Division Expriences Saint Aubin - BP 48 91192 Gif sur Yvette Cedex, France Tel. +33-1 69 35 97 15 E-Mail: hinsen||cnrs-orleans.fr ---------------------------------------------------------------------E-mail to subscribers: <a class="moz-txt-link-abbreviated" href="mailto:CHEMISTRY+/-ccl.net">CHEMISTRY+/-ccl.net</a> or use: <a class="moz-txt-link-freetext" href="http://www.ccl.net/cgi-bin/ccl/send_ccl_message">http://www.ccl.net/cgi-bin/ccl/send_ccl_message</a> E-mail to administrators: <a class="moz-txt-link-abbreviated" href="mailto:CHEMISTRY-REQUEST+/-ccl.net">CHEMISTRY-REQUEST+/-ccl.net</a> or use <a class="moz-txt-link-freetext" href="http://www.ccl.net/cgi-bin/ccl/send_ccl_message">http://www.ccl.net/cgi-bin/ccl/send_ccl_message</a<a class="moz-txt-link-freetext" href="http://www.ccl.net/chemistry/sub_unsub.shtml">http://www.ccl.net/chemistry/sub_unsub.shtml</a> Before posting, check wait time at: <a class="moz-txt-link-freetext" href="http://www.ccl.net">http://www.ccl.net</a> Job: <a class="moz-txt-link-freetext" href="http://www.ccl.net/jobs">http://www.ccl.net/jobs</a> Conferences: <a class="moz-txt-link-freetext" href="http://server.ccl.net/chemistry/announcements/conferences/">http://server.ccl.net/chemistry/announcements/conferences/</a> Search Messages: <a class="moz-txt-link-freetext" href="http://www.ccl.net/htdig">http://www.ccl.net/htdig</a> (login: ccl, Password: search)<a class="moz-txt-link-freetext" href="http://www.ccl.net/spammers.txt">http://www.ccl.net/spammers.txt</a> RTFI: <a class="moz-txt-link-freetext" href="http://www.ccl.net/chemistry/aboutccl/instructions/">http://www.ccl.net/chemistry/aboutccl/instructions/</a</pre> </blockquote> </body> </html> --------------060707080203000104070305--