Performance of Relative Binding Free Energy Calculations on an Automatically Generated Dataset of Halogen-Deshalogen Matched Molecular Pairs.

In this study, we generated a matched molecular pair dataset of halogen/deshalogen compounds with reliable binding affinity data and structural binding mode information from public databases. The workflow includes automated system preparation and setup of free energy perturbation relative binding free energy calculations. We demonstrate the suitability of these datasets to investigate the performance of molecular mechanics force fields and molecular simulation algorithms for the purpose of in silico affinity predictions in lead optimization. Our datasets of a total of 115 matched molecular pairs show highly accurate binding free energy predictions with an average error of <1 kcal/mol despite the semi-automated calculation scheme. We quantify the accuracy of the optimized potential for liquid simulations (OPLS) force field to predict the effect of halogen addition to compounds, a commonly employed chemical modification in the design of drug-like molecules.

MM/GBSA binding energy prediction on the PDBbind data set: successes, failures, and directions for further improvement.

We validate an automated implementation of a combined Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) method (VSGB 2.0 energy model) on a large and diverse selection of protein-ligand complexes (855 complexes). Although this data set is diverse with respect to both protein families and ligands, after carefully removing flawed structures, a significant correlation (R(2) = 0.63) between calculated and experimental binding affinities is obtained. Consistent explanations for "outlier" complexes are found. Visual analysis of the crystal structures and recourse to the original literature reveal that neglect of explicit solvent, ligand strain, and entropy contribute to the under- and overestimation of computed affinities. The limits of the Molecular Mechanics/Implicit Solvent approach to accurately estimate protein-ligand binding affinities is discussed as is the influence of the quality of protein-ligand complexes on computed free energy binding values.

Optimization and validation of a docking-scoring protocol; application to virtual screening for COX-2 inhibitors.

To exploit available structural information about the cyclooxygenase enzyme for the virtual screening of large chemical libraries, a docking-scoring protocol was tuned and validated. The screening accuracy was assessed using a series of known inhibitors and a set of diverse a priori inactive compounds that was seeded with known active ligands. The major parameters of the DOCK algorithm were investigated. A large improvement of the results was obtained on tweaking some of them. The generated complexes were rescored using six scoring functions. In this way, the striking importance of this step was demonstrated, as well as the good performances of DOCK energy and SCORE for this target. The results were further improved via a consensus approach. As a first application, a subset of a large compound library was screened using this protocol. Among the compounds that were selected for biological testing, a third of them turned out to have a significant enzyme inhibition.

2D QSAR consensus prediction for high-throughput virtual screening. An application to COX-2 inhibition modeling and screening of the NCI database.

Using classification (SOM, LVQ, Binary, Decision Tree) and regression algorithms (PLS, BRANN, k-NN, Linear), this paper details the building of eight 2D-QSAR models from a 266 COX-2 inhibitor training set. The predictive performances of these eight models were subsequently compared using an 88 COX-2 inhibitor test set. Each ligand is described by 52 2D descriptors expressed as van der Waals Surface Areas (P_VSA) and its COX-2 binding IC50. One of our best predictive models is the neural network model (BRANN), which is able to select a subset, from the 88 ligand test set, that contains 94% COX-2 active inhibitors (pIC50>7.5) and detects 71% of all the actives. We then introduce a QSAR consensus prediction protocol that is shown to be more predictive than any single QSAR model: our C3 consensus approach is able to select a subset from the 88 ligand test set that contains 94% active inhibitors and 83% of all the actives. The 2D QSAR consensus protocol was finally applied to the high-throughput virtual screening of the NCI database, containing 193,477 organic compounds.