Falcipain 2 and 3 as malarial drug targets: deciphering the effects of missense mutations and identification of allosteric modulators via computational approaches
Authors:   , Okeke Chiamaka Jessica 
Published: June 2024

Malaria, caused by an obligate unicellular protozoan parasite of the genus Plasmodium, is a disease of global health importance that remains a major cause of morbidity and mortality in developing countries. The World Health Organization (WHO) reported nearly 247 million malaria cases in 2021, causing 619,000 deaths, the vast majority ascribed to pregnant women and young children in sub-Saharan Africa. A critical component of malaria mitigation and elimination efforts worldwide is antimalarial drugs. However, resistance to available antimalarial drugs jeopardizes the treatment, prevention, and eradication of the disease. The recent emergence and spread of resistance to artemisinin (ART), the currently recommended first-line antimalarial drug, emphasizes the need to understand the resistance mechanism and apply this knowledge in developing new drugs that are effective against malaria. An insight into ART's mechanism of action indicates that ferrous iron (Fe2+) or heme, released when hemoglobin is degraded, cleaves the endoperoxide bridge. As a result, free radicals are formed, which alkylate many intracellular targets and result in plasmodial proteopathy. Aside from the existing evidence that mutations in the Kelch 13 protein propeller domain affect ART sensitivity and clearance rate by Plasmodium falciparum (Pf) parasites, recent investigations raise the possibility that additional target loci may be involved, and these include a nonsense (S69stop) and four missense variants (K255R, N257E, T343P, and D345G) in falcipain 2 (FP-2) protein. FP-2 and falcipain 3 (FP-3) are cysteine proteases responsible for hydrolyzing hemoglobin in the host erythrocytic cycle, a key virulence factor for malaria parasite growth and metabolism. Due to the obligatory nature of the hemoglobin degradation process, both proteases have become potential antimalarial drug targets attracting attention in recent years for the development of blood-stage antimalarial drugs. The alteration of the expression profile of FP-2 and FP-3 through gene manipulation approaches (knockout) or compound inhibition assays, respectively, induced parasites with swollen food vacuoles due to the accumulation of undegraded hemoglobin. Furthermore, missense mutations in FP-2 confer parasites with decreased ART sensitivity, probably due to altered enzyme efficiency and momentary decreased hemoglobin degradation. Hence, understanding how these mutations affect FP-2 (including those implicated in ART resistance) and FP-3 is imperative to finding potentially effective inhibitors. The first aim of this thesis is to characterize the effects of missense mutations on the partial zymogen complex and the catalytic domain of FP-2 and FP-3 using a range of computational approaches and tools such as homology modeling, molecular dynamics (MD) simulations, comparative essential dynamics, dynamic residue network (DRN) analysis, weighted residue contact map analysis, amongst others. The Pf genomic resource database (PlasmoDB) identified 41 missense mutations located in the partial zymogen and catalytic domains of FP-2 and FP-3. Using structure-based tools, six putative allosteric pockets were identified in FP-2 and FP-3. The effect of mutations on the whole protein, the central core, binding pocket residues and allosteric pockets was evaluated. The accurate 3D homology models of the WT and mutants were calculated. MD simulations were performed on the various systems as a quick starting point. MD simulations have provided a cornerstone for establishing numerous computational tools for describing changes arising from mutations, ligand binding, and environmental changes such as pH and temperature. Post-MD analysis was performed in two stages viz global and local analysis. Global analysis via radius of gyration (Rg) and comparative essential dynamic analysis revealed the conformational variability associated with all mutations. In the catalytic domain of FP-2, the presence of M245I mutation triggered the formation of a cryptic pocket via an exclusive mechanism involving the fusion of pockets 2 and 6. This striking observation was also detected in the partial zymogen complex of FP-2 and induced by A159V, M245I and E249A mutations. A similar observation was uncovered in the presence of A422T mutation in the catalytic domain of FP-3. Local DRN and contact map analyses identified conserved inter-residue interaction changes on important communication networks. This study brings a novel understanding of the effects of missense mutations in FP-2 and FP-3 and provides important insight which may help discover new anti-hemoglobinase drugs. The second aim is the identification of potential allosteric ligands against the WT and mutant systems of FP-2 and FP-3 using various computational tools. Of the six potential allosteric pockets identified in FP-2 and FP-3, pocket 1 was evaluated by SiteMap as the most druggable in both proteins. This pipeline was implemented to screen pocket 1 of FP-2 and FP-3 against 2089 repositionable compounds obtained from the DrugBank database. In order to ensure selectivity and specificity to the Plasmodium protein, the human homologs (Cat K and Cat L) were screened, and compounds binding to these proteins were exempted from further analysis. Subsequently, eight compounds (DB00128, DB00312, DB00766, DB00951, DB02893, DB03754, DB13972, and DB14159) were identified as potential allosteric hits for FP-2 and five (DB00853, DB00951, DB01613, DB04173 and DB09419) for FP-3. These compounds were subjected to MD simulation and post-MD trajectory analysis to ascertain their stability in their respective protein structures. The effects of the stable compounds on the WT and mutant systems of FP-2 and FP-3 were then evaluated using DRN analysis. Attention has recently been drawn towards identifying novel allosteric compounds targeting FP-2 and FP-3; hence this study explores the potential allosteric inhibitory mechanisms in the presence and absence of mutations in FP-2 and FP-3. Overall, the results presented in this thesis provide (i) an understanding of the role mutations in the partial zymogen complex play in the activation of the active enzyme, (ii) an insight into the possible allosteric mechanisms induced by mutations on the active enzymes, and (iii) a computational pipeline for the development of novel allosteric modulators for malaria inhibition studies.