Cysteine protease inhibition by nitrile-based inhibitors: a computational study

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Targeting SmCB1: Perspectives and Insights to Design Antischistosomal Drugs

Neglected tropical diseases (NTDs) are prevalent in tropical and subtropical countries, and schistosomiasis is among the most relevant diseases worldwide. In addition, one of the two biggest problems in developing drugs against this disease is related to drug resistance, which promotes the demand to develop new drug candidates for this purpose. Thus, one of the drug targets most explored, Schistosoma mansoni Cathepsin B1 (SmCB1 or Sm31), provides new opportunities in drug development due to its essential functions for the parasite's survival. In this way, here, the latest developments in drug design studies targeting SmCB1 were approached, focusing on the most promising analogs of nitrile, vinyl sulphones, and peptidomimetics. Thus, it was shown that despite being a disease known since ancient times, it remains prevalent throughout the world, with high mortality rates. The therapeutic arsenal of antischistosomal drugs (ASD) consists only of praziquantel, which is widely used for this purpose and has several advantages, such as efficacy and safety. However, it has limitations, such as the impossibility of acting on the immature worm and exploring new targets to overcome these limitations. SmCB1 shows its potential as a cysteine protease with a catalytic triad consisting of Cys100, His270, and Asn290. Thus, design studies of new inhibitors focus on their catalytic mechanism for designing new analogs. In fact, nitrile and sulfonamide analogs show the most significant potential in drug development, showing that these chemical groups can be better exploited in drug discovery against schistosomiasis. We hope this manuscript guides the authors in searching for promising new antischistosomal drugs.

Discrete Dynamics of Warhead Modulation on Covalent Inhibition of Oxyr: A QM/MM Study

The bacterial transcriptional factor OxyR, a peroxide sensor conserved in bacterial virulence pathways, has the capability to exhibit exceptional reactivity toward hydrogen peroxide (H₂O₂). H₂O₂ is essential for oxidizing cysteine thiolates to maintain cellular redox homeostasis and is dispensable for bacterial growth that can potentially mitigate drug resistance, thus underlining OxyR as a valuable target. We employ quantum mechanics/molecular mechanics (QM/MM) umbrella sampling (US) simulations at the DFTB3/MM level of theory and propose a reaction mechanism with four potential covalent inhibitors. The potential of mean force reveals the direct role of intrinsic reactivity of inhibitors, for instance, benzothiophenes and modified experimental inhibitors with methyl oxo-enoate warhead-activated carbonyl samples in the first step of reaction, which shed light on the significance of proton transfer indispensable for full inhibition, whereas the nitrile inhibitor undergoes a stepwise mechanism with a small proton-transfer energy barrier and lower imaginary frequencies that materialize instantly after nucleophilic attack. To unveil the molecular determinants of respective binding affinities, transition states along the reaction path are optimized and characterized with B3LYP 6-31+G(d,p). Furthermore, the post-simulation analysis indicates the catalytic triad (His130/Cys199/Thr129), thermodynamically favored for inhibition, which restricts water molecules from acting as the potential source of protonation/deprotonation. This study thus serves as a preamble to add variation in the proposed structures and unveils the impact of functional groups lying in warheads that modulate the kinetics of proton transfer, which will certainly aid to design more selective and efficient irreversible inhibitors of OxyR.

An update on the discovery and development of reversible covalent inhibitors

Small molecule drugs that covalently bind irreversibly to their target proteins have several advantages over conventional reversible inhibitors. They include increased duration of action, less-frequent drug dosing, reduced pharmacokinetic sensitivity, and the potential to target intractable shallow binding sites. Despite these advantages, the key challenges of irreversible covalent drugs are their potential for off-target toxicities and immunogenicity risks. Incorporating reversibility into covalent drugs would lead to less off-target toxicity by forming reversible adducts with off-target proteins and thus reducing the risk of idiosyncratic toxicities caused by the permanent modification of proteins, which leads to higher levels of potential haptens. Herein, we systematically review electrophilic warheads employed during the development of reversible covalent drugs. We hope the structural insights of electrophilic warheads would provide helpful information to medicinal chemists and aid in designing covalent drugs with better on-target selectivity and improved safety.

Graphical Abstract

Dipeptide-Derived Alkynes as Potent and Selective Irreversible Inhibitors of Cysteine Cathepsins

The potential of designing irreversible alkyne-based inhibitors of cysteine cathepsins by isoelectronic replacement in reversibly acting potent peptide nitriles was explored. The synthesis of the dipeptide alkynes was developed with special emphasis on stereochemically homogeneous products obtained in the Gilbert-Seyferth homologation for C≡C bond formation. Twenty-three dipeptide alkynes and 12 analogous nitriles were synthesized and investigated for their inhibition of cathepsins B, L, S, and K. Numerous combinations of residues at positions P1 and P2 as well as terminal acyl groups allowed for the derivation of extensive structure-activity relationships, which were rationalized by computational covalent docking for selected examples. The determined inactivation constants of the alkynes at the target enzymes span a range of >3 orders of magnitude (3-10 133 M^-1 s^-1). Notably, the selectivity profiles of alkynes do not necessarily reflect those of the nitriles. Inhibitory activity at the cellular level was demonstrated for selected compounds.

The pyridazine heterocycle in molecular recognition and drug discovery

The pyridazine ring is endowed with unique physicochemical properties, characterized by weak basicity, a high dipole moment that subtends π-π stacking interactions and robust, dual hydrogen-bonding capacity that can be of importance in drug-target interactions. These properties contribute to unique applications in molecular recognition while the inherent polarity, low cytochrome P450 inhibitory effects and potential to reduce interaction of a molecule with the cardiac hERG potassium channel add additional value in drug discovery and development. The recent approvals of the gonadotropin-releasing hormone receptor antagonist relugolix (24) and the allosteric tyrosine kinase 2 inhibitor deucravacitinib (25) represent the first examples of FDA-approved drugs that incorporate a pyridazine ring. In this review, the properties of the pyridazine ring are summarized in comparison to the other azines and its potential in drug discovery is illustrated through vignettes that explore applications that take advantage of the inherent physicochemical properties as an approach to solving challenges associated with candidate optimization.

Graphical Abstract

Targeting BAP1 with small compound inhibitor for colon cancer treatment

BRCA1-associated protein-1 (BAP1) is a ubiquitin C-terminal hydrolase domain-containing deubiquitinase. The gene encoding BAP1 is mutated in various human cancers, including mesothelioma, uveal melanoma and renal cell carcinoma. BAP1 plays roles in many cancer-related cellular functions, including cell proliferation, cell death, and nuclear processes crucial for genome stability, such as DNA repair and replication. While these findings suggest that BAP1 functions as a tumor suppressor, recent data also suggest that BAP1 might play tumor-promoting roles in certain cancers, such as breast cancer and hematopoietic malignancies. Here, we show that BAP1 is upregulated in colon cancer cells and tissues and that BAP1 depletion reduces colon cancer cell proliferation and tumor growth. BAP1 contributes to colon cancer cell proliferation by accelerating DNA replication and suppressing replication stress and concomitant apoptosis. A recently identified BAP1 inhibitor, TG2-179-1, which seems to covalently bind to the active site of BAP1, exhibits potent cytotoxic activity against colon cancer cells, with half-maximal inhibitory concentrations of less than 10 μM, and inhibits colon tumor growth. TG2-179-1 exerts cytotoxic activity by targeting BAP1, leading to defective replication and increased apoptosis. This work therefore shows that BAP1 acts oncogenically in colon cancer and is a potential therapeutic target for this cancer. Our work also suggests that TG2-179-1 can be developed as a potential therapeutic agent for colon cancer.

Mechanistic investigation of SARS-CoV-2 main protease to accelerate design of covalent inhibitors

Targeted covalent inhibition represents one possible strategy to block the function of SARS-CoV-2 Main Protease (M^PRO), an enzyme that plays a critical role in the replication of the novel SARS-CoV-2. Toward the design of covalent inhibitors, we built a covalent inhibitor dataset using deep learning models followed by high throughput virtual screening of these candidates against M^PRO. Two top-ranking inhibitors were selected for mechanistic investigations-one with an activated ester warhead that has a piperazine core and the other with an acrylamide warhead. Specifically, we performed a detailed analysis of the free energetics of covalent inhibition by hybrid quantum mechanics/molecular mechanics simulations. Cleavage of a fragment of the non-structured protein (NSP) from the SARS-CoV-2 genome was also simulated for reference. Simulations show that both candidates form more stable enzyme-inhibitor (E-I) complexes than the chosen NSP. It was found that both the NSP fragment and the activated ester inhibitor react with CYS145 of M^PRO in a concerted manner, whereas the acrylamide inhibitor follows a stepwise mechanism. Most importantly, the reversible reaction and the subsequent hydrolysis reaction from E-I complexes are less probable when compared to the reactions with an NSP fragment, showing promise for these candidates to be the base for efficient M^PRO inhibitors.

The Se-S Bond Formation in the Covalent Inhibition Mechanism of SARS-CoV-2 Main Protease by Ebselen-like Inhibitors: A Computational Study

The inhibition mechanism of the main protease (Mpro) of SARS-CoV-2 by ebselen (EBS) and its analog with a hydroxyl group at position 2 of the benzisoselenazol-3(2H)-one ring (EBS-OH) was studied by using a density functional level of theory. Preliminary molecular dynamics simulations on the apo form of Mpro were performed taking into account both the hydrogen donor and acceptor natures of the Nδ and Nε of His41, a member of the catalytic dyad. The potential energy surfaces for the formation of the Se-S covalent bond mediated by EBS and EBS-OH on Mpro are discussed in detail. The EBS-OH shows a distinctive behavior with respect to EBS in the formation of the noncovalent complex. Due to the presence of canonical H-bonds and noncanonical ones involving less electronegative atoms, such as sulfur and selenium, the influence on the energy barriers and reaction energy of the Minnesota hybrid meta-GGA functionals M06, M06-2X and M08HX, and the more recent range-separated hybrid functional wB97X were also considered. The knowledge of the inhibition mechanism of Mpro by the small protease inhibitors EBS or EBS-OH can enlarge the possibilities for designing more potent and selective inhibitor-based drugs to be used in combination with other antiviral therapies.

A Nut for Every Bolt: Subunit-Selective Inhibitors of the Immunoproteasome and Their Therapeutic Potential

At the heart of the ubiquitin-proteasome system, the 20S proteasome core particle (CP) breaks down the majority of intracellular proteins tagged for destruction. Thereby, the CP controls many cellular processes including cell cycle progression and cell signalling. Inhibitors of the CP can suppress these essential biological pathways, resulting in cytotoxicity, an effect that is beneficial for the treatment of certain blood cancer patients. During the last decade, several preclinical studies demonstrated that selective inhibition of the immunoproteasome (iCP), one of several CP variants in mammals, suppresses autoimmune diseases without inducing toxic side effects. These promising findings led to the identification of natural and synthetic iCP inhibitors with distinct chemical structures, varying potency and subunit selectivity. This review presents the most prominent iCP inhibitors with respect to possible scientific and medicinal applications, and discloses recent trends towards pan-immunoproteasome reactive inhibitors that cumulated in phase II clinical trials of the lead compound KZR-616 for chronic inflammations.

An insight into the interaction between α-ketoamide- based inhibitor and coronavirus main protease: A detailed in silico study

The search for therapeutic drugs that can neutralize the effects of COVID-2019 (SARS-CoV-2) infection is the main focus of current research. The coronavirus main protease (M_pro) is an attractive target for anti-coronavirus drug design. Further, α-ketoamide is proved to be very effective as a reversible covalent-inhibitor against cysteine proteases. Herein, we report on the non-covalent to the covalent adduct formation mechanism of α-ketoamide-based inhibitor with the enzyme active site amino acids by QM/SQM model (QM = quantum mechanical, SQM = semi-empirical QM). To uncover the mechanism, we focused on two approaches: a concerted and a stepwise fashion. The concerted pathway proceeds via deprotonation of the thiol of cysteine (here, Cys₁₄₅ SγH) and simultaneous reversible nucleophilic attack of sulfur onto the α-ketoamide warhead. In this work, we propose three plausible concerted pathways. On the contrary, in a traditional two-stage pathway, the first step is proton transfer from Cys₁₄₅ SγH to His₄₁ Nδ forming an ion pair, and consecutively, in the second step, the thiolate ion attacks the α-keto group to form a thiohemiketal. In this reaction, we find that the stability of the tetrahedral intermediate oxyanion/hydroxyl group plays an important role. Moreover, as the α-keto group has two faces Si or Re for the nucleophilic attack, we considered both possibilities of attack leading to S- and R-thiohemiketal. We computed the structural, electronic, and energetic parameters of all stationary points including transition states via ONIOM and pure DFT method. Additionally, to characterize covalent, weak noncovalent interaction (NCI) and hydrogen-bonds, we applied NCI-reduced density gradient (NCI-RDG) methods along with Bader's Quantum Theory of Atoms-in-Molecules (QTAIM) and natural bonding orbital (NBO) analysis.

Azanitrile Inhibitors of the SmCB1 Protease Target Are Lethal to Schistosoma mansoni: Structural and Mechanistic Insights into Chemotype Reactivity

Azapeptide nitriles are postulated to reversibly covalently react with the active-site cysteine residue of cysteine proteases and form isothiosemicarbazide adducts. We investigated the interaction of azadipeptide nitriles with the cathepsin B1 drug target (SmCB1) from Schistosoma mansoni, a pathogen that causes the global neglected disease schistosomiasis. Azadipeptide nitriles were superior inhibitors of SmCB1 over their parent carba analogs. We determined the crystal structure of SmCB1 in complex with an azadipeptide nitrile and analyzed the reaction mechanism using quantum chemical calculations. The data demonstrate that azadipeptide nitriles, in contrast to their carba counterparts, undergo a change from E- to Z-configuration upon binding, which gives rise to a highly favorable energy profile of noncovalent and covalent complex formation. Finally, azadipeptide nitriles were considerably more lethal than their carba analogs against the schistosome pathogen in culture, supporting the further development of this chemotype as a treatment for schistosomiasis.

Exploring the Catalytic Reaction of Cysteine Proteases

Cysteine proteases play a major role in many life processes and are the target of key drugs. The reaction mechanism of these enzymes is a complex process, which involves several steps that are divided into two main groups: acylation and deacylation. In this work, we studied the energy profile for the acylation and a part of the deacylation reaction of three different enzymes, cruzain, papain, and the Q19A-mutated papain with the benzyloxycarbonyl-phenylalanylarginine-4-methylcoumaryl-7-amide (CBZ-FR-AMC) substrate. The calculations were performed using the EVB and PDLD/S-LRA methods. The overall agreement between the calculated and observed results is encouraging and indicates that we captured the correct reaction mechanism. Finally, our finding indicates that the minimum of the reaction profile, between the acylation and deacylation steps, should provide an excellent state for the binding of covalent inhibitors.

Advances in Sustainable Catalysis: A Computational Perspective

The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.

Experimental study and computational modelling of cruzain cysteine protease inhibition by dipeptidyl nitriles

Chagas disease affects millions of people in Latin America. This disease is caused by the protozoan parasite Trypanossoma cruzi. The cysteine protease cruzain is a key enzyme for the survival and propagation of this parasite lifecycle. Nitrile-based inhibitors are efficient inhibitors of cruzain that bind by forming a covalent bond with this enzyme. Here, three nitrile-based inhibitors dubbed Neq0409, Neq0410 and Neq0570 were synthesized, and the thermodynamic profile of the bimolecular interaction with cruzain was determined using isothermal titration calorimetry (ITC). The result suggests the inhibition process is enthalpy driven, with a detrimental contribution of entropy. In addition, we have used hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) and Molecular Dynamics (MD) simulations to investigate the reaction mechanism of reversible covalent modification of cruzain by Neq0409, Neq0410 and Neq0570. The computed free energy profile shows that the nucleophilic attack of Cys25 on the carbon C1 of inhibitiors and the proton transfer from His162 to N1 of the dipeptidyl nitrile inhibitor take place in a single step. The calculated free energy of the inhibiton reaction is in agreement with covalent experimental binding. Altogether, the results reported here suggests that nitrile-based inhibitors are good candidates for the development of reversible covalent inhibitors of cruzain and other cysteine proteases.

Ca2+ binding induced sequential allosteric activation of sortase A: An example for ion-triggered conformational selection

The allosteric activation of the intrinsically disordered enzyme Staphylococcus aureus sortase A is initiated via binding of a Ca²⁺ ion. Although Ca²⁺ binding was shown to initiate structural changes inducing disorder-to-order transitions, the details of the allosteric activation mechanism remain elusive. We performed long-term molecular dynamics simulations of sortase A without (3 simulations of 1.6 μs) and with bound Ca²⁺ (simulations of 1.6 μs, 1.8 μs, and 2.5 μs). Our results show that Ca²⁺ binding causes not only ordering of the disordered β6/β7 loop of the protein, but also modulates hinge motions in the dynamic β7/β8 loop, which is important for the catalytic activity of the enzyme. Cation binding triggers signal transmission from the Ca²⁺ binding site to the dynamic β7/β8 loop via the repetitive folding/unfolding of short helical stretches of the disordered β6/β7 loop. These correlated structural rearrangements lead to several distinct conformational states of the binding groove, which show optimal binding features for the sorting signal motif and feature binding energies up to 20 kcal/mol more favorable than observed for the sortase A without Ca²⁺. The presented results indicate a highly correlated, conformational selection-based activation mechanism of the enzyme triggered by cation binding. They also demonstrate the importance of the dynamics of intrinsically disordered regions for allosteric regulation.

Density functional theory based probe of the affinity interaction of saccharide ligands with extra-cellular sialic acid residues

Changes in glycosylation pattern leads to malignant transformations among the cells. In combination with upregulated actions of sialyltransferases, it ultimately leads to differential expression of sialic acid (SA) at cell surface. Given its negative charge and localization to extracellular domain, SA has been exploited for the development of targeted theranostics using approaches, such as, cationization and appending recognition saccharides on carrier surface. In this study, we have performed quantum mechanical calculations based on density functional theory (DFT) to study the interaction of saccharides with extracellular SA. Gradient-corrected DFT with the three parameter function (B3) was utilized for the calculation of Lee-Yang-Parr (LYP) correlation function. Atomic charge, vibrational frequencies and energy of the optimized structures were calculated through B3LYP. Our calculations demonstrate a stronger galactose-sialic acid interaction at tumour-relevant low pH and hyperthermic condition. These results support the application of pH responsive delivery vehicles and targeted hyperthermic chemotherapy for eradicating solid tumour deposits. These studies, conducted a priori, can guide the formulation scientists over appropriate choice of ligands and their applications in the design of 'smart' theranostic tools.

Theoretical studies on the interaction between the nitrile-based inhibitors and the catalytic triad of Cathepsin K

Computational studies on the interaction of novel inhibitor compounds with the Cathepsin K protease have been performed to study the inhibition properties of the inhibitor compounds. The quantum chemical calculations have been performed to analyze the molecular geometries, structural stability, reactivity, nature of interaction, and the charge transfer properties using B3LYP level of theory by implementing 6-311g(d,p) basis set. The calculated C-S and N-H…N bond lengths of the inhibitor-triad complexes are found to agree well with the previous literature results. The chemical reactivity of the inhibitors and catalytic triad are analyzed through frontier molecular orbital analysis and found that the inhibitors are subjected to nucleophilic attack by the catalytic triad. The nature of inhibition of the inhibitor compounds is examined using the quantum theory of Atoms in Molecules analysis and found to be partially covalent. The NBO stabilization energy for the Cys - inhibitor are found to be most stable than the other interactions. The molecular dynamic simulations were performed to study the influence of dynamic of the active site on the QM results. The many body decomposition interaction energy calculated for the final results of MD simulation reveals that the dynamic of the active site induces significant changes in the interaction energy and occupancy of H-bonds plays a major role in the stabilizing the active site inhibitor interactions. The present study reveals that the inhibitor compounds can inhibit the proteolytic activity of the proteases on binding with the catalytic active site.

Chinese herb medicine against Sortase A catalyzed transformations, a key role in gram-positive bacterial infection progress

Many Gram-positive bacteria can anchor their surface proteins to the cell wall peptidoglycan covalently by a common mechanism with Sortase A (SrtA), thus escaping from the host's identification of immune cells. SrtA can complete this anchoring process by cleaving LPXTG motif conserved among these surface proteins and thus these proteins anchor on the cell wall. Moreover, those SrtA mutants lose this capability to anchor these relative proteins, with these bacteria no longer infectious. Therefore, SrtA inhibitors can be promising anti-infective agents to cure bacterial infections. Chinese herb medicines (CHMs) (chosen from Science Citation Index) have exhibited inhibition on SrtA of Gram-positive pathogens irreversibly or reversibly. In general, CHMs are likely to have important long-term impact as new antibacterial compounds and sought after by academia and the pharmaceutical industry. This review mainly focuses on SrtA inhibitors from CHMs and the potential inhibiting mechanism related to chemical structures of compounds in CHMs.

Methyl-methoxylpyrrolinone and flavinium nucleus binding signatures on falcipain-2 active site

Following the increasing reports of human toxicity and plasmodium resistance to artemisinin and its derivatives, falcipain-2 (FP-2) is now emerging as the choice antimalarial drug target. Coincidentally, FP-2 is the in vivo target of naturally occurring, therapeutically safe flavonoids (stenopalustroside, myricetin, and fisetin) and symplostatin (symplostatin 4) compounds known to exhibit potent in vitro and in vivo antiplasmodial actions. Here, the structural bases for their inhibitory actions have been studied using molecular dynamics simulation. Myricetin and fisetin act as proton transfer tunnel breakers by inserting between His174 and Cys42, which are key active site residues of FP-2, stenopalustroside inhibits the polarization of His174 by Asn173; a major preparatory step for Cys42/His174 proton transfer process. The roles of flavonoids are favored by T-shaped pi-pi interactions with His174. Symplostatin 4 inserts its methyl-methoxylpyrrolinone moiety into the active site where its proton acceptor function prepares Cys42 for nucleophilic attack on the Michael α,β-unsaturated bonds on its 4(S)-amino-2(E)-pentenoate moiety. Further analyses of the structures identified a unique bridge formed on FP-2 active site groove by stenopalustroside and symplostatin 4 during interaction with the sub-site I of FP-2, whereas fisetin preferentially interacts with sub-site II and myricetin interacts with sub-site III residues. Ultimately, symplostatin-4, myricetin, and fisetin were better than stenopalustroside at trapping FP-2 in its inactive state as revealed by comparative RSMD plots with X-ray structures of FP-2 co-crystallized with inhibitors. Comparative estimates of free energy of binding using the Molecular Mechanics-Poisson Boltzmann Surface Area (MMPBSA) method suggested that His174 protonation may further enhance stenopalustroside-FP-2 interaction. The unique binding signatures of the ligands within the FP-2 active site groove and its sub-sites may explain the subtle differences in their IC50 values and their mechanism of inhibition.