Влияние структурных свойств глобулярных белков на их механическую стабильность и термостабильность

Глякина, Анна Владимировна

Бесплатный автореферат и диссертация по биологии на тему
Влияние структурных свойств глобулярных белков на их механическую стабильность и термостабильность
ВАК РФ 03.01.02, Биофизика

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

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

бу

Глякнна Анна Владимировна

ВЛИЯНИЕ СТРУКТУРНЫХ СВОЙСТВ ГЛОБУЛЯРНЫХ БЕЛКОВ НА ИХ МЕХАНИЧЕСКУЮ СТАБИЛЬНОСТЬ И ТЕРМОСТАБИЛЬНОСТЬ

03.01.02 - биофизика

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

З О МАЙ 2013

Москва-2013

005060549

Работа выполнена в лаборатории молекулярной динамики Федерального государственного бюджетного учреждения науки Института математических проблем биологии Российской академии наук

Научные руководители: кандидат физико-математических наук,

доцент

Балабаев Николай Кириллович

доктор физико-математических наук Галзитская Оксана Валериановна

Официальные оппоненты: доктор физико-математических наук,

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

Шайтан Константин Вольдемарович

доктор биологических наук, ведущий научный сотрудник лаборатории молекулярной физиологии клетки ФГБУН Института биофизики клетки РАН Сивожелезов Виктор Семенович

Ведущая организация: Федеральное государственное бюджетное

учреждение науки Институт химической физики им. H.H. Семенова Российской академии наук

Защита состоится «13» июня 2013 г. в 14 час. 00 мин. на заседании Диссертационного совета Д 501.001.96 при Московском государственном университете имени М.В. Ломоносова по адресу: 119234, г. Москва, Ленинские горы, д. 1, стр. 12, биологический факультет, аудитория 389.

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

Автореферат разослан: " 6 " MQ.3 2013 г.

Ученый секретарь Диссертационного совета доктор биологических наук

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

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

Разворачивание отдельных белков растяжением их за концы стало возможным благодаря развитию атомно-силовой микроскопии. В конце 1986 года Г. Биннингом был разработан атомно-силовой микроскоп, открывший новые перспективы в изучении процесса сворачивания белков. Эксперименты по разворачиванию белков с помощью атомно-силового микроскопа обычно проводятся на многодоменных белках, либо для цепочки «сшитых» одинаковых или разных белков.

Компьютерные методы моделирования молекулярной динамики позволяют на атомном уровне исследовать процессы механического

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

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

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

Для достижения данной цели были поставлены следующие задачи:

1. Провести моделирование методом молекулярной динамики с использованием явной модели растворителя процесс силового разворачивания похожих по структуре, но различающихся по аминокислотной последовательности иммуноглобулинсвязывающих доменов белков Ь и в (эти белки являются хорошими модельными объектами и достаточно широко изучены экспериментально различными методами).

2. Выявить зависимость характеристик механического разворачивания белков от величины силы и скорости растяжения.

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

4. На примере белка L изучить влияние точечных аминокислотных замен на его механические свойства.

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

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

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

Апробация работы. Материалы диссертации были представлены на следующих конференциях: 11-я, 12-я, 13-я, 14-я Международная Пущинская школа-конференция молодых ученых «Биология - наука XXI века» (Пущино, Россия, 2007, 2008, 2009, 2010); И, III Международная конференция «Математическая биология и биоинформатика» (Пущино, Россия, 2008, 2010); IV Российский симпозиум «Белки и пептиды» (Казань, Россия, 2009); V Санкт-Петербургская конференция молодых учёных «Современные проблемы науки о полимерах» (Санкт-Петербург, Россия, 2009); VII

Национальная конференция «Рентгеновское, Синхротронное излучения, Нейтроны и Электроны для исследования наносистем и материалов. Нано-Био-Инфо-Когнитивные технологии РСНЭ-НБИК 2009» (Москва, Россия, 2009); XV Симпозиум по межмолекулярному взаимодействию и конформациям молекул (Петрозаводск, Россия, 2010); V Всероссийская Каргинская конференция «Полимеры-2010» (Москва, Россия, 2010); конференция молодых ученых «Современные проблемы теоретической физики» (Киев, Украина, 2010); XVIII Международная Конференция «Математика. Компьютер. Образование» (Пущино, Россия, 2011); XIX Российский национальный конгресс «Человек и лекарство» (Москва, Россия, 2012); IV Съезд биофизиков России (Нижний Новгород, Россия, 2012).

Публикации. Основные результаты диссертации представлены в 24 печатных работах, в том числе в 10 статьях в реферируемых научных изданиях.

Структура диссертации. Диссертация изложена на 103 страницах машинописного текста и состоит из введения, 6 глав (обзора литературы, изложения метода молекулярной динамики, обсуждения полученных результатов), благодарностей, выводов, и списка цитируемой литературы, включающего 109 ссылок. Работу иллюстрируют 26 рисунков и 12 таблиц.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ

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

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

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

Третья глава посвящена описанию объектов исследования, а также описанию постановки и проведения молекулярно-динамических экспериментов.

Четвертая глава посвящена исследованию механического разворачивания нативных белков Ь и С. Методом молекулярной динамики с использованием явной модели растворителя исследовались процессы механического разворачивания белков под действием приложенной внешней силы. С белками проводилось два типа вычислительных экспериментов: в одном случае к концам белка прикладывалась постоянная внешняя сила (разворачивание постоянной силой); в другом случае, концы белков равномерно удалялись относительно друг друга в противоположенные стороны (разворачивание постоянной скоростью).

В качестве объектов исследования были взяты два иммуноглобулинсвязывающих домена белков ЬиО (далее для краткости -белок Ь и белок в соответственно). Белок Ь (название файла в базе белковых структур - 2рИ) состоит из 63 аминокислотных остатков (остатки с 16-го по 78-й по последовательности) и содержит 951 атом. Белок О (название файла в базе белковых структур - ^Ь) содержит 56 аминокислотных остатков (853 атома). Оба этих белка состоят из двух р-шпилек, расположенных на концах белка (Ы- и С-шпильки), и а-спирали между ними. Белки имеют одинаковые пространственные структуры (среднеквадратичное отклонение выравненных Са-атомов этих белков составляет 1.38 А), но различные аминокислотные последовательности (рис. 1).

4 si Т-ЙС ЗЇТ 4- * XA ♦ + -i-GE+T DAK +Ї+

1-рФ и sue----s<

Рис. 1. а) Пространственное выравнивание структур белков L и G; б) упрощенное представление белков L и G, которое в дальнейшем используется при анализе полученных данных (цифрами обозначены концы элементов вторичной структуры); в) выравнивание аминокислотных последовательностей белков L и G.

Стартовые координаты атомов белков брались из банка белковых структур (Protein Data Bank). Белок помещайся в параллелепипед, заполненный молекулами воды. Вода, которая стерически перекрывалась с атомами белка, удалялась. Таким образом, были получены гидратированные белки L (1771 молекул воды) и G (1844 молекул воды). В дальнейшем, вся система белок + вода помещалась в сфероцилиндр с непроницаемыми отталкивающими стенками достаточно больших размеров, чтобы не влиять на динамику белка и прилегающую к нему воду, и в то же время не позволять испарившимся молекулам воды уходить на «бесконечность», возвращая их в область моделирования. В начальный момент времени всем атомам стартовой конфигурации системы приписывались случайные скорости, и проводилась молекулярная динамика при фиксированных положениях концевых атомов белка. Для каждого из белков проводилась серия из нескольких таких расчетов в независимых случайных столкновительных средах. Отрелаксированные таким образом в течение нескольких десятков пикосекунд системы являлись стартовыми для последующих вычислительных экспериментов.

Полный перечень проделанных вычислительньгх экспериментов приведен в таблице 1.

Нативные иммуноглобулинсвязывающие домены белков Ь и О

Число траекторий

Р=согкЛ, Т=350 К

Б, пН Белок Ь Белок С

600 10 10

700 10 10

800 10 10

900 10 10

1000 10 10

Т=соп51:, Р=1050 лН

Т,К Белок Ь Белок G

320 5 5

350 10 10

400 5 5

\'=со1к1 Т=350 К

V, А/пс Белок I, Белок в

0.0050 4 4

0.0625 10 10

0.1250 10 10

Белок с аминокислотными заменами

(49 одиночных замен на аланин и замена всех аминокислотных остатков на

аланин кроме глицинов)

у=0.05 А/пс, Т=350 К

150

(по 3 независимых вычислительных эксперимента при каждой

аминокислотной замене)

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

Графики зависимости силы от расстояния между концами белков приведены на рис. 2. Исследование этих зависимостей показало, что все полученные траектории можно поделить на три группы: траектории, на которых появляется один (в 6 случаях для белка Ь и в 8 для белка О), два (в 12 случаях для белка Ь и в 13 случаях для белка О) и три (в 6 случаях для белка ЬивЗ случаях для белка в) пика силы. Причем два пика силы как для белка Ь, так и для белка в наблюдаются в половине случаев не зависимо от величины скорости растяжения. Наличие нескольких пиков силы на

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

2100.14007000-

а БелокL

ШШ

JP* , - «

"Н~ГТГ"ГП~П 30 40 50 ео 70 80

Г КС, А

2100. 1400 -7000

б Белок G

1 I 1 I 1 I 1 I 1 I 20 30 40 50 60 70

Гыс, А

Рис. 2. Зависимость силы реакции F, действующей на концы белка, от расстояния между концами белка înc при различных скоростях растяжения (v=0.1250 А/пс - черные кривые, 0.0625 А/пс - светло-серые кривые, 0.0050 А/пс - серые кривые): а) белок L; б) белок G. Для скоростей 0.1250 и 0.0625 А/пс представлено по 10 траекторий для каждого белка, а для скорости 0.0050 А/пс - по 4 траектории.

Экспериментальные работы по растяжению этих белков в атомно-силовой микроскопии и теоретические работы, в которых используются грубые модели белков либо же растворитель учитывается неявно, указывают на то, что эти белки сворачиваются/разворачиваются по простому одностадийному механизму [Brockwell et al., 2005, Biophys. J., 89, 506; Cao and Li, 2007, Nat. Mater., 6, 109; Cao and Li, 2008, ./. Mol. Biol., 375, 316]. В опубликованных недавно экспериментальных работах [Waldauer et al., 2008, HFSP J„ 2, 388; Morrone et al., 2011, Biophys. J., 101, 2053] no разворачиванию белков L и G под воздействием денатуранта было показано, что этот процесс не описывается простой двухстадийной моделью разворачивания.

Наличие нескольких пиков силы на траектории не зависит от величины скорости растяжения. При каждой из использованных в данной работе скоростей растяжения есть траектории, на которых появляются как один, так и два, и три пика силы. От величины скорости растяжения зависит только высота пика силы, и естественно, что чем больше величина скорости растяжения, тем выше пик силы. Это объясняется тем, что при более высокой

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

На рис. За приведен график зависимости средней максимальной силы Fmax, которая возникает при таком механическом разворачивании белков, от величины скорости растяжения. Из этого графика видно, что при двух самых высоких скоростях оба белка характеризуются близкими значениями величины силы Fmax. Различия становятся заметны только для наименьшей из использованных скоростей растяжения. Для того чтобы развернуть белок О потребовалась большая сила, чем та, которая возникает при растяжении белка L. Из этого следует, что белок G - более механически стабилен, чем белок L. Это утверждение согласуется с экспериментальными данными из которых следует, что для разворачивания белка L надо приложить силу не менее 136 пН, а для белка G - 180 riH [Brockwell et al., 2005, Biophys. J., 89, 506; Cao and Li, 2007, Nat. Mater., 6, 109].

v, Á/пс F, nH

Рис. 3. Зависимости а) средней максимальной силы Ртох, возникающей при растяжении белков с постоянной скоростью; б) среднего времени начала разворачивания белков от величины силы растяжения.

В экспериментах при растяжении белков за концы с постоянной силой, следили за тем, как меняется расстояние между концами белка с течением времени. На рис. 4 приведены полученные зависимости для белков Ь и О и наборов независимых траекторий их разворачивания при различных величинах сил растяжения. Из этого рисунка видно, что при общей

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

О 10000 20000 о 400 800 1200 1600 2000 0 200 400 »00 800

t, ПС t, ПС t, ПС

Рис. 4. Изменение расстояния Гыс между концами белка от времени при заданной силе растяжения: а) 600 пН; б) 800 пН; в) 1000 пН. Черные кривые относятся к белку L, серые -к белку G. Разные кривые одного и того же типа относятся к различным (независимым) траекториям белка. Для каждой силы растяжения представлено по 10 траекторий для каждого белка.

f = воо пН

F= 1000 пН

На рис. 36 приведены зависимости среднего времени начала разворачивания белков Ь и в от величины силы растяжения. Этим временем считалось время, при котором расстояние между К- и С-концами белков начинает резко возрастать. Из рис. 36 видно, что при силе растяжения Р=600 пН для разворачивания белка в требуется порядка 8 не, а для белка Ь — 3 не. А из этого следует то, что белок в является более механически стабильным по сравнению с белком Ь в диапазоне малых сил растяжения (600 и 700 пН).

В таблице 2 приведены значения средних времен начала разворачивания белков Ь и в при силе растяжения Р=1050 пН при температурах Т=320, 350 и 400 К. Как и следовало ожидать, с повышением температуры время разворачивания для обоих белков уменьшается. Из таблицы 2 видно, что белок Ь разворачивается примерно в 1.5 раза быстрее по сравнению с белком в. Из этого опять же следует то, что белок в с повышением температуры остается более механически стабильным по сравнению с белком Ь.

Таблица 2. Среднее время начала разворачивания белков Ь и С при силе растяжения Р=1050 пН.

Температура, К 320 350 400

Белок Ь, время (пс) 224±Э6 (5 траекторий) 86±12 (10 траекторий) 71±22 (5 траекторий)

Белок в, время (пс) 558±261 (5 траекторий) 115±24 (10 траекторий) 10б±27 (5 траекторий)

Анализ изменения числа контактов между различными элементами вторичной структуры при механическом растяжении белков позволил выявить три возможных способа разворачивания белков Ь и Б. Во всех случаях белки сначала распадаются на две структурные единицы. В первом случае этими структурными единицами являются [Ы-шпилька + а-спираль] и [С-шпилька]. В этом случае первой разворачивается С-концевая р-шпилька, за ней а-спираль, и последней — К[-концевая р-шпилька. Во втором случае этими структурными единицами являются [Ы-шпилька] и [С-шпилька], а а-спираль не примыкает ни к одному из этих блоков. В данном случае первой разворачивается а-спираль, за ней С-концевая р-шпилька, и последней — Ы-концевая р-шпилька. И в третьем случае этими структурными единицами являются [Ы-шпилька] и [а-спираль + С-шпилька]. В этом случае первой разворачивается И-концевая р-шпилька, а за ней а-спираль и С-концевая Р-шпилька. Причем первый способ наблюдается наиболее часто (в 50 и 48 случаях из 74 для белков Ьиб соответственно), а третий способ наблюдался редко (как для белка Ь (в 7 случаях), так и для белка О (всего лишь в 2 случаях)) только при самых маленьких скоростях и силах растяжения.

Таким образом, было обнаружено, что при больших силах (Т > 800 пН) и скоростях растяжения (V > 0.0625 А/пс) средние времена и средние максимальные силы, которые требуются для разворачивания белков Ь и в близки. С уменьшением силы и скорости растяжения для разворачивания белка О в среднем требуется большее время и большая сила, чем для белка I,. А это говорит о том, что белок в — более механически стабилен по

сравнению с белком L в диапазоне малых скоростей и сил растяжения. Но в тоже время величина прикладываемой силы не изменяет пути механического разворачивания белков L и G.

Данное исследование показало, что и в белке L и в белке G самым механически стабильным элементом вторичной структуры является N-концевая р-шпилька. В экспериментах по разворачиванию этих двух белков под воздействием денатуранта было показано, что в ядро сворачивания белка L входит N-концевая Р-шпилька [Kim et al., 2000, J. Mol. Biol., 298, 971], a белка G - С-концевая р-шпилька [McCallister et al., 2000, Nat. Struct. Biol., 7, 669]. Таким образом, экспериментальные данные по разворачиванию этих белков под воздействием денатуранта, и теоретические данные по моделированию разворачивания этих белков под действием силы совпадают для белка L, а для белка G - нет. Экспериментально было показано, что пути механического и температурного разворачивания, либо разворачивания денатурантом не всегда совпадают.

Пятая глава работы посвящена изучению механического разворачивания белка L с точечными аминокислотными заменами. Для этого помимо нативной структуры белка L, были также получены 49 новых структур, отличающихся от нативной заменой одного из аминокислотных остатков на аланин. Затем моделировался процесс разворачивания полученных структур белков путем растяжения их за концы с постоянной скоростью (v=0.05 А/пс, Т=350 К). Для каждой аминокислотной замены было проведено по три независимых вычислительных эксперимента с различными начальными данными.

Анализ механического разворачивания белка L с точечными аминокислотными заменами показал, что есть замены, которые как понижают (в 6 случаях) максимальную силу реакции F, действующей на концы белка (пик силы), так и повышают ее (в 4 случаях). Все 6 замен, понижающих пик силы, располагаются во вторичной структуре белка. А из 4 замен, которые повышают пик силы, 3 располагаются во вторичной

структуре, и 1 в петле. Остальные 39 аминокислотных замен никак не повлияли на эту характеристику.

Так же была вычислена работа внешней силы по первичному разрушению структуры белка. Эта работа вычислялась как площадь под кривой зависимости силы реакции, действующей на концы белка, от расстояния между концами белка. Было показано, что эта работа в 19 случаях в среднем понижалась на 13%, по сравнению с нативной структурой. В остальных 30 случаях она не отличалась от работы, которую надо затратить на первичное разрушение нативного белка Ь. Из 19 замен, приводящих к понижению механической стабильности белка, 11 располагалось во вторичной структуре, а 8 в петлях. Полученный результат может служить указанием на то, что механическая стабильность белковой глобулы определяется не только мотивом и компактностью укладки вторичных структур белка, но и особенностями и составом аминокислот в соединяющих их нерегулярных участках (петлях) полипептидной цепи.

Средняя максимальная сила разворачивания для нативного белка Ь составляет 1312±87 пН (у=0.05 А/пс, Т=350 К). Для белка, в котором все аминокислотные остатки, кроме глицинов, заменены на аланины, становится равной 956±47 пН, т.е. в среднем понижается на 300 пН. Но при этом характер механического разворачивания такого белка не отличается от нативного. Этот белок разворачивается по наиболее вероятному пути, который характерен для нативного белка: первой разворачивается С-концевая Р-шпилька, за ней а-спираль, и последней - М-копцевая р-шпилька. При этом на кривой зависимости силы от расстояния между концами белка наблюдаются два пика силы, что говорит о появлении короткоживущих промежуточных состояний. Таким образом, механическая стабильность белков определяется взаимодействиями между аминокислотными остатками, а путь разворачивания зависит от топологии белка.

Известно, что белок G более термодинамически стабилен по сравнению с белком L (белок G разворачивается при 89°С [Stone M.J. et al., 2001, J. Am. Chem. Soc., 123, 185], а белок L - при 75°С [Gillespie В. et al., 2000, PNAS, 97, 12014]). И это может объясняться тем, что у белка G большее число межатомных контактов для внешних аминокислотных остатков, чем для белка L. Такая же картина наблюдается для белков из термофильных организмов, по сравнению с белками из мезофильных организмов. Анализ структурных различий между белками из термофильных и мезофильных организмов рассматривался в следующей главе.

Шестая глава работы посвящена изучению структурных факторов, ответственных за стабильность белков из термофильных организмов. Для этого была построена база, которая включала в себя 373 пары пространственно выравненных структур белков из термофильных и мезофильных организмов (http://bioinfo.protres.ru/resources/termo meso base.html). Термофилы - организмы, живущие при высоких температурах (50-80°С), мезофилы - организмы, живущие при нормальной температуре. В дальнейшем в тексте для краткости белки из термофильных организмов были названы просто термофилами, а из мезофильных - мезофилами. Исследуемая база белков была отобрана по следующим критериям: многодоменные белки были разделены на домены; длина доменов не превышала 400 аминокислотных остатков; разница длин белков в паре не превышала 10% от длины белка меньшего размера; домены, состоящие из двух или более отдаленных участков цепи, были исключены.

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

Типы остатков Доля, % Контактное расстояние 6 À Контактное расстошше 8 А

термофилы мезофилы термофилы мезофилы

Все 100 82.1±0.4 80.7±0.4 219.0±1.1 215.1±1.1

Внутренние 13 108.6±1.0 109.2±0.9 289.5±2.8 290.6±2.5

Внешние 46 61.0±0.2 59.5±0.2 163.6±0.6 159.4±0.6

Для данной базы доля аминокислотных остатков, вовлеченных в различные типы вторичной структуры, в термофильных и мезофильных белках одинакова и равна 80%. Из литературных данных известно, что термофильные белки более компактны [Robinson-Rechavi et al., 2006, J. Mol. Biol, 356, 547]. Это может быть следствием увеличения числа контактов между аминокислотными остатками у термофилов по сравнению с мезофиллами. Чтобы проверить это утверждение, было подсчитано число межатомных контактов на остаток для 373 пар белков термофил-мезофил при контактных расстояниях 6 и 8 А (считалось, что два атома находятся в контакте друг с другом, если расстояние между ними менее 6 А (8 А), межатомные контакты между соседними по цепи остатками в расчет не принимались). В термофильных белках наблюдается повышенное число межатомных контактов на остаток по сравнению с мезофильными (разница 1.4±0.4 для 6 А и 3.9±1.1 для 8 А, таблица 3). Чтобы определить за счёт каких аминокислотных остатков наблюдается такое повышение, в каждом белке аминокислотные остатки были разделены на промежуточные, внутренние и внешние. Остаток считался внутренним, если площадь его доступной для растворителя поверхности равнялась нулю; а внешним - если соответствующая площадь была больше 25% от максимальной площади поверхности остатка, доступной растворителю. Остальные аминокислотные остатки считались промежуточными [Fukuchi and Nishikawa, 2001, J. Mol. Biol., 309, 835]. Оказалось, что термофильные и мезофильные белки не отличаются по доли в них внутренних, внешних и промежуточных аминокислотных остатков: 13, 46 и 41%, соответственно. Было подсчитано число межатомных контактов на остаток для внутренних и внешних аминокислотных остатков термофильных и мезофильных белков (таблица 3). Внутренние аминокислотные остатки термофильных и мезофильных белков не отличаются по числу межатомных контактов на остаток. Это свидетельствует о том, что внутренние части термофильных и мезофильных белков одинаково плотно упакованы, а внешние аминокислотные остатки

термофильных белков имеют повышенное число межатомных контактов на остаток по сравнению с мезофильными (разница 1.5±0.2 для 6 А и 4.1±0.6 для

8 А).

Рис. 5. Доля каждого типа аминокислотных остатков, входящих в состав 373 пары термофильных и мезофильных белков: а) внутренние аминокислотные остатки; б) внешние аминокилотные остатки.

Также был проанализирован аминокислотный состав внутренних и внешних остатков термофильных и мезофильных белков. Из рис. 5а видно, что внутренние остатки термофильных и мезофильных белков в среднем не отличаются по своему аминокислотному составу. Среди внешних остатков термофильных белков наблюдается повышенное содержание по сравнению с мезофильными гомологами таких аминокислот, как лизин (разница 0.019±0.003), аргинин (0.017±0.003) и глутаминовая кислота (0.039±0.003). И, наоборот, среди внешних остатков мезофильных белков наблюдается повышенное содержание по сравнению с термофильными гомологами аланина (0.020±0.002), аспарагиновой кислоты (0.009±0.002), аспарагина (0.008±0.002), глутамина (0.022±0.002), треонина (0.011±0.002), серина (0.012±0.002) и гистидина (0.007±0.001) (рис. 56).

выводы

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

2. Показано, что на пути механического разворачивания белков Ь и в появляются короткоживущие промежуточные состояния.

3. Установлено наличие трех путей разворачивания белков Ь и & Во всех случаях нативные структуры белков сначала распадаются на две структурные единицы. В первом случае этими структурными единицами являются [Ы-шпилька + а-спираль] и [С-шпилька], во втором — [Ы-шпилька] и [С-шпилька] (а-спираль не примыкает ни к одному из этих блоков) и в третьем — [Ы-шпилька] и [а-спираль + С-шпилька]. В первом случае сначала разрушается С-концевая р-шпилька, во-втором - а-спираль, а в третьем - № концевая р-шпилька. Первый из сценариев разворачивания имеет большую вероятность (он встречается в 50 и 48 случаях из 74 для белков Ь и в соответственно), чем два других для обоих белков.

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

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

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи:

1. Глякина А.В., Лобанов М.Ю., Галзитская О.В. Поиск структурных факторов, ответственных за стабильность белков из термофильных организмов // Молекулярная биология. 2007. Т. 41. №4. С. 681-687.

2. Glyakina A.V., Garbuzynskiy S.O., Lobanov M.Yu., Galzitskaya O.V. Different Packing of External Residues Can Explain Differences in the Thermostability of Proteins from Thermophilic and Mesophilic Organisms // Bioinformatics. 2007. V. 23. Кг 17. P. 2231-2238.

3. Глякина A.B., Галзитская O.B., Балабаев H.K. Молекулярная динамика разворачивания белковых глобул под действием внешних сил // Вестник Тверского государственного университета. Серия Химия. 2008. Т. 68. №8. С. 109-119.

4. Глякина А.В., Балабаев Н.К., Галзитская О.В. Сравнение переходных состояний иммуноглобулинсвязываюгцих доменов белков L и G, полученных при моделировании разворачивания под действием внешних сил и в экспериментах под действием денэтуранта И Биохимия. 2009. Т. 74. С. 389-403.

5. Glyakina A.V., Balabaev N.K., Galzitskaya O.V. Mechanical unfolding of proteins L and G with constant force: similarities and differences H J. Chem. Phys. 2009. V. 131. №4. P. 045102.

6. Glyakina A.V., Balabaev NX, Galzitskaya O.V. Multiple unfolding intermediates obtained by molecular dynamics simulations under stretching for immunoglobulin-binding domain of protein G// The OpenBiochem. J. 2009. V. 3. P. 66-77.

7. Glyakina A.V., Balabaev NX, Galzitskaya O.V. Two-, Three-, and Four-State Events Occur in the Mechanical Unfolding of Small Protein L Using Molecular Dynamics Simulations // Protein and Peptide Utters. 2010. V. 17. №1. P. 92-103.

8. Глякина A.B., Галзитская O.B., Балабаев H.K. Исследование механических свойств иммуноглобулинсвязываюгцих доменов белков L и G методом молекулярной динамики // Компьютерные исследования tf моделирование. 2010. Т. 2. №1. С. 73-81.

9. Glyakina A.V., Galzitskaya O.V. Л comparative analysis of folding pathways of thermophilic and mesophilic proteins by Monte Carlo simulations II J Bioinform Comput Biol. 2010. V. 8. №3. P. 395-411.

10. Mamonova T.B., Glyakina A.V., Kumikova M.G., Galzitskaya O.V. Flexibility and mobility in mesophilic and thermophilic homologous proteins from molecular dynamics and FoldUnfoId method IIJ Bioinform Comput Biol. 2010. V. 8. № 3. P. 377-394.

Тезисы докладов:

11. Глякина A.B., Галзитская О.В., Балабаев Н.К. (2007) Молекулярная динамика разворачивания белковой глобулы под действием внешних сил // 11-я Международная Пущинская школа-конференция молодых ученых «Биология - наука XXI века». Россия, Пущино 2007. С. 56.

12. Глякина A.B., Балабаев Н.К., Галзитская О.В. Сравнение переходных состояний белков L и G, возникающих при моделировании их силового разворачивания и в экспериментах под действием денатурата // II Международная конференция «Математическая биология и биоинформатика». Россия, Пущино 2008. С. 147-148.

13. Глякина A.B., Галзитская О.В., Балабаев Н.К. Исследование путей механического разворачивания белков L и G // 12-я Международная Пущинская школа-конференция молодых ученых «Биология - наукаXXI века». Россия, Пущино 2008. С. 335.

14. Глякина A.B., Галзитская О.В., Балабаев Н.К. Влияние точечных аминокислотных замен на механические свойства иммуноглобулинсвязывающего домена белка L // IV Российский симпозиум «Белки и пептиды». Россия, Казань 2009. С. 314.

15. Глякина A.B., Галзитская О.В,, Балабаев Н.К Сравнение механических характеристик белков L и G И 13-я Международная Пущинская школа-конференция молодых ученых «Биология - наука XXI века». Россия, Пущино 2009. С. 258.

16. Глякина A.B., Галзитская О.В., Балабаев Н.К. Исследование механических свойств белков L и G // V Санкт-Петербургская конференция молодых учёных «Современные проблемы науки о полимерах». Россия, Санкт-Петербург 2009. С. 74.

17. Глякина A.B., Галзитская О.В., Балабаев Н.К. (2010) Влияние точечных аминокислотных замен на механические свойства иммуноглобулинсвязывающего домена белка L // VII Национальная конференция «Рентгеновское, Синхротронное излучения, Нейтроны и Электроны для исследования наносистем и материалов. Нано-Био-Инфо-Когнитивные технологии РСНЭ-НБИК 2009». Россия, Москва 2010. С. 439.

18. Глякина A.B., Галзитская О.В., Балабаев Н.К. Влияние точечных аминокислотных замен на механические свойства иммуноглобулин-связывающего домена белка L // 14-я Международная Пущинская школа-конференция молодых ученых «Биология - наука XXI века». Россия, Пущино 2010. С. 282.

19. Глякина A.B., Галзитская О.В., Балабаев Н.К. Исследование механических свойств белков методом молекулярной динамики // XV Симпозиум по межмолекулярному взаимодействию и конформациям молекул. Россия, Петрозаводск 2010. С. 94.

20. Глякина A.B., Галзитская О.В., Балабаев H.K. Сравнение механических свойств белков L и G. // V Всероссийская Каргипская конференция «Полимеры-2010». Россия, Москва 2010. 358_1.

21. Глякина A.B., Галзитская О.В., Балабаев Н.К. Молекулярно-динамическое моделирование механического разворачивания белков L и G // III Международная конференция «Математическая биология и биоинформатика». Россия, Пущино 2010. С. 177.

22. Глякина A.B., Галзитская О.В., Балабаев Н.К. Исследование механических свойств белков методом молекулярной динамики // Конференция молодых ученых «Современные проблемы теоретической физики». Украина, Киев 2010. С. 43.

23. Глякина A.B., Галзитская О.В., Балабаев Н.К. Исследование механических свойств белков методом молекулярной динамики // XVIII Международная Конференция «Математика Компьютер. Образование». Россия, Пущино 2011. С. 61.

24. Глякина A.B., Балабаев НХ, Галзитская О.В. Влияние структурных свойств белков на их механическую стабильность и термостабильность И IV Съезд биофизиков России. Россия, Нижний Новгород 2012. С. 74.

Подписано в печать:

26.04.2013

Заказ № 8442 Тираж - 100 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 www.autoreferat.ru

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

Федеральное государственное бюджетное учреждение науки ИНСТИТУТ МАТЕМАТИЧЕСКИХ ПРОБЛЕМ БИОЛОГИИ РОССИЙСКОЙ АКАДЕМИИ НАУК

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

04201357101

УДК 577.3

Глякина Анна Владимировна

ВЛИЯНИЕ СТРУКТУРНЫХ СВОЙСТВ ГЛОБУЛЯРНЫХ БЕЛКОВ НА ИХ МЕХАНИЧЕСКУЮ СТАБИЛЬНОСТЬ И ТЕРМОСТАБИЛЬНОСТЬ

Специальность: 03.01.02 - Биофизика

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

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

кандидат физико-математических наук Балабаев Николай Кириллович

доктор физико-математических наук Галзитская Оксана Валериановна

Москва-2013

Введение..................................................................................................................4

Глава 1. Обзор литературы.................................................................................8

1.1. Разворачивание белков в атомно-силовой микроскопии.........................8

1.2. Влияние силы на свободно-энергетический барьер разворачивания белка....................................................................................................................11

1.3. Анализ данных по механическому разворачиванию белков..................14

1.4. Ядро сворачивания белка. Ф-величины....................................................17

1.5. Экспериментальное и теоретическое исследование процессов механического разворачивания различных белков........................................19

1.5.1. Взаимосвязь между топологией белка и его механической стабильностью................................................................................................19

1.5.2. Влияние ориентации элементов вторичной структуры относительно направления приложения силы на механическую стабильность белка.... 26

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

1.5.4. Влияние денатуранта на механическую стабильность белка..........28

1.5.5. Промежуточные состояния на пути механического разворачивания белка.................................................................................................................29

1.5.6. Зависимость величины максимальной силы разворачивания от длины белка.....................................................................................................35

1.6. Термостабильность белков........................................................................35

Глава 2. Моделирование молекулярной динамики.....................................41

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

2.2. Т1РЗР модель воды......................................................................................44

2.3. Поддержание постоянной температуры и давления...............................46

2.4. Алгоритм вычисления взаимодействий валентно не связанных частиц (метод сканирования по пространству)...........................................................48

2.5. Численное интегрирование уравнений движения...................................50

Глава 3. Постановка вычислительных экспериментов..............................52

3.1. Объекты исследования...............................................................................52

3.1.1. Нативные иммуноглобулинсвязывающие домены белков Ь и О.... 52

3.1.2. Иммуноглобулинсвязывающий домен белка Ь с введёнными аминокислотными заменами.........................................................................53

3.2. Проведение молекулярно-динамических экспериментов......................54

3.3. Характеристики исследуемых моделей белков, рассчитываемые из полученных траекторий....................................................................................55

Глава 4. Моделирование механического разворачивания нативных иммуноглобулинсвязывающих доменов белков Ьив...............................57

4.1. Постоянная скорость (у=сопз1;)..................................................................57

4.2. Постоянная сила (Р=сопз1).........................................................................60

4.3. Пути механического разворачивания белков Ьив................................62

4.4. Расчет ядер сворачивания белков..............................................................68

Глава 5. Моделирование механического разворачивания белка Ь с

точечными аминокислотными заменами......................................................77

Глава 6. Поиск структурных факторов, ответственных за стабильность белков из термофильных организмов.............................................................82

6.1. Критерии отбора белков в базу.................................................................82

6.2. Структурные характеристики, рассчитываемые для белков из термофильных и мезофильных организмов....................................................84

ВЫВОДЫ..............................................................................................................88

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

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

Введение

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

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

Разворачивание отдельных белков растяжением их за концы стало возможным благодаря развитию атомно-силовой микроскопии (АСМ), которая открыла новые перспективы в изучении процесса сворачивания белков. Эксперименты по разворачиванию белков с помощью атомно-силового микроскопа обычно проводятся на многодоменных белках, либо для цепочки «сшитых» одинаковых или разных белков.

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

Экспериментально, при помощи АСМ при растяжении белков было обнаружено, что (3-структурные белки выдерживают существенно большие нагрузки, чем а-спиральные или а/(3 структурные белки. Считается, что упаковка элементов вторичной структуры белка является критическим фактором в определении его механической стабильности и термостабильности. Так, несколько водородных связей между соседними р-участками в р-листе стабилизируют структуру значительно сильнее, чем такое же количество гидрофобных контактов между спиралями в а-спиральных белках. Однако, не существует простого соотношения между структурой белка и его механической стабильностью и термостабильностью. Даже небольшие изменения в аминокислотной последовательности могут существенно повлиять на механические свойства белка и термостабильность. Направление приложения силы также является критическим фактором в определении механической стабильности. Детальное исследование путей силового разворачивания белка необходимо для понимания молекулярной природы механической стабильности белков.

Для достижения данной цели были поставлены следующие задачи:

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

4. На примере белка Ь изучить влияние точечных аминокислотных замен на его механические свойства.

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

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

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

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

1.1. Разворачивание белков в атомно-силовой микроскопии

Экспериментально разворачивание отдельных белков растяжением их за концы стало возможным благодаря развитию атомно-силовой микроскопии. В конце 1986 года Г. Биннингом с коллегами был разработан атомно-силовой микроскоп (АСМ), который в отличие от сканирующего туннельного микроскопа позволил изучать материалы, не проводящие электрический ток. Принцип работы АСМ основан на регистрации силового взаимодействия между поверхностью исследуемого образца и зондом. В качестве зонда используется наноразмерная игла, располагающаяся на конце упругой консоли, называемой кантилевером. Сила, действующая на зонд со стороны поверхности, приводит к изгибу кантилевера. Появление возвышенностей или впадин под остриём приводит к изменению силы, действующей на зонд, а значит, и изменению величины изгиба кантилевера. Таким образом, регистрируя величину изгиба, можно сделать вывод о рельефе поверхности. Под силами, действующими между зондом и образцом, в первую очередь подразумевают силы Ван-дер-Ваальса, которые сначала являются силами притяжения, а при дальнейшем сближении переходят в силы отталкивания [1,2].

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

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

образцом [3].

фото детектор

образец

лазер i ^ кантилевер

^подложка

о о

см

20 нм

Рис. 1.1. а) Схематическое изображение основных элементов атомно-силового микроскопа; б) характерная кривая зависимости силы от величины растяжения, получаемая в экспериментах с помощью атомно-силовой микроскопии. (Рисунок адаптирован из работы [4]).

Чувствительность АСМ позволила проводить эксперименты по изучению механических свойств белков. Эти эксперименты обычно проводятся на многодоменных белках, либо же в цепочку «сшивают» несколько одинаковых (или разных) белков, и затем эту цепочку растягивают.

До начала эксперимента цепочке из белков позволяют адсорбироваться на подложке. Затем иглой кантилевера цепляют за один ее конец и начинают растягивать. С увеличением расстояния между иглой кантилевера и

подложкой, растягиваемая цепочка белков генерирует возвращающую силу, которая вызывает изгиб кантилевера, а это в свою очередь вызывает отклонение лазерного луча, которое и регистрируется фотодетектором (рис. 1.1а). Такая система позволяет производить пространственные манипуляции менее чем на нанометр и может измерять силы в несколько пиконьютонов [4].

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

Для получения структурной интерпретации такой «пилообразной» зависимости силы от величины растяжения используются различные модельные подходы, в которых полноатомная (или грубая) модель [5, 6] белка растягивается за концы, и при этом используются различные способы учета растворителя (явно либо неявно) [7, 8].

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

наблюдается более сложное поведение с возможным существованием промежуточных состояний в процессе механического разворачивания [9-15].

Для разворачивания белков в АСМ используется два режима: режим постоянной скорости и режим постоянной силы. В первом случае подложка отодвигается от кантилевера с постоянной скоростью порядка 102-105 А/с. Во втором случае измеряемая сила сравнивается с заданным значением. Если эти силы различны, то генерируется сигнал «ошибки», который подается на усилитель, связанный с подложкой. Подложка изменяет свое положение посредством механизма обратной связи (время ответа порядка 20 миллисекунд), тем самым устраняется сигнал «ошибки» и поддерживается постоянное заданное значение силы [1].

1.2. Влияние силы на свободно-энергетический барьер разворачивания белка

Сворачивание белка описывается теорией переходного состояния. Из этой теории следует, что нативное и развернутое состояния белка разделены барьером свободной энергии. А скорости перехода через этот барьер имеют следующий вид:

ки = vk ехр

TS-N

kf = v/rexp

квТ J

КТ j

(1.1)

(1.2)

Где ки и к/ - константы скоростей разворачивания и сворачивания, соответственно; и - активационные энергии разворачивания и

сворачивания, соответственно; кв - константа Больцмана; Т - температура

в Кельвинах, V - частота колебаний; к - коэффициент перехода. Значение

предэкспоненциального множителя ук неизвестно для процесса сворачивания белка, но предполагается, что оно существенно ниже, чем для маленьких молекул [16].

те

Р(хи+хг)

Рис. 1.2. Свободно-энергетические барьеры разворачивания белка в отсутствии (серая кривая) и присутствии (черная кривая) дополнительной внешней силы. Под действием силы ^ барьер свободной энергии понижается на величину Рхи. Это стабилизирует переходное (ТБ) и развернутое (О) состояния, но не влияет на нативное (М) состояние. хи и х/ - расстояния до переходного состояния вдоль координаты

реакции для путей разворачивания и сворачивания соответственно. (Рисунок адаптирован из работы [1]).

Независимо от того, каким способом растягивается белок (в режиме постоянной силы, либо постоянной скорости), воздействие силы на белок интерпретируется использованием одной и той же теории (рис. 1.2). В отсутствии силы самое наименьшее по энергии состояние белка - это нативное состояние. Для того чтобы достигнуть развернутого состояния, белок должен преодолеть барьер свободной энергии. Высота этого барьера определяет скорость, с которой белок развернется вследствие тепловых флуктуаций. Под действием дополнительной внешней силы ^ барьер

свободной энергии понижается. Это увеличивает вероятность белка перейти в развернутое состояние. Экспоненциальный рост константы скорости разворачивания под действием внешней силы ки (р) описывается уравнением (1.3):

f Л ^ г-1 Л

K(F) = vk&IV - AGys-v FXu = к° exp

kBT

\квТ j

(1.3)

Где к° - константа скорости разворачивания при нулевой силе.

При растяжении белка в режиме постоянной скорости, сила (для механически стабильн�

Информация о работе