Huan
Ma
a,
Klaudia
Szeler
b,
Shina C. L.
Kamerlin
*b and
Mikael
Widersten
*a
aDepartment of Chemistry – BMC, Uppsala University, Box 576, SE-751 23 Uppsala, Sweden. E-mail: mikael.widersten@kemi.uu.SE
bDepartment of Cell and Molecular Biology, Uppsala University, Box 596, SE-751 24, Uppsala, Sweden. E-mail: kamerlin@icm.uu.SE
First published on 17th November 2015
DERA, 2-deoxyribose-5-phosphate aldolase, catalyzes the retro-aldol cleavage of 2-deoxy-ribose-5-phosphate (dR5P) into glyceraldehyde-3-phosphate (G3P) and acetaldehyde in a branch of the pentose phosphate pathway. In addition to the physiological reaction, DERA also catalyzes the reverse addition reaction and, hence, is an interesting candidate for bio-catalysis of carbo-ligation reactions, which are central to synthetic chemistry. An obstacle to overcome for this enzyme to become a truly useful biocatalyst, however, is to relax the very strict dependency of this enzyme on phosphorylated substrates. We have studied herein the role of the non-canonical phosphate-binding site of this enzyme, consisting of Ser238 and Ser239, by site-directed and site-saturation mutagenesis, coupled to kinetic analysis of mutants. In addition, we have performed molecular dynamics simulations on the wild-type and four mutant enzymes, to analyse how mutations at this phosphate-binding site may affect the protein structure and dynamics. Further examination of the S239P mutant revealed that this variant increases the enthalpy change at the transition state, relative to the wild-type enzyme, but concomitant loss in entropy causes an overall relative loss in the TS free energy change. This entropy loss, as measured by the temperature dependence of catalysed rates, was mirrored in both a drastic loss in dynamics of the enzyme, which contributes to phosphate binding, as well as an overall loss in anti-correlated motions distributed over the entire protein. Our combined data suggests that the degree of anticorrelated motions within the DERA structure is coupled to catalytic efficiency in the DERA-catalyzed retro-aldol cleavage reaction, and can be manipulated for engineering purposes.
More than 30 different aldolases have been described to date, which are in turn classified into two major classes according to their catalytic mechanisms.3 Class I aldolases such as the enzyme of interest to this work activate the donor molecule by replacing the carbonyl oxygen with a nitrogen from an active-site Lys residue, thus forming an imine Schiff base as a reaction intermediate (Fig. 1). Note that the stereochemistry of the newly formed stereocentre is controlled by the enzyme regardless of the structure of the reactant substrate, which facilitates predictions of the product's stereo-configuration. These enzymes have attracted great attention as potential environmentally friendly alternatives for the catalysis of carboligation reactions, and much effort is being invested into understanding their structure–activity relationships, and their putative roles as greener catalysts in organic synthesis.4
Escherichia coli 2-deoxyribose-5-phosphate aldolase (DERA, E.C. 4.1.2.4) is the only known Class I acetaldehyde dependent aldolase, and one out of only two known aldolases that catalyze the addition of two aldehydes.5 DERA is a 28 kDa monomeric protein consisting of 259 amino acid residues. It is a key enzyme in the pentose phosphate pathway of E. coli, as it catalyzes the reversible cleavage of D-2-deoxyribose-5-phosphate (dR5P) to D-glyceraldehyde-3-phosphate (G3P) and acetaldehyde. DERA displays a typical TIM (α/β)8 barrel fold, which is a common structural feature among the members of the Class I aldolase family (Fig. 2A).7 In DERA, the nucleophilic ε-amine of Lys167 located on the β6 strand, attacks the aldehyde functionality of the incoming substrate (acetaldehyde or dR5P), forming the corresponding imine Schiff base. If the substrate is dR5P, the reaction will proceed as outlined in Fig. 1 with C–C bond breaking and ultimate release of the products G3P and acetaldehyde. The rate-limiting step of the aldol cleavage has not been pinpointed, but is assumed to involve C–C bond breaking, which is in turn facilitated by the rearrangement of the formed carbanion to the stabilized enamine.
Fig. 2 (A) (α/β)8 TIM barrel structure of DERA, highlighting the active-site Lys167 that forms the Schiff base with the incoming aldehyde (dR5P or acetaldehyde). (B) Close-up of the DERA active site, with the Lys167–dR5P complex shown in a stick representation. The side-chains of Ser238 (two conformers in the crystal structure) and Ser239 and Lys172 are shown in sticks. The “plus” indicates the location of the N-terminus of α-helix 8. Image created with Chimera6 using the atomic coordinates in PDB ID 1JCJ.7,8 |
DERA is strictly dependent on phosphorylated acceptor substrates, which is illustrated by the drastic loss of activity when comparing G3P and D-glyceraldehyde as acceptor substrates.5 The main phosphate binding residues are Lys172 (via a bridging water molecule), as well as direct side-chain interactions with Ser238, direct peptide backbone interactions with Ser238 and Gly205, and water bridge interactions with the backbone moieties of Gly171, Val206, Gly236 and Ser239 (Fig. 2B). The Ser238–Ser239 motif located close to the C-terminal is an atypical phosphate-binding motif, and is not conserved among other Class I aldolases.9 Finally, a putative conformational change has been proposed to take place in this region upon substrate binding.9
To investigate the importance of the phosphate binding site in DERA, we conducted mutagenesis of the Ser238/Ser239 pair by two approaches, firstly replacing these serine residues for prolines, and secondly introducing 11 other substitutions at either position using the NDT codon set. The resulting mutants were screened for catalysis of the reaction shown in Fig. 3. The rationale for the Pro replacements was to slightly move the position of the peptide backbone and thereby affect both side-chain interactions and backbone interactions with the phosphate group of the substrate. The shift in the position of the backbone was also intended to slightly move the N-terminal end of α-helix 8, thereby decreasing its proposed additional contribution from its positive dipole moment that can otherwise attract the phosphate group9 (Fig. 2B). The resulting S238P and S239P single mutants and the S238P/S239P double mutant were subsequently analysed for retro-aldolase activity towards dR5P, and we also examined the temperature dependence of the wild-type and S239P variants of DERA, as well as performing molecular dynamics simulations on all mutants.
Fig. 3 Reaction used to screen for S238X/S239X DERA variants with apparent activity towards an aryl-substituted aldehyde. Asterisks indicate 14C-label. |
Recent years have seen significant interest in understanding potential links between enzyme dynamics and correlated motions with both catalysis10–17 and also enzyme promiscuity and functional evolution.18,19 Following from this, there has also been discussion of the relevance of manipulating such dynamics in artificial enzyme design.20–22 Our combined experimental and computational analysis of the dynamical behaviour of these DERA mutants strongly suggests a role for coupled motions and entropic changes in driving the catalytic activity of this enzyme, and that the detrimental changes in the activities of the mutants we studied can be linked to both a loss in correlated motions as well as changes in activation entropies. Such dynamical changes can, in turn, be used to modulate the activity of this biocatalytically important enzyme.
Enzyme | k cat (s−1) | k cat/KMb (s−1 M−1) | K M (M) |
---|---|---|---|
a Reactions were assayed at pH 8 and at 30 °C. b Calculated from the total concentration of dR5P. | |||
Wild type | 15 ± 0.4 | (1.3 ± 0.07) × 105 | (9.6 ± 0.8) × 10−5 |
S238P | — | <0.1 | — |
S239P | 11 ± 1 | (3.6 ± 0.06) × 103 | (3.0 ± 0.6) × 10−3 |
S238P/S239P | — | <0.1 | — |
S238I/S239I | 1.4 ± 0.2 | (1.5 ± 0.06) × 102 | (9.6 ± 1.5) × 10−3 |
Interpreting the precise origin of this change is complex, as in a multi-step reaction such as the DERA catalysed aldol cleavage, KM is a composite constant that is built up of the rates for formation and decay of both the Michaelis complex of E·dR5P in both its furanose and open-chain forms, and also that of the formation and breakdown of the covalent Schiff base intermediate. Hence, destabilization of any of these enzyme-substrate species will increase the value of KdR5PM. Additionally, as the perturbed pKa of the ammonium group of Lys167 ensures the presence of the nucleophilic amine at neutral pH, it is essential for catalytic activity. Therefore, if the acidity of this group was to be impaired by the S239P mutation, formation of the aldimine would also be impaired, which would also contribute to an increase in the value of KdR5PM. The pH dependencies of kcat and kcat/KM do not indicate this to be the case, however; the pH-rate profile of the S239P mutant is unchanged from that of wild-type DERA within the pH region assayed (Fig. S1†).
There has been significant discussion about the role of entropy in enzyme catalysis (e.g.ref. 23–26). Therefore, the relative effects of the S239P mutation on activation enthalpies and entropies were also determined from the respective temperature dependencies of the catalysed rates (Fig. 4). As can be seen from Table 1, this mutation results in a ∼40-fold reduction of kcat/KM compared to wild-type DERA. The thermodynamic data suggest that the main negative effect on catalysis caused by the S239P mutation is primarily a relative loss of activation entropy (Table 2). That is, the relative effects on the enthalpies are favourable, and of similar magnitude, for both the E·dR5P/E-dR5P → E + product(s) reaction, which would be represented by kcat, and the E + dR5P → E + product(s) reaction, which is reflected in kcat/KM. The observation that the mutation results in relative decreases in activation enthalpies is counter-intuitive, as one might assume that the phosphate interactions in the wild-type enzyme would contribute also to TS stabilization. We note, however, that since the values presented in Fig. 4 and Table 1 represent a double difference, it is not clear whether the main effect is on the transition state or the ground state, and it could equally be a ground state as a transition state effect.
Fig. 4 Modified Eyring plots of the effect by the S239P mutation on catalytic rates. The ratios of kcat (circles) and kcat/KM (squares) are plotted as a function of temperature. The slopes are described by −ΔΔH‡/RT, and the intercepts by ΔΔS‡/R. The extracted values of the differences in enthalpic and entropic parameters are given in Table 2. |
ΔΔH‡a (kcal mol−1) | TΔΔS‡a,b (kcal mol−1) | ΔΔG‡c (kcal mol−1) | |
---|---|---|---|
a Standard errors from linear least square regression. b Calculated at 30 °C. c Calculated at 30 °C from −RTln(rateS239P/ratewt). The ΔΔ values represent the difference in each parameter between wild-type and mutant enzymes. | |||
k cat | −2.9 ± 0.5 | −3.1 ± 0.5 | 0.19 ± 0.05 |
k cat/KM | −4.1 ± 1 | −6.2± 1 | 2.2 ± 0.02 |
To further test the necessity of the presence of a bona fide phosphate binding site for the catalytic activity of DERA, as well as the possibility of adapting DERA into also accepting non-phosphorylated acceptor aldehydes, we performed limited saturation mutagenesis using the NDT codon set to replace Ser238 and Ser239. The resulting variants were assayed for their ability to incorporate 14C-labeled acetaldehyde into the acceptor phenyl acetaldehyde (Fig. 3). The assay is not quantitative but useful for screening purposes, and the reaction illustrated in Fig. 3 was specifically chosen based on both commercial availability of the arylsubstituted aldehyde, and the current industrial interest in producing biocatalysts that can catalyse aldol addition.27,28 Fig. S2† shows a typical result from a screen of aldol addition activities in bacterial lysates from DERA-variant expressing cells. Sequence analysis of the scored hits revealed that Ser was no longer over-represented in the hits with apparent aldol addition activity. Rather, Ile displayed the highest over-representation at position 238, followed by Ser in second and Val in third places (Fig. S3†). The order of preferred residues at position 239 was Leu > Ile > His. Ser was absent at position 239 in all picked hit variants. Since these activity screens had targeted the aldol addition of acetaldehyde and phenyl acetaldehyde it was not obvious whether these DERA variants would retain activity with phosphorylated substrates such as dR5P. To test this, an Ile238/Ile239 variant was purified and analyzed for activity with dR5P. This mutant did indeed exhibit reasonable activity also with this phosphorylated substrate (Table 1), which clearly demonstrates that the Ser238/Ser239 pair is not critical for catalytic activity in and of itself.
Finally, in order to further explore the effect of mutations at positions 238 and 239 on overall structural stability, we performed molecular dynamics (MD) simulations of both wild-type DERA, and the S238P, S239P, S238P/S239P and S238I/S239I mutants in their substrate free forms. The simulation setup was as described in the Experimental section, and each variant was equilibrated for a total of 300 ns over 3 replicas with different initial velocities (100 ns each), leading to a total of 1.5 μs of cumulative simulation time for all variants. Fig. S4† shows the root mean square deviations of the backbone atoms for each variant, averaged over 3 replicas, demonstrating the convergence of the simulations. Following from this, we calculated root mean square fluctuations (RMSF) along the MD trajectories for each individual variant (Fig. 5), as well as calculating dynamic cross correlation maps (DCCM) from the MD trajectories, using the Bio3D package.29,30
From the RMSF plot, it can be seen that the enzyme is fairly flexible, with each peak on the RMSF plots corresponding to a different helix in DERA (see Fig. 5 for a colour coded comparison of the different helices of DERA and each peak on the RMSF plot), thus corresponding to “breathing motions” of this enzyme. As can be seen from these plots, the introduction of all mutants dampens the flexibility of the enzyme in several regions, and, in particular, the S239P mutant shows less flexibility around residues 25–75, 170–180 and 220–230 (see Fig. 5 for colouring). Following from this, a side-by-side comparison of the DCCM plots for the wild-type and S239P type mutant forms of DERA (Fig. 6) show that this substitution leads to a substantial loss in, in particular, anti-correlated motions. Interestingly, if one then extends this to placing all variants shown in Table 1 side by side, and ranking them according to their measured kcat/KM values, it can be seen that this loss of anti-correlated motions follows a clear trend across the mutants, with the variants with the lowest kcat/KM values (i.e. S238P and the S238P/S239P and S238I/S239I double mutants) showing much lower correlated motions (Fig. S5†). This is in good agreement with previous work that suggested the presence of coupled networks of residues in enzymes that play a role in enzyme function (e.g.ref. 10–17 among many others), and in particular work by Bruice which argues for anti-correlated motions as a driving force in catalysis.31 Therefore, the introduction of the mutations at positions 238 and 239 appears to alter the enzyme's dynamics, rigidifying it, which will in turn lead it to sample fewer catalytically active conformations and impair both substrate binding and subsequent TS stabilization, as shown also in the trends in kcat/KM in Table 1.
Fig. 6 Dynamic cross correlation maps (DCCM) of (A) wild-type DERA and (B) the S239P mutant. Corresponding DCCM plots for all mutants can be found in Fig. S5.† This figure was generated using Bio3D.29,30 |
DERA's requirement for phosphorylated substrates can be due to different reasons. (1) The obvious, and most trivial reason is to function as a negatively charged anchor binding to a complementary binding site, thereby increasing the affinity for the substrate in the ground state and throughout the reaction. This could be argued as being particularly important for small substrate molecules such as G3P. Any contribution to substrate selectivity, however, is questionable since the relatively strong binding interactions between the protein and phosphate moiety may overshadow the weaker (e.g. hydrogen bonding) interactions that steers substrate selectivity, including stereoselectivity. (2) Extensive experimental work by Richard and co-workers on other enzymes acting on phosphorylated substrates, which include orotidine monophosphate decarboxylase32–34 and triosephosphate isomerase,35–37 shows that remote protein–phosphodianion interactions can, in these cases, be crucial for transition state stabilization, in part by promoting the formation of (and stabilizing of) active high-energy conformations of the protein. (3) In the current case of DERA catalysed aldol cleavage of dR5P, the additional ring opening of the predominantly existing furanose form (99.8%38) of the aldose substrate to the reactive open-chain aldehyde has to be considered (the concentration of the aldehyde in solution is expected to be even lower than 0.2% of the total dR5P, since it will be further equilibrated with the hydrated gem-diol). The phosphate group on the furanose may in this case act as an affinity tag facilitating imperfect binding to the active site in which subsequent ring opening is catalysed. The ring opening can be facilitated both by acid catalysis but also by allowing for new binding interactions once the sugar is in the open-chain form, and the aldehyde thus formed is trapped as the imine intermediate after attack by Lys167.
For comparison, note that phosphoglucose isomerase, which catalyses the isomerization of glucose-1-phosphate and fructose-1-phosphate faces the same problem. This enzyme acts on the open chain of the aldose or ketose substrates, while these sugars are present in the respective pyranose and furanose forms in aqueous solution, to an even higher degree than in dR5P. Crystal structures have been solved in which both the cyclized and the open-chain form of the substrates can be observed within the same active site. In this enzyme, a His residue has been suggested to be the acid catalyst for sugar ring opening once bound within the active site.39 There is no residue in DERA that corresponds to the catalytic His residue present in phosphoglucose isomerase, but rather, there is a strategically positioned water molecule aligned in a hydrogen-bonding network involving the O4 of the open-chain Schiff base and the ammonium of Lys172 (PDB ID: 1JCJ7,8), suggesting possible acid functionality. The uncatalyzed ring-opening of furanoses is comparably slow, with rates well below 1 s−1.40,41 This is at least 15-fold below the value of kcat for cleavage of dR5P by DERA (Table 1). Hence, catalysed ring-opening must exceed at least the value of the turnover number. Phosphorylation, however, considerably increases the ring-opening rate (40–100 s−1for R5P, depending on anomer).42 The ring-opening of dR5P is therefore not expected to be rate-limiting to any extent. The aldehyde thus formed is trapped as the imine intermediate after attack by Lys167.
Purification of DERA variants was accomplished through nickel IMAC, as previously described.43 The purified protein was stored in 50 mM triethanolamine buffer, pH 8.0. The rates of catalysed aldol cleavage of dR5P were assayed as described in ref. 45, and steady state kinetic parameters were determined after fitting the Michaelis–Menten equation to the data using SIMFIT (http://www.simfit.org.uk). In addition, catalyzed rates, under saturating (kcat) or unsaturating (kcat/KM) conditions, in the presence of wild-type DERA or the S239P mutant, were determined at different temperatures (15 °C, 20 °C, 25 °C, 30 °C, 38 °C and 45 °C). The natural logarithms of the ratios of the determined rates were plotted as functions of the inverse of T to allow for extraction of the changes in transition state enthalpy and entropy (eqn (1)), caused by the S239P mutation.
(1) |
The steady state kinetic parameters for wild-type DERA, and the S239P and S238I/S239I mutants were also determined at pH values of 7, 8 and 9. The results are presented in Fig. S1 and Table S2.†
Footnote |
† Electronic supplementary information (ESI) available: Further experimental and computational data. See DOI: 10.1039/c5sc03666f |
This journal is © The Royal Society of Chemistry 2016 |