{"id":"10cab823-9dc0-414f-855b-6a43caff3e03","arxiv_id":"1908.07309","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A step-like template in a mass-conserving reaction-diffusion system localizes a new protein peak to the template edge when the two subdomain nullclines intersect at a point where only one has negative slope.","lead":"This paper shows how a protein pattern can serve as a spatial template, making a second protein accumulate at the template edge, using a simple two-component reaction-diffusion model. It gives a geometric rule to predict when this edge-sensing behavior occurs, with implications for cell division, wound healing, and synthetic biology.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Edge-sensing criterion rests on an untested regional-instability size condition; sufficiency of the nullcline-crossing rule is not established.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing spot: the regional-instability mechanism, which converts the homogeneous-domain nullcline-slope criterion into a prediction of edge-peak formation, is heuristic rather than rigorously derived. My read of the full text confirms that the paper itself flags this: Sec. II C calls the argument heuristic, and the Discussion states that an analytic approach using singular perturbation theory 'may help to cast our heuristic explanation of the localized eigenfunction, based on the concept of regional instability, into a more rigorous argument.' The concern is therefore not manufactured; it is the acknowledged gap between the geometric construction and the dynamical claim. The gap matters because the abstract and introduction present the nullcline-crossing condition as a predictive criterion for edge sensing, while the derivation in Sec. II C establishes at most a necessary condition and only for the specific nullcline arrangement studied. The numerical continuation and the Appendix F phase diagram provide genuine supporting evidence, including an independent phenomenological Cdc42 model, so the central claim is plausible and likely correct in the demonstrated regimes. But the sufficiency of the criterion, and the quantitative size condition in footnote 8, have not been tested against a parameter regime designed to violate them. My proposed check directly probes that untested condition: by shrinking the laterally unstable region below the supposed shortest growing wavelength while preserving the nullcline-crossing geometry, one can determine whether the geometric criterion is sufficient or whether an extra condition is needed. The conditional verdict remains appropriate: the paper makes a credible, well-illustrated mechanistic claim with numerical support, but the most load-bearing explanatory step is not yet backed by a derivation or a dedicated numerical test.","tokens_in":25838,"tokens_out":6615,"duration_ms":70435,"concrete_test":"Choose the Fig. 4 reaction kinetics and vary Dm/Dc (or k_fb) so that the nullcline intersection remains on the negative-slope branch but the size of the laterally unstable region, 2(xE-x0), becomes smaller than the shortest growing wavelength pi/q_min computed from the homogeneous dispersion relation of subdomain A at density n0 (Eq. 8 and footnote 8). At each parameter value, run numerical continuation to determine whether the stable edge-peak branch still exists and compute the leading Sturm-Liouville eigenfunction at the base-state saddle-node. If the edge-peak branch persists with a delocalized or non-edge eigenfunction, the regional-instability size condition is not necessary; if the branch disappears, the published geometric criterion is incomplete without a size caveat.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that edge sensing occurs when the subdomain nullclines intersect where one slope is steeper than -Dm/Dc is supported by the heuristic that a laterally unstable subregion around the base state's inflection point x0, inside subdomain A, grows into an edge peak. The load-bearing step is the transfer of the homogeneous-domain instability criterion (Eq. 8) to a finite region adjacent to the template edge. Three conditions are implicitly assumed, and none is derived: (i) the local equilibrium at x0 with density n0 is unstable with respect to modes fitting in the region; (ii) the region is large enough, 2(xE-x0) > pi/q_min (footnote 8); (iii) the coupling to subdomain B acts only as a perturbation and does not change the instability threshold. The paper explicitly labels this regional-instability argument heuristic and defers a rigorous treatment to future work (Sec. II C; Discussion). Since the main text states only a necessary condition ('A necessary condition to trigger a regional instability at the template edge...'), but the abstract and conclusions treat the nullcline geometry as predicting the operational regime, sufficiency of the criterion is not established. The numerical eigenfunction at one saddle-node (Fig. 5) and the Appendix F phase diagram are consistent, but they do not rule out nullcline arrangements satisfying the crossing condition where the unstable region is too small or the B-side reservoir stabilizes the edge mode.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies two-component mass-conserving reaction-diffusion systems on a one-dimensional domain in the presence of a step-like template that changes the reaction kinetics in two subdomains. Building on the authors' local equilibria theory, it constructs monotonic 'base states' from flux-balance and total-turnover balance (Sec. II A), predicts their saddle-node bifurcation, and shows that non-monotonic patterns with a peak at the template edge or at the system boundary arise (Sec. II B). The proposed mechanism is a regional mass-redistribution instability localized near the template edge, summarized by a geometric nullcline-crossing criterion (Sec. II C, Fig. 6). The results are supported by numerical continuation, finite-element simulations, and a parameter scan of a Cdc42 model (Appendix F), and the moving-template section explores pinning, depinning, and suppression.","tokens_in":26038,"tokens_out":7430,"duration_ms":77333,"significance":"If correct, the paper offers a simple and biologically plausible design rule for edge sensing in mass-conserving reaction-diffusion systems and demonstrates it in a second model. The analysis is transparent: the flux-balance construction is compared with numerical continuation, the heuristic nature of the regional-instability argument is acknowledged, and the geometric criterion is falsifiable. The main limitation is that the central criterion is derived only as a necessary condition, while several claims are phrased as predictions of an operational regime; this gap is acknowledged but not quantitatively resolved.","major_comments":[{"comment":"The main text states, in Sec. II C, that the nullcline-crossing condition is 'a necessary condition to trigger a regional instability at the template edge,' but the Abstract and Discussion present the same condition as predicting when edge sensing is operational. These are different logical claims. Sufficiency would require the three regional-instability assumptions introduced in Sec. II C and footnote 8—local instability of the equilibrium at x0, existence of a region larger than the shortest growing mode, and treatment of subdomain B as a weak reservoir—none of which is derived. Because the criterion is the central claim, either a numerical test of these assumptions in the main model (e.g., varying xE or the nullcline shapes to cross the size threshold) or an explicit restriction of the predictive claim to the models of Fig. 4 and Appendix F should be added.","section":"Sec. II C; Abstract; Discussion"},{"comment":"The size condition 2(xE - x0) > pi/q_min(n0) is never checked numerically or analytically, although it is load-bearing for the transition from base state to edge peak. The single eigenfunction shown at one saddle-node (Fig. 5) and the Appendix F phase diagram cannot rule out nullcline arrangements or parameter values where the unstable region is too small or where the mass-reservoir coupling to subdomain B stabilizes the edge mode. I recommend computing q_min from the dispersion relation of the local equilibrium and testing whether edge sensing is lost when the inequality is violated.","section":"Sec. II C, footnote 8"},{"comment":"The text claims, in Sec. II B, that the edge peak exists 'only when the two nullclines intersect at a point where only one of them has negative slope.' The paper analyzes this for one nullcline arrangement in detail and for one non-edge-sensing arrangement in Appendix B; a general proof for arbitrary nullcline arrangements is not given. Since the flux-balance construction is explicitly restricted to the nullcline shapes of Fig. 2(d), the statement should be qualified accordingly, or the scope of the claim should be narrowed.","section":"Sec. II B and Appendix B"}],"minor_comments":[{"comment":"The criterion is sometimes phrased as 'one nullcline has negative slope,' but the exact condition is that the slope is steeper than -Dm/Dc; the Discussion gives the correct version. Please make the precise condition consistent throughout, or state explicitly that the figures are drawn in the Dc >> Dm limit.","section":"Fig. 6 and Sec. II C"},{"comment":"The notation q_min(n0) for the shortest growing mode should be defined in terms of the unstable band of the dispersion relation, and the factor of pi rather than 2pi in the wavelength condition should be clarified.","section":"Footnote 8"},{"comment":"The label 'Non-montomic turnover' should be 'Non-monotonic turnover.'","section":"Fig. 3"},{"comment":"The monotonicity argument for the relation eta0(nbar) is compressed; writing the derivative of Eq. (11) with respect to eta0 explicitly would make the sign argument easier to verify.","section":"Sec. II A"},{"comment":"Several phase-space sketches are schematic; indicating which panels are actual numerical data, as is done in Fig. 8, would help the reader distinguish exact results from illustrations.","section":"Figs. 3-5"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid extension of the authors' local-equilibria program, and the novelty of the heterogeneous-template extension is clearly presented. My main reservation is framing: the word 'predict' in the abstract is stronger than the derived necessary condition. If the authors add a numerical test of the regional-instability size condition or soften the predictive claims, I would support acceptance. The overlap with the authors' prior work is properly credited."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper is a solid, genuinely useful extension of the local equilibria framework. The new piece is bringing step-like heterogeneity into the flux-balance construction: base states, their saddle-node bifurcation, and then the interpretation of the edge peak as a regional mass-redistribution instability. The geometric criterion (nullclines intersect where one slope is steeper than -Dm/Dc) is simple, checkable, and supported both by numerical continuation in the main model and by a Cdc42 model in Appendix F. I also credit the authors for clearly labelling the regional-instability argument as heuristic and for acknowledging that the analytic construction relies on large-domain approximations.\n\nThe soft spot is exactly the one the stress-test flags. The main text states the crossing condition as a necessary condition; the abstract and conclusions treat it as predicting the operational regime, i.e. sufficient. The load-bearing step—transferring the homogeneous nullcline-slope instability to a finite subregion near the edge, with the size condition in footnote 8 and the B-side acting only as a perturbation—is not derived. The numerics show consistency, not sufficiency. There could well be nullcline arrangements that satisfy the crossing condition but where the unstable region is too small or the B-side reservoir stabilizes the mode. This is not a fatal flaw, because the authors are upfront about the heuristic status, but it is a real gap between the abstract's claim and the analysis presented. I would ask them to either soften the language or run a systematic numerical sweep over the size condition and nullcline arrangements.\n\nMinor issues: no code or data are provided, so the continuation and eigenfunction plots are not independently checkable. The small-domain behaviour is acknowledged and shown in Appendix E. The citation pattern is fine; building on Refs 27 and 28 is justified, and the prior jump-type heterogeneity literature is cited.\n\nWho this is for: anyone modeling mass-conserving reaction-diffusion systems with spatial templates, especially in cell polarization and synthetic biology. The criterion is a useful design rule. It deserves a serious referee. My recommendation: send to peer review, with the expectation of a revision that fixes the necessary/sufficient mismatch and ideally includes a more systematic test of the regional-instability size condition.","headline":"A useful extension of local equilibria theory to step-like templates, with a geometric edge-sensing criterion that is well-supported numerically but whose sufficiency claim outruns the proof.","tokens_in":26602,"tokens_out":2889,"would_cite":true,"duration_ms":29511,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A step-like protein template can localize a downstream protein peak to its edge, and a simple nullcline-intersection rule predicts when.","keywords":["mass-conserving reaction-diffusion","edge sensing","local equilibria theory","flux-balance construction","regional instability","reactive nullcline","protein pattern formation","step-like template"],"falsifier":"Run the two-subdomain system with the template chosen so the nullclines intersect only where both slopes are positive or both negative, sweep the average mass adiabatically through the base-state bifurcation, and look at where the first peak appears: the criterion predicts a boundary peak, never an edge peak. A stable edge peak in that parameter regime would falsify the geometric rule. A second check is geometric: make the unstable patch near the edge narrower than the shortest growing mode; the regional-instability assumption says edge sensing should disappear, so an edge peak there would also be evidence against the mechanism.","tokens_in":25607,"feed_emoji":"📍","tokens_out":8832,"duration_ms":81795,"temperature":0.7,"pith_summary":"This paper asks how a cell's downstream proteins can \"sense\" the edge of an already-established protein pattern, as in Cdc42 polarization next to a bud scar or actin-ring formation around a PIP3 domain. Working with a two-component mass-conserving reaction-diffusion model, the authors extend local equilibria theory—which treats the domain as a chain of diffusively coupled well-mixed compartments—to a domain split by a step-like template with different reaction kinetics on each side. They show that the template can trigger a regional mass-redistribution instability near its edge, and that this instability grows into a stationary protein peak sitting on the edge. They reduce the condition for this edge sensing to a geometric rule: the two subdomains' reaction-balance curves (nullclines) must cross at a point where one curve is steeper than the diffusion-flux-balance line. Because the rule depends only on nullcline shape, it gives a testable design principle for synthetic pattern-forming systems and a way to screen molecular models for edge-sensing competence.","feed_headline":"A step-shaped template puts the protein peak right on its edge","feed_subtitle":"A geometric nullcline-crossing rule predicts which cells can localize downstream patterns to a template's edge.","key_machinery":"The engine is the phase-portrait construction of local equilibria theory: each subdomain's reaction kinetics is encoded in its reactive nullcline $f(m,c)=0$, and any stationary profile must lie in a single flux-balance subspace $c+(D_m/D_c)m=\\eta_0$ (a line whose slope is the diffusion ratio). Within this construction, the criterion $\\chi(\\bar n)<-D_m/D_c$ for lateral instability of a homogeneous state becomes a geometric condition on nullcline slope. The paper generalizes the construction to a step template by overlaying two nullclines, splitting non-monotonic profiles at extrema into monotonic segments, fixing $\\eta_0$ by total-turnover balance in the mesa segment, and locating where the profile crosses a laterally unstable nullcline branch. The object that carries the whole argument is the nullcline-intersection rule: edge sensing operates exactly when the two subdomains' nullclines intersect on a branch whose slope exceeds the flux-balance slope.","core_discovery":"The central claim is that a step-like template localizes pattern formation through a regional mass-redistribution instability. In the stationary monotonic \"base state\" connecting the two subdomains' plateaus at the template edge, the profile can develop an inflection point inside subdomain A; as the average mass is increased this state disappears in a saddle-node bifurcation, and the neutral eigenfunction at the bifurcation is concentrated at the template edge. The phase portrait explains this localization: the profile crosses a branch of the A-nullcline whose slope is steeper than the flux-balance slope $-D_m/D_c$, the same condition that makes a homogeneous steady state laterally unstable. Applying that instability criterion locally to the region around the template edge predicts that a peak accumulates there. The resulting geometric criterion—the two nullclines must intersect at a point where one has negative slope steeper than $-D_m/D_c$—divides edge-sensing from boundary-peak behavior, and the paper verifies it against numerical continuation and against the phenomenological Cdc42 model. A slowly moving template edge keeps the peak pinned; fast pulling depins it, fast pushing suppresses it.","pith_inferences":["Beyond the paper: a smooth (non-step) template edge should weaken the regional instability as the gradient flattens; the paper's own heuristic suggests a threshold steepness below which edge sensing disappears, which one could test by varying edge width.","Beyond the paper: because the criterion is purely geometric, it should transfer to other mass-conserving two-component systems—not only protein attachment-detachment—provided the nullcline and flux-balance slopes are defined; testing it in a reconstituted protein or DNA reaction-diffusion system would be a direct check.","Beyond the paper: the moving-edge results imply that a downstream peak can track an upstream pattern only up to a velocity set by cytosolic diffusion and peak width; in a cell with moving landmarks this sets a bound on when polarization can follow a cue."],"forward_implications":["In the regime where the edge-sensing criterion holds, an adiabatic increase of total mass carries the system from the low-mass base state through the saddle-node bifurcation into a stable peak at the template edge; the same sweep produces a boundary peak when the criterion fails.","The geometric criterion is quantitative enough to screen models: for the phenomenological Cdc42 model it predicts that the template must raise both attachment and detachment rates in one subdomain, matching the dual GEF/GAP role of Abr.","A peak pinned to a slowly moving template edge follows the edge; above a critical edge velocity the two directions behave differently—pulling depins the peak while pushing suppresses it.","Edge sensing persists beyond adiabatic mass changes as long as mass inflow is slower than diffusive transport across the domain; a phase diagram in $(L^2/D_c,\\kappa_s^{-1})$ delimits the regime.","Patterns with multiple peaks are unstable and coarsen into one of the two stable single-peak states, so the template edge acts as a selection mechanism between otherwise competing boundary-localized states."],"supporting_citations":[{"why":"Introduces local equilibria theory for mass-conserving reaction-diffusion systems, the phase-portrait framework this paper extends.","marker":"[27]"},{"why":"Provides the flux-balance construction, the nullcline-slope instability criterion (Eq. 8), and the homogeneous-domain phase portraits the edge-sensing analysis builds on.","marker":"[28]"},{"why":"Phenomenological Cdc42 model used in Appendix F to demonstrate that the geometric edge-sensing criterion predicts the observed regime.","marker":"[4]"},{"why":"Supplies the numerical continuation methods used to compute the bifurcation diagrams that verify the flux-balance construction.","marker":"[42]"}],"fun_headline_variants":["Mass redistribution pins protein peak to template edge","Nullcline geometry puts protein peak at template edge","Step template directs protein peak to its edge","Edge localization via mass redistribution instability","Protein peak forms at template edge due to nullcline slope"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The edge-peak prediction rests on assuming that a patch of the region next to the template edge can be treated as a small isolated homogeneous system, becoming unstable when its local reaction-balance curve is steeper than the diffusion-ratio line and the patch is larger than the shortest growing wavelength; the paper itself flags this regional-instability argument as heuristic, with a rigorous asymptotic derivation left for future work.","fun_headline_variants_meta":{"raw":{"variants":["Mass redistribution pins protein peak to template edge","Nullcline geometry puts protein peak at template edge","Step template directs protein peak to its edge","Edge localization via mass redistribution instability","Protein peak forms at template edge due to nullcline slope"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000714,"raw_usage":{"total_tokens":3277,"prompt_tokens":1081,"completion_tokens":2196,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":697,"completion_tokens_details":{"reasoning_tokens":2128}},"tokens_in":697,"tokens_out":2196,"duration_ms":17179,"temperature":1.0,"reasoning_tokens":2128,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:20:48.243562+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the two-subdomain system with the template chosen so the nullclines intersect only where both slopes are positive or both negative, sweep the average mass adiabatically through the base-state bifurcation, and look at where the first peak appears: the criterion predicts a boundary peak, never an edge peak. A stable edge peak in that parameter regime would falsify the geometric rule. A second check is geometric: make the unstable patch near the edge narrower than the shortest growing mode; the regional-instability assumption says edge sensing should disappear, so an edge peak there would also be evidence against the mechanism.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces local equilibria theory for mass-conserving reaction-diffusion systems, the phase-portrait framework this paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the flux-balance construction, the nullcline-slope instability criterion (Eq. 8), and the homogeneous-domain phase portraits the edge-sensing analysis builds on."},{"cited_title":"First, a peak forms at the domain boundary of subdomain A, as in the 16 Regional lateral instability (a) (b) (c) (d) FIG","cited_arxiv_id":null,"evidence_quote":"Phenomenological Cdc42 model used in Appendix F to demonstrate that the geometric edge-sensing criterion predicts the observed regime."},{"cited_title":"Doelman, P","cited_arxiv_id":null,"evidence_quote":"Supplies the numerical continuation methods used to compute the bifurcation diagrams that verify the flux-balance construction."}],"review_version":1}