On the compressibility and poroelasticity of human and murine skin

Poisson function and volume ratio of soft anisotropic materials under large deformations

Accurate transverse deformation measurements are required for the estimation of the Poisson function and volume ratio. In this study, pure silicone and soft composite specimens were subjected to uniaxial tension, and the digital image correlation method was used to measure longitudinal and in- and out-of-plane transverse stretches. To minimize the effects of measurement errors on parameter estimation, the measured transverse stretches were defined in terms of the longitudinal stretch using a new formulation based on Poisson's ratios and two stretch-dependent parameters. From this formulation, Poisson functions and volume ratio for soft materials under large deformations were obtained. The results showed that pure silicone can be considered isotropic and nearly incompressible under large deformations, as expected. In contrast, Poisson's ratio of silicone reinforced with extensible fabric can exceed classical bounds, including negative value (auxetic behavior). The incompressibility assumption can be employed for describing the stress-stretch curve of pure silicone, while volume ratios are required for soft composites. Data of human skin, aortic wall, and annulus fibrosus from the literature were selected and analyzed. Except for the aortic wall, which can be considered nearly incompressible, the studied soft tissues must be regarded as compressible. All tissues presented anisotropic behavior.

Variations in fluid chemical potential induce fibroblast mechano-response in 3D hydrogels

Mechanical deformation of skin creates variations in fluid chemical potential, leading to local changes in hydrostatic and osmotic pressure, whose effects on mechanobiology remain poorly understood. To study these effects, we investigate the specific influences of hydrostatic and osmotic pressure on primary human dermal fibroblasts in three-dimensional hydrogel culture models. Cyclic hydrostatic pressure and hyperosmotic stress enhanced the percentage of cells expressing the proliferation marker Ki67 in both collagen and PEG-based hydrogels. Osmotic pressure also activated the p38 MAPK stress response pathway and increased the expression of the osmoresponsive genes PRSS35 and NFAT5. When cells were cultured in two-dimension (2D), no change in proliferation was observed with either hydrostatic or osmotic pressure. Furthermore, basal, and osmotic pressure-induced expression of osmoresponsive genes differed in 2D culture versus 3D hydrogels, highlighting the role of dimensionality in skin cell mechanotransduction and stressing the importance of 3D tissue-like models that better replicate in vivo conditions. Overall, these results indicate that fluid chemical potential changes affect dermal fibroblast mechanobiology, which has implications for skin function and for tissue regeneration strategies.

Quantification of age-related changes in the structure and mechanical function of skin with multiscale imaging

The mechanical properties of skin change during aging but the relationships between structure and mechanical function remain poorly understood. Previous work has shown that young skin exhibits a substantial decrease in tissue volume, a large macro-scale Poisson's ratio, and an increase in micro-scale collagen fiber alignment during mechanical stretch. In this study, label-free multiphoton microscopy was used to quantify how the microstructure and fiber kinematics of aged mouse skin affect its mechanical function. In an unloaded state, aged skin was found to have less collagen alignment and more non-enzymatic collagen fiber crosslinks. Skin samples were then loaded in uniaxial tension and aged skin exhibited a lower mechanical stiffness compared to young skin. Aged tissue also demonstrated less volume reduction and a lower macro-scale Poisson's ratio at 10% uniaxial strain, but not at 20% strain. The magnitude of 3D fiber realignment in the direction of loading was not different between age groups, and the amount of realignment in young and aged skin was less than expected based on theoretical fiber kinematics affine to the local deformation. These findings provide key insights on how the collagen fiber microstructure changes with age, and how those changes affect the mechanical function of skin, findings which may help guide wound healing or anti-aging treatments.

Biomechanical and biochemical changes in murine skin during development and aging

Aging leads to biochemical and biomechanical changes in skin, with biological and functional consequences. Despite extensive literature on skin aging, there is a lack of studies which investigate the maturation of the tissue and connect the microscopic changes in the skin to its macroscopic biomechanical behavior as it evolves over time. The present work addresses this knowledge gap using multiscale characterization of skin in a murine model considering newborn, adult and aged mice. Monotonic uniaxial loading, tension relaxation with change of bath, and loading to failure tests were performed on murine skin samples from different age groups, complemented by inflation experiments and atomic force microscopy indentation measurements. In parallel, skin samples were characterized using histological and biochemical techniques to assess tissue morphology, collagen organization, as well as collagen content and cross-linking. We show that 1-week-old skin differs across nearly all measured parameters from adult skin, showing reduced strain stiffening and tensile strength, a thinner dermis, lower collagen content and altered crosslinking patterns. Surprisingly, adult and aged skin were similar across most biomechanical parameters in the physiologic loading range, while aged skin had lower tensile strength and lower stiffening behavior at large force values. This correlates with altered collagen content and cross-links. Based on a computational model, differences in mechanocoupled stimuli in the skin of the different age groups were calculated, pointing to a potential biological significance of the age-induced biomechanical changes in regulating the local biophysical environment of dermal cells. STATEMENT OF SIGNIFICANCE: Skin microstructure and the emerging mechanical properties change with age, leading to biological, functional and health-related consequences. Despite extensive literature on skin aging, only very limited quantitative data are available on microstructural changes and the corresponding macroscopic biomechanical behavior as they evolve over time. This work provides a wide-range multiscale mechanical characterization of skin of newborn, adult and aged mice, and quantifies microstructural correlations in tissue morphology, collagen content, organization and cross-linking. Remarkably, aged skin retained normal hydration and normal biomechanical function in the physiological loading range but showed significantly reduced properties at super-physiological loading. Our data show that age-related microstructural differences have a profound effect not only on tissue-level properties but also on the cell-level biophysical environment.

Second harmonic generation microscopy, biaxial mechanical tests and fiber dispersion models in human skin biomechanics

Advanced numerical simulations of the mechanical behavior of human skin require thorough calibration of the material's constitutive models based on experimental ex vivo mechanical tests along with images of tissue microstructure for a variety of biomedical applications. In this work, a total of 14 human healthy skin samples and 4 additional scarred skin samples were experimentally analyzed to gain deep insights into the biomechanics of human skin. In particular, second harmonic generation (SHG) microscopy was used to extract detailed images of the distribution of collagen fibers, which were subsequently processed using a three-dimensional Fourier transform-based method recently proposed by the authors to quantify the distribution of fiber orientations. Mechanical tests under both biaxial and uniaxial loading were performed to calibrate the relevant mechanical parameters of two widely used constitutive models of soft fiber-reinforced biological tissues that account for non-symmetrical fiber dispersion. The calibration of the models allowed us to identify correlations between the mechanical parameters of the constitutive models considered. STATEMENT OF SIGNIFICANCE: Constitutive models for soft collagenous tissues can accurately reproduce the complex nonlinear and anisotropic mechanical behavior of skin. However, a comprehensive analysis of both microstructural and mechanical parameters is still missing for human skin. In this study, these parameters are determined by combining biaxial mechanical tests and SHG stacks of collagen fibers on ex vivo healthy human skin samples. The constitutive parameters are provided for two widely used hyperelastic models and enable accurate characterization of skin mechanical behavior for advanced numerical simulations.

Investigation of dynamic responses of skin simulant against fragment impact through experiments and concurrent computational modeling

Perforation of the skin by fragment impact is a key determinant of the severity of an injury and incapacitation during modern asymmetric warfare. Computational models validated against experimental data are thus desired for simulating the responses of a skin simulant against fragment impact. Toward this end, experiments and concurrent computational modeling were used to investigate the dynamic responses of the skin simulant against fragment impact. Fragment simulating projectiles (FSPs) of masses 1.10 g and 2.79 g were considered herein, and the responses of the skin simulant were investigated in terms of the threshold velocity, energy density, peak displacement, and failure mechanisms. The results illustrate numerous salient aspects. The skin simulant failure involved cavity shearing followed by elastic hole enlargement, and these results were sensitive to the strain rate. The best agreement between the simulated and experimental results was achieved when the input stress-strain curves to the simulation were based on the full spectrum of strain rates. When a single stress-strain curve corresponding to a specific strain rate was used as the input, the threshold velocity and peak displacement of the skin simulant were either underpredicted or overpredicted depending on the strain rate considered. The threshold velocity was also sensitive to the input failure strain; here, the best agreement was obtained when the failure strain was based on the theoretical limiting strain. When the FSP materials were changed to plastics, the threshold velocities increased by up to 33%; however, the energy densities and generated stresses exceeded the contusion and laceration thresholds of the skin.

Modeling drug transport and absorption in subcutaneous injection of monoclonal antibodies: Impact of tissue deformation, devices, and physiology

The pharmaceutical industry has experienced a remarkable increase in the use of subcutaneous injection of monoclonal antibodies (mAbs), attributed mainly to its advantages in reducing healthcare-related costs and enhancing patient compliance. Despite this growth, there is a limited understanding of how tissue mechanics, physiological parameters, and different injection devices and techniques influence the transport and absorption of the drug. In this work, we propose a high-fidelity computational model to study drug transport and absorption during and after subcutaneous injection of mAbs. Our numerical model includes large-deformation mechanics, fluid flow, drug transport, and blood and lymphatic uptake. Through this computational framework, we analyze the tissue material responses, plume dynamics, and drug absorption. We analyze different devices, injection techniques, and physiological parameters such as BMI, flow rate, and injection depth. Finally, we compare our numerical results against the experimental data from the literature.

A quadriphasic mechanical model of the human dermis

The present study investigates the multiphasic nature of the mechanical behavior of human dermis. Motivated by experimental observations and by consideration of its composition, a quadriphasic model of the dermis is proposed, distinguishing solid matrix components, interstitial fluid and charged constituents moving within the fluid, i.e., anions and cations. Compression and tensile experiments with and without change of osmolarity of the bath are performed to characterize the chemo-mechanical coupling in the dermis. Model parameters are determined through inverse analysis. The computations predict a dominant role of the permeability in the determination of the temporal evolution of the mechanical response of the tissue. In line with the previous studies on other tissues, the analysis shows that an ideal model based on Donnan's equilibrium overestimates the osmotic pressure in skin for the case of very dilute solutions. The quadriphasic model is applied to predict changes in dermal cell environment and therefore alterations in what is called the "mechanome," associated with skin stretch. The simulations indicate that skin deformation causes a variation in several local variables, including in particular the electric field associated with a deformation-induced non-homogeneous distribution of fixed charges.

Mechanical Models of Collagen Networks for Understanding Changes in the Failure Properties of Aging Skin

Skin undergoes mechanical alterations due to changes in the composition and structure of the collagenous dermis with aging. Previous studies have conflicting findings, with both increased and decreased stiffness reported for aging skin. The underlying structure-function relationships that drive age-related changes are complex and difficult to study individually. One potential contributor to these variations is the accumulation of nonenzymatic crosslinks within collagen fibers, which affect dermal collagen remodeling and mechanical properties. Specifically, these crosslinks make individual fibers stiffer in their plastic loading region and lead to increased fragmentation of the collagenous network. To better understand the influence of these changes, we investigated the impact of nonenzymatic crosslink changes on the dermal microstructure using discrete fiber networks representative of the dermal microstructure. Our findings suggest that stiffening the plastic region of collagen's mechanical response has minimal effects on network-level stiffness and failure stresses. Conversely, simulating fragmentation through a loss of connectivity substantially reduces network stiffness and failure stress, while increasing stretch ratios at failure.

Sustained Physiological Stretch Induces Abdominal Skin Growth in Pregnancy

Supraphysiological stretches are exploited in skin expanders to induce tissue growth for autologous implants. As pregnancy is associated with large levels of sustained stretch, we investigated whether skin growth occurs in pregnancy. Therefore, we combined a mechanical model of skin and the observations from suction experiments on several body locations of five pregnant women at different gestational ages. The measurements show a continuous increase in stiffness, with the largest change observed during the last trimester. A comparison with numerical simulations indicates that the measured increase in skin stiffness is far below the level expected for the corresponding deformation of abdominal skin. A new set of simulations accounting for growth could rationalize all observations. The predicted amount of tissue growth corresponds to approximately 40% area increase before delivery. The results of the simulations also offered the opportunity to investigate the biophysical cues present in abdominal skin along gestation and to compare them with those arising in skin expanders. Alterations of the skin mechanome were quantified, including tissue stiffness, hydrostatic and osmotic pressure of the interstitial fluid, its flow velocity and electrical potential. The comparison between pregnancy and skin expansion highlights similarities as well as differences possibly influencing growth and remodeling.

Poroelastic Characterization and Modeling of Subcutaneous Tissue Under Confined Compression

Subcutaneous tissue mechanics are important for drug delivery. Yet, even though this material is poroelastic, its mechanical characterization has focused on its hyperelastic response. Moreover, advancement in subcutaneous drug delivery requires effective tissue mimics such as hydrogels for which similar gaps of poroelastic data exist. Porcine subcutaneous samples and gelatin hydrogels were tested under confined compression at physiological conditions and strain rates of 0.01%/s in 5% strain steps with 2600 s of stress relaxation between loading steps. Force-time data were used in an inverse finite element approach to obtain material parameters. Tissues and gels were modeled as porous neo-Hookean materials with properties specified via shear modulus, effective solid volume fraction, initial hydraulic permeability, permeability exponent, and normalized viscous relaxation moduli. The constitutive model was implemented into an isogeometric analysis (IGA) framework to study subcutaneous injection. Subcutaneous tissue exhibited an initial spike in stress due to compression of the solid and fluid pressure buildup, with rapid relaxation explained by fluid drainage, and longer time-scale relaxation explained by viscous dissipation. The inferred parameters aligned with the ranges reported in the literature. Hydraulic permeability, the most important parameter for drug delivery, was in the rangek 0 ∈ [ 0.142 , 0.203 ] mm4 /(N s). With these parameters, IGA simulations showed peak stresses due to a 1-mL injection to reach 48.8 kPa at the site of injection, decaying after drug volume disperses into the tissue. The poro-hyper-viscoelastic neo-Hookean model captures the confined compression response of subcutaneous tissue and gelatin hydrogels. IGA implementation enables predictive simulations of drug delivery.

Geometry of adipocyte packing in subcutaneous tissue contributes to nonlinear tissue properties captured through a Gaussian process surrogate model

Subcutaneous tissue mechanical response is governed by the geometry and mechanical properties at the microscale and drives physiological and clinical processes such as drug delivery. Even though adipocyte packing is known to change with age, disease, and from one individual to another, the link between the geometry of the packing and the overall mechanical response of adipose tissue remains poorly understood. Here we create 1200 periodic representative volume elements (RVEs) that sample the possible space of Laguerre packings describing adipose tissue. RVE mechanics are modeled under tri-axial loading. Equilibrium configuration of RVEs is solved by minimizing an energetic potential that includes volume change contributions from adipocyte expansion, and area change contributions from collagen foam stretching. The resulting mechanical response across all RVE samples is interpolated with the aid of a Gaussian process (GP), revealing how the microscale geometry dictates the overall RVE mechanics. For example, increase in adipocyte size and increase in sphericity lead to adipose tissue softening. We showcase the use of the homogenized model in finite element simulations of drug injection by implementing a Blatz-Ko model, informed by the GP, as a custom material in the popular open-source package FEBio. These simulations show how microscale geometry can lead to vastly different injection dynamics even if the constituent parameters are held constant, highlighting the importance of characterizing individual's adipose tissue structure in the development of personalized therapies.

The role of interface geometry and appendages on the mesoscale mechanics of the skin

The skin is the largest organ in the human body and serves various functions, including mechanical protection and mechanosensation. Yet, even though skin's biomechanics are attributed to two main layers-epidermis and dermis-computational models have often treated this tissue as a thin homogeneous material or, when considering multiple layers, have ignored the most prominent heterogeneities of skin seen at the mesoscale. Here, we create finite element models of representative volume elements (RVEs) of skin, including the three-dimensional variation of the interface between the epidermis and dermis as well as considering the presence of hair follicles. The sinusoidal interface, which approximates the anatomical features known as Rete ridges, does not affect the homogenized mechanical response of the RVE but contributes to stress concentration, particularly at the valleys of the Rete ridges. The stress profile is three-dimensional due to the skin's anisotropy, leading to high-stress bands connecting the valleys of the Rete ridges through one type of saddle point. The peaks of the Rete ridges and the other class of saddle points of the sinusoidal surface form a second set of low-stress bands under equi-biaxial loading. Another prominent feature of the heterogeneous stress pattern is a switch in the stress jump across the interface, which becomes lower with respect to the flat interface at increasing deformations. These features are seen in both tension and shear loading. The RVE with the hair follicle showed strains concentrating at the epidermis adjacent to the hair follicle, the epithelial tissue surrounding the hair right below the epidermis, and the bulb or base region of the hair follicle. The regions of strain concentration near the hair follicle in equi-biaxial and shear loading align with the presence of distinct mechanoreceptors in the skin, except for the bulb or base region. This study highlights the importance of skin heterogeneities, particularly its potential mechanophysiological role in the sense of touch and the prevention of skin delamination.

Cracks in tensile-contracting and tensile-dilating poroelastic materials

Fibrous gels such as cartilage, blood clots, and carbon-nanotube-based sponges with absorbed oils suffer a reduction in volume by the expulsion of liquid under uniaxial tension, and this directly affects crack-tip fields and energy release rates. A continuum model is formulated for isotropic fibrous gels that exhibit a range of behaviors from volume increasing to volume decreasing in uniaxial tension by changing the ratio of two material parameters. The motion of liquid in the pores of such gels is modeled using poroelasticity. The direction of liquid fluxes around cracks is shown to depend on whether the gel locally increases or decreases in volume. The energy release rate for cracks is computed using a surface-independent integral and it is shown to have two contributions - one from the stresses in the solid network, and another from the flow of liquid. The contribution to the integral from liquid permeation tends to be negative when the gel exhibits volume decrease, which effectively is a crack shielding mechanism.

Radial matrix constraint influences tissue contraction and promotes maturation of bi-layered skin equivalents

Human skin equivalents (HSEs) serve as important tools for mechanistic studies with human skin cells, drug discovery, pre-clinical applications in the field of tissue engineering and for skin transplantation on skin defects. Besides the cellular and extracellular matrix (ECM) components used for HSEs, physical constraints applied on the scaffold during HSEs maturation influence tissue organization, functionality, and homogeneity. In this study, we introduce a 3D-printed culture insert that exposes bi-layered HSEs to a static radial constraint through matrix adhesion. We examine the effect of various diameters of the ring-shaped culture insert on the HSE's characteristics and compare them to state-of-the-art unconstrained and planar constrained HSEs. We show that radial matrix constraint of HSEs regulates tissue contraction, promotes fibroblast and matrix organization that is similar to human skin in vivo and improves keratinocyte differentiation, epidermal stratification, and basement membrane formation depending on the culture insert diameter. Together, these data demonstrate that the degree of HSE's contraction is an important design consideration in skin tissue engineering. Therefore, this study can help to mimic various in vivo skin conditions and to increase the control of relevant tissue properties.

Computational modeling of the effect of skin pinch and stretch on subcutaneous injection of monoclonal antibodies using autoinjector devices

Subcutaneous injection of monoclonal antibodies (mAbs) has experienced unprecedented growth in the pharmaceutical industry due to its benefits in patient compliance and cost-effectiveness. However, the impact of different injection techniques and autoinjector devices on the drug's transport and uptake is poorly understood. Here, we develop a biphasic large-deformation chemomechanical model that accounts for the components of the extracellular matrix that govern solid deformation and fluid flow within the subcutaneous tissue: interstitial fluid, collagen fibers and negatively charged proteoglycan aggregates. We use this model to build a high-fidelity representation of a virtual patient performing a subcutaneous injection of mAbs. We analyze the impact of the pinch and stretch methods on the injection dynamics and the use of different handheld autoinjector devices. The results suggest that autoinjector base plates with a larger device-skin contact area cause significantly lower tissue mechanical stress, fluid pressure and fluid velocity during the injection process. Our simulations indicate that the stretch technique presents a higher risk of intramuscular injection for autoinjectors with a relatively long needle insertion depth.

Multiscale mechanical analysis of the elastic modulus of skin

The mechanical properties of the skin determine tissue function and regulate dermal cell behavior. Yet measuring these properties remains challenging, as evidenced by the large range of elastic moduli reported in the literature-from below one kPa to hundreds of MPa. Here, we reconcile these disparate results by dedicated experiments at both tissue and cellular length scales and by computational models considering the multiscale and multiphasic tissue structure. At the macroscopic tissue length scale, the collective behavior of the collagen fiber network under tension provides functional tissue stiffness, and its properties determine the corresponding elastic modulus (100-200 kPa). The compliant microscale environment (0.1-10 kPa), probed by atomic force microscopy, arises from the ground matrix without engaging the collagen fiber network. Our analysis indicates that indentation-based elasticity measurements, although probing tissue properties at the cell-relevant length scale, do not assess the deformation mechanisms activated by dermal cells when exerting traction forces on the extracellular matrix. Using dermal-equivalent collagen hydrogels, we demonstrate that indentation measurements of tissue stiffness do not correlate with the behavior of embedded dermal fibroblasts. These results provide a deeper understanding of tissue mechanics across length scales with important implications for skin mechanobiology and tissue engineering. STATEMENT OF SIGNIFICANCE: Measuring the mechanical properties of the skin is essential for understanding dermal cell mechanobiology and designing tissue-engineered skin substitutes. However, previous results reported for the elastic modulus of skin vary by six orders of magnitude. We show that two distinct deformation mechanisms, related to the tension-compression nonlinearity of the collagen fiber network, can explain the large variations in elastic moduli. Furthermore, we show that microscale indentation, which is frequently used to assess the stiffness perceived by cells, fails to engage the fiber network, and therefore cannot predict the behavior of dermal fibroblasts in stiffness-tunable fibrous hydrogels. This has important implications for how to measure and interpret the mechanical properties of soft tissues across length scales.

A quantitative comparison of devices for in vivo biomechanical characterization of human skin

Non-invasive skin characterization devices are emerging as a valuable tool in clinical skin research. In recent years, the range of available experimental techniques and methods used to determine the biomechanical properties of skin has increased considerably. Although a substantial amount of work has been devoted to assessing the working principle of macroscopic skin characterization devices individually, a rationalization and comparison between them is still lacking. This motivated the present study, which aimed to characterize and compare three commonly used working principles: suction, dynamic shear loading, and indentation. A synthetic model system with tunable mechanical properties was used to assess the three devices, and the results rationalized based on corresponding finite element models. In vivo measurements were performed on healthy volunteers to investigate the capability of differentiating the biomechanical properties of skin at different body locations, and to assess the intra- and inter-rater reliability of each device. The present comparative analysis indicates that the analyzed functional principles perceive the stiffness of human skin differently, with relevant implications for the interpretation of the respective measurement results.

Control of hydrostatic pressure and osmotic stress in 3D cell culture for mechanobiological studies

Hydrostatic pressure (HP) and osmotic stress (OS) play an important role in various biological processes, such as cell proliferation and differentiation. In contrast to canonical mechanical signals transmitted through the anchoring points of the cells with the extracellular matrix, the physical and molecular mechanisms that transduce HP and OS into cellular functions remain elusive. Three-dimensional cell cultures show great promise to replicate physiologically relevant signals in well-defined host bioreactors with the goal of shedding light on hidden aspects of the mechanobiology of HP and OS. This review starts by introducing prevalent mechanisms for the generation of HP and OS signals in biological tissues that are subject to pathophysiological mechanical loading. We then revisit various mechanisms in the mechanotransduction of HP and OS, and describe the current state of the art in bioreactors and biomaterials for the control of the corresponding physical signals.

Multi-Fidelity Gaussian Process Surrogate Modeling of Pediatric Tissue Expansion

Growth of skin in response to stretch is the basis for tissue expansion (TE), a procedure to gain new skin area for reconstruction of large defects. Unfortunately, complications and suboptimal outcomes persist because TE is planned and executed based on physician's experience and trial and error instead of predictive quantitative tools. Recently, we calibrated computational models of TE to a porcine animal model of tissue expansion, showing that skin growth is proportional to stretch with a characteristic time constant. Here, we use our calibrated model to predict skin growth in cases of pediatric reconstruction. Available from the clinical setting are the expander shapes and inflation protocols. We create low fidelity semi-analytical models and finite element models for each of the clinical cases. To account for uncertainty in the response expected from translating the models from the animal experiments to the pediatric population, we create multifidelity Gaussian process surrogates to propagate uncertainty in the mechanical properties and the biological response. Predictions with uncertainty for the clinical setting are essential to bridge our knowledge from the large animal experiments to guide and improve the treatment of pediatric patients. Future calibration of the model with patient-specific data-such as estimation of mechanical properties and area growth in the operating room-will change the standard for planning and execution of TE protocols.

Activation of Nrf2 in fibroblasts promotes a skin aging phenotype via an Nrf2-miRNA-collagen axis

Aging is associated with progressive skin fragility and a tendency to tear, which can lead to severe clinical complications. The transcription factor NRF2 is a key regulator of the cellular antioxidant response, and pharmacological NRF2 activation is a promising strategy for the prevention of age-related diseases. Using a combination of molecular and cellular biology, histology, imaging and biomechanical studies we show, however, that constitutive genetic activation of Nrf2 in fibroblasts of mice suppresses collagen and elastin expression, resulting in reduced skin strength as seen in aged mice. Mechanistically, the "aging matrisome" results in part from direct Nrf2-mediated overexpression of a network of microRNAs that target mRNAs of major skin collagens and other matrix components. Bioinformatics and functional studies revealed high NRF2 activity in aged human fibroblasts in 3D skin equivalents and human skin biopsies, highlighting the translational relevance of the functional mouse data. Together, these results identify activated NRF2 as a promoter of age-related molecular and biomechanical skin features.

Multiscale Computational Model Predicts Mouse Skin Kinematics Under Tensile Loading

Skin is a complex tissue whose biomechanical properties are generally understood in terms of an incompressible material whose microstructure undergoes affine deformations. A growing number of experiments, however, have demonstrated that skin has a high Poisson's ratio, substantially decreases in volume during uniaxial tensile loading, and demonstrates collagen fiber kinematics that are not affine with local deformation. In order to better understand the mechanical basis for these properties, we constructed multiscale mechanical models (MSM) of mouse skin based on microstructural multiphoton microscopy imaging of the dermal microstructure acquired during mechanical testing. Three models that spanned the cases of highly aligned, moderately aligned, and nearly random fiber networks were examined and compared to the data acquired from uniaxially stretched skin. Our results demonstrate that MSMs consisting of networks of matched fiber organization can predict the biomechanical behavior of mouse skin, including the large decrease in tissue volume and nonaffine fiber kinematics observed under uniaxial tension.

Biaxial mechanics of thermally denaturing skin - Part 2: Modeling

Understanding the response of skin to superphysiological temperatures is critical to the diagnosis and prognosis of thermal injuries, and to the development of temperature-based medical therapeutics. Unfortunately, this understanding has been hindered by our incomplete knowledge about the nonlinear coupling between skin temperature and its mechanics. In Part I of this study we experimentally demonstrated a complex interdependence of time, temperature, direction, and load in skin's response to superphysiological temperatures. In Part II of our study, we test two different models of skin's thermo-mechanics to explain our observations. In both models we assume that skin's response to superphysiological temperatures is governed by the denaturation of its highly collageneous microstructure. Thus, we capture skin's native mechanics via a microstructurally-motivated strain energy function which includes probability distributions for collagen fiber orientation and waviness. In the first model, we capture skin's response to superphysiological temperatures as a transition between two states that link the kinetics of collagen fiber denaturation to fiber coiling and to the transformation of each fiber's constitutive behavior from purely elastic to viscoelastic. In the second model, we capture skin's response to superphysiological temperatures instead via three states in which a sequence of two reactions link the kinetics of collagen fiber denaturation to fiber coiling, followed by a state of fiber damage. Given the success of both models in qualitatively and quantitatively capturing our observations, we expect that our work will provide guidance for future experiments that could probe each model's assumptions toward a better understanding of skin's coupled thermo-mechanics and that our work will be used to guide the engineering design of heat treatment therapies. STATEMENT OF SIGNIFICANCE: Quantifying and modeling skin thermo-mechanics is critical to our understanding of skin physiology, pathophysiology, as well as heat-based treatments. This work addresses a lack of theoretical and computational models of the coupled thermo-mechanics of skin. Our model accounts for skin microstructure through modeling the probability of fiber orientation and fiber stress-free states. Denaturing induces changes in the stress-free configuration of collagen, as well as changes in fiber stiffness and viscoelastic properties. We propose two competing models that fit all of our experimental observations. These models will enable future developments of thermal-therapeutics, prevention and management of skin thermal injuries, and set a foundation for improved mechanistic models of skin thermo-mechanics.

Effect of collagen fibre orientation on the Poisson's ratio and stress relaxation of skin: an ex vivo and in vivo study

During surgical treatment skin undergoes extensive deformation, hence it must be able to withstand large mechanical stresses without damage. Therefore, understanding the mechanical properties of skin becomes important. A detailed investigation on the relationship between the three-dimensional deformation response of skin and its microstructure is conducted in the current study. This study also discloses the underlying science of skin viscoelasticity. Deformation response of skin is captured using digital image correlation, whereas micro-CT, scanning electron microscopy and atomic force microscopy are used for microstructure analysis. Skin shows a large lateral contraction and expansion (auxeticity) when stretched parallel and perpendicular to the skin tension lines, respectively. Large lateral contraction is a result of fluid exudation from the tissue, while large rotation of the stiff collagen fibres in the loading direction explains the skin auxeticity. During stress relaxation, lateral contraction and fluid effluxion from skin reveal that tissue volume loss is the intrinsic science of skin viscoelasticity. Furthermore, the results obtained from in vivo study on human skin show the relevance of the ex vivo study to physiological conditions and stretching of the skin during its treatments.

Systems of conductive skin for power transfer in clinical applications

The primary aim of this article is to review the clinical challenges related to the supply of power in implanted left ventricular assist devices (LVADs) by means of transcutaneous drivelines. In effect of that, we present the preventive measures and post-operative protocols that are regularly employed to address the leading problem of driveline infections. Due to the lack of reliable wireless solutions for power transfer in LVADs, the development of new driveline configurations remains at the forefront of different strategies that aim to power LVADs in a less destructive manner. To this end, skin damage and breach formation around transcutaneous LVAD drivelines represent key challenges before improving the current standard of care. For this reason, we assess recent strategies on the surface functionalization of LVAD drivelines, which aim to limit the incidence of driveline infection by directing the responses of the skin tissue. Moreover, we propose a class of power transfer systems that could leverage the ability of skin tissue to effectively heal short diameter wounds. In this direction, we employed a novel method to generate thin conductive wires of controllable surface topography with the potential to minimize skin disruption and eliminate the problem of driveline infections. Our initial results suggest the viability of the small diameter wires for the investigation of new power transfer systems for LVADs. Overall, this review uniquely compiles a diverse number of topics with the aim to instigate new research ventures on the design of power transfer systems for IMDs, and specifically LVADs.

Improving reconstructive surgery design using Gaussian process surrogates to capture material behavior uncertainty

To produce functional, aesthetically natural results, reconstructive surgeries must be planned to minimize stress as excessive loads near wounds have been shown to produce pathological scarring and other complications (Gurtner et al., 2011). Presently, stress cannot easily be measured in the operating room. Consequently, surgeons rely on intuition and experience (Paul et al., 2016; Buchanan et al., 2016). Predictive computational tools are ideal candidates for surgery planning. Finite element (FE) simulations have shown promise in predicting stress fields on large skin patches and in complex cases, helping to identify potential regions of complication. Unfortunately, these simulations are computationally expensive and deterministic (Lee et al., 2018a). However, running a few, well selected FE simulations allows us to create Gaussian process (GP) surrogate models of local cutaneous flaps that are computationally efficient and able to predict stress and strain for arbitrary material parameters. Here, we create GP surrogates for the advancement, rotation, and transposition flaps. We then use the predictive capability of these surrogates to perform a global sensitivity analysis, ultimately showing that fiber direction has the most significant impact on strain field variations. We then perform an optimization to determine the optimal fiber direction for each flap for three different objectives driven by clinical guidelines (Leedy et al., 2005; Rohrer and Bhatia, 2005). While material properties are not controlled by the surgeon and are actually a source of uncertainty, the surgeon can in fact control the orientation of the flap with respect to the skin's relaxed tension lines, which are associated with the underlying fiber orientation (Borges, 1984). Therefore, fiber direction is the only material parameter that can be optimized clinically. The optimization task relies on the efficiency of the GP surrogates to calculate the expected cost of different strategies when the uncertainty of other material parameters is included. We propose optimal flap orientations for the three cost functions and that can help in reducing stress resulting from the surgery and ultimately reduce complications associated with excessive mechanical loading near wounds.

Plateau-Rayleigh instability in a soft viscoelastic material

A soft cylindrical interface endowed with surface tension can be unstable to wavy undulations. This is known as the Plateau-Rayleigh instability (PRI) and for solids the instability is governed by the competition between elasticity and capillarity. A dynamic stability analysis is performed for the cases of a soft (i) cylinder and (ii) cylindrical cavity assuming the material is viscoelastic with power-law rheology. The governing equations are made time-independent through the Laplace transform from which a solution is constructed using displacement potentials. The dispersion relationships are then derived, which depend upon the dimensionless elastocapillary number, solid Deborah number, and compressibility number, and the static stability limit, critical disturbance, and maximum growth rate are computed. This dynamic analysis recovers previous literature results in the appropriate limits. Elasticity stabilizes and compressibility destabilizes the PRI. For an incompressible material, viscoelasticity does not affect stability but does decrease the growth rate and shift the critical wavenumber to lower values. The critical wavenumber shows a more complex dependence upon compressibility for the cylinder but exhibits a predictable trend for the cylindrical cavity.

Visco- and poroelastic contributions of the zona pellucida to the mechanical response of oocytes

Probing mechanical properties of cells has been identified as a means to infer information on their current state, e.g. with respect to diseases or differentiation. Oocytes have gained particular interest, since mechanical parameters are considered potential indicators of the success of in vitro fertilisation procedures. Established tests provide the structural response of the oocyte resulting from the material properties of the cell's components and their disposition. Based on dedicated experiments and numerical simulations, we here provide novel insights on the origin of this response. In particular, polarised light microscopy is used to characterise the anisotropy of the zona pellucida, the outermost layer of the oocyte composed of glycoproteins. This information is combined with data on volumetric changes and the force measured in relaxation/cyclic, compression/indentation experiments to calibrate a multi-phasic hyper-viscoelastic model through inverse finite element analysis. These simulations capture the oocyte's overall force response, the distinct volume changes observed in the zona pellucida, and the structural alterations interpreted as a realignment of the glycoproteins with applied load. The analysis reveals the presence of two distinct timescales, roughly separated by three orders of magnitude, and associated with a rapid outflow of fluid across the external boundaries and a long-term, progressive relaxation of the glycoproteins, respectively. The new results allow breaking the overall response down into the contributions from fluid transport and the mechanical properties of the zona pellucida and ooplasm. In addition to the gain in fundamental knowledge, the outcome of this study may therefore serve an improved interpretation of the data obtained with current methods for mechanical oocyte characterisation.

Mechanical stimulation induces rapid fibroblast proliferation and accelerates the early maturation of human skin substitutes

The clinical treatment of large, full-thickness skin injuries with tissue-engineered autologous dermo-epidermal skin substitutes is an emerging alternative to split-thickness skin grafting. However, their production requires about one month of in vitro cell and tissue culture, which is a significant drawback for the treatment of patients with severe skin defects. With the aim to reduce the production time, we developed a new dynamic bioreactor setup that applies cyclic biaxial tension to collagen hydrogels for skin tissue engineering. By reliably controlling the time history of mechanical loading, the dynamic culturing results in a three-fold increase in collagen hydrogel stiffness and stimulates the embedded fibroblasts to enter the cell cycle. As a result, the number of fibroblasts is increased by 75% compared to under corresponding static culturing. Enhanced fibroblast proliferation promotes expression of dermal extracellular matrix proteins, keratinocyte proliferation, and the early establishment of the epidermis. The time required for early tissue maturation can therefore be reduced by one week. Analysis of the separate effects of cyclic loading, matrix stiffening, and interstitial fluid flow indicates that cyclic deformation is the dominant biophysical factor determining fibroblast proliferation, while tissue stiffening plays a lesser role. Local differences in the direction of deformation (in-plane equibiaxial vs. uniaxial strain) influence fibroblast orientation but not proliferation, nor the resulting tissue properties. Importantly, dynamic culturing does not activate fibroblast differentiation into myofibroblasts. The present work demonstrates that control of mechanobiological cues can be very effective in driving cell response toward a shorter production time for human skin substitutes.

Three-Dimensional Quantification of Collagen Microstructure During Tensile Mechanical Loading of Skin

Skin is a heterogeneous tissue that can undergo substantial structural and functional changes with age, disease, or following injury. Understanding how these changes impact the mechanical properties of skin requires three-dimensional (3D) quantification of the tissue microstructure and its kinematics. The goal of this study was to quantify these structure-function relationships via second harmonic generation (SHG) microscopy of mouse skin under tensile mechanical loading. Tissue deformation at the macro- and micro-scale was quantified, and a substantial decrease in tissue volume and a large Poisson’s ratio was detected with stretch, indicating the skin differs substantially from the hyperelastic material models historically used to explain its behavior. Additionally, the relative amount of measured strain did not significantly change between length scales, suggesting that the collagen fiber network is uniformly distributing applied strains. Analysis of undeformed collagen fiber organization and volume fraction revealed a length scale dependency for both metrics. 3D analysis of SHG volumes also showed that collagen fiber alignment increased in the direction of stretch, but fiber volume fraction did not change. Interestingly, 3D fiber kinematics was found to have a non-affine relationship with tissue deformation, and an affine transformation of the micro-scale fiber network overestimates the amount of fiber realignment. This result, along with the other outcomes, highlights the importance of accurate, scale-matched 3D experimental measurements when developing multi-scale models of skin mechanical function.

A biphasic multilayer computational model of human skin

The present study investigates the layer-specific mechanical behavior of human skin. Motivated by skin’s histology, a biphasic model is proposed which differentiates between epidermis, papillary and reticular dermis, and hypodermis. Inverse analysis of ex vivo tensile and in vivo suction experiments yields mechanical parameters for each layer and predicts a stiff reticular dermis and successively softer papillary dermis, epidermis and hypodermis. Layer-specific analysis of simulations underlines the dominating role of the reticular dermis in tensile loading. Furthermore, it shows that the observed out-of-plane deflection in ex vivo tensile tests is a direct consequence of the layered structure of skin. In in vivo suction experiments, the softer upper layers strongly influence the mechanical response, whose dissipative part is determined by interstitial fluid redistribution within the tissue. Magnetic resonance imaging-based visualization of skin deformation in suction experiments confirms the deformation pattern predicted by the multilayer model, showing a consistent decrease in dermal thickness for large probe opening diameters.

On the Reliability of Suction Measurements for Skin Characterization

in vivo skin characterization methods were shown to be useful in the detection of microstructural alterations of the dermis due to skin diseases. Specifically, the diagnostic potential of skin suction has been widely explored, yet measurement uncertainties prevented so far its application in clinical assessment. In this work, we analyze specific factors influencing the reliability of suction measurements. We recently proposed a novel suction device, called Nimble, addressing the limitations of existing instruments, and applied it in clinical trials quantifying mechanical differences between healthy skin and scars. Measurements were performed with the commercial device Cutometer and with the new device. A set of new suction measurements was carried out on scar tissue and healthy skin, and FE-based inverse analysis was applied to determine corresponding parameters of a hyperelastic-viscoelastic material model. FE simulations were used to rationalize differences between suction protocols and to analyze specific factors influencing the measurement procedure. Tissue stiffness obtained from Cutometer measurements was significantly higher compared to the one from Nimble measurements, which was shown to be associated with the higher deformation levels in the Cutometer and the nonlinear mechanical response of skin. The effect of the contact force exerted on skin during suction measurements was quantified, along with an analysis of the effectiveness of a corresponding correction procedure. Parametric studies demonstrated the inherently higher sensitivity of displacement- over load-controlled suction measurements, thus rationalizing the superior ability of the Nimble to distinguish between tissues.

How cosmetic tightening products modulate the biomechanics and morphology of human skin

The active and passive mechanical behavior of a cosmetic tightening product for skin anti-aging is investigated based on a wide range of in vivo and in vitro measurements. The experimental data are used to inform a numerical model of the attained cosmetic effect, which is then implemented in a commercial finite-element framework and used to analyze the mechanisms that regulate the biomechanical interaction between the native tissue and the tightening film. Such a film reduces wrinkles and enhances skin consistency by increasing its stiffness by 48-107% and reducing inelastic, non-recoverable deformations (-47%). The substrate deformability influences both the extent of tightening and the reduction of wrinkle amplitude. The present findings allow, for the first time, to rationalize the mechanisms of action of cosmetic products with a tightening action and provide quantitative evidence for further optimization of this fascinating class of biomaterials.

Extracting the elasticity of the human skin in microscale and in-vivo from atomic force microscopy experiments using viscoelastic models

Detecting mechanical properties of the intact skin in-vivo leads to a novel quantitative method to diagnose skin diseases and to monitor skin conditions in clinical settings. Current research and clinical methods that detect skin mechanics have major limitations. The in-vitro experiments are done in non-physiological conditions and in-vivo clinical methods measurer unwanted mechanics of underneath fat and muscle tissues but report the measurement as skin mechanics. An ideal skin mechanics should be captured at skin scale (i.e., micron-scale) and in-vivo. However, extreme challenges of capturing the in-vivo skin mechanics in micron-scale including skin motion due to heart beep, breathing and movement of the subject, has hindered measurement of skin mechanics in-vivo.This study for the first time captures micro-scale mechanics (elasticity and viscoelasticity) of top layers of skin (i.e., the stratum corneum (SC) and stratum granulosum (SG)) in-vivo. In this study, the relevant literature is reviewed and Atomic Force Microscopy (AFM) was used to capture force-indentation curves on the fingertip skin of four human subjects at a high indentation speed of 40 μm/s. The skin of the same subject were tested in-vitro at 10 different indentation speeds ranging from 0.125 to 40 μm/s by AFM. This study extracts the in-vivo elasticity of SC and SG by detecting time-dependency of tested tissue using a fractional viscoelastic standard linear model developed for indentation. The in-vivo elasticity of SC and SG were smaller in females and in-vitro elasticity were higher than that of in-vivo results. The results were consistent with previous observations.

The bioengineering theory of the key modes of action of a cyanoacrylate liquid skin protectant

The objective of this article is to formulate a new bioengineering theoretical framework for modelling the biomechanical efficacy of cyanoacrylate skin protectants, with specific focus on the Marathon technology (Medline Industries, Inc., Northfield, Illinois) and its modes of action. This work details the bioengineering and mathematical formulations of the theory, which is based on the classic engineering theories of flexural stiffness of coated elements and deformation friction. Based on the relevant skin anatomy and physiology, this paper demonstrates: (a) the contribution of the polymerised cyanoacrylate coating to flexural skin stiffness, which facilitates protection from non-axial (eg, compressive) localised mechanical forces; and (b) the contribution of the aforementioned coating to reduction in frictional forces and surface shear stresses applied by contacting objects such as medical devices. The present theoretical framework establishes that application of the cyanoacrylate coating provides considerable biomechanical protection to skin and subdermally, by shielding skin from both compressive and frictional (shearing) forces. Moreover, these analyses indicate that the prophylactic effects of the studied cyanoacrylate coating become particularly strong where the skin is thin or fragile (typically less than ~0.7 mm thick), which is characteristic to old age, post-neural injuries, neuromuscular diseases, and in disuse-induced tissue atrophy conditions.

Activin-mediated alterations of the fibroblast transcriptome and matrisome control the biomechanical properties of skin wounds

Matrix deposition is essential for wound repair, but when excessive, leads to hypertrophic scars and fibrosis. The factors that control matrix deposition in skin wounds have only partially been identified and the consequences of matrix alterations for the mechanical properties of wounds are largely unknown. Here, we report how a single diffusible factor, activin A, affects the healing process across scales. Bioinformatics analysis of wound fibroblast transcriptome data combined with biochemical and histopathological analyses of wounds and functional in vitro studies identify that activin promotes pro-fibrotic gene expression signatures and processes, including glycoprotein and proteoglycan biosynthesis, collagen deposition, and altered collagen cross-linking. As a consequence, activin strongly reduces the wound and scar deformability, as identified by a non-invasive in vivo method for biomechanical analysis. These results provide mechanistic insight into the roles of activin in wound repair and fibrosis and identify the functional consequences of alterations in the wound matrisome at the biomechanical level.

The relationship between histopathology, gene expression, and biochemical and mechanical properties of wounds is largely unknown. Here, the authors show that activin A alters wound healing at multiple levels by promoting pro-fibrotic gene expression and matrix deposition, thereby affecting biomechanical properties of skin wounds.

The regional-dependent biaxial behavior of young and aged mouse skin: A detailed histomechanical characterization, residual strain analysis, and constitutive model

Skin fulfills several vital functions, many of which are dependent on its mechanical properties. Therefore, as mice have become an invaluable model for skin research, determining murine skin's mechanical properties is important. Specifically, skin's mechanical properties are important for functional tests as well as for prognostic and diagnostic purposes. Additionally, computational simulations of skin behavior are becoming commonplace, rendering accurate models of murine skin's constitutive behavior necessary. To date, our knowledge of mouse skin mechanics shows significant gaps. For example, there are no comprehensive reports correlating skin's mechanical properties with region, age, and direction. Moreover, mouse skin's residual strain behavior has not been reported on. In our current work, we set out to fill these gaps. Based on histology, 2-photon microscopy, and planar biaxial testing, while accurately tracking various reference configurations, we report on differences in gross structure, microstructural organization, and constitutive response of skin, and cast those properties into a versatile Fung-type hyperelastic constitutive law for three reference configurations. Our data is the most comprehensive report contrasting the mechanical properties of young (12 weeks) and aged (52 weeks) mouse skin and will, thus, be valuable to basic science as control data, and provide accurate constitutive laws for mouse skin modeling. STATEMENT OF SIGNIFICANCE: Our findings are significant as they fill several gaps in our understanding of mouse skin mechanics. This is particularly important as mouse skin is becoming a frequent and critical model of human skin for cosmetic and medical science. Specifically, we quantified how mechanical properties of mice skin vary with age, with location, and with direction. Additionally, we cast our findings into constitutive models that can be used by others for predictive computer simulations of skin behavior.