Next Article in Journal
When Social Networks Meet D2D Communications: A Survey
Previous Article in Journal
Cortical Network Response to Acupuncture and the Effect of the Hegu Point: An fNIRS Study
Previous Article in Special Issue
General Architecture for Development of Virtual Coaches for Healthy Habits Monitoring and Encouragement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chain Modeling of Molecular Communications for Body Area Network

1
School of Communication and Information Engineering, Chongqing University of Posts and Telecommunica-Tions, Chongqing 400065, China
2
Key Laboratory of Optical Communication and Networks in Chongqing, Chongqing 400065, China
3
Key Laboratory of Ubiquitous Sensing and Networking in Chongqing, Chongqing 400065, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(2), 395; https://doi.org/10.3390/s19020395
Submission received: 30 October 2018 / Revised: 26 December 2018 / Accepted: 15 January 2019 / Published: 18 January 2019
(This article belongs to the Special Issue Wireless Body Area Networks and Connected Health)

Abstract

:
Molecular communications provide an attractive opportunity to precisely regulate biological signaling in nano-medicine applications of body area networks. In this paper, we utilize molecular communication tools to interpret how neural signals are generated in response to external stimuli. First, we propose a chain model of molecular communication system by considering three types of biological signaling through different communication media. Second, communication models of hormonal signaling, Ca 2 + signaling and neural signaling are developed based on existing knowledge. Third, an amplify-and-forward relaying mechanism is proposed to connect different types of signaling. Simulation results demonstrate that the proposed communication system facilitates the information exchange between the neural system and nano-machines, and suggests that proper adjustment can optimize the communication system performance.

1. Introduction

Rapid progress of nano-technology enables the manufacturing of nano-machines for medical applications of body area network (BAN) [1]. Thus far, a bio-inspired manufacturing method has attracted considerable attention. It assembles biological molecules or atoms into a micro- to nano-scale unit based on the functional structure of nature cells. To finish a complex task, the cooperation of nano-machines is indispensable due to their limited size and power. Molecular communication (MC) is a promising technology for the cooperation of nano-machines, exploiting biological molecules as intra-body information carriers [2]. Design of MC is inspired by mechanisms of biological activities in nature (e.g., hormone communication between ants and neural signaling inside human body), which can be promoted by the experience from diverse traditional communication systems [3,4]. Promising MC applications include health monitoring systems, drug delivery systems and bio-sensor networks [5].
MC could be divided into wireless and wired types based on different properties of biological mechanisms. Existing research efforts generally focus on a single MC type and lack the comprehensive understanding of the complex biological activities in human bodies. Biologically speaking, the human body is a complex system composed of various biological substances including organs, tissues, cells, molecules, etc. It is widely accepted that, even for a simple biological activity, different substances interact with each other through exchange of power and substance. Moreover, processes of biological activities usually involve multiple types of biological signaling, which should be investigated jointly rather than independently. From the aspect of communication, proper integration of bio-systems promote the implementations of complex functionalities. For example, neurons and blood vessels form neurovascular unit, which performs more complex functions via communication pathways, including maintaining the normal activities of neural system, repairing the damaged neurons with the nutrition from the blood, and regulating the vasodilatation of blood vessels [6].
In this paper, we develop a chain model of MC system, based on biological cells, and motivated by different chemical stimuli methods [7]. The model considers three types of MC signaling: hormonal signaling, Ca 2 + signaling and neural signaling. Hormonal signaling is based on the diffusion of hormones in fluid medium. Ca 2 + signaling is based on the oscillation of Ca 2 + ions in astrocytes. Neural signaling is based on the spike train transportation of action potential in neurons. The proposed model involves a chain of communications, in which hormonal signals are considered as the input, triggering oscillation of Ca 2 + signals. As the relay, Ca 2 + signals further induce neural signals as the output. The proposed chain model relies on different bio-sensors in biological cells. For example, hormones are observed by G-protein bio-sensors in astrocytes (see Section 5), and neural signals are collected by AMPA or NMDA bio-sensors in neurons (see Section 6).
This paper is developed based on our earlier work [8], where an initial hybrid model of MC is proposed including diffusive and neural channels. In that work, nano-machines are grouped to function as the sender, relay and receiver. The communication model in [8] includes high biological plurality and therefore is reliable. Novel contributions of this paper are listed as follows.
  • First, we develop a chain model of molecular communication based on biological signaling. The proposed model considers the biological interactions among hormone, Ca 2 + and neural signals, which is more general than the model of [8].
  • Second, we propose an implementable amplify-and-forward relaying mechanism instead of decode-and-forward relaying as in [8]. In addition, multiple astrocytes are utilized instead of a nano-machine, elevating the reliability of relaying from Ca 2 + signaling to neural signaling.
  • Third, based on the work in [8], we examined the relations between communication performance and more parameters of the proposed model. We also found that source coding is efficient in improving the communication performance, which may provide a guidance for nano-machine design.

2. Related Work

We examined the literature on MC related to this work. Diffusion-based MC has attracted great attention due to its universality in biology. In [9], a mathematical model of communication is developed based on Brownian motion, and the channel capacity is theoretically deduced. In [10], a joint optimization strategy is developed to improve the performance of the two-hop diffusion-based MC system. In [11], different constant composition codes are designed for channel detection at the receiver without channel state information.
Another potential type of MC is Ca 2 + signaling. In [12], a linear channel model for intra/inter-cellular Ca 2 + MC system is developed, where the channel gain and communication delay are analyzed. In [13], Ca 2 + signaling models for different types of cells are introduced and compared. In addition, various communication performance indicators are deduced for these models. In [14], a channel switching functionality is designed using Ca 2 + signaling in three biological cells, and different gating properties of gap junction channels are considered.
Neural signaling has been abundantly studied based on the knowledge of neural science. What concerns us most is the communication model in neural networks, and there are innumerable studies focusing on this area. In [15], a multi-input-single-output model of synaptic channel in hippocampus area of the brain is developed, and the communication performance is analyzed using Shannon theory. In [16], the capacity of the axonal channel is derived, where parameters are optimized to improve the performance. In [17], a fast algorithm is proposed for performance analysis of artificial and biological neural systems.

3. Overview of the Communication System

The proposed MC system is illustrated in Figure 1. The proposed communication system is motivated by the experience of multiple traditional communication systems [18,19,20,21,22,23]. Specifically, the communication system is modeled as a chain, which is composed of a sender, a receiver and three types of communication media. The sender and receiver are nano-machines with different functions, respectively, denoted by Tx and Rx. Tx is assumed to release hormones in response to external stimuli. Tx can be, for example, engineered endocrine cells [2], which are able to release hormones under control. Rx is a bio-inspired or electronic nanomachine that is able to detect neural signals.
The first type of communication medium is the fluid, where hormones propagate around the astrocytes. There are two types of propagation mechanisms, determined by types and properties of the hormones and fluid medium. First, hormones perform approximative random walk towards any direction with low energy consumption in the passive diffusion mechanism. For example, steroid hormones propagate in the brain via lipid mediator and passive diffusion [24]. Second, hormones drift towards specific directions with high energy consumption in active diffusion. For example, thyroid hormones permeate blood–brain barrier by a slow flux rate in bloodstream via active diffusion [25].
The second type of communication medium is astrocytes, which are specific cells connected with neurons. Astrocytes are closely related to neural activities, such as physically supporting structures of neurons, providing nutrition for neurons, and promoting the growth of neurons [26]. Astrocytes can promote the generation of neural signals through adjustment of Ca 2 + oscillation, thereby serving as nature bridges between another biological system and neural system. Reportedly, astrocytes can respond to hormones (e.g., norepinephrine) in generation of Ca 2 + signals [27]. Based on this knowledge, astrocytes function as the perfect relay in our proposed model, connecting hormonal signaling with neural signaling.
The third type of communication medium is neuron, which processes and transmits information through the human body. Neural communication is an efficient and reliable biological signaling, exploiting both electrical and chemical methods. Different types of neurons respond to different external stimuli (e.g., electricity, chemical and light) [28].
As shown in Figure 1, the communication system is developed to interpret the neural signaling evoked by chemical stimuli. We assume that the proposed biological communication system transmits information by encoding the presence or absence of signals in time scale. Moreover, we assume that time-slotted mechanism works for the proposed communication system, based on the biological clocks guaranteeing the synchronization of biological cells [29].
The communication process is described as follows. Initially, Tx releases hormones in a controlled pattern, responding to external stimuli. Then, hormones diffuse in the fluid medium, some of which are absorbed by the astrocytes, inducing the oscillation of Ca 2 + signaling. Oscillated Ca 2 + ions in astrocytes flow into the neurons and trigger neural signaling. Finally, the neural signals are detected and decoded by Rx.

4. Hormonal Signaling Model

In Tx, a list of binary information is encoded with the number of released hormone molecules per time slot. One bit is transmitted in a symbol period, denoted by T. It is found that secretion pattern of hormones is usually with burst and not sustaining [30]. Accordingly, we apply a simple coding for MC, i.e., on–off key (OOK) with square pulses. During a time of T, a certain number of hormone molecules Q 0 are released into fluid medium, if Tx transmits bit “1”; no hormone molecules are released into fluid medium, if Tx transmits bit “0”. Let R H ( t ) be the released molecules during the ith slot, the expression of OOK coding is given by,
R H ( t ) = Q 0 δ ( t i T ) bit = 1 0 bit = 0
where δ ( t ) is the Dirac’s delta function of t. Hormone molecules propagate with random or directed walks, motivated by the passive or active diffusion mentioned in Section 3. Plenty of factors impact the distance and velocity of diffusion processes, including molecules diameter, medium type, temperature and pressure. For simplicity, we only consider passive diffusion of hormones, i.e., hormones performing random walks. Concentration of the hormone molecules presents a decreasing gradient with the distance, which could be expressed by the well-known Fick’s second law,
d C H ( x , t ) d t = D H 2 C H ( x , t ) + R H ( t )
where 2 is the Laplace operator, indicating the sum of the second derivatives of C H ( x , t ) . R H ( t ) is the releasing rate of molecules at Tx in time scale. D H is the diffusion coefficient of hormones, indicating the diffusion rate. Let x S be the position of Tx, and x A denote the position of an astrocyte. By solving Equation (2), the concentration accumulated at the position of x A is given by,
C H ( x A , t ) = R H ( t ) · g ( d S , A , t ) = R H ( t ) ( 4 π D H t ) 3 2 e x p ( d S , A 2 4 D H t )
where g ( x , t ) is the Green function of hormones [9], expressing the concentration distribution of hormone molecules in space and time. d S , A denotes the distance between Tx and an astrocyte, given by d S , A = | x A x S | . A single hormonal molecule could be released by Tx in the n S th slot and absorbed by an astrocyte in the n A th slot. The corresponding probability is given by [31],
u n = n T ( n + 1 ) T g ( d S , A , t ) t m H ( v ) d v d t
where n = n A n S , and m H ( v ) denotes the probability function distribution of the molecular life expectancy, given by m H ( v ) = ι e ι v with mean of 1 / ι .

5. Calcium Signaling Model

Astrocytes are special star-shaped glial cells, physically close to neurons, and play a crucial role in regulation of neural activities. It is observed that astrocytes propagate Ca 2 + waves in response to external stimuli, which further induce and impact the behaviors of neural signaling [32]. In this section, we explain how Ca 2 + signals are generated and propagate in astrocytes.
We introduce a mathematical model of Ca 2 + signaling in astrocytes according to experimental observations [27]. The dynamic of Ca 2 + signaling in an astrocyte is shown in Figure 2. One major organelle that stores Ca 2 + ions is endoplasmic reticulum (ER). With an external stimulus in astrocyte, Inositol 1,4,5-trisphosphate (IP 3 ) molecules are generated by G-protein coupled receptors (GR), changing the permeability of IP 3 -sensitive Ca 2 + receptors (IP 3 R) of the ER and induce the Ca 2 + fluxes. Three major Ca 2 + fluxes are considered, namely channel fluxes released from the ER to the cytoplasm via receptors of IP 3 ( U c h a n ), uptake fluxes from the cytoplasm to the ER by Ca 2 + pumps ( U p u m p ), and leak fluxes from from the ER to the cytoplasm ( U l e a k ). Induced Ca 2 + fluxes incur the oscillation of Ca 2 + , and patterns of Ca 2 + waves are positively related to the concentration of external stimuli. Let C C Y and C E R denote the Ca 2 + concentrations in cytoplasm and ER. Variable h is employed to denote the Ca 2 + channel inactivation gate. On the time scale, the introduced model of Ca 2 + signaling is expressed as a set of equations,
d C C Y d t = a 1 ( U c h a n + U l e a k ) U p u m p
d C E R d t = U c h a n + U l e a k 1 a 1 · U p u m p
d h d t = a 3 h a 2
where,
U c h a n = v c h a n · ( C I P 3 C I P 3 + b 1 ) 3 ( C C Y C C Y + b 2 ) 3 h 3 ( C E R C C Y )
U l e a k = v l e a k · ( C E R C C Y )
U p u m p = v p u m p · C C Y 2 C C Y 2 + b 3 2
a 2 = 1 b 6 ( b f + C C Y )
a 3 = b f b f + C C Y
b f = b 5 C I P 3 + b 1 C I P 3 + b 4
In the above equations, a 1 denotes the ratio of ER to cytoplasm, a 2 denotes the time constant of the inactivation gate, and a 3 denotes the steady value of the inactivation gate. v c h a n , v l e a k and v p u m p denote the maximum flux rates of U c h a n , U l e a k and U p u m p , respectively. b 1 b 6 are constants of the model, and C I P 3 denotes the concentration of IP 3 .
On the spatial scale, the propagation of Ca 2 + ions in astrocytes is based on the diffusion of Ca 2 + [12]. ER is modeled as a point source of releasing Ca 2 + . Ca 2 + ions around the boundary of astrocytes can flow into the connected neurons. For any position x inside an astrocyte, the dynamic of the Ca 2 + concentration via diffusion can be expressed by,
d C C Y ( x , t ) d t = D C a Δ C C Y ( x , t ) + R C a ( t )
where D C a denotes diffusion coefficient of Ca 2 + smf R C a ( t ) denotes the variation rate of Ca 2 + concentration, determined by Equations (5)–(13). Equation (14) can be solved similarly to Equation (3).

6. Neural Signaling Model

Neural communication is a fast biological signaling, and its velocity varies due to different neuron types. Neural communication is also biologically reliable, as the neurons are able to eliminate errors in communications. In this section, we introduce a signaling model across two neighboring neurons, including the processes of neural firing, axonal transmission, gap junctional transmission, and postsynaptic response.
One neuron is normally composed of different number of dendrites, axons, somas, and axonal terminals. In a communication, neural signals are first generated by ions through neural firing process in soma, transmitted along the axon fiber in the pattern of electrical impulses. Then, the electrical impulses trigger the release of neurotransmitters from the vesicles located in axon terminals. The neurotransmitters propagate to the neighboring neurons, until they are absorbed by the dendrites. Finally, the postsynaptic response is induced that might trigger firing of neighboring neurons.

6.1. Amplify-and-Forward Neural Firing

The generation of neural signals is called neural firing, which is an all-or-none process. Various types of ions including Ca 2 + are released by astrocytes and flow into neurons. During the process, local potentials of neuron membrane are gradually elevated and aggregated into a global potential. If the global potential exceeds a threshold, neural signals are generated. Otherwise, it fails to generate neural signals since the global potential falls quickly.
Neural firing is a relaying process, in which information carrier varies from Ca 2 + to neural signals. In our model, we apply amplify-and-forward (AF) relaying in neural firing process. AF relaying is based on the concentration-mediated mechanism in astrocyte–neuron communications [33]. We assume that local neural signals are positively elevated according to the concentration of Ca 2 + without decoding.
To generate neural signals, there must be adequate Ca 2 + ions flowing into neurons. The number of astrocytes should be sufficient to generate the required intensities of Ca 2 + flows. The local potentials are elevated by Ca 2 + from connected astrocytes, which cooperate for achieving neural firing. This process is expressed by the well-known Hodgkin–Huxley model [34], given by,
V N ( t ) = V r + i M A t i ν ( s ) A N , i ( t s ) d s
where M A is the total number of astrocytes contributing to the neural firing, V N ( t ) denotes the global potential on the neural membrane during the firing process, V r is the resting global potential without any stimuli, ν is the linear response of the neural membrane to an input pulse of Ca 2 + flow, A N , i ( t ) denotes the intensity of Ca 2 + flow released by the ith astrocyte, injected into the neuron at time t i , and A N , i ( t ) is proportional to Ca 2 + concentration of astrocyte, expressed by A N , i ( t ) C C Y ( t ) . The action potential spikes are typical neural signal S ( t ) , generated when V N ( t ) is strong enough to exceed firing threshold θ 1 ,
S ( t ) = j φ ( t t j )
where φ ( t ) denotes the waveform of action potential spike, and t j denotes the time that the spike occurs. Based on experimental observations, S ( t ) follows or approaches Poisson distribution in response to random intensities of Ca 2 + flows [35]. We employ λ ( t ) to denote the firing rate, which is determined by not only the Ca 2 + flows, but also the inherent capability of neurons. The distribution of S ( t ) could be impacted by the noise or some other stochastic factors.

6.2. Axonal Transmission

In the process of axonal transmission, action potential spikes transmit along the axon until they arrive in the axon terminals. During the process, various ions (e.g., K + , N a + and C l ) exchange between the inside and outside of the axon fiber. The axonal transmission process could be modeled by the cable theory, expressed by the quantitative equation in [36],
γ m γ l d 2 V N ( x , t ) d x 2 = γ c γ m d V N ( x , t ) d t + V N ( x , t )
where γ m , γ l and γ c denote the membrane resistance, longitudinal resistance and capacitances, respectively. V N ( x , t ) denotes the potential of axonal fiber at time t and location x. The boundary condition of the axonal transmission is given by,
V N ( x , 0 ) = V 0 ( x ) V N ( 0 , t ) = 0
The above equations can be solved using the Fourier method [37], namely V N ( x , t ) is calculated as,
V N ( x , t ) = e x p ( t γ m γ c ) i = 1 ζ i e x p [ π 2 i 2 t L N 2 r l γ c ] c o s ( π · i · x L N )
where L N is the length of the axon fiber. ζ i denotes the Fourier coefficient, following the boundary conditions in Equation (18), and can be obtained by,
ζ i = 2 L N 0 L N V 0 ( x ) c o s ( π · i · x L N ) d x

6.3. Gap Junctional Transmission

6.3.1. Vesicle Release

We show the process of gap junctional transmission in Figure 3. Action potential spikes arrive in the axonal terminals, triggering vesicle movements and inducing vesicle-membrane fusion. As a result, neurotransmitters stored inside vesicles are released into the gap junction. The release process of neurotransmitters is stochastic, related to the vesicle-membrane distance [38]. Therefore, we exploit a pool model to classify vesicles into readily releasable pool (RRP) and reserve pool (RP), respectively. The vesicles of RP are difficult to release, since they are far from the membrane, whereas the vesicles of RRP are easy to release, since they are close to the membrane.
Let N R R P and N R P , respectively, denote the sizes of RRP and RP. The vesicles of RP are assumed to fill the RP with rate of τ f , if N R R P is reduced due to the release process. The average firing rate of vesicles is denoted by λ ¯ , which is calculated as the integral of firing rate λ ( t ) during a time slot,
λ ¯ = 0 T λ ( t ) d t T
The release probability of a vesicle is related to λ ¯ and N R R P , given by,
P r e l e = 1 e x p ( λ ¯ N R R P )
In local scale, P r e l e can be approximated as the linear relation with λ ¯ , i.e., P r e l e = λ ¯ N R R P . Due to the stochastic behavior of neurotransmitters, it is possible to cause errors during the vesicle release process. For binary information, the error probabilities of bit “1” and “0” are, respectively, calculated by,
P e _ n e u r o n 1 ( 1 ) = 1 ( P r e l e ) k i
P e _ n e u r o n 1 ( 0 ) = P r e l e k i
where k i is the number of action potential spikes during the ith slot. Equations (23) and (24) indicate that error transmission of binary information is caused by miss or redundant release of vesicles.

6.3.2. Queueing Model of Neurotransmitters

The neurotransmitters diffuse in the gap junctions until absorbed by the dendrite receptors of neighboring neurons. Due to the stochastic behavior of the vesicle release process, the absorbed time of neurotransmitters is random. The diffusion of neurotransmitters is slow; this process is normally ignored since gap junctions are extremely short (several nanometers). However, neurotransmitters can be absorbed with delay due to the limited absorbing ability of dendrites [39]. Generally, neurotransmitters gather around the receptors of dendrites and wait to be absorbed. To study this phenomenon, we use the queuing theory to calculate the delay of gap junctional transmission.
In our model, the neurotransmitters are defined as the batched customers, which are independent with identical distributions. We adopt the M x | G | 1 | K model, where M x is the absorbing time of these batched neurotransmitters. The time intervals of vesicle release follow the negative exponential distribution, i.e., f ( t i ) = σ λ b e x p ( σ λ b t i ) ( 0 < σ < 1 ) . μ b and τ b , respectively, denote the batch size and batch service rate, satisfying λ b = μ b τ b . K denotes the capacity of the receptor. G denotes absorbing time t s of the receptor with a distribution G ( t ) , which is given by 0 < t s = 1 μ = 0 t d G ( t ) < . The distribution of the batch sizes is expressed as ϵ ( x ) = ( 1 σ ) x 1 σ . Based on the result of [40], the blocking probability is given by,
P b l o = p T ( 0 ) τ b ( 1 σ ) 0 τ b ( 1 σ ) + σ + t = 1 N 1 p T ( t ) ( 1 σ ) ( N t ) + p T ( N )
where p T ( k ) is the probability that neurotransmitters are absorbed by k customers. The blocking phenomenon is more serious for bigger λ b . Based on the previous derivations, the average delay is calculated as,
ω = k = 1 N 1 r = i i = 1 N k p T ( t ) ϵ ( r ) 0 g 0 ( k + i 1 , t ) d t ϵ ( 1 P b l o ) + r = i i = 1 N p T ( 0 ) ϵ ( r ) 0 g 0 ( i 1 , t ) d t ϵ ( 1 P b l o )
where g 0 ( x , t ) is the probability that customers are served at time t.

6.4. Postsynaptic Response

The process of postsynaptic response includes the formation of postsynaptic potential and the decoding of neural signals.

6.4.1. Postsynaptic Potential

The neurotransmitters are absorbed by dendrites of neighboring neurons via various receptors, such as the AMPA and NMDA receptors. Neurotransmitters accumulate on the dendrite membrane, namely postsynaptic terminals. During this process, the postsynaptic potential of dendrite membrane is elevated to generate the waveform of output signals, which could be expressed by an alpha function [41],
ϝ ( t ) = G m a x t t p e x p ( 1 t t p )
where G m a x is the maximum amplitude of the potential, and t p is the corresponding time to peak.

6.4.2. Neural Decoding

Decoding process of the neural signal includes two steps: (1) estimation of the action potential spikes; and (2) extraction of the binary information.
In the first step, a waveform of output signal is denoted by y ( t ) , which could be calculated as a convolution of the action potential spikes and ϝ ( t ) ,
y ( t ) = q i δ ( t t i ) · ϝ ( t )
where q denotes a stochastic variable in the signal amplification, following the Gamma distribution [42]. To extract the binary information of y ( t ) , we first sample y ( t ) into discrete signal y ( n ) with interval τ s ,
y ( n ) = j = N s a m / 2 N s a m / 2 y ( n j τ s )
where N s a m is the number of samplings in waveform of y ( t ) , which depends on G m a x and t p . Let t p be the origin point; N s a m 2 sampling points are determined on both sides of the origin point, where k τ s < min ( t j t j 1 ) is satisfied. The sampling interval τ s should be properly regulated to alleviate the distortion of sampling signals. We adopt threshold ϑ 2 to determine whether a spike exists or not at the position of sampling. If n = 1 k y ( n ) k < ϑ 2 , a spike exists. Accordingly, we obtain a pattern of estimated spikes, denoted by y ( t ) .
In the second step, we extract the binary information by examining the distribution of y ( t ) . The Kolmogorov–Smirnov test [43] is employed to verify whether the distribution of y ( t ) follows the Poisson distribution in time slot T. We say y ( t ) follows the Poisson distribution if the following condition is satisfied,
i ϰ i > K α
where variable K α meets condition P ( K K α ) = 1 α , K denotes the Kolmogorov distribution, and ϰ i is the Kolmogorov–Smirnov statistic of the Poisson process during the ith slot. We obtain λ using the Maximum Likelihood Estimation (MLE). Threshold ϑ 3 is employed for decoding, which is smaller than λ . For the ith bit, if the condition in Equation (30) is satisfied and estimated λ is higher than ϑ 3 , the decoded information is bit “1”. Otherwise, the decoded information is bit “0”. Further, we deduce the transmission error probabilities of decoding for both bit “1” and “0”,
P e _ n e u r o n 2 ( 0 ) = P ( i ϰ i > K α ) P ( λ ϑ 3 )
P e _ n e u r o n 2 ( 1 ) = P ( i ϰ i K α ) + P ( i ϰ i > K α ) P ( λ < ϑ 3 )
According to Equations (23), (24), (31) and (32), the error probability of transmitting one bit neural signal is expressed as,
P n e u r o n ( 1 | 0 ) = P ( b i t = 0 ) P e _ n e u r o n 1 ( 0 ) ( 1 P e _ n e u r o n 2 ( 1 ) ) + p ( b i t = 0 ) ( 1 P e _ n e u r o n 1 ( 0 ) ) P e _ n e u r o n 2 ( 0 )
P n e u r o n ( 0 | 1 ) = P ( b i t = 1 ) P e _ n e u r o n 1 ( 1 ) ( 1 P e _ n e u r o n 2 ( 0 ) ) + p ( b i t = 0 ) ( 1 P e _ n e u r o n 1 ( 1 ) ) P e _ n e u r o n 2 ( 1 )

7. Channel Capacity and Transmission Delay

In this section, we deduce the expressions of channel capacity and transmission delay for the proposed communication system.

7.1. Channel Capacity

Let X i , Y i and Z i , respectively, represent the input, relay and output signal, indicating the diffusive concentration of hormone, oscillating concentration of Ca 2 + , and decoded action potential spikes. To deduce the channel capacity, we define the following equations,
P i ( 1 , 1 ) = P [ Y i = 0 | X i = 0 ]
P i ( 1 , 2 ) = P [ Y i = 1 | X i = 1 ]
P i ( 2 , 1 ) = P [ Y i = 0 | X i = 0 ] P [ Z i = 0 | Y i = 0 ]
P i ( 2 , 2 ) = P [ Y i = 1 | X i = 1 ] P [ Z i = 1 | Y i = 1 ]
P i ( 2 , 3 ) = P [ Y i = 0 | X i = 1 ] P [ Z i = 1 | Y i = 0 ]
P i ( 2 , 4 ) = P [ Y i = 1 | X i = 0 ] P [ Z i = 0 | Y i = 1 ]
In the above equations, P i ( 1 , 1 ) and P i ( 1 , 2 ) are related to the probability of transmitting bit “1”, denoted by p. According to the conclusions in [31], the expressions of P i ( 1 , 1 ) and P i ( 1 , 2 ) are, respectively, given by,
P i ( 1 ) = i = 1 N b 1 ( 1 p u i )
Q i ( 1 ) = 1 ( 1 u 0 ) i = 1 N b 1 ( 1 p u i )
where u i is the probability of Equation (4), and N b is the total number of time slots. The expressions of P [ Z i = 1 | Y i = 0 ] and P [ Z i = 0 | Y i = 1 ] can be calculated according to Equations (33) and (34). Based on information theory, the channel capacity is expressed as,
C = m a x i = 1 N b I ( X i ; Z i ) N b
We assume that, under the condition of Y, X is uncorrelated to Z. During the ith time slot, the mutual information is given by,
I ( X i ; Z i ) = I ( X i ; Y i Z i ) I ( X i ; Y i / Z i ) = I ( X i ; Y i ) I ( X i ; Y i / Z i )
The mutual information between X and Y is further calculated as,
I ( X i ; Y i ) = H ( Y i ) H ( Y i | X i ) = Ψ ( ( 1 p ) P i ( 1 , 1 ) + p ( 1 P i ( 1 , 2 ) ) ) p Ψ ( P i ( 1 , 2 ) ) ( 1 p ) Ψ ( P i ( 1 , 1 ) )
where the Ψ ( x ) = x l o g x ( 1 x ) l o g ( 1 x ) . The second item of Equation (44) is further expanded as,
I ( X i ; Y i / Z i ) = H ( X i | Z i ) H ( X i | Y i Z i ) = H ( Y i Z i | X i ) + H ( Y i Z i ) + H ( Z i | X i ) H ( Z i )

7.2. Delay

Transmission delay of the proposed communication system is denoted by Ω , which is calculated as the delay summation of hormonal signaling Ω H , Ca 2 + signaling Ω C a and neural signaling Ω N ,
Ω = Ω H + Ω C a + Ω N
The estimation of Ω H is related to the emission, propagation and reception processes of hormones. For simplicity, we only consider the propagation process in channel for the transmission delay calculation. The transmission delay is related to two factors on the scales of time and space. The first factor is the transmission distance on the spatial scale. One feasible method is to estimate based on the probability distribution function along propagation distance d S , A , given by,
f Ω H ( d S , A ) = 0 d S , A g ( x , t ) d x = e r f ( d S , A 2 4 D H Ω H ) 8 π D H Ω H
where e r f ( x ) is the error function, expressed as e r f ( x ) = 2 π 0 x e u 2 d u . Another factor determining transmission delay is the variation rate of molecules at the sender on the time scale. Due to the fixed number of molecules Q 0 in one time slot, the variation rate of molecules is determined by the length of each time slot. The average delay of hormonal diffusion is calculated as,
Ω H ( d S , A ) = 0 T Ω H f ( Ω H ) d Ω H
The estimation of Ω C a is based on the time difference, which starts with the generation of Ca 2 + waves in connected astrocytes, ends with the absorbing of Ca 2 + by a neuron. Similar with the estimation of Taynnan et al. [13], the variational concentrations of Ca 2 + in the astrocytes are assumed to be equal with the absorbed concentrations of Ca 2 + during an increasing delay of Ω C a [13], i.e.,
i = 1 M A 0 T Δ ( C C Y i ( t ) ) d t = 0 T + Ω C a Δ ( C N ( t ) ) d t
where Δ ( x ) denotes the variation of x during a time period, C C Y i ( t ) denotes the Ca 2 + concentration in cytoplasm of the ith connected astrocyte with the neuron, and C N ( t ) is the variational concentration of Ca 2 + absorbed by the neuron.
In the calculation of the neural signaling delay, we mainly consider the contribution of gap junctional transmission. Based on the proposed queuing model of gap junctional transmission, the average delay of the vesicle release process is approximated as the average waiting time ω of the queueing model. The delay of axonal transmission Ω A x is calculated based on Equation (19). Let N N be the number of neurons, and the average delay between two neighboring neurons is expressed as,
Ω N = N N ( ω + Ω A x )

8. Performance Evaluation

A computer simulation was employed to examine the performance of proposed communication system. The control variate method was exploited since there are too many parameters in the proposed system. We verified the relations between some important parameters and the performance indicators (i.e., channel capacity and transmission delay).

8.1. Simulation Design

Simulation parameters were grouped according to hormonal, Ca 2 + and neural signaling, respectively.
Parameters of the hormonal signaling are listed in [9,31,44]: diffusion coefficient D H was fixed at 4.8 μ m 2 /ms, number of released molecules per bit Q 0 varied from 1000 to 10,000, the distance between Tx and astrocytes were randomly generated within 1–20 μ m, the length of time slot T varied from 1 s to 20 s, molecule life expectancy ι was set to 0.2, and probability p of transmitting bit “1” varied between 0 and 1.
Parameters of the Ca 2 + signaling are listed in [27]: ratio of ER to cytoplasm a 1 was 0.185, the maximum released flux rate from ER to cytoplasm v c h a n was 6 s 1 , the maximum leaked flux rate from ER to cytoplasm v l e a k was 0.11 s 1 , the maximum pumped flux rate from cytoplasm to ER v p u m p was 0.9 μ Ms 1 , b 1 was 0.13 μ M, b 2 was 0.08234 μ M, b 3 was 0.1 μ M, b 4 was 0.9434 μ M, b 5 was 1.049 μ M, b 6 was 0.2 μ M 1 s 1 , diffusion coefficient D C a was 20 μ 2 s 1 , the average propagation distance of Ca 2 + signaling d E R , B was 10 μ m, and the length of time slot T varied from 1 to 20 s.
Parameters of neural signaling is listed in [15,36,38]: firing rate λ ¯ of neural signaling varied from 15 to 30 Hz, the average release probability of vesicles p v was set to 0.4, pool size of readily release pool N R R P varied from 3 to 10, the pool size of reserve pool N R P was set to 60, and the average filling time τ f varied from 10 to 40 ms. In axonal transmission, membrane resistance γ m was set as 64.1 M Ω , longitudinal resistance γ l was set as 8 M Ω , and capacitances γ c was set to 1 mF/cm 2 , the average length of axons L N was set to 43.2 μ m, and the length of time slot T varied from 1 to 20 s.

8.2. Behaviors of Signals

In the simulation, we present how various signals behave in 40 s to transmit 4 bits of binary information. Figure 4a,d,g presents hormonal signals, Figure 4b,e,h presents Ca 2 + signals, and Figure 4c,f,i presents neural signals. The results reveal that the hormonal and Ca 2 + signals propagate with low frequency, and the neural signals propagate with high frequency.
We analyzed the chain models by comparing a group of different signals. For example, the signals in Figure 4d,e,f transmit “1011”. Hormone concentrations generate three waves during 0–10 s and 20–40 s. Accordingly, oscillation of Ca 2 + concentration is triggered to generate periodic waves. However, oscillation frequency of 20–30 s is smaller than that of 30–40 s, because a part of hormones during 20–30 s still remain in 30–40 s (i.e., memory effect), increasing the hormonal concentration in 30–40 s and promoting the oscillation of Ca 2 + . Similarly, a faster oscillation of Ca 2 + generates more action potential spikes during 30–40 s, compared with those in 20–30 s.
We analyzed the behaviors of various signals by sending different bits. Comparing Figure 4a,d,g, we see that concentrations of hormones are elevated by the percentage of sending bit “1”. Similarly, the frequency of Ca 2 + oscillation is also related to percentage of sending bit “1”. In Figure 4h, there are totally four peaks, and the intervals between different peaks are decreasing, since the frequency of Ca 2 + oscillation is regulated. Moreover, Ca 2 + oscillation stops after 30 s, when hormonal concentration is too high, and the concentration of Ca 2 + stabilizes at a high level. In addition, the number of action potential spikes seems positively related to the oscillation frequency of Ca 2 + .

8.3. Channel Capacity

Channel capacity of the proposed communication system is jointly determined by the three types of signaling. Our target was to check the relations between channel parameters and the channel capacity.
Figure 5 shows how channel capacity C changes with different Q 0 , p and λ . Apparently, C is elevated by increasing number of released molecules Q 0 , because that astrocytes can easily absorb hormones due to the increase of Q 0 and the error probability of hormonal signaling is decreased. In addition, C reaches its peak if p is closer to 0.5. By comparing Figure 5a–c, we note that increase of firing rate λ of neural signaling restrains C, because a larger λ increases the error probability of vesicle release process.
Figure 6 shows how channel capacity C changes with different M A , p, N R R P and τ f . Obviously, the optimal p of maximizing C is not 0.5, which is different with traditional wireless communication. Thereby, the source coding of MC should be particularly designed. Increasing the number of astrocytes M A also elevates C, since the relaying capability of Ca 2 + signaling is enhanced. Moreover, Figure 6b reveals that C is larger with a larger readily release pool N R R P and a shorter filling time τ f .

8.4. Transmission Delay

We checked the relations between various channel parameters and the transmission delay. Figure 7a presents how transmission delay of hormonal signaling Ω H changes with different T and d S , A . Apparently, increasing the length of time slot T results in the growing Ω H , because releasing rate is decreased by increasing T and a fixed number of molecules Q 0 is released for each bit. The increasing rate of Ω H degrades fast when T is large. Therefore, we deduce that it is important to choose a proper T in the proposed model. Besides, we note that Ω H increases with growing distance d S , A between astrocyte and Tx, since the molecular concentration degrades significantly on the spatial scale.
Figure 7b presents how Ω C a changes with different T and M A . Similar with Ω H , the increased length of time slot T signifies the increased Ω C a . We guess that the variation of T impacts the hormonal signaling and Ca 2 + signaling. Moreover, a larger M A results in a smaller Ω C a , because the proposed system could transmit more Ca 2 + in a unit time, with the stronger relaying ability.
Figure 7c presents how the average waiting time ω in neural signaling changes with different N R R P , τ f and λ ¯ . The values of ω is much smaller than Ω C a and Ω H . We see that ω increases with firing rate λ , because firing rate promotes the vesicle release process, more neuro-transmitters are blocked in gap junctions, and ω increases. In addition, ω grows with the increasing size of readily release pool N R R P and the decreasing filling time τ f , since vesicle release process is promoted by those parameters.

9. Conclusions

In this paper, we design a chain model for molecular communication to understand and interpret neural signaling. The proposed model contains three types of signals, where Ca 2 + signals of astrocytes function as the relay of communication. Due to this advantage, the proposed model becomes more implementable. Behaviors of the different signals demonstrate their close relations with each other. Simulation results validate the effectiveness of the proposed chain model, and also reveal that proper source encoding helps improve the system performance.
Our work provides a philosophy of exploiting complex biological reactions for communication engineering in body area network. Future work will address the expansion of chain modeling and molecular communications in more complex networks, and explore more practical and complex biological activities.

Author Contributions

Conceptualization, P.H. and H.L.; methodology, P.H.; software, X.H.; validation, X.H.; formal analysis, P.H.; investigation, H.L.; resources, P.H. and H.L.; data curation, X.H.; writing—original draft preparation, P.H.; writing—review and editing, P.H. and H.L.; visualization, H.L.; supervision, P.H.; project administration, P.H.; funding acquisition, P.H.

Funding

This work was partially supported by the National Natural Science Foundation of China (61871062, 61471102, 61771082, and 61801065), and Program for Innovation Team Building at Institutions of Higher Education in Chongqing (CXTDX201601020).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cavallari, R.; Martelli, F.; Rosini, R.; Buratti, C.; Verdone, R. A survey on wireless body area networks: Technologies and design challenges. IEEE Commun. Surv. Tutor. 2014, 16, 1635–1657. [Google Scholar] [CrossRef]
  2. Nakano, T.; Eckford, A.W.; Haraguchi, T. Molecular Communication; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  3. Wu, D.; Yan, J.; Wang, H.; Wu, D.; Wang, R. Social attribute aware incentive mechanism for device-to-device video distribution. IEEE Trans. Multimed. 2017, 19, 1908–1920. [Google Scholar] [CrossRef]
  4. Wu, D.; Yan, J.; Wang, H.; Wu, D.; Wang, R. Dynamic trust relationships aware data privacy protection in mobile crowd-sensing. IEEE Internet Things J. 2018, 5, 2958–2970. [Google Scholar] [CrossRef]
  5. Qi, H.; Peng, Y.; Gao, Q.; Zhang, C. Applications of nanomaterials in electrogenerated chemiluminescence biosensors. Sensors 2009, 9, 674–695. [Google Scholar] [CrossRef] [PubMed]
  6. Lecrux, C.; Hamel, E. Nanonetworks: The neurovascular unit in brain function and disease. Acta Physiol. 2011, 203, 47–59. [Google Scholar] [CrossRef] [PubMed]
  7. Christian, C.A.; Mobley, J.L.; Moenter, S.M. Diurnal and estradiol-dependent changes in gonadotropin-releasing hormone neuron firing activity. Proc. Natl. Acad. Sci. USA 2005, 102, 15682–15687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. He, P.; Mao, Y.; Liu, Q.; Zhai, H. A diffusion-neuron hybrid channel for molecular communication. In Proceedings of the International Symposium on Intelligent Signal Processing and Communication Systems, Xiamen, China, 6–9 November 2017; pp. 855–860. [Google Scholar]
  9. Liu, Q.; Yang, K. Channel capacity analysis of a diffusion-based molecular communication system with ligand receptors. Int. J. Commun. Syst. 2015, 28, 1508–1520. [Google Scholar] [CrossRef]
  10. Satish, K.T.; Tadi, R.T.R.; Prabhat, K.U.; Daniel, B.D.C. Joint Optimization of Molecular Resource Allocation and Relay Positioning in Diffusive Nanonetworks. IEEE Access 2018, 66, 67681–67687. [Google Scholar]
  11. Jamali, V.; Ahmadzadeh, A.; Farsad, N.; Schober, R. Constant-composition codes for maximum likelihood detection without csi in diffusive molecular communications. IEEE Trans. Commun. 2018, 66, 1981–1995. [Google Scholar] [CrossRef]
  12. Bicen, A.O.; Akyildiz, I.F.; Balasubramaniam, S.; Koucheryavy, Y. Linear channel modeling and error analysis for intra/inter cellular Ca2+ molecular communication. IEEE Trans. NanoBiosci. 2016, 15, 488–498. [Google Scholar] [CrossRef]
  13. Taynnan, B.M.; Balasubramaniam, S.; Jennings, B. Comparative end-to-end analysis of Ca2+ signaling-based molecular communication in biological tissues. IEEE Trans. Commun. 2015, 63, 5128–5142. [Google Scholar] [CrossRef]
  14. He, P.; Nakano, T.; Mao, Y.; Lio, P.; Liu, Q.; Yang, K. Stochastic Channel Switching of Frequency-encoded Signals in Molecular Communication Networks. IEEE Commun. Lett. 2018, 22, 332–335. [Google Scholar] [CrossRef]
  15. Malak, D.; Akan, O. A Communication Theoretical Analysis of Synaptic Multiple-Access Channel in Hippocampal-Cortical Neurons. IEEE Trans. Commun. 2013, 61, 2457–2467. [Google Scholar] [CrossRef]
  16. Aghababaiyan, K.; Shah-Mansouri, V.; Maham, B. Axonal channel capacity in neuro-spike communication. IEEE Trans. NanoBiosci. 2018, 17, 78–87. [Google Scholar] [CrossRef] [PubMed]
  17. Bilgin, B.A.; Akan, O.B. A fast algorithm for analysis of molecular communication in artificial synapse. IEEE Trans. NanoBiosci. 2017, 16, 408–417. [Google Scholar] [CrossRef] [PubMed]
  18. Wu, D.; Liu, Q.; Wang, H.; Yang, Q.; Wang, R. Cache Less for More: Exploiting Cooperative Video Caching and Delivery in D2D Communications. IEEE Trans. Multimed. 2018. [Google Scholar] [CrossRef]
  19. Li, Y.; Liu, J.; Cao, B.; Wang, C. Joint Optimization of Radio and Virtual Machine Resources with Uncertain User Demands in Mobile Cloud Computing. IEEE Trans. Multimed. 2017, 20, 2427–2438. [Google Scholar] [CrossRef]
  20. Wu, D.; Zhang, F.; Wang, H.; Wang, R. Security-oriented opportunistic data forwarding in mobile social networks. Future Gener. Comput. Syst. 2017, 87, 803–815. [Google Scholar] [CrossRef]
  21. Li, Y.; Liang, Y.; Liu, Q.; Wang, H. Resources Allocation in Multi-Cell D2D Communications for Internet of Things. IEEE Internet Things J. 2018, 5, 4100–4108. [Google Scholar] [CrossRef]
  22. Wang, R.; Yan, J.; Wu, D.; Wang, H.; Yang, Q. Knowledge-centric edge computing based on virtualized d2d communication systems. IEEE Commun. Mag. 2018, 56, 32–38. [Google Scholar] [CrossRef]
  23. Wu, D.; Shi, H.; Wang, H.; Wang, R.; Fang, H. A Feature based Learning System for Internet of Things Applications. IEEE Internet Things J. 2018. [Google Scholar] [CrossRef]
  24. Pardridge, W.M.; Mietus, L.J. Transport of steroid hormones through the rat blood-brain barrier. Endocrinology 1979, 105, 605–612. [Google Scholar] [CrossRef] [PubMed]
  25. Villalobos, C.; Lucia, N.; Javier, G. Anterior pituitary thyrotropes are multifunctional cells. Am. J. Physiol. Endocrinol. Metab. 2004, 287, 1166–1170. [Google Scholar] [CrossRef] [PubMed]
  26. Reynolds, B.A.; Weiss, S. Generation of neurons and astrocytes from isolated cells of the adult mammalian central nervous system. Science 1992, 255, 1707–1710. [Google Scholar] [CrossRef] [PubMed]
  27. Roth, B.J.; Yagodin, S.V.; Holtzclaw, L.; Russell, J.T. A mathematical model of agonist-induced propagation of calcium waves in astrocytes. Cell Calcium 1995, 87, 53–64. [Google Scholar] [CrossRef]
  28. Flavio, F. Network Neuroscience; Elsevier Academic Press: Cambridge, MA, USA, 2016. [Google Scholar]
  29. Palmer, J.D.; Edmunds, L.N. Cellular and molecular bases of biological clocks. Bioscience 1988, 38, 161–162. [Google Scholar] [CrossRef]
  30. Kolesnick, R.N.; Gershengorn, M.C. Direct evidence that burst but not sustained secretion of prolactin stimulated by thyrotropin-releasing hormone is dependent on elevation of cytoplasmic calcium. J. Biol. Chem. 1985, 260, 5217–5220. [Google Scholar]
  31. Nakano, T.; Okaie, Y.; Liu, J.Q. Channel Model and Capacity Analysis of Molecular Communication with Brownian Motion. IEEE Commun. Lett. 2012, 16, 797–800. [Google Scholar] [CrossRef]
  32. Di, G.A.; Barbi, M.; Chillemi, S.; Alloisio, S.; Nobile, M. Calcium signalling in astrocytes and modulation of neural activity. Biosystems 2007, 89, 74–83. [Google Scholar]
  33. Navarrete, M.; Araque, A. Endocannabinoids mediate neuron-astrocyte communication. Neuron 2008, 57, 883–893. [Google Scholar] [CrossRef]
  34. Hodgkin, A.L.; Huxley, A.F. A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 1952, 117, 500–544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Pfister, J.P.; Toyoizumi, T.; Barber, D.; Gerstner, W. Optimal spike-timing-dependent plasticity for precise action potential firing in supervised learning. Neural Comput. 2006, 18, 1318–1348. [Google Scholar] [CrossRef] [PubMed]
  36. Rall, W. Optimal Cable theory for dendritic neurons. In Methods in Neuronal Modeling; MIT Press: Cambridge, MA, USA, 1989; pp. 9–92. [Google Scholar]
  37. Sreekumaran, E.; Ramakrishna, T.; Madhav, T.R.; Anandh, D.; Prabhu, B.M.; Sulekha, S.; Bindu, P.N.; Raju, T.R. Loss of dendritic connectivity in CA1, CA2, and CA3 neurons in hippocampus in rat under aluminum toxicity: Antidotal effect of pyridoxine. Brain Res. Bull. 2003, 59, 421–427. [Google Scholar] [CrossRef]
  38. Matveev, V.; Wang, X.J. Implications of all-or-none synaptic transmission and short-term depression beyond vesicle depletion: A computational study. Neuroscience 2000, 24, 1575–1588. [Google Scholar] [CrossRef]
  39. Pyle, J.L.; Kavalali, E.T.; Piedras-Rentería, E.S.; Tsien, R.W. Rapid reuse of readily releasable pool vesicles at hippocampal synapses. Neuron 2000, 28, 221–231. [Google Scholar] [CrossRef]
  40. Li, W.; Fretwell, R.J.; Kouvatsos, D.D. Performance Analysis of Queues with Batch Poisson Arrival and GE Service. In Proceedings of the 2011 7th International Conference on Wireless Communications, Networking and Mobile Computing, Wuhan, China, 23–25 September 2011; pp. 1033–1036. [Google Scholar]
  41. Rall, W. Distinguishing theoretical synaptic potentials computed for different soma-dendritic distributions of synaptic input. Neurophysiology 1967, 30, 1138–1168. [Google Scholar] [CrossRef] [PubMed]
  42. Bekkers, J.M. Origin of variability in quantal size in cultured hippocampal neurons and hippocampal slices. Neurobiology 1967, 87, 5359–5362. [Google Scholar] [CrossRef]
  43. Massey, J.F.J. The Kolmogorov-Smirnov test for goodness of fit. J. Am. Stat. Assoc. 1951, 46, 68–78. [Google Scholar] [CrossRef]
  44. Bogazzi, F.; Bartalena, L.; Brogioni, S.; Burelli, A.; Manetti, L.; Tanda, M.L.; Gasperi, M.; Martino, E. Thyroid vascularity and blood flow are not dependent on serum thyroid hormone levels: Studies in vivo by color flow Doppler sonography. Eur. J. Endocrinol. 1999, 140, 452–456. [Google Scholar] [CrossRef]
Figure 1. Proposed chain of MC system containing hormonal, Ca 2 + and neural signaling.
Figure 1. Proposed chain of MC system containing hormonal, Ca 2 + and neural signaling.
Sensors 19 00395 g001
Figure 2. Dynamic of Ca 2 + signals in an astrocyte.
Figure 2. Dynamic of Ca 2 + signals in an astrocyte.
Sensors 19 00395 g002
Figure 3. Diagram of gap junctional transmission.
Figure 3. Diagram of gap junctional transmission.
Sensors 19 00395 g003
Figure 4. Behavior of signals: (a) hormonal signals for “1001”; (b) Ca 2 + signals for “1001”; (c) neural signals for “1001”; (d) hormonal signals for “1011”; (e) Ca 2 + signals for “1011”; (f) neural signals for “1011”; (g) hormonal signals for “1111”; (h) Ca 2 + signals for “1111”; and (i) neural signals for “1111”.
Figure 4. Behavior of signals: (a) hormonal signals for “1001”; (b) Ca 2 + signals for “1001”; (c) neural signals for “1001”; (d) hormonal signals for “1011”; (e) Ca 2 + signals for “1011”; (f) neural signals for “1011”; (g) hormonal signals for “1111”; (h) Ca 2 + signals for “1111”; and (i) neural signals for “1111”.
Sensors 19 00395 g004
Figure 5. Channel capacity vs. neural firing rate, number of released molecules and probability of transmitting “1”. N R R P = 8 , τ f = 30 ms, M A = 3 : (a) λ ¯ = 20; (b) λ ¯ = 25; and (c) λ ¯ = 30.
Figure 5. Channel capacity vs. neural firing rate, number of released molecules and probability of transmitting “1”. N R R P = 8 , τ f = 30 ms, M A = 3 : (a) λ ¯ = 20; (b) λ ¯ = 25; and (c) λ ¯ = 30.
Sensors 19 00395 g005
Figure 6. Channel capacity versus: (a) astrocytes numbers; and (b) vesicle pool models.
Figure 6. Channel capacity versus: (a) astrocytes numbers; and (b) vesicle pool models.
Sensors 19 00395 g006
Figure 7. Transmission delay: (a) delay of hormonal signaling vs. time slot length and distance; (b) delay of Ca 2 + signaling vs. number of supporting astrocytes and time slot length; and (c) delay of neural signaling vs. pool models.
Figure 7. Transmission delay: (a) delay of hormonal signaling vs. time slot length and distance; (b) delay of Ca 2 + signaling vs. number of supporting astrocytes and time slot length; and (c) delay of neural signaling vs. pool models.
Sensors 19 00395 g007

Share and Cite

MDPI and ACS Style

He, P.; Han, X.; Liu, H. Chain Modeling of Molecular Communications for Body Area Network. Sensors 2019, 19, 395. https://doi.org/10.3390/s19020395

AMA Style

He P, Han X, Liu H. Chain Modeling of Molecular Communications for Body Area Network. Sensors. 2019; 19(2):395. https://doi.org/10.3390/s19020395

Chicago/Turabian Style

He, Peng, Xiaojuan Han, and Hanyong Liu. 2019. "Chain Modeling of Molecular Communications for Body Area Network" Sensors 19, no. 2: 395. https://doi.org/10.3390/s19020395

APA Style

He, P., Han, X., & Liu, H. (2019). Chain Modeling of Molecular Communications for Body Area Network. Sensors, 19(2), 395. https://doi.org/10.3390/s19020395

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

Article Metrics

Back to TopTop