A burn wound is a complex systemic disease at multiple levels. Current knowledge of scar formation after burn injury has come from traditional biological and clinical studies. These are normally focused on just a small part of the entire process, which has limited our ability to sufficiently understand the underlying mechanisms and to predict systems behaviour. Scar formation after burn injury is a result of a complex biological system-wound healing. It is a part of a larger whole. In this self-organising system, many components form networks of interactions with each other. These networks of interactions are typically non-linear and change their states dynamically, responding to the environment and showing emergent long-term behaviour. How molecular and cellular data relate to clinical phenomena, especially regarding effective therapies of burn wounds to achieve minimal scarring, is difficult to unravel and comprehend. Complexity science can help bridge this gap by integrating small parts into a larger whole, such that relevant biological mechanisms and data are combined in a computational model to better understand the complexity of the entire biological system. A better understanding of the complex biological system of post-burn scar formation could bring research and treatment regimens to the next level. The aim of this review/position paper is to create more awareness of complexity in scar formation after burn injury by describing the basic principles of complexity science and its potential for burn care professionals. Declaration of interest: The authors gratefully acknowledge financial support from the Dutch Burns Foundation (DBF), Beverwijk, project 19.105. The authors have no conflicts of interest to declare.

@en

A continuum hypothesis-based model is developed for the simulation of the contraction of burns in order to gain new insights into which elements of the healing response might have a substantial influence on this process. Tissue is modeled as a neo-Hookean solid. Furthermore, (myo)fibroblasts, collagen molecules, and a generic signaling molecule are selected as model components. An overview of the custom-made numerical algorithm is presented. Subsequently, good agreement is demonstrated with respect to variability in the evolution of the surface area of burns over time between the outcomes of computer simulations and measurements obtained in an experimental study. In the model this variability is caused by varying the values for some of its parameters simultaneously. A factorial design combined with a regression analysis are used to quantify the individual contributions of these parameter value variations to the dispersion in the surface area of healing burns. The analysis shows that almost all variability in the surface area can be explained by variability in the value for the myofibroblast apoptosis rate and, to a lesser extent, the value for the collagen molecule secretion rate. This suggests that most of the variability in the evolution of the surface area of burns over time in the experimental study might be attributed to variability in these two rates. Finally, a probabilistic analysis is used in order to investigate in more detail the effect of variability in the values for the two rates on the healing process. Results of this analysis are presented and discussed.@en

A continuum hypothesis-based, biomechanical model is presented for the simulation of the collagen bundle distribution-dependent contraction and subsequent retraction of healing dermal wounds that cover a large surface area. Since wound contraction mainly takes place in the dermal layer of the skin, solely a portion of this layer is included explicitly into the model. This portion of dermal layer is modeled as a heterogeneous, orthotropic continuous solid with bulk mechanical properties that are locally dependent on both the local concentration and the local geometrical arrangement of the collagen bundles. With respect to the dynamic regulation of the geometrical arrangement of the collagen bundles, it is assumed that a portion of the collagen molecules are deposited and reoriented in the direction of movement of (myo)fibroblasts. The remainder of the newly secreted collagen molecules are deposited by ratio in the direction of the present collagen bundles. Simulation results show that the distribution of the collagen bundles influences the evolution over time of both the shape of the wounded area and the degree of overall contraction of the wounded area. Interestingly, these effects are solely a consequence of alterations in the initial overall distribution of the collagen bundles, and not a consequence of alterations in the evolution over time of the different cell densities and concentrations of the modeled constituents. In accordance with experimental observations, simulation results show furthermore that ultimately the majority of the collagen molecules ends up permanently oriented toward the center of the wound and in the plane that runs parallel to the surface of the skin.@en

A continuum hypothesis-based model is presented for the simulation of the formation and the subsequent regression of hypertrophic scar tissue after dermal wounding. Solely the dermal layer of the skin is modeled explicitly and it is modeled as a heterogeneous, isotropic and compressible neo-Hookean solid. With respect to the constituents of the dermal layer, the following components are selected as primary model components: fibroblasts, myofibroblasts, a generic signaling molecule and collagen molecules. A good match with respect to the evolution of the thickness of the dermal layer of scars between the outcomes of simulations and clinical measurements on hypertrophic scars at different time points after injury in human subjects is demonstrated. Interestingly, the comparison between the outcomes of the simulations and the clinical measurements demonstrates that a relatively high apoptosis rate of myofibroblasts results in scar tissue that behaves more like normal scar tissue with respect to the evolution of the thickness of the tissue over time, while a relatively low apoptosis rate results in scar tissue that behaves like hypertrophic scar tissue with respect to the evolution of the thickness of the tissue over time. Our ultimate goal is to construct models with which the properties of newly generated tissues that form during wound healing can be predicted with a high degree of certainty. The development of the presented model is considered by us as a step toward their construction.@en