{"id":"4d3de9ec-aace-4688-b374-a27b2a545e15","arxiv_id":"2607.12766","paper_version":2,"verdict":"REJECT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"Entropic depletion forces from crowding molecules are modeled as the initiator of receptor-mediated endocytosis, yielding phase diagrams and an optimal virus radius of 30–60 nm that reduces to prior predictions when the new terms are dropped.","lead":"Endocytosis models usually assume the virus is already touching the cell; this paper says tiny surrounding molecules should push it there first, then the cell's squishy interior slows the wrapping. It uses physics equations to map when a particle can be fully swallowed, claiming the best size for viruses like HIV is about 50 nanometres.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Phase-2 kinetics and all 'solubility' phase boundaries rest on a constant-force Hertz creep solution applied to a depth-dependent force; Eq. (26) even gives t(0)<0, so the log-domain thresholds are not a consequence of the Onsager variational principle.","rationale":"Read in good faith: the paper identifies a real initiation gap and proposes a plausible entropic mechanism, and some standalone checks (e.g., the Asakura–Oosawa limit in Eq. (39)) are legitimate. But the central quantitative engine is not sound as written. The Onsager variational principle is invoked, yet Phase 2 is never derived from it; the creep-compliance route is a constant-load Green's function. Eq. (26)'s t(0)<0 is not a minor notation issue—it shows the formula is the inverse of a load-step response, not the time to reach depth h starting from contact. The 'solubility condition' is merely the log-domain constraint of that flawed inversion. Since all phase boundaries and thresholds derive from it, the headline predictions are unsecured. I did not rely on the reader's dimensional objection to Eq. (3); on inspection Eq. (3) is dimensionally homogeneous. The optimal-size relation (Eq. (36)) also does not validate the new physics because it is obtained by dropping depletion and cytoskeleton contributions. A revision could re-derive Phase 2 with a proper convolution or state-dependent friction and test against the current predictions; until then, rejection is appropriate.","tokens_in":22523,"tokens_out":9160,"duration_ms":95390,"concrete_test":"Take the standard linear solid compliance Φ(t) from Eq. (11) and solve the linear-viscoelastic contact problem for a prescribed force F(h(t)) given by Eq. (24): h^{3/2}(t) = (3/(8R^{1/2})) ∫_0^t Φ(t−s) dF(h(s)), with the appropriate step discontinuity at contact, numerically for R=20–100 nm and the parameters in Table II. Compute the marginal engulfment boundary and tc(R). Compare with Eq. (27) and Figs. 4, 7, and 8. If the corrected boundary differs by more than a few percent—or if complete engulfment occurs where Eq. (30) forbids it—the solubility condition is an artifact.","verdict_should_be":"REJECT","load_bearing_attack":"Eq. (22) and Eq. (13) adopt a^3(t)=3RFΦ(t)/8, the Lee–Radok solution for a step load constant in time, and h=a²/R. But the driving force in the wrapping phase is not constant: Eq. (24) contains −2γπh and −(√R/D)h^{3/2}. For a linear viscoelastic half-space with a time-varying load, the correct relation is a Boltzmann superposition integral, not this algebraic substitution. Inserting F(h) into the constant-load formula and solving for t gives Eq. (26). That formula is not the solution of the Onsager equation ζ(h)h˙=F(h); no Rayleigh dissipation function for Phase 2 is ever minimized. The unphysical character is visible immediately: at h=0, Eq. (26) yields t(0)=−τ ln(1+Ec/Ev)<0, whereas the wrapping process starts at t=0. The paper then turns the domain of this inverted step-load formula (Eqs. 28–30) into an 'Onsager solubility condition.' The lower bound is only the requirement that a constant force equal to the final F(2R) would reach h=2R at infinite time; it is neither necessary nor sufficient for the actual h-dependent force history. Every downstream quantitative result—amin (Eq. 31), ζmin (Eq. 32), size window (Eq. 34), critical stiffness (Eq. 33), and Figs. 7–8—is built on this unidentified assumption. Eq. (36) does not rescue the paper because it is derived after dropping the depletion and cytoskeleton terms, so the claimed match to HIV-1 is independent of the new physics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a unified Onsager-variational model of receptor-mediated endocytosis in which crowding-induced depletion forces initiate particle-membrane contact, and viscoelastic creep of the cytoskeleton controls wrapping kinetics. It derives a critical crowder concentration for initiation, a minimum ligand density, an engulfable size window, a critical cell stiffness, and an optimal virus radius claimed to match HIV-1. The main quantitative output is a kinetic phase diagram built from the 'Onsager solubility condition' applied to an expression for the complete engulfment time.","tokens_in":22940,"tokens_out":5883,"duration_ms":55742,"significance":"If sound, the model would fill an acknowledged gap by explaining how a virus and membrane come into proximity before specific binding, and it would make falsifiable predictions about ligand density, particle size, and cell stiffness. The paper also recovers the Asakura-Oosawa flat-surface depletion result in an asymptotic limit, which is a useful consistency check. However, the central Phase-2 kinetic derivation is not the solution of the stated Onsager equation, and the phase boundaries are obtained from the domain of a logarithm rather than from a variational principle. These defects are load-bearing: Eqs. (31)-(34) and Figs. 7-8 rest on them. The claimed optimal-size/HIV match is derived after removing the new entropic and cytoskeleton terms, so it does not validate the proposed mechanism.","major_comments":[{"comment":"Eqs. (22)-(26) adopt the Lee-Radok step-load solution a^3(t)=(3/8)RFΦ(t) for a constant force, but the driving force in Eq. (24) depends explicitly on h through −2γπh and −(√R/D)h^{3/2}. For a time-varying load the correct viscoelastic contact relation is a Boltzmann superposition integral, not the algebraic substitution F(h). Consequently Eq. (26) is not the solution of ζ(h)h˙=F(h), and no Rayleigh dissipation function for Phase 2 is actually minimized. The unphysical nature is visible at h=0: Eq. (26) gives t(0)=−τ ln(1+Ec/Ev)<0, whereas wrapping must start at t=0.","section":"§II.F-H and Appendix F"},{"comment":"The 'Onsager solubility condition' is simply the statement that the argument of the logarithm in Eq. (27) is positive and not larger than 1. It is the existence domain of the algebraic formula, not a consequence of the variational principle. The lower bound F(2R)>16√2 EcEvR^2/[3(Ec+Ev)] is the condition that a constant force equal to the final force F(2R) would reach h=2R at infinite time; it is neither necessary nor sufficient for the actual depth-dependent force history. Every downstream result—amin in Eq. (31), ζmin in Eq. (32), the quartic size window in Eq. (34), and the phase diagrams in Figs. 7-8—inherits this unidentified assumption.","section":"§III.A and Eqs. (28)-(30)"},{"comment":"The critical concentration formula is dimensionally inconsistent. In Eq. (3)/Eq. (B15), ccrit is a number density with dimensions L^{-3}, while the right-hand side evaluates to dimensions L^{-4} (the bracket 4πκ/R+2πγR has dimensions of force, J/L, not energy). Thus the threshold concentration cannot be correct as written, and the claimed initiation condition is not established.","section":"Eq. (3) and Appendix B.3"},{"comment":"The numerical evaluation of the Phase-1 time is wrong, and the accompanying physical explanation contradicts the formula. Eq. (A11) with Table II gives t1 ≈ 6×10^{-5} s, not 1.0×10^{-9} s as stated in Eq. (A13). Moreover, Appendix A.5 says the depletion force 'becomes singular near contact' and that the 'force diverges as the depletion volume derivative diverges,' but the sub-interval force F1 in Eq. (A9) is constant in h, and V1 is linear in h. This changes tapproach in Eq. (8) and weakens the time-scale argument in §II.C.","section":"Appendix A and §II.C"},{"comment":"The optimal radius Ropt=√(6κ/(a−2γ)) is derived after explicitly dropping the entropic and cytoskeleton terms (F(2R)≈2πRa−4πκ/R−4γπR). The claimed match to HIV-1 therefore does not test the new crowding or viscoelastic physics; it reduces to a previously known membrane-binding balance. This should be acknowledged, and the statement in §IV.C that the match 'validates' the present model is unsupported.","section":"Appendix I and Eq. (36)"}],"minor_comments":[{"comment":"Table II lists τ=ηm/Em=2×10^3 s, but the text in §IV.C estimates τ≈800 s using ηc=4000 kPa·s and Ec=5 kPa. The notation Em vs Ec is also inconsistent between the table and the main text, and the unused parameter δ=5 nm could be removed or described.","section":"Table II and §IV.C"},{"comment":"The caption states that 'the optimal size depends on the ligand-receptor binding energy density,' but the curves are plotted for different Ec at fixed parameters; the optimum is the same for all curves, as the text correctly notes. The caption should be revised.","section":"Fig. 4 caption"},{"comment":"Ref. [58] lacks author names and a full journal identification, and Ref. [52] duplicates Ref. [3]. A few equations use '−' unary spacing ambiguously (e.g., Eq. (A11)), making signs hard to parse.","section":"References"}],"recommendation":"reject","confidential_remarks":"The manuscript contains a promising idea—entropy-driven initiation—and one correct asymptotic limit, but the central kinetic and phase-boundary derivations are not valid as written. The Phase-2 result is not a solution of the stated Onsager equation, the 'solubility condition' is a log-domain artifact, and Eq. (3) has a dimensional inconsistency. Because the phase diagrams, thresholds, and size window all rest on these points, the errors are not local and cannot be fixed within the present scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's what you should know: the qualitative mechanism is genuinely new and worth a look, but the quantitative derivation as written does not hold together. The central phase boundaries and the HIV-size match are not supported by the equations.\n\nWhat's new: the paper identifies a real gap—most endocytosis models start with the particle already touching the membrane—and proposes that depletion forces from crowding agents drive the initial approach. That's a legitimate question, and the two-phase Onsager framework is a reasonable way to frame it. The asymptotic recovery of the Asakura-Oosawa result in the flat-surface limit is a nice algebraic check, and the paper is honest about several limitations.\n\nBut the soft spots are load-bearing, not cosmetic. Eq. (3), the headline critical concentration, is dimensionally inconsistent: c is a number density (m^-3) but the right-hand side has units of m^-4. The numerical value of t1 given in Eq. (A13) is about four orders of magnitude off from what the paper's own formula (A11) gives, and Appendix A.5 says the depletion force 'diverges near contact' even though the derivative of V1 is constant. More seriously, the entire wrapping-time calculation uses the constant-force Lee–Radok solution for a step load, then inserts the depth-dependent force F(h) into it. That is not the solution of the Onsager equation; the correct viscoelastic relation is a Boltzmann superposition integral. The problem shows up immediately: Eq. (26) gives t(0) < 0 at h=0. The 'Onsager solubility condition' used for all the phase boundaries is just the positivity of the argument of that log, not a consequence of the variational principle. And the optimal-size formula in Eq. (36) is derived after dropping the entropic and cytoskeleton terms, so the claimed match to HIV-1 is the old receptor-diffusion result, not a test of the new mechanism.\n\nWhere does that leave it? The qualitative idea deserves a serious referee, because if the kinetics were done correctly, this could be a useful contribution to the endocytosis modeling literature. But the paper in its current form is not close to publishable. I'd send it out rather than desk reject—the mechanism is worth evaluating—but with a clear expectation of major revision: fix the dimensional error, redo Phase 2 with a proper Boltzmann integral, and rederive the phase boundaries from the actual Onsager equation. As it stands, I wouldn't cite it or base any experimental design on its numbers.","headline":"A plausible new mechanism (depletion-driven initiation) buried under load-bearing math errors—worth referee attention, not acceptance as written.","tokens_in":23404,"tokens_out":4751,"would_cite":false,"duration_ms":40270,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Depletion forces from molecular crowding can initiate receptor-mediated endocytosis, and the model's kinetic phase diagram puts the optimal virus radius near HIV-1's ~50 nm.","keywords":["receptor-mediated endocytosis","depletion forces","molecular crowding","Onsager variational principle","viscoelastic cytoskeleton","engulfment phase diagram","optimal virus size","HIV-1"],"falsifier":"Measure uptake of well-characterized nanoparticles as a function of crowder concentration in a cellular or synthetic system: if particles below the predicted Rmin (roughly 22 nm) or above Rmax are internalized at the stated ligand densities, or if uptake persists when crowder concentration falls below the predicted critical value, the phase-boundary predictions are contradicted. A simpler variant: in a vesicle-only assay, check whether increasing crowder concentration produces membrane adhesion even in the complete absence of ligand–receptor binding.","tokens_in":22386,"feed_emoji":"🦠","tokens_out":4572,"duration_ms":46206,"temperature":0.7,"pith_summary":"The paper claims that the long-ignored initiation step of receptor-mediated endocytosis—what brings a virus or nanoparticle close enough to the membrane for ligand–receptor binding—is driven by entropic depletion forces from small crowding molecules. It builds a continuum model on the Onsager variational principle, with engulfment depth as the single generalized coordinate, and combines depletion attraction, ligand–receptor binding, membrane bending and tension, and viscoelastic cytoskeletal resistance. The model yields a kinetic phase diagram with quantitative thresholds: a critical crowder concentration for initiation, a minimum ligand density for complete engulfment, a finite particle-size window, and an optimal radius near 50 nm that shrinks as binding energy increases. If correct, this gives a variational foundation that connects molecular crowding, cell mechanics, and viral size, and explains why HIV-1 sits near the predicted optimum.","feed_headline":"Crowding entropy initiates endocytosis; optimal virus size: 50 nm","feed_subtitle":"Kinetic phase diagram from the Onsager principle explains how crowding gets viruses to the membrane and why HIV-1 sits at the sweet spot.","key_machinery":"The load-bearing object is the Onsager variational principle—a rule that the system's evolution minimizes the sum of free-energy change and dissipated power—applied to a single generalized coordinate, the engulfment depth h(t). The free energy E(h) adds four terms: depletion energy, which is linear in h and gives a constant driving force; ligand–receptor binding energy; membrane deformation from the Helfrich–Canham Hamiltonian; and cytoskeleton deformation via a viscoelastic Hertz contact built with the elastic–viscoelastic correspondence principle. Minimizing the total Rayleigh action gives the kinetic law ζ(h) dh/dt = F(h), and the Onsager solubility condition—the requirement that the loga","core_discovery":"The central claim is that endocytosis does not need pre-existing virus–membrane contact: crowding agents exert a depletion force that drives the particle to the membrane within about 7×10^-5 seconds, solving the initiation problem, after which ligand–receptor binding sustains wrapping against viscoelastic cytoskeletal resistance. From the free energy landscape the paper derives closed-form predictions: a critical crowder concentration for initiation (Eq. 3), a minimum ligand density for complete engulfment (Eqs. 31–32), a quartic equation that defines the engulfable size window (Eq. 34), and an optimal radius Ropt = sqrt(6κ/(a-2γ)) that decreases with binding energy density, giving roughly 4","pith_inferences":["If depletion forces really initiate uptake, then raising the concentration of inert crowding agents (such as dextran or Ficoll) in culture medium should be able to restore or accelerate uptake for particles that otherwise fail—a testable prediction the paper does not state explicitly.","Because the optimal radius is independent of cytoskeletal mechanics while the engulfment time is not, cells could modulate uptake speed without changing their size preference, which suggests a separation of control knobs for drug-delivery design.","The flat-surface limit recovering the Asakura–Oosawa result suggests the entropic driving force is generic, so the same initiation mechanism may apply to other membrane-adhesion events, such as immune-cell phagocytosis or uptake of lipid nanoparticles, beyond viruses."],"forward_implications":["Below a critical crowder concentration, initiation fails: even with high receptor density, no ligand–receptor encounter occurs, so uptake cannot start.","Complete engulfment requires ligand density above a minimum value; below it, the driving force cannot overcome membrane bending and cytoskeletal resistance, so the particle is only partially wrapped.","Only particles within a finite size window (roughly 20 nm to about 100 nm for typical parameters) can be fully internalized; smaller ones are blocked by bending cost, larger ones by cytoskeletal deformation, and engulfment time diverges at the lower bound.","The optimal particle radius is fixed by membrane constants and binding energy density, not by cell stiffness, and it decreases from about 60 nm to about 30 nm as binding strengthens; HIV-1's roughly 50 nm radius sits near the predicted optimum for its ligand–receptor system.","Stiffer cells slow engulfment, narrow the size window, and impose a critical stiffness beyond which no particle size can be completely engulfed."],"fun_headline_variants":["Crowding entropy starts endocytosis; HIV-1 size is optimal","Depletion force initiates uptake: no pre-contact needed","Onsager model predicts HIV-1's 50 nm sweet spot","Entropy drives virus to membrane; binding sustains wrap","Kinetic phase diagram: crowding gets viruses in"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The phase boundaries are obtained by requiring the argument of a logarithm in the wrapping-time formula to be positive and finite—a mathematical existence condition—and the paper treats this as a physical threshold for engulfment; if that interpretive step is not legitimate, the predicted thresholds (critical concentration, minimum ligand density, size window, critical stiffness) do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Crowding entropy starts endocytosis; HIV-1 size is optimal","Depletion force initiates uptake: no pre-contact needed","Onsager model predicts HIV-1's 50 nm sweet spot","Entropy drives virus to membrane; binding sustains wrap","Kinetic phase diagram: crowding gets viruses in"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000178,"raw_usage":{"total_tokens":1182,"prompt_tokens":839,"completion_tokens":343,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":259}},"tokens_in":583,"tokens_out":343,"duration_ms":4646,"temperature":1.0,"reasoning_tokens":259,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T06:20:01.481718+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure uptake of well-characterized nanoparticles as a function of crowder concentration in a cellular or synthetic system: if particles below the predicted Rmin (roughly 22 nm) or above Rmax are internalized at the stated ligand densities, or if uptake persists when crowder concentration falls below the predicted critical value, the phase-boundary predictions are contradicted. A simpler variant: in a vesicle-only assay, check whether increasing crowder concentration produces membrane adhesion even in the complete absence of ligand–receptor binding.","supporting_citations":[],"review_version":2}