Impact Factor 2.536 | CiteScore 4.8
More on impact ›


Front. Comput. Neurosci., 22 January 2019 |

A General Model of Ion Passive Transmembrane Transport Based on Ionic Concentration

  • School of Mathematics, South China University of Technology, Guangzhou, China

Current mainstream neural computing is based on the electricity model proposed by Hodgkin and Huxley in 1952, the core of which is ion passive transmembrane transport controlled by ion channels. However, studies on the evolutionary history of ion channels have shown that some neuronal ion channels predate the neurons. Thus, to deepen our understanding of neuronal activities, ion channel models should be applied to other cells. Expanding the scope of electrophysiological experiments from nerve to muscle, animal to plant, and metazoa to protozoa, has lead the discovery of a number of ion channels. Moreover, the properties of these newly discovered ion channels are too complex to be described by current common models. Hence this paper has presented a convenient method for estimating the distribution of ions under an electric field and established a general ionic concentration-based model of ion passive transmembrane transport that is simple but capable of explaining and simulating the complex phenomena of patch clamp experiments, is applicable to different ion channels in different cells of different species, and conforms to the current general understanding of ion channels. Finally, we designed a series of mathematical experiments, which we have compared with the results of typical electrophysiological experiments conducted on plant cells, oocytes, myocytes, cardiomyocytes, and neurocytes, to verify the model.

1. Introduction

It is well-established that ion channels and pumps generate ion transmembrane transport. The classic models of passive transport mechanisms were derived from the study of neurons, which were considered to be an electric capacitor. The most famous and earliest model, which depended on this idea, was designed by Hodgkin and Huxley (1952). However, in the biological sense, neurons are far from being an electric capacitor. In fact, electrophysiological phenomena are common in all types of cells; currents and potentials have been recorded in myocytes [e.g., denervated frog muscle fibers (Neher and Sakmann, 1976)] and plant cells [e.g., xylem parenchyma of barley roots (Wegner and Raschke, 1994)]. It is inappropriate to think of these cells as capacitors. The so-called process of cell charging, in fact, is the cationic entry into the cell, or the anion leaving; the discharge process is opposite. Thus, a cell is more like an ion container than a capacitor.

If cells are not thought of as capacitors, how should membrane potential and membrane current be read? In fact, membrane potential is a measure of electric charge density, generated by intracellular and extracellular ion distribution. When the internal positive and negative charge are not equal, owing to the electric field force generated by the internal ions which is perpendicular to the membrane, the charged ions inside and outside the cell will approximate to or deviate from the membrane surface, resulting in the non-uniformed distribution of ions.

In addition to the electric fields generated by the internal ions, people can stimulate neurons using an outer artificial electric field. Magnetic stimulation can also stimulate neurons (Walsh and Rushworth, 1999). For example, placing a transcranial magnetic stimulation coil above the skull over a region of interest, when a changing electric current flows through the coil, an electromagnetic field is created (Polson et al., 1982). According to Faraday's law, this induces an electric field that can stimulate cortical neurons (Pashut et al., 2011). However, experimental studies have shown that neurons are insensitive to transverse outer electric field stimulation (Reilly, 1989). Therefore, when we focus on transmembrane ion transport which is perpendicular to the membrane, we can ignore the outer electromagnetic field and consider only the electric field generated by the ions.

For these reasons, although we cannot use the Nernst equation to calculate the membrane potential directly using the ionic concentration before estimating the ion density on the membrane surface, we may still borrow some ideas from the study of the electrical double layer. The electrical double layer is an array of charged particles and oriented dipoles that is thought to exist at every interface (Grahame, 1947). Extending this idea to both sides of the membrane, we can estimate the ion density on the membrane surface. Then the membrane potential measured experimentally, can be calculated numerically as a function of the ion density on the membrane surface. Membrane current is the sum of the charges carried by the transmembrane ion flux. The stimulation current, in fact, generates a redox reaction, and creates an ionic flow to stimulate the membrane.

The classic ion channel model (Hodgkin and Huxley, 1952) was constructed based on the measurement of current-voltage relationships in the membrane of the giant axons of the Loligo (a genus of squid) (Hodgkin et al., 1952) using a voltage clamp for voltage-gated ion channels. Voltage-gated ion channels have been studied thoroughly because it is convenient to measure the electric current while controlling the voltage with a patch clamp. It is generally acknowledged that voltage-gated ion channels are responsible for the electrical activity in a variety of animal cell types; however, voltage-gated ion channels are also found in prokaryotes and diverse range of eukaryotes, and probably exist in all life forms. Furthermore, multicellular animals evolved more than 650 million years ago, whereas nervous systems and muscles evolved much later. For instance, it has been shown that voltage-gated sodium channels predate neurons (Zakon, 2012), like many other neuronal channels and receptors. Therefore, although the voltage-gated ion channels play a major role in the electrical excitability of nerve cells, they have been in existence before nerve system appeared. It is generally believed that voltage-gated ion channels evolved from a prokaryote 2-TM channel. Following the addition of four more domains, an early, ligand-gated, 6-TM protein gave rise to the voltage-gated potassium family and extant ligand-gated potassium channels. Then, following two rounds of gene duplication, this ligand-gated potassium channel formed a four-domain 6-TM channel, such as calcium and sodium channels (Anderson and Greenberg, 2001).

Thus, it is likely that ion channels were created for other purposes. For example, voltage-gated potassium channels often play a role in the depolarization of nerve cells; yet this channel appeared three billion years ago in bacteria and occurs in all organisms (Zakon, 2012). In plants, potassium is mainly absorbed by the roots, which serves an essential role as an osmoticum. Moreover, the potassium channels of root cells, which mediate potassium uptake, play an important role in regulating ion concentration (Hirsch et al., 1998). Some scholars have tried to model potassium uptake by cells via potassium channels (Dreyer, 2017), which is valuable. Therefore, it is not possible to grasp the nature of voltage-gated ion channels by means of electrical signal generation and conduction alone.

Hence, we should realize that the model of ion passive transmembrane transport, which is not universally applicable to different ion channels in the different cells of different species, is insufficient to describe the nature of neural activities. In this paper, we propose a general model.

2. Foundation

2.1. Operators

In this paper, we use some operation symbols for sequences to simplify the expressions: (ai, j)·(bi, j) = (ai, j·bi, j), (ai)·(bi, j) = (bi, j)·(ai) = (ai·bi, j), (ai)·(bi) = (ai·bi); (ai)·(bi)=(aibi), (ai)·b=(aib); (ai,j)(bi,j)=(ai,jbi,j), (ai)(bi)=(aibi); (ai)=iai, (ai,j)=(iai,j).

2.2. Ionic Flux

Electrical signals of nerve are mainly generated by the flow of charged ions. This paper mainly discusses transmembrane ion flow. Thus, spatially, we only consider the direction perpendicular to the membrane surface. The ion flow considered in this paper is mainly caused by diffusion and electromigration,

J=Jd+Je,    (1)

where ȷ = (ȷi), ȷi is the flux density of ion i, abbreviated as ionic flux in this paper. Ions diffuse from high concentration to low, and it is well-known that

Jd=Dcx,    (2)

where D = (Di), Di is the diffusion coefficient of ion i, and c = (ci), ci is the concentration of ion i. The uncommon operational symbols in this paper are explained above. In solution, the velocity of ionic electromigration is generally considered to be proportional to the electric field intensity ϵ, so the ionic flux of electromigration

je=ϵ · czμ,    (3)

where μ = μi, μi=DiFRT is experiential, R is the molar gas constant, F = eNA is the Faraday constant, T is the temperature (with a default of 290 K in this paper), and z = (zi), zi is the valence of the ion i. It is difficult to measure ionic flux directly, but easy to measure the current generated by the ions flowing and the redox reaction. The current density (hereafter referred to as current for short) generated by ionic flux is

ı=F·z·j.    (4)

Therefore, it is easy to calculate the current from ionic flux, but difficult to calculate ionic flux from the current. When the cell membrane is referred to as being given a stimulation current ıs, the current is not directly stimulating the membrane, but is involved in the redox reaction on the surface of the electrode, which generates an ionic flux that we refer to as stimulation ionic flux in this paper. There are several common electrodes, for the sake of convenience, here we choose silver electrodes (Ag/AgCl), on whose surface, the redox reaction is

AgCl+eAg+Cl.    (5)

The model would be too complicated if we considered the pipette solution. So, ignoring the interspace of the pipette, from Equations 4 and 5, the stimulation ionic flux can be estimated by

jsi={1F·ls(i=Cl)0(iCl)    (6)

where ıs is the intensity of stimulus current.

2.3. Ionic Concentration

How ionic flux is estimated by ionic concentration and electric field has been given above. Since, in this case, the electric field is generated by the ions, the ion flux can be estimated by the ionic concentration distribution. A membrane, in any shape, separates a solution into two regions: internal and external. Thus, the morphology of the membrane affects the distribution of ions intracellular and extracellular. There are various shapes of cells, such as sphere, polyhedra, spindle and cylinder. For convenience, the shapes are generally simplified, usually to cylindrical and spherical shapes.

2.3.1. Cylinder

It is tractable for mathematical analysis to reduce the complex branching structure of a dendritic tree to a simple cylinder (Holmes, 2015). In 1962, Rall showed that dendritic morphology could be reduced to a simple cylinder (Rall, 1962). This method was then applied to the study of axons. Roll the membrane into a cylinder, and let the inner surface of the membrane face inward and the outer surface face outward. Assume that the cylinder is standard, of which the length is long enough, and the ion distribution depends only on the distance to the axis of the cylinder. Then the rate of change of ionic concentration with time depends on ion flux,

tc=(1x+x)J,    (7)

where x is the distance to the axis of the cylinder, at x = r is the membrane, if for any x the ionic distribution is uniform, ignoring the thickness of the membrane. If ignore the outer electromagnetic field, then the electric field is generated by the ions and is perpendicular to the membrane surface,

ϵ(y)=zyε(y)·F0yc·x dx,    (8)

where ε is the permittivity.

2.3.2. Sphere

Since isolated cells in cultures or suspensions often have a spherical shape (Drasdo and Hohme, 2005), modeling a single cell as a sphere is also common, especially for electrophysiological experiments (Debruin and Krassowska, 1999). Packaged the membrane into a sphere, and also let the inner surface of the membrane face inward and the outer surface face outward. Assume that the sphere is standard, and the ion distribution depends only on the distance to the center of the sphere. Then the rate of change of ion concentration with time depends on ion flux,

tc=(2x+x)J,    (9)

where x is the distance to the center of the sphere, at x = r is the membrane, if for any x the distribution of ions is uniform, ignoring the thickness of the membrane. If we ignore the outer electromagnetic field, then the electric field is generated by the ions and is perpendicular to the membrane surface,

ϵ(y)=zy2ε(y)·F0yc·x2 dx.    (10)

2.4. Ionic Density on Surface

To simulate whole cell patch clamp method, focusing on transmembrane transport, the sphere is more appropriate.

2.4.1. Internal

The concentration of the solution in the experiment is, in general, the average concentration. The average internal ion concentration can be written as

Ainc=3r30rc·x2 dx,    (11)

and Ainci is traditionally written as [i]in. Transmembrane ion flow makes a small but significant change in the average internal ionic concentration that

ddtAinc=3r1·u,    (12)

where u = −ȷx = r is the transmembrane ion influx (then u < 0 means outflow, as a salute to the classical model). In practical application, the change of average internal ion concentration can be directly calculated by transmembrane ion flux, instead of integrating. Set

Qc=r20rc·x2 dx    (13)

to describe the total amount of intracellular ions. Then the internal charge (per unit area) Q can be calculated by internal ionic concentration,

Q=F·zQc,    (14)

of which the electric field derived from Equation 10 is

ϵ(r)=Fε1·zQc.    (15)

Additionally, in practical application, the change in internal charge Q can be directly calculated by transmembrane ion flux u,

ddtQ=F·uz,    (16)

instead of integrating. Thus, the ion density on the inner surface of the membrane, referred to as internal ionic density for short, can be calculated by

Ainc=r2rλr​​​​cx2dx,    (17)

where λ = (λi), λi is the diameter of the ion i (see Appendix).

2.4.2. External

For [0, r) is the internal space, (r, ∞) is external. The average external ion concentration is

Aexc=limy3y3r3ryc·x2 dx,    (18)

and Aexci is traditionally written as [i]ex, which is constant when the external space is large enough. At the same time, the ion density on the outer surface of the membrane, named external ionic density for short, should be

Texc=1r2rr+λcx2dx.    (19)

2.4.3. Estimating

Borrowing from the theory of electrical double layer (Grahame, 1947), the internal (external) solution could be thought to consist of two layers of ions: one of adsorbed ions and the other consisting of an ionic atmosphere. Thus, the membrane in contact with both the internal and external solution actually forms four layers. If the net charge is tiny, |z Qc|≪ΣQc, it can be derived from the equations in Ionic Concentration so that the net charge is very close to the membrane and is regarded as being adsorbed on the membrane surface, while the diffuse layer is relatively stable. Therefore, ȷ|x-r|>λiȷx=r0 or ȷ≡0, as

cicixϵFziRT,    (20)

for |xr| < λi. Thus, in the actual operation, a simple and convenient arithmetic can be used to estimate the membrane ionic density:

TincλAinc(1χ z Q c ·λz)(1+χz Q c · λz),    (21)
TexcλAexc(1χ z Q c ·λz)(1+χz Q c · λz),    (22)


χ=F22RTεm.    (23)

After measuring the membrane permittivity εm = ε(r), the membrane ionic density can be estimated.

2.5. Membrane Potential

In contact with the internal and external solution, the redox reaction on the electrode produces a potential difference, that is the membrane potential which is measured. The Nernst Equation is usually used to estimate the voltage, when the electrostatic force is absent, as ϵ = 0. However, it is also common that the internal ions almost inevitably generate an electric field force, which causes the ions to attach or leave the membrane surface. Actually, the membrane potential depends on the membrane ionic density,

Vm=V(TinTex)c.    (24)


V(ab)  =RTFz· ln (ab).    (25)

Therefore, after estimating the ionic density on the membrane, we could calculate the membrane potential difference under the electric field.

It can be seen that dependence on membrane potential is essentially a special case of dependence on membrane ion density.

In the whole cell patch clamp experiment, the pipette is connected to the intracellular fluid. Then, the difference between the stimulus voltage and the membrane potential produces a stimulus current, which is difficult to accurately calculate; however, using an approximation method, and ignoring the spatial difference on the membrane, it can be estimated as

ıs=VsVmRs·S,    (26)

where Rs is the resistance of the pipette filled with solution, Vs is the electrode voltage, and S is the superficial area of the cell. If the location of the patch clamp is taken into account, the point on the membrane close to the hole pierced by the patch clamp is much more stimulated. With Equation 6, the ionic flux of stimulation can be estimated.

3. Ion Channel Model

Using the analysis above, one can see that the membrane potential is actually a measured value of the membrane ion density. Therefore, voltage-sensitive ion channels are essentially sensitive to the ionic density on the membrane surface. To maintain the ion concentration, ion channels are selective and sensitive to ions. Especially considering the evolutionary history of ion channels, we speculate that ion channels have a general functional mechanism. Thus, it is necessary to establish a general ion channel model to describe the passive transport of ions across cell membranes. Here we set three basic assumptions based on the discussion above:

Assumption 1. The primary role of ion channels is to maintain specific intracellular ionic concentration.

Assumption 2. Ion channels can only detect ions on the inner surface of the membrane and the ions flowing through it.

Assumption 3. Ion channels can filter ions, but the flow of ions is passive.

The transmembrane passive ionic flux is related to the ionic density difference and the selective permeability of ion channels. In a membrane, considering channels of the type h as a whole (called channel h), the ionic flux depends on its selective permeability and the ionic density difference,

ζ=ρDc,    (27)

where ρ = (ρh, i), ρh, i is the selective permeability of channel h for ion i, ζ = (ζh, i), ζh, i is the influx of ion i through channel h and

D=TinTex.    (28)

Thus, the transmembrane ionic influx is

u=ζ+Js    (29)

if ion pump is absent, where ȷs is the stimulation ion influx. As generally understood, a channel can be abstracted as a combination of filter and gate. The filter that selects the ions is either conductive or inactivated, and the gate is either closed or open. Therefore, it can be designed that

ρ=fgϱ,    (30)

where f = (fh, i), fh, i describes the filter of channel h to ion i, gh describes the gate of channel h, and ϱ = (ϱh, i), ϱh, i is the standard selective permeability of channel h to ion i.

3.1. Sensor

Generally, the intracellular ion concentration is stable. It is assumed (Assumption 2) that the ion channel can only detect the concentration of ions on the inner surface of the membrane. It is reasonable to assume that each type of ion channel has a specific internal ion density at which it is most comfortable. Assumptions 1 and 2 require that the ion channels are sensitive to the specific ions. For each channel h, the internal ionic density detected, (named detected ionic density in this paper written as σh), depends on the internal ionic density and the sensitivity sh, i to ion i. Assuming that the dependency is linear, then

σ=Tinc·s,    (31)

and mark σ¯ for optimum, named optimum density. In the same way, with the sensitivity s, channel h detects the ionic flux by

ψh=ζh·sh,    (32)

named detected ionic flux in the same manner. We also know that the sensitivity s does not have to be proportional to the permeability ϱ. To be specific, channel h may be sensitive to ion i but allows ions j(≠i) to pass through. To facilitate the modeling below, we have designed two virtual indicators to describe the tendency of detected ionic density to deviate from the optimal density for ion channels — ϑ to point out the difference between detected ionic density and optimal density, and ϕ to point out the influence of detected ionic flux on density. According to usual practice, ϑ points to 0 for the optimum, as σ/σ¯=1; −π/2 and π/2 for the worst, as σ/σ¯=0,; ϕ points to 0 for the optimum, as ψ/σ = 0. To satisfy these conditions, it is natural to state that

ϑ=2atanσσ¯π2,    (33)


ϕ=τψσ¯,    (34)

where τ = (τh), τh is the empirical time constant of channel h. The larger τh is, the more stringent the restriction of gate on ion flow is.

3.2. Gate

The gate of the ion channel is a switch that controls the passive transportation of the intracellular and extracellular ions. According to the basic assumptions, an ion channel whose primary role is to maintain the intracellular ionic concentration, can only detect ionic density on the inner surface of the membrane and the flow of ions through it. When the detected ionic density is optimal, the ion channel h needs to keep this state and the gate should be closed. When the detected ion density deviates from the optimal density, σhσ¯hϑh>0, the ion channel needs to transfer ions. Given that this transportation is passive, the result may be better or worse. Thus, the opening and closing of the gate should be controlled by two parts,

g=vw.    (35)

One is to explore and the status is measured by vh ∈ [0, 1]. As assumed before, vh → 0 for ϑh = 0, vh → 1 for ϑh=±π2. Then we can obtain geometrically that

ddtv=α(v¯v),    (36)

where the limiting value

v¯=1cos·ηϑ,    (37)

α = (αh), αh describes the agility of exploring of channel h; and η = (ηh), ηh ∈ (0, ∞) is the tolerance of channel h to ionic density deviation, describes the degree of which the channel allows the ion density to deviate from the optimal value. The value of ηh is empirical and uniformly set as 1 in this paper for each h. The other is to judge and is measured by wh ∈ [0, 1]. While vh > 0:

• if σh approximates to the optimum, as ϕhh < 0, the judgement is to keep the gate open;

• if σh deviates from the optimum, as ϕh/ ϑh > 0, the judgement is to close the gate;

• otherwise, as ϕhh = 0, for random is the vibration of ions, wh → 1/2 from the viewpoint of probability.

The execution of judgment still needs a process, and the bigger vh is, the more urgent the process is. Thus, we state that

ddtw=β(w¯w)·(v+δ),    (38)

of which the objective is

w¯=1πacotϕϑ    (39)

to satisfy the above conditions about w, where β = (βh), βh describes the agility of judgment of channel h, δ = (δh), δh is positive, empirical, and constant on h, then βhδh(w¯h-wh) is the rate of change of channel h while vh = 0. If the detected ionic densities of all the channels on a membrane are optimal or very close to the optimum, then the membrane is in what is called a resting state.

3.3. Filter

It is generally acknowledged that ion channels filter the ions, like having a filter. When an ion channel is continuously stimulated, the filter changes the permeability of the ions. For example, application of a specific stimulus opens the intracellular gate of a potassium channel (activation), yielding a transient period of ion passive conduction until the selectivity filter spontaneously undergoes a conformational change toward a non-conductive state (inactivation). Removal of the stimulus closes the gate and allows the filter to interconvert back to its conductive conformation (recovery) (Ostmeyer et al., 2013). In addition, there is a remarkable correlation between the degree of gate opening and the conformation and ion occupancy of the selectivity filter (Cuello et al., 2010). For instance, the hERG channel opens when it is stimulated specifically and then gradually inactivates, but still allows the sodium to pass through in its specific inactive state (Gang and Zhang, 2006). In conclusion, we may consider that the filter occupied by the ion flow may gradually inactivate under continuous stimulation. It is inferred that f¯h,i decreases with the stimulus intensity (measured by ϑh2) and detected ion flux (measured by ϕh2). Following the description of the filter state above, set f¯h,i=1 at ϕh = 0 and ϑh = 0. The model of filter can then be designed as

ddtf=γ(f¯f),    (40)
f¯=exp(κϕ2mϑ2n),    (41)

γ = (γh), κ = (κh, i), where γh and κh,i+, mh, i and nh, i ∈ ℕ are experiential, γh describes the agility of the filter, and κh, i is the inactivation coefficient of channel h occupied by ion i.

4. Results

In Foundation, we have introduced how ionic flux affects ion concentration, and how to calculate membrane potential from ionic concentration and membrane current from passive transmembrane ionic flux. We then constructed a ion channel model, which can be used to simulate how the ion distribution affects the passive transmembrane ion flux. Thus, we may compare the numerical simulation results with the measured experimental data to validate the model and discuss its nature.

Although cells are different in shape and size, given our focus on the properties of membranes and channels, this does not affect the analysis of the model to unify the cells into a sphere with radius r = 10 μm, especially for ion channels.

4.1. Estimation of Capacitance

Here we use the method outlined in Foundation to calculate the membrane capacitance and compare it with the actual measurement. The detailed method is to calculate the membrane potential and internal charge first and then to estimate the membrane capacitance by Cm = dQ/dVm. To simulate a real environment and close all the channels, given common (AexAin)c and ıs, the relationship between intracellular electric charge Q and membrane potential Vm will be obtained, with which the membrane capacitance Cm with the potential can then be calculated. In reality, the relative permittivity of cytomembrane εr is generally between 1.4 and 20 (Zhang et al., 2014). Let εr = 4, for Vm/mV ∈ (−100, 100), the estimated membrane capacitance Cm1 μF/cm2, is consistent with the credible measured values of giant axons (Hodgkin et al., 1952); see Figure 1.


Figure 1. Numerical simulation of membrane capacitance. Relative permittivity (as default in this paper): εr = 4. Extracellular solution: (i). Intracellular solution: (ii). It shows that the membrane capacitance is not symmetrized to Vm = 0 mV — a phenomenon is common in all kinds of media (Hamelin et al., 1983).

4.2. Gate Properties

This subsection mainly analyses the nature of the inner gate in the model, which can be verified by the comparison of digital simulations and biological experiments.

4.2.1. Potassium Channel From Plant Root Cell

First simulate a plant root cell potassium channel by our model, of which the phenomenon is seen in plant root potassium channel electrophysiological experiments: if the external solution is rich in potassium ions, the hyperpolarization of wild-type root cell produces an inward current, which disappears when extracellular K+ is replaced by Cs+ (Hirsch et al., 1998). The numerical simulation is shown in Figure 2. Based on the properties of the model, in resting state:

• Holding Vs > V, the intracellular cations move to the inner surface of the membrane and the extracellular cations leave the outer surface. Then the inner surface potassium density is higher than the density in resting state, σ>σ¯v¯>0, and the potassium channel will be open. As an experimental environment in general, the inner surface potassium density is greater than the outer surface, ζ < 0, so w¯>1/2, and the gate will stay open.

• Holding Vs < V, the intracellular cations leave the inner surface of membrane and the extracellular cations move to the outer surface. Then the inner surface potassium density is lower than the density in resting state, σ<σ¯v¯>0, and the potassium channel will also be open.

– If the extracellular potassium ion is scarce, [K+]ex[K+]in0, the intracellular potassium will be transported passively to the outside of the cell, then w¯0, and the gate of the potassium channel will be closed to stop the potassium outflow.

– Otherwise, if the external potassium concentration is high enough, ϑ < 0 and ϕ>0 ⇒ w¯>1/2, and the potassium channel will stay open.


Figure 2. Numerical simulation current I(A/m2) of potassium channels in wild-type root cells. Channel list: Kv♣. Patch clamp: Rs = 2.4 MΩ; Vs = −40 mV for t/ms < 0, Vs = −200 to 120 mV for t/ms>0. Intracellular solution: (iv). Extracellular solution: (v) for (A), (vi) for (B). The points in (C) are the values of current recorded at t = 1s, black for (A), white for (B).

4.2.2. Sodium Channel and Refractory Period

In common patch clamp experiments, the intracellular solution is rich in potassium ions, and the sodium ion concentration is extremely low. Therefore, the sodium channel is a good contrast with the potassium channel for discussing the properties of the gate. In sodium channel experiments, the external solution is generally rich in sodium ions, [Na+]in[Na+]ex. In such an environment, when stimulating the membrane, the sodium channel opens briefly and then enters the refractory period. This phenomenon can be seen in the early experiment on the giant axon of the Loligo (Hodgkin and Huxley, 1952), of which the numerical simulation result is shown in Figure 3. After the stimulus is removed, the sodium channel will gradually recover. This phenomenon has been well documented in the study of excitable membranes (Hille, 2001), of which the numerical simulation result by our model is shown in Figure 4. By the model's calculation, in resting state:

1. Given Vs > V, the sodium ion density on the inner surface is higher than the resting state, σ>σ¯vv¯>0, and the sodium gate should be open as the channel is activated. However, when the external sodium ion concentration is much greater than the internal, as the gate is opening, the sodium ions will flow backward, ww¯0, and the gate will close again. In the case that α is large enough, such as α/β>0.1, the gate opens extremely fast. In the time between the opening and closing of the gate, there is a short but significant sodium ion influx. As ww¯0, the sodium gate tends to be closed, the channel inactivates, and the refractory period commences.

2. Once the stimulus is eliminated, w¯1, the channel will gradually recover. Subsequent w≈1 means the end of refractory period.


Figure 3. Numerical simulation from activation to inactivation of sodium channel. (A) shows the membrane current, of which the unit is A/m2, generated by sodium, potassium, and chloride ion flux; (B) the sodium influx, the unit is μmol/s/m2. Channel list: Nav♠, Kv♠, Clv♠. Patch clamp: Rs = 2.4 MΩ; Vs = −70 mV for t/ms ≤ 0, Vs = −70 mV to 30 mV for t/ms>0. Extracellular solution: (i). Intracellular solution: (ii).


Figure 4. Numerical simulation of recovery from refractory period. Patch clamp: Vs = −70 mV for t/ms ∈ [0, Δt], Vs = 30 mV for t/ms∉[0, Δt], Δt = 0–25. The other settings are the same as in Figure 3. With the Δt increasing, the sodium channel gradually revives to the resting state as the end of the refractory period.

4.3. Filter Properties

It was easy to infer from the concentration of the external solution that the outward current initially found on nerve cells was generated by the potassium ion flow. Based on the axon experiments, it was believed that the potassium channels would not be closed before the depolarization stimulus was removed. However, in the subsequent experiments, transient outward potassium current was observed.

4.3.1. Transient Outward Potassium Channel

The typical phenomena of transient outward current can be seen in experiments conducted on neural somata (Connor and Stevens, 1971) and cardiomyocytes (Akar et al., 2004), of which the numerical simulation results by our model are shown in Figures 5, 6.


Figure 5. Numerical simulation of transient outward potassium current. (A) Activation and inactivation curves. (B) Steady-state inactivation curves. Channel list: Kv4.3. Patch clamp: Rs = 2.4 MΩ; Vs = −70 mV for t/ms ≤ 0, Vs = −100–40 mV for t/ms>0. Extracellular solution: (i). Intracellular solution: (ii).


Figure 6. Recovery from inactivation of transient outward current by numerical simulation. Vs = 50mV for t/ms ∈ (−∞, 0)∪(Δt, ∞), Vs = −70mV for t/ms ∈ [0, Δt], Δt = 0 to 200. The other settings are the same as in Figure 5.

As we have mentioned above, the conformation of the selectivity filter is correlated to the degree of gate opening, and more specifically, depends on the detected stimulation intensity φ and the occupation ϕ by the ions. κ determines the type of ion channel inactivation. Based on the model calculation, in the environment of [K+]ex/[K+]in0, given β≫α≫κ:

• In the case of Vs < V, σ<σ¯. For the reason that SincK>SexcK ([K+]ex/[K+]in0), as the gate is opening, the potassium ion outflow leads to ww¯0, and the gate will be closed. For β ≫ α ≫ κ, the gate can be considered to always be closed so that the filter state remains constant.

• In the case of Vs>V, the gate will be open and will expand for the duration of the stimulation. For α ≫ κ, the change of f lags behind g. The potassium current will first increase significantly and then decrease. The peak increases with the stimulation. See Figure 5. After the stimulus is removed, the gate will be closed again and the inactive filter will gradually recover to the resting state. The longer the interval between two stimuli is, the closer the second peak is to the first; see Figure 6.

4.3.2. Delayed Rectifier Potassium Channel

As we have said earlier, κ determines the type of ion channel inactivation. In contrast to the transient outward potassium channel, is the rectifier potassium channel. The typical phenomena can be observed in an experiment conducted on hERG (Ficker et al., 1998), of which the numerical simulation by our model are shown in Figures 5, 6. By the model calculation, given κ ≫ α:

1. Holding Vs>V, the channel activates slowly on moderate membrane depolarization and exhibits limited outward current at a more positive membrane depolarization. After the membrane current is relatively stable, on return to negative membrane potentials, removal of C-type inactivation results in an initial strengthening in tail currents with slow deactivation kinetics, for dfdt-dgdt. The stronger the previous depolarization, the stronger the tail current. See Figure 7.

2. At a certain point in time after the outward current is relatively stable, the switch of Vs from depolarization to hyperpolarization of Vs<V will generate an instant inward current. The longer the depolarization lasts, the stronger the instant inward current is. See Figure 8.


Figure 7. Numerical simulation current of hERG. Channel list: hERG. Patch clamp: Rs = 2.4MΩ, Vs = −70mV for t/ms ≤ 0, Vs = −50mV to 70mV for t/ms ∈ (0, 300], Vs = −50mV for t/ms>300. Extracellular solution: (i). Intracellular solution: (ii). The black points are recorded at t = 0.3s, the white points are the peaks of the tail currents in (0.3, 0.5) of t/s.


Figure 8. Activation current of hERG by numerical simulation. Vs = −70mV for t/ms ≤ 0, Vs = 50mV for t/ms ∈ (0, Δt), Vs = −150mV for t/ms≥Δt, Δt = 0 to 400. The other settings are the same as in Figure 7.

Transient outward and delayed rectifier potassium channels are common in cardiomyocytes, of which the numerical simulation shows a general applicability to our model. We have analyzed the relationship between the filter properties and the types of ion channel inactivation within some typical cases.

4.4. Special Sensitivity

There are some cases we have given above in which sensitivity and permeability correspond incompletely (ϱh||ϱh||sh||sh||). To be more extreme, it is an absolute necessity to consider whether there are some channels particularly sensitive to certain ions but very permeable to others. An affirmative answer to this question may indicate a dual dependence on specific ion concentration and voltage of these channels. In reality, calcium-activated potassium channels are common and sensitive to calcium ions but permeable to potassium. Since all the channels mentioned above are sensitive and transparent to cations, we are naturally led to ask whether there are some ion channels sensitive to cations but permeable to anions. The calcium-activated chloride channel (CACC) is a good example.

4.4.1. Calcium-Activated Chloride Channel

Let us assume that a certain type of channel can be very sensitive to calcium ions but permeable to chlorine, and let [Ca2+]ex[Ca2+]in0.

• If [Ca2+]in is low enough to satisfy σσ¯ in the initial state, then w=w¯0g=g¯0, and the gate will remain almost closed under any stimulation.

• If [Ca2+]in is moderate, σσ¯g¯0 for Vs = V, the channel is almost closed for Vs = V, but open for VsV because of σσ¯ under this circumstance.

• If [Ca2+]in is high enough to satisfy σσ¯ for Vs = V, easily derived g¯>0 in this circumstance, the channel is already open before the stimulus.

Therefore, this type of chloride channel has a dual dependence on membrane potential and intracellular calcium concentration, of which the peculiarities, shown in Figures 9,10, are very consistent with the CACC (the samples were taken from xenopus oocytes) that play fundamental roles in physiological processes in many tissues (Kuruma and Hartzell, 2000).


Figure 9. Numerical simulation activation of CACC currents Channel list: CACC. Patch clamp: Rs = 4.0MΩ, Vs = 0mV for t/s < 0.1, −120mV to 120mV for t/s ∈ (0.1, 1], −120mV for t/s>1. Intracellular and Extracellular solution: (vii), using Ca-EGTA to control [Ca2+]in. The unit of I is A/m2, and the unit of t is ms.


Figure 10. Numerical simulation steady state current-voltage relationships of CACC. The values of electric current density recorded at t = 1s in Figure 9, Vs from −120mV to 120mV, white for 70 nmol/L, yellow for 600 nmol/L, blue for 1μmol/L, black for 2μmol/L.

4.5. Combined Effect of Channels

It is well-documented that there is more than one type of ion channel in a cell membrane. In our model, the ion channels are responsible only for the density of specific ions on the inner surface of the membrane. The state change of the ion channel depends on the ionic density on the inner surface of the membrane, such that there is no direct correlation among the ion channels. However, because the opening of one ion channel affects the distribution of ion concentrations in the membrane and thus affects the states of other ion channels, the ion channels on the same membrane interact with each other via ions, and produce a variety of complex electrophysiological phenomena, such as the nerve impulse, which is a result of the combined effect of the ion channels. In general, in experimental environments, extracellular fluid is rich in sodium, [Na+]in<[Na+]ex, intracellular fluid is rich in potassium, [K+]in>[K+]ex, and both are rich in chlorine, which is also the main electrode reactant. Thus, from the numerical model standpoint, the combination of sodium, potassium, and chlorine channels should be relatively sufficient; this is exactly the general structure observed in neurons.

4.5.1. Nerve Impulse

There exist several simple methods to observe changes in action potentials. From the perspective of our model, let us analyse the process of neural excitation.

Holding Vs < V, based on the analysis above, it can be concluded that the sodium and chlorine channels are open, whereas the potassium channel is almost closed. As the voltage clamp is released, the closure of the sodium and chlorine channels will take a short time. It takes a transient stimulus to produce a potential spike (see Figure 11), which is generated conjointly by three different channels (see Figure 12). The corresponding classic experimental phenomena were recorded in an earlier paper (Hodgkin and Huxley, 1952).


Figure 11. Numerical simulation nerve action potential of anode break excitation following sudden cessation of external current. Channel list: Nav♠, Kv♠, Clv♠. Stimulus current: Is/A/m2=-0.6 to −0.9 for t/ms ∈ [−50, 0], 0 for t/ms>0. The other settings were the same as in Figure 3.


Figure 12. Numerical simulation transmembrane ionic flux of anode break excitation following sudden cessation of external current in the nerve. The transmembrane ionic fluxes (μmol/s/m2) of sodium and potassium with time (ms) are shown. The settings were the same as in Figure 11.

Given a depolarized stable Is, as mentioned, the sodium ion channel will be open briefly. In this process, the membrane potential will increase sharply and the potassium and chloride ion channels will also be opened to make the membrane potential fall back down again.

• If the depolarizing direct current stimulus remains at an appropriate intensity, the channels may be excited repeatedly, forming periodic spiking (Ram-rez Piscina and Sancho, 2018), of which the numerical stimulation is shown in Figures 13, 14.

• If the alternating current stimulation is given with appropriate frequency and intensity, then bursting is generated (Figure 15).


Figure 13. Numerical nerve action potential of periodic spiking. The constant direct current stimulation: Is=0.1A/m2. The other settings were the same as in Figure 3.


Figure 14. Numerical simulation transmembrane ionic flux density of periodic spiking. The numerical ion fluxes (μmol/s/m2) of sodium, potassium and chloride with time (ms) are shown in order. The settings were the same as in Figure 13.


Figure 15. Numerical simulation nerve action potential of bursting. The alternating current stimulation: Is=0.2A/m2·sin20πt/s. The other settings were the same as in Figure 13.

The parameters of the sodium, potassium, and chloride ion channels in the model determine the frequency of periodic spiking. By comparing Figure 13 with 14, it is clear that the oscillations of ionic fluxes are much more complex than the membrane potential.

5. Discussion

In this paper we attempt to establish a general model of ion passive transmembrane transport based on ion concentration that conforms to the general understanding of ion channel today. The model is simple but capable of explaining and simulating the complex phenomena of patch clamp experiments, and is universally applicable to different ion channels in the different cells of different species.

In experiments, it has been found that many kinds of stimulation can produce electrophysiological phenomena, and the action of ion channels also seems to depend on a variety of external factors. These phenomena make the properties of ion channels so difficult to understand that some people believe that different models must be built to explain the different phenomena.

For example, a model based on the cooperative activation of sodium channels that reproduces the observed dynamics of action potential initiation to simulate the features of the initiation dynamics of cortical neuron action potentials has been proposed (Bjorn et al., 2006). However, it was shown that these features are to be expected from the model of Hodgkin and Huxley, if the spatial geometry and spike initiation properties of cortical neurons are taken into account. It is therefore unnecessary to invoke exotic channel-gating properties as an explanation (McCormick et al., 2007).

We believe that building different models for different phenomena, without explaining the commonality by a general model, will not help us further our understanding. It has been demonstrated countless times that the closer a model gets to the truth, the simpler it becomes. When we combine the factors which affect the ion channels, it is natural to establish a model based on membrane surface ion density. In this way, all the ion channels in this paper have been described by this general model, of which the numerical simulation results are consistent with the electrophysiological phenomena, universal to plant cells, oocytes, myocytes, cardiomyocytes, and neurocytes. We chose these experiments not only based on cell types, but also because their phenomena are representative.

These results suggest that when ion channels are responsible for intracellular ionic concentration, they are sufficient to produce a variety of electrophysiological phenomena.

6. Vista

6.1. Pump

Without the action of the pump, the intracellular ion concentration that deviated from the resting state due to passive transmembrane transport would not completely return to initial levels. A common model for the pumps has been built in our laboratory, though not yet validated by precise experimental data. The combination of ion channels and pumps constitute the complete transmembrane ion transport, which we will implement in the future. Once this is completed, our model will explain the process of intracellular ion concentration formation and maintenance.

6.2. Conduction

In this paper, the ionic flux is perpendicular to the membrane surface. To add the ionic flux parallel to the membrane surface, the complete picture of ionic flux will be calculated, especially for dendrites and axons. The early models only focus on the internal and transmembrane current or ionic flux of cells, generally. However, it is not difficult to deduced that an external ionic flux also exists. Further, in dendrites and axons, the direction of ion flow can be divided into transmembrane and axial and the axial ion flow is sensitive to the axial electric field. So the effect of the outer electromagnetic field is a very important subject, which needs to be further studied. Also, our model may extend to the connection between cells; for example, connections between nerve and nerve (synapse) and nerve and muscle (neuromuscular junction).

6.3. Information

Ion channels are sensitive to specific ions and control the passive transmembrane transport of ions. Therefore, the ionic flux and concentration can be regarded as the signals of cell state. From ionic density to voltage and ionic flux to electric current, there is a loss of information, which is hindering the exploration of the human brain. The artificial neural network theory based on the classical model is also inherently flawed. Our future research will try to distinguish and understand the ion signals of cells. On that basis, the establishment of a new artificial neural network model may contribute to the development of artificial intelligence.

Author Contributions

VW: Conceptualization, software and writing; SL: Supervision.


This work was supported by the National Natural Science Foundation of China, 11572127 and 11872183 (to SL).

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.


Akar, F. G., Wu, R. C., Deschenes, I., Armoundas, A. A., Rd, P. V., Houser, S. R., et al. (2004). Phenotypic differences in transient outward k+ current of human and canine ventricular myocytes: insights into molecular composition of ventricular ito. Am. J. Physiol. Heart Circ. Physiol. 286:H602. doi: 10.1152/ajpheart.00673.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, P. A. V., and Greenberg, R. M. (2001). Phylogeny of ion channels: clues to structure and function. Comparat. Biochem. Physiol. B Biochem. Mol. Biol. 129, 17–28. doi: 10.1016/s1096-4959(01)00376-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Bjorn, N., Fred, W., and Maxim, V. (2006). Unique features of action potential initiation in cortical neurons. Nature 440, 1060–1063. doi: 10.1038/nature04610

CrossRef Full Text | Google Scholar

Connor, J. A., and Stevens, C. F. (1971). Voltage clamp studies of a transient outward membrane current in gastropod neural somata. J. Physiol. 213, 21–30.

PubMed Abstract | Google Scholar

Cuello, L. G., Jogini, V., Cortes, D. M., and Perozo, E. (2010). Structural mechanism of c-type inactivation in k+ channels. Nature 466:203. doi: 10.1038/nature09153

PubMed Abstract | CrossRef Full Text | Google Scholar

Debruin, K. A., and Krassowska, W. (1999). Modeling electroporation in a single cell. i. effects of field strength and rest potential. Biophys. J. 77:1213.

PubMed Abstract | Google Scholar

Drasdo, D., and Hohme, S. (2005). A single-cell-based model of tumor growth in vitro: monolayers and spheroids. Phys. Biol. 2, 133–147. doi: 10.1088/1478-3975/2/3/001

PubMed Abstract | CrossRef Full Text | Google Scholar

Dreyer, I. (2017). Plant potassium channels are in general dual affinity uptake systems. Aims Biophys. 4, 90–106. doi: 10.3934/biophy.2017.1.90

CrossRef Full Text | Google Scholar

Ficker, E., Jarolimek, W., Kiehn, J., Baumann, A., and Brown, A. M. (1998). Molecular determinants of dofetilide block of herg K+ channels. Circ. Res. 82:386.

PubMed Abstract | Google Scholar

Gang, H., and Zhang, S. (2006). Na+ permeation and block of herg potassium channels. J. Gen. Physiol. 128, 55–71. doi: 10.1085/jgp.200609500

PubMed Abstract | CrossRef Full Text | Google Scholar

Grahame, D. C. (1947). The electrical double layer and the theory of electrocapillarity. Chem. Rev. 41, 441–501.

PubMed Abstract | Google Scholar

Hamelin, A., Vitanov, T., Sevastyanov, E., and Popov, A. (1983). The electrochemical double layer on sp metal single crystals: the current status of data. J. Electr. Chem. 145, 225–264.

Google Scholar

Hille, B. (2001). Ion Channels of Excitable Membranes. Sunderland: Sinauer Associates, Inc.

Google Scholar

Hirsch, R. E., Lewis, B. D., Spalding, E. P., and Sussman, M. R. (1998). A role for the akt1 potassium channel in plant nutrition. Science 280, 918–921.

PubMed Abstract | Google Scholar

Hodgkin, A. L., Huxley, A., and Katz, B. (1952). Measurement of current-voltage relations in the membrane of the giant axon of loligo. J. Physiol. 116:424.

PubMed Abstract | Google Scholar

Hodgkin, A. L., and Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117, 500–544.

PubMed Abstract | Google Scholar

Hodgkin, A. L., and Huxley, A. F. (1952). Currents carried by sodium and potassium ions through the membrane of the giant axon of loligo. J. Physiol. 116, 449–472.

PubMed Abstract | Google Scholar

Holmes, W. R. (2015). Equivalent Cylinder Model (Rall). New York, NY: Springer.

Google Scholar

Kuruma, A., and Hartzell, H. C. (2000). Bimodal control of a Ca2+-activated Cl channel by different Ca2+ signals. J. Gen. Physiol. 115:59. doi: 10.1085/jgp.115.1.59

CrossRef Full Text | Google Scholar

McCormick, D. A., Shu, Y., and Yu, Y. (2007). Hodgkin and huxley model—still standing? Nature 445, E1–E2. doi: 10.1038/nature05523

PubMed Abstract | CrossRef Full Text | Google Scholar

Neher, E., and Sakmann, B. (1976). Single-channel currents recorded from membrane of denervated frog muscle fibres. Nature 260, 799–802.

PubMed Abstract | Google Scholar

Ostmeyer, J., Chakrapani, S., Pan, A. C., Perozo, E., and Roux, B. (2013). Recovery from slow inactivation in K+ channels is controlled by water molecules. Biophys. J. 501, 121–124. doi: 10.1038/nature12395

PubMed Abstract | CrossRef Full Text | Google Scholar

Pashut, T., Wolfus, S., Friedman, A., Lavidor, M., Bargad, I., Yeshurun, Y., et al. (2011). Mechanisms of magnetic stimulation of central nervous system neurons. PLoS Comput. Biol. 7:e1002022. doi: 10.1371/journal.pcbi.1002022

PubMed Abstract | CrossRef Full Text | Google Scholar

Polson, M. J. R., Barker, A. T., and Freeston, I. L. (1982). Stimulation of nerve trunks with time-varying magnetic fields. Med. Biol. Eng. Comput. 20, 243–244.

PubMed Abstract | Google Scholar

Rall, W. (1962) Electrophysiology of a dendritic neuron model. Biophys. J. 2, 145–167.

PubMed Abstract | Google Scholar

Ram-rez Piscina, L., and Sancho, J. M. (2018). Periodic spiking by a pair of ionic channels. Physica A 505, 345–354. doi: 10.1016/j.physa.2018.03.075

CrossRef Full Text | Google Scholar

Reilly, J. P. (1989). Peripheral nerve stimulation by induced electric currents: exposure to time-varying magnetic fields. Med. Biol. Eng. Comput. 27, 101–110.

PubMed Abstract | Google Scholar

Walsh, V., and Rushworth, M. A (1999). primer of magnetic stimulation as a tool for neuropsychology. Neuropsychologia 37, 125–135.

PubMed Abstract | Google Scholar

Wegner, L. H., and Raschke, K. (1994). Ion channels in the xylem parenchyma of barley roots (a procedure to isolate protoplasts from this tissue and a patch-clamp exploration of salt passageways into xylem vessels). Plant Physiol. 105, 799–813.

PubMed Abstract | Google Scholar

Zhang, X. H., Zhang, Z. X., Li, Y. Z., Ma, H. B., Zhang, H. W., and Liu, W. F. (2014). Biologiacal effects of pulsed electric field based on multi-shelled dielectric model of spherical cell. Trans. Chin. Soc. Agric. Machinery 45, 253–258. doi: 10.6041/j.issn.1000-1298.2014.04.040

CrossRef Full Text | Google Scholar

Zakon, H. H. (2012). Adaptive evolution of voltage-gated sodium channels: the first 800 million years. Proc. Natl. Acad. Sci. U.S.A. 109(Suppl. 1), 10619–10625. doi: 10.1073/pnas.1201884109

PubMed Abstract | CrossRef Full Text | Google Scholar


A.1. Solution of Experiments

The solution of numerical experiments in this paper (in mmol/L): (i) NaCl 139, NaOH 1, KCl 5, CsCl 1, Ca-EGTA 1, MgCl2 1; (ii) NaCl 9, NaOH 1, KCl 120, CsCl 1, Ca-EGTA 1, MgCl2 2; (iii) NaCl 9, NaOH 1, KCl 120, CsCl 1, Mg-EGTA 1, MgCl2 2; (iv) NaOH 1, KCl 40, CsCl 1, Ca-EGTA 1, MgCl2 1, Sorbitol 90; (v) NaOH 1, KCl 1, CsCl 30, Ca-EGTA 1, MgCl2 1, Sorbitol 120; (vi) NaOH 1, KCl 30, CsCl 1, Ca-EGTA 1, MgCl2 1, Sorbitol 120; (vii) NMDG-Cl 139, NMDG-OH 1, KCl 5, CsCl 1, Ca-EGTA 0.001, MgCl2 1.

A.2. Ionic Radius

The ionic radiuses in this paper (pm): Na+ 102, Mg2+ 72, Cl 181, K+ 138, Ca2+ 100, Cs+ 167, OH 137.

A.3. Parameters of the Channel Model

All the parameters of the channel model in this paper are recorded in Tables A1A4 in Appendix. The parameters unassigned in appendix, defaults to 0.


Table 1. Filter state parameters.


Table 2. Properties of gate and sensor.


Table 3. The values of sensitivity parameter sh, i.


Table 4. Standard selective permeability ϱh,i(ms-1) of filter.

Keywords: membrane, ion channel, gate, filter, ionic concentration and flux, electrophysiology, neurosciences

Citation: Wang VQ and Liu S (2019) A General Model of Ion Passive Transmembrane Transport Based on Ionic Concentration. Front. Comput. Neurosci. 12:110. doi: 10.3389/fncom.2018.00110

Received: 28 August 2018; Accepted: 24 December 2018;
Published: 22 January 2019.

Edited by:

Yu-Guo Yu, Fudan University, China

Reviewed by:

Rukiye Uzun, Bülent Ecevit University, Turkey
Jun Ma, Lanzhou University of Technology, China

Copyright © 2019 Wang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Liu Shenquan,