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
- 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
The characterization of GTP Cyclohydrolase I and 6-Pyruvoyl Tetrahydropterin Synthase enzymes as potential anti-malarial drug targets
- Khairallah, Afrah Yousif Huseein
- Authors: Khairallah, Afrah Yousif Huseein
- Date: 2022-04-08
- Subjects: Antimalarials , Plasmodium falciparum , Malaria Chemotherapy , Malaria Africa , Drug resistance , Drug development , Molecular dynamics
- Language: English
- Type: Doctoral thesis , text
- Identifier: http://hdl.handle.net/10962/233784 , vital:50127 , DOI 10.21504/10962/233784
- Description: Malaria remains a public health problem and a high burden of disease, especially in developing countries. The unicellular protozoan malaria parasite of the genus Plasmodium infects about a quarter of a billion people annually, with an estimated 409 000 death cases. The majority of malaria cases occurred in Africa; hence, the region is regarded as endemic for malaria. Global efforts to eradicate the disease led to a decrease in morbidity and mortality rates. However, an enormous burden of malaria infection remains, and it cannot go unnoticed. Countries with limited resources are more affected by the disease, mainly on its public health and socio-economic development, due to many factors besides malaria itself, such as lack of access to adequate, affordable treatments and preventative regimes. Furthermore, the current antimalarial drugs are losing their efficacy because of parasite drug resistance. The emerged drug resistance has reduced the drug efficacy in clearing the parasite from the host system, causing prolonged illness and a higher risk of death. Therefore, the emerged antimalarial drug resistance has hindered the global efforts for malaria control and elimination and established an urgent need for new treatment strategies. When the resistance against classical antimalarial drugs emerged, the class of antifolate antimalarial medicines became the most common alternative. The antifolate antimalarial drugs target the malaria parasite de novo folate biosynthesis pathway by limiting folate derivates, which are essential for the parasite cell growth and survival. Yet again, the malaria parasite developed resistance against the available antifolate drugs, rendering the drugs ineffective in many cases. Given the previous success in targeting the malaria parasite de novo folate biosynthesis pathway, alternative enzymes within this pathway stand as good targets and can be explored to develop new antifolate drugs with novel mechanisms of action. The primary focus of this thesis is to contribute to the existing and growing knowledge of antimalarial drug discovery. The study aims to characterise the malaria parasite de novo folate synthesis pathway enzymes guanosine-5'-triphosphate (GTP) cyclohydrolase I (GCH1) and 6-pyruvoyl tetrahydropterin synthase (PTPS) as alternative drug targets for malaria treatment by using computational approaches. Further, discover new allosteric drug targeting sites within the two enzymes' 3D structures for future drug design and discovery. Sequence and structural analysis were carried out to characterise and pinpoint the two enzymes' unique sequence and structure-based features. From the analyses, key sequence and structure differences were identified between the malaria parasite enzymes relative to their human homolog; the identified sites can aid significantly in designing and developing new antimalarial antifolate drugs with good selectivity toward the parasites’ enzymes. GCH1 and PTPS contain a catalytically essential metal ion in their active site; therefore, force field parameters were needed to study their active sites accurately during all-atom molecular dynamic simulations (MD). The force field parameters were derived through quantum mechanics potential energy surface scans of the metals bonded terms and evaluated via all-atom MD simulations. Proteins structural dynamics is imperative for many biological processes; thus, it is essential to consider the structural dynamics of proteins whilst understanding their function. In this regard, the normal mode analysis (NMA) approach based on the elastic network model (ENM) was employed to study the intrinsic dynamics and conformations changes of GCH1 and PTPS enzymes. The NMA disclosed essential structural information about the protein’s intrinsic dynamics and mechanism of allosteric modulation of their binding properties, further highlighting regions that govern their conformational changes. The analysis also disclosed hotspot residues that are crucial for the proteins' fold stability and function. The NMA was further combined with sequence motif results and showed that conserved residues of GCH1 and PTPS were located within the identified key structural sites modulating the proteins' conformational rearrangement. The characterized structural features and hotspot residues were regarded as potential allosteric sites of important value for the design and development of allosteric drugs. Both GCH1 and PTPS enzymes have never been targeted before and can provide an excellent opportunity to overcome the antimalarial antifolate drug resistance problem. The data presented in this thesis contribute to the understanding of the sequence, structure, and global dynamics of both GCH1 and PTPS, further disclose potential allosteric drug targeting sites and unique structural features of both enzymes that can establish a solid starting point for drug design and development of new antimalarial drugs of a novel mechanism of actions. Lastly, the reported force field parameters will be of value for MD simulations for future in-silico drug discovery studies involving the two enzymes and other enzymes with the same Zn2+ binding motifs and coordination environments. The impact of this research can facilitate the discovery of new effective antimalarial medicines with novel mechanisms of action. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-04-08
- Authors: Khairallah, Afrah Yousif Huseein
- Date: 2022-04-08
- Subjects: Antimalarials , Plasmodium falciparum , Malaria Chemotherapy , Malaria Africa , Drug resistance , Drug development , Molecular dynamics
- Language: English
- Type: Doctoral thesis , text
- Identifier: http://hdl.handle.net/10962/233784 , vital:50127 , DOI 10.21504/10962/233784
- Description: Malaria remains a public health problem and a high burden of disease, especially in developing countries. The unicellular protozoan malaria parasite of the genus Plasmodium infects about a quarter of a billion people annually, with an estimated 409 000 death cases. The majority of malaria cases occurred in Africa; hence, the region is regarded as endemic for malaria. Global efforts to eradicate the disease led to a decrease in morbidity and mortality rates. However, an enormous burden of malaria infection remains, and it cannot go unnoticed. Countries with limited resources are more affected by the disease, mainly on its public health and socio-economic development, due to many factors besides malaria itself, such as lack of access to adequate, affordable treatments and preventative regimes. Furthermore, the current antimalarial drugs are losing their efficacy because of parasite drug resistance. The emerged drug resistance has reduced the drug efficacy in clearing the parasite from the host system, causing prolonged illness and a higher risk of death. Therefore, the emerged antimalarial drug resistance has hindered the global efforts for malaria control and elimination and established an urgent need for new treatment strategies. When the resistance against classical antimalarial drugs emerged, the class of antifolate antimalarial medicines became the most common alternative. The antifolate antimalarial drugs target the malaria parasite de novo folate biosynthesis pathway by limiting folate derivates, which are essential for the parasite cell growth and survival. Yet again, the malaria parasite developed resistance against the available antifolate drugs, rendering the drugs ineffective in many cases. Given the previous success in targeting the malaria parasite de novo folate biosynthesis pathway, alternative enzymes within this pathway stand as good targets and can be explored to develop new antifolate drugs with novel mechanisms of action. The primary focus of this thesis is to contribute to the existing and growing knowledge of antimalarial drug discovery. The study aims to characterise the malaria parasite de novo folate synthesis pathway enzymes guanosine-5'-triphosphate (GTP) cyclohydrolase I (GCH1) and 6-pyruvoyl tetrahydropterin synthase (PTPS) as alternative drug targets for malaria treatment by using computational approaches. Further, discover new allosteric drug targeting sites within the two enzymes' 3D structures for future drug design and discovery. Sequence and structural analysis were carried out to characterise and pinpoint the two enzymes' unique sequence and structure-based features. From the analyses, key sequence and structure differences were identified between the malaria parasite enzymes relative to their human homolog; the identified sites can aid significantly in designing and developing new antimalarial antifolate drugs with good selectivity toward the parasites’ enzymes. GCH1 and PTPS contain a catalytically essential metal ion in their active site; therefore, force field parameters were needed to study their active sites accurately during all-atom molecular dynamic simulations (MD). The force field parameters were derived through quantum mechanics potential energy surface scans of the metals bonded terms and evaluated via all-atom MD simulations. Proteins structural dynamics is imperative for many biological processes; thus, it is essential to consider the structural dynamics of proteins whilst understanding their function. In this regard, the normal mode analysis (NMA) approach based on the elastic network model (ENM) was employed to study the intrinsic dynamics and conformations changes of GCH1 and PTPS enzymes. The NMA disclosed essential structural information about the protein’s intrinsic dynamics and mechanism of allosteric modulation of their binding properties, further highlighting regions that govern their conformational changes. The analysis also disclosed hotspot residues that are crucial for the proteins' fold stability and function. The NMA was further combined with sequence motif results and showed that conserved residues of GCH1 and PTPS were located within the identified key structural sites modulating the proteins' conformational rearrangement. The characterized structural features and hotspot residues were regarded as potential allosteric sites of important value for the design and development of allosteric drugs. Both GCH1 and PTPS enzymes have never been targeted before and can provide an excellent opportunity to overcome the antimalarial antifolate drug resistance problem. The data presented in this thesis contribute to the understanding of the sequence, structure, and global dynamics of both GCH1 and PTPS, further disclose potential allosteric drug targeting sites and unique structural features of both enzymes that can establish a solid starting point for drug design and development of new antimalarial drugs of a novel mechanism of actions. Lastly, the reported force field parameters will be of value for MD simulations for future in-silico drug discovery studies involving the two enzymes and other enzymes with the same Zn2+ binding motifs and coordination environments. The impact of this research can facilitate the discovery of new effective antimalarial medicines with novel mechanisms of action. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-04-08
- «
- ‹
- 1
- ›
- »