
==== Front
J Phys Chem B
J Phys Chem B
jp
jpcbfk
The Journal of Physical Chemistry. B
1520-6106
1520-5207
American Chemical Society

39167050
10.1021/acs.jpcb.4c03047
Article
Oxidative Phosphorylation Does Not Violate the Second Law of Thermodynamics
https://orcid.org/0000-0003-3349-2876
Silverstein Todd P. *
Department of Chemistry (emeritus), Willamette University, Salem, Oregon 97301,United States
* Email: tsilvers@willamette.edu.
21 08 2024
05 09 2024
128 35 84488458
08 05 2024
15 08 2024
03 07 2024
© 2024 The Author. Published by American Chemical Society
2024
The Author
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

In a recent series of papers, James W. Lee reported that mitochondrial oxidative phosphorylation violates the second law of thermodynamics and that it is allowed to do so because it is a “Type-B” process that features lateral and longitudinal membrane asymmetry. We show here that these contentions are based on problematic interpretations of the literature. More reliable values of ΔGredox and ΔGATP synthesis show that the second law is not violated. More recent reports on the structures of the redox-driven proton pumps (Complexes I, III, and IV) suggest that longitudinal membrane asymmetry does not exist. Finally, Lee’s predictions for the concentration of protons localized at the P-side surface of the bioenergetic membrane are likely to be much too high due to several errors; thus, his predicted high values of ΔpHsurface that violate the second law are likely to be wrong. There is currently no strong experimental or theoretical evidence to support the contention that oxidative phosphorylation violates the second law of thermodynamics.

document-id-old-9jp4c03047
document-id-new-14jp4c03047
ccc-price
==== Body
pmcIntroduction

In a recent series of papers, James Lee reported that oxidative phosphorylation in mitochondria violates the second law of thermodynamics in two distinct ways.1−5 Here, we wish to examine carefully the evidence and assumptions underlying this bold claim. We begin with a brief introduction to oxidative phosphorylation, which is currently understood in light of the chemiosmotic theory, for which Peter Mitchell won the 1978 Nobel Prize.6 Chemiosmotic ATP synthesis is driven by a transmembrane difference of the proton electrochemical potential, ΔμH+. Initially, ΔμH+ is established by redox-driven proton pumping by the electron transfer chain coupling sites. There are three thermodynamic components to this system: (i) energy input from spontaneous redox reactions catalyzed by the coupling sites (e.g., oxidation of NADH, FADH2); (ii) a high-energy intermediate, the proton electrochemical gradient, ΔμH+; and (iii) energy output via proton-driven ATP synthesis.

ΔμH+ comprises an electrical component, the transmembrane potential difference (Δψ), and a chemical component, the proton concentration gradient (ΔpH), as given by eq 1:1

For our purposes, we will define the proton transport reaction in eq 1 as import from the outside (P) to the inside (N):7 Δψ ≡ ψN – ψP is invariably inside-negative, while ΔpH ≡ pHN – pHP is generally positive (i.e., inside-alkaline, except for alkaliphilic bacteria). Thus, ΔμH+(P → N) is invariably negative, proton import is spontaneous, and the free energy from proton influx can be used to drive nonspontaneous ATP synthesis.

Mitchell and many of the early proponents of the chemiosmotic theory insisted that the relevant proton concentrations were in the bulk aqueous phases on either side of the bioenergetic membrane: ΔpHbulk ≡ pHN,bulk – pHP,bulk. Studies on chloroplast thylakoid membranes, whose ATP synthase utilizes a high H+/ATP coupling ratio (4.67 = 14/3), supported this contention.8 However, as early as 1979, questions arose as to whether the energy in the bulk phase ΔpH was sufficient to drive ATP synthesis.9,10 Studies on alkaliphilic bacteria, which thrive at external pH values well above their cytoplasmic pH (making ΔpHbulk negative, thus detracting from ΔμH+), showed that oxidative phosphorylation readily occurred even when the bulk phase ΔμH+ was insufficient to drive ATP synthesis.11−14 Similar concerns were raised regarding low-potential bacteria and mitochondria, which utilize lower H+/ATP coupling factors of 3.3 and 2.7.8 Recently, high-precision narrowly focused measurements of mitochondrial pH have shown that ΔpHbulk across the mitochondrial F1F0 ATP synthase (≈0.07 to 0.32 pH units15−21) is even lower than previously reported (0.85–1.1 pH units). It has been suggested that this low ΔpH may yield a ΔμH+ that is insufficient to drive mitochondrial ATP synthesis.15,22

To explain these thermodynamically problematic results, many bioenergeticists hypothesized that proton concentration at the outer (P) surface of the bioenergetic membrane exceeded that in the bulk phase; in other words, the bulk and surface phases were maintained out of equilibrium, because protons localized at the membrane surface were not easily released into the bulk phase. This hypothesis (ΔpHsurface > ΔpHbulk) has been dubbed “localized chemiosmosis”. The lack of equilibration has been explained by two quite different models: Lee’s protonic capacitor model1,23 and the potential well/barrier model.24−27

In the past few years, James W. Lee has published a series of papers1−5 suggesting that oxidative phosphorylation violates the second law of thermodynamics in two different ways: (1) The energy required for ATP synthesis (final output) exceeds that supplied by redox reactions (initial input, dotted red B1 line in Figure 1), and also, (2) the energy in the surface-localized high-energy intermediate ΔμH+,TELP exceeds the redox energy input (dotted red B2 line in Figure 1). According to Lee, because of asymmetric structural aspects of the bioenergetic membrane, environmental heat energy can be harnessed to help synthesize ATP. Lee refers to this as a “thermotrophic” feature of a “type-B” process, which is not constrained by the second law.

Figure 1 Free energy profile of redox-driven proton pumping and proton-driven ATP synthesis. L.C. = localized chemiosmosis: the bulk ΔμH+ is insufficient to drive ATP synthesis, requiring a higher, localized ΔpHsurface. NADH and FADH2 oxidation yield ≈ −4.6 kcal/mol H+, which is identical to the energy required for proton-flow driven ATP synthesis in the matrix. B = type-B processes that violate the 2nd law of thermodynamics because the redox energy input is exceeded by (B1) ATP synthesis energy (output) and/or (B2) ΔμH+,TELP due to high ΔpHsurface established by a “protonic capacitor”.

A free energy profile depicts the three thermodynamic components of oxidative phosphorylation (Figure 1). The left-most reaction, NADH (or FADH2) oxidation, releases a maximum of ≈4.6 kcal per mole of protons pumped (Supporting Information, section IIA).8,22 Protons flowing back across the membrane through the F0 oligomer supply the 4.6 kcal/mol H+ (=12.3 kcal/mol ATP ÷ 2.67 H+/ATP) required for ATP synthesis in the matrix (Supporting Information, section IIB). The overall maximum thermodynamic efficiencya of the process in the matrix is ≈100%; for ATP exported to the cytoplasm, the maximum thermodynamic efficiency is 91% (15.4 kcal/mol ATP ÷ 3.67 H+/ATP ÷ 4.6 kcal/mol H+(redox); Supporting Information, section IIB). All processes obey the second law of thermodynamics: redox energy input ≥ ΔμH+,LC ≥ ATP synthesis energy output.

In the free energy profile (Figure 1), the need for localized chemiosmosis (L.C.) is depicted by the black dashed line segment: Here, the bulk-to-bulk ΔμH+ is insufficient to drive ATP synthesis. Lee’s two type-B processes are depicted by the two red dashed line segments, wherein redox energy input is exceeded by B1, ATP synthesis energy, and/or B2, the localized ΔμH+,TELP.

Claiming that oxidative phosphorylation violates the second law of thermodynamics is a bold contention. Oxidative phosphorylation, which first appeared in bacteria over 2 billion years ago, is a ubiquitous metabolic process in eukaryotic mitochondria and prokaryotic cells. It seems surprising at first glance that such a universal process could violate such a widely embraced thermodynamic law. We have therefore undertaken to examine the evidence supporting Lee’s claim.

Results and Discussion

B1: ATP Synthesis Requires More Energy than Redox Reactions Can Supply

In making this claim, Lee reported2,3,28 ΔμH+ = −5.26 kcal/mol H+ available from redox-driven proton pumping, but +5.85 kcal/mol H+ required for proton-driven ATP synthesis.b Thus, the per proton energy input exceeds output, and the second law of thermodynamicsc is violated. It turns out, however, that both of these numbers are too high. The value of −5.26 kcal/mol H+ available from redox-driven proton pumping is derived from standard redox potentials and thus applies to the reactions carried out when all reactant and product concentrations (except for H+) are 1 M. We have reported previously that under typical steady state concentrations of aqueous NADH, NAD+, FAD, FADH2, and O2, ΔμH+ = −4.45 and −4.68 kcal/mol H+ for oxidation of NADH and FADH2, respectively.8 (See the Supporting Information, section IIA, for further details.)

Lee calculated +5.85 kcal/mol H+ required for proton-driven ATP synthesis by dividing ΔGp = +15.6 kcal/mol ATP from Cockrell et al.29 by the H+/ATP coupling factor of 2.67 for the F1F0-catalyzed ATP synthesis in the mitochondrial matrix. However, Cockrell et al.’s value of ΔGp was not in fact measured in the matrix of normally functioning intact state 3 mitochondria (see the Supporting Information, section IIB). A comprehensive meta-analysis of reported ΔGp values (section IIC, Supporting Information) shows that its value in the mitochondrial matrix, where ATP is synthesized on the F1 complex, is 12.3 ± 0.8 kcal/mol ATP. Dividing this by 2.67 H+/ATP, we get ΔμH+ = +4.6 kcal/mol H+ flowing through the F1F0 ATP synthase. Note that this value is 21% lower than that used by Lee (5.85 kcal/mol H+) and is approximately equal to the available redox free energy (−4.6 kcal/mol H+). Considering ATP that is synthesized in the matrix and exported to the cytoplasm, we have ΔGp,cytopl. = 15.4 kcal/mol ÷ 3.67 H+/ATP = 4.2 kcal/mol H+ (see the Supporting Information, section IIB). Thus, the second law is not violated: Spontaneous redox reactions supply enough free energy to drive ATP synthesis.

B2: The Surface-Localized Proton Gradient Exceeds Redox Energy Input

Calculating Lee’s Surface-Localized Proton Chemical Potential, ΔμH+,P,L

Considering spontaneous proton influx (P → N) through the F0 complex that drives ATP synthesis, Lee distinguished between the chemical potential differences for three distinct proton transport processes (Figure 2): (1) from P bulk to N bulk phase, ΔμH+,bulk(P → N), which is given by Mitchell’s eq 1; (2) from P bulk to the P surface-localized TELP layer, ΔμH+,P,L; and (3 + 4) the total for localized (TELP) chemiosmosis from the P surface to the N bulk phase, ΔμH+,LC(P → N).

Figure 2 Proton transport steps in Lee’s model of transmembrane proton transport.

Since Pbulk → Nbulk transport (ΔμH+,bulk(P → N)) is equivalent to the sum of Pbulk → PTELP (ΔμH+,P,L) and PTELP → Nbulk (ΔμH+,LC(P → N)), we have2

Lee derived the following definition for ΔμH+,P,L:3

According to Lee, the proton concentration in the P surface layer = [H+]P,bulk + [H+]P,L, where [H+]P,L is the concentration of “excess” protons localized in the P surface TELP layer due to the TELP protonic capacitor electrostatic effect. From eq 3, we see that as long as [H+]P,L > 0, ΔμH+,P,L is positive; furthermore, as discussed above, due to redox-driven proton pumping, ΔμH+,bulk(P → N) is negative. Thus, we see from eq 2 that ΔμH+,LC(P → N) is more negative than ΔμH+,bulk(P → N), with the difference determined by the magnitude of ΔμH+,P,L: That is, more free energy is released in localized chemiosmosis than is available in the bulk phase gradient alone.

There are several serious problems with the thermodynamic and structural biological assumptions used by Lee in his derivations of eqs 2 and 3 above. These are discussed in detail in section III of the Supporting Information. However, for the purposes of our discussion here, we will accept the validity of these equations. In order to calculate ΔμH+,P,L (and thus the net ΔμH+,LC(P → N)), Lee used the following electrostatic model to calculate the concentration of “excess” localized protons, [H+]P,L.

Lee’s TELP Protonic Capacitor Equations

Since his first paper30 on this subject in 2012, Lee has characterized the localization of protons at the membrane surface using his transmembrane electrostatically localized protons (TELP) hypothesis. In this model, Lee assumes that a bioenergetic membrane that supports a nonzero transmembrane potential functions like a capacitor in an electrical circuit. Protons, because they diffuse much faster than other ions,d behave like electrons in a circuit. Two TELP surface layers that stretch 1 nm from either side of the low-dielectric membrane interior out toward the aqueous phase serve as the two plates of the “protonic capacitor”. Beyond the 1 nm-thick TELP layer lies the bulk phase. Excess protons at the P surface, and hydroxide anions at the N surface, are held within the TELP layers due to strong electrostatic attraction.e Using electrostatics equations describing electronic capacitors, Lee derived:4

where C = capacitance, A = area, l = thickness of the surface layer, and F = Faraday’s constant. Lee used l = 1 nm, as a “reasonable” thickness, with no explanation, and C/A = 13.2 mF/m2 = 1.32 μF/cm2, which was measured across leukemia cell plasma membranes.32 The multiplying factor, C/(A·l·F), then comes out to be 0.0136(8) M/V.

There are problems with Lee’s selected values for both C/A and l. C/A values are clearly higher for cancer cells (1.2–2.1 μF/cm2) than for normal cells and organelles (0.5–1.2 μF/cm2, Supporting Information, section V). Thus, Lee’s choice of 1.32 μF/cm2 is undoubtedly too high. The value reported33 for neuronal plasma membranes, 0.94 ± 0.20 μF/cm2 (range: 0.8–1.3 μF/cm2), seems like a reasonable choice for bioenergetic membranes. Furthermore, Lee’s choice of l = 1 nm supposedly represents the thickness of the bilayer headgroup region along with a monolayer of surface water at the headgroup/bulk interface (Lee, personal communication). Given that the lipid headgroup region is 1.1 ± 0.2 nm thick,34−37 and proton diffusion at the membrane surface involves 1–2 water monolayers at the headgroup/bulk interface,38,27,26,25 which would be about 0.4 nm thick, a better choice for the thickness of the TELP layer would be 1.5 nm. Using l = 1.5 nm and C/A = 0.94 μF/cm2, the multiplying factor, C/(A·l·F), comes out to be 0.0062(7) M/V, which is less than half of the value used by Lee.

Lee has pointed out that the actual [H+]P,L is lower than [H+]cap due to cations in solution exchanging with protons in the TELP surface layer. If the equilibrium constant for the Miz+ cation replacing protons in the L layer (H+L + Miz+bulk → H+bulk + Miz+L) is Kexchg,i, then considering n different cations in solution, i = 1 to n, Lee derived:5

Lee used literature values for bulk concentrations of Na+, K+, Mg2+, Ca2+, Zn2+, Fe2+, and other cations. He used Kexchg values of 5.07(10–8) for Na+, 6.93(10–8) for K+, ≈6(10–8) for NH4+, and 2.1(10–7) for Mg2+ and all other divalent cations.f Using these values together with literature values of mitochondrial P side bulk concentrations of the cations, Lee calculated the denominator in eq 5, dubbed the “exchange factor,” to be 1.29. There are reasons to believe that this value is orders of magnitude too low (see footnote f and Table 1).

Table 1 Key Thermodynamic Parameters Calculated from TELP Equations, Given pHN,bulk = 7.35 and pHP,bulk = 7.25a

Δψ(mV)	[H+]P,L, (mM)	ΔμH+,bulk(P → N)(kcal/mol H±)	ΔμH+,P,L(kcal/mol H±)	ΔμH+,LC(P → N)(kcal/mol H±)	
–81	8.65	–2.0	7.5	–9.5	
–123	13.1	–3.0	7.8	–10.74	
–159	17.0	–3.8	7.9	–11.7	
–180	19.2	–4.3	8.0	–12.3	
–81*	1.7(10–7)	–2.0	0.144	–2.15	
–123*	2.6(10–7)	–3.0	0.145	–3.1	
–159*	3.3(10–7)	–3.8	0.1455	–3.9(5)	
–180*	3.7(10–7)	–4.3	0.145	–4.4(4)	
a The exchange factor (eq 5 denominator) is 1.29 in the top four rows (following Lee) and 6.6(107) in the bottom four starred rows.

Lee used eq 5 to calculate the concentration of “excess” protons in the surface-localized TELP layer, [H+]P,L, as a function of Δψ, and then, using eqs 2 and 3, he calculated ΔμH+,LC(P → N). The range of Δψ that he employed was generally 50 to 200 mV, especially citing low literature values ranging from 56 to 114 mV. However, in vivo measurements suggesting a very low Δψ (−81 to −120 mV) have been challenged as experimental artifacts (Supporting Information, section IV).

In the Supporting Information (section IV), we present an extensive meta-analysis of reported Δψ values measured in mitochondria. The least negative reliable reported values are −123 mV, the most negative = −180 mV, and the average = −159 ± 16 mV (inside negative; Supporting Information, section IV). Lee selected bulk phase pH values of 7.25 on the P side and 7.35 on the N side; using eqs 2 and 5, we can calculate TELP predictions for key thermodynamic parameters under low-, average-, and high-potential conditions. Using Lee’s eq 5 denominator of 1.29 (Table 1, top four rows), the localized ΔμH+,P,L comprises most (65–79%) of the net ΔμH+,LC(P → N).

Furthermore, as Lee has noted, even at the most negative observed Δψ value, ΔμH+,bulk(P → N) is less negative than the 4.6 kcal/mol H+ required for ATP synthesis; on the other hand, the total ΔμH+,LC(P → N) far exceeds the required ATP synthesis energy even for the least negative Δψ. This is the crux of Lee’s TELP hypothesis, namely, that TELP-predicted values of ΔμH+,LC(P → N) exceed the required ATP synthesis energy for all observed values of Δψ. In his recent papers,1−4 Lee has pointed out that ΔμH+,P,L (and therefore, ΔμH+,LC(P → N)) exceeds the available redox energy input of −4.6 kcal/mol H+; hence, TELP predicts violation of the second law of thermodynamics.

There are two problems with the ΔμH+,P,L and ΔμH+,LC(P → N) values in the top four rows of Table 1: (i) They are predicted values and have not been observed experimentally, and (ii) they are based on an ion exchange factor of 1.29, which in turn is based on ion exchange Keq values that have been shown to be thermodynamically untenable40,39 and are undoubtedly orders of magnitude too low.

One obvious problem with ion exchange Keq values of ≈10–8 is that the sizes of the hydrated cations do not differ greatly: radii in Å = 2.80/H+, 3.31/K+, 3.58/Na+, 4.12/Ca2+, 4.28/Mg2+.41 Considering only electrostatic forces, one would expect membrane surface affinity from the strongest to weakest to be H+> Ca2+ > Mg2+ > K+ > Na+. However, based on ionic size and charge, the differences should be less than 2 orders of magnitude, not the 7 or 8 orders of magnitude that Lee uses. Indeed, using more reasonable Kexchg values of 0.001 for Na+, 0.0015 for K+, and 0.006 for Mg2+, the denominator in eq 5, 6.6(107), is more than 50 million times larger than 1.29. This in turn gives predicted [H+]P,L values 50 million-fold smaller than those calculated by Lee (Table 1). The minuscule localized ΔμH+,P,L = 0.14 kcal/mol yield ΔμH+,LC(P → N) values that are barely larger than ΔμH+,bulk(P → N); this, in turn, means that Lee’s calculated ΔμH+,LC(P → N) is insufficient to drive ATP synthesis (i.e., ≤4.6 kcal/mol H+).

Thus, all three of the parameters that Lee uses in eqs 4 and 5 are suspect: His C/A value of 1.32 μF/cm2 is too high by ≈45%, his arbitrarily selected TELP layer thickness of 1 nm is too low by ≈33%, and his M+/H+ ion exchange equilibrium constants are several orders of magnitude too low. This makes Lee’s calculated values of [H+]P,L and TELP/surface-localized free energy highly suspect. Using reasonable ion exchange equilibrium constant values of ≈10–3 gives very low (≈0.3 μM) localized [H+]P,L; this yields a localized ΔμH+,LC(P → N) that does not violate the second law of thermodynamics but is also insufficient to drive ATP synthesis. On the other hand, using Saeed and Lee’s untenable ion exchange equilibrium constant values of ≈10–8 gives high (≈10 mM) localized [H+]P,L; this yields a localized ΔμH+,LC(P → N) that is sufficient to drive ATP synthesis but violates the second law. Lee’s response has been to embrace this second law violation, proposing that due to structural and topological asymmetries in the bioenergetic membrane, the establishment of the P surface-localized proton layer is a so-called “Type-B” process that is allowed to violate the second law of thermodynamics.1−4 He claimed, for example, “that transmembrane electrostatically localized excess protons at the liquid water–membrane interface can isothermally utilize their molecular thermal motions to do useful work in driving ATP synthesis at the biophysical molecular scale”.1

TELP Experimental Support

Measurements of Surface pH

Using his TELP hypothesis, Lee has predicted low values of pHP,surface (≈2). We list in Table S3 (Supporting Information, section VI) experimentally measured values of surface pH within 1 nm of the water/hydrophobic interface measured or calculated at eight different interfaces; the median surface pH for all 18 reported values is 5.4 ± 0.4 (±one standard deviation). Based on the range of reported values, we may conclude that the surface pH lies between ≈4.5 and 6.5. This is dramatically less acidic than the pH ≈ 2 values predicted by Lee’s TELP hypothesis; thus, literature reports do not support the predictions of the TELP hypothesis.

Lee has compared his TELP-predicted ΔμH+,LC(P → N) or protonmotive force (pmf) values to two experimental results reported in the literature: the doubling time of alkaliphilic Bacillus pseudofirmus OF4 bacteria13 and the ATP efflux rate through the mitochondrial adenine nucleotide translocase (ANT).18

TELP and Alkaliphilic Bacterial Cell Doubling

Lee supported the bioenergetic significance of his TELP-predicted total protonmotive force (pmf) in alkaliphilic bacteria, writing the following: “the overall pattern of the total pmf [= -ΔμH+,LC(P → N)/F] amazingly matched with the B. pseudofirmus OF4... cell population growth doubling time... This indicates that we really are on something [sic] real”.1 However, we can see from the actual data plotted in in Figure 3 that the cell doubling time (black diamonds) remains stable (at 30–40 min) up to an external pH of 10.6; above that pH, doubling time rises exponentially toward infinity (i.e., no measurable doubling time above pH ≈ 11.5). On the other hand, total pmf (blue circles) is stable (at 470 mV), only up to external pH 8.5; above that pH, total pmf declines gradually and then linearly (pH ≥ 10.5). Thus, Lee’s claim of an “amazing” match between the pH-dependence of total pmf and cell doubling time is not supported.

Figure 3 TELP-predicted total pmf (bulk + localized; blue circles, left-hand y-axis), pmf required for ATP synthesis (red squares, left-hand y-axis), and B. pseudofirmus OF4 cell doubling time (black diamonds, right-hand y-axis) as a function of external pH. Calculated total pmf from ref (1); ADP phosphorylation potential (in mV) and OF4 cell doubling time from ref (13). Linear decline in total pmf above pH 10.5: slope = −174 mV/pH unit; intercept = 2123 mV; R2 = 0.992.

Another problematic observation is that at pH 10.8, the cell doubling time is more than twice that at pH 10.6, and yet the TELP-calculated total pmf is more than twice that required for ATP synthesis (red squares, Figure 3). In fact, TELP-calculated total pmf exceeds that required for ATP synthesis up to the highest pH tested. Thus, the relationship between the TELP-calculated total pmf and that required for ATP synthesis is completely decoupled from the cell doubling time. In other words, the dependence of cell doubling time on external pH does not support the bioenergetic significance of the TELP total pmf.

It is interesting to note that regarding cell doubling, the rate (e.g., per hour) is a more appropriate parameter to consider than the time. Also, one would expect the cell doubling rate to be more sensitive to cytoplasmic pH than to external pH. Accordingly, Figure 4 shows that the cell doubling rate falls with the deprotonation of three (n = 2.9 ± 0.4) titratable acidic groups of pKa = 8.51 ± 0.03. We have shown previously22 that the ATP synthesis rate in starved/respiration-inhibited Bacillus firmus OF4 cells12 energized solely by K+/valinomycin-induced Δψ ≈ −180 mV also falls with the deprotonation of titratable groups, pKa = 8.6 to 8.8. Guffanti and Krulwich hypothesized12 that this pH dependence of ATP synthesis rate was due to a protonation-dependent gating residue. Perhaps similar (or identical?) gating residues control cell doubling rate as well.

Figure 4 B. pseudofirmus OF4 cell doubling rate as a function of cytoplasmic pH. Data are adapted from ref (13). Points are fit by nonlinear regression (dotted curve) to the equation for a pH titration; best-fit values are pKa = 8.51 ± 0.03; maximum rate = 1.95 ± 0.07 h–1; n = 2.9 ± 0.4; and R2 = 0.991.

TELP and Mitochondrial ATP Synthesis and Efflux

Lee supported the bioenergetic significance of his TELP-predicted total pmf in mitochondria, writing the following: “the observed pattern of the ATP efflux rate [through the ANT], which decreases as mitochondrial membrane potential (Δψ) is reduced, is generally also in agreement with the pattern of the total pmf [predicted from TELP equations]”.3 Although it is true that both the TELP-predicted total pmf and the ANT-mediated ATP efflux rate (Figure 5) decline roughly linearly as Δψ increases (i.e., gets less negative), there are some important differences.

Figure 5 ATP-ADP steady-state exchange rate mediated by the ANT with Δψ clamped at various matrix pH values (6.8–7.9). Data from Figure 7C of ref (18). For Δψ ≥ −104 mV, the efflux rate falls to a constant value of 0.05 ± 0.04 mM ATP/min/mg protein; for Δψ ≤ −104 mV, the efflux rate declines linearly with increasing (i.e., less negative) Δψ: slope = −6.9 ± 0.8 μM ATP/min/mg protein per mV; x-intercept = −88 ± 18 mV; R2 = 0.66.

There are two different types of mitochondrial ATP synthesis: oxidative phosphorylation is Δψ-dependent, whereas substrate level phosphorylation is Δψ-independent. As pointed out by Lee, substrate level phosphorylation generally supplies <10% of ATP under optimal aerobic conditions. This can be seen in Figure 5, where the blue dotted line at 0.05 ± 0.04 mM ATP/min/mg protein shows the substrate level phosphorylation at low potential (Δψ ≥ −104 mV, six points), and the red dotted line shows oxidative phosphorylation increasing linearly as Δψ gets more negative, from −105 to −155 mV. The important thing to note is that oxidative phosphorylation reaches zero at −95 mV, where the blue and red dotted lines intersect in Figure 5. However, Lee’s Figure 1 shows that the TELP-predicted total pmf greatly exceeds the pmf required for ATP synthesis even at Δψ = −50 mV. Thus, according to Lee’s TELP, ATP synthesis via oxidative phosphorylation could occur at Δψ = −50 mV, but the baseline ATP efflux rate for Δψ ≥ −95 mV reported by Chinopoulos et al.18 shows that oxidative phosphorylation ceases when Δψ is above −95 mV.

ΔG, ΔS, and ΔH for Forming the TELP Surface Layer

In our discussion above, we introduced Lee’s thermodynamic parameter ΔμH+,P,L ≡ 2.3RTlog(), which covered proton transport from the P, bulk phase to the P surface-TELP phase; this parameter is positive for [H+]P,L > 0. Lee referred to this as ΔGL, but because his reported ΔGL values are negative,2 it is clear that he changed the transport direction: Lee’s negative ΔGL is for spontaneous proton transport from the P surface-TELP phase (high [H+]P,L) to the P, bulk phase (low [H+]P,bulk).

As with any other chemical process, we can break up free energy into enthalpy and entropy (eq 6):6

In a normal concentration cell, only solute concentration differs between the two phases; solvation is identical. Thus, we can approximate ΔHL ≈ 0, and7

So for the “L” reaction, proton diffusion from the P surface to bulk phase (high to low concentration), ΔSL will be positive and can be calculated from the equation8

Unfortunately, Lee again reversed his signs here, writing ΔGL ≈ TΔSL and calculating negative values for ΔSL. He wrote,2 “This study now, for the first time, numerically shows that transmembrane electrostatic proton localization (Type-B process) represents a negative entropy event with a local protonic entropy change (ΔSL) in a range from −95 to −110 J/K·mol.” Thus, using standard concentration cell thermodynamics, Lee has calculated an estimate of the entropy loss for proton transport from the low concentration P, bulk phase to the high concentration P, surface-TELP phase.

However, we must recall that the approximation sign in eq 8 stems from our assumption that ΔHL ≈ 0, i.e., proton hydration is essentially identical in the surface and bulk phases. But recall that the structure and dielectric permittivity of the surface water layer differs dramatically from that in the bulk phase. Hence, ΔHL ≈ 0 is a poor assumption. The parameter that Lee calculates as ΔSL, −ΔGL/T, is actually not a purely entropic term; rather, we can derive from eq 6 that it is actually ΔSL – ΔHL/T (eq 9):9

Given the high dielectric permittivity of bulk phase water (ε ≈ 80) relative to membrane surface water (ε ≈ 10),25,27,38,42 one would expect proton diffusion from bulk to surface to be endothermic. Thus, considering eq 9, the −ΔGL/T values tabulated by Lee are most likely more negative than the actual ΔSL for the proton diffusion process.

Redox-Attainable P Surface pH Values

Localized chemiosmosis for spontaneous proton flow through F1F0 to drive ATP synthesis can be described by eq 10, which is a modification of Mitchell’s original bulk phase equation:10

Here, Δψ ≡ ψN – ψP is negative-inside, and ΔpHLC, defined as pHN,surface – pHP,surface, is generally positive. Equation 10 can be rearranged to give11

Assuming that ΔμH+,LC(P → N) is equivalent to that available from the redox energy input = −4.6 kcal/mol H+, and using pHN = 7.41 (adjacent to the F1 complex) as reported by Rieger et al.,15 and T = 37 °C, eq 11 becomes12

Using eq 12, we can calculate the minimum pH at the P surface that is attainable from redox energy input as a function of the Δψ maintained. These minimum attainable pHP,surf values range from 7.1 for high-potential mitochondria (−180 mV) to 6.2 for low-potential mitochondria (−123 mV, see the Supporting Information, section IV). In agreement with these calculated pH minima, Cherepanov et al.42 predicted that attaining pHP,surf. below 6–6.5 would be highly unlikely (A. Mulkidjanian, person communication).

Type-B Processes and the Evidence for Membrane Asymmetry

Lee has proposed that oxidative phosphorylation violates the second law in two different ways: B1, in which redox energy input is exceeded by the requisite ATP synthesis energy output; and B2, in which redox energy input is exceeded by the energy in the TELP surface-localized proton gradient. We showed in the first section above that Lee’s B1 conclusion is based on ATP synthesis energy calculations that are too high due to two different errors. Correcting these errors shows that there is no B1-type violation of the second law: Redox energy supplies 4.6 kcal/mol H+, and ATP synthesis requires 4.6 kcal/mol H+ in the matrix or 4.2 kcal/mol H+ in the cytoplasm.

On the other hand, Lee’s TELP-predicted gradients (pHP,L ≈ 2) do present a B2-type violation of the second law, as long as one accepts his thermodynamically untenable ion exchange equilibrium constants (Table 1). As we mentioned above, Lee has argued in his recent papers that oxidative phosphorylation is an isothermal engine that can use thermal motion to do work, i.e., synthesize ATP: The establishment of the TELP P-surface localized proton layer is a type-B process that, due to structural and topological asymmetries in the bioenergetic membrane, is allowed to violate the second law.1−4

The literature on Type-B processes is an interesting subfield of physics and chemistry. Its critics claim that the second law of thermodynamics is a universal law, and experiments that demonstrate Type-B processes that seem to violate the law are either poorly designed or poorly carried out.43 On the other hand, supporters claim that their results are being censored from the peer-reviewed literature due to a superficial bias in favor of the universality of the second law, even though physicists in the 19th century already understood that the second law is not in fact universal.44,45 This is a fascinating controversy, but its details lie well beyond the scope of this paper. Here, we wish to examine the evidence in support of Lee’s contention that bioenergetic membrane asymmetry supports a type-B process.

Lee has proposed two distinct types of type-B supportive membrane asymmetries: lateral and longitudinal. Lateral asymmetry in the plane of the mitochondrial inner membrane has been strongly supported in the literature for 20 years or more. Three distinct regions of the mitochondrial inner membrane have been identified: highly curved cristae rims, cristae flat membranes, and flat inner boundary membranes (IBM, adjacent to the intermembrane space and across from the mitochondrial outer membrane). These regions feature distinct protein composition, Δψ, and ΔpH, which are summarized in Table 2.

Table 2 Examples of Inhomogeneity in the Mitochondrial Inner Membrane and Cristae

 	Cristae rims	Cristae flats	IBM	ref.	
enriched proteins	F1F0 dimers (proton sinks)	Complexes I, III, and IV (Proton sources)	 	(46, 47)	
Δψ	 	–165 ± 10 mV	–155 ± 8 mV	(48)	
pHPa	7.28 ± 0.24	6.88 ± 0.15	 	(15)	
 	6.88 ± 0.09	 	(49)	
7.20 ± 0.10	7.13 ± 0.17	6.85 ± 0.18	(17)	
pHNa	7.41 ± 0.20	7.73 ± 0.24	 	(15)	
 	7.78 ± 0.11	(49)	
a These values are measured with proton-sensitive GPF proteins fused to mitochondrial inner membrane proteins (Complexes IV and V). Although the GFP proton inlet sites are close to the membrane surface, it has been argued28,22 that they are not close enough to the headgroup/fatty acid interface to sample the surface layer pH.

The proton sources (redox-driven proton pumps, complexes I, III, and IV) are found primarily in the flat cristae membranes, whereas the proton sinks (dimeric F1F0 ATP synthases) are found primarily at the highly curved cristae rims.47 On the P side, protons are enriched at the cristae flat membranes and at the IBM compared to the cristae rims (ΔpH ≈ 0.35); on the N side, the difference is reversed, with proton enrichment at the cristae rims (ΔpH ≈ −0.35). Interestingly, this leads to a distinct difference in ΔpH, which is quite small at the cristae rims (≈ 0.13, ref (15)), and more substantial along the remainder of the cristae membrane (ΔpH ≈ 0.9). This makes sense when one considers that the F1F0 ATP synthase at the cristae rims consumes protons at the P surface and releases them on the N side.

Back in 2008, Strauss et al. reported50 that ribbons of F1F0 dimers helped to form cristae rims by “enforcing a strong local curvature on the membrane”. This curvature in turn caused an increase in charge density: Using Finite Element Modeling based on the Poisson equation to simulate electric field strength, they reported that for a crista axial ratio of 10:1, “the charge density on the curved membrane surface is up to 3.5 times higher than on planar membranes at the same constant membrane potential”. This would lead to proton enrichment at the crista rims, with pHflat – pHrim ≈ 0.5. Strauss et al. concluded that the purpose of mitochondrial cristae is not only to increase the surface area of inner membrane available to host oxidative phosphorylation complexes but also to focus protons at the F0 proton inlet, increasing ΔpH and the driving force and rate of ATP synthesis.50

Over a decade later, Lee published a paper23 reaching similar conclusions.g Using an unpublished set of equations51 (from an undergraduate physics problem set) to calculate surface charge density along ellipsoid and disk surfaces, Lee calculated that a 10:1 crista axial ratio would yield proton enrichment at the crista rim of 1 pH unit. He did not explain the discrepancy that his result was twice that of Strauss et al. Nevertheless, the evidence supporting lateral asymmetry along the mitochondrial inner membrane and crista membrane seems strong.

Such is not the case for Lee’s contention of longitudinal asymmetry perpendicular to the plane of the membrane. Based on a schematic representation (essentially a cartoon drawing) published by Dudkina et al.52 in 2020, Lee concluded23 that the redox-driven proton pumps (complexes I, III, and IV) release protons into the bulk aqueous phase of the cristae lumen (P side, Figure 6), whereas the proton inlet in the F0 complex of the ATP synthase is right at the surface of the crista rim membrane. This leads to a “longitudinal” asymmetry on the P side between proton release into the bulk phase (path i on the left side of Figure 6) vs uptake from the surface layer (id in Figure 6). According to Lee, this longitudinal asymmetry figures prominently in the ability of surface proton gradients to violate the second law of thermodynamics (type-B process).

Figure 6 Longitudinal asymmetry in the mitochondrial inner membrane. Mitchellian bulk-to-bulk proton transport is symmetric: paths i → ia → ii → iia; Lee H+ transport is longitudinally asymmetric: i → ib → id → ii → iia; most recent structural results suggest that H+ transport is symmetric, along the surface: i → ic → id → ii → iic → iid.

The problem is that this conclusion is based on a cartoon drawing in which proton release and uptake are depicted in the original paper simply by arrows: The actual sites of proton release are not identified; hence, one cannot conclude from this simple cartoon that the pumps release protons on the P-side into the bulk phase as opposed to the membrane surface. Since Dudkina’s 2010 paper, much better structural information on Complex I (NADH-UQ oxidoreductase) has been published.53−55 Figures in these papers make it clear that although the membrane arm of Complex I does protrude a few nm from the P surface of the crista membrane, the actual proton release sites lie almost flush with the membrane surface.

Similarly, although Complex III (UQ-cyt c oxidoreductase) protrudes tens of nanometers from the P surface, the Qo/P ubiquinone binding site is leveled with the edge of the low dielectric interior of the membrane, near the headgroup region. Because Complex III utilizes a Q cycle mechanism of proton pumping, protons are released on the P side when UQH2 is oxidized at the Qo/P site. Although the proton pathway leading from this site to the exterior of the protein is not yet fully characterized, extended H-bonded water chains have been found connecting each UQ binding site to the nearest membrane surface water layer (Jiapeng Zhu, personal communication). Thus, protons are likely to be released into the surface layer, not the bulk phase.

Finally, for Complex IV (cytochrome c oxidase), recent structural results have been reported for the three paths delineated for transmembrane proton pumping: the H-, D-, and K-channels.56−59 Although the protein surface of Complex IV extends several nanometers beyond both of the membrane surfaces, the proton inlet and outlet sites are generally within 0.1–0.5 nm of the lipid headgroup surface (i.e., within 1.2–1.6 nm of the low dielectric fatty acid/headgroup interface; see Supporting Information, section VII). Thus, they lie within the two water monolayers adjacent to the membrane surface.

In summary, Lee provided no evidence beyond a simple 2010 cartoon figure to support his contention of longitudinal asymmetry in proton source/sink machinery on the P side of the mitochondrial crista membrane. Recent high-resolution structural results suggest that such longitudinal asymmetry does not in fact exist: Pumped protons are most likely released into the same P side surface phase from which they are consumed by the F1F0 ATP synthase.60,61 They are also taken up from the same N side surface phase into which they are released by the F1F0 ATP synthase. We conclude from these recent structural results for Complexes I, III, IV, and V that proton flux is longitudinally symmetric, occurring along both membrane surfaces, as depicted by the red pathway arrows in Figure 6: i → ic → id → ii → iic → iid.

Conclusions

Redox energy input to oxidative phosphorylation equals or exceeds ATP synthesis energy output; hence, there is no violation of the second law of thermodynamics. For all mitochondria, from high potential (−180 mV) to low potential (−123 mV), the minimum pHP,surface attainable from redox energy input (pH 7.1 or 6.2, respectively) is sufficient to support ATP synthesis.

Given the physical and structural differences between bulk and surface water, in order to obtain the thermodynamic activity of surface protons, [H+]surface must be multiplied by its activity coefficient, γ(H+)surface, which is currently unknown but is probably significantly different from 1.

Lee’s TELP equations include several errors, including C/A = 1.32 vs 0.9 μF/cm2, surface layer thickness l = 1.0 vs 1.5 nm, and most significantly, ion exchange Keq values ≈ 10–8 vs 10–3. Thus, his predicted total pmf values are likely to be much too high.

Experimental corroboration of TELP predictions is lacking:Reported pHP,surface values range from 4.5 to 6.5; none are close to 2.

The dependence of alkaliphilic bacterial cell doubling rate on pH titrates with a pKa = 8.51 ± 0.03; this sigmoidal decrease with pH differs significantly from the pH dependence of the TELP-predicted total pmf.

The linear decline of mitochondrial ATP efflux rate (and ATP synthesis) with Δψ, to a steady low baseline value at ≥ −95 mV, cannot be explained by dependence on the TELP-predicted total pmf; this latter value exceeds that required for ATP synthesis even at −50 mV.

Lee’s calculation of ΔSL, a “novel” negative entropy event, contains several sign errors; more importantly, the value he calculates is actually −ΔGL/T (= ΔSL – ΔHL/T), which is more negative than ΔSL because ΔHL is positive.

Recent structural advances show that the longitudinal asymmetry claimed by Lee probably does not exist: Proton pumps most likely deliver protons into the same P-side surface layer from which they are drawn by the F1F0 ATP synthase.

There is no evidence that oxidative phosphorylation violates the second law of thermodynamics. Sufficient energy for ATP synthesis can be supplied by an enhanced surface ΔpH, i.e., the difference pHP,surface < pHP,bulk envisioned by standard localized chemiosmosis, which posits a lack of bulk/surface proton equilibration due to a potential barrier between the two phases.22,24

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcb.4c03047.Section I: Lee’s TELP model: implications of the first and second laws of thermodynamics; Section II: energy supplied by redox reactions is enough to drive ATP synthesis: |ΔGredox| ≥ ΔGp; Section III: capacitance cannot be used to calculate surface free energy or pmf; Section IV: meta-analysis (1969–2021) of mitochondrial transmembrane potentials; Section V: literature values for biological membrane-specific capacitance; Section VI: experimental measurements of pH at the membrane surface; Section VII: proton uptake and release sites in Complex IV (cytochrome c oxidase) (PDF)

Supplementary Material

jp4c03047_si_001.pdf

The author declares no competing financial interest.

Acknowledgments

I wish to thank Stephanie Schaertel for her insight into the first and second laws of thermodynamics and for her help with this discussion in the Supporting Information. I also thank James W. Lee for many hours of discussion over the years, and Alexey Nikulov for sharing his thoughts and literature on Type-B processes. No external funding supported this work.

a Thermodynamic efficiencies (100% × ΔGinput/|ΔGoutput|) are calculated here as 100% × ΔGATPsynth/|ΔGredox|. These are the maximum theoretical thermodynamic efficiencies, in that we are ignoring here all detracting processes like proton leak, slip, uncoupling, variable stoichiometry, and other “mechanical” limitations of the enzymes involved.

b Lee reported these values as pmf = 228 mV and 254 mV, respectively, on p. 8 of ref (3).

c For a discussion of potential first law violations, see Section I of the Supporting Information.

d Diffusion coefficients,31 in nm2/ns, are 9.3 for H+, 5.3 for OH–, 2.0 for K+ and Cl–, 1.3 for Na+, 1.2 for HCO3–. Thus, H+ diffuses ≈5 to 8-fold faster than K+, Na+, Cl–, and HCO3–. Is this difference really enough to prioritize protons as the charge carrier in the bioenergetic “capacitor”?

e We have shown22 that the energy of electrostatic attraction is actually insufficient to hold H+/OH– ion pairs at the surface. Over time, entropy drives spontaneous release into the bulk phase, dismantling the protonic capacitor.

f These values have been shown to be thermodynamically untenable39 and are likely to be orders of magnitude too low.

g However, one of Lee’s conclusions, “the formation of cristae creates more mitochondrial inner membrane surface area and thus more protonic capacitance for transmembrane-electrostatically localized proton energy storage”, is bioenergetically irrelevant. Lee’s TELP equations show that [H+]L depends not on C, but on C/A, which is an area-normalized parameter and thus remains the same with and without cristae. Although the number of surface protons increases with cristae, so does the volume of the surface layer, leaving [H+]L, and thus the pmf driving force, unchanged.
==== Refs
References

Lee J. W. Isothermal Environmental Heat Energy Utilization by Transmembrane Electrostatically Localized Protons at the Liquid–Membrane Interface. ACS omega 2020, 5 (28 ), 17385–17395. 10.1021/acsomega.0c01768.32715223
Lee J. W. Energy Renewal: Isothermal Utilization of Environmental Heat Energy with Asymmetric Structures. Entropy 2021, 23 (6 ), 665 10.3390/e23060665.34070431
Lee J. W. Mitochondrial Energetics with Transmembrane Electrostatically Localized Protons: Do We Have a Thermotrophic Feature?. Sci. Rep. 2021, 11 (1 ), 14575 10.1038/s41598-021-93853-x.34272427
Lee J. W. Type-B Energetic Processes and Their Associated Scientific Implications. J. Sci. Explor. 2022, 36 (3 ), 484–495. 10.31275/20222517.
Lee J. W. Thermotrophy Exploratory Study. J. Sci. Explor. 2023, 37 (1 ), 5–16. 10.31275/20232655.
Mitchell P. Nobel Prize in chemistry. https://www.nobelprize.org/prizes/chemistry/1978/ceremony-speech/ (accessed 2023–11–16).
Azzone G. ; Crofts A. ; Dilley R. ; Dutton P. L. ; Junge W. ; Krulwich T. A. ; Rottenberg H. ; Silverstein T. P. ; Skulachev V. P. ; Wikström M. ; et al. Transmembrane Measurements across Bioenergetic Membranes. Biochim. Biophys. Acta - Bioenerg. 1993, 1183 (1 ), 1–3. 10.1016/0005-2728(93)90002-W.
Silverstein T. P. An Exploration of How the Thermodynamic Efficiency of Bioenergetic Membrane Systems Varies with C-Subunit Stoichiometry of F1F0-ATP Synthases. J. Bioenerg. Biomembr. 2014, 46 (3 ), 229–241. 10.1007/s10863-014-9547-y.24706236
Kell D. B. On the Functional Proton Current Pathway of Electron Transport Phosphorylation: An Electrodic View. Biochim. Biophys. Acta - Rev. Bioenerg. 1979, 549 (1 ), 55–99. 10.1016/0304-4173(79)90018-1.
Ferguson S. J. Fully Delocalised Chemiosmotic or Localised Proton Flow Pathways in Energy Coupling?: A Scrutiny of Experimental Evidence. Biochim. Biophys. Acta - Rev. Bioenerg. 1985, 811 (1 ), 47–95. 10.1016/0304-4173(85)90005-9.
Krulwich T. A. ; Guffanti A. A. Proton-Coupled Bioenergetic Processes in Extremely Alkaliphilic Bacteria. J. Bioenerg. Biomembr. 1992, 24 (6 ), 587–599. 10.1007/BF00762351.1334072
Guffanti A. A. ; Krulwich T. A. Features of Apparent Nonchemiosmotic Energization of Oxidative Phosphorylation by Alkaliphilic Bacillus Firmus OF4. J. Biol. Chem. 1992, 267 (14 ), 9580–9588. 10.1016/S0021-9258(19)50130-2.1577797
Sturr M. G. ; Guffanti A. A. ; Krulwich T. A. Growth and Bioenergetics of Alkaliphilic Bacillus Firmus OF4 in Continuous Culture at High pH. J. Bacteriol. 1994, 176 (11 ), 3111–3116. 10.1128/jb.176.11.3111-3116.1994.8195065
Krulwich T. A. ; Gilmour R. ; Hicks D. B. ; Guffanti A. A. ; Ito M. Energetics of Alkaliphilic Bacillus Species: Physiology and Molecules. Adv. Microb. Phsyiol. 1998, 40 , 401–438. 10.1016/S0065-2911(08)60136-8.
Rieger B. ; Arroum T. ; Borowski M. ; Villalta J. ; Busch K. B. Mitochondrial F 1 F O ATP Synthase Determines the Local Proton Motive Force at Cristae Rims. EMBO Rep. 2021, 22 (12 ), e52727 10.15252/embr.202152727.34595823
Rieger B. ; Junge W. ; Busch K. B. Lateral pH Gradient between OXPHOS Complex IV and F 0 F 1 ATP-Synthase in Folded Mitochondrial Membranes. Nat. Commun. 2014, 5 (1 ), 3103 10.1038/ncomms4103.24476986
Toth A. ; Meyrat A. ; Stoldt S. ; Santiago R. ; Wenzel D. ; Jakobs S. ; von Ballmoos C. ; Ott M. Kinetic Coupling of the Respiratory Chain with ATP Synthase, but Not Proton Gradients, Drives ATP Production in Cristae Membranes. Proc. Natl. Acad. Sci. U.S.A. 2020, 117 (5 ), 2412–2421. 10.1073/pnas.1917968117.31964824
Chinopoulos C. ; Vajda S. ; Csanády L. ; Mándi M. ; Mathe K. ; Adam-Vizi V. A Novel Kinetic Assay of Mitochondrial ATP-ADP Exchange Rate Mediated by the ANT. Biophys. J. 2009, 96 (6 ), 2490–2504. 10.1016/j.bpj.2008.12.3915.19289073
Kim N. ; Ripple M. O. ; Springett R. Measurement of the Mitochondrial Membrane Potential and pH Gradient from the Redox Poise of the Hemes of the Bc1 Complex. Biophys. J. 2012, 102 (5 ), 1194–1203. 10.1016/j.bpj.2012.02.003.22404942
Ripple M. O. ; Kim N. ; Springett R. Mammalian Complex I Pumps 4 Protons per 2 Electrons at High and Physiological Proton Motive Force in Living Cells. J. Biol. Chem. 2013, 288 (8 ), 5374–5380. 10.1074/jbc.M112.438945.23306206
Tyrakis P. A. ; Yurkovich M. E. ; Sciacovelli M. ; Papachristou E. K. ; Bridges H. R. ; Gaude E. ; Schreiner A. ; D’Santos C. ; Hirst J. ; Hernandez-Fernaud J. ; Springett R. ; Griffiths J. R. ; Frezza C. Fumarate Hydratase Loss Causes Combined Respiratory Chain Defects. Cell Rep. 2017, 21 (4 ), 1036–1047. 10.1016/j.celrep.2017.09.092.29069586
Silverstein T. P. Is Localized Chemiosmosis Necessary in Mitochondria? Is Lee’s Protonic Capacitor Hypothesis a Valid Model?. Mitochondrial Commun. 2024, 2 , 48–57. 10.1016/j.mitoco.2024.06.001.
Lee J. W. Protonic Capacitor: Elucidating the Biological Significance of Mitochondrial Cristae Formation. Sci. Rep. 2020, 10 (1 ), 10304 10.1038/s41598-020-66203-6.32601276
Knyazev D. G. ; Silverstein T. P. ; Brescia S. ; Maznichenko A. ; Pohl P. A New Theory about Interfacial Proton Diffusion Revisited: The Commonly Accepted Laws of Electrostatics and Diffusion Prevail. Biomolecules 2023, 13 (11 ), 1641 10.3390/biom13111641.38002323
Nguyen T. H. ; Zhang C. ; Weichselbaum E. ; Knyazev D. G. ; Pohl P. ; Carloni P. Interfacial Water Molecules at Biological Membranes: Structural Features and Role for Lateral Proton Diffusion. PLoS One 2018, 13 (2 ), e0193454 10.1371/journal.pone.0193454.29474432
Weichselbaum E. ; Österbauer M. ; Knyazev D. G. ; Batishchev O. V. ; Akimov S. A. ; Hai Nguyen T. ; Zhang C. ; Knör G. ; Agmon N. ; Carloni P. ; Pohl P. ; et al. Origin of Proton Affinity to Membrane/Water Interfaces. Sci. Rep. 2017, 7 (1 ), 4553 10.1038/s41598-017-04675-9.28674402
Mulkidjanian A. Y. ; Heberle J. ; Cherepanov D. A. Protons@ Interfaces: Implications for Biological Energy Conversion. Biochim. Biophys. Acta - Bioenerg. 2006, 1757 (8 ), 913–930. 10.1016/j.bbabio.2006.02.015.
Lee J. W. TELP Theory: Elucidating the Major Observations of Rieger et al. 2021 in Mitochondria. Mitoch. Commun. 2023, 1 , 62–72. 10.1016/j.mitoco.2023.09.001.
Cockrell R. S. ; Harris E. J. ; Pressman B. C. Energetics of Potassium Transport in Mitochondria Induced by Valinomycin. Biochemistry 1966, 5 (7 ), 2326–2335. 10.1021/bi00871a022.5959455
Lee J. W. Proton-Electrostatics Hypothesis for Localized Proton Coupling Bioenergetics. Bioenergetics 2012, 1 (104 ), 1–8. 10.4172/2167-7662.1000104.
Aqion. Table of Diffusion Coefficients, 2020. https://www.aqion.de/site/diffusion-coefficients (accessed 2024–07–02).
Zheng Y. ; Shojaei-Baghini E. ; Wang C. ; Sun Y. Microfluidic Characterization of Specific Membrane Capacitance and Cytoplasm Conductivity of Single Cells. Biosens. Bioelectr. 2013, 42 , 496–502. 10.1016/j.bios.2012.10.081.
Gentet L. J. ; Stuart G. J. ; Clements J. D. Direct Measurement of Specific Membrane Capacitance in Neurons. Biophysical journal 2000, 79 (1 ), 314–320. 10.1016/S0006-3495(00)76293-X.10866957
Lewis B. A. ; Engelman D. M. Lipid Bilayer Thickness Varies Linearly with Acyl Chain Length in Fluid Phosphatidylcholine Vesicles. J. Mol. Biol. 1983, 166 (2 ), 211–217. 10.1016/S0022-2836(83)80007-2.6854644
Lindahl E. ; Edholm O. Mesoscopic Undulations and Thickness Fluctuations in Lipid Bilayers from Molecular Dynamics Simulations. Biophys. J. 2000, 79 (1 ), 426–433. 10.1016/S0006-3495(00)76304-1.10866968
Woodka A. C. ; Butler P. D. ; Porcar L. ; Farago B. ; Nagao M. Lipid Bilayers and Membrane Dynamics: Insight into Thickness Fluctuations. Phys. Rev. Lett. 2012, 109 (5 ), 058102 10.1103/PhysRevLett.109.058102.23006210
Frallicciardi J. ; Melcr J. ; Siginou P. ; Marrink S. J. ; Poolman B. Membrane Thickness, Lipid Phase and Sterol Type Are Determining Factors in the Permeability of Membranes to Small Solutes. Nat. Commun. 2022, 13 (1 ), 1605 10.1038/s41467-022-29272-x.35338137
Tocanne J.-F. ; Teissié J. Ionization of Phospholipids and Phospholipid-Supported Interfacial Lateral Diffusion of Protons in Membrane Model Systems. Biochim. Biophys. Acta - Rev. Biomembr. 1990, 1031 (1 ), 111–142. 10.1016/0304-4157(90)90005-W.
Silverstein T. P. Serious Problems with Saeed & Lee’s “Experimental Determination of Proton-Cation Exchange Equilibrium Constants at Water-Membrane Interface.. Water 2024, 13 , 103–114. 10.14294/WATER.2024.1.
Silverstein T. P. Lee’s″ Transmembrane Electrostatically-Localized Proton″ Model Does NOT Offer a Better Understanding of Neuronal Transmembrane Potentials. J. Neurophysiol. 2023, 130 , 123–127. 10.1152/jn.00173.2023.37314087
Zhong C. ; Deng Y. ; Hu W. ; Qiao J. ; Zhang L. ; Zhang J. A Review of Electrolyte Materials and Compositions for Electrochemical Supercapacitors. Chem. Soc. Rev. 2015, 44 (21 ), 7484–7539. 10.1039/C5CS00303B.26050756
Cherepanov D. A. ; Feniouk B. A. ; Junge W. ; Mulkidjanian A. Y. Low Dielectric Permittivity of Water at the Membrane Interface: Effect on the Energy Coupling Mechanism in Biological Membranes. Biophys. J. 2003, 85 (2 ), 1307–1316. 10.1016/S0006-3495(03)74565-2.12885673
Kostic M. M. The Second Law and Entropy Misconceptions Demystified. Entropy 2020, 22 , 648 10.3390/e22060648.33286420
(accessed 2024–01–15).

Nikulov A. Symmetry Breaking and the Law of Entropy Increase. Symmetry: Cult. Sci. 2005, 16 (1 ), 47–69.
Sheehan D. P. Maxwell Zombies: Conjuring the Thermodynamic Undead: Modern Descendants of Maxwell’s Celebrated Demon Defy Standard Thinking about Heat Engines and May Offer a Path to a New Kind of Energy Generation. Am. Sci. 2018, 106 (4 ), 234–242. 10.1511/2018.106.4.234.
Davies K. M. ; Strauss M. ; Daum B. ; Kief J. H. ; Osiewacz H. D. ; Rycovska A. ; Zickermann V. ; Kühlbrandt W. Macromolecular Organization of ATP Synthase and Complex I in Whole Mitochondria. Proc. Natl. Acad. Sci. U.S.A. 2011, 108 (34 ), 14121–14126. 10.1073/pnas.1103621108.21836051
Kondadi A. K. ; Anand R. ; Reichert A. S. Cristae Membrane Dynamics–a Paradigm Change. Trends Cell Biol. 2020, 30 (12 ), 923–936. 10.1016/j.tcb.2020.08.008.32978040
Wolf D. M. ; Segawa M. ; Kondadi A. K. ; Anand R. ; Bailey S. T. ; Reichert A. S. ; Van Der Bliek A. M. ; Shackelford D. B. ; Liesa M. ; Shirihai O. S. Individual Cristae within the Same Mitochondrion Display Different Membrane Potentials and Are Functionally Independent. EMBO J. 2019, 38 (22 ), e101056 10.15252/embj.2018101056.31609012
Porcelli A. M. ; Ghelli A. ; Zanna C. ; Pinton P. ; Rizzuto R. ; Rugolo M. pH Difference across the Outer Mitochondrial Membrane Measured with a Green Fluorescent Protein Mutant. Biochem. Biophys. Res. Commun. 2005, 326 (4 ), 799–804. 10.1016/j.bbrc.2004.11.105.15607740
Strauss M. ; Hofhaus G. ; Schröder R. R. ; Kühlbrandt W. Dimer Ribbons of ATP Synthase Shape the Inner Mitochondrial Membrane. EMBO J. 2008, 27 (7 ), 1154–1160. 10.1038/emboj.2008.35.18323778
McDonald K. T. Conducting Ellipsoid and Circular Disk, 2002. http://kirkmcd.princeton.edu/examples/ellipsoid.pdf (accessed 2024–07–02).
Dudkina N. V. ; Kouřil R. ; Peters K. ; Braun H.-P. ; Boekema E. J. Structure and Function of Mitochondrial Supercomplexes. Biochim. Biophys. Acta - Bioenerg. 2010, 1797 (6–7 ), 664–670. 10.1016/j.bbabio.2009.12.013.
Baradaran R. ; Berrisford J. M. ; Minhas G. S. ; Sazanov L. A. Crystal Structure of the Entire Respiratory Complex I. Nature 2013, 494 (7438 ), 443–448. 10.1038/nature11871.23417064
Sazanov L. A. A Giant Molecular Proton Pump: Structure and Mechanism of Respiratory Complex I. Nat. Rev. Mol. Cell Biol. 2015, 16 (6 ), 375–388. 10.1038/nrm3997.25991374
Parey K. ; Lasham J. ; Mills D. J. ; Djurabekova A. ; Haapanen O. ; Yoga E. G. ; Xie H. ; Kühlbrandt W. ; Sharma V. ; Zickermann V. ; et al. High-Resolution Structure and Dynamics of Mitochondrial Complex I—Insights into the Proton Pumping Mechanism. Sci. Adv. 2021, 7 (46 ), eabj3221 10.1126/sciadv.abj3221.34767441
Michel H. ; Behr J. ; Harrenga A. ; Kannt A. CYTOCHROME C OXIDASE: Structure and Spectroscopy. Annu. Rev. Biophys. Biomol. Struct. 1998, 27 (1 ), 329–356. 10.1146/annurev.biophys.27.1.329.9646871
Brzezinski P. ; Gennis R. B. Cytochrome c Oxidase: Exciting Progress and Remaining Mysteries. Journal of bioenergetics and biomembranes 2008, 40 (5 ), 521–531. 10.1007/s10863-008-9181-7.18975062
Yoshikawa S. ; Shimada A. Reaction Mechanism of Cytochrome c Oxidase. Chem. Rev. 2015, 115 (4 ), 1936–1989. 10.1021/cr500266a.25603498
Noodleman L. ; Götz A. W. ; Han Du W. G. ; Hunsicker-Wang L. Reaction Pathways, Proton Transfer, and Proton Pumping in Ba3 Class Cytochrome c Oxidase: Perspectives from DFT Quantum Chemistry and Molecular Dynamics. Front. Chem. 2023, 11 , 1186022 10.3389/fchem.2023.1186022.38188931
Sobti M. ; Smits C. ; Wong A. S. ; Ishmukhametov R. ; Stock D. ; Sandin S. ; Stewart A. G. Cryo-EM Structures of the Autoinhibited E. Coli ATP Synthase in Three Rotational States. Elife 2016, 5 , e21598 10.7554/eLife.21598.28001127
Hahn A. ; Vonck J. ; Mills D. J. ; Meier T. ; Kühlbrandt W. Structure, Mechanism, and Regulation of the Chloroplast ATP Synthase. Science 2018, 360 (6389 ), eaat4318 10.1126/science.aat4318.29748256
