El análisis de elementos finitos es una herramienta de uso frecuente para investigar el rendimiento mecánico de estructuras bajo carga. Aquí aplicamos su uso para modelar la biomecánica de la mandíbula del pez cebra.
Artículo de método
El análisis de elementos finitos es una herramienta de uso frecuente para investigar el rendimiento mecánico de estructuras bajo carga. Aquí aplicamos su uso para modelar la biomecánica de la mandíbula del pez cebra.
La morfogénesis esquelética ocurre a través de comportamientos celulares estrictamente regulados durante el desarrollo; muchos tipos de células alteran su comportamiento en respuesta a la tensión mecánica. Las articulaciones esqueléticas están sometidas a cargas mecánicas dinámicas. El análisis de elementos finitos (FEA) es un método computacional, utilizado con frecuencia en ingeniería, que puede predecir cómo responderá un material o estructura a la entrada mecánica. Al dividir un sistema completo (en este caso, el esqueleto de la mandíbula del pez cebra) en una malla de "elementos finitos" más pequeños, el FEA se puede utilizar para calcular la respuesta mecánica de la estructura a las cargas externas. Los resultados se pueden visualizar de muchas maneras, incluso como un "mapa de calor" que muestra la posición de las deformaciones principales máxima y mínima (una deformación principal positiva indica tensión, mientras que una negativa indica compresión. El máximo y el mínimo se refieren a la mayor y menor cepa). Estos se pueden utilizar para identificar qué regiones de la mandíbula y, por lo tanto, qué células probablemente estén bajo cargas de tensión o compresión particularmente altas durante el movimiento de la mandíbula y, por lo tanto, se pueden usar para identificar relaciones entre la deformación mecánica y el comportamiento de la célula. Este protocolo describe los pasos para generar modelos de elementos finitos a partir de datos de imágenes confocales en el sistema musculoesquelético, utilizando la mandíbula inferior del pez cebra como ejemplo práctico. El protocolo guía al lector a través de una serie de pasos: 1) tinción de los componentes musculoesqueléticos, 2) obtención de imágenes de los componentes musculoesqueléticos, 3) construcción de una superficie tridimensional (3D), 4) generación de una malla de elementos finitos, 5) resolución del FEA y finalmente 6) validación de los resultados mediante la comparación con los desplazamientos reales observados en los movimientos de la mandíbula del pez.
Modelización de elementos finitos (FE) es una técnica de ingeniería que computacionalmente puede calcular y asignar la magnitud y la localización de las cepas que actúan sobre una estructura 1. El modelo consta de la estructura 3D, representada por una malla de "elementos finitos", y el resultado final del análisis se rige por una serie de factores incluyendo la estructura y número de elementos de la malla, la magnitud y la ubicación de la mecánica cargas y las propiedades del material. Propiedades de los materiales describen ciertos aspectos del comportamiento de un material bajo un determinado tipo de carga; el módulo de Young (E) describe la elasticidad del material, mientras que la relación de Poisson describe la disminución proporcional en la anchura de un material a su longitud cuando se estira una muestra. modelado FE se puede utilizar para calcular una variedad de variables tales como el desplazamiento, el estrés, la presión y la tensión que actúa sobre el modelo teniendo en cuenta los datos de entrada único acerca de la estructura '; S forma, ubicación y magnitud de las cargas y de las propiedades de los materiales específicos.
Modelado de FE es ampliamente utilizado en la ingeniería 2 y 3 cada vez más para ortopédicos y paleontológicos aplicaciones 4. En el desarrollo se conocen fuerzas biomecánicas para actuar como un estímulo en muchas células para activar respuestas de las células 5-8 y es útil para predecir tanto las posiciones relativas y las magnitudes de los estímulos mecánicos dentro de desarrollo de sistemas de órganos, sin embargo, en la actualidad modelado FE ha sido poco utilizado para el desarrollo de pez cebra.
Tanto el cartílago y el hueso se ha demostrado que ser materiales mechanosensitive. Por ejemplo, la compresión in vitro se ha encontrado para activar las vías condrogénicas, mientras que la tensión se ha demostrado que es necesario para la formación de hueso 9. análisis FE (FEA) se ha explotado para modelar cepas que actúa en muestras biológicas, incluyendo las que actúan sobre los elementos esqueléticos durante hueso formación 10. Otras aplicaciones de desarrollo incluyen su uso para predecir la forma de una articulación después de que ha sido expuesto a fuerzas biomecánicas teóricos 11,12 y para mostrar el patrón de las cepas presentes durante polluelo rodilla morfogénesis conjunta 8.
Este protocolo está dirigido a compartir la experiencia de la generación de superficies de 3 dimensiones, mallas y modelos de elementos finitos de imágenes confocales con miras a la comprensión de la mecánica de los tejidos en desarrollo. También ponen de manifiesto formas de validar los modelos FE pesar de la captura de información conjunta desplazamiento real en vivo. Mientras que utilizamos la mandíbula de pez cebra como un ejemplo de las mismas técnicas podrían utilizarse en cualquier sistema biológico pequeño para el que la información en 3D sobre la estructura del sistema musculoesquelético se puede obtener mediante formación de imágenes confocal o multifotónica.
Todos los pasos dentro del protocolo siguen la Universidad de cuidado de los animales de Bristol y directrices para el bienestar y los del Ministerio del Interior del Reino Unido.
1. La visualización de la anatomía musculoesquelética
NOTA: Para visualizar la forma de los elementos esqueléticos, para cuantificar el músculo y para identificar la ubicación exacta de las inserciones musculares, immunostain (sección 1.1) de pescado a la edad apropiada para la miosina esquelética (que revela muscular) y colágeno tipo II (para visualizar cartílago). Alternativamente, visualizar la anatomía musculoesquelético usando líneas indicadoras fluorescentes transgénicos tales como el colágeno COL2A1 reportero a1: mCherry 13,14 para visualizar el cartílago y la miosina lenta cadena pesada smyhc reportero: GFP 15 para visualizar la posición de las inserciones musculares (Sección 1.2).
líneas alternativas que marcan el cartílago y músculo podrían funcionar igualmente bien.
2. Generación de una superficie 3D
3. Cálculo de las fuerzas de los músculos para ser utilizado en el modelo FE
4. Generación de una malla
5. Elementos Finitos ModEL Construcción
6. Validación de la mandíbula Deformación / Desplazamiento Distancias
La inmunotinción para muscular (Figura 1A) y el cartílago (Figura 1B) o formación de imágenes de reporteros transgénicos (Figura 1C) permite que la estructura 3D de la mandíbula para ser visualizada, junto con la musculatura asociada. Por formación de imágenes a alta resolución que era posible construir un modelo que captura tanto la forma tridimensional de la mandíbula (Figura 2) y la ubicación y la colocación de las cargas (Figura 3). La utilización de los desplazamientos in vivo visto a través de la captura de vídeo de alta velocidad (Figura 4) se verificó que el rango de movimiento en el modelo estaba dentro de un rango realista.
Los modelos FE Una vez ejecutado se pueden utilizar para mostrar una variedad de datos, como el estrés (Figura 5A), mínimo y máximo de deformaciones principales (Figura 5B - K). Estos resultados son de tres dim ensional lo que el modelo se puede ampliar para ver patrones finos de detalle (Figura 5E, 5I) hace girar para obtener puntos de vista relevantes (Figura 5F, 5G, 5J, 5K) y digital seccionadas (Figura 5E ', 5E' ', 5I', 5I '') para mostrar cómo los patrones de tensión, deformación o cambio de presión en todo el modelo. También es posible extraer los datos cuantitativos a partir del modelo (no mostrado). Al verificar el modelo y el uso de las propiedades de los materiales más precisos, cargas y malla forma el modelo FE puede ser utilizado para explorar la mejor estimación del entorno mecánico experimentado por las células durante esa ventana de desarrollo. Los resultados del modelo se pueden comparar directamente con los cambios en el comportamiento celular y la expresión génica 20.
re 1 "src =" / files / ftp_upload / 54811 / 54811fig1.jpg "/>
Figura 1:. Imágenes representativas de los elementos del aparato locomotor de la mandíbula inferior del pez cebra a 5 dpf pilas confocal representativos de la mandíbula inferior de las larvas 5dpf muestra con toda anterior al inicio (A) La inmunotinción para A4.1025 que tiñe toda la miosina esquelética (B) la inmunotinción para el colágeno tipo II que marca todo el cartílago (C) Pila de una larva viva expresión de los reporteros COL2A1 transgénicos: mCherry marcado cartílago (rojo) y smyhc: GFP músculo lento (verde). IA: intermandibularis anterior, PH: transportador hioideo, AM: mandibular aductor, Hawai: hioideo inferior, Hawai: hioideo superior, CH: esternohioideo, MC: cartílago de Meckel, PQ: palatoquadrate, CH:. Ceratohial Haga clic aquí para ver una versión más grande de esta cifra.


Figura 3: mallas representativo que muestra las limitaciones y los vectores de fuerza mallas representativos para una larva 5 dpf para (A) de cierre y la boca.(B) abertura de la boca. Los puntos blancos indican los lugares donde se ve limitada el modelo y en el que las dimensiones (por ejemplo, x e y o x, y y z). Las líneas blancas indican las posiciones de músculos, con el vector de la fuerza muscular indicado por las flechas blancas. Rojo muestra el cartílago y la interzona amarillo. Esta cifra ha sido modificado a partir de material complementario publicado previamente en Brunt et al. 15. Haga clic aquí para ver una versión más grande de esta figura.

Figura 4:. Pruebas de sensibilidad FE-modelo que simula el desplazamiento de la mandíbula en el pez cebra 5dpf para diferentes cartílago y interzonales módulos de Young. La mandíbula de desplazamiento (abierto a cerrado en micras) está marcado en la mandíbula; grabada utilizando la tecla de color. Cada modelo (A - L) tiene una combinación diferente de cartílago (c = 1,1, 3,1, o 6,1 MPa) o 1 MPa propiedades) interzona (i = de 0,25, 0,5, 0,75, o. Horizontal flecha negro resalta el valor del desplazamiento de la mandíbula en la punta del cartílago de Meckel (representado por la flecha negro vertical). M y N fotogramas de vídeos de 5 dpf larvas que presenten mínimo, es decir, la mandíbula cerrada (M) y el máximo, es decir, la mandíbula completamente abierta (N) con los dos superpuesta (O) - línea blanca en O representa el desplazamiento (de 43 micras). En este caso las propiedades del cartílago relativa de 1,1 con una interzona de 0,25 (A) la mejor partido de los desplazamientos observados en peces vivos (O). Sobre placas de aluminio de esta figura se han publicado previamente en Brunt et al. 15. Haga clic aquí para conocer el vers más grandesión de esta figura.

Figura 5: Los datos representativos de los modelos de simulación FE FE-modelo de todos los músculos aplicados en un 5 dpf larva (A - C).. (A) de Von Mises (EMaxmin) (B) cepa principal mínimo (E Min. P, μɛ) (C) la tensión principal máxima (E máx. P., μɛ). simulación con EF-modelo de cepas máximos y mínimos principales durante la apertura de la mandíbula. (D - K): la tensión principal máxima (. E Max P., μɛ) en (D) Muestra de la mandíbula ventral y (E) vista conjunta ventral (E) muestra la ubicación de las secciones proximal-distal a través del cartílago de Meckel conjunta y la interzona en (e ') y (e' '), respectivamente. (F): lateralla mandíbula vista. (G): vista lateral conjunta. (H - K): la tensión principal mínima (. E Min P, μɛ) en (H) vista ventral y la mandíbula (I) vista conjunta ventral. (I) muestra la ubicación de las secciones proximal-distal a través del cartílago articular de Meckel y en interzona (I ') y (I' ') respectivamente (J): Vista lateral de la mandíbula. (K): vista lateral conjunta. Esta cifra ha sido publicado previamente en Brunt et al. 15. Haga clic aquí para ver una versión más grande de esta figura.
| Número de fibras musculares | Área de la fibra muscular (m 2) | área de grupo muscular (micras2) | Fuerza (N) | |
| 5 dpf intermandibularis anterior | 5 | 23.8 | 119 | 4.76e-6 |
| 5 dpf transportador hioideo | 6 | 23.8 | 142,8 | 5.71e-6 |
| 5 dpf aductor mandibular | 9 | 23.8 | 214,2 | 8.57e-6 |
Tabla 1: La cuantificación del músculo. Calculado promedio fuerzas musculares para el Intermandibularis anterior, aductor mandibular y prolongador hioideo a 5 dpf utilizando 40 nN / m 2 (valor por unidad de superficie tomado de referencia 17). (Lorga et al., 2011) (n = 3).
Modelos de elementos finitos se han utilizado para relacionar las zonas de elementos esqueléticos que están bajo tensión con las que experimenta la formación de hueso 10, así como para mapear las áreas bajo tensión durante la osificación endocondral y la morfogénesis conjunta 8,12,21. Otros estudios también han sido capaces de aplicar los modelos teóricos de crecimiento para replicar los cambios durante el desarrollo conjunto 11,12. Aquí se muestra el protocolo para la construcción de modelos de FE para un sistema relativamente simple, la mandíbula de pez cebra 20. A diferencia de otros métodos de recopilación de imágenes en bruto para los modelos FE, tales como la TC 22, la imagen confocal de líneas transgénicas de pez cebra o inmunotinción permite múltiples tejidos para ser estudiados. Puede, por lo tanto, proporcionar información directa sobre los puntos de unión de músculos en relación con el cartílago. Entre los modelos de vertebrados pez cebra son particularmente susceptibles a la manipulación genética y farmacológica. La generación de modelos de FE para el pez cebracartílago craneofacial ahora se abre la posibilidad de un mayor estudio de la interacción entre la biomecánica y la genética en la morfogénesis conjunta.
Hay una serie de pasos críticos para el proceso de creación de un modelo de elementos finitos; la primera es la generación de una representación tridimensional precisa del sistema. Esto requiere formación de imágenes a una resolución lo suficientemente alta como para definir claramente los límites. Tenga en cuenta que incluso con imágenes de alta resolución para hacer una buena superficie uno puede tener que suavizar algunas regiones. Otro paso crítico es la definición de la colocación correcta de la carga y las limitaciones correctos. Un modelo insuficientemente limitada fallará para resolver y la colocación incorrecta de las cargas provocará un movimiento anormal.
Algunos de procesamiento de los datos en bruto (Figura 2) es necesario como una superficie generada a partir de los datos en bruto sería difícil de malla (Figura 2B). Filtramos los datos usando un filtro de Gauss (Figura 2C ) y que llevan a cabo algunos de suavizado manual de las curvas para producir un conjunto de contornos limpios que se pueden transformar en una superficie 3D. El exceso de suavizado puede producir una superficie "fundido" que ha perdido muchas de sus características. Elegir el tamaño del elemento correcta es un proceso iterativo como elegir demasiado pequeño un tamaño de elemento demasiado grande crea una malla que es computacionalmente intensivas. Sin embargo, elegir demasiado grande un tamaño de elemento producirá una malla que no recapitular la forma correcta de la estructura. La malla correcta tenía el tamaño más pequeño elemento que capturó la forma correcta de la mandíbula y se reunieron en una solución correcta, verifica por medio del desplazamiento de la mandíbula. También puede ser necesario modificar las propiedades del material o de los cálculos de carga para emular mejor el desplazamiento correcto como diferentes edades y especies tendrán propiedades sustancialmente diferentes.
Es importante recordar que siempre hay limitaciones a un modelo hipotético y unassumptions hizo correr modelos FE. Cuando sólo el modelado de uno o un pequeño número de muestras es crítica para asegurar que una muestra representativa se elige ya que es probable que sean pequeñas variaciones entre los individuos. Ya que sólo se incluyeron algunos de los elementos de mandíbula y los músculos, el modelo es una versión simplificada del sistema músculo-esquelético craneofacial pez cebra. Por lo tanto, las restricciones tuvieron que ser posicionado para tener en cuenta cuando los elementos de mordaza modelados conectarían con el resto del cráneo y el modelo fueron limitadas artificialmente en el centro para fijarlo en el "espacio". Esta restricción artificial no tuvo ningún impacto sobre la interpretación extraída de los modelos como el ceratohial en sí no fue analizado. La inclusión de más de la estructura craneofacial, especialmente otros músculos de la mandíbula de apertura como los sternohyals y su cartílago adjunta 23, podría haber añadido al modelo, pero las limitaciones incluyen la capacidad de los modelos más grandes para funcionar en el software de elementos finitos.
s = "jove_content"> Otra limitación es que no hemos modelado inserción de ligamentos, aunque esto podría lograrse mediante la inserción de los resortes 8. Otra suposición hecha en este caso fue que el modelo se comportaría de forma lineal. Las magnitudes de las tensiones en los modelos eran comparables a los de los modelos publicados y se aplicaron a células in vitro 10,24, con cepas de estar por debajo de 3500 y por encima de -5,000 μɛ aparte de puntos de restricción y de unión de músculos. Por lo tanto, las cepas en las regiones relevantes del modelo se consideraron dentro de un rango aceptable para un modelo lineal. El cartílago no se comporta como un material enteramente lineal y previamente se ha modelado como un material poroelástico, lo que permitió el análisis del comportamiento del fluido en el modelo de 25. La difusión de los puntos de unión de músculos entre un grupo de nodos locales distribuiría los picos de fuerza y representar con mayor precisión la inserción del músculo de ciertos músculos.
ent "> El uso de FE permite una evaluación de las deformaciones y tensiones que actúan sobre una estructura. Como una técnica que se utiliza con frecuencia en muchas disciplinas de ciencias biológicas, incluyendo la ortopedia, la paleontología y más recientemente la biología del desarrollo. A continuación se describe cómo construir FE para el pez cebra la mandíbula inferior. en el futuro estos modelos podría extenderse a mirar a toda la mandíbula, que incluye el paladar. técnicas similares se podrían utilizar para modelar la biomecánica de la columna vertebral en el pescado, que hasta la fecha la mayoría han sido estudiados por medios cinemáticos.Los autores no tienen nada que revelar.
Algunos datos de las Figuras 3-5 se ha reproducido de J.Biomech, 48 (12), Brunt et al., Modelización de elementos finitos predice cambios de forma conjunta y comportamiento de las células debido a la pérdida de la tensión muscular en el desarrollo de la mandíbula, 3112-22. de 2015, con permiso de Elsevier 15.
LHB fue financiado por el programa de doctorado de la célula Wellcome Trust dinámico; KAR fue financiado por la subvención del proyecto MRC MR / L002566 / 1 (otorgado a EJR y CLH) y CLH fue financiado por Aruk subvención 19479. También nos gustaría dar las gracias a la instalación de Wolfson Bioimagen para el asesoramiento de imagen.
| Nombre | Empresa | Número de catálogo | Comentarios |
|---|---|---|---|
| Coll2 | Abcam | ab34712 | Anticuerpo de colágeno tipo II: tiñe todo el cartílago |
| A4.1025 / MF20 | Banco de hibridomas de estudios del desarrollo | A4.1025 | Anticuerpo de misoína esquelética: marca todo el músculo esquelético |
| Agarosa de bajo punto de fusión | Sigma | A9414-5G | Para el montaje de pez cebra |
| MS222 (3-aminobenzoato de etilo metanosulfonato ) | Sigma | E10521-10G | Para hacer anestésico |
| Tripsina | Fisher | T/3760/48 | permeabilización de muestras |
| Dylight 488 Ratón IgG | Thermofisher | 35502 | Anticuerpo secundario |
| Dylight 550 Conejo IgG | Thermofisher | 84541 | Anticuerpo secundario |
| SP8/SP5 o SPE confocal | Leica | Para la imagen | |
| Software de captura LAS Leica | Software de imagen | Leica | |
| Aviso (versión 7.0.0) | FEI Visualization Science Group | (Sección 2) | |
| Hypermesh parte del paquete Hyperworks (versión 10)Software de generación de | modelos Altair Engineering | ||
| Abaqus (versión 6.14) | SIMULIA | FE (Sección 5.7-5.8) |