Next Article in Journal
Radially Excited AdS5 Black Holes in Einstein–Maxwell–Chern–Simons Theory
Previous Article in Journal
Fluctuation of Information Entropy Measures in Cell Image
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biological Aging and Life Span Based on Entropy Stress via Organ and Mitochondrial Metabolic Loading

Department of Mechanical Engineering, Texas A&M University, College Station, TX 77843-3123, USA
*
Author to whom correspondence should be addressed.
Entropy 2017, 19(10), 566; https://doi.org/10.3390/e19100566
Submission received: 7 July 2017 / Revised: 21 September 2017 / Accepted: 7 October 2017 / Published: 23 October 2017
(This article belongs to the Section Thermodynamics)

Abstract

:
The energy for sustaining life is released through the oxidation of glucose, fats, and proteins. A part of the energy released within each cell is stored as chemical energy of Adenosine Tri-Phosphate molecules, which is essential for performing life-sustaining functions, while the remainder is released as heat in order to maintain isothermal state of the body. Earlier literature introduced the availability concepts from thermodynamics, related the specific irreversibility and entropy generation rates to metabolic efficiency and energy release rate of organ k, computed whole body specific entropy generation rate of whole body at any given age as a sum of entropy generation within four vital organs Brain, Heart, Kidney, Liver (BHKL) with 5th organ being the rest of organs (R5) and estimated the life span using an upper limit on lifetime entropy generated per unit mass of body, σM,life. The organ entropy stress expressed in terms of lifetime specific entropy generated per unit mass of body organs (kJ/(K kg of organ k)) was used to rank organs and heart ranked highest while liver ranked lowest. The present work includes the effects of (1) two additional organs: adipose tissue (AT) and skeletal muscles (SM) which are of importance to athletes; (2) proportions of nutrients oxidized which affects blood temperature and metabolic efficiencies; (3) conversion of the entropy stress from organ/cellular level to mitochondrial level; and (4) use these parameters as metabolism-based biomarkers for quantifying the biological aging process in reaching the limit of σM,life. Based on the 7-organ model and Elia constants for organ metabolic rates for a male of 84 kg steady mass and using basic and derived allometric constants of organs, the lifetime energy expenditure is estimated to be 2725 MJ/kg body mass while lifetime entropy generated is 6050 kJ/(K kg body mass) with contributions of 190; 1835.0; 610; 290; 700; 1470 and 95 kJ/K contributed by AT-BHKL-SM-R7 to 1 kg body mass over life time. The corresponding life time entropy stresses of organs are: 1.2; 60.5; 110.5; 110.5; 50.5; 3.5; 3.0 MJ/K per kg organ mass. Thus, among vital organs highest stress is for heart and kidney and lowest stress is for liver. The 5-organ model (BHKL and R5) also shows similar ranking. Based on mitochondrial volume and 5-organ model, the entropy stresses of organs expressed in kJ/K per cm3 of Mito volume are: 12,670; 5465; 2855; 4730 kJ/cm3 of Mito for BHKL indicating brain to be highly stressed and liver to be least stressed. Thus, the organ entropy stress ranking based on unit volume of mitochondria within an organ (kJ/(K cm3 of Mito of organ k)) differs from entropy stress based on unit mass of organ. Based on metabolic loading, the brains of athletes already under extreme mitochondrial stress and facing reduced metabolic efficiency under concussion are subjected to more increased stress. In the absence of non-intrusive measurements for estimating organ-based metabolic rates which can serve as metabolism-based biomarkers for biological aging (BA) of whole body, alternate methods are suggested for estimating the biological aging rate.

1. Introduction and Rationale

Living Biological systems (BS) are open systems exchanging mass, energy and entropy with the surroundings. The required energy for sustaining the life is released through the oxidation of glucose (carbohydrates, CH), fat (F) and proteins (P) (called macro-nutrients) carried by the blood stream to all the organs (a group of tissues which are self-contained and are specialized in performing a particular function) and a part of the chemical energy of nutrients is converted into another form of chemical energy, which is stored within Adenosine Tri-Phosphate (ATP) molecules, called the energy currency of the body and which can be readily used to deliver work (W, organized energy; e.g., lifting weights and placing them on shelve resulting in energy transfer in the form of work, which causes potential energy gain of weight) and/or supply cellular energy for endothermic biological reactions, active transport of species against adverse gradient (e.g., transport of glucose in intestines) and for repair/reproduce of cells for sustaining life functions while the remainder of the chemical energy is released as heat (Q, random energy) to overcome heat loss from the body and to maintain the normal body temperature.
One may define aging as an indication of decrease in physiological life sustaining functions, which leads to an increase in age related mortalities along with a decrease in reproductive functions [1]. Weinert et al. [2] presented a comprehensive review on aging and listed multiple trajectories for aging: evolutionary [3], molecular, cellular and whole body (system) levels including the environmental factors. Medvedev [4] listed more than 300 hypotheses on aging. The evolutionary concept outlined in [3] deals with genetic variants which cause Alzheimerʹs disease are found to be less in people with longer lifespans (i.e., slower rate of aging), indicating that Darwin’s natural selection at work allowing better ones (without genetic mutants) to survive. Lee et al. [5] studied the regulatory effects of nutrients on aging and suggested that the nutrient composition may be altered rather than using CR for improving life span. Kirkwood’s theory on aging falls under evolutionary and cellular [6] trajectories and mostly relies on copying errors. The ROS theory falls under cellular mechanism. Under ROS theory, the aging is attributed to accumulated damage to the cells of various organs via the production of radical oxygen species (ROS) during oxidation process/electron transport chain, resulting in impaired functions of cells within organs in repairing and replacing dead cells and eventually result in its inability to overcome the adverse environmental factors. It is cautioned that there are many other factors and interventions that modulate speed of aging but not necessarily related to ROS production [7]. Thus, aging is a “multi-factor” process.
Rubner’s rate of living theory (ROL) assumes a fixed amount of specific metabolic energy release. qM,life (≈ 835 MJ/kg body mass, a first law hypothesis, “live fast, die young” [2]) over life span irrespective of metabolic efficiency (i.e., proportional to ATP production per unit mole of oxygen consumption) for every living organism and falls under systemic trajectory. The biological aging rate (BAR) is based on how fast or slow one reaches qM,life. Living systems are characteristics of systems being far from equilibrium and hence must generate entropy since birth/conception. Thus Silva and Annamalai [8,9] adopted the availability concepts from thermodynamics, showing that availability/exergetic efficiency in thermodynamics is almost same as metabolic efficiency (ηM), proposed the second law hypothesis on biological aging with a fixed amount of life span entropy generation/irreversibility, σM,life (10,000 kJ/(K kg body mass)) for the whole body. Such a hypothesis is also consistent with Andresan’s concept on constant total lifetime entropy production per unit body mass for all BS [10] except that Silva’s hypothesis includes the effect of metabolic efficiency. Unlike Rubner’s hypothesis, the energy release along with metabolic efficiency (i.e., ATP production efficiency) affects the entropy generation rate. Thus, BAR is based on how fast or slow one reaches σM,life. The σM,life depends on the heat part, “QM” (= MJ of energy released as heat/kg body mass) of specific energy release rate, (SERR) or called as specific metabolic rate (SMR) in biology.
Apart from thermodynamics governing internal biological processes, there is a perpetual outflow of energy in the form of heat loss, Q and hence disposal of entropy generated within the whole body in the form of heat to the environment. In layman’s language, the entropy generation within a system is a measure of how things can go wrong irreversibly [11].
There is striking similarity between the field of combustion science, which deals with oxidation at high temperature (order of 1200–1500 °C, without use of catalysts) [12], and the field of metabolism in BS where oxidation at low temperature (37 °C) is aided by enzymes. Thus, currently efforts are under way to translate, modify and link the extensive results from fields of thermodynamics and combustion science to the field of biology particularly in the following areas:
  • First and second Law and availability analyses of oxidation of nutrients: CH, F and P and their mixtures and the relation to life span of BS.
  • Adiabatic temperature rise in combustion of fuels [12] vs. maximum temperature rise of blood leaving the organs either from basic principles of nutrient oxidation or from allometric laws.
  • Oxidation of a single carbon particle and carbon particle temperature vs. oxidation of nutrients in cells and single cell temperature.
  • Adoption of literature on SERR from oxygen deficient carbon cloud in engineering [7] to specific energy release rate SERRk (W/(kg of organ k)) or specific metabolic rate of organ k (SMRk) (k = Brain, heart, kidney, etc.), and bridging of gap between data on body mass independent SMRk data of Elia [8] (but still depends on type of organ or its enzyme) with those of body mass-dependent SMRk data of Wang [13] and Singer’s data [10] in biology. The SMRk will be a function of size of organ or mass of organ k when one of the following conditions are satisfied: (1) diffusion of oxygen dominates for large organ (i.e., there are a large number of cells within organ or abundant active enzymes are present in mitochondria); (2) the energy release rate is kinetically limited with first order (Michaelis Menten (MM)) kinetics (enzymes are limited or not active) SMRk is independent of organ mass under zero order kinetics or when organ size is small which leads to Elia constant in biology [13].
  • Deduction of allometric laws/constant for organs and the whole-body Kleiber’s law using Combustion Science and Thermodynamics [9,14,15,16].
The present work concerns the extension of availability analyses to mitochondrial level (Item 1), estimation of maximum possible temperature rise of blood from each organ (Item 2; i.e., non-ATP process) and presents methods of determining BAR. A brief literature review is presented summarizing previous work, followed by objectives and methodology adopted, a summary of results, and basic and derived allometric constants in tabular forms, quantitative results for ranking entropy stresses of organs and finally methods of estimating BAR.

2. Literature Review and Background

2.1. Second Law, Irreversibility Rate ( I ˙ ) and Entropy Generation Rate ( σ ˙ )

Entropy generation occurs within the completely warm body (with control surface of control volume (CV) boundary just outside skin) due to internal and external irreversibilities. The internal irreversibilities include: (1) cells of organ k (σcl,m,k, J/(K kg of organ k) where vast number of irreversible chemical reaction and oxidation reactions occurs; (2) Interstitial fluid around cells where property gradients exist (O2, CO2, glucose, temperature); (3) Skeletal muscles where ATP is used to move muscle fibers against stationary muscles involving frictional process and heart muscles due to pumping action; and (4) pressure losses in circulation system, frictional/duct losses in breathing in and out etc. The external irreversibility occurs due to the property gradients between the skin of body and surroundings; they also influence internal irreversibilities; e.g., increased heat loss to colder ambience requires increased energy release rate ( q · ) thus increasing internal irreversibility. Of all the entropy generation sources, the major contribution to whole body entropy generation comes from internal irreversibilities where almost 2/3 of ERR is released as heat ( Q · ) . If control surface of CV is selected just inside the skin, the body is almost isothermal at body temperature TB. The life-sustaining gradients (e.g., ∆g, Gibbs function difference required for reactions) within the body disappear when a biological system dies and decays into its non-living form. At this juncture or near death, rate of entropy generation within body approaches almost zero and cumulative specific entropy generated by BS, σM,life (J/(K kg body mass)) over life span reaches a maximum. This should not be confused with Smax principle for isolated system in engineering (i.e., fixed internal energy (U), volume (V) and mass (m) or fixed enthalpy (H), pressure (P), m), since BS does not have fixed U, V and m and the BS is not an isolated system. Rather it operates at fixed T and P. There were two approaches in estimating specific entropy generation of whole body (σM,life, kJ/(K kg body)):
  • Homogeneous where one considers whole body as homogeneous system and
  • Heterogeneous where the metabolic rate of whole body is considered to be a sum of metabolic rates of internal organs with each organ serving as system.
The two approaches are used in biology for estimating specific metabolic energy release (SMR, kJ/kg body) and specific entropy generation ( σ M , kJ/(K kg body)).

2.2. Homogeneous Approach

The thermodynamics literature present the entropy balance equation for an open engineering system (where mass crosses the boundary) for whole body of mass mB as [12]:
dS dt   =   Q · T b   + m · i s i m · e s e   + σ ·
where Q · (random energy) is the rate of heat transfer across the boundary at boundary temperature Tb. The term dS / dt represents the entropy accumulation rate within body, m · i s i input of entropy through mass intake (air, nutrients through mouth), and m · e s e is entropy leaving through exhaust of waste products (CO2, H2O through nose).
Hershey and Wang [17] used second law and homogeneous approach and estimated average specific entropy generation rate as σ · M ( t ) ( = σ · ( t ) / m B ,   W / K   per   kg   body   mass ) at given age t. The specific entropy generation over human life span (σM(t), J/kg body mass K) is given as:
σ M ( t ) [ J kg   bodymass   K ]   =   t birth t σ · M ( t ) dt
In estimating, σ · M ( t ) they assumed that all the energy released ( q · ) via oxidation cause entropy generation and found that the advection terms m · i s i m · e s e is only about 2.5% of total entropy generation rate. Their work suggests that it is not the total entropy during lifespan, but the rate of change on entropy production that define the senile death. i.e., δσ/dt (kJ/K per kg body mass) reaches zero.
Instead of entropy balance equation in estimating σ · ( t ) , Silva and Annamalai [9] adopted the availability balance equation for whole body as though it is homogeneous system:
d ( E cv T 0 S ) dt   =   n , Inlet m · n , i Ψ n , i     n , exit m · n , e Ψ n , e W · I · , E = U + KE + PE = H PV + KE + PE
where “n” represents the nutrients n (e.g., n = CH, F and P) U, internal energy, KE, kinetic energy of whole body (e.g., running) and PE, potential energy (e.g., climbing stairs), ψn,i = hn,i T0sn,i, stream availability (as per EU definition) of nutrient “n” at inlet, (e.g., availability in through intake of food intake), T0 ambient temperature), ψn,e = hn,e T0sn,e, availability exiting (e.g., through waste products n = CO2, H2O, etc.), W · all forms of work crossing the boundary, and I · irreversibility rate. It is noted that exergy (as per EU definition; but called availability in US literature. EU notation is used here) is defined as ψ n , i ψ n , i , 0 where dead state availability is defined as ψ n , i , 0 = h n , i , 0 T 0 s n , i , 0 . The oxidative process/metabolism in BS involves energy release, q · and a part of q · is used in ATP production ( W · ATP ) and remainder as heat Q · (random energy) which is energy lost during electron transport and oxidative phosphorylation. The term W · in engineering thermal systems does not contribute to entropy generation and they replaced W · in Equation (3) by W · ATP (with assumption of reversible conversion of ADP to ATP) for all internal organs. When applied to human body the BS is almost isothermal within whole body at TB and hence Ti = Te = TB, ψn = hn, T0sn = hn, TBsn = gn, Gibbs function (J/kg) of nutrient “n”. It is the “Q” which contributes a major fraction of causes entropy generation rather than q · as assumed by Hershey and Wang [17]. Ignoring KE and PE (e.g., resting). Equation (3) becomes:
dG dt = n , Inlet m · n , i g n , i n , exit m · n , e g n , e + P dV dt W · ATP I ·
where G = H − TB S and I · = T B σ · .
Ignoring P dV/dt (e.g., breathing work) and assuming quasi-steady state on daily basis,
0 = n m · n , reacted | Δ G c , n | W · ATP I · ,
Typically g n , e < g n , i and hence Δ G c , n < 0 due to irreversible oxidation process involving species “n”; but above equation uses the absolute value of Δ G c , n . By setting I · = 0 in above equation one can show that total optimum work rate is equal to sum of optimum work of each macro-nutrient
W · ATP , opt = n m · n , reacted | Δ G c , n | ,   whole   body ,   n = CH , F   and   P
where m · n , reacted , the consumption rate of nutrients which could be less than m · n , i and unreacted Gibbs function of nutrients at inlet and exit cancel out. Silva and Annamalai [8] found that the availability efficiency (ηAvail = work delivered/optimum work) used in thermodynamics for an isothermal system is same as metabolic efficiency in biology. Thus, metabolic efficiency ( η Met )
η Avail , n = η Met , n = W ATP , n W ATP , n , opt = W ATP , n | Δ G c , n |
where the ∆Gc,n is maximum possible ATP work from nutrient “n”. The irreversibility rate of a system of interest is equal to sum of irreversibility of each macro-nutrient:
I · = W · ATP , opt W · ATP = T B σ · = n m · n , reacted | Δ G c , n | ( 1 η Met , n )
Since the engineering studies for several fuels revealed that the ratio of change in Gibbs function to heating values of fuels/nutrients, Δ G ° c ,   n /HVn ≈ Cg, a constant ≈ 0.95 − 1.05 [8,16] for most nutrients. Hence
| Δ G c , n | ( 1 η Met , n ) C g HV n ( 1 η Met , n )
where HVn is the heating value of nutrient n. Thus, after dividing by mB, Equation (6) becomes
i ( kJ S kg   body   mass   K ) = I · m B = T B σ · M = n m · n , reacted m B HV n ( 1 η Met , n )
where the term m · n , reacted HV n , i is the energy release rate, SERR (see Section 2.2) and m · n , reacted HV n η Met , n is work part of energy used by body for pumping blood, work for breathing and for repair and creation of new cells thus maintaining the life functions. Aoki demonstrated that the heat part { HV n ( 1 η Met , n ) } plays major role in estimation of entropy generation of reacting system [16,18] while convective part of entropy change is negligible. The irreversibility rate “i” is a function of age “t” since energy consumption changes with age t. The inclusion of metabolic efficiency in the entropy generation concept reveals that lower ATP production (i.e., lower amount of work) should result in higher Q, and higher temperature rise around cells. The increased temperature rise (due to higher Q) can also cause increased ROS and along with decreased ATP energy available for repair and creation of new cells including the damage cause by environmental factors, and accumulation of damaged cells causing rapid aging and reduced life span. This is consistent with hypothesis of Kirkwood [6] who discounted theory of gene programming.

2.3. Heterogeneous Approach

Wang [19] adopted heterogeneous approach where the exponents in the completely body-allometric law for energy release rate (ERR) can be derived by summing of the metabolic rates of five organs: vital organs such as Brain (B), Heart (H), Kidney (K), and Liver and the rest of body (R5). This is termed as Wang-5 model. Later they included two additional organs: skeletal muscle (SM) and adipose tissue AT (which is fatty tissue packed between organs) which is termed as Wang-7 model [20] and improved the accuracy on prediction of whole body allometric law exponents. Annamalai and Silva adopted similar heterogeneous approach for entropy generation of whole body as a sum of entropy generation of the BHKLR for an isothermal body [8]. Thus, Equation (7) at organ level for k = B, H, K, L, R becomes
T B { K } σ · k { Watts K } q · k ( 1 η Met , k ) ,   k = BHKLR , η Met , k = n η Met , n q · k , n q · k ; n = CH , F   and   P
where η Met , k an energy is weighted average metabolic efficiency for whole organ k. Major nutrient (almost 100%) for brain is CH with higher metabolic efficiency and unity respiratory Quotient (RQ) while a mix of fat (about 70%) with lower metabolic efficiency and RQ = 0.7 and CH (about 20%) with higher metabolic efficiency is used in heart [21]; thus the organs may have different metabolic efficiencies. The term q · k , is energy release rate from organ k due to oxidation of all nutrients within organ k, and q · k η Met , k is energy used in ATP production by all nutrients that do not cause entropy generation just as work delivered in engineering device does not cause entropy generation. Dividing by organ mass mk, the specific entropy generation rate is obtained as:
σ · k , m { Watts K   kg   of   k } q · k , m { Watts kg   of   k } ( 1 η Met , k ) T B { K } ,   k = BHKLR ,
Since the number of cells per unit mass is approximately constant for various organs, the specific entropy stress at cell level is proportional to σ · k , m . One uses allometric laws for organs to estimate q · k , m and known η Met , k in order to estimate the entropy generation rate of organ k.

2.4. Allometric Laws-Organs

S M R k ( q · k , m )   Controversy exists on the relation between the SMRk and body mass mB. Elia’s results show that SMRk remains constant regardless of body mass (i.e., regardless of age). Wang’s data fit for six mammalian species ranging from 0.5 kg (rabbit) to 70 kg (human) yields:
SMR k = q · k , m ( W kg   of   organ   k ) = e k m B f k
Thus, for Elia, fk = 0 while for Wang fk ≠ 0. Elia’s constant SMRk data presumes that the rates are kinetically limited indicating that SMR k exp ( E R ¯ T k ) where E, activation energy (J/kmol) for oxidation, R ¯ universal gas constant and Tk, organ temperature ≈ TB, body temperature in K for a specified organ k. Since body temperature remains constant through growth of organ k, SMRk must remain constant. However the extension of recent engineering literature on oxygen-deficient combustion to BS seems to indicate partial diffusion control of O2 and hence SMRk may depend on the size of organ or organ mass (mk). These results explain why fk < 0 [14] and confirm Singer’s data [22] on SMRk indicating oxygen deficiency for large organ.
The kth organ mass is given by the allometric law:
m k = c k m B d k
where ck and dk are allometric constants at organ level. Further the previous work [8,23] converted the “body mass”-based allometric relation for each organ into “organ mass”-based allometric relation:
SMR k = q · k , m = ( e k c k f k d k ) m k f k d k + 1
The exponents in Equation (12) satisfy in inequality relation: 1 / 3 < ( f k d k + 1 ) < 0 for most organs (see areas of research 4 and 5) [14,15].
Entropy Generation Rate: Equation (10) was used in Equation (9) to obtain specific organ entropy generation rate, σ · k , m ( t ) as a function of age t; then integrating over age “t” the organ-based specific entropy generation σ k , m is given as
σ k , m ( t ) [ J kg   organ   mass   K   ] = t birth t σ · k , m ( t ) dt ,   k = B , H , K , L   ( vital   organs )   and   R
and letting t = tlife, σm,life can be obtained [16]. Using heterogeneous approach, the whole body specific entropy generation ( σ M ( t ) ) as a function of age t is given as:
σ M ( t ) [ J kg   organ   mass   K   ] = t birth t k σ · k , m ( t ) ( m k ( t ) m B ( t ) ) dt ,   k = B , H , K , L   ( vital   organs )   and   R
where mB(t) change with age t. Setting t = tlife + tbirth in Equations (13) and (14), σk,m,life and σM,life were obtained; the earlier studies by Annamalai and Silva [8] revealed that heart has highest specific entropy generation (Hm,life) while liver (Lm,life) had the lowest.
Mitochondria: The oxidation occurs within the mitochondria of the cell. Mitochondria are known to be the main source of ROS due to electron transfer process. The ROS (also called as ROM, radical oxygen metabolites) within mitochondria is believed to be primary cause of mitochondrial DNA (mtDNA) damage in Alzheimer’s disease patients via leak of electrons in the electron transport chain process [24]. If one assumes that mitochondrial densities (Mito volume % in each cell) are same in all organ, then ranking of vital organ based on entropy stress per cell must follow the same scaling as entropy stress per unit mitochondrial volume. However, allometric laws reveal that Mito volume fraction (vfMito) vary from organ to organ [25] and is given by allometric law:
vf Mito = p k m B q k
Gemme et al. [24] indicated that there is better correlation between mitochondrial metabolic rate and mitochondrial free radical production rate and longevity (e.g., birds) than with the overall metabolic rate. Indeed, the mitochondrial rate of free radical production seems to have a much stronger correlation with maximum longevity. Thus, the organ stress based on mitochondrial metabolic loading could be different compared to average metabolic loading of organ.

2.5. Biological and Chronological Aging

The second law approach relies on entropy generation as one of the biomarkers for BAR. Improved metabolic efficiency (i.e., reduced amount of heat transferred to blood or lower temperature rise in blood) slows down the aging process and thus extends life span. Section 5.4 will show that temperature rise in blood when metabolic efficiency is reduced consistent with experimental data. Section 5.5 will show slower aging with increased metabolic efficiency validating 2nd law hypothesis.
The chronological aging (CA) is measured in terms of calendar years and most retirement planning are based on CA. Using the recommended nutrition data vs. age (calendar years) provided by Food and Nutrition board (FNB), Silva and Annamalai [9] estimated the cumulative entropy generation vs. calendar years (curves AM for Male and AF for female, Figure 1) and the slope will be denoted as calendar aging rate (CAR). If consumption is more than the amount recommended for indicated age group by FNB, the SMR will be much higher and hence cumulative entropy generated is faster and thus biological aging rate (BAR) is faster than CAR. One may be at age 50 but feels like age 80 (line AB, BAR > CAR) while others may be at age 80 but feels like 50’s (line AC, BAR < CAR). Thus, Milevsky and his co-researchers recommend using BAR for retirement planning [26]. The biomarkers for BAR could be (1) energy release from unit mass of the body (qM(t)) which follows Rubner’s hypothesis (first law, homogeneous approach, ordinate in Figure 1 replaced by energy release in kJ/kg body mass); (2) energy release from unit mass of the vital organ (qk,m(t)), first law, heterogeneous approach); (3) specific entropy generation from unit mass of the body, (σM(t)) which follows Silva and Annamalai’s hypothesis (Figure 1, second law, homogeneous approach); and (4) specific entropy generation from unit mass of the vital organ (σk,m(t)), second law, heterogeneous approach). Thus, the BAR of organs differ from each other. Smaller animals undergo aging at much faster rate compared to large animals due to high Calorie consumption per unit mass of body (i.e., heart has to pump more blood per min or heart beats much faster in clock time for smaller animals). Similarly, animal spending more time outdoor ages faster than animal in door.
The biological aging rate of organs can be used to indicate how fast the organs approach chronic conditions. It has been reported that about 91% of older adults are above one chronic condition and 73% have 2 chronic conditions [27]. The top four chronic conditions that cause 2/3 of deaths are: heart disease, cancer, stroke (brain), and diabetes (kidney/liver). The organ-based second law approach predicts the top two stressed organs are heart and brain [8]. In the current work, the metabolic and entropy stress are extended to mitochondrial level to investigate whether same order is followed.

3. Objectives and Methodology

3.1. Objectives

The overall goal of current aging research is to provide biomarkers for biological aging rate (BAR) using the second law approach. The objectives of the present study are as follows: (1) Couple BAR of whole body to organ aging using heterogeneous approach; (2) Extend the 5 organ (B-H-K-L-R) analysis of BS to 7 organ analysis by including the two additional organs: adipose tissue (AT) and skeletal muscle (SM) which are of importance to athletes; (3) Extend analysis to mitochondrial level; (4) Study the effects of proportions of nutrients being oxidized, (CH:F:P) which affect metabolic efficiencies as well as maximum temperature rise of blood and which relates to entropy generation; and (5) Compare the ranking of organs based on the lifetime specific entropy generation σ m , k , life (kJ/(K kg of organ k)) vs. mitochondrial volume-based specific entropy generation σ m , k , life (kJ/(K cm3 of Mito within organ k)).

3.2. Methodology

In order to achieve the stated objectives, the following methodology is adopted:
  • Extend the 5 organ analysis (B-H-K-L-R5) to 7 organs/compartments (AT-B-H-K-L-R7-SM) and tissues (SM and AT) and residual mass (R7).
  • Calculate the specific metabolic rate, q · k . m (W/kg of organ k) and entropy generation rate within organ k, σ · k , m (W/(K kg of organ k)) (Equations (12) and (13) respectively, Section 2) as function of age(t).
  • Using allometric laws for mitochondrial density (Equation (15)), estimate Mito volume-based specific entropy generation rate per unit volume of mitochondria ( σ · v , k W/cm3 of Mito within organ k).
  • Since body mass is not constant, then use body mass growth data with age and the SMRk relation in terms of body mass (Equation (8), Section 2) to estimate growth correction factor, GCF.
  • Using GCF, rank the organs based on mitochondrial metabolic loading (MJ in organ k over life/cm3 of Mito) and entropy stress (MJ in organ k over life/(K cm3 of Mito)).

4. Modeling and Relations for Specific Energy Release and Entropy Generation of Organs

This section presents brief derivations and explicit results obtained for several parameters of interest a few of which could serve as biomarkers for biological aging. Using the literature review section (Equation (9)),
σ · k , m q · k , m ( 1 η Met , k ) T B { K } ,   k = BHKLR ,
The lifetime energy release per unit mass of organ k ( q k , m , life ) over life span is given as:
q k , m , life ( kJ kg   of   organ   mass   k ) = b birth t st , 1 q · k , m ( t ) dt + t st , 1 t st , 2 q · k , m ( t ) dt + t st , 2 t life + t birth q · k , m ( t ) dt ,   where   k = B , H , K , L , AT , SM , R 7
where q · k , m ( t ) is a function of body mass mB(t) as given by allometric law (Equation (10)). The life span tlife is typically defined from birth to death; so t life = t death t birth .

4.1. Growth Correction Factor (GCF)

The relation for mB can be found from data given in [28] and was curve fitted using
m B ( t ) = m B , st ( t t st ) z
where “t” is age in years from day of conception so that mB = 0 at t = 0, t st is age required from conception to achieve steady body mass mB,st and as such organ masses which are related to body mass remain constant after tst. It is noted that Wang et al. analyzed an existing database for organ weights of humans and found that the mass almost remains constant for most vital organs after 21 years [20]. The present analysis uses UK data for body mass vs. age, which is presented in Figure 2. The curve fit yields “z” when tst,1 = 20 years for female and tst,1 = 30 for male.
Period-I: for tbirth < t < tst,1, z = 0.75 (mass increase period), tst,1 = 20 years for female and tst,1 = 30 for male;
Period-II: for tst,1 < t < tst,2, z = 0 (steady mass period);
Period-III: for tst,2 < t < tdeath, z = −0.75 (mass decrease period).

4.2. Lifetime Specific Energy Release of Organ

Letting t* = t/tdeath where tdeath = tlife + tbirth, using allometric laws q · k , m ( t ) given in terms of body mass mB(t) (Equation (10)) and Equation (17) for variation of mB(t) with age, Equation (16) is converted into non-dimensional form for qk,m,life as given below:
q k , m , life t death e k m B , st f k = q k , m , life t death e k m B , st f k ( k J kg   of   organ   mass k ) = t birth t st , 1 t z f k d t + t st , 1 t st , 2 t z f k d t + t st , 2 1 t z f k d t
The numerator of left hand side is actual specific metabolic energy over whole life span while denominator is specific metabolic energy over whole life span in case body mass remains constant. Thus, the right hand represents the growth correction factor (GCFq) for lifetime metabolic energy since body mass is not constant. Integrating above equation
q k , m , life t life e k m B , st f k = t st , 1 ( 1 { t birth t st , 1 } g k + 1 ) g k + 1 ( t st , 2 t st , 1 ) +   t st , 2 [ | 1 t st , 2 | g k + 1 1 ] ( g k + 1 )
where gk = z fk when computing life span metabolic energy released by unit mass of organ k. Ignoring period III, and hence t st , 2 = 1 .
q k , m , life t life e k m B , st f k = GCF q ( t st , 1 , t birth , g k ) = t st , 1 ( 1 { t birth t st , 1 } g k + 1 ) g k + 1 + ( 1 t st , 1 )
where GCFq, represents growth correction factor for q (metabolic energy). When fk = 0 (Elia data), gk = 0. GCFk = 1. For GCF derivations for other parameters (e.g., lifetime energy contribution by organ k to unit body mass), the “gk” definition changes as shown in column 4 of Table 1. See Annamalai and Silva [8] and Nanda [29] for more details. The summary of algebraic relations for other biomarkers of interest for BAR is presented in Table 1.

4.3. Mitochondrial Metabolic Loading

There are a large number of mitochondria present within the cells of an organ and occupy certain fraction of volume of cells (vfMito) called mitochondrial density. Allometric law for vfMito is given as
vf Mito = p k m B q k
Nutrients and oxygen enter the cells; but the oxidation occurs within mitochondria while glycolysis occurs outside the mitochondria (Figure 3). The oxidation produces significant amount of ATP and hence ATP due to less efficient anaerobic process (glycolysis) is neglected. The organ volume is a sum of cell volume and interstitial fluid volume. The organ metabolic rate and entropy generation rate are converted to rates at mitochondrial level using allometric law for mitochondrial volume fraction:
q · k , v , Mito ( W cm 3   of   Mito   of   organ   k ) = q · k , m ( W g   of   organ   k ) ρ k ( g   of   organ   k cm 3   of   organ   k ) 1 v f cell , k ( cm 3   of   cells   in   k cm 3   of   organ   k ) 1 v f Mito ( cm 3   of   Mito   k cm 3   of   cell   k ) = ρ k e k m B f k v f cell P k m B q k = ( ρ k e k v f cell P k ) m B f k q k
where vfcell,k volume fraction of cells = nk (π dcell3/6), nk, number of cells per cm3 of organ k, dcell, cell dia. For life span metabolic energy per unit volume of mitochondria gk = z (fk − qk) in Equation (18).
Previous literature provided details on conversion of allometric laws from body mass mB to organ mass k for various parameters of interest [23]. If cell density is equal to organ density, the allometric law for metabolic rate per unit mitochondrial volume in terms of vfMito,k is given as:
q · k , Mito , vol W cm 3   Mito   organ   k = r k ( vf Mito , k ) t k , r k = e k ρ k p k   f k / q k vf cell , k , t k = f k q k q k = f k q k 1
If f k q k > 1 , then tk > 0, Mito volume-based metabolic rate (W/m3 of Mito within k) will increase with increase in vfMito,k but rate of increase will be different for different organs depending on f k q k .

4.4. Maximum Possible Blood Temperature

The ROS production rate is extremely sensitive to temperature and is function of body temperature, oxygen concentration and activation energy. A part of the energy released is converted into work (ATP) and part is released as heat Q which controls the temperature rise and is also related to entropy generation; higher is Q, higher is temperature rise ∆T (of the order of 1 °C) and higher is σ. Thus, it is of interest to estimate maximum possible blood temperature rise estimated under zero work (i.e., no ATP production). Using conservation of energy for blood entering an organ and leaving an organ, the maximum rise in blood temperature is given as:
( T exit T in ) k = Y O 2 in OEF k HV O 2 ( 1 η Met , k ) c p
where Y O 2 in , oxygen mass fraction in blood entering organ (about 300 ppm), cp, specific heat of blood, OEF, oxygen extraction fraction, HV O 2 , energy released per g of O 2 consumed and Y O 2 in OEF k HV O 2 is energy released per g blood. If V · Bl , k is blood flow rate into organ k in volume units, then term Y O 2 in OEF k HV O 2   ρ Bl   V · Bl , k ρ Bl m k is SMRk, which is given by allometric law for the organ. Thus, knowledge of OEF is not necessary.
( T exit T in ) k =   Energy   released   by   nutrients   oxidation   per   kg   blood   ×   fraction   into   heat heat   capacity   per   kg   blood
( T exit T in ) k = SMR k m k ( 1 η Met ) V · Bl , k ρ Bl c p
Using the allometric law for blood flow rates to organs V · Bl , k ( mL   to   organ   k s ) = h k m B ( kg ) i k and allometry for SMRk, the allometric relation for temperature rise is derived as
Δ T = ( T exit T in ) ,   = N k m B L k ,   N k = ( c k e k h k ρ bl c p ) , L k = f k + d k i k
Note that body mass changes with age and as such temperature rise and hence ROS production rates of organs will change with age. An approximate relation between cell death rate and/or ROS production rate, body temperature (TB) and k, bl (entropy generation within organ k, kJ/(K kg blood flow into organ) is derived in [29]; it shows a linear relation as presented below:
d [ n cell , k ] dt ,   cells   under   disfunction cm 3   s = d [ ROS ] k dt A ( 1 + E RT B Δ T k T B ) exp ( E RT B ) = A ( 1 + E R ¯ T B ( σ k , bl c p ) ) exp ( E R ¯ T B )
where E activation energy for production of ROS, R ¯ universal gas constant (kJ/K kmol), ( σ k , bl c p ) = SMR k m k ( 1 η Met ) V · Bl , k ρ Bl T B c p , and A, pre-exponential factor; higher is heat part of energy released (which is higher when metabolic efficiency is low e.g., brain after concussion [30]), higher is temperature rise and higher is production rate of ROS.

5. Results and Discussion

Data input for thermal and chemical properties of nutrient are presented in Table 2 followed by basic allometric constants for Wang-5 organ and Wang-7 organ models where 5 and 7 denote number of organs considered in the model (Table 3). The derived allometric constants and blood temperature rise are summarized in Table 4. The results on entropy stress at organ and mitochondrial levels are then presented and compared. Finally, brief procedures are presented for use by medical professional in qualitative evaluation of the biological aging based on Calories/energy released.

5.1. Nutrient Properties and Body Mass Growth Data

The thermal and chemical properties of CH, F and P are presented in Table 2. It is seen that the HV O 2 is almost constant and ∆G°/∆H° is close to unity for most nutrients. Here mst = 84 kg, z = 0.75, tst,1 = 24 years and tst,2 = tdeath = 75 years for quantitative estimates (Equation (17)).

5.2. Basic Allometric Constants

The following basic allometric constants were used: (A) organ mass vs. body mass (mB) based on Wang-5 model and body mass-dependent Wang-5 SMRK data, (B) organ mass vs. body mass (mB) for Wang-7 organ with Elia’s body mass independent constant SMRk; and (C) blood flow rate, and mitochondrial volume vs. body mass.
  • Wang-5-Wang SMRk or SMRk (5 Organ Model): Here the whole body is divided into 5 organs: B, H, K, L, and R5 where R denotes residual mass and R5 denotes residual mass in the 5-organ model. Table 3 lists the basic allometric constants. The in vivo measurements of O2 concentration differences between artery and vein, and blood flow rate measurement provided data on organ metabolic rates. Using data for 6 species (Rat, rabbit, cat, dog, human-1, human-2) (0.48–70 kg), Wang et al. presented allometric relation for body mass-dependent SMRk and organ mass in terms of body mass and the basic allometric constants are listed (within parentheses) in Table 3 [19].
  • Wang 7–Elia-SMRk Model: The allometric constants were presented by Wang et al. considering a body with 14 organs: four higher metabolically active vital organs: BHKL, two intermediate organs, AT and SM and eight low level metabolic organs with equal metabolic intensity: adrenal, blood, gut, lung, skin, skeleton, spleen and thyroid. For the current work, the whole body is divided into 7 organs: BHKL, AT, SM and R7 (R7: Residual mass for 7-organ model). The allometric constants for SMR of R7 are different from R5 due to the following: body mass independent SMRk throughout growth [13] i.e., constant irrespective of size/mass of selected vital organ, R7 mass is different and the R7 includes 8 sub-organs: adrenal, blood, gut, lung, skin, skeleton, spleen and thyroid. The allometric constant for R7 are curve fitted in terms of body knowing the constants for 8 sub-organs. Table 3 lists the basic allometric constants for Wang-7 organ model.
  • Blood Flow Rate and Mitochondrial Volume Fraction (MiV): The basic allometric constants for blood flow rate to each organ (mL/s per kg organ k), and average mitochondrial volume fraction (vfk,Mito, cm3 of Mito/cm3 of cells within organ k), etc are given in terms of body mass (mB, kg) (Table 3).

5.3. Derived Allometric Constants

The body mass (mB, kg)-based allometric relations for (1) SMRk or q · k , m , q · k , V , Mito (W/cm3 of Mito for organ k); (2) blood flow rate into organ k (mL/s per kg organ k); (3) Mitochondrial volume fraction (vfk,Mito, cm3 of Mito/cm3 of cells within organ k); (4) oxygen extraction fraction OEFk (= oxygen used for energy release/oxygen supplied) in each organ k; and (5) maximum temperature rise (T − TB, °C) for blood flow around an organ k can be re-derived in terms of organ mass (mk, kg) and/or body mass mB and/or vfk,Mito. Derivations for the derived allometric relations are omitted here; for derivation details refer to thesis of Nanda [29]. Table 4 summarizes the derived allometric relations and lists the constants.

5.4. Maximum Possible Blood Temperature

Sacher et al. assumed that life span is a function of 4 variables: body mass (mB) and temperature (TB), specific metabolic rate of body and brain mass [36]. Since metabolic rate is related to mB, and body temperature (TB) is affected by the blood temperature, life span relation could be reduced to two variables: brain and body mass [37] if TB is constant. Since brain mass is related to body mass [19,30], then life span may be further reduced to a single variable of either body mass or metabolic rate. It is noted that these expressions for life span do not seem to account for varying body mass with age or growth correction factor (GCF) while Annamalai and Silva [8] accounted for GCF and as well as the effect of metabolic efficiency (which indirectly affects body temperature by affecting “Q” and hence Sacher’s relation showing effect of body temperature). Whichever organ’s venous blood temperature is high, it faces more oxidative damage. Using the relations presented in Section 4 and derived allometric constants listed in Table 4, the maximum temperature rise of vital organs were computed for 84 kg human (mB,st, metabolic efficiency = 0) and they are tabulated in Table 5. Except for brain data for Wang 7-Elia model, the temperature rise in almost all organs is less than 0.5 °C. Even though brain constitutes only 2–3% body mass, it consumes 20–25% of body’s total oxygen or total energy release indicating high specific energy release rate (W/kg brain mass) of body. Brain temperature in humans is normally 0.5–1 °C warmer compared to temperature around spinal cord of body [38]. Note that actual rise in temperature is about 70% of maximum possible rise when metabolic efficiency of about 30%. A concussion in brain causes metabolic efficiency to be about 15% (estimated from ATP data of [30] which shows ATP decreasing from 2320 nmol/g to 1095 nmol/g wet within 4 days after concussion) and according to current analysis the increased heat part of energy release will cause increased brain temperature with brain injury (i.e., due to reduced metabolic efficiency) is confirmed with data of increase in brain temperature by 1–2 °C [38,39]. Further change in temperature by more than 1 °C affects neurological functions of rat °C with metabolic rate increasing by 7% per °C rise in temperature. ROS relations in combustion science show that the ROS production rates are exponential functions of temperature (Equation (24)) and proportional to the concentrations of fuel and oxygen thus the decreased metabolic efficiency results in more ROS damage. The following affects the temperature rise:
  • body mass, which affects the blood flow rate and metabolic rates of an organ
  • blood flow rates and OEF during resting period exercise (during exercise the metabolic rate may increase to 480 W, almost 6 times that of normal resting metabolic rate which requires skin blood flow rate of almost 1.9 L/min for 1.8 m2 area person. Core to skin temperature difference is of the order of 5 °C with skin temperature around 30–31 °C) [40]
  • metabolic efficiency
  • nutrients oxidized and concussion in brain which affect metabolic efficiency
Lk (exponent for temperature rise allometry, Equation (23)) ≈ 0 for brain and kidney and hence temperature rise is fixed regardless of age. For heart Lk > 0 (Table 3), indicating increase in temperature rise with body mass which is an indication of increased ROS production while for liver Lk < 0 indicating decrease in temperature rise with body mass. If metabolic efficiency is reduced then there is further temperature rise.

5.5. Results Based on Organ

Values for cumulative entropy generation (σk,ife, MJ/K per kg of k) for each organ over life span were generated for both Wang-5 and Wang-7 models.

5.5.1. Metabolic Efficiency

In estimating the overall metabolic efficiency of all organs, the % energy released by each nutrient is used. With assumption of 55:30:15, CH:F:P (mass %), the energy release % is computed as 36.1%, 49.5% and 14.4% via CH, F and P [9] and hence η Met =   0.361 × 38.2 +   0.495 × 32.2 + 0.144 × 0.104   =   0.31 . Salin et al. [41] reports that there is variation of P/O ratio within an individual and amongst individuals due to diet, temperature and environmental factors. The current work presumes that it is affected only by diet, i.e., composition % of CH:F:P metabolized.

5.5.2. Growth Correction Factor (GCF)

With increase in body mass, the organ masses change along with SMRk (Wang-5 model) where fk is non-zero (Equation (10)); the exponent is typically negative and varies depending on size of an organ (theoretical rationale given in Miller [15]) and fk = 0 when Elia model is used. Recall that body mass is a function of age and metabolic rate increases with increase in body mass and reaches a steady value when body mass reaches a steady value. Equation (18) presents a relation for the growth correction factor GCF. A chart has been constructed for GCF vs. non-dimensional age with gk (= fk z) as parameter. Figure 4 shows the results for growth correction factor, GCF; GCF = 1 if mass of body remains constant (z = 0) irrespective of age or when fk = 0 (Elia Constant).
Last Column in Table 1 illustrates the procedure with an example for estimating lifetime entropy generation for the human heart.

5.5.3. Effects of Diet Composition

Figure 5 shows the results for organ entropy stress for the Wang-7 model as a function of various nutrient % metabolized in organs. The allometry-based EXCEL program for 7 organs was checked by summing up lifetime entropy contribution by each organ to unit body mass and the allometry-based sum was estimated to be 6050 kJ/(K kg body mass), which is slightly less than the estimated value by Silva and Annamalai (Figure 1) with use of dietary intake recommended by Food and Nutrition board. It can be inferred from Figure 5 that the lower metabolic efficiency of 55:30:15 diet (ηMet = 0.31) compared to pure CH (ηMet = 0.38) results in the higher values for Q, entropy generation rate (or ROS rate), and higher lifetime entropy generation for prescribed life span. Almost pure glucose (>90%) is used in brain with high ATP production rate accompanied by lesser heat release ( Q · ) and hence lesser cell temperature rise and lesser entropy generation rate per kg of tissue. The ATP or higher metabolic efficiency was shown to slow the process of aging [38], which seems to suggest the ability of ATP to restore/repair the damaged cell or create new cells with full functional characteristics. Increased ATP production or increased metabolic efficiency [35] (i.e., reduced entropy generation rate, slower BAR Figure 1) when flies are subjected to near Infra-red light and life span increased from 3 weeks (controls) life span increased 11–12 weeks; further activities increased. Edman, et al. [36] has shown that lifespan improves with a low-protein diet since low protein concentration implies higher CH and improving metabolic efficiency.
Higher metabolic efficiency (ηM) represents better use of nutrients in maintaining the healthy state of organs within BS. Further in support of the current hypothesis, pigeons were found to live ten times longer (35 years) than rats (4 years) [42], even though both species have the same mass and hence similar metabolic rates. It is noted that the effects of nutrients on aging and life span are studied strictly from the point of view calories released and metabolic efficiency. Typically, glucose concentration is 1000 ppm (or g per million g of blood in engineering terminology) while oxygen is only 300 ppm indicating rich nutrient:oxygen mixture. However the presence of increased % of glucose in blood stream affects the filter system (kidneys) and its effects on entropy stress of kidney and other organ functions are not studied. While the oxygen extraction fraction (OEF) is only 0.10–0.15 for healthy kidneys (0.5% of body mass), the diabetes requires additional metabolic input (more O2 consumption resulting in more ROS) in order to perform more filtration of solute [43] thus accelerating the entropy production rate or aging for kidneys and these factors are not included in the model.

5.5.4. Effects of 5 Organ vs. 7 Organ Models

Figure 6 shows a comparison of heart normalized entropy stress for the two models: Wang-5 and Wang-7 organ Model with selected diet (CH:F:P = 55:30:15). The Wang-5 include BHKL-R5 along with body mass-dependent SMRk while the Wang-7 Model includes 7 organs AT-BHKL-R7-SM with Elia’s body mass independent SMRk. Note that Wang-7 with Elia data shows almost same entropy stress as Wang-5 model.

5.6. Results Based on Mitochondria

5.6.1. SMRk (W/kg Organ k) Variation with Mito Volume

Typically, metabolism is modeled using Michaelis and Menten (MM) kinetics, which are analogous to Langmuir Hinshelwood (LH) kinetics on carbon combustion except for the difference that MM is based on volume of reacting system while LH is based on surface kinetics. If temperature of body is fixed, the kinetics controlled rate in W/cm3 of Mitochondria is a function of temperature and hence it is typically assumed that the consumption rate of O2 per unit volume of organ remains constant in most organs as presumed in Krogh cylinder model for metabolism. However the oxidation occurs in mitochondrion of the cell and as such the O2 consumption rate per unit volume of Mitochondrion is constant and higher the MiV fraction within the cell, higher is the SMRk. Thus, a plot on the SMRk (energy release rate per unit mass of organ) vs. mitochondrial volume % is presented in Figure 7 for both Wang-5 and Wang-7 models. Brain has the lowest Mito volume % (5.4%) while heart has the highest Mito volume % (15.5%). It is seen that higher the mitochondrial volume % in organ, higher the specific metabolic rate organ (W per kg of organ). If SMRk (W/kg organ k) is proportional to (MiV%)n, then least square straight line fit for vital organs yields the slopes as n = 0.83 (i.e., n < 1) with Wang-5 and 0.67 for Wang7-Elia but both indicates positive slopes. Note that Elia’s SMRk used in Wang-7 model is body mass independent but depends on type of organ, its enzyme and mitochondrial volume fraction. It is apparent from both models that SMRk does not increase in direct proportion to MiV fraction probably due to organ-dependent enzyme activities and/or oxygen transport rate limitation [14,15], i.e., slopes are less than 1. As MiV % increased consumption rate of O2 increases (Watts) which require faster diffusion rate from capillaries to the cells. If one compares these results to carbon combustion in engineering literature, diffusion control is a strong possibility particularly when mitochondrial volume is increased.

5.6.2. SMRk,V,Mito (W/cm3 of Mito) Variation with Mito Volume

In order to investigate whether SMR k , V , Mito (W/cm3) is same for all organs, and to ascertain the degree of damage to DNA within mitochondria, which is proportional to metabolic loading, SMR k , V , Mito quantitative estimation of SMR k , V , Mito is required for all organs. If mitochondrial volume fraction, MiVfk, number of cells per unit volume (nk), and density of organ k are known, then one can convert SMRk to SMRk,V,Mito (Watts/cm3 of Mito of k). Whether one use 7 organ or 5 organ data, the Mito volume % are same since it is a function of body mass only. While the plots in Figure 7 show that SMRk increases with the mitochondrial volume % in organ, the plots in Figure 8 shows that SMRK,v,Mito decreases with increase in MiV %. Since SMRk in Figure 7 are proportional to (MiVf)n where n < 1 and since SMR k , V , Mito W cm 3   of   Mito   in   k is proportional to (MiVf)n/MiVf ∝ (MiVf)n−1, then (n − 1) < 0 (See Figure 7); then n − 1= −0.17 for Wang 5 and −0.33 for Wang-7 while power law fit in Figure 8 indicates the corresponding slopes as −0.153 and −0.315.The R2 is very low since SMRk was divided by mitochondrial volume and uncertainty in Mito volume % is magnified.
It has been shown in earlier literature that SMRk is low particularly for large organs due to problems with oxygen accessibility [22], since capillaries must service a large number of cells in a large organ. In engineering literature detailed models called group combustion have been developed on the effect of oxygen deficiency on specific energy release rate (W/kg) vs. radius of spherical fuel clouds [12]; increasing size or mass of cloud (which is similar to increasing size of organ in biology) results in reduced energy release rate per unit mass of cloud. When this theory is extended to organs it shows decreasing SMRk with increasing size of organ [14,15]. To summarize, the SMRk is not only affected by mitochondrion volume %, and enzyme reactivity within mitochondrion of cell but also by transport rate of O2 and oxygen accessibility to all cells.
In order to ascertain the degree of damage to organ k relative to heart, the heart normalized metabolic intensity of mitochondrial intensity ratio (MIR) is defined as
MIR = q · Mito , k ( t ) q · Mito , H ( t ) = ρ k ρ H e k p H p k m B Q H Q k ,   relative   Mito DNA   damage
which assumes that all organs have same number density of cells and same cell size (Vcell,k = Vcell,H). Figure 9 shows the results. The net energy release rate per mitochondrial volume of the organs differ due to the difference in mitochondrial volume density (MiV) for each organ.
From Figure 10 it can be inferred that life span entropy generated per unit volume of mitochondrion is highest in the brain, followed by the heart, liver, and kidney. This result is vital for football athletes, who are subject to multiple concussions to the brain, which further reduce the metabolic efficiency of the organ [30,38], thereby increasing blood temperature and hence entropy generation i.e., brain is biologically aging faster for the athlete with repeated concussion. Once again, it is important to note that the ranking of the organs based on entropy generation changes when the analysis is based on MiV %. As the mitochondria is considered the engine of the cell, this model is more accurate to quantify organ failure as compared to whole organ model.

5.6.3. SMRk,V,Mito Variation of Organs with Body Mass

The Figure 7, Figure 8 and Figure 9 present results for SMRk,V,Mito of 84 kg human. Using Equation (19) and data in Table 3 for vital organs, SMRk,V,Mito can be obtained as function of body mass. Table 3 indicates that (fk − qk) < 0 indicating increased SMRk,V,Mito for brain and liver of smaller species compared to humans and vice versa, and (fk − qk) ≈ 0 for kidney and heart, indicating almost constant SMRk,V,Mito irrespective of body mass of BS.

6. Suggested Procedure for Estimation of Biological Aging

How these results could be used in the Medical Field in estimating metabolism-based biological aging rate (BAR)? Based on energy release rate per unit mass, there are two possible methodology for determining the biological aging rate (BAR) and entropy generation rate:
Method I: Organ Allometry Method: The determination of BARk of vital organ k (W/kg organ k) and the total BARM of organs per kg body mass (W/kg body mass) vs. age require non-intrusive measurements of oxygen consumption rate of each organ and mass of each organ. In absence of non-intrusive measurements, the body mass (mB(t)) vs. age(t) and the allometric laws for each organ in terms of body mass outlined here can be used to estimate metabolic rate of each organ (t) (W per kg organ mass) and the whole body q · M ( t ) (W per kg mass of body), vs. age may be used as biomarkers. With assumption of constant metabolic efficiency, σ · M ( t ) (W/(K per kg mass of organ)) and σ · M ( t ) (W/(K per kg mass of body)) vs. age(t) can also be determined.
Method II: Overall Body Calorimetry Method: Overall specific “Calories” release rate q · M ( t ) and entropy generation rate, σ · M ( t ) vs. age can be estimated from measured body mass and wet nasal gas analyses assuming dominant nutrients metabolized are CH and fat.
Steps in Organ Allometry (Method-I, heterogeneous): The procedure is given below for monitoring the biomarkers: q · k , m ( t ) and σ · k , m ( t ) (W/ (K kg body mass)} and σ M ( t ) (kJ/ (K kg body mass)) vs. age(t). Here step by step procedure is given with an example for heart of a human with steady mass of 84 kg (mB,st).
A. Energy Release Rates
  • First assume as though the person has constant body mass or steady mass (for e.g., the UK data reveals mB,st = 84 kg). Estimate kth organ mass using allometric relation and known ck, and dk (Table 3); thus for heart, k = H, mH,st = 0.48 kg when mB,st = 84 kg.
  • Estimate specific energy release rate (SERRk) or specific metabolic rate (SMRk) (W/kg organ k) during steady body mass period with known ek and fk given by allometric law:
    q · k , m , st ( W kg   organ   k ) = e k m B , st f k , q · H , m , st = 25.3 W kg   heart ,   ( k = H )
  • Multiplying tlife one get GJ per kg of heart over life span
    q · k , m , st = q · k , m , st t life   where   t life   in   s = 75 years × 3.15 × 10 7 s year = 2.37 × 10 9 s
    q H , m , st , life = 25.3 W kg   heart × 2.37 × 10 9 s × 1   GJ 10 9   J = 60   J kg   heart
    See last column in Table 4. Similarly estimate q · k , m , st and q k , m , st , life for all organs.
  • Estimate the contribution rate by each organ to unit mass of body:
    q · k , M , st W   by   k kg   body = q · k , m , st m k , st m B , st
    So with k = H, the contribution rate by heart to unit body mass is estimated as
    q · k , M , st = 25.3   W × 0.48   kg   heart 84   kg = 0.145   W
    q ˙ k , M , s t = 25.3   W ×   0.48   k g   h e a r t 84   k g = 0.145   W   c o n t r i u t i o n   r a t e   b y   h e a r t   t o   u n i t   m a s s   o f   b o d y
    Multiplying tlife one get GJ contributed by heart over life span to each unit body mass
    q H , M , st = 0.145 W   by   heart kg   body   mass × 2.37 × 10 9   s 10 9   s GJ = 0.35 GJ   by   heart kg   body
    Similarly estimate q · k , M , st and qk,M,st,life for all organs.
  • Add contributions by all organs to each unit body mass under steady period and check the total with estimate from Kleiber’s law. In Method A, the biomarkers which determine BAR are empirical organ allometric coefficients, steady mass and the metabolic rate (First law, heterogeneous) which determines the slope of biological aging curve in Figure 1.
B. Entropy Generation Rates
6.
Assuming constant average metabolic efficiency ( η Met ) for all organs, estimate heat part of SERRk
Q · k , m , st = q · k , m , st ( 1 η Met ) , Q · H , m , st = 25.3 W kg   heart ( 1 0.31 ) = 17.5 W kg   heart   with   η Met = 0.31
and estimate corresponding entropy generation rate from organ k as
σ · k , m , st = Q · k , m , st T B ,   T B   in   K = 310   K   and   for   heart
σ · H , m , st = 17.5 310 = 0.0565 W kg   heart   K
7.
Multiplying tlife (= 2.37 × 10 9 s over 75 years) one get GJ of heat per kg of heart over life span as
Q k , m , s t , l i f e = q ˙ k , m , s t ( 1 η M e t ) , Q H , m , s t , l i f e = 17.5 W k g   h e a r t × 2.37 × 10 9 s × 1   G J 10 9   J = 41.5   G J k g   h e a r t
and corresponding
σ H , m , s t = σ ˙ H , m , s t t l i f e = 41.5   G J 310   K   k g   h e a r t × 1000   M J G   J = 134   M J k g   h e a r t   K
8.
Following part A, the contribution by each organ to unit body mass are given below:
Q ˙ k , M , s t = q ˙ k , m , s t   ( 1 η M e t )   m k m B , s t , Q ˙ H , M , s t = 0.1 W   b y   h e a r t k g   b o d y , Q H , M , s t , l i f e = 0.24 G J   b y   h e a r t k g   b o d y
σ H , m , s t , l i f e = 0.78 M J   b y   h e a r t K   k g   b o d y   m a s s
9.
Add entropy contributions by all organs to each unit body mass under steady period to determine lifetime entropy generation per unit body mass. In Method B, the biomarkers which determine BAR are empirical organ allometric coefficients, steady mass, entropy generation rate (second law, heterogeneous)which determines the slope of biological aging curve in Figure 1.
10.
If Elia’s model is used where specific organ metabolic rate does not change with body mass (e.g., from birth to death), all the numbers in rate form per unit mass of organ will not change or growth correction factor = 1 for all rate form based on unit mass of organ. However, for Wang-5, the SMRk change with age since body mass does not remains constant from birth to death and as such GCF are necessary. Further organ mass changes with age (due to change with body mass with age) and hence their contribution to unit body mass will change even for Elia model; as such estimates for GCF are necessary when contribution to unit body mass is required for both Elia and Wang-5 models. Thus, one must use charts presented in Figure 4. When period of standard mass is short compared to life span, GCF ≈ 1.0.
Method II: Overall Body Calorimetry (Second Law, Homogeneous)
In Method II, the biomarkers that determine BAR are body mass, breathing rate and nasal exhaust analyses, which determine the whole body specific energy release rate (SERRM) and metabolic efficiency as function of age in determining the biological aging curve in Figure 1. In addition, one may use empirical organ allometry and sum over all organs to estimate SMRM and check with measured SMRM.
  • Monitor the patient’s breathing rate V ˙ B r e a t h e (L/s), mass of body (mB(t)) and CO2, O2 and H2O % in nasal exhaust gases vs. age.
  • If metabolism is dominated by CH (nutrient 1) and F (nutrient 2), then using nasal exhaust analysis, get RQ using Figure 11 [45]:
    R Q = ( X CO 2 , e ( N 2 / O 2 ) i X N 2 , e X O 2 , e ( N 2 / O 2 ) i )
    where Xk, mole (or volume) fraction of gaseous species k in dry exhaust.
  • Determine heat fraction contributed by CH and fat. Since RQ of blended nutrient mixture is given as
    R Q = H F 1 R Q 1 + ( 1 H F 1 ) R Q 2
    where HF1, ratio of energy released by nutrient 1 (CH) to energy released by the blend (mixture of CH and fat). Then for RQ1 ≠ RQ2, the above equation yields heat fraction faction of CH, (nutr 1) and fat (nutr 2)
    H F 1 = (   R Q   R Q 2 ) {   R Q 1 R Q 2 } ,    R Q 2 R Q 1
    H F 2 = (   R Q   R Q 1 ) {   R Q 1 R Q 2 } ,    R Q 2 R Q 1
    where HF2 represent the ratio of energy released by metabolism of fat (i.e., “Calories” released by fat) to total energy released by blend. RQ1 for glucose = 1, RQ2 for fat is 0.7. The mole fraction of CH in nutrient mixture can also be estimated if needed.
    X C H = M o l e F r C H = ( 23   R Q 16 ) ( 17   R Q 10 )
    Knowing molecular weights, mass fraction of CH can be determined.
  • Estimate oxygen extraction fraction (OEF) from inspired air:
    O E F = 1 X O 2 , e X N 2 , e ( N 2 / O 2 ) i ,
  • Estimate the overall metabolic efficiency of the body vs. age:
    η Met =   { Mole   Fr   CH   η Met , CH   Δ G CH + ( 1 Mole   Fr   CH ) η Met , F Δ G F } { M o l e   F r   C H   Δ G CH + ( 1 M o l e   Fr   CH ) Δ G F } η Met   H F 1 η Met , CH + ( 1 H F 2 ) η Met , F
    Figure 12 shows variation of metabolic efficiency and glucose mass % metabolized with RQ if it is assumed that dominant energy release occurs via metabolism of a mixture of CH and F. If RQ increases with age, the more glucose metabolized and less fat metabolized indicating accumulation of “fat” mass if “fat” consumption is known.
  • Using nasal exhaust analysis, get oxygen extraction fraction (OEF), determine specific metabolic rate (SMRM)
    S M R M = S M t R   { W a t t s } m B { k g } = V ˙ B r e a t h e ( L s )   ρ a i r ( k g L ) 0.23 O E F H V O 2 ( J k g   O 2 ) m B
    where oxygen mass fraction in air is assumed to be 0.23. Further specific entropy generation rate vs. age is given as
    σ ˙ M ( t ) ( k W K   k g   ) =   S M R M ( t ) ( 1 η M e t ( t ) ) T B { K } ,
    As such BAR cane be compared with CAR as shown in Figure 1. Equation (27) accounts for variation of metabolic rate and metabolic efficiency with age while the results reported in Figure 9 and Figure 10 are based on constant metabolic efficiency.
For e.g., if measured CO2, H2O and O2 % are 3.6, 6.2 and 15.6 [16], then using atom balance, RQ = 0.85, OEF = 0.21, ηMet = 0.352, CH % (mole): 80%; CH % (mass): 74%; rest is fat. ηMet = 0.37; If breathing rate is 360 LPH, then metabolic rate is 83 W and one can check whether the SMR of 83 W is close to value based on Kleiber’s law.
Note on Biological Aging and Lean Body Mass (LBM)
An average male has 76–82% of LBM while female has 69–75%.The LBM % decreases with age [46]. It is LBM which results metabolism. Thus when fat % is increased, the decreased mass of LBM must increase its specific metabolic rate so that heat loss from fat is compensated and fat mass can be maintained at body temperature. Hence the previous numbers under methods I and II must be divided by fraction of LBM in the body in order to estimate BAR. Hence biological aging rate is accelerated when LBM % decreases. Thus an additional measurement during office visit could be LBM %. Recently CDC reported that people of age 65 and older must exercise at least 2.5 h of exercise per week and it improves lean body mass % and mitochondrial functions [47]. Further exercise results in increasing open % of capillaries, serves more number cells and increasing blood flow under same pressure difference across the organs; it results in more ATP per unit volume of organ with more ability to repair and create new cells.

7. Summary and Conclusions

  • While Rubner’s relies on 1st law analysis for estimating the aging rate (irrespective of metabolic efficiency), the current work uses 2nd law analysis for estimating entropy generation per unit mass of organ and per unit volume of mitochondria, which can be used as one of the biomarkers for predicting the biological aging rate.
  • The inclusion of entropy generation for adipose tissue and skeletal muscle has minimal effect on whole body specific entropy generation rate.
  • Increased % of glucose diet seems to increase life span while increased protein % shortens life span due to low metabolic efficiency of proteins.
  • The increase in SMRV (per unit volume of organ) is less than proportional to increase in mitochondrial density (Mitochondrial volume %, MiV) within the cell.
  • The heart is the most stressed organ whether one adopts 5- or 7-organ model for resting humans when organ-based entropy stress model is used.
  • Ranking of entropy stress based on unit volume of mitochondria differs from entropy stress ranking based on unit mass of organ. At a mitochondrial level, the brain has the highest stress and liver has the lowest loading among all vital organs. Particularly athletes with constant concussions may be subject to more brain anomalies near the end of their lifetimes due to low metabolic efficiency with concussion.
  • While ATP produced via oxidation of nutrients supplies energy for repair and replacement of cells (synthesis) through oxidation, such an oxidation path also serves as the “destroyer” (decomposition) of cells through generation of ROS and hence involves delicate balance.
  • With currently available diagnostic tools for measurements of nasal exhaust gas composition, respiration rate and body mass vs. age, the whole body specific energy release and entropy generation (which accounts for whole body average metabolic efficiency) methodologies can be used to provide one of the biomarkers for tracking biological aging.

Acknowledgments

The authors wish to thank (1) three anonymous reviewers for bringing attention to review articles in this area, a few more pertinent references and many useful discussions which have been incorporated in the manuscript; and (2) Jothi Sukkawala MD for bringing attention to recent aging research of Mayo Clinic [47]. The principal author close to retirement recently started translating results from combustion and thermodynamics to biology and hence wish to acknowledge the funding provided through Paul Pepper Professorship and the James J. Cain Chair in Mechanical Engineering at Texas A&M University, which partially supported the graduate student(s) and the open access fees for the publication of this manuscript.

Author Contributions

Kalyan Annamalai conceived most of the basic concepts, formulated the governing equations, obtained the solutions for several parameters of interest, developed the Excel-based programs for use in obtaining quantitative data, prepared the rebuttal to three reviewers and finally revised the paper with supporting data from literature. Arnab Nanda collected literature and ran the various parametric cases as part of manuscript and wrote the initial draft of paper.

Conflicts of Interest

The authors declare no conflict of interest.

Symbols and Acronyms

ATPAdenosine Tri-Phosphate
ATAdipose Tissue
BBrain
BABiological Aging
BARBiological Aging Rate
BSBiological Systems
CSpecific Heat Capacity
CAChronological Aging
CARChronological Aging Rate
CHCarbohydrates
CRCalorie Restriction
CVControl Volume
DNADeoxyribonucleic Acid
ERREnergy Release Rate
EUEuropean Union
FFat
FNBFood and Nutrition Board
GCFGrowth Correction Factor
HHeart
HVHeating Value
HHVHigher Heating Value
KKidney
KEKinetic Energy
LLiver
LBMLean Body Mass
LHLangmuir Hinshelwood
MiVMitochondrial Volume Density
MMMichaelis Menten
OEFOxygen Extraction Factor
pPressure
PProtein
PEPotential Energy
RRest of Organs
ROLRate of Living theory
ROSRadical Oxygen Species
SMSkeletal Muscles
UInternal Energy
VVolume
vfMitoVolume Fraction of Mitochondria
WWork

Nomenclature

hEnthalpy, kJ/kg
IIrreversibility, kJ
I ˙ Irreversibility rate, kJ/s
MMass, kg
mBBody Mass
mkMass of organ, k
m ˙ nMass flow rate of nutrient n in organ k
m ˙ O 2 , n , k Consumption rate of oxygen by nutrient n in organ k
QHeat, kJ
Q ˙ Heat transfer rate due to metabolic heat release q k ˙ at organ k, kJ/s
q ˙ k,nSpecific metabolic energy release rate from organ k per unit mass of organ k
q ˙ k,MEnergy release rate of organ k contributed to the unit mass of body
SEntropy, kJ/K
sSpecific Entropy, kJ/kg K
TBBody temperature, K
tTime or age
tstTime to reach steady weight
UInternal Energy
WKWork delivered by metabolism at organ k
Δ G c ° Gibbs Free Energy for Combustion
Δ G M ° Gibbs Free Energy for metabolism (with ATP production)

Greek Symbols

ηMetabolic efficiency
σEntropy generation, kJ/K
σM,kEntropy contribution to unit mass of body by whole organ k
σ ˙ MEntropy generation rate per unit body mass (W/kg of body mass K)
σ ˙ m,kSpecific entropy generation rate of organ k (W/K-kg of k)
ψStream availability, kJ/kg
νO2,nStoichiometric oxygen mass per unit mass of nutrient n
ηn,kMetabolic efficiency of nutrient n in organ k

General Notes

A bar (-) on top of any property indicates its specific property per kmole of substance
A dot (.) on top of any property indicates its time rate of change

Appendix A. Nasal Gas Analyses

The readers are referred to more general derivations provided in [45] involving any C–H–O nutrients/fuels. Here derivations are provided for specific nutrients: glucose (CH) and fat (F). Consider a mixture of x moles of carbohydrate (CH) represented by C6H12O6 and (1 − x) moles of fat represented by C16H32O2. When this mixture undergoes stoichiometric oxidation,
x   C 6 H 12 O 6 + ( 1 x )   C 6 H 32 O 6 + ( 23 17 x )   O 2 ( 16 10 x )   CO 2 + ( 16 10 x )   H 2 O
The RQ factor along with “x” for a fuel can be determined from known CO2 % and O2 % in nasal exhaust gas. The general methodology is to formulate the following reaction equation. Assuming complete combustion, using atom conservation for CH: F fuel mixture;
x   C 6 H 12 O 6 + ( 1 x )   C 16 H 32 O 2 + a   O 2 + b   N 2 c   CO 2 + d   H 2 O + e   O 2 + f   N 2
Equation (A2) shows that there are 7 unknowns for C–H–O fuel: x, a, b, c, d, e, and f. Thus 7 equations are needed to determine the unknowns. Four equations are obtained using atom balance of C, H, O, and N. The three additional equations are generated as follows. (1) b/a in dry air = 79/21 = 3.76; (2) dry CO2 mole fraction in products, XCO2 = c/Ndry where Ndry = c + e + f; and (3) dry O2 mole fraction in products, XO2 = e/Ndry. Thus all 7 unknowns can be solved including “x”. In addition A: F ratio (= (b + a)), RQ (= (16 – 10x)/(23 − 17x)), and, equivalence ratio, φ (= O2 used or required in combustion/O2 supplied, same as oxygen extraction fraction, OEF for lean mixture or φ < 1). Derivations are skipped (See [39] for more details and any general C–H–O fuel). Charts are presented in the main body for limited parameters.

References

  1. Flatt, T. A New Definition of Aging? Front. Genet. 2012, 3, 148. [Google Scholar] [CrossRef] [PubMed]
  2. Weinert, B.T.; Timiras, P.S. Physiology of Aging Invited Review: Theories of aging. J. Appl. Physiol. 2003, 95, 1706–1716. [Google Scholar] [CrossRef] [PubMed]
  3. Mostafavi, H.; Berisa, T.; Day, F.R.; Perry, J.R.B.; Przeworski, M.; Pickrell, J.K. Identifying genetic variants that affect viability in large cohorts. PLoS Biol. 2017, 15, 9. [Google Scholar] [CrossRef] [PubMed]
  4. Medvedev, Z.A. An attempt at a rational classification of theories of ageing. J. Biol. Rev. Camb. Philos. Soc. 1990, 65, 375–398. [Google Scholar] [CrossRef]
  5. Lee, D.; Hwang, W.; Artan, M.; Jeong, D.-E.; Lee, S.-J. Effects of nutritional components on aging. Aging Cell 2015, 14, 8–16. [Google Scholar] [CrossRef] [PubMed]
  6. Kirkwood, T. Understanding the Odd Science of Aging. Cell 2005, 120, 437–447. [Google Scholar] [CrossRef] [PubMed]
  7. Guachalla, L.M.; Ju, Z.; Koziel, R.; von Figura, G.; Song, Z.; Fusser, M.; Epe, B.; Jansen-Durr, P.; Rudolph, K.L. Sod2 haploinsufficiency does not accelerate aging of telomere dysfunctional mice. Aging 2009, 1, 303–315. [Google Scholar] [CrossRef] [PubMed]
  8. Silva, C.; Annamalai, K. Entropy Generation and Human Aging: Lifespan Entropy and Effect of Physical Activity Level. J. Entropy 2008, 10, 100–123. [Google Scholar] [CrossRef]
  9. Silva, C.; Annamalai, K. Entropy generation and human aging: Lifespan entropy effect of diet composition and caloric restriction diets. J. Thermodyn. 2009, 2009, 1–10. [Google Scholar] [CrossRef]
  10. Andresen, B.; Shiner, J.S.; Uehlinger, D.E. Allometric scaling and maximum efficiency in physiological eigen time. PNAS 2002, 99, 5822–5825. [Google Scholar] [CrossRef] [PubMed]
  11. Pinker, S. The Second Law of Thermodynamics, Life History: Why Things Fall Apart. The Wall Street Journal, 31 December 2016–1 January 2017.
  12. Annamalai, K.; Puri, I.K. Combustion Science and Engineering; CRC Press, Taylor and Francis: Boca Raton, FL, USA, 2006; ISBN 9780849320712. [Google Scholar]
  13. Elia, M. Organ and tissue contribution to metabolic rate. In Energy Metabolism: Tissue Determinants and Cellular Corollaries; Kinney, J.M., Tucker, H.N., Eds.; Raven Press, Ltd.: New York, NY, USA, 1992; pp. 61–79. [Google Scholar]
  14. Annamalai, K.; Miller, J.A. Link between O2 Deficient Metabolism in Organs and Group Combustion in Engineering. In Proceedings of the 10th U.S. National Combustion Meeting, College Park, MD, USA, 23–26 April 2017. [Google Scholar]
  15. Miller, J. Oxygen Deficient Metabolism in Organs: A link to Combustion Science. Master’s Thesis, Texas A&M University, College Station, TX, USA, December 2014. [Google Scholar]
  16. Annamalai, K.; Silva, C. Thermal Sciences and Biological Systems, Chapter 14. In Advanced Thermodynamics Engineering, 2nd ed.; Annamalai, K., Puri, I., Jog, M., Eds.; CRC Press: Boca Raton, FL, USA, 2011; p. 1096. ISBN 9781439805725. [Google Scholar]
  17. Hershey, D.; Wang, H. A New Age-Scale for Humans; Lexington Books: New York, NY, USA, 1980. [Google Scholar]
  18. Aoki, I. Entropy production in human lifespan: A thermodynamical measure for aging. AGE 1994, 1, 29–31. [Google Scholar] [CrossRef]
  19. Wang, Z.; Zhang, J.; Ying, Z.; Heymsfield, S. Organ-Tissue Level Model of Resting Energy Expenditure across Mammals: New Insights into Kleiber’s law. ISRN Zool. 2012, 2012, 1–9. [Google Scholar] [CrossRef]
  20. Wang, Z.; Ying, Z.; Bosy-Westphal, A.; Zhang, J.; Schautz, B.; Later, W. Specific metabolic rates of major organs and tissues across adulthood: Evaluation by mechanistic model of resting energy expenditure. Am. J. Clin. Nutr. 2010, 92, 1369–1377. [Google Scholar] [CrossRef] [PubMed]
  21. Kummitha, C.M.; Kalhan, S.C.; Saidel, G.M.; Lai, N. Relating tissue/organ energy expenditure to metabolic fluxes in mouse and human: Experimental data integrated with mathematical modeling. Physiol. Rep. 2014, 2, e12159. [Google Scholar] [CrossRef] [PubMed]
  22. Singer, D. Size Relationship of Metabolic Rate: Oxygen Availability as the Missing Link between Structure and Function. Thermochim. Acta 2006, 446, 20–28. [Google Scholar] [CrossRef]
  23. Annamalai, K.; Silva, C. Entropy Stress and Scaling of Vital Organs over Life Span Based on Allometric Laws. J. Entropy 2012, 14, 2550–2577. [Google Scholar] [CrossRef]
  24. Gemma, C.; Vila, J.; Bachstetter, A.; Bickford, P. Oxidative Stress and the Aging Brain: From Theory to Prevention. Chapter 15. In Brain Aging: Models, Methods, and Mechanisms; CRC Press/Taylor & Francis: Boca Raton, FL, USA, 2007. [Google Scholar]
  25. Else, P.L.; Hulbert, A.J. Mammals: An allometric study of metabolism at tissue and mitochondrial level. Am. Physciol. Soc. 1985, 248, R415–R421. [Google Scholar]
  26. Huaxiong, H.; Milevsky, M.A.; Salisbury, T.S. Retirement Spending and Biological Age. Available online: https://ssrn.com/abstract=2918055 (accessed on 22 October 2017).
  27. Making Sure Our Longer Lives Are Healthy Ones. Available online: https://www.forbes.com/sites/nextavenue/2017/02/06/making-sure-our-longer-lives-are-healthy-ones/#5e2a662e28f6 (accessed on 22 October 2017).
  28. Summary Report of a Peer Involvement Workshop on the Development of an Exposure Factors Handbook for the Aging. Available online: https://cfpub.epa.gov/ncea/efp/recordisplay.cfm?deid=171923 (accessed on 22 October 2017).
  29. Nanda, N. Entropy Stress Based on Organ and Mitochondrial Metabolic Loading. Master’s Thesis, Mechanical Engineering, Texas A&M University, College Station, TX, USA, May 2017. [Google Scholar]
  30. Vagnozzi, R.; Tavazzi, B.; Signoretti, S.; Amorini, A.M.; Belli, A. Temporal window of metabolic brain vulnerability to concussions: Mitochondrial-related impairment—Part I. Nuerosurgery 2007, 61, 379–388. [Google Scholar] [CrossRef] [PubMed]
  31. Lee, C.P.; Gu, Q.; Xiong, Y.; Mitchell, R.A.; Ernster, L. P/O ratios reassessed: mitochondrial P/O ratios consistently exceed 1.5 with succinate and 2.5 with NAD-linked substrates. FASEB J. 1996, 10, 345–350. [Google Scholar] [PubMed]
  32. Skeletal muscle mass and distribution in 468 men and women aged 18–88 yr. Available online: http://jap.physiology.org/content/89/1/81.short (accessed on 22 October 2017).
  33. Wenger, B.C. The Regulation of Body Temperature. Chapter 29. Temperature Regulation and Exercise Physiology. Available online: http://www.academia.edu/4940904/The_Regulation_of_Body_Temperature (accessed on 22 October 2017).
  34. Hoppeler, H.; Lindstedt, S.L.; Claassen, H.; Taylor, C.R.; Mathieu, O.; Weibel, E.R. Scaling mitochondrial volume in heart to body mass. Respir. Physiol. 1984, 55, 131–137. [Google Scholar] [CrossRef]
  35. Abrahams, S.; Greenwald, L.; Stetson, D.L. Contribution of renal medullary mitochondrial density to urinary concentrating ability in mammals. Am. J. Physiol. 1991, 261, R719–R726. [Google Scholar] [PubMed]
  36. Porter, R.K.; Brand, M.D. Causes of differences in respiration rate of hepatocytes from mammals of different body mass. Am. J. Physiol. 1995, 269, R1213–R1224. [Google Scholar] [PubMed]
  37. Davies, M. On Body Size and Tissue Respiration. J. Cell. Compar. Physiol. 2005, 57, 135–147. [Google Scholar] [CrossRef]
  38. Sacher, G. Longevity, Aging, and Death: An Evolutionary Perspective. Gerontologist 1978, 18, 112–120. [Google Scholar] [CrossRef] [PubMed]
  39. Calder, W. Scaling of Physiological Processes in Homeothermic Animal. Ann. Rev. Physiol. 1981, 43, 301–322. [Google Scholar] [CrossRef] [PubMed]
  40. Segolene, M.; Fanny, V.; Thomas, G. Brain Temperature: Physiology and Pathophysiology after Brain Injury Review Article. Anesthesiol. Res. Pract. 2012, 2012, 1–13. [Google Scholar]
  41. Brain Core Temperature of Patients with Mild Traumatic Brain Injury as Assessed by DWI-Thermometry. Available online: https://link.springer.com/article/10.1007/s00234-014-1384-5 (accessed on 22 October 2017).
  42. Variation in the Link between Oxygen Consumption and ATP Production, and Its Relevance for Animal Performance. Available online: http://rspb.royalsocietypublishing.org/content/282/1812/20151028 (accessed on 22 October 2017).
  43. Flyod, R.A.; West, M.; Hensley, K. Oxidative Biochemical Markers Cluse to Understanding Aging in Long Lived Species. Exp. Gerontol. 2001, 36, 619–640. [Google Scholar] [CrossRef]
  44. Hansell, P.; William, J.W.; Roland, C.B.; Fredrik, P. Determinants of kidney oxygen consumption and their relationship to tissue oxygen tension in diabetes and hypertension. Clin. Exp. Pharm. Physiol. 2013, 40, 123–137. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, Z.; Heshka, S.; Gallaghar, D.; Boozer, C.; Kotler, D.P.; Heymsfield, S.B. Resting energy expenditure-fat-free mass relationship: New insights provided by body composition modeling. Am. J. Physiol. Endocrinol. Metab. 2000, 279, E539–E545. [Google Scholar] [PubMed]
  46. Annamalai, K.; Siva, T.; Ranjan, D. Ranking Renewable and Fossil Fuels on Global Warming Potential Using Respiratory Quotient (RQ) Concept. J. Combust. 2017, in press. [Google Scholar]
  47. Enhanced Protein Translation Underlies Improved Metabolic and Physical Adaptations to Different Exercise Training Modes in Young and Old Humans. Available online: http://www.cell.com/cell-metabolism/fulltext/S1550-4131(17)30099-2?_returnURL=http%3A%2F%2Flinkinghub.elsevier.com%2Fretrieve%2Fpii%2FS1550413117300992%3Fshowall%3Dtrue (accessed on 22 October 2017).
Figure 1. Entropy generation vs. aging in calendar tears (Chronological aging, CA based on recommended diet by Food and Nutrition Board, [8] for male and female. Curve AF: Female, Curve AM: Male. Within the upper limit of about 11,000 kJ per kg body mass, life span of man and woman are 70 (line AM) and 78 years (line AF). The slopes of AM and AF represent chronological aging rates (CAR) of male and female. With higher SERRM (kJ/kg body mass; uncontrolled diet, e.g., line AB), the biological aging rate, BAR > CAR (e.g., Reduced metabolic efficiency); under diet control (e.g., line AC, BAR < CAR [8]; (e.g., Increased metabolic efficiency). The BA lines are qualitative. Note that the BAR of cats, dogs, ants etc. are much faster compared to BAR of humans due to high SMR of whole body or high surface area to volume ratio. For Rubner’s hypothesis based on first law, the ordinate could be replaced by specific metabolic energy; lines AB and AC refer to energy (calories) consumption per unit mass of body with line AC referring to “Calorific restriction”, and AB refers to higher “calorific consumption”. Further line AB may also correspond to energy released from smaller species indicating faster BAR for smaller species. Under similar metabolic efficiency throughout life span both the 1st and 2nd law hypotheses yield similar results.
Figure 1. Entropy generation vs. aging in calendar tears (Chronological aging, CA based on recommended diet by Food and Nutrition Board, [8] for male and female. Curve AF: Female, Curve AM: Male. Within the upper limit of about 11,000 kJ per kg body mass, life span of man and woman are 70 (line AM) and 78 years (line AF). The slopes of AM and AF represent chronological aging rates (CAR) of male and female. With higher SERRM (kJ/kg body mass; uncontrolled diet, e.g., line AB), the biological aging rate, BAR > CAR (e.g., Reduced metabolic efficiency); under diet control (e.g., line AC, BAR < CAR [8]; (e.g., Increased metabolic efficiency). The BA lines are qualitative. Note that the BAR of cats, dogs, ants etc. are much faster compared to BAR of humans due to high SMR of whole body or high surface area to volume ratio. For Rubner’s hypothesis based on first law, the ordinate could be replaced by specific metabolic energy; lines AB and AC refer to energy (calories) consumption per unit mass of body with line AC referring to “Calorific restriction”, and AB refers to higher “calorific consumption”. Further line AB may also correspond to energy released from smaller species indicating faster BAR for smaller species. Under similar metabolic efficiency throughout life span both the 1st and 2nd law hypotheses yield similar results.
Entropy 19 00566 g001
Figure 2. Population-Weighted Mean Weight for NHANES3 Subjects of Different Ages (2.6 months to 96 years) (adopted from [28] and modified). The body mass, mB for humans is not constant and grows with age (t, since conception, Period I) until reaching steady mass mB,st at t = tst,1, remains almost constant until t = tst,2 (Period II) along with the rate of entropy production fixed and then it decreases tst,2 < t < tdeath (III). Current analysis ignores period III. Curve fit for growth for age 0 to 75 years: 0 < t < tst,1: m B = m B , st ( t / t st ) z , Male: mB,st: 84 kg, tst,1 = 24 years, tst,2 =75 years, z = 0.75 [8]; female: mB,st = 70 kg, tst,1 = 20 years, tst,2 = 75 years, z = 0.75.
Figure 2. Population-Weighted Mean Weight for NHANES3 Subjects of Different Ages (2.6 months to 96 years) (adopted from [28] and modified). The body mass, mB for humans is not constant and grows with age (t, since conception, Period I) until reaching steady mass mB,st at t = tst,1, remains almost constant until t = tst,2 (Period II) along with the rate of entropy production fixed and then it decreases tst,2 < t < tdeath (III). Current analysis ignores period III. Curve fit for growth for age 0 to 75 years: 0 < t < tst,1: m B = m B , st ( t / t st ) z , Male: mB,st: 84 kg, tst,1 = 24 years, tst,2 =75 years, z = 0.75 [8]; female: mB,st = 70 kg, tst,1 = 20 years, tst,2 = 75 years, z = 0.75.
Entropy 19 00566 g002
Figure 3. Schematic diagram of the thermodynamic system for estimation of metabolic loading based on mitochondrial volume; a large number of mitochondria within cell of given size results in increased mitochondrial volume (MiV) %.
Figure 3. Schematic diagram of the thermodynamic system for estimation of metabolic loading based on mitochondrial volume; a large number of mitochondria within cell of given size results in increased mitochondrial volume (MiV) %.
Entropy 19 00566 g003
Figure 4. Effect of Body mass growth on growth correction factor GCFk. For period I (see Figure 2), mB(t) = mB,st (t/tst)z and period II, mB(t) = mB,st a, z = 0, gk = 0; (1) For lifetime specific organ energy release or entropy generated: gk = z fk specific organ allometry (non-isometric, fk = 0) on the “F” factor. gk = z fk; z > 0 indicates body mass growth and z < 0 indicates decrease in body mass (period II); fk = 0 for constant SMRK (e.g., Elia data); Wang-5 data indicates 0 < fk < −0.27; in authors’ opinion 0 < fk < −1/3. Thus, mostly gk < 0 during body mass growth period and gk > 0 during body mass decrease period for most organs since fk < 0. For period I, mB(t)= mB,st (t/tst)z and period II, mB(t) = mB,st a, z = 0, gk = 0.; (2) gk = z (fk + dk – 1) in chart for lifetime energy release and lifetime entropy contribution by organ k to each unit mass of body; gk = dk − 1 in case fk = 0; (3) gk = z (fk + dk) for lifetime energy release and lifetime entropy contribution by organ k to whole body; gk = dk in case fk = 0 (Elia constants); (4) gk = z (fk − qk) for lifetime energy release and lifetime entropy contribution by organ k per cm3 of Mito.
Figure 4. Effect of Body mass growth on growth correction factor GCFk. For period I (see Figure 2), mB(t) = mB,st (t/tst)z and period II, mB(t) = mB,st a, z = 0, gk = 0; (1) For lifetime specific organ energy release or entropy generated: gk = z fk specific organ allometry (non-isometric, fk = 0) on the “F” factor. gk = z fk; z > 0 indicates body mass growth and z < 0 indicates decrease in body mass (period II); fk = 0 for constant SMRK (e.g., Elia data); Wang-5 data indicates 0 < fk < −0.27; in authors’ opinion 0 < fk < −1/3. Thus, mostly gk < 0 during body mass growth period and gk > 0 during body mass decrease period for most organs since fk < 0. For period I, mB(t)= mB,st (t/tst)z and period II, mB(t) = mB,st a, z = 0, gk = 0.; (2) gk = z (fk + dk – 1) in chart for lifetime energy release and lifetime entropy contribution by organ k to each unit mass of body; gk = dk − 1 in case fk = 0; (3) gk = z (fk + dk) for lifetime energy release and lifetime entropy contribution by organ k to whole body; gk = dk in case fk = 0 (Elia constants); (4) gk = z (fk − qk) for lifetime energy release and lifetime entropy contribution by organ k per cm3 of Mito.
Entropy 19 00566 g004
Figure 5. Lifespan Entropy Stress of organ MJ/K per kg organ mass—Wang-7 organ (1) CH:F:P = 55:30:15 Mass % (35:47:18, Heat %), RQ = 0.82, Met Eff = 0.312; (2) CH:F = 55:45 Mass % (33:67, Heat %), RQ = 0.79, Met Eff = 0.34; (3) CH = 100 Mass % (100% Heat %), RQ = 1.0, Met Eff = 0.38; Entropy stress is lower with pure CH since less Q is liberated; all organs are assumed to have same Met Eff.
Figure 5. Lifespan Entropy Stress of organ MJ/K per kg organ mass—Wang-7 organ (1) CH:F:P = 55:30:15 Mass % (35:47:18, Heat %), RQ = 0.82, Met Eff = 0.312; (2) CH:F = 55:45 Mass % (33:67, Heat %), RQ = 0.79, Met Eff = 0.34; (3) CH = 100 Mass % (100% Heat %), RQ = 1.0, Met Eff = 0.38; Entropy stress is lower with pure CH since less Q is liberated; all organs are assumed to have same Met Eff.
Entropy 19 00566 g005
Figure 6. Comparison between Wang-5 and Wang-7 Models for heart normalized Entropy stress. The 7 organ data had high brain mass (2.4 kg) vs. brain mass of 0.32 kg for 5 organ data. CH:F:P = 55:30:15, RQ = 0.82, Met Eff = 0.31, mst = 84 kg. Lifetime specific entropy generation for heart is 135 MJ (K kg heart) for Wang-5 and 110.62 MJ/(K kg heart) for Wang-7 essentially due to similar constant metabolic rate for Elia mode for both Kidney and heart. The heart is the most entropy stressed organ for Wang-5 and both heart and kidney are most stressed for Wang-7 model while AT and SM are least stressed. Thus, cells with higher metabolic rate per unit mass (e.g., heart) leads to faster biological aging of organs. Vital organs account for 4–6% body mass but account for metabolic rate 15–40 times that of SM and 5–100 times that of AT [44] indicating that vital organs age much faster compared to other organs.
Figure 6. Comparison between Wang-5 and Wang-7 Models for heart normalized Entropy stress. The 7 organ data had high brain mass (2.4 kg) vs. brain mass of 0.32 kg for 5 organ data. CH:F:P = 55:30:15, RQ = 0.82, Met Eff = 0.31, mst = 84 kg. Lifetime specific entropy generation for heart is 135 MJ (K kg heart) for Wang-5 and 110.62 MJ/(K kg heart) for Wang-7 essentially due to similar constant metabolic rate for Elia mode for both Kidney and heart. The heart is the most entropy stressed organ for Wang-5 and both heart and kidney are most stressed for Wang-7 model while AT and SM are least stressed. Thus, cells with higher metabolic rate per unit mass (e.g., heart) leads to faster biological aging of organs. Vital organs account for 4–6% body mass but account for metabolic rate 15–40 times that of SM and 5–100 times that of AT [44] indicating that vital organs age much faster compared to other organs.
Entropy 19 00566 g006
Figure 7. Relation between mitochondrial volume % ( MiV 100 ) and specific metabolic release rate (W/kg) for humans-84 kg person; simplified analysis and kinetics controlled metabolic rate indicates linear dependency of SMRk on Mitochondrial volume %. SMRk (W/kg organ k) is proportional to (MiV %)n with n = 0.83 for Wang-5 and 0.67 for Wang-7. Note n < 1.
Figure 7. Relation between mitochondrial volume % ( MiV 100 ) and specific metabolic release rate (W/kg) for humans-84 kg person; simplified analysis and kinetics controlled metabolic rate indicates linear dependency of SMRk on Mitochondrial volume %. SMRk (W/kg organ k) is proportional to (MiV %)n with n = 0.83 for Wang-5 and 0.67 for Wang-7. Note n < 1.
Entropy 19 00566 g007
Figure 8. The Volume-based Specific Metabolic Rate of Organ k (W/cm3 of Mito) vs. Mito volume %. It decreases with increased mitochondrion volume %, possibly due to diffusion limitations. SMRk,V (W/cm3 of Mito) is approximately proportional to (MiVf)n−1 (See Figure 7); then n − 1 = −0.17 for Wang 5 and −0.33 for Wang-7 while Figure indicates corresponding slopes as −0.153 and −0.315.The R2 is very low since SMRK was divided by mitochondrial volume and uncertainty in Mito volume % is magnified.
Figure 8. The Volume-based Specific Metabolic Rate of Organ k (W/cm3 of Mito) vs. Mito volume %. It decreases with increased mitochondrion volume %, possibly due to diffusion limitations. SMRk,V (W/cm3 of Mito) is approximately proportional to (MiVf)n−1 (See Figure 7); then n − 1 = −0.17 for Wang 5 and −0.33 for Wang-7 while Figure indicates corresponding slopes as −0.153 and −0.315.The R2 is very low since SMRK was divided by mitochondrial volume and uncertainty in Mito volume % is magnified.
Entropy 19 00566 g008
Figure 9. Comparison of Heart normalized Metabolic Stress based on unit organ mass (W/kg) and unit Mitochondrial volume (W/cm3 of Mito), 84 kg person; 1.26 W/cm3 of Mito (or 1260 W/kg of Mito) in heart, and 25.3 W/kg of heart. Wang7-Elia Model.
Figure 9. Comparison of Heart normalized Metabolic Stress based on unit organ mass (W/kg) and unit Mitochondrial volume (W/cm3 of Mito), 84 kg person; 1.26 W/cm3 of Mito (or 1260 W/kg of Mito) in heart, and 25.3 W/kg of heart. Wang7-Elia Model.
Entropy 19 00566 g009
Figure 10. Lifespan Entropy Loading in Mitochondria given in MJ life of Mito in k/(K cm3 of Mito); Wang7-Elia Model.
Figure 10. Lifespan Entropy Loading in Mitochondria given in MJ life of Mito in k/(K cm3 of Mito); Wang7-Elia Model.
Entropy 19 00566 g010
Figure 11. Determination of RQ (primary Y axis) and Oxygen Extraction Fraction (OEF, secondary Y (same as equivalence ratio) using nasal dry exhaust gas analysis. Valid for any C–H–O fuel under complete combustion (products CO2, H2O, O2 and N2) or metabolism of nutrients. Solid lines: RQ vs. XCO2 with oxygen mole fraction as a parameter; dashed lines: OEF vs. XCO2 with oxygen mole fraction as a parameter. As CO2 % in exhaust increase, it implies increasing production of CO2 or RQ for given O2 %. At given CO2 %, increase in O2 % in exhaust implies decreasing amount of O2 consumed or increasing RQ. The equivalence ratio for lean mixtures in engineering (O2 required to completely burn given amount of fuel/O2 supplied for same amount of fuel) is same as ratio of oxygen used for metabolism to oxygen supplied or oxygen extraction fraction in biology. The lean and rich mixture concept for fuel: air mixtures in engineering system does not apply to typical rich mixtures of nutrients (CH and F) and oxygen in blood stream (Figure adopted from [45]).
Figure 11. Determination of RQ (primary Y axis) and Oxygen Extraction Fraction (OEF, secondary Y (same as equivalence ratio) using nasal dry exhaust gas analysis. Valid for any C–H–O fuel under complete combustion (products CO2, H2O, O2 and N2) or metabolism of nutrients. Solid lines: RQ vs. XCO2 with oxygen mole fraction as a parameter; dashed lines: OEF vs. XCO2 with oxygen mole fraction as a parameter. As CO2 % in exhaust increase, it implies increasing production of CO2 or RQ for given O2 %. At given CO2 %, increase in O2 % in exhaust implies decreasing amount of O2 consumed or increasing RQ. The equivalence ratio for lean mixtures in engineering (O2 required to completely burn given amount of fuel/O2 supplied for same amount of fuel) is same as ratio of oxygen used for metabolism to oxygen supplied or oxygen extraction fraction in biology. The lean and rich mixture concept for fuel: air mixtures in engineering system does not apply to typical rich mixtures of nutrients (CH and F) and oxygen in blood stream (Figure adopted from [45]).
Entropy 19 00566 g011
Figure 12. Variation of Metabolic Efficiency and glucose mass % metabolized for a mixture of CH and F vs. RQ. CH: Glucose, F: Fat. Solid Line: Met Eff vs. RQ. Dashed Line: glucose mass % metabolized vs. RQ.
Figure 12. Variation of Metabolic Efficiency and glucose mass % metabolized for a mixture of CH and F vs. RQ. CH: Glucose, F: Fat. Solid Line: Met Eff vs. RQ. Dashed Line: glucose mass % metabolized vs. RQ.
Entropy 19 00566 g012
Table 1. Summary of Results Lifetime Entropy Generation of Whole body (Annamalai & Silva [23]).
Table 1. Summary of Results Lifetime Entropy Generation of Whole body (Annamalai & Silva [23]).
#Desired VariableNon-Dimensional Form of Desired VariableGrowth Correction Factor, GCF Right Side of Equation (18); Note 2Example for Heart, (Wang 5 Data SERRK Body Mass Depend; mBst = 84 kg, 75 years, tst,1 = 24 years, RQmix = 0.80, ηmet = 0.312, mH,std = 0.006 × 840.98 = 0.46 kg, z = 0.75; tst2 = tlife = 75 years (Note 1))
1Lifetime Specific metabolic energy release by organ k ( kJ k g   o f   o r g a n   m a s s   k ) q k , m , l i f e t d e a t h   q ˙ k , m , s t   or   q k , m , l i f e t l i f e   e k   m B , s t f k gk in Figure 4 = z fk; gk = 0 with Elia dataz = 0.75, fH = −0.12; gH = zfH = 0.75 × (−0.12) = −0.09, GCF from chart = 1.015 (Wang SERRk), = 1 for Elia; q ˙ H.m,st = 21 W/kg heart; qH,m,Life,st = 60.GJ/kg heart, qH,m,Life = 60 × 1.015 = 60.8 GJ/kg heart
2Lifetime Specific Entropy Generation of organ k, ( k J K   k g   o f   o r g a n   m a s s   k ) T σ k , m , l i f e t d e a t h   q ˙ k , m , s t ( 1 η )   or   T B σ k , m , l i f e t d e a t h   e k   m B , s t   f k ( 1 η k ) gk in Figure 4= z fkGCF = 1.015, t l i f e   e k   m B , s t f k ( 1 η k ) / T B , σH,m,st = 134 MJ/(K kg heart) σH,m,s = 134 × 1.015= 135 MJ/(K kg heart)
3Lifetime metabolic energy contribution by organ k to unit mass of body, ( k J   b y   k   k g   o f   b o d y   m a s s   ) q k , M , l i f e   t d e a t h   q ˙ k , M , s t   or   q k , m , l i f e ( m k , s t q k , m , s t , l i f e / m B , s t ) gk in Figure 4 = z(fk + dk − 1)gH = z(fH + dH − 1)=0.75 × (−0.12 + 0.98 − 1) = −0.105, GCF = 1.02 q H , M , Life , st = t l i f e   q ˙ k , M , s t = 60 × 0.46/84 = 0.33 GJ by heart to 1 kg body mass, qH,M,life = 0.33 × 1.02 = 0.336 GJ to 1 kg body mass
4Lifetime Entropy Generation contribution by organ k to unit mass of body ( k J   b y   k K   k g   o f   b o d y m a s s ) T σ k , M , l i f e   t d e a t h   q ˙ k , M , s t ( 1 η )     or   T σ k , m , l i f e t d e a t h   c k e k   m B , s t f k + d k 1 ( 1 η k ) gk in Figure 4 = z(fk + dk − 1)GCF = 1.02, σH,m,Life,st = 0.33 × 1000 × (1 − 0.31)/310= 0.73 MJ/K by heart to 1 kg body mass; σH,m,Life = 0.73 × 1.02 = 0.745 MJ/K by heart to 1 kg body mass
5Lifetime metabolic energy contribution by organ k to whole body (qk,life) q k , l i f e t d e a t h   q ˙ k , m , s t   or   q k , m , l i f e   ( m k , s t q k , m , s t , l i f e / m B , s t ) gk in Figure 4 = z(fk + dk)gH = z(fk + dk) = 0.75 × (−0.12 + 0.98) = 0.645, GCF = 0.874, qH,Life,st = 0.33 × 84 = 27.72 GJ st by heart to whole body qH,Life, = 27.72 GJ × 0.874 = 24.1 GJ by heart to whole body
6Lifetime Entropy Generation contribution by organ k to whole body T σ k , l i f e t d e a t h   q ˙ k , s t ( 1 η )   or   T σ k , m , l i f e t d e a t h   c k e k   m B , s t f k + d k ( 1 η k ) gk in Figure 4 = z(fk + dk)σH,Life,st = 27.72*(1 − 0.31) × 1000/310 = 61.34 MJ/K by heart to whole body at std mass; σH,Life,st = 61 × 34 × 0.874 = 53.6 MJ/K to whole body over life span
7Lifetime Specific metabolic energy release per m3 mito within organ k ( k J m 3   o f   M i t o   w i t h i n   o r g a n   k ) q k , V , M i t o , l i f e t d e a t h   q ˙ k , V , s t , M i t o   or   q k , V , M i t o , l i f e t d e a t h   ( 1 v f c e l l )   ( ρ k q ˙ k , m , s t p k m B , s t q k ) gk in Figure 4 = z (fk − qk)fH = −0.12 + 0.044 = −0.076; GCF = 1.015
Note 1: The Nallometric coefficients presume no AT, no SM; coefficients are from [23]. Note 2: When fk = 0, F = 1; one obtains Elia’s results; q ˙ k , m , s t   or SBMRk are: 0.581, 11.622, 21.307, 21.307, 9.685, 0.630, 0.581 W/kg of k, where F is called growth correction factor.
Table 2. Properties of Nutrients. RQ mix = ( ν CO 2 , CH X CH + ν CO 2 , F X F + ν CO 2 , P X CH , P ) / ( ν O 2 , CH X CH + ν O 2 , F X F + ν O 2 , P X CH , P ).
Table 2. Properties of Nutrients. RQ mix = ( ν CO 2 , CH X CH + ν CO 2 , F X F + ν CO 2 , P X CH , P ) / ( ν O 2 , CH X CH + ν O 2 , F X F + ν O 2 , P X CH , P ).
NutrientsM, kg/kmolSt. O2 kg/kg Nutrient (Mole per Mole Nutrient)HHV MJ/kg (MJ/kmol) = |ΔHc°|RQ H H V O 2 M J/kg of O2 Δ G ¯ ° c MJ/kmol Δ G ¯ ° M MJ/kmol | Δ G ¯ ° M | | Δ H ¯ ° c | kmol O2 Consumed per kmol ATPMetabolic Efficiency (%)
Glucose(CH), C6H12O61801.066 (6 mole O2/mole CH)15.630 (2813 MJ/kmol)114.65−2895-1.030.15838.2
Fat (F), C16H32O22562.869 (23 moles O2/mol F)39.125 (10,015 MJ/kmol)0.713.60−9840, −9800−67150.980.21732.2
Protein(P), C4.57H9.03N1.27O2.25S0.0461191.53 (5.70 moles O2/mole empirical P)22.79 (2712 MJ/kmol)0.80 0.8214.87−2665-0.98-10.4
Mixture, CH:F:P, 55:30:15 (Mass %)1821.68 (9.56 moles O2 per mole mix with empirical P)23.75 (4328 MJ/kmol)0.8214.1464,319,894-0.999-31.23
Note 1: These are energy ratios. Typically P/O ratio = 2.5–3.0 for NAD linked substrate and 1.5–2.0 for succinates [31].
Table 3. Basic Allometric Constants-Wang-7 Elia (Wang-5 models in parenthesis except as indicated). q ˙ k ,   m ( W k g   of   organ   k ) =   e k   m B f k , m k   ( k g   o r g a n   m a s s )   =   c k   m B d k , V ˙ B l , k ( mL   to   organ   k s ) =    h k   m B ( kg ) i k , v f M i t o = p k m B q k . For q ˙ k,m,st for 85 kg human; pk qk values from Else, P L and Hulbert A J for mammals ranging from 18–2067 g. Fat Free mass = adipose tissue free mass = body mass − adipose mass [32].
Table 3. Basic Allometric Constants-Wang-7 Elia (Wang-5 models in parenthesis except as indicated). q ˙ k ,   m ( W k g   of   organ   k ) =   e k   m B f k , m k   ( k g   o r g a n   m a s s )   =   c k   m B d k , V ˙ B l , k ( mL   to   organ   k s ) =    h k   m B ( kg ) i k , v f M i t o = p k m B q k . For q ˙ k,m,st for 85 kg human; pk qk values from Else, P L and Hulbert A J for mammals ranging from 18–2067 g. Fat Free mass = adipose tissue free mass = body mass − adipose mass [32].
OrganDensity g/ccck, kg [33]dkek W/kgfkhk (Note 4)ik (Note 4)pk (Note 3)qk (Note 3) q ˙ k , m , s t ( W k g   o f   k )
-ρk m k ( kg ) = c k m B ( t ) d k ,   m B   in   kg q ˙ k , m ( t ) ( W k g   o f   o r g a n   k ) =   e k   m B   f k V ˙ B l , k ( mL   to   organ   k s ) =   h k   m B ( k g )   i k ( MiVf ) = p k   m B q k -
Adipose Tisue-7 organ [20]0.9–0.92 [20]0.0753 (-)1.19 (-)0.220.0--- (-)- (-)-
Brain1.036 [20]0.1025 (0.011)0.71 (0.76)11.62 (21.62)0.0 (−0.14) −0.07 [34]0.1010.7040.0538 (0.0538) [25]0.00 (0.00) [25]0.0442
Heart (1.06) [20]1.060.006 (0.006)0.98 (0.98)21.3 (43.113)0.0 (−0.12)0.2240.7140.282 [34] (0.227)−0.044 [34] (−0.09)0.146
Kidney (1.05) [20]1.050.0089 (0.007)0.71 (0.85)21.3 (33.414)0.0 (−0.08)0.6030.7740.562 [35] (0.144)−0.056 [35] (−0.14)0.108
Liver (1.06) [20]1.060.0491 (0.0330)0.70 (0.87)9.7 (33.113)0.0 (−0.27) −0.24 [25]0.1560.8560.18 [36] (0.112)−0.09 [36] (−0.13)0.186
Residual Mass (RM)-Wang 2000-7 organ mass- Elia-Mod (SMRK for AT, BHKL SM and R with Wang)-0.3296 (0.939)1.01 (1.01)0.58 (1.446)0.0 (−0.17)--- (-)- (-)(0.192)
Skeletal Muscle [20,25]1.040.4683 (-)0.99 (-)0.63 (-)0.0 (-)0.769 (0.769)0.737 (0.737)0.0438 (0.0438) [25]−0.09 (−0.09) [25]-
Note 1: Elia values are: BHLK, SM, AT, RM: 11.62, 21.3, 9.7, 21.3, 0.63, 0.22, 0.58 W/kg [20]. Adipose Tissue; 21% body mass ref male (70 kg); 33% in ref female; wt: 12% of body wt; mB = 63 kg, adipose: 63 × 0.12 = 7.56 kg; 0.453 L of O2/h or 0.453 × 14,500 = 2.6 W or 2.6/7.56 = 0.33 W/kg corrected 0.165; SM: 0.492 × 0.63 kg; 3 L of O2/h = 0.43 W/kg [25,37]. Note 2: the allometric constants for mitochondrial volume fraction were selected from [34,35,36] for heart, liver and kidney respectively, for BS mass ranging from 0.02 kg to 250 kg. For kidney, the kidney cells are assumed to follow same allometric relation as reported for modular TAL. However for brain mitochondrial volume fraction for BS with wide ranging mass are not available. Hence Else and Hulbert data were selected for MiV generated for mass ranging from 0.02 to 2 kg. Since exponent for qk is zero (Brain), indicating body mass is independent of MiVf, due to which it can be applied to BS with higher mass. Note 3: Values in (…) parentheses are based on Else Hubert [25] data for MiVf allometry for mammals in mass range 0.018 kg mouse to 2.067 kg rabbit. Else and Hulbert data indicate body mass independent MiVf of cells within brain and hence assumed valid for humans of 84 kg. Reference [27], for Hepatocytes of 0.020 kg mice to 200 kg horse, 26 for Cardiocytes of 0.003 kg shrew to 920 kg cattle and 28 for renal medullary cells of 0.011 kg bats to 400 kg horse. Note 4: Here constants remain same for Wang-5 or Wang-7.
Table 4. Derived Allometric Relations-Wang 7 Elia Mod (Values in parentheses are for 5 organ data). nk = 1.75 × 108 cells/cm3 (assumed), vfcell = cell volume/organ volume = (nk π dcell3/6).
Table 4. Derived Allometric Relations-Wang 7 Elia Mod (Values in parentheses are for 5 organ data). nk = 1.75 × 108 cells/cm3 (assumed), vfcell = cell volume/organ volume = (nk π dcell3/6).
Organ kEkFkHkIkBkOkPkQkJkLkNk
- q ˙ k , m ,   W k g   o r g a n   k =   E k   m k F k   E k = ( e k c k f k d k ) , F k =   ( f k d k ) Note 1, Note 2 V ˙ B l , k , m , m L   o f   b l o o d s   k g   o r g a n   k =   H k m B I k , H k = h k / c k I k = i k d k ,   Note 1 q · k , m ( t ) ( W   of k kg organ   k ) = B k ( Mivf k ) O k ,   B k = e k p k f k q k ,   O k = e k p k f k q k , q k C , B k = e k p k f k , O k = f k , q k = 0 q · k , v , Mito ( W cm 3   of   Mito   organ   k ) = P k m B Q k ,   P k = ρ k e k vf cell p k ,   Q k = f k q k Note 1 O E F = J k m B L k J k = ( c k e k h k   ρ b l Y O 2 i n   H V O 2 ) , L k = f k + d k i k Note 1 ( T e x i t T i n ) ,   ° C = N k m B L k   N k = J k Y O 2 i n H V O 2 c p or   N k = ( c k e k h k   ρ b l c p ) Note 1
AT0.22 (-)0.0 (-)- (-)−1.19 (-)0.22 (-)0.0 (-)- (-)- (-)- (-)- (-)- (-)
Brain11.62 (9.42)0.0 (−0.184)0.985 (9.182)−0.006 (−0.056)11.62 (14.36)0.0 (−1400) Note 22.44 (4.544)0.0 (−0.14)2.623 (0.524)0.006 (−0.084)3.090 (0.617)
Heart21.31 (23.04)0.0 (−0.122)37.333 (37.333)−0.266 (−0.266)22.53 (39.16)0.0 (2.72)0.874 (1.769)0.044 (−0.076)0.127 (0.257)0.266 (0.146)0.150 (0.303)
Kidney21.3 (20.95)0.0 (−0.094)65.753 (86.143)0.064 (−0.076)21.99 (32.95)0.0 (1.43)0.434 (0.681)0.056 (−0.024)0.07 (0.086)−0.064 (−0.004)0.0824 (0.102)
Liver9.69 (11.49)0.0 (−0.310)3.177 (4.727)0.156 (−0.014)11.31 (24.32)0.0 (3.00)0.623 (2.128)0.09 (−0.18)0.678 (1.558)−0.156 (−0.256)0.799 (1.836)
Residual0.58 (1.43)0.0 (−0.168)- (-)−0.994 (−1.01)- (-)- (-)- (-)- (-)- (-)- (-)- (-)
SM0.63 (-)0.0 (-)1.643 (-)−0.253 (-)0.84 (-)0.0 (-)- (-)- (-)- (-)- (-)- (-)
Note: 1: ck, dk, ek and fk for AT and SM are not defined for 5 organ model. Note 2: When allometric laws are given say for mk and q ˙ k , m in terms of body mass, mB, then on conversion to organ mass-based allometry; Fk → ∞ when dk → 0; but still q ˙ k , m will tend to correct limit. For this case, one may set dk to a, low value say 10−4 in getting appropriate numbers. Similar approaches are adopted for other parameters.
Table 5. Selected Quantitative Results for Steady state (84 kg person) and over Life Span for 5 organ (in parenthesis) and 7 organ data with metabolic efficiency = 0.31; however for blood temperature, metabolic efficiency set to zero to get maximum possible values.
Table 5. Selected Quantitative Results for Steady state (84 kg person) and over Life Span for 5 organ (in parenthesis) and 7 organ data with metabolic efficiency = 0.31; however for blood temperature, metabolic efficiency set to zero to get maximum possible values.
OrganMax Temp Rise, COrgan Mass, kg SteadyMiVf-steady (Note 3) q ˙ k , m , s t (W/kg Organ Mass)-Steady q ˙ k , V , M i t o , s t W/cm3 of Mito (Note 1)Normalized Energy Release Rate W/cm3 of Mito in k/W per cm3 of Heart-Steadyk (W by k) q k , M ,   l i f e MJ to One kg Body) σ k , M ,   l i f e (kJ/K) to One kg Body Mass σ k , m ,   l i f e MJ/K per kg Organ MassHeart Normalized Entropy Stress during Lifetime (Note 1)
AT-14.620 (-)-0.22- 3.216 (NA)86.4 (NA)191.8 (NA)1.14 (-)1.0% (-)
Brain3.174 (0.425)2.382 (0.319)0.0538 (0.0377)11.62 (11.63)2.44 (2.44)2.3 (1.93)27.680 (3.71)825.8 (113.5)1832.8 (252.0)60.32 (62.25)54.5% (46.1%)
Heart0.486 (0.578)0.461 (0.461)0.232 (0.283)21.31 (25.33)1.06 (1.26)1 (1)9.829 (7.09)274.8 (335.6)609.9 (744.8)110.6 (135.00)100% (100%)
Kidney0.062 (0.099)0.207 (0.303)0.439 (0.204)21.31 (23.44)0.556 (0.612)0.524 (.485)4.406 (7.09)131.4 (208.2)291.7 (462.1)110.6 (123.79)100% (91.7%)
Liver0.400 (0.590)1.092 (1.558)0.121 (0.155)9.69 (10.01)0.928 (0.959)0.873 (0.76)10.578 (15.60)316.4 (480.4)702.3 (1066.2)50.30 (55.37)45.5% (41.0%)
Residual Mass-26.617 (82.45)-0.63 (0.68)-- (-)15.438 (56.17)430.4 (1620.6)955.2 (3596.9)3.01 (3673.4)2.7% (2.7%)
SM0.308 (-)37.608 (-)0.055 (0.0)0.58 (0.58)- (-)0.455 (0.455)23.693 (56.17)661.1 (NA)1467.3 (NA)3.27 (-)3% (-)
Sum-82.99 (85.09) (Note 1)----94.8 (94.2)2726.5 (2758.4)6051.0 (6121.9)--
Note 1: assumed n = 1.75 × 108 cells/cm3. Note 2: Average cell volume Vavg,cell = 1000 μm3 = 10−9 cm3, ∆T for blood around heart (Tout − Tin) through blood = 0.58 C; ∆T for Mito of heart in C. Note 3: MiVf is under estimated since it is body mass-based and body mass is 84 kg.

Share and Cite

MDPI and ACS Style

Annamalai, K.; Nanda, A. Biological Aging and Life Span Based on Entropy Stress via Organ and Mitochondrial Metabolic Loading. Entropy 2017, 19, 566. https://doi.org/10.3390/e19100566

AMA Style

Annamalai K, Nanda A. Biological Aging and Life Span Based on Entropy Stress via Organ and Mitochondrial Metabolic Loading. Entropy. 2017; 19(10):566. https://doi.org/10.3390/e19100566

Chicago/Turabian Style

Annamalai, Kalyan, and Arnab Nanda. 2017. "Biological Aging and Life Span Based on Entropy Stress via Organ and Mitochondrial Metabolic Loading" Entropy 19, no. 10: 566. https://doi.org/10.3390/e19100566

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop