Analysis of Cell Signal Transduction Based on Kullback–Leibler Divergence: Channel Capacity and Conservation of Its Production Rate during Cascade

Abstract

Kullback–Leibler divergence (KLD) is a type of extended mutual entropy, which is used as a measure of information gain when transferring from a prior distribution to a posterior distribution. In this study, KLD is applied to the thermodynamic analysis of cell signal transduction cascade and serves an alternative to mutual entropy. When KLD is minimized, the divergence is given by the ratio of the prior selection probability of the signaling molecule to the posterior selection probability. Moreover, the information gain during the entire channel is shown to be adequately described by average KLD production rate. Thus, this approach provides a framework for the quantitative analysis of signal transduction. Moreover, the proposed approach can identify an effective cascade for a signaling network.

1. Introduction

Kullback–Leibler divergence (KLD) is a type of generalized entropy or information quantity. It was introduced by Solomon Kullback and Richard A. Leibler, who discussed information source coding theory for information transmission efficiency [1]. At present, KLD finds diverse applications, including the imaging analytical field [2,3], hydrodynamics [4], clinical laboratory tests including electrocardiogram [5,6], network analysis [7,8] for biological applications [9], cellular biology [10], evaluating the bioequivalence of formulations of a drug [11], and experimental design for clinical study [12]. In this study, KLD is applied to analyze the cell information transmission, signal transduction, mediated by the cellular biochemical reaction. In particular, the proposed approach using Bayesian statistics [13,14], which is based on KLD, is expected to provide a novel theoretical framework [7].

Previous studies have discussed cell signal transduction through pathways from a viewpoint of similarity in the thermodynamic process that produces entropy [15]. Luo et al. [16] analyzed the heat production during carbohydrate metabolism and estimated the relationship between energy consumption and biological information from a biological metabolism perspective. Moreover, several research studies have applied variable concepts of entropy. In biologic genome informatics, a type of expanded Shannon entropy, such as the local Shannon-Jayne entropy, is utilized for analyzing the correlation between a set of gene expressions [17,18,19,20,21,22]. Teschendorff et al. [20] introduced a stochastic matrix, whose components are the normalized probabilities of the gene expressions in individual samples; the signal entropy rate was obtained using the matrix. In the network, the maximum entropy rate is determined by the adjacency matrix of the network [20]. For a multi-cell system, the network approach is useful for understanding biological signal transduction behaviors [23,24,25]. Recently, single cell entropy was introduced as the analytical basis of the variable phenotypic or genotypic state of a single cell, which is based on the assembly framework of statistical physics [26]. In another study, using a non-Markovian approach, the mutual entropy between the stimulus and the response in a biological system was considered in a sensory system, wherein the past trajectories are utilized to add useful information to the present state. In a simulation conducted by Becker et al., E. coli was used to reliably predict the concentration changes of environmental chemokines for chemotaxis [27].

Previously, the authors reported analyses of biological signal transduction based on information thermodynamics. In these studies, a theoretical framework was developed on Shannon entropy. However, the objective of the present study is to evaluate signal transduction efficiency of KLD in individual steps on the actual biochemical reaction kinetics [28,29]. KLD is introduced for analysis of signal transduction in reference to fluctuation theorem (FT) [30,31,32,33].

Signal Cascade Model

The signal events can be modeled as a cascade of modification and/or demodification cycle reactions of proteins in a cell that are named signaling molecules. Equation (1) presents a signal cascade model [34,35]. Here, suffixes m and j represent the number of cascades and the step number, respectively. In this model, the signaling molecule at step 1 of cascade m, denoted by X_m1, induces the modification of the X_m2 into X_m2* by binding the signal mediated molecule A such as adenosine triphosphate (ATP). Subsequently, X_m2 activates X_m3 in the same manner. In this way, the signaling molecule at the (j − 1)-th step of cascade m, denoted as X_mj−1, induces the modification of X_mj into X_mj*. As the opposite orientation of signal, demodification of X_mj* into X_mj occurs, at the −(j − 1)-th step of cascade m, and the pre-stimulation steady state is subsequently recovered [34]:

$$\begin{array}{l}{X}_{m1}+X_{m2}+A\leftrightarrow X_{m1}+X_{m2}*\hspace{0.17em}\hspace{0.17em}:1^{\mathrm{st}}\&-1^{\mathrm{st}}\mathrm{step}\\ \\ X_{mj-1}*+X_{mj}+A\leftrightarrow X_{mj-1}*+X_{mj}* :j-1^{\mathrm{th}}\&-{(j-1)}^{\mathrm{th}}\mathrm{step}\\ \\ X_{m,n-1}*+X_{mn}+A\leftrightarrow X_{m,n-1}*+X_{mn}*:n-1^{\mathrm{th}}\&-{(n-1)}^{\mathrm{th}}\mathrm{step}\end{array}$$

(1)

In the above model, the subscript m represents the total number of the cascade.

We introduce a priori (prior) selection probability of signaling molecule for the analysis. Here, q_mj, which represents the selection probability of inactive X_mj used in the j-th step in cascade m (forward direction), takes the form of the j-th molecule. On the other hand, q_mj*, which represents the selection probability of active X_mj*, is used in the −j-th step for cascade m (backward direction), as follows:

$$\begin{array}{l}{q}_{mj}=X_{mj}^0/X_m\\ q_{mj}*=X_{mj}^0*/X_m\end{array}$$

(2)

where,

$${\displaystyle \sum _{j=1}^n\left(q_{mj}+q_{mj}*\right)}=1$$

(3)

Here, X_m indicates the total concentration of signaling molecules in cascade m:

$$X_m={\displaystyle \sum _{j=1}^n\left(X_{mj}^0+X_{mj}^0*\right)}$$

(4)

and the total concentration of active and inactive signaling molecules is given by:

$${\displaystyle \sum _{j=1}^n{X}_{mj}^0}=X_m$$

(5)

and

$${\displaystyle \sum _{j=1}^n{X}_{mj}^0}*=X_m*$$

(6)

The total duration of cascade m, τ_m⁰, which indicates the sum of forward and backward cascades comprising a set of signaling molecules, is determined by:

$${\tau }_m^0={\displaystyle \sum _{j=1}^n\left(X_m^0{\tau }_{mj}^0-X_m^0*{\tau }_{-mj}^0\right)}$$

(7)

In Equations (2), (5), (6) and (7), the total duration was determined using the probabilities q_mj and q_mj*:

$${\tau }_m^0=X_m{\displaystyle \sum _{j=1}^n\left(q_{mj}{\tau }_{mj}^0-q_{mj}*{\tau }_{-mj}^0\right)}$$

(8)

The suffix 0 represents the prior state. Here, the duration, as forward τ_mj⁰ and backward τ_−mj⁰, are defined as shown in Figure 1. The Positive and negative values are assigned to τ_mj⁰ and τ_−mj⁰ corresponding to the direction of the step in the m cascade [34,35]. In the above equations, τ_mj⁰ represents the duration corresponding to positive code length in which the active molecule X_mj* increases in concentration. On the other hand, τ_−mj⁰ represents the duration corresponding to negative code length in which the active molecule X_mj* decreases in concentration. In this manner, the duration of individual step j-th can be represented as τ_mj⁰ − τ_−mj⁰.

2. Results

2.1. A Prior Probability Distribution of Signaling Molecules

Here, the author hypothesizes that the selection of signaling molecules is equal a priori. In our previous studies [34,35], Shannon’s entropy H_m for m-th cascade was demonstrated. Using Equations (3), (5) and (6), entropy H_m⁰ at a priori (prior) state can be represented as:

$$H_m^0=-X_m\left({\displaystyle \sum _{j=1}^n{q}_{mj}\log q_{mj}+}{\displaystyle \sum _{j=1}^n{q}_{mj}*\log q_{mj}*}\right)$$

(9)

To maximize H_m⁰, using non-determined parameters α_m⁰, and β_m⁰, in reference to the constraints established by Equations (3) and (8), let us introduce a function G.

$$\begin{array}{l}G\left(q_{m1},q_{m2},\cdots q_{mn};q_{m1}*,q_{m2}*,\cdots q_{mn}*;X_m\right)\\ =H_m^0-{\alpha }_m^0{\displaystyle \sum _{j=1}^n\left(q_{mj}+q_{mj}*\right)}-{\beta }_m^0{\tau }_m^0\\ =H_m^0-{\alpha }_m^0{\displaystyle \sum _{j=1}^n\left(q_{mj}+q_{mj}*\right)}-{\beta }_m^0{X}_m{\displaystyle \sum _{j=1}^n(q_{mj}{\tau }_{mj}^0-q_{mj}*{\tau }_{-mj}^0})\end{array}$$

(10)

Then, we have

$$\frac{\partial G}{\partial q_{mj}}=-X_m\left(\log q_{mj}-{\beta }_m^0{\tau }_{mj}^0\right)-{\alpha }_m^0-X_m$$

(11)

$$\frac{\partial G}{\partial q_{mj}*}=-X_m\left(\log q_{mj}+{\beta }_m^0{\tau }_{-mj}^0\right)-{\alpha }_m^0-X_m$$

(12)

$$\frac{\partial G}{\partial X_m}=-\left({\displaystyle \sum _{j=1}^n{q}_{mj}\log q_{mj}+q_{mj}*\log q_{mj}*}\right)-{\beta }_m^0{\displaystyle \sum _{j=1}^n\left(q_{mj}{\tau }_{mj}^0-q_{mj}*{\tau }_{-mj}^0\right)}$$

(13)

To maximize G, the right sides of Equations (11)–(13) are equated to zero.

$$-\log q_{mj}={\beta }_m^0{\tau }_{mj}^0({\tau }_{mj}^0>0)$$

(14)

$$-\log q_{mj}*=-{\beta }_m^0{\tau }_{-mj}^0({\tau }_{-mj}^0<0)$$

(15)

As indicated above, Equations (14) and (15) imply an important result; the coefficient β_m⁰ is independent of the step number j. Therefore, Equations (14) and (15) will be utilized as a prior probability distribution later. Therefore, from Equations (9), (14) and (15) the author has:

$$H_m^0=X_m{\beta }_m^0{\tau }_m^0$$

(16)

2.2. Average Entropy Production Rate in a Signal Cascade

Next, the kinetics were investigated using q_m (j|j − 1), which is the transitional probability of the j-th given (j − 1)-th step, and v_m (j|j − 1), which is the transitional rate of the j-th step in a forward signaling direction. Given j-th step, q_m (j − 1|j) is the transitional probability of the (j − 1)-th step given step j-th step. Similarly, given j-th step, v_m (j − 1|j) is the transitional rate of the (j − 1)-th step in a backward signaling direction in a given cascade. The cell system remains at detailed balance around the steady state, the homeostasis, as follows:

$$q_m\left(j|j-1\right)v_m^0\left(j|j-1\right)=q_m\left(j-1|j\right)v_m^0\left(j-1|j\right)$$

(17)

Therefore, the author has:

$$\log \frac{q_m\left(j-1|j\right)}{q_m\left(j|j-1\right)}=\log \frac{v_m^0\left(j|j-1\right)}{v_m^0\left(j-1|j\right)}$$

(18)

Using kinetic coefficients for (j − 1)-th step and the reverse −(j − 1)-th step, k_m,j−1 and k_m,−(j−1) in (1), the right side of (18) is given

$$\log \frac{q_m\left(j-1|j\right)}{q_m\left(j|j-1\right)}=\log \frac{v_m^0\left(j|j-1\right)}{v_m^0\left(j-1|j\right)}=\log \frac{k_{m,j-1}{X}_{m,j-1}^{0 *}{X}_{mj}^{0}A}{k_{m,-(j-1)}{X}_{m,j-1}^{0 *}{X}_{mj}^0*}=\log \frac{k_{m,j-1}{X}_{mj}^{0}A}{k_{m,-(j-1)}{X}_{mj}^0*}$$

(19)

When the change of X_m,j−1* is negligible during signal transduction, relative to the fluctuation of X_mj*, we have:

$$\log \frac{q_m\left(j-1|j\right)}{q_m\left(j|j-1\right)}=\log \frac{k_{m,j-1}{q}_{mj}A}{k_{m,-(j-1)}{q}_{mj}*}$$

(20)

In above, we used Equation (2). Dividing the both sides of above by τ_mj − τ_−mj and taking the limit, the variables q_mj and q_mj* remain in the right side:

$$\begin{array}{l}\underset{{\tau }_{mj}^0-{\tau }_{-mj}^0\to \infty }{\lim }\frac{1}{{\tau }_{mj}^0-{\tau }_{-mj}^0}\log \frac{q_m\left(j-1|j\right)}{q_m\left(j|j-1\right)}\\ =\underset{{\tau }_{mj}^0-{\tau }_{-mj}^0\to \infty }{\lim }\frac{1}{{\tau }_{mj}^0-{\tau }_{-mj}^0}\log \frac{q_{mj}}{q_{mj}*}\end{array}$$

(21)

Using Equations (14) and (15), the author has:

$$\underset{{\tau }_{mj}^0-{\tau }_{-mj}^0\to \infty }{\lim }\frac{1}{{\tau }_{mj}^0-{\tau }_{-mj}^0}\log \frac{q_m\left(j-1|j\right)}{q_m\left(j|j-1\right)}=-{\beta }_m^0\left(\frac{{\tau }_{mj}^0+{\tau }_{-mj}^0}{{\tau }_{mj}^0-{\tau }_{-mj}^0}\right)~-{\beta }_m^0$$

(22)

$$\underset{\left|-{\tau }_{mj}^0+{\tau }_{-mj}^0\right|->\infty }{\lim }\frac{1}{{\tau }_{mj}^0-{\tau }_{-mj}^0}\log \frac{q_m\left(j|j-1\right)}{q_m\left(j-1|j\right)}={\beta }_m^0\left(\frac{{\tau }_{mj}^0+{\tau }_{-mj}^0}{\left|-{\tau }_{mj}^0+{\tau }_{-mj}^0\right|}\right)~{\beta }_m^0$$

(23)

In above Equations (22) and (23), |τ_−mj⁰| is sufficiently longer than τ_mj⁰, according to experimental studies (Figure 1) [36,37,38,39,40,41,42,43,44,45,46,]. Here, using an arbitrary time parameter t, the average entropy production rate (AEPR), 〈ζ_mi〉 and 〈ζ_−mi〉 are defined during signal transduction for τ_mj⁰ − τ_−mj⁰ and |τ_−mj⁰ − τ_mj⁰|. respectively for m cascade and reverse cascade −m.

$$\langle {\zeta }_{mj}\rangle \triangleq \frac{1}{{\tau }_{mj}^0-{\tau }_{-mj}^0}{\displaystyle {\int }_0^{{\tau }_{mj}^0-{\tau }_{-mj}^0}{\zeta }_{mj}\left(t_{mj}\right)}d{t}_{mj}$$

(24)

$$\langle {\zeta }_{-mj}\rangle \triangleq \frac{1}{\left|{\tau }_{mj}^0-{\tau }_{-mj}^0\right|}{\displaystyle {\int }_0^{\left|{\tau }_{mj}^0-{\tau }_{-mj}^0\right|}{\zeta }_{-mj}\left(t_{-mj}\right)}d{t}_{-mj}$$

(25)

The fluctuation theorem (FT) states that the right sides of Equations (22)–(25) are equal to AEPR

$$\underset{{\tau }_{mj}^0-{\tau }_{-mj}^0\to \infty }{\lim }\frac{1}{{\tau }_{mj}^0-{\tau }_{-mj}^0}\log \frac{q_m\left(j|j-1\right)}{q_m\left(j-1|j\right)}=\langle {\zeta }_{mj}\rangle$$

(26)

$$\underset{\left|-{\tau }_{mj}^0+{\tau }_{-mj}^0\right|->\infty }{\lim }\frac{1}{\left|{\tau }_{mj}-{\tau }_{-mj}\right|}\log \frac{q_m\left(j-1|j\right)}{q_m\left(j|j-1\right)}=\langle {\zeta }_{-mj}\rangle$$

(27)

and

$${\beta }_m^0=-\langle {\zeta }_{mj}\rangle =\langle {\zeta }_{-mj}\rangle$$

(28)

Here, β_m⁰ has the dimension of entropy production rate and AEPRs are independent of the step number. Subsequently, AEPRs are redefined using Equations (14), as follows:

$$-\langle {\zeta }_{-mj}\rangle \triangleq \langle {\zeta }_m\rangle =-\frac{\log q_{mj}}{{\tau }_{mj}^0}$$

(29)

Notably, Equations (14), (15) and (28) indicate that AEPR is consistent during signal cascade. Here, the channel capacity is given by AEPR.

$$-\log q_{mj}=\langle {\zeta }_m\rangle {\tau }_{mj}^0$$

(30)

$$\log q_{mj}*=\langle {\zeta }_m\rangle {\tau }_{-mj}^0$$

(31)

In previous our studies, a simple formulation was proposed between selection probability q_mj and duration τ_mj, using an arbitrary parameter, ζ, which was independent of step numbers [34,35].

2.3. Multinomial Distribution with Population Distribution

KLD was used as a measure of information gain when obtaining a posteriori (posterior) distribution from the prior distribution in Equations (30) and (31) to a posterior distribution. Let probability q_mj be the prior distribution. Therefore, the uncertainty reduces:

$$\mathrm{Uncertainty}=-\log p_{mj}-\left(-\log q_{mj}\right)$$

(32)

We define information D_m, as the KLD of signal events in cascade m,

$$D_m\left(p_m\Vert q_m\right)=-X_m{\displaystyle \sum _{j=1}^n\left(p_{mj}\log \frac{p_{mj}}{q_{mj}}+p_{mj}*\log \frac{p_{mj}*}{q_{mj}*}\right)}$$

(33)

The above equation represents KLD, which indicates the average value of the information obtained from data I with respect to p_mj and p_mj*. Consequently, KLD is known as information gain. The maximum likelihood estimation is thought to be an estimation method that empirically minimizes KLD.

Posterior probabilities are defined:

$${\displaystyle \sum _{j=1}^n\left[p_{mj}+p_{mj}*\right]}=1$$

(34)

and

$$\begin{array}{l}{p}_{mj}=X_{mj}/X_m\\ p_{mj}*=X_{mj}*/X_m\end{array}$$

(35)

In addition,

$${\tau }_m =X_m{\displaystyle \sum _{j=1}^n\left(q_{mj}{\tau }_{mj}-q_{mj}*{\tau }_{-mj}\right)}$$

(36)

Therefore, when considering the signal transduction occurs under a certain given condition, the probability is transformed from q_mj* into p_mj* with minimum KLD under the given condition.

To minimize D_m (p_mj||q_mj) using non-determined parameters α_m, and β_m in reference to the constraints established by Equations (34) and (36), a function L was introduced to apply Lagrange’s method to undetermined multipliers.

$$\begin{array}{l}L\left(p_{m1},p_{m2},\cdots p_{mn};p_{m1}*,p_{m2}*,\cdots p_{mn}*;X_m\right)\\ =D_m\left(p_m\Vert q_m\right)-{\alpha }_m{\displaystyle \sum _{j=1}^n\left(p_{mj}+p_{mj}*\right)}-{\beta }_m{\tau }_m\end{array}$$

(37)

Here, the differences between α_m and α_m⁰ and β_m⁰ and β_m are indicative of the signaling. Subsequently,

$$\frac{\partial L}{\partial p_{mj}}=-X_m\left(\log \frac{p_{mj}}{q_{mj}}+{\beta }_m{\tau }_{mj}\right)-{\alpha }_m-X_m$$

(38)

$$\frac{\partial L}{\partial p_{mj}*}=-X_m\left(\log \frac{p_{mj}*}{q_{mj}*}-{\beta }_m{\tau }_{-mj}\right)-{\alpha }_m-X_m$$

(39)

$$\frac{\partial L}{\partial X_m}=-\left({\displaystyle \sum _{j=1}^n{p}_{mj}\log \frac{p_{mj}}{q_{mj}}}+{\displaystyle \sum _{j=1}^n{p}_{mj}*\log \frac{p_{mj}*}{q_{mj}*}}\right)-{\beta }_m\left({\displaystyle \sum _{j=1}^n{\tau }_{mj}{p}_{mj}}-{\displaystyle \sum _{j=1}^n{\tau }_{-mj}{p}_{mj}*}\right)$$

(40)

For the minimization of L, the right hand sides of Equations (38)–(40) are equated to zero, as follows:

$$-\log \frac{q_{mj}}{p_{mj}}={\beta }_m{\tau }_{mj}({\tau }_{mj}>0)$$

(41)

$$-\log \frac{q_{mj}*}{p_{mj}*}=-{\beta }_m{\tau }_{-mj}({\tau }_{-mj}<0)$$

(42)

Therefore, from (41) and (42), the author has

$$-\log p_{mj}={\beta }_m{\tau }_{mj}+{\beta }_m^0{\tau }_{mj}$$

(43)

$$-\log p_{mj}*=-{\beta }_m{\tau }_{-mj}-{\beta }_m^0{\tau }_{-mj}$$

(44)

$${\alpha }_m=-X_m$$

(45)

Accordingly, from Equations (7), (43) and (44), KLD is given by:

$$\begin{array}{l}{D}_m\left(p_m\Vert q_m\right)=-X_m{\displaystyle \sum _{j=1}^n\left(p_{mj}\log \frac{p_{mj}}{q_{mj}}+p_{mj}*\log \frac{p_{mj}*}{q_{mj}*}\right)}\\ =X_m{\beta }_m{\displaystyle \sum _{j=1}^n\left(p_{mj}{\tau }_{mj}-p_{mj}*{\tau }_{-mj}\right)}\\ ={\beta }_m{\tau }_m\end{array}$$

(46)

And the author has:

$${\beta }_m=\frac{D_m\left(p_m\Vert q_m\right)}{{\tau }_m}={\delta }_m\left(p_m\Vert q_m\right)$$

(47)

Accordingly, β_m is equal to the average KLD production rate, δ(p_m‖q_m), during the signal transduction and is consistent during the entire signal cascade. Therefore,

$$-\log \frac{p_{mj}}{q_{mj}}=\frac{D_m\left(p_m\Vert q_m\right)}{{\tau }_m}{\tau }_{mj}({\tau }_{mj}>0)$$

(48)

$$-\log \frac{p_{mj}*}{q_{mj}*}=-\frac{D_m\left(p_m\Vert q_m\right)}{{\tau }_m}{\tau }_{-mj}({\tau }_{-mj}<0)$$

(49)

The author defined p_m(j|j − 1), which is the transitional probability of the j-th given (j − 1)-th step, and v_m(j|j − 1), which is the transitional rate of the j-th step in a forward signaling direction. In addition, given j-th step, p_m(j − 1|j) is the transitional probability of the (j − 1)-th step given step j-th step. Similarly, given the j-th step, v_m(j − 1|j) is the transitional rate of the (j − 1)-th step in a backward signaling direction in a given cascade.

$$p_m\left(j|j-1\right)v_m\left(j|j-1\right)=p_m\left(j-1|j\right)v_m\left(j-1|j\right)$$

(50)

Likewise from (18)–(20), dividing the both sides of above by τ_mj − τ_−mj and taking the limit, we have:

$$\begin{array}{l}\underset{{\tau }_{mj}-{\tau }_{-mj}\to \infty }{\lim }\frac{1}{{\tau }_{mj}-{\tau }_{-mj}}\log \frac{p_m\left(j-1|j\right)}{p_m\left(j|j-1\right)}\\ =\underset{{\tau }_{mj}-{\tau }_{-mj}\to \infty }{\lim }\frac{1}{{\tau }_{mj}-{\tau }_{-mj}}\log \frac{p_{mj}}{p_{mj}*}\end{array}$$

(51)

Further, Equations (47), (50) and (51) give

$$\begin{array}{l}\underset{{\tau }_{mj}-{\tau }_{-mj}\to \infty }{\lim }\frac{1}{{\tau }_{mj}-{\tau }_{-mj}}\log \frac{p_m\left(j-1|j\right)/q_m\left(j-1|j\right)}{p_m\left(j|j-1\right)/q_m\left(j|j-1\right)}\\ =\underset{{\tau }_{mj}-{\tau }_{-mj}\to \infty }{\lim }\frac{1}{{\tau }_{mj}-{\tau }_{-mj}}\frac{p_{mj}}{p_{mj}*}\frac{q_{mj}}{q_{mj}*}\\ ~{\delta }_m\left(p_m\Vert q_m\right)\end{array}$$

(52)

In above, the author used τ_mj << τ_−mj (Figure 1). Above Equation corresponds to the extended FT considering KLD, and the extended channel capacity can be defined as:

$$C_m=\underset{\tau \to \infty }{\lim }\frac{K{D}_m\left(p_m\Vert q_m\right)}{{\tau }_m}=\underset{\tau \to \infty }{\lim }K{\delta }_m\left(p_m\Vert q_m\right)$$

(53)

Here, K is an arbitrary constant. If entropy unit is used, K = k_B, Boltzmann’s constant. On the other hand, in information science, K is equivalent to log₂e.

3. Conclusions

Recently, theoretical analysis of the transduction capacity of biochemical signaling networks has greatly developed [47]. In this study, KLD, the average KLD production rate, and channel capacity based on the average KLD production rate were shown to be critical quantities that can be attributed to the entire signal transduction step. This KLD is a tool to estimate entropy production from stationary trajectories [48].

In this work, the author deduced simple but important relational formulae, (47)–(53). For a prior distribution probability q_mj, the author derived Equations (14) and (15) in association with FT and source-coding theory [28]. This method was introduced in our previous study of Tsallis entropy [29]. The theoretical framework in the current study is shown in Figure 2.

In this study, uniform distribution was not applied as a prior distribution q_mj before stimulus. As reported previously in [34,35], the selection probability of signal molecules, q_mj, can be described by simple formulae given by Equations (14) and (15), which illustrate the simple relationship between the logarithm of the probability and time elapsed between tentative modification and demodification. Subsequent to the stimulus, the selection probability of the signal molecules are transformed into a posterior distribution probability p_mj. It is likely that the formulation of signal transformation using KLD is more intuitive and easier to understand than using mutual entropy.

The Mitogen-activated Protein Kinase (MAPK) pathway is a multistep signal conversion step, in which the process of modification and demodification in the entire cascade can be understood as a repeated cycle reaction. It has been experimentally demonstrated that the demodification process is significantly longer than the modification process, as shown in τ_mj << τ_−mj (Figure 1), with the former requiring a few hours for completion, while the latter achieves completion in a few minutes. This asymmetry in the time course kinetics points to an important result pertaining the conservation of AEPR and average KLD production rate with reference to FT. In addition, it should be noted that in Equation (53), the channel capacity of the entire signal cascade is represented by KLD. Moreover, the author introduced average KLD production rate in Equation (47), and the production rate was found to be consistent during the whole cascade, which is an expanded form of AEPR conservation during the whole cascade, as reported previously by the author [35].

In the experimental studies [36,37,38,39,40,41,42,43,44,45,46], the ratio of modified active type of signaling molecules was determined by the immunoblot intensity or other corresponding data. Thus, it is possible to compute, in principle, the channel capacity of the entire signal cascade. This can be attributed to the fact that in most cases, extracellular substances, such as ligands, simultaneously promote the activation of multiple cascades. Therefore, the selection of specific ligands for the given cascade is essential and enables the measurement of the rigorous activity of the cascade. In future, such recombinant protein ligand will be required to quantify the signal cascade accurately.

Thus, KLD and average KLD production rate may be regarded as a critical attribution of cell signal cascade. The thermodynamic approach in this manuscript can provide a theoretical framework for the quantitative analysis of signal transduction.