Structure-based discovery of interleukin-33 inhibitors: a pharmacophore modelling, molecular docking, and molecular dynamics simulation approach

Toward Unveiling Putative Binding Sites of Interleukin-33: Insights from Mixed-Solvent Molecular Dynamics Simulations of the Interleukin-1 Family

The interleukin (IL)-1 family is a major proinflammatory cytokine family, ranging from the well-studied IL-1s to the most recently discovered IL-33. As a new focus, IL-33 has attracted extensive research for its crucial immunoregulatory roles, leading to the development of notable monoclonal antibodies as clinical candidates. Efforts to develop small molecules disrupting IL-33/ST2 interaction remain highly desired but encounter challenges due to the shallow and featureless interfaces. The information from relative cytokines has shown that traditional binding site identification methods still struggle in mapping cryptic sites, necessitating dynamic approaches to uncover druggable pockets on IL-33. Here, we employed mixed-solvent molecular dynamics (MixMD) simulations with diverse-property probes to map the hotspots of IL-33 and identify potential binding sites. The protocol was first validated using the known binding sites of two IL-1 family members and then applied to the structure of IL-33. Our simulations revealed several binding sites and proposed side-chain rearrangements essential for the binding of a known inhibitor, aligning well with experimental NMR findings. Further microsecond-time scale simulations of this IL-33-protein complex unveiled distinct binding modes with varying occurrences. These results could facilitate future efforts in developing ligands to target challenging flexible pockets of IL-33 and IL-1 family cytokines in general.

Identification of efflux pump inhibitors for Pseudomonas aeruginosa MexAB-OprM via ligand-based pharmacophores, 2D-QSAR, molecular docking, and molecular dynamics approaches

Efflux pumps have been reported as one of the significant mechanisms by which bacteria evade the effects of multiple antibiotics. The tripartite efflux pump MexAB-OprM in Pseudomonas aeruginosa is one of the most significant multidrug efflux systems due to its broad resistance to antibiotics such as chloramphenicol, fluoroquinolones, lipophilic β-lactam antibiotics, nalidixic acid, novobiocin, rifampicin, and tetracycline. A promising strategy to overcome this resistance mechanism is to combine antibiotics with efflux pump inhibitors (EPIs), which can increase their intracellular concentration to enhance their biological activities. Based on 143 EPIs with chemically diverse skeletons, the 3D pharmacophore and 2D-QSAR modelings were developed and used for the virtual screening on 9.2 million compounds including ZINC15, DrugBank, and Traditional Chinese Medicine databases to identify new EPIs. The molecular docking was also performed to evaluate the binding affinity of potential EPIs to the distal-binding pocket of MexB and resulted in 611 potential EPIs. The structure-activity relationship analyses suggested that nitrogen heterocyclic compounds, piperazine and pyridine scaffolds, and amide derivatives are the most favorable chemically features for MexAB inhibitory activities. The results from molecular dynamics analysis in 100 ns indicated that ZINC009296881 and ZINC009200074 were the most potential MexB inhibitors with strong binding affinity to the distal pocket and MM/GBSA ∆Gbind values of - 38.97 and - 30.19 kcal mol-1, respectively. The predicted pharmacokinetic properties and toxicity of these compounds indicated their potential oral drugs. Multistep virtual screening of EPIs for MexAB-OprM, efflux pump multidrug resistant of P. aeruginosa.

Targeting Olokizumab-Interleukin 6 interaction interface to discover novel IL-6 inhibitors

The IL-6/IL-6R or IL-6/GP130 protein-protein interactions play a significant role in controlling the development of chronic inflammatory diseases, such as rheumatoid arthritis, Castleman disease, psoriasis, and, most recently, COVID-19. Modulating or antagonizing protein-protein interactions of IL6 binding to its receptors by oral drugs promises similar efficacy to biological therapy in patients, namely monoclonal antibodies. In this study, we used a crystal structure of the Fab part of olokizumab in a complex with IL-6 (PDB ID: 4CNI) to uncover starting points for small molecule IL-6 antagonist discovery. Firstly, a structure‑based pharmacophore model of the protein active site cavity was generated to identify possible candidates, followed by virtual screening with a significant database Drugbank. After the docking protocol validation, a virtual screening by molecular docking was carried out and a total of 11 top hits were reported. Detailed analysis of the best scoring molecules was performed with ADME/T analysis and molecular dynamics simulation. Furthermore, the Molecular Mechanics-Generalized Born Surface Area (MM/GBSA) technique has been utilized to evaluate the free binding energy. Based on the finding, one newly obtained compound in this study, namely DB15187, may serve as a lead compound for the discovery of IL-6 inhibitors.Communicated by Ramaswamy H. Sarma.

Discovery of AcrAB-TolC pump inhibitors: Virtual screening and molecular dynamics simulation approach

AcrAB-TolC tripartite efflux pump, which belongs to the RND superfamily, is a main multi-drug efflux system of Escherichia coli (E. coli) because of the broad resistance on various antibiotics. With the discovering of efflux pump inhibitors (EPIs), a combination between these and antibiotics is one of the most promising therapies. Therefore, building a virtual screening model with prediction capacities for the efflux pump inhibitory activities of candidates from DrugBank and ZINC15 dataset, is one of the key goals of this project. Based on the database of 170 diverse chemical structures collected from 28 research journals, two 2D-QSAR models and a 3D-pharmacophore model have been performed. On the AcrB protein (PDB 4DX7), two binding sites have been discovered that match to the hydrophobic trap in the distal pocket and the switch loop in the proximal pocket. After virtual screening processes, twenty candidate AcrAB-TolC inhibitors have been subjected to molecular dynamics simulations, binding free energy calculations and ADMET predictions. The results indicate that three compounds namely DB09233, DB02581, and DB15224 are potential inhibitors with ΔGbind of -42.30 ± 4.58, -40.76 ± 7.30 and -31.06 ± 7.63 kcal.mol-1, respectively.Communicated by Ramaswamy H. Sarma.

Research progress of p38 as a new therapeutic target against morphine tolerance and the current status of therapy of morphine tolerance

With the development of the medical industry, new painkillers continue to appear in people's field of vision, but so far no painkiller can replace morphine. While morphine has a strong analgesic effect, it is also easy to produce pain sensitivity and tolerance. Due to the great inter-individual differences in patient responses, there are few clear instructions on how to optimise morphine administration regimens, which complicates clinicians' treatment strategies and limits the effectiveness of morphine in long-term pain therapy. P38MAPK is a key member of the MAPK family. Across recent years, it has been discovered that p38MAPK rises dramatically in a wide range of morphine tolerance animal models. Morphine tolerance can be reduced or reversed by inhibiting p38MAPK. However, the role and specific mechanism of p38MAPK are not clear. In this review, we synthesise the relevant findings, highlight the function and potential mechanism of p38MAPK in morphine tolerance, as well as the present status and efficacy of morphine tolerance therapy, and underline the future promise of p38MAPK targeted morphine tolerance treatment.

In silico approach to identify novel allosteric intracellular antagonist for blocking the interleukin-8/CXCR2 receptor signaling pathway

The role of interleukin-8 (IL-8) and its receptor CXCR2 in inflammatory responses and tumor development and progression has been well documented. Our study aims to discover novel compounds as CXCR2 antagonists to block the IL-8 signaling pathway using an in silico drug design. Herein, a structure-based pharmacophore model was developed based on the crystal structure of inactive CXCR2 in a complex with an allosteric inhibitor. This model was validated and refined, followed by virtual screening with the ZINC15 database. Subsequent molecular docking allows for predicting the best pose of a ligand inside a receptor binding site. We found that the 35 top-ranked hits exhibited docking scores from -30.81 to -25.28 kJ/mol and better interaction potential comparing the reference inhibitor. Analysis of ADME and toxicity properties revealed the efficacy and safety of the selected seven compounds. To validate the stability of the protein-ligand complex structure MD simulations approach has also been performed and confirmed via the critical parameters. The MD results explained that the CXCR2 receptor bound with two best-proposed molecules, including ZINC77105530 and ZINC93176465, was quite stable states as observed from low RMSD, RMSF, Rg, SASA values, and high occupancy of the interaction types. Finally, our data identified that these compounds play as potential inhibitors of IL-8 signaling pathways with the MM/GBSA binding free energies of -41.77 ± 6.45 kcal/mol and -38.84 ± 6.16 kcal/mol, respectively.Communicated by Ramaswamy H. Sarma.

The Effects of One-Point Mutation on the New Delhi Metallo Beta-Lactamase-1 Resistance toward Carbapenem Antibiotics and β-Lactamase Inhibitors: An In Silico Systematic Approach

Antibiotic resistance has been becoming more and more critical due to bacteria's evolving hydrolysis enzymes. The NDM-1 enzyme could hydrolyze not only carbapenems but also most of β-lactam's antibiotics and inhibitors. In fact, variant strains could impose a high impact on the resistance of bacteria producing NDM-1. Although previous studies showed the effect of some variants toward antibiotics and inhibitors binding, there has been no research systematically evaluating the effects of alternative one-point mutations on the hydrolysis capacity of NDM-1. This study aims to identify which mutants could increase or decrease the effectiveness of antibiotics and β-lactamase inhibitors toward bacteria. Firstly, 35 different variants with a high probability of emergence based on the PAM-1 matrix were constructed and then docked with 5 ligands, namely d-captopril, l-captopril, thiorphan, imipenem, and meropenem. The selected complexes underwent molecular dynamics simulation and free energy binding estimation, with the results showing that the substitutions at residues 122 and 124 most influenced the binding ability of NDM-1 toward inhibitors and antibiotics. The H122R mutant decreases the binding ability between d-captopril and NDM-1 and diminishes the effectiveness of this antibiotic toward Enterobacteriaceae. However, the H122R mutant has a contrary impact on thiorphan, which should be tested in vitro and in vivo in further experiments.

LY3041658/ interleukin-8 complex structure as targets for IL-8 small molecule inhibitors discovery using a combination of in silico methods

Since interleukin-8 (IL-8/CXCL8) and its receptor, CXCR1 and CXCR2, were known in the early 1990s, biological pathways related to these proteins were proven to have high clinical value in cancer and inflammatory/autoimmune conditions treatment. Recently, IL-8 has been identified as biomarker for severe COVID-19 patients and COVID-19 prognosis. Boyles et al. (mAbs 12 (2020), pp. 1831880) have published a high-resolution X-ray crystal structure of the LY3041658 Fab in a complex human CXCL8. They described the ability to bind to IL-8 and the blocking of IL-8/its receptors interaction by the LY3041658 monoclonal antibody. Therefore, the study has been designed to identify potential small molecules inhibiting interleukin-8 by targeting LY3041658/IL-8 complex structure using an in silico approach. A structure‑based pharmacophore and molecular docking models of the protein active site cavity were generated to identify possible candidates, followed by virtual screening with the ZINC database. ADME analysis of hit compounds was also conducted. Molecular dynamics simulations were then performed to survey the behaviour and stability of the ligand-protein complexes. Furthermore, the MM/PBSA technique has been utilized to evaluate the free binding energy. The final data confirmed that one newly obtained compound, ZINC21882765, may serve as the best potential inhibitor for IL-8.

Discovery of benzothiazole-based thiazolidinones as potential anti-inflammatory agents: anti-inflammatory activity, soybean lipoxygenase inhibition effect and molecular docking studies

Despite the greatest achievement in the development of anti-inflammatory agents in the last two decades, the current clinical drugs suffer from a variety of complications in community settings and hospital. There is still an urgent need to design novel molecules with better safety profile and with different molecular targets from those in current clinical use. The aim of this research was to discover a series of benzothiazole-based thiazolidinones with lipoxygenase (LOX) inhibitory activity as a mechanism of anti-inflammatory action. Carrageenan-induced mouse foot paw oedema assay was carried out to determine the anti-inflammatory activity, while LOX inhibition was examined through the conversion of sodium linoleate to 13-hydroperoxylinoleic acid. Molecular docking studies were performed using AutoDock 4.2 software. The anti-inflammatory activity of the title compounds was determined in a range of 18.4%-69.57%, where compound #3 was found to be the most potent (69.57%) and also to be more active than the reference drug indomethacin (47%). Moreover, compound #3 showed the highest LOX inhibitory activity with IC50 of 13 μM being less potent to that of the reference NDGA (IC50 = 1.3 μM). Compound #3 has been identified as lead compound for further modification in an attempt to improve anti-inflammatory and LOX inhibitory activities.

Structure-based 3D-Pharmacophore modeling to discover novel interleukin 6 inhibitors: An in silico screening, molecular dynamics simulations and binding free energy calculations

Interleukin 6 (IL-6) is a cytokine with various biological functions in immune regulation, hematopoiesis, and inflammation. Elevated IL-6 levels have been identified in several severe disorders such as sepsis, acute respiratory distress syndrome (ARDS), and most recently, COVID-19. The biological activity of IL-6 relies on interactions with its specific receptor, IL-6Rα, including the membrane-bound IL-6 receptor (mIL-6R) and the soluble IL-6 receptor (sIL-6R). Thus, inhibition of the interaction between these two proteins would be a potential treatment for IL-6 related diseases. To date, no orally available small-molecule drug has been approved. This study focuses on finding potential small molecules that can inhibit protein-protein interactions between IL-6 and its receptor IL-6Rα using its crystal structure (PDB ID: 5FUC). First, two pharmacophore models were constructed based on the interactions between key residues of IL-6 (Phe74, Phe78, Leu178, Arg179, Arg182) and IL-6Rα (Phe229, Tyr230, Glu277, Glu278, Phe279). A database of approximately 22 million compounds was screened using 3D-pharmacophore models, molecular docking models, and ADMET properties. By analyzing the interactive capability of successfully docked compounds with important amino acids, 12 potential ligands were selected for further analysis via molecular dynamics simulations. Based on the stability of the complexes, the high interactions rate of each ligand with the key residues of IL-6/IL-6Rα, and the low binding free energy calculation, two compounds ZINC83804241 and ZINC02997430, were identified as the most potential IL-6 inhibitor candidates. These results will pave the way for the design and optimization of more specific compounds to combat cytokine storm in severe coronavirus patients.

Caffeic and 3-(3,4-dihydroxyphenyl)glyceric acid derivatives as antimicrobial agent: biological evaluation and molecular docking studies

Herein we report the evaluation of the antimicrobial activity of some previously synthesized 3-(3,4-dihydroxyphenyl)glyceric acid in benzylated and in free 3,4 hydroxy groups in catechol moiety along with some caffeic and 3-(3,4-dihydroxyphenyl)glyceric acid amides using the microdilution method. The evaluation revealed that compounds showed in general moderate to low activity with MIC in range of 0.36-4.5 mg/mL. Compounds were also studied against three resistant bacteria strains MRSA (Methicillin-resistant Staphylococcus aureus), E. coli and P. aeruginosa. Seven out of ten compounds were more potent than reference drugs ampicillin and streptomycin against MRSA, while against another two resistant strains seven compounds showed low activity and the rest were inactive. Antifungal activity of the tested compounds was much better than antibacterial, with MIC in the range of 0.019-3.0 mg/mL. Compounds #7 and 15 showed good activity against all fungi tested, being more potent than ketoconazole and in some case even better than bifonazole used as reference drugs. Docking studies revealed that the most active compound #7 binds to the haem group of the enzyme in the same way as ketoconazole.

Identification of potential interleukin-8 inhibitors acting on the interactive site between chemokine and CXCR2 receptor: A computational approach

Interactions between interleukin (IL)-8 and its receptors, CXCR1, and CXCR2, serve crucial roles in inflammatory conditions and various types of cancers. Inhibition of this signaling pathway has been exploited as a promising strategy in treating these diseases. However, most studies only focused on the design of allosteric antagonists-bound receptors on the intracellular side of IL-8 receptors. Recently, the first cryo-EM structures of IL-8-CXCR2-Gi complexes have been solved, revealing the unique binding and activation modes of the endogenous chemokine IL-8. Hence, we set to identify small molecule inhibitors for IL-8 using critical protein-protein interaction between IL-8 and CXCR2 at the orthosteric binding site. The pharmacophore models and molecular docking screened compounds from DrugBank and NCI databases. The oral bioavailability of the top 23 ligands from the screening was then predicted by the SwissAMDE tool. Molecular dynamics simulation and free binding energy calculation were performed for the best compounds. The result indicated that DB14770, DB12121, and DB03916 could form strong interactions and stable protein-ligand complexes with IL-8. These three candidates are potential IL-8 inhibitors that can be further evaluated by in vitro experiments in the next stage.

Discovery of small molecular inhibitors for interleukin-33/ST2 protein-protein interaction: a virtual screening, molecular dynamics simulations and binding free energy calculations

The interleukin-1 receptor like ST2 has emerged as a potential drug discovery target since it was identified as the receptor of the novel cytokine IL-33, which is involved in many inflammatory and autoimmune diseases. For the treatment of such IL-33-related disorders, efforts have been made to discover molecules that can inhibit the protein–protein interactions (PPIs) between IL-33 and ST2, but to date no drug has been approved. Although several anti-ST2 antibodies have entered clinical trials, the exploration of small molecular inhibitors is highly sought-after because of its advantages in terms of oral bioavailability and manufacturing cost. The aim of this study was to discover ST2 receptor inhibitors based on its PPIs with IL-33 in crystal structure (PDB ID: 4KC3) using virtual screening tools with pharmacophore modeling and molecular docking. From an enormous chemical space ZINC, a potential series of compounds has been discovered with stronger binding affinities than the control compound from a previous study. Among them, four compounds strongly interacted with the key residues of the receptor and had a binding free energy <  − 20 kcal/mol. By intensive calculations using data from molecular dynamics simulations, ZINC59514725 was identified as the most potential candidate for ST2 receptor inhibitor in this study.

In Silico Design, Synthesis, and Biological Evaluation of Anticancer Arylsulfonamide Endowed with Anti-Telomerase Activity

Telomerase, a reverse transcriptase enzyme involved in DNA synthesis, has a tangible role in tumor progression. Several studies have evidenced telomerase as a promising target for developing cancer therapeutics. The main reason is due to the overexpression of telomerase in cancer cells (85–90%) compared with normal cells where it is almost unexpressed. In this paper, we used a structure-based approach to design potential inhibitors of the telomerase active site. The MYSHAPE (Molecular dYnamics SHared PharmacophorE) approach and docking were used to screen an in-house library of 126 arylsulfonamide derivatives. Promising compounds were synthesized using classical and green methods. Compound 2C revealed an interesting IC₅₀ (33 ± 4 µM) against the K-562 cell line compared with the known telomerase inhibitor BIBR1532 IC₅₀ (208 ± 11 µM) with an SI ~10 compared to the BALB/3-T3 cell line. A 100 ns MD simulation of 2C in the telomerase active site evidenced Phe494 as the key residue as well as in BIBR1532. Each moiety of compound 2C was involved in key interactions with some residues of the active site: Arg557, Ile550, and Gly553. Compound 2C, as an arylsulfonamide derivative, is an interesting hit compound that deserves further investigation in terms of optimization of its structure to obtain more active telomerase inhibitors