pith. sign in

arxiv: 2602.16171 · v2 · pith:QIQ6PRSQnew · submitted 2026-02-18 · ⚛️ physics.bio-ph

Self-Organized Bioelectricity via Collective Pump Alignment: Toward a Physical Origin of Chemiosmosis

Ryosuke Nishide , Kunihiko Kaneko This is my paper

Pith reviewed 2026-05-21 13:28 UTC · model grok-4.3

classification ⚛️ physics.bio-ph
keywords self-organized bioelectricitycollective pump alignmentchemiosmosisnonequilibrium phase transitionmembrane potentialion pumpsprotocells
0
0 comments X

The pith

Ion pumps spontaneously align through feedback with the electric fields they create, producing sustained membrane potentials.

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

The paper presents a minimal model in which ion pumps on a membrane interact through the potentials generated by their own directional transport. Transport creates an electrochemical gradient that orients neighboring pumps, establishing positive feedback that drives collective alignment. Simulations and mean-field calculations identify a nonequilibrium transition from random orientations with no net flow to an aligned state that maintains a membrane potential. The transition follows mean-field Ising critical behavior, yet the aligning field is produced internally by the transport process itself. If correct, the model supplies a physical route by which chemiosmotic coupling could arise in protocells without requiring pre-existing complex machinery.

Core claim

In the model, directional ion transport generates a membrane potential that biases pump orientation, leading to self-organized collective alignment. Numerical simulations and mean-field analysis reveal a nonequilibrium transition from a disordered state without net transport to a pump-alignment state with sustained membrane potentials. The critical behavior is consistent with the mean-field Ising universality class; however, the effective field is generated self-consistently by nonequilibrium ion transport. Protocell asymmetry can bias the polarity of the membrane potential.

What carries the argument

The self-consistent effective field produced by nonequilibrium ion transport that biases pump orientations and thereby drives collective alignment.

If this is right

  • Net directional ion transport emerges spontaneously from initially disordered pump configurations.
  • Sustained membrane potentials arise as a direct consequence of the aligned state.
  • The transition exhibits critical behavior belonging to the mean-field Ising class, driven by an internally generated field.
  • Asymmetry in a protocell can select the polarity of the resulting membrane potential.

Where Pith is reading between the lines

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

  • The same feedback could operate in early evolutionary settings before specialized chemiosmotic proteins evolved.
  • Synthetic membrane systems might be engineered to exploit similar self-alignment for controlled ion flow.
  • Extensions that include multiple pump types or varying membrane curvature could reveal additional stable regimes.

Load-bearing premise

The electrochemical potential generated by ion transport directly biases the orientation of individual pumps to produce positive feedback that favors collective alignment.

What would settle it

A simulation or experiment in which the bias linking membrane potential to pump orientation is removed or reversed, resulting in the absence of any sustained net transport or membrane potential, would falsify the mechanism.

Figures

Figures reproduced from arXiv: 2602.16171 by Kunihiko Kaneko, Ryosuke Nishide.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: , histograms (for 2000 samples) and time series (for one sample) of the pump alignment order parameter m = 1 NP ∑ i Pi = 2p − 1, (7) are plotted, where p = nP /NP . For small α (say 0.01), no pump alignment is observed. The distribu￾tion of m exhibits a single peak with zero mean and fluctuations around it ( [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a), as α or J increases, the random state (blue region) transitions to the aligned state (red region). Fig￾ure 3(b) shows that the flipping rate is enhanced near the transition region. Additionally, the order parameter ⟨|m|⟩ increases sharply and its variance shows a sharp peak near α = 0.1 at J = 5.5 (Figs. 3(c) and (d)). Hence, the change between the random and aligned states can be interpreted as a pha… view at source ↗
Figure 4
Figure 4. Figure 4: shows pump alignment remains biased inward even slightly beyond the transition line. This occurs because, for α close to αC , the fixed point q¯ = q− is close to zero, so fluctuations are strong and inward orientation is pref￾erentially observed. For sufficiently large α, alignment in both directions is realized. Furthermore, decreasing rV enhances the energy asymmetry and suppresses fluc￾tuations around q… view at source ↗
read the original abstract

Directional ion transport across membranes maintains living systems in nonequilibrium, which underlies chemiosmotic energy conversion. However, the physical origin of collectively organized ion transport in primitive cellular systems remains unclear. Here, we propose a minimal model in which ion pumps collectively align through feedback between ion transport and electrostatic interactions. In the model, directional ion transport generates a membrane potential, while the resulting electrochemical potential biases pump orientation, leading to self-organized collective alignment. Numerical simulations and mean-field analysis reveal a nonequilibrium transition from a disordered state without net transport to a pump-alignment state with sustained membrane potentials. The critical behavior is consistent with the mean-field Ising universality class; however, the effective field is generated self-consistently by nonequilibrium ion transport. We further show that protocell asymmetry can bias the polarity of the membrane potential. These results provide a generic self-organizing mechanism for the emergence of bioelectricity and a physical route toward chemiosmotic coupling in protocells.

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

1 major / 2 minor

Summary. The manuscript proposes a minimal physical model in which ion pumps in a membrane collectively align through a self-consistent feedback loop: directional ion transport generates a membrane potential, which electrostatically biases the orientation of the pumps, leading to sustained net transport. Numerical simulations and mean-field analysis demonstrate a nonequilibrium transition from a disordered state (no net transport) to an aligned state with stable membrane potentials, with critical behavior consistent with the mean-field Ising universality class. The model further examines the role of protocell asymmetry in determining polarity.

Significance. If the proposed mechanism holds, it provides a generic, physically grounded route to the emergence of organized bioelectricity and chemiosmotic coupling in protocells, without requiring pre-existing complex machinery. The mapping to Ising-like criticality offers testable predictions for collective behavior in membrane systems. The self-consistent generation of the effective field by nonequilibrium transport is a notable feature that reduces reliance on external parameters.

major comments (1)
  1. [Abstract and model description] Abstract and model description paragraph: The positive-feedback bias from the electrochemical potential to individual pump orientations is introduced by assumption rather than derived from the underlying Nernst-Planck or pump-cycle kinetics. This assumption is load-bearing for the self-consistent loop that produces the alignment transition; if the microscopic bias is weak, opposite in sign, or screened by local gradients, the disordered state may remain stable. An independent derivation or explicit sensitivity analysis to the bias functional would be required to support the central claim.
minor comments (2)
  1. [Numerical simulations] Provide explicit details on parameter choices, data exclusion criteria, and convergence checks for the numerical simulations to allow independent verification of the reported transition.
  2. [Mean-field analysis] Clarify the precise form of the self-consistent effective aligning field in the mean-field equations, including any mean-field approximations or closures used.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address the major comment below.

read point-by-point responses

  1. Referee: [Abstract and model description] Abstract and model description paragraph: The positive-feedback bias from the electrochemical potential to individual pump orientations is introduced by assumption rather than derived from the underlying Nernst-Planck or pump-cycle kinetics. This assumption is load-bearing for the self-consistent loop that produces the alignment transition; if the microscopic bias is weak, opposite in sign, or screened by local gradients, the disordered state may remain stable. An independent derivation or explicit sensitivity analysis to the bias functional would be required to support the central claim.

    Authors: We agree that the orientational bias is introduced phenomenologically in this minimal model rather than being derived from a detailed treatment of Nernst-Planck transport or the full pump reaction cycle. This choice is deliberate to isolate the self-consistent feedback loop as the essential physical mechanism. The functional form is motivated by the electrostatic torque that a transmembrane field exerts on an asymmetrically charged membrane protein, which generically favors orientations aligned with the field. To address the referee's concern about robustness, we will add an explicit sensitivity analysis in the revised manuscript. This analysis will vary both the magnitude of the bias coefficient and its functional dependence on the local electrochemical potential, demonstrating that the nonequilibrium alignment transition persists over a broad range of positive-feedback strengths. We note that a first-principles derivation from microscopic kinetics would require specifying the molecular charge distribution and conformational states of the pumps, which lies outside the scope of the present minimal model. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained within proposed model

full rationale

The paper explicitly proposes a minimal model in which the electrochemical bias on pump orientation is introduced as part of the model definition to enable positive feedback. Numerical simulations and mean-field analysis are then used to derive the nonequilibrium transition to collective alignment and sustained membrane potentials. This emergent transition is not equivalent by construction to the input assumptions; it is a dynamical consequence obtained from solving the model equations, analogous to standard self-consistent mean-field treatments. No load-bearing self-citations, fitted parameters renamed as predictions, or reductions of the central claim to tautological inputs are present. The analysis remains self-contained against the model's own dynamics and does not rely on external unverified premises for its core result.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The model rests on standard biophysical assumptions about ion pumps and membranes plus the specific feedback rule introduced in the paper; no new particles or forces are postulated.

axioms (1)
  • domain assumption The electrochemical potential generated by ion transport biases the orientation of individual pumps.
    This feedback rule is the core of the self-organization mechanism described in the abstract.

pith-pipeline@v0.9.0 · 5704 in / 1159 out tokens · 39264 ms · 2026-05-21T13:28:23.008346+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

46 extracted references · 46 canonical work pages

  1. [1]

    Schoepp-Cothenet, R

    B. Schoepp-Cothenet, R. van Lis, A, Atteia, F. Baymann, L. Capowiez, A.-L. Ducluzeau, S. Duval, F. ten Brink, M. J. Russell, and W. Nitschke, On the universal core of bioenergetics, Biophysica Acta (BBA)-Bioenergetics, 2, 79–93 (2013)

    work page 2013

  2. [2]

    Mitchell, Coupling of phosphorylation to electron and hydrogen transfer by a chemi-osmotic type of mechanism, Nature, 191, 144–148 (1961)

    P. Mitchell, Coupling of phosphorylation to electron and hydrogen transfer by a chemi-osmotic type of mechanism, Nature, 191, 144–148 (1961)

    work page 1961

  3. [3]

    Mitchell, Chemiosmotic coupling in oxidative and photosynthetic phosphorylation, Biological Reviews, 41, 445–501 (1966)

    P. Mitchell, Chemiosmotic coupling in oxidative and photosynthetic phosphorylation, Biological Reviews, 41, 445–501 (1966)

    work page 1966

  4. [4]

    A. M. Morelli, S. Ravera, D. Calzia, and I. Panfoli, An update of the chemiosmotic theory as suggested by pos- sible proton currents inside the coupling membrane, 9, 180221 (2019)

    work page 2019

  5. [5]

    M. C. Weiss, F. L. Sousa, N. Mrnjavac, S. Neukirchen, M. Roettger, S. Nelson-Sathi, and W. F. Martin, The physi- ology and habitat of the last universal common ancestor, Nature microbiology, 1, 1–8 (2016)

    work page 2016

  6. [6]

    E. R. R. Moody, S. Álvarez-Carretero, T. A. Mahen- drarajah, J. W. Clark, H. C. Betts, N. Dombrowski, L. L. Szánthó, R. A. Boyle, S. Daines, X. Chen, N. Lane, Z. Yang, G. A. Shields, G. J. Szöllősi, A. Spang, D. Pisani, T. A. Williams, T. M. Lenton, and P. C. J. Donoghue, The nature of the last universal common ancestor and its impact on the early Earth s...

    work page 2024

  7. [7]

    V. Sojo, A. Pomiankowsli, and N. Lane, A bioenergetic basis for membrane divergence in archaea and bacteria, PLoS biology, 12, e1001926 (2014)

    work page 2014

  8. [8]

    F. Yu, J. Fei, Y. Jia, T. Wang, W. F. Martin, and J. Li, Chemiosmotic ATP synthesis by minimal protocells, Cell Reports Physical Science, 6 (2025)

    work page 2025

  9. [9]

    Dyson, Origins of life, Cambridge University Press (1999)

    F. Dyson, Origins of life, Cambridge University Press (1999)

    work page 1999

  10. [10]

    Jain, and S

    S. Jain, and S. Krishna, A model for the emergence of cooperation, interdependence, and structure in evolving networks, Proceedings of the National Academy of Sci- ences, 98, 543–547 (2001)

    work page 2001

  11. [11]

    Blokhuis, D

    A. Blokhuis, D. Lacoste, and P. Nghe, Universal motifs and the diversity of autocatalytic systems, Proceedings of the National Academy of Sciences, 117, 25230–25236 (2020)

    work page 2020

  12. [12]

    K. Kaneko, On recursive production and evolvability of cells: catalytic reaction network approach, Geometric Structures of Phase Space in Multidimensional Chaos: Applications to Chemical Reaction Dynamics in Com- plex Systems, 130, 543–598 (2005)

    work page 2005

  13. [13]

    S. A. Kauffman, The origins of order: Self-organization and selection in evolution, In Spin glasses and biology, 61–100 (1992)

    work page 1992

  14. [14]

    Eigen, Selforganization of matter and the evolution of biological macromolecules, Naturwissenschaften, 58, 465–523 (1971)

    M. Eigen, Selforganization of matter and the evolution of biological macromolecules, Naturwissenschaften, 58, 465–523 (1971)

    work page 1971

  15. [15]

    Eigen and P

    M. Eigen and P. Schuster, The hypercycle: a principle of natural self-organization, Springer Science & Business Media (2012)

    work page 2012

  16. [16]

    Kaneko and T

    K. Kaneko and T. Yomo, On a kinetic origin of heredity: minority control in a replicating system with mutually catalytic molecules, Journal of theoretical biology, 214, 563–576 (2002)

    work page 2002

  17. [17]

    A. V. Tkachenko and S. Maslov, Spontaneous emergence of autocatalytic information-coding polymers, The Jour- nal of Chemical Physics, 143 (2015)

    work page 2015

  18. [18]

    M. J. Falk, L. Zhou, Y. J. Matsubara, K. Husain, J. W. Szostak, and A. Murugan, Suppression of errors in collectively coded information, arXiv:2508.21806

    work page arXiv

  19. [19]

    Segré, D

    D. Segré, D. Ben-Eli, D.W. Deamer, and D. Lancet, The lipid world, Origins of Life and Evolution of the Bio- sphere, 31, 119–145 (2001)

    work page 2001

  20. [20]

    Lancet, R

    D. Lancet, R. Zidovetzki, and O. Markovitch, Systems protobiology: origin of life in lipid catalytic networks, Journal of The Royal Society Interface 15, 20180159 (2018)

    work page 2018

  21. [21]

    Szathmáry, The origin of replicators and reproduc- ers, Philosophical Transactions of the Royal Society B: Biological Sciences, 361, 1761–1776 (2006)

    E. Szathmáry, The origin of replicators and reproduc- ers, Philosophical Transactions of the Royal Society B: Biological Sciences, 361, 1761–1776 (2006)

    work page 2006

  22. [22]

    Goldenfeld and C

    N. Goldenfeld and C. Woese, Life is physics: evolution as a collective phenomenon far from equilibrium, Annu. Rev. Condens. Matter Phys., 2, 375-399 (2011)

    work page 2011

  23. [23]

    Branscomb, T

    E. Branscomb, T. Biancalani, N. Goldenfeld, and M. Russell, Escapement mechanisms and the conversion of disequilibria; the engines of creation, Phys. Rep. 677, 1–60 (2017)

    work page 2017

  24. [24]

    J. P. Morth, B. P. Pedersen, M. J. Buch-Pedersen, J. P. Andersen, B. Vilsen, M. G. Palmgren, and P. Nis- sen, A structural overview of the plasma membrane Na+, K+-ATPase and H+-ATPase ion pumps, Nature reviews Molecular cell biology 12, 60–70 (2011)

    work page 2011

  25. [25]

    M. Dyla, M. Kjærgaard, H. Poulsen, and P. Nissen, P, Structure and mechanism of P-type ATPase ion pumps. Annual review of biochemistry, 89, 583-603 (2020)

    work page 2020

  26. [26]

    Calisto, F

    F. Calisto, F. M. Sousa, F. V. Sena, P. N. Refojo, and M. M. Pereira, Mechanisms of energy transduction by charge translocating membrane proteins, Chemical re- views, 121, 1804–1844 (2021)

    work page 2021

  27. [27]

    D. A. Beard, A biophysical model of the mitochondrial respiratory system and oxidative phosphorylation, PLoS computational biology 1, e36 (2005)

    work page 2005

  28. [28]

    Terradot, E

    G. Terradot, E. Krasnopeeva, P. S. Swain, and T. Pili- zota, Escherichia coli Maintains pH via the Membrane Potential, Physical Review X, 2, 043015 (2024)

    work page 2024

  29. [29]

    W. Lo, E. Krasnopeeva, and T. Pilizota, Bacterial elec- trophysiology, Annual Review of Biophysics, 53 (2024)

    work page 2024

  30. [30]

    K. D. Garlid, A. D. Beavis, and S. K. Ratkje, On the nature of ion leaks in energy-transducing membranes, Biochimica et Biophysica Acta (BBA)-Bioenergetics, 12, 109–120 (1989)

    work page 1989

  31. [31]

    Kondepudi and I

    D. Kondepudi and I. Prigogine, Modern thermodynam- ics: from heat engines to dissipative structures, John wi- 6 ley & sons (2014)

    work page 2014

  32. [32]

    Von Heijne, Membrane-protein topology, Nature re- views Molecular cell biology, 7, 909–918 (2006)

    G. Von Heijne, Membrane-protein topology, Nature re- views Molecular cell biology, 7, 909–918 (2006)

    work page 2006

  33. [33]

    Bogdanov, H

    M. Bogdanov, H. Vitrac, and W. Dowhan, Flip-flopping membrane proteins: how the charge balance rule governs dynamic membrane protein topology, Biogenesis of Fatty Acids, Lipids and Membranes, 1–28 (2018)

    work page 2018

  34. [34]

    See Supplemental Material for details of simulation method, calculation, and supplemental figures

    work page

  35. [35]

    For simplicity, we assume that ion transport is dominated by pumps, and neglect ion channels and leakage

    Ions are passively transported via membrane-embedded ion channels and leakage. For simplicity, we assume that ion transport is dominated by pumps, and neglect ion channels and leakage. While channels and leakage act to reduce ion concentration gradients, they do not qualita- tively affect our results

    work page

  36. [36]

    β, z, and e are inverse temperature, the valence of ion, and the elementary charge, respectively

    The reference chemical potential µ0 is common to both regions and can be absorbed into the kinetic prefactor. β, z, and e are inverse temperature, the valence of ion, and the elementary charge, respectively

    work page

  37. [37]

    For simplicity, here we neglect the backward reaction

    work page

  38. [38]

    For simplicity, we set γ to be constant

    work page

  39. [39]

    q = 0 corresponds to cI = cE = c0, indicating the ab- sence of a membrane potential

    work page

  40. [40]

    The orientation can reverse in response to alterations in the charge distribution of residues[ 33]

    The orientation of membrane proteins is characterized by the positive-inside rule[ 32], whereby residues with posi- tive charge face the cytoplasmic side, reflecting differ- ences in charge and electrostatic potential across the membrane. The orientation can reverse in response to alterations in the charge distribution of residues[ 33]

    work page

  41. [41]

    Since the pumps align bidirectionally, the state is char- acterized by the absolute value of m

    work page

  42. [42]

    Here we use a backward pumping rate, γb = 10 −5

    For small rV , q can take large values ( −1 ≤ q ≤ 1/rV ), so that the backward reaction becomes non-negligible. Here we use a backward pumping rate, γb = 10 −5

    work page

  43. [43]

    ( 14) represents a generalization of Eq

    Eq. ( 14) represents a generalization of Eq. ( 9)

    work page

  44. [44]

    This occurs when the chemical reactions driving the pumps are asymmetric between the interior and exterior

    work page

  45. [45]

    This results in a shift from the origin with a contribution of (ln rγ)/2, as in Eq

    rγ represents the ratio of the pumping strength of the interior to that of the exterior, where rγ ̸= 1 leads to a steady state ¯q ̸= 0. This results in a shift from the origin with a contribution of (ln rγ)/2, as in Eq. ( 40). In con- trast, the term with rV , even for rV ̸= 1, always keeps a steady-state solution ¯q = 0, which induces concentration- depe...

    work page

  46. [46]

    END MA TTER Membrane potential The membrane potential ∆ϕ arises from an imbalance of charges across the membrane, which has low ion per- meability and acts as a capacitor[ 29]

    Note that we here consider the large NP limit, treat p as a continuous variable, and use the mean-field value of q. END MA TTER Membrane potential The membrane potential ∆ϕ arises from an imbalance of charges across the membrane, which has low ion per- meability and acts as a capacitor[ 29]. It is defined as the difference between the net charges inside, ...

    work page