Design of inhibitors of SARS-CoV-2 papain-like protease deriving from GRL0617: Structure-activity relationships.

The unique and complex structure of papain-like protease (PLpro) of the SARS-CoV-2 virus represents a difficult challenge for antiviral development, yet it offers a compelling validated target for effective therapy of COVID-19. The surge in scientific interest in inhibiting this cysteine protease emerged after its demonstrated connection to the cytokine storm in patients with COVID-19 disease. Furthermore, the development of new inhibitors against PLpro may also be beneficial for the treatment of respiratory infections caused by emerging coronavirus variants of concern. This review article provides a comprehensive overview of PLpro inhibitors, focusing on the structural framework of the known inhibitor GRL0617 and its analogs. We categorize PLpro inhibitors on the basis of their structures and binding site: Glu167 containing site, BL2 groove, Val70Ub site, and Cys111 containing catalytic site. We summarize and evaluate the majority of GRL0617-like inhibitors synthesized so far, highlighting their published biochemical parameters, which reflect their efficacy. Published research has shown that strategic modifications to GRL0617, such as decorating the naphthalene ring, extending the aromatic amino group or the orthomethyl group, can substantially decrease the IC50 from micromolar up to nanomolar concentration range. Some advantageous modifications significantly enhance inhibitory activity, paving the way for the development of new potent compounds. Our review places special emphasis on structures that involve direct modifications to the GRL0617 scaffold, including piperidine carboxamides and modified benzylmethylnaphthylethanamines (Jun9 scaffold). All these compounds are believed to inhibit the proteolytic, deubiquitination, and deISGylation activity of PLpro, biochemical processes linked to the severe progression of COVID-19. Finally, we summarize the development efforts for SARS-CoV-2 PLpro inhibitors, in detailed structure-activity relationships diagrams. This aims to inform and inspire future research in the search for potent antiviral agents against PLpro of current and emerging coronavirus threats.

Development of Novel Peptidyl Nitriles Targeting Rhodesain and Falcipain-2 for the Treatment of Sleeping Sickness and Malaria.

In recent decades, neglected tropical diseases and poverty-related diseases have become a serious health problem worldwide. Among these pathologies, human African trypanosomiasis, and malaria present therapeutic problems due to the onset of resistance, toxicity problems and the limited spectrum of action. In this drug discovery process, rhodesain and falcipain-2, of Trypanosoma brucei rhodesiense and Plasmodium falciparum, are currently considered the most promising targets for the development of novel antitrypanosomal and antiplasmodial agents, respectively. Therefore, in our study we identified a novel lead-like compound, i.e., inhibitor 2b, which we proved to be active against both targets, with a Ki = 5.06 µM towards rhodesain and an IC50 = 40.43 µM against falcipain-2.

Substitution-Induced Mechanistic Switching in SNAr-Warheads for Cysteine Proteases.

The aim of this study was to investigate the transition from non-covalent reversible over covalent reversible to covalent irreversible inhibition of cysteine proteases by making delicate structural changes to the warhead scaffold. To this end, dipeptidic rhodesain inhibitors with different N-terminal electrophilic arenes as warheads relying on the SNAr mechanism were synthesized and investigated. Strong structure-activity relationships of the inhibition potency, the degree of covalency, and the reversibility of binding on the arene substitution pattern were found. The studies were complemented and substantiated by molecular docking and quantum-mechanical calculations of model systems. Furthermore, the improvement in the membrane permeability of peptide esters in comparison to their corresponding carboxylic acids was exemplified.

Structure-Aided Computational Design of Triazole-Based Targeted Covalent Inhibitors of Cruzipain.

Cruzipain (CZP), the major cysteine protease present in T. cruzi, the ethiological agent of Chagas disease, has attracted particular attention as a therapeutic target for the development of targeted covalent inhibitors (TCI). The vast chemical space associated with the enormous molecular diversity feasible to explore by means of modern synthetic approaches allows the design of CZP inhibitors capable of exhibiting not only an efficient enzyme inhibition but also an adequate translation to anti-T. cruzi activity. In this work, a computer-aided design strategy was developed to combinatorially construct and screen large libraries of 1,4-disubstituted 1,2,3-triazole analogues, further identifying a selected set of candidates for advancement towards synthetic and biological activity evaluation stages. In this way, a virtual molecular library comprising more than 75 thousand diverse and synthetically feasible analogues was studied by means of molecular docking and molecular dynamic simulations in the search of potential TCI of CZP, guiding the synthetic efforts towards a subset of 48 candidates. These were synthesized by applying a Cu(I)-catalyzed azide-alkyne cycloaddition (CuAAC) centered synthetic scheme, resulting in moderate to good yields and leading to the identification of 12 hits selectively inhibiting CZP activity with IC50 in the low micromolar range. Furthermore, four triazole derivatives showed good anti-T. cruzi inhibition when studied at 50 µM; and Ald-6 excelled for its high antitrypanocidal activity and low cytotoxicity, exhibiting complete in vitro biological activity translation from CZP to T. cruzi. Overall, not only Ald-6 merits further advancement to preclinical in vivo studies, but these findings also shed light on a valuable chemical space where molecular diversity might be explored in the search for efficient triazole-based antichagasic agents.

Structure-Based Optimization of Protease-Inhibitor Interactions to Enhance Specificity of Human Stefin-A against Falcipain-2 from the Plasmodium falciparum 3D7 Strain.

The emergence of resistance in Plasmodium falciparum to frontline artemisinin-based combination therapies has raised global concerns and emphasized the identification of new drug targets for malaria. Cysteine protease falcipain-2 (FP2), involved in host hemoglobin degradation and instrumental in parasite survival, has long been proposed as a promising malarial drug target. However, designing active-site-targeted small-molecule inhibitors of FP2 becomes challenging due to their off-target specificity toward highly homologous human cysteine cathepsins. The use of proteinaceous inhibitors, which have nonconserved exosite interactions and larger interface area, can effectively circumvent this problem. In this study, we report for the first time that human stefin-A (STFA) efficiently inhibits FP2 with Ki values in the nanomolar range. The FP2-STFA complex crystal structure, determined in this study, and sequence analyses identify a unique nonconserved exosite interaction, compared to human cathepsins. Designing a mutation Lys68 > Arg in STFA amplifies its selectivity garnering a 3.3-fold lower Ki value against FP2, and the crystal structure of the FP2-STFAK68R complex shows stronger electrostatic interaction between side-chains of Arg68 (STFAK68R) and Asp109 (FP2). Comparative structural analyses and molecular dynamics (MD) simulation studies of the complexes further confirm higher buried surface areas, better interaction energies for FP2-STFAK68R, and consistency of the newly developed electrostatic interaction (STFA-R68-FP2-D109) in the MD trajectory. The STFA-K68R mutant also shows higher Ki values against human cathepsin-L and stefin, a step toward eliminating off-target specificity. Hence, this work underlines the design of host-based proteinaceous inhibitors against FP2, with further optimization to render them more potent and selective.

Experimental and Computational Study of Aryl-thiosemicarbazones Inhibiting Cruzain Reveals Reversible Inhibition and a Stepwise Mechanism.

Trypanosoma cruzi is a parasite that infects about 6-7 million people worldwide, mostly in Latin America, causing Chagas disease. Cruzain, the main cysteine protease of T. cruzi, is a validated target for developing drug candidates for Chagas disease. Thiosemicarbazones are one of the most relevant warheads used in covalent inhibitors targeting cruzain. Despite its relevance, the mechanism of inhibition of cruzain by thiosemicarbazones is unknown. Here, we combined experiments and simulations to unveil the covalent inhibition mechanism of cruzain by a thiosemicarbazone-based inhibitor (compound 1). Additionally, we studied a semicarbazone (compound 2), which is structurally similar to compound 1 but does not inhibit cruzain. Assays confirmed the reversibility of inhibition by compound 1 and suggested a two-step mechanism of inhibition. The Ki was estimated to be 36.3 µM and Ki* to be 11.5 µM, suggesting the pre-covalent complex to be relevant for inhibition. Molecular dynamics simulations of compounds 1 and 2 with cruzain were used to propose putative binding modes for the ligands. One-dimensional (1D) quantum mechanics/molecular mechanics (QM/MM) potential of mean force (PMF) and gas-phase energies showed that the attack of Cys25-S- on the C═S or C═O bond yields a more stable intermediate than the attack on the C═N bond of the thiosemicarbazone/semicarbazone. Two-dimensional (2D) QM/MM PMF revealed a putative reaction mechanism for compound 1, involving the proton transfer to the ligand, followed by the Cys25-S- attack at C═S. The ΔG and energy barrier were estimated to be -1.4 and 11.7 kcal/mol, respectively. Overall, our results shed light on the inhibition mechanism of cruzain by thiosemicarbazones.

Molecular Features of CA-074 pH-Dependent Inhibition of Cathepsin B.

CA-074 is a selective inhibitor of cathepsin B, a lysosomal cysteine protease. CA-074 has been utilized in numerous studies to demonstrate the role of this protease in cellular and physiological functions. Cathepsin B in numerous human disease mechanisms involves its translocation from acidic lysosomes of pH 4.6 to neutral pH 7.2 of cellular locations, including the cytosol and extracellular environment. To gain in-depth knowledge of CA-074 inhibition under these different pH conditions, this study evaluated the molecular features, potency, and selectivity of CA-074 for cathepsin B inhibition under acidic and neutral pH conditions. This study demonstrated that CA-074 is most effective at inhibiting cathepsin B at an acidic pH of 4.6 with nM potency, which was more than 100-fold more potent than its inhibition at a neutral pH of 7.2. The pH-dependent inhibition of CA-074 was abolished by methylation of its C-terminal proline, indicating the requirement for the free C-terminal carboxyl group for pH-dependent inhibition. Under these acidic and neutral pH conditions, CA-074 maintained its specificity for cathepsin B over other cysteine cathepsins, displayed irreversible inhibition, and inhibited diverse cleavages of peptide substrates of cathepsin B assessed by profiling mass spectrometry. Molecular docking suggested that pH-dependent ionic interactions of the C-terminal carboxylate of CA-074 occur with His110 and His111 residues in the S2' subsite of the enzyme at pH 4.6, but these interactions differ at pH 7.2. While high levels of CA-074 or CA-074Me (converted by cellular esterases to CA-074) are used in biological studies to inhibit cathepsin B at both acidic and neutral pH locations, it is possible that adjusted levels of CA-074 or CA-074Me may be explored to differentially affect cathepsin B activity at these different pH values. Overall, the results of this study demonstrate the molecular, kinetic, and protease specificity features of CA-074 pH-dependent inhibition of cathepsin B.

Inhibitors of Coronavirus 3CL Proteases Protect Cells from Protease-Mediated Cytotoxicity.

We describe a mammalian cell-based assay to identify coronavirus 3CL protease (3CLpro) inhibitors. This assay is based on rescuing protease-mediated cytotoxicity and does not require live virus. By enabling the facile testing of compounds across a range of 15 distantly related coronavirus 3CLpro enzymes, we identified compounds with broad 3CLpro-inhibitory activity. We also adapted the assay for use in compound screening and in doing so uncovered additional severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) 3CLpro inhibitors. We observed strong concordance between data emerging from this assay and those obtained from live-virus testing. The reported approach democratizes the testing of 3CLpro inhibitors by developing a simplified method for identifying coronavirus 3CLpro inhibitors that can be used by the majority of laboratories, rather than the few with extensive biosafety infrastructure. We identified two lead compounds, GC376 and compound 4, with broad activity against all 3CL proteases tested, including 3CLpro enzymes from understudied zoonotic coronaviruses. IMPORTANCE Multiple coronavirus pandemics have occurred over the last 2 decades. This has highlighted a need to be proactive in the development of therapeutics that can be readily deployed in the case of future coronavirus pandemics. We developed and validated a simplified cell-based assay for the identification of chemical inhibitors of 3CL proteases encoded by a wide range of coronaviruses. This assay is reporter free, does not require specialized biocontainment, and is optimized for performance in high-throughput screening. By testing reported 3CL protease inhibitors against a large collection of 3CL proteases with variable sequence similarity, we identified compounds with broad activity against 3CL proteases and uncovered structural insights into features that contribute to their broad activity. Furthermore, we demonstrated that this assay is suitable for identifying chemical inhibitors of proteases from families other than 3CL proteases.

Lead optimization of cathepsin K inhibitors for the treatment of Osteoarthritis.

Cathepsin K (Cat K) is a cysteine protease involved in bone remodeling. In addition to its role in bone biology, Cat K is upregulated in osteoclasts, chondrocytes and synoviocytes in osteoarthritic (OA) disease states making it a potential therapeutic target for disease-modifying OA. Starting from a prior preclinical compound, MK-1256, lead optimization efforts were carried out in the search for potent Cat K inhibitors with improved selectivity profiles with an emphasis on cathepsin F. Herein, we report the SAR studies which led to the discovery of the highly selective oxazole compound 23, which was subsequently shown to inhibit cathepsin K in vivo as measured by reduced levels of urinary C-telopeptide of collagen type I in dog.

Assessment of Reversibility for Covalent Cysteine Protease Inhibitors Using Quantum Mechanics/Molecular Mechanics Free Energy Surfaces.

We have used molecular dynamics (MD) simulations with hybrid quantum mechanics/molecular mechanics (QM/MM) potentials to investigate the reaction mechanism for covalent inhibition of cathepsin K and assess the reversibility of inhibition. The computed free energy profiles suggest that a nucleophilic attack by the catalytic cysteine on the inhibitor warhead and proton transfer from the catalytic histidine occur in a concerted manner. The results indicate that the reaction is more strongly exergonic for the alkyne-based inhibitors, which bind irreversibly to cathepsin K, than for the nitrile-based inhibitor odanacatib, which binds reversibly. Gas-phase energies were also calculated for the addition of methanethiol to structural prototypes for a number of warheads of interest in cysteine protease inhibitor design in order to assess electrophilicity. The approaches presented in this study are particularly applicable to assessment of novel warheads, and computed transition state geometries can be incorporated into molecular models for covalent docking.

Tetracycline Derivatives Inhibit Plasmodial Cysteine Protease Falcipain-2 through Binding to a Distal Allosteric Site.

Allosteric inhibitors regulate enzyme activity from remote and usually specific pockets. As they promise an avenue for less toxic and safer drugs, the identification and characterization of allosteric inhibitors has gained great academic and biomedical interest in recent years. Research on falcipain-2 (FP-2), the major papain-like cysteine hemoglobinase of Plasmodium falciparum, might benefit from this strategy to overcome the low selectivity against human cathepsins shown by active site-directed inhibitors. Encouraged by our previous finding that methacycline inhibits FP-2 noncompetitively, here we assessed other five tetracycline derivatives against this target and characterized their inhibition mechanism. As previously shown for methacycline, tetracycline derivatives inhibited FP-2 in a noncompetitive fashion, with Ki values ranging from 121 to 190 µM. A possible binding to the S' side of the FP-2 active site, similar to that described by X-ray crystallography (PDB: 6SSZ) for the noncompetitive inhibitor E-chalcone 48 (EC48), was experimentally discarded by kinetic analysis using a large peptidyl substrate spanning the whole active site. By combining lengthy molecular dynamics (MD) simulations that allowed methacycline to diffuse from solution to different FP-2 surface regions and free energy calculations, we predicted the most likely binding mode of the ligand. Of note, the proposed binding pose explains the low differences in Ki values observed for the tested tetracycline derivatives and the calculated binding free energies match the experimental values. Overall, this study has implications for the design of novel allosteric inhibitors against FP-2 and sets the basis for further optimization of the tetracycline scaffold to produce more potent and selective inhibitors.

Neoechinulin A as a Promising SARS-CoV-2 Mpro Inhibitor: In Vitro and In Silico Study Showing the Ability of Simulations in Discerning Active from Inactive Enzyme Inhibitors.

The COVID-19 pandemic and its continuing emerging variants emphasize the need to discover appropriate treatment, where vaccines alone have failed to show complete protection against the new variants of the virus. Therefore, treatment of the infected cases is critical. This paper discusses the bio-guided isolation of three indole diketopiperazine alkaloids, neoechinulin A (1), echinulin (2), and eurocristatine (3), from the Red Sea-derived Aspergillus fumigatus MR2012. Neoechinulin A (1) exhibited a potent inhibitory effect against SARS-CoV-2 Mpro with IC50 value of 0.47 µM, which is comparable to the reference standard GC376. Despite the structural similarity between the three compounds, only 1 showed a promising effect. The mechanism of inhibition is discussed in light of a series of extensive molecular docking, classical and steered molecular dynamics simulation experiments. This paper sheds light on indole diketopiperazine alkaloids as a potential structural motif against SARS-CoV-2 Mpro. Additionally, it highlights the potential of different molecular docking and molecular dynamics simulation approaches in the discrimination between active and inactive structurally related Mpro inhibitors.

Development of Reduced Peptide Bond Pseudopeptide Michael Acceptors for the Treatment of Human African Trypanosomiasis.

Human African Trypanosomiasis (HAT) is an endemic protozoan disease widespread in the sub-Saharan region that is caused by T. b. gambiense and T. b. rhodesiense. The development of molecules targeting rhodesain, the main cysteine protease of T. b. rhodesiense, has led to a panel of inhibitors endowed with micro/sub-micromolar activity towards the protozoa. However, whilst impressive binding affinity against rhodesain has been observed, the limited selectivity towards the target still remains a hard challenge for the development of antitrypanosomal agents. In this paper, we report the synthesis, biological evaluation, as well as docking studies of a series of reduced peptide bond pseudopeptide Michael acceptors (SPR10-SPR19) as potential anti-HAT agents. The new molecules show Ki values in the low-micro/sub-micromolar range against rhodesain, coupled with k2nd values between 1314 and 6950 M-1 min-1. With a few exceptions, an appreciable selectivity over human cathepsin L was observed. In in vitro assays against T. b. brucei cultures, SPR16 and SPR18 exhibited single-digit micromolar activity against the protozoa, comparable to those reported for very potent rhodesain inhibitors, while no significant cytotoxicity up to 70 µM towards mammalian cells was observed. The discrepancy between rhodesain inhibition and the antitrypanosomal effect could suggest additional mechanisms of action. The biological characterization of peptide inhibitor SPR34 highlights the essential role played by the reduced bond for the antitrypanosomal effect. Overall, this series of molecules could represent the starting point for further investigations of reduced peptide bond-containing analogs as potential anti-HAT agents.

Pan-3C Protease Inhibitor Rupintrivir Binds SARS-CoV-2 Main Protease in a Unique Binding Mode.

Rupintrivir targets the 3C cysteine proteases of the picornaviridae family, which includes rhinoviruses and enteroviruses that cause a range of human diseases. Despite being a pan-3C protease inhibitor, rupintrivir activity is extremely weak against the homologous 3C-like protease of SARS-CoV-2. In this study, the crystal structures of rupintrivir were determined bound to enterovirus 68 (EV68) 3C protease and the 3C-like main protease (Mpro) from SARS-CoV-2. While the EV68 3C protease-rupintrivir structure was similar to previously determined complexes with other picornavirus 3C proteases, rupintrivir bound in a unique conformation to the active site of SARS-CoV-2 Mpro splitting the catalytic cysteine and histidine residues. This bifurcation of the catalytic dyad may provide a novel approach for inhibiting cysteine proteases.

Molecular docking of potential SARS-CoV-2 papain-like protease inhibitors.

SARS-CoV-2 papain-like protease is considered as an important potential target for anti-SARS-CoV-2 drug discovery due to its crucial roles in viral spread and innate immunity. Here, we have utilized an in silico molecular docking approach to identify the possible inhibitors of the SARS-CoV-2 papain-like protease, by screening 21 antiviral, antifungal and anticancer compounds. Among them, Neobavaisoflavone has the highest binding energy for SARS-CoV-2 papain-like protease. These molecules could bind near the SARS-CoV-2 papain-like protease crucial catalytic triad, ubiquitination and ISGylation residues: Trp106, Asn109, Cys111, Met208, Lys232, Pro247, Tyr268, Gln269, His272, Asp286 and Thr301. Because blocking the papain-like protease is an important strategy in fighting against viruses, these compounds might be promising candidates for therapeutic intervention against COVID-19.

Equilibrium between tri- and tetra-coordinate chalcogenuranes is critical for cysteine protease inhibition.

There have been significant advances in the biological use of hypervalent selenium and tellurium compounds as cysteine protease inhibitors. However, the full understanding of their reaction mechanisms for and cysteine proteases inhibition is still elusive. Kinetic studies suggest an irreversible inhibition mechanism, which was explained by forming a covalent bond between the enzyme sulfhydryl group and the chalcogen atom at its hypervalent state (+4). In this work, we performed a theoretical investigation using density functional theory to propose the active inhibitor form in an aqueous solution. To this end, we investigated chloride ligand exchange reactions by oxygen and sulfur nucleophiles on hypervalent selenium and tellurium compounds. All tetra- and tri-coordinated chalcogen compounds and distinct protonation states of the nucleophiles were considered, totaling 34 unique species, 7 nucleophiles, and 155 free energies reactions. We discovered that chloride is easily replaced by a nonprotonated nucleophile (SH- or OH- ) in R2 SeCl2 . We also found that tri-coordinate species are more stable than their tetra-coordinate counterparts, with selenoxide (R2 SeO) protonation being strongly exergonic in acid pH. The thermodynamic and kinetic results suggest that the protonated selenoxide (R2 SeOH+ ) is the most probable active chemical species in biological media. The computed energetic profiles paint a possible picture for selenuranes activity, with successive exergonic steps leading to a covalent inhibition of thiol-dependent enzymes, like cysteine proteases. A second pathway has also been uncovered, with a direct reaction to chalcogenonium cation (R2 SeCl+ ) as the inhibition step. Tellurium compounds showed similar trends but formed telluroxide in a pH-independent fashion.

A head-to-head comparison of the inhibitory activities of 15 peptidomimetic SARS-CoV-2 3CLpro inhibitors.

The COVID-19 pandemic caused by SARS-CoV-2 has created an unprecedented global health emergency. As of July 2021, only three antiviral therapies have been approved by the FDA for treating infected patients, highlighting the urgent need for more antiviral drugs. The SARS-CoV-2 3CL protease (3CLpro) is deemed an attractive drug target due to its essential role in viral polyprotein processing and pathogenesis. Indeed, a number of peptidomimetic 3CLpro inhibitors armed with electrophilic warheads have been reported by various research groups that can potentially be developed for treating COVID-19. However, it is currently impossible to compare their relative potencies due to the different assays employed. To solve this, we conducted a head-to-head comparison of fifteen reported peptidomimetic inhibitors in a standard FRET-based SARS-CoV-2 3CLpro inhibition assay to compare and identify potent inhibitors for development. Inhibitor design and the suitability of various warheads are also discussed.

Structure-activity relationship studies of 3-substituted pyrazoles as novel allosteric inhibitors of MALT1 protease.

We report the discovery of a novel series of 1,5-bisphenylpyrazoles as potent MALT1 inhibitors. Structure-activity relationship exploration of a hit compound led to a potent MALT1 inhibitor. Compound 33 showed strong activity against MALT1 (IC50: 0.49 µM), potent cellular activity (NF-κB inhibition and inhibition of IL2 production), and high selectivity against caspase-3, -8, and -9. The results of a kinetics study suggest that compound 33 is a non-competitive inhibitor of MALT1 protein.

SARS-CoV-2 main protease inhibition by compounds isolated from Luffa cylindrica using molecular docking.

In this study, chemical investigation of methanol extract of the air-dried fruits of Luffa cylindrica led to the identification of a new δ-valerolactone (1), along with sixteen known compounds (2-17). Their chemical structures including the absolute configuration were elucidated by extensive spectroscopic analysis and electronic circular dichroism analysis, as well as by comparison with those reported in the literature. For the first time in literature, we have examined the binding potential of the isolated compounds to highly conserved protein, Mpro of SARS-CoV-2 using the molecular docking technique. We found that the isolated saponins (14-17) bind to the substrate-binding pocket of SARS-CoV-2 Mpro with docking energy scores of -7.13, -7.29, -7.47, and -7.54 kcal.mol-1, respectively, along with binding abilities equivalent to an already claimed N3 protease inhibitor (-7.51 kcal.mol-1).