New approach to molecular docking and its application to virtual screening of chemical databases

Molecular Recognition of Moringa oleifera Active Compounds for Stunted Growth Prevention Using Network Pharmacology and Molecular Modeling Approach

In this study, network pharmacology was used to analyze the active compounds of Moringa oleifera as food supplements for stunted growth prevention. Thirty-eight important proteins were discovered that may be strongly related to stunting. Those proteins were uploaded to several online tool platforms in order to determine the shared genes' pathways. Six pathways were identified that may be correlated with human growth. Furthermore, ligands for molecular docking analysis were retrieved from the top 5 active substances discovered through experimental investigation. In the meantime, the first-degree rank based on the protein-protein interaction (PPI) topological analysis was utilized to choose albumin protein (ALB) as a receptor. Our docking results showed that every ligand binds to the receptors, indicating that they can bind to the binding site of the ALB protein to form a complex formation. Further, MD simulation was used to verify the stability of the ligand in complex with the protein in the TIP3P water model. Based on the validation parameters, our results suggested that all models achieved a stable phase along the simulation. Additionally, the MM-GBSA method was used to calculate the binding energies of all models. Ligands 2 and 4 have strong binding to the binding pocket of ALB, followed by ligands 3, 5, and 2, suggesting that those ligands could be promising food supplements that can be utilized for the prevention of stunted growth in children.

In Silico Investigations into the Selectivity of Psychoactive and New Psychoactive Substances in Monoamine Transporters

New psychoactive substances (NPS) are a group of compounds that mimic the effects of illicit substances. A range of NPS have been shown to interact with the three main classes of monoamine transporters (DAT, NET, and SERT) to differing extents, but it is unclear why these differences arise. To aid in understanding the differences in affinity between the classes of monoamine transporters, several in silico experiments were conducted. Docking experiments showed there was no direct correlation between a range of scoring functions and experimental activity, but Spearman ranking analysis showed a significant correlation (α = 0.1) for DAT, with the affinity ΔG (0.42), αHB (0.40), GoldScore (0.40), and PLP (0.41) scoring functions, and for DAT (0.38) and SERT (0.40) using a consensus scoring approach. Qualitative structure-activity relationship (QSAR) experiments resulted in the generation of robust and predictive three-descriptor models for SERT (r ² = 0.87, q ² = 0.8, and test set r ² = 0.74) and DAT (r ² = 0.68, q ² = 0.51, test set r ² = 0.63). Both QSAR models described similar characteristics for binding, i.e., rigid hydrophobic molecules with a biogenic amine moiety, and were not sufficient to facilitate a deeper understanding of differences in affinity between the monoamine transporters. This contextualizes the observed promiscuity for NPS between the isoforms and highlights the difficulty in the design and development of compounds that are isoform-selective.

Disease-Ligand Identification Based on Flexible Neural Tree

In order to screen the disease-related compounds of a traditional Chinese medicine prescription in network pharmacology research accurately, a new virtual screening method based on flexible neural tree (FNT) model, hybrid evolutionary method and negative sample selection algorithm is proposed. A novel hybrid evolutionary algorithm based on the Grammar-guided genetic programming and salp swarm algorithm is proposed to infer the optimal FNT. According to hypertension, diabetes, and Corona Virus Disease 2019, disease-related compounds are collected from the up-to-date literatures. The unrelated compounds are chosen by negative sample selection algorithm. ECFP6, MACCS, Macrocycle, and RDKit are utilized to numerically characterize the chemical structure of each compound collected, respectively. The experiment results show that our proposed method performs better than classical classifiers [Support Vector Machine (SVM), random forest (RF), AdaBoost, decision tree (DT), Gradient Boosting Decision Tree (GBDT), KNN, logic regression (LR), and Naive Bayes (NB)], up-to-date classifier (gcForest), and deep learning method (forgeNet) in terms of AUC, ROC, TPR, FPR, Precision, Specificity, and F1. MACCS method is suitable for the maximum number of classifiers. All methods perform poorly with ECFP6 molecular descriptor.

Cyanopropyl functionalized benzimidazolium salts and their silver N-heterocyclic carbene complexes: Synthesis, antimicrobial activity, and theoretical analysis

The reaction of N-substituted benzimidazole with 4-bromobutyronitrile gives the corresponding benzimidazolium salts as N-heterocyclic carbene (NHC) precursors. Silver(I) carbene complexes are synthesized by the reaction of the corresponding benzimidazolium salts with Ag₂ O in dichloromethane. These new NHC precursors and Ag-NHC complexes were characterized by spectroscopy techniques and also screened for their antibacterial activities against the standard bacterial strains Escherichia coli, Pseudomonas aeruginosa, Acinetobacter baumannii, Klebsiella pneumoniae, Staphylococcus aureus, methicillin-resistant Staphylococcus aureus, and Enterococcus faecalis, and the standard fungal strains Candida albicans and Candida glabrata, and promising results were achieved. The compounds were also analyzed by density functional theory (DFT)/time-dependent DFT and docking methods.

Molecular modelling on SARS-CoV-2 papain-like protease: an integrated study with homology modelling, molecular docking, and molecular dynamics simulations

SARS-CoV-2 PLpro was investigated as a therapeutic target for potent antiviral drugs due to its essential role in not only viral replication but also in regulating the inborn immune response. Several computational approaches, including homology modelling, molecular docking, and molecular dynamics (MD) studies, were employed to search for promising drugs in treating SARS-CoV-2. Eighty-one compounds, sub-structurally similar to the antiviral drug, were used as potential inhibitors of PLpro. From our results, three complexes containing the ligands with Pubchem IDs: 153012995, 12149203, and 123608715 showed lower binding energies than the control (Ritonavir), indicating that they may become promising inhibitors for PLpro. MD was performed in a water solvent to validate the stability of the three complexes. All complexes achieved stable structure during the simulation as no significant fluctuations were observed in the validation parameters. Moreover, the binding energy for each complex was estimated using the MM-GBSA method. Complex 1 was the most stable structure based on the lowest binding energy score and its structure remained in a similar cavity with the docket snapshot. Based on our studies, three ligands were assumed to be potential inhibitors. The ligand of complex 1 may become the most promising antiviral drug against SARS-CoV-2 targeting PLpro.

Herbonanoceuticals: A Novel Beginning in Drug Discovery and Therapeutics

The Indian pharmaceutical industry is the world’s second largest industry (by volume) that develops products and market drugs licensed for use as medications. Medicines manufactured in the modern era are associated with major controversies such as non–target specificity, resistance, repeated administration, immune rejection, and other adverse effects on the body. Thus, there is a great need to find drugs that do not raise the aforementioned issues. Nature is an excellent hub providing a diverse range of phytoconstituents that open the way to phototherapeutics, which need a scientific path to deliver the active elements in a supported way to increase patient compliance and reduce the need for repeated administration. To discover a novel phytochemical as a lead compound for a therapeutic purpose is a real challenge. In former times, drug discovery was a complex process, as it took several years to find a lead compound for use against a particular disease. Nowadays, however, virtual screening methods have been developed, which are target specific, time consuming, and cost effective. To avoid increased and repeated administration of a drug, nanosized drug delivery systems for herbal drugs have been developed to enhance the activity and overcome problems associated with synthetic medicines. This review summarizes three main fields: drug discovery, docking for drug design, and last—but not least—drug delivery systems. Nowadays, nanobased drug delivery systems are in demand for delivery of herbal medicines used for therapeutic purposes. Herbonanoceuticals—herbal drugs of a nanosize—have better remedial value and fewer detrimental effects than modern medicines. Therefore, herbonanoceuticals can be a boon in the field of therapeutics.

Using a Consensus Docking Approach to Predict Adverse Drug Reactions in Combination Drug Therapies for Gulf War Illness

Gulf War Illness (GWI) is a chronic multisymptom illness characterized by fatigue, musculoskeletal pain, and gastrointestinal and cognitive dysfunction believed to stem from chemical exposures during the 1990⁻1991 Persian Gulf War. There are currently no treatments; however, previous studies have predicted a putative multi-intervention treatment composed of inhibiting Th1 immune cytokines followed by inhibition of the glucocorticoid receptor (GCR) to treat GWI. These predictions suggest the use of specific monoclonal antibodies or suramin to target interleukin-2 and tumor necrosis factor α , followed by mifepristone to inhibit the GCR. In addition to this putative treatment strategy, there exist a variety of medications that target GWI symptomatology. As pharmaceuticals are promiscuous molecules, binding to multiple sites beyond their intended targets, leading to off-target interactions, it is key to ensure that none of these medications interfere with the proposed treatment avenue. Here, we used the drug docking programs AutoDock 4.2, AutoDock Vina, and Schrödinger's Glide to assess the potential off-target immune and hormone interactions of 43 FDA-approved drugs commonly used to treat GWI symptoms in order to determine their putative polypharmacology and minimize adverse drug effects in a combined pharmaceutical treatment. Several of these FDA-approved drugs were predicted to be novel binders of immune and hormonal targets, suggesting caution for their use in the proposed GWI treatment strategy symptoms.

Synthesis, X-ray Single Crystal Structure, Molecular Docking and DFT Computations on N-[(1E)-1-(2H-1,3-Benzodioxol-5-yl)-3-(1H-imidazol-1-yl)propylidene]-hydroxylamine: A New Potential Antifungal Agent Precursor

Mycoses are serious health problem, especially in immunocompromised individuals. A new imidazole-bearing compound containing an oxime functionality was synthesized and characterized with different spectroscopic techniques to be used for the preparation of new antifungal agents. The stereochemistry of the oxime double bond was unequivocally determined via the single crystal X-ray technique. The title compound 4, C₁₃H₁₃N₃O₃·C₃H₈O, crystallizes in the monoclinic space group P2₁with a = 9.0963(3) Å, b = 14.7244(6) Å, c = 10.7035(4) Å, β = 94.298 (3)°, V = 1429.57(9) ų, Z = 2. The molecules were packed in the crystal structure by eight intermolecular hydrogen bond interactions. A comprehensive spectral analysis of the title molecule 4 has been performed based on the scaled quantum mechanical (SQM) force field obtained by density-functional theory (DFT) calculations. A molecular docking study illustrated the binding mode of the title compound 4 into its target protein. The preliminary antifungal activity of the title compound 4 was determined using a broth microdilution assay.

Low potency toxins reveal dense interaction networks in metabolism

BACKGROUND

The chemicals of metabolism are constructed of a small set of atoms and bonds. This may be because chemical structures outside the chemical space in which life operates are incompatible with biochemistry, or because mechanisms to make or utilize such excluded structures has not evolved. In this paper I address the extent to which biochemistry is restricted to a small fraction of the chemical space of possible chemicals, a restricted subset that I call Biochemical Space. I explore evidence that this restriction is at least in part due to selection again specific structures, and suggest a mechanism by which this occurs.

RESULTS

Chemicals that contain structures that our outside Biochemical Space (UnBiological groups) are more likely to be toxic to a wide range of organisms, even though they have no specifically toxic groups and no obvious mechanism of toxicity. This correlation of UnBiological with toxicity is stronger for low potency (millimolar) toxins. I relate this to the observation that most chemicals interact with many biological structures at low millimolar toxicity. I hypothesise that life has to select its components not only to have a specific set of functions but also to avoid interactions with all the other components of life that might degrade their function.

CONCLUSIONS

The chemistry of life has to form a dense, self-consistent network of chemical structures, and cannot easily be arbitrarily extended. The toxicity of arbitrary chemicals is a reflection of the disruption to that network occasioned by trying to insert a chemical into it without also selecting all the other components to tolerate that chemical. This suggests new ways to test for the toxicity of chemicals, and that engineering organisms to make high concentrations of materials such as chemical precursors or fuels may require more substantial engineering than just of the synthetic pathways involved.

Consensus docking: improving the reliability of docking in a virtual screening context

Structure-based virtual screening relies on scoring the predicted binding modes of compounds docked into the target. Because the accuracy of this scoring relies on the accuracy of the docking, methods that increase docking accuracy are valuable. Here, we present a relatively straightforward method for improving the probability of identifying accurately docked poses. The method is similar in concept to consensus scoring schemes, which have been shown to increase ranking power and thus hit rates, but combines information about predicted binding modes rather than predicted binding affinities. The pose prediction success rate of each docking program alone was found in this trial to be 55% for Autodock, 58% for DOCK, and 64% for Vina. By using more than one docking program to predict the binding pose, correct poses were identified in 82% or more of cases, a significant improvement. In a virtual screen, these more reliably posed compounds can be preferentially advanced to subsequent scoring stages to improve hit rates. Consensus docking can be easily introduced into established structure-based virtual screening methodologies.

Overview of high-throughput screening

High-throughput screening (HTS) is a key process used in drug discovery to identify hits from compound libraries that may become leads for medicinal chemistry optimization. This updated overview discusses the utilization of compound libraries, compounds derived from combinatorial and parallel synthesis campaigns and natural product sources; creation of mother and daughter plates; and compound storage, handling, and bar coding in HTS. The unit also presents an overview of established and emerging assay technologies (i.e., time-resolved fluorescence, fluorescence polarization, fluorescence-correlation spectroscopy, functional whole cell assays, and high-content assays) and their integration in automation hardware and IT systems. This revised unit provides updated descriptions of state-of-the-art instrumentation and technologies in this rapidly changing environment. The section on assay methodologies now also covers enzyme complementation assays and methods for high-throughput screening of ion channel activities. Finally, a section on criteria for assay robustness is included discussing the Z'-factor, which is now a widely accepted criterion for evaluation and validation of high throughput screening assays.

Integrating medicinal chemistry, organic/combinatorial chemistry, and computational chemistry for the discovery of selective estrogen receptor modulators with Forecaster, a novel platform for drug discovery

As part of a large medicinal chemistry program, we wish to develop novel selective estrogen receptor modulators (SERMs) as potential breast cancer treatments using a combination of experimental and computational approaches. However, one of the remaining difficulties nowadays is to fully integrate computational (i.e., virtual, theoretical) and medicinal (i.e., experimental, intuitive) chemistry to take advantage of the full potential of both. For this purpose, we have developed a Web-based platform, Forecaster, and a number of programs (e.g., Prepare, React, Select) with the aim of combining computational chemistry and medicinal chemistry expertise to facilitate drug discovery and development and more specifically to integrate synthesis into computer-aided drug design. In our quest for potent SERMs, this platform was used to build virtual combinatorial libraries, filter and extract a highly diverse library from the NCI database, and dock them to the estrogen receptor (ER), with all of these steps being fully automated by computational chemists for use by medicinal chemists. As a result, virtual screening of a diverse library seeded with active compounds followed by a search for analogs yielded an enrichment factor of 129, with 98% of the seeded active compounds recovered, while the screening of a designed virtual combinatorial library including known actives yielded an area under the receiver operating characteristic (AU-ROC) of 0.78. The lead optimization proved less successful, further demonstrating the challenge to simulate structure activity relationship studies.

POSE SCALING: GEOMETRICAL ASSESSMENT OF LIGAND BINDING POSES

A descriptor, the pose scaling factor, is proposed to quantitatively evaluate the geometrical match between a ligand and a target binding site. The pose scaling factor can be used to readily rank results of target-based in silico database screening or docking on large numbers of compounds. Such an approach will be of utility in the development and refinement of docking algorithms.

A Discovery Funnel for Nucleic Acid Binding Drug Candidates

Computational approaches are becoming increasingly popular for the discovery of drug candidates against a target of interest. Proteins have historically been the primary targets of many virtual screening efforts. While in silico screens targeting proteins has proven successful, other classes of targets, in particular DNA, remain largely unexplored using virtual screening methods. With the realization of the functional importance of many non-cannonical DNA structures such as G-quadruplexes, increased efforts are underway to discover new small molecules that can bind selectively to DNA structures. Here, we describe efforts to build an integrated in silico and in vitro platform for discovering compounds that may bind to a chosen DNA target. Millions of compounds are initially screened in silico for selective binding to a particular structure and ranked to identify several hundred best hits. An important element of our strategy is the inclusion of an array of possible competing structures in the in silico screen. The best hundred or so hits are validated experimentally for binding to the actual target structure by a high-throughput 96-well thermal denaturation assay to yield the top ten candidates. Finally, these most promising candidates are thoroughly characterized for binding to their DNA target by rigorous biophysical methods, including isothermal titration calorimetry, differential scanning calorimetry, spectroscopy and competition dialysis.This platform was validated using quadruplex DNA as a target and a newly discovered quadruplex binding compound with possible anti-cancer activity was discovered. Some considerations when embarking on virtual screening and in silico experiments are also discussed.

Design of new plasmepsin inhibitors: a virtual high throughput screening approach on the EGEE grid

Though different species of the genus Plasmodium may be responsible for malaria, the variant caused by P. falciparum is often very dangerous and even fatal if untreated. Hemoglobin degradation is one of the key metabolic processes for the survival of the Plasmodium parasite in its host. Plasmepsins, a family of aspartic proteases encoded by the Plasmodium genome, play a prominent role in host hemoglobin cleavage. In this paper we demonstrate the use of virtual screening, in particular molecular docking, employed at a very large scale to identify novel inhibitors for plasmepsins II and IV. A large grid infrastructure, the EGEE grid, was used to address the problem of large computation resources required for docking hundreds of thousands of chemical compounds on different plasmepsin targets of P. falciparum. A large compound library of about 1 million chemical compounds was docked on 5 different targets of plasmepsins using two different docking software, namely FlexX and AutoDock. Several strategies were employed to analyze the results of this virtual screening approach including docking scores, ideal binding modes, and interactions to key residues of the protein. Three different classes of structures with thiourea, diphenylurea, and guanidino scaffolds were identified to be promising hits. While the identification of diphenylurea compounds is in accordance with the literature and thus provides a sort of "positive control", the identification of novel compounds with a guanidino scaffold proves that high throughput docking can be effectively used to identify novel potential inhibitors of P. falciparum plasmepsins. Thus, with the work presented here, we do not only demonstrate the relevance of computational grids in drug discovery but also identify several promising small molecules which have the potential to serve as candidate inhibitors for P. falciparum plasmepsins. With the use of the EGEE grid infrastructure for the virtual screening campaign against the malaria causing parasite P. falciparum we have demonstrated that resource sharing on an eScience infrastructure such as EGEE provides a new model for doing collaborative research to fight diseases of the poor.

Discovery of novel DNA gyrase inhibitors by high-throughput virtual screening

The bacterial type II topoisomerases DNA gyrase and topoisomerase IV are validated targets for clinically useful quinolone antimicrobial drugs. A significant limitation to widely utilized quinolone inhibitors is the emergence of drug-resistant bacteria due to an altered DNA gyrase. To address this problem, we have used structure-based molecular docking to identify novel drug-like small molecules that target sites distinct from those targeted by quinolone inhibitors. A chemical ligand database containing approximately 140,000 small molecules (molecular weight, <500) was molecularly docked onto two sites of Escherichia coli DNA gyrase targeting (i) a previously unexplored structural pocket formed at the dimer interface of subunit A and (ii) a small region of the ATP binding pocket on subunit B overlapping the site targeted by coumarin and cyclothialidine drugs. This approach identified several small-molecule compounds that inhibited the DNA supercoiling activity of purified E. coli DNA gyrase. These compounds are structurally unrelated to previously identified gyrase inhibitors and represent potential scaffolds for the optimization of novel antibacterial agents that act on fluoroquinolone-resistant strains.

Scoring functions and enrichment: a case study on Hsp90

Background

The need for fast and accurate scoring functions has been driven by the increased use of in silico virtual screening twinned with high-throughput screening as a method to rapidly identify potential candidates in the early stages of drug development. We examine the ability of some the most common scoring functions (GOLD, ChemScore, DOCK, PMF, BLEEP and Consensus) to discriminate correctly and efficiently between active and non-active compounds among a library of ~3,600 diverse decoy compounds in a virtual screening experiment against heat shock protein 90 (Hsp90).

Results

Firstly, we investigated two ranking methodologies, GOLD_rank and BestScore_rank. GOLD_rank is based on ranks generated using GOLD. The various scoring functions, GOLD, ChemScore, DOCK, PMF, BLEEP and Consensus, are applied to the pose ranked number one by GOLD for that ligand. BestScore_rank uses multiple poses for each ligand and independently chooses the best ranked pose of the ligand according to each different scoring function. Secondly, we considered the effect of introducing the Thr184 hydrogen bond tether to guide the docking process towards a particular solution, and its effect on enrichment. Thirdly, we considered normalisation to account for the known bias of scoring functions to select larger molecules. All the scoring functions gave fairly similar enrichments, with the exception of PMF which was consistently the poorest performer. In most cases, GOLD was marginally the best performing individual function; the Consensus score usually performed similarly to the best single scoring function. Our best results were obtained using the Thr184 tether in combination with the BestScore_rank protocol and normalisation for molecular weight. For that particular combination, DOCK was the best individual function; DOCK recovered 90% of the actives in the top 10% of the ranked list; Consensus similarly recovered 89% of the actives in its top 10%.

Conclusion

Overall, we demonstrate the validity of virtual screening as a method for identifying new leads from a pool of ligands with similar physicochemical properties and we believe that the outcome of this study provides useful insight into the setting up of a suitable docking and scoring protocol, resulting in enrichment of 'target active' compounds.

A Common Reference Framework for Analyzing/Comparing Proteins and Ligands. Fingerprints for Ligands And Proteins (FLAP): Theory and Application

A fast new algorithm (Fingerprints for Ligands And Proteins or FLAP) able to describe small molecules and protein structures using a common reference framework of four-point pharmacophore fingerprints and a molecular-cavity shape is described in detail. The procedure starts by using the GRID force field to calculate molecular interaction fields, which are then used to identify particular target locations where an energetic interaction with small molecular features would be very favorable. The target points thus calculated are then used by FLAP to build all possible four-point pharmacophores present in the given target site. A related approach can be applied to small molecules, using directly the GRID atom types to identify pharmacophoric features, and this complementary description of the target and ligand then leads to several novel applications. FLAP can be used for selectivity studies or similarity analyses in order to compare macromolecules without superposing them. Protein families can be compared and clustered into target classes, without bias from previous knowledge and without requiring protein superposition, alignment, or knowledge-based comparison. FLAP can be used effectively for ligand-based virtual screening and structure-based virtual screening, with the pharmacophore molecular recognition. Finally, the new method can calculate descriptors for chemometric analysis and can initiate a docking procedure. This paper presents the background to the new procedure and includes case studies illustrating several relevant applications of the new approach.

PSI-DOCK: towards highly efficient and accurate flexible ligand docking

We have developed a new docking method, Pose-Sensitive Inclined (PSI)-DOCK, for flexible ligand docking. An improved SCORE function has been developed and used in PSI-DOCK for binding free energy evaluation. The improved SCORE function was able to reproduce the absolute binding free energies of a training set of 200 protein-ligand complexes with a correlation coefficient of 0.788 and a standard error of 8.13 kJ/mol. For ligand binding pose exploration, a unique searching strategy was designed in PSI-DOCK. In the first step, a tabu-enhanced genetic algorithm with a rapid shape-complementary scoring function is used to roughly explore and store potential binding poses of the ligand. Then, these predicted binding poses are optimized and compete against each other by using a genetic algorithm with the accurate SCORE function to determine the binding pose with the lowest docking energy. The PSI-DOCK 1.0 program is highly efficient in identifying the experimental binding pose. For a test dataset of 194 complexes, PSI-DOCK 1.0 achieved a 67% success rate (RMSD < 2.0 A) for only one run and a 74% success rate for 10 runs. PSI-DOCK can also predict the docking binding free energy with high accuracy. For a test set of 64 complexes, the correlation between the experimentally observed binding free energies and the docking binding free energies for 64 complexes is r = 0.777 with a standard deviation of 7.96 kJ/mol. Moreover, compared with other docking methods, PSI-DOCK 1.0 is extremely easy to use and requires minimum docking preparations. There is no requirement for the users to add hydrogen atoms to proteins because all protein hydrogen atoms and the flexibility of the terminal protein atoms are intrinsically taken into account in PSI-DOCK. There is also no requirement for the users to calculate partial atomic charges because PSI-DOCK does not calculate an electrostatic energy term. These features are not only convenient for the users but also help to avoid the influence of different preparation methods.

Screening Drug-Like Compounds by Docking to Homology Models: A Systematic Study

In the absence of an experimentally solved structure, a homology model of a protein target can be used instead for virtual screening of drug candidates by docking and scoring. This approach poses a number of questions regarding the choice of the template to use in constructing the model, the accuracy of the screening results, and the importance of allowing for protein flexibility. The present study addresses such questions with compound screening calculations for multiple homology models of five drug targets. A central result is that docking to homology models frequently yields enrichments of known ligands as good as that obtained by docking to a crystal structure of the actual target protein. Interestingly, however, standard measures of the similarity of the template used to build the homology model to the targeted protein show little correlation with the effectiveness of the screening calculations, and docking to the template itself often is as successful as docking to the corresponding homology model. Treating key side chains as mobile produces a modest improvement in the results. The reasons for these sometimes unexpected results, and their implications for future methodologic development, are discussed.

Unsupervised guided docking of covalently bound ligands

An approach for docking covalently bound ligands in protein enzymes or receptors was implemented in MacDOCK, a similarity-driven docking program based on DOCK 4.0. This approach was tested with a small number of covalent ligand-protein structures, using both native and non-native protein structures. In all cases, MacDOCK was able to generate orientations consistent with the known covalent binding mode of these complexes, with a performance similar to that of other docking programs. This method was also applied to search for known covalent thrombin inhibitors in a medium-sized molecular database (ca. 11,000 compounds). Detection of functional groups suitable for covalent docking was carried out automatically. A significant enrichment in known active molecules in the first 5% of the database was obtained, showing that MacDOCK can be used efficiently for the virtual screening of covalently bound ligands.

Virtual docking approaches to protein kinase B inhibition

We examined some in silico approaches to identify Akt (protein kinase B) inhibitors. Experimental validation of selected compounds was achieved using a fluorescence-based enzymatic assay and a substrate phosphorylation assay involving the protein GSK-3. We report on success and failure obtained by using several strategies including FlexX, GOLD, and CSCORE, where the 100-200 top-scoring compounds from a 50000-compound library were experimentally tested. This study led to the identification of low micromolar Akt1 inhibitors.

Docking and scoring in virtual screening for drug discovery: methods and applications

Computational approaches that 'dock' small molecules into the structures of macromolecular targets and 'score' their potential complementarity to binding sites are widely used in hit identification and lead optimization. Indeed, there are now a number of drugs whose development was heavily influenced by or based on structure-based design and screening strategies, such as HIV protease inhibitors. Nevertheless, there remain significant challenges in the application of these approaches, in particular in relation to current scoring schemes. Here, we review key concepts and specific features of small-molecule-protein docking methods, highlight selected applications and discuss recent advances that aim to address the acknowledged limitations of established approaches.

Validation of automated docking programs for docking and database screening against RNA drug targets

The increasing awareness of the essential role of RNA in controlling viral replication and in bacterial protein synthesis emphasizes the potential of ribonucleoproteins as targets for developing new antibacterial and antiviral drugs. RNA forms well defined three-dimensional structures with clefts and binding pockets reminiscent of the active sites of proteins. Furthermore, it precedes proteins in the translation pathway; inhibiting the function of a single RNA molecule would result in inhibition of multiple proteins. Thus, small molecules that bind RNA specifically would combine the advantages of antisense and RNAi strategies with the much more favorable medicinal chemistry of small-molecule therapeutics. The discovery of small-molecule inhibitors of RNA with attractive pharmacological potential would be facilitated if we had available effective computational tools of structure-based drug design. Here, we systematically test automated docking tools developed for proteins using existing three-dimensional structures of RNA-small molecule complexes. The results show that the native structures can generally be reproduced to within 2.5 angstroms more than 50-60% of the time. For more than half of the test complexes, the native ligand ranked among the top 10% compounds in a database-scoring test. Through this work, we provide parameters for the validated application of automated docking tools to the discovery of new inhibitors of RNA function.

Identification of Non-Phosphate-Containing Small Molecular Weight Inhibitors of the Tyrosine Kinase p56 Lck SH2 Domain via in Silico Screening against the pY + 3 Binding Site

The protein p56 lymphoid T cell tyrosine kinase (Lck) is predominantly expressed in T lymphocytes where it plays a critical role in T-cell-mediated immune response. Lck participates in phosphotyrosine-dependent protein-protein interactions through its modular binding unit, the Src homology-2 (SH2) domain. Accordingly, virtual screening methods combined with experimental assays were used to identify small molecular weight nonpeptidic compounds that block Lck SH2 domain-dependent interactions. Virtual screening included scoring normalization procedures and postdocking structural clustering that is shown to facilitate the selection of active compounds. By targeting the well-defined hydrophobic binding pocket known to impart specificity on Lck-protein interactions (i.e., pY + 3 site), inhibitors of the Lck SH2 domain were discovered that omit the phosphotyrosine (pY) or related moieties. The 34 out of 196 computationally selected compounds were shown to inhibit Lck SH2 domain association with phosphorylated immunoreceptor tyrosine based activation motifs peptide. Twenty-four of the active compounds were further tested for their ability to modulate biological function. Thirteen of these compounds showed inhibitory activity in mixed lymphocyte culture assay. Fluorescence titration experiments on four of these active compounds further verified their binding to the SH2 domain. Because of their simple chemical structures, these small organic compounds have the potential to act as lead compounds for the development of novel immunosuppressant drugs.

Virtual Screening Using Protein−Ligand Docking: Avoiding Artificial Enrichment

This study addresses a number of topical issues around the use of protein-ligand docking in virtual screening. We show that, for the validation of such methods, it is key to use focused libraries (containing compounds with one-dimensional properties, similar to the actives), rather than "random" or "drug-like" libraries to test the actives against. We also show that, to obtain good enrichments, the docking program needs to produce reliable binding modes. We demonstrate how pharmacophores can be used to guide the dockings and improve enrichments, and we compare the performance of three consensus-ranking protocols against ranking based on individual scoring functions. Finally, we show that protein-ligand docking can be an effective aid in the screening for weak, fragment-like binders, which has rapidly become a popular strategy for hit identification. All results presented are based on carefully constructed virtual screening experiments against four targets, using the protein-ligand docking program GOLD.

FDS: flexible ligand and receptor docking with a continuum solvent model and soft-core energy function

The docking of flexible small molecule ligands to large flexible protein targets is addressed in this article using a two-stage simulation-based method. The methodology presented is a hybrid approach where the first component is a dock of the ligand to the protein binding site, based on deriving sets of simultaneously satisfied intermolecular hydrogen bonds using graph theory and a recursive distance geometry algorithm. The output structures are reduced in number by cluster analysis based on distance similarities. These structures are submitted to a modified Monte Carlo algorithm using the AMBER-AA molecular mechanics force field with the Generalized Born/Surface Area (GB/SA) continuum model. This solvent model is not only less expensive than an explicit representation, but also yields increased sampling. Sampling is also increased using a rotamer library to direct some of the protein side-chain movements along with large dihedral moves. Finally, a softening function for the nonbonded force field terms is used, enabling the potential energy function to be slowly turned on throughout the course of the simulation. The docking procedure is optimized, and the results are presented for a single complex of the arabinose binding protein. It was found that for a rigid receptor model, the X-ray binding geometry was reproduced and uniquely identified based on the associated potential energy. However, when side-chain flexibility was included, although the X-ray structure was identified, it was one of three possible binding geometries that were energetically indistinguishable. These results suggest that on relaxing the constraint on receptor flexibility, the docking energy hypersurface changes from being funnel-like to rugged. A further 14 complexes were then examined using the optimized protocol. For each complex the docking methodology was tested for a fully flexible ligand, both with and without protein side-chain flexibility. For the rigid protein docking, 13 out of the 15 test cases were able to find the experimental binding mode; this number was reduced to 11 for the flexible protein docking. However, of these 11, in the majority of cases the experimental binding mode was not uniquely identified, but was present in a cluster of low energy structures that were energetically indistinguishable. These results not only support the presence of a rugged docking energy hypersurface, but also suggest that it may be necessary to consider the possibility of more than one binding conformation during ligand optimization.

A web-based platform for virtual screening

A fully integrated, web-based, virtual screening platform has been developed to allow rapid virtual screening of large numbers of compounds. ORACLE is used to store information at all stages of the process. The system includes a large database of historical compounds from high throughput screenings (HTS) chemical suppliers, ATLAS, containing over 3.1 million unique compounds with their associated physiochemical properties (ClogP, MW, etc.). The database can be screened using a web-based interface to produce compound subsets for virtual screening or virtual library (VL) enumeration. In order to carry out the latter task within ORACLE a reaction data cartridge has been developed. Virtual libraries can be enumerated rapidly using the web-based interface to the cartridge. The compound subsets can be seamlessly submitted for virtual screening experiments, and the results can be viewed via another web-based interface allowing ad hoc querying of the virtual screening data stored in ORACLE.

Analysis and optimization of structure-based virtual screening protocols. 2. Examination of docked ligand orientation sampling methodology: mapping a pharmacophore for success

An important element of any structure-based virtual screening (SVS) technique is the method used to orient the ligands in the target active site. This has been a somewhat overlooked issue in recent SVS validation studies, with the assumption being made that the performance of an algorithm for a given set of orientation sampling settings will be representative for the general behavior of said technique. Here, we analyze five different SVS targets using a variety of sampling paradigms within the DOCK, GOLD and PROMETHEUS programs over a data set of approximately 10,000 noise compounds, combined with data sets containing multiple active compounds. These sets have been broken down by chemotype, with chemotype hit rate used to provide a measure of enrichment with a potentially improved relevance to real world SVS experiments. The variability in enrichment results produced by different sampling paradigms is illustrated, as is the utility of using pharmacophores to constrain sampling to regions that reflect known structural biology. The difference in results when comparing chemotype with compound hit rates is also highlighted.

Improved protein-ligand docking using GOLD

The Chemscore function was implemented as a scoring function for the protein-ligand docking program GOLD, and its performance compared to the original Goldscore function and two consensus docking protocols, "Goldscore-CS" and "Chemscore-GS," in terms of docking accuracy, prediction of binding affinities, and speed. In the "Goldscore-CS" protocol, dockings produced with the Goldscore function are scored and ranked with the Chemscore function; in the "Chemscore-GS" protocol, dockings produced with the Chemscore function are scored and ranked with the Goldscore function. Comparisons were made for a "clean" set of 224 protein-ligand complexes, and for two subsets of this set, one for which the ligands are "drug-like," the other for which they are "fragment-like." For "drug-like" and "fragment-like" ligands, the docking accuracies obtained with Chemscore and Goldscore functions are similar. For larger ligands, Goldscore gives superior results. Docking with the Chemscore function is up to three times faster than docking with the Goldscore function. Both combined docking protocols give significant improvements in docking accuracy over the use of the Goldscore or Chemscore function alone. "Goldscore-CS" gives success rates of up to 81% (top-ranked GOLD solution within 2.0 A of the experimental binding mode) for the "clean list," but at the cost of long search times. For most virtual screening applications, "Chemscore-GS" seems optimal; search settings that give docking speeds of around 0.25-1.3 min/compound have success rates of about 78% for "drug-like" compounds and 85% for "fragment-like" compounds. In terms of producing binding energy estimates, the Goldscore function appears to perform better than the Chemscore function and the two consensus protocols, particularly for faster search settings. Even at docking speeds of around 1-2 min/compound, the Goldscore function predicts binding energies with a standard deviation of approximately 10.5 kJ/mol.

Analysis and optimization of structure-based virtual screening protocols (1): exploration of ligand conformational sampling techniques

Ligand conformational flexibility has long been recognized as an important issue in virtual screening (VS). To this end, a number of different methodologies have been adapted to tackle the problem. Many of said techniques were originally designed for ligand derived pharmacophore screens, but have subsequently been fashioned for application within structure-based virtual screening (SVS). A popular adaptation is the pre-calculation of diverse ligand conformations for subsequent docking in target active sites. In this paper, we study a number of the software programs currently being used in conformer generation, analyzing their ability to regenerate known ligand binding conformations. The implications of these studies are discussed, from the perspective of VS in general and SVS in particular.

Analysis and optimization of structure-based virtual screening protocols. (3). New methods and old problems in scoring function design

Scoring function research remains a primary focus of current structure-based virtual screening (SVS) technology development. Here, we present an alternative method for scoring function design that attempts to combine crystallographic structural information with data derived from directly within SVS calculations. The technique utilizes a genetic algorithm (GA) to optimize functions based on binding property data derived from multiple virtual screening calculations. These calculations are undertaken on protein data bank (PDB) complex active sites using ligands of known binding mode in conjunction with "noise" compounds. The advantages of such an approach are that the function does not rely on assay data and that it can potentially use the "noise" binding data to recognize the sub-optimal docking interactions inherent in SVS calculations. Initial efforts in technique exploration using DOCK are presented, with comparisons made to existing DOCK scoring functions. An analysis of the problems inherent to scoring function development is also made, including issues in dataset creation and limitations in descriptor utility when viewed from the perspective of docking mode resolution. The future directions such studies might take are also discussed in detail.

Virtual screening of virtual libraries

Virtual screening of virtual libraries (VSVL) is a rapidly changing area of research. Great efforts are being made to produce better algorithms, selection methods and infrastructure. Yet, the number of successful examples in the literature is not impressive, although the quality of work certainly is high. Why is this? One reason is that these methods tend to be applied at the lead generation stage and therefore there is a large lead-time before successful examples appear in the literature. However, any computational chemist would confirm that these methods are successful and there exists a glut of start-up companies specialising in virtual screening. Moreover, the scientific community would not be focussing so much attention on this area if it were not yielding results. Even so, the paucity of literature data is certainly a hindrance to the development of better methods. The VSVL process is unique within the discovery process, in that it is the only method that can screen the > 10(30) genuinely novel molecules out there. Already, some VSVL methods are evaluating 10(13) compounds, a capacity that high throughput screening can only dream of. There is a huge potential advantage for the company that develops efficient and effective methods, for lead generation, lead hopping and optimization of both potency and ADME properties. To do this, it requires more than the software, it requires confidence to exploit the methodology, to commit synthesis on the basis of it, and to build this approach into the medicinal chemistry strategy. It is a fact that these tools remain quite daunting for the majority of scientists working at the bench. The routine use of these methods is not simply a matter of education and training. Integration of these methods into accessible and robust end user software, without dilution of the science, must be a priority. We have reached a coincidence, where several technologies have the required level of maturity predictive computational chemistry methods, algorithms that manage the combinatorial explosion, high throughput crystallography and ADME measurements and the massive increase in computational horsepower from distributed computing. The author is confident that the synergy of these technologies will bring great benefit to the industry, with more efficient production of higher quality clinical candidates. The future is bright. The future is virtual!

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 approximately 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.

A new test set for validating predictions of protein-ligand interaction

We present a large test set of protein-ligand complexes for the purpose of validating algorithms that rely on the prediction of protein-ligand interactions. The set consists of 305 complexes with protonation states assigned by manual inspection. The following checks have been carried out to identify unsuitable entries in this set: (1) assessing the involvement of crystallographically related protein units in ligand binding; (2) identification of bad clashes between protein side chains and ligand; and (3) assessment of structural errors, and/or inconsistency of ligand placement with crystal structure electron density. In addition, the set has been pruned to assure diversity in terms of protein-ligand structures, and subsets are supplied for different protein-structure resolution ranges. A classification of the set by protein type is available. As an illustration, validation results are shown for GOLD and SuperStar. GOLD is a program that performs flexible protein-ligand docking, and SuperStar is used for the prediction of favorable interaction sites in proteins. The new CCDC/Astex test set is freely available to the scientific community (http://www.ccdc.cam.ac.uk).

Structure-based virtual screening: an overview

Enormous advances in genomics have resulted in a large increase in the number of potential therapeutic targets that are available for investigation. This growth in potential targets has increased the demand for reliable target validation, as well as technologies that can identify rapidly several quality lead candidates. Virtual screening, and in particular receptor-based virtual screening, has emerged as a reliable, inexpensive method for identifying leads. Although still an evolving method, advances in computational techniques have enabled virtual screening to have a positive impact on the discovery process. Here, the current strengths and weaknesses of the technology are discussed, and emphasis is placed on aspects of the work-flow of a virtual screening campaign, from preparation through to post-screening analysis.

A review of protein-small molecule docking methods

The binding of small molecule ligands to large protein targets is central to numerous biological processes. The accurate prediction of the binding modes between the ligand and protein, (the docking problem) is of fundamental importance in modern structure-based drug design. An overview of current docking techniques is presented with a description of applications including single docking experiments and the virtual screening of databases.

Modified AutoDock for accurate docking of protein kinase inhibitors

Protein kinases are an important class of enzymes controlling virtually all cellular signaling pathways. Consequently, selective inhibitors of protein kinases have attracted significant interest as potential new drugs for many diseases. Computational methods, including molecular docking, have increasingly been used in the inhibitor design process [1]. We have considered several docking packages in order to strengthen our kinase inhibitor work with computational capabilities. In our experience, AutoDock offered a reasonable combination of accuracy and speed, as opposed to methods that specialize either in fast database searches or detailed and computationally intensive calculations. However, AutoDock did not perform well in cases where extensive hydrophobic contacts were involved, such as docking of SB203580 to its target protein kinase p38. Another shortcoming was a hydrogen bonding energy function, which underestimated the attraction component and, thus, did not allow for sufficiently accurate modeling of the key hydrogen bonds in the kinase-inhibitor complexes. We have modified the parameter set used to model hydrogen bonds, which increased the accuracy of AutoDock and appeared to be generally applicable to many kinase-inhibitor pairs without customization. Binding to largely hydrophobic sites, such as the active site of p38, was significantly improved by introducing a correction factor selectively affecting only carbon and hydrogen energy grids, thus, providing an effective, although approximate, treatment of solvation.

Protein-ligand recognition using spherical harmonic molecular surfaces: towards a fast and efficient filter for large virtual throughput screening

Molecular surfaces are important because surface-shape complementarity is often a necessary condition in protein-ligand interactions and docking studies. We have previously described a fast and efficient method to obtain triangulated surface-meshes by topologically mapping ellipsoids on molecular surfaces. In this paper, we present an extension of our work to spherical harmonic surfaces in order to approximate molecular surfaces of both ligands and receptor-cavities and to easily check the surface-shape complementarity. The method consists of (1) finding lobes and holes on both ligand and cavity surfaces using contour maps of radius functions with spherical harmonic expansions, (2) superposing the surfaces around a given binding site by minimizing the distance between their respective expansion coefficients. This docking procedure capabilities was demonstrated by application to 35 protein-ligand complexes of known crystal structures. The method can also be easily and efficiently used as a filter to detect in a large conformational sampling the possible conformations presenting good complementarity with the receptor site, and being, therefore, good candidates for further more elaborate docking studies. This "virtual screening" was demonstrated on the platelet thrombin receptor.