L’analyse par éléments finis est un outil fréquemment utilisé pour étudier les performances mécaniques des structures sous charge. Ici, nous appliquons son utilisation à la modélisation de la biomécanique de la mâchoire du poisson-zèbre.
Article de méthode
L’analyse par éléments finis est un outil fréquemment utilisé pour étudier les performances mécaniques des structures sous charge. Ici, nous appliquons son utilisation à la modélisation de la biomécanique de la mâchoire du poisson-zèbre.
La morphogenèse squelettique se produit par des comportements cellulaires étroitement régulés au cours du développement ; de nombreux types de cellules modifient leur comportement en réponse à une contrainte mécanique. Les articulations squelettiques sont soumises à des charges mécaniques dynamiques. L’analyse par éléments finis (FEA) est une méthode de calcul, fréquemment utilisée en ingénierie, qui permet de prédire comment un matériau ou une structure réagira à une entrée mécanique. En divisant un système entier (dans ce cas, le squelette de la mâchoire du poisson zèbre) en un maillage d’éléments finis plus petits, l’analyse par éléments finis peut être utilisée pour calculer la réponse mécanique de la structure aux charges externes. Les résultats peuvent être visualisés de plusieurs façons, y compris sous la forme d’une « carte thermique » montrant la position des déformations principales maximales et minimales (une déformation principale positive indique une tension tandis qu’une déformation négative indique une compression). Le maximum et le minimum se réfèrent à la plus grande et à la plus petite souche). Ceux-ci peuvent être utilisés pour identifier les régions de la mâchoire et donc les cellules susceptibles d’être soumises à des charges de tension ou de compression particulièrement élevées lors du mouvement de la mâchoire et peuvent donc être utilisés pour identifier les relations entre la contrainte mécanique et le comportement des cellules. Ce protocole décrit les étapes de génération de modèles d’éléments finis à partir de données d’images confocales sur le système musculo-squelettique, en utilisant la mâchoire inférieure du poisson-zèbre comme exemple pratique. Le protocole guide le lecteur à travers une série d’étapes : 1) coloration des composants musculo-squelettiques, 2) imagerie des composants musculo-squelettiques, 3) construction d’une surface tridimensionnelle (3D), 4) génération d’un maillage d’éléments finis, 5) résolution de l’analyse par éléments finis et enfin 6) validation des résultats par comparaison avec les déplacements réels observés dans les mouvements de la mâchoire du poisson.
Finite Element (FE) la modélisation est une technique d'ingénierie qui peut informatiquement calculer et cartographier l'ampleur et la localisation des souches agissant sur une structure 1. Le modèle est constitué de la structure en 3D, représentée par un maillage de "éléments finis", et le résultat de l'analyse finale est régi par un certain nombre de facteurs, dont la structure et le nombre des éléments du maillage, l'ampleur et la localisation de la mécanique les charges et les propriétés des matériaux. Les propriétés des matériaux décrivent certains aspects du comportement d'un matériau dans un type donné de charge; Le module de Young (E) décrit l'élasticité du matériau tandis que le coefficient de Poisson décrit la diminution proportionnelle de la largeur d'un matériau à sa longueur lorsqu'un échantillon est étiré. la modélisation FE peut être utilisé pour calculer une variété de variables, y compris le déplacement, le stress, la pression et la contrainte agissant sur le modèle en prenant en compte les données d'entrée uniques sur la structure '; Forme s, l'emplacement et l'ampleur des charges et les propriétés des matériaux spécifiques.
La modélisation FE est largement utilisé dans l' ingénierie 2 et de plus en plus pour des applications orthopédiques 3 et paléontologiques 4. Dans le développement des forces biomécaniques sont connus pour agir comme un stimulus dans de nombreuses cellules pour activer les réponses des cellules 5-8 et il est utile de prévoir deux positions relatives et des amplitudes des stimuli mécaniques à l' intérieur le développement de systèmes d'organes, cependant, actuellement de modélisation FE a été peu utilisé pour le développement du poisson zèbre.
Deux cartilages et des os se sont avérés être des matières mécanosensibles. Par exemple, une compression in vitro a été trouvé pour activer les voies chondrogéniques, tandis que la tension a été démontré que nécessaire pour la formation d'os neuf. FE analyse (FEA) a été exploitée pour modéliser les souches agissant sur des échantillons biologiques, y compris ceux qui agissent sur les éléments du squelette lors de l'os formation 10. D' autres applications de développement comprennent son utilisation pour prédire la forme d'un joint après qu'il a été exposé à des forces biomécaniques théoriques 11,12 et de montrer le modèle des souches présentes au cours de poussin genou morphogenèse joint 8.
Ce protocole vise à partager l'expérience de générer des surfaces en 3 dimensions, mailles et modèles d'éléments finis à partir d'images confocale en vue de comprendre la mécanique des tissus en développement. Nous montrons aussi des moyens de valider les modèles FE si la capture réelle information de déplacement conjoint in vivo. Alors que nous utilisons la mâchoire zebrafish en tant qu'exemple les mêmes techniques peuvent être utilisées sur un petit système biologique pour lequel des informations sur la structure 3D du système musculo-squelettique peut être obtenu par imagerie confocale ou multiphotonique.
Toutes les étapes du protocole suivent l'Université du soin des animaux de Bristol et les lignes directrices de bien-être et ceux de la Home Office du Royaume-Uni.
1. Visualisation de l'appareil locomoteur Anatomie
NOTE: Pour visualiser la forme des éléments du squelette, de quantifier le muscle et d'identifier l'emplacement exact des attaches musculaires, immunocoloration (section 1.1) du poisson à l'âge approprié pour la myosine squelettique (qui révèle le muscle) et collagène de type II (pour visualiser cartilage). Vous pouvez également visualiser l'anatomie musculo - squelettique en utilisant des lignes de rapporteurs fluorescents transgéniques tels que le collagène a1 reporter COL2A1: mCherry 13,14 pour visualiser le cartilage et la lente chaîne lourde de myosine reporter smyhc: GFP 15 pour visualiser la position des attaches musculaires (section 1.2).
lignes alternatives qui marquent le cartilage et les muscles peuvent fonctionner aussi bien.
2. Création d'une surface 3D
3. Calcul des Forces musculaires à être utilisé dans le modèle FE
4. Génération d'un maillage
5. Finite Element Model Construction
6. Validation des Jaw Déformation / Déplacement Distances
Immunocoloration pour le muscle (figure 1A) et le cartilage (figure 1B) ou l' imagerie de reporters transgéniques (figure 1C) permet à la structure 3D de la mâchoire à visualiser, ainsi que la masse musculaire associée. En imagerie à haute résolution , il était possible de construire un modèle qui reflète à la fois la forme en trois dimensions de la mâchoire (Figure 2) et l'emplacement et l' emplacement des charges (Figure 3). Utilisant des déplacements in vivo vus à travers la capture vidéo haute vitesse (Figure 4) nous avons vérifié que l'amplitude du mouvement dans le modèle était dans une fourchette réaliste.
Les modèles FE une fois l' exécution peuvent être utilisés pour afficher une série de données, tels que le stress (figure 5A), minimum et déformation principale maximale (figure 5B - K). Ces résultats sont trois dim ensional de sorte que le modèle peut être agrandie pour voir des motifs fins de détail (Figure 5E, 5I) tourné pour obtenir des points de vue pertinents (Figure 5F, 5G, 5J, 5K) et numériquement sectionnés (Figure 5E ', 5E' ', 5I', 5I '') pour montrer comment les modèles de contrainte, déformation ou changement de pression dans l'ensemble du modèle. Il est également possible d'extraire des données quantitatives à partir du modèle (non représenté). En vérifiant le modèle et en utilisant les propriétés des matériaux les plus précis, les charges les mailles et la forme du modèle FE peut être utilisé pour explorer la meilleure estimation de l'environnement mécanique vécue par les cellules pendant cette fenêtre de développement. Les résultats du modèle peuvent être directement comparés aux changements dans le comportement cellulaire et l' expression des gènes 20.
re 1 "src =" / files / ftp_upload / 54811 / 54811fig1.jpg "/>
Figure 1:. Des images représentatives des éléments musculo - squelettiques de la mâchoire inférieure zebrafish à 5 dpf piles confocaux représentatifs de la mâchoire inférieure des larves 5dpf tous montré avec antérieure vers le haut (A) immunocoloration pour A4.1025 qui colore tout myosine squelettique (B) immunocoloration pour collagène de type II qui marque tout le cartilage (C) Stack à partir d' une larve vivant exprimant les journalistes transgéniques COL2A1: GFP muscle lent (vert): cartilage (rouge) et smyhc marquage mCherry. IA: intermandibularis antérieure, PH: protractor hyoïdien, AM: adductor mandibulae, HI: hyoïdien inférieure, HI: hyoïdien supérieure, CH: sternohyoideus, MC: cartilage de Meckel, PQ: Palatoquadrate, CH:. Ceratohyal S'il vous plaît cliquer ici pour voir une version plus grande de ce chiffre.


Figure 3: maillages représentatifs montrant les contraintes et les vecteurs de force maillages représentatifs pour une larve de 5 dpf pour (A) la fermeture de la bouche et.(B) ouverture de la bouche. Les points blancs représentent les endroits où le modèle est contraint et dans lequel les dimensions (par exemple, x et y ou x, y et z). des lignes blanches indiquent les positions des muscles, avec le vecteur de la force musculaire représenté par des flèches blanches. Rouge montre cartilage et jaune l'interzone. Ce chiffre a été modifié depuis le matériel supplémentaire précédemment publié dans Brunt et al. 15. S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 4:. Tests de sensibilité FE-modèle simulant la mâchoire déplacement dans zebrafish 5dpf pour le cartilage différent et des modules de Young interzones. Mâchoires déplacement (ouvert à fermé en um) est marqué sur la mâchoire; enregistré à l'aide de la touche de couleur. Chaque modèle (A - L) a une combinaison différente de cartilage (c = 1,1, 3,1 ou 6,1 MPa) ou interzone (i = de 0,25, 0,5, 0,75, ou 1 MPa) propriétés. Horizontal flèche noire met en évidence la valeur de la mâchoire de déplacement à l'extrémité du cartilage de Meckel (représenté par la flèche noire verticale). M et N images fixes à partir de vidéos 5 dpf larves montrant au minimum, à savoir, la mâchoire fermée (M) et maximum, à savoir, mâchoire complètement ouvert (N) avec les deux superposée (O) - ligne blanche sur O représente le déplacement (de 43 um). Dans ce cas , par rapport propriétés du cartilage de 1,1 avec une interzone de 0,25 (A) meilleur match les déplacements observés dans les poissons vivants (O). Panneaux AL de cette figure ont été publiées antérieurement dans Brunt et al. 15. S'il vous plaît cliquer ici pour voir un plus grand version de ce chiffre.

Figure 5: Les données représentatives des modèles EF simulation FE-modèle de tous les muscles appliqués dans un 5 dpf larve (A - C).. (A) Von Mises (EMaxmin) (B) de déformation minimum principal (E min. P, μɛ) (C) déformation principale maximale (E max. P., μɛ). simulation FE-modèle de souches principales maximales et minimales lors de l'ouverture de la mâchoire. (D - K): déformation principale maximale (. E Max P., μɛ) en (D) vue de la mâchoire ventral et (E) ventral vue commune (E) montre l' emplacement des sections proximale-distale à travers le cartilage de Meckel articulaire et l'interzone (E ') et (E'), respectivement. (F): latéralmâchoire vue. (G): vue latérale commune. (H - K): la souche principale minimale (. E Min P, μɛ) dans (H) vue de la mâchoire et ventrale (I) ventrale de vision commune. (I) montre l' emplacement des sections proximale-distale à travers l'articulation du cartilage de Meckel et interzones en (I ') et (I' ') respectivement (J): vue latérale de la mâchoire. (K): vue latérale commune. Ce chiffre a été publié précédemment dans Brunt et al. 15. S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.
| Nombre de fibres musculaires | Zone de fibres musculaires (pm 2) | zone de groupe musculaire (pm2) | Force de (N) | |
| 5 dpf intermandibularis antérieure | 5 | 23,8 | 119 | 4.76e-6 |
| 5 dpf protractor hyoïdien | 6 | 23,8 | 142,8 | 5.71e-6 |
| 5 dpf adductor mandibulae | 9 | 23,8 | 214,2 | 8.57e-6 |
Tableau 1: Muscle quantification. Calculées forces musculaires moyennes pour la Intermandibularis Anterior, Adductor mandibule et Protractor hyoïdien à 5 dpf utilisant 40 nN / um 2 (valeur par unité de surface prise de la référence 17). (Lorga et al., 2011) (n = 3).
Modèles d'éléments finis ont été utilisés pour relier les zones d'éléments squelettiques qui sont mis à rude épreuve avec ceux subissant la formation d'os 10, ainsi que pour cartographier les zones sous contrainte pendant l' ossification endochondrale et morphogenèse commune 8,12,21. D' autres études ont également été en mesure d'appliquer des modèles théoriques de croissance pour répliquer les modifications au cours du développement conjoint 11,12. Ici , nous montrons le protocole pour la construction de modèles FE pour un système relativement simple, la mâchoire zebrafish 20. Contrairement à d' autres méthodes de collecte des images brutes pour les modèles FE, comme la tomodensitométrie 22, imagerie confocale des lignées transgéniques ou zebrafish immunocolorées permet de multiples tissus à étudier. Il peut, par conséquent, fournir des informations directes sur les points d'attache des muscles par rapport au cartilage. Parmi les modèles vertébrés zebrafish sont particulièrement favorables à la manipulation génétique et pharmacologique. La génération de modèles FE pour le poisson zèbrecartilage crânio ouvre maintenant la possibilité de poursuivre l'étude de l'interaction entre la biomécanique et de la génétique dans la morphogenèse commune.
Il y a un certain nombre d'étapes critiques du processus de création d'un modèle FE; la première est la génération d'une représentation tridimensionnelle précise du système. Cela nécessite l'imagerie à résolution suffisamment élevée pour définir clairement les limites. Notez que même avec des images haute résolution pour faire une bonne surface ne peut avoir pour lisser certaines régions. Une autre étape cruciale consiste à définir le placement correct de la charge et des contraintes correctes. Un modèle insuffisamment contraint échouera à résoudre et placement incorrect des charges provoque des mouvements anormaux.
Une partie de traitement des données brutes (figure 2) est nécessaire en tant que surface générée à partir des données brutes serait difficile à engrener (figure 2B). Nous filtré les données en utilisant un filtre gaussien (figure 2C ) et nous avons effectué un certain lissage manuel des courbes pour produire un ensemble de lignes propres qui peuvent être convertis en une surface 3D. Trop lissage peut produire une surface "fondu" qui a perdu beaucoup de ses caractéristiques. Le choix de la taille de l'élément correct est un processus itératif que le choix de trop petite taille de l'élément crée un trop grand maillage qui est de calcul intensif. Toutefois, en choisissant une trop grande taille de l'élément va produire un maillage qui ne parvient pas à récapituler la forme correcte de la structure. L'exactitude maillage avait la plus petite taille de l'élément qui a capturé la forme correcte de la mâchoire et convergés sur une solution correcte, vérifiée en utilisant le déplacement de la mâchoire. Il peut également être nécessaire de modifier les propriétés des matériaux ou des calculs de charge pour mieux émuler le bon déplacement de différents âges et les espèces ont des propriétés sensiblement différentes.
Il est important de se rappeler qu'il ya toujours des limites à un modèle hypothétique et unssumptions fait pour exécuter des modèles FE. Lorsque seule la modélisation ou un petit nombre d'échantillons, il est essentiel de veiller à ce qu'un échantillon représentatif est choisi comme il est probable que de petites variations entre les individus. Comme seule une partie des éléments de mâchoire et les muscles ont été inclus, le modèle est une version simplifiée du système musculo-squelettique crânio poisson zèbre. Par conséquent, les contraintes devaient être positionnées pour expliquer où les éléments de mâchoire modélisés devraient se connecter avec le reste du crâne et le modèle a été contraint artificiellement dans le centre de le fixer dans un «espace». Cette contrainte artificielle n'a pas d'incidence sur l'interprétation tirée des modèles comme ceratohyal se n'a pas été analysé. L'inclusion de plus de la structure craniofaciale, en particulier les autres muscles de la mâchoire d' ouverture tels que les sternohyals et son cartilage joint 23, aurait pu ajouter au modèle, mais les limites comprennent la capacité des plus grands modèles de fonctionner dans le logiciel Finite Element.
s = "jove_content"> Une autre limitation est que nous n'avons pas modélisé l' insertion du ligament, bien que cela pourrait être réalisé par l'insertion de ressorts 8. Une autre hypothèse faite dans cette affaire était que le modèle se comporterait de façon linéaire. Les amplitudes des tensions sur les modèles étaient comparables à ceux des modèles publiés et appliqués à des cellules in vitro 10,24, avec les souches étant au- dessous et au- dessus 3,500 -5,000 μɛ en dehors des points de contrainte et de fixation du muscle. Par conséquent, les souches dans les régions concernées du modèle ont été considérées dans une gamme acceptable pour un modèle linéaire. Cartilages ne se comporte pas tout à fait comme un matériau linéaire et a déjà été modélisé en tant que matériau poroélastique, ce qui a permis d' analyser le comportement du fluide dans le modèle 25. Répandre les points d'attache des muscles parmi un groupe de noeuds locaux distribuerait les forces de pointe et représentent plus précisément l'insertion du muscle pour certains muscles.
ent "> L'utilisation de FE permet une évaluation des déformations et des contraintes agissant sur une structure. En tant que technique, il est fréquemment utilisé dans de nombreuses disciplines de sciences biologiques, y compris l'orthopédie, la paléontologie et plus récemment la biologie du développement. Nous décrivons ici comment construire FEs pour le poisson zèbre mâchoire inférieure. dans l'avenir, ces modèles pourrait être étendue à regarder l'ensemble mâchoire, y compris le palais. Des techniques similaires peuvent être utilisées pour modéliser la biomécanique de la colonne vertébrale chez les poissons, ce qui à ce jour ont surtout été étudiés par des moyens cinématiques.Les auteurs n'ont rien à dévoiler.
Quelques données sur les figures 3-5 a été tiré à part de J.Biomech, 48 (12), Brunt et al., La modélisation par éléments finis prédit des changements de forme conjointe et le comportement des cellules due à la perte de la souche musculaire dans le développement de la mâchoire, 3112-22. 2015, avec la permission de Elsevier 15.
LHB a été financée par le programme Cell PhD Wellcome Trust dynamique; KAR a été financé par la MRC subvention de projet MR / L002566 / 1 (attribué à EJR et CLH) et CLH a été financée par Aruk subvention 19479. Nous tenons également à remercier l'établissement Wolfson Bioimaging pour obtenir des conseils d'imagerie.
| Nom | Entreprise | Numéro de catalogue | Commentaires |
|---|---|---|---|
| Coll2 | Abcam | ab34712 | Anticorps anti-collagène de type II - colore tout le cartilage |
| A4.1025 / MF20 | Études développementales banque d’hybridomes | A4.1025 | Anticorps mysoin squelettique - marque tous les muscles squelettiques ; |
| Agarose à bas point de fusion | Sigma ; | A9414-5G | Pour le montage du poisson-zèbre |
| MS222 (Ethyl 3-aminobenzoate methanesulfonate ) | Sigma | E10521-10G | Pour faire un anesthésique |
| Trypsine | Fisher | T/3760/48 | perméabilisation de l’échantillon |
| Dylight 488 IgG de souris | Thermofisher | 35502 | Anticorps secondaire |
| Dylight 550 IgG de lapin | Thermofisher ; | 84541 | Anticorps secondaire |
| SP8/SP5 ou SPE confocal | Leica ; | Pour l’imagerie ; | |
| LAS Logiciel de capture Leica | Logiciel d’imagerie | Leica | |
| Aviso (version 7.0.0) | Logiciel d’analyse d’images 3D | FEI Visualization Science Group | |
| Hypermesh faisant partie du package Hyperworks (version 10) | Logiciel de génération de modèles Altair Engineering | ||
| Abaqus (version 6.14) | SIMULIA | FE (Section 5.7 et 5.8) |