Abstract: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 Novel Modeling and Simulation Approach for the Hindered Mobility of Charged Particles in Biological Hydrogels

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

Mathematical Modeling of Microscale Biology: Ion Pairing, Dielectric Decrement, and Born Energy in Glycosaminoglycan Brushes

Establishment of a Physical Model for Solute Diffusion in Hydrogel: Understanding the Diffusion of Proteins in Poly(sulfobetaine Methacrylate) Hydrogel

Active Diffusion of Self-Propelled Particles in Semi-Flexible Polymer Networks

Asymptotically consistent and computationally efficient modeling of short-ranged molecular interactions between curved slender fibers undergoing large 3D deformations

A computational model of cell polarization and motility coupling mechanics and biochemistry

Multiscale Modeling of Polymer Gels-Chemo-Electric Model versus Discrete Element Model

A Computational Model for Molecular Interactions Between Curved Slender Fibers Undergoing Large 3D Deformations With a Focus on Electrostatic, van der Waals and Repulsive Steric Forces

A moving grid finite element method applied to a mechanobiochemical model for 3D cell migration

Diffusiophoresis of a highly charged rigid colloid in a hydrogel incorporating ion steric interactions

Modeling of spatiotemporal dynamics of ligand-coated particle flow in targeted drug delivery processes

Modeling swelling of pH-responsive microgels: Theory and simulations

Soft nanofluidics governing minority ion exclusion in charged hydrogels.

Sorption and Spatial Distribution of Protein Globules in Charged Hydrogel Particles

Modeling Gel Swelling Equilibrium in Mean-Field: From explicit Models to Poisson-Boltzmann

An Extended Filament Based Lamellipodium Model Produces Various Moving Cell Shapes in the Presence of Chemotactic Signals

Phase Field Model Simulations of Hydrogel Dynamics under Chemical Stimulation

Collective migration under hydrodynamic interactions -- a computational approach

Modeling the Photoinduced Reconfiguration and Directed Motion of Polymer Gels

Modeling and Simulation of the Swelling Behavior of pH-Stimulus-Responsive Hydrogels