Search

Continuing advances in genomics, molecular and cellular mechanobiology and immunobiology, including transcriptomics and proteomics, and biomechanics increasingly reveal the complexity underlying native tissue and organ structure and function. Identifying methods to repair, regenerate, or replace vital tissues and organs remains one of the greatest challenges of modern biomedical engineering, one that deserves our very best effort. Notwithstanding the continuing need for improving standard methods of investigation, including cell, organoid, and tissue culture, biomaterials development and fabrication, animal models, and clinical research, it is increasingly evident that modern computational methods should play increasingly greater roles in advancing the basic science, bioengineering, and clinical application of regenerative medicine. This brief review focuses on the development and application of computational models of tissue and organ mechanobiology and mechanics for purposes of designing tissue engineered constructs and understanding their development in vitro and in situ. Although the basic approaches are general, for illustrative purposes we describe two recent examples from cardiovascular medicine—tissue engineered heart valves (TEHVs) and tissue engineered vascular grafts (TEVGs)—to highlight current methods of approach as well as continuing needs.

Blood vessels grow and remodel in response to mechanical stimuli. Many computational models capture this process phenomenologically, by assuming stress homeostasis, but this approach cannot unravel the underlying cellular mechanisms. Mechano-sensitive Notch signaling is well-known to be key in vascular development and homeostasis. Here, we present a multiscale framework coupling a constrained mixture model, capturing the mechanics and turnover of arterial constituents, to a cell–cell signaling model, describing Notch signaling dynamics among vascular smooth muscle cells (SMCs) as influenced by mechanical stimuli. Tissue turnover was regulated by both Notch activity, informed by in vitro data, and a phenomenological contribution, accounting for mechanisms other than Notch. This novel framework predicted changes in wall thickness and arterial composition in response to hypertension similar to previous in vivo data. The simulations suggested that Notch contributes to arterial growth in hypertension mainly by promoting SMC proliferation, while other mechanisms are needed to fully capture remodeling. The results also indicated that interventions to Notch, such as external Jagged ligands, can alter both the geometry and composition of hypertensive vessels, especially in the short term. Overall, our model enables a deeper analysis of the role of Notch and Notch interventions in arterial growth and remodeling and could be adopted to investigate therapeutic strategies and optimize vascular regeneration protocols.

Rebuilding the kidney in the context of tissue engineering offers a major challenge as the organ is structurally complex and has a high variety of specific functions. Recreation of kidney function is inherently connected to the formation of tubules since the functional subunit of the kidney, the nephron, is based on tubular structures. In vivo, tubulogenesis culminates in a perfectly shaped, patterned, and functional renal tubule via different morphogenic processes that depend on delicately orchestrated chemical, physical, and mechanical interactions between cells and between cells and their microenvironment. This review summarizes the current understanding of the role of the microenvironment in the morphogenic processes involved in in vivo renal tubulogenesis. We highlight the current state-of-the-art of renal tubular engineering and provide a view on the design elements that can be extracted from these studies. Next, we discuss how computational modeling can aid in specifying and identifying design parameters and provide directions on how these design parameters can be incorporated in biomaterials for the purpose of engineering renal tubulogenesis. Finally, we propose that a step-by-step reciprocal interaction between understanding and engineering is necessary to effectively guide renal tubulogenesis. Impact statement Tubular tissue engineering lies at the foundation of regenerating kidney tissue function, as the functional subunit of the kidney, the nephron, is based on tubular structures. Guiding renal tubulogenesis toward functional renal tubules requires in-depth knowledge of the developmental processes that lead to the formation of native tubules as well as engineering approaches to steer these processes. In this study, we review the role of the microenvironment in the developmental processes that lead to functional renal tubules and give directions how this knowledge can be harnessed for biomaterial-based tubular engineering using computational models.

In recent years, the mathematical and computational sciences have developed novel methodologies and insights that can aid in designing advanced bioreactors, microfluidic setups or organ-on-chip devices, in optimizing culture conditions, or predicting long-term behavior of engineered tissues in vivo. In this review, we introduce the concept of computational models and how they can be integrated in an interdisciplinary workflow for Tissue Engineering and Regenerative Medicine (TERM). We specifically aim this review of general concepts and examples at experimental scientists with little or no computational modeling experience. We also describe the contribution of computational models in understanding TERM processes and in advancing the TERM field by providing novel insights. Impact StatementAlthough in recent years the use of mathematical and computational sciences has increased in the Tissue Engineering and Regenerative Medicine (TERM) field, we believe that a further integration of experimental and computational approaches has a huge potential for advancing the field due to the ability of models to explain and predict experimental results and efficiently optimize TERM product and process designs. By providing an overview of existing computational models, how they have contributed to the field, as well as a future perspective, this review represents an important step to help realize TERM's ultimate goal: a cure instead of care.

Continuing advances in genomics, molecular and cellular mechanobiology and immunobiology, including transcriptomics and proteomics, and biomechanics increasingly reveal the complexity underlying native tissue and organ structure and function. Identifying methods to repair, regenerate, or replace vital tissues and organs remains one of the greatest challenges of modern biomedical engineering, one that deserves our very best effort. Notwithstanding the continuing need for improving standard methods of investigation, including cell, organoid, and tissue culture, biomaterials development and fabrication, animal models, and clinical research, it is increasingly evident that modern computational methods should play increasingly greater roles in advancing the basic science, bioengineering, and clinical application of regenerative medicine. This brief review focuses on the development and application of computational models of tissue and organ mechanobiology and mechanics for purposes of designing tissue engineered constructs and understanding their development in vitro and in situ. Although the basic approaches are general, for illustrative purposes we describe two recent examples from cardiovascular medicine-tissue engineered heart valves (TEHVs) and tissue engineered vascular grafts (TEVGs)-to highlight current methods of approach as well as continuing needs.

Understanding cardiovascular growth and remodeling (G&R) is fundamental for designing robust cardiovascular tissue engineering strategies, which enable synthetic or biological scaffolds to transform into healthy living tissues after implantation. Computational modeling, particularly when integrated with experimental research, is key for advancing our understanding, predicting the in vivo evolution of engineered tissues, and efficiently optimizing scaffold designs. As cells are ultimately the drivers of G&R and known to change their behavior in response to mechanical cues, increasing efforts are currently undertaken to capture (mechano-mediated) cell behavior in computational models. In this selective review, we highlight some recent examples that are relevant in the context of cardiovascular tissue engineering and discuss the current and future biological and computational challenges for modeling cell-mediated G&R.

Prestretch is observed in many soft biological tissues, directly influencing the mechanical behavior of the tissue in question. The development of this prestretch occurs through complex growth and remodeling phenomena, which yet remain to be elucidated. In the present study it was investigated whether local cell-mediated traction forces can explain the development of global anisotropic tissue prestretch in the mitral valve. Towards this end, a model predicting actin stress fiber-generated traction forces was implemented in a finite element framework of the mitral valve. The overall predicted magnitude of prestretch induced valvular contraction after release of in vivo boundary constraints was in good agreement with data reported on valvular retraction after excision from the heart. Next, by using a systematic variation of model parameters and structural properties, a more anisotropic prestretch development in the valve could be obtained, which was also similar to physiological values. In conclusion, this study shows that cell-generated traction forces could explain prestretch magnitude and anisotropy in the mitral valve.

Collagen fibre degradation is a strain-dependent process, whereby the magnitude of experienced strain dictates the rate of enzymatic cleavage. Studies have identified conflicting findings as to whether strain inhibits or enhances collagen degradation, which may be explained by the tissue type and tissue scale investigated, as well as the strain range considered. The aim of this study is to identify, for the first time, the strain-dependent degradation response of intact arterial vessels experiencing physiological pressures and apply these findings to a computational model to better understand degenerative arterial diseases, such as aneurysms. To achieve this, a series of quasi-static pressure inflation experiments were carried out on intact arteries in the presence of purified bacterial collagenase at physiologically relevant pressures to investigate collagen matrix degradation in the vascular wall. A complementary computational model was developed to explore the complex role of pressure, non-collagenous matrix contribution, and collagen fibre crimp in the ultimate degradation response of the vessel. Pressure induced inflation-degradation results identified an increased rate of vessel expansion and reduced time to failure with increasing pressure in the vessels. Interestingly, our computational model was able to capture this same response, including the elevated rates of degradation which occur at low pressures. These findings highlight the critical role of strain in collagen degradation, particularly in cases of arterial disease, such as aneurysm formation, whereby structural integrity may be compromised.

Valvular heart disease is a major cause of morbidity and mortality worldwide. Surgical valve repair or replacement has been the standard of care for patients with valvular heart disease for many decades, but transcatheter heart valve therapy has revolutionized the field in the past 15 years. However, despite the tremendous technical evolution of transcatheter heart valves, to date, the clinically available heart valve prostheses for surgical and transcatheter replacement have considerable limitations. The design of next-generation tissue-engineered heart valves (TEHVs) with repair, remodelling and regenerative capacity can address these limitations, and TEHVs could become a promising therapeutic alternative for patients with valvular disease. In this Review, we present a comprehensive overview of current clinically adopted heart valve replacement options, with a focus on transcatheter prostheses. We discuss the various concepts of heart valve tissue engineering underlying the design of next-generation TEHVs, focusing on off-the-shelf technologies. We also summarize the latest preclinical and clinical evidence for the use of these TEHVs and describe the current scientific, regulatory and clinical challenges associated with the safe and broad clinical translation of this technology.

Adherent cells are generally able to reorient in response to cyclic strain. In three-dimensional tissues, however, extracellular collagen can affect this cellular response. In this study, a computational model able to predict the combined effects of mechanical stimuli and collagen on cellular (re)orientation was developed. In particular, a recently proposed computational model (which only accounts for mechanical stimuli) was extended by considering two hypotheses on how collagen influences cellular (re)orientation: collagen contributes to cell alignment by providing topographical cues (contact guidance); or collagen causes a spatial obstruction for cellular reorientation (steric hindrance). In addition, we developed an evolution law to predict cell-induced collagen realignment. The hypotheses were tested by simulating bi- or uniaxially constrained cell-populated collagen gels with different collagen densities, subjected to immediate or delayed uniaxial cyclic strain with varying strain amplitudes. The simulation outcomes are in agreement with previous experimental reports. Taken together, our computational approach is a promising tool to understand and predict the remodeling of collagenous tissues, such as native or tissue-engineered arteries and heart valves.

Tissue growth and remodeling are known to govern mechanical homeostasis in biological tissue, but their relative contributions to homeostasis remain unclear. Here, we use mechanical models, fueled by experimental findings, to demonstrate that growth and remodeling have different effects on heart valve stretch homeostasis during physiological postnatal development. Two developmental stages were considered: early-stage (from infant to adolescent) and late-stage (from adolescent to adult) development. Our models indicated that growth and remodeling play opposing roles in preserving tissue stretch and with time. During early-stage development, excessive tissue stretch was decreased by tissue growth and increased by remodeling. In contrast, during late-stage development tissue stretch was decreased by remodeling and increased by growth. Our findings contribute to an improved understanding of native heart valve adaptation throughout life, and are highly relevant for the development of tissue-engineered heart valves.

The intermediate filament (IF) cytoskeleton has been proposed to regulate morphogenic processes by integrating the cell fate signaling machinery with mechanical cues. Signaling between endothelial cells (ECs) and vascular smooth muscle cells (VSMCs) through the Notch pathway regulates arterial remodeling in response to changes in blood flow. Here we show that the IF-protein vimentin regulates Notch signaling strength and arterial remodeling in response to hemodynamic forces. Vimentin is important for Notch transactivation by ECs and vimentin knockout mice (VimKO) display disrupted VSMC differentiation and adverse remodeling in aortic explants and in vivo. Shear stress increases Jagged1 levels and Notch activation in a vimentin-dependent manner. Shear stress induces phosphorylation of vimentin at serine 38 and phosphorylated vimentin interacts with Jagged1 and increases Notch activation potential. Reduced Jagged1-Notch transactivation strength disrupts lateral signal induction through the arterial wall leading to adverse remodeling. Taken together we demonstrate that vimentin forms a central part of a mechanochemical transduction pathway that regulates multilayer communication and structural homeostasis of the arterial wall.

Understanding collagen and stress fiber remodeling is essential for the development of engineered tissues with good functionality. These processes are complex, highly interrelated, and occur over different time scales. As a result, excessive computational costs are required to computationally predict the final organization of these fibers in response to dynamic mechanical conditions. In this study, an analytical approximation of a stress fiber remodeling evolution law was derived. A comparison of the developed technique with the direct numerical integration of the evolution law showed relatively small differences in results, and the proposed method is one to two orders of magnitude faster.

Summary In the presence of anisotropic biochemical or topographical patterns, cells tend to align in the direction of these cues—a widely reported phenomenon known as “contact guidance.” To investigate the origins of contact guidance, here, we created substrates micropatterned with parallel lines of fibronectin with dimensions spanning multiple orders of magnitude. Quantitative morphometric analysis of our experimental data reveals two regimes of contact guidance governed by the length scale of the cues that cannot be explained by enforced alignment of focal adhesions. Adopting computational simulations of cell remodeling on inhomogeneous substrates based on a statistical mechanics framework for living cells, we show that contact guidance emerges from anisotropic cell shape fluctuation and “gap avoidance,” i.e., the energetic penalty of cell adhesions on non-adhesive gaps. Our findings therefore point to general biophysical mechanisms underlying cellular contact guidance, without the necessity of invoking specific molecular pathways., Graphical Abstract, Highlights Cells exhibit an optimum alignment (contact guidance) to substrate micropatterns Contact guidance does not require constrained alignment of focal adhesions Statistical mechanics model reveals the energetic role of non-adhesive gaps The non-adhesive gap size determines the degree of cell alignment, Buskermolen et al. report that cells align to substrate micropatterns in a highly length-scale-dependent manner. This alignment emerges from energetically favorable cell elongation and “gap avoidance”: minimization of cellular adhesions on non-adhesive gaps. This finding points to biophysically driven mechanisms underlying cellular contact guidance without necessarily invoking specific molecular pathways.

The hemodynamic functionality of heart valves strongly depends on the distribution of collagen fibers, which are their main load-bearing constituents. It is known that collagen networks remodel in response to mechanical stimuli. Yet, the complex interplay between external load and collagen remodeling is poorly understood. In this study, we adopted a computational approach to simulate collagen remodeling occurring in native fetal and pediatric heart valves. The computational model accounted for several biological phenomena: cellular (re)orientation in response to both mechanical stimuli and topographical cues provided by collagen fibers; collagen deposition and traction forces along the main cellular direction; collagen degradation decreasing with stretch; and cell-mediated collagen prestretch. Importantly, the computational results were well in agreement with previous experimental data for all simulated heart valves. Simulations performed by varying some of the computational parameters suggest that cellular forces and (re)orientation in response to mechanical stimuli may be fundamental mechanisms for the emergence of the circumferential collagen alignment usually observed in native heart valves. On the other hand, the tendency of cells to coalign with collagen fibers is essential to maintain and reinforce that circumferential alignment during development. Statement of Significance: The hemodynamic functionality of heart valves is strongly influenced by the alignment of load-bearing collagen fibers. Currently, the mechanisms that are responsible for the development of the circumferential collagen alignment in native heart valves are not fully understood. In the present study, cell-mediated remodeling of native human heart valves during early development was computationally simulated to understand the impact of individual mechanisms on collagen alignment. Our simulations successfully predicted the degree of collagen alignment observed in native fetal and pediatric semilunar valves. The computational results suggest that the circumferential collagen alignment arises from cell traction and cellular (re)orientation in response to mechanical stimuli, and with increasing age is reinforced by the tendency of cells to co-align with pre-existing collagen fibers.

Collagen fibre remodelling is a strain dependent process which is stimulated by the degradation of existing collagen. To date, literature has focussed on strain dependent degradation of pure collagen or structurally simple collagenous tissues, often overlooking degradation within more complex, heterogenous soft tissues. The aim of this study is to identify, for the first time, the strain dependent degradation behaviour and mechanical factors influencing collagen degradation in arterial tissue using a combined experimental and numerical approach. To achieve this, structural analysis was carried out using small angle light scattering to determine the fibre level response due to strain induced degradation. Next, strain dependent degradation rates were determined from stress relaxation experiments in the presence of crude and purified collagenase to determine the tissue level degradation response. Finally, a 1D theoretical model was developed, incorporating matrix stiffness and a gradient of collagen fibre crimp to decouple the mechanism behind strain dependent arterial degradation. SALS structural analysis identified a strain mediated degradation response in arterial tissue at the fibre level not dissimilar to that found in literature for pure collagen. Interestingly, two distinctly different strain mediated degradation responses were identified experimentally at the tissue level, not seen in other collagenous tissues. Our model was able to accurately predict these experimental findings, but only once the load bearing matrix, its degradation response and the gradient of collagen fibre crimp across the arterial wall were incorporated. These findings highlight the critical role that the various tissue constituents play in the degradation response of arterial tissue. Statement of Significance Collagen fibre architecture is the dominant load bearing component of arterial tissue. Remodelling of this architecture is a strain dependent process stimulated by the degradation of existing collagen. Despite this, degradation of arterial tissue and in particular, arterial collagen, is not fully understood or studied. In the current study, we identified for the first time, the strain dependent degradation response of arterial tissue, which has not been observed in other collagenous tissues in literature. We hypothesised that this unique degradation response was due to the complex structure observed in arterial tissue. Based on this hypothesis, we developed a novel numerical model capable of explaining this unique degradation response which may provide critical insights into disease development and aid in the design of interventional medical devices.

Contractile stress generation by adherent cells is largely determined by the interplay of forces within their cytoskeleton. It is known that actin stress fibers, connected to focal adhesions, provide contractile stress generation, while microtubules and intermediate filaments provide cells compressive stiffness. Recent studies have shown the importance of the interplay between the stress fibers and the intermediate filament vimentin. Therefore, the effect of the interplay between the stress fibers and vimentin on stress generation was quantified in this study. We hypothesized that net stress generation comprises the stress fiber contraction combined with the vimentin resistance. We expected an increased net stress in vimentin knockout (VimKO) mouse embryonic fibroblasts (MEFs) compared to their wild-type (vimentin wild-type (VimWT)) counterparts, due to the decreased resistance against stress fiber contractility. To test this, the net stress generation by VimKO and VimWT MEFs was determined using the thin film method combined with sample-specific finite element modeling. Additionally, focal adhesion and stress fiber organization were examined via immunofluorescent staining. Net stress generation of VimKO MEFs was three-fold higher compared to VimWT MEFs. No differences in focal adhesion size or stress fiber organization and orientation were found between the two cell types. This suggests that the increased net stress generation in VimKO MEFs was caused by the absence of the resistance that vimentin provides against stress fiber contraction. Taken together, these data suggest that vimentin resists the stress fiber contractility, as hypothesized, thus indicating the importance of vimentin in regulating cellular stress generation by adherent cells.

Hemodynamic forces and Notch signaling are both known as key regulators of arterial remodeling and homeostasis. However, how these two factors integrate in vascular morphogenesis and homeostasis is unclear. Here, we combined experiments and modeling to evaluate the impact of the integration of mechanics and Notch signaling on vascular homeostasis. Vascular smooth muscle cells (VSMCs) were cyclically stretched on flexible membranes, as quantified via video tracking, demonstrating that the expression of Jagged1, Notch3, and target genes was down-regulated with strain. The data were incorporated in a computational framework of Notch signaling in the vascular wall, where the mechanical load was defined by the vascular geometry and blood pressure. Upon increasing wall thickness, the model predicted a switch-type behavior of the Notch signaling state with a steep transition of synthetic toward contractile VSMCs at a certain transition thickness. These thicknesses varied per investigated arterial location and were in good agreement with human anatomical data, thereby suggesting that the Notch response to hemodynamics plays an important role in the establishment of vascular homeostasis., Proceedings of the National Academy of Sciences of the United States of America, 115 (16), ISSN:0027-8424, ISSN:1091-6490

Recent studies on decellularized tissue engineered heart valves (DTEHVs) showed rapid host cell repopulation and increased valvular insufficiency developing over time, associated with leaflet shortening. A possible explanation for this result was found using computational simulations, which revealed radial leaflet compression in the original valvular geometry when subjected to physiological pressure conditions. Therefore, an improved geometry was suggested to enable radial leaflet extension to counteract for host cell mediated retraction. In this study, we propose a solution to impose this new geometry by using a constraining bioreactor insert during culture. Human cell based DTEHVs (n = 5) were produced as such, resulting in an enlarged coaptation area and profound belly curvature. Extracellular matrix was homogeneously distributed, with circumferential collagen alignment in the coaptation region and global tissue anisotropy. Based on in vitro functionality experiments, these DTEHVs showed competent hydrodynamic functionality under physiological pulmonary conditions and were fatigue resistant, with stable functionality up to 16 weeks in vivo simulation. Based on implemented mechanical data, our computational models revealed a considerable decrease in radial tissue compression with the obtained geometrical adjustments. Therefore, these improved DTEHV are expected to be less prone to host cell mediated leaflet retraction and will remain competent after implantation.

One of the hallmarks of biological soft tissues is their capacity to grow and remodel in response to changes in their environment. Although it is well-accepted that these processes occur at least partly to maintain a mechanical homeostasis, it remains unclear which mechanical constituent(s) determine(s) mechanical homeostasis. In the current study a nondestructive mechanical test and a two-step inverse analysis method were developed and validated to nondestructively estimate the mechanical properties of biological tissue during tissue culture. Nondestructive mechanical testing was achieved by performing an inflation test on tissues that were cultured inside a bioreactor, while the tissue displacement and thickness were nondestructively measured using ultrasound. The material parameters were estimated by an inverse finite element scheme, which was preceded by an analytical estimation step to rapidly obtain an initial estimate that already approximated the final solution. The efficiency and accuracy of the two-step inverse method was demonstrated on virtual experiments of several material types with known parameters. PDMS samples were used to demonstrate the method's feasibility, where it was shown that the proposed method yielded similar results to tensile testing. Finally, the method was applied to estimate the material properties of tissue-engineered constructs. Via this method, the evolution of mechanical properties during tissue growth and remodeling can now be monitored in a well-controlled system. The outcomes can be used to determine various mechanical constituents and to assess their contribution to mechanical homeostasis.

Tissue growth and remodeling are essential processes that should ensure long-term functionality of tissue-engineered (TE) constructs. Even though it is widely recognized that these processes strongly depend on mechanical stimuli, the underlying mechanisms of mechanically induced growth and remodeling are only partially understood. It is generally accepted that cells sense mechanical changes and respond by altering their surroundings, by means of extracellular matrix growth and remodeling, in an attempt to return to a certain preferred mechanical homeostatic state. However, the exact mechanical cues that trigger cells to synthesize and remodel their environment remain unclear. To identify the driving mechanical stimuli of these processes, it is critical to be able to temporarily follow the mechanical state of developing tissues under physiological loading conditions. Therefore, a novel "versatile tissue growth and remodeling" (Vertigro) bioreactor was developed that is capable of tissue culture and mechanical stimulation for a prolonged time period, while simultaneously performing mechanical testing. The Vertigro's unique two-chamber design allows easy, sterile handling of circular 3D TE constructs in a dedicated culture chamber, while a separate pressure chamber facilitates a pressure-driven dynamic loading regime during culture. As a proof-of-concept, temporal changes in the mechanical state of cultured tissues were quantified using nondestructive mechanical testing by means of a classical bulge test, in which the tissue displacement was tracked using ultrasound imaging. To demonstrate the successful development of the bioreactor system, compositional, structural, and geometrical changes were qualitatively and quantitatively assessed using a series of standard analysis techniques. With this bioreactor and associated mechanical analysis technique, a powerful toolbox has been developed to quantitatively study and identify the driving mechanical stimuli of engineered tissue growth and remodeling.

Valvular heart disease is a major cause of morbidity and mortality worldwide. Current heart valve prostheses have considerable clinical limitations due to their artificial, nonliving nature without regenerative capacity. To overcome these limitations, heart valve tissue engineering (TE) aiming to develop living, native-like heart valves with self-repair, remodeling, and regeneration capacity has been suggested as next-generation technology. A major roadblock to clinically relevant, safe, and robust TE solutions has been the high complexity and variability inherent to bioengineering approaches that rely on cell-driven tissue remodeling. For heart valve TE, this has limited long-term performance in vivo because of uncontrolled tissue remodeling phenomena, such as valve leaflet shortening, which often translates into valve failure regardless of the bioengineering methodology used to develop the implant. We tested the hypothesis that integration of a computationally inspired heart valve design into our TE methodologies could guide tissue remodeling toward long-term functionality in tissue-engineered heart valves (TEHVs). In a clinically and regulatory relevant sheep model, TEHVs implanted as pulmonary valve replacements using minimally invasive techniques were monitored for 1 year via multimodal in vivo imaging and comprehensive tissue remodeling assessments. TEHVs exhibited good preserved long-term in vivo performance and remodeling comparable to native heart valves, as predicted by and consistent with computational modeling. TEHV failure could be predicted for nonphysiological pressure loading. Beyond previous studies, this work suggests the relevance of an integrated in silico, in vitro, and in vivo bioengineering approach as a basis for the safe and efficient clinical translation of TEHVs.

Tendinopathy is a widespread and unresolved clinical challenge, in which associated pain and hampered mobility present a major cause for work-related disability. Tendinopathy associates with a change from a healthy tissue with aligned extracellular matrix (ECM) and highly polarized cells that are connected head-to-tail, towards a diseased tissue with a disorganized ECM and randomly distributed cells, scar-like features that are commonly attributed to poor innate regenerative capacity of the tissue. A fundamental clinical dilemma with this scarring process is whether treatment strategies should focus on healing the affected (disorganized) tissue or strengthen the remaining healthy (anisotropic) tissue. The question was thus asked whether the intrinsic remodeling capacity of tendon-derived cells depends on the organization of the 3D extracellular matrix (isotropic vs anisotropic). Progress in this field is hampered by the lack of suitable in vitro tissue platforms. We aimed at filling this critical gap by creating and exploiting a next generation tissue platform that mimics aspects of the tendon scarring process; cellular response to a gradient in tissue organization from isotropic (scarred/non-aligned) to highly anisotropic (unscarred/aligned) was studied, as was a transient change from isotropic towards highly anisotropic. Strikingly, cells residing in an ‘unscarred’ anisotropic tissue indicated superior remodeling capacity (increased gene expression levels of collagen, matrix metalloproteinases MMPs, tissue inhibitors of MMPs), when compared to their ‘scarred’ isotropic counterparts. A numerical model then supported the hypothesis that cellular remodeling capacity may correlate to cellular alignment strength. This in turn may have improved cellular communication, and could thus relate to the more pronounced connexin43 gap junctions observed in anisotropic tissues. In conclusion, increased tissue anisotropy was observed to enhance the cellular potential for functional remodeling of the matrix. This may explain the poor regenerative capacity of tenocytes in chronic tendinopathy, where the pathological process has resulted in ECM disorganization. Additionally, it lends support to treatment strategies that focus on strengthening the remaining healthy tissue, rather than regenerating scarred tissue., Matrix Biology, 65, ISSN:0945-053X, ISSN:1569-1802

Collagen is the main load-bearing component of many soft tissues and has a large influence on the mechanical behavior of tissues when exposed to mechanical loading. Therefore, it is important to increase our understanding of collagen remodeling in soft tissues to understand the mechanisms behind pathologies and to control the development of the collagen network in engineered tissues. In the present study, a constitutive model was developed by coupling a recently developed model describing the orientation and contractile stresses exerted by cells in response to mechanical stimuli to physically motivated collagen remodeling laws. In addition, cell-mediated contraction of the collagen fibers was included as a mechanism for tissue compaction. The model appeared to be successful in predicting a range of experimental observations, which are (1) the change in transition stretch of periosteum after remodeling at different applied stretches, (2) the compaction and alignment of collagen fibers in tissue-engineered strips, (3) the fiber alignment in cruciform gels with different arm widths, and (4) the alignment of collagen fibers in engineered vascular grafts. Moreover, by changing the boundary conditions, the model was able to predict a helical architecture in the vascular graft without assuming the presence of two helical fiber families a priori. Ultimately, this model may help to increase our understanding of collagen remodeling in physiological and pathological conditions, and it may provide a tool for determining the optimal experimental conditions for obtaining native-like collagen architectures in engineered tissues.

To build anisotropic, mechanically functioning tissue, it is essential to understand how cells orient in response to mechanical stimuli. Therefore, a computational model was developed which predicts cell orientation, based on the actin stress fiber distribution inside the cell. In the model, the stress fiber distribution evolves dynamically according to the following: (1) Stress fibers contain polymerized actin. The total amount of depolymerized plus polymerized actin is constant. (2) Stress fibers apply tension to their environment. This active tension is maximal when strain rate and absolute strain are zero and reduces with increasing shortening rate and absolute strain. (3) A high active fiber stress in a direction leads to a large amount of fibers in this direction. (4) The cell is attached to a substrate; all fiber stresses are homogenized into a total cell stress, which is in equilibrium with substrate stress. This model predicts that on a substrate of anisotropic stiffness, fibers align in the stiffest direction. Under cyclic strain when the cellular environment is so stiff that no compaction occurs (1 MPa), the model predicts strain avoidance, which is more pronounced with increasing strain frequency or amplitude. Under cyclic strain when the cellular environment is so soft that cells can compact it (10 kPa), the model predicts a preference for the cyclically strained compared to the compacting direction. These model predictions all agree with experimental evidence. For the first time, a computational model predicts cell orientation in response to this range of mechanical stimuli using a single set of parameters. © 2013 Springer-Verlag Berlin Heidelberg.

Deep tissue injury (DTI) is a severe form of pressure ulcer where tissue damage starts in deep tissues underneath intact skin. Tissue deformation may play an important role in the aetiology, which can be investigated using an experimental–numerical approach. Recently, an animal-specific finite element model has been developed to simulate experiments in which muscle tissue was compressed with an indenter. In this study, the material behaviour and boundary conditions were adapted to improve the agreement between model and experiment and to investigate the influence of these adaptations on the predicted strain distribution. The use of a highly nonlinear material law and including friction between the indenter and the muscle both improved the quality of the model and considerably influenced the estimated strain distribution. With the improved model, the required sample size to detect significant differences between loading conditions can be diminished, which is clearly relevant in experiments involving animals.

Pressure ulcers are localized areas of soft tissue breakdown due to mechanical loading. Susceptible individuals are subjected to pressure relief strategies to prevent long loading periods. Therefore, ischemia-reperfusion injury may play an important role in the etiology of pressure ulcers. To investigate the inter-relation between postischemic perfusion and changes in skeletal muscle integrity, the hindlimbs of Brown Norway rats were subjected to 4-h ischemia followed by 2-h reperfusion. Dynamic contrast-enhanced MRI was used to examine perfusion, and changes in skeletal muscle integrity were monitored with T2-weighted MRI. The dynamic contrast-enhanced MRI data showed a heterogeneous postischemic profile in the hindlimb, consisting of areas with increased contrast enhancement (14-76% of the hindlimb) and regions with no-reflow (5-77%). For T2, a gradual increase in the complete leg was observed during the 4-h ischemic period (from 34 to 41 msec). During the reperfusion phase, a heterogeneous distribution of T2 was observed. Areas with increased contrast enhancement were associated with a decrease in T2 (to 38 msec) toward preischemic levels, whereas no-reflow areas exhibited a further increase in T2 (to 42 msec). These results show that reperfusion after prolonged ischemia may not be complete, thereby continuing the ischemic condition and aggravating tissue damage.

In situ cardiovascular tissue-engineering can potentially address the shortcomings of the current replacement therapies, in particular, their inability to grow and remodel. In native tissues, it is widely accepted that physiological growth and remodelling occur to maintain a homeostatic mechanical state to conserve its function, regardless of changes in the mechanical environment. A similar homeostatic state should be reached for tissue-engineered (TE) prostheses to ensure proper functioning. For in situ tissue-engineering approaches obtaining such a state greatly relies on the initial scaffold design parameters. In this study, it is investigated if the simple scaffold design parameter initial thickness, influences the emergence of a mechanical and geometrical equilibrium state in in vitro TE constructs, which resemble thin cardiovascular tissues such as heart valves and arteries. Towards this end, two sample groups with different initial thicknesses of myofibroblast-seeded polycaprolactone-bisurea constructs were cultured for three weeks under dynamic loading conditions, while tracking geometrical and mechanical changes temporally using non-destructive ultrasound imaging. A mechanical equilibrium was reached in both groups, although at different magnitudes of the investigated mechanical quantities. Interestingly, a geometrically stable state was only established in the thicker constructs, while the thinner constructs’ length continuously increased. This demonstrates that reaching geometrical and mechanical stability in TE constructs is highly dependent on functional scaffold design.

For pressure ulcer prevention an ambitious goal would be the establishment of a mechanical threshold for tissue damage. In the past, several researchers have sought to establish such a threshold often involving the loading time. However, they have not resulted in a unique reliable value that could be used in practice. This limitation is probably due to the focus on interface pressure. The objective of this paper is to clarify to an audience with no conventional background in mechanics, why interface pressure is not the appropriate parameter to define a damage threshold, whereas internal local deformations (strains) may prove more suitable. The paper reveals that it may be possible to identify a damage threshold for healthy skeletal muscle tissue based on local internal deformations.

Sustained mechanical loading of skeletal muscle may result in the development of a severe type of pressure ulcer, referred to as deep tissue injury. Recently it was shown that the diffusion of large molecules (10-150 kDa) is impaired during deformation of tissue-engineered skeletal muscle, suggesting a role for impaired diffusion in the aetiology of deep tissue injury. However, the influence of deformation on diffusion of smaller molecules on its aetiology is less clear. This motivated the present study designed to investigate the influence of deformation of skeletal muscle on the diffusion of water, which can be measured with diffusion tensor magnetic resonance imaging (MRI). It could be predicted that this approach will provide valuable information on the diffusion of small molecules. Additionally the relationship between muscle temperature and diffusion was investigated. During deformation of the tibialis anterior a decrease of the apparent diffusion coefficient (ADC) was observed (7.2±3.9%). The use of a finite element model showed that no correlation existed between the maximum shear strain and the decrease of the ADC. The ADC in the uncompressed gastrocnemius muscle decreased with 5.9±3.7%. In an additional experiment a clear correlation was obtained between the decrease of the ADC and the relative temperature change of skeletal muscle tissue as measured by MRI. Taken together, it was concluded that (1) the decreased diffusion of water was not a direct effect of tissue deformation and (2) that it is likely that the observed decreased ADC during deformation was a result of a decreased muscle temperature. The present study therefore provides evidence that diffusion of small molecules, particularly oxygen and carbon dioxide, is not impaired during deformation of skeletal muscle tissue. © 2009 Elsevier Ltd. All rights reserved.

Continuous flow ventricular assist devices (VADs) for mechanical circulatory support (MCS) are generally smaller and believed to be more reliable than pulsatile VADs. However, regarding continuous flow, there are concerns about the decreased pulsatility and ventricular unloading. Moreover, pulsatile VADs offer a wider range in control strategies. For this reason, we used a computer model to evaluate whether pulsatile operation of a continuous flow VAD would be more beneficial than the standard constant pump speed. The computer model describes the left and right ventricle with one-fiber heart contraction models, and the systemic, pulmonary, and coronary circulation with lumped parameter hemodynamical models, while the heart rate is regulated with a baroreflex model. With this computer model, both normal and heart failure hemodynamics were simulated. A HeartMate II left ventricular assist device model was connected to this model, and both constant speed and pulsatile support were simulated. Pulsatile support did not solve the decreased pulsatility issue, but it did improve perfusion (cardiac index and coronary flow) and unloading (stroke work and heart rate) compared with constant speed. Also, pulsatile support would be beneficial for developing control strategies, as it offers more options to adjust assist device settings to the patient's needs. Because the mathematical model used in this study can simulate different assist device settings, it can play a valuable role in developing mechanical circulatory support control strategies. © 2009, International Center for Artificial Organs and Transplantation and Wiley Periodicals, Inc.

The paper describes the current views on the cause of a sub-class of pressure ulcers known as pressure induced deep tissue injury (DTI). A multi-scale approach was adopted using model systems ranging from single cells in culture, tissue engineered muscle to animal studies with small animals. This has led to a clear understanding on two damage mechanisms associated with the development of DTI. Direct deformation results from high, but physiologically relevant, strains and is a process that leads to the first signs of cell damage within minutes. Ischaemic damage is caused by occlusion of blood vessels, but takes several hours to develop. The paper ends with some clinical consequences.

In many tissue engineering applications it is essential to understand how cells orient under the influence of their mechanical environment. In vitro engineered models are used to investigate the orientation of F-actin stress fibers inside cells. One such in vitro model [1] consists of a mixture of cells, collagen and matrigel, that is constrained by an array of silicone posts (Figure 1). We have recently developed a computational model to describe the orientation of stress fibers in response to their mechanical environment [2]. In the present study, this computational model is extended to 3D and used to simulate cell behavior in the mentioned in vitro model. This improves our understanding of how stress fibers orient in response to the mechanical environment and aids in optimizing the use of the in vitro model.Copyright © 2013 by ASME

Tissue engineering represents a promising technique to overcome the limitations of the current valve replacements, since it allows for creating living autologous heart valves that have the potential to grow and remodel. However, also this approach still faces a number of challenges. One particular problem is regurgitation, caused by cell-mediated tissue retraction or the mismatch in geometrical and material properties between tissue-engineered heart valves (TEHVs) and their native counterparts. The goal of the present study was to assess the influence of valve geometry and tissue anisotropy on the deformation profile and closed configuration of TEHVs. To achieve this aim, a range of finite element models incorporating different valve shapes was developed, and the constitutive behavior of the tissue was modeled using an established computational framework, where the degree of anisotropy was varied between values representative of TEHVs and native valves. The results of this study suggest that valve geometry and tissue anisotropy are both important to maximize the radial strains and thereby the coaptation area. Additionally, the minimum degree of anisotropy that is required to obtain positive radial strains was shown to depend on the valve shape and the pressure to which the valves are exposed. Exposure to pulmonary diastolic pressure only yielded positive radial strains if the anisotropy was comparable to the native situation, whereas considerably less anisotropy was required if the valves were exposed to aortic diastolic pressure.

Deep pressure ulcers represent a major problem for individuals with spinal cord injury (SCI), with the initial damage often hidden underneath intact skin. Accordingly, early detection is difficult and treatment is problematic. In the present study, circulatory levels of biomarkers for muscle damage were investigated to explore their potential in the early detection of deep pressure ulcers. Baseline concentrations of creatine kinase, myoglobin (Mb), heart-type fatty acid binding protein (H-FABP), and C-reactive protein (CRP) were measured in small groups of nondisabled (age 39-66 yr) subjects and subjects with SCI (age 40-68 yr, American Spinal Injury Association grade A B, level of injury thoracic 11 to lumbar 3) over a period of 5 days. Each subject exhibited a unique concentration profile fer all markers, although some correlations were observed; for example, Mb and H-FABP were correlated for both subject groups. No significant differences were found in marker concentrations between the two subject groups, although a trend toward higher CRP levels was observed in the SCI subjects. Furthermore, one SCI subject with a category II pressure ulcer exhibited higher H-FABP and CRP concentrations than all other subjects. Because the variations in each of the marker concentrations were smaller than the predicted increases after pressure ulcers, this combination of plasma markers may prove appropriate for the early detection of deep pressure ulcers.

Sustained mechanical loading of soft tissues covering bony prominences may lead to degeneration of skeletal muscle tissue. This can result in a condition termed deep tissue injury (DTI), a severe kind of pressure ulcer that initiates in deep tissue layers, and progresses towards the skin. Previously, we have provided evidence that in a controlled animal model, deformation is the main trigger for damage within a 2 h loading period [1,2]. Recently, we also showed that ischemia and reperfusion may contribute to the damage process during prolonged loading [3]. In the present study, we investigated the relative effects of deformation, ischemia, and reperfusion on the temporal and spatial damage process of skeletal muscle tissue during a 6 h period using magnetic resonance imaging (MRI) techniques.

Mechanical loading of soft tissues covering bony prominences can cause skeletal muscle damage, ultimately resulting in a severe pressure ulcer termed deep tissue injury (DTI). Deformation plays an important role in the aetiology of DTI. Therefore, it is essential to minimise internal muscle deformations in subjects at risk of DTI. As an example, spinal cord-injured (SCI) individuals exhibit structural changes leading to a decrease in muscle thickness and stiffness, which subsequently increase the tissue deformations. In the present study, an animal-specific finite element model, where the geometry and boundary conditions were derived from magnetic resonance images, was developed. It was used to investigate the internal deformations in the muscle, fat and skin layers of the porcine buttocks during loading. The model indicated the presence of large deformations in both the muscle and the fat layers, with maximum shear strains up to 0.65 in muscle tissue and 0.63 in fat. Furthermore, a sensitivity analysis showed that the tissue deformations depend considerably on the relative stiffness values of the different tissues. For example, a change in muscle stiffness had a large effect on the muscle deformations. A 50% decrease in stiffness caused an increase in maximum shear strain from 0.65 to 0.99, whereas a 50% increase in stiffness resulted in a decrease in maximum shear strain from 0.65 to 0.49. These results indicate the importance of restoring tissue properties after SCI, with the use of, for example, electrical stimulation, to prevent the development of DTI.

Deep tissue injury (DTI) is a severe form of pressure ulcer where tissue damage starts in deep tissues underneath intact skin. In the present study, the contributions of deformation, ischemia, and reperfusion to skeletal muscle damage development were examined in a rat model during a 6-h period. Magnetic resonance imaging (MRI) was used to study perfusion (contrast-enhanced MRI) and tissue integrity (T2-weighted MRI). The levels of tissue deformation were estimated using finite element models. Complete ischemia caused a gradual homogeneous increase in T2 (∼20% during the 6-h period). The effect of reperfusion on T2 was highly variable, depending on the anatomical location. In experiments involving deformation, inevitably associated with partial ischemia, a variable T2 increase (17–66% during the 6-h period) was observed reflecting the significant variation in deformation (with two-dimensional strain energies of 0.60–1.51 J/mm) and ischemia (50.8–99.8% of the leg) between experiments. These results imply that deformation, ischemia, and reperfusion all contribute to the damage process during prolonged loading, although their importance varies with time. The critical deformation threshold and period of ischemia that cause muscle damage will certainly vary between individuals. These variations are related to intrinsic factors, such as pathological state, which partly explain the individual susceptibility to the development of DTI and highlight the need for regular assessments of individual subjects.