Crystal structures of the histo-aspartic protease (HAP) from Plasmodium falciparum

Post-Translational Modifications of Proteins of Malaria Parasites during the Life Cycle

Post-translational modifications (PTMs) are essential for regulating protein functions, influencing various fundamental processes in eukaryotes. These include, but are not limited to, cell signaling, protein trafficking, the epigenetic control of gene expression, and control of the cell cycle, as well as cell proliferation, differentiation, and interactions between cells. In this review, we discuss protein PTMs that play a key role in the malaria parasite biology and its pathogenesis. Phosphorylation, acetylation, methylation, lipidation and lipoxidation, glycosylation, ubiquitination and sumoylation, nitrosylation and glutathionylation, all of which occur in malarial parasites, are reviewed. We provide information regarding the biological significance of these modifications along all phases of the complex life cycle of Plasmodium spp. Importantly, not only the parasite, but also the host and vector protein PTMs are often crucial for parasite growth and development. In addition to metabolic regulations, protein PTMs can result in epitopes that are able to elicit both innate and adaptive immune responses of the host or vector. We discuss some existing and prospective results from antimalarial drug discovery trials that target various PTM-related processes in the parasite or host.

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.

Potentials of Terpenoids as Inhibitors of Multiple Plasmodium falciparum Protein Drug Targets

PURPOSE

The resistance of parasite to readily affordable antimalarial drugs, the high cost of currently potent drugs, and the resistance of vector mosquitoes to insecticides threaten the possibility of malaria eradication in malaria endemic areas. Due to the fact that quinine and artemisinin were isolated from plants sources, researchers have been encouraged to search for new antimalarials from medicinal plants. This is especially the case in Africa where a large percentage of the population depends on medicinal plant to treat malaria and other ailments.

METHOD

In this study, we evaluated previously characterized Plasmodium-cidal compounds obtained from the African flora to identify their likely biochemical targets, for an insight into their possible antimalarial chemotherapy. Molecular docking study was first conducted, after which remarkable compounds were submitted for molecular dynamic (MD) simulations studies.

RESULTS

From a total of 38 Plasmodium-cidal compounds docked with confirmed Plasmodium falciparum protein drug targets [plasmepsin II (PMII), histo-aspartic protein (HAP) and falcipain-2 (FP)], two pentacyclic triterpene, cucurbitacin B and 3 beta-O-acetyl oleanolic acid showed high binding affinity relative to artesunate. This implies their capacity to inhibit the three selected P. falciparum target proteins, and consequently, antimalarial potential. From the MD simulations studies and binding free energy outcomes, results confirmed that the two compounds are stable in complex with the selected antimalarial targets; they also showed excellent binding affinities during the 100 ns simulation.

CONCLUSION

These results showed that cucurbitacin B and 3 beta-O-acetyl oleanolic acid are potent antimalarials and should be considered for further studies.

Computational insights for predicting the binding and selectivity of peptidomimetic plasmepsin IV inhibitors against cathepsin D

Plasmepsins (Plms) are aspartic proteases involved in the degradation of human hemoglobin by P. falciparum and are essential for the survival and growth of the parasite. Therefore, Plm enzymes are reported as an important antimalarial drug target. Herein, we have applied molecular docking, molecular dynamics (MD) simulations, and binding free energy with the Linear Interaction Energy (LIE) approach to investigate the binding of peptidomimetic PlmIV inhibitors with a particular focus on understanding their selectivity against the human Asp protease cathepsin D (CatD). The residual decomposition analysis results suggest that amino acid differences in the subsite S3 of PlmIV and CatD are responsible for the higher selectivity of the 5a inhibitor. These findings yield excellent agreement with experimental binding data and provide new details regarding van der Waals and electrostatic interactions of subsite residues as well as structural properties of the PlmIV and CatD systems.

Phenylacetyl pepstatin inhibitors of aspartyl proteases from Streptomyces varsoviensis

A new pepstatin with a phenylacetyl group, pepstatin Pa (1), and its methyl ester (2) were isolated from Streptomyces varsoviensis DSM 40346. Their structures were determined by high-resolution mass spectrometry and nuclear magnetic resonance techniques. The absolute configuration was determined using the Marfey's method. Both pentapeptide products are inhibitors of pepsin and cathepsin D. Interestingly, the bacterial genome contains no biosynthetic gene cluster for the new pepstatin, suggesting an extrachromosomal origin of the biosynthetic genes.

Evaluation of antiplasmodial activity in silico and in vitro of N-acylhydrazone derivatives

N-acylhydrazones are considered privileged structures in medicinal chemistry, being part of antimicrobial compounds (for example). In this study we show the activity of N-acylhydrazone compounds, namely AH1, AH2, AH4, AH5 in in vitro tests against the chloroquine-resistant strain of Plasmodium falciparum (W2) and against WI26 VA-4 human cell lines. All compounds showed low cytotoxicity (LC₅₀ > 100 µM). The AH5 compound was the most active against Plasmodium falciparum, with an IC₅₀ value of 0.07 μM. AH4 and AH5 were selected among the tested compounds for molecular docking calculations to elucidate possible targets involved in their mechanism of action and the SwissADME analysis to predict their pharmacokinetic profile. The AH5 compound showed affinity for 12 targets with low selectivity, while the AH4 compound had greater affinity for only one target (3PHC). These compounds met Lipinski's standards in the ADME in silico tests, indicating good bioavailability results. These results demonstrate that these N-acylhydrazone compounds are good candidates for future preclinical studies against malaria.

Methods Used to Investigate the Plasmodium falciparum Digestive Vacuole

Plasmodium falciparum malaria remains a global health problem as parasites continue to develop resistance to all antimalarials in use. Infection causes clinical symptoms during the intra-erythrocytic stage of the lifecycle where the parasite infects and replicates within red blood cells (RBC). During this stage, P. falciparum digests the main constituent of the RBC, hemoglobin, in a specialized acidic compartment termed the digestive vacuole (DV), a process essential for survival. Many therapeutics in use target one or multiple aspects of the DV, with chloroquine and its derivatives, as well as artemisinin, having mechanisms of action within this organelle. In order to better understand how current therapeutics and those under development target DV processes, techniques used to investigate the DV are paramount. This review outlines the involvement of the DV in therapeutics currently in use and focuses on the range of techniques that are currently utilized to study this organelle including microscopy, biochemical analysis, genetic approaches and metabolomic studies. Importantly, continued development and application of these techniques will aid in our understanding of the DV and in the development of new therapeutics or therapeutic partners for the future.

A proteomic glimpse into the effect of antimalarial drugs on Plasmodium falciparum proteome towards highlighting possible therapeutic targets

There is no effective vaccine against malaria; therefore, chemotherapy is to date the only choice to fight against this infectious disease. However, there is growing evidences of drug-resistance mechanisms in malaria treatments. Therefore, the identification of new drug targets is an urgent need for the clinical management of the disease. Proteomic approaches offer the chance of determining the effects of antimalarial drugs on the proteome of Plasmodium parasites. Accordingly, we reviewed the effects of antimalarial drugs on the Plasmodium falciparum proteome pointing out the relevance of several proteins as possible drug targets in malaria treatment. In addition, some of the P. falciparum stage-specific altered proteins and parasite-host interactions might play important roles in pathogenicity, survival, invasion and metabolic pathways and thus serve as potential sources of drug targets. In this review, we have identified several proteins, including thioredoxin reductase, helicases, peptidyl-prolyl cis-trans isomerase, endoplasmic reticulum-resident calcium-binding protein, choline/ethanolamine phosphotransferase, purine nucleoside phosphorylase, apical membrane antigen 1, glutamate dehydrogenase, hypoxanthine guanine phosphoribosyl transferase, heat shock protein 70x, knob-associated histidine-rich protein and erythrocyte membrane protein 1, as promising antimalarial drugs targets. Overall, proteomic approaches are able to partially facilitate finding possible drug targets. However, the integration of other 'omics' and specific pharmaceutical techniques with proteomics may increase the therapeutic properties of the critical proteins identified in the P. falciparum proteome.

Malaria parasite plasmepsins: More than just plain old degradative pepsins

Plasmepsins are a group of diverse aspartic proteases in the malaria parasite Plasmodium Their functions are strikingly multifaceted, ranging from hemoglobin degradation to secretory organelle protein processing for egress, invasion, and effector export. Some, particularly the digestive vacuole plasmepsins, have been extensively characterized, whereas others, such as the transmission-stage plasmepsins, are minimally understood. Some (e.g. plasmepsin V) have exquisite cleavage sequence specificity; others are fairly promiscuous. Some have canonical pepsin-like aspartic protease features, whereas others have unusual attributes, including the nepenthesin loop of plasmepsin V and a histidine in place of a catalytic aspartate in plasmepsin III. We have learned much about the functioning of these enzymes, but more remains to be discovered about their cellular roles and even their mechanisms of action. Their importance in many key aspects of parasite biology makes them intriguing targets for antimalarial chemotherapy. Further consideration of their characteristics suggests that some are more viable drug targets than others. Indeed, inhibitors of invasion and egress offer hope for a desperately needed new drug to combat this nefarious organism.

Activation mechanism of plasmepsins, pepsin-like aspartic proteases from Plasmodium, follows a unique trans-activation pathway

Plasmodium parasites that cause malaria produce plasmepsins (PMs), pepsin-like aspartic proteases that are important antimalarial drug targets due to their role in host hemoglobin degradation. The enzymes are synthesized as inactive zymogens (pro-PMs), and the mechanism of their conversion to the active, mature forms has not been clearly elucidated. Our structural investigations of vacuolar pro-PMs with truncated prosegment (pro-tPMs) reveal that the formation of the S-shaped dimer is their innate property. Further structural studies, biochemical analysis, and molecular dynamics simulations indicate that disruption of the Tyr-Asp loop (121p-4), coordinated with the movement of the loop L1 (237-247) and helix H2 (101p-113p), is responsible for the extension of the pro-mature region (harboring the cleavage site). Consequently, under acidic pH conditions, these structural changes result in the dissociation of the dimers to monomers and the protonation of the residues in the prosegment prompts its unfolding. Subsequently, we demonstrated that the active site of the monomeric pro-tPMs with the unfolded prosegment is accessible for peptide substrate binding; in contrast, the active site is blocked in folded prosegment form of pro-tPMs. Thus, we propose a novel mechanism of auto-activation of vacuolar pro-tPMs that under acidic conditions can form a catalytically competent active site. One monomer cleaves the prosegment of the other one through a trans-activation process, resulting in formation of mature enzyme. As a result, once a mature enzyme is generated, it leads to the complete conversion of all the inactive pro-tPMs to their mature form. DATABASE: Atomic coordinates and structure factors have been submitted in the Protein Data Bank (PDB) under the PDB IDs 6KUB, 6KUC, and 6KUD.

Molecular docking and dynamics studies on novel benzene sulfonamide substituted pyrazole-pyrazoline analogues as potent inhibitors of Plasmodium falciparum Histo aspartic protease

Malaria is the major health issue in African, Asian and Mediterranean regions of the world. Due to the emerging resistance by the parasites and mosquitoes for the current medications and insecticides, respectively, the malaria free human world can be attained only by the novel design and development of new anti-malarial drugs. Hence, we attempted to carry out in silico screening of benzene sulfonamide substituted pyrazole-pyrazoline series against Histo aspartic protease. Our results reveal that the 65% of the data set with the free binding energy in the range of -11.58 to -11.21 kcal/mol, which is categorized as 'high scoring'. Ligands are docked with the catalytic residues Asp 215, Ser 75, Thr 33 and Ala 217, respectively. Molecular dynamic simulation study of free enzyme and the enzyme complex with 4-(5-(4-methoxyphenyl)-1'phenyl-3'-(p-tolyl)-3,4-1'H,2H-[3,4'-bipyrazol]-2-yl)benezenesulfonamide indicated structural stability. The trajectory analysis of complex reveals that the HAP-ligand complex is more stable than the free HAP. We are of the opinion that our results will be useful for designing potential anti-malarial compounds. AbbreviationsADTauto dock toolsBSPPbenzene sulfonamide substituted pyrazole-pyrazolineCQchloroquineHAPhisto aspartic proteaseKKelvinMDmolecular dockingMM/PBSAmolecular mechanics/Poisson Boltzmann surface areaNVTnormal volume and temperatureNPTnormal pressure and temperatureNsnanosecondsPDBprotein data bank.pdbprogram data base formatP. falciparumPlasmodium falciparumPspicosecondsPMsplasmepsinsP. vivaxPlasmodium vivaxRgradius of gyrationRMSDroot mean square deviationRMSFroot mean square fluctuationWHOWorld Health OrganizationCommunicated by Ramaswamy H. Sarma.

Acquired Removability of Aspartic Protease Inhibitors by Direct Biotinylation

Protease inhibitors are used as both research tools and therapeutics. Many of these inhibitors consist of substrate amino acid sequence-derived structure with a transition state mimic to interact with the active site of the protease, suppressing enzymatic activity. However, once they bind, macrodilution or protein denaturation is required to remove them, limiting their usage. In this study, we describe a removable protease inhibitor, which is a directly biotinylated analogue to control the activities of HIV-1 protease and human cathepsin D. In the substrate cleavage assay, we observed that the nanomolar inhibitory activities were lost upon the addition of streptavidin, while the enzymatic activities sufficiently recovered. HIV-1 protease mixed with the removable inhibitor, avoiding autolysis, was still active to be detected by adding streptavidin after one year at room temperature. We also observed that the inhibitor was an effective eluent for the simple detection of the activity of proteases purified from human serum and cells. These results demonstrate that direct biotinylation of protease inhibitors could be a novel method for controlling the enzymatic activity from OFF to ON. We proposed the phenomenon that binding equilibrium of inhibitor was shifted from protease to streptavidin with higher affinity, named "inhibitor stripping action by affinity competition", or ISAAC. We anticipate that ISAAC could be applicable for preservatives of proteases and activity-based diagnosis of protease related diseases. Furthermore, removable inhibitor to be designed for targeted proteases changing the inhibitor structure may elucidate enzymatic activity in intrinsic form with natural modifications from various biological samples.

FAM105A/OTULINL Is a Pseudodeubiquitinase of the OTU-Class that Localizes to the ER Membrane

Pseudoenzymes have been identified across a diverse range of enzyme classes and fulfill important cellular functions. Examples of pseudoenzymes exist within ubiquitin conjugating and deubiquitinase (DUB) protein families. Here we characterize FAM105A/OTULINL, the only putative pseudodeubiquitinase of the ovarian tumor protease (OTU domain) family in humans. The crystal structure of FAM105A revealed that the OTU domain possesses structural deficiencies in both active site and substrate-binding infrastructure predicted to impair normal DUB function. We confirmed the absence of catalytic function against all ubiquitin linkages and an inability of FAM105A to bind ubiquitin compared with catalytically active FAM105B/OTULIN. FAM105A co-localized with KDEL markers and Lamin B1 at the endoplasmic reticulum (ER) and nuclear envelope, respectively. Accordingly, the FAM105A interactome exhibited significant enrichment in proteins localized to the ER/outer nuclear, Golgi and vesicular membranes. In light of undetectable deubiquitinase activity, we posit that FAM105A/OTULINL functions through its ability to mediate protein-protein interactions.

Overexpression of plasmepsin II and plasmepsin III does not directly cause reduction in Plasmodium falciparum sensitivity to artesunate, chloroquine and piperaquine

Artemisinin derivatives and their partner drugs in artemisinin combination therapies (ACTs) have played a pivotal role in global malaria mortality reduction during the last two decades. The loss of artemisinin efficacy due to evolving drug-resistant parasites could become a serious global health threat. Dihydroartemisinin-piperaquine is a well tolerated and generally highly effective ACT. The implementation of a partner drug in ACTs is critical in the control of emerging artemisinin resistance. Even though artemisinin is highly effective in parasite clearance, it is labile in the human body. A partner drug is necessary for killing the remaining parasites when the pulses of artemisinin have ceased. A population of Plasmodium falciparum parasites in Cambodia and adjacent countries has become resistant to piperaquine. Increased copy number of the genes encoding the haemoglobinases Plasmepsin II and Plasmepsin III has been linked with piperaquine resistance by genome-wide association studies and in clinical trials, leading to the use of increased plasmepsin II/plasmepsin III copy number as a molecular marker for piperaquine resistance. Here we demonstrate that overexpression of plasmepsin II and plasmepsin III in the 3D7 genetic background failed to change the susceptibility of P. falciparum to artemisinin, chloroquine and piperaquine by both a standard dose-response analysis and a piperaquine survival assay. Whilst plasmepsin copy number polymorphism is currently implemented as a molecular surveillance resistance marker, further studies to discover the molecular basis of piperaquine resistance and potential epistatic interactions are needed.

A Perspective on Thiazolidinone Scaffold Development as a New Therapeutic Strategy for Toxoplasmosis

Toxoplasma gondii is one of the most successful parasites due to its ability to infect a wide variety of warm-blooded animals. It is estimated that one-third of the world's population is latently infected. The generic therapy for toxoplasmosis has been a combination of antifolates such as pyrimethamine or trimethoprim with either sulfadiazine or antibiotics such as clindamycin with a combination with leucovorin to prevent hematologic toxicity. This therapy shows limitations such as drug intolerance, low bioavailability or drug resistance by the parasite. There is a need for the development of new molecules with the capacity to block any stage of the parasite's life cycle in humans or in a different type of hosts. Heterocyclic compounds are promissory drugs due to its reported biological activity; for this reason, thiazolidinone and its derivatives are presented as a new alternative not only for its inhibitory activity against the parasite but also for its high selectivity-level with high therapeutic index. Thiazolidinones are an important scaffold known to be associated with anticancer, antibacterial, antifungal, antiviral, antioxidant, and antidiabetic activities. The molecule possesses an imidazole ring that has been described as an antiprotozoal agent with antiparasitic properties and less toxicity. Thiazolidinone derivatives have been reportedly as building blocks in organic chemistry and as scaffolds for drug discovery. Here we present a perspective of how structural modifications of the thiazolidinone core could generate new compounds with high anti-parasitic effect and less toxic results.

Deciphering the mechanism of potent peptidomimetic inhibitors targeting plasmepsins - biochemical and structural insights

Malaria is a deadly disease killing worldwide hundreds of thousands people each year and the responsible parasite has acquired resistance to the available drug combinations. The four vacuolar plasmepsins (PMs) in Plasmodium falciparum involved in hemoglobin (Hb) catabolism represent promising targets to combat drug resistance. High antimalarial activities can be achieved by developing a single drug that would simultaneously target all the vacuolar PMs. We have demonstrated for the first time the use of soluble recombinant plasmepsin II (PMII) for structure-guided drug discovery with KNI inhibitors. Compounds used in this study (KNI-10742, 10743, 10395, 10333, and 10343) exhibit nanomolar inhibition against PMII and are also effective in blocking the activities of PMI and PMIV with the low nanomolar K_i values. The high-resolution crystal structures of PMII-KNI inhibitor complexes reveal interesting features modulating their differential potency. Important individual characteristics of the inhibitors and their importance for potency have been established. The alkylamino analog, KNI-10743, shows intrinsic flexibility at the P2 position that potentiates its interactions with Asp132, Leu133, and Ser134. The phenylacetyl tripeptides, KNI-10333 and KNI-10343, accommodate different ρ-substituents at the P3 phenylacetyl ring that determine the orientation of the ring, thus creating novel hydrogen-bonding contacts. KNI-10743 and KNI-10333 possess significant antimalarial activity, block Hb degradation inside the food vacuole, and show no cytotoxicity on human cells; thus, they can be considered as promising candidates for further optimization. Based on our structural data, novel KNI derivatives with improved antimalarial activity could be designed for potential clinical use. DATABASE: Structural data are available in the PDB under the accession numbers 5YIE, 5YIB, 5YID, 5YIC, and 5YIA.

In vitro and in silico studies of naphthoquinones and peptidomimetics toward Plasmodium falciparum plasmepsin V

Plasmodium proteases play both regulatory and effector roles in essential biological processes in this important pathogen and have long been investigated as drug targets. Plasmepsin V from P. falciparum (PfPMV) is an essential protease that processes proteins for export into the host erythrocyte and is a focus of ongoing drug development efforts. In the present study, recombinant protein production, inhibition assays, binding studies as well as molecular docking and molecular dynamics simulation studies were used to investigate the mode of binding of a PEXEL-based peptidomimetic and naphthoquinone compounds to PfPMV. Consistent with our previous study, refolded PfPMVs were produced with functional characteristics similar to the soluble counterpart. Naphthoquinone compounds inhibited PfPMV activity by 50% at 50 μM but did not affect pepsin activity. The IC₅₀ values of compounds 31 and 37 against PfPMV were 22.25 and 68.94 μM, respectively. Molecular dynamics simulations revealed that PEXEL peptide interacted with PfPMV active site residues via electrostatic interactions while naphthoquinone binding preferred van der Waal interactions. P₁'-Ser of the PfEMP2 substrate formed an additional H-bond with Asp365 promoting the catalytic efficiency. Additionally, the effect of metal ions on the secondary structure of PfPMV was examined. Our results confirmed that Hg²⁺ ions reversibly induced the changes in secondary structure of the protein whereas Fe³⁺ ions induced irreversibly. No change was observed in the presence of Ca²⁺ ions. Overall, the results here suggested that naphthoquinone derivatives may represent another source of antimalarial inhibitors targeting aspartic proteases but further chemical modifications are required.

The Potential of Secondary Metabolites from Plants as Drugs or Leads against Protozoan Neglected Diseases-Part III: In-Silico Molecular Docking Investigations

Malaria, leishmaniasis, Chagas disease, and human African trypanosomiasis continue to cause considerable suffering and death in developing countries. Current treatment options for these parasitic protozoal diseases generally have severe side effects, may be ineffective or unavailable, and resistance is emerging. There is a constant need to discover new chemotherapeutic agents for these parasitic infections, and natural products continue to serve as a potential source. This review presents molecular docking studies of potential phytochemicals that target key protein targets in Leishmania spp., Trypanosoma spp., and Plasmodium spp.

Understanding the structural basis of substrate recognition by Plasmodium falciparum plasmepsin V to aid in the design of potent inhibitors

Plasmodium falciparum plasmepsin V (PfPMV) is an essential aspartic protease required for parasite survival, thus, considered as a potential drug target. This study reports the first detailed structural analysis and molecular dynamics simulation of PfPMV as an apoenzyme and its complexes with the substrate PEXEL as well as with the inhibitor saquinavir. The presence of pro-peptide in PfPMV may not structurally hinder the formation of a functionally competent catalytic active site. The structure of PfPMV-PEXEL complex shows that the unique positions of Glu179 and Gln222 are responsible for providing the specificity of PEXEL substrate with arginine at P3 position. The structural analysis also reveals that the S4 binding pocket in PfPMV is occupied by Ile94, Ala98, Phe370 and Tyr472, and therefore, does not allow binding of pepstatin, a potent inhibitor of most pepsin-like aspartic proteases. Among the screened inhibitors, the HIV-1 protease inhibitors and KNI compounds have higher binding affinities for PfPMV with saquinavir having the highest value. The presence of a flexible group at P2 and a bulky hydrophobic group at P3 position of the inhibitor is preferred in the PfPMV substrate binding pocket. Results from the present study will aid in the design of potent inhibitors of PMV.

Comparative genome-wide analysis and evolutionary history of haemoglobin-processing and haem detoxification enzymes in malarial parasites

Background

Malaria parasites have evolved a series of intricate mechanisms to survive and propagate within host red blood cells. Intra-erythrocytic parasitism requires these organisms to digest haemoglobin and detoxify iron-bound haem. These tasks are executed by haemoglobin-specific proteases and haem biocrystallization factors that are components of a large multi-subunit complex. Since haemoglobin processing machineries are functionally and genetically linked to the modes of action and resistance mechanisms of several anti-malarial drugs, an understanding of their evolutionary history is important for drug development and drug resistance prevention.

Methods

Maximum likelihood trees of genetic repertoires encoding haemoglobin processing machineries within Plasmodium species, and with the representatives of Apicomplexan species with various host tropisms, were created. Genetic variants were mapped onto existing three-dimensional structures. Genome-wide single nucleotide polymorphism data were used to analyse the selective pressure and the effect of these mutations at the structural level.

Results

Recent expansions in the falcipain and plasmepsin repertoires are unique to human malaria parasites especially in the Plasmodium falciparum and P. reichenowi lineage. Expansion of haemoglobin-specific plasmepsins occurred after the separation event of Plasmodium species, but the other members of the plasmepsin family were evolutionarily conserved with one copy for each sub-group in every Apicomplexan species. Haemoglobin-specific falcipains are separated from invasion-related falcipain, and their expansions within one specific locus arose independently in both P. falciparum and P. vivax lineages. Gene conversion between P. falciparum falcipain 2A and 2B was observed in artemisinin-resistant strains. Comparison between the numbers of non-synonymous and synonymous mutations suggests a strong selective pressure at falcipain and plasmepsin genes. The locations of amino acid changes from non-synonymous mutations mapped onto protein structures revealed clusters of amino acid residues in close proximity or near the active sites of proteases.

Conclusion

A high degree of polymorphism at the haemoglobin processing genes implicates an imposition of selective pressure. The identification in recent years of functional redundancy of haemoglobin-specific proteases makes them less appealing as potential drug targets, but their expansions, especially in the human malaria parasite lineages, unequivocally point toward their functional significance during the independent and repetitive adaptation events in malaria parasite evolutionary history.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1097-9) contains supplementary material, which is available to authorized users.

Enzymatic Characterization of Recombinant Food Vacuole Plasmepsin 4 from the Rodent Malaria Parasite Plasmodium berghei

The rodent malaria parasite Plasmodium berghei is a practical model organism for experimental studies of human malaria. Plasmepsins are a class of aspartic proteinase isoforms that exert multiple pathological effects in malaria parasites. Plasmepsins residing in the food vacuole (FV) of the parasite hydrolyze hemoglobin in red blood cells. In this study, we cloned PbPM4, the FV plasmepsin gene of P. berghei that encoded an N-terminally truncated pro-segment and the mature enzyme from genomic DNA. We over-expressed this PbPM4 zymogen as inclusion bodies (IB) in Escherichia coli, and purified the protein following in vitro IB refolding. Auto-maturation of the PbPM4 zymogen to mature enzyme was carried out at pH 4.5, 5.0, and 5.5. Interestingly, we found that the PbPM4 zymogen exhibited catalytic activity regardless of the presence of the pro-segment. We determined the optimal catalytic conditions for PbPM4 and studied enzyme kinetics on substrates and inhibitors of aspartic proteinases. Using combinatorial chemistry-based peptide libraries, we studied the active site preferences of PbPM4 at subsites S1, S2, S3, S1’, S2’ and S3’. Based on these results, we designed and synthesized a selective peptidomimetic compound and tested its inhibition of PbPM4, seven FV plasmepsins from human malaria parasites, and human cathepsin D (hcatD). We showed that this compound exhibited a >10-fold selectivity to PbPM4 and human malaria parasite plasmepsin 4 orthologs versus hcatD. Data from this study furthesr our understanding of enzymatic characteristics of the plasmepsin family and provides leads for anti-malarial drug design.

Flap flexibility amongst plasmepsins I, II, III, IV, and V: Sequence, structural, and molecular dynamics analyses

Herein, for the first time, we comparatively report the opening and closing of apo plasmepsin I - V. Plasmepsins belong the aspartic protease family of enzymes, and are expressed during the various stages of the P. falciparum lifecycle, the species responsible for the most lethal and virulent malaria to infect humans. Plasmepsin I, II, IV and HAP degrade hemoglobin from infected red blood cells, whereas plasmepsin V transport proteins crucial to the survival of the malaria parasite across the endoplasmic reticulum. Flap-structures covering the active site of aspartic proteases (such as HIV protease) are crucial to the conformational flexibility and dynamics of the protein, and ultimately control the binding landscape. The flap-structure in plasmepsins is made up of a flip tip in the N-terminal lying perpendicular to the active site, adjacent to the flexible loop region in the C-terminal. Using molecular dynamics, we propose three parameters to better describe the opening and closing of the flap-structure in apo plasmepsins. Namely, the distance, d1, between the flap tip and the flexible region; the dihedral angle, ϕ, to account for the twisting motion; and the TriCα angle, θ1. Simulations have shown that as the flap-structure twists, the flap and flexible region move apart opening the active site, or move toward each other closing the active site. The data from our study indicate that of all the plasmepsins investigated in the present study, Plm IV and V display the highest conformational flexibility and are more dynamic structures versus Plm I, II, and HAP.

Flap dynamics of plasmepsin proteases: insight into proposed parameters and molecular dynamics

Herein, for the first time, we report the flap opening and closing in Plasmepsin proteases – plasmepsin II (PlmII) was used as a prototype model.

From crystal to compound: structure-based antimalarial drug discovery

Despite a century of control and eradication campaigns, malaria remains one of the world's most devastating diseases. Our once-powerful therapeutic weapons are losing the war against the Plasmodium parasite, whose ability to rapidly develop and spread drug resistance hamper past and present malaria-control efforts. Finding new and effective treatments for malaria is now a top global health priority, fuelling an increase in funding and promoting open-source collaborations between researchers and pharmaceutical consortia around the world. The result of this is rapid advances in drug discovery approaches and technologies, with three major methods for antimalarial drug development emerging: (i) chemistry-based, (ii) target-based, and (iii) cell-based. Common to all three of these approaches is the unique ability of structural biology to inform and accelerate drug development. Where possible, SBDD (structure-based drug discovery) is a foundation for antimalarial drug development programmes, and has been invaluable to the development of a number of current pre-clinical and clinical candidates. However, as we expand our understanding of the malarial life cycle and mechanisms of resistance development, SBDD as a field must continue to evolve in order to develop compounds that adhere to the ideal characteristics for novel antimalarial therapeutics and to avoid high attrition rates pre- and post-clinic. In the present review, we aim to examine the contribution that SBDD has made to current antimalarial drug development efforts, covering hit discovery to lead optimization and prevention of parasite resistance. Finally, the potential for structural biology, particularly high-throughput structural genomics programmes, to identify future targets for drug discovery are discussed.

Disulfide linkages in Plasmodium falciparum plasmepsin-i are essential elements for its processing activity and multi-milligram recombinant production yield

Plasmodium falciparum plasmepsin-I (PM-I) has been considered a potential drug target for the parasite that causes fatal malaria in human. Determination of PM-I structures for rational design of its inhibitors is hindered by the difficulty in obtaining large quantity of soluble enzyme. Nearly all attempts for its heterologous expression in Escherichia coli result in the production of insoluble proteins in both semi-pro-PM-I and its truncated form, and thus require protein refolding. Moreover, the yields of purified, soluble PM-I from all reported studies are very limited. Exclusion of truncated semi-pro-PM-I expression in E. coli C41(DE3) is herein reported. We also show that the low preparation yield of purified semi-pro-PM-I with autoprocessing ability is mainly a result of structural instability of the refolded enzyme in acidic conditions due to incomplete formation of disulfide linkages. Upon formation of at least one of the two natural disulfide bonds, nearly all of the refolded semi-pro-PM-I could be activated to its mature form. A significantly improved yield of 10 mg of semi-pro-PM-I per liter of culture, which resulted in 6-8 mg of the mature PM-I, was routinely obtained using this strategy.

Crystal structure of a putative aspartic proteinase domain of the Mycobacterium tuberculosis cell surface antigen PE_PGRS16

We report the crystal structure of the first prokaryotic aspartic proteinase-like domain identified in the genome of Mycobacterium tuberculosis. A search in the genomes of Mycobacterium species showed that the C-terminal domains of some of the PE family proteins contain two classic DT/SG motifs of aspartic proteinases with a low overall sequence similarity to HIV proteinase. The three-dimensional structure of one of them, Rv0977 (PE_PGRS16) of M. tuberculosis revealed the characteristic pepsin-fold and catalytic site architecture. However, the active site was completely blocked by the N-terminal His-tag. Surprisingly, the enzyme was found to be inactive even after the removal of the N-terminal His-tag. A comparison of the structure with pepsins showed significant differences in the critical substrate binding residues and in the flap tyrosine conformation that could contribute to the lack of proteolytic activity of Rv0977.

Molecular modeling and prediction of binding mode and relative binding affinity of Art-Qui-OH with P. falciparum Histo-Aspartic Protease (HAP)

The relative binding affinity in terms of ΔΔG (bind-cald) value of the antimalarial compound artemisinin-quinine hybrid is primarily derived and is discussed in this article with reference to the ΔG (bind-cald) values of two known inhibitors Pepstatin-A and KNI-10006 complexed with HAP enzyme. The ΔG (bind-cald) value for KNI-10006 and Pepstatin-A is -14.10 kcal/mol and -13.09 kcal/mol respectively. The MM-GB/SA scoring results in the relative binding energy (ΔΔG (bind-cald)) of the hybrid molecule with respect to Pepstatin-A as 2.43 kcal/mol and 3.44 kcal/mol against KNI-10006. The overall binding mode of Art-Qui-OH resembles that of Pepstatin-A binding in HAP active site. We suggest here that the ΔΔG (bind-cald) value & proposed binding mode of the Art-Qui-OH for HAP enzyme should be considered for further structure-based drug design effort.

Structural studies of vacuolar plasmepsins

Plasmepsins (PMs) are pepsin-like aspartic proteases present in different species of parasite Plasmodium. Four Plasmodium spp. (P. vivax, P. ovale, P. malariae, and the most lethal P. falciparum) are mainly responsible for causing human malaria that affects millions worldwide. Due to the complexity and rate of parasite mutation coupled with regional variations, and the emergence of P. falciparum strains which are resistant to antimalarial agents such as chloroquine and sulfadoxine/pyrimethamine, there is constant pressure to find new and lasting chemotherapeutic drug therapies. Since many proteases represent therapeutic targets and PMs have been shown to play an important role in the survival of parasite, these enzymes have recently been identified as promising targets for the development of novel antimalarial drugs. The genome of P. falciparum encodes 10 PMs (PMI, PMII, PMIV-X and histo-aspartic protease (HAP)), 4 of which (PMI, PMII, PMIV and HAP) reside within the food vacuole, are directly involved in degradation of human hemoglobin, and share 50-79% amino acid sequence identity. This review focuses on structural studies of only these four enzymes, including their orthologs in other Plasmodium spp.. Almost all original crystallographic studies were performed with PMII, but more recent work on PMIV, PMI, and HAP resulted in a more complete picture of the structure-function relationship of vacuolar PMs. Many structures of inhibitor complexes of vacuolar plasmepsins, as well as their zymogens, have been reported in the last 15 years. Information gained by such studies will be helpful for the development of better inhibitors that could become a new class of potent antimalarial drugs. This article is part of a Special Issue entitled: Proteolysis 50 years after the discovery of lysosome.

Structural insights into the activation and inhibition of histo-aspartic protease from Plasmodium falciparum

Histo-aspartic protease (HAP) from Plasmodium falciparum is a promising target for the development of novel antimalarial drugs. The sequence of HAP is highly similar to those of pepsin-like aspartic proteases, but one of the two catalytic aspartates, Asp32, is replaced with histidine. Crystal structures of the truncated zymogen of HAP and of the complex of the mature enzyme with inhibitor KNI-10395 have been determined at 2.1 and 2.5 Å resolution, respectively. As in other proplasmepsins, the propeptide of the zymogen interacts with the C-terminal domain of the enzyme, forcing the N- and C-terminal domains apart, thereby separating His32 and Asp215 and preventing formation of the mature active site. In the inhibitor complex, the enzyme forms a tight domain-swapped dimer, not previously seen in any aspartic proteases. The inhibitor is found in an unprecedented conformation resembling the letter U, stabilized by two intramolecular hydrogen bonds. Surprisingly, the location and conformation of the inhibitor are similar to those of the fragment of helix 2 comprising residues 34p-38p in the prosegments of the zymogens of gastric aspartic proteases; a corresponding helix assumes a vastly different orientation in proplasmepsins. Each inhibitor molecule is in contact with two molecules of HAP, interacting with the carboxylate group of the catalytic Asp215 of one HAP protomer through a water molecule, while also making a direct hydrogen bond to Glu278A' of the other protomer. A comparison of the shifts in the positions of the catalytic residues in the inhibitor complex presented here with those published previously gives further hints regarding the enzymatic mechanism of HAP.

Computational perspectives into plasmepsins structure-function relationship: implications to inhibitors design

The development of efficient and selective antimalariais remains a challenge for the pharmaceutical industry. The aspartic proteases plasmepsins, whose inhibition leads to parasite death, are classified as targets for the design of potent drugs. Combinatorial synthesis is currently being used to generate inhibitor libraries for these enzymes, and together with computational methodologies have been demonstrated capable for the selection of lead compounds. The high structural flexibility of plasmepsins, revealed by their X-ray structures and molecular dynamics simulations, made even more complicated the prediction of putative binding modes, and therefore, the use of common computational tools, like docking and free-energy calculations. In this review, we revised the computational strategies utilized so far, for the structure-function relationship studies concerning the plasmepsin family, with special focus on the recent advances in the improvement of the linear interaction estimation (LIE) method, which is one of the most successful methodologies in the evaluation of plasmepsin-inhibitor binding affinity.

Crystal structures of the free and inhibited forms of plasmepsin I (PMI) from Plasmodium falciparum

Plasmepsin I (PMI) is one of the four vacuolar pepsin-like proteases responsible for hemoglobin degradation by the malarial parasite Plasmodium falciparum, and the only one with no crystal structure reported to date. Due to substantial functional redundancy of these enzymes, lack of inhibition of even a single plasmepsin can defeat efforts in creating effective antiparasitic agents. We have now solved crystal structures of the recombinant PMI as apoenzyme and in complex with the potent peptidic inhibitor, KNI-10006, at the resolution of 2.4 and 3.1Å, respectively. The apoenzyme crystallized in the orthorhombic space group P2(1)2(1)2(1) with two molecules in the asymmetric unit and the structure has been refined to the final R-factor of 20.7%. The KNI-10006 bound enzyme crystallized in the tetragonal space group P4(3) with four molecules in the asymmetric unit and the structure has been refined to the final R-factor of 21.1%. In the PMI-KNI-10006 complex, the inhibitors were bound identically to all four enzyme molecules, with the opposite directionality of the main chain of KNI-10006 relative to the direction of the enzyme substrates. Such a mode of binding of inhibitors containing an allophenylnorstatine-dimethylthioproline insert in the P1-P1' positions, previously reported in a complex with PMIV, demonstrates the importance of satisfying the requirements for the proper positioning of the functional groups in the mechanism-based inhibitors towards the catalytic machinery of aspartic proteases, as opposed to binding driven solely by the specificity of the individual enzymes. A comparison of the structure of the PMI-KNI-10006 complex with the structures of other vacuolar plasmepsins identified the important differences between them and may help in the design of specific inhibitors targeting the individual enzymes.

Characterization of the monomer-dimer equilibrium of recombinant histo-aspartic protease from Plasmodium falciparum

Histo-aspartic protease (HAP) from Plasmodium falciparum is an intriguing aspartic protease due to its unique structure. Our previous study reported the first recombinant expression of soluble HAP, in its truncated form (lys77p-Leu328) (p denotes prosegment), as a thioredoxin (Trx) fusion protein Trx-tHAP. The present study found that the recombinant Trx-tHAP fusion protein aggregated during purification which could be prevented through the addition of 0.2% CHAPS. Trx-tHAP fusion protein was processed into a mature form of tHAP (mtHAP) by both autoactivation, and activation with either enterokinase or plasmepsin II. Using gel filtration chromatography as well as sedimentation velocity and equilibrium ultracentrifugation, it was shown that the recombinant mtHAP exists in a dynamic monomer-dimer equilibrium with an increasing dissociation constant in the presence of CHAPS. Enzymatic activity data indicated that HAP was most active as a monomer. The dominant monomeric form showed a K(m) of 2.0 microM and a turnover number, k(cat), of 0.036s(-1) using the internally quenched fluorescent synthetic peptide substrate EDANS-CO-CH(2)-CH(2)-CO-Ala-Leu-Glu-Arg-Met-Phe-Leu-Ser-Phe-Pro-Dap-(DABCYL)-OH (2837b) at pH 5.2.