Computing intrinsic aqueous solubility of crystalline organic molecules

Research output: Contribution to conferenceSpeech

Abstract

Accurate computational methods to predict the solubility of crystalline organic molecules in water are highly sought after in many fields of the biomolecular sciences. For example, predictions of solubility are used in the agrochemical and pharmaceutical industries to assess the environmental fate of potential pollutants and the bioavailability of de novo designed drugs, respectively. While a number of statistical Quantitative Structure-Property Relationship (QSPR) solubility models have been built by many groups, including ourselves, development of a method based on chemical theory and molecular simulation has proved challenging. Successful construction of a fully theoretical method would bring both greater insight and systematic improvability, with physically meaningful component effects becoming quantifiable.
Here we show that reasonably accurate values for the intrinsic aqueous solubility of druglike organic molecules are obtainable from theory. We combine sublimation free energies calculated using the crystal lattice minimisation program DMACRYS with hydration free energies calculated using the 3D Reference Interaction Site Model (3DRISM) of the Integral Equation Theory of Molecular Liquids (IET). Intrinsic solubilities of 25 diverse crystalline druglike molecules are computed by our method, obtaining R=0.85 and RMSE =1.45 log10 S units. This is considerably better than results from implicit continuum solvent models such as the polarizable continuum model (PCM). We can fully computationally describe the thermodynamics of the transfer of a druglike molecule from the crystalline solid phase via the gas phase to dilute aqueous solution. Although the 3D RISM free energy functional that we use contains some parameters, our approach is not parameterized against experimentally measured solubilities. On-going work to improve the accuracy of the calculated sublimation thermodynamic parameters will also be discussed.

Conference

ConferenceInternational Conference on the Properties of Water and Steam (ICPWS16)
CountryUnited Kingdom
CityLondon
Period1/09/135/09/13

Fingerprint

solubility
molecules
free energy
sublimation
continuums
bioavailability
thermodynamics
crystal lattices
contaminants
hydration
integral equations
solid phases
drugs
industries
vapor phases
aqueous solutions
optimization
liquids
predictions
water

Keywords

  • solubility
  • solvation
  • druglike
  • pharmaceutical
  • in silico screening
  • physico-chemical properties

Cite this

Palmer, D., & Fedorov, M. (2013). Computing intrinsic aqueous solubility of crystalline organic molecules. International Conference on the Properties of Water and Steam (ICPWS16), London, United Kingdom.
Palmer, David ; Fedorov, Maxim. / Computing intrinsic aqueous solubility of crystalline organic molecules. International Conference on the Properties of Water and Steam (ICPWS16), London, United Kingdom.
@conference{e35d159f20d44a28aaaf2c2bce28ae56,
title = "Computing intrinsic aqueous solubility of crystalline organic molecules",
abstract = "Accurate computational methods to predict the solubility of crystalline organic molecules in water are highly sought after in many fields of the biomolecular sciences. For example, predictions of solubility are used in the agrochemical and pharmaceutical industries to assess the environmental fate of potential pollutants and the bioavailability of de novo designed drugs, respectively. While a number of statistical Quantitative Structure-Property Relationship (QSPR) solubility models have been built by many groups, including ourselves, development of a method based on chemical theory and molecular simulation has proved challenging. Successful construction of a fully theoretical method would bring both greater insight and systematic improvability, with physically meaningful component effects becoming quantifiable.Here we show that reasonably accurate values for the intrinsic aqueous solubility of druglike organic molecules are obtainable from theory. We combine sublimation free energies calculated using the crystal lattice minimisation program DMACRYS with hydration free energies calculated using the 3D Reference Interaction Site Model (3DRISM) of the Integral Equation Theory of Molecular Liquids (IET). Intrinsic solubilities of 25 diverse crystalline druglike molecules are computed by our method, obtaining R=0.85 and RMSE =1.45 log10 S units. This is considerably better than results from implicit continuum solvent models such as the polarizable continuum model (PCM). We can fully computationally describe the thermodynamics of the transfer of a druglike molecule from the crystalline solid phase via the gas phase to dilute aqueous solution. Although the 3D RISM free energy functional that we use contains some parameters, our approach is not parameterized against experimentally measured solubilities. On-going work to improve the accuracy of the calculated sublimation thermodynamic parameters will also be discussed.",
keywords = "solubility, solvation, druglike, pharmaceutical, in silico screening, physico-chemical properties",
author = "David Palmer and Maxim Fedorov",
year = "2013",
month = "9",
day = "4",
language = "English",
note = "International Conference on the Properties of Water and Steam (ICPWS16) ; Conference date: 01-09-2013 Through 05-09-2013",

}

Palmer, D & Fedorov, M 2013, 'Computing intrinsic aqueous solubility of crystalline organic molecules' International Conference on the Properties of Water and Steam (ICPWS16), London, United Kingdom, 1/09/13 - 5/09/13, .

Computing intrinsic aqueous solubility of crystalline organic molecules. / Palmer, David; Fedorov, Maxim.

2013. International Conference on the Properties of Water and Steam (ICPWS16), London, United Kingdom.

Research output: Contribution to conferenceSpeech

TY - CONF

T1 - Computing intrinsic aqueous solubility of crystalline organic molecules

AU - Palmer, David

AU - Fedorov, Maxim

PY - 2013/9/4

Y1 - 2013/9/4

N2 - Accurate computational methods to predict the solubility of crystalline organic molecules in water are highly sought after in many fields of the biomolecular sciences. For example, predictions of solubility are used in the agrochemical and pharmaceutical industries to assess the environmental fate of potential pollutants and the bioavailability of de novo designed drugs, respectively. While a number of statistical Quantitative Structure-Property Relationship (QSPR) solubility models have been built by many groups, including ourselves, development of a method based on chemical theory and molecular simulation has proved challenging. Successful construction of a fully theoretical method would bring both greater insight and systematic improvability, with physically meaningful component effects becoming quantifiable.Here we show that reasonably accurate values for the intrinsic aqueous solubility of druglike organic molecules are obtainable from theory. We combine sublimation free energies calculated using the crystal lattice minimisation program DMACRYS with hydration free energies calculated using the 3D Reference Interaction Site Model (3DRISM) of the Integral Equation Theory of Molecular Liquids (IET). Intrinsic solubilities of 25 diverse crystalline druglike molecules are computed by our method, obtaining R=0.85 and RMSE =1.45 log10 S units. This is considerably better than results from implicit continuum solvent models such as the polarizable continuum model (PCM). We can fully computationally describe the thermodynamics of the transfer of a druglike molecule from the crystalline solid phase via the gas phase to dilute aqueous solution. Although the 3D RISM free energy functional that we use contains some parameters, our approach is not parameterized against experimentally measured solubilities. On-going work to improve the accuracy of the calculated sublimation thermodynamic parameters will also be discussed.

AB - Accurate computational methods to predict the solubility of crystalline organic molecules in water are highly sought after in many fields of the biomolecular sciences. For example, predictions of solubility are used in the agrochemical and pharmaceutical industries to assess the environmental fate of potential pollutants and the bioavailability of de novo designed drugs, respectively. While a number of statistical Quantitative Structure-Property Relationship (QSPR) solubility models have been built by many groups, including ourselves, development of a method based on chemical theory and molecular simulation has proved challenging. Successful construction of a fully theoretical method would bring both greater insight and systematic improvability, with physically meaningful component effects becoming quantifiable.Here we show that reasonably accurate values for the intrinsic aqueous solubility of druglike organic molecules are obtainable from theory. We combine sublimation free energies calculated using the crystal lattice minimisation program DMACRYS with hydration free energies calculated using the 3D Reference Interaction Site Model (3DRISM) of the Integral Equation Theory of Molecular Liquids (IET). Intrinsic solubilities of 25 diverse crystalline druglike molecules are computed by our method, obtaining R=0.85 and RMSE =1.45 log10 S units. This is considerably better than results from implicit continuum solvent models such as the polarizable continuum model (PCM). We can fully computationally describe the thermodynamics of the transfer of a druglike molecule from the crystalline solid phase via the gas phase to dilute aqueous solution. Although the 3D RISM free energy functional that we use contains some parameters, our approach is not parameterized against experimentally measured solubilities. On-going work to improve the accuracy of the calculated sublimation thermodynamic parameters will also be discussed.

KW - solubility

KW - solvation

KW - druglike

KW - pharmaceutical

KW - in silico screening

KW - physico-chemical properties

UR - http://icpws16.imeche.org/

M3 - Speech

ER -

Palmer D, Fedorov M. Computing intrinsic aqueous solubility of crystalline organic molecules. 2013. International Conference on the Properties of Water and Steam (ICPWS16), London, United Kingdom.