Multiscale Computational Model Predicts Mouse Skin Kinematics Under Tensile Loading

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.

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.

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.

A Systems Approach to Biomechanics, Mechanobiology, and Biotransport

The human body represents a collection of interacting systems that range in scale from nanometers to meters. Investigations from a systems perspective focus on how the parts work together to enact changes across spatial scales, and further our understanding of how systems function and fail. Here, we highlight systems approaches presented at the 2022 Summer Biomechanics, Bio-engineering, and Biotransport Conference in the areas of solid mechanics; fluid mechanics; tissue and cellular engineering; biotransport; and design, dynamics, and rehabilitation; and biomechanics education. Systems approaches are yielding new insights into human biology by leveraging state-of-the-art tools, which could ultimately lead to more informed design of therapies and medical devices for preventing and treating disease as well as rehabilitating patients using strategies that are uniquely optimized for each patient. Educational approaches can also be designed to foster a foundation of systems-level thinking.

Elastic fibers define embryonic tissue stiffness to enable buckling morphogenesis of the small intestine

During embryonic development, tissues must possess precise material properties to ensure that cell-generated forces give rise to the stereotyped morphologies of developing organs. However, the question of how material properties are established and regulated during development remains understudied. Here, we aim to address these broader questions through the study of intestinal looping, a process by which the initially straight intestinal tube buckles into loops, permitting ordered packing within the body cavity. Looping results from elongation of the tube against the constraint of an attached tissue, the dorsal mesentery, which is elastically stretched by the elongating tube to nearly triple its length. This elastic energy storage allows the mesentery to provide stable compressive forces that ultimately buckle the tube into loops. Beginning with a transcriptomic analysis of the mesentery, we identified widespread upregulation of extracellular matrix related genes during looping, including genes related to elastic fiber deposition. Combining molecular and mechanical analyses, we conclude that elastin confers tensile stiffness to the mesentery, enabling its mechanical role in organizing the developing small intestine. These results shed light on the role of elastin as a driver of morphogenesis that extends beyond its more established role in resisting cyclic deformation in adult tissues.

The non-affine fiber network solver: A multiscale fiber network material model for finite-element analysis

Multiscale mechanical models in biomaterials research have largely relied on simplifying the microstructure in order to make large-scale simulations tractable. The microscale simplifications often rely on approximations of the constituent distributions and assumptions on the deformation of the constituents. Of particular interest in biomechanics are fiber embedded materials, where simplified fiber distributions and assumed affinity in the fiber deformation greatly influence the mechanical behavior. The consequences of these assumptions are problematic when dealing with microscale mechanical phenomena such as cellular mechanotransduction in growth and remodeling, and fiber-level failure events during tissue failure. In this work, we propose a technique for coupling non-affine network models to finite element solvers, allowing for simulation of discrete microstructural phenomena within macroscopically complex geometries. The developed plugin is readily available as an open-source library for use with the bio-focused finite element software FEBio, and the description of the implementation allows for the adaptation to other finite element solvers.