On the Use of the Metric rm 2 as an Effective Tool for Validation of QSAR Models in Computational Drug Design and Predictive Toxicology

Multivariate QSAR, similarity search and ADMET studies based in a set of methylamine derivatives described as dopamine transporter inhibitors

The dopamine transporter (DAT), responsible for the regulation of dopaminergic neurotransmission, is implicated in the etiology of several neuropsychiatric disorders which, in turn, have contributed to high rates of disability and numerous deaths in recent years, significantly impacting the global health system. Although the research for new drugs for the treatment of neuropsychiatric disorders has evolved in recent years, the availability of DAT-selective drugs that do not generate the same psychostimulant effects observed in drugs of abuse remains scarce. Therefore, we performed a QSAR study based on a dataset of 36 methylamine derivatives described as DAT inhibitors. The model was obtained based only in descriptors derived from 2D structures, and it was validated and generated satisfactory results considering the metrics used for internal and external validation. Subsequently, a virtual screening step also based on 2D similarity was performed, where it was possible to identify a total of 1157 compounds. After a series of reductions of the set using toxicity filters, applicability domain evaluation, and pharmacokinetic properties in silico assessment, seven hit compounds were selected as the most promising to be used, in future studies, as new scaffolds for the development of new DAT inhibitors.

Exploration of 2D and 3D-QSAR analysis and docking studies for novel dihydropteridone derivatives as promising therapeutic agents targeting glioblastoma

Background: Dihydropteridone derivatives represent a novel class of PLK1 inhibitors, exhibiting promising anticancer activity and potential as chemotherapeutic drugs for glioblastoma. Objective: The aim of this study is to develop 2D and 3D-QSAR models to validate the anticancer activity of dihydropteridone derivatives and identify optimal structural characteristics for the design of new therapeutic agents. Methods: The Heuristic method (HM) was employed to construct a 2D-linear QSAR model, while the gene expression programming (GEP) algorithm was utilized to develop a 2D-nonlinear QSAR model. Additionally, the CoMSIA approach was introduced to investigate the impact of drug structure on activity. A total of 200 novel anti-glioma dihydropteridone compounds were designed, and their activity levels were predicted using chemical descriptors and molecular field maps. The compounds with the highest activity were subjected to molecular docking to confirm their binding affinity. Results: Within the analytical purview, the coefficient of determination (R2) for the HM linear model is elucidated at 0.6682, accompanied by an R2 cv of 0.5669 and a residual sum of squares (S2) of 0.0199. The GEP nonlinear model delineates coefficients of determination for the training and validation sets at 0.79 and 0.76, respectively. Empirical modeling outcomes underscore the preeminence of the 3D-QSAR model, succeeded by the GEP nonlinear model, whilst the HM linear model manifested suboptimal efficacy. The 3D paradigm evinced an exemplary fit, characterized by formidable Q2 (0.628) and R2 (0.928) values, complemented by an impressive F-value (12.194) and a minimized standard error of estimate (SEE) at 0.160. The most significant molecular descriptor in the 2D model, which included six descriptors, was identified as "Min exchange energy for a C-N bond" (MECN). By combining the MECN descriptor with the hydrophobic field, suggestions for the creation of novel medications were generated. This led to the identification of compound 21E.153, a novel dihydropteridone derivative, which exhibited outstanding antitumor properties and docking capabilities. Conclusion: The development of 2D and 3D-QSAR models, along with the innovative integration of contour maps and molecular descriptors, offer novel concepts and techniques for the design of glioblastoma chemotherapeutic agents.

Quantitative structure-activity relationship and machine learning studies of 2-thiazolylhydrazone derivatives with anti-Cryptococcus neoformans activity

Cryptococcus neoformans is a fungus responsible for infections in humans with a significant number of cases in immunosuppressed patients, mainly in underdeveloped countries. In this context, the thiazolylhydrazones are a promising class of compounds with activity against C. neoformans. The understanding of the structure-activity relationship of these derivatives could lead to the design of robust compounds that could be promising drug candidates for fungal infections. Specifically, modern techniques such as 4D-QSAR and machine learning methods were employed in this work to generate two QSAR models (one 2D and one 4D) with high predictive power (r² for the test set equals to 0.934 and 0.831, respectively), and one random forest classification model was reported with Matthews correlation coefficient equals to 1 and 0.62 for internal and external validations, respectively. The physicochemical interpretation of selected models, indicated the importance of aliphatic substituents at the hydrazone moiety to antifungal activity, corroborating experimental data.Communicated by Ramaswamy H. Sarma.

Molecular Docking and QSAR Studies of Coumarin Derivatives as NMT Inhibitors: Simple Structural Features as Potential Modulators of Antifungal Activity

Background:

Treatments of fungal diseases, including Candidiasis, remain not up to scratch in spite of the mounting catalog of synthetic antifungal agents. These have served as the impetus for investigating new antifungal agents based on natural products. Consequently, genetic algorithm-multiple linear regression (GA-MLR) based QSAR (Quantitative Structure-Activity Relationship) studies of coumarin analogues along with molecular docking were carried out.

Methods:

Coumarin analogues with their MIC values were used to generate the training and test sets of compounds for QSAR models development; the analogues were also docked into the binding pocket of NMT (MyristoylCoA: protein N-myristoyltransferase).

Results and Discussion:

The statistical parameters for internal and external validation of QSAR analysis (R2 = 0.830, Q2 = 0.758, R2Pred = 0.610 and R2m overall = 0.683 ), Y Randomization, Ridge trace, VIF, tolerance and model criteria of Golbraikh and Tropsha data illustrate the robustness of the best proposed QSAR model. Most of the analogues bind to the electrostatic, hydrophobic clamp and display hydrogen bonding with amino acid residues of NMT. Interestingly, the most active coumarin analogue (MolDock score of -189.257) was docked deeply within the binding pocket of NMT, thereby displaying hydrogen bonding with Tyr107, Leu451, Leu450, Gln226, Cys393 and Leu394 amino acid residues.

Conclusion:

The combinations of descriptors from various descriptor subsets in QSAR analysis have highlighted the role of atomic properties such as polarizability and atomic van der Waals volume to explain the inhibitory activity. The models and related information may pave the way for important insight into the designing of putative NMT inhibitors for Candida albicans.

Lipoxygenase Inhibition Activity of Coumarin Derivatives-QSAR and Molecular Docking Study

Lipoxygenases (LOXs) are a family of enzymes found in plants, mammals, and microorganisms. In animals and plants, the enzyme has the capability for the peroxidation of unsaturated fatty acids. Although LOXs participate in the plant defense system, the enzyme’s metabolites can have numerous negative effects on human health. Therefore, many types of research are searching for compounds that can inhibit LOXs. The best quantitative structure–activity relationship (QSAR) model was obtained using a Genetic Algorithm (GA). Molecular docking was performed with iGEMDOCK. The inhibition of lipoxygenase was in the range of 7.1 to 96.6%, and the inhibition of lipid peroxidation was 7.0–91.0%. Among the synthesized compounds, the strongest inhibitor of soybean LOX-3 (96.6%) was found to be 3-benzoyl-7-(benzyloxy)-2H-chromen-2-one. A lipid peroxidation inhibition of 91.0% was achieved with ethyl 7-methoxy-2-oxo-2H-chromene-3-carboxylate. The docking scores for the soybean LOX-3 and human 5-LOX also indicated that this compound has the best affinity for these LOX enzymes. The best multiple linear QSAR model contains the atom-centered fragment descriptors C-06, RDF035p, and HATS8p. QSAR and molecular docking studies elucidated the structural features important for the enhanced inhibitory activity of the most active compounds, such as the presence of the benzoyl ring at the 3-position of coumarin’s core. Compounds with benzoyl substituents are promising candidates as potent lipoxygenase inhibitors.

Sugars and Sweeteners: Structure, Properties and In Silico Modeling

Several studies report the effects of excessive use of sugars and sweeteners in the diet. These include obesity, cardiac diseases, diabetes, and even lymphomas, leukemias, cancers of the bladder and brain, chronic fatigue syndrome, Parkinson's disease, Alzheimer's disease, multiple sclerosis, autism, and systemic lupus. On the other hand, each sugar and sweetener has a distinct metabolic assimilation process, and its chemical structure plays an important role in this process. Several scientific papers present the biological effects of the sugars and sweeteners in relation to their chemical structure. One important issue dealing with the sugars is the degree of similarity in their structures, focusing mostly on optical isomerism. Finding and developing new sugars and sweeteners with desired properties is an emerging research area, in which in silico approaches play an important role.

QSAR study of N-substituted oseltamivir derivatives as potent avian influenza virus H5N1 inhibitors using quantum chemical descriptors and statistical methods

In silico modelling studies were executed on thirty two N-substituted oseltamivir derivatives as inhibitors of influenza virus H5N1.

Neonicotinoid insecticide design: molecular docking, multiple chemometric approaches, and toxicity relationship with Cowpea aphids

Neonicotinoids are the fastest-growing class of insecticides successfully applied in plant protection, human and animal health care. The significant resistance increases led to the urgent need for alternative new neonicotinoids, with improved insecticidal activity. We performed molecular docking to describe a common binding mode of neonicotinoids into the nicotinic acetylcholine receptor, and to select the appropriate conformations to derive models. These were further used in a QSAR study employing both linear and nonlinear approaches to model the inhibitory activity against the Cowpea aphids. Linear modeling was performed by multiple linear regression and partial least squares and nonlinear modeling by artificial neural networks and support vector machine methods. The OECD principles were considered for QSAR models validation. Robust models with predictive power were found for neonicotinoid diverse structures. Based on our QSAR and docking outcomes, five new insecticides were predicted, according to the model applicability domain, the ligand efficiencies, and the binding mode. Therefore, the developed models can be confidently used for the prediction of the insecticidal activity of new chemicals, saving a substantial amount of time and money and, also, contributing to the chemical risk assessment.

An Approach to Identify New Insecticides Against Myzus Persicae. In silico Study Based on Linear and Non-linear Regression Techniques

Neonicotinoids are known to have high insecticidal potency, low mammalian toxicity and relatively tough activity for the development of resistance against aphids. A series of guadipyr insecticides, active against Myzus persicae was engaged in silico studies, based on Multiple Linear Regression (MLR), Partial Least Squares regression (PLS), Artificial Neural Networks (ANN), Support Vector Machine (SVM) and Pharmacophore modeling. Robust and predictive models were built using correlations between the insecticidal profile, expressed by experimental pLC₅₀ values, and molecular descriptors, calculated from the energy optimized structures. Four new potential insecticides active against Myzus persicae and their predicted pLC₅₀ toxicity values were reported for the first time. The models presented here can be used as an approach in the screening and prioritization of chemicals in a scientific and regulatory frame and for toxicity prediction.

In silico study toward the identification of new and safe potential inhibitors of photosynthetic electron transport

To address the rising global demand for food, it is necessary to search for new herbicides that can control resistant weeds. We performed a 2D-quantitative structure-activity relationship (QSAR) study to predict compounds with photosynthesis-inhibitory activity. A data set of 44 compounds (quinolines and naphthalenes), which are described as photosynthetic electron transport (PET) inhibitors, was used. The obtained model was approved in internal and external validation tests. 2D Similarity-based virtual screening was performed and 64 compounds were selected from the ZINC database. By using the VEGA QSAR software, 48 compounds were shown to have potential toxic effects (mutagenicity and carcinogenicity). Therefore, the model was also tested using a set of 16 molecules obtained by a similarity search of the ZINC database. Six compounds showed good predicted inhibition of PET. The obtained model shows potential utility in the design of new PET inhibitors, and the hit compounds found by virtual screening are novel bicyclic scaffolds of this class.

Deep Eutectic Solvents as Convenient Media for Synthesis of Novel Coumarinyl Schiff Bases and Their QSAR Studies

Deep eutectic solvents, as green and environmentally friendly media, were utilized in the synthesis of novel coumarinyl Schiff bases. Novel derivatives were synthesized from 2-((4-methyl-2-oxo-2H-chromen-7-yl)oxy)acetohydrazide and corresponding aldehyde in choline chloride:malonic acid (1:1) based deep eutectic solvent. In these reactions, deep eutectic solvent acted as a solvent and catalyst as well. Novel Schiff bases were synthesized in high yields (65–75%) with no need for further purification, and their structures were confirmed by mass spectra, ¹H and ¹³C NMR. Furthermore, their antioxidant activity was determined and compared to antioxidant activity of previously synthesized derivatives, thus investigating their structure–activity relationship utilizing quantitative structure-activity relationship QSAR studies. Calculation of molecular descriptors has been performed by DRAGON software. The best QSAR model (R_tr = 0.636; R_ext = 0.709) obtained with three descriptors (MATS3m, Mor22u, Hy) implies that the pairs of atoms higher mass at the path length 3, three-dimensional arrangement of atoms at scattering parameter s = 21 Å⁻¹, and higher number of hydrophilic groups (-OH, -NH) enhanced antioxidant activity. Electrostatic potential surface of the most active compounds showed possible regions for donation of electrons to 1,1-diphenyl-2-picryhydrazyl (DPPH) radicals.

A best​ comprehension about the toxicity of phenylsulfonyl carboxylates in Vibrio fischeri using quantitative structure activity/property relationship methods

Aromatic sulfones comprise a class of chemicals used in agrochemical and pharmaceutical industries and as floatation and extractant agents in petrochemical and metallurgy industries. In this study, new QSA(P)R studies were carried out to predict the toxicity against Vibrio fischeri of a set of 52 aromatic sulfones. The same approach was used to evaluate the relationship between these endpoint and the water solubility, another important environmental endpoint. The study resulted in models of good statistical quality and mechanistic interpretation with a possible correlation between the two endpoints, but the toxic effect is also likely to depend on other physicochemical properties. The use of the PLS2, a method not commonly used in QSA(P)R studies, also produced models of greater reliability, and the relationship between the two endpoints was reinforced to some degree. These results are useful for better understanding the process by which these compounds exert their environmental toxicity, thus aiding in the development of industrially useful compounds with less potential environmental damage.

Pharmacophore generation, atom-based 3D-QSAR, HQSAR and activity cliff analyses of benzothiazine and deazaxanthine derivatives as dual A2A antagonists/MAO‑B inhibitors

Dual inhibition of A_2A and MAO-B is an emerging strategy in neurodegenerative diseases, such as Alzheimer's disease (AD) and Parkinson's disease (PD). In this study, atom-based three-dimensional quantitative structure-activity relationship (3D-QSAR) and hologram quantitative structure-activity relationship (HQSAR) models were generated with benzothiazine and deazaxanthine derivatives. Based on activity against A_2A and MAO-B, two statistically significant 3D-QSAR models (r² = 0.96, q² = 0.76 and r² = 0.91, q² = 0.63) and HQSAR models (r² = 0.93, q² = 0.68 and r² = 0.97, q² = 0.58) were developed. In an activity cliff analysis, structural outliers were identified by calculating the Mahalanobis distance for a pair of compounds with A_2A and MAO-B inhibitory activities. The generated 3D-QSAR and HQSAR models, activity cliff analysis, molecular docking and dynamic studies for dual target protein inhibitors provide key structural scaffolds that serve as building blocks in designing drug-like molecules for neurodegenerative diseases.

Natural products as leads in schistosome drug discovery

Schistosomiasis is a neglected parasitic tropical disease that claims around 200,000 human lives every year. Praziquantel (PZQ), the only drug recommended by the World Health Organization for the treatment and control of human schistosomiasis, is now facing the threat of drug resistance, indicating the urgent need for new effective compounds to treat this disease. Therefore, globally, there is renewed interest in natural products (NPs) as a starting point for drug discovery and development for schistosomiasis. Recent advances in genomics, proteomics, bioinformatics, and cheminformatics have brought about unprecedented opportunities for the rapid and more cost-effective discovery of new bioactive compounds against neglected tropical diseases. This review highlights the main contributions that NP drug discovery and development have made in the treatment of schistosomiasis and it discusses how integration with virtual screening (VS) strategies may contribute to accelerating the development of new schistosomidal leads, especially through the identification of unexplored, biologically active chemical scaffolds and structural optimization of NPs with previously established activity.

Insight into the binding mode and the structural features of the pyrimidine derivatives as human A2A adenosine receptor antagonists

The interaction of 278 monocyclic and bicyclic pyrimidine derivatives with human A2A adenosine receptor (AR) was investigated by employing molecular dynamics, thermodynamic analysis and three-dimensional quantitative structure-activity relationship (3D-QSAR) approaches. The binding analysis reveals that the pyrimidine derivatives are anchored in TM2, 3, 5, 6 and 7 of A2A AR by the aromatic stacking and hydrogen bonding interactions. The key residues involving Phe168, Glu169, and Asn253 stabilize the monocyclic and bicyclic cores of inhibitors. The thermodynamic analysis by molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) approach also confirms the reasonableness of the binding modes. In addition, the ligand-/receptor-based comparative molecular similarity indices analysis (CoMSIA) models of high statistical significance were generated and the resulting contour maps correlate well with the structural features of the antagonists essential for high A2A AR affinity. A minor/bulky group with negative charge at C2/C6 of pyrimidine ring respectively enhances the activity for all these pyrimidine derivatives. Particularly, the higher electron density of the ring in the bicyclic derivatives, the more potent the antagonists. The obatined results might be helpful in rational design of novel candidate of A2A adenosine receptor antagonist for treatment of Parkinson's disease.

Preliminary studies on model development for rodent toxicity and its interspecies correlation with aquatic toxicities of pharmaceuticals

Environmental toxicity due to pharmaceuticals has been an issue of serious concern for long time. Development of chemometric models with reliable predictive power has been considered as an effective tool for the design of new drug agents with reduced or without ecotoxic potential. Considering a higher degree of similarity in genetic homology towards drug receptor with mammals, we have used a dataset of 194 compounds with reported rodent, fish, daphnia and algae toxicity data for extrapolation of their toxicity towards humans. Allowing for rodents as the most surrogate to human physiology, attempts have also been made to develop interspecies correlation models keeping rodent toxicity as dependent variable so that any drug without reported rodent toxicity can be predicted using fish, daphnia or algae toxicity data which can be consequently extrapolated to human toxicity. The developed models have been subjected to multiple validation strategies. Acceptable results have been obtained in both cases of direct and interspecies extrapolation quantitative structure-activity relationship models.

QSAR models for inhibitors of physiological impact of Escherichia coli that leads to diarrhea

Quantitative structure - activity relationships (QSARs) developed to evaluate percentage of inhibition of STa-stimulated (Escherichia coli) cGMP accumulation in T84 cells are calculated by the Monte Carlo method. This endpoint represents a measure of biological activity of a substance against diarrhea. Statistical quality of the developed models is quite good. The approach is tested using three random splits of data into the training and test sets. The statistical characteristics for three splits are the following: (1) n=20, r(2)=0.7208, q(2)=0.6583, s=16.9, F=46 (training set); n=11, r(2)=0.8986, s=14.6 (test set); (2) n=19, r(2)=0.6689, q(2)=0.5683, s=17.6, F=34 (training set); n=12, r(2)=0.8998, s=12.1 (test set); and (3) n=20, r(2)=0.7141, q(2)=0.6525, s=14.7, F=45 (training set); n=11, r(2)=0.8858, s=19.5 (test set). Based on the proposed here models hypothetical compounds which can be useful agents against diarrhea are suggested.