PDB, Jmol, PyMol�и работа с трехмерными�структурами
План
2
PDB файл
3
HEADER ELECTRON TRANSFER 19-APR-98 1BAW
TITLE PLASTOCYANIN FROM PHORMIDIUM LAMINOSUM
…
ATOM 1 N GLU A 1 35.039 18.146 24.600 1.00 73.96 N
ATOM 2 CA GLU A 1 34.182 18.513 25.759 1.00 73.96 C
ATOM 3 C GLU A 1 33.126 19.536 25.368 1.00 73.96 C
ATOM 4 O GLU A 1 33.124 20.038 24.244 1.00 87.91 O
ATOM 5 CB GLU A 1 35.048 19.052 26.891 1.00 87.91 C
ATOM 6 CG GLU A 1 36.075 18.056 27.376 1.00 87.91 C
ATOM 7 CD GLU A 1 37.009 18.647 28.405 1.00 87.91 C
ATOM 8 OE1 GLU A 1 36.536 19.012 29.502 1.00 87.91 O
ATOM 9 OE2 GLU A 1 38.220 18.744 28.120 1.00 87.91 O
ATOM 10 N THR A 2 32.242 19.858 26.307 1.00 27.58 N
ATOM 11 CA THR A 2 31.164 20.807 26.063 1.00 27.58 C
ATOM 12 C THR A 2 30.988 21.697 27.283 1.00 27.58 C
ATOM 13 O THR A 2 30.583 21.234 28.347 1.00 31.52 O
ATOM 14 CB THR A 2 29.834 20.077 25.761 1.00 31.52 C
ATOM 15 OG1 THR A 2 29.997 19.235 24.612 1.00 31.52 O
ATOM 16 CG2 THR A 2 28.726 21.078 25.487 1.00 31.52 C
ATOM 17 N PHE A 3 31.345 22.968 27.128 1.00 27.51 N
ATOM 18 CA PHE A 3 31.235 23.948 28.204 1.00 27.51 C
ATOM 19 C PHE A 3 29.925 24.713 28.083 1.00 27.51 C
Источники данных
4
5
Сайт Protein Data Bank https://www.rcsb.org/
6
Advanced – хороший поиск …
Advanced search
7
.. не только по ключевым словам, но и по наличию лигандов,
методу решения структуры и многому другому
Смотрим на сайте PDB!
Страница одной структуры
8
Structural viewers
9
Structural viewers
10
Им вообще все эти программы надо изучить :)
Про установку PyMol под Windows читайте https://pymolwiki.org/index.php/Windows_Install
PyMol (красивый и умелый)
11
PyMol (красивый и умелый)
12
Jmol/JSmol
13
Далее описан язык, которым надо пользоваться для подготовки изображений в Jmol/JSmol.
Язык �Jmol�(введение)
14
Мы видим картинку,
а Jmol знает про
множества атомов
Язык �Jmol�(load)
15
Формат команды
load <имя файла>
load =<PDB-код>
Назначение
Загрузить новый файл или скачать запись из PDB
Примечания
Язык �Jmol�
16
Все остальные команды
могут сделать что-то,�но никто, кроме вас, не может знать,�с какой частью молекулы это делать
поэтому используйте команды
restrict и select
Язык �Jmol�(select)
17
Формат команды
select <выражение>
Назначение
Назначить новое выделение
Примечания
Язык �Jmol�(простые выражения)
18
at | acidic | acyclic | aliphatic |
alpha | amino | aromatic | backbone |
basic | bonded | buried | cg |
charged | cyclic | cystine | helix |
hetero | hydrophobic | ions | large |
ligand | medium | neutral | nucleic |
polar | protein | purine | pyrimidine |
selected | sheet | sidechain | small |
solvent | surface | turn | water |
Язык �Jmol�(restrict)
19
Формат команды
restrict <выражение>
Назначение
Назначить новое выделение и спрятать изображение всего остального
Примечания
Язык �Jmol�(что мы видим)
20
Формат команды
<cпособ> [off|<число>|only]
Назначение
Нарисовать выбранным способом
Примеры способов изображения
21
wireframe
backbone
trace
cpk
spacefill
cartoons
strands
Язык �Jmol�(способы изображения)
Язык �Jmol�(color)
22
Формат команды
color <цвет>
Назначение
Покрасить выделенное в нужный цвет
Что такое <цвет>
23
Примечания
Немного об RGB цветах
Язык �Jmol�(color)
Язык �Jmol�(center)
24
Формат команды
center <выражение>
Назначение
Центрирует изображение по заданному множеству атомов
Можно, например, писать�center selected
или�center :a
Язык �Jmol�(zoom)
25
Формат команды
zoom <размер>
Назначение
Изменяет размер
Примечание: изменяется размер всего изображения, �а не только выделенного фрагмента
Язык �Jmol�(script)
26
Формат команды
script <имя файла>
Назначение
Выполнить команды, записанные в файле.
Примечания
Язык �Jmol�(пример)
27
load C:\PDB\1ncx.pdb — загрузить файл
restrict none — удалить всё изображение
select all — выделить всё (мы указываем программе, что далее будем работать со всеми атомами из файла)�1433 atoms selected — ответ программы (но изображения пока нет!)
cartoons — изобразить в виде ленточной модели
color group — раскрасить по номеру остатка
Язык �Jmol�(атомы)
28
Atom expressions
ser25:a.ca – Cα атом серина 25 цепочки A�25:a – все атомы остатка номер 25 цепочки A�25 – все атомы из остатков номер 25 всех цепочек �:a – все атомы цепочки A�ser – все атомы всех серинов
Можно использовать знаки * и ?
ser25.c? – все углероды 25-х серинов �ser*:a.ca – Cα атомы всех серинов цепочки :a
Язык �Jmol�(and,or,not)
29
Любые части выражений
можно комбинировать с помощью
логических операций и скобок
Все азоты, не входящие в состав белка
nitrogen and not protein
Все кислороды, входящие в состав белка �или нуклеиновых кислот
(protein or nucleic) and oxygen
:a and :c� ничего выделено не будет
Аденины в цепочках A и С
(:a or :c) and a
Язык �Jmol�(within)
30
Формат выражения
within(<расстояние>,выражение)
Назначение
В командах select, restrict и define определяет атомы, расположенные на расстоянии не более заданного от атомов выражения.
Примечание: расстояния надо указывать, используя точку (не «5», а «5.0»)
31
Примеры
Выделить атомы, расположенные
на расстоянии не более 3Å от DNA
select within(3.0,dna)
на расстоянии не более 2,5Å от DNA
select within(2.5,dna) and not dna
на расстоянии не более 3Å от белка (только атомы DNA)
select within(3.0,protein) and dna
остатки, имеющие хотя бы один атом
на расстоянии не более 2,5Å от DNA
define nears within(2.5,dna) and not dna
select within(group, nears)
Язык �Jmol�(within)
32
Формат команды
define <имя> <выражение>
Назначение
Назвать группу атомов своим именем
Зачем это надо?
Часто бывает необходимо построить очень
сложное выражение. Обычно его можно написать в одну
строчку, но неудобно. Кроме того, если потребуется его
использовать еще раз – то define позволяет просто
сослаться на ранее определенное имя
Язык �Jmol�(define)
Язык �Jmol�(пример)
33
load c:\PDB\1apl.pdb
restrict none
define needdna dna and within(18.0,:c)
define need needdna or :c
select need
center selected
select :c and backbone
wireframe 80
select needdna
wireframe
select needdna and C and not backbone
wireframe 80
wireframe
define cw :c and within(3.0,needdna)
select cw
select cw and not backbone
cpk 150
cpk off
define mainc :c and (136,175,182,184-186)
select mainc
wireframe 120
select needdna and within(5.0,mainc)
wireframe 120
Этот пример надо набить в�Jmol’е и объяснить
�И про PyMol
34