Next Article in Journal
Stimuli-Responsive Gels
Next Article in Special Issue
Advancements in Gel Science—A Special Issue in the Memory of Toyoichi Tanaka
Previous Article in Journal
Hydrogels for Cell Delivery
Previous Article in Special Issue
Dynamic Behaviors of Solvent Molecules Restricted in Poly (Acryl Amide) Gels Analyzed by Dielectric and Diffusion NMR Spectroscopy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Heterogeneous Gelation Dynamics and Their Application to Blood Coagulation

1
Division of Molecular Science, Graduate School of Science and Technology, Gunma University, Kiryu, Gunma 376-8515, Japan
2
Division of Pure and Applied Science, Graduate School of Science and Technology, Gunma University, Kiryu, Gunma 376-8515, Japan
*
Author to whom correspondence should be addressed.
Submission received: 11 May 2018 / Revised: 3 July 2018 / Accepted: 5 July 2018 / Published: 9 July 2018

Abstract

:
We present a scaling model based on a moving boundary picture to describe heterogeneous gelation dynamics. The dynamics of gelation induced by different gelation mechanisms is expressed by the scaled equation for the time taken for development of the gel layer with a few kinetic coefficients characterizing the system. The physical meaning obtained by the analysis for a simple boundary condition from the standpoint of the phase transition shows that the time development of the gelation layer depends on whether the dynamics of the order parameter expressing the gelation of the polymer solution is fast or slow compared with the diffusion of the gelators in the heterogeneous gelation. The analytical method is used to understand the coagulation of blood from various animals. An experiment using systems with plasma coagulation occurring at interfaces with calcium chloride solution and with packed erythrocytes is performed to provide the data for model fitting and it is clarified that a few key kinetic coefficients in plasma coagulation can be estimated from the analysis of gelation dynamics.

Graphical Abstract

1. Introduction

Commonly used gels are prepared using chemical reagents while mixing, by lowering the temperature and by the irradiation of high-energy electromagnetic waves, such as an electron beam and UV, where the cross-linking sites of polymer networks are dispersed randomly in gels, resulting in a macroscopically isotropic and homogeneous structure. This is because the maximum entropy of easily deformed polymer segments between cross-linking sites is obtained when the number of conformational states is maximum. Therefore, to prepare anisotropic gels some asymmetry, such as an asymmetric external force on the gel, should be required. To date, most studies performed on gelation have been on homogeneous systems [1,2,3]. In industry, however, gels are frequently prepared at an interface between two phases. Following the classification of chemical reactions into homogeneous and heterogeneous reactions, we refer to the gelation occurring in a homogeneous phase as homogeneous gelation and the gelation involving reactions at interfaces between different phases as heterogeneous gelation. In 1954, Thiele reported the fabrication of anisotropic beads by dripping an alginate solution in a bath of calcium chloride aqueous solution [4]. Thiele and coworkers prepared anisotropic channel-like alginate gels induced by diffusion of various multivalent cations [5]. These are among the first examples of heterogeneous gelation, on which many studies have now been reported [6,7,8,9,10,11,12,13,14,15,16,17]. The most important gels yielded by heterogeneous gelation are gels in parts of the human body. They are anisotropic, with the anisotropy often having a gradient, such as blood vessel walls [18]. Although the mechanism of biosynthesis is rather complex, most biotissues are yielded from a surface or a point and the growth is directional [19]. Therefore, understanding heterogeneous gelation is important from both industrial and biological viewpoints.
Among the systems that undergo heterogeneous gelation, polymer solutions in contact with solutions containing gelators such as cross-linkers have many promising applications. As a typical example of an anisotropic gel whose structure has been studied is a cylindrical gel prepared by the dialysis of sodium alginate aqueous solution in calcium chloride aqueous solution [4,5,6,7,8,9,10,20]. Maki et al., performed a small-angle X-ray scattering study on the structure of the alginate gel and showed that the polymer orientation is perpendicular to the flow of the gelators [20]. They showed that the gel is highly anisotropic near the surface (dialysis tube) but becomes isotropic toward the center. It has been shown for many combinations of polymers and gelators that the orientation is always similar to that observed for the alginate [14,15,21]. By changing the boundary conditions, we can prepare anisotropic gels with various shapes such as spherical gels and fiber-like gels [8]. The birefringence distribution of the gels is consistent with the orientation of polymer molecules perpendicular to the flow of the gelators [8]. In these systems, polymer molecules cannot permeate into the gel and the gelators selectively diffuse into the already yielded gel, resulting in the cross-linking of polymers at the reaction front. Since the polymers should have a high affinity with the gelators, they are in contact with the interface over a large contact area, which results in molecular orientation. Study of the microrheological events occurring on polymer molecules in the boundary layer at the sol-gel interface is challenging, although such studies should be performed in the near future.
For applications, the present approach to yielding anisotropy has an advantage of gelation in a self-organized manner under mild conditions. Several elaborate gels having unique properties have been developed. Gong’s group at Hokkaido University used a double network gelation involving the contact of synthetic polymer solutions with gelator solutions to prepare a variety of tough anisotropic artificial substitutes such as artificial cartilage [22,23]. Furusawa et al. prepared epithelial lumen-like engineered tissues by the contact of collagen solutions with buffer solutions having a high pH due to spinodal decomposition [24]. Konno et al., prepared a multilayer cylindrical gel consisting of an inner isotropic gel layer and an outer highly anisotropic gel layer [25,26]. The release of drugs from gel particles having such multilayers may exhibit different behavior from the common exponential behavior observed for homogeneous particles and used for stepwise or constant drug release, similarly to capsules or spheres having layer-by-layer wall membranes [27].
On the other hand, the results of studies on the dynamics of heterogeneous gelation have not been organized and almost no applications of such gelation dynamics have been reported, although basic data on the gelation dynamics for simple polymer systems that undergo gelation in contact with gelator solutions have been accumulated [28,29,30,31,32,33]. Measurable physical quantities characterizing the gelation process are the gel fraction and the molecular weight distribution in the sol fraction in homogeneous gelation [1], whereas they are the gel volume and the degree of orientation of the gel in heterogeneous gelation [8]. In most cases, gel layers are formed at the interface between the polymer solution and the gelator solution and extended into the polymer solution while keeping the sol-gel interface clear [8]. Yamamoto et al., developed an irreversible thermodynamic theory focusing on the motion of the sol-gel boundary (moving boundary (MB) picture) [28,29].
In this article, we attempt to reorganize the experimental results on the basis of the MB picture, focusing on the role of the gelator. Then we propose a way of classifying the previously reported systems into several types as a first step toward understanding the characteristics of heterogeneous gelation dynamics. The gelation dynamics are expressed using universal scaling equations corresponding to the gelation mechanism, such as gelation induced by the inflow of cross-linkers (Case 1), by solvent exchange (Case 2), by the inflow of catalysts (Case 3), by the exchange of solutes having very different diffusion constants (Case 4) and by nucleation at low supersaturation (Case 5) (see Section 2.3 for details). In Section 3, we discuss the physical meaning of isotropic and anisotropic gelation from the standpoint of phase separation. By fitting data for gelation dynamics whose gelation mechanism is unknown to the universal scaling equation for Cases 1–5, we can identify the type of the gelation of the system. Furthermore, we can determine kinetic coefficients from the fitting parameters. The second aim of this article is to apply the analytical method to one of the most important biomedical processes, blood coagulation. Although some properties of the blood of patients can be obtained from conventional biochemical tests and the time required for blood coagulation, we have few means of estimating kinetic properties in blood coagulation, which is part of the missing link between blood properties and disorder in blood coagulation. Analysis of the dynamics of blood coagulation by fitting the data to the universal scaling equation enables us to extract information on kinetic coefficients in the process that cannot be obtained by static measurements. In Section 4, we describe the application of the analytical method for gelation dynamics to two model systems of blood coagulation to determine several key kinetic coefficients relating to blood coagulation.

2. Experimental Results on Gelation Dynamics and Classification of Systems by Moving Boundary Picture

We propose a generalized model of heterogeneous gelation and show the observed gelation dynamics of chitosan solution [30] as a typical example in Section 2.1. Then we demonstrate a theoretical analysis of the dynamics based on the MB picture in Section 2.2. In Section 2.3 we discuss the factors which determine the characteristics of the dynamics in Section 2.2 and propose a way of classifying of various systems according to the key factors.

2.1. Gelation Dynamics of Chitosan Solution in Contact with Solutions with High pH

The gel growth behaviors from the interface of polymer solutions with various types of gelator solutions have been observed. In Figure 1 we show an illustration of one of the simplest cases, the one-dimensional growth of a gel in a polymer solution cell in contact with a gelator solution bath (left-hand side of the cell). This model was proposed for the analysis of gelation dynamics of chitosan solution induced by a change of pH. Figure 2 shows a typical observed time course of the gel layer thickness X(t) induced at the interface of chitosan solution with NaOH solution [25]. Chitosan solution is soluble at a low pH and forms a gel at a neutral pH by hydrogen bonding. Under the geometry shown in Figure 1, the part of the chitosan solution where the pH changes from a low pH to a neutral pH is transformed to a gel. The initial process appears to be expressed by square-root behavior X ~ t , as shown in the inset of Figure 2, suggesting diffusion-limited dynamics. Here we modify the model shown in Figure 1 and generalize it by including various cases in which A and B have different roles other than those of NaOH solution and chitosan solution, respectively. For example, the outflow of B may or may not be involved in gelation and gelation can occur with or without the consumption of A in the generalized model. Similar initial behaviors have, however, been observed for different types of gelator solutions under various geometries, whereas the late-stage behaviors were different from each other [28,29,30].

2.2. Moving Boundary Picture for Gelation Dynamics

The MB picture proposed by Yamamoto et al. [28,29] describes the gelation dynamics in terms of the sol-gel interface motion caused by the inflow of the gelator. The idea of the MB picture will be briefly explained through the gelation dynamics of chitosan solution. The dynamics introduced in Section 2.1 is analyzed as follows. In the system, the chitosan solution is enclosed in an L (length) × W (width) × D (depth) rectangular-solid space as shown in Figure 1 and is immersed in NaOH solution. The chitosan solution is in contact with the NaOH solution at a W × D side surface (the contact surface) [30].
The motion of the boundary x = X ( t ) is derived from the following assumptions.
(A) The sodium ions flow into the sol part and the acetate ions flow out from the sol part to the NaOH solution through the gel layer. The neutralization caused by the flows instantly results in the cross-linking of the inner chitosan solution to produce a new gel layer.
(B) The gel layer does not capture the inflow sodium-ions and the outflow acetate-ions by acting as a sink and the flows change so slowly that they can be considered to be in a steady state; all the gelators inflowing from the gelator solution arrive at the inner polymer solution through the gel layer to realize a steady state.
(C) The NaOH solution, gel and inner sol are in local thermodynamic equilibrium at the boundaries.
The concentrations of sodium ions in the NaOH solution, in the inner sol, in the gel layer, at the interface between the gel and the NaOH solution at the position x = x in the gel layer and at the interface between the gel and the inner sol in the gel layer are respectively denoted by ρ s , ρ 0 , ρ ( x ) , ρ s and ρ 0 . The concentrations of acetate ions at the corresponding positions are respectively denoted by C s , C 0 , C ( x ) , C s and C 0 . The chemical potentials of the sodium ion in the NaOH solution, in the inner sol and in the gel layer are respectively denoted by μ N a O H s ( ρ s ) , μ N a O H 0 ( ρ 0 ) and μ N a O H ( ρ ) and those of the acetate ion are denoted by μ A c s ( C s ) , μ A c 0 ( C 0 ) and μ A c ( C ) .
The flux of the inflow sodium-ion in the gel layer is expressed as
j N a O H ( x ) = j N a O H ( x ) e x
where j N a O H is the sodium ion flux density and is expressed in terms of the concentration ρ ( x ) and velocity v N a O H ( x ) of the sodium ion as j N a O H ( x ) = ρ ( x ) v N a O H ( x ) . The unit vector along the x -axis is denoted by e x . According to Fick’s law, the velocity of the sodium ion at x is given by
v N a O H ( x ) = k N a O H μ N a O H ( ρ ( x ) ) x
where k N a O H is the mobility of the sodium ion. Therefore, we have
j N a O H ( x ) = k N a O H ρ ( x ) μ N a O H ( ρ ( x ) ) x
In a similar way, we have the flux of the outflow acetate-ion as
j A c ( x ) = j A c ( x ) e x
with
j A c ( x ) = k A c C ( x ) μ A c ( C ( x ) ) x
In the above, k A c is the mobility of the acetate-ion.
Let a new gel layer thickness d X be produced during period d t by neutralization caused by the sodium ion inflow. Assumption (A) gives the following relationship:
1 ρ G ( C 0 ( t ) ) j N a o H ( X ) W D d t = W D d X
where ρ G ( C 0 ) is the number of sodium ions required to neutralize a unit volume of the inner chitosan solution with acetate ion concentration C 0 and is reasonably assumed to be proportional to C 0 ,
ρ G ( C 0 ) = α C 0
where α is a positive constant. Note that C 0 is a function of the immersion time t . Hence, we have the time development equation for the gel thickness:
d X d t = 1 α C 0 ( t ) j N a o H ( X ( t ) )
Assumption (B) requires the following relationships:
div j N a O H = j N a O H x = k N a O H x [ ρ ( x ) μ N a O H ( ρ ( x ) ) x ] = 0
div j A c = j A c x = k A c x [ C ( x ) μ A c ( C ( x ) ) x ] = 0
Integrating the differential Equations (9) and (10), we have the steady-state flows of sodium ions and acetate ions:
j N a O H = k N a O H g N a O H ( ρ s ) g N a O H ( ρ 0 ) X
j A c = k A c g A c ( C s ) g A c ( C 0 ) X
where g N a O H and g A c are the pressures of sodium ions and acetate ions respectively given by
g N a O H ( ρ ) = μ N a O H ( ρ ) ρ f N a O H ( ρ )
g A c ( C ) = μ A c ( C ) C f A c ( C )
Here, the free energies f N a O H and f A c per unit volume satisfy the relationships μ N a O H ( ρ ) = f N a O H ( ρ ) / ρ and μ A c ( ρ ) = f A c ( ρ ) / ρ .
Assumption (C) gives the chemical potential balance
{ μ N a O H ( ρ s ) = μ N a O H s ( ρ s ) μ N a O H ( ρ 0 ) = μ N a O H 0 ( ρ 0 ) μ A c ( C s ) = μ A c s ( C s ) μ A c ( C 0 ) = μ A c 0 ( C 0 )
Since the concentration of acetic acid is not so large, the chemical potential balance can be rewritten by assuming continuity of the concentration, C s C s and C 0 C 0 . In the immersion solution, the concentration of acetic acid is very small, C s 0 , because the volume of the immersion solution is very large. Hence, C s C s 0 . Therefore, the flux of the acetate-ion flow is given by
j A c = k A c g A c ( 0 ) g A c ( C 0 ) X
In the dilute limit of acetate ions in the gel, the pressure is expressed by g A c ( C ) = k B T C and the flux is given by
j A c = β C 0 X
where β = k A c k B T . The acetate ion concentration in the sol part C 0 decreases with increasing immersion time t since acetate ions flow out from the sol part and the time development of the acetate ion concentration is given by
V ( X ) d C 0 ( t ) d t = W D j A c
where V ( X ) = W D ( L X ) is the volume of the sol part. Using Equations (15) and (16), we have
d C 0 ( t ) d t = β C 0 ( t ) L X ( t )
Treating the sodium ion flow in the same manner as the acetate ion flow, we have
j N a O H = γ ρ s X
where γ = k N a O H k B T .
Solving the simultaneous equations, Equations (8), (19) and (20), with the initial condition C 0 ( 0 ) = C i n and X ( 0 ) = 0 and introducing the scaled variables X ˜ = X / L and t ˜ = t / L 2 , we have the scaled equation
ζ ˜ ( X ˜ , β K i n ) = K i n t ˜
with
ζ ˜ ( X ˜ , β K i n ) = 0 X ˜ u ˜ 1 β K i n ln ( 1 u ˜ ) d u ˜
and
K i n = γ α ρ s C i n
From the expansion ζ ˜ ( X ˜ , β K i n ) = 1 2 X ˜ 2 + O ( X ˜ 3 ) around X ˜ = 0 , we derive the initial-stage ( X ˜ 0 ) behavior as
1 2 X ˜ 2 ( t ˜ ) = K i n t ˜
This equation indicates square-root behavior in which the gel thickness increases proportionally to the square of the immersion time,
X = 2 K i n t
The square-root behavior is a characteristic feature of the diffusion-limited dynamics.
The experimental results shown in Figure 2 are analyzed using Equations (21) and (22). In the analysis, K i n and β are the fitting parameters. The results are plotted according to the equations in Figure 3. The slope Kin is proportional to the NaOH concentration in the immersion solution, inversely proportional to the acetic acid concentration in the chitosan solution and independent of the chitosan concentration, as predicted by Equation (23) [30] (not shown). Therefore, the time course of the gel thickness is fully explained by the MB picture. The observed gelation behavior is expressed by the scaled Equation (21) with system-dependent coefficients K i n and β . Thus, we can obtain information on the kinetic coefficients γ and kAc of the system from the fitting parameters K i n and β by comparing the experimental data with the theoretical equation.

2.3. Classification of Gelation Dynamics

In the gelation of chitosan solution, there are three types of characters having different roles. The first character is chitosan molecules, which are the element polymers constituting the gel. The second character is sodium ions, which are the gelator and are consumed to produce a gel layer. The third character is acetate ions, which are the gelation inhibitor.
Focusing on the gelator and the gelation inhibitor, we classify the gelation system into several types depending on the gelation mechanism. The simplest case is Case 1, in which cross-linking occurs simply by the inflow of cross-linkers as gelators; the gelation inhibitors are absent and the gelators are consumed to produce a gel layer. In Figure 1, the gelator A is involved in gelation whereas B is not involved in gelation. An example of Case 1 is a system where sodium alginate aqueous solution is in contact with calcium chloride aqueous solution (alginate/Ca2+ system) [4,5,6,7,8,9,10,20], in which calcium ions are the gelators. In Case 2, gelators flow into the element polymer solution and the gelation inhibitors flow out. Examples of Case 2 are polymers such as chitosan and collagen that undergo gelation via the formation of hydrogen bonds induced by a change in pH resulting from contact with a high pH solution. When chitosan or collagen in acetic acid solution with a low pH comes in contact with an aqueous solution of NaOH with a high or medium pH, anisotropic gels are prepared [15,30]. The gelator may be a catalyst, as in Case 3, making it necessary to consider the repeated use of the catalyst. In Case 3, although the inhibitors are absent, as in Case 1, the gelators are not consumed in the gelation. An example of Case 3 is the system where gelatin aqueous solution is in contact with transglutaminase aqueous solution [34].
Note that transient viscoelastic change occurs without any reactions when both solutes in liquid phases have very different diffusion constants [35]. We define this case as Case 4. In Case 4, the inhibitors are absent, as in Case 1. However, the gelator makes no links between the polymers. An example of Case 4 is the system where a high-molecular-weight DNA aqueous solution sandwiched between a pair of cover glasses is immersed in a low-molecular-weight DNA aqueous solution. In this case, the initial inflow of low-molecular-weight DNA causes the high-molecular-weight DNA solution to transiently form an anisotropic gel-like substance owning to the excluded volume effect and then the high-molecular-weight DNA diffuses to the immersion solution, finally resulting in a homogeneous solution. The gelation dynamics of the above-mentioned typical diffusion-limited systems were derived in Section 2.2 and are summarized in Table 1.
Case 5 is similar to the crystal growth from solutions at low supersaturation. Assume that the free energy of a polymer solution has a double minimum with different values. The lower and higher values correspond to gel phase and sol phase, respectively. If a gel region is nucleated by contact with a solid phase, then the polymer solution gels and the front line of the gel moves forward. In Case 5, neither gelators nor gelation inhibitors are present. The gelation is not diffusion-limited in this case. An example of Case 5 is the system where plasma is in contact with packed erythrocytes, which is discussed in detail in Section 4.2. The cases given above are not exhaustive but are those for which examples are reported [28,29,30,31,32,33,34,35].

3. Theoretical View of Heterogeneous Gelation from the Standpoint of Phase Transition Dynamics

Here we regard gelation as a phase transition and discuss the gelation dynamics in the context of the phase transition dynamics [36,37,38] for the polymer solution in contact with the gelator solution. Let us discuss the gelation system corresponding to Case 1 as the simplest case. The state of the polymer solution depends on the concentration ρ of gelators in the polymer solution. The state of the polymer solution is described in terms of the order parameter ψ , which expresses the degree of gelation; ψ > 0 when the polymer solution is in the gel state and ψ = 0 when it is in the sol state. Then, the local free energy of the polymer solution is a function of ρ and ψ ,
f = f ( ρ , ψ )
Let us consider a one-dimensional gel growth system, as shown in Figure 1. The free-energy functional expressing the whole of the polymer solution is given by
F ( { ρ } , { ψ } ) = [ 1 2 κ ρ ( ρ x ) 2 + 1 2 κ ψ ( ψ x ) 2 + f ( ρ ( x ) , ψ ( x ) ) ] d x
where the positive constants κ ρ and κ ψ express the increase in free energy with respect to the space inhomogeneity of the gelator concentration and the degree of gelation, respectively. Since the concentration of gelators is a conserved quantity, the time development equation for the gelator is given by
ρ ( x , t ) t = div [ Γ ρ δ F δ ρ ( x , t ) ]
where t denotes the elapsed time and Γ ρ is the kinetic coefficient for ρ . Since ψ is an order parameter and is not a conserved quantity, the time development equation for ψ is given by
ψ ( x , t ) t = Γ ψ δ F δ ψ ( x , t )
where Γ ψ is the kinetic coefficient for ψ . The boundary conditions
μ G ( ρ ( 0 , t ) , ψ ( 0 , t ) ) = μ S
is imposed, where μ G ( ρ ( 0 , t ) , ψ ( 0 , t ) ) and μ S are, respectively, the chemical potential of the gelator in the polymer solution at the liquid-liquid contact interface and that in the gelator solution. In terms of F , μ G is generally given by
μ G ( ρ ( x , t ) , ψ ( x , t ) ) = δ F δ ρ ( x , t )
One of the most characteristic features of the heterogeneous gelation process is that the boundary between the sol phase and gel phase is macroscopically clear. Thus, the gel growth dynamics can be visualized by the MB picture. The dynamics shown in Table 1 was derived by means of the MB picture.
Assumptions (A)–(C) in Section 2.2 for the gelation of chitosan solution (Case 2) is modified for Case 1 as follows [28] since the outflow ions are not involved in gelation in Case 1.
(A’) All the inflow gelators arriving at the inner polymer solution instantly cross-link the polymers to produce a new gel layer.
(B’) The gel layer does not capture the inflow gelators by acting as a sink and the gelator flow changes so slowly that it can be considered to be in a steady state; all the gelators inflowing from the gelator solution arrive at the inner polymer solution through the gel layer to realize a steady state.
(C’) The gelator solution, gel and inner polymer solution are in local thermodynamic equilibrium at the boundary.
According to the above three assumptions, in terms of the concentration ρ ( x ) and chemical potential μ G ( ρ ( x ) ) of the gelator in the gel layer, the gelator flow density j ( x ) is given by
j ( x ) = k ρ ( x ) μ G ( ρ ( x ) ) x
where k is the mobility of the gelator and we modify Equations (8), (19) and (20) governing the time development equations for the sol-gel boundary X ( t ) as [22]
{ d X ( t ) d t = j ( X ( t ) ) ρ G j ( X ) = k k B T ρ s X
where ρ G is the number of gelators required to produce a unit volume of gel. Note that ρ G is constant and the equation corresponding to Equation (18) is absent in Case 1 since gelation inhibitors are absent.
Using the initial condition X ( 0 ) = 0 , we obtain the solution of Equation (33) as
X ( t ) = 2 K t
where K = k k B T . For Case 1, the square root behavior occurs during the entire period, although the behavior only occurs during the initial stage for Case 2.
Let us compare the above results with the phase transition dynamics given by Equations (28) and (29). The MB picture does not refer to the time development of the degree of gelation. Assumption (A’), however, requires that the time development of ψ is very rapid and an equilibrium value that depends on the gelator density ρ is quickly achieved. Thus, the degree of gelation ψ is expressed as
ψ ( x , t ) ψ 0 Θ ( X ( t ) x )
where for simplicity we let the “level function” Θ ( x ) be given by the step function defined by
Θ ( x ) = { 0 x < 0 1 x > 0
and ψ 0 is a positive constant satisfying the equation
f ( ρ G , ψ 0 ) ψ 0 = 0
The gelator concentration in the MB picture indicates the concentration δ ρ ( x ) of the gelator dissolved in the solvent, which is given by
δ ρ ( x ) = ρ ( x ) ρ G
Therefore, in the context of phase transition dynamics, Equation (32) should be rewritten as
j ( x ) = k δ ρ ( x ) μ G ( ρ ( x ) ) x
Equation (39) requires the relationship
Γ ρ = Γ ρ ( ρ , ψ 0 ) = k δ ρ
for the kinetic coefficient Γ ρ = Γ ρ ( ρ , ψ ) . Assumption (B’) requires that the gelator flow can be regarded as a steady flow. The time dependence of the gelator concentration ρ is very weak and appears only through the time dependence of the boundary X ( t ) ; the boundary condition, that is, Assumption (C’), determines the weak time dependence. Hence, the gelator concentration is given by
ρ ( x , t ) ( ρ G + δ ρ ( x ; X ( t ) ) ) Θ ( X ( t ) x )
The chemical potential μ G ( x ) in Equations (32) and (39) is given by
μ G ( x ) = f ( ρ G + δ ρ ( x ) , ψ 0 ) δ ρ ( x )
To effectively use the description of the phase transition dynamics given by Equations (28) and (29), it is necessary to clarify the function forms of f ( ρ , ψ ) and Γ ρ and the value of the kinetic constant Γ ψ statistical-mechanically and/or experimentally. This requirement will be discussed as a future issue. Here, we discuss a problem independent of the details of f ( ρ , ψ ) and Γ ρ . Equation (33), expressing the gelation dynamics, is valid when the dynamics of ψ ( x , t ) is very fast compared with the gelator diffusion. What kind of dynamics is obtained in the opposite case, in which the dynamics of ψ ( x , t ) is very slow? The answer is obtained from the common phase transition dynamics. In this case, Equations (28) and (29) are rewritten as
f ( ρ e q , 0 ) ρ e q = 0
ψ ( x , t ) t = Γ ψ δ F ψ δ ψ ( x , t )
with
F ψ = [ 1 2 κ ψ ( ψ x ) 2 + f ( ρ e q , ψ ) ] d x
The boundary condition is given by
{ ψ ( 0 , t ) = ψ e q lim x ψ ( x , t ) = 0
where the equilibrium value of the degree of gelation ψ e q is obtained from the equation
f ( ρ e q , ψ e q ) ψ e q = 0
The dynamics described by Equation (44) is limited by the free energy and is regarded as the relaxation from an unstable state to a stable state minimizing the free energy.
When f ( ρ e q , ψ ) is a double-well-type function (both the sol phase ( ψ = 0 ) and the gel phase ( ψ = ψ e q ) locally minimize the local free energy f ( ρ e q , ψ ) ) and the constant κ ψ is small, the solution of Equation (44) is approximately [37,39]
ψ ( x , t ) ψ e q Θ ( X ( t ) x )
In this case, the gelation dynamics is also visualized by the motion of the sol-gel boundary. Therefore, to pay attention to the motion of the sol-gel boundary is valid. We could call the moving boundary picture the analysis method in which we pay attention to the sol-gel boundary motion even in the case of the energy-limited dynamics. The time dependence of X ( t ) , however, differs from the square-root behavior. The dynamics is expected to effectively visualize the gelation dynamics for Case 5.
Finally, we discuss how the anisotropy is taken into account. The presence of a gelator concentration gradient ρ ( x ) / x during gelation is one of the characteristics of heterogeneous gelation. Therefore, we assume that the local free energy depends on the concentration gradient and introduce the order parameter ϕ expressing the degree of anisotropy of the polymer solution, where ϕ > 0 when the polymer solution is anisotropic and ϕ = 0 when it is isotropic. The local free energy is given by
f = f ( ρ , ρ x , ψ , ϕ )
The free-energy functional is given by
F ( ρ , ρ x , ψ , ϕ ) = [ 1 2 κ ρ ( ρ x ) 2 + 1 2 κ ψ ( ψ x ) 2 + 1 2 κ ϕ ( ϕ x ) 2 + f ( ρ , ρ x , ψ , ϕ ) ] d x
where the positive constant κ ϕ expresses the increase in free energy with respect to the space inhomogeneity of the anisotropy. In addition to Equations (28) and (29), we introduce the time development equation for ϕ as
ϕ ( x , t ) t = Γ ϕ δ F δ ϕ ( x , t )
where Γ ϕ is the kinetic coefficient for ϕ .
Let us discuss the diffusion-limited gelation dynamics with the anisotropy expressed by the MB picture. The gelation dynamics of curdlan solution [28] is in this category. In this case, the time development of the anisotropy is expressed as
ϕ ( x , t ) ϕ 0 Θ ( X ( t ) x )
where
f ( ρ G , δ ρ ( X ) X , 0 , ϕ 0 ) ϕ 0 = 0 .
Equation (53) requires the local free energy to reach a local minimum at a finite value of the degree of anisotropy. Note that global minimization of the local free energy at a finite ϕ is not required. At a concentration gradient | δ ρ ( X ) X | larger than a threshold value, the local free energy is expected to have such a local minimum. The increase in the degree of gelation roughly coincides with the increase in anisotropy. The kinetic coefficient Γ ϕ rapidly drops to zero with increasing ψ since the gelation significantly interferes with the polymer motion in the polymer solution. Thus, the degree of anisotropy is fixed at a finite value ϕ 0 .

4. Application of Scaled Gelation Dynamics to Analysis of Blood Coagulation

One of the most interesting applications of the analysis of heterogeneous gelation dynamics is blood coagulation (gelation). Blood consists of about 45% blood corpuscles, which are mainly erythrocytes and the remaining liquid component, plasma. Fibrinogens, which comprise 7% of blood protein, are the main component of coagulants (gels). Blood coagulation is triggered by the contact of plasma with coagulant factors at cell surfaces or blood vessel walls, as illustrated in Figure 4. Therefore, the main process of macroscopic blood coagulation is the gelation of plasma induced by the contact of plasma and coagulant factors at a cell surface. With a sufficient amount of calcium ions, after complex cascade reactions, the key protein thrombins are activated and fibrinogens are hydrolyzed to form fibrins and then protofibrils, which are the building blocks of the coagulant. The biochemical reaction cascade in blood coagulation has been established by considering a number of simplified homogeneous systems [40]. However, the research on the dynamic aspects of blood coagulation by considering heterogeneous systems, such as a thrombin solution/plasma contact system, began relatively recently [41,42]. Determination of the kinetic coefficients of blood coagulation should help provide the missing link between the biochemical properties and disorders in blood coagulation. The kinetic coefficients can potentially be used as direct indicators of blood coagulability in diagnosis.
It is very difficult to completely replace an in vivo system by an in vitro system for blood testing, because when we draw blood from a body we must add some anticoagulants such as citrates or heparin to the blood. Therefore, it is important to study key systems in which part of the blood coagulation process can be reproduced to obtain useful information on blood. Here we discuss two systems of interest, citrated plasma in contact with calcium chloride aqueous solution (plasma/Ca+) and citrated plasma in contact with packed erythrocytes (plasma/packed erythrocytes).

4.1. Analysis of Citrated Plasma/Ca+ Contact System

In the clinical test to assess blood coagulability, a trigger protein or trigger lipid and calcium ions are mixed with citrated plasma and the clotting time is measured and compared with the standard value. Let us consider heterogeneous gelation as a control clinical testing system; this system consists of citrated plasma in contact with a calcium ion solution separated by a dialysis membrane, as illustrated in Figure 1. There is an inflow of calcium ions and an outflow of citrate ions through the dialysis membrane. When the calcium ion concentration in the plasma exceeds the threshold, gelation occurs. If we replace citrate ions, calcium ions and fibrinogen by acetate ions, sodium ions and chitosan, respectively, this system corresponds to Case 2 [43]. Therefore, the measurement of the gel layer thickness X as a function of time after the liquid-liquid contact t can be compared with the theoretical equation, Equation (T2-1) in Table 1, of
ζ ˜ ( X ˜ , β K i n ) = 0 X ˜ u ˜ 1 β K i n ln ( 1 u ˜ ) d u ˜ = K i n t ˜
where
X ˜ = X / L , t ˜ = t / L 2
and L is the effective length of the cell. Since the blood coagulation dynamics is expressed by the scaled form of Equation (54), the analytical results based on the equation do not depend on the experimental tools. This is an advantage of this method of extracting valuable information on blood.
Recently the entire time course of the gel thickness was compared with the theory for bovine blood and was found to be well expressed by Equation (54) [43]. However, since the blood coagulation behavior of animals of the even-hoof class, such as cattle, is known to be considerably different from that of other animals such as human, swine and horses [44], in this paper we performed the corresponding experiment using blood from a horse, one of the animals of the odd-hoof class. The time course was well expressed by Equation (54) as shown in Figure 5. The scaled time lag k used for fitting was close to the clotting time, which was the time required for macroscopic gelation by mixing (homogeneous gelation), used in clinical testing. Therefore, the analysis is also applicable for the blood of animals other than even-hoof animals. Note that the observed data could only be fit to Case 2 among the five cases. Calcium ions are known to play various roles in blood. They bind plasma proteins that are both involved and not involved in coagulation, each with a different binding constant. For example, fibrinogen has several strong and weak binding sites with calcium ions. Some enzymes are only activated in the presence of a sufficient amount of calcium ions. Thus, blood coagulation is generally complex. The present system is regarded as a simplified one to relate the observed gelation dynamics to Case 2. Here, we further expand the expressions for parameters β and Kin. From the definition, we have
β = α β 0 C i n ρ G 0 + α C i n
β 0 = k C I T k B T
K i n = γ ρ s ρ G 0 + α C i n
γ = kCakBT
where Cin is the initial concentration of citrate in the plasma, ρs is the initial concentration of calcium ions in the calcium chloride aqueous solution, γ and β0 are the diffusion constants of calcium ions and citrate ions in the gel, respectively and ρ G 0 is the calcium ion concentration required for the gelation of a unit volume of plasma. From Equations (56)–(59), we have
C i n = ρ G 0 + γ ρ s 0 K i n
and
C i n = ρ G 0 + β 0 C i n β
Therefore, by measuring the citrate concentration dependence of the parameters β and Kin, we can determine the three physically defined quantities γ, β0 and ρG0 from the plot using Equations (60a) and (60b). According to preliminary experiments, this analysis is also valid for human blood and the kinetic coefficients depends on subjects. Therefore, these parameters could be used for diagnosis as new indicators reflecting blood coagulability.

4.2. Analysis of Plasma/Packed Erythrocyte Contact System

Until recently, erythrocytes were regarded as being unrelated to blood coagulation. However, in the past few decades, several pieces of indirect evidence that erythrocytes have an active function in thrombosis and hemostasis have been found [45]. In particular, in the case of blood flow stagnation, such as in deep vein thrombosis, also called economy-class syndrome, the erythrocyte surface is suspected to be one of the factors that trigger blood coagulation [46]. Deep vein thrombosis is a blood clot that develops within a deep vein in the body, usually in the leg. When part of the blood clot breaks off, enters the bloodstream and blocks one of the blood vessels in the lungs, it leads to complications such as pulmonary embolism [47]. It has been assumed that protein C or protein S, both of which inhibit the coagulation cascade in the normal state, is inactivated, resulting in blood coagulation [48]. Kaibara et al. proposed a new pathway involving the activation of an intrinsic pathway in the coagulation cascade by the protein on the erythrocyte membrane, erythroelastase [46]. Although few continuous studies have been carried out on this pathway, the system involving the contact of erythrocytes and plasma seems to be appropriate for studying whether or not this scheme is consistent with or directly related to deep vein thrombosis. Because of the deformability of erythrocytes, we can easily prepare packed erythrocytes with a flat surface by centrifugation. In this study, we performed a study on gelation dynamics induced in swine plasma in contact with packed erythrocytes. Figure 6 shows the time course of the thickness of a gel grown from a packed erythrocyte surface. After a time lag, the gel thickness increased proportionally to time. The slope was roughly constant, irrespective of the calcium concentration, which means that the gelation behavior is the same above a threshold calcium concentration after gelation is initiated. According to the mechanism proposed by Kaibara et al., a coagulation protein called factor IX is activated on the surface of erythrocytes [46]. Then the coagulation cascade should be successively activated and positive feedback may yield a locally activated space with saturated fibrin molecules similarly to in vivo [49]. The fibrin molecules form a network from the surface of the packed erythrocytes, similarly to in the case of crystal growth, as illustrated in Figure 7. This situation is similar to Case 5. Choosing a quartic function of ψ in the function f ( ρ e q , ψ ) and the boundary condition Equation (46), we obtain the solution of Equation (48) according to Chan [50]. The position of the sol-gel interface is expressed by the linear function of time
X ( t ) = v ( t t 0 )
where t 0 is the lag time required so that the boundary condition is satisfied and v is a positive constant depending on the free-energy difference f ( ρ e q , ψ e q ) f ( ρ e q , 0 ) . Therefore, the observed linear increase in X with time t shown in Figure 6 is explained by Equation (61) and the proportionality coefficient is related to the free-energy difference between the initial state and the equilibrium state. On the other hand, such a gel layer was hardly observed on hardened erythrocytes modified by glutaraldehyde whose plasma membrane was inactivated. Thus, the experimental results are consistent with the mechanism proposed by Kaibara et al. [46]. The rate of increase in the gel layer thickness provides information on the coagulability of the examined blood through the proportionality coefficient v. According to a preliminary experiment, this approach is valid for human blood. The kinetic coefficient determined by the observed heterogeneous gelation dynamics may be related to the coagulability of the blood of subjects and can potentially be used as an indicator of deep vein thrombosis.

5. Conclusions

The gelation dynamics induced by the contact of polymer and gelator solutions was classified into several types and expressed by a system-independent scaling equation and system-dependent coefficients. From the viewpoint of the phase transition, it was clarified that the time taken for the development of the gelation layer depends on whether the dynamics of the order parameter expressing the gelation of the polymer solution is fast or slow compared with the diffusion of the gelators. As a straightforward application of one type of gelation dynamics (Case 2), in vitro blood coagulation was analyzed to find kinetic coefficients characterizing the dynamics. New information on blood from human subjects can be extracted by the analysis of gel growth behavior. Heterogeneous gelation is involved in many industrial processes such as microencapsulation by an insolubilization reaction to prepare artificial eggs and drug delivery carriers. Information on the extent of the reaction is important for controlling the retention and release capability of the core substance. For spherical entities, however, it is difficult to trace the gelation dynamics. The theoretical analysis in Section 2 may also be applicable to these materials to enabling the gel layer thickness to be estimated at a desired time.

6. Materials and Method

To study the plasma/Ca+ contact system, horse blood aseptically collected in equal volumes of Alsever’s solution was purchased from Nippon Bio-test Laboratories Inc. The blood was centrifuged at 1600× g for 10 min at 25 °C to obtain platelet-free plasma (PFP). The plasma was poured in a cell (10 mm in height, 3 mm in width and 20 mm depth) made of poly(methylmethacrylate) from the bottom to the upper end (inlet). Then the inlet was sealed with a dialysis membrane (Sanko Jun-yaku Co., Ltd., Tokyo, Japan) to confine plasma proteins such as fibrinogen in the plasma. The cell was immersed in 40 mL of calcium chloride aqueous solution with various calcium chloride concentrations in the range of 10–100 mM. The motion of the boundary was observed and recorded with a digital camera at 25 °C.
To study the plasma/packed RBC contact system, swine blood provided by Gunma Meat Wholesale Market Co. Ltd. was citrated with 3.2% sodium citrate in a polypropylene tube. The suspension was centrifuged at 1600× g for 10 min at 25 °C to separate the PFP and packed erythrocytes. Hardened erythrocytes were prepared by mixing 20% erythrocyte suspension with 1.6% glutaraldehyde aqueous solution with 40 times the amount of erythrocyte suspension and gently stirring it for 72 h. Then the precipitates were rinsed with MilliQ water five times and physiological saline solution twice. After placing 2 mL of the PFP with the desired concentration of calcium ions on 1 mL of packed erythrocytes in a collagen-coated polystyrene tube, we observed the growth of a turbid layer from the interface over time at 25 °C. The volume of the turbid layer was estimated by image analysis and the gel volume was measured directly; these two volumes were found to be equivalent. Therefore, the growth of the gel layer was traced by observing the thickness of the turbid layer.

Author Contributions

T.D. designed the experimental study and T.Y. built the theory. T.D. and T.Y. wrote the paper.

Funding

This work was partly supported by JSPS KAKENHI Grant Numbers JP24350115, JP15K05241 and JP16H04031.

Acknowledgments

The authors are grateful to Yoshiharu Toyama for his helpful discussion and Ryuta Kurasawa and Hiroki Shinoda at Gunma University for their help with the blood coagulation experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Flory, P.J. Principles of Polymer Chemistry; Cornell University Press: Ithaca, NY, USA, 1953. [Google Scholar]
  2. Weiss, R.G.; Terech, P. Molecular Gels Materials with Self-Assembled Fibrillar Networks; Springer: Berlin, Germany, 2006. [Google Scholar]
  3. Djabourov, M.; Nishinari, K.; Ross-Murphy, S.B. Physical Gels from Biological and Synthetic Polymers; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  4. Thiele, H. Ordered coagulation and gel formation. Discuss. Faraday Soc. 1954, 18, 294–314. [Google Scholar] [CrossRef]
  5. Thiele, H.; Hallich, K. Kapillarstrukturen in ionotropen Gelen. Kolloid Z. 1957, 151, 1–12. [Google Scholar]
  6. Mørch, Y.A.; Donati, I.; Strand, B.L.; Skjåk-Bræk, G. Effect of Ca2+, Ba2+ and Sr2+ on alginate microbeads. Biomacromolecules 2006, 7, 1471–1480. [Google Scholar] [CrossRef] [PubMed]
  7. Thumbs, J.; Kohler, H.H. Capillaries in alginate gel as an example of dissipative structure formation. Chem. Phys. 1996, 208, 9–24. [Google Scholar] [CrossRef]
  8. Dobashi, T.; Yamamoto, T. Anisotropic gel formation induced by dialysis. In Encyclopedia of Biocolloid and Biointerface Science; Ohshima, H., Ed.; Wiley: Hoboken, NJ, USA, 2016. [Google Scholar]
  9. Schuster, E.; Sott, K.; Ström, A.; Altskär, A.; Smisdom, N.; Gebäck, L.N.; Hermansson, A.M. Interplay between flow and diffusion in capillary alginate hydrogels. Soft Matter 2016, 12, 3897–3907. [Google Scholar] [CrossRef] [PubMed]
  10. Despang, F.; Dittrich, R.; Gelinsky, M. Novel biomaterials with parallel aligned pore channels by directed ionotropic gelation of alginate: Mimicking the anisotropic structure of bone tissue. In Advances in Biomimetics; George, A., Ed.; InTech: London, UK, 2011; Chapter 17; pp. 349–372. [Google Scholar]
  11. Dobashi, T.; Nobe, M.; Yoshihara, H.; Yamamoto, T.; Konno, A. Liquid crystalline gel with refractive index gradient of curdlan. Langmuir 2004, 20, 6530–6534. [Google Scholar] [CrossRef] [PubMed]
  12. Narita, T.; Tokita, M. Liesegang pattern formation in κ-carrageenan gel. Langmuir 2006, 22, 349–352. [Google Scholar] [CrossRef] [PubMed]
  13. Rivas-Araiza, R.; Alcouffe, P.; Rochas, C.; Montembault, A.; David, L. Micron range morphology of physical chitosan hydrogels. Langmuir 2010, 26, 17495–17504. [Google Scholar] [CrossRef] [PubMed]
  14. Mredha, M.T.I.; Zhang, X.; Nonoyama, T.; Nakajima, T.; Kurokawa, T.; Takagi, Y.; Gong, J.P. Swim bladder collagen forms hydrogel with macroscopic superstructure by diffusion induced fast gelation. J. Mater. Chem. B 2015, 3, 7658–7666. [Google Scholar] [CrossRef] [Green Version]
  15. Furusawa, K.; Sato, S.; Masumoto, J.; Hanazaki, Y.; Maki, Y.; Dobashi, T.; Yamamoto, T.; Fukui, A.; Sasaki, N. Studies on the formation mechanism and the structure of the anisotropic collagen gel prepared by dialysis-induced anisotropic gelation. Biomacromolecules 2012, 13, 29–30. [Google Scholar] [CrossRef] [PubMed]
  16. Mikkelsen, A.; Elgsaeter, A. Density ditribtuion of calcium-induced alginate gels. A numerical study. Biopolymers 1994, 3, 17–41. [Google Scholar]
  17. Miyamoto, Y.; Kaysser, W.A.; Rabin, B.H.; Kawasaki, A.; Ford, R.G. (Eds.) Functionally Graded Materials: Design, Processing and Applications; Springer: Berlin, Germany, 1999. [Google Scholar]
  18. Oka, S. Cardiovascular Hemorheology; Cambridge University Press: Cambridge, UK, 1981. [Google Scholar]
  19. Alberts, B.; Johnson, A.; Lewis, J.; Raff, M.; Roberts, K.; Walter, P. Molecular Biology of the Cell, 5th ed.; Garland Science: New York, NY, USA, 2008. [Google Scholar]
  20. Maki, Y.; Ito, K.; Hosoya, N.; Yoneyama, C.; Furusawa, K.; Yamamoto, T.; Dobashi, T.; Sugimoto, Y.; Wakabayashi, K. Anisotropic structure of calcium-induced alginate gels by optical and small-angle X-ray scattering measurements. Biomacromolecules 2011, 12, 2145–2152. [Google Scholar] [CrossRef] [PubMed]
  21. Maki, Y.; Furusawa, K.; Yasuraoka, S.; Okamura, H.; Hosoya, N.; Sunaga, M.; Dobashi, T.; Sugimoto, Y.; Wakabayashi, K. Universality and specificity in molecular orientation in anisotropic diffusion method. Carbohydr. Polym. 2014, 108, 118–126. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, W.; Furukawa, H.; Gong, J.P. Highly extensible double-network gels with self-assembling anisotropic structure. Adv. Mater. 2008, 20, 4499–4503. [Google Scholar] [CrossRef]
  23. Mredha, M.T.I.; Kitamura, N.; Nonoyama, T.; Wada, S.; Goto, K.; Zhang, X.; Nakajima, T.; Kurokawa, T.; Takagi, Y.; Yasuda, K.; et al. Anisotropic tough double network hydrogel from fish collagen and its spontaneous in vivo bonding to bone. Biomaterials 2017, 132, 85–95. [Google Scholar] [CrossRef] [PubMed]
  24. Furusawa, K.; Mizutani, T.; Machino, H.; Yahata, S.; Fukui, A.; Sasaki, N. Application of multichannel collagen gels in construction of epithelial lumen-like engineered tissues. ACS Biomater. Sci. Eng. 2015, 1, 539–548. [Google Scholar] [CrossRef]
  25. Konno, A.; Tsubouchi, M. Gel formation of curdlan. Kinran Tanki Daigaku Kenkyushi 1998, 29, 89–95. [Google Scholar]
  26. Maki, Y.; Furusawa, K.; Dobashi, T.; Sugimoto, Y.; Wakabayashi, K. Small-angle X-ray and light scattering analysis of multi-layered curdlan gels prepared by a diffusion method. Carbohydr. Polym. 2017, 155, 136–145. [Google Scholar] [CrossRef] [PubMed]
  27. Ariga, K.; Lvov, Y.M.; Kawakami, K.; Ji, Q.; Hill, J.P. Layer-by-layer self-assembled shells for drug delivery. Adv. Drug Deliv. Rev. 2011, 63, 762–771. [Google Scholar] [CrossRef] [PubMed]
  28. Nobe, M.; Dobashi, T.; Yamamoto, T. Dynamics in dialysis process for liquid crystalline gel formation. Langmuir 2005, 21, 8155–8160. [Google Scholar] [CrossRef] [PubMed]
  29. Yamamoto, T.; Tomita, N.; Maki, Y.; Dobashi, T. Dynamics in the process of formation of anisotropic chitosan hydrogel. J. Phys. Chem. B 2010, 114, 10002–10009. [Google Scholar] [CrossRef] [PubMed]
  30. Rokugawa, I.; Tomita, N.; Dobashi, T.; Yamamoto, T. One-dimensional growth of hydrogel by a contact of chitosan solution with high-pH solution. Soft Mater. 2014, 12, 36–41. [Google Scholar] [CrossRef]
  31. Skjåk-Bræk, G.; Grasdalen, H.; Smisrød, O. Inhomogeneous polysaccharide ionic gels. Carbohydr. Polym. 1989, 10, 31–54. [Google Scholar] [CrossRef]
  32. Furusawa, K.; Minamisawa, Y.; Dobashi, T.; Yamamoto, T. Dynamics of liquid crystalline gelation of DNA. J. Phys. Chem. B 2007, 111, 14423–14430. [Google Scholar] [CrossRef] [PubMed]
  33. Furusawa, K.; Narazaki, Y.; Tomita, N.; Dobashi, T.; Sasaki, N.; Yamamoto, T. Effect of pH on anisotropic gelation of DNA induced by aluminum cations. J. Phys. Chem. B 2010, 114, 13923–13932. [Google Scholar] [CrossRef] [PubMed]
  34. Yamamoto, T.; Kakinoki, K.; Maki, Y.; Dobashi, T. Gelation and orientation dynamics of protein solution induced by enzyme solution. In Proceedings of the JPS 2016 Autumn Meeting, Kanazawa, Japan, 13 September 2016. 13aBE3. [Google Scholar]
  35. Dobashi, T.; Maki, Y.; Furusawa, K.; Yamamoto, T. Anisotropic Structure Formation Induced by Liquid-Liquid Phase Contact and Diffusion. In Proceedings of the IUMRS-ICAM 2017, Kyoto, Japan, 31 August 2017. C6-I31-003. [Google Scholar]
  36. Hohenberg, P.C.; Halperin, B.I. Theory of dynamic critical phenomena. Rev. Mod. Phys. 1977, 49, 437–479. [Google Scholar] [CrossRef]
  37. Gunton, J.D.; Miguel, M.S.; Sahni, P.S. The Dynamics of First-order Phase Transitions. In Phase Transitions and Critical Phenomena, 8th ed.; Domb, C., Lebowitz, J.L., Eds.; Academic Press: London, UK; New York, NY, USA, 1983; pp. 267–466. [Google Scholar]
  38. Hohenberg, P.C.; Krekhov, A.P. An introduction to the Ginzburg-Landau theory of phase transitions and nonequilibirm patterns. Phys. Rep. 2015, 572, 1–42. [Google Scholar] [CrossRef]
  39. Allen, S.M.; Cahn, J.W. A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening. Acta Metall. 1979, 27, 1085–1095. [Google Scholar] [CrossRef]
  40. Voet, D.; Voet, J.G. Fundamentals of Biochemistry, Life at the Molecular Level, 4th ed.; Wiley: Hoboken, NJ, USA, 2013. [Google Scholar]
  41. Ataullakhanov, F.I.; Guria, G.T.; Sarbash, V.I.; Volkova, R.I. Spatiotemporal dynamics of clotting and pattern formation in human blood. Biochim. Biophy. Acta 1998, 1425, 453–468. [Google Scholar] [CrossRef]
  42. Zhalyalov, A.S.; Panteleev, M.A.; Gracheva, M.A.; Ataullakhanov, F.I.; Shibeko, A.M. Co-ordinated spatial propagation of blood plasma clotting and fibrinolytic fronts. PLoS ONE 2017, 12, e0180668. [Google Scholar] [CrossRef] [PubMed]
  43. Shida, N.; Kurasawa, R.; Maki, Y.; Toyama, Y.; Dobashi, T.; Yamamoto, T. Coagulation of plasma induced by a contact with calcium chloride solution. Soft Matter 2016, 12, 9471–9476. [Google Scholar] [CrossRef] [PubMed]
  44. Kaibara, M.; Shinozaki, T.; Kita, R.; Iwata, H.; Ujiie, H.; Sasaki, K.; Li, J.Y.; Sawasaki, T.; Ogawa, H. Analysis of coagulation of blood in different animal species with special reference to procoagulant activity of red blood cell. J. Jpn. Soc. Biorheol. 2006, 20, 35–43. [Google Scholar]
  45. Litvinov, I.; Weisel, J.W. Role of red blood cells in haemostasis and thrombosis. ISBT Sci. Ser. 2017, 12, 176–183. [Google Scholar] [CrossRef] [PubMed]
  46. Iwata, I.; Kaibara, M.; Dohmae, N.; Takio, K.; Himeno, R.; Kawakami, S. Purification, identification and characterization of elastase on erythrocyte membrane as factor IX-activating enzyme. Biochem. Biophys. Res. Commun. 2004, 316, 65–70. [Google Scholar] [CrossRef] [PubMed]
  47. Signorelli, S.S.; Ferrante, M.; Gaudio, A.; Fiore, V. Deep vein thrombosis related to environment. Mol. Med. Rep. 2017, 15, 3445–3448. [Google Scholar] [CrossRef] [PubMed]
  48. Tapson, V.F. Acute pulmonary embolism. N. Eng. J. Med. 2008, 358, 1037–1052. [Google Scholar] [CrossRef] [PubMed]
  49. Ieko, M. Characteristics of various anticoagulants and evaluation methods for risk of bleeding and thrombosis. Jpn. J. Electrocardiol. 2014, 34, 149–156. [Google Scholar]
  50. Chan, S.-K. Steady-state kinetics of diffusion less first order phase transition. J. Chem. Phys. 1977, 67, 5755–5762. [Google Scholar] [CrossRef]
Figure 1. Illustration of one-dimensional gel layer growth from the contact surface between polymer solution enclosed in an L (length) × W (width) × D (depth) rectangular cell in an immersing gelator solution (left-hand side of the cell) induced by inflow of A and outflow of B. The x -axis is chosen to be perpendicular to the contact surface and is oriented in the direction from the immersing gelator solution to the polymer solution. The origin of the x -axis is chosen at the contact interface. X(t) denotes the gel layer thickness at immersion time t. In the gelation of chitosan solution, the gelator solution and polymer solution are NaOH aqueous solution and chitosan in acetic acid aqueous solution, respectively.
Figure 1. Illustration of one-dimensional gel layer growth from the contact surface between polymer solution enclosed in an L (length) × W (width) × D (depth) rectangular cell in an immersing gelator solution (left-hand side of the cell) induced by inflow of A and outflow of B. The x -axis is chosen to be perpendicular to the contact surface and is oriented in the direction from the immersing gelator solution to the polymer solution. The origin of the x -axis is chosen at the contact interface. X(t) denotes the gel layer thickness at immersion time t. In the gelation of chitosan solution, the gelator solution and polymer solution are NaOH aqueous solution and chitosan in acetic acid aqueous solution, respectively.
Gels 04 00059 g001
Figure 2. Time course of gel thickness X observed for one-dimensional growth of chitosan gel induced in 2 wt. % chitosan in 2 wt. % acetic acid aqueous solution in contact with sodium hydroxide aqueous solution with concentrations of 0.1 M (circles), 0.3 M (squares), 0.5 M (upward triangles) and 1 M (downward triangles) [25]. The inset shows the X2 vs. t plot.
Figure 2. Time course of gel thickness X observed for one-dimensional growth of chitosan gel induced in 2 wt. % chitosan in 2 wt. % acetic acid aqueous solution in contact with sodium hydroxide aqueous solution with concentrations of 0.1 M (circles), 0.3 M (squares), 0.5 M (upward triangles) and 1 M (downward triangles) [25]. The inset shows the X2 vs. t plot.
Gels 04 00059 g002
Figure 3. Time course of gel thickness obtained using the function ζ given by Equation (22) [30] for the data in Figure 2. ζ is proportional to t ˜ , as predicted by Equation (21).
Figure 3. Time course of gel thickness obtained using the function ζ given by Equation (22) [30] for the data in Figure 2. ζ is proportional to t ˜ , as predicted by Equation (21).
Gels 04 00059 g003
Figure 4. Illustration of blood coagulation at blood vessel wall.
Figure 4. Illustration of blood coagulation at blood vessel wall.
Gels 04 00059 g004
Figure 5. Gelation dynamics expressed by Equation (54) for horse plasma in contact with calcium chloride with concentrations of 10 mM (△), 30 mM (□), 50 mM (◊) and 100 mM (○).
Figure 5. Gelation dynamics expressed by Equation (54) for horse plasma in contact with calcium chloride with concentrations of 10 mM (△), 30 mM (□), 50 mM (◊) and 100 mM (○).
Gels 04 00059 g005
Figure 6. Time course of gel layer thickness observed on packed swine erythrocytes in contact with plasma with calcium chloride concentrations of 1.63 mM (○), 2.50 mM (△), 4.59 mM (□) and 6.68 mM (×).
Figure 6. Time course of gel layer thickness observed on packed swine erythrocytes in contact with plasma with calcium chloride concentrations of 1.63 mM (○), 2.50 mM (△), 4.59 mM (□) and 6.68 mM (×).
Gels 04 00059 g006
Figure 7. Scheme of coagulation induced by contact of plasma with packed erythrocytes.
Figure 7. Scheme of coagulation induced by contact of plasma with packed erythrocytes.
Gels 04 00059 g007
Table 1. Types of gelation and time development equation.
Table 1. Types of gelation and time development equation.
Case #Geometry *Time Development of Gel ThicknessSystem-Dependent ParametersEquation #Ref.
1One-dimensional x = 2 K t K = k k B T ρ S ρ G T1-1
1Cylindrically symmetrical y ˜ 1 2 ( 1 x ˜ ) 2   ln ( 1 x ˜ ) 1 4 x ˜ 2 + 1 2 x ˜ = K t ˜ K = k k B T ρ S ρ G T1-2[26]
2One-dimensional ζ ˜ ( X ˜ , β K i n ) 0 X ˜ u ˜ 1 β K i n ln ( 1 u ˜ ) d u ˜ = K i n t ˜ K i n = k N a O H k B T α ρ s C i n ,   β = k A c k B T T2-1[25]
2Cylindrically symmetrical z ˜ ( x ˜ ; 2 β K i n ) 0 x ˜ ( 1 u ˜ ) ln 1 1 u ˜ 1 2 β K i n ln ( 1 u ˜ ) d u ˜ = K i n t ˜ K i n = k N a O H k B T α ρ s C i n ,   β = k A c k B T T2-2[28]
3One-dimensional x = 2 K t K = 2 k c t k B T T-3
4One-dimensional x = 2 K t (initial process) K = k k B T ρ S ρ G T1-1
* One dimensional geometry corresponds to Figure 1. Cylindrically symmetrical geometry refers the case such as dialysis of polymer solutions to gelator solutions. x: gel thickness; x ˜ : gel thickness scaled by the radius of cylinder or the effective length of the cell; t: time; t ˜ : time scaled by square of the radius of cylinder or of the effective length of the cell; ρ G : the number of gelator required to produce a unit volume of gel; ρ S : concentration of gelator or acetate ion; C i n : initial concentration of sodium ion; α : a positive numerical factor; k : mobility of gelator; k N a O H : mobility of sodium ion; k A c : mobility of acetate ion; k c t : mobility of catalyst; k B : Boltzmann constant; T: absolute temperature.

Share and Cite

MDPI and ACS Style

Dobashi, T.; Yamamoto, T. Analysis of Heterogeneous Gelation Dynamics and Their Application to Blood Coagulation. Gels 2018, 4, 59. https://doi.org/10.3390/gels4030059

AMA Style

Dobashi T, Yamamoto T. Analysis of Heterogeneous Gelation Dynamics and Their Application to Blood Coagulation. Gels. 2018; 4(3):59. https://doi.org/10.3390/gels4030059

Chicago/Turabian Style

Dobashi, Toshiaki, and Takao Yamamoto. 2018. "Analysis of Heterogeneous Gelation Dynamics and Their Application to Blood Coagulation" Gels 4, no. 3: 59. https://doi.org/10.3390/gels4030059

APA Style

Dobashi, T., & Yamamoto, T. (2018). Analysis of Heterogeneous Gelation Dynamics and Their Application to Blood Coagulation. Gels, 4(3), 59. https://doi.org/10.3390/gels4030059

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