Generation of Quantum Configurational Ensembles Using Approximate Potentials.

Conformational analysis is of paramount importance in drug design: it is crucial to determine pharmacological properties, understand molecular recognition processes, and characterize the conformations of ligands when unbound. Molecular Mechanics (MM) simulation methods, such as Monte Carlo (MC) and molecular dynamics (MD), are usually employed to generate ensembles of structures due to their ability to extensively sample the conformational space of molecules. The accuracy of these MM-based schemes strongly depends on the functional form of the force field (FF) and its parametrization, components that often hinder their performance. High-level methods, such as ab initio MD, provide reliable structural information but are still too computationally expensive to allow for extensive sampling. Therefore, to overcome these limitations, we present a multilevel MC method that is capable of generating quantum configurational ensembles while keeping the computational cost at a minimum. We show that FF reparametrization is an efficient route to generate FFs that reproduce QM results more closely, which, in turn, can be used as low-cost models to achieve the gold standard QM accuracy. We demonstrate that the MC acceptance rate is strongly correlated with various phase space overlap measurements and that it constitutes a robust metric to evaluate the similarity between the MM and QM levels of theory. As a more advanced application, we present a self-parametrizing version of the algorithm, which combines sampling and FF parametrization in one scheme, and apply the methodology to generate the QM/MM distribution of a ligand in aqueous solution.

Fragment-Guided Discovery of Pyrazole Carboxylic Acid Inhibitors of the Kelch-like ECH-Associated Protein 1: Nuclear Factor Erythroid 2 Related Factor 2 (KEAP1:NRF2) Protein-Protein Interaction.

The NRF2-mediated cytoprotective response is central to cellular homoeostasis, and there is increasing interest in developing small-molecule activators of this pathway as therapeutics for diseases involving chronic oxidative stress. The protein KEAP1, which regulates NRF2, is a key point for pharmacological intervention, and we recently described the use of fragment-based drug discovery to develop a tool compound that directly disrupts the protein-protein interaction between NRF2 and KEAP1. We now present the identification of a second, chemically distinct series of KEAP1 inhibitors, which provided an alternative chemotype for lead optimization. Pharmacophoric information from our original fragment screen was used to identify new hit matter through database searching and to evolve this into a new lead with high target affinity and cell-based activity. We highlight how knowledge obtained from fragment-based approaches can be used to focus additional screening campaigns in order to de-risk projects through the rapid identification of novel chemical series.

ParaMol: A Package for Automatic Parameterization of Molecular Mechanics Force Fields.

The ensemble of structures generated by molecular mechanics (MM) simulations is determined by the functional form of the force field employed and its parameterization. For a given functional form, the quality of the parameterization is crucial and will determine how accurately we can compute observable properties from simulations. While accurate force field parameterizations are available for biomolecules, such as proteins or DNA, the parameterization of new molecules, such as drug candidates, is particularly challenging as these may involve functional groups and interactions for which accurate parameters may not be available. Here, in an effort to address this problem, we present ParaMol, a Python package that has a special focus on the parameterization of bonded and nonbonded terms of druglike molecules by fitting to ab initio data. We demonstrate the software by deriving bonded terms' parameters of three widely known drug molecules, viz. aspirin, caffeine, and a norfloxacin analogue, for which we show that, within the constraints of the functional form, the methodologies implemented in ParaMol are able to derive near-ideal parameters. Additionally, we illustrate the best practices to follow when employing specific parameterization routes. We also determine the sensitivity of different fitting data sets, such as relaxed dihedral scans and configurational ensembles, to the parameterization procedure, and discuss the features of the various weighting methods available to weight configurations. Owing to ParaMol's capabilities, we propose that this software can be introduced as a routine step in the protocol normally employed to parameterize druglike molecules for MM simulations.

Practical High-Quality Electrostatic Potential Surfaces for Drug Discovery Using a Graph-Convolutional Deep Neural Network.

Inspecting protein and ligand electrostatic potential (ESP) surfaces in order to optimize electrostatic complementarity is a key activity in drug design. These ESP surfaces need to reflect the true electrostatic nature of the molecules, which typically means time-consuming high-level quantum mechanics (QM) calculations are required. For interactive design much faster alternative methods are required. Here, we present a graph convolutional deep neural network (DNN) model, trained on ESP surfaces derived from high quality QM calculations, that generates ESP surfaces for ligands in a fraction of a second. Additionally, we describe a method for constructing fast QM-trained ESP surfaces for proteins. We show that the DNN model generates ESP surfaces that are in good agreement with QM and that the ESP values correlate well with experimental properties relevant to medicinal chemistry. We believe that these high-quality, interactive ESP surfaces form a powerful tool for driving drug discovery programs forward. The trained model and associated code are available from https://github.com/AstexUK/ESP_DNN.

The Fragment Network: A Chemistry Recommendation Engine Built Using a Graph Database.

The hit validation stage of a fragment-based drug discovery campaign involves probing the SAR around one or more fragment hits. This often requires a search for similar compounds in a corporate collection or from commercial suppliers. The Fragment Network is a graph database that allows a user to efficiently search chemical space around a compound of interest. The result set is chemically intuitive, naturally grouped by substitution pattern and meaningfully sorted according to the number of observations of each transformation in medicinal chemistry databases. This paper describes the algorithms used to construct and search the Fragment Network and provides examples of how it may be used in a drug discovery context.

Predicting "Hot" and "Warm" Spots for Fragment Binding.

Computational fragment mapping methods aim to predict hotspots on protein surfaces where small fragments will bind. Such methods are popular for druggability assessment as well as structure-based design. However, to date researchers developing or using such tools have had no clear way of assessing the performance of these methods. Here, we introduce the first diverse, high quality validation set for computational fragment mapping. The set contains 52 diverse examples of fragment binding "hot" and "warm" spots from the Protein Data Bank (PDB). Additionally, we describe PLImap, a novel protocol for fragment mapping based on the Protein-Ligand Informatics force field (PLIff). We evaluate PLImap against the new fragment mapping test set, and compare its performance to that of simple shape-based algorithms and fragment docking using GOLD. PLImap is made publicly available from https://bitbucket.org/AstexUK/pli .

Protein-Ligand Informatics Force Field (PLIff): Toward a Fully Knowledge Driven "Force Field" for Biomolecular Interactions.

The Protein Data Bank (PDB) contains a wealth of data on nonbonded biomolecular interactions. If this information could be distilled down to nonbonded interaction potentials, these would have some key advantages over standard force fields. However, there are some important outstanding issues to address in order to do this successfully. This paper introduces the protein-ligand informatics "force field", PLIff, which begins to address these key challenges ( https://bitbucket.org/AstexUK/pli ). As a result of their knowledge-based nature, the next-generation nonbonded potentials that make up PLIff automatically capture a wide range of interaction types, including special interactions that are often poorly described by standard force fields. We illustrate how PLIff may be used in structure-based design applications, including interaction fields, fragment mapping, and protein-ligand docking. PLIff performs at least as well as state-of-the art scoring functions in terms of pose predictions and ranking compounds in a virtual screening context.

Detection of secondary binding sites in proteins using fragment screening.

Proteins need to be tightly regulated as they control biological processes in most normal cellular functions. The precise mechanisms of regulation are rarely completely understood but can involve binding of endogenous ligands and/or partner proteins at specific locations on a protein that can modulate function. Often, these additional secondary binding sites appear separate to the primary binding site, which, for example for an enzyme, may bind a substrate. In previous work, we have uncovered several examples in which secondary binding sites were discovered on proteins using fragment screening approaches. In each case, we were able to establish that the newly identified secondary binding site was biologically relevant as it was able to modulate function by the binding of a small molecule. In this study, we investigate how often secondary binding sites are located on proteins by analyzing 24 protein targets for which we have performed a fragment screen using X-ray crystallography. Our analysis shows that, surprisingly, the majority of proteins contain secondary binding sites based on their ability to bind fragments. Furthermore, sequence analysis of these previously unknown sites indicate high conservation, which suggests that they may have a biological function, perhaps via an allosteric mechanism. Comparing the physicochemical properties of the secondary sites with known primary ligand binding sites also shows broad similarities indicating that many of the secondary sites may be druggable in nature with small molecules that could provide new opportunities to modulate potential therapeutic targets.

Potential and limitations of ensemble docking.

A major problem in structure-based virtual screening applications is the appropriate selection of a single or even multiple protein structures to be used in the virtual screening process. A priori it is unknown which protein structure(s) will perform best in a virtual screening experiment. We investigated the performance of ensemble docking, as a function of ensemble size, for eight targets of pharmaceutical interest. Starting from single protein structure docking results, for each ensemble size up to 500,000 combinations of protein structures were generated, and, for each ensemble, pose prediction and virtual screening results were derived. Comparison of single to multiple protein structure results suggests improvements when looking at the performance of the worst and the average over all single protein structures to the performance of the worst and average over all protein ensembles of size two or greater, respectively. We identified several key factors affecting ensemble docking performance, including the sampling accuracy of the docking algorithm, the choice of the scoring function, and the similarity of database ligands to the cocrystallized ligands of ligand-bound protein structures in an ensemble. Due to these factors, the prospective selection of optimum ensembles is a challenging task, shown by a reassessment of published ensemble selection protocols.

Experiences in fragment-based drug discovery.

Fragment-based drug discovery (FBDD) has become established in both industry and academia as an alternative approach to high-throughput screening for the generation of chemical leads for drug targets. In FBDD, specialised detection methods are used to identify small chemical compounds (fragments) that bind to the drug target, and structural biology is usually employed to establish their binding mode and to facilitate their optimisation. In this article, we present three recent and successful case histories in FBDD. We then re-examine the key concepts and challenges of FBDD with particular emphasis on recent literature and our own experience from a substantial number of FBDD applications. Our opinion is that careful application of FBDD is living up to its promise of delivering high quality leads with good physical properties and that in future many drug molecules will be derived from fragment-based approaches.

Docking performance of fragments and druglike compounds.

This paper addresses two questions of key interest to researchers working with protein-ligand docking methods: (i) Why is there such a large variation in docking performance between different test sets reported in the literature? (ii) Are fragments more difficult to dock than druglike compounds? To answer these, we construct a test set of in-house X-ray structures of protein-ligand complexes from drug discovery projects, half of which contain fragment ligands, the other half druglike ligands. We find that a key factor affecting docking performance is ligand efficiency (LE). High LE compounds are significantly easier to dock than low LE compounds, which we believe could explain the differences observed between test sets reported in the literature. There is no significant difference in docking performance between fragments and druglike compounds, but the reasons why dockings fail appear to be different.

Protein-ligand docking against non-native protein conformers.

In the validation of protein-ligand docking protocols, performance is mostly measured against native protein conformers, i.e. each ligand is docked into the protein conformation from the structure that contained that ligand. In real-life applications, however, ligands are docked against non-native conformations of the protein, i.e. the apo structure or a structure of a different protein-ligand complex. Here, we have constructed an extensive test set for assessing docking performance against non-native protein conformations. This new test set is built on the Astex Diverse Set (which we recently constructed for assessing native docking performance) and contains 1112 non-native structures for 65 drug targets. Using the protein-ligand docking program GOLD, the Astex Diverse Set and the new Astex Non-native Set, we established that, whereas docking performance (top-ranked solution within 2 A rmsd of the experimental binding mode) is approximately 80% for native docking, this drops to 61% for non-native docking. A similar drop-off is observed for sampling performance (any solution within 2 A): 91% for native docking vs 72% for non-native docking. No significant differences were observed between docking performance against apo and nonapo structures. We found that, whereas small variations in protein conformation are generally tolerated by our rigid docking protocol, larger protein movements result in a catastrophic drop-off in performance. Some docking performance and nearly all sampling performance can be recovered by considering dockings produced against a small number of non-native structures simultaneously. Docking against non-native structures of complexes containing ligands that are similar to the docked ligand also significantly improves both docking performance and sampling performance.

Identification of inhibitors of protein kinase B using fragment-based lead discovery.

Using fragment-based screening techniques, 5-methyl-4-phenyl-1H-pyrazole (IC50 80 microM) was identified as a novel, low molecular weight inhibitor of protein kinase B (PKB). Herein we describe the rapid elaboration of highly potent and ligand efficient analogues using a fragment growing approach. Iterative structure-based design was supported by protein-ligand structure determinations using a PKA-PKB "chimera" and a final protein-ligand structure of a lead compound in PKBbeta itself.

A structural comparison of inhibitor binding to PKB, PKA and PKA-PKB chimera.

Although the crystal structure of the anti-cancer target protein kinase B (PKBbeta/Akt-2) has been useful in guiding inhibitor design, the closely related kinase PKA has generally been used as a structural mimic due to its facile crystallization with a range of ligands. The use of PKB-inhibitor crystallography would bring important benefits, including a more rigorous understanding of factors dictating PKA/PKB selectivity, and the opportunity to validate the utility of PKA-based surrogates. We present a "back-soaking" method for obtaining PKBbeta-ligand crystal structures, and provide a structural comparison of inhibitor binding to PKB, PKA, and PKA-PKB chimera. One inhibitor presented here exhibits no PKB/PKA selectivity, and the compound adopts a similar binding mode in all three systems. By contrast, the PKB-selective inhibitor A-443654 adopts a conformation in PKB and PKA-PKB that differs from that with PKA. We provide a structural explanation for this difference, and highlight the ability of PKA-PKB to mimic the true PKB binding mode in this case.

Diverse, high-quality test set for the validation of protein-ligand docking performance.

A procedure for analyzing and classifying publicly available crystal structures has been developed. It has been used to identify high-resolution protein-ligand complexes that can be assessed by reconstructing the electron density for the ligand using the deposited structure factors. The complexes have been clustered according to the protein sequences, and clusters have been discarded if they do not represent proteins thought to be of direct interest to the pharmaceutical or agrochemical industry. Rules have been used to exclude complexes containing non-drug-like ligands. One complex from each cluster has been selected where a structure of sufficient quality was available. The final Astex diverse set contains 85 diverse, relevant protein-ligand complexes, which have been prepared in a format suitable for docking and are to be made freely available to the entire research community (http://www.ccdc.cam.ac.uk). The performance of the docking program GOLD against the new set is assessed using a variety of protocols. Relatively unbiased protocols give success rates of approximately 80% for redocking into native structures, but it is possible to get success rates of over 90% with some protocols.

Protein-Ligand Complexes: Computation of the Relative Free Energy of Different Scaffolds and Binding Modes.

A methodology for the calculation of the free energy difference between a pair of molecules of arbitrary topology is proposed. The protocol relies on a dual-topology paradigm, a softening of the intermolecular interactions, and a constraint that prevents the perturbed molecules from drifting away from each other at the end states. The equivalence and the performance of the methodology against a single-topology approach are demonstrated on a pair of harmonic oscillators, the calculation of the relative solvation free energy of ethane and methanol, and the relative binding free energy of two congeneric inhibitors of cyclooxygenase 2. The stability of two alternative binding modes of an inhibitor of cyclin-dependent kinase 2 is then investigated. Finally, the relative binding free energy of two structurally different inhibitors of cyclin-dependent kinase 2 is calculated. The proposed methodology allows the study of a range of problems that are beyond the reach of traditional relative free energy calculation protocols and should prove useful in drug design studies.

Protein-ligand binding affinity predictions by implicit solvent simulations: a tool for lead optimization?

Continuum electrostatics is combined with rigorous free-energy calculations in an effort to deliver a reliable and efficient method for in silico lead optimization. The methodology is tested by calculation of the relative binding free energies of a set of inhibitors of neuraminidase, cyclooxygenase2, and cyclin-dependent kinase 2. The calculated free energies are compared to the results obtained with explicit solvent simulations and empirical scoring functions. For cyclooxygenase2, deficiencies in the continuum electrostatics theory are identified and corrected with a modified simulation protocol. For neuraminidase, it is shown that a continuum representation of the solvent leads to markedly different protein-ligand interactions compared to the explicit solvent simulations, and a reconciliation of the two protocols is problematic. Cyclin-dependent kinase 2 proves more challenging, and none of the methods employed in this study yield high quality predictions. Despite the differences observed, for these systems, the use of an implicit solvent framework to predict the ranking of congeneric inhibitors to a protein is shown to be faster, as accurate or more accurate than the explicit solvent protocol, and superior to empirical scoring schemes.

Automated protein-ligand crystallography for structure-based drug design.

An approach to automate protein-ligand crystallography is presented, with the aim of increasing the number of structures available to structure-based drug design. The methods we propose deal with the automatic interpretation of diffraction data for targets with known protein structures, and provide easy access to the results. Central to the system is a novel procedure that fully automates the placement of ligands into electron density maps. Automation provides an objective way to structure solution, whereas manual placement can be rather subjective, especially for data of low to medium resolution. Ligands are placed by docking into electron density, whilst taking care of protein-ligand interactions. The ligand fitting procedure has been validated on both public domain and in-house examples. Some of the latter deal with cocktails of low-molecular weight compounds, as used in fragment-based drug discovery by crystallography. For such library-screening experiments we show that the method can automatically identify which of the compounds from a cocktail is bound.