An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves

The viscoelastic behavior of soft materials, especially cells and tissues, has been extensively investigated due to its importance in many biological and physiological processes that take place during development and even disease.*
Many techniques are used to quantify the mechanical properties of cells, among them micropipette aspiration, optical stretching, deformability cytometry and atomic force microscopy (AFM).*

The AFM, in particular, is still nowadays one of the most popular methods due to its conformity with various material types and geometries and the rather simple analysis process of the material properties.*

For a typical AFM indentation measurement, an AFM cantilever, with a distinct AFM tip shape, moves toward the sample with a predefined velocity and indents it until a prescribed force is reached. The AFM cantilever then moves upwards while detaching from the sample. The deflection and displacement signals of the AFM cantilever are processed further to extract the mechanical properties of the sample. Generally, a Hertzian model is fitted to the approach part of the force-indentation curves to quantify the apparent Young’s modulus.*

When applying the Hertzian model, few assumptions need to be considered, such as the material being homogeneous, isotropic, and linearly elastic. *

Cells and tissues, however, show not only elastic but also viscous behavior that is evident from the hysteresis between the approach and retraction segments of the force-indentation curve. Consequently, assessing this viscoelastic behavior is imperative for understanding the complex nature of biological matter.*

A number of studies utilized AFM to measure the viscoelastic properties of cells in both time and frequency domains.*

Ideally, to investigate the whole range of the viscoelastic behavior one needs to probe the material for a long time and observe its response or apply oscillatory signals and evaluate its phase lag. These approaches require the user to alter the probing method and add several steps to account for the time-dependent drift or the effect of the hydrodynamic drag of the surrounding medium. On top of that, in many of studies, the biological materials were probed with a linear approach followed by immediate retraction. The force-indentation curves from these studies were used to evaluate the apparent elastic modulus of the probed material using the standard Hertzian model. However, additional information concerning energy dissipation can still be extracted from the same curves to evaluate the viscoelasticity of the material.*

In the article “An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves”, Shada Abuhattum, Dominic Mokbel, Paul Müller, Despina Soteriou, Jochen Guck and Sebastian Aland propose a new fitting model to extract the viscoelastic properties of soft materials from AFM force-indentation curves. *

To construct the explicit relation of force and indentation, the authors first use a generalization of Maxwell and Kelvin-Voigt models to describe soft materials, and numerically simulate the indentation of such material with a spherical indenter. *

Shada Abuhattum et al. show that the proposed Kelvin-Voigt-Maxwell (KVM) model adequately captures the force-indentation curves of materials having different mechanical characteristics. *

Based on the simulation results, Shada Abuhattum et al. further propose an explicit force-indentation relation to be fitted to the force-indentation curves. This explicit relation simplifies the association of the mechanical properties with physically meaningful components and processes.
Finally, the authors apply the fitting model to a number of samples, including poroelastic and viscoelastic hydrogels as well as HeLa cells in two different cell cycle phases, interphase and mitotic. *

Shada Abuhattum et al. demonstrate that the distinct nature of the hydrogels, arising from the different crosslinking mechanisms, can be described with the fitting model. For the HeLa cells, the mitotic cells had a higher apparent elasticity and a lower apparent viscosity, implying a stiffer actin cortex and a diluted cytoplasm protein concentration, when compared with interphase cells.*

Their findings demonstrate that the proposed model can reliably extract viscoelastic properties from conventional force-indentation curves. Moreover, the model is able to assess the contribution of the different elastic and viscous elements, and thus allows a direct comparison between the viscoelastic nature of different materials.*

AFM measurements were preformed using a commercially available Atomic Force Microscope. To indent the samples, NanoWorld Pyrex-Nitride tipless AFM cantilevers PNP-TR-TL with a nominal spring constant of 0.08 mN/m were modified by gluing 5 μm diameter polystyrene beads to the underside of the AFM cantilevers using two component glue.*

The AFM cantilevers were calibrated prior to each experiment using the thermal noise method and their accurate spring constant ranged between 0.047-0.059 mN/m. For PAAm and agarose hydrogels, the AFM cantilever was lowered with a constant velocity (5, 10, or 15 μm/s) toward the surface of the sample until a force of 2 nN for agarose and 4 nN for PAAm was reached. These force set points accounted for an indentation in the range of 0.5–1 μm. For HeLa cells, the AFM cantilever was lowered with a constant velocity of 2 μm/s and the cells were indented until a force of 2 nN was reached, which accounted for an indentation depth in the range of 0.5–1.5 μm.*

Graphical abstract for the article “An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves” by Shada Abuhattum, Dominic Mokbel, Paul Müller, Despina Soteriou, Jochen Guck and Sebastian Aland consisting of 4 squares. showing a symbol for numerical simulations in the top left square, an arrow points to the bottom left square showing a graph and a formula as symbols for fitting algorithm a further arrow points to the bottom right square symbolizing the extraction of viscoelastic properties. The pictures in this square show on the left a drawing of the end of a tipless AFM cantilever on which a sphere is glued pressing on a cell, on the right of this picture there is another picture showing the end of a tipless AFM cantilever on which a sphere is glued pressing on a sphere or bead, underneath a graph symbolizing the mechanical properties of hydrogels is shown. Above this square on the top right a graph with a symbol for the mechanical behavior of the indented material is shown.NanoWorld Pyrex-Nitride tipless AFM cantilevers PNP-TR-TL with a nominal spring constant of 0.08 mN/m were modified by gluing 5 μm diameter polystyrene beads to the underside of the AFM cantilevers using two component glue were used for the atomic force microscopy indentation measurements described in the cited article.
Graphical abstract for the article “An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves” by Shada Abuhattum at al. 2022. NanoWorld Pyrex-Nitride tipless AFM cantilevers PNP-TR-TL with a nominal spring constant of 0.08 mN/m were modified by gluing 5 μm diameter polystyrene beads to the underside of the AFM cantilevers using two component glue were used for the atomic force microscopy indentation measurements described in the cited article.

*Shada Abuhattum, Dominic Mokbel, Paul Müller, Despina Soteriou, Jochen Guck and Sebastian Aland
An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves
iScience, Volume 25, ISSUE 4, 104016, April 15, 2022
DOI: https://doi.org/10.1016/j.isci.2022.104016

The article “An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves” by Shada Abuhattum, Dominic Mokbel, Paul Müller, Despina Soteriou, Jochen Guck and Sebastian Aland is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third-party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/.

Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer

Flux-closure structures, vortices/antivortices, skyrmions, and merons in oxides, metals and polymers represent non-trivial topologies in which a local polar/magnetic order undergoes quasi-continuous spatial variations in a host crystal lattice. These structures are now extensively studied due to emergent functionalities, but the application of electrical/mechanical fields has so far only served to destroy the polar topologies of interest. *

Topology created by quasi-continuous spatial variations of a local polarization direction represents an exotic state of matter, but field-driven manipulation has been hitherto limited to creation and destruction. *In the article “Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer” Mengfan Guo, Erxiang Xu, Houbing Huang, Changqing Guo, Hetian Chen, Shulin Chen, Shan He, Le Zhou, Jing Ma, Zhonghui Shen, Ben Xu, Di Yi, Peng Gao, Ce-Wen Nan, Neil. D. Mathur and Yang Shen report that relatively small electric or mechanical fields can drive the non-volatile rotation of polar spirals in discretized microregions of the relaxor ferroelectric polymer poly(vinylidene fluoride-ran-trifluoroethylene).*

These polar spirals arise from the asymmetric Coulomb interaction between vertically aligned helical polymer chains, and can be rotated in-plane through various angles with robust retention. *

Given also that their manipulation of topological order can be detected via infrared absorption, Mengfan Guo et al.’s work suggests a new direction for the application of complex materials. *

Mengfan Guo et al. produced a 100-nm-thick monolayer of face-on lamellae with vertically aligned polymer chains by melt-recrystallizing spin-coated thin films of P(VDF-TrFE).

The resulting melt-recrystallized thin film of the relaxor ferroelectric polymer was characterized by the authors using a commercial atomic force microscope for in-plane piezo-response force microscopy (IP-PFM).

NanoWorld Platinum Iridium coated Arrow-CONTPt AFM probes (typical resonant frequency: 14 kHz, typical force constant: 0.2 N/m, typical AFM tip radius 25 nm) were used for the in-plane (IP) PFM tests and the PFM lithography tests.

For piezo-response force microscopy (PFM) imaging, Vector mode was used where AFM tips were modulated at around 240 kHz for IP imaging, with the AC voltage set at 2 V. The images obtained by Vector Mode were double checked by using dual AC resonance tracking (DART) mode and the patterns could be reproduced. *

For angle-resolved IP-PFM tests, the rotation of sample was controlled by a protractor. To ensure identical position was imaged after rotating the sample, the authors made cross-scratches as a mark on the sample surface in advance. This method was applied to locate the scanning position in other situations if Mengfan Guo et al. had to move the sample in between the scanning probe microscopic studies.

For electric-field-induced manipulations using PFM lithography, the DC voltage on AFM tip was previously edited in the software. The scan speed was set at 1.95 Hz and no AC voltage was applied during the scanning. The DC voltage was divided by film thickness (100 nm) to obtain the electric field value. And an electric field with downward direction is defined with a positive sign.

For stress-induced manipulations, the deflection value of the PFM cantilever, which is a signal from photodetector, was preset to control the stress/force applied onto the sample. The difference in deflection value between a pressed AFM cantilever and a free AFM cantilever reflects how hard the AFM tip and sample surface are pressed to each other.*

To obtain the force value F, Mengfan Guo et al. first calibrated the AFM tips by the thermal noise method, and obtain the inverse optical lever sensitivity (InvOLS) and the spring constant k of the AFM tips.

The authors also conducted a polarization analysis based on their PFM measurements. *

To obtain the nominal toroidal order evaluated by the local curvature, the obtained IP-PFM amplitude image was firstly divided into 33 × 33 arrays, and each region was then subjected to a recognition of potential domain walls and measurement of an averaged curvature radius. *

To obtain polarization maps, angle-resolved IP-PFM images were first aligned to correct spatial distortion in nanoscale measurement. Positions with specific morphological characteristics were selected as reference points to determine the coordinate. After the correction, improved angle-resolved IP-PFM phase images would be divided into 64 × 64 arrays for deriving polarization maps. *

Fig. 1 from Mengfan Guo et al. (2024) “Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer”:Observation of a microregion containing an in-plane polar spiral. a Morphology of a melt-recrystallized thin film of the relaxor ferroelectric polymer. The scale bar is 2 μm. IP-PFM phase (b) and amplitude (c) images of the same area in a exhibiting concentric ring-shaped domains in curly stripe domains. d Distribution of domain wall curvatures in the same area in a–c evidencing nominal toroidal order. It is assumed that the local polarization is parallel to the nearest domain wall so that larger curvature (denoted red) reflects stronger toroidal order. IP-PFM phase images of identical concentric ring-shaped domains with the axis along vertical (e) and horizontal (f) measurement directions. The scale bar is 0.3 μm. The curl (g) and the divergence (h) of local polarization in the same area as e and f, revealing the polar spiral topology. i Schematic stereoscopic view of a CCW polar spiral, arrows represent regions of polarization. The red/blue arrows denote the polar source/sink that spirals in/out. The white arrows represent Néel rotation along the radial direction, as shown in more detail via the inset. NanoWorld Platinum Iridium coated Arrow-CONTPt AFM probes (typical resonant frequency: 14 kHz, typical force constant: 0.2 N/m, typical AFM tip radius 25 nm) were used for the in-plane (IP) PFM tests and the PFM lithography tests.
Fig. 1 from Mengfan Guo et al. (2024) “Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer”:
Observation of a microregion containing an in-plane polar spiral.
a Morphology of a melt-recrystallized thin film of the relaxor ferroelectric polymer. The scale bar is 2 μm. IP-PFM phase (b) and amplitude (c) images of the same area in a exhibiting concentric ring-shaped domains in curly stripe domains. d Distribution of domain wall curvatures in the same area in a–c evidencing nominal toroidal order. It is assumed that the local polarization is parallel to the nearest domain wall so that larger curvature (denoted red) reflects stronger toroidal order. IP-PFM phase images of identical concentric ring-shaped domains with the axis along vertical (e) and horizontal (f) measurement directions. The scale bar is 0.3 μm. The curl (g) and the divergence (h) of local polarization in the same area as e and f, revealing the polar spiral topology. i Schematic stereoscopic view of a CCW polar spiral, arrows represent regions of polarization. The red/blue arrows denote the polar source/sink that spirals in/out. The white arrows represent Néel rotation along the radial direction, as shown in more detail via the inset.

*Mengfan Guo, Erxiang Xu, Houbing Huang, Changqing Guo, Hetian Chen, Shulin Chen, Shan He, Le Zhou, Jing Ma, Zhonghui Shen, Ben Xu, Di Yi, Peng Gao, Ce-Wen Nan, Neil. D. Mathur and Yang Shen
Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer
Nature Communications volume 15, Article number: 348 (2024)
DOI: https://doi.org/10.1038/s41467-023-44395-5

The article “Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer” by Mengfan Guo, Erxiang Xu, Houbing Huang, Changqing Guo, Hetian Chen, Shulin Chen, Shan He, Le Zhou, Jing Ma, Zhonghui Shen, Ben Xu, Di Yi, Peng Gao, Ce-Wen Nan, Neil. D. Mathur and Yang Shen is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third-party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/.

Stiffness Mediated-Mechanosensations of Airway Smooth Muscle Cells on Linear Stiffness Gradient Hydrogels

Airflow limitation in obstructive airway disease is characterized by narrowing of the airway lumen from excessive contraction of airway smooth muscle (ASM) and remodeling of the airway wall which includes changes in the extracellular matrix (ECM) of the ASM layer.*

Previous studies on human airway smooth muscle cells ( hASMC ) have independently assessed the influence of extracellular matrix (ECM) proteins on substrates of supra-physiological stiffnesses, such as tissue culture plastic or glass.*

While the influence of discrete substrate stiffness on hASMC behavior has been examined, manipulation of both substrate stiffness and ECM proteins simultaneously (as expected in disease) has not been extensively modeled in vitro.*

In the article “Stiffness Mediated-Mechanosensation of Airway Smooth Muscle Cells on Linear Stiffness Gradient Hydrogels” Yong Hwee Tan, Kimberley C. W. Wang, Ian L. Chin, Rowan W. Sanderson, Jiayue Li, Brendan F. Kennedy, Peter B. Noble and Yu Suk Choi highlight the interplay and complexities between stiffness and ECM protein type on hASMC mechanosensation, relevant to airway remodelling in obstructive airway diseases.*

The authors first determined a physiological range of ASM layer stiffness using a porcine airway and used these empirical recordings to inform the fabrication of a linear stiffness gradient platform coated with different ECM proteins.*

Using this linear stiffness gradient platform, Yong Hwee Tan et al. profiled hASMC morphology, contractile function with alpha-smooth muscle actin (αSMA) and mechanisms of mechanosensation, specifically with nuclear translocalization of Yes-associated protein (YAP) and lamin-A expression.*

Yong Hwee Tan et al.’s assessment of hASMC mechanosensation utilized an innovative hydrogel platform delivering a linear stiffness gradient to understand stiffness-mediated cell behavior with an ECM substrate for cellular adhesion. *

The employment of a stiffness gradient that was designed after empirical measurements performed on ex vivo ASM tissue, enabled the presentation of physiologically relevant stiffnesses to study hASMC behavior.*

Using this platform, the authors of the article found that hASMC mechanosense underlying mechanical cues more than the types of proteins they are anchored to by screening hASMC morphology, contractile phenotype, and mechanomarker expression, with a few exceptions.*

While the authors acknowledge that the findings from their study were done using cells from only one donor they still think that their study provides a proof of concept for the relevance of hASMC mechanosensation to ECM stiffness, and is another step in the right direction for understanding the pathophysiological impact of airway remodeling in obstructive diseases and exploring potential avenues for improving therapy through greater fidelity of in vitro platforms that include key concepts of mechanosensation. *

Yong Hwee Tan et al. wanted to use the same method which is used to assess hydrogel stiffness, namely atomic force microscopy (AFM), to measure ASM stiffness.*

However, nanoscale measurements of ASM strips by AFM proved to be difficult due to an uneven tissue surface after de-epithelialization (Figure S1C, Supporting Information of the cited article), resulting in false force triggering.

To validate the translation of stiffness values measured from macroscale compression (ASM strips) to nanoscale indentation (AFM on hydrogels), Yong Hwee Tan et al. fabricated additional hydrogels of four different stiffnesses using well-characterized polyacrylamide and compared the stiffness of hydrogels measured by uniaxial compression tester and atomic force microscopy (Figure S2A, Supporting Information of the cited article).

The nanoscale stiffness of hydrogels was assessed using an atomic force microscope (AFM) with NanoWorld triangular Pyrex-Nitride PNP-TR AFM probes (the longer AFM cantilever beam – CB 2 – with 200 µm length was used).

These AFM tips probed hydrogels immersed in 1 × PBS with 2 nN, an approach velocity of 2 µm s−1 and a retraction velocity of 10 µm s−1.

Young’s modulus was determined from linear portions of contact-generated force curves using a custom-written code in Igor Pro.

All probe indentations were made in triplicate and averaged for a stiffness measurement in kilopascals (kPa).

An example force curve is shown in Figure S2B, Supporting Information of the cited article. Validation of a linear stiffness gradient was achieved with eight indentations on the hydrogel, 2 mm away from both edges of the hydrogel and at 1 mm intervals along the stiffness gradient axis. Measurements were plotted against displacement from the hydrogel edge (soft to stiff) (Figure 2B of the cited article).

Figure 2 from Yong Hwee Tan et al. 2024 “Stiffness Mediated-Mechanosensation of Airway Smooth Muscle Cells on Linear Stiffness Gradient Hydrogels”:Linear stiffness gradient hydrogel fabrication. A) A schematic of a two-step polymerization process. i) 120 µL of mixed polyacrylamide (PA) solution (% acrylamide + % bis-acrylamide) was added to the primary mold and left to polymerize under ii) a methacrylated coverslip for 20 min. iii) Wedge-shaped 1° gel was removed and flipped for placement of a iv) secondary mold before v) addition of a second 120 µL PA solution and polymerized under a vi) dichlorodimethylsilane-coated coverslip for 20 min. vii) Removal of coverslip and mold completes the fabrication of bi-layered stiffness gradient hydrogel. viii) the dotted arrow indicating the direction of gradient and atomic force microscopy (AFM) measurement. B) Young's moduli gradient measured by AFM. Twelve hydrogels were selected (one gel per batch) and assessed for stiffness, yielding a gradient of 4.0 kPa mm−1, with a range of 1.7 ± 1.2 to 29.6 ± 4.3 kPa (R2 = 0.998, n = 8). Data are presented as mean ± SD. NanoWorld triangular Pyrex-Nitride PNP-TR AFM probes were used to assess the stiffness of the hydrogels with atomic force microscopy.
Figure 2 from Yong Hwee Tan et al. 2024 “Stiffness Mediated-Mechanosensation of Airway Smooth Muscle Cells on Linear Stiffness Gradient Hydrogels”:
Linear stiffness gradient hydrogel fabrication. A) A schematic of a two-step polymerization process. i) 120 µL of mixed polyacrylamide (PA) solution (% acrylamide + % bis-acrylamide) was added to the primary mold and left to polymerize under ii) a methacrylated coverslip for 20 min. iii) Wedge-shaped 1° gel was removed and flipped for placement of a iv) secondary mold before v) addition of a second 120 µL PA solution and polymerized under a vi) dichlorodimethylsilane-coated coverslip for 20 min. vii) Removal of coverslip and mold completes the fabrication of bi-layered stiffness gradient hydrogel. viii) the dotted arrow indicating the direction of gradient and atomic force microscopy (AFM) measurement. B) Young’s moduli gradient measured by AFM. Twelve hydrogels were selected (one gel per batch) and assessed for stiffness, yielding a gradient of 4.0 kPa mm−1, with a range of 1.7 ± 1.2 to 29.6 ± 4.3 kPa (R2 = 0.998, n = 8). Data are presented as mean ± SD.

 

Figure S2 from Yong Hwee Tan et al. 2024 “Stiffness Mediated-Mechanosensation of Airway Smooth Muscle Cells on Linear Stiffness Gradient Hydrogels”:(A) Correlation of Young’s modulus from macroscale stiffness (UCT) assessment with nanoindentation (AFM), was conducted using cylindrical PA hydrogels of different Acrylamide %/Bis-acrylamide % derived from Tse and Engler [47] 10 %/0.06 %, 10 %/0.1 %, 10 %/0.15 % and 10 %/0.3 % (Linear regression, P < 0.0001, R2 = 0.9288, n = 4). Data are presented as mean  SEM. (B) An example force curve from atomic force microscopy with an approach velocity of 2 μm/s, until a 2 nN trigger force was registered, and retraction of indenter at 10 μm/s. NanoWorld triangular Pyrex-Nitride PNP-TR AFM probes were used to assess the stiffness of the hydrogels with atomic force microscopy.
Figure S2 from Yong Hwee Tan et al. 2024 “Stiffness Mediated-Mechanosensation of Airway Smooth Muscle Cells on Linear Stiffness Gradient Hydrogels”:
(A) Correlation of Young’s modulus from macroscale stiffness (UCT) assessment with nanoindentation (AFM), was conducted using cylindrical PA hydrogels of different Acrylamide %/Bis-acrylamide % derived from Tse and Engler [47] 10 %/0.06 %, 10 %/0.1 %, 10 %/0.15 % and 10 %/0.3 % (Linear regression, P < 0.0001, R2 = 0.9288, n = 4). Data are presented as mean  SEM. (B) An example force curve from atomic force microscopy with an approach velocity of 2 μm/s, until a 2 nN trigger force was registered, and retraction of indenter at 10 μm/s.
*Yong Hwee Tan, Kimberley C. W. Wang, Ian L. Chin, Rowan W. Sanderson, Jiayue Li, Brendan F. Kennedy, Peter B. Noble and Yu Suk Choi
Stiffness Mediated-Mechanosensation of Airway Smooth Muscle Cells on Linear Stiffness Gradient Hydrogels
Advanced Healthcare Materials 2024, 2304254
DOI: https://doi.org/10.1002/adhm.202304254

The article “Stiffness Mediated-Mechanosensation of Airway Smooth Muscle Cells on Linear Stiffness Gradient Hydrogels” by Yong Hwee Tan, Kimberley C. W. Wang, Ian L. Chin, Rowan W. Sanderson, Jiayue Li, Brendan F. Kennedy, Peter B. Noble and Yu Suk Choi is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third-party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/.