Mitochondrial energetics with transmembrane electrostatically localized protons: do we have a thermotrophic feature?

Transmembrane electrostatically localized protons (TELP) theory has been recently recognized as an important addition over the classic Mitchell’s chemiosmosis; thus, the proton motive force (pmf) is largely contributed from TELP near the membrane. As an extension to this theory, a novel phenomenon of mitochondrial thermotrophic function is now characterized by biophysical analyses of pmf in relation to the TELP concentrations at the liquid-membrane interface. This leads to the conclusion that the oxidative phosphorylation also utilizes environmental heat energy associated with the thermal kinetic energy (kBT) of TELP in mitochondria. The local pmf is now calculated to be in a range from 300 to 340 mV while the classic pmf (which underestimates the total pmf) is in a range from 60 to 210 mV in relation to a range of membrane potentials from 50 to 200 mV. Depending on TELP concentrations in mitochondria, this thermotrophic function raises pmf significantly by a factor of 2.6 to sixfold over the classic pmf. Therefore, mitochondria are capable of effectively utilizing the environmental heat energy with TELP for the synthesis of ATP, i.e., it can lock heat energy into the chemical form of energy for cellular functions.

Mitochondrial energetics with transmembrane electrostatically localized protons: do we have a thermotrophic feature?

James Weifu Lee
Transmembrane electrostatically localized protons (TELP) theory has been recently recognized as an important addition over the classic Mitchell's chemiosmosis; thus, the proton motive force (pmf) is largely contributed from TELP near the membrane. As an extension to this theory, a novel phenomenon of mitochondrial thermotrophic function is now characterized by biophysical analyses of pmf in relation to the TELP concentrations at the liquid-membrane interface. This leads to the conclusion that the oxidative phosphorylation also utilizes environmental heat energy associated with the thermal kinetic energy (k B T) of TELP in mitochondria. The local pmf is now calculated to be in a range from 300 to 340 mV while the classic pmf (which underestimates the total pmf) is in a range from 60 to 210 mV in relation to a range of membrane potentials from 50 to 200 mV. Depending on TELP concentrations in mitochondria, this thermotrophic function raises pmf significantly by a factor of 2.6 to sixfold over the classic pmf. Therefore, mitochondria are capable of effectively utilizing the environmental heat energy with TELP for the synthesis of ATP, i.e., it can lock heat energy into the chemical form of energy for cellular functions.
In the past, there was a common belief that the living organisms on Earth could only utilize light energy and/ or chemical energy, but not the environmental heat energy. Consequently, the life on Earth has been classified as two types based on their sources of energy: phototrophs and chemotrophs. Recently through bioenergetics elucidation studies with a new transmembrane electrostatic proton localization theory [1][2][3][4][5][6][7] now called also as the transmembrane electrostatically localized protons (TELP) theory, Lee discovered that certain biosystems such as alkalophilic bacteria Bacillus pseuodofirmus are capable of utilizing environmental heat energy isothermally with TELP to help drive ATP synthesis [8][9][10][11][12][13] . This discovery indicated that the protonic bioenergetic systems may have a thermotrophic function that is able to isothermally generate significant amounts of Gibbs free energy from environmental heat (dissipated-heat energy) [8][9][10][11][12][13][14] . It naturally raises a fundamentally important question: Do the mitochondria-powered chemotrophs have a thermotrophic function featured with isothermal environmental heat energy utilization as well? The answer to this scientific question is now positive [10][11][12][13]15 .
Recent studies 2,3,5,7,8,13,[19][20][21][22] showed that this Mitchellian textbook pmf equation (Eq. 1) could not really explain the bioenergetics in many biological systems. The most obvious evidences that clearly invalidate Mitchell www.nature.com/scientificreports/ ATP synthesis in biological systems [63][64][65][66] and thus proposed an alternative theory of energy transduction 67 . The mitochondrial membrane potentials in living cells are now known to be mostly about 56 mV, 105 ± 0.9 mV and 81 ± 0.7 mV 68 , 91 ± 11 mV and 81 ± 13 mV 69 , and also 114 mV 70 and 123 mV 71 . That is, Mitchell's chemiosmotic theory and its associated textbook pmf equation (Eq. 1) is not able to explain the protonic energetics even in mitochondria! This energetic problem of paramount importance will be addressed here in this article employing the TELP theory 1,2,[4][5][6][7]9,13 with new discovery in analyzing the protonic motive force (pmf) including the local pmf. Note, pmf can be translated to Gibbs free energy change G according to a simple relation ( �G = −Fpmf ) with the Faraday constant ( F ). Therefore, the study reported here may represent not only new progress in protonic bioenergetics, but also the latest discovery on a naturally occurring novel energy phenomenon: the thermotrophic function featured as isothermal environmental heat energy utilization with TELP to do useful work in driving ATP synthesis that we now know to occur in the mitochondria of our body as well.

Methods
Newly formulated protonic motive force equations with TELP. The bioenergetic process for oxidative phosphorylation (ATP synthesis) in mitochondria comprises the steps of creating an excess number of protons on the intermembrane space/cristae side of the mitochondrial inner membrane accompanied by a corresponding number of hydroxyl ions on the matrix side 5,13 typically, through the respiratory redox-driven electron-transport-coupled proton pumps across the mitochondrial inner membrane 7 . According to the TELP theory 1,2,[4][5][6][7]9,13,22 , the excess positively charged protons in the aqueous liquid on the intermembrane space/cristae side of the mitochondrial inner membrane will electrostatically become localized at the liquid-membrane interface, attracting an equal number of excess negatively charged hydroxyl anions to the other side (matrix) of the mitochondrial inner membrane to form a "protons-membrane-anions capacitor structure". This theory is built on the fundamental understanding that liquid water can act as a protonic conductor, which is well in line with the knowledge that protons quickly transfer among water molecules by the "hops and turns" mechanism first outlined by Grotthuss [72][73][74][75] . Therefore, a newly formulated equation for the protonic motive force (pmf) across a biomembrane considering TELP was introduced more recently through the author's latest publications 5,7,13,22 as Here �ψ is the "membrane potential from the p-side to the n-side" as defined by Mitchell 76,77 ,Nicholls and Ferguson 25,26 ; H + L is the TELP concentration at the liquid-membrane interface on the positive (p) side of the membrane; H + pB is the "proton concentration in the bulk aqueous p-phase" (intermembrane space in the case of mitochondria); and H + nB is the "proton concentration in the bulk liquid n-phase" (matrix in mitochondria) 5 . The first two terms of Eq. (2) comprise the "Mitchellian bulk phase-to-bulk phase proton electrochemical potential gradients" that we now call as the "classic" pmf, equivalent to that of Eq. (1); whereas the last term accounts for the "local" pmf from TELP at the liquid-membrane interface 7,13 .
For a protonic capacitor, as reported previously 2,5,7,13,22,51 , the ideal TELP concentration [H + L ] 0 on the positive (p) side of the membrane is a function of the transmembrane potential �ψ as expressed in the following equation: where C/S is "the specific membrane capacitance per unit surface area", l is "the thickness of the localized proton layer" 22 .
In mitochondria, non-proton cations in the bulk aqueous liquid phase can exchange with TELP at the liquidmembrane interface and thereby reduce TELP concentration. To account for the exchanging effect, this study employed the following equation developed recently also by Lee for the steady-state TELP concentration H + L in consideration of cation-proton exchange with each of the cation species M i+ pB of the bulk liquid p-phase at the equilibrium state as reported in the recent publications 2,5,7,13,22,51 , Here M i+ pB is the non-proton cation concentration in the bulk liquid p-phase, and K Pi is the equilibrium constant for the cation to exchange with TELP. The equilibrium constant K Pi is defined (and can be experimentally measured) as the ratio of the delocalized proton concentration H + pB to the cation concentration M i+ pB in the bulk aqueous p-phase when the cation-proton exchanging process reaches the midpoint at an equilibrium state where the steady-state TELP concentration H + L is equal to the localized cation concentration M i+ pL at the liquidmembrane interface 3 .
Note, each of the physical quantities appearing in Eqs. (2−4) may, in principle, be determined through experimental measurements. There are "no freely adjustable parameters". The calculations with Eq. (3) here used " C/S = 13.2 mf/m 2 as an averaged membrane capacitance based on measured experimental data" 78 and " l = 1 nm as a reasonable thickness of TELP layer" as previously explained in Ref. 2,5,7,13,19,22 . In the work reported in the next section, the TELP bioenergetics analysis with Eqs. (2−4) is extended to mitochondria not only by making a better treatment for cation-proton exchange 3,49 , but also by incorporating this treatment to calculate the total pmf (including both the "classic" and "local" pmf values) for the measured values for �ψ , H + pB , and H + nB using the well-documented animal mitochondria experimental data 61 .

Cation exchange reduction factors on TELP.
In the experimental study of "mitochondrial ATP-ADP exchange" by Chinopoulos et al. 61 , "mitochondrial membrane potential �ψ was measured in a range from 60 to 160 mV using fluorescence quenching with a cationic dye owing to its accumulation in energized mitochondria". The Chinopoulos et al. experimental results demonstrated the synthesis of ATP that was measured as the "ATP efflux rate" at a membrane potential �ψ as low as 60 and 80 mV. Furthermore, their experimental work 61 also demonstrated that "there is essentially no or little bulk-phase pH difference" between the matrix and the intermembrane space: the "∆pH max is only ~ 0.11". That is, under the given reaction medium pH 7.25 ( pH pB ), mitochondria matrix pH was about 7.35 ( pH nB ) during the state three. Another independent study 79 also consistently showed that the mitochondria matrix pH is about 7.3, which essentially is identical to that of the cytosol. These experimental observations are well corroborated with the Lee team's "experimental results from a biomimetic anode water-membrane-water cathode system where the bulk-phase liquid pH in the anode liquid chamber was observed to be essentially about the same as that in the cathode liquid chamber before and after energization by excess protons at one side of the membrane and excess hydroxyl anions at the other side" 3,8,19,49 . Therefore, the "measured experimental parameters (data) of the reaction medium pH 7.25 ( pH pB ) and mitochondria matrix pH 7.35 ( pH nB ) during the state three as reported by Chinopoulos et al. " 61 are employed in the biophysical analyses here as previously reported 5,7 . The concentrations of the bulk liquid-phase cations Na + , K + , and Mg 2+ in the mitochondria are listed in Table 1 (adapted from Ref. 5 ). As recently reported 5,7 , Table 1 presents the value for each of the cation-proton exchange constants K Pi for Na + , K + , and Mg 2+ used in the calculation. More importantly, it shows the "calculated cation exchange reduction factors of TELP concentration for pH pB = 7.25, which was the liquid medium pH used in the mitochondrial membrane potential ( �ψ ) measurement experiment" 61 . It also shows that "the total product of these factors (total cation exchange reduction factor in the denominator of Eq. (4)) is 1.29, which is close to one, indicating a relatively small role of cation-proton exchange in modulating TELP concentration at the liquid-membrane interface in mitochondria" 5,7 .

Results and discussions
Contribution from TELP to mitochondrial pmf. Transmembrane electrostatically localized protons (TELP) significantly contribute to the total protonic motive force (pmf) that can be translated to Gibbs free energy change ( �G = −Fpmf ) . Table 2 lists the values of mitochondrial pmf and the associated properties calculated as a function of transmembrane electrical potential difference �ψ using Eqs. (2−4) under the given reaction medium pH 7.25 ( pH pB ), mitochondria matrix pH 7.35 ( pH nB ) and taking cation-proton exchange into account as described above. The calculated pmf as a function of the transmembrane potential �ψ is displayed in Fig. 1, in comparison with the classic and local pmf contributions. These results showed that TELP dominantly contribute to the overall strength of pmf.
A pmf of at least 156 mV is required to overcome the mitochondrial phosphorylation potential of − 416 mV (− 9.6 × 4.187 ÷ F ) as calculated from the Gibbs free energy change (ΔG ATP ) of + 9.6 kcal/mol reported in Ref. 48,80,81 for ATP synthesis through oxidative phosphorylation with a proton-to-ATP ratio of 8/3 (416 mV/2.67 = 156 mV). The proton-to-ATP ratio of 8/3 is consistent with the animal mitochondrial F 0 F 1 -ATP synthase structure, which has 3 catalytic sites for ATP synthesis, driven by a flow of 8 protons per revolution through the 8 c-subunits in its rotary molecular machinery 60,[82][83][84] .
In consideration of an additional proton being used in the transport of Pi and ADP into, and of ATP out of, the mitochondrial matrix, the overall H + /ATP stoichiometry is increased by 1 to account for the additional proton (2.67 + 1 = 3.67). If the mitochondrial phosphorylation potential is conservatively assumed as − 416 mV Table 1. The total product (1.29) of cation exchange reduction factors calculated from experimental cation concentrations in the bulk liquid reaction medium for mitochondria (as reported in Ref. 5,7,61 ), their associated cation-proton exchange equilibrium constants, and cation exchange reduction factors for TELP concentration at the liquid-membrane interface with the bulk liquid reaction medium pH pB = 7.25. Adapted from Ref. 5 .  85 , the critical ATP hydrolysis free energy is approximately − 63.5 kJ mol −1 , which may be translated to a phosphorylation potential of − 658 mV (− 63.5 kJ mol −1 /F ). Accordingly, the minimally required pmf for ATP synthesis and transport together should be 179 mV (658 mV/3.67) in a functional animal cell. Table 2. Mitochondrial protonic motive force (pmf) and the associated properties including local pmf, calculated as a function of membrane potential �ψ using Eqs. (2−4) based on the measured properties ( pH pB , pH nB , �ψ) with the known reaction medium compositions of ref. 61 . The cation concentrations, protoncation exchange equilibrium constants and cation exchange reduction factor (1.29) are from Table 1; and the temperature T = 310 K. The "local" pmf is the last term in Eq. (2), while the first two terms of Eq. (2) give the "classic" Mitchellian pmf. Adapted and updated from Ref. 5 .   www.nature.com/scientificreports/ Note, the phosphorylation potential for ATP synthesis used by Slater 81 for his 1967 evaluation of the Mitchellian chemiosmotic hypothesis is + 15.6 kcal/mol that was measured by Cockrell et al. 80 in isolated rat liver mitochondria. Remarkably, this phosphorylation potential of + 15.6 kcal/mol (translating to 65.3 kJ mol −1 ) for ATP synthesis is quite close to the magnitude of the critical free energy − 63.5 kJ mol −1 for ATP hydrolysis in a functional animal heart cell obtained by Wu et al. 85 . According to this phosphorylation potential of + 15.6 kcal/ mol (translating to − 65.3 kJ mol −1 ÷ F = − 677 mV), the minimally required pmf for ATP synthesis in mitochondria should be 254 mV (677 mV/2.67) while the minimally required pmf for ATP synthesis and transport out of mitochondria would be 184 mV (677 mV/3.67).
To give the most benefits of doubts for the Mitchellian chemiosmotic theory and its associated textbook pmf equation (Eq. 1), we will use the minimally required pmf of 156 mV for ATP synthesis based on the conservative value of + 9.6 kcal/mol (translating to 40.2 kJ mol −1 ) 48 for the phosphorylation potential of ATP synthesis in the following analysis. As shown in Fig. 1, the total pmf including the contribution from TELP is well above the minimally required value of 156 mV to synthesize ATP, while the Mitchellian pmf is below this minimally required pmf at any of the �ψ values below 150 mV. This result indicates that the classic Mitchellian chemiosmotic theory could hardly explain the energetics even in mitochondria as noticed previously also by an independent study 48 . The mitochondria membrane potentials in a range from 40 to 140 mV have been determined in vivo in living human fibroblast cells with the techniques of fluorescence microscopic imaging and deconvolution 68 . It is now quite clear that the classic pmf alone cannot explain how the living cells are able to synthesize ATP and grow; In contrast, the synthesis of ATP and cell growth can now be well explained by the total pmf (Fig. 1) as calculated according to the TELP theory 1,2,4,5,7,9,13 .
Based on the in vivo measurement by Zhang et al. 68 using "fluorescent dye tetramethylrhodamine methyl ester" (which equilibrates between mitochondria and cytosol in living cells), the distribution of the mitochondrial membrane potentials in fibroblasts was determined to be mostly about 56 mV when the mean cytosolic dye fluorescent intensity was used as a threshold. With a somewhat higher threshold ("mean + SD") in their analysis 68 , the mitochondrial membrane potential ( �ψ ) values in fibroblasts and "N2a" cells were reported to be 105 ± 0.9 mV and 81 ± 0.7 mV, respectively. At any of these reported mitochondrial membrane potential ( �ψ ) values (56 mV, 105 ± 0.9 mV, and 81 ± 0.7 mV), the classic pmf value as presented in Fig. 1 is below the minimally required pmf value of 156 mV to drive ATP synthesis. That is, the "∆pH max is only ~ 0.11" observed by Chinopoulos et al. 61 can be translated to a delocalized pmf of no more than 6.8 mV; with these low mitochondrial membrane potential ( �ψ ) values, the total classic pmf values are no more than 63 mV, 112 ± 0.9 mV, and 88 ± 0.7 mV, which are all below the minimally required pmf of 156 mV to drive ATP synthesis. Therefore, the Mitchellian chemiosmotic theory and its pmf equation (Eq. 1) clearly failed to explain how the mitochondria can make ATP through oxidative phosphorylation to support the cell growth here.
An independent experimental study by Akhmedov et al. 86 showed that "the resting matrix pH was only 0.19 pH units higher than the cytosolic pH", which is remarkably similar to the observation of "∆pH max is only ~ 0.11" by Chinopoulos et al. 61 . If this ∆pH of 0.19 (translating to a delocalized pmf of 12 mV) is used to calculate pmf (with the mitochondrial membrane potential ( �ψ ) values of 56 mV, 105 ± 0.9 mV, and 81 ± 0.7 mV), the total classic pmf values would then be 68 mV, 117 ± 0.9 mV, and 93 ± 0.7 mV, all of which are also below the minimally required pmf value of 156 mV to drive ATP synthesis. Again, the Mitchellian chemiosmotic theory and its pmf equation (Eq. 1) clearly failed to explain how the mitochondria can make ATP through oxidative phosphorylation to support the cell growth here.
Akhmedov et al. 86 further showed that "during glucose stimulation (increasing the glucose concentration from 2.5 to 16.7 mM), mitochondrial matrix pH in the insulin-secreting INS-1E cells increased from pH 7.25 ± 0.04 to 7.78 ± 0.02" and "fifteen minutes after initiation of the glucose response, the ∆pH increased 3.5-fold to 0.66". If this high glucose-stimulated pH of 0.66 (translating to a delocalized pmf of 40 mV) is used with the mitochondrial membrane potential ( �ψ ) values (56 mV, 105 ± 0.9 mV, and 81 ± 0.7 mV), the total classic pmf values would then be 106 mV, 145 ± 0.9 mV, and 121 ± 0.7 mV; all of them are still below the minimally required pmf value of 156 mV to drive ATP synthesis. Therefore, once again, the classic Mitchellian chemiosmotic theory and its pmf equation (Eq. 1) again failed to explain how the mitochondria can make ATP to support the cell growth.
These findings are well corroborated with the mysterious problem previously noticed by Silverstein (2014) as a "thermodynamic efficiency of 113%" in mitochondria at a membrane potential of around 80 mV 48 . According to the classic pmf equation (Eq. 1), to avoid the "impossibly high efficiency (> 100%)" for mitochondria, one would have to "adjust" the bulk-phase "ΔpH (in-out)" to an arbitrary value of at least " + 2.5". However, it is now quite clear that the bulk-phase ΔpH (in-out) is nearly zero: "ΔpH max is only ~ 0.11" based on the modern experimental measurements 61 and modeling analysis of mitochondria 87 . The observed bulk-phase ΔpH of nearly zero in mitochondria 61 is also corroborated with the prediction from the transmembrane electrostatic proton localization theory 2,5,7-9 with the understating that mitochondrial inner membrane is rather impermeable to ions 88,89 . Another independent study using a pH-sensitive GFP 90 has recently also showed "that the intracristae lumen does not provide a reservoir for substrate protons for ATP synthesis" indicating "kinetic coupling of the respiratory chain with ATP synthase, but not proton gradients, drives ATP production in cristae membranes".
Note, the data in Table 2 and Fig. 1 show that the classic pmf value would be above the minimally required pmf of 156 mV only if the mitochondrial membrane potential is above 150 mV. In some of the recent literatures 91 , the bulk-liquid phase pH across the mitochondrial inner membrane was believed to be typically less than 0.5 unit, such as "the resting matrix pH (~ 7.6) and ΔpH m (~ 0.45) of HeLa cells at 37 °C" reported by Poburko et al. 92 . In that case, only when the mitochondrial membrane potential is above 130 mV, the classic pmf value would be above the minimally required pmf (156 mV). This is an interesting feature, since it can probably explain why it took so long in the last 50 years for the scientific community to fully realize the problems with the Mitchellian chemiosmotic theory including its pmf equation (Eq. 1). For example, like in the parable of "the blind men and an elephant" where a blind man touching the side of an elephant body may passionately feel the elephant "is a www.nature.com/scientificreports/ wall", some of the readers especially those who may have been fully occupied by the classic Mitchellian chemiosmotic theory looking at these classic pmf data points above the 156 mV level ( Table 2 and Fig. 1) such as a classic pmf value from 186 to 196 mV at a membrane potential ( �ψ ) from 180 to 190 mV previously reported in isolated mitochondria 93,94 could still feel hard to judge whether the Mitchellian chemiosmotic theory really has deficiencies or not. We now know that, even at the high membrane potentials ( �ψ ) such as 190 and 200 mV, its associated classic pmf value 196 and 206 mV as shown in Table 2 actually are still below the minimally required pmf of 254 mV (677 mV/2.67) for ATP synthesis in mitochondria, according to the phosphorylation potential of + 15.6 kcal/ mol (translating to − 677 mV) used by Slater 81 .
Anyhow, a good theory should be able to explain the bioenergetics with a full range of scientific data, for example, in the entire range of mitochondrial membrane potential ( �ψ ) values from 50 up to 200 mV, not just at the high membrane potential above 150 mV. Furthermore, the in vivo mitochondrial membrane potential ( �ψ) 70,95,96 values are mostly below 150 mV including 56 mV, 105 ± 0.9 mV, and 81 ± 0.7 mV reported by Zhang et al. 68 69 , and also 114 mV and 123 mV measured in swine and human respectfully using an improved PET-based method by Alpert et al. 70 and by Pelletier-Galarneau et al. 71 . We now understand that the classic Mitchellian chemiosmotic theory and its pmf equation (Eq. 1) cannot explain the mitochondrial energetics of oxidative phosphorylation in living cells because it fatally misses to account for the pmf contribution from TELP 5,7,13 .
As shown in Table 2 and Fig. 1, the local pmf from TELP was calculated to be a range from 306 to 343 mV, which is a function of the membrane potential ( �ψ ) in a range from 50 up to 200 mV. The total pmf is the sum of the classic pmf and local pmf, which in this case is in a range from 362 to 549 mV. The use of the local pmf (and the total pmf) can now easily explain how mitochondria are able to synthesize ATP since they are all well above the minimally required pmf of 156 mV at any of the mitochondrial membrane potential ( �ψ ) values in the entire range from 50 to 200 mV. Thus, the newly formulated pmf expression in Eqs. (2−4) consistently provides an excellent elucidation for the energetics of oxidative phosphorylation in mitochondria without requiring any arbitrary adjustment in the number of the bulk-phase "ΔpH (in-out)" that the previous study 48 with the classic Mitchellian doctrine had to require.
This new elucidation for the energetics in mitochondria can now also be well corroborated with the experimental results from the measurements of ATP production as a function of membrane potential ( �ψ ) during the state 3 after the addition of ADP in isolated rat mitochondria under oxidative phosphorylation experimental conditions 61 . Figure 2 presents the ATP production rate measured as the steady-state ATP efflux rate mediated through the adenine nucleotide translocase as a function of mitochondrial membrane potential ( �ψ ) manipulated through titration using uncoupler SF 6847 during the state 3 to various �ψ values in a range from 60 to 160 mV at various matrix pH (pHi) values in isolated rat mitochondria energized with 5 mM K-glutamate and 5 mM K-malate in the presence of dissolved O 2 , as reported previously in great details 61 . As shown in Fig. 2, the oxidative phosphorylation ATP production through the activity of mitochondrial ATP synthase is indeed still going on even in the low membrane potential range (60-150 mV) as predicted by the newly calculated pmf data ( Table 2 and Fig. 1) using the newly formulated pmf equations (Eqs. 2−4) where the total pmf values owning to the effect of TELP are still well above the minimally required pmf of 156 mV for ATP synthesis. Furthermore, the observed pattern of the ATP efflux rate which decreases as mitochondrial membrane potential ( �ψ ) is reduced (Fig. 2) is generally also in agreement with the pattern of the total pmf including the local pmf that is a function of the membrane potential ( �ψ ) as shown in Fig. 1.   Figure 2. ATP production rate in isolated mitochondria experimentally measured as the steady-state ATP efflux rate mediated through the adenine nucleotide translocase (ANT) with mitochondrial membrane potential ( �ψ ) manipulated to various values in a range from 60 to 160 mV at various matrix pH (pHi) values in isolated rat mitochondria energized with "5 mM K-glutamate and 5 mM K-malate" in the presence of dissolved O 2 . Adapted from Chinopoulos et al. 61 Fig. 2, based on the known metabolic pathways [97][98][99][100][101] , the conversion of a fed glutamate molecule to α-ketoglutarate and finally to oxaloacetate through part of the tricarboxylic acid (TCA) cycle in mitochondria could produce 1 FADH 2 plus 3 NADH (1 NADH from glutamate to α-ketoglutarate and 2 NADH from α-ketoglutarate to oxaloacetate) and an ATP by substrate-level phosphorylation (at the step from succinyl-CoA to succinate), while the use of a fed malate molecule to oxaloacetate with the TCA cycle can produce at least 1 NADH. The use of 4 NADH and 1 FADH 2 through the mitochondrial oxidative phosphorylation process could synthesize and transport 12.5 ATP [(4 NADH × 10 protons/NADH + 1 FADH × 6 protons/ FADH 2 ) / (3.67 protons/ATP)] out of mitochondria via the adenine nucleotide translocase. That is, the measured steady-state ATP efflux rate (Fig. 2) apparently is contributed mainly by the ATP production from the oxidative phosphorylation process. Merely a small fraction (1/13.5 × 100% = 7.4%) of the ATP production could be from the substrate-level phosphorylation of the TCA cycle, although the percentage from the substrate-level phosphorylation could be higher than 7.4% when the mitochondrial membrane potential ( �ψ ) is reduced to an extremely low level (near or below 60 mV) by titration using large amounts of uncoupler SF 6847.
Therefore, the newly calculated pmf data including the local pmf (Table 2 and Fig. 1) in conjunction with the experimentally measured ATP production (Fig. 2) have now clearly showed that ATP synthesis through oxidative phosphorylation by mitochondrial ATP synthase can indeed occur at a low membrane potential ( �ψ ) anywhere in a range from about 60 to 150 mV. This finding is remarkably in line with the independent observations of mitochondrial membrane potentials in living cells being mostly about 56 mV, 105 ± 0.9 mV and 81 ± 0.7 mV 68 , 91 ± 11 mV and 81 ± 13 mV 69 , and also 114 mV 70 and 123 mV 71 where apparently significant amounts of ATP are synthesized at such low mitochondrial membrane potentials to support the growth and activities of the living cells. That is, this new understanding may have major scientific implications with significant insights in explaining how ATP can be synthesized through mitochondrial ATP synthase with the effect of TELP at such remarkably low mitochondrial membrane potentials ( �ψ about 56 mV, 105 ± 0.9 mV and 81 ± 0.7; 91 ± 11 mV and 81 ± 13 mV, and also 114 mV and 123 mV) in the living cells 68-71 . TELP utilizing mitochondrial environmental heat energy for ATP synthesis: a thermotrophic feature? As presented in Table 2 and Fig. 1, all the total pmf values calculated through Eqs. (2−4) for mitochondria are well above the minimally required pmf value for ATP synthesis; Thus, the ATP synthesis in mitochondria can now be well explained for the entire range of membrane potential ( �ψ ) from 50 to 200 mV. However, it is surprising that many of these pmf values (such as 362, 369, 391, 404, 430 and 455 mV) are significantly larger than the "redox potential energy pmf upper limit (228 mV)" that could be maximally accounted for by "the redox-driven protonic pump system" based on the redox potential difference (about 1140 mV) "between the electron donor NADH (E m,7 = − 320 mV) to the terminal electron acceptor O 2 (E m,7 = + 820 mV)" 26 as analyzed previously in regarding to the redox potential chemical energy limit of 228 mV 13 .
Note, this redox potential chemical energy upper limit for the pmf of 228 mV is still smaller than the minimally required pmf of 254 mV (677 mV/2.67) for ATP synthesis in mitochondria, according to the mitochondrial phosphorylation potential of + 15.6 kcal/mol (translating to − 677 mV) 81 . Therefore, it also indicates that there must be another disparate energy mechanism (which we now know is the thermotrophic feature associated with TELP herein) to explain the mitochondrial energetics of oxidative phosphorylation.
How could the total transmembrane pmf exceed the mitochondrial chemical energy upper limit (228 mV)? We now understand that mitochondria can isothermally utilize the environmental heat energy (also known as the temperature-dependent molecular thermal motion kinetic energy) associated with TELP in "driving the synthesis of ATP from ADP and Pi through F 0 F 1 -ATP synthase". Therefore, mitochondria represent not only a chemotrophic system but also have a significant thermotrophic feature.
As recently reported 13 , the formation of a TELP layer apparently constitutes some kind of "negative entropy effect" 8,9,13 that "brings the excess protons to the mouths of the pmf users (F 0 F 1 -ATP synthase) where the protons can isothermally utilize their molecular thermal motions (protonic thermal kinetic energy k B T) possibly including their random and chaotic Brownian motions to push through the doors of F 0 F 1 -ATP synthase in driving ATP synthesis". That is, although the thermal energy-associated protonic translational motions (kinetic energy) are random and chaotic in all directions, "a localized proton at the water-membrane interface" has a much higher probability to chaotically hit through the mouth (F 0 protonic channel) of F 0 F 1 -ATP synthase in driving the "F 0 rotary molecular machinery for ATP synthesis" (thus the protonic thermal motion kinetic energy may be utilized) than "a delocalized proton in the bulk liquid phase that is far away from the protonic users". As shown with the third term in Eq. (2) "for the local pmf ", the thermotrophic function featured as "the utilization of protonic thermal kinetic energy k B T" is essentially expressed as "RT (= k B T ⋅ N A ) which equals to the product of the Boltzmann constant k B , the mitochondrial temperature T and the Avogadro constant N A ".
The "delocalized protons" in the mitochondrial cristae bulk liquid phase "are quite far away from the membrane surface"; Thus, the random and chaotic thermal motions of the delocalized protons in the bulk liquid phase are "not within the striking distance for them to hit into the F 0 F 1 -ATP synthase protonic channel to drive the rotary molecular machinery for ATP synthesis". As discussed also in the recent publication 13 , "the delocalized protons could also do the work when they are at the liquid-membrane interface near the F 0 F 1 -ATP synthase". Therefore, "the thermal energy factor RT (= k B T ⋅ N A )" is also in the second term of Eq. (2). However, the value of log 10 H + pB / H + nB was nearly zero in the case of mitochondria; consequently, the delocalized protons in mitochondria did not significantly contribute to the thermotrophic feature in driving ATP synthesis here.
As pointed out in the recent publication 13 , "the protonic bioenergetics systems operate widely in nearly all organisms known today". Therefore, this special thermotrophic function associated with TELP "has occurred probably for billions of years already on Earth". As presented in Fig. 1 www.nature.com/scientificreports/ function isothermally utilizing environmental heat energy associated with TELP. That is, the following local pmf equation has fundamental scientific significance in relation to the thermotrophic function featured as the "isothermal environmental heat energy utilization" with TELP.
According to Eq. (5), this special protonic thermotrophic function is mathematically related to the ratio ( H + L / H + pB ) of TELP concentration H + L at the liquid-membrane interface to the bulk-phase proton concentration H + pB in the intermembrane space/crista space at the same p-side of the mitochondrial inner membrane, which is quite surprising.
When the TELP-associated thermotrophic activities utilize mitochondrial environmental heat energy (k B T) in driving the molecular turbine of F 0 F 1 -ATP synthase for the synthesis of ATP from ADP and Pi , as discussed also in the recent publication 13 "a fraction of the environmental heat (k B T) energy may consequently be locked into the chemical form of energy in ATP molecules; and it would thus result in a small drop in the environmental temperature theoretically because of the TELP-associated isothermal environmental heat utilization". However, in mitochondria and the cells, "there are many other processes (including the glycolysis, tricarboxylic acid cycle, and the redox-driven proton-pumping electron transport activities as well as the ATP utilization processes such as ATP hydrolysis) releasing heat energy, which could mask the thermotrophic function that features as the isothermal environmental heat energy utilization process". Therefore, the energetic phenomenon in mitochondria (and the cells) may likely represent an interconnected mixture of both the chemotrophic processes and the thermotrophic processes. This subtle complexity could probably also explain why it took so long for the human beings on Earth to finally understand and figure it out here.
Based on the phosphorylation potential of + 15.6 kcal/mol (translating to − 65.3 kJ mol −1 ÷ F = − 677 mV) previously used by Slater 81 that is remarkably close to the magnitude of the critical free energy (− 63.5 kJ mol −1 ) for ATP hydrolysis in a functional animal heart cell obtained by Wu et al. 85 , the energy efficiency for the utilization of total pmf (including local pmf) in driving the synthesis of ATP can now be estimated. For example, according to the data in Table 2, at a mitochondrial membrane potential of 100 mV where the total pmf is 430 mV (including the local pmf of 324 mV), the energy efficiency for the utilization of total pmf in driving the synthesis of ATP is now estimated to be about 60% (677 × 100%/(2.67 × 430)), which thermodynamically appears to be a quite reasonable energy conversion efficiency. This also indicates that a substantial portion of the total pmf energy including the local pmf energy from TELP associated with the isothermal utilization of the mitochondrial environmental heat-related molecular motion kinetic energy (k B T) can indeed be locked to the ATP chemical energy.
Therefore, the discovery of thermotrophic feature isothermally utilizing environmental heat energy reported here "may have profound scientific and practical implications in bettering the fundamental understanding of bioenergetics and energy renewal" as recently pointed out for sustainable development on Earth 8 . With the new knowledge learned from this study, it may be possible to have the benefits of "mimicking this biophysical molecular scale process to create a new way in producing useful energy by isothermally utilizing the environmental heat energy from the ambient environment" 8,9,102 . It may also have scientific implications on the questions of who we are and how life began on Earth in terms of protonic bioenergetics. We now know that "water serves not only as a solvent and a substrate, but also as a protonic conductor" in living organisms including mitochondria 5,7 . The thermotrophic function isothermally utilizing environmental heat energy discovered herein was previously believed to be impossible for centuries. We now understand it is water as protonic conductor and consequently the formation of a protonic membrane capacitor with TELP that makes this thermotrophic function possible.
Conclusion with the discovery. The classic Mitchellian chemiosmotic theory and its textbook pmf equation (Eq. 1) cannot fully explain the energetics even in mitochondria. The TELP theory 1,2,4,5,7,9,13 with our updated pmf equations (Eqs. 2−5) can now fully elucidate the energetics of "oxidative phosphorylation in mitochondria" for the entire range of membrane potential ( �ψ ) from 50 to 200 mV ( Table 2 and Fig. 1). Through this study, we have now identified a novel thermotrophic function featured as "the isothermal utilization of environmental heat energy associated with the thermal motion kinetic energy of TELP in driving the synthesis of ATP" in mitochondria.
It is now, for the first time, numerically showed that the local pmf associated with TELP isothermally utilizing their environmental thermal motion kinetic energy represents about 75% of the total pmf which is the sum of the classic pmf and local pmf. The local pmf is now calculated, for the first time, to be in a range from 300 to 340 mV while the classic pmf is in a range from 60 to 210 mV in mitochondria with membrane potentials in a range from 50 to 200 mV, respectively. Apparently, a substantial fraction of the "environmental heat energy (k B T)" associated with TELP is "locked into the chemical form of energy in ATP molecules". Therefore, it is now quite clear that the low-grade environmental heat energy (i.e., the thermal motion kinetic energy of TELP) associated with the 37 °C human body temperature can be isothermally utilized by mitochondria "to perform useful work driving the synthesis of ATP". As shown in an example with a transmembrane potential of 100 mV, mitochondria obtain as much as 324 mV (translating to the Gibbs free energy change of − 31 kJ/mol) of local pmf through isothermally utilizing environmental heat (protonic thermal motion kinetic energy), which surprisingly represents 75% of the total pmf (430 mV, translating to the Gibbs free energy change of − 41 kJ/ mol); and only 25% is from the classic Mitchellian pmf component (106 mV, translating to the Gibbs free energy of − 10 kJ/mol). Consequently, we human as a mitochondria-powered organism is not only a chemotroph but also has a significant thermotrophic feature isothermally utilizing environmental heat energy in our human body environment to do work such as ATP synthesis.  www.nature.com/scientificreports/ Thermotropy of cells utilizing environmental heat energy was previously thought to be impossible for centuries. We now understand it is water as protonic conductor and thus the formation of protonic membrane capacitor that makes this thermotrophic feature possible. The energy efficiency for the utilization of total pmf in "driving the synthesis of ATP" is estimated to be about 60%, which again indicates a considerable fraction of the mitochondrial "environmental heat energy" associated with the thermal motion kinetics energy (k B T) of TELP being "locked into the chemical form of energy in ATP molecules". Therefore, this discovery may have profound scientific implications.
Quite clearly, one of the immediate scientific implications from this discovery is that we human as a mitochondria-powered organism is not only a chemotroph, but now understandably also carries a significant thermotrophic feature!

Data availability
All data generated or analyzed during this study are included in this published article and its Supplementary Information (xls efile of excel spreadsheet showing how the protonic motive forces were calculated) available online at the journal website.