Remodelling of continuously distributed collagen fibres in soft connective tissues

Extracellular matrix remodelling plays an essential role in tissue engineering of load-bearing structures. The goal of this study is to model changes in collagen fibre content and orientation in soft connective tissues due to mechanical stimuli. A theory is presented describing the mechanical condition within the tissue and accounting for the effects of collagen fibre alignment and changes in fibre content. A fibre orientation tensor is defined to represent the continuous distribution of collagen fibre directions. A constitutive model is introduced to relate the fibre configuration to the macroscopic stress within the material. The constitutive model is extended with a structural parameter, the fibre volume fraction, to account for the amount of fibres present within the material. It is hypothesised that collagen fibre reorientation is induced by macroscopic deformations and the amount of collagen fibres is assumed to increase with the mean fibre stretch. The capabilities of the model are demonstrated by considering remodelling within a biaxially stretched cube. The model is then applied to analyse remodelling within a closed stented aortic heart valve. The computed preferred fibre orientation runs from commissure to commissure and resembles the fibre directions in the native aortic valve.     

The effect of collagen reinforcement in the behaviour of the temporomandibular joint disc

In this paper, the influence of collagen fibres in the behaviour of the temporomandibular joint disc is studied. A three-dimensional finite element model of the joint is developed from a set of medical images. The model comprises the mandible, part of the cranium and both temporomandibular joints. Joints have been considered to be composed of the articular discs and the temporomandibular ligaments. A fibre-reinforced porohyperelastic model was used to study the response under clenching of the fibrocartilage that composes the articular disc. This was divided in an intermediate zone, and two bands, an anterior and other posterior, in order to define the orientation of collagen fibres. The study demonstrates that the introduction of collagen fibres in the biphasic behaviour of the articular disc implies for a prescribed displacement not only an increase of the pressurization in the tissue, but also higher stresses in the anterior and posterior bands, as well as in the lateral zone of the disc. Thus, modelling the disc as an isotropic solid matrix leads in this case to an overestimation of the stresses in the intermediate zone, an underestimation of the pore pressure in this area, and an underestimation of the stresses in the rest of the disc.     

Remodeling of the collagen fiber architecture due to compaction in small vessels under tissue engineered conditions

Mechanical loading protocols in tissue engineering (TE) aim to improve the deposition of a properly organized collagen fiber network. In addition to collagen remodeling, these conditioning protocols can result in tissue compaction. Tissue compaction is beneficial to tissue collagen alignment, yet it may lead to a loss of functionality of the TE construct due to changes in geometry after culture. Here, a mathematical model is presented to relate the changes in collagen architecture to the local compaction within a TE small blood vessel, assuming that under static conditions, compaction is the main factor responsible for collagen fiber organization. An existing structurally based model is extended to incorporate volumetric tissue compaction. Subsequently, the model is applied to describe the collagen architecture of TE constructs under either strain based or stress based stimulus functions. Our computations indicate that stress based simulations result in a helical collagen fiber distribution along the vessel wall. The helix pitch angle increases from a circumferential direction in the inner wall, over about 45 deg in the middle vessel layer, to a longitudinal direction in the outer wall. These results are consistent with experimental data from TE small diameter blood vessels. In addition, our results suggest a stress dependent remodeling of the collagen, suggesting that cell traction is responsible for collagen orientation. These findings may be of value to design improved mechanical conditioning protocols to optimize the collagen architecture in engineered tissues.     

Mechanical stimulation to stimulate formation of a physiological collagen architecture in tissue-engineered cartilage: a numerical study

The load-bearing capacity of today's tissue-engineered (TE) cartilage is insufficient. The arcade-like collagen network in native cartilage plays an important role in its load-bearing properties. Inducing the formation of such collagen architecture in engineered cartilage can, therefore, enhance mechanical properties of TE cartilage. Considering the well-defined relationship between tensile strains and collagen alignment in the literature, we assume that cues for inducing this orientation should come from mechanical loading. In this study, strain fields prescribed by loading conditions of unconfined compression, sliding indentation and a novel loading regime of compression-sliding indentation are numerically evaluated to assess the probability that these would trigger a physiological collagen architecture. Results suggest that sliding indentation is likely to stimulate the formation of an appropriate superficial zone with parallel fibres. Adding lateral compression may stimulate the formation of a deep zone with perpendicularly aligned fibres. These insights may be used to improve loading conditions for cartilage tissue engineering.     

The effect of tissue-engineered cartilage biomechanical and biochemical properties on its post-implantation mechanical behavior

The insufficient load-bearing capacity of today’s tissue-engineered (TE) cartilage limits its clinical application. Focus has been on engineering cartilage with enhanced mechanical stiffness by reproducing native biochemical compositions. More recently, depth dependency of the biochemical content and the collagen network architecture has gained interest. However, it is unknown whether the mechanical performance of TE cartilage would benefit more from higher content of biochemical compositions or from achieving an appropriate collagen organization. Furthermore, the relative synthesis rate of collagen and proteoglycans during the TE process may affect implant performance. Such insights would assist tissue engineers to focus on those aspects that are most important. The aim of the present study is therefore to elucidate the relative importance of implant ground substance stiffness, collagen content, and collagen architecture of the implant, as well as the synthesis rate of the biochemical constituents for the post-implantation mechanical behavior of the implant. We approach this by computing the post-implantation mechanical conditions using a composition-based fibril-reinforced poro-viscoelastic swelling model of the medial tibia plateau. Results show that adverse implant composition and ultrastructure may lead to post-implantation excessive mechanical loads, with collagen orientation being the most critical variable. In addition, we predict that a faster synthesis rate of proteoglycans compared to that of collagen during TE culture may result in excessive loads on collagen fibers post-implantation. This indicates that even with similar final contents, constructs may behave differently depending on their development. Considering these aspects may help to engineer TE cartilage implants with improved survival rates.     

A physically motivated constitutive model for cell-mediated compaction and collagen remodeling in soft tissues

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.     

Importance of collagen orientation and depth-dependent fixed charge densities of cartilage on mechanical behavior of chondrocytes

The collagen network and proteoglycan matrix of articular cartilage are thought to play an important role in controlling the stresses and strains in and around chondrocytes, in regulating the biosynthesis of the solid matrix, and consequently in maintaining the health of diarthrodial joints. Understanding the detailed effects of the mechanical environment of chondrocytes on cell behavior is therefore essential for the study of the development, adaptation, and degeneration of articular cartilage. Recent progress in macroscopic models has improved our understanding of depth-dependent properties of cartilage. However, none of the previous works considered the effect of realistic collagen orientation or depth-dependent negative charges in microscopic models of chondrocyte mechanics. The aim of this study was to investigate the effects of the collagen network and fixed charge densities of cartilage on the mechanical environment of the chondrocytes in a depth-dependent manner. We developed an anisotropic, inhomogeneous, microstructural fibril-reinforced finite element model of articular cartilage for application in unconfined compression. The model consisted of the extracellular matrix and chondrocytes located in the superficial, middle, and deep zones. Chondrocytes were surrounded by a pericellular matrix and were assumed spherical prior to tissue swelling and load application. Material properties of the chondrocytes, pericellular matrix, and extracellular matrix were obtained from the literature. The loading protocol included a free swelling step followed by a stress-relaxation step. Results from traditional isotropic and transversely isotropic biphasic models were used for comparison with predictions from the current model. In the superficial zone, cell shapes changed from rounded to elliptic after free swelling. The stresses and strains as well as fluid flow in cells were greatly affected by the modulus of the collagen network. The fixed charge density of the chondrocytes, pericellular matrix, and extracellular matrix primarily affected the aspect ratios (height/width) and the solid matrix stresses of cells. The mechanical responses of the cells were strongly location and time dependent. The current model highlights that the collagen orientation and the depth-dependent negative fixed charge densities of articular cartilage have a great effect in modulating the mechanical environment in the vicinity of chondrocytes, and it provides an important improvement over earlier models in describing the possible pathways from loading of articular cartilage to the mechanical and biological responses of chondrocytes.     

Structural three-dimensional constitutive law for the passive myocardium

A three-dimensional constitutive law is proposed for the myocardium. Its formulation is based on a structural approach in which the total strain energy of the tissue is the sum of the strain energies of its constituents: the muscle fibers, the collagen fibers and the fluid matrix which embeds them. The ensuing material law expresses the specific structural and mechanical properties of the tissue, namely, the spatial orientation of the comprising fibers, their waviness in the unstressed state and their stress-strain behavior when stretched. Having assumed specific functional forms for the distribution of the fibers spatial orientation and waviness, the results of biaxial mechanical tests serve for the estimation of the material constants appearing in the constitutive equations. A very good fit is obtained between the measured and the calculated stresses, indicating the suitability of the proposed model for describing the mechanical behavior of the passive myocardium. Moreover, the results provide general conclusions concerning the structural basis for the tissue overall mechanical properties, the main of which is that the collagen matrix, though comprising a relatively small fraction of the whole tissue volume, is the dominant component accounting for its stiffness.     

The Evolution of Collagen Fiber Orientation in Engineered Cardiovascular Tissues Visualized by Diffusion Tensor Imaging

The collagen architecture is the major determinant of the function and mechanical behavior of cardiovascular tissues. In order to engineer a functional and load-bearing cardiovascular tissue with a structure that mimics the native tissue to meet in vivo mechanical demands, a complete understanding of the collagen orientation mechanism is required. Several methods have been used to visualize collagen architecture in tissue-engineered (TE) constructs, but they either have a limited imaging depth or have a complicated set up. In this study, Diffusion Tensor Imaging (DTI) is explored as a fast and reliable method to visualize collagen arrangement, and Confocal Laser Scanning Microscopy (CLSM) was used as a validation technique. Uniaxially constrained TE strips were cultured for 2 days, 10 days, 3 and 6 weeks to investigate the evolution of the collagen orientation with time. Moreover, a comparison of the collagen orientation in high and low aspect ratio (length/width) TE constructs was made with both methods. Both methods showed similar fiber orientation in TE constructs. Collagen fibers in the high aspect ratio samples were mostly aligned in the constrained direction, while the collagen fibers in low aspect ratio strips were mainly oriented in the oblique direction. The orientation changed to the oblique direction by extending culture time and could also be visualized. DTI captured the collagen orientation differences between low and high aspect ratio samples and with time. Therefore, it can be used as a fast, non-destructive and reliable tool to study the evolution of the collagen orientation in TE constructs.     

Contribution of postnatal collagen reorientation to depth-dependent mechanical properties of articular cartilage

The collagen fibril network is an important factor for the depth-dependent mechanical behaviour of adult articular cartilage (AC). Recent studies show that collagen orientation is parallel to the articular surface throughout the tissue depth in perinatal animals, and that the collagen orientations transform to a depth-dependent arcade-like structure in adult animals. Current understanding on the mechanobiology of postnatal AC development is incomplete. In the current paper, we investigate the contribution of collagen fibril orientation changes to the depth-dependent mechanical properties of AC. We use a composition-based finite element model to simulate in a 1-D confined compression geometry the effects of ten different collagen orientation patterns that were measured in developing sheep. In initial postnatal life, AC is mostly subject to growth and we observe only small changes in depth-dependent mechanical behaviour. Functional adaptation of depth-dependent mechanical behaviour of AC takes place in the second half of life before puberty. Changes in fibril orientation alone increase cartilage stiffness during development through the modulation of swelling strains and osmotic pressures. Changes in stiffness are most pronounced for small stresses and for cartilage adjacent to the bone. We hypothesize that postnatal changes in collagen fibril orientation induce mechanical effects that in turn promote these changes. We further hypothesize that a part of the depth-dependent postnatal increase in collagen content in literature is initiated by the depth-dependent postnatal increase in fibril strain due to collagen fibril reorientation.     

Mechanical behavior of collagen-fibrin co-gels reflects transition from series to parallel interactions with increasing collagen content

Fibrin and collagen, biopolymers occurring naturally in the body, are biomaterials commonly-used as scaffolds for tissue engineering. How collagen and fibrin interact to confer macroscopic mechanical properties in collagen-fibrin composite systems remains poorly understood. In this study, we formulated collagen-fibrin co-gels at different collagen-to-fibrin ratios to observe changes in the overall mechanical behavior and microstructure. A modeling framework of a two-network system was developed by modifying our micro-scale model, considering two forms of interaction between the networks: (a) two interpenetrating but noninteracting networks ("parallel"), and (b) a single network consisting of randomly alternating collagen and fibrin fibrils ("series"). Mechanical testing of our gels show that collagen-fibrin co-gels exhibit intermediate properties (UTS, strain at failure, tangent modulus) compared to those of pure collagen and fibrin. The comparison with model predictions show that the parallel and series model cases provide upper and lower bounds, respectively, for the experimental data, suggesting that a combination of such interactions exists between the collagen and fibrin in co-gels. A transition from the series model to the parallel model occurs with increasing collagen content, with the series model best describing predominantly fibrin co-gels, and the parallel model best describing predominantly collagen co-gels.     

Influence of substrate stiffness on circulating progenitor cell fate

In situ vascular tissue engineering (TE) aims at regenerating vessels using implanted synthetic scaffolds. An envisioned strategy is to capture and differentiate progenitor cells from the bloodstream into the porous scaffold to initiate tissue formation. Among these cells are the endothelial colonies forming cells (ECFCs) that can differentiate into endothelial cells and transdifferentiate into smooth muscle cells under biochemical stimulation. The influence of mechanical stimulation is unknown, but relevant for in situ vascular TE because the cells perceive a change in mechanical environment when captured inside the scaffold, where they are shielded from blood flow induced shear stresses. Here we investigate the effects of substrate stiffness as one of the environmental mechanical cues to control ECFC fate within scaffolds. ECFCs were seeded on soft (3.58±0.90 kPa), intermediate (21.59±2.91 kPa), and stiff (93.75±18.36 kPa) fibronectin-coated polyacrylamide gels, as well as on glass controls, and compared to peripheral blood mononuclear cells (PBMC). Cell behavior was analyzed in terms of adhesion (vinculin staining), proliferation (BrdU), phenotype (CD31, αSMA staining, and flow cytometry), and collagen production (col I, III, and IV). While ECFCs adhesion and proliferation increased with substrate stiffness, no change in phenotype was observed. The cells produced no collagen type I, but abundant amounts of collagen type III and IV, albeit in a stiffness-dependent organization. PBMCs did not adhere to the gels, but they did adhere to glass, where they expressed CD31 and collagen type III. Addition mechanical cues, such as cyclic strains, should be studied to further investigate the effect of the mechanical environment on captured circulating cells for in situ TE purposes.     

Introducing mesoscopic information into constitutive equations for arterial walls

We propose a new elastic constitutive law for arterial tissue in which the limiting polymeric chain extensibility of both collagen and elastin fibres is accounted for. The elastic strain-energy function is separated additively into two parts: an isotropic contribution associated with the matrix (incorporating the elastin fibre network) and an anisotropic one associated with the collagen fibres. Information on the limiting extensibility in each case provides some mesoscopic input into the model. The (logarithm-based) model is compared with the Fung-Demiray exponential model and certain other recently proposed models. Some aspects of the elastic response under extension and inflation of a thin-walled circular cylindrical tube (the artery) are then examined and compared with the corresponding response of a rubber-like tube. We point out that our model, when both isotropic and anisotropic terms are included, can be developed to accommodate changing mechanical properties associated with degradation of the elastin and collagen by considering the material constants that define the limit of chain extensibility to evolve in time.     

Multiscale mechanobiology: Coupling models of adhesion kinetics and nonlinear tissue mechanics

The mechanical behavior of tissues at the macroscale is tightly coupled to cellular activity at the microscale. Dermal wound healing is a prominent example of a complex system in which multiscale mechanics regulate restoration of tissue form and function. In cutaneous wound healing, a fibrin matrix is populated by fibroblasts migrating in from a surrounding tissue made mostly out of collagen. Fibroblasts both respond to mechanical cues, such as fiber alignment and stiffness, as well as exert active stresses needed for wound closure.Here, we develop a multiscale model with a two-way coupling between a microscale cell adhesion model and a macroscale tissue mechanics model. Starting from the well-known model of adhesion kinetics proposed by Bell, we extend the formulation to account for nonlinear mechanics of fibrin and collagen and show how this nonlinear response naturally captures stretch-driven mechanosensing. We then embed the new nonlinear adhesion model into a custom finite element implementation of tissue mechanical equilibrium. Strains and stresses at the tissue level are coupled with the solution of the microscale adhesion model at each integration point of the finite element mesh. In addition, solution of the adhesion model is coupled with the active contractile stress of the cell population. The multiscale model successfully captures the mechanical response of biopolymer fibers and gels, contractile stresses generated by fibroblasts, and stress-strain contours observed during wound healing. We anticipate that this framework will not only increase our understanding of how mechanical cues guide cellular behavior in cutaneous wound healing, but will also be helpful in the study of mechanobiology, growth, and remodeling in other tissues.     

Collagen fiber alignment does not explain mechanical anisotropy in fibroblast populated collagen gels

Many load-bearing soft tissues exhibit mechanical anisotropy. In order to understand the behavior of natural tissues and to create tissue engineered replacements, quantitative relationships must be developed between the tissue structures and their mechanical behavior. We used a novel collagen gel system to test the hypothesis that collagen fiber alignment is the primary mechanism for the mechanical anisotropy we have reported in structurally anisotropic gels. Loading constraints applied during culture were used to control the structural organization of the collagen fibers of fibroblast populated collagen gels. Gels constrained uniaxially during culture developed fiber alignment and a high degree of mechanical anisotropy, while gels constrained biaxially remained isotropic with randomly distributed collagen fibers. We hypothesized that the mechanical anisotropy that developed in these gels was due primarily to collagen fiber orientation. We tested this hypothesis using two mathematical models that incorporated measured collagen fiber orientations: a structural continuum model that assumes affine fiber kinematics and a network model that allows for nonaffine fiber kinematics. Collagen fiber mechanical properties were determined by fitting biaxial mechanical test data from isotropic collagen gels. The fiber properties of each isotropic gel were then used to predict the biaxial mechanical behavior of paired anisotropic gels. Both models accurately described the isotropic collagen gel behavior. However, the structural continuum model dramatically underestimated the level of mechanical anisotropy in aligned collagen gels despite incorporation of measured fiber orientations; when estimated remodeling-induced changes in collagen fiber length were included, the continuum model slightly overestimated mechanical anisotropy. The network model provided the closest match to experimental data from aligned collagen gels, but still did not fully explain the observed mechanics. Two different modeling approaches showed that the level of collagen fiber alignment in our uniaxially constrained gels cannot explain the high degree of mechanical anisotropy observed in these gels. Our modeling results suggest that remodeling-induced redistribution of collagen fiber lengths, nonaffine fiber kinematics, or some combination of these effects must also be considered in order to explain the dramatic mechanical anisotropy observed in this collagen gel model system.     

A structural constitutive model considering angular dispersion and waviness of collagen fibres of rabbit facial veins

Background  

Structural constitutive models of vascular wall integrate information on composition and structural arrangements of tissue. In blood vessels, collagen fibres are arranged in coiled and wavy bundles and the individual collagen fibres have a deviation from their mean orientation. A complete structural constitutive model for vascular wall should incorporate both waviness and orientational distribution of fibres. We have previously developed a model, for passive properties of vascular wall, which considers the waviness of collagen fibres. However, to our knowledge there is no structural model of vascular wall which integrates both these features.     

Fibre Orientation in Human Fetal Heart and Ventricular Mechanics : A Small Perturbation Analysis

The study of the topological organisation of myocardial cells is a basic requirement for understanding the mechanical design of the normal and pathological heart. Anatomical observations show that cardiac muscle tissue has a highly specialized architecture. We have made new quantitative measurements of fibre orientation through the heart wall by means of polarized light analysis on some thick sections of human fetal heart embedded in a resin and polymerized. A small perturbation method to find an equilibrium solution in a cylindrical left ventricular (LV) geometry with fibres running on toroidal shells of revolution is used to investigate the mechanical behaviour of three human fetal hearts (FH) of 14, 20 and 33 weeks of gestational age. The results of fibre strains and stresses presented for end-systolic state show significant differences when compared to results of the cylindrical geometry with regular helicoidal fibres running on cylindrical surfaces. Moreover, the toroidal shells of revolution explain shear stresses and strains in the transverse plane which also exist in the adult heart.     

Affine versus non-affine fibril kinematics in collagen networks: theoretical studies of network behavior

The microstructure of tissues and tissue equivalents (TEs) plays a critical role in determining the mechanical properties thereof. One of the key challenges in constitutive modeling of TEs is incorporating the kinematics at both the macroscopic and the microscopic scale. Models of fibrous microstructure commonly assume fibrils to move homogeneously, that is affine with the macroscopic deformation. While intuitive for situations of fibril-matrix load transfer, the relevance of the affine assumption is less clear when primary load transfer is from fibril to fibril. The microstructure of TEs is a hydrated network of collagen fibrils, making its microstructural kinematics an open question. Numerical simulation of uniaxial extensile behavior in planar TE networks was performed with fibril kinematics dictated by the network model and by the affine model. The average fibril orientation evolved similarly with strain for both models. The individual fibril kinematics, however, were markedly different. There was no correlation between fibril strain and orientation in the network model, and fibril strains were contained by extensive reorientation. As a result, the macroscopic stress given by the network model was roughly threefold lower than the affine model. Also, the network model showed a toe region, where fibril reorientation precluded the development of significant fibril strain. We conclude that network fibril kinematics are not governed by affine principles, an important consideration in the understanding of tissue and TE mechanics, especially when load bearing is primarily by an interconnected fibril network.     

A computational model for understanding the micro-mechanics of collagen fiber network in the tunica adventitia

Abdominal aortic aneurysm is a prevalent cardiovascular disease with high mortality rates. The mechanical response of the arterial wall relies on the organizational and structural behavior of its microstructural components, and thus, a detailed understanding of the microscopic mechanical response of the arterial wall layers at loads ranging up to rupture is necessary to improve diagnostic techniques and possibly treatments. Following the common notion that adventitia is the ultimate barrier at loads close to rupture, in the present study, a finite element model of adventitial collagen network was developed to study the mechanical state at the fiber level under uniaxial loading. Image stacks of the rabbit carotid adventitial tissue at rest and under uniaxial tension obtained using multi-photon microscopy were used in this study, as well as the force–displacement curves obtained from previously published experiments. Morphological parameters like fiber orientation distribution, waviness, and volume fraction were extracted for one sample from the confocal image stacks. An inverse random sampling approach combined with a random walk algorithm was employed to reconstruct the collagen network for numerical simulation. The model was then verified using experimental stress–stretch curves. The model shows the remarkable capacity of collagen fibers to uncrimp and reorient in the loading direction. These results further show that at high stretches, collagen network behaves in a highly non-affine manner, which was quantified for each sample. A comprehensive parameter study to understand the relationship between structural parameters and their influence on mechanical behavior is presented. Through this study, the model was used to conclude important structure–function relationships that control the mechanical response. Our results also show that at loads close to rupture, the probability of failure occurring at the fiber level is up to 2%. Uncertainties in usually employed rupture risk indicators and the stochastic nature of the event of rupture combined with limited knowledge on the microscopic determinants motivate the development of such an analysis. Moreover, this study will advance the study of coupling microscopic mechanisms to rupture of the artery as a whole.