{"id":"620dd53a-c66a-4395-b168-fe3d28bf9310","arxiv_id":"2608.05130","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Coarse-grained molecular dynamics shows that cellulose-like glycan chains slow water's rearrangement into tetrahedral ice-like order at 180 K, supporting a tetrahedrality-based mechanism for antifreeze behavior.","lead":"This simulation study asks how sugar chains called glycans stop water from freezing. It reports that cellulose-like chains disturb the tetrahedral arrangement of nearby water, keeping that water liquid far below the normal freezing point, and could guide the design of sustainable antifreeze materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The CG cellulose model is unvalidated and its large LJ sigma likely creates the low-q hydration layer as an excluded-volume artifact; without a control, the central mechanism is unproven.","rationale":"The paper's strongest claim is that glycans prevent freezing by disrupting tetrahedral water structure. The evidence chain is: q dip near surface (Fig. 2), slower C_q(t) decay at 180 K near CG8 (Fig. 3), and chain-length independence. Every link in this chain uses the q order parameter, so anything that biases q near the solute is load-bearing. The reader identified the coarse-grained model as the weakest assumption; I agree and sharpen it. The Methods do not validate the CG model: no parameterization table, no comparison with atomistic cellulose or M3B, no hydration free energy or contact-angle check. The solute-water interaction is a bare 12-6 LJ on a one-site mW water model, which has no hydrogen-bond donor/acceptor sites; calling epsilon=1.0 kcal/mol 'hydrophilic' is an assertion, not a result. With sigma=4.5 Å (about 1.9 times the mW water sigma), the excluded volume alone will deplete the first water shell. Under those conditions, water molecules with fewer than four neighbors within 3.1 Å are assigned q=0 (§2.2). That assignment can produce exactly the low-q 'disrupted' layer shown in Fig. 2b without any specific glycan chemistry. Similarly, slower C_q(t) at 180 K may just reflect a larger liquid-like region pinned by a large inert obstacle, not an active antifreeze mechanism. The paper reports no ice-fraction comparison, and the two control simulations needed to separate excluded-volume from hydrophilicity are absent. I therefore do not reject the manuscript; the hypothesis is plausible and the prior ab initio work provides some motivation. But the central claim is conditional on a model validation that is missing. The concrete test—a repulsive-bead control plus a small-sigma control—would settle whether the observed disruption is chemistry or geometry. If the paper passes this test and reports ice fractions, the CONDITIONAL verdict could be upgraded.","tokens_in":5604,"tokens_out":10608,"duration_ms":134850,"concrete_test":"Re-run the CG8+water system at 180 K with two control beads: (i) a purely repulsive bead of identical sigma=4.5 Å (epsilon=0) and (ii) a smaller solute sigma=2.8 Å with epsilon=1.0 kcal/mol. Compare the near-surface q(r) profile and C_q(t) for both controls against the published model. If the purely repulsive or small-sigma control reproduces the low-q layer and slowed C_q(t), the disruption signature is an excluded-volume artifact of the unvalidated bead size, and the central claim is not supported; if the q profiles differ substantially from both controls, the specific glycan parameters are implicated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that glycans inhibit freezing by disrupting water's tetrahedral network—rests entirely on the q distributions and q autocorrelations of §2.2–2.3. These simulations use a coarse-grained cellulose model whose only solute–water coupling is an isotropic 12-6 Lennard-Jones term with epsilon=1.0 kcal/mol and sigma=4.5 Å (Methods). No validation of this mapping is reported: the paper asserts 'we ensured that the cellulose-water interaction is hydrophilic' without comparing to atomistic cellulose or to the M3B force field it invokes. This matters because mW water has no hydrogen-bonding sites; solute hydrophilicity can only be encoded in the two-body epsilon, and 1.0 kcal/mol is small compared with the mW water-water two-body epsilon (~6.2 kcal/mol). Meanwhile sigma=4.5 Å is roughly twice the mW water sigma (~2.39 Å), so each bead excludes a large volume. Near such a big repulsive core, water is necessarily under-coordinated; the paper even assigns q=0 to water molecules with fewer than four neighbors within 3.1 Å (§2.2). The observed low-q, 'disrupted' layer within 6 Å of the chain may therefore be a geometric excluded-volume artifact of the bead size rather than a chemical property of glycans. Because this layer is the only direct microstructural evidence for the proposed mechanism, the model fidelity is the load-bearing assumption; it is asserted, not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports coarse-grained molecular dynamics simulations of cellulose-type glycan chains in mW water at 300 K and 180 K, using the tetrahedrality order parameter q to characterize hydration water. The authors find that water within a few angstroms of the glycan chain has lower and more variable tetrahedrality than bulk water, and that the q autocorrelation decays more slowly in the CG8+water system than in pure water at 180 K. From this they conclude that glycans prevent water from freezing by disrupting the tetrahedral rearrangement of water, and propose tetrahedrality as a design principle for cellulose-based antifreeze materials, framed as validation of their earlier ab initio hypothesis.","tokens_in":5883,"tokens_out":5175,"duration_ms":65158,"significance":"If the central claim is correct, the paper would establish a simple, physically interpretable descriptor—hydration-water tetrahedrality—for the antifreeze activity of cellulose-type glycans, with potential implications for sustainable antifreeze design. The work uses an established order parameter (Errington–Debenedetti q), a widely used water model (mW), and publicly available LAMMPS, and it connects to the authors' prior density-functional study of cellulose–ice interfaces. These are genuine strengths. However, the evidence presented is indirect: there is no direct ice-fraction or nucleation measurement, the coarse-grained model is not validated, and the main figures lack statistical characterization. The significance is therefore conditional on substantial additional support.","major_comments":[{"comment":"The claim that glycans 'prevent water from freezing' is a phase-behavior statement, but the manuscript reports no direct freezing observable. The simulations are cooled from 300 K to 180 K at 0.1 K/ns with only 25 ns of production, and no ice fraction, critical nucleus, or freezing/melting temperature is reported. In §2.2 the ice phase is identified by chill+ and removed before the q analysis, which means the remaining water distribution is by construction not ice-like. The C_q(t) decay at 180 K in Fig. 3 is a dynamical correlation, not a measure of the presence or absence of ice. The authors should report ice fraction as a function of temperature or time, or characterize the state of the system (liquid, supercooled, or partially frozen) over the production run, to support the freezing-suppression claim.","section":"Abstract, §2.3, Fig. 3"},{"comment":"The coarse-grained glycan model is not validated for the property it is used to measure. The glycan–water interaction is a single isotropic 12-6 Lennard-Jones term with ε=1.0 kcal/mol and σ=4.5 Å, and the text only asserts that 'we ensured that the cellulose-water interaction is hydrophilic' without comparing with atomistic cellulose, with the M3B force field invoked in ref. [7], or with any experimental hydration data. Because mW water has no explicit hydrogen-bonding sites, solute hydrophilicity is encoded entirely in this two-body term; a well depth of 1.0 kcal/mol is small compared with the mW water–water two-body term, and a sigma of 4.5 Å is roughly twice the mW water sigma. The low-q hydration layer within 6 Å of the chain may therefore be a geometric excluded-volume artifact of a large repulsive bead rather than a chemical property of glycans. A validation study and a control simulation with a non-hydrophilic or hard-sphere solute of the same size are needed to separate these effects.","section":"§4, Materials and Methods"},{"comment":"The assignment of q=0 to every water molecule that has fewer than four neighbors within the 3.1 Å cutoff directly biases the main observable in exactly the region where the model's large excluded volume acts. Near a bead with σ=4.5 Å, many water molecules will be undercoordinated by this criterion, so the 'disrupted' low-q layer is partly definitional. The paper should report the fraction of q=0 molecules as a function of distance from the chain, and should recompute the tetrahedrality using the four nearest neighbors regardless of an absolute cutoff, to distinguish genuine tetrahedral disruption from simple geometric undercoordination.","section":"§2.2, Eq. (1)"},{"comment":"The central quantitative claims are not supported by uncertainty estimates or replicate information. Figure 2 shows color-coded single snapshots with no error bars, and the conclusion that the q distribution 'remains fairly the same' for CG4, CG6, and CG8 is made without a quantitative comparison of distributions. Figure 3 compares C_q(t) decays with no confidence intervals, no number of independent runs, and no statement of equilibration or convergence. The authors should report averages and standard errors over at least several independent trajectories, and provide a direct statistical test for the chain-length dependence and for the difference between CG8+water and pure water at 180 K.","section":"Figures 2 and 3, §2.2–2.3"}],"minor_comments":[{"comment":"The first sentence repeats a phrase: 'Antifreeze materials prevent ice-formation by disrupting the ice-formation by binding to certain ice-planes.' This appears to be a copyediting error.","section":"Abstract"},{"comment":"The color scale for q, the exact definition of the distance cutoff (per-atom minimum distance to any chain bead?), and the number of water molecules included after chill+ removal should be stated in the caption; currently the reader cannot tell whether blue domains are q=0 by cutoff or genuinely low-q water.","section":"Fig. 2"},{"comment":"The sentence 'Here, the variance⟨q(0)⟩ 2 is computed over all water particles over the entire trajectory' is not standard notation; the denominator in Eq. (2) is ⟨q(0)^2⟩−⟨q(0)⟩^2, and the text should say the variance of q(0) is computed, not 'the variance ⟨q(0)⟩^2'.","section":"Eq. (2)"},{"comment":"The bonded interaction parameters (harmonic bond and angle constants) are not given, despite being part of the model introduced in ref. [4]. All force-field parameters should be listed to make the simulations reproducible.","section":"§4, Materials and Methods"},{"comment":"The cooling rate of 0.1 K/ns is very fast relative to typical ice nucleation timescales for mW water; the authors should discuss whether the 180 K state is equilibrated or a glassy/vitreous state, since this directly affects the interpretation of the q autocorrelation.","section":"§4, Materials and Methods"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is heavily self-referential, repeatedly framing the work as validation of the authors' previous hypothesis (ref. [4]) and using the same tetrahedrality concept and a model described as introduced in that paper. This is not equation-level circularity, but the novelty is essentially a simulation follow-up; the broad 'design principle' claim in the conclusion goes beyond what the data demonstrate. An editor may wish to ensure that the advocacy tone is appropriately tempered after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The specific observation here is new: CG cellulose-type glycan chains in mW water slow the decay of the tetrahedrality autocorrelation at 180 K compared with pure water, and the near-chain water is less tetrahedral than bulk. That is a clean, reproducible simulation result, and connecting it to the antifreeze-protein hydration literature is a sensible extension of the authors' prior ab initio work. The chain-length independence is also a useful negative result.\n\nThe central claim, though — that glycans prevent freezing by disrupting tetrahedral structure — goes beyond what the data show. There is no direct ice-fraction or nucleation measurement, only q distributions and one autocorrelation at one temperature. The q=0 assignment for under-coordinated water biases the distribution near the chain, and the CG bead has sigma=4.5 Å, roughly twice the mW water sigma, so the low-q layer may be an excluded-volume artifact rather than a chemical property of glycans. The authors state they “ensured” hydrophilicity, but provide no validation against atomistic cellulose or the M3B model they invoke. Figures 2 and 3 have no error bars or replicate information, so it is hard to tell whether the 180 K difference is robust.\n\nThe paper is not circular in the equation-level sense, and the self-citation to prior work is not itself a problem. The interpretive framework is self-referential, but the q measurements are independent. What is missing is a control: a hydrophobic bead of the same size, or an atomistic benchmark, would distinguish chemical disruption from geometric exclusion.\n\nThe writing is clear and honest about scope; the abstract has a copyedit issue (“disrupting the ice-formation by disrupting the ice-formation”), but that is minor. The paper deserves a serious referee because the proposed design descriptor is testable and the simulations are cheap enough to strengthen with controls. As it stands, the mechanism is plausible, not demonstrated. I would send it out, but expect heavy revision: add replicates and error bars, include a direct freezing metric (e.g., ice fraction or critical nucleus size), and validate the CG mapping against an atomistic model or hydrophobic control.","headline":"Plausible mechanism, unproven: the q autocorrelation result is new and useful, but the excluded-volume artifact and missing direct freezing metric make the central claim premature.","tokens_in":6449,"tokens_out":2506,"would_cite":false,"duration_ms":31952,"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":"This paper claims that cellulose-type glycans keep nearby water from freezing by preventing it from rearranging into a highly tetrahedral, ice-like structure, making tetrahedrality the key measure for antifreeze materials.","keywords":["tetrahedrality","antifreeze","cellulose","glycan","hydration water","molecular dynamics","mW water","ice"],"falsifier":"A fully atomistic molecular dynamics simulation of cellulose in an atomistic water model cooled under the same protocol, computing the same q distribution and C_q(t), would settle the claim: if hydration water in the first shell reaches bulk-like tetrahedrality and its autocorrelation decays as fast as in pure water at 180 K, the coarse-grained result is a model artifact.","tokens_in":1627,"feed_emoji":"❄️","tokens_out":3165,"duration_ms":153512,"temperature":0.7,"pith_summary":"The paper sets out to explain why cellulose-type glycans show antifreeze activity by simulating their hydration water with a coarse-grained water model built around tetrahedral order. Using the tetrahedrality parameter q, it claims that water within roughly 6 Å of the glycan surface is prevented from rearranging into the highly tetrahedral ice structure at 180 K, while water slightly farther away appears ice-like. The paper reports that this effect is local and essentially independent of chain length, and that the time autocorrelation of q decays more slowly around a glycan chain than in pure water at 180 K. The conclusion is that disrupting tetrahedral linkage in water is the central mechanism of glycan antifreeze behavior and a design principle for cellulose-based antifreeze materials.","feed_headline":"Glycans suppress ice by stalling water's tetrahedral rearrangement","feed_subtitle":"A simulation links glycan antifreeze activity to disruption of water's tetrahedral order, offering a screening metric for new materials.","key_machinery":"The key object is the tetrahedrality order parameter q, a number between 0 and 1 that measures how close a water molecule's four nearest neighbors are to a perfect tetrahedron. The simulations pair this with the coarse-grained mW water model, which is constructed to reproduce the tetrahedral ordering of water and its freezing behavior, and a three-bead (A/B/C) coarse-grained representation of each glucose ring. The glycan beads interact with water through a Lennard-Jones potential tuned to be hydrophilic, with hand-set parameters of epsilon = 1.0 kcal/mol and sigma = 4.5 Å. The authors compute spatial q distributions around the chain after removing ice with an ice-identification routine, and the time autocorrelation C_q(t), which reports how long water retains its tetrahedral state; the slow decay of C_q(t) at 180 K near the glycan is the direct evidence that rearrangement to tetrahedral ice is inhibited.","core_discovery":"The central claim is that the degree of tetrahedrality of hydration water is the intrinsic measure of glycan antifreeze behavior. In simulations, water near a cellulose-type glycan chain at 180 K does not reorganize into the tetrahedral coordination characteristic of ice; the q parameter shows disrupted order within the first hydration shell and ice-like order beyond it, and the q autocorrelation decays more slowly than in pure water. Since the q distribution does not change as the chain grows from four to eight repeat units, the effect is localized at the glycan–water interface. The authors interpret this as validation of their earlier ab initio finding that cellulose binds ice planes through tetrahedral coordination, and conclude that suppressing tetrahedral rearrangement is the design principle for cellulose-based antifreeze materials.","pith_inferences":["If tetrahedrality is the controlling descriptor, the same q-based analysis could rank other polysaccharides and synthetic hydrophilic polymers by antifreeze potency without simulating ice growth, a screening use the paper does not explicitly develop.","The chain-length independence hints that tethering cellulose-like beads to a surface or nanoparticle could confer local frost resistance, extending the antifreeze principle beyond dissolved chains.","A direct experimental test would be to measure low-temperature water reorientation near cellulose with ultrafast infrared or broadband dielectric spectroscopy; slower-than-bulk reorientation in the first hydration shell would corroborate the simulated tetrahedral arrest."],"forward_implications":["Since the q distribution is nearly unchanged from four to eight cellulose repeat units, short cellulose oligomers should show the same antifreeze character as longer chains, and the effect is confined to the glycan–water interface.","At 180 K the presence of the glycan slows the decay of C_q(t) relative to pure water, so the hydration layer retains liquid-like mobility and cannot complete the tetrahedral rearrangement needed for freezing.","Tetrahedrality of hydration water can serve as a screening metric for designing cellulose-based antifreeze additives, because it captures the molecular mechanism without simulating full ice growth.","Together with the earlier ab initio work, the results imply that tetrahedral coordination governs both ends of the process: cellulose recognition of ice planes and disruption of tetrahedral order in nearby water."],"supporting_citations":[{"why":"It is the preceding ab initio study that this work extends; it established that cellulose binds ice planes through tetrahedral coordination.","marker":"[4]"},{"why":"It supplies the coarse-grained carbohydrate bead representation that motivates the A/B/C glycan model used here.","marker":"[7]"},{"why":"It provides the molecular-dynamics precedent that hydration water near an antifreeze protein's ice-binding face has enhanced tetrahedrality.","marker":"[10]"},{"why":"It defines the tetrahedrality order parameter q used to measure water structure in this work.","marker":"[12]"},{"why":"It supplies the ice-identification routine used to remove the ice phase before computing tetrahedrality distributions.","marker":"[13]"},{"why":"It supplies the tetrahedrality autocorrelation function C_q(t) used to quantify the timescale of tetrahedral rearrangement.","marker":"[14]"},{"why":"It provides the coarse-grained mW water model, which is constructed to reproduce tetrahedral water structure and ice formation.","marker":"[17]"}],"fun_headline_variants":["Tetrahedral ordering of water is key to glycan antifreeze","Glycan antifreeze hinges on water's tetrahedral stall","Tetrahedrality metric reveals glycan ice-blocking power","Water tetrahedral disruption drives glycan antifreeze"],"cache_read_input_tokens":8448,"weakest_assumption_plain":"The load-bearing premise is that the coarse-grained bead representation of cellulose, with its hand-set hydrophilic water-bead interactions, preserves the way real cellulose disrupts the tetrahedral ordering of water; if that mapping is wrong, the q-based conclusions describe the model rather than real glycans.","fun_headline_variants_meta":{"raw":{"variants":["Tetrahedral ordering of water is key to glycan antifreeze","Glycan antifreeze hinges on water's tetrahedral stall","Tetrahedrality metric reveals glycan ice-blocking power","Water tetrahedral disruption drives glycan antifreeze"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000422,"raw_usage":{"total_tokens":2129,"prompt_tokens":868,"completion_tokens":1261,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":1193}},"tokens_in":484,"tokens_out":1261,"duration_ms":11825,"temperature":1.0,"reasoning_tokens":1193,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T04:34:22.264874+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A fully atomistic molecular dynamics simulation of cellulose in an atomistic water model cooled under the same protocol, computing the same q distribution and C_q(t), would settle the claim: if hydration water in the first shell reaches bulk-like tetrahedrality and its autocorrelation decays as fast as in pure water at 180 K, the coarse-grained result is a model artifact.","supporting_citations":[{"cited_title":"Molecular Systems Design & Engineering11, 107–117 (2026) https://doi.org/10.1039/D5ME00137D","cited_arxiv_id":null,"evidence_quote":"It is the preceding ab initio study that this work extends; it established that cellulose binds ice planes through tetrahedral coordination."},{"cited_title":"The Jour- nal of Physical Chemistry B108(4), 1414–1427 (2004) https://doi.org/10.1021/ jp0354752","cited_arxiv_id":null,"evidence_quote":"It supplies the coarse-grained carbohydrate bead representation that motivates the A/B/C glycan model used here."},{"cited_title":"Journal of the American Chemical Society130(39), 13066–13073 (2008) https://doi.org/10","cited_arxiv_id":null,"evidence_quote":"It provides the molecular-dynamics precedent that hydration water near an antifreeze protein's ice-binding face has enhanced tetrahedrality."},{"cited_title":"The Journal of Physical Chemistry B119(29), 9369–9376 (2015) https://doi.org/10.1021/ jp510289t","cited_arxiv_id":null,"evidence_quote":"It supplies the ice-identification routine used to remove the ice phase before computing tetrahedrality distributions."},{"cited_title":"Pro- ceedings of the National Academy of Sciences106(52), 22130–22134 (2009) https://doi.org/10.1073/pnas.0911094106","cited_arxiv_id":null,"evidence_quote":"It supplies the tetrahedrality autocorrelation function C_q(t) used to quantify the timescale of tetrahedral rearrangement."}],"review_version":1}