On the effects of residual stress in microindentation tests of soft tissue structures

Microindentation methods are commonly used to determine material properties of soft tissues at the cell or even sub-cellular level. In determining properties from force-displacement (FD) data, it is often assumed that the tissue is initially a stress-free, homogeneous, linear elastic half-space. Residual stress, however, can strongly influence such results. In this paper, we present a new microindentation method for determining both elastic properties and residual stress in soft tissues that, to a first approximation, can be regarded as a pre-stressed layer embedded in or adhered to an underlying relatively soft, elastic foundation. The effects of residual stress are shown using two linear elastic models that approximate specific biological structures. The first model is an axially loaded beam on a relatively soft, elastic foundation (i.e., stress-fiber embedded in cytoplasm), while the second is a radially loaded plate on a foundation (e.g., cell membrane or epithelium). To illustrate our method, we use a nonlinear finite element (FE) model and experimental FD and surface contour data to find elastic properties and residual stress in the early embryonic chick heart, which, in the region near the indenter tip, is approximated as an isotropic circular plate under tension on a foundation. It is shown that the deformation of the surface in a microindentation test can be used along with FD data to estimate material properties, as well as residual stress, in soft tissue structures that can be regarded as a plate under tension on an elastic foundation. This method may not be as useful, however, for structures that behave as a beam on a foundation.     

Robust strategies for automated AFM force curve analysis-II: adhesion-influenced indentation of soft, elastic materials

In the first of this two-part discourse on the extraction of elastic properties from atomic force microscopy (AFM) data, a scheme for automating the analysis of force-distance curves was introduced and experimentally validated for the Hertzian (i.e., linearly elastic and noninteractive probe-sample pairs) indentation of soft, inhomogeneous materials. In the presence of probe-sample adhesive interactions, which are common especially during retraction of the rigid tip from soft materials, the Hertzian models are no longer adequate. A number of theories (e.g., Johnson-Kendall-Roberts and Derjaguin-Muller-Toporov), covering the full range of sample compliance relative to adhesive force and tip radius, are available for analysis of such data. We incorporated Pietrement and Troyon's approximation (2000, "General Equations Describing Elastic Indentation Depth and Normal Contact Stiffness Versus Load," J. Colloid Interface Sci., 226(1), pp. 166-171) of the Maugis-Dugdale model into the automated procedure. The scheme developed for the processing of Hertzian data was extended to allow for adhesive contact by applying the Pietrement-Troyon equation. Retraction force-displacement data from the indentation of polyvinyl alcohol gels were processed using the customized software. Many of the retraction curves exhibited strong adhesive interactions that were absent in extension. We compared the values of Young's modulus extracted from the retraction data to the values obtained from the extension data and from macroscopic uniaxial compression tests. Application of adhesive contact models and the automated scheme to the retraction curves yielded average values of Young's modulus close to those obtained with Hertzian models for the extension curves. The Pietrement-Troyon equation provided a good fit to the data as indicated by small values of the mean-square error. The Maugis-Dugdale theory is capable of accurately modeling adhesive contact between a rigid spherical indenter and a soft, elastic sample. Pietrement and Troyon's empirical equation greatly simplifies the theory and renders it compatible with the general automation strategies that we developed for Hertzian analysis. Our comprehensive algorithm for automated extraction of Young's moduli from AFM indentation data has been expanded to recognize the presence of either adhesive or Hertzian behavior and apply the appropriate contact model.     

Assessing Micromechanical Properties of Cells with Atomic Force Microscopy: Importance of the Contact Point

Mechanical properties are obtainable from atomic force microscopy (AFM) indentation force–depth curves, which are calculated from relationships between tip deflection and cantilever position, i.e. deflection curves. Indentation depth is the difference between tip deflections on a rigid and a soft material for the same amount of cantilever advancement, after contact is made. Since the contact point cannot be unequivocally identified from experimental data, there is some uncertainty in estimating material properties. Using simulations, this study examines some important issues related to the influence of contact point identification on estimated material properties. Simulations for linear materials using a typical stiffness for an AFM cantilever demonstrate that certain portions of the post-contact region of deflection curves for soft and very stiff materials can be approximated by quadratic and linear functions, respectively. Based on these findings, we first develop and verify an objective, automatic method to identify the contact point for materials with linear properties. We then assess the effect of misidentifying the contact point, with and without noise. If the contact point is missed by <50 nm, material properties for small indentations are erroneous but the error decreases asymptotically beyond 200 nm of indentation and the correct estimate of material stiffness is obtained. If the contact point is missed by >100 nm, however, the true material properties cannot be estimated accurately. Noise adds to uncertainty in material properties at small indentations but the combined effect of missing the contact point and noise is dominated by the former. Even though the algorithm was developed for linear materials, it is also suitable for certain nonlinear materials making it more generally applicable.     

A mechanics model of the compression of cells with finite initial contact area

The effect of intercellular bonding on the stress-strain behavior of soft plant tissue is considered. In our mechanical model, a conglomerate of identical cells is arranged in a regular array. Each cell is pressurized and bonded across flat contact areas with adjacent cells in the direction of the applied load. The cell wall is a finitely-deformed mechanical membrane bounding an incompressible fluid (the cytoplasm). A nonlinear elastic constitutive law is presented that describes data for apple parenchyma. Results show that intercellular bonding has a strong effect on the macroscopic properties of the whole tissue. A larger intercellular contact area increases tissue stiffness and magnifies the effect of initial turgor pressure on tissue stiffness.     

Frictional contact mechanics methods for soft materials: application to tracking breast cancers

Mammography is currently the most widely used screening and diagnostic tool for breast cancer. Because X-ray images are 2D projections of a 3D object, it is not trivial to localise features identified in mammogram pairs within the breast volume. Furthermore, mammograms represent highly deformed configurations of the breast under compression, thus the tumour localisation process relies on the clinician's experience. Biomechanical models of the breast undergoing mammographic compressions have been developed to overcome this limitation. In this study, we present the development of a modelling framework that implements Coulomb's frictional law with a finite element analysis using a C(1)-continuous Hermite mesh. We compared two methods of this contact mechanics implementation: the penalty method, and the augmented Lagrangian method, the latter of which is more accurate but computationally more expensive compared to the former. Simulation results were compared with experimental data from a soft silicon gel phantom in order to evaluate the modelling accuracy of each method. Both methods resulted in surface-deformation root-mean-square errors of less than 2mm, whilst the maximum internal marker prediction error was less than 3mm when simulating two mammographic-like compressions. Simulation results were confirmed using the augmented Lagrangian method, which provided similar accuracy. We conclude that contact mechanics on soft elastic materials using the penalty method with an appropriate choice of the penalty parameters provides sufficient accuracy (with contact constraints suitably enforced), and may thus be useful for tracking breast tumours between clinical images.     

Dehomogenized Elastic Properties of Heterogeneous Layered Materials in AFM Indentation Experiments

Atomic force microscopy (AFM) is used to study mechanical properties of biological materials at submicron length scales. However, such samples are often structurally heterogeneous even at the local level, with different regions having distinct mechanical properties. Physical or chemical disruption can isolate individual structural elements but may alter the properties being measured. Therefore, to determine the micromechanical properties of intact heterogeneous multilayered samples indented by AFM, we propose the Hybrid Eshelby Decomposition (HED) analysis, which combines a modified homogenization theory and finite element modeling to extract layer-specific elastic moduli of composite structures from single indentations, utilizing knowledge of the component distribution to achieve solution uniqueness. Using finite element model-simulated indentation of layered samples with micron-scale thickness dimensions, biologically relevant elastic properties for incompressible soft tissues, and layer-specific heterogeneity of an order of magnitude or less, HED analysis recovered the prescribed modulus values typically within 10% error. Experimental validation using bilayer spin-coated polydimethylsiloxane samples also yielded self-consistent layer-specific modulus values whether arranged as stiff layer on soft substrate or soft layer on stiff substrate. We further examined a biophysical application by characterizing layer-specific microelastic properties of full-thickness mouse aortic wall tissue, demonstrating that the HED-extracted modulus of the tunica media was more than fivefold stiffer than the intima and not significantly different from direct indentation of exposed media tissue. Our results show that the elastic properties of surface and subsurface layers of microscale synthetic and biological samples can be simultaneously extracted from the composite material response to AFM indentation. HED analysis offers a robust approach to studying regional micromechanics of heterogeneous multilayered samples without destructively separating individual components before testing.     

Ultrasonic model and system for measurement of corneal biomechanical properties and validation on phantoms

Non-invasive measurement of biomechanical properties of corneas may provide important information for ocular disease management and therapeutic procedures. An ultrasonic non-destructive evaluation method with a wave propagation model was developed to determine corneal biomechanical properties in vivo. In this study, we tested the feasibility of the approach in differentiating the mechanical properties of soft contact lenses as corneal phantoms. Three material types of soft contact lenses (six samples in each group) were measured using a broadband ultrasound transducer. The ultrasonic reflections from the contact lenses were recorded by a 500MHz/8-bit digitizer, and displayed and processed by a PC. A reference signal was recorded to compute the normalized power spectra using Fast Fourier Transformation. An inverse algorithm based on least-squares minimization was used to reconstruct three parameters of the contact lenses: density, thickness, and elastic constants lambda+2micro. The thickness of each sample was verified using an electronic thickness gauge, and the averaged density for each type of lenses was verified using Archimedes' principle and manufacturer's report. Our results demonstrated that the ultrasonic system was able to differentiate the elastic properties of the three types of the soft contact lenses with statistical significance (P-value<0.001). The reconstructed thicknesses and densities agreed well with the independent measurements. Our studies on corneal phantoms indicated that the ultrasonic system was sensitive and accurate in measuring the material properties of cornea-like structures. It is important to optimize the system for in vivo measurements.     

A Mathematical Formulation for 3D Quasi-Static Multibody Models of Diarthrodial Joints

This study describes a genera] set of equations for quasi-static analysis of three-dimensional multibody systems, with a particular emphasis on modeling of diarthrodial joints. The model includes articular contact, muscle forces, tendons and tendon pulleys, ligaments, and the wrapping of soft tissue structures around bone and cartilage surfaces. The general set of equations governing this problem are derived using a consistent notation for all types of links, which can be converted conveniently into efficient computer codes. The computational efficiency of the model is enhanced by the use of analytical Jacobians, particularly in the analysis of articular surface contact and wrapping of soft tissue structures around bone and cartilage surfaces. The usefulness of the multibody model is demonstrated by modeling the patellofemoral joint of six cadaver knees, using cadaver-specific data for the articular surface and bone geometries, as well as tendon and ligament insertions and muscle lines of actions. Good accuracy was observed when comparing the model patellar kinematic predictions to experimental data (mean ± stand, dev. error in translation: 0.63 ± 1.19 mm, 0.10 ± 0.71 mm, -0.29 ± 0.84 mm along medial, proximal, and anterior directions, respectively; in rotation: -1.41 ± 1.71°, 0.27 ±2.38°, -1.13 ± 1.83° in flexion, tilt and rotation, respectively). The accuracy which can be achieved with this type of model, and the computational efficiency of the algorithm employed in this study may serve in many applications such as computer-aided surgical planning, and real-time computer-assisted surgery in the operating room.     

A Mathematical Formulation for 3D Quasi-Static Multibody Models of Diarthrodial Joints

This study describes a general set of equations for quasi-static analysis of three-dimensional multibody systems, with a particular emphasis on modeling of diarthrodial joints. The model includes articular contact, muscle forces, tendons and tendon pulleys, ligaments, and the wrapping of soft tissue structures around bone and cartilage surfaces. The general set of equations governing this problem are derived using a consistent notation for all types of links, which can be converted conveniently into efficient computer codes. The computational efficiency of the model is enhanced by the use of analytical Jacobians, particularly in the analysis of articular surface contact and wrapping of soft tissue structures around bone and cartilage surfaces. The usefulness of the multibody model is demonstrated by modeling the patellofemoral joint of six cadaver knees, using cadaver-specific data for the articular surface and bone geometries, as well as tendon and ligament insertions and muscle lines of actions. Good accuracy was observed when comparing the model patellar kinematic predictions to experimental data (mean +/- stand. dev. error in translation: 0.63 +/- 1.19 mm, 0.10 +/- 0.71 mm, -0.29 +/- 0.84 mm along medial, proximal, and anterior directions, respectively; in rotation: -1.41 +/- 1.71 degrees, 0.27 +/- 2.38 degrees, -1.13 +/- 1.83 degrees in flexion, tilt and rotation, respectively). The accuracy which can be achieved with this type of model, and the computational efficiency of the algorithm employed in this study may serve in many applications such as computer-aided surgical planning, and real-time computer-assisted surgery in the operating room.     

An analytical model of joint contact

The stress distribution in the region of contact between a layered elastic sphere and a layered elastic cavity is determined using an analytical model to stimulate contact of articulating joints. The purpose is to use the solution to analyze the effects of cartilage thickness and stiffness, bone stiffness and joint curvature on the resulting stress field, and investigate the possibility of cracking of the material due to tensile and shear stresses. Vertical cracking of cartilage as well as horizontal splitting at the cartilage-calcified cartilage interface has been observed in osteoarthritic joints. The current results indicate that for a given system (material properties mu and nu constant), the stress distribution is a function of the ratio of contact radius to layer thickness (a/h), and while tensile stresses are seen to occur only when a/h is small, tensile strain is observed for all a/h values. Significant shear stresses are observed at the cartilage-bone interface. Softening of cartilage results in an increase in a/h, and a decrease in maximum normal stress. Cartilage thinning increases a/h and the maximum contact stress, while thickening has the opposite effect. A reduction in the indenting radius reduces a/h and increases the maximum normal stress. Bone softening is seen to have negligible effect on the resulting contact parameters and stress distribution.     

How the stiffness of meniscal attachments and meniscal material properties affect tibio-femoral contact pressure computed using a validated finite element model of the human knee joint

In an effort to prevent degeneration of articular cartilage associated with meniscectomies, both meniscal allografts and synthetic replacements are subjects of current interest and investigation. The objectives of the current study were to (1) determine whether a transversely isotropic, linearly elastic, homogeneous material model of the meniscal tissue is necessary to achieve a normal contact pressure distribution on the tibial plateau, (2) determine which material and boundary condition (attachments) parameters affect the contact pressure distribution most strongly, and (3) set tolerances on these parameters to restore the contact pressure distribution to within a specified error. To satisfy these objectives, a finite element model of the tibio-femoral joint of a human cadaveric knee (including both menisci) was used to study the contact pressure distribution on the tibial plateau. To validate the model, the contact pressure distribution on the tibial plateau was measured experimentally in the same knee used to create the model. Within physiologically reasonable bounds on five material parameters and four attachment parameters associated with a meniscal replacement, an optimization was performed under 1200 N of compressive load on the set of nine parameters to minimize the difference between the experimental and model results. The error between the experimental and model contact variables was minimized to 5.4%. The contact pressure distribution of the tibial plateau was sensitive to the circumferential modulus, axial/radial modulus, and horn stiffness, but relatively insensitive to the remaining six parameters. Consequently, both the circumferential and axial/radial moduli are important determinants of the contact pressure distribution, and hence should be matched in the design and/or selection of meniscal replacements. In addition, during surgical implantation of a meniscal replacement, the horns should be attached with high stiffness bone plugs, and the attachments of the transverse ligament and deep medial collateral ligament should be restored to minimize changes in the contact pressure distribution, and thereby possibly prevent the degradation of articular cartilage.     

High-performance cation-exchange chromatography of proteins

Approximately one-third of all proteins reported in the literature have a pI sufficiently high to be resolved by cation-exchange chromatography. This paper reports the preparation and use of new high-performance polymeric-bonded-phase cation-exchange columns. Starting from a very stable, covalently bonded polyamide coating on microparticulate silica, simple derivatization produces a versatile cation-exchange material useful for separations traditionally performed on classical carboxymethylated soft gel supports. Column behavior was monitored using chymotrypsinogen, cytochrome c, and lysozyme as standards. The polymeric bonded phase was stable to pH 2.5 and exhibits enhanced selectivity for proteins due to a slight hydrophobic character of the matrix. Several separations of biological interest that demonstrate the utility of these small cation-exchange columns for modern biochemical separations are shown.     

Mechanics of inactive swelling and bursting of porate pollen grains

The structure of pollen grains, which is typically characterized by soft apertures in an otherwise stiff exine shell, guides their response to changes in the humidity of the environment. These changes can lead to desiccation of the grain and its infolding but also to excessive swelling of the grain and even its bursting. Here we use an elastic model to explore the mechanics of pollen grain swelling and the role of soft, circular apertures (pores) in this process. Small, circular apertures typically occur in airborne and allergenic pollen grains so that the bursting of such grains is important in the context of human health. We identify and quantify a mechanical weakness of the pores, which are prone to rapid inflation when the grain swells to a critical extent. The inflation occurs as a sudden transition and may induce bursting of the grain and release of its content. This process crucially depends on the size of the pores and their softness. Our results provide insight into the inactive part of the mechanical response of pollen grains to hydration when they land on a stigma as well as bursting of airborne pollen grains during changes in air humidity.     

Breakdown of Methods for Phasing and Imputation in the Presence of Double Genotype Sharing

In genome-wide association studies, results have been improved through imputation of a denser marker set based on reference haplotypes and phasing of the genotype data. To better handle very large sets of reference haplotypes, pre-phasing with only study individuals has been suggested. We present a possible problem which is aggravated when pre-phasing strategies are used, and suggest a modification avoiding the resulting issues with application to the MaCH tool, although the underlying problem is not specific to that tool.We evaluate the effectiveness of our remedy to a subset of Hapmap data, comparing the original version of MaCH and our modified approach. Improvements are demonstrated on the original data (phase switch error rate decreasing by 10%), but the differences are more pronounced in cases where the data is augmented to represent the presence of closely related individuals, especially when siblings are present (30% reduction in switch error rate in the presence of children, 47% reduction in the presence of siblings).The main conclusion of this investigation is that existing statistical methods for phasing and imputation of unrelated individuals might give results of sub-par quality if a subset of study individuals nonetheless are related. As the populations collected for general genome-wide association studies grow in size, including relatives might become more common. If a general GWAS framework for unrelated individuals would be employed on datasets with some related individuals, such as including familial data or material from domesticated animals, caution should also be taken regarding the quality of haplotypes.Our modification to MaCH is available on request and straightforward to implement. We hope that this mode, if found to be of use, could be integrated as an option in future standard distributions of MaCH.     

Nonlinear viscoelastic, thermodynamically consistent, models for biological soft tissue

The mechanical behavior of most biological soft tissue is nonlinear viscoelastic rather than elastic. Many of the models previously proposed for soft tissue involve ad hoc systems of springs and dashpots or require measurement of time-dependent constitutive coefficient functions. The model proposed here is a system of evolution differential equations, which are determined by the long-term behavior of the material as represented by an energy function of the type used for elasticity. The necessary empirical data is time independent and therefore easier to obtain. These evolution equations, which represent non-equilibrium, transient responses such as creep, stress relaxation, or variable loading, are derived from a maximum energy dissipation principle, which supplements the second law of thermodynamics. The evolution model can represent both creep and stress relaxation, depending on the choice of control variables, because of the assumption that a unique long-term manifold exists for both processes. It succeeds, with one set of material constants, in reproducing the loading-unloading hysteresis for soft tissue. The models are thermodynamically consistent so that, given data, they may be extended to the temperature-dependent behavior of biological tissue, such as the change in temperature during uniaxial loading. The Holzapfel et al. three-dimensional two-layer elastic model for healthy artery tissue is shown to generate evolution equations by this construction for biaxial loading of a flat specimen. A simplified version of the Shah-Humphrey model for the elastodynamical behavior of a saccular aneurysm is extended to viscoelastic behavior.     

Mathematical correlation between biomaterial and cellular parameters--critical reflection of statistics

For mathematical modelling of the biomaterial-cell contact, it is necessary to find both parameters characterizing physical and chemical properties of the material surface and also such describing the reaction of the adhering cells. Only those material and cell parameters that correlate with each other are applicable to model this contact mathematically. Only few papers are dealing with this special problem. The aim of this paper is to present results of physical/chemical and biological investigations made on differently modified rough titanium implant surfaces in order to find out only the correlating parameters. Furthermore we discuss several ways to apply statistical methods to the correlation problem. Only few ones of all investigated parameters both on material and on cellular side were applicable for correlation. For example we found in our studies that fractal structure parameter topothesy has influence on the spreading behaviour of the osteoblastic cells. However the value of the correlation coefficient and its statistical significance heavily depend on the method of averaging the available data. Especially the biological data (spreading area) were afflicted with relatively high error up to 30%. Averaging of this data masks the true facts. That is why the correlation coefficient considerably decreases if the biological parameters are not averaged. On the other hand, the statistical reliability increases due to the higher number of investigated cases. Critical error discussion is necessary in statistical correlation between material and biological parameters. Often the results are heavily influenced by the statistical handling of data, especially if only few data are available. May be that new unconventional methods like bootstrap method can show a way out of this dilemma.     

Ordinary least squares regression is indicated for studies of allometry

When it comes to fitting simple allometric slopes through measurement data, evolutionary biologists have been torn between regression methods. On the one hand, there is the ordinary least squares (OLS) regression, which is commonly used across many disciplines of biology to fit lines through data, but which has a reputation for underestimating slopes when measurement error is present. On the other hand, there is the reduced major axis (RMA) regression, which is often recommended as a substitute for OLS regression in studies of allometry, but which has several weaknesses of its own. Here, we review statistical theory as it applies to evolutionary biology and studies of allometry. We point out that the concerns that arise from measurement error for OLS regression are small and straightforward to deal with, whereas RMA has several key properties that make it unfit for use in the field of allometry. The recommended approach for researchers interested in allometry is to use OLS regression on measurements taken with low (but realistically achievable) measurement error. If measurement error is unavoidable and relatively large, it is preferable to correct for slope attenuation rather than to turn to RMA regression, or to take the expected amount of attenuation into account when interpreting the data.     

The effect of Glisson's capsule on the superficial elasticity measurements of the liver

In the past decade, novel tools for surgical planning and disease diagnosis have been developed to detect the liver's mechanical properties. Some tools utilize superficial indentation type measurements to determine the elasticity of the liver parenchyma and to assume material homogeneity. In fact, the liver is a soft tissue covered with a connective sheathing that is called Glisson's capsule. This article quantifies the effect of this capsule on the measured or "effective" elastic modulus obtained by indentation with a spherical geometry. Two sets of parametric computational studies in which the Glisson capsule thickness and elasticity were varied, demonstrated the relationship between the measured elastic response and the underlying parenchymal elastic response. Previously reported in vivo indentation data on the human liver were utilized to determine the elasticity of its parenchyma. The results indicated a linear relationship between the effective (measured) elastic response and the underlying parenchyma for the Glisson capsule thicknesses considered. When previously published human liver indentation data were analyzed, the measured elastic modulus was nearly 6.9% greater than the underlying parenchyma elastic modulus. Although the analyzed data were obtained from a single liver and yet displayed a significant variation, the Glisson capsule elasticity induced a significant but systematic error as well. The Glisson capsule thickness error was negligible for capsule parameters associated with a normal liver. Based on this work, an emphasis on the Glisson capsule's contribution to the mechanical response of the liver would enhance the clinical potential of indentation-based novel tools for liver care.     

An improved parameter estimation and comparison for soft tissue constitutive models containing an exponential function

Motivated by the well-known result that stiffness of soft tissue is proportional to the stress, many of the constitutive laws for soft tissues contain an exponential function. In this work, we analyze properties of the exponential function and how it affects the estimation and comparison of elastic parameters for soft tissues. In particular, we find that as a consequence of the exponential function there are lines of high covariance in the elastic parameter space. As a result, one can have widely varying mechanical parameters defining the tissue stiffness but similar effective stress–strain responses. Drawing from elementary algebra, we propose simple changes in the norm and the parameter space, which significantly improve the convergence of parameter estimation and robustness in the presence of noise. More importantly, we demonstrate that these changes improve the conditioning of the problem and provide a more robust solution in the case of heterogeneous material by reducing the chances of getting trapped in a local minima. Based upon the new insight, we also propose a transformed parameter space which will allow for rational parameter comparison and avoid misleading conclusions regarding soft tissue mechanics.