include rare tautomers in output

User 9ba4ed1585

18-08-2014 09:03:31

Hi,


I'm preparing a set of ligands for docking, I use majormicrospecies to convert smiles into mol2, then standardizer and corina, is there any way to include the rare protonation forms when using majormicro (e.g. only use form that have an abundancy bigger than 10%)? I believe that in some cases the rare tautomer is actually the active drug form,


best


Tomek

ChemAxon d51151248d

19-08-2014 11:46:16

Dear Tomek, 


can you send us an example file so we can investigate your problem?


Thanks, 


Daniel

User 9ba4ed1585

20-08-2014 17:40:36

Ok, here are three example ligands in SMILEs format, they have two major tautomeric forms/protonation states, one of which is dominant, but I would like the cxcalc majormicrospecies two give me every form which has an abundance higher than 20%, not only the major form,


best 


Tomek

ChemAxon d51151248d

21-08-2014 12:37:53

Dear Tomek, 


could please attach the SMILES codes you mentioned?


Thanks, 


Daniel

User 9ba4ed1585

21-08-2014 13:32:24

Sorry, for some strange reason the forum doesnt want to attach my file, where here are the smiles:


Fc1ccccc1NC(=O)CCN1CCC(CC1)Cc1ccccc1 CHEMBL101133


O=C(Nc1ccccc1[N+](=O)[O-])CCN1CCN(CC1)Cc1ccccc1 CHEMBL101122


O=C(Nc1ccccc1CC)CCN1CCN(CC1)Cc1ccccc1 CHEMBL101088


 

ChemAxon d51151248d

22-08-2014 12:18:12

Dear Tomek, 


you can use the pKa plugin in MarvinSketch to generate the microspecies distribution  for your molecule.From here you can see e.g. that molecule #1 has 4, while molecule #2 has 6 microspecies alltogether. You can also see the microspecies with relevant distribution at a given pH.


In cxcalc you can geneate all microspecies with their distribtion at a given pH (default is 7.4) using the 


cxcalc microspeciesdistribution test.mol

command. In this case the microspecies distribution is stored in a MSDISTR[pH=...]property. You can use the


cxcalc -h 

command to get a full list of calculations.


BR, 


Daniel