Finding pharmaceutically relevant target conformations from an arbitrary set of protein conformations remains a challenge in structure-based virtual screening (SBVS). The growth in the number of available conformations, either experimentally determined or computationally derived, obscures the situation further. While the inflated conformation space potentially contains viable druggable targets, the increase of conformational complexity, as a consequence, poses a selection problem. To address this challenge, we took advantage of machine learning methods, namely an over-sampling and a binary classification procedure, and present a novel method to select druggable receptor conformations. Specifically, we trained a binary classifier on a set of nuclear receptor conformations, wherein each conformation was labeled with an enrichment measure for a corresponding SBVS. The classifier enabled us to formulate suggestions and identify enriching SBVS targets for six of seven nuclear receptors. Further, the classifier can be extended to other proteins of interest simply by feeding new training data sets to the classifier. Our work, thus, provides a methodology to identify pharmaceutically interesting receptor conformations for nuclear receptors and other drug targets.
Ensemble docking can be a successful virtual screening technique that addresses the innate conformational heterogeneity of macromolecular drug targets. Yet, lacking a method to identify a subset of conformational states that effectively segregates active and inactive small molecules, ensemble docking may result in the recommendation of a large number of false positives. Here, three knowledge-based methods that construct structural ensembles for virtual screening are presented. Each method selects ensembles by optimizing an objective function calculated using the receiver operating characteristic (ROC) curve: either the area under the ROC curve (AUC) or a ROC enrichment factor (EF). As the number of receptor conformations, N, becomes large, the methods differ in their asymptotic scaling. Given a set of small molecules with known activities and a collection of target conformations, the most resource intense method is guaranteed to find the optimal ensemble but scales as O(2(N)). A recursive approximation to the optimal solution scales as O(N(2)), and a more severe approximation leads to a faster method that scales linearly, O(N). The techniques are generally applicable to any system, and we demonstrate their effectiveness on the androgen nuclear hormone receptor (AR), cyclin-dependent kinase 2 (CDK2), and the peroxisome proliferator-activated receptor δ (PPAR-δ) drug targets. Conformations that consisted of a crystal structure and molecular dynamics simulation cluster centroids were used to form AR and CDK2 ensembles. Multiple available crystal structures were used to form PPAR-δ ensembles. For each target, we show that the three methods perform similarly to one another on both the training and test sets.
In this Viewpoint, we provide a commentary on the impact of the Journal of Chemical Information and Modeling Special Issue on Women in Computational Chemistry published in May 2019 and the feedback we received.
Recent outbreaks of highly pathogenic and occasional drug-resistant influenza strains have highlighted the need to develop novel anti-influenza therapeutics. Here, we report computational and experimental efforts to identify influenza neuraminidase inhibitors from among the 3000 natural compounds in the Malaysian-Plants Natural-Product (NADI) database. These 3000 compounds were first docked into the neuraminidase active site. The five plants with the largest number of top predicted ligands were selected for experimental evaluation. Twelve specific compounds isolated from these five plants were shown to inhibit neuraminidase, including two compounds with IC50 values less than 92 μM. Furthermore, four of the 12 isolated compounds had also been identified in the top 100 compounds from the virtual screen. Together, these results suggest an effective new approach for identifying bioactive plant species that will further the identification of new pharmacologically active compounds from diverse natural-product resources.