1
$\begingroup$

How to predict the half-maximal inhibitory concentration (IC50) of my ligands.

I'm currently using Schrodinger suite software and I don't know how to predict the IC50?

$\endgroup$
1
  • $\begingroup$ It is usually helpful to do a little research with a search engine or to give more information about what exactly you need. See e.g. schrodinger.com/kb/936 $\endgroup$ Sep 23, 2020 at 17:45

1 Answer 1

2
$\begingroup$

I don't think it is possible just from the docking scores. Docking scores are designed to give you an estimate of the binding affinity (entalphy), but they do not take entropic contributions into account very well. The best you can do are molecular dynamics simulations to calculate the free energy of binding. Which can be computationally quite costly e.g. with GROMACS or CHARMM.

You can estimate the IC50 with molecular dynamics simulations as well, but that is not straightforward. E.g. with Umbrella simulations. An example for a potassium channel and a very simple ligand can be found here.

Another approach would be to build a predictive model. If you have a library of compounds at hand you could use the chemical structures and e.g. the docking scores to estimate the IC50 based on statistics or machine learning. But biology is very messy in general, though if you get a representative dataset with very similar compounds to the one that you have in your dataset then you might be successful. You have to be sure that the binding site and pose are correct though.

$\endgroup$

Not the answer you're looking for? Browse other questions tagged or ask your own question.