Hydrophobic Molecular Similarity from MST Fractional Contributions to the Octanol/water Partition Coefficient

Summary The use of a recently proposed hydrophobic similarity index for the alignment of molecules and the prediction of their differences in biological activity is described. The hydrophobic similarity index exploits atomic contributions to the octanol/water transfer free energy, which are evaluated by means of the fractional partitioning scheme developed within the framework of the Miertus-Scrocco-Tomasi continuum model. Those contributions are used to define global and local measures of hydrophobic similarity. The suitability of this computational strategy is examined for two series of compounds (ACAT inhibitors and 5-HT₃ receptor agonists), which are aligned to maximize the global hydrophobic similarity using a Monte Carlo-simulated protocol. Indeed, the concept of local hydrophobic similarity is used to explore structure–activity relationships in a series of COX-2 inhibitors. Inspection of the 3D distribution of hydrophobic/hydrophilic contributions in the aligned molecules is valuable to identify regions of very similar hydrophobicity, which can define pharmacophoric recognition patterns. Moreover, low similar regions permit to identify structural elements that modulate the differences in activity between molecules. Finally, the quantitative relationships found between the pharmacological activity and the hydrophobic similarity index points out that not only the global hydrophobicity, but its 3D distribution, is important to gain insight into the activity of molecules. J.M.M. and S.P. have contributed equally to this study.     

Energy decomposition in molecular complexes: implications for the treatment of polarization in molecular simulations

This study examines the contribution of electrostatic and polarization to the interaction energy in a variety of molecular complexes. The results obtained from the Kitaura-Morokuma (KM) energy decomposition analysis at the HF/6-31G(d) level indicate that, for intermolecular distances around the equilibrium geometries, the polarization energy can be determined as the addition of the polarization energies of interacting blocks, as the mixed polarization term is typically negligible. Comparison of KM and QM/MM results shows that the electrostatic energy determined in the KM method is underestimated (in absolute value) by QM/MM methods. The reason of such underestimation can be attributed to the simplified representation of treating the interaction between overlapping charge distribution by the interaction of a QM molecule with a set of point charges. Nevertheless, the polarization energies calculated by KM and QM/MM methods are in close agreement. Finally, a consistent, automated strategy to derive charge distributions that include implicitly polarization effects in pairwise, additive force fields is presented. The strategy relies in the simultaneous fitting of electrostatic and polarization energies computed by placing a suitable perturbing particle at selected points around the molecule. The suitability of these charges to describe molecular interactions is discussed.     

Accuracy of free energies of hydration for organic molecules from 6-31g*-derived partial charges

Absolute free energies of hydration have been computed for 13 diverse organic molecules using partial charges derived from ab initio 6-31G* wave functions. Both Mulliken charges and charges fit to the electrostatic potential surface (EPS) were considered in conjunction with OPLS Lennard–Jones parameters for the organic molecules and the TIP4P model of water. Monte Carlo simulations with statistical perturbation theory yielded relative free energies of hydration. These were converted to absolute quantities through perturbations to reference molecules for which absolute free energies of hydration had been obtained previously in TIP4P water. The average errors in the computed absolute free energies of hydration are 1.1 kcal/mol for the 6-31G* EPS charges and 4.0 kcal/mol for the Mulliken charges. For the EPS charges, the largest individual errors are under 2 kcal/mol except for acetamide, in which case the error is 3.7 kcal/mol. The hydrogen bonding between the organic solutes and water has also been characterized. © John Wiley & Sons, Inc.