Aging in a Relativistic Biological Space-Time

1. Introduction

The connection between one's chronological age and biological age is something that we all perceive. In a sense, it is the difference between the age you “feel” and the age you are. Some people look “young” for their age, while some become frail earlier than others (Ness et al., 2013). Molecular “clocks” and markers of “biological age" can change throughout one's lifetime. In some cases, the rate of change of a person's biological age is greater than the rate of change of their chronological age. For instance, some cancers have been shown to increase biological aging, which can be described as an increase in the speed of biological clocks (Horvath, 2013).

From a mathematical perspective, biological dynamics as they relate to aging are frequently modeled as periodic or oscillating “clocks” describing, for example, circadian rhythms (Klerman and Hilaire, 2007). When studied in isolation, biological clocks can be described and predicted with periodic functions, which may speed up, slow down, or even stop and start over the course of a persons' lifetime. The more complex case of multiple integrated biological clocks can be modeled with coupled oscillators and dynamical systems theory (Shiju and Sriram, 2017). However, a fundamental assumption of these approaches is that the rate of change of these biological clocks are measured with respect to a linear passage of time and that the rate is independent from the biological space in which the biological clocks operate.

Here we investigate a mathematical model of biological space-time which includes the effects of time dilation and contraction resulting from accelerated or decelerated biological clocks, which may provide a new theoretical foundation and perspective on rates of aging. The principle assumption of our model is that a dynamic biological process may be represented as the motion of a (massless) point along a trajectory on a manifold. We then investigate the consequence of time dilation and contraction in terms of accelerating or decelerating motion of the point along the trajectory, and relate these concepts to rates of biological aging, with a particular focus on aging in cancer. In particular, we put forth the hypothesis that a biological space-time may be used to model aging from an arbitrary biological viewpoint relative to a common frame of reference. A consequence of this hypothesis is that precipitous (faster than chronological time) or protracted (slower) aging may be modeled as relativistic corrections of dilation or contraction of time along the trajectory on the manifold in which the biological aging process occurs.

To the best of our knowledge, only a few groups have proposed similar concepts. Bailly and colleagues (Bailly et al., 2011; Longo and Montévil, 2014) have proposed a mathematical definition of “biological time" as a means to model biological rhythms and periodic biological processes. However, they do not define a biological “space" nor include the possibility of the dilation or contraction of time. Systems biology pioneer Denis Noble has proposed a theory of biological relativity, which asserts that there is no “privileged level of causation" in biology (Noble, 2012). Nobel's theory contends that biological processes occur on many scales in space and in time, and that these scales are coupled to each other and should not be separated; that no single scale is responsible for the dynamics of the whole. Noble's theory implicitly couples space and time, but does not include a definition of biological space or concept of dilation or contraction of time. Consequently, our work is among the first to introduce both the mathematical interpretation and specific concept of a relativistic biological space-time to be used to study rates of biological aging.

This manuscript is structured as follows: first we describe the mathematical objects which we use to define biological space-time. After defining the meaning of relativistic dynamics in this context, we show how relativistic corrections of dilation and contraction of time can be used to model precipitous or protracted aging. We then show examples of how these principles may be used to model aspects of the aging process with a particular emphasis on aging in cancer. We discuss the implications of this theory, including criteria that may be used to test or falsify the theory and suggest novel biological quantities that may improve our interpretation of biological aging.

2. Biological Space-Time

Biology, and biological processes, are measured and observed in our conventional notion and understanding of physical space. Cells, tissues, and organisms move and change in a physical space that we can measure with length and time scales in conventional units. However, we may also consider the functional, or phenotype space in which biological processes can be represented as locations in the space. We refer to movement in a biological space as a sequence of locations in the space that form a trajectory. These general concepts have been used to characterize biological states such as hematopoietic differentiation, where 2- or 3-dimensional representations of biological space are constructed with dimension reduction techniques applied to high-dimensional single cell RNA-sequencing data (Mojtahedi et al., 2016; Nestorowa et al., 2016; Rizvi et al., 2017). The idea of the relativity of aging is to apply the special relativity machinery to provide a rigorously defined mathematical framework to represent biological dynamics as trajectories on manifolds moving at different acceleration rates relative to a common frame of reference. In other words, the difference in aging between two different people will be explained in terms of different dynamics in biological space-time.

2.1. A Manifold 𝔐 and Submanifolds 𝔐_i

In order to provide a conceptual picture of our mathematical framework, we imagine the space related to a biological process identified by the index i to be a smooth manifold ${{M}}_i$. A manifold is a mathematical object that is locally, but not necessarily globally, Euclidean. The canonical example of a manifold is the Earth: locally flat, globally round. If a point on Earth represents an unique locus, a point on ${{M}}_i$ represents a unique possible configuration or a state of a biological process. The evolution of the biological process in time, is then represented by a collection of points which form a curve, or trajectory, on the manifold ${{M}}_i$. If the biological process does not change in time, the trajectory degenerates to a point, with the consequent lack of motion in the biological space. We then identify the i-th biological space-time manifold by mathematically combining the time component and the spatial components of the i-th biological process.

Given a subset U_i of a topological space (i.e., a set in which at each point it is possible to associate a neighborhood) ${{M}}_i$, a d_i-dimensional chart is an injective (one-to-one) function ${\phi }_i:U_i\subset {{M}}_i\to {ℝ}^{d_i}$. A point q_i on ${{M}}_i$ is identified by a set of d_i spatial coordinates $(q_i^1,q_i^2,\dots q_i^{d_i})$. An atlas on ${{M}}_i$ is the set ${{A}}_i=\left\{{\phi }_{i,\alpha }:U_{i,\alpha }\to {ℝ}^{d_i}\right\}$ for some finite values of the index α, where the union U_{i, 1}∪U_{i, 2}∪…  is the whole space ${{M}}_i$. A space ${{M}}_i$ equipped with an atlas ${{A}}_i$ is a d_i-dimensional differential manifold. For a more detailed discussion see (Tu, 2011).

A trajectory, or curve, γ_i on ${{M}}_i$ is a smooth map:

$\begin{array}{ll}{\gamma }_i:\hfill & {{I}}_i\subset {ℝ}^{+}\to {{M}}_i\hfill \\ \hfill & t\mapsto {\gamma }_i(t)=q_i(t)=(q_i^1(t),q_i^2(t),\dots ,,q_i^{d_i}(t)).\hfill & (1)\end{array}$

Here $t\in {{I}}_i\subset {ℝ}^{+}$ is a parameter that is interpreted as the time variable associated to the manifold ${{M}}_i$. The curve γ_i represents the time evolution of the i-th biological process starting from its beginning to its end (Figure 1). If the point q_i(t) does not change for all $t\in {{I}}_i\subset {ℝ}^{+}$, the trajectory degenerates to a point. We stress the fact that the variable t is a parameter by which is possible to parametrize the curve γ_i(t).

FIGURE 1

Figure 1. The trajectory γ_i on the manifold ${{M}}_i$ represents the evolution of the i-th biological process which starts at q_i(a_i) and ends at q_i(b_i), where a_i, $b_i\in {ℝ}^{+}$ . The atlas ${{A}}_i$ provides the connection with the Euclidean space ℝ² identified by the two unit vectors e_x and e_y.

Given the manifold ℝ⁺, we now define the i-th biological space-time manifold 𝔐_i as the (d_i + 1)-dimensional Lorentzian manifold given by the Cartesian product (the set of all ordered pairs) (see Tu, 2011):

$\begin{array}{ll}{\mathfrak{M}}_i={ℝ}^{+}\times {{M}}_i=\left\{(t,q_i(t))|t\in {ℝ}^{+},q_i(t)\in {{M}}_i\right\}.& (2)\end{array}$

A point ${{Q}}_i$ on the resulting manifold 𝔐_i is then identified by a set of d_i + 1 coordinates $(t,q_i^1(t),q_i^2(t),\dots ,q_i^{d_i}(t))$. The invariant square of an infinitesimal line element between the points ${{Q}}_i$ and ${{Q}}_i+d{{Q}}_i$, referred to as space-time interval, is then evaluated by:

$\begin{array}{ll}d{s}_i^2={\displaystyle \sum _{\mu ,\nu =0}^{d_i}{g}_{i,\mu \nu }}d{{Q}}_i^{\mu }d{{Q}}_i^{\nu },& (3)\end{array}$

where g_{i, μν} is the metric tensor of the i-th biological manifold whose entries are functions of the local coordinates, $(t,q_i^1(t),q_i^2(t),\dots ,q_i^{d_i}(t))$. Since the metric tensor is an intrinsic property of the manifold, its provides information on the geometry of the manifold and vice-versa. It is therefore a key quantity which characterizes the biological space-time and it can be evaluated once the structure of the manifold is known or hypothesized. What we emphasize is that the notion of distance is not necessarily the usual Euclidean distance and therefore, when dealing with biological processes, distances can be measured only once the entries of the metric tensor are known.

Since a large number of biological processes take place in the body of an individual, we assume the existence of several manifolds which may be indexed ${{M}}_i$ and associated with individual or aggregate biological processes. We then consider all such biological processes by constructing the following manifold:

$\begin{array}{ll}{M}={{M}}_1\times {{M}}_2\times \dots \times {{M}}_N,\text{dim}({M})={\displaystyle \sum _{i=1}^N{d}_i}=m.& (4)\end{array}$

A point q on ${M}$ is then identified by the set of coordinates (q₁, q₂, …, q_N) where each $q_i\in {{M}}_i$, i = 1, 2, …, N and it represents a possible configuration, or state, in which the human body can be found. Following the definition given by Equation (1), a trajectory Γ on ${M}$ is a smooth map:

$\begin{array}{ll}\begin{array}{ll}\Gamma :\hfill & {I}\subset {ℝ}^{+}\to {M}\hfill \\ \hfill & t\mapsto \gamma (t)=q(t)=(q_1(t),q_2(t),\dots ,q_N(t)),\hfill \end{array}& (5)\end{array}$

where $t\in {I}\subset {ℝ}^{+}$ is a parameter that is interpreted as the time variable associated to the manifold ${M}$ and the curve Γ represents the entire life of a person from its birth, to its death. The connection between the curve Γ on ${M}$ and the curve γ_i on ${{M}}_i$ will be clarified later, when the projection map will be introduced.

In analogy with Equation (2), the biological space-time 𝔐 for all biological processes occurring in a body of an individual is then defined by the m + 1-dimensional Lorentzian manifold given by:

$\begin{array}{ll}\mathfrak{M}={ℝ}^{+}\times {M}=\left\{(t,q(t))|t\in {ℝ}^{+},q(t)\in {M}\right\}.& (6)\end{array}$

We consider the Cartesian product $𝔐={ℝ}^{+}\times {M}$ instead of 𝔐 = 𝔐₁ × 𝔐₂ × ⋯ × 𝔐_N because in the latter case it is unclear how all time variables, one per each submanifold 𝔐_i, will combine together and define the time on the resulting manifold 𝔐. A point ${Q}$ on the manifold 𝔐 is then identified by a set of coordinates as follows:

$\begin{array}{ll}{Q}=(\stackrel{\in {ℝ}^{+}}{\overbrace{t}},\underset{\in {{M}}_1}{\underbrace{q_1^1,q_1^2,\dots q_1^{d_1}}},\underset{\in {{M}}_2}{\underbrace{q_2^1,q_2^2,\dots q_2^{d_2}}},\dots ,\underset{\in {{M}}_N}{\underbrace{q_N^1,q_N^2,\dots q_N^{d_N}}}).& (7)\end{array}$

The idea for which time flows at different rates for different biological processes is now mathematically modeled by introducing a set of projection maps Π_i such that:

$\begin{array}{ll}\begin{array}{ll}{\Pi }_i:\hfill & {ℝ}^{+}\times {M}\to {ℝ}^{+}\times {{M}}_i\hfill \\ \hfill & {Q}\mapsto {{Q}}_i={\Pi }_i({Q})=({\tau }_i(t)=t_i,{\pi }_i(q)=q_i),\hfill \end{array}& (8)\end{array}$

where the spatial part ${\pi }_i:{M}\to {{M}}_i$ is the canonical projection map which maps the curve Γ on ${M}$ onto the curve γ_i on ${{M}}_i$ while the temporal part ${\tau }_i:{ℝ}^{+}\to {ℝ}^{+}$ provides the connection between the time t on the manifold 𝔐 and the time t_i on the i-th manifold 𝔐_i (Figure 2).

FIGURE 2

Figure 2. A torus is shown as an example of a possible biological space-time manifold that may be decomposed into submanifolds. The time ℝ⁺ is represented by the vertical arrows. The torus is decomposed into the two circles (submanifolds) $M_1=S_1^1$ and $M_2=S_2^1$. The point ${{Q}}_0$ at time t₀ is mapped onto the points ${{Q}}_{1,0}$ and ${{Q}}_{2,0}$ while the point ${{Q}}_1$ at time t₁ is mapped onto the points ${{Q}}_{1,1}$ and ${{Q}}_{2,1}$ in the two space-time submanifolds ${𝔐}_1={ℝ}^{+}\times S_1^1$ and ${𝔐}_2={ℝ}^{+}\times S_2^1$. The trajectory on the torus is then mapped onto two different trajectories on the two circles from the initial point (green circle) to the final point (red circle).

We emphasize that the relationship between t and t_i is not necessary linear. In fact, in the next section we will show that the time measured in a particular manifold depends upon the acceleration of the particle along a trajectory and it will produce a non-linear relation between t and t_i.

2.1.1. Example: Decomposing the Torus

We now provide a simple example which can be easily visualized in a three-dimensional space. We consider a manifold, the torus embedded in the three dimensional space ℝ³, which can be decomposed into two circles, the two submanifolds. In this particular example, this is the only possible decomposition (Tu, 2011).

The torus, embedded in ℝ³, is defined by the following set of equations:

$\begin{array}{ll}\{\begin{array}{l}x(\theta ,\varphi )=(R+r\text{ cos }\theta )\text{ cos }\varphi \\ y(\theta ,\varphi )=(R+r\text{ cos }\theta )\text{ sin }\varphi \\ z(\theta ,\varphi )=r\text{ sin }\varphi \end{array}& (9)\end{array}$

and is decomposed into the two submanifolds $S_1^1$ and $S_2^1$ given by:

$\begin{array}{ll}{S}_1^1=\left\{r{e}^{i\theta }:\theta \in [0,2\pi )\right\},S_2^1=\left\{R{e}^{i\varphi }:\varphi \in [0,2\pi )\right\},& (10)\end{array}$

where, θ identifies the poloidal direction (along the orange vertical circle in Figure 2), ϕ identifies the toroidal direction (along the blue horizontal circle in Figure 2), R the major radius (distance from the center of the torus), and r < R is the minor radius of the torus (radius of the tube). These two circles can represent two submanifolds related to two different biological processes, and their Cartesian product produces the torus. Therefore, the two biological space-times associated to this particular example, and hence to these two biological processes, will be given ${𝔐}_1={ℝ}^{+}\times S_1^1$ and ${𝔐}_2={ℝ}^{+}\times S_2^1$. A trajectory Γ on the torus, identified by the equations:

$\begin{array}{ll}\{\begin{array}{l}x(t)=(R+r\text{ cos }(2\pi nt))\text{ cos }(2\pi mt)\\ y(t)=(R+r\text{ cos }(2\pi nt))\text{ sin }(2\pi mt)\\ z(t)=r\text{ sin }(2\pi mt)\end{array}& (11)\end{array}$

where t ∈ ℝ⁺, is then decomposed onto two trajectories:

$\begin{array}{l}{s}_1(t)=\left\{r{e}^{2\pi nit}:t\in [0,1)\subset {ℝ}^{+}\right\},\\ s_2(t)=\left\{R{e}^{2\pi mit}:t\in [0,1)\subset {ℝ}^{+}\right\},& (12)\end{array}$

where n and m are the winding numbers associated with the poloidal and toroidal directions, respectively. The two submanifolds (two circles) have two different trajectories defined by Equation (12). However, when their Cartesian product is considered, the resulting manifold (the torus) has a single unique trajectory on it defined by Equation (11). In this particular example, the magnitude of the velocity and the acceleration of the projected motion on the two submanifolds are given by:

$\begin{array}{ll}\{\begin{array}{l}{v}_1=2n\pi r\\ v_2=2m\pi R\end{array}\text{and}\{\begin{array}{l}{a}_1=4n^2{\pi }^{2}r\\ a_2=4m^2{\pi }^{2}R\end{array}& (13)\end{array}$

which implies a₂ > a₁ for the same value of n and m. The position, the velocity and the acceleration on the torus are projected onto the two subspaces in which the dynamics are identified by different values for the position, velocity and acceleration. If in Equation (11) we set m = 0, the trajectory reduces to a circle around the poloidal direction. In this case, Equation (13) become:

$\begin{array}{ll}\{\begin{array}{l}{v}_1=2n\pi r\\ v_2=0\end{array}\text{and}\{\begin{array}{l}{a}_1=4n^2{\pi }^{2}r\\ a_2=0\end{array}& (14)\end{array}$

This example shows that there can be motion on the torus and on one submanifold, but an absence of motion on the other submanifold.

The example of the torus illustrates how trajectories on submanifolds (e.g., representations of different biological dynamics) with different accelerations may be combined and interpreted as a single trajectory on a larger manifold, and vice versa. It must be noted that in this decomposition, the time t in the submanifolds is a parameter along the curve and not the biological time. In fact, as we will see in the next section, because the accelerations in the two submanifolds are different we expect time to flow at different rates in each submanifold.

3. Dilation and Contraction of Time

The fact that the Maxwell's equations are invariant under the Lorentz transformation implies that any inertial observer will measure light moving at the same constant speed c. Moreover, the Lorentz transformations define an invariant quantity, called proper time, which represents a space-time interval which assumes the same value for any inertial (non-accelerating) observer. In physics, the presence of fundamental equations such Maxwell's equations, helps us in understanding the type of coordinate transformations that leave the equations unchanged and can be used to define invariant quantities. Here, we deal with a biological space-time for which the existence of fundamental equations, coordinate transformations, and corresponding invariant quantities, are unknown.

Although many investigators have posed the question of whether or not governing laws exist in biology (Ruse, 1970; Brandon, 1997; Wagner, 2017), we contend this question has yet to be conclusively answered. Rather, we believe a more fundamental question is to investigate the nature of the space-time in which any biological phenomenon occurs. In the context of Special Relativity, it has been shown by Levy-Leblond (1976) that if the space is homogeneous and isotropic, the Lorentz transformations are characterized by a parameter with the dimensions of a velocity which is an intrinsic property of the space. Moreover, the meaning of this parameter is related to the maximum velocity at which information can travel or be transmitted (Brunner et al., 2004) and this velocity turns out to be the speed of light c.

3.1. Biological Invariants and Information

Although information in the context of biological processes can not be easily defined (Gatenby and Frieden, 2007), it can be related to the efficiency in the conversion of energy into order (Frieden and Gatenby, 2011). In our biological context, assuming our biological space-time to be homogeneous and isotropic, there will be a parameter analogous to c, denoted by $\tilde{c}$, which can be related to the rate of conversion of energy into order. We stress the fact that, although we are not sure about the meaning of information in this particular context, we assume that the meaning of $\tilde{c}$ is to provide an upper limit for the information to be transmitted. This upper limit is what we believe is an invariant quantity in the biological space-time.

As described in section 2.1, a point $q\in {M}$ defines the state of a biological process. In what follows we assume the manifold to be a flat plane and to be constructed such that two points q(t) and q(t + dt) represent the information of the system at time t and t + dt, respectively (see Figure 3). The existence of a value $\tilde{c}$ defines a circular region Ω, centered at q(t) with boundary ∂Ω given by the circumference with radius $r=\tilde{c}dt$, such that $\left|\right|q(t)-q(t+dt)\left|\right|<\tilde{c}dt$. In other words, the value of $\tilde{c}$ defines a front of information which propagates in time along the manifold (see Figure 3). In the case of a generic manifold ${M}$, we will have a generic region Ω bounded by a curve ∂Ω which will depend on the geometry of the manifold. In this case the distance can be evaluated once the metric tensor g_μν is known. Moreover, the existence of such a parameter defines the non-relativistic and the relativistic limits in this context. The former occurs for velocities $v<<\tilde{c}$ while the latter occurs for velocities $v~\tilde{c}$ (see Figure 3, top and bottom left). In analogy with the Special Relativity case, we refer to relativistic corrections, as those corrections which need to be made in order to correctly describe and characterize the dynamics of the moving particle with velocity $v~\tilde{c}$. We want to clarify the fact that at this stage we do not have a procedure to build a manifold for a given biological process, rather, we assume the existence of the manifold.

FIGURE 3

Figure 3. Top left. Representation of the information front ∂Ω propagating at speed $v=\tilde{c}$ from the center of the disk Ω on the flat plane ${M}$. The distance traveled by a point in the interval of time dt is given by ||q(t + dt) − q(t)|| and it will be always less than $\tilde{c}dt$ (dashed yellow line), where $\tilde{c}$ is the maximum speed of information. Bottom left. The analogous of the relativistic corrections are needed if we want to fully characterize the dynamics of a particle which is moving at speed $v~\tilde{c}$ Right. In the case of a generic manifold ${M}$ the information front ∂Ω will be the boundary of a region Ω.

3.2. Special Relativity and Biological Processes

Special relativity is universally recognized as a theory which describes properties of the ordinary space-time in which physical phenomena occur. In the case of a biological process, unlike the physical phenomenon, the biological space-time is not known a priori and it needs to be mathematically constructed. In this section we assume the existence of such a space and we consider the motion of two frames of reference whose coordinates are related by a set of coordinate transformations. In particular we investigate the dilation and contraction of biological time resulting from accelerated motions. This approach is used to model and to explore different rates of aging.

As shown by Levy-Leblond (1976), Lorentz transformations are the only set of transformations of coordinates in a homogeneous and isotropic space which also preserve the group structure and the causality condition. Therefore, assuming our space to be homogeneous and isotropic, we assume their validity also in a biological space-time. In this case, instead of the speed of light c, the Lorentz transformations will depend on the biological invariant $\tilde{c}$.

We consider two frames of reference S and S′ in the flat Minkowski (the space given by considering the time component and the ordinary three-dimensional space) space-time M⁴ = ℝ × ℝ³. In general, the motion of a particle is along a curvilinear trajectory, and the two frames of reference can be in a relative motion like the one represented in Figure 4 on the left. However, for simplicity, we assume a relative motion with a constant relative velocity v along the x-direction. We denote with (t, x, y, z) and (t′, x′, y′, z′) an event recorded by S and S′, respectively. These two events are linked by the Lorentz transformations

$\begin{array}{ll}\begin{array}{l}{x}^{\prime }=\gamma (x-vt),y^{\prime }=y,z^{\prime }=z,t^{\prime }=\gamma \left(t-\frac{vx}{{\tilde{c}}^2}\right),\hfill \\ \gamma =\frac{1}{\sqrt{1-v^2/{\tilde{c}}^2}},\hfill \end{array}& (15)\end{array}$

where γ > 1 is the Lorentz factor and $\tilde{c}$ is the biological invariant previously introduced. It is straightforward to show that these transformations leave invariant the space-time interval $d{\tau }^2={\tilde{c}}^{2}d{t}^2-d{x}^2-d{y}^2-d{z}^2$ which we will call biological proper time. The invariance of the biological proper time will lead, also for biological systems, to the well known phenomenon of time dilation:

$\begin{array}{ll}d{t}^{\prime }=\gamma dt.& (16)\end{array}$

The time interval dt′ measured from an observer in motion is dilated by the factor γ and we say that a moving clock runs slower. It has to be remarked that this effect is reciprocal when both frames of reference are inertial (there is a relative motion at constant speed between them). In that case, no preferred frame of reference exists and each observer can argue that the other one is moving. For this reason, we do consider accelerating frames of reference in which this ambiguity is solved. We now consider the particular case in which a point is moving along a straight line oriented along the x-axis with a constant acceleration $a_x^{\prime }$ (Figure 4 on the right). We assume S′ to be located at the position of the particle and therefore, its velocity is the same as the velocity of the particle and it increases in time. In this particular situation, it can be shown, (see Appendix A), that the time t′ of the accelerated particle, and hence in the accelerated frame of reference, is given by:

$\begin{array}{ll}{t}^{\prime }=\frac{\tilde{c}}{{a^{\prime }}_x}\text{arcsinh}\left(\frac{{a^{\prime }}_{x}t}{\tilde{c}}\right).& (17)\end{array}$

This equation tells us how much time t′ elapses in the accelerated frame of reference S′, in terms of the time t elapsed in the inertial frame of reference S. This function represents the time component τ_i(t) of the projection map Π_i introduced in Equation (8). In the case of an accelerating frame of reference, the phenomenon of time dilation is not reciprocal and the observer in the frame of reference S′ who experiences the constant acceleration cannot argue that the other observer is moving and he is still. Therefore, the time t′ of the accelerating frame of reference measured from the observer at rest in the frame of reference S appears to run slower according to Equation (17). In what follows we will apply this idea to few biological processes always assuming that the time dilation is described by Equation (17).

FIGURE 4

Figure 4. On the left) The particle (purple dot) is moving along a trajectory with velocity u and acceleration a. On the right) The frame of reference S′ is moving at constant velocity v away from S: the components of the velocity and the acceleration of the particle in these two frames of references are related by Equations (A.1) and (A.2).

3.3. Biological Age

Many markers have been used to define biological age. In our modeling framework, biological time is defined as the rate at which biological processes take place, as measured against chronological time. Therefore, our model characterizes one's biological age by the degree of dilation or contraction of time resulting from the acceleration or deceleration of biological processes that are associated with biological age.

One marker of biological age is the methlyation state of the genome, composed of varying degrees of hyper- or hypomethlyated states (Bocklandt et al., 2011; Horvath, 2013). The “epigenome" is a component of gene regulation and has been proposed to contain information relevant to the overall state of the biological system (Jenkinson et al., 2017). Linear age-related epigenetic drift has been associated with cancer incidence (Curtius et al., 2016), and statistical methods have been proposed to identify deviation from time linearity in epigenetic aging, although without a theoretical rationale to describe—or predict—the nature of the non-linearity (Snir et al., 2016).

In order to explain the process of changing methylation states with age, and the use of methylation state as a surrogate marker of biological age with our model, we first consider the trajectory on a manifold which describes the methylation process of a living person and we assume that this trajectory is mapped onto a straight line in the Minkowski space. We now consider one frame of reference S at rest in the origin of the trajectory and one frame of reference S′ that is moving at constant acceleration $a_x^{\prime }$ for which the time t′ is given by Equation (17). We then assume that the frame of reference slows down to an acceleration $a_x^{″}<a_x^{\prime }$ for which the time is t″ and it is given by Equation (17). We define the rate of methylation dM/dt′ and dM/dt″ when the frame of reference moves with acceleration equal to $a_x^{\prime }$ and $a_x^{″}$, respectively and we have:

$\begin{array}{ll}\frac{dM}{d{t}^{\prime }}=\frac{dM}{dt}\frac{dt}{d{t}^{\prime }},\frac{dM}{d{t}^{″}}=\frac{dM}{dt}\frac{dt}{d{t}^{″}}.& (18)\end{array}$

If $a_x^{\prime }>a_x^{″}$ it can be shown that dt′ / dt < dt″ / dt, and therefore, we conclude that:

$\begin{array}{ll}\frac{dM}{d{t}^{\prime }}>\frac{dM}{d{t}^{″}}.& (19)\end{array}$

The slowing down of the rate of demethylation is therefore interpreted as a decelerating frame of reference. We note that this model does not assume any specific functional form for the methylation trajectory. In fact, during the lifetime of a person, Equation (18) permits the acceleration or deceleration of age-related changes to the methylation state, which can be accentuated or modified in the context of cancer.

4. Examples of Space-Time Dynamics in Biology

In this section we connect our model to different rates of aging between individuals, and discuss examples of precipitous or protracted aging within an individual.

4.1. Aging

We define aging as the functional and structural decline of an organism, resulting in an increasing risk of disease, impairment and mortality over time. At the molecular and cellular level several hallmarks of aging have been proposed to define common characteristics of aging in mammals (López-Otín et al., 2013). These hallmarks include: genomic instability, telomere attrition, epigenetic alterations, loss of proteostasis, altered nutrient-sensing, mitochondrial dysfunction, cellular senescence, stem cell exhaustion and altered cell-cell communication. These processes are interdependent and influenced by cell microenvironmental cues.

Ultimately, the rate of age-related decline varies depending on how genetic variation, environmental exposure and lifestyle factors impact these mechanisms. Consequently, age, when measured chronologically, is often not a reliable indicator of the body's rate of decline or physiological breakdown. Over the years, the idea of quantifying the “biological age" based on biomarkers for cellular and systemic changes that accompany the aging process have been explored (Levine, 2013), although without a rigorous mathematical treatment of biological age. Major efforts have been made to dissect the individual contribution of each of these hallmarks and factors on aging, but the major challenge remains how to determine how their interconnectedness as a whole impacts the aging dynamics. Here we connect concepts in our mathematical model to biological and chronological aging. We note that for the sake of clarity, we will refer to aging as precipitous (faster than chronological aging) or protracted (slower) in order to clearly associate the terms “accelerated" and “decelerated" to quantities in our modeling framework.

As we derived in the previous section, an unambiguous time dilation effect requires an accelerating frame of reference. We now imagine a range of accelerations representing a distribution within a population between the values a_min and a_max (green region in Figure 5 on the left). The dashed and dotted curves corresponds to Equation (17) by using the values a_min and a_max, respectively. Any acceleration within this range will correspond to a curve in the green region. An interval of time Δt on the x-axes corresponds to a reduced interval of time Δt′ on the y-axis and therefore, the ticking of the second clock is slowed down. The existence of a distribution, or range of values of acceleration is necessary to capture both features of time dilation, (green region in Figure 5 on the left) and time contraction, when the value of the acceleration is less than a_min (pink region). Moreover, values of the acceleration much larger than a_max can be used to define the arrest of biological time such as in a hibernating state or cryogenic freezing (blue region). The case of zero acceleration is given by Equation (16) and is not considered in this model context.

FIGURE 5

Figure 5. On the left) For an accelerating frame of reference a time interval Δt = 1 corresponds to shorter Δt′ with an consequent dilation of time. The green region represents a distribution of accelerations associated with an average range of accelerations that is used as a common reference for accelerated and decelerated processes. Smaller values of the acceleration are related to time contraction (pink region, $a_x^{\prime }<a_{min}^{\prime }$) while large values are associated with time dilation (blue region, $a_x^{\prime }>a_{max}^{\prime }$). On the right) In the mental illustration of the lifetime of Bob, a deceleration will bend Bob's trajectory (black line with red diamonds markers) out of the green region with a corresponding contraction of time.

4.1.1. Alice and Bob

We illustrate the contraction and dilation of time in biological space-time (i.e., on a submanifold) and how they can be related to the presence of a disease, we consider the following thought experiment involving Alice and Bob, two individuals with different rates of aging. For simplicity, we assume that their dynamics occur on a flat torus ${{T}}^2$ defined as the Cartesian ℝ² plane under the identifications (x, y) ~ (x + 2π, y) ~ (x, y + 2π). The opposite edges x, x + 2π and y, y + 2π of the domain are identified and therefore, they must be interpreted as the same point (see red and blue edges in Figure 6). In other words, the trajectory that comes out from the right edge of the plane, will come in from the left edge and the same applies from the upper and lower part of the plane. In this case, the manifold is trivially mapped into ℝ² with the identity map, and hence the dynamics on the manifold corresponds to the dynamics in the Minkowski space-time.

FIGURE 6

Figure 6. The flat torus ${{T}}^2$ in which the dynamics of Alice $(S_A^{\prime })$ and Bob $(S_B^{\prime })$ occur along the straight dotted line. The green and the purple dots represent their birth and their actual position on the manifold. The effect of time dilation is more evident for Alice whose frame of reference moves at higher acceleration.

The point on the manifold corresponding to the birth of Alice and Bob is indicated in Figure 6 by the green circle in which we assume to place a frame of reference S identified by the three unit vectors (e_x, e_y, e_z). The dynamics of Alice is represented by the motion of a purple point in which we imagine to place a frame of reference $S_A^{\prime }$. We also assume that the point is moving along a straight line with constant acceleration given by $a_A^{\prime }=a_A^{\prime }{e}_x^{\prime }$. The same conditions apply to Bob, but in this case the frame of reference is $S_B^{\prime }$ and the particle is moving with a constant acceleration $a_B^{\prime }=a_B^{\prime }{e}_x^{\prime }$ (see Figure 6).

We now can apply the results obtained in the previous section to infer that the time in $S_A^{\prime }$ and $S_B^{\prime }$ are given by Equation (17):

$\begin{array}{ll}{t^{\prime }}_A=\frac{\tilde{c}}{{a^{\prime }}_A}\text{arcsinh}\left(\frac{{a^{\prime }}_{A}t}{\tilde{c}}\right),{t^{\prime }}_B=\frac{\tilde{c}}{{a^{\prime }}_B}\text{arcsinh}\left(\frac{{a^{\prime }}_{B}t}{\tilde{c}}\right).& (20)\end{array}$

Assuming the acceleration of Alice to lie within the range [a_min, a_max], her proper time $t_A^{\prime }$ lies in the green region (see Figure 5 on the right). On the other hand, Bob's frame of reference is moving with acceleration $a_B^{\prime }<a_A^{\prime }$ and its proper time $t_B^{\prime }>t_A^{\prime }$ deviates from the green region toward an older biological age in the pink region. With respect to Alice, Bob is moving with a slower acceleration and hence, his clock is running faster, resulting in a time contraction with respect to Alice's clock. If instead we assume $a_B^{\prime }>a_A^{\prime }$, Bob's proper time $t_B^{\prime }$ will be larger than $t_A^{\prime }$ and it deviates from the green toward the blue region. In this case Bob's clock will run slow with a consequent more evident effect of time dilation. Therefore, different rates of aging could be generally compared between individuals by determining their accelerations as compared to a frame of reference at rest.

We now consider only Bob and we assume his frame of reference $S_B^{\prime }$ to have an acceleration $a_1^{\prime }$ within the range $\left[a_{min}^{\prime },a_{max}^{\prime }\right]$. Suddenly, his acceleration is decreasing to the value $a_2^{\prime }$ and we therefore want to know what will happen to its proper time t′. In particular we hypothesize the acceleration to change as follows:

${a^{\prime }}_B=\{\begin{array}{ll}{a^{\prime }}_1\hfill & t<t_1\hfill \\ \alpha +\beta \text{ cos }(\pi (t_1-t)/(t_2-t_1))\hfill & t_1\le t\le t_2\hfill \\ {a^{\prime }}_2\hfill & t>t_2\hfill \end{array}$

where $\alpha =(a_1^{\prime }+a_2^{\prime })/2$ and $\beta =(a_1^{\prime }-a_2^{\prime })/2$. By solving Equations (A.5) and (17) we obtain the time t′ for the decelerating frame of reference shown in Figure 5 on the right: the line marked with green circles and the line with violet squared markers represent the trajectory corresponding to the accelerations $a_1^{\prime }$ and $a_2^{\prime }$, respectively. However, due to the deceleration, the line marked with red diamonds will not follow the line with the green circles but it will move outside the green region. The deceleration of Bob's frame of reference is reflected to a contraction of time while a larger time dilation is a consequence of a frame of reference with an increasing acceleration.

4.2. Cancer as Inducer of Precipitous Aging

Major pathologies, such as cancer, diabetes, cardiovascular disorders and neurodegenerative diseases have an impact on aging. Cancer and chemotherapy in particular are known to accelerate the aging process (Alfano et al., 2017). As described above, aging involves multiple complex changes at molecular and cellular level that lead to decline in physiologic reserve capacity across virtually all organ systems (Seals et al., 2016). In cancer, accumulation of cellular damage aggravates the hallmarks of aging and accelerates aging. In cancer patients, the trajectory of decline worsens not only because of the direct physiologic insult inflicted by cancer, but also because of the injuries induced by anti-cancer therapies, such as chemotherapy or radiation, to systems that maintain physiologic reserves (Koelwyn et al., 2012; Hurria et al., 2016). In this context, there is a dose-dependent effect, whereby the more intensive the treatment of cancer or more vulnerable and frail the physiologic state, the steeper the decline in physiologic reserves (Henderson et al., 2014; Kohanski et al., 2016). Below, we use these two variables, intensity of therapy and frailty, to illustrate our model of accelerated aging.

It is not surprising that the two cancer populations most affected by precipitous aging caused by cancer are survivors of childhood cancer who are typically exposed to intensive multi-agent therapy at a young age (Henderson et al., 2014; Ness et al., 2015) and adult patients who undergo hematopoietic cell transplantation (HCT) for refractory hematologic malignancies (Arora et al., 2016). These two scenarios can be explained by our model as follows. The early onset of advanced biological age in childhood cancers corresponds to a rapid contraction of time that deviates away from the green region, similar to Bob's modified trajectory, but with an earlier onset (see Figure 5 on the right). The case of allogeneic HCT, in which case the transplanted cells come from another person, corresponds to a discontinuity in biological space-time. The recipient's trajectory on the hematopoietic manifold and/or the biological age is reset to that of the donor. If the biological ages of the donor and recipient of the HCT do not agree, the recipient will experience an abrupt and persistent change in the rate of biological aging corresponding to the discontinuity in both biological space and time.

Frailty, characterized by a cluster of measurements of physical states, is the best described measure of aging in a population, and identifies individuals who are highly vulnerable to adverse health outcomes and premature mortality. Although frailty is not a perfect corollary with biological age, it is a measure of abnormal aging at a population level. In long-term survivors of childhood cancer [median age 33 years (range 18–50 years)], the prevalence of frailty (Rockwood and Mitnitski, 2011) has been shown to be comparable to that reported among adults greater than 65 years of age (Ness et al., 2013). Comparable high rates have been reported in adult survivors of HCT (Arora et al., 2016). In fact, frail HCT recipients have a three-fold increased risk of subsequent mortality compared with the non-frail counterparts, (Arora et al., 2016) which is similar to the downstream consequences of frailty seen in the general population such as adverse health outcomes (Rockwood and Mitnitski, 2011), and early mortality (Fried et al., 2001; Hogan et al., 2003). From the lens of our model, frailty can be associated with a biological age and therefore defined as a threshold in time ($t_f^{\prime }$) so that the onset of frailty in the cancer population $(t_c,t_f^{\prime })$ will appear sooner relative to the general population $(t_g,t_f^{\prime })$ where t_c < t_g.

4.2.1. Chronic Inflammation

The normal process of aging is associated with chronic low grade inflammation and with cumulative oxidative stress, independently of disease. Inflammation and oxidative stress are critical responses in host defense and injury repair and are essential for normal body functions. However, with advanced age there is a loss of sensitivity in the injury-repair cycle leading to persistent chronic inflammation, and a natural decline in the endogenous anti-oxidant capacity leading to cumulative oxidative stress (Mittal et al., 2014). These two process are interdependent, and contribute to the hallmarks of aging, influencing telomere length, mitochondrial function, epigenetics and stem cell self-renewal.

Chronic inflammation and oxidative stress are also common underlying factors of age-associated diseases. Inflammation, particularly chronic low-grade inflammation, has been found to contribute to the initiation and progression of multiple age-related pathologies such as type II Diabetes, Alzheimer's disease, cardiovascular disease and cancer (Mantovani et al., 2008). In addition, chronic diseases, and in particular cancer, elicit and promote an inflammatory tumor microenvironment that increases cancer fitness. Therefore, chronic inflammation has been suggested to underlie and accelerate biological aging (Fougère et al., 2017), in particular when associated with disease. The concept of “inflamm-aging,” inflammation-associated aging, can be used to provide a systemic perspective of biological aging in the human population (Franceschi et al., 2007).

On the level of tissue homeostasis, the best characterized example is hematopoietic aging associated with chronic inflammatory signaling in the bone marrow microenvironment. In fact, the “age" of a young hematopoietic stem cell can be “reprogrammed" when transplanted into an aged or inflammatory environment (Kovtonyuk et al., 2016), highlighting the impact of "inflamm-aging" and the plasticity of molecular clocks. Aging is associated with clonal hematopoiesis and accumulation of mutations in hematopoietic progenitors (Steensma et al., 2015), likely due to the underlying inflammation in the aged bone marrow niche. Indeed, an inflamed bone marrow can induce pre-leukemic conditions in mice similar to those occurring in elderly patients (Wang et al., 2014; Dong et al., 2016). Our model can be applied to measure how microenvironmental cues, such as inflammation, impact aging of hematopoietic cells.

4.3. Protracted Aging

In contrast to precipitous biological aging which corresponds to the contraction of biological time, protracted aging corresponds to the dilation of biological time. This can be illustrated by prolonged periods of near zero biological activity, for instance in freezing conditions or hibernation which is part of a continuum of biological and metabolic states (van Breukelen and Martin, 2015). How can this situation be explained in a relativistic framework? We imagine a cryogenic-manifold mapped onto a Minkowski space. If in this space the frame of reference is moving with a high acceleration in the direction of motion then the interval of time in its frame of reference tends to zero: a frozen cell corresponds to a frame of reference moving at very large acceleration relative to its unfrozen state on this particular biological manifold. Figure 5 shows the limiting case of a near perfectly frozen (cryogenic) biological process (blue region), which would correspond to an extreme dilation of time, or equivalently, a very large relativistic correction.

5. Discussion

Here we have investigated a mathematical model of biological aging. We define a biological space-time by mathematically combining manifolds and submanifolds, and apply the principles of relativity to compare different rates of biological aging. We illustrate the concepts of precipitous and protracted aging as relativistic corrections of biological time with a mental illustration comparing the lifespans of two individuals, Bob and Alice. This analogy provides a framework to compare the rates of aging between individuals by determining their rate of acceleration as compared to a common frame of reference.

A critical component of our theory is the construction of “biological space-time,” hence a submanifold that represents a biological process. Following the work of (Levy-Leblond, 1976) we believe that the “biological space-time" should have the properties of being homogeneous and isotropic. If we also require that the transformations between frames of reference in this space preserve the group composition law and causality, the transformations are characterized by a parameter $\tilde{c}$, that has the dimensions of a velocity, which is an intrinsic property of the space (Levy-Leblond, 1976). In analogy with the case of the Lorentz transformation for which this parameter is the speed of light c and it represents the upper limit for the information to be transmitted (Brunner et al., 2004), we believe that $\tilde{c}$ represents the maximum speed at which information can be transmitted in the biological space-time. However, it might be possible that a space which satisfies the above constraints can not be built and hence, our theory would not apply.

On the other hand, given our ability to determine the value of $\tilde{c}$ and to represent dynamics along a trajectory on a manifold, our theory can be falsified or verified by considering the following cases. The first case occurs if the particle on the manifold moves with a constant velocity v. In this case, no accelerations are present and hence our theory based on accelerations does not apply. The second case takes place if the particle has an acceleration but its velocity $v<<\tilde{c}$. The particle is in a non-relativistic regime and hence, relativistic corrections are not needed. On the other hand, if accelerations are present and if $v~\tilde{c}$, we are in the relativistic regime and hence, relativistic corrections must be taken into account if we want to properly describe the dynamics on the manifold. We note that we construct this theory with the intention of using it to better understand biological aging. We therefore postulate that a $\tilde{c}$ may exist—or be assumed—on the same scale of observed velocities from human aging data so that we may use the theory to better study the data.

Once the submanifolds are constructed, then it is always possible to combine them considering their Cartesian product defined by Equation (4). Given our ability to construct these submanifolds, the comparison between the dynamics along different trajectories can be done if and only if they do belong to the same submanifold. So for example, if we are interested in comparing different rates of aging of two individuals with respect to the methylation status, we will need to construct the methylation status manifold, evaluate the trajectories of their methylation status in time, and analyze the dynamics on each of them. In particular, the comparison between the accelerations along both trajectories may help in understanding the difference between the rate of aging of the two individuals. In this sense, the importance is the different dynamics along each trajectory and not the fact that there are two different trajectories.

Other groups have proposed a framework for biological time to explain biological rhythms and other oscillating or biological processes that repeat periodically (Bailly et al., 2011). Our biological space-time model is more general, and is aimed towards explaining non-linear biological phenomena that are not necessarily periodic or oscillating, although these can be considered special cases of more generalized trajectories. In complement to Noble's hypothesis that there exists no privileged level of causation in biology (Noble, 2012), we suggest that the multiple scales of biology may be interpreted as trajectories on distinct manifolds that may be combined and coupled to each other. Although generalizations of relativity and space-time have been proposed, here we adapt these mathematical structures in order to interpret trajectories on manifolds to represent biological processes (O'Neill, 1983). We want to note that, as described in section 2.1.1, there can be submanifolds in which the velocity and the acceleration are identically zero. However, this does not imply an absence of motion on the other submanifolds (the other circle $S_2^1$) or on the principal manifold (the torus). In this case, it will be impossible to draw any conclusion about the possible difference between chronological and biological time for the process that takes place in the manifold in which the velocity, and hence the acceleration, is zero.

A novel concept which naturally follows from our approach is the notion of force in a biological context. A biological force may be a generalization of a physical force in a biological context which changes the dynamics of a particle on a manifold. We may hypothesize then, that accelerated dynamics on the manifold may be the effect of biological forces acting on the particle. For example, cancer and chemotherapy may be interpreted as forces which combine and have a resultant which may change the dynamics of the point on a manifold with a consequent impact on the rate of aging.

The shape of the manifolds which characterize biological space-time also affect our notion of distance, which requires the knowledge of the metric tensor of the manifold. Although a manifold is defined to be locally Euclidean, two biological processes or objects that are sufficiently different from each other (in space-time) may not be measured by using the conventional definition of Euclidean distance. The degree to which we must redefine our notion of a metric in the biological space-time may also affect our ability to interpret long-time dynamical behaviors of biological systems and whether or not relativistic corrections are required in order to properly analyze these dynamics. We therefore hypothesize that the underlying geometry of the manifold, or equivalently, the elements and structure of the metric tensor, may provide valuable insights into the biological process.

6. Summary

In summary, we have proposed a mathematical framework and criteria that may be used to define and construct a biological space-time. We use this as a tool to model and study different rates of biological aging based on the concept of relativistic corrections of time due to the acceleration or deceleration of biological dynamics relative to a common frame of reference. We discuss some examples of biological processes that illustrate these concepts, provide criteria that may be used to test or falsify the theory, and discuss implications and novel hypotheses that are generated by this model.

Author Contributions

RR and DM: Conception, design, and approval of final manuscript. All authors: Manuscript writing and editing.

Funding

Research reported in this publication was supported by the National Cancer Institute of the National Institutes of Health under award number P30CA033572, R01CA205247, and R01CA178387. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Leo D. Wang for a critical reading of the manuscript, and Jacob G. Scott for discussions that inspired and motivated this work. DM would like to thank G. Bocca for useful discussions.

References

Alfano, C. M., Peng, J., Andridge, R. R., Lindgren, M. E., Povoski, S. P., Lipari, A. M., et al. (2017). Inflammatory cytokines and comorbidity development in breast cancer survivors versus noncancer controls: evidence for accelerated aging? J. Clin. Oncol. 35, 149–156. doi: 10.1200/JCO.2016.67.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

Arora, M., Sun, C.L., and Ness, K.K. (2016). Physiologic frailty in nonelderly hematopoietic cell transplantation patients: results from the bone marrow transplant survivor study. JAMA Oncol. 2, 1277–1286. doi: 10.1001/jamaoncol.2016.0855

PubMed Abstract | CrossRef Full Text | Google Scholar

Bailly, F., Longo, G., and Montevil, M. (2011). A 2-dimensional geometry for biological time. Prog. Biophys. Mol. Biol. 106, 474–484. doi: 10.1016/j.pbiomolbio.2011.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Bocklandt, S., Lin, W., Sehl, M. E., Sanchez, F. J., Sinsheimer, J. S., Horvath, S., et al. (2011). Epigenetic predictor of age. PLoS ONE 6:e14821. doi: 10.1371/journal.pone.0014821

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandon, R. N. (1997). Does biology have laws? The experimental evidence. Philos. Sci. 64, S444–S457. doi: 10.1086/392621

CrossRef Full Text | Google Scholar

Brunner, N., Scarani, V., Wegmüller, M., Legré, M., and Gisin, N. (2004). Direct measurement of superluminal group velocity and signal velocity in an optical fiber. Phys. Rev. Lett. 93:203902. doi: 10.1103/PhysRevLett.93.203902

PubMed Abstract | CrossRef Full Text | Google Scholar

Curtius, K., Wong, C.-J., Hazelton, W. D., Kaz, A. M., Chak, A., Willis, J. E., et al. (2016). A molecular clock infers heterogeneous tissue age among patients with Barrett's esophagus. PLoS Comput. Biol. 12:e1004919. doi: 10.1371/journal.pcbi.1004919

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, L., Yu, W. M., Zheng, H., Loh, M. L., Bunting, S. T., Pauly, M., et al. (2016). Leukaemogenic effects of Ptpn11 activating mutations in the stem cell microenvironment. Nature 539, 304–308. doi: 10.1038/nature20131

PubMed Abstract | CrossRef Full Text | Google Scholar

Fougère, B., Boulanger, E., Nourhashémi, F., Guyonnet, S., and Cesari, M. (2017). Chronic inflammation: accelerator of biological aging. J. Gerontol. Series A Biol. Sci. Med. Sci. 72, 1218–1225. doi: 10.1093/gerona/glw240

PubMed Abstract | CrossRef Full Text | Google Scholar

Franceschi, C., Capri, M., Monti, D., Giunta, S., Olivieri, F., Sevini, F., et al. (2007). Inflammaging and anti-inflammaging: a systemic perspective on aging and longevity emerged from studies in humans. Mech. Ageing Dev. 128, 92–105. doi: 10.1016/j.mad.2006.11.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Fried, L. P., Tangen, C. M., Walston, J., Newman, A. B., Hirsch, C., Gottdiener, J., et al. (2001). Frailty in older adults evidence for a phenotype. J. Gerontol. Series A 56, M146–M157. doi: 10.1093/gerona/56.3.M146

PubMed Abstract | CrossRef Full Text | Google Scholar

Frieden, B. R., and Gatenby, R. A. (2011). Information dynamics in living systems: prokaryotes, eukaryotes, and cancer. PLoS ONE 6:e22085. doi: 10.1371/journal.pone.0022085

PubMed Abstract | CrossRef Full Text | Google Scholar

Gatenby, R. A., and Frieden, B. R. (2007). Information theory in living systems, methods, applications, and challenges. Bull. Math. Biol. 69, 635–657. doi: 10.1007/s11538-006-9141-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Henderson, T. O., Ness, K. K., and Cohen, H. J. (2014). Accelerated aging among cancer survivors: from pediatrics to geriatrics. Am. Soc. Clin. Oncol. Educ. Book. 2014:e423. doi: 10.14694/EdBook_AM.2014.34.e423

CrossRef Full Text | Google Scholar

Hogan, D. B., MacKnight, C., and Bergman, H. (2003). Models, definitions, and criteria of frailty. Aging Clin. Exp. Res. 15(3 Suppl.), 1–29. doi: 10.1016/B978-012369391-4/50051-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Horvath, S. (2013). DNA methylation age of human tissues and cell types. Genome Biol. 14:R115. doi: 10.1186/gb-2013-14-10-r115

PubMed Abstract | CrossRef Full Text | Google Scholar

Hurria, A., Jones, L., and Muss, H. B. (2016). Cancer treatment as an accelerated aging process: assessment, biomarkers, and interventions. Am. Soc. Clin. Oncol. Educ. Book 35, e516–e522. doi: 10.14694/EDBK_156160

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenkinson, G., Pujadas, E., Goutsias, J., and Feinberg, A. P. (2017). Potential energy landscapes reveal the information-theoretic nature of the epigenome. Nat. Genet. 49:059766. doi: 10.1038/ng.3811

CrossRef Full Text | Google Scholar

Klerman, E. B., and Hilaire, M. S. (2007). Review: on mathematical modeling of circadian rhythms, performance, and alertness. J. Biol. Rhythms 22, 91–102. doi: 10.1177/0748730407299200

PubMed Abstract | CrossRef Full Text | Google Scholar

Koelwyn, G. J., Khouri, M., Mackey, J. R., Douglas, P. S., and Jones, L. W. (2012). Running on empty: cardiovascular reserve capacity and late effects of therapy in cancer survivorship. J. Clin. Oncol. 30, 4458–4461. doi: 10.1200/JCO.2012.44.0891

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohanski, R. A., Deeks, S. G., Gravekamp, C., Halter, J. B., High, K., Hurria, A., et al. (2016). Reverse geroscience: how does exposure to early diseases accelerate the age-related decline in health? Ann. N. Y. Acad. Sci. 1386, 30–44. doi: 10.1111/nyas.13297

PubMed Abstract | CrossRef Full Text | Google Scholar

Kovtonyuk, L. V., Fritsch, K., Feng, X., Manz, M. G., and Takizawa, H. (2016). Inflamm-aging of hematopoiesis, hematopoietic stem cells, and the bone marrow microenvironment. Front. Immunol. 7:502. doi: 10.3389/fimmu.2016.00502

PubMed Abstract | CrossRef Full Text | Google Scholar

Levine, M. E. (2013). Modeling the rate of senescence: can estimated biological age predict mortality more accurately than chronological age? J. Gerontol. Series A Biol. Sci. Med. Sci. 68, 667–674. doi: 10.1093/gerona/gls233

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy-Leblond, J.-M. (1976). One more derivation of the Lorentz transformation. Am. J. Phys. 44, 271–277. doi: 10.1119/1.10490

CrossRef Full Text | Google Scholar

Longo, G., and Montévil, M. (2014). Perspectives on Organisms. Berlin; Heidelberg: Springer.

Google Scholar

López-Otín, C., Blasco, M. A., Partridge, L., Serrano, M., and Kroemer, G. (2013). The hallmarks of aging. Cell 153, 1194–1217. doi: 10.1016/j.cell.2013.05.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Mantovani, A., Allavena, P., Sica, A., and Balkwill, F. (2008). Cancer-related inflammation. Nature 454, 436–444. doi: 10.1038/nature07205

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittal, M., Siddiqui, M. R., Tran, K., Reddy, S. P., and Malik, A. B. (2014). Reactive oxygen species in inflammation and tissue injury. Antiox. Redox Signal. 20, 1126–1167. doi: 10.1089/ars.2012.5149

PubMed Abstract | CrossRef Full Text | Google Scholar

Mojtahedi, M., Skupin, A., Zhou, J., Castaño, I. G., Leong-Quong, R. Y. Y., Chang, H., et al. (2016). Cell fate decision as high-dimensional critical state transition. PLoS Biol. 14:e2000640. doi: 10.1371/journal.pbio.2000640

PubMed Abstract | CrossRef Full Text | Google Scholar

Ness, K. K., Armstrong, G. T., Kundu, M., Wilson, C. L., Tchkonia, T., and Kirkland, J. L. (2015). Frailty in childhood cancer survivors. Cancer 121, 1540–1547. doi: 10.1002/cncr.29211

PubMed Abstract | CrossRef Full Text | Google Scholar

Ness, K. K., Krull, K. R., Jones, K. E., Mulrooney, D. A., Armstrong, G. T., Green, D. M., et al. (2013). Physiologic frailty as a sign of accelerated aging among adult survivors of childhood cancer: a report from the st jude lifetime cohort study. J. Clin. Oncol. 31, 4496–4503. doi: 10.1200/JCO.2013.52.2268

PubMed Abstract | CrossRef Full Text | Google Scholar

Nestorowa, S., Hamey, F. K., Pijuan Sala, B., Diamanti, E., Shepherd, M., Laurenti, E., et al. (2016). A single-cell resolution map of mouse hematopoietic stem and progenitor cell differentiation. Blood 128, 20–32. doi: 10.1182/blood-2016-05-716480

PubMed Abstract | CrossRef Full Text | Google Scholar

Noble, D. (2012). A theory of biological relativity: no privileged level of causation. Interface Focus 2, 55–64. doi: 10.1098/rsfs.2011.0067

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Neill, B. (1983). Semi-Riemannian Geometry with Applications to Relativity, 1st Edn. Los Angeles, CA: Academic Press.

Google Scholar

Rizvi, A. H., Camara, P. G., Kandror, E. K., Roberts, T. J., Schieren, I., Maniatis, T., and Rabadan, R. (2017). Single-cell topological RNA-seq analysis reveals insights into cellular differentiation and development. Nat. Biotechnol. 35, 551–560. doi: 10.1038/nbt.3854

PubMed Abstract | CrossRef Full Text | Google Scholar

Rockwood, K., and Mitnitski, A. (2011). Frailty defined by deficit accumulation and geriatric medicine defined by frailty. Clin. Geriatric Med. 27, 17 – 26. doi: 10.1016/j.cger.2010.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruse, M. E. (1970). Are there laws in biology? Aust. J. Philos. 48, 234–246. doi: 10.1080/00048407012341201

CrossRef Full Text | Google Scholar

Seals, D. R., Justice, J. N., and LaRocca, T. J. (2016). Physiological geroscience: targeting function to increase healthspan and achieve optimal longevity. J. Physiol. 594, 2001–2024. doi: 10.1113/jphysiol.2014.282665

PubMed Abstract | CrossRef Full Text | Google Scholar

Shiju, S., and Sriram, K. (2017). Hypothesis driven single cell dual oscillator mathematical model of circadian rhythms. PLoS ONE 12:e0177197. doi: 10.1371/journal.pone.0177197

PubMed Abstract | CrossRef Full Text | Google Scholar

Snir, S., VonHoldt, B. M., and Pellegrini, M. (2016). A statistical framework to identify deviation from time linearity in epigenetic aging. PLoS Comput. Biol. 12:e1005183. doi: 10.1371/journal.pcbi.1005183

PubMed Abstract | CrossRef Full Text | Google Scholar

Steensma, D. P., Bejar, R., Jaiswal, S., Lindsley, R. C., Sekeres, M. A., Hasserjian, R. P., et al. (2015). Clonal hematopoiesis of indeterminate potential and its distinction from myelodysplastic syndromes. Blood 126, 9–16. doi: 10.1182/blood-2015-03-631747

PubMed Abstract | CrossRef Full Text | Google Scholar

Tu, L. W. (2011). An Introduction to Manifolds. New York, NY: Springer.

van Breukelen, F., and Martin, S. L. (2015). The hibernation continuum: physiological and molecular aspects of metabolic plasticity in mammals. Physiology 30, 273–281. doi: 10.1152/physiol.00010.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, G. P. (2017). The power of negative [theoretical] results. Proc. Natl. Acad. Sci. U.S.A. 114, 12851–12852. doi: 10.1073/pnas.1718862114

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhang, H., Rodriguez, S., Cao, L., Parish, J., Mumaw, C., et al. (2014). Notch-dependent repression of miR-155 in the bone marrow niche regulates hematopoiesis in an NF-κB-dependent manner. Cell Stem Cell 15, 51–65. doi: 10.1016/j.stem.2014.04.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Appendix A

By differentiating Equation (15) we obtain the transformations for the velocity

$\begin{array}{ll}\begin{array}{l}{u^{\prime }}_x=\frac{d{x}^{\prime }}{d{t}^{\prime }}=\frac{dx-vdt}{dt-vdx/{\tilde{c}}^2}=\frac{u_x-v}{1-v{u}_x/{\tilde{c}}^2},\hfill \\ {u^{\prime }}_y=\frac{d{y}^{\prime }}{d{t}^{\prime }}=\frac{u_y}{1-v{u}_x/{\tilde{c}}^2},\hfill \\ {u^{\prime }}_z=\frac{d{z}^{\prime }}{d{t}^{\prime }}=\frac{u_z}{1-v{u}_x/{\tilde{c}}^2}.\hfill \end{array}& (\mathrm{A.1})\end{array}$

where u_x = dx/dt, u_x = dx/dt and u_x = dx/dt are the velocities along the corresponding direction. An additional differentiation leads to the transformations of the acceleration

$\begin{array}{ll}\begin{array}{l}{a^{\prime }}_x=\frac{d{u^{\prime }}_x}{d{t}^{\prime }}=\frac{d{u^{\prime }}_x}{dt}\frac{dt}{d{t}^{\prime }}=\frac{a_x\left(1-v{u}_x/{\tilde{c}}^2\right)+(u_x-v)\left(v{a}_x/{\tilde{c}}^2\right)}{\gamma {(1-v{u}_x/{\tilde{c}}^2)}^3}\hfill \\ =\frac{a_x(1-v^2/c^2)}{\gamma {(1-v{u}_x/{\tilde{c}}^2)}^3}=\frac{a_x}{{\gamma }^3{(1-v{u}_x/{\tilde{c}}^2)}^3}\hfill \\ {a^{\prime }}_y=\frac{a_y}{{\gamma }^2{(1-v{u}_x/{\tilde{c}}^2)}^2}-\frac{u_{y}v{a}_x/{\tilde{c}}^2}{\gamma {(1-v{u}_x/{\tilde{c}}^2)}^3},\hfill \\ {a^{\prime }}_z=\frac{a_z}{{\gamma }^2{(1-v{u}_x/{\tilde{c}}^2)}^2}-\frac{u_{y}v{a}_x/{\tilde{c}}^2}{\gamma {(1-v{u}_x/{\tilde{c}}^2)}^3},\hfill \end{array}& (\mathrm{A.2})\end{array}$

where we have evaluated dt′/dt by considering the last of Equation (15) and we used the definition of the Lorentz factor γ. The sets of Equation (A.1) and (A.2) provide the connection between the velocity and the acceleration of an object measured by an observer in S and another observer in S′. We now consider the particular case in which a point is moving along a straight line oriented along the x-axis with a constant acceleration $a_x^{\prime }$ (Figure 4 on the right). We assume S′ to be located at the position of the particle and therefore, its velocity is the same as the velocity of the particle and it increases in time. Hence, by setting v = u_x the x-component of the acceleration $a_x^{\prime }$ becomes

$\begin{array}{ll}{a^{\prime }}_x=\frac{a_x}{{\gamma }^3{(1-u_x^2/{\tilde{c}}^2)}^3}=\frac{a_x{(1-u_x^2/{\tilde{c}}^2)}^{3/2}}{{(1-u_x^2/{\tilde{c}}^2)}^3}=\frac{a_x}{{(1-u_x^2/{\tilde{c}}^2)}^{3/2}},& (\mathrm{A.3})\end{array}$

and the acceleration in S is given by

$\begin{array}{ll}{a}_x={a^{\prime }}_x{(1-u_x^2/{\tilde{c}}^2)}^{3/2}.& (\mathrm{A.4})\end{array}$

By integrating the above equation with respect to time and by imposing v(0) = 0 we obtain the velocity

$\begin{array}{ll}{u}_x=\frac{{a^{\prime }}_{x}t}{\sqrt{1+{a^{\prime }}_x{}^2{t}^2/{\tilde{c}}^2}}.& (\mathrm{A.5})\end{array}$

By integrating once again with respect to time and by setting $x(0)={\tilde{c}}^2/a_x^{\prime }$ we obtain the position

$\begin{array}{ll}x(t)=\frac{{\tilde{c}}^2}{{a^{\prime }}_x}\sqrt{1+\frac{{a^{\prime }}_x^2{t}^2}{{\tilde{c}}^2}}.& (\mathrm{A.6})\end{array}$

The time t′ of the accelerated particle, and hence in the accelerated frame of reference, is given by integrating Equation (16)

$\begin{array}{ll}\begin{array}{l}{t}^{\prime }={\displaystyle {\int }_0^t\sqrt{1-u_x^2(s)/{\tilde{c}}^2}}ds={\displaystyle {\int }_0^t\frac{ds}{\sqrt{1+{({a^{\prime }}_{x}s/\tilde{c})}^2}}}\hfill \\ =\frac{\tilde{c}}{{a^{\prime }}_x}\text{arcsinh}\left(\frac{{a^{\prime }}_{x}t}{\tilde{c}}\right).\hfill \end{array}& (\mathrm{A.7})\end{array}$