Computational search for nature-derived dual-action inhibitors of HIV-1 reverse transcriptase and integrase: a potential strategy to mitigate drug resistance progression
- Authors: Mwiinga, Luyando
- Date: 2024-10-11
- Subjects: HIV (Viruses) , Reverse transcriptase , Antiretroviral agents , RDKit , Drug resistance , Docking
- Language: English
- Type: Academic theses , Master's theses , text
- Identifier: http://hdl.handle.net/10962/463930 , vital:76458
- Description: Human immunodeficiency virus Type 1 (HIV-1) is a devastating viral infection affecting millions worldwide and presents significant challenges in treatment and management. In 2022, approximately 39 million people were living with HIV with Sub-Saharan Africa having two thirds of these infections. Devastatingly, there were approximately 300 000 HIV/AIDS related deaths in Sub-Saharan Africa alone in 2022 alone. Antiretroviral therapy (ART) which is fundamental for HIV treatment, comprises of a combination of drugs such as nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTs), protease inhibitors (PIs) and integrase strand transfer inhibitors (INSTIs). However, although 28.7 million people out of the estimated 38.4 million people living with HIV in 2021 were receiving ART, the emergence of drug-resistant strains further complicates treatment efforts, highlighting the need for novel therapeutic approaches. This study aimed to address the challenges raised by drug resistance and significant side effects by identifying potential dual inhibitors against HIV-1 Reverse Transcriptase (RT) and Integrase (IN) using in silico techniques. RT RNase H and IN were chosen as targets for their shared dependency on Mg2+ ions within their active sites, which are crucial for catalytic activity. The selection of dual inhibitors was motivated by the fact that the virus would need to replicate at two points simultaneously to develop resistance, making it less likely. The objectives of this study included the creation of a natural derivative compound library using RDKit with the aid of SciFinder, utilizing (-)-epigallocatechin-3-O-gallate (EGCG), because of its dual inhibitory effects against RT and IN, as indicated by a study conducted by Sanna et al. 2019. The natural derivatives were chosen to take advantage of their chemical diversity and to explore potential novel therapeutic options for combating HIV drug resistance. The compound library created comprised of 125 203 compounds. Then docking studies were conducted to assess proteinligand binding. After the correlation of the RT and IN docking studies, 288 compounds were filtered to have potential dual inhibitory activity. Then quantitative estimation of druggability (QED) analysis identified three compounds with superior properties compared to EGCG and FDAapproved drug raltegravir (RAL). Molecular docking simulations revealed interactions between the inhibitors and the key active site residues of RT and IN, along with the chelation of at least one 3 Mg2+, suggesting the potential for enzymatic disruption. Furthermore, molecular dynamic (MD) simulations were then conducted to assess protein-ligand system behavior, through RMSD and RMSF analysis. The RMSD analysis uncovered instability in the IN-Sci30703 complex, leading to its exclusion as a potential dual action inhibitor. RMSF analysis for IN showed that all the inhibitors had the ability to limit the flexibility of the catalytic loop which is essential for catalytic activity. Therefore, further in vitro studies are required to evaluate the effectiveness of the remaining two EGCG derivatives (Sci33211 and Sci48919) in inhibiting RT and IN through the chelation of at least one Mg2+ ion to determine if they have superior dual inhibitory effects compared to EGCG. This study adds to the ongoing efforts to develop effective strategies against HIV-1 drug resistance and emphasizes the importance of continued research in this field. , Thesis (MSc) -- Faculty of Science, Biochemistry, Microbiology & Bioinformatics, 2024
- Full Text:
- Date Issued: 2024-10-11
- Authors: Mwiinga, Luyando
- Date: 2024-10-11
- Subjects: HIV (Viruses) , Reverse transcriptase , Antiretroviral agents , RDKit , Drug resistance , Docking
- Language: English
- Type: Academic theses , Master's theses , text
- Identifier: http://hdl.handle.net/10962/463930 , vital:76458
- Description: Human immunodeficiency virus Type 1 (HIV-1) is a devastating viral infection affecting millions worldwide and presents significant challenges in treatment and management. In 2022, approximately 39 million people were living with HIV with Sub-Saharan Africa having two thirds of these infections. Devastatingly, there were approximately 300 000 HIV/AIDS related deaths in Sub-Saharan Africa alone in 2022 alone. Antiretroviral therapy (ART) which is fundamental for HIV treatment, comprises of a combination of drugs such as nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTs), protease inhibitors (PIs) and integrase strand transfer inhibitors (INSTIs). However, although 28.7 million people out of the estimated 38.4 million people living with HIV in 2021 were receiving ART, the emergence of drug-resistant strains further complicates treatment efforts, highlighting the need for novel therapeutic approaches. This study aimed to address the challenges raised by drug resistance and significant side effects by identifying potential dual inhibitors against HIV-1 Reverse Transcriptase (RT) and Integrase (IN) using in silico techniques. RT RNase H and IN were chosen as targets for their shared dependency on Mg2+ ions within their active sites, which are crucial for catalytic activity. The selection of dual inhibitors was motivated by the fact that the virus would need to replicate at two points simultaneously to develop resistance, making it less likely. The objectives of this study included the creation of a natural derivative compound library using RDKit with the aid of SciFinder, utilizing (-)-epigallocatechin-3-O-gallate (EGCG), because of its dual inhibitory effects against RT and IN, as indicated by a study conducted by Sanna et al. 2019. The natural derivatives were chosen to take advantage of their chemical diversity and to explore potential novel therapeutic options for combating HIV drug resistance. The compound library created comprised of 125 203 compounds. Then docking studies were conducted to assess proteinligand binding. After the correlation of the RT and IN docking studies, 288 compounds were filtered to have potential dual inhibitory activity. Then quantitative estimation of druggability (QED) analysis identified three compounds with superior properties compared to EGCG and FDAapproved drug raltegravir (RAL). Molecular docking simulations revealed interactions between the inhibitors and the key active site residues of RT and IN, along with the chelation of at least one 3 Mg2+, suggesting the potential for enzymatic disruption. Furthermore, molecular dynamic (MD) simulations were then conducted to assess protein-ligand system behavior, through RMSD and RMSF analysis. The RMSD analysis uncovered instability in the IN-Sci30703 complex, leading to its exclusion as a potential dual action inhibitor. RMSF analysis for IN showed that all the inhibitors had the ability to limit the flexibility of the catalytic loop which is essential for catalytic activity. Therefore, further in vitro studies are required to evaluate the effectiveness of the remaining two EGCG derivatives (Sci33211 and Sci48919) in inhibiting RT and IN through the chelation of at least one Mg2+ ion to determine if they have superior dual inhibitory effects compared to EGCG. This study adds to the ongoing efforts to develop effective strategies against HIV-1 drug resistance and emphasizes the importance of continued research in this field. , Thesis (MSc) -- Faculty of Science, Biochemistry, Microbiology & Bioinformatics, 2024
- Full Text:
- Date Issued: 2024-10-11
Towards a possible future solution against Multidrug Resistance: An in silico exploration of the Multidrug and Toxic compound Extrusion (MATE) transporter proteins as potential antimicrobial drug targets
- Authors: Damji, Amira Mahamood
- Date: 2024-04-04
- Subjects: Multidrug resistance , Multidrug and toxic compound extrusion family, eukaryotic , Docking , Molecular dynamics , Drug development , Transmembrane protein
- Language: English
- Type: Academic theses , Master's theses , text
- Identifier: http://hdl.handle.net/10962/435009 , vital:73123
- Description: The rise of multidrug resistance (MDR) has become a pressing global issue, hindering the treatment of cancers and infectious diseases, and imposing a burden on healthcare systems and the economy. The Multidrug and Toxic compound Extrusion (MATE) superfamily of membrane efflux transporters is one of the key players contributing to MDR due to their ability to export a wide range of cationic and hydrophilic xenobiotics, including treatment drugs, from cells, diminishing their efficacy. Targeting MATE transporters holds great promise in achieving some cellular control over MDR, but first, a deeper understanding of their structure-function-dynamics link is required. This study aimed to explore the MATE transporters as potential antimicrobial drug targets using a two-fold in silico approach. First, virtual screening of compounds from the South African Natural Compounds Database (SANCDB) was performed to identify prospective lead inhibitory compounds against the MATE transporters using molecular docking, and top hits were selected based on their binding energy and interaction with the active site on the N-lobe of the protein. Second, to investigate the molecular-level dynamics of their extrusion mechanism, the MATE transporter structures were embedded in a POPC membrane bilayer using the CHARMM-GUI online tool and then subjected to MD simulations for 100 ns with the CHARMM 36m force field using GROMACS. The resulting trajectories were evaluated using three standard metrics – RMSD, RMSF, and Rg; significant global structural changes were observed and key functional regions in both membrane- and non-membrane transmembrane (TM) segments were identified, containing more dynamic and flexible residues than other regions. Furthermore, the MATE transporters showed more of a loosely-packed structure, providing flexibility to allow for conformational switching during their substrate-transport cycle, which is typical for proteins whose secondary structures are composed of all α-helices. The scope of this study lied in the preliminary stages of the computer-aided drug design process, and provided insights that can be used to guide the development of strategies aimed at regulating or inhibiting the function of the MATE transporters, offering a possible future solution to the growing challenge of MDR. , Thesis (MSc) -- Faculty of Science, Biochemistry and Microbiology, 2024
- Full Text:
- Date Issued: 2024-04-04
- Authors: Damji, Amira Mahamood
- Date: 2024-04-04
- Subjects: Multidrug resistance , Multidrug and toxic compound extrusion family, eukaryotic , Docking , Molecular dynamics , Drug development , Transmembrane protein
- Language: English
- Type: Academic theses , Master's theses , text
- Identifier: http://hdl.handle.net/10962/435009 , vital:73123
- Description: The rise of multidrug resistance (MDR) has become a pressing global issue, hindering the treatment of cancers and infectious diseases, and imposing a burden on healthcare systems and the economy. The Multidrug and Toxic compound Extrusion (MATE) superfamily of membrane efflux transporters is one of the key players contributing to MDR due to their ability to export a wide range of cationic and hydrophilic xenobiotics, including treatment drugs, from cells, diminishing their efficacy. Targeting MATE transporters holds great promise in achieving some cellular control over MDR, but first, a deeper understanding of their structure-function-dynamics link is required. This study aimed to explore the MATE transporters as potential antimicrobial drug targets using a two-fold in silico approach. First, virtual screening of compounds from the South African Natural Compounds Database (SANCDB) was performed to identify prospective lead inhibitory compounds against the MATE transporters using molecular docking, and top hits were selected based on their binding energy and interaction with the active site on the N-lobe of the protein. Second, to investigate the molecular-level dynamics of their extrusion mechanism, the MATE transporter structures were embedded in a POPC membrane bilayer using the CHARMM-GUI online tool and then subjected to MD simulations for 100 ns with the CHARMM 36m force field using GROMACS. The resulting trajectories were evaluated using three standard metrics – RMSD, RMSF, and Rg; significant global structural changes were observed and key functional regions in both membrane- and non-membrane transmembrane (TM) segments were identified, containing more dynamic and flexible residues than other regions. Furthermore, the MATE transporters showed more of a loosely-packed structure, providing flexibility to allow for conformational switching during their substrate-transport cycle, which is typical for proteins whose secondary structures are composed of all α-helices. The scope of this study lied in the preliminary stages of the computer-aided drug design process, and provided insights that can be used to guide the development of strategies aimed at regulating or inhibiting the function of the MATE transporters, offering a possible future solution to the growing challenge of MDR. , Thesis (MSc) -- Faculty of Science, Biochemistry and Microbiology, 2024
- Full Text:
- Date Issued: 2024-04-04
An in-silico study of the type II NADH: Quinone Oxidoreductase (ndh2). A new anti-malaria drug target
- Authors: Baye, Bertha Cinthia
- Date: 2022-10-14
- Subjects: Malaria , Plasmodium , Molecular dynamics , Computer simulation , Quinone , Antimalarials , Molecules Models , Docking , Drugs Computer-aided design
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/365633 , vital:65767 , DOI https://doi.org/10.21504/10962/365633
- Description: Malaria is caused by Plasmodium parasites, spread to people through the bites of infected female Anopheles mosquitoes. This study focuses on all 5 (Plasmodium falciparum, Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale and Plasmodium vivax) parasites that cause malaria in humans. Africa is a developing continent, and it is the most affected with an estimation of 90% of more than 400 000 malaria-related deaths reported by the World Health Organization (WHO) report in 2020, in which 61% of that number are children under the ages of five. Malaria resistance was initially observed in early 1986 and with the progression of time anti-malarial drug resistance has only increased. As a result, there is a need to study the malarial proteins mechanism of action and identify alternative treatment strategies for this disease. Type II NADH: quinone oxidoreductase (NDH2) is a monotopic protein that catalyses the electron transfer from NADH to quinone via FAD without a proton-pumping activity, and functions as an initial enzyme, either in addition to or as an alternative to proton-pumping NADH dehydrogenase (complex I) in the respiratory chain of bacteria, archaea, and fungal and plant mitochondrial. The structures for the Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale and Plasmodium vivax were modelled from the crystal structure of Plasmodium falciparum (5JWA). Compounds from the South African natural compounds database (SANCDB) were docked against both the NDH2 crystal structure and modelled structures. By performing in silico screening the study aimed to find potential compounds that might interrupt the electron transfer to quinone therefore disturbing the enzyme‟s function and thereby possibly eliminating the plasmodium parasite. CHARMM-GUI was used to create the membrane (since this work is with membrane-bound proteins) and to orient the protein on the membrane using OPM server guidelines, the interface produced GROMACS topology files that were used in molecular dynamics simulations. Molecular dynamics simulations were performed in the Centre for high performance computing (CHPC) cluster under the CHEM0802 project and the trajectories produced were further analysed. In this work not only were hit compounds from SANCDB identified, but also differences in behaviour across species and in the presence or absence of the membrane were described. This highlights the need to include the correct protein environment when studying these systems. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-10-14
- Authors: Baye, Bertha Cinthia
- Date: 2022-10-14
- Subjects: Malaria , Plasmodium , Molecular dynamics , Computer simulation , Quinone , Antimalarials , Molecules Models , Docking , Drugs Computer-aided design
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/365633 , vital:65767 , DOI https://doi.org/10.21504/10962/365633
- Description: Malaria is caused by Plasmodium parasites, spread to people through the bites of infected female Anopheles mosquitoes. This study focuses on all 5 (Plasmodium falciparum, Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale and Plasmodium vivax) parasites that cause malaria in humans. Africa is a developing continent, and it is the most affected with an estimation of 90% of more than 400 000 malaria-related deaths reported by the World Health Organization (WHO) report in 2020, in which 61% of that number are children under the ages of five. Malaria resistance was initially observed in early 1986 and with the progression of time anti-malarial drug resistance has only increased. As a result, there is a need to study the malarial proteins mechanism of action and identify alternative treatment strategies for this disease. Type II NADH: quinone oxidoreductase (NDH2) is a monotopic protein that catalyses the electron transfer from NADH to quinone via FAD without a proton-pumping activity, and functions as an initial enzyme, either in addition to or as an alternative to proton-pumping NADH dehydrogenase (complex I) in the respiratory chain of bacteria, archaea, and fungal and plant mitochondrial. The structures for the Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale and Plasmodium vivax were modelled from the crystal structure of Plasmodium falciparum (5JWA). Compounds from the South African natural compounds database (SANCDB) were docked against both the NDH2 crystal structure and modelled structures. By performing in silico screening the study aimed to find potential compounds that might interrupt the electron transfer to quinone therefore disturbing the enzyme‟s function and thereby possibly eliminating the plasmodium parasite. CHARMM-GUI was used to create the membrane (since this work is with membrane-bound proteins) and to orient the protein on the membrane using OPM server guidelines, the interface produced GROMACS topology files that were used in molecular dynamics simulations. Molecular dynamics simulations were performed in the Centre for high performance computing (CHPC) cluster under the CHEM0802 project and the trajectories produced were further analysed. In this work not only were hit compounds from SANCDB identified, but also differences in behaviour across species and in the presence or absence of the membrane were described. This highlights the need to include the correct protein environment when studying these systems. , Thesis (PhD) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-10-14
In silico substrate binding profiling for SARS-COV-2 main protease (mpro) using hexapeptide substrates
- Authors: Zabo, Sophakama
- Date: 2022-10-14
- Subjects: COVID-19 (Disease) , Peptides , Chymotrypsin like , Chymotrypsin , Proteases , Proteolytic enzymes
- Language: English
- Type: Academic theses , Master's theses , text
- Identifier: http://hdl.handle.net/10962/365566 , vital:65760
- Description: COVID-19, as a disease resulting from SARS-CoV-2 infection, and a pandemic has had a devastating effect on the world. There are limited effective measures that control the spread and treatment of COVID-19 illness. The homodimeric cysteine main protease (Mpro) is crucial to the life cycle of the virus, as it cleaves the large polyproteins 1a and 1ab into matured, functional non-structural proteins. The Mpro exhibits high degrees of conservation in sequence, structure and specificity across coronavirus species, making it an ideal drug target. The Mpro substrate-binding profiles remain, despite the resolution of its recognition sequence and cleavage points (Leu-Gln↓(Ser/Ala/Gly)). In this study, a series of hexapeptide sequences containing the appropriate recognition sequence and cleavage points were generated and screened against the Mpro to study these binding profiles, and to further be the basis for efficiency-driven drug design. A multi-conformer hexapeptide substrate library comprising optimised 81000 models of 810 unique sequences was generated using RDKit within the context of python. Terminal capping with ACE and NMe was effected using SMILES and SMARTS matching. Multiple hexapeptides were complexed with chain B of crystallographic Mpro (PDS ID: 6XHM), following the validation of chain B for this purpose using AutoDock Vina at high levels of exhaustiveness (480). The resulting Vina scores ranged between -8.7 and -7.0 kcal.mol-1, and the reproducibility of best poses was validated through redocking. Ligand efficiency indices were calculated to identify substrate residues with high binding efficiency at their respective positions, revealing Val (P3), Ala (P1′); and Gly and Ala (P2′ and P3′) as leading efficient binders. Binding efficiencies were lowered by molecular weight. Substrate recognition was assessed by mapping of binding subsites, and Mpro specificity was evaluated through the resolution of intermolecular interaction at the binding interface. Molecular dynamics simulations for 20 ns were performed to assess the stability and behaviour of 132 Mpro systems complexed with KLQ*** substrates. Principal component analysis (PCA), was performed to assess II protein motions and conformational changes during the simulations. A strategy was formulated to classify and evaluate relations in the Mpro PCA motions, revealing four main clades of similarity. Similarity within a clade (Group 2) and dissimilarity between clades were confirmed. Trajectory visualisation revealed complex stability, substrate unbinding and dimer dissociation for various Mpro systems. , Thesis (MSc) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-10-14
- Authors: Zabo, Sophakama
- Date: 2022-10-14
- Subjects: COVID-19 (Disease) , Peptides , Chymotrypsin like , Chymotrypsin , Proteases , Proteolytic enzymes
- Language: English
- Type: Academic theses , Master's theses , text
- Identifier: http://hdl.handle.net/10962/365566 , vital:65760
- Description: COVID-19, as a disease resulting from SARS-CoV-2 infection, and a pandemic has had a devastating effect on the world. There are limited effective measures that control the spread and treatment of COVID-19 illness. The homodimeric cysteine main protease (Mpro) is crucial to the life cycle of the virus, as it cleaves the large polyproteins 1a and 1ab into matured, functional non-structural proteins. The Mpro exhibits high degrees of conservation in sequence, structure and specificity across coronavirus species, making it an ideal drug target. The Mpro substrate-binding profiles remain, despite the resolution of its recognition sequence and cleavage points (Leu-Gln↓(Ser/Ala/Gly)). In this study, a series of hexapeptide sequences containing the appropriate recognition sequence and cleavage points were generated and screened against the Mpro to study these binding profiles, and to further be the basis for efficiency-driven drug design. A multi-conformer hexapeptide substrate library comprising optimised 81000 models of 810 unique sequences was generated using RDKit within the context of python. Terminal capping with ACE and NMe was effected using SMILES and SMARTS matching. Multiple hexapeptides were complexed with chain B of crystallographic Mpro (PDS ID: 6XHM), following the validation of chain B for this purpose using AutoDock Vina at high levels of exhaustiveness (480). The resulting Vina scores ranged between -8.7 and -7.0 kcal.mol-1, and the reproducibility of best poses was validated through redocking. Ligand efficiency indices were calculated to identify substrate residues with high binding efficiency at their respective positions, revealing Val (P3), Ala (P1′); and Gly and Ala (P2′ and P3′) as leading efficient binders. Binding efficiencies were lowered by molecular weight. Substrate recognition was assessed by mapping of binding subsites, and Mpro specificity was evaluated through the resolution of intermolecular interaction at the binding interface. Molecular dynamics simulations for 20 ns were performed to assess the stability and behaviour of 132 Mpro systems complexed with KLQ*** substrates. Principal component analysis (PCA), was performed to assess II protein motions and conformational changes during the simulations. A strategy was formulated to classify and evaluate relations in the Mpro PCA motions, revealing four main clades of similarity. Similarity within a clade (Group 2) and dissimilarity between clades were confirmed. Trajectory visualisation revealed complex stability, substrate unbinding and dimer dissociation for various Mpro systems. , Thesis (MSc) -- Faculty of Science, Biochemistry and Microbiology, 2022
- Full Text:
- Date Issued: 2022-10-14
Large and multi scale mechanistic modeling of Diels-Alder reactions
- Authors: Isamura, Bienfait Kabuyaya
- Date: 2022-04-06
- Subjects: Computational chemistry , Diels-Alder reaction , Python (Computer program language) , Reaction force theory , Fullerenes , Diolefins , AMADAR (Automated workflow for Mechanistic Analysis of Diels-Alder Reactions , ONIOM
- Language: English
- Type: Master's thesis , text
- Identifier: http://hdl.handle.net/10962/232317 , vital:49981
- Description: The [4+2] cycloaddition reaction between conjugated dienes and substituted alkenes is known as the Diels-Alder (DA) reaction, in honor of two German chemists, Otto Diels and Kurt Alder, who first reported this marvelous chemical transformation. The DA reaction is one of the most popular reactions in organic chemistry, allowing for the regio- and stereospecific establishment of six-membered rings with up to four stereogenic centers. This pericyclic reaction has found many applications in areas as diverse as natural products chemistry, polymer chemistry, and agrochemistry. Over the past decades, the mechanism of the Diels-Alder (DA) reaction has been the subject of numerous studies, dealing with questions as diverse as the mechanistic pathway, the synchronicity, the use of catalysts, the effect of solvents and salts, etc. On the other hand, as an example, fullerenes (and particularly [60] fullerene) have been found to act as good dienophiles in DA reactions to the extent that many functionalized fullerenes with interesting applications are still synthesized by reacting C60 with dienes. However, despite the very abundant literature about the mechanism of the DA reaction, some pertinent questions have been still pending, including, without being restricted to, the prediction of transition state (TS) geometries and the modeling of DA reactions involving large systems, such as those of C60 fullerene. It must be emphasized that TSs are not easy to predict and the main reason is that many existing algorithms require that the search is initiated from a good starting point (guess TS), which must be very similar to the actual TS. This problem is even more difficult when many TSs are to be located as may be the case in large-scale studies. Moreover, due to the large size of the C60 molecule, the usage of accurate high-level computational methods in the investigation of its reactivity towards dienes is computationally costly, implying the need to find the best threshold between accuracy and computational cost. Therefore, the present study was carried out to contribute to solving the problems of large-scale prediction of DA transition state geometries and the multi-scale modeling of C60 fullerene DA reactions. To address the first problem (large-scale prediction of TSs), we have developed a python program named “AMADAR”, which predicts an unlimited number of DA transition states, using only the SMILES strings of the cycloadducts. AMADAR is customizable and allows for the description of intramolecular DA reactions as well as systems resulting in competing paths. In addition, The AMADAR tool contains two separate modules that perform reaction force analyses and atomic decomposition of energy derivatives from the predicted Intrinsic Reaction Coordinates (IRC) paths. The performance of AMADAR was assessed using 2000 DA cycloadducts and showed a success rate of ~ 95%. Most of the errors were due to basis set inconsistencies or convergence issues that we are still working on. Furthermore, a set of 150 IRC paths generated by the AMADAR program were analyzed to get insight into the (a)synchronicity of DA reactions. This investigation confirmed that the reaction force constant 𝜅 (second derivatives of the system energy with respect to the reaction coordinate) was a good indicator of synchronicity in DA reactions. A close inspection of the profile of 𝜅 has enabled us to propose an alternative classification of DA reactions based on their synchronicity degree, in terms of (quasi)-synchronous, moderate asynchronous, asynchronous, and likely two-steps DA reactions. Natural population analyses seemed to indicate that the global maximum of the reaction force constant could be identified with the formation of all the bonds in the reaction site. Finally, the atomic resolution of energy derivatives suggested that the mechanism of the DA reaction involves two inner elementary processes associated with the formation of each C-C bond. A striking mechanistic difference between synchronous and asynchronous DA reactions emerging from this study is that, in asynchronous reactions, the driving and retarding forces are mainly caused by the fast and slow-forming bonds (elementary process) respectively, while in the case of synchronous ones both elementary processes retard and drive the process concomitantly and equivalently. Regarding the DA reaction of C60 fullerene that was considered to illustrate the problem of multiscale modeling, we have constructed 12 ONIOM2 and 10 ONIOM3 models combining five semi-empirical methods (AM1, PM3, PM3MM, PDDG, PM6) and the LDA(SVWN) functional in conjunction with the B3LYP/6-31G(d) level. Then, their accuracy and efficiency were assessed in comparison with the pure B3LYP/6-31G(d) level considering first the DA reaction between C60 and cyclopentadiene whose experimental data are available. Further, different DFT functionals were employed in place of the B3LYP functional to describe the higher-layer of the best ONIOM partition, and the results obtained were compared to experimental data. At this step, the ONIOM2(M06-2X/6-31 G(d): SVWN/STO-3G) model, where the higher layer encompasses the diene and pyracyclene portion of C60, was found to provide the best tradeoff between accuracy and cost, with respect to experimental data. This model showed errors lower than 2.6 and 2.0 kcal/mol for the estimation of the activation and reaction enthalpies respectively. We have also demonstrated, by comparing several ONIOM2(DFT/6-31G(d): SVWN/STO-3G) models, the importance of dispersion corrections in the accurate estimation of reaction and activation energies. Finally, we have considered a set of 21 dienes, including anthracene, 1,3-butadiene, 1,3-cyclopentadiene, furan, thiophene, selenothiophene, pyrrole and their mono-cyano and hydroxyl derivatives to get insight into the DA reaction of C60 using the best ONIOM2(M06-2X/6-31 G(d): SVWN/STO-3G) model. For a given diene and its derivatives, the analysis of frontier molecular orbitals provides a consistent explanation for the substituent effect on the activation barrier. It revealed that electron-donating (withdrawing) groups such as -OH (–CN) cut down on the activation barrier of the reaction by lowering (extending) of the HOMOdiene – LUMOC60 gap and consequently enhancing (weakening) the interaction between the two reactants. Further, the decomposition of the activation energy into the strain and interaction components suggested that, for a given diene, electron-donating groups (here –OH) diminish the height of the activation barrier not only by favoring the attractive interaction between the diene and C60, but also by reducing the strain energy of the system; the opposite effect is observed for electron-withdrawing groups (here –CN). In contrast with some previous findings on typical DA reactions, we could not infer any general rule applicable to the entire dataset for the prediction of activation energies because the latter does not correlate well with either of the TS polarity, electrophilicity of the diene, or the reaction energy. , Thesis (MSc) -- Faculty of Science, Chemistry, 2022
- Full Text:
- Date Issued: 2022-04-06
- Authors: Isamura, Bienfait Kabuyaya
- Date: 2022-04-06
- Subjects: Computational chemistry , Diels-Alder reaction , Python (Computer program language) , Reaction force theory , Fullerenes , Diolefins , AMADAR (Automated workflow for Mechanistic Analysis of Diels-Alder Reactions , ONIOM
- Language: English
- Type: Master's thesis , text
- Identifier: http://hdl.handle.net/10962/232317 , vital:49981
- Description: The [4+2] cycloaddition reaction between conjugated dienes and substituted alkenes is known as the Diels-Alder (DA) reaction, in honor of two German chemists, Otto Diels and Kurt Alder, who first reported this marvelous chemical transformation. The DA reaction is one of the most popular reactions in organic chemistry, allowing for the regio- and stereospecific establishment of six-membered rings with up to four stereogenic centers. This pericyclic reaction has found many applications in areas as diverse as natural products chemistry, polymer chemistry, and agrochemistry. Over the past decades, the mechanism of the Diels-Alder (DA) reaction has been the subject of numerous studies, dealing with questions as diverse as the mechanistic pathway, the synchronicity, the use of catalysts, the effect of solvents and salts, etc. On the other hand, as an example, fullerenes (and particularly [60] fullerene) have been found to act as good dienophiles in DA reactions to the extent that many functionalized fullerenes with interesting applications are still synthesized by reacting C60 with dienes. However, despite the very abundant literature about the mechanism of the DA reaction, some pertinent questions have been still pending, including, without being restricted to, the prediction of transition state (TS) geometries and the modeling of DA reactions involving large systems, such as those of C60 fullerene. It must be emphasized that TSs are not easy to predict and the main reason is that many existing algorithms require that the search is initiated from a good starting point (guess TS), which must be very similar to the actual TS. This problem is even more difficult when many TSs are to be located as may be the case in large-scale studies. Moreover, due to the large size of the C60 molecule, the usage of accurate high-level computational methods in the investigation of its reactivity towards dienes is computationally costly, implying the need to find the best threshold between accuracy and computational cost. Therefore, the present study was carried out to contribute to solving the problems of large-scale prediction of DA transition state geometries and the multi-scale modeling of C60 fullerene DA reactions. To address the first problem (large-scale prediction of TSs), we have developed a python program named “AMADAR”, which predicts an unlimited number of DA transition states, using only the SMILES strings of the cycloadducts. AMADAR is customizable and allows for the description of intramolecular DA reactions as well as systems resulting in competing paths. In addition, The AMADAR tool contains two separate modules that perform reaction force analyses and atomic decomposition of energy derivatives from the predicted Intrinsic Reaction Coordinates (IRC) paths. The performance of AMADAR was assessed using 2000 DA cycloadducts and showed a success rate of ~ 95%. Most of the errors were due to basis set inconsistencies or convergence issues that we are still working on. Furthermore, a set of 150 IRC paths generated by the AMADAR program were analyzed to get insight into the (a)synchronicity of DA reactions. This investigation confirmed that the reaction force constant 𝜅 (second derivatives of the system energy with respect to the reaction coordinate) was a good indicator of synchronicity in DA reactions. A close inspection of the profile of 𝜅 has enabled us to propose an alternative classification of DA reactions based on their synchronicity degree, in terms of (quasi)-synchronous, moderate asynchronous, asynchronous, and likely two-steps DA reactions. Natural population analyses seemed to indicate that the global maximum of the reaction force constant could be identified with the formation of all the bonds in the reaction site. Finally, the atomic resolution of energy derivatives suggested that the mechanism of the DA reaction involves two inner elementary processes associated with the formation of each C-C bond. A striking mechanistic difference between synchronous and asynchronous DA reactions emerging from this study is that, in asynchronous reactions, the driving and retarding forces are mainly caused by the fast and slow-forming bonds (elementary process) respectively, while in the case of synchronous ones both elementary processes retard and drive the process concomitantly and equivalently. Regarding the DA reaction of C60 fullerene that was considered to illustrate the problem of multiscale modeling, we have constructed 12 ONIOM2 and 10 ONIOM3 models combining five semi-empirical methods (AM1, PM3, PM3MM, PDDG, PM6) and the LDA(SVWN) functional in conjunction with the B3LYP/6-31G(d) level. Then, their accuracy and efficiency were assessed in comparison with the pure B3LYP/6-31G(d) level considering first the DA reaction between C60 and cyclopentadiene whose experimental data are available. Further, different DFT functionals were employed in place of the B3LYP functional to describe the higher-layer of the best ONIOM partition, and the results obtained were compared to experimental data. At this step, the ONIOM2(M06-2X/6-31 G(d): SVWN/STO-3G) model, where the higher layer encompasses the diene and pyracyclene portion of C60, was found to provide the best tradeoff between accuracy and cost, with respect to experimental data. This model showed errors lower than 2.6 and 2.0 kcal/mol for the estimation of the activation and reaction enthalpies respectively. We have also demonstrated, by comparing several ONIOM2(DFT/6-31G(d): SVWN/STO-3G) models, the importance of dispersion corrections in the accurate estimation of reaction and activation energies. Finally, we have considered a set of 21 dienes, including anthracene, 1,3-butadiene, 1,3-cyclopentadiene, furan, thiophene, selenothiophene, pyrrole and their mono-cyano and hydroxyl derivatives to get insight into the DA reaction of C60 using the best ONIOM2(M06-2X/6-31 G(d): SVWN/STO-3G) model. For a given diene and its derivatives, the analysis of frontier molecular orbitals provides a consistent explanation for the substituent effect on the activation barrier. It revealed that electron-donating (withdrawing) groups such as -OH (–CN) cut down on the activation barrier of the reaction by lowering (extending) of the HOMOdiene – LUMOC60 gap and consequently enhancing (weakening) the interaction between the two reactants. Further, the decomposition of the activation energy into the strain and interaction components suggested that, for a given diene, electron-donating groups (here –OH) diminish the height of the activation barrier not only by favoring the attractive interaction between the diene and C60, but also by reducing the strain energy of the system; the opposite effect is observed for electron-withdrawing groups (here –CN). In contrast with some previous findings on typical DA reactions, we could not infer any general rule applicable to the entire dataset for the prediction of activation energies because the latter does not correlate well with either of the TS polarity, electrophilicity of the diene, or the reaction energy. , Thesis (MSc) -- Faculty of Science, Chemistry, 2022
- Full Text:
- Date Issued: 2022-04-06
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
- 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
Activity of diverse chalcones against several targets: statistical analysis of a high-throughput virtual screen of a custom chalcone library
- Authors: Sarron, Arthur F D
- Date: 2020
- Subjects: Acetophenone , Benzaldehyde , Ketones , Pyruvate kinase , Drug development , Aromatic compounds , Heat shock proteins
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/116028 , vital:34291
- Description: Chalcone family molecules are well known to have therapeutic proprieties (anti-inflammatory, anti-microbial or anti-cancer, etc). However the mechanism of action in some cases is not well known. A virtual library of this family of compounds was constructed using custom scripts, based on the aldol condensation, and this library was modified further to analogues by expansion of the α,β-unsaturated ketone linker. Acetophenone and benzaldehyde derivatives which are available and purchasable were used as a base to design the chalcone virtual library. 8063 chalcones were constructed and geometrically optimized with Gaussian 09. Their physicochemical characteristics linked to the Lipinski rules were analyzed with Knime and CDK. The entire library was after docked against several targets including HIV-1 integrase, MRSA pyruvate kinase, HSP90, COX-1, COX-2, ALR2, MAOA, MAOB, acetylcholinesterase, butyrylcholinesterase and PLA2. With the exception of MAOA, which does not have a crystal structure ligand, all dockings were validated by redocking the original ligand provided by the literature. These targets are known in the literature to be inhibited by chalcone-derivatives. However, specificity of the particular known chalcone inhibitors to the particular targets is not known. To this end the performance of the generated chalcone library against the list of targets was of interest. The binding energy of ligand-protein complexes was generally good across the library. Statistical analysis including principal component analysis and hierarchical clustering analysis were made in order to investigate for any physical/chemical characteristics which might explain what chalcone features affect the binding energy of the ligand-protein complexes. The spherical polar coordinates defining the orientation of the binding poses were also calculated and used in the statistical analysis. The statistical analysis has allowed us to hypothesize the importance of these radial distances and the polar angles of key atoms in the chalcones in binding to the pyruvate kinase crystal structure. This was validated by the docking of another small library of compound models in which the α,β-unsaturated ketone chain of the chalcone was replaced by incrementally longer conjugated chains. Further studies on the chalcones themselves reveal rotameric systems in both cis and trans-configurations (which may impact binding), and also studied was the effect of Topliss-based modification and its impact of binding to HSP90. Molecular dynamics confirmed good binding of identified chalcone hits.
- Full Text:
- Date Issued: 2020
- Authors: Sarron, Arthur F D
- Date: 2020
- Subjects: Acetophenone , Benzaldehyde , Ketones , Pyruvate kinase , Drug development , Aromatic compounds , Heat shock proteins
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/116028 , vital:34291
- Description: Chalcone family molecules are well known to have therapeutic proprieties (anti-inflammatory, anti-microbial or anti-cancer, etc). However the mechanism of action in some cases is not well known. A virtual library of this family of compounds was constructed using custom scripts, based on the aldol condensation, and this library was modified further to analogues by expansion of the α,β-unsaturated ketone linker. Acetophenone and benzaldehyde derivatives which are available and purchasable were used as a base to design the chalcone virtual library. 8063 chalcones were constructed and geometrically optimized with Gaussian 09. Their physicochemical characteristics linked to the Lipinski rules were analyzed with Knime and CDK. The entire library was after docked against several targets including HIV-1 integrase, MRSA pyruvate kinase, HSP90, COX-1, COX-2, ALR2, MAOA, MAOB, acetylcholinesterase, butyrylcholinesterase and PLA2. With the exception of MAOA, which does not have a crystal structure ligand, all dockings were validated by redocking the original ligand provided by the literature. These targets are known in the literature to be inhibited by chalcone-derivatives. However, specificity of the particular known chalcone inhibitors to the particular targets is not known. To this end the performance of the generated chalcone library against the list of targets was of interest. The binding energy of ligand-protein complexes was generally good across the library. Statistical analysis including principal component analysis and hierarchical clustering analysis were made in order to investigate for any physical/chemical characteristics which might explain what chalcone features affect the binding energy of the ligand-protein complexes. The spherical polar coordinates defining the orientation of the binding poses were also calculated and used in the statistical analysis. The statistical analysis has allowed us to hypothesize the importance of these radial distances and the polar angles of key atoms in the chalcones in binding to the pyruvate kinase crystal structure. This was validated by the docking of another small library of compound models in which the α,β-unsaturated ketone chain of the chalcone was replaced by incrementally longer conjugated chains. Further studies on the chalcones themselves reveal rotameric systems in both cis and trans-configurations (which may impact binding), and also studied was the effect of Topliss-based modification and its impact of binding to HSP90. Molecular dynamics confirmed good binding of identified chalcone hits.
- Full Text:
- Date Issued: 2020
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
Application of computational methods in elucidating the isomerization step in the biosynthesis of coumarins
- Authors: Tshiwawa, Tendamudzimu
- Date: 2019
- Subjects: Coumarins , Isomerization , Biosynthesis , Organic compounds -- Synthesis , Cinnamic acid
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/67646 , vital:29124
- Description: The identity of the enzyme(s) responsible for the biosynthetic transformation of cinnamic acid derivatives to important, naturally occurring coumarins has yet to be established. This study constitutes a high-level theoretical analysis of the possibility that a recently reported molecular mechanism of the synthesis of coumarins from Baylis-Hillman adducts, may provide a viable model for three critical phases in the biosynthetic pathway Particular attention has been given to the first of these phases: i) E→Z isomerisation of the cinnamic acid precursor; ii) Cyclisation (lactonisation) to the hemi-acetal intermediate; and ii) Dehydration to afford the coumarin derivative. In order to accomplish this analysis, an enzyme capable, theoretically, of effecting this E→Z isomerisation required identification, and its potential involvement in the transformation mechanism explored. Combined Molecular Mechanics and high-level Quantum Mechanical/DFT calculations were used to access complementary models of appropriate complexes and relevant processes within the enzyme active sites of a range of eleven Chalcone Isomerase (CHI) enzyme candidates, the structures of which were downloaded from the Protein Data Bank. Detailed B3LYP/6-31+G(d,p) calculations have provided pictures of the relative populations of conformations within the ensemble of conformations available at normal temperatures. Conformations of several protonation states of cinnamic acid derivatives have been studied in this way, and the results obtained showed that coupled protonation and deprotonation of (E)-o-coumaric acid provides a viable approach to achieve the E→Z isomerization. In silico docking of the B3LYP/6-31+G(d,p) optimized (E)-o-coumaric acid derivatives in the active sites of each of the candidate CHI enzymes (CHI) revealed that (E)-o-coumaric acid fits well within the active sites of Medicago Sativa CHI crystallographic structures with 1FM8 showing best potential for not only accommodating (E)-o-coumaric acid , but also providing appropriate protein active site residues to effect the simultaneous protonation and deprotonation of the substrate , two residues being optimally placed to facilitate these critical processes. Further exploration of the chemical properties and qualities of selected CHI enzymes, undertaken using High Throughput Virtual Screening (HTVS), confirmed 1FM8 as a viable choice for further studies of the enzyme-catalysed E→Z isomerization of (E)-o-coumaric acid. A molecular dynamics study, performed to further evaluate the evolution of (E)-o-coumaric acid in the CHI active site over time, showed that the ligand in the 1FM8 active site is not only stable, but also that the desired protein-ligand interactions persist throughout the simulation period to facilitate the E→Z isomerization. An integrated molecular orbital and molecular mechanics (ONIOM) study of the 1FM8-(E)-o-coumaric acid complex, involving the direct protonation and deprotonation of the ligand by protein residues; has provided a plausible mechanism for the E → Z isomerization of (E)-o-coumaric acid within the 1FM8 active site; a transition state complex (with an activation energy of ca. 50 kCal.mol-1) has been located and its connection with both the (E)- and (Z)-o-coumaric acid isomer has been confirmed by Intrinsic Reaction Coordinate (IRC) calculations. More realistic models of the 1FM8-(E)-o-coumaric acid complex, with the inclusion of water solvent molecules, have been obtained at both the QM/MM and adaptive QM/MM levels which simulate the dynamic active site at the QM level. The results indicate that the simultaneous protonation and deprotonation of (E)-o-coumaric acid within the CHI enzyme is a water-mediated process – a conclusion consistent with similar reported processes. Visual inspection of the 1FM8-(Z)-o-coumaric acid complex reveals both the necessary orientation of the phenolic and carboxylic acid moieties of the (Z)-o-coumaric acid and the presence of appropriate, proximal active site residues with the potential to permit catalysis of the subsequent lactonisation and dehydration steps required to generate coumarin.
- Full Text:
- Date Issued: 2019
- Authors: Tshiwawa, Tendamudzimu
- Date: 2019
- Subjects: Coumarins , Isomerization , Biosynthesis , Organic compounds -- Synthesis , Cinnamic acid
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/67646 , vital:29124
- Description: The identity of the enzyme(s) responsible for the biosynthetic transformation of cinnamic acid derivatives to important, naturally occurring coumarins has yet to be established. This study constitutes a high-level theoretical analysis of the possibility that a recently reported molecular mechanism of the synthesis of coumarins from Baylis-Hillman adducts, may provide a viable model for three critical phases in the biosynthetic pathway Particular attention has been given to the first of these phases: i) E→Z isomerisation of the cinnamic acid precursor; ii) Cyclisation (lactonisation) to the hemi-acetal intermediate; and ii) Dehydration to afford the coumarin derivative. In order to accomplish this analysis, an enzyme capable, theoretically, of effecting this E→Z isomerisation required identification, and its potential involvement in the transformation mechanism explored. Combined Molecular Mechanics and high-level Quantum Mechanical/DFT calculations were used to access complementary models of appropriate complexes and relevant processes within the enzyme active sites of a range of eleven Chalcone Isomerase (CHI) enzyme candidates, the structures of which were downloaded from the Protein Data Bank. Detailed B3LYP/6-31+G(d,p) calculations have provided pictures of the relative populations of conformations within the ensemble of conformations available at normal temperatures. Conformations of several protonation states of cinnamic acid derivatives have been studied in this way, and the results obtained showed that coupled protonation and deprotonation of (E)-o-coumaric acid provides a viable approach to achieve the E→Z isomerization. In silico docking of the B3LYP/6-31+G(d,p) optimized (E)-o-coumaric acid derivatives in the active sites of each of the candidate CHI enzymes (CHI) revealed that (E)-o-coumaric acid fits well within the active sites of Medicago Sativa CHI crystallographic structures with 1FM8 showing best potential for not only accommodating (E)-o-coumaric acid , but also providing appropriate protein active site residues to effect the simultaneous protonation and deprotonation of the substrate , two residues being optimally placed to facilitate these critical processes. Further exploration of the chemical properties and qualities of selected CHI enzymes, undertaken using High Throughput Virtual Screening (HTVS), confirmed 1FM8 as a viable choice for further studies of the enzyme-catalysed E→Z isomerization of (E)-o-coumaric acid. A molecular dynamics study, performed to further evaluate the evolution of (E)-o-coumaric acid in the CHI active site over time, showed that the ligand in the 1FM8 active site is not only stable, but also that the desired protein-ligand interactions persist throughout the simulation period to facilitate the E→Z isomerization. An integrated molecular orbital and molecular mechanics (ONIOM) study of the 1FM8-(E)-o-coumaric acid complex, involving the direct protonation and deprotonation of the ligand by protein residues; has provided a plausible mechanism for the E → Z isomerization of (E)-o-coumaric acid within the 1FM8 active site; a transition state complex (with an activation energy of ca. 50 kCal.mol-1) has been located and its connection with both the (E)- and (Z)-o-coumaric acid isomer has been confirmed by Intrinsic Reaction Coordinate (IRC) calculations. More realistic models of the 1FM8-(E)-o-coumaric acid complex, with the inclusion of water solvent molecules, have been obtained at both the QM/MM and adaptive QM/MM levels which simulate the dynamic active site at the QM level. The results indicate that the simultaneous protonation and deprotonation of (E)-o-coumaric acid within the CHI enzyme is a water-mediated process – a conclusion consistent with similar reported processes. Visual inspection of the 1FM8-(Z)-o-coumaric acid complex reveals both the necessary orientation of the phenolic and carboxylic acid moieties of the (Z)-o-coumaric acid and the presence of appropriate, proximal active site residues with the potential to permit catalysis of the subsequent lactonisation and dehydration steps required to generate coumarin.
- Full Text:
- Date Issued: 2019
Enumeration, conformation sampling and population of libraries of peptide macrocycles for the search of chemotherapeutic cardioprotection agents
- Authors: Sigauke, Lester Takunda
- Date: 2019
- Subjects: Peptides -- Synthesis , Macrocyclic compounds , Drug development , Drug discovery , Cardiovascular system -- Diseases -- Prevention , Proteins -- Synthesis
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/116056 , vital:34293
- Description: Peptides are uniquely endowed with features that allow them to perturb previously difficult to drug biomolecular targets. Peptide macrocycles in particular have seen a flurry of recent interest due to their enhanced bioavailability, tunability and specificity. Although these properties make them attractive hit-candidates in early stage drug discovery, knowing which peptides to pursue is non‐trivial due to the magnitude of the peptide sequence space. Computational screening approaches show promise in their ability to address the size of this search space but suffer from their inability to accurately interrogate the conformational landscape of peptide macrocycles. We developed an in‐silico compound enumerator that was tasked with populating a conformationally laden peptide virtual library. This library was then used in the search for cardio‐protective agents (that may be administered, reducing tissue damage during reperfusion after ischemia (heart attacks)). Our enumerator successfully generated a library of 15.2 billion compounds, requiring the use of compression algorithms, conformational sampling protocols and management of aggregated compute resources in the context of a local cluster. In the absence of experimental biophysical data, we performed biased sampling during alchemical molecular dynamics simulations in order to observe cyclophilin‐D perturbation by cyclosporine A and its mitochondrial targeted analogue. Reliable intermediate state averaging through a WHAM analysis of the biased dynamic pulling simulations confirmed that the cardio‐protective activity of cyclosporine A was due to its mitochondrial targeting. Paralleltempered solution molecular dynamics in combination with efficient clustering isolated the essential dynamics of a cyclic peptide scaffold. The rapid enumeration of skeletons from these essential dynamics gave rise to a conformation laden virtual library of all the 15.2 Billion unique cyclic peptides (given the limits on peptide sequence imposed). Analysis of this library showed the exact extent of physicochemical properties covered, relative to the bare scaffold precursor. Molecular docking of a subset of the virtual library against cyclophilin‐D showed significant improvements in affinity to the target (relative to cyclosporine A). The conformation laden virtual library, accessed by our methodology, provided derivatives that were able to make many interactions per peptide with the cyclophilin‐D target. Machine learning methods showed promise in the training of Support Vector Machines for synthetic feasibility prediction for this library. The synergy between enumeration and conformational sampling greatly improves the performance of this library during virtual screening, even when only a subset is used.
- Full Text:
- Date Issued: 2019
- Authors: Sigauke, Lester Takunda
- Date: 2019
- Subjects: Peptides -- Synthesis , Macrocyclic compounds , Drug development , Drug discovery , Cardiovascular system -- Diseases -- Prevention , Proteins -- Synthesis
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/116056 , vital:34293
- Description: Peptides are uniquely endowed with features that allow them to perturb previously difficult to drug biomolecular targets. Peptide macrocycles in particular have seen a flurry of recent interest due to their enhanced bioavailability, tunability and specificity. Although these properties make them attractive hit-candidates in early stage drug discovery, knowing which peptides to pursue is non‐trivial due to the magnitude of the peptide sequence space. Computational screening approaches show promise in their ability to address the size of this search space but suffer from their inability to accurately interrogate the conformational landscape of peptide macrocycles. We developed an in‐silico compound enumerator that was tasked with populating a conformationally laden peptide virtual library. This library was then used in the search for cardio‐protective agents (that may be administered, reducing tissue damage during reperfusion after ischemia (heart attacks)). Our enumerator successfully generated a library of 15.2 billion compounds, requiring the use of compression algorithms, conformational sampling protocols and management of aggregated compute resources in the context of a local cluster. In the absence of experimental biophysical data, we performed biased sampling during alchemical molecular dynamics simulations in order to observe cyclophilin‐D perturbation by cyclosporine A and its mitochondrial targeted analogue. Reliable intermediate state averaging through a WHAM analysis of the biased dynamic pulling simulations confirmed that the cardio‐protective activity of cyclosporine A was due to its mitochondrial targeting. Paralleltempered solution molecular dynamics in combination with efficient clustering isolated the essential dynamics of a cyclic peptide scaffold. The rapid enumeration of skeletons from these essential dynamics gave rise to a conformation laden virtual library of all the 15.2 Billion unique cyclic peptides (given the limits on peptide sequence imposed). Analysis of this library showed the exact extent of physicochemical properties covered, relative to the bare scaffold precursor. Molecular docking of a subset of the virtual library against cyclophilin‐D showed significant improvements in affinity to the target (relative to cyclosporine A). The conformation laden virtual library, accessed by our methodology, provided derivatives that were able to make many interactions per peptide with the cyclophilin‐D target. Machine learning methods showed promise in the training of Support Vector Machines for synthetic feasibility prediction for this library. The synergy between enumeration and conformational sampling greatly improves the performance of this library during virtual screening, even when only a subset is used.
- Full Text:
- Date Issued: 2019
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
The investigation of type-specific features of the copper coordinating AA9 proteins and their effect on the interaction with crystalline cellulose using molecular dynamics studies
- Authors: Moses, Vuyani
- Date: 2018
- Subjects: Copper proteins , Cellulose , Molecular dynamics , Cellulose -- Biodegradation , Bioinformatics
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/58327 , vital:27230
- Description: AA9 proteins are metallo-enzymes which are crucial for the early stages of cellulose degradation. AA9 proteins have been suggested to cleave glycosidic bonds linking cellulose through the use of their Cu2+ coordinating active site. AA9 proteins possess different regioselectivities depending on the resulting cleavage they form and as result, are grouped accordingly. Type 1 AA9 proteins cleave the C1 carbon of cellulose while Type 2 AA9 proteins cleave the C4 carbon and Type 3 AA9 proteins cleave either C1 or C4 carbons. The steric congestion of the AA9 active site has been proposed to be a contributor to the observed regioselectivity. As such, a bioinformatics characterisation of type-specific sequence and structural features was performed. Initially AA9 protein sequences were obtained from the Pfam database and multiple sequence alignment was performed. The sequences were phylogenetically characterised and sequences were grouped into their respective types and sub-groups were identified. A selection analysis was performed on AA9 LPMO types to determine the selective pressure acting on AA9 protein residues. Motif discovery was then performed to identify conserved sequence motifs in AA9 proteins. Once type-specific sequence features were identified structural mapping was performed to assess possible effects on substrate interaction. Physicochemical property analysis was also performed to assess biochemical differences between AA9 LPMO types. Molecular dynamics (MD) simulations were then employed to dynamically assess the consequences of the discovered type-specific features on AA9-cellulose interaction. Due to the absence of AA9 specific force field parameters MD simulations were not readily applicable. As a result, Potential Energy Surface (PES) scans were performed to evaluate the force field parameters for the AA9 active site using the PM6 semi empirical approach and least squares fitting. A Type 1 AA9 active site was constructed from the crystal structure 4B5Q, encompassing only the Cu2+ coordinating residues, the Cu2+ ion and two water residues. Due to the similarity in AA9 active sites, the Type force field parameters were validated on all three AA9 LPMO types. Two MD simulations for each AA9 LPMO types were conducted using two separate Lennard-Jones parameter sets. Once completed, the MD trajectories were analysed for various features including the RMSD, RMSF, radius of gyration, coordination during simulation, hydrogen bonding, secondary structure conservation and overall protein movement. Force field parameters were successfully evaluated and validated for AA9 proteins. MD simulations of AA9 proteins were able to reveal the presence of unique type-specific binding modes of AA9 active sites to cellulose. These binding modes were characterised by the presence of unique type-specific loops which were present in Type 2 and 3 AA9 proteins but not in Type 1 AA9 proteins. The loops were found to result in steric congestion that affects how the Cu2+ ion interacts with cellulose. As a result, Cu2+ binding to cellulose was observed for Type 1 and not Type 2 and 3 AA9 proteins. In this study force field parameters have been evaluated for the Type 1 active site of AA9 proteins and this parameters were evaluated on all three types and binding. Future work will focus on identifying the nature of the reactive oxygen species and performing QM/MM calculations to elucidate the reactive mechanism of all three AA9 LPMO types.
- Full Text:
- Date Issued: 2018
- Authors: Moses, Vuyani
- Date: 2018
- Subjects: Copper proteins , Cellulose , Molecular dynamics , Cellulose -- Biodegradation , Bioinformatics
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/58327 , vital:27230
- Description: AA9 proteins are metallo-enzymes which are crucial for the early stages of cellulose degradation. AA9 proteins have been suggested to cleave glycosidic bonds linking cellulose through the use of their Cu2+ coordinating active site. AA9 proteins possess different regioselectivities depending on the resulting cleavage they form and as result, are grouped accordingly. Type 1 AA9 proteins cleave the C1 carbon of cellulose while Type 2 AA9 proteins cleave the C4 carbon and Type 3 AA9 proteins cleave either C1 or C4 carbons. The steric congestion of the AA9 active site has been proposed to be a contributor to the observed regioselectivity. As such, a bioinformatics characterisation of type-specific sequence and structural features was performed. Initially AA9 protein sequences were obtained from the Pfam database and multiple sequence alignment was performed. The sequences were phylogenetically characterised and sequences were grouped into their respective types and sub-groups were identified. A selection analysis was performed on AA9 LPMO types to determine the selective pressure acting on AA9 protein residues. Motif discovery was then performed to identify conserved sequence motifs in AA9 proteins. Once type-specific sequence features were identified structural mapping was performed to assess possible effects on substrate interaction. Physicochemical property analysis was also performed to assess biochemical differences between AA9 LPMO types. Molecular dynamics (MD) simulations were then employed to dynamically assess the consequences of the discovered type-specific features on AA9-cellulose interaction. Due to the absence of AA9 specific force field parameters MD simulations were not readily applicable. As a result, Potential Energy Surface (PES) scans were performed to evaluate the force field parameters for the AA9 active site using the PM6 semi empirical approach and least squares fitting. A Type 1 AA9 active site was constructed from the crystal structure 4B5Q, encompassing only the Cu2+ coordinating residues, the Cu2+ ion and two water residues. Due to the similarity in AA9 active sites, the Type force field parameters were validated on all three AA9 LPMO types. Two MD simulations for each AA9 LPMO types were conducted using two separate Lennard-Jones parameter sets. Once completed, the MD trajectories were analysed for various features including the RMSD, RMSF, radius of gyration, coordination during simulation, hydrogen bonding, secondary structure conservation and overall protein movement. Force field parameters were successfully evaluated and validated for AA9 proteins. MD simulations of AA9 proteins were able to reveal the presence of unique type-specific binding modes of AA9 active sites to cellulose. These binding modes were characterised by the presence of unique type-specific loops which were present in Type 2 and 3 AA9 proteins but not in Type 1 AA9 proteins. The loops were found to result in steric congestion that affects how the Cu2+ ion interacts with cellulose. As a result, Cu2+ binding to cellulose was observed for Type 1 and not Type 2 and 3 AA9 proteins. In this study force field parameters have been evaluated for the Type 1 active site of AA9 proteins and this parameters were evaluated on all three types and binding. Future work will focus on identifying the nature of the reactive oxygen species and performing QM/MM calculations to elucidate the reactive mechanism of all three AA9 LPMO types.
- Full Text:
- Date Issued: 2018
Synthesis and biolgical screening of potential plasmodium falciparum DXR inhibitors
- Authors: Adeyemi, Christiana Modupe
- Date: 2017-04
- Subjects: Plasmodium falciparum , Enzyme inhibitors , Malaria , Antimalarials , Drug development , Malaria -- Chemotherapy , Isopentenoids -- Synthesis , Fosmidomycin , 1-Deoxy-D-xylulose 5-phosphate
- Language: English
- Type: Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/61790 , vital:28060
- Description: The non-mevalonate isoprenoid pathway, also known as the 1-deoxy-D-xylulose-5- phosphate DXP pathway, is absent in humans, but present in the anopheles mosquito responsible for the transmission of malaria. DXP reductoisomerase - a key enzyme in the DXP pathway in Plasmodium falciparum (PfDXR) has been identified as a target for the design of novel anti-malarial drugs. Fosmidomycin and its acetyl analogue (FR900098) are known to be inhibitors of PfDXR and, in this study, synthetic variations of the fosmidomycin scaffold have led to four series of novel analogues. Particular attention has been centred on the introduction of various substituted benzyl groups in each of these series in order to occupy a recently discovered vacant pocket in the PfDXR active-site and thus enhance ligand-enzyme binding. In the process 160 ligands and precursors have been prepared, no less than 119 of them novel. Fistly, a series of C-benzylated phosphonate esters and phosphonic acids were synthesised, in which the fosmidomycin hydroxamate Mg2+- coordinating moiety was replaced by an amide funtionality and the number of methylene groups in the “hydrophobic patch” between the phosphonate and the hydroxamate moiety was decreased from two to one. Several approaches were explored for this series, the most successful involving reaction of 3- substituted anilines with a-bromo propanoic acid in the presence of the coupling agent 1,1'- carbonyldiimidazole (CDI), followed by Michaelis-Arbuzov phosphonation using triethyl phosphite. Reaction of the resulting chiral phosphonate esters with bromotrimethylsilane gave the corresponding phosphonic acids in good yields. In order to obviate chirality issues, a second series of potential “reverse” fosmidomycin analogues was synthesised by replacing the methylene group adjacent to the the phosphonate moiety with a nitrogen atom. Deprotonation, alkylation and phosphorylation of various amines gave diethyl #-benzylphosphoramidate ester intermediate. Aza-Michael addition of these intermediates, followed by hydrolysis gave the corresponding carboxylic acids which could be reacted with different hydroxylamine hydrochloride derivatives to afford the novel hydroxamic acid derivatives in good yields. Thirdly, a series of a novel #-benzylated phosphoramidate derivatives were prepared as aza- FR900098 analogues. Alkylation of different amines using bromoacetalde-hyde diethylacetal gave a series of N-benzyl-2,2-diethoxyethylamine compounds, which were then elaborated via a futher six steps to afford novel #-benzylated phosphoramidate derivatives. Finally, in order to ensure syn-orientation of the donor atoms in the Mg - coordinating group and, at the same time, introduce conformational constraints in the ligand, the hydrophobic patch and the hydroxamate moiety were replaced by cyclic systems. Several approaches towards the synthesis of such conformationally constrained phosphoramidate analogues from maleic anhydride led to the unexpected isolation of an unprecedented acyclic furfuryl compound, and 1H NMR and DFT-level theoretical studies have been initiated to explore the reaction sequence. A series of #-benzylated phosphoramidate derivatives containing dihydroxy aromatic rings (as the conformationally constrained groups) to replace the hydroxamate moiety, were successfully obtained in six steps from the starting material, 3,4-dihydroxylbenzaldehyde. While in vitro assays have been conducted on all of the synthesised compounds, and some of the ligands show promising anti-malarial inhibitory activity - most especially the conformationally constrained cyclic #-benzylated phosphoramidate series. Interestingly, a number of these compounds has also shown activity against T.brucei - the causative agent of sleeping sickness. In silico docking studies of selected compounds has revealed the capacity of some of the ligands to bind effectively in the PfDXR active-site with the newly introduced benzyl group occupying the adjacent vacant pocket. The physico-chemical properties of these ligands were also explored in order to predict the oral-bioavailability. Most of the ligands obeyed the Lipinski rule of 5, while QSAR methods have been used in an attempt to correlate structural variations and calculated molecular properties with the bioassay data. , Thesis (PhD) -- Faculty of Science, Chemistry, 2017
- Full Text:
- Date Issued: 2017-04
- Authors: Adeyemi, Christiana Modupe
- Date: 2017-04
- Subjects: Plasmodium falciparum , Enzyme inhibitors , Malaria , Antimalarials , Drug development , Malaria -- Chemotherapy , Isopentenoids -- Synthesis , Fosmidomycin , 1-Deoxy-D-xylulose 5-phosphate
- Language: English
- Type: Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/61790 , vital:28060
- Description: The non-mevalonate isoprenoid pathway, also known as the 1-deoxy-D-xylulose-5- phosphate DXP pathway, is absent in humans, but present in the anopheles mosquito responsible for the transmission of malaria. DXP reductoisomerase - a key enzyme in the DXP pathway in Plasmodium falciparum (PfDXR) has been identified as a target for the design of novel anti-malarial drugs. Fosmidomycin and its acetyl analogue (FR900098) are known to be inhibitors of PfDXR and, in this study, synthetic variations of the fosmidomycin scaffold have led to four series of novel analogues. Particular attention has been centred on the introduction of various substituted benzyl groups in each of these series in order to occupy a recently discovered vacant pocket in the PfDXR active-site and thus enhance ligand-enzyme binding. In the process 160 ligands and precursors have been prepared, no less than 119 of them novel. Fistly, a series of C-benzylated phosphonate esters and phosphonic acids were synthesised, in which the fosmidomycin hydroxamate Mg2+- coordinating moiety was replaced by an amide funtionality and the number of methylene groups in the “hydrophobic patch” between the phosphonate and the hydroxamate moiety was decreased from two to one. Several approaches were explored for this series, the most successful involving reaction of 3- substituted anilines with a-bromo propanoic acid in the presence of the coupling agent 1,1'- carbonyldiimidazole (CDI), followed by Michaelis-Arbuzov phosphonation using triethyl phosphite. Reaction of the resulting chiral phosphonate esters with bromotrimethylsilane gave the corresponding phosphonic acids in good yields. In order to obviate chirality issues, a second series of potential “reverse” fosmidomycin analogues was synthesised by replacing the methylene group adjacent to the the phosphonate moiety with a nitrogen atom. Deprotonation, alkylation and phosphorylation of various amines gave diethyl #-benzylphosphoramidate ester intermediate. Aza-Michael addition of these intermediates, followed by hydrolysis gave the corresponding carboxylic acids which could be reacted with different hydroxylamine hydrochloride derivatives to afford the novel hydroxamic acid derivatives in good yields. Thirdly, a series of a novel #-benzylated phosphoramidate derivatives were prepared as aza- FR900098 analogues. Alkylation of different amines using bromoacetalde-hyde diethylacetal gave a series of N-benzyl-2,2-diethoxyethylamine compounds, which were then elaborated via a futher six steps to afford novel #-benzylated phosphoramidate derivatives. Finally, in order to ensure syn-orientation of the donor atoms in the Mg - coordinating group and, at the same time, introduce conformational constraints in the ligand, the hydrophobic patch and the hydroxamate moiety were replaced by cyclic systems. Several approaches towards the synthesis of such conformationally constrained phosphoramidate analogues from maleic anhydride led to the unexpected isolation of an unprecedented acyclic furfuryl compound, and 1H NMR and DFT-level theoretical studies have been initiated to explore the reaction sequence. A series of #-benzylated phosphoramidate derivatives containing dihydroxy aromatic rings (as the conformationally constrained groups) to replace the hydroxamate moiety, were successfully obtained in six steps from the starting material, 3,4-dihydroxylbenzaldehyde. While in vitro assays have been conducted on all of the synthesised compounds, and some of the ligands show promising anti-malarial inhibitory activity - most especially the conformationally constrained cyclic #-benzylated phosphoramidate series. Interestingly, a number of these compounds has also shown activity against T.brucei - the causative agent of sleeping sickness. In silico docking studies of selected compounds has revealed the capacity of some of the ligands to bind effectively in the PfDXR active-site with the newly introduced benzyl group occupying the adjacent vacant pocket. The physico-chemical properties of these ligands were also explored in order to predict the oral-bioavailability. Most of the ligands obeyed the Lipinski rule of 5, while QSAR methods have been used in an attempt to correlate structural variations and calculated molecular properties with the bioassay data. , Thesis (PhD) -- Faculty of Science, Chemistry, 2017
- Full Text:
- Date Issued: 2017-04
Combined in silico approaches towards the identification of novel malarial cysteine protease inhibitors
- Authors: Musyoka, Thommas Mutemi
- Date: 2017
- Language: English
- Type: Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/4488 , vital:20679
- Description: Malaria an infectious disease caused by a group of parasitic organisms of the Plasmodium genus remains a severe public health problem in Africa, South America and parts of Asia. The leading causes for the persistence of malaria are the emergence of drug resistance to common antimalarial drugs, lack of effective vaccines and the inadequate control of mosquito vectors. Worryingly, accumulating evidence shows that the parasite has developed resistant to the current first-line treatment based on artemisinin. Hence, the identification and characterization of novel drug targets and drugs with unique mode of action remains an urgent priority. The successful sequencing and assembly of genomes from several Plasmodium species has opened an opportune window for the identification of new drug targets. Cysteine proteases are one of the major drug targets to be identified so far. The use of cysteine protease inhibitors coupled with gene manipulation studies has defined specific and putative roles of cysteine proteases which include hemoglobin degradation, erythrocyte rupture, immune evasion and erythrocyte invasion, steps which are central for the completion of the Plasmodium parasite life cycle. In an aim to discover potential novel antimalarials, this thesis focussed on falcipains (FPs), a group of four papain-like cysteine proteases from Plasmodium falciparum. Two of these enzymes, FP-2 and FP-3 are the major hemoglobinases and have been validated as drug targets. For the successful elimination of malaria, drugs must be safe and target both human and wild Plasmodium infective forms. Thus, an incipient aim was to identify protein homologs of these two proteases from other Plasmodium species and the host (human). From BLASTP analysis, up to 16 FP-2 and FP-3 homologs were identified (13 plasmodial proteases and 3 human cathepsins). Using in silico characterization approaches, the intra and inter group sequence, structural, phylogenetic and physicochemical differences were determined. To extend previous work (MSc student) involving docking studies on the identified proteins using known FP-2 and FP-3 inhibitors, a South African natural compound and its ZINC analogs, molecular dynamics and binding free energy studies were performed to determine the stabilities and quantification of the strength of interactions between the different protein-ligand complexes. From the results, key structural elements that regulate the binding and selectivity of non-peptidic compounds onto the different proteins were deciphered. Interaction fingerprints and energy decomposition analysis identified key residues and energetic terms that are central for effective ligand binding. This research presents novel insight essential for the structure-based molecular drug design of more potent antimalarial drugs.
- Full Text:
- Date Issued: 2017
- Authors: Musyoka, Thommas Mutemi
- Date: 2017
- Language: English
- Type: Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/4488 , vital:20679
- Description: Malaria an infectious disease caused by a group of parasitic organisms of the Plasmodium genus remains a severe public health problem in Africa, South America and parts of Asia. The leading causes for the persistence of malaria are the emergence of drug resistance to common antimalarial drugs, lack of effective vaccines and the inadequate control of mosquito vectors. Worryingly, accumulating evidence shows that the parasite has developed resistant to the current first-line treatment based on artemisinin. Hence, the identification and characterization of novel drug targets and drugs with unique mode of action remains an urgent priority. The successful sequencing and assembly of genomes from several Plasmodium species has opened an opportune window for the identification of new drug targets. Cysteine proteases are one of the major drug targets to be identified so far. The use of cysteine protease inhibitors coupled with gene manipulation studies has defined specific and putative roles of cysteine proteases which include hemoglobin degradation, erythrocyte rupture, immune evasion and erythrocyte invasion, steps which are central for the completion of the Plasmodium parasite life cycle. In an aim to discover potential novel antimalarials, this thesis focussed on falcipains (FPs), a group of four papain-like cysteine proteases from Plasmodium falciparum. Two of these enzymes, FP-2 and FP-3 are the major hemoglobinases and have been validated as drug targets. For the successful elimination of malaria, drugs must be safe and target both human and wild Plasmodium infective forms. Thus, an incipient aim was to identify protein homologs of these two proteases from other Plasmodium species and the host (human). From BLASTP analysis, up to 16 FP-2 and FP-3 homologs were identified (13 plasmodial proteases and 3 human cathepsins). Using in silico characterization approaches, the intra and inter group sequence, structural, phylogenetic and physicochemical differences were determined. To extend previous work (MSc student) involving docking studies on the identified proteins using known FP-2 and FP-3 inhibitors, a South African natural compound and its ZINC analogs, molecular dynamics and binding free energy studies were performed to determine the stabilities and quantification of the strength of interactions between the different protein-ligand complexes. From the results, key structural elements that regulate the binding and selectivity of non-peptidic compounds onto the different proteins were deciphered. Interaction fingerprints and energy decomposition analysis identified key residues and energetic terms that are central for effective ligand binding. This research presents novel insight essential for the structure-based molecular drug design of more potent antimalarial drugs.
- Full Text:
- Date Issued: 2017
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
Phenomenology: preconceptions and experiences of non-chemists at Rhodes University using milk paint
- Authors: Kelly, Kelvin Leigh
- Date: 2017
- Subjects: Phenomenology , Art and science , Casein , Paint , Chemistry -- Study and teaching , Science -- Study and teaching -- Philosophy
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/37942 , vital:24711
- Description: There exists an ever-increasing crisis in science education where students experience disinterest because of an inability to grasp true understanding of scientific subjects, and therefore there should be a call to increase the research of phenomenology in combination with science education. A rebalance and paradigm shift in the focus of the modes of teaching could result in a great improvement in the learning, comprehension, and intellectual self-confidence of students interested in the sciences. To study this, three research questions were established: How is chemistry perceived by non-chemists; what is the experience of the participants’ during the chemistry practical in a laboratory and; do the participants’ perspectives about chemistry change during the experience. The performed study consisted of a chemistry practical, two art works and, in some cases, an interview. Nine participants were asked to create the art under specific instructions of points of focus, namely their preconceptions prior to the practical (Artwork 1) and their lived experience during the practical (Artwork 2). Participants’ artworks were examined using methods of visual semiotics and classical art analysis techniques, looking at line, shape, and colour choice. The iterative analysis of the interviews from participants 1, 2, 7, and 9 coded with ATLAS.ti 7 software, led to the emergence of themes that constitute the core of the participants’ experience. This phenomenological study presents a path to engage the non-chemist with processes taking place in the laboratory by using ‘Kitchen Chemistry’ and illustrates how a phenomenological engagement with chemistry can make the subject more applicable to the general population of non-chemists.
- Full Text:
- Date Issued: 2017
- Authors: Kelly, Kelvin Leigh
- Date: 2017
- Subjects: Phenomenology , Art and science , Casein , Paint , Chemistry -- Study and teaching , Science -- Study and teaching -- Philosophy
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/37942 , vital:24711
- Description: There exists an ever-increasing crisis in science education where students experience disinterest because of an inability to grasp true understanding of scientific subjects, and therefore there should be a call to increase the research of phenomenology in combination with science education. A rebalance and paradigm shift in the focus of the modes of teaching could result in a great improvement in the learning, comprehension, and intellectual self-confidence of students interested in the sciences. To study this, three research questions were established: How is chemistry perceived by non-chemists; what is the experience of the participants’ during the chemistry practical in a laboratory and; do the participants’ perspectives about chemistry change during the experience. The performed study consisted of a chemistry practical, two art works and, in some cases, an interview. Nine participants were asked to create the art under specific instructions of points of focus, namely their preconceptions prior to the practical (Artwork 1) and their lived experience during the practical (Artwork 2). Participants’ artworks were examined using methods of visual semiotics and classical art analysis techniques, looking at line, shape, and colour choice. The iterative analysis of the interviews from participants 1, 2, 7, and 9 coded with ATLAS.ti 7 software, led to the emergence of themes that constitute the core of the participants’ experience. This phenomenological study presents a path to engage the non-chemist with processes taking place in the laboratory by using ‘Kitchen Chemistry’ and illustrates how a phenomenological engagement with chemistry can make the subject more applicable to the general population of non-chemists.
- Full Text:
- Date Issued: 2017
Physical organic studies of substituted norbornyl systems: aspects of mechanisms and chirality
- Authors: Singh, Alicia
- Date: 2017
- Language: English
- Type: Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/50558 , vital:25999
- Description: Fenchone and camphor are essential natural products that are available optically pure and contribute to the chiral pool in asymmetric synthesis. Further, they are both derivatives of norbornane, a structure that undergoes a remarkable diversity of rearrangements in acidic conditions. This work explores two aspects of the camphor/fenchone derived systems. Firstly an attempt to clarify rearrangement mechanisms on a camphor system successfully via deuterium labelling and unsuccessfully via derivatization of fenchone (with rearrangement) to produce other 13C-labelled camphor substitutions, has resulted in confirmation of a theoretically proposed, highly concerted Wagner-Meerwein, 6,2 - hydride shift, Wagner-Meerwein rearrangement in competition with an associated 2,3-methide shift. Kinetics and activation parameters for many steps have been resolved in this rearrangement of the deuterium labelled camphor-derived tosylate system to two pairs of isotopomers. Further, the kinetics and formation of an unexpected pair of dimers encountered during the scheme for 13C labelling are investigated in detail. These dimers (forming during the initial stages of the synthetic scheme) are unusual in that they are not expected rotamers of each other, but diastereomers resulting from chirality of a sulfur atom in a sulfite moiety. A feasible mechanism of formation that matches the kinetics has been proposed in this unexpectedly complex system, and thermodynamic parameters have been determined. The second aspect of substituted norbornyl systems pertains to their chirality, and the influence of this chirality on reaction mixtures, with an aim to identify novel chiral micellar catalysts for use in heterogeneous reaction mixtures. Headway has been made towards the synthesis of the appropriate surfactants to be used in the construction of these micelles, but extensive molecular dynamics simulations have illustrated the feasibility of forming the stable chiral micelles in a dual-solvent system, and detail precisely the influence of chirality on surrounding media. These studies add important physical organic data as well as show the immense possibilities pertaining to substituted norbornane systems.
- Full Text:
- Date Issued: 2017
- Authors: Singh, Alicia
- Date: 2017
- Language: English
- Type: Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/50558 , vital:25999
- Description: Fenchone and camphor are essential natural products that are available optically pure and contribute to the chiral pool in asymmetric synthesis. Further, they are both derivatives of norbornane, a structure that undergoes a remarkable diversity of rearrangements in acidic conditions. This work explores two aspects of the camphor/fenchone derived systems. Firstly an attempt to clarify rearrangement mechanisms on a camphor system successfully via deuterium labelling and unsuccessfully via derivatization of fenchone (with rearrangement) to produce other 13C-labelled camphor substitutions, has resulted in confirmation of a theoretically proposed, highly concerted Wagner-Meerwein, 6,2 - hydride shift, Wagner-Meerwein rearrangement in competition with an associated 2,3-methide shift. Kinetics and activation parameters for many steps have been resolved in this rearrangement of the deuterium labelled camphor-derived tosylate system to two pairs of isotopomers. Further, the kinetics and formation of an unexpected pair of dimers encountered during the scheme for 13C labelling are investigated in detail. These dimers (forming during the initial stages of the synthetic scheme) are unusual in that they are not expected rotamers of each other, but diastereomers resulting from chirality of a sulfur atom in a sulfite moiety. A feasible mechanism of formation that matches the kinetics has been proposed in this unexpectedly complex system, and thermodynamic parameters have been determined. The second aspect of substituted norbornyl systems pertains to their chirality, and the influence of this chirality on reaction mixtures, with an aim to identify novel chiral micellar catalysts for use in heterogeneous reaction mixtures. Headway has been made towards the synthesis of the appropriate surfactants to be used in the construction of these micelles, but extensive molecular dynamics simulations have illustrated the feasibility of forming the stable chiral micelles in a dual-solvent system, and detail precisely the influence of chirality on surrounding media. These studies add important physical organic data as well as show the immense possibilities pertaining to substituted norbornane systems.
- Full Text:
- Date Issued: 2017
Structural analysis of proteases from South African HIV-1 (subtype C) patients undergoing Lopinavir treatment, using comparative modeling, ligand-docking and molecular dynamics
- Authors: Sheik-Amamuddy, Olivier
- Date: 2017
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/4931 , vital:20744
- Description: HIV is regarded as one of the most devastating infectious diseases of the last few decades, and has a high prevalence in South Africa, subtype C being the most common. Palliative measures used to fight HIV involve the use various types of inhibitors, including the use of HIV protease inhibitors. Representatives from this class of inhibitors are gradually losing their efficacy due to development of resistance mutations from HIV-1. In this study, compounds from the South African Natural Compound Database (SANCDB) were screened against HIV-1 protease models generated from protease protein sequences belonging to 11 South African HIV patients before and after treatment with Lopinavir. The effect of Lopinavir on the alteration of drug-binding affinity before and after treatment is investigated by molecular docking of the protease against other FDA-approved drugs and detection of mutation types using the HIVdb tool. A network representation of hydrogen bonding between docked ligands and their receptor proteases has been developed and a profiling method of visualizing receptor-ligand docking energies at the local level is presented. Four potential HIV-1 protease inhibitors were identified from the list of 599 natural compounds on the basis of receptor conformation and binding free energy. Ligand stabilities were monitored by 20ns molecular dynamics runs using the GROMACS software.
- Full Text:
- Date Issued: 2017
- Authors: Sheik-Amamuddy, Olivier
- Date: 2017
- Language: English
- Type: Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/4931 , vital:20744
- Description: HIV is regarded as one of the most devastating infectious diseases of the last few decades, and has a high prevalence in South Africa, subtype C being the most common. Palliative measures used to fight HIV involve the use various types of inhibitors, including the use of HIV protease inhibitors. Representatives from this class of inhibitors are gradually losing their efficacy due to development of resistance mutations from HIV-1. In this study, compounds from the South African Natural Compound Database (SANCDB) were screened against HIV-1 protease models generated from protease protein sequences belonging to 11 South African HIV patients before and after treatment with Lopinavir. The effect of Lopinavir on the alteration of drug-binding affinity before and after treatment is investigated by molecular docking of the protease against other FDA-approved drugs and detection of mutation types using the HIVdb tool. A network representation of hydrogen bonding between docked ligands and their receptor proteases has been developed and a profiling method of visualizing receptor-ligand docking energies at the local level is presented. Four potential HIV-1 protease inhibitors were identified from the list of 599 natural compounds on the basis of receptor conformation and binding free energy. Ligand stabilities were monitored by 20ns molecular dynamics runs using the GROMACS software.
- Full Text:
- Date Issued: 2017