{"id":"f6d305e0-ad80-4b01-84ec-5c276d130afb","arxiv_id":"2608.13016","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Cyclic switching of membrane-bound proteins between two curvature states produces blinking domains whose size oscillates in place, even on tensionless membranes.","lead":"This simulation paper shows that proteins that bind to a membrane in two different curved shapes can make curved membrane patches blink, growing and shrinking in place. The finding matters because it predicts a new type of oscillating pattern that can appear on floppy, tensionless membranes where static curved domains would normally be unstable.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Blinking may be an artifact of the MC flip-attempt interval tau_MC=0.01tau; the paper never sweeps this unphysical parameter.","rationale":"The paper's central observation—blinking domains that persist at zero surface tension—is presented as a property of the nonequilibrium cycling model. I read the simulation protocol carefully. The membrane dynamics are integrated with a leapfrog scheme at dt=0.002 tau, while the three-state flips are attempted once per particle every tau_MC=0.01 tau with the Metropolis rule Eq. (6). The paper sweeps chemical potentials and surface tension but never varies tau_MC. This matters because the blinking mechanism is explicitly a race: convex s=2 domains must grow to a visible size and then convert to s=1 before bud formation (Sec. III B and III C). The rate of that conversion is proportional to the number of MC attempts per unit time (and to the acceptance probability). Thus the period and even the existence of BD depend on tau_MC. At gamma=0, the alternative to blinking is vesiculation (VES), so if tau_MC is raised by a factor of 10, one would expect the s=2 domains to have more time to pinch off. At smaller tau_MC, fast cycling may prevent domain growth. Since tau_MC is an unphysical simulation parameter (it does not correspond directly to a known protein binding rate), the headline 'blinking domains can form even in tensionless membranes' is, as stated, a statement about a particular simulation protocol, not about the physical system unless the result is robust to this parameter. The reader's weakest assumption concerned the memoryless, static-chemical-potential implementation; my concern is more specific and operational: the MC sweep interval is a dial that controls the very race that produces blinking, and the paper does not tell us how the phase diagram changes as this dial is turned. I therefore recommend keeping the CONDITIONAL verdict, with the explicit condition that the authors report a tau_MC sweep (e.g., 0.001-0.1 tau) at the representative BD points and at gamma=0 and gamma=1. The paper otherwise shows good internal consistency: the equilibrium reference reproduces the known HD phase, and the BD snapshots and time traces are compelling. This is a single missing robustness check, not a demonstration of error.","tokens_in":10025,"tokens_out":11086,"duration_ms":116651,"concrete_test":"Fix mu01=2, mu02=8 at gamma=0 and gamma=1. Rerun the same N-gamma-T simulations with tau_MC = 0.001, 0.003, 0.03, and 0.1 tau, keeping all other parameters fixed. Classify the resulting steady states using the paper's f3co>0.2 criterion and by visual inspection of snapshots/movies. If BD (oscillating convex domains) persists for all tau_MC values spanning two decades, the finding is robust to the flip-attempt rate. If BD appears only near 0.01 tau, the central claim is an artifact of the MC sweep interval and the verdict should be reconsidered.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim (blinking domains at gamma=1 and gamma=0) rests on a dynamical balance between the rate of protein state changes and the rate of membrane shape relaxation. In the model, the state-flip rate is set by two arbitrary simulation parameters: the MC attempt interval tau_MC=0.01tau and the acceptance rule in Eq. (6). The paper varies only the chemical potentials mu01, mu02, mu12 and the surface tension gamma; it never varies tau_MC. The BD mechanism described in Sec. III B requires that an s=2 domain grows to a finite size and then converts to s=1 before budding; this is a race between domain growth (set by membrane relaxation, diffusion, and line tension) and the s=2->s=1 flip rate. If tau_MC were larger, the s=2 domains would have more time to grow and could form buds (VES) instead of blinking; if tau_MC were smaller, the system might stay homogeneous or exhibit FD. Thus the existence and boundaries of the BD phase are controlled by an unphysical protocol parameter, and the paper provides no evidence that the chosen value represents a physical binding/unbinding rate. Because the abstract's headline result ('blinking domains can form even in tensionless membranes') is demonstrated only for tau_MC=0.01tau, the claim is conditional on this arbitrary choice.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports meshless membrane simulations of a three-state active Potts-like model in which membrane particles can be unbound (s=0) or bound in two curvature-inducing states (s=1 and s=2 with spontaneous curvatures C0=0.05 and 0.1, respectively). State changes are driven by imposed chemical potentials mu01, mu02, and mu12, so the protein cycle s=0 -> 2 -> 1 -> 0 operates out of equilibrium. At positive surface tension (gamma=1) the authors find, besides the equilibrium hexagonal domain (HD) phase, a blinking domain (BD) mode in which convex s=2 domains grow, convert to the lower-curvature s=1 state, and shrink before reforming at nearly the same location, as well as flat moving-domain (FD) and vesicle (VES) regimes. The central claim is that BD is a standing-wave-like pattern that can form even at zero surface tension, where equilibrium hexagonal domains are unstable. The paper includes phase diagrams in the (mu01, mu02) plane for gamma=1 and gamma=0, time traces of state densities, cluster-size measures, and representative snapshots and movies.","tokens_in":10345,"tokens_out":5152,"duration_ms":53997,"significance":"If the results hold, the paper identifies a new nonequilibrium membrane patterning mode: cyclic protein binding/unbinding with two curvature states produces stationary blinking domains, in contrast to equilibrium HD domains that require positive surface tension. This extends earlier lattice and off-lattice active Potts studies by coupling the state dynamics to membrane deformation and thermal fluctuations, and it offers a concrete mechanism for standing-wave-like membrane patterns. The evidence is based on direct simulation output: snapshots, density time traces, cluster measures, and multiple independent runs. The main weakness is that the dynamical phase boundaries and the robustness of BD at gamma=0 rest on an unexplored simulation protocol parameter (the MC flip interval tau_MC) and on hand-tuned classification thresholds, so the quantitative phase diagram is not yet established with the same confidence as the existence of BD at representative parameter points.","major_comments":[{"comment":"The state-flip attempt interval tau_MC = 0.01 tau is introduced in Sec. II and used in all simulations, but it is never varied. The BD mechanism described in Sec. III B relies on a race between s=2 domain growth (governed by membrane relaxation, diffusion, and line tension) and the s=2 -> s=1 flip rate. With a larger tau_MC, s=2 domains would have more time to grow and could bud into vesicles instead of blinking; with a smaller tau_MC, the system might stay homogeneous or enter the FD regime. Since tau_MC is an arbitrary protocol parameter rather than a derived physical binding/unbinding rate, the headline claim that blinking domains can form even at gamma=0 is conditional on this single choice. A sweep of tau_MC, or an explicit mapping from tau_MC to a physical rate for the chemical cycle, is needed to establish that BD is not an artifact of this choice.","section":"Sec. II, Eq. (6), and Sec. III B-C"},{"comment":"The dynamic phase boundaries are assigned using hand-tuned thresholds: BD is defined by f3co > 0.2, where three-state coexistence requires phi_k > phi_th for all states, with phi_th = 0.005 for most conditions and phi_th = 0.01-0.03 in other cases, and no statistical error bars are shown for the phase boundaries. Because these thresholds are chosen ad hoc, the location and, in narrow regions, even the existence of the BD phase could depend on the classification rule. The snapshots and time traces convincingly demonstrate BD at representative points, but the quantitative phase diagrams and the claim that BD is a 'robust' dynamic mode would be considerably strengthened by a sensitivity analysis with respect to the thresholds and by reporting classification results across a range of thresholds.","section":"Sec. III B, Figs. 3 and 8"}],"minor_comments":[{"comment":"The criterion for three-state coexistence uses the condition phi_k > phi_th, but phi_k is not explicitly defined in the text; it should be stated that these are the instantaneous densities of the three states used in the time-fraction calculation f3co.","section":"Sec. III B, after Fig. 7"},{"comment":"The manuscript states that statistical errors are calculated from three or more independent runs, but it does not report the number of runs used for each phase diagram point or show error bars on the phase boundaries; a brief statement of run counts and measurement uncertainties would improve reproducibility.","section":"Sec. II, statistical errors"},{"comment":"The y-axis label of Fig. 2(e) appears to read 'HDE0', which is likely a typographical error; it should read 'E0' or be otherwise clarified.","section":"Fig. 2(e)"},{"comment":"The expression for the entropy production rate, (mu02 - mu01 - mu12) q_f, is stated without derivation or reference; a short derivation or a citation to the relevant result would help the reader connect this quantity to the fluctuation theorem.","section":"Sec. III B, entropy production"}],"recommendation":"major_revision","confidential_remarks":"This is a technically competent simulation study from an experienced group, and the central phenomenon is interesting and plausibly correct. The main substantive gap is the missing sweep of tau_MC, which controls the rate of the chemical cycle relative to membrane relaxation; without it, the robustness of the 'tensionless blinking' claim is not fully established. The phase-threshold issue is secondary but should be addressed. I recommend major revision rather than rejection because the concerns are fixable by additional simulations and a robustness analysis, and the paper's contribution would then be solid."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: the paper reports a new dynamic membrane pattern—blinking domains—where convex domains of high-curvature proteins grow, flip to a lower-curvature state, shrink, and reform in place. The central evidence is direct: snapshots, time traces of densities, cluster-size measures, and multiple independent runs. The pattern is credible, and the tensionless-membrane case is the strongest claim, since equilibrium theory says hexagonal domains need positive tension.\n\nWhat is new: the BD mode is a standing wave in domain size without ballistic motion, and the paper makes a reasonable case that it has not appeared in earlier lattice or off-lattice active Potts models. The mechanism is transparent: cyclic chemical potentials set up a race between domain growth and state conversion, and membrane bending pins the domain location. The model is the author's established meshless membrane approach, and reusing the same parameter set as prior work helps consistency. The zero-tension result is a genuine nonequilibrium addition.\n\nThe main soft spot is the stress-test concern. The MC flip-attempt interval tau_MC=0.01tau is never varied. The BD phase is a dynamical balance between flip rate and membrane relaxation; at a slower flip rate the s=2 domains could grow into buds, and at a faster rate the system might stay homogeneous or show flat domains. The paper demonstrates BD at this one ratio but does not show the BD region is robust to tau_MC. That does not invalidate the simulation at that parameter set, but it makes the headline claim conditional on a choice the paper does not justify as physical. A sweep of tau_MC over a decade would settle it. Minor soft spots: the phase diagrams use hand-tuned thresholds (f3co>0.2, phi_th=0.005 or 0.01–0.03) and no error bars on boundaries, so the phase boundaries are approximate; and no code or data are shipped, which is typical but limits independent reproduction.\n\nThis paper is for active membrane physics, nonequilibrium soft matter, and people modeling protein binding/unbinding on deformable surfaces. It is a solid within-subfield advance, not a field reshaped. A serious referee should ask for the tau_MC sweep and clearer classification, but the core observation is likely to hold. I would take it at peer review.","headline":"Blinking domains are a genuinely new nonequilibrium membrane pattern with direct evidence, but the unswept flip-attempt interval makes the headline claim conditional.","tokens_in":10854,"tokens_out":3944,"would_cite":true,"duration_ms":38048,"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":"Cyclic binding and unbinding of curvature-inducing proteins in two states makes membrane domains blink—growing, shrinking, and reforming in place—even at zero surface tension.","keywords":["blinking domains","active Potts model","curvature-inducing proteins","nonequilibrium membrane patterns","meshless membrane simulation","spontaneous curvature","standing waves"],"falsifier":"Run the same two-state curvature cycle with state-flip rates that depend on local protein density or on instantaneous local curvature, for example a Gillespie-type binding kinetics with concentration-dependent rates, and check whether the blinking domain phase at $γ=0$ survives; if it disappears or becomes a traveling wave, the memoryless fixed-chemical-potential flip rule is essential to the reported pattern. In an experiment, fluorescently labeled curvature-inducing proteins with two switchable conformations on tensionless vesicles under an ATP-driven binding cycle should show stationary oscillating domains if the claim is right.","tokens_in":9828,"feed_emoji":"🫧","tokens_out":6070,"duration_ms":59094,"temperature":0.7,"pith_summary":"This paper establishes that a membrane driven out of equilibrium by the cyclic binding and unbinding of curvature-inducing proteins can support \"blinking domains\": patches that grow, shrink, and regrow at the same location without moving away. The proteins have two bound states with different spontaneous curvatures ($C_0=0.05$ and $0.1$), and the chemical-potential imbalance drives a cycle $s=0\\to2\\to1\\to0$. Convex domains of the higher-curvature state grow, then switch to the lower-curvature state and shrink, repeating indefinitely. The result matters because it shows a nonequilibrium standing-wave pattern on a deformable membrane that does not require surface tension, whereas equilibrium hexagonal domains are stable only under positive tension.","feed_headline":"Protein cycles make membrane domains blink","feed_subtitle":"Simulations show two-state curvature-inducing proteins drive domains to grow, shrink, and reform at the same spot.","key_machinery":"The central object is a meshless membrane model in which each membrane particle is a binding site with three states: unbound ($s=0$) and two bound states ($s=1$ and $s=2$) with spontaneous curvatures $C_0=0.05$ and $0.1$, respectively. State changes are attempted as single-particle Metropolis flips with acceptance rate $p_{\\rm acpt}=\\min[1,\\exp(\\pm(\\Delta U-μ_{\\alpha\\beta})/k_BT)]$, where $μ_{\\alpha\\beta}$ is the chemical potential between states; the drive is a fixed chemical-potential imbalance $μ_{02}-μ_{01}-μ_{12}>0$ that forces the cyclic direction $s=0\\to2\\to1\\to0$. The load-bearing identity is the relation between bound-state curvature and membrane bending: higher-curvature $s=2$ domains are stabilized by bending energy under tension, whereas switching to lower-curvature $s=1$ makes them unstable, producing the growth-then-shrink blink cycle. Bending rigidity and the saddle-splay modulus both change on binding, coupling protein distribution to local membrane shape and keeping the domain location pinned.","core_discovery":"The central claim is that cyclic protein binding and unbinding creates a new nonequilibrium pattern, blinking domains, in which the domain size oscillates in place, and this pattern survives in tensionless membranes where equilibrium convex domains are unstable. At $γ=1$, with $μ_{01}\\simeq2$ and $μ_{02}\\simeq8$, convex $s=2$ domains grow toward their stable size, then switch to $s=1$; the resulting $s=1$ domains are unstable against the unbound state and shrink, after which the cycle restarts. Membrane bending holds the domains at nearly fixed locations, so the pattern behaves like a standing wave rather than a traveling or diffusing pattern. In thermal equilibrium the same model gives hexagonally ordered convex domains only under positive surface tension, but the blinking mode appears even at $γ=0$, before domains can grow into buds or vesicles. The paper argues this makes blinking a robust pattern on deformable membranes and a minimal off-lattice analogue of standing-wave concentration dynamics.","pith_inferences":["If the blinking mechanism extends to cooperative or concentration-dependent binding, it could provide a motor-free way for cells to create oscillating protein patches at fixed membrane sites, for instance during endocytic or signaling cluster turnover.","The role of the saddle-splay change on binding is implicit in the model but not isolated; a testable extension would compare blinking with particles that change only spontaneous curvature and not $\\bar{\\kappa}$, to see whether Gaussian-curvature sensing contributes to pinning.","The two-state cycle is the minimal conformational switch; adding a third state or competing cycles in the deformable membrane may turn the standing blink into traveling waves, analogous to the $q>3$ lattice results the paper cites.","A direct experimental test would use two switchable curvature-generating domains on a tensionless giant vesicle with an externally cycled chemical fuel; the paper predicts stationary oscillating domains rather than uniform budding."],"forward_implications":["Blinking domains form at zero surface tension, so nonequilibrium protein cycling can organize membrane curvature patterns where equilibrium phase separation would fail.","Because the blinking domains stay in place, the pattern is a standing-wave-type concentration oscillation, distinct from the traveling or diffusing domains seen in earlier lattice and off-lattice active Potts models.","The blink frequency and the fraction of the $s=1$ state increase with the chemical-potential imbalance, so the oscillation period is controlled by the drive strength.","The steady state dissipates energy at rate $(μ_{02}-μ_{01}-μ_{12})q_f$, tying the pattern directly to entropy production.","At too large or too small $μ_{02}$ in tensionless membranes the cycle instead buds off vesicles, so blinking occupies a finite window of drive parameters."],"supporting_citations":[{"why":"Provides the previous off-lattice active Potts membrane model with the same parameter set and biphasic domain dynamics that this paper extends.","marker":"[51]"},{"why":"Establishes the equilibrium hexagonal convex-domain phase that is the reference state for the blinking and other nonequilibrium modes.","marker":"[42]"},{"why":"Supplies the meshless membrane model potential used throughout the simulation.","marker":"[58]"},{"why":"Defines the relation between the model bending parameter and spontaneous curvature, $C_0=C_{bd}/2\\sigma$.","marker":"[59]"},{"why":"Supplies the membrane model property that the saddle-splay modulus is proportional to bending rigidity, used when binding changes rigidity.","marker":"[60]"},{"why":"Provides the three-state active Potts lattice model whose blinking and propagating domains motivate the off-lattice implementation here.","marker":"[52]"}],"fun_headline_variants":["Protein cycles blink membrane domains at fixed spots","Blinking domains arise from protein binding cycles","Cyclic protein binding makes membrane domains blink","Binding-unbinding cycles drive membrane domain blinking","Protein binding cycles create blinking membrane domains"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the nonequilibrium drive can be represented by fixed chemical potentials acting on memoryless, single-particle state flips, so that binding and unbinding kinetics are assumed to depend only on the instantaneous local energy difference and not on protein concentration, cooperativity, or shape changes during the flip itself.","fun_headline_variants_meta":{"raw":{"variants":["Protein cycles blink membrane domains at fixed spots","Blinking domains arise from protein binding cycles","Cyclic protein binding makes membrane domains blink","Binding-unbinding cycles drive membrane domain blinking","Protein binding cycles create blinking membrane domains"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000519,"raw_usage":{"total_tokens":2475,"prompt_tokens":868,"completion_tokens":1607,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":1542}},"tokens_in":484,"tokens_out":1607,"duration_ms":11993,"temperature":1.0,"reasoning_tokens":1542,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:29:54.076017+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same two-state curvature cycle with state-flip rates that depend on local protein density or on instantaneous local curvature, for example a Gillespie-type binding kinetics with concentration-dependent rates, and check whether the blinking domain phase at $γ=0$ survives; if it disappears or becomes a traveling wave, the memoryless fixed-chemical-potential flip rule is essential to the reported pattern. In an experiment, fluorescently labeled curvature-inducing proteins with two switchable conformations on tensionless vesicles under an ATP-driven binding cycle should show stationary oscillating domains if the claim is right.","supporting_citations":[{"cited_title":"Noguchi \\ and\\ author G","cited_arxiv_id":null,"evidence_quote":"Supplies the meshless membrane model potential used throughout the simulation."},{"cited_title":"Shiba \\ and\\ author H","cited_arxiv_id":null,"evidence_quote":"Defines the relation between the model bending parameter and spontaneous curvature, $C_0=C_{bd}/2\\sigma$."},{"cited_title":"Noguchi ,\\ 10.1063/1.5113646 journal journal J.\\ Chem.\\ Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the membrane model property that the saddle-splay modulus is proportional to bending rigidity, used when binding changes rigidity."},{"cited_title":"Noguchi , author F","cited_arxiv_id":null,"evidence_quote":"Provides the three-state active Potts lattice model whose blinking and propagating domains motivate the off-lattice implementation here."}],"review_version":1}