Sequence, structure, dynamics, and substrate specificity analyses of bacterial Glycoside Hydrolase 1 enzymes from several activities
- Authors: Veldman, Wayde Michael
- Date: 2022-04-08
- Subjects: Glycosidases , Bioinformatics , Molecular dynamics , Ligands (Biochemistry) , Enzymes , Ligand binding (Biochemistry) , Sequence alignment (Bioinformatics) , Structural bioinformatics
- Language: English
- Type: Doctoral thesis , text
- Identifier: http://hdl.handle.net/10962/233805 , vital:50129 , DOI 10.21504/10962/233810
- Description: Glycoside hydrolase 1 (GH1) enzymes are a ubiquitous family of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. Despite their conserved catalytic domain, these enzymes have many different enzyme activities and/or substrate specificities as a change of only a few residues in the active site can alter their function. Most GH1 active site residues are situated in loop regions, and it is known that enzymes are more likely to develop new functions (broad specificity) if they possess an active site with a high proportion of loops. Furthermore, the GH1 active site consists of several subsites and cooperative binding makes the binding affinity of sites difficult to measure because the properties of one subsite are influenced by the binding of the other subsites. Extensive knowledge of protein-ligand interactions is critical to the comprehension of biology at the molecular level. However, the structural determinants and molecular details of GH1 ligand specificity and affinity are very broad, highly complex, not well understood, and therefore still need to be clarified. The aim of this study was to computationally characterise the activity of three newly solved GH1 crystallographic structures sent to us by our collaborators, and to provide evidence for their ligand-binding specificities. In addition, the differences in structural and biochemical contributions to enzyme specificity and/or function between different GH1 activities/enzymes was assessed, and the sequence/structure/function relationship of several activities of GH1 enzymes was analysed and compared. To accomplish the research aims, sequence analyses involving sequence identity, phylogenetics, and motif discovery were performed. As protein structure is more conserved than sequence, the discovered motifs were mapped to 3D structures for structural analysis and comparisons. To obtain information on enzyme mechanism or mode of action, as well as structure-function relationship, computational methods such as docking, molecular dynamics, binding free energy calculations, and essential dynamics were implemented. These computational approaches can provide information on the active site, binding residues, protein-ligand interactions, binding affinity, conformational change, and most structural or dynamic elements that play a role in enzyme function. The three new structures received from our collaborators are the first GH1 crystallographic structures from Bacillus licheniformis ever determined. As phospho-glycoside compounds were unavailable for purchase for use in activity assays, and as the active sites of the structures were absent of ligand, in silico docking and MD simulations were performed to provide evidence for their GH1 activities and substrate specificities. First though, the amino acid sequences of all known characterised bacterial GH1 enzymes were retrieved from the CAZy database and compared to the sequences of the three new B. licheniformis crystallographic structures which provided evidence of the putative 6Pβ-glucosidase activity of enzyme BlBglH, and dual 6Pβ-glucosidase/6Pβ-galactosidase (dual-phospho) activity of enzymes BlBglB and BlBglC. As all three enzymes were determined to be putative 6Pβ-glycosidase activity enzymes, much of the thesis focused on the overall analysis and comparison of the 6Pβ-glucosidase, 6Pβ-galactosidase, and dual-phospho activities that make up the 6Pβ-glycosidases. The 6Pβ-glycosidase active site residues were identified through consensus of binding interactions using all known 6Pβ-glycosidase PDB structures complexed complete ligand substrates. With regards to the 6Pβ-glucosidase activity, it was found that the L8b loop is longer and forms extra interactions with the L8a loop likely leading to increased L8 loop rigidity which would prevent the displacement of residue Ala423 ensuring a steric clash with galactoconfigured ligands and may engender substrate specificity for gluco-configured ligands only. Also, during molecular dynamics simulations using enzyme BlBglH (6Pβ-glucosidase activity), it was revealed that the favourable binding of substrate stabilises the loops that surround and make up the enzyme active site. Using the BlBglC (dual-phospho activity) enzyme structure with either galacto- (PNP6Pgal) or gluco-configured (PNP6Pglc) ligands, MD simulations in triplicate revealed important details of the broad specificity of dual-phospho activity enzymes. The ligand O4 hydroxyl position is the only difference between PNP6Pgal and PNP6Pgal, and it was found that residues Gln23 and Trp433 bind strongly to the ligand O3 hydroxyl group in the PNP6Pgal-enzyme complex, but to the ligand O4 hydroxyl group in the PNP6Pglc-enzyme complex. Also, His124 formed many hydrogen bonds with the PNP6Pgal O3 hydroxyl group but had none with PNP6Pglc. Alternatively, residues Tyr173, Tyr301, Gln302 and Thr321 formed hydrogen bonds with PNP6Pglc but not PNP6Pgal. Lastly, using multiple 3D structures from various GH1 activities, a large network of conserved interactions between active site residues (and other important residues) was uncovered, which most likely stabilise the loop regions that contain these residues, helping to retain their positions needed for binding molecules. Alternatively, there exists several differing residue-residue interactions when comparing each of the activities which could contribute towards individual activity substrate specificity by causing slightly different overall structure and malleability of the active site. Altogether, the findings in this thesis shed light on the function, mechanisms, dynamics, and ligand-binding of GH1 enzymes – particularly of the 6Pβ-glycosidase activities. , 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
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
Analysis of bacterial Mur amide ligase enzymes for the identification of inhibitory compounds by in silico methods
- Authors: Chamboko, Chiratidzo Respina
- Date: 2020
- Subjects: Mur amide ligases , Ligases , Ligand binding (Biochemistry) , Antibacterial agents
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/163430 , vital:41036
- Description: An increased emergence of resistant pathogenic bacterial strains over the years has resulted in many people dying of untreatable infections. This has become one of the most critical global public health problems, as resistant strains are complicating treatment of infectious diseases, increasing human morbidity, mortality, and health care costs. A very limited amount of effective antibiotics is currently available, but the development of novel classes of antibacterial agents is becoming a priority. Mur amide ligases are enzymes that have been identified as potentially good targets for antibiotics, as they are uniquely found in bacteria. They are responsible for the formation of peptide bonds in a growing peptidoglycan structure for bacterial cell walls. The current work presented here focused on characterizing these Mur amide ligase enzymes and obtaining inhibitory compounds that could potentially be of use in drug discovery of antibacterial agents. To do this, multiple sequence alignment, motif analysis and phylogenetic tree constructions were carried out, followed by docking studies and molecular dynamic simulations. Prior to docking, homology modelling of missing residues in the MurF structure (PDB 1GG4) was performed. Characterization results revealed the Mur amide ligase enzymes contained defined conservation in limited regions, that ultimately mapped towards the central domain responsible for ATP binding (presence of a conserved GKT motif). Further analysis of results further unraveled the unique patterns observed within each group of the family of enzymes. As a result of these findings, docking studies were carried out on each Mur amide ligase structure. At most, two ligands were identified to be sufficiently inhibiting each Mur amide ligase. The ligands obtained were SANC00574 and SANC00575 for MurC, SANC00290 and SANC00438 for MurD, SANC00290 and SANC00525 for MurE and SANC00290 and SANC00434 for MurF. The two best ligands identified for each enzyme had docked in the active site of their respective proteins, passed Lipinski’s rule of five and had substantially low binding energies. Molecular dynamic simulations were then performed to analyze the behavior of the proteins and protein-ligand complexes, to confirm the lead compounds as good inhibitors of the Mur amide ligases. In the case of MurC, MurD and MurE complexes, the identified ligands clearly impacted the behavior of the protein, as the ligand bound proteins became more compact and stable, while flexibility decreased. There was however an opposite effect on MurF complexes, that resulted in identified inhibitors being discarded. As a potential next step, in vivo and in vitro experiments can be performed with identified ligands from this research, to further support the information presented.
- Full Text:
- Date Issued: 2020
The creation and validation of aptamers binding to murine 3T3-L1 Preadipocytes: preliminary implications for controlled cellular attachment, differentiation and cell fate
- Authors: Rubidge, Mark Lourens
- Date: 2017
- Subjects: Oligonucleotides , Fat cells , Stem cells , Ligand binding (Biochemistry) , Fluorimetry
- Language: English
- Type: Master's theses , text
- Identifier: http://hdl.handle.net/10962/65247 , vital:28714
- Description: The controlled seeding of a variety of stem cells in vitro has been reported to alter the patterns of their subsequent differentiation. This has been attributed to the control of the surface microenvironment onto which adherent stem cells are cultured, especially control of the proximal density of neighbouring cells. Simultaneously, advances in the generation of aptamers - synthetic ligand molecules developed using in vitro selection techniques targeting complex molecules - have aided in the production of molecules capable of selectively binding to a variety of commercial stem cell lines. Combining the aforementioned research fields, the project reported in this thesis aimed to generate DNA-based aptamers capable of assisting with the selective binding of murine 3T3- L1 preadipocytes to a solid surface. This was performed with a view to, eventually, control the seeding densities of the adherent preadipocytes on the surface of the tissue culture dish in subsequent researchers. In the process of meeting this goal, several optimisations of the in vitro process by which aptamers binding to cells are generated (Cell-SELEX) were performed: an analysis into a variety of methods used for the removal of the single stranded aptamer candidate sequences attached to the surface of 3T3-L1 preadipocytes, a comparison of methods for the generation of single-stranded aptamer sequences from double-stranded DNA template molecules and a method for quantifying the removed ssDNA from the cell surface. Their use is further reported in this work. Initially, it was determined that a fluorimetric evaluation of the unbound single stranded DNA was the optimum technique to use to evaluate the relative amounts of aptamer DNA binding to target cells during cell-SELEX; this arose from the release of DNA, and other cell lysate contaminates, which interfered UV/ Vis quantification. The evaluation into different methods of ssDNA removal from the cell surface showed that although trypsinisation of the cells demonstrated the highest level of aptamer detachment (quantified by fluorimetry), there is a decrease the number of potential targets that aptamers could attach to. The most common method for detaching bound DNA aptamer molecules from cellular targets reported in literature, the use of high temperatures, was selected for cell-SELEX to increase the variability in potential target sites on the cell surface. Using techniques optimised in this work, fluorescently-tagged single-stranded oligonucleotide aptamers were later generated with a positive selection pressure to bind to the surface of the 3T3-L1 preadipocytes, but not to their differentiated adipocyte counterparts. After eight cycles of cell-SELEX, fluorescent spectroscopic analysis depicted a 74 % binding retention of the selection pool in the positive preadipocyte selection pool, as opposed to a 0.69 % binding of sequences to the negative differentiated preadipocytes. Following the isolation and identification of candidate sequences, seven separate sequences were identified as being successfully generated from the selection process. Bioinformatic characterisation of these placed sequenced aptamer candidates into two separate families, that were then analysed in opposition to each for their binding affinity toward each other. Using fluorescently-tagged sequences, the binding selectivity of the generated aptamers was validated using both epifluorescent microscopy and confocal microscopy. At this stage, an aptamer sequence selected from prior in-house research to serve as a negative control also demonstrated significant binding to the extracellular matrix of both preadipocytes and mature adipocytes. 5’-thiolated aptamer sequences were used to form self-assembled monolayers on the electrode surfaces of the impedimetric Roche xCELLigence Real-Time Cell Analysis. The use of aptamer sequences to capture the seeded preadipocytes, demonstrated a slight increase in the extent of binding of the preadipocytes to the gold electrode surface and produced some preliminary indications of alterations to the pattern and rate of subsequent differentiation in the preadipocytes. This provides preliminary evidence that aptamers developed to bind specifically to a stem cell line in vitro show potential to be used as to capture said cell when cast in a self- assembled monolayer assembly. This provides a future opportunity to control the seeding densities of the cells in vitro. The effects of cellular differentiation at a set of predefined cellular densities can be demonstrated on a desired stem cell line. , Thesis (MSc) -- Faculty of Faculty of Science, Biotechnology Innovation Centre, 2017
- Full Text:
- Date Issued: 2017
In silico analysis of human Hsp90 for the identification of novel anti-cancer drug target sites and natural compound inhibitors
- Authors: Penkler, David Lawrence
- Date: 2015
- Subjects: Heat shock proteins , Cancer -- Treatment , Molecular chaperones , Homeostasis , Carcinogenesis , Chemotherapy , Ligand binding (Biochemistry) , Protein-protein interactions
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: vital:4162 , http://hdl.handle.net/10962/d1018938
- Description: The 90-KDa heat shock protein (Hsp90) is part of the molecular chaperone family, and as such it is involved in the regulation of protein homeostasis within cells. Specifically, Hsp90 aids in the folding of nascent proteins and re-folding of denatured proteins. It also plays an important role in the prevention of protein aggregation. Hsp90’s functionality is attributed to its several staged, multi-conformational ATPase cycle, in which associated client proteins are bound and released. Hsp90 is known to be associated with a wide array of client proteins, some of which are thought to be involved in multiple oncogenic processes. Indeed Hsp90 is known to be directly involved in perpetuating the stability and function of multiple mutated, chimeric and over-expressed signalling proteins that are known to promote the growth and survival of cancer cells. Hsp90 inhibitors are thus thought to be promising therapeutic agents for cancer treatment. A lack of a 3D structure of human Hsp90 however has restricted Hsp90 inhibitor development in large to in vivo investigations. This study, aims to investigate and calculate hypothetical homology models of the full human Hsp90 protein, and to probe these structural models for novel drug target sites using several in silico techniques. A multi-template homology modelling methodology was developed and in conjunction with protein-protein docking techniques, two functionally important human Hsp90 structural models were calculated; the nucleotide free “v-like” open and nucleotide bound closed conformations. Based on the conservation of ligand binding, virtual screening experiments conducted on both models using 316 natural compounds indigenous to South Africa, revealed three novel putative target sites. Two binding pockets in close association with important Hsp90-Hop interaction residues and a single binding pocket on the dimerization interface in the C-terminal domain. Targeted molecular docking experiments at these sites revealed two compounds (721395-11-5 and 264624-39-7) as putative inhibitors, both showing strong binding affinities for at least one of the three investigated target sites. Furthermore both compounds were found to only violate one Lipinski’s rules, suggesting their potential as candidates for further drug development. The combined work described here provides a putative platform for the development of next generation inhibitors of human Hsp90.
- Full Text:
- Date Issued: 2015