Prospective Evaluation of Free Energy Calculations for the Prioritization of Cathepsin L Inhibitors.

Improving the binding affinity of a chemical series by systematically probing one of its exit vectors is a medicinal chemistry activity that can benefit from molecular modeling input. Herein, we compare the effectiveness of four approaches in prioritizing building blocks with better potency: selection by a medicinal chemist, manual modeling, docking followed by manual filtering, and free energy calculations (FEP). Our study focused on identifying novel substituents for the apolar S2 pocket of cathepsin L and was conducted entirely in a prospective manner with synthesis and activity determination of 36 novel compounds. We found that FEP selected compounds with improved affinity for 8 out of 10 picks compared to 1 out of 10 for the other approaches. From this result and other additional analyses, we conclude that FEP can be a useful approach to guide this type of medicinal chemistry optimization once it has been validated for the system under consideration.

[1]  John M. Barnard,et al.  Chemical Similarity Searching , 1998, J. Chem. Inf. Comput. Sci..

[2]  Matthias Rarey,et al.  Torsion Library Reloaded: A New Version of Expert-Derived SMARTS Rules for Assessing Conformations of Small Molecules , 2016, J. Chem. Inf. Model..

[3]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[4]  Daniel Cappel,et al.  Accurate Binding Free Energy Predictions in Fragment Optimization , 2015, J. Chem. Inf. Model..

[5]  F. Diederich,et al.  Potent and Selective Inhibition of Cysteine Proteases from Plasmodium falciparum and Trypanosoma brucei , 2010, ChemMedChem.

[6]  Jennifer L. Knight,et al.  Modeling Local Structural Rearrangements Using FEP/REST: Application to Relative Binding Affinity Predictions of CDK2 Inhibitors. , 2013, Journal of chemical theory and computation.

[7]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[8]  Y. Cheng,et al.  Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction. , 1973, Biochemical pharmacology.

[9]  Holger Gohlke,et al.  Binding Free Energy Calculations for Lead Optimization: Assessment of Their Accuracy in an Industrial Drug Design Context. , 2014, Journal of chemical theory and computation.

[10]  Jennifer L. Knight,et al.  Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field. , 2015, Journal of the American Chemical Society.

[11]  Shuai Liu,et al.  D3R grand challenge 2015: Evaluation of protein–ligand pose and affinity predictions , 2016, Journal of Computer-Aided Molecular Design.

[12]  P. Leeson,et al.  The influence of drug-like concepts on decision-making in medicinal chemistry , 2007, Nature Reviews Drug Discovery.

[13]  Atli Thorarensen,et al.  Imidazotriazines: Spleen Tyrosine Kinase (Syk) Inhibitors Identified by Free‐Energy Perturbation (FEP) , 2016, ChemMedChem.

[14]  Robert Abel,et al.  Leveraging Data Fusion Strategies in Multireceptor Lead Optimization MM/GBSA End-Point Methods. , 2014, Journal of chemical theory and computation.

[15]  David Rogers,et al.  Extended-Connectivity Fingerprints , 2010, J. Chem. Inf. Model..

[16]  Jennifer L. Knight,et al.  OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. , 2016, Journal of chemical theory and computation.

[17]  Emilio Gallicchio,et al.  Advances in all atom sampling methods for modeling protein-ligand binding affinities. , 2011, Current opinion in structural biology.

[18]  Paul R. Gerber,et al.  MAB, a generally applicable molecular force field for structure modelling in medicinal chemistry , 1995, J. Comput. Aided Mol. Des..

[19]  W Patrick Walters,et al.  A detailed comparison of current docking and scoring methods on systems of pharmaceutical relevance , 2004, Proteins.

[20]  Martin Stahl,et al.  Small Molecule Conformational Preferences Derived from Crystal Structure Data. A Medicinal Chemistry Focused Analysis , 2008, J. Chem. Inf. Model..

[21]  S. Buchwald,et al.  Mild Palladium-Catalyzed Cyanation of (Hetero)aryl Halides and Triflates in Aqueous Media , 2015, Organic letters.

[22]  F. Diederich,et al.  Optimization of Triazine Nitriles as Rhodesain Inhibitors: Structure–Activity Relationships, Bioisosteric Imidazopyridine Nitriles, and X‐ray Crystal Structure Analysis with Human Cathepsin L , 2013, ChemMedChem.

[23]  Matthias Rarey,et al.  Torsion angle preferences in druglike chemical space: a comprehensive guide. , 2013, Journal of medicinal chemistry.

[24]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[25]  Robert Abel,et al.  Accelerating drug discovery through tight integration of expert molecular design and predictive scoring. , 2017, Current opinion in structural biology.

[26]  R. Glen,et al.  Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation. , 1995, Journal of molecular biology.

[27]  A. Hopkins,et al.  Ligand efficiency: a useful metric for lead selection. , 2004, Drug discovery today.

[28]  F. Diederich,et al.  Tuning and predicting biological affinity: aryl nitriles as cysteine protease inhibitors. , 2012, Organic & biomolecular chemistry.

[29]  H. Berendsen,et al.  THERMODYNAMICS OF CAVITY FORMATION IN WATER - A MOLECULAR-DYNAMICS STUDY , 1982 .

[30]  W. L. Jorgensen,et al.  Monte Carlo simulation of differences in free energies of hydration , 1985 .

[31]  Dariusz Plewczynski,et al.  Can we trust docking results? Evaluation of seven commonly used programs on PDBbind database , 2011, J. Comput. Chem..

[32]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[33]  Laura Pérez-Benito,et al.  Application of Free Energy Perturbation for the Design of BACE1 Inhibitors , 2016, J. Chem. Inf. Model..

[34]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[35]  Martin Stahl,et al.  Rationalizing Tight Ligand Binding through Cooperative Interaction Networks , 2011, J. Chem. Inf. Model..