How accurate are continuum solvation models for drug-like molecules?

Jacob Kongsted, Pär Söderhjelm, Ulf Ryde

Research output: Contribution to journalArticlepeer-review

65 Citations (SciVal)


We have estimated the hydration free energy for 20 neutral drug-like molecules, as well as for three series of 6-11 inhibitors to avidin, factor Xa, and galectin-3 with four different continuum solvent approaches (the polarised continuum method the Langevin dipole method, the finite-difference solution of the Poisson equation, and the generalised Born method), and several variants of each, giving in total 24 different methods. All four types of methods have been thoroughly calibrated for a number of experimentally known small organic molecules with a mean absolute deviation (MAD) of 1-6 kJ/mol for neutral molecules and 4-30 kJ/mol for ions. However, for the drug-like molecules, the accuracy seems to be appreciably worse. The reason for this is that drug-like molecules are more polar than small organic molecules and that the uncertainty of the methods is proportional to the size of the solvation energy. Therefore, the accuracy of continuum solvation methods should be discussed in relative, rather than absolute, terms. In fact, the mean unsigned relative deviations of the best solvation methods, 0.09 for neutral and 0.05 for ionic molecules, correspond to 2-20 kJ/mol absolute error for the drug-like molecules in this investigation, or 2-3,000 in terms of binding constants. Fortunately, the accuracy of all methods can be improved if only relative energies within a series of inhibitors are considered, especially if all of them have the same net charge. Then, all except two methods give MADs of 2-5 kJ/mol (corresponding to an uncertainty of a factor of 2-7 in the binding constant). Interestingly, the generalised Born methods typically give better results than the Poison-Boltzmann methods.
Original languageEnglish
Pages (from-to)395-409
JournalJournal of Computer-Aided Molecular Design
Issue number7
Publication statusPublished - 2009

Bibliographical note

The information about affiliations in this record was updated in December 2015.
The record was previously connected to the following departments: Theoretical Chemistry (S) (011001039)

Subject classification (UKÄ)

  • Theoretical Chemistry


  • Free energy of hydration
  • Continuum methods
  • PCM
  • Langevin dipoles
  • Generalised Born
  • FDPB
  • Solvation energy


Dive into the research topics of 'How accurate are continuum solvation models for drug-like molecules?'. Together they form a unique fingerprint.

Cite this