Файл: Дашевский, В. Г. Конформации органических молекул.pdf

ных олигопептидов, является проблема многих минимумов, и в этом смысле олигопептиды являются «маленькими белками». Дей­ ствительно, если каждый аминокислотный остаток («пептидная единица») может иметь три конформации R, В и L, то в декапепти­ де возможно уже З10 конформаций. В циклических пептидах бла­ годаря условиям замыкания число возможностей несколько мень­ ше; гораздо меньше оно для симметричных циклических пептидов, и потому в ряде случаев поиск оптимальной структуры становится технически выполнимым.

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

В работе [102] было экспериментально установлено, что молекула этого «шестизвенного» цикла может иметь, в зависимости от поляр­ ности среды, две конформации — «неполярную» N и «полярную» Р. Расчеты конформаций однозначно предсказывают эти две струк­ туры и помогают объяснить большой экспериментальный материал по этому соединению (в частности, спектры ЯМР).

Конформационные карты Ы-ацетил-К-метил-Двалина и диметиламида О-ацетил-ДДоксиизовалериановой кислоты, модели­ рующих два фрагмента энниатина В, приведены на рис. 8.12 а и б соответственно. Как видим, первый фрагмент может иметь четы­ ре конформации приблизительно равной энергии, тогда как вто­ рой фрагмент — только одну. Для поиска оптимальных конфор­ маций всей молекулы энниатина В вовсе не обязательно миними­ зировать потенциальную функцию по многим переменным — достаточно собрать скелетные модели с разрешенными конформация­ ми (комбинации из четырех конформаций первого фрагмента и од­ ной конформации второго фрагмента). Рассмотрение моделей пока­ зывает, что молекула действительно может иметь две конформа­ ции, если шестизвенный цикл замкнут. Одна из них — «поляр­ ная» форма — показана ниже вместе с ионом К+, входящим в по­ лость молекулы и образующим координационные связи с атомами кислорода.

384

О— атомы углерода 0 —атомы кислорода

Если ионы калия отсутствуют, молекула энниатина В перехо­ дит в другую конформацию — «неполярную». Таким образом, рас­ четы и экспериментальные данные (на последних мы не остана­ вливаемся) показывают, что в основе механизма переноса ионов этими молекулами лежат конформационные перестройки.

Рис. 8.12. Конформационные карты Ы-ацетил-Ы-метил-/.-валина (а) и диметиламида 0-ацетил-0£-оксиизовалериановой кислоты (б) — моделей фраг­ ментов молекулы энниатина В

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

L-Val—L-Orn—L-Leu—D-Phe —L-Pro

L-Pro—D-Phe—L-Leu—L-Om — L-Val

Рентгеноструктурное исследование этой молекулы [103] пока­ зало наличие оси симметрии второго порядка. Хотя полностью расшифровать структуру не удалось, авторы этого исследования предложили пространственную модель грамицидина С. Десятиле­ тием позже Ликвори [104] и Шерага 113, 14] нашли пространствен­ ные модели, удовлетворяющие критерию допустимых контактов, причем Шерага делал выбор из 282 различных конформаций, ми­ нимизируя потенциальную функцию. Однако все три предложен­ ные структуры не удовлетворяли некоторым косвенным экспери­ ментальным данным [105]. В связи с этим было рассмотрено еще несколько моделей пространственной структуры грамицидина С [106, 107J. И все же до настоящего времени конформацию этой мо­ лекулы пока еще нельзя считать окончательно выясненной, хотя наиболее вероятной представляется структура, стабилизируемая двумя водородными связями между группами С = 0 и N—Н остат­ ков валина и лейцина [1081.

Итак, для предсказания конформации нерегулярного олиго­ пептида требуется найти абсолютный минимум функции многих переменных, что само по себе является весьма трудной задачей (см. раздел 5 гл. 2). Далее, надо иметь в виду, что в выборе опти­ мальной структуры играет роль не только глубина ямы, но и ее форма. В самом деле, оптимальная структура определяется ми­ нимумом свободной энергии, а не энтальпии, а форма ямы дает энтропийную часть свободной энергии.

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

4. ПРОСТРАНСТВЕННАЯ СТРУКТУРА ГЛОБУЛЯРНЫХ БЕЛКОВ

Конформации дипептидных фрагментов

вбелках

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

386

добный анализ был проведен Брантом и Шиммелем [109], Рамачандраном 116], Шерага [17] и другими авторами.

На рис. 8.13 показаны конформационные карты метиламида N-ацетилаланина с нанесенными на них точками, соответствую­ щими углам вращения в миоглобине кашалота [ПО, 111] и лизо­ циме яичного белка [112, 113]. Понятно, что для подобного анали­ за аланин является наиболее подходящим остатком, поскольку он моделирует все другие остатки, содержащие атом Ср. Карты

Рис. 8.13. Конформационные карты метиламида Ы-ацетил-/.-аланина с на­ несенными на них точками, соответствующими углам <р и ф в неспиральных участках миглобина кашалота (а) и лизоцима яичного белка (б). Кружки соответствуют остаткам глицина.

на рис. 8. 13 построены с потенциалами Дашевского без учета элек­ тростатических взаимодействий и водородных связей. Действи­ тельно, эти взаимодействия несущественны для понимания кон­ формаций дипептидных фрагментов в белках: электростатические взаимодействия сказываются лишь на относительной стабильно­ сти конформеров и в значительной степени зависят от природы бо­ кового радикала R, а водородные связи, стабилизирующие свер­ нутые формы дипептидов, встречаются в глобулярных белках край­ не редко; они возможны лишь в участках, ответственных за пово­ роты полипептидной цепи. Наконец, следует заметить, что точки, соответствующие спиральным участкам белков, очевидно, не имеет смысла наносить на карты, поскольку все они должны оказаться в благоприятных областях.

Как видно из рис. 8.13а, все точки миоглобина попадают в об­ ласти низких энергий — до 2 ккал/моль. В лизоциме, как было указано в работе [108], 18 остатков оказываются в неблагоприят­ ных областях. Из них пять остатков глицина (кружки на

рис. 8.136), которые на карте метиламида N-ацетилглицина на­ ходились бы в пределах контура 1 ккал/моль. Далее, конформа­ ции большого числа остатков попадают в область перешейка меж­ ду формами R и В, т. е. эти конформации проигрывают оптималь­ ным примерно 1,5 ккал/моль. Однако конформационный анализ трипептидных фрагментов П14, 1151 показывает, что подобные конформации могут стабилизироваться за счет образования водо­ родных связей между первым и третьим пептидными звеньями. При этом вместо семичленных «циклических» структур, характер­ ных для свернутых форм дипептидов, здесь уже возникают деся­ тичленные «циклы». Но если семичленные «циклы» в белках крайне редки, то трипептидные десятичленные циклы удовлетво­ ряют всем стереохимическим условиям и действительно часто наблюдаются.

Таким образом, конформационная карта метиламида N-ацетил- аланина превосходно объясняет распределение точек (ф, ф) в глобулярных белках. Однако конформационный анализ таких ма­ лых фрагментов, как дипептиды, полезен для объяснения конфор­ маций пептидных звеньев в белках, но явно недостаточен для их предсказания. Теперь мы обратим внимание на те соображения и конформационные расчеты, которые могут быть положены в осно­ ву предсказаний.

Уже давно высказывалось мнение, что первичная структура белка — аминокислотный состав — определяет его пространст­ венную структуру [116, 117). Поскольку за последнее десятилетие структура пяти белков — гемоглобина [1181, миоглобина [ПО, 1111, лизоцима [112, 113], а-химотрипсина [119] и рибонуклеазы А [120] — стала известна детально, а еще для нескольких белков рентгеноструктурный анализ дал достаточное разрешение для то­ го, чтобы увидеть основные черты вторичной и третичной струк­ тур, многие исследователи пытаются найти корреляцию между аминокислотным кодом (составом и последовательностью амино­ кислотных остатков) и конформационным кодом (пространствен­ ной структурой).

Поиски корреляций связаны с тем обстоятельством, что полное предсказание пространственной структуры белка по его амино­ кислотной последовательности сопряжено с неопределимыми в на­ стоящее время трудностями. Поэтому представляется полезным выявить некоторые детали пространственной структуры белка, ос­ нованные на его аминокислотном коде. Большой интерес пред­ ставляет теоретическая оценка процента спиральности, а также выделение спиральных и неспиральных участков. Можно надеять­ ся, что в будущем такой подход приведет к предсказанию конформационной структуры белков с большим содержанием а- или |3- спиральных участков. Тогда задачу предсказания структуры бел­ ка можно разбить на два этапа: на первом этапе определяются спиральные и неспиральные участки, а затем подбираются такие взаимные расположения спиральных и неспиральных участков.

388

которые дают минимум энергии для белка со сформировавшейся вторичной структурой.

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

Статистический подход

В некоторых работах [121—125] на основании результатов экспе­ римента или статистической обработки данных по простран­ ственной структуре белков были выделены «спиральные» и «антиспиральные» аминокислотные остатки. Например, в работе [121] было показано, что остатки пролина, рассеянные случайным образом по полипептиду, полностью разрушают спираль, если со­ держание пролина достигнет 8%; в работе [122] на основании экспериментальных данных были выделены антиспиральные остат­ ки — Ser, Tre, Val, Не и Cys. В работе 11241из статистических рас­ четов миоглобина кашалота, а- и [5-цепей гемоглобина лошади и лизоцима яичного белка были выделены в качестве спиральных остатков остатки Ala, Leu, Val, а в качестве антцспиральны-х — остатки Arg, Asp—NH2, Pro. О. Б. Птицын 1125] провел статисти­ ческие расчеты по 6 белкам, т. е. приблизительно по 1000 амино­ кислотным остаткам: помимо тех белков, которые рассматривал Кук, были включены рибонуклеаза и а-химотрипсин. В резуль­ тате были выделены спиральные остатки Leu, Ala, Glu и антиспи­ ральные Tre, Asp—NH2, причем «выбрасывание» из расчета двух гомологичных белков — глобинов — не повлияло на результаты.

Имеется ряд работ, в которых были введены более сложные критерии для нахождения спиральных и неспиральных участков в белках. Итальянские авторы [126] вводят для каждого остатка (учитывая особенности окружения) так называемые спиральные потенциалы. Статистика, учитывающая способность как отдель­ ных остатков, так и их комбинаций образовывать спирали, разви­ та в работах [127—130]. Например, Протеро [127] выдвинул пра­ вило, согласно которому любой участок нолипептидной цепи из пяти аминокислот будет спиральным, если три остатка в нем —

Ala, Val, Leu или Glu.

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

389