 Research
 Open Access
 Published:
Bioengineering thermodynamics of biological cells
Theoretical Biology and Medical Modelling volume 12, Article number: 29 (2015)
Abstract
Background
Cells are open complex thermodynamic systems. They can be also regarded as complex engines that execute a series of chemical reactions. Energy transformations, thermoelectrochemical processes and transports phenomena can occur across the cells membranes. Moreover, cells can also actively modify their behaviours in relation to changes in their environment.
Methods
Different thermoelectrobiochemical behaviours occur between health and disease states. But, all the living systems waste heat, which is no more than the result of their internal irreversibility. This heat is dissipated into the environment. But, this wasted heat represent also a sort of information, which outflows from the cell toward its environment, completely accessible to any observer.
Results
The analysis of irreversibility related to this wasted heat can represent a new approach to study the behaviour of the cells themselves and to control their behaviours. So, this approach allows us to consider the living systems as black boxes and analyze only the inflows and outflows and their changes in relation to the modification of the environment. Therefore, information on the systems can be obtained by analyzing the changes in the cell heat wasted in relation to external perturbations.
Conclusions
The bioengineering thermodynamics bases are summarized and used to analyse possible controls of the calls behaviours based on the control of the ions fluxes across the cells membranes.
Background
Nature, from a physical, biological, chemical and mathematical point of view, is a complex system, while from an engineering point of view, it is the “first” engineer! In particular, cells can be modelled as an adaptive thermal and chemical engines which convert energy in one form to another by coupling metabolic and chemical reactions with transport processes [1–5], by consuming irreversibly [6–8] free energy for thermal and chemical processes, transport of matter, energy and ions.
Energy is a thermodynamic property of any system in relation to a reference state, which changes during any process, while its total amount remains constant in relation to the universe, considering it as the system together with its environment. In cells, many processes such as replication, transcription and translation need to convert molecular binding energy, chemical bond hydrolysis and electromagnetic gradients into mechanical work, related to conformational changes and displacements [9]. The biomechanical analysis of DNA has pointed out the connections among forces, thermodynamics, nanomechanical and electromagnetic behaviour of biological structures and kinetics [10].
Engineering thermodynamics is the science which studies both energy and its best use in relation to the available energy resources with particular regards to energy conversion, including power production, refrigeration and relationships among the properties of matter, including also living matter. So, engineering thermodynamics can be introduced in the mechanobiological and system biological approach in order to improve these sciences by analysing the biosystems also from a thermal point of view: a new engineering science could be considered, the bioengineering thermodynamics. Indeed, the first law of thermodynamics expresses the conservation of energy, while the second law states that entropy continuously increases for the system and its environment and introduces the statistical and informational meaning of global quantities [11–14].
In this paper we develop the bioengineering thermodynamic of biological cells, with particular regards to possible control of the cells growth by a control of the ions transport across the cell membrane. To do so, we consider that cells spontaneously exchange heat, and this heat is related to their biochemical and biophysical behaviour. This wasted heat represents the interaction between the cell and its environment, a sort of “spontaneous communication” towards environment. This interaction is fundamental to developing a thermodynamic study of the cell. Indeed, cells are too complex to understand the contribution of each process to the global result, and the study of cells as black boxes allows us to simplify the analysis by considering only the inflow and outflow balances [15]. Moreover, it is easier to have access to the cell environment than to the living cell itself. These considerations allow us to introduce the bases of the bioengineering thermodynamic approach introduced in the study of the cells:

1.
An open irreversible real linear or nonlinear system is considered;

2.
Each process has a finite lifetime τ;

3.
What happens in each instant in the range [0,τ] cannot be known, but what has happened after time τ (the result of the process) is wellknown (at least it is sufficient to wait and observe): local equilibrium is not necessarily required;

4.
The balance equations are balance of fluxes of energy, mass and ions.
The fundamental quantity used in this analysis is the global entropy [16, 17], related to systems changes, highlighted as the only effective criterion for spontaneity of change in any system, with particular regards to the entropy variation due to irreversibility, named entropy generation [18], which is the result of the global effect of the entropy variation

1.
due to the interaction with the environment

2.
within the system itself.
The introduction of entropy generation comes from the need to avoid inequalities: entropy results as a state function, so nothing is really produced or generated. Therefore, entropy is nothing more than a parameter characterising the thermodynamic state, and the term due to irreversibility, S _{ g }, measures how far the system is from the state that will be attained in a reversible way [12]. It is always S _{ g } ≥ 0.
Recently, it has been highlighted that any effect in Nature is always the consequence of the dynamic balances of the interactions between the real systems and their environments [12] and the real systems evolution is always related to the decrease of their free energy, in the least time [19–21]. So, bioengineering thermodynamics is based just on two fundamental concepts of physics: interactions and flows. The result is the analytical formulation of flowbased analysis in thermodynamics, which can play the role of a “rallying point” of the different modelling approach to biosystems. Indeed, if we consider natural systems we can highlight that they are always open systems, which means that they can exchange heat and mass with their environment. So, the interaction with the environment is a fundamental concept for the thermodynamic analysis.
We consider the environment as a thermostat and the system, together with its environment, is an adiabatic closed system [18]. But, for an adiabatic close system, the total entropy, defined as:
it always increases, as a consequence of the second law [18]. In relation (1) dS is the variation of the total entropy elementary, d _{ e } S is the entropy variation for interaction between the open system considered and its environment, and d _{ i } S is the entropy variation due to irreversibility, such that:
Now, we can write the relation (1) as [22]:
where Q is the heat flow, T is the temperature, V is the volume, t is the time and ṡ _{ g } is the density of the entropy generation rate. Now, we consider that the stationary states of the open system correspond to the equilibrium states of the adiabatic closed system. Considering the system together with its environment, we are analyzing an adiabatic closed system, so the entropy variation for the volume considered is maximum at the equilibrium [23]:
and
This last relation allows us to state that the flows between the open system and its environment cause the entropy generation rate density, so the interaction between system and environment is responsible of irreversibility. But, we cannot state if the cause of changes is the change of the entropy inside the cell or the fluxes across the cell membrane. We can only highlight the relation between changes and fluxes, but this approach doesn’t allow us to establish if are the fluxes to cause entropy changes or if entropy changes causes fluxes.
Now, considering that the entropy generation rate density can be written as [22]:
where J _{ k } is the flow of the kth quantity involved in the process considered and X _{ k } is the related thermodynamic force. Now, considering that:
the relation (5) becomes:
in agreement with Le Chatelier’s principle [24], for which any change in concentration, temperature, volume, or pressure generates a readjustment of the system in opposition to the effects of the applied changes in order to establish a new equilibrium, or stationary state. It follows that the fundamental imperative of Nature is to consume free energy in least time. Any readjustment of the state of the system can be obtained only by generating fluxes of free energy which entail any process where the system evolves from one state to another.
Results and discussion
The existence of bioelectric signalling among most cell types suggests a wide field of applicability of these electromagnetical signals. Here, we provide bioengineering thermodynamic theory that suggest how to explain the effects of energy, mass and ionic flows across cell membranes and, consequently, to control the cell behaviour by a control of ion fluxes.
Living cells are separated from their environment by the lipid bilayer membrane, which presents a different concentration of specific ion species on both sides. As a consequence, a charge separation across the membrane is generated by the electrodiffusion of ions down their electrochemical gradient. These ions move into a negative (inside the cell) membrane potential of around −70 to −100 mV. The hydrophobic component of the lipid bilayers behaves as a capacitor dielectric, which maintains the ionic gradients across the membrane; in some instances, the action of ATPdriven ionic pumps supports this effect by separating the charges. The cell function is regulated by the membrane proteins, sensitive to electric field; indeed, changes in the electric field are transduced into a conformational change that accomplishes the function of the membrane protein with consequences for the regulation of cell functions. The charged species, their arrangements, the local field strength, charges and dipoles disposition and movements can vary with the result of changing the electric field which is tranduced into a conformational change related to the protein functions themselves [32].
These considerations suggest that control and regulation of the membrane’s electric field could represent a new approach to therapies against diseases such as cancer. To understand how to control the fluxes across the membrane we consider the concentration of the ions on the opposite sides of the membrane [33]:
where c is the molar concentration of the chemical species, R is the universal constant of gas, T is the temperature and Φ is the electric potential energy. As a consequence of this concentration difference the cell can move the ions, and change the pH inside and outside its membrane. The ion drift velocity v _{ drift } across the cell membrane can be obtained by using the classical kinetic theory [34] as:
where Ze is the electric charge of the ion, m is the ion mass, ϕ is the electric potential across the membrane, d is the length of the membrane and τ _{ drift } is the mean time between two collisions [33]:
where σ is the electric conductivity. Consequently, an electric current I occurs for each ion i = H^{+}, Na^{+}, K^{+}, Ca^{2+}, Cl^{−}, Mg^{2+}, etc.:
where A is the mean surface area of the membrane. Now, considering the equivalent RC electric circuit for a membrane it is possible to state that the resonant frequency for such a circuit results in (2πRC)^{−1}, where R is the electric resistivity for the ion considered and C is the membrane capacity.
It follows, that if we want to control the crossmembrane flux we must impact the current. The easier physical way to interact with a current is to use an electromagnetic wave of the resonant frequency for the membrane, in relation to the ion considered, with its amplitude being related to the entropy generation as just obtained in Ref. [25–30].
In Figs. 1 and 2, it is represented an example of this kind of control. Figure 1 represents the natural behaviour of cell requirement of energy to grow. Figure 2 represents the cell requirement of energy by cell to grow when they are inside an electromagnetic field. It represents the ratio between the variation in percentage of the energy used by a cancer in a magnetic field (50 μT, 40 Hz) respect the energy used by a cancer outside of the field, related to the energy used by the cancer outside the field vs the growth of the cancer in terms of volume growth (ratio between the cancer volume during the cancer growth and the initial volume). It has been obtained by using the entropy generation approach described in the following section on methods. It is possible to highlight how the different ions have different effect. The positive ions determines a decreasing of the energy used while the negative ion increase it. So the positive ions determine an opposition to the growth. The more effective ion is Ca^{2+}. It means that a control of calcium ion can determine a control of the volume growth of a cancer. Here, the control is suggested by the use of an electromagnetic field. The field induces in the cell a greater use of energy to obtain the same growth.
Conclusions
Life is an organisational and thermodynamic process that tends towards the maximum conversion of available energy. The biochemical reactions produce or consume external metabolites, and they connect internal metabolites, in constant concentrations in the cells at their steady states. To do so, the cell must exchange energy and matter through its membrane. The fundamental phenomena used by cells to reach their optimality consist of a redistributing of the flow patterns through their metabolic network.
By using the bioengineering thermodynamics, it has been highlighted how the different ions have different effect on the use of energy by the cell to grow. To do so, a control of the cells behaviours is introduced. Here, an electromagnetic field is used as a control system, but other field could be used. Cells inside and outside an electromagnetic field have been considered. The positive ions determines a decreasing of the energy used by the cancer, such that the cancer cannot grow as outside the field. On the other hand, the negative ion increase the use of energy. It means that a control of ions can determine a control of the volume growth of a cancer. This result can be extended to all the molecular fluxed across the cell membrane, obtaining a possible bioengineering thermodynamic approach to control the cancer growth.
Methods
The approach previously used is based on the following considerations:

1.
The energy lost by a system is gained by the environment, consequently, the information lost by the system is gained by the environment: here the problem is to codify this information;

2.
The environment is completely accessible by any observer, so it is easy to collect data on the lost energy of any system;

3.
The flows cause entropy generation variations, consequently we can evaluate the entropy generation to obtain information to the flows, even when we are unable to evaluate the flows themselves;

4.
The entropy generation is a global quantity, so we can obtain global information on the cells, but from a biomedical point of view just the global cells behaviour is the useful information.
Biological systems are very interesting because they are able to adapt to the variation of environmental conditions; indeed, cells attain their “optimal” performance by a selection process driven by their environmental interactions. The resultant effect is a redistribution of energy, ions and mass flows in their metabolic network, by using regulatory proteins.
The bioengineering thermodynamic approach to biological systems consists in the analysis of the biological optimization process realized by Nature. It is no more than the classical and engineering thermodynamic analysis of the steadystate flux distribution, which, for a cell, are no more than the metabolic flows. So, starting from Equation (1) and considering the second law for the open systems [18]:
where Q is the heat exchanged, T is the temperature of the thermal source, s is the specific entropy, G is the mass flow and τ is the lifetime of the process. But, for any open system, the entropy balance in a local form results [22]:
where s = S/m, is the specific entropy, S is the entropy, σ is the entropy production density, v is the specific volume, Q is the heat flow, ẋ_{ i } is the relative velocity in relation to the centre of mass reference, and ẋ_{ B } is the centre of mass velocity. Now, considering that [22]:
where s is the specific entropy, u is the internal specific energy, v is the specific volume, p is the pressure, μ _{ i } are the chemical potentials, c _{ i } is the concentrations, T is the temperature, d/dt = ∂/∂t + ẋ_{ B } ⋅ ∇, q is the heat per unit mass, Π = P – p I with Π total pressure tensor, p hydrostatic pressure and I identity matrix of which the elements are I _{ jk } = δ _{ jk } = 1 if j = k and 0 in the other cases, a:b = Σ_{ ij } a _{ ij } b _{ ji } is the product between two tensors a and b, J _{ k } = ρ _{ k } (ẋ _{i} − ẋ_{ B }) is the diffusion flows and F _{ k } are the forces, J _{ j } is the chemical reaction rate of the jth chemical reaction and ν _{ ij } are quantities such that if they are divided by the molecular mass of the ith component they are proportional to the stoichiometric coefficients. Now, introducing the electrochemical affinity Ã = A + Z Δϕ related also to pH variation and the electric field variation, with A_{ j } = Σ_{ k } ν _{ kj } μ _{ j }, Z the electric charge per unit mass, ϕ the electrostatic potential, the relation (AA) holds [25–28]:

1.
S _{ g,tf } is the entropy generation due to the thermal flux driven by temperature difference;

2.
S _{ g,dc } is the entropy generation due to the diffusion current driven by chemical potential gradients, with \( \tilde{\mu} \) = μ + Z ϕ electrochemical potential, μ chemical potential;

3.
S _{ g,vg } is the entropy generation due to the velocity gradient coupled with viscous stress;

4.
S _{ g,cr } is the entropy generation due to the chemical reaction rate driven by affinity, always positive;

5.
S _{ g,de } is the entropy generation due to the dissipation due to work by interaction with the environment;
and τ _{ i }, i ∈ [1, 5], are the lifetimes of any process and the volume of the cell is evaluated by a characteristic length, in transport phenomena usually considered the diameter of the cell approximated as the diameter of a sphere L = (6V/π)^{1/3} = 2r, with r being the cell radius;

1.
the mean environmental temperature can be assumed as T _{0} = 310 K and the mean cell temperature has been estimated to be T _{0} + ΔT. The quantity ΔT would be experimentally evaluated for different cells lines in relation to their metabolism;

2.
the internal energy density results in u = 3.95 × 10^{7} Jm^{−3}, being calculated as the ratio between the ATP energy, U = 3 × 10^{−7} J and the mean value of the cell inside the human body, V = 7600 μm^{3}. It must be emphasized that this is an approximation because the cell volume inside the human body is in the range of 200–15000 μm^{3};

3.
the thermal molecular mean velocity inside the cytoplasm is considered to be = 5 × 10^{−5} m s^{−1};

4.
the membrane volume is calculated with \( {\mathrm{V}}_{\mathrm{m}}=\frac{4}{3}{\uppi \mathrm{r}}^3\hbox{} \frac{4}{3}\uppi {\left({\mathrm{r}\hbox{} \mathrm{d}}_{\mathrm{e}}\right)}^3=\frac{4}{3}{\uppi \mathrm{r}}^3\hbox{} \frac{4}{3}\pi {\left(r\hbox{} 0.2r\right)}^3=0.992V \) being d _{ e } = 0.2 r;

5.
the chemical potential gradient can be approximated through the ratio between the mean value of the chemical potential μ = 1.20 × 10^{−9} J kg^{−1} and the membrane length d _{ m } = 0.01 μm, with the mean density being ρ = 1000 kg m^{−3};

6.
the viscosity is taken to be 6.91 × 10^{−3} N s m^{−2};

7.
η ~ 2.07 × 10^{−3} N s m^{−2} at 30 °C;

8.
ẋ_{ B } is set as 3.0 × 10^{−6} m s^{−1};

9.
τ _{1} is the time related to the thermal flow driven by temperature difference. It can be assessed considering that the time constant of the thermal transient for heat conduction is τ _{ cv } ≈ pcV/(hA) with ρ ≈ 1000 kg m^{−3} density, V the cell volume, A the external cell surface, c ≈ 4186 J kg^{−1}K^{−1} specific heath, and h the convection heat transfer coefficient evaluated as: h ≈ 0.023Re ^{0.8} Pr ^{0.35} λ/L,where λ ≈ 0.6 W m^{−1}K^{−1} of heat conductibility, L the characteristic dimension of the cell (here we have considered the diameter), Re ≈ 0.2 the Reynolds number and Pr ≈ 7 the Prandtl number. The process would have occurred in a time τ _{1} ≈ 5 τ _{ cv }. For human cells this value can be considered in the range 15–269 ms;

10.
τ _{2} is the time related to the diffusion current driven by chemical potential gradients. It can be evaluated as τ _{2} ≈ d/D, with d = 0.01 μm, i.e., the length of the membrane, and D being the diffusion coefficient. Considering that the diffusion coefficient of glucose is approximately 10^{−9} m^{2}s^{−1} it follows that τ _{2} ≈ 10 s;

11.
τ _{3} is the time related to the velocity gradient coupled with viscous stress. This time can be evaluated as the propagating time of a mechanical wave on the surface of the cell \( {\tau}_3\approx \frac{2\pi r}{c} \) with c ~ 1540 m s^{−1} the sound velocity, considered to be the same in biological tissue;

12.
τ _{4} is the time related to the chemical reaction rate driven by affinity and it can be evaluated considering the magnitude order of a chemical reaction in a cell (~10^{−7} mol s^{−1}l^{−1}). Moreover, we consider that the moles number is proportional to the density of the chemical species (for glucose 1540 kg m^{−3}) and the volume of the cell itself. It follows that this time is in the range 17–1283 ns;

13.
τ _{5} is the time related to the dissipation due to work by interaction with the external forces. It depends on the interaction considered;

14.
L is a characteristic length, introduced as usually done in transport phenomena.
An experiment has been developed to obtain also a direct proof [31]. Therefore the spontaneous heat exchanged by the cell represents the interaction or the spontaneous communication between the cell and its environment. The proposed thermodynamic theory predicts that the temperature difference between cells with distinct metabolic characteristics can be amplified by an altered interaction with the external environment, due to the entropy generation term related to the interaction of the system with the external fields. The experiments carried out on cells exposed to low frequency electromagnetic waves consolidate the thermodynamic approach. Indeed, through infrared thermography an adimensional number, maned thermal dispersion index, was evaluated. This adimensional number represents the inability of the cells to fit their thermal power to environmental changes. Primary fibroblasts display a high dispersion index, with a maximal value of 800 % vs NIH3T3 immortalized line, which means that the primary fibroblasts adjust more efficiently their thermal production or dissipation than the NIH3T3 cells. This significant difference implies that, when exposed to selected environmental conditions, transformed cells dissipate heat more slowly than their normal counterpart. The results of this experimental approach demonstrate that selecting environmental conditions it is possible to appreciate distinct cellular phenotypes; these differences can be evaluated by thermal dispersion patterns measured by infrared thermography. The experiment confirmed the bioengineering thermodynamics theoretical results.
The results obtained can be improved by considering other approach to bioengineering thermodynamics devoted to the study of organization in living systems and by linking each others. Indeed, evolution over the long term requires a constant generation of new alternative forms, a biological behavior named mutation [35]. The cooperative effect of mutation and selection consists in different processes on different time scales:

1.
Microevolution: changes within a natural populations, in the composition of populations due to mutation and natural selection. It occurs on a time scale of generations, it represents an adaptive change;

2.
Macroevolution: changes between populations, in the composition of lineages due to speciation and extinction. It occurs in a geological time scale.
Macroevolution has its origin in microevolution as the result of natural selection acting on genotypic and phenotypic variation [35–38]. These natural processes can be described by introducing mathematical models, based on two thermodynamic actions [35, 39]:

1.
The acquisition of resources from the external environment and its conversion into energy storage;

2.
The transformation of the metabolic energy into useful work.
The bases of these processes is the interaction between biosystem and environment [40]. This brings to non equilibrium states, and the mathematical formalisms developed to the biosystems analysis was the dynamical systems, based on the studies of Bowen [41], Ruelle [42] and Sinai [43], who provided new perspectives in the analysis of far from equilibrium systems by the discovery of certain connections between nonequilibrium statistical mechanics and the ergodic theory of dynamical systems. In this context the fundamental concept is the entropy and just this concept represents the link between the dynamical systems approach and the thermodynamic approach here developed. Indeed, following Ruelle [44], considering a classical system with isokinetic time evolution described by the equation:
with
p∈R ^{N} and q∈R ^{N} momentum and position respectively, ξ a nongradient time independent force, m mass and (−α p) the isokinetic thermostat mathematical expression with α defined as:
so that [44]:
Under these conditions Ruelle defined the entropy increment as [44]
with δ _{ τ } F is a timedependent small perturbation of F, \( {\rho}_{\boldsymbol{\upxi}}\left({\nabla}_{\mathbf{x}}\left(d\mathbf{x}\circ {f}_{{\boldsymbol{\upxi}}_0}^{t\tau}\right)\cdot {\delta}_{\tau}\mathbf{F}\left(\mathbf{x}\right)\right) \) probability distribution, f _{ ξ } ^{t − τ} the solution of the equation (18) at the time tτ corresponding to the initial conditions ξ _{0}, Φ(x) = (N  1)α. Then, Denbigh [18, 45] expressed the fundamental processes of living systems, introducing an entropy approach:
where dS is the total entropy elementary variation, dS _{int} is the entropy elementary production within the system due to its metabolism of ingested exergy and dS _{ext} is the entropy exchange with the environment. Entropy is a path independent state function, and the overall reaction entropy ΔS _{ R } can be evaluated by the macroscopic reaction stoichiometry between external metabolites:
where ΔS _{ i } = − c ln p _{ i } is the entropy of reaction, s _{ li } = (h _{ li } – g _{ li })/T, with h _{ li } molar enthalpy and g _{ li } Gibbs molar energy, are the molar entropies of the k reactants and products, ν _{ l } are the stoichiometry coefficients, p _{ i } is the probability of the ith mode and c is a constant related to the numbers of elementary modes and on their reaction entropies. It represents the state of the fully evolved metabolic network [46]. When the living systems increase in organization, they increase their entropy and, far from equilibrium, they have a high exergy content [47]; indeed, considering two systems with the same mass and the same chemical composition, the one, that has a large amount of organization, has also higher exergy content. During their evolution, the living systems, and also ecosystems, increase their structure in organization, which is a working information useful for resilience and integrity, and also their efficiency in converting exergy to entropy, in order to reduce the applied exergy gradient, while their internal entropic state continue to decrease [48, 49]. Then, while dS _{int} is always positive defined (dS _{int} ≥ 0), dS _{ext} can have any sign.
The inner entropy generation rate σ is defined as the local first time derivative of the [50] internal component of the entropy:
If the irreversible processes are sufficiently slow, the Gibbs equation can be applied to any subsystem [50]:
and the entropy can be expressed in terms of fluxes J _{ i } and conjugated generalized forces X _{ i } [50]:
The nonequilibrium stationary states, which are the states whose variables are independent of time, play a fundamental role in the irreversible processes. After a characteristic time, the system achieves the equilibrium if no restraints are imposed on it, while if a number of constant restraints are imposed, a steady state is attained [50]. In any steady state the total entropy is independent of time, consequently:
but
and it is possible to argue that the entropy generation rate in a stationary system must be compensated by the liberation of entropy to the surroundings. This means also that nonequilibrium steady states cannot occur in isolated systems because these last systems do not allow exchange of entropy between the systems and the surroundings [8]. Prigogine proved that [51–53]:
On the use of the Prigogine’s results there is little doubt that a mature organism may reached a stationary state; indeed, the homeostasis of all self regulating systems is interpreted as tendency to return from a perturbed state to that of highest stability compatible with biological constraints [50].
Moreover, considering an irreversible and open system, it is composed by N elementary volumes. Every ith element of this system is located by a position vector x _{ i }, it has a velocity ẋ_{ i }, a mass m _{ i } and a momentum p _{ i } = m _{ i }ẋ_{ i }. The total mass of the system is m = ∑_{ i } m _{ i } and its density is ρ = m/V with V = ∑_{ i } V _{ i } total volume. The position of the centre of mass is x _{ B } and its velocity results ẋ_{ B } = ∑_{ i } m _{ i }ẋ_{ i }/m, while the diffusion velocity is u _{ i } = ẋ_{ i } − ẋ_{ B }. The total mass of the system is conserved, so the following relation \( \dot{\rho}+\rho \nabla \cdot {\dot{\mathbf{x}}}_B=0 \) is satisfied together with its local expression \( {\dot{\rho}}_i+\rho \nabla \cdot {\dot{\mathbf{x}}}_i=\rho {\varXi}_i \), related to the density of the ith elementary volume of density ρ _{ i } and a source Ξ generated by matter transfer, chemical reactions or thermodynamic transformations. For an open system, as just described in macroscopic way, the equation of the entropy balance is [22]:
where s = S/m, is the specific entropy, S entropy, σ the density of the entropy generation rate, v the specific volume and J _{ S } is the entropic flux defined as:
with Q heat flux.
Any dynamical state of this system can be described by the 3N canonical coordinates {x _{ i } ∈ R ^{3}, i ∈ [1,N]} and their conjugate momenta {p _{ i } ∈ R ^{3}, i ∈ [1,N]}. The 6N − dimensional space spanned by{(p _{ i },x _{ i }), i ∈ [1,N]} is the phase space Ω of the open system considered. Any point q _{ i } = (p _{ i },x _{ i }), q _{ i } ∈ R ^{6N} in the phase space Ω, represents a state of the entire N − elements system [54]. Any family {ξ(t), t ∈ R} is called stochastic process in the phase space Ω and it can be represented by a family of equivalent classes of random variables ξ(t) on Ω, {γ(σ(t)) : t ∈ R}. The point function γ(q(t)) is called trajectory of the stochastic process ξ(t): a description of a physical system in terms of a trajectory of a stochastic process corresponds to a point dynamics, while its description in terms of equivalent classes of trajectories and their associated probability measure corresponds to an ensemble dynamics [55]. So it is considered a nonequilibrium system moving in the Ωspace between two states, which are in two elementary cells of a given partition of the phase space. We use the concept of path of classical mechanics: if the motion of the system is regular, or if the phase manifold has positive or zero Riemannian curvature, there will be only a fine bundle of paths which track each other between the initial and the final cells [13]. For a system in chaotic motion, or when the Riemannian curvature of the phase manifold is negative, two points indistinguishable in the initial cell can separate from each other exponentially [54]. Then, between two given phase cells, there may be many possible paths γ _{ k }, k ∈ [1,ω] with ω number of all the paths, with different travelling time t _{ γk } of the system and different probability p _{ γk } for the system to take the path k, called path probability distribution [56–59]. It is considered an ensemble of a large number L of identical systems, all moving in the phase space from two cells with ω possible paths, and L _{ k } systems travelling on the path γ _{ k }. The probability p _{ γk } that the system take the path γ _{ k } is defined as usual by p _{ γk } = L _{ k }/L. If ω _{ k } = 1 then p _{ γk } = 1. By definition, p _{ γk } is the transition probability from the two states considered. These trajectories must be the paths minimizing action according to the principle of least action [54]. Since 1962, Jaynes argued that Gibbs’ formalism of equilibrium statistical mechanics could be generalised in a statistical inference theory for nonequilibrium systems [60]. Jaynes developed the nonequilibrium statistical mechanics for the stationary state constraint on the basis of maximum entropy; his approach consists of maximising the path Shannon information entropy written for the path, S _{ I } = − Σ_{ γ } p _{ γ } ln p _{ γ }, with respect to p _{ γ } of the path γ, with the probability subject to the actual constraints. According to Shannon, ‘the information entropy is the logarithm of the number of the outcomes i with nonnegligible probability p _{ i }’, while in ‘nonequilibrium statistical mechanics it is the logarithm of the number of microscopic phasespace paths γ having nonnegligible probability p _{ γ }’ [60]. Jaynes’ approach consists of finding the ‘most probable macroscopic path realised by the greater number of microscopic paths compatible with the imposed constrained’ [60], in analogy with the Boltzmann microstate counting: ‘paths rather then states are the central objects of interest in nonequilibrium systems, because of the presence of nonzero macroscopic fluxes whose statistical description requires considering the underlying microscopic behaviour over time’ [60] which implies that ‘the macroscopic behaviour is reproducible under given constraints’ and it is ‘characteristic of each of the great number of microscopic paths compatible with those constraints’ [60]. Following this approach and these considerations, the statistical expression of the entropy generation has been written as [56–59]:
It can be also interpreted as the missing information necessary for predicting which path a system of the ensemble takes during the transition from a state to another.
In the theory of probability the stochastic order is introduced. Two random variables X and Y are in stochastic order if there exists a random variable Z and functions ψ _{1} and ψ _{2} such that X = ψ _{1}(Z) and Y = ψ _{2}(Z), with ψ _{1}(Z) ≤ ψ _{2}(Z) [61]. Now, the set of paths {γ _{ k }, k ∈ [1,ω]} is considered, with ω number of all the paths between two thermodynamic states, represented by two points in the phase space. It is possible to define a stochastic order among the paths, saying that a path γ _{ i } is stochastically smaller than a path γ _{ j } if its probability p _{ γi } is smaller that the probability of the other path, p _{ γj } [13]:
So, the probability of a path can be expressed in term of the first order differential of the entropy generation respect to the probability itself, as follows [13]:
But, in the analysis of the complex systems, it was highlighted how chaotic and fractal behaviour are very widespread in nature: any numerical evaluation based on accessible states in the phase space is incomplete because of the rejected, singular or inaccessible points [56]. The basis of the incomplete information is that a part of information on complex system may not be completely accessible. The consequence is that irreversibility is the physical model by which thermodynamic phenomena can be completely described [54]. The related information is incomplete because, for complex systems, it occurs that ∑ _{ j = 1} ^{ν} p _{ j } = θ ≤ 1, with ν number of accessible or accountable states, smaller of the total number of states, and θ incompleteness of the treatment and linked to the nature of the system, consequence of the partial knowledge of the dynamics or of the inaccessible states of the system itself [54]. Nonstatistical mechanics replaces the complete probability normalisation by:
with ϑ incompleteness parameter such that ϑ = 1 if the probability distribution is complete. It can be related to the incompleteness θ by the following relation [13]:
The phase space cells, which represent the stationary states, was proven to form a subset of all the cells on which the evolution acts as a onecycle permutation: this kind of ergodicity has been defined ergodicity for irreversibility [54]. Moreover, in nonequilibrium transformation, the volume of the phase space was proven to contract indefinitely [54]. Recently, considering the expression for the probability p _{ γi } of a path γ _{ i } and the statistical results on the entropy generation [54], it was proven that [13]:
which means that the paths are statistically ordered. The stochastic order of the path proves that the evolution of the biosystems is related to their irreversibility and the quantity useful to evaluate the allowed paths and their probability is the entropy generation. Consequently, a link between the bioengineering thermodynamic approach proposed and the dynamical system approach is obtained.
References
 1.
Demirel Y, Sandler SI. Thermodynamics and bioenergetics. Biophys Chem. 2002;97:87–111.
 2.
Toussaint O, Schneider ED. The thermodynamic and evolution of complexity in biological systems. Comp Biochem Physiol A. 1998;120:3–9.
 3.
Caplan SR, Essig A. Bioenergetics and Linear Nonequilibrium Thermodynamics, The Steady State. Cambridge: Harvard University Press; 1983.
 4.
Lucia U. Entropy generation approach to cell systems. Physica A. 2014;406:1–11.
 5.
Lucia U. Molecular machine as chemicalthermodynamic devices. Chem Phys Lett. 2013;556:242–4.
 6.
Salerian AJ, Saleri NG. Cooling core body temperature may slow down neurodegeneration. CNS Spectr. 2008;13:227–9.
 7.
Katchalsky A, Curran PF. Nonequilibrium thermodynamics in biophysics. Cambridge: Harvard University Press; 1967.
 8.
Lucia U. Irreversibility in biophysical and biochemical engineering. Physica A. 2012;391:5997–6007.
 9.
Bustamante C, Chemla YR, Forde NR, Izhaky D. Mechanical processes in biochemistry. Annu Rev Biochem. 2004;73:705–48.
 10.
Hudspeth A, Choe Y, Mehta A, Martin P. Putting ion channels to work: mechanoelectrical transduction, adaptation and amplication by hair cells. Proc Natl Acad Sci. 2000;97(22):11765.
 11.
Lucia U. Irreversibility, entropy and incomplete information. Physica A. 2009;388(19):4025–33.
 12.
Lucia U. (2013b) Stationary open systems: A brief review on contemporary theories on irreversibility. Physica A. 2013;392(5):1051–62.
 13.
Lucia U. Thermodynamic paths and stochastic order in open systems. Physica A. 2013;392(18):3912–9.
 14.
Lucia U. Thermodynamics and cancer stationary states. Physica A. 2013;392(17):3648–53.
 15.
Lucia U. Transport processes in biological systems: tumoral cells and human brain. Physica A. 2014;393:327–36.
 16.
Denbigh KG. Note on entropy, disorder and disorganization. Brit J Phil Sci. 1989;40:323–32.
 17.
Denbigh KG. The many faces of irreversibility. Brit J Phil Sci. 1989;40:501–18.
 18.
Bejan A. Advance engineering thermodynamics. New York: John Wiley; 2006.
 19.
Lucia U. Some considerations on molecular machines and Loschmidt paradox. Chem Phys Lett. 2015;623:98–100.
 20.
Lucia U. A link between nano and classical thermodynamics: dissipation analysis (the entropy generation approach in nanothermodynamics). Entropy. 2015;17(2):1309–28.
 21.
Lucia U. Entropy production and generation: clarity from nanosystems considerations. Chem Phys Lett. 2015;629:87–90.
 22.
de Groot SG, Mazur P. Nonequilibrium thermodynamics. Amsterdam: NorthHolland Publishing; 1984.
 23.
Zemansky MW. Heat and thermodynamics. New York: McGrawHill; 1966.
 24.
Atkins PW. The elements of physical chemistry. 3rd ed. Oxford: Oxford University Press; 1993.
 25.
Lucia U. Entropy generation and cell growth with comments for a thermodynamic anticancer approach. Physica A. 2014;406:107–18.
 26.
Lucia U. Thermodynamic approach to nanoproperties of cell membrane. Physica A. 2014;407:185–91.
 27.
Lucia U. Transport processes and irreversible thermodynamics analysis in tumoral systems. Physica A. 2014;410:380–90.
 28.
Lucia U. The gouystodola theorem in bioenergetic analysis of living systems (Irreversibility in bioenergetics of living systems). Energies. 2014;7:5717–39.
 29.
Lucia U, Ponzetto A, Deisboeck TS. A thermophysical analysis of the proton pump vacuolarATPase: the constructal approach. Sci Rep. 2014;4:6763.
 30.
Lucia U, Ponzetto A, Deisboeck TS. A thermodynamic approach to the ‘mitosis/apoptosis’ ratio in cancer. Physica A. 2015;436:246–55.
 31.
Lucia U, Grazzini G, Montrucchio B, Grisolia G, Borchiellini R, Gervino G, et al. Constructal thermodynamics combined with infrared experiments to evaluate temperature differences in cells. Sci Rep. 2015;5:11587.
 32.
Tuszynski JA, Kurzynski M. Introduction to molecular biophysics. Raton: Taylor & Francis; 2003.
 33.
Newman J. Physics of the life sciences. Berlin: Springer; 2008.
 34.
Kittel C, Kroemer H. Thermal physics. London: W.H. Freeman and Company; 1980.
 35.
Demetrius LA, Gundlach VM. Directionality theory and the entropic principle in natural selection. Entropy. 2014;16:5428–522.
 36.
Mayr E. What evolution is. New York: Weidenfeld and Nicholson; 2002.
 37.
Mayr E. Evolution and the diversity of life. Harvard: Harvard University Press; 1976.
 38.
Demetrius L. Directionality principles in thermodynamics and evolution. Proc Natl Acad Sci U S A. 1997;94:3491–8.
 39.
Lehninger A. Bioenergetics. New York: WA Benjamin Inc.; 1965.
 40.
Lewontin RC. Gene, organism and environment. In: Bendall DS, editor. Evolution from molecules to men. Cambridge: Cambridge University Press; 1983. p. 273–85.
 41.
Bowen R. Equilibrium states and the ergodic theory of anosov diffeomorphisms. Vol. 470 lecture notes in math. New York: Springer; 1975.
 42.
Ruelle D. Thermodynamic formalism. Vol. 5 encyclopedia of mathematics and its applications. Reading: AddisonWesley; 1978.
 43.
Sinai YG. Gibbs measures in ergodic theory. Russ Math Surv. 1972;4:21–69.
 44.
Ruelle D. Extending the definition of entropy to nonequilibrium steady states. Proc Natl Acad Sci U S A. 2003;100:3054–8.
 45.
Denbigh KG. Note on entropy, disorder and disorganization. Brit J Phil Sci. 1989;40:323–32.
 46.
Sandler SI, Orbey H. On the thermodynamics of microbialgrowth processes. Biotech Bioeng. 1991;38:697–718.
 47.
Swenson R. Emergent attractors and the law of maximum entropy prpduction: foundations to a theory of general evolution. Systems Research. 1989;6:187–97.
 48.
Brooks DR, Collier J, Maurer BA, Smith JDH, Wiley EO. Entropy and information in evolving biological systems. Biology & Philosophy. 1989;4:407–32.
 49.
Günther F, Folke C. Characteristics of nested living systems. J Biological Systems. 1993;1:257–74.
 50.
Katchalsky A, Kedem O. Thermodynamics of flow processes in biological systems. Biophys J. 1962;2:53–78.
 51.
Glansdorff P, Prigogine I. Thermodynamic theory of structure, stability and fluctuations. New York: John Wiley & Sons; 1971.
 52.
Prigogine I. Etude Thermodynamique des Phénomènes Irrèversibles. Liège: Desoer; 1947.
 53.
Prigogine I. Introduction to thermodynamics of irreversible processes. New York: Interscience; 1961.
 54.
Lucia U. Probability, ergodicity, irreversibility and dynamical systems. Proc Royal Soc A. 2008;464:1089–184.
 55.
Primas H. Basic elements and problems of probability theory. J Sci Explor. 1999;13(4):579–613.
 56.
Lucia U. Statistical approach of the irreversible entropy variation. Physica A. 2008;387(14):3454–60.
 57.
Lucia U. Irreversibility, entropy and incomplete information. Physica A. 2009;388:4025–33.
 58.
Lucia U. Maximum entropy generation and κ − exponential model. Physica A. 2010;389:4558–63.
 59.
Lucia U. Irreversibility entropy variation and the problem of the trend to equilibrium. Physica A. 2007;376:289–92.
 60.
Dewar R. Information theory explanation of the fluctuation theorem, maximum entropy production and selforganized criticality in nonequilibrium stationary states. J Phys A: Math Gen. 2003;36:631–41.
 61.
Shaked M. Stochastic orders. New York: Springer; 2006.
Acknowledgments
The author must thank prof. Antonio Ponzetto (Università di Torino) for his scientific support in biomedical interpretation of the thermodynamic results and Prof. Bartolomeo Montrucchio (Politecnico di Torino) for his support in the design of the experimental devices.
Author information
Additional information
Competing interests
The author declare that he has no competing interests.
Authors’ contributions
UL has developed the thermodynamics approach and written the paper.
Authors’ information
Umberto Lucia, MSc in Physics at Turin University and PhD in Energetics at Florence University, teaches Engineering thermodynamics at Biomedical Engineering School at Politecnico di Torino. He is member of the Energy Department and develops researches in thermodynamics biosystems, and classical, statistical and quantum thermodynamics of complex systems. After twenty years spent to study the relations between irreversibility and steady states in open systems, in the last three years he developed the bioengineering thermodynamics applied to biosystems, biological cells, cancer and non equilibrium states. This approach to biosystems has been experimentally confirmed. Consequently, he continue to improve the theory developed in order to obtain a full thermodynamic approach useful to predict the behaviour of the biological cells with particular interest to cancer and its thermodynamic control.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Lucia, U. Bioengineering thermodynamics of biological cells. Theor Biol Med Model 12, 29 (2015). https://doi.org/10.1186/s129760150024z
Received:
Accepted:
Published:
Keywords
 Entropy generation
 Exergy
 Irreversibility
 Living systems
 Medicine and biochemistry thermodynamics
 Membrane transport