1
|
Shaping the Future of Obesity Treatment: In Silico Multi-Modeling of IP6K1 Inhibitors for Obesity and Metabolic Dysfunction. Pharmaceuticals (Basel) 2024; 17:263. [PMID: 38399478 PMCID: PMC10891520 DOI: 10.3390/ph17020263] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Revised: 02/14/2024] [Accepted: 02/16/2024] [Indexed: 02/25/2024] Open
Abstract
Recent research has uncovered a promising approach to addressing the growing global health concern of obesity and related disorders. The inhibition of inositol hexakisphosphate kinase 1 (IP6K1) has emerged as a potential therapeutic strategy. This study employs multiple ligand-based in silico modeling techniques to investigate the structural requirements for benzisoxazole derivatives as IP6K1 inhibitors. Firstly, we developed linear 2D Quantitative Structure-Activity Relationship (2D-QSAR) models to ensure both their mechanistic interpretability and predictive accuracy. Then, ligand-based pharmacophore modeling was performed to identify the essential features responsible for the compounds' high activity. To gain insights into the 3D requirements for enhanced potency against the IP6K1 enzyme, we employed multiple alignment techniques to set up 3D-QSAR models. Given the absence of an available X-ray crystal structure for IP6K1, a reliable homology model for the enzyme was developed and structurally validated in order to perform structure-based analyses on the selected dataset compounds. Finally, molecular dynamic simulations, using the docked poses of these compounds, provided further insights. Our findings consistently supported the mechanistic interpretations derived from both ligand-based and structure-based analyses. This study offers valuable guidance on the design of novel IP6K1 inhibitors. Importantly, our work exclusively relies on non-commercial software packages, ensuring accessibility for reproducing the reported models.
Collapse
|
2
|
In Silico Modeling and Structural Analysis of Soluble Epoxide Hydrolase Inhibitors for Enhanced Therapeutic Design. Molecules 2023; 28:6379. [PMID: 37687207 PMCID: PMC10490281 DOI: 10.3390/molecules28176379] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Revised: 08/17/2023] [Accepted: 08/28/2023] [Indexed: 09/10/2023] Open
Abstract
Human soluble epoxide hydrolase (sEH), a dual-functioning homodimeric enzyme with hydrolase and phosphatase activities, is known for its pivotal role in the hydrolysis of epoxyeicosatrienoic acids. Inhibitors targeting sEH have shown promising potential in the treatment of various life-threatening diseases. In this study, we employed a range of in silico modeling approaches to investigate a diverse dataset of structurally distinct sEH inhibitors. Our primary aim was to develop predictive and validated models while gaining insights into the structural requirements necessary for achieving higher inhibitory potential. To accomplish this, we initially calculated molecular descriptors using nine different descriptor-calculating tools, coupled with stochastic and non-stochastic feature selection strategies, to identify the most statistically significant linear 2D-QSAR model. The resulting model highlighted the critical roles played by topological characteristics, 2D pharmacophore features, and specific physicochemical properties in enhancing inhibitory potential. In addition to conventional 2D-QSAR modeling, we implemented the Transformer-CNN methodology to develop QSAR models, enabling us to obtain structural interpretations based on the Layer-wise Relevance Propagation (LRP) algorithm. Moreover, a comprehensive 3D-QSAR analysis provided additional insights into the structural requirements of these compounds as potent sEH inhibitors. To validate the findings from the QSAR modeling studies, we performed molecular dynamics (MD) simulations using selected compounds from the dataset. The simulation results offered crucial insights into receptor-ligand interactions, supporting the predictions obtained from the QSAR models. Collectively, our work serves as an essential guideline for the rational design of novel sEH inhibitors with enhanced therapeutic potential. Importantly, all the in silico studies were performed using open-access tools to ensure reproducibility and accessibility.
Collapse
|
3
|
Designing multi-target drugs for the treatment of major depressive disorder. Expert Opin Drug Discov 2023; 18:643-658. [PMID: 37183604 DOI: 10.1080/17460441.2023.2214361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
INTRODUCTION Major depressive disorders (MDD) pose major health burdens globally. Currently available medications have their limitations due to serious adverse effects, long latency periods as well as resistance. Considering the highly complicated pathological nature of this disorder, it has been suggested that multitarget drugs or multi-target-directed ligands (MTDLs) may provide long-term therapeutic solutions for the treatment of MDD. AREAS COVERED In the current review, recent lead design and lead modification strategies have been covered. Important investigations reported in the last ten years (2013-2022) for the pre-clinical development of MTDLs (through synthetic medicinal chemistry and biological evaluation) for the treatment of MDD were discussed as case studies to focus on the recent design strategies. The discussions are categorized based on the pharmacological targets. On the basis of these important case studies, the challenges involved in different design strategies were discussed in detail. EXPERT OPINION Even though large variations were observed in the selection of pharmacological targets, some potential biological targets (NMDA, melatonin receptors) are required to be explored extensively for the design of MTDLs. Similarly, apart from structure activity relationship (SAR), in silico techniques such as multitasking cheminformatic modelling, molecular dynamics simulation and virtual screening should be exploited to a greater extent.
Collapse
|
4
|
In silico characterization of aryl benzoyl hydrazide derivatives as potential inhibitors of RdRp enzyme of H5N1 influenza virus. Front Pharmacol 2022; 13:1004255. [PMID: 36225563 PMCID: PMC9548590 DOI: 10.3389/fphar.2022.1004255] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Accepted: 08/24/2022] [Indexed: 11/13/2022] Open
Abstract
RNA-dependent RNA polymerase (RdRp) is a potential therapeutic target for the discovery of novel antiviral agents for the treatment of life-threatening infections caused by newly emerged strains of the influenza virus. Being one of the most conserved enzymes among RNA viruses, RdRp and its inhibitors require further investigations to design novel antiviral agents. In this work, we systematically investigated the structural requirements for antiviral properties of some recently reported aryl benzoyl hydrazide derivatives through a range of in silico tools such as 2D-quantitative structure-activity relationship (2D-QSAR), 3D-QSAR, structure-based pharmacophore modeling, molecular docking and molecular dynamics simulations. The 2D-QSAR models developed in the current work achieved high statistical reliability and simultaneously afforded in-depth mechanistic interpretability towards structural requirements. The structure-based pharmacophore model developed with the docked conformation of one of the most potent compounds with the RdRp protein of H5N1 influenza strain was utilized for developing a 3D-QSAR model with satisfactory statistical quality validating both the docking and the pharmacophore modeling methodologies performed in this work. However, it is the atom-based alignment of the compounds that afforded the most statistically reliable 3D-QSAR model, the results of which provided mechanistic interpretations consistent with the 2D-QSAR results. Additionally, molecular dynamics simulations performed with the apoprotein as well as the docked complex of RdRp revealed the dynamic stability of the ligand at the proposed binding site of the receptor. At the same time, it also supported the mechanistic interpretations drawn from 2D-, 3D-QSAR and pharmacophore modeling. The present study, performed mostly with open-source tools and webservers, returns important guidelines for research aimed at the future design and development of novel anti-viral agents against various RNA viruses like influenza virus, human immunodeficiency virus-1, hepatitis C virus, corona virus, and so forth.
Collapse
|
5
|
Predicting the Surface Tension of Deep Eutectic Solvents: A Step Forward in the Use of Greener Solvents. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27154896. [PMID: 35956845 PMCID: PMC9370217 DOI: 10.3390/molecules27154896] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 07/28/2022] [Accepted: 07/28/2022] [Indexed: 11/16/2022]
Abstract
Deep eutectic solvents (DES) are an important class of green solvents that have been developed as an alternative to toxic solvents. However, the large-scale industrial application of DESs requires fine-tuning their physicochemical properties. Among others, surface tension is one of such properties that have to be considered while designing novel DESs. In this work, we present the results of a detailed evaluation of Quantitative Structure-Property Relationships (QSPR) modeling efforts designed to predict the surface tension of DESs, following the Organization for Economic Co-operation and Development (OECD) guidelines. The data set used comprises a large number of structurally diverse binary DESs and the models were built systematically through rigorous validation methods, including ‘mixtures-out’- and ‘compounds-out’-based data splitting. The most predictive individual QSPR model found is shown to be statistically robust, besides providing valuable information about the structural and physicochemical features responsible for the surface tension of DESs. Furthermore, the intelligent consensus prediction strategy applied to multiple predictive models led to consensus models with similar statistical robustness to the individual QSPR model. The benefits of the present work stand out also from its reproducibility since it relies on fully specified computational procedures and on publicly available tools. Finally, our results not only guide the future design and screening of novel DESs with a desirable surface tension but also lays out strategies for efficiently setting up silico-based models for binary mixtures.
Collapse
|
6
|
Structural behavior of monomer of SARS-CoV-2 spike protein during initial stage of adsorption on graphene. MATERIALS TODAY. CHEMISTRY 2021; 22:100572. [PMID: 34485782 PMCID: PMC8405511 DOI: 10.1016/j.mtchem.2021.100572] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2021] [Revised: 08/22/2021] [Accepted: 08/26/2021] [Indexed: 05/14/2023]
Abstract
Spike glycoprotein of the SARS-CoV-2 virus and its structure play a crucial role in the infections of cells containing angiotensin-converting enzyme 2 (ACE2) as well as in the interactions of this virus with surfaces. Protection against viruses and often even their deactivation is one of the great varieties of graphene applications. The structural changes of the non-glycosylated monomer of the spike glycoprotein trimer (denoted as S-protein in this work) triggered by its adsorption onto graphene at the initial stage are investigated by means of atomistic molecular dynamics simulations. The adsorption of the S-protein happens readily during the first 10 ns. The shape of the S-protein becomes more prolate during the adsorption, but this trend, albeit less pronounced, is observed also for the freely relaxing S-protein in water. The receptor-binding domain (RBD) of the free and adsorbed S-protein manifests itself as the most rigid fragment of the whole S-protein. The adsorption even enhances the rigidity of the whole S-protein as well as its subunits. Only one residue of the RBD involved in the specific interactions with ACE2 during the cell infection is involved in the direct contact of the adsorbed S-protein with the graphene. The new intramolecular hydrogen bonds formed during the S-protein adsorption replace the S-protein-water hydrogen bonds; this trend, although less apparent, is observed also during the relaxation of the free S-protein in water. In the initial phase, the secondary structure of the RBD fragment specifically interacting with ACE2 receptor is not affected during the S-protein adsorption onto the graphene.
Collapse
|
7
|
Multi-Target QSAR Approaches for Modeling Protein Inhibitors. Simultaneous Prediction of Activities Against Biomacromolecules Present in Gram-Negative Bacteria. Curr Top Med Chem 2016; 15:1801-13. [PMID: 25961517 DOI: 10.2174/1568026615666150506144814] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Revised: 02/25/2015] [Accepted: 03/24/2015] [Indexed: 11/22/2022]
Abstract
Drug discovery is aimed at finding therapeutic agents for the treatment of many diverse diseases and infections. However, this is a very slow an expensive process, and for this reason, in silico approaches are needed to rationalize the search for new molecular entities with desired biological profiles. Models focused on quantitative structure-activity relationships (QSAR) have constituted useful complementary tools in medicinal chemistry, allowing the virtual predictions of dissimilar pharmacological activities of compounds. In the last 10 years, multi-target (mt) QSAR models have been reported, representing great advances with respect to those models generated from classical approaches. Thus, mt- QSAR models can simultaneously predict activities against different biological targets (proteins, microorganisms, cell lines, etc.) by using large and heterogeneous datasets of chemicals. The present review is devoted to discuss the most promising mt-QSAR models, particularly those developed for the prediction of protein inhibitors. We also report the first multi-tasking QSAR (mtk-QSAR) model for simultaneous prediction of inhibitors against biomacromolecules (specifically proteins) present in Gram-negative bacteria. This model allowed us to consider both different proteins and multiple experimental conditions under which the inhibitory activities of the chemicals were determined. The mtk-QSAR model exhibited accuracies higher than 98% in both training and prediction sets, also displaying a very good performance in the classification of active and inactive cases that depended on the specific elements of the experimental conditions. The physicochemical interpretations of the molecular descriptors were also analyzed, providing important insights regarding the molecular patterns associated with the appearance/enhancement of the inhibitory potency.
Collapse
|
8
|
First Multitarget Chemo-Bioinformatic Model To Enable the Discovery of Antibacterial Peptides against Multiple Gram-Positive Pathogens. J Chem Inf Model 2016; 56:588-98. [PMID: 26960000 DOI: 10.1021/acs.jcim.5b00630] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Antimicrobial peptides (AMPs) have emerged as promising therapeutic alternatives to fight against the diverse infections caused by different pathogenic microorganisms. In this context, theoretical approaches in bioinformatics have paved the way toward the creation of several in silico models capable of predicting antimicrobial activities of peptides. All current models have several significant handicaps, which prevent the efficient search for highly active AMPs. Here, we introduce the first multitarget (mt) chemo-bioinformatic model devoted to performing alignment-free prediction of antibacterial activity of peptides against multiple Gram-positive bacterial strains. The model was constructed from a data set containing 2488 cases of AMPs sequences assayed against at least 1 out of 50 Gram-positive bacterial strains. This mt-chemo-bioinformatic model displayed percentages of correct classification higher than 90.00% in both training and prediction (test) sets. For the first time, two computational approaches derived from basic concepts in genetics and molecular biology were applied, allowing the calculations of the relative contributions of any amino acid (in a defined position) to the antibacterial activity of an AMP and depending on the bacterial strain used in the biological assay. The present mt-chemo-bioinformatic model constitutes a powerful tool to enable the discovery of potent and versatile AMPs.
Collapse
|
9
|
A general ANN-based multitasking model for the discovery of potent and safer antibacterial agents. Methods Mol Biol 2015; 1260:45-64. [PMID: 25502375 DOI: 10.1007/978-1-4939-2239-0_4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Bacteria have been one of the world's most dangerous and deadliest pathogens for mankind, nowadays giving rise to significant public health concerns. Given the prevalence of these microbial pathogens and their increasing resistance to existing antibiotics, there is a pressing need for new antibacterial drugs. However, development of a successful drug is a complex, costly, and time-consuming process. Quantitative Structure-Activity Relationships (QSAR)-based approaches are valuable tools for shortening the time of lead compound identification but also for focusing and limiting time-costly synthetic activities and in vitro/vivo evaluations. QSAR-based approaches, supported by powerful statistical techniques such as artificial neural networks (ANNs), have evolved to the point of integrating dissimilar types of chemical and biological data. This chapter reports an overview of the current research and potential applications of QSAR modeling tools toward the rational design of more efficient antibacterial agents. Particular emphasis is given to the setup of multitasking models along with ANNs aimed at jointly predicting different antibacterial activities and safety profiles of drugs/chemicals under diverse experimental conditions.
Collapse
|
10
|
|
11
|
Current tendencies in antimicrobial research: medicinal chemistry of antibacterial agents and advances in the use of computational methodologies. Curr Top Med Chem 2013; 13:3011-2. [PMID: 24200364 DOI: 10.2174/15680266113136660216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
12
|
Role of ligand-based drug design methodologies toward the discovery of new anti- Alzheimer agents: futures perspectives in Fragment-Based Ligand Design. Curr Med Chem 2012; 19:1635-45. [PMID: 22376033 DOI: 10.2174/092986712799945058] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2011] [Revised: 01/23/2012] [Accepted: 01/23/2012] [Indexed: 11/22/2022]
Abstract
Alzheimer's disease (AD), a degenerative disease affecting the brain, is the single most common source of dementia in adults. The cause and the progression of AD still remains a mystery among medical experts. As a result, a cure has not yet been discovered, even after decade's worth of research that started since 1906, when the disease was first identified. Despite the efforts of the scientific community, several of the biological receptors associated with AD have not been sufficiently studied to date, limiting in turn the design of new and more potent anti-AD agents. Thus, the search for new drug candidates as inhibitors of different targets associated with AD constitutes an essential part towards the discovery of new and more efficient anti-AD therapies. The present work is focused on the role of the Ligand-Based Drug Design (LBDD) methodologies which have been applied for the elucidation of new molecular entities with high inhibitory activity against targets related with AD. Particular emphasis is given also to the current state of fragment-based ligand approaches as alternatives of the Fragment-Based Drug Discovery (FBDD) methodologies. Finally, several guidelines are offered to show how the use of fragment-based descriptors can be determinant for the design of multi-target inhibitors of proteins associated with AD.
Collapse
|
13
|
Quantitative structure-activity relationship modelling of the carcinogenic risk of nitroso compounds using regression analysis and the TOPS-MODE approach. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2010; 21:277-304. [PMID: 20544552 DOI: 10.1080/10629361003773930] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
Worldwide, legislative and governmental efforts are focusing on establishing simple screening tools for identifying those chemicals most likely to cause adverse effects without experimentally testing all chemicals of regulatory concern. This is because even the most basic biological testing of compounds of concern, apart from requiring a huge number of test animals, would be neither resource nor time effective. Thus, alternative approaches such as the one proposed here, quantitative structure-activity relationship (QSAR) modelling, are increasingly being used for identifying the potential health hazards and subsequent regulation of new industrial chemicals. This paper follows up on our earlier work that demonstrated the use of the TOPological Substructural MOlecular DEsign (TOPS-MODE) approach to QSAR modelling for predictions of the carcinogenic potency of nitroso compounds. The data set comprises 56 nitroso compounds which have been bio-assayed in female rats and administered by the oral water route. The QSAR model was able to account for about 81% of the variance in the experimental activity and exhibited good cross-validation statistics. A reasonable interpretation of the TOPS-MODE descriptors was achieved by means of bond contributions, which in turn afforded the recognition of structural alerts (SAs) regarding carcinogenicity. A comparison of the SAs obtained from different data sets showed that experimental factors, such as the sex and the oral administration route, exert a major influence on the carcinogenicity of nitroso compounds. The present and previous QSAR models combined together provide a reliable tool for estimating the carcinogenic potency of yet untested nitroso compounds and they should allow the identification of SAs, which can be used as the basis of prediction systems for the rodent carcinogenicity of these compounds.
Collapse
|
14
|
Phenolic acid derivatives with potential anticancer properties--a structure-activity relationship study. Part 1: methyl, propyl and octyl esters of caffeic and gallic acids. Bioorg Med Chem 2005; 12:3581-9. [PMID: 15186842 DOI: 10.1016/j.bmc.2004.04.026] [Citation(s) in RCA: 207] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2004] [Accepted: 04/19/2004] [Indexed: 11/25/2022]
Abstract
The antiproliferative and cytotoxic properties of polyphenolic acid derivatives, structurally related with the natural models caffeic and gallic acids, have been tested in human cervix adenocarcinoma cells (HeLa). Simultaneous structural information was obtained for these compounds through theoretical ab initio methods. This study was conducted for the following esters: methyl caffeate (MC, 1), propyl caffeate (PC, 2), octyl caffeate (OC, 3), methyl gallate (MG, 4), propyl gallate (PG, 5) and octyl gallate (OG, 6). A significant growth-inhibition effect was assessed for some of these compounds, clearly dependent on their structural characteristics. Marked structure-activity relationships (SARs)--namely the number of hydroxyl ring substituents--were found to rule the biological effect of such systems.
Collapse
|
15
|
Synthesis and QSAR study of the anticancer activity of some novel indane carbocyclic nucleosides. Bioorg Med Chem 2003; 11:4999-5006. [PMID: 14604662 DOI: 10.1016/j.bmc.2003.09.005] [Citation(s) in RCA: 49] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A set of 14 indane carbocyclic nucleosides were synthesized and experimentally assayed for their inhibitory effects in the proliferation of murine leukemia (L1210/0) and human T-lymphocyte (Molt4/C8, CEM/0) cells. The compounds have promising inhibitory activity judging from the IC(50) values obtained for all these cellular lines. Multiple linear regression analysis was then applied to build up consistent QSAR models based on quantum mechanics-derived molecular descriptors. The derived models reproduce well the experimental data of both three cells (r(2) >/=0.90), display a good predictive power and are, above all, easily interpretable. They show that frontier-orbital energies and hydrophobicity are mainly responsible for the activity of the synthesized compounds and also, suggest similar mechanisms of action. The final QSAR-models involve only two descriptors: the lowest unoccupied molecular orbital energy and the solvent accessible-hydrophobic surface area, but describe a sound correlation between predicted and experimental activity data (r(2)=0.931, r(2)=0.936 and r(2)=0.931 for the cells L1210/0, Molt4/C8 and CEM/0, respectively).
Collapse
|