{"id":"10eb47a2-2eaf-41f3-ba88-51e06707ad79","arxiv_id":"2507.17434","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Vesicle deformability suppresses full engulfment when the vesicle size is comparable to the bendo-capillary length, while geometry alone dictates wrapping for larger vesicles.","lead":"Researchers measured how small vesicles get wrapped and swallowed by larger ones, using microscope images and computer models. They found that whether full wrapping happens depends mainly on geometry for larger vesicles, but for smaller ones the vesicle's deformability can block complete engulfment.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equilibrium assumption untested for the deformability-suppression claim: the single R/L=2.9 partially wrapped pair may be kinetically trapped rather than a stable thermodynamic state.","rationale":"The reader's weakest assumption is exactly the load-bearing premise: energy-minimized simulations predict the measured morphologies, implying experimental states sit at global free-energy minima with fixed volume and area. My reading sharpens this into a concrete kinetic-trapping concern targeted at the deformability-suppression claim. The paper's own evidence for that claim is thin: one partially wrapped outlier at R/L=2.9 versus a fully wrapped point at R/L=7.1, with no reported error bars or replicate counts, and the time-lapse data show only forward wrapping. Since the neck that forms near full wrapping is a known topological barrier, the partial state could be metastable. Energy curves in Supporting Fig. S4 may already address this, but the paper does not report starting-configuration dependence or barrier-crossing tests, so the concern remains live. This is an addressable, falsifiable issue rather than a fatal flaw; it warrants the same CONDITIONAL verdict the reader gave. I did not identify a stronger objection: the geometry-regime analysis (Eq. 7) is internally consistent, the L_exp calibration via curvature matching is reasonable, and the azo-PC area-modulation experiments, while relying on a 3% area estimate, are supportive rather than load-bearing for the deformability claim.","tokens_in":18129,"tokens_out":2994,"duration_ms":37169,"concrete_test":"Re-run the Surface Evolver minimization for the exact dimensionless parameters of the R/L=2.9 partially wrapped pair, but start from an initially fully wrapped configuration (small vesicle fully enclosed, neck treated as in Supporting Fig. S6). If the global minimum is partially wrapped, the system should spontaneously unwrap to that state; if it remains fully wrapped, the landscape has a barrier and kinetic trapping is a viable alternative explanation. Complement this with an experimental reversibility check: after observing full wrapping, reduce adhesion strength (e.g., by diluting polyacrylamide or changing temperature) and monitor whether the vesicle unwraps on experimental timescales. If it does not, the equilibrium assumption underlying Fig. 6 is violated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that small-vesicle deformability inhibits full engulfment when R_small/L≈1 rests on the equilibrium phase boundaries in Fig. 6 and on two experimental points at νγ=0.93 (R/L=7.1 fully wrapped vs R/L=2.9 partially wrapped). The paper does not establish that the R/L=2.9 pair is at a thermodynamic minimum rather than a metastable state. Surface Evolver is a local gradient-descent minimizer from an initial guess, and the time series in Fig. 4B only shows wrapping progressing forward; no spontaneous unwrapping is shown. If the partially wrapped state is trapped behind the catenoidal-neck barrier, the observed suppression is a kinetic effect, not the equilibrium deformability effect claimed. The geometry-dominated result (Eq. 7) would survive, but the headline claim that deformability 'dictates' the shallow-to-full transition would be undercut. Moreover, the experimental support for Fig. 6 is essentially one partially wrapped pair compared with one fully wrapped pair at a different R/L; if those pairs also differ in ν_small, φ, or ν_large (the text says 'nearly identical' but does not tabulate all four parameters), the attribution to R/L/deformability is confounded. Error bars and replicate counts are not provided for these two points.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript combines confocal fluorescence imaging of GUV–GUV engulfment with Surface Evolver continuum simulations to map wrapping morphologies as functions of volume ratio, reduced volumes, and the ratio of small-vesicle size to the bendo-capillary length. The authors derive an analytic geometric condition for full engulfment in the R_small/L >> 1 regime (Eq. 7), calibrate L_exp by matching peak mean curvature, and construct state diagrams for endocytic and exocytic geometries. They report a geometry-dominated regime and a deformability-dominated regime near R_small/L ~ 1, where lower small-vesicle reduced volume suppresses complete wrapping despite sufficient excess area.","tokens_in":18305,"tokens_out":12233,"duration_ms":115729,"significance":"If the central claim holds, this is a valuable direct experimental test of long-standing theoretical predictions that soft-object deformability inhibits full engulfment, with implications for endocytosis, viral entry, and drug delivery. The paper's strengths include quantitative 3D morphometry, the clean analytic formula in Eq. 7, treatment of both endo- and exocytic geometries, and a photo-switchable area-expansion control. The main evidential weakness is that the headline deformability-suppression result rests on two experimental points and an untested equilibrium assumption, so the definitive status of that claim depends on additional experiments and analysis.","major_comments":[{"comment":"The central claim that deformability suppresses full wrapping when R_small/L ≈ 1 assumes each observed pair is at the global minimum of the Helfrich free energy in Eq. 8. Surface Evolver is a local gradient-descent minimizer from an initial guess, and the time series in Fig. 4B only shows progressive wrapping; no spontaneous unwrapping or repeated cycling is reported. A partially wrapped pair at R/L = 2.9 could be kinetically trapped behind the catenoidal-neck barrier rather than representing the equilibrium phase boundary in Fig. 6. Please test reversibility directly, for example by lowering the polymer concentration or switching azo-PC back from cis to trans and showing that partially and fully wrapped states interconvert, and/or by computing and reporting energy barriers from near-full initial configurations. Without such a test, the geometry-dominated result (Eq. 7, Fig. 5) remains secure, but the deformability-suppression claim is not fully supported.","section":"Modelling the vesicle-vesicle engulfment; Fig. 4B; Fig. 6"},{"comment":"The experimental evidence for Fig. 6 consists of two data points at νγ = 0.93 with R/L = 7.1 (fully wrapped) and R/L = 2.9 (partially wrapped). This single comparison carries the headline claim, but the manuscript does not tabulate all four dimensionless parameters for these two pairs, does not provide replicate counts or experimental uncertainties, and does not propagate the ±0.10 µm uncertainty in L_exp to R/L. If the two pairs differ appreciably in φ, ν_small, or ν_large, or if the R/L = 2.9 point has an uncertainty reaching the phase boundary, the attribution to R/L and deformability is confounded. Please provide a table with exact values and errors for all parameters of these pairs, and ideally add more data points in the R/L ≈ 1–5 range at fixed ν_small, with error bars on wrapping fractions as well.","section":"Fig. 6 and the 'nearly identical parameters' comparison after Fig. 5"}],"minor_comments":[{"comment":"The wrapping fractions in Fig. 3C are shown without error bars or replicate counts; please state the number of independent vesicle pairs per condition and the measurement uncertainty, or acknowledge the absence of replicates in the figure caption.","section":"Fig. 3C"},{"comment":"As rendered in the arXiv text, Eq. 2 appears to define ν_i as the cube root of (4πV_i/A_i^{3/2}), which is not the standard reduced volume and would be inconsistent with Eq. 7. If the intended expression is ν_i = 3√(4π) V_i / A_i^{3/2}, please typeset it unambiguously.","section":"Eq. 2"},{"comment":"The availability statement says only that data are present in the paper and/or Supporting Information; please state whether Surface Evolver input files and analysis scripts are available, as this would strengthen reproducibility.","section":"Data and materials availability"},{"comment":"The estimated ~3% membrane area increase is inferred rather than measured directly; please clarify whether this is an upper bound and how the uncertainty affects the inferred ν_large values in Fig. 5B.","section":"Azo-PC area expansion"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good fit for the journal and the modeling framework is sound. My main concern is overinterpretation of a two-point experimental comparison in Fig. 6 and the untested equilibrium assumption; both are addressable with additional experiments or a clearly framed barrier analysis. I see no issues with citation practice or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a good paper, worth engaging. The experiments are careful, and the combination of 3D confocal reconstruction, Surface Evolver energy minimization, and a simple geometric transition condition (Eq. 7) gives a coherent picture of when GUVs fully engulf smaller GUVs under depletion adhesion. The curvature-matching calibration of L_exp from peak mean curvature is a genuinely useful idea, and the resulting L_exp = 0.61 ± 0.10 µm is credible. The state diagrams in Figs. 4-6 organize the data well.\n\nWhat's new: the experimental demonstration that for R_small/L >> 1 the full-wrapping transition is set by the effective reduced volume νγ (Eq. 7), and that the data collapse onto that boundary despite a wide range of R_small/L from ~1.9 to 19. The azo-PC area-expansion control experiments are a nice touch, even if the estimated 3% area change is rough. The exocytic vs endocytic comparison is also valuable.\n\nSoft spots: the central claim about deformability suppressing full engulfment near R_small/L ~ 1 rests on exactly two experimental points, both at νγ = 0.93, one fully wrapped at R/L = 7.1 and one partially wrapped at R/L = 2.9. That's a thin basis for the headline, and the text frankly says the outlier is 'close to the boundary' of the geometric regime. The stress-test worry that the partially wrapped state could be kinetically trapped is real but not fatal in my reading: the time series in Fig. 4B shows the system progressing to full wrapping, and the simulations reproduce the shape sequence, so the equilibrium assumption isn't obviously violated. Still, the authors don't test reversibility, and one pair is not enough to pin the deformability effect. I'd ask for more data points, error bars on the wrapping fractions, and a stability check (e.g., starting simulation from the fully wrapped state and seeing if it unwraps) before I'd take the suppression claim as established.\n\nThe calibration of L_exp by matching peak curvature to the same model used to predict the phase boundaries is a potential circularity, but it's really just parameter fitting; the wrapping fractions and state diagrams are independent outputs, so I don't see it as a flaw. The gravity explanation for large-vesicle shape mismatches is plausible but not tested directly.\n\nBottom line: the geometric transition (Eq. 7) is solid and likely to be cited; the deformability suppression is promising but under-supported. This deserves peer review and will probably survive with revisions. I'd bring it to our group meeting and cite it.","headline":"Careful experiment/simulation study with a credible geometric wrapping transition; the deformability-suppression headline is promising but rests on two points, so treat as conditional.","tokens_in":18905,"tokens_out":2924,"would_cite":true,"duration_ms":29627,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A small vesicle's deformability can suppress its complete engulfment by a larger membrane even when enough membrane area exists, according to experiments and continuum simulations.","keywords":["membrane wrapping","giant unilamellar vesicles","depletion adhesion","bendo-capillary length","reduced volume","Helfrich free energy","endocytic engulfment","vesicle deformability"],"falsifier":"Take a vesicle pair with $\\nu_{\\rm small} \\approx 0.8$ and $R_{\\rm small}/L \\approx 1$ that lies just below the predicted full-wrapping boundary, slowly raise the adhesion strength until it wraps, then slowly lower it again; if the pair stays fully wrapped below the predicted boundary, or if partial wrapping persists where the boundary predicts full wrapping, the observed states are kinetically trapped rather than equilibrium ones.","tokens_in":17867,"feed_emoji":"🫧","tokens_out":10410,"duration_ms":91248,"temperature":0.7,"pith_summary":"This paper asks what controls whether a small lipid vesicle is swallowed by a larger one, and shows that the answer depends on whether the small vesicle is much larger or about the same size as the bendo-capillary length, the scale at which membrane bending and adhesion energies balance. In the geometry-dominated regime, full engulfment is decided by whether the large vesicle has enough excess membrane area, captured by the effective reduced volume $\\nu_\\gamma$ from Eq.~7. Near $R_{\\rm small}/L \\approx 1$, that geometric criterion is no longer enough: a more deformable small vesicle, one with lower reduced volume $\\nu_{\\rm small}$, needs stronger adhesion or a larger size ratio to become fully wrapped. The support comes from matching three-dimensional confocal reconstructions of real giant unilamellar vesicle pairs to energy-minimized continuum simulations, with the bendo-capillary length calibrated by curvature matching. If the paper is right, the same energy balance explains both endocytic and exocytic uptake and gives a unified framework for soft cargo internalization in synthetic and biological settings.","feed_headline":"Floppy vesicles resist full engulfment even with excess membrane area","feed_subtitle":"A single length-scale ratio decides whether geometry or softness sets the uptake outcome.","key_machinery":"The central object is the bendo-capillary length $L = \\sqrt{\\kappa/w}$, the scale at which membrane bending and adhesion energies balance, together with the reduced volume $\\nu_i$ that measures how much excess membrane area a vesicle has and hence how deformable it is. The study organizes the problem around the ratio $R_{\\rm small}/L$ and derives an effective reduced volume $\\nu_\\gamma$ (Eq.~7) for the large vesicle once it has wrapped the small one; the condition $\\nu_\\gamma \\le 1$ is the geometric criterion for full engulfment. Energy-minimized continuum simulations of the Helfrich free energy (Eq.~8), with constant membrane area and volume enforced by Lagrange multipliers, supply the equilibrium morphologies and phase boundaries that the experiments are compared against.","core_discovery":"The paper establishes a two-regime picture of vesicle-vesicle engulfment. When $R_{\\rm small}/L \\gg 1$, adhesion dominates over bending and the partial-to-full wrapping transition is governed by geometry alone: full engulfment is possible precisely when the effective reduced volume of the large vesicle after wrapping, $\\nu_\\gamma = (1+\\phi)/(\\nu_{\\rm large}^{-2/3} - \\phi^{2/3}\\nu_{\\rm small}^{-2/3})^{3/2}$, is at most 1. When $R_{\\rm small}/L \\approx 1$, full engulfment can be suppressed even if $\\nu_\\gamma < 1$, and the suppression grows as the small vesicle becomes more deformable: the transition curves require increasing adhesion strength, equivalently larger $R_{\\rm small}/L$, as $\\nu_{\\rm small}$ drops from 0.99 to 0.80. Experiments and simulations agree on wrapping fractions and morphologies across volume ratios from 0.001 to 0.8 and size ratios from 1.6 to 19, and the measured bendo-capillary length is $L_{\\rm exp} = 0.61 \\pm 0.10\\,\\mu\\rm m$. The paper therefore claims that deformability is not a minor correction but a decisive control parameter in the crossover regime, and that the same energetic balance explains both endocytic and exocytic engulfment.","pith_inferences":["An untested consequence of the equilibrium picture is reversibility: if adhesion is slowly removed, a fully wrapped vesicle should unwrap along the same path. The paper only shows progressive wrapping, so a direct unwrapping experiment would show whether the phase boundaries are thermodynamic or partly kinetic.","The single outlier near the geometric boundary suggests that bending energy can inhibit full engulfment even when Eq.~7 says it is possible; mapping the transition systematically at intermediate $R_{\\rm small}/L$, roughly 2 to 10, would quantify how sharp the crossover really is.","The geometric criterion $\\nu_\\gamma \\le 1$ treats the small vesicle through its volume and area only, so extending the framework to non-spherical or multi-domain cargo would require an effective shape parameter beyond the reduced volume."],"forward_implications":["In the geometry-dominated regime, the outcome of engulfment can be predicted from three measured quantities, the volume ratio and the two reduced volumes, without needing a precise value of the adhesion strength or bending rigidity.","Near $R_{\\rm small}/L \\approx 1$, modest changes in the deformability of the cargo vesicle shift the adhesion strength required for full uptake, so softness is a practical control parameter for engulfment.","Increasing the excess membrane area of the large vesicle, for example by light-triggered area expansion, drives a pair from shallow to deep to full wrapping and offers an external switch for the process.","Because the same energy balance reproduces both endocytic and exocytic morphologies, the framework should transfer to uptake of soft carriers in synthetic-cell and drug-delivery contexts.","The cargo vesicle changes from oblate to prolate as wrapping deepens, so wrapping fraction and cargo shape co-evolve instead of the cargo staying rigid."],"supporting_citations":[{"why":"Supplies the experimental and theoretical setup for depletion-induced engulfment of particles by GUVs that this study adapts to vesicle-vesicle pairs.","marker":"[14]"},{"why":"Provides the theoretical prediction that wrapping transitions for nanoparticles depend on geometry and that bending relaxation near $\\nu_\\gamma=1$ lowers the required adhesion.","marker":"[19]"},{"why":"Predicts that deformability of soft particles inhibits full engulfment, the effect this paper tests experimentally and locates in the $R_{\\rm small}/L \\approx 1$ regime.","marker":"[21]"},{"why":"Models membrane interactions of non-spherical elastic particles and supports the deformability-dependent shape remodeling seen in the simulations.","marker":"[24]"},{"why":"Derives the mutual remodeling of interacting droplets and vesicles and underlies the geometric excess-area criterion expressed in Eq.~7.","marker":"[25]"},{"why":"Gives the depletion interaction expression $E_{\\rm ad}=wA_c$ used to define the adhesion strength.","marker":"[31]"},{"why":"Provides the equilibrium shape phase diagram for vesicles with given reduced volume, used to interpret the free and wrapped shapes of the vesicles.","marker":"[34]"},{"why":"Explains the gravity-induced shape deviations observed for the largest vesicles, allowing those cases to be set aside from the energy-minimization comparison.","marker":"[41]"}],"fun_headline_variants":["Deformability, not area, decides vesicle engulfment fate","Soft vesicles stall full wrapping when size matches bend length","Geometry vs softness: one length scale sets engulfment outcome","Vesicle uptake: deformability suppresses full engulfment","Beyond geometry: vesicle softness gates full wrapping"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison rests on the assumption that each measured vesicle pair has relaxed to the global minimum of the Helfrich free energy at fixed volume and area, so the computed equilibrium phase boundaries describe the experimental morphologies; the paper does not demonstrate that wrapping is reversible or that the states are not kinetically trapped by the narrow neck that forms near full wrapping.","fun_headline_variants_meta":{"raw":{"variants":["Deformability, not area, decides vesicle engulfment fate","Soft vesicles stall full wrapping when size matches bend length","Geometry vs softness: one length scale sets engulfment outcome","Vesicle uptake: deformability suppresses full engulfment","Beyond geometry: vesicle softness gates full wrapping"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000172,"raw_usage":{"total_tokens":1322,"prompt_tokens":1039,"completion_tokens":283,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":655,"completion_tokens_details":{"reasoning_tokens":201}},"tokens_in":655,"tokens_out":283,"duration_ms":3833,"temperature":1.0,"reasoning_tokens":201,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:48:08.717289+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a vesicle pair with $\\nu_{\\rm small} \\approx 0.8$ and $R_{\\rm small}/L \\approx 1$ that lies just below the predicted full-wrapping boundary, slowly raise the adhesion strength until it wraps, then slowly lower it again; if the pair stays fully wrapped below the predicted boundary, or if partial wrapping persists where the boundary predicts full wrapping, the observed states are kinetically trapped rather than equilibrium ones.","supporting_citations":[{"cited_title":"van der Ham, J","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental and theoretical setup for depletion-induced engulfment of particles by GUVs that this study adapts to vesicle-vesicle pairs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the theoretical prediction that wrapping transitions for nanoparticles depend on geometry and that bending relaxation near $\\nu_\\gamma=1$ lowers the required adhesion."},{"cited_title":"Yi and H","cited_arxiv_id":null,"evidence_quote":"Predicts that deformability of soft particles inhibits full engulfment, the effect this paper tests experimentally and locates in the $R_{\\rm small}/L \\approx 1$ regime."},{"cited_title":"Midya, T","cited_arxiv_id":null,"evidence_quote":"Models membrane interactions of non-spherical elastic particles and supports the deformability-dependent shape remodeling seen in the simulations."},{"cited_title":"Satarifard and R","cited_arxiv_id":null,"evidence_quote":"Derives the mutual remodeling of interacting droplets and vesicles and underlies the geometric excess-area criterion expressed in Eq.~7."},{"cited_title":"Asakura and F","cited_arxiv_id":null,"evidence_quote":"Gives the depletion interaction expression $E_{\\rm ad}=wA_c$ used to define the adhesion strength."},{"cited_title":"Seifert, K","cited_arxiv_id":null,"evidence_quote":"Provides the equilibrium shape phase diagram for vesicles with given reduced volume, used to interpret the free and wrapped shapes of the vesicles."},{"cited_title":"Kraus, U","cited_arxiv_id":null,"evidence_quote":"Explains the gravity-induced shape deviations observed for the largest vesicles, allowing those cases to be set aside from the energy-minimization comparison."}],"review_version":1}