Waiting
Login processing...

Trial ends in Request Full Access Tell Your Colleague About Jove
JoVE Journal Bioengineering
Click here for the English version

Bioengineering

Характеристика сложных систем по постановке эксперимента подхода: Переходный Белок Expression табачными в качестве примера

Published: January 31, 2014 doi: 10.3791/51216
Automatic Translation
1, 1,2
1Institute for Molecular Biotechnology, RWTH Aachen University, 2Institute for Molecular Biology and Applied Ecology, Fraunhofer Gesellschaft

Summary

Опишем планирование эксперимента подхода, который может использоваться, чтобы определить и смоделировать влияние регуляторных элементов трансгенных, рост и развитие растений, параметров и условий инкубации на временной экспрессии моноклональных антител и репортерных белков в растениях.

Abstract

Растения обеспечивают множество преимуществ для производства биофармацевтических препаратов, включая низкие затраты, масштабируемость и безопасность. Кратковременная экспрессия предлагает дополнительное преимущество короткого развития и производства раз, но уровни экспрессии могут существенно различаться между партиями таким образом, послуживших основанием для регуляторных проблем в контексте надлежащей производственной практики. Мы использовали дизайн эксперименты (DOE) подхода для определения влияния основных факторов, таких как регуляторных элементов в конструкции экспрессии, рост растений и параметров развития, а также условия инкубации в течение выражения, по изменчивости экспрессии между партиями. Мы протестировали растения, выражающие модель анти-ВИЧ моноклональное антитело (2G12) и флуоресцентный маркер белок (DsRed). Обсуждается обоснование выбора определенных свойств модели и выявить его потенциальные ограничения. Общий подход может быть легко переведены в другие проблемы, поскольку принципы модели Авновь широко применимы: основанной на знаниях выбор параметров, снижение сложности путем разделения исходную задачу на более мелкие модули, программное обеспечение наведением установки оптимальных комбинаций эксперимента и поэтапного проектирования увеличения. Таким образом, методология не только полезно для характеристики экспрессии белка в растениях, но и для исследования других сложных систем, не имеющих механистический описание. Прогнозные уравнения, описывающие взаимосвязь между параметрами можно использовать для установления механистические модели для других сложных систем.

Introduction

Производство биофармацевтических белков в растениях является выгодным, поскольку растения являются недорогими расти, платформа может быть расширен только путем выращивания больше растений, и патогены человека не в состоянии воспроизвести 1,2. Переходные стратегии экспрессии на основе, например, на инфильтрацией листьев с Agrobacterium tumefaciens обеспечивает дополнительные преимущества, так как время между точкой доставки ДНК и доставки очищенного продукта снижается с нескольких лет до менее чем за 2 месяца 3. Переходный выражение используется также для функционального анализа, например, для тестирования генов для их способность к комплементации с потерей функции мутантов или для исследования белковых взаимодействий 4-6. Тем не менее, переходные уровни экспрессии как правило, показывают больший разброс от партии к партии, чем уровни экспрессии в трансгенных растениях 7-9. Это снижает вероятность того, что биофармацевтических производственных процессов на основе временной экспрессии Wiбудете утверждены в контексте надлежащей производственной практики (GMP), потому что воспроизводимость является атрибутом критическим качество и подлежит оценке риска 10. Такое изменение может также маскировать любые взаимодействия, что исследователи намерены провести расследование. Таким образом, мы намерены выявить основные факторы, которые влияют переходные уровни экспрессии в растениях и построить качественную количественную модель прогнозирования.

(OFAT) подход один фактор-на-а-время часто используется для характеристики воздействия (эффект) определенных параметров (факторов) на исход (ответ) эксперимента 11. Но это не является оптимальным, так как отдельные тесты (работает) в ходе расследования (эксперимент) будет приведена в соответствие, как жемчужины на нитке через потенциальный области, натянутого на факторы, которые проверены (дизайн помещения). Охват дизайн помещения и, следовательно, от степени информации, полученной из экспериментанизкий, как показано на рисунке 1A 12. Кроме того, взаимозависимости между различными факторами (коэффициент взаимодействия) может оставаться скрытым что приводит к плохим моделей и / или прогнозирования ложной оптимума, как показано на фиг.1В 13.

Недостатки, описанные выше, можно избежать с помощью планирование экспериментов (DOE) подход, в котором работает эксперимента разбросаны более равномерно в течение дизайна помещения, а это означает, что более чем один фактор варьируется между двумя отрезками 14. Существуют специализированные проекты для смесей, скрининг факторов (факторные планы) и количественная оценка воздействия фактора на ответах (методы поверхности отклика, RSM S) 15. Кроме того, RSMS может быть реализована как центральных композиционных конструкций, но также может быть достигнуто эффективно с помощью специализированного программного обеспечения, которое можно применить различные критерии для отбора трасс. Например, так называемый D-optimalitКритерий у выберет работает так, чтобы минимизировать ошибку в коэффициентах полученной модели, в то время как критерий IV-оптимальность выбирает работает, что достижения минимально дисперсию предсказания на протяжении всего дизайна помещения 15,16. РСМ мы описываем здесь позволяет точную количественную оценку временной экспрессии белка в растениях, но он может быть легко переведены в любой системе с участием нескольких (~ 5-8) числовые коэффициенты (например, температура, время, концентрация) и несколько (~ 2 - 4) категориальные факторы (например, промоутер, цвет), в котором механистический описание отсутствует или слишком сложны для моделирования.

Подход DoE возникла в сельскохозяйственных наук, но распространился на другие области, потому что это может быть передана любой ситуации, где это полезно уменьшить количество трасс, необходимых для получения достоверных данных и генерировать описательные модели для сложных процессов. Это в свою очередь привело к включению DoE в «Руководства дляПромышленность, Q8 (R2) Фармацевтическая разработка ", опубликованном Международной конференции по гармонизации технических требований к регистрации лекарственных препаратов для человека (ICH) 17. Министерство энергетики теперь широко используется в научных исследованиях и промышленности 18. Тем не менее, следует соблюдать осторожность во время планирование и проведение эксперимента, потому что выбора неправильного степень полинома для модели множественного линейного регрессионного (базовая модель) можно ввести необходимость в дополнительных трасс правильно моделировать все эффекты фактора. Кроме того, повреждены или отсутствуют данные генерировать неправильные модели и недостатки прогнозы, и даже может предотвратить любую модель здания попытку, как описано в секциях протокола и дискуссионных 18. В разделе протокола, мы первоначально установлен наиболее важные этапы планирования для эксперимента РСМ-основе, а затем объяснить дизайн, основанный на DoE Программное обеспечение DesignExpert v8.1. Однако подобные конструкции могут быть построены с другим программным обеспечением includiнг JMP, Modde и STATISTICA. Экспериментальные процедуры следуют инструкции для анализа и оценки данных.

Рисунок 1
Рисунок 1. Сравнение OFAT и НОО. А. Последовательное изменение одного фактора, в то время (OFAT) в эксперименте (черный, красный и синий кружки) достигает малый охват дизайн помещения (заштрихованные области). В противоположность этому, изменение более чем один фактор за раз, используя конструкцию экспериментов (DOE) Стратегии (зеленые кружки) усиливает охват и, следовательно, точность результирующих моделей. B. Предвзятым дизайн пространства охват означает, что OFAT эксперименты (черные кружки) могут и не выявить оптимальные эксплуатационные регионы (красный) и предсказать суб-оптимальные решения (большой черный круг), в то время как Министерство энергетики стратегическиэс (черные звезды), скорее всего, определить предпочтительные условия (большой черный Star).

Subscription Required. Please recommend JoVE to your librarian.

Protocol

1. Планирование стратегии DoE

  1. Определите соответствующие факторы и ответы для включения в проект.
    1. Определите один или несколько ответов для измерения. Здесь были использованы уровни экспрессии 2G12 и DsRed (мкг / мл), в том числе минимальном допустимом разницы рассматривается как соответствующие (10 и 20 мкг / мл, соответственно) и приблизительное значение для предполагаемого стандартного отклонения системы (4 и 8 мкг / мл, соответственно) на основе предыдущих экспериментов.
    2. Используйте имеющуюся литературу, данные из предыдущих экспериментов или специализированных конструкций скрининга (например, факторного плана, см. введение), чтобы выбрать значимые факторы, воздействие которых на ответах будет количественно (табл. 1) 7,8,19,20.
    3. Выделяют типы фактор (цифровые или категориального) и выбрать диапазон, в рамках которой цифровые факторы будут изменяться в ходе расследования DoE (табл. 1).
    4. Определить численныйс факторами, для которых непрерывное изменение трудно реализовать.
      1. Не следует использовать непрерывное изменение факторы, которые не могут быть отрегулированы точно определенном уровне. Например, температура инкубации обычно можно контролировать только в пределах ± 2 ° C, поэтому непрерывное изменение, используя значения, такие как 27,2 ° С, 25,9 ° С, и 29,3 ° С должна следует избегать.
      2. Вместо этого, выберите ряд дискретных уровней этих факторов (табл. 1). Число уровней должно соответствовать ожидаемой базовой модели (шаг 1.2). Это важно только для оптимальных конструкций, потому центральных композиционных всегда иметь дискретные уровни факторов.
    5. Связать ли категорическое фактором является номинальной, то есть никакой подразумеваемый порядок (например, различные производители), или порядковый и, таким образом дискретную числовой параметр (например, различные листья на заводе).
    6. Выберите уровни для обоих типов категориальных факторов. </ Li>
  2. Выберите полезную базовую модель.
    1. На основании предварительных экспериментов и литературы, предвидеть связь между каждого фактора и реакции, а также фактор взаимодействия и реакции. Например, линейный (чем больше, тем лучше), квадратичной (один оптимальный или нелинейное увеличение / уменьшение) или кубический (перекос оптимальный).
    2. Убедитесь, что количество уровней для дискретных числовых факторов п + 1, с п будучи многочлен степень связи между фактором и ответ. Например, если температура, как ожидается, имеют квадратичную действие на ответ, п = 2 и, следовательно, по крайней мере, три уровня температуры должны быть исследованы, чтобы достичь подгонку к квадратичного эффекта.
  3. Определить долю пространства дизайна (FDS), для которых дисперсия предсказание должно быть ниже определенного порога (альфа уровень). FDS должна быть> 0,95 (охватывают более 95% от дизайна помещения) для достижения модели, которые дают надежные прогнозы на протяжениивесь дизайн помещения 21-23.
    Примечание: порог 0,05 соответствует 5% уровне значимости (тип я ошибка) и является типичным значение. Увеличение этого значения приведет к сокращению числа экспериментов, необходимых для стратегии DoE но одновременно повысит вероятность того, что значительные последствия будут пропустили и что оценки их влияния на ответ будет неправильным.

Таблица 1
Таблица 1. Факторы, влияющие на временную экспрессию белка в табака, включая вариации колеблется во Доу. Факторы выделены жирным шрифтом были только включены в дизайн для экспериментов, описанных в разделе "A описательной модели для накопления DsRed во временной экспрессии с использованием различных промотор / 5'UTRs", тогда как факторы, выделенные курсивом были включены только в дизайне для "Оптимизация incubatioN условия и схемы урожай для производства моноклональных антител в растениях с использованием временной экспрессии ".

Рисунок 2
Рисунок 2. Министерство энергетики процесс планирования. Факторы со значительным воздействием на ответ исследуемого выбираются на основе имеющихся данных. Тогда атрибуты фактор (например, цифровые), диапазоны и уровни присваиваются. Предыдущий знания и эксперименты используются для определения подходящего базовой модели. Требования предсказательной силой определяются на основе применения / целью окончательной модели. Собранные данные затем могут быть переданы в соответствующее программное обеспечение Доу.

2. Настройка RSM в DesignExpert

  1. Начните DesignExpert и выберите пункт "New Design". В окне "поверхности отклика", выберите &# 34; Оптимальное "и введите количество числовых и категориальных факторов, выбранных в разделе 1.1).
  2. Введите имена фактор, блоки, типы, подтипы, количество уровней / границ уровня и значения для уровней для всех факторов в соответствующих полях.
  3. Продолжить на следующую страницу и под "Поиск" Варианты выбора "Бест", а также желаемый "Об оптимальности" критерий, здесь "D-оптимальный".
  4. В меню "Правка модель ...", выбрать модель, которая содержит факторы и взаимодействия ожидаемые в разделе 1.2). Если вы не уверены, выберите полную модель, которая охватывает все факторы и взаимодействия для определенной степени полинома, здесь "Квадратичные". Обратите внимание, что выбирая полную модель (содержащий все взаимодействия) может увеличить количество экспериментов, необходимых для DoE.
  5. Выберите число "блоков". Здесь, один блок был использован, так как все растения из той же партии, все бактерии культивировали одновременно иВсе инъекции были проведены в тот же день одним и тем же оператором. Используйте более одного блока, если более одного партия растений используется или несколько операторов обработки инъекции.
  6. Баланс дизайн
    1. Включить "Force категорический баланс" для того, чтобы распределить работает эксперимента поровну между категориальных уровней факторов, например, различные промоутеры, которые проверяются.
      Примечание: это может уменьшить оптимальность проектирования до такой степени, что DesignExpert сообщит как% значения в зависимости от точек, выбранных по алгоритму оптимальности.
    2. Пересчитать дизайн несколько раз с помощью реализован случайный алгоритм семян, чтобы минимизировать снижение оптимальности.
      Примечание: потери в оптимальности может быть в диапазоне от ~ 3-40% с использованием тех же входных данных, но изменение количества повторов может быть использован для поддержания баланса категориальный.
  7. Программное обеспечение будет предложить значение для числа «точек типовых" на основе базовой модели, Здесь 70 трасс. Отрегулируйте количество показов для "повторах" и "Для оценки отсутствие подгонки", чтобы убедиться, что он> 5% из числа «точек типовых" для каждого, здесь 10 трасс каждый.
  8. Перейдите к следующей странице выберите количество ответов и введите их имена и единиц. Дополнительные ответы также могут быть добавлены в любое время при оценке данных без негативного воздействия на конструкции.
  9. Продолжить, чтобы начать алгоритм, расчета уровней фактором для прогонов DoE. Если выбранный базовая модель не соответствует количество уровней в течение определенного фактора, выводится уведомление и расчет не начнется. В этом случае либо:
    1. Увеличение количества уровней для этого фактора, или
    2. Снимите условия в базовой модели, что не может быть рассчитана исходя из текущего количества уровней.
    3. Например, если есть два уровня для фактора «Время инкубации" т (т.е. 2 г и 5 г) в тон в настоящее время конструкция и квадратичной базовая модель в том числе термина т 2 был выбран, либо добавить третий уровень для т (т.е. 8 г) или удалите фактор T 2 из модели.
  10. Сохраните таблицу, содержащую оптимальные сочетания факторов, которые отображаются сразу завершения расчета.
  11. В "Дизайн" узла выберите "оценка" суб-узел и перейдите на вкладку "Графики".
    1. В "инструмент Графики" выбрать "FDS" и в поле "FDS Graph" выберите "Pred" в качестве типа ошибки. Затем введите минимального обнаруживаемого разницу, а также приближенное значение для предполагаемого стандартного отклонения системы, определенной в разделе 1.1.1) в качестве значений "г" и "е" соответственно (здесь 20 и 8 мкг / мл для DsRed). Также введите уровень альфа (приемлемым% отсутствует существенное влияние), подходящую для применения, как правило, 0,05 (5%).
    2. Если это не так (как найти здесь, FDS был 1%, как показано на рисунке 3А), вернуться к «Дизайн» узла и на панели задач, выберите "Средства разработки", затем "Augment дизайн ..." и выберите " Дополните ".
    3. Выберите те же "Поиск" и "Об оптимальности" критерии, как раньше. Также выбрал тот же "Изменить модель" и "форс-категориальный баланс" Настройки. Если дизайн увеличение выполняется до любого эксперимента (как это было сделано здесь), а затем изменить "Положите работает в блок" в положение "Блок 1".
    4. В разделе "прогонов" введите количество дополнительных «точек модели", сохраняя, по крайней мере 5% "повторах" и "Для оценки отсутствиеиз вписаться "в общем дизайне. Здесь добавлены 100 дополнительных Опыты и никаких" Реплицирует "и" Для оценки отсутствие подгонки "не были включены.
    5. После расчета завершена, пересмотреть график FDS, как описано выше. Если FDS все еще неудовлетворительное, повторите дизайн увеличение как описано выше. Здесь FDS был на 100% после увеличения и таким образом не далее увеличение не требовалось (рис. 3В).

Рисунок 3
Рисунок 3. Сравнение FDS участков. . Министерство энергетики, состоящий из 90 серий производит недостаточное FDS всего в 1% для стандартной ошибки прогноза, используя квадратное базовой модели в сочетании со значениями для минимального обнаруживаемого разницы (20 мкг / мл) и оцеповязана стандартное отклонение системы (8 мкг / мл). B. Увеличение НОО в общей сложности 210 трасс достигается 100% FDS и плоскую кривую с указанием равномерное точность модели во всем дизайн помещения.

3. Клонирование и анализ экспрессии Кассеты

  1. Выращивают кишечной палочки в 5 мл LB среды (10 г / л триптона, 5 г / л дрожжевого экстракта, 170 мМ NaCl, 50 мг / л ампициллина; для LB-агаром включают 15 г / л агара) при 37 ° С в течение 6 -8 час или на ночь на орбитальном шейкере при 160 оборотах в минуту. ВНИМАНИЕ: ампициллин является вредным веществом. Привить среду LB с или 50 мкл жидкой культуры или передаче клеток из колоний, выросших на чашках с использованием пипетки.
  2. Для очистки плазмидной ДНК ПСО и других производных PPAM (GenBank AY027531), из Е. Штамм E.coli К12 DH5a, следуйте инструкциям в комплекте руководства очистки ДНК 24. ВНИМАНИЕ: очищение кон содержит вредных химических веществ (см. выше пособия для деталей). Последовательность плазмиду описано на фиг.4 и Покупатель и соавт. 20.
  3. Определить концентрацию ДНК в очищенном элюата путем измерения оптической плотности образца 2 мкл при 260 нм в NanoDrop устройства.
  4. Подтверждение личности очищенной ДНК плазмиды расщеплением рестриктаз (ВИЭ).
    1. Выберите УЭ в соответствии с последовательностью плазмиды, так что уникальные и различимые узоры фрагмент генерируются. Следуйте рекомендациям производителя по условиям пищеварения, таких как реакционного объема, времени, температуры и концентрации стабилизатора бычьего сывороточного альбумина. В зависимости от чувствительности устройство анализа, используют 50-500 нг плазмидной ДНК в пищеварении.
    2. Подготовка 0,8-2,0% гели для разделения фрагментов ДНК в агарозном кипячением в Трис-борат-ЭДТА (КЭ) буфера (90 мМ Трис, 90 мМ бората, 2 мМ ЭДТА, рН 8,0). Тон больше ожидаемых фрагментов, тем меньше агарозы должны быть использованы в геле.
    3. Добавьте 5 мкл пять-кратного буфера для образцов (5x SABU, 0,1% (вес / объем) бромфеноловый синий, 0.1% (вес / объем) ксиленцианол, 10% (вес / объем) глицерина в КЭ) в 50 мкл RE-обработанный образец ДНК и отделить фрагменты электрофорезом в агарозном геле при 100 В в течение ~ 40 мин или до четкого разделения фрагментов достигается. На каждый гель, включают в себя полосу, содержащую размер ДНК лестница маркеры для сравнения, например, 2-3 мкл 1-кб лестнице.
  5. Замените омега 5'UTR в PSO с одним из трех других 5'UTRs.
    1. Отпустите омега 5'UTR последовательность из ~ 4 мкг очищенной ДНК PSO обработкой 20 единиц Eco RI-HF и сержантского Я-HF ВИЭ в NEBuffer 4 при 37 ° С в течение ~ 60 мин. Тогда отделить фрагменты, как описано в разделах 3.4.2 и 3.4.3.
    2. Изолировать больший фрагмент (PS, "становым хребтом") из агарозного геля с помощью гель дополнительныйводств Набор по изготовитель Руководство 25 и определения концентрации очищенной ДНК, как описано в разделе 3.3). ВНИМАНИЕ: Комплект для слива гель содержит вредных химических веществ, см. в руководстве по эксплуатации.
    3. Изолировать CHS, CHS-ЛПХ и TL 5'UTRs из подходящих векторов доноров обработкой ~ 10 мкг каждого вектора с 20 единиц Eco RI-HF и Nco Я-HF ВИЭ в NEBuffer 4 при 37 ° С в течение ~ 60 мин. Тогда отделить фрагменты, как описано в разделах 3.4.2 и 3.4.3 и очистить меньшую 5'UTR-содержащий фрагмент, как описано в разделе 3.5.2.
    4. Лигируют изолированные очищенные 5'UTRs в отдельных аликвотах в изолированной очищенной линейная пСм вектора при 25 ° С в течение 5 мин в соответствии с рекомендациями 26 изготовителя. Использование ~ 50 нг векторной ДНК и три-кратный молярный избыток 5'UTR ДНК в общем объеме 20 мкл.
  6. Трансформация E палочка клетки с гое рекомбинантные плазмиды 26,27.
    1. Добавить ~ 10 нг (~ 4-5 мкл) лигирования смеси (см. раздел 3.5.4) до 50 мкл RbCl компетентным Е. палочка и аккуратно перемешать, а затем инкубировать в течение 30-60 мин на льду. Тепловой удар в течение 1,5 мин при 42 ° С и охладить на льду в течение 5-30 мин.
    2. Добавить 950 мкл антибиотической свободной LB среды и инкубируют в течение 1 ч при 37 ° С и 160 оборотах в минуту. Для отбора трансформантов, распространение 50 и 100 мкл культуры на чашках с агаром LB, содержащих ампициллин и инкубировать в течение ~ 16-20 часов при 37 ° С.
    3. Привить 5-10 отдельных колоний, представляющих каждый промоутер-5'UTR перевязки, описанный в разделе 3.1 в 5-мл аликвотами LB среде, содержащей ампициллин. Тогда очистить ДНК плазмиды и его идентификации, как описано в разделах 3.2 до 3,4.
  7. Замените 35SS промотор с промоутером н.у.к в каждой из четырех 5'UTR плазмид использованием ВИЭ возрастанию I и Eco RI в NEBuffER 4 при 37 ° С в течение 1 ч, как описано в разделах 3.5 и 3.6.
  8. Представьте каждой из восьми полученных плазмид в А. tumefaciens штамм GV3101: pMP90RK путем электропорации 28.
    1. Добавить ~ 500 нг очищенной плазмидной ДНК в 50 мкл компетентного А. tumefaciens клетки на льду. Аккуратно перемешайте и передать в prechilled 0,2 см кювету для электропорации. Убедитесь, что смесь на дне кюветы и не содержит пузырей.
    2. Импульсные клетки на 2,5 кВ в течение 5 мс и подтверждающие интенсивность и продолжительность импульса.
    3. Не допускать повышения концентрации соли в образце ДНК или неправильной подготовки компетентного А. tumefaciens клетки, поскольку это может привести к высоких ионных токов, вызывающих мгновенное испарение клеточной суспензии, что значительно снижает эффективность трансформации.
    4. Добавить 950 мкл YEB среде без антибиотиков (выдержка 5 г / л говядина, 1 г / л дрожжевого экстракта, 5 г / л пептоне, 5 г / л сахарозы, 2 мМ MgSO 4, рН 7,0), осторожно перемешать и передавать непосредственно в стерильный 1,5 мл реакционной трубки. Инкубировать в течение 2-4 ч при 26-28 ° С и 160 оборотах в минуту.
    5. Распространение 1-2 мкл на YEB агаровых чашках, содержащих антибиотики (50 мг / л   карбенициллина, 25 мг / л канамицина, 25 мг / л рифампицина) для отбора трансформантов. ВНИМАНИЕ: рифампицин является токсичным веществом.
    6. Инокулировать три 5 мл аликвоты YEB среде, содержащей антибиотики с клетками из отдельных колоний, представляющих каждый промотор-5'-UTR лигирование и инкубировать в течение 48-72 ч при 26-28 ° С и 160 оборотах в минуту.
  9. Подтвердите успех А. tumefaciens трансформации.
    1. Передача 2 мкл аликвоты каждого из раздела 3.8.6 отделить 48 мкл аликвоты ПЦР Mastermix (2 мкл каждого из 10 мкМ праймера складе (400 нМ конечная концентрация каждого праймера),1 мкл 10 мМ смеси дНТФ (200 мкМ конечную концентрацию каждого дезоксинуклеозид-трифосфата), 5 мкл 10-кратного Expand High Fidelity буфера с 15 мМ MgCl 2, 0,75 мкл Expand High Fidelity смеси ферментов, и 39,25 мкл стерильной дистиллированной воды).
    2. Amplify кассете экспрессии каждой плазмиды с использованием соответствующих праймеров (Fwd: 5'-CCT CAG GAA GAG CAA TAC-3 ', связывающий 1026 нуклеотидов в обратном направлении от промотора; REV: 5'-AAG CCA CGA GTA CAC AAC-3', связывание в полиаденилирования сайте 35S) при соответствующих условиях ПЦР. Здесь, 94 ° C была использована для первоначальной денатурации с последующими 30 циклами 94 ° С денатурации в течение 15 сек, 51 ° C отжиг в течение 30 сек и 72 ° С для удлинения 120 сек, а конечная стадия удлинения при 72 ° С в течение 8 минимум
    3. Определите размер продуктов ПЦР с помощью электрофореза в агарозном геле, как описано в разделе 3.4.3.
  10. Подготовьте А. tumefaciens запасы глицерина посмешивание 500 мкл 50% (об / об) стерильный глицерин с 500 мкл А. tumefaciens культур из раздела 3.8.6, для которых преобразование было подтверждено (раздел 3.9.3). не хранить запасы глицерине при -80 ° С до дальнейшего использования.
  11. Вычислить энергию складные различных мРНК с использованием RNAfold веб-сервер 29 и включают в себя нуклеотидную последовательность от сайтом начала транскрипции до первого 50 п.н. кодирующей области.
    1. В разделе "Fold алгоритмы и основные параметры" выберите "минимум свободной энергии (МКЭ) и функцию секционирования" и "избежать изолированных пар оснований" варианты.
    2. В разделе "Дополнительные параметры складные" выберите "болтались энергии с обеих сторон спирали в любом случае", "параметры РНК (Тернер модели, 2004)" и изменить "масштабировать энергетических параметров для данной температуре (C)" до 25 ° С.
    3. В разделе "Опции вывода" выбрать все варианты.

Рисунок 4
Рисунок 4. Промотором и 5 'UTR вариантов. Кассеты экспрессии были получены путем ступенчатой ​​обмена 5'-UTR, в результате чего четыре комбинации с промотором CaMV 35SS с последующим замещением этого промотора с последовательностью NOS получают четыре дополнительные варианты и в общей сложности восемь различных комбинаций промотор / 5'UTR.

4. Растениеводство

  1. Подготовьте 0,1%-ный раствор удобрения Ferty2 Мега в деионизированной воде при рН 5,9.
  2. Подготовка 10 х 10 х 8 см Rockwool блоков обширной промывки деионизованной водой для удаления остаточных химических веществ и, наконец, равновесие с удобрением.
  3. Семя растения табака путем размещения 1-2 семена табака на каждого блока минеральной ваты последующим кратким заподлицо с удобрением будьте осторожны, избегайте мытья семена далеко.
  4. Прорастают и культивировать растения табака в течение 42 дней в теплице при 25/22 ° C температура днем / ночью, относительной влажности 70% и 16-часовой фотопериода (180 ммоль с -1 м -2; λ = 400-700 нм) . В течение этого фотопериода, орошения с удобрением в течение 15 мин каждый час в гидропоники культуры.

5. Переходные Экспрессия белка

  1. Подготовьте А. tumefaciens для инъекций в листьях.
    1. Инокулировать 5-50 мл YEB среде, содержащей антибиотики с 1% А. tumefaciens крио-исходной культуры и инкубировать при 27 ° С, пока ОП <суб> 600 нм достигает 5,0 (~ 48-72 ч в зависимости от объема и типа судна).
    2. Развести А. tumefaciens культура водой и 2 раза инфильтрации среды (4,3 г / л Мурасиге и Скуга солей (рН 5,6), 5 г / л сахарозы, 1,8 г / л глюкозы, 100 мМ ацетосирингона) в соответствии с ОП 600 нм, необходимое для инъекций. Подтвердите OD 600 нм непосредственно перед инъекцией. Примечание: OD 600 нм 1.0 соответствует ~ 1,43 ± 0,12 х 10 9 колониеобразующих единиц на мл.
  2. Введите А. tumefaciens суспензии в листьях.
    1. Выберите и маркировать noncotyledon листья и позиции по этому вопросу рассматриваться (например, в соответствии со стратегией DOE).
    2. Не вводят различные A. tumefaciens решения в том же межреберной области. Вместо этого используйте противоположные разделы по каждой стороне среднего вены оси. Если более чем два различных решения должны быть введены в том же положении использовать additionaл завод.
    3. Тщательно встряхните А. tumefaciens решение ресуспендируйте любые клетки, которые, возможно, собралась с подготовки разбавление и аспирации в 1-мл шприц.
    4. Аккуратно поцарапать эпидермис на место, предназначенное для инъекций с кончика пипетки или подобного содействия приток в А. tumefaciens решение. Избегайте разрыва листовой пластинки, делая так.
    5. Удерживая шприц перпендикулярно к листовой пластинки касаясь барр к межреберной области, подлежащей обработке, и осторожно нажмите на выход (без иглой) на нижней стороне листьев. Нажмите верхнюю сторону листа осторожно в то же время предотвратить листовой пластинки от сдвига или разрыва.
    6. Осторожно надавите на поршень шприца. А. Раствор tumefaciens войдет в межклеточные пространства внутри листовой пластинки, как показано на обработанных участках, возникающих темно-зеленый и влажной. Повторяйте эту процедуру в нескольких местах, пока вся intercostaл поле проникли с А. tumefaciens. Тогда переходите к следующему межреберной области.
    7. Убедитесь, что шприц остается перпендикулярно к листу. Наклон шприц может привести к бактериальной суспензии на бить струей под высоким давлением.
    8. Если разные А. tumefaciens используют растворы (например, для тестирования различных промоторов), удалить любой избыток раствора, остающегося на нижней стороне листа от первой инъекции с использованием бумажное полотенце или подобный перед нанесением следующего раствора.
  3. После инфильтрации инкубации растений и отбора проб.
    1. Подготовьте фитотроне для обработанных растений ~ 24 часа в сутки перед инъекцией, чтобы позволить температуры и влажности уравновешивание на уровнях, необходимых для DoE.
    2. После инъекции, передавать растения в фитотроне и разместить их в лотки, достаточно больших для орошения водой до конца инкубационного периода, определяемой Министерством энергетики.
    3. Создание 16-часовом фотoperiod с использованием шести Osram холодные белые 36 Вт Люминесцентные лампы на 0,7 м 2 (75 ммоль с -1 м -2; λ = 400-700 нм). Предотвратить растения от затенения друг друга, ограничивая количество растений до шести за 0,7 м 2.
    4. Перед выборки, убедитесь, что правильно межреберная поле была выбрана путем сравнения этикетку добавленную в п. 5.2.1 и план Доу.
    5. Используйте пробки буром для удаления 4-5 дисков листьев из обработанных межреберных полей на должности и времени, указанные в DoE. Не снимайте весь лист с завода во время отбора проб. Стабилизировать лист с ручным бумажным полотенцем, удаляя диски, чтобы предотвратить разрыв.
    6. Определяют массу каждого образца и поместить его в 1,5 мл пластиковую пробирку реакции с именем образца и массы. Хранить образцы при -20 ° С или -80 ° С до белка количественного определения. Этот процесс может быть приостановлен на этой стадии в течение нескольких месяцев в зависимости от стабильности образцаи температура хранения.

6. Белок Количественное

  1. Извлечение белков из образцов листовых дисков.
    1. Добавить 3 мл экстракционного буфера (50 мМ фосфата натрия, 500 мМ хлорида натрия, рН 8,0) на мг массы образца и измельчить дисков листьев в реакционной трубке с помощью электрического пестиком до тех пор, никаких больших фрагментов остаются. Избегайте перегрева образца.
    2. Удалить диспергированных твердых частиц центрифугированием образца в два раза при 16000 х г в течение 20 мин при 4 ° С. Передача супернатант в чистую 1,5 мл реакционной трубки после каждого шага, не нарушая гранул.
    3. После центрифугирования, процесс может быть приостановлена ​​путем замораживания растительных экстрактов при температуре -20 ° С или -80 ° С в течение нескольких месяцев в зависимости от стабильности образца и температуры хранения. Убедитесь, что цикл при замораживании-оттаивании не влияет на концентрацию белка-мишени (ей).
  2. Измерьте флуоресценции DsRed.
    1. Подготовьте три технических повторов каждого образца в черные 96-луночные планшеты половине площади (50 мкл экстракта на лунку). Избегайте образования пузырьков во время пипетки.
    2. Использование набор из шести разведений (0, 25, 75, 125, 175 и 225 мкг / мл) стандартного DsRed в 96-луночный планшет для генерации эталонной гильзы кривой. Подготовка разведений в PBS и хранить их при 4 ° С для использования в течение 3 месяцев.
    3. Измерьте флуоресценции дважды, последовательно, в ридере 96-а, снабженной возбуждения 530/25 нм и 590/35 фильтров выбросов нм.
    4. Для каждого образца, в среднем флуоресценцию над ними гласит и три технических повторов и вычтите значение записанного для пустой управления, содержащий 0 мкг / мл DsRed. Кроме того вычесть эту величину из считывает из стандартных разведений и использовать эти пустые коррекцией значения для линейной регрессии получая калибровочную кривую (через начало координат ординат).
    5. Использование наклон эталонной кривой для преобразования Fluorescence, измеренна в образцах в DsRed концентрациях. При необходимости разбавить образцы так, что значение находится в интервале концентраций стандартов и учитывать этот фактор разбавления в последующих расчетах.
  3. Определите концентрацию 2G12.
    1. Подготовка поверхностного плазмонного резонанса (SPR) устройство для измерения концентрации антител путем связывания белок А к активированной поверхности одной клетки потока. Используйте другую ячейку потока в качестве эталона путем инактивации поверхность без сцепления белка 30,31.
    2. Развести растительные экстракты 1:20 в SPR рабочем буфере (10 мМ HEPES рН 7,4, 3 мМ ЭДТА, 150 мМ NaCl, 0,05% об / об Tween-20) и измерить единиц отклика (RU) связывания антитела с белком А в три технические повтора каждого образца в конце инъекции 90 мкл (180 сек при 30 мкл / мин). Вычтите RU оцениваются по клетке ссылки потока (без белка A) от РУ измеренного в экспериментальной этой ячейки (Белок поверхность).
    3. Измерение стандарт 585 нг / мл через каждые 2G12 10-15 образцов и вычесть величину, измеренную в клетках ссылки потока, как описано выше. Использование среднюю RU этих стандартов для вычисления линейной калибровочной кривой через начало координат ординат.
    4. Рассчитать концентрацию 2G12 в образцах на основе разведения 1:20 и наклон эталонной кривой.
    5. Проверьте сильный неспецифическое связывание с опорной ячейки, которые могут привести к повреждению измерение. Также проверьте стандарты 2G12 сохранить приблизительно постоянные значения (<5% вариации) на протяжении анализа, как более сильному изменению отражает старение белка поверхности.

7. Анализ и оценка данных

  1. Вручную анализировать ответы, наблюдаемые в экспериментах экспрессии для (I) экстремальных значений (высоких или низких), указывающий на ошибки измерения; (II) необычных результатов, например комбинаций факторов, которые генныхоценить неожиданную реакцию, указывающие обмена данными, и (III) отсутствующих значений.
    Примечание: использовать эти меры для предотвращения в построение моделей на основе ошибочных данных.
  2. Передача данных анализируемые реагирования (здесь Концентрации белка) в "Дизайн" узла DesignExpert. Убедитесь, что ответные данные правильно назначены для соответствующих установок факторов. Существует обширный раздел помощь доступна в DesignExpert программного обеспечения, которая охватывает аспекты рассматриваются ниже.
  3. В «Анализ» узла, выберите ответ анализируемый и первоначально выберите "нет" на вкладке трансформации.
    Примечание: преобразование рекомендуется для мин / макс соотношениях ответа больше, чем 10. Самое полезное преобразование может быть получено из Box-Cox-участка на вкладке "Диагностика" в разделе "Диагностика" в "инструмент Диагностика", описанной ниже (раздел 7.8). Вход 10 преобразование часто необходимо.
    4. (например взаимодействий двухфакторных, квадратичных эффектов). Программное обеспечение будет предложить начальную модель, основанную на его значении.
  4. На вкладке "типового", начальная модель выбирается заранее на основе «Установить резюме" результатов. Используйте автоматизированном режиме редактировать эту модель:
    1. Выберите "порядок процесс", который является одним заказа выше предложенной модели, например, если предложенная модель является "2FI" (два фактором взаимодействия) выберите "квадратичной".
    2. Выберите «назад» в поле «Выбор», которая будет многократно удалить незначительное термины из модели после перейти к вкладке "ANOVA» на основе значения «Альфа из», который изначально должен быть 0,100.
    3. На вопрос, с помощью программного обеспечения, всегда автоматически корректировать модель иерархии, потому что это necessarу генерировать надежные модели 32,33.
  5. На вкладке "ANOVA", исследовать предложенную модель и включенные факторы. При необходимости вручную удалить любые факторы с р-значениями выше установленного порога (здесь 0,05, что соответствует 5% уровне значимости) или те, которые вряд ли на основе механистического рассмотрения перейти обратно на вкладке "Модель", изменяя "Выбор" "Ручной" и устранения соответствующих коэффициентов из модели.
  6. Вернуться во вкладку "ANOVA", оценить обновленный модели в терминах р-значения "Модели" (низкое значение является желательным, поскольку это свидетельствует о значимости) и "Отсутствие подгонки" (высокое значение желательно потому что это указывает на отсутствие значения), а также значения для "R-квадрат", "Скорректированный R-квадрат" и "Прогнозируемая R-квадрат" (значения> 0,80 желательны для всех трех значений).
    1. Сравните эти значения FOг различных моделей, включая / исключая факторов на нескольких уровнях значимости.
      Примечание: это может быть полезно, чтобы дублировать столбец реагирования в «Дизайн» узла и выполнять каждый анализ индивидуально или экспортировать всю таблицу "ANOVA" в другую программу, например, электронную таблицу.
  7. Перейдите к закладке "Диагностика", чтобы подтвердить качество модели и выявления потенциальных выбросов в наборе данных, которые имеют сильное влияние на модели путем изучения все вкладки в «Диагностика Tool" ("Диагностика" и раздел "Влияние") .
    1. В разделе "Диагностика" в "инструмент Диагностика", выполните следующие действия:
      1. Убедитесь, что точки беспорядочно в пределах в "Разности и предсказанных" диаграммы. Здесь шеврона образное распределение указывает на необходимость преобразования данных.
      2. Обеспечить точки беспорядочно внутри границ"Остатки против Прогнозируемая" диаграммы. Здесь шеврона образное распределение указывает на необходимость преобразования данных.
      3. Исследовать ли точки в "Остатки против перспективе" схемы разброса случайно в пределах. Образец (например, "лестница") указывает на тенденцию в здании данных и блок может быть использован, чтобы противодействовать его влияние на модели.
      4. Ищите прямой диагональной линии, указывающей идеальную модель в "Прогнозируемая и фактическое" заговора. Чем больше рассеяние данных вокруг диагонали, тем меньше точность модели.
      5. Проверьте перекрытием на ровном месте (лучший) и зеленый (фактические) вертикальные линии в Box-Cox-участка с указанием преобразование данных. Если зеленая линия не соответствует синюю (вне указанном промежутке красными вертикальными линиями) вернуться к разделу 7.3 и улучшить модели здания круглый выбрав преобразование данных, предложенную Box-Cox-участка. Опять дублирование реколонка Отклик может быть полезно для сравнения различных моделей.
      6. Убедитесь в том, что точки в «Остатки по сравнению фактор" диаграмм (есть один график для каждого фактора модели) рассеивают случайно в пределах. Кривизна в распределении данных может указывать на коэффициент отсутствующий в модели.
    2. В разделе "Влияние" из "инструмента диагностики", убедитесь, что есть случайные рассеяние данных в рамках в "Внешне стьюдентизированной остатков", "плечо", "DFFITS" и "DFBETAS" участки и равномерно мало распространен без экстремальные значения в "Кука Расстояние" участка.
  8. На вкладке "Модель графиков", визуализировать оценочную модель. Для ограниченного числа числовых факторов (например, 3) поверхность отклика ("3D поверхность") представление полезно для оценки Optima / характеристики вручную.
    Примечание: поверхности реагирования только проиллюстрировать влияние два факТЗ на ответ под следствием. Эффект дополнительного фактора на ответ раскрывается путем изменения его значения (числовые факторы) или уровень (категориальные факторы) в окне «Факторы инструмента". Кроме того, факторы могут быть отнесены к оси участка щелкнув правой кнопкой мыши их в окне "Факторы инструмента" и выбрав нужный независимой переменной ось.
    1. Манипулировать уровни факторов и назначить их на графике ординат, используя "Факторы функции".
    2. Экспорт графики с помощью "Экспорт график в файл ..." команду на вкладке "Файл".
  9. Используйте "Численное" суб-узел в "оптимизации" узла Для оптимизации чувствительности численно (свернуть, развернуть, в диапазоне) в зависимости от модели факторов, к которым, в свою очередь определенные ограничения могут быть применены (например, ограничения и вес) через Вкладка "Критерии".
    1. Рассчитать и изучить численные решения в "Solutions"Вкладка на основе материалов, представленных на вкладке «критерии».
    2. Экспорт этих решений с другим программным обеспечением (например, электронные таблицы) для дальнейшего анализа (например, гистограмм) выявление установки факторов, связанных с высокой или низкой значений отклика.
      Примечание: это полезно, если более трех числовые факторы исследуются и 3D представление трудно.
  10. Используйте "Точка предсказания" суб-узел, чтобы предсказать ответ для конкретных установок факторов.
    1. Используйте этот прогноз для всех настроек должны быть оценены (здесь, все промотор и 5'UTR комбинации, а также все листья и времени инкубации или листьев позиции и температуру инкубации).
    2. Экспорт предсказанных значений и собрать их в массиве данных, которые могут использоваться для описания временной экспрессии в растениях табака, например, путем вычисления выражения для конкретной комбинации промотор-5'UTR в определенное время, усредненное по всем позиции лист или ACROсс все листья растения.
    3. Нормализовать данные экспрессии на основе масс различных листьев, получая определенный уровень экспрессии (мкг г -1) или удельную скорость экспрессии (мкг г -1 ч -1).
    4. С другой стороны, экспорт коэффициенты "Финал уравнения" из нижней части вкладки "ANOVA" в электронную таблицу и умножить их с массивом установок факторов, получая тот же массив данных.
      Примечание: этот массив может содержать только значения факторов изнутри начальной дизайна помещения, так как оборудована модель уравнение не подходит для экстраполяции.
  11. Провести подтверждение эксперимент для модели точек, которые представляют большой интерес.
    1. Выполните выражение переходный белка в условиях, выбранных для "предсказания Пойнт" (раздел 7.11). Например использовать одинаковые температуры и листья и т.д.
    2. Определить концентрации белка для этого переходного выражения еXperiment, как описано выше (раздел 6) и сравнить их с значений, предсказываемых моделью.
    3. Проверьте, попадает ли средняя концентрация белка в эксперименте подтверждения в модели поверхности отклика в интервале прогнозирования, полученной во время предсказания точки. Матч подтверждает предсказательной силой модели. Несоответствие указывает, может потребоваться низкое качество модели и дополнительные трассы.
    4. Примечание: завод партии к партии изменчивость расширяет интервал прогнозирования и может девальвировать подтверждения экспериментов, выполненных с другой партией растений. Однако образцы, которые не были использованы для создания модели, но которые возникли из той же партии растений как образцы содержали в модели может служить в качестве эффективного контроля.

Subscription Required. Please recommend JoVE to your librarian.

Representative Results

Описательная модель для накопления DsRed во временной экспрессии с использованием различных промоутеров и 5'UTRs

DsRed флуоресценции в экстрактах листьев используется для указания уровня экспрессии рекомбинантного белка и, таким образом, был использован в качестве ответа в стратегии Доу. Минимальная определяемая разница мы рассматривали отношение было 20 мкг / мл и оценивается стандартное отклонение системы было 8 мкг / мл на основе первоначальных экспериментов. Факторы, включенные в переходном модели выражения были отобраны на основе литературы 7,8 данных и наших предыдущих результатов 9. Исследуемые диапазоны также были выбраны в соответствии с этими данными (табл. 1). По крайней мере, три уровня были выбраны для всех дискретных числовых факторов позволяет рассчитывать квадратичной базовой модели. D-оптимальный алгоритм выбора был выбран для отбора НОО работает, чтобы получить наиболее точные оценки длякоэффициенты регрессионной модели. Конструкция изначально предложено DesignExpert состоял из 90 трасс, но FDS было недостаточно для достижения 1% стандартную ошибку предсказания (рис. 3А). D-оптимальный увеличение конструкции в общей сложности 210 трасс решены эту проблему и привело к FDS 100% с более равномерной точности прогнозирования через дизайн пространстве, обозначенном плоской кривой (рис. 3В).

Концентрации DsRed были определены для всех серий 210 и данные были преобразованы Log 10. Факторы Модель были выбраны автоматизированной обратной отбора из кубического модели с альфа уровне 0.100. Это привело к значительному модели (табл. 2) с незначительными отсутствие подгонки и высокие значения для множества коэффициентов корреляции (табл. 2). Р-значение всех модельных факторов был <0,05 и, таким образом, не было необходимости дальнейшего ручного манипулирования модели. Модель содержитсятри-фактор взаимодействия, которые не были частью первоначального квадратичной базовой модели (табл. 2). Переоценка графике FDS используя все факторы, включенные в окончательной модели прогнозирования показали, что FDS для стандартной ошибки прогноза существенно не уменьшилась, включив дополнительные взаимодействий трех факторов (FDS = 0,99).

Средства диагностики качества модель в DesignExpert указано, что преобразование данных была полезной и не было никаких недостающие факторы в модели, потому что нормальный участок остатков показал линейное поведение и никаких конкретных шаблон не наблюдалось остатков против прогнозируемого участка (5А и 5Б ). Там также не было тенденции на протяжении всего эксперимента, чтобы показать скрытую зависящих от времени переменной (рис. 5C). Вместо этого предсказания модели были в очень хорошем согласии с наблюдаемым DsRed флуоресценции (рис. 5D). Мы тherefore предполагается, что выбранная модель была полезна для прогнозирования временную экспрессию DsRed в не-семядолей листьев табака приводимых в другой промотор / 5'-UTR комбинаций во время после инфильтрации инкубационного периода продолжительностью 8 дней. Мы также выбрали искусственную модель линейной регрессии без преобразования данных для иллюстрации последствий неправильного выбора фактора и трансформации (рис. 6). Очевидно, что нормальный участок остатков отклоняется от ожидаемого линейного поведения (фиг.6А) и есть "V"-образной кривой в остатках против предсказано участок вместо случайного разброса (фиг.6В). Кроме того, остатки против запуска участок подчеркивает два крайних значения (рис. 6в), в то время как прогнозы бедны для обоих, малых и больших значений, отклоняясь от оптимальной модели (по диагонали) (рис. 6D).

Комбинации 5'UTR с CaMV 35SS промотор привело к усилению экспрессии DsRed, чем комбинации с NOS промотор (фиг. 5A и 5B), как сообщалось ранее 34 и соответствующие факторы были признаны значимыми в модели DoE (табл. 2). Модель также предсказывает, что лист возраст является важным фактором (табл. 2) с уровнями экспрессии, показанных выше в молодых листьях (рис. 7В и 7D), который был в хорошем согласии с нашими предыдущими выводами 19 и других людей 7,8. Прогрессирование накопления DsRed в листьях не было линейным или экспоненциальным но затем сигмоидальную кривую в течение 8 дней после инфильтрации инкубации в зависимости от модели. Интересно, что не было никакой фиксированной порядок 5'UTRs через соответствующие выражения DsRed. Таким образом, «сила» в 5'-UTR зависела от сопутствующей промоутера. Маловероятно, что этовзаимозависимость между промотором и 5'-UTR бы были выявлены с помощью OFAT эксперимент. Прогнозная модель также указал, что некоторые пары промотор / 5'UTR комбинаций (например. CaMV 35SS/CHS и н.у.к. / CHS (7А и 7В), CaMV 35SS/omega и NOS / CHS или CaMV 35SS/CHS и CaMV 35SS / TL) привело к сбалансированной уровней экспрессии, отличаясь менее чем на 30% по сравнению с определенном соотношении во всех листьев и времени инкубации> 2 дня (например,. для CaMV 35SS/omega и н.у.к / CHS соотношение было ~ 8,0) 20. Такой сбалансированный выражение было бы полезно для экспрессии мультимерных белков, таких как IgA с определенной стехиометрией 35,36.

Рисунок 5
Рисунок 5. Graphical оценка качества модели. . Нормальный вероятностный график внутренне стьюдентизированной остатков указывает нормальное распределение, потому что остатки следовать линии. B. Внутренне стьюдентизированные остатки разбрасывать вокруг нулевого значения для всех диапазонов прогнозируемой реакции (флуоресценции) с указанием соответствующего преобразования данных. C. Внутренне стьюдентизированные остатки разбрасывать вокруг нулевого значения в течение всего хода эксперимента, показывая отсутствие скрытых временных эффектов. D. Прогнозируемая и фактические значения ответа находятся в хорошем согласии как все точки лежат близко к диагонали (идеальной модели).

Рисунок 6
Рисунок 6. Диагностика с указанием низкое качество модели. B. Внутренне стьюдентизированные остатки показывают "V"-образный распределение указывает неуместным преобразования данных. C. Внутренне стьюдентизированные остатки рассеивают не вокруг нулевого значения в течение всего хода эксперимента, но проявляют тенденцию к положительным значениям. Кроме того, два крайних значения могут быть найдены. D. Прогнозируемая и фактические значения отклика плохо согласуются как наибольшее количество очков отклоняться от диагонали (идеальной модели).

Рисунок 7
Рисунок 7. Ответ поверхностей для временной экспрессии DsRed в листьях табака. . Чонг> неопределенные / CHS в лист 2; B. CaMV 35SS/CHS в лист 2; C. CaMV 35SS/TL в лист 6, D. CaMV 35SS/CHS в лист 6. Концентрации DsRed в течение инкубационного периода увеличилась в сигмоидальной образом. Уровни экспрессии были ниже в старых листьев (например, лист 2 в А и В) по сравнению с молодыми листьями (например листьев 6 в С и D) и для NOS (A) по сравнению с промотором CaMV 35SS (B). 5'UTR также имели значительное влияние на DsRed выражения (например, TL (C) по сравнению с CHS (D)), но сила выражение зависит от сопровождающего промоутера.

"> DF
Источник Сумма квадратов Средний квадрат F-значение р-значение
Модель 174.85 95 1.84 182.12 <0,0001
Позиция на листе (А) 0.16 1 0.16 16.06 0.0001
Время инкубации (В) 112.46 1 112.46 11128,22 <0,0001
Лист возраст (С) 15.24 7 2.18 215.39 <0,0001
Организатор (D) 23.49 1 23.49 2324.29 <0,0001
5'-UTR (E) 0.93 3 0.31 30.61 <0,0001
Переменный ток 0.24 7 0.034 3.38 0,0026
До нашей эры 1.46 7 0.21 20.64 <0,0001
BD 2.27 1 2.27 224.51 <0,0001
BE 0.87 3 0,29 28.74 <0,0001
Компакт-диск 0,29 7 0.042 4.11 0.0005
CE 0.43 21 0.021 2.04 0,0093
Делавэр 0.48 3 0.16 15.73 <0,0001
В2 6.34 1 6.34 627,29 <0,0001
BCD 0.95 7 0,14 13.42 <0,0001
До н.э. 0,39 21 0.019 0.0230
BDE 0.16 3 0.054 5.37 0.0017
B2D 1.49 1 1.49 147.93 <0,0001
Остаточный 1.15 114 0.010
Отсутствие подгонки 1.08 104 0.010 1.45 0,2669
Чистая ошибка 0.072 10 7.153E-003
Общая Корреляция 176.01 209
Коэффициенты корреляции Значение
R-Squared 0.9935
Скорректированный R-Squared 0.9880
Прогнозируемая R-Squared 0.9770

Таблица 2. Факторы, включенные в прогностическаямодель для временной экспрессии DsRed в табачных листьев используя различные комбинации промотор / 5 'UTR. Три-фактор взаимодействия выделены жирным шрифтом.

Оптимизация условий инкубации и схемы сбора урожая для производства моноклональных антител в растениях по временной экспрессии

Подход DoE был также использован для оптимизации температуры инкубации, бактериальный OD600nm для схем лист инфильтрации, возрастных растений и сбора урожая, для одновременного производства моноклональных ВИЧ-нейтрализующих антител 2G12 и DsRed табаком 19. Схема урожай определили, какие были использованы листья для извлечения белка, например шести верхних листьев. Поэтому мы установили прогностическую модель для выражения каждого белка (2G12 и DsRed) в растениях в разном возрасте (стадия молодого = рано бутон = урожая на 40 дней после посева; старый = ступень зрелый бутон = урожая на 47 дней после см.динь). Эти четыре модели были затем оценены и общепринятую модель установлено, что включен каждый фактор оказывается значимым в отдельных моделях. Затем мы подтвердили, что модель консенсуса по-прежнему хорошее представление всех начальных наборов данных (табл. 3). Затем консенсус модель использовалась для определения оптимальных температуры инкубации (~ 22 ° C для 2G12 и ~ 25 ° C для DsRed) и бактериальный OD 600 нм (~ 1.0) для обоих белков. Эти оптимальные условия были затем использованы для прогнозирования концентрации белка во всех листьев (1-8) и листовых позиций (1-4) в молодых и старых растений. Профили концентрации были объединены с данными биомассы 19, получая абсолютное количество белка-мишени, полученной из различных листьев и возраста растений (фиг.8А). Абсолютные суммы белковые затем были сопоставлены с соответствующими вниз по течению затрат позволяет нам проводить анализ затрат и выгод для обработки каждого возраста лист / растений. Это показало, чтомолодые растения были выгодны для временной экспрессии, так как белки достиг более высокие концентрации в течение коротких периодов роста, несмотря на более низкой общей биомассы по сравнению с старых растений (Цифры 8В и 8С). Мы также обнаружили, что обработка всех листьев растений из старых был дороже, чем отбрасывания листьев 1-3 и увеличения количества растений в пакете вместо (рис. 8d). Таким образом, модели DoE основе подходят не только для обозначения заключительный шаг эксперимента, но и для комбинации с другими данными, чтобы облегчить более сложные аспекты анализа процессов. Серия взаимосвязанных моделей, охватывающих весь процесс производства на биофармацевтической белка может быть конечной целью.

Рисунок 8
Рисунок 8. Интегрируя ион биомассы и белковых данных о концентрации в результате в абсолютных выход белка и производственных затрат. . Лист биомассы и концентрация белка цель не развиваться в коллинеарном образом, в результате предвзятого накопления 2G12 в молодых листьях. B. То же самое количество DsRed и ~ 65% 2G12 было найдено у молодых растений по сравнению с старых несмотря на ~ 50% меньше среднего биомассы. С. Это отражало более высокую удельную экспрессию белка в молодых растений. D. Повышенная специфическая экспрессия переведены на снижение издержек производства для обоих белков, потому что меньше биомассы необходимы для обработки требует меньше парниковых пространства, меньше расходных материалов (фильтры, например.), И менее последующего оборудования.

пт; ширина: 48pt "ширина =" 64 "> возрастные завод [ДПС] т ">
40 47
Целевой белок [-] DsRed 2G12 DsRed 2G12
Модель [-] Первоначальный Консенсус Первоначальный Консенсус Первоначальный Консенсус Первоначальный Консенсус
R-квадрат [-] 0.9829 0,9577 0,9321 0,9099 0,9436 0,9403 0,8826 0.8782
Скорректированный R-квадрат [-] 0,9711 0.9480 0.9059 0,8893 0.9362 0,9350 0.8716 0.8674
Прогнозируемая R-квадрат [-] 0.9510 0.9336 0.8272 0,8587 0,9254 0,9282 0,8554 0,8516

Таблица 3. Сравнение коэффициентов корреляции для начальных RSМ модели и окончательный консенсус модель DsRed и 2G12 выражения в растениях табака.

Subscription Required. Please recommend JoVE to your librarian.

Discussion

Каждый эксперимент требует тщательного планирования, поскольку ресурсы часто дефицитным и дорогим. Это особенно верно для стратегий Доу, так как ошибки на этапе планирования (например, выбора базовой модели, которая не покрывает все значимые взаимодействия факторов) может существенно уменьшить предсказательной силой получаемых моделей и тем самым обесценить весь эксперимент. Тем не менее, эти ошибки можно легко избежать путем следующих основных процедур.

Соображения при планировании DoE

Во-первых, важно, чтобы выбрать ответ, который подходит для оценки результатов (положительный или отрицательный) для каждого Doe Run. Например: концентрация белка-мишени, определяемой поглощения УФ 37 полезно оценить активность определенной комбинации промотор / 5'UTR. Тем не менее, функциональная оценка белка (например. По флуоресценции в случае DsRed или связывания с белком А в КАСе 2G12) еще лучше, потому что это также включает в себя качество белка аспект. Реакции также могут быть значения, рассчитанные из двух или более отдельных параметров, таких как количество электроэнергии НП в ферментационных экспериментов (в расчете от входной мощности, скорости вращения и диаметра мешалки). Ответ должен дать значения с высокой точностью (т.е. минимальному стандартное отклонение) и точности (т.е. ноль или несколько мешающих параметров) по улучшению качества модели. Поэтому Подходящие устройства и способы обнаружения должна быть оптимизирована до эксперимента, если необходимо.

Факторы, влияющие на реакцию, должны быть определены до того, как Министерство энергетики реализуется, потому что РСМ используется, чтобы определить точное воздействие факторов, а не выбирать их. Выбор факторов со значительным влиянием на ответ может быть достигнуто путем изучения литературных данных, но полные или дробные факторные стратегии DoE являются более эффективными. Эти фильтры эксперименты должнытакже можно использовать для определения диапазона для каждого фактора в котором значимые результаты могут быть достигнуты, например, в биологических системах фактор "температура" часто ограничивается диапазоне 10-40 ° С. Дискретных уровней должны быть определены в этих пределах для числовых факторов, которые технически трудно контролировать, например, фитотроне температуры могут колебаться в пределах ± 2,0 ° С, поэтому нецелесообразно исследовать разницу температур ± 1,0 ° С. Практические ограничения могут также вступают в игру, например, с помощью> 25 различных концентраций конкретного соединения. Тем не менее, число уровней должно быть по крайней мере один больше, чем ожидалось степенной порядок этого фактора в базовой модели, например, если многочлен порядок температуре п = 2, то по крайней мере три уровня температуры должны быть расследованы. Использование базовых моделей высокого степенной порядок (кубической или выше) является вариантом, если многочлен порядок фактор является неопределенным.Тем не менее, это может существенно увеличить количество трасс, необходимых для DoE.

Для сложных систем с многими факторами и высоких порядков полиномиальных это может привести к конструкции, которые не могут быть вычислены эффективно. В таких случаях может быть целесообразным разделить проблему на две или более конструкций следственных различные подмножества факторов. Тем не менее, существенные факторы исключены из таких «сплит-конструкций" должны быть тщательно доводят до фиксированных значений для предотвращения искажения результатов. Каждый сплит-дизайн должен быть повторен при различных настройках для опущенных факторов выявить взаимодействия между ними и факторов, включенных в проект. Это выгодно, чтобы определить факторы, оказавшая наибольшее влияние на ответ и включить их во всех конструкций.

Критические аспекты во DoE оценки

FDS должна быть первый параметр оценивается в недавно рассчитанной DoE. Оператор должен подтвердить, что THэ DoE обеспечивает необходимую охват дизайн помещения для требуемого минимального обнаруживаемого различия и стандартное отклонение ответа. Если это не так, то конструкция должна быть увеличена с помощью соответствующего числа серий и повторяет, пока требования не будут выполнены. Процедура для расчета Doe работает можно инициализировать несколько раз, чтобы выбрать проекты с минимальными потерями оптимальности, если баланс между категориальных факторов включена.

Преобразование данных должно осуществляться, если значения реагирования занимать более одной порядок. Вход 10 преобразование часто целесообразно но сюжет Бокса-Кокса будет предложить трансформацию, которая даст наиболее стабильную модель. Преобразование предложил в этом участке может измениться после того, как окончательная модель была выбрана, и поэтому должны быть перепроверены. Только значимые факторы (р-значение <0,05) должны быть включены в модели поверхности отклика для того, чтобы гарантировать ее интеллектуального рАуэр. Исключения могут быть сделаны для поддержания иерархии модели 32,33 или включить факторы, которые, как известно, оказывают влияние на основе предыдущих результатов. В последнем случае, желательно, чтобы исследовать, почему такие важные факторы, не были признаны значимыми в текущей Доу. Отсутствие подгонки параметр может быть использован в сочетании с R-квадрат величины (табл. 2 и 3) для быстрого обнаружения плохих моделей / припадков, тогда как скорректированный R-квадрат указывает, были ли выбраны слишком много / Незначительные факторы и предсказал R-квадрат является показателем предсказательной силы модели. Тем не менее, остаточные участки, предоставленные в разделе диагностическом DesignExpert еще более важны для оценки качества модели, поскольку они позволяют обнаруживать тенденции в наборе данных, свидетельствующих о наличии скрытых факторов, а также идентификацию экстремальных значений. Последние имеют выше среднего влияние на посадку и исказить модель. Таких значений может исходить от некорректных данных (например перестановкой значений), которые должны быть исправлены, или из неудачных экспериментов, которые должны быть повторены. Если экстремальное значение оказывается "реальный" в таких экспериментах повторения желательно добавить несколько трасс в непосредственной близости от этого значения, чтобы улучшить качество модели в пределах этой области дизайна помещения.

Поправки к конструкции, которые несовершенны или слишком маленький

Независимо от почему дополнительные пробеги могут потребоваться (например, отсутствуют данные, удалось трасс, подлинные экстремальные значения и регионов, непредвиденные взаимодействия факторов высокого порядка), эти трассы должны быть добавлены к существующей конструкции в упорядоченной манере с помощью "Дизайн увеличение" инструмент ввести новые трассы в отдельном блоке. Это позволит облегчить использование предыдущего набора данных при выборе наиболее значимые дополнительные рейсы. Результат будет модель, которая подходит к данным и позволяет надежнымпрогнозирование будущих результатов.

Subscription Required. Please recommend JoVE to your librarian.

Disclosures

Пошлина за публикацию частично спонсируется компаниями Statease, Inc (США) и Statcon (Германия), которые не были вовлечены в участие в подготовке рукописи или ответственность за его содержание.

Acknowledgments

Авторы благодарны д-р Томас Радемахер для обеспечения растений вектор экспрессии PPAM и Ибрахима аль Amedi для выращивания табачные растения, используемые в данном исследовании. Мы хотели бы поблагодарить д-ра Ричарда М. Twyman за его помощь в редактировании рукописи. Эта работа была частично финансируется Европейским исследовательским советом Расширенный Грант "Future-Фарма", число предложений 269110 и Фраунгофера Zukunftsstiftung (Fraunhofer Будущее Фонд).

Materials

Name Company Catalog Number Comments
Design-Expert(R) 8 Stat-Ease, Inc. n.a. DoE software
Tryptone Carl Roth GmbH 8952.2 Media component
Yeast extract Carl Roth GmbH 2363.2 Media component
Sodium chloride Carl Roth GmbH P029.2 Media component
Ampicillin Carl Roth GmbH K029.2 Antibiotic
Agar-Agar Carl Roth GmbH 5210.2 Media component
Escherichia coli K12 DH5a Life Technologies 18263-012 Microorganism
pPAM GenBank AY027531 Cloning/expression vector; 
NucleoSpin Plasmid  MACHEREY-NAGEL GmbH 740588.250 Plasmid DNA isolation kit
NucleoSpin Gel and PCR Clean-up MACHEREY-NAGEL GmbH 740609.250 Plasmid DNA purification kit
NanoDrop 2000 Thermo Scientific n.a. Spectrophotometer
NcoI New England Biolabs Inc. R3193L Restrictionendonuclease
EcoRI New England Biolabs Inc. R3101L Restrictionendonuclease
AscI New England Biolabs Inc. R0558L Restrictionendonuclease
NEB 4 New England Biolabs Inc. B7004S Restrictionendonuclease buffer
TRIS Carl Roth GmbH 4855.3 Media component
Disodium tetraborate Carl Roth GmbH 4403.3 Media component
EDTA Carl Roth GmbH 8040.2 Media component
Agarose Carl Roth GmbH 6352.4 Media component
Bromophenol blue Carl Roth GmbH A512.1 Color indicator
Xylene cyanol Carl Roth GmbH A513.1 Color indicator
Glycerol Carl Roth GmbH 7530.2 Media component
Mini-Sub Cell GT Cell BioRad 170-4406 Gel electrophoresis chamber
Agrobacterium tumefaciens strain GV3101:pMP90RK DSMZ 12365 Microorganism
Electroporator 2510 Eppendorf 4307000.658 Electroporator
Beef extract Carl Roth GmbH X975.2 Media component
Peptone Carl Roth GmbH 2365.2 Media component
Sucrose Carl Roth GmbH 4621.2 Media component
Magnesium sulfate Carl Roth GmbH 0261.3 Media component
Carbenicillin Carl Roth GmbH 6344.2 Antibiotic
Kanamycin Carl Roth GmbH T832.3 Antibiotic
Rifampicin Carl Roth GmbH 4163.2 Antibiotic
FWD primer Eurofins MWG Operon n.a. CCT CAG GAA GAG CAA TAC
REV primer Eurofins MWG Operon n.a. CCA AAG CGA GTA CAC AAC
2720 Thermal cycler Applied Biosystems 4359659 Thermocycler
RNAfold webserver University of Vienna n.a. Software
Ferty 2 Mega Kammlott 5.220072 Fertilizer
Grodan Rockwool Cubes 10 x10 cm Grodan n.a. Rockwool block
Greenhouse n.a. n.a. For plant cultivation
Phytotron Ilka Zell n.a. For plant cultivation
Omnifix-F Solo B. Braun 6064204 Syringe
Murashige and Skoog salts Duchefa M 0222.0010 Media component
Glucose Carl Roth GmbH 6780.2 Media component
Acetosyringone Sigma-Aldrich D134406-5G Phytohormon analogon
 BioPhotometer plus Eppendorf 6132 000.008 Photometer
Osram cool white 36 W Osram 4930440 Light source
Disodium phosphate Carl Roth GmbH 4984.3 Media component
Centrifuge 5415D Eppendorf 5424 000.410 Centrifuge
Forma -86C ULT freezer ThermoFisher 88400 Freezer
Synergy HT BioTek SIAFRT Fluorescence plate reader
Biacore T200 GE Healthcare n.a. SPR device
Protein A Life Technologies 10-1006 Antibody binding protein
HEPES Carl Roth GmbH 9105.3 Media component
Tween-20 Carl Roth GmbH 9127.3 Media component
2G12 antibody Polymun AB002 Reference antibody

DOWNLOAD MATERIALS LIST

References

  1. Fischer, R., Emans, N. Molecular farming of pharmaceutical proteins. Transgenic research. 9, 277-299 (2000).
  2. Commandeur, U., Twyman, R. M., Fischer, R. The biosafety of molecular farming in plants. AgBiotechNet. 5, 9 (2003).
  3. Shoji, Y., et al. A plant-based system for rapid production of influenza vaccine antigens. Influenza Other Resp. 6, 204-210 (2012).
  4. Goodin, M. M., Zaitlin, D., Naidu, R. A., Lommel, S. A. Nicotiana benthamiana: Its history and future as a model for plant-pathogen interactions. Mol Plant Microbe In. 21, 1015-1026 (2008).
  5. Berg, R. H., Beachy, R. N. Fluorescent protein applications in plants. Method Cell Biol. 85, 153 (2008).
  6. Chung, S. M., Vaidya, M., Tzfira, T. Agrobacterium is not alone: gene transfer to plants by viruses and other bacteria. Trends in plant science. 11, 1-4 (2006).
  7. Sheludko, Y. V., Sindarovska, Y. R., Gerasymenko, I. M., Bannikova, M. A., Kuchuk, N. V. Comparison of several Nicotiana species as hosts for high-scale Agrobacterium-mediated transient expression. Biotechnology and Bioengineering. 96, 608-614 (2007).
  8. Wydro, M., Kozubek, E., Lehmann, P. Optimization of transient Agrobacterium-mediated gene expression system in leaves of Nicotiana benthamiana. Acta Biochimica Polonica. 53, 289-298 (2006).
  9. Buyel, J. F., Fischer, R. Processing heterogeneous biomass: Overcoming the hurdles in model building. Bioengineered. 4, (2013).
  10. Fischer, R., Schillberg, S., Hellwig, S., Twyman, R. M., Drossard, J. GMP issues for recombinant plant-derived pharmaceutical proteins. Biotechnol Adv. 30, 434-439 (2012).
  11. Daniel, C. One-at-a-time plans. Journal of the American Statistical Association. 68, 353-360 (1973).
  12. Czitrom, V. One-Factor-at-a-Time versus Designed Experiments The American Statistician. 53, 6 (1999).
  13. Anderson, M. J., Kraber, S. L. Keys to successful designed experiments. ASQ - The global voice of quality. 6, 6 (1999).
  14. Montgomery, D. C. Design and Analysis of Experiments. , John Wiley & Sons Inc. (2007).
  15. Myers, R. H., Montgomery, D. C., Anderson-Cook, C. M. Response Surface Methodology: Process and Product Optimization Using Designed Experiments. , Wiley. (2009).
  16. Piepel, G. F. Programs for generating extreme vertices and centroids of linearly constrained experimental regions. J Qual Technol. 20, 15 (1988).
  17. FDA. , U.S. Department of Health and Human Services. Rockville, MD, USA. (2009).
  18. Shivhare, M., McCreath, G. Practical Considerations for DoE Implementation in Quality By Design. BioProcess International. 8, 9 (2010).
  19. Buyel, J. F., Fischer, R. Predictive models for transient protein expression in tobacco (Nicotiana tabacum L.) can optimize process time, yield, and downstream costs. Biotechnology and bioengineering. 109, 2575-2588 (2012).
  20. Buyel, J. F., Kaever, T., Buyel, J. J., Fischer, R. Predictive models for the accumulation of a fluorescent marker protein in tobacco leaves according to the promoter/5'UTR combination. Biotechnology and bioengineering. 110, 471-482 (2013).
  21. Anderson, M. J., Whitcomb, P. J. DOE Simplified: Practical Tools for Effective Experimentation. , Productivity Inc. (2000).
  22. Anderson, M. J., Whitcomb, P. J. Response Surface Methods Simplified. , Productivity Press. (2005).
  23. De Gryze, S., Langhans, I., Vandebroek, M. Using the correct intervals for prediction: A tutorial on tolerance intervals for ordinary least-squares regression. Chemometr Intell Lab. 87, 147-154 (2007).
  24. Plasmid DNA purification User manual. , MACHEREY-NAGEL. Düren. (2012).
  25. PCR clean-up Gel extraction User manual. , MACHEREY-NAGEL. Düren. (2012).
  26. Quick Ligation Protocol. 4, NewEnglandBiolabs. Ipswich. (2009).
  27. Inoue, H., Nojima, H., Okayama, H. High-Efficiency Transformation of Escherichia-Coli with Plasmids. Gene. 96, 23-28 (1990).
  28. Main, G. D., Reynolds, S., Gartland, J. S. Electroporation protocols for Agrobacterium. Methods in Molecular Biology. 44, 405-412 (1995).
  29. Gruber, A. R., Lorenz, R., Bernhart, S. H., Neubock, R., Hofacker, I. L. The Vienna RNA websuite. Nucleic acids research. 36, 70-74 (2008).
  30. Howell, S., Kenmore, M., Kirkland, M., Badley, R. A. High-density immobilization of an antibody fragment to a carboxymethylated dextran-linked biosensor surface. J Mol Recognit. 11, 200-203 (1998).
  31. Newcombe, A. R., et al. Evaluation of a biosensor assay to quantify polyclonal IgG in ovine serum used for the production of biotherapeutic antibody fragments. Process Biochem. 41, 842-847 (2006).
  32. Peixoto, J. L. Hierarchical Variable Selection in Polynomial Regression-Models. Am Stat. 41, 311-313 (1987).
  33. Peixoto, J. L. A Property of Well-Formulated Polynomial Regression-Models. Am Stat. 44, 26-30 (1990).
  34. Sanders, P. R., Winter, J. A., Barnason, A. R., Rogers, S. G., Fraley, R. T. Comparison of cauliflower mosaic virus 35S and nopaline synthase promoters in transgenic plants. Nucleic acids research. 15, 1543-1558 (1987).
  35. Ma, J. K. C., et al. Generation and Assembly of Secretory Antibodies in Plants. Science. 268, 716-719 (1995).
  36. Wycoff, K. L. Secretory IgA antibodies from plants. Curr Pharm Design. 11, 2429-2437 (2005).
  37. Pace, C. N., Vajdos, F., Fee, L., Grimsley, G., Gray, T. How to measure and predict the molar absorption coefficient of a protein. Protein Sci. 4, 2411-2423 (1995).

Tags

Биоинженерия выпуск 83 дизайн экспериментов (DOE) выражение переходный белка растительного происхождения биофармацевтических промоутер 5'UTR люминесцентные репортер белка построение модели условия инкубации моноклональное антитело
Характеристика сложных систем по постановке эксперимента подхода: Переходный Белок Expression табачными в качестве примера
Play Video
PDF DOI DOWNLOAD MATERIALS LIST

Cite this Article

Buyel, J. F., Fischer, R.… More

Buyel, J. F., Fischer, R. Characterization of Complex Systems Using the Design of Experiments Approach: Transient Protein Expression in Tobacco as a Case Study. J. Vis. Exp. (83), e51216, doi:10.3791/51216 (2014).

Less
Copy Citation Download Citation Reprints and Permissions
View Video

Get cutting-edge science videos from JoVE sent straight to your inbox every month.

Waiting X
Simple Hit Counter