In silico characterization of missense mutations in infectious diseases: case studies of tuberculosis and COVID-19
- Authors: Barozi, Victor
- Date: 2023-10-13
- Subjects: Microbial mutation , COVID-19 (Disease) , Drug resistance in microorganisms , Antitubercular agents , Tuberculosis , Molecular dynamics , Single nucleotide polymorphisms
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/431626 , vital:72791 , DOI 10.21504/10962/431626
- Description: One of the greatest challenges facing modern medicine and the global public health today is antimicrobial drug resistance (AMR). This “silent pandemic,” as coined by the world health organization (WHO), is steadily increasing with an estimated 4.95 million mortalities attributed to AMR in 2019, 1.27 million of which were directly linked to AMR. Some of the contributors to AMR include self-prescription, drug overuse, sub-optimal drug prescriptions by health workers, and inaccessibility to drugs, especially in remote areas, which leads to poor adherence. The situation is aggravated by the upsurge of new zoonotic infections like the coronavirus disease 2019, which present unique challenges and take the bulk of resources hence stunting the fight against AMR. Quite alarming still is our current antimicrobial arsenal, which hasn’t had any novel antimicrobial drug discovery/addition, of a new class, since the 1980s. This puts a burden on the existing broad-spectrum antimicrobial drugs which are already struggling against multi-drug resistant strains like multi-drug resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB). Besides the search for new antimicrobial agents, the other avenue for addressing AMR is studying drug resistance mechanisms, especially single nucleotide polymorphisms (SNPs), that change drug target characteristics. With the advancement of computational power and data storage resources, computational approaches can be applied in mutational studies to provide insight into the drug resistance mechanisms with an aim to inform future drug design and development. Therefore, in the first part of this thesis, we employ integrative in silico approaches, including 3D structure modeling, molecular dynamic (MD) simulations, comparative essential dynamics (ED), and protein network analysis approaches i.e., dynamic residue network (DRN) analysis to decipher drug resistance mechanisms in tuberculosis (TB). This involved an investigation of the drug resistance mutations in the catalase-peroxidase (KatG) and pyrazinamidase (MtPncA) enzymes which are responsible for activation of TB first-line drugs; Isoniazid (INH) and Pyrazinamide (PZA), respectively. In the case of KatG, eleven high confidence (HC) KatG mutations associated with a high prevalence of phenotypic INH resistance were identified and their 3D structures modeled before subjecting them to MD simulations. Global analysis showed an unstable KatG structure and active site environment in the mutants compared to the wildtype. Active site dynamics in the mutants compromised cofactor (heme) interactions resulting in less bonds/interactions compared to the wildtype. Given the importance of the heme, reduced interactions affect enzyme function. Trajectory analysis also showed asymmetric protomer behavior both in the wildtype and mutant systems. DRN analysis identified the KatG dimerization domain and C-terminal domain as functionally important and influential in the enzyme function as per betweenness centrality and eigenvector centrality distribution. In the case of the MtPncA enzyme, our main focus was on understanding the MtPncA binding ability of Nicotinamide (an analogue of PZA) in comparison to PZA, especially in the presence of 82 resistance conferring MtPncA mutations. Like in KatG, the mutant structures were modeled and subjected to MD simulations and analysis. Interestingly, more MtPncA mutants favored NAM interactions compared to PZA i.e., 34 MtPncA mutants steadily coordinated NAM compared to 21 in the case of PZA. Trajectory and ligand interaction analysis showed how increased active site lid loop dynamics affect the NAM binding, especially in the systems with the active site mutations i.e., H51Y, W68R, C72R, L82R, K96N, L159N, and L159R. This led to fewer protein-ligand interactions and eventually ligand ejection. Network analysis further identified the protein core, metal binding site (MBS), and substrate binding site as the most important regions of the enzyme. Furthermore, the degree of centrality analysis showed how specific MtPncA mutations i.e., C14H, F17D, and T412P, interrupt intra-protein communication from the MtPncA core to the MBS, affecting enzyme activity. The analysis of KatG and MtPncA enzyme mutations not only identified the effects of mutations on enzyme behaviour and communication, but also established a framework of computational approaches that can be used for mutational studies in any protein. Besides AMR, the continued encroachment of wildlife habitats due to population growth has exposed humans to wildlife pathogens leading to zoonotic diseases, a recent example being coronavirus disease 2019 (COVID-19). In the second part of the thesis, the established computational approaches in Part 1, were employed to investigate the changes in inter-protein interactions and communication patterns between the severe acute respiratory coronavirus 2 (SARS-CoV-2) with the human host receptor protein (ACE2: angiotensin-converting enzyme 2) consequent to mutations in the SARS-CoV-2 receptor binding domain (RBD). Here, the focus was on RBD mutations of the Omicron sub-lineages. We identified four Omicron-sub lineages with RBD mutations i.e., BA.1, BA.2, BA.3 and BA.4. Each sub-lineage mutations were modeled into RBD structure in complex with the hACE2. MD analysis of the RBD-hACE2 complex highlighted how the RBD mutations change the conformational flexibility of both the RBD and hACE2 compared to the wildtype (WT). Furthermore, DRN analysis identified novel allosteric paths composed of residues with high betweenness and eigenvector centralities linking the RBD to the hACE2 in both the wildtype and mutant systems. Interestingly, these paths were modified with the progression of Omicron sub-lineages, highlighting how the virus evolution affects protein interaction. Lastly, the effect of mutations on S RBD and hACE2 interaction was investigated from the hACE2 perspective by focusing on mutations in the hACE2 protein. Here, naturally occurring hACE2 polymorphisms in African populations i.e., S19P, K26R, M82I, K341R, N546D, and D597Q, were identified and their effects on RBD-hACE2 interactions investigated in presence of the Omicron BA.4/5 RBD mutations. The hACE2 polymorphisms subtly affected the complex dynamics; however, RBD-hACE2 interaction analysis showed that hACE2 mutations effect the complex formation and interaction. Here, the K26R mutation favored RBD-hACE2 interactions, whereas S19P resulted in fewer inter-protein interactions than the reference system. The M82I mutation resulted in a higher RBD-hACE2 binding energy compared to the wildtype meaning that the mutation might not favor RBD binding to the hACE2. On the other hand, K341R had the most RBD-hACE2 interactions suggesting that it probably favors RBD binding to the hACE2. N546D and D597Q had diminutive differences to the reference system. Interestingly, the network of high betweenness centrality residues linking the two proteins, as seen in the previous paragraph, were maintained/modified in presence of hACE2 mutations. HACE2 mutations also changed the enzyme network patterns resulting in a concentration of high eigenvector centrality residues around the zinc-binding and active site region, ultimately influencing the enzyme functionality. Altogether, the thesis highlights fundamental structural and network changes consequent to mutations both in TB and COVID-19 proteins of interest using in silico approaches. These approaches not only provide a new context on impact of mutations in TB and COVID target proteins, but also presents a framework that be implemented in other protein mutation studies. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2023
- Full Text:
- Date Issued: 2023-10-13
- Authors: Barozi, Victor
- Date: 2023-10-13
- Subjects: Microbial mutation , COVID-19 (Disease) , Drug resistance in microorganisms , Antitubercular agents , Tuberculosis , Molecular dynamics , Single nucleotide polymorphisms
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/431626 , vital:72791 , DOI 10.21504/10962/431626
- Description: One of the greatest challenges facing modern medicine and the global public health today is antimicrobial drug resistance (AMR). This “silent pandemic,” as coined by the world health organization (WHO), is steadily increasing with an estimated 4.95 million mortalities attributed to AMR in 2019, 1.27 million of which were directly linked to AMR. Some of the contributors to AMR include self-prescription, drug overuse, sub-optimal drug prescriptions by health workers, and inaccessibility to drugs, especially in remote areas, which leads to poor adherence. The situation is aggravated by the upsurge of new zoonotic infections like the coronavirus disease 2019, which present unique challenges and take the bulk of resources hence stunting the fight against AMR. Quite alarming still is our current antimicrobial arsenal, which hasn’t had any novel antimicrobial drug discovery/addition, of a new class, since the 1980s. This puts a burden on the existing broad-spectrum antimicrobial drugs which are already struggling against multi-drug resistant strains like multi-drug resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB). Besides the search for new antimicrobial agents, the other avenue for addressing AMR is studying drug resistance mechanisms, especially single nucleotide polymorphisms (SNPs), that change drug target characteristics. With the advancement of computational power and data storage resources, computational approaches can be applied in mutational studies to provide insight into the drug resistance mechanisms with an aim to inform future drug design and development. Therefore, in the first part of this thesis, we employ integrative in silico approaches, including 3D structure modeling, molecular dynamic (MD) simulations, comparative essential dynamics (ED), and protein network analysis approaches i.e., dynamic residue network (DRN) analysis to decipher drug resistance mechanisms in tuberculosis (TB). This involved an investigation of the drug resistance mutations in the catalase-peroxidase (KatG) and pyrazinamidase (MtPncA) enzymes which are responsible for activation of TB first-line drugs; Isoniazid (INH) and Pyrazinamide (PZA), respectively. In the case of KatG, eleven high confidence (HC) KatG mutations associated with a high prevalence of phenotypic INH resistance were identified and their 3D structures modeled before subjecting them to MD simulations. Global analysis showed an unstable KatG structure and active site environment in the mutants compared to the wildtype. Active site dynamics in the mutants compromised cofactor (heme) interactions resulting in less bonds/interactions compared to the wildtype. Given the importance of the heme, reduced interactions affect enzyme function. Trajectory analysis also showed asymmetric protomer behavior both in the wildtype and mutant systems. DRN analysis identified the KatG dimerization domain and C-terminal domain as functionally important and influential in the enzyme function as per betweenness centrality and eigenvector centrality distribution. In the case of the MtPncA enzyme, our main focus was on understanding the MtPncA binding ability of Nicotinamide (an analogue of PZA) in comparison to PZA, especially in the presence of 82 resistance conferring MtPncA mutations. Like in KatG, the mutant structures were modeled and subjected to MD simulations and analysis. Interestingly, more MtPncA mutants favored NAM interactions compared to PZA i.e., 34 MtPncA mutants steadily coordinated NAM compared to 21 in the case of PZA. Trajectory and ligand interaction analysis showed how increased active site lid loop dynamics affect the NAM binding, especially in the systems with the active site mutations i.e., H51Y, W68R, C72R, L82R, K96N, L159N, and L159R. This led to fewer protein-ligand interactions and eventually ligand ejection. Network analysis further identified the protein core, metal binding site (MBS), and substrate binding site as the most important regions of the enzyme. Furthermore, the degree of centrality analysis showed how specific MtPncA mutations i.e., C14H, F17D, and T412P, interrupt intra-protein communication from the MtPncA core to the MBS, affecting enzyme activity. The analysis of KatG and MtPncA enzyme mutations not only identified the effects of mutations on enzyme behaviour and communication, but also established a framework of computational approaches that can be used for mutational studies in any protein. Besides AMR, the continued encroachment of wildlife habitats due to population growth has exposed humans to wildlife pathogens leading to zoonotic diseases, a recent example being coronavirus disease 2019 (COVID-19). In the second part of the thesis, the established computational approaches in Part 1, were employed to investigate the changes in inter-protein interactions and communication patterns between the severe acute respiratory coronavirus 2 (SARS-CoV-2) with the human host receptor protein (ACE2: angiotensin-converting enzyme 2) consequent to mutations in the SARS-CoV-2 receptor binding domain (RBD). Here, the focus was on RBD mutations of the Omicron sub-lineages. We identified four Omicron-sub lineages with RBD mutations i.e., BA.1, BA.2, BA.3 and BA.4. Each sub-lineage mutations were modeled into RBD structure in complex with the hACE2. MD analysis of the RBD-hACE2 complex highlighted how the RBD mutations change the conformational flexibility of both the RBD and hACE2 compared to the wildtype (WT). Furthermore, DRN analysis identified novel allosteric paths composed of residues with high betweenness and eigenvector centralities linking the RBD to the hACE2 in both the wildtype and mutant systems. Interestingly, these paths were modified with the progression of Omicron sub-lineages, highlighting how the virus evolution affects protein interaction. Lastly, the effect of mutations on S RBD and hACE2 interaction was investigated from the hACE2 perspective by focusing on mutations in the hACE2 protein. Here, naturally occurring hACE2 polymorphisms in African populations i.e., S19P, K26R, M82I, K341R, N546D, and D597Q, were identified and their effects on RBD-hACE2 interactions investigated in presence of the Omicron BA.4/5 RBD mutations. The hACE2 polymorphisms subtly affected the complex dynamics; however, RBD-hACE2 interaction analysis showed that hACE2 mutations effect the complex formation and interaction. Here, the K26R mutation favored RBD-hACE2 interactions, whereas S19P resulted in fewer inter-protein interactions than the reference system. The M82I mutation resulted in a higher RBD-hACE2 binding energy compared to the wildtype meaning that the mutation might not favor RBD binding to the hACE2. On the other hand, K341R had the most RBD-hACE2 interactions suggesting that it probably favors RBD binding to the hACE2. N546D and D597Q had diminutive differences to the reference system. Interestingly, the network of high betweenness centrality residues linking the two proteins, as seen in the previous paragraph, were maintained/modified in presence of hACE2 mutations. HACE2 mutations also changed the enzyme network patterns resulting in a concentration of high eigenvector centrality residues around the zinc-binding and active site region, ultimately influencing the enzyme functionality. Altogether, the thesis highlights fundamental structural and network changes consequent to mutations both in TB and COVID-19 proteins of interest using in silico approaches. These approaches not only provide a new context on impact of mutations in TB and COVID target proteins, but also presents a framework that be implemented in other protein mutation studies. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2023
- Full Text:
- Date Issued: 2023-10-13
In silico identification of selective novel hits against the active site of wild type mycobacterium tuberculosis pyrazinamidase and its mutants
- Authors: Gowo, Prudence
- Date: 2021-04
- Subjects: Mycobacterium tuberculosis , Pyrazinamide , Multidrug resistance , Antitubercular agents , Molecular dynamics , Hydrogen bonding , Ligand binding (Biochemistry) , Dynamic Residue Network
- Language: English
- Type: thesis , text , Masters , MSc
- Identifier: http://hdl.handle.net/10962/178007 , vital:42898
- Description: The World Health Organization declared Tuberculosis a global health emergency and has set a goal to eradicate it by 2035. However, effective treatment and control of the disease is being hindered by the emerging Multi-Drug Resistant and Extensively Drug Resistant strains on the most effective first line prodrug, Pyrazinamide (PZA). Studies have shown that the main cause of PZA resistance is due to mutations in the pncA gene that codes for the target protein Pyrazinamidase (PZase). Therefore, this study aimed to identify novel drug compounds that bind to the active site of wild type PZase and study the dynamics of these potential anti-TB drugs in the mutant systems of PZase. This approach will aid in identifying drugs that may be repurposed for TB therapy and/or designed to counteract PZA resistance. This was achieved by screening 2089 DrugBank compounds against the whole wild type (WT) PZase protein in molecular docking using AutoDOCK4.2. Compound screening based on docking binding energy, hydrogen bonds, molecular weight and active site proximity identified 47 compounds meeting all the set selection criteria. The stability of these compounds were analysed in Molecular Dynamic (MD) simulations and were further studied in PZase mutant systems of A3P, A134V, A146V, D8G, D49A, D49G, D63G, H51P, H137R, L85R, L116R, Q10P, R140S, T61P, V139M and Y103S. Generally, mutant-ligand systems displayed little deviation from the WT systems. The compound systems remained compact, with less fluctuations and more hydrogen bond interactions throughout the simulation (DB00255, DB00655, DB00672, DB00782, DB00977, DB01196, DB04573, DB06414, DB08981, DB11181, DB11760, DB13867, DB13952). From this research study, potential drugs that may be repurposed for TB therapy were identified. Majority of these drugs are currently used in the treatment of hypertension, menopause disorders and inflammation. To further understand the mutant-ligand dynamic systems, calculations such as Dynamic Residue Network (DRN) may be done. Also, the bioactivity of these drugs on Mycobacterium tuberculosis may be studied in wet laboratory, to understand their clinical impart in vivo experiments. , Thesis (MSc) -- Faculty of Science, Biochemistry and Microbiology, 2021
- Full Text:
- Date Issued: 2021-04
- Authors: Gowo, Prudence
- Date: 2021-04
- Subjects: Mycobacterium tuberculosis , Pyrazinamide , Multidrug resistance , Antitubercular agents , Molecular dynamics , Hydrogen bonding , Ligand binding (Biochemistry) , Dynamic Residue Network
- Language: English
- Type: thesis , text , Masters , MSc
- Identifier: http://hdl.handle.net/10962/178007 , vital:42898
- Description: The World Health Organization declared Tuberculosis a global health emergency and has set a goal to eradicate it by 2035. However, effective treatment and control of the disease is being hindered by the emerging Multi-Drug Resistant and Extensively Drug Resistant strains on the most effective first line prodrug, Pyrazinamide (PZA). Studies have shown that the main cause of PZA resistance is due to mutations in the pncA gene that codes for the target protein Pyrazinamidase (PZase). Therefore, this study aimed to identify novel drug compounds that bind to the active site of wild type PZase and study the dynamics of these potential anti-TB drugs in the mutant systems of PZase. This approach will aid in identifying drugs that may be repurposed for TB therapy and/or designed to counteract PZA resistance. This was achieved by screening 2089 DrugBank compounds against the whole wild type (WT) PZase protein in molecular docking using AutoDOCK4.2. Compound screening based on docking binding energy, hydrogen bonds, molecular weight and active site proximity identified 47 compounds meeting all the set selection criteria. The stability of these compounds were analysed in Molecular Dynamic (MD) simulations and were further studied in PZase mutant systems of A3P, A134V, A146V, D8G, D49A, D49G, D63G, H51P, H137R, L85R, L116R, Q10P, R140S, T61P, V139M and Y103S. Generally, mutant-ligand systems displayed little deviation from the WT systems. The compound systems remained compact, with less fluctuations and more hydrogen bond interactions throughout the simulation (DB00255, DB00655, DB00672, DB00782, DB00977, DB01196, DB04573, DB06414, DB08981, DB11181, DB11760, DB13867, DB13952). From this research study, potential drugs that may be repurposed for TB therapy were identified. Majority of these drugs are currently used in the treatment of hypertension, menopause disorders and inflammation. To further understand the mutant-ligand dynamic systems, calculations such as Dynamic Residue Network (DRN) may be done. Also, the bioactivity of these drugs on Mycobacterium tuberculosis may be studied in wet laboratory, to understand their clinical impart in vivo experiments. , Thesis (MSc) -- Faculty of Science, Biochemistry and Microbiology, 2021
- Full Text:
- Date Issued: 2021-04
- «
- ‹
- 1
- ›
- »