Comparison of analytical and inverse finite element approaches to estimate cell viscoelastic properties by micropipette aspiration

The viscoelastic properties of cells are important in predicting cell deformation under mechanical loading and may reflect cell phenotype or pathological transition. Previous studies have demonstrated that viscoelastic parameters estimated by finite element (FE) analyses of micropipette aspiration (MA) data differ from those estimated by the analytical half-space model. However, it is unclear whether these differences are statistically significant, as previous studies have been based on average cell properties or parametric analyses that do not reflect the inherent experimental and biological variability of real experimental data. To determine whether cell material parameters estimated by the half-space model are significantly different from those predicted by the FE method, we implemented an inverse FE method to estimate the viscoelastic parameters of a population of primary porcine aortic valve interstitial cells tested by MA. We found that inherent differences between the analytical and inverse FE estimation methods resulted in statistically significant differences in individual cell properties. However, in cases with small pipette to cell radius ratios and short loading periods, model-dependent differences were masked by experimental and cell-to-cell variability. Analytical models that account for finite cell-size and loading rate further relaxed the experimental conditions for which accurate cell material parameter estimates could be obtained. These data provide practical guidelines for analysis of MA data that account for the wide range of conditions encountered in typical experiments.     

Flow in an arteriosclerotic blood vessel with rigid permeable wall

This paper concerns the fluid-mechanical study of the effects of permeability of the wall of an arteriosclerotic blood vessel by idealizing the tissue space as a porous medium bounding the blood vessel and the arteriosclerotic blood vessel as a constricted axisymetric tube of slowly but arbitrarily varying cross-secton. An analytical solution in the general case is obtained by perturbation technique and several important particular geometries of constriction are discussed as special cases. Numerical results for the effects of permeability on the wall shear stress and filtration velocity are shown graphically.     

Viscoelastic properties of cell walls of single living plant cells determined by dynamic nanoindentation

Plant development results from controlled cell divisions, structural modifications, and reorganizations of the cell wall. Thereby, regulation of cell wall behaviour takes place at multiple length scales involving compositional and architectural aspects in addition to various developmental and/or environmental factors. The physical properties of the primary wall are largely determined by the nature of the complex polymer network, which exhibits time-dependent behaviour representative of viscoelastic materials. Here, a dynamic nanoindentation technique is used to measure the time-dependent response and the viscoelastic behaviour of the cell wall in single living cells at a micron or sub-micron scale. With this approach, significant changes in storage (stiffness) and loss (loss of energy) moduli are captured among the tested cells. The results reveal hitherto unknown differences in the viscoelastic parameters of the walls of same-age similarly positioned cells of the Arabidopsis ecotypes (Col 0 and Ws 2). The technique is also shown to be sensitive enough to detect changes in cell wall properties in cells deficient in the activity of the chromatin modifier ATX1. Extensive computational modelling of the experimental measurements (i.e. modelling the cell as a viscoelastic pressure vessel) is used to analyse the influence of the wall thickness, as well as the turgor pressure, at the positions of our measurements. By combining the nanoDMA technique with finite element simulations quantifiable measurements of the viscoelastic properties of plant cell walls are achieved. Such techniques are expected to find broader applications in quantifying the influence of genetic, biological, and environmental factors on the nanoscale mechanical properties of the cell wall.     

Pulse wave velocity as a diagnostic index: the pitfalls of tethering versus stiffening of the arterial wall

Pulse wave velocity (PWV) is often used as a clinical index of aging, vascular disease, or age related hypertension. This practice is based on the assumption that a higher wave speed indicates vascular stiffening. This assumption is well grounded in the physics of pulsatile flow of an incompressible fluid where it is fully established that a pulse wave travels faster in a tube of stiffer wall, the wave speed becoming infinite in the mathematical limit of a rigid wall. However, in this paper we point out that the physical principal of higher pulse wave velocity in a stiffer tube is strictly valid only when the wall is free from outside constraints, which in the physiological setting is present in the form of tethering of the vessel wall. The use of PWV as an index of arterial stiffening may thus lose its validity if tethering is involved. A solution of the problem of vessel wall mechanics as they arise from the physiological pulsatile flow problem is presented for the purpose of resolving this issue. The vessel wall is considered to have finite thickness with or without tethering and with a range of mechanical properties ranging from viscoelastic to stiff. The results show that, indeed, while the wave speed becomes infinite in the mathematical limit of a rigid free wall, the opposite actually happens if the vessel wall is tethered. Here the wave speed actually diminishes as the degree of tethering increases. This dichotomy in the effects of tethering versus stiffening of the arterial wall may clearly lead to error in the interpretation of PWV as an index of vessel wall stiffness. In particular, a normal value of PWV may lead to the conclusion that vessel wall stiffening is absent while this value may in fact have been lowered by tethering. In other words, the diagnostic test may lead to a false negative diagnosis. Our results indicate that the reason for which PWV is lower in a tethered wall compared with that in a free wall of the same stiffness is that the radial movements of the wall are greatly reduced by tethering. More precisely, the results show that PWV depends strongly on the ratio of radial to axial displacements and that this ratio is much lower in a tethered wall than it is in a free wall of the same stiffness.     

Comparative study of viscoelastic arterial wall models in nonlinear one-dimensional finite element simulations of blood flow

It is well known that blood vessels exhibit viscoelastic properties, which are modeled in the literature with different mathematical forms and experimental bases. The wide range of existing viscoelastic wall models may produce significantly different blood flow, pressure, and vessel deformation solutions in cardiovascular simulations. In this paper, we present a novel comparative study of two different viscoelastic wall models in nonlinear one-dimensional (1D) simulations of blood flow. The viscoelastic models are from papers by Holenstein et al. in 1980 (model V1) and Valdez-Jasso et al. in 2009 (model V2). The static elastic or zero-frequency responses of both models are chosen to be identical. The nonlinear 1D blood flow equations incorporating wall viscoelasticity are solved using a space-time finite element method and the implementation is verified with the Method of Manufactured Solutions. Simulation results using models V1, V2 and the common static elastic model are compared in three application examples: (i) wave propagation study in an idealized vessel with reflection-free outflow boundary condition; (ii) carotid artery model with nonperiodic boundary conditions; and (iii) subject-specific abdominal aorta model under rest and simulated lower limb exercise conditions. In the wave propagation study the damping and wave speed were largest for model V2 and lowest for the elastic model. In the carotid and abdominal aorta studies the most significant differences between wall models were observed in the hysteresis (pressure-area) loops, which were larger for V2 than V1, indicating that V2 is a more dissipative model. The cross-sectional area oscillations over the cardiac cycle were smaller for the viscoelastic models compared to the elastic model. In the abdominal aorta study, differences between constitutive models were more pronounced under exercise conditions than at rest. Inlet pressure pulse for model V1 was larger than the pulse for V2 and the elastic model in the exercise case. In this paper, we have successfully implemented and verified two viscoelastic wall models in a nonlinear 1D finite element blood flow solver and analyzed differences between these models in various idealized and physiological simulations, including exercise. The computational model of blood flow presented here can be utilized in further studies of the cardiovascular system incorporating viscoelastic wall properties.     

A study on the non-linear flow of blood through arteries

The pulsatile flow of blood through arteries is investigated in this paper by treating the blood vessel as a thin-walled anisotropic, non-linearly viscoelastic, incompressible circular cylindrical shell; nonlinearities of the flow of blood are also paid due consideration. The displacement components of the vessel wall are obtained from the equations of equilibrium which have been linearized by employing the principle of superimposition of a small deformation on a state of known finite deformation. The influence of the wall deformation on the flow properties of blood, has been accounted for by considering suitably formulated continuity conditions. A finitedifference scheme is employed for solving the flow equations together with the boundary and initial conditions by using the locally measured values of pressure and pressure gradient. Numerical results obtained for the velocity profile of blood flowing in a canine middle descending thoracic aorta have been presented through figures.     

Is arterial wall-strain stiffening an additional process responsible for atherosclerosis in coronary bifurcations?: an in vivo study based on dynamic CT and MRI

Coronary bifurcations represent specific regions of the arterial tree that are susceptible to atherosclerotic lesions. While the effects of vessel compliance, curvature, pulsatile blood flow, and cardiac motion on coronary endothelial shear stress have been widely explored, the effects of myocardial contraction on arterial wall stress/strain (WS/S) and vessel stiffness distributions remain unclear. Local increase of vessel stiffness resulting from wall-strain stiffening phenomenon (a local process due to the nonlinear mechanical properties of the arterial wall) may be critical in the development of atherosclerotic lesions. Therefore, the aim of this study was to quantify WS/S and stiffness in coronary bifurcations and to investigate correlations with plaque sites. Anatomic coronary geometry and cardiac motion were generated based on both computed tomography and MRI examinations of eight patients with minimal coronary disease. Computational structural analyses using the finite element method were subsequently performed, and spatial luminal arterial wall stretch (LW(Stretch)) and stiffness (LW(Stiff)) distributions in the left main coronary bifurcations were calculated. Our results show that all plaque sites were concomitantly subject to high LW(Stretch) and high LW(Stiff), with mean amplitudes of 34.7 ± 1.6% and 442.4 ± 113.0 kPa, respectively. The mean LW(Stiff) amplitude was found slightly greater at the plaque sites on the left main coronary artery (mean value: 482.2 ± 88.1 kPa) compared with those computed on the left anterior descending and left circumflex coronary arteries (416.3 ± 61.5 and 428.7 ± 181.8 kPa, respectively). These findings suggest that local wall stiffness plays a role in the initiation of atherosclerotic lesions.     

Oscillatory flow of blood in a stenosed artery

In order to understand the abnormal flow conditions of blood in a locally constricted blood vessel, the analytical results are obtained for the oscillatory flow of blood which behaves as a Newtonian fluid. It is here assumed that the surface roughness is cosine-shaped and the maximum height of the roughness is very small compared with the radius of the unconstricted tube. Numerical solutions are presented for the instantaneous flow rate, resistive impedance, wall shear stress and phase lag.     

A mechanism for arteriolar remodeling based on maintenance of smooth muscle cell activation

Structural adaptation in arterioles is part of normal vascular physiology but is also seen in disease states such as hypertension. Smooth muscle cell (SMC) activation has been shown to be central to microvascular remodeling. We hypothesize that, in a remodeling process driven by SMC activation, stress sensitivity of the vascular wall is a key element in the process of achieving a stable vascular structure. We address whether the adaptive changes in arterioles under different conditions can arise through a common mechanism: remodeling in a stress-sensitive wall driven by a shift in SMC activation. We present a simple dynamic model and show that structural remodeling of the vessel radius by rearrangement of the wall material around a lumen of a different diameter and driven by differences in SMC activation can lead to vascular structures similar to those observed experimentally under various conditions. The change in structure simultaneously leads to uniform levels of circumferential wall stress and wall strain, despite differences in transmural pressure. A simulated vasoconstriction caused by increased SMC activation leads to inward remodeling, whereas outward remodeling follows relaxation of the vascular wall. The results are independent of the specific myogenic properties of the vessel. The simulated results are robust in the face of parameter changes and, hence, may be generalized to vessels from different vascular beds.     

A method for measuring linearly viscoelastic properties of human tympanic membrane using nanoindentation

A viscoelastic nanoindentation technique was developed to measure both in-plane and through-thickness viscoelastic properties of human tympanic membrane (TM). For measurement of in-plane Young's relaxation modulus, the TM sample was clamped on a circular hole and a nanoindenter tip was used to apply a concentrated force at the center of the TM sample. In this setup, the resistance to nanoindentation displacement can be considered due primarily to the in-plane stiffness. The load-displacement curve obtained was used along with finite element analysis to determine the in-plane viscoelastic properties of TM. For measurements of Young's relaxation modulus in the through-thickness (out-of-plane) direction, the TM sample was placed on a relatively rigid solid substrate and nanoindentation was made on the sample surface. In this latter setup, the resistance to nanoindentation displacement arises primarily due to out-of-plane stiffness. The load-displacement curve obtained in this manner was used to determine the out-of-plane relaxation modulus using the method appropriate for viscoelastic materials. From our sample tests, we obtained the steady-state values for in-plane moduli as approximately 17.4 MPa and approximately 19.0 MPa for posterior and anterior portions of TM samples, respectively, and the value for through-thickness modulus as approximately 6.0 MPa for both posterior and anterior TM samples. Using this technique, the local out-of-plane viscoelastic modulus can be determined for different locations over the entire TM, and the in-plane properties can be determined for different quadrants of the TM.     

Cability of the neuro-muscular spatial relationships in the wall of the blood vessels]

Electronmicroscopic studies have been made on the distribution of the nervous fibers in the wall of the constricted and dilated auricular artery of the rabbit. It was demonstrated that during the dilatation of the vessel, the dimensions of the neuro-muscular cleft decrease. In the constricted artery, 50% of the nervous fibers are found in the adventitial zone which has a thickness of 4.6 mu. In the dilated vessel this zone decreases up to 1.8 mu. It is suggested that this mechanism is responsible for the regulation of the blood vessel tone.     

The effect of blood viscoelasticity on pulsatile flow in stationary and axially moving tubes

An analytical solution for pulsatile flow of a generalized Maxwell fluid in straight rigid tubes, with and without axial vessel motion, has been used to calculate the effect of blood viscoelasticity on velocity profiles and shear stress in flows representative of those in the large arteries. Measured bulk flow rate Q waveforms were used as starting points in the calculations for the aorta and femoral arteries, from which axial pressure gradient ▿P waves were derived that would reproduce the starting Q waves for viscoelastic flow. The ▿P waves were then used to calculate velocity profiles for both viscoelastic and purely viscous flow. For the coronary artery, published ▿P and axial vessel acceleration waveforms were used in a similar procedure to determine the separate and combined influences of viscoelasticity and vessel motion.Differences in local velocities, comparing viscous flow to viscoelastic flow, were in all cases less than about 2% of the peak local velocity. Differences in peak wall shear stress were less than about 3%.In the coronary artery, wall shear stress differences between viscous and viscoelastic flow were small, regardless of whether axial vessel motion was included. The shape of the wall shear stress waveform and its difference, however, changed dramatically between the stationary and moving vessel cases. The peaks in wall shear stress difference corresponded with large temporal gradients in the combined driving force for the flow.     

Viscoelastic properties of relaxed papillary muscle under physiological hypertrophy

Viscoelastic properties of relaxed rat papillary muscles in physiological hypertrophy (intensive swimming for 5 weeks) have been obtained. It has been ascertained that viscoelastic properties of hypertrophied muscles are not significantly distinguished from those of control papillary muscles. A three-dimensional model of myocardial fascicle has been verified in compliance with experimental data of biomechanical tests of hypertrophied muscles. Elastic and viscous parameters of structural elements of the model negligibly differ from the parameters of the model of a control muscle. It is shown that physiological hypertrophy has a slight influence on viscoelastic properties of papillary muscles.     

Loss of stability-a new model for stress relaxation in plant cell walls

This study addresses the mechanism of wall stress relaxation in growing plant cells. The current viscoelastic model of cell wall relaxation, which dates from the work of Preston, Cleland, Lockhart, and others in the 1960s, has serious shortcomings. It has been shown however that the theory of loss of stability (LOS) can be applied to materials in tension, leading to the conclusion that the relaxation of stresses in the walls of any pressure vessel is rigorously modeled using LOS. We propose that LOS also provides a more appropriate and versatile model of stress relaxation in growing plant cells. We argue that when treated as a manifestation of LOS, the regulation of cell turgor has a rigorous and demonstrable basis in the geometrical and physical properties of the cell wall and the cell's ability to import water. Thus plant cell growth can be regarded as an inherently self-limiting process, tunable by biochemical or structural means. Lastly, despite the current limitations of our model, we apply direct measurement of elastic modulus, wall thickness and cell radius obtained from cylindrical Chara corallina cells to generate an initial calculation of critical pressures in a hypothetical spherical cell with the same material properties.     

Platelet Mechanosensing of Collagen Matrices

During vascular injury, platelets adhere to exposed subendothelial proteins, such as collagen, on the blood vessel walls to trigger clot formation. Although the biochemical signalings of platelet-collagen interactions have been well characterized, little is known about the role microenvironmental biomechanical properties, such as vascular wall stiffness, may have on clot formation. To that end, we investigated how substrates of varying stiffness conjugated with the same concentration of Type I collagen affect platelet adhesion, spreading, and activation. Using collagen-conjugated polyacrylamide (PA) gels of different stiffnesses, we observed that platelets do in fact mechanotransduce the stiffness cues of collagen substrates, manifesting in increased platelet spreading on stiffer substrates. In addition, increasing substrate stiffness also increases phosphatidylserine exposure, a key aspect of platelet activation that initiates coagulation on the platelet surface. Mechanistically, these collagen substrate stiffness effects are mediated by extracellular calcium levels and actomyosin pathways driven by myosin light chain kinase but not Rho-associated protein kinase. Overall, our results improve our understanding of how the mechanics of different tissues and stroma affect clot formation, what role the increased vessel wall stiffness in atherosclerosis may directly have on thrombosis leading to heart attacks and strokes, and how age-related increased vessel wall stiffness affects hemostasis and thrombosis.     

Effect of initial stresses on the wave propagation in arteries

A theoretical analysis for the problem of wave propagation in arteries is presented. Blood is treated as a Newtonian, viscous incompressible fluid. On the basis of information derived from experimental investigations on the mechanical properties of wall tissues, the arterial wall is considered to be nonlinearly viscoelastic and orthotropic. The analysis is carried out for a cylindrical artery, under the purview of membrane theory, by taking account the effect of initial stresses. The motion of the wall and that of the fluid are assumed to be axisymmetric. Particular emphasis has been paid to the propagation of small amplitude harmonic waves whose wavelength is large compared to the radius of the vessel. By employing the equations of motion of the fluid and those for the wall, together with the equations of continuity, a frequency equation is derived by exploiting the conditions of continuity of the velocity of the arterial wall and that of blood on the endosteal surface of the wall. In order to illustrate the validity of the derived analytical expressions a quantitative analysis is made for the variations of the phase velocities as well as the transmission coefficient with frequency corresponding to different transmural pressures which relate to different initial stresses. Computational results indicate that phase velocities increase with the increase of transmural pressures.     

Loss of stability: a new look at the physics of cell wall behavior during plant cell growth

In this article we investigate aspects of turgor-driven plant cell growth within the framework of a model derived from the Eulerian concept of instability. In particular we explore the relationship between cell geometry and cell turgor pressure by extending loss of stability theory to encompass cylindrical cells. Beginning with an analysis of the three-dimensional stress and strain of a cylindrical pressure vessel, we demonstrate that loss of stability is the inevitable result of gradually increasing internal pressure in a cylindrical cell. The turgor pressure predictions based on this model differ from the more traditional viscoelastic or creep-based models in that they incorporate both cell geometry and wall mechanical properties in a single term. To confirm our predicted working turgor pressures, we obtained wall dimensions, elastic moduli, and turgor pressures of sequential internodal cells of intact Chara corallina plants by direct measurement. The results show that turgor pressure predictions based on loss of stability theory fall within the expected physiological range of turgor pressures for this plant. We also studied the effect of varying wall Poisson's ratio nu on extension growth in living cells, showing that while increasing elastic modulus has an understandably negative effect on wall expansion, increasing Poisson's ratio would be expected to accelerate wall expansion.     

A viscoelastic model for use in predicting arterial pulse waves

In nonlinear mathematical models of the arterial circulation, the viscoelasticity of the vessel walls has generally been neglected or only taken into account in a highly approximate manner. A new method is proposed to simulate the nonlinear viscoelastic properties of the wall material with the aid of a convolution integral of the creep function and the pressure history. With this simulation it is possible to properly describe the measured characteristics of arterial viscoelasticity. Moreover, it is utilized in a mathematical model of arterial pulse propagation to study the influence of the internal wall friction on the shape, amplitude and mean value of pressure and flow pulses. The corresponding predictions are in much better agreement with in-vivo measurements, especially for the distal part of the circulation, than those obtained without viscoelasticity.     

Analysis of some effects of aortic smooth muscle on<Emphasis Type="Italic">in vivo</Emphasis> aortic pressure curves

When aortic pressure curves were predicted previously on the basis of a newly developed model of visco-elastic properties of the aorta, it was necessary to use published viscoelastic constants. These were usually obtained from longitudinal strips of blood vessels long removed from the animal, and therefore probably containing deteriorated smooth muscle. The predicted curves had the same form as actual tracings, substantiating the analysis somewhat, but the pressure levels were low. These low levels, if due to inadequate visco-elastic constants, could be attributed to the use of longitudinal rather than circumferential segments as well as to the use of segments with deteriorated muscle. The present analysis uses data collected by the author testing circumferential viscoelastic properties of fourteen different aortic regions in a way suggested by the author's model of an aortic wall. Moreover, the constants were measured on segments containing muscle relaxed by EDTA solutions and on similar segments containing muscle contracted by neosynephrine. These visco-elastic constants were used in the author's nonlinear differential equation of motion of the aortic wallin vivo to predictin vivo pressure curves. The predicted curves were low in any given aortic region if relaxed constants were used, but at normal levels with contracted constants. In fact, pressure curves predicted using constants obtained from aortic segments containing contracted muscle resembled actual tracings in form and pressure levels. Even the observed variations in the form of the systolic pressure curve down the aorta were predicted by this analysis.     

Biaxial elastic material properties of porcine coronary media and adventitia

The importance of mechanical stresses and strains has become well recognized in vascular physiology and pathology. To compute the stress and strain on the various components of the vessel wall, we must know the constitutive equations for the different layers of the vessel wall. The objective of the present study is to determine the constitutive equation of the coronary artery treated as a two-layer composite: intima-media and adventitial layers. Twelve hearts were obtained from a local slaughterhouse, and the right coronary artery and left anterior descending artery were dissected free from the myocardium. The vessel wall was initially mechanically tested biaxially (inflation and axial extension) as a whole (intact wall) and subsequently as intima-media or adventitial layer. A Fung-type exponential strain energy function was used to curve fit the experimental data for the intact wall and individual layers for the right coronary artery and left anterior descending artery. Two methods were used for the determination of material constants, including the Marquardt-Levenberg nonlinear least squares method and the genetic algorithm method. Our results show that there were no statistically significant differences in the material constants obtained from the two methods and that either set of elastic constants results in good fit of the data. Furthermore, at an in vivo value of axial stretch ratio, we find that the stiffness is as follows: intima-media > intact > adventitia. These results underscore the composite nature of coronary arteries with different material properties in each layer. The present results are necessary for analysis of coronary artery mechanics and to provide a fundamental understanding of vessel physiology.