Cruzain Inhibitors Research Articles

Combining Charge Density Analysis with Machine Learning Tools To Investigate the Cruzain Inhibition Mechanism.

Trypanosoma cruzi, a flagellate protozoan parasite, is responsible for Chagas disease. The parasite major cysteine protease, cruzain (Cz), plays a vital role at every stage of its life cycle and the active-site region of the enzyme, similar to those of other members of the papain superfamily, is well characterized. Taking advantage of structural information available in public databases about Cz bound to known covalent inhibitors, along with their corresponding activity annotations, in this work, we performed a deep analysis of the molecular interactions at the Cz binding cleft, in order to investigate the enzyme inhibition mechanism. Our toolbox for performing this study consisted of the charge density topological analysis of the complexes to extract the molecular interactions and machine learning classification models to relate the interactions with biological activity. More precisely, such a combination was useful for the classification of molecular interactions as “active-like” or “inactive-like” according to whether they are prevalent in the most active or less active complexes, respectively. Further analysis of interactions with the help of unsupervised learning tools also allowed the understanding of how these interactions come into play together to trigger the enzyme into a particular conformational state. Most active inhibitors induce some conformational changes within the enzyme that lead to an overall better fit of the inhibitor into the binding cleft. Curiously, some of these conformational changes can be considered as a hallmark of the substrate recognition event, which means that most active inhibitors are likely recognized by the enzyme as if they were its own substrate so that the catalytic machinery is arranged as if it is about to break the substrate scissile bond. Overall, these results contribute to a better understanding of the enzyme inhibition mechanism. Moreover, the information about main interactions extracted through this work is already being used in our lab to guide docking solutions in ongoing prospective virtual screening campaigns to search for novel noncovalent cruzain inhibitors.

Synthesis and structure-activity relationship of nitrile-based cruzain inhibitors incorporating a trifluoroethylamine-based P2 amide replacement

The structure-activity relationship for nitrile-based cruzain inhibitors incorporating a P2 amide replacement based on trifluoroethylamine was explored by deconstruction of a published series of inhibitors. It was demonstrated that the P3 biphenyl substituent present in the published inhibitor structures could be truncated to phenyl with only a small loss of affinity. The effects of inverting the configuration of the P2 amide replacement and linking a benzyl substituent at P1 were observed to be strongly nonadditive. We show that plotting affinity against molecular size provides a means to visualize both the molecular size efficiency of structural transformations and the nonadditivity in the structure-activity relationship. We also show how the relationship between affinity and lipophilicity, measured by high-performance liquid chromatography with an immobilized artificial membrane stationary phase, may be used to normalize affinity with respect to lipophilicity.

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 Ki=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.

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.

Cruzain structures: apocruzain and cruzain bound to S-methyl thiomethanesulfonate and implications for drug design.

Chagas disease, which is caused by Trypanosoma cruzi, affects more than six million people worldwide. Cruzain is the major cysteine protease involved in the survival of this parasite. Here, the expression, purification and crystallization of this enzyme are reported. The cruzain crystals diffracted to 1.2 Å resolution, yielding two novel cruzain structures: apocruzain and cruzain bound to the reversible covalent inhibitor S-methyl thiomethanesulfonate. Mass-spectrometric experiments confirmed the presence of a methylthiol group attached to the catalytic cysteine. Comparison of these structures with previously published structures indicates the rigidity of the cruzain structure. These results provide further structural information about the enzyme and may help in new in silico studies to identify or optimize novel prototypes of cruzain inhibitors.

Cruzain Inhibitors as Prominent Molecules with The Potential to become Drug Candidates against Chagas Disease

Cruzain Inhibitors as Prominent Molecules with The Potential to become Drug Candidates against Chagas Disease

Structure-Based Virtual Screening and In Vitro Evaluation of New Trypanosoma cruzi Cruzain Inhibitors.

Chagas disease (CD), or American trypanosomiasis, causes more than 10,000 deaths per year in the Americas. Current medical therapy for CD has low efficacy in the chronic phase of the disease and serious adverse effects; therefore, it is necessary to search for new pharmacological treatments. In this work, the ZINC15 database was filtered using the N-acylhydrazone moiety and a subsequent structure-based virtual screening was performed using the cruzain enzyme of Trypanosoma cruzi to predict new potential cruzain inhibitors. After a rational selection process, four compounds, Z2 (ZINC9873043), Z3 (ZINC9870651), Z5 (ZINC9715287), and Z6 (ZINC9861447), were chosen to evaluate their in vitro trypanocidal activity and enzyme inhibition. Compound Z5 showed the best trypanocidal activity against epimatigote (IC50 = 36.26 ± 9.9 μM) and trypomastigote (IC50 = 166.21 ± 14.5 μM and 185.1 ± 8.5 μM on NINOA and INC-5 strains, respectively) forms of Trypanosoma cruzi. In addition, Z5 showed a better inhibitory effect on Trypanosoma cruzi proteases than S1 (STK552090, 8-chloro-N-(3-morpholinopropyl)-5H-pyrimido[5,4-b]-indol-4-amine), a known cruzain inhibitor. This study encourages the use of computational tools for the rational search for trypanocidal drugs.

Carbamoyl Imidazoles As Potent, Reversible and Competitive Cruzain Inhibitors with <i>in vitro</i> and <i>in vivo</i> Trypanocidal Activity: A Structure-Based Drug Design Approach

A virtual screening conducted with nearly 4,000,000 compounds from lead-like and fragment-like subsets enabled the identification of a small-molecule inhibitor (1) of the Trypanosoma cruzi cruzain enzyme, a validated drug target for Chagas disease. Subsequent comprehensive structure-based drug design and structure-activity relationship studies led to the discovery of carbamoyl imidazoles as potent, reversible and competitive cruzain inhibitors. The most potent carbamoyl imidazole inhibitor (45) exhibited high affinity with a Ki value of 20 nM, presenting both in vitro and in vivo activity against T. cruzi. Furthermore, the most promising compounds reduced parasite burden in vivo and showed no toxicity at a dose of 100 mg/kg. These carbamoyl imidazoles are structurally attractive, nonpeptidic, and easy to prepare and synthetically modify. Finally, these results further advance our understanding of the noncovalent mode of inhibition of this pharmaceutically relevant enzyme, building strong foundations for drug discovery efforts.

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

The protozoan cysteine proteases cruzain in Trypanosoma cruzi and rhodesain in Trypanosoma brucei are therapeutic targets for Chagas disease and Human African Trypanosomiasis (HAT), respectively. A benzimidazole series was previously characterized as potent noncovalent competitive cruzain and rhodesain inhibitors with activity against trypanosomes. Common structure-activity relationships (SAR) trends and structural modifications leading to selectivity against each enzyme were described. However, some of these trends could not be understood based on the reported binding mode of lead compound 1. Therefore, we employed microsecond molecular dynamics simulations and free energy calculations to understand qualitative SAR trends and to quantitatively recapitulate them. Simulations revealed the most stable protein-ligand interactions and provided insights concerning enzyme selectivity. Calculated relative binding free energies of compound 1 analogs exhibited deviations of 1.1 and 2.2 kcal/mol from the experimental values for cruzain and rhodesain, respectively. These data encourage prospective thermodynamic integration (TI) studies to optimize this series and facilitate the prioritization of compounds for synthesis.

Molecular description of α-keto-based inhibitors of cruzain with activity against Chagas disease combining 3D-QSAR studies and molecular dynamics.

In this work, a group of α-keto-based inhibitors of the cruzain enzyme with anti-chagas activity was selected for a three-dimensional quantitative structure-activity relationship study (3D-QSAR) combined with molecular dynamics (MD). Firstly, statistical models based on Partial Least Square (PLS) regression were developed employing comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) descriptors. Validation parameters (q2 and r2 )for the models were, respectively, 0.910 and 0.997 (CoMFA) and 0.913 and 0.992 (CoMSIA). In addition, external validation for the models using a test group revealed r2pred =0.728 (CoMFA) and 0.971 (CoMSIA). The most relevant aspect in this study was the generation of molecular fields in both favorable and unfavorable regions based on the models developed. These fields are important to interpret modifications necessary to enhance the biological activities of the inhibitors. This analysis was restricted considering the inhibitors in a fixed conformation, not interacting with their target, the cruzain enzyme. Then, MD was employed taking into account important variables such as time and temperature. MD helped describe the behavior of the inhibitors and their properties showed similar results as those generated by QSAR-3D study.

Leveraging the cruzain S3 subsite to increase affinity for reversible covalent inhibitors

Cruzain is the major cysteine protease of Trypanosoma cruzi, the etiological agent of Chagas disease. Reversible covalent cruzain inhibitors can block the steps of cell differentiation in the parasite and kill the organism. To this end, the description of how inhibitors modified at the P2/P3 positions lead to analogs with greater cruzain affinity to the S2/S3 subsites is of fundamental importance. Albeit many efforts are being employed in the characterization of the interaction processes with S2 subsite, little is known about the cruzain S3 subsite. In this work, we show a brief but consistent study to identify favorable substitutions in P3 of dipeptidyl nitriles that increase cruzain affinity. Using molecular dynamics simulations, we have identified some dipeptidyl nitrile analogs with modifications at P3 position that had higher cruzain inhibition than the original unsubstituted compound. A matched molecular pair analysis shows the importance of including a chlorine atom in the P3-meta position. The modifications implemented in P3 are confirmed when profiling the thermodynamic parameters via isothermal titration calorimetry. The classical enthalpy-entropy compensation phenomenon, in which enthalpy changes are counterbalanced by entropy results in a small modification of ΔG. The inclusion of the chlorine atom in the P3-meta position results in the highest reduction of the detrimental entropic contribution observed in P3.

Investigation of the binding mode of a novel cruzain inhibitor by docking, molecular dynamics, ab initio and MM/PBSA calculations.

Chagas disease remains a major health problem in South America, and throughout the world. The two drugs clinically available for its treatment have limited efficacy and cause serious adverse effects. Cruzain is an established therapeutic target of Trypanosoma cruzi, the protozoan that causes Chagas disease. Our group recently identified a competitive cruzain inhibitor (compound 1) with an IC50 = 15 µM that is also more synthetically accessible than the previously reported lead, compound 2. Prior studies, however, did not propose a binding mode for compound 1, hindering understanding of the structure-activity relationship and optimization. Here, the cruzain binding mode of compound 1 was investigated using docking, molecular dynamics (MD) simulations with ab initio derived parameters, ab initio calculations, and MM/PBSA. Two ligand protonation states and four binding poses were evaluated. A careful ligand parameterization method was employed to derive more physically meaningful parameters than those obtained by automated tools. The poses of unprotonated 1 were unstable in MD, showing large conformational changes and diffusing away from the binding site, whereas the protonated form showed higher stability and interaction with negatively charged residues Asp161 and Cys25. MM/PBSA also suggested that these two residues contribute favorably to binding of compound 1. By combining results from MD, ab initio calculations, and MM/PBSA, a binding mode of 1 is proposed. The results also provide insights for further optimization of 1, an interesting lead compound for the development of new cruzain inhibitors.

Molecular Docking Studies Applied to a Dataset of Cruzain Inhibitors

Chagas' disease is one of the main causes of heart failure in developing countries. The disadvantages of current therapy include the undesirable side-effects, resistance, and therapeutic adhesion. The development of new efficient and safe drugs is, therefore, an issue of extreme importance. In order to gain a better understanding of how the compounds interact with the target, computational methods are essential. In this theoretical study, we report a docking protocol applied to a dataset of 173 cruzain inhibitors with IC50 values of less than 10 µM, belonging 16 different chemical classes. A preliminary analysis was performed, where the best protein structure for the study was identified. The enzyme was validated by redocking and a fingerprint graph for the ligand-enzyme interactions was generated, allowing the identification of the main amino acid residues related to the activity. Additionally, a larger cluster was generated, allowing the visualization of the orientation of the compounds and providing binding information for the different classes of compounds as well as their interaction in the cruzain active site. Amino acid residues other than those known as the catalytic triad (Gly23, Cys25, and Gly65) were identified, for example, Gln19 and Asp158. This provides a better insight into the mode of interaction of various cruzain inhibitors, which show IC50 values in the nanomolar range but which do not interact with the triad. These findings can help researchers to find new cruzain inhibitors for use in the fight against the Chagas disease.

A comparative study of warheads for design of cysteine protease inhibitors

The effects on potency of cruzain inhibition of replacing a nitrile group with alternative warheads were explored. The oxime was almost an order of magnitude more potent than the corresponding nitrile and has the potential to provide access to the prime side of the catalytic site. Dipeptide aldehydes and azadipeptide nitriles were found to be two orders of magnitude more potent cruzain inhibitors than the corresponding dipeptide nitriles although potency differences were modulated by substitution at P1 and P3. Replacement of the α methylene of a dipeptide aldehyde with cyclopropane led to a loss of potency of almost three orders of magnitude. The vinyl esters and amides that were characterized as reversible inhibitors were less potent than the corresponding nitrile by between one and two orders of magnitude.

Development of a pharmacophore for cruzain using oxadiazoles as virtual molecular probes: quantitative structure-activity relationship studies.

Chagas's is a neglected tropical disease caused by the protozoan parasite Trypanosoma cruzi. According to the World Health Organization, 7million people are infected worldwide leading to 7000 deaths per year. Drugs available, nifurtimox and benzimidazole, are limited due to low efficacy and high toxicity. As a validated target, cruzain represents a major front in drug discovery attempts for Chagas disease. Herein, we describe the development of 2D QSAR ([Formula: see text] = 0.81) and a 3D-QSAR-based pharmacophore ([Formula: see text] = 0.82) from a series of non-covalent cruzain inhibitors represented mostly by oxadiazoles (lead compound, IC50 = 200nM). Both models allowed us to map key intermolecular interactions in S1', S2 and S3 cruzain sub-sites (including halogen bond and C‒H/π). To probe the predictive capacity of obtained models, inhibitors available in the literature from different classes displaying a range of scaffolds were evaluate achieving mean absolute deviation of 0.33 and 0.51 for 2D and 3D models, respectively. CoMFA revealed an unexplored region where addition of bulky substituents to produce new compounds in the series could be beneficial to improve biological activity.

Identification of Structure-Stabilizing Interactions in Enzymes: A Novel Mechanism to Impact Enzyme Activity.

Cruzain, a cysteine protease in the cathepsin family, is pivotal to the life-cycle of Trypanosoma cruzi, the etiological agent in Chagas disease. Current inhibitors of cruzain suffer from drawbacks involving gastrointestinal and neurological side effects and as a result have spurred the search for alternative anti-trypanocidals. Through sequence alignment studies and intra-residue interaction analysis of the pro-protein of cruzain (pro-cruzain), we have identified a host of non-active site residues that are conserved among the cathepsins. We hypothesize that these conserved amino acids play a critical role in structure-stabilizing interactions among the cathepsins and are therefore crucial for eventually gaining protease activity. As predicted, mutation of selected conserved non-active site amino-acid candidates in cruzain resulted in a compromised structural stability and a corresponding loss in enzymatic activity relative to wild-type enzyme. By advancing the discovery of novel, non-active-site-based targets to arrest enzymatic activity our results potentially open the field of alternative inhibitor design. The advantages of defining such a non-active-site inhibitor design space is discussed.

Repositioning FDA Drugs as Potential Cruzain Inhibitors from Trypanosoma cruzi: Virtual Screening, In Vitro and In Vivo Studies.

Chagas disease (CD) is a neglected disease caused by the parasite Trypanosoma cruzi, which affects underdeveloped countries. The current drugs of choice are nifurtimox and benznidazole, but both have severe adverse effects and less effectivity in chronic infections; therefore, the need to discover new drugs is essential. A computer-guided drug repositioning method was applied to identify potential FDA drugs (approved and withdrawn) as cruzain (Cz) inhibitors and trypanocidal effects were confirmed by in vitro and in vivo studies. 3180 FDA drugs were virtually screened using a structure-based approach. From a first molecular docking analysis, a set of 33 compounds with the best binding energies were selected. Subsequent consensus affinity binding, ligand amino acid contact clustering analysis, and ranked position were used to choose four known pharmacological compounds to be tested in vitro. Mouse blood samples infected with trypomastigotes from INC-5 and NINOA strains were used to test the trypanocidal effect of four selected compounds. Among these drugs, one fibrate antilipemic (etofyllin clofibrate) and three β-lactam antibiotics (piperacillin, cefoperazone, and flucloxacillin) showed better trypanocidal effects (LC50 range 15.8–26.1 μg/mL) in comparison with benznidazole and nifurtimox (LC50 range 33.1–46.7 μg/mL). A short-term in vivo evaluation of these compounds showed a reduction of parasitemia in infected mice (range 90–60%) at 6 h, but this was low compared to benznidazole (50%). This work suggests that four known FDA drugs could be used to design and obtain new trypanocidal agents.

Molecular modeling and structure-activity relationships for a series of benzimidazole derivatives as cruzain inhibitors.

Chagas disease is endemic in Latin America and no effective treatment is available. Efforts in drug research have focused on several enzymes from Trypanosoma cruzi, among which cruzain is a validated pharmacological target. Chemometric analyses were performed on the data set using the hologram quantitative structure-activity relationship, comparative molecular field analysis and comparative molecular similarity index analysis methods. Docking simulations were executed using the crystallographic structure of cruzain in complex with a benzimidazole inhibitor. The top-scoring enzyme-inhibitor complexes were selected for the development of the 3D quantitative structure-activity relationship (QSAR) models and to assess the inhibitor binding modes and intermolecular interactions. Benzimidazole derivatives as cruzain inhibitors were used in molecular docking and QSAR studies. Significant statistical indicators were obtained, and the best models demonstrated high predictive ability for an external test set (r 2pred = 0.65, 0.94 and 0.82 for hologram QSAR, comparative molecular field analysis and comparative molecular similarity index analysis, respectively). Additionally, the graphical information of the chemometric analyses demonstrated substantial complementarity with the enzyme-binding site. These results demonstrate the relevance of the QSAR models to guide the design of structurally related benzimidazole derivatives with improved potency.

Synthesis and biological evaluation of potential inhibitors of the cysteine proteases cruzain and rhodesain designed by molecular simplification

Analogues of 8-chloro-N-(3-morpholinopropyl)-5H-pyrimido[5,4-b]indol-4-amine 1, a known cruzain inhibitor, were synthesized using a molecular simplification strategy. Five series of analogues were obtained: indole, pyrimidine, quinoline, aniline and pyrrole derivatives. The activity of the compounds was evaluated against the enzymes cruzain and rhodesain as well as against Trypanosoma cruzi amastigote and trypomastigote forms. The 4-aminoquinoline derivatives showed promising activity against both enzymes, with IC50 values ranging from 15 to 125µM. These derivatives were selective inhibitors for the parasitic proteases, being unable to inhibit mammalian cathepsins B and S. The most active compound against cruzain (compound 5a; IC50=15µM) is considerably more synthetically accessible than 1, while retaining its ligand efficiency. As observed for the original lead, compound 5a was shown to be a competitive enzyme inhibitor. In addition, it was also active against T. cruzi (IC50=67.7µM). Interestingly, the pyrimidine derivative 4b, although inactive in enzymatic assays, was highly active against T. cruzi (IC50=3.1µM) with remarkable selectivity index (SI=128) compared to uninfected fibroblasts. Both 5a and 4b exhibit drug-like physicochemical properties and are predicted to have a favorable ADME profile, therefore having great potential as candidates for lead optimization in the search for new drugs to treat Chagas disease.

Synthesis, Biological Evaluation and Molecular Docking of New Benzenesulfonylhydrazone as Potential anti-Trypanosoma cruzi Agents

Chagas disease is a public health problem caused by Trypanosoma cruzi. Cruzain is a pharmacological target for designing a new drug against this parasite. Hydrazone and Nacylhydrazone derivatives have been traditionally associated as potential Cruzain inhibitors. Additionally, benzenesulfonyl derivatives show trypanocidal activity. Therefore, in this study, the combination of both structures has been taken into account for drug design. Seven benzenesulfonylhydrazone (BS-H) and seven N-propionyl benzenesulfonylhydrazone (BS-NAH) derivatives were synthetized and elucidated by infrared spectroscopy, nuclear magnetic resonance, and elemental analysis. All compounds were evaluated biologically in vitro against two strains of Trypanosoma cruzi (NINOA and INC-5), which are endemic in Mexico, and compared with the reference drugs nifurtimox and benznidazole. In order to gain insight into the putative molecular origin of the trypanocidal properties of these derivatives, docking studies were carried out with Cruzain. Compounds 4 and 6 (BS-H) and 10, 12-14 (BS-NAH) showed the best biological activity against NINOA and INC-5 strains, respectively. Compound 13 was the most potent trypanocidal compound showing a LC50 of 0.06 µM against INC-5 strain. However, compound 4 showed the best activity against both strains (LC50 <30 µM). Theoretical binding modes obtained suggested covalent binding that could explain their biological activity. Benzenesulfonyl and N-propionyl benzenesulfonyl hydrazone derivatives are good options for developing new trypanocidal agents. Particularly, compound 4 could be considered a lead compound.