|Home | About | Journals | Submit | Contact Us | Français|
The structural basis for allosteric regulation of phenylalanine hydroxylase (PAH), whose dysfunction causes phenylketonuria (PKU), is poorly understood. A new morpheein model for PAH allostery is proposed to consist of a dissociative equilibrium between two architecturally different tetramers whose interconversion requires a ~90° rotation between the PAH catalytic and regulatory domains, the latter of which contains an ACT domain. This unprecedented model is supported by in vitro data on purified full length rat and human PAH. The conformational change is both predicted to and shown to render the tetramers chromatographically separable using ion exchange methods. One novel aspect of the activated tetramer model is an allosteric phenylalanine binding site at the inter-subunit interface of ACT domains. Amino acid ligand-stabilized ACT domain dimerization follows the multimerization and ligand binding behavior of ACT domains present in other proteins in the PDB. Spectroscopic, chromatographic, and electrophoretic methods demonstrate a PAH equilibrium consisting of two architecturally distinct tetramers as well as dimers. We postulate that PKU-associated mutations may shift the PAH quaternary structure equilibrium in favor of the low activity assemblies. Pharmacological chaperones that stabilize the ACT:ACT interface can potentially provide PKU patients with a novel small molecule therapeutic.
Diminished activity of human phenylalanine hydroxylase (PAH), which catalyzes the iron, dioxygen, and BH4 dependent oxidation of phenylalanine (Phe) to tyrosine, is the most common inborn error of metabolism; the most severe cases result in phenylketonuria (PKU) . Despite the success of dietary therapy in preventing severe PKU-associated mental deficits, patients continue to struggle with lifelong control of blood Phe levels. Poor control can compromise behavior, learning, and executive function [2–3]. Thus, alternative therapies are sought, including small molecule therapeutics [1, 4]. Pharmacological chaperones that stabilize PAH structure are gaining support as a long term therapeutic approach [5–12]. The cofactor BH4, which is an active site directed pharmacological chaperone, is in clinical use with some success . Although biological and chemical chaperones are generally viewed as assisting in all aspects of protein folding, a focus on quaternary structure dynamics provides a different perspective that is applicable to PAH. Pharmacological chaperones can function allosterically (e.g. not at the enzyme active site) to assist a protein in achieving and/or maintaining its most active structure. PAH is known to be an allosteric enzyme where allostery involves major conformational changes and an equilibrium of multimers including dimers and tetramers. Although the structural basis for allosteric activation by Phe is poorly understood, PAH was recently identified as a putative morpheein, and the current work builds on this hypothesis [14–15]. Unlike previous views of PAH multimerization equilibria, the morpheein model contains two different PAH tetramers separated by a large kinetic barrier that requires multimer dissociation, conformational change in the dissociated state, and reassociation to an alternate tetramer. Herein we 1) detail the rationale behind developing a morpheein model for PAH allostery, 2) present protein structure models for components of a PAH quaternary structure equilibrium, and 3) provide experimental evidence on purified full length rat and human PAH consistent with the proposed equilibrium. The current work sets the stage for our long-term goal, which is to develop allosteric pharmacological chaperones to expand the PKU patient population responsive to small molecule therapy.
A contemporary explanation for PAH dysfunction causing PKU identifies the disorder as one of PAH protein stability and/or folding [10, 12, 16–18]. For multimeric proteins, where the final step in the folding process is multimer assembly, perturbation of an equilibrium of functionally distinct alternative multimers provides both a mechanism for allostery as well as a mechanism for protein dysfunction [19–20]. Extensive literature describes mammalian PAH as an allosteric enzyme [21–22]. At low Phe levels the protein has basal activity which allows retention of sufficient Phe to support protein biosynthesis. At high Phe levels protein activation is accompanied by a large conformational change which dramatically changes various characteristics of the protein, including the fluorescence spectrum and the susceptibility to proteolytic digestion [23–25]. PAH allostery is associated with a quaternary structure equilibrium consisting of dimers and tetramers with some propensity to form larger aggregates [21, 23, 26–29]. Phe serves both as substrate and as allosteric activator; allosteric activation by Phe has been shown to draw an equilibrium of tetramers and dimers toward the tetramers [28, 30]. A structural basis for allosteric Phe binding causing stabilization of a specific, fully active conformation of the tetramer is proposed herein.
The new model for allosteric regulation of PAH, its relationship to PKU, and the promise of developing pharmacological chaperones that function allosterically follows our recently reviewed work with the prototype morpheein porphobilinogen synthase (PBGS) whose dysfunction causes the inborn error of metabolism known as ALAD porphyria [31–32]. Four important corollaries exist between the behavior of PBGS and the model proposed for PAH.
The goals of this study are to prepare and support a model for the quaternary structure equilibrium that governs PAH activity and to develop methods that can evaluate the position of the PAH quaternary structure equilibrium in a variety of physiologically relevant environments. The developed methods will allow future investigation of the effect of PKU-associated variants on the position of the PAH quaternary structure equilibrium. The developed protein structure models will assist in pursuit of allosteric pharmacologic chaperones that promote PAH to achieve and maintain its active multimeric assembly.
Both human PAH tetramer models are based on the human PAH structure 2PAH, which is a tetramer of a truncated form of PAH that contains only the catalytic and multimerization domains. A human regulatory domain homology model (residues 19–117) was prepared from the rat PAH crystal structure 1PHZ, which is of a truncated rat PAH containing only the regulatory and catalytic domains. Modeling was performed using our program MolIDE, which includes side chain rotamer optimization [49–50]. The low activity human PAH tetramer model was prepared using the program CHIMERA [51–52] to align the catalytic domains of the rat structure 1PHZ with the human structure 2PAH and the human regulatory domain model. Loop modeling (amino acids 113–121) was used to attach a regulatory domain to each of the human PAH subunits to yield our structure of 4mer*, using the program YASARA . Following loop modeling, loop optimization included energy minimization by steepest descent methods. A final optimization of side chain rotamers was done in the context of the full multimer using SCWRL4 . Preparation of the high activity human 4mer model started with the preparation of a dimer of the human PAH regulatory domain homology model, prepared by aligning two regulatory domains with the ACT domain dimer seen in crystal structure 1PSD. Following energy minimization and side chain optimization, the human PAH regulatory domain dimers were manually positioned to optimize the connection to the catalytic domains in the region of amino acid 118 of two catalytic domains with amino acid 117 of two halves of the regulatory domain homology model dimer. Loop modeling followed by side chain optimization and energy minimization completed the preparation of the human PAH 4mer model. Similar methods were used to prepare models of rat PAH 4mer* and 4mer.
The expression construct for a the MBP fusion protein of human PAH was a kind gift of Dr. A. Martinez and protein purification followed standard procedures on an amylose column. Expression constructs for full length rat and human PAH were a kind gift of Prof. P. Fitzpatrick. All gifted plamids were amplified and the sequence was confirmed to be as published. Expression procedures included small variations on the procedures used for expression of mammalian PBGS . Expression typically was carried out at reduced temperature (15 – 25 °C), Fe(II) was added at the same time as IPTG (at various concentrations), and in some cases glycerol was included as a chemical chaperone. None of these variations had a significant impact on protein expression levels. Purification of rat PAH followed the procedure of Shiman wherein crude cell extract was passed over a 75 ml Phenyl Sepharose-CL4B column in the presence of phenylalanine, followed by extensive washes with a variety of buffered solvents, and elution of the protein following removal of Phe from the elution buffer, which contains 15% glycerol, all carried out at 4 °C . This affinity purification is based on a major conformational change that has long been known to accompany the process wherein Phe activates mammalian PAH . Although this purification works well for rat PAH, human PAH, which expressed poorly, emerged from the column requiring significant further purification, which was carried out by ion exchange chromatography (IEC) on a 1 ml HiTrap Q column in 30 mM Tris-HCl, pH 7.4 using a KCl gradient (0–0.4 M) at room temperature.
Analytical SEC (100 μl sample size) was carried out on a 24 ml Superdex 200 column at 0.5 ml/min using 30 mM Tris pH 7.4, 150 mM KCl buffer with or without various concentrations of Phe. Analytical IEC was carried out on a 1 ml HiTrap Q column at a flow rate of 0.5 ml/min using 30 M Tris pH 7.4 with a KCl gradient from 20 mM to 0.4M with a variety of additives (see results).
Electrophoretic methods used the GE PhastSystem using 12.5% acrylamide gels and either native or SDS buffer strips. Western blotting on both native and denaturing gels used the WesternBreeze Chromogenic Western Blot Immunodetection Kit, WB7103 (Invitrogen). The primary antibody was MAB5278 (Millipore,10 ug/blot), which recognizes an epitope in the catalytic domain that is common to human tryptophan hydroxylase, tyrosine hydroxylase, and PAH.
Rat PAH was subjected to limited proteolysis by trypsin (Roche sequencing grade). Trypsin stock solutions were prepared in 1 mM HCl, and the total amount of HCl added to each reaction was normalized. Digestion reactions were carried out at room temperature in phosphate buffered saline pH 7.4 with PAH held constant at 0.2 mg/mL. Ratios of trypsin:PAH were varied from 1:5 to 1:500, in the presence or absence of 1 mM Phe. At various times, aliquots of the digestion were removed and the digestion was terminated with 10 mM PMSF prior to resolution of the digestion products by SDS-PAGE.
PAH activity determination used a continuous fluorescence assay monitoring the production of tyrosine  adapted to single microcuvette format for use with a Photon Technologies International fluorescence spectrophotometer equipped with a xenon flash lamp. Fluorescence was monitored at λEx = 275 nm and λEm = 305 nm, with excitation and emission slit widths set to 2 nm. The standard assay (200 μL final volume) contained 20 mM Na-Hepes pH 7.3, 10 μM ferrous ammonium sulfate, 2 mM dithiothreitol, 1 mg/mL catalase (to yield 2000–5000 U/mL), 1 mM Phe, 75 μM BH4, and 1–5 μg PAH at 30 °C. For experiments monitoring “order-of-addition” effects, the assay components were allowed to equilibrate at 30 °C for 5 minutes (rat PAH) or 15 minutes (human PAH) prior to starting the reaction by addition of the final assay reaction component. Quantification of tyrosine was performed by constructing standard curves of tyrosine under assay solution conditions in the presence of the relevant concentrations of BH4, which can quench tyrosine fluorescence .
Human and rat PAH were diluted to 0.2 mg/mL in 20 mM Na-Hepes pH 7.3, 2 mM DTT, 10 μM ferrous ammonium sulfate and 200 mM KCl in the presence or absence of 1 mM Phe. Emission spectra were recorded at 25 °C with λEx = 295 nm, and the excitation and emission slits set to 2 and 5 nm, respectively.
A protein structure model for the full length human PAH (Fig 1b, coordinates included in Supplementary Information) was prepared by combining 1) a homology model of the regulatory domain drawn from the crystal structure of a truncated rat PAH containing both the regulatory and catalytic domains (PDBid 1PHZ) with 2) the crystal structure of a truncated human PAH containing both the catalytic and multimerization domains (PDBid 2PAH). Similar straightforward combined models of full length tetrameric PAH have been presented previously : they provide the standard representation of tetrameric PAH . In these models the most N-terminal ordered region partially blocks access to the Fe-containing enzyme active site (Fig 1b). It is well documented that removal of the regulatory domain renders the protein constitutively active; this form is insensitive to Phe activation and has a specific activity comparable to the Phe-activated form of the full length protein . Thus the tetramer shown in Fig 1b is assigned to the full length protein assembly with low or basal activity, which we designate 4mer*. The conformational change that differentiates the low activity from the high activity multimers is posited to disallow this inhibitory interaction between the regulatory and catalytic domains. A similar model for rat PAH 4mer*, where the multimerization domain is homology modeled, is included in Supplementary Information.
Formulating a structural basis for PAH allostery draws on the comprehensive review of Kaufman , which describes extensive biochemical, biophysical, and kinetic data on full-length human and rat PAH isolated from mammalian tissue. Although not previously interpreted in this framework, published data are largely consistent with a dissociative allosteric model as detailed for the morpheein model of allostery [14–15, 19, 59–60]. The proposed dissociative mechanism for allosteric regulation of mammalian PAH consists of an equilibrium of two structurally and functionally distinct PAH tetramers (4mer*, which is low activity, and 4mer, which is high activity) where the spatial relationship between the regulatory and catalytic domains changes dramatically. The required conformational change is proposed to be sterically hindered in the tetrameric states, and to occur at the level of a dimer, hypothetically through the pathway illustrated in Fig 1b & c. The structure of the activated tetramer is conceived to include a protein-protein interaction, which is not present in the low activity tetramer, and which provides a hypothetical binding site for an allosteric Phe. Allosteric Phe binding is supported by published Phe binding stoichiometry studies [61–62]. Based on the behavior of ACT domains, described in more detail below, allosteric Phe binding is proposed to occur at an interaction interface between the regulatory domains of two subunits (Fig 1c, coordinates included in Supplementary Information). This proposition is consistent with a recent biophysical characterization of the isolated rat PAH regulatory domain, which was unexpectedly discovered to be both dimeric and capable of binding Phe . Consistent with our model, early work on the kinetic characteristics of dimeric and tetrameric mammalian PAH yielded the hypothesis that an allosteric Phe binding site is present in a tetramer and absent in the dimer . According to the proposed equilibrium, 4mer* 2mer* 2mer 4mer, allosteric Phe binding would stabilize the 4mer-specific interaction and draw the equilibrium toward 4mer. The active sites of 4mer, as modeled, have unfettered access to substrates (Fig 1c). Although the data presented herein do not address whether dissociation to the dimer is an essential component of the interconversion between the low activity tetramer (4mer*) and the high activity tetramer (4mer), we provide evidence for inclusion of a PAH dimer in the equilibrium as well as evidence for additional species in the equilibration of 4mer* with 4mer. A significant unknown in modeling the pathway between the tetramers is an uncertainty as to the conformation of the most N-terminal portion of the regulatory domain in the active 4mer (see below). Hydrogen/deuterium exchange studies show this region to be flexible both before and after Phe activation . Indirectly supporting a role for dimers in this pathway, the same studies show significant flexibility in the C-terminal helical portion of the multimerization domain both in the presence and absence of the allosteric Phe.
A large portion of the PAH regulatory domain contains a sequence motif known as an ACT domain, named for three allosterically regulated enzymes containing this sequence (aspartate kinase, chorismate mutase, and TyrA (a.k.a. prephenate dehydrogenase)) [33, 35]. Modeling the structure of 4mer draws on the structure(s) of ACT domains in the PDB  as analyzed by ProtCID and PDBfam [65–66] as well as a recent review on the role of ACT domains in the regulation of amino acid and nucleoside metabolism . Analysis of ACT domains in the PDB reveals that ACT domains commonly form dimers, as illustrated in Figure 2, where the ACT:ACT interface is often seen bound to a regulatory amino acid ligand. Publications on ACT domains often cite PAH as an outlier that does not include this interface [35–36]. However, applying this multimerization motif to the commonly accepted PAH tetramer structure, 4mer*, would require a ~90° rotation of the regulatory domain relative to the catalytic domain, not unlike the ~114° rotation seen in an ACT-like regulatory domain of 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase (DAH7PS) (see discussion) . In building the high activity or allosterically activated PAH tetramer model, the initial step was to dimerize the human PAH regulatory domain homology model by forming an interface between the ACT domain segments. The ACT domain dimer seen in the crystal structure of phosphoglycerate dehydrogenase (PDBid 1PSD, Fig 2b) was used as a template for orienting the PAH ACT domains. The resultant regulatory domain dimer(s) were attached to the human PAH tetramer (PDB code 2PAH) at the hinge region that connects the regulatory and catalytic domains to yield the 4mer model illustrated in Figure 1c. One aspect of the 4mer model which has low confidence is the conformation of the N-terminal 33 amino acids and its orientation relative to the ACT domain; the model simply retains these structural elements from the rat PAH crystal structure. The coordinates for the human PAH 4mer model and a similarly prepared rat PAH 4mer model are both included in Supplementary Information.
The models of 4mer* and 4mer are consistent with the well documented behaviors of PAH. For example, fluorescence ascribed to Trp120 and tryptic susceptibility in a highly basic region at the C-terminal region of the regulatory domain are two characteristics well documented to differentiate low activity PAH from Phe-activated PAH [23, 68–69]. Figure 3a illustrates that these regions of PAH are in different environments in 4mer relative to 4mer* consistent with published data. Using purified rat and human PAH, described below, we confirmed the 10 nm red shift in the PAH fluorescence spectra upon addition of sufficient Phe known to provide maximum activity stimulation (~1 mM Phe) . We also show that the susceptibility of purified rat PAH to limited tryptic digestion is a function of whether or not Phe is present (Fig 3b). Under the conditions shown, 1 mM Phe provides nearly total protection from limited tryptic digestion, which has previously been shown to cleave within the basic region at the C-terminal end of the regulatory domain . For rat PAH (Fig 3a) the surface accessibility of basic residues in this region decrease significantly when the conformation changes from the low activity to the high activity tetramers (e.g. for Arg111, Lys113, and Arg123, modeled surface accessibilities decrease to 69, 25, and 79 percent respectively). Hydrogen/deuterium exchange studies indicate that the entire interface between the regulatory and catalytic domains is more solvent exposed in the allosterically activated form, consistent with the 4mer (more exposed) and 4mer* (less exposed) models . The proposed models also provide a structural explanation for some published data on disease-associated PAH variants. For example, ex post facto inspection of 4mer reveals that Gly46, Ala47, Thr63, His64, Ile65, and Arg68 correspond to the amino acid ligand binding sites seen in ACT domain dimers (compare Figs 2b, c, d, and e), and would thus be predicted to be in the vicinity of the allosteric Phe binding site on 4mer. Notably, and in strong support of our model for ACT domain dimerization forming the allosteric Phe binding site, in vitro characterization of PAH fusion proteins for the PKU-associated variants G46S, A47T, T63P/H64N, I65T, and R68S showed that these mutations eliminate Phe binding .
Although the proposed allosteric and structural models are consistent with published data on mammalian PAH, a unique aspect of the morpheein model is an equilibrium of two tetrameric assemblies separated by an unusual kinetic barrier which is proposed to include multimer dissociation, conformational change in the dissociated state, and assembly to an alternate multimer structure. In particular, the novel 4mer model suggests that the alternate tetramers can be separated using ion exchange chromatography (IEC); Fig 1d shows the enhanced anionic character of the 4mer surface relative to the 4mer* surface. IEC has been useful in distinguishing between the PBGS hexamer and octamer and in evaluating how disease associated variants shift the PBGS quaternary structure equilibrium . We foresee IEC as a major tool in the future evaluation of disease-associated PAH variants. Using such classic biochemical techniques to establish an equilibrium of slowly exchanging protein quaternary structure isoforms can be straightforward if one presupposes the existence of such an equilibrium in the interpretation of experimental data. In addition to IEC, particularly useful and generally accessible experimental methods are native polyacrylamide gel electrophoresis (PAGE) and size exclusion chromatography (SEC). For an equilibrium of homo-multimers, native PAGE separates on the basis of multimerization stoichiometry because all multimers of a given protein have the same charge/mass ratio. SEC provides an independent estimate of multimerization stoichiometry (size). Both native PAGE and SEC can deviate from ideal behavior due to differences in shape, and such evidence is presented for PAH to support tetramers of significantly different shape (see below). Although these techniques are relatively low resolution, as described below, together they present compelling data in favor of a PAH quaternary structure equilibrium that consists of dimers and at least two architecturally different tetramers.
Experimental data has been obtained on mammalian PAH, heterologously expressed in Escherichia coli, from three different expression vectors. Because the human protein is reported to be unstable, most of the published work utilizes either a fusion protein or fragments lacking the regulatory domain. Below we describe studies of human PAH purified with a cleavable N-terminal maltose binding protein (MBP) as well as studies with full length human PAH expressed and purified without a tag. Fortunately, full length rat PAH is reported to be more tractable and expression, purification, and characterization of the tag-less rat protein informed our approach to working with tag-less full length human PAH. Our long-range goal of understanding how PKU-associated variants perturb a human PAH quaternary structure equilibrium demands establishing a baseline for the full length wild type human protein.
Initially human PAH was expressed as the MBP-PAH fusion and was found to largely follow published behavior , which is consistent with the proposed allosteric model. However, in the modeled fragile protein quaternary structure equilibrium, where the long-range plan is to evaluate disease-associated variants for their effect on said equilibrium, it is not reasonable to expect the MBP-PAH fusion to reflect a physiologically relevant equilibrium representative of full length PAH. For example, an N-terminal His-tag on DAH7PS has been shown to lock the position of its ACT-like regulatory domain and prevent its rotation . Various conditions used to cleave MBP from the fusion protein resulted in either incomplete digestion or undesired nonspecific proteolysis yielding PAH dimers and tetramers that were heteromeric in size, surface charge, and charge/mass ratio, and thus unacceptable for the proposed methodology.
Full length (tag-less) rat and human PAH were expressed from published pET constructs with superior expression observed for the rat protein . Following the Phenyl Sepharose affinity method of Shiman  the yield for rat PAH was ~12 mg at >90% purity (SDS PAGE in Fig 4a) from 9.1 g of E. coli cell paste with a specific activity of 7900 nmol min−1 mg−1 when addition of PAH is used to initiate activity. Expression of full length human PAH in E. coli used a variety of protocols, all of which provided relatively poor yield. Consequently 47.7 g of E. coli cell paste was used in the classic Shiman single column purification protocol. Unlike the behavior of rat PAH, human PAH activity eluted very slowly from the Phenyl Sepharose CL-4B (PS) column and the resulting PS pool measured 670 ml. This pool was concentrated to ~15 ml at ~ 1 mg/ml where SDS PAGE showed it to be of poor purity (Fig 4b) with a specific activity of only 1200 nmol min−1 mg−1. Nevertheless, an SDS Western blot yielded a single band for PAH (Fig 4b). A portion of the PS purified human PAH was further purified by IEC, described in more detail below, which yielded ~ 1mg of pure protein (Fig. 4b).
The PS purified rat PAH was used to illustrate dramatic differences in PAH catalyzed reaction rates during assays that vary only in the order-of-addition of reaction components; this kinetic behavior is one characteristic of morpheeins . As illustrated in Figure 5, when the PAH assay is initiated by the addition of BH4, the initial rate is rapid and decays as O2 becomes limiting. When the assay is initiated by the addition of Phe, the initial rate is quite slow. An intermediate rate is seen when both Phe and BH4 are used to initiate the assay. One interpretation of these anomalies is that Phe stabilizes 4mer, yielding high initial activity, while BH4 (in the absence of Phe) stabilizes a form with low or basal activity, perhaps 4mer*. IEC purified human PAH showed nearly identical order of addition effects (not shown).
Considerable information on the behavior of PAH was derived from chromatographic and electrophoretic data. Native PAGE (Fig 6a) and SEC (Fig 6b) are consistent with the freshly purified rat PAH protein (at neutral pH with 15% glycerol) being predominantly tetrameric with a small amount of dimer; the dimeric component is more apparent in the Western blot relative to the Coomassie stained gel. Close inspection of the native PAGE results suggests two different rat PAH tetramers of different architecture. Close inspection of the SEC results suggests that rat PAH defrosted from −80 °C at 37 °C may contain more than one tetrameric component. This is not seen when the protein is defrosted at 25 °C, which is the temperature later used to defrost all other PAH samples. For comparison, a Coomassie stained SDS gel of freshly purified rat PAH, shown in Fig 4a, indicates a single band at the expected 53 kDa. To address whether or not IEC could resolve components of the quaternary structure equilibrium, purified rat PAH (~100 μg) was subjected to analytical IEC in a neutral pH buffer with a KCl gradient (Fig 6c) in the presence and absence of Phe. In the absence of Phe, the chromatogram showed predominantly one peak centered at 256 mM KCl. These same buffer conditions used for SEC analysis showed predominantly tetramer with a small amount of dimer (Fig 6b). Inclusion of 1 mM Phe in the IEC elution buffer, which is sufficient to fully activate rat PAH, dramatically altered the elution profile such that the main peak is centered at 294 mM KCl. The observed increased retention of rat PAH upon exposure to 1 mM Phe is consistent with the surface charge predicted for 4mer relative to 4mer* in rat PAH tetramer models (see Fig 1d). SEC (Figure 6b) indicates that 1 mM Phe does not dramatically alter the size of tetrameric PAH, although the precise elution of the PAH tetramer is slightly earlier in the presence of 1 mM Phe, consistent with the notion that 4mer* and 4mer have somewhat different shapes. The key demonstration that rat PAH exists in an equilibrium of multiple quaternary structure isoforms is its behavior on the analytical IEC column when 50 μM Phe is included in the elution buffers, but not the PAH sample; 50 μM Phe is sub-saturating for the allosteric activation of PAH. Under these conditions, rat PAH separates into both different tetramers, as well as additional forms (Fig 6c). At 50 μM Phe, rat PAH runs as an equilibrium of the form eluting at ~260 mM KCl (presumably 4mer*), the form eluting at ~300 mM KCl (presumably 4mer), as well as additional components that elute at higher KCl (~380 mM). The appearance of additional distinct chromatographically separable components is consistent with the participation of dimeric forms in the equilibration between alternate PAH tetramer structures.
Using knowledge gained from the IEC behavior of rat PAH, the low-purity PS-purified human PAH was further purified by IEC, yielding two pools of activity (Fig 7a). The first peak, well resolved from the others, appeared pure by Coomassie stained SDS PAGE (Fig 4b) and of high specific activity. The fractions were brought to 15% glycerol and stored as the central fraction (Pool A, 6800 nmoles tyrosine min−1 mg−1), and the edge fractions (Pool A1, 6100 nmoles tyrosine min−1 mg−1). The region of disperse low activity (Pool B) coeluted with other proteins and was of lower specific activity (700 nmoles tyrosine min−1 mg−1). Coomassie stained native PAGE and native Western blots of these Pools are illustrated in Fig 7b, which shows that although Pool A is pure PAH (see Fig 4b), it contains tetramer, dimer, and some higher order PAH multimers that barely enter the gel (see Discussion).
In order to investigate whether the IEC-resolved peaks of human PAH are in equilibrium, both Pools A1 and B were dialyzed against the PS elution buffer, and individually re-chromatographed by IEC. In case of Pool A1, rechromatography resulted in several peaks, notably one at 85 minutes. In the case of Pool B, rechromatography resulted in the reemergence of the original Pool B peaks plus a small peak at 74 min. The combination of these analyses supports the conclusion that the form of human PAH that elutes at 74 min is in equilibrium with a form of human PAH that elutes at 85 min. SDS PAGE (not shown) indicates that the protein subunits do not change molecular mass during this equilibration.
IEC purified human PAH (Pool A) was used to evaluate the effect of Phe on IEC mobility (Fig 7c). In this experiment 100 μl aliquots of protein were diluted to 1 ml in the presence of Phe at 0, 50 μM, and 1 mM, followed by a 1 hour incubation prior to IEC; the results showed the main peak of human PAH migrating at 214 mM KCl in the absence of Phe and at 247 mM KCl in the presence of 1 mM Phe. At 50 μM Phe a broad peak of intermediate mobility is seen. Without the incubation, the chromatographic changes were attenuated. Although interpretation of human PAH IEC chromatograms is not as clear-cut as a similar study done with rat PAH (Fig 6c), the data indicate that the Phe induced conformational change causes the protein to be further retained on an anion exchange resin.
To demonstrate that Phe stabilizes a PAH tetramer relative to the dimer(s), the predominantly tetrameric rat PAH (~ 1 mg/ml, Fig 6a), was incubated at ~ 0.2 mg/ml in the presence and absence of ligands. Under these conditions, time points (not shown) indicate that the protein slowly equilibrates between tetramer and dimer, yielding about a 50/50 mixture after 24 h at 25 °C (Fig 8). Phe at 1 mM provided complete stabilization of a tetramer, presumably 4mer; Phe at 50 μM did not appreciably stabilize the tetramer; BH4 at 75 μM appears to promote formation of a dimer; and this effect is modulated in the presence of both 75 μM BH4 and 50 μM Phe.
We demonstrate for the first time that mammalian PAH exists in an equilibrium consisting of tetramers of two distinct architectures and that one or more dimeric forms of PAH are part of the equilibrium. The published behavior of mammalian PAH is consistent with a dissociative allosteric mechanism [21–22], which is approximated by a morpheein model of PAH allostery consisting of dimers and tetramers in the equilibrium 4mer* 2mer* 2mer 4mer; 4mer is the most active assembly and has a distinct allosteric binding site for Phe that is not present in the other components of the equilibrium. 4mer stabilization is proposed as the mechanism wherein Phe acts as an allosteric activator of mammalian PAH. In 4mer the N-terminal portion of the regulatory domain cannot sample the active-site-blocking conformation illustrated for 4mer* and seen in the rat PAH two-domain crystal structure. Although the current work does not disallow a direct equilibration between the tetramers, e.g. 4mer* 4mer, which would complete a thermodynamic box, we propose this interconversion to be sterically forbidden. We are developing the necessary tools to address this question directly as we have done for PBGS [59, 73]. A key to this demonstration is the isolation of a heterologously expressed human PAH variant whose quaternary structure equilibrium is dramatically shifted, as we found in the F12L variant of PBGS. Such variants are expected to emerge from the planned characterization of select disease associated PAH proteins.
In formulating further hypotheses about the PAH quaternary structure equilibrium, we look to the behavior of human PBGS, where a morpheein model for allostery is 6mer* 2mer* 2mer 8mer, and the high activity 8mer is a fragile assembly . The human PBGS protein assembly defaults to the low activity 6mer* in response to a variety of specific point mutations, to basic pH, and to the absence of active site ligands . In fact the basic arm of the bell-shaped pH activity profile is attributed to a shift in the equilibrium from predominantly high activity 8mer to predominantly low activity 6mer* . Specifically, all of the point mutations that encode variants associated with the inborn error of metabolism ALAD porphyria result in protein with an increased propensity to sample the hexameric assembly . A variety of small molecule ligands have been described as acting at the active site as well as allosterically to modulate the human PBGS quaternary structure equilibrium and thus modulate PBGS activity [43–44, 47–48]. A plethora of literature on mammalian PAH suggests comparable behavior wherein a variety of factors such as pH and ligand binding are reported to affect enzyme activity and protein stability. Here we consider that the reported inability of some disease associated PAH variants to achieve an active conformation may reflect an altered quaternary structure equilibrium rather than a failure to properly fold the individual subunits. Accumulation of the low activity components of the 4mer* 2mer* 2mer 4mer equilibrium is not an accumulation of improperly folded PAH; it is an accumulation of protein in a physiologically meaningful state with basal activity. Although we cannot predict the percentage of PKU-associated mutations that shift the PAH quaternary structure equilibrium, based on our experience with small molecule modulation of the PBGS quaternary structure equilibrium, we are confident that small molecules will be found that effectively shift the equilibrium toward the active 4mer.
It is well documented that PAH has a tendency to “aggregate”, or to form higher order multimers, and we see evidence for this in our native Western blots (e.g. Fig 7b). The formation of 4mer from two 2mers, stabilized by ACT domain dimerization, forms a closed system without the ability to further multimerize. This is illustrated schematically in Figure 9a. However, one can model the observed higher order multimerization through the interaction between ACT domains of the 4mer* and 2mer* components of the PAH quaternary structure equilibrium. The association of ACT domains between 4mer* and/or 2mer* components always leaves an exposed ACT domain interaction surface that can further multimerize (Fig 9b). It is unknown if these higher order multimers exist in vivo or if they have any meaningful place in PAH allostery. Other multimeric allosteric metabolic enzymes such as inosine monophosphate dehydrogenase and CTP synthase have been shown to form higher order multimers that appear to function meaningfully in aspects of cell cycle control .
DAH7PS (E.C. 188.8.131.52) catalyzes the first step in de novo aromatic amino acid biosynthesis via the shikimate pathway. DAH7PS, like PAH and PBGS, has evolved a phylogenetically diverse repertoire of mechanisms for allosteric regulation; for DAH7PS the regulation appears to be a classic feedback inhibition by tyrosine or phenylalanine. T. maritima DAH7PS is a homotetramer in which the structural basis for feedback inhibition was recently revealed to involve a large rotation (~114°) of an ACT-like regulatory domain relative to a catalytic domain . Similar to the proposed allosteric regulation of PAH, the DAH7PS domain reorientation serves to gate access to the enzyme active site; this is not unlike the PBGS case wherein allosteric regulation modulates the opening/closing of an active site lid. In all three cases, the fundamental difference between the low activity state and the high activity state lies in the modulation of active site access, not in any repositioning of active site residues. Strikingly similar to the proposed PAH regulatory model, allosteric regulation of T. maritima DAH7PS involves dimerization of the ACT-like domain and tyrosine is seen bound at the domain interface (Fig S1). Strikingly dissimilar is the dimerization mode (Fig S1). However, unlike PBGS and PAH, there is no evidence to date that allosteric regulation T. maritima DAH7PS involves tetramer dissociation.
The authors gratefully acknowledge support from the Developmental Therapeutics program at the Fox Chase Cancer Center, from the National Institutes of Health, National Cancer Institute Comprehensive Cancer Center grant P30CA006927, and the National Institute of General Medical Sciences R01GM84453 (RLD). We acknowledge contributions from the FCCC Molecular Modeling Facility, the Glasswashing Facility, and the DNA Sequencing Facility.
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.