Finite element implementation of mechanochemical phenomena in neutral deformable porous media under finite deformation

Biological soft tissues and cells may be subjected to mechanical as well as chemical (osmotic) loading under their natural physiological environment or various experimental conditions. The interaction of mechanical and chemical effects may be very significant under some of these conditions, yet the highly nonlinear nature of the set of governing equations describing these mechanisms poses a challenge for the modeling of such phenomena. This study formulated and implemented a finite element algorithm for analyzing mechanochemical events in neutral deformable porous media under finite deformation. The algorithm employed the framework of mixture theory to model the porous permeable solid matrix and interstitial fluid, where the fluid consists of a mixture of solvent and solute. A special emphasis was placed on solute-solid matrix interactions, such as solute exclusion from a fraction of the matrix pore space (solubility) and frictional momentum exchange that produces solute hindrance and pumping under certain dynamic loading conditions. The finite element formulation implemented full coupling of mechanical and chemical effects, providing a framework where material properties and response functions may depend on solid matrix strain as well as solute concentration. The implementation was validated using selected canonical problems for which analytical or alternative numerical solutions exist. This finite element code includes a number of unique features that enhance the modeling of mechanochemical phenomena in biological tissues. The code is available in the public domain, open source finite element program FEBio (http:∕∕mrl.sci.utah.edu∕software).     

Convection and Diffusion in Charged Hydrated Soft Tissues: A Mixture Theory Approach

The extracellular matrix of cartilage is a charged porous fibrous material. Transport phenomena in such a medium are very complex. In this study, solute diffusive flux and convective flux in porous fibrous media were investigated using a continuum mixture theory approach. The intrinsic diffusion coefficient of solute in the mixture was defined and its relation to drag coefficients was presented. The effect of mechanical loading on solute diffusion in cartilage under unconfined compression with a frictionless boundary condition was analyzed numerically using the model developed. Both strain-dependent hydraulic permeability and diffusivity were considered. Analyses and results show that (1) In porous media, the convective velocity for each solute phase is different. (2) The solute convection in tissue is governed by the relative convective velocity (i.e., relative to solid velocity). (3) Under the assumption that all the frictional interactions among solutes are negligible, the relative convective velocity for α-solute phase is equal to the relative solvent velocity multiplied by its convective coefficient (H ^α) which is also known as the hindrance factor in the literature. The relationship between the convective coefficient and the relative diffusivity of solute is presented. (4) Solute concentration profile within the cartilage sample depends on the phase of dynamic compression.     

Loss of translational entropy in molecular associations

Molecular associations in solution are opposed by the loss of entropy (DeltaS) that results from the restriction of motion of each component in the complex. Theoretical estimates of DeltaS are essential for rationalizing binding affinities, as well as for calculating entropic contribution to enzyme catalysis. Recently a statistical-mechanical framework has been proposed for estimating efficiently the translational entropy loss (DeltaS(trsl)), while taking explicitly into account the complex intermolecular interactions between the solute and the solvent. This framework relates the translational entropy of a solute in solution to its "free volume," defined as the volume accessible to the center of mass of the solute in the presence of the solvent and calculated by using an extension of the cell model (CM) for condensed phases. The translational entropy of pure water, estimated with the CM algorithm, shows good agreement with the experimental information. The free volume of various solutes in water, calculated within the CM by using molecular dynamics simulations with explicit solvent, displays a strong correlation with the solutes' polar and total surface areas. This correlation is used to propose a parameterization that can be used to calculate routinely the translational entropy of a solute in water. We also applied the CM formalism to calculate the free volume and translational entropy loss (DeltaS(trsl)) on binding of benzene to a cavity in a mutant T4-lysozyme. Our results agree with previously published estimates of the binding of benzene to this mutant T4-lysozyme. These and other considerations suggest that the cell model is a simple yet efficient theoretical framework to evaluate the translational entropy loss on molecular association in solution.     

Modeling of neutral solute transport in a dynamically loaded porous permeable gel: implications for articular cartilage biosynthesis and tissue engineering

A primary mechanism of solute transport in articular cartilage is believed to occur through passive diffusion across the articular surface, but cyclical loading has been shown experimentally to enhance the transport of large solutes. The objective of this study is to examine the effect of dynamic loading within a theoretical context, and to investigate the circumstances under which convective transport induced by dynamic loading might supplement diffusive transport. The theory of incompressible mixtures was used to model the tissue (gel) as a mixture of a gel solid matrix (extracellular matrix/scaffold), and two fluid phases (interstitial fluid solvent and neutral solute), to solve the problem of solute transport through the lateral surface of a cylindrical sample loaded dynamically in unconfined compression with frictionless impermeable platens in a bathing solution containing an excess of solute. The resulting equations are governed by nondimensional parameters, the most significant of which are the ratio of the diffusive velocity of the interstitial fluid in the gel to the solute diffusivity in the gel (Rg), the ratio of actual to ideal solute diffusive velocities inside the gel (Rd), the ratio of loading frequency to the characteristic frequency of the gel (f), and the compressive strain amplitude (epsilon 0). Results show that when Rg > 1, Rd < 1, and f > 1, dynamic loading can significantly enhance solute transport into the gel, and that this effect is enhanced as epsilon 0 increases. Based on representative material properties of cartilage and agarose gels, and diffusivities of various solutes in these gels, it is found that the ranges Rg > 1, Rd < 1, correspond to large solutes, whereas f > 1 is in the range of physiological loading frequencies. These theoretical predictions are thus in agreement with the limited experimental data available in the literature. The results of this study apply to any porous hydrated tissue or material, and it is therefore plausible to hypothesize that dynamic loading may serve to enhance solute transport in a variety of physiological processes.     

Computational modeling of chemical reactions and interstitial growth and remodeling involving charged solutes and solid-bound molecules

Mechanobiological processes are rooted in mechanics and chemistry, and such processes may be modeled in a framework that couples their governing equations starting from fundamental principles. In many biological applications, the reactants and products of chemical reactions may be electrically charged, and these charge effects may produce driving forces and constraints that significantly influence outcomes. In this study, a novel formulation and computational implementation are presented for modeling chemical reactions in biological tissues that involve charged solutes and solid-bound molecules within a deformable porous hydrated solid matrix, coupling mechanics with chemistry while accounting for electric charges. The deposition or removal of solid-bound molecules contributes to the growth and remodeling of the solid matrix; in particular, volumetric growth may be driven by Donnan osmotic swelling, resulting from charged molecular species fixed to the solid matrix. This formulation incorporates the state of strain as a state variable in the production rate of chemical reactions, explicitly tying chemistry with mechanics for the purpose of modeling mechanobiology. To achieve these objectives, this treatment identifies the specific theoretical and computational challenges faced in modeling complex systems of interacting neutral and charged constituents while accommodating any number of simultaneous reactions where reactants and products may be modeled explicitly or implicitly. Several finite element verification problems are shown to agree with closed-form analytical solutions. An illustrative tissue engineering analysis demonstrates tissue growth and swelling resulting from the deposition of chondroitin sulfate, a charged solid-bound molecular species. This implementation is released in the open-source program FEBio (www.febio.org). The availability of this framework may be particularly beneficial to optimizing tissue engineering culture systems by examining the influence of nutrient availability on the evolution of inhomogeneous tissue composition and mechanical properties, the evolution of construct dimensions with growth, the influence of solute and solid matrix electric charge on the transport of cytokines, the influence of binding kinetics on transport, the influence of loading on binding kinetics, and the differential growth response to dynamically loaded versus free-swelling culture conditions.     

A mixture theory model of fluid and solute transport in the microvasculature of normal and malignant tissues. I. Theory

In order to better understand the mechanisms governing transport of drugs, nanoparticle-based treatments, and therapeutic biomolecules, and the role of the various physiological parameters, a number of mathematical models have previously been proposed. The limitations of the existing transport models indicate the need for a comprehensive model that includes transport in the vessel lumen, the vessel wall, and the interstitial space and considers the effects of the solute concentration on fluid flow. In this study, a general model to describe the transient distribution of fluid and multiple solutes at the microvascular level was developed using mixture theory. The model captures the experimentally observed dependence of the hydraulic permeability coefficient of the capillary wall on the concentration of solutes present in the capillary wall and the surrounding tissue. Additionally, the model demonstrates that transport phenomena across the capillary wall and in the interstitium are related to the solute concentration as well as the hydrostatic pressure. The model is used in a companion paper to examine fluid and solute transport for the simplified case of an axisymmetric geometry with no solid deformation or interconversion of mass.     

Effect of N-trimethyl chitosan enhancing the dissolution properties of the lipophilic drug cyclosporin A

The aim was to evaluate the influence of N-trimethyl chitosan chloride (TMC) as a carrier for solid dispersion on the dissolution of poorly water-soluble drugs. In this study, we used cyclosporin A(CyA) as a model drug and TMC as a carrier. The effect of various formulation and process variables including TMC-to-CyA mixing weight ratio, weigh molecular(M_w) of TMC and methods used to disperse CyA along with the TMC on the drug dissolution was investigated. The nature of CyA dispersed in the matrix was studied by powder X-ray diffractometry (PXRD), diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and dissolution rate analyses. It was proved that all solid mixtures of CyA with TMCs showed a significantly rapid dissolution rate compared to pure drug and physical mixture. The greater the TMC content the higher the drug dissolution was, up to a maximum corresponding to a polymer: drug ratio of 3:1. The lower the M_w of TMC, the more important the polymer effect was. The dissolution of CyA was remarkably improved by the solid dispersion. The drug dissolution enhancement was attributed to the decreased drug crystallinity and size and polymer wetting effect. There was no significant difference in the efficiency of improving the drug dissolution between the solid dispersions prepared by solvent dispersing and by co-grinding. It was suggested that the TMC with a lower molecular weight is a useful carrier for solid dispersion.     

The physical origin of the low solubility of nonpolar solutes in water

Elementary but general statistical-mechanical relations are derived that relate the thermodynamic properties of the dissolution process to those of the pure solvent. A number of conclusions are drawn from qualitative arguments that these relations suggest. These include the following: (1) The low solubility of nonpolar solutes in water arises not from the fact that water molecules can form hydrogen bonds, but rather from the fact that they are small in size. (2) The large entropy decrease attending the transfer of an inert solute from a nonaqueous solvent to water is largely due to the decrease in entropy of the nonaqueous solvent as the solvent–solvent interaction is restored on removal of the solute from it. (3) It is improper to use values of thermodynamic quantities obtained from small-molecule transfer studies for those that involve macromolecular folding and interaction.     

Validation of theoretical framework explaining active solute uptake in dynamically loaded porous media

Solute transport in biological tissues is a fundamental process necessary for cell metabolism. In connective soft tissues, such as articular cartilage, cells are embedded within a dense extracellular matrix that hinders the transport of solutes. However, according to a recent theoretical study (Mauck et al., 2003, J. Biomech. Eng. 125, 602–614), the convective motion of a dynamically loaded porous solid matrix can also impart momentum to solutes, pumping them into the tissue and giving rise to concentrations which exceed those achived under passive diffusion alone. In this study, the theoretical predictions of this model are verified against experimental measurements. The mechanical and transport properties of an agarose–dextran model system were characterized from independent measurements and substituted into the theory to predict solute uptake or desorption under dynamic mechanical loading for various agarose concentrations and dextran molecular weights, as well as different boundary and initial conditions. In every tested case, agreement was observed between experiments and theoretical predictions as assessed by coefficients of determination ranging from R²=0.61 to 0.95. These results provide strong support for the hypothesis that dynamic loading of a deformable porous tissue can produce active transport of solutes via a pumping mechanisms mediated by momentum exchange between the solute and solid matrix.     

Modeling of heat and mass transfer for solid state fermentation process in tray bioreactor

A mathematical model is developed to simulate oxygen consumption, heat generation and cell growth in solid state fermentation (SSF). The fungal growth on the solid substrate particles results in the increase of the cell film thickness around the particles. The model incorporates this increase in the biofilm size which leads to decrease in the porosity of the substrate bed and diffusivity of oxygen in the bed. The model also takes into account the effect of steric hindrance limitations in SSF. The growth of cells around single particle and resulting expansion of biofilm around the particle is analyzed for simplified zero and first order oxygen consumption kinetics. Under conditions of zero order kinetics, the model predicts upper limit on cell density. The model simulations for packed bed of solid particles in tray bioreactor show distinct limitations on growth due to simultaneous heat and mass transport phenomena accompanying solid state fermentation process. The extent of limitation due to heat and/or mass transport phenomena is analyzed during different stages of fermentation. It is expected that the model will lead to better understanding of the transport processes in SSF, and therefore, will assist in optimal design of bioreactors for SSF.     

A growth mixture theory for cartilage with application to growth-related experiments on cartilage explants

In this paper, we present a growth mixture model for cartilage. The main features of this model are illustrated in a simple equilibrium boundary-value problem that is chosen to illustrate how a mechanical theory of cartilage growth may be applied to growth-related experiments on cartilage explants. The cartilage growth mixture model describes the independent growth of the proteoglycan and collagen constituents due to volumetric mass deposition, which leads to the remodeling of the composition and the mechanical properties of the solid matrix. The model developed here also describes how the material constants of the collagen constituent depend on a scalar parameter that may change over time (e.g., crosslink density); this leads to a remodeling of the structural and mechanical properties of the collagen constituent. The equilibrium boundary-value problem that describes the changes observed in cartilage explants harvested at different stages of a growth or a degenerative process is formulated. This boundary-value problem is solved using existing experimental data for developing bovine cartilage explants harvested at three developmental stages. The solution of the boundary-value problem in conjunction with existing experimental data suggest the types of experimental studies that need to be conducted in the future to determine model parameters and to further refine the model.     

A mixture theory for charged-hydrated soft tissues containing multi-electrolytes: passive transport and swelling behaviors.

A new mixture theory was developed to model the mechano-electrochemical behaviors of charged-hydrated soft tissues containing multi-electrolytes. The mixture is composed of n + 2 constituents (1 charged solid phase, 1 noncharged solvent phase, and n ion species). Results from this theory show that three types of force are involved in the transport of ions and solvent through such materials: (1) a mechanochemical force (including hydraulic and osmotic pressures); (2) an electrochemical force; and (3) an electrical force. Our results also show that three types of material coefficients are required to characterize the transport rates of these ions and solvent: (1) a hydraulic permeability; (2) mechano-electrochemical coupling coefficients; and (3) an ionic conductance matrix. Specifically, we derived the fundamental governing relationships between these forces and material coefficients to describe such mechano-electrochemical transduction effects as streaming potential, streaming current, diffusion (membrane) potential, electro-osmosis, and anomalous (negative) osmosis. As an example, we showed that the well-known formula for the resting cell membrane potential (Hodgkin and Huxley, 1952a, b) could be derived using our new n + 2 mixture model (a generalized triphasic theory). In general, the n + 2 mixture theory is consistent with and subsumes all previous theories pertaining to specific aspects of charged-hydrated tissues. In addition, our results provided the stress, strain, and fluid velocity fields within a tissue of finite thickness during a one-dimensional steady diffusion process. Numerical results were provided for the exchange of Na+ and Ca++ through the tissue. These numerical results support our hypothesis that tissue fixed charge density (CF) plays a significant role in modulating kinetics of ions and solvent transport through charged-hydrated soft tissues.     

A novel technique for measuring solute diffusivities in entrapment matrices used in immobilization

A novel technique has been developed for measuring effective solute diffusivities in entrapment matrices used for cell immobilization. In this technique radiotracers were used to measure effective diffusivities and equilibrium partition coefficients of the solute between the liquid and solid matrix. Ca-alginate was used in this study, because it is one of the most commonly employed matrices for the immobilization of microbial, plant and mammalian cells. The experimental apparatus consisted of a single spherical Ca-alginate bead which was attached to a rotating rod and immersed in water containing C(14)-glucose. The rotational speed of the spherical bead was controlled and resulted in excellent mixing, and negligible external film mass transfer resistance, which allowed the measurement of true effective solute diffusivity within the solid matrix. The rates of C(14)-glucose diffusion within the Ca-alginate sphere were measured using a scintillation spectrometer. A mathematical model of unsteady-state diffusion in a sphere was used with appropriate boundary conditions, and the effective diffusivity of glucose was found from the best fit of the experimental data using a computer regression analysis method. Using 2% (w/v) Ca-alginate beads in this new radiotracer technique the effective diffusivity and partition coefficient of glucose were found to be 6.62 x 10(-10) m(2)/s and 0.98, respectively. The accuracy, advantages, and simplicity of this new method for diffusivity measurements are also compared to other existing methods.     

Development of a system for simultaneous dissolution studies and magnetic resonance imaging of water transport in hydrodynamically balanced systems: A technical note

Summary and Conclusions  Without disturbing the observed processes, the MRI methods combined with the dissolution studies provide insight into the phenomena occurring when the dosage form comes into contact with aqueous fluids. The MR images allow one to observe the solvent penetration into the hydrophilic matrix and the hydrogel formation. The data obtained in the MRI studies complement information obtained from the dissolution studies. The analysis of MR images may support the explanation of differences in the drug-releasing or floating properties of HBS. Published: February 23, 2007     

The Relative Concentration of Solids in the Nucleolus, Nucleus, and Cytoplasm of the Developing Nerve Cell of the Chick

Growing and differentiating nerve cells of the fifth cranial ganglion of the chick embryo were studied by several means. During the period of 70 hours to 11 days of incubation (Hamburger-Hamilton stages 19 to 37) average cell mass increased more than 4.5 times while cells changed from relatively undifferentiated neuroblasts to morphologically characteristic nerve cells with long processes. By making simplifying assumptions about thickness of nucleus and nucleolus, relative to cytoplasmic thickness, it was possible to calculate solute concentration of nucleus and nucleolus relative to that of the cytoplasm from measurements of optical retardations through living cells. Differences in relative solute concentration were observed in nucleolus, cytoplasm, and nucleoplasm in the approximate ratio 1.2:1.0:0.8, respectively. The ratio remained essentially constant during the growth period examined despite the fact that the cell components grow at markedly different rates. This suggests that solid concentrations are physical characteristics of nucleus, nucleolus, and cytoplasm which are maintained even during rapid growth and differentiation. By cytochemical means it was demonstrated that mass increase in the nucleus is not associated with increase in deoxyribonucleic acid. Both ribonucleic acid and protein are in greater concentration in nucleolus and cytoplasm than in the nucleoplasm. Electron microscopy shows interruptions in the nuclear envelope as well as an approximately even distribution of electron density in nucleus and cytoplasm. It is pointed out that consistent differences in solid concentration can exist on either side of the nuclear envelope even though it contains "pores." Implications of these data are discussed.     

Interface Concentrated‐Confinement Suppressing Cathode Dissolution in Water‐in‐Salt Electrolyte

Mass dissolution is one main problems for cathodes in aqueous electrolytes due to the strong polarity of water molecules. In principle, mass dissolution is a thermodynamically favorable process as determined by the Gibbs free energy. However, in real situations, dissolution kinetics, which include viscosity, dissolving mass mobility, and interface properties, are also a critical factor influencing the dissolution rate. Both thermodynamic and kinetic dissolving factors can be regulated by the ratio of salt to solvent in the electrolyte. In this study, concentration‐controlled cathode dissolution is investigated in a susceptible Na₃V₂(PO₄)₃ cathode whose time‐, cycle‐, and state‐of‐charge‐dependent dissolubility are evaluated by multiple electrochemical and chemical methods. It is verified that the super‐highly concentrated water‐in‐salt electrolyte has a high viscosity, low vanadium ion diffusion, low polarity of solvated water, and scarce solute?water dissolving surfaces. These factors significantly lower the thermodynamic‐controlled solubility and the dissolving kinetics via time and physical space local mass interfacial confinement, thereby inducing a new mechanism of interface concentrated‐confinement which improves the cycling stability in real aqueous rechargeable sodium‐ion batteries.     

Implicit solvent models

Implicit solvent models for biomolecular simulations are reviewed and their underlying statistical mechanical basis is discussed. The fundamental quantity that implicit models seek to approximate is the solute potential of mean force, which determines the statistical weight of solute conformations, and which is obtained by averaging over the solvent degrees of freedom. It is possible to express the total free energy as the reversible work performed in two successive steps. First, the solute is inserted in the solvent with zero atomic partial charges; second, the atomic partial charges of the solute are switched from zero to their full values. Consequently, the total solvation free energy corresponds to a sum of non-polar and electrostatic contributions. These two contributions are often approximated by simple geometrical models (such as solvent exposed area models) and by macroscopic continuum electrostatics, respectively. One powerful route is to approximate the average solvent density distribution around the solute, i.e. the solute-solvent density correlation functions, as in statistical mechanical integral equations. Recent progress with semi-analytical approximations makes continuum electrostatics treatments very efficient. Still more efficient are fully empirical, knowledge-based models, whose relation to explicit solvent treatments is not fully resolved, however. Continuum models that treat both solute and solvent as dielectric continua are also discussed, and the relation between the solute fluctuations and its macroscopic dielectric constant(s) clarified.     

Fractionation of macromolecules in an alternating transverse electric field: simulation of the method

An electric field of alternating polarity applied in a direction transverse to the direction of solute transport is used as the basis of a method for the separation of biological macromolecules. The method derives directly from the ability of an electric field to induce movement of a charged macromolecule and from the physics of laminar fluid flow; no adsorptive immobile phase component is involved.

The method is simulated by computer for the case of solute molecules in a solvent flowing through a narrow chamber of recta generates an electric field orthogonal to the direction of solvent flow. Solute molecules repetitively traverse the solvent channel at rates determined by their electrophoretic mobility. During the transit across the channel, solute molecules are transported in the direction of solvent flow; at the channel wall, solvent velocity is negligible and solute transport is limited to that provided by transient diffusion into a mobile solvent zone. Molecules of different intrinsic electrophoretic mobility are separated.

The computer model was used to illustrate the process and to demonstrate the ‘tunability’ of the method as a function of the oscillation frequency and voltage wave form. Because of this tunability, a single instrument can function as the equivalent of several different chromatographic systems. Because fractionation is effected by direct physicochemical phenomena rather than via interaction with chromatographic sites, variations in fractionation results arising from formation of polymers for gel electrophoresis, packing of chromatography columns, or deterioration of columns with use are avoided. This method may be of particular use for the purification of nucleic acid fragments and for the analysis of protei: nucleic acid interactions.     

High shear mixing granulation of ibuprofen and β-cyclodextrin: Effects of process variables on ibuprofen dissolution

The aims of the study were to evaluate the effect of high shear mixer (HSM) granulation process parameters and scale-up on wet mass consistency and granulation characteristics. A mixer torque rheometer (MTR) was employed to evaluate the granulating solvents used (water, isopropanol, and 1:1 vol/vol mixture of both) based on the wet mass consistency. Gral 25 and mini-HSM were used for the granulation. The MTR study showed that the water significantly enhanced the beta-cyclodextrin (βCD) binding tendency and the strength of liquid bridges formed between the particles, whereas the isopropanol/water mixture yielded more suitable agglomerates. Mini-HSM granulation with the isopropanol/water mixture (1:1 vol/vol) showed a reduction in the extent of torque value rise by increasing the impeller speed as a result of more breakdown of agglomerates than coalescence. In contrast, increasing the impeller speed of the Gral 25 resulted in higher torque readings, larger granule size, and consequently, slower dissolution. This was due to a remarkable rise in temperature during Gral granulation that reduced the isopropanol/water ratio in the granulating solvent as a result of evaporation and consequently increased the βCD binding strength. In general, the HSM granulation retarded ibuprofen dissolution compared with the physical mixture because of densification and agglomeration. However, a successful HSM granulation scale-up was not achieved due to the difference in the solvent mixture’s effect from 1 scale to the other.