A Spontaneous In Situ Thiol-Ene Crosslinking Hydrogel with Thermo-Responsive Mechanical Properties

The thermo-responsive behavior of Poly(N-isopropylacrylamide) makes it an ideal candidate to easily embed cells and allows the polymer mixture to be injected. However, P(NiPAAm) hydrogels possess minor mechanical properties. To increase the mechanical properties, a covalent bond is introduced into the P(NIPAAm) network through a biocompatible thiol-ene click-reaction by mixing two polymer solutions. Co-polymers with variable thiol or acrylate groups to thermo-responsive co-monomer ratios, ranging from 1% to 10%, were synthesized. Precise control of the crosslink density allowed customization of the hydrogel's mechanical properties to match different tissue stiffness levels. Increasing the temperature of the hydrogel above its transition temperature of 31 °C induced the formation of additional physical interactions. These additional interactions both further increased the stiffness of the material and impacted its relaxation behavior. The developed optimized hydrogels reach stiffnesses more than ten times higher compared to the state of the art using similar polymers. Furthermore, when adding cells to the precursor polymer solutions, homogeneous thermo-responsive hydrogels with good cell viability were created upon mixing. In future work, the influence of the mechanical micro-environment on the cell's behavior can be studied in vitro in a continuous manner by changing the incubation temperature.

SGABU computational platform for multiscale modeling: Bridging the gap between education and research

BACKGROUND AND OBJECTIVE

In accordance with the latest aspirations in the field of bioengineering, there is a need to create a web accessible, but powerful cloud computational platform that combines datasets and multiscale models related to bone modeling, cancer, cardiovascular diseases and tissue engineering. The SGABU platform may become a powerful information system for research and education that can integrate data, extract information, and facilitate knowledge exchange with the goal of creating and developing appropriate computing pipelines to provide accurate and comprehensive biological information from the molecular to organ level.

METHODS

The datasets integrated into the platform are obtained from experimental and/or clinical studies and are mainly in tabular or image file format, including metadata. The implementation of multiscale models, is an ambitious effort of the platform to capture phenomena at different length scales, described using partial and ordinary differential equations, which are solved numerically on complex geometries with the use of the finite element method. The majority of the SGABU platform's simulation pipelines are provided as Common Workflow Language (CWL) workflows. Each of them requires creating a CWL implementation on the backend and a user-friendly interface using standard web technologies. Platform is available at https://sgabu-test.unic.kg.ac.rs/login.

RESULTS

The main dashboard of the SGABU platform is divided into sections for each field of research, each one of which includes a subsection of datasets and multiscale models. The datasets can be presented in a simple form as tabular data, or using technologies such as Plotly.js for 2D plot interactivity, Kitware Paraview Glance for 3D view. Regarding the models, the usage of Docker containerization for packing the individual tools and CWL orchestration for describing inputs with validation forms and outputs with tabular views for output visualization, interactive diagrams, 3D views and animations.

CONCLUSIONS

In practice, the structure of SGABU platform means that any of the integrated workflows can work equally well on any other bioengineering platform. The key advantage of the SGABU platform over similar efforts is its versatility offered with the use of modern, modular, and extensible technology for various levels of architecture.

Synthetic fibrous hydrogels as a platform to decipher cell-matrix mechanical interactions

Cells continuously sense external forces from their microenvironment, the extracellular matrix (ECM). In turn, they generate contractile forces, which stiffen and remodel this matrix. Although this bidirectional mechanical exchange is crucial for many cell functions, it remains poorly understood. Key challenges are that the majority of available matrices for such studies, either natural or synthetic, are difficult to control or lack biological relevance. Here, we use a synthetic, yet highly biomimetic hydrogel based on polyisocyanide (PIC) polymers to investigate the effects of the fibrous architecture and the nonlinear mechanics on cell-matrix interactions. Live-cell rheology was combined with advanced microscopy-based approaches to understand the mechanisms behind cell-induced matrix stiffening and plastic remodeling. We demonstrate how cell-mediated fiber remodeling and the propagation of fiber displacements are modulated by adjusting the biological and mechanical properties of this material. Moreover, we validate the biological relevance of our results by demonstrating that cellular tractions in PIC gels develop analogously to those in the natural ECM. This study highlights the potential of PIC gels to disentangle complex bidirectional cell-matrix interactions and to improve the design of materials for mechanobiology studies.

Mechanical properties-translucency-microstructure relationships in commercial monolayer and multilayer monolithic zirconia ceramics

OBJECTIVES

To evaluate the phase composition, microstructure, optical properties and mechanical properties of eight commercially available multilayer and monolayer monolithic dental zirconias.

METHODS

Five commercial 3Y-TZP (GC ST, GC HT [GC, Tokyo Japan]; Katana ML, Katana HT [Kuraray Noritake] and Lava Plus [3M Oral Care]) and three Y-PSZ (Katana STML, Katana UTML [Kuraray Noritake]; GC UHT [GC, Tokyo Japan]) zirconia ceramic grades were cut in plate-shaped specimens, sintered according to the manufacturer's instructions and mirror polished. The zirconia chemical composition was determined using X-ray fluorescence (XRF), phase composition was characterized using X-ray diffraction (XRD), while the grain size was measured using scanning electron microscopy (SEM). The translucency Parameter (TP) and Contrast Ratio (CR) were measured with a spectrophotometer (n = 10/group). The indentation fracture toughness (n = 10), Vickers hardness (n = 10) and biaxial strength (n = 20) of the sintered ceramics were assessed. The stress distribution during biaxial testing was assessed by Finite element analysis (FEA). Statistical analysis involved one-way ANOVA and post-hoc Tukey's HSD test and Pearson correlation test (α = 0.05).

RESULTS

FEA showed that the stress distribution in plate shape specimens was the same as for disks, rationalizing the use of plates for biaxial strength testing. As expected, higher quantities of Y₂O₃ were related to a higher cubic ZrO₂ phase content and lower tetragonality t-ZrO₂, which improved translucency but diminished flexural strength and toughness. While there was no significant correlation between grain size and other material properties, addition of pigments to the zirconia grade statistically negatively affected hardness.

CONCLUSION

Even though an improvement in strength and translucency could be recorded for the last Y-TZP generation, future research still needs to strive for combined improvement of optical properties and mechanical reliability of zirconia ceramics.

Patterned dextran ester films as a tailorable cell culture platform

The elucidation of cell-surface interactions and the development of model platforms to help uncover their underlying mechanisms remains vital to the design of effective biomaterials. To this end, dextran palmitates with varying degrees of substitution were synthesised with a multipurpose functionality: an ability to modulate surface energy through surface chemistry, and an ideal thermal behaviour for patterning. Herein, dextran palmitate films are produced by spin coating, and patterned by thermal nanoimprint lithography with nano-to-microscale topographies. These films of moderately hydrophobic polysaccharide esters with low nanoscale roughness performed as well as fibronectin coatings in the culture of bovine aortic endothelial cells. Upon patterning, they display distinct regions of roughness, restricting cell adhesion to the smoothest surfaces, while guiding multicellular arrangements in the patterned topographies. The development of biomaterial interfaces through topochemical fabrication such as this could prove useful in understanding protein and cell-surface interactions.

Fibrodysplasia Ossificans Progressiva: What Have We Achieved and Where Are We Now? Follow-up to the 2015 Lorentz Workshop

Fibrodysplasia ossificans progressiva (FOP) is an ultra-rare progressive genetic disease effecting one in a million individuals. During their life, patients with FOP progressively develop bone in the soft tissues resulting in increasing immobility and early death. A mutation in the ACVR1 gene was identified as the causative mutation of FOP in 2006. After this, the pathophysiology of FOP has been further elucidated through the efforts of research groups worldwide. In 2015, a workshop was held to gather these groups and discuss the new challenges in FOP research. Here we present an overview and update on these topics.

Simulating flow induced migration in vascular remodelling

Shear stress induces directed endothelial cell (EC) migration in blood vessels leading to vessel diameter increase and induction of vascular maturation. Other factors, such as EC elongation and interaction between ECs and non-vascular areas are also important. Computational models have previously been used to study collective cell migration. These models can be used to predict EC migration and its effect on vascular remodelling during embryogenesis. We combined live time-lapse imaging of the remodelling vasculature of the quail embryo yolk sac with flow quantification using a combination of micro-Particle Image Velocimetry and computational fluid dynamics. We then used the flow and remodelling data to inform a model of EC migration during remodelling. To obtain the relation between shear stress and velocity in vitro for EC cells, we developed a flow chamber to assess how confluent sheets of ECs migrate in response to shear stress. Using these data as an input, we developed a multiphase, self-propelled particles (SPP) model where individual agents are driven to migrate based on the level of shear stress while maintaining appropriate spatial relationship to nearby agents. These agents elongate, interact with each other, and with avascular agents at each time-step of the model. We compared predicted vascular shape to real vascular shape after 4 hours from our time-lapse movies and performed sensitivity analysis on the various model parameters. Our model shows that shear stress has the largest effect on the remodelling process. Importantly, however, elongation played an especially important part in remodelling. This model provides a powerful tool to study the input of different biological processes on remodelling.

Shear stress is known to play a leading role in endothelial cell (EC) migration and hence, vascular remodelling. Vascular remodelling is, however, more complicated than only EC migration. To achieve a better understanding of this process, we developed a computational model in which, shear stress mediated EC migration has the leading role and other factors, such as avascular regions and EC elongation, are also accounted for. We have tested this model for different vessel shapes during remodelling and could study the role that each of these factors play in remodelling. This model gives us the possibility of addition of other factors such as biochemical signals and angiogenesis which will help us in the study of vascular remodelling in both development and disease.

Actomyosin-dependent invasion of endothelial sprouts in collagen

During sprouting angiogenesis-the growth of blood vessels from the existing vasculature-endothelial cells (ECs) adopt an elongated invasive form and exert forces at cell-cell and cell-matrix interaction sites. These cell shape changes and cellular tractions require extensive reorganizations of the actomyosin network. However, the respective roles of actin and myosin for endothelial sprouting are not fully elucidated. In this study, we further investigate these roles by treating 2D-migrating and 3D-sprouting ECs with chemical compounds targeting either myosin or actin. These treatments affected the endothelial cytoskeleton drastically and reduced the invasive response in a compound-specific manner; pointing toward a tight control of the actin and myosin activity during sprouting. Clusters in the data further illustrate that endothelial sprout morphology is sensitive to the in vitro model mechanical microenvironment and directs future research toward mechanical substrate guidance as a strategy for promoting engineered tissue vascularization. In summary, our results add to a growing corpus of research highlighting a key role of the cytoskeleton for sprouting angiogenesis.

Intercellular Adhesion Stiffness Moderates Cell Decoupling as a Function of Substrate Stiffness

The interplay between cell-cell and cell-substrate interactions is complex yet necessary for the formation and healthy functioning of tissues. The same mechanosensing mechanisms used by the cell to sense its extracellular matrix also play a role in intercellular interactions. We used the discrete element method to develop a computational model of a deformable cell that includes subcellular components responsible for mechanosensing. We modeled a three-dimensional cell pair on a patterned (two-dimensional) substrate, a simple laboratory setup to study intercellular interactions. We explicitly modeled focal adhesions and adherens junctions. These mechanosensing adhesions matured, becoming stabilized by force. We also modeled contractile stress fibers that bind the discrete adhesions. The mechanosensing fibers strengthened upon stalling. Traction exerted on the substrate was used to generate traction maps (along the cell-substrate interface). These simulated maps are compared to experimental maps obtained via traction force microscopy. The model recreates the dependence on substrate stiffness of the tractions' spatial distribution, contractile moment of the cell pair, intercellular force, and number of focal adhesions. It also recreates the phenomenon of cell decoupling, in which cells exert forces separately when substrate stiffness increases. More importantly, the model provides viable molecular explanations for decoupling: mechanosensing mechanisms are responsible for competition between different fiber-adhesion configurations present in the cell pair. The point at which an increasing substrate stiffness becomes as high as that of the cell-cell interface is the tipping point at which configurations that favor cell-substrate adhesion dominate over those favoring cell-cell adhesion. This competition is responsible for decoupling.

Modeling of Mechanosensing Mechanisms Reveals Distinct Cell Migration Modes to Emerge From Combinations of Substrate Stiffness and Adhesion Receptor-Ligand Affinity

Mesenchymal cell migration is an integral process in development and healing. The process is regulated by both mechanical and biochemical properties. Mechanical properties of the environment are sensed through mechanosensing, which consists of molecular responses mediated by mechanical signals. We developed a computational model of a deformable 3D cell on a flat substrate using discrete element modeling. The cell is polarized in a single direction and thus moves along the long axis of the substrate. By modeling discrete focal adhesions and stress fibers, we implement two mechanosensing mechanisms: focal adhesion stabilization by force and stress fiber strengthening upon contraction stalling. Two substrate-associated properties, substrate (ligand) stiffness and adhesion receptor–ligand affinity (in the form of focal adhesion disassembly rate), were varied for different model setups in which the mechanosensing mechanisms are set as active or inactive. Cell displacement, focal adhesion number, and cellular traction were quantified and tracked in time. We found that varying substrate stiffness (a mechanical property) and adhesion receptor–ligand affinity (a biochemical property) simultaneously dictate the mode in which cells migrate; cells either move in a smooth manner reminiscent of keratocytes or in a cyclical manner reminiscent of epithelial cells. Mechanosensing mechanisms are responsible for the range of conditions in which a cell adopts a particular migration mode. Stress fiber strengthening, specifically, is responsible for cyclical migration due to build-up of enough force to elicit rupture of focal adhesions and retraction of the cellular rear. Together, both mechanisms explain bimodal dependence of cell migration on substrate stiffness observed in the literature.

Matrix deformations around angiogenic sprouts correlate to sprout dynamics and suggest pulling activity

Angiogenesis is the formation of new blood vessels from the pre-existing vasculature. It is essential for normal tissue growth and regeneration, and also plays a key role in many diseases [Carmeliet in Nat Med 9:653-660, 2003]. Cytoskeletal components have been shown to be important for angiogenic sprout initiation and maintenance [Kniazeva and Putnam in Am J Physiol 297:C179-C187, 2009] as well as endothelial cell shape control during invasion [Elliott et al. in Nat Cell Biol 17:137-147, 2015]. The exact nature of cytoskeleton-mediated forces for sprout initiation and progression, however, remains poorly understood. Questions on the importance of tip cell pulling versus stalk cell pushing are to a large extent unanswered, which among others has to do with the difficulty of quantifying and resolving those forces in time and space. We developed methods based on time-lapse confocal microscopy and image processing-further termed 4D displacement microscopy-to acquire detailed, spatially and temporally resolved extracellular matrix (ECM) deformations, indicative of cell-ECM mechanical interactions around invading sprouts. We demonstrate that matrix deformations dependent on actin-mediated force generation are spatio-temporally correlated with sprout morphological dynamics. Furthermore, sprout tips were found to exert radially pulling forces on the extracellular matrix, which were quantified by means of a computational model of collagen ECM mechanics. Protrusions from extending sprouts mostly increase their pulling forces, while retracting protrusions mainly reduce their pulling forces. Displacement microscopy analysis further unveiled a characteristic dipole-like deformation pattern along the sprout direction that was consistent among seemingly very different sprout shapes-with oppositely oriented displacements at sprout tip versus sprout base and a transition zone of negligible displacements in between. These results demonstrate that sprout-ECM interactions are dominated by pulling forces and underline the key role of tip cell pulling for sprouting angiogenesis.

The role of actin protrusion dynamics in cell migration through a degradable viscoelastic extracellular matrix: Insights from a computational model

Actin protrusion dynamics plays an important role in the regulation of three-dimensional (3D) cell migration. Cells form protrusions that adhere to the surrounding extracellular matrix (ECM), mechanically probe the ECM and contract in order to displace the cell body. This results in cell migration that can be directed by the mechanical anisotropy of the ECM. However, the subcellular processes that regulate protrusion dynamics in 3D cell migration are difficult to investigate experimentally and therefore not well understood. Here, we present a computational model of cell migration through a degradable viscoelastic ECM. This model is a 2D representation of 3D cell migration. The cell is modeled as an active deformable object that captures the viscoelastic behavior of the actin cortex and the subcellular processes underlying 3D cell migration. The ECM is regarded as a viscoelastic material, with or without anisotropy due to fibrillar strain stiffening, and modeled by means of the meshless Lagrangian smoothed particle hydrodynamics (SPH) method. ECM degradation is captured by local fluidization of the material and permits cell migration through the ECM. We demonstrate that changes in ECM stiffness and cell strength affect cell migration and are accompanied by changes in number, lifetime and length of protrusions. Interestingly, directly changing the total protrusion number or the average lifetime or length of protrusions does not affect cell migration. A stochastic variability in protrusion lifetime proves to be enough to explain differences in cell migration velocity. Force-dependent adhesion disassembly does not result in faster migration, but can make migration more efficient. We also demonstrate that when a number of simultaneous protrusions is enforced, the optimal number of simultaneous protrusions is one or two, depending on ECM anisotropy. Together, the model provides non-trivial new insights in the role of protrusions in 3D cell migration and can be a valuable contribution to increase the understanding of 3D cell migration mechanics.

Fast quantitative time lapse displacement imaging of endothelial cell invasion

In order to unravel rapid mechano-chemical feedback mechanisms in sprouting angiogenesis, we combine selective plane illumination microscopy (SPIM) and tailored image registration algorithms - further referred to as SPIM-based displacement microscopy - with an in vitro model of angiogenesis. SPIM successfully tackles the problem of imaging large volumes while upholding the spatial resolution required for the analysis of matrix displacements at a subcellular level. Applied to in vitro angiogenic sprouts, this unique methodological combination relates subcellular activity - minute to second time scale growing and retracting of protrusions - of a multicellular systems to the surrounding matrix deformations with an exceptional temporal resolution of 1 minute for a stack with multiple sprouts simultaneously or every 4 seconds for a single sprout, which is 20 times faster than with a conventional confocal setup. Our study reveals collective but non-synchronised, non-continuous activity of adjacent sprouting cells along with correlations between matrix deformations and protrusion dynamics.

Effect of ultrasound on bone fracture healing: A computational mechanobioregulatory model

Bone healing process is a complicated phenomenon regulated by biochemical and mechanical signals. Experimental studies have shown that ultrasound (US) accelerates bone ossification and has a multiple influence on cell differentiation and angiogenesis. In a recent work of the authors, a bioregulatory model for providing bone-healing predictions was addressed, taking into account for the first time the salutary effect of US on the involved angiogenesis. In the present work, a mechanobioregulatory model of bone solidification under the US presence incorporating also the mechanical environment on the regeneration process, which is known to affect cellular processes, is presented. An iterative procedure is adopted, where the finite element method is employed to compute the mechanical stimuli at the linear elastic phases of the poroelastic callus region and a coupled system of partial differential equations to simulate the enhancement by the US cell angiogenesis process and thus the oxygen concentration in the fractured area. Numerical simulations with and without the presence of US that illustrate the influence of progenitor cells' origin in the healing pattern and the healing rate and simultaneously demonstrate the salutary effect of US on bone repair are presented and discussed.

Combustion-derived particles inhibit in vitro human lung fibroblast-mediated matrix remodeling

BACKGROUND

The continuously growing human exposure to combustion-derived particles (CDPs) drives in depth investigation of the involved complex toxicological mechanisms of those particles. The current study evaluated the hypothesis that CDPs could affect cell-induced remodeling of the extracellular matrix due to their underlying toxicological mechanisms. The effects of two ultrafine and one fine form of CDPs on human lung fibroblasts (MRC-5 cell line) were investigated, both in 2D cell culture and in 3D collagen type I hydrogels. A multi-parametric analysis was employed.

RESULTS

In vitro dynamic 3D analysis of collagen matrices showed that matrix displacement fields induced by human lung fibroblasts are disturbed when exposed to carbonaceous particles, resulting in inhibition of matrix remodeling. In depth analysis using general toxicological assays revealed that a plausible explanation comprises a cascade of numerous detrimental effects evoked by the carbon particles, including oxidative stress, mitochondrial damage and energy storage depletion. Also, ultrafine particles revealed stronger toxicological and inhibitory effects compared to their larger counterparts. The inhibitory effects can be almost fully restored when treating the impaired cells with antioxidants like vitamin C.

CONCLUSIONS

The unraveled in vitro pathway, by which ultrafine particles alter the fibroblasts' vital role of matrix remodeling, extends our knowledge about the contribution of these biologically active particles in impaired lung tissue repair mechanisms, and development and exacerbation of chronic lung diseases. The new insights may even pave the way to precautionary actions. The results provide justification for toxicological assessments to include mechanism-linked assays besides the traditional in vitro toxicological screening assays.

Full L1-regularized Traction Force Microscopy over whole cells

Background

Traction Force Microscopy (TFM) is a widespread technique to estimate the tractions that cells exert on the surrounding substrate. To recover the tractions, it is necessary to solve an inverse problem, which is ill-posed and needs regularization to make the solution stable. The typical regularization scheme is given by the minimization of a cost functional, which is divided in two terms: the error present in the data or data fidelity term; and the regularization or penalty term. The classical approach is to use zero-order Tikhonov or L₂-regularization, which uses the L₂-norm for both terms in the cost function. Recently, some studies have demonstrated an improved performance using L₁-regularization (L₁-norm in the penalty term) related to an increase in the spatial resolution and sensitivity of the recovered traction field. In this manuscript, we present a comparison between the previous two regularization schemes (relying in the L₂-norm for the data fidelity term) and the full L₁-regularization (using the L₁-norm for both terms in the cost function) for synthetic and real data.

Results

Our results reveal that L₁-regularizations give an improved spatial resolution (more important for full L₁-regularization) and a reduction in the background noise with respect to the classical zero-order Tikhonov regularization. In addition, we present an approximation, which makes feasible the recovery of cellular tractions over whole cells on typical full-size microscope images when working in the spatial domain.

Conclusions

The proposed full L₁-regularization improves the sensitivity to recover small stress footprints. Moreover, the proposed method has been validated to work on full-field microscopy images of real cells, what certainly demonstrates it is a promising tool for biological applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1771-0) contains supplementary material, which is available to authorized users.

3D full-field quantification of cell-induced large deformations in fibrillar biomaterials by combining non-rigid image registration with label-free second harmonic generation

To advance our current understanding of cell-matrix mechanics and its importance for biomaterials development, advanced three-dimensional (3D) measurement techniques are necessary. Cell-induced deformations of the surrounding matrix are commonly derived from the displacement of embedded fiducial markers, as part of traction force microscopy (TFM) procedures. However, these fluorescent markers may alter the mechanical properties of the matrix or can be taken up by the embedded cells, and therefore influence cellular behavior and fate. In addition, the currently developed methods for calculating cell-induced deformations are generally limited to relatively small deformations, with displacement magnitudes and strains typically of the order of a few microns and less than 10% respectively. Yet, large, complex deformation fields can be expected from cells exerting tractions in fibrillar biomaterials, like collagen. To circumvent these hurdles, we present a technique for the 3D full-field quantification of large cell-generated deformations in collagen, without the need of fiducial markers. We applied non-rigid, Free Form Deformation (FFD)-based image registration to compute full-field displacements induced by MRC-5 human lung fibroblasts in a collagen type I hydrogel by solely relying on second harmonic generation (SHG) from the collagen fibrils. By executing comparative experiments, we show that comparable displacement fields can be derived from both fibrils and fluorescent beads. SHG-based fibril imaging can circumvent all described disadvantages of using fiducial markers. This approach allows measuring 3D full-field deformations under large displacement (of the order of 10 μm) and strain regimes (up to 40%). As such, it holds great promise for the study of large cell-induced deformations as an inherent component of cell-biomaterial interactions and cell-mediated biomaterial remodeling.

Fibrin structural and diffusional analysis suggests that fibers are permeable to solute transport

Fibrin hydrogels are promising carrier materials in tissue engineering. They are biocompatible and easy to prepare, they can bind growth factors and they can be prepared from a patient's own blood. While fibrin structure and mechanics have been extensively studied, not much is known about the relation between structure and diffusivity of solutes within the network. This is particularly relevant for solutes with a size similar to that of growth factors. A novel methodological approach has been used in this study to retrieve quantitative structural characteristics of fibrin hydrogels, by combining two complementary techniques, namely confocal fluorescence microscopy with a fiber extraction algorithm and turbidity measurements. Bulk rheological measurements were conducted to determine the impact of fibrin hydrogel structure on mechanical properties. From these measurements it can be concluded that variations in the fibrin hydrogel structure have a large impact on the rheological response of the hydrogels (up to two orders of magnitude difference in storage modulus) but only a moderate influence on the diffusivity of dextran solutes (up to 25% difference). By analyzing the diffusivity measurements by means of the Ogston diffusion model we further provide evidence that individual fibrin fibers can be semi-permeable to solute transport, depending on the average distance between individual protofibrils. This can be important for reducing mass transport limitations, for modulating fibrinolysis and for growth factor binding, which are all relevant for tissue engineering.

STATEMENT OF SIGNIFICANCE

Fibrin is a natural biopolymer that has drawn much interest as a biomimetic carrier in tissue engineering applications. We hereby use a novel combined approach for the structural characterization of fibrin networks based on optical microscopy and light scattering methods that can also be applied to other fibrillar hydrogels, like collagen. Furthermore, our findings on the relation between solute transport and fibrin structural properties can lead to the optimized design of fibrin hydrogel constructs for controlled release applications. Finally, we provide new evidence for the fact that fibrin fibers may be permeable for solutes with a molecular weight comparable to that of growth factors. This finding may open new avenues for tailoring mass transport properties of fibrin carriers.

Computational modeling of bone fracture non-unions: four clinically relevant case studies

The human skeleton has a remarkable regeneration capacity. Nevertheless, 5–10 % of the bone fractures fails to heal and develops into a non-union which is a challenging orthopedic complication requiring complex and expensive treatment. This review paper will discuss four different computational models, each capturing a particular clinical case of non-union: non-union induced by reaming of the marrow canal and periosteal stripping, non-union due to a large interfragmentary gap, non-union due to a genetic disorder [i.e. NF1 related congenital pseudoarthrosis of the tibia (CPT)] and non-union due to mechanical overload. Together, the four computational models are able to capture the etiology of a wide range of fracture non-union types and design novel treatment strategies thereof. Further research is required to corroborate the computational models in both animal and human settings and translate them from bench to bed side.

Free Form Deformation-Based Image Registration Improves Accuracy of Traction Force Microscopy

Traction Force Microscopy (TFM) is a widespread method used to recover cellular tractions from the deformation that they cause in their surrounding substrate. Particle Image Velocimetry (PIV) is commonly used to quantify the substrate’s deformations, due to its simplicity and efficiency. However, PIV relies on a block-matching scheme that easily underestimates the deformations. This is especially relevant in the case of large, locally non-uniform deformations as those usually found in the vicinity of a cell’s adhesions to the substrate. To overcome these limitations, we formulate the calculation of the deformation of the substrate in TFM as a non-rigid image registration process that warps the image of the unstressed material to match the image of the stressed one. In particular, we propose to use a B-spline -based Free Form Deformation (FFD) algorithm that uses a connected deformable mesh to model a wide range of flexible deformations caused by cellular tractions. Our FFD approach is validated in 3D fields using synthetic (simulated) data as well as with experimental data obtained using isolated endothelial cells lying on a deformable, polyacrylamide substrate. Our results show that FFD outperforms PIV providing a deformation field that allows a better recovery of the magnitude and orientation of tractions. Together, these results demonstrate the added value of the FFD algorithm for improving the accuracy of traction recovery.

A mechano-regulatory model for bone healing predictions under the influence of ultrasound

The bone healing process involves a sequence of cellular action and interaction, regulated by biochemical and mechanical signals. Experimental studies have shown that ultrasound accelerates bone solidification and enhances the underlying healing mechanisms. An integrated computational model is presented for deriving predictions of bone healing under the presence of ultrasound.

Size does matter: an integrative in vivo-in silico approach for the treatment of critical size bone defects

Although bone has a unique restorative capacity, i.e., it has the potential to heal scarlessly, the conditions for spontaneous bone healing are not always present, leading to a delayed union or a non-union. In this work, we use an integrative in vivo - in silico approach to investigate the occurrence of non-unions, as well as to design possible treatment strategies thereof. The gap size of the domain geometry of a previously published mathematical model was enlarged in order to study the complex interplay of blood vessel formation, oxygen supply, growth factors and cell proliferation on the final healing outcome in large bone defects. The multiscale oxygen model was not only able to capture the essential aspects of in vivo non-unions, it also assisted in understanding the underlying mechanisms of action, i.e., the delayed vascularization of the central callus region resulted in harsh hypoxic conditions, cell death and finally disrupted bone healing. Inspired by the importance of a timely vascularization, as well as by the limited biological potential of the fracture hematoma, the influence of the host environment on the bone healing process in critical size defects was explored further. Moreover, dependent on the host environment, several treatment strategies were designed and tested for effectiveness. A qualitative correspondence between the predicted outcomes of certain treatment strategies and experimental observations was obtained, clearly illustrating the model's potential. In conclusion, the results of this study demonstrate that due to the complex non-linear dynamics of blood vessel formation, oxygen supply, growth factor production and cell proliferation and the interactions thereof with the host environment, an integrative in silico-in vivo approach is a crucial tool to further unravel the occurrence and treatments of challenging critical sized bone defects.

In 5–10% of fracture patients, the bone fractures do not heal in the normal healing period (delayed healing) or do not heal at all (non-union). In order to investigate the causes of impaired healing and design potential treatment strategies, we have used a combined experimental and computational approach. More specifically, large bone defects were analyzed in mouse models and simulated by a previously published computational model. After showing that the predictions of the computational model match the observations of the experimental model, we have used the computational model to investigate the underlying mechanisms of action. In particular, the results indicated that the new blood vessels do not reach the central fracture zone in time due to the large defect size, which leads to insufficient oxygen delivery, increased cell death and disrupted bone healing. The healing, however, could be rescued by adequate blood vessel ingrowth from the overlying soft tissues. Moreover, potential treatment strategies were designed based on the influence of these soft tissues. In conclusion, this study demonstrates the potential of a combined experimental and computational approach to contribute to the understanding of pathological processes like the impaired bone regeneration in large bone defects and design future treatments thereof.

Reporter cell activity within hydrogel constructs quantified from oxygen-independent bioluminescence

By providing a three-dimensional (3D) support to cells, hydrogels offer a more relevant in vivo tissue-like environment as compared to two-dimensional cell cultures. Hydrogels can be applied as screening platforms to investigate in 3D the role of biochemical and biophysical cues on cell behaviour using bioluminescent reporter cells. Gradients in oxygen concentration that result from the interplay between molecular transport and cell metabolism can however cause substantial variability in the observed bioluminescent reporter cell activity. To assess the influence of these oxygen gradients on the emitted bioluminescence for various hydrogel geometries, a combined experimental and modelling approach was implemented. We show that the applied model is able to predict oxygen gradient independent bioluminescent intensities which correlate better to the experimentally determined viable cell numbers, as compared to the experimentally measured bioluminescent intensities. By analysis of the bioluminescence reaction dynamics we obtained a quantitative description of cellular oxygen metabolism within the hydrogel, which was validated by direct measurements of oxygen concentration within the hydrogel. Bioluminescence peak intensities can therefore be used as a quantitative measurement of reporter cell activity within a hydrogel, but an unambiguous interpretation of these intensities requires a compensation for the influence of cell-induced oxygen gradients on the luciferase activity.

A causal relation between bioluminescence and oxygen to quantify the cell niche

Bioluminescence imaging assays have become a widely integrated technique to quantify effectiveness of cell-based therapies by monitoring fate and survival of transplanted cells. To date these assays are still largely qualitative and often erroneous due to the complexity and dynamics of local micro-environments (niches) in which the cells reside. Here, we report, using a combined experimental and computational approach, on oxygen that besides being a critical niche component responsible for cellular energy metabolism and cell-fate commitment, also serves a primary role in regulating bioluminescent light kinetics. We demonstrate the potential of an oxygen dependent Michaelis-Menten relation in quantifying intrinsic bioluminescence intensities by resolving cell-associated oxygen gradients from bioluminescent light that is emitted from three-dimensional (3D) cell-seeded hydrogels. Furthermore, the experimental and computational data indicate a strong causal relation of oxygen concentration with emitted bioluminescence intensities. Altogether our approach demonstrates the importance of oxygen to evolve towards quantitative bioluminescence and holds great potential for future microscale measurement of oxygen tension in an easily accessible manner.

Deciphering mechanical regulation of chondrogenesis in fibrin-polyurethane composite scaffolds enriched with human mesenchymal stem cells: a dual computational and experimental approach

Fibrin-polyurethane composite scaffolds support chondrogenesis of human mesenchymal stem cells (hMSCs) derived from bone marrow and due to their robust mechanical properties allow mechanical loading in dynamic bioreactors, which has been shown to increase the chondrogenic differentiation of MSCs through the transforming growth factor beta pathway. The aim of this study was to use the finite element method, mechanical testing, and dynamic in vitro cell culture experiments on hMSC-enriched fibrin-polyurethane composite scaffolds to quantitatively decipher the mechanoregulation of chondrogenesis within these constructs. The study identified compressive principal strains as the key regulator of chondrogenesis in the constructs. Although dynamic uniaxial compression did not induce chondrogenesis, multiaxial loading by combined application of dynamic compression and interfacial shear induced significant chondrogenesis at locations where all the three principal strains were compressive and had a minimum magnitude of 10%. In contrast, no direct correlation was identified between the level of pore fluid velocity and chondrogenesis. Due to the high permeability of the constructs, the pore fluid pressures could not be increased sufficiently by mechanical loading, and instead, chondrogenesis was induced by triaxial compressive deformations of the matrix with a minimum magnitude of 10%. Thus, it can be concluded that dynamic triaxial compressive deformations of the matrix is sufficient to induce chondrogenesis in a threshold-dependent manner, even where the pore fluid pressure is negligible.

Analysis of initial cell spreading using mechanistic contact formulations for a deformable cell model

Adhesion governs to a large extent the mechanical interaction between a cell and its microenvironment. As initial cell spreading is purely adhesion driven, understanding this phenomenon leads to profound insight in both cell adhesion and cell-substrate interaction. It has been found that across a wide variety of cell types, initial spreading behavior universally follows the same power laws. The simplest cell type providing this scaling of the radius of the spreading area with time are modified red blood cells (RBCs), whose elastic responses are well characterized. Using a mechanistic description of the contact interaction between a cell and its substrate in combination with a deformable RBC model, we are now able to investigate in detail the mechanisms behind this universal power law. The presented model suggests that the initial slope of the spreading curve with time results from a purely geometrical effect facilitated mainly by dissipation upon contact. Later on, the spreading rate decreases due to increasing tension and dissipation in the cell's cortex as the cell spreads more and more. To reproduce this observed initial spreading, no irreversible deformations are required. Since the model created in this effort is extensible to more complex cell types and can cope with arbitrarily shaped, smooth mechanical microenvironments of the cells, it can be useful for a wide range of investigations where forces at the cell boundary play a decisive role.

How cells spread on a newly encountered surface is an important issue, since it hints at how cells interact physically with the specific material in general. It has been shown before that many cell types have very similar early spreading behavior. This observation has been linked to the mechanical nature of the phenomenon, during which a cell cannot yet react by changing its structure and behavior. Understanding in detail how this passive spreading occurs, and what clues a cell may later respond to is the goal of this work. At the same time, the model we develop here should be very valuable for more complex situations of interacting cells, since it is able to reproduce the purely mechanical response in detail. We find that spreading is limited mainly by energy dissipation upon contact and later dissipation in the cell's cortex and that no irreversible deformation occurs during the spreading of red blood cells on an adhesive surface.

Fluorescent oxygen sensitive microbead incorporation for measuring oxygen tension in cell aggregates

Molecular oxygen is a main regulator of various cell functions. Imaging methods designed as screening tools for fast, in situ, 3D and non-interfering measurement of oxygen tension in the cellular microenvironment would serve great purpose in identifying and monitoring this vital and pivotal signalling molecule. We describe the use of dual luminophore oxygen sensitive microbeads to measure absolute oxygen concentrations in cellular aggregates. Stable microbead integration, a prerequisite for their practical application, was ensured by a site-specific delivery method that is based on the interactions between streptavidin and biotin. The spatial stability introduced by this method allowed for long term measurements of oxygen tension without interfering with the cell aggregation process. By making multiple calibration experiments we further demonstrated the potential of these sensors to measure local oxygen tension in optically dense cellular environments.

MOSAIC: a multiscale model of osteogenesis and sprouting angiogenesis with lateral inhibition of endothelial cells

The healing of a fracture depends largely on the development of a new blood vessel network (angiogenesis) in the callus. During angiogenesis tip cells lead the developing sprout in response to extracellular signals, amongst which vascular endothelial growth factor (VEGF) is critical. In order to ensure a correct development of the vasculature, the balance between stalk and tip cell phenotypes must be tightly controlled, which is primarily achieved by the Dll4-Notch1 signaling pathway. This study presents a novel multiscale model of osteogenesis and sprouting angiogenesis, incorporating lateral inhibition of endothelial cells (further denoted MOSAIC model) through Dll4-Notch1 signaling, and applies it to fracture healing. The MOSAIC model correctly predicted the bone regeneration process and recapitulated many experimentally observed aspects of tip cell selection: the salt and pepper pattern seen for cell fates, an increased tip cell density due to the loss of Dll4 and an excessive number of tip cells in high VEGF environments. When VEGF concentration was even further increased, the MOSAIC model predicted the absence of a vascular network and fracture healing, thereby leading to a non-union, which is a direct consequence of the mutual inhibition of neighboring cells through Dll4-Notch1 signaling. This result was not retrieved for a more phenomenological model that only considers extracellular signals for tip cell migration, which illustrates the importance of implementing the actual signaling pathway rather than phenomenological rules. Finally, the MOSAIC model demonstrated the importance of a proper criterion for tip cell selection and the need for experimental data to further explore this. In conclusion, this study demonstrates that the MOSAIC model creates enhanced capabilities for investigating the influence of molecular mechanisms on angiogenesis and its relation to bone formation in a more mechanistic way and across different time and spatial scales.

Modelling the effect of repositioning on the evolution of skeletal muscle damage in deep tissue injury

Deep tissue injury (DTI) is a localized area of tissue necrosis that originates in the subcutaneous layers under an intact skin and tends to develop when soft tissue is compressed for a prolonged period of time. In clinical practice, DTI is particularly common in bedridden patients and remains a serious issue in todays health care. Repositioning is generally considered to be an effective preventive measure of pressure ulcers. However, limited experimental research and no computational studies have been undertaken on this method. In this study, a methodology was developed to evaluate the influence of different repositioning intervals on the location, size and severity of DTI in bedridden patients. The spatiotemporal evolution of compressive stresses and skeletal muscle viability during the first 48 h of DTI onset was simulated for repositioning schemes in which a patient is turned every 2, 3, 4 or 6 h. The model was able to reproduce important experimental findings, including the morphology and location of DTI in human patients as well as the discrepancy between the internal tissue loads and the contact pressure at the interface with the environment. In addition, the model indicated that the severity and size of DTI were reduced by shortening the repositioning intervals. In conclusion, the computational framework presented in this study provides a promising modelling approach that can help to objectively select the appropriate repositioning scheme that is effective and efficient in the prevention of DTI.

Relating the chondrocyte gene network to growth plate morphology: from genes to phenotype

During endochondral ossification, chondrocyte growth and differentiation is controlled by many local signalling pathways. Due to crosstalks and feedback mechanisms, these interwoven pathways display a network like structure. In this study, a large-scale literature based logical model of the growth plate network was developed. The network is able to capture the different states (resting, proliferating and hypertrophic) that chondrocytes go through as they progress within the growth plate. In a first corroboration step, the effect of mutations in various signalling pathways of the growth plate network was investigated.

Occurrence and treatment of bone atrophic non-unions investigated by an integrative approach

Recently developed atrophic non-union models are a good representation of the clinical situation in which many non-unions develop. Based on previous experimental studies with these atrophic non-union models, it was hypothesized that in order to obtain successful fracture healing, blood vessels, growth factors, and (proliferative) precursor cells all need to be present in the callus at the same time. This study uses a combined in vivo-in silico approach to investigate these different aspects (vasculature, growth factors, cell proliferation). The mathematical model, initially developed for the study of normal fracture healing, is able to capture essential aspects of the in vivo atrophic non-union model despite a number of deviations that are mainly due to simplifications in the in silico model. The mathematical model is subsequently used to test possible treatment strategies for atrophic non-unions (i.e. cell transplant at post-osteotomy, week 3). Preliminary in vivo experiments corroborate the numerical predictions. Finally, the mathematical model is applied to explain experimental observations and identify potentially crucial steps in the treatments and can thereby be used to optimize experimental and clinical studies in this area. This study demonstrates the potential of the combined in silico-in vivo approach and its clinical implications for the early treatment of patients with problematic fractures.

In light of the ageing population, the occurrence of bone fractures is expected to rise substantially in the near future. In 5 to 10% of these cases, the healing process does not succeed in repairing the bone, leading to the formation of delayed unions or even non-unions. In this study we used a combination of an animal model mimicking a clinical non-union situation and a mathematical model developed for normal fracture healing to investigate both the causes of non-union formation and potential therapeutic strategies that can be applied to restart the healing process. After showing that the mathematical model is able to simulate key aspects of the non-union formation, we have used it to investigate several treatment strategies. One of these strategies, the treatment of a non-union involving a transplantation of cells from the bone marrow to the fracture site, was also tested in a pilot animal experiment. Both the simulations and the experiments showed the formation of a bony union between the fractured bone ends. In addition, we used the mathematical model to explain some unexpected experimental observations. This study demonstrates the added value of using a combination of mathematical modelling and experimental research as well the potential of using cell transplantation for the treatment of non-unions.

The influence of Young's modulus of loaded implants on bone remodeling: an experimental and numerical study in the goat knee

The aim of this study was to examine the influence of the Young's modulus of the implant material on the bone remodeling in a loaded condition. A combined animal experimental and computational study was set up. The animal experimental group comprised of 16 Saanen goats, each receiving one titanium implant (Young's modulus 110 GPa) and one high-density polyethylene (HDPE) implant (Young's modulus 1 GPa) in the left femoral condyle. Both types of implants received a titanium coating of 100 nm thickness. The implants protruded in the knee joint space and were directly weight bearing. The first group of eight goats was sacrificed after 6 weeks of loading and the second group of eight goats after 6 months of loading. The 16 femoral condyles with the 32 implants were prepared for microfocus computed tomography (micro-CT) scanning and histological sectioning. Three-dimensional trabecular bone parameters were calculated on the micro-CT images for the zones neck, middle, and apex of the implant. The percent of bone contact with the implant was measured on longitudinal histological sections. An axisymmetric finite element (FE) model was created to compare peri-implant bone strains and relative motion between a titanium and a HDPE implant for the experimental loading condition, and to assess the influence of different bone-implant interface (contact) conditions. From the statistical analysis of the 3D bone parameters, the difference between the titanium and HDPE implants was not significantly different (p > 0.05) between the zones (neck, middle, and apex) for both groups of goats. The implants could be considered in their entirety. After 6 weeks of loading, the PE implant presented lower connectivity and smaller marrow spaces in the circular region of 0-500 microm. In the region 500-1500 microm more bone volume was present for the PE implant. After 6 months, the PE implants showed more bone volume and thicker trabeculae than the titanium implants for the entire length of the implant. This effect was already present in the smallest region of interest, 0-500 microm. After 6 months more fibrous encapsulation was found around titanium implants. FE results demonstrated a substantial influence of the interface conditions on peri-implant strains and relative motion. For interface conditions that were representative for the early postoperative situation (involving press-fit and friction), differences in peri-implant bone strain distributions between titanium and HDPE could be related to the experimentally observed differences in amounts of bone and fibrous encapsulation. In contrast, differences in relative motion did not seem to play a role. Both the experimental and computational results suggest that implant stiffness can affect the peri-implant tissue response, which may be related to differences in peri-implant strains.

The remodeling of cardiovascular bioprostheses under influence of stem cell homing signal pathways

Optimizing current heart valve replacement strategies by creating living prostheses is a necessity to alleviate complications with current bioprosthetic devices such as calcification and degeneration. Regenerative medicine, mostly in vitro tissue engineering, is the forerunner of this optimization search, yet here we show the functionality of an in vivo alternative making use of 2 homing axes for stem cells. In rats we studied the signaling pathways of stem cells on implanted bioprosthetic tissue (photooxidized bovine pericardium (POP)), by gene and protein expression analysis. We found that SDF-1alpha/CXCR4 and FN/VLA4 homing axes play a role. When we implanted vascular grafts impregnated with SDF-1alpha and/or FN as carotid artery interpositions, primitive cells were attracted from the circulation. Next, bioprosthetic heart valves, constructed from POP impregnated with SDF-1alpha and/or FN, were implanted in pulmonary position. As shown by CD90, CD34 and CD117 immunofluorescent staining they became completely recellularized after 5 months, had a normal function and biomechanical properties and specifically the combination of SDF-1alpha and FN had an optimal valve-cell phenotype.