Toward a universal model to calculate the solvation thermodynamics of druglike molecules: the importance of new experimental databases.

We demonstrate that a new free energy functional in the integral equation theory of molecular liquids gives accurate calculations of hydration thermodynamics for druglike molecules. The functional provides an improved description of excluded volume effects by incorporating two free coefficients. When the values of these coefficients are obtained from experimental data for simple organic molecules, the hydration free energies of an external test set of druglike molecules can be calculated with an accuracy of about 1 kcal/mol. The 3D RISM/UC method proposed here is easily implemented using existing computational software and allows in silico screening of the solvation thermodynamics of potential pharmaceutical molecules at significantly lower computational expense than explicit solvent simulations.

[1]  John B. O. Mitchell,et al.  Predicting intrinsic aqueous solubility by a thermodynamic cycle. , 2008, Molecular pharmaceutics.

[2]  Fumio Hirata,et al.  Potential of Mean Force between Two Molecular Ions in a Polar Molecular Solvent: A Study by the Three-Dimensional Reference Interaction Site Model , 1999 .

[3]  A. Bauer-Brandl,et al.  Thermodynamic properties of flufenamic and niflumic acids--specific and non-specific interactions in solution and in crystal lattices, mechanism of solvation, partitioning and distribution. , 2007, Journal of pharmaceutical and biomedical analysis.

[4]  David Chandler,et al.  Free energy functions in the extended RISM approximation , 1985 .

[5]  Maxim V. Fedorov,et al.  Multigrid solver for the reference interaction site model of molecular liquids theory , 2011, J. Comput. Chem..

[6]  Jean-François Truchon,et al.  Predictions of hydration free energies from continuum solvent with solute polarizable models: the SAMPL2 blind challenge , 2010, J. Comput. Aided Mol. Des..

[7]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[8]  Donald G Truhlar,et al.  SM6:  A Density Functional Theory Continuum Solvation Model for Calculating Aqueous Solvation Free Energies of Neutrals, Ions, and Solute-Water Clusters. , 2005, Journal of chemical theory and computation.

[9]  B. Roux,et al.  Solvation Free Energy of Polar and Nonpolar Molecules in Water: An Extended Interaction Site Integral Equation Theory in Three Dimensions , 2000 .

[10]  C. Cramer,et al.  Self-Consistent Reaction Field Model for Aqueous and Nonaqueous Solutions Based on Accurate Polarized Partial Charges. , 2007, Journal of chemical theory and computation.

[11]  Fumio Hirata,et al.  Solvation thermodynamics of protein studied by the 3D-RISM theory , 2004 .

[12]  Maxim V Fedorov,et al.  Accurate calculations of the hydration free energies of druglike molecules using the reference interaction site model. , 2010, The Journal of chemical physics.

[13]  M. Kinoshita,et al.  Morphometric approach to the solvation free energy of complex molecules. , 2006, Physical review letters.

[14]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[15]  F. J. Luque,et al.  Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. , 2000, Chemical reviews.

[16]  Maxim V Fedorov,et al.  An accurate prediction of hydration free energies by combination of molecular integral equations theory with structural descriptors. , 2010, The journal of physical chemistry. B.

[17]  Miguel Jorge,et al.  1-Octanol/Water Partition Coefficients of n-Alkanes from Molecular Simulations of Absolute Solvation Free Energies. , 2009, Journal of chemical theory and computation.

[18]  J. Guthrie,et al.  A blind challenge for computational solvation free energies: introduction and overview. , 2009, The journal of physical chemistry. B.

[19]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[20]  Seiichiro Ten-no,et al.  Comparative study on solvation free energy expressions in reference interaction site model integral equation theory. , 2005, The journal of physical chemistry. B.

[21]  Andriy Kovalenko,et al.  Association thermodynamics and conformational stability of beta-sheet amyloid beta(17-42) oligomers: effects of E22Q (Dutch) mutation and charge neutralization. , 2010, Biophysical journal.

[22]  Andriy Kovalenko,et al.  An MM/3D-RISM approach for ligand binding affinities. , 2010, The journal of physical chemistry. B.

[23]  Fumio Hirata,et al.  Solution of three‐dimensional reference interaction site model and hypernetted chain equations for simple point charge water by modified method of direct inversion in iterative subspace , 1999 .

[24]  Yuichi Harano,et al.  Theoretical study for partial molar volume of amino acids and polypeptides by the three-dimensional reference interaction site model , 2001 .

[25]  A. Bauer-Brandl,et al.  Influence of Position and Size of Substituents on the Mechanism of Partitioning: A Thermodynamic Study on Acetaminophens, Hydroxybenzoic Acids, and Parabens , 2008, AAPS PharmSciTech.

[26]  J. Andrew Grant,et al.  SAMPL2 and continuum modeling , 2010, J. Comput. Aided Mol. Des..

[27]  Maxim V Fedorov,et al.  Combination of RISM and Cheminformatics for Efficient Predictions of Hydration Free Energy of Polyfragment Molecules: Application to a Set of Organic Pollutants. , 2011, Journal of chemical theory and computation.

[28]  F. Hirata,et al.  Theoretical study for volume changes associated with the helix-coil transition of peptides. , 2001, Biopolymers.

[29]  István Kolossváry,et al.  Low Mode Search . An Efficient , Automated Computational Method for Conformational Analysis : Application to Cyclic and Acyclic Alkanes and Cyclic Peptides , 1997 .

[30]  Carlos Simmerling,et al.  Three-dimensional molecular theory of solvation coupled with molecular dynamics in Amber. , 2010, Journal of chemical theory and computation.

[31]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[32]  B. Montgomery Pettitt,et al.  A dielectrically consistent interaction site theory for solvent—electrolyte mixtures , 1992 .

[33]  D. Blankschtein,et al.  Liquid-state theory of hydrocarbon-water systems: application to methane, ethane, and propane , 1992 .

[34]  Norio Yoshida,et al.  Selective ion-binding by protein probed with the 3D-RISM theory. , 2006, Journal of the American Chemical Society.

[35]  Florian Nigsch,et al.  Why Are Some Properties More Difficult To Predict than Others? A Study of QSPR Models of Solubility, Melting Point, and Log P , 2008, J. Chem. Inf. Model..

[36]  A. Bauer-Brandl,et al.  Towards an understanding of the molecular mechanism of solvation of drug molecules: a thermodynamic approach by crystal lattice energy, sublimation, and solubility exemplified by hydroxybenzoic acids. , 2006, Journal of pharmaceutical sciences.

[37]  Yuko Okamoto,et al.  First-Principle Determination of Peptide Conformations in Solvents: Combination of Monte Carlo Simulated Annealing and RISM Theory , 1998 .

[38]  S. Ten-no,et al.  Assessment of free energy expressions in RISM integral equation theory: theoretical prediction of partition coefficients revisited , 2010 .

[39]  David L Mobley,et al.  Predicting small-molecule solvation free energies: an informal blind test for computational chemistry. , 2008, Journal of medicinal chemistry.

[40]  D. Chandler,et al.  Excess electrons in simple fluids. I. General equilibrium theory for classical hard sphere solvents , 1984 .

[41]  A. Bauer-Brandl,et al.  Energetic aspects of diclofenac acid in crystal modifications and in solutions--mechanism of solvation, partitioning and distribution. , 2007, Journal of pharmaceutical sciences.

[42]  A. Bauer-Brandl,et al.  Solvation of drugs as a key for understanding partitioning and passive transport exemplified by NSAIDs. , 2004, Current drug delivery.

[43]  W. L. Jorgensen The Many Roles of Computation in Drug Discovery , 2004, Science.

[44]  Fumio Hirata,et al.  An extended rism equation for molecular polar fluids , 1981 .

[45]  G. Perlovich,et al.  Thermodynamic and structural study of tolfenamic acid polymorphs. , 2009, Journal of pharmaceutical and biomedical analysis.

[46]  Norio Yoshida,et al.  Theoretical study of CO escaping pathway in myoglobin with the 3D-RISM theory. , 2009, Journal of the American Chemical Society.

[47]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[48]  G. Maggiora,et al.  Solvation thermodynamics of polar molecules in aqueous solution by the XRISM method , 1993 .

[49]  Jason Crain,et al.  Improved estimates for hydration free energy obtained by the reference interaction site model , 2007 .

[50]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[51]  Fumio Hirata,et al.  Hydration free energy of hydrophobic solutes studied by a reference interaction site model with a repulsive bridge correction and a thermodynamic perturbation method , 2000 .

[52]  Donald G Truhlar,et al.  Charge Model 4 and Intramolecular Charge Polarization. , 2007, Journal of chemical theory and computation.

[53]  Anthony Nicholls,et al.  The SAMPL2 blind prediction challenge: introduction and overview , 2010, J. Comput. Aided Mol. Des..

[54]  Arieh Ben-Naim,et al.  Standard thermodynamics of transfer. Uses and misuses , 1978 .

[55]  Fumio Hirata,et al.  Molecular Theory of Solvation , 2004 .

[56]  Maxim V Fedorov,et al.  Towards a universal method for calculating hydration free energies: a 3D reference interaction site model with partial molar volume correction , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[57]  Andriy Kovalenko,et al.  Calculation of local water densities in biological systems: a comparison of molecular dynamics simulations and the 3D-RISM-KH molecular theory of solvation. , 2011, The journal of physical chemistry. B.

[58]  Donald G Truhlar,et al.  Performance of SM6, SM8, and SMD on the SAMPL1 test set for the prediction of small-molecule solvation free energies. , 2009, The journal of physical chemistry. B.

[59]  M. Kinoshita Roles of translational motion of water molecules in sustaining life. , 2009, Frontiers in bioscience.

[60]  F. A. Neugebauer,et al.  Electrochemical oxidation and structural changes of 5,6-dihydrobenzo[c]cinnolines , 1996 .

[61]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach , 2000 .

[62]  Fumio Hirata,et al.  Water molecules in a protein cavity detected by a statistical-mechanical theory. , 2005, Journal of the American Chemical Society.

[63]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[64]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[65]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[66]  David Chandler,et al.  Optimized Cluster Expansions for Classical Fluids. II. Theory of Molecular Liquids , 1972 .

[67]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[68]  Jan H. Jensen,et al.  Prediction and rationalization of protein pKa values using QM and QM/MM methods. , 2005, The journal of physical chemistry. A.

[69]  István Kolossváry,et al.  Low‐mode conformational search elucidated: Application to C39H80 and flexible docking of 9‐deazaguanine inhibitors into PNP , 1999 .

[70]  A. Kovalenko,et al.  Microtubule stability studied by three-dimensional molecular theory of solvation. , 2007, Biophysical journal.

[71]  Seiichiro Ten-no,et al.  Free energy of solvation for the reference interaction site model: Critical comparison of expressions , 2001 .

[72]  W. L. Jorgensen,et al.  Prediction of drug solubility from structure. , 2002, Advanced drug delivery reviews.

[73]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[74]  A. D. J. Haymet,et al.  Integral equation theory for uncharged liquids: The Lennard‐Jones fluid and the bridge function , 1995 .

[75]  Masahiro Kinoshita,et al.  Peptide Conformations in Alcohol and Water: Analyses by the Reference Interaction Site Model Theory , 2000 .

[76]  Fumio Hirata,et al.  Ligand mapping on protein surfaces by the 3D-RISM theory: toward computational fragment-based drug design. , 2009, Journal of the American Chemical Society.

[77]  Arieh Ben-Naim,et al.  Solvation thermodynamics of nonionic solutes , 1984 .

[78]  A. Bauer-Brandl,et al.  Extent and mechanism of solvation and partitioning of isomers of substituted benzoic acids: a thermodynamic study in the solid state and in solution. , 2008, Journal of pharmaceutical sciences.