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

Quantum mechanical-based strategies in drug discovery: Finding the pace to new challenges in drug design

The expansion of the chemical space to tangible libraries containing billions of synthesizable molecules opens exciting opportunities for drug discovery, but also challenges the power of computer-aided drug design to prioritize the best candidates. This directly hits quantum mechanics (QM) methods, which provide chemically accurate properties, but subject to small-sized systems. Preserving accuracy while optimizing the computational cost is at the heart of many efforts to develop high-quality, efficient QM-based strategies, reflected in refined algorithms and computational approaches. The design of QM-tailored physics-based force fields and the coupling of QM with machine learning, in conjunction with the computing performance of supercomputing resources, will enhance the ability to use these methods in drug discovery. The challenge is formidable, but we will undoubtedly see impressive advances that will define a new era.

Investigation of steric hindrance effect on the interactions between four alkaloids and HSA by isothermal titration calorimetry and molecular docking

The binding of four alkaloids with human serum albumin (HSA) was investigated by isothermal titration calorimetry (ITC), spectroscopy and molecular docking techniques. The findings demonstrated that theophylline or caffeine can bind to HAS, respectively. The number of binding sites and binding constants are obtained. The binding mode is a static quenching process. The effects of steric hindrance, temperature, salt concentration and buffer solution on the binding indicated that theophylline and HSA have higher binding affinity than caffeine. The fluorescence and ITC results showed that the interaction between HSA and theophylline or caffeine is an entropy-driven spontaneous exothermic process. The hydrophobic force was the primary driving factor. The experimental results were consistent with the molecular docking data. Based on the molecular structures of the four alkaloids, steric hindrance might be a major factor in the binding between HSA and these four alkaloids. This study elucidates the mechanism of interactions between four alkaloids and HSA.

Discovery of Novel Inhibitors of Cruzain Cysteine Protease of Trypanosoma cruzi

Chagas disease (CD) is a parasitic disease endemic in several developing countries. According to the World Health Organization, approximately 6-8 million people worldwide are inflicted by CD. The scarcity of new drugs, mainly for the chronic phase, is the main reason for treatment limitation in CD. Therefore, there is an urgent need to discover new targets for which new therapeutical agents could be developed. Cruzain cysteine protease (CCP) is a promising alternative because this enzyme exhibits pleiotropic effects by acting as a virulence factor, modulating host immune cells, and interacting with host cells. This systematic review was conducted to discover new compounds that act as cruzain inhibitors, and their effects in vitro were studied through enzymatic assays and molecular docking. Additionally, the advances and perspectives of these inhibitors are discussed. These findings are expected to contribute to medicinal chemistry in view of the design of new, safe, and efficacious inhibitors against Trypanosoma cruzi CCP detected in the last decade (2013-2022) to provide scaffolds for further optimization, aiming toward the discovery of new drugs.

The Power of Molecular Dynamics Simulations and Their Applications to Discover Cysteine Protease Inhibitors

A large family of enzymes with the function of hydrolyzing peptide bonds, called peptidases or cysteine proteases (CPs), are divided into three categories according to the peptide chain involved. CPs catalyze the hydrolysis of amide, ester, thiol ester, and thioester peptide bonds. They can be divided into several groups, such as papain-like (CA), viral chymotrypsin-like CPs (CB), papainlike endopeptidases of RNA viruses (CC), legumain-type caspases (CD), and showing active residues of His, Glu/Asp, Gln, Cys (CE). The catalytic mechanism of CPs is the essential cysteine residue present in the active site. These mechanisms are often studied through computational methods that provide new information about the catalytic mechanism and identify inhibitors. The role of computational methods during drug design and development stages is increasing. Methods in Computer-Aided Drug Design (CADD) accelerate the discovery process, increase the chances of selecting more promising molecules for experimental studies, and can identify critical mechanisms involved in the pathophysiology and molecular pathways of action. Molecular dynamics (MD) simulations are essential in any drug discovery program due to their high capacity for simulating a physiological environment capable of unveiling significant inhibition mechanisms of new compounds against target proteins, especially CPs. Here, a brief approach will be shown on MD simulations and how the studies were applied to identify inhibitors or critical information against cysteine protease from several microorganisms, such as Trypanosoma cruzi (cruzain), Trypanosoma brucei (rhodesain), Plasmodium spp. (falcipain), and SARS-CoV-2 (Mpro). We hope the readers will gain new insights and use our study as a guide for potential compound identifications using MD simulations.

Landscape of In Silico Tools for Modeling Covalent Modification of Proteins: A Review on Computational Covalent Drug Discovery

Covalent drug discovery has been a challenging research area given the struggle of finding a sweet balance between selectivity and reactivity for these drugs, the lack of which often leads to off-target activities and hence undesirable side effects. However, there has been a resurgence in covalent drug design following the success of several covalent drugs such as boceprevir (2011), ibrutinib (2013), neratinib (2017), dacomitinib (2018), zanubrutinib (2019), and many others. Design of covalent drugs includes many crucial factors, where "evaluation of the binding affinity" and "a detailed mechanistic understanding on covalent inhibition" are at the top of the list. Well-defined experimental techniques are available to elucidate these factors; however, often they are expensive and/or time-consuming and hence not suitable for high throughput screens. Recent developments in in silico methods provide promise in this direction. In this report, we review a set of recent publications that focused on developing and/or implementing novel in silico techniques in "Computational Covalent Drug Discovery (CCDD)". We also discuss the advantages and disadvantages of these approaches along with what improvements are required to make it a great tool in medicinal chemistry in the near future.

Covalent Inhibitors for Neglected Diseases: An Exploration of Novel Therapeutic Options

Neglected diseases, primarily found in tropical regions of the world, present a significant challenge for impoverished populations. Currently, there are 20 diseases considered neglected, which greatly impact the health of affected populations and result in difficult-to-control social and economic consequences. Unfortunately, for the majority of these diseases, there are few or no drugs available for patient treatment, and the few drugs that do exist often lack adequate safety and efficacy. As a result, there is a pressing need to discover and design new drugs to address these neglected diseases. This requires the identification of different targets and interactions to be studied. In recent years, there has been a growing focus on studying enzyme covalent inhibitors as a potential treatment for neglected diseases. In this review, we will explore examples of how these inhibitors have been used to target Human African Trypanosomiasis, Chagas disease, and Malaria, highlighting some of the most promising results so far. Ultimately, this review aims to inspire medicinal chemists to pursue the development of new drug candidates for these neglected diseases, and to encourage greater investment in research in this area.

Multiscale Simulations of the Covalent Inhibition of the SARS-CoV-2 Main Protease: Four Compounds and Three Reaction Mechanisms

We report the results of computational modeling of the reactions of the SARS-CoV-2 main protease (MPro) with four potential covalent inhibitors. Two of them, carmofur and nirmatrelvir, have shown experimentally the ability to inhibit MPro. Two other compounds, X77A and X77C, were designed computationally in this work. They were derived from the structure of X77, a non-covalent inhibitor forming a tight surface complex with MPro. We modified the X77 structure by introducing warheads capable of reacting with the catalytic cysteine residue in the MPro active site. The reaction mechanisms of the four molecules with MPro were investigated by quantum mechanics/molecular mechanics (QM/MM) simulations. The results show that all four compounds form covalent adducts with the catalytic cysteine Cys 145 of MPro. From the chemical perspective, the reactions of these four molecules with MPro follow three distinct mechanisms. The reactions are initiated by a nucleophilic attack of the thiolate group of the deprotonated cysteine residue from the catalytic dyad Cys145-His41 of MPro. In the case of carmofur and X77A, the covalent binding of the thiolate to the ligand is accompanied by the formation of the fluoro-uracil leaving group. The reaction with X77C follows the nucleophilic aromatic substitution SNAr mechanism. The reaction of MPro with nirmatrelvir (which has a reactive nitrile group) leads to the formation of a covalent thioimidate adduct with the thiolate of the Cys145 residue in the enzyme active site. Our results contribute to the ongoing search for efficient inhibitors of the SARS-CoV-2 enzymes.

Impact of the Recognition Part of Dipeptidyl Nitroalkene Compounds on the Inhibition Mechanism of Cysteine Proteases Cruzain and Cathepsin L

Cysteine proteases (CPs) are an important class of enzymes, many of which are responsible for several human diseases. For instance, cruzain of protozoan parasite Trypanosoma cruzi is responsible for the Chagas disease, while the role of human cathepsin L is associated with some cancers or is a potential target for the treatment of COVID-19. However, despite paramount work carried out during the past years, the compounds that have been proposed so far show limited inhibitory action against these enzymes. We present a study of proposed covalent inhibitors of these two CPs, cruzain and cathepsin L, based on the design, synthesis, kinetic measurements, and QM/MM computational simulations on dipeptidyl nitroalkene compounds. The experimentally determined inhibition data, together with the analysis and the predicted inhibition constants derived from the free energy landscape of the full inhibition process, allowed describing the impact of the recognition part of these compounds and, in particular, the modifications on the P2 site. The designed compounds and, in particular, the one with a bulky group (Trp) at the P2 site show promising in vitro inhibition activities against cruzain and cathepsin L for use as a starting lead compound in the development of drugs with medical applications for the treatment of human diseases and future designs.

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 K_i was estimated to be 36.3 μM and K_i* 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.

Nitriles: an attractive approach to the development of covalent inhibitors

Nitriles have broad applications in medicinal chemistry, with more than 60 small molecule drugs on the market containing the cyano functional group. In addition to the well-known noncovalent interactions that nitriles can perform with macromolecular targets, they are also known to improve drug candidates' pharmacokinetic profiles. Moreover, the cyano group can be used as an electrophilic warhead to covalently bind an inhibitor to a target of interest, forming a covalent adduct, a strategy that can present benefits over noncovalent inhibitors. This approach has gained much notoriety in recent years, mainly with diabetes and COVID-19-approved drugs. Nevertheless, the application of nitriles in covalent ligands is not restricted to it being the reactive center, as it can also be employed to convert irreversible inhibitors into reversible ones, a promising strategy for kinase inhibition and protein degradation. In this review, we introduce and discuss the roles of the cyano group in covalent inhibitors, how to tune its reactivity and the possibility of achieving selectivity only by replacing the warhead. Finally, we provide an overview of nitrile-based covalent compounds in approved drugs and inhibitors recently described in the literature.

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.

2D-QSAR, docking, molecular dynamics, studies of PF-07321332 analogues to identify alternative inhibitors against 3CLpro enzyme in SARS-CoV disease

Given the results of the Pfizer-developed inhibitor PF-07321332 in the treatment of the SARS-Covid-19 epidemic, we aimed to identify potential alternatives to this compound by utilizing various methods; we developed 2 D-QSAR models to predict the therapeutic activity of 78 analogues of PF-07321332, three statistical learning techniques including (MLP-ANN), (SVR), and (MLR) were exploited. Various validation approaches were applied to the three models developed following the use of five most relevant descriptors. The study of the characteristics of these descriptors proved that the inhibitory activity of PF-07321332 analogues is specifically affected by the structure of the molecule, its polarizability, and by the hydrogen bonds. The best model, named MLP-ANN (with a 5-3-1 architecture), was selected on the basis of the following statistical parameters: r2 = 0.922, Q2 = 0.921. In addition, we performed a molecular docking and a molecular dynamics analysis of these compounds. The obtained results confirm that compound 8 can be a good alternative for compound PF-07321332.Communicated by Ramaswamy H. Sarma.

Combinatorial Effect of Dietary Oregano Extracts and 3,4,5-Trihydroxy Benzoic Acid on Growth Performance and Elimination of Coccidiosis in Broiler Chickens

We aimed to compare the combinatorial effect of 3,4,5-trihydroxybenzoic acid (THB) and oregano extracts (OE) with THB alone on the growth performance and elimination of deleterious effects in coccidiosis-infected broilers. A total of 210 one-day-old broilers were randomly assigned to one of five dietary treatments, with six replicates each, for 35 days. Dietary treatments were: 1) non-challenged, non-treated (NC); 2) challenged, non-treated (PC); 3) PC+ Salinomycin (0.05 g/kg; AB); 4) PC+THB (0.1 g/kg; THB); and 5) PC+THB+OE (0.1 g/kg; COM). On day 14, all groups except for NC were challenged with a 10-fold dose of Livacox® T anticoccidial vaccine to induce mild coccidiosis. All treatments significantly improved (P<0.05) body weight, average daily gain, and average daily feed intake, compared to PC, on days 21, 28, and 35. However, all treatments significantly reduced (P<0.05) the feed conversion ratio of PC by more than 14.60% on day 35, 11.76% during growing period, and 10.36% through the entire period. Broilers receiving anticoccidial treatments had 54.23% and 51.86% lower lesion scores (P<0.05) at 4 and 7 days post-infection, respectively, compared to PC. Additionally, the villus height of COM was significantly longer (P < 0.05) than that of THB. Although the molecular action of COM remains unclear, OE addition to THB reduced the shedding of oocysts better than THB alone (P<0.05, 9-11 days post-infection). Most importantly, COM effectively minimized the mortality of challenged birds from as high as 11.90% (PC) to 0%, a level similar to NC and AB, while THB maintained a mortality of 2.38%. In conclusion, the anticoccidial effect of THB can be enhanced by the addition of OE for better animal performance and the elimination of deleterious effects from coccidiosis-infected broilers for 35 days.

Mechanistic Modeling of Monoglyceride Lipase Covalent Modification Elucidates the Role of Leaving Group Expulsion and Discriminates Inhibitors with High and Low Potency

Inhibition of monoglyceride lipase (MGL), also known as monoacylglycerol lipase (MAGL), has emerged as a promising approach for treating neurological diseases. To gain useful insights in the design of agents with balanced potency and reactivity, we investigated the mechanism of MGL carbamoylation by the reference triazole urea SAR629 (IC₅₀ = 0.2 nM) and two recently described inhibitors featuring a pyrazole (IC₅₀ = 1800 nM) or a 4-cyanopyrazole (IC₅₀ = 8 nM) leaving group (LG), using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. Opposite to what was found for substrate 2-arachidonoyl-sn-glycerol (2-AG), covalent modification of MGL by azole ureas is controlled by LG expulsion. Simulations indicated that changes in the electronic structure of the LG greatly affect reaction energetics with triazole and 4-cyanopyrazole inhibitors following a more accessible carbamoylation path compared to the unsubstituted pyrazole derivative. The computational protocol provided reaction barriers able to discriminate between MGL inhibitors with different potencies. These results highlight how QM/MM simulations can contribute to elucidating structure–activity relationships and provide insights for the design of covalent inhibitors.

Covalent Reversible Inhibitors of Cysteine Proteases Containing the Nitrile Warhead: Recent Advancement in the Field of Viral and Parasitic Diseases

In the field of drug discovery, the nitrile group is well represented among drugs and biologically active compounds. It can form both non-covalent and covalent interactions with diverse biological targets, and it is amenable as an electrophilic warhead for covalent inhibition. The main advantage of the nitrile group as a warhead is mainly due to its milder electrophilic character relative to other more reactive groups (e.g., -CHO), reducing the possibility of unwanted reactions that would hinder the development of safe drugs, coupled to the ease of installation through different synthetic approaches. The covalent inhibition is a well-assessed design approach for serine, threonine, and cysteine protease inhibitors. The mechanism of hydrolysis of these enzymes involves the formation of a covalent acyl intermediate, and this mechanism can be exploited by introducing electrophilic warheads in order to mimic this covalent intermediate. Due to the relevant role played by the cysteine protease in the survival and replication of infective agents, spanning from viruses to protozoan parasites, we will review the most relevant and recent examples of protease inhibitors presenting a nitrile group that have been introduced to form or to facilitate the formation of a covalent bond with the catalytic cysteine active site residue.

Characterization of the non-covalent interaction between the PF-07321332 inhibitor and the SARS-CoV-2 main protease

We have studied the non-covalent interaction between PF-07321332 and SARS-CoV-2 main protease at the atomic level using a computational approach based on extensive molecular dynamics simulations with explicit solvent. PF-07321332, whose chemical structure has been recently disclosed, is a promising oral antiviral clinical candidate with well-established anti-SARS-CoV-2 activity in vitro. The drug, currently in phase III clinical trials in combination with ritonavir, relies on the electrophilic attack of a nitrile warhead to the catalytic cysteine of the protease. Nonbonded interaction between the inhibitor and the residues of the binding pocket, as well as with water molecules on the protein surface, have been characterized using two different force fields and the two possible protonation states of the main protease catalytic dyad HIS41-CYS145. When the catalytic dyad is in the neutral state, the non-covalent binding is likely to be stronger. Molecular dynamics simulations seems to lend support for an inhibitory mechanism in two steps: a first non-covalent addition with the dyad in neutral form and then the formation of the thiolate-imidazolium ion pair and the ligand relocation for finalising the electrophilic attack.

Predicting the Relative Binding Affinity for Reversible Covalent Inhibitors by Free Energy Perturbation Calculations

Covalent inhibitors are assuming central importance in drug discovery projects, especially in this pandemic scenario. Many research groups have focused their attention on inhibiting viral proteases or human proteases such as cathepsin L (hCatL). The inhibition of these critical enzymes may impair viral replication. However, molecular modeling of covalent ligands is challenging since covalent and noncovalent ligand-bound states must be considered in the binding process. In this work, we evaluated the suitability of free energy perturbation (FEP) calculations as a tool for predicting the binding affinity of reversible covalent inhibitors of hCatL. Our strategy relies on the relative free energy calculated for both covalent and noncovalent complexes and the free energy changes have been compared with experimental data for eight nitrile-based inhibitors, including three new inhibitors of hCatL. Our results demonstrate that the covalent complex can be employed to properly rank the inhibitors. Nevertheless, a comparison of the free energy changes in both noncovalent and covalent states is valuable to interpret the effect triggered by the formation of the covalent bond on the interactions played by functional groups distant from the warhead. Overall, FEP can be employed as a powerful predictor tool in developing and understanding the activity of reversible covalent inhibitors.

New functionalized 6-thienylpyrimidine-5-carbonitriles as antiproliferative agents against human breast cancer cells

6-Thienylpyrimidine-5-carbonitrile derivatives were synthesized and screened for their in vitro antiproliferative activities against two human breast cancer cell lines in comparison to 5-fluorouracil as a reference. Compounds 2, 3a-c, and 6b evolved as the most active congeners against both cell lines, while others showed selectivity for only one cell line. Compound 2 exerted its effect through inhibition of the epidermal growth factor receptor (EGFR), while 6b showed less aromatase inhibitory activity than letrozole. The rest of the tested compounds did not show significant inhibition, and it can be assumed that they exert their antiproliferative activity through different target mechanisms. In addition, caspase-9 protein activation assays, cell cycle analysis using flow cytometry, and annexin V-fluorescein isothiocyanate-propidium iodide (FITC/PI) dual staining assays were performed for the most active compounds. All the tested compounds were found to be potent pyrimidine derivatives able to initiate apoptosis in MCF-7 and MDA-MB-231 cells.

Catalytic Enantioselective Strecker Reaction of Isatin-Derived N-Unsubstituted Ketimines

A catalytic enantioselective Strecker reaction of isatin-derived N-unsubstituted ketimines directly afforded the N-unprotected α-aminonitriles with a tetrasubstituted carbon stereocenter in up to 99% ee without requiring protection/deprotection steps. One-pot Strecker reactions from the parent carbonyl compounds were also realized with comparable yields and enantioselectivities. Direct transformations of the N-unprotected α-aminonitrile products streamlined the synthesis of unnatural amino acid derivatives and achieved the shortest one-pot stereoselective routes to a biologically active compound reported to date.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

QSAR Modeling for Multi-Target Drug Discovery: Designing Simultaneous Inhibitors of Proteins in Diverse Pathogenic Parasites

Parasitic diseases remain as unresolved health issues worldwide. While for some parasites the treatments involve drug combinations with serious side effects, for others, chemical therapies are inefficient due to the emergence of drug resistance. This urges the search for novel antiparasitic agents able to act through multiple mechanisms of action. Here, we report the first multi-target model based on quantitative structure-activity relationships and a multilayer perceptron neural network (mt-QSAR-MLP) to virtually design and predict versatile inhibitors of proteins involved in the survival and/or infectivity of different pathogenic parasites. The mt-QSAR-MLP model exhibited high accuracy (>80%) in both training and test sets for the classification/prediction of protein inhibitors. Several fragments were directly extracted from the physicochemical and structural interpretations of the molecular descriptors in the mt-QSAR-MLP model. Such interpretations enabled the generation of four molecules that were predicted as multi-target inhibitors against at least three of the five parasitic proteins reported here with two of the molecules being predicted to inhibit all the proteins. Docking calculations converged with the mt-QSAR-MLP model regarding the multi-target profile of the designed molecules. The designed molecules exhibited drug-like properties, complying with Lipinski’s rule of five, as well as Ghose’s filter and Veber’s guidelines.

Identification of tectoridin as the inhibitor of FTO by isothermal titration calorimetric and spectroscopic methods

The characteristics of binding between tectoridin and the fat mass and obesity-associated protein were investigated.

Design, synthesis and stepwise optimization of nitrile-based inhibitors of cathepsins B and L

Human cathepsin B (CatB) is an important biological target in cancer therapy. In this work, we performed a knowledge-based design approach and the synthesis of a new set of 19 peptide-like nitrile-based cathepsin inhibitors. Reported compounds were assayed against a panel of human cysteine proteases: CatB, CatL, CatK, and CatS. Three compounds (7h, 7i, and 7j) displayed nanomolar inhibition of CatB and selectivity over CatK and CatL. The selectivity was achieved by using the combination of a para biphenyl ring at P3, halogenated phenylalanine in P2 and Thr-O-Bz group at P1. Likewise, compounds 7i and 7j showed selective CatB inhibition among the panel of enzymes studied. We have also described a successful example of bioisosteric replacement of the amide bond for a sulfonamide one [7e → 6b], where we observed an increase in affinity and selectivity for CatB while lowering the compound lipophilicity (ilogP). Our knowledge-based design approach and the respective structure-activity relationships provide insights into the specific ligand-target interactions for therapeutically relevant cathepsins.

Update on relevant trypanosome peptidases: Validated targets and future challenges

Trypanosoma cruzi, the agent of the American Trypanosomiasis, Chagas disease, and Trypanosoma brucei gambiense and Trypanosoma brucei rhodesiense, the agents of Sleeping sickness (Human African Trypanosomiasis, HAT), as well as Trypanosoma brucei brucei, the agent of the cattle disease nagana, contain cysteine, serine, threonine, aspartyl and metallo peptidases. The most abundant among these enzymes are the cysteine proteases from the Clan CA, the Cathepsin L-like cruzipain and rhodesain, and the Cathepsin B-like enzymes, which have essential roles in the parasites and thus are potential targets for chemotherapy. In addition, several other proteases, present in one or both parasites, have been characterized, and some of them are also promising candidates for the developing of new drugs. Recently, new inhibitors, with good selectivity for the parasite proteasomes, have been described and are very promising as lead compounds for the development of new therapies for these neglected diseases. This article is part of a Special Issue entitled: "Play and interplay of proteases in health and disease".

Crystal structure of Leishmania mexicana cysteine protease B in complex with a high-affinity azadipeptide nitrile inhibitor

Leishmania mexicana is an obligate intracellular protozoan parasite that causes the cutaneous form of leishmaniasis affecting South America and Mexico. The cysteine protease LmCPB is essential for the virulence of the parasite and therefore, it is an appealing target for antiparasitic therapy. A library of nitrile-based cysteine protease inhibitors was screened against LmCPB to develop a treatment of cutaneous leishmaniasis. Several compounds are sufficiently high-affinity LmCPB inhibitors to serve both as starting points for drug discovery projects and as probes for target validation. A 1.4 Å X ray crystal structure, the first to be reported for LmCPB, was determined for the complex of this enzyme covalently bound to an azadipeptide nitrile ligand. Mapping the structure-activity relationships for LmCPB inhibition revealed superadditive effects for two pairs of structural transformations. Therefore, this work advances our understanding of azadipeptidyl and dipeptidyl nitrile structure-activity relationships for LmCPB structure-based inhibitor design. We also tested the same series of inhibitors on related cysteine proteases cathepsin L and Trypanosoma cruzi cruzain. The modulation of these mammalian and protozoan proteases represents a new framework for targeting papain-like cysteine proteases.

Discovery of new pyrimidine-5-carbonitrile derivatives as anticancer agents targeting EGFRWT and EGFRT790M

A new series of pyrimidine-5-carbonitrile derivatives has been designed as ATP mimicking tyrosine kinase inhibitors of the epidermal growth factor receptor (EGFR). These compounds were synthesized and evaluated for their in vitro cytotoxic activities against a panel of four human tumor cell lines, namely colorectal carcinoma (HCT-116), hepatocellular carcinoma (HepG-2), breast cancer (MCF-7), and non-small cell lung cancer cells (A549). Five of the synthesized compounds, 11a, 11b, 12b, 15b and 16a, were found to exhibit moderate antiproliferative activity against the tested cell lines and were more active than the EGFR inhibitor erlotinib. In particular, compound 11b showed 4.5- to 8.4-fold erlotinib activity against HCT-116, HepG-2, MCF-7, and A549 cells with IC50 values of 3.37, 3.04, 4.14, and 2.4 μM respectively. Moreover, the most cytotoxic compounds that showed promising IC50 values against the four cancer cell lines were subjected to further investigation for their kinase inhibitory activities against EGFRWT and EGFRT790M using homogeneous time resolved fluorescence (HTRF) assay. Compound 11b was also found to be the most active compound against both EGFRWT and mutant EGFRT790M, exhibiting IC50 values of 0.09 and 4.03 μM, respectively. The cell cycle and apoptosis analyses revealed that compound 11b can arrest the cell cycle at the G2/M phase and induce significant apoptotic effects in HCT-116, HepG-2, and MCF-7 cells. Additionally, compound 11b upregulated the level of caspase-3 by 6.5 fold in HepG-2 when compared with the control. Finally, molecular docking studies were carried out to examine the binding mode of the synthesized compounds against the proposed targets; EGFRWT and EGFRT790M. Additional in silico ADMET studies were performed to explore drug-likeness properties.

Synthesis and matched molecular pair analysis of covalent reversible inhibitors of the cysteine protease CPB

Cysteine protease B (CPB) can be targeted by reversible covalent inhibitors that could serve as antileishmanial compounds. Here, sixteen dipeptidyl nitrile derivatives were synthesized, tested against CPB, and analyzed using matched molecular pairs to determine the effects of stereochemistry and p-phenyl substitution on enzyme inhibition. The compound (S)-2-(((S)-1-(4-bromophenyl)-2,2,2-trifluoroethyl)amino)-N-(1-cyanocyclopropyl)-3-phenylpropanamide (5) was the most potent CPB inhibitor (pKi = 6.82), which was also selective for human cathepsin B (pKi < 5). The inversion of the stereochemistry from S to R was more detrimental to potency when placed at the P2 position than at P3. The p-Br derivatives were more potent than the p-CH₃ and p-OCH₃ derivatives, probably due to intermolecular interactions with the S3 subsite.

On the intrinsic reactivity of highly potent trypanocidal cruzain inhibitors

The cysteine protease cruzipain is considered to be a validated target for therapeutic intervention in the treatment of Chagas disease. Hence, peptidomimetic cruzipain inhibitors having a reactive group (known as warhead) are subject to continuous studies to discover novel antichagasic compounds. Here, we evaluated how different warheads for a set of structurally similar related compounds could inhibit the activity of cruzipain and, ultimately, their trypanocidal effect. We first investigated in silico the intrinsic reactivity of these compounds by applying the Fukui index to correlate it with the enzymatic affinity. Then, we evaluated their potency against T. cruzi (Y and Tulahuen strains), which revealed the reversible cruzain inhibitor Neq0656 as a better trypanocidal agent (EC^Y.strain ₅₀ = 0.1 μM; SI = 58.4) than the current drug benznidazole (EC^Y.strain ₅₀ = 5.1 μM; SI > 19.6). We also measured the half-life time by HPLC analysis of three lead compounds in the presence of glutathione and cysteine to experimentally assess their intrinsic reactivity. Results clearly illustrated the reactivity trend for the warheads (azanitrile > aldehyde > nitrile), where the aldehyde displayed an intermediate intrinsic reactivity. Therefore, the aldehyde bearing peptidomimetic compounds should be subject for in-depth evaluation in the drug discovery process.

Optimization strategy of single-digit nanomolar cross-class inhibitors of mammalian and protozoa cysteine proteases

Cysteine proteases (CPs) are involved in a myriad of actions that include not only protein degradation, but also play an essential biological role in infectious and systemic diseases such as cancer. CPs also act as biomarkers and can be reached by active-based probes for diagnostic and mechanistic purposes that are critical in health and disease. In this paper, we present the modulation of a CP panel of parasites and mammals (Trypanosoma cruzi cruzain, LmCPB, CatK, CatL and CatS), whose inhibition by nitrile peptidomimetics allowed the identification of specificity and selectivity for a given CP. The activity cliffs identified at the CP inhibition level are useful for retrieving trends through multiple structure-activity relationships. For two of the cruzain inhibitors (10g and 4e), both enthalpy and entropy are favourable to Gibbs binding energy, thus overcoming enthalpy-entropy compensation (EEC). Group contribution of individual molecular modification through changes in enthalpy and entropy results in a separate partition on the relative differences of Gibbs binding energy (ΔΔG). Overall, this study highlights the role of CPs in polypharmacology and multi-target screening, which represents an imperative trend in the actual drug discovery effort.

Assessment of the Cruzain Cysteine Protease Reversible and Irreversible Covalent Inhibition Mechanism

Reversible and irreversible covalent ligands are advanced cysteine protease inhibitors in the drug development pipeline. K777 is an irreversible inhibitor of cruzain, a necessary enzyme for the survival of the Trypanosoma cruzi (T. cruzi) parasite, the causative agent of Chagas disease. Despite their importance, irreversible covalent inhibitors are still often avoided due to the risk of adverse effects. Herein, we replaced the K777 vinyl sulfone group with a nitrile moiety to obtain a reversible covalent inhibitor (Neq0682) of cysteine protease. Then, we used advanced experimental and computational techniques to explore details of the inhibition mechanism of cruzain by reversible and irreversible inhibitors. The isothermal titration calorimetry (ITC) analysis shows that inhibition of cruzain by an irreversible inhibitor is thermodynamically more favorable than by a reversible one. The hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) simulations were used to explore the mechanism of the reaction inhibition of cruzain by K777 and Neq0682. The calculated free energy profiles show that the Cys25 nucleophilic attack and His162 proton transfer occur in a single step for a reversible inhibitor and two steps for an irreversible covalent inhibitor. The hybrid QM/MM calculated free energies for the inhibition reaction correspond to -26.7 and -5.9 kcal mol^-1 for K777 and Neq0682 at the MP2/MM level, respectively. These results indicate that the ΔG of the reaction is very negative for the process involving K777, consequently, the covalent adduct cannot revert to a noncovalent protein-ligand complex, and its binding tends to be irreversible. Overall, the present study provides insights into a covalent inhibition mechanism of cysteine proteases.

Evaluating QM/MM Free Energy Surfaces for Ranking Cysteine Protease Covalent Inhibitors

One tactic for cysteine protease inhibition is to form a covalent bond between an electrophilic atom of the inhibitor and the thiol of the catalytic cysteine. In this study, we evaluate the reaction free energy obtained from a hybrid quantum mechanical/molecular mechanical (QM/MM) free energy profile as a predictor of affinity for reversible, covalent inhibitors of rhodesain. We demonstrate that the reaction free energy calculated with the PM6/MM potential is in agreement with the experimental data and suggest that the free energy profile for covalent bond formation in a protein environment may be a useful tool for the inhibitor design.

Quantum Mechanics/Molecular Mechanics Studies of the Mechanism of Cysteine Proteases Inhibition by Dipeptidyl Nitroalkenes

In this work a computational study of the mechanism of inhibition of cruzain, rhodesain, and cathepsin L cysteine proteases by the dipeptidyl nitroalkene Cbz-Phe-Ala-CH=CH-NO₂ has been carried out by means of molecular dynamics simulations with hybrid QM/MM potentials. The free-energy surfaces confirmed that the inhibition takes place by the formation of a covalent bond between the protein and the β-carbon atom of the inhibitor. According to the results, the tested inhibitor should be a much more efficient inhibitor of cruzain than of rhodesain, and little activity would be expected against cathepsin L, in total correspondence with the available experimental data. The origin of these differences may lie in the different stabilizing electrostatic interactions established between the inhibitor and the residues of the active site and S2 pocket of these enzymes. These results may be useful for the rational design of new dipeptidyl nitroalkenes with higher and more selective inhibitory activity against cysteine proteases.

Predicting the affinity of halogenated reversible covalent inhibitors through relative binding free energy

Nitrile reversible covalent inhibitors of human cathepsin L (hCatL) bind covalently to the side chain of the catalytic Cys25 residue in the S1 pocket to form thioimidates. Predicting the binding of reversible covalent inhibitors is essential for their practical application in drug design. In this report, five nitrile-based inhibitors coded Neq0570, Neq0710, Neq0802, Neq0803 and Neq0804 had their hCatL inhibition constants, K_i, determined. These analogs of the prototypical Neq0570 are halogenated reversible covalent inhibitors of hCatL, which bear a halogen atom in the meta position of the P3 benzyl ring that can form a halogen bond with the Gly61 of the hCatL. To describe halogen bonding interaction in an inhibitor-hCatL complex, we applied an extra point (EP) of charge to represent the anisotropic distribution of charge on the iodine, bromine and chlorine atoms. Besides, we have used alchemical free energy calculations for evaluating the overall relative binding free energies of these inhibitors using a two-state binding model: noncovalent and covalent bond states. Our results show that free energy perturbation (FEP) can predict the hCatL binding affinities of halogenated reversible covalent inhibitors in close agreement with experiments.

Ligand-induced conformational selection predicts the selectivity of cysteine protease inhibitors

Cruzain, a cysteine protease of Trypanosoma cruzi, is a validated target for the treatment of Chagas disease. Due to its high similarity in three-dimensional structure with human cathepsins and their sequence identity above 70% in the active site regions, identifying potent but selective cruzain inhibitors with low side effects on the host organism represents a significant challenge. Here a panel of nitrile ligands with varying potencies against cathepsin K, cathepsin L and cruzain, are studied by molecular dynamics simulations as both non-covalent and covalent complexes. Principal component analysis (PCA), identifies and quantifies patterns of ligand-induced conformational selection that enable the construction of a decision tree which can predict with high confidence a low-nanomolar inhibitor of each of three proteins, and determine the selectivity for one against others.