Proofreading through spatial gradients

Version of Record

Accepted for publication after peer review and revision.

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record published: January 18, 2021 (This version)
Accepted: December 24, 2020
Accepted Manuscript published: December 24, 2020 (Go to version)
Received: June 26, 2020

1. Of interest
Substrate evaporation drives collective construction in termites

Giulio Facchini, Alann Rathery ... Andrea Perna

Research Article Apr 10, 2024
Further reading

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Appendix 2
Appendix 3
Appendix 4
Appendix 5
Appendix 6
Data availability
References
Article and author information
Metrics

Abstract

Key enzymatic processes use the nonequilibrium error correction mechanism called kinetic proofreading to enhance their specificity. The applicability of traditional proofreading schemes, however, is limited because they typically require dedicated structural features in the enzyme, such as a nucleotide hydrolysis site or multiple intermediate conformations. Here, we explore an alternative conceptual mechanism that achieves error correction by having substrate binding and subsequent product formation occur at distinct physical locations. The time taken by the enzyme–substrate complex to diffuse from one location to another is leveraged to discard wrong substrates. This mechanism does not have the typical structural requirements, making it easier to overlook in experiments. We discuss how the length scales of molecular gradients dictate proofreading performance, and quantify the limitations imposed by realistic diffusion and reaction rates. Our work broadens the applicability of kinetic proofreading and sets the stage for studying spatial gradients as a possible route to specificity.

Introduction

The nonequilibrium mechanism called kinetic proofreading (Hopfield, 1974; Ninio, 1975) is used for reducing the error rates of many biochemical processes important for cell function (e.g. DNA replication [Kunkel, 2004], transcription [Sydow and Cramer, 2009], translation [Rodnina and Wintermeyer, 2001; Ieong et al., 2016], signal transduction [Swain and Siggia, 2002], or pathogen recognition [McKeithan, 1995; Goldstein et al., 2004; Cui and Mehta, 2018]). Proofreading mechanisms operate by inducing a delay between substrate binding and product formation via intermediate states for the enzyme–substrate complex. Such a delay gives the enzyme multiple chances to release the wrong substrate after initial binding, allowing far lower error rates than what one would expect solely from the binding energy difference between right and wrong substrates.

Traditional proofreading schemes require dedicated molecular features such as an exonuclease pocket in DNA polymerases (Kunkel, 2004) or multiple phosphorylation sites on T-cell receptors (McKeithan, 1995; Goldstein et al., 2004); such features create intermediate states that delay product formation (Figure 1a) and thus allow proofreading. Additionally, since proofreading is an active nonequilibrium process often involving near–irreversible reactions, the enzyme typically needs to have an ATP or GTP hydrolysis site to enable the use of energy supplies of the cell (Yamane and Hopfield, 1977; Rodnina and Wintermeyer, 2001). Due to such stringent structural requirements, the number of confirmed proofreading enzymes is relatively small. Furthermore, generic enzymes without such dedicated features are assumed to not have active error correction available to them.

Figure 1

Download asset Open asset

Error correction schemes that operate by delaying product formation.

(a) The traditional proofreading scheme with multiple biochemically distinct intermediates, transitions between which are typically accompanied by energy–consuming reactions. The T-cell activation mechanism with successive phosphorylation events is used for demonstration (McKeithan, 1995; Cui and Mehta, 2018). (b) The spatial proofreading scheme where the delay between binding and catalysis is created by constraining these events to distinct physical locations. The wavy arrows stand for the diffusive motion of the complex. Binding events primarily take place on the length scale $λ_{_{S}}$ of substrate localization.

In this work, we propose an alternative scheme where the delay between initial substrate binding and product formation steps is achieved by separating these events in space. If substrates are spatially localized and product formation is favorable only in a region of low substrate concentration where an activating effector is present then the time taken by the enzyme–substrate complex to travel from one location to the other can be used to discard the wrong substrates, which are assumed to unbind from the enzyme more readily than the right substrates (Figure 1b). When this delay is longer than substrate unbinding time scales, very low error rates of product formation can be achieved, allowing this spatial proofreading scheme to outperform biochemical mechanisms with a finite number of proofreading steps.

In contrast to traditional proofreading, the nonequilibrium mechanism here does not require any direct energy consumption by the enzyme or substrate itself (e.g. through ATP hydrolysis). This liberates the enzyme from any proofreading-specific molecular features; indeed, any ‘equilibrium’ enzyme with a localized effector can proofread using our scheme if appropriate concentration gradients of the substrates or enzymes are set up. In this way, the energetic and structural requirements of proofreading can be outsourced from the enzyme and substrate to the gradient maintaining mechanism. It also means that spatial proofreading is easy to overlook in experiments, and that the fidelity of reconstituted reactions in vitro could be lower than the fidelity in vivo.

The lack of reliance on structure makes spatial proofreading more adaptable. We study how tuning the length scale of concentration gradients can trade off error rate against speed and energy consumption on the fly. In contrast, traditional proofreading schemes rely on nucleotide chemical potentials, for example, the out of equilibrium [ATP]/[ADP] ratio in the cell, and cannot modulate their operation without broader physiological disruptions.

Our proposed scheme can be leveraged for specificity if appropriate concentration gradients are set. Such gradients arise in multiple cellular contexts (e.g. near the nucleus, the plasma membrane, the Golgi apparatus, the endoplasmic reticulum [ER], kinetochores, microtubules [Bivona et al., 2003; Caudron et al., 2005; Kholodenko, 2006]) and several gradient-forming mechanisms have been discussed in the literature (Wu et al., 2018; Kholodenko, 2006; Kholodenko, 2003). We conclude our analysis of spatial proofreading by quantifying its limitations as set by realistic reaction rates and gradient formation mechanisms, and discuss examples from the literature, including the localization of mRNAs in polarised cells, and the non-vesicular transport of lipids in eukaryotic cells, in which this mechanism might be in play. Our work motivates a detailed investigation of spatial structures and compartmentalization in living cells as possible delay mechanisms for proofreading enzymatic reactions.

Results

Slow transport of enzymatic complex enables proofreading

Our proposed scheme is based on spatially separating substrate binding and product formation events for the enzyme (Figure 1b). Such a setting arises naturally if substrates are spatially localized by having concentration gradients in a cellular compartment. Similarly, an effector needed for product formation (e.g. through allosteric activation) may have a spatial concentration gradient localized elsewhere in that compartment. To keep our model simple, we assume that the right (R) and wrong (W) substrates have identical concentration gradients of length scale $λ_{_{S}}$ but that the effector is entirely localized to one end of the compartment, for example via membrane tethering. In Appendix 4, we extend our study of model performance to the scenario where the two substrates have different localization length scales.

We model our system using coupled reaction–diffusion equations for the substrate-bound (‘ES’ with $S = R, W$ ) and free (‘E’) enzyme densities, namely,

\frac{\partial ρ_{_{ER}}}{\partial t} = D \frac{\partial^{2} ρ_{_{ER}}}{\partial x^{2}} - k_{off}^{R} ρ_{_{ER}} + k_{on} ρ_{_{R}} ρ_{_{E}},

\frac{\partial ρ_{_{EW}}}{\partial t} = D \frac{\partial^{2} ρ_{_{EW}}}{\partial x^{2}} - k_{off}^{W} ρ_{_{EW}} + k_{on} ρ_{_{W}} ρ_{_{E}},

\frac{\partial ρ_{_{E}}}{\partial t} = D \frac{\partial^{2} ρ_{_{E}}}{\partial x^{2}} + \sum_{S = R, W} k_{off}^{S} ρ_{_{ES}} - \sum_{S = R, W} k_{on} ρ_{_{S}} ρ_{_{E}} .

Here, D is the enzyme diffusion constant, $k_{on}$ and $k_{off}^{S}$ (with $k_{off}^{W} > k_{off}^{R}$ ) are the substrate binding and unbinding rates, respectively, and $ρ_{_{S}} (x) \sim e^{- x / λ_{_{S}}}$ is the spatially localized substrate concentration profile which we take to be exponentially decaying, which is often the case for profiles created by cellular gradient formation mechanisms (Driever and Nüsslein-Volhard, 1988; Brown and Kholodenko, 1999). We limit our discussion to this one-dimensional setting of the system, though our treatment can be generalized to two and three dimensions in a straightforward way.

The above model does not explicitly account for several effects relevant to living cells, such as depletion of substrates or distinct diffusion rates for the free and substrate-bound enzymes. More importantly, it does not account for the mechanism of substrate gradient formation. We analyze a biochemically detailed model with this latter feature and experimentally constrained parameters later in the paper. Here, we proceed with the minimal model above for explanatory purposes. To identify the key determinants of the model’s performance, we assume throughout our analysis that the amount of substrates is sufficiently low that the enzymes are mostly free with a roughly uniform profile (i.e. $ρ_{_{E}} \approx constant$ ). This assumption makes Equations (1-3) linear and allows us to solve them analytically at steady state. We demonstrate in Appendix 5 that proofreading is, in fact, most effective under this assumption and discuss the consequences of having high substrate amounts on the performance of the scheme.

In our simplified picture, enzyme activation and catalysis take place upon reaching the right boundary at a rate r that is identical for both substrates. Therefore, the density of substrate–bound enzymes at the right boundary can be taken as a proxy for the rate of product formation $v_{S}$ , since

\begin{matrix} v_{S} = r ρ_{_{ES}} (L), \end{matrix}

where L is the size of the compartment. In order to keep the analytical results concise and intuitive, we perform our main analyses under the assumption that catalysis is slow, mirroring the study of traditional proofreading schemes (Hopfield, 1974). In Appendix 3, we derive the precise conditions under which this treatment is valid, and generalize our analysis to arbitrary catalysis rates.

To demonstrate the proofreading capacity of the model, we first analyze the limiting case where substrates are localized to the left end of the compartment ( $λ_{_{S}} \to 0$ ). In this limit, the fidelity $η$ , defined as the number of right products formed per single wrong product, becomes

\begin{matrix} η = \frac{v_{R}}{v_{W}} = \sqrt{η_{eq}} \frac{\sinh (\sqrt{τ_{D} k_{off}^{W}})}{\sinh (\sqrt{τ_{D} k_{off}^{R}})}, \end{matrix}

where $η_{eq} = k_{off}^{W} / k_{off}^{R}$ is the equilibrium fidelity, and $τ_{D} = L^{2} / D$ is the characteristic time scale of diffusion across the compartment (see Appendix 1 for the derivation).

Equation 5 is plotted in Figure 2 for a family of different parameter values. As can be seen, when diffusion is fast (small $τ_{D}$ ), fidelity converges to its equilibrium value and proofreading is lost ( $η \approx \sqrt{η_{eq}} \times \sqrt{τ_{D} k_{off}^{W} / τ_{D} k_{off}^{R}} = η_{eq}$ ). Conversely, when diffusion is slow (large $τ_{D}$ ), the enzyme undergoes multiple rounds of binding a substrate at the left end and unbinding midway until it manages to diffuse across the whole compartment as a complex and form a product. These rounds serve as ‘futile cycles’ that endow the system with proofreading. In this regime, fidelity scales as

η \sim e^{(\sqrt{k_{off}^{W}} - \sqrt{k_{off}^{R}}) \sqrt{τ_{D}}} .

Figure 2

Download asset Open asset

Dependence of fidelity on the diffusion time scale in the limit of very high substrate localization.

Individual curves were made for different choices of $k_{off}^{W}$ (varied in the $[10 - 100] k_{off}^{R}$ range). $τ_{off}^{R} = 1 / k_{off}^{R}$ is the unbinding time scale of right substrates, kept fixed in the study. Fidelity values corresponding to integer degrees of proofreading in a traditional sense ( $η / η_{eq} = η_{eq}^{n}$ , $n = 1, 2, 3, \dots$ ) are marked as circles. Dominant processes in the two limiting regimes are highlighted in red in the schematics shown as insets.

To get further insights, we introduce an effective number of extra biochemical intermediates (n) that a traditional proofreading scheme would need to have in order to yield the same fidelity, that is $η / η_{eq} = η_{eq}^{n}$ . We calculate this number as (see Appendix 1)

\begin{matrix} n \approx \frac{\sqrt{τ_{D} k_{off}^{W}}}{\ln η_{eq}} . \end{matrix}

Notably, since $τ_{D} \sim L^{2}$ , the result above suggests a linear relationship between the effective number of proofreading realizations and the compartment size ( $n \sim L$ ). In addition, because the right-hand side of Equation 7 is an increasing function of $k_{off}^{W}$ , the proofreading efficiency of the scheme rises with larger differences in substrate off-rates (Figure 2) – a feature that ‘hard–wired’ traditional proofreading schemes with a fixed number of proofreading steps lack.

Navigating the speed–fidelity trade-off

As is inherent to all proofreading schemes, the fidelity enhancement described earlier comes at a cost of reduced product formation speed. This reduction, in our case, happens because of increased delays in diffusive transport. Here, we explore the resulting speed–fidelity trade-off and its different regimes by varying two of the model parameters: diffusion time scale $τ_{D}$ and the substrate localization length scale $λ_{_{S}}$ .

Speed and fidelity for different sampled values of $τ_{D}$ and $λ_{_{S}}$ are depicted in Figure 3a. As can be seen, for a fixed $τ_{D}$ , the reduction of $λ_{_{S}}$ can trade off fidelity against speed. This trade-off is intuitive; with tighter substrate localization, the complexes are formed closer to the left boundary. Hence, a smaller fraction of complexes reach the activation region, reducing reaction speed. The Pareto-optimal front of the trade-off over the whole parameter space, shown as a red curve on the plot, is reached in the limit of ideal substrate localization ( $λ_{_{S}} \to 0$ ). Varying the diffusion time scale allows one to navigate this optimal trade-off curve and access different performance regimes.

Figure 3

Download asset Open asset

Speed–fidelity trade-off and consequences of having weak substrate gradients.

(a) Speed and fidelity evaluated for sampled values of the diffusion time scale ( $τ_{D}$ ) and substrate localization length scale ( $λ_{_{S}}$ ). Here, $v_{eq} \sim 1 / k_{off}^{R}$ is the speed in the equilibrium limit of a uniform substrate profile ( $λ_{_{S}} \to \infty$ ). The red line corresponds to the Pareto-optimal front and is reached in the high substrate localization limit. The example speed–fidelity trade-off illustrated through the black dotted curve is obtained for $τ_{D} \approx 20 τ_{off}^{R}$ . (b) Density profiles of wrong (EW) and right (ER) complexes in three qualitatively different performance regimes. The normalization factor $ρ_{_{ES}}^{eq}$ corresponds to the equilibrium complex densities. (c) Fidelity as a function of diffusion time scale for different choices of $λ_{_{S}}$ (varied in the $[0.04, 0.4] L$ range). The dashed line corresponds to the ideal substrate localization limit ( $λ_{_{S}} \to 0$ ). Inset: Fidelity as a function of $L / λ_{_{S}}$ for a fixed $τ_{D}$ . Shaded area indicates the range where the bulk of fidelity enhancement takes place. Equilibrium fidelity $η_{eq} = 10$ was used in generating all the panels.

Specifically, if the diffusion time scale is fast compared with the time scales of substrate unbinding (i.e. $τ_{D} ≪ 1 / k_{off}^{R}, 1 / k_{off}^{W}$ ), then both right and wrong complexes that form near the left boundary arrive at the activation region with high probability, resulting in high speeds, although at the expense of error–prone product formation (Figure 3b, top). In the opposite limit of slow diffusion, both types of complexes have exponentially low densities at the activation region, but due to the difference in substrate off-rates, production is highly accurate (Figure 3b, bottom). There also exists an intermediate regime where a significant fraction of right complexes reach the activation region while the vast majority of wrong complexes do not (Figure 3b, middle). As a result, an advantageous trade-off is achieved where a moderate decrease in the production rate yields high fidelity enhancement – a feature that was also identified in multi-step traditional proofreading models (Murugan et al., 2012).

In Appendix 3, we also study this trade-off caused by varying the catalysis rate $r$ . Briefly, we find that when all other parameters are fixed, increasing $r$ trades off fidelity against speed in a linear fashion, with the ratio of highest and lowest fidelity values falling in the $[\sqrt{η_{eq}}, η_{eq}]$ range. The Pareto–optimal front of the trade-off, however, monotonically shifts toward the higher speed region, suggesting that faster catalysis is, in fact, more favorable if the diffusion time scale $τ_{D}$ can be adjusted accordingly (see Appendix 3 for details).

We saw in Figure 3a that in the case of ideal substrate localization, the slowdown of diffusive transport necessarily reduced the production rate and increased the fidelity. The latter part of this statement, however, breaks down when substrate gradients are weak. Indeed, fidelity exhibits a non-monotonic response to tuning $τ_{D}$ when the substrate gradient length scale $λ_{_{S}}$ is non-zero (Figure 3c). The reason for the eventual decay in fidelity is the fact that with slower diffusion (larger $τ_{D}$ ), substrate binding and unbinding events take place more locally and therefore, the right and wrong complex profiles start to resemble the substrate profile itself, which does not discriminate between the two substrate kinds. We show in Appendix 1 that the optimal diffusion time scale can be roughly approximated as $τ_{D}^{*} / τ_{off}^{R} \approx η_{eq}^{- 1} {(L / λ_{_{S}})}^{2}$ , which increases monotonically with $L / λ_{_{S}}$ , consistent with the shifting peaks in Figure 3c.

Not surprisingly, the error–correcting capacity of the scheme improves with better substrate localization (lower $λ_{_{S}}$ ). For a fixed $τ_{D}$ , the bulk of this improvement takes place when $L / λ_{_{S}}$ is tuned in a range set by the two key dimensionless numbers of the model, namely, $\sqrt{τ_{D} k_{off}^{R}}$ and $\sqrt{τ_{D} k_{off}^{W}}$ (Figure 3c, inset). In Appendix 1, we provide an analytical justification for this result. Taken together, these parametric studies uncover the operational principles of the spatial proofreading scheme and demonstrate how the speed–fidelity trade-off could be dynamically navigated as needed by tuning the key time and length scales of the model.

Energy dissipation and limits of proofreading performance

A hallmark signature of proofreading is that it is a nonequilibrium mechanism with an associated free energy cost. In our scheme, the enzyme itself is not directly involved in any energy-consuming reactions, such as hydrolysis. Instead, the free energy cost comes from maintaining the spatial gradient of substrates, which the enzymatic reaction tends to homogenize by releasing bound substrates in regions of low substrate concentration. As the activating effectors are assumed to be tethered at $x = L$ , they do not require dissipation to remain localized.

While mechanisms of substrate gradient maintenance may differ in their energetic efficiency, there exists a thermodynamically dictated minimum energy that any such mechanism must dissipate per unit time. We calculate this minimum power P as

\begin{matrix} P & = \sum_{S = {R,W}} \int_{0}^{L} j_{_{S}} (x) μ (x) d x . \end{matrix}

Here $j_{_{S}} (x) = k_{on} ρ_{_{S}} (x) ρ_{_{E}} - k_{off}^{S} ρ_{_{ES}} (x)$ is the net local binding flux of substrate ‘S’, and $μ (x)$ is the local chemical potential (see Appendix 2.1 for details). For substrates with an exponentially decaying profile considered here, the chemical potential is given by

\begin{matrix} μ (x) = μ (0) + k_{B} T \ln \frac{ρ_{_{S}} (x)}{ρ_{_{S}} (0)} = μ (0) - k_{B} T \frac{x}{λ_{_{S}}}, \end{matrix}

where $k_{B} T$ is the thermal energy scale. Notably, the chemical potential difference across the compartment, which serves as an effective driving force for the scheme, is set by the inverse of the nondimensionalized substrate localization length scale, namely,

\begin{matrix} β Δ μ = \frac{L}{λ_{_{S}}}, \end{matrix}

where $β^{- 1} = k_{B} T$ . This driving force is zero for a uniform substrate profile ( $λ_{_{S}} \to \infty$ ) and increases with tighter localization (lower $λ_{_{S}}$ ), as intuitively expected.

We used Equation 8 to study the relationship between dissipation and fidelity enhancement as we tuned $Δ μ$ for different choices of the diffusion time scale $τ_{D}$ . As can be seen in Figure 4, power rises with increasing fidelity, diverging when fidelity reaches its asymptotic maximum given by Equation 5 in the large $Δ μ$ limit. For the bulk of each curve, power scales as the logarithm of fidelity, suggesting that a linear increase in dissipation can yield an exponential reduction in error. Notably, such a scaling relationship has also been calculated in the context of E. coli chemoreceptor adaptation (Lan et al., 2012). In particular, it was shown that the adaptation error decreases exponentially with energy dissipated through multiple methylation–demethylation cycles which are used to stabilize the activity state of the receptor. Analogies in the cost-performance trade-off across these functionally distinct mechanisms contribute to the search for overarching thermodynamic themes underlying cellular information processing (Lan et al., 2012; Lan and Tu, 2013; Horowitz et al., 2017; Sartori and Pigolotti, 2015).

Figure 4

Download asset Open asset

Power–fidelity relationship when tuning the effective driving force $Δ μ$ for different choices of the diffusion time scale $τ_{D}$ .

$J_{b i n d} = k_{o n} ρ_{E} \int ρ_{S} (x) d x$ is the integrated rate of substrate binding. The red line indicates the large dissipation limit of fidelity given by Equation 5. The circles indicate the $Δ μ$ range specified in Equation 11 for different $τ_{D}$ choices. For sufficiently large $τ_{D}$ values, the cost per binding event approaches $β \sqrt{η_{eq}}$ at the end of this range (see Appendix 2.1 for details). In making this plot, $η_{eq} = 10$ was used.

The logarithmic scaling is achieved in our model when the driving force is in a range where most of the fidelity enhancement takes place, namely,

β Δ μ \in [\sqrt{τ_{D} k_{off}^{R}}, \sqrt{τ_{D} k_{off}^{W}}] .

At the end of this range, the cost per substrate binding event approaches $\sqrt{η_{eq}}$ in $k_{B} T$ units (see Appendix 2.1 for details). And beyond the range, additional error correction is attained at an increasingly higher cost.

Note that the power computed here does not include the baseline cost of creating the substrate gradient, which, for instance, would depend on the substrate diffusion constant. We only account for the additional cost to be paid due to the operation of the proofreading scheme which works to homogenize this substrate gradient. The baseline cost in our case is analogous to the work that ATP synthase needs to perform to maintain a nonequilibrium [ATP]/[ADP] ratio in the cell, whereas our calculated power is analogous to the rate of ATP hydrolysis by a traditional proofreading enzyme. We discuss these two classes of dissipation in greater detail in Appendix 2.3.

Just as the cellular chemical potential of ATP or GTP imposes a thermodynamic upper bound on the fidelity enhancement by any proofreading mechanism (Qian, 2006), the effective driving force $Δ μ$ imposes a similar constraint for the spatial proofreading model. This thermodynamic limit depends only on the available chemical potential and is equal to $e^{β Δ μ}$ . This limit can be approached very closely by our model, which for $Δ μ ≳ 1$ achieves the exponential enhancement with an additional linear prefactor, namely, ${(η / η_{eq})}^{max} \approx e^{β Δ μ} / β Δ μ$ (see Appendix 2.2). Such scaling behavior was theoretically accessible only to infinite-state traditional proofreading schemes (Qian, 2006; Ehrenberg and Blomberg, 1980). This offers a view of spatial proofreading as a procession of the enzyme through an infinite series of spatial filters and suggests that, from the perspective of peak error reduction capacity, our model outperforms the finite-state schemes.

Proofreading by biochemically plausible intracellular gradients

Our discussion of the minimal model thus far was not aimed at a particular biochemical system and thus did not involve the use of realistic reaction rates and diffusion constants typically seen in living cells. Furthermore, we did not account for the possibility of substrate diffusion, as well as for the homogenization of substrate concentration gradients due to enzymatic reactions, and have thereby abstracted away the gradient maintaining mechanism. The quantitative inspection of such mechanisms is important for understanding the constraints on spatial proofreading in realistic settings.

Here, we investigate proofreading based on a widely applicable mechanism for creating gradients by the spatial separation of two opposing enzymes (Stelling and Kholodenko, 2009; Bivona et al., 2003; Brown and Kholodenko, 1999). Consider a protein $S$ that in its free state is phosphorylated by a membrane-bound kinase and dephosphorylated by a delocalized cytoplasmic phosphatase, as shown in Figure 5a. This setup will naturally create a gradient of the active form of protein ( $S^{*}$ ), with the gradient length scale controlled by the rate of phosphatase activity $k_{p}$ ( $S^{*} \overset{k_{p}}{\to} S$ ). Such mechanisms are known to create gradients of the active forms of MEK and ERK (Kholodenko, 2006), of GTPases such as Ran (with GEF and GAP [Kalab et al., 2002] playing the role of kinase and phosphatase, respectively), of cAMP (Kholodenko, 2006) and of stathmin oncoprotein 18 (Op18) (Bastiaens et al., 2006; Niethammer et al., 2004) near the plasma membrane, the Golgi apparatus, the ER, kinetochores and other places.

Figure 5

Download asset Open asset

Proofreading based on substrate gradients formed by spatially separated kinases and phosphatases.

(a) The active form $S^{*}$ of many proteins exhibits gradients because kinases that phosphorylate $S$ are anchored to a membrane while phosphatases can diffuse in the cytoplasm (Kholodenko, 2006). An enzyme can exploit the resulting spatial gradient for proofreading. (b) At low enzyme activity (i.e. low $k_{on} ρ_{_{E}}$ ), the gradient of $S^{*}$ is successfully maintained, allowing for proofreading. The upper dashed line corresponds to the peak fidelity when the substrate profile is exponential. At high enzyme activity (large $k_{on} ρ_{_{E}}$ ), the dephosphorylation with rate $k_{p} = 5$ s^-1 is no longer sufficient to maintain the gradient and proofreading is lost. (c) Profiles of right substrates for different choices of enzyme activity. Numbers indicate $k_{on} ρ_{_{E}}$ in s^-1 units. The black line shows an exponential substrate profile with a length scale $λ_{_{S}} = \sqrt{D / k_{p}} \sim 0.5$ μm.

We test the proofreading power of such gradients, assuming experimentally constrained biophysical parameters for the gradient forming mechanism. Specifically, we consider an enzyme $E$ that acts on the active forms of cognate ( $R^{*}$ ) and non-cognate ( $W^{*}$ ) substrates which have off-rates 0.1 s^-1 and 1 s^-1, respectively (hence, $η_{eq} = 10$ ). These off-rates are consistent with typical values for substrates proofread by cellular signaling systems (Cui and Mehta, 2018; Gascoigne et al., 2001). The kinases and phosphatases in our setup act identically on right and wrong substrates. We consider a dephosphorylation rate constant $k_{p} = 5$ s^-1 that falls in the range 0.1−100 s^-1 reported for different phosphatases (Brown and Kholodenko, 1999; Kholodenko et al., 2000; Todd et al., 1999), and a cytosolic diffusion constant $D = 1$ μm²/s for all proteins in this model. With this setup, exponential gradients of length scale ∼0.5 μm are formed for $R^{*}$ and $W^{*}$ . We evaluate the proofreading and energetic performance of the model in a compartment of size $L = 10$ μm – a typical length scale in eukaryotic cells (see Appendix 6 for details).

Although not cost-efficient, this setup achieves proofreading in a wide range of regimes. Specifically, it is most effective when the enzyme–substrate binding is slow, in which case the exponential substrate profile is maintained and the system attains the fidelity predicted by our earlier explanatory model (Figure 5b). The system’s proofreading capacity is retained if the first–order on-rate is raised up to $k_{on} ρ_{_{E}} \sim 10$ s^-1, where around 10-fold increase in fidelity is still possible. If the binding rate constant ( $k_{on}$ ) or the enzyme’s expression level ( $ρ_{_{E}}$ ) is any higher, then enzymatic reactions overwhelm the ability of the kinase/phosphatase system to keep the active forms of substrates sufficiently localized (Figure 5c) and proofreading is lost. Overall, this model suggests that enzymes can work at reasonable binding rates and still proofread, when accounting for an experimentally characterized gradient maintaining mechanism.

Discussion

We have outlined a way for enzymatic reactions to proofread and improve specificity by exploiting spatial concentration gradients of substrates. Like the classic model, our proposed spatial proofreading scheme is based on a time delay; but unlike the classic model, here the delay is due to spatial transport rather than transitions through biochemical intermediates. Consequently, the enzyme is liberated from the stringent structural requirements imposed by traditional proofreading, such as multiple intermediate conformations and hydrolysis sites for energy coupling. Instead, our scheme exploits the free energy supplied by active mechanisms that maintain spatial structures.

The decoupling of the two crucial features of proofreading – time delay and free energy dissipation – allows the cell to tune proofreading on the fly. For instance, all proofreading schemes offer fidelity at the expense of reaction speed and energy. For traditional schemes, navigating this trade-off is not always feasible, as it needs to involve structural changes via mutations or modulation of the [ATP]/[ADP] ratio which can cause collateral effects on the rest of the cell. In contrast, the spatial proofreading scheme is more adaptable to the changing conditions and needs of the cell. The scheme can prioritize speed in one context, and fidelity in another, simply by tuning the length scale of intracellular gradients (e.g. through the regulation of the phosphotase or free enzyme concentration in the scheme discussed earlier).

On the other hand, this modular decoupling can complicate the experimental identification of proofreading enzymes and the interpretation of their fidelity. Here, the enzymes need not be endowed with the structural and biochemical properties typically sought for in a proofreading enzyme. At the same time, any attempt to reconstitute enzymatic activity in a well-mixed, in vitro assay, will show poor fidelity compared to in vivo measurements, even when all necessary molecular players are present in vitro. Therefore, more care is required in studies of cellular information processing mechanisms that hijack a distant source of free energy compared to the case where the relevant energy consumption is local and easier to link causally to function.

While we focused on spatially localized substrates and delocalized enzymes, our framework would apply equally well to other scenarios, like one with a spatially localized enzyme (or its active form [Kalab et al., 2002; Nalbant et al., 2004]) and effector with delocalized substrates, an example of which would be an alternative version of the scheme in Figure 5a where the target of the kinase/phosphatase activity is changed from substrates to enzymes. Our framework can also be extended to signaling cascades, where slightly different phosphatase activities can result in magnified concentration ratios of two competing signaling molecules at the spatial location of the next cascade step (Roy and Cyert, 2009; Bauman and Scott, 2002; Kholodenko, 2006).

The spatial gradients needed for the operation of our model can be created and maintained through multiple mechanisms in the cell, ranging from the kinase/phosphatase system modeled here, to the passive diffusion of substrates/ligands combined with active degradation (e.g. Bicoid and other developmental morphogens), to active transport processes combined with diffusion. A particularly simple implementation of our scheme is via compartmentalization – substrates and effectors are localized in two spatially separated compartments with the enzyme–substrate complex having to travel from one to another to complete the reaction.

Many molecular localization pathways involving the naturally compartmentalized parts of the cell require high substrate selectivity and are therefore potential candidates for the implementation of spatial proofreading. For example, in polarized, asymmetric cells (e.g. budding yeast or neuronal cells) gene expression often needs to be spatially regulated (Parton et al., 2014; Martin and Ephrussi, 2009). Such regulation is achieved with designated ribonucleoproteins that bind specific mRNAs near the cell nucleus, perform a biased random walk to the mRNA localization site and deliver them for translation. During transport, mRNAs are protected from ribosome binding and when they unbind, they are subject to degradation which would prevent rebinding events at intermediate locations. Another example process is the non-vesicular transport of lipids between the membrane–bound domains of the cells (e.g. the ER, mitochondria, the Golgi apparatus, or the plasma membrane). This transport mechanism is mediated by lipid-transfer proteins that bind lipids on the donor membrane, diffuse to the acceptor membrane and upon interacting with it, undergo a conformational change, delivering the ‘cargo’ (Lev, 2010). Although the higher proximity of the two membranes is thought to enhance the transport efficiency, it would be interesting to study the optimality of the inter-membrane distance in the context of fidelity–transport efficiency trade-off, given the fact that some of the lipid-transfer proteins are known to exhibit specificity for their cognate substrates.

Our scheme may also be applicable as a quality control mechanism in protein secretion pathways (Ellgaard and Helenius, 2003; Arvan et al., 2002), in high-fidelity targeting of membrane proteins mediated by signal recognition particles (Rao et al., 2016; Chio et al., 2017), as well as in selective glycosylation reactions in the Golgi apparatus (Jaiman and Thattai, 2020). Lastly, considering the recent advances in generating synthetic morphogen patterns in multicellular organisms (Toda et al., 2020; Stapornwongkul et al., 2020), spatial proofreading could also be employed in pathways acting on engineered protein gradients. Experimental investigations of these processes in light of our work will reveal the extent to which spatial transport promotes specificity.

In conclusion, we have analyzed the role played by spatial structures in endowing enzymatic reactions with kinetic proofreading. Simply by spatially segregating substrate binding from catalysis, enzymes can enhance their specificity. This suggests that enzymatic reactions may acquire de novo proofreading capabilities by coupling to pre-existing spatial gradients in the cell.

Materials and methods

Detailed derivations of the analytical results presented in the main text along with additional studies on our model are included in the Appendices. In addition, Python scripts and Jupyter notebooks used to generate all the plots in the main text and Appendices are included as Supplementary files.

Appendix 1

Analytical calculations of the complex density profile and fidelity

We begin this section by deriving an analytical expression for the density profile of substrate-bound enzymes ( $ρ_{_{ES}} (x)$ ) in the case where the $ρ (x) \approx constant$ assumption holds. Based on this result, we then obtain expressions for fidelity in low, high, and intermediate substrate localization regimes. We reserve the studies of speed and fidelity in the general case of a nonuniform free enzyme profile to Appendix 5.

1. Derivation of the complex density profile $ρ_{_{ES}} (x)$

The ordinary differential equation (ODE) that defines the steady state profile of substrate-bound enzymes is

\begin{matrix} \underset{diffusion}{\underset{⏟}{D \frac{d^{2} ρ_{_{ES}}}{d x^{2}}}} - \underset{unbinding}{\underset{⏟}{k_{off}^{S} ρ_{_{ES}} (x)}} + \underset{binding}{\underset{⏟}{k_{on} ρ_{_{S}} (0) e^{- x / λ_{_{S}}} ρ_{_{E}} (x)}} = 0 . \end{matrix}

Here, $ρ_{_{S}} (0)$ is the substrate density at the leftmost boundary, whose value can be calculated from the condition that the total number of free substrates is $S_{total}$ , namely,

\begin{array}{ll} S_{total} & = \int_{x = 0}^{L} ρ_{_{S}} (0) e^{- x / λ_{_{S}}} d x \\ = ρ_{_{S}} (0) λ_{_{S}} (1 - e^{- L / λ_{_{S}}}) \Rightarrow \end{array}

ρ_{S} (0) = \frac{S_{t o t a l}}{λ_{S} (1 - e^{- L / λ_{S}})} .

In the limit of low substrate amounts where the approximation $ρ_{_{E}} (x) \approx constant$ is valid, Equation S1 represents a linear nonhomogeneous ODE. Hence, its solution can be written as

\begin{matrix} ρ_{_{ES}} (x) = ρ_{_{ES}}^{(h)} (x) + ρ_{_{ES}}^{(p)} (x), \end{matrix}

where $ρ_{_{ES}}^{(h)} (x)$ is the general solution to the corresponding homogeneous equation, while $ρ_{_{ES}}^{(p)} (x)$ is a particular solution.

Looking for solutions of the form $C e^{- x / λ}$ for the homogeneous part, we find

\begin{matrix} C (\frac{D}{λ^{2}} - k_{off}^{S}) e^{- x / λ} = 0 . \end{matrix}

The two possible roots for $λ$ are $\pm \sqrt{D / k_{off}^{S}}$ . Calling the positive root $λ_{_{ES}}$ , which represents the mean distance traveled by the substrate–bound enzyme before releasing the substrate, we can write the general solution to the homogeneous part of Equation S1 as

\begin{matrix} ρ_{_{ES}}^{(h)} (x) = C_{1} e^{- x / λ_{_{ES}}} + C_{2} e^{x / λ_{_{ES}}}, \end{matrix}

where C₁ and C₂ are constants which will be determined from the boundary conditions.

Since the nonhomogeneous part of Equation S1 is a scaled exponential, we look for a particular solution of the same functional form, namely, $ρ_{_{ES}}^{(p)} (x) = C_{p} e^{- x / λ_{_{S}}}$ . Substituting this form into the ODE, we obtain

\begin{matrix} C_{p} (\frac{D}{λ_{_{S}}^{2}} - k_{off}^{S}) e^{- x / λ_{_{S}}} = - k_{on} ρ_{_{S}} (0) e^{- x / λ_{_{S}}} ρ_{_{E}} . \end{matrix}

The constant coefficient $C_{p}$ can then be found as

\begin{aligned} C_{p} & = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S} - \frac{D}{λ_{_{S}}^{2}}} = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S} (1 - \frac{D / k_{off}^{S}}{λ_{_{S}}^{2}})} \\ = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S} (1 - \frac{λ_{_{ES}}^{2}}{λ_{_{S}}^{2}})}, \end{aligned}

where we have used the equality $λ_{_{ES}} = \sqrt{D / k_{off}^{S}}$ .

Now, to find the unknown coefficients C₁ and C₂, we impose the no-flux boundary conditions for the density $ρ_{E S} (x)$ at the left and right boundaries of the compartment, namely,

\frac{d ρ_{_{E S}}}{d x} |_{x = 0} = - \frac{C_{1}}{λ_{_{E S}}} + \frac{C_{2}}{λ_{_{E S}}} - \frac{C_{p}}{λ_{_{S}}} = 0,

\frac{d ρ_{_{E S}}}{d x} |_{x = L} = - \frac{C_{1}}{λ_{_{E S}}} e^{- \frac{L}{λ_{_{E S}}}} + \frac{C_{2}}{λ_{_{E S}}} e^{\frac{L}{λ_{_{E S}}}} - \frac{C_{p}}{λ_{_{S}}} e^{- \frac{L}{λ_{_{S}}}} = 0.

Note that we did not take into account the product formation flux at the rightmost boundary when writing Equation S10 in order to simplify our calculations. This is justified in the limit of slow catalysis – an assumption that we make in our treatment. The above system of two equations can then be solved for C₁ and C₂, yielding

C_{1} = - \frac{λ_{E S}}{2 λ_{S}} \frac{e^{L / λ_{E S}} - e^{- L / λ_{S}}}{\sinh (L / λ_{E S})} C_{p},

C_{2} = \frac{λ_{E S}}{2 λ_{S}} \frac{e^{- L / λ_{S}} - e^{- L / λ_{E S}}}{\sinh (L / λ_{E S})} C_{p} .

With the constant coefficients known, we obtain the general solution for the complex profile as

\begin{array}{ll} ρ_{_{ES}} (x) & = C_{1} e^{- x / λ_{_{ES}}} + C_{2} e^{x / λ_{_{ES}}} + C_{p} e^{- x / λ_{_{S}}} \\ = C_{p} (\frac{λ_{_{ES}}}{λ_{_{S}} \sinh (L / λ_{_{ES}})} [- \frac{e^{(L - x) / λ_{_{ES}}} + e^{(x - L) / λ_{_{ES}}}}{2} + \frac{e^{- x / λ_{_{ES}}} + e^{x / λ_{_{ES}}}}{2} e^{- L / λ_{_{S}}}] + e^{- x / λ_{_{S}}}) \\ = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S} (1 - λ_{_{ES}}^{2} / λ_{_{S}}^{2})} (\frac{λ_{_{ES}}}{λ_{_{S}} \sinh (L / λ_{_{ES}})} [- \cosh (\frac{L - x}{λ_{_{ES}}}) + \cosh (\frac{x}{λ_{_{ES}}}) e^{- L / λ_{_{S}}}] + e^{- x / λ_{_{S}}}) \\ = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S} (1 - λ_{_{ES}}^{2} / λ_{_{S}}^{2})} (\frac{λ_{_{ES}}}{λ_{_{S}} \sinh (L / λ_{_{ES}})} [- \cosh (\frac{L - x}{λ_{_{ES}}}) + \cosh (\frac{x}{λ_{_{ES}}}) e^{- L / λ_{_{S}}}] + e^{- x / λ_{_{S}}}) . \end{array}

2. Density profile in low and high substrate localization regimes

If substrate localization is very poor ( $λ_{_{S}} ≫ L$ ), the substrate distribution will be uniform ( $ρ_{_{S}} (x) = {\bar{ρ}}_{_{S}} = S_{total} / L$ ), resulting in a similarly flat profile of enzyme–substrate complexes with their density $ρ_{_{ES}}^{\infty}$ given by

\begin{matrix} ρ_{_{ES}}^{\infty} & = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S}} \\ = \frac{k_{on} {\bar{ρ}}_{_{S}} ρ_{_{E}}}{k_{off}^{S}} . \end{matrix}

This is the expected equilibrium result where the complex concentration is inversely proportional to the dissociation constant ( $k_{off}^{S} / k_{on}$ ).

In the opposite limit where the substrates are highly localized ( $λ_{_{S}} ≪ λ_{_{ES}}, L$ and $ρ_{_{S}} (0) \approx S_{total} / λ_{_{S}}$ from Equation S3), the complex density profile simplifies into

\begin{aligned} ρ_{_{ES}} (x) & \approx \frac{k_{on} S_{total} ρ_{_{E}}}{k_{off}^{S} λ_{_{S}} (- λ_{_{ES}}^{2} / λ_{_{S}}^{2})} (- \frac{λ_{_{ES}}}{λ_{_{S}} \sinh (L / λ_{_{ES}})} \cosh (\frac{L - x}{λ_{_{ES}}})) \\ = \frac{k_{on} S_{total} ρ_{_{E}}}{k_{off}^{S} L} \frac{L / λ_{_{ES}}}{\sinh (L / λ_{_{ES}})} \cosh (\frac{L - x}{λ_{_{ES}}}) \\ = ρ_{_{ES}}^{\infty} \times \frac{L / λ_{_{ES}}}{\sinh (L / λ_{_{ES}})} \cosh (\frac{L - x}{λ_{_{ES}}}) . \end{aligned}

The x-dependence through the $\cosh (\cdot)$ function suggests that the complex density is the highest at the leftmost boundary and lowest at the rightmost boundary, with the degree of complex localization dictated by the length scale parameter $λ_{_{ES}}$ . Notably, this localization of complexes does not alter their total number, since the average complex density is conserved, that is,

\begin{aligned} ⟨ ρ_{_{ES}} ⟩ & = \int_{0}^{L} ρ_{_{ES}} (x) d x \\ = ρ_{_{ES}}^{\infty} \times \frac{L / λ_{_{ES}}}{\sinh (L / λ_{_{ES}})} \times \frac{1}{L} \int_{0}^{L} \cosh (\frac{L - x}{λ_{_{ES}}}) d x \\ = ρ_{_{ES}}^{\infty} \times \frac{L / λ_{_{ES}}}{\sinh (L / λ_{_{ES}})} \times \frac{λ_{_{ES}}}{L} \sinh (L / λ_{_{ES}}) \\ = ρ_{_{ES}}^{\infty} . \end{aligned}

Equation S15 for the complex profile can be alternatively written in terms of the diffusion time scale $τ_{D} = L^{2} / D$ and the substrate off-rate $k_{off}^{S}$ . Noting that $L / λ_{_{ES}} = \sqrt{L^{2} k_{off}^{S} / D} = \sqrt{τ_{D} k_{off}^{S}}$ and introducing a dimensionless coordinate $\tilde{x} = x / L$ , we find

\begin{matrix} ρ_{_{ES}} (x) = ρ_{_{ES}}^{\infty} \times \frac{\sqrt{τ_{D} k_{off}^{S}}}{\sinh (\sqrt{τ_{D} k_{off}^{S}})} \cosh (\sqrt{τ_{D} k_{off}^{S}} (1 - \tilde{x})) . \end{matrix}

The above equation is what was used for generating the plots in Figure 3b of the main text for different choices of the diffusion time scale.

3. Fidelity in low and high substrate localization regimes

Let us now evaluate the fidelity of the model in the two limiting regimes discussed earlier. In the poor substrate localization case, which corresponds to an equilibrium setting, the fidelity can be found from Equation S14 as

η_{eq} = \frac{r ρ_{E R}^{\infty}}{r ρ_{E W}^{\infty}} = \frac{k_{off}^{W}}{k_{off}^{R}},

where we have employed the assumption about the right and wrong substrates having identical density profiles. This is the expected result for equilibrium discrimination where no advantage is taken of the system’s spatial structure.

In the regime with high substrate localization, the enzyme–substrate complexes have a nonuniform spatial distribution. What matters for product formation is the complex density at the rightmost boundary ( $\tilde{x} = 1$ ), which we obtain from Equation S17 as

\begin{matrix} ρ_{_{ES}} (L) = ρ_{_{ES}}^{\infty} \times \frac{\sqrt{τ_{D} k_{off}^{S}}}{\sinh (\sqrt{τ_{D} k_{off}^{S}})} . \end{matrix}

Substituting the above expression written for right and wrong complexes into the definition of fidelity, we find

\begin{array}{ll} η & = \frac{r ρ_{_{ER}} (L)}{r ρ_{_{EW}} (L)} \\ = η_{eq} \times \sqrt{\frac{k_{off}^{R}}{k_{off}^{W}}} \frac{\sinh (\sqrt{τ_{D} k_{off}^{W}})}{\sinh (\sqrt{τ_{D} k_{off}^{R}})} \\ = \sqrt{η_{eq}} \frac{\sinh (\sqrt{τ_{D} k_{off}^{W}})}{\sinh (\sqrt{τ_{D} k_{off}^{R}})} . \end{array}

This is the result reported in Equation 5 of the main text. To gain more intuition about it and draw parallels with traditional kinetic proofreading, let us consider the limit of long diffusion time scales where proofreading is the most effective. In this limit, the hyperbolic sine functions above can be approximated as $\sinh (\sqrt{τ_{D} k_{off}^{S}}) \approx 0.5 e^{\sqrt{τ_{D} k_{off}^{S}}}$ , simplifying the fidelity expression into

\begin{aligned} η & = \sqrt{η_{eq}} \frac{e^{\sqrt{τ_{D} k_{off}^{W}}}}{e^{\sqrt{τ_{D} k_{off}^{R}}}} \\ = \sqrt{η_{eq}} e^{\sqrt{τ_{D} k_{off}^{W}} - \sqrt{τ_{D} k_{off}^{R}}} \\ = \sqrt{η_{eq}} e^{\sqrt{τ_{D} k_{off}^{R}} (\sqrt{η_{eq}} - 1)}, \end{aligned}

where we have used the definition of equilibrium fidelity (Equation S18). In traditional proofreading, a scheme with n proofreading realizations can yield a maximum fidelity of $η / η_{eq} = η_{eq}^{n}$ . The value of n for the original Hopfield model, for instance, is 1. It would be informative to also know the effective parameter n for the spatial proofreading model. Dividing Equation S21 by $η_{eq}$ , we find

\begin{aligned} \frac{η}{η_{eq}} = \frac{1}{\sqrt{η_{eq}}} e^{\sqrt{τ_{D} k_{off}^{R}} (\sqrt{η_{eq}} - 1)} & = η_{eq}^{n}, \\ e^{\sqrt{τ_{D} k_{off}^{R}} (\sqrt{η_{eq}} - 1)} & = η_{eq}^{n + \frac{1}{2}}, \\ \sqrt{τ_{D} k_{off}^{R}} (\sqrt{η_{eq}} - 1) & = (n + \frac{1}{2}) \ln η_{eq} \Rightarrow \\ n + \frac{1}{2} & = \frac{\sqrt{η_{eq}} - 1}{\ln η_{eq}} \sqrt{τ_{D} k_{off}^{R}} . \end{aligned}

This exact result can be simplified into an approximate form when diffusion is slow and $η_{eq} ≫ 1$ , yielding the expression reported in Equation 7 of the main text, namely,

\begin{array}{ll} n & \approx \frac{\sqrt{η_{eq}} \sqrt{τ_{D} k_{off}^{R}}}{\ln η_{eq}} \\ = \frac{\sqrt{τ_{D} k_{off}^{W}}}{\ln η_{eq}} . \end{array}

4. Fidelity in an intermediate substrate localization regime

The generic expression for complex density at the rightmost boundary ( $x = L$ ) can be written using Equation S13 as

ρ_{_{ES}} (L) = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S} (1 - λ_{_{ES}}^{2} / λ_{_{S}}^{2})} (\frac{λ_{_{ES}}}{λ_{_{S}} \sinh (L / λ_{_{ES}})} [\cosh (\frac{L}{λ_{_{ES}}}) e^{- L / λ_{_{S}}} - 1] + e^{- L / λ_{_{S}}}) .

For the system to proofread, substrates need to be sufficiently localized ( $λ_{_{S}} < L$ ) and diffusion needs to be sufficiently slow ( $τ_{D} k_{off}^{S} > 1$ or, $λ_{_{ES}} < L$ ). Under these conditions, the substrate profile can be approximated using Equation S3 as $ρ_{_{S}} (x) \approx λ_{_{S}}^{- 1} S_{total} e^{- x / λ_{_{S}}}$ , while the hyperbolic sine and cosine functions used above can be approximated as $\sinh (L / λ_{_{ES}}) \approx \cosh (L / λ_{_{ES}}) \approx 0.5 e^{L / λ_{_{ES}}}$ . With these approximations, the complex density expression simplifies into

\begin{aligned} ρ_{_{ES}} (L) & = \frac{k_{on} S_{total} ρ_{_{E}}}{k_{off}^{S} λ_{_{S}} (1 - λ_{_{ES}}^{2} / λ_{_{S}}^{2})} (\frac{λ_{_{ES}}}{λ_{_{S}}} [e^{- L / λ_{_{S}}} - 2 e^{- L / λ_{_{ES}}}] + e^{- L / λ_{_{S}}}) \\ = \frac{k_{on} S_{total} ρ_{_{E}}}{k_{off}^{S} (λ_{_{S}}^{2} - λ_{_{ES}}^{2})} ((λ_{_{S}} + λ_{_{ES}}) e^{- L / λ_{_{S}}} - 2 λ_{_{ES}} e^{- L / λ_{_{ES}}}) . \end{aligned}

Now, depending on how $λ_{_{S}}$ compares with $λ_{_{ES}}$ , there can be two qualitatively different regimes for the complex density, namely,

ρ_{_{ES}} (L) = ρ_{_{ES}}^{\infty} \times {\begin{matrix} \frac{2 L}{λ_{_{ES}}} e^{- L / λ_{_{ES}}}, & i f λ_{S} ≪ λ_{E S} (L / λ_{S} ≫ \sqrt{τ_{D} k_{o f f}^{S}}) \\ \frac{L}{λ_{S}} e^{- L / λ_{S}}, & i f λ_{E S} ≪ λ_{S} (\sqrt{τ_{D} k_{o f f}^{S}} ≫ L / λ_{S}) \end{matrix}

where we used the equilibrium complex density $ρ_{_{ES}}^{\infty}$ defined in Equation S14.

Notably, the first regime effectively corresponds to the case of ideal substrate localization where complex density is independent of the precise value of $λ_{_{S}}$ . The dimensionless number $\sqrt{τ_{D} k_{off}^{S}}$ sets the scale for the minimum $L / λ_{_{S}}$ value beyond which ideal localization can be assumed. Conversely, the second regime corresponds to the case where the distance traveled by a complex before dissociating is so short that the complex profile is dictated by the substrate profile itself. Because of that, the complex density reduction from its equilibrium limit is independent of the precise values of $τ_{D}$ and $k_{off}^{S}$ , as long as the condition $λ_{_{ES}} ≪ λ_{_{S}}$ is met.

The scheme yields its highest fidelity when both right and wrong complex densities are in the first regime (ideal localization). When both densities are in the second regime, fidelity is reduced down to its equilibrium value $η_{eq}$ (Appendix 1—table 1). The transition between these two extremes happens when the density profiles of right and wrong complexes fall under different regimes. Fidelity can be navigated in the transition zone by tuning the substrate gradient length scale $λ_{_{S}}$ . This is demonstrated in Appendix 1—figure 1 for three different choices of $η_{eq}$ . In all three cases, the dimensionless numbers $\sqrt{τ_{D} k_{off}^{R}}$ and $\sqrt{τ_{D} k_{off}^{W}}$ set the approximate range in which the bulk of fidelity enhancement occurs, as stated in the main text.

Appendix 1—table 1

Fidelity of the scheme in different regimes of right and wrong complex densities.

The upper-right cell is empty because the two conditions on $λ_{_{S}}$ cannot be simultaneously met, since $λ_{_{ER}} > λ_{_{EW}}$ by construction (follows from $k_{off}^{R} < k_{off}^{W}$ ).

	$λ_{_{S}} ≪ λ_{_{ER}}$	$λ_{_{S}} ≫ λ_{_{ER}}$
$λ_{_{S}} ≪ λ_{_{EW}}$	$\frac{λ_{_{EW}}}{λ_{_{ER}}} e^{L (λ_{_{EW}}^{- 1} - λ_{_{ER}}^{- 1})}$	-
$λ_{_{S}} ≫ λ_{_{EW}}$	$\frac{2 λ_{_{S}}}{λ_{_{ER}}} e^{L (λ_{_{S}}^{- 1} - λ_{_{ER}}^{- 1})}$	$η_{eq}$

Appendix 1—figure 1

Download asset Open asset

The effective number of proofreading realizations ( $n_{eff}$ ) as a function of $L / λ_{_{S}}$ .

The shaded region represents the range of $L / λ_{_{S}}$ values set by the key dimensionless numbers $\sqrt{τ_{D} k_{off}^{R}}$ and $\sqrt{τ_{D} k_{off}^{W}}$ . $τ_{D}$ values chosen for the demonstration were 60, 40, and 20 (in $1 / k_{off}^{R}$ units) for the three different choices of $η_{eq}$ , respectively.

5. Optimal diffusion time scale for maximum fidelity

Figure 3c of the main text illustrated the non-monotonic dependence of fidelity on the diffusion time scale $τ_{D}$ for different fixed values of $λ_{_{S}}$ . Here, we further explore this feature by asking what sets the optimal $τ_{D}$ . To gain analytical insights, we focus on the case where the system can proofread, which, as we argued in the previous section, happens when $λ_{_{S}}, λ_{_{ES}} < L$ . Under this condition, we identified two qualitatively different regimes of complex density reduction (Equation S26). Namely, we found that for sufficiently fast diffusion the system acted as if the substrates were localized ideally, whereas for sufficiently slow diffusion the complex density reduction was dictated solely by $λ_{_{S}}$ and did not discriminate between the two substrate kinds. These two limiting behaviors are indeed reflected in Figure 3c where in the low $τ_{D}$ limit (fast diffusion) the family of curves matches the dotted ideal localization curve, while in the high $τ_{D}$ limit (slow diffusion) all curves decay to 1, corresponding to the loss of error correction.

An intuitive approach for identifying the optimal $τ_{D}$ is to slow down diffusion up to the point where the density of wrong complexes at $x = L$ approaches a plateau and effectively stops decreasing. Going past this threshold would only reduce the density of right complexes at $x = L$ and thereby, reduce the fidelity. We know from Equation S26 that plateauing for wrong complexes happens when $λ_{_{EW}} ≪ λ_{_{S}}$ (equivalently, $\sqrt{τ_{D} k_{off}^{W}} ≫ L / λ_{_{S}}$ ). Hence, our first guess for the optimal diffusion time scale $τ_{D}^{*}$ is

\begin{aligned} (S27) & τ_{D}^{*} k_{off}^{W} & \sim {(\frac{L}{λ_{_{S}}})}^{2} \Rightarrow \\ (S28) & τ_{D}^{*} k_{off}^{R} & \sim \frac{k_{off}^{R}}{k_{off}^{W}} {(\frac{L}{λ_{_{S}}})}^{2} \Rightarrow \\ (S29) & τ_{D}^{*} / τ_{off}^{R} & \sim \frac{1}{η_{eq}} {(\frac{L}{λ_{_{S}}})}^{2} . \end{aligned}

To test the soundness of this expression, we compared its predictions to the optimal $τ_{D}$ values in Figure 3 that were identified numerically for different choices of $λ_{_{S}}$ . The results of the comparison are shown in Appendix 1—figure 2. As can be seen, for sufficiently high degrees of substrate localization ( $L / λ_{_{S}}$ ), the prediction of Equation S29 provides a good approximation of the true optimum. However, it is apparent that the prediction consistently underestimates the true $τ_{D}^{*}$ , which was expected since plateauing of $ρ_{_{EW}} (L)$ happens not under equality but a strict inequality condition (i.e. $\sqrt{τ_{D}^{*} k_{off}^{W}} ≫ L / λ_{_{S}}$ ). Because an exact analytical expression for $τ_{D}^{*}$ is not available, we performed different approximations to the fidelity formula and found an empirical correction term for our earlier estimate given by $2 (L / λ_{_{S}}) / \sqrt{η_{eq}}$ . The prediction for $τ_{D}^{*}$ with the correction term is now accurate starting a much lower value of $L / λ_{_{S}}$ , corresponding to a regime where the system proofreads once ( $n_{eff} \approx 1$ ). Overall, these analytical results provide good initial guesses for $τ_{D}^{*}$ which should be refined using a numerical approach for a higher accuracy.

Appendix 1—figure 2

Download asset Open asset

Optimal diffusion time scale for different choices of $λ_{_{S}}$ .

Blue dots represent the exact values obtained numerically for the data in Figure 3c. Dashed and solid lines represent the analytical estimates with and without the correction term. Vertical lines correspond to those values of $L / λ_{_{S}}$ that yield an integer number of effective proofreading realizations.

Appendix 2

Energetics of the scheme

We start this section by deriving an analytical expression for the minimum dissipated power, which was used in making Figure 4 of the main text. Then, we calculate the upper limit on fidelity enhancement available to our model for a finite substrate gradient length scale and compare this limit with the fundamental thermodynamic bound. We end the section by providing an estimate for the baseline cost of setting up gradients and compare this cost with the maintenance cost reported in the main text. Similar to our treatment of Appendix 1, here too our calculations are based on the $ρ_{_{E}} \approx constant$ assumption to allow for intuitive analytical results.

1. Derivation of the minimum dissipated power

As stated in the main text, we calculate the minimum rate of energy dissipation necessary for maintaining the substrate profiles as

\begin{matrix} P = \sum_{S=R,W} \int_{0}^{L} j_{_{S}} (x) μ (x) d x, \end{matrix}

where $j_{_{S}} (x) = k_{on} ρ_{_{S}} (x) ρ_{_{E}} - k_{off}^{S} ρ_{_{ES}} (x)$ is the net local substrate binding flux and $μ (x) = μ (0) + k_{B} T \ln ρ_{_{S}} (x) / ρ_{_{S}} (0) = μ (0) - k_{B} T \cdot x / λ_{_{S}}$ is the local chemical potential.

Our choice for the expression of power at steady state is motivated by that fact that the enzyme transport is passive and therefore, energy needs to be spent only on counteracting the local binding/unbinding events that tend to homogenize the substrate profile. To demonstrate the validity of our proposed expression more formally, we invoke the standard approaches for calculating power (Hill, 1977; Zhang et al., 2012). In particular, for a system that is described through discrete states with transition rates $k_{i \to j}$ between them, the rate of energy dissipation at steady state is given by

\begin{matrix} P = k_{B} T \sum_{i > j} (J_{i \to j} - J_{j \to i}) \ln \frac{k_{i \to j}}{k_{j \to i}}, \end{matrix}

where $J_{i \to j}$ is the flux from state i into state j. We note here that a similar expression for the rate of total entropy production involves a $\ln (J_{i \to j} / J_{j \to i})$ term (statistical forces) instead of the $\ln (k_{i \to j} / k_{j \to i})$ term (deterministic driving forces). At steady state, however, these two expressions are mathematically equivalent (Zhang et al., 2012). Our choice for Equation S31 stems from the better physical intuition that it provides in our context.

So far, the description of our system has been in terms of continuous density functions. To apply Equation S31 for calculating power, we consider the discrete-state representation of enzyme dynamics shown in Appendix 2—figure 1. There, space is discretized into intervals of size $δ x$ and diffusion is represented through jumps between neighboring sites with a rate $D / δ x^{2}$ . What keeps the system out of equilibrium is the spatially varying substrate profile $ρ_{_{S}} (x)$ .

Appendix 2—figure 1

Download asset Open asset

Discrete-state representation of diffusive transport and substrate binding/unbinding events.

Transparent clusters of different numbers of substrates illustrate the spatial variation of substrate concentration.

Because forward and backward diffusive transitions have identical rates, according to Equation S31 they will not contribute to energy dissipation (since $\ln (1) = 0$ ). The contribution from the remaining substrate binding/unbinding events can then be written as

\begin{matrix} P = k_{B} T \sum_{S=R,W} \sum_{i} (k_{on} ρ_{_{S}} (x_{i}) \times δ n_{i}^{E} - k_{off}^{S} \times δ n_{i}^{ES}) \ln \frac{k_{on} ρ_{_{S}} (x_{i})}{k_{off}^{S}}, \end{matrix}

where $δ n_{i}^{E} = ρ_{_{E}} δ x$ and $δ n_{i}^{ES} = ρ_{_{ES}} (x_{i}) δ x$ are the numbers of free and substrate–bound enzymes, respectively, in the $[x_{i}, x_{i} + δ x]$ interval. In the limit of a large number of discrete spatial intervals, the sum over i in Equation S32 can be rewritten as an integral over the coordinate x, namely,

\begin{matrix} P = k_{B} T \sum_{S=R,W} \int_{x = 0}^{\infty} \underset{j_{_{S}} (x)}{\underset{⏟}{(k_{on} ρ_{_{S}} (x) ρ_{_{E}} - k_{off}^{S} ρ_{_{ES}} (x))}} \ln \frac{k_{on} ρ_{_{S}} (x)}{k_{off}^{S}} d x . \end{matrix}

Comparing the form of Equation S33 to that of Equation S30 (with $μ (x)$ substituted), one can notice a difference in the terms that multiply $j_{_{S}} (x)$ . Specifically, in Equation S30 we have $μ (x) = μ (0) - k_{B} T \ln ρ_{_{S}} (0) + k_{B} T \ln ρ_{_{S}} (x)$ while the corresponding term in Equation S33 is $k_{B} T \ln (k_{on} / k_{off}^{S}) + k_{B} T \ln ρ_{_{S}} (x)$ . The difference between them, however, is in the parts that do not depend on x, while the spatially varying parts (namely, the $k_{B} T \ln ρ_{_{S}} (x)$ contributions) are identical. Now, since the number of bound complexes is constant at steady state, we have $\int_{0}^{\infty} j_{_{S}} (x) d x = 0$ . This means that the x-independent parts discussed earlier all integrate to zero, making the power estimates by Equation S30 and Equation S33 identical, thereby justifying our proposed expression.

To estimate power, we substitute the analytical expression for $ρ_{_{ES}} (x)$ found earlier (Equation S13) into $j_{_{S}} (x)$ and performing a somewhat cumbersome integral, obtain

\begin{matrix} β P = J_{bind} \sum_{S=R,W} \frac{1}{1 - λ_{_{S}}^{2} / λ_{_{ES}}^{2}} (\frac{λ_{_{ES}}}{λ_{_{S}}} \frac{\tanh (L / 2 λ_{_{ES}})}{\tanh (L / 2 λ_{_{S}})} - 1), \end{matrix}

where $β^{- 1} = k_{B} T$ , and $J_{bind} = k_{on} S_{total} ρ_{_{E}}$ is the net binding rate of each substrate. Figure 4 in the main text was made using this expression for power.

To get additional insights about this result, let us consider the case where substrates are highly localized ( $λ_{_{S}} ≪ L$ ) and diffusion is slow ( $λ_{_{ES}} ≪ L$ ) – conditions needed for effective proofreading. Under these conditions, the hyperbolic tangent terms become 1 and the expression for the power expenditure simplifies into

\begin{matrix} β P = J_{bind} \sum_{S=R,W} \frac{λ_{_{ES}}^{2}}{λ_{_{S}} (λ_{_{ES}} + λ_{_{S}})} . \end{matrix}

The monotonic increase of power with $λ_{_{ES}}$ suggests that energy is primarily spent on maintaining the concentration gradient of right substrates. This is not surprising, since typically right complexes travel a much greater distance into the low concentration region of the compartment before releasing the bound substrate (i.e. $λ_{_{ER}} ≫ λ_{_{EW}}$ ). Therefore, neglecting the contribution from wrong substrates and considering the range of $λ_{_{S}}$ values where the bulk of power–fidelity trade-off takes place ( $λ_{_{ER}} > λ_{_{S}} > λ_{_{EW}}$ ), we further simplify the power expression into

\begin{matrix} β P \approx \frac{J_{bind} λ_{_{ER}}}{λ_{_{S}}} = \frac{J_{bind} \cdot β Δ μ}{\sqrt{τ_{D} k_{off}^{R}}}, \end{matrix}

where we used the identities $β Δ μ = L / λ_{_{S}}$ and $λ_{_{ER}} = L / \sqrt{τ_{D} k_{off}^{R}}$ . This simple linear relation suggests that in order to maintain the exponential substrate profile, the minimum energy spent per substrate binding event should be at least $P / J_{bind} \approx k_{B} T \cdot λ_{_{ER}} / λ_{_{S}} > 1 k_{B} T$ (since $λ_{_{ER}} > λ_{_{S}}$ ).

We can also use Equation S36 to estimate the minimum dissipation per substrate binding event at $λ_{_{S}} \approx λ_{_{EW}}$ where the logarithmic power–fidelity scaling regime ends (see Figure 4 of the main text). Substituting the value of $λ_{_{S}}$ , we obtain $β P / J_{bind} \approx (λ_{_{ER}} / λ_{_{EW}}) = \sqrt{η_{eq}}$ , which is the result illustrated in Figure 4.

2. Limits on fidelity enhancement

The error reduction capacity of the spatial proofreading scheme improves with a greater difference in substrate off-rates, as was demonstrated in Figure 2 of the main text. At the same time, Figure 3c showed that the finite length scale of substrate localization (or, finite driving force) sets an upper limit on fidelity enhancement for substrates with fixed off-rates. It is therefore of interest to consider these two features together to find the absolute limit on fidelity enhancement available to our model and then compare it with the fundamental bound set by thermodynamics.

Intuitively, fidelity will be enhanced the most if the density of right complexes does not decay across the compartment, while that of wrong complexes decays maximally. The first condition can be met if diffusion is fast or if the unbinding rate of right substrates is low, in which case we have

\begin{matrix} ρ_{_{ER}} (L) \approx ρ_{_{ER}}^{\infty}, \end{matrix}

where $ρ_{E R}^{\infty}$ is the equilibrium density of right complexes. Conversely, when the unbinding rate of wrong substrates is very large, the density of wrong complexes is maximally reduced at the rightmost boundary and can be obtained from Equation S24 by taking the $λ_{_{ES}} \to 0$ limit, namely,

\begin{aligned} ρ_{_{EW}} (L) & \approx \frac{k_{on} ρ_{_{E}} ρ_{_{S}} (0) e^{- L / λ_{_{S}}}}{k_{off}^{W}} = \frac{k_{on} ρ_{_{E}} S_{total} e^{- L / λ_{_{S}}}}{λ_{_{S}} (1 - e^{- L / λ_{_{S}}}) k_{off}^{W}} \\ = \frac{k_{on} ρ_{_{E}} S_{total}}{k_{off}^{W} L} \times \frac{L e^{- L / λ_{_{S}}}}{λ_{_{S}} (1 - e^{- L / λ_{_{S}}})} \\ = ρ_{_{EW}}^{\infty} \times \frac{β Δ μ e^{- β Δ μ}}{1 - e^{- β Δ μ}} . \end{aligned}

Here, $ρ_{_{EW}}^{\infty}$ is the equilibrium density of wrong complexes, and $β Δ μ = L / λ_{_{S}}$ is the effective driving force of the scheme. Taking the ratio of Equations S37 and S38. Limits on fidelity enhancement, we obtain the largest fidelity enhancement of the scheme for the given driving force, namely,

η^{max} = \frac{ρ_{_{ER}} (L)}{ρ_{_{EW}} (L)} = \underset{η_{eq}}{\underset{⏟}{\frac{ρ_{_{ER}}^{\infty}}{ρ_{_{EW}}^{\infty}}}} \times \frac{e^{β Δ μ} - 1}{β Δ μ} \Rightarrow

{(η / η_{eq})}^{max} = (e^{β Δ μ} - 1) / β Δ μ .

When $β Δ μ ≳ 1$ (or, $λ_{_{S}} ≲ L)$ , the limit above gets further simplified into

\begin{matrix} {(η / η_{eq})}^{max} \approx e^{β Δ μ} / β Δ μ . \end{matrix}

Appendix 2—figure 2

Download asset Open asset

Fidelity enhancement as a function of the effective driving force for varying choices of $k_{off}^{W}$ .

The red dashed line indicates the thermodynamic bound given by $e^{β Δ μ}$ . The black dashed line corresponds to the model’s upper limit on fidelity enhancement given by Equation S40.

Now, thermodynamics imposes an upper bound on fidelity enhancement by any proofreading scheme operating with a finite chemical potential $Δ μ$ . This bound is equal to $e^{β Δ μ}$ and is reached when the entire chemical potential is used to increase the free energy difference between right and wrong substrates (Qian, 2006). Comparing it with the result in Equation S41, we can see that fidelity enhancement in the spatial proofreading model has the same exponential scaling term, but with an additional linear factor. Since the dominant contribution comes from the exponential term (as captured also in Appendix 2—figure 2), we can claim that our proposed model can operate very close to the fundamental thermodynamic limit.

3. Energetic cost to setup a concentration gradient

Earlier in the section, we calculated the rate at which energy needs to be dissipated to counteract the homogenizing effect that enzyme activity has on the substrate gradient. In addition to this cost, however, there is also a baseline cost for setting up a gradient in the absence of any enzyme. Here, we calculate this cost in the case where the gradient formation mechanism needs to work against diffusion that tends to flatten the substrate profile.

As before, we consider an exponentially decaying substrate gradient with a decay length scale $λ_{_{S}}$ and a total number of substrates $S_{total}$ . We write the minimum power $P_{D}$ required for counteracting the diffusion of substrates as

\begin{matrix} P_{D} = - \int_{0}^{L} J_{D} (x) μ^{'} (x) d x, \end{matrix}

where $J_{D} = - D_{_{S}} \nabla ρ_{_{S}} (x)$ is the diffusive flux, with $D_{_{S}}$ being the substrate diffusion constant. The rationale for writing this form is that diffusion moves substrates from a higher chemical potential region into a neighboring lower chemical potential region. The gradient maintaining mechanism would need to spend at least this chemical potential difference ( $δ μ = - μ^{'} (x) δ x$ ) per each substrate diffusing a distance $δ x$ down the chemical potential gradient. Adding up the contribution from all local neighborhoods with a local diffusive flux $J_{D} (x)$ results in Equation S42.

Now, substituting $ρ_{_{S}} (x) \sim e^{- x / λ_{_{S}}}$ for the substrate profile and $μ (x) = μ (0) + k_{B} T \ln (ρ_{_{S}} (x) / ρ_{_{S}} (0))$ for the chemical potential, we obtain

\begin{matrix} β P_{D} & = \int_{0}^{L} D_{_{S}} ρ_{_{S}}^{'} (x) {(\ln ρ_{_{S}} (x))}^{'} d x \\ = D_{_{S}} \int_{0}^{L} \frac{{(ρ_{_{S}}^{'} (x))}^{2}}{ρ_{_{S}} (x)} d x \\ = D_{_{S}} \int_{0}^{L} \frac{ρ_{_{S}} (x)}{λ_{_{S}}^{2}} d x \\ = \frac{D_{_{S}} S_{total}}{λ_{_{S}}^{2}}, \end{matrix}

where in the third step we used the relation $ρ_{_{S}}^{'} (x) = - ρ_{_{S}} (x) / λ_{_{S}}$ . This suggests that the minimum dissipated power required for setting up an exponential gradient increases quadratically with decreasing localization length scale $λ_{_{S}}$ .

It is informative to also make a comparison between this result and the earlier calculated minimum dissipation needed to counteract the enzyme’s homogenizing activity. Recall that when substrates were sufficiently localized and when diffusion was sufficiently slow, proofreading power could be approximated as (Equation S35)

\begin{matrix} β P \approx J_{bind} \frac{λ_{_{ES}}^{2}}{λ_{_{S}} (λ_{_{ES}} + λ_{_{S}})}, \end{matrix}

where $J_{bind} = k_{on} S_{total} ρ_{_{E}}$ is the total substrate binding flux. Using the identities $λ_{_{ES}} = \sqrt{D / k_{off}^{S}}$ and $K_{d}^{S} = k_{off}^{S} / k_{on}$ , we can calculate the ratio of the proofreading power to baseline power as

\begin{matrix} \frac{P}{P_{D}} & = \frac{k_{on} S_{total} ρ_{_{E}} λ_{_{ES}}^{2}}{D_{_{S}} S_{total}} \times \frac{λ_{_{S}}^{2}}{λ_{_{S}} (λ_{_{ES}} + λ_{_{S}})} \\ = \frac{D}{D_{_{S}}} \times \frac{ρ_{_{E}}}{K_{d}^{S}} \times \frac{λ_{_{S}} / λ_{_{ES}}}{1 + λ_{_{S}} / λ_{_{ES}}} . \end{matrix}

Presuming for simplicity that the enzyme and substrate diffusion constants are the same, we see that two factors determine the power ratio: (1) the amount of free enzyme in the system ( $ρ_{_{E}} / K_{d}^{S}$ ) and (2) the substrate localization length scale relative to the characteristic length scale of complex diffusion ( $λ_{_{S}} / λ_{_{ES}}$ ). Now, recall that the proofreading cost is spent largely on counteracting the homogenizing activity of the enzyme on right substrates (Appendix 2.1) and that the bulk of fidelity enhancement takes place when $λ_{_{S}} ≲ λ_{_{ER}}$ (Appendix 1.4). Therefore, when tuning $λ_{_{S}}$ down, initially the power ratio would only depend on the amount of free enzyme in the system ( $ρ_{_{E}} / K_{d}^{S}$ ) and then, with tighter substrate localization, the relative contribution of the proofreading power would start to decrease.

In the end, we would like to note that spatial gradients can also be set up using an external potential without a continuous dissipation of energy. In an in vivo setting, gravity can give rise to spatial structures in oocytes (Feric and Brangwynne, 2013), while in an in vitro setting, electric fields can create gradients and power the transport of the complex (Hansen et al., 2017). We leave the investigations of such alternative strategies to future work.

Appendix 3

Studies on the effect of catalysis on the model performance

In Appendix 1, we considered the rate of catalysis at the right boundary to be very small for the analytical simplicity of our derivations. This resulted in expressions for fidelity that were independent of the rate of catalysis r and allowed us to use the complex density at the right boundary as a proxy for speed. In this section, we relax this assumption and explore the consequences of having non-negligible catalysis rates on the model’s fidelity and on the speed–fidelity trade-off.

1. Derivation of the complex density profile $ρ_{_{ES}} (x)$

Accounting for catalysis in our model should be done through a boundary condition for the complex density equation (Equation S1). Earlier, we imposed a no-flux boundary condition at $x = L$ under the slow catalysis assumption. With non-negligible catalysis, this assumption is no longer valid, and the boundary condition is modified into

- D \frac{d ρ_{_{ES}}}{d x} |_{x = L} = \underset{catalysis flux}{\underset{⏟}{r ρ_{_{ES}} (L)}} .

Recall from Equations S4, S6 and S8 that the general solution for the complex profile had the form

\begin{aligned} (S47) & ρ_{E S} (x) & = C_{1} e^{- x / λ_{_{ES}}} + C_{2} e^{x / λ_{_{ES}}} + C_{p} e^{- x / λ_{_{S}}}, where \\ (S48) & C_{p} & = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S} (1 - \frac{λ_{_{ES}}^{2}}{λ_{_{S}}^{2}})} . \end{aligned}

Imposing the no-flux boundary condition at $x = 0$ allows us to eliminate one of the integration constants, namely,

\begin{aligned} (S49) & - D \frac{d ρ_{_{ES}}}{d x} |_{x = 0} & = - D (- \frac{C_{1}}{λ_{_{ES}}} + \frac{C_{2}}{λ_{_{ES}}} - \frac{C_{p}}{λ_{_{S}}}) = 0 \Rightarrow \\ (S50) & C_{2} & = C_{1} + \frac{λ_{_{ES}}}{λ_{_{S}}} C_{p} \Rightarrow \\ ρ_{_{ES}} (x) & = C_{1} (e^{- x / λ_{_{ES}}} + e^{x / λ_{_{ES}}}) + \frac{C_{p}}{λ_{_{S}}} (λ_{_{ES}} e^{x / λ_{_{ES}}} + λ_{_{S}} e^{- x / λ_{_{S}}}) \\ (S51) & = 2 C_{1} \cosh (x / λ_{_{ES}}) + \frac{C_{p}}{λ_{_{S}}} (λ_{_{ES}} e^{x / λ_{_{ES}}} + λ_{_{S}} e^{- x / λ_{_{S}}}) . \end{aligned}

Next, we impose the new boundary condition at $x = L$ (Equation S46), which yields

\begin{aligned} - D (\frac{2 C_{1}}{λ_{_{ES}}} \sinh (L / λ_{_{ES}}) + \frac{C_{p}}{λ_{_{S}}} (e^{L / λ_{_{ES}}} - e^{- L / λ_{_{S}}})) \\ = r (2 C_{1} \cosh (L / λ_{_{ES}}) + \frac{C_{p}}{λ_{_{S}}} (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}})) \\ \Rightarrow 2 C_{1} \sinh (L / λ_{_{ES}}) + \frac{C_{p}}{λ_{_{S}}} λ_{_{ES}} (e^{L / λ_{_{ES}}} - e^{- L / λ_{_{S}}}) \\ = - \underset{ε}{\underset{⏟}{\frac{λ_{_{ES}} r}{D}}} (2 C_{1} \cosh (L / λ_{_{ES}}) + \frac{C_{p}}{λ_{_{S}}} (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}})) . \end{aligned}

Note that we have introduced the dimensionless variable $ε$ , which, as will see later, will define the extent to which the presence of catalysis affects the fidelity. For convenience, here we write different equivalent forms for $ε$ as

\begin{matrix} ε = \frac{λ_{_{ES}} r}{D} = \frac{r}{\sqrt{D k_{off}^{S}}} = \frac{r}{L k_{off}^{S}} \sqrt{τ_{D} k_{off}^{S}} . \end{matrix}

Solving for the remaining unknown coefficient C₁ in Equation S52, we find

\begin{matrix} C_{1} = - \frac{C_{p}}{2 λ_{_{S}}} \frac{λ_{_{ES}} (e^{L / λ_{_{ES}}} - e^{- L / λ_{_{S}}}) + ε (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}})}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} . \end{matrix}

Lastly, we substitute this result for C₁ into Equation S51 and obtain a general expression for the complex density profile as

ρ_{E S} (x) = - \frac{C_{p}}{λ_{_{S}}} \frac{λ_{_{ES}} (e^{L / λ_{_{ES}}} - e^{- L / λ_{_{S}}}) + ε (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}})}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} \cosh (x / λ_{_{ES}}) + \frac{C_{p}}{λ_{_{S}}} (λ_{_{ES}} e^{x / λ_{_{ES}}} + λ_{_{S}} e^{- x / λ_{_{S}}}) .

One can show in a straightforward way that this result reduces to Equation S13 in the $ε \to 0$ limit.

2. Effects on fidelity in low and high substrate localization regimes

Accounting for the catalysis flux has made the general expression for the complex density profile even more incomprehensible. In order to gain insights about the qualitative as well as quantitative changes introduced by catalysis, we will focus on two characteristic limits of substrate localization – uniform substrate profile ( $λ_{_{S}} \to \infty$ ) and ideal substrate localization ( $λ_{_{S}} \to 0$ ).

2.1. Uniform substrate profile

In this case, no mechanism for localizing substrates is in play. Let us start off by evaluating the coefficient $C_{p}$ (Equation S48) in the $λ_{_{S}} \to \infty$ limit. Recalling from Equation S3 that $ρ_{_{S}} (0) = S_{total} / (λ_{_{S}} (1 - e^{- L / λ_{_{S}}})$ ), we find

\begin{aligned} (S56) & ρ_{_{S}} (0) & \approx \frac{S_{total}}{L} \Rightarrow \\ C_{p} & \approx \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{k_{off}^{S}} \\ \approx \frac{k_{on} S_{total} ρ_{_{E}}}{L k_{off}^{S}} \\ (S57) & = \frac{J_{bind}}{L k_{off}^{S}}, \end{aligned}

where $J_{bind} = k_{on} S_{total} ρ_{_{E}}$ is the total substrate binding flux.

Substituting the expression for $C_{p}$ into Equation S55 and eliminating all the terms that vanish upon taking the $λ_{_{S}} \to \infty$ limit, we obtain

\begin{aligned} ρ_{_{ES}} (x) & \approx C_{p} (1 - \frac{ε \cosh (x / λ_{_{ES}})}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})}) \\ = \frac{J_{bind}}{L k_{off}^{S}} \times \frac{\sinh (L / λ_{_{ES}}) + ε (\cosh (L / λ_{_{ES}}) - \cosh (x / λ_{_{ES}}))}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} . \end{aligned}

Ultimately, we are interested in knowing the rate of product formation defined via $v_{S} = r ρ_{_{ES}} (L)$ . We therefore evaluate the complex density at $x = L$ and multiply it by r, which yields

\begin{aligned} v_{S} = r ρ_{_{ES}} (L) & = J_{bind} \times (\frac{r}{L k_{off}^{S}}) \times \frac{\sinh (L / λ_{_{ES}})}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} \\ = J_{bind} \times (\frac{r}{L k_{off}^{S}}) \times \frac{\tanh (L / λ_{_{ES}})}{\tanh (L / λ_{_{ES}}) + ε} \\ \equiv J_{bind} \times (\frac{r}{L k_{off}^{S}}) \times \frac{\tanh (\sqrt{τ_{D} k_{off}^{S}})}{\tanh (\sqrt{τ_{D} k_{off}^{S}}) + ε}, \end{aligned}

where in the last step we wrote an equivalent expression using the $L / λ_{_{ES}} = \sqrt{τ_{D} k_{off}^{S}}$ identity. To analyze this result further, we will consider two limiting cases.

Case 1: Fast diffusion ( $\sqrt{τ_{D} k_{𝐨𝐟𝐟}^{𝐒}} ≪ 1$ ). If diffusion is fast, we can approximate the hyperbolic tangent functions as the arguments themselves (i.e. $\tanh (z) \approx z$ for $z ≪ 1$ ). Then, using the last form of $ε$ in Equation S53, we simplify the expression for speed as

\begin{aligned} v_{S} & \approx J_{bind} \times (\frac{r}{L k_{off}^{S}}) \times \frac{\sqrt{τ_{D} k_{off}^{S}}}{\sqrt{τ_{D} k_{off}^{S}} + \frac{r}{L k_{off}^{S}} \sqrt{τ_{D} k_{off}^{S}}} \\ (S60) & = J_{bind} \times (\frac{r}{L k_{off}^{S}}) \times \frac{1}{1 + \frac{r}{L k_{off}^{S}}} \Rightarrow \\ (S61) & v_{S} & = J_{bind} \times \frac{\tilde{r}}{k_{off}^{S} + \tilde{r}}, where \\ (S62) & \tilde{r} & = r / L . \end{aligned}

This is an intuitive result, suggesting that an enzyme that diffuses fast acts like a standard Michaelis–Menten enzyme with an effective catalysis rate $\tilde{r}$ . For such an enzyme, the probability of catalysis for a bound substrate is $\tilde{r} / (k_{off}^{S} + \tilde{r})$ . Multiplying this probability by the net substrate binding flux yields the expression for speed in Equation S61.

Fidelity of the model in this fast diffusion setting can be written as

η = \frac{v_{R}}{v_{W}} = \frac{k_{off}^{W} + \tilde{r}}{k_{off}^{R} + \tilde{r}} .

In the limit where catalysis is very slow ( $\tilde{r} ≪ k_{off}^{R}$ ), the equilibrium fidelity given by the ratio of off-rates is recovered. And in the opposite limit of very fast catalysis ( $\tilde{r} ≫ k_{off}^{W}$ ), the discriminatory capacity of the enzyme disappears altogether (Appendix 3—figure 1a).

Appendix 3—figure 1

Download asset Open asset

Dependence of fidelity on the catalysis rate in the case where the substrate profile is uniform.

(a) Fast diffusion setting ( $\sqrt{τ_{D} k_{off}^{R}} ≪ 1$ ). The highest fidelity reduction is a factor of $η_{eq}$ . (b) Slow diffusion setting ( $\sqrt{τ_{D} k_{off}^{R}} ≳ 1$ ). The highest fidelity reduction is a factor of $\sqrt{η_{eq}}$ . In both cases, $η_{e q} = 10$ was used.

Case 2: Slow diffusion ( $\sqrt{τ_{D} k_{𝐨𝐟𝐟}^{𝐒}} ≳ 1$ ). A more interesting case is when diffusion is slow. Now, the hyperbolic tangent functions in Equation S59 are approximately 1, allowing us to simplify the expression for speed into

\begin{array}{ll} v_{S} & = J_{bind} \times (\frac{r}{L k_{off}^{S}}) \times \frac{1}{1 + \frac{r}{L k_{off}^{S}} \sqrt{τ_{D} k_{off}^{S}}} \\ = J_{bind} \times \frac{\tilde{r}}{k_{off}^{S} + \tilde{r} \sqrt{τ_{D} k_{off}^{S}}} . \end{array}

Drawing an analogy between the above result and Equation S61, one can notice the presence of an extra $\sqrt{τ_{D} k_{off}^{S}}$ factor for $\tilde{r}$ in the denominator.

Evaluating the speeds of right and wrong product formation, we can write fidelity in this slow diffusion setting as

\begin{matrix} η = \frac{v_{R}}{v_{W}} = \frac{k_{off}^{W} + \tilde{r} \sqrt{τ_{D} k_{off}^{W}}}{k_{off}^{R} + \tilde{r} \sqrt{τ_{D} k_{off}^{R}}} . \end{matrix}

Like the fast diffusion case, when catalysis is very slow ( $\tilde{r} ≪ \sqrt{k_{off}^{R} / τ_{D}}$ or, equivalently, $r ≪ \sqrt{D k_{off}^{R}}$ ), the equilibrium fidelity is recovered. Unlike the fast diffusion case, however, if catalysis is very fast ( $r ≫ \sqrt{D k_{off}^{W}}$ ), the enzyme partly preserves its discriminatory capacity (Appendix 3—figure 1b). In this limit, a fidelity equal to the square root of the equilibrium fidelity is still attainable, namely,

\begin{matrix} η \approx \frac{\sqrt{k_{off}^{W}}}{\sqrt{k_{off}^{R}}} = \sqrt{η_{eq}} . \end{matrix}

This unexpected result suggests a potential advantage of localizing fast catalytic reactions instead of having them occur in a well–mixed solution.

2.2. Ideal substrate localization

We next consider the effect of catalysis on model fidelity in the ideal substrate localization limit ( $λ_{_{S}} \to 0$ ). We begin by evaluating the $C_{p} / λ_{_{S}}$ ratio that appears in the density profile expression (Equation S55). Using Equations S48 and Equations S3, we find

\begin{aligned} (S67) & ρ_{_{S}} (0) & \approx \frac{S_{total}}{λ_{_{S}}} \\ \frac{C_{p}}{λ_{_{S}}} & = \frac{k_{on} ρ_{_{S}} (0) ρ_{_{E}}}{λ_{_{S}} k_{off}^{S} (1 - λ_{_{ES}}^{2} / λ_{_{S}}^{2})} \\ \approx \frac{k_{on} S_{total} ρ_{_{E}}}{- k_{off}^{S} λ_{_{ES}}^{2}} \\ (S68) & = - \frac{J_{bind}}{D}, \end{aligned}

where in the last step we invoked the identities $λ_{_{ES}}^{2} = D / k_{off}^{S}$ and $J_{bind} = k_{on} S_{total} ρ_{_{E}}$ . We then substitute our result for $C_{p} / λ_{_{S}}$ into Equation S55 and simplify the complex density expression into

\begin{array}{ll} ρ_{_{ES}} (x) & = \frac{J_{bind}}{D} \times λ_{_{ES}} (\frac{e^{L / λ_{_{ES}}} + ε e^{L / λ_{_{ES}}}}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} \cosh (x / λ_{_{ES}}) - e^{x / λ_{_{ES}}}) \\ = J_{bind} \times \frac{λ_{_{ES}}}{D} \frac{\cosh ((L - x) / λ_{_{ES}}) + ε \sinh ((L - x) / λ_{_{ES}})}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} . \end{array}

To obtain the speed, we evaluate $ρ_{_{ES}} (x)$ at the right boundary ( $x = L$ ) and multiply it by r, namely,

\begin{array}{ll} v_{S} = r ρ_{_{ES}} (L) & = J_{bind} \underset{ε}{\underset{⏟}{\frac{λ_{_{ES}} r}{D}}} \frac{1}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} \\ = J_{bind} \times \frac{ε}{\sinh (\sqrt{τ_{D} k_{off}^{S}}) + ε \cosh (\sqrt{τ_{D} k_{off}^{S}})} . \end{array}

To evaluate the effect of catalysis further, we again consider two special limits – those of fast and slow diffusion.

Case 1: Fast diffusion ( $\sqrt{τ_{D} k_{𝐨𝐟𝐟}^{𝐒}} ≪ 1$ ). In this limit, the hyperbolic sine function can be approximated by its argument (i.e. $\sinh (z) \approx z$ for $z ≪ 1$ ), while the hyperbolic cosine function is approximately 1. Making these approximations and substituting the expression for $ε$ , we obtain

\begin{array}{ll} v_{S} & \approx J_{bind} \times \frac{\frac{r}{L k_{off}^{S}} \sqrt{τ_{D} k_{off}^{S}}}{\sqrt{τ_{D} k_{off}^{S}} + \frac{r}{L k_{off}^{S}} \sqrt{τ_{D} k_{off}^{S}}} \\ = J_{bind} \times \frac{\frac{r}{L k_{off}^{S}}}{1 + \frac{r}{L k_{off}^{S}}} \\ = J_{bind} \times \frac{\tilde{r}}{k_{off}^{S} + \tilde{r}} . \end{array}

This result is identical to what we found in the fast diffusion limit for the $λ_{_{S}} \to \infty$ setting (Equation S61), which is reasonable, since the location of substrate binding is irrelevant if diffusion is very fast (Appendix 3—figure 2a).

Appendix 3—figure 2

Download asset Open asset

Fidelity as a function of the catalysis rate in an ideal substrate localization setting.

(a) Fast diffusion case, where the behavior of the system is identical to that in Appendix 3—figure 1a. (b) Slow diffusion case where efficient proofreading is achieved. Catalysis can reduce the fidelity by up to a factor of $\sqrt{η_{eq}}$ . In both cases, $η_{e q} = 10$ was used.

Case 2: Slow diffusion ( $\sqrt{τ_{D} k_{𝐨𝐟𝐟}^{𝐒}} ≫ 1$ ). In this limit, the hyperbolic sine and cosine functions can be approximated as exponentials with a $1 / 2$ prefactor, simplifying the expression of speed into

\begin{matrix} v_{S} \approx J_{bind} \times \frac{2 ε}{1 + ε} e^{- \sqrt{τ_{D} k_{off}^{S}}} . \end{matrix}

Recalling the identity $ε = r / \sqrt{D k_{off}^{S}}$ (note that $ε$ depends on the substrate kind), we evaluate the speed for right and wrong product formation and, dividing them, obtain the fidelity as

\begin{matrix} η = \frac{v_{R}}{v_{W}} & = \frac{1 + r / \sqrt{D k_{off}^{W}}}{1 + r / \sqrt{D k_{off}^{R}}} \times \frac{\sqrt{k_{off}^{W}}}{\sqrt{k_{off}^{R}}} e^{\sqrt{τ_{D} k_{off}^{W}} - \sqrt{τ_{D} k_{off}^{R}}} \\ = \frac{1 + r / \sqrt{D k_{off}^{W}}}{1 + r / \sqrt{D k_{off}^{R}}} \times \sqrt{η_{eq}} e^{\sqrt{τ_{D} k_{off}^{R}} (\sqrt{η_{eq}} - 1)} . \end{matrix}

In the case where catalysis is slow ( $r ≪ \sqrt{D k_{off}^{R}}$ ), the first term in the fidelity expression becomes approximately 1, and the our earlier result obtained with no account of catalysis is recovered (Equation S21). In the opposite limit of fast catalysis ( $r ≫ \sqrt{D k_{off}^{W}}$ ), the first term is no longer 1, and we find

\begin{matrix} η & \approx \underset{1}{\underset{⏟}{\sqrt{\frac{k_{off}^{R}}{k_{off}^{W}}} \sqrt{η_{eq}}}} e^{\sqrt{τ_{D} k_{off}^{R}} (\sqrt{η_{eq}} - 1)} \\ = e^{\sqrt{τ_{D} k_{off}^{R}} (\sqrt{η_{eq}} - 1)} . \end{matrix}

As we can see, fast catalysis in the slow diffusion regime reduces the fidelity by $\sqrt{η_{eq}}$ or, equivalently, reduces the effective number of proofreading realizations by one half, without affecting the exponential amplification term (Appendix 3—figure 2b).

To conclude, our study demonstrated the expected reduction of fidelity with increasing catalysis rate. In the case of fast diffusion, up to a factor of $η_{eq}$ reduction is possible, as is the case for the original (Hopfield, 1974; Wong et al., 2018). In the case of slow diffusion, however, the cap on the amount of reduction is decreased down to $\sqrt{η_{eq}}$ . The advantage of this feature is most notable in the limit of a non-localized (i.e. uniform) substrate profile and fast catalysis where a diffusing enzyme is still capable of discriminating between substrates. This behavior would not be possible for a Michaelis–Menten enzyme in a well-mixed solution.

3. Effects on the speed–fidelity trade-off

In Figure 3a of the main text we explored the speed–fidelity trade-off in the slow catalysis limit. This trade-off arose in response to tuning the substrate localization length scale ( $λ_{_{S}}$ ) and the diffusion time scale ( $τ_{D}$ ). Here, we explore the changes to this trade-off behavior in the case where the effects of catalysis are not negligible. For concreteness, we focus on alterations to the Pareto front of the trade-off achieved in the $λ_{_{S}} \to 0$ limit.

Appendix 3—figure 3a compares the Pareto fronts in the cases of slow and fast catalysis limits. In each case, speed is normalized by the corresponding effective Michaelis–Menten speed that is reached in the fast diffusion limit and is given by $v_{_{MM}} = J_{bind} \times \tilde{r} / (k_{off}^{R} + \tilde{r})$ , where $\tilde{r} = r / L$ . One can notice a shift of the fast catalysis front toward the low-fidelity region, which was expected since earlier we observed the complete loss of substrate discrimination when diffusion and catalysis were both fast (Appendix 3—figure 2a).

Appendix 3—figure 3

Download asset Open asset

Pareto front of the speed–fidelity trade-off at different levels of catalytic activity.

(a) Cases of slow and fast catalysis limits, with the y-axis for speed normalized to the [0,1] interval. (b) Family of Pareto fronts for different choices of the catalysis rate. Speed on the y-axis is reported relative to the substrate binding flux $J_{bind}$ .

Appendix 3—figure 3a may leave an impression that faster catalysis leads to a less favorable speed–fidelity trade-off. Note, however, that the speed $v_{_{MM}} (\tilde{r})$ used to normalize the y-axis is itself a function of the catalysis rate and penalizes the fast catalysis case more than its slow counterpart. To eliminate this ambiguity, we plotted a family of Pareto fronts for increasing values of the catalysis rate but this time normalizing the y-axis by the r-independent quantity $J_{bind}$ (Appendix 3—figure 3b). As can be seen, faster catalysis in fact improves the speed–fidelity trade-off, meaning that in order to maximize fidelity at a given speed level, the best strategy would be to increase the catalysis rate and correspondingly slow down the diffusion.

A trade-off between speed and fidelity also arises in response to the sole alteration of the catalysis rate, while keeping the rest of the model parameters fixed. To explore this trade-off for an arbitrary fixed choice of $λ_{_{S}}$ and $τ_{D}$ , we begin by evaluating speed from Equation S55, namely,

\begin{array}{ll} v_{S} & = r ρ_{_{ES}} (L) \\ = r \times (- \frac{C_{p}}{λ_{_{S}}} \frac{λ_{_{ES}} (e^{L / λ_{_{ES}}} - e^{- L / λ_{_{S}}}) + ε (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}})}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} \cosh (L / λ_{_{ES}}) + \frac{C_{p}}{λ_{_{S}}} (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}})) \\ = r \times \frac{C_{p}}{λ_{_{S}}} (\frac{- λ_{_{ES}} (e^{L / λ_{_{ES}}} - e^{- L / λ_{_{S}}}) \times \cosh (L / λ_{_{ES}}) + (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}}) \times \sinh (L / λ_{_{ES}})}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} + \\ ε \times \underset{= 0}{\underset{⏟}{\frac{- (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}}) \times \cosh (L / λ_{_{ES}}) + (λ_{_{ES}} e^{L / λ_{_{ES}}} + λ_{_{S}} e^{- L / λ_{_{S}}}) \times \cosh (L / λ_{_{ES}})}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})}}}) \\ = r \times \frac{C_{p}}{λ_{_{S}}} \frac{(λ_{_{S}} \sinh (L / λ_{_{ES}}) + λ_{_{ES}} \cosh (L / λ_{_{ES}})) e^{- L / λ_{_{S}}} - λ_{_{ES}}}{\sinh (L / λ_{_{ES}}) + ε \cosh (L / λ_{_{ES}})} \\ = \frac{a_{_{S}} r}{1 + b_{_{S}} r} . \end{array}

In the last step, we introduced coefficients $a_{_{S}}$ and $b_{_{S}}$ that are independent from r, and used the fact that $ε \sim r$ .

Now, using the definition of fidelity and the result obtained above, we can write

\begin{aligned} η = \frac{v_{R}}{v_{W}} = \frac{a_{_{R}}}{a_{_{W}}} \frac{1 + b_{_{W}} r}{1 + b_{_{R}} r} . \end{aligned}

Notice that the ratio $a_{_{R}} / a_{_{W}} \equiv η_{0}$ is the fidelity in the limit of very slow catalysis ( $r \to 0$ ). Substituting it, we write

\begin{aligned} η & = η_{0} (\frac{1 + b_{_{R}} r - (b_{_{R}} - b_{_{W}}) r}{1 + b_{_{R}} r}) \\ (S77) & = η_{0} (1 - (b_{_{R}} - b_{_{W}}) \times \underset{v_{R} / a_{_{R}}}{\underset{⏟}{\frac{r}{1 + b_{_{R}} r}}}) \Rightarrow \\ (S78) & η & = η_{0} (1 - \frac{Δ b}{a_{_{R}}} v_{R}), \end{aligned}

where $Δ b = b_{_{R}} - b_{_{W}}$ . Recalling that $ε = λ_{_{ES}} r / D$ and noting the function form of the denominator in Equation S75, one can show that $b_{_{S}} = D^{- 1} λ_{_{ES}} / \tanh (L / λ_{_{ES}})$ . This is an increasing function of $λ_{_{ES}}$ and hence, a decreasing function of $k_{off}^{S}$ , implying that $Δ b > 0$ .

With this condition in mind, we can see from Equation S78 that speed and fidelity are anticorrelated with a linear slope when tuning the catalysis rate, unlike the more sophisticated trade-off relations when tuning the other model parameters. The peak fidelity $η_{0}$ is attained in the limit of vanishing speed. And conversely, speed is the highest when fidelity is the lowest for the given fixed values of $λ_{_{S}}$ and $τ_{D}$ (Appendix 3—figure 4).

Appendix 3—figure 4

Download asset Open asset

Linear trade-off between speed and fidelity when tuning the rate of catalysis.

$η_{min}$ is the fidelity in the fast catalysis limit and is up to $η_{eq}$ lower than $η_{0}$ (based on the results of the previous section). Linear scale is used for both axes.

Overall, our result illustrates the simple speed–fidelity trade-off that can be navigated by tuning the catalysis rate. This, for instance, can be achieved by changing the concentration of effectors that activate the enzyme for catalysis.

Appendix 4

Proofreading for substrates with different localization conditions

Following the original treatment by Hopfield, 1974, we have performed the studies of our model under the assumption that discrimination between right and wrong substrates is solely based on their off–rates ( $k_{off}^{W} > k_{off}^{R}$ ). Although this is often the signature difference between substrates, in a cellular setting substrate discrimination may occur through other factors also. For example, substrates may be present at different amounts or they may have non-identical on–rates. These differences, however, have a multiplicative effect on the fidelity (i.e. $η \sim (k_{on}^{R} [R]) / (k_{on}^{W} [W])$ ) and do not highlight the proofreading capacity of a particular model.

Unlike these two features, differences in the degree to which right and wrong substrates are localized can have a non-trivial effect on the proofreading performance. In this Appendix, we generalize our study of the model fidelity to cases where right and wrong substrates have unequal localization length scales $λ_{_{R}}$ and $λ_{_{W}}$ , respectively.

1. Limiting cases

We start off by exploring the limiting cases first. From the earlier derived Equation S14 and Equation S15, we know that the complex density at $x = L$ in very low ( $λ_{_{S}} ≫ L$ ) and very high ( $λ_{_{S}} ≪ L$ ) substrate localization regimes is given by

\begin{aligned} (S79) & ρ_{_{ES}}^{\infty} & = \frac{k_{on} {\bar{ρ}}_{_{S}} ρ_{_{E}}}{k_{off}^{S}} and \\ (S80) & ρ_{_{ES}}^{ideal} (L) & = ρ_{_{ES}}^{\infty} \times \frac{L / λ_{_{ES}}}{\sinh (L / λ_{_{ES}})}, \end{aligned}

respectively. Note that the complex density in the ideal localization case is necessarily lower than that in the case of a uniform profile, since the inequality $\sinh (L / λ_{_{ES}}) > L / λ_{_{ES}}$ holds for all choices of $λ_{_{ES}}$ . If $λ_{_{R}}$ and $λ_{_{W}}$ are not constrained to be equal, then the highest fidelity for a given $τ_{D}$ will be attained when the right substrates are distributed uniformly while the wrong substrates are highly localized ( $λ_{_{R}} ≫ L$ and $λ_{_{W}} ≪ L$ , respectively). We obtain the fidelity in this case as

\begin{aligned} η^{max} & = \frac{ρ_{_{ER}}^{\infty}}{ρ_{_{EW}}^{ideal} (L)} \\ = \frac{ρ_{_{ER}}^{\infty}}{ρ_{_{EW}}^{\infty}} \times \frac{\sinh (L / λ_{_{EW}})}{L / λ_{_{EW}}} \\ (S81) & = η_{eq} \times \frac{\sinh (L / λ_{_{EW}})}{L / λ_{_{EW}}} \Rightarrow \\ (S82) & \frac{η^{max}}{η_{eq}} & = \frac{\sinh (L / λ_{_{EW}})}{L / λ_{_{EW}}} \equiv \frac{\sinh (\sqrt{τ_{D} k_{off}^{W}})}{\sqrt{τ_{D} k_{off}^{W}}} . \end{aligned}

Notably, this result for maximum fidelity enhancement is independent of $k_{off}^{R}$ . Furthermore, it exceeds the ideal localization fidelity reported in the main text (Equation 5, derived in the $λ_{_{S}} \to 0$ limit), which was expected since now the right complexes on average travel a shorter distance to reach the activation site than the wrong complexes.

In the opposite scenario where the wrong substrates are uniformly distributed and the right ones are highly localized ( $λ_{_{R}} ≪ L$ and $λ_{_{W}} ≫ L$ , respectively), the system attains its lowest fidelity for a given $τ_{D}$ , namely,

\begin{aligned} η^{min} & = \frac{ρ_{_{ER}}^{ideal} (L)}{ρ_{_{EW}}^{\infty}} = \frac{ρ_{_{ER}}^{\infty}}{ρ_{_{EW}}^{\infty}} \times \frac{L / λ_{_{ER}}}{\sinh (L / λ_{_{ER}})} \\ (S83) & = η_{eq} \times \frac{L / λ_{_{ER}}}{\sinh (L / λ_{_{ER}})} \Rightarrow \\ (S84) & \frac{η^{min}}{η_{eq}} & = \frac{L / λ_{_{ER}}}{\sinh (L / λ_{_{ER}})} = \frac{\sqrt{τ_{D} k_{off}^{R}}}{\sinh (\sqrt{τ_{D} k_{off}^{R}})} . \end{aligned}

Since $L / λ_{_{ER}} < \sinh (L / λ_{_{ER}})$ , the lowest fidelity is less than the equilibrium fidelity itself ( $η^{min} < η_{eq}$ ), suggesting that the enzyme may in fact do anti-proofreading (Murugan et al., 2014) if the wrong substrates are generally closer to the catalytic site.

2. Intermediate levels of substrate localization

In Figure 3 inset as well as in Appendix 1.4, we explored the dependence of fidelity on the substrate localization length scale $λ_{_{S}}$ when it was the same for the two substrate kinds. Here, we expand this study to the case where this constraint is relaxed.

In particular, using Equation S24, we calculate complex densities and corresponding fidelity values as a function of $λ_{_{R}}$ for different fixed choices of the length scale ratio $λ_{_{R}} / λ_{_{W}}$ . The results of the study are captured in Appendix 4—figure 1. In the special case where the two length scales are equal ( $λ_{_{R}}$ = $λ_{_{W}}$ , solid black line), fidelity exhibits a monotonic depends on $L / λ_{_{R}}$ , and in the limit of ideal localization (very large $L / λ_{_{R}}$ ) the result in Equation 5 of the main text is recovered.

When $λ_{_{R}} \neq λ_{_{W}}$ , the dependence of fidelity on $L / λ_{_{R}}$ is no longer monotonic. If right substrates are more localized than the wrong ones (red curves), then the fidelity curves have a minimum where the enzyme does anti-proofreading (i.e. $η < η_{eq}$ ). The proofreading portion of the curves (when $η > η_{eq}$ ) is shifted to the right, suggesting that much higher substrate localization is needed for the enzyme to proofread.

The opposite case is when the right substrates have a shallower gradient than the wrong ones (blue curves). The fidelity curves are now shifted to the left and have a peak that is greater than the large $L / λ_{_{R}}$ limit of fidelity. This means that there is an optimal degree of substrate localization, going beyond which makes the model performance worse in terms of both error correction and energy consumption.

Appendix 4—figure 1

Download asset Open asset

Fidelity as a function of $L / λ_{_{R}}$ for different choices of the ratio $λ_{_{R}} / λ_{_{W}}$ .

The solid black line corresponds to the earlier studied regime where substrates had identical localization length scales. The blue curves represent the cases where $λ_{_{R}} > λ_{_{W}}$ , while the red curves represent the cases where $λ_{_{R}} < λ_{_{W}}$ . Numbers next to the curves correspond to the $λ_{_{R}} / λ_{_{W}}$ ratios used for generating them. Expressions for the highest and lowest fidelity values, as well as the fidelity expression in the limit where both substrates are highly localized are shown on the right side of the figure. $τ_{D} = 40 τ_{off}^{R}$ and $η_{eq} = 10$ were used for demonstration.

Over the course of its diffusive transport, a bound enzyme is more likely to deposit a right substrate in a substrate-depleted region than a wrong one, because right substrates stay attached to the enzyme for a longer time. Therefore, if the gradient-maintaining mechanism does not discriminate between substrates (which we assume is the case for the kinase/phosphatase-based one), then it will be easier for it to maintain the wrong ones localized since they tend to get deposited closer to the localization site (see Appendix 6—figure 1c as an example). This means that in a realistic setting the spatial organization of substrates is more likely to be in the advantageous blue region of Appendix 4—figure 1 where $λ_{_{R}} > λ_{_{W}}$ , facilitating the realization of spatial proofreading.

Appendix 5

Studies on the validity of the uniform free enzyme profile assumption

In our treatment of the model so far, we have assumed for mathematical convenience that free enzymes are in excess, which suggested the approximation $ρ_{_{E}} (x) \approx constant$ . Example enzyme density profiles shown in Appendix 5—figure 1, however, demonstrate that this assumption does not hold in general. Specifically, there is a depletion of free enzymes near the substrate localization site and abundance near the catalysis site. Because of this depletion at the leftmost edge, we expect a reduction in speed in comparison with our earlier treatment where a flat profile was assumed. In addition, if substrates have a weak gradient, we expect the fidelity to also be reduced, since more enzymes will bind substrates at intermediate positions, reducing the average travel distance to the catalytic site. In what follows, we discuss in greater detail the consequences of having a nonuniform free enzyme distribution on the model performance.

Appendix 5—figure 1

Download asset Open asset

Example profiles of free and substrate-bound enzymes.

Enzyme profiles are normalized so that the sum of areas under the curves is unity. The substrate profile (rescaled on the y-axis) is shown in transparent gray.

1. Effects that relaxing the $ρ_{_{E}} (x) \approx constant$ assumption has on the Pareto front

We begin by studying the effects of relaxing the uniform free enzyme profile assumption on the Pareto front of the speed–fidelity trade-off (Figure 3a of the main text). This front is reached in the ideal substrate localization limit ( $λ_{_{S}} \to 0$ ). Though in general enzyme profiles need to be obtained using numerical methods due to the nonlinearity of reaction–diffusion equations, in this particular limit ( $λ_{_{S}} \to 0$ ) an analytical solution is available. To obtain it, we write the reaction–diffusion equations in the bulk region of space as

\begin{aligned} (S85) & \frac{\partial ρ_{_{ER}}}{\partial t} & = D \frac{\partial^{2} ρ_{_{ER}}}{\partial x^{2}} - k_{off}^{R} ρ_{_{ER}} \\ (S86) & \frac{\partial ρ_{_{EW}}}{\partial t} & = D \frac{\partial^{2} ρ_{_{EW}}}{\partial x^{2}} - k_{off}^{W} ρ_{_{EW}} \\ (S87) & \frac{\partial ρ_{_{E}}}{\partial t} & = D \frac{\partial^{2} ρ_{_{E}}}{\partial x^{2}} + \sum_{S = R, W} k_{off}^{S} ρ_{_{ES}} . \end{aligned}

Substrate binding reactions did not enter the above equations, as they occur at the leftmost boundary only. They are instead accounted for via boundary conditions, which read

\begin{aligned} (S88) & - D \frac{\partial ρ_{E R}}{\partial x} |_{x = 0} & = k_{o n} S_{t o a l} ρ_{E} (0), \\ (S89) & - D \frac{\partial ρ_{E W}}{\partial x} |_{x = 0} & = k_{o n} S_{t o t a l} ρ_{E} (0), \\ (S90) & - D \frac{\partial ρ_{E}}{\partial x} |_{x = 0} & = - 2 k_{o n} S_{t o t a l} ρ_{E} (0), \end{aligned}

where $S_{total}$ is the total amount of free substrate of each kind concentrated at $x = 0$ .

Relating local enzyme concentrations

Considering the system at steady state, we add Equations S85-S87 and obtain

0 = D \frac{d^{2} ρ_{E R}}{d x^{2}} + D \frac{d^{2} ρ_{E W}}{d x^{2}} + D \frac{d^{2} ρ_{E}}{d x^{2}},

where we replaced the partial derivatives with total derivative since the profiles are time-independent. Dividing Equation S91 by D and integrating once, we find

\frac{d ρ_{_{ER}}}{d x} + \frac{d ρ_{_{EW}}}{d x} + \frac{d ρ_{_{E}}}{d x} = C_{1} .

The above relation must hold for arbitrary position x. Choosing $x = 0$ and noting that from Equations S88-S90 the sum of fluxes should be zero, we can claim that $C_{1} = 0$ . Integrating for the second time, we obtain

ρ_{E R} (x) + ρ_{E W} (x) + ρ_{E} (x) = C_{2},

where C₂ is now a different constant. To find it, we perform an integral for the last time across the entire compartment, namely,

\int_{0}^{L} (ρ_{E R} (x) + ρ_{E W} (x) + ρ_{E} (x)) d x = E_{t o t a l} = C_{2} L .

Here, we introduced the parameter $E_{total}$ as the total number of enzymes in the system (in free or bound forms). The constant C₂, which we will rename into $ρ_{0}$ , is then the average enzyme density, that is,

\begin{matrix} ρ_{0} = E_{total} / L . \end{matrix}

Substituting this result into Equation S93, we find an insightful relation between free and bound enzyme densities at an arbitrary position, namely,

ρ_{E} (x) = ρ_{0} - ρ_{E R} (x) - ρ_{E W} (x) .

This relation suggests that whenever the local concentration of bound enzymes is high, the local concentration of free enzymes should be correspondingly low, as we see reflected in the profiles of Appendix 5—figure 1.

Deriving the fidelity expression

Next, we consider Equations S85 and S86 separately at steady state, written in the form

D \frac{d^{2} ρ_{_{ES}}}{d x^{2}} - k_{off}^{S} ρ_{_{ES}} = 0.

The general solution to this ODE reads

ρ_{E S} (x) = C_{1}^{S} e^{- x / λ_{E S}} + C_{2}^{S} e^{x / λ_{E S}},

where $λ_{_{ES}} = \sqrt{D / k_{off}^{S}}$ , and $C_{1}^{S}$ and $C_{2}^{S}$ (S = R,W) are constants which are different for right and wrong complexes. The no-flux boundary condition at $x = L$ can be used to relate these constants and simplify the complex profile expression, namely,

\begin{aligned} (S99) & - D \frac{d ρ_{_{ES}} (x)}{d x} |_{x = L} & = - \frac{D}{λ_{_{ES}}} (- C_{1}^{S} e^{- L / λ_{_{ES}}} + C_{2}^{S} e^{L / λ_{_{ES}}}) = 0 \Rightarrow \\ (S100) & C_{2}^{S} & = e^{- 2 L / λ_{_{ES}}} C_{1}^{S} \Rightarrow \\ ρ_{_{ES}} (x) & = C_{1}^{S} e^{- x / λ_{_{ES}}} + C_{1}^{S} e^{- 2 L / λ_{_{ES}}} e^{x / λ_{_{ES}}} \\ = 2 C_{1}^{S} e^{- L / λ_{_{ES}}} \cosh (\frac{L - x}{λ_{_{ES}}}) \\ (S101) & = {\tilde{C}}_{1}^{S} \cosh (\frac{L - x}{λ_{_{ES}}}), \end{aligned}

where ${\tilde{C}}_{1}^{S} = 2 C_{1}^{S} e^{- L / λ_{_{ES}}}$ is a new constant coefficient introduced for convenience.

Now, the fidelity of the scheme is the ratio of right and wrong complex densities at $x = L$ . Using the result above, the fidelity can be written as

\begin{aligned} η = \frac{ρ_{_{ER}} (L)}{ρ_{_{EW}} (L)} = \frac{{\tilde{C}}_{1}^{R}}{{\tilde{C}}_{1}^{W}} . \end{aligned}

The ratio of these constant coefficients can be obtained by noting that the diffusive fluxes of right and wrong complexes at $x = 0$ are identical (from Equations S38 and S38), that is,

\begin{aligned} (S103) & - D \frac{\partial ρ_{_{ER}}}{\partial x} |_{x = 0} & = - D \frac{\partial ρ_{_{EW}}}{\partial x} |_{x = 0} \Rightarrow \\ (S104) & {\tilde{C}}_{1}^{R} \times \frac{\sinh (L / λ_{_{ER}})}{λ_{_{ER}}} & = {\tilde{C}}_{1}^{W} \times \frac{\sinh (L / λ_{_{EW}})}{λ_{_{EW}}} \Rightarrow \\ (S105) & \frac{{\tilde{C}}_{1}^{R}}{{\tilde{C}}_{1}^{W}} & = \frac{λ_{_{ER}}}{λ_{_{EW}}} \frac{\sinh (L / λ_{_{EW}})}{\sinh (L / λ_{_{ER}})} . \end{aligned}

Substituting this result into Equation S102, and recalling the equality $L / λ_{_{ES}} = \sqrt{τ_{D} k_{off}^{S}}$ , we obtain

\begin{matrix} η = \frac{\sqrt{τ_{D} k_{off}^{W}}}{\sqrt{τ_{D} k_{off}^{R}}} \frac{\sinh (\sqrt{τ_{D} k_{off}^{W}})}{\sinh (\sqrt{τ_{D} k_{off}^{R}})} = \sqrt{η_{eq}} \frac{\sinh (\sqrt{τ_{D} k_{off}^{W}})}{\sinh (\sqrt{τ_{D} k_{off}^{R}})} . \end{matrix}

This expression is identical to that in Equation S20 which was derived under the $ρ_{_{E}} (x) \approx constant$ assumption, suggesting that when substrates are highly localized, the shape of the free enzyme profile does not dictate the fidelity.

Deriving the speed expression

To keep the expression of speed compact while still illustrating the key consequences of relaxing the $ρ (x) \approx constant$ assumption, we will assume moving forward that the density of wrong complexes is much lower than that of the right complexes, that is, $ρ_{_{EW}} (x) ≪ ρ_{_{ER}} (x)$ . This assumption holds as long as the right and wrong complexes have sufficiently different off-rates. To see why it is the case, note that the ratio $ρ_{_{EW}} (x) / ρ_{_{ER}} (x)$ is the highest at $x = 0$ . We therefore calculate an upper bound for the ratio using Equation S101 and Equation S105 as

\begin{aligned} \frac{ρ_{E W} (x)}{ρ_{E R} (x)} < \frac{ρ_{E W} (0)}{ρ_{E R} (0)} = \frac{λ_{E W}}{λ_{E R}} \frac{\tanh (L / λ_{E R})}{\tanh (L / λ_{E W})} < \frac{λ_{E W}}{λ_{E R}} = \sqrt{\frac{k_{off}^{R}}{k_{off}^{W}}} = \frac{1}{\sqrt{η_{eq}}} . \end{aligned}

As long as $η_{eq} ≳ 10$ , it is fair to assume that the right complexes greatly outnumber the wrong ones, which allows us to approximate the free enzyme density from Equation S96 as $ρ_{_{E}} (x) \approx ρ_{0} - ρ_{_{ER}} (x)$ .

The specification of the right complex density profile requires the knowledge of the unknown coefficient ${\tilde{C}}_{1}^{R}$ . To find this coefficient, we use the boundary condition in Equation S88 and the approximation $ρ_{_{E}} (x) \approx ρ_{0} - ρ_{_{ER}} (x)$ to write

\begin{aligned} (S108) & D \frac{{\tilde{C}}_{1}^{R}}{λ_{_{ER}}} \sinh (L / λ_{_{ER}}) & = k_{on} S_{total} (ρ_{0} - {\tilde{C}}_{1}^{R} \cosh (L / λ_{_{ER}})) \Rightarrow \\ {\tilde{C}}_{1}^{R} & = \frac{k_{on} S_{total} ρ_{0}}{\frac{D}{λ_{_{ER}}} \sinh (L / λ_{_{ER}}) + k_{on} S_{total} \cosh (L / λ_{_{ER}})} \\ = \frac{k_{on} S_{total} ρ_{0}}{λ_{_{ER}} k_{off}^{R} \sinh (L / λ_{_{ER}}) + k_{on} S_{total} \cosh (L / λ_{_{ER}})} \\ (S109) & = ρ_{0} \times \frac{\frac{k_{on} S_{total}}{k_{off}^{R} L}}{1 + \frac{L}{λ_{_{ER}}} \frac{\cosh (L / λ_{_{ER}})}{\sinh (L / λ_{_{ER}})} \frac{k_{on} S_{total}}{k_{off}^{R} L}} \times \frac{L / λ_{_{ER}}}{\sinh (L / λ_{_{ER}})} . \end{aligned}

With the constant coefficient known, the right complex density then becomes

ρ_{_{ER}} (x) = ρ_{0} \times \frac{\frac{{\bar{ρ}}_{_{S}}}{K_{d}^{R}}}{1 + \frac{L}{λ_{_{ER}}} \frac{\cosh (L / λ_{_{ER}})}{\sinh (L / λ_{_{ER}})} \frac{{\bar{ρ}}_{_{S}}}{K_{d}^{R}}} \times \frac{L / λ_{_{ER}}}{\sinh (L / λ_{_{ER}})} \cosh (\frac{L - x}{λ_{_{ER}}}),

where we used the definitions of the mean substrate density ${\bar{ρ}}_{_{S}} = S_{total} / L$ and the dissociation constant $K_{d}^{R} = k_{off}^{R} / k_{on}$ .

To enable a direct parallel between this general treatment and the earlier one with the $ρ_{_{E}} (x) \approx constant$ approximation, let us introduce $ρ_{_{ER}}^{\infty}$ as the uniform right complex density when diffusion is very fast ( $λ_{_{ER}} ≫ L$ ) and calculate it from Equation S110 as

ρ_{_{ER}}^{\infty} = ρ_{0} \times \frac{\frac{{\bar{ρ}}_{_{S}}}{K_{d}^{R}}}{1 + \frac{{\bar{ρ}}_{_{S}}}{K_{d}^{R}}} .

Now, using the $ρ_{_{ER}}^{\infty}$ expression, we rewrite Equation S110 as

\begin{array}{ll} ρ_{_{ER}} (x) & = \frac{1 + \frac{{\bar{ρ}}_{_{S}}}{K_{d}^{R}}}{1 + \frac{L}{λ_{_{ER}}} \frac{\cosh (L / λ_{_{ER}})}{\sinh (L / λ_{_{ER}})} \frac{{\bar{ρ}}_{_{S}}}{K_{d}^{R}}} \times ρ_{_{ER}}^{\infty} \times \frac{L / λ_{_{ER}}}{\sinh (L / λ_{_{ER}})} \cosh (\frac{L - x}{λ_{_{ER}}}) \\ = \frac{1 + \frac{{\bar{ρ}}_{_{S}}}{K_{d}^{R}}}{1 + \underset{γ}{\underset{⏟}{\frac{L}{λ_{_{ER}}} \frac{\cosh (L / λ_{_{ER}})}{\sinh (L / λ_{_{ER}})}}} \frac{{\bar{ρ}}_{_{S}}}{K_{d}^{R}}} \times ρ_{_{ER}}^{const} (x), \end{array}

where $ρ_{_{ER}}^{const} (x)$ is the complex density obtained under the $ρ_{_{E}} (x) \approx constant$ assumption (Equation S15). The extra factor that appears on front does not exceed 1 since $γ \geq 1$ , indicating a reduction in speed, as we anticipated in our more qualitative discussion at the beginning of the section. The presence of the extra factor suggests two possibilities for the approximation to hold true; first, $γ \approx 1$ which happens when $λ_{_{ER}} ≳ L$ or when the right complex does not decay noticeably across the compartment, and second, when $γ > 1$ and ${\bar{ρ}}_{_{S}} ≪ γ^{- 1} K_{d}^{R}$ , which is when right complexes do decay but their fraction is low compared with free enzymes because of low substrate concentration.

Let us demonstrate the last statement more explicitly. Specifically, let us show that the validity of the approximation $ρ_{_{E}} (x) \approx constant$ is indeed linked directly to the fraction of bound enzymes. To that end, we evaluate $ρ_{_{E}} (0) / ρ_{_{E}} (L)$ as a metric that quantifies the degree to which $ρ_{_{E}} (x) \approx constant$ holds. If there is a large depletion of free enzymes near the substrate-binding site, then the metric will be significantly less than 1; conversely, if the free enzyme profile is practically flat, then the metric will be close to 1. Invoking the relation $ρ_{_{E}} (x) \approx ρ_{0} - ρ_{_{ER}} (x)$ and using our result for the complex density (Equation S110) as well as the definition of $γ$ in Equation S112, we evaluate this metric as

\begin{array}{ll} \frac{ρ_{_{E}} (0)}{ρ_{_{E}} (L)} & \approx \frac{ρ_{0} - ρ_{_{ER}} (0)}{ρ_{0} - ρ_{_{ER}} (L)} \\ = \frac{1 - \frac{γ {\bar{ρ}}_{_{S}} / K_{d}^{R}}{1 + γ {\bar{ρ}}_{_{S}} / K_{d}^{R}}}{1 - \frac{γ {\bar{ρ}}_{_{S}} / K_{d}^{R}}{\cosh (L / λ_{_{ER}}) (1 + γ {\bar{ρ}}_{_{S}} / K_{d}^{R})}} \\ = \frac{1}{1 + (1 - \frac{1}{\cosh (L / λ_{_{ER}})}) γ {\bar{ρ}}_{_{S}} / K_{d}^{R}} . \end{array}

Next, we calculate the fraction of bound enzymes $p_{bound}$ from Equation S110 as

\begin{array}{ll} p_{bound} & \approx E_{total}^{- 1} \int_{0}^{L} ρ_{_{ER}} (x) d x \\ = \frac{ρ_{0} L}{E_{total}} \frac{{\bar{ρ}}_{_{S}} / K_{d}^{R}}{1 + γ {\bar{ρ}}_{_{S}} / K_{d}^{R}} \\ = \frac{{\bar{ρ}}_{_{S}} / K_{d}^{R}}{1 + γ {\bar{ρ}}_{_{S}} / K_{d}^{R}} . \end{array}

Note that $γ^{- 1}$ emerges as the highest fraction of bound enzymes ( $p_{bound}^{max}$ ) reached in the large substrate concentration limit.

To link the metric $ρ_{_{E}} (0) / ρ_{_{E}} (L)$ to the fraction of bound enzymes, we express ${\bar{ρ}}_{_{S}} / K_{d}^{R}$ in terms of $p_{bound}$ and substitute it into Equation S113, namely,

\begin{aligned} (S115) & {\bar{ρ}}_{_{S}} / K_{d}^{R} & = \frac{p_{bound}}{1 - γ p_{bound}} \Rightarrow \\ \frac{ρ_{_{E}} (0)}{ρ_{_{E}} (L)} & = \frac{1}{1 + (1 - \frac{1}{\cosh (L / λ_{_{ER}})}) \frac{γ p_{bound}}{1 - γ p_{bound}}} \\ = \frac{1 - γ p_{bound}}{(1 - γ p_{bound}) + (1 - \frac{1}{\cosh (L / λ_{_{ER}})}) γ p_{bound}} \\ = \frac{1 - γ p_{bound}}{1 - γ p_{bound} / \cosh (L / λ_{_{ER}})} \\ (S116) & = \frac{p_{bound}^{max} - p_{bound}}{p_{bound}^{max} - p_{bound} / \cosh (L / λ_{_{ER}})} . \end{aligned}

Now, when the complexes do not decay appreciably across the compartment ( $λ_{_{ER}} ≳ L$ and thus, $\cosh (L / λ_{_{ER}}) \approx 1$ ), the metric becomes roughly equal to 1, suggesting that the free enzyme profile is practically flat. A more interesting case is when the complexes do decay ( $λ_{_{ER}} < L$ ), as in Appendix 5—figure 1. In this case, applying the condition $\cosh (L / λ_{_{ER}}) ≫ 1$ , we find

\frac{ρ_{E} (0)}{ρ_{E} (L)} \approx 1 - \frac{p_{b o u n d}}{p_{b o u n d}^{m a x}} .

The anti-correlation between the $ρ_{_{E}} (0) / ρ_{_{E}} (L)$ and $p_{bound}$ in the above result demonstrates that the degree to which the approximation $ρ_{_{E}} (x) \approx constant$ is violated is indeed dictated by the fraction of bound enzymes.

Pareto front shift

The previous calculations showed that in the ideal substrate localization limit relaxing the $ρ (x) \approx constant$ assumption keeps the fidelity the same while the speed gets reduced. And this reduction is greater for higher substrate concentrations. We therefore expect a shift in the Pareto front when going to the high substrate concentration limit, as is illustrated in Appendix 5—figure 2a. To get more intuition about the effect of this shift caused by tuning the amount of substrates, we consider the effective number of proofreading realizations at half-maximum speed (n₅₀) and study how this number changes as a function of the fraction of enzymes bound ( $p_{bound}$ ), which increases monotonically with $S_{total}$ as suggested by Equation S114. Appendix 5—figure 2b shows this dependence. As can be seen, n₅₀ reduces roughly linearly with $p_{bound}$ ; for example, if 10% of the enzymes are bound, then a 10% reduction in n₅₀ is expected. This suggests that as long as the fraction of bound enzymes is low, our findings related to the Pareto front made under the $ρ_{_{E}} \approx constant$ assumption will generally hold true.

Appendix 5—figure 2

Download asset Open asset

Consequences of relaxing the $ρ_{_{E}} (x) \approx constant$ assumption on the Pareto front.

(a) Pareto fronts in the low and high substrate concentration limits. (b) Reduction in the effective number of proofreading realizations at half-maximum speed as a function of the fraction of enzymes bound. $η_{eq} = 10$ was used in making the plots.

2. Effects that relaxing the $ρ_{_{E}} (x) \approx constant$ assumption has on fidelity in a weak substrate gradient setting

In this section, we study how accounting for the spatial distribution of free enzymes affects our results on the model’s fidelity in the setting where substrates have a finite localization length scale $λ_{_{S}}$ . In this setting, Equations (1–3) (in the main text) describing the system’s dynamics become a system of nonlinear equations, which we solve at steady state using numerical methods.

An example curve of how fidelity changes with tuning diffusion time scale in a finite $λ_{_{S}}$ setting is shown in Appendix 5—figure 3. As expected, the nonuniform free enzyme profile leads to a reduction in fidelity. This reduction is not significant when diffusion is relatively fast as in that case the free enzyme profile manages to flatten out rapidly. The reduction is not significant also in the very slow diffusion limit where binding events that lead to production primarily take place in the proximity of the activation region and hence, the nonuniform profile of free enzymes across the compartment has little impact on fidelity. The greatest reduction happens at intermediate diffusion time scales; in particular, when the system achieves its peak fidelity.

Appendix 5—figure 3

Download asset Open asset

Fidelity as a function of diffusion time scale calculated with and without making the $ρ_{_{E}} (x) \approx constant$ approximation.

The total number of free substrates is chosen so that ${\bar{ρ}}_{_{S}} / K_{d}^{R} = 3$ . The substrate localization length scale used for generating the solid curves is $λ_{_{S}} / L = 0.04$ .

To quantify the extent of this highest reduction, we calculated the peak value of the effective number of proofreading realizations ( $n_{max}$ ) for different free substrate amounts which regulate the fraction of bound enzymes ( $p_{bound}$ ). The results obtained for different choices of $λ_{_{S}}$ are summarized in Appendix 5—figure 4. As can be seen, for the high substrate localization case ( $λ_{_{S}} / L = 0.04$ ), there is a roughly linear dependence between $n_{max}$ and $p_{bound}$ . The initial decrease in $n_{max}$ with growing $p_{bound}$ is even slower when substrates are less tightly localized ( $λ_{_{S}} / L = 0.10, 0.30$ ).

Appendix 5—figure 4

Download asset Open asset

Reduction in the peak effective number of proofreading realizations as a function of $p_{bound}$ .

$n_{m a x}^{l o w}$ represents the peak value of $n_{eff}$ in the limit of low substrate concentration (the maximum of the solid blue curve in Appendix 5—figure 3).

Taken together, these results suggest that if the substrate concentration is low enough to leave most of the enzymes unbound, then our proposed scheme will proofread efficiently. And this requirement on substrate amount will be further relaxed if diffusion is fast, or if substrates are not very tightly localized.

Appendix 6

Proofreading on a kinase/phosphatase-induced gradient

In this section, we introduce the mathematical modeling setup for the kinase/phosphatase-based gradient formation scheme and describe how its fidelity is calculated numerically. In the end, we discuss the energetics of setting up the substrate concentration gradient and link our calculations to the lower bounds on energy cost obtained earlier in Appendix 2.

1. Setup and estimation of fidelity

In the analysis thus far, we have imposed a gradient of free substrates and analyzed the proofreading capability of an enzyme acting on this gradient. In a living cell, gradients themselves are maintained by active cellular processes. However, the action of the enzyme – that is, binding a substrate in one spatial location, diffusing away, and releasing the substrate elsewhere – can destroy the gradient, and thereby lead to a loss of proofreading. Here, we analyze the consequences of free substrate depletion and gradient flattening caused by the enzyme.

We model the formation of a substrate gradient by a combination of localized activation and delocalized deactivation. We suppose that substrates can exist in phosphorylated or dephosphorylated forms, and that only the phosphorylated form is capable of binding to the enzyme. The substrates are phosphorylated by a kinase with rate $k_{kin} = 0.2$ s⁻¹, and dephosphorylated by a phosphatase with rate $k_{p} = 5$ s⁻¹. Crucially, we assume that phosphatases are found everywhere in the domain of size $L \sim 10$ μm (a typical length scale in a eukaryotic cell), while kinases are localized to one end of the domain (at $x = 0$ ), as may occur naturally if kinases are bound to one of the membranes enclosing the domain.

The minimal dynamics of phosphorylated substrates and enzyme–substrate complexes is then given by

\begin{array}{ll} \frac{\partial ρ_{S}}{\partial t} & = D \nabla^{2} ρ_{S} - k_{b} ρ_{S} + k_{o f f}^{S} ρ_{E S} - k_{p} ρ_{S}, \\ \frac{\partial ρ_{E S}}{\partial t} & = D \nabla^{2} ρ_{E S} + k_{b} ρ_{S} - k_{o f f}^{S} ρ_{E S}, \end{array}

augmented by the boundary conditions

\begin{array}{ll} Substrate phosphorylation: - D \nabla ρ_{_{S}} |_{x = 0} = k_{kin}, \\ No-flux: - D \nabla ρ_{S} |_{x = L} = - D \nabla ρ_{E S} |_{x = L} = - D \nabla ρ_{E S} |_{x = 0} = 0. \end{array}

Here, we have supposed that the densities of free enzymes, dephosphorylated substrates, and phosphatases are fixed and uniform, and have absorbed them into the relevant rate constants ( $k_{b} = k_{on} ρ_{_{E}}$ , $k_{kin}$ , and $k_{p}$ , respectively). For simplicity, we have also assumed that the free substrates and enzyme–substrate complexes have the same diffusion coefficient $D = 1$ μm²/s. We note that accounting for distinct diffusivities of phosphorylated and unphosphorylated substrate forms (Kholodenko, 2009) would affect the speed, while accounting for the slower diffusion of the enzyme–substrate complex would alter the estimates of both speed and fidelity of the model. One or several of these effects can be considered when studying a specific biological system where these microscopic details are known.

We numerically solve Equations S118 and S119 at steady state to obtain the concentration profiles. First, the equations of dynamics are made dimensionless by settings units of length and time by L ( $\bar{x} = x / L$ ) and $τ_{D} \equiv L^{2} / D$ ( $\bar{t} = t / τ_{D}$ ), respectively. At steady state, the dimensionless equations read

\begin{aligned} {\bar{\nabla}}^{2} {\bar{ρ}}_{S} & = ({\bar{k}}_{b} + {\bar{k}}_{p}) {\bar{ρ}}_{S} - {\bar{k}}_{o f f}^{S} {\bar{ρ}}_{E S}, \\ {\bar{\nabla}}^{2} {\bar{ρ}}_{E S} & = - {\bar{k}}_{b} {\bar{ρ}}_{S} + k_{o f f}^{S} {\bar{ρ}}_{E S}, \end{aligned}

with boundary conditions

\begin{aligned} \bar{\nabla} {\bar{ρ}}_{S} |_{\bar{x} = 0} & = - {\bar{k}}_{kin}, \\ \bar{\nabla} {\bar{ρ}}_{S} |_{\bar{x} = 1} & = \bar{\nabla} {\bar{ρ}}_{E S} |_{\bar{x} = 1} = \bar{\nabla} {\bar{ρ}}_{E S} |_{\bar{x} = 0} = 0, \end{aligned}

where concentrations have been rescaled as $\bar{ρ} = ρ L$ , and kinetic rates as $\bar{k} = k τ_{D}$ .

We discretize the steady state equations on a grid with spacing $Δ \bar{x} = 0.01$ , approximating the second derivative as

{\bar{\nabla}}^{2} \bar{ρ} \approx \frac{1}{Δ {\bar{x}}^{2}} (\bar{ρ} (\bar{x} + Δ \bar{x}) + \bar{ρ} (\bar{x} - Δ \bar{x}) - 2 \bar{ρ} (\bar{x})) .

This is ill-defined at the boundaries $\bar{x} = 0$ and $\bar{x} = 1$ , which is addressed by incorporating the boundary conditions. For illustration, consider the left boundary, $\bar{x} = 0$ , and suppose that our domain included also a point at $\bar{x} = - Δ \bar{x}$ . Then, we could approximate the boundary condition ${\bar{\nabla} {\bar{ρ}}_{_{S}} |}_{\bar{x} = 0} = - {\bar{k}}_{kin}$ by a centred difference scheme, and solve out for the fictional point at $\bar{x} = - Δ \bar{x}$ , namely,

\begin{array}{ll} \bar{\nabla} {\bar{ρ}}_{S} |_{\bar{x} = 0} = - {\bar{k}}_{kin} \\ \Rightarrow \frac{1}{2 Δ \bar{x}} ({\bar{ρ}}_{S} (Δ \bar{x}) - {\bar{ρ}}_{S} (- Δ \bar{x})) = - {\bar{k}}_{kin} \\ \Rightarrow {\bar{ρ}}_{S} (- Δ \bar{x}) = {\bar{ρ}}_{S} (Δ \bar{x}) + 2 Δ \bar{x} {\bar{k}}_{kin}, \end{array}

which, when inserted into Equation S122, specifies ${\bar{\nabla}}^{2} {\bar{ρ}}_{_{S}}$ at $\bar{x} = 0$ , that is,

{\bar{\nabla}}^{2} {\bar{ρ}}_{_{S}} |_{\bar{x} = 0} = \frac{1}{Δ {\bar{x}}^{2}} (2 {\bar{ρ}}_{_{S}} (Δ \bar{x}) - 2 {\bar{ρ}}_{_{S}} (0)) + \frac{2}{Δ \bar{x}} {\bar{k}}_{kin} .

For the boundary at the right ( $\bar{x} = 1$ ) as well as for the boundary conditions for ${\bar{ρ}}_{_{ES}}$ , we similarly implement no-flux boundary conditions. After discretizing, Equation S120 can then be written in a matrix form as

\begin{array}{ll} \overset{M_{S}}{\overset{⏞}{(\frac{1}{Δ {\bar{x}}^{2}} (\begin{matrix} - 2 & 2 & 0 & \dots & 0 \\ 1 & - 2 & 1 & \dots & 0 \\ ⋮ & ⋮ & ⋮ & ⋱ & ⋮ \\ 0 & \dots & 1 & - 2 & 1 \\ 0 & 0 & \dots & 1 & - 1 \end{matrix}) - ({\bar{k}}_{b} + {\bar{k}}_{p}) I)}} {\vec{ρ}}_{S} = - {\bar{k}}_{o f f}^{S} {\vec{ρ}}_{ES} + \overset{\vec{b}}{\overset{⏞}{(\begin{matrix} - \frac{2}{Δ \bar{x}} {\bar{k}}_{kin} \\ 0 \\ ⋮ \\ 0 \\ 0 \end{matrix})}}, \\ \underset{M_{ES}}{\underset{⏟}{(\frac{1}{Δ {\bar{x}}^{2}} (\begin{matrix} - 1 & 1 & 0 & \dots & 0 \\ 1 & - 2 & 1 & \dots & 0 \\ ⋮ & ⋮ & ⋮ & ⋱ & ⋮ \\ 0 & \dots & 1 & - 2 & 1 \\ 0 & 0 & \dots & 1 & - 1 \end{matrix}) - {\bar{k}}_{o f f}^{S} I)}} {\vec{ρ}}_{ES} = - {\bar{k}}_{b} {\vec{ρ}}_{S}, \end{array}

where ${\vec{ρ}}_{S}$ , ${\vec{ρ}}_{ES}$ are column vectors of the nondimensionalized concentration profiles evaluated at the spatial grid points, that is, ${[\bar{ρ} (0), \bar{ρ} (Δ \bar{x}), \dots]}^{T}$ . Solving these matrix equations yields

\begin{aligned} {\vec{ρ}}_{S} & = {(M_{S} - {\bar{k}}_{o f f}^{S} {\bar{k}}_{b} M_{ES}^{- 1})}^{- 1} \vec{b}, \\ {\vec{ρ}}_{E S} & = - {\bar{k}}_{b} {(M_{S} M_{ES} - {\bar{k}}_{o f f}^{S} {\bar{k}}_{b} I)}^{- 1} \vec{b} . \end{aligned}

We compute Equation S125 numerically for two substrates: a cognate (‘R’) and a non-cognate (‘W’), which differ in their off-rates ( $k_{off}^{R} = 0.1 s^{- 1}$ and $k_{off}^{W} = 1 s^{- 1}$ , respectively). Having the density profiles, the fidelity of the model becomes $η \approx {\bar{ρ}}_{_{ER}} (\bar{x} = 1) / {\bar{ρ}}_{_{EW}} (\bar{x} = 1)$ . We calculate the fidelity for different choices of the first–order rate of enzyme–substrate binding ( $k_{b} = k_{on} ρ_{_{E}}$ ); this may be thought of as varying the concentration of free enzyme in the cell. The results are shown in Figure 5 of the main text.

2. Energy dissipation

In Appendices 2.1 and 2.3, we estimated lower bounds on the minimum power that needs to be dissipated in order to counter the homogenizing effect that enzyme activity and substrate diffusion respectively have on localized substrate profiles. Here, we calculate the energy dissipation required to run the kinase/phosphatase-based mechanism and compare it with these lower bounds estimated earlier.

Let us assume that phosphorylation and dephosphorylation reactions by kinases and phosphatases are nearly irreversible with associated free energy costs of $Δ ε_{kin}$ and $Δ ε_{phosph}$ per reaction, respectively. The net rate at which active substrates get dephosphorylated is $k_{p} S_{phosphorylated}$ and it needs to be identical to the net phosphorylation rate of inactive substrates in order for $S_{phosphorylated}$ to remain constant. With the costs of each reaction known, we can write the rate of energy dissipation $P_{k/p}$ as

P_{k/p} = k_{p} S_{phosphorylated} (Δ ε_{kin} + Δ ε_{phosph}) .

To gain analytical intuition, we first consider the case where the enzyme activity is very low, so that the kinase/phosphatase–based mechanism maintains an exponential profile of active substrates with a decay length scale $λ_{_{S}} = \sqrt{D_{_{S}} / k_{p}}$ . Expressing the rate of phosphorylation in terms of $λ_{_{S}}$ and $D_{_{S}}$ (i.e., $k_{p} = D_{_{S}} / λ_{_{S}}^{2}$ ), and substituting it into Equation S126, we obtain

\begin{aligned} P_{k/p} = \frac{D_{_{S}} S_{phosphorylated}}{λ_{_{S}}^{2}} (Δ ε_{kin} + Δ ε_{phosph}) . \end{aligned}

Comparing this result with the lower dissipation bound found earlier (Equation S43), we can note the presence of an extra factor $β (Δ ε_{kin} + Δ ε_{phosph})$ . Since the free energy consumption during ATP hydrolysis is $\sim 10 k_{B} T$ , we can say that the power dissipated by the kinase/phosphatase system for setting up an exponential gradient surpasses the lower limit necessary for counteracting diffusion roughly by an order of magnitude.

Next, we explore the energetics of the kinase/phosphatase-based mechanism in the context of the power–fidelity trade-off. Our study of the trade-off in Figure 4 of the main text was performed under the assumption that substrate profiles were exponentially decaying in the entire spatial domain. In Appendix 6—figure 1a, we show the trade-off curves obtained under this assumption and compare them with the trade-off curve for the kinase/phosphatase-based mechanism that arises in response to changing the substrate localization by tuning $k_{p}$ . As can be seen, the predicted lower bound (sum of the minimum powers needed to counteract the enzyme action and substrate diffusion) is roughly an order of magnitude lower than the total dissipation of the mechanism, and this difference increases with higher fidelity.

Note, however, that the assumption about an exponential substrate localization is not generally valid for the kinase/phosphatase-based mechanism because substrates can be deposited in low–concentration regions and not get immediately dephosphorylated (Appendix 6—figure 1c). We therefore refine our lower bounds on the dissipated power by estimating them numerically using their generic definitions, namely, Equation S30 for counteracting enzymatic action, and Equation S42 for counteracting substrate diffusion. These refined estimates suggest a factor of ∼10 difference between the total cost and its lower bound consistently across a wide region of the trade-off curve. This means that substrate gradient maintenance through practically irreversible phosphorylation and dephosphorylation reactions has low energetic efficiency for doing spatial proofreading, which, however, may be sustainable depending on the energy budget of the cell.

Appendix 6—figure 1

Download asset Open asset

Energetic performance of the kinase/phosphatase–based mechanism.

(a) Total dissipation and calculated lower bounds under the assumption of exponential substrate localization. (b) Total dissipation and lower bounds estimated without assuming exponential substrate profiles. In both (a) and (b) $k_{b} = 1 s^{- 1}$ and $Δ ε_{kin} = Δ ε_{phosph} = 10 k_{B} T$ were used. (c) Example profiles of right and wrong substrates for the physiologically relevant dephosphorylation rate $k_{p} = 5 s^{- 1}$ . Exponential decay of the substrate profile with the predicted length scale $λ_{_{S}} = \sqrt{D_{_{S}} / k_{p}}$ holds in the first $\sim 2$ μm of the compartment.

Data availability

All scripts used to generate the data for making the plots are provided in supporting files.

References

(2002) Secretory pathway quality control operating in Golgi, plasmalemmal, and endosomal systems
Traffic 3:771–780.
https://doi.org/10.1034/j.1600-0854.2002.31102.x
- PubMed
- Google Scholar
(2006) Gradients in the self-organization of the mitotic spindle
Trends in Cell Biology 16:125–134.
https://doi.org/10.1016/j.tcb.2006.01.005
- PubMed
- Google Scholar
1. Bauman AL
2. Scott JD
(2002) Kinase- and phosphatase-anchoring proteins: harnessing the dynamic duo
Nature Cell Biology 4:E203–E206.
https://doi.org/10.1038/ncb0802-e203
- PubMed
- Google Scholar
1. Bivona TG
2. Pérez De Castro I
3. Ahearn IM
4. Grana TM
5. Chiu VK
6. Lockyer PJ
7. Cullen PJ
8. Pellicer A
9. Cox AD
10. Philips MR
(2003) Phospholipase Cγ activates Ras on the Golgi apparatus by means of RasGRP1
Nature 424:694–698.
https://doi.org/10.1038/nature01806
- PubMed
- Google Scholar
1. Brown GC
2. Kholodenko BN
(1999) Spatial gradients of cellular phospho-proteins
FEBS Letters 457:452–454.
https://doi.org/10.1016/S0014-5793(99)01058-3
- PubMed
- Google Scholar
(2005) Spatial coordination of spindle assembly by chromosome-mediated signaling gradients
Science 309:1373–1376.
https://doi.org/10.1126/science.1115964
- PubMed
- Google Scholar
1. Chio US
2. Cho H
3. Shan SO
(2017) Mechanisms of tail-anchored membrane protein targeting and insertion
Annual Review of Cell and Developmental Biology 33:417–438.
https://doi.org/10.1146/annurev-cellbio-100616-060839
- PubMed
- Google Scholar
1. Cui W
2. Mehta P
(2018) Identifying feasible operating regimes for early T-cell recognition: the speed, energy, accuracy trade-off in kinetic proofreading and adaptive sorting
PLOS ONE 13:e0202331.
https://doi.org/10.1371/journal.pone.0202331
- PubMed
- Google Scholar
1. Driever W
2. Nüsslein-Volhard C
(1988) A gradient of bicoid protein in Drosophila embryos
Cell 54:83–93.
https://doi.org/10.1016/0092-8674(88)90182-1
- PubMed
- Google Scholar
1. Ehrenberg M
2. Blomberg C
(1980) Thermodynamic constraints on kinetic proofreading in biosynthetic pathways
Biophysical Journal 31:333–358.
https://doi.org/10.1016/S0006-3495(80)85063-6
- PubMed
- Google Scholar
1. Ellgaard L
2. Helenius A
(2003) Quality control in the endoplasmic reticulum
Nature Reviews Molecular Cell Biology 4:181–191.
https://doi.org/10.1038/nrm1052
- PubMed
- Google Scholar
1. Feric M
2. Brangwynne CP
(2013) A nuclear F-actin scaffold stabilizes ribonucleoprotein droplets against gravity in large cells
Nature Cell Biology 15:1253–1259.
https://doi.org/10.1038/ncb2830
- PubMed
- Google Scholar
(2001) T-cell receptor binding kinetics in T-cell development and activation
Expert Reviews in Molecular Medicine 3:1–17.
https://doi.org/10.1017/S1462399401002502
- Google Scholar
(2004) Mathematical and computational models of immune-receptor signalling
Nature Reviews Immunology 4:445–456.
https://doi.org/10.1038/nri1374
- PubMed
- Google Scholar
1. Hansen CH
2. Yang D
3. Koussa MA
4. Wong WP
(2017) Nanoswitch-linked immunosorbent assay (NLISA) for fast, sensitive, and specific protein detection
PNAS 114:10367–10372.
https://doi.org/10.1073/pnas.1708148114
- PubMed
- Google Scholar
Book
1. Hill TL
(1977) Free Energy Transduction in Biology: The Steady State Kinetic and Thermodynamic Formalism
Academic Press.
https://doi.org/10.1016/B978-0-12-348250-1.X5001-3
- Google Scholar
1. Hopfield JJ
(1974) Kinetic proofreading: a new mechanism for reducing errors in biosynthetic processes requiring high specificity
PNAS 71:4135–4139.
https://doi.org/10.1073/pnas.71.10.4135
- PubMed
- Google Scholar
(2017) Minimum energetic cost to maintain a target nonequilibrium state
Physical Review E 95:042102.
https://doi.org/10.1103/PhysRevE.95.042102
- PubMed
- Google Scholar
1. Ieong KW
2. Uzun Ü
3. Selmer M
4. Ehrenberg M
(2016) Two proofreading steps amplify the accuracy of genetic code translation
PNAS 113:13744–13749.
https://doi.org/10.1073/pnas.1610917113
- PubMed
- Google Scholar
1. Jaiman A
2. Thattai M
(2020) Golgi compartments enable controlled biomolecular assembly using promiscuous enzymes
eLife 9:e49573.
https://doi.org/10.7554/eLife.49573
- Google Scholar
1. Kalab P
2. Weis K
3. Heald R
(2002) Visualization of a Ran-GTP gradient in interphase and mitotic Xenopus egg extracts
Science 295:2452–2456.
https://doi.org/10.1126/science.1068798
- PubMed
- Google Scholar
(2000) Diffusion control of protein phosphorylation in signal transduction pathways
Biochemical Journal 350:901–907.
https://doi.org/10.1042/bj3500901
- Google Scholar
1. Kholodenko BN
(2003) Four-dimensional organization of protein kinase signaling cascades: the roles of diffusion, endocytosis and molecular motors
Journal of Experimental Biology 206:2073–2082.
https://doi.org/10.1242/jeb.00298
- Google Scholar
1. Kholodenko BN
(2006) Cell-signalling dynamics in time and space
Nature Reviews Molecular Cell Biology 7:165–176.
https://doi.org/10.1038/nrm1838
- PubMed
- Google Scholar
1. Kholodenko BN
(2009) Spatially distributed cell signalling
FEBS Letters 583:4006–4012.
https://doi.org/10.1016/j.febslet.2009.09.045
- Google Scholar
1. Kunkel TA
(2004) DNA replication fidelity
Journal of Biological Chemistry 279:16895–16898.
https://doi.org/10.1074/jbc.R400006200
- Google Scholar
1. Lan G
2. Sartori P
3. Neumann S
4. Sourjik V
5. Tu Y
(2012) The energy-speed-accuracy tradeoff in sensory adaptation
Nature Physics 8:422–428.
https://doi.org/10.1038/nphys2276
- PubMed
- Google Scholar
1. Lan G
2. Tu Y
(2013) The cost of sensitive response and accurate adaptation in networks with an incoherent type-1 feed-forward loop
Journal of the Royal Society Interface 10:20130489.
https://doi.org/10.1098/rsif.2013.0489
- Google Scholar
1. Lev S
(2010) Non-vesicular lipid transport by lipid-transfer proteins and beyond
Nature Reviews Molecular Cell Biology 11:739–750.
https://doi.org/10.1038/nrm2971
- PubMed
- Google Scholar
1. Martin KC
2. Ephrussi A
(2009) mRNA localization: gene expression in the spatial dimension
Cell 136:719–730.
https://doi.org/10.1016/j.cell.2009.01.044
- PubMed
- Google Scholar
1. McKeithan TW
(1995) Kinetic proofreading in T-cell receptor signal transduction
PNAS 92:5042–5046.
https://doi.org/10.1073/pnas.92.11.5042
- PubMed
- Google Scholar
(2012) Speed, dissipation, and error in kinetic proofreading
PNAS 109:12034–12039.
https://doi.org/10.1073/pnas.1119911109
- PubMed
- Google Scholar
(2014) Discriminatory proofreading regimes in nonequilibrium systems
Physical Review X 4:021016.
https://doi.org/10.1103/PhysRevX.4.021016
- Google Scholar
(2004) Activation of endogenous Cdc42 visualized in living cells
Science 305:1615–1619.
https://doi.org/10.1126/science.1100367
- PubMed
- Google Scholar
(2004) Stathmin-tubulin interaction gradients in motile and mitotic cells
Science 303:1862–1866.
https://doi.org/10.1126/science.1094108
- PubMed
- Google Scholar
1. Ninio J
(1975) Kinetic amplification of enzyme discrimination
Biochimie 57:587–595.
https://doi.org/10.1016/S0300-9084(75)80139-8
- PubMed
- Google Scholar
1. Parton RM
2. Davidson A
3. Davis I
4. Weil TT
(2014) Subcellular mRNA localisation at a glance
Journal of Cell Science 127:2127–2133.
https://doi.org/10.1242/jcs.114272
- PubMed
- Google Scholar
1. Qian H
(2006) Reducing intrinsic biochemical noise in cells and its thermodynamic limit
Journal of Molecular Biology 362:387–392.
https://doi.org/10.1016/j.jmb.2006.07.068
- PubMed
- Google Scholar
1. Rao M
2. Okreglak V
3. Chio US
4. Cho H
5. Walter P
6. Shan SO
(2016) Multiple selection filters ensure accurate tail-anchored membrane protein targeting
eLife 5:e21301.
https://doi.org/10.7554/eLife.21301
- PubMed
- Google Scholar
1. Rodnina MV
2. Wintermeyer W
(2001) Fidelity of aminoacyl-tRNA selection on the ribosome: kinetic and structural mechanisms
Annual Review of Biochemistry 70:415–435.
https://doi.org/10.1146/annurev.biochem.70.1.415
- PubMed
- Google Scholar
1. Roy J
2. Cyert MS
(2009) Cracking the phosphatase code: docking interactions determine substrate specificity
Science Signaling 2:re9.
https://doi.org/10.1126/scisignal.2100re9
- PubMed
- Google Scholar
1. Sartori P
2. Pigolotti S
(2015) Thermodynamics of error correction
Physical Review X 5:041039.
https://doi.org/10.1103/PhysRevX.5.041039
- Google Scholar
(2020) Patterning and growth control in vivo by an engineered GFP gradient
Science 370:321–327.
https://doi.org/10.1126/science.abb8205
- PubMed
- Google Scholar
1. Stelling J
2. Kholodenko BN
(2009) Signaling cascades as cellular devices for spatial computations
Journal of Mathematical Biology 58:35–55.
https://doi.org/10.1007/s00285-008-0162-6
- PubMed
- Google Scholar
1. Swain PS
2. Siggia ED
(2002) The role of proofreading in signal transduction specificity
Biophysical Journal 82:2928–2933.
https://doi.org/10.1016/S0006-3495(02)75633-6
- PubMed
- Google Scholar
1. Sydow JF
2. Cramer P
(2009) RNA polymerase fidelity and transcriptional proofreading
Current Opinion in Structural Biology 19:732–739.
https://doi.org/10.1016/j.sbi.2009.10.009
- PubMed
- Google Scholar
1. Toda S
2. McKeithan WL
3. Hakkinen TJ
4. Lopez P
5. Klein OD
6. Lim WA
(2020) Engineering synthetic morphogen systems that can program multicellular patterning
Science 370:327–331.
https://doi.org/10.1126/science.abc0033
- PubMed
- Google Scholar
(1999) Extracellular regulated kinases (ERK) 1 and ERK2 are authentic substrates for the dual-specificity protein-tyrosine phosphatase VHR. A novel role in down-regulating the ERK pathway
The Journal of Biological Chemistry 274:13271–13280.
https://doi.org/10.1074/jbc.274.19.13271
- PubMed
- Google Scholar
(2018) Energy-speed-accuracy relation in complex networks for biological discrimination
Physical Review E 98:012420.
https://doi.org/10.1103/PhysRevE.98.012420
- PubMed
- Google Scholar
1. Wu Y
2. Han B
3. Li Y
4. Munro E
5. Odde DJ
6. Griffin EE
(2018) Rapid diffusion-state switching underlies stable cytoplasmic gradients in the Caenorhabditis elegans zygote
PNAS 115:E8440–E8449.
https://doi.org/10.1073/pnas.1722162115
- PubMed
- Google Scholar
1. Yamane T
2. Hopfield JJ
(1977) Experimental evidence for kinetic proofreading in the aminoacylation of tRNA by synthetase
PNAS 74:2246–2250.
https://doi.org/10.1073/pnas.74.6.2246
- PubMed
- Google Scholar
1. Zhang X-J
2. Qian H
3. Qian M
(2012) Stochastic theory of nonequilibrium steady states and its applications. Part I
Physics Reports 510:1–86.
https://doi.org/10.1016/j.physrep.2011.09.002
- Google Scholar

Article and author information

Author details

Vahe Galstyan

Biochemistry and Molecular Biophysics Option, California Institute of Technology, Pasadena, United States

Contribution
Conceptualization, Formal analysis, Investigation, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-7073-9175
Kabir Husain

Department of Physics and the James Franck Institute, University of Chicago, Chicago, United States

Contribution
Conceptualization, Formal analysis, Investigation, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared
Fangzhou Xiao

Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, United States

Contribution
Conceptualization, Formal analysis, Investigation, Methodology, Writing - review and editing

Competing interests
No competing interests declared
Arvind Murugan

Department of Physics and the James Franck Institute, University of Chicago, Chicago, United States

Contribution
Conceptualization, Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing

For correspondence
amurugan@uchicago.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-5464-917X
Rob Phillips
1. Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, United States
2. Department of Physics, California Institute of Technology, Pasadena, United States
Contribution
Conceptualization, Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing

For correspondence
phillips@pboc.caltech.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-3082-2809

Funding

James S. McDonnell Foundation

Kabir Husain

Simons Foundation

Arvind Murugan

John Templeton Foundation

Rob Phillips

National Institute of General Medical Sciences

Rob Phillips

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Anatoly Kolomeisky, Shu-ou Shan and Erik Winfree for insightful discussions, Soichi Hirokawa and Avi Flamholz for providing useful feedback on the manuscript. We also thank Alexander Grosberg whose idea of a compartmentalized ‘rotary demon’ motivated the development of our model. This work was supported by the NIH Grant 1R35 GM118043-01, the John Templeton Foundation Grants 51250 and 60973 (to RP), a James S. McDonnell Foundation postdoctoral fellowship (to KH), and the Simons Foundation (AM).

Version history

Received: June 26, 2020
Accepted: December 24, 2020
Accepted Manuscript published: December 24, 2020 (version 1)
Version of Record published: January 18, 2021 (version 2)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

5,538

views
475

downloads
14

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Article PDF

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Vahe Galstyan
Kabir Husain
Fangzhou Xiao
Arvind Murugan
Rob Phillips

(2020)

Proofreading through spatial gradients

eLife 9:e60415.

https://doi.org/10.7554/eLife.60415

Share this article

Cite this article

Error correction schemes that operate by delaying product formation.

Dependence of fidelity on the diffusion time scale in the limit of very high substrate localization.

Speed–fidelity trade-off and consequences of having weak substrate gradients.

Power–fidelity relationship when tuning the effective driving force Δ⁢μ for different choices of the diffusion time scale τD.

Proofreading based on substrate gradients formed by spatially separated kinases and phosphatases.

Fidelity of the scheme in different regimes of right and wrong complex densities.

The effective number of proofreading realizations (neff) as a function of L/λS.

Optimal diffusion time scale for different choices of λS.

Discrete-state representation of diffusive transport and substrate binding/unbinding events.

Fidelity enhancement as a function of the effective driving force for varying choices of koffW.

Dependence of fidelity on the catalysis rate in the case where the substrate profile is uniform.

Fidelity as a function of the catalysis rate in an ideal substrate localization setting.

Pareto front of the speed–fidelity trade-off at different levels of catalytic activity.

Linear trade-off between speed and fidelity when tuning the rate of catalysis.

Fidelity as a function of L/λR for different choices of the ratio λR/λW.

Example profiles of free and substrate-bound enzymes.

Consequences of relaxing the ρE⁢(x)≈constant assumption on the Pareto front.

Fidelity as a function of diffusion time scale calculated with and without making the ρE⁢(x)≈constant approximation.

Reduction in the peak effective number of proofreading realizations as a function of pbound.

Energetic performance of the kinase/phosphatase–based mechanism.

Author details

Vahe Galstyan

Contribution

Competing interests

Kabir Husain

Contribution

Competing interests

Fangzhou Xiao

Contribution

Competing interests

Arvind Murugan

Contribution

For correspondence

Competing interests

Rob Phillips

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

Power–fidelity relationship when tuning the effective driving force $Δ μ$ for different choices of the diffusion time scale $τ_{D}$ .

The effective number of proofreading realizations ( $n_{eff}$ ) as a function of $L / λ_{_{S}}$ .

Optimal diffusion time scale for different choices of $λ_{_{S}}$ .

Fidelity enhancement as a function of the effective driving force for varying choices of $k_{off}^{W}$ .

Fidelity as a function of $L / λ_{_{R}}$ for different choices of the ratio $λ_{_{R}} / λ_{_{W}}$ .

Consequences of relaxing the $ρ_{_{E}} (x) \approx constant$ assumption on the Pareto front.

Fidelity as a function of diffusion time scale calculated with and without making the $ρ_{_{E}} (x) \approx constant$ approximation.

Reduction in the peak effective number of proofreading realizations as a function of $p_{bound}$ .