A mechanochemical 3D continuum model for smooth muscle contraction under finite strains.

This paper presents a modelling framework in which the mechanochemical properties of smooth muscle cells may be studied. The activation of smooth muscles is considered in a three-dimensional continuum model which is key to realistically capture the function of hollow organs such as blood vessels. On the basis of a general thermodynamical framework the mechanical and chemical phases are specialized in order to quantify the coupled mechanochemical process. A free-energy function is proposed as the sum of a mechanical energy stored in the passive tissue, a coupling between the mechanical and chemical kinetics and an energy related purely to the chemical kinetics and the calcium ion concentration. For the chemical phase it is shown that the cross-bridge model of Hai and Murphy [1988. Am. J. Physiol. Cell Physiol. 254, C99-C106] is included in the developed evolution law as a special case. In order to show the specific features and the potential of the proposed continuum model a uniaxial extension test of a tissue strip is analysed in detail and the related kinematics and stress-stretch relations are derived. Parameter studies point to coupling phenomena; in particular the tissue response is analysed in terms of the calcium ion level. The model for smooth muscle contraction may significantly contribute to current modelling efforts of smooth muscle tissue responses.

In vivo estimation of the contribution of elastin and collagen to the mechanical properties in the human abdominal aorta: effect of age and sex.

The mechanical properties of the aorta affect cardiac function and are related to cardiovascular morbidity/mortality. This study was designed to evaluate the isotropic (mainly elastin, elastin(iso)) and anisotropic (mainly collagen, collagen(ani)) material parameters within the human aorta in vivo. Thirty healthy men and women in three different age categories (23-30, 41-54, and 67-72 yr) were included. A novel mechanical model was used to identify the mechanical properties and the strain field with aid of simultaneously recorded pressure and radius in the abdominal aorta. The magnitudes of the material parameters relating to both the stiffness of elastin(iso) and collagen(ani) were in agreement with earlier in vitro studies. The load-bearing fraction attributed to collagen(ani) oscillated from 10 to 30% between diastolic and systolic pressures during the cardiac cycle. With age, stiffness of elastin(iso) increased in men, despite the decrease in elastin content that has been found due to elastolysis. Furthermore, an increase in stiffness of collagen(ani) at high physiological pressure was found. This might be due to increased glycation, as well as changed isoforms of collagen in the aortic wall with age. A marked sex difference was observed, with a much less age-related effect, both on elastin(iso) and collagen(ani) stiffness in women. Possible factors of importance could be the effect of sex hormones, as well as differing collagen isoforms, between the sexes.

Smooth muscle contraction: mechanochemical formulation for homogeneous finite strains.

Chemical kinetics of smooth muscle contraction affect mechanical properties of organs that function under finite strains. In an effort to gain further insight into organ physiology, we formulate a mechanochemical finite strain model by considering the interaction between mechanical and biochemical components of cell function during activation. We propose a new constitutive framework and use a mechanochemical device that consists of two parallel elements: (i) spring for the cell stiffness; (ii) contractile element for the sarcomere. We use a multiplicative decomposition of cell elongation into filament contraction and cross-bridge deformation, and suggest that the free energy be a function of stretches, four variables (free unphosphorylated myosin, phosphorylated cross-bridges, phosphorylated and dephosphorylated cross-bridges attached to actin), chemical state variable driven by Ca2+-concentration, and temperature. The derived constitutive laws are thermodynamically consistent. Assuming isothermal conditions, we specialize the mechanical phase such that we recover the linear model of Yang et al. [2003a. The myogenic response in isolated rat cerebrovascular arteries: smooth muscle cell. Med. Eng. Phys. 25, 691-709]. The chemical phase is also specialized so that the linearized chemical evolution law leads to the four-state model of Hai and Murphy [1988. Cross-bridge phosphorylation and regulation of latch state in smooth muscle. Am. J. Physiol. 254, C99-C106]. One numerical example shows typical mechanochemical effects and the efficiency of the proposed approach. We discuss related parameter identification, and illustrate the dependence of muscle contraction (Ca2+-concentration) on active stress and related stretch. Mechanochemical models of this kind serve the mathematical basis for analyzing coupled processes such as the dependency of tissue properties on the chemical kinetics of smooth muscle.

Modeling initial strain distribution in soft tissues with application to arteries.

A general theory for computing and identifying the stress field in a residually stressed tissue is presented in this paper. The theory is based on the assumption that a stress free state is obtained by letting each point deform independently of its adjacent points. This local unloading represents an initial strain, and can be described by a tangent map. When experimental data is at hand in a specific situation, the initial strain field may be identified by stating a nonlinear minimization problem where this data is fitted to its corresponding model response. To illustrate the potential of such a method for identifying initial strain fields, the application to an in vivo pressure-radius measurement for a human aorta is presented. The result shows that the initial strain is inconsistent with the strain obtained with the opening-angle-method. This indicates that the opening-angle-method has a too restrictive residual strain parameterization, in this case.

Aorta in vivo parameter identification using an axial force constraint.

It was shown in a previous study by Stålhand et al. (2004) that both material and residual strain parameters for an artery can be identified noninvasively from an in vivo clinical pressure-diameter measurement. The only constraints placed on the model parameters in this previous study was a set of simple box constraints. More advanced constraints can also be utilized, however. These constraints restrict the model parameters implicitly by demanding the state of the artery to behave in a specified way. It has been observed in vitro that the axial force is nearly invariant to the pressure at the physiological operation point. In this paper, we study the possibility to include this behaviour as a constraint in the parameter optimization. The method is tested on an in vivo obtained pressure-diameter cycle for a 24-year-old human. Presented results show that the constrained parameter identification procedure proposed here can be used to obtain good results, and we believe that it may be applied to account for other observed behaviours as well.

Towards in vivo aorta material identification and stress estimation.

This paper addresses the problem of constructing a mechanical model for the abdominal aorta and calibrating its parameters to in vivo measurable data. The aorta is modeled as a pseudoelastic, thick-walled, orthotropic, residually stressed cylindrical tube, subjected to an internal pressure. The model parameters are determined by stating a minimization problem for the model pressure and computing the optimal solution by a minimization algorithm. The data used in this study is in vivo pressure-diameter data for the abdominal aorta of a 24-year-old man. The results show that the axial, circumferential and radial stresses have magnitudes in the span 0 to 180 kPa. Furthermore, the results show that it is possible to determine model parameters directly from in vivo measurable data. In particular, the parameters describing the residual stress distribution can be obtained without interventional procedures.