Falcipain-2 inhibition by suramin and suramin analogues

Plasmodium falciparum proteases as new drug targets with special focus on metalloproteases

Malaria poses a significant global health threat particularly due to the prevalence of Plasmodium falciparum infection. With the emergence of parasite resistance to existing drugs including the recently discovered artemisinin, ongoing research seeks novel therapeutic avenues within the malaria parasite. Proteases are promising drug targets due to their essential roles in parasite biology, including hemoglobin digestion, merozoite invasion, and egress. While exploring the genomic landscape of Plasmodium falciparum, it has been revealed that there are 92 predicted proteases, with only approximately 14 of them having been characterized. These proteases are further distributed among 26 families grouped into five clans: aspartic proteases, cysteine proteases, metalloproteases, serine proteases, and threonine proteases. Focus on metalloprotease class shows further role in organelle processing for mitochondria and apicoplasts suggesting the potential of metalloproteases as viable drug targets. Holistic understanding of the parasite intricate life cycle and identification of potential drug targets are essential for developing effective therapeutic strategies against malaria and mitigating its devastating global impact.

Three Decades of Targeting Falcipains to Develop Antiplasmodial Agents: What have we Learned and What can be Done Next?

Malaria is a devastating infectious disease that affects large swathes of human populations across the planet's tropical regions. It is caused by parasites of the genus Plasmodium, with Plasmodium falciparum being responsible for the most lethal form of the disease. During the intraerythrocytic stage in the human hosts, malaria parasites multiply and degrade hemoglobin (Hb) using a battery of proteases, which include two cysteine proteases, falcipains 2 and 3 (FP-2 and FP-3). Due to their role as major hemoglobinases, FP-2 and FP-3 have been targeted in studies aiming to discover new antimalarials and numerous inhibitors with activity against these enzymes, and parasites in culture have been identified. Nonetheless, cross-inhibition of human cysteine cathepsins remains a serious hurdle to overcome for these compounds to be used clinically. In this article, we have reviewed key functional and structural properties of FP-2/3 and described different compound series reported as inhibitors of these proteases during decades of active research in the field. Special attention is also paid to the wide range of computer-aided drug design (CADD) techniques successfully applied to discover new active compounds. Finally, we provide guidelines that, in our understanding, will help advance the rational discovery of new FP-2/3 inhibitors.

Recent applications of computational methods to allosteric drug discovery

Interest in exploiting allosteric sites for the development of new therapeutics has grown considerably over the last two decades. The chief driving force behind the interest in allostery for drug discovery stems from the fact that in comparison to orthosteric sites, allosteric sites are less conserved across a protein family, thereby offering greater opportunity for selectivity and ultimately tolerability. While there is significant overlap between structure-based drug design for orthosteric and allosteric sites, allosteric sites offer additional challenges mostly involving the need to better understand protein flexibility and its relationship to protein function. Here we examine the extent to which structure-based drug design is impacting allosteric drug design by highlighting several targets across a variety of target classes.

Tetracycline Derivatives Inhibit Plasmodial Cysteine Protease Falcipain-2 through Binding to a Distal Allosteric Site

Allosteric inhibitors regulate enzyme activity from remote and usually specific pockets. As they promise an avenue for less toxic and safer drugs, the identification and characterization of allosteric inhibitors has gained great academic and biomedical interest in recent years. Research on falcipain-2 (FP-2), the major papain-like cysteine hemoglobinase of Plasmodium falciparum, might benefit from this strategy to overcome the low selectivity against human cathepsins shown by active site-directed inhibitors. Encouraged by our previous finding that methacycline inhibits FP-2 noncompetitively, here we assessed other five tetracycline derivatives against this target and characterized their inhibition mechanism. As previously shown for methacycline, tetracycline derivatives inhibited FP-2 in a noncompetitive fashion, with K_i values ranging from 121 to 190 μM. A possible binding to the S' side of the FP-2 active site, similar to that described by X-ray crystallography (PDB: 6SSZ) for the noncompetitive inhibitor E-chalcone 48 (EC48), was experimentally discarded by kinetic analysis using a large peptidyl substrate spanning the whole active site. By combining lengthy molecular dynamics (MD) simulations that allowed methacycline to diffuse from solution to different FP-2 surface regions and free energy calculations, we predicted the most likely binding mode of the ligand. Of note, the proposed binding pose explains the low differences in K_i values observed for the tested tetracycline derivatives and the calculated binding free energies match the experimental values. Overall, this study has implications for the design of novel allosteric inhibitors against FP-2 and sets the basis for further optimization of the tetracycline scaffold to produce more potent and selective inhibitors.

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

Targeted Redesign of Suramin Analogs for Novel Antimicrobial Lead Development

The emergence of new viral infections and drug-resistant bacteria urgently necessitates expedient therapeutic development. Repurposing and redesign of existing drugs against different targets are one potential way in which to accelerate this process. Suramin was initially developed as a successful antiparasitic drug but has also shown promising antiviral and antibacterial activities. However, due to its high conformational flexibility and negative charge, suramin is considered quite promiscuous toward positively charged sites within nucleic acid binding proteins. Although some suramin analogs have been developed against specific targets, only limited structure-activity relationship studies were performed, and virtual screening has yet to be used to identify more specific inhibitor(s) based on its scaffold. Using available structures, we investigated suramin's target diversity, confirming that suramin preferentially binds to protein pockets that are both positively charged and enriched in aromatic or leucine residues. Further, suramin's high conformational flexibility allows adaptation to structurally diverse binding surfaces. From this platform, we developed a framework for structure- and docking-guided elaboration of suramin analog scaffolds using virtual screening of suramin and heparin analogs against a panel of diverse therapeutically relevant viral and bacterial protein targets. Use of this new framework to design potentially specific suramin analogs is exemplified using the SARS-CoV-2 RNA-dependent RNA polymerase and nucleocapsid protein, identifying leads that might inhibit a wide range of coronaviruses. The approach presented here establishes a computational framework for designing suramin analogs against different bacterial and viral targets and repurposing existing drugs for more specific inhibitory activity.

Repurposing Drugs to Fight Hepatic Malaria Parasites

Malaria remains one of the most prevalent infectious diseases worldwide, primarily affecting some of the most vulnerable populations around the globe. Despite achievements in the treatment of this devastating disease, there is still an urgent need for the discovery of new drugs that tackle infection by Plasmodium parasites. However, de novo drug development is a costly and time-consuming process. An alternative strategy is to evaluate the anti-plasmodial activity of compounds that are already approved for other purposes, an approach known as drug repurposing. Here, we will review efforts to assess the anti-plasmodial activity of existing drugs, with an emphasis on the obligatory and clinically silent liver stage of infection. We will also review the current knowledge on the classes of compounds that might be therapeutically relevant against Plasmodium in the context of other communicable diseases that are prevalent in regions where malaria is endemic. Repositioning existing compounds may constitute a faster solution to the current gap of prophylactic and therapeutic drugs that act on Plasmodium parasites, overall contributing to the global effort of malaria eradication.

100 Years of Suramin

Suramin is 100 years old and is still being used to treat the first stage of acute human sleeping sickness, caused by Trypanosoma brucei rhodesiense Suramin is a multifunctional molecule with a wide array of potential applications, from parasitic and viral diseases to cancer, snakebite, and autism. Suramin is also an enigmatic molecule: What are its targets? How does it get into cells in the first place? Here, we provide an overview of the many different candidate targets of suramin and discuss its modes of action and routes of cellular uptake. We reason that, once the polypharmacology of suramin is understood at the molecular level, new, more specific, and less toxic molecules can be identified for the numerous potential applications of suramin.

Identification of Plasmodium dipeptidyl aminopeptidase allosteric inhibitors by high throughput screening

Dipeptidyl aminopeptidases (DPAPs) are cysteine proteases that cleave dipeptides from the N-terminus of protein substrates and have been shown to play important roles in many pathologies including parasitic diseases such as malaria, toxoplasmosis and Chagas's disease. Inhibitors of the mammalian homologue cathepsin C have been used in clinical trials as potential drugs to treat chronic inflammatory disorders, thus proving that these enzymes are druggable. In Plasmodium species, DPAPs play important functions at different stages of parasite development, thus making them potential antimalarial targets. Most DPAP inhibitors developed to date are peptide-based or peptidomimetic competitive inhibitors. Here, we used a high throughput screening approach to identify novel inhibitor scaffolds that block the activity of Plasmodium falciparum DPAP1. Most of the hits identified in this screen also inhibit Plasmodium falciparum DPAP3, cathepsin C, and to a lesser extent other malarial clan CA proteases, indicating that these might be general DPAP inhibitors. Interestingly, our mechanism of inhibition studies indicate that most hits are allosteric inhibitors, which opens a completely new strategy to inhibit these enzymes, study their biological function, and potentially develop new inhibitors as starting points for drug development.

South African Abietane Diterpenoids and Their Analogs as Potential Antimalarials: Novel Insights from Hybrid Computational Approaches

The hemoglobin degradation process in Plasmodium parasites is vital for nutrient acquisition required for their growth and proliferation. In P. falciparum, falcipains (FP-2 and FP-3) are the major hemoglobinases, and remain attractive antimalarial drug targets. Other Plasmodium species also possess highly homologous proteins to FP-2 and FP-3. Although several inhibitors have been designed against these proteins, none has been commercialized due to associated toxicity on human cathepsins (Cat-K, Cat-L and Cat-S). Despite the two enzyme groups sharing a common structural fold and catalytic mechanism, distinct active site variations have been identified, and can be exploited for drug development. Here, we utilize in silico approaches to screen 628 compounds from the South African natural sources to identify potential hits that can selectively inhibit the plasmodial proteases. Using docking studies, seven abietane diterpenoids, binding strongly to the plasmodial proteases, and three additional analogs from PubChem were identified. Important residues involved in ligand stabilization were identified for all potential hits through binding pose analysis and their energetic contribution determined by binding free energy calculations. The identified compounds present important scaffolds that could be further developed as plasmodial protease inhibitors. Previous laboratory assays showed the effect of the seven diterpenoids as antimalarials. Here, for the first time, we demonstrate that their possible mechanism of action could be by interacting with falcipains and their plasmodial homologs. Dynamic residue network (DRN) analysis on the plasmodial proteases identified functionally important residues, including a region with high betweenness centrality, which had previously been proposed as a potential allosteric site in FP-2.

In silico Guided Drug Repurposing: Discovery of New Competitive and Non-competitive Inhibitors of Falcipain-2

Malaria is among the leading causes of death worldwide. The emergence of Plasmodium falciparum resistant strains with reduced sensitivity to the first line combination therapy and suboptimal responses to insecticides used for Anopheles vector management have led to renewed interest in novel therapeutic options. Here, we report the development and validation of an ensemble of ligand-based computational models capable of identifying falcipain-2 inhibitors, and their subsequent application in the virtual screening of DrugBank and Sweetlead libraries. Among four hits submitted to enzymatic assays, two (odanacatib, an abandoned investigational treatment for osteoporosis and bone metastasis, and the antibiotic methacycline) confirmed inhibitory effects on falcipain-2, with Ki of 98.2 nM and 84.4 μM. Interestingly, Methacycline proved to be a non-competitive inhibitor (α = 1.42) of falcipain-2. The effects of both hits on falcipain-2 hemoglobinase activity and on the development of P. falciparum were also studied.

Dissecting a novel allosteric mechanism of cruzain: A computer-aided approach

Trypanosoma cruzi is the causative agent of Chagas disease, a neglected infection affecting millions of people in tropical regions. There are several chemotherapeutic agents for the treatment of this disease, but most of them are highly toxic and generate resistance. Currently, the development of allosteric inhibitors constitutes a promising research field, since it can improve the accessibility to more selective and less toxic medicines. To date, the allosteric drugs prediction is a state-of-the-art topic in rational structure-based computational design. In this work, a simulation strategy was developed for computational discovery of allosteric inhibitors, and it was applied to cruzain, a promising target and the major cysteine protease of T. cruzi. Molecular dynamics simulations, binding free energy calculations and network-based modelling of residue interactions were combined to characterize and compare molecular distinctive features of the apo form and the cruzain-allosteric inhibitor complexes. By using geometry-based criteria on trajectory snapshots, we predicted two main allosteric sites suitable for drug targeting. The results suggest dissimilar mechanisms exerted by the same allosteric site when binding different potential allosteric inhibitors. Finally, we identified the residues involved in suboptimal paths linking the identified site and the orthosteric site. The present study constitutes the first approximation to the design of cruzain allosteric inhibitors and may serve for future pharmacological intervention. Here, no major effects on active site structure were observed due to compound binding (modification of distance and angles between catalytic residues), which indicates that allosteric regulation in cruzain might be mediated via alterations of its dynamical properties similarly to allosteric inhibition of human cathepsin K (HCatK). The current findings are particularly relevant for the design of allosteric modulators of papain-like cysteine proteases.

Synthesis of pyrazole acrylic acid based oxadiazole and amide derivatives as antimalarial and anticancer agents

Depravity of malaria in terms of morbidity and mortality in human beings makes it a major health issue in tropical and subtropical areas of the globe. Drug counterfeiting and non-adherence to the treatment regimen have significantly contributed to development and spread of multidrug resistance that has highlighted the need for development of novel and more efficient antimalarial drugs. Complexity associated with cancer disease and prevalence of diversified cell populations vindicates highly specific treatment options for treatment of cancer. Resistance to these anticancer agents has posed a great hindrance in successful treatment of cancer. Pondering this ongoing situation, it was speculated to develop novel compounds targeting malaria and cancer. Moving on the same aisle, we synthesized pyrazole acrylic acid based oxadiazole and amide derivatives using multi-step reaction pathways (6a-x; 6a'-h'). Schizont maturation inhibition assay was employed to determine antimalarial potential. Compound 6v emerged as the most potent antimalarial agent targeting falcipain-2 enzyme. Anticancer activity was done using sulforhodamine B assay. Compounds 6b' and 6g' demonstrated promising results against all the tested cell lines. Further, Microscopic view clearly indicated formation of apoptotic bodies, chromatin condensation, shrinkage of cells and bleb formation. Validation of the results was achieved using molecular docking studies. From the obtained results, it was observed that cyclization (oxadiazole) favored antimalarial activity while non-cyclized compounds (amides) emerged as better anticancer agents.

Structural and functional attributes of malaria parasite diadenosine tetraphosphate hydrolase

Malaria symptoms are driven by periodic multiplication cycles of Plasmodium parasites in human red blood corpuscles (RBCs). Malaria infection still accounts for ~600,000 annual deaths, and hence discovery of both new drug targets and drugs remains vital. In the present study, we have investigated the malaria parasite enzyme diadenosine tetraphosphate (Ap4A) hydrolase that regulates levels of signalling molecules like Ap4A by hydrolyzing them to ATP and AMP. We have tracked the spatial distribution of parasitic Ap4A hydrolase in infected RBCs, and reveal its unusual localization on the infected RBC membrane in subpopulation of infected cells. Interestingly, enzyme activity assays reveal an interaction between Ap4A hydrolase and the parasite growth inhibitor suramin. We also present a high resolution crystal structure of Ap4A hydrolase in apo- and sulphate- bound state, where the sulphate resides in the enzyme active site by mimicking the phosphate of substrates like Ap4A. The unexpected infected erythrocyte localization of the parasitic Ap4A hydrolase hints at a possible role of this enzyme in purinerigic signaling. In addition, atomic structure of Ap4A hydrolase provides insights for selective drug targeting.

Allosteric regulation of the Plasmodium falciparum cysteine protease falcipain-2 by heme

During the erythrocytic cycle of Plasmodium falciparum malaria parasites break down host hemoglobin, resulting in the release of free heme (ferriprotoporphyrin IX). Heme is a generator of free radicals that cause oxidative stress, but it is detoxified by crystallization into hemozoin inside the food vacuole. We evaluated the interaction of heme and heme analogues with falcipain-2, a P. falciparum food vacuole cysteine protease that plays a key role in hemoglobin digestion. Heme bound to falcipain-2 with a 1:1 stoichiometry, and heme inhibited falcipain-2 activity against both human hemoglobin and chromogenic peptide substrates through a noncompetitive-like mechanism. A series of porphyrin analogues was screened for inhibition of falcipain-2, demonstrating a minor contribution of iron to heme-falcipain-2 interaction, and revealing dependence on both propionic and vinyl groups for inhibition of falcipain-2 by heme. Docking and molecular dynamics simulation unveiled a novel, inducible heme-binding moiety in falcipain-2 adjacent to the catalytic site. Kinetic data suggested that the noncompetitive-like inhibition was substrate inhibition induced by heme. Collectively these data suggest that binding of heme to falcipain-2 may limit the accumulation of free heme in the parasite food vacuole, providing a means of heme detoxification in addition to crystallization into hemozoin.

DOCLASP - Docking ligands to target proteins using spatial and electrostatic congruence extracted from a known holoenzyme and applying simple geometrical transformations

The ability to accurately and effectively predict the interaction between proteins and small drug-like compounds has long intrigued researchers for pedagogic, humanitarian and economic reasons. Protein docking methods (AutoDock, GOLD, DOCK, FlexX and Glide to name a few) rank a large number of possible conformations of protein-ligand complexes using fast algorithms. Previously, it has been shown that structural congruence leading to the same enzymatic function necessitates the congruence of electrostatic properties (CLASP). The current work presents a methodology for docking a ligand into a target protein, provided that there is at least one known holoenzyme with ligand bound - DOCLASP (Docking using CLASP). The contact points of the ligand in the holoenzyme defines a motif, which is used to query the target enzyme using CLASP. If there are significant matches, the holoenzyme and the target protein are superimposed based on congruent atoms. The same linear and rotational transformations are also applied to the ligand, thus creating a unified coordinate framework having the holoenzyme, the ligand and the target enzyme. In the current work, the dipeptidyl peptidase-IV inhibitor vildagliptin was docked to the PI-PLC structure complexed with myo-inositol using DOCLASP. Also, corroboration of the docking of phenylthiourea to the modelled structure of polyphenol oxidase (JrPPO1) from walnut is provided based on the subsequently solved structure of JrPPO1 (PDBid:5CE9). Analysis of the binding of the antitrypanosomial drug suramin to nine non-homologous proteins in the PDB database shows a diverse set of binding motifs, and multiple binding sites in the phospholipase A2-likeproteins from the Bothrops genus of pitvipers. The conformational changes in the suramin molecule on binding highlights the challenges in docking flexible ligands into an already 'plastic' binding site. Thus, DOCLASP presents a method for 'soft docking' ligands to proteins with low computational requirements.

DOCLASP - Docking ligands to target proteins using spatial and electrostatic congruence extracted from a known holoenzyme and applying simple geometrical transformations

The ability to accurately and effectively predict the interaction between proteins and small drug-like compounds has long intrigued researchers for pedagogic, humanitarian and economic reasons. Protein docking methods (AutoDock, GOLD, DOCK, FlexX and Glide to name a few) rank a large number of possible conformations of protein-ligand complexes using fast algorithms. Previously, it has been shown that structural congruence leading to the same enzymatic function necessitates the congruence of electrostatic properties (CLASP). The current work presents a methodology for docking a ligand into a target protein, provided that there is at least one known holoenzyme with ligand bound - DOCLASP (Docking using CLASP). The contact points of the ligand in the holoenzyme defines a motif, which is used to query the target enzyme using CLASP. If there are significant matches, the holoenzyme and the target protein are superimposed based on congruent atoms. The same linear and rotational transformations are also applied to the ligand, thus creating a unified coordinate framework having the holoenzyme, the ligand and the target enzyme. In the current work, the dipeptidyl peptidase-IV inhibitor vildagliptin was docked to the PI-PLC structure complexed with myo-inositol using DOCLASP. Also, corroboration of the docking of phenylthiourea to the modelled structure of polyphenol oxidase (JrPPO1) from walnut is provided based on the subsequently solved structure of JrPPO1 (PDBid:5CE9). Analysis of the binding of the antitrypanosomial drug suramin to nine non-homologous proteins in the PDB database shows a diverse set of binding motifs, and multiple binding sites in the phospholipase A2-likeproteins from the Bothrops genus of pitvipers. The conformational changes in the suramin molecule on binding highlights the challenges in docking flexible ligands into an already ’plastic’ binding site. Thus, DOCLASP presents a method for ’soft docking’ ligands to proteins with low computational requirements.

Synthetic compounds from an in house library as inhibitors of falcipain-2 from Plasmodium falciparum

Falcipain-2 (FP-2) is a key cysteine protease from the malaria parasite Plasmodium falciparum. Many previous studies have identified FP-2 inhibitors; however, none has yet met the criteria for an antimalarial drug candidate. In this work, we assayed an in-house library of non-peptidic organic compounds, including (E)-chalcones, (E)-N'-benzylidene-benzohydrazides and alkyl-esters of gallic acid, and assessed the activity toward FP-2 and their mechanisms of inhibition. The (E)-chalcones 48, 54 and 66 showed the lowest IC50 values (8.5 ± 0.8 µM, 9.5 ± 0.2 µM and 4.9 ± 1.3 µM, respectively). The best inhibitor (compound 66) demonstrated non-competitive inhibition, and using mass spectrometry and fluorescence spectroscopy assays, we suggest a potential allosteric site for the interaction of this compound, located between the catalytic site and the hemoglobin binding arm in FP-2. We combined structural biology tools and mass spectrometry to characterize the inhibition mechanisms of novel compounds targeting FP-2.

DOCLASP - Docking ligands to target proteins using spatial and electrostatic congruence extracted from a known holoenzyme and applying simple geometrical transformations

The ability to accurately and effectively predict the interaction between proteins and small drug-like compounds has long intrigued researchers for pedagogic, humanitarian and economic reasons. Protein docking methods (AutoDock, GOLD, DOCK, FlexX and Glide to name a few) rank a large number of possible conformations of protein-ligand complexes using fast algorithms. Previously, it has been shown that structural congruence leading to the same enzymatic function necessitates the congruence of electrostatic properties (CLASP). The current work presents a methodology for docking a ligand into a target protein, provided that there is at least one known holoenzyme with ligand bound - DOCLASP (Docking using CLASP). The contact points of the ligand in the holoenzyme defines a motif, which is used to query the target enzyme using CLASP. If there are significant matches, the holoenzyme and the target protein are superimposed based on congruent atoms. The same linear and rotational transformations are also applied to the ligand, thus creating a unified coordinate framework having the holoenzyme, the ligand and the target enzyme. In the current work, the dipeptidyl peptidase-IV inhibitor vildagliptin was docked to the PI-PLC structure complexed with myo-inositol using DOCLASP. Also, corroboration of the docking of phenylthiourea to the modelled structure of polyphenol oxidase (JrPPO1) from walnut is provided based on the subsequently solved structure of JrPPO1 (PDBid:5CE9). Analysis of the binding of the antitrypanosomial drug suramin to nine non-homologous proteins in the PDB database shows a diverse set of binding motifs, and multiple binding sites in the phospholipase A2-likeproteins from the Bothrops genus of pitvipers. The conformational changes in the suramin molecule on binding highlights the challenges in docking flexible ligands into an already 'plastic' binding site. Thus, DOCLASP presents a method for 'soft docking' ligands to proteins with low computational requirements.