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
Quinolone-Pyrazinamide Derivatives: synthesis, characterisation, in silico ADME analysis and in vitro biological evaluation against Mycobacterium tuberculosis
- Authors: Rukweza, Kudakwashe Gerald
- Date: 2023-10-13
- Subjects: Quinolone antibacterial agents , Mycobacterium tuberculosis , Antitubercular agents , Tuberculosis Chemotherapy , Drug resistance , Moxifloxacin , Isoniazid
- Language: English
- Type: Academic theses , Master's theses , text
- Identifier: http://hdl.handle.net/10962/390901 , vital:68596
- Description: Tuberculosis is one of the leading causes of death worldwide caused by an infectious species, Mycobacterium tuberculosis (Mtb). Some of the factors that contribute to the prevalence of this disease include the complexity of diagnosis, prolonged period of therapy, side effects associated with current TB drugs, the prevalence of resistance against the current treatment options and a high incidence of co-infection with HIV/AIDS. Thus, there is a need for new alternative drugs to provide safer and shorter treatment therapy options that are not susceptible to the development of drug resistance. In this project, we focus our attention on the quinolone pharmacophore. Quinolones are currently used as alternative options in the treatment of resistant strains of Mtb. Previous work pertaining to quinolone-isoniazid hybrid compounds showed promising in vitro activity against the H37Rv strain of Mtb and served as the inspiration to pursue this project. The journey commenced with the synthesis of quinolone-pyrazinamide hybrid compounds (Figure 3.1). These compounds were synthesised, through the attachment of the quinolone and the pyrazinamide entity through a hydrazine linker. The synthesised compounds were purified, and their structural identity confirmed using common spectroscopic techniques including 1H and 13C NMR, infra-red (IR) and mass spectrometry. In vitro biological assays were performed by testing for the activity against the H37RvMA strain of Mtb. The bioassays were performed in triplicates to ensure the accuracy of the results. Moxifloxacin and isoniazid were tested as control compounds. Finally, the resultant compounds were profiled in silico for physicochemical and ADMET properties using open access software SwissADME. All the synthesised compounds 3.8a-f showed no activity against H37RvMA. In most cases, the resulting compounds showed minimal to no activity (MICs ≥ 57.3 μM) in all three media. During the in vitro studies, the compounds showed significant precipitation in the media over time suggesting poor aqueous solubility. The SwissADME analysis of these compounds indicated poor solubility in aqueous media, which is likely linked to their molecular size and complexity. Despite poor aqueous solubility, compounds 3.8a-f showed acceptable physicochemical properties and ADME parameters. No PAINs (Pan-assay interference compounds) were observed. Minimal to no interaction with CYP enzymes were predicted. Most of the compounds were compatible with the Lipinski’s rules of five. , Thesis (MSc) -- Faculty of Science, Pharmacy, 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
Design, formulation and evaluation of liposomes co-loaded with human serum albumin and rifampicin
- Authors: Bapolisi, Alain Murhimalika
- Date: 2020
- Subjects: Liposomes , Rifampin , Antitubercular agents , Serum albumin , Albumins , Tuberculosis -- Treatment
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/163179 , vital:41016
- Description: Tuberculosis (TB) is a devastating infectious disease caused by Mycobacterium tuberculosis and is the leading cause of death from a single infectious agent. The high morbidity and mortality rates of TB are partly due to factors such as the lengthy regimen (of 6–24 months), the development of drug resistance, and the pathogen location within the macrophages. These, with poor physiochemical properties of existing drugs hamper the effectiveness of the treatment despite the existence of potent antibiotics such as Rifampicin (Rif). Hydrophobicity plagues many drugs, including Rif, which are then particularly affected due to inherently poor intracellular availability. Novel drug delivery approaches are therefore needed in order to optimize the cytotoxic potential of said antitubercular drugs. To improve the bioavailability of hydrophobic drugs, numerous delivery strategies have been developed. Amongst these, the coordination of cytotoxic drugs to therapeutic proteins have shown some success for improved efficacy in the management of illnesses including infectious diseases. Of therapeutic proteins, Human Serum Albumin (HSA) is an attractive drug carrier with interesting benefits such as low immunogenicity, antioxidant properties and improving cellular uptake of drugs through HSA-specific binding sites which are expressed on most cells including macrophages, where M. tuberculosis often resides. Hence, combination of Rif to HSA (Rif-HSA) seems a promising approach for improved intracellular delivery of Rif. However, the in vivo stability of colloidal protein-based therapeutics is mostly challenging and an effective vehicle is needed to control the biological fate of such conjugates. Liposomes seem to be appropriate carriers for the Rif-HSA complex due to their reputable applicability for encapsulating diverse materials (i.e., hydrophobic and hydrophilic compounds or small and complex molecules) and preventing chemical and biological degradation of the cargo. Therefore, the main objective of this study was to simultaneously encapsulate Rif and HSA in liposomes, which, to the best of our knowledge, has not been done before. The dual liposomes (Rif-HSA-lip) were made by a modified “Reverse Phase Evaporation” method (REV), following a Design of Experiments (DOE) approach to determine which factors impact the formulation. In addition, liposomes were made from crude soybean lecithin (CSL), rather than expensive and highly purified lipids. iv The liposomes were fully characterised, and the encapsulation efficiency (î) was monitored using high-performance liquid chromatography (HPLC). The results were correlated with factors such as organic and aqueous phase composition, as well as the in vitro release profile of Rif. Transmission electron microscopy (TEM) results confirmed the formation of spherical dual liposomes nanoparticles of roughly 200 nm. Dynamic light scattering (DLS) and Zeta potential measurements showed a negative charge (<–45 mV) and with satisfactory polydispersity (PDI<0.5). HSA dramatically improved the aqueous solubility of Rif (from1.9 mg/ml in water to around 4.3 mg/ml in HSA 10% solution) mainly due to Rif-HSA hydrophobic interactions. This resulted in a good î of almost 60% for Rif, despite the presence of bulky HSA in the lipid bilayer. These details were confirmed using proton nuclear magnetic resonance (1H NMR) and Fourier-transform infrared spectroscopy (FTIR). Furthermore, energy dispersive X-ray (EDX) and DLS data suggested the presence of HSA poking out on the surface of liposomes, which is encouraging for potential targeted delivery in the future. The in vitro release studies also depicted a substantial improvement in the diffusion of Rif in dual liposomes versus free Rif, from 65% after 12 hours for free Rif to 95% after only 5 hours for Rif- HSA-lip. Finally, stability studies conducted over 30 days at room temperature, showed that the freeze-dried formulations of Rif-HSA-lip exhibited good shelf stability over liposomes with no HSA. This study represents an illustrative example of co-loading of antibiotics and proteins into liposomes, which could encourage further development of novel nanoparticulate tools for the effective management of both drug-susceptible and -resistant infectious diseases such as TB.
- Full Text:
- Date Issued: 2020