Targeting cryptic-orthosteric site of PD-L1 for inhibitor identification using structure-guided approach

Progress in small-molecule inhibitors targeting PD-L1

PD-L1 is a transmembrane protein overexpressed by tumor cells. It binds to PD-1 on the surface of T-cells, suppresses T-cell activity and hinders the immune response against cancer. Clinically, several monoclonal antibodies targeting PD-1/PD-L1 have achieved significant success in cancer immunotherapy. Nevertheless, their disadvantages, such as unchecked immune responses, high cost and long half-life, stimulated pharmacologists to develop small-molecule inhibitors targeting PD-1/PD-L1. After a batch of excellent inhibitors with a biphenyl core structure were firstly reported by BMS, more and more researchers focused on small-molecule inhibitors targeting PD-L1 rather than PD-1. Numerous small-molecule inhibitors were extensively designed and synthesized in the past few years. In this paper, the structural characteristics of PD-L1 and complexes of PD-L1 with its inhibitors are elaborated and small molecule inhibitors developed in the last decade are summarized as well. This paper aims to provide insights into further designing and synthesis of small molecule inhibitors targeting PD-L1.

Hotspot site microenvironment in the deubiquitinase OTUB1 drives its stability and aggregation

Lewy bodies (LB) are aberrant protein accumulations observed in the brain cells of individuals affected by Parkinson's disease (PD). A comprehensive analysis of LB proteome identified over a hundred proteins, many co-enriched with α-synuclein, a major constituent of LB. Within this context, OTUB1, a deubiquitinase detected in LB, exhibits amyloidogenic properties, yet the mechanisms underlying its aggregation remain elusive. In this study, we identify two critical sites in OTUB1-namely, positions 133 and 173-that significantly impact its amyloid aggregation. Substituting alanine at position 133 and lysine at position 173 enhances both thermodynamic and kinetic stability, effectively preventing amyloid aggregation. Remarkably, lysine at position 173 demonstrates the highest stability without compromising enzymatic activity. The increased stability and inhibition of amyloid aggregation are attributed mainly to the changes in the specific microenvironment at the hotspot. In our exploration of the in-vivo co-occurrence of α-synuclein and OTUB1 in LB, we observed a synergistic modulation of each other's aggregation. Collectively, our study unveils the molecular determinants influencing OTUB1 aggregation, shedding light on the role of specific residues in modulating aggregation kinetics and structural transition. These findings contribute valuable insights into the complex interplay of amino acid properties and protein aggregation, with potential implications for understanding broader aspects of protein folding and aggregation phenomena.

Elevated glycosylation of CD36 in platelets is a risk factor for oxLDL-mediated platelet activation in type 2 diabetes

Platelet activation and related cardiovascular complications are the hallmarks of type 2 diabetes (T2D). We investigated the mechanism of platelet activation in T2D using MS-based identification of differentially expressed platelet proteins with a focus on glycosylated forms. Glycosylation is considered one of the common post-translational modifications in T2D, and N/O-linked glycosylation of glycoproteins (GPs)/integrins is known to play crucial roles in platelet activation. Our platelet proteome data revealed elevated levels of GPs GPIbα, GPIIbIIIa, GPIV (CD36), GPV and integrins in T2D patients. T2D platelets had elevated N-linked glycosylation of CD36 at asparagine (Asn)408,417 . Enrichment analysis revealed a close association of glycosylated CD36 with thrombospondin-1, fibrinogen and SERPINA1 in T2D platelets. The glycosylation of CD36 has previously been reported to increase cellular uptake of long-chain fatty acids. Our in silico molecular docking data also showed a favorable binding of cholesterol with glycosylated Asn417 CD36 compared to the non-glycosylated form. We further investigated the CD36:LDL cholesterol axis in T2D. Elevated levels of oxidized-low density lipoprotein (oxLDL) were found to cause significant platelet activation via CD36-mediated stimulation of Lyn-JNK signaling. Sulfo-N-succinimidyl oleate, an inhibitor of CD36, effectively inhibited oxLDL-mediated platelet activation and adhesion in vitro. Our study suggests increased glycosylation of CD36 in T2D platelets as a potential route for oxLDL-mediated platelet activation. The oxLDL:CD36 axis may thus be exploited as a prospective target to develop therapeutics against thrombosis in T2D.

Harnessing the druggability at orthosteric and allosteric sites of PD-1 for small molecule discovery by an integrated in silico pipeline

The PD-1/PD-L1 interaction is a promising target for small molecule inhibitors in cancer immunotherapy, but targeting this interface has been challenging. While efforts have been made to identify compounds that target the orthosteric sites, no reports have explored the potential of small molecules to target the allosteric region of PD-1. Therefore, our study aims to establish a pipeline to identify small molecules that can effectively bind to either the orthosteric or allosteric pockets of PD-1. We categorized the PD-1 interface into two hot-spot zones (P-and N-zones) based on extensive analysis of its structural, dynamical, and energetic properties. These zones correspond to the orthosteric and allosteric PPI sites, respectively, targeted by monoclonal antibodies. We used a guided virtual screening workflow to identify hits from ∼7 million compounds library, which were then clustered based on structural similarity and assessed by interaction fingerprinting. The selective and diverse chemical representatives were subjected to MD simulations and binding energetics calculations to filter out false positives and identify actual binders. Binding poses metadynamics calculations confirmed the stability of the final hits in the pocket. This study emphasizes the need for an integrated pipeline that uses molecular dynamics simulations and binding energetics to identify potential binders for the dynamic PD-1/PD-L1 interface, due to the lack of small molecule co-crystals. Only a few potential binders were discovered from a large pool of molecules targeting both the allosteric and orthosteric zones. Our results suggest that the allosteric site has more potential than the orthosteric site for inhibitor design. The identified "computational hits" hold potential as starting points for in vitro evaluations followed by hit-to-lead optimization. Overall, this study represents an effort to establish a computational pipeline for exploring and enriching both the allosteric and orthosteric sites of PPI interfaces, "a tough but indispensable nut to crack".

Understanding the role of conserved proline and serine residues in the SARS-CoV-2 spike cleavage sites in the virus entry, fusion, and infectivity

The spike (S) glycoprotein of the SARS-CoV-2 virus binds to the host cell receptor and promotes the virus's entry into the target host cell. This interaction is primed by host cell proteases like furin and TMPRSS2, which act at the S1/S2 and S2´ cleavage sites, respectively. Both cleavage sites have serine or proline residues flanking either the single or polybasic region and were found to be conserved in coronaviruses. Unravelling the effects of these conserved residues on the virus entry and infectivity might facilitate the development of novel therapeutics. Here, we have investigated the role of the conserved serine and proline residues in the SARS-CoV-2 spike mediated entry, fusogenicity, and viral infectivity by using the HIV-1/spike-based pseudovirus system. A conserved serine residue mutation to alanine (S2´S-A) at the S2´ cleavage site resulted in the complete loss of spike cleavage. Exogenous treatment with trypsin or overexpression of TMPRSS2 protease could not rescue the loss of spike cleavage and biological activity. The S2´S-A mutant showed no significant responses against E-64d, TMPRSS2 or other relevant inhibitors. Taken together, serine at the S2´ site in the spike protein was indispensable for spike protein cleavage and virus infectivity. Thus, novel interventions targeting the conserved serine at the S2´ cleavage site should be explored to reduce severe disease caused by SARS-CoV-2-and novel emerging variants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03749-y.

A multitier virtual screening of antagonists targeting PD-1/PD-L1 interface for the management of triple-negative breast cancer

Immunotherapies are promising therapeutic options for the management of triple-negative breast cancer because of its high mutation rate and genomic instability. Of note, the blockade of the immune checkpoint protein PD-1 and its ligand PD-L1 has been proven to be an efficient and potent strategy to combat triple-negative breast cancer. To date, various anti-PD-1/anti-PD-L1 antibodies have been approved. However, the intrinsic constraints of these therapeutic antibodies significantly limit their application, making small molecules a potentially significant option for PD-1/PD-L1 inhibition. In light of this, the current study aims to use a high-throughput virtual screening technique to identify potential repurposed candidates as PD-L1 inhibitors. Thus, the present study explored binding efficiency of 2509 FDA-approved compounds retrieved from the drug bank database against PD-L1 protein. The binding affinity of the compounds was determined using the glide XP docking programme. Furthermore, prime-MM/GBSA, DFT calculations, and RF score were used to precisely re-score the binding free energy of the docked complexes. In addition, the ADME and toxicity profiles for the lead compounds were also examined to address PK/PD characteristics. Altogether, the screening process identified three molecules, namely DB01238, DB06016 and DB01167 as potential therapeutics for the PD-L1 protein. To conclude, a molecular dynamic simulation of 100 ns was run to characterise the stability and inhibitory action of the three lead compounds. The results from the simulation study confirm the robust structural and thermodynamic stability of DB01238 than other investigated molecules. Thus, our findings hypothesize that DB01238 could serve as potential PD-L1 inhibitor in the near future for triple-negative breast cancer patients.

Non-oncology drug (meticrane) shows anti-cancer ability in synergy with epigenetic inhibitors and appears to be involved passively in targeting cancer cells

Emerging evidence suggests that chemotherapeutic agents and targeted anticancer drugs have serious side effects on the healthy cells/tissues of the patient. To overcome this, the use of non-oncology drugs as potential cancer therapies has been gaining momentum. Herein, we investigated one non-oncology drug named meticrane (a thiazide diuretic used to treat essential hypertension), which has been reported to indescribably improve the therapeutic efficacy of anti-CTLA4 in mice with AB1 HA tumors. In our hypothesis-driven study, we tested anti-cancer potential meticrane in hematological malignance (leukemia and multiple myeloma) and liver cancer cell lines. Our analysis showed that: 1) Meticrane induced alteration in the cell viability and proliferation in leukemia cells (Jurkat and K562 cells) and liver cancer (SK-hep-1), however, no evidence of apoptosis was detectable. 2) Meticrane showed additive/synergistic effects with epigenetic inhibitors (DNMT1/5AC, HDACs/CUDC-101 and HDAC6/ACY1215). 3) A genome-wide transcriptional analysis showed that meticrane treatment induces changes in the expression of genes associated with non-cancer associated pathways. Of importance, differentially expressed genes showed favorable correlation with the survival-related genes in the cancer genome. 4) We also performed molecular docking analysis and found considerable binding affinity scores of meticrane against PD-L1, TIM-3, CD73, and HDACs. Additionally, we tested its suitability for immunotherapy against cancers, but meticrane showed no response to the cytotoxicity of cytokine-induced killer (CIK) cells. To our knowledge, our study is the first attempt to identify and experimentally confirm the anti-cancer potential of meticrane, being also the first to test the suitability of any non-oncology drug in CIK cell therapy. Beyond that, we have expressed some concerns confronted during testing meticrane that also apply to other non-oncology drugs when considered for future clinical or preclinical purposes. Taken together, meticrane is involved in some anticancer pathways that are passively targeting cancer cells and may be considered as compatible with epigenetic inhibitors.

Intradermal Immunization of Soluble Influenza HA Derived from a Lethal Virus Induces High Magnitude and Breadth of Antibody Responses and Provides Complete Protection In Vivo

Immunogens mimicking the native-like structure of surface-exposed viral antigens are considered promising vaccine candidates. Influenza viruses are important zoonotic respiratory viruses with high pandemic potential. Recombinant soluble hemagglutinin (HA) glycoprotein-based protein subunit vaccines against Influenza have been shown to induce protective efficacy when administered intramuscularly. Here, we have expressed a recombinant soluble trimeric HA protein in Expi 293F cells and purified the protein derived from the Inf A/Guangdong-Maonan/ SWL1536/2019 virus which was found to be highly virulent in the mouse. The trimeric HA protein was found to be in the oligomeric state, highly stable, and the efficacy study in the BALB/c mouse challenge model through intradermal immunization with the prime-boost regimen conferred complete protection against a high lethal dose of homologous and mouse-adapted InfA/PR8 virus challenge. Furthermore, the immunogen induced high hemagglutinin inhibition (HI) titers and showed cross-protection against other Inf A and Inf B subtypes. The results are promising and warrant trimeric HA as a suitable vaccine candidate.

Computational Approaches Drive Developments in Immune-Oncology Therapies for PD-1/PD-L1 Immune Checkpoint Inhibitors

Computational approaches in immune-oncology therapies focus on using data-driven methods to identify potential immune targets and develop novel drug candidates. In particular, the search for PD-1/PD-L1 immune checkpoint inhibitors (ICIs) has enlivened the field, leveraging the use of cheminformatics and bioinformatics tools to analyze large datasets of molecules, gene expression and protein-protein interactions. Up to now, there is still an unmet clinical need for improved ICIs and reliable predictive biomarkers. In this review, we highlight the computational methodologies applied to discovering and developing PD-1/PD-L1 ICIs for improved cancer immunotherapies with a greater focus in the last five years. The use of computer-aided drug design structure- and ligand-based virtual screening processes, molecular docking, homology modeling and molecular dynamics simulations methodologies essential for successful drug discovery campaigns focusing on antibodies, peptides or small-molecule ICIs are addressed. A list of recent databases and web tools used in the context of cancer and immunotherapy has been compilated and made available, namely regarding a general scope, cancer and immunology. In summary, computational approaches have become valuable tools for discovering and developing ICIs. Despite significant progress, there is still a need for improved ICIs and biomarkers, and recent databases and web tools have been compiled to aid in this pursuit.

Discovery of non-nucleoside oxindole derivatives as potent inhibitors against dengue RNA-dependent RNA polymerase

A series of thiazole linked Oxindole-5-Sulfonamide (OSA) derivatives were designed as inhibitors of RNA-dependent RNA polymerase (RdRp) activity of Dengue virus. These were synthesized and then evaluated for their efficacy in ex-vivo virus replication assay using human cell lines. Among 20 primary compounds in the series, OSA-15 was identified as a hit. A series of analogues were synthesized by replacing the difluoro benzyl group of OSA-15 with different substituted benzyl groups. The efficacy of OSA-15derivatives was less than that of the parent compound, except OSA-15-17, which has shown improved efficacy than OSA-15. The further optimization was carried out by adding dimethyl (DM) groups to both the sulfonamide and oxindole NH's to produce OSA-15-DM and OSA-15-17-DM. These two compounds were showing no detectable cytotoxicity and the latter was more efficacious. Further, both these compounds were tested for inhibition in all the serotypes of the Dengue virus using an ex-vivo assay. The EC₅₀ of OSA-15-17-DM was observed in a low micromolar range between 2.5 and 5.0 µg/ml. Computation docking and molecular dynamics simulation studies confirmed the binding of identified hits to DENV RdRp. OSA15-17-DM blocks the RNA entrance and elongation site for their biological activity with high binding affinity. Overall, the identified oxindole derivatives are novel compounds that can inhibit Dengue replication, working as non-nucleoside inhibitors (NNI) to explore as anti-viral RdRp activity.

Deciphering the Structural Determinants Critical in Attaining the FXR Partial Agonism

Elucidation of structural determinants is pivotal for structure-based drug discovery. The Farnesoid X receptor (FXR) is a proven target for NASH; however, its full agonism causes certain clinical complications. Therefore, partial agonism (PA) appears as a viable alternative for improved therapeutics. Since the agonist and PA both share the same binding site, i.e., ligand-binding pocket (LBP), which is highly dynamic and has synergy with the substrate binding site, the selective designing of PA is challenging. The identification of structural and conformational determinants is critical for PA compared with an agonist. Furthermore, the mechanism by which PA modulates the structural dynamics of FXR at the residue level, a prerequisite for PA designing, is still elusive. Here, by using ∼4.5 μs of MD simulations and residue-wise communication network analysis, we identified the structural regions which are flexible with PA but frozen with an agonist. Also, the network analysis identified the considerable changes between an agonist and PA in biologically essential zones of FXR such as helix H10/H11 and loop L:H11/H12, which lead to the modulation of synergy between LBP and the substrate binding site. Furthermore, the thermodynamic profiling suggested the methionine residues, mainly M³²⁸, M³⁶⁵, and M⁴⁵⁰, seem to be responsible for the recruitment of PA. The other residues I³⁵⁷, Y³⁶¹, L⁴⁶⁵, F³⁰⁸, Q³¹⁶, and K³²¹ are also identified, exclusively interacting with PA. This study offers novel structural and mechanistic insights that are critical for FXR targeted drug discovery for PA designing.

konnect2prot: a web application to explore the protein properties in a functional protein-protein interaction network

MOTIVATION

The regulation of proteins governs the biological processes and functions and, therefore, the organisms' phenotype. So there is an unmet need for a systematic tool for identifying the proteins that play a crucial role in information processing in a protein-protein interaction (PPI) network. However, the current protein databases and web servers still lag behind to provide an end-to-end pipeline that can leverage the topological understanding of a context-specific PPI network to identify the influential spreaders. Addressing this, we developed a web application, 'konnect2prot' (k2p), which can generate context-specific directional PPI network from the input proteins and detect their biological and topological importance in the network.

RESULTS

We pooled together a large amount of ontological knowledge, parsed it down into a functional network, and gained insight into the molecular underpinnings of the disease development by creating a one-stop junction for PPI data. k2p contains both local and global information about a protein, such as protein class, disease mutations, ligands and PDB structure, enriched processes and pathways, multi-disease interactome and hubs and bottlenecks in the directional network. It also identifies spreaders in the network and maps them to disease hallmarks to determine whether they can affect the disease state or not.

AVAILABILITY AND IMPLEMENTATION

konnect2prot is freely accessible using the link https://konnect2prot.thsti.in. The code repository is https://github.com/samrat-lab/k2p_bioinfo-2022.

Characterizing (un)binding mechanism of USP7 inhibitors to unravel the cause of enhanced binding potencies at allosteric checkpoint

The ability to predict the intricate mechanistic behavior of ligands and associated structural determinants during protein-ligand (un)binding is of great practical importance in drug discovery. Ubiquitin specific protease-7 (USP7) is a newly emerging attractive cancer therapeutic target with bound allosteric inhibitors. However, none of the inhibitors have reached clinical trials, allowing opportunities to examine every aspect of allosteric modulation. The crystallographic insights reveal that these inhibitors have common properties such as chemical scaffolds, binding site and interaction fingerprinting. However, they still possess a broader range of binding potencies, ranging from 22 nM to 1,300 nM. Hence, it becomes more critical to decipher the structural determinants guiding the enhanced binding potency of the inhibitors. In this regard, we elucidated the atomic-level insights from both interacting partners, that is, protein-ligand perspective, and established the structure-activity link between USP7 inhibitors by using classical and advanced molecular dynamics simulations combined with linear interaction energy and molecular mechanics-Poisson Boltzmann surface area. We revealed the inhibitor potency differences by examining the contributions of chemical moieties and USP7 residues, the involvement of water-mediated interactions, and the thermodynamic landscape alterations. Additionally, the dissociation profiles aided in the establishment of a correlation between experimental potencies and structural determinants. Our study demonstrates the critical role of blocking loop 1 in allosteric inhibition and enhanced binding affinity. Comprehensively, our findings provide a constructive expansion of experimental outcomes and show the basis for varying binding potency using in-silico approaches. We expect this atomistic approach to be useful for effective drug design.

A Single Nucleotide Mixture Enhances the Antitumor Activity of Molecular-Targeted Drugs Against Hepatocellular Carcinoma

New strategies for molecular-targeted drug therapy for advanced hepatocellular carcinoma (HCC) ignore the contribution of the nutritional status of patients and nutritional support to improve physical status and immunity. We aimed to elucidate the role of a single nucleotide mixture (SNM) in the anti-tumor therapy of HCC, and to explore the importance of a SNM as adjuvant therapy for HCC. Compared with a lipid emulsion (commonly used nutritional supplement for HCC patients), the SNM could not induce metabolic abnormalities in HCC cells (Warburg effect), and did not affect expression of metabolic abnormality-related factors in HCC cells. The SNM could also attenuate the lymphocyte injury induced by antitumor drugs in vitro and in vivo, and promote the recruitment and survival of lymphocytes in HCC tissues. Using HCC models in SCID (server combined immune-deficiency) mice or BalB/c mice, the SNM had anti-tumor activity, and could significantly upregulate the antitumor activity of molecular-targeted drugs (tyrosine-kinase inhibitors [TKI] and immune-checkpoint inhibitors [ICI]) against HCC. We employed research models in vivo and in vitro to reveal the anti-tumor activity of the SNM on HCC. Our findings expand understanding of the SNM and contribute to HCC (especially nutritional support) therapy.