Molecular dynamics simulations of substrate channeling through an α–β barrel protein

Abstract Steered molecular dynamics simulations are used to probe the energetics of substrate channeling in an enzyme regulating histidine biosynthesis, imidazole glycerol phosphate synthase (IGPS). IGPS is a multidomain globular protein complex: the glutaminase domain hydrolyzes glutamine to form glutamate and ammonia, and is docked to the cyclase domain, a (β/α)8 barrel protein that completes the ring formation of imidazole glycerol phosphate. Recently, it has been suggested that this protein exploits its barrel structure to channel ammonia from one remote active-site to the other. The current work includes both domains, their substrates, ammonia, and explicit solvent. Compared to the apo-complex, the inclusion of substrates does indeed affect the barrier to ammonia entry into the channel as well its transport through the barrel. Based on bioinformatic data, we suggest an “open-gate” mechanism that has a low barrier to ammonia entry. We also perform the first systematic investigation of interface water molecules near the channel gate and argue that the optimum number of water molecules inside the channel is one.

[1]  A. Mattevi,et al.  Cross-talk and ammonia channeling between active centers in the unexpected domain arrangement of glutamate synthase. , 2000, Structure.

[2]  Janet L. Smith,et al.  The crystal structure of GMP synthetase reveals a novel catalytic triad and is a structural paradigm for two enzyme families , 1996, Nature Structural Biology.

[3]  A. McDermott,et al.  The time scale of the catalytic loop motion in triosephosphate isomerase. , 2001, Journal of molecular biology.

[4]  J. L. Smith,et al.  Coupled formation of an amidotransferase interdomain ammonia channel and a phosphoribosyltransferase active site. , 1997, Biochemistry.

[5]  K. Schulten,et al.  Free energy calculation from steered molecular dynamics simulations using Jarzynski's equality , 2003 .

[6]  G. Hummer,et al.  Free energy reconstruction from nonequilibrium single-molecule pulling experiments , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[7]  B. Roux,et al.  Free energy profiles for H+ conduction along hydrogen-bonded chains of water molecules. , 1998, Biophysical journal.

[8]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[9]  P. Babbitt,et al.  Divergent evolution of enzymatic function: mechanistically diverse superfamilies and functionally distinct suprafamilies. , 2001, Annual review of biochemistry.

[10]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[11]  V. Davisson,et al.  Imidazole glycerol phosphate synthase: the glutamine amidotransferase in histidine biosynthesis. , 1993, Biochemistry.

[12]  Cathy H. Wu,et al.  UniProt: the Universal Protein knowledgebase , 2004, Nucleic Acids Res..

[13]  C. Orengo,et al.  One fold with many functions: the evolutionary relationships between TIM barrel families based on their sequences, structures and functions. , 2002, Journal of molecular biology.

[14]  C. Matthews,et al.  Proline replacements and the simplification of the complex, parallel channel folding mechanism for the alpha subunit of Trp synthase, a TIM barrel protein. , 2003, Journal of molecular biology.

[15]  R. Callender,et al.  Active site loop motion in triosephosphate isomerase: T-jump relaxation spectroscopy of thermal activation. , 2003, Biochemistry.

[16]  R. Wierenga,et al.  The TIM‐barrel fold: a versatile framework for efficient enzymes , 2001, FEBS letters.

[17]  Sridar V. Chittur,et al.  Crystal Structure of Imidazole Glycerol Phosphate Synthase: A Tunnel through a (β/α)8 Barrel Joins Two Active Sites , 2001 .

[18]  L. Tong,et al.  Solution-state NMR investigations of triosephosphate isomerase active site loop motion: ligand release in relation to active site loop dynamics. , 2001, Journal of molecular biology.

[19]  S. Schuster,et al.  Three-dimensional structure of Escherichia coli asparagine synthetase B: a short journey from substrate to product. , 1999, Biochemistry.

[20]  V Jo Davisson,et al.  Substrate-induced changes in the ammonia channel for imidazole glycerol phosphate synthase. , 2003, Biochemistry.

[21]  C. Jarzynski Equilibrium free-energy differences from nonequilibrium measurements: A master-equation approach , 1997, cond-mat/9707325.

[22]  K. Schulten,et al.  Calculating potentials of mean force from steered molecular dynamics simulations. , 2004, The Journal of chemical physics.

[23]  H. Dyson,et al.  Structural basis for Hif-1α/CBP recognition in the cellular hypoxic response , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[24]  R. Sterner,et al.  Imidazole glycerol phosphate synthase from Thermotoga maritima. Quaternary structure, steady-state kinetics, and reaction mechanism of the bienzyme complex. , 2001, The Journal of biological chemistry.

[25]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[26]  J Hermans,et al.  Hydrophilicity of cavities in proteins , 1996, Proteins.

[27]  Perforation of the tunnel wall in carbamoyl phosphate synthetase derails the passage of ammonia between sequential active sites. , 2004 .

[28]  M Wilmanns,et al.  Structural evidence for evolution of the beta/alpha barrel scaffold by gene duplication and fusion. , 2000, Science.

[29]  H M Holden,et al.  The small subunit of carbamoyl phosphate synthetase: snapshots along the reaction pathway. , 1999, Biochemistry.

[30]  Peter G Wolynes,et al.  Role of water mediated interactions in protein-protein recognition landscapes. , 2003, Journal of the American Chemical Society.

[31]  V Jo Davisson,et al.  Toward understanding the mechanism of the complex cyclization reaction catalyzed by imidazole glycerolphosphate synthase: crystal structures of a ternary complex and the free enzyme. , 2003, Biochemistry.

[32]  Matthias Wilmanns,et al.  Structural evidence for ammonia tunneling across the (beta alpha)(8) barrel of the imidazole glycerol phosphate synthase bienzyme complex. , 2002, Structure.

[33]  F. Raushel,et al.  The amidotransferase family of enzymes: molecular machines for the production and delivery of ammonia. , 1999, Biochemistry.

[34]  I. Tinoco,et al.  Equilibrium Information from Nonequilibrium Measurements in an Experimental Test of Jarzynski's Equality , 2002, Science.

[35]  C. Jarzynski Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[36]  Rommie E. Amaro,et al.  Developing an energy landscape for the novel function of a (β/α)8 barrel: Ammonia conduction through HisF , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[37]  F. Ritort,et al.  Bias and error in estimates of equilibrium free-energy differences from nonequilibrium measurements , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[38]  M Wilmanns,et al.  Directed evolution of a (beta alpha)8-barrel enzyme to catalyze related reactions in two different metabolic pathways. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[39]  A Teplyakov,et al.  Channeling of ammonia in glucosamine-6-phosphate synthase. , 2001, Journal of molecular biology.

[40]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[41]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[42]  K. Schulten,et al.  Energetics of glycerol conduction through aquaglyceroporin GlpF , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[43]  Structure of HisF, a histidine biosynthetic protein from Pyrobaculum aerophilum. , 2001, Acta crystallographica. Section D, Biological crystallography.

[44]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[45]  I. Rayment,et al.  Structure of carbamoyl phosphate synthetase: a journey of 96 A from substrate to product. , 1997, Biochemistry.

[46]  K. Schulten,et al.  Control of the Selectivity of the Aquaporin Water Channel Family by Global Orientational Tuning , 2002, Science.

[47]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[48]  L. Tong,et al.  Optimal alignment for enzymatic proton transfer: Structure of the Michaelis complex of triosephosphate isomerase at 1.2-Å resolution , 2002, Proceedings of the National Academy of Sciences of the United States of America.