Model development and numerical simulation of electric-stimulus-responsive hydrogels subject to an externally applied electric field

Based on a multi-phasic mixture theory with consideration of ionic diffusion and convection, a multi-physic model, called the multi-effect-coupling electric-stimulus (MECe) model, is developed for simulation of responsive behavior of the electric-sensitive hydrogels when they are immersed into a bathing solution subject to an externally applied electric field. In the developed model, with chemo-electro-mechanical coupling effects, the convection-diffusion equations for concentration distribution of diffusive ions incorporate the influence of electric potential. The electroneutrality condition is replaced by the Poisson equation for distribution of electric potential. The steady and transient analyses of hydrogel deformation are easily carried out by the continuity and momentum equations of the mixture phase. Further, the computational domain of the present model covers both the hydrogel and the surrounding solution. In order to solve the present mathematical model consisting of multi-field coupled nonlinear partial differential governing equations, a hierarchical iteration technique is proposed and a meshless Hermite-Cloud method (HCM) is employed. The steady-state simulation of the electric-stimulus responsive hydrogel is numerically conducted when it is subjected to an externally applied electric field. The hydrogel deformation and the ionic concentrations as well as electric potentials of both the hydrogel and external solution are investigated. The parameter influences on the swelling behaviors of the hydrogel are also discussed in detail. The simulating results are in good agreement with the experimental data and they validate the presently developed model.     

A network thermodynamic method for numerical solution of the Nernst-Planck and Poisson equation system with application to ionic transport through membranes

Simple techniques of network thermodynamics are used to obtain the numerical solution of the Nernst-Planck and Poisson equation system. A network model for a particular physical situation, namely ionic transport through a thin membrane with simultaneous diffusion, convection and electric current, is proposed. Concentration and electric field profiles across the membrane, as well as diffusion potential, have been simulated using the electric circuit simulation program, SPICE. The method is quite general and extremely efficient, permitting treatments of multi-ion systems whatever the boundary and experimental conditions may be.     

Adsorption of ions at charged sites and phase boundary potentials

The Poisson equation is used to determine for the equilibrium case the space profile of the electric potential at the boundary of two phases, one of which contains anionic sites able to adsorb selectively available cations. The solution has given an expression for the effective width of the space profile of the phase boundary potential in terms of the concentrations and the adsorption characteristics of the ions. It is shown that the sign of one (lumped) parameter determines the sign of the potential. It is also shown that a value of 66 Å can be obtained for the effective width of the space profile of the phase boundary potential. The implications concerning the nature of the “membrane” are briefly discussed.     

Electroviscosity of polyelectrolyte solutions

A theoretical expression for the electroviscous effect in polyelectrolyte solutions, caused by the distortion of counterion-distribution and counterion flow around a polyion under a velocity gradient of solvent flow, was obtained to elucidate the characteristic behaviour of the viscosity of highly charged polyelectrolyte solutions observed at low salt concentration. The derivation of the theory was performed on the basis of the Navier-Stokes-Onsager equation, Poisson equation, and diffusion equations for low molecular ions by the use of a cell model (free-volume model) for a polyion. Energy dissipation was obtained without directly solving these equations. It was found that the derived expression of viscosity explained the experimental results satisfactorily, and that the streaming potential effect caused by the counterion flow played an essential role in the increase in viscosity of polyelectrolyte solutions at finite polymer concentration and low salt concentration ranges.     

Hybrid MD-Nernst Planck model of α-hemolysin conductance properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an α-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson–Nernst–Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K⁺ and Cl^?) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1?M KCl solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5–7 times in comparison to bulk values. Significant statistical variations (17–45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240?mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius ～9?Å with two constriction blocks where the radius is reduced to ～6?Å. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the α-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.     

Isoelectric focusing using non-amphoteric buffers in free solution: I. Determination of stable concentration profiles

For large-scale separations of proteins, the use of simple non-amphoteric buffers in free solution and in multicompartment electrolyzers seems promising for industrial applications. The stabilization of a pH profile with this type of buffer requires the strict observation of two conditions: choice of an adequate buffer; stationary profiles of concentrations. During electrolysis in free solution, the ions of the buffer are displaced across the compartments by migration and by diffusion. To keep a stationary composition, the inflow and outflow of all individual ionic species through each compartment must be identical. At high current, diffusion may be neglected against migration and the ionic flows will be identical if the transport number of each ion is constant at each location within the cell. In these conditions, stationary compositions will be independent of the electric current. This condition of constant transport numbers implies the use of profiles of buffer concentrations different from those published up to now. The new equations for these profiles of concentrations are given in the present paper. The constant migration of the ions must be compensated in the end compartments of the isoelectric focusing cell to provide a stable steady state. Two methods are proposed in the literature: the buffer renewal method and the external recycling method (rheoelectrolysis). Here modified buffer renewal method is proposed. Using stationary mass balances, analytical equations are given to calculate the flows and the composition of the solutions to be recycled or added. Using these equations and the profiles of concentrations to keep constant transport numbers, it is demonstrated that only a renewal of the buffers in the end compartments may lead to stable pH profiles and thus to valid conditions of separation.     

Intrinsic electric fields and proton diffusion in immobilized protein membranes. Effects of electrolytes and buffers.

We present a theory for proton diffusion through an immobilized protein membrane perfused with an electrolyte and a buffer. Using a Nernst-Planck equation for each species and assuming local charge neutrality, we obtain two coupled nonlinear diffusion equations with new diffusion coefficients dependent on the concentration of all species, the diffusion constants or mobilities of the buffers and salts, the pH-derivative of the titration curves of the mobile buffer and the immobilized protein, and the derivative with respect to ionic strength of the protein titration curve. Transient time scales are locally pH-dependent because of protonation-deprotonation reactions with the fixed protein and are ionic strength-dependent because salts provide charge carriers to shield internal electric fields. Intrinsic electric fields arise proportional to the gradient of an "effective" charge concentration. The field may reverse locally if buffer concentrations are large (greater to or equal to 0.1 M) and if the diffusivity of the electrolyte species is sufficiently small. The "ideal" electrolyte case (where each species has the same diffusivity) reduces to a simple form. We apply these theoretical considerations to membranes composed of papain and bovine serum albumin (BSA) and show that intrinsic electric fields greatly enhance the mobility of protons when the ionic strength of the salts is smaller than 0.1 M. These results are consistent with experiments where pH changes are observed to depend strongly on buffer, salt, and proton concentrations in baths adjacent to the membranes.     

A network model for biofilm development in Escherichia coli K-12

Proteins in any solution with a pH value that differs from their isoelectric point exert both an electric Donnan effect (DE) and colloid osmotic pressure. While the former alters the distribution of ions, the latter forces water diffusion. In cells with highly Cl^--permeable membranes, the resting potential is more dependent on the cytoplasmic pH value, which alters the Donnan effect of cell proteins, than on the current action of Na/K pumps. Any weak (positive or negative) electric disturbances of their resting potential are quickly corrected by chloride shifts. In many excitable cells, the spreading of action potentials is mediated through fast, voltage-gated sodium channels. Tissue cells share similar concentrations of cytoplasmic proteins and almost the same exposure to the interstitial fluid (IF) chloride concentration. The consequence is that similar intra- and extra-cellular chloride concentrations make these cells share the same Nernst value for Cl^-. Further extrapolation indicates that cells with the same chloride Nernst value and high chloride permeability should have similar resting membrane potentials, more negative than -80 mV. Fast sodium channels require potassium levels >20 times higher inside the cell than around it, while the concentration of Cl^- ions needs to be >20 times higher outside the cell. When osmotic forces, electroneutrality and other ions are all taken into account, the overall osmolarity needs to be near 280 to 300 mosm/L to reach the required resting potential in excitable cells. High plasma protein concentrations keep the IF chloride concentration stable, which is important in keeping the resting membrane potential similar in all chloride-permeable cells. Probable consequences of this concept for neuron excitability, erythrocyte membrane permeability and several features of circulation design are briefly discussed.     

The theoretical distributions and diffusivities of small ions in chondroitin sulphate and hyaluronate

The electrostatic interactions between polyionic glycosaminoglycans and small mobile ions are investigated using the Poisson-Boltzmann equation and a rod-in-cell model of the polyelectrolyte. Calculations are made for the range of polyelectrolyte concentrations and buffer compositions for which measurements of ion distributions and diffusivities are reported in a companion paper (Maroudas et al., Biophys. Chem. 32 (1988) 257). We conclude that the distribution of mobile ions is largely determined by the 'far-field' potential and is adequately described by the Poisson-Boltzmann theory and also by more approximate theories such as ideal Donnan or 'condensation' theory. The measured variations in cation diffusivities, particularly the increase in diffusivity with increasing matrix concentration at low ionic strengths, are predicted qualitatively using an approximate diffusion theory together with the calculated potential fields. However, the same theory applied to anion diffusion gives qualitatively wrong results.     

Irreversible thermodynamic analysis of electrophoretic light scattering experiments

A fluctuation theory for electrolyte solutions is developed based on the coupling between the equations of nonequilibrium thermodynamics and the Poisson equation. The resulting fluctuation theory is applied to the analysis of electrophoretic light scattering. It is shown that in a binary electrolyte solution (two ionic species), the Doppler shift is not determined by the electrical mobility of either ion, but depends instead on the rate of change of transference number with salt concentration. In addition the ionic relaxation time is shown to be proportional to the conductivity of the solution.     

Modeling and characterization of glucose-sensitive hydrogel: Effect of Young's modulus

A multiphysics model concerning the diffusion and enzyme reaction simultaneously is developed in this paper to characterize the equilibrium behavior of the glucose-sensitive hydrogel, which is called the multi-effect-coupling glucose-stimulus (MECglu) model. The responsive behavior of the hydrogel in the chemo-electro-mechanical coupled energy domains is modeled by the nonlinear coupled partial differential equations. They include the Nernst–Planck equations for the diffusion of mobile species and the enzyme reaction catalyzed by the glucose oxidase and the catalase, the Poisson equation for electric potential, and the mechanical equilibrium equation for finite deformation of the glucose-oxidase-loaded pH-sensitive hydrogel. Numerical simulations demonstrate that the MECglu model can consist well with the published experiment for the practical physiological glucose concentration ranging from 0 to 16.5 mM (300 mg/ml). The effect of Young's modulus of the hydrogel is investigated on the distributive concentrations of reacting and diffusive species and the deformation of the glucose-sensitive hydrogels.     

Flux, coupling, and selectivity in ionic channels of one conformation.

Ions crossing biological membranes are described as a concentration of charge flowing through a selective open channel of one conformation and analyzed by a combination of Poisson and Nernst-Planck equations and boundary conditions, called the PNP theory for short. The ion fluxes in this theory interact much as ion fluxes interact in biological channels and mediated transporters, provided the theoretical channel contains permanent charge and has selectivity created by (electro-chemical) resistance at its ends. Interaction occurs because the flux of different ionic species depends on the same electric field. That electric field is a variable, changing with experimental conditions because the screening (i.e., shielding) of the permanent charge within the channel changes with experimental conditions. For example, the screening of charge and the shape of the electric field depend on the concentration of all ionic species on both sides of the channel. As experimental interventions vary the screening, the electric field varies, and thus the flux of each ionic species varies conjointly, and is, in that sense, coupled. Interdependence and interaction are the rule, independence is the exception, in this channel.     

Immobilized enzymes: Electrokinetic effects on reaction rates under external diffusion

The rates of reactions catalyzed by enzymes immobilized on a nonporous solid surface have been computed employing a Nernst film model. The Nernst-Planck equations for the transport of the charged substrate and product species in the film and the Poisson equation for the distribution of electrical potential are solved numerically with the appropriate boundary conditions. The electrical charge at the surface is assumed to arise from the dissociation equilibria of the acidic and basic surface groups of the enzyme. The pH at the surface affects both the surface charge as well as the intrinsic kinetics of the enzyme-catalyzed reaction. Factors which determine the pH at the surface include the pH in the bulk solution and the release of H(+) ions in the enzyme-catalyzed reaction. The latter causes a lowering of pH at the surface, causing the reaction rate to differ from that computed assuming an equilibrium distribution of electrical potential. Another kind of nonequilibrium contribution is caused by unequal charges or diffusivities of the substrate and products, which results in a diffusion potential being set up. Two moduli are introduced to evaluate the significance of the reaction-generated lowering of pH and the diffusion potential effect. The effect of changing various parameters, e.g., reaction rate constant, substrate concentration, enzyme concentration, pH, etc., on the overall reaction rate are studied.     

Single file ion diffusion and the Schmoluchowski equation

The Schmoluchowski equation is introduced into the problem of single file ion diffusion in a channel. The ions mutually interact due to coulomb repulsion and are also subject to a single ion potential due to the channel. The positions of the ions are represented by a continuous co-ordinate. The problem is reduced to the solution of a pair of transfer integral equations. The resistivity of finite and infinite channels is calculated for various dielectric constants and mean ionic separations. The ionic density for finite channels is also calculated. The results clearly demonstrate that strong coulomb interaction leads to a co-operative motion of the ions across channels.     

The effect of translocating cylindrical particles on the ionic current through a nanopore

The electric field-induced translocation of cylindrical particles through nanopores with circular cross sections is studied theoretically. The coupled Nernst-Planck equations (multi-ion model, MIM) for the concentration fields of the ions in solution and the Stokes equation for the flow field are solved simultaneously. The predictions of the multi-ion model are compared with the predictions of two simplified models based on the Poisson-Boltzmann equation (PBM) and the Smoluchowski's slip velocity (SVM). The concentration field, the ionic current though the pore, and the particle's velocity are computed as functions of the particle's size, location, and electric charge; the pore's size and electric charge; the electric field intensity; and the bulk solution's concentration. In qualitative agreement with experimental data, the MIM predicts that, depending on the bulk solution's concentration, the translocating particle may either block or enhance the ionic current. When the thickness of the electric double layer is relatively large, the PBM and SVM predictions do not agree with the MIM predictions. The limitations of the PBM and SVM are delineated. The theoretical predictions are compared with and used to explain experimental data pertaining to the translocation of DNA molecules through nanopores.     

Sensitivity analysis of the Poisson Nernst–Planck equations: a finite element approximation for the sensitive analysis of an electrodiffusion model

Biological structures exhibiting electric potential fluctuations such as neuron and neural structures with complex geometries are modelled using an electrodiffusion or Poisson Nernst–Planck system of equations. These structures typically depend upon several parameters displaying a large degree of variation or that cannot be precisely inferred experimentally. It is crucial to understand how the mathematical model (and resulting simulations) depend on specific values of these parameters. Here we develop a rigorous approach based on the sensitivity equation for the electrodiffusion model. To illustrate the proposed methodology, we investigate the sensitivity of the electrical response of a node of Ranvier with respect to ionic diffusion coefficients and the membrane dielectric permittivity.

     

Electroosmotic Flow in Nanoscale Parallel-plate Channels: Molecular Simulation Study and Comparison with Classical Poisson–Boltzmann Theory

Molecular dynamics simulations have been carried out for simple electrolyte systems to study the electrokinetically driven osmotic flow in parallel-plate channels of widths ～10–120?nm. The results are compared with the classical theory predictions based on the solution to the Poisson–Boltzmann equation. We find that despite some of the limitations in the Poisson–Boltzmann equation, such as assumption of the Boltzmann distribution for the ions, the classical theory captures the general trend of the variations of the osmotic flow with channel width, as characterized by the mobility of the fluid in channels between ～10 and 120?nm at moderate to low ion concentration. At moderate concentration (corresponding to relatively low surface potential), the classical theory is almost quantitative. The theory and simulation show more disagreement at low concentration, primarily caused by the high surface potential where the assumption of Boltzmann distribution becomes inaccurate. We discuss the limitations of the Poisson–Boltzmann equation as applied to the nanoscale channels.     

Approximation method for electrolyte diffusion through a fixed charge medium

The diffusion of a binary salt through media containing uniformly dispersed immobile ions is formally equivalent to the diffusion of a nonelectrolyte with a concentration-dependent diffusion coefficient. This permits the application of a quasi-steady-state approximation, developed in a previous paper (Bull. Math. Biophysics,21, 19–32, 1959), to problems of transient diffusion. Free diffusion across an interface involving a discontinuity in immobile ion concentration is considered and expressions for the salt concentration and electrical potential at any point are developed. The dependence of electrical potential on the fixed charge density, ionic mobilities, as well as on the boundary concentrations is illustrated by a numerical example. This particular example is suggested by experiments on connective tissue.