Mining the PDB for Tractable Cases Where X-ray Crystallography Combined with Fragment Screens Can Be Used to Systematically Design Protein-Protein Inhibitors: Two Test Cases Illustrated by IL1β-IL1R and p38α-TAB1 Complexes

Ligandability Assessment of IL-1β by Integrated Hit Identification Approaches

Human interleukin-1β (IL-1β) is a pro-inflammatory cytokine that plays a critical role in the regulation of the immune response and the development of various inflammatory diseases. In this publication, we disclose our efforts toward the discovery of IL-1β binders that interfere with IL-1β signaling. To this end, several technologies were used in parallel, including fragment-based screening (FBS), DNA-encoded library (DEL) technology, peptide discovery platform (PDP), and virtual screening. The utilization of distinct technologies resulted in the identification of new chemical entities exploiting three different sites on IL-1β, all of them also inhibiting the interaction with the IL-1R1 receptor. Moreover, we identified lysine 103 of IL-1β as a target residue suitable for the development of covalent, low-molecular-weight IL-1β antagonists.

An integrated strategy of UPLC-Q-TOF-MS analysis, network pharmacology, and molecular docking to explore the chemical constituents and mechanism of Zixue Powder against febrile seizures

Febrile seizures (FS) are the most common type of seizures for children. As a commonly used representative cold formula for resuscitation, Zixue Powder (ZP) has shown great efficacy for the treatment of FS in clinic, while its active ingredients and underlying mechanism remain largely unclear. This study aimed to preliminarily elucidate the material basis of ZP and the potential mechanism for the treatment of FS through ultra-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry (UPLC-Q-TOF-MS), network pharmacology, and molecular docking. UPLC-Q-TOF-MS was firstly applied to characterize the ingredients in ZP, followed by network pharmacology to explore the potential bioactive ingredients and pathways of ZP against FS. Furthermore, molecular docking technique was employed to verify the binding affinity between the screened active ingredients and targets. As a result, 75 ingredients were identified, containing flavonoids, chromogenic ketones, triterpenes and their saponins, organic acids, etc. Through the current study, we focused on 13 potential active ingredients and 14 key potential anti-FS targets of ZP, such as IL6, STAT3, TNF, and MMP9. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that inflammatory response, EGFR tyrosine kinase inhibitor resistance, AGE-RAGE signaling pathway in diabetic complications, and neuroactive ligand-receptor interaction were the main anti-FS signaling pathways. Licochalcones A and B, 26-deoxycimicifugoside, and hederagenin were screened as the main potential active ingredients by molecular docking. In conclusion, this study provides an effective in-depth investigation of the chemical composition, potential bioactive components, and possible anti-FS mechanism of ZP, which lays the foundation for pharmacodynamic studies and clinical applications of ZP.

Isolation, Identification, Anti-Inflammatory, and In Silico Analysis of New Lignans from the Resin of Ferula sinkiangensis

Ferula sinkiangensis K. M. Shen (Apiaceae) is distributed in arid desert areas of Xinjiang, and its resin is a traditional Chinese medicine to treat gastrointestinal digestive diseases. To explore bioactive components from F. sinkiangensis, three new lignans and thirteen known components were isolated. The structural elucidation of the components was established utilizing spectroscopic analyses together with ECD calculations. Griess reaction results indicated new compounds 1 and 2 significantly decreased NO production in LPS-stimulated RAW 264.7 macrophages, and ELISA results indicated that they effectively attenuated LPS-induced inflammation by inhibiting TNF-α, IL-1β, and IL-6 expressions. The in silico approach confirmed that compound 1 docked into the receptors with strong binding energies of -5.84~-10.79 kcal/mol. In addition, compound 6 inhibited the proliferation of AGS gastric cancer cells with IC50 values of 15.2 μM by suppressing the cell migration and invasion. This study disclosed that F. sinkiangensis might be a promising potential resource for bioactive components.

Discovery of a selective and biologically active low-molecular weight antagonist of human interleukin-1β

Human interleukin-1β (hIL-1β) is a pro-inflammatory cytokine involved in many diseases. While hIL-1β directed antibodies have shown clinical benefit, an orally available low-molecular weight antagonist is still elusive, limiting the applications of hIL-1β-directed therapies. Here we describe the discovery of a low-molecular weight hIL-1β antagonist that blocks the interaction with the IL-1R1 receptor. Starting from a low affinity fragment-based screening hit 1, structure-based optimization resulted in a compound (S)-2 that binds and antagonizes hIL-1β with single-digit micromolar activity in biophysical, biochemical, and cellular assays. X-ray analysis reveals an allosteric mode of action that involves a hitherto unknown binding site in hIL-1β encompassing two loops involved in hIL-1R1/hIL-1β interactions. We show that residues of this binding site are part of a conformationally excited state of the mature cytokine. The compound antagonizes hIL-1β function in cells, including primary human fibroblasts, demonstrating the relevance of this discovery for future development of hIL-1β directed therapeutics.

Combining structural and coevolution information to unveil allosteric sites

Understanding allosteric regulation in biomolecules is of great interest to pharmaceutical research and computational methods emerged during the last decades to characterize allosteric coupling. However, the prediction of allosteric sites in a protein structure remains a challenging task. Here, we integrate local binding site information, coevolutionary information, and information on dynamic allostery into a structure-based three-parameter model to identify potentially hidden allosteric sites in ensembles of protein structures with orthosteric ligands. When tested on five allosteric proteins (LFA-1, p38-α, GR, MAT2A, and BCKDK), the model successfully ranked all known allosteric pockets in the top three positions. Finally, we identified a novel druggable site in MAT2A confirmed by X-ray crystallography and SPR and a hitherto unknown druggable allosteric site in BCKDK validated by biochemical and X-ray crystallography analyses. Our model can be applied in drug discovery to identify allosteric pockets.

Characterization of p38α autophosphorylation inhibitors that target the non-canonical activation pathway

p38α is a versatile protein kinase that can control numerous processes and plays important roles in the cellular responses to stress. Dysregulation of p38α signaling has been linked to several diseases including inflammation, immune disorders and cancer, suggesting that targeting p38α could be therapeutically beneficial. Over the last two decades, numerous p38α inhibitors have been developed, which showed promising effects in pre-clinical studies but results from clinical trials have been disappointing, fueling the interest in the generation of alternative mechanisms of p38α modulation. Here, we report the in silico identification of compounds that we refer to as non-canonical p38α inhibitors (NC-p38i). By combining biochemical and structural analyses, we show that NC-p38i efficiently inhibit p38α autophosphorylation but weakly affect the activity of the canonical pathway. Our results demonstrate how the structural plasticity of p38α can be leveraged to develop therapeutic opportunities targeting a subset of the functions regulated by this pathway.

Discovery of Chalcone-Based Hybrid Structures as High Affinity and Site-Specific Inhibitors against SARS-CoV-2: A Comprehensive Structural Analysis Based on Various Host-Based and Viral Targets

Previous studies indicated that natural-based chalcones have significant inhibitory effects on the coronavirus enzymes 3CLpro and PLpro as well as modulation of some host-based antiviral targets (HBATs). In this study, a comprehensive computational and structural study was performed to investigate the affinity of our compound library consisting of 757 chalcone-based structures (CHA-1 to CHA-757) for inhibiting the 3CLpro and PLpro enzymes and against twelve selected host-based targets. Our results indicated that CHA-12 (VUF 4819) is the most potent and multi-target inhibitor in our chemical library over all viral and host-based targets. Correspondingly, CHA-384 and its congeners containing ureide moieties were found to be potent and selective 3CLpro inhibitors, and benzotriazole moiety in CHA-37 was found to be a main fragment for inhibiting the 3CLpro and PLpro. Surprisingly, our results indicate that the ureide and sulfonamide moieties are integral fragments for the optimum 3CLpro inhibition while occupying the S1 and S3 subsites, which is fully consistent with recent reports on the site-specific 3CLpro inhibitors. Finding the multi-target inhibitor CHA-12, previously reported as an LTD4 antagonist for the treatment of inflammatory pulmonary diseases, prompted us to suggest it as a concomitant agent for relieving respiratory symptoms and suppressing COVID-19 infection.

'Chemistry at the speed of sound': automated 1536-well nanoscale synthesis of 16 scaffolds in parallel

Screening of large and diverse libraries is the 'bread and butter' in the first phase of the discovery of novel drugs. However, maintenance and periodic renewal of high-quality large compound collections pose considerable logistic, environmental and monetary problems. Here, we exercise an alternative, the 'on-the-fly' synthesis of large and diverse libraries on a nanoscale in a highly automated fashion. For the first time, we show the feasibility of the synthesis of a large library based on 16 different chemistries in parallel on several 384-well plates using the acoustic dispensing ejection (ADE) technology platform. In contrast to combinatorial chemistry, we produced 16 scaffolds at the same time and in a sparse matrix fashion, and each compound was produced by a random combination of diverse large building blocks. The synthesis, analytics, resynthesis of selected compounds, and chemoinformatic analysis of the library are described. The advantages of the herein described automated nanoscale synthesis approach include great library diversity, absence of library storage logistics, superior economics, speed of synthesis by automation, increased safety, and hence sustainable chemistry.

Molecular Mechanics Simulations and Experimental Investigation of the Effect of Tadalafil on Various Inflammatory Pain Mediators

PURPOSE

Tadalafil's exact analgesic mechanism is still unclear. The current study aimed to elucidate this mechanism in an inflammatory pain model.

METHODS

Computer-assisted simulation docking experiments were carried out to assess the binding of tadalafil to different ligands. The anti-inflammatory and analgesic effects of tadalafil were evaluated using formalin-induced paw edema and a von Frey filament test, respectively. The plantar paw of the mice was then dissected to quantify iNOS, nNOS, COX-2, TNFα, IL1, and IL10 gene expression levels using a real-time polymerase chain reaction. iNOS, TNFα, and COX-2 inhibition was reassessed in vitro using the ELISA technique. One-way analysis of variance followed by post hoc Tukey test or t-test was used to compare the means.

RESULTS

Docking analysis showed a superior binding score of tadalafil to COX-2, iNOS, IL-1, and TNF-α compared to that of indomethacin and morphine and a similar binding score to nNOS and IL-10 relative to that of indomethacin. In the in vivo study, tadalafil, after an hour of formalin administration, inhibited significantly paw edema, similar to indomethacin. Furthermore, it significantly increased the withdrawal force in the von Frey filament test as compared to the negative control, which was similar to the effect observed with indomethacin and morphine. The RT-PCR revealed that tadalafil reduced significantly the iNOS, COX-2, and TNF-α gene expressions but had no effect on nNOS, IL 1, and IL10. In vitro ELISA tests confirmed the inhibition of iNOS, COX-2, and TNF-α.

CONCLUSION

Tadalafil probably exerts its analgesic effect through the simultaneous inhibition of iNOS, COX-2, and TNF-α, which is not the case with other nonsteroidal anti-inflammatory drugs. Nevertheless, further studies are required to confirm its mechanism.

TAB1 binding induced p38α conformation change: an accelerated molecular dynamics simulation study

p38α mitogen-activated protein kinase (MAPK) undergoes autophosphorylation induced by the binding of TGFβ-activated kinase 1 binding protein 1 (TAB1) in myocardial ischemia. Investigation of the conformational transformations in p38α triggered by TAB1 binding is motivated by the need to find selective p38α activation inhibitors to treat myocardial ischemia. Herein, the conformational transformations of p38α were studied via all-atom accelerated molecular dynamics simulations and principal component analysis. With the binding of TAB1, the conformational changes of p38α auto-activation were characterized by the movement of the activation loop (A-loop) away from the αG helix toward the αF, αE helixes and L16-loop. In addition, a diverse intermediate state with an extensional and phosphorylated A-loop different from the transition intermediate state was explored. The conformational changes, including the A-loop alpha-structure breaking and the stronger hydrogen bond network formation, are accompanied by the extension of the A-loop and more intramolecular interactions in p38α. TAB1 correlates with other regions of p38α that are distal from the TAB1-binding site, including the A-loop, αC helix, and L16-loop, which regulates the intramolecular correlation of p38α. And, the phosphorylation further enhances the correlations between the A-loop and the other regions of p38α. The correlation results imply the regulation process of p38α conformational transformations. These findings will improve our understanding of the autophosphorylation of kinase and facilitate the development of selective inhibitors for the treatment of ischemic injury.

Jiedu Huoxue Decoction for Cytokine Storm and Thrombosis in Severe COVID-19: A Combined Bioinformatics and Computational Chemistry Approach

Jiedu Huoxue Decoction (JHD), a recommended traditional prescription for patients with severe COVID-19, has appeared in the treatment protocols in China. Based on bioinformatics and computational chemistry methods, including molecular docking, molecular dynamics (MD) simulation, and Molecular Mechanics Generalized Born Surface Area (MM/GBSA) calculation, we aimed to reveal the mechanism of JHD in treating severe COVID-19. The compounds in JHD were obtained and screened on TCMSP, SwissADME, and ADMETLab platforms. The compound targets were obtained from TCMSP and STITCH, while COVID-19 targets were obtained from Genecards and NCBI. The protein-protein interaction network was constructed by using STRING. Gene Ontology (GO) and KEGG enrichment were performed with ClueGO and R language. AutoDock vina was employed for molecular docking. 100 ns MD simulation of the optimal docking complex was carried out with AmberTools 20. A total of 84 compounds and 29 potential targets of JHD for COVID-19 were collected. The key phytochemicals included quercetin, luteolin, β-sitosterol, puerarin, stigmasterol, kaempferol, and wogonin, which could regulate the immune system. The hub genes included IL6, IL10, VEGFA, IL1B, CCL2, HMOX1, DPP4, and ACE2. ACE2 and DPP4 were related to SARS-CoV-2 entering cells. GO and KEGG analysis showed that JHD could intervene in cytokine storm and endothelial proliferation and migration related to thrombosis. The molecular docking, 100 ns MD simulation, and MM/GBSA calculation confirmed that targets enriched in the COVID-19 pathway had high affinities with related compounds, and the conformations of the puerarin-ACE2, quercetin-EGFR, luteolin-EGFR, and quercetin-IL1B complexes were stable. In a word, JHD could treat COVID-19 by intervening in cytokine storm, thrombosis, and the entry of SARS-CoV-2, while regulating the immune system. These mechanisms were consistent with JHD's therapeutic concept of “detoxification” and “promoting blood circulation and removing blood stasis” in treating COVID-19. The research provides a theoretical basis for the development and application of JHD.

Computational Cosolvent Mapping Analysis Leads to Identify Salicylic Acid Analogs as Weak Inhibitors of ST2 and IL33 Binding

Cytokine signaling initiated by the binding of the cytokine receptors to cytokines plays important roles in immune regulation and diseases. Structurally, cytokine receptors interact with cytokines via an extensive, rugged interface that represents a challenge in inhibitor development. Our computational analysis has previously indicated that butyric acid, mimicking acidic residues, preferentially binds to sites in ST2 (Stimulation-2) that interact with acidic residues of IL33, the endogenous cytokine for ST2. To investigate if a charged group in small molecules facilitates ligand binding to ST2, we developed a biochemical homogeneous time resolved fluorescence assay to determine the inhibition of ST2/IL33 binding by five molecules containing an aromatic ring and a charged group. Three molecules, including niacin, salicylic acid, and benzamidine, exhibit inhibition activities at millimolar concentrations. We further employed the computational cosolvent mapping analysis to identify a shared mode of interaction between niacin, salicylic acid, and ST2. The mode of interaction was further confirmed by four analogous compounds that exhibited similar or improved activities. Our study provided the evidence of inhibition of ST2 and IL33 binding by salicylic acid and analogs. The results suggest that biological activity of salicylic acid may be partly mediated through modulating extracellular cytokine receptors and cytokine interaction.

Combining High-Throughput Synthesis and High-Throughput Protein Crystallography for Accelerated Hit Identification

Protein crystallography (PX) is widely used to drive advanced stages of drug optimization or to discover medicinal chemistry starting points by fragment soaking. However, recent progress in PX could allow for a more integrated role into early drug discovery. Here, we demonstrate for the first time the interplay of high throughput synthesis and high throughput PX. We describe a practical multicomponent reaction approach to acrylamides and -esters from diverse building blocks suitable for mmol scale synthesis on 96-well format and on a high-throughput nanoscale format in a highly automated fashion. High-throughput PX of our libraries efficiently yielded potent covalent inhibitors of the main protease of the COVID-19 causing agent, SARS-CoV-2. Our results demonstrate, that the marriage of in situ HT synthesis of (covalent) libraires and HT PX has the potential to accelerate hit finding and to provide meaningful strategies for medicinal chemistry projects.

Binding mechanism of inhibitors to p38α MAP kinase deciphered by using multiple replica Gaussian accelerated molecular dynamics and calculations of binding free energies

The p38α MAP Kinase has been an important target of drug design for treatment of inflammatory diseases and cancers. This work applies multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations and the molecular mechanics generalized Born surface area (MM-GBSA) method to probe the binding mechanism of inhibitors L51, R24 and 1AU to p38α. Dynamics analyses show that inhibitor bindings exert significant effect on conformational changes of the active helix α2 and the conserved DFG loop. The rank of binding free energies calculated with MM-GBSA not only agrees well with that determined by the experimental IC50 values but also suggests that mutual compensation between the enthalpy and entropy changes can improve binding of inhibitors to p38α. The analyses of free energy landscapes indicate that the L51, R24 and 1AU bound p38α display a DFG-out conformation. The residue-based free energy decomposition method is used to evaluate contributions of separate residues to the inhibitor-p38α binding and the results imply that residues V30, V38, L74, L75, I84, T106, H107, L108, M109, L167, F169 and D168 can be utilized as efficient targets of potent inhibitors toward p38α.

Atypical p38 Signaling, Activation, and Implications for Disease

The mitogen-activated protein kinase (MAPK) p38 is an essential family of kinases, regulating responses to environmental stress and inflammation. There is an ever-increasing plethora of physiological and pathophysiological conditions attributed to p38 activity, ranging from cell division and embryonic development to the control of a multitude of diseases including retinal, cardiovascular, and neurodegenerative diseases, diabetes, and cancer. Despite the decades of intense investigation, a viable therapeutic approach to disrupt p38 signaling remains elusive. A growing body of evidence supports the pathological significance of an understudied atypical p38 signaling pathway. Atypical p38 signaling is driven by a direct interaction between the adaptor protein TAB1 and p38α, driving p38 autophosphorylation independent from the classical MKK3 and MKK6 pathways. Unlike the classical MKK3/6 signaling pathway, atypical signaling is selective for just p38α, and at present has only been characterized during pathophysiological stimulation. Recent studies have linked atypical signaling to dermal and vascular inflammation, myocardial ischemia, cancer metastasis, diabetes, complications during pregnancy, and bacterial and viral infections. Additional studies are required to fully understand how, when, where, and why atypical p38 signaling is induced. Furthermore, the development of selective TAB1-p38 inhibitors represents an exciting new opportunity to selectively inhibit pathological p38 signaling in a wide array of diseases.

Diversity and versatility of p38 kinase signalling in health and disease

The ability of cells to deal with different types of stressful situations in a precise and coordinated manner is key for survival and involves various signalling networks. Over the past 25 years, p38 kinases — in particular, p38α — have been implicated in the cellular response to stress at many levels. These span from environmental and intracellular stresses, such as hyperosmolarity, oxidative stress or DNA damage, to physiological situations that involve important cellular changes such as differentiation. Given that p38α controls a plethora of functions, dysregulation of this pathway has been linked to diseases such as inflammation, immune disorders or cancer, suggesting the possibility that targeting p38α could be of therapeutic interest. In this Review, we discuss the organization of this signalling pathway focusing on the diversity of p38α substrates, their mechanisms and their links to particular cellular functions. We then address how the different cellular responses can be generated depending on the signal received and the cell type, and highlight the roles of this kinase in human physiology and in pathological contexts.

p38α — the best-characterized member of the p38 kinase family — is a key mediator of cellular stress responses. p38α is activated by a plethora of signals and functions through a multitude of substrates to regulate different cellular behaviours. Understanding context-dependent p38α signalling provides important insights into p38α roles in physiology and pathology.

A New Crystal Form of the SARS-CoV-2 Receptor Binding Domain: CR3022 Complex-An Ideal Target for In-Crystal Fragment Screening of the ACE2 Binding Site Surface

In-crystal fragment screening is a powerful tool to chemically probe the surfaces used by proteins to interact, and identify the chemical space worth exploring to design protein-protein inhibitors. A crucial prerequisite is the identification of a crystal form where the target area is exposed and accessible to be probed by fragments. Here we report a crystal form of the SARS-CoV-2 Receptor Binding Domain in complex with the CR3022 antibody where the ACE2 binding site on the Receptor Binding Domain is exposed and accessible. This crystal form of the complex is a valuable tool to develop antiviral molecules that could act by blocking the virus entry in cells.