Abstract
In the present study, we investigated the genetic diversity of Plasmodium vivax metacaspase 1 (PvMCA1) catalytic domain in two municipalities of the main malaria hotspot in Brazil, i.e., the Juruá Valley, and observed complete sequence identity among all P. vivax field isolates and the Sal-1 reference strain. Analysis of PvMCA1 catalytic domain in different P. vivax genomic sequences publicly available also revealed a high degree of conservation worldwide, with very few amino acid substitutions that were not related to putative histidine and cysteine catalytic residues, whose involvement with the active site of protease was herein predicted by molecular modeling. The genetic conservation presented by PvMCA1 may contribute to its eligibility as a druggable target candidate in vivax malaria.
Key words:
P. vivax; metacaspase; genetic diversity
The emergence and spread of drug-resistant parasites constitute major obstacles to malaria control in the world,11. Antony HA, Parija SC. Antimalarial drug resistance: an overview. Trop Parasitol. 2016; 6(1): 30-41. where Plasmodium falciparum and P. vivax are the most important species causing malaria.22. WHO - World Health Organization. World Malaria Report 2019. Geneva: WHO; 2019. Available from: https://www.who.int/publications/i/item/9789241565721. Throughout the history of antimalarial therapy, P. falciparum has been notorious for its capacity to develop resistance to diverse antimalarial drugs, including chloroquine, sulfadoxine/pyrimethamine and, more recently, artemisinin in Southeast Asia.11. Antony HA, Parija SC. Antimalarial drug resistance: an overview. Trop Parasitol. 2016; 6(1): 30-41.,33. Dondorp AM, Smithuis FM, Woodrow C, von Seidlein L. How to contain artemisinin- and multidrug-resistant falciparum malaria. Trends Parasitol. 2017; 33(5): 353-63. Nevertheless, emergence of drug-resistant P. vivax has also been documented since the late 1980s in Papua New Guinea and,44. Renapurkar DM, Pradhan VR, Sutar NK, Deshmukh RA, Pandit CH, Marathe SN. Micro test for assaying sensitivity of Plasmodium vivax in vitro. Chemotherapy. 1989; 35(3): 160-3.,55. Rieckmann KH, Davis DR, Hutton DC. Plasmodium vivax resistance to chloroquine? Lancet. 1989; 2(8673): 1183-4. currently, therapeutic failure of chloroquine, which combined with primaquine comprise the first-line treatment for vivax malaria, is reported in many vivax-endemic countries,66. Price RN, von Seidlein L, Valecha N, Nosten F, Baird JK, White NJ. Global extent of chloroquine-resistant Plasmodium vivax: a systematic review and meta-analysis. Lancet Infect Dis. 2014; 14(10): 982-91.,77. Dayananda KK, Achur RN, Gowda DC. Epidemiology, drug resistance, and pathophysiology of Plasmodium vivax malaria. J Vector Borne Dis. 2018; 55(1): 1-8. while the magnitude of primaquine tolerance remains largely unknown.88. Thomas D, Tazerouni H, Sundararaj KGS, Cooper JC. Therapeutic failure of primaquine and need for new medicines in radical cure of Plasmodium vivax. Acta Trop. 2016; 160: 35-8. These facts highlight the urgent need for the development of novel compounds targeting plasmodial molecular pathways with mechanisms of action different from the classical antimalarial drugs.
In this scenario, proteases playing significant roles in parasite survival have been for long time considered molecular targets for antimalarial drug development and,99. Levy MR, Chou SC. Inhibition of macromolecular synthesis in the malarial parasites by inhibitors of proteolytic enzymes. Experientia. 1975; 31(1): 52-4.,1010. Rosenthal PJ. Proteases of malaria parasites: new targets for chemotherapy. Emerg Infect Dis. 1998; 4(1): 49-57. more recently, members of metacaspase family, whose activity was already suggested to be implicated in both growth and cell death of Plasmodium and other protozoan parasites, have emerged as novel candidates.1111. Kumar B, Mohammad T, Amaduddin, Hussain A, Islam A, Ahmad F, et al. Targeting metacaspase-3 from Plasmodium falciparum towards antimalarial therapy: a combined approach of in-silico and in-vitro investigation. J Biomol Struct Dyn. 2021; 39(2): 421-30.,1212. Vandana, Shankar S, Prasad KM, Kashif M, Kalia I, Rai R, et al. A nonpeptidyl molecule modulates apoptosis-like cell death by inhibiting P. falciparum metacaspase-2. Biochem J. 2020; 477(7): 1323-44. Metacaspases are cysteine proteases belonging to the C14 family that are found in genome of protists, fungi, algae and plants. They present structural similarity to metazoan caspases, both having a conserved His-Cys catalytic dyad in the large (p20) domain, but different substrate specificity.1313. Klemencic M, Funk C. Structural and functional diversity of caspase homologues in non-metazoan organisms. Protoplasma. 2018; 255(1): 387-97. In Plasmodium genomes, three metacaspases (MCA1-3) occur1414. Le Chat L, Sinden RE, Dessens JT. The role of metacaspase 1 in Plasmodium berghei development and apoptosis. Mol Biochem Parasitol. 2007; 153(1): 41-7. and a study performed with field isolates from Mauritania, Sudan and Oman showed that P. vivax metacaspase 1 (PvMCA1) can present single nucleotide polymorphisms in the putative His-Cys catalytic residues,1515. Sow F, Bonnot G, Ahmed BR, Diagana SM, Kebe H, Koita M, et al. Genetic diversity of Plasmodium vivax metacaspase 1 and Plasmodium vivax multi-drug resistance 1 genes of field isolates from Mauritania, Sudan and Oman. Malar J. 2017; 16(1): 61. which in some extent can be a factor limiting the eligibility of PvMCA1 as a novel drug target for vivax malaria. In the present study, therefore, we investigated the genetic diversity of the PvMCA1 catalytic domain in a P. vivax population from the main malaria hotspot in Brazil - a country where this plasmodial species is highly prevalent (~90%).22. WHO - World Health Organization. World Malaria Report 2019. Geneva: WHO; 2019. Available from: https://www.who.int/publications/i/item/9789241565721.,1616. Ferreira MU, Castro MC. Challenges for malaria elimination in Brazil. Malar J. 2016; 15(1): 284.
Study area and sample collection - P. vivax parasite isolates were collected from two municipalities [Cruzeiro do Sul (CZS; 7º37′51″S, 72º40′12″W) and Mâncio Lima (ML; 07º36′50″S, 72º53′4″W)] situated in the Juruá Valley, Northwest of the Acre State, Brazilian Amazon, from June to August 2016 and 2018. Covering together a surface area of 12,597 km2 (CZS: 7,925 km2; ML: 4,672 km2), these municipalities are among the endemic areas with the highest annual parasite incidence (API) (positive blood slides per 1,000 inhabitants) in Brazil, presenting an epidemiological profile with sustained high transmission over the last decade. For reference, the Brazilian Ministry of Health considers high risk areas those with API ≥ 50, and Mâncio Lima and Cruzeiro do Sul registered APIs of 436.4 and 231.9 in 2016 and 422.8 and 147.5 in 2018, respectively. In these areas, eighty-three P. vivax mono-infected individuals living in different localities and presenting uncomplicated malaria were diagnosed through microscopic examination of Giemsa-stained thick blood smears and, then, a single blood sample was collected in EDTA tubes from each individual before initiation of treatment. Samples were centrifuged at 350 g for 10 min to remove plasma and cell pellet was preserved at -20ºC in glycerolyte solution (1:2). Written informed consent was obtained from all donors and the study was reviewed and approved by the Oswaldo Cruz Foundation Ethical Committee and the National Ethical Committee of Brazil (CEP-FIOCRUZ CAAE 46084015.1.0000.5248). All P. vivax infections were later confirmed by polymerase chain reaction (PCR) assay, as previously described.1717. Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, Rosario VE, et al. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993; 61(2): 315-20.
DNA extraction and PCR amplification - parasite DNA was extracted from cryopreserved blood samples by QIAamp DNA blood midi kit (QIAgen), following the manufacturer’s instructions, and stored at -20ºC until use. Gene segment coding for the catalytic domain of PvMCA1 was amplified by standard PCR method using a pair of specific primers (forward, 5′-CATGGAAACAAAAAAAAGG-3′; reverse, 5′-CGAAAACTCCATATCTTTGC-3′), as previously described.1515. Sow F, Bonnot G, Ahmed BR, Diagana SM, Kebe H, Koita M, et al. Genetic diversity of Plasmodium vivax metacaspase 1 and Plasmodium vivax multi-drug resistance 1 genes of field isolates from Mauritania, Sudan and Oman. Malar J. 2017; 16(1): 61. PCR was performed in a Veriti 96-well Thermal Cycler (Applied Biosystems) with a total volume of 25 μL reaction mixture containing 3 μL genomic DNA, 10 pmol/μL each primer, 2.5 unit AmpliTaqTM Gold DNA polymerase, 3 mM MgCl2 and 2 μL of 10X PCR buffer. The following conditions were used: 35 cycles at 95ºC for 30 s, 56ºC for 30 s and 72ºC for 2 min. Amplified products were size-fractionated by electrophoresis within 2% agarose gel (Sigma) containing 0.5 μg/mL of ethidium bromide and, then, were visualised by ultraviolet illumination.
Sequencing and polymorphism analysis - PCR products were purified by Wizard® SV Gel and PCR Clean-Up System (Promega) following the manufacturer’s instructions and sequenced in both directions using above-mentioned primers. The sequencing reaction was performed in duplicate according BigDye™ Terminator v3.1 Cycle Sequencing Kit Applied Biosystems™ using 75-100 ng of the purified PCR products and, then, the obtained products were read on a 3730 xl DNA Analyser (Applied Biosystems). Forward and reverse sequences were compared and checked for quality by using SeqMan v.7.0.0 of the DNASTAR software package (Lasergen, Madison, WI, USA), with default parameters, following manual inspection of chromatograms to eliminate ambiguous bases. A minimum quality score of 20 (base call accuracy ≥ 99%) was considered. Alignment of edited sequences was performed in MEGA 7 using Clustal X2 algorithm to identify polymorphisms relative to PvMCA1 sequence from El Salvador reference strain (Sal-1; PlasmoDB: PVX_114725). Additionally, for worldwide analysis, 112 nucleotide sequences for the complete catalytic domain of PvMCA1, as predicted using Pfam database (https://pfam.xfam.org/), were recovered from PlasmoDB and GenBank databases by BLAST and aligned to search for genetic variability. Metadata obtained were used to identify geographical distribution of the isolates recovered from these databases.
Molecular modeling - The sequence of G0ZIA8 (Sal-1; UniProtKB) was submitted to Blast, in the search of homologous proteins. The template (4AFR, Trypanosoma brucei metacaspase) found in PDB share 38.4% identity with G0ZIA8. After alignment with blast, the coverage is 32%. The e-value was 7e-35. The structural model of G0ZIA8 was constructed with MODELLER (https://salilab.org/modeller/). Fifty models were constructed and the best one was selected according to the DOPE score. The electrostatic surface potential was calculated by APBS software. The structural analysis was performed with Pymol.
To investigate the genetic diversity of PvMCA1, eighty-three wild isolates of P. vivax from Juruá Valley in the Brazilian Amazon were submitted to DNA sequencing for gene segment coding the consensus His372-Cys428 catalytic dyad (13). In contrast to work by Sow et al.,1515. Sow F, Bonnot G, Ahmed BR, Diagana SM, Kebe H, Koita M, et al. Genetic diversity of Plasmodium vivax metacaspase 1 and Plasmodium vivax multi-drug resistance 1 genes of field isolates from Mauritania, Sudan and Oman. Malar J. 2017; 16(1): 61. in which the majority of isolates (24/28, 85.7%) presented amino acid substitution in both catalytic dyad and an upstream cysteine residue (Cys305) also inserted in the peptidase C14 domain, all the isolates herein studied showed complete nucleotide sequence identity each other as well as compared to P. vivax Sal-1 reference, with conservation of all the amino acid residues successfully sequenced (A288 to K446) (Fig. 1), supporting that the proteolytic activity of PvMCA1 may be critical to the parasite.
multiple alignment of PvMCA1 peptidase domain from Brazilian Amazon field isolates and Plasmodium vivax strains from different endemic countries around the world. Deduced amino acid sequences of Peptidase_C14 domain of PvMCA1 were obtained from 83 P. vivax isolates collected in two municipalities of the Juruá Valley (CZS: Cruzeiro do Sul and; Mâncio Lima: ML) and were, then, compared against sequences deduced from 112 P. vivax nucleotide genomic sequences available in GenBank and PlasmoDB, using Sal-1 as reference strain. Selected sequences are shown and the complete analysis is available in Supplementary data. The canonical His372-Cys428 catalytic dyad and the substitutions of amino acid residues are indicated by (*) and (#), respectively. (↓) indicates the adjacent cysteine residue (Cys427) with possible involvement in protease activity. CZS/ML represents all field isolates studied, since complete identity was observed. Monkey-adapted strains are represented by Sal1, BrazilI, MauritaniaI and IndiaVII; clinical isolates by Brazil32, PvC01, PvT01 and PvP01 and; wild ape isolates by PvSY42, PvSY56, Pv104 and Pvl10. (-): indicates non-determined amino acid residues.
Since functional disruption of the single metacaspase (MCA1) in Saccharomyces cerevisiae by mutagenesis was already shown to induce in vitro drug-resistance as a result of refractoriness to apoptosis in the yeasts,1818. Biswas C, Zuo X, Chen SCA, Schibeci SD, Forwood JK, Jolliffe KA, et al. Functional disruption of yeast metacaspase, Mca1, leads to miltefosine resistance and inability to mediate miltefosine-induced apoptotic effects. Fungal Genet Biol. 2014; 67: 71-81. one possibility is that polymorphism of PvMCA1 could be related to chemoresistance phenomena reported in vivax malaria, whose occurrence is reported to be infrequent in Juruá Valley1919. de Oliveira TC, Corder RM, Early A, Rodrigues PT, Ladeia-Andrade S, Alves JMP, et al. Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil. PLoS Negl Trop Dis. 2020; 14(10): e0008808. and that was not characterised in our cross-sectional study. Indeed, metacaspases have been implicated in apoptosis-like cell-death pathways in different microorganisms2020. Tsiatsiani L, Van Breusegem F, Gallois P, Zavialov A, Lam E, Bozhkov PV. Metacaspases. Cell Death Differ. 2011; 18(8): 1279-88. and drug-resistant strains of P. falciparum seem to be insensitive to choroquine-induced apoptosis2121. Picot S, Burnod J, Bracchi V, Chumpitazi BF, Ambroise-Thomas P. Apoptosis related to chloroquine sensitivity of the human malaria parasite Plasmodium falciparum. Trans R Soc Trop Med Hyg. 1997; 91(5): 590-1. as well as to anti-cell death effect promoted by inhibition of metacaspase activity.2222. Meslin B, Beavogui AH, Fasel N, Picot S. Plasmodium falciparum metacaspase PfMCA-1 triggers a z-VAD-fmk inhibitable protease to promote cell death. PLoS One. 2011; 6(8): e23867. Nevertheless, our preliminary study with recurrent infection in five P. vivax patients receiving appropriate antimalarial therapy (chloroquine/primaquine) at a Reference Center for malaria in Manaus city, Amazonas state, does not indicate a relation between PvMCA1 genetic diversity and parasite chemoresistance, as no mutation in the protease catalytic domain was detected in these patients (data not shown). Even so, the conclusion of the survey with a larger number of samples appropriately characterised in terms of parasite drug-sensitivity is needed to better clarify this issue.
In order to know the variability of PvMCA1 worldwide, we extended our analysis to P. vivax isolates from different geographical regions whose genomic sequences are publicly available and, additionally, evaluated complete peptidase C14 domain (K286-S497), as predicted using Pfam database of protein families. Nucleotide sequence analysis of more than one hundred isolates from Americas, Africa, Asia and Indonesia, which included clinical and monkey-adapted strains as well as ape-infecting strains from Cameroon and Gabon (Supplementary data), revealed only four non-synonymous polymorphisms at position 310 (Leu to Val), 360 (Gln to Glu), 457 (Val to Met) and 479 (Lys to Arg) (Fig. 1). Three of them (L310V, Q360E, V457M) were restricted to wild ape parasites, with L310V substitution occurring markedly in all eight isolates analysed, while K479R substitution was the unique variation observed in the clinical and monkey-adapted isolates (Fig. 1 and Supplementary data). It is noteworthy that substitutions were not related to cysteine and histidine residues; even in the ape-infecting P. vivax strains, whose remarkable genetic diversity comparing to human strains suggests distinct demographic histories.2323. Loy DE, Plenderleith LJ, Sundararaman SA, Liu W, Gruszczyk J, Chen YJ, et al. Evolutionary history of human Plasmodium vivax revealed by genome-wide analyses of related ape parasites. Proc Natl Acad Sci USA. 2018; 115(36): E8450-9. These data indicate, therefore, that the peptidase domain of PvMCA1 presents evolutionary conservation of amino acid residues that are putatively involved in metacaspase activity, as described across diverse taxa.2424. Klemencic M, Asplund-Samuelsson J, Dolinar M, Funk C. Phylogenetic distribution and diversity of bacterial pseudo-orthocaspases underline their putative role in photosynthesis. Front Plant Sci. 2019; 10: 293.,2525. Kumar B, Verma S, Kashif M, Sharma R, Atul, Dixit R, et al. Metacaspase-3 of Plasmodium falciparum: an atypical trypsin-like serine protease. Int J Biol Macromol. 2019; 138: 309-20.
Similar results showing few polymorphic nucleotides with high conservation of essential amino acid residues required to catalytic site formation have previously been reported for other P. vivax proteases, such as vivapain-1, -2 and -3, plasmepsin-4 and PvSERA-4.2626. Na BK, Shenai BR, Sijwali PS, Choe Y, Pandey KC, Singh A, et al. Identification and biochemical characterization of vivapains, cysteine proteases of the malaria parasite Plasmodium vivax. Biochem J. 2004; 378(Pt 2): 529-38.,2727. Rawat M, Vijay S, Gupta Y, Dixit R, Tiwari PK, Sharma A. Sequence homology and structural analysis of plasmepsin 4 isolated from Indian Plasmodium vivax isolates. Infect Genet Evol. 2011; 11(5): 924-33.,2828. Rahul CN, Shiva Krishna K, Pawar AP, Bai M, Kumar V, Phadke S, et al. Genetic and structural characterization of PvSERA4: potential implication as therapeutic target for Plasmodium vivax malaria. J Biomol Struct Dyn. 2014; 32(4): 580-90. Curiously, the polymorphisms reported by Sow et al.1515. Sow F, Bonnot G, Ahmed BR, Diagana SM, Kebe H, Koita M, et al. Genetic diversity of Plasmodium vivax metacaspase 1 and Plasmodium vivax multi-drug resistance 1 genes of field isolates from Mauritania, Sudan and Oman. Malar J. 2017; 16(1): 61. in PvMCA1 peptidase domain were restricted to histidine and cysteine residues, including the putative H372-C428 catalytic dyad, which could impact protease activity. Abrogation of proteolytic activity can indeed be achieved by site-specific mutagenesis of histidine and cysteine residues forming consensus catalytic dyad in metacaspases of T. brucei (TbMCA2), Leishmania major (LmjMCA) and Candida albicans (CaMca1).2929. González IJ, Desponds C, Schaff C, Mottram JC, Fasel N. Leishmania major metacaspase can replace yeast metacaspase in programmed cell death and has arginine-specific cysteine peptidase activity. Int J Parasitol. 2007; 37(2): 161-72.,3030. Moss CX, Westrop GD, Juliano L, Coombs GH, Mottram JC. Metacaspase 2 of Trypanosoma brucei is a calcium-dependent cysteine peptidase active without processing. FEBS Lett. 2007; 581(29): 5635-9.,3131. Jeong JH, Lee SE, Kim J. Mutational analysis of metacaspase CaMca1 and decapping activator Edc3 in the pathogenicity of Candida albicans. Fungal Genet Biol. 2016; 97: 18-23. But in other metacaspases this dyad seems not to be essential for catalytic activity, such as in LdMCA1 of L. donovani,3232. Lee N, Gannavaram S, Selvapandiyan A, Debrabant A. Characterization of metacaspases with trypsin-like activity and their putative role in programmed cell death in the protozoan parasite Leishmania. Eukaryot Cell. 2007; 6(10): 1745-57. and a secondary catalytic cysteine has been evidenced in T. congolense (TcoMCA5) as well as in the plant species Arabidopsis thaliana (AtMC9) and Triticum aestivum (TaeMCAII).3333. Belenghi B, Romero-Puertas MC, Vercammen D, Brackenier A, Inzé D, Delledonne M, et al. Metacaspase activity of Arabidopsis thaliana is regulated by S-nitrosylation of a critical cysteine residue. J Biol Chem. 2007; 282(2): 1352-8.,3434. Piszczek E, Dudkiewicz M, Mielecki M. Biochemical and bioinformatic characterization of type II metacaspase protein (TaeMCAII) from wheat. Plant Mol Biol Report. 2012; 30(6): 1338-47.,3535. Eyssen LEA, Coetzer THT. Expression, purification and characterisation of Trypanosoma congolense metacaspase 5 (TcoMCA5) - a potential drug target for animal African trypanosomiasis. Protein Expr Purif. 2019; 164: 105465.
This raises the possibility of alternative amino acid residues in PvMCA1 compensating for the mutations previously reported by Sow and colleagues.1515. Sow F, Bonnot G, Ahmed BR, Diagana SM, Kebe H, Koita M, et al. Genetic diversity of Plasmodium vivax metacaspase 1 and Plasmodium vivax multi-drug resistance 1 genes of field isolates from Mauritania, Sudan and Oman. Malar J. 2017; 16(1): 61. For instance, in MCA1 of P. berghei, P. chabaudi and P. yoelii the predicted catalytic cysteine occurs immediately before of the consensus position, where a substitute residue (proline) is present.1414. Le Chat L, Sinden RE, Dessens JT. The role of metacaspase 1 in Plasmodium berghei development and apoptosis. Mol Biochem Parasitol. 2007; 153(1): 41-7. The same is true for T. brucei MCA4 (TbMCA4), in which the activity of such cysteine was experimentally demonstrated.3636. Szallies A, Kubata BK, Duszenko M. A metacaspase of Trypanosoma brucei causes loss of respiration competence and clonal death in the yeast Saccharomyces cerevisiae. FEBS Lett. 2002; 517(1-3): 144-50. Interestingly, PvMCA1 (Fig. 1; Cys427) and other metacaspases, including MCA1 of P. knowlesi and P. gallinaceum, present a second conserved cysteine adjacently preceding the consensus catalytic cysteine1111. Kumar B, Mohammad T, Amaduddin, Hussain A, Islam A, Ahmad F, et al. Targeting metacaspase-3 from Plasmodium falciparum towards antimalarial therapy: a combined approach of in-silico and in-vitro investigation. J Biomol Struct Dyn. 2021; 39(2): 421-30.,1414. Le Chat L, Sinden RE, Dessens JT. The role of metacaspase 1 in Plasmodium berghei development and apoptosis. Mol Biochem Parasitol. 2007; 153(1): 41-7. and the involvement in proteolytic activity has been shown as well.2929. González IJ, Desponds C, Schaff C, Mottram JC, Fasel N. Leishmania major metacaspase can replace yeast metacaspase in programmed cell death and has arginine-specific cysteine peptidase activity. Int J Parasitol. 2007; 37(2): 161-72.,3030. Moss CX, Westrop GD, Juliano L, Coombs GH, Mottram JC. Metacaspase 2 of Trypanosoma brucei is a calcium-dependent cysteine peptidase active without processing. FEBS Lett. 2007; 581(29): 5635-9.,3535. Eyssen LEA, Coetzer THT. Expression, purification and characterisation of Trypanosoma congolense metacaspase 5 (TcoMCA5) - a potential drug target for animal African trypanosomiasis. Protein Expr Purif. 2019; 164: 105465. Notably, substitution of this alternative cysteine is not found in any of the P. vivax isolates analysed until now (Cys427), indicating by its high degree of conservation a possible role in PvMCA1 activity.
We then predicted the three-dimensional structure of PvMCA1 catalytic domain by homology modeling to examine if C427 and additional histidine residues could, in some extent, participate in catalytic site formation. As shown in Fig. 2, the area harboring the canonical catalytic dyad (His372-Cys428) presented typical electrostatic potential of metacaspase catalytic pocket (Fig. 2B), which is endowed with negative charges and confers specificity toward substrates containing basic arginine or lysine residues,1313. Klemencic M, Funk C. Structural and functional diversity of caspase homologues in non-metazoan organisms. Protoplasma. 2018; 255(1): 387-97.,2424. Klemencic M, Asplund-Samuelsson J, Dolinar M, Funk C. Phylogenetic distribution and diversity of bacterial pseudo-orthocaspases underline their putative role in photosynthesis. Front Plant Sci. 2019; 10: 293. while Cys427 was located in a more positively charged area (Fig. 2B). On the other hand, the proximity of both Cys427 (7.8 Å) and Cys428 (6.1 Å) to His372 residue was shown to be consistent with the active site reported for some cysteine proteases (Fig. 2A), such as legumains, caspases and metacaspases.3434. Piszczek E, Dudkiewicz M, Mielecki M. Biochemical and bioinformatic characterization of type II metacaspase protein (TaeMCAII) from wheat. Plant Mol Biol Report. 2012; 30(6): 1338-47.,3737. Wong AHH, Yan C, Shi Y. Crystal structure of the yeast metacaspase Yca1. J Biol Chem. 2012; 287(35): 29251-9.,3838. Elsässer B, Zauner FB, Messner J, Soh WT, Dall E, Brandstette H. Distinct roles of catalytic cysteine and histidine in the protease and ligase mechanisms of human legumain as revealed by DFT-based QM/MM simulations. ACS Catal. 2017; 7(9): 5585-93. In the single metacaspase of S. cerevisiae (YCA1), for instance, the distance between catalytic residues (Cys-His) was about 9.2 Å, while the nearest histidine (5.8 Å) was not implicated in the proteolytic process,3535. Eyssen LEA, Coetzer THT. Expression, purification and characterisation of Trypanosoma congolense metacaspase 5 (TcoMCA5) - a potential drug target for animal African trypanosomiasis. Protein Expr Purif. 2019; 164: 105465. supporting that Cys427 could participate in the activity of PvMCA1. In this context, an additional histidine residue (His480) was also identified in the vicinity of the canonical catalytic dyad but, unless a conformational change occurs before catalytic site activity, His480 does not seem to constitute a partner residue for Cys427 or Cys 428, since it is >13.0 Å away (Fig. 2A).
prediction of catalytic dyad of PvMCA1 by molecular modeling. (A) Three-dimensional structure of PvMCA1 catalytic domain was modeled using comparative modeling and the distances (in Å) between the putative amino acid residues participating in catalytic dyad formation were estimated considering the Nε atom of histidine (H) imidazole ring and Sγ atom of cysteine (C), which are crucial for the proteolytic process. (B) Electrostatic surface potential analysis of PvMCA1 catalytic domain: negatively charged regions are shown in red and positively regions in blue.
In summary, although mutations in the putative catalytic dyad of PvMCA1 have previously been reported in field isolates from Mauritania, Sudan and Oman, a highly conserved peptidase domain of PvMCA1 was remarkably observed in Brazilian Amazon as well as between P. vivax isolates from different endemic countries around the world, which can contribute to eligibility of this metacaspase as a druggable target candidate in vivax malaria. Nevertheless, further studies aiming biochemical and functional characterisation of PvMCA1, together to identification of amino acid residues that effectively participate in the proteolytic process are still required to elucidate the role of PvMCA1 in P. vivax biology.
ACKNOWLEDGEMENTS
To DNA Sequencing facility of Fiocruz-RJ (P01-001) for the support in sequencing assays.
REFERENCES
-
1Antony HA, Parija SC. Antimalarial drug resistance: an overview. Trop Parasitol. 2016; 6(1): 30-41.
-
2WHO - World Health Organization. World Malaria Report 2019. Geneva: WHO; 2019. Available from: https://www.who.int/publications/i/item/9789241565721.
-
3Dondorp AM, Smithuis FM, Woodrow C, von Seidlein L. How to contain artemisinin- and multidrug-resistant falciparum malaria. Trends Parasitol. 2017; 33(5): 353-63.
-
4Renapurkar DM, Pradhan VR, Sutar NK, Deshmukh RA, Pandit CH, Marathe SN. Micro test for assaying sensitivity of Plasmodium vivax in vitro. Chemotherapy. 1989; 35(3): 160-3.
-
5Rieckmann KH, Davis DR, Hutton DC. Plasmodium vivax resistance to chloroquine? Lancet. 1989; 2(8673): 1183-4.
-
6Price RN, von Seidlein L, Valecha N, Nosten F, Baird JK, White NJ. Global extent of chloroquine-resistant Plasmodium vivax: a systematic review and meta-analysis. Lancet Infect Dis. 2014; 14(10): 982-91.
-
7Dayananda KK, Achur RN, Gowda DC. Epidemiology, drug resistance, and pathophysiology of Plasmodium vivax malaria. J Vector Borne Dis. 2018; 55(1): 1-8.
-
8Thomas D, Tazerouni H, Sundararaj KGS, Cooper JC. Therapeutic failure of primaquine and need for new medicines in radical cure of Plasmodium vivax. Acta Trop. 2016; 160: 35-8.
-
9Levy MR, Chou SC. Inhibition of macromolecular synthesis in the malarial parasites by inhibitors of proteolytic enzymes. Experientia. 1975; 31(1): 52-4.
-
10Rosenthal PJ. Proteases of malaria parasites: new targets for chemotherapy. Emerg Infect Dis. 1998; 4(1): 49-57.
-
11Kumar B, Mohammad T, Amaduddin, Hussain A, Islam A, Ahmad F, et al. Targeting metacaspase-3 from Plasmodium falciparum towards antimalarial therapy: a combined approach of in-silico and in-vitro investigation. J Biomol Struct Dyn. 2021; 39(2): 421-30.
-
12Vandana, Shankar S, Prasad KM, Kashif M, Kalia I, Rai R, et al. A nonpeptidyl molecule modulates apoptosis-like cell death by inhibiting P. falciparum metacaspase-2. Biochem J. 2020; 477(7): 1323-44.
-
13Klemencic M, Funk C. Structural and functional diversity of caspase homologues in non-metazoan organisms. Protoplasma. 2018; 255(1): 387-97.
-
14Le Chat L, Sinden RE, Dessens JT. The role of metacaspase 1 in Plasmodium berghei development and apoptosis. Mol Biochem Parasitol. 2007; 153(1): 41-7.
-
15Sow F, Bonnot G, Ahmed BR, Diagana SM, Kebe H, Koita M, et al. Genetic diversity of Plasmodium vivax metacaspase 1 and Plasmodium vivax multi-drug resistance 1 genes of field isolates from Mauritania, Sudan and Oman. Malar J. 2017; 16(1): 61.
-
16Ferreira MU, Castro MC. Challenges for malaria elimination in Brazil. Malar J. 2016; 15(1): 284.
-
17Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, Rosario VE, et al. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993; 61(2): 315-20.
-
18Biswas C, Zuo X, Chen SCA, Schibeci SD, Forwood JK, Jolliffe KA, et al. Functional disruption of yeast metacaspase, Mca1, leads to miltefosine resistance and inability to mediate miltefosine-induced apoptotic effects. Fungal Genet Biol. 2014; 67: 71-81.
-
19de Oliveira TC, Corder RM, Early A, Rodrigues PT, Ladeia-Andrade S, Alves JMP, et al. Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil. PLoS Negl Trop Dis. 2020; 14(10): e0008808.
-
20Tsiatsiani L, Van Breusegem F, Gallois P, Zavialov A, Lam E, Bozhkov PV. Metacaspases. Cell Death Differ. 2011; 18(8): 1279-88.
-
21Picot S, Burnod J, Bracchi V, Chumpitazi BF, Ambroise-Thomas P. Apoptosis related to chloroquine sensitivity of the human malaria parasite Plasmodium falciparum. Trans R Soc Trop Med Hyg. 1997; 91(5): 590-1.
-
22Meslin B, Beavogui AH, Fasel N, Picot S. Plasmodium falciparum metacaspase PfMCA-1 triggers a z-VAD-fmk inhibitable protease to promote cell death. PLoS One. 2011; 6(8): e23867.
-
23Loy DE, Plenderleith LJ, Sundararaman SA, Liu W, Gruszczyk J, Chen YJ, et al. Evolutionary history of human Plasmodium vivax revealed by genome-wide analyses of related ape parasites. Proc Natl Acad Sci USA. 2018; 115(36): E8450-9.
-
24Klemencic M, Asplund-Samuelsson J, Dolinar M, Funk C. Phylogenetic distribution and diversity of bacterial pseudo-orthocaspases underline their putative role in photosynthesis. Front Plant Sci. 2019; 10: 293.
-
25Kumar B, Verma S, Kashif M, Sharma R, Atul, Dixit R, et al. Metacaspase-3 of Plasmodium falciparum: an atypical trypsin-like serine protease. Int J Biol Macromol. 2019; 138: 309-20.
-
26Na BK, Shenai BR, Sijwali PS, Choe Y, Pandey KC, Singh A, et al. Identification and biochemical characterization of vivapains, cysteine proteases of the malaria parasite Plasmodium vivax. Biochem J. 2004; 378(Pt 2): 529-38.
-
27Rawat M, Vijay S, Gupta Y, Dixit R, Tiwari PK, Sharma A. Sequence homology and structural analysis of plasmepsin 4 isolated from Indian Plasmodium vivax isolates. Infect Genet Evol. 2011; 11(5): 924-33.
-
28Rahul CN, Shiva Krishna K, Pawar AP, Bai M, Kumar V, Phadke S, et al. Genetic and structural characterization of PvSERA4: potential implication as therapeutic target for Plasmodium vivax malaria. J Biomol Struct Dyn. 2014; 32(4): 580-90.
-
29González IJ, Desponds C, Schaff C, Mottram JC, Fasel N. Leishmania major metacaspase can replace yeast metacaspase in programmed cell death and has arginine-specific cysteine peptidase activity. Int J Parasitol. 2007; 37(2): 161-72.
-
30Moss CX, Westrop GD, Juliano L, Coombs GH, Mottram JC. Metacaspase 2 of Trypanosoma brucei is a calcium-dependent cysteine peptidase active without processing. FEBS Lett. 2007; 581(29): 5635-9.
-
31Jeong JH, Lee SE, Kim J. Mutational analysis of metacaspase CaMca1 and decapping activator Edc3 in the pathogenicity of Candida albicans. Fungal Genet Biol. 2016; 97: 18-23.
-
32Lee N, Gannavaram S, Selvapandiyan A, Debrabant A. Characterization of metacaspases with trypsin-like activity and their putative role in programmed cell death in the protozoan parasite Leishmania. Eukaryot Cell. 2007; 6(10): 1745-57.
-
33Belenghi B, Romero-Puertas MC, Vercammen D, Brackenier A, Inzé D, Delledonne M, et al. Metacaspase activity of Arabidopsis thaliana is regulated by S-nitrosylation of a critical cysteine residue. J Biol Chem. 2007; 282(2): 1352-8.
-
34Piszczek E, Dudkiewicz M, Mielecki M. Biochemical and bioinformatic characterization of type II metacaspase protein (TaeMCAII) from wheat. Plant Mol Biol Report. 2012; 30(6): 1338-47.
-
35Eyssen LEA, Coetzer THT. Expression, purification and characterisation of Trypanosoma congolense metacaspase 5 (TcoMCA5) - a potential drug target for animal African trypanosomiasis. Protein Expr Purif. 2019; 164: 105465.
-
36Szallies A, Kubata BK, Duszenko M. A metacaspase of Trypanosoma brucei causes loss of respiration competence and clonal death in the yeast Saccharomyces cerevisiae. FEBS Lett. 2002; 517(1-3): 144-50.
-
37Wong AHH, Yan C, Shi Y. Crystal structure of the yeast metacaspase Yca1. J Biol Chem. 2012; 287(35): 29251-9.
-
38Elsässer B, Zauner FB, Messner J, Soh WT, Dall E, Brandstette H. Distinct roles of catalytic cysteine and histidine in the protease and ligase mechanisms of human legumain as revealed by DFT-based QM/MM simulations. ACS Catal. 2017; 7(9): 5585-93.
-
Financial support: IOC (FIOCRUZ), FAPERJ, CAPES. MVGL and CTDR receive a Productivity Research Fellowship from the CNPq, Brazil; CTDR is a “Cientista do Nosso Estado” from FAPERJ.
Publication Dates
-
Publication in this collection
31 May 2021 -
Date of issue
2021
History
-
Received
13 Nov 2020 -
Accepted
04 May 2021