Deformation and fracture of echinoderm collagen networks

Collagen networks provide the main structural component of most tissues and represent an important ingredient for bio-mimetic materials for bio-medical applications. Here we study the mechanical properties of stiff collagen networks derived from three different echinoderms and show that they exhibit non-linear stiffening followed by brittle fracture. The disordered nature of the network leads to strong sample-to-sample fluctuations in elasticity and fracture strength. We perform numerical simulations of a three dimensional model for the deformation of a cross-linked elastic fibril network which is able to reproduce the macroscopic features of the experimental results and provide insights into the internal mechanics of stiff collagen networks. Our numerical model provides an avenue for the design of collagen membranes with tunable mechanical properties.


Introduction
Living systems provide a formidable source of inspiration for the design and synthesis of new classes of materials with a potentially wide range of medical and non-medical applications [1,2].The structural and mechanical properties of animals mostly relies on collagen [3], one of the main protein constituents of the extracellular matrix in connective tissues and an essential component of mammalian skin, bones and tendons.To achieve remarkable mechanical properties, long collagen molecules self-assemble into hierarchical structures by forming fibrils that can then arrange into bundles, known as fibers, and extended elastic networks.The mechanics of fiber networks has received a wide attention both experimentally [4,5,6,7,8] and computationally [9,10,11,12,13,14,15] because of its relevance for applications but also as a model system to understand soft and disordered matter in general.
Most of the experimental studies focus on shear deformation of very soft collagen gels which are typically extremely soft at small deformations and become stiffer as the deformation proceeds [16].This peculiar non-linear strain stiffening, observed widely in soft collagen gels under shear [4,5,6,7,8] but also present in stiffer skin tissues under tension [17], has important mechanical and physiological implications since it provides at the same time stability for large deformations and facilitates soft remodeling at small deformations.The precise origin of non-linear stiffening is still under scrutiny.In biopolymer networks this general behavior is usually attributed to a crossover between bending dominated and stretching dominated elasticity [18] or to the strain induced reorientation of the fibers [4,12].The application of these ideas to collagen is still debated [8] and it is not clear if the deformation is dominated by the non-linear response of individual fibers or by collective properties of the network [19].
Marine-derived biomaterials in general -and marine collagens in particular -are promising materials for applications to tissue engineering.The mutable collagenous tissues (MCTs) of echinoderms, widespread in all five extant echinoderm classes, display striking passive mechanical properties (viscosity, tensile strength, and stiffness) that can be actively transformed within seconds by active mechanisms controlled by the nervous system [20].Hence, in addition to the mechanical functions usually associated with conventional collagenous structures (i.e.energy storage, transmission and dissipation), MCTs allow for the detachment of appendages or body parts in response to disease, trauma or predator attack [20] and for the energy-sparing maintenance of posture [21].As in mammalian connective tissues, most mutable collagenous structures are formed by parallel aggregations of collagen fibrils attached to proteoglycans [22,23].Collagen fibers are usually surrounded by an elastic network further improving the structural and mechanical properties of the tissue [20].
Besides their physiological relevance, MCTs have been recently proposed as a promising source of collagen for applications to regenerative medicine [24,23].In particular, collagen barrier-membranes of mammalian origin are typically used to facilitate proper tissue regrowth in anatomically separated compartments during guided tissue regeneration [25].Marine collagen and MCTs in particular could provide an innovative and safer alternative to mammalian ones [26,24,27].Moreover, other new intriguing applications could be derived from the use of collagen derived from echinoderms: for example to investigate cellular contractility as a biomechanics-related readout for phenotypic drug assessment [28] In this paper, we characterize the elastic tensile deformation and failure of cross-linked collagen networks derived from three echinoderms (Fig. 1): sea urchin (Paracentrotus lividus), starfish (Echinaster sepositus) and sea cucumber (Holothuria tubulosa).We analyze the stiffening behavior of the networks and characterize sample-to-sample fluctuations in the elastic modulus and tensile strength.We compare our results with numerical simulations for a three dimensional model for cross-linked elastic fibrils with parameters chosen to mimic the microstructure of our samples as revealed by SEM microscopy.Our model is similar in spirit to recent three dimensional models of cross-linked fiber networks [9,29,11,14,13] but its features and parameters are adapted to the stiff collagen networks we study experimentally.In particular, we explicitly endow each cross-linker with bending, shear and torsional rigidities [13].Simulations of the model highlight the importance of cross-linking in determining the mechanical properties of these networks, setting the bases for a computational guided design for this kind of bio-inspired materials.2

Experimental animals
Adult specimens of the starfish Echinaster sepositus, the sea urchin Paracentrotus lividus and the sea cucumber Holothuria tubulosa were collected by scuba divers in Paraggi (Marine Protected Area of Portofino, Ligurian Sea, Italy) and quickly transferred to the Department of Biosciences (University of Milan) in order to be immediately dissected.The aboral body walls (ABW) of the starfish, the peristomial membrane (PM) of the sea urchins and the whole body walls (BW) of the sea cucumbers were dissected and stored at -20 • C for the subsequent collagen extraction protocol.

Echinoderm collagen extraction
Sea urchin and starfish collagen were extracted from the PM and ABW, respectively, as previously described [24] with only slight modifications for the latter (see below).Briefly, both starfish ABW and sea urchin PM were dissected in small pieces, rinsed in artificial sea water (ASW), left in a hypotonic buffer (10 mM Tris, 0,1% EDTA) for 12 h (RT) and then in a decellularizing solution (10 mM Tris, 0,1% SDS) for 12 h (RT).After several washings in PBS, samples were placed in a disaggregating solution (0.5 M NaCl, 0,1 M Tris-HCl pH 8,0, 0,1 M β-mercaptoethanol, 0,05 M EDTA-Na).Starfish samples underwent an additional washing step in citric acid (1 mM pH 3-4) between the decellularizing and the disaggregating solutions.All these steps were performed under stirring (rotating) conditions.The obtained collagen suspension was filtered and dialyzed against 0,5 M EDTA-Na solution (pH 8,0) for 3 h (RT) and against dH2O overnight (RT).Sea cucumber BW was dissected in small pieces, placed in PBS and gentamicin (40 µg/mL) and left in stirring condition (RT) for at least 5 days until an homogeneous collagen suspension was obtained, which was then subsequently filtered.The suspensions obtained from the three experimental models were stored at -20 • C until membrane preparation.

Echinoderm-derived collagen membrane production
Membranes of the three different collagen types were prepared as previously described for sea urchin collagen matrices [24].The collagen suspensions were centrifuged for 10 minutes at 50x g and then for 20 minutes at 2000x g.The pellet was re-suspended in autoclaved dH2O (final concentration 2mg/mL) and 8,7mL of this suspension were placed in 5cm x 3 cm rubber silicone molds and left dry overnight at +37 • C. The so obtained collagen membranes were then immersed in a EDC/NHS cross-linker solution (EDC 30 mM/NHS 15 mM in MES buffer 100 mM pH 5,5) for 4 h and then washed with PBS, dH2O and ethanol 70%.

Scanning Electron Microscopy
Substrates were fixed with 2% glutaraldehyde in 0.1 M cacodylate buffer for 2 h at 4 • C and post-fixed with 1% osmium tetroxide in 0.1 M sodium cacodylate buffer (2 h, room temperature).After several washings with dH 2 O, they were dehydrated in an increasing ethanol scale and treated with a series of solutions of HMDS (Hexamethyldisilazane) and ethanol in different proportions (1:3, 1:1, 3:1 and 100% HMDS), mounted on stubs, covered with pure gold (Agar SEM Auto Sputter, Stansted, UK).The samples are observed under a scanning electron microscope (SEM) (LEO-1430, Zeiss, Oberkochen, Germany).Examples of the observed network structure is reported in Fig. 1ace.Images were analyzed using the ImageJ plugin Diameter-J [30] to extract diameter distributions (see Fig. 1bdf).3

Tensile tests
The collagen membranes were cut into narrow strips by pressing them under an assembly of evenly spaced razor blades.The ends of the collagen strips were then glued between pieces of rough emery paper to facilitate clamping.The length and width of the exposed collagen surface was 10 mm and 2.25 mm, respectively.Before testing, the thickness of each sample was measured using a LGF-01100L-B transmission-type photo-electric linear encode (Mitutoyo) with 100 nm resolution.Measurements were taken from three positions and averaged.The average thickness was in the range 4.6−18.6µm.We did not find any systematic variation of the mechanical properties with the thickness of the sample.The tensile tests were performed under controlled temperature (21-22 • C) and humidity (64-69%) with a Kammrath-Weiss B.1205.A tensile tester equipped with displacement gauge and DDS3 Controller (Kammrath & Weiss GmbH, Dortmund, Germany) and a 100 N load cell.After clamping samples into place, Leibovitz L-15 liquid was added on the surface in order to wet them completely.Additional droplets of liquid were added every two minutes during the tests to keep the samples from drying.All samples were tested under uniaxial tension until failure, with a constant strain rate (8.35 • 10 −4 s −1 ).Fig. 2 reports a representative example of a broken sample visualized under SEM microscopy.We observe a very smooth fracture surface, indicative of brittle failure.

Experimental results
The individual stress-strain curves obtained experimentally for n equivalent samples are reported in Fig. 3a for sea urchin (n = 14), in Fig. 3b for starfish (n = 20) and in Fig. 3c for sea cucumber (n = 7).The curves display a wide variability between sample to sample.The shape of each curve is, however, similar in all cases and showing a characteristic stiffening behavior with an initial non-linear regime, followed by a linear regime that ends in abrupt failure, in qualitative agreement with previous experiments on softer mammalian collagen networks [31].To reveal the typical behavior of the system, we compute the average stress strain curve by defining strain bins of size ∆ = 0.01 and computing the average stress in each bin (Fig. 3).Using the same procedure, we estimate the strain-dependent Young modulus E( ) as the average slope of the stress-strain curves in each of the strain bins.The result, plotted in Fig. 3d, shows a marked linear region where E( ) increases until it reaches a plateau E * .
To better visualize the large variations observed from sample to sample, we report in Fig. 4a the individual values of the Young moduli E * and tensile strengths σ c , defined as the peak stresses of each curve.As reported for other natural and synthetic materials [2], strength and stiffness are correlated, scaling in our case as σ c 0.12E * .The sample to sample fluctuations can be quantified by measuring the cumulative distributions of the Young moduli P(E * ) and of tensile strengths, P(σ c ), reported in Fig. 4b and 4c, respectively.The functional form of P(σ c ) is well described by the classical Weibull law P(x) = 1 − exp(−(x/x 0 ) k ) (Fig. 4d) as commonly observed in the fracture of disordered media [32].Our results indicate that sea urchin and sea cucumber collagen networks display comparable mechanical properties, while starfish derived collagen membranes are slightly stronger and stiffer.

Model
We construct a three dimensional model of the deformation of cross-linked elastic fibrils.Collagen fibrils are discretized as bead-spring polymers of diameter d with stretching and bending stiffness.The rest length of each discretization segment is a and the elastic energy of a fibril of L nodes can be written as where k f = πEd 2 /4a is the stretching stiffness, K = 3πEd 4 /64a is the bending stiffness, and θ j is the angle between two fibril segments.The fibrils are assumed to have no torsional resistance.The Young modulus of an individual fibril is chosen to be E = 10 GPa, as measured in AFM tests in dry conditions [33].To simulate fibril-fibril mutual hindrance, fibril nodes also interact by a purely repulsive Lennard-Jones (LJ) potential U(r) = 4 [(σ/r) 12 − (σ/r) 6 ] for r < σ and zero otherwise.The parameter σ is chosen to represent the typical fibril diameter (e.g.σ = 75nm).
To simulate cross-linking, we randomly place links between neighboring fibril nodes.In particular, fibrils are cross-linked with probability p when two nodes are closer than a distance b.In order to provide the cross-linkers with stretching, bending and torsional stiffness, we model each link by a set of four springs (see Fig. 5a).Considering two fibrils, A and B, a link involves four nodes, two on fibril A, i − 1 and i + 1 and two on fibril B, j − 1 and j + 1, which are connected by four springs i ± 1 → j ± 1 (Fig. 5a).We notice that internal stresses are present in real samples [34], but are not easily quantified in our case.Hence, for simplicity, we assume that initially there is no stress in the system and set the rest length of each of the four link connections to their initial distance b α i, j , with α = 1, ...4.The cross-linking energy between fibril i and j is given by dr is the actual length of the link.In most of the simulations, we assume that both fibrils and cross-linking bonds are linearly elastic and do not consider fracture events or the formation of new cross-linkers.We also perform simulations where the cross-linkers display non-linear elastic behavior as recently shown in Ref. [14].In particular, the cross-linker k l stiffness in Eq. 2 is increased from k l,1 = 13.3N/m to a value k l,2 = 1330N/m when the strain in the link is larger than a threshold 1 .
Most simulations are performed with N = 1500 fibrils randomly arranged in a 7.5 × 30 × 3,75 µm 3 block with open boundary conditions (see Fig. 5b for a typical starting configuration).We also perform additional simulations with N = 750 and N = 12000 fibrils, with all the dimensions of the simulation box doubled in the latter case.Other geometrical parameters, estimated on the basis of the SEM micrographs, are reported in Table 1.We deform the sample by holding a set of nodes on one edge of the sample fixed and moving nodes at the other edge at constant velocity.The applied strain-rate for all simulations is chosen to be ˙ yy = 10s −1 .As it is common in molecular dynamics simulations, we can only simulate strain rates that are much larger than the experimental ones.The simulations are performed with a linear damping term (F = −γv), with damping coefficient γ = 2 × 10 −14 kg/s and the mass of each element is set to M = 10 −18 kg, consistent with the estimate of Ref. [3] for a smaller volume element.The model is implemented in the LAMMPS molecular dynamics simulator package [35] and the integration time step is set to dt = 6 × 10 −8 s.

Simulations
We perform numerical simulations of the network model using the parameters reported in Table 1.The evolution of the deformation process is reported in Fig. 5c-d  observe that fibrils are first reoriented in the initial stage and then the sample becomes much thinner.This is quantified in Fig. 6 displaying the evolution with strain of the distribution of the angle θ formed by the tensile axis and fibril segments.Fig. 6 shows that as the sample deforms the typical value of θ decreases, as shown by the left-ward shift of the peak of the distributions.We checked that the effect is more marked when the sample size is larger.
From the simulations, we notice that the amount of cross-linking plays a major role in the mechanical response.This is illustrated in Fig. 7a where we show that the stress-strain curves change substantially by varying the number of cross-linkers N link .The same curves display features that are in qualitative agreement with experiments (Fig. 3).For each curve we compute the local slope that we fit with an exponential function of the strain E( ) = E 0 + (E * − E 0 )(1 − exp(− / 0 )), which allows to estimate the steady state Young modulus E * (Fig. 3b).Fig. 7cd shows that the resulting Young modulus E * increases with the number of cross-linkers N link and also depends on their stiffness k link .The results reported in Fig. 7c correspond to a set of 750 fibrils, while those in Fig. 7d correspond to 1500 fibrils.Finally, we investigate how the elastic energy is distributed in the system as it is loaded.In Fig. 8ab, we report the distribution of elastic energies E el of fibrils (Eq. 1) and cross-linkers (Eq.2).Both distributions decay exponentially with a decay constant that simply scales with the loading strain energy (e.g.E el ∝ 2 ) as illustrated in Fig. 8c.Finally Fig. 8d shows the fraction of elastic energy carried by fibrils.As the sample is loaded a part of the load carried by the cross-linkers is transferred to the fibrils, in agreement with the observation of a reorientation of the latter.We notice that the non-linearity of the simulated stress-strain curves and of the resulting Young moduli is less pronounced than in experiments (compare Fig. 3b with Fig. 7b).This suggests that geometrical effects due to fibril arrangement and reorientation are not enough to fully account for the experimentally observed stiffening.To overcome this problem, we explicitly include nonlinear effects in the cross-linker elasticity following Ref.[14].Fig. 9a shows the stress-strain curves obtained for two values of the linear cross-linker spring stiffnesses (i.e.k l = 13.3N/m and k l = 1330N/m).The curves are then compared with simulations with non-linear cross-linker springs which interpolate between the two linear spring stiffness values.The resulting stressstrain curves displays slopes that interpolates between the curves obtained with linear crosslinker springs.The local Young moduli, shown in Fig. 9b, now display a more pronounced non-linearity that is closer to experimental results.

Conclusions
In this paper, we investigate experimentally and computationally the mechanical properties of collagen networks.We focus on echinoderm-derived collagen networks, an alternative to the most widely used collagen of mammalian origin, and show that their mechanical properties are very similar to the latter.We notice that for the three types of collagen networks the resulting mechanical properties are similar, with Young moduli in the 100-500 MPa range and tensile strength in the 20-50MPa range, values that are remarkably close to what is measured in the human skin [36].This feature is possibly relevant in order to use marine collagen in biomedical applications [24,23].Slight quantitative differences exist for starfish collagen that appears to be stronger and stiffer than the other two.This could be attributed to different cross-linking properties due to differences in the sample preparation and extraction.Furthermore, as in other fiber networks studied in the past [4,5,6,7,8], our samples also display non-linear stiffening and sample-to-sample strength fluctuations typical of systems with a heterogeneous microstructure [32].
To clarify the microscopic origin of the experimentally observed behavior, we develop a computational model for the deformation of cross-linked elastic fibril networks.The model allows to reproduce some key features of the experimental results such as the non-linear stiffening behavior.If we use linear-elastic cross-linkers, we find non-linear stress-strain curves induced by a reorientation of the fibrils.The resulting non-linear strengthening is, however, weaker than the one measured experimentally suggesting that other mechanisms playing a role.When we include non-linear elastic cross-linkers in the model, as suggested by Ref. [14], we obtain a more pronounced stiffening.Additional sources of non-linearity may come from non-linear elastic deformation of the fibril themselves [3] and from internal stresses due to sample preparation [34].Furthermore, our model highlights the crucial role played by the amount of cross-linking in determining the mechanical response of the sample.This suggests an interpretation from the experimentally measured sample-to-sample fluctuations in the stress strain curves displaying large variations (around 50%) in the Young modulus.We have performed simulations of the model using the same parameters but different initial arrangements for the fibrils.When we do this we find much smaller variations of the resulting Young modulus (around 5%).A possible explanation is that the experimental variations reflect internal fluctuations in the cross-linking probability or in their stiffness, quantities that we have held constant in the simulations but have shown (Fig. 7) to have a crucial influence on the resulting Young modulus.
Our modeling strategy allows us to obtain an agreement with experimental results and suggests mechanisms for the observed non-linear elasticity.The current version of the computational model is presently not able to simulate fracture, since simply breaking overstretched cross-linkers or fibrils leads to discontinuities that are difficult to deal with in the present framework.A possible solution would involve quasi-static approaches where the network deformation is obtained by energy minimization as for other models for fracture [32].These types of models pave the way for a computationally based design of bio-inspired materials where the mechanical proper-ties can first be explored in silico to guide the production of materials with desired mechanical properties.

Figure 1 :
Figure 1: Echinoderm collagen networks.The fibril structure of collagen networks as observed by SEM are reported in panels a),c),e) for sea urchin (Paracentrotus lividus) (inset of b, photo by Federico Betti) starfish (Echinaster sepositus, (inset of d, photo by Tato Grasso, CC 3.0) and sea cucumber (Holothuria tubulosa (inset of f, photo by Roberto Pillon, CC 3.0).The corresponding diameter distributions are reported in panels b),d),f). 2µm

Figure 2 :Figure 3 :
Figure 2: Fracture surface of a collagen membrane.A fractured sea urchin collagen membrane as observed by SEM.Notice the very smooth fracture surface typical of brittle fracture.

Figure 4 :Figure 5 :
Figure 4: Stiffness and strength display sample-to-sample fluctuations.a)We observe considerable sample-to-sample variations in the Young modulus E * and the tensile strength σ c which appear to be correlated (correlation coefficient r = 0.73, linear least square fit yields a slope b = 0.12 ± 0.02 and intercept a = 6 ± 9 MPa).We also compute the cumulative distributions b) for Young moduli P(E * ) and c) tensile strengths P(σ c ). d) The strength distribution can be reasonably fit by the Weibull law.We observe no significant difference between sea urchin and sea cucumber derived membranes while membranes made out of starfish collagen appear to be stiffer and stronger.

Figure 6 :Figure 7 :Figure 8 :Figure 9 :
Figure 6: Fibrils reorient during deformation.The distribution of fibril angles θ as a function of strain for N = 12000 fibrils.The data indicate that the angles between between the loading axis and the fibrils decrease with strain.

Table 1 :
Parameters used in the numerical simulations in SI units.Notice that the parameter ranges (e.g. in k l ) represent the range that has been explored in different simulations.