Авторефераты по всем темам  >>  Авторефераты по биологии  

РОССИЙСКАЯ АКАДЕМИЯ НАУК

ИНСТИТУТ ТЕОРЕТИЧЕСКОЙ И ЭКСПЕРИМЕНТАЛЬНОЙ БИОФИЗИКИ

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

БАШАРОВ МАХМУД АШЫГ-ОГЛЫ

ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ БАРЬЕРОВ ВНУТРЕННЕГО ВРАЩЕНИЯ И МЕХАНИЗМА СВОРАЧИВАНИЯ БЕЛКОВ

03.00.02 - биофизика

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

доктора физико-математических наук

Пущино - 2009

Работа выполнена в Институте теоретической и экспериментальной биофизики РАН

Научный консультант: доктор биологических наук Аллахвердиев С.И.

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

доктор химических наук, профессор Шишков А.В.

доктор физико-математических наук, профессор Сергеев Н.М.

доктор физико-математических наук, профессор Туманян В.Г.

Ведущая организация: Институт биофизики клетки РАН

Защита состоится У___Ф _____________ 20__ г. в ________ на заседании Диссертационного совета Д 501.001.96  при Московском государственном университете по адресу: 119991, Москва, ГСП-1 Воробьевы горы, МГУ, Биологический факультет

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

Автореферат разослан У___Ф _____________ 20__ г.

Ученый секретарь Диссертационного совета

доктор биологических наук, профессор _____________  Кренделева Т.Е.

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

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

Решением этой проблемы интенсивно занимаются во многих научных коллективах мира исходя из двух концепций. Согласно одной из них белки сворачиваются в процессе биосинтеза, по мере выхода полипептидной цепи из рибосомы, котрансляционно (рис.1). Концепция была обоснована тем, что синтез белковой цепи на рибосоме происходит с N-конца к C-концу направленно и подтверждается результатами многочисленных экспериментов in vivo и на белоксинтезирующих модельных системах. Однако неизвестно, как могут сворачиваться белки котрансляционно.

Рис.1. Ко- (слева) и посттрансляционное (справа) сворачивание белков

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

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

Итак, решение проблемы сворачивания белков требует ответить на обозначенные выше три фундаментальных вопроса: 1) сворачиваются ли белки ко- или посттрансляционно, т.е., в процессе, или после синтеза полипептидной цепи на рибосоме, 2) каков механизм сворачивания белков, 3) каково объяснение быстрого сворачивания белков?

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

Gстаб. = Hцепь - T Sцепь + Gгидрат. (*)

где Hцепь- энтальпия цепи, которая представляет собой внутримолекулярные и межмолекулярные взаимодействия полипептидной цепи и растворителя суммарно, Sцепь- энтропия цепи, Gгидрат.- свободная энергия гидратации (сольватации) цепи, которая включает в себе также свободную энергию растворителя, Т- температура. Предложены различные формулы и рецепты для вычисления каждого из входящих в выражение (*) вкладов. Эти рецепты, однако, далеки от совершенства. В частности, давно известно, что популярные функционалы, построенные на основе потенциалов невалентного взаимодействия типа У6-ехрФ или У6-12Ф и торсионных потенциалов, используемые для оценок потенциальной энергии полипептидов, страдают от очевидных серьезных недостатков. Поэтому является актуальным создание простых и надежных расчетных методов для реалистической оценки потенциальной энергии, внутримолекулярных и межмолекулярных взаимодействий полипептидов.

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

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

- Создание надежного расчетного метода для конформационного анализа и изучения взаимодействия полипептидов и других органических молекул - модифицированного метода фрагмент-фрагментных взаимодействий (ФФВ);

- Создание алгоритмов и вычислительных программ для: построения молекулярной модели олигопептидов, расчетов поверхности потенциальной энергии (ППЭ) олигопептидов и отдельных аминокислотных остатков в них, обработки данных банка белковых структур и вычисления конформации полипептидных цепей из атомных координат белков, статистического анализа конформаций аминокислотных остатков в трехмерных структурах белков, построения конформационных карт;

- Расчетно-теоретическое исследование ППЭ олигопептидов аминокислот методом ФФВ и оценка величины БВВ вокруг связей NC и CC полипептидной цепи белков;

- Анализ конформаций аминокислотных остатков в трехмерных структурах белков из банка белковых данных;

- Систематизация литературных данных по химически синтезированным белкам и выяснение их отношение к проблеме сворачивания белков;

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

- Выяснение степени адекватности ко- и посттрансляционного подходов изучению механизма сворачивания белков;

- Выдвижение возможного механизма сворачивания белков и предъявление подтверждающих этот механизм данных.

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

Разработан надежный расчетный метод для оценок барьеров внутреннего вращения (БВВ) в белках - модифицированный метод фрагмент-фрагментных взаимодействий (ФФВ) для конформационного анализа и изучения взаимодействия полипептидов и других органических молекул.

Установлено, методом ФФВ, что величины БВВ вокруг связей NC и CC основной цепи полипептидов, которые традиционно считаются незначительными, около 0,7 ккал/моль для обеих величин, достаточно значительны и составляют соответственно около 16 ккал/моль - что сопоставимо с транс-цис барьером пептидной связи, и 6 ккал/моль.

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

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

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

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

Научно-практическая значимость полученных результатов. Разработан расчетный метод для конформационного анализа и изучения взаимодействия органических молекул и полипептидов. Метод может быть использован для изучения структуру, конформаций, энергию взаимодействия, дипольные моменты различных молекул и молекулярных комплексов. Даны реалистические оценки величины БВВ вокруг связей NC и CC основной цепи полипептидов, что может способствовать развитию структурной протеомики и сыграть важную роль в исследованиях проблемы сворачивания и пространственной структуры белков. Обоснована концепция котрансляционного сворачивания белков, что существенно облегчает целенаправленный поиск эффективных путей для реалистического описания соотношения между аминокислотной последовательностью и пространственной структурой белков. Предложен механизм сворачивания белков, что позволит разработать адекватные алгоритмы и методы для предсказания пространственной структуры пептидов и белков на основе известной первичной структуры. Полученные по теме диссертации результаты используются в учебных курсах по молекулярной биофизике в ряде университетах США, по меньшей мере, и могут быть использованы в учебных курсах других вузов по молекулярной биофизике, биохимии, молекулярной биологии. Они могут найти активное применение в биотехнологических разработках и молекулярной медицине, при проектировании и создании новых эффективных лекарственных препаратов пептидной природы.

Апробация работы. Результаты работы были представлены на: VI и VII Всесоюзных симпозиумах по межмолекулярному взаимодействию и конформациям молекул (Вильнюс, 1982; Рига, 1986), I Всесоюзном биофизическом съезде (Москва, 1982), IХ Всесоюзном совещании по квантовой химии (Иваново, 1985), I международной школе-конференции молодых ученых стран-членов СЭВ по биофизике (Братислава, 1985), Российской научной конференции с участием зарубежных ученых Математические модели нелинейных возбуждений, переноса, динамики, управления в конденсированных системах и других средах (Тверь, 1994), Международной научной конференции Математические модели нелинейных возбуждений, переноса, динамики, управления в конденсированных системах и других средах (Тверь, 1996), International Symposium УProtein Structure, Stability and Folding. Fundamental and Medical AspectsФ (Moscow, 1998), 3-rd International Conference on Molecular Structural Biology (Vienna, 1999).

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

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

Объем, структура, общая характеристика диссертации. Диссертация состоит из введения, обзора литературы, результатов исследования и их обсуждения, выводов, списка опубликованных по теме диссертации и цитированных работ. Объем диссертации составляет 203 стр. и включает 39 рисунков, 37 таблиц и 498 ссылок. Во введении представлены актуальность проблемы и краткий анализ известных подходов к ее решению, цель исследования, отражена научная новизна работы. В обзоре литературы представлены успехи, достигнутые в исследованиях проблемы сворачивания белков. Раздел результаты и обсуждение состоит из трех частей, включающих 10 глав. В первой части рассматриваются основы гипотезы случайного клубка применительно к белкам, где внимание сосредоточено на обсуждении БВВ в белках. Для оценок БВВ в белках предлагается новая улучшенная версия расчетного метода ФФВ. Приведены результаты расчетов поверхности потенциальной энергии олигопептидов ряда аминокислот, полученные с применением метода ФФВ. Значительное место уделено результатам обработки данных о трехмерных структурах белков из банка белковых структур и анализу конформаций аминокислотных остатков встречаемых в этих структурах. Вторая часть посвящена анализу концепции посттрансляционного сворачивания белков. Результаты указывают на недостаточную обоснованность адекватности посттрансляционного подхода решению проблемы сворачивания белков. В третьей части вкратце прослежен биогенез белковой цепи от начала ее биосинтеза на рибосоме и до завершения сворачивания в разных водных компартментах клетки. Имеющиеся данные позволили обосновать, почему белки должны сворачиваться котрансляционно, и предложить возможный механизм сворачивания белков. Приведены данные, известные из литературы экспериментальные, и собственные, полученные в расчетах олигопептидов, которые поддерживают предлагаемую схему сворачивания белков.

СОДЕРЖАНИЕ РАБОТЫ

ОБЗОР ЛИТЕРАТУРЫ. Обзор занимает первые две главы работы. Первая посвящена результатам экспериментальных и теоретических исследований проблемы сворачивания белков на основе посттрансляционного подхода; приведено краткое описание популярных моделей сворачивания. Во второй главе представлены результаты исследований биосинтеза и сворачивания белков in vivo и на синтезирующих белок модельных системах, которые подтверждают котрансляционное сворачивание белков.

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

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

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

Анализу традиционных и выяснению реалистических оценок БВВ в белках посвящена первая часть работы, состоящая из четырех глав (гл. 3 - 6).

Часть первая. Барьеры внутреннего вращения (БВВ) в белках

Глава 3. Общепринятое представление о БВВ в белках. Исследователи, занимающихся проблемой сворачивания белков считают, что БВВ вокруг единичных связей NC и CC главной полипептидной цепи и большинства связей боковых цепей аминокислотных остатков белков незначительны, а вращение вокруг этих связей происходит почти свободно. Предположение о незначительности БВВ в белках было сделано в 1960-е годы, из-за невозможности определить их экспериментально. Рядом авторов было предложено, что они должны быть близкими к экспериментально известным малым значениям БВВ вокруг аналогичных связей в малых органических молекулах, в частности, различных производных амидной группы (Scheraga, AdPhysOrgChem,6,103,1968; Ramachandran&Sasisekharan,AdProtChem,23, 283,1968). Рекомендованные значения БВВ вокруг связей NC и CC главной цепи для обеих величин составляют в среднем около 0,7 ккал/моль. Именно эти, сравнимые с кТ при нормальной температуре, значения цитируют обычно при обсуждении величины БВВ в белках (Tanford, AdProtChem, 23,121,1968; Creighton, BiochJ,270,1990; Smith etal,Fold&Des,1,R95,1996; Fitzkee&Rose, PNAS, 101,12497,2004; Финкельштейн и Птицын Физика белка 2005). На предположении о незначительности БВВ построены и популярные модели (например, спиновое стекло, решеточные) и подходы к решению проблемы сворачивания (напр. использующие приближение случайной энергии) белков.

Глава 4. БВВ в белках значительно выше общепринятых для них значений: расчетно-теоретическое обоснование. Как известно, достаточно точные значения БВВ можно получить квантово-химическими ab initio методами. Расчеты поверхности потенциальной энергии (ППЭ) дипептидов методами ab initio проводятся регулярно с начала 1970-х. Что касается олигопептидов, рассчитать их методами ab initio невозможно, и поэтому нужны надежные и простые расчетные методы.

К концу 1970-х Головановым, Соболевым и Волькенштейном были разработаны новая потенциальная функция и алгоритмы с целью создания надежного расчетного метода для оценок внутри- и межмолекулярных взаимодействий молекул (ЖСХ:20,26,1981; ДАН СССР:259,877,1981). На основе этих разработок нами был создан метод связь-связевых взаимодействий для конформационного анализа и изучения взаимодействия органических молекул (Башаров и др. 1984аЦ1984е). Впоследствии нами были созданы и улучшенные версии метода: метод фрагмент-фрагментных взаимодействий (ФФВ) (Башаров и др.1989а, 1989б) и модифицированный метод ФФВ (Башаров 1995а). Последняя версия метода - она представлена в главе 6 - является одним из основных результатов диссертации. Метод воспроизводит достаточно адекватно ППЭ дипептидов, полученные квантово-химическими ab initio методами (Башаров 1995б, 1996), а также позволяет рассчитать ППЭ олигопептидов. При оценке величины БВВ в белках будем использовать результаты наших расчетов ППЭ ди- и олигопептидов ряда аминокислот методом ФФВ и расчетов дипептидов методами ab initio, сосредоточенных на БВВ вокруг связей NC и CC главной полипептидной цепи белков.

ППЭ дипептидов. Дипептиды глицина и аланина. Карты двумерной ППЭ глицинового и аланинового дипептидов, полученные методом ФФВ, представлены на рис.2, а относительные энергии некоторых конформаций и значения барьеров перехода между ними - в табл.1. Согласно этим данным, относительные энергии популярных конформаций аминокислотных остатков в белках и барьеры перехода между ними значительны. В частности, БВВ вокруг связи NC (C7(R)-αL; =0,=90) составляет не менее 16 ккал/моль, а вокруг связи CC (C7(R)-αR; = -90,=0) - не менее 5 ккал/моль.

Дипептиды других аминокислот. Нами были рассчитаны полные ППЭ Е (ϕ, ψ, χi) дипептидов серина, валина, фенилаланина и тирозина, где χi - углы вращения вокруг связей боковых цепей аминокислот (Башаров 1995б, 1996; неопубликованные результаты данной работы). Углы ϕ, ψ и χi варьировали в интервале от Ц180o до 180o, ϕ, ψ с шагом 10o, а i - 30o. Число рассмотренных конформаций каждого дипептида составило 37*37*13*n, где n число вращательных степеней боковой цепи аминокислотного остатка. Известны и расчеты ППЭ дипептидов ряда аминокислот методами ab initio.

Из результатов наших и ab initio расчетов следует, что ППЭ рассмотренных дипептидов в определенных конформациях их боковых радикалов мало отличаются от ППЭ дипептида аланина. Заметная разница лишь в том, что доступные области конформационного пространства значительно сужается с увеличением размера боковых цепей. Важным из результатов является то, что БВВ вокруг связей NCα и CαC рассмотренных дипептидов сравнимы с соответствующими величинами БВВ в дипептидах глицина и аланина (составляют в среднем около 17 и 7 ккал/моль; см. далее).

Рис.2. ППЭ дипептидов глицина (слева) и аланина (справа) по методу ФФВ. Вертикальное направление представляет вращение вокруг связи CC по углу , а горизонтальное - NC по углу . Крестиками обозначены минимумы. Цифры на эквипотенциальных линиях - относительная энергия в ккал/моль. Энергия областей заштрихованных параллельными линиями - от 15 до 25 ккал/моль, дважды заштрихованных - более 25 ккал/моль (Башаров 1995б).

Табл.1. Относительные энергии некоторых конформаций дипептидов глицина и аланина (энергия в ккал/моль, углы в градусах)

Конформация |  Относительная энергия

| Глицин  Аланин

ϕ ψ |  [1] [2]  [3] [4] [*]  [1] [2]  [5]  [*]

C5 180  180  0.0 0.0 0.0 0,8 0,7 1,5 2,5 1,4  3,4

C7(R)  -90  60  1.4 1.0 3.3 0,0 0,0 0,0 0,0  0,0 0,0

αR -60 -50  9.6 4.5 8.9 4,7 6,5 9,3 3,0  6,7 8,1

C7(R)-αR -90  0  7.5 5.5 5.5 7,9 9,5  11,0 4,9 7,9

C7(R)-αL  0 90 19.2  16.0  21.7  17,8 19,7  16,0  17,2

0  -90  25,5 18 24,5

C7(L) 60 -110 11,9 14  11,0

α(L)  60  50 10,6  5,0  8,4

Ab initio: [1]-Hiller&Robson,JTheorBiol,76,83,1979; [2]-Peters&Peters,JMolStr,85,107, 1981; [3]-Wright&Borkman,JPhysChem,86,3856,1982; [4]-Schaferetal,JChemPhys,72, 1439,1982; [5]-Scarsdaleetal,JACS,105,3438,1983; [*]- Метод ФФВ: Башаров (1995б).

ППЭ глицина и аланина в олигопептидах. Нами были рассчитаны ППЭ олигопептидов глицина и аланина, от ди- до декапептида в их спиральных конформациях, а также серединного 5-го и концевого 9-го остатков в декапептидах глицина и аланина в конформациях правой α-спирали (ϕ=-600,ψ=-400), близкой к вытянутой (ϕ=-1400,ψ=1500) и С7(ϕ=-800,ψ=800).

Величины БВВ вокруг связей NCα и CαC исследованных нами аминокислот в олигопептидах более наглядно отражены на кривых на рис.3 и 4 соответственно. Кривые представляют, приблизительно, одномерные профили предпочтительного пути конформационных изменений (изменения потенциальной энергии) остатков в пространстве Е (ϕ,ψ): на рис.3 в зависимости от угла ϕ (U (ϕ)) и на рис.4 в зависимости от угла ψ (U (ψ)), по всему интервалу их значений (Ц1800 ϕ, ψ +1800). Эти пути проходят по двум оврагам ППЭ, по углу ϕ приблизительно через седловую точку ϕ = 00,ψ = 900, и по углу ψ - ϕ = -900, ψ = 00, как это видно из карт на рис.2. Кривые построены на основе ППЭ соответствующих остатков следующим образом. На рис.3 на каждой кривой величина энергии при угле ϕi соответствует минимуму потенциальной энергии на ППЭ данного остатка в зависимости от угла ψ в интервале 0 ψ 180 при фиксированном ϕi. Аналогично, на рис.4 на каждой кривой величина энергии при ψi соответствует минимуму энергии на ППЭ в зависимости от угла ϕ в интервале -180 ϕ 0 при фиксированном ψi. Таким образом, кривые на рис.3 представляют потенциальную энергию вращения вокруг связи NCα, а на рис.4 - вокруг связи CαC соответствующих аминокислотных остатков в полипептидной цепи. Нулевой энергии на кривых на обоих рисунках соответствует предпочтительная конформация.

На кривых U (ϕ) локальному минимуму в правой части соответствует конформация левой -спирали, а величина потенциальной энергии при максимуме с координатой ϕ ~ 0 представляет величину БВВ вокруг связи NCα (рис.3). Аналогичным образом, на кривых U (ψ) локальному минимуму в правой части соответствует конформация правой -спирали, а величина потенциальной энергии при максимуме с координатой ψ ~ 0 представляет величину БВВ вокруг связи CαC (рис.4).

Общий вид карты ППЭ и расположение минимумов и седловых точек на ней, а также предпочтительные пути конформационных изменений по углам ϕ и ψ , описанные выше, в основном характерны также для других аминокислотных остатков, кроме пролина, в определенных конформациях их боковых радикалов. Утверждение известно из результатов расчетов ППЭ дипептидов с использованием эмпирических потенциальных функций, а также ряда не рассмотренных нами аминокислот методами ab initio.

Величины БВВ вокруг связей NCα и CαC аминокислот в дипептидах и олигопептидах, полученные в наших расчетах, представлены в таблицах в правой части рис.3 и 4. Для дипептидов приведены также величины БВВ, оцененные из имеющихся в литературе карт ППЭ, рассчитанные квантово-химическими ab initio методами.

Согласно представленным данным, БВВ вокруг связи NCα в дипептидах составляет 16 - 22 ккал/моль, в среднем 18,5 ккал/моль, по данным квантово-химических расчетов, и 15,1 - 18,2 ккал/моль, в среднем 16,7 ккал/моль, по результатам наших расчетов (рис3).

	Дипептид  Барьер  Ссылка Gly 18,2  [*] Ala 17,6  [*] Val 16,0  [*] >15 ,<20  [1] Ser 15,1 [*] >15, <20  [2] Phe  16,4  [*] Tyr  16,7  [*] Thr  >15 ,<20 [3]   Ile >15 ,<20 [1]   (a)
	n-Ala  Высота барьера [*] 2 17,6 3 16,0 4 15,4 5 15,1 6 14,8 7 14,7 8 14,5 9 14,4 10 14,4   (б)
	Пептид Высота барьера [*]   Gly9C7  16,7   Gly9al  17,4   Ala9C5  16,6   Ala9al  16,7   Gly5C7  16,2   Ala5C7  16,1   Gly5C5  16,4   Ala5C5  16,0   (в)

Рис.3. Предпочтительные пути конформационных изменений по углу ϕ (U (ϕ)) и величина БВВ вокруг связи NCα ряда аминокислот в (а)- дипептидах, (б)- олигопептидах аланина, (с)- 5-го и 9-го остатков в декапептидах Gly и Ala в конформациях С7, (al) и С5, полученные по методу ФФВ [*].Энергия в ккал/моль. Данные аb initio дипептидов: Gly и Ala см. табл.1; [1]ЦPeters&Peters JMolStr,88,157, 1982; [2]ЦPeters&Peters JMolStr,90,305,1982; [3]- Peters&Peters JMolStr,90,305,1982.

	Дипептид Барьер  Ссылка   Gly  7,6  [*]   Ala  7,4  [*]   Val  8,0  [*]   >6, <10  [1] Ser  5,3  [*]   6,3 [2] Phe 7,1 [*] Tyr 7,0 [*] Thr >6, <10  [3] Ile >6, <10  [1] (a)
	_________________________________________ n-Ala  Высота барьера [*] 2 7,4 3 6,3 4 5,8 5 5,5 6 5,3 7 5,2 8 5,1 9 5,0 10 4,9 ( б)
	Пептид  Высота барьера [*]   Gly9C7  7,9   Gly9al  6,7   Ala9C5  7,4   Ala9al  5,8   Gly5C7  7,0   Ala5C7  6,3   Gly5C5  7,3   Ala5C5  6,8   (в)

Рис.4. Предпочтительные пути конформационных изменений по углу ψ (U (ψ)) и величина БВВ вокруг связи CαC ряда аминокислот в (а)- дипептидах, (б)- олигопептидах аланина, (с)- 5-го и 9-го остатков в декапептидах Gly и Ala в конформациях С7, (al) и С5, полученные по методу ФФВ [*].Энергия в ккал/моль. Данные аb initio дипептидов: Gly и Ala см. табл.1; [1]ЦPeters&Peters JMolStr,88,157, 1982; [2]ЦPeters&Peters JMolStr,90,305,1982; [3]- Peters&Peters JMolStr,90,305,1982.

В олигопептидах глицина и аланина величина БВВ вокруг связи NC постепенно уменьшается с увеличением числа остатков, с 17,2 ккал/моль в дипептиде до 15 ккал/моль в пентапептиде. Она остается почти неизменной с октапептида, составляя 14,4 ккал/моль. Уменьшение барьера обусловлено малым увеличением энергии стабилизации кулоновского вклада с увеличением числа остатков. Основным составляющим барьера является энергия обменного отталкивания занятых молекулярных орбиталей. Для серединного 5-го и концевого 9-го остатков в трех избранных конформациях декапептидов глицина и аланина величина БВВ вокруг связи NCα составляет 15 - 16,8 ккал/моль, в среднем 16 ккал/моль.

Величина БВВ вокруг связи CC в дипептидах составляет 5 - 10 ккал/моль, 7,5 ккал/моль в среднем, по данным квантово-химических расчетов, и 5,3 - 8 ккал/моль, в среднем 6,6 ккал/моль, по результатам наших расчетов (рис.4). В олигопептидах аланина величина барьера несколько уменьшается с увеличением числа остатков, с 7,4 ккал/моль в дипептиде до 5,3 ккал/моль в пентапептиде. Уменьшение барьера обусловлено уменьшением вклада электростатического отталкивания. Величина барьера остается почти неизменной с гексапептида, составляя около 5 ккал/моль. Для серединного 5-го и концевого 9-го остатков в конформациях правой α-спирали (ϕ=-500,ψ=-500), близкой к С5 (ϕ=-1400,ψ=1500) и С7(ϕ=-800,ψ=800) декапептидов глицина и аланина величина БВВ вокруг связи CαC составляет 5,3 - 7,8 ккал/моль, в среднем 6,5 ккал/моль. Основным составляющим барьера во всех случаях является электростатическое отталкивание.

Влияние изменения молекулярной геометрии на ППЭ пептидов. Расчеты ab initio дипептидов с оптимизацией молекулярной геометрии показывают, что как относительные энергии отдельных конформеров, так и высота барьеров между ними претерпевают заметные изменения за счет значительной деформации молекулярной геометрии (Hiller&Robson,JTheorBiol, 1979,76,83; Head-Gordon et al,JACS,1991,113,5989). Однако, согласно данным белковой кристаллографии и кристаллохимии (Engh&Huber, ActCryst,1991,A47, 392), а также полученным нами из анализа молекулярной геометрии олигопептидов, заметные деформации в плотноупакованной и напряженной атомной системе, какой является полипептидная цепь, приводят к недопустимому сближению валентно-несвязанных атомов.

Итоги расчетов ППЭ олигопептидов. БВВ вокруг связей NCα и CαC главной цепи полипептидов. На основании представленных выше результатов мы заключаем, что величина БВВ вокруг связей NCα и CαC в полипептидной цепи белков составляет, приблизительно, около 16 ккал/моль и 6 ккал/моль соответственно. Надежность предлагаемых значений может быть оправдана надежностью методов расчета, с использованием которых они получены. Надежность широкого базиса квантово-химических ab initio методов, безусловно, вполне приемлема. Надежность метода ФФВ проверена на многочисленных органических молекулах различного класса, а также дипептидах. Несомненно, рассчитанная величина БВВ зависит в определенной степени от геометрии молекулы - мы использовали стандартные геометрии пептидов. Предлагаемые значения БВВ вокруг связей NCα и CαC, естественно, не претендуют быть абсолютно точными, но они однозначно обозначают порядок этих величин и, по нашему мнению, близки к действительным их значениям. Они значительно выше традиционно считаемых для них общепринятых - около 0,7 ккал/моль для обеих величин.

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

Теоретические предпосылки. Обратим внимание на предпочтительные пути конформационных изменений ряда аминокислотных остатков по углу при вращении вокруг связи NC, представленные на рис.3. Назовем отрицательной конформацию остатка с отрицательным значением угла , и положительной с положительным значением угла . Из рис.3 видно, что для неглициновых остатков практически запрещен переход из отрицательных конформаций в положительные через линию =180о. В тоже время глициновый остаток может переходить через эту линию в положительную конформацию и обратно беспрепятственно. Видно также, что возможные переходы неглициновых остатков между доступными отрицательными и положительными конформациями могут осуществляться через линию =0о путем преодоления NC барьера, который по нашим оценкам составляет около 16 ккал/моль. Величина достаточно значительна и почти сопоставима с хорошо известным значением барьера между транс- и цис- конформациями пептидной связи (около 18 ккал/моль), что должна делать практически нереализуемым переход неглициновых остатков из отрицательных конформаций в положительные, если принимать во внимание тот факт, что на рибосоме полипептидная цепь наиболее вероятно синтезируется в правой -спиральной конформации (Lim&Spirin,1985, JMB,188,565).

Обратим внимание и на локальный минимум потенциальной энергии в области положительных конформаций (рис.3), которому соответствует конформация левой -спирали. Он выше глобального минимума на ~8 ккал/моль в дипептидах, на ~4 ккал/моль в тетрапептиде и на ~1 ккал/моль в декапептиде. Поэтому приобретение неглициновыми остатками этой конформации не только затруднено кинетически из-за высоких барьеров, но и невыгодно остаткам находиться в этой конформации термодинамически, как следует из оценки по формуле Больцмана-Аррениуса.

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

Использовали данные банка о трехмерных структурах белков (PDB; версия 1990г.). Из банка отобрали файлы, которые содержали координаты атомов расшифрованных при разрешении 2Е и лучше белковых структур (всего 185). Белки принадлежали различным классам или типам по функциональным и семейственным отношениям: транспорта электрона, транспорта кислорода, медьсодержащие, нуклеазы, токсины, флавопротеины, ферросеропротеины, иммуноглобулины, лизоцимы, кальцийсвязывающие, ингибиторы, металлогидролазы, протеиназы, лиазы, зимогены, комплексы ферментов, и др. Из файлов PDB выбрали координаты атомов и вычислили длины валентных связей, валентные и двугранные углы для дальнейшей обработки. Обработку данных PDB и необходимые вычисления проводили с использованием собственных алгоритмов и программ.

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

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

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

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

Учитывая эти результаты, из 185 структур мы отобрали: наиболее точную структуру среди различных структур того же белка, структуру одного представителя различных лигандных, функциональных и мутантных состояний того же белка, а также белка выделенного из различных источников, структуру одной субъединицы из идентичных субъединиц белка. Число отобранных структур для дальнейшего анализа составило 81.

Анализ распределения конформаций аминокислотных остатков в белковых структурах. Общее число аминокислотных остатков в 81 структуре (без учета концевых остатков белковых цепей) составило 14706. В соответствии со сделанным из анализа ППЭ пептидов выводом, конформации глициновых и неглициновых остатков анализировали отдельно. Соответствующие карты распределения представлены на рис.5.

Рис.5. Карты распределения конформаций глициновых (слева) и неглициновых (справа) остатков в 81 белковой структуре (Башаров 1997).

В 81 структуре из 1371 глициновых остатков в положительной конформации находятся 798 (58,3%). Число неглициновых остатков и их относительная встречаемость в положительных конформациях в каждой структуре приведены в табл.4. Суммарное число неглициновых остатков во всех структурах составило 13335, из которых в положительных конформациях находятся 507 (3,8%). В 30-и структурах доля неглициновых остатков в положительных конформациях составляет не более 2%, в 12-и структурах - не более  1%. В семи структурах неглициновые остатки в положительной конформации не встречаются (табл.4).

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

Белковые структуры с серьезными ошибками. Из представленных в табл.4 структур выберем те, в которых доля неглициновых остатков в положительных конформациях составляет 7,6% и более (т.е. в два и более раза больше их средней встречаемости). Таких оказалось девять: 1HDS (разрешение 1,98Е, относительная встречаемость неглициновых остатков в положительных конформациях 8,1%), 3FAB (2Е, 12,5%), 1ACX (2Е, 9,8%), 1NXB (1,38Е, 10,9%), 3WGA (1,8Е, 8,7%), 1RNS (2Е, 8,5%), 2SNS (1,5Е, 13,1%), 2SOD (2Е, 8,1%) и 4LDG (2Е, 10,9%). Их анализировали более детально.

Табл.4. Встречаемости неглициновых остатков в положительной конформации в 81 структуре

а б в  г  д  е  а б в г д е  а б в  г  д  е

1ECA 134 123  0 0.0 1 2LHB  147 141  1 0.7 0 4HHB  283 263  2 0.8 1

1HDS 564 520 42  8.1  30  1LH1  151 145  2 1.4 0 1HMQ  111 106  2 1.9 1

1MBO  151 141  2 1.4 0  1CCR  109 98  1 1.0 0 2CCY 125 115 2 1.7 0

3CYT 101  89  0 0.0 0 2C2C 110 102  3 2.9 0  351C 80  73  1 1.4 0

2B5C 83  77  2  2.6 0 2CDV 105  95  1 1.1 0 1AZA  127  115  3 2.6 0

1PAZ 118 110  4 3.6 1 1PCY  97  87 0 0.0 2 4FD1 104 101 7  6.9 0

1FDX  52  48  3 6.3  0 3RXN  50  45 2 4.4 1 1HIP  83  78  2 2.6 0

4FXN  136 122  5 4.1 0 1FX1 145 127  5 3.9 0 1CPV 106  98 4  4.1 0

1ICB 73 68  0 0.0 0 5TNC 159 146  2 1.4 0  1FB4  441  401  17 4.2 3

3FAB  424 392  49  12.5  2 1REI  105  97  3  3.1 2  2RHE 112  99  4 4.0 0

1INS 47 44  0  0.0 1  1PPT 34 33  1  3.0 0  2OVO  54  50  2  4.0 1

4PTI 56 50 1 2.0 0 2CI2 63  61  2 3.3 0  1ACX  105 92 9  9.8 0

1MLT 24 22 0 0.0 0  2EBX  60 55 2 3.6 0 1SN3 63  54 2  3.7 1

1NXB 60 55  6  10.9 4  2CNA  265 246 10 4.1 0  3WGA  169  127 11 8.7 0

1UBQ 74 69  3 4.3  0  1CRN  44 40 0 0.0 0 2PAB 112  105  3  2.8 0

1GCR  172 159  9 5.7 0  2WRP  102 97 1 1.0 0 6LYZ 127  115  6  5.2 0

1LZ1 128 117  7 6.0 0  2LZM 162 151  2  1.3 0 2SNS  139 130 17  13.1  1

1RSM  122 119  1 0.8 2  1RNS 120 117 10 8.5 2 1RNT  102  90 3  3.3  2

5CPA 305 282  5  1.8 3  3TLN  314 278 16 5.8 1 9PAP  210 182  3  1.6 1

2ALP 196 165  4 2.4 1  2APR 323 276 3 1.1 2 2ACT  216 188  7  3.7 1

2APP 321 281  6 2.1 2  5CHA 231 209  8 3.8 0 2CGA  243 220  7  3.2 0

3RP2 222 204 8 3.9 1  3EST  238 213  5 2.3 0 2PKA  228 206  7  3.4 2

2PRK 277  244 8 3.3 1  1TON 224 204  7 3.4 1 2SGA  179 147  4  2.7 1

1TGN  220 195  7 3.6 0  1CSE  333 295  5 1.7 2  2CAB  254 238  4 1.7 2

2CTS 435 402  5 1.2 0  1GD1 331 308  7 2.3  0  1GP1  180 166 3 1.8 2

2CYP  291 267  7 2.6 0  3GRS 459  418  12 2.9 2  3DFR  160 150 3 2.0 2

2CPP 403 378  4 1.1  3 4LDG 327 303  33  10.9  0  2SOD  596 496 40 8.1 0

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

Анализ расстояний между валентно несвязанными атомами. Вычисляли эти расстояния и сравнивали их с уменьшенными на 0,5Е значениями стандартных наименьших контактных расстояний между соответствующими атомами. Межатомные расстояния, которые оказывались меньше выбранных для сравнения значений, считали аномально короткими. В результате анализа мы обнаружили 45 аномально коротких расстояний в структуре 1HDS, 233 - FAB, 21 - 1ACX, 11 Ц1NXB, 4 - 1RNS, 92 - 2SOD, 62 - 4LDG. В структурах 2SNS и 3WGA аномально короткие межатомные расстояния обнаружено не было.

Анализ валентных связей и валентных углов. Вычисленные длины валентных связей и валентные углы главной цепи белков сравнивали с их стандартными значениями. Считали их аномальными, если валентные связи отклонялись от стандартных более чем на 0,35Е, а углы более чем на 20о. Было выявлено в структуре 1HDS 215 аномальных валентных связей, а 1NXB и 1RNS по две. В 1HDS 126 валентных связей главной цепи отклоняются от стандартных более чем на 0,4Е, а 42 - более чем на 0,5Е. В структуре 1HDS имеется 480 аномальных валентных углов, 2SOD - 38, 3FAB - 2, NXB - 18, RNS - 5, 4LDG - 6. В 1HDS 160 валентных углов отклоняются от соответствующих стандартных более чем на 30o, а 71 - более чем на 40o. В структуре 2SOD таких углов имеются соответственно 10 и 3.

Сравнение полярности угла эквивалентных остатков в идентичных субъединицах 1HDS и 2SOD. Анализ показал, что в четырех идентичных субъединицах 2SOD из 40 неглициновых остатков в положительных конформациях только восемь составляют эквивалентную пару, по два остатка в каждой субъединице. В 1HDS из 45 остатков в положительных конформациях в двух идентичных парах субъединиц только четыре оказались эквивалентными (по одному остатку в каждой идентичной паре).

Цис-пептидные связи в 1HDS. Среди представленных в табл.4 структура 1HDS (гемоглобин лани) единственная, в которой в цис- конформации много пептидных связей (30). В тоже время в структуре других гемоглобинов - человека (4HHB) и лошади (2MHB), несмотря на высокую гомологию в первичной структуре, имеется всего одна цис- пептидная связь.

Хиральность остатков. Результаты анализа показали, что остатки Phe-33 первой, Val-112 второй и Asp-116 четвертой субъединицы 1HDS, Arg-1, Lys-51 и Glu-56 1NXB, как ни странно, находятся в D-изомерной форме.

Представленные выше данные свидетельствуют о наличии серьезных ошибок в структурах 1HDS, 1NXB и 2SOD. Относительно остальных шести структур можно отметить следующее. По заметкам авторов расшифровавших структуру 3FAB, она определена неточно и может содержать ошибки (Saul et al,1978,JBiolChem,253,585). В структуре 4LDG мы обнаружили заметное количество аномально коротких межатомных расстояний и аномальных валентных углов (см. выше). Сравнение этой структуры лактатдегидрогеназы с опубликованной в литературе новой уточненной структурой показало, что . вместо 34 неглициновых остатков в положительных конформациях в старой структуре, в новой их всего 16. Структура нуклеазы 2SNS уточнена при достаточно высоком (1,54Е) разрешении и в ней неглициновых остатков в положительных конформациях 17. Таких остатков в новой структуре оказалось всего 7. При анализе литературы мы обнаружили, что структура 1RNS несовершенна также: по заметке авторов структура лишена точности. Что касается 1ACX, помимо указанных выше ошибок, авторы отмечают, что координаты атомов получены в промежуточном этапе уточнения.

Итак, результаты анализа девяти белковых структур, в которых относительно много неглициновых остатков в положительных конформациях, показывают, что восемь из них (1HDS, 3FAB, 1ACX, 1NXB, 1RNS, 2SNS, 2SOD, 4LDG) содержат ошибки. Эти структуры мы исключили из дальнейшего анализа. Что касается девятой структуры - 3WGA, из 11 неглициновых остатков в положительных конформациях восемь по два находятся в четырех гомологичных доменах молекулы. Таким образом, наличие в белковой структуре относительно много неглициновых остатков в положительной конформации оказалось простым показателем несовершенства этой структуры. Неглициновые остатки в положительных конформациях в остальных 73 структурах анализировали также.

Отрицательные конформации неглициновых остатков в белках. Карты распределения конформаций глициновых и неглициновых остатков в 73 структурах представлены на рис.5. Сравнивая эти карты с аналогичными картами на рис.4 для 81 структуры визуально, можно убедиться в том, что они заметно УчистыеФ: распределение конформаций остатков в отдельных областях более компактное и резкое. Доля неглициновых остатков в положительных конформациях в 73 структурах составила 2,7% - 301 остатка из 11230, по сравнению с 3,8% в 81 структуре. Встречаемости каждого из 20 канонических аминокислотных остатков в положительных конформациях в 81 и 73 структурах приведены в табл.5.

Рис. 5. Карты Рамачандрана глициновых (а) и неглициновых (б) остатков в 73 белковых структурах (Башаров 1997).

Далее мы анализировали каждого из 301 неглицинового остатка в положительной конформации отдельно, путем инспекции их молекулярной геометрии на корректность. Длины валентных связей и валентные углы основной цепи остатков сравнивали со стандартными общепринятыми значениями соответствующих величин. Остатки, для которых рассчитанные из структурных данных наблюдаемые значения отклонялись от стандартных более чем на 0,2Е для валентных связей, 7o для угла NCC и 5o для остальных валентных углов, и пептидных связей от планарности более чем на 15o, классифицировали как с некорректной молекулярной геометрией. Таких остатков оказалось 46. Конформация 13 других остатков, по меньшей мере, попадала в запрещенные области конформационной карты.

Табл.5. Встречаемости 20 канонических аминокислотных остатков в положительных конформациях в 81 и 73 белковых структурах

а б  в

1  2 3 1 2 3  1  2  3

Trp 218 1 0,46 196 0  0,0  196  0  0,0

Pro 654 4 0,61 548 2  0,36  546  0  0,0

Val 1109 12 1,08 895 5  0,56  890  0  0,0

Ile 672  6 0,89 593  4 0,67 591 2  0,34

Thr 962 23 2,39  795  7 0,88 788  0  0,0

Leu 1115  26 2,33 920  12 1,30 911 3  0,33

Phe 533  8 1,50 454  7 1,54 450 3  0,67

Met  228  5 2,19 200  4 2,00 199  3 1,51

Glu 713  20 2,81 630  13 2,06  622 5  0,80

Lys  882 38 4,31 716  15 2,09  710 9  1,27

Arg  485 18  3,71  422  9  2,13 418 5  1,20

Ala  1247  41  3,29  1053  23  2,18 1039  9 0,87

Cys  365 11 3,01 315  7  2,22  311  3 0,96

Ser 1179 60 5,09 988  26  2,63  976  14  1,43

Gln  548 25 4,56 456  14  3,07  446  4  0,90

Tyr  536 22 4,10 481  16 3,33 472  7  1,48

His  375 23 6,13 281  11 3,91 276  6  2,17

Asp 801  51  6,37 689  33 4,79  675  19 2,81

Asn 716  113 15,78  596  93 15,60  551  48 8,71

13335  507 3,8 11233 301  2,68  11067 140 1,26

Gly 1371  798 58,20 1141 688  60,30  1141 688 60,30

Примечание. а- в 81 структуре, б- в 73 структурах, в- в 73 структурах без учета остатков, чьи координаты определены некорректно. Столбцы 1-3: общее число остатков, число остатков в положительных конформациях, относительное число остатков в положительных конформациях. - неглициновые остатки суммарно.

В результате анализа литературы мы обнаружили, что координаты атомов 111 остатков определены плохо: или их не видны на карте электронной плотности, или они имеют слабую электронную плотность, или обладают высоким температурным фактором (табл.5, столбец в). Таким образом, из 301 остатка атомные координаты 170 остатков нуждаются в уточнении. Если исключить эти остатки из расчета, то относительное число неглициновых остатков в положительных конформациях составило бы 1,2%.

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

Глава 6. Модифицированный метод фрагмент-фрагментных взаимодействий (ФФВ) для конформационного анализа и изучения взаимодействия олигопептидов и других биоорганических молекул. Оценка потенциальной энергии молекул. В 1970-е гг. в работах Голованова, Соболева и Волькенштейна (ДАН СССР: 1977,236,1140; 1977,237,375; ЖСХ,1979, 20,26;Мол.биол,1979,13,633) были высказаны соображения, которые можно было использовать для создания простого расчетного метода для конформационного анализа и изучения взаимодействия биоорганических молекул. Предполагается, что молекулы составлены из фрагментов, в качестве которых рассматриваются не атомы, а локализованные единичные связи (ЛЕС) и неподеленные электронные пары (НП). Энергия взаимодействия между двумя молекулами А и В (межмолекулярное взаимодействие - ММВ), или двумя группами А и В одной молекулы (внутримолекулярное взаимодействие - ВМВ) определяется суммарной энергией взаимодействия их фрагментов а и b  - Eab:

= (1),

Eab оценивается по выражению: Eab = Eкул. + Eобм.отт. + Eпер.зар. + Eдисп. , (2)

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

Eкул. = qi qj (Ri j)-1 ,  Eобм.отт. = K1 (SЗaЗb)2 (Rab)-3 ,

  Eпер.зар. = - K2 [(SЗаCb)2 + (SСаСb)2], Eдисп. = - K3 (Rab)-6  (3),

где qi  , qj - заряды фрагментов a и b, Rij - расстояние между зарядовыми центрами. SЗaЗb - интеграл перекрывания занятых молекулярных орбиталей (МО) фрагментов a и b, SЗаCb - интеграл перекрывания занятой МО фрагмента a со свободной МО фрагмента b, SCаЗb - интеграл перекрывания свободной МО фрагмента a с занятой МО фрагмента b. Rab - расстояние между центрами дисперсионных сил фрагментов, K1, K2, и K3 - параметры.

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

Дисперсионный вклад оценивается на основе формулы Лондона:

Е (D) = 1.51IaIb (Ia  + Ib )-1 αaαb (Rab)-6 , по выражению Едисп. = f(Rab)*E(D), (4)

где Ia  и Ib - потенциалы ионизации, а αa и αb - поляризуемости фрагментов a и b, Rab - расстояние между центрами фрагментов, f(Rab) - функция, медленно убывающая с расстоянием. Объединение выражений (1) - (4) дает общий вид функционала для вычисления энергии взаимодействия двух молекул А и В или двух групп А и В одной молекулы:

EAB = =  (qι qj (Rι j)-1 -1.51*f (Rab)*IaIb (Ia+Ib )-1αaαb (Rab)-6  - K2 [(SЗаCb)2 +(SСаСb)2] + K1 (SЗaЗb)2 (Rab)-3  (5)

Для вычисления EAB  по выражению (5), необходимо: смоделировать фрагменты (единичные, двойные и сопряженные связи, НП), выбрать характеристики фрагментов, уточнить выражения для вкладов и вид функционала, определить параметры.

На основе изложенного выше алгоритма нами был создан метод связь-связевых взаимодействий (Башаров и др.1984аЦ1984е) и его улучшенная версия - метод ФФВ (Башаров и др. 1989а,1989б). Метод ФФВ позволил реалистически оценить ВМВ и ММВ в малых молекулах. Однако он имел по меньшей мере два серьезных недостатка: а) ВМВ и ММВ оценивались независимо с использованием двух различных наборов параметров К1 и К2; б) для расчетов ВМВ был введен дополнительный параметр смягчения кулоновского взаимодействия НП, кратных и сопряженных связей. Кроме того, метод давал завышенные значения равновесных межмолекулярных расстояний для некоторых комплексов. Эти недостатки оставляли под сомнением применимость метода для расчетов длинных  цепных молекул - в том числе и полипептидов, в силу того, что таких молекулах взаимодействие отдаленных вдоль цепи фрагментов можно рассматривать как межмолекулярное, а взаимодействие близких - как внутримолекулярное.

Модифицированный метод ФФВ. Вклады оценивали по формулам (3) - (4), энергию взаимодействия - по выражению (5). Модели фрагментов использовали те же, что и в методе ФФВ. В отличие от предыдущей версии метода, дисперсионный вклад оценили иначе. Известно, что формула Лондона дает завышенные значения энергии взаимодействия в средних расстояниях. Этот эффект учитывается медленно убывающей с расстоянием функцией f(R). В предыдущей версии метода ее мы определяли как: f(Rab) = 1, если Rab > 4.5Е, а в расстояниях Rab < 4.5Е дисперсионный вклад считали постоянным и равным его значению при Rab=4.5Е: Едисп. = ЕD(Rab) = ED(Rab=4.5Е). В новой версии функцию f(R) мы определяли из расчетов тестовых молекул в виде: f(Rab) = (2Rab/(R0 + Rab))2, где R0 = Rab если Rab > 4.5Е, R0  = 4.5Е если 2Е < Rab < 4.5Е, R0  = 4.5Е, Rab = 2.0Е если Rab < 2.0Е. В предыдущей версии метода потенциал ионизации всех фрагментов мы считали одинаковым и равным 13,6эв. В новой версии каждый фрагмент имел собственный потенциал ионизации, которого мы определяли из потенциалов ионизации простых тестовых молекул. Далее нами были выбраны или определены, из расчетов тестовых молекул, новые характеристики и параметры ЛЕС, НП, кратных и сопряженных связей: величины и распределение зарядов, расположение орбитальных центров, зарядов и дисперсионных сил. Из расчетов тестовых молекул был найден один набор параметров K1 и K2 для оценок как ВМВ так и ММВ (вместо двух наборов в предыдущей версии). Подробно метод описан в диссертации и опубликован (Башаров 1995а).

Метод ФФВ и внутримолекулярное взаимодействие (ВМВ). Рассчитали БВВ и дипольные моменты значительного числа тестовых молекул (углеводороды насыщенные и с двойными и сопряженными связями, эфиры, альдегиды, кетоны, органические кислоты, аминокислоты, амиды, пептиды), для которых были известны экспериментальные или полученные в строгих ab initio расчетах значения соответствующих величин. Во всех рассмотренных примерах точно воспроизводятся предпочтительные и метастабильные конформации, БВВ и дипольные моменты молекул. Некоторые результаты приведены в таблицах 6 и 7.

Метод ФФВ и межмолекулярное взаимодействие (ММВ). В расчетах мы выясняли предпочтительную конфигурацию, равновесное расстояние, энергию связывания и дипольные моменты димеров и комплексов органических молекул: неполярных, полярных и неполярных с полярными. С физико-химической точки зрения, исследованные классы комплексов интересны тем, что они являются простыми моделями гидрофобных и гидрофильных взаимодействий в белках. Результаты расчетов комплексов мы не приводим; они представлены в диссертации в трех таблицах и опубликованы в печати (Башаров 1995а). Анализ отдельных вкладов в общую энергию показал, что предпочтительная конфигурация комплексов неполярных молекул формируется в основном за счет дисперсионного вклада. Определяющим вкладом в стабилизацию комплексов неполярных молекул с полярными, а также димеров и комплексов полярных молекул является кулоновский.

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

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

Таблица 6. БВВ (ккал/моль) и дипольные моменты ( дебай) насыщенных молекул

Молекула	БВВ			μ
	Экспер.	ab initio	ФФВ	Экспер.	ФФВ
CH3 ЦCH3	2.93	2.6 - 3.32	2.97
CH3 ЦC2H5	3.33	3.45	3.60
CH3 ЦCH2CH2CH3 т-г	3.4± 0.4	3.58	4.13
C2H5 - C2H5  г-г	6.1 - 7.4	5.72	8.38
Δ(г-г)	0.6	1.13	0.74
CH3 - CH(CH3)2	3.9±0.7	-	4.53
CH3 - C(CH3)3	4.5±0.2		5.40
CH3 - NH2	1.96	2.13	2.19	1.33	1.62
CH3 - NHCH3	3.28	3.62	3.66	1.03	1.40
CH3 - N(CH3)2	4.41	-	4.73	0.63	1.19
CH3 - CH2NH2	3.74		3.57		1.51
CH3 CH2- NH2 т-г	-	2.23	2.72		1.51
Δ(т-г)		0.52	0.44		1.51
CH3 - OH	1.09	1.06	1.34	1.71	1.60
CH3 - OCH3	2.72	2.98	3.17	1.30	1.01
CH3 - CH2OH			2.84		1.54
CH3 CH2 - OH  т-г		1.34	1.68		1.54
Δ(т-г)		0.63	0.49		1.54
CH3 - CH2 OCH3	3.3		3.11		1.22

Примечание. т-г - транс-гош; Δ(т-г) - разность энергий транс-гош. г-г - гошЦгош.

Таблица 7. БВВ (ккал/моль) в молекулах с кратными и сопряженными связями

Молекула	БВВ			μ
	Экспер	ab initio	ФФВ	экспер.	ФФВ
CH3 ЦCHCH2	1.98	1.55	1.89
CH3 ЦC(CH3)CH2	2.21	1.70	1.99
Цис- CH3 ЦCНCHCH3	0.731	0.42	0.57
Транс- CH3 ЦCНCHCH3	1.95	1.54	1.80
CH3 ЦC(CH)3CHCH2	2.03		2.24
CH3 ЦC6H5	0.014		0.02
Орто- CH3 ЦC6H4CH3	1.49		1.46
CH3 ЦCH=O	1.16	1.09	0.93	2.69	2.12
CH3 ЦC(CH3)=O	0.78		0.97	2.90	2.22
CH3 ЦCH=С=O	1.18		1.50	1.79	2.09
CH3 ЦC(CH3)=С=O	2.07		1.93	1.94	2.12
CH3 ЦCH=NCH3	1.64		1.90		1.74
CH3 ЦCH=N-OH -транс	1.83		1.98		0.94
CH3 ЦOCHO	1.19	3.68	1.39	1.77	1.70
CH3 ЦCOOH	0.48	1.1	0.51	1.70	1.62
CH3 ЦCONH2		0.15 - 1.16	1.14	3.4 - 3.9	3.76
CH3 ЦCONHCH3		0.06 - 1.10	0.94	3.6 - 4.2	3.78
CH3CONHЦCH3		-.73 - -.03	0.92	3.6 - 3.9	3.78
CH3 ЦNHCHO	0.1 - 1.2	0.12 - 0.86	0.82	3.4 - 3.9	3.66

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

Н- водородные связи в пептидах. Методом ФФВ нами была исследована возможность образования специфических Н водородных связей в комплексах амидов, пептидах и белках. В образовании таких связей в качестве акцептора протона выступают связи (пептидные связи основной цепи полипептидов и сопряженные связи боковой цепи ароматических и полярных аминокислотных остатков), а донором - атомы водорода при атомах N, O и C. В исследованиях на модельных системах нами было установлено, что энергия Н- связей может составлять до 2,5 ккал/моль (Башаров и др. 1986б, 1987б), и поэтому эти связи могут сыграть существенную роль в формировании и стабилизации пространственной структуры белков. На основании нашего предсказания Н- связи были обнаружены в структуре миоглобина, а впоследствии в трехмерных структурах многих белков.

С использованием метода ФФВ нами были рассчитаны ППЭ дипептидов глицина, аланина, валина, фенилаланина, серина и тирозина, гомологичных олигопептидов (от ди- до декапептида) глицина и аланина. Полученные результаты были представлены в главе 4.

Часть вторая. Степень адекватности посттрансляционного подхода решению проблемы сворачивания белков

Глава 7. Степень обоснованности гипотезы случайного клубка применительно к полипептидной цепи. Обычно в исследованиях по выяснению механизма сворачивания белков объектом исследования рассматривают состояние случайного клубка полипептидных цепей белков. Насколько обосновано такое рассмотрение? Выяснению этого вопроса посвящена настоящая глава работы.

Гипотеза случайного клубка и БВВ. В предыдущей части работы из расчетов ППЭ олигопептидов мы выяснили, что БВВ вокруг связей NC и CC главной цепи полипептидов значительны и составляют около 16 ккал/моль и 6 ккал/моль соответственно. Наличие высоких энергетических барьеров дает основание считать полипептидную цепь белка достаточно жесткоцепной макромолекулой. Понятие же случайный клубок, как известно, используется для описания гибкоцепных линейных макромолекул, в которых возможны свободные вращения. Это понятие было применено к белкам на основании предположения о незначительности БВВ в полипептидных цепях.

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

Анализ показал, что полученные в многочисленных экспериментах в течение около последних двух десятилетий результаты нового поколения по денатурации белков противоположны результатам классических экспериментов. Согласно новым результатам, почти все исследованные белки остаются достаточно компактными и обладают значительной остаточной вторичной и третичной структурами при любых, даже сильно денатурирующих условиях. Т.е., любое денатурированное состояние белка далеко от состояния случайного клубка. Такие результаты в частности были получены на белках фосфоглицераткиназа, иммуноглобулин-связывающий домен белка L, барстар, субтилизиновый ингибитор, цитохром c и его мутанты, рибонуклеаза А и ее S- пептид, ДНК связывающий домен 434-репрессора, FK506 связывающий белок, барназа, лизоцим, шаперонин GroEL, папаин, стафилококковая нуклеаза, химотрипсиновый игибитор-2, субтилизин, лактамаза, апомиоглобин, парвальбумин, -амилаза, лактальбумин, трипсиновый ингибитор, рибонуклеаза Т1, и многие др.

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

Глава 8. Синтетические белки и проблема сворачивания белков. Возможность получения синтетических аналогов нативных белков путем химического синтеза традиционно считают одним из серьезных доводов в пользу посттрансляционного сворачивания белков. Она, как полагают, также отвергает котрансляционное сворачивание белков, поскольку химический синтез белковой цепи обычно проводится в противоположном к биосинтезу направлении, от C-конца к N-концу. Мы впервые провели анализ данных о химически синтезированных белках. Результаты позволили нам дать достаточно адекватное объяснение - почему в экспериментах удавалось получить синтетические аналоги белков, вкратце следующим образом.

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

Глава 9. Эксперименты денатурации-ренатурации и проблема сворачивания белков. Денатурация и физические свойства белков. Эксперименты по денатурации и ренатурации традиционно используют как удобный способ для выяснения механизма сворачивания белков. Обычно эти эксперименты проводятся следующим образом. Объектом для исследования используют нативный белок. Сначала внешним воздействием (химические реагенты, температура, рН и др.) нарушают нативную конформацию - белок денатурируют. Далее воздействие убирают и следят за восстанавливаемостью той же исходной нативной структуры белковой молекулы во времени. Такой прием - обычный способ изучения физических свойств исследуемого объекта (твердых тел, полимеров и составляющих их молекул), таких, как упругость, эластичность, гибкость, деформация и т.д. Способность макромолекулы к конформационным изменениям тоже является важным ее физическим свойством. С этой позиции можно оценивать и процесс денатурации-ренатурации белков. При денатурации белок претерпевает конформационные изменения под внешним воздействием, преодолевая энергетические барьеры по вращательным степеням свободы. Глубина денатурации (степень нарушения нативной структуры) белка, очевидно, может зависеть как от гибкости (жесткости) самой полипептидной цепи, так и от характера и величины внешнего воздействия. Зависимость конформационных изменений белка от величины внешнего воздействия может быть представлена следующим образом. Если структурно-функциональные характеристики белка не претерпевают заметных изменений под воздействием внешней силы, то имеет место устойчивость (стабильность) - одно из важных физических свойств нативной структуры. В терминах конформации молекул это может означать, что полипептидная цепь не преодолевает барьеры по вращательным степеням свободы, а совершает крутильные колебания в пределах соответствующих нативной конформации потенциальных ям. Далее, с увеличением внешнего воздействия белок начинает денатурировать, преодолевая энергетические барьеры по вращательным степеням свободы, и переходить в отличное от нативного конформационное состояние. При этом могут проявиться, в зависимости от величины внешнего воздействия, свойства эластичности и невозвратной денатурации - другие физические характеристики белка. Показателем того или иного свойства может служить время спонтанной релаксации белка из денатурированного в нативное состояние после удаления денатурирующего воздействия. Область эластичности нативной структуры может быть идентифицирована "биологически приемлемым" релаксационным временем, протяженность которого составляет, допустим, до нескольких минут в лучшем случае. Если нативная конформация белка не восстанавливается за это время, значит, белок переведен в область необратимой денатурации. Проявление всех трех типов свойств белками в экспериментах по денатурации и ренатурации, известно, явление обычное.

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

Времена ренатурации денатурированных и in vivo сворачивания белков. Считается, что ренатурация денатурированного белка и сворачивание того же белка in vivo происходят за одинаковое время, которое сравнимо со временем биосинтеза полипептидной цепи. Это мнение, однако, и  противоречивое и неверное. Прежде всего, тем, что, согласно посттрансляционной концепции, полипептидная цепь начинает сворачиваться после синтеза и полного выхода из рибосомы. Поэтому нет никакой взаимосвязи между временами биосинтеза полипептидной цепи на рибосоме и сворачиванием этой цепи вне рибосомы. Процессы биосинтеза и сворачивания различны по своей природе: первый - это сопряженный химический и требующий энергию механохимический, тогда как второй - спонтанный физический. Оба процесса происходят независимо друг от друга и в различных местах, разделенных расстоянием между пептидилтрансферазным центром и местом выхода цепи из рибосомы (70 - 100Е), по меньшей мере. Биосинтез белка занимает время, равное произведению среднего периода элонгации цепи на один аминокислотных остаток (около 0,1 с/остаток) на число остатков в среднем. О времени же сворачивания белков достоверно известно лишь то, что белки проявляют активность сразу после выхода из рибосомы и даже раньше. Что же касается быстрого спонтанного восстановления нативного состояния денатурированных белков в ренатурационных экспериментах, это явление нерядовое. Во многих случаях регенерировать нативной конформации белков после денатурации задача достаточно проблематичная.

Белки сворачиваются, как традиционно считают, за времена 10-4с - 10-2с, а вторичные и компактные структуры в белках, как известно, формируются за времена 10-4с - 10-2с. В то время как наблюдаемое время для наращивания полипептидной цепи на один аминокислотный остаток рибосомой составляет в среднем 0,1с (см. третью часть). Из сравнения этих масштабов следует, что времена формирования вторичных структур и компактных состояний в полипептидной цепи, а также предполагаемые времена для сворачивания белков в целом, меньше времени изменения химической структуры полипептидной цепи на один аминокислотный остаток при биосинтезе. Именно поэтому мы считаем, что достаточно проблематично выяснить механизм формирования нативной конформации белка в исследованиях, которые используют целые белки или полноценной длины полипептидные цепи белков, как это делают в экспериментах по денатурации-ренатурации. Как было отмечено выше, полипептидные цепи денатурированных и ново синтезированных белков не имеют вид случайного клубка. Поэтому не имеет твердую основу утверждение, что полипептидная цепь сворачиваться в нативную структуру из этого состояния, а процесс сворачивания белка может быть представлен как перехода клубок - глобула.

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

В последние годы появляется все большее число данных, которые свидетельствуют о том, что функционально активное состояние многих белков не является наиболее стабильным термодинамически (Baker&Agard, Biochem:1994,33,7505; Thomas etal,TIBS:1995,20,456; Baker,NatStrBiol :1998,5,1021; Berkenpas etal,EMBO J:1995,14,2969; Sohl etal,Nature:1998,392,817). Объяснением этих наблюдений может быть то, что из-за высоких энергетических барьеров молекула белка может находиться в кинетически стабильном и термодинамически метастабильном состоянии, соответствующей нативному белку, достаточно долго, не имея возможности приобрести соответствующую глобальному минимуму свободной энергии конформацию. Высокие энергетические барьеры в полипептидной цепи, как мы установили (глава 4), имеет место.

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

Часть третья. Котрансляционное постепенное сворачивание белков

Результаты множественных экспериментов in vivo и на модельных системах свидетельствуют о котрансляционном, во время биосинтеза сворачивании белков (введение, глава 1). Принимая во внимание эти результаты, а также представленный во второй части анализ степени адекватности посттрансляционного подхода решению проблемы сворачивания белков, мы пытались найти обоснование, почему белки должны сворачиваться котрансляционно и выяснить возможный механизм процесса сворачивания на этой основе. Чтобы решить эти задачи, мы полагали, что необходимо проследить биогенез белков от начала их биосинтеза на рибосоме и до завершения сворачивания. При этом учитывать следующие, по меньшей мере, факторы: 1) механизм, способ и скорость импорта полипептидных цепей из рибосомы в цитоплазму для сворачивания, 2) возможные изменения в химической структуре - аминокислотной последовательности - и пространственной структуре полипептидной цепи, которые она может претерпевать в цитоплазме с начала выхода из рибосомы и до завершения сворачивания, 3) силы, обуславливающие сворачивание белков, 4) времена сворачивания белков и белковых фрагментов.

Глава 10. Сворачивание белка около рибосомы. Почему белки должны сворачиваться котрансляционно? Как известно, белковая цепь строится в пептидилтрансферазном центре рибосомы с N-конца, по одному аминокислотному остатку. Синтезируемая цепь продвигается далее по узкому рибосомному каналу и котрансляционно выходит из рибосомы в цитоплазму (рис.6).

Известно, что рибосомный канал заполнен водой, а его стенки обладают редгезивным свойством, - полипептидная цепь не прилипает к стенкам канала. Установлены размеры канала, которые составляют, по разным данным, по длине 50 - 100Е, и в диаметре 15 - 28Е. Известно также, что полипептидная цепь продвигается по каналу и выходит из него без резких изгибов в развернутом виде.

Рис.6. Схематическое представление синтезирующей полипептид рибосомы. (А)- узнающий антикодон участок, (Р)- пептидилтрансферазный центр, (Ex)-место выхода цепи из рибосомы. (Адаптирован из: Спирин, Молекулярная  биология том 1 (1986).)

Обратим внимание на скорость выхода полипептидной цепи из рибосомы, т.е. скорость элонгации цепи на рибосоме. Средние скорости элонгации полипептидных цепей различных белков определены в разных клетках и при различных условиях в многочисленных экспериментах in vivo и на модельных системах. Определены также абсолютные скорости трансляции отдельных кодонов м-РНК в E.coli. Наибольшая скорость составила 21,6 кодонов/с, следовательно, наименьшее время для трансляции одного кодона - 0,046 секунды. Учитывая это значение, можно даже предполагать, что наименьший период трансляции кодонов - периодичность поступления остатков в цитоплазму, не быстрее 0,01 с/остаток (т.е., меньше наименьшего наблюдаемого периода почти в пять раз). Справедливость нашего предположения обоснована (Basharov,2003). Можно посмотреть и цитированные там литературу, где обсуждается взаимосвязь скорости и надежности трансляции м-РНК на рибосоме, а также принимать во внимание данные о том, что в водной среде максимальная скорость питаемого нуклеозид-трифосфатом клеточного мотора не быстрее 102/сек (Vale,1996, JCellBiol:135,291; Kinosita et al,1998,Cell:93,21).

Далее сосредоточимся на данных о временах сворачивания белков и белковых фрагментов. Хорошо известно, что характеристические времена конформационных изменений молекул составляет 10-12 - 10-7с, -спирали и -листы в полипептидах формируются за 10-7 - 10-4с (Eigen,1968,ChimPhys:65,53), а компактные состояния за 10-5 - 10-2с (Eaton et al,1997,CurrOpinStrBiol:7,10). Сворачивание малых однодоменных белков, как полагают, происходит 10-4 - 10-2с (Ptitsyn,1991,FEBS Lett,285,176; Sosnick etal,1997,PNAS:94,8545). Белки сворачиваются очень быстро также in vivo (Kolb et al, JBC:275,16597).

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

Примерный механизм сворачивания белков: котрансляционное постепенное формирование нативной конформации. Как было отмечено выше: а) аминокислотные остатки синтезируемой полипептидной цепи поступают из рибосомы в цитоплазму для сворачивания с N-концевого, последовательно друг за другом и постепенно с периодичностью не быстрее 0,01с; б) времена сворачивания белковых фрагментов и белков меньше этого времени. Опираясь лишь на эти данные, и учитывая наблюдения, о том, что рибосома и синтезируемая полипептидная цепь не образуют ассоциаты ни при синтезе цепи и ни после его завершения, можно сформулировать следующую простую и вполне очевидную примерную схему процесса сворачивания белка в непосредственной близости рибосомы.

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

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

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

- На всех этапах процесса, формирование структуры N-концевой части цепи, а также на последнем этапе - нативной структуры белка, происходит не дольше 0,01с. Общее время процесса сворачивания белка имеет такой же порядок, что и суммарное время синтеза цепи на рибосоме.

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

Анализ известных моделей о возможном механизме сворачивания белков. В данном параграфе диссертации представлен анализ известных моделей о возможном механизме котрансляционного сворачивания белков и указаны их недостатки. Результаты анализа опубликованы (Башаров 2000б).

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

Полипептидные цепи предшественников белков пересекают одну или несколько мембран, чтобы попасть в компартменты сворачивания. (Пути транспорта определяет расположенный впереди N- конца цепи прекурсор, состоящий из одного или нескольких расщепляемых сигнальных пептидов.) Цепи поступают в компартмент через узкий гидрофильный белковый канал в мембране диаметром 20 - 26Е. Канал является частью системы импорта белков, транслокона. Все транслоконы сильно похожи (рис.7). Они состоят из нескольких мембрано-встроенных, -причаленных или -ассоциированных гетероолигомерных белков, которые формируют четыре функционально разные системы: узнавания и причаливания сигнального пептида, проводящий белок узкий гидрофильный трансмембранный канал, узнавания и процессинга сигнальной последовательности, транслокационный мотор.

Транслокация полипептидных цепей в различные компартменты сворачивания происходит по одинаковому механизму. Предшественники полипептиды импортируются в компартменты в развернутом виде NЦконцом (рис.7). Сигнальный пептид отщепляется сигнальной пептидазой котранслокационно, как только NЦконец цепи поступает в компартмент. Любая другая реакция процессинга полипептидной цепи в компартменте происходит котранслокационно также (напр., гликозилирование остатков Asn в люмене ЭР). Импорт полипептидов в компартмент осуществляется посредством транслокационного мотора, за счет энергии АТФ или ГТФ.

Рис.7. Система и механизм импорта полипептидов в митохондрии (а) и хлоропласты (б) (Haucke&Schatz,1997,TICB,7,103), бактерии (E.coli) и эукариот (S.cerevisiae) (Pohlschroder et al,1997,Cell,91,563). Импорт белков из цитоплазмы осуществляется транслоконом внешней (OM) и внутренней (IM) мембран митохондрий (Tom, Tim) и хлоропласта (Oep/lap, Iep/lap). По мере синтеза пробелки могут ассоциировать в цитоплазме с шаперонами (Hsp70, MSF), которые доставляют их к рецепторным белкам на внешней мембране. SRP - сигнал-узнающая частица. SecYEG и Sec61 формируют проводящий белок канал в мембране. Hsp70,mHsp70,sHsp100,SecA,Kar2 - транслокационные моторы.

Итак, согласно известным данным молекулярной клеточной биологии, изложенным выше вкратце, транслокация полипептидных цепей белков в отдаленные от рибосомы компартменты клетки происходит определенным,  одинаковым как из рибосомы, образом. Цепи поступают в эти компартменты N-концом в развернутом виде из узкого белкового канала в мембране, как из канала рибосомы в цитоплазму (рис.8).

Выясним теперь скорость транслокации полипептидных цепей через мембрану в компартменты, которая почти не исследована. Представляется разумным, если предположим, что транслокация цепей через мембраны осуществляется при скоростях не быстрее высокой скорости элонгации полипептидной цепи на рибосоме, т.е. с периодичностью 0,01 с/остаток (см. предыдущий з). Предположение справедливо в случае котрансляционной транслокации полипептидной цепи в люмен ЭР, по меньшей мере, т.е. когда рибосома причалена к мембране ЭР и непосредственно синтезирует цепь прямо в люмен через транслокационный канал мембраны (рис.8б).

Рис.8.Полипептидные цепи поступают в различные компартменты из узкого канала рибосомы (а) или транслокона в мембране (b, с)  и сворачиваются там во время транслокации, ко- (а, b), или посттрансляционно (с) (Basharov,2003).

Предположение поддерживается оценочными экспериментальными данными: от 20 до 30 остатков переходит через мембрану E.coli в периплазму за время от нескольких сек. до одной мин. (Economou etal,1995,Cell: 83,1171), максимальный период транслокации цепи через внутреннюю мембрану митохондрий посредством транслокона ~0,7 с/остаток (Schwartz & Matouschek, 1999,PNAS:96,13086). Предположение оправдано, если учесть уже упомянутые данные о том, что в водной среде максимальная скорость клеточного мотора не быстрее 102/сек, и тем, что рибосома и транслоконы обладают общими структурно-функциональными особенностями в части выполнения транслокационных функций. Рибосома сама является транслокационной машиной и состоит из тех же функциональных элементов, что и транслоконы, за исключением системы процессинга сигнального пептида.

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

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

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

Первое положение экспериментальный факт, установленный клеточными биологами (Nissen etal,2000,Science:289,920; Kramer etal,2001,IntJBioch&CellBiol, 33,541). В этой связи мы обращаем особое внимание на экранирующую роль гидратного слоя вблизи полярной белковой поверхности, около которой сворачивается белок. Второе положение хорошо известные данные о временах формирования вторичных и компактных структур в полипептидной цепи. Оно подтверждается также результатами экспериментов in vivo, согласно которым белки проявляют активность сразу после выхода из рибосомы, или даже пока продолжается синтез полипептидной цепи.

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

Динамичность конформации в зависимости от длины полипептидной цепи. Экспериментальная поддержка модели. Конформационный анализ десяти различной длины фрагментов полипептидной цепи химотрипсинового ингибитора (ХИ2; 64 остатков) проведен методами КД, ЯМР и др. (Prat Gay et al,1995a,JMolBiol.254,968; 1995b,PNAS,92,3683). Нативная структура ХИ2 составлена из одной -спирали (остатки 12-24), шести листов (3-5, 8-11, 28-34, 45-51, 55-58, 60-64) и четырех поворотов (5-8, 25-28, 35-42, 52-54). Фрагменты состояли из 5 (ХИ2-5), 13 (ХИ2-13), 25 (ХИ2-25), 28 (ХИ2-28), 40 (ХИ2-40), 50 (ХИ2-50), 53 (ХИ2-53), 60 (ХИ2-60), 62 (ХИ2-62) и 63 (ХИ2-63) N-концевых остатков. Пептиды ХИ2-5 и ХИ2-13 проявляют признаки поворота, тогда как в нативном белке остатки 1-13 образуют два листа и один поворот. Пептиды ХИ2-25 и ХИ2-28 неупорядоченно компактны, тогда как в нативном ХИ2 остатки 12 - 24 образуют стабильную -спираль. Пептид ХИ2-40 обладает компактной структурой с гидрофобным кластером, а ХИ2-49 компактнее, но нативнооподобные взаимодействия в них отсутствуют. Фрагмент ХИ2-53 обладает хорошо сформированным мини ядром, а ХИ2-60 свернут сильнее и более компактен. В формировании нативноподобной структуры ХИ2 определяющую роль играют пять С-концевых остатков. Эти результаты свидетельствуют о том, что пространственная структура цепи ХИ2 при ее постепенном удлинении динамична - вторичная и третичная структуры полипептида постепенно меняются по мере его роста.

В исследованиях зависимости структуры семи фрагментов барназы от длины биофизическими методами установлено, что вторичная и третичная структуры фрагментов формируются вместе и динамичны: компактность неуклонно растет в фрагментах с длиной от 22 до 68 остатков по мере удлинения цепи, потом резко падает в фрагменте с длиной 79 остатков, а затем снова растет (Neira&Fersht,1999, JMolBiol:287,421).

Результаты исследований зависимости структуры N-концевых фрагментов апомиоглобина с длиной 36, 77 и 119 остатков представлены на рис.9, где также приведены известные схемы сворачивания этого белка.

Рис.9. Схемы сворачивания апо-Мб. а-ренатурация из денатурирован-ного состояния (Jennings&Wright, 1993,Sci:262,892;ЯМР); b- гипотети-ческий механизм на основе посттр-ансляционного подхода: цепь имеет вид случайного клубка пока не синтезирована полностью; c-другой гипотетический механизм: вторич-ная структура формируется по мере синтеза соответствующего сегмента цепи; d- механизм, предложенный на основе ЯМР исследований (Chow et al., 2003).

Основной вывод этих исследований такой же, как и в случае с белком ХИ-2: вторичная и третичная структуры динамичны и значительно меняются с удлинением цепи. Содержание элементов вторичной структуры исследованных фрагментов и нативного апомиоглобина приведены в табл.8. Сравнивая данные соответствующих столбцов панелей а и d рис.9 и табл.8 можно убедиться и в том, что пути формирования пространственной структуры апомиоглобина из денатурированного мочевиной состояния и при постепенном росте его полипептидной цепи отличаются кардинально.

Таблица 8. Содержание элементов вторичной структуры в апомиоглобине (апо-Мб) и его фрагментах

Фрагмент апо-Мб	Содержание элементов (в долях)
				+ t	rc
(1 - 36) апо-Мб	0.04		0.40		0.62	0.34
(1 - 77) апо-Мб	0.16		0.29		0.51	0.33
(1 - 119) апо-Мб	0.30		0.18		0.40	0.30
(1 - 153) апо-Мб (нативный)	0.50		0.04		0.28	0.21

Примечание. Символы , , +t и rc соответственно обозначают: спирали, листы, листы плюс повороты, неупорядоченные структуры (random coil).

Зависимость стабильности спиральных конформаций олигопептидов от длины цепи. Динамическое изменение конформации полипептидной цепи в зависимости от ее длины следует из наших расчетов ППЭ олигопептидов глицина и аланина методом ФФВ. В качестве примера в табл.9 приведены относительные энергии двух конформеров последовательных олигопептидов этих аминокислот. Конформеры соответствуют популярным конформациям: полностью вытянутой С5 (=±180o, =±180o) которая, как считают, наиболее удобно представляет конформацию статистического клубка, и -спиральной (=-60o, =-40o).

Таблица 9. Относительные энергии  двух конформаций олигопептидов глицина и аланина с длиной от двух до десяти остатков (в ккал/моль)

n- пептид	Gly				Ala
	С5(=180o,=180o)		(=-60o,=-40o)		С5(=180o,=180o)		(=-60o,=-40o)
	Е пол.	Е кул.	Е пол.	Е кул.	Е пол.	Е кул.	Е пол.	Е кул.
2	0,8	0,8	6,7	6,4	3,9	0,0	6,7	6,4
3	0,6	0,6	5,2	5,2	3,8	0,4	5,7	5,1
4	0,7	0,6	3,8	4,0	3,8	0,4	4,1	4,0
5	0,6	0,6	2,9	3,4	3,8	0,3	3,1	3,1
6	0,6	0,6	2,3	2,9	3,8	0,3	3,1	3,3
7	0,6	0,6	1,8	2,3	3,8	0,3	2,4	2,7
8	0,6	0,6	1,4	2,0	3,8	0,3	1,4	2,0
9	0,6	0,6	1,1	1,7	3,8	0,3	1,1	1,8
10	0,6	0,6	0,9	1,5	3,8	0,3	0,9	1,6

В таблице приведены также величины кулоновского вклада в общую потенциальную энергию. Из представленных в табл.9 данных следует, что в олигопептидах аланина от дипептида до пентапептида предпочтительной является полностью вытянутая С5 конформация, а начиная с пентапептида предпочтительной становится -спиральная. В дипептиде глицина -спиральная конформация невыраженная, энергия ее стабилизации растет с увеличением числа остатков, с 6,7 ккал/моль в дипептиде до 0,6 ккал/моль в декапептиде, и судя по тенденции, она становится предпочтительной начиная с одиннадцати пептида.

Заключение. Изменение химической структуры полипептидной цепи белка на один аминокислотный остаток при биосинтезе происходит на порядки медленнее формирования элементов вторичной структуры и компактных состояний в развернутой цепи, а также предполагаемого для сворачивания белков времени. Именно поэтому белки должны сворачиваться котрансляционно, в процессе биосинтеза. Сворачивание белка динамический процесс формирования пространственной структуры постепенно удлиняющегося с N-конца полипептида, происходящий последовательными этапами по мере его роста - изменения химической структуры, в процессе биосинтеза первичной структуры. На каждом этапе этого процесса, протяженность которого не дольше 0,01 сек, формируются, одновременно и сопряженно, вторичная и третичная структуры синтезированного N-концевого сегмента, которые могут изменяться на следующем этапе, а на последнем этапе - нативной структуры полной нативной последовательности полипептидной цепи белка.

Основные выводы

I. Впервые проведено фундаментальное исследование по выяснению реалистических оценок барьеров внутреннего вращения (БВВ) вокруг связей NCa и CaC главной цепи полипептидов.

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

2. Расчетами полной поверхности потенциальной энергии аминокислот глицина и аланина методом ФФВ установлено, что относительные энергии предпочтительной и равновесных конформаций аминокислот и величины потенциальных барьеров между этими конформациями составляют не более 3 и 5 ккал/моль соответственно.

3. Расчетами комплексов амидов и пептидов методом ФФВ установлено, что в белках возможны образования специфических Н-π водородных связей с энергией стабилизации до 2,5 ккал/моль, которые могут сыграть существенную роль в формировании и стабилизации пространственной структуры белков; наличие таких связей в структурах белков было обнаружено впоследствии кристаллографами.

4. Впервые даны реалистические оценки барьеров внутреннего вращения вокруг связей NCa и CaC главной цепи полипептидов. Исследованиями методом ФФВ полной поверхности потенциальной энергии ряда олигопептидов (глицина и аланина от ди- до декапептида, дипептидов серина, валина, фенилаланина и тирозина, серединного пятого и концевого девятого остатков в декапептидах глицина и аланина) установлено, что барьер вращения вокруг связей NCa составляет около 16 ккал/моль - что сопоставимо с транс-цис барьером пептидной связи, и вокруг связей CaC - около 6 ккал/моль. Эти величины значительно выше общепринятых традиционных значений около 0,7 ккал/моль для обоих барьеров.

5. Из-за высоких барьеров вращения вокруг связей NCa и CaC в белках аминокислотные остатки обладают ограниченной способностью к конформационным изменениям, и возможностью находиться лишь в нескольких дискретных конформациях. Маловероятно, что неглициновые остатки могут приобрести конформацию с положительным значением угла вращения вокруг связи NC из-за высоких барьеров вращения вокруг связей NC, так как синтез белка на рибосоме происходит в конформации правой a-спирали. Стабильность нативной структуры и способности к ренатурации белков, избирательное каталитическое действие ферментов, а также эффективность сворачивания белков, скорее всего, обусловлены высокими потенциальными барьерами вращения вокруг связей NCa и CaC полипептидной цепи.

6. Установлено, - из анализа данных банка белковых структур о кристаллических структурах белков (PDB), что изменение функционального состояния, условий кристаллизации, кристаллизационной среды, кристаллической формы, а также мутации белка не приводит к изменению полярности конформации по углу неглициновых остатков. В наилучшим образом расшифрованных белковых структурах неглициновые остатки в положительных конформациях встречаются редко. Наличие в структуре белка заметного числа неглициновых остатков в положительной конформации указывает на присутствие локальных неточностей в расшифровке структуры, что может служить простым критерием оценки качества при расшифровке структуры белков традиционными методами рентгеновской кристаллографии и ЯМР. Эти результаты надежное косвенное экспериментальное подтверждение наличия достаточно высокого БВВ вокруг связей NCa в белках, на основе  которых также может быть создан новый вычислительный подход для дополнительного уточнения белковой структуры.

II Полипептидная цепь белка вряд ли может приобрести состояние случайного клубка из-за высоких БВВ вокруг связей NCa и CaC так как понятие случайный клубок приемлемо, по определению, для свободносочлененных линейных макромолекул, в которых возможны свободные вращения. Белки не приобретают вид случайного клубка в сильно денатурирующих условиях в отличие от традиционных представлений, а обладают значительной упорядоченной структурой - общая закономерность, выявленная из систематизации результатов нового поколения точных экспериментов по денатурации белков. Вследствие этих данных, по-видимому, является недостаточно обоснованным представление процесса сворачивания белка как переход клубок - глобула.

III. Сравнительный анализ известных данных биофизики, физической химии и клеточной биологии в области протеомики свидетельствует о том, что удлинение белковой цепи на один аминокислотный остаток при биосинтезе происходит на порядки медленнее формирования элементов вторичной структуры и компактных состояний в развернутом полипептиде и предполагаемого для сворачивания белков времени. На этом основании, белки должны сворачиваться в процессе биосинтеза - котрансляционно - в соответствии со следующей схемой. Сворачивание белка происходит последовательными этапами, по мере выхода аминокислотных остатков синтезируемой полипептидной цепи из рибосомы. На каждом этапе формируется конкретная конформация N-концевой части цепи, выступающей из рибосомы, которая может измениться на следующем этапе. На всех этапах процесса, формирование структуры N-концевой части, а на последнем этапе - нативной конформации белка происходит по времени не более 0,01 секунды. Предложенный механизм подтверждается рядом биофизических экспериментов, а его универсальность - данными клеточной биологии.

Список основных публикаций

1. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Ермаков Г.Л., Научитель В.В., Соболев В.М. (1984а). Связь-связевые взаимодействия. I. Простое соотношение для оценок связь-связевого взаимодействия. Ж. структ. химии 25, N1, 31-35.

2. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Соболев В.М. (1984б). Связь-связевые взаимодействия. II. Барьеры внутреннего вращения в насыщенных органических молекулах. Ж. структ. химии 25, N1, 36-41.

3. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Соболев В.М. (1984в). Связь-связевые взаимодействия. III. Барьеры внутреннего вращения в ненасыщенных углеводородах. Ж. структ. химии 25, N2, 3-8.

4. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Ермаков Г.Л., Соболев В.М. (1984г). Связь-связевые взаимодействия. IV. Барьеры внутреннего вращения в альдегидах, кетонах, сложных эфирах, органических кислотах, амидах и аминокислотах. Ж. структ. химии 25, N2, 9-13.

5. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Ермаков Г.Л., Научитель В.В., Соболев В.М. (1984д). Связь-связевые взаимодействия. V. Димеры молекул водорода и углеводородов различных классов. Ж. структ. химии 25, N3, 17-22.

6. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Научитель В.В., Соболев В.М. (1986а). Роль электростатических взаимодействий в стабилизации Цспиралей пептидов. Докл. АН СССР 286, 1261-1264.

7. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Лазарев Ю.А., Соболев В.М. (1986б). Н-связи в амидах и пептидах. Докл. АН СССР 287, 211-215.

8. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Гренадер А.К., Ермаков Г.Л. (1986в). Конформационный анализ и корреляционное соотношение структураЦактивность некоторых потенциально активных антиаритмиков лидокаинового ряда. Биофизика 31, 741-746.

9. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Научитель ВВ, Соболев В.М. (1987а). Стабильность Цспиральной структуры олигопептидов. Молекулярная биология 21, 743-749.

10. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Научитель В.В., Соболев В.М. (1987б). Роль Н-водородных связей в стабилизации спиральных структур олигопептидов. Молекулярная биология 21, 1339-1345.

11. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Научитель В.В., Соболев В.М. (1987в). Роль Н-водородных связей и электростатических взаимодействий в стабилизации различных структур пептидов. В сб: Межмолекулярное взаимодействие и конформации молекул. Ред. Петропавлов Н.Н., Зоркий П.М., Пущино, 96-104.

12. Башаров М.А., Волькенштейн М.В., Голованов И.Б. (1988). Влияние боковых радикалов на структуру пептидов. Поверхности потенциальной энергии модельных дипептидов аланина и фенилаланина. Молекулярная биология 22, 1217-1225.

13. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Научитель В.В., Соболев В.М. (1989а). Метод фрагмент-фрагментных взаимодействий. I. Краткое описание и документация метода для оценок межмолекулярных взаимодействий. Ж. обшей химии 59, 489-502.

14. Башаров М.А., Волькенштейн М.В., Голованов И.Б., Научитель В.В., Соболев В.М. (1989б). Метод фрагмент-фрагментных взаимодействий. II. Оценки внутримолекулярных взаимодействий. Ж. обшей химии 59, 503-512.

15. Башаров М.А. (1990). Некоторый взгляд на формирование пространственной структуры белков. Молекулярная биология 24, 1274-1282.

16. Башаров М.А. (1995а). Модифицированный метод фрагмент-фрагментных взаимодействий для изучения внутри- и межмолекулярных взаимодействий. Журнал физической химии 69, 1422-1426.

17. Башаров М.А. (1995б). Поверхности потенциальной энергии пептидов и конформации аминокислотных остатков в белках. Предварительные результаты обработки банка белковых структур. Биофизика 40, 260-273.

18. Башаров М.А. (1996). Роль стартовой конформации в формировании пространственной структуры белка. Отрицательные конформации неглициновых остатков в белковых структурах. Докл.Акад.наук 348, 115-118.

19. Башаров М.А. (1997а). Распределение конформаций аминокислотных остатков, встречающихся в трехмерных структурах белков. Анализ неглициновых остатков, находящихся в УположительныхФ конформациях. Биофизика 42, 753-764.

20. Башаров М.А. (1997б). Роль шаперонинов в сворачивании белков. Новая модель архитектуры строения GroEL/GroES комплекса. Биохимия 62, 489-498.

21. Башаров М.А. (2000а). Концепция посттрансляционного сворачивания белков: насколько она реалистична? Биохимия 65, 1400-1408.

22. Башаров М.А. (2000б). Котрансляционное сворачивание белков. Биохимия 65, 1639-1644.

23. Башаров М.А. (2002). Имеют ли синтетические белки отношение к проблеме сворачивания белков? Биофизика 47, 989-995.

24. Basharov, M.A. (2002). Protein folding: Chemically synthesized proteins. In: Recent Research Developments in Protein Folding, Stability and Design. M. Gromiha, S. Selveraj, Eds. (Research Signpost, Trivandrum, India), pp. 167-175.

25. Basharov, M.A. (2003). Protein folding. J. Cell. Mol. Med. 7, 223-237.

26. Basharov, M.A., Allakhverdiev, S.I. (2008). Expedience of protein folding modeling during progressive elongation of polypeptide chain. The Open Struct.Biology J. 2, 31-32.

27. Basharov, M.A. (2008). Protein folding: Requirement for simulations on the basis of sequential growth of polypeptides. Biochemistry: An Indian J. 2, 63-65.

  Авторефераты по всем темам  >>  Авторефераты по биологии