# Bioengineering thermodynamics of biological cells

- Umberto Lucia
^{1}Email authorView ORCID ID profile

**12**:29

https://doi.org/10.1186/s12976-015-0024-z

© Lucia. 2015

**Received: **6 August 2015

**Accepted: **11 November 2015

**Published: **1 December 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, thermo-electro-chemical 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 thermo-electro-biochemical 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.

### Keywords

Entropy generation Exergy Irreversibility Living systems Medicine and biochemistry thermodynamics Membrane transport## 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, nano-mechanical 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].

- 1.
An open irreversible real linear or non-linear 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 well-known (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.

- 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 flow-based 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.

*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:

**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]:

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.

**J**

_{ k }is the flow of the

*k*-th quantity involved in the process considered and

**X**

_{ k }is the related thermodynamic force. Now, considering that:

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 electro-magnetical 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 electro-diffusion 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 ATP-driven 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].

*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:

*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]:

*σ*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 cross-membrane 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].

^{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

- 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.

*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]:

*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]:

*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

*j*-th chemical reaction and

*ν*

_{ ij }are quantities such that if they are divided by the molecular mass of the

*i*-th component they are proportional to the stoichiometric coefficients. Now, introducing the electro-chemical 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;

*τ*

_{ 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*= (6

*V*/

*π*)

^{1/3}= 2

*r*, 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.023*Re*^{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.

- 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.

- 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.

**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:

*δ*

_{ τ }

**F**is a time-dependent 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:

*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 *i*-th 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.

*σ*is defined as the local first time derivative of the [50] internal component of the entropy:

*J*

_{ i }and conjugated generalized forces

*X*

_{ i }[50]:

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].

*N*elementary volumes. Every

*i*-th 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

*i*-th 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]:

*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.

*N*canonical coordinates {

**x**

_{ i }∈

*R*

^{3},

*i*∈ [1,

*N*]} and their conjugate momenta {

**p**

_{ i }∈

*R*

^{3},

*i*∈ [1,

*N*]}. The 6

*N*− 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 non-equilibrium 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 non-equilibrium systems [60]. Jaynes developed the non-equilibrium 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 non-negligible probability

*p*

_{ i }’, while in ‘non-equilibrium statistical mechanics it is the logarithm of the number of microscopic phase-space paths

*γ*having non-negligible 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 non-equilibrium systems, because of the presence of non-zero 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.

*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]:

_{ 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]. Non-statistical mechanics replaces the complete probability normalisation by:

*ϑ*incompleteness parameter such that

*ϑ*= 1 if the probability distribution is complete. It can be related to the incompleteness

*θ*by the following relation [13]:

*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 bio-systems 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.

## Declarations

### 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.

**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.

## Authors’ Affiliations

## References

- Demirel Y, Sandler SI. Thermodynamics and bioenergetics. Biophys Chem. 2002;97:87–111.View ArticlePubMedGoogle Scholar
- Toussaint O, Schneider ED. The thermodynamic and evolution of complexity in biological systems. Comp Biochem Physiol A. 1998;120:3–9.View ArticleGoogle Scholar
- Caplan SR, Essig A. Bioenergetics and Linear Nonequilibrium Thermodynamics, The Steady State. Cambridge: Harvard University Press; 1983.View ArticleGoogle Scholar
- Lucia U. Entropy generation approach to cell systems. Physica A. 2014;406:1–11.View ArticleGoogle Scholar
- Lucia U. Molecular machine as chemical-thermodynamic devices. Chem Phys Lett. 2013;556:242–4.View ArticleGoogle Scholar
- Salerian AJ, Saleri NG. Cooling core body temperature may slow down neurodegeneration. CNS Spectr. 2008;13:227–9.PubMedGoogle Scholar
- Katchalsky A, Curran PF. Nonequilibrium thermodynamics in biophysics. Cambridge: Harvard University Press; 1967.Google Scholar
- Lucia U. Irreversibility in biophysical and biochemical engineering. Physica A. 2012;391:5997–6007.View ArticleGoogle Scholar
- Bustamante C, Chemla YR, Forde NR, Izhaky D. Mechanical processes in biochemistry. Annu Rev Biochem. 2004;73:705–48.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Lucia U. Irreversibility, entropy and incomplete information. Physica A. 2009;388(19):4025–33.View ArticleGoogle Scholar
- Lucia U. (2013b) Stationary open systems: A brief review on contemporary theories on irreversibility. Physica A. 2013;392(5):1051–62.View ArticleGoogle Scholar
- Lucia U. Thermodynamic paths and stochastic order in open systems. Physica A. 2013;392(18):3912–9.View ArticleGoogle Scholar
- Lucia U. Thermodynamics and cancer stationary states. Physica A. 2013;392(17):3648–53.View ArticleGoogle Scholar
- Lucia U. Transport processes in biological systems: tumoral cells and human brain. Physica A. 2014;393:327–36.View ArticleGoogle Scholar
- Denbigh KG. Note on entropy, disorder and disorganization. Brit J Phil Sci. 1989;40:323–32.View ArticleGoogle Scholar
- Denbigh KG. The many faces of irreversibility. Brit J Phil Sci. 1989;40:501–18.View ArticleGoogle Scholar
- Bejan A. Advance engineering thermodynamics. New York: John Wiley; 2006.Google Scholar
- Lucia U. Some considerations on molecular machines and Loschmidt paradox. Chem Phys Lett. 2015;623:98–100.View ArticleGoogle Scholar
- Lucia U. A link between nano- and classical thermodynamics: dissipation analysis (the entropy generation approach in nano-thermodynamics). Entropy. 2015;17(2):1309–28.View ArticleGoogle Scholar
- Lucia U. Entropy production and generation: clarity from nanosystems considerations. Chem Phys Lett. 2015;629:87–90.View ArticleGoogle Scholar
- de Groot SG, Mazur P. Non-equilibrium thermodynamics. Amsterdam: North-Holland Publishing; 1984.Google Scholar
- Zemansky MW. Heat and thermodynamics. New York: McGraw-Hill; 1966.Google Scholar
- Atkins PW. The elements of physical chemistry. 3rd ed. Oxford: Oxford University Press; 1993.Google Scholar
- Lucia U. Entropy generation and cell growth with comments for a thermodynamic anticancer approach. Physica A. 2014;406:107–18.View ArticleGoogle Scholar
- Lucia U. Thermodynamic approach to nano-properties of cell membrane. Physica A. 2014;407:185–91.View ArticleGoogle Scholar
- Lucia U. Transport processes and irreversible thermodynamics analysis in tumoral systems. Physica A. 2014;410:380–90.View ArticleGoogle Scholar
- Lucia U. The gouy-stodola theorem in bioenergetic analysis of living systems (Irreversibility in bioenergetics of living systems). Energies. 2014;7:5717–39.View ArticleGoogle Scholar
- Lucia U, Ponzetto A, Deisboeck TS. A thermo-physical analysis of the proton pump vacuolar-ATPase: the constructal approach. Sci Rep. 2014;4:6763.PubMed CentralView ArticlePubMedGoogle Scholar
- Lucia U, Ponzetto A, Deisboeck TS. A thermodynamic approach to the ‘mitosis/apoptosis’ ratio in cancer. Physica A. 2015;436:246–55.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Tuszynski JA, Kurzynski M. Introduction to molecular biophysics. Raton: Taylor & Francis; 2003.View ArticleGoogle Scholar
- Newman J. Physics of the life sciences. Berlin: Springer; 2008.View ArticleGoogle Scholar
- Kittel C, Kroemer H. Thermal physics. London: W.H. Freeman and Company; 1980.Google Scholar
- Demetrius LA, Gundlach VM. Directionality theory and the entropic principle in natural selection. Entropy. 2014;16:5428–522.View ArticleGoogle Scholar
- Mayr E. What evolution is. New York: Weidenfeld and Nicholson; 2002.Google Scholar
- Mayr E. Evolution and the diversity of life. Harvard: Harvard University Press; 1976.Google Scholar
- Demetrius L. Directionality principles in thermodynamics and evolution. Proc Natl Acad Sci U S A. 1997;94:3491–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Lehninger A. Bioenergetics. New York: WA Benjamin Inc.; 1965.Google Scholar
- Lewontin RC. Gene, organism and environment. In: Bendall DS, editor. Evolution from molecules to men. Cambridge: Cambridge University Press; 1983. p. 273–85.Google Scholar
- Bowen R. Equilibrium states and the ergodic theory of anosov diffeomorphisms. Vol. 470 lecture notes in math. New York: Springer; 1975.Google Scholar
- Ruelle D. Thermodynamic formalism. Vol. 5 encyclopedia of mathematics and its applications. Reading: Addison-Wesley; 1978.Google Scholar
- Sinai YG. Gibbs measures in ergodic theory. Russ Math Surv. 1972;4:21–69.View ArticleGoogle Scholar
- Ruelle D. Extending the definition of entropy to nonequilibrium steady states. Proc Natl Acad Sci U S A. 2003;100:3054–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Denbigh KG. Note on entropy, disorder and disorganization. Brit J Phil Sci. 1989;40:323–32.View ArticleGoogle Scholar
- Sandler SI, Orbey H. On the thermodynamics of microbial-growth processes. Biotech Bioeng. 1991;38:697–718.View ArticleGoogle Scholar
- Swenson R. Emergent attractors and the law of maximum entropy prpduction: foundations to a theory of general evolution. Systems Research. 1989;6:187–97.View ArticleGoogle Scholar
- Brooks DR, Collier J, Maurer BA, Smith JDH, Wiley EO. Entropy and information in evolving biological systems. Biology & Philosophy. 1989;4:407–32.View ArticleGoogle Scholar
- Günther F, Folke C. Characteristics of nested living systems. J Biological Systems. 1993;1:257–74.View ArticleGoogle Scholar
- Katchalsky A, Kedem O. Thermodynamics of flow processes in biological systems. Biophys J. 1962;2:53–78.PubMed CentralView ArticlePubMedGoogle Scholar
- Glansdorff P, Prigogine I. Thermodynamic theory of structure, stability and fluctuations. New York: John Wiley & Sons; 1971.Google Scholar
- Prigogine I. Etude Thermodynamique des Phénomènes Irrèversibles. Liège: Desoer; 1947.Google Scholar
- Prigogine I. Introduction to thermodynamics of irreversible processes. New York: Interscience; 1961.Google Scholar
- Lucia U. Probability, ergodicity, irreversibility and dynamical systems. Proc Royal Soc A. 2008;464:1089–184.View ArticleGoogle Scholar
- Primas H. Basic elements and problems of probability theory. J Sci Explor. 1999;13(4):579–613.Google Scholar
- Lucia U. Statistical approach of the irreversible entropy variation. Physica A. 2008;387(14):3454–60.View ArticleGoogle Scholar
- Lucia U. Irreversibility, entropy and incomplete information. Physica A. 2009;388:4025–33.View ArticleGoogle Scholar
- Lucia U. Maximum entropy generation and
*κ*− exponential model. Physica A. 2010;389:4558–63.View ArticleGoogle Scholar - Lucia U. Irreversibility entropy variation and the problem of the trend to equilibrium. Physica A. 2007;376:289–92.View ArticleGoogle Scholar
- Dewar R. Information theory explanation of the fluctuation theorem, maximum entropy production and self-organized criticality in non-equilibrium stationary states. J Phys A: Math Gen. 2003;36:631–41.View ArticleGoogle Scholar
- Shaked M. Stochastic orders. New York: Springer; 2006.Google Scholar