The influence of hydrogen bonding on partition coefficients.

This Perspective explores how consideration of hydrogen bonding can be used to both predict and better understand partition coefficients. It is shown how polarity of both compounds and substructures can be estimated from measured alkane/water partition coefficients. When polarity is defined in this manner, hydrogen bond donors are typically less polar than hydrogen bond acceptors. Analysis of alkane/water partition coefficients in conjunction with molecular electrostatic potential calculations suggests that aromatic chloro substituents may be less lipophilic than is generally believed and that some of the effect of chloro-substitution stems from making the aromatic π-cloud less available to hydrogen bond donors. Relationships between polarity and calculated hydrogen bond basicity are derived for aromatic nitrogen and carbonyl oxygen. Aligned hydrogen bond acceptors appear to present special challenges for prediction of alkane/water partition coefficients and this may reflect 'frustration' of solvation resulting from overlapping hydration spheres. It is also shown how calculated hydrogen bond basicity can be used to model the effect of aromatic aza-substitution on octanol/water partition coefficients.

ClogP(alk): a method for predicting alkane/water partition coefficient.

Alkane/water partition coefficients (P(alk)) are less familiar to the molecular design community than their 1-octanol/water equivalents and access to both data and prediction tools is much more limited. A method for predicting alkane/water partition coefficient from molecular structure is introduced. The basis for the ClogP(alk) model is the strong (R² = 0.987) relationship between alkane/water partition coefficient and molecular surface area (MSA) that was observed for saturated hydrocarbons. The model treats a molecule as a perturbation of a saturated hydrocarbon molecule with the same MSA and uses increments defined for functional groups to quantify the extent to which logP(alk) is perturbed by the introduction each functional group. Interactions between functional groups, such as intramolecular hydrogen bonds are also parameterized within a perturbation framework. The functional groups and interactions between them are specified substructurally in a transparent and reproducible manner using SMARTS notation. The ClogP(alk) model was parameterized using data measured for structurally prototypical compounds that dominate the literature on alkane/water partition coefficients and then validated using an external test set of 100 alkane/water logP measurements, the majority of which were for drugs.

Automated molecule editing in molecular design.

The ability to modify chemical structures in an automated and controlled manner is useful in molecular design. This Perspective introduces the MUDO molecule editor and shows how automated molecule editing can be used to standardize structures, enumerate tautomeric and ionization states, identify matched molecular pairs. Unlike its predecessor Leatherface, MUDO can also process 3D structures and this capability can be used to link non-covalently docked ligands to proteins.

Discovery of novel Trypanosoma cruzi glyceraldehyde-3-phosphate dehydrogenase inhibitors.

Based on its essential role in the life cycle of Trypanosoma cruzi, the glycolytic enzyme glyceraldehyde 3-phosphate dehydrogenase (GAPDH) has been considered a promising target for the development of novel chemotherapeutic agents for the treatment of Chagas' disease. In the course of our research program to discover novel inhibitors of this trypanosomatid enzyme, we have explored a combination of structure and ligand-based virtual screening techniques as a complementary approach to a biochemical screening of natural products using a standard biochemical assay. Seven natural products, including anacardic acids, flavonoid derivatives, and one glucosylxanthone were identified as novel inhibitors of T. cruzi GAPDH. Promiscuous inhibition induced by nonspecific aggregation has been discarded as specific inhibition was not reversed or affected in all cases in the presence of Triton X-100, demonstrating the ability of the assay to find authentic inhibitors of the enzyme. The structural diversity of this series of promising natural products is of special interest in drug design, and should therefore be useful in future medicinal chemistry efforts aimed at the development of new GAPDH inhibitors having increased potency.

Hydrogen Bond Basicity Prediction for Medicinal Chemistry Design.

Hydrogen bonding is discussed in the context of medicinal chemistry design. Minimized molecular electrostatic potential (Vmin) is shown to be an effective predictor of hydrogen bond basicity (pKBHX), and predictive models are presented for a number of hydrogen bond acceptor types relevant to medicinal chemistry. The problems posed by the presence of nonequivalent hydrogen bond acceptor sites in molecular structures are addressed by using nonlinear regression to fit measured pKBHX to calculated Vmin. Predictions are made for hydrogen bond basicity of fluorine in situations where relevant experimental measurements are not available. It is shown how predicted pKBHX can be used to provide insight into the nature of bioisosterism and to profile heterocycles. Examples of pKBHX prediction for molecular structures with multiple, nonequivalent hydrogen bond acceptors are presented.

Integration of methods in cheminformatics and biocalorimetry for the design of trypanosomatid enzyme inhibitors.

BACKGROUND: The enzyme gapdh, which acts in the glycolytic pathway, is seen as a potential target for pharmaceutical intervention of chagas disease. RESULTS: Herein, we report the discovery of new Trypanosoma cruzi GAPDH (TcGAPDH) inhibitors from target- and ligand-based virtual screening protocols using isothermal titration calorimetry (ITC) and molecular dynamics. Molecular dynamics simulations were used to gain insight on the binding poses of newly identified inhibitors acting at the TcGAPDH substrate (G3P) site. CONCLUSION: Nequimed125, the most potent inhibitor to act upon TcGAPDH so far, which sits on the G3P site without any contact with the co-factor (NAD(+)) site, underpins the result obtained by ITC that it is a G3P-competitive inhibitor. Molecular dynamics simulation provides biding poses of TcGAPDH inhibitors that correlate with mechanisms of inhibition observed by ITC. Overall, a new class of dihydroindole compounds that act upon TcGAPDH through a competitive mechanism of inhibition as proven by ITC measurements also kills T. cruzi.