Next Article in Journal
Study of the Chemical Recovery and Selectivity against U in the Radiochemical Separation of Th with Tri-n-butyl Phosphate by Varying the Proportion of Xylene and HCl Concentration
Previous Article in Journal
Trivalent Disulfide Unit-Masked System Efficiently Delivers Large Oligonucleotide
Previous Article in Special Issue
Participation of Oxidative Stress in the Activity of Compounds Isolated from Eleutherine plicata Herb
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Juan Pablo Cerutti
1,2,
Lucas Abreu Diniz
3,
Viviane Corrêa Santos
3,
Salomé Catalina Vilchez Larrea
4,
Guillermo Daniel Alonso
4,
Rafaela Salgado Ferreira
3,
Wim Dehaen
2 and
Mario Alfredo Quevedo
1,*
1
Unidad de Investigación y Desarrollo en Tecnología Farmacéutica (UNITEFA-CONICET), Facultad de Ciencias Químicas, Universidad Nacional de Córdoba (FCQ-UNC), Haya de la Torre y Medina Allende, Córdoba 5000, Argentina
2
Sustainable Chemistry for Metals and Molecules, Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
3
Departamento de Bioquímica e Imunologia, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Av. Antônio Carlos 6627, Belo Horizonte 31270-901, Brazil
4
Instituto de Investigaciones en Ingeniería Genética y Biología Molecular (INGEBI-CONICET), Vuelta de Obligado 2490, Ciudad de Buenos Aires 1428, Argentina
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(17), 4224; https://doi.org/10.3390/molecules29174224
Submission received: 3 July 2024 / Revised: 29 August 2024 / Accepted: 3 September 2024 / Published: 5 September 2024

Abstract

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

1. Introduction

Chagas disease (CD), also known as American trypanosomiasis, is defined by the WHO as a potentially life-threatening infectious systemic disease caused by the protozoan parasite Trypanosoma cruzi (T. cruzi) [1,2,3]. Approximately 6.5 million individuals globally suffer from CD, with an annual incidence of 30,000 new cases and up to 13,000 deaths per year [1,2,3,4,5,6]. Despite its significant health burden, CD is one of the thirteen most neglected tropical diseases and the third most common parasitic disease worldwide, just after malaria and schistosomiasis [1,5], highlighting the need for policies capable of addressing the public health challenges faced by the affected regions and countries [2,3].
The development of new drug candidates for the treatment of CD constitutes a very active scientific topic among the academic research community [7]; however, to date, only two nonspecific drugs approved about 50 years ago still constitute the unique treatment options: nifurtimox (1940, NFX) and benznidazole (1974, BZN). Despite the efficacy of NFX and BZN in primarily addressing acute stages of the disease, the severity and high number of adverse effects associated with their clinical use significantly limits their pharmacotherapeutic application, ultimately resulting in poor patient adherence and drug resistance.
By the end of the 1970s, a group of cysteine proteases (CP) were identified as crucial components of the T. cruzi biological machinery; with one representative attracting special interest as a potential druggable target. This CP was named Cruzipain (CZP) by Cazzulo et al. [8], and is nowadays known to represent the main CP involved in the T. cruzi lifecycle. Its catalytic domain consists of a single polypeptide chain of 215 amino acid (AA) residues, forming α -helices and antiparallel β -sheets that fold in two distinct subdomains and delineate the active site located within the interface (Figure 1a) [9,10]. It includes a catalytic triad (Cys25, His162 and Asn182) responsible for the enzyme hydrolytic activity [11,12] and four subsites identified as S1’, S1, S2 and S3 [8,9,10] (Figure 1b). Detailed structural and topological exploration has been reported for CZP substrates and inhibitors interacting within these four subsites, with different extents of structure–activity relationships being reported [8,9,10,13,14,15,16,17,18,19,20,21].
A broad diversity of CZP inhibitors has been reported to date, including both non-covalent and covalent inhibitors resulting either from drug repurposing or de novo synthesis [16,17,18,23,24,25,26,27,28,29,30,31,32,33,34,35]. In recent years, increased scientific interest has been observed towards the development of a specific class of covalent inhibitors identified as targeted covalent inhibitors (TCI) [36,37], in which both a maximized complementarity and modulated reactivity with the binding site operate in a synergistic way towards safely and effectively inhibiting the corresponding target. Both of these features can be efficiently optimized by applying modern bioinformatic approaches and technologies related to structure-based drug design (SBDD). Several examples of clinically useful TCI have been reported in the literature [36,37,38,39]; however, to the best of our knowledge, none targeting CZP are currently under clinical development to date, while most of the preclinical candidates lack adequate translation of inhibitory activity from CZP to T. cruzi [40].
Most CZP inhibitors present a peptide-like structure that confers them with unfavorable physicochemical properties and suboptimal biopharmaceutic performance, such as a low oral bioavailability, poor stability and high rate of first pass metabolism [41,42,43]. Particular interest has been placed in strategies involving bioisosteric peptide bond replacement to further enhance these aspects. One example is the report by Brak et al. in 2008 [44], which developed a lead CZP-inhibitor based on the 1,4-disubstituted 1,2,3-triazole scaffold. Although this compound exhibited a high CZP inhibitory potency (Figure 2, Ts-370, IC50CZP = 370 nM), its bioactivity was not efficiently translated into T. cruzi inhibition, as reflected by its 14-fold lower antitrypanosomal potency (IC50Tc = 5.1 μ M). This lack of translation constitutes a recognized hurdle for the development of effective trypanocidal agents and has also been reported for the covalent inhibitor K777 (Figure 2), which stands out for its high CZP inhibitory potency (IC50CZP = 2 nM) but translates into IC50Tc values of ≈10 μ M. Montanari et al. recently reported a carefully designed machine-learning-based study aimed to assess the chemical space shared by candidates combining high CZP and T. cruzi inhibitory activity [40]. Interestingly, the authors concluded that only a small chemical space associated with CZP inhibition intersects with high T. cruzi bioactivity. This limited shared bioactive chemical space is further narrowed if selectivity requirements for CZP inhibition relative to homologous CP, such as human L-cathepsin (hCatL) [10], are taken into consideration. In this context, carefully designed protocols capable of combinatorially screening CZP inhibitors covering a wide chemical space are required for the advancement in the preclinical identification of TCI meriting further progress as antichagasic agents.
In this work, a SBDD approach was applied to combinatorially generate molecular diverse libraries of 1,4-disubstituted 1,2,3-triazole derivatives, further explored by means of molecular modeling (i.e., molecular docking and molecular dynamics) as potential TCI of CZP with maximized complementarity with the corresponding CZP active site. In this search, the atomistic details and known structure–activity relationships with respect to CZP S1’, S1, S2 and S3 subsites were also comprehensively explored.
Overall, a virtual library containing 78,540 diverse drug-like 1,4-disubstituted 1,2,3-triazole analogues was exhaustively studied in silico, further identifying 22 promising candidates exhibiting a binding mode compatible with high CZP inhibitory activity. By applying a versatile Cu(I)-catalyzed azide-alkyne cicloaddition (CuAAC) centered synthetic scheme, the 22 positive binder candidates were obtained and further complemented with a set of 23 moderate to poor performing binders included as inhibitor decoys. The bioactivity of this set of triazole-derivatives was explored in detail, including in vitro measurements of inhibitory activity against CZP and T. cruzi in infection cell assays, as well as their selectivity over hCatL and cytotoxicity over human cells.

2. Results and Discussion

2.1. Virtual High Throughput Screening

2.1.1. Development of an SBDD Screening Workflow

In order to validate the virtual high throughput screening (vHTS) conditions for hit identification, a set of 20 triazole-based CZP inhibitors bearing a¸ tetrafluorophenoxymethyl ketone (4FPMK) warhead (WH) was used as a training set [46] (Figure S1 and Table S1). The reported crystallographic structure of Ts-370 (PDB 3IUT [22]) was considered as reference for the identification of the corresponding bioactive binding poses (see Experimental section for technical details).
Among the diverse docking approaches available, in this study, a classical non-covalent docking protocol was applied, with no bias (i.e., covalent docking or grid based) being included with respect to the orientation of the corresponding WH. Thus, it was possible to evaluate the overall binding pose clusters resulting from the orientation of R1, R2 and R3 within each subsite (S1, S2, S3), further analyzing the resulting reactive pose in the context of the landscape of alternative non-reactive orientations.
The molecular docking search parameters were systematically explored until bioactive poses were successfully identified not only for Ts-370 (Figure S2, RMSD value < 1 Å respect to the crystal) but also for the remaining 19 analogues. The final pose selection criteria to identify hits included a combination of the resulting docking score, cluster size and satisfaction of four H-bond contacts within the CZP active site (i.e., Gln19, Ser61, Gly65, Gly66). As shown in Figure S3a, a clear trend of more negative docking score values was found for the reference compounds exhibiting lower IC50, confirming the predictive power of the molecular docking protocol to identify bioactive poses and sort them according to their inhibitory potency.
To describe the time-dependent conformational evolution of the identified docked poses for the training set, they were subjected to explicit solvent classical molecular dynamics (MM-MD) simulations at 300 K, with free-energy MMPBSA analyses performed on the resulting simulated trajectories. All encounter complexes showed appropriate stability during the simulated time (100 ns); and, as can be seen in Figure S3b, the results obtained also allowed an accurate correlation between the ligand potencies and their predicted Δ G b i n d . Additionally, the electrostatic (EEL) and van der Waals (vdW) interaction patterns analyzed for each ligand–CZP complex (Figures S4 and S5) were in agreement with the reported structure–activity relationships (SAR) of S1’, S2 and S3 subsites.
Correspondingly, the molecular-docking-centered search protocol showed adequate efficiency to screen triazole-based CZP inhibitors, effectively discriminating the most potent candidates.

2.1.2. Construction of the Working Chemical Space

The chemical space associated with the molecular diversity to be explored was based on a synthetic scheme designed to obtain 1,4-disubstituted 1,2,3-triazoles (Figure 3), which allows the combinatorial exploration of R1, R2 and R3 substitutions. In this respect, a Cu(I)-catalyzed 1,3-dipolar cycloaddition (CuAAC [47,48,49,50,51]) was considered as the core of the scheme, allowing the combinatorial coupling of libraries of azides and alkynes.
For the synthesis of azides (Az), an ISA.H2SO4-mediated diazotransfer reaction [52,53] to transform α -amino esters into their corresponding azido-analogues was considered; while for the alkyne building blocks, the approach was based on MW-assisted A3 coupling reactions [54,55,56], resulting in a set of terminal propargylamines (PA). Additionally, two alternative WH groups were explored: esters (Es) and aldehydes (Ald), with the former being obtained from α -azido esters via CuAAC, while the latter were derived from the resulting Es by applying a DIBAL-based reduction [57]. Further details regarding the corresponding synthetic procedures are provided in Section 2.2 (Synthesis of selected candidates and biological evaluation against CZP), Section 3.2 (Synthetic procedures) and Supplementary Materials.
Once the synthetic protocols were defined, and considering the great diversity of potential starting materials that could be employed, massive virtual screening libraries were constructed by enumerating the presented chemical reactions and evaluated using the vHTS method developed for the identification of promising candidates. These libraries were built and managed using TidyScreen (https://github.com/alfredoq/TidyScreen, accessed on 1 March 2024), an open-source Python-based in-house developed framework that leverages a structured SQL database to streamline the organization, combinatorial synthesis and molecular docking execution of large screening projects (see Section 3.1.1 for further details). Thus, the database of organic molecules eMolecules [58] was used as a building block source for the combinatorial construction of massive virtual libraries of triazole derivatives. In order to narrow down the size of the libraries to drug-like candidates, those exhibiting favorable properties (i.e., MW < 700 Da, ClogP between 0 and 6, less than 5 HBD and 10 HBA groups, chiral centers < 3, and a max. of 10 rotatable bonds) were considered for further exploration.
Overall, a set of 1517 AA, 6022 aldehydes and 4037 amines were retained for in silico synthesis. Considering the large number of building block combinations (>300 billion possibilities), a strategy applying an initial stepwise comprehensive regional virtual screening to identify the most promising R1, R2 and R3 substituents preceded a final combinatorial analysis of the passing results, as described in detail in the following sections.

2.1.3. Virtual Screening on R 1

The chemical space on R1 was screened in silico evaluating a set of 1517 AA, while keeping constant the cyclopentyl, pyrazolopyrimidine and 4FPMK scaffolds of Ts-3 as R2, R3 and WH, respectively. The set of derivatives generated were subjected to molecular docking studies, with results being classified based on the hit identification criteria defined as part of the protocol validation stage (i.e., lowest energy and most populated cluster, complying with the pharmacophoric intermolecular interactions). This selection procedure is graphically presented in (Figure 4), with the upper left region corresponding to derivatives identified as preferred R1 substituents. From the 35 promising R1 (Figure S6), it is worth highlighting the identification of Arg, Lys, His and Phe as promising starting materials, exhibiting even lower docked energies compared to the reference compound Ts-3 (Figure 4), consistent with the known preference of CZP for basic and/or aromatic AA side chains oriented towards S1. In addition, a wide range of physicochemical properties is observed, with scaffolds ranging from totally apolar to ionizable (cationic); highlighting the chemical space versatility that can be explored in the search for adequate biological activity translation.
As was found for the training set during the vHTS protocol validation, most of the ligands oriented their R1 substituent towards the bulk of the solvent. A notable exception was the ligand derived from Lys (R1Lys-4FPMK, Figure S7a), which exhibited more favorable docking energies compared to the rest of the candidates (−16.64 kcal·mol−1 vs. −12.17 kcal·mol−1 for the second with the lowest energy, Figure 4), placing it as the top-ranked derivative resulting from the docking campaign. Visual inspection of its binding mode evidenced an extra contact point not typically identified for compounds in the training set. Specifically, the methylcarboxylate side chain of Asp161, located at the S1’/S2 interface of the CZP binding site, is able to interact with the ionized amine moiety of the Lys side chain strongly favored by the flexibility, size and length of the n-butylamino group (Figure S7b).

2.1.4. Virtual Screening on R2

The molecular diversity generated on the R2 substitution resulted from 6022 selected aldehydes, which were screened by maintaining the R1, R3 and WH moieties of Ts-3. Upon construction of the library and applying the docking procedures, 51 promising substructures were identified (Figures S8 and S9). Among the privileged R2 moieties, two main groups can be distinguished: aromatic and aliphatic substituents, with a reduced amount of heteroatomic scaffolds, in accordance with the highly recognized hydrophobicity of S2 [8,13,14,15,16,17,18,19]. It is worth noting the outstanding result achieved with the derivative bearing a cyclohexyl group, not only in terms of energy but also in the population of the docked cluster (cHex, Figure S9).

2.1.5. Virtual Screening on R 3

Screening on R3 was accomplished by considering 4037 amine building blocks as alternatives of the pyrazolopyrimidin-5-ylmethanamine moiety present in Ts-3, resulting in a total of 22 scaffolds being identified as privileged based on their complementarity within the S3 subsite (Figures S10 and S11). Most of them exhibit at least one H-bond interaction with the hydroxymethyl side chain of Ser61, while no aromatic amine passed the bioactive pose filtering criteria. These observations highlight the role of Ser61 in ligand binding, emphasizing the need for some degree of flexibility to optimize this interaction. The findings are consistent with the already reported preference of S3 for bulky and/or positively charged aromatic groups [8,10,17,20], which confers binding selectivity for CZP over human-cathepsins [17,18]. It is worth pointing out the remarkable molecular docking performance of the derivative bearing a benzodioxole group as R3 (BDX, Figure S11), complying with the requirements of bulkiness and H-bond interaction with Ser61, which, to date, has not been reported in CZP inhibitors.

2.1.6. Combinatorial Screening of Favored R 1 , R 2 and R 3 with Inclusion of Selected WH

In a last step, a virtual combinatorial library resulting from 35 AA (R1), 51 aldehydes (R2), 22 amines (R3) and 2 reference WH (aldehyde, Ald and methylester, Es) was constructed to systematically explore the best overall combination of substituents, leading to a library of 78,540 candidates. Upon applying the docking analyses, as expected, most of them passed the interaction-based post-docking hit-selection criteria (n = 77,024) with medium to high cluster sizes (31–70%, Figure 5).
Remarkably, about 65% of the new candidates generated from the combinatorial synthesis of promising fragments showed more favorable docking energies (−25 to −17 kcal·mol−1) than both the training set (Section 2.1.1, −6.9 to −11.3 kcal·mol−1) and those compounds obtained during the regional search on R1, R2 and R3 (Section 2.1.3, Section 2.1.4 and Section 2.1.5; −17 to −5 kcal·mol−1). A certain preference could be identified for derivatives bearing cationic, polar and aromatic R1 groups (Lys, Arg, Ser, Phe, Tyr and His), combined with cHex as R2 and BDX as R3 (Figure 5a, green box, top left). The identification of Lys derivatives (i.e., Ald-6 and Es-6) as promising candidates, outstanding over their analogues mainly with respect to their docking energies, should once again be underlined. Additionally, no significant differences in the resulting docked energies were observed between compounds bearing an Ald or Es as WH; suggesting that any difference in inhibitory activity may be governed by the reactivity of the corresponding WH.
Hence, considering all the commercially available starting materials included as part of the screening campaign, a chemical workspace (n = 77,024, Figure S12) was identified in silico to obtain triazole derivatives with good drug-likeness and the potential to act as antichagasic agents; leading to synthetic efforts as disclosed in the following section.

2.2. Synthesis of Selected Candidates and Biological Evaluation against CZP

The synthetic scheme shown in Figure 6 was followed to obtain the corresponding triazole derivatives.
Based on the results obtained from the in silico screening, the main focus was placed on varying the R1 substituents, using AA classified as optimal (n = 6; Arg, Lys, His, Ser, Tyr, Phe), regular (n = 5; Cys, Met, Trp, Leu and Val) and suboptimal (n = 4; Ala, Thr, Glu, Phg). The first two groups successfully passed the massive virtual screening stages, but with significantly different docking energy and cluster size values, while the AA classified as suboptimal did not overcome any stage of the vHTS. In addition, the c-Hex scaffold was kept constant, since it was identified as a highly privileged R2 fragment, while two alternative moieties were used to decorate R3: BDX (capable of establishing H-bond interactions with Ser61) and THQ (with no HBA capability).
As stated in the previous section (Figure 3), the 15 azide building blocks (Az-1 to Az-15) were obtained via an ISA.H2SO4-mediated diazotransfer reaction adapting the previously reported conditions [52,53], with yields ranging from 16 to 95% (see Section 3.2.2 for further details). PA-1 was obtained by a copper(I)-catalyzed microwave-assisted A3-coupling reaction [54] followed by KHO-mediated desilylation in methanol. In contrast, the copper(I)-catalyzed decarboxylative A3 coupling reaction reported by Ermolat’ev et al. [55] was employed to obtain PA-2, allowing the terminal propargylamine to be yielded in one single step. Details of the reaction conditions and corresponding characterization of the building blocks are described in Section 3.2 Synthetic procedures and Supplementary Materials (pp. 24–28).
Overall, a set of 45 chemically diverse 1,4-disubstituted 1,2,3-triazoles were synthesized via CuAAC (Figure S12) and further subjected to CZP bioassays, as presented in Table 1.
As can be seen, good to excellent synthetic yields (28–90%) were achieved for all combinations of R1 moieties with c-Hex and BDX or THQ (i.e., PA-1 and PA-2), including an Es WH, resulting in 24 triazole derivatives (Es-1 to Es-24). The bioactivity screening against CZP was performed at 100 μ M using E-64 as positive control, a recognized specific covalent irreversible inhibitor of cystein proteases [59]. As expected, E-64 inhibited CZP in 98% (±2%) at 1 μ M after 10 min of preincubation, thus validating the experimental setup.
Among the 24 Es derivatives evaluated, Es-10, Es-14 and Es-15 inhibited the enzyme in at least 47%; with Es-15 excelling for its IC50CZP of 6.8 ± 0.3 μ M. The target-specific mechanism of Es-15 was confirmed upon observing unchanged CZP inhibition percentages at various concentrations of Triton X-100 and after preincubation with BSA, thus discarding colloidal aggregation [60,61,62] (Table S2). Moreover, a time-dependent inhibitory mechanism was also ruled out for Es-15 after observing unaltered percentages of CZP activity inhibition upon assays without pre-incubation (81 ± 2%, Table S3). Thus, to determine the inhibitory mechanism of Es-15, enzyme initial velocities at different substrate and compound concentrations were monitored. From the Michaelis–Menten graph (Figure 7a), the V m a x , K M a p p and α values were determined. An α value equal to 8.3 indicated a competitive mechanism of inhibition; while a positive slope was determined upon evaluating the K M a p p at varying compound concentrations (Figure 7b). The Lineweaver–Burk plot reveals that V m a x decreases when Es-15 concentration increases, as an increment in the intercept value on the ordinate can be observed (Figure 7c). Therefore, the increased K M a p p and decreased V m a x values upon higher Es-15 concentration suggest a mixed mechanism of inhibition, with a calculated K i of 3 ± 1 μ M (general linear F-test, p-value = 0.0127).
The aforementioned observations suggest that the ester group is not acting as an electrophilic center capable of being nucleophilically attacked by Cys25, which might be the main reason for the lack of inhibitory potency observed for the Es derivatives, including R1, identified as highly favored by in silico screening (Figure 5).
The corresponding Ald-1 to Ald-24 derivatives were synthesized from their ester-based precursors with fair to good yields (Table 1). It was not possible to synthesize Ald analogues of Es-14, Es-15 and Es-24 by applying the DIBAL-mediated reduction, since the crude reaction mixtures were found to be riddled with by-products. Among the 13 compounds bearing BDX as R3 (Ald-1 to Ald-13), those bearing R1 scaffolds identified as optimal (Ald-1, Ald-4, Ald-6, Ald-10, Ald-13) inhibited CZP in a range of 64 to 87% at 100 μ M; while those derived from moieties predicted as regular (Ald-2, Ald-5, Ald-7, Ald-12) or suboptimal (Ald-3, Ald-8, Ald-9, Ald-11) R1 inhibited the enzyme activity to a lower extent (8 to 64%). Among the best performers, Ald-6 (R1: Lys) and Ald-10 (R1: Ser) stood out as the most potent ones, with IC50 of 3.3 ± 0.3 μ M and 3.4 ± 0.5 μ M, respectively. Correspondingly, most Ald-based triazoles bearing THQ as R3 (Ald-16 to Ald-23) exhibited inhibition values ranging from 8% to 55% at 100 μ M and were considered inactive, with the exception of Ald-20 (R1: Ser), which presented 83% of CZP inhibition and an IC50 of 7 ± 4 μ M. From the results stated, the marked difference in inhibitory activity of BDX derivatives with respect to THQ-derived ones reinforces the relevance of Ser61 as a hotspot at S3, favoring appropriate CZP recognition.
After ruling out nonspecific enzyme inhibition caused by colloidal aggregation (Table S2), a time-dependent mechanism was verified for Ald-6 and Ald-10 based on the significant reduction in their inhibitory percentages without prior incubation with CZP (Table S3). Their modes of action were further examined through jump dilution assays (JDA, Figure 8a,b) [63,64,65,66].
From these studies, it was notable that Ald-6 and Ald-10 progress curves exhibited slopes equal to about 12% and 24% of that of the control (DMSO), respectively, showing a homologous behavior to that observed for the widely known covalent irreversible inhibitor E-64, which selectively targets Cys25 [63]. Altogether, the presented mechanistic studies for Ald-6 and Ald-10 strongly support the involvement of a covalent irreversible inhibition of the catalytic residue of CZP, which is in agreement with reports on tight-binders bearing an aldehyde WH [67,68].

2.3. Structure–Activity Relationships (SAR)

From the analysis of the biological and computational results for the set of triazole derivatives evaluated, some conclusions regarding their SAR are worth noting.
The first very important conclusion refers to the efficacy of the bioisosteric replacement of the amide bond by the 1,4-disubstituted 1,2,3-triazole group. In particular, the triazole ring effectively mimics the classical peptide bond of the CZP substrates, complying with stable H-bond interactions with Gly65 and Gly66, as observed in MD simulations (Figure S13); which are confirmed as pharmacophoric contacts.
Although the involvement of the oxyanion hole and the catalytic triad for the mechanistic processes of CZP is well known, their relevance for the formation of the encounter complexes becomes evident here. Gln19 also appears as a pharmacophoric contact, acting as both an H-bond acceptor and donor residue whose electrostatic contribution to the stabilization of these compounds is remarkable (Figure S13). In fact, together with Ser61 (S3), it is responsible for the low inhibitory activity found for Ald derivatives containing THQ as R1.
Similarly, for Ald-10, one of the most potent covalent inhibitors identified (IC50: 3.4 ± 0.5 μ M), an interaction pattern characterized by contact with His162, a member of the CZP catalytic triad, was found. Although the R1 hydroxymethyl substituent of this Ser derivative is small and lacks basic character, its interaction with the imidazolium group of His162 strongly favors the formation of the encounter complex, while the analysis of the docking and MM-MD results revealed a stable H-bond interaction throughout the simulation between the -OH of the R1 and His162 residue (Figures S13 and S14). This interaction also seems to favor the stabilization of the corresponding WH by orienting and positioning it closer to the catalytic AA Cys25. Indeed, despite not possessing a favorable anchor point at R3, Ald-20, also derived from Ser at R1, showed moderate inhibitory activity (IC50: 7 ± 4 μ M). Although these results are not evident from the computed Δ G b i n d values (Figure S15), the preferential orientation and stabilization of the WH could imply an improvement in ligand reactivity, an effect that cannot be assessed by classical computational design strategies. This hypothesis is also supported by the null activity of Ald-2, Ald-8 and Ald-11, three Ald-10 analogues of similar size but lacking relevant contacts in R1 with His162.
Furthermore, Ser64 and Asp161 were identified as being responsible for the CZP preference for inhibitors bearing basic and/or aromatic groups oriented towards S1/S1’, with three compounds including a basic side chain on R1 among the most active ones.
The docking pose and dynamic behavior of Ald-6 (IC50: 3.3 ± 0.3 μ M, Figure 9), synthesized from the AA Lys to introduce the R1 substituent, always excelled over the rest of the analogues, with the flexible and elongated n-butylamino group able to optimize electrostatic interactions with Asp161 (S1’) when complexed to CZP (Figure S13). As discussed above, this AA is located at a distant edge of the binding site. This encouraging result was also reflected in classical molecular dynamics studies complemented by binding free energy analysis, as this complex stood out from the other aldehyde derivatives due to its favorable Δ G b i n d , with a difference of about 7 kcal·mol−1 compared to the second best Ald (Figure S15), mainly attributed to an EEL contribution.
The importance of Asp161 is also reflected for Es-15 (Figure S16) and Ald-4 (Figure S17), since the encounter complex formation during the molecular docking studies is favored by the presence of this anionic AA. However, the larger size and lower flexibility of the propylguanidinium group of Es-15 and the lower extension of the methylimidazole of Ald-4 prevent this contact from remaining constant over time, as observed in the MM-MD. The presence of an additional contact point with Ser64 seems to be responsible for the high inhibitory activity of Es-15 (Figure S13), leading to an optimized encounter complex that results in the non-covalent inhibition of the enzyme. On the other hand, inhibitors including an acid moiety as R1 were not active at all, with Ald-3 constituting a representative example, highlighting the role of Asp161 for ligand recognition.
Finally, Phe, Trp and Tyr derivatives (Ald-1, Ald-12 and Ald-13) also showed promising results, although with potencies 12.4, 13.8 and 3.9 times lower than the most active analogues (Ald-6 and Ald-10). The size and bulk of their R1 groups seems sufficient to establish a possible π –anion interaction with the Asp161 carboxylate [69] (Figure S18a), although the energetic contribution of this contact is not accounted in the computed docked energies since it is not parameterized in the scoring function used. More detailed tests involving quantum studies and/or hybrid QM/MM-MD studies are required to further explore this hypothesis. In the particular case of Ald-13, its higher potency compared to the other two R1-aromatic derivatives seems to be justified by the additional anchor point with Ser64, which also reinforces the role of this AA for the stabilization of CZP:inhibitor complexes within S1 (Figure S18b).
When the bioactivity of the studied compounds was plotted in the context of the explored chemical space [70] (Figure 10), one main region containing the highest proportion of active CZP inhibitors was found. As can be seen in its magnification (Figure 10, right), a significant number of combinations with expected high potency have not been synthetically produced as part of this initial research effort, although the methodological implementation should be straightforward based on the principle of similarity [70,71,72,73,74].

2.4. Selectivity Assays: Screening against hCatL

In view of the homology between CZP and hCatL, it is crucial to evaluate the selectivity of CZP inhibitors with respect to the potential unwanted hCatL inhibition, since this may directly affect the in vivo safety profile of antichagasic agents. In this respect, the inhibitory activity against hCatL of the most potent triazole-based CZP inhibitors was evaluated at 100 μ M in triplicate. As for the CZP assays, E-64 was used as positive control.
All the evaluated compounds inhibited hCatL activity in less than 26% at 100 μ M (Table 2). Therefore, Es-14, Es-15, Ald-1, Ald-4, Ald-5, Ald-6, Ald-10, Ald-12, Ald-13 and Ald-20 can be considered selective CZP inhibitors.

2.5. Anti T. cruzi Bioactivity

Given the significance of attaining a good translation of the CZP inhibitory potency to the corresponding trypanocidal effect, the most potent CZP inhibitors were evaluated on a T. cruzi infection model in Vero cells, allowing the analysis against both the trypomastigote (extracellular) and amastigote (intracellular) forms. In addition, to further assess the safety of these compounds, cell viability was also measured. Initially, all inhibitors were tested in triplicate at 50 μ M, quantifying the percentage of remaining infection after measuring β -galactosidase activity (Figure 11).
Es-15, Ald-4, Ald-6, Ald-10 and Ald-20 exhibited acceptable anti-T. cruzi/cytotoxicity balance; among which Ald-6 and Ald-10 stood out for reducing 97% of infection at the evaluated concentration (Figure 11a). Interestingly, these derivatives also exhibited the highest CZP inhibition, specifically those derived from Lys (Ald-6) and Ser (Ald-10) as R1 substituents, as predicted by the in silico studies. These two derivatives also showed low cytotoxicity (Figure 11b) in contrast to the behavior observed for Es-14 and Ald-2, which, despite their high inhibitory potency, showed significantly reduced cell viability.
In a second screening stage, Ald-6 and Ald-10 were further assayed at 10 μ M. Ald-10 only reduced T. cruzi viability by 38%, while Ald-6 excelled for its excellent antitrypanosomal activity, with only 7% of parasite survival after 72 h of incubation.
In the final evaluation stage, Ald-6 was assayed at 1 μ M, finding a 27% reduction in T. cruzi infection; calculating an IC50Tc of 3.2 ± 0.4 μ M (Figure S20). Notably, the mean IC50 values determined for both CZP (IC50CZP = 3.3 ± 0.1 μ M) and T. cruzi are statistically equivalent (ANOVA P-value: 0.6563), evidencing a complete bioactivity translation, with no biopharmaceutical limitations at the cellular level. In this respect, Ald-6 exhibited a translation ratio close to 1 (IC50Tc/IC50CZP ratio: 0.97 ± 0.16) and, as observed in Figure 10, it occupies a different region of the computed chemical space compared to the rest of triazole derivatives assayed, thus defining a region of particular interest for the further design of analogues.

3. Materials and Methods

3.1. Computer-Aided Studies

3.1.1. Construction of Chemical Libraries

The workflow for the generation of three-dimensional structures of each compound to be evaluated is initiated with a list of SMILES accompanied by an ID or identification code. The SMILES of the 20 triazole derivatives included in the test set were generated using MarvinSketch [75]. Massive libraries of potential new inhibitors were generated employing the Python-based package TidyScreen (https://github.com/alfredoq/TidyScreen, accessed on 1 March 2024). Initially, the SMILES lists of the corresponding building blocks (aldehydes, amines and AA) were extracted from eMolecules [58]. By using the RDKit module [76], the lists were filtered according to criteria based on drug-like properties and excluding those compounds containing metals, non-biocompatible elements or those capable of interfering with the synthetic procedures, such as duplicated reactive groups leading to side products in the different stages of synthesis. In this way, the lists were used to synthesize the triazole derivatives employing the RDKit’s Chem module by means of SMARTS-based reactions [76], combining the different filtered substances by replicating the planned reactions to be executed for their synthesis (diazotransfer, A3 coupling, CuAAC and DIBAL reduction). Once the triazoles were generated, they were subjected to a second round of drug-like-based filtration. From the resulting lists, the chiral centers present in each compound were determined, with the corresponding diastereomers being computed. In-house developed scripts allowed the selection of only those isomers of interest or, alternatively, all of them in case no stereoselectivity was restricted. In the case of the test set, only those exhibiting an (S,S) configuration were retained. Considering the planned synthetic protocol, the conformation of the chiral α -centre at the N1-triazole position was kept constant, according to the configuration of the α -AA used for their synthesis, varying only that of the C α at C4. Subsequently, using Dimorphite-DL 1.3 [77], the ionization state of each ligand was evaluated on the basis of an estimated pKa and a defined pH range. Based on this criterion, pH = 5.5–7.4 was set, and all predominant ionization equilibrium forms in this range were retained and considered. Next, RDKit was used to generate a maximum of 100 different conformations, further applying a minimization protocol by considering the Merck Molecular Force Field (MMFF94) [78]. The resulting conformers were ranked on the basis of their computed energy, retaining the lowest energy 3D conformation to initiate ligand preparation procedures. Finally, each ligand was parameterized, generating a coordinate file in the PDBQT format, as required by the docking engine employed (AutoDock suite) [79]. The database containing the building blocks and triazole derivatives evaluated are available for download at https://github.com/alfredoq/Triazole_based_TCI_2024_paper (accessed on 14 March 2024), and a step-by-step tutorial to build and reproduce the libraries is detailed in section “Preparation of a Virtual Chemical Library for vHTS campaigns” at https://github.com/alfredoq/TidyScreen (accessed on 1 March 2024).

3.1.2. CZP Three-Dimensional Structure

The CZP structure deposited under the code PDB ID 3IUT [22] was used as receptor template and was subjected to optimization procedures and further parameterized for molecular docking procedures. In this respect, in a first instance, the missing disulfide bonds were added (Cys22-Cys63, Cys56-Cys101 and Cys155-Cys203), and the structures of the ligands and additives present in the crystallographies were deleted. Considering that many studies claim that the Cys25-His162 pair of the catalytic triad would be mostly in their ionized states, both structures were parameterized so as to generate the corresponding ion pair Cys25S-His162H+. The whole preprocessing was performed using ADT v4.2 [79], which also included the addition of polar hydrogen atoms and assignment of Gasteiger charges to the macromolecular structure. The processed protein structure was finally saved as a PDBQT file and stored for later use.

3.1.3. Molecular Docking

The AutoDockGPU 4.2.6 [80] software was employed to accomplish molecular docking studies, running 100 docking runs using a Lamarckian genetic algorithm (LGA) to sample the ligands’ conformational space. The corresponding affinity grids were generated using a 40 × 42 × 60 Å box centered on Cys25 and further extended in order to include the entire active site of CZP. A grid spacing of 0.375 Å was used to calculate the affinity points. The docking results were filtered by extracting the lowest energy and highest cluster enrichment docking poses satisfying the considered pharmacophoric contacts: H-bond interactions (default distance cutoff: 3.7 Å) with Gln19, Ser61, Gly65 and Gly66. Clustering of the resulting docked poses was performed using the built-in functionality of AutoDock-GPU, allowing an RMSD tolerance of 2.0 Å, handling symmetry in the corresponding calculation and saving all computed clusters.

3.1.4. MM-MD and Free Energy Analysis

Molecular dynamic (MM-MD) simulations were performed under explicit solvent conditions using the AMBER22 software package. Ligand and receptor parameters were assigned from the GAFF2 and ff14SB force fields, respectively; with the system preparation and parametrization being performed using the tLeap module of AMBER22 package [81]. The simulation trajectories were obtained by applying periodic boundary conditions, constructing octagonal boxes of explicitly pre-equilibrated TIP3P water molecules, with a minimum distance of 16 Å between the edges of the box and the solute. To obtain the corresponding trajectories, a previously validated MD workflow was applied using the AMBER22 pmemd.cuda module. First, a double-stage minimization step, in which the water molecules were initially minimized with the restricted solute, followed by a minimization of the whole system (10,000 steps for each stage), was carried out. Next, a heating phase was applied, thermostatizing the system from 0 to 300 K for 500 ps, followed by an equilibration phase for an additional 5 ns. Finally, the equilibrated systems were subjected to production runs (100 ns) obtained under constant pressure and temperature conditions, using a time step of 2 fs. The convergence of the simulation was checked by structural and energetic inspection using the cpptraj module of AMBER22, complemented with in-house Python scripts developed exclusively for this purpose. Afterwards, global and per-residue energy analyses of each MM-MD simulation were performed using the MMPBSA.py script [82,83,84], which were applied over the entire trajectory (100 ns), with individual snapshots sampled every 10 frames.

3.1.5. Molecular Visualization

Snapshot visualization of docking poses, structural properties and protein–ligand complexes were calculated and represented using VMD v1.9 software [85]. Most of the 3D-depictions included throughout the document were also rendered from VMD. Schematic diagrams of intermolecular interactions were generated from the results obtained with ProLIF [86] and LigPlot [87].

3.1.6. Computational Infrastructure

All the in silico studies and analyses were carried out using computational resources provided by the MedChem Lab, Departamento de Ciencias Farmacéuticas, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba (FCQ-UNC, Córdoba, Argentina) headed by Prof. Mario A. Quevedo; and at the Centro de Computación de Alto Desempeño (CCAD, UNC [88]), Sistema Nacional de Computación de Alto Desempeño, of the Ministerio de Ciencia, Tecnología e Innovación (SNCAD, MinCyT), particularly in the Mendieta and Serafín clusters.

3.2. Synthetic Procedures

3.2.1. Chemistry

The synthetic procedures were carried out at the Laboratory of Organic Synthesis, Sustainable Chemistry for Metals and Molecules Division, Department of Chemistry, KU Leuven (Belgium), headed by Dr Wim Dehaen. All chemicals were purchased from Acros Organics (Geel, Belgium), Alfa Aesar (Karlsruhe, Germany), Fluorochem (Glossop, UK), Merck (Darmstadt, Germany) and TCI Europe (Zwijndrecht, Belgium) and used as is without further purification. Dry diethyl ether, THF and toluene were obtained using an M-Braun SPS-800 system. Moisture sensitive reactions were carried out under nitrogen or argon atmosphere using flame dried glassware. A CEM-Discover microwave reactor was employed for the MW-assisted procedures. Reaction conversion was monitored via TLC analysis using MilliporeSigma™ Silica Gel 60 F254 Coated Aluminum-Backed TLC Sheets (Sigma, St. Louis, MO, USA). Column chromatography was performed via standard column chromatography, employing 70–230 mesh silica 60 (Acros) as the stationary phase. 1H NMR spectra were recorded on Bruker Avance 300 (300 MHz) or Bruker Avance 400 (400 MHz) spectrometers. The 13C NMR spectra were recorded on a Bruker Avance 400 (101 MHz working frequency, Bruker, Billerica, MA, USA). The NMR samples were dissolved in a corresponding deuterated solvent, and chemical shifts ( δ ) were reported in parts per million (ppm) and referenced relative to the deuterated solvent signal [89]. The FTIR-spectra were recorded on a Bruker Vertex 70 spectrometer (Bruker, Billerica, MA, USA), and OPUS 8.5 software was used to analyze the recorded spectra. The high-resolution mass spectra were acquired on a quadrupole orthogonal acceleration time-of-flight mass spectrometer (Synapt G2 HDMS, Waters, Milford, MA, USA). The samples were infused at 3 μ L/min and spectra were obtained in positive or negative ionization mode with a resolution of 15,000 (FWHM—full width at half maximum) using leucine enkephalin as a lock mass. The full characterization data of intermediates and evaluated compounds are available within the Supplementary Materials.

3.2.2. Azides (Az)

Imidazole-1-sulfonyl azide hydrogen sulfate (ISA·H2SO4) was synthesized according to Method B reported by Potter et al. [52] To a cooled and stirred suspension of sodium azide (20 mmol) in dry EtOAc (20 mL) under a N2 atmosphere, sulfuryl chloride (20 mmol) was added dropwise over 10min, after which the mixture was warmed to room temperature and stirred for 17 h. After recooling at 0 °C, imidazole (40 mmol) was added. The thick suspension was stirred for 3 h at 0 °C and then basified with saturated NaHCO3. The organic layer was washed, dried, cooled (0 °C) and placed under a N2 atmosphere. H2SO4 (conc., 20 mmol) was added dropwise over the course of 5 min, after which the resulting mixture was gradually warmed to room temperature. After 45 min, the precipitate formed was filtered and dried to yield pure ISA.H2SO4 (16.5 mmol, 4.5 g, 83%).
General procedure A for the synthesis of α -azido esters (Az-1 to Az-15). This method was adapted from the one reported by Goddard-Borger and Stick [53]. To a round-bottom flask equipped with a magnetic stirring bar, the corresponding α -amino ester (1 eq), K2CO3 (3.4 eq), CuSO4·5H2O (0.01 eq) and ISA.H2SO4 (1.2 eq) were added and suspended in methanol (5 mL). The reaction mixture was stirred overnight at room temperature for the specified time (14–15 h) and concentrated under reduced pressure. The crude was resuspended in EtOAc and extracted twice with 1 M HCl, washed with water, dried over MgSO4, filtered and concentrated in vacuo. Flash chromatography was carried out when necessary to afford the desired product.

3.2.3. CuAAC

General procedure B for the synthesis of 1,4-disubstituted 1,2,3-triazole derivatives (Es-1 to Es-24). This method was adapted from the one reported by Zsabka et al. [90] In a round-bottom flask equipped with a magnetic stirring bar, the corresponding alkyne (1 eq) and azide (1 eq) were dissolved in DCM (0.05 M). A solution of sodium ascorbate (2 eq) in H2O (0.1 M) was added as a second phase; followed by the dropwise addition of CuSO4·5H2O (0.02 eq) in H2O (5 mM) and DIPEA (3.5 eq). The resulting mixture was stirred vigorously at room temperature until full conversion, as indicated by TLC. The resulting mixture was extracted three times with H2O:CHCl3, and the combined organic layers were dried over MgSO4. The solvents were evaporated under reduced pressure, and the residue was purified (when necessary) to afford the corresponding 1,4-disubstituted 1,2,3-triazole derivative.

3.2.4. Aldehydes

General procedure C for the synthesis of aldehyde derivatives (Ald 1-23). This method was adapted from the one reported by Wood et al. [57]. To an oven-dried round-bottom flask equipped with a magnetic stirring bar, the ester derivative (1 eq) was added and dissolved in dry THF (0.2 M) under inert atmosphere. The solution was cooled in a dry ice/acetone cooling bath, followed by the dropwise addition of diisobutylaluminium hydride (DIBAL, 3 eq) in THF. After full conversion, as indicated by TLC (2–5 h), the temperature was raised to 0 °C, and the reaction was quenched by the slow addition of acetic acid (15 eq). A saturated aqueous potassium sodium tartrate solution was added, and the reaction mixture was stirred for 20 min. The resulting mixture was extracted five times with EtOAc, and the combined organic layers were dried over MgSO4. The solvents were evaporated under reduced pressure, and the residue was purified by column chromatography to afford the corresponding aldehyde derivative.

3.3. Biological Assays

3.3.1. Screening against CZP

The previously reported methods were followed to evaluate the in vitro activity against CZP [91,92,93,94]. Allison Doak and Prof. Brian Shoichet (University of California San Francisco, San Francisco, CA, USA) generously provided recombinant CZP (cruzain). The enzyme activity was measured by monitoring the cleavage of the fluorogenic substrate Z-FR-AMC at 25 °C. Unless stated otherwise, assays were performed using 2.5 μ M of substrate ( K M = 0.5 ± 0.1 μ M) and 0.5 nM CZP in a 0.1 M sodium acetate buffer, pH 5.5, 0.01% Triton X-100 and 2 mM β -mercaptoethanol. DMSO and 1 μ M E-64 were employed as negative and positive controls, respectively, in all assays. Fluorescence was monitored at 340/440 nm (excitation/emission) over time in a Biotec 87 Synergy 2 fluorimeter at the Multiuser Laboratory of the Biochemistry and Immunology Department (UFMG, Belo Horizonte, Brazil). All assays were performed in triplicate in two independent experiments.
Initial screening. For the initial screening against CZP, compounds at 100 μ M were tested against CZP with and without a 10min pre-incubation with the enzyme. The fluorescence was monitored for 10 min at 5 s intervals. DMSO and E-64 were employed in all assays as negative and positive controls, respectively. Initial velocities of the negative control (DMSO) were compared to those in the presence of compounds to calculate the percentage of inhibition. The reported values correspond to the mean and standard error of the mean (SEM) of experiments performed in triplicate in two independent experiments.
IC50 determination. Compounds in at least eight distinct compound concentrations (400 μ M, 100 μ M, 25 μ M, 6.25 μ M, 1.5 μ M, 0.4 μ M, 0.1 μ M and 0.025 μ M) were incubated with CZP for 10 min, and once the substrate was added, the fluorescence was monitored for 10 min at 5 s intervals. IC50 curves were obtained from a non-linear fit of the inhibition values versus the logarithmic concentration of the compounds. The reported IC50 values are the mean and SEM of two independent measurements performed in triplicate.
Aggregation Assays. To evaluate detergent sensitivity, CZP was pre-incubated with compounds at their IC50 for 10 min at room temperature in varying concentrations of Triton X-100 (0%, 0.1% and 0.01%). Additionally, compounds at their IC50 were incubated for 10min with BSA (4 mg/mL, Sigma-Aldrich, St. Louis, MO, USA) followed by 10 min incubation in the presence of CZP, and upon the addition of substrate, the enzyme activity was monitored. For comparison, the four conditions were assayed in parallel. Aggregation was excluded when the difference in inhibition percentages between the different conditions was less than 20%.
Jump dilution assay. For CZP at 100-fold, its final assay concentration was incubated with the inhibitor(s) at 10-fold its respective IC50 value for 30 min in a volume of 2 μ L in a 96-well plate. This mixture was diluted 100-fold with an assay buffer containing a 5 μ M Z-FR-AMC to a final volume of 200 μ L, resulting in a standard concentration of enzyme and 0.1 times the IC50 value of inhibitor. Fluorescence intensities were monitored continuously for 2 h.
Mechanism of inhibition. To determine the mechanism of inhibition and K i values of Es-15, CZP activity was monitored at a minimum of six different substrate concentrations (0.3–20.0 μ M) and five inhibitor concentrations (from 0.25× to 2× IC50) and in the absence of the compound. The results were calculated by nonlinear (Michaelis–Menten) and linear regression (Lineweaver–Burk plots). The effect of the inhibitor concentration on the K M a p p and its K i was evaluated by the general linear F-test.

3.3.2. CatL Screening

The recombinant human cathepsin L (hCatL) (Sigma-Aldrich) was activated prior to the assays. First, 10 ng/ μ L of the enzyme was incubated in the activation buffer (50 mM MES, pH 6, 1 mM EDTA, 5mM DTT) for 30 min at 37 °C. hCatL activity was measured based on the cleavage of the fluorescent substrate Z-FR-AMC in a Biotec 87 Synergy 2 fluorimeter at the Multiuser Laboratory of the Biochemistry and Immunology Department (UFMG, Brazil). The assays were executed in 96-well round bottom black plates in 50 mM MES buffer, pH 6, containing 1 mM EDTA, 5 mM DTT, 0.001% Triton-X100, 5 pg/ μ L enzyme and 25 μ M Z-FR-AMC at 37 °C. The triazole-based derivatives were tested against hCatL, at 100 μ M, with a 10 min pre-incubation with the enzyme. The initial velocities of the DMSO control were compared to those in the presence of the compounds. Assays were executed in triplicate in two independent experiments, and enzyme activity was measured for 20 min.

3.3.3. T. cruzi Inhibition Assay

To perform the assay, the drug screening method reported by Buckner et al. was adapted [95,96]. Vero cells were seeded at 1 × 104 per well on 96-well culture plates in 100 μ L MEM-5% FBS. Tulahuen trypomastigotes overexpressing the E. coli β -galactosidase protein were distributed in triplicates on the plate. After a 2 h incubation at 37 °C, the culture media were extracted and replaced with fresh MEM-5% FBS with 1 μ M, 10 μ M or 50 μ M solutions of the triazole-based inhibitors, and the plate was incubated another 72 h at 37 °C. Ninety-six hours post infection, the cell culture media were removed and cells and intracellular amastigotes were lysed in 100 μ L lysis buffer (25 mM Tris pH 8, 2 mM EDTA, 2 mM DTT, 1% Triton X-100, 10% glycerol in ultrapure MQ water) for 10 min at 37 °C. Then, 100 μ L 2× reaction buffer (200 mM sodium phosphate pH 8, 2 mM MgCl2, 100 mM 2-mercaptoethanol and 1.33 mg.mL−1 o-nitrophenyl- β -galactoside (ONPG)) was added and the reaction was allowed to proceed until a yellow color developed (1–2 h at 37 °C). The absorbance at 420 nm was measured in a Synergy HTX multi-mode microplate reader (Biotek Instruments, Winooski, VT, USA) and normalized to the value obtained for the infection in the absence of inhibitors. To determine its IC50, Ald-6 was evaluated at nine different concentrations (0.098, 0.195, 0.39, 0.78, 1.56 3.13 6.25, 12.5 and 50.0 μ M), following the same protocol. Each infection condition and inhibitor were tested in triplicates in four independent experiments.

3.3.4. Cytotoxicity Assay

To evaluate the viability of Vero cells, the Alamar blue method was followed [96]. First, 104 Vero cells/well in 100 μ L MEM-10% FBS were seeded in 96-well culture plates. After 24 h, 50 μ M inhibitor solution was added to cell monolayers and incubated for 96 h. After incubation, resazurin solution (final concentration 10 μ g·mL−1) was added to each well, and fluorescence was measured after 2 h in a Synergy HTX multi-mode microplate reader (Biotek Instruments, Winooski, VT, USA) using the 530–560 nm emission and 590 nm excitation filters. Each condition was tested in triplicate and in at least three independent experiments.

3.3.5. Statistical Analysis

All statistical analyses were performed with GraphPad Prism 6 and complemented with Python-based scripts.

4. Conclusions

The present work has been conducted in the context of CD, which remains a neglected disease requiring continuous scientific efforts to develop potent and safe drugs. The validation of CZP as a therapeutic target for T. cruzi and the renewed interest in TCI as pharmacotherapeutic alternatives stand out. Building on the previous research, we have further demonstrated the potential of the triazole scaffold to develop peptidomimetic antichagasic agents targeting CZP. The versatility of triazole chemistry has enabled us to explore a vast chemical space through in silico methodologies, with good predictions of CZP inhibition, which were largely validated through complementary experimental efforts.
In this study, we elucidated a comprehensive structure–activity relationship (SAR) concerning the binding mode of potent and selective CZP inhibitors. Our results confirm that the chemical space for effective bioactivity translation from CZP to T. cruzi is indeed narrow but compatible with adequate selectivity over hCatL.
Carefully designed screening campaigns, coupled with robust criteria for sampling the compound regions to be experimentally explored, enhance the chances of identifying candidates with the desired trypanocidal effect. Following this approach, we identified a potent, safe and selective antichagasic candidate (Ald-6) that merits further advancement to more complex preclinical bioactivity assessments, such as in vivo models. Additionally, we delineated a confined chemical space where the likelihood of identifying clinically useful CZP TCI is maximized.
Overall, in agreement with the recent findings, we emphasize the critical need to include cell-based T. cruzi infection assays alongside classical CZP inhibition evaluations in the preclinical development pipeline of antichagasic agents.
As ongoing work related to the presented research efforts, additional WH groups are envisioned that not only may provide an adequate reactivity for Cys25 labeling but may also contribute to the specific recognition of the inhibitor within the CZP catalytic site. In this respect, the WH groups reported in the present work (i.e., Es and Ald) are feasible to be further modified towards reactive groups that provide additional intermolecular recognition contributions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29174224/s1, Additional experimental details and findings represented or summarized in figures and tables. Synthesis details, characterization, and 1H NMR and 13C NMR spectra for all precursors and triazole derivatives described and evaluated. Reference [97] is cited in the Supplementary Material.

Author Contributions

Conceptualization, J.P.C., W.D. and M.A.Q.; methodology, J.P.C., G.D.A., R.S.F., W.D. and M.A.Q.; software, J.P.C. and M.A.Q.; validation, J.P.C., L.A.D., V.C.S. and S.C.V.L.; formal analysis, J.P.C., L.A.D., V.C.S. and S.C.V.L.; investigation, J.P.C., L.A.D., V.C.S. and S.C.V.L.; resources, G.D.A., R.S.F., W.D. and M.A.Q.; data curation, J.P.C.; writing—original draft preparation, J.P.C.; writing—review and editing, L.A.D., V.C.S., S.C.V.L., G.D.A., R.S.F., W.D. and M.A.Q.; visualization, J.P.C.; supervision, G.D.A., R.S.F., W.D. and M.A.Q.; project administration, W.D. and M.A.Q.; funding acquisition, W.D. and M.A.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by grants from Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), PIP-CONICET 2021-2023 Nbr. 11220200100776CO; Secretaría de Ciencia y Técnica (SeCyT-UNC 2023), Prj. Nbr: 336-202301-00159CB; Agencia Nacional de Promoción Científica y Tecnología, FONCYT PICT-2021-I-A-00487; the Research Foundation Flanders (FWO), I002720N and I001920N; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) 310197/2021-0 and Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG) grants BPD-00076-22, APQ-00789-22 and RED-00096-22.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The SMILES of the 20 4FPMK derivatives [46] of the training set were generated using MarvinSketch (https://docs.chemaxon.com/display/lts-europium/marvinsketch.md). The SMILES lists of building blocks were obtained from the eMolecules database (https://www.emolecules.com/). The RDKit (https://www.rdkit.org/) package was employed to construct and filter the libraries of evaluated triazole derivatives, and Dimorphite-DL 1.3 (https://pypi.org/project/dimorphite-dl/) to set their corresponding ionization state. The meeko package (https://pypi.org/project/meeko/) was used to obtain PDBQT files as required for docking studies by AutoDock-GPU (https://github.com/ccsb-scripps/AutoDock-GPU), which was employed at each instance of molecular docking. The protein structure of CZP was retrieved from the Protein Data Bank (PDB ID: 3IUT, https://www.rcsb.org/structure/3IUT, accessed on 28 August 2021). All protein–ligand complex topologies were prepared using the AmberTools22 package as freeware (https://ambermd.org/AmberTools.php), while molecular dynamics simulations were computed using Amber22 (https://ambermd.org/AmberMD.php) through the non-commercial program. Binding free energies were computed using the MMPBSA.py script (https://ambermd.org/). For the visualization and analyses of in silico results, VMD v1.9 software (https://www.ks.uiuc.edu/Research/vmd/), ProLIF (https://github.com/chemosim-lab/ProLIF) and LigPlot (https://www.ebi.ac.uk/thornton-srv/software/LigPlus/) were employed. The spectroscopic data were analyzed by TopSpin 3.6.5 (https://www.bruker.com/en/products-and-solutions/mr/nmr-software/topspin.html), and GraphPad Prism 6 was used for statistical analyses related to biological assays (https://www.graphpad.com/scientific-software/prism/). Most input data and parameters are available for download at https://github.com/alfredoq/Triazole_based_TCI_2024_paper (accessed on 1 March 2024). The full MM-MD trajectories, additional scripts and data are available from the authors upon request.

Acknowledgments

We thank Allison Doak and Brian Shoichet (UCSF) for gently providing recombinant cruzain; and Bart Van Huffel, Gert Steurs and Wim De Borggraeve for technical support with NMR. The authors kindly acknowledge the CCAD-UNC (SNCAD-MinCyT, Argentina) for granting access to supercomputing infrastructure.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
Aldaldehyde
BDX1,3-benzodioxole
CDChagas disease
CPCysteine proteases
CZPCruzipain
Esester
hCatLhuman cathepsin L
TCITargeted covalent inhibitors
THQ1,2,3,4-tetrahydroisoquinoline
vHTSvirtual high throughput screening
WHwarhead
4FPMKtetrafluorophenoxymethyl ketone

References

  1. World Health Organization. Chagas Disease. Available online: https://www.who.int/health-topics/chagas-disease#tab=tab_1 (accessed on 8 August 2023).
  2. World Health Organization. Accelerating Work to Overcome the Global Impact of Neglected Tropical Diseases: A Roadmap for Implementation; Technical report; World Health Organization: Geneva, Switzerland, 2012. [Google Scholar]
  3. World Health Organization. Ending the Neglect to Attain the Sustainable Development Goals: A Sustainability Framework for Action Against Neglected Tropical Diseases 2021–2030; Technical report; World Health Organization: Geneva, Switzerland, 2021. [Google Scholar]
  4. Hochberg, N.S.; Montgomery, S.P. Chagas Disease. Ann. Intern. Med. 2023, 176, ITC17–ITC32. [Google Scholar] [CrossRef]
  5. Morillo, C.; Marin-Neto, J.; Avezum, A.; Sosa-Estani, S.; Rassi, A., Jr.; Rosas, F.; Villena, E.; Quiroz, R.; Bonilla, R.; Britto, C.; et al. Randomized trial of benznidazole for chronic Chagas’ cardiomyopathy. N. Engl. J. Med. 2015, 373, 1295–1306. [Google Scholar] [CrossRef] [PubMed]
  6. Pérez-Molina, J.A.; Molina, I. Chagas disease. Lancet 2018, 391, 82–94. [Google Scholar] [CrossRef]
  7. Pathak, S.; Bhardwaj, M.; Agrawal, N.; Bhardwaj, A. A comprehensive review on potential candidates for the treatment of chagas disease. Chem. Biol. Drug Des. 2023, 102, 587–605. [Google Scholar] [CrossRef] [PubMed]
  8. Cazzulo, J.J. Proteinases of Trypanosoma cruzi: Potential targets for the chemotherapy of Chagas disease. Curr. Top. Med. Chem. 2002, 2, 1261–1271. [Google Scholar] [CrossRef] [PubMed]
  9. Cazzulo, J.; Stoka, V.; Turk, V. Cruzipain, the major cysteine proteinase from the protozoan parasite Trypanosoma cruzi. Biol. Chem. Hoppe Seyler 1997, 378, 1–10. [Google Scholar]
  10. Del Nery, E.; Juliano, M.; Meldal, M.; Svendsen, I.; Scharfstein, J.; Walmsley, A.; Juliano, L. Characterization of the substrate specificity of the major cysteine protease (cruzipain) from Trypanosoma cruzi using a portion-mixing combinatorial library and fluorogenic peptides. Biochem. J. 1997, 323, 427–433. [Google Scholar] [CrossRef]
  11. Hayes, D.M.; Kollman, P.A. A Comparison of the Energetics of Proton Transfer in the Serine and Cysteine “Charge Relay” Systems and the Role of the Protein Electrostatic Potential on the Proton Transfer Energetics. In Proceedings of the Catalysis in Chemistry and Biochemistry Theory and Experiment: Proceedings of the Twelfth Jerusalem Symposium on Quantum Chemistry and Biochemistry, Jerusalem, Israel, 2–4 April 1979; Springer: Berlin/Heidelberg, Germany, 1979; pp. 77–90. [Google Scholar]
  12. Duschak, V.G.; Couto, A.S. Cruzipain, the major cysteine protease of Trypanosoma cruzi: A sulfated glycoprotein antigen as relevant candidate for vaccine development and drug target. A review. Curr. Med. Chem. 2009, 16, 3174–3202. [Google Scholar] [CrossRef]
  13. Cazzulo, J.; Couso, R.; Raimondi, A.; Wernstedt, C.; Hellman, U. Further characterization and partial amino acid sequence of a cysteine proteinase from Trypanosoma cruzi. Mol. Biochem. Parasitol. 1989, 33, 33–41. [Google Scholar] [CrossRef]
  14. Judice, W.; Cezari, M.; Lima, A.; Scharfstein, J.; Chagas, J.; Tersariol, I.; Juliano, M.; Juliano, L. Comparison of the specificity, stability and individual rate constants with respective activation parameters for the peptidase activity of cruzipain and its recombinant form, cruzain, from Trypanosoma cruzi. Eur. J. Biochem. 2001, 268, 6578–6586. [Google Scholar] [CrossRef]
  15. Robinson, M.W.; Dalton, J.P. Cysteine Proteases of Pathogenic Organisms; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011; Volume 712, ISBN 978-1489979421. [Google Scholar]
  16. Choe, Y.; Brinen, L.; Price, M.; Engel, J.; Lange, M.; Grisostomi, C.; Weston, S.; Pallai, P.; Cheng, H.; Hardy, L.; et al. Development of α-keto-based inhibitors of cruzain, a cysteine protease implicated in Chagas disease. Bioorganic Med. Chem. 2005, 13, 2141–2156. [Google Scholar] [CrossRef] [PubMed]
  17. Beaulieu, C.; Isabel, E.; Fortier, A.; Massé, F.; Mellon, C.; Méthot, N.; Ndao, M.; Nicoll-Griffith, D.; Lee, D.; Park, H.; et al. Identification of potent and reversible cruzipain inhibitors for the treatment of Chagas disease. Bioorganic Med. Chem. Lett. 2010, 20, 7444–7449. [Google Scholar] [CrossRef]
  18. Nicoll-Griffith, D.A. Use of cysteine-reactive small molecules in drug discovery for trypanosomal disease. Expert Opin. Drug Discov. 2012, 7, 353–366. [Google Scholar] [CrossRef] [PubMed]
  19. Bryant, C.; Kerr, I.; Debnath, M.; Ang, K.; Ratnam, J.; Ferreira, R.; Jaishankar, P.; Zhao, D.; Arkin, M.; McKerrow, J.; et al. Novel non-peptidic vinylsulfones targeting the S2 and S3 subsites of parasite cysteine proteases. Bioorganic Med. Chem. Lett. 2009, 19, 6218–6221. [Google Scholar] [CrossRef]
  20. Santos, V.C.; Ferreira, R.S. Computational approaches towards the discovery and optimisation of cruzain inhibitors. Memórias Inst. Oswaldo Cruz 2022, 117, e210385. [Google Scholar] [CrossRef]
  21. Santos, V.; Oliveira, A.; Campos, A.; Reis-Cunha, J.; Bartholomeu, D.; Teixeira, S.; Lima, A.; Ferreira, R. The gene repertoire of the main cysteine protease of Trypanosoma cruzi, cruzipain, reveals four sub-types with distinct active sites. Sci. Rep. 2021, 11, 18231. [Google Scholar] [CrossRef] [PubMed]
  22. Research Collaboratory for Structural Bioinformatics (RCSB). RCSB PDB-3IUT: Crystal Structure of Cruzipain in Complex with a Vinyl Sulfone Inhibitor. Available online: https://www.rcsb.org/structure/3IUT (accessed on 28 August 2021).
  23. Pauli, I.; Rezende, C., Jr.; Slafer, B.; Dessoy, M.; de Souza, M.; Ferreira, L.; Adjanohun, A.; Ferreira, R.; Magalhães, L.; Krogh, R.; et al. Multiparameter Optimization of Trypanocidal Cruzain Inhibitors With In Vivo Activity and Favorable Pharmacokinetics. Front. Pharmacol. 2022, 12, 774069. [Google Scholar] [CrossRef]
  24. Ferreira, R.; Dessoy, M.; Pauli, I.; Souza, M.; Krogh, R.; Sales, A.; Oliva, G.; Dias, L.; Andricopulo, A. Synthesis, biological evaluation, and structure–activity relationships of potent noncovalent and nonpeptidic cruzain inhibitors as anti-Trypanosoma cruzi agents. J. Med. Chem. 2014, 57, 2380–2392. [Google Scholar] [CrossRef]
  25. Ferreira, R.; Pauli, I.; Sampaio, T.; de Souza, M.; Ferreira, L.; Magalhaes, L.; Rezende, C., Jr.; Ferreira, R.; Krogh, R.; Dias, L.; et al. Structure-based and molecular modeling studies for the discovery of cyclic imides as reversible cruzain inhibitors with potent anti-Trypanosoma cruzi activity. Front. Chem. 2019, 7, 798. [Google Scholar] [CrossRef] [PubMed]
  26. Reis, C.; Souza, H.; Leme, R.; Castelo-Branco, F.; Fernandes, T.; Boechat, N.; Dias, L.; Hoelz, L. Study of the dynamic behavior of the cruzain enzyme in free and complexed forms with competitive and noncovalent benzimidazole inhibitors. J. Biomol. Struct. Dyn. 2023, 41, 4368–4382. [Google Scholar] [CrossRef]
  27. Li, L.; Chenna, B.; Yang, K.; Cole, T.; Goodall, Z.; Giardini, M.; Moghadamchargari, Z.; Hernandez, E.; Gomez, J.; Calvet, C.; et al. Self-masked aldehyde inhibitors: A novel strategy for inhibiting cysteine proteases. J. Med. Chem. 2021, 64, 11267–11287. [Google Scholar] [CrossRef] [PubMed]
  28. Ndao, M.; Beaulieu, C.; Black, W.; Isabel, E.; Vasquez-Camargo, F.; Nath-Chowdhury, M.; Massé, F.; Mellon, C.; Methot, N.; Nicoll-Griffith, D. Reversible cysteine protease inhibitors show promise for a Chagas disease cure. Antimicrob. Agents Chemother. 2014, 58, 1167–1178. [Google Scholar] [CrossRef]
  29. Fonseca, N.; da Cruz, L.; da Silva Villela, F.; do Nascimento Pereira, G.; de Siqueira-Neto, J.; Kellar, D.; Suzuki, B.; Ray, D.; de Souza, T.; Alves, R.; et al. Synthesis of a sugar-based thiosemicarbazone series and structure-activity relationship versus the parasite cysteine proteases rhodesain, cruzain, and Schistosoma mansoni cathepsin B1. Antimicrob. Agents Chemother. 2015, 59, 2666–2677. [Google Scholar] [CrossRef]
  30. Du, X.; Guo, C.; Hansell, E.; Doyle, P.; Caffrey, C.; Holler, T.; McKerrow, J.; Cohen, F. Synthesis and structure-activity relationship study of potent trypanocidal thiosemicarbazone inhibitors of the trypanosomal cysteine protease cruzain. J. Med. Chem. 2002, 45, 2695–2707. [Google Scholar] [CrossRef] [PubMed]
  31. Chipeleme, A.; Gut, J.; Rosenthal, P.; Chibale, K. Synthesis and biological evaluation of phenolic Mannich bases of benzaldehyde and (thio) semicarbazone derivatives against the cysteine protease falcipain-2 and a chloroquine resistant strain of Plasmodium falciparum. Bioorganic Med. Chem. 2007, 15, 273–282. [Google Scholar] [CrossRef] [PubMed]
  32. Engel, J.; Doyle, P.; Hsieh, I.; McKerrow, J. Cysteine protease inhibitors cure an experimental Trypanosoma cruzi infection. J. Exp. Med. 1998, 188, 725–734. [Google Scholar] [CrossRef]
  33. McKerrow, J.; Doyle, P.; Engel, J.; Podust, L.; Robertson, S.; Ferreira, R.; Saxton, T.; Arkin, M.; Kerr, I.; Brinen, L.; et al. Two approaches to discovering and developing new drugs for Chagas disease. Memórias Inst. Oswaldo Cruz 2009, 104, 263–269. [Google Scholar] [CrossRef]
  34. Kerr, I.; Lee, J.; Farady, C.; Marion, R.; Rickert, M.; Sajid, M.; Pandey, K.; Caffrey, C.; Legac, J.; Hansell, E.; et al. Vinyl sulfones as antiparasitic agents and a structural basis for drug design. J. Biol. Chem. 2009, 284, 25697–25703. [Google Scholar] [CrossRef] [PubMed]
  35. Silva, J.; Cianni, L.; Araujo, D.; Batista, P.; de Vita, D.; Rosini, F.; Leitão, A.; Lameira, J.; Montanari, C. Assessment of the cruzain cysteine protease reversible and irreversible covalent inhibition mechanism. J. Chem. Inf. Model. 2020, 60, 1666–1677. [Google Scholar] [CrossRef]
  36. Lonsdale, R.; Ward, R.A. Structure-based design of targeted covalent inhibitors. Chem. Soc. Rev. 2018, 47, 3816–3830. [Google Scholar] [CrossRef]
  37. Singh, J. The ascension of targeted covalent inhibitors. J. Med. Chem. 2022, 65, 5886–5901. [Google Scholar] [CrossRef] [PubMed]
  38. Baillie, T.A. Targeted covalent inhibitors for drug design. Angew. Chem. Int. Ed. 2016, 55, 13408–13421. [Google Scholar] [CrossRef]
  39. Gehringer, M.; Laufer, S.A. Emerging and re-emerging warheads for targeted covalent inhibitors: Applications in medicinal chemistry and chemical biology. J. Med. Chem. 2018, 62, 5673–5724. [Google Scholar] [CrossRef]
  40. Lameiro, R.F.; Montanari, C.A. Investigating the lack of translation from cruzain inhibition to Trypanosoma cruzi activity with machine learning and chemical space analyses. ChemMedChem 2023, 18, e202200434. [Google Scholar] [CrossRef]
  41. Craik, D.; Fairlie, D.; Liras, S.; Price, D. The future of peptide-based drugs. Chem. Biol. Drug Des. 2013, 81, 136–147. [Google Scholar] [CrossRef]
  42. Hamman, J.; Enslin, G.; Kotzé, A. Oral delivery of peptide drugs: Barriers and developments. BioDrugs 2005, 19, 165–177. [Google Scholar] [CrossRef] [PubMed]
  43. Tran, H.; ElSayed, M.E. Progress and limitations of oral peptide delivery as a potentially transformative therapy. Expert Opin. Drug Deliv. 2022, 19, 163–178. [Google Scholar] [CrossRef] [PubMed]
  44. Brak, K.; Doyle, P.; McKerrow, J.; Ellman, J. Identification of a new class of nonpeptidic inhibitors of cruzain. J. Am. Chem. Soc. 2008, 130, 6404–6410. [Google Scholar] [CrossRef] [PubMed]
  45. Brak, K.; Kerr, I.; Barrett, K.; Fuchi, N.; Debnath, M.; Ang, K.; Engel, J.; McKerrow, J.; Doyle, P.; Brinen, L.; et al. Nonpeptidic tetrafluorophenoxymethyl ketone cruzain inhibitors as promising new leads for Chagas disease chemotherapy. J. Med. Chem. 2010, 53, 1763–1773. [Google Scholar] [CrossRef]
  46. Neitz, R.; Bryant, C.; Chen, S.; Gut, J.; Caselli, E.; Ponce, S.; Chowdhury, S.; Xu, H.; Arkin, M.; Ellman, J.; et al. Tetrafluorophenoxymethyl ketone cruzain inhibitors with improved pharmacokinetic properties as therapeutic leads for Chagas’ disease. Bioorganic Med. Chem. Lett. 2015, 25, 4834–4837. [Google Scholar] [CrossRef]
  47. Huisgen, R. 1,3-Dipolare Cycloadditionen Rückschau und Ausblick. Angew. Chem. 1963, 75, 604–637. [Google Scholar] [CrossRef]
  48. Huisgen, R. Kinetik und Mechanismus 1,3-Dipolarer Cycloadditionen. Angew. Chem. 1963, 75, 742–754. [Google Scholar] [CrossRef]
  49. Agalave, S.; Maujan, S.; Pore, V. Click chemistry: 1,2,3-triazoles as pharmacophores. Chem. Asian J. 2011, 6, 2696–2718. [Google Scholar] [CrossRef] [PubMed]
  50. Rostovtsev, V.; Green, L.; Fokin, V.; Sharpless, K. A stepwise huisgen cycloaddition process: Copper (I)-catalyzed regioselective “ligation” of azides and terminal alkynes. Angew. Chem. 2002, 114, 2708–2711. [Google Scholar] [CrossRef]
  51. Tornøe, C.; Christensen, C.; Meldal, M. Peptidotriazoles on solid phase:[1,2,3]-triazoles by regiospecific copper (I)-catalyzed 1,3-dipolar cycloadditions of terminal alkynes to azides. J. Org. Chem. 2002, 67, 3057–3064. [Google Scholar] [CrossRef]
  52. Potter, G.; Jayson, G.; Miller, G.; Gardiner, J. An updated synthesis of the diazo-transfer reagent imidazole-1-sulfonyl azide hydrogen sulfate. J. Org. Chem. 2016, 81, 3443–3446. [Google Scholar] [CrossRef]
  53. Goddard-Borger, E.D.; Stick, R.V. An efficient, inexpensive, and shelf-stable diazotransfer reagent: Imidazole-1-sulfonyl azide hydrochloride. Org. Lett. 2007, 9, 3797–3800. [Google Scholar] [CrossRef]
  54. v. Wachenfeldt, H.; Paulsen, F.; Sundin, A.; Strand, D. Synthesis of Substituted Oxazoles from N-Benzyl Propargyl Amines and Acid Chlorides. Eur. J. Org. Chem. 2013, 2013, 4578–4585. [Google Scholar] [CrossRef]
  55. Ermolat’ev, D.; Feng, H.; Song, G.; Van der Eycken, E. Copper (I)-Catalyzed Decarboxylative Coupling of Propiolic Acids with Secondary Amines and Aldehydes. Eur. J. Org. Chem. 2014, 2014, 5346–5350. [Google Scholar] [CrossRef]
  56. Li, J.; Li, L.; Vessally, E. Decarboxylative A3-coupling reactions: An overview. J. Chin. Chem. Soc. 2021, 68, 13–26. [Google Scholar] [CrossRef]
  57. Wood, W.; Patterson, A.; Tsuruoka, H.; Jain, R.; Ellman, J. Substrate Activity Screening: A Fragment-Based Method for the Rapid Identification of Nonpeptidic Protease Inhibitors. J. Am. Chem. Soc. 2005, 127, 15521–15527. [Google Scholar] [CrossRef]
  58. eMolecules. Available online: https://www.emolecules.com/ (accessed on 30 August 2023).
  59. Silva, L.R.; Guimaraes, A.S.; do Nascimento, J.; do Santos Nascimento, I.J.; da Silva, E.B.; McKerrow, J.H.; Cardoso, S.H.; da Silva-Junior, E.F. Computer-aided design of 1, 4-naphthoquinone-based inhibitors targeting cruzain and rhodesain cysteine proteases. Bioorganic Med. Chem. 2021, 41, 116213. [Google Scholar] [CrossRef] [PubMed]
  60. Shoichet, B.K. Screening in a spirit haunted world. Drug Discov. Today 2006, 11, 607–615. [Google Scholar] [CrossRef] [PubMed]
  61. Shoichet, B.K. Interpreting steep dose-response curves in early inhibitor discovery. J. Med. Chem. 2006, 49, 7274–7277. [Google Scholar] [CrossRef]
  62. Jadhav, A.; Ferreira, R.; Klumpp, C.; Mott, B.; Austin, C.; Inglese, J.; Thomas, C.; Maloney, D.; Shoichet, B.; Simeonov, A. Quantitative Analyses of Aggregation, Autofluorescence, and Reactivity Artifacts in a Screen for Inhibitors of a Thiol Protease. J. Med. Chem. 2010, 53, 37–51. [Google Scholar] [CrossRef]
  63. Martins, L.; de Oliveira, R.; Lameira, J.; Ferreira, R. Experimental and Computational Study of Aryl-thiosemicarbazones Inhibiting Cruzain Reveals Reversible Inhibition and a Stepwise Mechanism. J. Chem. Inf. Model. 2023, 63, 1506–1520. [Google Scholar] [CrossRef]
  64. Xu, Y.; Chigan, J.; Li, J.; Ding, H.; Sun, L.; Liu, L.; Hu, Z.; Yang, K. Hydroxamate and thiosemicarbazone: Two highly promising scaffolds for the development of SARS-CoV-2 antivirals. Bioorganic Chem. 2022, 124, 105799. [Google Scholar] [CrossRef] [PubMed]
  65. Copeland, R.; Basavapathruni, A.; Moyer, M.; Scott, M. Impact of enzyme concentration and residence time on apparent activity recovery in jump dilution analysis. Anal. Biochem. 2011, 416, 206–210. [Google Scholar] [CrossRef]
  66. Morrison, J.F.; Walsh, C.T. The behavior and significance of slow-binding enzyme inhibitors. Adv. Enzymol. Relat. Areas Mol. Biol. 1988, 61, 201–301. [Google Scholar]
  67. Zhu, L.; George, S.; Schmidt, M.; Al-Gharabli, S.; Rademann, J.; Hilgenfeld, R. Peptide aldehyde inhibitors challenge the substrate specificity of the SARS-coronavirus main protease. Antivir. Res. 2011, 92, 204–212. [Google Scholar] [CrossRef] [PubMed]
  68. Dax, C.; Coinçon, M.; Sygusch, J.; Blonski, C. Hydroxynaphthaldehyde phosphate derivatives as potent covalent Schiff base inhibitors of fructose-1,6-bisphosphate aldolase. Biochemistry 2005, 44, 5430–5443. [Google Scholar] [CrossRef] [PubMed]
  69. Philip, V.; Harris, J.; Adams, R.; Nguyen, D.; Spiers, J.; Baudry, J.; Howell, E.; Hinde, R. A survey of aspartate- phenylalanine and glutamate- phenylalanine interactions in the protein data bank: Searching for anion- π pairs. Biochemistry 2011, 50, 2939–2950. [Google Scholar] [CrossRef] [PubMed]
  70. Cihan Sorkun, M.; Mullaj, D.; Koelman, J.; Er, S. ChemPlot, a python library for chemical space visualization. Chemistry-Methods 2022, 2, e202200005. [Google Scholar] [CrossRef]
  71. Medina-Franco, J.; Martínez-Mayorga, K.; Giulianotti, M.; Houghten, R.; Pinilla, C. Visualization of the chemical space in drug discovery. Curr. Comput.-Aided Drug Des. 2008, 4, 322–333. [Google Scholar] [CrossRef]
  72. Reymond, J.; Van Deursen, R.; Blum, L.; Ruddigkeit, L. Chemical space as a source for new drugs. MedChemComm 2010, 1, 30–38. [Google Scholar] [CrossRef]
  73. Lipinski, C.; Hopkins, A. Navigating chemical space for biology and medicine. Nature 2004, 432, 855–861. [Google Scholar] [CrossRef]
  74. Oprea, T.I. Chemical space navigation in lead discovery. Curr. Opin. Chem. Biol. 2002, 6, 384–389. [Google Scholar] [CrossRef]
  75. Csizmadia, P. MarvinSketch and MarvinView: Molecule Applets for the World Wide Web. In Proceedings of the 3rd International Electronic Conference on Synthetic Organic Chemistry, Basel, Switzerland, 1–30 November 1999. [Google Scholar]
  76. Landrum, G. RDKit: A software suite for cheminformatics, computational chemistry, and predictive modeling. Greg Landrum 2013, 8, 31. [Google Scholar]
  77. Ropp, P.; Kaminsky, J.; Yablonski, S.; Durrant, J. Dimorphite-DL: An open-source program for enumerating the ionization states of drug-like small molecules. J. Cheminform. 2019, 11, 1–8. [Google Scholar] [CrossRef]
  78. Halgren, T.A. Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J. Comput. Chem. 1996, 17, 490–519. [Google Scholar] [CrossRef]
  79. Morris, G.; Huey, R.; Lindstrom, W.; Sanner, M.; Belew, R.; Goodsell, D.; Olson, A. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef]
  80. Santos-Martins, D.; Solis-Vasquez, L.; Tillack, A.; Sanner, M.; Koch, A.; Forli, S. Accelerating AutoDock4 with GPUs and Gradient-Based Local Search. J. Chem. Theory Comput. 2021, 17, 1060–1073. [Google Scholar] [CrossRef]
  81. Case, D.; Duke, R.; Walker, R.; Skrynnikov, N.; Cheatham, T., III; Mikhailovskii, O.; Simmerling, C.; Xue, Y.; Roitberg, A.; Izmailov, S.; et al. Amber 2022; University of California: San Francisco, CA, USA, 2022. [Google Scholar]
  82. Kollman, P.; Massova, I.; Reyes, C.; Kuhn, B.; Huo, S.; Chong, L.; Lee, M.; Lee, T.; Duan, Y.; Wang, W.; et al. Calculating structures and free energies of complex molecules: Combining molecular mechanics and continuum models. Accounts Chem. Res. 2000, 33, 889–897. [Google Scholar] [CrossRef] [PubMed]
  83. Kuhn, B.; Gerber, P.; Schulz-Gasch, T.; Stahl, M. Validation and use of the MM-PBSA approach for drug discovery. J. Med. Chem. 2005, 48, 4040–4048. [Google Scholar] [CrossRef] [PubMed]
  84. Miller III, B.; McGee, T., Jr.; Swails, J.; Homeyer, N.; Gohlke, H.; Roitberg, A. MMPBSA.py: An efficient program for end-state free energy calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [CrossRef]
  85. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef] [PubMed]
  86. Cédric, B.; Sébastien, F. ProLIF: A library to encode molecular interactions as fingerprints. J. Cheminformatics 2021, 13, 72. [Google Scholar]
  87. Wallace, A.C.; Laskowski, R.A.; Thornton, J.M. LIGPLOT: A program to generate schematic diagrams of protein-ligand interactions. Protein Eng. Des. Sel. 1995, 8, 127–134. [Google Scholar] [CrossRef]
  88. Centro de Computación de Alto Desempeño de la Universidad Nacional de Córdoba. Available online: https://ccad.unc.edu.ar/ (accessed on 1 March 2024).
  89. Fulmer, G.; Miller, A.; Sherden, N.; Gottlieb, H.; Nudelman, A.; Stoltz, B.; Bercaw, J.; Goldberg, K. NMR chemical shifts of trace impurities: Common laboratory solvents, organics, and gases in deuterated solvents relevant to the organometallic chemist. Organometallics 2010, 29, 2176–2179. [Google Scholar] [CrossRef]
  90. Zsabka, P.; Opsomer, T.; Van Hecke, K.; Dehaen, W.; Wilden, A.; Modolo, G.; Verwerft, M.; Binnemans, K.; Cardinaels, T. Solvent Extraction Studies for the Separation of Trivalent Actinides from Lanthanides with a Triazole-functionalized 1,10-phenanthroline Extractant. Solvent Extr. Ion Exch. 2020, 38, 719–734. [Google Scholar] [CrossRef]
  91. do Valle Moreira, T.; Martins, L.; Diniz, L.; Bernardes, T.; de Oliveira, R.; Ferreira, R. Screening the Pathogen Box to Discover and Characterize New Cruzain and Tbr CatL Inhibitors. Pathogens 2023, 12, 251. [Google Scholar] [CrossRef] [PubMed]
  92. Barbosa da Silva, E.; Rocha, D.; Fortes, I.; Yang, W.; Monti, L.; Siqueira-Neto, J.; Caffrey, C.; McKerrow, J.; Andrade, S.; Ferreira, R. Structure-based optimization of quinazolines as cruzain and Tbr CATL inhibitors. J. Med. Chem. 2021, 64, 13054–13071. [Google Scholar] [CrossRef] [PubMed]
  93. Pereira, G.; da Silva, E.; Braga, S.; Leite, P.; Martins, L.; Vieira, R.; Soh, W.; Soh, W.; Villela, F.; Costa, F.; et al. Discovery and characterization of trypanocidal cysteine protease inhibitors from the ‘malaria box’. Eur. J. Med. Chem. 2019, 179, 765–778. [Google Scholar] [CrossRef] [PubMed]
  94. Boudreau, P.; Miller, B.; McCall, L.; Almaliti, J.; Reher, R.; Hirata, K.; Le, T.; Siqueira-Neto, J.; Hook, V.; Gerwick, W. Design of gallinamide A analogs as potent inhibitors of the cysteine proteases human cathepsin L and Trypanosoma cruzi cruzain. J. Med. Chem. 2019, 62, 9026–9044. [Google Scholar] [CrossRef]
  95. Buckner, F.; Verlinde, C.; La Flamme, A.; Van Voorhis, W. Efficient technique for screening drugs for activity against Trypanosoma cruzi using parasites expressing beta-galactosidase. Antimicrob. Agents Chemother. 1996, 40, 2592–2597. [Google Scholar] [CrossRef]
  96. Lafon-Hughes, L.; Villamil, S.; Larrea, S. Tankyrase inhibitors hinder Trypanosoma cruzi infection by altering host-cell signalling pathways. Parasitology 2021, 148, 1680–1690. [Google Scholar] [CrossRef]
  97. Chen, Z.; Sun, T.; Qing, G. cAMP-modulated biomimetic ionic nanochannels based on a smart polymer. J. Mater. Chem. B 2019, 7, 3710–3715. [Google Scholar] [CrossRef]
Figure 1. (a) 3D structure of CZP (PDB 3IUT [22]). (b) 3D structure of CZP active site, highlighting the S1’ (red), S1 (green), S2 (blue), S3 (purple) and catalytic (yellow) subsites. The catalytic triad and oxyanion hole are represented in licorice.
Figure 1. (a) 3D structure of CZP (PDB 3IUT [22]). (b) 3D structure of CZP active site, highlighting the S1’ (red), S1 (green), S2 (blue), S3 (purple) and catalytic (yellow) subsites. The catalytic triad and oxyanion hole are represented in licorice.
Molecules 29 04224 g001
Figure 2. Peptidomimetic (Ts-370) and peptide-based (K777) promising CZP inhibitors reported in literature, with poor translation to T. cruzi inhibitory activity [32,45].
Figure 2. Peptidomimetic (Ts-370) and peptide-based (K777) promising CZP inhibitors reported in literature, with poor translation to T. cruzi inhibitory activity [32,45].
Molecules 29 04224 g002
Figure 3. Synthetic core scheme for the synthesis of 1,4-disubstituted 1,2,3-triazole-based potential inhibitors of CZP.
Figure 3. Synthetic core scheme for the synthesis of 1,4-disubstituted 1,2,3-triazole-based potential inhibitors of CZP.
Molecules 29 04224 g003
Figure 4. Summary of the docking-based vHTS results focused on R1 with CZP.
Figure 4. Summary of the docking-based vHTS results focused on R1 with CZP.
Molecules 29 04224 g004
Figure 5. (a) Summary of the docking-based vHTS results for the 1,2,3-triazoles obtained by combinatorial synthesis. Clusters of optimal (green box, upper right) and less favored (fuchsia box, bottom right) synthesized compounds are highlighted. (b) Chemical structure of the most promising triazole-derivatives identified by means of vHTS and synthesized. (c) Docking pose of Ald-6 (licorice), highlighting its interaction with Asp161 (CPK).
Figure 5. (a) Summary of the docking-based vHTS results for the 1,2,3-triazoles obtained by combinatorial synthesis. Clusters of optimal (green box, upper right) and less favored (fuchsia box, bottom right) synthesized compounds are highlighted. (b) Chemical structure of the most promising triazole-derivatives identified by means of vHTS and synthesized. (c) Docking pose of Ald-6 (licorice), highlighting its interaction with Asp161 (CPK).
Molecules 29 04224 g005
Figure 6. Synthetic core scheme for the synthesis of Es- and Ald-derived 1,4-disubstituted 1,2,3-triazoles. * The yields of Ald-15, Ald-16 and Ald-24 are not included in the range since they could not be isolated.
Figure 6. Synthetic core scheme for the synthesis of Es- and Ald-derived 1,4-disubstituted 1,2,3-triazoles. * The yields of Ald-15, Ald-16 and Ald-24 are not included in the range since they could not be isolated.
Molecules 29 04224 g006
Figure 7. Mechanism of CZP inhibition by Es-15, evaluated in triplicate. (a) Michaelis–Menten plot. The curves correspond to the fitting of a mixed-model of inhibition to data. (b) Plot of the K m a p p for varying Es-15 concentrations (c) Lineweaver–Burk reciprocal plot.
Figure 7. Mechanism of CZP inhibition by Es-15, evaluated in triplicate. (a) Michaelis–Menten plot. The curves correspond to the fitting of a mixed-model of inhibition to data. (b) Plot of the K m a p p for varying Es-15 concentrations (c) Lineweaver–Burk reciprocal plot.
Molecules 29 04224 g007
Figure 8. Jump dilution assay results for (a) Ald-6 and (b) Ald-10. The inhibitors were incubated at 10 × IC50 with 100 × [CZP] for 30 min, followed by 100× dilution and 2 h monitoring. The ligand behavior during the first 15 min of the assay is depicted in the magnifications on the top left. The mean and SEM for three replicates are shown.
Figure 8. Jump dilution assay results for (a) Ald-6 and (b) Ald-10. The inhibitors were incubated at 10 × IC50 with 100 × [CZP] for 30 min, followed by 100× dilution and 2 h monitoring. The ligand behavior during the first 15 min of the assay is depicted in the magnifications on the top left. The mean and SEM for three replicates are shown.
Molecules 29 04224 g008
Figure 9. (a) Interaction pattern identified for Ald-6:CZP. The alternative poses of the R1 substituent are represented (gray). (b) MM-MD snapshots of Ald-6 (licorice), highlighting the alternative poses of its R1 when interacting with Ser64 and Asp161 (CPK). (c) Distance plots of R1-n-butylamino of Ald-6 with CZP-Asp161 and Ser64.
Figure 9. (a) Interaction pattern identified for Ald-6:CZP. The alternative poses of the R1 substituent are represented (gray). (b) MM-MD snapshots of Ald-6 (licorice), highlighting the alternative poses of its R1 when interacting with Ser64 and Asp161 (CPK). (c) Distance plots of R1-n-butylamino of Ald-6 with CZP-Asp161 and Ser64.
Molecules 29 04224 g009
Figure 10. Chemical space of the triazole-based derivatives feasible to be synthesized with the biologically evaluated ones colored according to their IC50 against CZP. The subplot on the right shows a magnification of the most populated cluster.
Figure 10. Chemical space of the triazole-based derivatives feasible to be synthesized with the biologically evaluated ones colored according to their IC50 against CZP. The subplot on the right shows a magnification of the most populated cluster.
Molecules 29 04224 g010
Figure 11. (a) Mean ± SD of T. cruzi viability quantified in triplicate ([Inh.] = 50 μ M) 96 h after cell seeding (72 h p.i.) by β -gal activity assay and expressed as % of remaining infection. ** [Ald-4] = 20 μ M. (b) Percentage of uninfected Vero cell viability after 96h incubation with triazole-based CZP inhibitors (50 μ M), quantified by Alamar blue assay.
Figure 11. (a) Mean ± SD of T. cruzi viability quantified in triplicate ([Inh.] = 50 μ M) 96 h after cell seeding (72 h p.i.) by β -gal activity assay and expressed as % of remaining infection. ** [Ald-4] = 20 μ M. (b) Percentage of uninfected Vero cell viability after 96h incubation with triazole-based CZP inhibitors (50 μ M), quantified by Alamar blue assay.
Molecules 29 04224 g011
Table 1. Structural, synthetic and enzymatic data of triazole-based Es and Ald derivatives.
Table 1. Structural, synthetic and enzymatic data of triazole-based Es and Ald derivatives.
R3R1Es aYield c% CZP inh. eIC50 ( μ M) f Ald a Yield d % CZP inh. eIC50 ( μ M) f
BDXPheEs-1b89%15 ± 2n.d.Ald-136%64 ± 541 ± 4
ValEs-285%9 ± 7n.d.Ald-253%55 ± 3n.d.
GluEs-350%8 ± 7n.d.Ald-321%8 ± 2n.d.
HisEs-489%5 ± 4n.d.Ald-442%83 ± 25 ± 2
LeuEs-590%5 ± 2n.d.Ald-532%64 ± 27 ± 2
LysEs-648%9 ± 2n.d.Ald-651%81 ± 23.3 ± 0.3
MetEs-7b84%9 ± 8n.d.Ald-759%42 ± 7n.d.
AlaEs-884%8 ± 7n.d.Ald-838%8 ± 4n.d.
PhgEs-982%6 ± 5n.d.Ald-937%11 ± 3n.d.
SerEs-1087%47 ± 6n.d.Ald-1056%87 ± 13.4 ± 0.5
ThrEs-1165%0 ± 2n.d.Ald-1144%58 ± 4n.d.
TrpEs-1268%15 ± 2n.d.Ald-1258%56 ± 745.4 ± 0.3
TyrEs-1351%8 ± 1n.d.Ald-1343%83 ± 213 ± 5
CysEs-1473%53 ± 4n.d.Ald-14---
ArgEs-1543%86 ± 16.8 ± 0.3Ald-15---
THQPheEs-16b76%13 ± 6n.d.Ald-1655%34 ± 5n.d.
HisEs-1754%6 ± 3n.d.Ald-1756%38 ± 6n.d.
MetEs-18b60%0 ± 6n.d.Ald-1843%16 ± 4n.d.
PhgEs-1951%0 ± 3n.d.Ald-1949%3 ± 10n.d.
SerEs-2077%5 ± 4n.d.Ald-2038%83 ± 37 ± 4
ThrEs-2171%0 ± 7n.d.Ald-2137%8 ± 7n.d.
TrpEs-2258%1 ± 6n.d.Ald-2250%55 ± 9n.d.
TyrEs-2351%1 ± 14n.d.Ald-2344%20 ± 4n.d.
ArgEs-2428%6 ± 1n.d.Ald-24---
a The R2 for all of them is c-Hex. b Es-derivatives obtained from L-AA-OEt, while the remaining were synthesized from L-AA-OMe. Yields of c CuAAC and d DIBAL-mediated reduction steps. e The initial screening against CZP was performed in triplicate at 100 μM (% CZP inhib.) with 10’ pre-incubation of the enzyme in the presence of the compounds. f IC50 are the mean ± SEM of two independent experiments. n.d.: not determined.
Table 2. Activity profiles of 1,4-disubstituted 1,2,3-triazole derivatives against CZP and hCatL at 100 μ M after 10 min of incubation. * E-64: positive control.
Table 2. Activity profiles of 1,4-disubstituted 1,2,3-triazole derivatives against CZP and hCatL at 100 μ M after 10 min of incubation. * E-64: positive control.
Compd ID% CZP inhib.% hCatL inhib.
Es-1453 ± 40 ± 26
Es-1586 ± 126 ± 16
Ald-164 ± 50 ± 21
Ald-483 ± 210 ± 3
Ald-564 ± 224 ± 9
Ald-681 ± 226 ± 5
Ald-1087 ± 17 ± 2
Ald-1256 ± 70 ± 16
Ald-1383 ± 20 ± 22
Ald-2083 ± 30 ± 39
E-64 *98 ± 297 ± 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cerutti, J.P.; Diniz, L.A.; Santos, V.C.; Vilchez Larrea, S.C.; Alonso, G.D.; Ferreira, R.S.; Dehaen, W.; Quevedo, M.A. Structure-Aided Computational Design of Triazole-Based Targeted Covalent Inhibitors of Cruzipain. Molecules 2024, 29, 4224. https://doi.org/10.3390/molecules29174224

AMA Style

Cerutti JP, Diniz LA, Santos VC, Vilchez Larrea SC, Alonso GD, Ferreira RS, Dehaen W, Quevedo MA. Structure-Aided Computational Design of Triazole-Based Targeted Covalent Inhibitors of Cruzipain. Molecules. 2024; 29(17):4224. https://doi.org/10.3390/molecules29174224

Chicago/Turabian Style

Cerutti, Juan Pablo, Lucas Abreu Diniz, Viviane Corrêa Santos, Salomé Catalina Vilchez Larrea, Guillermo Daniel Alonso, Rafaela Salgado Ferreira, Wim Dehaen, and Mario Alfredo Quevedo. 2024. "Structure-Aided Computational Design of Triazole-Based Targeted Covalent Inhibitors of Cruzipain" Molecules 29, no. 17: 4224. https://doi.org/10.3390/molecules29174224

APA Style

Cerutti, J. P., Diniz, L. A., Santos, V. C., Vilchez Larrea, S. C., Alonso, G. D., Ferreira, R. S., Dehaen, W., & Quevedo, M. A. (2024). Structure-Aided Computational Design of Triazole-Based Targeted Covalent Inhibitors of Cruzipain. Molecules, 29(17), 4224. https://doi.org/10.3390/molecules29174224

Article Metrics

Back to TopTop