Investigating metalloproteinases MMP-2 and MMP-9 mechanosensitivity to feedback loops involved in the regulation of in vitro angiogenesis by endogenous mechanical stresses.

Angiogenesis is a complex morphogenetic process regulated by growth factors, but also by the force balance between endothelial cells (EC) traction stresses and extracellular matrix (ECM) viscoelastic resistance. Studies conducted with in vitro angiogenesis assays demonstrated that decreasing ECM stiffness triggers an angiogenic switch that promotes organization of EC into tubular cords or pseudo-capillaries. Thus, mechano-sensitivity of EC with regard to proteases secretion, and notably matrix metalloproteinases (MMPs), should likely play a pivotal role in this switching mechanism. While most studies analysing strain regulation of MMPs used cell cultured on stretched membranes, this work focuses on MMP expression during self-assembly of EC into capillary-like structures within fibrin gels, i.e. on conditions that mimics more closely the in vivo cellular mechanical microenvironment. The activity of MMP-2 and MMP-9, two MMPs that have a pivotal role in capillaries formation, has been monitored in pace with the progressive elongation of EAhy926 cells that takes place during the emergence of cellular cords. We found an increase of the zymogen proMMP-2 that correlates with the initial stages of EC cords formation. However, MMP-2 was not detected. ProMMP-9 secretion decreased, with levels of MMP-9 kept at a rather low value. In order to analyse more precisely the observed differences of EAhy926 response on fibrin and plastic substrates, we proposed a theoretical model of the mechano-regulation of proMMP-2 activation in the presence of type 2 tissue inhibitor of MMPs (TIMP-2). Using association/dissociation rates experimentally reported for this enzymatic network, the model adequately describes the synergism of proMMP-2 and TIMP-2 strain activation during pseudo-capillary morphogenesis. All together, these results provide a first step toward a systems biology approach of angiogenesis mechano-regulation by cell-generated extracellular stresses and strains.

Elucidating atherosclerotic vulnerable plaque rupture by modeling cross substitution of ApoE-/- mouse and human plaque components stiffnesses.

The structure of mouse atherosclerotic lesions may differ from that of humans, and mouse atherosclerotic plaques do not rupture except in some specific locations such as the brachiocephalic artery. Recently, our group was the first to observe that the amplitudes of in vivo stresses in ApoE-/- mouse aortic atherosclerotic lesions were much lower and differed from those found in a previous work performed on human lesions. In this previous preliminary work, we hypothesized that the plaque mechanical properties (MP) may in turn be responsible for such species differences. However, the limited number of human samples used in our previous comparative study was relevant but not sufficient to broadly validate such hypothesis. Therefore, in this study, we propose an original finite element strategy that reconstructs the in vivo stress/strain (IVS/S) distributions in ApoE-/- artherosclerotic vessels based on cross substitution of ApoE-/- mouse and human plaque components stiffnesses and including residual stress/strain (RS/S). Our results: (1) showed that including RS/S decreases by a factor 2 the amplitude of maximal IVS/S, and more importantly, (2) demonstrated that the MP of the ApoE-/- plaque constituents are mainly responsible for the low level-compared with human-of intraplaque stress in ApoE-/- mouse aortic atherosclerotic lesions (8.36 ± 2.63 kPa vs. 182.25 ± 55.88 kPa for human). Our study highlights that such differences in the distribution and amplitude of vessel wall stress might be one key feature for explaining for the difference in lesion stability between human coronary and mouse aortic lesions.

Robustness analysis and behavior discrimination in enzymatic reaction networks.

Characterizing the behavior and robustness of enzymatic networks with numerous variables and unknown parameter values is a major challenge in biology, especially when some enzymes have counter-intuitive properties or switch-like behavior between activation and inhibition. In this paper, we propose new methodological and tool-supported contributions, based on the intuitive formalism of temporal logic, to express in a rigorous manner arbitrarily complex dynamical properties. Our multi-step analysis allows efficient sampling of the parameter space in order to define feasible regions in which the model exhibits imposed or experimentally observed behaviors. In a first step, an algorithmic methodology involving sensitivity analysis is conducted to determine bifurcation thresholds for a limited number of model parameters or initial conditions. In a second step, this boundary detection is supplemented by a global robustness analysis, based on quasi-Monte Carlo approach that takes into account all model parameters. We apply this method to a well-documented enzymatic reaction network describing collagen proteolysis by matrix metalloproteinase MMP2 and membrane type 1 metalloproteinase (MT1-MMP) in the presence of tissue inhibitor of metalloproteinase TIMP2. For this model, our method provides an extended analysis and quantification of network robustness toward paradoxical TIMP2 switching activity between activation or inhibition of MMP2 production. Further implication of our approach is illustrated by demonstrating and analyzing the possible existence of oscillatory behaviors when considering an extended open configuration of the enzymatic network. Notably, we construct bifurcation diagrams that specify key parameters values controlling the co-existence of stable steady and non-steady oscillatory proteolytic dynamics.

Is arterial wall-strain stiffening an additional process responsible for atherosclerosis in coronary bifurcations?: an in vivo study based on dynamic CT and MRI.

Coronary bifurcations represent specific regions of the arterial tree that are susceptible to atherosclerotic lesions. While the effects of vessel compliance, curvature, pulsatile blood flow, and cardiac motion on coronary endothelial shear stress have been widely explored, the effects of myocardial contraction on arterial wall stress/strain (WS/S) and vessel stiffness distributions remain unclear. Local increase of vessel stiffness resulting from wall-strain stiffening phenomenon (a local process due to the nonlinear mechanical properties of the arterial wall) may be critical in the development of atherosclerotic lesions. Therefore, the aim of this study was to quantify WS/S and stiffness in coronary bifurcations and to investigate correlations with plaque sites. Anatomic coronary geometry and cardiac motion were generated based on both computed tomography and MRI examinations of eight patients with minimal coronary disease. Computational structural analyses using the finite element method were subsequently performed, and spatial luminal arterial wall stretch (LW(Stretch)) and stiffness (LW(Stiff)) distributions in the left main coronary bifurcations were calculated. Our results show that all plaque sites were concomitantly subject to high LW(Stretch) and high LW(Stiff), with mean amplitudes of 34.7 ± 1.6% and 442.4 ± 113.0 kPa, respectively. The mean LW(Stiff) amplitude was found slightly greater at the plaque sites on the left main coronary artery (mean value: 482.2 ± 88.1 kPa) compared with those computed on the left anterior descending and left circumflex coronary arteries (416.3 ± 61.5 and 428.7 ± 181.8 kPa, respectively). These findings suggest that local wall stiffness plays a role in the initiation of atherosclerotic lesions.

Assessing low levels of mechanical stress in aortic atherosclerotic lesions from apolipoprotein E-/- mice--brief report.

OBJECTIVE: Despite the fact that mechanical stresses are well recognized as key determinants for atherosclerotic plaque rupture, very little is known about stress amplitude and distribution in atherosclerotic lesions, even in the standard apolipoprotein E (apoE)-/- mouse model of atherosclerosis. Our objectives were to combine immunohistology, atomic force microscopy measurements, and finite element computational analysis for the accurate quantification of stress amplitude and distribution in apoE-/- mouse aortic atherosclerotic lesions. METHODS AND RESULTS: Residual stresses and strains were released by radially cutting aortic arch segments from 7- to 30-week-old pathological apoE-/- (n=25) and healthy control mice (n=20). Immunohistology, atomic force microscopy, and biomechanical modeling taking into account regional residual stresses and strains were performed. Maximum stress values were observed in the normal arterial wall (276±71 kPa), whereas low values (<20 kPa) were observed in all plaque areas. Stress distribution was not correlated to macrophage infiltration. CONCLUSIONS: Low mechanical stress amplitude was observed in apoE-/- mouse aortic atherosclerotic lesions. This original study provides a basis for further investigations aimed at determining whether low stress levels are responsible for the apparently higher stability of murine aortic atherosclerotic lesions.

Mapping elasticity moduli of atherosclerotic plaque in situ via atomic force microscopy.

Several studies have suggested that evolving mechanical stresses and strains drive atherosclerotic plaque development and vulnerability. Especially, stress distribution in the plaque fibrous capsule is an important determinant for the risk of vulnerable plaque rupture. Knowledge of the stiffness of atherosclerotic plaque components is therefore of critical importance. In this work, force mapping experiments using atomic force microscopy (AFM) were conducted in apolipoprotein E-deficient (ApoE(-/-)) mouse, which represents the most widely used experimental model for studying mechanisms underlying the development of atherosclerotic lesions. To obtain the elastic material properties of fibrous caps and lipidic cores of atherosclerotic plaques, serial cross-sections of aortic arch lesions were probed at different sites. Atherosclerotic plaque sub-structures were subdivided into cellular fibrotic, hypocellular fibrotic and lipidic rich areas according to histological staining. Hertz's contact mechanics were used to determine elasticity (Young's) moduli that were related to the underlying histological plaque structure. Cellular fibrotic regions exhibit a mean Young modulus of 10.4±5.7kPa. Hypocellular fibrous caps were almost six-times stiffer, with average modulus value of 59.4±47.4kPa, locally rising up to ∼250kPa. Lipid rich areas exhibit a rather large range of Young's moduli, with average value of 5.5±3.5kPa. Such precise quantification of plaque stiffness heterogeneity will allow investigators to have prospectively a better monitoring of atherosclerotic disease evolution, including arterial wall remodeling and plaque rupture, in response to mechanical constraints imposed by vascular shear stress and blood pressure.

On the potential of a new IVUS elasticity modulus imaging approach for detecting vulnerable atherosclerotic coronary plaques: in vitro vessel phantom study.

Peak cap stress amplitude is recognized as a good indicator of vulnerable plaque (VP) rupture. However, such stress evaluation strongly relies on a precise, but still lacking, knowledge of the mechanical properties exhibited by the plaque components. As a first response to this limitation, our group recently developed, in a previous theoretical study, an original approach, called iMOD (imaging modulography), which reconstructs elasticity maps (or modulograms) of atheroma plaques from the estimation of strain fields. In the present in vitro experimental study, conducted on polyvinyl alcohol cryogel arterial phantoms, we investigate the benefit of coupling the iMOD procedure with the acquisition of intravascular ultrasound (IVUS) measurements for detection of VP. Our results show that the combined iMOD-IVUS strategy: (1) successfully detected and quantified soft inclusion contours with high positive predictive and sensitivity values of 89.7 ± 3.9% and 81.5 ± 8.8%, respectively, (2) estimated reasonably cap thicknesses larger than ∼300 µm, but underestimated thinner caps, and (3) quantified satisfactorily Young's modulus of hard medium (mean value of 109.7 ± 23.7 kPa instead of 145.4 ± 31.8 kPa), but overestimated the stiffness of soft inclusions (mean Young`s moduli of 31.4 ± 9.7 kPa instead of 17.6 ± 3.4 kPa). All together, these results demonstrate a promising benefit of the new iMOD-IVUS clinical imaging method for in vivo VP detection.

Functionalization of matrices by cyclically stretched osteoblasts through matrix targeting of VEGF.

Vascular endothelial growth factor A (VEGF) plays a central role in load-induced bone gain. We previously showed that increasing cyclic stretch frequency from 0.05 to 5 Hz induce parallel increased in entrapment of VEGF (mVEGF) into osteoblast secreted extracellular matrix. We ask in this study if mVEGF could be protective against apoptotic signals and biologically active in vitro on endothelial cell migration as well as in vivo on angiogenesis. We established that mechanically-induced VEGF entrapment using stretched silicone membrane was saturable after 3 exposures at high frequency stretches (5Hz). We found that mVEGF stimulates microvascular cells migration and enhanced angiogenesis more importantly than VEGF 165 controls suggesting the absence of potent anti-angiogenic factors in our functionalized matrices. Indeed we found that the anti-angiogenic factors, tissue inhibitor of metalloproteinase (TIMP2) and pigment epithelium-derived factor (PEDF) were specifically downregulated for 5 Hz stretch and that the release of these potent factors was increased for low frequency of stretch (0.05 Hz). This study qualifies high frequency cyclic stretch as an interesting approach for surfaces activation of deformable biomaterials.

An integrated formulation of anisotropic force-calcium relations driving spatio-temporal contractions of cardiac myocytes.

Isolated cardiac myocytes exhibit spontaneous patterns of rhythmic contraction, driven by intracellular calcium waves. In order to study the coupling between spatio-temporal calcium dynamics and cell contraction in large deformation regimes, a new strain-energy function, describing the influence of sarcomere length on the calcium-dependent generation of active intracellular stresses, is proposed. This strain-energy function includes anisotropic passive and active contributions that were first validated separately from experimental stress-strain curves and stress-sarcomere length curves, respectively. An extended validation of this formulation was then conducted by considering this strain-energy function as the core of an integrated mechano-chemical three-dimensional model of cardiac myocyte contraction, where autocatalytic intracellular calcium dynamics were described by a representative two-variable model able to generate realistic intracellular calcium waves similar to those observed experimentally. Finite-element simulations of the three-dimensional cell model, conducted for different intracellular locations of triggering calcium sparks, explained very satisfactorily, both qualitatively and quantitatively, the contraction patterns of cardiac myocytes observed by time-lapse videomicroscopy. This integrative approach of the mechano-chemical couplings driving cardiac myocyte contraction provides a comprehensive framework for analysing active stress regulation and associated mechano-transduction processes that contribute to the efficiency of cardiac cell contractility in both physiological and pathological contexts.

Nonlinear elastic properties of polyacrylamide gels: implications for quantification of cellular forces.

Because of their tunable mechanical properties, polyacrylamide gels (PAG) are frequently used for studying cell adhesion and migratory responses to extracellular substrate stiffness. Since these responses are known to heavily depend on the tensional balance between cell contractility and substrate mechanical resistance, a precise knowledge of PAG's mechanical properties becomes quite crucial. Using the micropipette aspiration technique, we first exhibited the nonlinear elastic behavior of PAG and then successfully modeled it by an original strain-energy function. This function depends on the Poisson's ratio and on two material parameters, which have been explicitly related to acrylamide and bis-acrylamide concentrations. Implications of these results have been highlighted with regard to traction force microscopy experiments where cellular force quantification is derived from displacements of beads embedded in PAG. We found that considering PAG as a linear elastic medium tends to significantly underestimate traction forces for substrate displacements larger than 2 mum. Interestingly, we also showed that in the range of cellular force amplitude and PAG stiffness currently used in cell traction force experiments, finite size effects become critical for PAG substrate thickness below 60 mum. Thus, our improved characterization of PAG nonlinear mechanical properties through a new constitutive law could have significant impact onto biological experimentations where such extracellular substrates experience large strains.

Vulnerable atherosclerotic plaque elasticity reconstruction based on a segmentation-driven optimization procedure using strain measurements: theoretical framework.

It is now recognized that prediction of the vulnerable coronary plaque rupture requires not only an accurate quantification of fibrous cap thickness and necrotic core morphology but also a precise knowledge of the mechanical properties of plaque components. Indeed, such knowledge would allow a precise evaluation of the peak cap-stress amplitude, which is known to be a good biomechanical predictor of plaque rupture. Several studies have been performed to reconstruct a Young's modulus map from strain elastograms. It seems that the main issue for improving such methods does not rely on the optimization algorithm itself, but rather on preconditioning requiring the best estimation of the plaque components' contours. The present theoretical study was therefore designed to develop: 1) a preconditioning model to extract the plaque morphology in order to initiate the optimization process, and 2) an approach combining a dynamic segmentation method with an optimization procedure to highlight the modulogram of the atherosclerotic plaque. This methodology, based on the continuum mechanics theory prescribing the strain field, was successfully applied to seven intravascular ultrasound coronary lesion morphologies. The reconstructed cap thickness, necrotic core area, calcium area, and the Young's moduli of the calcium, necrotic core, and fibrosis were obtained with mean relative errors of 12%, 4% and 1%, 43%, 32%, and 2%, respectively.

Necrotic core thickness and positive arterial remodeling index: emergent biomechanical factors for evaluating the risk of plaque rupture.

Fibrous cap thickness is often considered as diagnostic of the degree of plaque instability. Necrotic core area (Core(area)) and the arterial remodeling index (Remod(index)), on the other hand, are difficult to use as clinical morphological indexes: literature data show a wide dispersion of Core(area) thresholds above which plaque becomes unstable. Although histopathology shows a strong correlation between Core(area) and Remod(index), it remains unclear how these interact and affect peak cap stress (Cap(stress)), a known predictor of rupture. The aim of this study was to investigate the change in plaque vulnerability as a function of necrotic core size and plaque morphology. Cap(stress) value was calculated on 5,500 idealized atherosclerotic vessel models that had the original feature of mimicking the positive arterial remodeling process described by Glagov. Twenty-four nonruptured plaques acquired by intravascular ultrasound on patients were used to test the performance of the associated idealized morphological models. Taking advantage of the extensive simulations, we investigated the effects of anatomical plaque features on Cap(stress). It was found that: 1) at the early stages of positive remodeling, lesions were more prone to rupture, which could explain the progression and growth of clinically silent plaques and 2) in addition to cap thickness, necrotic core thickness, rather than area, was critical in determining plaque stability. This study demonstrates that plaque instability is to be viewed not as a consequence of fibrous cap thickness alone but rather as a combination of cap thickness, necrotic core thickness, and the arterial remodeling index.

A computational model of cell migration coupling the growth of focal adhesions with oscillatory cell protrusions.

Cell migration is a highly integrated process where actin turnover, actomyosin contractility, and adhesion dynamics are all closely linked. In this paper, we propose a computational model investigating the coupling of these fundamental processes within the context of spontaneous (i.e. unstimulated) cell migration. In the unstimulated cell, membrane oscillations originating from the interaction between passive hydrostatic pressure and contractility are sufficient to lead to the formation of adhesion spots. Cell contractility then leads to the maturation of these adhesion spots into focal adhesions. Due to active actin polymerization, which reinforces protrusion at the leading edge, the traction force required for cell translocation can be generated. Computational simulations first show that the model hypotheses allow one to reproduce the main features of fibroblast cell migration and established results on the biphasic aspect of the cell speed as a function of adhesion strength. The model also demonstrates that certain temporal parameters, such as the adhesion proteins recycling time and adhesion lifetimes, influence cell motion patterns, particularly cell speed and persistence of the direction of migration. This study provides some elements, which allow a better understanding of spontaneous cell migration and enables a first glance at how an individual cell would potentially react once exposed to a stimulus.

The motility of normal and cancer cells in response to the combined influence of the substrate rigidity and anisotropic microstructure.

Cell adhesion and migration are strongly influenced by extracellular matrix (ECM) architecture and rigidity, but little is known about the concomitant influence of such environmental signals to cell responses, especially when considering cells of similar origin and morphology, but exhibiting a normal or cancerous phenotype. Using micropatterned polydimethylsiloxane substrates (PDMS) with tunable stiffness (500 kPa, 750 kPa, 2000 kPa) and topography (lines, pillars or unpatterned), we systematically analyse the differential response of normal (3T3) and cancer (SaI/N) fibroblastic cells. Our results demonstrate that both cells exhibit differential morphology and motility responses to changes in substrate rigidity and microtopography. 3T3 polarisation and spreading are influenced by substrate microtopography and rigidity. The cells exhibit a persistent type of migration, which depends on the substrate anisotropy. In contrast, the dynamic of SaI/N spreading is strongly modified by the substrate topography but not by substrate rigidity. SaI/N morphology and migration seem to escape from extracellular cues: the cells exhibit uncorrelated migration trajectories and a large dispersion of their migration speed, which increases with substrate rigidity.

Estimation of cell Young's modulus of adherent cells probed by optical and magnetic tweezers: influence of cell thickness and bead immersion.

A precise characterization of cell elastic properties is crucial for understanding the mechanisms by which cells sense mechanical stimuli and how these factors alter cellular functions. Optical and magnetic tweezers are micromanipulation techniques which are widely used for quantifying the stiffness of adherent cells from their response to an external force applied on a bead partially embedded within the cell cortex. However, the relationships between imposed external force and resulting bead translation or rotation obtained from these experimental techniques only characterize the apparent cell stiffness. Indeed, the value of the estimated apparent cell stiffness integrates the effect of different geometrical parameters, the most important being the bead embedding angle 2gamma, bead radius R, and cell height h. In this paper, a three-dimensional finite element analysis was used to compute the cell mechanical response to applied force in tweezer experiments and to explicit the correcting functions which have to be used in order to infer the intrinsic cell Young's modulus from the apparent elasticity modulus. Our analysis, performed for an extensive set of values of gamma, h, and R, shows that the most relevant parameters for computing the correcting functions are the embedding half angle gamma and the ratio h(u)/2R, where h(u) is the under bead cell thickness. This paper provides original analytical expressions of these correcting functions as well as the critical values of the cell thickness below which corrections of the apparent modulus are necessary to get an accurate value of cell Young's modulus. Moreover, considering these results and taking benefit of previous results obtained on the estimation of cell Young's modulus of adherent cells probed by magnetic twisting cytometry (MTC) (Ohayon, J., and Tracqui, P., 2005, Ann. Biomed. Eng., 33, pp. 131-141), we were able to clarify and to solve the still unexplained discrepancies reported between estimations of elasticity modulus performed on the same cell type and probed with MTC and optical tweezers (OT). More generally, this study may strengthen the applicability of optical and magnetic tweezers techniques by insuring a more precise estimation of the intrinsic cell Young's modulus (CYM).

Influence of residual stress/strain on the biomechanical stability of vulnerable coronary plaques: potential impact for evaluating the risk of plaque rupture.

In a vulnerable plaque (VP), rupture often occurs at a site of high stress within the cap. It is also known that vessels do not become free of stress when all external loads are removed. Previous studies have shown that such residual stress/strain (RS/S) tends to make the stress distribution more uniform throughout the media of a normal artery. However, the influence of RS/S on the wall stress distribution in pathological coronaries remains unclear. The aim of this study was to investigate the effects of RS/S on the biomechanical stability of VPs. RS/S patterns were studied ex vivo in six human vulnerable coronary plaque samples. Because the existence of RS/S can only be assessed by releasing it, the opening angle technique was the experimental approach used to study the geometrical opening configurations of the diseased arteries, producing an arterial wall in a near-zero stress state. Reciprocally, these opening geometries were used in finite element simulations to reconstruct the RS/S distributions in closed arteries. It was found that the RS/S 1) is not negligible, 2) dramatically affects the physiological peak stress amplitude in the thin fibrous cap, 3) spotlights some new high stress areas, and 4) could be a landmark of the lipid core's developmental process within a VP. This study demonstrates that plaque rupture is not to be viewed as a consequence of intravascular pressure alone, but rather of a subtle combination of external loading and intraplaque RS/S.

An extended relationship for the characterization of Young's modulus and Poisson's ratio of tunable polyacrylamide gels.

Substrates with tunable mechanical properties are crucial for the study of cellular processes, and polyacrylamide gels (PAGs) are frequently used in this context. Several experimental techniques have been proposed to obtain the mechanical properties of PAGs. However, the range of the considered Poisson's ratio values remains quite large and no attempt has been made to propose an analytical relationship allowing the estimation of PAG Young's modulus when both bis-acrylamide and acrylamide concentrations are known. In order to complete the actual knowledge on the mechanical properties of PAGs, we took benefit of our original method based on the micropipette aspiration technique (Boudou et al., J. Biomech. 2006) for characterizing gels made with concentrations in the range 0.02% < or =[Bis]< or =0.20% and 3% < or =[Acry]< or =10%. We found that the PAGs Young's modulus varies nonlinearly with the acrylamide amount. Moreover, our study validates the quasi-incompressibility hypothesis usually made in studies using PAGs (mean Poisson's ratio of 0.480+/-0.012). More generally, and in agreement with data published by other groups, we propose an original nonlinear mathematical relationship allowing the computation of Young's modulus of PAG for any given acrylamide and bis-acrylamide amounts taken in the range of values we considered.

QxDB: a generic database to support mathematical modelling in biology.

QxDB (quantitative x-modelling database) is a web-based generic database package designed especially to house quantitative and structural information. Its development was motivated by the need for centralized access to such results for development of mathematical models, but its usefulness extends to the general research community of both modellers and experimentalists. Written in PHP (Hyper Preprocessor) and MYSQL, the database is easily adapted to new fields of research and ported to Apache-based web servers. Unlike most existing databases, experimental and observational results curated in QxDB are supplemented by comments from the experts who contribute input to the database, giving their evaluations of experimental techniques, breadth of validity of results, experimental conditions, and the like, thus providing the visitor with a basis for gauging the quality (or appropriateness) of each item for his/her needs. QxDB can be easily customized by adapting the contents of the database table containing the descriptors that characterize each data record according to an informal ontology of the research domain. We will illustrate this adaptability of QxDB by presenting two examples, the first dealing with modelling in oncology and the second with mechanical properties of cells and tissues.

Spatiotemporal dynamics of actin-rich adhesion microdomains: influence of substrate flexibility.

In this study we analyse the formation and dynamics of specific actin-rich structures called podosomes. Podosomes are very dynamic punctual adhesion sites tightly linked to the actin cytoskeleton. Mechanical properties of substrates are emerging as important physical modulators of anchorage-dependent processes involved in the cellular response. We investigate the influence of substrate flexibility on the dynamic properties of podosomes. We used mouse NIH-3T3 fibroblasts, transfected with GFP-actin and cultured on polyacrylamide collagen-coated substrates of varying stiffness. Static and dynamic features of cell morphologies associated with an optical flow analysis of the dynamics of podosomes revealed that: (1) they have constant structural properties, i.e. their shape factor and width do not change with the substrate flexibility; (2) the lifespan of podosomes and mean minimum distance between them depend on the substrate flexibility; (3) there is a variation in the displacement speed of the rosette of podosomes. Moreover, the rosettes sometimes appear as periodically emergent F-actin structures, which suggests that a two-level self-organisation process may drive first, the formation of clusters of podosomes and second, the organisation of these clusters into oscillating rings. Such dynamic features give new perspectives regarding the potential function of podosomes as mechanosensory structures.