Chemical space sampling by different scoring functions and crystal structures

Leveraging Data Fusion Strategies in Multireceptor Lead Optimization MM/GBSA End-Point Methods

Accurate and efficient affinity calculations are critical to enhancing the contribution of in silico modeling during the lead optimization phase of a drug discovery campaign. Here, we present a large-scale study of the efficacy of data fusion strategies to leverage results from end-point MM/GBSA calculations in multiple receptors to identify potent inhibitors among an ensemble of congeneric ligands. The retrospective analysis of 13 congeneric ligand series curated from publicly available data across seven biological targets demonstrates that in 90% of the individual receptor structures MM/GBSA scores successfully identify subsets of inhibitors that are more potent than a random selection, and data fusion strategies that combine MM/GBSA scores from each of the receptors significantly increase the robustness of the predictions. Among nine different data fusion metrics based on consensus scores or receptor rankings, the SumZScore (i.e., converting MM/GBSA scores into standardized Z-Scores within a receptor and computing the sum of the Z-Scores for a given ligand across the ensemble of receptors) is found to be a robust and physically meaningful metric for combining results across multiple receptors. Perhaps most surprisingly, even with relatively low to modest overall correlations between SumZScore and experimental binding affinities, SumZScore tends to reliably prioritize subsets of inhibitors that are at least as potent as those that are prioritized from a "best" single receptor identified from known compounds within the congeneric series.

Incorporation of protein flexibility and conformational energy penalties in docking screens to improve ligand discovery

Proteins fluctuate between alternative conformations, which presents a challenge for ligand discovery because such flexibility is difficult to treat computationally owing to problems with conformational sampling and energy weighting. Here we describe a flexible docking method that samples and weights protein conformations using experimentally derived conformations as a guide. The crystallographically refined occupancies of these conformations, which are observable in an apo receptor structure, define energy penalties for docking. In a large prospective library screen, we identified new ligands that target specific receptor conformations of a cavity in cytochrome c peroxidase, and we confirm both ligand pose and associated receptor conformation predictions by crystallography. The inclusion of receptor flexibility led to ligands with new chemotypes and physical properties. By exploiting experimental measures of loop and side-chain flexibility, this method can be extended to the discovery of new ligands for hundreds of targets in the Protein Data Bank for which similar experimental information is available.

Receptor-ligand molecular docking

Docking methodology aims to predict the experimental binding modes and affinities of small molecules within the binding site of particular receptor targets and is currently used as a standard computational tool in drug design for lead compound optimisation and in virtual screening studies to find novel biologically active molecules. The basic tools of a docking methodology include a search algorithm and an energy scoring function for generating and evaluating ligand poses. In this review, we present the search algorithms and scoring functions most commonly used in current molecular docking methods that focus on protein-ligand applications. We summarise the main topics and recent computational and methodological advances in protein-ligand docking. Protein flexibility, multiple ligand binding modes and the free-energy landscape profile for binding affinity prediction are important and interconnected challenges to be overcome by further methodological developments in the docking field.

Discovery of inhibitors to block interactions of HIV-1 integrase with human LEDGF/p75 via structure-based virtual screening and bioassays

This study aims to identify inhibitors that bind at the interface of HIV-1 integrase (IN) and human LEDGF/p75, which represents a novel target for anti-HIV therapy. To date, only a few such inhibitors have been reported. Here structure-based virtual screening was performed to search for the inhibitors from an in-house library of natural products and their derivatives. Among the 38 compounds selected by our strategy, 18 hits were discovered. The two most potent inhibitors showed IC(50) values at 0.32 and 0.26 μM, respectively. Three compounds were subsequently selected for anti-HIV assays, among which (E)-3-(2-chlorophenyl)-1-(2,4-dihydroxyphenyl)prop-2-en-1-one (NPD170) showed the highest antiviral activity (EC(50) = 1.81 μM). The antiviral mechanism of these compounds was further explored, and the results validated that the compounds interrupted the binding of transfected IN to endogenous LEDGF/p75. These findings could be helpful for anti-HIV drug discovery.

Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and prime/MM-GBSA calculation

Developing chemicals that inhibit checkpoint kinase 1 (Chk1) is a promising adjuvant therapeutic to improve the efficacy and selectivity of DNA-targeting agents. Reliable prediction of binding-free energy and binding affinity of Chk1 inhibitors can provide a guide for rational drug design. In this study, multiple docking strategies and Prime/Molecular Mechanics Generalized Born Surface Area (Prime/MM-GBSA) calculation were applied to predict the binding mode and free energy for a series of benzoisoquinolinones as Chk1 inhibitors. Reliable docking results were obtained using induced-fit docking and quantum mechanics/molecular mechanics (QM/MM) docking, which showed superior performance on both ligand binding pose and docking score accuracy to the rigid-receptor docking. Then, the Prime/MM-GBSA method based on the docking complex was used to predict the binding-free energy. The combined use of QM/MM docking and Prime/MM-GBSA method could give a high correlation between the predicted binding-free energy and experimentally determined pIC(50) . The molecular docking combined with Prime/MM-GBSA simulation can not only be used to rapidly and accurately predict the binding-free energy of novel Chk1 inhibitors but also provide a novel strategy for lead discovery and optimization targeting Chk1.

Recent trends and observations in the design of high-quality screening collections

The design of a high-quality screening collection is of utmost importance for the early drug-discovery process and provides, in combination with high-quality assay systems, the foundation of future discoveries. Herein, we review recent trends and observations to successfully expand the access to bioactive chemical space, including the feedback from hit assessment interviews of high-throughput screening campaigns; recent successes with chemogenomics target family approaches, the identification of new relevant target/domain families, diversity-oriented synthesis and new emerging compound classes, and non-classical approaches, such as fragment-based screening and DNA-encoded chemical libraries. The role of in silico library design approaches are emphasized.

The role of computational methods in the identification of bioactive compounds

Computational methods play an ever increasing role in lead finding. A vast repertoire of molecular design and virtual screening methods emerged in the past two decades and are today routinely used. There is increasing awareness that there is no single best computational protocol and correspondingly there is a shift recommending the combination of complementary methods. A promising trend for the application of computational methods in lead finding is to take advantage of the vast amounts of HTS (High Throughput Screening) data to allow lead assessment by detailed systems-based data analysis, especially for phenotypic screens where the identification of compound-target pairs is the primary goal. Herein, we review trends and provide examples of successful applications of computational methods in lead finding.