pith. sign in

arxiv: 2603.15170 · v2 · submitted 2026-03-16 · ❄️ cond-mat.soft · physics.bio-ph

Cracking donuts and sorting lipids: Geometry controls archaeal membrane stability and lipid organization

Felix Frey , Miguel Amaral , An{\dj}ela \v{S}ari\'c This is my paper

Pith reviewed 2026-05-15 10:21 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.bio-ph
keywords archaeal membranesbolalipidsmembrane curvaturelipid sortingmolecular dynamicsmembrane stabilitytoroidal vesiclescurvature coupling
0
0 comments X

The pith

Curvature sorts lipids and destabilizes bolalipid monolayers in archaeal membranes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper uses simulations of toroidal vesicles to examine how mean and Gaussian curvatures influence the stability and lipid arrangement in archaeal membranes containing both bilayer lipids and bolalipids. Bilayer membranes remain intact across all curvature combinations tested, while bolalipid monolayer membranes tend to reshape or break apart. Lipids that mimic bilayers, including u-shaped bolalipids, collect in zones of high mean curvature regardless of Gaussian curvature. The work positions this geometry-driven lipid segregation as a basic physical process underlying how archaeal membranes can remodel during cellular activities.

Core claim

Coarse-grained molecular dynamics simulations of toroidal vesicles show that soft bilayer membranes sustain all imposed curvatures, whereas rigid bolalipid monolayer membranes transition to other shapes or rupture. Bilayer-mimicking u-shaped bolalipids and bilayer lipids accumulate preferentially in regions of high mean curvature, independent of Gaussian curvature. The simulations thereby identify curvature-composition coupling as the physical signature of archaeal membrane remodeling.

What carries the argument

Curvature-composition coupling arising from the conformational flexibility of bolalipids (straight for monolayers, u-shaped for bilayers) and their differential response to mean curvature in mixed membranes.

If this is right

  • Bilayer lipids and u-shaped bolalipids will concentrate at sites of high mean curvature during membrane budding, fission, or vesicle formation.
  • Monolayer bolalipid membranes will require specific overall shapes or additional stabilizers to avoid rupture or shape transitions under bending.
  • The ratio of straight to u-shaped bolalipids will determine overall membrane stability as curvature increases.
  • Gaussian curvature will play little role in driving lipid sorting compared with mean curvature.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This mechanism could allow archaea to locally tune membrane thickness or fluidity purely through geometry during division without needing dedicated protein machinery.
  • Synthetic vesicles built from similar lipid mixtures might be designed to self-sort components into specific curvature zones for targeted release or shape control.
  • Live-cell imaging of archaea with fluorescently labeled bolalipids under varying osmotic or mechanical stresses could directly test whether the simulated sorting occurs in vivo.

Load-bearing premise

The coarse-grained molecular dynamics model faithfully captures the conformational flexibility, packing, and mechanical response of archaeal bolalipids and bilayer lipids under imposed mean and Gaussian curvatures.

What would settle it

Experimental measurements on real archaeal membranes or reconstituted systems under controlled toroidal curvatures showing no preferential accumulation of bilayer lipids at high mean curvature regions, or showing equal stability between monolayers and bilayers, would falsify the central claim.

Figures

Figures reproduced from arXiv: 2603.15170 by An{\dj}ela \v{S}ari\'c, Felix Frey, Miguel Amaral.

Figure 1
Figure 1. Figure 1: Computational model, toroidal vesicle geometry and curvature distribution on toroidal vesicles. (A, left) [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Bolalipid rigidity controls stability of toroidal vesicles composed of bolalipids. (A) Shape diagram of pure [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Bilayer lipid fraction controls stability of toroidal vesicles of archaeal mixture membranes composed of bilayer [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Curvature controls lipid distribution in pure bolalipid and archaeal mixture membranes. (A) For [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
read the original abstract

Cells are defined by lipid membranes that differ in their structure across the tree of life. While the membranes of most bacteria and eukaryotes consist of single-headed bilayer lipids, the membranes of archaea are composed of mixtures of single-headed bilayer lipids and double-headed bolalipids. Archaeal bolalipids can adopt straight or u-shaped conformations, enabling them - together with bilayer lipids - to control whether membranes form bilayer or monolayer structures. Yet, the physical principles governing archaeal membranes remain largely unexplored, especially how membrane structure couples to externally imposed curvature during membrane remodeling. Here, we perform coarse-grained molecular dynamics simulations of toroidal vesicles to systematically probe the effects of all relevant combinations of mean and Gaussian curvatures on shape stability and lipid organization. We find that soft bilayer membranes can sustain all curvatures induced, whereas rigid bolalipid monolayer membranes either transition to different vesicle shapes or rupture. Bilayer-mimicking u-shaped bolalipids and bilayer lipids are spatially accumulated in regions of high mean membrane curvature independent of Gaussian curvature. Our work identifies curvature-composition coupling as a physical signature of archaeal membrane remodeling.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports coarse-grained molecular dynamics simulations of toroidal vesicles formed from mixtures of archaeal bolalipids (capable of straight or u-shaped conformations) and bilayer lipids. It examines the effects of all combinations of mean and Gaussian curvature on vesicle stability and lateral lipid organization, finding that soft bilayer membranes remain stable under imposed curvatures while rigid bolalipid monolayers either reshape or rupture, and that u-shaped bolalipids together with bilayer lipids accumulate preferentially in regions of high mean curvature independent of Gaussian curvature. The central claim is that curvature-composition coupling constitutes a physical signature of archaeal membrane remodeling.

Significance. If the coarse-grained model is shown to be reliable, the work provides a systematic computational exploration of how geometry couples to composition in mixed archaeal membranes, which is relevant to understanding membrane stability in extremophiles and to biophysical mechanisms of remodeling. The toroidal geometry approach to independently vary mean and Gaussian curvature is a clear methodological strength.

major comments (2)
  1. [Methods] Methods section on coarse-grained model parameterization: no direct validation is reported for the free-energy barrier between straight and u-shaped bolalipid conformations, nor for their packing density and bending rigidity under simultaneous nonzero mean and Gaussian curvature. This mapping is load-bearing for the stability and sorting claims.
  2. [Results] Results on lipid sorting (e.g., the reported accumulation of u-shaped bolalipids at high mean curvature): the manuscript does not present error estimates, convergence checks across independent runs, or controls for finite-size effects on the toroidal vesicles, leaving the robustness of the curvature-independent sorting result unclear.
minor comments (2)
  1. [Figures] Figure captions for the toroidal vesicle snapshots could explicitly state the imposed curvature values and the simulation time window over which the configurations are averaged.
  2. [Abstract] The abstract states that 'all relevant combinations of mean and Gaussian curvatures' were probed but does not list the specific numerical ranges or the number of independent trajectories per condition.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive assessment of our work and for highlighting the need to strengthen the validation of the coarse-grained model and the statistical robustness of the lipid-sorting results. We address each major comment below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Methods] Methods section on coarse-grained model parameterization: no direct validation is reported for the free-energy barrier between straight and u-shaped bolalipid conformations, nor for their packing density and bending rigidity under simultaneous nonzero mean and Gaussian curvature. This mapping is load-bearing for the stability and sorting claims.

    Authors: We agree that explicit validation of the conformational free-energy barrier and the mechanical properties under combined curvatures would strengthen the manuscript. The coarse-grained model was parameterized from all-atom reference simulations and experimental packing data for archaeal lipids, but direct comparisons were not included in the original submission. In the revised version we will add a dedicated Methods subsection together with supplementary figures that report (i) the free-energy profile for the straight-to-u-shaped transition obtained from umbrella sampling and (ii) bending-rigidity estimates extracted from fluctuation spectra on both flat bilayers and toroidal geometries. These additions will directly support the stability and sorting claims. revision: yes

  2. Referee: [Results] Results on lipid sorting (e.g., the reported accumulation of u-shaped bolalipids at high mean curvature): the manuscript does not present error estimates, convergence checks across independent runs, or controls for finite-size effects on the toroidal vesicles, leaving the robustness of the curvature-independent sorting result unclear.

    Authors: We acknowledge that the original manuscript omitted quantitative error estimates and convergence diagnostics for the composition profiles. In the revision we will include (i) standard deviations computed from at least three independent trajectories for each curvature combination, (ii) time-series plots demonstrating that the sorting metrics have converged within the production runs, and (iii) a finite-size control comparing composition profiles on tori whose major radius differs by a factor of two. These additions will confirm that the preferential accumulation of u-shaped bolalipids and bilayer lipids at high mean curvature is independent of Gaussian curvature and robust to system size. revision: yes

Circularity Check

0 steps flagged

No circularity: direct simulation outputs without derivation or fitting

full rationale

The manuscript reports results exclusively from coarse-grained molecular dynamics simulations of toroidal vesicles under imposed curvatures. No mathematical derivation chain exists, no parameters are fitted to the target lipid sorting or stability data, and no self-citation is invoked to justify a uniqueness theorem or ansatz that would reduce the central claim to its own inputs. Curvature-composition coupling emerges as an observed outcome of the model dynamics rather than being defined or predicted by construction from prior results within the paper.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The study rests on the validity of a standard coarse-grained lipid model whose interaction parameters and conformational rules are taken from prior literature; no new entities are introduced.

free parameters (1)
  • Coarse-grained lipid interaction parameters
    Standard CG force-field parameters that control head-group repulsion, tail packing, and bolalipid bending stiffness; typically tuned to experimental bilayer properties.
axioms (1)
  • domain assumption The chosen coarse-grained representation faithfully reproduces the mechanical response and conformational ensemble of real archaeal bolalipids under curvature
    Invoked throughout the simulation protocol and necessary for all stability and sorting conclusions.

pith-pipeline@v0.9.0 · 5506 in / 1260 out tokens · 63415 ms · 2026-05-15T10:21:29.545841+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    (G) Normalized Gaussian curvature on the Clifford torus. (H) Toroidal vesicles for pure bolalipid membranes atk bola = 0k BT(top), atk bola = 1k BT(center) and a mixture membrane with bilayer content off bi = 0.7 and bolalipids atk bola = 2k BT. mixture membranes, both u-shaped bolalipids and bi- layer lipids accumulate in regions of high mean cur- vature...

    work page

  2. [2]

    4B) andf bi (Fig

    and the analysis ofu f (Fig. 4B) andf bi (Fig. 4D and Fig. S7). Simulationρ(geometry)ρ(lipid composition) Bolalipid torusρ= 1.39ρ= 1.48 Bilayer torusf bi = 0.7ρ= 1.36ρ= 1.30 Bilayer torusf bi = 0.8ρ= 1.38ρ= 1.20 Bilayer torusf bi = 0.9ρ= 1.36ρ= 1.31 enrichment of u-shaped bolalipids. To address this, we focus on mixture membranes with bilayer lipid frac- ...

    work page arXiv

  3. [3]

    Alberts, A

    B. Alberts, A. Johnson, J. Lewis, D. Morgan, M. Raff, K. Roberts, and P. Walter,Molecular biology of the cell: Sixth edition(Garland science, New York, 2015)

    work page 2015

  4. [4]

    S. B. Gould, Curr. Biol.28, R381 (2018)

    work page 2018

  5. [5]

    Albers and B

    S.-V. Albers and B. H. Meyer, Nat. Rev. Microbiol. 9, 414 (2011)

    work page 2011

  6. [6]

    D. P. Brownholland, G. S. Longo, A. V. Struts, M. J. Justice, I. Szleifer, H. I. Petrache, M. F. Brown, and D. H. Thompson, Biophysical Journal97, 2700 (2009)

    work page 2009

  7. [7]

    Bhattacharya, I

    A. Bhattacharya, I. D. Falk, F. R. Moss, T. M. Weiss, K. N. Tran, N. Z. Burns, and S. G. Boxer, Chem. Sci. 15, 14273 (2024)

    work page 2024

  8. [8]

    Tarrason Risa, F

    G. Tarrason Risa, F. Hurtig, S. Bray, A. E. Hafner, L. Harker-Kirschneck, P. Faull, C. Davis, D. Pap- atziamou, D. R. Mutavchiev, C. Fan, L. Meneguello, A. Arashiro Pulschen, G. Dey, S. Culley, M. Kilkenny, D. P. Souza, L. Pellegrini, R. A. M. de Bruin, R. Henriques, A. P. Snijders, A. ˇSari´ c, A.-C. Lind˚ as, N. P. Robinson, and B. Baum, Science369, eaa...

    work page 2020

  9. [9]

    J. Liu, N. Soler, A. Gorlas, V. Cvirkaite-Krupovic, M. Krupovic, and P. Forterre, microLife2, uqab007 (2021)

    work page 2021

  10. [10]

    Warner and B

    J. Warner and B. O’Shaughnessy, Biophys. J.103, 689 (2012)

    work page 2012

  11. [11]

    R. K. W. Spencer, I. Santos-P´ erez, I. Rodr´ ıguez- Renovales, J. M. Martinez Galvez, A. V. Shnyrova, and M. M¨ uller, Nat. Commun.15, 2793 (2024)

    work page 2024

  12. [12]

    Mu˜ noz-Basagoiti, F

    M. Mu˜ noz-Basagoiti, F. Frey, B. Meadowcroft, M. Amaral, A. Prada, and A. ˇSari´ c, Soft Matter21, 7736 (2025)

    work page 2025

  13. [13]

    Seifert, Adv

    U. Seifert, Adv. Phys.46, 13 (1997)

    work page 1997

  14. [14]

    Deserno, Chem

    M. Deserno, Chem. Phys. Lipids.185, 11 (2015)

    work page 2015

  15. [15]

    H. T. McMahon and J. L. Gallop, Nature438, 590 (2005)

    work page 2005

  16. [16]

    Zimmerberg and M

    J. Zimmerberg and M. M. Kozlov, Nat. Rev. Mol. Cell Biol.7, 9 (2006)

    work page 2006

  17. [17]

    Frey and T

    F. Frey and T. Idema, Soft Matter17, 3533 (2021)

    work page 2021

  18. [18]

    P. L.-G. Chong, Chem. Phys. Lipids163, 253 (2010)

    work page 2010

  19. [19]

    Tourte, P

    M. Tourte, P. Schaeffer, V. Grossi, and P. M. Oger, Front. Microbiol.11, 526 (2020)

    work page 2020

  20. [20]

    I. R. Cooke and M. Deserno, J. Chem. Phys.123, 224710 (2005)

    work page 2005

  21. [21]

    Amaral, F

    M. Amaral, F. Frey, X. Jiang, B. Baum, and A. ˇSari´ c, eLife14, RP105432 (2025)

    work page 2025

  22. [22]

    Ou-Yang Zhong-can, Phys. Rev. A41, 4517 (1990)

    work page 1990

  23. [23]

    M. Hu, J. J. Briguglio, and M. Deserno, Biophys. J. 102, 1403 (2012)

    work page 2012

  24. [24]

    Wohlert, W

    J. Wohlert, W. K. den Otter, O. Edholm, and W. J. Briels, J. Chem. Phys.124, 154905 (2006)

    work page 2006

  25. [25]

    I. R. Cooke and M. Deserno, Biophysical Journal91, 487 (2006)

    work page 2006

  26. [26]

    Derganc, Phys

    J. Derganc, Phys. Biol.4, 317 (2007)

    work page 2007

  27. [27]

    Callan-Jones, B

    A. Callan-Jones, B. Sorre, and P. Bassereau, Cold Spring Harb. Perspect. Biol.3, a004648 (2011)

    work page 2011

  28. [28]

    K¨ onig, R

    M. K¨ onig, R. De Vries, F. Gr¨ unewald, S. Marrink, and W. Pezeshkian, bioRxiv 10.1101/2023.12.15.571845 (2023)

    work page doi:10.1101/2023.12.15.571845 2023

  29. [29]

    Sorre, A

    B. Sorre, A. Callan-Jones, J.-B. Manneville, P. Nas- soy, J.-F. Joanny, J. Prost, B. Goud, and P. Bassereau, Proc. Natl. Acad. Sci. U.S.A.106, 5622 (2009)

    work page 2009

  30. [30]

    Baoukina, H

    S. Baoukina, H. I. Ing´ olfsson, S. J. Marrink, and D. P. Tieleman, Advcd Theory and Sims1, 1800034 (2018)

    work page 2018

  31. [31]

    Mutz and D

    M. Mutz and D. Bensimon, Phys. Rev. A43, 4525 (1991)

    work page 1991

  32. [32]

    Tian and T

    A. Tian and T. Baumgart, Biophys. J.96, 2676 (2009)

    work page 2009

  33. [33]

    M. M. Kamal, D. Mills, M. Grzybek, and J. Howard, Proc. Natl. Acad. Sci. U.S.A.106, 22245 (2009)

    work page 2009

  34. [34]

    Beltr´ an-Heredia, F.-C

    E. Beltr´ an-Heredia, F.-C. Tsai, S. Salinas-Almaguer, F. J. Cao, P. Bassereau, and F. Monroy, Commun Biol 2, 225 (2019)

    work page 2019

  35. [35]

    P¨ ohnl, M

    M. P¨ ohnl, M. F. W. Trollmann, and R. A. B¨ ockmann, Nat. Commun.14, 8038 (2023)

    work page 2023

  36. [36]

    K. S. Makarova, M. Y. Galperin, and E. V. Koonin, Biochimie118, 302 (2015)

    work page 2015

  37. [37]

    Verch` ere, W.-L

    A. Verch` ere, W.-L. Ou, B. Ploier, T. Morizumi, M. A. Goren, P. B¨ utikofer, O. P. Ernst, G. Khelashvili, and A. K. Menon, Sci. Rep.7, 9522 (2017)

    work page 2017

  38. [38]

    A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, Comput. Phys. Commun.271, 108171 (2022)

    work page 2022

  39. [39]

    Stukowski, Model Simul

    A. Stukowski, Model Simul. Mat. Sci. Eng.18, 015012 (2009). 1 Supplementary material Cracking donuts and sorting lipids: geometry controls archaeal membrane stability and lipid organization Felix Frey, Miguel Amaral, Andela ˇSari´ c I. COMPUTATIONAL MODEL A. Coarse-grained particle based model of bilayer and bolalipids To model bilayer and bolalipids we f...

    work page 2009

  40. [40]

    In the model, a bilayer lipid consists of one head bead and two tail beads (Fig

    and our extension to bolalipids [19]. In the model, a bilayer lipid consists of one head bead and two tail beads (Fig. 1B). A bolalipid is represented by joining two bilayer lipids along their tails, so it consists of two head beads that are connected by four tail beads (Fig. 1B). As a result, bilayer and bolalipids contain identical head and tail groups....

    work page 2000

  41. [41]

    Depending on whether we simulate a bolalipid or mixture membrane, we either place one bolalipid or two bilayer lipids lined up in radial direction at each sampled position. In Table I we summarize the number of bolalipids (# bolalipids), bilayer lipids (# bilayer lipids), the total number of lipids (# lipids) and the lipid bead number (# beads) used for d...

    work page