{"id":"d3583f16-b072-47cf-874f-a98fb8c26273","arxiv_id":"2501.04862","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Mitochondrial fusion and fission give rise to a quasi-particle-like collective degree of freedom that acts as a low-pass filter, so cells ignore brief stress fluctuations and commit to apoptosis only for persistent stress.","lead":"Cells commit to death by apoptosis based on stress signals, and this paper claims that mitochondria's constant fusion and splitting acts like an electrical low-pass filter that smooths out brief stresses and responds to long-lasting ones. The authors derive a general theory of such compartmentalized systems and test it with microscopy of dying cells.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main-text Eq. (2) has a sign opposite to the supplement's own derivation: adiabatic elimination of s from Eqs. (S50)-(S56) and the mechanical analogue give x_dot = F + gamma F'', not F - gamma F''.","rationale":"The reader's weakest_assumption on the S26 mean-field flux is legitimate, but the more immediately decisive issue is that the main-text equations contradict the supplement's own derivation. This matters because the paper's central claim is literally Eq. (2); if the sign is wrong, the claimed suppression/amplification and shifted fixed points do not follow from the stated equation. The sign error is concrete and testable, and it cannot be dismissed as a different convention because the Langevin equation defines F as the drift velocity and gamma is positive. I verified the supplement's derivation path: the exact Fokker-Planck mean is d_t<c> = <F>, so the Taylor expansion gives +gamma F''; fusion conserves the mean; the median then inherits a +Lambda s term; and the adiabatic elimination gives the plus sign. The three-mass analogue confirms this. The biological low-pass filter may still be true because the full stochastic simulations (Fig. S11, Fig. 5f) do not rely on the printed Eq. (2), but the paper must fix Eqs. (1) and (2) and re-run the qualitative analysis before the central claim can be accepted. This strengthens the reader's CONDITIONAL verdict; it does not change it.","tokens_in":53151,"tokens_out":14293,"duration_ms":151502,"concrete_test":"Re-derive Eq. (2) from Eqs. (S50) and (S53) exactly: write s = <c> - m, insert the exact mean equation d_t<c> = <F> and the median equation d_t m = F(m) + Lambda s, keep all terms, and identify the sign of the surviving gamma F'' term; compare with the three-mass analogue of Section 2.7 by simulating the overdamped system with a force with F'' > 0 and measuring the steady-state sign of s. If the sign is positive, main-text Eq. (2) is wrong and the graphical analysis and low-pass filter argument must be re-run with the corrected equation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing defect is an internal sign contradiction in the central equations. The supplement derives the quasi-particle equation with a plus sign: the mean obeys d_t<c> = F(<c>) + (gamma/2) F''(<c>) (S50), the skew obeys d_t s = -Lambda s + (gamma/2) F''(m) (S54b), and hence d_t m = F + (gamma/2) F'' (S56). The mechanical analogue in Section 2.7 gives the same sign: in steady state the spring displacement s is proportional to +F''. The main text instead states d_t s = -Lambda s - gamma F'' and Eq. (2), x_dot = F(x) - gamma F''(x). Since s enters x_dot additively through Lambda s, the two signs predict opposite curvature-induced drift whenever F'' != 0. The sign cannot be absorbed by a redefinition of gamma, which is defined positive and proportional to D Lambda^{-1}. Main-text Eq. (1) carries a related plus-sign inconsistency with the standard Fokker-Planck form used in the supplement (S15). Even if the mean-field fusion-flux approximation S26 is granted, the advertised universal equation of motion is not the one that follows from the paper's own derivation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a general theory for dynamically compartmentalized stochastic systems. Starting from a population-balance description with chemical reaction noise and Smoluchowski-type fusion/fragmentation, it derives a localized quasi-particle degree of freedom and reduces the multi-scale dynamics to a single equation of motion, xdot = F(x) - gamma F''(x) (main-text Eq. (2)). The framework is then applied to Bax-mediated apoptosis: dynamic mitochondrial fusion/fission is predicted to localize the distribution of mitochondrial Bax concentrations and to create a kinetic low-pass filter that suppresses fast stress fluctuations while facilitating responses to persistent stimuli. The experimental part reports that inhibiting mitochondrial fission with M1 increases cell-to-mitochondrion variability in Bax localization and changes the time course of apoptosis induction in a manner qualitatively consistent with a low-pass filter. The supplement contains the full derivation, a mechanical analogue, and stochastic simulations.","tokens_in":1856,"tokens_out":1973,"duration_ms":125609,"significance":"If the central reduction is correct, the paper's theoretical contribution is substantial: it offers a rare example in which a high-dimensional multi-scale stochastic system is reduced to a one-dimensional equation with a curvature-dependent drift, and it connects that reduction to a concrete biological decision. The authors deserve credit for deriving the quasi-particle equation from microdynamics rather than postulating it, for anchoring parameters to independent measurements (Bax translocation time, MOMP rates, mitochondrial fusion/fission rates) rather than fitting the response curves, and for supporting the predicted sigmoidal response with full stochastic simulations. The experimental design with and without the M1 fission inhibitor is a direct test of the predicted effect. However, the printed central equation is internally inconsistent with the supplement's own derivation, so the significance can be assessed only after the sign and normalization issues are resolved.","major_comments":[{"comment":"The central equations of motion are internally inconsistent. The main text states xdot = F(x) + Lambda s and sdot = -Lambda s - gamma F''(x), leading to Eq. (2), xdot = F(x) - gamma F''(x). The supplement derives instead sdot = -Lambda s + (gamma/2) F'' (Eq. S54b) and d_t m = F + (gamma/2) F'' (Eq. S56); the mechanical analogue in §2.7 gives sdot = -k s + l0^2 F'' and hence a steady-state s proportional to +F''. Since s enters xdot additively through +Lambda s, the two signs predict opposite curvature-induced drift whenever F'' is nonzero. The sign cannot be absorbed by redefining gamma, which is defined positive and proportional to D Lambda^{-1}. The authors must identify the correct sign and propagate it through Eq. (2), Fig. 3, and all qualitative predictions; the factor of two between S54b/S56 and the main text should also be reconciled.","section":"Main text Eq. (2); Supplemental Theory §2.6–2.7"},{"comment":"As printed, Eq. (1) is not consistent with the Fokker-Planck equation used in the supplement. The deterministic term appears as +∂c[(F + ∂cΦ) f], whereas Supplemental Eq. (S15) has the standard form ∂t f = -∂c(F f) + ∂c((∂c D)^T f). With the displayed sign, the drift direction is reversed relative to the microdynamics, and the stationary distribution of the noninteracting problem would be incorrect. If Φ is meant to be an effective potential, the force should appear as -∂cΦ rather than +∂cΦ. Please correct Eq. (1) and verify that the localization argument and the identification of the quasi-particle are consistent with Eq. (S15).","section":"Main text Eq. (1)"},{"comment":"The central closure, J_fus,sr(c) ≈ (µ/2) N f(c,t)(⟨c⟩ - c), rests on a sequence of strong approximations: mean-field factorization f(c1,c2,t) = f(c1,t) f(c2,t), delta-distributed equal-size fragmentation, and a triangle linearization of the angular integral. The supplement itself states that the short-range approximation is accurate only near the mean and fails in the tails (Fig. S1d), and the finite-fragmentation extension further assumes f(c,v,t) = f(c,t) m(v,t). Because Eq. (2) and the universal claim are built on this closure, the authors should provide a quantitative test in the biological regime, for example by comparing the exact Smoluchowski flux with Eq. (S26) for the parameters of Figs. S10–S11.","section":"Supplemental Theory §2.4, Eq. (S26)"},{"comment":"The derivation is stated for weak noise D(c) → 0, but gamma is defined as proportional to D Lambda^{-1} and as the width of the distribution f. In the strict limit D → 0, gamma → 0 and the curvature correction in Eq. (2) vanishes, leaving xdot = F(x). The intended asymptotic ordering should be stated precisely (for example, small but finite noise with gamma held at the stationary variance), so that Eq. (2) is derived in a regime where its extra term is nonzero.","section":"Main text, paragraph preceding Eq. (2)"}],"minor_comments":[{"comment":"The abbreviation hiSPC appears in the Results but the Methods use hiPSC; please make the abbreviation consistent.","section":"Results, BAX immunofluorescence"},{"comment":"Two occurrences of Mann-Withney should read Mann-Whitney (the test is the two-sided Mann-Whitney U test).","section":"Results, Fig. 4 caption and text"},{"comment":"The phrase cells not preptreated should read cells not pretreated.","section":"Results, second experimental paragraph"},{"comment":"The sentence for times much longer than the time scales much longer than Λ^{-1} contains a redundant timescale condition and should be edited to a single clear condition.","section":"Results, paragraph before Eq. (2)"},{"comment":"In the sentence describing Eq. (S26), the reference Eq. (S26 is missing a closing parenthesis.","section":"Supplemental Theory §2.4"},{"comment":"The caption says Block dots are fixed points but should say Black dots; the main text also refers to Fig. 3c, while Fig. 3 appears to have only panels a and b.","section":"Fig. 3"},{"comment":"The text calls the semi-analytic prediction parameter-free after fixing gamma from an independent simulation; this wording is misleading and should be replaced, for example, with with gamma determined from an independent simulation and no parameters fit to the response data.","section":"Supplemental Theory §3.3, near Fig. S11"},{"comment":"The main text refers to an exponent characterizing the suppression, while the figure plots the ratio tau_out/tau_in; the axis labels and the text should be aligned with the definition of responsiveness R in the supplement.","section":"Fig. 5f and Supplemental §3.4"}],"recommendation":"major_revision","confidential_remarks":"The sign inconsistency between the main text and the supplement is serious enough that I cannot recommend acceptance in the current form; however, the supplement's own derivation points to a plus sign, so a consistent correction appears achievable within the scope of a revision. I would also ask the authors to check whether the corrected sign changes the graphical analysis in Fig. 3 and the interpretation of the low-pass filter. The novelty relative to the authors' related preprint (arXiv:2312.09307) and the thesis companion [1] should be clarified in the response letter, since the supplement already discloses a close textual overlap."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one for the idea, not for the equations as printed. The central equation of motion, Eq. (2) in the main text, has a sign that contradicts the supplement's own derivation. The supplement gets d_t m = F + (γ/2) F'' (S56) and the mechanical analogue gives s proportional to +F''; the main text writes d_t s = -Λs - γF'' and x_dot = F - γF''. That is not a typo you can absorb by redefining γ, since γ is positive and set by D/Λ. Eq. (1) also has a drift sign opposite to the standard Fokker-Planck form used in S15. If the sign is wrong, the predicted shift of fixed points and the low-pass filter mechanism are at least quantitatively off, and possibly reversed for realistic potentials.\n\nWhat is genuinely new and worth credit: the paper translates the authors' earlier population-balance theory of compartmentalized stochastic systems into a concrete apoptosis story, predicts a kinetic low-pass filter, and tests it with two experiments (Bax variance under M1, cell death timing). The experiments are qualitative but they are not the whole argument. The parameters are anchored to external measurements and an independent simulation, not fit to the death curves. The supplement is transparent that it condenses prior work (arXiv:2312.09307 and ref [1]), which is fine.\n\nSoft spots beyond the sign: the fusion-flux approximation (S26) relies on mean-field factorization, equal-size fragmentation, and a triangle linearization of the angular integral. Mitochondrial size distributions are broad and fragmentation is not equal-sized; the authors acknowledge some of this in the supplement, but they do not test the approximation against a realistic size distribution in the main text. The experimental part underreports: the 24 h Bax comparison is non-significant and omitted from the results narrative. No code or raw data are shipped, despite the claim that code is available on request; for a theory paper with simulations that is a solvable but real deficiency.\n\nNet: this deserves a serious referee, but the referee should insist on fixing the sign convention and re-running the semi-analytic predictions. As printed, the universal equation advertised in the abstract is not the one derived. Still, the idea that organelle dynamics implement a low-pass filter in apoptosis is plausible and the experiments, once fully reported, are worth having in the literature.\n\nMy recommendation: send to peer review, but with a clear request to correct the sign and deposit code/data. I would not cite it in its current form.","headline":"Sign contradiction in the central equation, but the apoptosis low-pass filter idea is worth taking seriously.","tokens_in":54022,"tokens_out":2676,"would_cite":false,"duration_ms":26955,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Dynamic compartmentalization reduces stochastic cell-death machinery to a single quasi-particle degree of freedom whose kinetics act as a low-pass filter for stress signals.","keywords":["compartmentalization","quasi-particle kinetics","apoptosis","mitochondrial dynamics","stochastic resetting","kinetic low-pass filter","cell fate decision"],"falsifier":"Titrate the mitochondrial fusion-fission rate — with graded doses of the fission inhibitor M1 or Drp1 knockdown — while keeping stress fixed, and measure the variance of membrane-bound Bax across mitochondria within single cells. The theory predicts variance scales as $D/\\mu$ (inverse fusion rate); if the variance stays flat or grows as the fusion rate rises, the mean-field quasi-particle picture is refuted. A complementary test measures the responsiveness ratio (output-to-input correlation time) under a stress signal with known correlation time: the theory predicts this ratio falls monotonically with fusion rate and approaches one only for very slow signals.","tokens_in":52899,"feed_emoji":"⚖️","tokens_out":15725,"duration_ms":137772,"temperature":0.7,"pith_summary":"Cells decide whether to die inside a noisy environment: stress signals fluctuate, and the proteins that execute apoptosis (Bax on mitochondrial membranes) accumulate stochastically on mitochondria that constantly fuse and split. This paper claims that the combination collapses into a single emergent degree of freedom — a quasi-particle whose position in concentration space moves under an effective force $F(x) - \\gamma F''(x)$, where $F$ is the bare biochemical force and the second-derivative term comes from the compartment dynamics. If that is right, compartmentalization is not a passive container but an active filter: the quasi-particle kinetics suppress responses to fast stress fluctuations and amplify slow, persistent ones, converting the bistable Bax switch into a death decision that requires sustained stress. The authors support the claim with fluorescence-microscopy experiments comparing cells with dynamic mitochondria to cells whose fission is inhibited, together with stochastic simulations whose parameters are fixed by measured timescales. The scope is general: any organelle-associated decision whose compartment turnover mixes with reaction timescales — not just apoptosis — would inherit this filtering behavior.","feed_headline":"Cell death collapses to one quasi-particle that filters stress","feed_subtitle":"Dynamic mitochondria suppress fast noise and amplify persistent stress, so cells only die when the signal is real.","key_machinery":"The quasi-particle: a localized mode of the concentration distribution $f(c,t)$ that forms because fusion followed by fragmentation produces an effective mean-field flux $J \\approx (\\mu/2)\\, N f(c,t)(\\langle c \\rangle - c)$, pulling every compartment's concentration toward the ensemble mean against the dispersive noise. Its position $x$ is the geometric median of $f$; its internal deformation $s = \\langle c \\rangle - x$ relaxes at the fusion rate $\\Lambda$ and is driven by local curvature, $\\dot{s} = -\\Lambda s - \\gamma F''(x)$; eliminating the fast variable gives $\\dot{x} = F(x) - \\gamma F''(x)$. The mechanical intuition is three overdamped point masses coupled by springs: the central mass ($x$) feels the discrete second derivative of the force field from its neighbors, which is exactly the third derivative of the potential. The $\\gamma F''(x)$ term carries the argument: it explains why compartmentalized steady states differ from the biochemistry's fixed points, why the response to a step stimulus is sigmoidal, and why fast fluctuations in the stress signal are filtered out with a cutoff set by the fusion-fission rate.","core_discovery":"The central claim is that the time evolution of the whole multi-scale stochastic system is effectively described by a single degree of freedom, $x$, the geometric median of the distribution of compartment concentrations, which obeys $\\dot{x} = F(x) - \\gamma F''(x)$ (main-text Eq. 2). The unusual term $\\gamma F''(x)$ — the local curvature of the bare force weighted by the ensemble variance over the fusion rate — is the fingerprint of compartmentalization: it shifts and deepens the effective potential, removes metastable fixed points of the bare biochemistry, and produces response kinetics that a well-mixed description cannot produce. Applied to apoptosis, the authors write $F$ for the bistable accumulation of membrane-bound Bax, fix $\\gamma$ from published mitochondrial fusion-fission rates and Bax translocation timescales, and predict and observe that the mitochondrial ensemble responds sigmoidally to weak apoptotic stimuli: the response is suppressed on short timescales and facilitated on long ones. They further show that for a fluctuating stress signal $\\eta(t)$, dynamic mitochondria act as a kinetic low-pass filter — the ratio of output to input correlation time drops as the fusion-fission rate rises — so cells distinguish slow, biologically relevant stress from fast, irrelevant fluctuations. The experimental demonstrations are comparative: fission-inhibited (M1-treated) cells lose the localization of mitochondrial Bax concentrations and show the unfiltered, stochastic switching behavior.","pith_inferences":["The paper does not directly test the quantitative prediction that the quasi-particle variance scales as the inverse fusion rate ($\\sigma^2 \\propto D/\\mu$); a graded titration of fusion activity while imaging single cells would turn the qualitative localization result into a quantitative test of the mean-field picture.","Because the filter operates at the ensemble level, cells with very few mitochondria should filter less and die more sporadically under fluctuating stress — a prediction testable in engineered cells with low mitochondrial content.","If the filter is a tunable information-processing device, mitochondrial dynamics may be under selection for decision reliability rather than energetics alone; comparing noise-filtering performance across cell types with different intrinsic fusion-fission rates would probe that hypothesis."],"forward_implications":["Under weak apoptotic stimuli the mitochondrial ensemble responds sigmoidally — suppressed at short times and facilitated at long times — so the kinetics of cell death, not just its biochemistry, determines whether a cell dies.","Dynamic mitochondria localize mitochondrial Bax concentrations in concentration space; blocking fission with M1 destroys this localization and restores independent stochastic switching of individual mitochondria, observed as a high-variability subpopulation.","The mitochondrial ensemble acts as a kinetic low-pass filter: faster fusion-fission suppresses fast stress fluctuations more strongly, so a cell can tune its stress cutoff by changing mitochondrial dynamics.","Steady states of the compartmentalized system need not coincide with fixed points of the isolated biochemistry, so inferences about cell-fate states drawn from well-mixed reaction networks can be wrong in real cells.","The quasi-particle reduction should carry over to other organelle-associated decisions — mTORC1 translocation to lysosomes, endosome maturation, mitochondrial respiration — whenever compartment turnover and reaction timescales mix."],"supporting_citations":[{"why":"Supplies the theoretical framework the paper extends: localization of probability densities in concentration space under compartment interactions, the basis of Eq. (1).","marker":"[14]"},{"why":"Provides the BCL-2 family interaction biochemistry from which the paper derives the bistable Bax force term $F(x)$.","marker":"[18]"},{"why":"Establishes membrane binding by tBid and ordered Bax permeabilization, the pore-formation biology the model captures.","marker":"[19]"},{"why":"Reviews BCL-2 family mechanisms in mitochondrial apoptosis, grounding the MOMP and cytochrome-c release readout.","marker":"[20]"},{"why":"Supplies the timescale of mitochondrial membrane permeability waves during apoptosis, used to argue compartment and reaction timescales mix.","marker":"[27]"},{"why":"Characterizes the mitochondrial fusion-fission cycle kinetics that fix the compartment-dynamics regime ($\\Lambda$) in the theory.","marker":"[31]"},{"why":"Provides the experimentally measured Bax-accumulation timescale (about 15 minutes) used to calibrate the effective model.","marker":"[36]"},{"why":"Introduces the M1 fission inhibitor, the pharmacological perturbation that defines the static-mitochondria control in both experiments.","marker":"[37]"}],"fun_headline_variants":["Cell death shrinks to one quasi-particle that filters noise","Quasi-particle kinetics turn cell death into a low-pass filter","Dynamic mitochondria as quasi-particles control the death switch","One quasi-particle emerges to filter signals and decide apoptosis","Cell suicide collapses to a single quasi-particle stress filter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reduction depends on treating mitochondrial fusion and rapid fragmentation as a single mean-field averaging step that pulls each mitochondrion's Bax concentration toward the ensemble average; if broad mitochondrial size distributions and finite fission rates break that averaging, the quasi-particle equation and the predicted low-pass filter no longer follow.","fun_headline_variants_meta":{"raw":{"variants":["Cell death shrinks to one quasi-particle that filters noise","Quasi-particle kinetics turn cell death into a low-pass filter","Dynamic mitochondria as quasi-particles control the death switch","One quasi-particle emerges to filter signals and decide apoptosis","Cell suicide collapses to a single quasi-particle stress filter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000796,"raw_usage":{"total_tokens":3528,"prompt_tokens":991,"completion_tokens":2537,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":607,"completion_tokens_details":{"reasoning_tokens":2458}},"tokens_in":607,"tokens_out":2537,"duration_ms":19023,"temperature":1.0,"reasoning_tokens":2458,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:24:43.292513+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Titrate the mitochondrial fusion-fission rate — with graded doses of the fission inhibitor M1 or Drp1 knockdown — while keeping stress fixed, and measure the variance of membrane-bound Bax across mitochondria within single cells. The theory predicts variance scales as $D/\\mu$ (inverse fusion rate); if the variance stays flat or grows as the fusion rate rises, the mean-field quasi-particle picture is refuted. A complementary test measures the responsiveness ratio (output-to-input correlation time) under a stress signal with known correlation time: the theory predicts this ratio falls monotonically with fusion rate and approaches one only for very slow signals.","supporting_citations":[{"cited_title":"Stochastic reso- nance,","cited_arxiv_id":null,"evidence_quote":"Supplies the theoretical framework the paper extends: localization of probability densities in concentration space under compartment interactions, the basis of Eq. (1)."},{"cited_title":"Emergent behaviour in multi-particle systems with non-local interactions,","cited_arxiv_id":null,"evidence_quote":"Provides the BCL-2 family interaction biochemistry from which the paper derives the bistable Bax force term $F(x)$."},{"cited_title":"Mitochondria as multifaceted regulators of cell death,","cited_arxiv_id":null,"evidence_quote":"Establishes membrane binding by tBid and ordered Bax permeabilization, the pore-formation biology the model captures."},{"cited_title":"BCL-2 family proteins: chang- ing partners in the dance towards death,","cited_arxiv_id":null,"evidence_quote":"Reviews BCL-2 family mechanisms in mitochondrial apoptosis, grounding the MOMP and cytochrome-c release readout."},{"cited_title":"Bcl-XL Inhibits Membrane Per- meabilization by Competing with Bax,","cited_arxiv_id":null,"evidence_quote":"Supplies the timescale of mitochondrial membrane permeability waves during apoptosis, used to argue compartment and reaction timescales mix."},{"cited_title":"Emer- gence of the Mitochondrial Reticulum from Fission and Fusion Dynamics,","cited_arxiv_id":null,"evidence_quote":"Characterizes the mitochondrial fusion-fission cycle kinetics that fix the compartment-dynamics regime ($\\Lambda$) in the theory."},{"cited_title":"Live-cell imaging to mea- sure BAX recruitment kinetics to mitochondria during apoptosis,","cited_arxiv_id":null,"evidence_quote":"Provides the experimentally measured Bax-accumulation timescale (about 15 minutes) used to calibrate the effective model."},{"cited_title":"Bcl-xL Retrotranslocates Bax from the Mitochondria into the Cytosol,","cited_arxiv_id":null,"evidence_quote":"Introduces the M1 fission inhibitor, the pharmacological perturbation that defines the static-mitochondria control in both experiments."}],"review_version":1}