Скачать презентацию
Идет загрузка презентации. Пожалуйста, подождите
Презентация была опубликована 11 лет назад пользователемscc.acad.bg
1 Сравнительный анализ пространственных структур белков 1. Оценка качества пространственной структуры (по G.Kleywegt, 2000)
2 Построение пространственной модели белка в результате рентгеноструктурного анализа кристалла не есть точная наука!
3 Модель зависит от ее авторов: их опыта самокритичности мотивации используемого программного обеспечения Ошибки в модели почти неизбежны!
4 При высоком разрешении (
5 Остальные 5% модели зависят от того, какие программы использовались как моделировался B-фактор (есть свобода выбора) допускались ли альтернативные конформации моделировались ли водороды
6 Продолжение какие сгущения электронной плотности интерпретировались как молекулы растворителя (воды), а какие – как шум как отнеслись к некристаллографической симметрии: посчитали ли 2е или более копий одной и той же молекулы белка в ассиметрической ячейке –идентичными или –разрешили им отличаться?
7 При разрешении >2 Å (>90% РСА моделей в PDB) роль воображения авторов модели возрастает!
8 Задачи, решаемые пользователем PDB Выбор хорошей (лучшей) модели из нескольких моделей одного и того же или нескольких родственных белков Интегральная оценка качества и, следовательно, степень доверия модели Выявление маргинальных групп атомов Маргинальные атомы в модели могут быть либо ошибками расшифровки (отщепенцами), либо особенными атомами белка (белыми воронами)
9 Основные индикаторы маргиналов 1.Пространственный R-фактор (Jones et al., 1991) 2.Карта Рамачандрана (Ramachandran&Ramachandran, 1965) 3.Ротамеровость боковой цепи (Jones et al., 1991) 4.Комфортность окружения атомов (Vriend&Sanders, 1993), в частности, анализ водородных связей 5.Инверсии петидной цепи от C α до следующего C α (pep-flip; Jones et al., 1991, Kleywegt&Jones, 1998)
10 высокий B-фактор ( >40) странный коэффициент заполнения (напр. 0!?) длины связей углы торсионные углы (не) планарность (неправильная) хиральность (С α, для Thr, Ile – еще и С β ) Другие индикаторы маргиналов
11 Продолжение недопустимое сближение несвязанных атомов существенное различие B-факторов ковалентно связанных атомов существенное различие конформаций в разных копиях одного и того же белка из ассиметрической ячейки
12 Пространственный R-фактор (Real Space R, RSR) Характеризует насколько модель атома (-ов) соответствует экспериментальной электронной плотности Сумма берется по узлам пространственной решетки в окружении атома (атомов) Хорошие значения: RSR20%
13 Продолжение ρ эксп восстанавливается по фазам из модели и интенсивностям рефлексов из эксперимента (часто берут взвешенные интенсивности: эксперимент+модель – отсюда разные карты электронной плотности для одной модели) Карты электронной плотности моделей, для которых в PDB есть файл структурных факторов, доступны на сайте Electron Density Server (EDS)
14 Вариант RSR – коэфициент корреляции между ρ эксп и ρ модели Как и RSR, вычисляется по узлам пространственной решетки Не зависит от значений ρ эксп, а зависит от согласованности изменений ρ эксп и ρ модели Coeff.corr.
15 Еще один вариант RSR – Z-score Для вычисления Z остатка (напр. Ala57) его RSR сравнивается со средним RSR для того же типа остатков (Ala) по выборке из PDB с примерно таким же разрешением (напр Å)
16 Карта Рамачандрана характеризует насколько конформация плипептидной цепи остатка (псевдо- )энергетически выгодна -180° 0° +180° φ ψ +180° 0° -180°
17 Концепция псевдоэнергии Идея: то, что наблюдается чаще – более энергетически выгодно Реализация: –строим гистограмму интересующего нас параметра по подходящей выборке из PDB –приближаем ее разумной функцией (например, гауссовой) –считаем, что полученный график отражает энергию в зависимости от параметра
18 81782 остатка на карте Рамачандрана Beta (B) Left helix (L) Alpha(A)
19 Области на карте Рамачандрана В хорошей модели >90% остатков, не считая Gly, Pro, находятся в предпочитаемой области Классификация областей (PROCHECK): - предпочитаемая (A,B,L) - разрешенная (a,b,l,p) - допустимая (~a,~b,~l,~p) - запрещенная 1CNR, разрешение 1.05
20 Карта какого остатка? PROLINE
21 Где Trp, где Ser? TrpSer
22 Анализ водородных связей В моделях встречается инверсия боковых цепей His, Asn, Gln Asn His Gln
23 Инверсия (?) в Asn51 гомеодомена 2 Атом ND2 Атом OD1 Атом N6 – донор H Атом N7 – акцептор H ? ? 1 2 A103 Asn51 Еще в 36 структурах гомеодоменов – так же, как в 1й; еще в 2х – как во второй
24 Анализ молекул воды. Пример из модели 1CBS HOH376 HOH375 Leu28.CD 3.9Å 4.3Å Может ли HOH376 фиксироваться в одинаковых точках во всех ячейках кристалла!??? Думаю, что нет
25 Интегральная оценка комфортности окружения остатка В программе WhatCheck рассчитывается Z-score для комфортности окружения каждой боковой цепи Маргиналы – Z-score < -5 Более показательны участки цепи с низким Z, для их обнаружения строится сглаженный график зависимости Z от номера остатка
26 Статистика конформаций боковых цепей остатков Боковые цепи имеют от 0 (Gly, Ala) до 4х (Lys, Arg) степеней свободы вращения вокруг ковалентных связей Обозначаются χ 1,..., χ 4
27 Распределение угла χ 1 в моделях PDB Всего в выборке было остатков
28 Карта углов χ 1 и χ 2 для Leu Ось X: chi_1 (0-360°) Ось Y: chi_2 (0-360°) Всего 6638 остатков Изображены линии уровня плотности числа остатков
29 Ротамеры - это боковые цепи в типичных для данного типа остатков конформациях У каждого типа остатков (Leu, Trp, Arg и т.п.) свое число ротамеров SwissPDBviewer знает ротамеры и может их показать
30 Примеры грубых ошибок в записях PDB Кристалл не того белка !!! Полностью ошибочная модель (!) Правильно найдены, но неправильно соединены элементы вторичной структуры Участок цепи вписан в электронную плотность с ошибкой
31 Janin, 1990: в лаборатории очищали и кристаллизовали ретинат-связывающий белок; при РСА обнаружилось, что получился лизоцим …
32 Две структуры фотоактивного желтого белка рецептора фототаксиса Ectothiorhdospira halophila: почти ничего общего! 1PHY, 2.4 Å McRee et al PHY, 1.4 Å Borgstahl et al. 1995
33 Хлормуконат циклоизомераза из Alcaligenes eutrophus: совмещенные структуры 1chr и 2chr 1CHR: Hoier et al., 1993 Разрешение 3.00 Å 2CHR: Kleywegt et al., 1996 Разрешение 3.00 Å
34 Наложение 1CHR_A и 2СHR
35 Gln20 Ser22 Ala5 Lys16 Gly40 1CHR и 2CHR: пример сдвига рамки при расшифровке
36 Выравнивание последовательностей 1CHR и 2CHR по близости C_alpha атомов при наложении структур
37 Мог ли пользователь заподозрить ошибку в структуре 1CHR до появления 2CHR?
38 Server EDS: Real-space R-value vs Residue for 1chr
39 RSR: пространственный R-фактор для всех остатков структуры 1CHR
40 1СHR: коэффициент корреляции для RSR
41 То же – для всех остатков
42 EDS server: Z-score vs Residue for 1chr Z-score=(RSR- )/sigma
43 Z-score для всех остатков структуры 1CHR
44 Real-space R-value vs Residue for 2chr
45 RSR для всех оснований из 2chr
46 2СHR: коэффициент корреляции для RSR
47 То же – для всех остатков
48 EDS server: Z-score vs Residue for 2chr
49 Z-score для RSR для всех остатков структуры 2CHR
50 Packing environmen for 1CHR: Second generation quality Z-score plot smoothed over a 10 residue window (PDBCheck)
51 1CHR: Arg35 оказался в гидрофобном кармане 2CHR: Arg35 образует одну водородную связь
52 Наложение участка 5-47 двух цепочек A и B из структуры 1CHR
55 1CHR 2CHR
56 Сравнение структур 1CHR и 2CHR Statistic1CHR 2CHRComments Crystallo graphic R-value Not very informative Free R- value ???0.264Cannot compare
57 Продолжение Statistic1CHR 2CHRComments WHAT IF Z-score Ramachandran plot appearance (for 3.0 Å structures) CHR worse than average WHAT IF Z-score rotamer quality (for 3.0 Å structures) CHR worse than average; 2CHR better than average WHAT IF Z-score backbone conformation quality (for 3.0 Å structures) CHR worse than average; 2CHR better than average
58 Сервисы и программы PDBsum PDB PDBCheck (программа WhatCheck из пакета WhatIf) PDBReport EDS (RSR во всех видах, файлы с электронными плотностями)
59 Продолжение Ramachandran plots for all types of residues ( Rotamers for all types of side chains (
Еще похожие презентации в нашем архиве:
© 2024 MyShared Inc.
All rights reserved.