Understanding the underlying resistance mechanism of Mycobacterium tuberculosis against Rifampicin by analyzing mutant DNA - directed RNA polymerase proteins via bioinformatics approaches
- Authors: Monama, Mokgerwa Zacharia
- Date: 2020
- Subjects: Mycobacterium tuberculosis , Rifampin , Drug resistance , Homology (Biology) , Tuberculosis -- Chemotherapy
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/167508 , vital:41487
- Description: Tuberculosis or TB is an airborne disease caused by the non-motile bacilli, Mycobacterium tuberculosis (MTB). There are two main forms of TB, namely, latent TB or LTB, asymptomatic and non-contagious version which according to the World Health Organization (WHO) is estimated to afflict over a third of the world’s population; and active TB or ATB, a symptomatic and contagious version which continues to spread, affecting millions worldwide. With the already high reported prevalence of TB, the emergence of drug-resistant strains has prompted the development of novel approaches to enhance the efficacy of known drugs and a desperate search for novel compounds to combat MTB infections. It was for this very purpose that this study was conducted. A look into the resistance mechanism of Rifampicin (Rifampin or RIF), one of the more potent first-line drugs, might prove beneficial in predicting the consequence of an introduced mutation (which usually occur as single nucleotide polymorphisms or SNPs) and perhaps even overcome it using appropriate therapeutic interventions that improve RIF’s efficacy. To accomplish this task, models of acceptable quality were generated for the WT and clinically relevant, RIF resistance conferring, SNPs occurring at codon positions D516, H526 and S531 (E .coli numbering system) using MODELLER. The models were accordingly ranked using GA341 and z-DOPE score, and subsequently validated with QMEAN, PROCHECK and VERIFY3D. MD simulations spanning 100 ns were run for RIF-bound (complex) and RIF-free (holo) DNA-directed RNA polymerase (DDRP) protein systems for the WT and SNP mutants using GROMACS. The MD frames were analyzed using RMSD, Rg and RMSF. For further analysis, MD-TASK was used to analyze the calculated dynamic residue networks (DRNs) from the generated MD frames, determining both change in average shortest path (ΔL) and betweenness centrality (ΔBC). The RMSD analysis revealed that all of the SNP complex models displayed a level instability higher than that of the WT complex. A majority of the SNP complex models were also observed to have similar compactness to the WT holo when looking at the calculated Rg. The RMSF results also hinted towards possible physiological consequences of the mutations (generally referred to as a fitness cost) highlighted by the increased fluctuations of the zinc-binding domain and the MTB SI α helical coiled coil. For the first time, to the knowledge of the authors, DRN analysis was employed for the DDRP protein for both holo and complex systems, revealing insightful information about the residues that play a key role in the change in distance between residue pairs along with residues that play an essential role in protein communication within the calculated RIN. Overall, the data supported the conclusions drawn by a recent study that only concentrated on RIF-resistance in rpoB models which suggested that the binding pocket for the SNP models may result in the changed coordination of RIF which may be the main contributor to its impaired efficacy.
- Full Text:
- Date Issued: 2020
- Authors: Monama, Mokgerwa Zacharia
- Date: 2020
- Subjects: Mycobacterium tuberculosis , Rifampin , Drug resistance , Homology (Biology) , Tuberculosis -- Chemotherapy
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/167508 , vital:41487
- Description: Tuberculosis or TB is an airborne disease caused by the non-motile bacilli, Mycobacterium tuberculosis (MTB). There are two main forms of TB, namely, latent TB or LTB, asymptomatic and non-contagious version which according to the World Health Organization (WHO) is estimated to afflict over a third of the world’s population; and active TB or ATB, a symptomatic and contagious version which continues to spread, affecting millions worldwide. With the already high reported prevalence of TB, the emergence of drug-resistant strains has prompted the development of novel approaches to enhance the efficacy of known drugs and a desperate search for novel compounds to combat MTB infections. It was for this very purpose that this study was conducted. A look into the resistance mechanism of Rifampicin (Rifampin or RIF), one of the more potent first-line drugs, might prove beneficial in predicting the consequence of an introduced mutation (which usually occur as single nucleotide polymorphisms or SNPs) and perhaps even overcome it using appropriate therapeutic interventions that improve RIF’s efficacy. To accomplish this task, models of acceptable quality were generated for the WT and clinically relevant, RIF resistance conferring, SNPs occurring at codon positions D516, H526 and S531 (E .coli numbering system) using MODELLER. The models were accordingly ranked using GA341 and z-DOPE score, and subsequently validated with QMEAN, PROCHECK and VERIFY3D. MD simulations spanning 100 ns were run for RIF-bound (complex) and RIF-free (holo) DNA-directed RNA polymerase (DDRP) protein systems for the WT and SNP mutants using GROMACS. The MD frames were analyzed using RMSD, Rg and RMSF. For further analysis, MD-TASK was used to analyze the calculated dynamic residue networks (DRNs) from the generated MD frames, determining both change in average shortest path (ΔL) and betweenness centrality (ΔBC). The RMSD analysis revealed that all of the SNP complex models displayed a level instability higher than that of the WT complex. A majority of the SNP complex models were also observed to have similar compactness to the WT holo when looking at the calculated Rg. The RMSF results also hinted towards possible physiological consequences of the mutations (generally referred to as a fitness cost) highlighted by the increased fluctuations of the zinc-binding domain and the MTB SI α helical coiled coil. For the first time, to the knowledge of the authors, DRN analysis was employed for the DDRP protein for both holo and complex systems, revealing insightful information about the residues that play a key role in the change in distance between residue pairs along with residues that play an essential role in protein communication within the calculated RIN. Overall, the data supported the conclusions drawn by a recent study that only concentrated on RIF-resistance in rpoB models which suggested that the binding pocket for the SNP models may result in the changed coordination of RIF which may be the main contributor to its impaired efficacy.
- Full Text:
- Date Issued: 2020
In-silico analysis of Plasmodium falciparum Hop protein and its interactions with Hsp70 and Hsp90
- Authors: Clitheroe, Crystal-Leigh
- Date: 2013
- Subjects: Plasmodium falciparum , Heat shock proteins , Molecular chaperones , Homology (Biology) , Protein-protein interactions , Malaria -- Chemotherapy
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: vital:3896 , http://hdl.handle.net/10962/d1003819 , Plasmodium falciparum , Heat shock proteins , Molecular chaperones , Homology (Biology) , Protein-protein interactions , Malaria -- Chemotherapy
- Description: A lessor understood co-chaperone, the Hsp70/Hsp90 organising protein (Hop), has been found to play an important role in modulating the activity and co-interaction of two essential chaperones; Hsp90 and Hsp70. The best understood aspects of Hop so far indicate that residues in the concave surfaces of the three tetratricopeptide repeat (TPR) domains in the protein bind selectively to the C-terminal motifs of Hsp70 and Hsp90. Recent research suggests that P. falciparum Hop (PfHop), PfHsp90 and PfHsp70 do interact and form complex in the P. falciparum trophozooite and are overexpressed in this infective stage. However, there has been almost no computational research on malarial Hop protein in complex with other malarial Hsps.The current work has focussed on several aspects of the in-silico characterisation of PfHop, including an in-depth multiple sequence alignment and phylogenetic analysis of the protein; which showed that Hop is very well conserved across a wide range of available phyla (four Kingdoms, 60 species). Homology modelling was employed to predict several protein structures for these interactions in P. falciparum, as well as predict structures of the relevant TPR domains of Human Hop (HsHop) in complex with its own Hsp90 and Hsp70 C-terminal peptide partners for comparison. Protein complex interaction analyses indicate that concave TPR sites bound to the C-terminal motifs of partner proteins are very similar in both species, due to the excellent conservation of the TPR domain’s “double carboxylate binding clamp”. Motif analysis was combined with phylogenetic trees and structure mapping in novel ways to attain more information on the evolutionary conservation of important structural and functional sites on Hop. Alternative sites of interaction between Hop TPR2 and Hsp90’s M and C domains are distinctly less well conserved between the two species, but still important to complex formation, making this a likely interaction site for selective drug targeting. Binding and interaction energies for all modelled complexes have been calculated; indicating that all HsHop TPR domains have higher affinities for their respective C-terminal partners than do their P. falciparum counterparts. An alternate motif corresponding to the C-terminal motif of PfHsp70-x (exported to the infected erythrocyte cytosol) in complex with both human and malarial TPR1 and TPR2B domains was analysed, and these studies suggest that the human TPR domains have a higher affinity for this motif than do the respective PfHop TPR domains. This may indicate potential for a cross species protein interaction to take place, as PfHop is not transported to the human erythrocyte cytosol.
- Full Text:
- Date Issued: 2013
- Authors: Clitheroe, Crystal-Leigh
- Date: 2013
- Subjects: Plasmodium falciparum , Heat shock proteins , Molecular chaperones , Homology (Biology) , Protein-protein interactions , Malaria -- Chemotherapy
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: vital:3896 , http://hdl.handle.net/10962/d1003819 , Plasmodium falciparum , Heat shock proteins , Molecular chaperones , Homology (Biology) , Protein-protein interactions , Malaria -- Chemotherapy
- Description: A lessor understood co-chaperone, the Hsp70/Hsp90 organising protein (Hop), has been found to play an important role in modulating the activity and co-interaction of two essential chaperones; Hsp90 and Hsp70. The best understood aspects of Hop so far indicate that residues in the concave surfaces of the three tetratricopeptide repeat (TPR) domains in the protein bind selectively to the C-terminal motifs of Hsp70 and Hsp90. Recent research suggests that P. falciparum Hop (PfHop), PfHsp90 and PfHsp70 do interact and form complex in the P. falciparum trophozooite and are overexpressed in this infective stage. However, there has been almost no computational research on malarial Hop protein in complex with other malarial Hsps.The current work has focussed on several aspects of the in-silico characterisation of PfHop, including an in-depth multiple sequence alignment and phylogenetic analysis of the protein; which showed that Hop is very well conserved across a wide range of available phyla (four Kingdoms, 60 species). Homology modelling was employed to predict several protein structures for these interactions in P. falciparum, as well as predict structures of the relevant TPR domains of Human Hop (HsHop) in complex with its own Hsp90 and Hsp70 C-terminal peptide partners for comparison. Protein complex interaction analyses indicate that concave TPR sites bound to the C-terminal motifs of partner proteins are very similar in both species, due to the excellent conservation of the TPR domain’s “double carboxylate binding clamp”. Motif analysis was combined with phylogenetic trees and structure mapping in novel ways to attain more information on the evolutionary conservation of important structural and functional sites on Hop. Alternative sites of interaction between Hop TPR2 and Hsp90’s M and C domains are distinctly less well conserved between the two species, but still important to complex formation, making this a likely interaction site for selective drug targeting. Binding and interaction energies for all modelled complexes have been calculated; indicating that all HsHop TPR domains have higher affinities for their respective C-terminal partners than do their P. falciparum counterparts. An alternate motif corresponding to the C-terminal motif of PfHsp70-x (exported to the infected erythrocyte cytosol) in complex with both human and malarial TPR1 and TPR2B domains was analysed, and these studies suggest that the human TPR domains have a higher affinity for this motif than do the respective PfHop TPR domains. This may indicate potential for a cross species protein interaction to take place, as PfHop is not transported to the human erythrocyte cytosol.
- Full Text:
- Date Issued: 2013
The spatial evolution of the chemotaxis proteins of the Bacillus subtilis group
- Yssel, Anna Elizabeth Johanna
- Authors: Yssel, Anna Elizabeth Johanna
- Date: 2011
- Subjects: Chemotaxis , Bacillus subtilis , Bacillus (Bacteria) , Homology (Biology) , Plants -- Microbiology
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: vital:4027 , http://hdl.handle.net/10962/d1004087 , Chemotaxis , Bacillus subtilis , Bacillus (Bacteria) , Homology (Biology) , Plants -- Microbiology
- Description: The aim of this work was to study spatial evolution of the chemotaxis proteins of a group of plant-associated soil-dwelling bacteria vernacularly referred to as the B. subtilis group. This was achieved by creating homology models for the chemotaxis proteins if a suitable template was available, and by analysing the selective forces (positive, purifying or neutral) acting upon the chemotaxis proteins. Chemotaxis is the phenomenon in which bacteria direct their movement towards more favourable conditions, and is critical for processes such as obtaining nutrients, escaping toxic compounds, host colonization and bio-film formation. Members of the B. subtilis group exhibit different preferences for certain host plants, and it is therefore feasible that their chemotactic machinery are fine-tuned to respond optimally to the conditions of the various niches that the strains inhabit. Homology models were inferred for the plant growth promoting B. amyloliquefaciens FZB42 proteins CheB, CheC, CheD, CheR, CheW and CheY. The interactions between: CheC-CheD, the P1 and P2 domains of CheA with CheY and CheB, and the P4 and P5 domains of CheA with CheW were also modelled. The hydrophobic interactions contributing to intra- and inter-protein contacts were analysed. The models of the interactions between CheB and the various domains of CheA are of particular interest, because to date no structures have been solved that show an interaction between a histidine kinase (such as CheA) and a multidomain response regulator (such as CheB). Furthermore, evidence that phospho-CheB may inhibit the formation of phospho-CheY by competitively binding to the P2 domain of CheA is also presented. Proteins were analysed to determine if individual amino acid sites are under positive, neutral or purifying selection. The Methyl Accepting Chemotaxis Proteins (MCPs), CheA and CheV were also analyzed, but due to a lack of suitable templates, no homology models were constructed. Site-specific positive and purifying selection were estimated by comparing the ratios of non-synonymous to synonymous substitutions at each site in the sequences for the chemotaxis proteins as well as for the receptors McpA, McpB, and McpC. Homology models were coloured according to intensity of selective forces. It was found that the chemotaxis proteins of member of the B. subtilis group are under strong evolutionary constraints, hence it is unlikely that positive selection in these proteins are responsible for the differences in habitat preference that these organism exhibit.
- Full Text:
- Date Issued: 2011
- Authors: Yssel, Anna Elizabeth Johanna
- Date: 2011
- Subjects: Chemotaxis , Bacillus subtilis , Bacillus (Bacteria) , Homology (Biology) , Plants -- Microbiology
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: vital:4027 , http://hdl.handle.net/10962/d1004087 , Chemotaxis , Bacillus subtilis , Bacillus (Bacteria) , Homology (Biology) , Plants -- Microbiology
- Description: The aim of this work was to study spatial evolution of the chemotaxis proteins of a group of plant-associated soil-dwelling bacteria vernacularly referred to as the B. subtilis group. This was achieved by creating homology models for the chemotaxis proteins if a suitable template was available, and by analysing the selective forces (positive, purifying or neutral) acting upon the chemotaxis proteins. Chemotaxis is the phenomenon in which bacteria direct their movement towards more favourable conditions, and is critical for processes such as obtaining nutrients, escaping toxic compounds, host colonization and bio-film formation. Members of the B. subtilis group exhibit different preferences for certain host plants, and it is therefore feasible that their chemotactic machinery are fine-tuned to respond optimally to the conditions of the various niches that the strains inhabit. Homology models were inferred for the plant growth promoting B. amyloliquefaciens FZB42 proteins CheB, CheC, CheD, CheR, CheW and CheY. The interactions between: CheC-CheD, the P1 and P2 domains of CheA with CheY and CheB, and the P4 and P5 domains of CheA with CheW were also modelled. The hydrophobic interactions contributing to intra- and inter-protein contacts were analysed. The models of the interactions between CheB and the various domains of CheA are of particular interest, because to date no structures have been solved that show an interaction between a histidine kinase (such as CheA) and a multidomain response regulator (such as CheB). Furthermore, evidence that phospho-CheB may inhibit the formation of phospho-CheY by competitively binding to the P2 domain of CheA is also presented. Proteins were analysed to determine if individual amino acid sites are under positive, neutral or purifying selection. The Methyl Accepting Chemotaxis Proteins (MCPs), CheA and CheV were also analyzed, but due to a lack of suitable templates, no homology models were constructed. Site-specific positive and purifying selection were estimated by comparing the ratios of non-synonymous to synonymous substitutions at each site in the sequences for the chemotaxis proteins as well as for the receptors McpA, McpB, and McpC. Homology models were coloured according to intensity of selective forces. It was found that the chemotaxis proteins of member of the B. subtilis group are under strong evolutionary constraints, hence it is unlikely that positive selection in these proteins are responsible for the differences in habitat preference that these organism exhibit.
- Full Text:
- Date Issued: 2011
- «
- ‹
- 1
- ›
- »