NON-EQUILIBRIUM THERMODYNAMICS OF HETEROGENEOUS GROWING BIOSYSTEMS Natalya Kizilova Department of Theoretical and Applied Mechanics Kharkov National University Ukraine
Outline 1. Biological growth: definition, types, properties 2. Experiments with growing plant materials (leaves) 3. Experiment-based mathematical model of growing continuum. Parameter identification. 4. Biological growth in tissue engineering. Experimental technologies and models. 5. A mixture model of the inhomogeneous growing tissue. Application to the tissue growth in the degradable scaffold 6. Conclusions
Growth = irreversible changes in the mass (volume, size) of an object provided by new mass accumulation Tissues=cells + extracellular solid matter + interstitial liquid Plant cells = immovable cells + rigid cellular walls Animal cells = movable (migrating) cells + extracellular solids and liquids I: cell growth and divisions II: extracellular matter production and self-assembling
Biosystems are • open TD systems with are in permanent mass and energy exchange with environment (circulatory, respiratory, excretory systems; outer and internal surfaces) • in permanent non-equilibrium (NE) state working against equilibrium; supporting non-zero gradients and corresponding fluxes; exhibiting complex cross-related phenomena • non-uniform systems (cell types, gradient fields) at permanent dynamical loading (gravity, muscle contractions, flow oscillations, electric impulses) • active systems (parameter-dependent properties; local chemical and mechanical + central nervous and humoral systems) • optimal systems possessing maximal performance at given conditions (minimal energy expenses/entropy production)
Growth types:
Surface growth • Mass accumulation/resorbtion at external surfaces • Coupling of dissolution-crystallization • Driven by c , ,... a e • Features: growth anisotropy; non- uniformity • TD consideration: solidification fronts • Examples: bones, skull, tree trunks, branches, shoots
Inner growth (remodeling) • Mass increase/decrease in each point • Non-zero stress field • Examples: plant leaves and roots, inner organs, tumors • Features: anisotropic growth; residual stresses
Volume growth • Mass increase/decrease in each point • Non-zero stress field • Examples: plant leaves and roots, inner organs, tumors • Features: anisotropic growth; residual stresses
Experimental study of plant leaf growth at zero stress conditions v a ( t ) r r r v 0 v ( v , v ) r v a ( t ) r r r v a ( t ) r
Experimental study of plant leaf growth at mechanical restrictions
Leaf blade deflection and boundary angle measurements
Experiment-based conclusions: • Extraction/compression stimulates/oppresses growth in the corresponding direction • Growth rate at zero-stress conditions is a function of time and concentrations of growth factors/regulators • Growth rate at nonzero-stress conditions is a function of stress tensor components • New material accumulates according to principals of the stress tensor providing the lightweight design • Stress-induced elongation of cells (endothelial cells in vessel wall, skeletal muscle cells, conducting vessels)
Mathematical modeling of growing continua div( v) q ˆ ˆ ˆ 1 ˆ ˆ ˆ e A(t) B (E) d / dt t 1 v v ˆ div 0 i k ˆ e 2 x x k i 0, v 0 n 1 2 v x A xx 2 2 x 2 A A A yy xy xx 2 v 2 2 x y y y x A yy y v v y x 2 A xy y x 2 2 2 A 2 2 A A jj ij ii 2 , i , j 1 , 2 , 3 A , A ii ij 2 2 2 x x x x x x x i j i j i j i
Growth viscosity tensor, Beltrami-Michell equations, growth problem formulation 2 2 e A B ( B B B ) ( B B B ) ii ii ij jj ik kk ji ii jj jj jk kk ik ik iklm lm 2 2 x x B B B 0 0 0 11 12 13 j i B B B 0 0 0 2 21 22 23 ( B ) mm ij B B B 0 0 0 2 31 32 33 B x x i j 0 0 0 B 0 0 44 , i , j , k 1 , 2 , 3 ; m i 3 0 0 0 0 B 0 55 0 0 0 0 0 B 66 3 v 1 v 1 v j m i A A mm ij x b x x 2B x x i m j pp j i m 1 1 v v i k A F 0 ik i x 2B x x k qq k i div F 0 ˆ i,j,k 1 ,2,3, q 9 i k, p 9 i j B B B 11 12 13 * b b B B B det b v 0 ik 21 22 23 ik n B B B 31 32 33
Conclusions • In spite of different shape, size, physiology, evolutionary age, etc… the narrow limits for growth parameters have been found * ~ 0.03 – 0.05 MPa * A ~ 0.5 3 mm / day * 1 B ~ 0.1 1(Pa s) • Transportation systems have the same principles of design (dependences between the lengths, diameters, branching angles, drianage areas) which corresponds to the model of optimal pipeline providing homogenous flow delivery at minimum energy expences.
Biological growth in tissue engineering
Successful laboratory and clinical reports on tissue engineering of: blood and lymphatic vessels [Shin’oka T., et al, 2001] heart valves [Sodian R., et al, 2000] cardiac tissue [Carrier R.L., et al, 1999] bone and cartilage [Vacanti C.A., et al, 1994] tendon [Cao D., et al, 2006] skin [Parenteau N.L., et al, 1991] liver [Kim T.H., et al, 2000] stomach [Maemura T., et al, 2003] intestine [Choi R.S., et al, 1998] bladder [Oberpenning F., et al, 1998] skeletal muscle [Geris L., et al, 2001] nerves [Fansa H., et al, 2003]
3D tissue and organ printing
Polymer and metal scaffolds with regular structure • Role of geometry (strength, lightweight design, porosity, shape of pores, adequate pore sizes for easy penetration of the growing cells/structures) • Role of material (biocompatibility and non-toxicity; controlled degradation kinetics corresponding to the new tissue formation).
Diffusion models of growth C J C C t b J Cb b t J D b b b J D C f(b) b C C b b D D ( ,b, ,R, , (T)) F C C R 6 R
Particle dynamic growth models ( j) 2 ( j) ( j) ( j) (k) d r dr r r dr ( j) ( j) rm m k (r r ) 6 R k a S r 2 n dt dt dt ( j) (k) r r k j adhesion drag repulsion random walks
A multi-phase model of the growing inhomogeneous tissue Solid phases: 1 – cells of different types 2 - vessel walls, connective tissues, airways 3 – extracellular matrix Liquid phases: 4 – intracellular liquids 5 – extracellular (tissue) liquid 6 – delivering liquid Components: 1 - nutrition (glucose, O 2 , …) 2 – growth factors, …
Mass balance equations div v t dC divJ M k dt J C (v v)
Momentum balance equations k v k j kj k k k ( v v ) p R M f j j t k k M v k k (R M ) 0
Energy balance equations E k j k k ( v E ) Q v R f N W j j k t N E k (v R N W ) 0 k
Additional equations (active movement, structure formation, aggregation, …) (v ) (C ,v ,...) (C ,v ,...) t n 1 div(nv ) G t
Recommend
More recommend