Computational studies in human African trypanosomiasis
- Authors: Muronzi, Tendai
- Date: 2023-10-13
- Subjects: African trypanosomiasis , Apolipoprotein L1 , Docking , Protein-protein interactions , Homology modeling , Tetrahydrofolate dehydrogenase , Pteridine reductase
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/431883 , vital:72812 , DOI 10.21504/10962/431885
- Description: Human African trypanosomiasis (HAT) is a neglected tropical disease (NTD) caused by two subspecies of the parasite, namely Trypanosoma brucei (Tb) gambiense (g-HAT) and rhodesiense (r-HAT). HAT is endemic in sub-Saharan countries, where the parasite transmission vectors, tsetse flies, breed. An estimated 70 million people remain at risk of contracting the disease, where the infection is classified as acute or chronic for g-HAT and r-HAT, respectively, with both forms ending in fatal meningoencephalitis when left untreated. Both g-HAT and r-HAT are responsible for widespread fatal epidemics throughout sub-Saharan African history, resulting from the complex molecular interplay between trypanosomes and humans through unique, innate immunity evasion mechanisms. Of interest, the Tbr subspecies expresses a serum resistance-associated protein (SRA), which binds to human serum lytic factor, apolipoprotein L1 (ApoL1), nullifying any trypanocidal activity. In response, ApoL1 (G1 and G2) variants found in humans of sub-Saharan African lineage have been cited for conferring resistance to the r-HAT infection in an interaction that is not fully elucidated In the event of successful infection, current HAT chemotherapeutics are plagued with complexity of administration, poor efficacy, toxicity, and potential drug resistance, highlighting a need for improved approaches. The parasite folate pathway provides a strategic target for alternative anti-trypanosomal drug development as trypanosomatids are folate auxotrophs, requiring host folate for growth and survival. Validated drug targets pteridine reductase (TbPTR1) and dihydrofolate reductase (TbDHFR) are essential for salvaging cofactors folate and folate biopterin crucial to parasite survival, making them viable targets for anti-folate investigation. The overall aims of this thesis were to a) provide insights into the molecular and dynamic basis of the SRA and ApoL1 interplay in HAT infection and b) identify safer and more efficient anti-folate anti-trypanosomal drug alternatives through in silico approaches. To achieve our first aim, in silico structure prediction was applied to generate 3D models of ApoL1 C-terminal variants G0, G1, G1G/M, G2 and G1G2, and four SRA variants retrieved from the NCBI database. The SRA and ApoL1 structures were inspected dynamically to identify the effect of the variants through molecular dynamics (MD) simulations. Dynamic residue network (DRN) analysis of MD trajectories was fundamental in identifying residues playing a vital role in the intramolecular communication of both proteins in the presence of mutations. Protein-protein docking was then applied to calculate plausible SRA-ApoL1 C-terminal wild-type complex structures to further elucidate the nature of SRA-mediated infection. Through MD simulations, twelve SRA-ApoL1 dimeric structures were narrowed down from five to two energetically sound complexes. The two feasible SRA-ApoL1 complexes (1 and 2) exhibited favourable communication observed through DRN analysis, including the retaining key communication residues identified in prior monomer DRN calculations. ApoL1 C-terminal variants were additionally incorporated into SRA-ApoL1 complexes 1 and 2 for further complex dynamics analysis This investigation into the nature of SRA-ApoL1 binding resulted in five primary outcomes: 1) highlighting the intramolecular effects ApoL1 variants have on the stability of the protein, 2) the identification of crucial SRA and ApoL1 communication residues in both monomeric or dimeric form, 3) the isolation of feasible SRA-ApoL1 complexes determined through global and local structural analyses 4) identification of residues crucial to the complex formation and maintenance of SRA-ApoL1, overlapping with those identified in (1), and 5) the minimal dissociative role of the G1 mutations in the complex, but compounding effect of the G2 deletion mutation. Computational modelling and drug repurposing were employed to achieve the thesis's second aim as they drastically cut down the costs involved in drug discovery and provide a more time-efficient screening method through numerous drug candidates. Using high throughput virtual screening, a subset of 2089 approved DrugBank compounds were screened against TbPTR1. The outputs were filtered to 24 viable compounds in 54 binding poses using binding energy and molecular interactions. Through subsequent MD simulations of 200ns, thirteen potential hit compounds were identified. The resultant hit compounds were subjected to further blind docking against TbDHFR and molecular dynamics to identify compounds with the potential for dual inhibition. The filtered subset was also tested in in vitro single concentration and dose-response bioassays to assess inhibitory properties against Trypanosoma brucei, complementing in silico findings. Post-molecular dynamics, four compounds exhibited high stabilities and molecular interactions with both TbPTR1 and TbDHFR, with two presenting favourable results in the in vitro assays. Three compounds additionally shared common structural moieties. In all, the in silico repurposing highlighted drugs characterised by favourable interactions and stabilities in TbPTR1, thus providing (1) a framework for further studies investigating anti-folate HAT compounds and (2) modulatory scaffolds based on identified moieties that can be used for the design of safe anti-folate trypanosomal drugs. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2023
- Full Text:
- Date Issued: 2023-10-13
An in-silico study of the type II NADH: Quinone Oxidoreductase (ndh2). A new anti-malaria drug target
- Authors: Baye, Bertha Cinthia
- Date: 2022-10-14
- Subjects: Malaria , Plasmodium , Molecular dynamics , Computer simulation , Quinone , Antimalarials , Molecules Models , Docking , Drugs Computer-aided design
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/365633 , vital:65767 , DOI https://doi.org/10.21504/10962/365633
- Description: Malaria is caused by Plasmodium parasites, spread to people through the bites of infected female Anopheles mosquitoes. This study focuses on all 5 (Plasmodium falciparum, Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale and Plasmodium vivax) parasites that cause malaria in humans. Africa is a developing continent, and it is the most affected with an estimation of 90% of more than 400 000 malaria-related deaths reported by the World Health Organization (WHO) report in 2020, in which 61% of that number are children under the ages of five. Malaria resistance was initially observed in early 1986 and with the progression of time anti-malarial drug resistance has only increased. As a result, there is a need to study the malarial proteins mechanism of action and identify alternative treatment strategies for this disease. Type II NADH: quinone oxidoreductase (NDH2) is a monotopic protein that catalyses the electron transfer from NADH to quinone via FAD without a proton-pumping activity, and functions as an initial enzyme, either in addition to or as an alternative to proton-pumping NADH dehydrogenase (complex I) in the respiratory chain of bacteria, archaea, and fungal and plant mitochondrial. The structures for the Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale and Plasmodium vivax were modelled from the crystal structure of Plasmodium falciparum (5JWA). Compounds from the South African natural compounds database (SANCDB) were docked against both the NDH2 crystal structure and modelled structures. By performing in silico screening the study aimed to find potential compounds that might interrupt the electron transfer to quinone therefore disturbing the enzyme‟s function and thereby possibly eliminating the plasmodium parasite. CHARMM-GUI was used to create the membrane (since this work is with membrane-bound proteins) and to orient the protein on the membrane using OPM server guidelines, the interface produced GROMACS topology files that were used in molecular dynamics simulations. Molecular dynamics simulations were performed in the Centre for high performance computing (CHPC) cluster under the CHEM0802 project and the trajectories produced were further analysed. In this work not only were hit compounds from SANCDB identified, but also differences in behaviour across species and in the presence or absence of the membrane were described. This highlights the need to include the correct protein environment when studying these systems. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-10-14
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
Identification of selective novel hits against Mycobacterium tuberculosis KasA potential allosteric sites using bioinformatics approaches
- Authors: Hare, Fadzayi Faith
- Date: 2022-10-14
- Subjects: Tuberculosis , Docking , Molecules Models , Virtual screening , Multidrug-resistant tuberculosis , Fatty acids Synthesis , Drugs Design
- Language: English
- Type: Academic theses , Master's theses , text
- Identifier: http://hdl.handle.net/10962/362842 , vital:65367
- Description: Tuberculosis (TB) is a global health threat that has led to approximately 1.5 million deaths annually. According to the World Health Organization (WHO), TB is among the top ten deadly diseases and is the leading cause of death due to a single infectious agent. The main challenge in the effective treatment and control of TB is the ongoing emergence of resistant strains of Mycobacterium tuberculosis (Mtb) which lead to multi-drug resistant (MDR) and extensive-drug resistant (XDR) TB. Hence, the identification and characterization of novel drug targets and drugs that modulate the activity of the pathogen are an urgent priority. The current situation even necessitates the reengineering or repurposing of drugs in order to achieve effective control. The β-ketoacyl-acyl carrier protein synthase I (KasA) of Mycobacterium tuberculosis is an essential enzyme in the mycobacterial fatty acid synthesis (FAS-II) pathway and is believed to be a promising target for drug discovery in TB. It is one of the five main proteins of the FAS-II pathway and catalyzes a key condensation reaction in the synthesis of meromycolate chains, the precursors of mycolic acids involved in cell wall formation. Although this protein has been extensively studied, little research has been devoted to the allosteric inhibition of potential drug compounds. The main aim of this research was to identify the allosteric sites on the protein that could be involved in the inhibition of substrate binding activities and novel drug compounds that bind to these sites by use of in-silico approaches. The bioinformatics approaches used in this study were divided into four main objectives namely identification of KasA homolog sequences, sequence analysis and protein characterization, allosteric site search and lastly virtual screening of DrugBank compounds via molecular docking. Fifteen homolog sequences were identified from the BLASTP analysis and were derived from bacteria, fungi and mammals. In order to discover important residues and regions within the KasA proteins, sequence alignment, motif analysis and phylogenetic studies were performed using Mtb KasA as a reference. Sequence alignment revealed conserved residues in all KasA proteins that have functional importance such as the catalytic triad residues (Cys171, His311 and His345). Motif analysis identified 18 highly conserved motifs within the KasA proteins with structural and functional roles. In addition, motifs unique to the Mtb KasA protein were also identified and explored for inhibitor drug design purposes. Phylogenetic analysis of the homolog sequences showed a distinct clustering of prokaryotes and eukaryotes. A distinctive clustering was also observed for species belonging to the same genus. Since the mechanism of action of most drugs involves the active site, allosteric site search was conducted on Mtb KasA and the human homolog protein using a combination of pocket detection algorithms with the aim of identifying sites that could be utilized in allosteric modulator drug discovery. This was followed by the virtual screening of 2089 FDA approved DrugBank compounds against the entire protein surfaces of Mtb KasA and Hsmt KasA, performed via molecular docking using AutoDock Vina. Screening of the compounds was based on the binding energies, with more focus on identifying ligands that bound exclusively to the acyl-binding tunnel of Mtb KasA. This reduced the data set to 27 promising drug compounds with a relatively high binding affinity for Mtb KasA, however, further experiments need to be performed to validate this result. Among these compounds were DB08889, DB06755, DB09270, DB11226, DB00392, DB12278, DB08936, DB00781, DB13720 and DB00392, which displayed relatively low binding energies for Mtb KasA when compared to the human homolog protein. , Thesis (MSc) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-10-14
Synthesis, In-Silico molecular modelling and biological studies of 1,4-Dihydroxyanthraquinone and its derivatives
- Authors: Kisula, Lydia Mboje
- Date: 2022-10-14
- Subjects: Computer simulation , Molecules Models , Dihydroxyanthraquinone , Trypanosomiasis , Leishmaniasis , Docking
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/365828 , vital:65793 , DOI https://doi.org/10.21504/10962/365828
- Description: This current study of investigation reports on the synthesis of 1,4-dihydroxyanthraquinone and its derivatives on explorations of their medicinal potential. The study initially aimed to synthesize an analogue of a natural anthraquinone, 1,3,6-trihydroxy-7-((S)-1- hydroxyethyl)anthracene-9,10-dione 5 using Friedel-Crafts acylation of phthalic anhydride and a benzene derivative. Synthetic transformation of anacardic acid 63, obtained as a by- product of the cashew industry successfully afforded 4-ethoxyisobenzofuran-1,3-dione 89. However, when attempted to couple 4-ethoxyisobenzofuran-1,3-dione 89 with 2- hydroxyacetophenone 91 in a Friedel-Crafts acylation manner to form 2-acetyl-1,8- dihydroxyanthracene-9,10-dione 87 the reaction did not work efficiently. A simple derivative of benzene which is; benzene-1,4-diol 102 was reacted instead with 3-ethoxyphthalic acid 71 and isobenzofuran-1,3-dione 96 to form 1,4,5-trihydroxy anthraquinone 72 and 1,4- dihydroxyanthraquinone 42, respectively. A modified Marschalk reaction was then used to introduce the hydroxyl alkyl group to 1,4-dihydroxy anthraquinone 42, which allowed further elaboration of the hydroxyl-substituent in moderate to good yields (22-80%). A molecular docking study was performed using Schrödinger software to predict the binding affinity of the test compounds to the target protein trypanothione reductase (PDB ID: 6BU7). An in-vitro screening of 1,4-dihydroxyanthraquinone derivatives and some selected precursors for antitrypanosomal, antiplasmodial, antibacterial, and cytotoxicity activities produced encouraging results. Derivatives of anacardic acid and cardanol from CNSL were found to have moderate activity against trypanosomes with no activity against Plasmodium falciparum. Almost 63% of synthesized 1,4-dihydroxyanthraquinone derivatives displayed activity against trypanosomes. The in-vitro evaluation and the in silico molecular docking studies revealed that 1,4-dihydroxyanthraquinone derivatives can be potential drug-like candidates active against T.brucei parasites (IC50 = 0.70-1.20 μM). Only four 1,4- iv dihydroxyanthraquinone derivatives with thiosemicarbazone, chloride, pyrrole, and diethanolamine functionality displayed activity against Plasmodium falciparum (IC50 = 3.17- 14.36 μM). In-vitro evaluated of test compounds against antibacterial screen and cytotoxicity effects significantly showed that 2-hydroxy-6-pentadecylbenzoic acid 63a and 2-((2- chlorophenyl)(piperazin-1-yl) methyl)-1,4-dihydroxyanthracene-9,10-dione 78 have potency against Staphylococcus aureus and reduced the viability of the cells below 20% at an initial concentration of 50 μg/mL. Only 1,4-dihydroxyanthraquinone derivatives with thiosemicarbazone 76, piperazine 78, and diethanolamine 80 motifs were active against HeLa cells and reduced the viability of cells below 20% at a concentration of 50 μg/mL. In conclusion, this current reported study has generated useful knowlege on the applicability of the agro-waste CNSL as an agent active against trypanosomiasis but also as a low-cost starting material to synthesize hydroxy anthraquinones. The study has further given an overview to the understanding of the medicinal value 1,4-dihydroxyanthraquinone derivatives as promising candidates towards developing drugs suitable for treating neglected tropical diseases particularly trypanosomiasis. , Thesis (PhD) -- Faculty of Science, Chemistry, 2022
- Full Text:
- Date Issued: 2022-10-14