{"id":"a37ec7d4-a38a-4f53-aedd-9719d1ab6f0c","arxiv_id":"2502.06582","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A model shows that a proton gradient across a spherical vesicle with pH-sensitive bending can deform it into a crista-like shape with flat poles and a curved equator.","lead":"This paper models a mitochondrial crista as a spherical vesicle whose membrane bends in response to a proton gradient flowing from the poles to the equator. The model predicts that such a gradient alone can deform the sphere into a crista-like shape with flat protein-rich zones and a highly curved ATP synthase zone, and it maps which membrane parameters give a well-functioning crista.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fitted proton profile in Eq. (2) is inconsistent with the cited experimental pH values and has the wrong sign, so the driving force in Eq. (17) may be reversed.","rationale":"The central claim requires that the prescribed proton concentration h(θ) drives the deformation in the correct direction. The paper's own text says h decreases from RC (θ0) to ATP-S (π/2), consistent with the cited pH measurements (6.4 near RC, 7.1 near ATP-S). However, Eq. (2) with the stated values of α0 and α1 gives the opposite behavior: h is negative near RC and positive near ATP-S, and the fit does not match either experimental pH value. Because the active stress in the shape equation enters through δξ=βh, this sign error reverses the active force and can invert the predicted morphology and shift the optimal parameter region. This is an internal inconsistency, not merely a parameter chosen by hand, so it is more elementary than the reader's concern about deformation feedback. Since the error is correctable by re-fitting Eq. (2) and recomputing the shapes, a conditional verdict remains appropriate; the condition should explicitly require that correction. I therefore keep the reader's CONDITIONAL verdict unchanged, but for a different and more specific reason.","tokens_in":15226,"tokens_out":12647,"duration_ms":107495,"concrete_test":"Re-fit α0 and α1 in Eq. (2) directly to the experimental values h(θ0)=10−6.4−10−7 and h(π/2)=10−7.1−10−7, using θ0≈0.1 (or the actual r/R). Then recompute u(θ) from Eqs. (17)–(26), the score map in Fig. 5, and the shapes in Figs. 4 and 6. If the optimal zone and the flat/curved morphology disappear or invert, the central claim depends on an incorrect proton profile.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (2) with the caption values α0=−6.0×10−8 mol L−1 and α1=7.9×10−8 mol L−1 does not reproduce the two experimental pH values it is claimed to fit. At the ATP-synthase location θ=π/2, h(π/2)=α1=7.9×10−8, giving pH=−log10(h0+h)≈6.75 for h0=10−7, not the cited 7.1 (ref 29). At the RC location θ0≈0.1, h≈−1.0×10−7, making h0+h≈0 and pH undefined. The text states that h decreases monotonically from RC to ATP-S; with Eq. (2) that requires α0>0, not α0<0. Since the active stress enters through δξ=βh in Eq. (17), this sign/calibration error reverses the direction of the active force. The predicted flat-pole/curved-equator morphology and the optimal region (A=6, ξ=−1) in Fig. 5 are therefore not established for a proton gradient flowing from RC to ATP-S; they may be artifacts of a reversed gradient.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models a mitochondrion crista as a spherical vesicle whose membrane has a pH-dependent spontaneous curvature. A proton concentration field h(θ) is taken to be purely diffusive between the respiratory complexes at the poles and ATP synthase at the equator, and the resulting linearized Helfrich shape equation is solved in the small-deformation regime. The authors define a scalar functionality score S based on geometric features (equator size, proton path length, flatness near the RC) and map the (A, ξ) parameter space to identify a 'well-functioning crista' region, with an optimum at (A=6, ξ=−1).","tokens_in":15377,"tokens_out":4749,"duration_ms":94808,"significance":"The manuscript provides a detailed derivation of the linearized shape equation for a closed spherical vesicle with a position-dependent spontaneous curvature, including explicit treatments of the degenerate cases A=0 and A=2 and error bounds for the truncation of the kernel functions. This is a useful methodological contribution for active-membrane models. If the pH-gradient calibration and the sign of the driving force were correct, the model would offer a minimal mechanism by which a proton gradient alone could produce inhomogeneous curvature with flat polar regions and a curved equator. However, the central quantitative claims are currently compromised by an internal inconsistency in the proton-field fit, and the 'well-functioning crista' conclusion rests on an ad hoc score. The paper is therefore promising but not yet established.","major_comments":[{"comment":"The fitted proton profile is internally inconsistent with the stated experimental pH values and has the wrong sign. With α0=−6.0×10−8 mol L−1 and α1=7.9×10−8 mol L−1, Eq. (2) gives h(π/2)=α1=7.9×10−8, so with h0=10−7 mol L−1 the pH at the ATP-S location is −log10(1.79×10−7)≈6.75, not 7.1 as claimed; at the RC location θ0≈0.1, h≈−1.0×10−7, making h0+h≈0 and pH undefined. Moreover, since dh/dθ=−α0/sinθ, the stated negative α0 makes h increase from RC to ATP-S, contradicting the text's statement that h decreases monotonically from RC to ATP-S and the pH gradient direction cited from Ref. 29. Because the active stress enters through δξ=βh in Eq. (17), this sign error reverses the direction of the driving force. The predicted flat-pole/curved-equator morphology and the optimum in Fig. 5 are therefore not established for a proton gradient flowing from RC to ATP-S.","section":"Section II, Eq. (2) and Fig. 2"},{"comment":"The functionality score S = S1 − S2 − S3 uses equal weights of 1 for the three geometric contributions, with no biophysical justification for these weights or for the chosen normalizations, and no sensitivity analysis with respect to the weighting. The paper's central claim that shapes in the identified parameter region 'correspond to a well-functioning crista' depends entirely on this score. Since S1, S2, and S3 are precisely the geometric features the model is designed to produce, the maximization of S in parameter space is partly self-referential: the score rewards the exact outputs of the mechanism being tested. The authors should either derive the score from an independent biophysical argument or demonstrate robustness of the optimal region to reasonable weight variations.","section":"Section IV, Eq. (34)"},{"comment":"The coupling parameter β = R C1^0 is chosen by hand so that the variation of the natural curvature is 'of the order of 10% of the curvature of the system' (caption of Fig. 4). No experimental or simulation constraint is given for β, and the paper does not test whether the predicted shapes persist for smaller β. Since the deformation is linearly proportional to β through Eq. (17) (and the particular solution in Eq. (22)), the existence of the observed flat-pole/curved-equator morphology is conditional on this arbitrary choice. A sensitivity analysis over β, or a measurement-based estimate, is required to support the claim that the proton gradient suffices to deform the vesicle.","section":"Section III, after Eq. (19)"},{"comment":"The abstract states that the phase diagram is 'compared to experimental measurements,' but the paper does not make a quantitative comparison to experimental crista shapes or pH measurements. The only experimental inputs are the pH values used (incorrectly, as noted above) and parameter values cited for R, κ, and σ0. The Discussion's statement that the optimum gives C0≈−0.01 nm 'in agreement with recent simulations' is not a comparison to measurements. To support the claimed comparison, the authors should plot the model's deformed shapes against experimental images or quantitative morphometric data for cristae, or at minimum against the in vitro vesicle data of Refs. 14 and 15.","section":"Abstract and Discussion"}],"minor_comments":[{"comment":"The symbol h0 is written as '10 × 10−7 mol L−1', which is ambiguous: it could mean 10−7 or 10−6 mol L−1. Since the pH calibration depends on h0, this must be stated unambiguously.","section":"Section II"},{"comment":"The angle γ is used in the boundary condition but is defined only in the caption of Fig. 1. Define γ explicitly in the main text when it first appears.","section":"Section III, Eq. (28)"},{"comment":"The volume constraint δV = ∫ sin(θ)dθ u(θ) = 0 is missing the factor 2πR^3 that appears in the full volume integral. While this factor cancels in the constraint, the equation as written is dimensionally inconsistent with δV.","section":"Section IV, Eq. (31)"},{"comment":"The caption states that h(θ) decreases monotonically from RC to ATP-S, but for the given α0 and α1, dh/dθ>0, so it increases. This is the same sign issue as in the main text; the caption and text should be corrected together.","section":"Fig. 2 caption"},{"comment":"The third row of the matrix M in Eq. (D4) appears garbled: it reads '1 cot( γ/2) 1 cot( γ/2) 1 cot( γ/2)' and is missing apparent entries for the C0,0 and C0,1 columns. Please check the typesetting.","section":"Appendix D, matrix M"},{"comment":"The expression for N(A) contains 'ln ϵ + ln(1−cos²θ0) − ln|bm0|' inside a fraction; verify that the division by cos(θ0) applies as intended and that the bracketing in the floor expression is correct. Also, the product in Eq. (23) for n=0 is an empty product; define it as 1 for clarity.","section":"Appendix E, Eq. (E4)"}],"recommendation":"major_revision","confidential_remarks":"The reader's report and the stress-test note highlight a specific sign/calibration error in Eq. (2) that I have verified: for the stated α0 and α1, the proton concentration increases from RC to ATP-S, reversing the direction of the active driving force. This is a load-bearing error, but it is correctable in principle (e.g., by fixing the sign of α0 and recalibrating against Ref. 29). The ad hoc nature of the functionality score is a further concern that will require either a derivation or a robustness analysis. I recommend major revision rather than rejection because the derivation methodology is sound and the central mechanism remains plausible once the calibration is fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The stress-test note holds up, and it hits the load-bearing part of the paper. Eq. (2) with the reported α0 = −6×10⁻⁸ mol/L and α1 = 7.9×10⁻⁸ mol/L gives a proton concentration that increases from the RC pole (θ≈0.1) to the ATP-synthase equator (θ=π/2), exactly opposite to the decrease the text describes and the biology requires. At the RC, h is about −1×10⁻⁷, so the total proton concentration nearly vanishes and pH is undefined; at the equator h is positive, giving pH ≈ 6.75 rather than the cited 7.1. So the fitted profile is both inconsistent with the experimental values and reversed relative to the physical gradient. Since the active stress enters as βh, the driving force in Eq. (17) is flipped. The flat-pole/curved-equator morphology and the optimal region near (A=6, ξ=−1) are computed for a gradient flowing from the equator to the pole, not from RC to ATP-S. The central biological claim is therefore not established for the real system.\n\nThe derivation itself is the strong part. The stress-tensor formalism is applied consistently, the linearized shape equation is correct, and the appendices give explicit solutions for the degenerate A=0 and A=2 cases, plus a bound on the series truncation. Extending the group's earlier flat/cylindrical active-membrane work to a closed spherical vesicle with protein-mimicking boundary conditions is a legitimate step. The authors also honestly acknowledge the missing cardiolipin-sorting feedback.\n\nThe functionality score is the next soft spot. S is built from the exact features the model is designed to produce—large equator, short proton path, flat RC zone—so maximizing S in parameter space is partly circular. The weights and cutoffs are arbitrary, and the claimed comparison to experimental measurements is qualitative rather than a quantitative overlay. The coupling β is chosen by hand to give a 10% curvature change, with no sensitivity analysis.\n\nBecause the calibration error is correctable and the mathematical core is sound, this deserves referee time despite the flawed input. But I would not cite it in its current form, and the reading group will get more from it as a case study in how a sign error can invert a conclusion than as a reliable result. My recommendation: send to peer review, and require the authors to refit h to ref. 29, verify the sign of α0, and show what changes in Fig. 5 under the correct gradient.","headline":"Careful derivation, but the proton-concentration fit that drives the deformation is reversed, so the crista-like shape and the functional phase diagram are not established for a real RC-to-ATP-S gradient.","tokens_in":15993,"tokens_out":7962,"would_cite":false,"duration_ms":67832,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A purely diffusive proton gradient from pole to equator can deform a spherical, pH-sensitive vesicle into a crista-like shape with flat zones at the respiratory complexes and strongly curved zones at the ATP synthase, and the model's…","keywords":["mitochondrial cristae","pH gradient","spontaneous curvature","Helfrich Hamiltonian","active membrane","vesicle deformation","proton diffusion","cardiolipin"],"falsifier":"Measure the deformation profile of a cardiolipin-containing giant unilamellar vesicle under a controlled pole-to-equator pH gradient and compare it with Eq. (26) using independently measured values of the bending modulus, surface tension, and pH–curvature coupling: if the vesicle does not develop a flatter polar region and a more curved equator as the gradient is increased, or if the measured proton profile near the deformed membrane deviates strongly from the Laplace solution of Eq. (1), the central claim fails.","tokens_in":14935,"feed_emoji":"🫧","tokens_out":10108,"duration_ms":81328,"temperature":0.7,"pith_summary":"Mitochondrial cristae are folds of the inner membrane where ATP is made, and their shape—flat zones hosting respiratory complexes and tightly curved zones hosting ATP synthase—is part of how they work. This paper asks whether the proton current that drives ATP synthesis can itself sculpt that shape. The authors model a crista as a spherical vesicle whose membrane has a spontaneous curvature that depends linearly on the local proton concentration, with the protons diffusing passively from the poles to the equator. In the small-deformation regime they solve the resulting shape equation and find a parameter region in which the vesicle flattens at the poles and becomes highly curved at the equator, exactly the morphology of a functioning crista. If correct, the result shows that no protein scaffolding is needed: a pure diffusive pH gradient plus a pH-sensitive lipid is sufficient to organize the crista geometry that maximizes ATP output.","feed_headline":"A proton gradient alone can sculpt a vesicle into a working crista","feed_subtitle":"Pure diffusion of protons creates flat protein zones and curved ATP-synthase zones, matching crista anatomy.","key_machinery":"The machinery is a pH-dependent Helfrich Hamiltonian, $\\mathcal{H} = \\int_S d^2X \\,[\\tfrac{1}{2}\\kappa(C-C_0[h])^2 + \\sigma_0] - PV$, with the spontaneous curvature coupled linearly to the local proton concentration, $C_0(h)=C_0^0 + C_1^0 h(\\theta)$. The proton field $h(\\theta)$ is prescribed by the spherical Laplace equation $\\Delta_S h=0$ between the respiratory-complex boundary at $\\theta=\\theta_0$ and the ATP-synthase boundary at $\\theta=\\pi/2$, giving $h(\\theta)=-\\alpha_0 \\ln(\\tan(\\theta/2))+\\alpha_1$ fitted to measured crista pH values. The normal force balance reduces to a fourth-order linear ODE for the radial deformation $u(\\theta)$; its solution is the sum of a particular solution proportional to $h(\\theta)$ and kernel modes built from the functions $C_{A,b}(\\theta)$, with coefficients fixed by five biological boundary conditions (flat RC patch, ATP-S opening angle, and force balances at both protein sites) plus volume conservation. A functionality score $S=S_1-S_2-S_3$, combining equatorial perimeter, proton path length, and pole flatness, classifies the resulting shapes and locates the well-functioning crista region in the $(A,\\xi)$ plane.","core_discovery":"On the paper's own terms, the central discovery is that a closed spherical membrane with a pH-dependent spontaneous curvature, driven only by a diffusive proton concentration field, can adopt shapes with strongly inhomogeneous curvature. The authors state in the Discussion that 'we show that such a system can adopt shapes with inhomogeneous curvature. In the relevant parameter space, we identify the zone leading to both highly curved in the ATP-S location and flat zones in the RC location. These shapes correspond to a well-functioning crista.' Concretely, the shape equation $(\\Delta_S+2)(\\Delta_S+A)u = -[\\Delta_S+2(\\xi-1)]\\beta h + \\delta p$, solved with boundary conditions encoding the flat respiratory-complex patch and the opening angle of ATP synthase, yields a deformation with near-zero mean curvature near the pole and enhanced curvature at the equator; the model's functionality score peaks at $A=6$, $\\xi=-1$, corresponding to a typical crista radius $R\\approx 100$ nm and a spontaneous curvature $C_0\\approx -0.01$ nm$^{-1}$.","pith_inferences":["The paper leaves open how cardiolipin sorting would feed back on the pH field; coupling diffusion to curvature would likely amplify the predicted flat-pole/curved-equator split, since cardiolipin concentrates in curved zones and raises local pH sensitivity.","The optimal operating point $A=6$ sits at an eigenvalue of the operator that controls spherical stability, so the model implicitly predicts that the best-functioning crista is near a shape instability; varying $A$ across this value in the same calculation would test how robust the well-functioning zone is.","A direct experiment the paper does not propose: reverse the proton source and sink, placing RC at the equator and ATP-S at the poles; the symmetry of the equations then predicts a shape with the flat and curved zones swapped, which would distinguish this passive-diffusion mechanism from protein-scaffolding models."],"forward_implications":["A purely diffusive proton field is sufficient to create the flat-pole/curved-equator anatomy of a crista; the respiratory complexes and ATP synthase act only as boundary conditions, not as scaffolds that hold the shape.","The model predicts an optimal operating point ($A=6$, $\\xi=-1$) that corresponds, for $R\\approx 100$ nm, to a spontaneous curvature $C_0\\approx -0.01$ nm$^{-1}$, a value compatible with molecular-dynamics estimates for cardiolipin-containing membranes.","The phase diagram of the functionality score delimits, in the space of reduced pressure and spontaneous curvature, which parameter combinations yield well-functioning cristae; the measured lateral pH gradient (pH 6.4 near the complexes, pH 7.1 near ATP synthase) is used as input, so the diagram can be tested against pH-sensitive vesicle experiments.","Without a proton flux the vesicle stays spherical; the proton field is the sole driver of deformation, so the model makes the testable prediction that turning off the pH gradient erases the crista-like shape."],"supporting_citations":[{"why":"Supplies the measured lateral pH gradient (pH 6.4 near the respiratory complexes, pH 7.1 near ATP synthase) used to fit the proton concentration h(θ).","marker":"29"},{"why":"Defines the elastic membrane energy (Helfrich Hamiltonian) that the paper modifies with a pH-dependent spontaneous curvature.","marker":"18"},{"why":"In vitro experiment showing a local pH gradient remodels cardiolipin-containing giant vesicles into crista-like invaginations; primary motivation for the pH–curvature coupling.","marker":"14"},{"why":"Shows pH-induced lipid packing variations in cardiolipin bilayers, identifying the physical driving force behind the crista-like shape instability.","marker":"15"},{"why":"Molecular-dynamics evidence that proton accumulation near a cardiolipin-containing membrane modifies its spontaneous curvature, supporting the sign and linear form of C0(h).","marker":"16"},{"why":"Earlier model of a membrane deformed by a surface proton gradient on a cylindrical geometry; the stress-tensor framework used here extends it to a closed sphere.","marker":"26"},{"why":"Derives membrane deformations driven by a surface gradient and identifies bulges and necks reminiscent of cristae; direct methodological predecessor.","marker":"27"},{"why":"Provides the second-variation stability equation for spherical membranes, whose homogeneous part the paper identifies with its no-flux limit.","marker":"19"},{"why":"Establishes the macromolecular organization of ATP synthase and complex I in cristae, fixing the biological locations of flat and curved zones used as boundary conditions.","marker":"2"},{"why":"Molecular-dynamics estimate of spontaneous curvature of cardiolipin-containing mitochondrial membranes; used to compare the optimal C0≈−0.01 nm⁻¹.","marker":"33"}],"fun_headline_variants":["Proton gradient alone sculpts a vesicle into a working crista","pH gradient deforms vesicle into crista-like shape","Proton flow shapes membrane into crista pattern","Diffusive protons alone craft flat and curved crista zones","One proton gradient yields functional crista geometry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the proton concentration is set by pure diffusion on the undeformed spherical surface and does not change when the membrane deforms, with the pH–curvature coupling strength β chosen by hand so that the induced spontaneous curvature is about 10% of the sphere's curvature; if the shape feeds back on the proton distribution, or if β is much weaker, the predicted flat-pole/curved-equator shape need not appear.","fun_headline_variants_meta":{"raw":{"variants":["Proton gradient alone sculpts a vesicle into a working crista","pH gradient deforms vesicle into crista-like shape","Proton flow shapes membrane into crista pattern","Diffusive protons alone craft flat and curved crista zones","One proton gradient yields functional crista geometry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000497,"raw_usage":{"total_tokens":2412,"prompt_tokens":897,"completion_tokens":1515,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":513,"completion_tokens_details":{"reasoning_tokens":1446}},"tokens_in":513,"tokens_out":1515,"duration_ms":10660,"temperature":1.0,"reasoning_tokens":1446,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T14:57:36.379218+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the deformation profile of a cardiolipin-containing giant unilamellar vesicle under a controlled pole-to-equator pH gradient and compare it with Eq. (26) using independently measured values of the bending modulus, surface tension, and pH–curvature coupling: if the vesicle does not develop a flatter polar region and a more curved equator as the gradient is increased, or if the measured proton profile near the deformed membrane deviates strongly from the Laplace solution of Eq. (1), the central claim fails.","supporting_citations":[{"cited_title":"Rieger, W","cited_arxiv_id":null,"evidence_quote":"Supplies the measured lateral pH gradient (pH 6.4 near the respiratory complexes, pH 7.1 near ATP synthase) used to fit the proton concentration h(θ)."},{"cited_title":"Helfrich","cited_arxiv_id":null,"evidence_quote":"Defines the elastic membrane energy (Helfrich Hamiltonian) that the paper modifies with a pH-dependent spontaneous curvature."},{"cited_title":"Khalifat, N","cited_arxiv_id":null,"evidence_quote":"In vitro experiment showing a local pH gradient remodels cardiolipin-containing giant vesicles into crista-like invaginations; primary motivation for the pH–curvature coupling."},{"cited_title":"Khalifat, J-B","cited_arxiv_id":null,"evidence_quote":"Shows pH-induced lipid packing variations in cardiolipin bilayers, identifying the physical driving force behind the crista-like shape instability."},{"cited_title":"Allolio and D","cited_arxiv_id":null,"evidence_quote":"Molecular-dynamics evidence that proton accumulation near a cardiolipin-containing membrane modifies its spontaneous curvature, supporting the sign and linear form of C0(h)."},{"cited_title":"Patil, S","cited_arxiv_id":null,"evidence_quote":"Earlier model of a membrane deformed by a surface proton gradient on a cylindrical geometry; the stress-tensor framework used here extends it to a closed sphere."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives membrane deformations driven by a surface gradient and identifies bulges and necks reminiscent of cristae; direct methodological predecessor."},{"cited_title":"Helfrich","cited_arxiv_id":null,"evidence_quote":"Provides the second-variation stability equation for spherical membranes, whose homogeneous part the paper identifies with its no-flux limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the macromolecular organization of ATP synthase and complex I in cristae, fixing the biological locations of flat and curved zones used as boundary conditions."},{"cited_title":"Konar, H","cited_arxiv_id":null,"evidence_quote":"Molecular-dynamics estimate of spontaneous curvature of cardiolipin-containing mitochondrial membranes; used to compare the optimal C0≈−0.01 nm⁻¹."}],"review_version":1}