Adrian Buganza Tepole

Skin growth and remodeling underlies health, disease, and treatments such as tissue expansion (TE). The mechanotransduction pathways in dermal fibroblasts are increasingly well characterized, and tissue-level growth has been described phenomenologically, but coupling between cell-level signaling and tissue-level growth remains poorly understood. We develop a dermal fibroblast signaling network through extensive literature data curation, comprising 151 reactions among 96 nodes. The inputs are mechanical stretch and eight ligands (TGFbeta, PDGF, FGF, IL1, IL6, TNF\alpha, AngII, ET1); outputs of interest span ECM-enzymes (proMMP1/2/9, MMP1/2/9), ECM proteins (CImRNA, collagen I, fibronectin), and fibroblast activity (alphaSMA, proliferation). Implemented as a logic-based ODE system, the network reproduces 82% of the calibration dataset and agrees with independent validation data. Sensitivity analysis reveals a tension-dependent regulation of signaling: at baseline tension, outputs are governed by many boosters (nodes that positively influence downstream targets) and one dominant brake, LATS1/2, whereas at high tension control consolidates and new, tension-specific regulators such as integrin (ITGB1) emerge. Multiple pathway axes converge on a few central regulators, producing pronounced crosstalk, most notably between TGFbeta and mechanical tension. Finally, linking the collagen outputs to a tissue-level growth formulation yields a bidirectional mechanical-biochemical coupling that reproduces tension-induced skin growth measured in a porcine TE model. This framework establishes a comprehensively calibrated dermal fibroblast signaling network coupled to tissue-level growth, opening opportunities for targeted TE interventions.

Pressure ulcers arise from sustained mechanical loading that impairs perfusion and damages skin tissues, yet the coupled mechanical and biological mechanisms of their formation and healing remain poorly characterized. We addressed this gap using a mouse model in which dorsal skin underwent 72 hours of magnet-induced ischemia followed by reperfusion, with tissue collected at 0, 3, 6, and 9 days and compared with baseline controls. From each mouse, we obtained paired samples from pressure ulcer and remote control (non-loaded) sites, mapped thickness by tissue profilometry, and performed equibiaxial testing with full-field digital image correlation and inverse finite element analysis to estimate regional material parameters. In parallel, we quantified CD31 ⁺ vasculature, F4/80 ⁺ macrophages, collagen content, and key cytokines. Pressure ulcer sites were compressed and thinner at Day 0, developed ulcers by Day 3, and continued to remodel through Day 9. Mechanical tests revealed heterogeneous strain fields with elevated deformation along ulcer borders, while remote control tissue deformed more homogeneously. These mechanical changes evolved alongside dynamic vessel and macrophage repopulation, increased collagen content at early time points, and cytokine upregulation within pressure ulcer tissue. Collectively, our data define the spatiotemporal co-evolution of tissue geometry, mechanics, collagen remodeling, and inflammation in pressure ulcers and provide a quantitative foundation for predictive mechanobiological models.

SUMMARY: Breast-conserving surgery (BCS; lumpectomy) is widely used to treat early-stage breast cancer, yet the complexity and patient-specific variability of postoperative cavity remodeling make healing trajectories and physical outcomes difficult to predict. Although inflammatory and vascular processes are central to these outcomes, mathematical models of tissue healing do not capture their coupled interactions or calibrate them against experimental data. Here, we extend a computational model of breast cavity healing to incorporate coupled inflammatory and vascular dynamics, including angiogenesis, oxygen transport, and inflammatory cell activity. The model is calibrated using preclinical porcine lumpectomy histology and literature data. Model parameters are inferred using a multi-task Gaussian process surrogate within a Bayesian inference frame-work to align predictions with experimental observations and quantify uncertainty. Thus, this work provides a mechanistically grounded framework for inflammatory and vascular remodeling during breast cavity healing and provides a foundation for patient-specific prediction of healing and physical outcomes following lumpectomy.

Angiogenesis drives the formation and remodeling of capillary networks throughout tissue repair, regulating the vascular environment that supports healing and tissue remodeling. Experimental characterization of these processes is commonly performed using CD31-stained histological tissue sections to quantify capillary surface density and morphology throughout healing. However, these measurements provide only two-dimensional characterization of an underlying three-dimensional (3D) vascular network, limiting direct estimation of volumetric capillary density and vascular architecture. To address this limitation, an experimentally informed framework was developed to generate representative 3D capillary networks, enabling estimation of volumetric capillary density from histologically quantified vascular measurements. CD31-stained histological sections obtained from a longitudinal porcine lumpectomy study were analyzed to quantify the percentage of CD31-positive area (%CD31+) and capillary morphology within healthy tissue and healing surgical cavities. Histologically quantified morphology distributions and literature-informed vascular branching characteristics were incorporated into a capillary network generation framework to construct representative 3D vascular networks. Capillary branches were iteratively generated within representative tissue volumes until virtual histological sections reproduced experimental %CD31+ measurements, enabling estimation of volumetric capillary density. Generated capillary networks demonstrated good agreement with experimentally characterized 3D vascular architecture, while simulated histological sections accurately reproduced experimentally quantified capillary counts and vascularization measurements. Application of the framework to the porcine lumpectomy dataset captured temporal changes in vascular remodeling throughout healing, revealing progressive increases in volumetric capillary density and vascular maturation. Collectively, this framework provides an experimentally informed methodology for relating histological vascular measurements to volumetric capillary density estimates, supporting future computational studies of angiogenesis and tissue repair.

Accurate in vivo characterization of skin mechanical properties is essential for diagnostics and treatment planning across dermatological and surgical applications. Existing noninvasive techniques are limited in capturing the nonlinear and anisotropic behavior of skin. In this work, we propose a Bayesian inference framework that leverages active membranes to induce desired deformations and infer patient-specific skin properties from a measured strain field. A finite element model of skin–membrane interaction, parameterized using the Holzapfel–Gasser–Ogden model, is used to generate strain field data under various membrane actuation conditions. To overcome the computational cost of repeated simulations required for Bayesian sampling, we construct a data-driven surrogate using principal component analysis for dimensionality reduction and Gaussian process regression for rapid evaluation. Our approach enables probabilistic inference of key skin parameters, including shear modulus, fiber stiffness, dispersion, and orientation. An advantage of the proposed method is that inference of skin biomechanics does not require direct force measurements; rather, the method relies on known properties of active membranes (which can be tested ahead of time). The method does require strain field measurements. Through synthetic studies, we demonstrate that our method accurately recovers most model parameters even under moderate levels of spatially correlated noise, and that multiframe or multimembrane observations significantly enhance identifiability. These results establish the potential of active membranes as a viable platform for noninvasive, in vivo skin biomechanics assessment.