Modeling the effect of collagen fibril alignment on ligament mechanical behavior

Age-dependent changes in collagen crosslinks reduce the mechanical toughness of human meniscus

The mechanical resilience of the knee meniscus is provided by a group of structural proteins in the extracellular matrix. Aging can alter the quantity and molecular structure of these proteins making the meniscus more susceptible to debilitating tears. In this study, we determined the effect of aging on the quantity of structural proteins and collagen crosslinks in human lateral meniscus, and examined whether the quantity of these molecules was predictive of tensile toughness (area under the stress-strain curve). Two age groups were tested: a young group under 40 and an older group over 65 years old. Using mass spectrometry, we quantified the abundance of proteins and collagen crosslinks in meniscal tissue that was adjacent to the dumbbell-shaped specimens used to measure uniaxial tensile toughness parallel or perpendicular to the circumferential fiber orientation. We found that the enzymatic collagen crosslink deoxypyridinoline had a significant positive correlation with toughness, and reductions in the quantity of this crosslink with aging were associated with a loss of toughness in the ground substance and fibers. The non-enzymatic collagen crosslink carboxymethyl-lysine increased in quantity with aging, and these increases corresponded to reductions in ground substance toughness. For the collagenous (Types I, II, IV, VI, VIII) and non-collagenous structural proteins (elastin, decorin, biglycan, prolargin) analyzed in this study, only the quantity of collagen VIII was predictive of toughness. This study provides valuable insights on the structure-function relationships of the human meniscus, and how aging causes structural adaptations that weaken the tissue's mechanical integrity.

Regional shear wave speeds track regional axial stress in nonuniformly loaded fibrous soft tissues

Ligaments and tendons undergo nonuniform deformation during movement. While deformations can be imaged, it remains challenging to use such information to infer regional tissue loading. Shear wave tensiometry is a promising noninvasive technique to gauge axial stress and is premised on a tensioned beam model. However, it is unknown whether tensiometry can predict regional stress in a nonuniformly loaded structure. The objectives of this study were to (1) determine whether regional shear wave speed tracks regional axial stress in nonuniformly loaded fibrous soft tissues, and (2) determine the sensitivity of regional axial stress and shear wave speed to nonuniform load distribution and fiber alignment. We created a representative set of 12,000 dynamic finite element models of a fibrous soft tissue with probabilistic variations in fiber alignment, stiffness, and aspect ratio. In each model, we applied a randomly selected nonuniform load distribution, and then excited a shear wave and tracked its regional propagation. We found that regional shear wave speed was an excellent predictor of the regional axial stress (RMSE = 0.57 MPa) and that the nature of the regional shear wave speed-stress relationship was consistent with a tensioned beam model (R2 = 0.99). Variations in nonuniform load distribution and fiber alignment did not substantially alter the wave speed-stress relationship, particularly at higher loads. Thus, these findings suggests that shear wave tensiometry could provide a quantitative estimate of regional tissue stress in ligaments and tendons.

Contribution of Elastic and Collagen Fibers to the Mechanical Behavior of Bovine Nuchal Ligament

Ligamentum nuchae is a highly elastic tissue commonly used to study the structure and mechanics of elastin. This study combines imaging, mechanical testing, and constitutive modeling to examine the structural organization of elastic and collagen fibers and their contributions to the nonlinear stress-strain behavior of the tissue. Rectangular samples of bovine ligamentum nuchae cut in both longitudinal and transverse directions were tested in uniaxial tension. Purified elastin samples were also obtained and tested. It was observed that the stress-stretch response of purified elastin tissue follows a similar curve as the intact tissue initially, but the intact tissue shows a significant stiffening behavior for stretches above 1.29 with collagen engagement. Multiphoton and histology images confirm the elastin-dominated bulk of ligamentum nuchae interspersed with small bundles of collagen fibrils and sporadic collagen-rich regions with cellular components and ground substance. A transversely isotropic constitutive model that considers the longitudinal organization of elastic and collagen fibers was developed to describe the mechanical behavior of both intact and purified elastin tissue under uniaxial tension. These findings shed light on the unique structural and mechanical roles of elastic and collagen fibers in tissue mechanics and may aid in future use of ligamentum nuchae in tissue grafting.

Computational approach for structure generation of anisotropic particles (CASGAP) with targeted distributions of particle design and orientational order

The macroscopic properties of materials are governed by their microscopic structure which depends on the materials' composition (i.e., building blocks) and processing conditions. In many classes of synthetic, bioinspired, or natural soft and/or nanomaterials, one can find structural anisotropy in the microscopic structure due to anisotropic building blocks and/or anisotropic domains formed through the processing conditions. Experimental characterization and complementary physics-based or data-driven modeling of materials' structural anisotropy are critical for understanding structure-property relationships and enabling targeted design of materials with desired macroscopic properties. In this pursuit, to interpret experimentally obtained characterization results (e.g., scattering profiles) of soft materials with structural anisotropy using data-driven computational approaches, there is a need for creating real space three-dimensional structures of the designer soft materials with realistic physical features (e.g., dispersity in building block sizes) and anisotropy (i.e., aspect ratios of the building blocks, their orientational and positional order). These real space structures can then be used to compute and complement experimentally obtained characterization results or be used as initial configurations for physics-based simulations/calculations that can then provide training data for machine learning models. To address this need, we present a new computational approach called CASGAP - Computational Approach for Structure Generation of Anisotropic Particles - for generating any desired three dimensional real-space structure of anisotropic building blocks (modeled as particles) adhering to target distributions of particle shape, size, and positional and orientational order.

Adaptation of Fibril-Reinforced Poroviscoelastic Properties in Rabbit Collateral Ligaments 8 Weeks After Anterior Cruciate Ligament Transection

Ligaments of the knee provide stability and prevent excessive motions of the joint. Rupture of the anterior cruciate ligament (ACL), a common sports injury, results in an altered loading environment for other tissues in the joint, likely leading to their mechanical adaptation. In the collateral ligaments, the patterns and mechanisms of biomechanical adaptation following ACL transection (ACLT) remain unknown. We aimed to characterize the adaptation of elastic and viscoelastic properties of the lateral and medial collateral ligaments eight weeks after ACLT. Unilateral ACLT was performed in six rabbits, and collateral ligaments were harvested from transected and contralateral knee joints after eight weeks, and from an intact control group (eight knees from four animals). The cross-sectional areas were measured with micro-computed tomography. Stepwise tensile stress-relaxation testing was conducted up to 6% final strain, and the elastic and viscoelastic properties were characterized with a fibril-reinforced poroviscoelastic material model. We found that the cross-sectional area of the collateral ligaments in the ACL transected knees increased, the nonlinear elastic collagen network modulus of the LCL decreased, and the amount of fast relaxation in the MCL decreased. Our results indicate that rupture of the ACL leads to an early adaptation of the elastic and viscoelastic properties of the collagen fibrillar network in the collateral ligaments. These adaptations may be important to consider when evaluating whole knee joint mechanics after ACL rupture, and the results aid in understanding the consequences of ACL rupture on other tissues.

Three-dimensional analysis of load-dependent changes in the orientation of dermal collagen fibers in human skin: A pilot study

The availability of quantitative structural data on the orientation of collagen fibers is of crucial importance for understanding the behavior of connective tissues. These fibers can be visualized using a variety of imaging techniques, including second harmonic generation (SHG) microscopy. However, characterization of the collagen network requires the accurate extraction of parameters from imaging data. To this end, several automated processes have been developed to identify the preferred orientation of collagen fibers. Common methods include fast Fourier transforms and curvelet transforms, but these tools are mostly used to infer a single preferred orientation. The purpose of this pilot study was to develop an easy procedure for comprehensively comparing multiple human skin samples with the goal of analyzing load-dependent changes via SHG microscopy. We created a 3D model based upon 2D image stacks that provide fiber orientation data perpendicular and parallel to the plane of the epidermis. The SHG images were analyzed by CurveAlign to obtain angle histogram plots containing information about the multiple fiber orientations in each single image. Subsequently, contour plots of the angle histogram intensities were created to provide a useful visual plotting method to clearly show the anomalies in the angle histograms in all samples. Our results provided additional details on how the collagen network carries a load. In fact, analysis of SHG images indicated that increased stretch was accompanied by an increase in the alignment of fibers in the loading direction. Moreover, these images demonstrated that more than one type of preferred orientation is present. In particular, the 3D network of fibers appears to have two preferred orientations in the planes both perpendicular and parallel to the plane of the epidermis.

A specialized protocol for mechanical testing of isolated networks of type II collagen

The mechanical responses of most soft biological tissues rely heavily on networks of collagen fibers, thus quantifying the mechanics of both individual collagen fibers and networks of these fibers advances understanding of biological tissues in health and disease. The mechanics of type I collagen are well-studied and quantified. Yet no data exist on the tensile mechanical responses of individual type II collagen fibers nor of isolated networks comprised of type II collagen. We aimed to establish methods to facilitate studies of networked and individual type II collagen fibers within the native networked structure, specifically to establish best practices for isolating and mechanically testing type II collagen networks in tension. We systematically investigated mechanical tests of networks of type II collagen undergoing uniaxial extension, and quantified ranges for each of the important variables to help ensure that the experiment itself does not affect the measured mechanical parameters. Specifically we determined both the specimen (establishing networks of isolated collagen, the footprint and thickness of the specimen) and the mechanical test (both the device and the strain rate) to establish a repeatable and practical protocol. Mechanical testing of isolated networks of type II collagen fibers leveraging this protocol will lead to better understanding of the mechanics both of these networks and of the individual fibers. Such understanding may aid in developing and testing therapeutics, understanding inter-constituent interactions (and their roles in bulk-tissue biomechanics), investigating mechanical/biochemical modifications to networked type II collagen, and proposing, calibrating, and validating constitutive models for finite element analyses.

Diaphragm muscle fibrosis involves changes in collagen organization with mechanical implications in Duchenne muscular dystrophy

In Duchenne muscular dystrophy (DMD), diaphragm muscle dysfunction results in respiratory insufficiency, a leading cause of death in patients. Increased muscle stiffness occurs with buildup of fibrotic tissue, characterized by excessive accumulation of extracellular matrix (ECM) components such as collagen, and prevents the diaphragm from achieving the excursion lengths required for respiration. However, changes in mechanical properties are not explained by collagen amount alone and we must consider the complex structure and mechanics of fibrotic tissue. The goals of our study were to 1) determine if and how collagen organization changes with the progression of DMD in diaphragm muscle tissue and 2) predict how collagen organization influences the mechanical properties of the ECM. We first visualized collagen structure with scanning electron microscopy (SEM) images and then developed an analysis framework to quantify collagen organization and generate image-based finite-element models. Image analysis revealed increased collagen fiber straightness and alignment in mdx over wild type (WT) at 3 mo (straightness: mdx = 0.976 ± 0.0108, WT = 0.887 ± 0.0309, alignment: mdx = 0.876 ± 0.0333, WT = 0.759 ± 0.0416) and 6 mo (straightness: mdx = 0.942 ± 0.0182, WT = 0.881 ± 0.0163, alignment: mdx = 0.840 ± 0.0315, WT = 0.759 ± 0.0368). Collagen fibers retained a transverse orientation relative to muscle fibers (70°-90°) in all groups. Mechanical models predicted an increase in the transverse relative to longitudinal (muscle fiber direction) stiffness, with stiffness ratio (transverse/longitudinal) increased in mdx over WT at 3 mo (mdx = 5.45 ± 2.04, WT = 1.97 ± 0.670) and 6 mo (mdx = 4.05 ± 0.985, WT = 1.96 ± 0.506). This study revealed changes in diaphragm ECM structure and mechanics during disease progression in the mdx muscular dystrophy mouse phenotype, highlighting the need to consider the role of collagen organization on diaphragm muscle function.NEW & NOTEWORTHY Scanning electron microscopy images of decellularized diaphragm muscle from WT and mdx, Duchenne muscular dystrophy model, mice revealed that collagen fibers in the epimysium are oriented transverse to muscle fibers, with age- and disease-dependent changes in collagen arrangement. Finite-element models generated from these images predicted that changes in collagen arrangement during disease progression influence the mechanical properties of the extracellular matrix. Thus, changes in collagen fiber-level structure are implicated on tissue-level properties during fibrosis.

Development and evaluation of ligament phantoms targeted for shear wave tensiometry

Developing a shear wave tensiometer capable of non-invasively measuring ligament tension holds promise for enhancing research and clinical assessments of ligament function. Such development would benefit from tunable test specimens fabricated from well-characterized and consistent materials. Although previous work found that yarn can replicate the mechanical behavior of collateral ligaments, it is not obvious whether yarn-based phantoms would be suitable for development of a shear wave tensiometer for measuring ligament tension. Accordingly, the primary objective of this study was to characterize the mechanical properties and shear wave speed - stress relationships of ligament phantoms fabricated from yarn and silicone, and compare these results to published data from biological ligaments. We measured the mechanical properties and shear wave speeds during axial loading in nine phantoms with systematically varied material properties. We performed a simple linear regression between shear wave speed squared and axial stress to determine the shear wave speed - stress relationship for each phantom. We found comparable elastic moduli, hysteresis, and shear wave speed squared - stress regression parameters between the phantoms and collateral ligaments. For example, the ranges of the coefficients of determination (R²) and slopes across the nine phantoms were 0.84-0.95, and 0.78-1.27 kPa/m²/s², respectively, which overlapped with the ranges found in a prior study in porcine collateral ligaments (0.84-0.996 and 0.34-1.18 kPa/m²/s², respectively). Additionally, the shear wave speed squared - stress regression parameters varied predictably with the density of the phantom and the shear modulus of the silicone. In summary, we found that yarn-based phantoms serve as mechanical analogs for ligaments (i.e., are ligament mimicking), and thus, should prove beneficial for investigations into ligament structure-function relationships and in the development of a shear wave tensiometer for measuring ligament tension.

Sensitivity of the shear wave speed-stress relationship to soft tissue material properties and fiber alignment

The use of shear wave propagation to noninvasively measure material properties and loading in tendons and ligaments is a growing area of interest in biomechanics. Prior models and experiments suggest that shear wave speed primarily depends on the apparent shear modulus (i.e., shear modulus accounting for contributions from all constituents) at low loads, and then increases with axial stress when axially loaded. However, differences in the magnitudes of shear wave speeds between ligaments and tendons, which have different substructures, suggest that the tissue's composition and fiber alignment may also affect shear wave propagation. Accordingly, the objectives of this study were to (1) characterize changes in the apparent shear modulus induced by variations in constitutive properties and fiber alignment, and (2) determine the sensitivity of the shear wave speed-stress relationship to variations in constitutive properties and fiber alignment. To enable systematic variations of both constitutive properties and fiber alignment, we developed a finite element model that represented an isotropic ground matrix with an embedded fiber distribution. Using this model, we performed dynamic simulations of shear wave propagation at axial strains from 0% to 10%. We characterized the shear wave speed-stress relationship using a simple linear regression between shear wave speed squared and axial stress, which is based on an analytical relationship derived from a tensioned beam model. We found that predicted shear wave speeds were both in-range with shear wave speeds in previous in vivo and ex vivo studies, and strongly correlated with the axial stress (R² = 0.99). The slope of the squared shear wave speed-axial stress relationship was highly sensitive to changes in tissue density. Both the intercept of this relationship and the apparent shear modulus were sensitive to both the shear modulus of the ground matrix and the stiffness of the fibers' toe-region when the fibers were less well-aligned to the loading direction. We also determined that the tensioned beam model overpredicted the axial tissue stress with increasing load when the model had less well-aligned fibers. This indicates that the shear wave speed increases likely in response to a load-dependent increase in the apparent shear modulus. Our findings suggest that researchers may need to consider both the material and structural properties (i.e., fiber alignment) of tendon and ligament when measuring shear wave speeds in pathological tissues or tissues with less well-aligned fibers.

A Hybrid Microstructural-Continuum Multiscale Approach for Modeling Hyperelastic Fibrous Soft Tissue

The heterogeneous, nonlinear, anisotropic material behavior of biological tissues makes precise definition of an accurate constitutive model difficult. One possible solution to this issue would be to define microstructural elements and perform fully coupled multiscale simulation. However, for complex geometries and loading scenarios, the computational costs of such simulations can be prohibitive. Ideally then, we should seek a method that contains microstructural detail, but leverages the speed of classical continuum-based finite-element (FE) modeling. In this work, we demonstrate the use of the Holzapfel-Gasser-Ogden (HGO) model [1, 2] to fit the behavior of microstructural network models. We show that Delaunay microstructural networks can be fit to the HGO strain energy function by calculating fiber network strain energy and average fiber stretch ratio. We then use the HGO constitutive model in a FE framework to improve the speed of our hybrid model, and demonstrate that this method, combined with a material property update scheme, can match a full multiscale simulation. This method gives us flexibility in defining complex FE simulations that would be impossible, or at least prohibitively time consuming, in multiscale simulation, while still accounting for microstructural heterogeneity.

Effect of age on the failure properties of human meniscus: High-speed strain mapping of tissue tears

The knee meniscus is a soft fibrous tissue with a high incidence of injury in older populations. The objective of this study was to determine the effect of age on the failure behavior of human knee meniscus when applying uniaxial tensile loads parallel or perpendicular to the primary circumferential fiber orientation. Two age groups were tested: under 40 and over 65 years old. We paired high-speed video with digital image correlation to quantify for the first time the planar strains occurring in the tear region at precise time points, including at ultimate tensile stress, when the tissue begins losing load-bearing capacity. On average, older meniscus specimens loaded parallel to the fiber axis had approximately one-third less ultimate tensile strain and absorbed 60% less energy to failure within the tear region than younger specimens (p < 0.05). Older specimens also had significantly reduced strength and material toughness when loaded perpendicular to the fibers (p < 0.05). These age-related changes indicate a loss of collagen fiber extensibility and weakening of the non-fibrous matrix with age. In addition, we found that when loaded perpendicular to the circumferential fibers, tears propagated near the planes of maximum tensile stress and strain. Whereas when loaded parallel to the circumferential fibers, tears propagated oblique to the loading axis, closer to the planes of maximum shear stress and strain. Our experimental results can assist the selection of valid failure criteria for meniscus, and provide insight into the effect of age on the failure mechanisms of soft fibrous tissue.

Three-dimensional collagen fiber mapping and tractography of human uterine tissue using OCT

Automatic quantification and visualization of 3-D collagen fiber architecture using Optical Coherence Tomography (OCT) has previously relied on polarization information and/or prior knowledge of tissue-specific fiber architecture. This study explores image processing, enhancement, segmentation, and detection algorithms to map 3-D collagen fiber architecture from OCT images alone. 3-D fiber mapping, histogram analysis, and 3-D tractography revealed fiber groupings and macro-organization previously unseen in uterine tissue samples. We applied our method on centimeter-scale mosaic OCT volumes of uterine tissue blocks from pregnant and non-pregnant specimens revealing a complex, patient-specific network of fibrous collagen and myocyte bundles.

General method for classification of fiber families in fiber-reinforced materials: application to in-vivo human skin images

Fiber structures play a major role for the function of fiber-reinforced materials such as biological tissue. An objective classification of the fiber orientations into fiber families is crucial to understand its mechanical properties. We introduce the Fiber Image Network Evaluation Algorithm (FINE algorithm) to classify and quantify the number of fiber families in scientific images. Each fiber family is characterized by an amplitude, a mean orientation, and a dispersion. A new alignment index giving the averaged fraction of aligned fibers is defined. The FINE algorithm is validated by realistic grayscale Monte-Carlo fiber images. We apply the algorithm to an in-vivo depth scan of second harmonic generation images of dermal collagen in human skin. The derived alignment index exhibits a crossover at a critical depth where two fiber families with a perpendicular orientation around the main tension line arise. This strongly suggests the presence of a transition from the papillary to the reticular dermis. Hence, the FINE algorithm provides a valuable tool for a reliable classification and a meaningful interpretation of in-vivo collagen fiber networks and general fiber reinforced materials.

Noise reduction and quantification of fiber orientations in greyscale images

Quantification of the angular orientation distribution of fibrous tissue structures in scientific images benefits from the Fourier image analysis to obtain quantitative information. Measurement uncertainties represent a major challenge and need to be considered by propagating them in order to determine an adaptive anisotropic Fourier filter. Our adaptive filter method (AF) is based on the maximum relative uncertainty δcut of the power spectrum as well as a weighted radial sum with weighting factor α. We use a Monte-Carlo simulation to obtain realistic greyscale images that include defined variations in fiber thickness, length, and angular dispersion as well as variations in noise. From this simulation the best agreement between predefined and derived angular orientation distribution is found for evaluation parameters δcut = 2.1% and α = 1.5. The resulting cumulative orientation distribution was modeled by a sigmoid function to obtain the mean angle and the fiber dispersion. A comparison to a state-of-the-art band-pass method revealed that the AF method is more suitable for the application on greyscale fiber images, since the error of the fiber dispersion significantly decreased from (33.9 ± 26.5)% to (13.2 ± 12.7)%. Both methods were found to accurately quantify the mean fiber orientation with an error of (1.9 ± 1.5)° and (2.3 ± 2.1)° in case of the AF and the band-pass method, respectively. We demonstrate that the AF method is able to accurately quantify the fiber orientation distribution in in vivo second-harmonic generation images of dermal collagen with a mean fiber orientation error of (6.0 ± 4.0)° and a dispersion error of (9.3 ± 12.1)%.

Mechanical Roles in Formation of Oriented Collagen Fibers

Collagen is a structural protein that is widely present in vertebrates, being usually distributed in tissues in the form of fibers. In living organisms, fibers are organized in different orientations in various tissues. As the structural base in connective tissue and load-bearing tissue, the orientation of collagen fibers plays an extremely important role in the mechanical properties and physiological and biochemical functions. The study on mechanics role in formation of oriented collagen fibers enables us to understand how discrete cells use limited molecular materials to create tissues with different structures, thereby promoting our understanding of the mechanism of tissue formation from scratch, from invisible to tangible. However, the current understanding of the mechanism of fiber orientation is still insufficient. In addition, existing fabrication methods of oriented fibers are varied and involve interdisciplinary study, and the achievements of each experiment are favorable to the construction and improvement of the fiber orientation theory. To this end, this review focuses on the preparation methods of oriented fibers and proposes a model explaining the formation process of oriented fibers in tendons based on the existing fiber theory. Impact statement As the structural base in connective tissue and load-bearing tissue, the orientation of collagen fibers plays an extremely important role in the mechanical properties and physiological and biochemical functions. However, the current understanding of the mechanism of fiber orientation is still insufficient, which is greatly responsible for the challenge of functional tissue repair and regeneration. Understanding the mechanism of fiber orientation can promote the successful application of fiber orientation scaffolds in tissue repair and regeneration, as well as providing an insight for the mechanism of tissue histomorphology.

Preparing an injectable hydrogel with sodium alginate and Type I collagen to create better MSCs growth microenvironment

In the past few decades, stem cell transplantation has been generally accepted as an effective method on the treatment of tissue and organ injury. However, the insufficient number of transplanted stem cells and low survival rate that caused by series of negative conditions limit the therapeutic effect. In this contribution, we developed an injectable hydrogel composed of sodium alginate (SA) and Type I collagen (ColI), as the tissue scaffold to create better growth microenvironment for the stem cells. Compared the traditional SA scaffold, the ColI/SA hydrogel inherits its biomimetic properties, and simultaneously has shorter gelation time which means less loss of the transplanted stem cells. The mesenchyma stem cell (MSC) culture experiments indicated that the ColI/SA hydrogel could prevent the MSC apoptosis and contributed to faster MSC proliferation. It is highlighted that this ColI/SA hydrogel may have potential application for tissue regeneration and organ repair as the stem cell scaffold.

High-speed collagen fiber modeling and orientation quantification for optical coherence tomography imaging

Quantifying collagen fiber architecture has clinical and scientific relevance across a variety of tissue types and adds functionality to otherwise largely qualitative imaging modalities. Optical coherence tomography (OCT) is uniquely suited for this task due to its ability to capture the collagen microstructure over larger fields of view than traditional microscopy. Existing image processing techniques for quantifying fiber architecture, while accurate and effective, are very slow for processing large datasets and tend to lack structural specificity. We describe here a computationally efficient method for quantifying and visualizing collagen fiber organization. The algorithm is demonstrated on swine atria, bovine anterior cruciate ligament, and human cervical tissue samples. Additionally, we show an improved performance for images with crimped fiber textures and low signal to noise when compared to similar methods.

Extracellular Matrix Expression and Production in Fibroblast-Collagen Gels: Towards an In Vitro Model for Ligament Wound Healing

Ligament wound healing involves the proliferation of a dense and disorganized fibrous matrix that slowly remodels into scar tissue at the injury site. This remodeling process does not fully restore the highly aligned collagen network that exists in native tissue, and consequently repaired ligament has decreased strength and durability. In order to identify treatments that stimulate collagen alignment and strengthen ligament repair, there is a need to develop in vitro models to study fibroblast activation during ligament wound healing. The objective of this study was to measure gene expression and matrix protein accumulation in fibroblast-collagen gels that were subjected to different static stress conditions (stress-free, biaxial stress, and uniaxial stress) for three time points (1, 2 or 3 weeks). By comparing our in vitro results to prior in vivo studies, we found that stress-free gels had time-dependent changes in gene expression (col3a1, TnC) corresponding to early scar formation, and biaxial stress gels had protein levels (collagen type III, decorin) corresponding to early scar formation. This is the first study to conduct a targeted evaluation of ligament healing biomarkers in fibroblast-collagen gels, and the results suggest that biomimetic in-vitro models of early scar formation should be initially cultured under biaxial stress conditions.