Ion mediated crosslink driven mucous swelling kinetics S. Sircara and A. J. Roberts School of Mathematical Sciences, University of Adelaide, SA 5005, Australia arXiv:1501.05003v1 [q-bio.TO] 20 Jan 2015 a Corresponding Author (email: sarthok.sircar@adelaide.edu.au) Abstract We present an experimentally guided, multi-phasic, multi-species ionic gel model to compare and make qualitative predictions on the rheology of mucus of healthy individuals (Wild Type) versus those infected with Cystic Fibrosis. The mixture theory consists of the mucus (polymer phase) and water (solvent phase) as well as several different ions: H+ , Na+ and Ca2+ . The model is linearized to study the hydration of spherically symmetric mucus gels and calibrated against the experimental data of mucus diffusivities. Near equilibrium, the linearized form of the equation describing the radial size of the gel, reduces to the well-known expression used in the kinetic theory of swelling hydrogels. Numerical studies reveal that the Donnan potential is the dominating mechanism driving the mucus swelling/deswelling transition. However, the altered swelling kinetics of the Cystic Fibrosis infected mucus is not merely governed by the hydroelectric composition of the swelling media, but also due to the altered movement of electrolytes as well as due to the defective properties of the mucin polymer network. Keywords: Donnan potential, Cystic Fibrosis, mucus diffusivity, polyelectrolyte gel 1 Introduction Mucus is a polyelectrolyte biogel that plays a critical role as a protective, exchange and transport medium in the digestive, respiratory and reproductive systems of humans and other vertebrates [1, 2]. Its swelling mechanism are of special interest because of its role in understanding a variety of diseases including cystic fibrosis (CF) [3, 4, 5]. The conformation of the long-chain, negatively charged mucus glycoproteins depends strongly on factors such as pH, ionic strength and ionic bath composition [6]. Mucin is present in secretory vesicles at very high concentrations where they are shielded primarily by a combination of divalent ions (e.g., Ca2+ ) [7]. Experiments show that the mucus gel may swell explosively, up to 600-fold, in times that are the order of a few seconds, a process not observed in hydrogel swelling [8, 9]. Experiments also confirm that this rapid and massive expansion of the mucus gel is driven by an exchange of calcium in the vesicle for a monovalent ion (e.g., Na+ ) in the extracellular environment [9]. This is because, calcium being divalent, must balance two negative charges rather than one. Hence, a divalent Ca2+ ion can act as a ‘cross-linker’ between two polymer strands, allowing much tighter condensation than when the negative charges of the network are shielded by monovalent ions (Fig. 1). Further, experiments demonstrate that the exocytosed mucin is recondensed if the calcium concentration of the ionic medium is increased sufficiently [10]. These observations indicate that the amount and the nature 1 of the salt dissolved in the solvent determines the initial and the equilibrium configuration of these gels. Experimental findings have reported that (unlike hydrogel swelling) a simple osmotic pressure difference of solute particles across the mucus gel cannot explain the massive and the explosive post-exocytotic swelling observations [11]. Mucus hydration results from a balance between osmotic forces, diffusional hindrance (of the ensemble of entangled polymer mesh of the mucus), polyionic charges of the mucin and the other fixed polyions entrapped within the gel’s matrix, the last two factors being the main idea behind the forces via Donnan potential [12]. Figure 1: Schematics showing how the The Donnan potential is not constant, but modulated by ion-displacement changes between calthe concentration of free cations and polycations and the cium ion (Ca2+ ) and a monovalent ion pH of the hydrating fluid [13]. Marriott et al. observed (e.g., Na+ ) leading to changes in the netthat a decreased monovalent ion concentration, or an inwork structure of the mucus matrix and creased Ca2+ concentration, in the airway surface liquid its eventual expansion. was a limiting factor in mucus hydration in CF-infected mucus [14]. Lisle indicated the role of defective processing of mucin polymers, via increased sulfation, in the abnormal mucus hydration in CF [15]. Increased sulfation and sialilation found in CF mucins [16] are important for mucus swelling, since the high affinity of sulfate residues for calcium drastically shifts the normal monovalent/divalent ion-exchange properties of mucins and their swelling characteristics following the exocytotic release [17]. In summary, besides the ionic composition of water on the surface of the airway, the defective processing of mucins and its abnormal ion exchange properties are essential factors to explain the characteristically deficient hydration and rheology of mucus in CF. However, the precise quantitative link between these elements and mucus hydration still remains unknown [18]. The theory to understand the swelling and deswelling of ionic gels has a long history beginning with the classical work of Flory [19, 20, 21] and Katchalsky [22] (see also [23, 24]). An early study to understand the kinetics (and not simply the equilibria) of swelling and deswelling, was done by Tanaka and colleagues [25, 26] who developed a kinetic theory of swelling gels viewing a gel as a linear elastic solid immersed in a viscous fluid. Although they neglected the motion of the fluid solvent, the model reasonably explained the swelling of a gel to its equilibrium volume fraction. Their work gave rise to the concept of gel diffusivity as a way to characterize the kinetics of swelling. The gel diffusivity is defined as D = L2/τ where L is the equilibrium size (length or radius) of a gel and τ is the time constant of exponential swelling toward the equilibrium size. They found that expansion of the gel was governed approximately by a diffusion equation with diffusion coefficient D. However, how the diffusivity of the gel, or more generally the kinetics of swelling (perhaps with large changes in volume fraction) is affected by the movement of the ions and defective properties of the mucus polymer network is poorly understood. Subsequent studies relaxed the assumption of linear elasticity by defining the force on the gel to be the functional derivative of the free energy function for the polymer mesh [27, 28, 29, 30, 31]. However, most of these works neglected the fluid flow that must accompany swelling. Wang et al. [32] added fluid flow by application of two-phase flow theory but considered only small polymer volume fractions and small gradients in the volume fraction. Durning and Morman [27] 2 also used continuity equations to describe the flow of solvent and solution in the gel, but used a diffusion approximation with a constant diffusion coefficient to determine the fluid motion. More recently, Wolgemuth et al. [13] extended these theories to study the swelling of a polyelectrolyte gel. However, these current state of the art models do not couple the binding/unbinding of the ions to the network micro-structure, an idea which is critical in describing the bulk mechanical properties of the polymer [33, 34]. The purpose of this paper is to provide a new comprehensive model detailing the swelling kinetics of mucin-like ionic gels and use it to calibrate the rheological data of spherically symmetric mucin granules, exocytosed from the goblet cells of CF-patients versus those released from the non-CF individuals. The novel feature of this model is the coupling of the dynamical motion of the swelling gel with the ionic binding to the polymer network. We use this model to show how the altered kinetics of the CF-infected mucus depends on the electro-chemical environment of the swelling media as well as the rheological properties of the polymer network. The next section presents the details of this model (§2), including equations of motion (§2.1) and chemistry of binding reactions (§2.2). Sections §2.3 and §2.4 outline the linearized analysis of spherically symmetric swelling gels and methods to estimate model parameters, respectively. The results pertaining to the swelling kinetics and the equilibrium configuration of these ionic gels under different chemical stimuli are presented in §3. We conclude with a brief discussion of the implication of these results on the pathophysiology of the infected mucus. 2 Multi-species, multi-phase mucus-gel model The polymer gel is modeled as a multi-component material, consisting of k different types of particles. Specifically, this material consists of solvent particles, polymers, and several small molecular ion species. The polymer is assumed to be made up of monomers (i.e., charge units), denoted as M, each of which carries a single negative charge. The positively charged ions in the solvent are Hydrogen (H+ ), Sodium (Na+ ) and Calcium (Ca2+ ). The negatively charged ions could include Hydronium (OH− ) and Chloride (Cl− ). Because the negatively charged ions are assumed to be not involved in any binding reactions with the gel, acting only as counterions to positive charges, we identify these ions by the name chloride. The binding reactions of the positively charged ions with the monomers are k k k k k−h k−n k−c k−x n c x h M− + H+ −− −− HM, M− + Na+ −− −− NaM, M− + Ca2+ −− −− MCa+ , M− + MCa+ −− −− M2 Ca, (1) where kC and k−C , for C = h, n, c, x, are the binding and the unbinding rates, respectively. We assume that all the binding sites/charge sites are identical although the binding affinities for the different ions are different. The species M2 Ca are cross-linked monomer pairs, and the species M− , MCa+ , NaM and HM are different monomer species, all of which move with the same polymer velocity. The ion species are freely diffusible, but because they are ions, their movement is restricted by the requirement to maintain electroneutrality. Finally, because a small amount of water dissociates into hydrogen and hydronium, we are guaranteed that there are always some positive and negative ions in the solvent. 3 2.1 Equations of motion and interface conditions Suppose we have some volume V of a mixture comprised of k types of particles (including polymer, solvent and ion species) each with particle density (number of particles per unit volume), nj (x, t) and particle volume, νj , moving with velocity, vj (x, t), j = 1, . . . , k. The subscripts (1, 2) denote the polymer and the solvent phase, respectively. We denote these phases with subcripts (p, s) respectively. The volume fractions for polymer and solvent are θp = νp np and θs = νs ns . Conservation of polymer implies ∂θp + ∇ · (vp θp ) = 0. (2) ∂t We assume that the other molecular species do not contribute significantly to the volume. Therefore, the volume fractions, θp + θs = 1. It follows from (2) and a similar conservation argument for solvent that ∇ · (θs vs + θp vp ) = 0. (3) The motion of the polymer and solvent phase of this multi-component mixture is governed by the Stokes equation for Newtonian fluid, which are ∇ · (θp σp (vp )) − ξ and ∇ · (θs σs (vs )) − ξ θs θp φp (vp − vs ) − ∇µp = 0, νs νp θs θs φp (vs − vp ) − νs νs θs φˆj ∇µj − ∇µs = 0, νs j≥3 (4) (5) η where σj (v) = 2j (∇v+∇vT )+λj I∇·v, are the viscous stresses, ηj > 0 and λj are the viscosities, ξ is the drag coefficients. φp = np/np +ns , φs = ns/np +ns , φj = nj/ i=1,2 ni for j ≥ 3 (assuming that the ions are dissolved in the solvent), are the polymer, solvent and ion species per total solvent particle fractions, respectively. The first and second terms in P D E (4, 5) represents the viscous forces and drag forces due to the friction between the two phases, respectively. The third term in P D E (5) represents the force due to the dissolved counter-ions in the solvent (osmotic effect), while the last terms in P D E (4, 5) are the chemical forces due to the respective phases. The chemical potentials are given by µp = kB T Mp + zm Φe + νp P, µs = kB T Ms + νs P, (6) where Mp and Ms represent the entropic contribution to the respective chemical potentials (described later in this section), kB is the Boltzmann constant, T is temperature, P is pressure, zm is the average charge per monomer (defined later in §2.2), Φe is the electric potential arising from the unbalanced charges. For the ion species, the particle volume, νj , is effectively zero so that µj = kB T (ln φj + 1 − 2σI ) + zj Φe , j ≥ 3, (7) where zj is the charge on the jth ionic species. In calculating particle fractions, φj , we assume that the particle density of the ions is insignificantly small compared to the particle density of polymer and solvent, j≥3 nj ns + np . The ion species satisfy the force balance ξj nj (vs − vj ) − nj ∇µj = 0, 4 j ≥ 3. (8) Finally, since there is a free moving-edge to the gel, on one side of which (inside the gel) θp− = θp , and on the other side of which (outside the gel) θp+ = 0, θs+ = 1, implying that there is no polymer outside the gel. The superscripts (+, −) denote regions inside and outside the gel, respectively. The interface conditions are σp (vp− )n = kB T νm and Mp− + zm Ψe + νp kB T νs σs (vs+ ) − σs (vs− ) n = P kB T n, Ms+ − Ms− − σI+ + σI− − νs (9) P kB T n, (10) where 1 ln φp + N 1 T0 χ 2 − 1 φs + φ + µp0 , N T 2 s 1 T0 χ 2 Ms = ln φs + 1 − φ + µs0 . φp − σI + N T 2 p Mp = (11) N is the number of monomers per polymer-chain and T0 is a reference temperature. The hydrostatic pressure is P − = P , P + = 0 and the normalized electrostatic potential (Ψe = Φe/kB T ) is, Ψ− e = Ψe , Ψ+ = 0. The normal to the free surface is denoted by n. The term σ in the solvent chemical I e nj/ns ) and represents osmotic potential (Eq. (11)) is the total ion particle fraction (σI = j≥3 pressure as characterized by van’t Hoff’s law. The quantities χ, µp0 , µs0 are the Flory interaction parameter, standard free energies for pure polymer and pure solvent respectively. These are related to the cross-link fraction, α (the fraction of monomers bound with calcium), χ = z( 1 + z + 2 z µs0 = − 2 . 2 µp0 = − 1 2) 4 −2 1− 1− 1 N 1 N + 1 − 1 α, 1 3 α, 2 (12) (i = 1, . . . , 4) are the nearest neighbor interaction energies parameters of the various monomermonomer and monomer-solvent pairs and z is the number of interaction sites on a lattice (i.e., coordination number) [35]. Eliminating P from Eqns. (9, 10) we find a single interface condition i σp (vp− ) − σs (vs− ) + σs (vs+ ) n = Σnet n, (13) where Σnet is the net swelling pressure at the interface, Mp− Ms− T0 µ0s zm Σnet σ− σ+ = − + + Ψe + I − I . kB T νm νs T νs νm νs νs (14) In the above equation, the fourth term represents swelling pressure coming from any electric charge on the monomers, referred as the Donnan pressure (zp Ψe is the corresponding Donnan potential) whereas the last two terms represents the osmotic swelling pressure coming from the difference between the concentrations of ions dissolved in the gel and those dissolved in the bath. Sircar et al. 5 [35] derived these equations of motion and the interface conditions using the standard variational arguments to minimize the rate of work dissipated within the polymer and the solvent. Unlike previous theories, this model accurately captures how the binding/unbinding of the dissolved ions influences the motion of the swelling gel. The chemistry of the dissolved ions and the charged polymer species is described next. 2.2 Ionization chemistry This section formulates a mathematical model to represent how the chemical species move and react. Let the concentrations per total volume of the polymer species be denoted by x = [M2 Ca], m = [M− ], v = [NaM], w = [MCa+ ] and y = [HM], with the total monomer concentration mT = m + 2x + w + v + y. (15) The concentrations per solvent volume of the ion species are denoted as c = [Ca2+ ], n = [Na+ ], h = [H+ ], and cl = [Cl− ]. With concentrations expressed in units of moles per liter, the relationship between ion particle fractions φj and concentrations cj is φj = νs NA cj , where NA is Avagadro’s number. To describe the chemical reactions, we use the law of mass action. Since all the monomer species are advected with the polymer velocity vp , the monomer species evolve according to ∂j + ∇ · (vp j) = Rj+ − Rj− , ∂t j = x, v, w, y (16) where Rj+ , Rj− are the forward and backward rates for the monomer binding reactions in Eqn. (1), respectively. The monomer concentration, m, is obtained from Eqn. (15) and mT = θp/νp NA . Under the assumption of fast chemistry, we set the right hand side of the P D E (16) to zero, which reduces into the following set of equations for each of the monomer species w= θs mc, Kc φ2s v= θs mn, Kn φ2s y= θs mh, Kh φ2s x= θs m2 c, 2 4 4Kc φs (17) where Kc = k−c/kc , Kn = k−n/kn , Kh = k−h/kh , Kx = k−x/kx = 4Kc . The above expressions (17), assume that the unbinding (dissociation) reactions are ionization reactions that require two “units” of solvent. We take the unbinding reaction rates to be k−C φ2s , for C = c, x, n, h and because calcium is a divalent ion, 2kx = kc and k−x = 2k−c . Similarly, under the assumption of fast diffusion and chemistry, the law of mass action for the ion species reduces to C = Cb e−zC Ψe (18) with symbol C = c, h, n, cl , zn = zh = 1, zc = 2 and zcl = −1. The subscript ‘b’ denotes the corresponding bath concentrations [35]. The electrostatic potential, Ψe , is determined by the electroneutrality constraint inside the gel, namely, (2c + n + h − cl )θs + zm mT = 0, (19) where zm is the average residual charge of the unbound monomers which depends on the amount of binding with ions, zm mT = w − m. (20) 6 Since both the electrostatic potential and polymer particle fraction are assumed to be zero outside the gel, electroneutrality in the bath requires that 2cb + nb + hb − clb = 0. (21) Finally, the crosslink fraction (or the fraction of monomers bound with calcium) is α = mxT , where Eqns. (15, 17) define the concentrations, mT , x, respectively. In summary, the dynamical motion of a freely swelling mucus-gel is modeled by the system of equations including the mass conservation P D E (2), total volume conservation P D E (3) together with force balance P D E S (4-5) and interface condition Eqn. (13), subject to the constraints Eqn. (15) (monomer conservation), Eqn. (18) (ion motion) and Eqn. (19) (electroneutrality). 2.3 Linearized analysis of spherically symmetric swelling gels We now consider the swelling kinetics of a radially symmetric mucus gel using a linearized analysis of the governing equations, described in the previous two sections. The volume fraction θp is nonzero on the domain 0 ≤ r < R(t), with R (maximum radius of the sphere) a function of time. Since only the finite solutions of the governing equations are of interest (in particular, finite solutions of the force balance, P D E S (4, 5)) inside this domain, we assume that vp = vs = 0 at r = 0. The volume conservation, P D E (3), implies the constraint θp vp +θs vs = 0 throughout the domain. Using this constraint, we reduce the two force balance equations into a single equation (by multiplying P D E (4) by volume fraction θs and P D E (5) by volume fraction θp and subtracting) ∂ ∂vp λp ∂ 2 + 2 r vp + θp ∇ · θs ηs ∂r r ∂r ∂r φp φs + − θp θs ∇ (Σnet ) = 0, νs νp θs ∇ · θp ηp + ξvp θp θp vp θs + λs ∂ r2 ∂r r2 θp vp θs (22) whereas the interface condition (Eqn. (13)) reduces to ∂ ∂vp + ηs ηp ∂r ∂r θp vp θs λs ∂ λp ∂ + 2 r 2 vp + 2 r ∂r r ∂r r2 θp vp θs = Σnet at r = R(t). (23) Note that ∇( νµmp − µνss ) = ∇(Σnet ), where Σnet is the interface swelling pressure defined in Eqn. (14). − + The velocities inside the gel are vp/s = vp/s and outside the gel are vp/s = 0. Next, we assume that ∗ the linearized velocity vp is small (the equilibrium velocity vp = 0) and that θp = θp∗ + δθp , where θp∗ is the equilibrium polymer volume fraction. Sircar et al. [35] gives details on the equilibrium solution. Since this is a moving boundary problem, it is convenient to map the domain 0 ≤ r < R(t) onto the fixed domain 0 ≤ y < 1 by making the change of variables r = R(τ )y and t = τ . Further, we seek space-time variable separated solutions of the form vp = a(τ )f1 (y) and δθp = a(τ )f2 (y). Under these assumptions the reduced force balance, P D E (22) is f1 + 2(ηe + 2λe )R 2λe ξR f1 + + 2 (ηe + λe )y (ηe + λe )y η e + λe φ∗s φ∗p + νp νs f1 = −θs∗ RΣθnet (θp∗ )f2 , (24) and the interface condition, Eqn. (23), is (ηe + λe )f1 + 2λe f1 = θs∗ Σθnet (θp∗ )f2 y 7 at y = 1. (25) The supercripts ( ) and (θ ) denote the derivative of the functions with respect to the variables y and θ, respectively. The net shear and bulk viscosities are given by ηe = θs∗ ηp + θp∗ ηs and λe = θs∗ λp + θp∗ λs . φ∗p , φ∗s are the equilibrium particle fractions for the polymer and solvent, respectively. For a finite solution of Eqn. (24), it is assumed that f1 (0) = 0. Further, we assume that f2 = f1 . These assumptions reduce Eqn. (24) into the homogeneous, spherical Bessel differential equation whose solutions have the form f1 = y γ Jn (βy), where γ = 1 2 − (ηe +2λe )R , ηe +λe β = ∗ ξR ( φs ηe +λe νp + φ∗p ) νs + θs∗ RΣθnet (θp∗ ) e |. The Bessel functions, Jn , are of the first kind of order n. Since there are and n = |γ 2 − ηe2λ +λe several solutions satisfying the condition f1 (0) = 0, we choose the solution with the lowest order n = nmin , which determines the equilibrium radius of the swelling gel, Rf = R(τ → ∞), Rf = 1 − 2 n2min + 2λe ηe + λe ηe + λe ηe + 2λe . (26) The variable-separable, linearized form of solution (linearized near equilibrium) for mass conservation, Eqn. (2), gives 1 ∂a 1 θp∗ 1 =− (y 2 f1 ) = − , (27) 2 a ∂τ f2 Rf y τch where τch is a constant (to be determined later). The first and the last part of Eqn. (27) implies that a(τ ) = e−τ /τch . Because the polymer is conserved (neither created nor destroyed) the velocity of the moving boundary must be the same as the gel velocity at the boundary, vp (y = 1), ∂R = vp (y = 1) = e−τ /τch f1 (1). ∂τ (28) The solution to the above equation is numerically computed via Matlab ODE solver ode15, with trivial initial conditions. In particular, we note that if R ≈ Rf , then f1 (1) ≈ constant and in this case Eqn. (28) can be solved exactly, i.e., R(τ ) = Rf (1 − e−τ /τch ). This expression is identical to the radial expansion of the linearized form given by Tanaka et. al [25] in their kinetic theory of swelling hydrogels. Finally, the time constant of gel-swelling towards the equilibrium size (Eqn. (27)), is f2 R f y 2 τch = ∗ 2 , (29) θp (y f1 ) y=1 which defines the diffusivity for spherical gels, D = Rf2/τch . 2.4 Parameter Estimation Two sets of data are used to calibrate the model for the kinetics of spherically symmetric swelling mucus gels. The first experiment, conducted by Verdugo et al., measures the diffusivity data (Rf2 versus τch values) of the swelling kinetics of exocytotic mucin granules from four (human) CFpatients and three healthy individuals at pH = 7.2, T = 37◦ C and [Ca2+ ]b = 1mM as well as [Ca2+ ]b = 2.5mM [8]. In the second experiment, performed by Kuver et al., kinetic data was collected from cultured gallbladder epithelial cells from wild-type and CF-infected mice at pH = 7.0, T = 37◦ C, [Ca2+ ]b = 4mM and [Na+ ]b = 140mM [36]. The microscopic composition of the 8 mucus in both of these experiments were reported identical. Fig. 2a and Fig. 2b presents the results from these two experiments, respectively. Table 1 lists the values of the parameters used in our numerical calculations. The constants in the model are the monomer particle volume, νp , the solvent particle volume, νs , the coordination number of the polymer lattice, z, and the nearest neighbor interaction energies, i (Eqn. (12)), shear and bulk viscosity coefficients of solvent (i.e., water), ηs , λs , respectively, at reference temperature. The undetermined constants are the binding affinities of the various cations with the gel, Kh , Kc , Kn (introduced in Eqn. (17)) and the polymer viscosity coefficients, µp , λp , and the drag coefficient, ξ (Eqn. (4)). Table 1: Constant parameters common to all the numerical simulations. The reference temperature for the viscosity coefficients and solubility parameters is fixed at T0 = 25◦ C. Constants Value Units Source Repeat unit per chain (N ) 266 – [7] −20 3 Molecular volume of mucus (νp ) 5 × 10 m [7] Molecular volume of water (νs ) 2 × 10−23 m3 [7] −4 Shear viscosity of water (ηs ) 8.88 × 10 Pa s [37] Bulk viscosity of water (λs ) 2.47 × 10−3 Pa s [37] 1/2 Hildebrand solubility (δp ) 1.0928 (α = 0), 1.3258 (α = 1) MPa [38] 1/2 Hildebrand solubility for water (δs ) 48.07 MPa [39] The radially symmetrically swelling mucus gel has a 3-D configuration, which suggests that we choose the coordination number, z = 6, mimicking the 3-D structure of the polymer lattice [19]. The standard free energies, kB T0 µ0p and kB T0 µ0s and the energy interaction parameters, i (Eqn. (12)) are found from the Hildebrand solubility data, δi [39]. The values for the solubility data for materials mimicking mucus glycoproteins are given in a study by Mimura [38]. The fully un-crosslinked (no Ca2+ binding) and fully crosslinked (Ca2+ bound) states are denoted by α = 0 and α = 1, respectively. The standard free energy is the energy of all the interactions between the molecule and its neighbors in a pure state that have to be disrupted to remove the molecule from the pure state. The relation between the standard free energies and the solubility parameters (Table 1) are −kB T0 µ0s = νs δs2 −kB T0 µ0p = νm δp2 , (30) where νs = 2 × 10−23 cm3 is the volume of one molecule of water at reference temperature, T0 = 298K. The negative sign in Eqn. (30) indicates that kB T0 µ0p , kB T0 µ0s < 0, since they are the interaction energies. Using the relations in Eqn. (30), these values are fixed at 1 = 4.84, 2 = 3.74, 3 = −13.70, 4 = 0. The reference temperature of the experiments was fixed at T0 = 25◦ C. The volume of a mucin oligomer/multimer chain is calculated by mathematically modeling a chain comprising N freely jointed cylindrical segments of length l and width d, where N = L/l is the number of Kuhn segments (or effective rigid cylindrical segments which determines the different conformations a chain can have), and L is the end-to-end length of the chain. For a gel forming mucin (e.g., human M U C 5 A C , used in the experimental data to calibrate our model), l ≈ 0.03 µm, L ≈ 8 µm and N = 8/0.03 ≈ 266 [6, 40]. The radius of gyration, Rg , (defined as the average 9 distance from all Kuhn segments to the center of mass of the chain) and the pervade volume, V , (i.e., the approximate spherical volume of a sphere with radius Rg ) is [23] Rg ≈ l0.8 × d0.2 × N 0.588 , 2.45 V ≈4× π × Rg3 . 3 (31) Substituting values, the volume V (= νp ) ≈ 0.05 µm3 . The undetermined parameters, namely the binding affinities Kh , Kn , Kc and the rheological coefficients ηp , λp , ξ are computed by minimizing a nonlinear least-square difference function between the experimental values of diffusivity data (Rf2 versus τch ) and the corresponding model output (Eqns. (26, 29)), implemented via the M AT L A B non-linear least-square minimization function lsqnonlin. These values are found as log10 (Kn ) = −2.27, log10 (Kh ) = −3.65, log10 (Kc ) = −3.12, ηp = 0.11, λp = 0.31, ξ = 1.97 (wild-type mucus) and log10 (Kn ) = −2.55, log10 (Kh ) = −3.98, log10 (Kc ) = −7.12, ηp = 1.02, λp = 2.87, ξ = 0.21 (CF-infected mucus). The closeness of fit between the model (highlighted by the solid lines) and the experimental data points is shown in Fig 2. The error bars represent the maximum and minimum deviation from the sample points and set at 5% margin of error. Notice the linear relationship between Rf2 and τch , predicted by the experiments and accurately captured by the model. 3 Results and discussion The main idea behind this work is to provide an objective comparison of the swelling properties of normal versus unhealthy mucus, immersed in an extracellular medium with chemically controlled composition. Well established results exist in literature which detail the relationship between the swelling rate, final size of the gel and the swelling time [25, 26]. However, these results do not explicitly outline how these and other physiochemical elements influencing the mucus hydration and rheology, depend on the nature of the swelling media. Hence, in §3.1, we explore the relationship between the radial size of the swelling mucus gel and the rheology of the mucus polymer, name the drag and the viscosity coefficients. The effect of the calcium and the sodium ions in the solvent on the equilibrium size of the mucus blob are detailed in §3.2 and §3.3, respectively. 3.1 Swelling kinetics Numerical simulations were performed to outline the differences between the swelling kinetics for WT and CF-infected mucus, by altering the electro-chemical composition of the swelling media. Fig. 3 presents the radius of the swelling gel, R(τ ), versus time, τ , for WT mucus with calcium bath concentrations, Cb = 1 mM and Cb = 2.5 mM (the ‘dash-dot’ and ‘solid’ curve, respectively), as well as for CF-infected mucus (highlighted by the ‘dotted’ and ‘dashed’ curve, respectively). The sodium ion concentration and the pH in the bath are fixed at Nb = 0 and pH = 7.2. These concentrations correspond to the in vitro conditions of Verdugo’s experiments [8]. In particular, note that the expansion of the CF-infected mucus (‘dashed’ and ‘dotted’ curves, Fig. 3) is appreciably slower than WT mucus (‘dash-dot’ and ‘solid’ curves, Fig. 3). This difference in swelling profiles is a consequence of small drag coefficient of the CF-infected mucus gels, relative ξ to the viscosity coefficients (i.e., the ratio δ = θ∗ (ηp +λp )+θ is small for CF-infected mucus). ∗ s p (ηs +λs ) Using the rheological parameters, estimated in §2.4 (i.e., ηp , λp , ξ), the solvent parameters listed in 10 (a) pH = 7.2, [Ca2+ ]b = 1 mM, [Na+ ]b = 0 (b) pH = 7.0, [Ca2+ ]b = 4 mM, [Na+ ]b = 140 mM Figure 2: Experimental data of the final size squared (Rf2 ) versus characteristic swelling time (τch ) for the two in vitro electro-chemical conditions mentioned in §2.4. For WT mucus the data points are represented by ( ) and for CF-infected mucus these points are shown by ( ). The temperature of the experiments was set at 37◦ C. The model output (i.e., solid lines) closely fits the experimental data points. Numerical simulations are done for the data represented by the points (i) and (ii). These are shown in Fig. 3. 11 Figure 3: Swelling radius, R(τ ), of a spherically symmetric mucin gel versus time, τ , with kinetics governed by Eqn. (28), for WT-mucus at calcium concentration Cb = 1.0 mM (dash-dot curve), CF-mucus at Cb = 1.0 mM (dotted curve), WT-mucus at Cb = 2.5 mM (solid curve), CF-mucus at Cb = 2.5 mM (dashed curve). The (Rf , τch ) values for the dash-dot and the dotted curve correspond to the data represented by points (i) and (ii) in Fig. 2a, respectively. Table 1 (i.e., ηs , λs ) and the equilibrium polymer volume fraction, θp∗ (i.e., by solving Eqns. (13, 19), values shown in Fig. 4a); we find that δ = 82.8655, 78.7264 for WT mucus for concentrations Cb = 1 mM and Cb = 2.5 mM, respectively while the corresponding values are δ = 0.0885, 0.0914 for CF-infected mucus, respectively. Gels with higher viscosities expand at a slower rate and this explains a relatively slow or insufficient hydration due to the defective rheology of the CF-infected mucus. Previous numerical studies by Sircar at al. corroborate these results [41]. The corresponding diffusivity data for the curves in Fig. 3 are given in Table 2. The diffusivity decreases by 33% for normal mucus and 88% for CF-infected mucus, when the Ca2+ concentration in the swelling medium is changed from 1 to 2.5 mM, respectively. Further, there is a 76% reduction in the diffusivity of CF-infected mucus versus the WT mucus in a solvent containing 1 mM Ca2+ , in concurrence with Verdugo’s swelling experiments [8]. A lower diffusivity of the defective mucus is a consequence of the slow rate of expansion of these gels (or a bigger time constant, τch ), which again lead to the conclusion that the defective rheology of mucus plays a crucial role in the eventual mucus hydration. Table 2: Diffusivity data for WT and CF-infected mucus at calcium bath concentrations Cb = 1 mM and Cb = 2.5 mM. Rf , τch and D are measured in µm, sec and cm2 /s, respectively. The sodium concentration and the pH in the bath are fixed at Nb = 0 and pH = 7.2. 1 mM Ca2+ 2.5 mM Ca2+ WT Rf = 2.41, τch = 0.2, D = 2.90 × 10−7 Rf = 2.37, τch = 0.29, D = 1.93 × 10−7 CF Rf = 2.37, τch = 0.79, D = 7.09 × 10−8 Rf = 2.37, τch = 6.58, D = 0.85 × 10−8 12 3.2 Equilibrium configuration versus bath [Ca2+ ] In the next two sections, we explore the role of the electrolytic composition of the swelling media on the equilibrium configuration of the healthy versus the diseased mucus. Fig. 4a,b,c highlight the equilibrium volume fraction, θp∗ (which determines the swelled/deswelled state of the gel), the Donnan potential and the crosslink fraction, at pH= 7.2 and pH= 5.0, versus the bath concentration of calcium, respectively. The emergence of three different swelling regions are observed in Fig. 4. At low bath concentrations of calcium (Cb < 10−14 M), the gels are immersed in a solvent which lacks cations to bind with the charged monomers. This leads to a highly ionized state (i.e., a high average negative charge per monomer) and consequently a higher Donnan potential (Fig. 4b). Sircar et al. has shown earlier that Donnan potential is the dominant mechanism driving the swelling of ionic gels [35]. Solvents with lower pH (swelling profiles represented by the ‘dash-dot’ and ‘dotted’ curves, Fig. 4a) can furnish more H+ ions, leading to a lower Donnan potential which translates to a relatively de-swelled state (i.e., comparing the θp∗ values at pH= 5 versus those at pH= 7.2). The hydration properties of the diseased versus healthy mucus are markedly different for small to intermediate calcium concentrations, 10−14 M<Cb < 10 M (Fig. 4). In particular, the diseased mucus (swelling profiles shown by the ‘dashed’ and ‘dotted’ curve) shows massive deswelling. The explanation for these swelling features in this region is that there is a complex interplay of ionization via Donnan potential (which favors swelling) and the energy gain due to increased crosslinking (which promotes deswelling) [35]. Defective mucus gels have higher calcium binding affinity (i.e., log10 Kc = −7.12 for CF-infected mucus versus log10 Kc = −3.12 for WT mucus, refer §2.4 where these binding affinity values are mentioned), This leads to a relatively de-ionized state in the defective gels (or a lower Donnan potential, Fig. 4b). Consequently, a de-swelled state (mediated by the de-ionization due to calcium crosslinks) dominates. At sufficiently high calcium ion concentrations (Cb > 10 M) the average charge on the monomer is positive rather than negative. The positive charge on the monomer is because in the divalent ion case, the binding of a monomer with an ion converts it from a negatively charged ion, M− , into a positively charged ion, MCa+ (Eqn. (1)). This has little effect on the overall swelling pressure since at high ion concentrations the charge on the gel is small relative to the overall number of available ions, and the Donnan potential, vanishes. As a consequence, at sufficiently high Cb the gel (both normal as well as the diseased type) behaves as if it is uncharged. In summary, an increased affinity of the CF-infected mucus to bind with calcium (in physiologically significant concentration range) leads to a highly de-swelled state of the gel. This defective binding property of the diseased mucus partly explains how the altered movement of the electrolytes leads to an insufficient hydration of the gel. 3.3 Equilibrium configuration versus bath [Na+ ] This section reports the salient features of the equilibrium state of the mucus gel immersed in a bath containing monovalent cations (i.e., Na+ , Fig. 5). Again, three different swelling transition regimes are seen from these figures: a de-swelled state (i.e., a relatively high polymer volume fraction, θp∗ ) is favored at low (Nb < 10−7 M) and high (Nb > 1 M) bath concentrations. A higher polymer volume fraction is because of a low Donnan potential in these concentration ranges, which is due to either the furnished H+ ions (at low sodium concentrations) or excessive Na+ ions (otherwise). 13 (a) Equilibrium polymer volume fraction, θp∗ (b) Donnan swelling pressure, zm Ψe 14 In small to intermediate sodium concentrations (10−7 M<Nb < 10−1 M, i.e., approximately in the region Nb ≈Hb ), there is an additional feature with increasing sodium concentration, namely, swelling either gradually (pH = 5.0 curves, Fig. 5a) or via phase-transition (pH = 7.2 curves, Fig. 5b) followed by de-swelling. The explanation for this swelling feature is a complicated interplay between ionization via hydrogen unbinding (which promotes swelling) and a de-ionization via sodium binding (which promotes de-swelling). However, the volume transition curves are qualitatively different for the WT and the CF-infected mucus, e.g., notice the hysteretic swelling transition for the WT-mucus (solid curve) versus a doublehysteretic curve for the defective mucus (dashed curve) at pH = 7.2. These differences, again, arise due to different binding affinities (log10 Kn = −2.27, log10 Kh = −3.65 for CF-infected mucus versus log10 Kn = −2.55, log10 Kh = −3.98 for WT mucus, values listed in §2.4). Thus, the differences in the ion-binding property of the mucus influences the not only the swelling-deswelling volume transitions but also the equilibrium state of the gel. 4 Conclusions This paper develops a new, comprehensive, multi-phase, multi-species model to quantify the swelling / deswelling mechanism for mucin gels. This model explains how the final configuration of these gels depends on complex interactions between competing effects that alters the gel ionization, that is the Donnan potential, changes in the bath concentration of ions and their corresponding binding affinity with mucus. Near equilibrium, the radial size of the swelling mucus gel reduces to the well known expression for hydrogel swelling [25]. The diffusivity, D (which is similar to the one defined by Tanaka’s hydrogel theory) accurately characterizes the swelling properties of the mucin network, including the effect of calcium binding on the equilibrium configuration of the gel. The reasons behind the experimental outcomes of deficient mucus hydration for CF-infected patients are examined. In particular, the defective swelling properties of the diseased mucus is the combined effect of an abnormal mucus rheology (i.e., an abnormally low drag relative to the polymer viscosity), an altered electrolytic composition of the extracellular medium (i.e., an increased calcium ion concentration in the airway surface liquid surrounding the mucus) as well as an altered movement of ions in the medium (i.e., a higher affinity of the CF-infected mucus to bind with calcium ion). All these factors favor a condensed state of the gel, resulting in a mucus gel with defective mucociliary transport contributing to the chronic airway infection found in CF patients. References [1] P. Verdugo, ‘Hydration kinetics of exocytosed mucins in cultured secretory cells of the rabbit trachea: a new model’, Mucus and mucosa, Ciba Foundation symposium 109, 212225 (1984). [2] P. Verdugo, ‘Polymer gel phase transition in condensation-decondensation of secretory products’, Adv. Polm. Sci., 110, 145-156 (1993). [3] J. Barasch, B. Kiss, A. Prince, L. Saiman, D. Gruenert and Q. A. Awqati, ‘Defective acidification of intracellular organelles in cystic fibrosis’, Nature, 352, 70-73 (1991). 15 (a) Equilibrium polymer volume fraction, θp∗ (b) Donnan swelling pressure, zm Ψe Figure 5: Variables in equilibrium versus sodium bath concentration, Nb (in mol/lt) for WT-mucus at pH = 7.2 (solid curve), CF-infected mucus at pH = 7.2 (dashed curve), WT-mucus at pH = 5.0 (dash-dot curve), CF-infected mucus at pH = 5.0 (dotted curve). The bath concentration of calcium and the temperature are fixed at Cb = 0 and T = 37◦ C, respectively. 16 [4] P. W. Cheng, T. F. Boat, K. Cranfill, J. R. Yankaskas and R. C. Boucher, ‘Increased Sulfonation of Glycoconjugates by Cultured Nasal Epithelial Cells from Patients with Cystic Fibrosis’, J. Clin. Invest., 84, 68-72 (1989). [5] R. Kuver, J. H. Klinkspoor, W. R. A. Osborne and S. P. Lee, ‘Mucous granule exocytosis and CFTR expression in gallbladder epithelium’, Glycobiology, 10(2), 149-157 (2000). [6] R. Bansil and B. S. Turner, ‘Mucin structure, aggregation, physiological functions and biomedical applications’, Curr. Opin. Colloid and Interface Sci., 11, 164-170 (2006). [7] J. P. Villar, ‘Mucin Granule Intraluminal Organization’, Am. J. Respir. Cell Mol. Biol., 36, 183-190 (2007). [8] P. Verdugo, Cilia, Mucus and Mucociliary Interactions, Chap 19, 167-189, Marcel Decker, New York (1998). [9] P. Verdugo, M. Aitken, L. Langley and M.J. Villalon, ‘Molecular mechanism of product storage and release in mucin secretion. II. The role of extracelluar Ca++ ’, Biorheology, 24, 625-633 (1987). [10] S. Baconnais, F. Delavoie, J.M. Zahm, M. Milliot, C. Terryn, N. Castillon, V. Banchet, J. Micheal, O. Danos, M. Merten, T. Chinet, K. Zierold, N. Bonnet, E. Puchelle and G. Balossier, ‘Abnormal ion content, hydration and granule expansion of the secretory granules from cystic fibrosis airway glandular cells’, Exp. Cell Res., 309, 296-304 (2005). [11] P. Y. Tam and P. Verdugo, ‘Control of mucus hydration as a Donnan equilibrium process’, Nature, 292, 340-342 (1981). [12] A. Katchalsky, S. Lifson and H. J. Eisenberg, ‘Equation of swelling for polyelectrolyte gels’, J. Polym. Sci., 7, 571-574 (1951). [13] C. W. Wolgemuth, A. Mogilner and G. Oster, ‘The hydration dynamics of polyelectrolyte gels with applications to cell motility and drug delivery’, Euro Biophys. J., 33, 146-158 (2004). [14] R. S. Crowther and C. Marriott, ‘Counter-ion binding to mucus glycoproteins’, J. Pharm. Pharmacol., 36, 21-26 (1984). [15] R. C. Lisle, ‘Increased expression of sulfated gp300 and acinar tissue pathology in pancreas of CFTR (-/-) mice’, Am. J. Physiol., 268, G717-G723 (1995). [16] K. V. Chase, D. S. Leahy, R. Martin, M. Carubelli, M. Flux and G. P. Sachdev, ‘Respiratory mucus secretions in patients with cystic fibrosis: relationship between levels of highly sulfated mucin components and severity of disease’, Clin. Chim. Acta, 132, 143-155 (1983). [17] J. Barash, B. Kiss, A. Prince, L. Saiman, D. Gruenert and Q. Al-Awqati, ‘Defective acidification of intracellular organelles in cystic fibrosis’, Nature, 352, 70-73 (1991). [18] T. Nguyen, W. C. Chin and P. Verdugo, ‘Role of Ca2+ /K+ ion exchange in intracellular storage and release of Ca2+ ’, Nature, 395, (1998). 17 [19] P. J. Flory, Principles of Polymer Chemistry, Cornell University Press, NY (1953). [20] Flory P.J. and Rehner J., ‘Statistical mechanics of crosslinked polymer networks I. rubberlike elasticity.’, J. Chem. Phys. 11, 512 (1943). [21] Flory P.J. and Rehner J., ‘Statistical mechanics of crosslinked polymer networks II. swelling.’, J. Chem. Phys. 11, 521 (1943). [22] Katchalsky A. and Michaeli I., ‘Polyelectrolyte gels in salt solutions..’, J. Polym. Sci. 15, 69-86 (1955). [23] Doi M. and Edwards S. F., ‘Theory of polymer dynamics’, Clarendon Press, Oxford, England (1986). [24] Flory P.J., ‘Statistical thermodynamics of random networks, Proc. R. Soc. Lond. A. 351, 351-380 (1976). [25] T. Tanaka and D. Fillmore, ‘Kinetics of swelling gels’, J. Chem. Phys., 70(03), 1214 (1979). [26] T. Tanaka, D. Fillmore, S. T. Sun, I. Nishio, G. Swislow and A. Shah, ‘Phase transition in Ionic Gels’, Phys. Rev. Lett., 45(20), 1636 (1980). [27] Durning C.J. and Morman K.N., ‘Nonlinear swelling of polymer gels.’, J. Chem. Phys. 98, 4275-4293 (1993). [28] Maskawa J., Takeuchi T., Maki K., Tsuijii K. and Tanaka T., ‘Theory and numerical calculation of pattern formation in shrinking gels’, J. Chem. Phys. 110, 10993-10999 (1999). [29] Onuki A., ‘Theory of pattern formation in gels: surface folding in highly compressible elastic bodies’, Phys. Rev. A, 39, 5932-5948 (1989). [30] Sekimoto K., Suematsu N. and Kawasaki K., ‘Spongelike domain structure in a twodimensional model gel undergoing volume-phase transition’, Phys. Rev. A, 39, 4912-4914 (1989). [31] Yamaue T., Taniguchi T. and Doi M., ‘Shrinking process of gels by stress-diffusion coupled dynamics’, Prog. Theor. Phys. Supp., 138, 416-417 (2000). [32] Wang C., Li Y. and Hu Z., ‘Swelling kinetics of polymer gels’, Macromolecules, 30, 47274732 (1997). [33] Calderer M.C., Chen H., Micek C. and Mori Y., A Dynamic Model of Polyelectrolyte Gels, SIAM J. Appl. Math. 73(1), 104-133 (2013). [34] English A.E., Mafe S., Manzanares J., Yu X., Grosberg A.Y. and Tanaka T.,‘Equilibrium swelling properties of polyampholytic hydrogels’, J. Chem. Phys.104 8713-8720 (1996). [35] S. Sircar, J. P. Keener and A.L. Fogelson, ‘The effect of divalent vs. monovalent ions on the swelling of Mucin-like polyelectrolyte gels: Governing equations and equilibrium analysis’, J. Chem. Phys., 138 014901 (2013). 18 [36] R. Kuver, T. Wong, J. H. Klinkspoor and S. P. Lee, ‘Absence of CFTR is associated with pleiotropic effects on mucins in mouse gallbladder epithelial cells’, Am. J. Physiol. Gastrointest Liver Physiol. 291, G1148-G1154 (2006). [37] M. J. Holmes, N. G. Parker and M. J. W. Povey, ‘Temperature dependence of bulk viscosity in water using acoustic spectroscopy’, J. Phys.: Conf. Ser., 269 012011 1-7 (2011). [38] Y. Mimura, Y. Inoue, N. J. Maeji and R. Chujo, ‘N.m.r. study on conformation of Ac-Thr(αGalNAc)-Ala-Ala-OMe as a model for mucin type glycoprotein’, Int. J. Peptide Protein Res. 34 363-368 (1989). [39] A. F. M. Barton, ‘Handbook of polymer-liquid interaction parameters and solubility parameters’, CRC press (1990). [40] P. Verdugo, ‘Goblet cells secretion and mucogenesis’, Annu. Rev. Physiol., 52 157-176 (1990). [41] J. P. Keener, S. Sircar and A. Fogelson ‘Kinetics of swelling gels’, SIAM J. Appl. Math. 71(3), 854-875 (2011). 19

© Copyright 2017 ExploreDoc