A thermo-physical analysis of the proton pump vacuolar-ATPase: the constructal approach.

Pumping protons across a membrane was a critical step at the origin of life on earth, and it is still performed in all living organisms, including in human cells. Proton pumping is paramount to keep normal cells alive, e.g. for lysosomal digestion and for preparing peptides for immune recognition, but it goes awry in cancer cells. They acidify their microenvironment hence membrane voltage is lowered, which in turn induces cell proliferation, a hallmark of cancer. Proton pumping is achieved by means of rotary motors, namely vacuolar ATPases (V-ATPase), which are present at many of the multiple cellular interfaces. Therefore, we undertook an examination of the thermodynamic properties of V-ATPases. The principal result is that the V-ATPase-mediated control of the cell membrane potential and the related and consequent environmental pH can potentially represent a valuable support strategy for anticancer therapies. A constructal theory approach is used as a new viewpoint to study how V-ATPase can be modulated for therapeutic purposes. In particular, V-ATPase can be regulated by using external fields, such as electromagnetic fields, and a theoretical approach has been introduced to quantify the appropriate field strength and frequency for this new adjuvant therapeutic strategy.

Introduction

A fundamental characteristic of intracellular membrane compartments is the difference between their luminal pH and the bulk cytoplasm pH. The vacuolar ATPase is the main mechanism responsible for this pH differential1. These proteins include a class of proton pumps structurally homologous to the F-ATPases that produce ATP by using the proton-motive force across the mitochondrial inner membrane234.

A vast amount of new information has been obtained on the structure, mechanics, and biochemistry of the F-ATPases. For example, the rotation of the F1 motor was proven to develop over 40 pN nm and advance in three steps per revolution, with the hydrolysis of one ATP for each step5. Moreover, the F0 motor was found to consist of 10 or 12 subunits with rotational symmetry6. The F0 motor counters the large F1 torque by generating an even larger torque in the opposite direction to synthesize ATP. To do so it converts the transmembrane proton-motive force into rotary motion.

A mechanochemical model for the V-ATPase was suggested by Grabe, Wang, and Oster2. This model allows us to predict proton pumping rates over a wide range of environmental conditions which proves useful to determine the acidification of organelles. The model is based on the hypothesis that the V-ATPase works under normal operating conditions and that ATP concentrations are sufficiently high so that hydrolysis is not rate limiting2. The V-ATPase structure is composed of a counter-rotating stator and a rotor. The membrane-inserted/transmembrane section V0 is affected by the hydrolysis of ATP in the V1-soluble headpiece. A two-channel model and a one-channel model have been suggested to explain the rotor-stator assemblies7. They differ in the protons' path through the enzyme and communicate with the cytoplasm through the protons bound to the rotor section2. However, some experiments on sodium V-ATPases seem to support the one channel model8.

The accepted model for the active transmembrane ion transport is the alternating access mechanism. Ions are bound tightly on the low concentration side of the membrane. A conformational change weakens their binding affinity by exposing them to the high concentration side; as a consequence, they dissociate. Then, the pump changes its conformation in order to begin the cycle again2.

These types of processes are usually described by mechanochemical approaches. But, recently, an applied thermodynamic approach has been developed to analyse the biophysical and biochemical behaviour of the biological systems910111213. This approach introduces a new method in biophysics and biochemistry: the differences stem from the conceptual bases of the approach itself11.

Cells are systems, as of yet too complex to fully identify and dissect the contribution of each component and the interactions among them. Consequently, it is difficult to model an ideal sequence inside the cell and to develop a well posed thermodynamic approach. As such, the cell is considered as a black box, which maintains communication with its environment. The subject of the approach is precisely this communication: it is characterized by heat and mass transfer between a well-known environment and an open system, the cell. So, our attention is focused on the spontaneous flows. In particular, cells discharge wasted heat. This heat is the consequence of the internal cellular biochemical reactions, which are irreversible processes. So, from an applied thermodynamic point of view the heat flows are the wasted heat for irreversible processes of an open system in a non-equilibrium state91011. Consequently, the analysis of this irreversibility is fundamental to study the possible biological states of the cell: normal or diseased (e.g., cancerous). Indeed, normal and cancerous (or otherwise diseased) cells will dissipate different heat because different biochemical paths occur910. In applied thermodynamics, irreversibility is studied by introducing the concept of entropy generation and its link to the membrane flows12, as proven just by the constructal law1213141516.

This paper will develop a constructal analysis of the V-ATPase and will suggest some bio-medical consequences, ultimately aimed at improving the presently availably anticancer therapies.

To do so, in Section 2 the thermodynamic approach will be summarised, while we develop in Section 3 its applicability for V-ATPase analysis and, in Section 4, we propose some biomedical hypotheses.

The approach to irreversibility

Cells are living systems. They grow and, at a proper time and depending on the tissues they belong to, each of them divides into two different daughter cells. At that time, the size of any cell can vary as well as their daughters' sizes, a phenomenon occurring within a particular range. So a cell's life is a cyclic process. Indeed, it begins with the emergence from cell separation, and it ends with the separation of the daughter cells. Cells are composed of:

The cellular membrane, which controls the mass and energy inflow and outflow;

The cytoplasm, an aqueous solution containing thousands of structures and a vast array of chemical species;

The organelles, specialised subunits suspended in the cytoplasm, each enclosed within a membrane separating it from the cytoplasm. They perform specific, specialized, functions;

A network of tubular structures maintaining cell form, and allowing directional transport such as microtubules, cilia, organizer.

Within cells, chemical reactions occur that produce energy and macromolecules, and that increase and modify cell volume and its form17.

From a thermodynamic point of view, cells are open and complex systems able to convert their energy in the most efficient way for transport of substances across their membranes. They behave in two distinctively different ways: evolving towards maximum disorder or maintaining a high degree of organization in space and time. To do so, they must couple metabolic and chemical reactions with transport processes across their borders18.

Any system in nature has shape and structure13. They are macroscopic, finite size, and recognizable as patterns. The previous classical thermodynamic analysis highlighted that the flows in cell systems are fundamental to evaluate the behaviour of the systems themselves. Consequently, the analysis of the flows through the cell membrane appears fundamental in the comprehension of the biophysical and biochemical mechanisms inside the cell12. But this kind of analysis is powerfully described by the constructal theory. Indeed, by referring to the constructal law, a living system presents two characteristics: it flows and it morphs freely toward configurations that allow all its currents to flow more easily over time13. Life and evolution are a physics phenomenon, and they belong in physics16. Constructal law is a new approach introduced in thermodynamics in order to explain optimal shapes of natural structures13141516192021. The fundamental bases of the Constructal law was expressed21 as follows: “For a finite-size flow system to persist in time (to live), its configuration must evolve in such a way that provides greater and greater access to the currents that flow through it.”.

But, in a cell, a part of the energy is lost as heat outflow and only the resulting products of biochemical processes are known, while any individual step is inaccessible17. So, a constructal approach can represent a powerful theoretical method to analyse cell behaviour. Indeed, constructal theory highlights the fundamental role that flows across the system's border play in any thermodynamic process. This can represent a new viewpoint in the analysis of the biochemical and biophysical behaviour of cells. Instead of studying the cell, a very complex system, we can now study how the cells exchange components and information with their environments, and the interactions between cells and environments, which consist of the flows across the cell membranes. Indeed, cells are so complex that it is very difficult to understand the single effect of a given cellular process in relation to the ‘global’ result for the cell. Consequently, the study of cells can be developed by introducing the black box model, and considering only the spontaneous cell flows. Therefore the spontaneous heat cell exchange represents the interaction or, here, spontaneous communication between the cell and its environment. Lastly, it is, of course, easier to access the environment than the living cell.

Therefore, we decided to analyse the heat and mass flows across the membrane. This is what following a constructal approach suggests. But, in the analysis of the cell membrane, in relation to the molecular motors involved in this process, we have no useful data to evaluate directly the flows. Consequently we use constructal theory as a new fundamental viewpoint but we need a related method to develop the calculations. Notably, the heat flow is the consequence of the irreversible processes within cells and this is easily developed by using the Gouy-Stodola theorem22; it considers only the work lost for irreversibility and the temperature of the system's environment. This constructal based approach is theoretically interesting because22:

It allows to obtain the physical conditions in which an open system persists in its stationary states;

This is a power description of complex phenomena because it allows to evaluate the global effects and their fluctuations around the stationary states;

It involves the definition of exergy flows, which are the flows of the available energy of the system.

In 1873, Gibbs introduced the available energy, today named exergy. It is the function that expresses the maximum useful work that a system can obtain in a thermodynamic equilibrium with its environment. The exergy lost or dissipated E, i.e. the available energy or work lost W, in an irreversible process, for us the heat emerging from the cell, can be obtained through the Gouy-Stodola theorem21: with T0 the environment temperature and S the entropy generation. By using this relation and evaluating any process across the membrane, always by using the well accepted applied thermodynamic relations21, the exergy flows across the membrane can be related to the entropy generation for a cell which has been recently obtained9101112 as: where:

S is the entropy generation due to the thermal flux driven by the temperature difference, in which τ1 is the lifetime of this process;

T is the temperature;

S is the entropy generation due to the diffusion current driven by the chemical potential gradients, in which τ2 is the lifetime of this process, μ is the chemical potential of the i-th species;

S is the entropy generation due to the velocity gradient coupled with the viscous stress, in which τ3 is the lifetime of this process;

S is the entropy generation due to the chemical reaction rate driven by the affinity, in which τ4 is the lifetime of this process, N is the number per unit time and volume of the i-th chemical reaction and is the affinity, evaluated as the variation of the standard Gibbs' free energy;

S is the entropy generation due to the dissipation that results from the interaction between external forces and the system, in which τ5 is the lifetime of this process, F is the force generated by the interaction with the external field and J is the associated flow.

Where the volume of the cell is defined as12 with δ1(t) and δ2(t) the long and the short axes dimensions of the cell, but δ1(t) and δ2(t) must be experimentally evaluated, so for a theoretical approach the diameter of the cell is approximated as the diameter of a sphere12: Consequently, the characteristic length of the cell results in:

r = L/2 is the cell radius

The mean environmental temperature can be assumed to be9101117 T0 = 310 K and the mean cell temperature has been estimated to be T0 + ΔT;

ΔT is the difference between the temperature inside the cell and that of its environment910111217. It has been evaluated as ΔT ≈ 0.4°C, but it would be different for each cell line and for each cell line it would have to be different between normal and cancerous (or more generally, diseased) states9101112;

The characteristic length171920 can be evaluated as L = 2 r;

The internal energy density u can be evaluated as the ratio between the cell's mean internal energy91011121718, considered the same as that of ATP, U = 3 × 10−7 J and the mean value of the cell inside the human body V = 7600 μm3, hence the cell volume in the human body being in the range 200–15000 μm3, so it results in u = 3.95 × 107 Jm−3;

The thermal molecular mean velocity inside the cytoplasm is considered to be917 ;

d can be assessed as17 d = 0.2 r;

The membrane volume is evaluated as9101112 with π = 3.14;

The chemical potential gradient91011121718 can be calculated as the ratio between the mean value of the chemical potential μ = 1.20 × 10−9 J kg−1 and the membrane length d = 0.01 μm, with the mean density being ρ = 1000 kg m−3;

The viscosity910111217 is taken to be 6.91 × 10−3 N s m−2;

η~ 2.07 × 10−3 N s m−2 at910111217 30°C;

is evaluated as910111217 3.0 × 10−6 m s−1;

F is the external field and J is the associated flow.

Constructal approach to V-ATPase

The basis of metabolism energetics consists in the generation and the hydrolysis of ATP, which occurs across a trans-membrane electromotive gradient. Indeed, this conversion of energy can be obtained by transitioning the electrochemical energy into the chemical energy of the terminal phosphoric anhydride bond of the ATP26. This can occur by the action of an enzyme, which works as a proton-pumping ATP synthetase. Here the proton pump V-ATPase will be discussed by introducing just the constructal approach.

As described in the introduction, the V-APTase works through a counter-rotating stator and a rotor mechanism. It hydrolyses ATP to obtain the required energy for its work. The fundamental reaction is: and, consequently, a H+ ion is pumped into the cell: where out means outside, in refers to inside and memb stands for across the membrane.

This proton-pumping can be modelled considering the membrane as an electric RC-circuit, while the V-ATPase can be modelled as a DC motor as represented in Figure 1. Indeed, the V-ATPase rotor can be considered to be the equivalent of a simple DC-motor rotor. The energy required by the rotor movement is generated by the energy conversion of the ATP-hydrolysis (7), while the stators rotate as gears dragged by the rotor itself. Moreover, all the DC motors convert electric energy into work with high efficiency (about 1), so, introducing this model for the rotor, the irreversibility results only in the gears. Consequently, the efficiency of the V-ATPase can be evaluated as2: where ΔG is the free energy variation due to the hydrolysis of a single ATP molecule (~21 k = 50 kJ mol−1, being k the Boltzmann constant and T the temperature), χ is the coupling ratio (χ = J/J, being J the proton flux and J the ATP hydrolysis rate) and ΔG is the free energy variation required to move the proton across the membrane2: with Δϕ being the membrane potential, R is the gas constant (8.314 J mol−1K−1), F is the Faraday constant (96.485 × 103 A s mol−1), and 2.3 ΔpΗ is the physiological concentration gradient. The coupling ratio χ is affected both by the pH gradients and by the membrane potential. The average rotation ν can be calculated as a function of the load τ as7: Under physiological conditions, this leads23 to 15–20 Hz.

Figure 1

The equivalent schema for the V-ATPase system.

It is a DC motor in series with a rectifier. The V-ATPase rotor can be considered to be the equivalent of a simple DC-motor rotor. The energy required by the rotor movement is generated by the energy conversion of the ATP-hydrolysis, while the stators rotate as gears dragged by the rotor itself. The DC motor converts electric energy into work with high efficiency (about 1). Consequently, the irreversibility results only in the gears.

From these results the work dissipated in wasted heat yields: and the related power lost in heat flow: As a consequence of the previous Section it is now possible to evaluate the entropy generation due to membrane flows as: and the related entropy generation rate as: with S being zero because the V-ATPase is not driven by temperature difference, and S being zero because the V-ATPase is not driven by the velocity gradient.

Considering relations (14) and (15) together with relations (12) and (13) it is now possible to argue that the entropy generation, and consequently the irreversibility, depend on:

The chemical potential at the membrane,

The affinity,

The electric potential at the membrane,

The H+/ATP rate,

The pH gradient, and

The working temperature.

But, as highlighted in the introduction, these quantities are characteristic quantity of the thermodynamic state of a cell, which is different between normal and cancerous cells. Consequently, these quantities are also different between normal and cancerous cells of the same cell line. So, for a cancer cell it follows: where the significant quantities are considered different from normal cells; to indicate it, a c symbol has been introduced. As a consequence, we get a different value of the entropy generation, and the variation between a cancer cell and a normal cell leads to:

Furthermore, if we are able to change the entropic behaviour of the tumor cell it is possible to compel the cancer to behave as it would be a normal cell. To do so, the component S related to the dissipation due to work by interaction with the external field can be introduced obtaining the entropy generation (always multiplied with the environmental temperature to obtain energy balances as previously done) as: and the related entropy generation rate: To obtain the required effect it must be:

The effect could be obtained by using2425262728293031:

catalysis32,

electric field interaction33,

electromagnetic or ultrasound waves,

molecular machines,

local inflow of nano-particles of ferro-fluids in interaction with a magnetic field, but this technique has yet to be designed,

and/or coupling some of these possible techniques.

Here, we want to highlight that this conceptual therapy is meant to become a supporting strategy to the currently applied anticancer therapies.

Considering the use of magnetic fields it is possible to argue that any such field would modify the rotation of the diseased cell towards the normal one. So, it must induce an electric field strong enough to obtain the normal rotation of the rotor. That means, the frequency of the electric field must be such that: and it would supposedly to be in the 0 ÷ 40 Hz range, which is twice the physiological one and this will not damage the normal cell(s). The torque can be related to the membrane potential and results in values of the order of2 10−21 Nm. Now, this leads us to consider that this torque can be obtained also as34: and that the relation between the electric field and a magnetic field for an electromagnetic wave is35: with μ being the magnetic permeability and ε representing the electric permittivity of the cell membrane.

Considerations

All types of ATPases, i.e., the A-ATPase of Archaea, the E-,F-P- and the V-ATPases are essential for life and all of them generate an electrochemical ion gradient across the membrane and hydrolyse or synthesize ATP. Also, from a structural point of view, they are similar; indeed, they are enzymatic complexes, which work as molecular rotary motors. In particular, V-ATPase plays an important role36 in receptor-mediated endocytosis1, in intracellular transport, and in the acidification of late endosomes37.

Metabolism is a cycle composed of a continuous sequence of oxidations and reductions. Considering our DC motor approach, the electromotive force is the free energy variation, ΔG, required to move the proton across the membrane represented by the Nernst equation (10). From this equation it possible to highlight that any variation of the electromotive force of the V-ATPase rotor can determine two possible consequences:

A change in the membrane potential, or

A change in the pH.

Of course, these effects can also occur together.

V-ATPase supplies energy to the membrane, regulates intracellular pH, and causes extracellular acidification or alkalinization. Its activity is fundamental for both vacuolar acidification in response to glucose metabolism and the regulation of pH homeostasis36. In particular, vacuolar acidification was found in the transport of lysosomal enzymes from the Golgi apparatus to the lysosomes38 and V-ATPase contributes to the homeostasis of cytoplasmic pH. Consequently, V-ATPase has a fundamental role in securing the microenvironment required for correct protein transport and membrane exchanges. As such, any breakdown in the V-ATPase can generate a different acid microenvironment. Indeed, V-ATPase breakdowns can cause lower pH in the microenvironment and this can be exploited by lysosomal enzymes of cancer cells that participate in the degradation of the extracellular matrix required for invasion and metastasis3940.

Four regulatory mechanisms are known36:

The regulation of pump density, useful to maintain their cytoplasmic and vacuolar pH stable;

The regulation of V1 and V0 domain association/dissociation;

The regulation of secretory activity, thereby maintaining the balance in the formation of bisulfite and binding efficiency between H+ and the pump;

The modifications of the membrane potential due to electrogenic force.

Using a thermodynamic approach, the aim of this paper is to highlight how external electromagnetic waves can modify the mechanical behaviour of the V-ATPase rotor. So, now, we must evaluate the intensity and the frequency of the applied external magnetic field to obtain some beneficial anticancer effects yet without endangering the patient. To do so, we consider that any breakdown in the rotor can be modelled as a difference Δτ in the torque in comparison with the normal value τ, such that: By using this last relation in equation (22) it follows that with E + ΔE denoting the electric field related to the breakdown of the V-ATPase rotor. Consequently, the effect of an external magnetic field must be to contrast this field variation, so that, considering relation (26), it follows: Now, considering a torque variation for the breakdown in the order of 0 ÷ 3 pN nm−1, so that the magnetic permittivity is about the value of the permittivity in the air, μ ~ 10−6 H m−1, and so that the electric permeability is around ε ~ 10−6 ÷ 10−5 F m−1, ~ 10−12 ÷ 10−11 s m−1 and that the electric field at the membrane is around 107 V m−1 it follows that the therapeutic magnetic field must be of the order of 10−5 ÷ 10−4 T, i.e., of the same order as the Earth's magnetic field.

Now, from relation (24) we can obtain the possible range of frequencies of the magnetic wave. Considering that the normal value of the torque is around 10 pN nm−1, we assume a torque variation for the breakdown of the order of 0 ÷ 3 pN nm−1; so, it is possible to obtain a frequency range of the order 0 ÷ 10 kHz for therapeutic use.

Moreover, this frequency can be confirmed by evaluating the time τ2 related to the diffusion current driven by chemical potential gradients. It can be calculated by considering the macroscopic phenomenon of diffusion across the membrane as: with d ≈ 10−8 m being the length of the membrane and D the diffusion coefficient. Considering that the diffusion coefficient of the ions across a membrane is around42 10−12 ÷ 10−8 m2s−1 it follows τ2 ≈ 10−4 ÷ 1 s, which is equivalent to a frequency ν = 1/τ2 of 0 ÷ 10 kHz.

Conclusions

Proteins play a fundamental role in ion transport because membranes represent a barrier to free diffusion of molecules.

In active transport an ion crosses the membrane against its electrochemical potential. The required energy for this process is obtained by the chemical energy released from hydrolysis of ATP, or pyrophosphate, or from the movement of a co-transported or coupled ion along its electrochemical gradient41. Indeed, the coupled downhill and uphill movement of ions is a common transport path through the membrane. In this context, the role of the pumps is fundamental, with particular regards to V-ATPase, i.e. the pump managing the H+ transport. Indeed, it moves positive charges into and from the cytoplasm through the hydrolysis of ATP, with the consequence of establishing a large membrane voltage (inside negative and outside positive) and, consequently, a pH gradient41 of about 400 mV for H+. The V-ATPase is composed of a membrane extrinsic and a membrane intrinsic sector, and couples catalysis of ATP hydrolysis to proton transport by a rotational mechanism42. V-ATPase is fundamental in the analysis of cell behavior because the H+ gradient established by this molecular pump is used to drive coupled active movements of other ions across the cell membrane. One example is the Cl− transport: indeed, Cl− is actively transported across the membrane because the membrane potential is more negative than the equilibrium potential for this ion.

The ion channels and transporters provide different permeability to distinct ions, such as Na+, K+, Ca2+, and Cl−. As a consequence of the asymmetry in these ion distributions, a membrane potential exists between the cytoplasm and the extracellular environment. It is expressed relative to the extracellular environment and a cell depolarizes if the membrane potential is relatively less negative, and vice versa43. The membrane potential can be calculated by using the Goldman–Hodgkin–Katz equation4445: where P is the permeability of the ion, [A] means concentration of the A-ion, R is the ideal gas constant (8.314 J mol−1K−1), T is the temperature, and F is the Faraday constant (96.485 × 103 A s mol−1). From relation (31) it is possible to state that the membrane potential can be changed by alterations in the conductance of one or more ions. In particular, from the previous considerations on the V-ATPase we can highlight that this pump can change the transport of H+ and, as a consequence of the Cl−-H+ coupled transport, it can change the Cl− transport. Consequently, both the membrane potential and the pH are changed by any alteration of the V-ATPase.

It is noteworthy that in the analysis of the mitotic activities in sarcoma cells, the membrane potential was found to undergo hyper-polarization before entering M phase. It suggests that the level of Δϕ is correlated with cell cycle progression. Moreover, membrane hyper-polarization was shown to block reversibly DNA synthesis and mitosis and to be correlated with the level of differentiation464748. Consequently, the membrane potential represents a fundamental quantity for the control of critical cell functions, particularly, with regards to proliferation, migration, and differentiation. To support this conclusion, some experimental evidence can be cited: First, direct in vitro and in vivo comparisons of the membrane potentials have highlighted that cancer cells are more depolarized in relation to normal cells; some evidence of this electrochemical behavior of the cancer cells can be summarized as follows43: between normal and cancerous breast cells49, hepatocytes and hepatocellular carcinoma cells50, normal and neoplastic adrenocortical tissues51, normal embryonic fibroblasts and fibrosarcoma52, benign and cancerous skin cells53, and between normal and cancerous ovarian tissue54.

Lastly, cell migration is controlled by the movement of ions and water43 in that an acidic environment furthers this phenomenon. This environmental pH is regulated by the H+ concentration, which is related to the V-ATPase functions. In addition, the membrane potential is considered an indirect factor of cell migration, strictly related to the electrical driving force for Ca2+ whereas a hyperpolarized membrane potential increases intracellular Ca2+ through the TRP channels; in contrast, membrane depolarization activates the Ca2+ channels55. Notably, migrating cells have a high intracellular Ca2+ concentration gradient56.

In summary, all these findings highlight that V-ATPase-mediated control of the cell membrane potential and that of the environmental pH can potentially represent a valuable support strategy for anticancer therapies. In here, a constructal approach has been used to study how V-ATPase can be modulated for therapeutic purposes. In particular, V-ATPase can be regulated by using external fields, such as electromagnetic fields, and we have introduced a theoretical approach to quantify the appropriate field strength and frequency for this new adjuvant therapeutic strategy. Here, in contrast with the usual applications of constructal theory575859, the flows have been assessed by evaluating the entropy generation because its variation is strictly related to the flows60 across the cell membrane.

Author Contributions

U.L. developed the thermodynamic approach. U.L., A.P. and T.S.D. contributed to the main manuscript text. All authors (U.L., A.P., T.S.D.) reviewed the manuscript.