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.

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.

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.

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.

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.