Ligand-specific scoring functions: improved ranking of docking solutions

Molecular docking is a powerful computational method that has been widely used in many biomolecular studies to predict geometry of a protein-ligand complex. However, while its conformational search algorithms are usually able to generate correct conformation of a ligand in the binding site, the scoring methods often fail to discriminate it among many false variants. We propose to treat this problem by applying more precise ligand-specific scoring filters to re-rank docking solutions. In this way specific features of interactions between protein and different types of compounds can be implicitly taken into account. New scoring functions were constructed including hydrogen bonds, hydrophobic and hydrophilic complementarity terms. These scoring functions also discriminate ligands by the size of the molecule, the total hydrophobicity, and the number of peptide bonds for peptide ligands. Weighting coefficients of the scoring functions were adjusted using a training set of 60 protein-ligand complexes. The proposed method was then tested on the results of docking obtained for an additional 70 complexes. In both cases the success rate was 5-8% better compared to the standard functions implemented in popular docking software.     

Ligand-specific scoring functions: improved ranking of docking solutions1

Molecular docking is a powerful computational method that has been widely used in many biomolecular studies to predict geometry of a protein-ligand complex. However, while its conformational search algorithms are usually able to generate correct conformation of a ligand in the binding site, the scoring methods often fail to discriminate it among many false variants. We propose to treat this problem by applying more precise ligand-specific scoring filters to re-rank docking solutions. In this way specific features of interactions between protein and different types of compounds can be implicitly taken into account. New scoring functions were constructed including hydrogen bonds, hydrophobic and hydrophilic complementarity terms. These scoring functions also discriminate ligands by the size of the molecule, the total hydrophobicity, and the number of peptide bonds for peptide ligands. Weighting coefficients of the scoring functions were adjusted using a training set of 60 protein–ligand complexes. The proposed method was then tested on the results of docking obtained for an additional 70 complexes. In both cases the success rate was 5–8% better compared to the standard functions implemented in popular docking software.     

Supervised consensus scoring for docking and virtual screening

Docking programs are widely used to discover novel ligands efficiently and can predict protein-ligand complex structures with reasonable accuracy and speed. However, there is an emerging demand for better performance from the scoring methods. Consensus scoring (CS) methods improve the performance by compensating for the deficiencies of each scoring function. However, conventional CS and existing scoring functions have the same problems, such as a lack of protein flexibility, inadequate treatment of salvation, and the simplistic nature of the energy function used. Although there are many problems in current scoring functions, we focus our attention on the incorporation of unbound ligand conformations. To address this problem, we propose supervised consensus scoring (SCS), which takes into account protein-ligand binding process using unbound ligand conformations with supervised learning. An evaluation of docking accuracy for 100 diverse protein-ligand complexes shows that SCS outperforms both CS and 11 scoring functions (PLP, F-Score, LigScore, DrugScore, LUDI, X-Score, AutoDock, PMF, G-Score, ChemScore, and D-score). The success rates of SCS range from 89% to 91% in the range of rmsd < 2 A, while those of CS range from 80% to 85%, and those of the scoring functions range from 26% to 76%. Moreover, we also introduce a method for judging whether a compound is active or inactive with the appropriate criterion for virtual screening. SCS performs quite well in docking accuracy and is presumably useful for screening large-scale compound databases before predicting binding affinity.     

GFscore: a general nonlinear consensus scoring function for high-throughput docking

Most of the recent published works in the field of docking and scoring protein/ligand complexes have focused on ranking true positives resulting from a Virtual Library Screening (VLS) through the use of a specified or consensus linear scoring function. In this work, we present a methodology to speed up the High Throughput Screening (HTS) process, by allowing focused screens or for hitlist triaging when a prohibitively large number of hits is identified in the primary screen, where we have extended the principle of consensus scoring in a nonlinear neural network manner. This led us to introduce a nonlinear Generalist scoring Function, GFscore, which was trained to discriminate true positives from false positives in a data set of diverse chemical compounds. This original Generalist scoring Function is a combination of the five scoring functions found in the CScore package from Tripos Inc. GFscore eliminates up to 75% of molecules, with a confidence rate of 90%. The final result is a Hit Enrichment in the list of molecules to investigate during a research campaign for biological active compounds where the remaining 25% of molecules would be sent to in vitro screening experiments. GFscore is therefore a powerful tool for the biologist, saving both time and money.     

Effects of protein conformation in docking: improved pose prediction through protein pocket adaptation

Computational methods for docking ligands have been shown to be remarkably dependent on precise protein conformation, where acceptable results in pose prediction have been generally possible only in the artificial case of re-docking a ligand into a protein binding site whose conformation was determined in the presence of the same ligand (the “cognate” docking problem). In such cases, on well curated protein/ligand complexes, accurate dockings can be returned as top-scoring over 75% of the time using tools such as Surflex-Dock. A critical application of docking in modeling for lead optimization requires accurate pose prediction for novel ligands, ranging from simple synthetic analogs to very different molecular scaffolds. Typical results for widely used programs in the “cross-docking case” (making use of a single fixed protein conformation) have rates closer to 20% success. By making use of protein conformations from multiple complexes, Surflex-Dock yields an average success rate of 61% across eight pharmaceutically relevant targets. Following docking, protein pocket adaptation and rescoring identifies single pose families that are correct an average of 67% of the time. Consideration of the best of two pose families (from alternate scoring regimes) yields a 75% mean success rate.     

Optimizing the signal-to-noise ratio of scoring functions for protein--ligand docking

Empirical scoring functions provide estimates of the free energy of protein-ligand binding in situations when atomic-scale simulations are intractable, for example, in virtual high-throughput screening. Currently, such scoring functions are often inaccurate, and further improvements are complicated by the lack of reliable training data, the complex interplay between scoring functions and docking algorithms, and an inconsistent statistical treatment of positive and negative training data. In comparison to various other performance measures of scoring functions, "analysis of variance" provides a well-behaved objective function for optimization, which focuses on the signal-to-noise ratio of ligand-decoy discrimination. In combination with a large database of ligands and decoys, an in situ optimization of scoring function parameters was able to generate improved, target-specific scoring functions for three different proteins of pharmaceutical interest: cyclin-dependent kinase 2, the estrogen receptor, and cyclooxygenase-2. Statistical analysis of the improvements observed in "receiver-operating characteristic" curves showed that the optimized scoring functions achieved a significantly (between p < 0.0001 and p < 0.05) higher enrichment of true ligands. A scaffold dependence of the resulting binding modes was observed, which is discussed in conjunction with the rigid receptor hypothesis commonly made in protein-ligand docking. In summary, the approach described here represents a well-adapted statistical method for setting up scoring functions.     

Drug efficiency indices for improvement of molecular docking scoring functions

A dataset of protein‐drug complexes with experimental binding energy and crystal structure were analyzed and the performance of different docking engines and scoring functions (as well as components of these) for predicting the free energy of binding and several ligand efficiency indices were compared. The aim was not to evaluate the best docking method, but to determine the effect of different efficiency indices on the experimental and predicted free energy. Some ligand efficiency indices, such as ΔG/W (Wiener index), ΔG/NoC (number of carbons), and ΔG/P (partition coefficient), improve the correlation between experimental and calculated values. This effect was shown to be valid across the different scoring functions and docking programs. It also removes the common bias of scoring functions in favor of larger ligands. For all scoring functions, the efficiency indices effectively normalize the free energy derived indices, to give values closer to experiment. Compound collection filtering can be done prior or after docking, using pharmacokinetic as well as pharmacodynamic profiles. Achieving these better correlations with experiment can improve the ability of docking scoring functions to predict active molecules in virtual screening. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Theory of docking scores and its application to a customizable scoring function

In general, the docking scoring tends to have a size dependence related to the ranking of compounds. In this paper, we describe a novel method of parameter optimization for docking scores which reduce the size dependence and can efficiently discriminate active compounds from chemical databases. This method is based on a simplified theoretical model of docking scores which enables us to utilize large amounts of data of known active and inactive compounds for a particular target without requiring large computational resources or a complicated procedure. This method is useful for making scoring functions for the identification of novel scaffolds using the knowledge of active compounds for a particular target or a customized scoring function for an interesting family of drug targets.     

Antagonist-perturbation mechanism for activation function-2 fixed motifs: active conformation and docking mode of retinoid X receptor antagonists

HX531, which contains a dibenzodiazepine skeleton, is one of the first retinoid X receptor (RXR) antagonists. Functioning via RXR-PPARγ heterodimer, this compound is receiving a lot of attention as a therapeutic drug candidate for diabetic disease controlling differentiation of adipose tissue. However, the active conformation of HX531 for RXRs is not well established. In the present study, quantum mechanics calculations and molecular mechanical docking simulations were carried out to precisely study the docking mode of HX531 with the human RXRα ligand-binding domain, as well as to provide a new approach to drug design using a structure-based perspective. It was suggested that HX531, which has the R configuration for the bent dibenzodiazepine plane together with the equatorial configuration for the N-methyl group attached to the nitrogen atom in the seven-membered diazepine ring, is a typical activation function-2 (AF-2) fixed motif perturbation type antagonist, which destabilizes the formation of AF-2 fixed motifs. On the other hand, the docking simulations supported the experimental result that LG100754 is an RXR homodimer antagonist and an RXR heterodimer agonist.     

Approaches to solving the rigid receptor problem by identifying a minimal set of flexible residues during ligand docking

BACKGROUND: Using fixed receptor sites derived from high-resolution crystal structures in structure-based drug design does not properly account for ligand-induced enzyme conformational change and imparts a bias into the discovery and design of novel ligands. We sought to facilitate the design of improved drug leads by defining residues most likely to change conformation, and then defining a minimal manifold of possible conformations of a target site for drug design based on a small number of identified flexible residues. RESULTS: The crystal structure of thymidylate synthase from an important pathogenic target Pneumocystis carinii (PcTS) bound to its substrate and the inhibitor, BW1843U89, is reported here and reveals a new conformation with respect to the structure of PcTS bound to substrate and the more conventional antifolate inhibitor, CB3717. We developed an algorithm for determining which residues provide 'soft spots' in the protein, regions where conformational adaptation suggests possible modifications for a drug lead that may yield higher affinity. Remodeling the active site of thymidylate synthase with new conformations for only three residues that were identified with this algorithm yields scores for ligands that are compatible with experimental kinetic data. CONCLUSIONS: Based on the examination of many protein/ligand complexes, we develop an algorithm (SOFTSPOTS) for identifying regions of a protein target that are more likely to accommodate plastically to regions of a drug molecule. Using these indicators we develop a second algorithm (PLASTIC) that provides a minimal manifold of possible conformations of a protein target for drug design, reducing the bias in structure-based drug design imparted by structures of enzymes co-crystallized with inhibitors.     

Homology model-based virtual screening for GPCR ligands using docking and target-biased scoring

The current study investigates the combination of two recently reported techniques for the improvement of homology model-based virtual screening for G-protein coupled receptor (GPCR) ligands. First, ligand-supported homology modeling was used to generate receptor models that were in agreement with mutagenesis data and structure-activity relationship information of the ligands. Second, interaction patterns from known ligands to the receptor were applied for scoring and rank ordering compounds from a virtual library using ligand-receptor interaction fingerprint-based similarity (IFS). Our approach was evaluated in retrospective virtual screening experiments for antagonists of the metabotropic glutamate receptor (mGluR) subtype 5. The results of our approach were compared to the results obtained by conventional scoring functions (Dock-Score, PMF-Score, Gold-Score, ChemScore, and FlexX-Score). The IFS lead to significantly higher enrichment rates, relative to the competing scoring functions. Though using a target-biased scoring approach, the results were not biased toward the chemical classes of the reference structures. Our results indicate that the presented approach has the potential to serve as a general setup for successful structure-based GPCR virtual screening.     

Evaluation and application of multiple scoring functions for a virtual screening experiment

In order to identify novel chemical classes of factor Xa inhibitors, five scoring functions (FlexX, DOCK, GOLD, ChemScore and PMF) were engaged to evaluate the multiple docking poses generated by FlexX. The compound collection was composed of confirmed potent factor Xa inhibitors and a subset of the LeadQuest screening compound library. Except for PMF the other four scoring functions succeeded in reproducing the crystal complex (PDB code: 1FAX). During virtual screening the highest hit rate (80%) was demonstrated by FlexX at an energy cutoff of -40 kJ/mol, which is about 40-fold over random screening (2.06%). Limited results suggest that presenting more poses of a single molecule to the scoring functions could deteriorate their enrichment factors. A series of promising scaffolds with favorable binding scores was retrieved from LeadQuest. Consensus scoring by pair-wise intersection failed to enrich the hit rate yielded by single scorings (i.e. FlexX). We note that reported successes of consensus scoring in hit rate enrichment could be artificial because their comparisons were based on a selected subset of single scoring and a markedly reduced subset of double or triple scoring. The findings presented in this report are based upon a single biological system and support further studies.     

Predicting conformations and orientations of guests within a water soluble host: a molecular docking approach

In this study we have examined conformations and orientations of guests within a water-soluble host known by the trivial name Octa Acid (OA). Docking program Vina, which was originally developed for screening drug-like molecules, has been used to identify binding modes and affinities of select guest molecules with OA. Hydrophobic guests were encapsulated into the nonpolar cavity of OA capsule owing to solvophobic interactions. Amphiphilic guests were bound by keeping the nonpolar part within the cavity of OA, while pointing the polar anionic group out of the cavity. All these results obtained from the docking study were in accord with experimental findings. The post-complexation attributes of the guests were regulated by available free space and the specific interactions between guest–OA pair, which led to unusual conformations and orientations. This study showed that scoring function available with Vina, which was derived from protein–ligand data set, could successfully predict post-complexed structural features of guests within OA, thus opening opportunities to modulate physical and chemical behavior of guest molecules.     

Application of semiempirical quantum chemical methods as a scoring function in docking

A crucial point in docking simulations is the scoring function used for estimation of the target-ligand interaction energy. The usual practice is to employ fast but simplified empirical scoring functions. Rigorous quantum chemical methods are too slow to screen virtual combinatorial libraries consisting of thousands of molecules, but they can be used in the final step of the simulations for assessing the results obtained. At this stage quantum chemical calculations can be performed only for the 10–100 top binders predicted by simplified scoring functions, and only using linear-scaling semiempirical quantum chemical methods such as MOZYME. The possibilities and potentialities of the quantum chemical methods for estimation of the binding affinities in docking simulations are a largely unexplored area, so the main goal of this study is a detailed evaluation of the potential and limitations of the MOZYME methodology for estimation of the target-ligand binding energies and its comparison with available experimental data.Proceedings of the 11th International Congress of Quantum Chemistry satellite meeting in honor of Jean-Louis Rivail     

Maximum common binding modes (MCBM): consensus docking scoring using multiple ligand information and interaction fingerprints

Improving the scoring functions for small molecule-protein docking is a highly challenging task in current computational drug design. Here we present a novel consensus scoring concept for the prediction of binding modes for multiple known active ligands. Similar ligands are generally believed to bind to their receptor in a similar fashion. The presumption of our approach was that the true binding modes of similar ligands should be more similar to each other compared to false positive binding modes. The number of conserved (consensus) interactions between similar ligands was used as a docking score. Patterns of interactions were modeled using ligand receptor interaction fingerprints. Our approach was evaluated for four different data sets of known cocrystal structures (CDK-2, dihydrofolate reductase, HIV-1 protease, and thrombin). Docking poses were generated with FlexX and rescored by our approach. For comparison the CScore scoring functions from Sybyl were used, and consensus scores were calculated thereof. Our approach performed better than individual scoring functions and was comparable to consensus scoring. Analysis of the distribution of docking poses by self-organizing maps (SOM) and interaction fingerprints confirmed that clusters of docking poses composed of multiple ligands were preferentially observed near the native binding mode. Being conceptually unrelated to commonly used docking scoring functions our approach provides a powerful method to complement and improve computational docking experiments.     

MiniMuDS: a new optimizer using knowledge-based potentials improves scoring of docking solutions

In small molecule docking, the scoring and ranking of generated conformations is an important, though still not a completely resolved problem. Rescoring schemes often improve the quality of the obtained rankings. It is known that a local optimization is essential before a valid rescore value can be calculated. Here, we present a method that improves rescoring results obtained with the DrugScore function due to a new optimization technique. The method implements a more sophisticated search algorithm compared to the classic local optimization procedures used in this context. We validated the proposed method on a set of 192 protein-ligand complexes. Results show substantial improvements compared to original docking results with success rates increased by up to 10% for top scored solutions below 2 ? root-mean-square deviation to the native state and up to 18% increase below 1 ?, respectively.     

WIN55212-2 docking to the CB1 cannabinoid receptor and multiple pathways for conformational induction

Key pharmacophoric elements for the (aminoalkyl)indole (AAI) CB1 cannabinoid receptor agonists are the aminoalkyl moiety, the lipophilic aroyl group, and the heterocyclic indole ring. In the present study, the docking space allowed for (R)-[2,3-dihydro-5-methyl-3-[(4-morpholinyl)methyl]pyrrolo[1,2,3-de]-1,4-benzoxazin-6-yl](1-naphthalenyl)methanone (WIN55212-2; 1) within the CB1 receptor was extensively explored by a docking approach that combines Monte Carlo (MC) and molecular dynamics (MD) simulations. The goals were to understand the key binding interactions of AAIs within the CB1 receptor and to examine the role of the ligand in inducing a receptor conformational change. From the findings of extensive SAR studies on the cannabinoid compounds and correlation between AAI binding affinity data and calculated binding energies, we proposed two alternative binding conformations, aroyl-up1 and aroyl-up2. These denote the directionality of the ligand naphthyl ring within the receptor upward with respect to the extracellular side. A comprehensive structural analysis of 1 demonstrated that the aroyl ring moiety could be important as the steric trigger for inducing CB1 receptor conformational change. Thus, it appears that aromatic-aromatic interactions are important not only for the binding of 1 but also for inducing receptor conformational change. It is possible that differences in the nature of the ligand binding could contribute to ligand-specific conformational changes in the receptor.     

Conformational and docking studies of acyl homoserine lactones as a robust method to investigate bioactive conformations

A method aiming at investigating possible bioactive conformations of acyl homoserine lactone (AHL) quorum sensing (QS) modulators is established. The method relies on the exhaustive conformational analysis of AHLs by varying torsion angles around the amide group then on the selection of the closest conformation to those known from co-crystallized XRD data of AHL-receptor complexes. These latter are then docked as rigid ligand within the receptor binding site, leading to interactions with binding site residues which are highly consistent as compared with the data arising from XRD studies. The method is first validated using AHLs for which XRD data of their complexes with their cognate receptor are available, then extended to examples for which the binding mode is still unknown.Three compounds were used to validate the method: hexanoyl homoserine lactone (HHL) as an example of autoinducer, 3-oxo-butanoyl homoserine lactone (OBHL), as a representative model of 3-oxo-AHLs, and 4-(4-chlorophenoxy)butanoyl homoserine lactone (CPOBHL) as an example of a QS inhibitor. The conformational analysis of these three compounds to their cognate protein (TraR, SdiA, LasR and CviR) provides the data which enable the next rigid docking step. Further rigid docking of the closest conformations compared to the known bioactive ones within the binding sites allows to recover the expected binding mode with high precision (atomic RMSD < 2 Å). This “conformational analysis/torsion angle filter/rigid ligand docking” method was then used for investigating three non-natural AHL-type QS inhibitors without known co-crystallized XRD structures, namely was 2-hexenoyl homoserine lactone (HenHL), 3-oxo-4-phenylbutanoyl homoserine lactone (OPBHL) and 3-(4-bromophenyl)propanoyl homoserine lactone (BPPHL). Results provide insights into their possible binding mode by identifying specific interactions with some key residues within the receptor binding site, allowing discussion of their biological activity.     

Distilling the essential features of a protein surface for improving protein-ligand docking,scoring, and virtual screening

For the successful identification and docking of new ligands to a protein target by virtual screening, the essential features of the protein and ligand surfaces must be captured and distilled in an efficient representation. Since the running time for docking increases exponentially with the number of points representing the protein and each ligand candidate, it is important to place these points where the best interactions can be made between the protein and the ligand. This definition of favorable points of interaction can also guide protein structure-based ligand design, which typically focuses on which chemical groups provide the most energetically favorable contacts. In this paper, we present an alternative method of protein template and ligand interaction point design that identifies the most favorable points for making hydrophobic and hydrogen–bond interactions by using a knowledge base. The knowledge-based protein and ligand representations have been incorporated in version 2.0 of SLIDE and resulted in dockings closer to the crystal structure orientations when screening a set of 57 known thrombin and glutathione S–transferase (GST) ligands against the apo structures of these proteins. There was also improved scoring enrichment of the dockings, meaning better differentiation between the chemically diverse known ligands and a 15,000-molecule dataset of randomly-chosen small organic molecules. This approach for identifying the most important points of interaction between proteins and their ligands can equally well be used in other docking and design techniques. While much recent effort has focused on improving scoring functions for protein-ligand docking, our results indicate that improving the representation of the chemistry of proteins and their ligands is another avenue that can lead to significant improvements in the identification, docking, and scoring of ligands.(These authors contributed equally to this work)     

Using correlated parameters for improved ranking of protein-protein docking decoys

A successful protein–protein docking study culminates in identification of decoys at top ranks with near‐native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near‐native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein–protein interface. In this work, we describe a novel method combining an interface‐size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein–protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having ≤10 Å backbone root mean square deviation from true binding geometry. Comparisons with other state‐of‐art methods confirm the improved ranking power of our method without the use of any experiment‐guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less‐difficult bound binary protein–protein complexes with ≤35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with ≤5Å backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011.