{"id":"efd2fb7e-5175-4aa0-b00b-7b1477d09d06","arxiv_id":"2412.00771","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A stochastic simulation of protein-based cell signaling shows that messenger size has efficiency optima set by tradeoffs between synthesis cost, diffusion, degradation, and receptor binding.","lead":"The paper builds a computer model of cells sending protein messages and asks what messenger size is best when synthesis, diffusion, breakdown, and receptor binding all depend on size. It finds that optimal sizes exist, with the best size depending on what you optimize, and argues this may explain why real signaling proteins cluster in a narrow mass range.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The optimal messenger size is produced by a hand-picked degradation scaling (KM=(Lp/2)^δ, δ≈2) that the text says was chosen to compensate other size penalties, so the biological claim stands or falls on an unvalidated exponent.","rationale":"I agree with the reader: the weakest point is the degradation scaling. Eqs. (13)-(14) are the compensation mechanism that makes larger proteins viable; all other size dependencies (synthesis λ=1/4, diffusion L^-0.33, binding L^-1.5) monotonically disfavor larger messengers. Without the steep L^-2 degradation advantage, the simulation has no optimum, as the authors state. The qualitative citations do not determine the functional form, and since no Lp→kDa mapping is provided, the claimed relevance to 8-14 kDa chemokines is not established. A further supporting issue is that the efficiency measures include a factor Lp bits per bound protein, which also favors larger messengers; however, the paper claims robustness to this choice under a smaller diffusion coefficient, so I do not make that the primary attack. The verdict remains CONDITIONAL: the proposed test would settle whether the concern lands. If peaks persist with measured degradation, the paper is much stronger; if they vanish or move outside the chemokine range, the biological claim should be revised.","tokens_in":15188,"tokens_out":10216,"duration_ms":102574,"concrete_test":"Using the same simulation code and remaining parameters, replace Eqs. (13)-(14) with pdeg(Lp) fitted to independent measurements of extracellular protein half-life versus molecular mass (e.g., human plasma protein turnover or recombinant protein stability datasets), mapping Lp to kDa as residue count × ~110 Da. Record the locations of the ηE, ηT, and ηN peaks. If the peaks disappear, shift outside 8-14 kDa, or survive only for δ≈2 imposed by hand, the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on Section III Eqs. (13)-(14): pdeg = Vmax/(KM + Pt), Vmax = 1/Ld, KM = (Lp/2)^δ, with δ≈2. The authors state that this 'rapidly decreasing' degradation probability was assumed 'to compensate the other processes which are suppressing the transmission of larger proteins' (synthesis rate/energy via λ=1/4, diffusion D∝L^-0.33, binding pon∝L^-1.5). Since δ=2 is steeper than the combined penalties, the survival advantage overwhelms them and creates the reported η peaks. The paper explicitly concedes that for δ≪2 and zoff≪10 no optimum appears. The cited biology (renal filtration, proteasomal stability) supports only a qualitative trend toward greater stability with size, not a power-law KM∝L^2 with exponent 2. Moreover, the simulation range Lp≈2-40 corresponds, if Lp is amino-acid number, to peptides below about 5 kDa, far from the 8-14 kDa chemokines the paper motivates. Thus the optimality is an artifact of a tuned parameter unless the degradation-size scaling is independently established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a stochastic simulation model of cell-to-cell chemical communication in which a messenger protein is synthesized in a transmitter cell, diffuses in a three-dimensional domain, is subject to size-dependent degradation, and binds to receptors on neighboring cells. Three signal shapes (step, exponential, power-law) are considered, and the authors define efficiency measures that quantify received information per unit of energy, time, and number of generated proteins, with each bound protein assumed to contribute Lp bits of information. The main result is that these efficiency measures display maxima at certain messenger lengths, with total time and number efficiencies peaking around Lp < 10 and max-time efficiencies peaking at Lp > 10, which the authors interpret as evidence for optimal messenger sizes and relate to the observed narrow mass range of chemokines.","tokens_in":15485,"tokens_out":3101,"duration_ms":33861,"significance":"If the underlying assumptions were independently justified, the paper would provide a conceptually appealing explanation for why secreted signaling proteins cluster in a narrow mass range. The strengths of the manuscript are its explicit, transparent simulation model, the use of three different signal shapes, and the large number of realizations (10^6) for the stochastic averages. However, the central biological claim is not yet supported: the main optimum is produced by a hand-chosen degradation scaling, and the efficiency measures assume an information-per-protein scaling that is not derived from the communication channel. The paper is honest about some of these limitations, but the abstract and conclusions overstate the biological implications.","major_comments":[{"comment":"The degradation probability is modeled as pdeg = Vmax/(KM + Pt) with KM = (Lp/2)^delta and delta = 2, explicitly chosen to compensate for the penalties that larger proteins face in synthesis, diffusion, and binding. The text itself concedes that for delta << 2 and z_off << 10 no optimum appears. This makes the appearance of the optimum a direct consequence of the chosen exponent, and the cited biology (renal filtration, proteasomal stability, glycosylation) supports only a qualitative trend toward greater stability with size, not a quantitative power law with exponent 2. The central claim therefore requires independent empirical determination of delta from degradation data, or at minimum a demonstration that the main qualitative results are robust across a biologically plausible range of delta values. As written, the simulation is internally consistent but the biological optimality result is not established.","section":"Section III, Eqs. (13)-(14)"},{"comment":"All efficiency measures are multiplied by Lp, reflecting the assumption that each bound protein carries Lp bits of information and that the effective number of informative sequences is e^{Lp}. This assumption is load-bearing because it introduces a strong preference for larger proteins that is not derived from the model's communication mechanism: receptor occupancy is a binary or concentration-related readout, and no coding scheme is specified by which a single bound protein would transmit Lp bits of sequence information. The statement in the text that the qualitative behavior persists for measures not proportional to Lp is relegated to the Supplementary Information and is not quantified in the main text. The authors should either derive the Lp information scaling from a concrete signaling/coding model or present the results without this factor as a central figure.","section":"Section IV, Eqs. (18)-(26)"},{"comment":"The paper motivates the study with the observation that chemokines have molecular masses in the 8-14 kDa range, but the simulations cover Lp values of roughly 2-40. If Lp is the number of amino acids, this corresponds to peptides below about 5 kDa, far from the chemokine range; if Lp is an arbitrary length scale, no mapping to molecular mass is provided. The paper should calibrate the protein length Lp to molecular mass and explicitly simulate the range corresponding to the chemokine mass window before claiming that the model explains the observed chemokine size distribution.","section":"Abstract and Section I, Fig. 1"}],"minor_comments":[{"comment":"The bottom panels are described as showing quantities 'when B(t) is maximal at Tmin', but they should refer to Tmax; the x-axis and the text in Section IV indicate these are maximum-time quantities.","section":"Fig. 7 caption"},{"comment":"The quantity Vmax = 1/Ld is introduced as a rate, but all quantities are stated to be dimensionless; please clarify the dimensional interpretation or consistently state that all variables are expressed in dimensionless units.","section":"Section III, Eq. (14)"},{"comment":"The claim in Section V that eta_E(T) decreases monotonically with Lp is not visible in the main text figure panels shown; please ensure the reported qualitative behavior matches the displayed data or specify the parameter regime.","section":"Section IV, Fig. 8"},{"comment":"The Supplementary Information is referenced for results with other parameter values and for efficiency measures not proportional to Lp, but the SI is not included in the manuscript; please ensure it is available to reviewers and that the key robustness claims are summarized in the main text.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper's main selling point is the claim that a simple physical model can explain the narrow mass range of signaling proteins. The largest risk is that the optimum is manufactured by the choice of delta = 2 in the degradation model and by the Lp-bit information assumption. The authors are transparent about some of these limitations, but the abstract and introduction present the result as a biological finding. A major revision that adds independent degradation data or systematically varies delta, and that either justifies or removes the Lp information scaling, is necessary before the paper can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely useful: it puts synthesis cost, diffusion, degradation, and receptor binding into one simulation and asks where messenger size optima might come from. The model is clearly described, the simulations are extensive, and the authors are honest that some compensation among size penalties is required for any optimum to appear. The distinction between efficiencies measured at first binding, peak binding, and over the whole signal is a nice structural insight, and the correlation analysis between free and bound proteins adds texture beyond the main efficiency plots.\n\nThe soft spots are real, and they sit exactly where the abstract overclaims. The degradation probability is pdeg = Vmax/(KM+Pt) with KM = (Lp/2)^δ and δ≈2. The text says this choice was made “to compensate the other processes which are suppressing the transmission of larger proteins,” and it concedes that for δ≪2 and small z_off no optimum appears. So the central peak is largely an output of the assumed degradation scaling, not a derivation of it. The cited biology (renal filtration, proteasomal stability) supports a qualitative trend toward greater stability with size, but not specifically KM ∝ Lp^2. That exponent is doing the work.\n\nThe second structural assumption is the information content per bound protein, taken as ∝ Lp bits. Multiplying every efficiency by Lp will pull optima toward larger sizes even if the underlying transport numbers are flat. The authors mention that some optima survive without that term for a smaller diffusion coefficient, but the Lp-factor is still a strong prior, not a measured property of chemokines.\n\nThere is also a mapping problem. The paper identifies Lp with amino-acid number, and the simulated range Lp=2–40 corresponds to peptides under roughly 5 kDa. The chemokines motivating the study are 8–14 kDa. The authors never connect their dimensionless protein length to molecular mass, so the quantitative link to the biological observation is missing.\n\nNone of this kills the paper as a modeling exercise. The authors explicitly frame it as a “working principle,” and the simulation is internally consistent. But the abstract states that optimal mass values exist for real signals, and that claim is not supported at the same standard as the math. The fix is achievable: validate the degradation-size scaling with independent data, map Lp to kDa, and either justify or drop the Lp-bit information assumption. I would send this to review, because the question is important and the model is a reasonable scaffold — but I would expect the revision to substantially reframe the biological conclusion.","headline":"A clean simulation of a size-dependent messenger tradeoff, but the central optimality is built in by a hand-picked degradation scaling and an assumed information-per-length term, so the chemokine claim is not yet earned.","tokens_in":15990,"tokens_out":1822,"would_cite":false,"duration_ms":20527,"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":"A physical tradeoff among synthesis, diffusion, degradation, and binding predicts optimal messenger protein sizes, explaining why chemokines cluster in a narrow mass range.","keywords":["cellular communication","chemokine size","messenger optimality","protein diffusion","extracellular degradation","receptor binding","stochastic simulation","information efficiency"],"falsifier":"Measure the extracellular degradation half-life of a size series of secreted proteins (for example 2, 4, 8, 16, and 32 kDa) in the same medium. If the half-life does not grow about as the square of protein length over that range, the degradation term that creates the optimal sizes in the model is absent, and rerunning the simulation with measured degradation rates should erase the predicted efficiency peaks.","tokens_in":15006,"feed_emoji":"🧬","tokens_out":12613,"duration_ms":97425,"temperature":0.7,"pith_summary":"The paper argues that the size of a secreted messenger protein is not a free parameter: it determines the energy and rate of synthesis, the speed of diffusion, the lifetime set by extracellular degradation, and the probability of binding a receptor. The authors put these four size-dependent effects into a stochastic simulation of a transmitter cell signaling to neighbors and measure information delivered per unit energy, per unit time, and per protein produced. They find that these efficiencies are non-monotonic in protein length, with time- and protein-number efficiencies peaking for small proteins ($L_p<10$) and peak-occupancy efficiencies favoring larger proteins ($L_p>10$). If this physical account is right, the observed clustering of human chemokines near 8–14 kDa is a generic design principle of chemical communication rather than a sequence-specific accident.","feed_headline":"Chemokine size is set by a physical cost tradeoff","feed_subtitle":"A four-way physical tradeoff puts peak signaling efficiency at protein lengths below and above 10, matching real chemokines.","key_machinery":"The argument is carried by a discrete-time stochastic simulation of synthesis, diffusion, degradation, and binding whose input terms scale with protein length $L_p$: synthesis energy cost $\\Delta E_s$ grows with $L_p$; the diffusion coefficient is $D \\propto L_p^{-0.33}$; the degradation probability is $p_{\\mathrm{deg}} = V_{\\max}/(K_M + P_t)$ with $K_M = (L_p/2)^{\\delta}$ and $\\delta \\approx 2$; and the binding probability is $p_{\\mathrm{on}} = (L_p/2)^{-\\beta}$ with $\\beta = 3/2$. The steep degradation-size relation is the balancing term that lets larger proteins remain competitive despite slower diffusion, costlier synthesis, and weaker binding. Performance is scored by the ratios $\\eta_E = L_p B / E$, $\\eta_T = L_p B / T$, and $\\eta_N = L_p B / N$, treating each bound protein as carrying $L_p$ bits of information.","core_discovery":"The central claim is that optimal messenger sizes emerge from a competition among processes that favor small proteins and one that favors large proteins. Small proteins are cheap to synthesize, diffuse quickly, and bind readily; large proteins survive longer in the extracellular space because the degradation probability is taken to fall steeply with size. Across three signal shapes (step, exponential, and power-law), the simulation produces efficiency maxima: total information per time and per generated protein peak around $L_p<10$, while information per energy at the moment of maximum receptor occupancy and the other peak-occupancy efficiencies peak at $L_p>10$, with the optimum shifting upward as signal duration or collision probability increases. The authors conclude that the narrow molecular-mass window of natural secreted messengers reflects this size–cost tradeoff.","pith_inferences":["Extending the paper's setting, if degradation rates differ by tissue, optimal messenger size should vary across tissues; this could be searched for in proteomic surveys of secreted proteins.","The paper's information measure counts each bound protein as carrying $L_p$ bits; a mutual-information or channel-capacity version could sharpen or shift the optima, since the information per protein is itself size-dependent.","The model leaves out protein–protein interactions and active transport; either could move real messenger sizes away from the diffusion-only optimum, so the narrow chemokine window may reflect extra constraints beyond passive diffusion.","A synthetic test would be to engineer two vesicle-based communication systems with different degradation-size scalings and check whether the transferred signal peaks at the predicted messenger sizes."],"forward_implications":["If the central claim is correct, the narrow mass range of native chemokines (8–14 kDa) is a consequence of the synthesis–diffusion–degradation–binding tradeoff, not a historical accident.","The model predicts two distinct preferred messenger sizes: small messengers ($L_p<10$) for overall time and protein-budget efficiency, and larger messengers ($L_p>10$) for efficiency at the moment of peak receptor occupancy.","Optimal messenger size is tunable: longer signal durations and higher collision probabilities shift the preferred size upward, so signal shape and receptor density can select different messenger sizes.","The same efficiency measures can serve as design rules for biomimetic chemical communication, where synthesis cost, degradation rate, and binding affinity can be adjusted to place the optimum where it is wanted.","Because total energy efficiency decreases monotonically with size in the simulations, energy-limited signaling environments should favor the smallest messengers, while time-limited environments should favor intermediate sizes."],"supporting_citations":[{"why":"Supplies the UniProt mass distribution of 204,088 human proteins showing chemokines confined to roughly 8–14 kDa, the empirical pattern the model seeks to explain.","marker":"[7]"},{"why":"Documents that chemokine molecular masses fall in a narrow 8–10 kDa window, the motivating observation.","marker":"[5, 6]"},{"why":"Provides the size-based Michaelis–Menten degradation model the simulation adopts for the degradation probability.","marker":"[43]"},{"why":"Supplies evidence that extracellular proteasomal degradation is slower for larger proteins, supporting the degradation-size trend.","marker":"[12]"},{"why":"Reports shorter lifetimes for smaller extracellular proteins like cytokines and growth factors, grounding the assumed lifetime-size relation.","marker":"[13]"},{"why":"Gives the Arrhenius-type synthesis relation from which the size-dependent energy cost is derived.","marker":"[36]"},{"why":"Relates binding energy to a mass-dependent force constant, yielding a binding ratio that falls with protein length.","marker":"[49]"},{"why":"Provides the standard result that larger proteins diffuse more slowly, supporting the diffusion-size scaling in the model.","marker":"[8, 9]"}],"fun_headline_variants":["Four-way cost tradeoff sets messenger size optimality","Chemokine size emerges from physical cost balances","Optimal messenger size from synthesis-survival tradeoff","Physical costs pin down optimal messenger size","Cost tradeoff explains messenger size sweet spots"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on the assumption that extracellular degradation becomes much slower for larger proteins, steeply enough (roughly as protein length squared) to compensate for the disadvantages of larger proteins in synthesis cost, diffusion, and binding; if real degradation does not scale this way, the predicted optimal sizes disappear.","fun_headline_variants_meta":{"raw":{"variants":["Four-way cost tradeoff sets messenger size optimality","Chemokine size emerges from physical cost balances","Optimal messenger size from synthesis-survival tradeoff","Physical costs pin down optimal messenger size","Cost tradeoff explains messenger size sweet spots"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000466,"raw_usage":{"total_tokens":2247,"prompt_tokens":790,"completion_tokens":1457,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":406,"completion_tokens_details":{"reasoning_tokens":1388}},"tokens_in":406,"tokens_out":1457,"duration_ms":14474,"temperature":1.0,"reasoning_tokens":1388,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:01:27.867199+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the extracellular degradation half-life of a size series of secreted proteins (for example 2, 4, 8, 16, and 32 kDa) in the same medium. If the half-life does not grow about as the square of protein length over that range, the degradation term that creates the optimal sizes in the model is absent, and rerunning the simulation with measured degradation rates should erase the predicted efficiency peaks.","supporting_citations":[{"cited_title":"Calfon, H","cited_arxiv_id":null,"evidence_quote":"Supplies the UniProt mass distribution of 204,088 human proteins showing chemokines confined to roughly 8–14 kDa, the empirical pattern the model seeks to explain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the size-based Michaelis–Menten degradation model the simulation adopts for the degradation probability."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies evidence that extracellular proteasomal degradation is slower for larger proteins, supporting the degradation-size trend."},{"cited_title":"Sigarjonsson et al., Scand","cited_arxiv_id":null,"evidence_quote":"Reports shorter lifetimes for smaller extracellular proteins like cytokines and growth factors, grounding the assumed lifetime-size relation."},{"cited_title":"Bostrom et al., J","cited_arxiv_id":null,"evidence_quote":"Gives the Arrhenius-type synthesis relation from which the size-dependent energy cost is derived."},{"cited_title":"Shabek et.al, Mol","cited_arxiv_id":null,"evidence_quote":"Relates binding energy to a mass-dependent force constant, yielding a binding ratio that falls with protein length."}],"review_version":1}