Skip to main content

IVSPlat 1.0: an integrated virtual screening platform with a molecular graphical interface

Abstract

Background

The virtual screening (VS) of lead compounds using molecular docking and pharmacophore detection is now an important tool in drug discovery. VS tasks typically require a combination of several software tools and a molecular graphics system. Thus, the integration of all the requisite tools in a single operating environment could reduce the complexity of running VS experiments. However, only a few freely available integrated software platforms have been developed.

Results

A free open-source platform, IVSPlat 1.0, was developed in this study for the management and automation of VS tasks. We integrated several VS-related programs into a molecular graphics system to provide a comprehensive platform for the solution of VS tasks based on molecular docking, pharmacophore detection, and a combination of both methods. This tool can be used to visualize intermediate and final results of the VS execution, while also providing a clustering tool for the analysis of VS results. A case study was conducted to demonstrate the applicability of this platform.

Conclusions

IVSPlat 1.0 provides a plug-in-based solution for the management, automation, and visualization of VS tasks. IVSPlat 1.0 is an open framework that allows the integration of extra software to extend its functionality and modified versions can be freely distributed. The open source code and documentation are available at http://kyc.nenu.edu.cn/IVSPlat/.

Background

The successful application of virtual screening in drug discovery means that medicinal chemists and pharmacologists are increasingly using this tool in drug discovery research [1]. Two general strategies are employed in virtual screening: (1) structure-based virtual screening (SBVS) methods for screening compound libraries where the three-dimensional (3D) structures of targets are available; and (2) ligand-based virtual screening (LBVS) methods for identifying potential hits in compound libraries, typically where 3D target structures are unknown [2]. SBVS is dependent on knowledge of the 3D structure of the target. The docking of a compound collection is tested with the target structure and a quantified interaction score is used to identify candidate lead compounds. Thus, SBVS is not dependent on the existence of known active compounds, which increases the prospects for identifying new active lead compounds [3]. Several docking programs have been developed for virtual screening since the initial development of UCSF Dock [4, 5], such as AutoDock [6], GOLD [7], and GLIDE [8, 9]. There have been reports of the successful identification of lead compounds using docking-based VS methods [10]. LBVS assumes that structurally similar compounds are likely to exhibit similar biological activities. LBVS makes use of the physicochemical, structural, and energetic properties of known active compounds when searching large compound libraries for related or novel chemical compounds [11]. The pharmacophore-based database searching technique is a widely used VS strategy [12], which relies on knowledge of the biological activity of multiple hits when identifying key features during a search. A pharmacophore is a spatial arrangement of features that allows a compound to interact with a target receptor at a specific binding mode. Thus, a pharmacophore can serve as a powerful model with versatile applications in rational drug design, such as virtual screening, de novo design, lead optimization, and ADME/Tox studies. Several algorithms and software tools have been developed for pharmacophore detection [13–19]. In practice, SBVS and LBVS can be used in combination or separately, depending on the specific aim. However, very few open-source platforms have integrated both approaches for virtual screening [20].

In this study, we developed an integrated virtual screening platform, IVSPlat 1.0, which is an easily operated PyMOL [21] plug-in that performs many practical virtual screening tasks derived from UCSF Dock 6.5 [22] and PharmaGist [18], while also providing the facility to analyze results. Molecular visualization can be used to highlight connectivity and important structural features in interacting molecules, which is crucial for drug design. IVSPlat 1.0 was developed to exploit PyMOL's exceptional molecular viewing capabilities, which are seamlessly integrated with other programs. Thus, IVSPlat 1.0 facilitates molecular docking and pharmacophore studies, which can be visualized in a PyMOL window. Furthermore, IVSPlat 1.0 is an open source framework, thereby allowing developers to add other new open source computational applications using a uniform graphical user interface (GUI). IVSPlat 1.0 was written in Python and its installation instructions and video demonstrations can be found on the IVSPlat website.

Implementation

Related technologies

IVSPlat 1.0 integrates the programs detailed in the following sections and it also provides a uniform interface for accessing these programs, which facilitates transparency for the user. Auxiliary software can also be incorporate provided that they are preinstalled, such as Openbabel [23], DMS [24], Sphgen [6], and MPICH2 [25].

PyMol

PyMol is one of the most popular molecular visualization programs. The core of the program is a full-featured Python [26] interpreter, which is extended by an OpenGL 3D display, a Tcl/Tk based GUI, and a PyMol API that facilitates plug-in creation. Several popular plug-in extensions have been developed in the field of molecular interactions. The APBS plug-in [27] provides an interface for the popular adaptive Poisson-Boltzmann solver (APBS) program, thereby giving easy access to electrostatics calculations, and visualization of potential energy surfaces and charge densities on protein surfaces. PyDeT [28] is a PyMOL plug-in for visualizing tessellations derived from the protein structure and the source protein. PyETV [29] is a PyMOL plug-in used for viewing, analyzing, and manipulating the prediction of evolutionarily important residues and sites in protein structures and their complexes. The Autodock/Vina plug-in [30] provides an interface between PyMOL and two popular docking programs, Autodock and Autodock Vina [31], thereby facilitating molecular docking, virtual screening, and binding site analysis via PyMOL.

UCSF Dock 6.5

UCSF Dock is one of the most frequently used molecular docking applications in VS studies [2]. The latest version of UCSF Dock is version 6.5 (UCSF Dock 6.5). The following new features have been added to UCSF Dock 6.5: 1) an anchor can be selected by specifying an atom in a fragment and the number of anchors used can be limited during multi-anchor docking; 2) the new scoring function (descriptor score) has been introduced, which includes a hydrogen bond term and footprint similarity scoring; and 3) PB/SA score has been generalized to some extent and efficiency improvements have been included that make docking, rather than rescoring, more tractable for nontrivial systems [22]. However, no GUI tool is currently available for automating VS tasks using UCSF Docks.

PharmaGist

IVSPlat 1.0 uses the freely available program PharmaGist [18] for LBVS. PharmaGist is a virtual screening method that is based on pharmacophore detection. Its workflow is as follows: 1) specify a set of active compounds; 2) generate candidate pharmacophores that are solved using an algorithm via multiple flexible alignment of the active compounds; and 3) perform pharmacophore-based virtual screening. PharmaGist is highly efficient and a typical run of up to 32 active compounds takes from seconds to a few minutes using a standard PC [18].

PaDEL-Descriptor

PaDEL [32] is an open source program that can calculate 801 descriptors (667 1D/2D descriptors and 134 3D descriptors) and 10 types of fingerprints. The descriptors and fingerprints are calculated using the Chemistry Development Kit [33] with additional descriptors and fingerprints. These additions include atom type electrotopological state descriptors, McGowan volume, molecular linear free energy relation descriptors, ring counts, counts of chemical substructures identified by Laggner, binary fingerprints, and counts of chemical substructures identified by Klekota and Roth. PubChem fingerprints options were used to produce molecular descriptors in our platform.

Hcluster

Hcluster [34] is a hierarchical clustering library developed in the Python language with the Scipy package [35]. This library provides Python functions for agglomerative clustering, including the generation of hierarchical clusters from distance matrices, computation of distance matrices from observation vectors, calculation of cluster statistics, linkage cutting to generate flat clusters, and cluster visualization with dendrograms. Hcluster also provides algorithms for clustering, including, linkage, single, complete, average, median, and Ward's method. Ward's method is the most adaptive for molecular clustering [36], so it was used as the default in our platform. However, users can still select their own preferred method. Hcluster outputs a dendrogram to illustrate the cluster results.

Architecture of IVSPlat 1.0

Figure 1 shows that the architecture of the platform consists of three layers: application layer, class interface layer, and software entity layer. The class interface layer includes a class library written in Python that constitutes the core of the platform. Several software programs are integrated with the molecular graphics system in the software entity layer, while various virtual screening services are provided in the application layer. Python is an interactive, object-oriented, extensible programming language with strong support for integration with other languages, so it was selected as the main programming tool to produce a uniform software environment for the simultaneous application of various scientific programs. The class library in the class interface layer has four components: GUI class, software interface class, results analysis class, and format reading class. The software interface class is the core, which consists of the Thread_run and IVSPlat classes. IVSPlat class is the main class for calling other classes in all the methods used in the application layer, including the thread distribution, conversion of file formats between the integrated programs, database operation, and the analysis of the VS results. The application layer is composed of a GUI and three functional modules. The GUI provides functions for visualizing the binding geometries and interactions between ligands (compounds) and receptors using PyMOL, while the three functional modules allow docking-based VS, pharmacophore-based VS, and VS results analysis. The platform is launched by running PyMOL and it allows the user to specify the runtime parameters. The main platform screen is shown in Figure 2(a), which includes nine functional items: configuration, grid, receptor, ligand, docking, analysis, pharmacophore, cluster, and database management.

Figure 1
figure 1

Architecture of the IVSPlat 1.0 VS platform.

Figure 2
figure 2

Main applications of the IVSPlat 1.0 VS platform. (a) Main screen of the IVSPlat 1.0 platform. (b) Preparation of receptor binding site(s). (c) Preparation of ligand database. (d) Definition of the grid box. (e) Docking-based VS panel with the docked results shown in a table. (f) Pharmacophore-based VS panel with the results listed in a table. (g) Clustering panel and clustering results displayed as a dendrogram. (h) Database management.

Molecular docking-based virtual screening

A molecular docking-based virtual screening task can be conducted starting with a target protein and a large, diverse compound library. The main steps are described in the following sections.

Running environment configuration

This module defines the location and the environment variables required for the execution of each program. To facilitate system maintenance and reconfiguration, all the correctly configured options can be saved to a text file, which can be manually modified as required.

Receptor and ligand preparation

Receptor preparation in UCSF Dock 6.5 includes a series of tasks, including, adding hydrogen atoms and charges, protonating amines, deleting solvent molecules [37], and localizing the receptor binding sites. These tasks are frequently accomplished by manual operation and inspection. In particular, the localization of binding site(s) is composed of three related operations: 1) generating the molecular surface of the receptor using the DMS program; 2) generating spheres surrounding the receptor using the Sphgen program; and 3) selecting a subset of spheres to represent the binding site(s) using the Sphere_selector program. In IVSPlat 1.0, the receptor functional module provides a workbench with an interactive GUI that automates the process of binding site preparation. Users are not required to consider the inputs and outputs of all the intermediate processes. The resulting binding site(s) can also be visualized and analyzed in a PyMOL window (Figure 2(b)).

During the preparation of ligands (compounds), UCSF Dock 6.5 and PharmaGist have to add charges and hydrogen atoms to ligands in a MOL2 format. It is recommended that all ligand structures used in VS experiments are obtained from the ZINC database [38]. ZINC is the most popular free database for virtual screening and it contains compounds that are commercially available. Moreover, the charges and hydrogen atoms have already been added to each compound in the ZINC compound database. Therefore, our platform does not require a functionality for adding charges and hydrogen atoms. All compounds can be inspected and selected via the PyMOL window in the ligand functional module, before a compound database for virtual screening is set up (Figure 2(c)). This module also converts other formats to MOL2 during ligand preparation by calling the Openbabel program.

Grid generation

UCSF Dock 6.5 requires a rectangular box to generate a grid file and the grid functional module provides a graphical environment for defining the centre of the box, adjusting the box size and position, and generating the final box file. PyMOL defined the box centre by selecting explicit coordinates or by specifying a PyMOL selection (e.g., a reference ligand) to calculate the mean coordinates of the atoms. The docking box can be displayed in the PyMOL window while the size and the exact position of the box can be adjusted depending on the user's requirements (Figure 2(d)).

Docking

Depending on a user's computer configuration, UCSF Dock 6.5 can be installed in several different running modes, including, GNU, GNU parallel, and GNU parallel pbsa [22]. UCSF Dock 6.5 also has a wealth of parameters for tailoring a particular VS service, including a choice of multiple scoring functions. UCSF Dock 6.5 typically runs from the command line in a standard Unix shell and it reads an input parameter file containing all the field/value options. We created a standard input file known as dock.in and a control panel is provided so users can conveniently adjust the running parameters in a file (Figure 2(e)). During docking experiments, the grid box size must be set high enough to allow all poses of the ligands to lie within the box, although excessively large grid boxes will waste computational resources. Thus, we placed an option in the control panel to facilitate trial and error, i.e., max_orientations. The primary score parameter option in the control panel helps users to select an appropriate scoring function. Furthermore, the parallel version of UCSF Dock 6.5 with a PB/SA scoring function can be employed if users install the OpenEye Toolkits [39], which are freely available to most academic and government institutions.

Analysis of results

A series of output files are created after the docking task is complete, including an anchor_and_grow_ranked.mol2 file containing the top scored ligands sorted using a score function value. Users can list all the data in the anchor_and_grow_ranked.mol2 file as a table in the analysis functional module, followed by loading the data in a Pymol window to visualize and inspect ligand poses. In general, the number of hits in a virtual screening task far exceeds the number of compounds that can be biologically tested, so a molecular clustering procedure is required to reduce the number of hits to a subset of compounds that best represents the overall result set. In IVSPlat 1.0, the Hcluster Python library and the PaDEL-Descriptor program are used to produce a molecular clustering module. The PaDEL-Descriptor program is used to calculate the molecular descriptors and fingerprints, while the Hcluster library is used to perform clustering tasks based on the descriptors and fingerprints. The clustering results are displayed in a hierarchical dendrogram (Figure 2(g)), which can be saved as a text file.

Pharmacophore-based virtual screening

IVSPlat 1.0 provides a pharmacophore-based virtual screening module via PharmaGist. The module is simple and easy to use, and users are only required to upload a ligand file and a compound database file as input, before setting an optional running parameter, such as the scoring weight of each feature. The module has two main components: generating pharmacophores and screening the compound database with the resulting pharmacophore. The method is similar to that used in molecular docking-based virtual screening. The pharmacophores can be loaded and analyzed in a PyMOL window before selecting the appropriate public pharmacophore (Figure 2(f)). An output such as *mol2_dbres.html (* means any string) is created after the VS task is complete, which contains each output compound's score and structure information. The screening results can be summarized as a table in score order and these data can be visualized in a PyMOL window (Figure 2(f)).

Combined virtual screening

Different VS approaches are characterized by their screening strategies, i.e., SBVS or LBVS, which can be used independently or in combination. Combined virtual screening is regarded as a novel VS protocol that efficiently increases hit rates during virtual screening, while reducing the computational costs [40–42]. IVSPlat 1.0 provides an environment for combined virtual screening using pharmacophore-based VS and molecular docking-based VS. Pharmacophore-based VS can be used as pre-filter during the virtual screening of a huge database, before performing molecular docking-based VS only with the top results obtained during the first step [40]. Pharmacophore-based VS and molecular docking-based VS can also be used in parallel to improve the enrichment factor (EF). An example is presented in the next section to demonstrate the applicability of this platform.

Database management

All the input and output files used in the current version of the IVSPlat platform are text files, including the compound database files (*.mol2). A suite of text access methods were developed to overcome compatibility issues in input/output files of the integrated software tools. A Python interface to SQLite (a lightweight SQL database) was developed, which is mainly used by IVSPlat 1.0 for the data management of VS output results (Figure 2(h)). This provides a base for further extension.

Application of the platform - a case study

Experimental dataset

The experimental dataset was taken from the DUD database [43, 44], which was designed for testing docking algorithms by providing challenging decoys. DUD contains 2950 active compounds with 40 targets. We used thrombin (PDB id: 1H8D, H chain) as the target protein in our experiments. Thrombin is a trypsin-like serine protease with a central role in thrombosis and haemostasis and its inhibitors are potentially pharmacologically important antithrombotics [45]. The X-ray crystal structures of thrombin are known and structures of many small molecule inhibitors complexed with thrombin have been identified [43–45]. We also used 72 compounds that are active against thrombin and we randomly selected 1000 compounds as decoys. We tested the applicability of the platform by conducting virtual screening of thrombin inhibitors using several screening strategies.

Screening for thrombin inhibitors using molecular docking-based VS

The target protein (thrombin, PDB id: 1H8D) was prepared using UCSF Chimera [38] by the following two steps: 1) remove L chain, I chain (inhibitor), and water molecules from the crystal structure of thrombin; and 2) add all hydrogen atoms and all the Gasteiger atom charges. The residues within 10.0 Ã… of the bound ligand were defined as the active site. To reduce the collision probability, the parameter min_anchor_size was set to 5. The pruning_conformer_score_cutoff was set at 100, the max_orientations was set at 1000, the maximum number of orientations to be scored was set at 1000, and the grid box size was set at 10. We treated all compounds as flexible using DOCK's anchor-and-grow algorithm while the grid energy scoring method was used to evaluate compounds. UCSF Dock 6.5 in parallel mode was used for the screening task (using a Dell precision T7400 workstation with two processors: Intel Xeon E5430 quad core). The total running time was c. 2 h. The results were analyzed in terms of EF. EF is a common metric used for comparing virtual screening results, which is defined as follows:

E F  =  H I T S s a m p l e d N s a m p l e d H I T S t o t a l N t o t a l
(1)

where Ntotal is the number of compounds in the database, HITStotal is the total number of known active compounds, Nsampled is the number of top compounds to be sampled, and HITSsampled is the number of known active compounds found in the Nsampled compounds.

Molecular docking-based VS provided good performance in the experiments. Figure 3 shows that the maximum EF was 14 with the top 1% of the database (Figure 3(a)), while the frequency of thrombin-active compounds in the top 1% of the total library was c. 15% (Figure 3(b)). The frequency of active compounds that were screened gradually increased as we increased the percentage of the database that was sampled. In contrast, the EFs gradually decreased and c. 90% of the active compounds was successfully screened from the top 50% of the database (Figure 3(b)).

Figure 3
figure 3

EF of docking-based VS and the corresponding percentage of active compounds. EF of docking-based VS. (b) The corresponding percentage of active compounds.

Screening thrombin inhibitors using pharmacophore-based VS

The current version of PharmaGist can only handle up to 32 compounds when generating pharmacophores, so we separated the 72 compounds into eight equal groups with nine compounds in each. We used the nine compounds in each group as seeds for constructing eight test datasets, by taking extra random samples from the other groups. Thus, eight test datasets were constructed and each dataset was again divided randomly into four small groups that contained 32, 16, 8, and 4 compounds respectively. The running parameters of PharmaGist were set as follows: the number of output pharmacophores was set at 2; the number of output results was set at 1000; while other parameters were set at the default values. The results for the first and fourth dataset are shown in Figure 4. Figure 4(a-1) shows that the 32-compound group used for generating pharmacophores performed better than the other cases. Its maximum enrichment factor was 12.5 with the top 1% from the database and it retrieved c. 90% of active compounds from the top 50% in the database (Figure 4(a-2)). The enrichment factor and the percentage of active compounds found in the experiments declined steadily as the number of compounds used for generating pharmacophores was decreased. The first dataset was the worst case in our experiments and Figure 4(b-1) and Figure 4(b-2) show that the maximum enrichment factor was 4.2 with the top 1% from the database (Figure 4(b-1)) while 47% of the active compounds were retrieved from the top 50% in the database (Figure 4(b-2)). The results with other datasets are shown in Supplementary Figure S1 (see Additional file 1). These data indicated that the performance of the pharmacophore-based VS generally tended to be influenced by the compounds used to generate the pharmacophores. Thus, the quality and number of the compounds selected were key factors in the pharmacophore-based VS.

Figure 4
figure 4

EF and corresponding percentages of active compounds with pharmacophore-based VS using the fourth and first datasets. (a-1) EF of pharmacophore-based VS with the fourth experimental dataset. (a-2) The corresponding percentages of active compounds with the fourth experimental dataset. (b-1) Enrichment factor when screening the first experimental dataset. (b-2) The corresponding percentages of active compounds when screening the first experimental dataset. The plot shows the results with four ligand groups: 32 ligands (blue); 16 ligands (red); 8 ligands (green); and 4 ligands (purple).

Combination of the two VS protocols

Pharmacophore-based VS is suitable when a set of known active compounds is available. The strength of the pharmacophore-based VS is its high speed. However, docking-based VS is more appropriate when the structure of a protein target is known. Compared with pharmacophore-based VS, docking-based VS is more suitable for searching for active compounds with novel scaffolds, although its running speed is generally slow. A combination of both VS protocols may provide a superior alternative to the single VS protocol when the structure of the target protein and a set of known active compounds are available. Table 1 shows the results obtained with a combination of both protocols with the fourth dataset. Figure 5 shows that 65 active compounds were successfully screened using pharmacophore- and docking- based VS. We found that 58 compounds were screened with combination of both, whereas 14 of the 72 active compounds were screened when using only one of the two protocols. Screening results for the remaining seven datasets are shown in Supplementary Tables S1-S7 and Figures S2-S8 (see Additional file 1). In summary, pharmacophore- and docking- based VS identified some common active compounds, but they also identified different individual active compounds, indicating the complementary effects of each method. Thus, a combination of both VS protocols is feasible.

Table 1 Screening results with the fourth dataset for the top 50% of the database using dock- and pharmacophore- based VS, and a combination of both
Figure 5
figure 5

Number of ligands identified by docking- and pharmacophore- based VS, and a combination of both, with the fourth dataset for the top 50% of the database. "Fail of both" indicates the number of ligands we failed to identify using both methods.

Pharmacophore-based methods can be used as a pre-filter during the screening of large databases to reduce data input with docking-based VS. Alternatively, the two methods can be used in parallel to improve the hit rate with VS.

Conclusions

IVSPlat 1.0 is a plug-in based software platform for the management of VS tasks and the analysis of VS results, which integrates pharmacophore-based VS, docking-based VS, a molecular graphics system, and the analysis of results in a single, user-friendly workbench. IVSPlat 1.0 is an open framework which provides the potential for further development. Work is ongoing to improve the platform and allows researchers to perform QSAR model-based VS, ADMET-based virtual screening, and knowledge-based virtual screening. Researchers can also use more auxiliary programs for workbench construction and results analysis. Another major feature of this project will be the development of online updates. The platform will be updated on the IVSPlat web page.

Availability and requirements

Project name: IVSPlat 1.0

Project home page: http://kyc.nenu.edu.cn/IVSPlat/

Operating system(s): Linux (Fedora 14)

Programming language: Python

Other requirements: PyMOL (v 1.4.1), Dock 6.5, Pharmagist, PaDEL-Descriptor, Hcluster (v 0.2.0), NumPy (v 1.6), DMS, MPICH2 and Openbabel (v 2.3)

License: GNU GPL

Any restrictions to use by non-academics: Contact authors

References

  1. Shoichet BK: Virtual screening of chemical libraries. Nature. 2004, 432: 862-865. 10.1038/nature03197.

    Article  CAS  Google Scholar 

  2. Lengauer T, Lemmen C, Rarey M, Zimmermann M: Novel technologies for virtual screening. Drug Discov Today. 2004, 9: 27-34. 10.1016/S1359-6446(04)02939-3.

    Article  CAS  Google Scholar 

  3. Kitchen DB, Decornez H, Furr JR, Bajorath J: Docking and scoring in virtual screening for drug discovery: methods and applications. Nat Rev Drug Discov. 2004, 3: 935-949. 10.1038/nrd1549.

    Article  CAS  Google Scholar 

  4. Tuccinardi T: Docking-based virtual screening: recent developments. Comb Chem High Throughput Screen. 2009, 12: 303-314. 10.2174/138620709787581666.

    Article  CAS  Google Scholar 

  5. Kuntz ID, Blaney JM, Oatley SJ, Langridge R, Ferrin TE: A geometric approach to macromolecule-ligand interactions. J Mol Biol. 1982, 161: 269-288. 10.1016/0022-2836(82)90153-X.

    Article  CAS  Google Scholar 

  6. Morris GM, Goodsell DS, Halliday DS, Huey R, Hart WE, Belew R, Olson AJ: Automated docking using a Lamarckian genetic algorithm and and empirical binding free energy function. J Comp Chem. 1998, 19: 1639-1662. 10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-B.

    Article  CAS  Google Scholar 

  7. Jones G, Willett P, Glen RC, Leach AR, Taylor R: Development and validation of a genetic algorithm for flexible docking. J Mol Bio. 1997, 267: 727-748. 10.1006/jmbi.1996.0897.

    Article  CAS  Google Scholar 

  8. Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, Repasky MP, Knoll EH, Shelley M, Perry JK, Shaw DE, Francis P, Shenkin PS: Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J Med Chem. 2004, 47: 1739-1749. 10.1021/jm0306430.

    Article  CAS  Google Scholar 

  9. Halgren TA, Murphy RB, Friesner RA, Beard HS, Frye LL, Pollard WT, Banks JL: Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J Med Chem. 2004, 47: 1750-1759. 10.1021/jm030644s.

    Article  CAS  Google Scholar 

  10. Sousa SF, Cerqueira NM, Fernandes PA, Ramos MJ: Virtual screening in drug design and development. Comb Chem High Throughput Screen. 2010, 13: 442-453. 10.2174/138620710791293001.

    Article  CAS  Google Scholar 

  11. Stahura FL, Bajorath J: New Methodologies for Ligand-Based Virtual Screening. Curr Pharm Des. 2005, 11: 1189-1202. 10.2174/1381612053507549.

    Article  CAS  Google Scholar 

  12. Khedkar SA, Malde AK, Coutinho EC, Srivastava S: Pharmacophore modeling in drug discovery and development: an overview. Med Chem. 2007, 3: 187-197. 10.2174/157340607780059521.

    Article  CAS  Google Scholar 

  13. Chen X, Rusinko A, Tropsha A, Young SS: Automated pharmacophore identification for large chemical data sets. J Chem Inf Comput Sci. 1999, 39: 887-896. 10.1021/ci990327n.

    Article  CAS  Google Scholar 

  14. Feng J, Sanil A, Young SS: PharmID: pharmacophore identification using Gibbs sampling. J Chem Inf Model. 2006, 46: 1352-1359. 10.1021/ci050427v.

    Article  CAS  Google Scholar 

  15. Richmond NJ, Abrams CA, Wolohan PR, Abrahamian E, Willett P, Clark RD: GALAHAD: 1. pharmacophore identification by hypermolecular alignment of ligands in 3D. J Comput Aided Mol Des. 2006, 20: 567-587. 10.1007/s10822-006-9082-y.

    Article  CAS  Google Scholar 

  16. Zhu F, Agrafiotis DK: Recursive distance partitioning algorithm for common pharmacophore identification. J Chem Inf Model. 2007, 47: 1619-1625. 10.1021/ci7000583.

    Article  CAS  Google Scholar 

  17. Schneidman-Duhovny D, Dror O, Inbar Y, Nussinov R, Wolfson HJ: Deterministic pharmacophore detection via multiple flexible alignment of drug-like molecules. J Comput Biol. 2008, 15: 737-754. 10.1089/cmb.2007.0130.

    Article  CAS  Google Scholar 

  18. Schneidman-Duhovny D, Dror O, Inbar Y, Nussinov R, Wolfson HJ: PharmaGist: a webserver for ligand-based pharmacophore detection. Nucleic Acids Res. 2008, 36: W223-228. 10.1093/nar/gkn187.

    Article  CAS  Google Scholar 

  19. Podolyan Y, Karypis G: Common pharmacophore identification using frequent clique detection algorithm. J Chem Inf Model. 2009, 49: 13-21. 10.1021/ci8002478.

    Article  CAS  Google Scholar 

  20. Cabrera AC, Gil-Redondo R, Perona A, Gago F, Morreale A: VSDMIP 1.5: an automated structure- and ligand-based virtual screening platform with a PyMOL graphical user interface. J Comput Aided Mol Des. 2011, 25: 813-824. 10.1007/s10822-011-9465-6.

    Article  CAS  Google Scholar 

  21. The PyMOL molecular graphics system. Available online: [http://www.pymol.org] (accessed on 12 October 2011)

  22. DOCK 6.5 Users Manual. Available online: [http://dock.compbio.ucsf.edu/DOCK_6/dock6_manual.htm] (accessed on 12 October 2011)

  23. The Openbabel programs. Available online: [http://openbabel.org] (accessed on 12 October 2011)

  24. Baker N, Sept D, Joseph S, Holst M, McCammon J: Electrostatics of nanosystems: application to microtubules and the ribosome. Proc Natl Acad Sci USA. 2001, 98: 10037-10041. 10.1073/pnas.181342398.

    Article  CAS  Google Scholar 

  25. The Mpich2 programs. Available online: [http://www.mcs.anl.gov/research/projects/mpich2/] (accessed on 12 October 2011)

  26. Sanner M: Python: a programming language for software integration and development. J Mol Graphics Mod. 1999, 17: 57-61.

    CAS  Google Scholar 

  27. Apbs plugin for pymol. Available online: [http://sourceforge.net/projects/pymolapbsplugin/] (accessed on 12 October 2011)

  28. Ordog R: PyDeT, a PyMOL plug-in for visualizing geometric concepts around proteins. Bioinformation. 2008, 2: 346-347.

    Article  Google Scholar 

  29. Lua RC, Lichtarge O: PyETV: a PyMOL evolutionary trace viewer to analyze functional site predictions in protein complexes. Bioinformatics. 2010, 26: 2981-2982. 10.1093/bioinformatics/btq566.

    Article  CAS  Google Scholar 

  30. Seeliger D, de Groot BL: Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J Comput Aided Mol Des. 2010, 24: 417-422. 10.1007/s10822-010-9352-6.

    Article  CAS  Google Scholar 

  31. Trott O, Olson A: AutoDock/Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010, 31: 455-461.

    CAS  Google Scholar 

  32. Yap CW: PaDEL-descriptor: an open source software to calculate molecular descriptors and fingerprints. J Comput Chem. 2011, 32: 1466-1474. 10.1002/jcc.21707.

    Article  CAS  Google Scholar 

  33. The Chemistry Development Kit (CDK). Available online: [http://sourceforge.net/apps/mediawiki/cdk/index.php?title=Main_Page] (accessed on 12 October 2011)

  34. Hcluster package of Scipy. Available online: [http://code.google.com/p/scipy-cluster/] (accessed on 12 October 2011)

  35. The Scipy package. Available online: [http://www.scipy.org] (accessed on 12 October 2011)

  36. Bocker A, Schneider G, Teckentrup A: NIPALSTREE: a new hierarchical clustering approach for large compound libraries and its application to virtual screening. J Chem Inf Model. 2006, 46: 2220-2229. 10.1021/ci050541d.

    Article  Google Scholar 

  37. Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem. 2004, 25: 1605-1612. 10.1002/jcc.20084.

    Article  CAS  Google Scholar 

  38. Irwin JJ, Shoichet BK: ZINC - a free database of commercially available compounds for virtual screening. J Chem Inf Model. 2005, 45: 177-182. 10.1021/ci049714+.

    Article  CAS  Google Scholar 

  39. OpenEye Toolkits. Available online: [http://www.eyesopen.com/] (accessed on 12 October 2011)

  40. Liu X, Xie H, Luo C, Tong L, Wang Y, Peng T, Ding J, Jiang H, Li H: Discovery and SAR of thiazolidine-2,4-dione analogues as insulin-like growth factor-1 receptor (IGF-1R) inhibitors via hierarchical virtual screening. J Med Chem. 2010, 53: 2661-2665. 10.1021/jm901798e.

    Article  CAS  Google Scholar 

  41. Ren JX, Li LL, Zheng RL, Xie HZ, Cao ZX, Feng S, Pan YL, Chen X, Wei YQ, Yang SY: Discovery of novel Pim-1 kinase inhibitors by a hierarchical multistage virtual screening approach based on SVM model, pharmacophore, and molecular docking. J Chem Inf Model. 2011, 51: 1364-1375. 10.1021/ci100464b.

    Article  CAS  Google Scholar 

  42. Vidovic D, Busby SA, Griffin PR, Schürer SC: A combined ligand- and structure-based virtual screening protocol identifies submicromolar PPARγ partial agonists. Chem Med Chem. 2011, 6: 94-103.

    Article  CAS  Google Scholar 

  43. Huang N, Shoichet BK, Irwin JJ: Benchmarking sets for molecular docking. J Med Chem. 2006, 49: 6789-6801. 10.1021/jm0608356.

    Article  CAS  Google Scholar 

  44. The DUD database. Available online: [http://dud.docking.org/] (accessed on 12 October 2011)

  45. Bovio JA, Smith SM, Gums JG: Dabigatran etexilate: a novel oral thrombin inhibitor for thromboembolic disease. Ann Pharmacother. 2011, 45: 603-614. 10.1345/aph.1P644.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Natural Science Foundation of China (No. 61172183), National Key New Drug Creation and Manufacturing Program of Ministry of Science and Technology(No. 2011ZX09401-305-03), the Fundamental Research Funds for the Central Universities, the Natural Science Foundation of Jilin Province (No. 20101503), the Grant of Jilin Province Science & Technology Committee, China(Nos. 20102203, 20110938 and 20110711) and the Scientific and Technical Project of Administration of Traditional Chinese Medicine of Jilin Province (Nos. 2010pt067 and 2011- zd16).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Yan Xin Huang, Jun Kong or Yu Xin Li.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

YXS and YXH conceived and designed the research. YXS, HYW, FLL and CF performed the research including data collection, test and analysis. LGS and ZQM carefully proofread the manuscript. YLB and JK suggested extension and modifications to the research. YXL supervised the whole research and revised the manuscript critically. All authors have read and approved the final manuscript.

Electronic supplementary material

13065_2011_342_MOESM1_ESM.DOC

Additional file 1: Supplementary experimental results. The file contains Supplementary Tables S1 to S7 and Supplementary Figures S1 to S8. (DOC 695 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This is an open access article distributed under the terms of the Creative Commons Attribution Noncommercial License ( https://creativecommons.org/licenses/by-nc/2.0 ), which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.

Reprints and permissions

About this article

Cite this article

Sun, Y.X., Huang, Y.X., Li, F.L. et al. IVSPlat 1.0: an integrated virtual screening platform with a molecular graphical interface. Chemistry Central Journal 6, 2 (2012). https://doi.org/10.1186/1752-153X-6-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1752-153X-6-2

Keywords