openmolecules.org

 
Home » DataWarrior » Functionality » Globularity (how do exisitng DW shape/ linearity calculations compare to globularity score?)
Re: Globularity [message #1264 is a reply to message #1262] Wed, 24 March 2021 17:13 Go to previous messageGo to previous message
thomas is currently offline  thomas
Messages: 715
Registered: June 2014
Senior Member
Hi Patrick,

the shape index is based on the 2D-graph of the molecule and looks at the those two atoms that have the longest distance in terms of bonds between them. This length and the atom count of the molecule are used to calculate the shape index.

The globularity is calculated from atoms coordinates in conformers. It is basically doing a singular value decomposition to rotate the molecule in the coordinate system such that its width (size in x-direction) is maximized. Then it is rotated around the x-axis until its height is maximized. After that globularity is calculated as the molecule size in z-direction devided by the size in x-direction.

Shape index and globularity have a slight negative correlation, but one can hardly used as substitute for the other. I have just added a globularity calculation to DataWarrior (dev version). The results correlate (0.83) with the Nature paper's glob value. This is not a perfect correlation. I use up to 128 non-energy minimized conformers and average the result. If I minimize conformers with MMFF94s+ or reduce their number and therefore overpopulate low energy conformers, then the correlation gets worse. Probably, the conformers used in the paper were random conformers with no bias on low energy and most likely not forcefield minimized.

Hope this helps,

Thomas

Report message to a moderator

[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: reaction enumeration reagent connection
Next Topic: missing entry tautobase
Goto Forum:
  

[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Nov 22 21:20:48 CET 2024

Total time taken to generate the page: 0.04198 seconds
openmolecules.org