Trypanosoma cruzi is a parasite that infects about 6-7 million people worldwide, mostly in Latin America, causing Chagas disease. Cruzain, the main cysteine protease of T. cruzi, is a validated target for developing drug candidates for Chagas disease. Thiosemicarbazones are one of the most relevant warheads used in covalent inhibitors targeting cruzain. Despite its relevance, the mechanism of inhibition of cruzain by thiosemicarbazones is unknown. Here, we combined experiments and simulations to unveil the covalent inhibition mechanism of cruzain by a thiosemicarbazone-based inhibitor (compound 1). Additionally, we studied a semicarbazone (compound 2), which is structurally similar to compound 1 but does not inhibit cruzain. Assays confirmed the reversibility of inhibition by compound 1 and suggested a two-step mechanism of inhibition. The Ki was estimated to be 36.3 μM and Ki* to be 11.5 μM, suggesting the pre-covalent complex to be relevant for inhibition. Molecular dynamics simulations of compounds 1 and 2 with cruzain were used to propose putative binding modes for the ligands. One-dimensional (1D) quantum mechanics/molecular mechanics (QM/MM) potential of mean force (PMF) and gas-phase energies showed that the attack of Cys25-S- on the C═S or C═O bond yields a more stable intermediate than the attack on the C═N bond of the thiosemicarbazone/semicarbazone. Two-dimensional (2D) QM/MM PMF revealed a putative reaction mechanism for compound 1, involving the proton transfer to the ligand, followed by the Cys25-S- attack at C═S. The ΔG and energy barrier were estimated to be -1.4 and 11.7 kcal/mol, respectively. Overall, our results shed light on the inhibition mechanism of cruzain by thiosemicarbazones.

The emergence of artemisinin-resistant variants of Plasmodium falciparum necessitates the urgent search for novel antimalarial drugs. In this regard, an in silico study to screen antimalarial drug candidates from a series of benzimidazole-thiosemicarbazone hybrid molecules with interesting antiplasmodial properties and explore their falcipain-2 (FP2) inhibitory potentials has been undertaken herein. FP2 is a key cysteine protease that degrades hemoglobin in Plasmodium falciparum and is an important biomolecular target in the development of antimalarial drugs. Pharmacokinetic properties, ADMET profiles, MM/GBSA-based binding free energies, reaction mechanisms, and associated barrier heights have been investigated. DFT, molecular dynamics simulation, molecular docking, and ONIOM methods were used. From the results obtained, four 4N-substituted derivatives of the hybrid molecule (E)-2-(1-(5-chloro-1H-benzo[d]imidazol-2-yl)ethylidene)hydrazine-1-carbothioamide (1A) denoted 1B, 1C, 1D, and 1E are drug-like and promising inhibitors of FP2, exhibiting remarkably small inhibitory constants (5.94 × 10-14 - 2.59 × 10-04 n M) and favorable binding free energies (-30.32 to -17.17 kcal/mol). Moreover, the ONIOM results have revealed that 1B and possibly 1C and 1D may act as covalent inhibitors of FP2. The rate-determining step of the thermodynamically favorable covalent binding mechanism occurs across a surmountable barrier height of 24.18 kcal/mol in water and 28.42 kcal/mol in diethyl ether. Our findings are useful for further experimental investigations on the antimalarial activities of the hybrid molecules studied.

Virtual screening a collection of ~ 25,000 ChemBridge molecule collection identified two nitrogenous heterocyclic molecules, 12 and 15, with potential dual inhibitory properties against trypanosomal cruzain and rhodesain cysteine proteases. Similarity search in DrugBank found the two virtual hits with novel chemical structures with unreported anti-trypanosomal activities. Investigations into the binding mechanism by molecular dynamics simulations for 100 ns revealed the molecules were able to occupy the binding sites and stabilise the protease complexes. Binding affinities calculated using the MM/PBSA method for the last 20 ns showed that the virtual hits have comparable binding affinities to other known inhibitors from literature suggesting both molecules as promising scaffolds with dual cruzain and rhodesain inhibition properties, i.e. 12 has predicted ΔGbind values of - 38.1 and - 38.2 kcal/mol to cruzain and rhodesain, respectively, and 15 has predicted ΔGbind values of - 34.4 and - 25.8 kcal/mol to rhodesain. Per residue binding free energy decomposition studies and visual inspection at 100 ns snapshots revealed hydrogen bonding and non-polar attractions with important amino acid residues that contributed to the ΔGbind values. The interactions are similar to those previously reported in the literature. The overall ADMET predictions for the two molecules were favourable for drug development with acceptable pharmacokinetic profiles and adequate oral bioavailability.

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

The development of cruzipain inhibitors represents one of the most attractive challenges in the search for drugs for the treatment of Chagas disease. A recombinant form of this enzyme, cruzain, has been crystallized with numerous inhibitors, excluding thiosemicarbazones. These compounds have been established as potent inhibitors of cruzain, although there is very little data in the literature of thiosemicarbazones tested on cruzipain. In this work, we present the results of the evaluation of eleven thiosemicarbazones on cruzipain, isolated from T. cruzi epimastigotes, six of them previously evaluated on cruzain. For these latter, we studied through computational methods, the mode of interaction with the active site of cruzain and the contribution of geometric parameters to the possible mechanism of action involved in the observed inhibition. Finally, from some geometric parameters analyzed on modeled TSC-cruzain complexes, a semi-quantitative relationship was established that could explain the inhibitory activity of thiosemicarbazones on cruzipain, the enzyme actually present in the parasite.

The COVID-19 pandemic and its continuing emerging variants emphasize the need to discover appropriate treatment, where vaccines alone have failed to show complete protection against the new variants of the virus. Therefore, treatment of the infected cases is critical. This paper discusses the bio-guided isolation of three indole diketopiperazine alkaloids, neoechinulin A (1), echinulin (2), and eurocristatine (3), from the Red Sea-derived Aspergillus fumigatus MR2012. Neoechinulin A (1) exhibited a potent inhibitory effect against SARS-CoV-2 Mpro with IC50 value of 0.47 μM, which is comparable to the reference standard GC376. Despite the structural similarity between the three compounds, only 1 showed a promising effect. The mechanism of inhibition is discussed in light of a series of extensive molecular docking, classical and steered molecular dynamics simulation experiments. This paper sheds light on indole diketopiperazine alkaloids as a potential structural motif against SARS-CoV-2 Mpro. Additionally, it highlights the potential of different molecular docking and molecular dynamics simulation approaches in the discrimination between active and inactive structurally related Mpro inhibitors.

Papain-Like Protease (PLpro) is a key protein for SARS-CoV-2 viral replication which is the cause of the emerging COVID-19 pandemic. Targeting PLpro can suppress viral replication and provide treatment options for COVID-19. Due to the dynamic nature of its binding site loop, PLpro multiple conformations were generated through a long-range 1 micro-second molecular dynamics (MD) simulation. Clustering the MD trajectory enabled us to extract representative structures for the conformational space generated. Adding to the MD representative structures, X-ray structures were involved in an ensemble docking approach to screen the FDA approved drugs for a drug repositioning endeavor. Guided by our recent benchmarking study of SARS-CoV-2 PLpro, FRED docking software was selected for such a virtual screening task. The results highlighted potential consensus binders to many of the MD clusters as well as the newly introduced X-ray structure of PLpro complexed with a small molecule. For instance, three drugs Benserazide, Dobutamine and Masoprocol showed a superior consensus enrichment against the PLpro conformations. Further MD simulations for these drugs complexed with PLpro suggested the superior stability and binding of dobutamine and masoprocol inside the binding site compared to Benserazide. Generally, this approach can facilitate identifying drugs for repositioning via targeting multiple conformations of a crucial target for the rapidly emerging COVID-19 pandemic.

Several derivatives of benzoic acid and semisynthetic alkyl gallates were investigated by an in silico approach to evaluate their potential antiviral activity against SARS-CoV-2 main protease. Molecular docking studies were used to predict their binding affinity and interactions with amino acids residues from the active binding site of SARS-CoV-2 main protease, compared to boceprevir. Deep structural insights and quantum chemical reactivity analysis according to Koopmans\' theorem, as a result of density functional theory (DFT) computations, are reported. Additionally, drug-likeness assessment in terms of Lipinski\'s and Weber\'s rules for pharmaceutical candidates, is provided. The outcomes of docking and key molecular descriptors and properties were forward analyzed by the statistical approach of principal component analysis (PCA) to identify the degree of their correlation. The obtained results suggest two promising candidates for future drug development to fight against the coronavirus infection.

The Ananas comosus stem extract is a complex mixture containing various cysteine ​​proteases of the C1A subfamily, such as bromelain and ananain. This mixture used for centuries in Chinese medicine, has several potential therapeutic applications as anti-cancer, anti-inflammatory and ecchymosis degradation agent. In the present work we determined the structures of bromelain and ananain, both in their free forms and in complex with the inhibitors E64 and TLCK. These structures combined with protease-substrate complexes modeling clearly identified the Glu68 as responsible for the high discrimination of bromelain in favor of substrates with positively charged residues at P2, and unveil the reasons for its weak inhibition by cystatins and E64. Our results with purified and fully active bromelain, ananain and papain show a strong reduction of cell proliferation with MDA-MB231 and A2058 cancer cell lines at a concentration of about 1 μM, control experiments clearly emphasizing the need for proteolytic activity. In contrast, while bromelain and ananain had a strong effect on the proliferation of the OCI-LY19 and HL-60 non-adherent cell lines, papain, the archetypal member of the C1A subfamily, had none. This indicates that, in this case, sequence/structure identity beyond the active site of bromelain and ananain is more important than substrate specificity.

Oseltamivir is a first-line antiviral drug, especially in primary hospitals. During the ongoing outbreak of coronavirus disease 2019 (COVID-19), most patients with COVID-19 who are symptomatic have used oseltamivir. Considering its popular and important role as an antiviral drug, it is necessary to evaluate oseltamivir in the treatment of COVID-19.
To evaluate the effect of oseltamivir against COVID-19.
Swiss-model was used to construct the structure of the N-terminal RNA-binding domain (NRBD) of the nucleoprotein (NC), papain-like protease (PLpro), and RNA-directed RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). TM-align program was performed to compare the structure of the viral proteins with the structure of the neuraminidase of influenza A. Molecular docking was used to analyze the theoretical possibility of effective binding of oseltamivir with the active centers of the viral proteins. In vitro study was used to evaluate the antiviral efficiency of oseltamivir against SARS-CoV-2. By clinical case analysis, we statistically evaluated whether the history of oseltamivir use influenced the progression of the disease.
The structures of NRBD, PLpro, and RdRp were built successfully. The results from TM-align suggested that the S protein, NRBD, 3C-like protease (3CLpro), PLPrO, and RdRp were structurally similar to the influenza A neuraminidase, with TM-scores of 0.30077, 0.19254, 0.28766, 0.30666, and 0.34047, respectively. Interestingly, the active center of 3CL pro was found to be similar to the active center from the neuraminidase of influenza A. Through an analysis of molecular docking, we discovered that oseltamivir carboxylic acid was more favorable to bind to the active site of 3CLpro effectively, but its inhibitory effect was not strong compared with the positive group. Finally, we used in vitro study and retrospective case analysis to verify our speculations. We found that oseltamivir is ineffective against SARS-CoV-2 in vitro study and the clinical use of oseltamivir did not improve the patients\' symptoms and signs and did not slow the disease progression.
We consider that oseltamivir isn\'t suitable for the treatment of COVID-19. During the outbreak of novel coronavirus, when oseltamivir is not effective for the patients after they take it, health workers should be highly vigilant about the possibility of COVID-19.