Автоматическое моделирование пространственной структуры антител и их комплексов с лигандом
Выполнил ученик 10Н класса Рустам Усманов
Научные руководители:
Олег Владимирович Колясников
Владимир Константинович Аржаник
Оглавление
1. Литературный обзор
1.1 Применение
Антитела применяются в разных областях науки - от проверки продуктов питания на гербициды до биоинженерии и микробиологии. Но все же главной отраслью их применения остается медицина. В этом случае их называют терапевтическими антителами
Терапевтические антитела - род биофармацевтических препаратов, которые нашли применение в различных областях медицины: от онкологии до лечения иммунных и воспалительных заболеваний. (Успешное развитие новых терапевтических антител требует понимания механизмов действия содержащих их лекарственных препаратов на организм.) Зная точную структуру специфических антител, можно разработать новые подходы к избавлению людей от этих недугов (взаимодействия антиген-антитело).
1.2 Строение
Антитела имеют характерную для белков структуру. Их первичная структура представлена последовательностью аминокислот (обычно 450-550), а вторичная- в -слоями (их направление указано стрелочками) и б-спиралями. Четвертичная структура напоминает букву Y.
Антитела состоят из двух лёгких цепей и двух тяжелых цепей.
Антитела распознают антигены, связываясь с определённым эпитопом - характерным фрагментом поверхности или линейной аминокислотной цепи антигена.
1. Fab- участок
2. Fc-участок
3. Тяжелая цепь (состоит из VH, CH1, шарнирного участка, CH2 и CH3 участков)
4. Легкая цепь (состоит из VL, CL участков)
5. Место крепления антигена
6. Шарнирный участок
Некоторые иммуноглобулиновые домены состоят из двух тяжелых цепей (красная и синяя) и двух легких цепей (зеленая и желтая) антитела. Иммуноглобулиновые домены состоят из от 7 (для константных доменов) и 9 (для вариабельных доменов) в-стрэндов.
Переменная часть антитела являются его V-регионом, а постоянная часть его С-регионом.
Некоторые части антитела имеют одни и те же функции. N-концевая часть антитела, например, содержит сайты, которые могут связывать два антигена. Эта область антитела называется Fab (fragment, antigen binding) region (фрагмент, связывающий антиген). Он состоит из одного постоянного и одного вариабельного домена от, соответственно, тяжелой и легкой цепи антитела.
С-концевая часть Fab играет большую роль и в регуляции активности иммунных клеток. За это отвечает область под названием Fc (кристаллизующийся фрагмент). Она состоит из двух тяжелых цепей, образующих два или три (в зависимости от класса антитела) константных домена. Таким образом, область Fc гарантирует, что каждое антитело генерирует соответствующий иммунный ответ на данный антиген, путем связывания с конкретным классом Fc рецепторов. Делая это, он является посредником различных физиологических эффектов, включая опознавание опсонированных (подлежащих уничтожению, лизису) частиц, лизис и очищение от токсинов жировых клеток, базофилов и эозинофилов.
Вариабельный домен (также упоминается как область Fv) является наиболее важным регионом для связывания с антигенами. В частности, вариабельные петли, три на каркасном участке легкой (VL) и три - на участке тяжелой (VH) цепи ответственны за связывание с антигеном. Эти петли называются областями определяющими комплементарность (CDR). Стоит заметить, что такое огромное разнообразие антител объясняется лишь различностью их антиген-связывающего участка, называемого эпитопом. Основные его части называются CDR - Complementary Determining Regions). Участки тяжелых и легких цепей, содержащих CDR, называют вариабельными.
1.3 Получение структур антител
Экспериментальные методы получения структуры антител, такие как рентгеновская кристаллография и ядерный магнитный резонанс, являются энергоемкими, трудоемкими и дорогостоящими. Полученные структуры заносятся в базу данных белков в специальном формате PDB (Protein DataBase) с координатами каждого атома молекулы. В то время как вычислительный метод предсказания структуры антител (т.е. их компьютерное моделирование) обеспечивает быстрый и недорогой путь для их получения структур. Хотел бы также подчеркнуть, что структуры некоторых антител не могут быть получены иным образом. Наиболее распространенный метод - моделирование по гомологии (автоматическое моделирование).
Одной из программ для работы с такими файлами является PYMOL, написанная на языке Python. Используется трехмерная декартова система координат, в которой положение каждого атома определяется его X, Y, Z координатами.
Остановимся подробнее на одном из этих серверов:
В Rosetta Antibody сервер (http://antibody.graylab.jhu.edu) прогнозирует пространственную структуру первичной последовательности антител в CDR, учитывая аминокислотные последовательности, соответствующие легким и тяжелым цепям. Сначала, сервер определяет и отображает наиболее похожие на специальный шаблон структуры (в отношении CDR-петель). Впоследствии наиболее схожие из них собираются в оптимизированную "грубую" модель, и сервер отображает эскиз структуры антитела и файл содержащий координаты трехмерной модели в пространстве.
Rosetta Antibody генерирует 2000 независимых структур, и сервер отображает изображения, файлы координат, и подробную скоринговую информацию для 10 самых правдоподобных моделей. Они позволяют пользователям использовать рациональное суждение в выборе оптимальной для них модели или использовать комплект как ансамбль для дальнейших исследований, таких, как докинг.
Следует отметить, что существующие серверные приложения не обеспечивают высокого разрешения уточнения структуры антитела и не учитывает термодинамических характеристик в процессе моделирования.
Итак, для пользователя работа описанных приложений выглядит так: вы вводите в специальную форму только аминокислотные последовательности легких и тяжелых цепей, имя антитела, имя пользователя, а на выходе получаете файл, содержащий:
1) Дату и время подачи документов, оформления и завершения.
2) Диаграммы показывающие границы CDR петель.
3) Таблицу, показывающую последовательности, наиболее гомологичные шаблону (VH-VL структур, CDR, L1, L2, L3, H1, H2 и H3).
4) Пространственную структуру(ы), приписываемую данному антителу.
Rosetta также разработала новый алгоритм "KIC", позволяющий предсказывать структуры антител на порядок аккуратнее (по слепому сравнению с эталоном - кристаллографической структурой) чем это делает устоявшийся алгоритм ROSIE.
Материалы и методы: Нами был использован сайт http://rosettadock.graylab.jhu.edu/ для моделирования пространственной модели структуры обоих вариабельных доменов 5 антител по гомологии. При этом исходным материалом для меня служили последовательности аминокислот первичной структуры из статьи Франека [1].
2. Результаты и обсуждение
антиген антитело генетический алгоритм
Для визуализации полученных моделей я использовал программу PyMol 1.5.0.3. на языке программирования Python.2.7.6. Были обнаружены сайты связывания, образованные всегда CDR-петлями тяжелых и легких цепей антитела.
Для удобства, на иллюстрациях атомы углерода тяжелой цепи выделены зеленым, легкой - голубым, кислорода и азота обеих цепей - синим и красным, а атомы сайта связывания - фиолетовым. Аминокислоты боковой цепи CDR петлей, непосредственно учавствующие в связывании лиганда обозначены палочками.
Ниже приведены иллюстрации сайтов связывания 5 антител, выполненные в PyMol:
1. E2/G2
|
Аминокислоты, формирующие сайт связывания |
Местоположения в антителе |
Характеристика сайта связывания |
|
|
Tyr H 95 |
CDR-H3 |
Сайт связывания представлен практически плоской поверхностью; сформирован совокупностью небольших гидрофобных и гидрофильных участков. |
|
|
Tyr L 91 |
CDR-L3 |
||
|
Tyr L 96 |
CDR-L3 |
1. 2. E2/B5
|
Аминокислоты, формирующие сайт связывания |
Местоположения в антителе |
Характеристика сайта связывания |
|
|
Val H 33 |
CDR-Н 1 |
На поверхности есть большие гидрофобные участки, однако, форма сайта, в целом, - плоская |
|
|
Tyr H 52 |
CDR-H2 |
||
|
Tyr H 95 |
CDR-H3 |
||
|
Tyr H 98 |
CDR-H3 |
||
|
Ser L 32 |
CDR-L1 |
||
|
His L 33 |
CDR-L1 |
||
|
Tyr L 91 |
CDR-L3 |
||
|
Tyr L 96 |
CDR-L3 |
3. E4/C2
|
Аминокислоты, формирующие сайт связывания |
Местоположения в антителе |
Характеристика сайта связывания |
|
|
Phe H 32 |
CDR-H1 |
Полость образована гидрофобными участками Tyr, Val, Leu, Ile, Ser, формирующие гидрофобное ядро. Небольшая часть полости образована гидрофильным азотом боковой цепи триптофана H47. Глутаминовая кислота и гидрофильные участки боковых цепей (тирозины L96, L91) формирует гидрофильный вход в полость. |
|
|
Ile H 33 |
CDR-H1 |
||
|
Tyr H 95 |
CDR-H3 |
||
|
Gly H 96 |
CDR-H3 |
||
|
Tyr H 97 |
CDR-H3 |
||
|
Ser L 32 |
CDR-L1 |
||
|
His L 43 |
Каркасные участки вар. домена |
||
|
Tyr L 91 |
CDR-L3 |
||
|
Tyr L 94 |
CDR-L3 |
||
|
Phe L 96 |
CDR-L3 |
4. B5/C3
|
Аминокислоты, формирующие сайт связывания |
Местоположения в антителе |
Характеристика сайта связывания |
|
|
Phe H 32 |
CDR-H1 |
С виду сайт связывания плоский, но оказалось, что внутри имеется полость. За счет теплового движения боковой цепи аминокислоты могут менять свое положение, освобождая путь к внутрь полости, однако, получившийся сайт связывания говорит о некоторых неточностях в моделировании комплекса с лигандом. Полость сформирована большими гидрофобными участками аминокислот Tyr, Phe, His. Также имеется 2 участка, сформированных положительным электрическим зарядом аргинина и отрицательным зарядом глютаминовой кислоты. |
|
|
Ile H 33 |
CDR-H1 |
||
|
Tyr H 95 |
CDR-H3 |
||
|
Tyr H 97 |
CDR-H3 |
||
|
His L 34 |
CDR-L1 |
||
|
Tyr L 91 |
CDR-L3 |
||
|
Tyr L 94 |
CDR-L3 |
||
|
Phe L 96 |
CDR-L3 |
1. 5. B7
|
Аминокислоты, формирующие сайт связывания |
Местоположения в антителе |
Характеристика сайта связывания |
|
|
Tyr H 32 |
CDR-H1 |
Полость образована гидрофобными участками Tyr, Val, Leu, формирующие гидрофобное ядро. Небольшая часть полости образована гидрофильной боковой цепью серина H 35 и азотом боковой цепи триптофана H 47. Глутаминовая кислота и гидрофильные участки боковых цепей: тирозин L 96, H52 и кислород тирозина H95 формирует гидрофильный вход в полость. |
|
|
Val H 33 |
CDR-H1 |
||
|
Ser H 35 |
CDR-H1 |
||
|
Val H 37 |
CDR-H1 |
||
|
Trp H 47 |
Каркасные участки вар. домена |
||
|
Glu H 50 |
CDR-H2 |
||
|
Tyr H 52 |
CDR-H3 |
||
|
Tyr H 95 |
CDR-H3 |
||
|
Tyr H 97 |
CDR-H3 |
||
|
His L 34 |
CDR-L1 |
||
|
Tyr L 49 |
CDR-L2 |
||
|
Leu L 89 |
CDR-L3 |
||
|
Tyr L 91 |
CDR-L3 |
||
|
Tyr L 96 |
CDR-L3 |
||
|
Phe L 98 |
CDR-L3 |