Modulating extracellular matrix (ECM) elasticity with fibrillar collagen offers great potential for regenerative medicine, drug discovery, and disease modeling by replicating in vivo mechanical signals. This enhances the understanding of cellular responses and fosters therapeutic innovation. However, precise ECM elasticity measurements are still lacking. This study couples time-resolved Brillouin spectroscopy and pulsed laser-induced Scholte wave generation. We measure how collagen fibrillation affects sound velocity and refractive index. These insights are advancing tissue engineering and cellular biomechanics.

Collagen, making up 30% of the body’s proteins, is prevalent in bones, muscles, blood, and skin, where it constitutes 75%. It is crucial for skin elasticity and serves as an ideal matrix for joints, tendons, and ligaments. Understanding collagen fibrillation kinetics is essential in tissue engineering and cell culture, where optimal conditions for cell growth and function are paramount. Manipulating the mechanical properties of the extracellular matrix (ECM) is crucial as ECM elasticity significantly influences cellular processes, such as adhesion, migration, proliferation, and differentiation. Fibrillar collagen, with its adjustable mechanical properties, stands out due to its ability to influence cellular behavior and tissue development. By controlling collagen concentration, cross-linking density, and fiber alignment, researchers can tailor ECM stiffness and elasticity to mimic various tissue microenvironments, enhancing the relevance of in vitro models.1–3 This modulation has significant implications for regenerative medicine, drug discovery, and disease modeling, aiding in understanding cellular responses to different conditions.

Despite its potential, ECM elasticity’s quantitative measurements remain incomplete. Brillouin spectroscopy, which is particularly well suited for investigations of semitransparent systems, has already been applied to the elastic characterization in numerous biological studies.1–3 In this article, we report on alternative experiments on the thermoelastic excitation of acoustic waves induced using a pulsed laser in a collagen solution. Simultaneous detection of longitudinal waves and Scholte waves at the collagen/solid interface enables us to quantitatively measure the effect of the fibrillation process on sound velocity and refractive index vs the fibrillation time, providing insights into the relationship between collagen fibrillation and ECM mechanical properties.

Collagen fibrillation was performed directly on the supports, which were 1 × 1 cm2 silicon wafers, for variable durations ranging from 30 minutes to 72 h. The collagen utilized was concentrated to 20 mg/ml, specifically SYMATESE collagen, which contained 17 mM acetic acid. For each sample, a mixture of 0.387 ml collagen, 0.018 ml sodium hydroxide, and 0.044 ml phosphatebuffered saline (PBS) was prepared. The mixture was stirred to ensure uniform distribution of sodium hydroxide and PBS within the collagen. Subsequently, the dish was resealed and placed in the culture room’s incubator for the required fibrillation duration.

Time resolved pump and probe experiments have been realized using a mode-locked Ti:sapphire (MAI-TAI Spectra) laser source operated at 800 nm with a pulse duration below 100 fs at the laser output and a repetition rate of 79.4 MHz. Synchronous interferometric detection is performed by modulating the pump beam at 1.8 MHz using an acousto-optic modulator. A Michelson scheme allowed measuring the perpendicular surface displacement. The pump and probe are first superimposed on the acoustic transducer to study longitudinal acoustic waves. To obtain information on surface or interface waves, the probe beam is swept around the pump excitation using a 4f configuration. In both configurations, we performed a two-color pump–probe experiment by doubling the frequency of the probe or pump beam. As the fibrillation process starts at the free surface of the collagen, we used a sufficiently thin metallic transducer, deposited on a sapphire or glass substrate, placed in contact with this free surface. The titanium or aluminum layer is thick enough to generate an acoustic burst and thin enough to allow the probe beam to follow the propagation of the acoustic pulse in the first few micrometers of the biomaterial layer.

First, let us consider longitudinal acoustic wave propagating perpendicular to the surface, using the experimental geometry described in Fig. 1(a). In this case, the advantage of sapphire, compared with a glass substrate for example, lies in the fact that the Brillouin signature of this transparent compound is in a much higher frequency range, around 100 GHz, than that expected in the collagen gel. In Fig. 1(b), reflectivity variations are measured with variable pump and probe delay. After subtraction of the photothermal background, a double Brillouin oscillation can be clearly observed. The high-frequency signature comes from the acoustic wave transmitted into the sapphire, while the low-frequency component corresponds to the collagen layer signature. In this case, the transducer should be pressed sufficiently on the collagen layer to ensure good acoustic transfer.

FIG. 1.

(a) Experimental geometry. (b) Time resolved measurement. Bottom: Raw data with a collagen layer with a fibrillation time of 24 h. Middle: Signal showing the Brillouin oscillation without the background. Top: Brillouin oscillation with a fibrillation time of 48 h.

FIG. 1.

(a) Experimental geometry. (b) Time resolved measurement. Bottom: Raw data with a collagen layer with a fibrillation time of 24 h. Middle: Signal showing the Brillouin oscillation without the background. Top: Brillouin oscillation with a fibrillation time of 48 h.

Close modal
The lower frequency component evolves as a function of fibrillation time. The evolution of Brillouin frequency vs fibrillation time is shown in Fig. 2. Over a 72 h period, the frequency changes by almost 26%. This frequency is linked to the layer’s elastic properties by the following formula:
(1)
where λ is the probe wavelength, here 400 nm; V is the longitudinal velocity; and n is the refractive index.
FIG. 2.

Variation of Brillouin frequency in collagen 20 mg/ml vs fibrillation time.

FIG. 2.

Variation of Brillouin frequency in collagen 20 mg/ml vs fibrillation time.

Close modal

Consequently, knowledge of the refractive index is mandatory to determine longitudinal velocity. This parameter varies considerably in the literature and is also likely to change with fibrillation time. To overcome this limitation, we used an alternative geometry to directly measure the wave velocity in the layer, independently of the optical properties. While our previous focus was on longitudinal volume waves, we are now concentrating on waves that may exist at the interface between the metal transducer and the collagen gel.

Following this, the probe beam is scanned around the pump epicenter to be sensitive to in-plane propagation in the collagen layer. The experimental setup is briefly illustrated in Fig. 3(a). The mapping of surface displacements when the collagen layer is simply replaced by air is shown in Fig. 3(b). A set of concentric circles corresponding to the Rayleigh waves emitted by the various pump pulses can be observed.4 

FIG. 3.

(a) Experimental geometry. The Al transducer thickness is chosen to have enough acoustic emission in the plane and enough probe transmission to detect interface waves. (b) Rayleigh wave mapping at the free surface. (c) Scholte wave mapping at the interface Al–collagen. A time-resolved measurement of the Scholte wave shows a central frequency of around 500 MHz.

FIG. 3.

(a) Experimental geometry. The Al transducer thickness is chosen to have enough acoustic emission in the plane and enough probe transmission to detect interface waves. (b) Rayleigh wave mapping at the free surface. (c) Scholte wave mapping at the interface Al–collagen. A time-resolved measurement of the Scholte wave shows a central frequency of around 500 MHz.

Close modal

The isotropic nature within the plane is associated with the isotropy of the glass substrate used here. Figure 3(c) shows a similar situation when the aluminum transducer is in contact with non-fibrillated collagen. In this case, we note the presence of similar features that clearly exhibited a much lower propagation velocity. By pointing out the position of these rings and knowing the repetition rate of our laser, 12.59 ns in this case, we obtain Fig. 4, enabling us to extract the propagation velocity of these different waves.

FIG. 4.

Position of the wave front at different times with wave velocity deduced. Black: data for a Rayleigh wave. Red: data for a Scholte wave in collagen. Blue: data for a Scholte wave in ethanol.

FIG. 4.

Position of the wave front at different times with wave velocity deduced. Black: data for a Rayleigh wave. Red: data for a Scholte wave in collagen. Blue: data for a Scholte wave in ethanol.

Close modal

For an air layer, a velocity of 3020 m/s is obtained, which is in perfect agreement with a Rayleigh wave propagating in a glass substrate.5 Here, the very thin aluminum layer has no influence on this propagation due to the low frequency of the surface acoustic wave involved.

In the case of ethanol or collagen, we obtain much lower velocities: 1118 m/s and 1507 m/s, respectively, indicating that we are dealing with a different type of wave. In this hard solid/fluid configuration,6,7 where the transducer has a much higher longitudinal velocity than the liquid, it is well known that the Rayleigh wave is transformed into a strongly attenuated leaky Rayleigh wave; additionally, a Scholte wave appears at the interface, capable of propagating over long distances. The velocity vs of this wave is given by the following equation:8 
(2)
where v and vt are the longitudinal and transverse velocities, respectively; ρ is the density; and indices 1 and 2 refer to solid and liquid, respectively.

For the glass–ethanol interface, the Scholte wave velocity is expected to be very close to the longitudinal wave velocity of ethanol, equal to 1150 m/s, which is in good agreement with the experimental speed deduced from Fig. 4, demonstrating that the concentric rings depicted in Fig. 3(b) may be associated with such Scholte waves. Scholte waves induced by laser excitation have already been demonstrated in the lower frequency range.9 Conversely, by introducing the Scholte velocity measured in a collagen solution into Eq. (2), we can deduce a longitudinal velocity of 1512 m s−1 for the non-fibrillated collagen using a density equal to 1. Then, considering this deduced longitudinal velocity, and the Brillouin frequency measured in Fig. 2, we obtain a value for refractive index at 400 nm equal to 1.35. This new approach, which enables the refractive index to be determined using purely acoustic measurements, sheds new light on issues relating to the elasticity of living cells, where this parameter has hitherto been neglected due to a lack of quantitative measurements.10–14 However, it should be pointed out here that this approach neglects possible viscoelastic effects that may affect the propagation velocity as a function of the frequency range considered, which here is at least one order of magnitude. As a result, the underestimation of velocity in the Brillouin frequency may lead to a slight overestimate in the optical index.

By reproducing the same acoustic mappings for fibrillated solutions, we can extract the evolution of longitudinal velocity and refractive index as a function of fibrillation time, without having to rely on the knowledge of additional parameters, except density. These evolutions are shown in Fig. 5. Consequently, the fluctuation in Brillouin frequency observed in Fig. 2 was, in fact, essentially linked to the variation in collagen elasticity, with the optical index varying by only around 3%. Now, if the density is assumed to be roughly constant, equal to 1, then the bulk modulus given by B = ρV2 may change by 40% varying from 2.5 GPa up to 3.5 GPa. It should, nevertheless, be kept in mind that the treatment proposed here is based on the assumption that the longitudinal velocities in the plane and perpendicular to the layer are equal, and potential anisotropies in the propagation properties are neglected.

FIG. 5.

Deduced longitudinal velocity (in black) and optical index (in red) vs fibrillation time.

FIG. 5.

Deduced longitudinal velocity (in black) and optical index (in red) vs fibrillation time.

Close modal

We have developed a new technique to determine the longitudinal velocities and refractive index in materials whose elasticity can be modified on demand, which is particularly relevant for biological systems such as biomaterials. This method has shown remarkable efficiency in studying cell growth, allowing for variations of up to 40% in the bulk elastic modulus. The application of this technique opens new perspectives for the analysis and engineering of biomaterials, providing precise and adaptable means to optimize mechanical properties according to the specific needs of biological studies. However, a major challenge remains: the precise determination of material’s density, a crucial step for complete characterization of its mechanical properties; interested readers may find some pointers in Ref. 13. Continuing research in this direction could not only enhance our understanding of biomaterials but also lead to significant advancements in their practical applications.

The authors have no conflicts to disclose.

A. Hamraoui: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal). O. Sénépart: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal). L. Belliard: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal).

The data that support the findings of this study are available within the article.

1.
B.
Cortese
,
I.
Elena Palamà
,
S.
D’Amone
, and
G.
Gigli
, “
Influence of electrotaxis on cell behaviour
,”
Integr. Biol.
6
(
9
),
817
830
(
2014
).
2.
K.
Hamid
,
D.
Cikes
,
M.
Samalova
,
L.
Schindler
,
L.-M.
Appel
,
M.
Urbanek
,
I.
Yudushkin
,
D.
Slade
,
W. J.
Weninger
,
A.
Peaucelle
,
J.
Penninger
, and
K.
Elsayad
, “
Brillouin light scattering anisotropy microscopy for imaging the viscoelastic anisotropy in living cells
,”
Nat. Photonics
18
(
3
),
276
285
(
2024
).
3.
K.
Elsayad
,
S.
Werner
,
M.
Gallemí
,
J.
Kong
,
E. R. S.
Guajardo
,
L.
Zhang
,
Y.
Jaillais
,
T.
Greb
, and
Y.
Belkhadir
, “
Mapping the subcellular mechanical properties of live cells in tissues with fluorescence emission–Brillouin imaging
,”
Sci. Signaling
9
(
435
),
rs5
(
2016
).
4.
D. H.
Hurley
,
O. B.
Wright
,
O.
Matsuda
,
T.
Suzuki
,
S.
Tamura
, and
Y.
Sugawara
, “
Time-resolved surface acoustic wave propagation across a single grain boundary
,”
Phys. Rev. B
73
,
125403
(
2006
).
5.
P. A.
Mante
,
J. F.
Robillard
, and
A.
Devos
, “
Complete thin film mechanical characterization using picosecond ultrasonics and nanostructured transducers: Experimental demonstration on SiO2
,”
Appl. Phys. Lett.
93
(
7
),
071909
(
2008
).
6.
C.
Glorieux
,
K.
Van De Rostyne
,
K.
Nelson
,
W.
Gao
,
W.
Lauriks
, and
J.
Thoen
, “
On the character of acoustic waves at the interface between hard and soft solids and liquids
,”
J. Acoust. Soc. Am.
110
(
3
),
1299
1306
(
2001
).
7.
V.
Gusev
,
C.
Desmet
,
W.
Lauriks
,
C.
Glorieux
, and
J.
Thoen
, “
Theory of Scholte, leaky Rayleigh, and lateral wave excitation via the laser-induced thermoelastic effect
,”
J. Acoust. Soc. Am.
100
(
3
),
1514
1528
(
1996
).
8.
L. M.
Brekhovskikh
,“
Waves in layered media
,” in
Applied Mathematics and Mechanics
(
Academic Press
,
1980
).
9.
C.
Desmet
,
V.
Gusev
,
W.
Lauriks
,
C.
Glorieux
, and
J.
Thoen
, “
Laser-induced thermoelastic excitation of Scholte waves
,”
Appl. Phys. Lett.
68
(
21
),
2939
2941
(
1996
).
10.
S.
Danworaphong
,
M.
Tomoda
,
Y.
Matsumoto
,
O.
Matsuda
,
T.
Ohashi
,
H.
Watanabe
,
M.
Nagayama
,
K.
Gohara
,
P. H.
Otsuka
, and
O. B.
Wright
, “
Three-dimensional imaging of biological cells with picosecond ultrasonics
,”
Appl. Phys. Lett.
106
(
16
),
163701
(
2015
).
11.
L.
Liu
,
A.
Viel
,
G.
Le Saux
,
L.
Plawinski
,
G.
Muggiolu
,
P.
Barberet
,
M.
Pereira
,
C.
Ayela
,
H.
Seznec
,
M.-C.
Durrieu
,
J.-M.
Olive
, and
B.
Audoin
, “
Remote imaging of single cell 3D morphology with ultrafast coherent phonons and their resonance harmonics
,”
Sci. Rep.
9
(
1
),
6409
(
2019
).
12.
F.
Pérez-Cota
,
R. J.
Smith
,
H. M.
Elsheikha
, and
M.
Clark
, “
New insights into the mechanical properties of Acanthamoeba castellanii cysts as revealed by phonon microscopy
,”
Biomed. Opt. Express
10
(
5
),
2399
(
2019
).
13.
A.
Viel
,
E.
Péronne
,
O.
Sénépart
,
L.
Becerra
,
C.
Legay
,
F.
Semprez
,
L.
Trichet
,
T.
Coradin
,
H.
Ahmed
, and
B.
Laurent
, “
Picosecond ultrasounds as elasticity probes in neuron-like cells models
,”
Appl. Phys. Lett.
115
(
21
),
213701
(
2019
).
14.
L.
Liu
,
L.
Plawinski
,
M.-C.
Durrieu
, and
A.
Bertrand
, “
Label-free multi-parametric imaging of single cells: Dual picosecond optoacoustic microscopy
,”
J. Biophot.
12
(
8
),
e201900045
(
2019
).