Rethinking the applicability domain analysis in QSAR models.

Notwithstanding the wide adoption of the OECD principles (or best practices) for QSAR modeling, disparities between in silico predictions and experimental results are frequent, suggesting that model predictions are often too optimistic. Of these OECD principles, the applicability domain (AD) estimation has been recognized in several reports in the literature to be one of the most challenging, implying that the actual reliability measures of model predictions are often unreliable. Applying tree-based error analysis workflows on 5 QSAR models reported in the literature and available in the QsarDB repository, i.e., androgen receptor bioactivity (agonists, antagonists, and binders, respectively) and membrane permeability (highest membrane permeability and the intrinsic permeability), we demonstrate that predictions erroneously tagged as reliable (AD prediction errors) overwhelmingly correspond to instances in subspaces (cohorts) with the highest prediction error rates, highlighting the inhomogeneity of the AD space. In this sense, we call for more stringent AD analysis guidelines which require the incorporation of model error analysis schemes, to provide critical insight on the reliability of underlying AD algorithms. Additionally, any selected AD method should be rigorously validated to demonstrate its suitability for the model space over which it is applied. These steps will ultimately contribute to more accurate estimations of the reliability of model predictions. Finally, error analysis may also be useful in "rational" model refinement in that data expansion efforts and model retraining are focused on cohorts with the highest error rates.

ElectroPredictor: An Application to Predict Mayr's Electrophilicity E through Implementation of an Ensemble Model Based on Machine Learning Algorithms.

Electrophilicity (E) is one of the most important parameters to understand the reactivity of an organic molecule. Although the theoretical electrophilicity index (ω) has been associated with E in a small homologous series, the use of w to predict E in a structurally heterogeneous set of compounds is not a trivial task. In this study, a robust ensemble model is created using Mayr's database of reactivity parameters. A combination of topological and quantum mechanical descriptors and different machine learning algorithms are employed for the model's development. The predictability of the model is assessed using different statistical parameters, and its validation is examined, including a training/test partition, an applicability domain, and a y-scrambling test. The global ensemble model presents a Q5-fold2 of 0.909 and a Qext2 of 0.912, demonstrating an excellent predictability performance of E values and showing that w is not a good descriptor for the prediction of E, especially for the case of neutral compounds. ElectroPredictor, a noncommercial Python application (https://github.com/mmoreno1/ElectroPredictor), is developed to predict E. QM9, a well-known large dataset containing 133885 neutral molecules, is used to perform a virtual screening (94.0% coverage). Finally, the 10 most electrophilic molecules are analyzed as possible new Mayr's electrophiles, which have not yet been experimentally tested. This study confirms the necessity to build an ensemble model using nonlinear machine learning algorithms, topographic descriptors, and separating molecules into charged and neutral compounds to predict E with precision.

Elucidating the Racemization Mechanism of Aliphatic and Aromatic Amino Acids by In Silico Tools.

The racemization of biomolecules in the active site can reduce the biological activity of drugs, and the mechanism involved in this process is still not fully comprehended. The present study investigates the impact of aromaticity on racemization using advanced theoretical techniques based on density functional theory. Calculations were performed at the ωb97xd/6-311++g(d,p) level of theory. A compelling explanation for the observed aromatic stabilization via resonance is put forward, involving a carbanion intermediate. The analysis, employing Hammett's parameters, convincingly supports the presence of a negative charge within the transition state of aromatic compounds. Moreover, the combined utilization of natural bond orbital (NBO) analysis and intrinsic reaction coordinate (IRC) calculations confirms the pronounced stabilization of electron distribution within the carbanion intermediate. To enhance our understanding of the racemization process, a thorough examination of the evolution of NBO charges and Wiberg bond indices (WBIs) at all points along the IRC profile is performed. This approach offers valuable insights into the synchronicity parameters governing the racemization reactions.

Effect of the Nucleophile's Nature on Chloroacetanilide Herbicides Cleavage Reaction Mechanism. A DFT Study.

In this study, the degradation mechanism of chloroacetanilide herbicides in the presence of four different nucleophiles, namely: Br-, I-, HS-, and S2O3-2, was theoretically evaluated using the dispersion-corrected hybrid functional wB97XD and the DGDZVP as a basis set. The comparison of computed activation energies with experimental data shows an excellent correlation (R2 = 0.98 for alachlor and 0.97 for propachlor). The results suggest that the best nucleophiles are those where a sulfur atom performs the nucleophilic attack, whereas the other species are less reactive. Furthermore, it was observed that the different R groups of chloroacetanilide herbicides have a negligible effect on the activation energy of the process. Further insights into the mechanism show that geometrical changes and electronic rearrangements contribute 60% and 40% of the activation energy, respectively. A deeper analysis of the reaction coordinate was conducted, employing the evolution chemical potential, hardness, and electrophilicity index, as well as the electronic flux. The charge analysis shows that the electron density of chlorine increases as the nucleophilic attack occurs. Finally, NBO analysis indicates that the nucleophilic substitution in chloroacetanilides is an asynchronous process with a late transition state for all models except for the case of the iodide attack, which occurs through an early transition state in the reaction.

In Silico Screening of the DrugBank Database to Search for Possible Drugs against SARS-CoV-2.

Coronavirus desease 2019 (COVID-19) is responsible for more than 1.80 M deaths worldwide. A Quantitative Structure-Activity Relationships (QSAR) model is developed based on experimental pIC50 values reported for a structurally diverse dataset. A robust model with only five descriptors is found, with values of R2 = 0.897, Q2LOO = 0.854, and Q2ext = 0.876 and complying with all the parameters established in the validation Tropsha's test. The analysis of the applicability domain (AD) reveals coverage of about 90% for the external test set. Docking and molecular dynamic analysis are performed on the three most relevant biological targets for SARS-CoV-2: main protease, papain-like protease, and RNA-dependent RNA polymerase. A screening of the DrugBank database is executed, predicting the pIC50 value of 6664 drugs, which are IN the AD of the model (coverage = 79%). Fifty-seven possible potent anti-COVID-19 candidates with pIC50 values > 6.6 are identified, and based on a pharmacophore modelling analysis, four compounds of this set can be suggested as potent candidates to be potential inhibitors of SARS-CoV-2. Finally, the biological activity of the compounds was related to the frontier molecular orbitals shapes.

Synthesis, In Vitro, and In Silico Analysis of the Antioxidative Activity of Dapsone Imine Derivatives.

Dapsone (DDS) is an antibacterial drug with well-known antioxidant properties. However, the antioxidant behavior of its derivatives has not been well explored. In the present work, the antioxidant activity of 10 dapsone derivatives 4-substituted was determined by an evaluation in two in vitro models (DPPH radical scavenging assay and ferric reducing antioxidant power). These imine derivatives 1-10 were obtained through condensation between DDS and the corresponding aromatic aldehydes 4-substuited. Three derivatives presented better results than DDS in the determination of DPPH (2, 9, and 10). Likewise, we have three compounds with better reducing activity than dapsone (4, 9, and 10). In order to be more insight, the redox process, a conceptual DFT analysis was carried out. Molecular descriptors such as electronic distribution, the total charge accepting/donating capacity (I/A), and the partial charge accepting/donating capacity (ω+/ω-) were calculated to analyze the relative donor-acceptor capacity through employing a donor acceptor map (DAM). The DFT calculation allowed us to establish a relationship between GAPHOMO-LUMO and DAM with the observed antioxidant effects. According to the results, we concluded that compounds 2 and 3 have the lowest Ra values, representing a good antioxidant behavior observed experimentally in DPPH radical capturing. On the other hand, derivatives 4, 9, and 10 display the best reducing capacity activity with the highest ω- and Rd values. Consequently, we propose these compounds as the best antireductants in our DDS imine derivative series.

Computational approaches for lead compound discovery in dipeptidyl peptidase-4 inhibition using machine learning and molecular dynamics techniques.

The prediction of possible lead compounds from already-known drugs that may present DPP-4 inhibition activity imply a advantage in the drug development in terms of time and cost to find alternative medicines for the treatment of Type 2 Diabetes Mellitus (T2DM). The inhibition of dipeptidyl peptidase-4 (DPP-4) has been one of the most explored strategies to develop potential drugs against this condition. A diverse dataset of molecules with known experimental inhibitory activity against DPP-4 was constructed and used to develop predictive models using different machine-learning algorithms. Model M36 is the most promising one based on the internal and external performance showing values of Q2CV = 0.813, and Q2EXT = 0.803. The applicability domain evaluation and Tropsha's analysis were conducted to validate M36, indicating its robustness and accuracy in predicting pIC50 values for organic molecules within the established domain. The physicochemical properties of the ligands, including electronegativity, polarizability, and van der Waals volume were relevant to predict the inhibition process. The model was then employed in the virtual screening of potential DPP4 inhibitors, finding 448 compounds from the DrugBank and 9 from DiaNat with potential inhibitory activity. Molecular docking and molecular dynamics simulations were used to get insight into the ligand-protein interaction. From the screening and the favorable molecular dynamic results, several compounds including Skimmin (pIC50 = 3.54, Binding energy = -8.86 kcal/mol), bergenin (pIC50 = 2.69, Binding energy = -13.90 kcal/mol), and DB07272 (pIC50 = 3.97, Binding energy = -25.28 kcal/mol) seem to be promising hits to be tested and optimized in the treatment of T2DM. This results imply a important reduction in cost and time on the application of this drugs because all the information about the its metabolism is already available.

Study of potential inhibition of the estrogen receptor α by cannabinoids using an in silico approach: Agonist vs antagonist mechanism.

Breast cancer is the main cancer type with more than 2.2 million cases in 2020, and is the principal cause of death in women; with 685000 deaths in 2020 worldwide. The estrogen receptor is involved at least in 70% of breast cancer diagnoses, and the agonist and antagonist properties of the drug in this receptor play a pivotal role in the control of this illness. This work evaluated the agonist and antagonist mechanisms of 30 cannabinoids by employing molecular docking and dynamic simulations. Compounds with docking scores < -8 kcal/mol were analyzed by molecular dynamic simulation at 300 ns, and relevant insights are given about the protein's structural changes, centered on the helicity in alpha-helices H3, H8, H11, and H12. Cannabicitran was the cannabinoid that presented the best relative binding-free energy (-34.96 kcal/mol), and based on rational modification, we found a new natural-based compound with relative binding-free energy (-44.83 kcal/mol) better than the controls hydroxytamoxifen and acolbifen. Structure modifications that could increase biological activity are suggested.

In Silico Searching for Alternative Lead Compounds to Treat Type 2 Diabetes through a QSAR and Molecular Dynamics Study.

Free fatty acid receptor 1 (FFA1) stimulates insulin secretion in pancreatic ß-cells. An advantage of therapies that target FFA1 is their reduced risk of hypoglycemia relative to common type 2 diabetes treatments. In this work, quantitative structure-activity relationship (QSAR) approach was used to construct models to identify possible FFA1 agonists by applying four different machine-learning algorithms. The best model (M2) meets the Tropsha's test requirements and has the statistics parameters R2 = 0.843, Q2CV = 0.785, and Q2ext = 0.855. Also, coverage of 100% of the test set based on the applicability domain analysis was obtained. Furthermore, a deep analysis based on the ADME predictions, molecular docking, and molecular dynamics simulations was performed. The lipophilicity and the residue interactions were used as relevant criteria for selecting a candidate from the screening of the DiaNat and DrugBank databases. Finally, the FDA-approved drugs bilastine, bromfenac, and fenofibric acid are suggested as potential and lead FFA1 agonists.

Searching glycolate oxidase inhibitors based on QSAR, molecular docking, and molecular dynamic simulation approaches.

Primary hyperoxaluria type 1 (PHT1) treatment is mainly focused on inhibiting the enzyme glycolate oxidase, which plays a pivotal role in the production of glyoxylate, which undergoes oxidation to produce oxalate. When the renal secretion capacity exceeds, calcium oxalate forms stones that accumulate in the kidneys. In this respect, detailed QSAR analysis, molecular docking, and dynamics simulations of a series of inhibitors containing glycolic, glyoxylic, and salicylic acid groups have been performed employing different regression machine learning techniques. Three robust models with less than 9 descriptors-based on a tenfold cross (Q2 CV) and external (Q2 EXT) validation-were found i.e., MLR1 (Q2 CV = 0.893, Q2 EXT = 0.897), RF1 (Q2 CV = 0.889, Q2 EXT = 0.907), and IBK1 (Q2 CV = 0.891, Q2 EXT = 0.907). An ensemble model was built by averaging the predicted pIC50 of the three models, obtaining a Q2 EXT = 0.933. Physicochemical properties such as charge, electronegativity, hardness, softness, van der Waals volume, and polarizability were considered as attributes to build the models. To get more insight into the potential biological activity of the compouds studied herein, docking and dynamic analysis were carried out, finding the hydrophobic and polar residues show important interactions with the ligands. A screening of the DrugBank database V.5.1.7 was performed, leading to the proposal of seven commercial drugs within the applicability domain of the models, that can be suggested as possible PHT1 treatment.