{"id":"c03df6e2-f158-4217-8699-4f93a5d68a9d","arxiv_id":"2506.14347","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Monte Carlo simulations show arc-shaped membrane components drive a pearling-to-cylinder transition as nematic alignment increases, while saddle-shaped components stabilize necks between convex bare membrane regions.","lead":"This paper uses computer simulations to show that banana-shaped and saddle-shaped curved proteins on a vesicle surface can bend it into tubes, pearls, and narrow necks. It maps when each shape appears based on how strongly the proteins align with each other and how much volume the vesicle keeps.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The saddle-CMC neck claim is the paper's headline mechanism, yet Section 4.2.3 concedes the relevant phase diagram is metastable and poorly reproducible; equilibrium support is missing.","rationale":"The reader's weakest assumption matches the most load-bearing gap I can identify. The pearling-to-cylinder transition for arc CMCs is comparatively well supported: Figure 3 shows convergence of energy and volume for selected arc cases, and SI 8.3.1 shows the transition from an initially empty cylinder, providing reproducible-looking evidence. The saddle neck stabilization, however, is the part of the abstract that is most novel and most directly tied to the 'mechanisms' claim, and it is the part the authors themselves flag as metastable and poorly reproducible. The finite MC length (500 sweeps) and single snapshots per parameter point cannot establish that the ensemble is equilibrated; the paper's own admission that phases coexist at the same v and w is evidence of multiple basins. This does not make the paper worthless—the observed metastable necks may still be biologically relevant—but it means the central equilibrium-phase claim is not yet established. A multi-seed, longer-time, initial-condition-swap test would settle whether the concern lands. The reader's CONDITIONAL verdict remains appropriate pending that check.","tokens_in":21076,"tokens_out":4375,"duration_ms":48272,"concrete_test":"Run replicate MC simulations at the neck-phase parameters of Figure 13: ρ=0.16 and ρ=0.5, Hm=0, Dm=0.98, w=3, and reduced volumes v=0.54 and 0.61. Use at least 10 independent random seeds, starting from both a sphere and a pre-equilibrated neck, and run 10x the default 500 sweeps. Compare final morphology occupancy and energy distributions. If the neck appears from both initial conditions with equal or lower energy than coexisting oblate/prolate states, the equilibrium claim is supported; if necks only persist when pre-formed, or have higher energy than coexisting states, the phase must be reclassified as metastable/kinetic and the abstract's stabilization claim weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that saddle-shaped CMCs 'stabilize necks between the convex regions of bare membrane' is a thermodynamic statement, but the only evidence is finite-length Monte Carlo trajectories. Section 4.2.3 states that the saddle-CMC phase diagram (Figure 12) 'features a high degree of metastability, resulting in approximate phases that often coexist for the same v and w' and that 'poor reproducibility of these transitions warrants further investigation.' The same section notes that the gyration-tensor eigenvalues 'cannot reliably distinguish between phases,' yet the boundaries in Figure 12 are drawn by hand from those eigenvalues. Convergence is demonstrated in Figure 3 only for three arc-CMC cases, not for the saddle neck or invaginated states. Without run-to-run statistics, longer equilibration checks, or tests of initial-condition independence, the neck phase could be a long-lived kinetic trap rather than an equilibrium morphology. Since the abstract's distinctive claim—neck stabilization by saddle CMCs—depends on this phase, the admitted metastability is a load-bearing gap, not a peripheral caveat.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript uses Monte Carlo simulations of triangulated vesicles with a deviatoric-elasticity bending energy to study how anisotropic curved membrane components (arc- and saddle-shaped), nematic inter-particle interactions, and volume constraints determine vesicle shape. It reports a pearling-to-cylinder transition for fully arc-covered vesicles, a shift of the oblate-to-prolate transition for half arc-covered vesicles, and neck stabilization by saddle-shaped CMCs. It also analyzes nematic topological defects and compares a fully covered saddle-CMC shape with an axisymmetric calculation. The main evidence is a set of steady-state phase diagrams in the reduced-volume versus interaction-strength plane.","tokens_in":21308,"tokens_out":6343,"duration_ms":63705,"significance":"If the reported phases are true equilibrium morphologies, the work offers a useful framework for anisotropic curvature sensing and membrane remodeling and has clear biological relevance (BAR domains, neck-stabilizing proteins). The manuscript is transparent about its simulation protocol, builds on a well-defined prior theory rather than fitting phase diagrams, and includes a valuable cross-check of defect positions against a fixed-shape axisymmetric calculation. The principal significance risk is that the headline saddle-CMC neck-stabilization claim rests on a phase diagram that the authors themselves label as highly metastable and poorly reproducible, so the thermodynamic reading of the neck phase is not yet supported.","major_comments":[{"comment":"The abstract's distinctive claim that saddle-shaped CMCs 'stabilize necks between the convex regions of bare membrane' is supported only by the neck phase in Figure 12, but the text states that this phase diagram 'features a high degree of metastability, resulting in approximate phases that often coexist for the same v and w' and that the gyration-tensor heatmaps 'cannot reliably distinguish between phases.' No run-to-run statistics, initial-condition independence tests, or convergence traces are reported for the neck or invaginated saddle-CMC states; the convergence checks in Figure 3 cover only three arc-CMC cases. A neck observed in a single finite-length trajectory is not evidence for equilibrium stabilization. Please provide, for representative neck, invaginated, and oblate saddle-CMC parameters: multiple independent runs from distinct initial shapes, longer runs with convergence analysis of energy, volume, and a shape order parameter, and a hysteresis test starting separately from necked and spherical configurations at the same parameters. Without such data, the neck-stabilization claim is not established.","section":"Section 4.2.3 and Figure 12"},{"comment":"The phase boundaries in Figures 9 and 12 are drawn by hand from gyration-tensor eigenvalues, but for the saddle-CMC case the text admits that these eigenvalues 'cannot reliably distinguish between phases,' and for the arc-CMC case the separation between boomerang, dumb-bell, and prolate phases is described as 'more qualitative than quantitative.' A hand-drawn boundary based on an order parameter that cannot resolve the phases is not reproducible and does not support the quantitative phase maps. Please define a reproducible classification protocol—for example, thresholds on asphericity and prolateness combined with cluster curvature statistics, or a documented clustering algorithm on shape descriptors—and report the resulting boundaries, including a statement of how robust they are to the chosen thresholds.","section":"Section 4.2.1, Figure 9; Section 4.2.3, Figure 12"},{"comment":"The manuscript reports one representative snapshot per phase point but does not state how many independent runs were performed for each parameter set or how representative the displayed shape is. In the metastable regions of Figure 12 this matters directly: a single run cannot establish which phases 'coexist' or their relative probabilities. Please report the number of independent runs per point, the fraction of runs that ended in each phase, and the criterion for selecting the displayed snapshot.","section":"Section 3"}],"minor_comments":[{"comment":"The symbol E2 is used for both the isotropic binding energy (Eq. 3) and the nematic interaction energy (Eqs. 4-6), which confuses the two distinct mechanisms; please use different symbols such as E_iso and E_nem.","section":"Section 2.3"},{"comment":"The sentence 'uses 8 through a dot product' appears to be a typo for 'uses Eq. (8) through a dot product'; please correct it.","section":"Section 2.5"},{"comment":"The caption contains 'stady-state', which should read 'steady-state'.","section":"Figure 12 caption"},{"comment":"The mismatch tensor is first defined in Section 2.1 as M = R Cm R^{-1} - C with an explicit rotation angle omega, but Eq. (23) presents the vertex-level mismatch directly as S minus a diagonal Cm; please clarify how the rotation angle omega enters the discretized computation.","section":"Section 8.2.2 and Eq. (23)"},{"comment":"The derivation from the Frank energy (Eq. 4) to the Lebwohl-Lasher and XY forms (Eqs. 5-6) should state explicitly which gradient terms are dropped and how the one-constant approximation sets the scale of w relative to k'_G and k'_c.","section":"Section 2.3"}],"recommendation":"major_revision","confidential_remarks":"This is a serious simulation study that fits the scope of cond-mat.soft, and the authors' candid admission of metastability is to their credit, but the central biological claim in the abstract is currently ahead of the evidence. The load-bearing gap is the lack of equilibrium support for the saddle-CMC neck phase; the requested reproducibility and equilibration controls can be supplied within a revision. I do not see a circularity problem, since the phase diagrams are simulation outputs rather than fits to the target phenomena."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper deserves a serious referee, with the expectation of revision. The genuinely new content is the systematic phase mapping in the v–w plane for partially covered vesicles: the pearling-to-cylinder transition for arc-shaped CMCs as nematic strength grows, a boomerang phase between oblate and prolate, and the local nematic order Si at transitions. The topological-defect analysis is clean and consistent with Gauss–Bonnet/Poincaré–Hopf, and the comparison to the axisymmetric calculation adds credibility. The paper is also honest about its limitations, which is worth respecting.\n\nThe soft spot is load-bearing, not cosmetic. The abstract's distinctive claim—that saddle-shaped CMCs stabilize necks between convex bare-membrane regions—is an equilibrium statement, but the evidence in Section 4.2.3 is a phase diagram the authors themselves say features \"a high degree of metastability\" and \"poor reproducibility of these transitions,\" with boundaries drawn by hand from gyration eigenvalues they acknowledge \"cannot reliably distinguish between phases.\" Convergence is demonstrated for three arc-shaped cases (Figure 3), not for the neck or invaginated states. Without run-to-run statistics, longer equilibration checks, or tests of initial-condition independence, the neck phase could be a long-lived kinetic trap rather than an equilibrium morphology. The arc-shaped results are on firmer ground: the pearling-to-cylinder transition is consistent across parameters and matches earlier analytic predictions. The saddle neck claim needs more work before it can be treated as equilibrium.\n\nMinor issues: the volume modulus kv is never stated, no code or data are provided, and the phase boundaries are drawn by hand. The novelty is incremental—the Hamiltonian and Monte Carlo scheme extend the authors' own prior framework—but the new phase diagrams are genuinely new and useful. The citation pattern looks appropriate; the heavy self-citation is mostly to the original deviatoric-curvature model, which is fair.\n\nRecommendation: send it to peer review. An editor should ask for code/data, multiple initial conditions, longer runs, and a more defensible definition of the saddle-CMC phase boundaries. If the neck phase survives that scrutiny, this becomes a solid reference for the membrane-remodeling community.","headline":"Useful systematic phase maps for arc-shaped CMCs; the saddle-CMC neck claim needs equilibration evidence before it can be trusted as equilibrium.","tokens_in":21880,"tokens_out":3358,"would_cite":false,"duration_ms":37024,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["87.16.Dg","87.10.Rt"],"model":"deepseek-v4-flash","headline":"Simulations show arc-shaped membrane components pass from pearls to a smooth cylinder as nematic alignment strengthens, while saddle-shaped components stabilize the necks between bare convex vesicle regions.","keywords":["membrane curvature","nematic ordering","vesicle morphology","topological defects","anisotropic membrane inclusions","Monte Carlo simulation","deviatoric curvature","reduced volume"],"falsifier":"Re-run the same reduced-volume and interaction-strength points from several unrelated starting shapes — a random sphere, a pre-formed pearl chain, a ready cylinder — and extend the runs well beyond the 500 Monte Carlo sweeps used here; if the reported neck and invaginated phases appear only for particular starting shapes or dissolve on longer runs, those phases are kinetic transients rather than equilibrium states.","tokens_in":20883,"feed_emoji":"🫧","tokens_out":12153,"duration_ms":108861,"temperature":0.7,"pith_summary":"The paper claims that two minimal ingredients — an anisotropic intrinsic curvature printed on membrane components and a tendency for neighboring components to align — are enough to account for the variety of shapes closed membranes take. In Monte Carlo simulation, vesicles fully covered with arc-shaped components pass from a chain of pearls to a smooth cylinder as nematic interaction strengthens, and vesicles half covered with saddle-shaped components collect those components in the necks between bare convex regions, holding the neck open. The paper maps these steady-state shapes in the plane of reduced volume and interaction strength, revealing oblate, prolate, boomerang, capped, invaginated, and neck phases, and shows that the in-plane nematic field develops topological defects at high-curvature sites, with total charge matching the Poincaré-Hopf theorem. A sympathetic reader would take this as a candidate minimal mechanism by which BAR-domain proteins and other anisotropic inclusions shape organelles and whole cells.","feed_headline":"Two protein shapes can steer vesicles into pearls, tubes, and necks","feed_subtitle":"Simulations show arc-shaped components pearl then tubulate as their alignment grows, while saddles clamp membrane necks.","key_machinery":"The engine of the model is the mismatch-tensor bending energy, written per area as $(2K_1+K_2)(H-H_m)^2 - K_2(D^2 - 2DD_m\\cos 2\\omega + D_m^2)$, which penalizes both the difference between the membrane's mean and deviatoric curvatures ($H$, $D$) and the component's intrinsic values ($H_m$, $D_m$), and the in-plane rotation angle $\\omega$ between the membrane principal frame and the component's principal frame. An arc-shaped component has $H_m = D_m > 0$ (one curved, one flat direction), while a saddle-shaped component has $H_m = 0$, $D_m > 0$ and is frustrated on convex surface. To this are added a nematic (Lebwohl-Lasher / XY-type) coupling between neighboring occupied vertices and a harmonic constraint on the reduced volume $\\bar{v} = 6\\sqrt{\\pi}V/A^{3/2}$. The combined energy is minimized by Monte Carlo moves on a dynamically triangulated, self-avoiding mesh with bond flips, and steady states are classified by the eigenvalues of the gyration tensor and by the local nematic order parameter $S_i = \\frac{1}{2}(3\\cos^2\\theta - 1)$.","core_discovery":"The central claim is that coupling the mismatch between a membrane's local curvature tensor and a component's intrinsic curvature tensor to the nematic alignment of neighboring components, under a fixed reduced volume, is sufficient to produce the steady-state phase repertoire reported here. For full coverage with arc-shaped components (nonzero intrinsic mean curvature $H_m$, zero deviator), no nematic interaction yields pearling — a string of near-spheres joined by thin necks in which local nematic order vanishes — and raising the interaction strength $w$ converts the pearls into smooth cylinders whose radius is set by the component curvature. For full coverage with saddle-shaped components ($H_m = 0$, nonzero deviator $D_m$), stronger $w$ first flattens the vesicle and then gives a bow-tie or protrusive shape with a $+1$ defect at the tip and two $+1/2$ defects at its base. On half-covered vesicles, arc components produce oblate, prolate, boomerang, mixed, and capped phases that push the oblate-prolate boundary to lower reduced volume than bare membranes, an entropy effect from components clustering on the oblate rim; saddle components instead stabilize neck phases in which the components aggregate in the negatively curved necks between bare convex caps, persisting even with no volume constraint. The topological defect charge sums to $+2$ throughout, and defects sit where Gaussian curvature is largest.","pith_inferences":["A testable prediction follows from the boomerang phase: anisotropic components should generically stabilize an asymmetric intermediate along the oblate-prolate path in experiments on giant unilamellar vesicles with adsorbed anisotropic proteins, something isotropic inclusions would not produce.","The paper's reported metastability leaves open that the neck and invaginated phases are kinetically trapped; a simulated-annealing extension that slowly ramps reduced volume or interaction strength could reveal whether the phases are separated by a genuine free-energy barrier.","The framework implies a design rule: a protein's intrinsic curvature signature, through the ratio $H_m/D_m$, should predict its morphological role — arc-like components tubulate while saddle-like components clamp necks — a classification that could be tested against measured intrinsic curvatures of BAR and F-BAR domain proteins."],"forward_implications":["Weak nematic coupling among arc-shaped components leaves a fully covered vesicle pearled, so initial curvature sensing requires only weak interactions, while stronger alignment drives full tubulation.","Saddle-shaped components stabilize the necks between bare convex regions at concentrations below 10% and even without volume constraints, offering a minimal mechanism for localizing fission sites and connecting tubules.","Because arc-shaped components shift the oblate-to-prolate transition to lower reduced volume purely through mixing entropy, anisotropic components can reshape a vesicle even with no direct component-component interaction.","The pearl radius and the neck width are set by the component's intrinsic mean and deviatoric curvatures, so the steady-state shape directly encodes the geometry of the components that made it."],"supporting_citations":[{"why":"supplies the deviatoric-curvature elasticity whose energy expression (Eq. 2) drives the pearling and neck physics.","marker":"[34]"},{"why":"provides the mesoscale curvature-computation method for triangulated membranes that the simulations adapt.","marker":"[29]"},{"why":"the CMC-CMC interaction model the paper adopts for neighboring anisotropic inclusions.","marker":"[39]"},{"why":"supplies the shape-operator estimation on the mesh used to evaluate the anisotropic bending energy each Monte Carlo step.","marker":"[44]"},{"why":"the triangulated, self-avoiding Monte Carlo membrane code the simulations are run with.","marker":"[41]"},{"why":"the bare-vesicle phase diagram that serves as the no-CMC baseline for the oblate-prolate comparison.","marker":"[47]"},{"why":"the axisymmetric nematic-shell calculation the paper compares against to confirm the saddle-CMC steady state and defect positions.","marker":"[53]"},{"why":"the topological-defect framework used to interpret where defects sit and why their charges sum to two.","marker":"[50]"}],"fun_headline_variants":["Arc proteins pearl vesicles, saddles pinch necks","Curved proteins shape vesicles: pearls, tubes, necks","How arc and saddle proteins drive vesicle shape phases","Nematic alignment turns vesicle pearls into cylinders","Curvature-nematic coupling explains vesicle shapes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the finite Monte Carlo runs reach steady states that represent equilibrium morphologies, when the paper itself reports high metastability and poor reproducibility for the saddle-shaped-CMC transitions, so the neck and invaginated phases could be kinetic artifacts rather than true equilibrium shapes.","fun_headline_variants_meta":{"raw":{"variants":["Arc proteins pearl vesicles, saddles pinch necks","Curved proteins shape vesicles: pearls, tubes, necks","How arc and saddle proteins drive vesicle shape phases","Nematic alignment turns vesicle pearls into cylinders","Curvature-nematic coupling explains vesicle shapes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1425,"prompt_tokens":1016,"completion_tokens":409,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":632,"completion_tokens_details":{"reasoning_tokens":334}},"tokens_in":632,"tokens_out":409,"duration_ms":4842,"temperature":1.0,"reasoning_tokens":334,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:17:47.117956+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the same reduced-volume and interaction-strength points from several unrelated starting shapes — a random sphere, a pre-formed pearl chain, a ready cylinder — and extend the runs well beyond the 500 Monte Carlo sweeps used here; if the reported neck and invaginated phases appear only for particular starting shapes or dissolve on longer runs, those phases are kinetic transients rather than equilibrium states.","supporting_citations":[{"cited_title":"On the role of membrane anisotropy in the beading transition of undulated tubular membrane structures","cited_arxiv_id":null,"evidence_quote":"supplies the deviatoric-curvature elasticity whose energy expression (Eq. 2) drives the pearling and neck physics."},{"cited_title":"Mesoscale computational stud- ies of membrane bilayer remodeling by curvature-inducing proteins","cited_arxiv_id":null,"evidence_quote":"provides the mesoscale curvature-computation method for triangulated membranes that the simulations adapt."},{"cited_title":"A review of mechanics- based mesoscopic membrane remodeling methods: capturing both the physics and the chemical diversity","cited_arxiv_id":null,"evidence_quote":"the CMC-CMC interaction model the paper adopts for neighboring anisotropic inclusions."},{"cited_title":"Monte Carlo simulations of fluid vesicles with in-plane orientational ordering","cited_arxiv_id":null,"evidence_quote":"supplies the shape-operator estimation on the mesh used to evaluate the anisotropic bending energy each Monte Carlo step."},{"cited_title":"Theoretical study of vesicle shapes driven by coupling curved proteins and active cytoskeletal forces","cited_arxiv_id":null,"evidence_quote":"the triangulated, self-avoiding Monte Carlo membrane code the simulations are run with."},{"cited_title":"Shape transformations of vesicles: Phase diagram for spontaneous-curvature and bilayer-coupling models","cited_arxiv_id":null,"evidence_quote":"the bare-vesicle phase diagram that serves as the no-CMC baseline for the oblate-prolate comparison."},{"cited_title":"Coupling of nematic in-plane orientational ordering and equilibrium shapes of closed flexible nematic shells","cited_arxiv_id":null,"evidence_quote":"the axisymmetric nematic-shell calculation the paper compares against to confirm the saddle-CMC steady state and defect positions."},{"cited_title":"Effective topological charge cancelation mechanism","cited_arxiv_id":null,"evidence_quote":"the topological-defect framework used to interpret where defects sit and why their charges sum to two."}],"review_version":1}