QSPR predicting the vapor pressure of pesticides into high/low volatility classes.

In this work, the vapor pressure of pesticides is employed as an indicator of their volatility potential. Quantitative Structure-Property Relationship models are established to predict the classification of compounds according to their volatility, into the high and low binary classes separated by the 1-mPa limit. A large dataset of 1005 structurally diverse pesticides with known experimental vapor pressure data at 20 °C is compiled from the publicly available Pesticide Properties DataBase (PPDB) and used for model development. The freely available PaDEL-Descriptor and ISIDA/Fragmentor molecular descriptor programs provide a large number of 19,947 non-conformational molecular descriptors that are analyzed through multivariable linear regressions and the Replacement Method technique. Through the selection of appropriate molecular descriptors of the substructure fragment type and the use of different standard classification metrics of model's quality, the classification of the structure-property relationship achieves acceptable results for discerning between the high and low volatility classes. Finally, an application of the obtained QSPR model is performed to predict the classes for 504 pesticides not having experimentally measured vapor pressures.

Resonance structure contributions, flexibility, and frontier molecular orbitals (HOMO-LUMO) of pelargonidin, cyanidin, and delphinidin throughout the conformational space: application to antioxidant and antimutagenic activities.

This research refers to the study and understanding of the conformational space of the positive-charged anthocyanidin structures in relation with the known chemical reactivities and bioactivities of these compounds. Therefore, the planar (P) and nonplanar (Z) conformers of the three hydroxylated anthocyanidins pelargonidin, cyanidin, and delphinidin were analyzed throughout the conformational space at the B3LYP/6-311 ++ G** level of theory. The outcome displayed eleven new conformers for pelargonidin, fifty-four for cyanidin, and thirty-one for delphinidin. Positive-charged quinoidal structures showed a significant statistical weight in the conformational space, thus coexisting simultaneously with other resonance structures, such that under certain reaction conditions, the anthocyanidins behave as positive-charged quinoidal structures instead of oxonium salts. The calculations of the permanent dipole moment and the polarizability showed relationships with the quantity and arrangement of hydroxyls in the structure. In addition, theoretical calculations were used to analyze the frontier molecular orbitals (HOMO-LUMO) of the three anthocyanidins. The novel conception of this work lies in the fact that dipole moment, polarizability, and HOMO-LUMO values were related to the reactivity/bioactivity of these three anthocyanidins. HOMO-LUMO energy gaps were useful to explain the antioxidant activity, while the percent atom contributions to HOMO were appropriate to demonstrate the antimutagenic activity as enzyme inhibitors, as well as the steric and electrostatic requirements to form the pharmacophore. Delphinidin was the strongest antioxidant anthocyanidin, and pelargonidin the best anthocyanidin with antimutagenic activity.

Dietary anthocyanins balance immune signs in osteoarthritis and obesity - update of human in vitro studies and clinical trials.

Anthocyanins are known to change ligand-receptor bindings, cell membrane permeability, and intracellular signaling pathways. The beneficial effects of dietary anthocyanins have been chronologically demonstrated in interventional and observational studies, including fourteen human chondrocyte studies and related cell culture assays, nineteen human clinical trials in osteoarthritis patients, seven in vivo obesity assays, nineteen in vitro assays in preadipocytes and related cells, and twenty-two clinical trials in overweight/obese subjects, which are critically discussed in this update. Strawberries, cherries, berries, pomegranate, tropical fruits, rosehip, purple rice, purple corn, red beans, and black soybean, together with cyanidin, delphinidin, malvidin, peonidin, some 3-O-glycosides, metabolites, and acylated anthocyanins from a potato cultivar have shown the best outcomes. The set of these five key tests and clinical trials, taken together, contributes to the understanding of the underlying mechanisms and pathways involved. Furthermore, this set shows the value of anthocyanins in counteracting the progression of osteoarthritis/obesity. The interplay between the inflammation of osteoarthritis and obesity, and the subsequent regulation/immunomodulation was performed through isolated and food anthocyanins. The antioxidant, anti-inflammatory, and immunomodulatory properties of anthocyanins explain the findings of the studies analyzed. However, further interventional studies should be conducted to finally establish the appropriate doses for anthocyanin supplementation, dose-response, and length of consumption, to include dietary recommendations for osteoarthritis/obese patients for preventive and management purposes.

In Silico Antiprotozoal Evaluation of 1,4-Naphthoquinone Derivatives against Chagas and Leishmaniasis Diseases Using QSAR, Molecular Docking, and ADME Approaches.

Chagas and leishmaniasis are two neglected diseases considered as public health problems worldwide, for which there is no effective, low-cost, and low-toxicity treatment for the host. Naphthoquinones are ligands with redox properties involved in oxidative biological processes with a wide variety of activities, including antiparasitic. In this work, in silico methods of quantitative structure-activity relationship (QSAR), molecular docking, and calculation of ADME (absorption, distribution, metabolism, and excretion) properties were used to evaluate naphthoquinone derivatives with unknown antiprotozoal activity. QSAR models were developed for predicting antiparasitic activity against Trypanosoma cruzi, Leishmania amazonensis, and Leishmania infatum, as well as the QSAR model for toxicity activity. Most of the evaluated ligands presented high antiparasitic activity. According to the docking results, the family of triazole derivatives presented the best affinity with the different macromolecular targets. The ADME results showed that most of the evaluated compounds present adequate conditions to be administered orally. Naphthoquinone derivatives show good biological activity results, depending on the substituents attached to the quinone ring, and perhaps the potential to be converted into drugs or starting molecules.

Ensemble learning application to discover new trypanothione synthetase inhibitors.

Trypanosomatid-caused diseases are among the neglected infectious diseases with the highest disease burden, affecting about 27 million people worldwide and, in particular, socio-economically vulnerable populations. Trypanothione synthetase (TryS) is considered one of the most attractive drug targets within the thiol-polyamine metabolism of typanosomatids, being unique, essential and druggable. Here, we have compiled a dataset of 401 T. brucei TryS inhibitors that includes compounds with inhibitory data reported in the literature, but also in-house acquired data. QSAR classifiers were derived and validated from such dataset, using publicly available and open-source software, thus assuring the portability of the obtained models. The performance and robustness of the resulting models were substantially improved through ensemble learning. The performance of the individual models and the model ensembles was further assessed through retrospective virtual screening campaigns. At last, as an application example, the chosen model-ensemble has been applied in a prospective virtual screening campaign on DrugBank 5.1.6 compound library. All the in-house scripts used in this study are available on request, whereas the dataset has been included as supplementary material.

A non-conformational QSAR study for plant-derived larvicides against Zika Aedes aegypti L. vector.

A set of 263 plant-derived compounds with larvicidal activity against Aedes aegypti L. (Diptera: Culicidae) vector is collected from the literature, and is studied by means of a non-conformational quantitative structure-activity relationships (QSAR) approach. The balanced subsets method (BSM) is employed to split the complete dataset into training, validation and test sets. From 26,775 freely available molecular descriptors, the most relevant structural features of compounds affecting the bioactivity are taken. The molecular descriptors are calculated through four different freewares, such as PaDEL, Mold2, EPI Suite and QuBiLs-MAS. The replacement method (RM) variable subset selection technique leads to the best linear regression models. A successful QSAR equation involves 7-conformation-independent molecular descriptors, fulfiling the evaluated internal (loo, l30%o, VIF and Y-randomization) and external (test set with Ntest = 65 compounds) validation criteria. The practical application of this QSAR model reveals promising predicted values for some natural compounds with unknown experimental larvicidal activity. Therefore, the present model constitutes the first one based on a large molecular set, being a useful computational tool for identifying and guiding the synthesis of new active molecules inspired by natural products.

QSAR studies of the antioxidant activity of anthocyanins.

Through experimental information available from antioxidant assays of seventeen anthocyanins, and six common anthocyanidins, quantitative structure-activity relationships (QSAR) have been established in the present work. The antioxidant bioactivity has been predicted in three different lipid environments: emulsified and bulk oil (methyl linoleate) (in vitro tests) at concentrations of 50 and 250 µM, and 50 µM of the inhibitor, respectively, and in human LDL (low-density lipoprotein; "bad cholesterol") (ex vivo test) at concentrations of 2.5, 10, and 25 µM of the inhibitor. Radical scavenging activity was predicted in the assay with the 1,1-diphenyl-2-picrylhydrazyl radical (DPPH·). The QSAR models developed for each test and concentration used allowed to obtain prospective information on the constitutional and topological molecular characteristics for anthocyanin/anthocyanidin compounds. Therefore, the antioxidant activity was predicted for twenty-one compounds with unknown experimental values, leading for some of them to a favorable predicted bioactivity.

Foodinformatics: Quantitative Structure-Property Relationship Modeling of Volatile Organic Compounds in Peppers.

The aim of this work was the foodinformatic (chemoinformatic) modeling of volatile organic compounds (VOCs) of different samples of peppers based on a quantitative structure-property relationship (QSPR) for the retention indices of 273 identified compounds. The experimental retention indices were measured by means of comprehensive two-dimensional gas chromatography combined with quadrupole-mass spectrometry (GC × GC/qMS) using the BPX5 and BP20 column coupled system. All the VOCs were represented by means of both conformation-independent molecular descriptors and molecular fingerprints calculated in the Dragon and PaDEL-Descriptor software. The dataset was divided into training, validation and test sets of molecules according to the Balanced Subsets Method (BSM). Subsequently, the V-WSP unsupervised variable reduction method was used to reduce the presence of multicollinearity, redundancy, and noise in the initial pool of 4,336 molecular descriptors and fingerprints. Using this method, a reduced pool of 1,664 was submitted to the supervised selection by means of the replacement method (RM) variable subset selection in order to define a four-descriptor model. The quality of the model was measured by means of the coefficient of determination and the root-mean-square deviation in fitting ( R train 2 = 0 . 879 and RMSD train = 72.1 ), validation ( R val 2 = 0 . 832 and RMSD val = 91.7 ), and prediction ( R test 2 = 0 . 915 and RMSD test = 55.4 ). The negligible differences among the parameters in the three sets indicate a stable and predictive QSPR model. This quantitative structure-activity relationship was developed keeping in mind the five principles defined by the Organization for Economic Co-operation and Development (OECD) to make it applicable. PRACTICAL APPLICATION: This predictive mathematical model developed from the retention indices of 273 volatile organic compounds (VOCs) detected in pepper samples could be useful for chromatographers working on the identification of other common VOCs in peppers or other foods by means of comprehensive two-dimensional gas chromatography combined with quadrupole-mass spectrometry (GC × GC/qMS) using a bi-dimensional stationary phase coupled system (BPX5 and BP20).

QSAR Study of Biologically Active Essential Oils against Beetles Infesting the Walnut in Catamarca, Argentina.

Essential oils from six species of aromatic plants collected in the Catamarca Province of Argentina were evaluated for their chemical composition and repellent and insecticidal activities against beetles of the genus Carpophilus (Coleoptera: Nitidulidae) and Oryzaephilus (Coleoptera: Silvanidae) that infest the local walnut production. Experimental data were analyzed using generalized estimating equations, with normal distribution and the identity link function. From the spectral information from the tested essential oils, we worked their molecular modeling as mixtures by developing mixture descriptors ( Dmix) that combined the molecular descriptor of each component in the mixture ( d i) and its relative concentration ( x i), i.e., Dmix = f( d i, x i). The application of chemoinformatic approaches determined that a combination of mixture descriptors related to molecular size, branchedness, charge distribution, and electronegativity were useful to explain the bioactivity profile against Carpophilus spp. and Oryzaephilus spp. The reported models were rigorously validated using stringent statistical parameters and essential oils reported with repellent activity against other beetle species from the Nitidulidae and Silvanidae families. This model confirmed each essential oil as a repellent with a comparable performance to the experimental reports.

Linear Regression QSAR Models for Polo-Like Kinase-1 Inhibitors.

A structurally diverse dataset of 530 polo-like kinase-1 (PLK1) inhibitors is compiled from the ChEMBL database and studied by means of a conformation-independent quantitative structure-activity relationship (QSAR) approach. A large number (26,761) of molecular descriptors are explored with the main intention of capturing the most relevant structural characteristics affecting the bioactivity. The structural descriptors are derived with different freeware, such as PaDEL, Mold², and QuBiLs-MAS; such descriptor software complements each other and improves the QSAR results. The best multivariable linear regression models are found with the replacement method variable subset selection technique. The balanced subsets method partitions the dataset into training, validation, and test sets. It is found that the proposed linear QSAR model improves previously reported models by leading to a simpler alternative structure-activity relationship.

Quantitative structure-activity relationship (QSAR) analysis of plant-derived compounds with larvicidal activity against Zika Aedes aegypti (Diptera: Culicidae) vector using freely available descriptors.

BACKGROUND: We have developed a quantitative structure-activity relationship (QSAR) model for predicting the larvicidal activity of 60 plant-derived molecules against Aedes aegypti L. (Diptera: Culicidae), a vector of several diseases such as dengue, yellow fever, chikungunya and Zika. The balanced subsets method (BSM) based on k-means cluster analysis (k-MCA) was employed to split the data set. The replacement method (RM) variable subset selection technique coupled with multivariable linear regression (MLR) proved to be successful for exploring 18 326 molecular descriptors and fingerprints calculated with PaDEL, Mold2 and EPI Suite open-source softwares. RESULTS: A robust QSAR model (Rtrain2=0.84, Strain = 0.20 and Rtest2=0.92, Stest = 0.23) involving five non-conformational descriptors was established. The model was validated and tested through the use of an external test set of compounds, the leave-one-out (LOO) and leave-more-out (LMO) cross-validation methods, Y-randomization and applicability domain (AD) analysis. CONCLUSION: The QSAR model surpasses previously published models based on geometrical descriptors, thereby representing a suitable tool for predicting larvicidal activity against the vector A. aegypti using a conformation-independent approach. © 2018 Society of Chemical Industry.

The conformation-independent QSPR approach for predicting the oxidation rate constant of water micropollutants.

In advanced water treatment processes, the degradation efficiency of contaminants depends on the reactivity of the hydroxyl radical toward a target micropollutant. The present study predicts the hydroxyl radical rate constant in water (k OH) for 118 emerging micropollutants, by means of quantitative structure-property relationships (QSPR). The conformation-independent QSPR approach is employed, together with a large number of 15,251 molecular descriptors derived with the PaDEL, Epi Suite, and Mold2 freewares. The best multivariable linear regression (MLR) models are found with the replacement method variable subset selection technique. The proposed five-descriptor model has the following statistics for the training set: [Formula: see text], RMS train = 0.21, while for the test set is [Formula: see text], RMS test = 0.11. This QSPR serves as a rational guide for predicting oxidation processes of micropollutants.

Conformation-Independent QSPR Approach for the Soil Sorption Coefficient of Heterogeneous Compounds.

We predict the soil sorption coefficient for a heterogeneous set of 643 organic non-ionic compounds by means of Quantitative Structure-Property Relationships (QSPR). A conformation-independent representation of the chemical structure is established. The 17,538 molecular descriptors derived with PaDEL and EPI Suite softwares are simultaneously analyzed through linear regressions obtained with the Replacement Method variable subset selection technique. The best predictive three-descriptors QSPR is developed on a reduced training set of 93 chemicals, having an acceptable predictive capability on 550 test set compounds. We also establish a model with a single optimal descriptor derived from CORAL freeware. The present approach compares fairly well with a previously reported one that uses Dragon descriptors.

Quantitative Structure-Antifungal Activity Relationships for cinnamate derivatives.

Quantitative Structure-Activity Relationships (QSAR) are established with the aim of analyzing the fungicidal activities of a set of 27 active cinnamate derivatives. The exploration of more than a thousand of constitutional, topological, geometrical and electronic molecular descriptors, which are calculated with Dragon software, leads to predictions of the growth inhibition on Pythium sp and Corticium rolfsii fungi species, in close agreement to the experimental values extracted from the literature. A set containing 21 new structurally related cinnamate compounds is prepared. The developed QSAR models are applied to predict the unknown fungicidal activity of this set, showing that cinnamates like 38, 28 and 42 are expected to be highly active for Pythium sp, while this is also predicted for 28 and 34 in C. rolfsii.

Amino acid profiles and quantitative structure-property relationship models as markers for Merlot and Torrontés wines.

Quantitative structure-property relationships (QSPRs) were applied to the aminograms obtained by HPLC in our laboratories for Torrontés and Merlot wines. Dragon theoretical descriptors were derived for a set of optimized amino acid structures with the purpose of establishing QSPR models. The statistical Replacement Method was used for designing the best multi-parametric linear regression models, which included structural features selected from a pool containing 1497 constitutional, topological, geometrical or electronic molecular descriptors. Predicted QSPR results were in good agreement with experimental amino acid profiles. The developed QSPR approach showed to be of practical value for distinguishing each wine varietal, and for calculating experimentally non-available amino acid concentrations of Torrontés and Merlot wines. It was also useful for assessing wine authenticity; the models were especially suitable for Merlot and Torrontés wines.

QSAR study for carcinogenicity in a large set of organic compounds.

In our continuing efforts to find out acceptable Absorption, Distribution, Metabolization, Elimination and Toxicity (ADMET) properties of organic compounds, we establish linear QSAR models for the carcinogenic potential prediction of 1464 compounds taken from the "Galvez data set", that include many marketed drugs. More than a thousand of geometry-independent molecular descriptors are simultaneously analyzed, obtained with the softwares E-Dragon and Recon. The variable subset selection method employed is the Replacement Method, and also the improved version Enhanced Replacement Method. The established models are properly validated through an external test set of compounds, and by means of the Leave-Group-Out Cross Validation method. In addition, we apply the Y-Randomization strategy and analyze the Applicability Domain of the developed model. Finally, we compare the results obtained in present study with the previous ones from the literature. The novelty of present work relies on the development of an alternative predictive structure-carcinogenicity relationship in a large heterogeneous set of organic compounds, by only using a reduced number of geometry independent molecular descriptors.

An integrated drug development approach applying topological descriptors.

We describe the opportunities posed by computer-assisted drug design in the light of two aspects of the current drug discovery scenario: the decline of innovation due to high attrition rates at clinical stage of development and the combinatorial explosion emerging from exponential growth of feasible small molecules and genome and proteome exploration. We present an overview of recent reports from our group in the field of rational drug development, by using topological descriptors (either alone, or in combination with different 3D approaches) and a diversity of modeling techniques such as Linear Discriminant Analysis and the Replacement Method. Modeling efforts aimed at the integrated prediction of several significant molecular properties in the field of drug discovery, such as pharmacological activity, aqueous solubility, human intestinal permeability and affinity to P-glycoprotein (ABCB1, MDR1) are reviewed. The suitability of conformation-independent descriptors to explore large chemical repositories is highlighted, as well as the opportunities posed by in silico guided drug repurposing.

Quantitative structure-activity relationships of mosquito larvicidal chalcone derivatives.

The mosquito larvicidal activities of a series of chalcones and some derivatives were subjected to a quantitative structure-activity relationship (QSAR) study, using more than a thousand constitutional, topological, geometrical, and electronic molecular descriptors calculated with Dragon software. The larvicidal activity values for 28 active compounds of the series were predicted, showing in general a good approximation to the experimental values found in the literature. Chalcones having one or both electron-rich rings showed high toxicity. However, the activity of chalcones was reduced by electron-withdrawing groups, and this was roughly diminished by derivatization of the carbonyl group. A set of six chalcones being structurally similar to some of the active ones, with a still unknown larvicidal activity, were prepared. Their activity values were predicted by applying the developed QSAR models, showing that two chalcones of such set, both 32 and 34, were expected to be highly active.

A comparative QSAR on 1,2,5-thiadiazolidin-3-one 1,1-dioxide compounds as selective inhibitors of human serine proteinases.

Selective inhibitors of target serine proteinases have a potential therapeutic role for the treatment of various inflammatory and related diseases. We develop a comparative quantitative structure-activity relationships based analysis on compounds embodying the 1,2,5-thiadiazolidin-3-one 1,1-dioxide scaffold. By means of classical Molecular Dynamics we obtain the conformation of each lowest-energy molecular structure from which we derive more than a thousand of structural descriptors necessary for building predictive QSAR models. We resort to two different modeling approaches with the purpose of testing the consistency of our results: (a) multivariable linear regressions based on the replacement method and forward stepwise regression, and (b) the calculation of flexible descriptors with the CORAL program. All the models are properly validated by means of standard procedures. The resulting QSAR models are supposed to be of great utility for the rational search and design (including synthesis and/or in vitro biochemical studies) of new effective non-peptidyl inhibitors of serine proteinases.

Advances in the replacement and enhanced replacement method in QSAR and QSPR theories.

The selection of an optimal set of molecular descriptors from a much greater pool of such regression variables is a crucial step in the development of QSAR and QSPR models. The aim of this work is to further improve this important selection process. For this reason three different alternatives for the initial steps of our recently developed enhanced replacement method (ERM) and replacement method (RM) are proposed. These approaches had previously proven to yield near optimal results with a much smaller number of linear regressions than the full search. The algorithms were tested on four different experimental data sets, formed by collections of 116, 200, 78, and 100 experimental records from different compounds and 1268, 1338, 1187, and 1306 molecular descriptors, respectively. The comparisons showed that one of the new alternatives further improves the ERM, which has shown to be superior to genetic algorithms for the selection of an optimal set of molecular descriptors from a much greater pool. The new proposed alternative also improves the simpler and the lower computational demand algorithm RM.