Understanding Structure–Activity Relationships for Trypanosomal Cysteine Protease Inhibitors by Simulations and Free Energy Calculations

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.

Molecular targets for Chagas disease: validation, challenges and lead compounds for widely exploited targets

INTRODUCTION

Chagas disease (CD) imposes social and economic burdens, yet the available treatments have limited efficacy in the disease's chronic phase and cause serious adverse effects. To address this challenge, target-based approaches are a possible strategy to develop new, safe, and active treatments for both phases of the disease.

AREAS COVERED

This review delves into target-based approaches applied to CD drug discovery, emphasizing the studies from the last five years. We highlight the proteins cruzain (CZ), trypanothione reductase (TR), sterol 14 α-demethylase (CPY51), iron superoxide dismutase (Fe-SOD), proteasome, cytochrome b (Cytb), and cleavage and polyadenylation specificity factor 3 (CPSF3), chosen based on their biological and chemical validation as drug targets. For each, we discuss its biological relevance and validation as a target, currently related challenges, and the status of the most promising inhibitors.

EXPERT OPINION

Target-based approaches toward developing potential CD therapeutics have yielded promising leads in recent years. We expect a significant advance in this field in the next decade, fueled by the new options for Trypanosoma cruzi genetic manipulation that arose in the past decade, combined with recent advances in computational chemistry and chemical biology.

Structure-based discovery of novel cruzain inhibitors with distinct trypanocidal activity profiles

Over 110 years after the first formal description of Chagas disease, the trypanocidal drugs thus far available have limited efficacy and several side effects. This encourages the search for novel treatments that inhibit T. cruzi targets. One of the most studied anti-T. cruzi targets is the cysteine protease cruzain; it is associated with metacyclogenesis, replication, and invasion of the host cells. We used computational techniques to identify novel molecular scaffolds that act as cruzain inhibitors. First, with a docking-based virtual screening, we identified compound 8, a competitive cruzain inhibitor with a K_i of 4.6 μM. Then, aided by molecular dynamics simulations, cheminformatics, and docking, we identified the analog compound 22 with a K_i of 27 μM. Surprisingly, despite sharing the same isoquinoline scaffold, compound 8 presented higher trypanocidal activity against the epimastigote forms, while compound 22, against the trypomastigotes and amastigotes. Taken together, compounds 8 and 22 represent a promising scaffold for further development of trypanocidal compounds as drug candidates for treating Chagas disease.

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.

Study of the dynamic behavior of the cruzain enzyme in free and complexed forms with competitive and noncovalent benzimidazole inhibitors

There are only two drugs for the treatment of Chagas disease, namely, nifurtimox and benznidazole, that can cause several adverse effects. Despite the effectiveness of these drugs in the disease's acute phase, they are not recognized as curative in the chronic phase, establishing the need for more effective treatment in all stages of the disease. Cruzain is an enzyme that plays a vital role in the life cycle of the etiologic agent, the protozoan Trypanosoma cruzi, being relevant as a therapeutic target in the planning of new drugs. Using molecular docking and dynamics simulations, we have investigated the structural and dynamic factors that can be involved in the enzyme inhibition process at the atomic-molecular level by benzimidazole compounds that are potent cruzain inhibitors with in vitro trypanocidal activity. The study suggests that these inhibitors bind cruzain through steric and hydrogen bonding interactions without altering its secondary structure content and protein compaction. Besides, we observed that these inhibitors decrease the correlation of movements between Cα-atoms of cruzain, increasing the number of atomic communities, mainly in the α-helix that presents the catalytic Cys25 residue. As expected, we also observed a correlation between the inhibitory activity of each inhibitor and their respective binding-free energies, reinforcing that the affinity of the complexes seems to be a relevant factor for enzymatic inhibition. Hence, the results presented in this work contribute to a better understanding of the cruzain enzyme inhibition mechanism through competitive and non-covalent inhibitors.Communicated by Ramaswamy H. Sarma.

Computational approaches towards the discovery and optimisation of cruzain inhibitors

The need to develop safer and more efficacious drugs to treat Chagas disease has motivated the search for cruzain inhibitors. Cruzain is the recombinant, truncated version of cruzipain, a cysteine protease from Trypanosoma cruzi with important roles during the parasite life cycle. Several computational techniques have been applied to discover and optimise cruzain inhibitors, providing a molecular basis to guide this process. Here, we review some of the most recent computational studies that provided important information for the design of cruzain inhibitors. Moreover, we highlight the diversity of applications of in silico techniques and their impact.

The gene repertoire of the main cysteine protease of Trypanosoma cruzi, cruzipain, reveals four sub-types with distinct active sites

Cruzipains are the main papain-like cysteine proteases of Trypanosoma cruzi, the protozoan parasite that causes Chagas disease. Encoded by a multigenic family, previous studies have estimated the presence of dozens of copies spread over multiple chromosomes in different parasite strains. Here, we describe the complete gene repertoire of cruzipain in three parasite strains, their genomic organization, and expression pattern throughout the parasite life cycle. Furthermore, we have analyzed primary sequence variations among distinct family members as well as structural differences between the main groups of cruzipains. Based on phylogenetic inferences and residue positions crucial for enzyme function and specificity, we propose the classification of cruzipains into two families (I and II), whose genes are distributed in two or three separate clusters in the parasite genome, according with the strain. Family I comprises nearly identical copies to the previously characterized cruzipain 1/cruzain, whereas Family II encompasses three structurally distinct sub-types, named cruzipain 2, cruzipain 3, and cruzipain 4. RNA-seq data derived from the CL Brener strain indicates that Family I genes are mainly expressed by epimastigotes, whereas trypomastigotes mainly express Family II genes. Significant differences in the active sites among the enzyme sub-types were also identified, which may play a role in their substrate selectivity and impact their inhibition by small molecules.

Elucidating the Dual Mode of Action of Dipeptidyl Enoates in the Inhibition of Rhodesain Cysteine Proteases

A computational study of the two possible inhibition mechanisms of rhodesain cysteine protease by the dipeptidyl enoate Cbz-Phe-Leu-CH=CH-CO₂ C₂ H₅ has been carried out by means of molecular dynamics simulations with hybrid QM/MM potentials. The low free energy barriers confirm that the Cys25 residue can attack both C_β and C1 atoms of the inhibitor, confirming a dual mode of action in the inhibition of the rhodesain by enoates. According to the results, the inhibition process through the Cys25 attack on the C_β atom of the inhibitor is an exergonic and irreversible process, while the inhibition process when Cys25 attacks on the C1 atom of the inhibitor is and exergonic but reversible process. The interactions between the inhibitor and rhodesain suggest that P2 is the most important fragment to consider in the design of new efficient inhibitors of rhodesain. These results may be useful for the design of new inhibitors of rhodesain and other related cysteine proteases based on dipeptidyl enoates scaffolds.

Shared Binding Mode of Perrottetinene and Tetrahydrocannabinol Diastereomers inside the CB1 Receptor May Incentivize Novel Medicinal Drug Design: Findings from an in Silico Assay

In recent years, therapeutic compounds derived from phytocannabinoids have brought renewed attention to the benefits they offer to ameliorate chronic disease symptoms. Among cannabinoids, tetrahydrocannabinol (THC) is a well-known component of the Cannabis plant, whose active principles have been studied through the years. Another psychoactive phytocannabinoid, derived from liverworts Radula, perrottetinene (PET), has created interest, especially as a pharmaceutical product and for its legal recreational use. Unfortunately, so far, the interaction mode of these compounds at the type 1 cannabinoid receptors (CB1R) binding site remains unknown, and no experimental three-dimensional structure in complex with THC or PET is available in the Protein Data Bank. Today, many computational methodologies can assist in this crusade and help unveil how these molecules bind, based on the already known pose of a structurally similar compound. In this work, we aim to elucidate the binding mode of THC and PET molecules in both cis and trans conformers, using a combination of several computational methodologies, including molecular docking, molecular dynamics, free energy calculations, and protein-energy network studies. We found that THC and PET interact similarly with the CB1R, in a different conformation depending on the considered diastereomer. We have observed that cis ligands adopted a half-chair conformation of the cycle ring containing the dimethyl group, assuming an axial or equatorial conformation producing a different induced fitting of the surrounding residues compared with trans ligands, with higher interaction energy than the trans conformer. For PET, we have seen that Trp-279 and Trp-356 have a marked influence on the binding. After binding, Trp-279 accommodates its side chain to better interact with the PET's terminal phenyl group, disturbing CB1R residues communication. The interaction with Trp-356 might impair the activation of CB1R and can influence the binding of PET as a partial agonist. Understanding the PET association with CB1R from a molecular perspective can offer a glimpse of preventing potential toxicological or recreational effects since it is an attractive lead for drug development with fewer side effects than trans-THC.

The role of imidazole and benzimidazole heterocycles in Chagas disease: A review

The haemoflagellate protozoan Trypanosoma cruzi (T. cruzi) is the causative agent of Chagas disease (CD), a potentially life-threatening disease. Little by little, remarkable progress has been achieved against CD, although it is still not enough. In the absence of effective chemotherapy, many research groups, organizations and pharmaceutical companies have focused their efforts on the search for compounds that could become viable drugs against CD. Within the wide variety of reported derivatives, this review summarizes and provides a global vision of the situation of those compounds that include broadly studied heterocycles in their structures due to their applications in medicinal chemistry: imidazole and benzimidazole rings. Therefore, the intention of this work is to present a compilation, as much as possible, of all the reported information, regarding these imidazole and benzimidazole derivatives against T. cruzi, as a starting point for future researchers in this field.

Profiling selectivity of chagasin mutants towards cysteine proteases cruzain or cathepsin L through molecular dynamics simulations

Chagasin, an endogenous cysteine protease inhibitor from Trypanosoma cruzi, can control the activity of the parasitic cruzain and its homologous human cathepsin L. While chagasin inhibits both enzymes with similar potency, mutations have different effects on binding to these enzymes. Mutants T31A and T31A/T32A bind well to cathepsin L, but their affinity for cruzain drops ∼40 to 140-fold. On the other hand, the mutant W93A binds well to cruzain, but it loses potency against cathepsin L. Here, we employed molecular dynamics simulations to understand the selectivity in inhibition of cruzain or cathepsin L by chagasin mutants W93A, T31A, and T31A/T32A. Our results allowed profiling the nonbonded interactions in the interfaces of each mutant with these cysteine proteases. Additionally, we observed differences in the binding conformation of the chagasin loops L2 and L6 of the W93A mutant, favoring interactions with cruzain and reducing interactions with cathepsin L. These differences are associated with a partial dissociation of the W93A-cathepsin L complex, providing a likely cause for the selectivity of the mutant W93A towards cruzain.Communicated by Ramaswamy H. Sarma.

Benzimidazole inhibitors of the major cysteine protease of Trypanosoma brucei

Aim: Limitations in available therapies for trypanosomiases indicate the need for improved medicines. Cysteine proteases cruzain and rhodesain are validated targets for treatment of Chagas disease and human African trypanosomiasis. Previous studies reported a benzimidazole series as potent cruzain inhibitors. Results & methodology: Considering the high similarity between these proteases, we evaluated 40 benzimidazoles against rhodesain. We describe their structure-activity relationships (SAR), revealing trends similar to those observed for cruzain and features that lead to enzyme selectivity. This series comprises noncovalent competitive inhibitors (best K_i = 0.21 μM against rhodesain) and micromolar activity against Trypanosoma brucei brucei. A cheminformatics analysis confirms scaffold novelty, and the inhibitors described have favorable predicted physicochemical properties. Conclusion: Our results support this series as a starting point for new human African trypanosomiasis medicines.

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.

Identification of potential cruzain inhibitors using de novo design, molecular docking and dynamics simulations studies

Communicated by Ramaswamy H. Sarma.

Discovery and characterization of trypanocidal cysteine protease inhibitors from the 'malaria box'

Chagas disease, Human African Trypanosomiasis, and schistosomiasis are neglected parasitic diseases for which new treatments are urgently needed. To identify new chemical leads, we screened the 400 compounds of the Open Access Malaria Box against the cysteine proteases, cruzain (Trypanosoma cruzi), rhodesain (Trypanosoma brucei) and SmCB1 (Schistosoma mansoni), which are therapeutic targets for these diseases. Whereas just three hits were observed for SmCB1, 70 compounds inhibited cruzain or rhodesain by at least 50% at 5 μM. Among those, 15 commercially available compounds were selected for confirmatory assays, given their potency, time-dependent inhibition profile and reported activity against parasites. Additional assays led to the confirmation of four novel classes of cruzain and rhodesain inhibitors, with potency in the low-to mid-micromolar range against enzymes and T. cruzi. Assays against mammalian cathepsins S and B revealed inhibitor selectivity for parasitic proteases. For the two competitive inhibitors identified (compounds 7 and 12), their binding mode was predicted by docking, providing a basis for structure-based optimization efforts. Compound 12 also acted directly against the trypomastigote and the intracellular amastigote forms of T. cruzi at 3 μM. Therefore, through a combination of experimental and computational approaches, we report promising hits for optimization in the development of new trypanocidal drugs.