Energetic analysis of fragment docking and application to structure-based pharmacophore hypothesis generation

We have developed a method that uses energetic analysis of structure-based fragment docking to elucidate key features for molecular recognition. This hybrid ligand- and structure-based methodology uses an atomic breakdown of the energy terms from the Glide XP scoring function to locate key pharmacophoric features from the docked fragments. First, we show that Glide accurately docks fragments, producing a root mean squared deviation (RMSD) of <1.0 Å for the top scoring pose to the native crystal structure. We then describe fragment-specific docking settings developed to generate poses that explore every pocket of a binding site while maintaining the docking accuracy of the top scoring pose. Next, we describe how the energy terms from the Glide XP scoring function are mapped onto pharmacophore sites from the docked fragments in order to rank their importance for binding. Using this energetic analysis we show that the most energetically favorable pharmacophore sites are consistent with features from known tight binding compounds. Finally, we describe a method to use the energetically selected sites from fragment docking to develop a pharmacophore hypothesis that can be used in virtual database screening to retrieve diverse compounds. We find that this method produces viable hypotheses that are consistent with known active compounds. In addition to retrieving diverse compounds that are not biased by the co-crystallized ligand, the method is able to recover known active compounds from a database screen, with an average enrichment of 8.1 in the top 1% of the database.     

Evaluating docking programs: keeping the playing field level

Over recent years many enrichment studies have been published which purport to rigorously compare the performance of two or more docking protocols. It has become clear however that such studies often have flaws within their methodologies, which cast doubt on the rigour of the conclusions. Setting up such comparisons is fraught with difficulties and no best mode of practice is available to guide the experimenter. Careful choice of structural models and ligands appropriate to those models is important. The protein structure should be representative for the target. In addition the set of active ligands selected should be appropriate to the structure in cases where different forms of the protein bind different classes of ligand. Binding site definition is also an area in which errors arise. Particular care is needed in deciding which crystallographic waters to retain and again this may be predicated by knowledge of the likely binding modes of the ligands making up the active ligand list. Geometric integrity of the ligand structures used is clearly important yet it is apparent that published sets of actives + decoys may contain sometimes high proportions of incorrect structures. Choice of protocol for docking and analysis needs careful consideration as many programs can be tweaked for optimum performance. Should studies be run using ‘black box’ protocols supplied by the software provider? Lastly, the correct method of analysis of enrichment studies is a much discussed topic at the moment. However currently promoted approaches do not consider a crucial aspect of a successful virtual screen, namely that a good structural diversity of hits be returned. Overall there is much to consider in the experimental design of enrichment studies. Hopefully this study will be of benefit in helping others plan such experiments.     

Pharmacophore modeling,virtual screening and molecular docking of ATPase inhibitors of HSP70

Heat shock protein 70 is an effective anticancer target as it influences many signaling pathways. Hence the study investigated the important pharmacophore feature required for ATPase inhibitors of HSP70 by generating a ligand based pharmacophore model followed by virtual based screening and subsequent validation by molecular docking in Discovery studio V4.0. The most extrapolative pharmacophore model (hypotheses 8) consisted of four hydrogen bond acceptors. Further validation by external test set prediction identified 200 hits from Mini Maybridge, Drug Diverse, SCPDB compounds and Phytochemicals. Consequently, the screened compounds were refined by rule of five, ADMET and molecular docking to retain the best competitive hits. Finally Phytochemical compounds Muricatetrocin B, Diacetylphiladelphicalactone C, Eleutheroside B and 5-(3-{[1-(benzylsulfonyl)piperidin-4-yl]amino}phenyl)- 4-bromo-3-(carboxymethoxy)thiophene-2-carboxylic acid were obtained as leads to inhibit the ATPase activity of HSP70 in our findings and thus can be proposed for further in vitro and in vivo evaluation.     

Improving database enrichment through ensemble docking

While it may seem intuitive that using an ensemble of multiple conformations of a receptor in structure-based virtual screening experiments would necessarily yield improved enrichment of actives relative to using just a single receptor, it turns out that at least in the p38 MAP kinase model system studied here, a very large majority of all possible ensembles do not yield improved enrichment of actives. However, there are combinations of receptor structures that do lead to improved enrichment results. We present here a method to select the ensembles that produce the best enrichments that does not rely on knowledge of active compounds or sophisticated analyses of the 3D receptor structures. In the system studied here, the small fraction of ensembles of up to 3 receptors that do yield good enrichments of actives were identified by selecting ensembles that have the best mean GlideScore for the top 1% of the docked ligands in a database screen of actives and drug-like "decoy" ligands. Ensembles of two receptors identified using this mean GlideScore metric generally outperform single receptors, while ensembles of three receptors identified using this metric consistently give optimal enrichment factors in which, for example, 40% of the known actives outrank all the other ligands in the database.     

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.     

Selective ATP competitive leads of CDK4: Discovery by 3D-QSAR pharmacophore mapping and molecular docking approach

The discovery of ATP competitive CDK4 inhibitors is an on-going challenging task in cancer therapy. Here, an attempt has been made to develop new leads targeting ATP binding site of CDK4 by applying 3D-QSAR pharmacophore mapping and molecular docking methods The outcome of 6 leads offers a significant contribution for selective CDK4 inhibition, since they show potential binding interactions with Val⁹⁶, Arg¹⁰¹, and Glu¹⁴⁴ residues of CDK4, that are unique and from other kinases. It is worth noting that there is a striking similarity in binding interactions of the leads and known CDK4 inhibitors, namely Abemaciclib, Palbociclib and Ribociclib. Further key features, including high dock score value, good predicted activity, scaffold diversity, and the acceptable ADME profile of leads, provide a great opportunity for the development of highly potent and selective ATP competitive inhibitors of CDK4.     

Virtual fragment screening: an exploration of various docking and scoring protocols for fragments using Glide

Fragment-based drug discovery approaches allow for a greater coverage of chemical space and generally produce high efficiency ligands. As such, virtual and experimental fragment screening are increasingly being coupled in an effort to identify new leads for specific therapeutic targets. Fragment docking is employed to create target-focussed subset of compounds for testing along side generic fragment libraries. The utility of the program Glide with various scoring schemes for fragment docking is discussed. Fragment docking results for two test cases, prostaglandin D2 synthase and DNA ligase, are presented and compared to experimental screening data. Self-docking, cross-docking, and enrichment studies are performed. For the enrichment runs, experimental data exists indicating that the docking decoys in fact do not inhibit the corresponding enzyme being examined. Results indicate that even for difficult test cases fragment docking can yield enrichments significantly better than random. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Flexibases: A way to enhance the use of molecular docking methods

Summary Specially expanded databases containing three-dimensional structures are created to enhance the utility of docking methods to find new leads, i.e., active compounds of pharmacological interest. The expansion is based on the automatic generation of a set of maximally dissimilar conformations. The ligand receptor system of methotrexate and dihydrofolate reductase is used to demonstrate the feasibility of creating flexibases and their utility in docking studies.     

Recommendations for evaluation of computational methods

The field of computational chemistry, particularly as applied to drug design, has become increasingly important in terms of the practical application of predictive modeling to pharmaceutical research and development. Tools for exploiting protein structures or sets of ligands known to bind particular targets can be used for binding-mode prediction, virtual screening, and prediction of activity. A serious weakness within the field is a lack of standards with respect to quantitative evaluation of methods, data set preparation, and data set sharing. Our goal should be to report new methods or comparative evaluations of methods in a manner that supports decision making for practical applications. Here we propose a modest beginning, with recommendations for requirements on statistical reporting, requirements for data sharing, and best practices for benchmark preparation and usage.     

Surflex-Dock 2.1: Robust performance from ligand energetic modeling,ring flexibility,and knowledge-based search

The Surflex flexible molecular docking method has been generalized and extended in two primary areas related to the search component of docking. First, incorporation of a small-molecule force-field extends the search into Cartesian coordinates constrained by internal ligand energetics. Whereas previous versions searched only the alignment and acyclic torsional space of the ligand, the new approach supports dynamic ring flexibility and all-atom optimization of docked ligand poses. Second, knowledge of well established molecular interactions between ligand fragments and a target protein can be directly exploited to guide the search process. This offers advantages in some cases over the search strategy where ligand alignment is guided solely by a “protomol” (a pre-computed molecular representation of an idealized ligand). Results are presented on both docking accuracy and screening utility using multiple publicly available benchmark data sets that place Surflex’s performance in the context of other molecular docking methods. In terms of docking accuracy, Surflex-Dock 2.1 performs as well as the best available methods. In the area of screening utility, Surflex’s performance is extremely robust, and it is clearly superior to other methods within the set of cases for which comparative data are available, with roughly double the screening enrichment performance.     

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     

Prediction of new chromene-based inhibitors of tubulin using structure-based virtual screening and molecular dynamics simulation methods

Multidrug resistance (MDR) is one of the serious problems in cancer research that causes failure in chemotherapy. Chromene-based compounds have been proven to be the novel anti-MDR agents for inhibiting proliferation of tumor cells through tubulin polymerization inhibition of by binding at the colchicine binding site. In this study, we screened a chromene-based database of small molecules using physicochemical, ADMET properties and molecular docking to identify potential hit compounds. In order to validate our hit compounds, molecular dynamics simulations and related analysis were carried out and the results suggest that our hit compounds (PubChem CIDs: 16814409, 17594471, 57367244 and 69899719) can prove to be potential inhibitors of tubulin. The in silico results show that the present hits, like colchicine, effectively suppressed the dynamic instability of microtubules and induced microtubule-depolymerization and cell cycle arrest.     

VSDMIP: virtual screening data management on an integrated platform

A novel software (VSDMIP) for the virtual screening (VS) of chemical libraries integrated within a MySQL relational database is presented. Two main features make VSDMIP clearly distinguishable from other existing computational tools: (i) its database, which stores not only ligand information but also the results from every step in the VS process, and (ii) its modular and pluggable architecture, which allows customization of the VS stages (such as the programs used for conformer generation or docking), through the definition of a detailed workflow employing user-configurable XML files. VSDMIP, therefore, facilitates the storage and retrieval of VS results, easily adapts to the specific requirements of each method and tool used in the experiments, and allows the comparison of different VS methodologies. To validate the usefulness of VSDMIP as an automated tool for carrying out VS several experiments were run on six protein targets (acetylcholinesterase, cyclin-dependent kinase 2, coagulation factor Xa, estrogen receptor alpha, p38 MAP kinase, and neuraminidase) using nine binary (actives/inactive) test sets. The performance of several VS configurations was evaluated by means of enrichment factors and receiver operating characteristic plots. ángel R. Ortiz deceased on May 5, 2008.     

Dihydropyrazole and dihydropyrrole structures based design of Kif15 inhibitors as novel therapeutic agents for cancer

Mitotic Kinesin motors, Eg5 and Kif15, have recently emerged as good targets for cancer as they play an inevitable role during mitosis. But, most of the Eg5 inhibitors were found ineffective when the cancer cells develop resistance to them by escalating the expression of Kif15 as alternative to Eg5. Therefore, the drugs that target Kif15 became necessary to be used either as a single or in combination with Eg5 inhibitors. The present study used 39 dihydropyrazole and 13 dihydropyrrole derivatives that were having in vitro inhibitory potential against kinesin motors to develop a common pharmacophore hypothesis AHRR and atom-based QSAR model. The model was used for virtual screening of ZINC database and the resultant hits were docked against Kif15. The four drug candidates with high docking score were examined for their activity and pharmacokinetic behaviour. Based on the results these drugs could be considered as lead candidates in further drug development for cancer.     

Discovery of potential inhibitor against human acetylcholinesterase: a molecular docking and molecular dynamics investigation

Alzheimer’s disease (AD) is a progressive neurodegenerative disease of central nervous system among elderly people. Human acetylcholinesterase (hAChE), an important enzyme in neuronal signaling, is responsible for the degradation of acetylcholine which in turn prevents the post synaptic signal transmissions. hAChE has been an attractive target of drug discovery for the search of therapeutics against AD. In the recent past hAChE has become hot target for the investigation of new potential therapeutics. We performed virtual screening of entire database against hAChE. Further, the extra precision molecular docking was carried out to refine the docking results and the best complex was passed for molecular dynamics simulations in order of understanding the hAChE dynamics and its behavior in complex with the ligand which corroborate the outcomes of virtual screening. This also provides binding free energy data that establishes the ligands efficiency for inhibiting hAChE. The computational findings discussed in this paper provide initial information of inhibitory effects of ligand, (drugbank entry DB00983), over hAChE.     

Synthesis,characterization, computational,excited state properties,wave function and molecular docking studies of (E)-1-(perfluorophenyl)-N-(p-tolyl) methanimine

The compound (E)-1-(perfluorophenyl)-N-(p-tolyl)methanimine (PFPT) was synthesized and characterized by the Infrared, UV–Visible, and NMR analysis. Using density functional theory, the current work is a set of theoretical studies on PFPT. The compound molecular structure and geometry were defined using DFT. Topological studies, like electron localized function, localized orbital locator, average localized ionization energy, and reduced density gradient studies, were done with the Multiwfn-3.8 to find the main binding areas and weak interactions in the molecule. Using the IEFPCM solvation model used to study the calculated UV–Visible spectrum, we used two different solvents. The HOMO-LUMO, MEP, and NLO properties were carried out by DFT/B3LYP/cc-pVDZ in the gas phase. The NBO calculations are used to study how charges move between and within the molecule and the stability of this molecule. A pharmacological analysis is done using an online tool like Swiss-ADME, to see if the molecule could be a potential drug candidate; this evaluation looks at the drug-likeness, ADME, and eco-friendly toxicity properties of the PFPT molecule. Auto-dock suite and Discovery studio Visualizer are used to do molecular docking against 2QFA protein.     

Discovery of potential Zika virus RNA polymerase inhibitors by docking-based virtual screening

Zika virus (ZIKV) infection has been associated with Guillain-Barre syndrome in adults and microcephaly in infants. The existence of insufficient structural data in most of the protein databases hinders the synthesis of anti-ZIKV pharmaceutics. In this work, we attempted to model the catalytic domain of the ZIKV RNA polymerase (RdRpC) along with a detailed assessment of conserved aspartates in ZIKV RdRpC palm domain as potential drug targets. The conserved and catalytically active aspartate residues present in the predicted RdRpC protein were virtually screened against a ZINC database for inhibitors, and the selected potential drug candidates were further filtered based on their ADMET profiles. One of the pharmacokinetically active compounds (Ligand 6) showed a remarkable docking profile against the strictly conserved aspartate residues of the RdRpC active site. We hypothesize that the Ligand 6 may form a potential drug candidate for RdRpC inhibition in the clinical treatment of ZIKV infection.     

Virtual Screening and Structure Generation Applied to Drug Design

The methods of computer-aided drug design can be divided into two categories according to whether or not the structures of receptors are known1, corresponding to two principal strategies: (1) searching the bio-active ligands against virtual combinatorial libraries and calculating the affinity energy between ligand and receptor by docking ; (2) QSAR and 3D-structure data-mining. 3D-QSAR method is now applied widely to drug discovery, but this method is generally limited to refine the structu…