Customizable de novo design strategies for DOCK: Application to HIVgp41 and other therapeutic targets

De novo design can be used to explore vast areas of chemical space in computational lead discovery. As a complement to virtual screening, from‐scratch construction of molecules is not limited to compounds in pre‐existing vendor catalogs. Here, we present an iterative fragment growth method, integrated into the program DOCK, in which new molecules are built using rules for allowable connections based on known molecules. The method leverages DOCK's advanced scoring and pruning approaches and users can define very specific criteria in terms of properties or features to customize growth toward a particular region of chemical space. The code was validated using three increasingly difficult classes of calculations: (1) Rebuilding known X‐ray ligands taken from 663 complexes using only their component parts (focused libraries), (2) construction of new ligands in 57 drug target sites using a library derived from ∼13M drug‐like compounds (generic libraries), and (3) application to a challenging protein‐protein interface on the viral drug target HIVgp41. The computational testing confirms that the de novo DOCK routines are robust and working as envisioned, and the compelling results highlight the potential utility for designing new molecules against a wide variety of important protein targets. © 2017 Wiley Periodicals, Inc.     

Comparative performance of several flexible docking programs and scoring functions: enrichment studies for a diverse set of pharmaceutically relevant targets

Virtual screening by molecular docking has become a widely used approach to lead discovery in the pharmaceutical industry when a high-resolution structure of the biological target of interest is available. The performance of three widely used docking programs (Glide, GOLD, and DOCK) for virtual database screening is studied when they are applied to the same protein target and ligand set. Comparisons of the docking programs and scoring functions using a large and diverse data set of pharmaceutically interesting targets and active compounds are carried out. We focus on the problem of docking and scoring flexible compounds which are sterically capable of docking into a rigid conformation of the receptor. The Glide XP methodology is shown to consistently yield enrichments superior to the two alternative methods, while GOLD outperforms DOCK on average. The study also shows that docking into multiple receptor structures can decrease the docking error in screening a diverse set of active compounds.     

Evaluation of docking calculations on X-ray structures using CONSENSUS-DOCK

We are participating in the challenge of identifying active compounds for target proteins using structure-based virtual screening (SBVS). We use an in-house customized docking program, CONSENSUS-DOCK, which is a customized version of the DOCK4 program in which three scoring functions (DOCK4, FlexX and PMF) and consensus scoring have been implemented. This paper compares the docking calculation results obtained using CONSENSUS-DOCK and DOCK4, and demonstrates that CONSENSUS-DOCK produces better results than DOCK4 for major X-ray structures obtained from the Protein Data Bank (PDB).     

SPOT‐Ligand: Fast and effective structure‐based virtual screening by binding homology search according to ligand and receptor similarity

Structure‐based virtual screening usually involves docking of a library of chemical compounds onto the functional pocket of the target receptor so as to discover novel classes of ligands. However, the overall success rate remains low and screening a large library is computationally intensive. An alternative to this “ab initio” approach is virtual screening by binding homology search. In this approach, potential ligands are predicted based on similar interaction pairs (similarity in receptors and ligands). SPOT‐Ligand is an approach that integrates ligand similarity by Tanimoto coefficient and receptor similarity by protein structure alignment program SPalign. The method was found to yield a consistent performance in DUD and DUD‐E docking benchmarks even if model structures were employed. It improves over docking methods (DOCK6 and AUTODOCK Vina) and has a performance comparable to or better than other binding‐homology methods (FINDsite and PoLi) with higher computational efficiency. The server is available at http://sparks-lab.org . © 2016 Wiley Periodicals, Inc.     

Computer‐aided Discovery of Potential Inhibitors for Transthyretin‐related Amyloidosis

Transthyretin (TTR) is a homotetrameric plasma protein associated with human amyloid diseases. Although Tafamidis has been recently approved for the treatment of TTR familial amyloid polyneuropathy (FAP), there is still a need for more effective drugs in the treatment of TTR amyloidosis diseases. In this study, a computer‐aided approach combining molecular docking, virtual screening and molecular dynamics (MD) simulations was employed to identify potent TTR amyloidosis inhibitors from National Cancer Institute (NCI), Maybridge and Asdi databases. A receptor‐specific scoring function was also developed using comparative binding energy (COMBINE) method to accurately predict the inhibitory activities for the selected compounds during virtual screening. The developed receptor‐specific scoring function demonstrated good predictive ability by yielding strong correlation coefficients between experimental activities and estimated activities for 32 training set and 9 test set compounds, respectively. Moreover, it was successfully applied to rank the candidate compounds from structure‐based virtual screening. Finally, three compounds (NSC220163, MFCD00276817 and SPB06319) were identified as potential leads, which exhibited higher predicted inhibitory activities and higher binding affinities in comparison to the Tafamidis. Our results further suggest that halogen bonding interaction plays a crucial role in stabilizing the TTR‐inhibitor complex. These results indicate that our computational approach could effectively discover more potent TTR amyloidosis inhibitors, which can be further validate by in vitro and in vivo biological tests.     

Comparison of two implementations of the incremental construction algorithm in flexible docking of thrombin inhibitors

A set of 32 known thrombin inhibitors representing different chemical classes has been used to evaluate the performance of two implementations of incremental construction algorithms for flexible molecular docking: DOCK 4.0 and FlexX 1.5. Both docking tools are able to dock 10–35% of our test set within 2 Å of their known, bound conformations using default sampling and scoring parameters. Although flexible docking with DOCK or FlexX is not able to reconstruct all native complexes, it does offer a significant improvement over rigid body docking of single, rule-based conformations, which is still often used for docking of large databases. Docking of sets of multiple conformers of each inhibitor, obtained with a novel protocol for diverse conformer generation and selection, yielded results comparable to those obtained by flexible docking. Chemical scoring, which is an empirically modified force field scoring method implemented in DOCK 4.0, outperforms both interaction energy scoring by DOCK and the Böhm scoring function used by FlexX in rigid and flexible docking of thrombin inhibitors. Our results indicate that for reliable docking of flexible ligands the selection of anchor fragments, conformational sampling and currently available scoring methods still require improvement.     

Comparative study of several algorithms for flexible ligand docking

We have performed a comparative assessment of several programs for flexible molecular docking: DOCK 4.0, FlexX 1.8, AutoDock 3.0, GOLD 1.2 and ICM 2.8. This was accomplished using two different studies: docking experiments on a data set of 37 protein-ligand complexes and screening a library containing 10,037 entries against 11 different proteins. The docking accuracy of the methods was judged based on the corresponding rank-one solutions. We have found that the fraction of molecules docked with acceptable accuracy is 0.47, 0.31, 0.35, 0.52 and 0.93 for, respectively, AutoDock, DOCK, FlexX, GOLD and ICM. Thus ICM provided the highest accuracy in ligand docking against these receptors. The results from the other programs are found to be less accurate and of approximately the same quality. A speed comparison demonstrated that FlexX was the fastest and AutoDock was the slowest among the tested docking programs. The database screening was performed using DOCK, FlexX and ICM. ICM was able to identify the original ligands within the top 1% of the total library in 17 cases. The corresponding number for DOCK and FlexX was 7 and 8, respectively. We have estimated that in virtual database screening, 50% of the potentially active compounds will be found among approximately 1.5% of the top scoring solutions found with ICM and among approximately 9% of the top scoring solutions produced by DOCK and FlexX.     

A shape- and chemistry-based docking method and its use in the design of HIV-1 protease inhibitors

Summary The program DOCK [1,2] has been used successfully to identify molecules which will bind to a specified receptor [3]. The original method ranks molecules based on their shape complementarity to the receptor site and relies on the chemist to bring the appropriate electrostatic or hydrogen bond properties into the molecular skeletons obtained in the search. This is useful when screening a small database of compounds, where it is not likely that molecules with both the correct shape and electrostatic properties will be found. As large databases are more likely to have redundant molecular shapes with a variety of functionality (e.g., members of a congeneric series), it would be useful to have a method which identifies molecules with both the correct shape and functionality. To this end we have modified the DOCK 1.0 method to target user-specified atom types to selected positions in the receptor site. The target sites can be chosen based on structural evidence, calculation or inspection. Targeted-DOCK improves the ability of the DOCK method to find the crystallographically determined binding mode of a ligand. Additionally, targeted-DOCK searches a database of small molecules at 100–1000 times the rate of DOCK 1.0, allowing more molecules to be screened and more sophisticated scoring schemes to be employed. Targeted-DOCK has been used successfully in the design of a novel non-peptide inhibitor of HIV-1 protease.     

The MPSim-Dock hierarchical docking algorithm: application to the eight trypsin inhibitor cocrystals

To help improve the accuracy of protein-ligand docking as a useful tool for drug discovery, we developed MPSim-Dock, which ensures a comprehensive sampling of diverse families of ligand conformations in the binding region followed by an enrichment of the good energy scoring families so that the energy scores of the sampled conformations can be reliably used to select the best conformation of the ligand. This combines elements of DOCK4.0 with molecular dynamics (MD) methods available in the software, MPSim. We test here the efficacy of MPSim-Dock to predict the 64 protein-ligand combinations formed by starting with eight trypsin cocrystals, and crossdocking the other seven ligands to each protein conformation. We consider this as a model for how well the method would work for one given target protein structure. Using as a criterion that the structures within 2 kcal/mol of the top scoring include a conformation within a coordinate root mean square (CRMS) of 1 A of the crystal structure, we find that 100% of the 64 cases are predicted correctly. This indicates that MPSim-Dock can be used reliably to identify strongly binding ligands, making it useful for virtual ligand screening.     

SAMPL4 & DOCK3.7: lessons for automated docking procedures

The SAMPL4 challenges were used to test current automated methods for solvation energy, virtual screening, pose and affinity prediction of the molecular docking pipeline DOCK 3.7. Additionally, first-order models of binding affinity were proposed as milestones for any method predicting binding affinity. Several important discoveries about the molecular docking software were made during the challenge: (1) Solvation energies of ligands were five-fold worse than any other method used in SAMPL4, including methods that were similarly fast, (2) HIV Integrase is a challenging target, but automated docking on the correct allosteric site performed well in terms of virtual screening and pose prediction (compared to other methods) but affinity prediction, as expected, was very poor, (3) Molecular docking grid sizes can be very important, serious errors were discovered with default settings that have been adjusted for all future work. Overall, lessons from SAMPL4 suggest many changes to molecular docking tools, not just DOCK 3.7, that could improve the state of the art. Future difficulties and projects will be discussed.     

Virtual high throughput screening in new lead identification

Drug discovery continues to be one of the greatest contemporary challenges and rational application of modelling approaches is the first important step to obtain lead compounds, which can be optimised further. Virtual high throughput screening (VHTS) is one of the efficient approaches to obtain lead structures for a given target. Strategic application of different screening filters like pharmacophore mapping, shape-based, ligand-based, molecular similarity etc., in combination with other drug design protocols provide invaluable insights in lead identification and optimization. Screening of large databases using these computational methods provides potential lead compounds, thus triggering a meaningful interplay between computations and experiments. In this review, we present a critical account on the relevance of molecular modelling approaches in general, lead optimization and virtual screening methods in particular for new lead identification. The importance of developing reliable scoring functions for non-bonded interactions has been highlighted, as it is an extremely important measure for the reliability of scoring function. The lead optimization and new lead design has also been illustrated with examples. The importance of employing a combination of general and target specific screening protocols has also been highlighted.     

A virtual high throughput screen for high affinity cytochrome P450cam substrates. Implications for in silico prediction of drug metabolism

Structure-based virtual screening techniques require reliable scoring functions to discriminate potential substrates effectively. In this study we compared the performance of GOLD, PMF, DOCK and FlexX scoring functions in FlexX flexible docking to cytochrome P450cam binding site. Crystal structures of protein-substrate complexes were most effectively reproduced by the FlexX/PMF method. On the other hand, the FlexX/GOLD approach provided the best correlation between experimental binding constants and predicted scores. Binding modes selected by the FlexX/PMF approach were rescored by GOLD to obtain a reliable measure of binding energetics. The effectiveness of the FlexX/PMF/GOLD method was demonstrated by the correct classification of 32 out of the 33 experimentally studied compounds and also in a virtual HTS test on a library of 10,000 compounds. Although almost all the available functions were developed to be general, our study on cytochrome P450cam substrates suggests that careful selection or even tailoring the scoring function might increase the prediction power of virtual screens significantly. The FlexX/PMF/GOLD methodology was tested on cytochrome P450 3A4 substrates and inhibitors. This preliminary study revealed that the combined function was able to recognise 334 out of the 345 compounds bound to 3A4.     

Simultaneous screening of four epidermal growth factor receptor antagonists from Curcuma longa via cell membrane chromatography online coupled with HPLC–MS

The epidermal growth factor receptors (EGFRs) are significant targets for screening active compounds. In this work, an analytical method was established for rapid screening, separation, and identification of EGFRs antagonists from Curcuma longa. Human embryonic kidney 293 cells with a steadily high expression of EGFRs were used to prepare the cell membrane stationary phase in a cell membrane chromatography model for screening active compounds. Separation and identification of the retention chromatographic peaks was achieved by HPLC–MS. The active sites, docking extents and inhibitory effects of the active compounds were also demonstrated. The screening result found that ar‐turmerone, curcumin, demethoxycurcumin, and bisdemethoxycurcumin from Curcuma longa could be active components in a similar manner to gefitinib. Biological trials showed that all of four compounds can inhibit EGFRs protein secretion and cell growth in a dose‐dependent manner, and downregulate the phosphorylation of EGFRs. This analytical method demonstrated fast and effective characteristics for screening, separation and identification of the active compounds from a complex system and should be useful for drug discovery with natural medicinal herbs.     

SHAFTS: a hybrid approach for 3D molecular similarity calculation. 1. Method and assessment of virtual screening

We developed a novel approach called SHAFTS (SHApe-FeaTure Similarity) for 3D molecular similarity calculation and ligand-based virtual screening. SHAFTS adopts a hybrid similarity metric combined with molecular shape and colored (labeled) chemistry groups annotated by pharmacophore features for 3D similarity calculation and ranking, which is designed to integrate the strength of pharmacophore matching and volumetric overlay approaches. A feature triplet hashing method is used for fast molecular alignment poses enumeration, and the optimal superposition between the target and the query molecules can be prioritized by calculating corresponding "hybrid similarities". SHAFTS is suitable for large-scale virtual screening with single or multiple bioactive compounds as the query "templates" regardless of whether corresponding experimentally determined conformations are available. Two public test sets (DUD and Jain's sets) including active and decoy molecules from a panel of useful drug targets were adopted to evaluate the virtual screening performance. SHAFTS outperformed several other widely used virtual screening methods in terms of enrichment of known active compounds as well as novel chemotypes, thereby indicating its robustness in hit compounds identification and potential of scaffold hopping in virtual screening.     

Evaluation of library ranking efficacy in virtual screening

We present the results of a comprehensive study in which we explored how the docking procedure affects the performance of a virtual screening approach. We used four docking engines and applied 10 scoring functions to the top-ranked docking solutions of seeded databases against six target proteins. The scores of the experimental poses were placed within the total set to assess whether the scoring function required an accurate pose to provide the appropriate rank for the seeded compounds. This method allows a direct comparison of library ranking efficacy. Our results indicate that the LigandFit/Ligscore1 and LigandFit/GOLD docking/scoring combinations, and to a lesser degree FlexX/FlexX, Glide/Ligscore1, DOCK/PMF (Tripos implementation), LigandFit1/Ligscore2 and LigandFit/PMF (Tripos implementation) were able to retrieve the highest number of actives at a 10% fraction of the database when all targets were looked upon collectively. We also show that the scoring functions rank the observed binding modes higher than the inaccurate poses provided that the experimental poses are available. This finding stresses the discriminatory ability of the scoring algorithms, when better poses are available, and suggests that the number of false positives can be lowered with conformers closer to bioactive ones.     

Combinational numeral fingerprint spectra of Glycyrrhiza and analysis of common peak ratio invariableness in HPLC

To overcome the instability of traditional Chinese medicine (TCM) fingerprint spectra (FPS) and to build up absolute identification standard of TCM, the construction method of combinational numeral fingerprint spectra (CNFPS) was set up. The analysis of invariableness based on CNFPS was carried out. It can be used as absolute quantitative standard to identify the difference of TCM samples. According to this method, the HPLC FPS of components extracted from Glycyrrhiza root samples by ethanol was analyzed. Stable results and accurate quality evaluation were obtained perfectly. The common peak ratio invariableness is an absolute standard in identifying TCM samples, which is better than the similarity method of the FPS.     

CCK1受体的同源模拟和分子对接研究

采用同源建模法对CCK1受体的三维结构进行了模拟,并采用分子动力学方法对模型进行修正和优化,再采用与训练集激动剂和拮抗剂分子对接的方法分别得到激动状态和拮抗状态CCK1受体的三维结构模型。得到的模型使用DOCK对接软件对训练集中的分子进行对接,所得结果与其实际活性拟合度较好,说明我们建立的激动和拮抗状态下的CCK1受体的三维结构模型比较合理,可以作为化合物虚拟筛选的模型对新化合物进行虚拟筛选。     

Pharmacophore alignment search tool: influence of scoring systems on text-based similarity searching

The text-based similarity searching method Pharmacophore Alignment Search Tool is grounded on pairwise comparisons of potential pharmacophoric points between a query and screening compounds. The underlying scoring matrix is of critical importance for successful virtual screening and hit retrieval from large compound libraries. Here, we compare three conceptually different computational methods for systematic deduction of scoring matrices: assignment-based, alignment-based, and stochastic optimization. All three methods resulted in optimized pharmacophore scoring matrices with significantly superior retrospective performance in comparison with simplistic scoring schemes. Computer-generated similarity matrices of pharmacophoric features turned out to agree well with a manually constructed matrix. We introduce the concept of position-specific scoring to text-based similarity searching so that knowledge about specific ligand-receptor binding patterns can be included and demonstrate its benefit for hit retrieval. The approach was also used for automated pharmacophore elucidation in agonists of peroxisome proliferator activated receptor gamma, successfully identifying key interactions for receptor activation.