1
|
Soni M, Srivastava G, Ramalingam K, Shakya AK, Siddiqi MI, Pratap JV. Identification of potent inhibitors for Leishmania donovani homoserine kinase: an integrated in silico and kinetic study. J Biomol Struct Dyn 2023:1-16. [PMID: 37962849 DOI: 10.1080/07391102.2023.2279279] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Accepted: 10/30/2023] [Indexed: 11/15/2023]
Abstract
Leishmaniasis is caused by ∼20 species of Leishmania that affects millions in endemic areas. Available therapies are not sufficient to effectively control the disease, cause severe side effects and eventually lead to drug resistance, making the discovery of novel therapeutic molecules an immediate need. Molecular target-based drug discovery, where the target is a defined molecular gene, protein or a mechanism, is a rationale driven approach for novel therapeutics. Humans obtain the essential amino acid such as threonine from dietary sources, while Leishmania synthesize it de-novo. Enzymes of the threonine biosynthesis pathway, including the rate limiting Homoserine kinase (HSK) which converts L-homoserine into ortho-phospho homoserine are thus attractive targets for rationale driven therapy. The absence of HSK in humans and its presence in Leishmania donovani enhances the opportunity to exploit HSK as a molecular target for anti-leishmanials therapeutic development. In this study, we utilize structure-based high throughput drug discovery (SBDD), followed by biochemical validation and identified two potential inhibitors (RH00038 and S02587) from Maybridge chemical library that targets L. donovani HSK. These two inhibitors effectively induced the mortality of Leishmania donovani in both amastigote and promastigote stages, with one of them being specific to parasite and twice as effective as the standard therapeutic molecule.
Collapse
Affiliation(s)
- Mohini Soni
- Biochemistry and Structural Biology Division, CSIR-Central Drug Research Institute, Lucknow, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, UP, India
| | - Gaurava Srivastava
- Biochemistry and Structural Biology Division, CSIR-Central Drug Research Institute, Lucknow, India
| | - Karthik Ramalingam
- Biochemistry and Structural Biology Division, CSIR-Central Drug Research Institute, Lucknow, India
| | - Anil Kumar Shakya
- Biochemistry and Structural Biology Division, CSIR-Central Drug Research Institute, Lucknow, India
| | - Mohammad Imran Siddiqi
- Biochemistry and Structural Biology Division, CSIR-Central Drug Research Institute, Lucknow, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, UP, India
| | - J Venkatesh Pratap
- Biochemistry and Structural Biology Division, CSIR-Central Drug Research Institute, Lucknow, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, UP, India
| |
Collapse
|
2
|
Meshram RJ, Bagul KT, Aouti SU, Shirsath AM, Duggal H, Gacche RN. Modeling and simulation study to identify threonine synthase as possible drug target in Leishmania major. Mol Divers 2020; 25:1679-1700. [PMID: 32737682 DOI: 10.1007/s11030-020-10129-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Accepted: 07/22/2020] [Indexed: 01/16/2023]
Abstract
Leishmaniasis is one of the most neglected tropical diseases that demand immediate attention to the identification of new drug targets and effective drug candidates. The present study demonstrates the possibility of using threonine synthase (TS) as a putative drug target in leishmaniasis disease management. We report the construction of an effective homology model of the enzyme that appears to be structurally as well as functionally well conserved. The 200 nanosecond molecular dynamics data on TS with and without pyridoxal phosphate (PLP) shed light on mechanistic details of PLP-induced conformational changes. Moreover, we address some important structural and dynamic interactions in the PLP binding region of TS that are in good agreement with previously speculated crystallographic estimations. Additionally, after screening more than 44,000 compounds, we propose 10 putative inhibitor candidates for TS based on virtual screening data and refined Molecular Mechanics Generalized Born Surface Area calculations. We expect that structural and functional dynamics data disclosed in this study will help initiate experimental endeavors toward establishing TS as an effective antileishmanial drug target.
Collapse
Affiliation(s)
- Rohan J Meshram
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India.
| | - Kamini T Bagul
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| | - Snehal U Aouti
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| | - Akshay M Shirsath
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| | - Harleen Duggal
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| | - Rajesh N Gacche
- Department of Biotechnology, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| |
Collapse
|
3
|
Meshram RJ, Shirsath A, Aouti S, Bagul K, Gacche RN. Molecular modeling and simulation study of homoserine kinase as an effective leishmanial drug target. J Mol Model 2020; 26:218. [PMID: 32720228 DOI: 10.1007/s00894-020-04473-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2019] [Accepted: 07/14/2020] [Indexed: 11/25/2022]
Abstract
Leishmaniasis is a tropical neglected disease that imposes major health concerns in many endemic countries worldwide and requires urgent attention to the identification of new drug targets as well as drug candidates. In the current study, we propose homoserine kinase (HSK) inhibition as a strategy to induce pathogen mortality via generating threonine deficiency. We introduce a homology-based molecular model of leishmanial HSK that appears to possess all conserved structural as well as functional features in the GHMP kinase family. Furthermore, 200 ns molecular dynamics data of the enzyme in open and closed state attempts to provide the mechanistic details involved in the substrate as well as phosphate binding to this enzyme. We discuss the structural and functional significance of movements involved in various loops (motif 1, 2, 3) and lips (upper and lower) in the transition of leishmanial HSK from closed to open state. Virtual screening data of more than 40,000 compounds from the present investigation tries to identify a few potential HSK inhibitors that possess important features to act as efficient HSK inhibitors. These compounds can be considered an effective starting point for the identification of novel drug-like scaffolds. We hope the structural wealth that is offered in this report will be utilized in designing competent experimental and therapeutic interventions for leishmaniasis management. Graphical abstract.
Collapse
Affiliation(s)
- Rohan J Meshram
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India.
| | - Akshay Shirsath
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| | - Snehal Aouti
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| | - Kamini Bagul
- Bioinformatics Centre, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| | - Rajesh N Gacche
- Department of Biotechnology, Savitribai Phule Pune University, Pune, Maharashtra, 411007, India
| |
Collapse
|
4
|
Chen CL, Mermoud JC, Paul LN, Steussy CN, Stauffacher CV. Mevalonate 5-diphosphate mediates ATP binding to the mevalonate diphosphate decarboxylase from the bacterial pathogen Enterococcus faecalis. J Biol Chem 2017; 292:21340-21351. [PMID: 29025876 PMCID: PMC5766736 DOI: 10.1074/jbc.m117.802223] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2017] [Revised: 09/23/2017] [Indexed: 01/07/2023] Open
Abstract
The mevalonate pathway produces isopentenyl diphosphate (IPP), a building block for polyisoprenoid synthesis, and is a crucial pathway for growth of the human bacterial pathogen Enterococcus faecalis The final enzyme in this pathway, mevalonate diphosphate decarboxylase (MDD), acts on mevalonate diphosphate (MVAPP) to produce IPP while consuming ATP. This essential enzyme has been suggested as a therapeutic target for the treatment of drug-resistant bacterial infections. Here, we report functional and structural studies on the mevalonate diphosphate decarboxylase from E. faecalis (MDDEF). The MDDEF crystal structure in complex with ATP (MDDEF-ATP) revealed that the phosphate-binding loop (amino acids 97-105) is not involved in ATP binding and that the phosphate tail of ATP in this structure is in an outward-facing position pointing away from the active site. This suggested that binding of MDDEF to MVAPP is necessary to guide ATP into a catalytically favorable position. Enzymology experiments show that the MDDEF performs a sequential ordered bi-substrate reaction with MVAPP as the first substrate, consistent with the isothermal titration calorimetry (ITC) experiments. On the basis of ITC results, we propose that this initial prerequisite binding of MVAPP enhances ATP binding. In summary, our findings reveal a substrate-induced substrate-binding event that occurs during the MDDEF-catalyzed reaction. The disengagement of the phosphate-binding loop concomitant with the alternative ATP-binding configuration may provide the structural basis for antimicrobial design against these pathogenic enterococci.
Collapse
Affiliation(s)
| | | | - Lake N. Paul
- the Biophysical Analysis Laboratory, Bindley Bioscience Center, Purdue University, West Lafayette, Indiana 47906
| | | | - Cynthia V. Stauffacher
- From the Department of Biological Sciences and ,Purdue University Center for Cancer Research, Purdue University, West Lafayette, Indiana 47907 and , To whom correspondence should be addressed:
Dept. of Biological Sciences, Purdue University, Hockmeyer Hall, Rm. 327, 240 South Martin Jischke Dr., West Lafayette, IN 47907. Tel.:
765-494-4937; Fax:
765-496-1189; E-mail:
| |
Collapse
|
5
|
Improved production of L-threonine in Escherichia coli by use of a DNA scaffold system. Appl Environ Microbiol 2012; 79:774-82. [PMID: 23160128 DOI: 10.1128/aem.02578-12] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Despite numerous approaches for the development of l-threonine-producing strains, strain development is still hampered by the intrinsic inefficiency of metabolic reactions caused by simple diffusion and random collisions of enzymes and metabolites. A scaffold system, which can promote the proximity of metabolic enzymes and increase the local concentration of intermediates, was reported to be one of the most promising solutions. Here, we report an improvement in l-threonine production in Escherichia coli using a DNA scaffold system, in which a zinc finger protein serves as an adapter for the site-specific binding of each enzyme involved in l-threonine production to a precisely ordered location on a DNA double helix to increase the proximity of enzymes and the local concentration of metabolites to maximize production. The optimized DNA scaffold system for l-threonine production significantly increased the efficiency of the threonine biosynthetic pathway in E. coli, substantially reducing the production time for l-threonine (by over 50%). In addition, this DNA scaffold system enhanced the growth rate of the host strain by reducing the intracellular concentration of toxic intermediates, such as homoserine. Our DNA scaffold system can be used as a platform technology for the construction and optimization of artificial metabolic pathways as well as for the production of many useful biomaterials.
Collapse
|
6
|
Biochemical characterization of pantoate kinase, a novel enzyme necessary for coenzyme A biosynthesis in the Archaea. J Bacteriol 2012; 194:5434-43. [PMID: 22865846 DOI: 10.1128/jb.06624-11] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Although bacteria and eukaryotes share a pathway for coenzyme A (CoA) biosynthesis, we previously clarified that most archaea utilize a distinct pathway for the conversion of pantoate to 4'-phosphopantothenate. Whereas bacteria/eukaryotes use pantothenate synthetase and pantothenate kinase (PanK), the hyperthermophilic archaeon Thermococcus kodakarensis utilizes two novel enzymes: pantoate kinase (PoK) and phosphopantothenate synthetase (PPS). Here, we report a detailed biochemical examination of PoK from T. kodakarensis. Kinetic analyses revealed that the PoK reaction displayed Michaelis-Menten kinetics toward ATP, whereas substrate inhibition was observed with pantoate. PoK activity was not affected by the addition of CoA/acetyl-CoA. Interestingly, PoK displayed broad nucleotide specificity and utilized ATP, GTP, UTP, and CTP with comparable k(cat)/K(m) values. Sequence alignment of 27 PoK homologs revealed seven conserved residues with reactive side chains, and variant proteins were constructed for each residue. Activity was not detected when mutations were introduced to Ser104, Glu134, and Asp143, suggesting that these residues play vital roles in PoK catalysis. Kinetic analysis of the other variant proteins, with mutations S28A, H131A, R155A, and T186A, indicated that all four residues are involved in pantoate recognition and that Arg155 and Thr186 play important roles in PoK catalysis. Gel filtration analyses of the variant proteins indicated that Thr186 is also involved in dimer assembly. A sequence comparison between PoK and other members of the GHMP kinase family suggests that Ser104 and Glu134 are involved in binding with phosphate and Mg(2+), respectively, while Asp143 is the base responsible for proton abstraction from the pantoate hydroxy group.
Collapse
|
7
|
Morneau DJK, Abouassaf E, Skanes JE, Aitken SM. Development of a continuous assay and steady-state characterization of Escherichia coli threonine synthase. Anal Biochem 2012; 423:78-85. [PMID: 22289691 DOI: 10.1016/j.ab.2012.01.004] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2011] [Revised: 01/02/2012] [Accepted: 01/09/2012] [Indexed: 10/14/2022]
Abstract
Threonine synthase (TS) catalyzes the hydrolysis of O-phospho-L-homoserine (OPHS) to produce L-threonine (L-Thr) and inorganic phosphate. Here, we report a simplified purification protocol for the OPHS substrate and a continuous, coupled-coupled, spectrophotometric TS assay. The sequential actions of threonine deaminase (TD) and hydroxyisocaproate dehydrogenase (HO-HxoDH) convert the L-Thr product of TS to α-ketobutyrate (α-KB) and then to 2-hydroxybutyrate, respectively, and are monitored as the decrease in absorbance at 340 nm resulting from the concomitant oxidation of β-nicotinamide adenine dinucleotide (NADH) to NAD(+) by HO-HxoDH. The effect of pH on the activities of Escherichia coli TD and Lactobacillus delbrueckii HO-HxoDH was determined to establish this continuous assay as suitable for steady-state characterization and to facilitate the optimization of coupling enzyme concentrations under different assay conditions to enable studies of TS across phyla. To validate this assay, TS from E. coli was characterized. The kinetic parameters (k(cat)=4s(-1) and K(m)=0.34 mM) and the pH optimum of 8.7, determined using the continuous assay, are consistent with values reported for this enzyme based on the discontinuous malachite green assay. The k(cat)/K(m)(OPHS) versus pH profile of E. coli TS is bell-shaped, and the apparent pK(a) values for the acidic and basic limbs are 7.1 and 10.4, respectively.
Collapse
|
8
|
Najdi TS, Yang CR, Shapiro BE, Hatfield GW, Mjolsness ED. APPLICATION OF A GENERALIZED MWC MODEL FOR THE MATHEMATICAL SIMULATION OF METABOLIC PATHWAYS REGULATED BY ALLOSTERIC ENZYMES. J Bioinform Comput Biol 2011; 4:335-55. [PMID: 16819787 DOI: 10.1142/s0219720006001862] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2005] [Accepted: 01/31/2006] [Indexed: 11/18/2022]
Abstract
In our effort to elucidate the systems biology of the model organism, Escherichia coli, we have developed a mathematical model that simulates the allosteric regulation for threonine biosynthesis pathway starting from aspartate. To achieve this goal, we used kMech, a Cellerator language extension that describes enzyme mechanisms for the mathematical modeling of metabolic pathways. These mechanisms are converted by Cellerator into ordinary differential equations (ODEs) solvable by Mathematica™. In this paper, we describe a more flexible model in Cellerator, which generalizes the Monod, Wyman, Changeux (MWC) model for enzyme allosteric regulation to allow for multiple substrate, activator and inhibitor binding sites. Furthermore, we have developed a model that describes the behavior of the bifunctional allosteric enzyme aspartate kinase I-homoserine dehydrogenase I (AKI-HDHI). This model predicts the partition of enzyme activities in the steady state which paves the way for a more generalized prediction of the behavior of bifunctional enzymes.
Collapse
Affiliation(s)
- Tarek S Najdi
- Institute for Genomics and Bioinformatics, Department of Microbiology and Molecular Genetics, College of Medicine, University of California, Irvine, Irvine, California 92697, USA.
| | | | | | | | | |
Collapse
|
9
|
De Pascale G, Griffiths EJ, Shakya T, Nazi I, Wright GD. Identification and Characterization of New Inhibitors of Fungal Homoserine Kinase. Chembiochem 2011; 12:1179-82. [DOI: 10.1002/cbic.201100121] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2011] [Indexed: 11/09/2022]
|
10
|
The role of the active site residues in human galactokinase: implications for the mechanisms of GHMP kinases. Bioorg Chem 2011; 39:120-6. [PMID: 21474160 DOI: 10.1016/j.bioorg.2011.03.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2011] [Revised: 03/01/2011] [Accepted: 03/02/2011] [Indexed: 11/20/2022]
Abstract
Galactokinase catalyses the phosphorylation of galactose at the expense of ATP. Like other members of the GHMP family of kinases it is postulated to function through an active site base mechanism in which Asp-186 abstracts a proton from galactose. This asparate residue was altered to alanine and to asparagine by site-directed mutagenesis of the corresponding gene. This resulted in variant enzyme with no detectable galactokinase activity. Alteration of Arg-37, which lies adjacent to Asp-186 and is postulated to assist the catalytic base, to lysine resulted in an active enzyme. However, alteration of this residue to glutamate abolished activity. All the variant enzymes, except the arginine to lysine substitution, were structurally unstable (as judged by native gel electrophoresis in the presence of urea) compared to the wild type. This suggests that the lack of activity results from this structural instability, in addition to any direct effects on the catalytic mechanism. Computational estimations of the pK(a) values of the arginine and aspartate residues, suggest that Arg-37 remains protonated throughout the catalytic cycle whereas Asp-186 has an abnormally high pK(a) value (7.18). Quantum mechanics/molecular mechanics (QM/MM) calculations suggest that Asp-186 moves closer to the galactose molecule during catalysis. The experimental and theoretical studies presented here argue for a mechanism in which the C(1)-OH bond in the sugar is weakened by the presence of Asp-186 thus facilitating nucleophilic attack by the oxygen atom on the γ-phosphorus of ATP.
Collapse
|
11
|
Stockbridge RB, Wolfenden R. The intrinsic reactivity of ATP and the catalytic proficiencies of kinases acting on glucose, N-acetylgalactosamine, and homoserine: a thermodynamic analysis. J Biol Chem 2009; 284:22747-57. [PMID: 19531469 DOI: 10.1074/jbc.m109.017806] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
To evaluate the rate enhancements produced by representative kinases and their thermodynamic basis, rate constants were determined as a function of changing temperature for 1) the spontaneous methanolysis of ATP and 2) reactions catalyzed by kinases to which different mechanisms of action have been ascribed. For each of these enzymes, the minor effects of changing viscosity indicate that k(cat)/K(m) is governed by the central chemical events in the enzyme-substrate complex rather than by enzyme-substrate encounter. Individual Arrhenius plots, obtained at intervals between pH 4.8 and 11.0, yielded Delta H(#) and T Delta S(#) for the nonenzymatic methanolysis of ATP(2-), ATP(3-), and ATP(4-) in the absence of Mg(2+). The addition of Mg(2+) led to partly compensating changes in Delta H(#) and T Delta S(#), accelerating the nonenzymatic methanolysis of ATP 11-fold at pH 7 and 25 degrees C. The rate enhancements produced by yeast hexokinase, homoserine kinase, and N-acetylgalactosamine kinase (obtained by comparison of their k(cat)/K(m) values in the presence of saturating phosphoryl acceptor with the second order rate constant for methanolysis of MgATP) ranged between 10(12)- and 10(14)-fold. Their nominal affinities for the altered substrates in the transition state were 2.1 x 10(-16) m for N-acetylgalactosamine kinase, 7.4 x 10(-17) m for homoserine kinase, and 6.4 x 10(-18) m for hexokinase. Compared with nonenzymatic phosphoryl transfer, all three kinases were found to produce major reductions in the entropy of activation, in accord with the likelihood that substrate juxtaposition and desolvation play prominent roles in their catalytic action.
Collapse
Affiliation(s)
- Randy B Stockbridge
- Department of Biochemistry and Biophysics, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | | |
Collapse
|
12
|
Fan C, Fromm HJ, Bobik TA. Kinetic and functional analysis of L-threonine kinase, the PduX enzyme of Salmonella enterica. J Biol Chem 2009; 284:20240-8. [PMID: 19509296 DOI: 10.1074/jbc.m109.027425] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The PduX enzyme of Salmonella enterica is an l-threonine kinase used for the de novo synthesis of coenzyme B(12) and the assimilation of cobyric acid. PduX with an N-terminal histidine tag (His(8)-PduX) was produced in Escherichiacoli and purified. The recombinant enzyme was soluble and active. Kinetic analysis indicated a steady-state Ordered Bi Bi complex mechanism in which ATP is the first substrate to bind. Based on a multiple sequence alignment of PduX homologues and other GHMP (galactokinase, homoserine kinase, mevalonate kinase, and phosphomevalonate kinase) family members, 14 PduX variants having changes at 10 conserved serine/threonine and aspartate/glutamate sites were constructed by site-directed mutagenesis. Each variant was produced in E. coli and purified. Comparison of the circular dichroism spectra and kinetic properties of the PduX variants with those of the wild-type enzyme indicated that Glu-24 and Asp-135 are needed for proper folding, Ser-99 and Glu-132 are used for ATP binding, and Ser-253 and Ser-255 are critical to l-threonine binding whereas Ser-100 is essential to catalysis, but its precise role is uncertain. The studies reported here are the first to investigate the kinetic and catalytic mechanisms of l-threonine kinase from any organism.
Collapse
Affiliation(s)
- Chenguang Fan
- Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, Iowa 50011, USA
| | | | | |
Collapse
|
13
|
Najdi TS, Yang CR, Shapiro BE, Hatfield GW, Mjolsness ED. Application of a generalized MWC model for the mathematical simulation of metabolic pathways regulated by allosteric enzymes. PROCEEDINGS. IEEE COMPUTATIONAL SYSTEMS BIOINFORMATICS CONFERENCE 2005:279-88. [PMID: 16447985 DOI: 10.1109/csb.2005.15] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
In our effort to elucidate the systems biology of the model organism, Escherichia coli, we have developed a mathematical model that simulates the allosteric regulation for threonine biosynthesis pathway starting from aspartate. To achieve this goal, we used kMech, a Cellerator language extension that describes enzyme mechanisms for the mathematical modeling of metabolic pathways. These mechanisms are converted by Cellerator into ordinary differential equations (ODEs) solvable by Mathematica. In this paper, we describe a more flexible model in Cellerator, which generalizes the Monod, Wyman, Changeux (MWC) model for enzyme allosteric regulation to allow for multiple substrate, activator and inhibitor binding sites. Furthermore, we have developed a model that describes the behavior of the bifunctional allosteric enzyme aspartate Kinase I-Homoserine Dehydrogenase I (AKI-HDHI). This model predicts the partition of enzyme activities in the steady state which paves a way for a more generalized prediction of the behavior of bifunctional enzymes.
Collapse
Affiliation(s)
- Tarek S Najdi
- Institute for Genomics and Bioinformatics, University of California at Irvine, Irvine, CA 92697, USA
| | | | | | | | | |
Collapse
|
14
|
Chassagnole C, Raïs B, Quentin E, Fell DA, Mazat JP. An integrated study of threonine-pathway enzyme kinetics in Escherichia coli. Biochem J 2001; 356:415-23. [PMID: 11368768 PMCID: PMC1221852 DOI: 10.1042/0264-6021:3560415] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
We have determined the kinetic parameters of the individual steps of the threonine pathway from aspartate in Escherichia coli under a single set of experimental conditions chosen to be physiologically relevant. Our aim was to summarize the kinetic behaviour of each enzyme in a single tractable equation that takes into account the effect of the products as competitive inhibitors of the substrates in the forward reaction and also, when appropriate (e.g. near-equilibrium reactions), as substrates of the reverse reactions. Co-operative feedback inhibition by threonine and lysine was also included as necessary. We derived the simplest rate equations that describe the salient features of the enzymes in the physiological range of metabolite concentrations in order to incorporate them ultimately into a complete model of the threonine pathway, able to predict quantitatively the behaviour of the pathway under natural or engineered conditions.
Collapse
Affiliation(s)
- C Chassagnole
- INSERM EMI 9929, University Victor Segalen Bordeaux 2, 146 rue Léo Saignat, 33076 Bordeaux, France
| | | | | | | | | |
Collapse
|
15
|
Zhou T, Daugherty M, Grishin NV, Osterman AL, Zhang H. Structure and mechanism of homoserine kinase: prototype for the GHMP kinase superfamily. Structure 2000; 8:1247-57. [PMID: 11188689 DOI: 10.1016/s0969-2126(00)00533-5] [Citation(s) in RCA: 102] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
BACKGROUND Homoserine kinase (HSK) catalyzes an important step in the threonine biosynthesis pathway. It belongs to a large yet unique class of small metabolite kinases, the GHMP kinase superfamily. Members in the GHMP superfamily participate in several essential metabolic pathways, such as amino acid biosynthesis, galactose metabolism, and the mevalonate pathway. RESULTS The crystal structure of HSK and its complex with ADP reveal a novel nucleotide binding fold. The N-terminal domain contains an unusual left-handed betaalphabeta unit, while the C-terminal domain has a central alpha-beta plait fold with an insertion of four helices. The phosphate binding loop in HSK is distinct from the classical P loops found in many ATP/GTP binding proteins. The bound ADP molecule adopts a rare syn conformation and is in the opposite orientation from those bound to the P loop-containing proteins. Inspection of the substrate binding cavity indicates several amino acid residues that are likely to be involved in substrate binding and catalysis. CONCLUSIONS The crystal structure of HSK is the first representative in the GHMP superfamily to have determined structure. It provides insight into the structure and nucleotide binding mechanism of not only the HSK family but also a variety of enzymes in the GHMP superfamily. Such enzymes include galactokinases, mevalonate kinases, phosphomevalonate kinases, mevalonate pyrophosphate decarboxylases, and several proteins of yet unknown functions.
Collapse
Affiliation(s)
- T Zhou
- Department of Biochemistry, University of Texas Southwestern Medical Center, Dallas 75390, USA
| | | | | | | | | |
Collapse
|
16
|
Lee M, Leustek T. Identification of the gene encoding homoserine kinase from Arabidopsis thaliana and characterization of the recombinant enzyme derived from the gene. Arch Biochem Biophys 1999; 372:135-42. [PMID: 10562426 DOI: 10.1006/abbi.1999.1481] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Homoserine kinase (EC 2.7.1.39) catalyzes the formation of O-phospho-l-homoserine, a branch point intermediate in the pathways for Met and Thr in plants. A genomic open reading frame located on the top arm of chromosome II and a corresponding cDNA have been identified from Arabidopsis thaliana that encode homoserine kinase. The HSK gene is composed of an 1113-bp continuous open reading frame that could produce a 38-kDa protein. The gene product has homology with homoserine kinase from bacteria and fungi. It contains a conserved motif, known as GHMP, found in a group of ATP-dependent metabolite kinases and thought to comprise the ATP binding site. The amino-terminal 50 amino acids of the HSK protein show features of a transit peptide for localization to plastids. Genomic blot analysis revealed that there is a single locus in A. thaliana to which the HSK cDNA hybridizes. The HSK protein expressed as a His-tagged construct in Escherichia coli shows a specific activity in an l-homoserine-dependent ADP synthesis assay of 3.09 +/- 0.25 micromol min(-1) mg(-1) protein at pH 8.5 and 37 degrees C. The apparent K(m) values are 0.40 mM for l-homoserine and 0.32 mM for Mg-ATP. Other hydroxylated compounds are not used as substrates. The enzyme requires 40 mM K(+) and 3 mM Mg(2+) for activity. It has an unusually high temperature optimum, yet it is very unstable, losing more than 80% of its activity after a single cycle of freeze-thawing. The HSK enzyme shows no significant regulation by amino acids in vitro.
Collapse
Affiliation(s)
- M Lee
- Biotechnology Center for Agriculture and the Environment and the Plant Science Department, Rutgers University, New Brunswick, New Jersey, 08901-8520, USA
| | | |
Collapse
|
17
|
Blickling S, Beisel HG, Bozic D, Knäblein J, Laber B, Huber R. Structure of dihydrodipicolinate synthase of Nicotiana sylvestris reveals novel quaternary structure. J Mol Biol 1997; 274:608-21. [PMID: 9417939 DOI: 10.1006/jmbi.1997.1393] [Citation(s) in RCA: 67] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
DHDPS is the first enzyme unique to the lysine biosynthetic pathway in plants and bacteria and catalyses the formation of (4S)-4-hydroxy-2,3,4,5-tetrahydro-(2S)-dipicolinic acid. It is feedback-regulated in plants by L-lysine. The crystal structure of Nicotiana sylvestris DHDPS with and without inhibitory lysine bound to the enzyme has been solved to a resolution of 2.8 A. The molecule is a homotetramer composed of a dimer of dimers. Comparison with the structure of Escherichia coli DHDPS showed a novel quaternary structure by a profound rearrangement of the dimers forming the tetramer. The crystal structure of the enzyme in the presence of L-lysine revealed substantial changes. These changes together with the novel quaternary structure provide a structural basis for the strong inhibition of plant DHDPS enzymes by L-lysine.
Collapse
Affiliation(s)
- S Blickling
- Max-Planck-Institut für Biochemie, Abteilung Strukturforschung, Martinsried, Germany
| | | | | | | | | | | |
Collapse
|
18
|
Laber B, Lindell SD, Pohlenz HD. Inactivation of Escherichia coli threonine synthase by DL-Z-2-amino-5-phosphono-3-pentenoic acid. Arch Microbiol 1994; 161:400-3. [PMID: 8042902 DOI: 10.1007/bf00288949] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The rhizocticines and plumbemicines are two groups of di- and tripeptid antibiotics thought to interfere with threonine or threonine-related metabolism. Z-2-amino-5-phosphono-3-pentenoic acid, the common unusual amino acid constituent of the rhizocticines and plumbemicines, was found to irreversibly inhibit Escherichia coli threonine synthase in a time-dependent reaction that followed pseudo-first order and saturation kinetics. These data provide evidence that the toxicity of the rhizocticines and plumbemicines is due to the inhibition of threonine synthase by Z-2-amino-5-phosphone-3-pentenoic acid, which is liberated by peptidases after uptake into the target cell. Additionally, methods for the purification of threonine synthase from an overproducing E. coli strain and for the enzymatic synthesis of L-homoserine phosphate are described.
Collapse
Affiliation(s)
- B Laber
- Schering AG, Berlin, Germany
| | | | | |
Collapse
|
19
|
Martin-Rendon E, Farfán MJ, Ramos C, Calderon IL. Isolation of a mutant allele that deregulates the threonine biosynthesis in Saccharomyces cerevisiae. Curr Genet 1993; 24:465-71. [PMID: 8299165 DOI: 10.1007/bf00351707] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
We have cloned the yeast allele HOM3-R2, that codes for a mutant aspartate kinase which is insensitive to feedback inhibition by threonine, by gap-repair. A strain carrying this allele in a multicopy plasmid, or integrated into the genome, accumulates 14-times and 8-times more threonine than the wild-type, respectively. The sequence of the mutant allele differs from that of the wild-type in a single base pair change, namely a G by an A, at position 1355 in the open reading frame. The fact that the presence of this mutant allele in a cell induces threonine overproduction points to aspartate kinase as the key enzyme in the regulation of threonine biosynthesis in yeast.
Collapse
Affiliation(s)
- E Martin-Rendon
- Departamento de Genética, Facultad de Biología, Universidad de Sevilla, Spain
| | | | | | | |
Collapse
|
20
|
Martínez-Force E, Benítez T. Regulation of aspartate-derived amino acid biosynthesis in the yeastSaccharomyces cerevisiae. Curr Microbiol 1993. [DOI: 10.1007/bf01576262] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
21
|
Ramos C, Calderon IL. Overproduction of threonine by Saccharomyces cerevisiae mutants resistant to hydroxynorvaline. Appl Environ Microbiol 1992; 58:1677-82. [PMID: 1622238 PMCID: PMC195657 DOI: 10.1128/aem.58.5.1677-1682.1992] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
In this work, we isolated and characterized mutants that overproduce threonine from Saccharomyces cerevisiae. The mutants were selected for resistance to the threonine analog alpha-amino-beta-hydroxynorvalerate (hydroxynorvaline), and, of these, the ones able to excrete threonine to the medium were chosen. The mutant strains produce between 15 and 30 times more threonine than the wild type does, and, to a lesser degree, they also accumulate isoleucine. Genetic and biochemical studies have revealed that the threonine overproduction is, in all cases studied, associated with the presence in the strain of a HOM3 allele coding for a mutant aspartate kinase that is totally or partially insensitive to feedback inhibition by threonine. This enzyme seems, therefore, to be crucial in the regulation of threonine biosynthesis in S. cerevisiae. The results obtained suggest that this strategy could be efficiently applied to the isolation of threonine-overproducing strains of yeasts other than S. cerevisiae, even those used industrially.
Collapse
Affiliation(s)
- C Ramos
- Departamento de Genética, Facultad de Biología, Universidad de Sevilla, Spain
| | | |
Collapse
|
22
|
Eikmanns BJ, Metzger M, Reinscheid D, Kircher M, Sahm H. Amplification of three threonine biosynthesis genes in Corynebacterium glutamicum and its influence on carbon flux in different strains. Appl Microbiol Biotechnol 1991; 34:617-22. [PMID: 1369320 DOI: 10.1007/bf00167910] [Citation(s) in RCA: 100] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
The hom-thrB operon (homoserine dehydrogenase/homoserine kinase) and the thrC gene (threonine synthase) of Corynebacterium glutamicum ATCC 13,032 and the homFBR (homoserine dehydrogenase resistant to feedback inhibition by threonine) alone as well as homFBR-thrB operon of C. glutamicum DM 368-3 were cloned separately and in combination in the Escherichia coli/C. glutamicum shuttle vector pEK0 and introduced into different corynebacterial strains. All recombinant strains showed 8- to 20-fold higher specific activities of homoserine dehydrogenase, homoserine kinase, and/or threonine synthase compared to the respective host. In wild-type C. glutamicum, amplification of the threonine genes did not result in secretion of threonine. In the lysine producer C. glutamicum DG 52-5 and in the lysine-plus-threonine producer C. glutamicum DM 368-3 overexpression of hom-thrB resulted in a notable shift of carbon flux from lysine to threonine whereas cloning of homFBR-thrB as well as of homFBR in C. glutamicum DM 368-3 led to a complete shift towards threonine or towards threonine and its precursor homoserine, respectively. Overexpression of thrC alone or in combination with that of homFBR and thrB had no effect on threonine or lysine formation in all recombinant strains tested.
Collapse
Affiliation(s)
- B J Eikmanns
- Institut für Biotechnologie 1, Forschungszentrums Jülich GmbH, Federal Republic of Germany
| | | | | | | | | |
Collapse
|
23
|
Ramos C, Delgado MA, Calderon IL. Inhibition by different amino acids of the aspartate kinase and the homoserine kinase of the yeast Saccharomyces cerevisiae. FEBS Lett 1991; 278:123-6. [PMID: 1847111 DOI: 10.1016/0014-5793(91)80098-n] [Citation(s) in RCA: 28] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
In this paper, we describe a simple method to measure the yeast homoserine kinase and aspartate kinase activities, independently but in the same extract. With this method, we have determined some kinetic parameters for the physiological substrates of both enzymes, and investigated the inhibition exerted by different amino acids on these activities. Of all natural amino acids tested, only threonine inhibits effectively both enzymatic activities, although to a different degree. We did not find the reported inhibition by L-homoserine over the aspartate kinase. Altogether the data point to the aspartate kinase and to the threonine as the key factors in the regulation of this route.
Collapse
Affiliation(s)
- C Ramos
- Departamento de Genética y Biotecnia, Facultad de Biología, Universidad de Sevilla, Spain
| | | | | |
Collapse
|
24
|
Finkelnburg B, Klemme JH. Homoserine kinase from the phototrophic bacteriumRhodospirillum rubrumis not sensitive to feedback inhibition by l-threonine. FEMS Microbiol Lett 1987. [DOI: 10.1111/j.1574-6968.1987.tb02522.x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
|
25
|
Shames SL, Ash DE, Wedler FC, Villafranca JJ. Interaction of aspartate and aspartate-derived antimetabolites with the enzymes of the threonine biosynthetic pathway of Escherichia coli. J Biol Chem 1984. [DOI: 10.1016/s0021-9258(17)42554-3] [Citation(s) in RCA: 30] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
|