Generation of a virtual library of terpenes using graph theory, and its application in exploration of the mechanisms of terpene biosynthesis
- Authors: Dendera, Washington
- Date: 2020
- Subjects: Terpenes , Plants -- Metabolism , Computational biology , Bioinformatics , Organic compounds -- Synthesis , Monoterpenes , Molecular biology -- Computer simulation
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/123453 , vital:35439
- Description: Terpenes form a large group of organic compounds which have proven to be of use to many living organisms being used by plants for metabolism (Pichersky and Gershenzon, 1934; McGarvey and Croteau, 1995; Gershenzon and Dudareva, 2007), defence or as a means to attract pollinators and also used by humans in medical, pharmaceutical and food industry (Bicas, Dionísio and Pastore, 2009; Marmulla and Harder, 2014; Kandi et al., 2015). Following on literature methods to generate chemical libraries using graph theoretic techniques, complete libraries of all possible terpene isomers have been constructed with the goal of construction of derivative libraries of possible carbocation intermediates which are important in the elucidation of mechanisms in the biosynthesis of terpenes. Virtual library generation of monoterpenes was first achieved by generating graphs of order 7, 8, 9 and 10 using the Nauty and Traces suite. These were screened and processed with a set of collated Python scripts written to recognize the graphs in text format and translate them to molecules, minimizing through Tinker whilst discarding graphs that violate chemistry laws. As a result of the computational time required only order 7 and order 10 graphs were processed. Out of the 873 graphs generated from order seven, 353 were converted to molecules and from the 11,7 million produced from order 10 half were processed resulting in the production of 442928 compounds (repeats included). For screening, 55 366 compounds were docked in the active site of limonene synthase; of these 2355 ligands had a good Vina docking score with a binding energy of between -7.0 and -7.4 kcal.mol-1. When these best docked molecules were overlaid in the active site a map of possible ligand positions within the active site of limonene synthase was traced out.
- Full Text:
- Date Issued: 2020
- Authors: Dendera, Washington
- Date: 2020
- Subjects: Terpenes , Plants -- Metabolism , Computational biology , Bioinformatics , Organic compounds -- Synthesis , Monoterpenes , Molecular biology -- Computer simulation
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/123453 , vital:35439
- Description: Terpenes form a large group of organic compounds which have proven to be of use to many living organisms being used by plants for metabolism (Pichersky and Gershenzon, 1934; McGarvey and Croteau, 1995; Gershenzon and Dudareva, 2007), defence or as a means to attract pollinators and also used by humans in medical, pharmaceutical and food industry (Bicas, Dionísio and Pastore, 2009; Marmulla and Harder, 2014; Kandi et al., 2015). Following on literature methods to generate chemical libraries using graph theoretic techniques, complete libraries of all possible terpene isomers have been constructed with the goal of construction of derivative libraries of possible carbocation intermediates which are important in the elucidation of mechanisms in the biosynthesis of terpenes. Virtual library generation of monoterpenes was first achieved by generating graphs of order 7, 8, 9 and 10 using the Nauty and Traces suite. These were screened and processed with a set of collated Python scripts written to recognize the graphs in text format and translate them to molecules, minimizing through Tinker whilst discarding graphs that violate chemistry laws. As a result of the computational time required only order 7 and order 10 graphs were processed. Out of the 873 graphs generated from order seven, 353 were converted to molecules and from the 11,7 million produced from order 10 half were processed resulting in the production of 442928 compounds (repeats included). For screening, 55 366 compounds were docked in the active site of limonene synthase; of these 2355 ligands had a good Vina docking score with a binding energy of between -7.0 and -7.4 kcal.mol-1. When these best docked molecules were overlaid in the active site a map of possible ligand positions within the active site of limonene synthase was traced out.
- Full Text:
- Date Issued: 2020
Prediction of mass spectra for natural products using an ab initio approach
- Authors: Novokoza, Yolanda
- Date: 2020
- Subjects: Molecular dynamics , Molecular dynamics -- Computer simulation , Mass spectroscopy , Electron impact ionization
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/167166 , vital:41443
- Description: Mass spectrometry (MS) is a technique that measures the fragmentation of molecules, dependent on the molecule’s chemical composition and structure, by first introducing a charge on the molecules. The instrument records the mass to charge ratio, but the energy from the ionization process causes the molecule to fragment. The resultant mass spectrum is highly indicative of not only the molecule analyzed, but also its chemical composition. MS is used in research and industry for both routine and research purposes. One such way to ionize molecules for MS is by bombarding the molecule with electrons which is the basis of electron impact mass spectrometry (EIMS). Although EIMS is widely used, prediction of electron impact mass spectra from first principles is a challenging problem due to a need to accurately determine the probability of different fragmentation pathways of a molecule. Ab initio molecular dynamics based methods are able to explore in an automatic fashion the energetically available fragmentation paths thus give reaction mechanisms in an unbiased way. The mass spectra of five molecules have been explored in work-flows leading to the prediction of mass spectra. These molecules include three natural products alpha-hispanolol, PFB oxime derivative and boronolide (for which experimental mass spectra were not available) and two compounds from the NIST database (for which experimental mass spectra were available). For each of these systems many random conformations were generated using the RDKit library. To all conformations random velocities were applied to each atom. Ab initio molecular dynamics was performed on each conformer, using these initial random velocities using CP2K software, at DFTB+ level at a variety of highly raised temperatures (to accelerate the formation of fragments) Fragmentation was monitored by iterating through all bonds, and identifying bond breakages during dynamics. Graph theoretical packages were used then to track distinct fragments generated. For each of these fragments, charges were determined from Mulliken analysis for all atoms on the fragment from the QM calculations and sum of atomic spin densities per fragment was also plotted. The fragment with the greatest charge (corresponding to the formation of a cation fragment) was taken for plotting on the mass spectrum. Finally, from the mass of the fragment and its elemental composition, the isotopic distribution for the fragment was determined, and this distribution was included by addition in to the mass spectrum. For all trajectories, the sum of all isotopic distributions determined the final mass spectrum.
- Full Text:
- Date Issued: 2020
- Authors: Novokoza, Yolanda
- Date: 2020
- Subjects: Molecular dynamics , Molecular dynamics -- Computer simulation , Mass spectroscopy , Electron impact ionization
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/167166 , vital:41443
- Description: Mass spectrometry (MS) is a technique that measures the fragmentation of molecules, dependent on the molecule’s chemical composition and structure, by first introducing a charge on the molecules. The instrument records the mass to charge ratio, but the energy from the ionization process causes the molecule to fragment. The resultant mass spectrum is highly indicative of not only the molecule analyzed, but also its chemical composition. MS is used in research and industry for both routine and research purposes. One such way to ionize molecules for MS is by bombarding the molecule with electrons which is the basis of electron impact mass spectrometry (EIMS). Although EIMS is widely used, prediction of electron impact mass spectra from first principles is a challenging problem due to a need to accurately determine the probability of different fragmentation pathways of a molecule. Ab initio molecular dynamics based methods are able to explore in an automatic fashion the energetically available fragmentation paths thus give reaction mechanisms in an unbiased way. The mass spectra of five molecules have been explored in work-flows leading to the prediction of mass spectra. These molecules include three natural products alpha-hispanolol, PFB oxime derivative and boronolide (for which experimental mass spectra were not available) and two compounds from the NIST database (for which experimental mass spectra were available). For each of these systems many random conformations were generated using the RDKit library. To all conformations random velocities were applied to each atom. Ab initio molecular dynamics was performed on each conformer, using these initial random velocities using CP2K software, at DFTB+ level at a variety of highly raised temperatures (to accelerate the formation of fragments) Fragmentation was monitored by iterating through all bonds, and identifying bond breakages during dynamics. Graph theoretical packages were used then to track distinct fragments generated. For each of these fragments, charges were determined from Mulliken analysis for all atoms on the fragment from the QM calculations and sum of atomic spin densities per fragment was also plotted. The fragment with the greatest charge (corresponding to the formation of a cation fragment) was taken for plotting on the mass spectrum. Finally, from the mass of the fragment and its elemental composition, the isotopic distribution for the fragment was determined, and this distribution was included by addition in to the mass spectrum. For all trajectories, the sum of all isotopic distributions determined the final mass spectrum.
- Full Text:
- Date Issued: 2020
In silico study of Plasmodium 1-deoxy-dxylulose 5-phosphate reductoisomerase (DXR) for identification of novel inhibitors from SANCDB
- Authors: Diallo, Bakary N'tji
- Date: 2018
- Subjects: Plasmodium 1-deoxy-dxylulose 5-phosphate reductoisomerase , Isoprenoids , Plasmodium , Antimalarials , Malaria -- Chemotherapy , Molecules -- Models , Molecular dynamics , South African Natural Compounds Database
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/64012 , vital:28523
- Description: Malaria remains a major health concern with a complex parasite constantly developing resistance to the different drugs introduced to treat it, threatening the efficacy of the current ACT treatment recommended by WHO (World Health Organization). Different antimalarial compounds with different mechanisms of action are ideal as this decreases chances of resistance occurring. Inhibiting DXR and consequently the MEP pathway is a good strategy to find a new antimalarial with a novel mode of action. From literature, all the enzymes of the MEP pathway have also been shown to be indispensable for the synthesis of isoprenoids. They have been validated as drug targets and the X-ray structure of each of the enzymes has been solved. DXR is a protein which catalyses the second step of the MEP pathway. There are currently 255 DXR inhibitors in the Binding Database (accessed November 2017) generally based on the fosmidomycin structural scaffold and thus often showing poor drug likeness properties. This study aims to research new DXR inhibitors using in silico techniques. We analysed the protein sequence and built 3D models in close and open conformations for the different Plasmodium sequences. Then SANCDB compounds were screened to identify new potential DXR inhibitors with new chemical scaffolds. Finally, the identified hits were submitted to molecular dynamics studies, preceded by a parameterization of the manganese atom in the protein active site.
- Full Text:
- Date Issued: 2018
- Authors: Diallo, Bakary N'tji
- Date: 2018
- Subjects: Plasmodium 1-deoxy-dxylulose 5-phosphate reductoisomerase , Isoprenoids , Plasmodium , Antimalarials , Malaria -- Chemotherapy , Molecules -- Models , Molecular dynamics , South African Natural Compounds Database
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/64012 , vital:28523
- Description: Malaria remains a major health concern with a complex parasite constantly developing resistance to the different drugs introduced to treat it, threatening the efficacy of the current ACT treatment recommended by WHO (World Health Organization). Different antimalarial compounds with different mechanisms of action are ideal as this decreases chances of resistance occurring. Inhibiting DXR and consequently the MEP pathway is a good strategy to find a new antimalarial with a novel mode of action. From literature, all the enzymes of the MEP pathway have also been shown to be indispensable for the synthesis of isoprenoids. They have been validated as drug targets and the X-ray structure of each of the enzymes has been solved. DXR is a protein which catalyses the second step of the MEP pathway. There are currently 255 DXR inhibitors in the Binding Database (accessed November 2017) generally based on the fosmidomycin structural scaffold and thus often showing poor drug likeness properties. This study aims to research new DXR inhibitors using in silico techniques. We analysed the protein sequence and built 3D models in close and open conformations for the different Plasmodium sequences. Then SANCDB compounds were screened to identify new potential DXR inhibitors with new chemical scaffolds. Finally, the identified hits were submitted to molecular dynamics studies, preceded by a parameterization of the manganese atom in the protein active site.
- Full Text:
- Date Issued: 2018
Investigating the influence of ring substitution on indole hydrogen bonding, with amino acids
- Authors: Nel, Donovan
- Date: 2018
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/63509 , vital:28426
- Description: Expected release date-April 2019
- Full Text: false
- Date Issued: 2018
- Authors: Nel, Donovan
- Date: 2018
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/63509 , vital:28426
- Description: Expected release date-April 2019
- Full Text: false
- Date Issued: 2018
In silico analysis of plasmodium falciparum Hsp70-x for potential binding sites and hits
- Authors: Amusengeri, Arnold
- Date: 2017
- Subjects: Uncatalogued
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/59136 , vital:27435
- Description: Restricted access-thesis embargoed for 1 year - release date April 2019
- Full Text:
- Date Issued: 2017
- Authors: Amusengeri, Arnold
- Date: 2017
- Subjects: Uncatalogued
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/59136 , vital:27435
- Description: Restricted access-thesis embargoed for 1 year - release date April 2019
- Full Text:
- Date Issued: 2017
- «
- ‹
- 1
- ›
- »