pith. sign in

arxiv: 2502.09798 · v1 · submitted 2025-02-13 · ⚛️ physics.bio-ph · cond-mat.soft

A tutorial for mesoscale computer simulations of lipid membranes: tether pulling, tubulation and fluctuation

Pith reviewed 2026-05-23 03:13 UTC · model grok-4.3

classification ⚛️ physics.bio-ph cond-mat.soft
keywords coarse-grained membrane simulationslipid membranestether pullingtubulationmembrane fluctuationsmesoscale modelingbiophysics simulations
0
0 comments X

The pith

Coarse-grained simulations of lipid membranes become accessible through tutorials on tether pulling and tubulation.

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

This paper introduces non-experts to coarse-grained membrane simulations at the mesoscale by first overviewing modeling approaches and offering guidance on developing such simulations. It then delivers a hands-on tutorial that applies three different membrane models to the specific cases of tether pulling, tubulation, and fluctuations. The work compares the models on scope and computational demands while supplying a public repository of ready-to-run tutorials to lower entry barriers for researchers in soft matter and biophysics.

Core claim

The paper establishes that a pedagogical tutorial using three representative coarse-grained membrane models can teach the practical application of mesoscale simulations to membrane deformations, with explicit discussion of each model's applicability range and resource requirements together with executable code examples.

What carries the argument

Three distinct coarse-grained membrane models applied to tether pulling, tubulation, and fluctuation scenarios, functioning as the central pedagogical examples with direct comparisons of their scope and cost.

If this is right

  • Researchers can immediately reproduce tether pulling and tubulation results using the provided models and code.
  • Model choice can be guided by matching the required deformation type to the model's documented scope and computational cost.
  • Fluctuation analysis becomes a practical entry point for studying membrane properties across different modeling resolutions.
  • The tutorial structure supports extension to additional membrane phenomena while maintaining the same comparative framework.

Where Pith is reading between the lines

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

  • The repository format could be adopted for tutorials on other soft-matter systems beyond membranes.
  • Direct comparison of the three models might highlight systematic differences in predicted deformation energies that warrant further analytical study.
  • Integration of the tutorial outputs with experimental membrane data could provide a low-cost validation route not explored in the paper.

Load-bearing premise

The three selected membrane models are representative of the main classes used in the field and the repository instructions are sufficient for a non-expert to execute and interpret the simulations without additional debugging.

What would settle it

A reader following the repository instructions who cannot successfully run and interpret the tether pulling or tubulation simulations without further external assistance would show that the tutorial does not achieve its intended accessibility.

Figures

Figures reproduced from arXiv: 2502.09798 by Adam Prada, An{\dj}ela \v{S}ari\'c, Billie Meadowcroft, Felix Frey, Maitane Mu\~noz-Basagoiti, Miguel Amaral.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Tether pulling simulation using a mesh model. (A) In a dynamically triangulated mesh, membrane fluidity is ensured [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (A) (top) Diagram of Cooke bilayer lipid and interaction matrix showing pair potential sketches. (bottom) Simulation [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (A) CPU time (hours) spent in simulating different membrane deformations using the membrane models in Section III. [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
read the original abstract

Lipid membranes and membrane deformations are a long-standing area of research in soft matter and biophysics. Computer simulations have complemented analytical and experimental approaches as one of the pillars in the field. However, setting up and using membrane simulations can come with barriers due to the multidisciplinary effort involved and the vast choice of existing simulations models. In this review, we introduce the non-expert reader to coarse-grained membrane simulations (CGMS) at the mesoscale. Firstly, we give a concise overview of the modelling approaches to study fluid membranes, together with guidance to more specialized references. Secondly, we provide a conceptual guide on how to develop CGMS. Lastly, we construct a hands-on tutorial on how to apply CGMS, by providing a pedagogical examination of tether pulling, tubulation and fluctuations with three different membrane models, and discussing them in terms of their scope and how resource-intensive they are. To ease the reader's venture into the field, we provide a repository with ready-to-run tutorials.

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

0 major / 1 minor

Summary. The manuscript is a tutorial/review introducing non-expert readers to coarse-grained mesoscale simulations (CGMS) of fluid lipid membranes. It first surveys modeling approaches for membranes, then offers a conceptual guide to developing such simulations, and finally provides hands-on pedagogical examples of tether pulling, tubulation, and fluctuation analysis performed with three distinct membrane models. The authors supply a public repository containing ready-to-run tutorials and discuss the scope and computational cost of each model.

Significance. If the provided models are representative and the repository instructions are executable without hidden barriers, the work would lower the entry threshold for researchers new to membrane simulations by combining conceptual orientation with immediately usable code. The explicit provision of a ready-to-run repository constitutes a concrete strength for reproducibility and pedagogy in a field where setup overhead is high.

minor comments (1)
  1. The abstract and introduction refer to 'three different membrane models' without naming them or their key distinguishing features in the opening sections; a brief table or explicit list early in the text would improve navigation for readers.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript and their recommendation to accept. We appreciate the recognition of the tutorial's value in providing both conceptual guidance and immediately usable code for non-experts in mesoscale membrane simulations.

Circularity Check

0 steps flagged

No significant circularity: tutorial with no derivations or predictions

full rationale

The manuscript is explicitly a tutorial and review that reuses three established membrane models for pedagogical demonstration of tether pulling, tubulation, and fluctuations. It contains no novel quantitative predictions, derivations, or fitted parameters that could reduce to inputs defined within the paper. The central contribution is instructional guidance plus a code repository; any self-citations are limited to background references for the pre-existing models and do not bear load on any claimed result. The paper is therefore self-contained against external benchmarks with no circular steps.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The work is a tutorial on established methods and therefore introduces no new free parameters, axioms, or invented entities; it relies on standard coarse-grained membrane models already present in the prior literature.

pith-pipeline@v0.9.0 · 5730 in / 1137 out tokens · 25026 ms · 2026-05-23T03:13:49.230574+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

180 extracted references · 180 canonical work pages

  1. [1]

    A detailed review on such models can be found in [123]

    Several-beads-per-lipid models Here we introduce models that represent each lipid as several (usually 3-7) particles. A detailed review on such models can be found in [123]. Several-beads-per-lipid models are generic, top-down models that do not address specific detail of the lipids. In these models, a lipid is rep- resented by one hydrophilic head partic...

  2. [2]

    Several similar one-bead mod- els deploying anisotropic potentials have been developed [114, 139–141]

    One-bead-per-lipid-patch models To simulate membranes at the largest possible scale with particle-based models, a high level of coarse- graining has been achieved by representing a patch of lipids by just one bead. Several similar one-bead mod- els deploying anisotropic potentials have been developed [114, 139–141]. One-bead models are unique in that they...

  3. [3]

    Us- ing typical numbers for bending rigidity κ and membrane tension γ, we expect tube radii below optical resolution, ranging from tens to hundreds of nanometres

    Simulation set-up Membrane tubes extruded by pulling a tether have been showed to follow predictions of the Helfrich ana- lytic theory, namely, a radii of Req = p κ/(2γ) and a required force for tether extrusion f = 2π√2κγ [43]. Us- ing typical numbers for bending rigidity κ and membrane tension γ, we expect tube radii below optical resolution, ranging fr...

  4. [4]

    Integrating the dynamics of the membrane TriLMP simulates the dynamics of a fluid membrane using an HMC approach, that is, by combining MC and MD stages throughout the integration. Specifically, a typical simulation run consists of a sequence of two stages: bond-swaps, which are determined by an MC al- gorithm (Metropolis criterion), and vertex moves, whi...

  5. [5]

    Protocol to measure the extrusion force Mimicking the optical tweezer experiments [9, 54], the desired membrane deformation can be achieved by teth- ering a bead to the membrane patch and pulling on it (Fig. 2A). In our tutorial, we use a bead that is tethered to the central vertex of the membrane through a har- monic bond. It is important to choose the e...

  6. [6]

    Results: analyzing the force-elongation profile The force-elongation profile during tether extrusion of a fluid membrane can be divided into several regimes that have been experimentally measured [54, 109]. At small elongations, one finds an elastic regime where the force required to deform the membrane is proportional to the elongation; here the membrane...

  7. [7]

    This set-up models a biological membrane as a single layer of particles interacting via a 2-body potential dependent on the relative distance and orientation of the particles (Fig

    The YLZ potential One of the most widely used 1-bead-per-lipid-patch models is the YLZ model; developed in 2010 and named after its authors[114]. This set-up models a biological membrane as a single layer of particles interacting via a 2-body potential dependent on the relative distance and orientation of the particles (Fig. 3A). When this potential is im...

  8. [8]

    The orientation vectors are all taken to point outside of the tube, normal to the surface of the cylinder

    Simulation set-up and protocol to measure the equilibrium radius To set up the initial state, we put YLZ particles on a regular trigonal lattice on the surface of a cylinder with the edge length roughly corresponding to the minimum of the radial potential ( ≈ 1σ, where σ sets the lengthscale of the system). The orientation vectors are all taken to point o...

  9. [9]

    Regulating membrane tension in a simulation The main conceptually difficult part of the simulation is the use of a barostat to control the tension of the mem- brane. As can be derived from the definition of a surface tension, the membrane tension of a tube can be related to its equilibrium radius and the force in the direction of 13 the tube axis exerted ...

  10. [10]

    3D is ln Req = x ln κ 2 − x ln γ (8) The simulations are consistent in the value of x (see Fig 3D), but slightly deviate from the theoretical value of 1/2 (see Fig

    Comparison with theory and computation of the bending rigidity κ Informed by the analytical formula for Req(γ), we plot ln Req as a function of ln γ and observe that these two quantities adhere extremely well to a power law, which can be represented as Req = κ 2γ x (7) The equation fitted to the data in Fig. 3D is ln Req = x ln κ 2 − x ln γ (8) The simula...

  11. [11]

    The results agree with the theory ex- tremely well

    Results: The effects of the tension γ and the bending rigidity κ The simulations were stable except under very high tensions and occasionally shedding a small fraction of the particles (which can be solved by adjusting barostat- ing parameters). The results agree with the theory ex- tremely well. For the same set of membrane parameters and the same membra...

  12. [12]

    k/idcjjicotU6JrXJsbY3MtVDbY=

    The three-beads-per-lipid model One of the most recent and widely used several-beads- per-lipid models is the Cooke membrane model [113], where the solvent is modelled implicitly. We refer the reader to its seminal paper[113] for specifics, and give only a brief overview here. 14 0.0 0.5 1.0 Simulation time / 𝜏×103 0 5 10 15 20Tube radius / 𝜎 1 0 2 4 6 Si...

  13. [13]

    We also used the same timestep of dt = 0.01τ

    Simulation set-up The seminal paper for the Cooke model [113] tuned two parameters, the tail interaction range w which we set to 1.5σ, and the temperature which we set up so thatkBT = 1ϵ; this guarantees that our system will be in the liquid phase. We also used the same timestep of dt = 0.01τ. To simulate this system, we pre-assemble a flat membrane in a ...

  14. [14]

    Analysis pipeline We generally follow the analysis provided in [113], in- formed by a latter work which focuses on fluctuation spectrum analysis [153]. Fig. 4C provides an overview of our analysis. First, we determine if the equilibration stage is complete by plotting together the 4 replicas. We verify qualitatively that the simulation has equilibrated wi...

  15. [15]

    4C, we show the membrane height fluctuation spectrum

    Results: analysing membrane fluctuations In Fig. 4C, we show the membrane height fluctuation spectrum. For the membrane parameters and the sim- ulation duration used, modes are well-sampled only if their wavenumber n ≥ 2, or equivalently, if their wave- length is smaller than half the length of the simulation box. We also observe that the initially linear...

  16. [16]

    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

  17. [17]

    Physics of Biological Membranes , ed. P. Bassereau and P. Sens, Springer International Publishing, Cham, 2018

  18. [18]

    Frey and T

    F. Frey and T. Idema, Soft Matter, 2021, 17, 3533–3549

  19. [19]

    Phillips, J

    R. Phillips, J. Kondev, J. Theriot and H. Garcia, Physi- cal biology of the cell, Garland Science, New York, 2012

  20. [20]

    Siggel, S

    M. Siggel, S. Kehl, K. Reuter, J. K¨ ofinger and G. Hum- mer, The Journal of Chemical Physics , 2022, 157, 174801

  21. [21]

    Dimova, Annu

    R. Dimova, Annu. Rev. Biophys., 2019, 48, 93–119

  22. [22]

    Idema and D

    T. Idema and D. J. Kraft, Curr. Opin. Colloid Interface Sci., 2019, 40, 58–69

  23. [23]

    Turlier and T

    H. Turlier and T. Betz, Annu. Rev. Condens. Matter Phys., 2019, 10, 213–232

  24. [24]

    Bassereau, B

    P. Bassereau, B. Sorre and A. L´ evy,Advances in Colloid and Interface Science, 2014, 208, 47–57

  25. [25]

    H. T. McMahon and J. L. Gallop, Nature, 2005, 438, 590–596

  26. [26]

    Lipowsky, Advanced Biology, 2022, 6, 2101020

    R. Lipowsky, Advanced Biology, 2022, 6, 2101020

  27. [27]

    Deserno, Chem

    M. Deserno, Chem. Phys. Lipids , 2015, 185, 11–45

  28. [28]

    Seifert, Adv

    U. Seifert, Adv. Phys., 1997, 46, 13–137

  29. [29]

    Lipowsky, Nature, 1991, 349, 475–481

    R. Lipowsky, Nature, 1991, 349, 475–481

  30. [30]

    Beiter and G

    J. Beiter and G. A. Voth, Curr. Opin. Cell Biol. , 2024, 87, 102831

  31. [31]

    L. Li, X. Huang and J. Wang, Giant, 2024, 18, 100263

  32. [32]

    A. K. Dasanna and D. A. Fedosov, Eur. Phys. J. Spec. Top., 2024

  33. [33]

    Kumar, S

    G. Kumar, S. C. Duggisetty and A. Srivastava, J Mem- brane Biol, 2022, 255, 757–777

  34. [34]

    K¨ uhlbrandt,Science, 2014, 343, 1443–1444

    W. K¨ uhlbrandt,Science, 2014, 343, 1443–1444

  35. [35]

    S. J. Sahl, S. W. Hell and S. Jakobs, Nat. Rev. Mol. Cell Biol., 2017, 18, 685–701

  36. [36]

    Pfitzner, V

    A.-K. Pfitzner, V. Mercier, X. Jiang, J. Moser von Filseck, B. Baum, A. ˇSari´ c and A. Roux, Cell, 2020, 182, 1140–1155

  37. [37]

    M. Mund, A. Tschanz, Y.-L. Wu, F. Frey, J. L. Mehl, M. Kaksonen, O. Avinoam, U. S. Schwarz and J. Ries, Journal of Cell Biology , 2023, 222, e202206038

  38. [38]

    Steink¨ uhler, R

    J. Steink¨ uhler, R. L. Knorr, Z. Zhao, T. Bhatia, S. M. Bartelt, S. Wegner, R. Dimova and R. Lipowsky, Nat. Commun., 2020, 11, 905

  39. [39]

    Azadbakht, B

    A. Azadbakht, B. Meadowcroft, T. Varkevisser, A.ˇSari´ c and D. J. Kraft, Nano Letters, 2023, 23, 4267–4273

  40. [40]

    Azadbakht, B

    A. Azadbakht, B. Meadowcroft, J. M´ ajek, A.ˇSari´ c and D. J. Kraft, Biophysical Journal, 2024, 123, 307–316

  41. [41]

    Willems, A

    V. Willems, A. Baron, D. A. M. Fernandez, G. Wolfis- berg, E. R. Dufresne, J.-C. Baret and L. Alvarez, Phase separation dependent active motion of Janus lipid vesicles, 2023, http://arxiv.org/abs/2311.00685, arXiv:2311.00685 [cond-mat, physics:physics]

  42. [42]

    S. J. Marrink, A. H. De Vries and D. P. Tieleman, Biochim. Biophys. Acta Biomembr. , 2009, 1788, 149– 168

  43. [43]

    ´Alvarez, J

    D. ´Alvarez, J. Sapia and S. Vanni, Current Opinion in Cell Biology, 2023, 83, 102205

  44. [44]

    A. H. Larsen, IJMS, 2022, 23, 8098

  45. [45]

    Noguchi, J

    H. Noguchi, J. Phys. Soc. Jpn. , 2009, 78, 041007

  46. [46]

    Guckenberger and S

    A. Guckenberger and S. Gekle, J. Phys.: Condens. Mat- ter, 2017, 29, 203001

  47. [47]

    Pezeshkian and S

    W. Pezeshkian and S. J. Marrink, Current Opinion in Cell Biology, 2021, 71, 103–111

  48. [48]

    Munoz-Basagoiti, F

    M. Munoz-Basagoiti, F. Frey, B. Meadowcroft, M. Ama- ral, A. Prada and A. Saric, Tutorials for mesoscale com- puter simulations of lipid membranes, https://github. com/Saric-Group/MembraneReviewTutorials

  49. [49]

    Helfrich and R.-M

    W. Helfrich and R.-M. Servuss, Il Nuovo Cimento D , 1984, 3, 137–151

  50. [50]

    Deserno, Phys

    M. Deserno, Phys. Rev. E , 2004, 69, 031903

  51. [51]

    Foret, The European Physical Journal E , 2014, 37, 42

    L. Foret, The European Physical Journal E , 2014, 37, 42

  52. [52]

    F. Frey, F. Ziebert and U. S. Schwarz, Phys. Rev. Lett., 2019, 122, 088102. 20

  53. [53]

    F. Frey, F. Ziebert and U. S. Schwarz, Phys. Rev. E , 2019, 100, 052403

  54. [54]

    Seifert, K

    U. Seifert, K. Berndl and R. Lipowsky, Phys. Rev. A , 1991, 44, 1182–1202

  55. [55]

    Christ, T

    S. Christ, T. Litschel, P. Schwille and R. Lipowsky, Soft Matter, 2021, 17, 319–330

  56. [56]

    Frey and T

    F. Frey and T. Idema, Phys. Rev. E, 2022, 106, 024401

  57. [57]

    A.-S. c. v. Smith, E. Sackmann and U. Seifert, Phys. Rev. Lett., 2004, 92, 208101

  58. [58]

    Der´ enyi, F

    I. Der´ enyi, F. J¨ ulicher and J. Prost,Phys. Rev. Lett. , 2002, 88, 238101

  59. [59]

    Seifert and R

    U. Seifert and R. Lipowsky, Phys. Rev. A , 1990, 42, 4768–4771

  60. [60]

    J¨ ulicher and R

    F. J¨ ulicher and R. Lipowsky, Phys. Rev. E , 1996, 53, 2670–2683

  61. [61]

    J¨ ulicher and R

    F. J¨ ulicher and R. Lipowsky,Phys. Rev. Lett., 1993, 70, 2964–2967

  62. [62]

    Dasgupta, T

    S. Dasgupta, T. Auth and G. Gompper, Soft Matter , 2013, 9, 5473–5482

  63. [63]

    Dasgupta, T

    S. Dasgupta, T. Auth and G. Gompper, Nano Lett. , 2014, 14, 687–693

  64. [64]

    Ni and G

    H. Ni and G. A. Papoian, J. Phys. Chem. B , 2021, 125, 10710–10719

  65. [65]

    H. R. Vutukuri, M. Hoore, C. Abaurrea-Velasco, L. van Buren, A. Dutto, T. Auth, D. A. Fedosov, G. Gompper and J. Vermant, Nature, 2020, 586, 52–56

  66. [66]

    P. Iyer, G. Gompper and D. A. Fedosov, Soft Matter , 2023, 19, 3436–3449

  67. [67]

    Pezeshkian and J

    W. Pezeshkian and J. H. Ipsen, Nat Commun, 2024, 15, 548

  68. [68]

    Weichsel and P

    J. Weichsel and P. L. Geissler, PLOS Computational Biology, 2016, 12, e1004982

  69. [69]

    Koster, A

    G. Koster, A. Cacciuto, I. Der´ enyi, D. Frenkel and M. Dogterom, Physical Review Letters , 2005, 94, 068101

  70. [70]

    K. A. Brakke, Experimental Mathematics, 1992, 1, 141– 165

  71. [71]

    C. Zhu, C. T. Lee and P. Rangamani, Biophys. Rep. , 2022, 2, 100062

  72. [72]

    D. A. Matoz-Fernandez, S. Li, M. O. de la Cruz and R. Sknepnek, PyMembrane: A flexible frame- work for efficient simulations of elastic and liquid membranes, 2023, http://arxiv.org/abs/2308.12754, arXiv:2308.12754 [cond-mat, physics:physics]

  73. [73]

    Dadunashvili and T

    G. Dadunashvili and T. Idema, flippy: User friendly and open source framework for lipid membrane sim- ulations, 2023, http://arxiv.org/abs/2303.12305, arXiv:2303.12305 [cond-mat, physics:physics]

  74. [74]

    T. S. Carpenter, C. A. L´ opez, C. Neale, C. Montour, H. I. Ing´ olfsson, F. Di Natale, F. C. Lightstone and S. Gnanakaran, Journal of Chemical Theory and Com- putation, 2018, 14, 6050–6062

  75. [75]

    R. S. Davis, P. B. Sunil Kumar, M. M. Sperotto and M. Laradji, The Journal of Physical Chemistry B , 2013, 117, 4072–4080

  76. [76]

    T. A. Wassenaar, H. I. Ing´ olfsson, R. A. B¨ ockmann, D. P. Tieleman and S. J. Marrink, Journal of Chemical Theory and Computation , 2015, 11, 2144–2155

  77. [77]

    S. Jo, J. B. Lim, J. B. Klauda and W. Im, Biophysical Journal, 2009, 97, 50–58

  78. [78]

    Y. Qi, H. I. Ing´ olfsson, X. Cheng, J. Lee, S. J. Marrink and W. Im, Journal of Chemical Theory and Computa- tion, 2015, 11, 4486–4494

  79. [79]

    S. O. Yesylevskyy, L. V. Sch¨ afer, D. Sengupta and S. J. Marrink, PLOS Computational Biology , 2010, 6, e1000810

  80. [80]

    Khalili-Araghi, J

    F. Khalili-Araghi, J. Gumbart, P.-C. Wen, M. So- tomayor, E. Tajkhorshid and K. Schulten, Current Opinion in Structural Biology , 2009, 19, 128–137

Showing first 80 references.