Covalent complexes of proteasome model with peptide aldehyde inhibitors MG132 and MG101: docking and molecular dynamics study

Synthesis, Characterization, and Biological Evaluation of Meldrum's Acid Derivatives: Dual Activity and Molecular Docking Study

In the presented study, eight novel Meldrum's acid derivatives containing various vanillic groups were synthesized. Vanillidene Meldrum's acid compounds were tested against different cancer cell lines and microbes. Out of nine, three showed very good biological activity against E. coli, and HeLa and A549 cell lines. It is shown that the O-alkyl substituted derivatives possessed better antimicrobial and anticancer activities in comparison with the O-acyl ones. The decyl substituted molecule (3i) has the highest activity against E. coli (MIC = 12.4 μM) and cancer cell lines (HeLa, A549, and LS174 = 15.7, 21.8, and 30.5 μM, respectively). The selectivity index of 3i is 4.8 (HeLa). The molecular docking study indicates that compound 3i showed good binding affinity to DNA, E. coli Gyrase B, and topoisomerase II beta. The covalent docking showed that 3i was a Michael acceptor for the nucleophiles Lys and Ser. The best Eb was noted for the topoisomerase II beta-LYS482-3i cluster.

Advances in covalent drug discovery

Covalent drugs have been used to treat diseases for more than a century, but tools that facilitate the rational design of covalent drugs have emerged more recently. The purposeful addition of reactive functional groups to existing ligands can enable potent and selective inhibition of target proteins, as demonstrated by the covalent epidermal growth factor receptor (EGFR) and Bruton's tyrosine kinase (BTK) inhibitors used to treat various cancers. Moreover, the identification of covalent ligands through 'electrophile-first' approaches has also led to the discovery of covalent drugs, such as covalent inhibitors for KRAS(G12C) and SARS-CoV-2 main protease. In particular, the discovery of KRAS(G12C) inhibitors validates the use of covalent screening technologies, which have become more powerful and widespread over the past decade. Chemoproteomics platforms have emerged to complement covalent ligand screening and assist in ligand discovery, selectivity profiling and target identification. This Review showcases covalent drug discovery milestones with emphasis on the lessons learned from these programmes and how an evolving toolbox of covalent drug discovery techniques facilitates success in this field.

Phage-Related Ribosomal Protease (Prp) of Staphylococcus aureus: In Vitro Michaelis-Menten Kinetics, Screening for Inhibitors, and Crystal Structure of a Covalent Inhibition Product Complex

Phage-related ribosomal proteases (Prps) are essential for the assembly and maturation of the ribosome in Firmicutes, including the human pathogens Staphylococcus aureus, Streptococcus pneumoniae, and Clostridium difficile. These bacterial proteases cleave off an N-terminal extension of a precursor of ribosomal protein L27, a processing step that is essential for the formation of functional ribosomes. This essential role of Prp in these pathogens has identified this protease as a potential antibiotic target. In this work, we determine the X-ray crystal structure of a covalent inhibition complex at 2.35 Å resolution, giving the first complete picture of the active site of a functional Prp. We also characterize the kinetic activity and screen for potential inhibitors of Prp. This work gives the most complete characterization of the structure and specificity of this novel class of proteases to date.

Potential drug discovery for COVID-19 treatment targeting Cathepsin L using a deep learning-based strategy

Cathepsin L (CTSL), a cysteine protease that can cleave and activate the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein, could be a promising therapeutic target for coronavirus disease 2019 (COVID-19). However, there is still no clinically available CTSL inhibitor that can be used. Here, we applied Chemprop, a newly trained directed-message passing deep neural network approach, to identify small molecules and FDA-approved drugs that can block CTSL activity to expand the discovery of CTSL inhibitors for drug development and repurposing for COVID-19. We found 5 molecules (Mg-132, Z-FA-FMK, leupeptin hemisulfate, Mg-101 and calpeptin) that were able to significantly inhibit the activity of CTSL in the nanomolar range and inhibit the infection of both pseudotype and live SARS-CoV-2. Notably, we discovered that daptomycin, an FDA-approved antibiotic, has a prominent CTSL inhibitory effect and can inhibit SARS-CoV-2 pseudovirus infection. Further, molecular docking calculation showed stable and robust binding of these compounds with CTSL. In conclusion, this study suggested for the first time that Chemprop is ideally suited to predict additional inhibitors of enzymes and revealed the noteworthy strategy for screening novel molecules and drugs for the treatment of COVID-19 and other diseases with unmet needs.

Repurposing FDA-approved Drugs Targeting SARS-CoV2 3CLpro: a study by applying Virtual Screening, Molecular Dynamics, MM-PBSA Calculations and Covalent Docking

Background:

Since the end of 2019, the etiologic agent SAR-CoV-2 responsible for one of the most significant epidemics in history has caused severe global economic, social, and health damages. The drug repurposing approach and application of Structure-based Drug Discovery (SBDD) using in silico techniques are increasingly frequent, leading to the identification of several molecules that may represent promising potential.

Method:

In this context, here we use in silico methods of virtual screening (VS), pharmacophore modeling (PM), and fragment-based drug design (FBDD), in addition to molecular dynamics (MD), molecular mechanics/Poisson-Boltzmann surface area (MM -PBSA) calculations, and covalent docking (CD) for the identification of potential treatments against SARS-CoV-2. We initially validated the docking protocol followed by VS in 1,613 FDA-approved drugs obtained from the ZINC database. Thus, we identified 15 top hits, of which three of them were selected for further simulations. In parallel, for the compounds with a fit score value ≤ of 30, we performed the FBDD protocol, where we designed 12 compounds

Result:

By applying a PM protocol in the ZINC database, we identified three promising drug candidates. Then, the 9 top hits were evaluated in simulations of MD, MM-PBSA, and CD. Subsequently, MD showed that all identified hits showed stability at the active site without significant changes in the protein's structural integrity, as evidenced by the RMSD, RMSF, Rg, SASA graphics. They also showed interactions with the catalytic dyad (His41 and Cys145) and other essential residues for activity (Glu166 and Gln189) and high affinity for MM-PBSA, with possible covalent inhibition mechanism.

Conclution:

Finally, our protocol helped identify potential compounds wherein ZINC896717 (Zafirlukast), ZINC1546066 (Erlotinib), and ZINC1554274 (Rilpivirine) were more promising and could be explored in vitro, in vivo, and clinical trials to prove their potential as antiviral agents.

Exploring EZH2-Proteasome Dual-Targeting Drug Discovery through a Computational Strategy to Fight Multiple Myeloma

Multiple myeloma is an incurable plasma cell neoplastic disease representing about 10-15% of all haematological malignancies diagnosed in developed countries. Proteasome is a key player in multiple myeloma and proteasome inhibitors are the current first-line of treatment. However, these are associated with limited clinical efficacy due to acquired resistance. One of the solutions to overcome this problem is a polypharmacology approach, namely combination therapy and multitargeting drugs. Several polypharmacology avenues are currently being explored. The simultaneous inhibition of EZH2 and Proteasome 20S remains to be investigated, despite the encouraging evidence of therapeutic synergy between the two. Therefore, we sought to bridge this gap by proposing a holistic in silico strategy to find new dual-target inhibitors. First, we assessed the characteristics of both pockets and compared the chemical space of EZH2 and Proteasome 20S inhibitors, to establish the feasibility of dual targeting. This was followed by molecular docking calculations performed on EZH2 and Proteasome 20S inhibitors from ChEMBL 25, from which we derived a predictive model to propose new EZH2 inhibitors among Proteasome 20S compounds, and vice versa, which yielded two dual-inhibitor hits. Complementarily, we built a machine learning QSAR model for each target but realised their application to our data is very limited as each dataset occupies a different region of chemical space. We finally proceeded with molecular dynamics simulations of the two docking hits against the two targets. Overall, we concluded that one of the hit compounds is particularly promising as a dual-inhibitor candidate exhibiting extensive hydrogen bonding with both targets. Furthermore, this work serves as a framework for how to rationally approach a dual-targeting drug discovery project, from the selection of the targets to the prediction of new hit compounds.

Virtual screening using covalent docking to find activators for G245S mutant p53

TP53 is the most mutated gene in all cancers. The mutant protein also accumulates in cells. The high frequency of p53 mutations makes the protein a promising target for anti-cancer therapy. Only a few molecules have been found, using in vitro screening, to reactivate the mutant protein. APR-246 is currently the most successful mutant p53 activator, which reactivates the transcriptional activity of p53 by covalently binding to C124 of the protein. We have recently created in silico models of G245S-mp53 in its apo and DNA-bound forms. In this paper we further report on our in silico screening for potential activators of G245S-mp53. We filtered the ZINC15 database (13 million compounds) to only include drug-like molecules with moderate to standard reactivity. Our filtered database of 130,000 compounds was screened using the DOCKTITE protocol in the Molecular Operating Environment software. We performed covalent docking at C124 of G245S-mp53 to identify potential activators of the mutant protein. The docked compounds were ranked using a consensus scoring approach. We also used ADMET Predictor™ to predict pharmacokinetics and the possible toxicities of the compounds. Our screening procedure has identified compounds, mostly thiosemicarbazones and halo-carbonyls, with the best potential as G245S-mp53 activators, which are described in this work. Based on its binding scores and ADMET risk score, compound 2 is likely to have the best potential as a G245S-mp53 activator compared to the other top hits.

Docking of Covalent Ligands: Challenges and Approaches

Covalent ligands have recently regained considerable attention in drug discovery. The rational design of such ligands, however, is still faced with particular challenges, mostly related to the fact that covalent bond formation is a quantum mechanical phenomenon which cannot adequately be handled by the force fields or empirical approaches typically used for noncovalent protein-ligand interactions. Although the necessity for quantum chemical approaches is clear, they cannot yet routinely be applied on large data sets of ligands or for a broader exploration of binding modes in docking calculations. On the other hand, technical solutions for performing docking calculations with covalent ligands are available, but their scope is normally quite limited. Scoring functions typically neglect the contribution from covalent bond formation completely. In this situation, the question arises how to approach covalent ligands and which methods to choose for their docking and design.

Proteasome properties of hemocytes differ between the whiteleg shrimp Penaeus vannamei and the brown shrimp Crangon crangon (Crustacea, Decapoda)

Crustaceans are intensively farmed in aquaculture facilities where they are vulnerable to parasites, bacteria, or viruses, often severely compromising the rearing success. The ubiquitin-proteasome system (UPS) is crucial for the maintenance of cellular integrity. Analogous to higher vertebrates, the UPS of crustaceans may also play an important role in stress resistance and pathogen defense. We studied the general properties of the proteasome system in the hemocytes of the whiteleg shrimp, Penaeus vannamei, and the European brown shrimp Crangon crangon. The 20S proteasome was the predominant proteasome population in the hemocytes of both species. The specific activities of the trypsin-like (Try-like), chymotrypsin-like (Chy-like), and caspase-like (Cas-like) enzymes of the shrimp proteasome differed between species. P. vannamei exhibited a higher ratio of Try-like to Chy-like activities and Cas-like to Chy-like activities than C. crangon. Notably, the Chy-like activity of P. vannamei showed substrate or product inhibition at concentrations of more than 25 mmol L-1. The K M values ranged from 0.072 mmol L-1 for the Try-like activity of P. vannamei to 0.309 mmol L-1 for the Cas-like activity of C. crangon. Inhibition of the proteasome of P. vannamei by proteasome inhibitors was stronger than in C. crangon. The pH profiles were similar in both species. The Try-like, Chy-like, and Cas-like sites showed the highest activities between pH 7.5 and 8.5. The proteasomes of both species were sensitive against repeated freezing and thawing losing ~80-90% of activity. This study forms the basis for future investigations on the shrimp response against infectious diseases, and the role of the UPS therein.

Design, synthesis and biological evaluation of novel non-peptide boronic acid derivatives as proteasome inhibitors

A novel series of non-peptide proteasome inhibitors bearing the 1, 4-naphthoquinone scaffold and boronic acid warhead was developed. In the biological evaluation on the chymotrypsin-like activity of human 20S proteasome, five compounds showed IC₅₀ values in the nanomolar range. Docking experiments into the yeast 20S proteasome rationalized their biological activities and allowed further optimization of this interesting class of inhibitors. Within the cellular proliferation inhibition assay and western blot analysis, compound 3e demonstrated excellent anti-proliferative activity against solid tumor cells and clear accumulation of ubiquitinated cellular proteins. Furthermore, in the microsomal stability assay compound 3e demonstrated much improved metabolic stability compared to bortezomib, emerging as a promising lead compound for further design of non-peptide proteasome inhibitors.

Computational Approaches for the Discovery of Human Proteasome Inhibitors: An Overview

Proteasome emerged as an important target in recent pharmacological research due to its pivotal role in degrading proteins in the cytoplasm and nucleus of eukaryotic cells, regulating a wide variety of cellular pathways, including cell growth and proliferation, apoptosis, DNA repair, transcription, immune response, and signaling processes. The last two decades witnessed intensive efforts to discover 20S proteasome inhibitors with significant chemical diversity and efficacy. To date, the US FDA approved to market three proteasome inhibitors: bortezomib, carfilzomib, and ixazomib. However new, safer and more efficient drugs are still required. Computer-aided drug discovery has long being used in drug discovery campaigns targeting the human proteasome. The aim of this review is to illustrate selected in silico methods like homology modeling, molecular docking, pharmacophore modeling, virtual screening, and combined methods that have been used in proteasome inhibitors discovery. Applications of these methods to proteasome inhibitors discovery will also be presented and discussed to raise improvements in this particular field.

Covalent Inhibition of HIV-1 Integrase by N-Succinimidyl Peptides

We present a new approach for the covalent inhibition of HIV-1 integrase (IN) by an LEDGF/p75-derived peptide modified with an N-terminal succinimide group. The covalent inhibition is mediated by direct binding of the succinimide to the amine group of a lysine residue in IN. The peptide serves as a specific recognition sequence for the target protein, while the succinimide serves as the binding moiety. The combination of a readily synthesizable peptide precursor with easy and efficient binding to the target protein makes this approach a promising new strategy for designing lead compounds.

Crystal structure of N-{N-[N-acetyl-(S)-leuc-yl]-(S)-leuc-yl}norleucinal (ALLN), an inhibitor of proteasome

The crystal structure of ALLN, the tripeptidic inhibitor of proteasomes, is solved from synchrotron diffraction data. An infinite β-sheet extended through the crystal is formed by symmetry-related oligopeptide mol­ecules in extended conformation.

The title compound, C₂₀H₃₇N₃O₄, also known by the acronym ALLN, is a tripeptidic inhibitor of the proteolytic activity of the proteasomes, enzyme complexes implicated in several neurodegenerative diseases and other disorders, including cancer. The crystal structure of ALLN, solved from synchrotron radiation diffraction data, revealed the mol­ecules in extended conformation of the backbone and engaging all peptide N and O atoms in inter­molecular hydrogen bonds forming an infinite anti­parallel β-sheet.

Theory and applications of covalent docking in drug discovery: merits and pitfalls

The present art of drug discovery and design of new drugs is based on suicidal irreversible inhibitors. Covalent inhibition is the strategy that is used to achieve irreversible inhibition. Irreversible inhibitors interact with their targets in a time-dependent fashion, and the reaction proceeds to completion rather than to equilibrium. Covalent inhibitors possessed some significant advantages over non-covalent inhibitors such as covalent warheads can target rare, non-conserved residue of a particular target protein and thus led to development of highly selective inhibitors, covalent inhibitors can be effective in targeting proteins with shallow binding cleavage which will led to development of novel inhibitors with increased potency than non-covalent inhibitors. Several computational approaches have been developed to simulate covalent interactions; however, this is still a challenging area to explore. Covalent molecular docking has been recently implemented in the computer-aided drug design workflows to describe covalent interactions between inhibitors and biological targets. In this review we highlight: (i) covalent interactions in biomolecular systems; (ii) the mathematical framework of covalent molecular docking; (iii) implementation of covalent docking protocol in drug design workflows; (iv) applications covalent docking: case studies and (v) shortcomings and future perspectives of covalent docking. To the best of our knowledge; this review is the first account that highlights different aspects of covalent docking with its merits and pitfalls. We believe that the method and applications highlighted in this study will help future efforts towards the design of irreversible inhibitors.

Discovery of novel covalent proteasome inhibitors through a combination of pharmacophore screening, covalent docking, and molecular dynamics simulations

The ubiquitin-proteasome pathway plays a pivotal role in the regulation of cellular protein processing and degradation. Proteasome inhibitors (PIs) have enormous potential to treat multiple myeloma, solid tumors, parasites, inflammation, and immune diseases, which is spurring the development of new types of PIs with enhanced efficacy, fewer side effects, and reduced drug resistance. Nevertheless, virtual screening for covalent PIs has rarely been reported because calculating the covalent binding energy is a challenging task. The aim of this study was to discover new covalent inhibitors of the 20S proteasome. The structures of PIs were manually divided into two parts: a noncovalent binding part resulting from virtual screening, and an epoxyketone group that was pre-selected as a covalent binding part. The SPECS database was screened by noncovalent docking and a pharmacophore model built with the 20S proteasome. After validating the covalent conjugation, 88 hits with epoxyketone were covalently docked into the 20S proteasome to analyze the intermolecular interactions. Four compounds were selected after multiple filtration and validations. Molecular dynamics simulations were performed to check the stability of the noncovalent and covalent docked ligand-enzyme complexes and investigate the interaction patterns of the screened inhibitors. Finally, two compounds with novel aromatic backbones, reasonable interactions, and stable covalent binding modes were retained. These compounds can serve as potential hits for further biological evaluation.

New C(4)- and C(1)-derivatives of furo[3,4-c]pyridine-3-ones and related compounds: evidence for site-specific inhibition of the constitutive proteasome and its immunoisoform

A set of 18 new C(4) and C(1) derivatives of nor-cerpegin (1,1-dimethyl furo[3,4-c]pyridine-3-one), 6 model compounds (γ- and δ-lactones) and 20 furo- or thieno[2,3-d]-pyrimidine-4-one related compounds were designed and synthesized. Each compound was assayed for inhibition of CT-L, T-L and PA proteolytic activities of 20S constitutive proteasome (c20S). Most performant compounds were also assayed on 20S immunoproteasome (i20S). Compound 10 with a benzylamino group at C(4) and dimethylated at C(1) of the furopyridine ring was the most efficient PA site-specific inhibitor of the c20S (IC50(cPA) of 600nM) without noticeable inhibition of the i20S PA site (iPA). In silico docking assays for 10 at the iPA catalytic site revealed the absence of poses normally observed for this compound and related ones at the constitutive PA site (cPA). The thieno[2,3-d]pyrimidine-4-one 40 was T-L site-specific with a mild inhibition of both c20S and i20S in vitro (IC50(cT-L) of 9.9μM and IC50(iT-L) of 6.7μM). In silico docking assays of 40 at T-L sites of c20S and i20S revealed almost identical first rank poses in the two types of sites with no possibility left for nucleophilic attack by Thr1 as observed for the fused furopyridine-3-one 10.

ROS inhibitor N-acetyl-L-cysteine antagonizes the activity of proteasome inhibitors

NAC (N-acetyl-L-cysteine) is commonly used to identify and test ROS (reactive oxygen species) inducers, and to inhibit ROS. In the present study, we identified inhibition of proteasome inhibitors as a novel activity of NAC. Both NAC and catalase, another known scavenger of ROS, similarly inhibited ROS levels and apoptosis associated with H₂O₂. However, only NAC, and not catalase or another ROS scavenger Trolox, was able to prevent effects linked to proteasome inhibition, such as protein stabilization, apoptosis and accumulation of ubiquitin conjugates. These observations suggest that NAC has a dual activity as an inhibitor of ROS and proteasome inhibitors. Recently, NAC was used as a ROS inhibitor to functionally characterize a novel anticancer compound, piperlongumine, leading to its description as a ROS inducer. In contrast, our own experiments showed that this compound depicts features of proteasome inhibitors including suppression of FOXM1 (Forkhead box protein M1), stabilization of cellular proteins, induction of ROS-independent apoptosis and enhanced accumulation of ubiquitin conjugates. In addition, NAC, but not catalase or Trolox, interfered with the activity of piperlongumine, further supporting that piperlongumine is a proteasome inhibitor. Most importantly, we showed that NAC, but not other ROS scavengers, directly binds to proteasome inhibitors. To our knowledge, NAC is the first known compound that directly interacts with and antagonizes the activity of proteasome inhibitors. Taken together, the findings of the present study suggest that, as a result of the dual nature of NAC, data interpretation might not be straightforward when NAC is utilized as an antioxidant to demonstrate ROS involvement in drug-induced apoptosis.

Microarray analysis of isolated human islet transcriptome in type 2 diabetes and the role of the ubiquitin-proteasome system in pancreatic beta cell dysfunction

To shed light on islet cell molecular phenotype in human type 2 diabetes (T2D), we studied the transcriptome of non-diabetic (ND) and T2D islets to then focus on the ubiquitin-proteasome system (UPS), the major protein degradation pathway. We assessed gene expression, amount of ubiquitinated proteins, proteasome activity, and the effects of proteasome inhibition and prolonged exposure to palmitate. Microarray analysis identified more than one thousand genes differently expressed in T2D islets, involved in many structures and functions, with consistent alterations of the UPS. Quantitative RT-PCR demonstrated downregulation of selected UPS genes in T2D islets and beta cell fractions, with greater ubiquitin accumulation and reduced proteasome activity. Chemically induced reduction of proteasome activity was associated with lower glucose-stimulated insulin secretion, which was partly reproduced by palmitate exposure. These results show the presence of many changes in islet transcriptome in T2D islets and underline the importance of the association between UPS alterations and beta cell dysfunction in human T2D.

C1 and N5 derivatives of cerpegin: synthesis of a new series based on structure-activity relationships to optimize their inhibitory effect on 20S proteasome

Thirty-two new derivatives of cerpegin (1,1,5-trimethylfuro[3,4-c]pyridine-3,4-dione) were designed and synthesized in high yield by a new method, combining several C(1) and N(5) substituents. All compounds were tested for their inhibitory effect on the CT-L, T-L and PA proteolytic activities of a purified mammalian 20S proteasome. Only one molecule inhibited both CT-L and PA activities. Sixteen molecules specifically inhibited PA at the micromolar range, out of which fourteen had IC50 values around 5 μM and two had IC50 values closer to 2 μM. Except in one case, neither calpain I nor cathepsin B was inhibited. In silico docking suggests a unique mode of binding of the most efficient compounds to the β1 catalytic site (PA activity) in relation to the chemical nature of C(1) substituents.

Docking-based virtual screening of covalently binding ligands: an orthogonal lead discovery approach

In pharmaceutical industry, lead discovery strategies and screening collections have been predominantly tailored to discover compounds that modulate target proteins through noncovalent interactions. Conversely, covalent linkage formation is an important mechanism for a quantity of successful drugs in the market, which are discovered in most cases by hindsight instead of systematical design. In this article, the implementation of a docking-based virtual screening workflow for the retrieval of covalent binders is presented considering human cathepsin K as a test case. By use of the docking conditions that led to the best enrichment of known actives, 44 candidate compounds with unknown activity on cathepsin K were finally selected for experimental evaluation. The most potent inhibitor, 4-(N-phenylanilino)-6-pyrrolidin-1-yl-1,3,5-triazine-2-carbonitrile (CP243522), showed a K(i) of 21 nM and was confirmed to have a covalent reversible mechanism of inhibition. The presented approach will have great potential in cases where covalent inhibition is the desired drug discovery strategy.