In silico thermodynamics stability change analysis involved in BH4 responsive mutations in phenylalanine hydroxylase: QM/MM and MD simulations analysis

The mammalian tetrahydrobiopterin (BH4)-dependent phenylalanine hydroxylases (PAH), involved in important metabolic pathways of phenylalanine, belong to non-heme iron-containing aromatic acid hydroxylases’ enzyme (AAH) family. AAHs utilize BH4 as protein co-factor and thus promote hydroxylation reactions of their substrates. Any alterations in BH4 -mediated AAH’s pathway or mutations in these enzymes are responsible for various disorders, and thus highlights the importance of mutational analysis to assess the effect on their biosynthetic pathways. Our present studies are aimed at single-site mutations in PAH that lead to thermodynamic stability change upon folding and further validation of designed non-reduced BH2 designed co-factors. We have presented single-site mutational analysis of PAH where single-site mutations have been identified from known literature. Further, in silico studies with the PAH, in silico mutant PAH, and crystallized known mutant A313T forms, involved QM/MM and Molecular Dynamics (MD) simulations analysis. The modified co-factor A showed high affinity with PAH and all mutant PAH with high G-score of −14.851. The best pose high affinity co-factor A subjected to QM/MM optimization which leads to square-pyramidal coordination of non-heme active site. The structural and energetic information obtained from the production phase of 20 ns MD simulation of co-factor-metalloprotein complex results helped to understand the binding mode and involvement of three molecules throughout the reaction pathways’ catalysis of PAH. The free energies of binding (dG) of A were found to be −68.181 kcal/mol and −72.249 for 1DMW and 1TDW for A313T mutant. Binding of Co-factor A do not perturb the coordination environment of iron at the active site which resides in 2-Histdine and 1-Glutamate triad, and may enhance the percentage response towards co-factor-mediated therapy.

[1]  G. K. Smith,et al.  Biosynthesis and metabolism of tetrahydrobiopterin and molybdopterin. , 1985, Annual review of biochemistry.

[2]  T. Sugimoto,et al.  Stereochemistry of biopterin cofactor and facile methods for the determination of the stereochemistry of a biologically active 5,6,7,8-tetrahydropterin. , 1985, Journal of biochemistry.

[3]  J. Haavik,et al.  Tyrosine hydroxylase and Parkinson's disease , 1998, Molecular Neurobiology.

[4]  R. Stevens,et al.  Correction of kinetic and stability defects by tetrahydrobiopterin in phenylketonuria patients with certain phenylalanine hydroxylase mutations. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[5]  P. Fitzpatrick.  Mechanism of aromatic amino acid hydroxylation. , 2003, Biochemistry.

[6]  S. Kaufman.  A new cofactor required for the enzymatic conversion of phenylalanine to tyrosine. , 1958, The Journal of biological chemistry.

[7]  Isopenicillin N Synthase: Mechanistic Studies , 1991 .

[8]  N. Holtzman,et al.  Phenylketonuria due to a deficiency of dihydropteridine reductase. , 1975, The New England journal of medicine.

[9]  S. Dahlén,et al.  Leukotrienes and lipoxins: structures, biosynthesis, and biological effects. , 1987, Science.

[10]  S. Gersting,et al.  Phenylketonuria as a model for protein misfolding diseases and for the development of next generation orphan drugs for patients with inborn errors of metabolism , 2010, Journal of Inherited Metabolic Disease.

[11]  A. Ponzone,et al.  Atypical phenylketonuria with “dihydrobiopterin synthetase” deficiency: Absence of phosphate-eliminating enzyme activity demonstrated in liver , 1985, European Journal of Pediatrics.

[12]  Philippe Bogaerts,et al.  Fast and accurate predictions of protein stability changes upon mutations using statistical potentials and neural networks: PoPMuSiC-2.0 , 2009, Bioinform..

[13]  N. Blau,et al.  The metabolic and molecular bases of tetrahydrobiopterin-responsive phenylalanine hydroxylase deficiency. , 2004, Molecular genetics and metabolism.

[14]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[15]  R. Ho,et al.  Dioxygen Activation by Enzymes with Mononuclear Non-Heme Iron Active Sites. , 1996, Chemical reviews.

[16]  L. Que,et al.  Dioxygen activation at mononuclear nonheme iron active sites: enzymes, models, and intermediates. , 2004, Chemical reviews.

[17]  J. Klinman,et al.  Oxygen-18 kinetic isotope effect studies of the tyrosine hydroxylase reaction: Evidence of rate limiting oxygen activation , 1998 .

[18]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[19]  S. Udenfriend.  B. TYROSINE HYDROXYLASE , 1966 .

[20]  M. W. van der Kamp,et al.  Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. , 2013, Biochemistry.

[21]  C. Townsend,et al.  Purification and characterization of clavaminate synthase from Streptomyces clavuligerus: an unusual oxidative enzyme in natural product biosynthesis. , 1990, Biochemistry.

[22]  R. Stevens,et al.  Crystal structure and site-specific mutagenesis of pterin-bound human phenylalanine hydroxylase. , 2000, Biochemistry.

[23]  N. Blau,et al.  GTP cyclohydrolase I deficiency, a new enzyme defect causing hyperphenylalaninemia with neopterin, biopterin, dopamine, and serotonin deficiencies and muscular hypotonia , 1984, European Journal of Pediatrics.

[24]  L. Que,et al.  Oxygen activating nonheme iron enzymes. , 1998, Current opinion in chemical biology.

[25]  Circular Dichroism and Magnetic Circular Dichroism Spectroscopy of the Catalytically Competent Ferrous Active Site of Phenylalanine Hydroxylase and Its Interaction with Pterin Cofactor , 1999 .

[26]  C. Townsend,et al.  Elucidation of the order of oxidations and identification of an intermediate in the multistep clavaminate synthase reaction. , 1991, Biochemistry.

[27]  D. Ballou,et al.  Purification and characterization of phthalate oxygenase and phthalate oxygenase reductase from Pseudomonas cepacia. , 1987, The Journal of biological chemistry.

[28]  Akash Khandelwal,et al.  A combination of docking, QM/MM methods, and MD simulation for binding affinity estimation of metalloprotein ligands. , 2005, Journal of medicinal chemistry.

[29]  Barry L. Stoddard,et al.  The 2.1-A resolution structure of iron superoxide dismutase from Pseudomonas ovalis. , 1990, Biochemistry.

[30]  T. J. Kappock,et al.  Pterin-Dependent Amino Acid Hydroxylases. , 1996, Chemical reviews.

[31]  S. Kaufman.  New tetrahydrobiopterin-dependent systems. , 1993, Annual review of nutrition.

[32]  A. Roscher,et al.  Tetrahydrobiopterin as an alternative treatment for mild phenylketonuria. , 2002, The New England journal of medicine.

[33]  Andrew L. Feig,et al.  Reactions of Non-Heme Iron(II) Centers with Dioxygen in Biology and Chemistry , 1994 .

[34]  J. Haavik,et al.  A structural approach into human tryptophan hydroxylase and its implications for the regulation of serotonin biosynthesis. , 2001, Current medicinal chemistry.

[35]  D. Grahame-Smith.  Tryptophan hydroxylation in brain. , 1964, Biochemical and biophysical research communications.

[36]  C. A. Nichol,et al.  Biopterin cofactor biosynthesis: GTP cyclohydrolase, neopterin and biopterin in tissues and body fluids of mammalian species. , 1984, Life sciences.

[37]  B. Eser,et al.  Single turnover kinetics of tryptophan hydroxylase: evidence for a new intermediate in the reaction of the aromatic amino acid hydroxylases. , 2010, Biochemistry.