Бесплатный автореферат и диссертация по биологии на тему

Влияние аминокислотных замен на спектральные свойства флуоресцентных белков на примере mRFP1

ВАК РФ 03.01.02, Биофизика

Автореферат диссертации по теме "Влияние аминокислотных замен на спектральные свойства флуоресцентных белков на примере mRFP1"

005015319

На правах рукописи

Вржещ Евгений Петрович

Влияние аминокислотных замен на спектральные свойства флуоресцентных белков на примере

т^Р!

03.01.02-биофизика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

1 2 мдр 2012

МОСКВА 2012

005015319

Работа выполнена на кафедре общей физики Физического факультета и Факультете биоинженерии и биоинформатики Московского Государственного Университета им. М.В. Ломоносова.

Научный руководитель:

Защита состоится 22 марта 2012 года в 14 часов 00 минут на заседании Диссертационного совета Д 501.001.96 при Московском Государственном Университете им. М.В. Ломоносова по адресу: 119991, Россия, г. Москва, Ленинские горы 1/12, МГУ, биологический факультет, кафедра биофизики, аудитория «Новая».

С диссертацией можно ознакомиться в библиотеке биологического факультета Московского Государственного Университета им. М.В. Ломоносова.

доктор физико-математических наук, профессор Салецкий Александр Михайлович заведующий кафедрой общей физики Физического факультета МГУ

Официальные оппоненты:

Савицкий Александр Павлович

доктор химических наук, профессор руководитель лаборатории физической биохимии Института биохимии им. А.Н. Баха РАН

Григоренко Белла Людвиговна

доктор физико-математических наук старший научный сотрудник Химического факультета МГУ

Ведущая организация:

Оренбургский Государственный Университет.

Автореферат разослан <<_> февраля 2012 года.

Ученый секретарь

Диссертационного совета Д 501.001.96 Доктор биологических наук

Страховская М.Г.

Общая характеристика работы

Актуальность темы

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

Флуоресцентные белки активно используются в исследовании прижизненных процессов в живых клетках, тканях и организмах, таких как активность работы промотеров, биохимические взаимодействия, локализации различных структур и др. Ген флуоресцентного белка с помощью стандартных биохимических методов вносится в генетический аппарат исследуемого объекта и дальнейшая экспрессия белка и его свечение не требуют наличия дополнительных кофакторов или специфических реагентов (кроме атмосферного кислорода). Набор существующих на данный момент флуоресцентных белков достаточно разнообразен: максимум эмиссии флуоресценции белков занимает широкую область спектра: от синей (менее 450 нм) до красной (708 нм) частей видимого излучения. Также различным флуоресцентным белкам свойственно большое разнообразие значений таких параметров, как квантовый выход флуоресценции, коэффициент экстинкции, фотостабильность и пр.

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

Цели и задачи исследования

Целью данного исследования является изучение влияния аминокислотных замен на структуру и свойства мутантов красного флуоресцентного белка mRFPl in silico и in vitro с использованием методов квантовой химии, молекулярной динамики, генной инженерии и спектрально-люминесцентной спектроскопии. Исходя из этого, были сформулированы следующие задачи исследования:

• создать топологию хромофора белка mRFPl как нового аминокислотного остатка в силовом поле и вычислить ее параметры;

• создать PDB-файлы для мутантов белка mRFPl, отредактировать файлы силового поля и создать топологии, соответствующие хромофорам мутантов;

• рассчитать равновесные трехмерные структуры мутантов белка mRFPl и провести анализ влияния замены аминокислотного остатка глутамина по 66 положению на геометрические параметры белка;

• методом генной инженерии создать плазмиды, содержащие гены 20 мутантных белков на основе белка mRFPl со всеми вариантами аминокислотных остатков по 66 положению, провести целевую наработку препаратов белка, определить оптические свойства полученных препаратов;

• провести анализ влияния замены аминокислотного остатка глутамина по 66 положению на оптические свойства белка и сопоставить с изменением геометрических параметров;

• методами нелинейной лазерной флуориметрии установить истинные фотофизические параметры отдельных молекул белков и концентрации различных форм белка в растворе.

Научная новизна результатов и практическая значимость работы

Разработан и опробован метод генно-инженерного внесения любых наперед заданных мутаций в ген красного мономерного флуоресцентного белка mRFPl в области хромофора, позволяющий проводить направленный мутагенез белков с целью скрининга свойств в зависимости от изменения аминокислотной последовательности. С помощью метода впервые получено девятнадцать уникальных новых мутантов белка mRFPl по 66 положению в виде концентрированных препаратов белка и определены их оптические свойства.

Впервые разработан метод описания и параметризации объектов типа красных флуоресцентных белков в силовом поле OPLS-AA (в частности, хромофора как отдельной

уникальной аминокислоты) для проведения молекулярно-динамических расчетов состояния и геометрии белков в растворе т яШсо.

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

Разработан метод, позволяющий определять фотофизические параметры индивидуальных форм хромофоров в составе смеси без ее разделения, что позволит более точно следить за характеристиками как при получении новых вариантов флуоресцентных белков, так и при применении их как маркеров процессов в живых системах.

Впервые определенные данным методом фотофизические характеристики белка тЫЛ и некоторых его мутантов показали, что истинный коэффициент экстинкции отдельных форм флуоресцентных белков может значительно превышать 100 ООО М"'см"'.

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

Положения, выносимые на защиту

1. предложен метод генно-инженерного внесения любых наперед заданных мутаций в ген красного мономерного флуоресцентного белка тЮЛ в области хромофора, позволяющий проводить направленный мутагенез с целью скрининга свойств в зависимости от изменения аминокислотной последовательности;

2. разработан метод описания и параметризации объектов типа красных флуоресцентных белков в силовом поле ОРЬБ-АА, позволяющий проводить молекулярно-динамические расчеты состояния и геометрии белков в растворе т яШсо;

3. выявлены эмпирические зависимости оптических свойств флуоресцентных белков от геометрических параметров аминокислот, составляющих и/или окружающих хромофор, что может применяться при планировании мутагенеза, направленного на получение новых флуоресцентных белков с улучшенными свойствами;

4. разработан метод, позволяющий определять фотофизические параметры индивидуальных форм хромофоров в составе смеси без ее разделения.

Апробация работы

Основные результаты работы были доложены на конференции «Альманах клинической медицины» в 2006 году, международных конференциях студентов, аспирантов и молодых ученых «Ломоносов», студенческих конференциях Факультета Биоинженерии и Биоинформатики МГУ в 2006-2007 годах.

Публикации

По результатам работы опубликовано 9 работ, из них статей в рецензируемых журналах, включенных в список ВАК РФ - 5, тезисов докладов и материалов конференций - 4.

Личный вклад автора

Все основные результаты работы получены лично диссертантом. Вклад диссертанта в работу является определяющим.

Структура диссертации

Диссертационная работа состоит из введения, пяти глав результатов, выводов и списка цитируемой литературы. Объем диссертации составляет 127 страниц, в том числе 35 рисунков, 21 таблица и 27 формул. Список цитируемой литературы включает 94 наименования.

Содержание работы

Во введении обоснована актуальность темы диссертации, ее научная новизна и практическая значимость, а также сформулирована цель и поставлены задачи исследования.

Первая глава содержит обзор опубликованных данных по теме диссертации. Описана история открытия флуоресцентных белков и основные природные источники. Приведена классификация белков по семействам с указаниями спектральных и структурных свойств. Описывается биологическая роль флуоресцентных белков в различных организмах. Проведен анализ данных о кристаллических структурах белков различных семейств. Описываются механизмы формирования (созревания) хромофоров зеленого СРР и красного ОвЯес! флуоресцентных белков, указывается различие в количестве стадий и промежуточных форм в механизмах созревания. Анализируется влияние аминокислотного окружения будущего хромофора на механизм созревания. Описаны флуоресцентные белки с различными олигомерными состояниями, на примере белка ОвИ-ес! показаны интерфейсы взаимодействия различных мономеров при образовании тетрамера. Описан механизм мутагенеза, приводящий к образованию красного мономерного флуоресцентного белка тЯБР1 из тетрамерного ОвЯес!. Описаны физические и оптические характеристики белков СРР, БвЯес! и т!1РР1. Проведен анализ данных о влиянии замен (как на природные, так и на искусственные аминокислотные остатки) по различным позициям хромофора на флуоресцентные свойства белка. В заключение обзора литературы формулируется важность понимания взаимосвязи между первичной структурой белков и оптическими свойствами.

Во второй главе диссертации описываются методы молекулярной динамики (МД), конкретная их реализация для пакета программ СЯОМАСБ, синтаксис и параметры силового поля ОРЬЭ-АА, Описывается применение метода МД для нахождения равновесных структур флуоресцентных белков в растворе, приведены все параметры для запуска МД-симуляций и анализа полученных траекторий. Использовались программный пакет ОЯОМАСЭ и силовое поле ОРЬБ-АА в виде его модификации ОРЬв-АА/ОзЛес!, содержащей параметры топологий всевозможных хромофоров красного типа как отдельных уникальных аминокислотных остатков. Указаны все стадии, уравнения и параметры расчета.

В третьей главе описывается создание топологий всех вариантов хромофоров для проведения молекулярно-динамических симуляций. Описаны методы и параметры квантово-механических (КМ) расчетов энергии модельной молекулы, отвечающей основной части хромофора, содержащей тс-электроны ароматических колец. Описано получение зарядов, равновесных длин связей и значений валентных и двугранных углов модельной молекулы.

Описано получение профилей энергии валентного и двух двугранных углов вокруг специфического участка -С_Ш_СА1, не похожего ни на одну из описанных в поле ОРЬБ-АА молекул и/или частей молекул. Приведены полные таблицы параметров силового поля для описания хромофора для проведения МД-симуляций. Приведено сравнение структур белка ОзЯеё, полученных в результате МД-симуляций и данных рентгено-структурного анализа. Показано, что структуры совпадают с точностью до теплового движения молекул, т.е. топология хромофора верна. Проведен анализ изменения геометрии для всех вариантов мутаций хромофора по первой позиции, проведен поиск корреляций изменений геометрических параметров и оптических свойств.

В четвертой главе диссертации описана процедура создания методами генной инженерии двадцати плазмид на основе гена белка тКРР1, содержащих гены всех возможных вариантов флуоресцентных белков с заменами по 66 позиции. Описаны методы наработки и очистки целевого продукта белков.

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

Основные результаты работы

Расчет параметров топологи хромофоров

Хромофор белка DsRed рассматривается как уникальный аминокислотный остаток в составе полипептидной цепи. Хромофор обладает сильной делокализацией электронной плотности, что значительно отличает его от любого соединения, ранее описанного в силовом поле OPLS-AA. Это обусловило необходимость получения оригинальных параметров для описания хромофора.

Для проведения ab initio КМ-расчетов использовался программный пакет GAMESS. Процедура геометрической оптимизации структуры молекулы проводилась ab initio с использованием ограниченного метода Хартри-Фока RHF, базисом волновых функций 6-31G**. Расчет проводился с применением теории функционала плотности DFT с использованием трехпараметрического функционала Бека B3LYP.

Для нахождения значений эффективных зарядов атомов, образующих хромофор, была использована программа R.E.D. Эта программа позволяет находить такие эффективные значения точечных зарядов атомов молекулы в рамках МД-представлений, что электростатическое поле, порождаемое такими зарядами, аппроксимирует электростатическое поле молекулы, рассчитанное квантово-механически.

Строение нейтральной формы хромофора и его связь с белком DsRed показаны на рисунке 1. Каждому атому, входящему в состав хромофора, было присвоено уникальное имя. Имена атомов состоят из символа элемента, а также буквенно-цифрового индекса положения атома в молекуле.

В качестве первого приближения использовалась топология хромофора белка DsRed, порождаемая типами атомов, входящих в хромофор. Типы атомов подбирали по аналогии с соединениями, описанными в поле OPLS-AA.

Рис. 1. Строение нейтральной формы хромофора DsRed в составе белковой цепи. Хромофор очерчен пунктиром, Nur - N-концевой фрагмент белковой цепи, CKr - С-концевой фрагмент белковой цепи. Указаны имена атомов для создания топологии.

Для атомов CZ, ОН и НН хромофора (рис. 1) в качестве аналогов для задания типа были выбраны атом углерода С при гидроксильном радикале фенола, а также атомы кислорода и водорода гидроксильного радикала фенола из имеющегося в поле OPLS-AA описания молекулы фенола. Для атомов CEI, СЕ2, CD1, CD2, CG2, HEI, НЕ2, HD1, HD2 и атома СВ2 хромофора в качестве аналогов были выбраны атомы углерода и водорода из описания молекулы бензола в поле OPLS-AA. Для атомов СА2 и Cl, N2, N3 хромофора в качестве аналогов были выбраны, соответственно, атомы С4, N7, N9 из описания аденина в поле OPLS-AA. Атому N1 был приписан тип атома N3 из описания в силовом поле OPLS-AA цитозина. Атомам С2 и 02 были приписаны типы атомов С6 и связанного с ним кислорода О, соответственно, из описания гуанина.

Для атомов хромофора, входящих в состав бокового радикала, соответствующего боковому радикалу аминокислотного остатка Gln66 незрелого белка DsRed в качестве аналогов были выбраны соответствующие атомы из описания глутамина в силовом поле OPLS-AA. Для атома САЗ хромофора белка DsRed в качестве аналога был выбран Са-атом из описания глицина в поле OPLS-AA.

Данный набор параметров является первым приближением топологии нейтральной формы хромофора белка DsRed.

В конечный вариант топологии хромофора белка DsRed из данного начального набора вошли параметры атомов для расчета потенциала Леннарда-Джонса и жесткости связей, 10

валентных и двугранных углов (за исключением углов -C_N1_CA1, -CA_-C_N1_CA1 и C_N1_CA ICI, жесткости которых были рассчитаны отдельно).

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

.ОН _СЕ1 ,HD1 / \ / X / IH CZ

,СЕ2 ^CG2 ^.НВ2 ^С02 СВ2

I II

HD2 СА2

/ \ ^.02 N2 С2^

НВ12 / \ / ,С1-N3

/ \ ^НЛЗ'

исмо -

1 /

-о=с-

I

^-СА--НА

t /Л

-НА1 \

-НА2

Рис. 2. Модельная молекула А.

Параметры атомов и межатомных взаимодействий бокового радикала хромофора, соответствующего боковому радикалу аминокислотного остатка Gin66 незрелого белка DsRed задавались равными параметрам глутамина из силового поля OPLS-AA. Поэтому, для упрощения расчетов, в модельную молекулу А не были включены атомы: CGI, HG11, HG12, CD3, 0Е1, NE1, НЕ11, НЕ12 (рис.7), а к атому СВ1 для насыщения валентности был добавлен атом водорода HB 13 (рис. 2).

Дополнительно для расчета параметров валентного и двугранных углов при атоме N1 в модельную молекулу А был включен фрагмент предшествующего хромофору остатка фенилаланина Phe65 (нумерация по аминокислотной последовательности белка DsRed), включающий атомы: -С, -О, -CA и -НА. К атому -CA были добавлены атомы водорода -НА1 и -НА2 (рис. 2) для насыщения свободных валентностей.

Для определения значений эффективных зарядов атомов, равновесных значений длин связей, валентных и двугранных углов квантово-механически была проведена геометрическая оптимизация структуры модельной молекулы А. В качестве исходной структуры модельной

молекулы А была использована структура хромофора в структуре белка РвЯес! (РОВ Ю ЮОХ). Данная структура модельной молекулы А была названа А100Х. В результате расчета была получена равновесная структура модельной молекулы А, названная А4, отвечающая состоянию молекулы А в энергетическом минимуме. Полученные значения длин связей и углов структуры А64, были использованы в качестве значений равновесных длин связей Ьф валентных углов 8Ц1т и двугранных углов в окончательном описании топологии хромофора белка ОэКес! в поле

ОР1.8-АА.

Для структуры Асч квантово-механически были рассчитаны эффективные значения точечных зарядов на атомах при помощи программы Я.Е.О. Полученные значения зарядов использовались в качестве уточненных зарядов атомов хромофора и были внесены в окончательный вариант топологии хромофора белка ОзЯес! в поле ОРЬЭ-АА.

Для расчета параметров валентного угла -С_Ш_СА1 была использована структура А54. В ней изменяли равновесное значение валентного угла -С_>П_СА1 (0° = 123.8°) на величину -5° до +5° с шагом в 1", оставляя значения всех остальных углов и связей неизменными. Тем самым, был получен набор из десяти неравновесных структур, обозначенных В* (Д8).

Для каждой такой структуры квантово-механически была проведена процедура геометрической оптимизации с фиксированным значением валентного угла 0. В результате был получен набор из десяти новых равновесных структур, с фиксированным значением валентного угла 0, обозначенных как В'4 (Л0).

Все структуры вместе со структурой А84 имеют заданное отклонение Д9 валентного угла 0 от его равновесного значения в пределах -5° до +5°, включая 0° (равновесное значение для структуры ЛСС|) и геометрически оптимизированную структуру всех остальных частей молекулы. Для полученных структур была вычислена энергия (рис. 3).

Для нахождения параметров потенциалов двугранных углов -СА_-С_Ш_СА1 и -С_Ш_СА1_С1 в качестве исходной структуры модельной молекулы А бьиа использована равновесная структура АК|, в которой изменяли значение соответствующего угла, в результате чего был получен набор неравновесных структур, в которых углы принимал значения от 0° до 360° с шагом 20°, проведена процедура геометрической КМ-оптимизации с фиксированным значением соответствующего двугранного угла и вычислена энергия (рис. 4).

2,5 -,

2,0 -

1,0 -

<1 >

0,5 -

0,0 -

—-1-.-1---1-.-1-.-1-.-1

-6 -4 -2 0 2 4 6

А9, градусы

Рис. 3. Зависимость энергии от значения валентного угла -С_Ы1_СА1. Квадратами показаны значения энергии, полученной в результате квантово-механического расчета, сплошная линия -аппроксимация гармоническим потенциалом. 60

150 200 250

угол, градусы

300

350

Рис. 4. Зависимость энергии от значений валентный углов -СА_-С_Ы1_СА1 и С_Ы1_СА1_С1. Квадратами и кругами показаны значения энергии по результатам квантово-механического расчета для структур с фиксацией углов -СА_-С_Ы1_СА1 и С_М1_СА1_С1, соответственно, сплошная линия - аппроксимация расчетных данных периодическим гармоническим потенциалом.

Полученные наборы точек были аппроксимированы кривыми, описывающими периодический потенциал изменения величины двугранного угла. Параметры потенциалов С„ (где п=1,2,..., 6) для углов -CA_-C_N1_CA1 и -C_N1_CA1_C1 вошли в окончательный вариант топологии хромофора белка DsRed.

Пользуясь стандартной схемой проведения МД расчета с использованием силового поля OPLS-AA/DsRed, был проведен МД расчет для системы, содержащей тетрамер белка DsRed, 12 510 молекул воды и четыре иона натрия на протяжении 4.2 нс и получена усредненная по траектории структура хромофора, обозначенная QYGMD,Averase.

Сравнение геометрии хромофора QYGMD'Averai!C со структурами, полученными РСА (PDB Ш 1ZGO, 1G7K, 1GGX) показало, что с точностью до теплового движения все структуры идентичны, что подтверждает правильность параметризации хромофора.

Создание топологии хромофоров мутантов mRFPl

Поскольку на момент проведения работы отсутствовали данные о расшифровке структур mRFPl и мутантных белков методами РСА или ЯМР, для проведения молекулярно-динамических расчетов в качестве начальных использовались структуры, построенные на основе структуры ближайшего гомолога, белка DsRed, путем внесения соответствующих аминокислотных замен с помощью программы SwissPDBViewer. Полученные структуры были оптимизированы методом минимизации энергии по алгоритму L-BFGS в программном пакете GROMACS.

МД-расчеты проводились на основе силового поля OPLS-AA/DsRed. Для создания описаний хромофоров мутантных белков требовалось внести изменения, соответствующие модифицированному участку хромофора. Параметры были извлечены из записей соответствующих аминокислот.

Анализ полученных в результате молекулярно-динамических расчетов траекторий также осуществлялся средствами программного пакета GROMACS.

Создание генов и наработка целевого продукта мутантов белка mRFPl

В работе использовался ген красного мономерного флуоресцентного белка mRFPl из плазмидного вектора pMT-mRFPl, любезно предоставленного «Fungal Genetics Stock Center», Университет Миссури, США.

Основу олигонуклеотид-направленного сайт-специфического мутагенеза составил метод «ник-трансляции», при котором участок векторной ДНК тем или иным способом получают в однотяжевой форме, связывают его с синтетическим олигонуклеотидом, спланированным так, что 14 .

он комплементарен выбранному участку ДНК за исключением небольшого фрагмента в области внесения мутации. Полученный комплекс достраивают ДНК-полимеразой и замыкают ДНК-лигазой в двутяжевую ковалентно-замкнугую кольцевую векторную ДНК с небольшим некомплементарным участком в районе планируемой мутации. После трансформации подходящих компетентных клеток и клонирования проводят поиск и идентификацию целевых мутантов.

Для проведения генно-инженерных модификаций была осуществлена оригинальная реализация метода «ник-трансляции» (рис. 5), для проведения которой недалеко от места введения мутации должен существовать уникальный сайт узнавания эндонуклеазы рестрикции.

сайт рестрикции

рсстрикчэза SnaB [

чш

жп

II

РП~П~0~

отжиг праймеров Р^

IV

Исходная двутяжевая кольцевая плазмидная ДНК была расщеплена в уникальном сайте эндонуклеазой рестрикции с образованием линеаризованной двутяжевой плазмидной ДНК с короткими однотяжевыми участками на концах («липкими концами») (стадия I на рис. 5).

Затем плазмидная ДНК была обработана экзонуклеазой III, которая с обоих концов линеаризованного вектора высвобождает достаточно протяженные однотяжевые участки (стадия II на рис. 5).

На освобожденные участки комплементарно были отожжены синтетические олигонуклеотиды: мутаген и адаптер (стадия III на рис. 5). Олигонуклеотид-мутаген имел фосфатную группу на 5'-конце, в то время как синтетический олигонуклеотид-адаптер не имел такой группы. Структура З'-концевых участков олигонуклеотидов была подобрана таким образом,

15

чтобы после гибридизации и обработки конструкции ДНК-лигазой концы вектора комплементарно соединились с восстановлением исходного базового сайта эндонуклеазы рестрикции (стадия IV на рис. 5).

На следующем этапе полученный гибрид был достроен фрагментом Кленова ДНК-полимеразы I Е.соН в присутствии смеси дезоксирибонуклиозидтрифосфатов и ДНК-лигазы (стадия V на рис. 5). Основной продукт этой реакции - кольцевой гетеродуплекс ДНК с разрывом в одном из тяжей (в области гибридизации с дефосфоршшрованным адаптером) и нарушением двойной спирали вследствие некомплементарности в области мутагенеза, причем целевая «мугантная» форма гена встроена в ковалентно-замкнутый тяж ДНК, не имеющий разрыва, а исходная структура гена едикого типа» встроена в тяж с искусственным разрывом, расположенным в направлении 3' —» 5' на относительно небольшом расстоянии от места мутагенеза.

Обработка полученной гетеродуплексной конструкции ДНК-полимеразой I Е.соН привела к реакции так называемой «ник-трансляции» (стадия VI на рис. 5). В ходе этой реакции фермент осуществляет деградацию цепи с разрывом с З'-конца и одновременную достройку с 5'-конца на основе комплементарной цепи. В качестве «затравки» используется нефосфорилированный олигонуклеотид-адаптер. После прохождения полимеразой района некомплементарности образуется ДНК, обе цепи которой в районе мутагенеза будут иметь запланированную структуру.

Полученной кольцевой ковалентно-замкнутой молекулой ДНК трансформировали клетки Е.соН. Выделенная из полученной культуры клеток плазмидная ДНК была подвергнута рестриктному анализу, структура гена была проверена методом прямого секвенирования. Отбор клонов проводился с использованием анализа методом ПЦР, рестриктного анализа и секвенирования участка, содержащего предполагаемые мутации. Плазмидными ДНК отобранных клонов трансформировали компетентные клетки Е.соН, которые инокулировали в питательную среду БУТ. Из выросшей биомассы клеток-суперпродуцентов (ночной культуры) после лизиса методом никель-афинной хроматографии были выделены экспрессируемые белки. Пробы белков были подвергнуты диализу для очистки. Концентрации полученных препаратов белков составили от 0.3 до 1.5 мг/мл.

Оптические свойства полученных мутантных белков

Для всех полученных препаратов белков были измерены спектры поглощения, возбуждения и испускания флуоресценции.

Двенадцать из полученных 20 белков (60%) обладают окраской и флуоресценцией (табл. /). Одиннадцать из них флуоресцируют в красной области, мутант с заменой на глицин обладает 16

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

№ Кодовое имя Замена Название аминокислоты Максимум возбуждения, нм Максимум ЭМИССИИ, 1IM V с, М"'см"'

1 плт-пт нет глутамин 584 607 0.27 41 800

2 А5 N аспарагин 570 604 0.17 6 900

3 41 А аланин 578 605 0.19 2 400

4 70 L лейцин 577 613 0.12 3 600

5 11т С цистеин 568 588 0.21 35 650

6 2258 S серии 561 579 0.21 32 700

7 А8 H гистидин 588 618 0.13 12 700

8 160 G глицин 504 515 0.025 16 200

9 14Э Е глутамат 577-588 610 0.34 9 000

10 13шп M метионин не окрашен не флуоресцирует

11 12тп К лизин не окрашен не флуоресцирует

12 51 I изолейцин не окрашен не флуоресцирует

13 11т Т треонин 560-570 589 0.024 25 000

14 5Р F фенилаланин 595 624 0.03 24 350

15 ЮУДО V валин не окрашен не флуоресцирует

16 3911 R аргинин не окрашен не флуоресцирует

17 13Р Р пролин не окрашен не флуоресцирует

18 19У Y тирозин 595 624 0.04 2 430

19 1Ш D аспартат не окрашен не флуоресцирует

20 зш W триптофан не окрашен не флуоресцирует

Таблица 1. Оптические характеристики тЕРР и мутантных белков с заменами в 66 положении аминокислотной последовательности. Приведены кодовые имена мутантов и соответствующие им замены, 4/ - квантовый выход флуоресценции, е - эффективный коэффициент экстинкции.

Анализ спектров поглощения и возбуждения

У белка тКРР1 и его мутантов с аминокислотными заменами 066А, (366Ь, 0665, С}66Н,

С>66С, С!66К в спектрах поглощения присутствуют две полосы поглощения: зеленая и красная

17

(табл. 2). Для всех белков, кроме белка С>66М, значения длин волн, соответствующих максимумам зеленой полосы, лежат в диапазоне 500-509 нм, в то время как для белка максимум

поглощения зеленой полосы находится при 525 нм. Эту полосу можно отнести к поглощению не полностью созревшей зеленой формы хромофора. Значения длин волн, соответствующих максимумам красной полосы, лежат в диапазоне 561-598 нм и различны для различных замен в 66 положении.

Зеленая форма Красная форма

Белок Максимум поглощения, нм М'см1 Максимум поглощения, нм М'см'1 Максимум возбуждения, нм Максимум эмиссии, мм V

mRFPl 503 34 600 584 41 800 584 607 0.27

Q66N 525 13 700 570 6 900 570 604 0.17

Q66H 504 35 600 588 12 700 588 618 0.13

Q66L 500 64 000 582 3 600 577 613 0.12

Q66C 505 36 400 568 35 650 568 588 0.21

Q66A 500 29 600 582 2 400 578 605 0.19

Q66S 506 24 900 562 32 700 561 579 0.21

Q66F 508 26 600 593 8 000 593 617 0.03

Q66T 509 плечо 561 25 000 560 - 570 589 0.024

Q66Y 508 1200 592 2 400 595 624 0.04

Q66E 503 57 000 588 9000 577-588 610 0.34

Таблица 2. Оптические свойства зеленых и красных форм хромофоров белка т№Р1 и мутантных белков, е - эффективный коэффициент экстинкции в максимуме поглощения, у> -квантовый выход флуоресценции.

Определены эффективные коэффициенты экстинкции и квантовые выходы флуоресценции при длинах волн, соответствующих максимумам зеленых и красных полос (табл. 2).

Положение максимумов поглощения зеленой полосы мутантов практически не зависит от заменяемого остатка, напротив, положения максимумов поглощения красной полосы мутантов, а также максимумов возбуждения и эмиссии флуоресценции зависит от заменяемого остатка.

Экспериментальные данные не противоречат представлениям, согласно которым в полученных препаратах имеется две формы созревших флуоресцентных белков - красная и зеленая.

Если предположить, что зеленая форма хромофора отличается от красной отсутствием

двойной связи между Са и N атомами 66 аминокислотного остатка, то собственно система сопряженных я-связей зеленой формы хромофора не будет в себя включать 66 аминокислотный остаток и можно предположить, что его замена не должна значительно влиять на спектральные свойства. В случае же «красного» хромофора мутантов дегидрирование связи между Са и N атомами вовлекает углеродный скелет 66 аминокислотного остатка в систему сопряжения и изменения в боковом радикале этого остатка могут приводить к «тонкой» настройке спектров поглощения и флуоресценции.

Выявлены корреляции между спектральными свойствами полученных белков (длинами волн максимумов возбуждения и эмиссии флуоресценции) и геометрическими параметрами мутируемых аминокислотных остатков (площадь и объем). С увеличением размеров остатка сдвигаются в более красную область положения максимумов возбуждения и эмиссии флуоресценции (рис. 6).

Для набора мутантов С?66А, (366Ц (}66$, С266Н, 066С, <266М, 066Т, 066У, (266Р, 066Е положения максимумов возбуждения флуоресценции коррелируют с объемом остатка по 66 положению с коэффициентом корреляции 0.66, положения максимумов эмиссии коррелируют с объемом со значением 0.63. Однако если из набора мутантов исключить неполярные аминокислотные остатки (аланин, лейцин, тирозин и фенилаланин), то в оставшемся ряду полярных аминокислотных остатков степени корреляций возрастают до 0.92 и 0.83, соответственно Примечательно, что в ряду неполярных замен коэффициенты корреляции низки -0.28 и 0.69, соответственно. Данное расщепление можно объяснить наличием двух факторов, оказывающих влияние на спектральные свойства:

1. геометрический размер заменяемой аминокислоты (приводящий к отдалению от хромофора части аминокислот окружения с увеличением размера остатка)

2. в зависимости от типа аминокислоты по 66 положению возможно образование водородных связей полярными аминокислотами с остатком С1и215, а так же изменение энергии системы в случае замены на алифатические аминокислотные остатки за счет гидрофобных взаимодействий.

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

о&ьем, А1

Рис. 6. Зависимости положений максимумов возбуждения и эмиссии флуоресценции от объемов боковых радикалов аминокислотных остатков, находящихся в первом положении хромофора для белка тКРР1 и его мутантов с заменами в 66 положении. Квадратами обозначены значения максимума возбуждения, кругами -максимума эмиссии.

объем, А3

Рис. 7, Зависимость положений максимумов возбуждения (квадраты) и эмиссии (круги) флуоресценции от объемов боковых радикалов аминокислотных остатков для белка йРР и его мутантов с заменами в 65 положении.

Анализ МД-структур мутантных белков и т!*РР1

Структуры, полученные в результате молекулярно-динамических расчетов, сохранили общую форму, характерную для флуоресцентных белков (р-бочонок, внутри которого проходит а-спираль с хромофором в ее центральной части). Хромофоры остались плоскими и при наложении хорошо совмещаются (см. рис. 8).

Рис. 8. Сравнение 5 А окружения хромофоров белков тККР1, Q66N, Q66C Q66A и <2665.

Анализ полученных в результате молекулярно-динамических симуляций средних по траектории структур белков показывает, что полярные боковые радикалы аминокислот в 66 положении могут образовывать водородные связи с глутамином 215. Образование выгодной по энергии водородной связи требует сближения соответствующих атомов. Однако такое сближение не может происходить за счет изменения положения 01и215, т.к. он является частью жесткого |3-бочонка, а нарушение конфигурации бочонка требует больших энергетических затрат. Хромофор является гораздо более подвижной структурой, поэтому осуществить образование необходимой водородной связи за счет изменения конформации хромофора энергетически выгоднее. Как показали результаты молекулярно-динамических расчетов, у разных мутантных белков

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

Определение оптических характеристик методом нелинейной флуориметрии

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

Определение соотношения концентраций .между зеленой и красной формами белков

Используя совместно метод нелинейной флуориметрии и классические методы, можно определить соотношение концентраций между флуоресцирующим и нефлуоресцирующим хромопротеинами красного флуоресцентного белка в растворе, решая следующую систему уравнений:

Г

Фо(Яех=570) «{Г- »а0

<532) X „ .. ~<532) = 2.3 • Б'532'/1

'^гееп ^гееп пге<1 ^геен = п0

индексы 532, 570 и red, green обозначают значения длины волны возбуждающего излучения и формы белка, соответственно;

где nred, ngreen - концентрации красной и зеленой форм белка в растворе [см"3], по - общая концентрация молекул белка, о - сечение поглощения хромофора, <trs - сечение комбинационного рассеяния (КР) воды, D - оптическая плотность раствора белка, Фо - флуоресцентный параметр в отсутствии насыщения флуоресценции, 1 - толщина слоя раствора при абсорбционных измерениях, также предполагается, что квантовый выход г) не зависит от длины волны возбуждения.

Значение длины волны возбуждения 570 нм было выбрано так, чтобы оптическая плотность D, измеренная с помощью спектрофотометра, определялась поглощением только флуоресцирующей красной формой белка, т.е. чтобы для всех образцов было выполнено D = Dred.

Полученное значение сечения позволяет определить индивидуальные значения

сечения поглощения а [см2] и коэффициента экстинкции г = 2.6Ю20а [М"'см~'] хромофора в области длин волн, где измеренные величины оптической плотности D определяются только поглощением красной формы белка, в частности, коэффициент экстинкции в максимуме полосы поглощения (см. таблицу 3).

Метод был опробован на примере трех белков, имеющих явно выраженные красную и зеленую полосы в спектрах поглощения: mRFPl, Q66S и Q66C (см. таблицу 3).

mRFPl Q66C Q66S

a«fd2), 10"16см2 3.4±0.6 3±0.5 2.2+0.5

Ю"'7см2 2.4±0.4 1.6+0.4 2.3+0.4

of™\ 10"16см2 6.5+1.1 5.1+0.8 3±0.5

мМ-'см"1 215+40 135+20 85+13

nred / По 0.26±0.06 0.17±0.06 0.34±0.06

ngreen / По 0.74+0.06 0.8310.06 0.66+0.06

nred / ngreen 1:3 1:5 1:2

Тз, НС 3+0.15 2.8+0.14 2.9+0.14

Г|Т 0±0.02 0±0.02 0.05+0.02

П 0.24±0.03 0.19±0.04 0.20+0.04

Таблица 3. Индивидуальные оптические характеристики белков mRFPl, Q66C и Q66S.

В максимуме полосы поглощения красной формы белка шЯРР1 (к ~ 584 нм): сечение поглощения хромофора о"^^ составляет (8.2+1.5)-10'16 [см2], что соответствует коэффициенту экстинкции равному 215+40 [мМ"'см"']. Полученное значение в 4 раза выше определенного ранее другими авторами. Такое различие связанно с тем, что в этих работах при расчете г использовали суммарную концентрацию белка (и поэтому определяли интегральный коэффициент экстинкции), а не парциальную. Как следствие, определение парциальной концентрации хромофоров позволило найти индивидуальные коэффициенты экстинкции двух форм белков. В общем случае, определение оптических характеристик ФБ на основе интегральных характеристик препарата является некорректным. Предложенный метод позволяет определять истинные характеристики флуоресцентных белков.

Как было указано ранее, положение максимумов поглощения зеленой полосы мутантов практически не зависит от типа заменяемого остатка. Напротив, положения максимумов поглощения красной полосы мутантов, а также максимумов возбуждения и эмиссии флуоресценции зависят от заменяемого остатка и хорошо коррелируют с объемом бокового радикала 66-ого аминокислотного остатка (рис. 6), причем с увеличением объема положения максимумов сдвигаются в длинноволновую область спектра. Аналогичная корреляция была выявлена и для истинных коэффициентов экстинкции (рис. 9).

80 90 100 110 120 130

0 3

объем,А

Рис. 9. Зависимость истинного коэффициента экстинкции от объемов аминокислотного остатка по 66 положению.

Результаты и выводы

1. Методами квантовой химии и молекулярной динамики в силовом поле ОРЦЗ-АА создана топология хромофора белка шЯРР1 как нового аминокислотного остатка и рассчитаны ее параметры. Проведенное молекулярно-динамическое моделирование структуры белка ОзЯе<1 (с аналогичным хромофором) и сравнение с имеющимися рентгено-структурными данными показало правильность созданной топологии.

2. Методами квантовой химии и молекулярной динамики в силовом поле OPL.S-.AA созданы топологии и рассчитаны параметры для вссх 20 возможных хромофоров с заменами по 66 позиции.

3. Созданы РБВ-файлы структур белка т1ШЧ и его мутантов.

4. Проведено молекулярно-динамическое моделирование пространственных структур белков в растворе. Рассчитаны равновесные структуры белков. Наложение средних по траектории структур показало, что значительных изменений в конформации хромофора и его окружения не произошло, однако имеются некоторые геометрические отличия, которые, видимо, отвечают за сдвиги в спектрах поглощения и флуоресценции у мутантов. В результате анализа структур хромофоров мутантов выяснено, что наиболее вариабельными участками являются область соединения фенольного и имидозалидонового колец хромофора, а также зона присоединения мутированного радикала.

5. Методом олигонуклеотид-направленного сайт-специфического мутагенеза гена мономерного красного флуоресцентного белка (тКРР1), клонированного в модифицированной плазмиде р(2ЕМ, получен набор мугантных генов с заменами кодона 01п66. Экспрессия в Е.соИ и последующая аффинная очистка позволили получить серию препаратов из 20 мугантных белков тЮТЧ. Проведена наработка и выделение целевого белка. Проведенный электрофорез проб показал наличие во всех белка с массой 26 кДа, что соответствует целевому продукту. Измерены оптические свойства препаратов. Из 20 белков окраской и способностью флуоресцировать в видимой области обладают 12. Выявлена корреляция между объемом бокового радикала а.о. по 66 положению и оптическими свойствами полученного белка.

6. Определены индивидуальные оптические характеристики белка тМР 1 и его мутантов ((}66С и С?668), что позволило определить концентрации «зеленой» и «красной» форм хромофоров в растворе. Выявлена зависимость индивидуального коэффициента экстинкции красной формы белка (в максимуме полосы поглощения) от объема аминокислоты по 66 положению, что может быть применено для прогнозирования свойств новых мутантов флуоресцентных белков и синтеза белков с наперед заданными свойствами.

Список работ, опубликованных по теме диссертации

Статьи

1. Банишев АА, Маслов ДВ, Ширшин ЕА, Фадеев ВВ, Вржещ ЕП, Дмитриеико ДВ, Друца ВЛ, Вржещ ПВ, Пащенко ВЗ. Метод определения индивидуальных оптических характеристик посттрансляционных флуоресцентных форм флуоресцентных белков с использованием нелинейной лазерной флуориметрии. Биофизика, 2007,52(5), 792-798.

2. Khrameeva ЕЕ, Drutsa VL, Vrzheshch ЕР, Dmitrienko DV, Vrzheshch PV. Mutants of monomeric redfluorescent protein mRFPl at residue 66: structure modeling by molecular dynamics and search for correlations with spectral properties Biochemistry, 2008, 73(10), 1085-1095.

3. Dmitrienko DV, Vrzheshch EP, Drutsa VL, Vrzheshch PV. Red fluorescent protein DsRed: parametrization of its chromophore as an amino acid residue for computer modeling in the OPLS-AA force field. Biochemistry, 2006, 71(10), 1133-1152.

4. Вржещ ЕП, Дмитриеико ДВ, Руданов ГС, Загидуллин ВЭ, Пащенко ВЗ, Разживин АП, Салецкий AM, Вржещ ПВ. Оптические свойства мономерного красного флуоресцентного белка mRFPl. Вестник Московского Университета, серия биология, 2008,3,21-24.

5. Banishev АА, Vrzheshch ЕР, Shirshin ЕА. Application of laser fluorimetry for determining the influence of a single amino-acid substitution on the individual photophysical parameters of a fluorescent form of a fluorescent protein mRFPl. Quantum Electronics, 2009, 39(3), 273-278.

Тезисы докладов на конференциях

1. Ширшин ЕА, Банишев АА, Вржещ ЕП. Фотофизические параметры и механизмы флуоресценции белка. Международная научная конференция студентов, аспирантов и молодых ученых «Ломоносов-2006», секция «Физика», Москва, 2006, Тезисы докладов, с. 157.

2. Фадеев ВВ, Вржещ ПВ, Маслов ДВ, Банишев АА, Вржещ ЕП, Дмитриенко ДВ, Ширшин ЕА. Флуоресцентные характеристики и молекулярные фотофизические параметры красного флуоресцентного белка mRFPl. Труды конференции «Альманах клинической медицины», издается МОНИКИ им. М.Ф. Владимирского, Москва, 2006. Том XII, с. 42.

3. Вржещ ЕП, Дмитриенко ДВ, Вржещ ПВ. Компьютерное моделирование трехмерных структур флуоресцентных белков. Труды конференции «Альманах клинической медицины», издается МОНИКИ им. М.Ф. Владимирского, Москва, 2006. Том XII, с. 51.

4. Цой О, Вржещ ЕП. Изучение физико-химических свойств мутантов красного флуоресцентного белка mRFPl. XIII международная конференция студентов, аспирантов и молодых ученых «Ломоносов-2006», Секция «Биоинженерия и биоинформатика», Москва, 2006, Тезисы докладов,

с. 10.

Заказ № 72-П/02/2012 Подписано в печать 10.02.2012 Тираж 110 экз. Усл. п.л.1,3

ООО "Цифровичок", тел. (495) 649-83-30 www.cfr.ru; e-mail:info@cfr.ru

Текст научной работыДиссертация по биологии, кандидата физико-математических наук, Вржещ, Евгений Петрович, Москва

61 12-1/586

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ИМЕНИ М.В. ЛОМОНОСОВА

ФИЗИЧЕСКИЙ ФАКУЛЬТЕТ КАФЕДРА ОБЩЕЙ ФИЗИКИ

На правах рукописи

Вржещ Евгений Петрович

Влияние аминокислотных замен на спектральные свойства флуоресцентных белков на примере

тЖ1

03.01.02 - биофизика Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель: д.ф-м.н., проф. Салецкий А.М.

МОСКВА 2012

Содержание

Содержание.................................................................................................................................................2

Список сокращений..................................................................................................................................6

Введение.......................................................................................................................................................8

Задачи исследования...............................................................................................................................10

1. Обзор литературы................................................................................................................................11

1.1. Флуоресцентные белки, природные источники и применение..............................................11

1.2. Биологическая роль ФБ..............................................................................................................16

1.3. Кристаллическая структура ФБ................................................................................................17

1.3.1. Кристаллическая структура вРР........................................................................................17

1.3.2. Кристаллическая структура ОзЯес!....................................................................................19

1.4. Формирование хромофоров белков..........................................................................................21

1.4.1. Формирование хромофора вРР..........................................................................................21

1.4.2. Формирование хромофора БэКес!......................................................................................23

1.4.3. Влияние окружения хромофора на созревание................................................................27

1.5. Олигомеризация флуоресцентных белков...............................................................................29

1.6. Получение тКРР1 - мономера БвКес!......................................................................................30

1.7. Физические характеристики белков ОБР и DsR.ec!..................................................................32

1.8. Спектральные характеристики ФБ...........................................................................................32

1.8.1. Спектральные характеристики вБР...................................................................................33

1.8.2. Спектральные характеристики DsR.ec!...............................................................................33

1.8.3. Спектральные характеристики пгЕШЧ..............................................................................34

1.8.4. Влияние а.о. в хромофоре на флуоресцентные свойства белка......................................35

1.9. Заключение..................................................................................................................................38

2. Молекулярная динамика...................................................................................................................39

2.1. Методы........................................................................................................................................39

2.1.1. Метод молекулярной динамики.........................................................................................39

2.1.2. Силовое поле ОРЬ8-АА......................................................................................................41

2.1.3. Программный пакет ОЮМАСЗ........................................................................................43

2.2. Схема проведения МД-расчетов...............................................................................................44

2.2.1. Исходная структура белка ВвЛес!......................................................................................44

2.2.2. Задание системы для МД-расчетов....................................................................................44

2.2.3. Минимизации энергии системы.......................................................-..................................44

2.2.4. Проведение МД-расчета с ограничением подвижности атомов белка..........................45

2.2.5. Получение МД-траектории системы.................................................................................46

2.3. Анализ МД-траектории..............................................................................................................47

2.3.1. Расчет RMSD атомов мономера DsRed.............................................................................47

2.3.2. Получение структуры хромофора, усредненной по траектории....................................47

2.3.3. Расчет RMSF атомов в МД-траектории и визуализация их подвижности.....................47

2.3.4. Получение распределений значений валентных и двугранных углов...........................48

2.4. Результаты...................................................................................................................................48

3. Квантовая механика...........................................................................................................................49

3.1. Методы........................................................................................................................................49

3.1.1. Метод квантовой механики................................................................................................49

3.1.2. Программный пакет GAMESS...........................................................................................49

3.1.3. Процедура геометрической оптимизации.........................................................................50

3.1.4. Расчет эффективных зарядов атомов хромофора.............................................................50

3.2. Создание топологии хромофора белка DsRed.........................................................................50

3.2.1. Первое приближение топологии хромофора DsRed........................................................52

3.2.2. Уточнение топологии хромофора белка DsRed................................................................53

3.2.2.1. Модельная молекула для КМ-расчетов......................................................................53

3.2.2.2. КМ-расчет равновесной структуры модельной молекулы.......................................54

3.2.2.3. Уточненные значения эффективных точечных зарядов...........................................54

3.2.2.4. Уточненные параметры связей....................................................................................58

3.2.2.5. Уточненные параметры валентных углов..................................................................59

3.2.2.6. Уточненные параметры двугранных углов................................................................62

3.2.2.7. Параметры валентного угла -CN1CA1...................................................................65

3.2.2.8. Параметры двугранных углов -CA-CN1CA1 и C_N1_CA1_C1.........................66

3.2.3. Окончательный вариант топологии, расширение силового поля OPLS-AA.................68

3.2.4. МД-траектория тетрамера белка DsRed в силовом поле OPLS-AA/DsRed...................68

3.2.5. Структура хромофора белка DsRed, усредненная по МД-траектории...........................69

3.2.5.1. Значения углов -C_N1_CA1, -C_N1_CA1_C1 и -CA1_-C_N1_CA1........................73

3.2.5.2. Средние значения и дисперсии для двугранных углов.............................................75

3.3. Результаты...................................................................................................................................77

3.3.1. Сравнение количественных характеристик структур хромофора..................................77

3.3.2. Анализ значений углов -C_N1_CA1, -C_N1_CA1_C1 и -CA1-CN1CA1..................79

3.3.3. Анализ значений двугранных углов хромофора..............................................................79

3.3.4. Создание топологии хромофоров мутантов белка mRFPl..............................................80

3.3.4.1. Моделирование структур мутантных белков и тКРР1............................................80

3.3.4.2. Структуры мутантных белков и тМЛ.....................................................................80

3.3.4.3. Поиск корреляций.........................................................................................................81

4. Генная инженерия...............................................................................................................................85

4.1. Материалы...................................................................................................................................85

4.1.1. Реактивы...............................................................................................................................85

4.1.2. Ферменты.............................................................................................................................85

4.1.3. Генно-инженерные конструкции.......................................................................................85

4.1.4. Буферные растворы.............................................................................................................86

4.1.5. Приготовление культуральных сред..................................................................................87

4.2. Методы........................................................................................................................................87

4.2.1. Создание экспрессионного вектора рРЕбО-тЯБРЬБ......................................................87

4.2.1.1. Переклонирование гена белка тШП в экспрессионный вектор р(}Е-60..............87

4.2.1.2. Создание вспомогательного сайта рестрикции 8паВ1..............................................88

4.2.2. Олигонуклеотид-направленный сайт-специфический мутагенез ДНК.........................88

4.2.2.1. Расщепление исходной кольцевой ДНК в сайте рестрикции 8паВ1........................90

4.2.2.2. Осаждение ДНК............................................................................................................90

4.2.2.3. Обработка линеаризованной плазмиды экзонуклеазой III.......................................91

4.2.2.4. Проверка степени гидролиза экзонуклеазой III.........................................................91

4.2.2.5. Внесение мутаций в ДНК гена тКРР1, кодирующую хромофор............................92

4.2.2.6. Фосфорилирование 5'-концов олигонуклеотидов-мутагенов (кинирование)........92

4.2.2.7. Получение гетеродуплексной ДНК............................................................................92

4.2.2.8. Проведение «ник-трансляции»....................................................................................93

4.2.3. Трансформация и клонирование........................................................................................93

4.2.4. Отбор клонов, содержащих целевые мутантные гены.....................................................93

4.2.4.1. Анализ методом ПЦР...................................................................................................94

4.2.4.2. Рестриктный анализ......................................................................................................94

4.2.4.3. Секвенирование по методу Сенгера...........................................................................94

4.2.5. Экспрессия белка.................................................................................................................95

4.2.5.1. Лизис клеток для выделения белка.............................................................................96

4.2.5.2. Аффинное выделение белка на колонке с №-ЭТА агарозой...................................96

4.2.5.3. Диализ белков...............................................................................................................96

4.2.5.4. Электрофорез белков....................................................................................................97

4.3. Результаты...................................................................................................................................98

5. Оптика...................................................................................................................................................99

5.1. Материалы...................................................................................................................................99

5.2. Методы .........................................................................................................................................99

5.2.1. Классические спектральные измерения............................................................................99

5.2.2. Временные измерения.........................................................................................................99

5.2.3. Лазерные спектральные измерения.................................................................................100

5.2.4. Метод нелинейной лазерной флуориметрии..................................................................100

5.2.5. Регистрация спектров оптической плотности................................................................102

5.2.6. Определение коэффициентов экстинкции......................................................................103

5.2.7. Регистрация спектров флуоресценции............................................................................103

5.2.8. Определение квантового выхода.....................................................................................104

5.3. Результаты.................................................................................................................................105

5.3.1. Оптические свойства белков............................................................................................105

5.3.2. Анализ спектров поглощения и возбуждения................................................................111

5.3.3. Определение оптических характеристик методом нелинейной флуориметрии.........113

5.3.4. Определение соотношения концентраций между зеленой и красной форм белков... 115

5.3.5. Обсуждение результатов...................................................................................................117

Основные результаты и выводы работы..........................................................................................120

Список литературы...............................................................................................................................121

Список сокращений

amajGFP

amFP486

asCP

ATP

BamHI

dATP

dCTP

dGTP

dNTP

DsRed

dTTP

EcoRI

EDTA

EGFP

eqFP611

EYFP

GAMESS

GFP GFPs

GROMACS

HcRed

Hindlll

Kaede

KFP

KikG mRFPl

Ncol

OPLS-AA

зеленый флуоресцентный белок из анемоны Anemonia majano циановый флуоресцентный белок из анемоны Anemonia majano красный флуоресцентный белок из актинии Anemonia sulcata АТФ, аденозинтрифосфат, adenosine-5 '-triphosphate эндонуклеаза рестрикции, сайт GAGATCC дАТФ, дезоксиаденинтрифосфат, deoxyadenosine triphosphate дСТФ, дезоксицитозинтрифосфат, deoxycytidine triphosphate дГТФ, дезоксигуанинтрифосфат, deoxyguanosine triphosphate дНТФ, дезоксинуклеозидтрифосфат, deoxyribonucleotide triphosphate (drFP583) красный флуоресцентный белок из коралла Discosoma striata дТТФ, дезокситиминтрифосфат, deoxythymidine triphosphate эндонуклеаза рестрикции, сайт GAAATTC

ЭДТУ, ethylenediaminetetraacetic acid, этилендиаминтетрауксусная кислота зеленый флуоресцентный белок, мутант GFP

красный флуоресцентный белок из анемоны Entacmaea quadricolor

желтый флуоресцентный белок, мутант GFP

General Atomic and Molecular Electronic Structure System,

http://www.msg.chem.iastate.edu/gamess

зеленый флуоресцентный белок из медузы Aequorea victoria

зеленые флуоресцентные белки (максимум поглощения 485-520 нм)

GROningen MAchine for Chemical Simulations, http://www.gromacs.org

красный флуоресцентный белок из анемоны Heteractis crispa

эндонуклеаза рестрикции, сайт AAAGCTT

зеленый флуоресцентный белок из коралла Trachyphyllia geoffroyi красный «разжигаемый» флуоресцентный белок, полученный с помощью сайт-направленного и случайного мутагенеза хромобелка asFP595 из актинии Anemonia sulcata

зеленый флуоресцентный белок из коралла Favia favus

мономерный красный флуоресцентный белок, полученный на основе DsRed путем направленного мутагенеза эндонуклеаза рестрикции, сайт С ACATGG

optimized potentials for liquid simulations all atom (силовое поле пакета GROMACS)

PDB PSA R.E.D. RFPs RMSD RMSF Rtms5 SDS SnaBI spdbv

ТАЕ-буфер «трис-ацетат буфер с EDTA» для электрофореза ДНК

TEMED tetramethylethylenediamine, тетраметилэтилендиамин

YFPs желтые флуоресцентные белки (максимум поглощения 540 нм)

zFP506 зеленый флуоресцентный белок из полипа Zoanthus sp.

zFP538 желтый флуоресцентный белок из полипа Zoanthus sp.

a.o. аминокислотный остаток

БСА BSA, бычий сывороточный альбумин, bovine serum albnumin

Да Da, Дальтон, а.е.м., атомная единица массы, Dalton

кДНК комплементарная ДНК

KM квантовая механика

KP комбинационное рассеяние

МД молекулярная динамика

ПААГ PAGE, полиакриламидный гель, polyacrylamide gel electrophoresis

ПЦР PCR, полимеразная цепная реакция, polymerase chain reaction

ПЭГ PEG, полиэтиленгликоль, polyethylene glycol

РСА X-ray, рентгено-структурный анализ

СОС сложное органическое соединение

ТЕ-буфер «трис-EDTA» буфер для работы с ДНК

ФБ FP, флуоресцентный белок, fluorescent protein

ЯМР NMR, ядерно-магнитный резонанс, nuclear magnetic resonance

Protein Data Bank, http://www.rcsb.org

ПСА, персульфат аммония, ammonium persulfate

RESP ESP charge Derive, http://www.u-picardie.ir/labo/lbpd/RED

оранжево-красные флуоресцентные белки (максимум поглощения > 570 нм)

root mean square deviation, среднеквадратичное отклонение

root mean square fluctuations, среднеквадратичные флуктуации

нефлуоресцентный хромопротеин из коралла Montipora efflorescens

sodium dodecyl sulfate, додецилсульфат натрия, лаурилсульфат натрия

эндонуклеаза рестрикции, сайт TACAGTA

Swiss-PDB Viewer Deep View, http://spdbv.vital-it.ch

Введение

Семейство флуоресцентных белков включает в себя способные испускать и поглощать излучение в видимой части спектра белки, обнаруженные у различных представителей типа Кишечнополостные (Сое1еп1ега1а). Флуоресцентные белки имеют сходную структуру, представленную 11-тяжевым (3-бочонком, внутри которого проходит а-спираль с хромофором в ее центральной части (рис. 1):

«бочонка», фиолетовым — участки а-спиралей. В центре бочонка располагается хромофор (выделен сипим цветом). Серым обозначены петлевые участки.

Хромофор - это гетерогруппа, образующаяся в результате авто каталитической посттрансляционной реакции между тремя смежными аминокислотами. Его наличие обеспечивает окраску белков и возможность флуоресцировать (нефлуоресцирующие формы носят название хромопротеинов). В настоящее время описаны представители данного семейства с различными свойствами, такими как быстрота созревания, стабильность, оптические свойства (спектры поглощения и флуоресценции; квантовый выход флуоресценции, фотовыцветание и пр.).

Долгие годы флуоресцентные белки изучались лишь узким кругом ученых как компоненты

биолюминесцентной системы [2]. Решающий прорыв в их практическом применении произошел

после того, как в 1992 году клонировали зеленый белок ОРР [3], [4] и было показано, что

8

экспрессия его гена в других организмах также приводит к синтезу флуоресцентного белка [2], [4], [5], [6]. Следовательно, ген содержит всю информацию, необходимую для посттрансляционного синтеза хромофора, и в этом процессе не задействованы какие-либо специфические ферменты [2].

После этого открытия интерес к GFP чрезвычайно возрос, и в настоящее время GFP и его мутанты, а также другие флуоресцентные белки являются наиболее широко используемыми маркерами [2].

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

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

Часто для решения научных задач требуется помечать разные клеточные структуры маркерами разных цветов, что с легкостью осуществимо при помощи флуоресцентных белков. Но для этого необходимо иметь белки, флуоресценция которых приходится на различные участки спектра. Поэтому в последние годы ведется активное создание новых мутантных форм с измененными свойствами.

В 1999 году целый ряд окрашенных белков, включающий красный белок DsRed,