Identification of novel compounds against Plasmodium falciparum Cytochrome bc1 Complex inhibiting the trans-membrane electron transfer pathway: an In Silico study
- Authors: Chebon, Lorna Jemosop
- Date: 2022-10-14
- Subjects: Malaria , Plasmodium falciparum , Molecular dynamics , Antimalarials , Molecules Models , Docking , Cytochromes , Drug resistance , Computer simulation , Drugs Computer-aided design , System analysis
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/365666 , vital:65774 , DOI https://doi.org/10.21504/10962/365666
- Description: Malaria continues to be a burden globally with a myriad of challenges deterring eradication efforts. With most antimalarials facing drug resistance, such as atovaquone (ATQ), alternative compounds that can withstand resistance are warranted. The Plasmodium falciparum cytochrome b (PfCytb), a subunit of P. falciparum cytochrome bc1 complex, is a validated drug target. Structurally, cytochrome b, cytochrome c1, and iron sulphur protein (ISP) subunits form the catalytic domain of the protein complex having heme bL, heme bH and iron-sulphur [2FE-2S] cluster cofactors. These cofactos have redox centres to aid in the electron transfer (ET) process. These subunits promote ET mainly through the enzyme’s ubiquinol oxidation (Qo) and ubiquinone reduction (Qi) processes in the catalytic domain. ATQ drug has been used in the prevention and treatment of uncomplicated malaria by targeting PfCytb protein. Once the mitochondrial transmembrane ET pathway is inhibited, it causes a collapse in its membrane potential. Previously reported ATQ drug resistance has been associated with the point mutations Y268C, Y268N and Y268S. Thus, in finding alternatives to the ATQ drug, this research aimed to: i) employ in silico approaches incorporating protein into phospholipid bilayer for the first time to understand the parasites’ resistance mechanism; ii) determine any sequence and structural differences that could be explored in drug design studies; and iii) screen for PfCytb-iron sulphur protein (Cytb-ISP) hit compounds from South African natural compound database (SANCDB) and Medicines for Malaria Venture (MMV) that can withstand the identified mutations. Using computational tools, comparative sequence and structural analyses were performed on the cytochrome b protein, where the ultimate focus was on P. falciparum cytochrome b and its human homolog. Through multiple sequence alignment, motif discovery and phylogeny, differences between P. falciparum and H. sapiens cytochrome b were identified. Protein modelling of both P. falciparum and H. sapiens cytochrome b - iron sulphur protein (PfCytb-ISP and HsCytb-ISP) was performed. Results showed that at the sequence level, there were few amino acid residue differences because the protein is highly conserved. Important to note is the four-residue deletion in Plasmodium spp. absent in the human homolog. Motif analysis discovered five unique motifs in P. falciparum cytochrome b protein which were mapped onto the predicted protein model. These motifs were not in regions of functional importance; hence their function is still unknown. At a structural level, the four-residue deletion was observed to alter the Qo substrate binding pocket as reported in previous studies and confirmed in this study. This deletion resulted in a 0.83 Å structural displacement. Also, there are currently no in silico studies that have performed experiments with P. falciparum cytochrome b protein incorporated into a phospholipid bilayer. Using 350 ns molecular dynamics (MD) simulations of the holo and ATQ-bound systems, the study highlighted the resistance mechanism of the parasite protein where the loss of active site residue-residue interactions was identified, all linked to the three mutations. The identified compromised interactions are likely to destabilise the protein’s function, specifically in the Qo substrate binding site. This showed the possible effect of mutations on ATQ drug activity, where all three mutations were reported to share a similar resistance mechanism. Thereafter, this research work utilised in silico approaches where both Qo active site and interface pocket were targeted by screening the South African natural compounds database (SANCDB) and Medicines for Malaria Venture (MMV) compounds to identify novel selective hits. SANCDB compounds are known for their structural complexity that preserves the potency of the drug molecule. Both SANCDB and MMV compounds have not been explored as inhibitors against the PfCytb drug target. Molecular docking, molecular dynamics (MD) simulations, principal component, and dynamic residue network (DRN; global and local) analyses were utilised to identify and confirm the potential selective inhibitors. Docking results identified compounds that bound selectively onto PfCytb-ISP with a binding energy ≤ -8.7 kcal/mol-1. Further, this work validated a total of eight potential selective compounds to inhibit PfCytb-ISP protein (Qo active site) not only in the wild-type but also in the presence of the point mutations Y268C, Y268N and Y268S. The selective binding of these hit compounds could be linked to the differences reported at sequence/residue level in chapter 3. DRN and residue contact map analyses of the eight compounds in holo and ligand-bound systems revealed reduced residue interactions and decreased protein communication. This suggests that the eight compounds show the possibility of inhibiting the parasite and disrupting important residue-residue interactions. Additionally, 13 selective compounds were identified to bind at the protein’s heterodimer interface, where global and local analysis confirmed their effect on active site residues (distal location) as well as on the communication network. Based on the sequence differences between PfCytb and the human homolog, these findings suggest these selective compounds as potential allosteric modulators of the parasite enzyme, which may serve as possible replacements of the already resistant ATQ drug. Therefore, these findings pave the way for further in vitro studies to establish their anti-plasmodial inhibition levels. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-10-14
- Authors: Chebon, Lorna Jemosop
- Date: 2022-10-14
- Subjects: Malaria , Plasmodium falciparum , Molecular dynamics , Antimalarials , Molecules Models , Docking , Cytochromes , Drug resistance , Computer simulation , Drugs Computer-aided design , System analysis
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/365666 , vital:65774 , DOI https://doi.org/10.21504/10962/365666
- Description: Malaria continues to be a burden globally with a myriad of challenges deterring eradication efforts. With most antimalarials facing drug resistance, such as atovaquone (ATQ), alternative compounds that can withstand resistance are warranted. The Plasmodium falciparum cytochrome b (PfCytb), a subunit of P. falciparum cytochrome bc1 complex, is a validated drug target. Structurally, cytochrome b, cytochrome c1, and iron sulphur protein (ISP) subunits form the catalytic domain of the protein complex having heme bL, heme bH and iron-sulphur [2FE-2S] cluster cofactors. These cofactos have redox centres to aid in the electron transfer (ET) process. These subunits promote ET mainly through the enzyme’s ubiquinol oxidation (Qo) and ubiquinone reduction (Qi) processes in the catalytic domain. ATQ drug has been used in the prevention and treatment of uncomplicated malaria by targeting PfCytb protein. Once the mitochondrial transmembrane ET pathway is inhibited, it causes a collapse in its membrane potential. Previously reported ATQ drug resistance has been associated with the point mutations Y268C, Y268N and Y268S. Thus, in finding alternatives to the ATQ drug, this research aimed to: i) employ in silico approaches incorporating protein into phospholipid bilayer for the first time to understand the parasites’ resistance mechanism; ii) determine any sequence and structural differences that could be explored in drug design studies; and iii) screen for PfCytb-iron sulphur protein (Cytb-ISP) hit compounds from South African natural compound database (SANCDB) and Medicines for Malaria Venture (MMV) that can withstand the identified mutations. Using computational tools, comparative sequence and structural analyses were performed on the cytochrome b protein, where the ultimate focus was on P. falciparum cytochrome b and its human homolog. Through multiple sequence alignment, motif discovery and phylogeny, differences between P. falciparum and H. sapiens cytochrome b were identified. Protein modelling of both P. falciparum and H. sapiens cytochrome b - iron sulphur protein (PfCytb-ISP and HsCytb-ISP) was performed. Results showed that at the sequence level, there were few amino acid residue differences because the protein is highly conserved. Important to note is the four-residue deletion in Plasmodium spp. absent in the human homolog. Motif analysis discovered five unique motifs in P. falciparum cytochrome b protein which were mapped onto the predicted protein model. These motifs were not in regions of functional importance; hence their function is still unknown. At a structural level, the four-residue deletion was observed to alter the Qo substrate binding pocket as reported in previous studies and confirmed in this study. This deletion resulted in a 0.83 Å structural displacement. Also, there are currently no in silico studies that have performed experiments with P. falciparum cytochrome b protein incorporated into a phospholipid bilayer. Using 350 ns molecular dynamics (MD) simulations of the holo and ATQ-bound systems, the study highlighted the resistance mechanism of the parasite protein where the loss of active site residue-residue interactions was identified, all linked to the three mutations. The identified compromised interactions are likely to destabilise the protein’s function, specifically in the Qo substrate binding site. This showed the possible effect of mutations on ATQ drug activity, where all three mutations were reported to share a similar resistance mechanism. Thereafter, this research work utilised in silico approaches where both Qo active site and interface pocket were targeted by screening the South African natural compounds database (SANCDB) and Medicines for Malaria Venture (MMV) compounds to identify novel selective hits. SANCDB compounds are known for their structural complexity that preserves the potency of the drug molecule. Both SANCDB and MMV compounds have not been explored as inhibitors against the PfCytb drug target. Molecular docking, molecular dynamics (MD) simulations, principal component, and dynamic residue network (DRN; global and local) analyses were utilised to identify and confirm the potential selective inhibitors. Docking results identified compounds that bound selectively onto PfCytb-ISP with a binding energy ≤ -8.7 kcal/mol-1. Further, this work validated a total of eight potential selective compounds to inhibit PfCytb-ISP protein (Qo active site) not only in the wild-type but also in the presence of the point mutations Y268C, Y268N and Y268S. The selective binding of these hit compounds could be linked to the differences reported at sequence/residue level in chapter 3. DRN and residue contact map analyses of the eight compounds in holo and ligand-bound systems revealed reduced residue interactions and decreased protein communication. This suggests that the eight compounds show the possibility of inhibiting the parasite and disrupting important residue-residue interactions. Additionally, 13 selective compounds were identified to bind at the protein’s heterodimer interface, where global and local analysis confirmed their effect on active site residues (distal location) as well as on the communication network. Based on the sequence differences between PfCytb and the human homolog, these findings suggest these selective compounds as potential allosteric modulators of the parasite enzyme, which may serve as possible replacements of the already resistant ATQ drug. Therefore, these findings pave the way for further in vitro studies to establish their anti-plasmodial inhibition levels. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-10-14
The characterization of GTP Cyclohydrolase I and 6-Pyruvoyl Tetrahydropterin Synthase enzymes as potential anti-malarial drug targets
- Khairallah, Afrah Yousif Huseein
- Authors: Khairallah, Afrah Yousif Huseein
- Date: 2022-04-08
- Subjects: Antimalarials , Plasmodium falciparum , Malaria Chemotherapy , Malaria Africa , Drug resistance , Drug development , Molecular dynamics
- Language: English
- Type: Doctoral thesis , text
- Identifier: http://hdl.handle.net/10962/233784 , vital:50127 , DOI 10.21504/10962/233784
- Description: Malaria remains a public health problem and a high burden of disease, especially in developing countries. The unicellular protozoan malaria parasite of the genus Plasmodium infects about a quarter of a billion people annually, with an estimated 409 000 death cases. The majority of malaria cases occurred in Africa; hence, the region is regarded as endemic for malaria. Global efforts to eradicate the disease led to a decrease in morbidity and mortality rates. However, an enormous burden of malaria infection remains, and it cannot go unnoticed. Countries with limited resources are more affected by the disease, mainly on its public health and socio-economic development, due to many factors besides malaria itself, such as lack of access to adequate, affordable treatments and preventative regimes. Furthermore, the current antimalarial drugs are losing their efficacy because of parasite drug resistance. The emerged drug resistance has reduced the drug efficacy in clearing the parasite from the host system, causing prolonged illness and a higher risk of death. Therefore, the emerged antimalarial drug resistance has hindered the global efforts for malaria control and elimination and established an urgent need for new treatment strategies. When the resistance against classical antimalarial drugs emerged, the class of antifolate antimalarial medicines became the most common alternative. The antifolate antimalarial drugs target the malaria parasite de novo folate biosynthesis pathway by limiting folate derivates, which are essential for the parasite cell growth and survival. Yet again, the malaria parasite developed resistance against the available antifolate drugs, rendering the drugs ineffective in many cases. Given the previous success in targeting the malaria parasite de novo folate biosynthesis pathway, alternative enzymes within this pathway stand as good targets and can be explored to develop new antifolate drugs with novel mechanisms of action. The primary focus of this thesis is to contribute to the existing and growing knowledge of antimalarial drug discovery. The study aims to characterise the malaria parasite de novo folate synthesis pathway enzymes guanosine-5'-triphosphate (GTP) cyclohydrolase I (GCH1) and 6-pyruvoyl tetrahydropterin synthase (PTPS) as alternative drug targets for malaria treatment by using computational approaches. Further, discover new allosteric drug targeting sites within the two enzymes' 3D structures for future drug design and discovery. Sequence and structural analysis were carried out to characterise and pinpoint the two enzymes' unique sequence and structure-based features. From the analyses, key sequence and structure differences were identified between the malaria parasite enzymes relative to their human homolog; the identified sites can aid significantly in designing and developing new antimalarial antifolate drugs with good selectivity toward the parasites’ enzymes. GCH1 and PTPS contain a catalytically essential metal ion in their active site; therefore, force field parameters were needed to study their active sites accurately during all-atom molecular dynamic simulations (MD). The force field parameters were derived through quantum mechanics potential energy surface scans of the metals bonded terms and evaluated via all-atom MD simulations. Proteins structural dynamics is imperative for many biological processes; thus, it is essential to consider the structural dynamics of proteins whilst understanding their function. In this regard, the normal mode analysis (NMA) approach based on the elastic network model (ENM) was employed to study the intrinsic dynamics and conformations changes of GCH1 and PTPS enzymes. The NMA disclosed essential structural information about the protein’s intrinsic dynamics and mechanism of allosteric modulation of their binding properties, further highlighting regions that govern their conformational changes. The analysis also disclosed hotspot residues that are crucial for the proteins' fold stability and function. The NMA was further combined with sequence motif results and showed that conserved residues of GCH1 and PTPS were located within the identified key structural sites modulating the proteins' conformational rearrangement. The characterized structural features and hotspot residues were regarded as potential allosteric sites of important value for the design and development of allosteric drugs. Both GCH1 and PTPS enzymes have never been targeted before and can provide an excellent opportunity to overcome the antimalarial antifolate drug resistance problem. The data presented in this thesis contribute to the understanding of the sequence, structure, and global dynamics of both GCH1 and PTPS, further disclose potential allosteric drug targeting sites and unique structural features of both enzymes that can establish a solid starting point for drug design and development of new antimalarial drugs of a novel mechanism of actions. Lastly, the reported force field parameters will be of value for MD simulations for future in-silico drug discovery studies involving the two enzymes and other enzymes with the same Zn2+ binding motifs and coordination environments. The impact of this research can facilitate the discovery of new effective antimalarial medicines with novel mechanisms of action. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-04-08
- Authors: Khairallah, Afrah Yousif Huseein
- Date: 2022-04-08
- Subjects: Antimalarials , Plasmodium falciparum , Malaria Chemotherapy , Malaria Africa , Drug resistance , Drug development , Molecular dynamics
- Language: English
- Type: Doctoral thesis , text
- Identifier: http://hdl.handle.net/10962/233784 , vital:50127 , DOI 10.21504/10962/233784
- Description: Malaria remains a public health problem and a high burden of disease, especially in developing countries. The unicellular protozoan malaria parasite of the genus Plasmodium infects about a quarter of a billion people annually, with an estimated 409 000 death cases. The majority of malaria cases occurred in Africa; hence, the region is regarded as endemic for malaria. Global efforts to eradicate the disease led to a decrease in morbidity and mortality rates. However, an enormous burden of malaria infection remains, and it cannot go unnoticed. Countries with limited resources are more affected by the disease, mainly on its public health and socio-economic development, due to many factors besides malaria itself, such as lack of access to adequate, affordable treatments and preventative regimes. Furthermore, the current antimalarial drugs are losing their efficacy because of parasite drug resistance. The emerged drug resistance has reduced the drug efficacy in clearing the parasite from the host system, causing prolonged illness and a higher risk of death. Therefore, the emerged antimalarial drug resistance has hindered the global efforts for malaria control and elimination and established an urgent need for new treatment strategies. When the resistance against classical antimalarial drugs emerged, the class of antifolate antimalarial medicines became the most common alternative. The antifolate antimalarial drugs target the malaria parasite de novo folate biosynthesis pathway by limiting folate derivates, which are essential for the parasite cell growth and survival. Yet again, the malaria parasite developed resistance against the available antifolate drugs, rendering the drugs ineffective in many cases. Given the previous success in targeting the malaria parasite de novo folate biosynthesis pathway, alternative enzymes within this pathway stand as good targets and can be explored to develop new antifolate drugs with novel mechanisms of action. The primary focus of this thesis is to contribute to the existing and growing knowledge of antimalarial drug discovery. The study aims to characterise the malaria parasite de novo folate synthesis pathway enzymes guanosine-5'-triphosphate (GTP) cyclohydrolase I (GCH1) and 6-pyruvoyl tetrahydropterin synthase (PTPS) as alternative drug targets for malaria treatment by using computational approaches. Further, discover new allosteric drug targeting sites within the two enzymes' 3D structures for future drug design and discovery. Sequence and structural analysis were carried out to characterise and pinpoint the two enzymes' unique sequence and structure-based features. From the analyses, key sequence and structure differences were identified between the malaria parasite enzymes relative to their human homolog; the identified sites can aid significantly in designing and developing new antimalarial antifolate drugs with good selectivity toward the parasites’ enzymes. GCH1 and PTPS contain a catalytically essential metal ion in their active site; therefore, force field parameters were needed to study their active sites accurately during all-atom molecular dynamic simulations (MD). The force field parameters were derived through quantum mechanics potential energy surface scans of the metals bonded terms and evaluated via all-atom MD simulations. Proteins structural dynamics is imperative for many biological processes; thus, it is essential to consider the structural dynamics of proteins whilst understanding their function. In this regard, the normal mode analysis (NMA) approach based on the elastic network model (ENM) was employed to study the intrinsic dynamics and conformations changes of GCH1 and PTPS enzymes. The NMA disclosed essential structural information about the protein’s intrinsic dynamics and mechanism of allosteric modulation of their binding properties, further highlighting regions that govern their conformational changes. The analysis also disclosed hotspot residues that are crucial for the proteins' fold stability and function. The NMA was further combined with sequence motif results and showed that conserved residues of GCH1 and PTPS were located within the identified key structural sites modulating the proteins' conformational rearrangement. The characterized structural features and hotspot residues were regarded as potential allosteric sites of important value for the design and development of allosteric drugs. Both GCH1 and PTPS enzymes have never been targeted before and can provide an excellent opportunity to overcome the antimalarial antifolate drug resistance problem. The data presented in this thesis contribute to the understanding of the sequence, structure, and global dynamics of both GCH1 and PTPS, further disclose potential allosteric drug targeting sites and unique structural features of both enzymes that can establish a solid starting point for drug design and development of new antimalarial drugs of a novel mechanism of actions. Lastly, the reported force field parameters will be of value for MD simulations for future in-silico drug discovery studies involving the two enzymes and other enzymes with the same Zn2+ binding motifs and coordination environments. The impact of this research can facilitate the discovery of new effective antimalarial medicines with novel mechanisms of action. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-04-08
Application of computer-aided drug design for identification of P. falciparum inhibitors
- Authors: Diallo, Bakary N’tji
- Date: 2021-10-29
- Subjects: Plasmodium falciparum , Malaria -- Chemotherapy , Molecular dynamics , Antimalarials , Cheminformatics , Drug development , Ligand binding (Biochemistry) , Plasmodium falciparum1-deoxy-D-xylulose-5-phosphate reductoisomerase (PfDXR) , South African Natural Compounds Database
- Language: English
- Type: Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/192798 , vital:45265 , 10.21504/10962/192798
- Description: Malaria is a millennia-old disease with the first recorded cases dating back to 2700 BC found in Chinese medical records, and later in other civilizations. It has claimed human lives to such an extent that there are a notable associated socio-economic consequences. Currently, according to the World Health Organization (WHO), Africa holds the highest disease burden with 94% of deaths and 82% of cases with P. falciparum having ~100% prevalence. Chemotherapy, such as artemisinin combination therapy, has been and continues to be the work horse in the fight against the disease, together with seasonal malaria chemoprevention and the use of insecticides. Natural products such as quinine and artemisinin are particularly important in terms of their antimalarial activity. The emphasis in current chemotherapy research is the need for time and cost-effective workflows focussed on new mechanisms of action (MoAs) covering the target candidate profiles (TCPs). Despite a decline in cases over the past decades with, countries increasingly becoming certified malaria free, a stalling trend has been observed in the past five years resulting in missing the 2020 Global Technical Strategy (GTS) milestones. With no effective vaccine, a reduction in funding, slower drug approval than resistance emergence from resistant and invasive vectors, and threats in diagnosis with the pfhrp2/3 gene deletion, malaria remains a major health concern. Motivated by these reasons, the primary aim of this work was a contribution to the antimalarial pipeline through in silico approaches focusing on P. falciparum. We first intended an exploration of malarial targets through a proteome scale screening on 36 targets using multiple metrics to account for the multi-objective nature of drug discovery. The continuous growth of structural data offers the ideal scenario for mining new MoAs covering antimalarials TCPs. This was combined with a repurposing strategy using a set of orally available FDA approved drugs. Further, use was made of time- and cost-effective strategies combining QVina-W efficiency metrics that integrate molecular properties, GRIM rescoring for molecular interactions and a hydrogen mass repartitioning (HMR) molecular dynamics (MD) scheme for accelerated development of antimalarials in the context of resistance. This pipeline further integrates a complex ranking for better drug-target selectivity, and normalization strategies to overcome docking scoring function bias. The different metrics, ranking, normalization strategies and their combinations were first assessed using their mean ranking error (MRE). A version combining all metrics was used to select 36 unique protein-ligand complexes, assessed in MD, with the final retention of 25. From the 16 in vitro tested hits of the 25, fingolimod, abiraterone, prazosin, and terazosin showed antiplasmodial activity with IC50 2.21, 3.37, 16.67 and 34.72 μM respectively and of these, only fingolimod was found to be not safe with respect to human cell viability. These compounds were predicted active on different molecular targets, abiraterone was predicted to interact with a putative liver-stage essential target, hence promising as a transmission-blocking agent. The pipeline had a promising 25% hit rate considering the proteome-scale and use of cost-effective approaches. Secondly, we focused on Plasmodium falciparum 1-deoxy-D-xylulose-5-phosphate reductoisomerase (PfDXR) using a more extensive screening pipeline to overcome some of the current in silico screening limitations. Starting from the ZINC lead-like library of ~3M, hierarchical ligand-based virtual screening (LBVS) and structure-based virtual screening (SBVS) approaches with molecular docking and re-scoring using eleven scoring functions (SFs) were used. Later ranking with an exponential consensus strategy was included. Selected hits were further assessed through Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA), advanced MD sampling in a ligand pulling simulations and (Weighted Histogram Analysis Method) WHAM analysis for umbrella sampling (US) to derive binding free energies. Four leads had better predicted affinities in US than LC5, a 280 nM potent PfDXR inhibitor with ZINC000050633276 showing a promising binding of -20.43 kcal/mol. As shown with fosmidomycin, DXR inhibition offers fast acting compounds fulfilling antimalarials TCP1. Yet, fosmidomycin has a high polarity causing its short half-life and hampering its clinical use. These leads scaffolds are different from fosmidomycin and hence may offer better pharmacokinetic and pharmacodynamic properties and may also be promising for lead optimization. A combined analysis of residues’ contributions to the free energy of binding in MM-PBSA and to steered molecular dynamics (SMD) Fmax indicated GLU233, CYS268, SER270, TRP296, and HIS341 as exploitable for compound optimization. Finally, we updated the SANCDB library with new NPs and their commercially available analogs as a solution to NP availability. The library is extended to 1005 compounds from its initial 600 compounds and the database is integrated to Mcule and Molport APIs for analogs automatic update. The new set may contribute to virtual screening and to antimalarials as the most effective ones have NP origin. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2021
- Full Text:
- Date Issued: 2021-10-29
- Authors: Diallo, Bakary N’tji
- Date: 2021-10-29
- Subjects: Plasmodium falciparum , Malaria -- Chemotherapy , Molecular dynamics , Antimalarials , Cheminformatics , Drug development , Ligand binding (Biochemistry) , Plasmodium falciparum1-deoxy-D-xylulose-5-phosphate reductoisomerase (PfDXR) , South African Natural Compounds Database
- Language: English
- Type: Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/192798 , vital:45265 , 10.21504/10962/192798
- Description: Malaria is a millennia-old disease with the first recorded cases dating back to 2700 BC found in Chinese medical records, and later in other civilizations. It has claimed human lives to such an extent that there are a notable associated socio-economic consequences. Currently, according to the World Health Organization (WHO), Africa holds the highest disease burden with 94% of deaths and 82% of cases with P. falciparum having ~100% prevalence. Chemotherapy, such as artemisinin combination therapy, has been and continues to be the work horse in the fight against the disease, together with seasonal malaria chemoprevention and the use of insecticides. Natural products such as quinine and artemisinin are particularly important in terms of their antimalarial activity. The emphasis in current chemotherapy research is the need for time and cost-effective workflows focussed on new mechanisms of action (MoAs) covering the target candidate profiles (TCPs). Despite a decline in cases over the past decades with, countries increasingly becoming certified malaria free, a stalling trend has been observed in the past five years resulting in missing the 2020 Global Technical Strategy (GTS) milestones. With no effective vaccine, a reduction in funding, slower drug approval than resistance emergence from resistant and invasive vectors, and threats in diagnosis with the pfhrp2/3 gene deletion, malaria remains a major health concern. Motivated by these reasons, the primary aim of this work was a contribution to the antimalarial pipeline through in silico approaches focusing on P. falciparum. We first intended an exploration of malarial targets through a proteome scale screening on 36 targets using multiple metrics to account for the multi-objective nature of drug discovery. The continuous growth of structural data offers the ideal scenario for mining new MoAs covering antimalarials TCPs. This was combined with a repurposing strategy using a set of orally available FDA approved drugs. Further, use was made of time- and cost-effective strategies combining QVina-W efficiency metrics that integrate molecular properties, GRIM rescoring for molecular interactions and a hydrogen mass repartitioning (HMR) molecular dynamics (MD) scheme for accelerated development of antimalarials in the context of resistance. This pipeline further integrates a complex ranking for better drug-target selectivity, and normalization strategies to overcome docking scoring function bias. The different metrics, ranking, normalization strategies and their combinations were first assessed using their mean ranking error (MRE). A version combining all metrics was used to select 36 unique protein-ligand complexes, assessed in MD, with the final retention of 25. From the 16 in vitro tested hits of the 25, fingolimod, abiraterone, prazosin, and terazosin showed antiplasmodial activity with IC50 2.21, 3.37, 16.67 and 34.72 μM respectively and of these, only fingolimod was found to be not safe with respect to human cell viability. These compounds were predicted active on different molecular targets, abiraterone was predicted to interact with a putative liver-stage essential target, hence promising as a transmission-blocking agent. The pipeline had a promising 25% hit rate considering the proteome-scale and use of cost-effective approaches. Secondly, we focused on Plasmodium falciparum 1-deoxy-D-xylulose-5-phosphate reductoisomerase (PfDXR) using a more extensive screening pipeline to overcome some of the current in silico screening limitations. Starting from the ZINC lead-like library of ~3M, hierarchical ligand-based virtual screening (LBVS) and structure-based virtual screening (SBVS) approaches with molecular docking and re-scoring using eleven scoring functions (SFs) were used. Later ranking with an exponential consensus strategy was included. Selected hits were further assessed through Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA), advanced MD sampling in a ligand pulling simulations and (Weighted Histogram Analysis Method) WHAM analysis for umbrella sampling (US) to derive binding free energies. Four leads had better predicted affinities in US than LC5, a 280 nM potent PfDXR inhibitor with ZINC000050633276 showing a promising binding of -20.43 kcal/mol. As shown with fosmidomycin, DXR inhibition offers fast acting compounds fulfilling antimalarials TCP1. Yet, fosmidomycin has a high polarity causing its short half-life and hampering its clinical use. These leads scaffolds are different from fosmidomycin and hence may offer better pharmacokinetic and pharmacodynamic properties and may also be promising for lead optimization. A combined analysis of residues’ contributions to the free energy of binding in MM-PBSA and to steered molecular dynamics (SMD) Fmax indicated GLU233, CYS268, SER270, TRP296, and HIS341 as exploitable for compound optimization. Finally, we updated the SANCDB library with new NPs and their commercially available analogs as a solution to NP availability. The library is extended to 1005 compounds from its initial 600 compounds and the database is integrated to Mcule and Molport APIs for analogs automatic update. The new set may contribute to virtual screening and to antimalarials as the most effective ones have NP origin. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2021
- Full Text:
- Date Issued: 2021-10-29
- «
- ‹
- 1
- ›
- »