Virtual Ligand Screening Using PL-PatchSurfer2, a Molecular Surface-Based Protein-Ligand Docking Method.

Virtual screening is a computational technique for predicting a potent binding compound for a receptor protein from a ligand library. It has been a widely used in the drug discovery field to reduce the efforts of medicinal chemists to find hit compounds by experiments.Here, we introduce our novel structure-based virtual screening program, PL-PatchSurfer, which uses molecular surface representation with the three-dimensional Zernike descriptors, which is an effective mathematical representation for identifying physicochemical complementarities between local surfaces of a target protein and a ligand. The advantage of the surface-patch description is its tolerance on a receptor and compound structure variation. PL-PatchSurfer2 achieves higher accuracy on apo form and computationally modeled receptor structures than conventional structure-based virtual screening programs. Thus, PL-PatchSurfer2 opens up an opportunity for targets that do not have their crystal structures. The program is provided as a stand-alone program at http://kiharalab.org/plps2 . We also provide files for two ligand libraries, ChEMBL and ZINC Drug-like.

PL-PatchSurfer2: Improved Local Surface Matching-Based Virtual Screening Method That Is Tolerant to Target and Ligand Structure Variation.

Virtual screening has become an indispensable procedure in drug discovery. Virtual screening methods can be classified into two categories: ligand-based and structure-based. While the former have advantages, including being quick to compute, in general they are relatively weak at discovering novel active compounds because they use known actives as references. On the other hand, structure-based methods have higher potential to find novel compounds because they directly predict the binding affinity of a ligand in a target binding pocket, albeit with substantially lower speed than ligand-based methods. Here we report a novel structure-based virtual screening method, PL-PatchSurfer2. In PL-PatchSurfer2, protein and ligand surfaces are represented by a set of overlapping local patches, each of which is represented by three-dimensional Zernike descriptors (3DZDs). By means of 3DZDs, the shapes and physicochemical complementarities of local surface regions of a pocket surface and a ligand molecule can be concisely and effectively computed. Compared with the previous version of the program, the performance of PL-PatchSurfer2 is substantially improved by the addition of two more features, atom-based hydrophobicity and hydrogen-bond acceptors and donors. Benchmark studies showed that PL-PatchSurfer2 performed better than or comparable to popular existing methods. Particularly, PL-PatchSurfer2 significantly outperformed existing methods when apo-form or template-based protein models were used for queries. The computational time of PL-PatchSurfer2 is about 20 times shorter than those of conventional structure-based methods. The PL-PatchSurfer2 program is available at http://www.kiharalab.org/plps2/ .

Three-dimensional compound comparison methods and their application in drug discovery.

Virtual screening has been widely used in the drug discovery process. Ligand-based virtual screening (LBVS) methods compare a library of compounds with a known active ligand. Two notable advantages of LBVS methods are that they do not require structural information of a target receptor and that they are faster than structure-based methods. LBVS methods can be classified based on the complexity of ligand structure information utilized: one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D). Unlike 1D and 2D methods, 3D methods can have enhanced performance since they treat the conformational flexibility of compounds. In this paper, a number of 3D methods will be reviewed. In addition, four representative 3D methods were benchmarked to understand their performance in virtual screening. Specifically, we tested overall performance in key aspects including the ability to find dissimilar active compounds, and computational speed.

Measured and predicted affinities of binding and relative potencies to activate the AhR of PAHs and their alkylated analogues.

Polycyclic aromatic hydrocarbons (PAHs) and their alkylated forms are important components of crude oil. Both groups of PAHs have been reported to cause dioxin-like responses, mediated by aryl hydrocarbon receptor (AhR). Thus, characterization of binding affinity to the AhR of unsubstituted or alkylated PAHs is important to understand the toxicological consequences of oil contamination on ecosystems. We investigated the potencies of major PAHs of crude oil, e.g., chrysene, phenanthrene and dibenzothiophene, and their alkylated forms (n=17) to upregulate expression of AhR-mediated processes by use of the H4IIE-luc transactivation bioassay. In addition, molecular descriptors of different AhR activation potencies among PAHs were investigated by use of computational molecular docking models. Based on responses of the H4IIE-luc in vitro assay, it was shown that potencies of PAHs were determined by alkylation in addition to the number and conformation of rings. Potencies of AhR-mediated processes were generally greater when a chrysene group was substituted, especially in 1-methyl-chrysene. Significant negative correlations were observed between the in vitro dioxin-like potency measured in H4IIE-luc cells and the binding distance estimated from the in silico modeling. The difference in relative potency for AhR activation observed among PAHs and their alkylated forms could be explained by differences among binding distances in the ligand binding domain of the AhR caused by alkylation. The docking model developed in the present study may have utility in predicting risks of environmental contaminants of which toxicities are mediated by AhR binding.