X ENS Physique B PC 2021Sujet, corrigé et rapport du jury

Mécanique de la cellule

Pas encore noté

Faisable en Sup

Faisable en sup ?

Faisable en partie en PCSI : 28 questions sur 48, parties III et IV en entier. Autres questions faisables en PCSI

Autres premières années : MP2I 27 · BCPST1 21 questions.

Détail par partie

Partie I : Questions préliminaires et ordres de grandeur0/32e année
Partie II : Suivi de particule unique et viscosité du cytoplasme0/172e année
Partie III : Expériences d'indentation et d'élasticité de la cellule14/14Faisable
Partie IV : Modèles visco-élastiques du cytoplasme cellulaire14/14Faisable

Explorer les annales

BanqueX-ENS Année2021 FilièrePC MatièrePhysique

Annales X-ENS PC Physique X-ENS PC 2021 PC Physique Physique 2021

Téléchargements

Télécharger le sujet (PDF)

Corrigé détaillé en ligne

Corrigé sur Discord

Rapport du jury Lire en ligne LaTeX Word

L'épreuve en chiffres

Moyenne 10,33 / 20 · écart-type 2,96 · 1 008 présents · où vous situez-vous ?

Moyenne: 10,33/ 20
Écart-type: 2,96
Présents: 1 008

Deux tiers des copies environ (moyenne ± écart-type)

Où je me situe ?

Votre note sur 20 à ce sujet, en conditions de concours.

Source : document officiel du concours. Notes publiées par le concours (après harmonisation le cas échéant). Courbe : estimation par une loi normale.

Corrigé détaillé

48 questions

Indice, stratégie puis corrigé complet pour chaque question, avec les remarques du jury.

48 questions corrigées une à une : un indice, puis la stratégie, puis le corrigé complet avec le résultat encadré et les remarques du jury. Chaque question est vérifiée avant publication.

Voir les 40 autres questions

Ouvrir le corrigé détaillé Comment sont produits nos corrigés

Description

Ce sujet de l'épreuve de Physique B (XEULC) aborde la mécanique de la matière vivante à travers l'étude approfondie des propriétés rhéologiques du cytoplasme et de la cellule eucaryote.

Ces sujets peuvent vous intéresser

Lecture du sujet en ligne

L'énoncé complet, avec les formules et les figures, sans ouvrir le PDF.

Sujet en PDF

ECOLE POLYTECHNIQUE - ESPCI ECOLES NORMALES SUPERIEURES

CONCOURS D'ADMI SSION 2021

MERCREDI 14 AVRIL 2021
08h00-12h00
FILIERE PC - Epreuve

n^\circ5

PHYSI QUE B (XEULC)

Durée : 4 heures
L'utilisation des calculatrices n'est pas
autorisée pour cette épreuve

Mécanique de la cellule

Le but de ce problème est d'étudier quelques propriétés mécaniques des cellules biologiques. D'un point de vue physique, une cellule est un objet 'mou' qui consomme de l'énergie sous forme d'ATP. Il s'agit donc d'un système hors d'équilibre thermodynamique. Nous nous intéresserons ici essentiellement à la caractérisation de la rhéologie cellulaire, c'est-à-dire aux propriétés mécaniques de la cellule. Le formulaire donne les valeurs des paramètres utiles pour les estimations d'ordre de grandeur et les applications numériques ainsi qu'un bref formulaire d'analyse vectorielle et une figure rappelant les caractéristiques d'une ellipse lors d'une étude harmonique.

Les deux dernières parties du problème (III et IV) sont indépendantes des deux premières (I et II). Des résultats intermédiaires au sein des différents paragraphes permettent de progresser dans le sujet.

Notations, formulaire et données numériques

Vitesse caractéristique de l'écoulement dans le cytoplasme : $u = 0, 1μ m.s^(- 1)$
Taille caractéristique d'une cellule : $L = 10μ m$
Viscosité dynamique du cytoplasme : $η_c = 1$ Pa.s $^(- 1)$
Masse volumique du cytoplasme : $ρ_c = 1, 00 \times 10^3 kg \cdot m^(- 3)$
Viscosité cinématique du cytoplasme : $ν_c = (η_c)/(ρ_c) = 1 \times 10^(- 3) m^2 \cdot s^(- 1)$
Rayon d'une bille de latex : $R_b = 1μ m$
Température : $T = 310 K$
Constante de Boltzmann : $k_B = 1, 4 \times 10^(- 23) J \cdot K^(- 1)$
Expression de l'opérateur laplacien scalaire $Δφ = div(grad^(\to-)φ)$ , où $φ$ est un champ scalaire

Figure 1 - Caractéristiques d'une ellipse pour une étude harmonique.

I Questions préliminaires et ordres de grandeur

La figure 2 schématise une cellule eucaryote

^1

. Le problème étudiera les propriétés mécaniques de la cellule entière ou de son cytoplasme. La taille caractéristique de la cellule considérée se situe entre 10 et

100μ m

. L'intérieur de la cellule, appelé cytoplasme, est un milieu riche en protéines, compartiments membranaires et filaments du cytosquelette. Du point de vue de la physique de la matière molle, le cytoplasme est un fluide complexe dont les propriétés mécaniques font l'objet du sujet. Les vitesses de déplacements d'objets à l'intérieur du cytoplasme, par exemple les vésicules de transport intracellulaire ou les filaments du cytosquelette qui polymérisent ou dépolymérisent, sont de l'ordre de

0, 1μ m \cdot s^(- 1)

Dans cette partie préliminaire, nous supposons que le cytoplasme est un fluide newtonien de masse volumique

ρ_c

. Nous utiliserons des arguments d'ordre de grandeur, ainsi que les grandeurs sans dimension définies en mécanique des fluides pour simplifier l'étude du cytoplasme.

Figure 2 - Schéma d'une cellule eucaryote. Les différents compartiments cellulaires évoqués dans l'énoncé sont représentés, ainsi que les échelles de longueurs utilisées.

On désigne par $v⃗(r⃗, t)$ le champ de vitesse du fluide intracellulaire en un point $M(r⃗)$ du cytoplasme et à l'instant $t$ , dans un repère $R = (Oxyz)$ associé au référentiel de la cellule, supposé galiléen.

Nous rappelons l'équation de Navier-Stokes pour un écoulement incompressible :

ρ_c d/(dt)v⃗ = f⃗ - grad^(\to-)p + η_c Δv⃗,

où

f⃗

désigne la résultante des forces volumiques extérieures exercée sur une particule de fluide du cytoplasme,

p

est la pression,

η_c

est la viscosité dynamique du cytoplasme. L'écoulement dans le cytoplasme est considéré comme incompressible. L'opérateur

Δ

représente l'opérateur laplacien. La dérivée particulaire

d/(dt)

a pour expression :

d/(dt) = \partial/(\partialt) + v⃗ \cdot grad^(\to-)

Afin de simplifier l'étude du cytoplasme, il est nécessaire de comparer les effets de différents phénomènes entre eux. Pour cela, nous introduisons des grandeurs sans dimension et différents temps caractéristiques associés à certains phénomènes pour identifier ceux qui sont déterminants dans le comportement mécanique du cytoplasme.

On définit le nombre de Reynolds, noté Re, comme le rapport sans dimension entre le terme convectif de la dérivée particulaire et les forces volumiques de viscosité intervenant dans l'équation de Navier-Stokes (ou, de façon équivalente, comme le rapport entre le temps caractéristique de diffusion de la quantité de mouvement et le temps caractéristique de convection).

Donner l'expression de Re pour l'écoulement dans le cytoplasme cellulaire en fonction de la dimension caractéristique de l'extension de l'écoulement $L$ , de la viscosité cinématique $ν_c$ et de la vitesse caractéristique $u$ . Donner un ordre de grandeur de Re et conclure sur le rôle du terme convectif dans les écoulements du cytoplasme à l'intérieur de la cellule.

Nous supposons que les forces volumiques extérieures et de pression exercées sur une particule de fluide dans le cytoplasme sont négligeables devant les forces volumiques de viscosité.

Pour estimer la contribution du terme inertiel (terme en

\partialv⃗/\partialt

dans l'expression de la dérivée particulaire), on définit un temps caractéristique

τ_(inertie)

tel que

|ρ_c(\partialv⃗)/(\partialt)| \sim ρ_c u/(τ_(inertie)).

On fait l'hypothèse que les forces volumiques de viscosité et le terme inertiel dans l'équation de Navier-Stokes sont du même ordre de grandeur. Exprimer $τ_(inertie)$ en fonction de $L, ρ_c$ et $η_c$ puis calculer numériquement $τ_(inertie)$ . En comparant $τ_(inertie)$ avec les temps caractéristiques (supérieurs à la milliseconde) mis en jeu au cours des expériences de rhéologie cellulaire que nous évoquerons dans la suite, justifier le fait que l'on conserve les forces volumiques de viscosité et le terme inertiel dans l'équation de Navier-Stokes.
Déduire des approximations et hypothèses précédentes l'équation qui régit la dynamique de l'écoulement cytoplasmique. Rappeler le nom de ce type d'équation et donner une interprétation du coefficient de viscosité cinématique $ν_c$ relativement à cette équation.

II Suivi de particule unique et viscosité du cytoplasme

On étudie dans cette partie les phénomènes de diffusion dans le cytoplasme cellulaire, considéré comme un fluide newtonien. Les paragraphes préliminaires (II.A, II.B et II.C) abordent quelques modèles simplifiés. Ils nous permettront, à partir d'une étude microscopique de la diffusion de particules d'un gaz, et celle d'une particule en mouvement brownien dans un fluide, de relier le coefficient de diffusion

D

d'une particule brownienne à la viscosité dynamique

η

d'un fluide. Cette relation, appelée relation d'Einstein, s'applique à un fluide newtonien. Elle est à la base d'une méthode de rhéologie passive permettant d'accéder à la viscosité

η_c

du cytoplasme à partir du suivi d'une particule unique (II.D).

II.A Approche microscopique de la diffusion dans un gaz

On considère un gaz monoatomique, dont les particules ont une masse

m

, un rayon

r_p

et une vitesse moyenne d'agitation thermique

v¯

(cette vitesse est la vitesse quadratique moyenne définie comme la racine carrée de la moyenne de la vitesse au carré soit

v¯ = \sqrt(⟨v^2⟩)

). Le nombre total de particules est noté

N

, le volume du gaz est noté

V

et sa température thermodynamique est notée

T

. Le gaz est considéré à l'équilibre thermique. On notera

k_B

la constante de Boltzmann,

n = N/V

le nombre de particules du gaz par unité de volume et

ρ = Nm/V = nm

sa masse volumique.

On modélise les particules du gaz comme des billes de rayon

r_p

, qui ne peuvent s'interpénétrer, et qui peuvent subir des collisions entre elles (ce modèle s'appelle le modèle des sphères dures). On note

τ_(col)

le temps de collision défini comme le temps moyen entre deux collisions successives. La distance moyenne parcourue par une particule entre deux collisions successives subies par la même particule est appelée libre parcours moyen, noté

ℓ

. Enfin, la section efficace de choc d'une particule,

σ_(eff) = 4πr_p^2

, est la surface plane perpendiculaire à la vitesse d'une particule définissant son extension spatiale.

Nous étudions dans ce paragraphe le transport de particules du gaz par diffusion au sein de ce gaz.
4. Relier le libre parcours moyen

ℓ

au temps de collision

τ_(col)

et à la vitesse

v¯

. Quel est l'ordre de grandeur du nombre de particules se trouvant en moyenne dans le volume

σ_(eff)ℓ

? En déduire que

ℓ = c/(nσ_(eff))

, où

c

est une constante numérique. On admettra que

c = 1/(\sqrt2)

. On pourra effectuer un schéma pour illustrer le raisonnement.

On s'intéresse ici à la diffusion des particules de gaz se déplaçant sous l'effet d'un gradient de nombre de particules par unité de volume $n(x)$ unidimensionnel. La vitesse des particules est supposée uniforme (indépendante de l'abscisse $x$ ) et égale à la vitesse quadratique moyenne du gaz, $v¯$ . On considère un volume $V = 2Sℓ$ de gaz situé entre les abscisses $x - ℓ$ et $x + ℓ$ (voir figure 3 ). Le nombre volumique de particules de gaz à l'abscisse $x - ℓ$ (resp. $x + ℓ$ ) est notée $n(x - ℓ)$ (resp. $n(x + ℓ)$ ). On note $Φ_(N, e)$ le flux de particules de gaz à l'abscisse $x - ℓ$ , par unité de surface $S$ , entrant dans le volume $V$ et $Φ_(N, s)$ le flux de particules de gaz à l'abscisse $x + ℓ$ , par unité de surface $S$ , sortant du volume $V$ .

On suppose que les particules se déplacent dans la direction

+ x

et que le nombre de particules par unité de volume

n(x)

varie lentement sur des distances de l'ordre de

ℓ

. S'il entre

dN_e

particules de gaz dans le volume

V

à l'abscisse

x - ℓ

entre les instants

t

t + dt

, alors

Φ_(N, e) = 1/S(dN_e)/(dt)

. De la même façon, s'il sort

dN_s

particules de gaz du volume

V

à l'abscisse

x + ℓ

entre les instants

t

t + dt

, alors

Φ_(N, s) = 1/S(dN_s)/(dt)

Figure 3 - Modèle à une dimension pour l'approche microscopique de la diffusion de particules dans un gaz.

Montrer qu'entre les dates $t$ et $t + dt, dN_e$ particules de gaz entrent dans le volume $V$ et $dN_s$ particules de gaz en sortent, avec $dN_e = Sv¯n(x - ℓ)dt$ et $dN_s = Sv¯n(x + ℓ)dt$ . À l'aide d'un bilan portant sur le nombre de particules dans le volume $V$ , en déduire que le flux de particules par unité de surface dans la direction $+ x$ défini par $Φ_(N, + x) = Φ_(N, e) - Φ_(N, s)$ s'écrit $Φ_(N, + x) = - 2v¯ℓ(dn)/(dx)$ .
En tenant compte des déplacements des particules dans toutes les directions de l'espace, on peut montrer que le flux de particules par unité de surface dans $V$ est $Φ_N = 1/6Φ_(N, x) = - 1/3v¯ℓ(dn)/(dx)$ . En identifiant la relation précédente à la loi de FiCK reliant le flux de particules au gradient de nombre de particules $n(x)$ , montrer que le coefficient de diffusion $D$ des particules de gaz s'écrit $D = 1/3v¯ℓ$ .
À l'aide du théorème d'équipartition de l'énergie : $1/2mv¯^2 = 3/2k_B T$ , montrer que le coefficient de diffusion $D$ s'écrit :

D = (k_B T)/mτ_(col).

II.B Mouvement brownien et modèle de marche au hasard

Nous nous intéressons dans la suite au mouvement brownien d'une particule dans un fluide sous l'effet de l'agitation thermique. Nous allons établir, pour une telle particule, la relation

D = A(ℓ^2)/(τ_(col))

, où

A

est un facteur numérique que nous déterminerons.

Pour décrire le mouvement d'une telle particule soumise à l'agitation thermique, on considère une marche aléatoire dans un espace à trois dimensions de

N

pas partant du point

A = P_0

fixe et arrivant au point

B = P_N

, dont la position dépend de la trajectoire suivie (paramétrée par les

N + 1

points

(P_0, P_1, \dots, P_i, \dots, P_N)

, avec

i = 0, 1, \dots, N)

. On note

d_i\to

le vecteur tel que

d_i\to = P_i P_(i + 1)^(\to-)

. La suite des vecteurs

d_i\to

, où l'indice

i

prend ici des valeurs entre 0 et

N - 1

, représente une trajectoire (une réalisation) suivie par la particule entre les points

A

B

fixés (voir figure 4). Les points

P_i

, appelés points de collision, sont les points où la particule change de direction. La durée de parcours entre chaque point de collision est en moyenne égale au temps de collision

τ_(col)

Le temps moyen mis par la particule pour aller de

A

B

est noté

t

Figure 4 - Modèle de marche aléatoire à

N

pas dans un espace à trois dimensions.

Nous adoptons ci-dessous quelques hypothèses et définissons certaines notations concernant la marche au hasard présentée ci-dessus et schématisée sur la figure 4.

⋆

La marche au hasard est isotrope : toutes les orientations des vecteurs

d_i\to

sont équiprobables.

La marche au hasard est homogène : chaque point de collision obéit à la même loi de probabilité.
On désigne par $⟨f⟩$ la moyenne d'une grandeur $f(t)$ sur tous les chemins possibles à $N$ pas, entre A et B. On admet les propriétés suivantes : $d/(dt)⟨f⟩ = ⟨(df)/(dt)⟩$ et $⟨fg⟩ = ⟨f⟩⟨g⟩$ si $f$ et $g$ sont deux grandeurs indépendantes.
On note $f¯$ , la moyenne temporelle de $f$ sur un chemin donné entre A et B . L'hypothèse dite d'ergodicité stipule $⟨f⟩ = f¯$ .
Nous n'étudions pas le détail des collisions et nous nous limitons aux caractéristiques de la marche aléatoire évoquées ci-dessus.

On désigne par

δ⃗ = \sum_(i = 0)^(N - 1)d_i\to = AB^(\to-)

un chemin possible pour aller de A à B en effectuant

N

pas. La marche aléatoire étant homogène et isotrope,

⟨δ⃗⟩ = 0\to

, ce qui signifie que B est situé sur une sphère
de centre A et de rayon moyen au carré

⟨‖AB^(\to-)‖^2⟩

.
8. Montrer que

⟨AP_(i + 1)^2⟩ = ⟨AP_i^2⟩ + ⟨d_i^2⟩ + ⟨2d⃗_i \cdot AP^(\to-)_i⟩

.
9. On note

θ_i

l'angle entre les vecteurs

AP^(\to-)_i

d⃗_i

(voir figure 4). Le libre parcours moyen

ℓ

est ici la valeur moyenne qu'aurait la composante du vecteur

d_i\to

dans l'une des trois directions de l'espace, repérées par les vecteurs

e_x^(\to-), e_y\to

e_z\to

, formant une base orthonormée. On a donc

⟨d_(i, x)^2⟩ = ⟨d_(i, y)^2⟩ = ⟨d_(i, z)^2⟩ = ℓ^2

. En considérant que les variables

‖d_i\to‖, ‖AP_i^(\to-)‖

θ_i

sont des grandeurs indépendantes, montrer que

⟨AP_(i + 1)^2⟩ = ⟨AP_i^2⟩ + 3ℓ^2

.
10. En déduire que

⟨‖δ⃗‖^2⟩ = 3Nℓ^2

.
11. Indiquer la relation entre le temps moyen

t

pour aller de A à B , le temps de collision

τ_(col)

et le nombre de pas de la marche aléatoire

N

. En déduire alors l'expression de

⟨‖δ⃗‖^2⟩

en fonction de

t, ℓ

τ_(col)

.
12. Le coefficient de diffusion de particules à trois dimensions étant défini par

‖δ⃗‖^2^- = 6Dt

, préciser l'expression de

D

en fonction du libre parcours moyen

ℓ

et du temps de collision

τ_(col)

. Comparer cette expression à celle de

D

, obtenue par l'approche microscopique de la diffusion de particules à la question 6 .

II.C Particule brownienne mésoscopique et modèle de Langevin

Nous nous intéressons maintenant à une particule brownienne de grandes dimensions devant les dimensions atomiques et se déplaçant dans un fluide à la température

T

. Nous utilisons le modèle de Langevin pour établir la relation d'Einstein, reliant le coefficient de diffusion

D

de cette particule brownienne et la viscosité dynamique du fluide. Notons que le coefficient de diffusion que nous obtiendrons dans cette situation ne possède pas le même statut que celui du paragraphe II.A, qui décrivait l'auto-diffusion dans un gaz.

On considère des particules sphériques de rayon

r_p

, de masse

m

, se déplaçant à la vitesse

V⃗

dans un liquide de viscosité dynamique

η

. La position d'une particule au point M est repérée par le vecteur

OM^(\to-)(t) = (x(t), y(t), z(t))

dans un système de coordonnées cartésiennes de repère (

O, e_x^(\to-), e_y\to, e_z\to

). Les particules sont soumises à :

une force de frottement fluide due au liquide, dont l'expression est donnée par la force de Stokes $F⃗_η = - 6πηr_p V⃗$ , avec $V⃗$ la vitesse d'une particule;
$⋆$ une force aléatoire $f⃗(t) = (f_x, f_y, f_z)$ telle que $⟨f_x⟩ = ⟨f_y⟩ = ⟨f_z⟩ = 0$ (où $⟨\dots⟩$ désigne la moyenne sur un grand nombre de particules). Cette moyenne est l'analogue de la moyenne sur les chemins possibles, définie précédemment dans le cadre de la marche aléatoire, et possède les mêmes propriétés.

On suppose que l'équilibre thermodynamique est atteint.
13. En utilisant le théorème d'équipartition de l'énergie, déterminer l'expression de

⟨x˙^2⟩

, où la notation

x˙

désigne la dérivée de

x(t)

par rapport au temps, composante de la vitesse dans la direction

e_x^(\to-)

.
14. On suppose que les variables

x

f_x

sont indépendantes. À partir du principe fondamental de la dynamique appliqué à une particule, établir l'équation différentielle vérifiée par

⟨xx˙⟩

.
15. Résoudre cette équation différentielle et montrer qu'en régime permanent, on obtient

⟨xx˙⟩ = k_B T/6πηr_p

. On introduira une durée caractéristique

τ

et on donnera un critère d'atteinte du régime permanent.
16. En supposant que

⟨x^2⟩(t = 0) = 0

et des durées longues devant

τ

, montrer que

⟨x^2⟩(t) = 2Dt

et établir la relation d'Einstein

D = (k_B T)/(6πηr_p).

II.D Rhéologie passive du cytoplasme cellulaire par suivi de particule unique

Une technique expérimentale de rhéologie passive pour mesurer la viscosité

η_c

du cytoplasme cellulaire consiste à suivre au microscope le déplacement d'une bille de taille micrométrique (rayon

R_b

), qui a été introduite dans le cytoplasme, au cours du temps (figure 5). La bille est visualisée à intervalles de temps réguliers par vidéomicroscopie. On note

Δt_(min)

l'intervalle de temps entre deux images (ou temps d'acquisition des images) et

N_(tot)

le nombre total d'images. Un algorithme de suivi de particule unique permet ensuite d'obtenir les coordonnées

x_b(t)

y_b(t)

du centre de masse de la bille et sa trajectoire en deux dimensions, à une résolution inférieure au pixel de la caméra utilisée. À partir de

x_b(t)

y_b(t)

, on calcule le déplacement quadratique moyen

δ^2(t) = ⟨x_b(t)^2 + y_b(t)^2⟩

, à la date

t = t_i = iΔt_(min)

où

i = 1\dotsN_(tot)

Figure 5 - Schéma d'une cellule eucaryote en présence d'une bille. Une bille de diamètre

2R_b

est placée à l'intérieur de la cellule. Elle permet d'étudier les actions mécaniques exercées dans le cytoplasme en jouant le rôle de particule-sonde de l'écoulement dans le cytoplasme.

En supposant que le mouvement de la bille dans le cytoplasme est purement diffusif (on notera $D$ son coefficient de diffusion), représenter graphiquement l'évolution temporelle de $δ^2(t)$ (de façon qualitative) et expliquer comment la mesure de $δ^2(t)$ permet d'accéder à la viscosité du cytoplasme $η_c$ .
On visualise le déplacement d'une bille de rayon $R_b = 1μ m$ à 6 images par seconde. Le déplacement quadratique moyen entre deux images est mesuré et on obtient 0,2 pixel. Sachant que les images obtenues sont calibrées à $0, 05μ m$ par pixel, en déduire une estimation numérique du coefficient de diffusion de la bille, puis de la viscosité du cytoplasme.

En pratique, le mouvement de la bille aux temps longs est rarement gouverné par le seul phénomène de diffusion. Certains constituants du cytoplasme peuvent avoir une influence sur le mouvement de la bille et induire une déviation par rapport au mouvement diffusif usuel.

Des assemblages moléculaires, appelés moteurs moléculaires, peuvent entraîner la bille à une vitesse moyenne de dérive $v¯_(mot)$ . Ce phénomène vient s'ajouter au phénomène de diffusion.
Les filaments du cytosquelette forment un réseau complexe dont la maille caractéristique est de l'ordre de 50 à 100 nm (Fig. 5). Le mouvement de la bille peut être restreint dans une zone de l'espace de rayon caractéristique $R_(cage)$ , à l'intérieur du cytoplasme.

Illustrer à l'aide d'un schéma la trajectoire de la bille en présence de moteurs moléculaires. On supposera que le mouvement de dérive à la vitesse $v¯_(mot)$ dû aux moteurs moléculaires se superpose au mouvement dû à la diffusion, et conduit à un comportement diffusif avec dérive. Établir l'équation donnant $δ^2(t)$ en fonction de $v¯_(mot), D$ et de $t$ . Vérifier que le mouvement est bien diffusif sans dérive aux temps courts. Représenter graphiquement, et de façon qualitative, l'évolution temporelle de $δ^2(t)$ en représentation $log$ -log.
Représenter à l'aide d'un schéma la trajectoire de la bille en situation de diffusion confinée dans une zone restreinte de l'espace. Montrer que l'expression $δ^2(t) = αR_(cage)^2[1 - exp(- 4Dt/R_(cage)^2)]$ , avec $α$ un facteur numérique, peut modéliser un comportement où le mouvement de la bille est restreint dans une cage de rayon caractéristique $R_(cage)$ . Vérifier que le mouvement est bien diffusif sans confinement aux temps courts. Représenter graphiquement, et de façon qualitative, l'évolution temporelle de $δ^2(t)$ en représentation $log$ -log.

III Expériences d'indentation et d'élasticité de la cellule

Dans cette partie nous allons nous intéresser aux déformations élastiques de la cellule. Pour cela, nous définissons ici les notions de déformation et de contrainte, notions centrales de la théorie de l'élasticité des matériaux. La déformation

ε(t)

du matériau est une grandeur sans dimension qui dépend de la géométrie du système physique modélisé. Différentes situations sont représentées sur la figure 6. En général, la déformation

ε

est égale à

ε = x/L

où

x

est la distance sur laquelle le matériau est déformé relativement à une longueur caractéristique

L

de l'échantillon. La contrainte

σ(t)

exercée sur le matériau est définie comme

σ = F/S

où

F

est la force qui provoque la déformation appliquée sur une surface

S

et homogène à une pression. La déformation et la contrainte sont des grandeurs algébriques.

Figure 6 - (a) : Déformation d'un matériau par étirement. (b) : Déformation d'un matériau par cisaillement. (a) Dans le cas d'un étirement uniaxial

ΔL = L(t) - L_0

sous l'effet d'une force verticale, la déformation est

ε = ΔL/L_0

et la contrainte est

σ = F/S

où la force d'étirement

F⃗

est normale à la surface

S

. (b) Dans le cas d'un cisaillement horizontal sous l'effet d'une force horizontale, le plan

z = 0

est fixe et la déformation est

ε = Δx/L

où

Δx

est le déplacement horizontal du matériau à une hauteur

z = L

dans le matériau. La contrainte est

σ = F/S

où la force de cisaillement

F⃗

est tangentielle à la surface

S

. On note

θ

l'angle d'inclinaison par rapport à la verticale.

III.A Déformations d'un solide élastique

Pour un matériau homogène, linéaire et purement élastique, la déformation est proportionnelle à la contrainte appliquée. On parle alors de solide de Hooke. La constante de proportionnalité entre la contrainte et la déformation est analogue à la constante de raideur

k

d'un ressort obéissant à la loi de Hooke

F = kΔx

où

Δx

est l'allongement du ressort et

F

la norme de la force de rappel du ressort. Nous étudions dans la suite deux cas particuliers simples de déformations d'un solide élastique pour comprendre quelques éléments de la théorie de l'élasticité.

III.A.a Étirement longitudinal

On considère un cylindre d'axe vertical (

Oz

), de longueur

L_0

et de rayon

r_0

. On étire le cylindre verticalement selon son axe. Sa longueur augmente et son rayon diminue alors, tels que

L = L_0 + ΔL

r = r_0 + Δr

(voir Figure 7).

Extension uniaxiale

Figure 7 - Cylindre en traction. On étire verticalement un cylindre d'axe vertical, de longueur

L_0

et de rayon

r_0

. Après étirement, sa longueur est

L = L_0 + ΔL

et son rayon est

r = r_0 + Δr

.
21. On suppose que le matériau constituant le cylindre est incompressible et donc que le volume du cylindre est conservé. En supposant les déformations petites devant les dimensions du cylindre

(|ΔL|≪L_0

|Δr|≪r_0)

, montrer qu'on a la relation

(Δr)/(r_0) = - ν(ΔL)/(L_0)

où

ν

est une constante dont on déterminera la valeur.

Pour un matériau élastique, linéaire et compressible, la relation $(Δr)/(r_0) = - ν(ΔL)/(L_0)$ est vérifiée et signifie que la déformation du solide dans la direction normale à la contrainte $ε_⊥ = (Δr)/(r_0)$ est proportionnelle à la déformation dans la direction parallèle à la contrainte $ε_‖ = (ΔL)/(L_0)$ de sorte que $ε_⊥ = - νε_‖$ . Le coefficient $ν = - (ε_⊥)/(ε_‖) = - (Δr/r_0)/(ΔL/L_0)$ est appelé coefficient de PoISSON du solide et tel que $- 1 \leq ν \leq 1/2$ .

Pour une contrainte unidirectionnelle normale à la surface

S

(figures 6 (a) et 7 ), la linéarité entre contrainte et déformation dans la direction de la contrainte s'écrit

σ = Eε_‖

(loi de Hooke) où la constante de proportionnalité

E

est appelée module de Young du solide.
22. Déterminer la dimension du module de Young

E

.
23. On considère un cylindre élastique, linéaire et compressible, de coefficient de PoISSON

ν

soumis à la même contrainte d'extension verticale. On suppose toujours que les déformations sont petites devant les dimensions du cylindre. Montrer que la variation relative de volume du cylindre

(ΔV)/(V_0)

est reliée à la contrainte appliquée

σ

par la relation

(ΔV)/(V_0) = (1 - 2ν)/Eσ

où

V_0

est le volume du cylindre avant extension.

III.A.b Contrainte de cisaillement

On considère un solide de forme parallélépipèdique soumis à une contrainte de cisaillement (Figure 6 (b)). Dans ce cas, avec les notations de la figure 6 (b), la contrainte est tangentielle à la surface

S

et la déformation dans la direction de la contrainte est

ε_‖ = Δx/L = tanθ

. Pour de petites déformations, la linéarité entre contrainte et déformation dans la direction de la contrainte se traduit par

σ = Gε_‖

(la loi de Hooke pour le cisaillement) où la constante de proportionnalité

G

est appelée module de cisaillement du solide. La théorie de l'élasticité montre que le module de cisaillement s'écrit en fonction du module de Young et du coefficient de Poisson selon l'expression

G = E/(2(1 + ν))

Pour une petite déformation $(|Δx|≪L)$ , montrer que le module de cisaillement s'écrit $G = F/(Sθ)$ , où $θ$ est l'angle d'inclinaison du solide par rapport à la verticale, défini sur la figure 6 (b).

Figure 8 - Cylindre sollicité en torsion. La surface supérieure d'un cylindre de longueur

L

et de rayon

R

subit une rotation d'angle

φ

par rapport à la surface inférieure du cylindre qui reste fixe. Le moment de torsion crée une contrainte de cisaillement dans le matériau.

Un exemple classique de cisaillement pur est donné par la torsion d'un fil cylindrique (fil de torsion). On considère un fil élastique assimilé à un cylindre de rayon $R$ , de longueur $L$ (avec $L≫R$ ), d'axe vertical ( $Oz$ ) soumis à un moment de torsion $M^(\to-)$ de direction parallèle à l'axe du cylindre. Le plan supérieur du cylindre tourne alors d'un angle $φ$ , supposé petit, par rapport au plan inférieur qui reste fixe (Figure 8). On note $e_z\to$ le vecteur unitaire vertical. On appelle ici O le point d'intersection entre l'axe du cylindre et le plan supérieur du cylindre.

Montrer qu'un élément de surface $dS$ du plan supérieur du cylindre situé à la distance $r$ de l'axe du cylindre subit une déformation $ε(r) = r/Lφ$ selon la direction orthoradiale.
À partir de la loi de Hooke pour le cisaillement et de la définition de la contrainte, déduire l'expression de la norme $dF$ de la force tangentielle $dF^(\to-)$ qui s'exerce sur l'élément de surface $dS$ en fonction de $G, r, L, φ$ , et $dS$ .
Exprimer le moment de torsion $M^(\to-)$ dû aux forces de cisaillement et montrer qu'il se met sous la forme $M^(\to-) = Cφe_z\to$ où la raideur du fil de torsion $C$ est une constante (appelée constante de torsion) qu'on déterminera en fonction de $G, R$ et $L$ .

III.B Indendation par pointe AFM (Atomic Force Microscope)

Une technique expérimentale permettant de sonder les propriétés mécaniques des cellules consiste à indenter (c'est-à-dire déformer localement) la cellule à l'aide d'une pointe de microscope à force atomique (AFM pour 'atomic force microscope'). Les figures 9 et 10 représentent, respectivement, le principe et les résultats d'une expérience d'indentation d'une cellule par une pointe AFM pyramidale. La géométrie de la pointe d'AFM joue un rôle important et détermine le comportement de la force

F(δ)

appliquée sur la cellule en fonction de l'indentation, notée

δ

. On réalise la déformation en restant dans le régime élastique de la cellule (déformation réversible). La pointe AFM est le plus souvent une pyramide à base carrée et, dans ce cas, la dépendance de la force

F(δ)

n'est plus linéaire en fonction de l'indentation, à la différence des exemples abordés au paragraphe III.A. La force varie alors en

δ^2

et s'exprime de la manière suivante :

F(δ) = CE^* δ^2

où

C

est une constante numérique dépendant de la pointe AFM et

E^* = E/(1 - ν^2)

est appelé module de Young effectif.

Lors de l'expérience d'indentation, la pointe AFM est déplacée verticalement selon l'axe (

Oz

) vers le bas grâce à un levier flexible dont la raideur est modélisée par une constante de raideur

k_(AFM)

et dont la hauteur

z

à son extrémité est imposée par le déplacement d'un système piézo-électrique.

On note

ζ

la distance verticale entre la pointe et l'extrémité du levier (aussi appelée déflexion du levier). Au début de l'expérience (

t = 0

), la base de la pointe affleure la cellule sans l'indenter, et l'extrémité du levier et le haut de la pointe sont alignés de sorte que

ζ(0) = z(0) = δ = 0

(voir figure 9). Lorsqu'on indente la cellule (

t > 0

), la force

F(δ)

et l'indentation

δ

augmentent. On suppose que la force exercée par le levier sur la pointe est celle d'un ressort de constante de raideur

k_(AFM)

et d'allongement

ζ

. On donne les caractéristiques suivantes de la pointe et du levier :

C = 0, 5

k_(AFM) = 0, 1 N \cdot m^(- 1)

L'indentation est effectuée à l'aide d'un déplacement du levier à vitesse constante et mesurée grâce à une méthode optique. Une photodiode permet de repérer la position d'un point lumineux issu de la réflexion d'un faisceau laser sur un miroir placé sur la pointe. On néglige les éventuels frottements exercés par le milieu extracellulaire sur la pointe et le levier dans cette expérience.
28. Expliquer, synthétiquement, comment on peut déduire la force

F

appliquée sur la cellule à partir de la mesure de la position de l'extrémité du levier et de la position de la pointe.

Figure 9 - Schéma de principe d'une expérience d'indentation d'une cellule par une pointe AFM. Les schémas représentent l'expérience vue de profil (axe vertical (

Oz

) dirigé vers le bas) à l'instant initial

t = 0

où le levier est horizontal et à un instant

t > 0

au cours de l'indentation où le levier est défléchi. Lorsque le levier est abaissé, la pointe indente la cellule et le levier fléchit. L'extrémité du levier est repérée par sa position verticale

z

. L'indentation

δ

correspond à la position verticale de la pointe. La déflexion du levier est notée

ζ

Figure 10 - Résultats d'une expérience d'indentation d'une cellule par une pointe AFM. Graphe de gauche : courbes force-indentation

F(δ)

pour deux cellules différentes (courbes 1 et 2 respectivement). Graphe de droite : expérience où la cellule est d'abord indentée jusqu'à

1μ m

environ en descendant le levier puis en remontant le levier.

La figure 10 (graphe de gauche) représente les courbes d'indentation $F(δ)$ pour deux cellules différentes correspondant aux courbes 1 et 2 respectivement. Déterminer en exploitant la figure quelle cellule sera la plus déformable. Donner une estimation numérique du module de Young effectif $E^*$ des deux cellules.
On effectue une expérience d'indentation jusqu'à $1μ m$ environ puis on retire la pointe (voir figure 10, graphe de droite). Proposer une explication permettant de comprendre le comportement quasi-statique observé et analyser le rôle éventuel de la dissipation visqueuse.

III.C Oscillations de la pointe AFM et réponse harmonique

On utilise le dispositif expérimental, décrit au paragraphe précédent (figure 9), dans lequel on fait osciller la pointe AFM pour indenter la cellule de manière périodique. L'évolution temporelle du
déplacement piézoélectrique du levier est

z(t) = z_0 + z_1 cos(ωt + φ_z)

où

z_0

est la position moyenne du levier,

z_1

l'amplitude des oscillations,

ω

la pulsation et

φ_z

la phase. On suppose que l'amplitude des oscillations

z_1

est petite devant la position moyenne

z_0

. Dans ces conditions, l'indentation

δ

et la force

F

sont également périodiques de pulsation

ω

et leurs évolutions sont données par :

δ(t) = δ_0 + δ_1 cos(ωt + φ_δ)

F(t) = F_0 + F_1 cos(ωt + φ_F)

. Nous rappelons que la pointe AFM est pyramidale et donc

F(δ) = CE^* δ^2

. Les valeurs moyennes de

F

δ

sont telles que

F(δ_0) = F_0

.
31. En linéarisant l'expression de

F(δ)

aux petites amplitudes d'oscillations, exprimer

F(δ) - F_0

en fonction de

δ - δ_0, C, E^*

δ_0

.
32. En utilisant l'expression du module de cisaillement

G = E/(2(1 + ν))

, vue au paragraphe III.A, et en exprimant

E^*

en fonction de

G

, déduire de la question précédente l'expression du module de cisaillement

G

de la cellule en fonction de

δ_0, C, ν, F - F_0

, et

δ - δ_0

On introduit les notations complexes suivantes

F_-(ω)

δ_-(ω)

associées à

F - F_0

δ - δ_0

telles que :

F_-(ω) = F_1 e^(i(ωt + φ_F)) et δ_-(ω) = δ_1 e^(i(ωt + φ_δ)), F_1, δ_1 \in ℝ^+.

Le module de cisaillement complexe

G_-(ω)

est défini par

G_-(ω) = ((1 - ν))/(4Cδ_0)(F_-(ω))/(δ_-(ω))

. On peut alors écrire une relation linéaire entre les grandeurs complexes

F_-(ω)

δ_-(ω)

, en introduisant une grandeur

Z_-(ω)

analogue à une impédance électrique, telle que :

Z_-(ω) = (F_-(ω))/(δ_-(ω)) = (4Cδ_0)/(1 - ν)G_-(ω).

Le module de cisaillement

G_-(ω)

comporte une partie réelle, notée

G^'(ω)

, et une partie imaginaire, notée

G^('')(ω)

Figure 11 - Réponse harmonique dans une expérience d'oscillations à faible amplitude de la pointe AFM. Courbes représentant la force

F

en fonction de l'indentation

δ

à des fréquences d'oscillations variant entre 0,1 et 100 Hz .

En interprétant l'expression $F_-(ω) = (Re(Z_-(ω)) + iωIm((Z_-(ω))/ω))δ_-(ω)$ comme une relation entre la force $F$ , l'indentation $δ$ et la dérivée de $δ$ par rapport au temps, proposer une interprétation physique des parties réelle et imaginaire de $G_-(ω)$ . Déterminer le déphasage $φ$ entre $F_-(ω)$ et $δ_-(ω)$ en fonction de $G^'(ω)$ et $G^('')(ω)$ et discuter les cas $φ = 0$ et $φ = π/2$ .
La figure 11 représente les résultats expérimentaux d'une expérience d'oscillations d'une pointe AFM pour une amplitude d'oscillation de 50 nm à différentes fréquences. Interpréter le tracé de la force $F$ en fonction de l'indentation $δ$ . Discuter notamment l'évolution de la direction du demi-grand axe des ellipses observées ainsi que leur largeur, en fonction de la fréquence des oscillations. On pourra s'appuyer sur la figure 1 en préambule du sujet.

IV Modèles visco-élastiques du cytoplasme cellulaire

La cellule et son cytoplasme ont à la fois des caractéristiques mécaniques d'un solide élastique (solide de Hooke) et d'un fluide visqueux (fluide newtonien). Ce sont donc des milieux visco-élastiques. Nous introduisons ici les modèles visco-élastiques les plus couramment utilisés en mécanique cellulaire. À partir de la déformation

ε(t)

du matériau et de la contrainte

σ(t)

appliquée au matériau (définies au paragraphe III et sur la figure 6), on définit le module de cisaillement par

G(t) = (σ(t))/(ε(t))

.
Pour caractériser un matériau visco-élastique, on utilise :

⋆

l'expérience dite de fluage où on impose une contrainte constante

σ(t) = σ_0 = C

te et on mesure la déformation en fonction du temps

ε(t)

;

l'expérience de dynamique oscillatoire où la déformation et la contrainte oscillent de façon harmonique à la pulsation $ω, ε(t) = ε_0 cos(ωt)$ et $σ(t) = σ_0 cos(ωt + φ)$ , où $φ$ est le déphasage entre la contrainte et la déformation.

Dans cette partie, nous étudierons d'abord les propriétés des deux éléments de base, le solide de Hooke (représenté par un élément élastique analogue à un ressort de module de Young

E

) et le fluide newtonien (représenté par un élément visqueux analogue à un amortisseur fluide de viscosité

η

) (paragraphes IV.A et IV.B) avant de les combiner dans le modèle visco-élastique de Kelvin-Voigt (paragraphe IV.C) qui associe un élément élastique et un élément visqueux en parallèle (Figure 12).

Figure 12 - Éléments de base pour la modélisation d'un matériau visco-élastique linéaire. De gauche à droite, sont représentés : un élément élastique (solide de Hooke de module de Young

E

), un élément visqueux (fluide newtonien de viscosité

η

) et une association en parallèle de ces deux éléments dans le modèle de Kelvin-Voigt.

IV.A Cas élémentaires : solide de Hooke et fluide newtonien

Nous rappelons qu'un solide de Hooke obéit à la loi du même nom : $σ = Eε$ .

Nous étudions la réponse du solide de Hooke dans une expérience de fluage. On considère un échelon de contrainte, qui a pour expression $σ(t) = σ_0 = C$ te pour toute date $t > 0$ . Déterminer l'expression de $ε(t)$ .

Un fluide newtonien est un matériau purement visqueux. Il est caractérisé par sa viscosité $η$ . On peut montrer que la contrainte $σ$ est proportionnelle au taux de déformation $ε˙ = (dε)/(dt)$ , telle que $σ = ηε˙$ .

Dans le cadre d'une expérience de fluage, on soumet le fluide newtonien à un échelon de contrainte ( $σ(t) = σ_0$ pour tout temps $t > 0$ ). Déterminer l'expression de $ε(t)$ , en supposant la déformation du fluide nulle à $t = 0$ .

IV.B Cas harmonique et aspects énergétiques

On considère maintenant un matériau visco-élastique soumis à une sollicitation où la déformation et la contrainte oscillent à la pulsation

ω

. On utilise la notation complexe :

ε_- = ε_0 e^(iωt)

σ_- = σ_0 e^(i(ωt + φ)) = σ_0 e^(iφ)e^(iωt)

, où

ε_0 \in ℝ

σ_0 \in ℝ

. Le module de cisaillement complexe est

G_- = σ_-/ε_- = G^' + iG^('')

, avec

G^' = Re(G_-)

G^('') = Im(G_-)

. On note

G_0 = |G_-|

le module de

G_-

.
37. Déterminer l'expression de

G_-

pour un solide de Hooke d'une part, puis pour un fluide newtonien d'autre part. En déduire

G^', G^(''), G_0

et le déphasage

φ

pour un solide de Hooke puis pour un fluide newtonien.
38. Définir le travail élémentaire par unité de volume

δw

fourni par une contrainte

σ

au cours d'une déformation infinitésimale

dε

.
39. Déterminer

δw

dans le cas où

ε(t) = ε_0 cos(ωt)

σ(t) = σ_0 cos(ωt + φ)

. On mettra l'expression de

δw

sous la forme suivante :

δw = δw_1 + δw_2

, avec

δw_1 = - ωσ_0 ε_0 cosφcosωtsinωt dt

δw_2 = ωσ_0 ε_0 sinφsin^2 ωt dt

.
40. Exprimer l'énergie volumique de déformation

w_T

fournie à l'unité de volume sur une période

T = (2π)/ω

en fonction de

ε_0, σ_0

φ

. En déduire

w_T

pour un solide de Hooke, puis pour un fluide newtonien. Proposer une interprétation de

w_T

pour un matériau visco-élastique linéaire.
41. Exprimer l'énergie volumique

w_(1, T/4)

, due à

δw_1

, fournie à l'unité de volume sur un quart de période

T/4

et définie par :

w_(1, T/4) = \int_0^(T/4)δw_1

. En déduire

w_(1, T/4)

pour un solide de Hooke et pour un fluide newtonien. Proposer une interprétation de

w_(1, T/4)

pour un matériau visco-élastique linéaire.
42. À partir des questions 40 et 41 , déterminer

tanφ

et en donner une signification à partir de grandeurs énergétiques pour un matériau visco-élastique linéaire.

IV.C Modèles de Kelvin-Voigt d'un fluide visco-élastique

Le modèle de Kelvin-Voigt est un des modèles les plus élémentaires décrivant le comportement d'un matériau visco-élastique linéaire. Il combine les propriétés d'un solide de Hooke et d'un fluide newtonien.

La figure 12 représente le modèle de Kelvin-Voigt sous la forme d'un ressort et d'un amortisseur associés en parallèle et avec les propriétés suivantes :

⋆

le ressort est un élément élastique de module de Young

E

, soumis à la contrainte

σ_(el)

, et subissant la déformation

ε_(el)

;

⋆

l'amortisseur est un élément visqueux de viscosité

η

, de contrainte

σ_(visc)

, et de déformation

ε_(visc)

.
43. Établir l'expression de la déformation

ε_(KV)

et de la contrainte

σ_(KV)

en fonction des paramètres du ressort et de l'amortisseur.
44. Établir l'équation différentielle vérifiée par la contrainte

σ_(KV)

et la déformation

ε_(KV)

du matériau visco-élastique. On fera apparaître un temps de relaxation

τ = η/E

.
45. Dans le cas d'une expérience de fluage (échelon de contrainte

σ_(KV)(t > 0) = σ_0

), déterminer l'évolution de la déformation

ε_(KV)(t)

en résolvant l'équation différentielle de la question 44 et en prenant comme condition initiale

ε_(KV)(t = 0) = 0

. Donner l'allure graphique de

σ_(KV)(t)

ε_(KV)(t)

.
46. Dans le cas d'une déformation et d'une contrainte harmoniques, établir l'équation algébrique, associée à l'équation différentielle de la question 44, à l'aide des notations complexes introduites au paragraphe IV.B. En déduire l'expression du module de cisaillement complexe

G_-_(KV)(ω)

, puis les expressions de

G_(KV)^'(ω), G_(KV)^('')(ω)

φ_(KV)(ω)

en fonction de

ω, E

τ

.
47. Donner l'allure graphique de

G_(KV)^'

G_(KV)^('')

. On fera apparaître une pulsation caractéristique

ω_c

et on analysera les régimes asymptotiques. Discuter le comportement aux temps longs (basses fréquences) et aux temps courts (hautes fréquences) du matériau visco-élastique dans le modèle de Kelvin-Voigt et comparer au solide élastique et au fluide newtonien.
48. À l'aide des résultats des expériences d'indentation par une pointe AFM en réponse harmonique présentés sur la figure 11 (partie III.C, question 34), discuter la validité du modèle de Kelvin-Voigt pour appréhender le comportement mécanique d'une cellule. On pourra s'appuyer sur l'évolution des caractéristiques des ellipses observées dans les expériences en fonction de la fréquence.

1. Cellule comportant un noyau et un cytoplasme.

Pas de description pour le moment

Commentaires• X ENS Physique B PC 2021

Connectez-vous pour participer aux discussions

Partagez vos avis, posez des questions et échangez avec la communauté

Chargement des commentaires...