Virtual screening (VS) has emerged in drug discovery as a powerful computational approach to screen large libraries of small molecules for new hits with desired properties that can then be tested experimentally. Similar to other computational approaches, VS intention is not to replace in vitro or in vivo assays, but to speed up the discovery process, to reduce the number of candidates to be tested experimentally, and to rationalize their choice. Moreover, VS has become very popular in pharmaceutical companies and academic organizations due to its time-, cost-, resources-, and labor-saving. Among the VS approaches, quantitative structure–activity relationship (QSAR) analysis is the most powerful method due to its high and fast throughput and good hit rate. As the first preliminary step of a QSAR model development, relevant chemogenomics data are collected from databases and the literature. Then, chemical descriptors are calculated on different levels of representation of molecular structure, ranging from 1D to nD, and then correlated with the biological property using machine learning techniques. Once developed and validated, QSAR models are applied to predict the biological property of novel compounds. Although the experimental testing of computational hits is not an inherent part of QSAR methodology, it is highly desired and should be performed as an ultimate validation of developed models. In this mini-review, we summarize and critically analyze the recent trends of QSAR-based VS in drug discovery and demonstrate successful applications in identifying perspective compounds with desired properties. Moreover, we provide some recommendations about the best practices for QSAR-based VS along with the future perspectives of this approach.

Schistosomiasis is a debilitating neglected tropical disease, caused by flatworms of Schistosoma genus. The treatment relies on a single drug, praziquantel (PZQ), making the discovery of new compounds extremely urgent. In this work, we integrated QSAR-based virtual screening (VS) of Schistosoma mansoni thioredoxin glutathione reductase (SmTGR) inhibitors and high content screening (HCS) aiming to discover new antischistosomal agents. Initially, binary QSAR models for inhibition of SmTGR were developed and validated using the Organization for Economic Co-operation and Development (OECD) guidance. Using these models, we prioritized 29 compounds for further testing in two HCS platforms based on image analysis of assay plates. Among them, 2-[2-(3-methyl-4-nitro-5-isoxazolyl)vinyl]pyridine and 2-(benzylsulfonyl)-1,3-benzothiazole, two compounds representing new chemical scaffolds have activity against schistosomula and adult worms at low micromolar concentrations and therefore represent promising antischistosomal hits for further hit-to-lead optimization.

The main protease (M^pro) of the SARS‐CoV‐2 has been proposed as one of the major drug targets for COVID‐19. We have identified the experimental data on the inhibitory activity of compounds tested against the closely related (96 % sequence identity, 100 % active site conservation) M^pro of SARS‐CoV. We developed QSAR models of these inhibitors and employed these models for virtual screening of all drugs in the DrugBank database. Similarity searching and molecular docking were explored in parallel, but docking failed to correctly discriminate between experimentally active and inactive compounds, so it was not relied upon for prospective virtual screening. Forty‐two compounds were identified by our models as consensus computational hits. Subsequent to our computational studies, NCATS reported the results of experimental screening of their drug collection in SARS‐CoV‐2 cytopathic effect assay (https://opendata.ncats.nih.gov/covid19/). Coincidentally, NCATS tested 11 of our 42 hits, and three of them, cenicriviroc (AC₅₀ of 8.9 μM), proglumetacin (tested twice independently, with AC₅₀ of 8.9 μM and 12.5 μM), and sufugolix (AC₅₀ 12.6 μM), were shown to be active. These observations support the value of our modeling approaches and models for guiding the experimental investigations of putative anti‐COVID‐19 drug candidates. All data and models used in this study are publicly available via Supplementary Materials, GitHub (https://github.com/alvesvm/sars‐cov‐mpro), and Chembench web portal (https://chembench.mml.unc.edu/).

Drug discovery is mostly guided by innovative and knowledge by the application of experimental and computational approaches. Quantitative structure-activity relationships (QSAR) have a critical task in the discovery and optimization of lead compounds, thereby contributing to the development of new chemical entities. 3D-QSAR methods use the information of the tridimensional molecular structure of ligands and can be applied to elucidate the relationships between 3D molecular interactions and their measured biological property, therefore, providing a rational approach for the development of new potential compounds. The purpose of this review is to provide a perspective of the utility of 3DQSAR approaches in drug design, focusing on progress, challenges and future orientations. The essential steps involved to generate reliable and predictive CoMFA models are discussed. Moreover, we present an example of application of a CoMFA study to derive 3D-QSAR models for a series of oxadiazoles inhibitors of Schistosoma mansoni Thioredoxin Glutathione Reductase (SmTGR).

Schistosomiasis is a neglected tropical disease that affects millions of people worldwide. Thioredoxin glutathione reductase of Schistosoma mansoni (SmTGR) is a validated drug target that plays a crucial role in the redox homeostasis of the parasite. We report the discovery of new chemical scaffolds against S. mansoni using a combi-QSAR approach followed by virtual screening of a commercial database and confirmation of top ranking compounds by in vitro experimental evaluation with automated imaging of schistosomula and adult worms. We constructed 2D and 3D quantitative structure–activity relationship (QSAR) models using a series of oxadiazoles-2-oxides reported in the literature as SmTGR inhibitors and combined the best models in a consensus QSAR model. This model was used for a virtual screening of Hit2Lead set of ChemBridge database and allowed the identification of ten new potential SmTGR inhibitors. Further experimental testing on both shistosomula and adult worms showed that 4-nitro-3,5-bis(1-nitro-1H-pyrazol-4-yl)-1H-pyrazole (LabMol-17) and 3-nitro-4-{[(4-nitro-1,2,5-oxadiazol-3-yl)oxy]methyl}-1,2,5-oxadiazole (LabMol-19), two compounds representing new chemical scaffolds, have high activity in both systems. These compounds will be the subjects for additional testing and, if necessary, modification to serve as new schistosomicidal agents.