|
|
||||||||
1 Department of Biomedicine, University of Bergen, 5009-Bergen, Norway
2 Department of Biochemistry, Universidade de Coimbra, 3000 Coimbra, Portugal
3 Departamento de Química Física, Facultad de Ciencias, 18071 Granada, Spain
Reprint requests to: Aurora Martínez, Department of Biomedicine, Jonas Lies vei 91, 5009-Bergen, Norway; e-mail: aurora.martinez{at}ibmb.uib.no; fax: 47-55586360.
(RECEIVED December 28, 2003; FINAL REVISION February 4, 2004; ACCEPTED February 9, 2004)
| Abstract |
|---|
|
|
|---|
Keywords: phosphorylation; phenylalanine hydroxylase; molecular dynamics simulations; conformational stability; electrostatic interactions; proteolysis
Abbreviations: BH4, (6R)-L-erythro-5,6,7,8-tetrahydrobiopterin CD, circular dichroism DSC, differential scanning calorimetry MD, molecular dynamics PAH, phenylalanine hydroxylase PKA, cyclic AMPdependent protein kinase wt-PAH, wild-type phenylalanine hydroxylase.
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.03595904.
| Introduction |
|---|
|
|
|---|
-helical content and resistance to tryptic proteolysis (Miranda et al. 2002). NMR spectroscopy has also shown that the mobile N-terminal region of PAH shows an increased ordered secondary structure upon phosphorylation (Horne et al. 2002). FRET measurements and molecular dynamics (MD) simulations performed on a 10-residue peptide around the MAP kinase substrate Ser-31 in the homologous enzyme tyrosine hydroxylase also indicate that phosphorylation causes the peptide backbone to adopt a compact structure away from the active site (Stultz et al. 2002). This change also appears to be initiated by the favorable electrostatic interaction established by the phosphoserine and the side chain of a neighboring arginine. Similarly, phosphorylation at Ser 16 of phospholamban by PKA promotes a coil-to-helix transition that modulates the structural coupling between the transmembrane and cytosolic domains of the protein, most probably through an electrostatic linkage between the phosphate group and the kinase-recognizing residue Arg13 (Li et al. 2003). Thus, it appears that for many proteins the regulatory effects exerted by phosphorylation rely on a conformational switch put forward by the formation of salt bridges between the phosphate group and adjacent positively charged residues. In this work we have performed MD simulations in the modeled structure of the full-length unphosphorylated and phosphorylated enzyme in order to investigate if the rearrangement of the N-terminal end brought about by the electrostatic interaction between the phosphate and Arg13 is built up long range to a global conformational change in the rest of the enzyme. The stabilizing and destabilizing interactions in the resulting structures were analyzed by our implementation of the Tanford-Kirkwood model (Ibarra-Molero et al. 1999a; Sanchez-Ruiz and Makhatadze 2001). This is admittedly a very simple model, but it has been shown to predict qualitatively the effect of charge-reversal and charge-deletion mutations on the stability of several proteins (Sanchez-Ruiz and Makhatadze 2001). In addition, the structural information obtained from the MD simulated conformers has been validated by differential scanning calorimetry (DSC), thermal-dependent circular dichroism (CD), and fluorescence spectroscopy using unphosphorylated and phosphorylated human PAH and the mutants S16D-PAH and S16E-PAH.
| Results |
|---|
|
|
|---|
atoms of Arg13 and the C
of Glu381 at the catalytic domain (Fig. 1A
|
|
Wi
values, due to interactions mainly with groups of the opposite charge (Fig. 3
Wq - q
values are negative for both enzyme forms, namely, 115 and 126 kJ/mole for the unphosphorylated and phosphorylated subunit PAH, respectively.
Wq - q
is considered to be an approximation of the chargecharge contributions to the unfolding Gibbs energies (
Gq - q), and our results indicate that the phosphorylated form is more stable. Nevertheless, a substantial contribution to the 
Wq - q
=16 kJ/mole per subunit between the phosphorylated and unphosphorylated form seems to arise mainly from the favorable interactions that Arg13 and Lys14 establish with the phosphate at Ser16 (Fig. 3B
|
|
-helical content and a concomitant decrease of the unordered structure in the N-terminal tail (Fig. 2
|
|
|
| Discussion |
|---|
|
|
|---|
-helical structure and an increased affinity for the substrate (Miranda et al. 2002). On the other hand, the cooperative binding of the substrate certainly induces large conformational changes affecting the secondary, tertiary, and quaternary structure of the enzyme (Kaufman 1993; Thórólfsson et al. 2003).
Several proteins, such as the botulinum neurotoxins, have been shown to increase their stability upon phosphorylation (Encinar et al. 1998; Blanes-Mira et al. 2001). Phosphory-lated PAH, however, does not show increased thermal stability, despite the estimated 
Wq - q
=16 kJ/mole per subunit between the phosphorylated and unphosphorylated form. This discrepancy may reflect the limitations of the simple electrostatic model used, such as the facts that denatured-state interactions are not taken into account and only chargecharge interactions are considered (i.e., charge-solvation effects, which may be important, are not included). This notwithstanding, our simple electrostatic calculations do seem to be consistent with the electrostatic effect of phosphorylation being local (the significant change upon phosphorylation in the calculated chargecharge interaction energy arises mainly from the interaction between the phosphate and the nearby residues Arg13 and Lys14). Thus, the calculations do not suggest large conformational changes upon phosphorylation and, in this sense, seem consistent with the similar thermal stability of unphosphorylated and phosphorylated PAH, and their similar rate of proteolytic degradation by chymotrypsin. In this context, the increased resistance toward limited tryptic proteolysis exhibited by the phosphorylated form appears noteworthy. Our earlier results have shown that phosphorylation protects toward the specific cleavage at the carboxyl site of Arg13 most probably due to the inability of the protease to bind to this recognition site when it is in turn interacting with the phosphoserine (Døskeland et al. 1996; Miranda et al. 2002). This inability seems to decelerate the overall tryptic proteolysis rate in the rest of the protein. It is tentative to speculate on the physiological implications for a specific protection of phosphorylation toward degradation and turnover in vivo by proteases with trypsin-like specificity. The protection against cellular degradation shown by both phosphorylated ATF2 transcription factor (Fuchs et al. 2000) and the p27 (Kip1) protein (Ishida et al. 2000) has been associated to a protection from ubiquitination and decrease of degradation by the ubiquitin-dependent proteasome pathway. Although it has been established that human PAH is a target for ubiquitination by ubiquitin-conjugating enzyme system isolated from rat liver (Døskeland and Flatmark 2001), nothing is known about the effects of phosphorylation on ubiquitination of the enzyme and vice versa. Other studies have linked the higher intracellular stability of phosphorylated proteins, for example,
II spectin (Nicolas et al. 2002) to its protection toward the nonlysosomal neutral thiol protease calpain with partial trypsin-like cleavage specificity (Yajima and Kawashima 2002). It seems plausible that the local specific conformational change induced by phosphorylation of Ser16 at the autoregulatory sequence has consequences on the turnover of the enzyme in vivo, in addition to the short-term activation (~2.4 increased catalytic efficiency for the basal activity of the enzyme) and the synergy with L-Phe activation. With respect to this activation, the role of the local conformational change increasing the accessibility of the active site might be reinforced by the favorable effect that phosphorylation has on the energy of the chargecharge interactions of Arg270 and His285, residues directly involved on the binding of the substrate. A direct effect of phosphorylation on His285 can also be a way to positively modulate the reactivity of the coordinated nonheme iron and activate the enzyme, which would explain the higher catalytic efficiency of the phosphorylated enzyme even prior to activation by incubation with L-Phe (Miranda et al. 2002). In this context, it is interesting to note that two of the three iron-coordinated water molecules in unphosphorylated resting PAH (PDB 2PHM
[PDB]
; Kobe et al. 1999) are not observed in the structure of the phosphorylated enzyme (PDB 1PHZ
[PDB]
), which resembles the iron-coordination in the substrate (and cofactor) bound ternary complexes of PAH (Andersen et al. 2002).
| Materials and methods |
|---|
|
|
|---|
Limited proteolysis by trypsin and chymotrypsin
For the limited proteolysis, either TPCK-treated trypsin from bovine pancreas (Sigma) or type II chymotrypsin (Sigma) was used at a protease/PAH ratio of 0 : 1, 0.002 : 1, 0.005 : 1, 0.01 : 1, and 0.1 : 1, in a final volume of 45 µL. The concentration of the various forms of PAH was 0.1 µg/µL. The reaction was performed for 20 min at 30°C and quenched with SDS denaturation buffer (3-min treatment at 95°C). The samples were then analyzed by SDS-PAGE (performed on 10% [w/v] acrylamide gels at 15 mA/ gel). The gels were stained by Coomassie brilliant blue, dried, scanned, and further analyzed using the software Deskscan II (Hewlett-Packard Co.) and Phoretix 1D Plus (Nonlinear Dynamics Ltd).
Circular dichroism
CD measurements were performed on a Jasco J-810 spectropolarimeter equipped with a Jasco Peltier 423S element, for temperature control. Purified samples of the enzyme were prepared on a degassed solution of 20 mM Na-phosphate and 0.15 M KF (pH 7.0) at a concentration of 10 µM. Equimolar amounts of ferrous ammonium sulfate per enzyme subunit were added in order to convert some iron-free apoenzyme present in the samples into holoenzyme (Chehin et al. 1998; Thórólfsson et al. 2002). Quartz cells of 1-mm path length were used. Thermal denaturation was monitored by following the changes in ellipticity at 222 nm, at a scan rate of 0.7 K/min in the temperature range 30°C to 70°C. Analysis of the data was performed by using the Standard Analysis program provided with the instrument.
Fluorescence measurements
Measurements were performed on a Perkin-Elmer LS-50B luminescence spectrometer with a constant temperature cell holder, using 0.5-cm path-length quartz cells. Enzyme samples (1 µM) were prepared in 20 mM Na-Hepes, 0.2 M NaCl (pH 7.0). The excitation wavelength was 295 nm, and the excitation and emission slits were three and five, respectively. All spectra were corrected for blank emission.
Differential scanning calorimetry
Measurements were performed on a MicroCal VP-DSC differential scanning calorimeter (MicroCal Inc.) with cell volumes of 0.5 mL at the indicated scan rates. A 100 mM Na-Hepes buffer, 0.1 M NaCl (pH 7.0) was used in all experiments. Calorimetric cells were kept under an excess pressure of 30 psi to prevent degassing during the scan. Purified tetrameric PAH enzymes (30 µM per subunit) with equimolar amounts of ferrous ammonium sulfate per subunit were used. Further details about the calorimetric experiment and the calorimetric data processing can be found elsewhere (Ibarra-Molero et al. 1999b; Thórólfsson et al. 2002).
MD simulations
MD simulations were performed in the dimeric form of the modeled full-length unphosphorylated and phosphorylated PAH at Ser16 (Miranda et al. 2002), previously prepared by molecular docking of the 18 N-terminal residues into the crystal structure of rat PAH (PDB codes 2PHM
[PDB]
and 1PHZ
[PDB]
; Kobe et al. 1999), which lacks structural information about this N-terminal region. The structures were solvated in a water box, using the TIP3P model for water molecules extended at least 6 Å in each direction from the solute (25,000 water molecules included), and Na+ and counter ions were added by using the LEAP module of the AMBER 7 program (parm94 force field included). The system was first energy-minimized and then heated for 10 psec to 300 K, and the Particle Mesh Ewald method was used to calculate long-range electrostatic interactions. Following the heating step, the system was maintained at 300 K, and the MD simulation was computed at 1.0-fsec intervals for 600 psec, with frames collected every 5.0 psec. The nonbonded interaction cutoff was set to 8.0 Å. From each of the frames of the unphosphorylated and phosphorylated structures, atomic coordinates were collected for comparison.
Chargecharge interactions calculations
Calculations of the energies of chargecharge interactions were carried by using our implementation of the Tanford-Kirkwood model with the solvent accessibility correction of Gurd, as we have previously described in detail (Ibarra-Molero et al. 1999a; Sanchez-Ruiz and Makhatadze 2001). Dielectric constants of 4.0 and 78.5 were used for the protein and the aqueous solvent, respectively. In this work we have focused on the energy due to the chargecharge interaction of group i with the rest of the ionizable groups in the protein
Wi
, which can be used to estimate the total chargecharge interaction energy in the proteins (
Wq - q
):
![]() | (1) |
| Acknowledgments |
|---|
The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
| References |
|---|
|
|
|---|
Blanes-Mira, C., Ibañez, C., Fernandez-Ballester, G., Planells-Cases, R., Perez-Paya, E., and Ferrer-Montiel, A. 2001. Thermal stabilization of the catalytic domain of botulinum neurotoxin E by phosphorylation of a single tyrosine residue. Biochemistry 40: 22342242.[CrossRef][Medline]
Chehin, R., Thórólfsson, M., Knappskog, P.M., Martínez, A., Flatmark, T., Arrondo, J.L., and Muga, A. 1998. Domain structure and stability of human phenylalanine hydroxylase inferred from infrared spectroscopy. FEBS Lett. 422: 225230.[CrossRef][Medline]
Døskeland, A.P. and Flatmark, T. 2001. Conjugation of phenylalanine hydroxylase with polyubiquitin chains catalysed by rat liver enzymes. Biochim. Biophys. Acta 1547: 379386.[Medline]
Døskeland, A.P., Martínez, A., Knappskog, P.M., and Flatmark, T. 1996. Phosphorylation of recombinant human phenylalanine hydroxylase: Effect on catalytic activity, substrate activation and protection against non- specific cleavage of the fusion protein by restriction protease. Biochem. J. 313: 409414.
Encinar, J.A., Fernandez, A., Ferragut, J.A., Gonzalez-Ros, J.M., DasGupta, B.R., Montal, M., and Ferrer-Montiel, A. 1998. Structural stabilization of botulinum neurotoxins by tyrosine phosphorylation. FEBS Lett. 429: 7882.[CrossRef][Medline]
Erlandsen, H. and Stevens, R.C. 1999. The structural basis of phenylketonuria. Mol. Genet. Metab. 68: 103125.[CrossRef][Medline]
Fitzpatrick, P.F. 2000. The aromatic amino acid hydroxylases. Adv. Enzymol. Relat. Areas Mol. Biol. 74: 235294.[CrossRef][Medline]
Flatmark, T. and Stevens, R.C. 1999. Structural insight into the aromatic amino acid hydroxylases and their disease-related mutant forms. Chem. Rev. 99: 21372160.[CrossRef][Medline]
Fuchs, S.Y., Tappin, I., and Ronai, Z. 2000. Stability of the ATF2 transcription factor is regulated by phosphorylation and dephosphorylation. J. Biol. Chem. 275: 1256012564.
Graves, J.D. and Krebs, E.G. 1999. Protein phosphorylation and signal transduction. Pharmacol. Ther. 82: 111121.[CrossRef][Medline]
Horne, J., Jennings, I.G., Teh, T., Gooley, P.R., and Kobe, B. 2002. Structural characterization of the N-terminal autoregulatory sequence of phenylalanine hydroxylase. Protein Sci. 11: 20412047.
Ibarra-Molero, B., Loladze, V.V., Makhatadze, G.I., and Sanchez-Ruiz, J.M. 1999a. Thermal versus guanidine-induced unfolding of ubiquitin: An analysis in terms of the contributions from chargecharge interactions to protein stability. Biochemistry 38: 81388149.[CrossRef][Medline]
Ibarra-Molero, B., Makhatadze, G.I., and Sanchez-Ruiz, J.M. 1999b. Cold denaturation of ubiquitin. Biochim. Biophys. Acta 1429: 384390.[CrossRef][Medline]
Ishida, N., Kitagawa, M., Hatakeyama, S., and Nakayama, K. 2000. Phosphorylation at serine 10, a major phosphorylation site of p27(Kip1), increases its protein stability. J. Biol. Chem. 275: 2514625154.
Johnson, L.N. 1992. Glycogen phosphorylase: Control by phosphorylation and allosteric effectors. FASEB J. 6: 22742282.[Abstract]
Johnson, L.N. and Lewis, R.J. 2001. Structural basis for control by phosphorylation. Chem. Rev. 101: 22092242.[CrossRef][Medline]
Kaufman, S. 1993. The phenylalanine hydroxylating system. Adv. Enzymol. Relat. Areas Mol. Biol. 67: 77264.[Medline]
Knappskog, P.M. and Haavik, J. 1995. Tryptophan fluorescence of human phenylalanine hydroxylase produced in Escherichia coli. Biochemistry 34: 1179011799.[CrossRef][Medline]
Knappskog, P.M. and Martínez, A. 1997. Effect of mutations at Cys237 on the activation state and activity of human phenylalanine hydroxylase. FEBS Lett. 409: 711.[Medline]
Kobe, B., Jennings, I.G., House, C.M., Michell, B.J., Goodwill, K.E., Santarsiero, B.D., Stevens, R.C., Cotton, R.G., and Kemp, B.E. 1999. Structural basis of autoregulation of phenylalanine hydroxylase. Nat. Struct. Biol. 6: 442448.[CrossRef][Medline]
Leach, A.R. and Kuntz, I.D. 1992. Conformational-analysis of flexible ligands in macromolecular receptor-sites. J. Comput. Chem. 13: 730748.[CrossRef]
Li, J., Bigelow, D.J., and Squier, T.C. 2003. Phosphorylation by cAMP-dependent protein kinase modulates the structural coupling between the transmembrane and cytosolic domains of phospholamban. Biochemistry 42: 1067410682.[CrossRef][Medline]
Lin, K., Rath, V.L., Dai, S.C., Fletterick, R.J., and Hwang, P.K. 1996. A protein phosphorylation switch at the conserved allosteric site in GP. Science 273: 15391542.[Abstract]
Martínez, A., Knappskog, P.M., Olafsdottir, S., Døskeland, A.P., Eiken, H.G., Svebak, R.M., Bozzini, M., Apold, J., and Flatmark, T. 1995. Expression of recombinant human phenylalanine hydroxylase as fusion protein in Escherichia coli circumvents proteolytic degradation by host cell proteases: Isolation and characterization of the wild-type enzyme. Biochem. J. 306: 589597.
Miranda, F.F., Teigen, K., Thórólfsson, M., Svebak, R.M., Knappskog, P.M., Flatmark, T., and Martínez, A. 2002. Phosphorylation and mutations of Ser(16) in human phenylalanine hydroxylase: Kinetic and structural effects. J. Biol. Chem. 277: 4093740943.
Nicolas, G., Fournier, C.M., Galand, C., Malbert-Colas, L., Bournier, O., Kroviarski, Y., Bourgeois, M., Camonis, J.H., Dhermy, D., Grandchamp, B., et al. 2002. Tyrosine phosphorylation regulates
II spectrin cleavage by calpain. Mol. Cell. Biol. 22: 35273536.
Phillips, R.S., Parniak, M.A., and Kaufman, S. 1984. Spectroscopic investigation of ligand interaction with hepatic phenylalanine hydroxylase: Evidence for a conformational change associated with activation. Biochemistry 23: 38363842.[CrossRef][Medline]
Russo, A.A., Jeffrey, P.D., and Pavletich, N.P. 1996. Structural basis of cyclin-dependent kinase activation by phosphorylation. Nat. Struct. Biol. 3: 696700.[CrossRef][Medline]
Sanchez-Ruiz, J.M. and Makhatadze, G.I. 2001. To charge or not to charge? Trends Biotechnol. 19: 132135.[CrossRef][Medline]
Scriver, C.R., Eisensmith, R.C., Woo, S.L.C., and Kaufman, S. 1995. The hyperphenylalaninemias in man and mouse. Annu. Rev. Genet. 28: 141165.[CrossRef]
Stultz, C.M., Levin, A.D., and Edelman, E.R. 2002. Phosphorylation-induced conformational changes in a mitogen-activated protein kinase substrate: Implications for tyrosine hydroxylase activation. J. Biol. Chem. 277: 4765347661.
Teigen, K., Frøystein, N.Å., and Martínez, A. 1999. The structural basis of the recognition of phenylalanine and pterin cofactors by phenylalanine hydroxylase: Implications for the catalytic mechanism. J. Mol. Biol. 294: 807823.[CrossRef][Medline]
Thórólfsson, M., Ibarra-Molero, B., Fojan, P., Petersen, S.B., Sanchez-Ruiz, J.M., and Martínez, A. 2002. L-Phenylalanine binding and domain organization in human phenylalanine hydroxylase: A differential scanning calorimetry study. Biochemistry 41: 75737585.[CrossRef][Medline]
Thórólfsson, M., Teigen, K., and Martínez, A. 2003. Activation of phenylalanine hydroxylase: Effect of substitutions at Arg68 and Cys237. Biochemistry 42: 34193428.[CrossRef][Medline]
Wretborn, M., Humble, E., Ragnarsson, U., and Engstrom, L. 1980. Amino acid sequence at the phosphorylated site of rat liver phenylalanine hydroxylase and phosphorylation of a corresponding synthetic peptide. Biochem. Biophys. Res. Commun. 93: 403408.[CrossRef][Medline]
Yajima, Y. and Kawashima, S. 2002. Calpain function in the differentiation of mesenchymal stem cells. Biol. Chem. 383: 757764.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |