Louis Y. P.
Luk
,
E. Joel
Loveridge
and
Rudolf K.
Allemann
*
School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, UK. E-mail: allemannrk@cardiff.ac.uk
First published on 9th April 2015
The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.
Dihydrofolate reductase (DHFR) has been widely used in studies of the relationship between enzyme structure, motions and catalysis.8–30 DHFR transfers a hydride from NADPH and a proton from the solvent to C-6 and N-5 of 7,8-dihydrofolate (DHF) to generate 5,6,7,8-tetrahydrofolate (THF) (Fig. 1). DHFR from Escherichia coli (EcDHFR) contains a number of mobile segments including the M20 (residues 9–24), FG (residues 116–132) and GH (residues 142–149) loops and switches between a closed and an occluded conformation during the catalytic cycle (Fig. 2).16,31 Upon binding of substrate and cofactor, the M20 loop closes over the active site by forming stable hydrogen bonds with the FG loop16 to create an optimal electrostatic environment conducive to hydride transfer.13,16 Once the products have formed, the M20 loop releases the nicotinamide ring of the oxidized cofactor and occludes part of the active site by forming an alternative hydrogen bonding pattern with the GH loop,16 which triggers the exchange of NADP+ and NADPH. Finally, the product THF is released from the active site in the rate-limiting step at pH 7 and the enzyme returns to the closed conformation and another reduced cofactor enters the active site.16
![]() | ||
Fig. 2 Cartoon representations of EcDHFR (1RX2),16 EcDHFR-N23PP/S148A (3QL0),24 MpDHFR (3IA5),54 TmDHFR (1D1G)39 and BsDHFR (1ZDR).35 NADP+ and folic acid serve as ligands in EcDHFR; the M20 (red), FG (yellow) and GH (orange) loops are highlighted. |
A number of DHFRs have been shown to follow similar catalytic cycles to EcDHFR, including DHFR from humans32 as well as those from bacteria including Lactobacillus casei (LcDHFR)33 and the psychrophile Moritella profunda (MpDHFR).34 In all these cases, and for DHFR from the thermophilic bacterium Geobacillus stearothermophilus (BsDHFR),35,36 the steady state turnover at pH 7 is limited by a physical step rather than the actual chemical step of hydride transfer. In contrast, catalytic turnover under steady-state conditions in DHFR from the hyperthermophilic bacterium Thermotoga maritima (TmDHFR) is at least partially limited by hydride transfer;37,38 TmDHFR has dramatically lower hydride transfer rate constants than other characterised chromosomal DHFRs.38 Unlike these DHFRs, TmDHFR relies on its dimeric structure for activity and with a melting temperature of 81 °C is the most thermostable DHFR characterised.39–41
The kinetic isotope effects (KIE) on the DHFR-catalysed reaction have been measured chiefly by two experimental techniques. For many DHFRs, the turnover number kcat at pH 7 reports mainly on product release14,32–34 and transient kinetic techniques such as stopped flow must be used to extract information about the chemical step of the catalytic cycle. These have the advantage that they can be performed under conditions where hydride transfer makes little or no contribution to kcat.42 We have made widespread use of single turnover stopped flow kinetics at pH 7 to study the DHFR-catalysed reaction. However, the observed single-turnover (or burst phase) stopped flow rate constants do not report exclusively on hydride transfer, but are complicated by typically small contributions from physical events.42 On the other hand, hydride transfer becomes increasingly rate limiting when the pH is elevated, so kcat at pH > 8.5 reports predominantly on hydride transfer with small contributions from physical steps.14 Competitive multiple-KIE measurements at elevated pH allow extraction of intrinsic KIEs, relatively free from kinetic complexity and reporting solely on the hydride transfer step, but such measurements are not ideal because the protonation state and conformational preference of the enzyme are altered, and the catalytic behaviour of the enzyme may not mirror that under physiological conditions42,43 and most computational investigations of the EcDHFR catalysed reaction have employed X-ray structures determined under conditions of neutral pH.18,21,44–46 The intrinsic KIE of EcDHFR has been estimated at various pH values showing that the degree of kinetic complexity depends on a number of factors, including reaction temperature and enzyme conformational behaviour as well as the protonation state of the enzyme.42 It has also been shown that the intrinsic KIE of EcDHFR at pH 7 is different to that at pH 9.47 Accordingly, single turnover KIEs under physiological and intrinsic KIEs at elevated pH must both be interpreted with caution.
![]() | ||
Fig. 3 Cartoon representation of EcDHFR (1RX2).24 Bound ligands are shown as sticks; the α-carbons of the residues discussed in the text are marked with orange spheres. |
Hydride transfer is simply not possible in the occluded conformation as the reactants are not sufficiently close to one another. More generally, the ability to form the occluded conformation does not affect the chemical step of the catalytic cycle directly, as shown both by the similarity of the single turnover rate constants of EcDHFR and MpDHFR at pH 712,48 and by the existence of EcDHFR variants such as EcDHFR-S148A that are incapable of forming an occluded conformation yet maintain wild-type-like single turnover rate constants.43,44 It is likely that only DHFRs with a Ser148 equivalent are capable of adopting an occluded conformation.43,44 However, even when large conformational changes are not present, protein motions can still be involved in progression through the catalytic cycle, as demonstrated by the changes in millisecond motions in different EcDHFR product complexes,31 all of which adopt the occluded conformation.16 In agreement with the observation that only certain DHFRs are likely to form an occluded conformation, and that this does not affect the overall turnover rates, it has also been shown that native-like protein equilibrium motions are not required for the function of beta-lactamases.55
As discussed elsewhere,56 such equilibrium motions are the consequence of the free energy surface and modulation of this surface by ligand exchange is well known. Enzyme motions clearly affect the electrostatic environment by changing the relative positions of the interacting atoms. However, electrostatic effects determine the form of the free energy surface, and motions of the enzyme simply reflect the change in interactions between the atoms as the enzyme moves across it; they do not actively alter the free energy surface. Therefore a change in the motions of the enzyme is always a consequence of a change to the electrostatics of the enzyme, not the other way round. This does not diminish the importance of studies of conformational and other equilibrium motions, as they greatly aid the understanding of the free energy landscape. The importance of the free energy landscape in tuning conformational changes across multiple steps of an enzyme's catalytic cycle is well known.57 Enzyme motions, or rather enzyme conformational flexibility, have also been suggested to be important for the evolution of new enzymatic functions.58
Binding of antifolates to EcDHFR dampens millisecond conformational motions,59 consistent with the EcDHFR catalytic cycle involving millisecond motions that gate ligand binding and release.31 Enzyme-inhibitor interactions will be more favourable than those obtained in alternative conformational states (placing the enzyme in what has been referred to as a ‘dynamic straightjacket’59), imposing changes to the free energy surface that cause these motions to be dampened. However, measuring the effect on the enzyme motions will potentially provide a quantitative measure of how strong an interaction needs to be to overcome conformational switching.
While bacterial DHFRs typically have a mobile M20 loop, mammalian ones are more rigid due to the presence of a proline-rich sequence at the end of that loop.24,28,29 Incorporation of a diproline repeat into the M20 loop of EcDHFR to form EcDHFR-N23PP (Fig. 3) prevented formation of the occluded conformation and caused a general loss of the millisecond motions found in the wild type Michaelis complex; it has therefore been proposed that such motions could be linked to the chemical step.24 A subsequent computational study suggested that the impaired catalytic activity results from changes to the reorganisation free energy of the reaction brought about by changes to the electrostatic preorganisation within the active site rather than directly from changes to the motional behaviour of the enzyme.8 An experimental investigation supported the dominant effect of changes to the activation entropy on the single turnover rate constant at pH 7.13 Although the temperature dependence of the intrinsic KIE differed considerably between the wild type and variant enzymes, which was interpreted as being due to changes in the values and distribution of the donor–acceptor distance,47 such changes are an equilibrium property of the enzyme, affected by changes to the active site electrostatics.
It was later shown that the extension of the M20 loop rather than the presence of proline per se was responsible for the loss of conformational flexibility,28 presumably by disrupting the ability of the M20 loop to form stabilising hydrogen bonds to the GH loop. Additionally, it was shown that an insertion in the folate-binding region, acquired earlier in evolutionary history than the proline-rich region, is necessary for high catalytic activity in mammalian DHFRs.29 EcDHFR-N23PP/G51PEKN (Fig. 3) has a comparable rate constant for hydride transfer to wild type human DHFR,29 and the temperature dependence of the intrinsic KIE at pH 9 is similar to that of EcDHFR.47 Hence, despite the similarities in protein architectures, and even where the kinetics of the chemical step are highly similar, the conformational landscapes of different DHFRs may be quite different.
The G121V mutation causes a change to the ground-state structure of EcDHFR including alterations to the active site itself.65,69–71 Evidence from NMR measurements indicated that EcDHFR-G121V adopts an occluded conformation in the ground state of the Michaelis complex72 and that conformational fluctuations sample a state other than the closed conformation.71 This is unsurprising as the isopropyl side chain introduced by the conversion of Gly to Val is likely to point toward the interior of the enzyme16 leading to unfavourable steric interactions with the M20 loop and destabilisation of the hydrogen bonds between the FG and M20 loops that are necessary for stabilisation of the catalytically competent closed conformation,65,73 as well as consequent disruption of correlated motions between these loops. Networks of coupled motions have been identified in TmDHFR74 and in DHFR from Bacillus subtilis,75 although differences between these and the network in EcDHFR were seen, likely reflecting differences in their conformational landscapes. Indeed, the G123V mutation in MpDHFR has little effect on the single turnover rate constant or its observed KIE at pH 7, consistent with the proposal that interactions between the FG and M20 loops do not play a significant role in MpDHFR catalysis.49
Binding of NADPH and the inhibitor methotrexate to EcDHFR-G121V forms a putative mimic of the transition state and causes the enzyme to adopt a closed conformation,26 confirming that EcDHFR-G121V is capable of forming this conformation for hydride transfer. However, the mutation leads to aberrant millisecond conformational switching of the M20 and FG loops,26 as would be expected when the closed conformation is strongly destabilised, suggesting that these motions are anti-catalytic and destabilise the optimum active site configuration.26 The EcDHFR-M42W complex with NADPH and methotrexate also forms a closed conformation and shows slower millisecond motion than the wild-type enzyme.76 Interestingly, the effects of both the G121V and M42W mutations on the ps–ns motions are small,26,76 although side chain methyl fluctuations are sensitive to the M42W mutation with many showing increased rigidity and a smaller number showing greater flexibility.76 Computational studies indicate that EcDHFR-G121V/M42W has reduced M20 loop conformational motions and reduced flexibility at the transition state, related to an increased entropic barrier. This is in contrast to EcDHFR-N23PP/S148A, where increased flexibility at the transition state77 and in the Michaelis complex73 was observed despite the loss of millisecond conformational motions.24
These results indicate a role for networks of coupled motions in ‘promoting’ the chemical step in DHFR. However, they do so through thermal motions that are directly affected by changes to the active site electrostatics brought about by the mutations, however minor. This again is an equilibrium property and a manifestation of changes to the free energy surface of the enzyme. As such these motions are not the same as ‘promoting motions’ or ‘promoting vibrations’ that have been proposed to modulate the reaction coordinate by actively reducing the barrier height and/or width.3,78–87 Indeed, computational studies of EcDHFR indicate that motions are mostly thermally dissipated before the chemical transformation takes place,8,27,45,68,77,88 and the influence of mutations at Gly121 in EcDHFR has been explained in terms of the effect on the free energy surface.21,89,90 Additionally, while statistical motions may be coupled over large distances, true dynamic correlations (i.e. involving non-statistical motions) only exist over very small distances and non-local, long-range dynamics therefore appear not to have a significant role in EcDHFR catalysis.10 Formation of an optimal ‘reaction-ready’ active site configuration could be interpreted either as an enzyme rearrangement following substrate binding or as part of the ascent on the free energy surface towards the transition state. In either case coupled motions and their proposed effect on hydride transfer are a manifestation of movement on the enzyme's energy surface rather than an effect in their own right; these motions do not ‘couple’ directly to hydride transfer.
Similar effects of active site volume have been proposed for enzymes such as morphinone reductase,92 soybean lipoxygenase93,94 and alcohol dehydrogenase95 and the effect of pressure on hydride transfer kinetics has been investigated and related to ‘promoting vibrations’ in the active site.96–98 Changes in active site volume and in pressure do however not cause true ‘dynamic’ effects, but alter the equilibrium conformational ensemble of the enzyme. It has been shown that pressure does not directly compress the free energy surface,99 and more generally that barrier compression would lead to reduced tunnelling and reduced KIEs,21,99 in contrast to recent claims that shorter donor–acceptor distances lead to reduced KIEs due to increased deuterium tunnelling.47
For TmDHFR, EcDHFR, MpDHFR and BsDHFR, viscosity had no effect on the single turnover rate constants at pH 7 and solvent composition did not affect the KIE.11,12,22,36 In contrast, the dielectric constant of the medium had a pronounced effect on both the single turnover and steady state rate constants at pH 7. Although no viscosity effect on the steady state rate constant was observed for TmDHFR,22 MpDHFR12 or BsDHFR,36 solvents with similar dielectric constants reveal a small, but consistent, effect from viscosity in the steady state for EcDHFR.11 This is consistent with the rate-determining step in the catalytic cycle of EcDHFR, but not those of the other three DHFRs, involving a substantial conformational change (vide supra).
We initially interpreted our results for TmDHFR in terms of an environmentally coupled tunnelling model of catalysis22 but it became clear that this model could not adequately explain the effect of co-solvents on catalysis by other DHFRs.11,12 The dominant effect of the dielectric constant provides strong support for an overriding role for electrostatics in controlling catalysis, as proposed for enzymes in general104,105 and DHFR specifically.21,53,90 Changes to the solvent composition will affect the free energy surface for the protein and so affect the reorganisation energy of the catalysed reaction. In contrast, the lack of a viscosity effect provides strong evidence against long-range motions on any timescale being involved in the chemical step.
The effect of protein isotope labelling is assessed by comparing the reaction rate constants for the light enzyme, kLE, and heavy enzyme, kHE, to give an ‘enzyme KIE’, kLE/kHE, where the reactants are not isotopically labelled. In contrast to primary hydrogen KIEs, where extensive literature exists for the interpretation of their temperature dependence,13,21,44,93,112–114 the effect of temperature on the enzyme KIE has not yet been fully analysed. To investigate dynamic (i.e. non-statistical) contributions to the hydride transfer rate constant ensemble-averaged variational transition-state theory (EA-VTST) calculations were performed that incorporate a correction to the quasi-classical rate constant, the transmission coefficient, which takes into account dynamic and tunnelling contributions to the reaction:115–117
Γ(T) = γ(T)·κ(T) |
To date, we have performed ‘heavy’ enzyme studies on four DHFRs, namely EcDHFR, EcDHFR-N23PP/S148A, the thermophilic BsDHFR and the hyperthermophilic TmDHFR. Under pH 7 steady state conditions, EcDHFR gave a small, temperature-dependent enzyme KIEcat (kLEcat/kHEcat) (Fig. 4),27 while for EcDHFR-N23PP/S148A kcat was not affected by isotope labelling.77 As the rate-limiting product release step in EcDHFR catalysis involves movement of the M20 and FG loops,31 whereas in EcDHFR-N23PP/S148A NADP+ release is rate limiting and is not coupled to a significant conformational change,24 a temperature-dependent enzyme KIEcat is likely an indication of the involvement of a conformational change in the rate-limiting step of the catalytic cycle.27,77 ‘Heavy’ BsDHFR gave a relatively large, temperature independent enzyme KIEcat of ∼2.6 (Fig. 4). As the temperature independence of KIEcat suggests no involvement of a large-scale conformational change in kcat, the measured kinetic difference was attributed to the inherent flexibility of the enzyme.111
![]() | ||
Fig. 4 Enzyme KIEs at pH 7 for EcDHFR (![]() ![]() ![]() ![]() |
At pH 7 the enzyme KIEH (kLEH/kHEH, where kH is the single turnover rate constant that reports predominantly on hydride transfer) for EcDHFR was slightly inverse at 5 °C but increased with temperature, passing through zero and increasing to give a small, normal value at 40 °C (Fig. 4).27 It was demonstrated that the data for EcDHFR could be fit to a single-conformer kinetic model, without the need to invoke multiple catalytically competent conformations.27,119,120 A weakly temperature dependent enzyme KIE was observed for hydride transfer catalysed by EcDHFR-N23PP/S148A, but with greater magnitude than for the wild type enzyme.77
EA-VTST calculations indicated that the tunnelling coefficient remains unchanged on isotope substitution and is not significantly different between the wild type and variant enzymes,27,77 in agreement with a previous computational study that found no role for ‘promoting motions’, including those in the network of coupled motions (vide supra), in driving barrier passage in EcDHFR.121 Instead, the observed difference in the single turnover rate constant is due to a small but significant difference between the recrossing coefficients in the ‘light’ and ‘heavy’ enzymes. It appears that the reaction in the ‘heavy’ enzymes is more prone to dynamic recrossing because the response from the environmental motions along the reaction coordinate is slowed.27,77 The non-productive recrossing trajectories were enhanced in the N23PP/S148A variant,77 suggesting that additional protein vibrations are incorporated into the transition state. Hence, although EcDHFR-N23PP/S148A is a ‘dynamic knock-out’ on the millisecond timescale, it is actually a ‘dynamic knock-in’ on the timescale of the chemical step, which contributes to the reduction in the efficiency of hydride transfer.77 While similar concepts have been discussed in previous computational studies,7,21 this report provided the first experimental evidence for the existence of ‘demoting vibrations’.77
The enzyme KIEH for the BsDHFR catalysed reaction showed biphasic behaviour, being relatively low (∼1.1) and temperature independent above 20 °C, but increasing sharply with decreasing temperature (Fig. 4).111 As observed previously,27,77 only the recrossing coefficients, not the tunnelling coefficients, were affected by isotope labeling,111 again in disagreement with the proposals of barrier modulation and the involvement of ‘promoting’ motions. The singly deuterium-labelled and doubly 15N and 13C-labelled BsDHFRs were also prepared to investigate the biphasic temperature dependence. Since the molecular weights of these enzymes increased by ∼5%, instead of ∼10% as observed for the triply 15N, 13C, 2H-labelled, fully ‘heavy’ BsDHFR, the enzyme KIEs measured with these partially labelled enzymes would reveal clues about the intensity of the dynamic coupling.111 Above 25 °C, the enzyme KIEH measured with these labelled enzymes were statistically the same as those with the fully ‘heavy’ enzyme, suggesting that dynamic coupling is limited.111 This also confirmed that the observed enzyme KIE is caused by alteration of the protein's dynamics rather than other mass-induced effects such as a change in van der Waals radii or hydrophobicity due to deuterium labelling.111 In contrast, at 5 °C the enzyme KIEH measured with the fully ‘heavy’ enzyme was noticeably higher than those of the singly and doubly labelled enzymes. The thermal melting temperatures of the ‘light’ and ‘heavy’ BsHDFRs also showed no significant difference.
In contrast to EcDHFR, EcDHFR-N23PP/S148A and BsDHFR, protein isotope labelling exerted no effect on the single-turnover rate constants for TmDHFR at pH 7 resulting in an enzyme KIEH of unity for all temperatures investigated (Fig. 4). Accordingly, mass-dependent protein dynamics do not appear to exert any detectable effect on the chemical step of TmDHFR catalysis. It may therefore be possible that TmDHFR has gained a small catalytic benefit by eliminating dynamic recrossing trajectories from the transition state dividing surface.110 It had previously been postulated that active site motions may be detrimental to catalysis and that this observation becomes more pronounced in enzymes from thermophilic species.122 Our studies however indicate that TmDHFR is able to overcome this, by eliminating dynamic coupling at high temperatures. Under steady-state conditions at pH 7 the enzyme KIEcat for TmDHFR remained at ∼1.35 from 15 °C to 65 °C, but increased sharply below 15 °C (Fig. 4). The temperature-independent region likely relates to the absence of significant conformational changes in the TmDHFR catalytic cycle, whereas the temperature-dependent region suggests a change in the conformational equilibrium favourable for reaction. It should be noted that wild type TmDHFR revealed subtle biphasic behaviour under pH 7 steady state conditions, but TmDHFR variants with disrupted dimeric interface did not show this behaviour.112,123 Hence, this hypothetical switch in conformational equilibrium could result from a change in inter-subunit interactions.110
Of further interest is the fact that KIEH increases with increasing temperature for EcDHFR27 (and to a lesser extent EcDHFR-N23PP/S148A77), is temperature independent (and negligible) for TmDHFR,110 but decreases with increasing temperature for BsDHFR.111 For BsDHFR, nine pairs of recrossing coefficients were calculated to investigate the temperature dependence of KIEH. The recrossing coefficient was shown to absorb all non-statistical effects; in particular, at low temperature it is noticeably lower in the ‘heavy’ enzyme than in the ‘light’ enzyme.111In silico studies also indicated that BsDHFR is more flexible on the nanosecond timescale,111 supporting earlier experimental studies.124–126 The reverse temperature dependence of the enzyme KIEH in BsDHFR relative to EcDHFR may be related to this intrinsic flexibility of BsDHFR.111 Heavy isotope labelling raises the enthalpic barrier for electrostatic preorganisation and the reorganisation energy for reorienting the substrates; consequently, the transition state experiences additional friction in the enzyme active site. For BsDHFR, this is more pronounced at low temperatures, because the enzyme lacks the energy required to adopt an ideal configuration, leading to an increased recrossing contribution in the ‘heavy’ enzyme and a strong enzyme KIEH. However, most of this enthalpic barrier can be overcome at physiological temperatures, allowing hydride transfer to proceed in an ideal electrostatic configuration with minimal dynamic effects. This was also evident from the activation parameters; the higher activation enthalpy (ΔH‡) in ‘heavy’ BsDHFR is accompanied by a smaller magnitude of the activation entropy (ΔS‡). The activation entropy that is contributed by the temperature-dependent recrossing coefficient (ΔS‡γ) can be expressed as:111
Recently, a complementary experimental analysis of ‘heavy’ EcDHFR was reported.30 Differences in thermal melting temperatures and binding kinetics of certain ligands were observed, leading to the conclusion that the conformational ensemble is altered on protein isotope labelling.30 Furthermore, the intrinsic primary KIE at pH 9 was found to be unchanged by enzyme isotopic substitution at temperatures above 20 °C, but to diverge at low temperatures, giving a ∼2-fold difference between the ‘light’ and ‘heavy’ enzymes at 5 °C.30 Based on a phenomenological modified Marcus model,67 it was suggested that the average distance between hydride donor and acceptor increases at low temperatures in the ‘heavy’ enzyme.30 However, all of these reports conclude that dynamic coupling is not a dominant factor in the DHFR-catalysed reaction.27,30,77,110,111 HIV-1 protease, purine nucleoside phosphorylase, alanine racemase and pentaerythritol tetranitrate reductase have also been used for protein isotope substitution studies.82,106–108 The kinetics were affected differently among these enzymes, implying that there is no universal relationship between the conformational and chemical coordinates. Nevertheless, these experiments demonstrate that protein isotope labelling in combination with computational studies, provide a pragmatic approach to the study of the role of protein dynamics in enzyme catalysis.
To define the direct relation between enzyme motions and the chemical step of catalysis is more challenging. Both BsDHFR111,124,125 and MpDHFR12 are more flexible than EcDHFR, yet at pH 7 the three enzymes have similar single turnover rate constants.12,36 EcDHFR-N23PP/S148A has reduced thermal motion on the millisecond timescale, which impedes the conformational changes required for optimal progression through the catalytic cycle.24 The increased dynamic coupling to the reaction on the femtosecond–picosecond timescale is detrimental to the chemical step as it increases the proportion of unproductive trajectories on the transition state dividing surface.77 Phenomenological models of enzyme catalysis that interpret the KIE in terms of a distribution of donor–acceptor distances cannot account for alternative explanations for these data.67 Dynamic coupling to the reaction coordinate is effectively minimised in TmDHFR, which should provide some benefit, but its catalytic performance is poor due to the dominance of other factors.110 Thermal motions prior to the chemical step may ‘promote’ it by providing an optimal reaction-ready active site configuration, but correlations between motions, and between motions and kinetics, must be interpreted with caution. Although this is most obvious when the timescales of the motions and the chemical events are different, it has recently been demonstrated that many motions in DHFR do not affect the chemical step despite occurring on similar timescales to it.129 Furthermore, even where motions are truly coupled to the reaction coordinate, it does not necessarily follow that these motions drive catalysis since they are controlled by the same free energy landscape.
In the majority of the cases discussed here, the motions are equilibrium fluctuations. ‘Well tuned’ equilibrium motions are a consequence of well-tuned electrostatics. However, this does not subtract from the importance of studies of protein motions in enzyme catalysis. Although statistical motions are a consequence of the free energy landscape of the enzyme rather than a cause of it, and non-statistical dynamics can be affected by statistical motions, their effects can provide valuable information about the physical effects of the free energy surface and changes to it.
It appears that efficient enzymes minimise dynamic coupling in their transition states. The chemical step is a fast step of catalysis and so selective pressure cannot easily act on it. Although most DHFRs give a measurable steady-state KIE at pH 7, the typically >10-fold difference between kcat and kH indicates that the small benefit obtained from eliminating dynamic coupling is unlikely sufficient for selective pressure to act on. Evolution-based arguments may therefore appear inappropriate on first examination. However, evolution does appear to have managed both to minimise dynamic coupling to the reaction coordinate and to maintain this favourable state. One possibility for how this occurred is that enzyme dynamics were optimised early in history, when the chemistry was rate limiting, and neutral drift has not overcome this even when physical steps became rate limiting. This seems an unsatisfactory explanation as it relies on the absence of neutral drift. A more likely alternative may be that coupling of fast dynamics to the reaction coordinate is tied to other factors that selective pressure can act on, such as conformational sampling or other equilibrium motions on millisecond timescales. As all motions are affected by the same free energy surface this seems reasonable, although the question remains how the local non-statistical dynamics, which appear to be uncoupled from motions in the surrounding environment,10,121 are affected by the free energy surface.
These findings lead to the question of why dynamic coupling exists, if it is detrimental to catalysis. Indeed, it is a relatively common phenomenon with a number of enzymes showing a measurable enzyme KIE at or near physiological conditions.27,30,82,106,130 At least in the case of DHFR catalysis, maximal catalytic power appears to be achieved by carefully tuning the involvement of protein motions. Conformational transitions couple to the physical steps, but an ideal reaction-ready active site configuration should provide a static environment, where hydride transfer can proceed efficiently. The dynamic effects observed in DHFR are likely ‘residual’ motions from the reorganisation needed to facilitate the charge transfer in forming the transition state species. As the system climbs the energy barrier towards the transition state, degrees of freedom are progressively lost. This constraint of the atomic positions is likely to translate to a loss of non-statistical motions. Dynamic coupling in DHFR is therefore an undesired by-product that has been imperfectly eliminated. It may be that since the chemical transformation catalysed by DHFR is relatively simple with few charge transfers involved, dynamic effects are automatically minimised so long as the enzyme samples an ideal reaction ready configuration. For more complex, multi-step chemical transformations, dynamic coupling may be less well suppressed. These findings may find uses in designing enzymes with new functionality. To test the validity of the current proposal, investigations of further enzyme systems will be required.
This journal is © the Owner Societies 2015 |