Search

In the United States, a pre-market regulatory submission for any medical device that comes into contact with either a patient or the clinical practitioner must include an adequate toxicity evaluation of chemical substances that can be released from the device during its intended use. These substances, also referred to as extractables and leachables, must be evaluated for their potential to induce sensitization/allergenicity, which traditionally has been done in animal assays such as the guinea pig maximization test (GPMT). However, advances in basic and applied science are continuously presenting opportunities to employ new approach methodologies, including computational methods which, when qualified, could replace animal testing methods to support regulatory submissions. Herein, we developed a new computational tool for rapid and accurate prediction of the GPMT outcome that we have named PreS/MD (predictor of sensitization for medical devices). To enable model development, we (1) collected, curated, and integrated the largest publicly available dataset for GPMT results; (2) succeeded in developing externally predictive (balanced accuracy of 70%–74% as evaluated by both 5-fold external cross-validation and testing of novel compounds) quantitative structure-activity relationships (QSAR) models for GPMT using machine learning algorithms, including deep learning; and (3) developed a publicly accessible web portal integrating PreS/MD models that can predict GPMT outcomes for any molecule of interest. We expect that PreS/MD will be used by both industry and regulatory scientists in medical device safety assessments and help replace, reduce, or refine the use of animals in toxicity testing. PreS/MD is freely available at https://presmd.mml.unc.edu/.

Background: Chemoinformatics has several applications in the field of drug design, helping to identify new compounds against a range of ailments. Among these are Leishmaniasis, effective treatments for which are currently limited. Objective: To construct new indole 2-aminothiophene molecules using computational tools and to test their effectiveness against Leishmania amazonensis (sp.). Methods: Based on the chemical structure of thiophene-indol hybrids, we built regression models and performed molecular docking, and used these data as bases for design of 92 new molecules with predicted pIC50 and molecular docking. Among these, six compounds were selected for the synthesis and to perform biological assays (leishmanicidal activity and cytotoxicity). Results: The prediction models and docking allowed inference of characteristics that could have positive influences on the leishmanicidal activity of the planned compounds. Six compounds were synthesized, one-third of which showed promising antileishmanial activities, with IC50 ranging from 2.16 and 2.97 μM (against promastigote forms) and 0.9 and 1.71 μM (against amastigote forms), with selectivity indexes (SI) of 52 and 75. Conclusion: These results demonstrate the ability of Quantitative Structure-Activity Relationship (QSAR)-based rational drug design to predict molecules with promising leishmanicidal potential, and confirming the potential of thiophene-indole hybrids as potential new leishmanial agents.

Modern chemical toxicology is facing a growing need to Reduce, Refine, and Replace animal tests (Russell 1959) for hazard identification. The most common type of animal assays for acute toxicity assessment of chemicals used as pesticides, pharmaceuticals, or in cosmetic products is known as a "6-pack" battery of tests, including three topical (skin sensitization, skin irritation and corrosion, and eye irritation and corrosion) and three systemic (acute oral toxicity, acute inhalation toxicity, and acute dermal toxicity) end points.We compiled, curated, and integrated, to the best of our knowledge, the largest publicly available data sets and developed an ensemble of quantitative structure-activity relationship (QSAR) models for all six end points. All models were validated according to the Organisation for Economic Co-operation and Development (OECD) QSAR principles, using data on compounds not included in the training sets.In addition to high internal accuracy assessed by cross-validation, all models demonstrated an external correct classification rate ranging from 70% to 77%. We established a publicly accessible Systemic and Topical chemical Toxicity (STopTox) web portal (https://stoptox.mml.unc.edu/) integrating all developed models for 6-pack assays.We developed STopTox, a comprehensive collection of computational models that can be used as an alternative to

Chemically induced toxicity is the leading cause of recent extinction of honey bees. In this regard, we developed an innovative artificial intelligence-based web app (BeeToxAI) for assessing the acute toxicity of chemicals to Apis mellifera. Initially, we developed and externally validated QSAR models for classification (external set accuracy ∼91%) through the combination of Random Forest and molecular fingerprints to predict the potential for chemicals to cause acute contact toxicity and acute oral toxicity to honey bees. Then, we developed and externally validated regression QSAR models ( R 2 = 0.75) using Feedforward Neural Networks (FNNs). Afterward, the best models were implemented in the publicly available BeeToxAI web app ( http://beetoxai.labmol.com.br/ ) . The outputs of BeeToxAI are: toxicity predictions with estimated confidence, applicability domain estimation, and color-coded maps of relative structure fragment contributions to toxicity. As an additional assessment of BeeToxAI performance, we collected an external set of pesticides with known bee toxicity that were not included in our modeling dataset. BeeToxAI classification models were able to predict four out of five pesticides correctly. The acute contact toxicity model correctly predicted all of the eight pesticides. Here we demonstrate that BeeToxAI can be used as a rapid new approach methodology for predicting acute toxicity of chemicals in honey bees.

Recently, we have introduced fast-calculated and reliable statistical criteria to estimate whether a predictive QSAR model can be built for a given chemical dataset. These modelability criteria were successfully applied to more than 100 datasets with a binary response variable. In this study, we have extended the modelability approach to datasets with a continuous response variable. The QSARome datasets with built k Nearest Neighbors (kNN) and Random Forest (RF) QSAR models for different IC50 values were used as the training set. Two modeling criteria were proposed based on our preliminary results: MODI_ssR2 and MODI_q2. Both criteria were based on the kNN approach without variable selection and optimization of the number of nearest neighbors. MODI_ssR2 was based on leave-group-out (20%-out or fivefold) external cross-validation. MODI_q2 values were calculated as q2 values for the entire dataset (using leave-one-out cross-validation). Modeling criteria were calculated for 1, 3 and 5 nearest neighbors. QSAR models were built using 34 datasets with different methodologies and molecular descriptors. The highest coefficients of determination were found between QSAR_R2 and the MODI_ssR2 criterion calculated for k = 5 nearest neighbors: 0.91, 0.90, and 0.90, respectively, for training, test, and both training and test sets together (entire dataset). Four additional external datasets were employed to validate our approach. Our results demonstrate that the modelability criteria can be used as a good estimate of the predictive power of QSAR models even prior to model building. They establish the lower bound of the predictive power of QSAR models in terms of QSAR_R2. Regression equations between QSAR_R2 and modelability criteria can also be used for estimation of the predictive power of QSAR models for datasets not modeled yet. The use of modelability criteria can eliminate the futile attempts to build QSAR models for “tough” datasets.

Parasitic diseases remain as unresolved health issues worldwide. While for some parasites the treatments involve drug combinations with serious side effects, for others, chemical therapies are inefficient due to the emergence of drug resistance. This urges the search for novel antiparasitic agents able to act through multiple mechanisms of action. Here, we report the first multi-target model based on quantitative structure-activity relationships and a multilayer perceptron neural network (mt-QSAR-MLP) to virtually design and predict versatile inhibitors of proteins involved in the survival and/or infectivity of different pathogenic parasites. The mt-QSAR-MLP model exhibited high accuracy (>80%) in both training and test sets for the classification/prediction of protein inhibitors. Several fragments were directly extracted from the physicochemical and structural interpretations of the molecular descriptors in the mt-QSAR-MLP model. Such interpretations enabled the generation of four molecules that were predicted as multi-target inhibitors against at least three of the five parasitic proteins reported here with two of the molecules being predicted to inhibit all the proteins. Docking calculations converged with the mt-QSAR-MLP model regarding the multi-target profile of the designed molecules. The designed molecules exhibited drug-like properties, complying with Lipinski’s rule of five, as well as Ghose’s filter and Veber’s guidelines.

We review the development and application of the Simplex approach for the solution of various QSAR/QSPR problems. The general concept of the simplex method and its varieties are described. The advantages of utilizing this methodology, especially for the interpretation of QSAR/QSPR models, are presented in comparison to other fragmentary methods of molecular structure representation. The utility of SiRMS is demonstrated not only in the standard QSAR/QSPR applications, but also for mixtures, polymers, materials, and other complex systems. In addition to many different types of biological activity (antiviral, antimicrobial, antitumor, psychotropic, analgesic, etc.), toxicity and bioavailability, the review examines the simulation of important properties, such as water solubility, lipophilicity, as well as luminescence, and thermodynamic properties (melting and boiling temperatures, critical parameters, etc.). This review focuses on the stereochemical description of molecules within the simplex approach and details the possibilities of universal molecular stereo-analysis and stereochemical configuration description, along with stereo-isomerization mechanism and molecular fragment “topography” identification.

Since 2009, animal testing for cosmetic products has been prohibited in Europe, and in 2016, US EPA announced their intent to modernize the so-called "6-pack" of acute toxicity tests (acute oral toxicity, acute dermal toxicity, acute inhalation toxicity, skin irritation and corrosion, eye irritation and corrosion, and skin sensitization) and expand acceptance of alternative methods to reduce animal testing of pesticides. We have compiled, curated, and integrated the largest publicly available dataset and developed an ensemble of QSAR models for all six endpoints. All models were validated according to the OECD QSAR principles and tested using newly identified data on compounds not included in the training sets. We have established a publicly accessible Systemic and Topical chemical Toxicity (STopTox) web portal (https://stoptox.mml.unc.edu/) integrating all developed models for “6-pack” assays. This portal can be used by scientists and regulators to identify putative toxicants or non-toxicants in chemical libraries of interest.

Pesticides are used to control and combat insects and pests in the agricultural sector, households, and public health programs. The frequent and disorderly use of these pesticides may lead to variety of undesired effects. Therefore, natural products have many advantages over to synthetic compounds to be used as insecticides. The goal of this study was to find natural products with insecticidal potential against Musca domestica and Mythimna separata. To achieve this goal, we developed predictive QSAR models using MuDRA, PLS, and RF approaches and performed virtual screening of 117 natural products. As a result of QSAR modeling, we formulated the recommendations regarding physico-chemical characteristics for promising compounds active against Musca domestica and Mythimna separata. Homology models were successfully built for both species and molecular docking of QSAR hits vs known insecticides allowed us to prioritize twenty-two compounds against Musca domestica and six against Mythimna separata. Our results suggest that pimarane diterpenes, abietanes diterpenes, dimeric diterpenes and scopadulane diterpenes obtained from aerial parts of species of the genus Calceolaria (Calceolariaceae: Scrophulariaceae) can be considered as potential insecticidal.

Safety assessment is an essential component of the regulatory acceptance of industrial chemicals. Previously, we have developed a model to predict the skin sensitization potential of chemicals for two assays, the human patch test and murine local lymph node assay, and implemented this model in a web portal. Here, we report on the substantially revised and expanded freely available web tool, Pred-Skin version 3.0. This up-to-date version of Pred-Skin incorporates multiple quantitative structure-activity relationship (QSAR) models developed with in vitro, in chemico, and mice and human in vivo data, integrated into a consensus naive Bayes model that predicts human effects. Individual QSAR models were generated using skin sensitization data derived from human repeat insult patch tests, human maximization tests, and mouse local lymph node assays. In addition, data for three validated alternative methods, the direct peptide reactivity assay, KeratinoSens, and the human cell line activation test, were employed as well. Models were developed using open-source tools and rigorously validated according to the best practices of QSAR modeling. Predictions obtained from these models were then used to build a naive Bayes model for predicting human skin sensitization with the following external prediction accuracy: correct classification rate (89%), sensitivity (94%), positive predicted value (91%), specificity (84%), and negative predicted value (89%). As an additional assessment of model performance, we identified 11 cosmetic ingredients known to cause skin sensitization but were not included in our training set, and nine of them were accurately predicted as sensitizers by our models. Pred-Skin can be used as a reliable alternative to animal tests for predicting human skin sensitization.

Background: Organocalcogens are a class of organic compounds obtained by the synthesis experiments to include S, Se, or Te. Among the elements that comprise this class, Se is characterized as an essential mineral and nutrient for humans. Se has been widely studied in many aspects. Organic synthesis of organoselenides is used for obtaining new potential drug candidates and may be highly beneficial from the use of computational approaches to reduce time and cost of the experiments. Thus, the goal of our study is to evaluate the computational approaches used in the organoselenides research from 1999 to 2019. Methods: A literature review was performed by searching the database “Web of Sciences”. Results: Most of the theoretical studies included structural elucidation or structure-property analysis. We also found research regarding molecular docking approaches and Quantitative Structure-Activity Relationship (QSAR) studies. Conclusion: Computational studies have been widely applied to organoselenides. They demonstrated promising results and resulted in reduced the cost of research, increased efficacy, and, ultimately, novel organoselenides with desired properties.

The goal of this study was to perform in silico identification of bioinsecticidal potential of 42 monoterpenes against

Quantitative Structure-Activity Relationships (QSAR) models are often seen as a “black box” because they are considered difficult to interpret. Meanwhile, qualitative approaches, e.g., structural alerts (SA) or read-across, provide mechanistic insight, which is preferred for regulatory purposes, but predictive accuracy of such approaches is often low. Herein, we introduce the Chemistry-Wide Association Study (CWAS) approach, a novel framework that both addresses such deficiencies and combines advantages of statistical QSAR and alert-based approaches. The CWAS framework consists of the following steps: (i) QSAR model building for an endpoint of interest; (ii) identification of key chemical features; (iii) determination of communities of such features disproportionately co-occurring more frequently in the active than in the inactive class; and (iv) assembling these communities to form larger (and not necessarily chemically connected) novel structural alerts with high specificity. As a proof-of-concept, we have applied CWAS to model Ames mutagenicity and Stevens-Johnson Syndrome (SJS). For the well-studied Ames mutagenicity dataset, we have identified 76 important individual fragments and assembled co-occurring fragments into SA both replicative of known as well as representing novel mutagenicity alerts. For the SJS dataset, we identified 29 important fragments and assembled co-occurring communities into SA including both known and novel alerts. In summary, we demonstrate that CWAS provides a new framework to interpret predictive QSAR models and derive refined structural alerts for more effective design and safety assessment of drugs and drug candidates.

New anti-tuberculosis (anti-TB) drugs are urgently needed to battle drug-resistant Mycobacterium tuberculosis strains and to shorten the current 6-12-month treatment regimen. In this work, we have continued the efforts to develop chalcone-based anti-TB compounds by using an in silico design and QSAR-driven approach. Initially, we developed SAR rules and binary QSAR models using literature data for targeted design of new heteroaryl chalcone compounds with anti-TB activity. Using these models, we prioritized 33 compounds for synthesis and biological evaluation. As a result, 10 heteroaryl chalcone compounds (4, 8, 9, 11, 13, 17-20, and 23) were found to exhibit nanomolar activity against replicating mycobacteria, low micromolar activity against nonreplicating bacteria, and nanomolar and micromolar against rifampin (RMP) and isoniazid (INH) monoresistant strains (rRMP and rINH) (

Chemically induced skin sensitization is a complex immunological disease with a profound impact on quality of life and working ability. Despite some progress in developing alternative methods for assessing the skin sensitization potential of chemical substances, there is no in vitro test that correlates well with human data. Computational QSAR models provide a rapid screening approach and contribute valuable information for the assessment of chemical toxicity. We describe the development of a freely accessible web-based and mobile application for the identification of potential skin sensitizers. The application is based on previously developed binary QSAR models of skin sensitization potential from human (109 compounds) and murine local lymph node assay (LLNA, 515 compounds) data with good external correct classification rate (0.70-0.81 and 0.72-0.84, respectively). We also included a multiclass skin sensitization potency model based on LLNA data (accuracy ranging between 0.73 and 0.76). When a user evaluates a compound in the web app, the outputs are (i) binary predictions of human and murine skin sensitization potential; (ii) multiclass prediction of murine skin sensitization; and (iii) probability maps illustrating the predicted contribution of chemical fragments. The app is the first tool available that incorporates quantitative structure-activity relationship (QSAR) models based on human data as well as multiclass models for LLNA. The Pred-Skin web app version 1.0 is freely available for the web, iOS, and Android (in development) at the LabMol web portal ( http://labmol.com.br/predskin/ ), in the Apple Store, and on Google Play, respectively. We will continuously update the app as new skin sensitization data and respective models become available.

Prediction of chemical bioactivity and physical properties has been one of the most important applications of statistical and more recently, machine learning and artificial intelligence methods in chemical sciences. This field of research, broadly known as quantitative structure-activity relationships (QSAR) modeling, has developed many important algorithms and has found a broad range of applications in physical organic and medicinal chemistry in the past 55+ years. This Perspective summarizes recent technological advances in QSAR modeling but it also highlights the applicability of algorithms, modeling methods, and validation practices developed in QSAR to a wide range of research areas outside of traditional QSAR boundaries including synthesis planning, nanotechnology, materials science, biomaterials, and clinical informatics. As modern research methods generate rapidly increasing amounts of data, the knowledge of robust data-driven modelling methods professed within the QSAR field can become essential for scientists working both within and outside of chemical research. We hope that this contribution highlighting the generalizable components of QSAR modeling will serve to address this challenge.

Data pre-processing that includes removal of descriptors with low variance is a standard first step in quantitative structure–activity relationship modeling. In this paper, we study low-variance descriptors and show that some of them contain significant amounts of useful information. In particular, we define the notion of extreme descriptors (those variables that have the same value for almost all compounds and only a few values that are different from the common median). We show that extreme descriptors can be helpful for activity prediction in a standard binary classification setting. Moreover, we demonstrate using two case studies (M2 muscarinic receptors and skin sensitization) that extreme descriptors can be used for the identification of possibly mislabeled compounds. Because of these previously unknown, but important, properties, extreme descriptors should be considered in quantitative structure–activity relationship modeling studies. Copyright © 2016 John Wiley & Sons, Ltd.

The main protease (Mpro) 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) Mpro 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 (AC50 of 8.9 μM), proglumetacin (tested twice independently, with AC50 of 8.9 μM and 12.5 μM), and sufugolix (AC50 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/).

The poor aqueous solubility of active pharmaceutical ingredients (APIs) places a limit on their therapeutic potential. Cyclodextrins (CDs) have been shown to improve the solubility of APIs, but the magnitude of the improvement depends on the structure of both the CDs and APIs. We have developed quantitative structure&ndash, property relationship (QSPR) models that predict the stability of the complexes formed by a popular poorly soluble antibiotic, cefuroxime axetil (CA) and different CDs. We applied this model to five CA&ndash, CD systems not included in the modeling set. Two out of three systems predicted to have poor stability and poor CA solubility, and both CA&ndash, CD systems predicted to have high stability and high CA solubility were confirmed experimentally. One of the CDs that significantly improved CA solubility, methyl-&beta, CD, is described here for the first time, and we propose this CD as a novel promising excipient. Computational approaches and models developed and validated in this study could help accelerate the development of multifunctional CDs-based formulations.

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.

Multiple approaches to quantitative structure-activity relationship (QSAR) modeling using various statistical or machine learning techniques and different types of chemical descriptors have been developed over the years. Oftentimes models are used in consensus to make more accurate predictions at the expense of model interpretation. We propose a simple, fast, and reliable method termed Multi-Descriptor Read Across (MuDRA) for developing both accurate and interpretable models. The method is conceptually related to the well-known kNN approach but uses different types of chemical descriptors simultaneously for similarity assessment. To benchmark the new method, we have built MuDRA models for six different end points (Ames mutagenicity, aquatic toxicity, hepatotoxicity, hERG liability, skin sensitization, and endocrine disruption) and compared the results with those generated with conventional consensus QSAR modeling. We find that models built with MuDRA show consistently high external accuracy similar to that of conventional QSAR models. However, MuDRA models excel in terms of transparency, interpretability, and computational efficiency. We posit that due to its methodological simplicity and reliable predictive accuracy, MuDRA provides a powerful alternative to a much more complex consensus QSAR modeling. MuDRA is implemented and freely available at the Chembench web portal ( https://chembench.mml.unc.edu/mudra ).

Background: Human health assessments synthesize human, animal, and mechanistic data to produce toxicity values that are key inputs to risk-based decision making. Traditional assessments are data-, ...

The Ebola virus (EBOV) causes severe human infection that lacks effective treatment. A recent screen identified a series of compounds that block EBOV-like particle entry into human cells. Using data from this screen, Quantitative Structure-Activity Relationship (QSAR) models were built and employed for virtual screening of a ~17 million compound library. Experimental testing of 102 hits yielded 14 compounds with IC(50) values under 10 μM, including several sub-micromolar inhibitors, and more than 10-fold selectivity against host cytotoxicity. These confirmed hits include FDA-approved drugs and clinical candidates with non-antiviral indications, as well as compounds with novel scaffolds and no previously known bioactivity. Five selected hits inhibited BSL-4 live-EBOV infection in a dose-dependent manner, including vindesine (0.34 μM). Additional studies of these novel anti-EBOV compounds revealed their mechanisms of action, including the inhibition of NPC1 protein, cathepsin B/L, and lysosomal function. Compounds identified in this study are among the most potent and well-characterized anti-EBOV inhibitors reported to date.

Malaria is a life-threatening infectious disease caused by parasites of the genus Plasmodium, affecting more than 200 million people worldwide every year and leading to about a half million deaths. Malaria parasites of humans have evolved resistance to all current antimalarial drugs, urging for the discovery of new effective compounds. Given that the inhibition of deoxyuridine triphosphatase of Plasmodium falciparum (PfdUTPase) induces wrong insertions in plasmodial DNA and consequently leading the parasite to death, this enzyme is considered an attractive antimalarial drug target. Using a combi-QSAR (quantitative structure-activity relationship) approach followed by virtual screening and in vitro experimental evaluation, we report herein the discovery of novel chemical scaffolds with in vitro potency against asexual blood stages of both P. falciparum multidrug-resistant and sensitive strains and against sporogonic development of P. berghei. We developed 2D- and 3D-QSAR models using a series of nucleosides reported in the literature as PfdUTPase inhibitors. The best models were combined in a consensus approach and used for virtual screening of the ChemBridge database, leading to the identification of five new virtual PfdUTPase inhibitors. Further in vitro testing on P. falciparum multidrug-resistant (W2) and sensitive (3D7) parasites showed that compounds LabMol-144 and LabMol-146 demonstrated fair activity against both strains and presented good selectivity versus mammalian cells. In addition, LabMol-144 showed good in vitro inhibition of P. berghei ookinete formation, demonstrating that hit-to-lead optimization based on this compound may also lead to new antimalarials with transmission blocking activity. publishersversion published

Chagas disease is a neglected tropical disease (NTD) caused by the protozoan parasite Trypanosoma cruzi and is primarily transmitted to humans by the feces of infected Triatominae insects during their blood meal. The disease affects 6–8 million people, mostly in Latin America countries, and kills more people in the region each year than any other parasite-born disease, including malaria. Moreover, patient numbers are currently increasing in non-endemic, developed countries, such as Australia, Japan, Canada, and the United States. The treatment is limited to one drug, benznidazole, which is only effective in the acute phase of the disease and is very toxic. Thus, there is an urgent need to develop new, safer, and effective drugs against the chronic phase of Chagas disease. Using a QSAR-based virtual screening followed by in vitro experimental evaluation, we report herein the identification of novel potent and selective hits against T. cruzi intracellular stage. We developed and validated binary QSAR models for prediction of anti-trypanosomal activity and cytotoxicity against mammalian cells using the best practices for QSAR modeling. These models were then used for virtual screening of a commercial database, leading to the identification of 39 virtual hits. Further in vitro assays showed that seven compounds were potent against intracellular T. cruzi at submicromolar concentrations (EC50 30). Furthermore, other six compounds were also inside the hit criteria for Chagas disease, which presented activity at low micromolar concentrations (EC50 15). Moreover, we performed a multi-parameter analysis for the comparison of tested compounds regarding their balance between potency, selectivity, and predicted ADMET properties. In the next studies, the most promising compounds will be submitted to additional in vitro and in vivo assays in acute model of Chagas disease, and can be further optimized for the development of new promising drug candidates against this important yet neglected disease.

Cheminformatics is an interdisciplinary field between chemistry and informatics, which has evolved considerably since its inception in the 1960s. Initially, the cheminformatics community dealt primarily with practical and technical aspects of chemical structure representation, manipulation, and processing, while modern research explores a new role: the exploration and interpretation of large chemical databases and the discovery of new compounds with desired activity and safety profiles. Despite the recent release of several hallmark reviews addressing methods and application of cheminformatics written in Portuguese, so far there are no scientific articles presenting cheminformatics research to the Brazilian scientific community yet. To address this gap, we aim to introduce the field of cheminformatics to both students and researchers in a simple and didactic way by narrating important historical facts and contextualizing information within the scope of various applications.

Severe adverse drug reactions (ADRs) are the fourth leading cause of fatality in the U.S. with more than 100,000 deaths per year. As up to 30% of all ADRs are believed to be caused by drug-drug interactions (DDIs), typically mediated by cytochrome P450s, possibilities to predict DDIs from existing knowledge are important. We collected data from public sources on 1485, 2628, 4371, and 27,966 possible DDIs mediated by four cytochrome P450 isoforms 1A2, 2C9, 2D6, and 3A4 for 55, 73, 94, and 237 drugs, respectively. For each of these data sets, we developed and validated QSAR models for the prediction of DDIs. As a unique feature of our approach, the interacting drug pairs were represented as binary chemical mixtures in a 1:1 ratio. We used two types of chemical descriptors: quantitative neighborhoods of atoms (QNA) and simplex descriptors. Radial basis functions with self-consistent regression (RBF-SCR) and random forest (RF) were utilized to build QSAR models predicting the likelihood of DDIs for any pair of drug molecules. Our models showed balanced accuracy of 72-79% for the external test sets with a coverage of 81.36-100% when a conservative threshold for the model's applicability domain was applied. We generated virtually all possible binary combinations of marketed drugs and employed our models to identify drug pairs predicted to be instances of DDI. More than 4500 of these predicted DDIs that were not found in our training sets were confirmed by data from the DrugBank database.

This study summarizes the results of our recent QSAR and QSPR investigations on prediction of numerous aspects of environmental behavior of nitro compounds. In this study, we applied the QSAR/QSPR models previously developed by our group for virtual screening of energetic compounds, their precursors and other compounds containing nitro groups. To make predictions on the environmental impact of nitro compounds, we analyzed the trends in the change of the experimentally obtained and QSAR/QSPR-predicted values of aqueous solubility, lipophilicity, Ames mutagenicity, bioavailability, blood–brain barrier penetration, aquatic toxicity on T. pyriformis and acute oral toxicity on rats as a function of chemical structure of nitro compounds. All the models were developed using simplex descriptors in combination with random forest (RF) modeling techniques. We interpreted the possible environmental impact (different toxicological properties) in terms of dividing considered nitro compounds based on hydrophobic and hydrophilic characteristics and in terms of the influence of their molecular fragments that promote and interfere with toxicity. In particular, we found that, in general, the presence of amide or tertiary amine groups leads to an increase in toxicity. Also, it was predicted that compounds containing a NO2 group in the para-position of a benzene ring are more toxic than meta-isomers, which, in turn, are more toxic than ortho-isomers. In general, we concluded that hydrophobic nitroaromatic compounds, especially the ones with electron-accepting substituents, halogens and amino groups, are the most environmentally hazardous.

Computational tools are recognized to provide high-quality predictions for the assessment of chemical toxicity. In the recent years, mobile devices have become ubiquitous, allowing for the development of innovative and useful models implemented as chemical software applications. Here, we will briefly discuss this recent uptick in the development of web-based and mobile applications for chemical problems, focusing on best practices, development, usage and interpretation. As an example, we also describe two innovative apps (Pred-hERG and Pred-Skin) for chemical toxicity prediction developed in our laboratory. These applications are based on predictive quantitative structure-activity relationships (QSAR) models developed using the largest publicly available datasets of structurally diverse compounds. The developed tools ensure both highly accurate predictions and easy interpretation of the models, allowing users to discriminate potential toxicants and to purpose structural modifications to design safer chemicals.

Prediction of chemical bioactivity and physical properties has been one of the most important applications of statistical and more recently, machine learning and artificial intelligence methods in chemical sciences. This field of research, broadly known as Quantitative Structure-Activity Relationships (QSAR) modeling, has developed many important algorithms and has found a broad range of applications in physical organic and medicinal chemistry in the past 55+ years. This Perspective summarizes recent technological advances in QSAR modeling. Importantly, it also highlights the applicability of algorithms, modeling methods, and validation practices developed in QSAR to a wide range of research areas beyond traditional QSAR fields. These fields include nanotechnology, materials science, biomaterials, clinical informatics, and others. As modern research methods generate rapidly increasing amounts of data, knowledge of robust data-driven modelling methods is becoming essential for scientists in many disciplines both within and outside of chemical research. We hope that this contribution will serve to address this challenge.

Several non-cardiovascular drugs have been withdrawn from the market due to their inhibition of hERG K+ channels that can potentially lead to severe heart arrhythmia and death. As hERG safety testing is a mandatory FDA-required procedure, there is a considerable interest for developing predictive computational tools to identify and filter out potential hERG blockers early in the drug discovery process. In this study, we aimed to generate predictive and well-characterized quantitative structure–activity relationship (QSAR) models for hERG blockage using the largest publicly available dataset of 11,958 compounds from the ChEMBL database. The models have been developed and validated according to OECD guidelines using four types of descriptors and four different machine-learning techniques. The classification accuracies discriminating blockers from non-blockers were as high as 0.83–0.93 on external set. Model interpretation revealed several SAR rules, which can guide structural optimization of some hERG blockers into non-blockers. We have also applied the generated models for screening the World Drug Index (WDI) database and identify putative hERG blockers and non-blockers among currently marketed drugs. The developed models can reliably identify blockers and non-blockers, which could be useful for the scientific community. A freely accessible web server has been developed allowing users to identify putative hERG blockers and non-blockers in chemical libraries of their interest (http://labmol.farmacia.ufg.br/predherg).

We introduce a simple MODelability Index (MODI) that estimates the feasibility of obtaining predictive QSAR models (Correct Classification Rate above 0.7) for a binary dataset of bioactive compounds. MODI is defined as an activity class-weighted ratio of the number of the nearest neighbor pairs of compounds with the same activity class versus the total number of pairs. The MODI values were calculated for more than 100 datasets and the threshold of 0.65 was found to separate non-modelable from the modelable datasets.

The relationship between structure and reproductive toxicity was studied in organic compounds. The effects of molecular structure on acute toxicity among 53 organic compounds were studied using 2D simplex representations of molecular structure and the Partial Least Squares method. Suitable QSAR (Quantitative Structure Activity Relationships) models were obtained. These models were used to identify the structural fragments increasing toxicity. In addition, the relative effects of a number of physicochemical factors on changes in reproductive toxicity were evaluated using QSAR models.

We report on the prediction accuracy of ligand-based (2D QSAR) and structure-based (MedusaDock) methods used both independently and in consensus for ranking the congeneric series of ligands binding to three protein targets (UK, ERK2, and CHK1) from the CSAR 2011 benchmark exercise. An ensemble of predictive QSAR models was developed using known binders of these three targets extracted from the publicly available ChEMBL database. Selected models were used to predict the binding affinity of CSAR compounds toward the corresponding targets and rank them accordingly; the overall ranking accuracy evaluated by Spearman correlation was as high as 0.78 for UK, 0.60 for ERK2, and 0.56 for CHK1, placing our predictions in the top 10% among all the participants. In parallel, MedusaDock, designed to predict reliable docking poses, was also used for ranking the CSAR ligands according to their docking scores; the resulting accuracy (Spearman correlation) for UK, ERK2, and CHK1 were 0.76, 0.31, and 0.26, respectively. In addition, performance of several consensus approaches combining MedusaDock- and QSAR-predicted ranks altogether has been explored; the best approach yielded Spearman correlation coefficients for UK, ERK2, and CHK1 of 0.82, 0.50, and 0.45, respectively. This study shows that (i) externally validated 2D QSAR models were capable of ranking CSAR ligands at least as accurately as more computationally intensive structure-based approaches used both by us and by other groups and (ii) ligand-based QSAR models can complement structure-based approaches by boosting the prediction performances when used in consensus.

The enormous increase in the amount of publicly available chemical genomics data and the growing emphasis on data sharing and open science mandates that cheminformaticians also make their models publicly available for broad use by the scientific community. Chembench is one of the first publicly accessible, integrated cheminformatics Web portals. It has been extensively used by researchers from different fields for curation, visualization, analysis, and modeling of chemogenomics data. Since its launch in 2008, Chembench has been accessed more than 1 million times by more than 5000 users from a total of 98 countries. We report on the recent updates and improvements that increase the simplicity of use, computational efficiency, accuracy, and accessibility of a broad range of tools and services for computer-assisted drug design and computational toxicology available on Chembench. Chembench remains freely accessible at https://chembench.mml.unc.edu.

We have applied Hierarchical QSAR Technology (HiT QSAR) to the prediction of antiviral effects of paired combinations of picornavirus replication inhibitors against poliovirus 1 (Mahoney) in vitro. The inhibition from all binary combinations of eight antivirals were investigated. Simplex representation of molecular structure (SiRMS) was used for the generation of molecular descriptors of both pure compounds and all dual mixture combinations. Predictive QSAR models were obtained using the partial least squares (PLS) method. Predictive power of the developed models was validated using eightfold external cross-validation (CV, $$Q_{\text{ext}}^{ 2}$$ = 0.67–0.93). Adequate models ( $$Q_{\text{ext}}^{ 2}$$ = 0.53–0.97) were obtained in the same way for predicting measured inhibitory concentrations at other levels (i.e., IC30, IC40, IC60, IC70). The usage of predicted values of these concentrations in the framework of the feature net (FN) approach led to an insignificant increase in the quality of the obtained QSAR models ( $$Q_{\text{ext}}^{ 2}$$ = 0.71–0.94). Developed QSAR models were analyzed and interpreted so that structural fragments and components of the combination promoting the antiviral activity were determined (e.g., 2-(4-methoxyphenyl)-4,5-dihydrooxazole or the combination of N-hydroxybenzimidoyl and 3-methylisoxasole). Then the resulting consensus model was used to predict novel potent combinations of drugs. Combinations of enviroxime with pleconaril, WIN52084, and rupintrivir and the mixture of rupintrivir with disoxaril were predicted to cause the most inhibition of poliovirus 1 replication. HiT QSAR proved itself as an adequate tool for QSAR analysis of mixtures and, although the method described here is suitable only for binary mixtures, it can be easily extended for more complex combinations.

Computational models have earned broad acceptance for assessing chemical toxicity during early stages of drug discovery or environmental safety assessment. The majority of publicly available QSAR toxicity models have been developed for datasets including mostly drugs or drug-like compounds. We have evaluated and compared chemical spaces occupied by cosmetics, drugs, and pesticides, and explored whether current computational models of toxicity endpoints can be universally applied to all these chemicals. Our analysis of the chemical space overlap and applicability domain (AD) of models built previously for twenty different toxicity endpoints showed that most of these models afforded high coverage (>90%) for all three classes of compounds analyzed herein. Only T. pyriformis models demonstrated lower coverage for drugs and pesticides (38% and 54%, respectively). These results show that, for the most part, historical QSAR models built with data available for different toxicity endpoints can be used for toxicity assessment of novel chemicals irrespective of the intended commercial use; however, the AD restriction is necessary to assure the expected prediction accuracy. Local models may need to be developed to capture chemicals that appear as outliers with respect to global models.

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.

There is a growing public concern about the lack of reproducibility of experimental data published in peer-reviewed scientific literature. Herein, we review the most recent alerts regarding experimental data quality and discuss initiatives taken thus far to address this problem, especially in the area of chemical genomics. Going beyond just acknowledging the issue, we propose a chemical and biological data curation workflow that relies on existing cheminformatics approaches to flag, and when appropriate, correct possibly erroneous entries in large chemogenomics data sets. We posit that the adherence to the best practices for data curation is important for both experimental scientists who generate primary data and deposit them in chemical genomics databases and computational researchers who rely on these data for model development.

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.

Prior to using a quantitative structure activity relationship (QSAR) model for external predictions, its predictive power should be established and validated. In the absence of a true external data set, the best way to validate the predictive ability of a model is to perform its statistical external validation. In statistical external validation, the overall data set is divided into training and test sets. Commonly, this splitting is performed using random division. Rational splitting methods can divide data sets into training and test sets in an intelligent fashion. The purpose of this study was to determine whether rational division methods lead to more predictive models compared to random division. A special data splitting procedure was used to facilitate the comparison between random and rational division methods. For each toxicity end point, the overall data set was divided into a modeling set (80% of the overall set) and an external evaluation set (20% of the overall set) using random division. The modeling set was then subdivided into a training set (80% of the modeling set) and a test set (20% of the modeling set) using rational division methods and by using random division. The Kennard-Stone, minimal test set dissimilarity, and sphere exclusion algorithms were used as the rational division methods. The hierarchical clustering, random forest, and k-nearest neighbor (kNN) methods were used to develop QSAR models based on the training sets. For kNN QSAR, multiple training and test sets were generated, and multiple QSAR models were built. The results of this study indicate that models based on rational division methods generate better statistical results for the test sets than models based on random division, but the predictive power of both types of models are comparable.

This review is devoted to the critical analysis of advantages and disadvantages of existing mixture descriptors and their usage in various QSAR/QSPR tasks. We describe good practices for the QSAR modeling of mixtures, data sources for mixtures, a discussion of various mixture descriptors and their application, recommendations about proper external validation specific for mixture QSAR modeling, and future perspectives of this field. The biggest problem in QSAR of mixtures is the lack of reliable data about the mixtures' properties. Various mixture descriptors are used for the modeling of different endpoints. However, these descriptors have certain disadvantages, such as applicability only to 1 : 1 binary mixtures, and additive nature. The field of QSAR of mixtures is still under development, and existing efforts could be considered as a foundation for future approaches and studies. The usage of non-additive mixture descriptors, which are sensitive to interaction effects, in combination with best practices of QSAR model development (e.g., thorough data collection and curation, rigorous external validation, etc.) will significantly improve the quality of QSAR studies of mixtures.

Quantitative structure-activity relationship (QSAR) modeling and toxicogenomics are typically used independently as predictive tools in toxicology. In this study, we evaluated the power of several statistical models for predicting drug hepatotoxicity in rats using different descriptors of drug molecules, namely, their chemical descriptors and toxicogenomics profiles. The records were taken from the Toxicogenomics Project rat liver microarray database containing information on 127 drugs ( http://toxico.nibio.go.jp/datalist.html ). The model end point was hepatotoxicity in the rat following 28 days of continuous exposure, established by liver histopathology and serum chemistry. First, we developed multiple conventional QSAR classification models using a comprehensive set of chemical descriptors and several classification methods (k nearest neighbor, support vector machines, random forests, and distance weighted discrimination). With chemical descriptors alone, external predictivity (correct classification rate, CCR) from 5-fold external cross-validation was 61%. Next, the same classification methods were employed to build models using only toxicogenomics data (24 h after a single exposure) treated as biological descriptors. The optimized models used only 85 selected toxicogenomics descriptors and had CCR as high as 76%. Finally, hybrid models combining both chemical descriptors and transcripts were developed; their CCRs were between 68 and 77%. Although the accuracy of hybrid models did not exceed that of the models based on toxicogenomics data alone, the use of both chemical and biological descriptors enriched the interpretation of the models. In addition to finding 85 transcripts that were predictive and highly relevant to the mechanisms of drug-induced liver injury, chemical structural alerts for hepatotoxicity were identified. These results suggest that concurrent exploration of the chemical features and acute treatment-induced changes in transcript levels will both enrich the mechanistic understanding of subchronic liver injury and afford models capable of accurate prediction of hepatotoxicity from chemical structure and short-term assay results.

The Hierarchical Technology for Quantitative Structure - Activity Relationships (HiT QSAR) was applied to 95 diverse nitroaromatic compounds (including some widely known explosives) tested for their toxicity (50% inhibition growth concentration, IGC50) against the ciliate Tetrahymena pyriformis. The dataset was divided into subsets according to putative mechanisms of toxicity. Classification and Regression Trees (CART) approach implemented within HiT QSAR has been used for prediction of mechanism of toxicity for new compounds. The resulting models were shown to have ~80% accuracy for external datasets indicating that the mechanistic dataset division was sensible. Then, Partial Least Squares (PLS) statistical approach was used for the development of 2D QSAR models. Validated PLS models were explored to (i) elucidate the effects of different substituents in nitroaromatic compounds on toxicity; (ii) differentiate compounds by probable mechanisms of toxicity based on their structural descriptors; (iii) analyze the role of various physical-chemical factors responsible for compounds’ toxicity. Models were interpreted in terms of molecular fragments promoting or interfering with toxicity. It was also shown that mutual influence of substituents in benzene ring plays the determining role in toxicity variation. Although chemical mechanism based models were statistically significant and externally predictive (R2ext=0.64 for the external set of 63 nitroaromatics identified after all calculations have been completed), they were also shown to have limited coverage (57% for modeling and 76% for external set).

Background: Antiviral drugs are urgently needed for the treatment of acute and chronic diseases caused by enteroviruses such as coxsackievirus B3 (CVB3). The main goal of this study is quantitative structure–activity relationship (QSAR) analysis of anti-CVB3 activity (clinical CVB3 isolate 97927 [log IC50, µM]) and investigation of the selectivity of 25 ([biphenyloxy]propyl)isoxazoles, followed by computer-aided design and virtual screening of novel active compounds. Discussion: The 2D QSAR obtained models are quite satisfactory (R2 = 0.84–0.99, Q2 = 0.76–0.92, R2 ext = 0.62–0.79). Compounds with high antiviral activity and selectivity have to contain 5-trifluoromethyl-[1,2,4]oxadiazole or 2,4-difluorophenyl fragments. Insertion of 2,5-dimethylbenzene, napthyl and especially biphenyl substituents into investigated compounds substantially decreases both their antiviral activity and selectivity. Several compounds were proposed as a result of design and virtual screening. A high level of activity of 2-methoxy-1-phenyl-1H-imidazo[4,5-c]pyridine (sm428) was confirmed experimentally. Conclusion: Simplex representation of molecular structure allows successful QSAR analysis of anti-CVB3 activity of ([biphenyloxy]propyl)isoxazole derivatives. Two possible ways of battling CVB3 are considered as a future perspective.

This review explores the application of the Simplex representation of molecular structure (SiRMS) QSAR approach in antiviral research. We provide an introduction to and description of SiRMS, its application in antiviral research and future directions of development of the Simplex approach and the whole QSAR field. In the Simplex approach every molecule is represented as a system of different simplexes (tetratomic fragments with fixed composition, structure, chirality and symmetry). The main advantages of SiRMS are consideration of the different physical–chemical properties of atoms, high adequacy and good interpretability of models obtained and clear procedures for molecular design. The reliability of developed QSAR models as predictive virtual screening tools and their ability to serve as the basis of directed drug design was validated by subsequent synthetic and biological experiments. The SiRMS approach is realized as the complex of the computer program ‘HiT QSAR’, which is available on request.

Molecular modelers and cheminformaticians typically analyze experimental data generated by other scientists. Consequently, when it comes to data accuracy, cheminformaticians are always at the mercy of data providers who may inadvertently publish (partially) erroneous data. Thus, dataset curation is crucial for any cheminformatics analysis such as similarity searching, clustering, QSAR modeling, virtual screening, etc., especially nowadays when the availability of chemical datasets in public domain has skyrocketed in recent years. Despite the obvious importance of this preliminary step in the computational analysis of any dataset, there appears to be no commonly accepted guidance or set of procedures for chemical data curation. The main objective of this paper is to emphasize the need for a standardized chemical data curation strategy that should be followed at the onset of any molecular modeling investigation. Herein, we discuss several simple but important steps for cleaning chemical records in a database including the removal of a fraction of the data that cannot be appropriately handled by conventional cheminformatics techniques. Such steps include the removal of inorganic and organometallic compounds, counterions, salts and mixtures; structure validation; ring aromatization; normalization of specific chemotypes; curation of tautomeric forms; and the deletion of duplicates. To emphasize the importance of data curation as a mandatory step in data analysis, we discuss several case studies where chemical curation of the original “raw” database enabled the successful modeling study (specifically, QSAR analysis) or resulted in a significant improvement of model's prediction accuracy. We also demonstrate that in some cases rigorously developed QSAR models could be even used to correct erroneous biological data associated with chemical compounds. We believe that good practices for curation of chemical records outlined in this paper will be of value to all scientists working in the fields of molecular modeling, cheminformatics, and QSAR studies.

This work is devoted to the application of the random forest approach to QSAR analysis of aquatic toxicity of chemical compounds tested on Tetrahymena pyriformis. The simplex representation of the molecular structure approach implemented in HiT QSAR Software was used for descriptors generation on a two-dimensional level. Adequate models based on simplex descriptors and the RF statistical approach were obtained on a modeling set of 644 compounds. Model predictivity was validated on two external test sets of 339 and 110 compounds. The high impact of lipophilicity and polarizability of investigated compounds on toxicity was determined. It was shown that RF models were tolerant for insertion of irrelevant descriptors as well as for randomization of some part of toxicity values that were representing a "noise". The fast procedure of optimization of the number of trees in the random forest has been proposed. The discussed RF model had comparable or better statistical characteristics than the corresponding PLS or KNN models.

In this article the Hierarchical QSAR technology (HiT QSAR) has been used for consensus QSAR modeling of Acetylcholinesterase (AChE) inhibition by various organophosphate compounds. Simplex representation of molecular structure (SiRMS) and Lattice model (LM) QSAR approaches have been used for descriptors' generation. Statistical models have been obtained by partial least squares (PLS) method. Various chiral organophosphates represented by their (R)- and (S)-isomers, racemic mixtures and achiral structures have been investigated. Successful consensus model (R2=0.978) based on fourteen best QSAR models obtained using different QSAR approaches and training sets for several levels and methods of molecular structure representation (2.5D, double 2.5D, and 3D) has been used for the prediction of AChE inhibition of new compounds. In order to avoid chance correlations 1000 rounds of Y-scrambling were performed for each selected model. Leverage and ellipsoid Applicability domain (DA) approaches have been used for additional estimation of the quality of prognosis. Molecular fragments enhancing and interfering with inhibitory activity have been determined. It was revealed that atom's individuality and stereochemistry of chiral surroundings of the asymmetric atom of phosphorus play a vital role in AChE inhibition. Thus, (R)-isomers are always less active than the (S) isomers and the racemate. HiT QSAR proves to be a powerful tool for investigation of activity of stereo selective interactions and can be used in such subsequent studies.

The influence of molecular structure of 4-thiazolidone derivatives (about 70 compounds) on their antimicrobial activity has been investigated by QSAR approach based on Simplex Representation of Molecular Structure (SiRMS). This method is characterized by the absence of molecular alignment problem as well as by taking into account the different atom characteristics (partial charge, lipophilicity, etc) and the possibility of drug design. The following micro-organisms were used as reference cultures: Candida albicans S, Citrobacter freundii, Klebsiella pneumoniae, Pseudomonas aeruginosa (R and S strains) and Staphylococcus aureus MSSA. Statistical characteristics of QSAR models obtained by Partial Least Squares (PLS) method are rather good (R2=0.843–0.989, Q2=0.679–0.864, R2test=0.744–0.943). Y-scrambling test shows the possibility of chance correlations for equations developed for C. freundii and K. pneumonia. These two cultures were excluded from subsequent analysis. The molecular fragments promoting/interfering with antimicrobial activity and degree of their influence on the activities have been determined. Thus, naphthalene-substituted fragment (irrespective to its location within the molecule) decreases the antimicrobial activities significantly. The requirements to molecular design were formulated; e.g. highly active compounds must contain indole fragment. In addition, the influence of heterocyclic system evolution on antimicrobial properties of 4-thiazolidones derivatives was demonstrated. The important role of electrostatic characteristics of their molecular structure was also shown. The compounds with the highest predicted values of antimicrobial activities against several cultures were chosen as a result of virtual screening.