User 53611d53ae
26-01-2012 20:36:09
Hello ChemAxon,
I am contacting you about a problem that we, at the bioinformatic research group at SRI International, find to be significant. I have spoken with people at ChEBI (Marcus Ennis) and UniProt (Anne Morgat) and they encounter the same hurdles.
When describing the stereochemistry of structures such as sterols, it is customary to draw several hydrogen bonds with stereo information (for example, look at the structure of cholesterol at wikipedia, which has 4 such hydrogens -see http://en.wikipedia.org/wiki/File:Cholesterol.svg
When the Marvin protonation plugin protonates such structures, it removes all of these hydrogens, and represents the stereo information that was defined by them by adding up and down directions to some of the wedges within the ring itself (see below). The structures end up looking so bad that we have stopped using Marvin for protonating any of the compounds that could be effected, and now perform the protonation of these compounds by hand. Similarly, curators at ChEBI manually fix every such compound following the analysis by Marvin.
Needless to say, these hydrogens all have pKa values that fall far beyond the relevant range, and really should not even be considered by the pKa module. Would it be possible to modify the protonation code so that it ignores hydrogens whose pKa falls outside the range of 1-14, and leaves those bonds as-is? This would be a great help to several of your "power users".
Thank you,
Ron Caspi
ChemAxon e08c317633
27-01-2012 10:51:06
rcaspi wrote: |
When the Marvin protonation plugin protonates such structures, it removes all of these hydrogens, and represents the stereo information that was defined by them by adding up and down directions to some of the wedges within the ring itself (see below). The structures end up looking so bad that we have stopped using Marvin for protonating any of the compounds that could be effected, and now perform the protonation of these compounds by hand. Similarly, curators at ChEBI manually fix every such compound following the analysis by Marvin.
|
Thanks for the feedback, we will consider implementing it.
rcaspi wrote: |
Needless to say, these hydrogens all have pKa values that fall far beyond the relevant range, and really should not even be considered by the pKa module. Would it be possible to modify the protonation code so that it ignores hydrogens whose pKa falls outside the range of 1-14, and leaves those bonds as-is? This would be a great help to several of your "power users".
|
pKa plugin already has options to set the pKa range to consider. The options are:
Min basic pKa: widens the calculation range because weak bases will have lower pKa values than the default -10.
Max acidic pKa: widens the calculation range because weak acids will have higher pKa values than the default 20.
For more see http://www.chemaxon.com/marvin/help/calculations/protonation.html#pka
Zsolt
User 851ac690a0
06-02-2012 14:32:33
Dear Ron,
Thanks for reporting this problem!
We will fix it asap.
Jozsi
User 53611d53ae
07-02-2012 00:43:21
Dear Jozsi,
Great. We are looking forward!
Take care,
Ron
User 53611d53ae
11-06-2012 20:57:53
Hi Jozsi,
looks like a fix didn't make it into beta 5.10. Do you have an idea when we could expect a fix for this problem?
Thanks,
Ron
User 851ac690a0
13-06-2012 08:56:21
Hi,
The fix will be available in 5.10.1. It will be released in July.
Jozsi
User 53611d53ae
31-05-2013 18:19:29
Jozsi,
Looks like this problem is still around in 6.0.
Ron
User 851ac690a0
31-05-2013 18:47:15
Hi,
I have tried the 6.0 version on the Chemaxon home site. See the attached figure. It seems to work now.
A new option was introduced "Keep explicit hydrogens". But unfortunately we managed to hide it from users.
Please switch on the "Keep explicit hydrogens" option. Because according to the default setting it is switched off (we will change the default).
Jozsi
User 53611d53ae
31-05-2013 21:52:36
Right, it works! I did not notice that option. Thanks!