A theoretical and non-destructive experimental approach for direct inclusion of measured collagen orientation and recruitment into mechanical models of the artery wall

Gradual collagen recruitment has been hypothesized as the underlying mechanism for the mechanical stiffening with increasing stress in arteries. In this work, we investigated this hypothesis in eight rabbit carotid arteries by directly measuring the distribution of collagen recruitment stretch under increasing circumferential loading using a custom uniaxial (UA) extension device combined with a multi-photon microscope (MPM). This approach allowed simultaneous mechanical testing and imaging of collagen fibers without traditional destructive fixation methods. Fiber recruitment was quantified from 3D rendered MPM images, and fiber orientation was measured in projected stacks of images. Collagen recruitment was observed to initiate at a finite strain, corresponding to a sharp increase in the measured mechanical stiffness, confirming the previous hypothesis and motivating the development of a new constitutive model to capture this response. Previous constitutive equations for the arterial wall have modeled the collagen contribution with either abrupt recruitment at zero strain, abrupt recruitment at finite strain or as gradual recruitment beginning at infinitesimal strain. Based on our experimental data, a new combined constitutive model was presented in which fiber recruitment begins at a finite strain with activation stretch represented by a probability distribution function. By directly including this recruitment data, the collagen contribution was modeled using a simple Neo-Hookean equation. As a result, only two phenomenological material constants were required from the fit to the stress stretch data. Three other models for the arterial wall were then compared with these results. The approach taken here was successful in combining stress-strain analysis with simultaneous microstructural imaging of collagen recruitment and orientation, providing a new approach by which underlying fiber architecture may be quantified and included in constitutive equations.     

A multiscale model for red blood cell mechanics

The objective of this article is the derivation of a continuum model for mechanics of red blood cells via multiscale analysis. On the microscopic level, we consider realistic discrete models in terms of energy functionals defined on networks/lattices. Using concepts of Γ-convergence, convergence results as well as explicit homogenisation formulae are derived. Based on a characterisation via energy functionals, appropriate macroscopic stress–strain relationships (constitutive equations) can be determined. Further, mechanical moduli of the derived macroscopic continuum model are directly related to microscopic moduli. As a test case we consider optical tweezers experiments, one of the most common experiments to study mechanical properties of cells. Our simulations of the derived continuum model are based on finite element methods and account explicitly for membrane mechanics and its coupling with bulk mechanics. Since the discretisation of the continuum model can be chosen freely, rather than it is given by the topology of the microscopic cytoskeletal network, the approach allows a significant reduction of computational efforts. Our approach is highly flexible and can be generalised to many other cell models, also including biochemical control.     

A 3D damage model for trabecular bone based on fabric tensors

Motivated by the role of damage in normal and pathological conditions of trabecular bone, a novel 3D constitutive law was developed that describes anisotropic elasticity and the rate-independent degradation in mechanical properties resulting from the growth of cracks or voids in the trabecular tissue. The theoretical model was formulated within the framework of continuum damage mechanics and based on two fabric tensors characterizing the local trabecular morphology. Experimental validation of the model was achieved by uniaxial and torsional testing of waisted bovine trabecular bone specimens. Strong correlations were found between cumulated permanent strain, reduction in elastic moduli and nonlinear postyield stress which support the hypothesis that these variables reflect the same underlying damage process.     

The micromechanical environment of intervertebral disc cells determined by a finite deformation,anisotropic, and biphasic finite element model

Cellular response to mechanical loading varies between the anatomic zones of the intervertebral disc. This difference may be related to differences in the structure and mechanics of both cells and extracellular matrix, which are expected to cause differences in the physical stimuli (such as pressure, stress, and strain) in the cellular micromechanical environment. In this study, a finite element model was developed that was capable of describing the cell micromechanical environment in the intervertebral disc. The model was capable of describing a number of important mechanical phenomena: flow-dependent viscoelasticity using the biphasic theory for soft tissues; finite deformation effects using a hyperelastic constitutive law for the solid phase; and material anisotropy by including a fiber-reinforced continuum law in the hyperelastic strain energy function. To construct accurate finite element meshes, the in situ geometry of IVD cells were measured experimentally using laser scanning confocal microscopy and three-dimensional reconstruction techniques. The model predicted that the cellular micromechanical environment varies dramatically between the anatomic zones, with larger cellular strains predicted in the anisotropic anulus fibrosus and transition zone compared to the isotropic nucleus pulposus. These results suggest that deformation related stimuli may dominate for anulus fibrosus and transition zone cells, while hydrostatic pressurization may dominate in the nucleus pulposus. Furthermore, the model predicted that micromechanical environment is strongly influenced by cell geometry, suggesting that the geometry of IVD cells in situ may be an adaptation to reduce cellular strains during tissue loading.     

A model for the nonlinear elastic response of large arteries

The nonlinear elastic response of large arteries subjected to finite deformations due to action of biaxial principal stresses, is described by simple constitutive equations. Generalized measures of strain and stress are introduced to account for material nonlinearity. This also ensures the existence of a strain energy density function. The orthotropic elastic response is described via quasi-linear relations between strains and stresses. One nonlinear parameter which defines the measures of strain and stress, and three elastic moduli are assumed to be constants. The lateral strain parameters (equivalent to Poisson's ratios in infinitesimal deformations) are deformation dependent. This dependence is defined by empirical relations developed via the incompressibility condition, and by the introduction of a fifth material parameter. The resulting constitutive model compares well with biaxial experimental data of canine carotid arteries.     

Kinematics framework optimized for deformation, growth, and remodeling in vascular organs

A basic tenant of constitutive theory is that phenomenological relations can be derivable from phenomenological behavior or material tests; and yet, conventional representation formulas, such as those of Rivlin and Fung, fail in this regard because of the choice of kinematical variables. Granted, with these representation formulas a particular constitutive relation may be guessed that fits data, but if the relation is non-unique and cannot be derived de novo from actual and/or hypothetical tests, then such a relation is indeterminable. The representation formula of Rivlin is indeterminable because of excessive covariance or coalignment in the kinematical variables. The representation formula of Fung is indeterminable because the incompressibility constraint is not utilized to reduce the kinematical variables a priori. The proposed kinematics framework succeeds in achieving determinability for hyperelastic materials because, primarily, the kinematical variables have minimal coalignment and dilatation and distortion are separated. Determinability is discussed and demonstrated in the context of hyperelasticity. However, any representation formula, whether it is for visco-elasticity or remodeling or etcetera, will be indeterminable when kinematical variables are highly coaligned and/or are subject to a non-reducible constraint. In other words, conventional kinematical frameworks are non-starters for experimentally determining constitutive representations for soft tissues. For the sake of determinability and/or validity of continuum models of vascular tissue, the proposed framework is needed. Moreover, this framework is optimized to simplify the balance equations for tubular structures.     

A chemo-mechanical constitutive model for transiently cross-linked actin networks and a theoretical assessment of their viscoelastic behaviour

Biological materials can undergo large deformations and also show viscoelastic behaviour. One such material is the network of actin filaments found in biological cells, giving the cell much of its mechanical stiffness. A theory for predicting the relaxation behaviour of actin networks cross-linked with the cross-linker α-actinin is proposed. The constitutive model is based on a continuum approach involving a neo-Hookean material model, modified in terms of concentration of chemically activated cross-links. The chemical model builds on work done by Spiros (Doctoral thesis, University of British Columbia, Vancouver, Canada, 1998) and has been modified to respond to mechanical stress experienced by the network. The deformation is split into a viscous and elastic part, and a thermodynamically motivated rate equation is assigned for the evolution of viscous deformation. The model predictions were evaluated for stress relaxation tests at different levels of strain and found to be in good agreement with experimental results for actin networks cross-linked with α-actinin.     

A theoretically-based experimental approach for identifying vascular constitutive relations

We employ a structurally-motivated phenomenological formulation to identify biomechanical experiments which can be used to determine a vascular constitutive relation directly from data. Large deformations, nonlinear material behavior, load-dependent anisotropy, material heterogeneity and incompressibility are accounted for in the analysis. For purposes of illustration, we outline a procedure for studying elastic arteries wherein the behavior of the media and adventitia is considered separately. This general approach for identifying vascular constitutive relations can be applied to any vessel or airway, however, and should provide certain advantages over previous microstructural or purely phenomenological formulations.     

A poroviscoelastic description of fibrin gels

The mechanical induction of specific cell phenotypes can only be properly controlled if the local stimuli applied to the cells are known as a function of the external applied loads. Finite element analysis of the cell carriers would be one method to calculate these local conditions. Furthermore, the constitutive model of the construct material should be able to describe mechanical events known to be responsible for cell stimulation, such as interstitial fluid flow. The aim of this study was to define a biphasic constitutive model for fibrin, a natural hydrogel often used for tissue engineering but not yet thoroughly characterized. Large strain poroelastic and poroviscoelastic constitutive equations were implemented into a finite element model of a fibrin gel. The parameter values for both formulations were found by either analytically solving equivalent low strain equations, or by optimizing directly the large strain equations based on experimental stress relaxation data. No poroelastic parameters that satisfactorily described the fibrin carrier behaviour could be found, suggesting that network viscoelasticity and fluid-flow time-dependent behaviour must be separately accounted for. It was demonstrated that fibrin can be described as a poroviscoelastic material, but a large strain characterization of the parameter values was necessary. The analytical resolution of the low strain poroviscoelastic equations was, however, accurate enough to serve as a reliable initial condition for further optimization of the parameter values with the large strain formulation.     

A transversely isotropic hyperelastic constitutive model of the PDL. Analytical and computational aspects

This study describes the development of a constitutive law for the modelling of the periodontal ligament (PDL) and its practical implementation into a commercial finite element code. The constitutive equations encompass the essential mechanical features of this biological soft tissue: non-linear behaviour, large deformations, anisotropy, distinct behaviour in tension and compression and the fibrous characteristics. The approach is based on the theory of continuum fibre-reinforced composites at finite strain where a compressible transversely isotropic hyperelastic strain energy function is defined. This strain energy density function is further split into volumetric and deviatoric contributions separating the bulk and shear responses of the material. Explicit expressions of the stress tensors in the material and spatial configurations are first established followed by original expressions of the elasticity tensors in the material and spatial configurations. As a simple application of the constitutive model, two finite element analyses simulating the mechanical behaviour of the PDL are performed. The results highlight the significance of integrating the fibrous architecture of the PDL as this feature is shown to be responsible for the complex strain distribution observed.     

A Transversely Isotropic Hyperelastic Constitutive Model of the PDL. Analytical and Computational Aspects

This study describes the development of a constitutive law for the modelling of the periodontal ligament (PDL) and its practical implementation into a commercial finite element code. The constitutive equations encompass the essential mechanical features of this biological soft tissue: non-linear behaviour, large deformations, anisotropy, distinct behaviour in tension and compression and the fibrous characteristics. The approach is based on the theory of continuum fibre-reinforced composites at finite strain where a compressible transversely isotropic hyperelastic strain energy function is defined. This strain energy density function is further split into volumetric and deviatoric contributions separating the bulk and shear responses of the material. Explicit expressions of the stress tensors in the material and spatial configurations are first established followed by original expressions of the elasticity tensors in the material and spatial configurations. As a simple application of the constitutive model, two finite element analyses simulating the mechanical behaviour of the PDL are performed. The results highlight the significance of integrating the fibrous architecture of the PDL as this feature is shown to be responsible for the complex strain distribution observed.     

Stress distributions and material properties determined in articular cartilage from MRI-based finite strains

The noninvasive measurement of finite strains in biomaterials and tissues by magnetic resonance imaging (MRI) enables mathematical estimates of stress distributions and material properties. Such methods allow for non-contact and patient-specific modeling in a manner not possible with traditional mechanical testing or finite element techniques. Here, we employed three constitutive (i.e. linear Hookean, and nonlinear Neo-Hookean and Mooney-Rivlin) relations with known loading conditions and MRI-based finite strains to estimate stress patterns and material properties in the articular cartilage of tibiofemoral joints. Displacement-encoded MRI was used to determine two-dimensional finite strains in juvenile porcine joints, and an iterative technique estimated stress distributions and material properties with defined constitutive relations. Stress distributions were consistent across all relations, although the stress magnitudes varied. Material properties for femoral and tibial cartilage were found to be consistent with those reported in literature. Further, the stress estimates from Hookean and Neo-Hookean, but not Mooney-Rivlin, relations agreed with finite element-based simulations. A nonlinear Neo-Hookean relation provided the most appropriate model for the characterization of complex and spatially dependent stresses using two-dimensional MRI-based finite strain. These results demonstrate the feasibility of a new and computationally efficient technique incorporating MRI-based deformation with mathematical modeling to non-invasively evaluate the mechanical behavior of biological tissues and materials.     

Rapid Growth of Cartilage Rudiments may Generate Perichondrial Structures by Mechanical Induction

Experimental and theoretical research suggest that mechanical stimuli may play a role in morphogenesis. We investigated whether theoretically predicted patterns of stress and strain generated during the growth of a skeletal condensation are similar to in vivo expression patterns of chondrogenic and osteogenic genes. The analysis showed that predicted patterns of compressive hydrostatic stress (pressure) correspond to the expression patterns of chondrogenic genes, and predicted patterns of tensile strain correspond to the expression patterns of osteogenic genes. Furthermore, the results of iterative application of the analysis suggest that stresses and strains generated by the growing condensation could promote the formation and refinement of stiff tissue surrounding the condensation, a prediction that is in agreement with an observed increase in collagen bundling surrounding the cartilage condensation, as indicated by picro-sirius red staining. These results are consistent with mechanical stimuli playing an inductive or maintenance role in the developing cartilage and associated perichondrium and bone collar. This theoretical analysis provides insight into the potential importance of mechanical stimuli during the growth of skeletogenic condensations.The authors J. H. Henderson, L. de la Fuente have contributed equally to this work.     

An Eulerian/XFEM formulation for the large deformation of cortical cell membrane

Most animal cells are surrounded by a thin layer of actin meshwork below their membrane, commonly known as the actin cortex (or cortical membrane). An increasing number of studies have highlighted the role of this structure in many cell functions including contraction and locomotion, but modelling has been limited by the fact that the membrane thickness (about 1?μm) is usually much smaller than the typical size of a cell (10-100?μm). To overcome theoretical and numerical issues resulting from this observation, we introduce in this paper a continuum formulation, based on surface elasticity, that views the cortex as an infinitely thin membrane that can resists tangential deformation. To accurately model the large deformations of cells, we introduced equilibrium equations and constitutive relations within the Eulerian viewpoint such that all quantities (stress, rate of deformation) lie in the current configuration. A solution procedure is then introduced based on a coupled extended finite element approach that enables a continuum solution to the boundary value problem in which discontinuities in both strain and displacement (due to cortical elasticity) are easily handled. We validate the approach by studying the effect of cortical elasticity on the deformation of a cell adhering on a stiff substrate and undergoing internal contraction. Results show very good prediction of the proposed method when compared with experimental observations and analytical solutions for simple cases. In particular, the model can be used to study how cell properties such as stiffness and contraction of both cytoskeleton and cortical membrane lead to variations in cell's surface curvature. These numerical results show that the proposed method can be used to gain critical insights into how the cortical membrane affects cell deformation and how it may be used as a means to determine a cell's mechanical properties by measuring curvatures of its membrane.     

Quantification of the mechanical behavior of carotid arteries from wild-type, dystrophin-deficient, and sarcoglycan-delta knockout mice

As patients with muscular dystrophy live longer because of improved clinical care, they will become increasingly susceptible to many of the cardiovascular diseases that affect the general population. There is, therefore, a pressing need to better understand both the biology and the mechanics of the arterial wall in these patients. In this paper, we use nonlinear constitutive relations to model, for the first time, the biaxial mechanical behavior of carotid arteries from two common mouse models of muscular dystrophy (dystrophin-deficient and sarcoglycan-delta null) and wild-type controls. It is shown that a structurally motivated four-fiber family stress-strain relation describes the passive behavior of all three genotypes better than does a commonly used phenomenological exponential model, and that a Rachev-Hayashi model describes the mechanical contribution of smooth muscle contraction under basal tone. Because structurally motivated constitutive relations can be extended easily to model adaptations to altered hemodynamics, results from this study represent an important step toward the ultimate goal of understanding better the mechanobiology and pathophysiology of arteries in muscular dystrophy.     

A new constitutive formulation for characterizing the mechanical behavior of soft tissues.

We present a new constitutive formulation that combines certain desirable features of two previously used approaches (phenomenological and microstructural). Specifically, we assume that certain soft tissues can be idealized as composed of various families of noninteracting fibers and a homogeneous matrix. Both the fibers and the matrix are assumed to follow the gross deformation. Within the usual framework of pseudoelasticity, incompressibility, homogeneity, and the continuum hypothesis, a pseudostrain-energy function (W) is proposed wherein W is expressed in terms of matrix and fibrous contributions. Unlike phenomenological approaches where a W is usually chosen in an ad hoc manner, the present approach can be used to postulate reasonable forms of W based on limited structural information and multiaxial stress-strain data. Illustrative applications of the theory are discussed for visceral pleura and myocardium. Concise structurally motivated constitutive relations result, wherein load-dependent anisotropy, nonlinear material behavior, finite deformations, and incompressibility are accounted for.     

Trabecular bone density and loading history: regulation of connective tissue biology by mechanical energy

The method of considering a single loading condition in the study of stress/morphology relationships in trabecular bone is expanded to include the multiple loading conditions experienced by bone in vivo. The bone daily loading histories are characterized in terms of stress magnitudes or cyclic strain energy density and the number of loading cycles. Relationships between local bone apparent density and loading history are developed which assume that bone mass is adjusted in response to strength or energy considerations. Three different bone maintenance criteria are described which are formulated based upon: (1) continuum model effective stress, (2) continuum model fatigue damage accumulation density, and (3) bone tissue strain energy density. These approaches can be applied to predict variations in apparent density within bone and among bones. We show that all three criteria have similar mathematical forms and may be related to the density (or concentration) of bone strain energy which is transferred (dissipated) in the mineralized tissue. The loading history and energy transfer concepts developed here can be applied to many different situations of growth, functional adaptation, injury, and aging of connective tissues.     

NF-kappaB responds to mechanical strains in osteoblast-like cells, and lighter strains create an NF-kappaB response more readily

This study was to examine the early responses of nuclear factor kappa B (NF-kappaB) to mechanical strains in MG-63. MG-63 cells were subjected to cyclic uniaxial compressive or tensile strain, produced by a four-point bending system, at 1000 microstrain or 4000 microstrain for 5 min, 15 min, 30 min and 1h, respectively. Control cells received the same treatment with no mechanical stress loading. Expression of NF-kappaB (p60) was measured by Western blotting. NF-kappaB responded rapidly to mechanical stimuli in MG-63 cells. NF-kappaB was activated by cyclic uniaxial stretch at 1000 microstrain while it was restrained under a compressive strain environment at 1000 microstrain (P<0.001). The effects reversed for tension and compression at 4000 microstrain (P<0.001). Furthermore, strains at 1000 microstrain affected NF-kappaB expression much easier than those at 4000 microstrain. This indicates that there may be different responding mechanisms or mechanotransduction pathways for different mechanical stimuli.     

A constitutive model for vascular tissue that integrates fibril, fiber and continuum levels with application to the isotropic and passive properties of the infrarenal aorta

A fundamental understanding of the mechanical properties of the extracellular matrix (ECM) is critically important to quantify the amount of macroscopic stress and/or strain transmitted to the cellular level of vascular tissue. Structural constitutive models integrate histological and mechanical information, and hence, allocate stress and strain to the different microstructural components of the vascular wall. The present work proposes a novel multi-scale structural constitutive model of passive vascular tissue, where collagen fibers are assembled by proteoglycan (PG) cross-linked collagen fibrils and reinforce an otherwise isotropic matrix material. Multiplicative kinematics account for the straightening and stretching of collagen fibrils, and an orientation density function captures the spatial organization of collagen fibers in the tissue. Mechanical and structural assumptions at the collagen fibril level define a piece-wise analytical stress-stretch response of collagen fibers, which in turn is integrated over the unit sphere to constitute the tissue's macroscopic mechanical properties. The proposed model displays the salient macroscopic features of vascular tissue, and employs the material and structural parameters of clear physical meaning. Likewise, the constitutive concept renders a highly efficient multi-scale structural approach that allows for the numerical analysis at the organ level. Model parameters were estimated from isotropic mean-population data of the normal and aneurysmatic aortic wall and used to predict in-vivo stress states of patient-specific vascular geometries, thought to demonstrate the robustness of the particular Finite Element (FE) implementation. The collagen fibril level of the multi-scale constitutive formulation provided an interface to integrate vascular wall biology and to account for collagen turnover.