RESUMO
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.
Assuntos
Dipeptidil Peptidase 4 , Inibidores da Dipeptidil Peptidase IV , Descoberta de Drogas , Aprendizado de Máquina , Simulação de Dinâmica Molecular , Inibidores da Dipeptidil Peptidase IV/química , Inibidores da Dipeptidil Peptidase IV/farmacologia , Dipeptidil Peptidase 4/metabolismo , Dipeptidil Peptidase 4/química , Humanos , Estrutura Molecular , Diabetes Mellitus Tipo 2/tratamento farmacológicoRESUMO
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.
Assuntos
Algoritmos , Relação Quantitativa Estrutura-Atividade , Reprodutibilidade dos TestesRESUMO
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.
Assuntos
Aminoácidos Aromáticos , Ligação de HidrogênioRESUMO
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.
Assuntos
Algoritmos , Aprendizado de Máquina , Bases de Dados FactuaisRESUMO
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.
Assuntos
Antioxidantes/síntese química , Dapsona/química , Iminas/síntese química , Antioxidantes/química , Antioxidantes/farmacologia , Simulação por Computador , Teoria da Densidade Funcional , Iminas/química , Iminas/farmacologia , Estrutura Molecular , Relação Estrutura-AtividadeRESUMO
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.
Assuntos
Acetamidas/química , Teoria da Densidade Funcional , Enxofre/químicaRESUMO
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.