Анализ методом конечных элементов — это часто используемый инструмент для исследования механических характеристик конструкций под нагрузкой. Здесь мы применяем его использование для моделирования биомеханики челюсти данио-рерио.
Методическая статья
Анализ методом конечных элементов — это часто используемый инструмент для исследования механических характеристик конструкций под нагрузкой. Здесь мы применяем его использование для моделирования биомеханики челюсти данио-рерио.
Скелетный морфогенез происходит через жестко регулируемое поведение клеток во время развития; многие типы клеток изменяют свое поведение в ответ на механическое напряжение. Скелетные суставы подвергаются динамической механической нагрузке. Анализ методом конечных элементов (МКЭ) — это вычислительный метод, часто используемый в инженерии, который позволяет предсказать, как материал или структура будут реагировать на механическое воздействие. Разделив всю систему (в данном случае скелет челюсти рыбки данио) на сетку более мелких «конечных элементов», МКЭ может быть использована для расчета механической реакции конструкции на внешние нагрузки. Результаты могут быть визуализированы различными способами, в том числе в виде «тепловой карты», показывающей положение максимальной и минимальной главных деформаций (положительная главная деформация указывает на растяжение, а отрицательная — на сжатие. Максимум и минимум относятся к наибольшему и наименьшему напряжению). Они могут быть использованы для определения того, какие области челюсти и, следовательно, какие клетки могут подвергаться особенно высоким растягивающим или сжимающим нагрузкам во время движения челюсти, и, следовательно, могут быть использованы для определения взаимосвязи между механическим напряжением и поведением клеток. В этом протоколе описываются шаги по созданию моделей конечных элементов на основе данных конфокального изображения опорно-двигательного аппарата на примере нижней челюсти рыбки данио-рерио. Протокол проводит читателя через ряд этапов: 1) окрашивание компонентов опорно-двигательного аппарата, 2) визуализация компонентов опорно-двигательного аппарата, 3) построение трехмерной (3D) поверхности, 4) создание сетки конечных элементов, 5) решение FEA и, наконец, 6) проверка результатов путем сравнения с реальными смещениями, наблюдаемыми в движениях челюсти рыбы.
Методом конечных элементов (FE) моделирование является инженерная техника , которая может вычислительно рассчитать и сопоставить величину и местоположение штаммов , действующих на конструкцию 1. Модель состоит из 3D-структуры, представленной сеткой "конечных элементов", и конечный результат анализа определяется целым рядом факторов, включая структуру и количество элементов в сетке, величины и расположения механического нагрузки и свойств материала. Существенные свойства описывают некоторые аспекты поведения материала при заданном типе нагрузки; модуль Юнга (Е) описывает эластичность материала в то время как коэффициент Пуассона описывает пропорциональное уменьшение ширины материала к его длине при растяжении образца. КЭ моделирование может быть использовано для вычисления различных переменных, включая перемещение, напряжение, давление и деформации, действующей на модели, принимая во внимание уникальные входные данные о структуре "; S форма, расположение и величина нагрузки и специфические свойства материала.
FE моделирование широко используется в технике 2 и все чаще для ортопедических 3 и палеонтологических приложений 4. В развитии биомеханические силы , как известно, выступают в качестве стимула во многих клетках , чтобы активировать клеточные ответы 5-8 и полезно для прогнозирования как относительные положения и величины механических стимулов в развивающихся систем органов, тем не менее, в настоящее время КЭ моделирование было мало используется для данио развития.
Оба хрящевой и костной было показано, что механочувствительных материалы. Например, сжатие в пробирке было обнаружено , чтобы активировать хондрогенное пути, в то время как напряжение было показано, что необходимо для образования костной ткани 9. FE анализ (FEA) было использовано для моделирования штаммов, действующие на биологических образцах, в том числе и те, которые действуют на скелетных элементов во время кости Foия 10. Другие разработки приложений включают его использование для прогнозирования форму сустава после того, как он подвергся воздействию теоретических биомеханических сил 11,12 и показать картину штаммов , присутствующих во время птенца коленного сустава морфогенеза 8.
Этот протокол направлен на обмен опытом создания 3-мерных поверхностей, сеток и конечно-элементных моделей из конфокальных изображений с целью понимания механики развивающихся тканей. Мы также покажем способы проверки моделей FE , хотя захватив реальное совместное перемещения информации в естественных условиях. В то время как мы используем данио челюсть как образец одни и те же методы могут быть использованы на любой малой биологической системы, для которой 3D информация о структуре костно-мышечной системы могут быть получены с помощью конфокальной или многофотонной томографии.
Все шаги в рамках протокола следуют Университет ухода за животными Бристоля и принципов социального обеспечения и тех, Министерство внутренних дел Великобритании.
1. Визуализация костно-мышечной анатомии
Примечание: Для того, чтобы визуализировать форму скелетных элементов, чтобы определить количество мышечной массы и определить точное расположение креплений мышц, immunostain (раздел 1.1) рыбы в соответствующем возрасте для скелетных миозина (который показывает мышцы) и тип II коллагена (для визуализации хрящ). В качестве альтернативы, визуализировать костно - мышечной анатомии , используя трансгенные флуоресцентные репортер линии , такие как коллаген a1 репортер Col2a1: mCherry 13,14 визуализировать хрящ и медленный тяжелой цепи миозина репортер smyhc: GFP 15 , чтобы визуализировать положение мышечных вложений (раздел 1.2).
Альтернативные линии, отмечающие хрящей и мышц может работать одинаково хорошо.
2. Создание 3D Surface
3. Расчет мышечных усилий, чтобы быть использованы в FE Model
4. Создание сетки
5. конечных элементов Modэль Строительство
6. Проверка Челюсти / Деформация смещения Расстояния
Иммуногистохимическое мышцы (рис 1А) и хрящевой ткани (рис 1б) или изображений трансгенных репортеров (фиг.1С) позволяет 3D структура челюсти , чтобы визуализировать, вместе с соответствующим мускулатуру. По визуализации при высоком разрешении можно было построить модель , которая как трехмерную форму челюсти (рисунок 2) и расположение и размещение грузов (рисунок 3). Использование в естественных условиях смещений видели при захвате видео высокой скорости (рисунок 4) , мы убедились , что диапазон движения в модели была в пределах реалистичного диапазона.
Модели FE После запуска может быть использован для отображения диапазона данных, таких как стресс (рисунок 5А), минимальной и максимальной главной деформации (фиг.5В - K). Эти результаты являются тремя тусклыми ensional поэтому модель можно увеличить , чтобы увидеть прекрасные образцы деталей (рис 5e, 5i) поворачивается , чтобы получить соответствующие взгляды (рис 5F, 5G, 5J, 5K) и в цифровом виде секционного (рис 5E ', 5E' ', 5i', 5i '') , чтобы показать , как паттерны напряжения, деформации или изменения давления по всей модели. Кроме того, можно извлечь количественные данные из модели (не показана). Подтвердив модели и используя самые точные свойства материалов, нагрузки и сетки форму модель FE может быть использован для изучения наилучшей оценки механической среды, испытываемых клеток во время этого окна развития. Результаты модели могут быть непосредственно по сравнению с изменениями в клеточном поведении и экспрессии генов 20.
Re 1 "SRC =" / файлы / ftp_upload / 54811 / 54811fig1.jpg "/>
Рисунок 1:. Характерные образы элементов опорно - двигательного аппарата данио нижней челюсти на 5 DPF Типичные конфокальной стеки нижней челюсти личинок 5dpf все показано с передним вверх (А) иммуногистохимическое A4.1025 , который окрашивает все скелетные миозин (B) иммуногистохимическое II типа коллагена , который помечает все хрящей (C) Стек из живой личинкой выражения трансгенные репортеров Col2a1: mCherry маркировка хрящ (красный) и smyhc: GFP медленный мышцы (зеленый). IA: intermandibularis передняя, PH: транспортир hyoideus, AM: аддуктор челюстями, HI: hyoideus хуже, HI: hyoideus превосходит, CH: sternohyoideus, MC: хрящ Меккелев, PQ: Palatoquadrate, CH:. Ceratohyal Пожалуйста , нажмите здесь , чтобы увидеть увеличенную версию этой фигуры.


Рисунок 3: Типичные сетки , показывающие ограничения и векторы силы Типичные сетки для 5 DPF личинкой для (A) и закрытия рта.(B) во рту открытие. Белые точки обозначают места , где ограничена модель и в которой размеры (например, х и у или х, у и г). Белые линии обозначают мышечные позиции, с вектором силы мышц, обозначенном белыми стрелками. Красный показывает хрящ и Желтый Interzone. Эта цифра была изменена из дополнительного материала , ранее опубликованного в Брант и др. 15. Пожалуйста , нажмите здесь , чтобы посмотреть увеличенную версию этой фигуры.

Рис . 4: тестирование чувствительности FE-модель имитации смещения челюсти в 5dpf данио для различных хрящей и межзоны модулей Юнга. Челюсть перемещение (открытого в закрытое в мкм) отмечен на челюсти; записанные с помощью ключа цвета. Каждая модель (А - Л) имеет различную комбинацию хряща (с = 1,1, 3,1 или 6,1 МПа) или межзоны (I = 0.25, 0.5, 0.75, или 1 МПа) свойствами. Горизонтальная черная стрелка указывает на значение смещения челюсти на кончик хряща Меккеля ( в лице вертикальной черной стрелкой). М и N кадры из видео , 5 DPF личинок показывая минимум, то есть, челюсти закрыты (М) и максимум, т.е. челюсть полностью открыт (N) с двумя наложенными (O) - белая линия на O представляет смещение (43 мкм). В этом случае относительная свойства хрящ 1.1 с Interzone 0,25 (А) лучший матч перемещения видели в живую рыбу (O). Панели AL этой фигуры были ранее опубликованы в Бранта и др. 15. Пожалуйста , нажмите здесь , чтобы увидеть увеличенное версион этой фигуры.

Рисунок 5: Типичные данные модели FE FE-имитационная модель всех мышц , применяемых в 5 денье личинка (А - С).. (A) Мизеса (EMaxmin) (B) Минимальная главной деформации (Е Мин. Р, μɛ) (C) Максимальная основная деформация (E Макс. P., μɛ). FE-имитационная модель максимальных и минимальных главных деформаций при открывании челюсти. (D - K): Максимальная главная деформация (. E Макс П., μɛ) в (D) вентральной челюсти и (E) вентральный совместный вид (Е) показывает расположение проксимальных-дистальных отделов через хряща Меккеля сустава и Interzone в (Е ') и (Е' '), соответственно. (F): боковаячелюсть вид. (G): боковой вид сустава. (H - K): Минимальная главной деформации (. E Min P, μɛ) в (H) зрения вентральной челюсти и (I) вентральный совместного просмотра. (I) показывает расположение проксимальной-дистальной участков через хрящ сустава Меккеля и межзоны в (I ') и (I' ') , соответственно , (J): вид сбоку в челюсть. (K): боковой вид сустава. Эта цифра была ранее опубликована в Бранта и др. 15. Пожалуйста , нажмите здесь , чтобы посмотреть увеличенную версию этой фигуры.
| Количество мышечных волокон | Область мышечного волокна (мкм 2) | Площадь группы мышц (мкм2) | Сила (N) | |
| 5 DPF intermandibularis передняя | 5 | 23,8 | 119 | 4.76e-6 |
| 5 DPF транспортир hyoideus | 6 | 23,8 | 142,8 | 5.71e-6 |
| 5 DPF аддуктор челюстями | 9 | 23,8 | 214,2 | 8.57e-6 |
Таблица 1: Мышцы Количественное. Вычисленной средней силы мышц для Intermandibularis Передней, Adductor челюстями и транспортир Hyoideus на 5 DPF с использованием 40 Nn / мкм 2 (стоимость в расчете на единицу площади , взятой из ссылки 17). (Lorga и др., 2011) (N = 3).
Конечно - элементных моделей были использованы для соотнесения областей скелетных элементов , которые находятся под напряжением с тех , которые подвергаются костеобразование 10, а также для отображения области , находящиеся под напряжением во время эндохондральной окостенение и совместного морфогенеза 8,12,21. Другие исследования также были в состоянии применить теоретические модели роста для репликации изменений во время совместной разработки 11,12. Здесь мы показываем протокол для построения моделей FE для относительно простой системы, данио челюсти 20. В отличие от альтернативных методов сбора необработанных изображений для моделей FE, таких как КТ 22, конфокальной микроскопии трансгенных линий или иммунологически данио позволяет нескольким ткани должны быть изучены. Он может, таким образом, обеспечивают прямую информацию о точках крепления мышц по отношению к хрящевой ткани. Среди позвоночных животных моделей данио особенно поддаются генетической и фармакологической манипуляции. Поколение моделей FE для даниокраниофациальные хрящ в настоящее время открывает возможность дальнейшего изучения взаимосвязи между биомеханики и генетики в совместном формообразования.
Есть целый ряд важных шагов в процессе создания модели FE; первый генерирует точный трехмерное представление системы. Это требует обработки изображений при достаточно высоком разрешении, чтобы четко определить границы. Обратите внимание, что даже с высокой разрешающей способностью, чтобы сделать хорошую поверхность, вероятно, придется сгладить некоторые регионы. Другим важным шагом является определение правильного размещения груза и правильных ограничений. Недостаточно ограничена модель будет не в состоянии решить и неправильное расположение нагрузок вызовет ненормальное движение.
Некоторые обработки необработанных данных (рисунок 2) необходимо в качестве поверхности , генерируется из необработанных данных , было бы трудно меш (Фигура 2В). Мы фильтруют данные с помощью гауссовского фильтра (рис 2C ) и мы провели некоторое ручное сглаживание кривых, чтобы произвести набор чистых контуров, которые могут быть преобразованы в 3D-поверхности. Слишком много сглаживания может произвести "топленое" поверхность, которая потеряла многие из его особенностей. Выбор правильного размера элемента является повторяющимся процессом, как выбор слишком мал размером элемент создает слишком большую сетку, которая является вычислительно интенсивным. Тем не менее, выбирая слишком большой размер элемента будет производить сетку материала, который не перепросматривать правильную форму структуры. Правильная сетка имела наименьший размер элемента, который захватил правильную форму челюсти и сходились на правильное решение, проверяется с помощью смещения челюсти. Он также может быть необходимо изменить свойства материала или расчетов нагрузки, чтобы лучше имитировать правильное смещение в разных возрастов и видов будут иметь существенно различные свойства.
Важно помнить, что всегда есть ограничения гипотетической модели иssumptions сделал для запуска моделей FE. При моделировании только один или небольшое число выборок, которые имеет решающее значение для обеспечения того, чтобы репрезентативная выборка выбирается так как, вероятно, будут небольшие вариации между людьми. Как только некоторые из элементов челюсти и мышцы были включены, модель представляет собой упрощенную версию данио краниофациальной системы опорно-двигательного аппарата. Таким образом, ограничения должны были быть установлены, чтобы учесть, где смоделированные элементы челюстной бы соединиться с остальной частью черепа и модель была искусственно ограничена в центре, чтобы зафиксировать его в "пространстве". Это искусственное ограничение не повлияло на интерпретации, проведенной из моделей как ceratohyal сам не был проанализирован. Включение большего количества краниофациальной структуры, особенно другие открытия мышцы челюсти , такие как sternohyals и подключенного к нему хряща 23, мог бы добавить к модели, но ограничения включают в себя способность больших моделей для работы в программном обеспечении конечных элементов.
s = "jove_content"> Другим ограничением является то, что мы не смоделирован вставки связки, хотя это может быть достигнуто путем введения пружин 8. Еще одно предположение, сделанное в этом случае было то, что модель будет вести себя линейно. Величины напряжений на моделях , были сопоставимы с теми , в опубликованных моделях и применены к клеткам в пробирке 10,24, штаммами быть ниже +3,500 и выше -5,000 μɛ кроме ограничений и прикрепления мышц точек. Таким образом, штаммы в соответствующих областях модели были признаны в пределах диапазона приемлемых для линейной модели. Хрящ не ведет себя полностью как линейный материал и ранее был смоделирован как пороупругой материал, который позволил анализировать поведение жидкости в модели 25. Распространение места прикрепления мышц среди кластера локальных узлов будет распределять пиковые силы и более точно отражают вставку мышц для определенных мышц.
лор "> Использование FE позволяет провести оценку деформаций и напряжений, действующих на конструкции. В качестве методики она часто используется во многих бионаучных дисциплин, включая ортопеда, палеонтологии и совсем недавно биологии развития. Здесь мы опишем, как построить FES для данио нижняя челюсть. в будущем эти модели могут быть расширены, чтобы смотреть на всю челюсть, включая вкус. Подобные методы могут быть использованы для моделирования позвоночника биомеханики в рыбе, которая до сих пор в основном были изучены с помощью кинематических средств.Авторы не имеют ничего раскрывать.
Некоторые данные на фиг 3-5 была перепечатана из J.Biomech, 48 (12), Брант и др., Моделирование методом конечных элементов предсказывает изменения в совместной форме и поведения клеток из - за потери мышечного напряжения в развитии челюсти, 3112-22. , 2015, с разрешения Elsevier 15.
LHB финансировался программой Wellcome Trust Dynamic Cell PhD; КАР финансировалось за счет гранта проекта MRC MR / L002566 / 1 (присужден и CLH Эйр) и CLH финансировалось за счет гранта Aruk 19479. Мы также хотели бы поблагодарить центр Wolfson биоимиджинга за советом к визуализации.
| Имя | Компания | Каталожный номер | Комментарии |
|---|---|---|---|
| Coll2 | Abcam | ab34712 | Антитела к коллагену типа II - окрашивает все хрящи |
| A4.1025 / MF20 | Исследования развития гибридома банка | A4.1025 | Скелетные антитела к мисоину - маркируют все скелетные мышцы |
| Малоплавкая агароза | Sigma | A9414-5G | Для монтажа данио-рерио |
| MS222 (этил-3-аминобензоатметасульфонат) | Sigma | E10521-10G | Для изготовления анестетика |
| Трипсин | Фишер | T/3760/48 | образец пермеабилизации |
| Dylight 488 Мышь IgG | Thermofisher | 35502 | Вторичное антитело |
| Dylight 550 Кролик IgG | Thermofisher | 84541 | Вторичное антитело |
| SP8/SP5 или конфокальный | SPE Leica | Для визуализации | |
| Программное обеспечение для захвата LAS Leica | , программное обеспечение для обработки изображений | Leica | |
| Aviso (версия 7.0.0) | Программное обеспечение для анализа 3D-изображений | FEI Visualization Science Group | |
| Hypermesh, входящая в состав пакета Hyperworks (версия 10) | Программное обеспечение для создания моделей FE от Altair Engineering | ||
| Abaqus (версия 6.14) | Программное обеспечение для анализа SIMULIA | FE (Раздел пункты 5.7-5.8) |