Approaches for machine learning intermolecular interaction energies and application to energy components from symmetry adapted perturbation theory.

Accurate prediction of intermolecular interaction energies is a fundamental challenge in electronic structure theory due to their subtle character and small magnitudes relative to total molecular energies. Symmetry adapted perturbation theory (SAPT) provides rigorous quantum mechanical means for computing such quantities directly and accurately, but for a computational cost of at least O(N5), where N is the number of atoms. Here, we report machine learned models of SAPT components with a computational cost that scales asymptotically linearly, O(N). We use modified multi-target Behler-Parrinello neural networks and specialized intermolecular symmetry functions to address the idiosyncrasies of the intermolecular problem, achieving 1.2 kcal mol-1 mean absolute errors on a test set of hydrogen bound complexes including structural data extracted from the Cambridge Structural Database and Protein Data Bank, spanning an interaction energy range of 20 kcal mol-1. Additionally, we recover accurate predictions of the physically meaningful SAPT component energies, of which dispersion and induction/polarization were the easiest to predict and electrostatics and exchange-repulsion are the most difficult.

Collaborating to improve the use of free-energy and other quantitative methods in drug discovery.

In May and August, 2016, several pharmaceutical companies convened to discuss and compare experiences with Free Energy Perturbation (FEP). This unusual synchronization of interest was prompted by Schrödinger's FEP+ implementation and offered the opportunity to share fresh studies with FEP and enable broader discussions on the topic. This article summarizes key conclusions of the meetings, including a path forward of actions for this group to aid the accelerated evaluation, application and development of free energy and related quantitative, structure-based design methods.

To measure is to know: an approach to CADD performance metrics.

With the financial and productivity challenges currently facing the pharmaceutical industry, there is constant pressure to justify resources and improve efficiency. With process-driven activities, understanding the contribution of these resources is reasonably straightforward. By contrast, measuring the contribution of knowledge workers is less obvious. Here, we present an impact-oriented approach to assessing the performance of an industrial computer-assisted drug design group. We discuss how these metrics are used to understand and optimize resource allocation in support of drug discovery programs.

Trends in kinase selectivity: insights for target class-focused library screening.

A kinome-wide selectivity screen of >20000 compounds with a rich representation of many structural classes has been completed. Analysis of the selectivity patterns for each class shows that a broad spectrum of structural scaffolds can achieve specificity for many kinase families. Kinase selectivity and potency are inversely correlated, a trend that is also found in a large set of kinase functional data. Although selective and nonselective compounds are mostly similar in their physicochemical characteristics, we identify specific features that are present more frequently in compounds that bind to many kinases. Our results support a scaffold-oriented approach for building compound collections to screen kinase targets.

A simple strategy for mitigating the effect of data variability on the identification of active chemotypes from high-throughput screening data.

Among the several goals of a high-throughput screening campaign is the identification of as many active chemotypes as possible for further evaluation. Often, however, the number of concentration response curves (e.g., IC(50)s or K(i)s) that can be collected following a primary screen is limited by practical constraints such as protein supply, screening workload, and so forth. One possible approach to this dilemma is to cluster the hits from the primary screen and sample only a few compounds from each cluster. This introduces the question as to how many compounds must be selected from a cluster to ensure that an active compound is identified, if it exists at all. This article seeks to address this question using a Monte Carlo simulation in which the dependence of the success of sampling is directly linked to screening data variability. Furthermore, the authors demonstrate that the use of replicated compounds in the screening collection can easily assess this variability and provide a priori guidance to the screener and chemist as to the extent of sampling required to maximize chemotype identification during the triage process. The individual steps of the Monte Carlo simulation provide insight into the correspondence between the percentage inhibition and eventual IC(50) curves.

Tryptamine and homotryptamine-based sulfonamides as potent and selective inhibitors of 15-lipoxygenase.

A series of inhibitors of mammalian 15-lipoxygenase based on tryptamine and homotryptamine scaffolds is described. Compounds with aryl substituents at C-2 of the indole core of tryptamine and homotryptamine sulfonamides (e.g., 37a-p) proved to be potent inhibitors of the isolated enzyme. Selected compounds also demonstrated desirable inhibition selectivities over isozymes 5- and P-12-LO.