Metabolic regulation, compartmentation and system limits

A pool size, a flux and a control coefficient are different

A metabolite pool is an amount or concentration present at a time. A flux is an amount converted per time, or a compatible concentration rate under defined volume conditions. A control coefficient describes how a system quantity responds to a small parameter perturbation in a specified state. These quantities cannot be inferred from one another without a model.

An enzyme may be regulated locally by substrate, product, partner binding or modifications. The effect on a whole network also depends on other reactions, transport and compartments. Calling one step “rate limiting” is not a universal substitute for calculating or measuring system-level control.

An open toy model shares control between steps

Invent an open linear model A ⇌ B → sink with clamped c_A=1.00 mmol L⁻¹ and an effective zero sink concentration. Assign v₁=k₁(c_A−c_B), v₂=k₂c_B, k₁=0.400 s⁻¹ and k₂=0.100 s⁻¹. Fixed volume and reservoir conditions are part of this concentration-rate model.

At steady B, equate the rates. Then c_B,ss=k₁c_A/(k₁+k₂)=0.800 mmol L⁻¹ and J=k₁k₂c_A/(k₁+k₂)=0.0800 mmol L⁻¹ s⁻¹. The intermediate is constant while both nonzero fluxes continue.

The differential flux-control coefficients C_i^J=∂ln J/∂ln k_i are k₂/(k₁+k₂)=0.200 for k₁ and k₁/(k₁+k₂)=0.800 for k₂. They sum to one in this particular model and scaling convention. They are not exact finite-change multipliers or universal physiological constants.

Doubling k₁ alone gives J=0.0888888889 mmol L⁻¹ s⁻¹ and B=0.888888889 mmol L⁻¹. Doubling both coefficients gives J=0.160 and leaves B=0.800. Uniform scaling and a single-step change therefore have different effects.

A bulk average can miss what the enzyme experiences

Invent two fixed compartments with volumes 1.00 and 9.00 μL and free substrate concentrations 10.0 and 0.100 mmol L⁻¹. Put the modeled enzyme only in compartment A. Total substrate amount is 10.0+0.900=10.9 nmol, since 1 mmol L⁻¹×1 μL=1 nmol.

The volume-weighted concentration is (1×10+9×0.100)/10=1.09 mmol L⁻¹. For an assigned local rate fraction v/V=S_local/(2.00+S_local), the enzyme in A sees fraction 10/(2+10)=0.833333333. Substituting the bulk mean would instead give 1.09/(2+1.09)=0.352750809.

Only an extra stipulated equilibration that equalizes unbound substrate without loss between the fixed volumes would make the enzyme see the mean concentration 1.09. Actual cells can have transport, binding and different local environments. The invented two-compartment numbers are not a measured cellular distribution.

A stable pool can support continuing throughput

For an independent amount-pool example, start with 3.00 μmol and supply constant input 8.00 μmol min⁻¹ with outputs 5.00 and 2.00 μmol min⁻¹ over four minutes. The balance is dn/dt=8−5−2=1.00 μmol min⁻¹, so the pool reaches 7.00 μmol. It is accumulating rather than steady.

Holding the same outputs while input becomes 7.00 μmol min⁻¹ gives zero accumulation. This does not make either output zero; a constant pool can have continuing throughput. The balance alone also does not identify a unique enzyme controlling the flux.

These fluxes are supplied and constant over the stated interval. A changing concentration that feeds back on rate, a growing volume or another outlet would require a different dynamical model.

A feedback exponent does not count sites

Specify a phenomenological decreasing response f(P)=1/[1+(P/P_½)^h], with assigned half-response concentration P_½=0.500 mmol L⁻¹, exponent h=2 and unmodified flux V=10.0 μmol min⁻¹. The variable P is a hypothetical feedback quantity, not an identified biological inhibitor.

At P=0.250, 0.500 and 1.00 mmol L⁻¹, response fractions are 0.800, 0.500 and 0.200, giving fluxes 8.00, 5.00 and 2.00 μmol min⁻¹. Doubling P from 0.500 to 1.00 changes flux from five to two, rather than halving it.

The exponent two describes the shape of the stipulated response. It does not prove two binding or catalytic sites, an allosteric mechanism or a stable network. P_½ is the defined half-response parameter for this function; it is not automatically a Michaelis K_m. The static graph alone does not settle dynamic regulation.

Regulation depends on the state and the boundary

Shared control, compartment access, pool accumulation and a feedback response describe different aspects of a system. A complete explanation must specify reservoirs, concentrations, active enzyme states, transport and the perturbation being considered. A model valid for one organism, compartment or state is not automatically universal.

Use independent observations to test the part of the model that supports the proposed function. The four original models here show why a local rate curve or stable metabolite level cannot establish every network claim. They contain no actual cellular control coefficients, intervention dose or physiological performance prediction.

Worked examples

Try it yourself

1. Double k₁ alone, then double both k₁/k₂. Give the steady flux and B in each open-model case.

   Show the answer and reasoning

   k₁ alone: J=0.0888888889, B=0.888888889. Both: J=0.160, B=0.800 in mmol L⁻¹ s⁻¹ and mmol L⁻¹ respectively.

   Use the exact steady formulas. Differential coefficients describe small sensitivities rather than exact doubling effects.

2. Suppose the two compartments now equalize free substrate without loss at fixed volumes. Find the enzyme’s local concentration and assigned rate fraction.

   Show the answer and reasoning

   1.09 mmol L⁻¹ and 0.352750809.

   Conserved 10.9 nmol divided by 10.0 μL gives the common concentration. This equality needs the new equilibration assumption.

3. With continuing outputs five and two μmol min⁻¹, what input stabilizes the pool? Does the balance identify one flux-control enzyme?

   Show the answer and reasoning

   7.00 μmol min⁻¹; the balance does not identify a unique control enzyme.

   Zero accumulation means matched total inflow/outflow, not zero flux or a measured parameter sensitivity.

4. For the feedback function, calculate the flux after P doubles from 0.500 to 1.00. Then find P for a ten-percent response.

   Show the answer and reasoning

   Flux 2.00 μmol min⁻¹, not half of five; ten-percent response at P=1.50 mmol L⁻¹.

   Solve 1+(P/0.500)²=10, so P/0.500=3. This static inversion proves no molecular site count or network stability.

What to carry forward

• Pools, fluxes and differential control coefficients are distinct.
• Control can be shared and state dependent.
• A volume-weighted mean may differ from local substrate.
• Stable pools can have nonzero throughput.
• A phenomenological feedback exponent does not count molecular sites.