{"id":"a34325a3-c456-4a41-904c-46b4372ea1fc","arxiv_id":"2501.13837","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For membrane-bound receptors, the paper derives spatial versions of Michaelis-Menten, substrate competition, competitive inhibition, and uncompetitive inhibition rate laws, showing well-mixed formulas overestimate influx when receptor trapping is strong.","lead":"Receptors on cell membranes bind diffusing molecules, but classical well-mixed enzyme kinetics ignores the spatial crowding of these receptors. This paper derives analytical corrections for four receptor reaction schemes and shows where the classical formulas overestimate uptake rates.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline V≪V regime is exactly where the discrete-receptor homogenization plus mean-field availability fraction is least secure; Eqs. (19)-(34) are not validated against discrete-receptor simulations, so the central quantitative claim is conditional.","rationale":"I checked the algebra: the root of the quadratic (12)-(14) is algebraically equivalent to Eq. (19) with φ in (20), and the competition and inhibition analogs follow the same pattern. So there is no internal inconsistency in the derivation from the homogenized PDE-ODE model. The vulnerability is the model itself. The homogeneous Robin condition is standard, but the extension to finite kinetics via u/u0 is not derived from a systematic homogenization of a surface with mixed absorbing, reflecting, and occupied patches; it is an ansatz. Because the paper's strongest biological claim is the large correction at κ≫1, and κ≫1 is the regime where receptor competition is strong, the unvalidated ansatz is load-bearing. A discrete simulation is a direct and feasible test. The reader identified the same assumption, and I agree with that identification. The paper remains a valuable theoretical contribution, but its central quantitative claims should be accepted only conditionally on the discrete-receptor check.","tokens_in":15051,"tokens_out":11984,"duration_ms":116276,"concrete_test":"Simulate the discrete-receptor problem: a sphere of radius R=1 μm with N=10^4 disk receptors of radius εR=10^-3 μm, substrate diffusivity D=500 μm²/s, binding rate ka=4εDR, catalytic rate kc=10^3/s, and fixed far-field concentration S, for S/K spanning 10^-2 to 10^2. Use Brownian dynamics with surface reactions or finite elements with mixed boundary conditions on the disks and irreversible ES→E+P, and measure the steady-state total influx. Compare to Eq. (19). Repeat with receptors placed in tight clusters (e.g., 100 clusters of 100 disks) to test the roughly-evenly-distributed assumption. If the discrete-simulation influx matches Eq. (19) to within about 10% at S/K=0.1 and κ≈3.2, the mean-field homogenization is validated; if the error grows as S/K decreases, the headline overestimate factor needs revision.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central result is the closed-form correction φ in Eq. (19) and its analogs, and the conclusion that V≪V for S/K not much larger than κ≫1. That conclusion rests on two successive approximations. First, Section 2.1 replaces discrete absorbing receptor disks by the uniform Robin condition D∂r c = (κD/R)c with κ=εN/π; this is standard homogenization valid for N≫1, ε²N/4≪1, and roughly uniform receptor placement. Second, Section 2.2 assumes finite receptor kinetics can be inserted by replacing the boundary condition with D∂r c = ka u c, i.e. by multiplying the fully available flux by the mean-field available fraction u/u0 while using the homogenized surface concentration c(R) in the ODEs. This second step is a closure assumption: it posits that the local flux to each available receptor is unchanged except for the overall fraction, ignoring spatial correlations in receptor occupation and any modification of the depletion field around a receptor when neighbors are occupied. The regime highlighted by the paper—κ≫1 and S/K not large—is precisely the regime of strong receptor competition, where these correlations are most likely to matter. The paper does not test the homogenized model against a discrete-receptor simulation or experiment, so the exact inequalities (21)-(23) and the numerical factors in Eqs. (19), (26), (30), and (34) are not established for that regime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops a spatial model for the steady-state influx of diffusing substrate molecules into membrane-bound receptors. Using boundary homogenization, the authors replace discrete receptor disks on a spherical cell by a uniform Robin boundary condition and couple the resulting bulk diffusion PDE to surface ODEs for receptor occupancy. For Michaelis-Menten kinetics, substrate competition, competitive inhibition, and uncompetitive inhibition, they derive closed-form reaction rates that have the same functional form as the classical well-mixed rates but with a concentration-dependent correction factor (φ or ψ) multiplying the half-saturation constant. They show that the spatial rate is always below the well-mixed rate, identify parameter regimes in which the well-mixed rate overestimates uptake, and illustrate the effect with biophysical parameters.","tokens_in":15232,"tokens_out":15933,"duration_ms":141430,"significance":"The derivations are systematic and the steady-state algebra is verifiable: each scenario reduces to a quadratic whose relevant root yields the stated rate formula, and no parameters are fitted to the target results. The paper provides a unified analytical treatment of four common receptor-kinetic schemes and makes falsifiable predictions that could be tested against discrete-receptor simulations or experiments. If the homogenized model is valid in the strong-competition regime, the correction factors in Eqs. (19), (26), (30), and (34) give a practical way to extend classical enzyme kinetics to membrane-bound receptors. The main caveats are that the parameter regime for the headline overestimate is stated too broadly and that the mean-field closure underlying the model is not numerically validated.","major_comments":[{"comment":"The parameter regime for the headline overestimate is misstated. The paper claims V ≪ Vbar if S/K is not much larger than κ and κ ≫ 1, but this is not sufficient. For example, if S/K = κ with κ ≫ 1, Eq. (20) gives φ ≈ sqrt(κ), so Vbar/V → 1 and there is no strong overestimate. A sufficient condition is more restrictive: roughly S/K ≪ κ, and for S/K ≫ 1 actually S/K ≪ sqrt(κ) is needed for the ratio Vbar/V to grow. Please correct the statements in the abstract, the Introduction, Eq. (23), and the analogous conditions for competitive inhibition in Eq. (31) and uncompetitive inhibition in Eq. (35).","section":"Section 1; Section 3.1, Eq. (23)"},{"comment":"The model relies on two successive approximations: boundary homogenization, which replaces discrete receptor disks by a uniform Robin condition, and a mean-field closure, which multiplies the boundary flux by the available-receptor fraction u/u0. This closure neglects spatial correlations between receptor occupancy and the local depletion field around each receptor. The regime highlighted in Eqs. (21)-(23), with κ ≫ 1 and S/K small, is precisely where receptor competition and such correlations are strongest. The manuscript does not test the homogenized PDE-ODE model against a discrete-receptor simulation or against the original mixed-boundary problem. Please either add a numerical validation (for example, Brownian dynamics with N finite-rate absorbing disks) or state explicitly that the exact correction factors in Eqs. (19), (26), (30), and (34) are conditional on this approximation and discuss the expected error.","section":"Section 2.2, Eq. (10)"}],"minor_comments":[{"comment":"The sentence 'an important avenue for future would is to develop' contains a grammatical error and should read 'would be to develop'.","section":"Section 4"},{"comment":"The substrate concentration axis label appears as '[7M]' in the manuscript; this should render as 'µM'.","section":"Figures 3 and 4"},{"comment":"The well-mixed case is described as receptors diffusing in three dimensions, but the association rate ka = 4εRD used for both the well-mixed and spatial rates is the rate for a disk on a reflecting plane. Please clarify that the comparison fixes the intrinsic receptor binding rate rather than using the three-dimensional diffusion-limited rate of a free disk in solution.","section":"Introduction"},{"comment":"The statement 'V ≈ Vbar if S/K ≫ κ or κ ≪ 1' could be made more precise by specifying the relevant limits, for example by writing Vbar/V → 1 as S/K → ∞ with κ fixed, or as κ → 0 with S/K fixed.","section":"Section 3.1, Eq. (22)"}],"recommendation":"major_revision","confidential_remarks":"The paper contains sound algebra and a useful framework for spatially corrected receptor kinetics. My main concerns are the overly broad asymptotic claim in Eq. (23) and the absence of numerical validation of the mean-field closure; both are addressable in revision. I would look favorably on a revised version that corrects these points."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know one thing before reading: the paper is careful and the derivations hold up. It extends the authors' earlier finite-receptor kinetics work to substrate competition, competitive inhibition, and uncompetitive inhibition, and it does so without fitting anything. The central message—that well-mixed Michaelis-Menten theory overestimates membrane-receptor influx, often by a factor up to 1+kappa—is real and analytically proven from the model's assumptions.\n\nWhat's genuinely new: the competition and inhibition formulas (Eqs. 26, 30, 34) and the explicit parameter regimes where the well-mixed approximation fails. The single-substrate MM result (Eq. 19) is a rederivation/extension of Ref. [20], but the paper is honest about that. The algebra is clean; each scenario reduces to a quadratic whose root gives the stated rate. The inequalities V < V and the asymptotic overestimate factor follow from the closed form. Good use of biophysical parameters, and the Monod-equation discussion is a nice hook.\n\nThe soft spot is exactly where the stress-test note lands. The headline regime—low S/K and large kappa—is where the boundary homogenization and the mean-field availability fraction are least secure. The paper assumes that the local flux to each receptor is captured by multiplying the fully available homogenized flux by u/u0. That ignores spatial correlations in receptor occupation and depletion, which are likely strongest when receptors compete hard for substrate. The paper presents no discrete-receptor simulation or experimental test. This is a genuine gap, but it is a gap in validation, not a fatal flaw. The qualitative overestimate is robust and intuitive; the exact numeric factors in Eqs. (19)-(34) are conditional in that regime. Also, for moderate kappa (~3, as in the E. coli example), \"V << V\" is more like a factor of 4 than an order of magnitude. The paper's wording overstates slightly, though the asymptotic statement is fine.\n\nMy read: this deserves serious peer review. A good referee should push for a discrete-receptor simulation in the high-kappa regime and a sensitivity check on the homogenization closure. But the analytical core is sound, the extensions are useful, and the writing is clear. This will be of interest to anyone modeling membrane receptor kinetics or deciding when to trust well-mixed approximations.","headline":"Solid, clearly written extension of boundary-homogenization kinetics to four receptor schemes; the math checks out, but the exact correction factors are not tested against discrete-receptor simulations in the regime where they matter most.","tokens_in":15826,"tokens_out":2255,"would_cite":true,"duration_ms":23034,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92C45","35B27","35Q92"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives a spatial correction showing that classical well-mixed Michaelis–Menten rates overestimate membrane-receptor uptake, with the gap largest at low substrate and high receptor trapping.","keywords":["membrane receptors","Michaelis-Menten kinetics","boundary homogenization","reaction-diffusion equations","substrate competition","competitive inhibition","uncompetitive inhibition","nutrient uptake kinetics"],"falsifier":"Run a three-dimensional Brownian-dynamics or finite-element simulation of the discrete problem—$N$ absorbing receptor disks of radius $\\varepsilon R$ on a sphere, reversible binding with catalytic conversion at rate $k_c$, far-field substrate $S$—and compare the measured steady-state influx with Eq. (19) across the regime $S/K \\lesssim \\kappa$, $\\kappa \\gg 1$; if the simulated influx matches the classical well-mixed rate rather than the corrected rate, the homogenization correction fails.","tokens_in":14756,"feed_emoji":"🧫","tokens_out":8250,"duration_ms":71709,"temperature":0.7,"pith_summary":"This paper tries to show that the textbook reaction-rate laws for membrane-bound receptors—Michaelis–Menten kinetics, substrate competition, competitive inhibition, and uncompetitive inhibition—are systematically too optimistic because they ignore the fact that receptors are packed on a two-dimensional surface. The authors replace discrete receptor disks with a homogenized boundary condition and solve the coupled diffusion-reaction problem at steady state. The result is that the classical rate laws keep their overall shape, but the half-saturation constant must be multiplied by a function $\\phi(S/K)$ that is always greater than 1. At low substrate or high receptor trapping $\\kappa$, the well-mixed rate can exceed the true spatial rate by up to a factor $1+\\kappa$; at high substrate, the correction disappears. If correct, this gives a parameter-based test for when spatial modeling matters and when the classical formulas can be trusted.","feed_headline":"Classical reaction rates overestimate membrane receptor uptake","feed_subtitle":"Receptors on a cell surface compete for substrate; the new correction can cut predicted uptake severalfold.","key_machinery":"The load-bearing object is the dimensionless trapping rate $\\kappa = \\varepsilon N/\\pi$, the effective rate at which the homogenized cell surface consumes diffusing molecules, together with the correction function $$\\$\\varphi$(x) = \\frac12\\left(1+\\kappa - x + \\sqrt{(1+\\kappa + x)^2 - 4\\kappa x}\\right).$$ The function $\\phi$ arises as the algebraic root of the steady-state boundary-closure equation and multiplies the classical half-saturation constant $K$; it encodes receptor competition on the surface. Because $\\phi$ decreases monotonically from $1+\\kappa$ at $x=0$ to $1$ at $x=\\infty$, it tells exactly when spatial correlations matter: for $x \\gg \\kappa$ the reaction is catalysis-limited and classical, while for $x \\lesssim \\kappa$ the arrival of substrate is rate-limiting and the surface geometry cuts the rate.","core_discovery":"For a spherical cell with $N$ receptors of radius $\\varepsilon R$, the paper's central discovery is that the steady-state molecular influx follows the classical Michaelis–Menten form $V = V_{\\max} S / (\\phi(S/K)K + S)$, with the same $V_{\\max}$ and $K$ as the well-mixed theory but a dimensionless half-saturation correction $\\phi(S/K) \\in (1, 1+\\kappa)$, where $\\kappa = \\varepsilon N/\\pi$ is the surface trapping rate. The correction comes from the positive root of a quadratic equation that couples bulk diffusion to surface receptor state equations. The paper proves that $V < \\overline{V}$ for every positive substrate concentration and trapping rate, and shows the gap is largest when $S/K$ is not much larger than $\\kappa$ and $\\kappa \\gg 1$: in the low-substrate limit, the relative error approaches $\\kappa$. The same structure repeats for two competing substrates, competitive inhibition, and uncompetitive inhibition, with the ratio $S/K$ replaced by the relevant combination of substrate and inhibitor concentrations.","pith_inferences":["A testable extension: simulate the discrete receptor problem with random versus clustered receptor placements; if clustering changes the flux beyond what the mean available-receptor fraction predicts, the homogenized correction would need arrangement-dependent modifications.","Beyond the paper, fitting classical Michaelis–Menten curves to spatial uptake data would produce a substrate-dependent apparent $K$; interpreting that dependence through $\\phi$ could let experimentalists estimate the effective trapping rate $\\kappa$ from dose–response data.","The analysis assumes isolated cells; for dense cell clusters, shielding between neighboring cells would likely compound the overestimate, making the well-mixed rate even less reliable than the single-cell correction indicates.","The same boundary-homogenization machinery should extend to non-spherical or corrugated membranes by redefining $\\kappa$ through local geometry, although the paper does not address that case."],"forward_implications":["Wherever receptors are dense on a small cell ($\\kappa \\gg 1$) and substrate is scarce ($S/K \\lesssim \\kappa$), classical Michaelis–Menten uptake can overestimate the true influx by up to a factor $1+\\kappa$; for the paper's illustrative bacterial parameters, $\\kappa \\approx 3.2$ and the overestimate appears below roughly 1 $\\mu$M substrate.","In the opposite limit, $S/K \\gg \\kappa$ or $\\kappa \\ll 1$, the classical formulas are safe: the correction $\\phi \\to 1$ and the well-mixed rate is recovered.","For substrate competition, the same correction applies with $x = S_1/K_1 + S_2/K_2$, so it is the combined substrate load that determines when spatial effects matter.","Both competitive and uncompetitive inhibitors shrink the parameter region in which the well-mixed rate badly overestimates, because inhibitor occupancy lowers the effective catalytic demand on the surface.","The derivation gives a mechanistic route to the empirical bacterial growth–nutrient uptake law: the uptake rate is the corrected Michaelis–Menten form, not the classical one."],"supporting_citations":[{"why":"Supplies the boundary-homogenization approximation that replaces discrete receptor disks by a uniform reactive boundary.","marker":"[19]"},{"why":"Extends that homogenization to finite receptor kinetics and yields the association rate $k_a = 4\\varepsilon D R$ used throughout.","marker":"[20]"},{"why":"Defines the dimensionless trapping rate $\\kappa = \\varepsilon N/\\pi$ that controls the correction.","marker":"[18]"},{"why":"Gives the classical well-mixed Michaelis–Menten rate that the spatial rate is compared against.","marker":"[7]"},{"why":"Gives the classical two-substrate competition formula that the spatial version modifies.","marker":"[8]"},{"why":"Gives the classical competitive and uncompetitive inhibition formulas that the spatial versions modify.","marker":"[9]"},{"why":"Fixes the illustrative receptor radius, cell radius, and related biophysical parameters used in the numerical comparison.","marker":"[24, 25]"}],"fun_headline_variants":["Spatial receptor model reveals classical rate overestimates","Well-mixed theory overstated membrane receptor uptake","Surface trapping corrects classic Michaelis-Menten rates","New model cuts predicted receptor uptake severalfold","Spatial effects fix overestimates in receptor kinetics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands on the homogenization step in which many small receptor disks are replaced by a uniform boundary condition with trapping rate $\\kappa = \\varepsilon N/\\pi$, combined with the assumption that occupied receptors only reduce the flux by the mean available fraction $u/u_0$.","fun_headline_variants_meta":{"raw":{"variants":["Spatial receptor model reveals classical rate overestimates","Well-mixed theory overstated membrane receptor uptake","Surface trapping corrects classic Michaelis-Menten rates","New model cuts predicted receptor uptake severalfold","Spatial effects fix overestimates in receptor kinetics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000818,"raw_usage":{"total_tokens":3551,"prompt_tokens":881,"completion_tokens":2670,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":2596}},"tokens_in":497,"tokens_out":2670,"duration_ms":19183,"temperature":1.0,"reasoning_tokens":2596,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:34:57.641893+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a three-dimensional Brownian-dynamics or finite-element simulation of the discrete problem—$N$ absorbing receptor disks of radius $\\varepsilon R$ on a sphere, reversible binding with catalytic conversion at rate $k_c$, far-field substrate $S$—and compare the measured steady-state influx with Eq. (19) across the regime $S/K \\lesssim \\kappa$, $\\kappa \\gg 1$; if the simulated influx matches the classical well-mixed rate rather than the corrected rate, the homogenization correction fails.","supporting_citations":[{"cited_title":"Shoup and A","cited_arxiv_id":null,"evidence_quote":"Supplies the boundary-homogenization approximation that replaces discrete receptor disks by a uniform reactive boundary."},{"cited_title":"Berg and E.M","cited_arxiv_id":null,"evidence_quote":"Extends that homogenization to finite receptor kinetics and yields the association rate $k_a = 4\\varepsilon D R$ used throughout."},{"cited_title":"Mechanism of action of memantine","cited_arxiv_id":null,"evidence_quote":"Defines the dimensionless trapping rate $\\kappa = \\varepsilon N/\\pi$ that controls the correction."},{"cited_title":"Johnson and Roger S","cited_arxiv_id":null,"evidence_quote":"Gives the classical well-mixed Michaelis–Menten rate that the spatial rate is compared against."},{"cited_title":"A Note on the Kinetics of Enzyme Action","cited_arxiv_id":null,"evidence_quote":"Gives the classical two-substrate competition formula that the spatial version modifies."},{"cited_title":"Pocklington and J","cited_arxiv_id":null,"evidence_quote":"Gives the classical competitive and uncompetitive inhibition formulas that the spatial versions modify."}],"review_version":1}