Technical Support Forum Index
Technical Support Forum
Access ChemAxon scientists and developers here. For registration and login issues contact website support.

Support Ticket System is replacing forum

This forum was converted into a searchable archive. You cannot add posts here any more. For support please use our new Ticket System.

Create your first ticket
Methology details of pKa calculator?
To watch this topic for replies  Register (enables digests) or give email address:
This topic is locked: you cannot edit posts or make replies.
Display posts from previous:   
    View previous topic :: View next topic    
Author Message
Ed

Joined: 01 Feb 2013
Posts: 5

View user's profile

Back to top
Link to postPosted: Wed Aug 17, 2016 2:26 pmPost subject: Methology details of pKa calculator? Reply with quote

Hello to all,

 

After having some communications with the ChemAxon support team about details regarding the pKa algorithm, it would be a good idea to make the questions public.  So, here goes - hope you can answer them!

-          Which partial-charge algorithm did you use?  It wouldn’t be Gasteiger-Marsili by any chance, if not an in-house ChemAxon algorithm?

-          Same for polarisabililty – how did you calculate it?

-          What are “structure-specific” increments and how do you define ionisable sites?  A ChemAxon poster mentioned some simple definitions of acidic and basic groups, but are these actually used in the regression model?  Surely if it’s site-specific then you’d need to take into account the functional group (as well as surrounding atoms to account for steric effects)?

-          What sort of compounds did you train the regression model on? 

-          How does this relate to cxtrain when a “local model” is trained?  This is very important as we have several unique chemical series which one won’t find in public datasets.

 

In addition, do you have any published (or private) examples of cxtrain-pKa being used?  My experience on our compounds showed little benefit, but I also think that I didn’t use cxtrain properly, yet I've found no evidence in the literature (yet) that shows cxtrain being used for pKa.  Seeing as it’s a regression model based on descriptors as well-respected as partial charges, I find it somewhat hard to believe that training did not work in our case...

 

Thanks in advance,

Ed.

Daniel
ChemAxon personnel
Joined: 22 Jan 2012
Posts: 315

View user's profile

Back to top
Link to postPosted: Thu Aug 18, 2016 2:30 pmPost subject: Reply with quote

Hi Ed,

Here are my answers for your questions:

  1. We use the Gasteiger partial charge calculation method, which we extended with some parameters. We had to heavily extend the original set of atom types of the Gasteiger method, and provide parameters for them. For example we can handle molecules with metal atoms, or ionic molecules, which the original method was uncapable of. We also had to tweak some of the calculations for delocalized structures. 

  2. We calculate polarizability using the methods in our reference, basically.

    https://docs.chemaxon.com/display/docs/Polarizability+Plugin      

    In this case we also had to extend the basic set of atom types, and in some cases partial charges have to included in the calculation.

  3.  Structure-specific increments are increments that are attributed to different structure-specific properties that should be taken into account during calculation, for example intra-molecular H-bonds.

    We define ionizable sites by the equilibrium between the conjugated acid and its base, so by the classical conjugated acid-base pair definition. Yes, we include/use in the regression model those simple acidic, basic groups.
     
  4.  We trained our regression model on literature data e.g. Perrin's book or drug compendiums.
     
  5.  When you train with an experimental dataset, ionizable sites are identified and a suitable regression equation is chosen for calculating pKa.
     
  6.  We only have some individual examples here, which shows how cxtrain is used, but not a full publication:

    https://docs.chemaxon.com/display/docs/Training+the+pKa+Plugin

  7.  Well cxtrain should work for your molecules as well, but without the settings and your data I can't tell more. 

I hope this helps, 

Daniel

 

 

This topic is locked: you cannot edit posts or make replies.
Page 1 of 1


To watch this topic for replies   Register (enables digests) or give email address  
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum
You cannot attach files in this forum
You cannot download files in this forum