Исследование вторичной структуры белка HutP_BACSU (PDB ID: 1WPV)

На всякий случай: 1WPV.pdb.

Общее описание вторичной структуры белка HutP_BACSU.

Изображение цепи A, полученное с помощью команд: ribbons 150, colour structure. α-спирали окрашены розовым, β-листы - желтым, повороты - видимо, серым (по Manual должны быть светло-голубыми), а остальные остатки - белым.

Данные, полученные командой structure:

Кол-во спиралей Кол-во β-тяжей Кол-во реверсивных повоторов (β-поворотов) Тип структуры
6 (т.к. для всей структуры RasMol писал 19, посмотрела в pdb-файле, оказалось что по 6 в A и B цепях, а в цепи C - 7) 6 11 α+β (во вторичной структуре встречаются как α-спирали, так и β-листы, причем они не чередуются регулярно, как в α/β белках).

Выбор элементов вторичной структуры, опредение множества атомов остова для каждого элемента.

Для исследования выбраны следующие элементы вторичной структуры:

  • самая длинная α-спираль 69-88 (TYR69:A-ARG88:A) - myhelix
  • β-тяж 99-109 (THR99:A-GLY109:A) и антипараллельные по отношению к нему β-тяжи 36-44 (LYS36:A-SER44:A) и 120-130 (TRP120:A-GLY130:A) - mysheet
  • петля, в которой полипептидная цепь круто разворачивается на 180°, 130-137 (GLY130:A-GLU137:A) - myturn
Соседних параллельных β-тяжей в моем белке нет, поэтому исследоваться будут антипараллельные.

Используемый для получения изображения скрипт: all.spt.



А вот пара параллельных тяжей 366-370 и 397-400 для цепи А белка YEAST ARGINYL-TRNA SYNTHETASE (PDB ID: 1BS2), а также скрипт к ней: other.spt. Иллюстрации для параллельных тяжей приводить не буду, можно посмотреть их с помощью скрипта.

Создание картинок, иллюстирующих свойств выбранной альфа-спирали.

Используемый скрипт: helix.spt.

Определение наиболее часто встречающихся Н-связей.

Общий вид записи: Н(i,j), где i - номер остатка, которому принадлежит атом кислорода, j- номер остатка, которому принадлежит атом азота. ...,H(6,10),H(7,11),...,H(10,14),...- некоторое множество встречающихся Н-связей. Видим, что j=i+4. Обобщенный паттерн: H(i,i+4).

Определение расположения атомов Cβ относительно спирали.

Расположив спираль (в данном случае ее часть, т.к. сама спираль довольно длинная) перепедикулярно плоскости экрана, видим, что Cβ атомы, под приблизительно равными по отношению к связям углами (близкими по значению к тетраэдрическим), поочередно обращены на внешнюю сторону спирали.

Создание иллюстраций свойств выбранного бета-тяжа.

Используемый скрипт: sheet.spt.

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

Определение наиболее часто встречающихся Н-связей.

Составим таблицу для водородных связей для двух антипараллельных β-тяжей, указав номера остатков, начиная с N-концов β-тяжей и принимая за i - первый остаток первого тяжа, за k - второго (начиная с N-концов β-тяжей):

Тяж 36-44 36, i 38, i+2 40, i+4 42, i+6 44, i+8
Тяж 99-109 108, k+8 106, k+6 104, k+4 102, k+2 100, k

Если n - число остатков от начала первого β-тяжа до начала второго, то таблицу можно представить в виде:

Тяж 36-44 36, i 38, i+2 40, i+4 42, i+6 44, i+8
Тяж 99-109 108, i+n+8 106, i+n+4 104, i+n 102, i+n-4 100, i+n-8
Паттерн (i,i+n+8) (i,i+n+4) (i,i+n) (i,i+n-4) (i,i+n-8)

Теперь определение расположения атомов Cβ относительно "плоскости" β-листа.

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

Создание изображения реверсивного поворота.

Выбранный Cβ-поворот: Gly130Ala131Pro132Ile133Lys134Gly135Leu136Glu137 (130-137:A).

Н-связь (136,131) наиболее вероятно способна "скрепить" поворот.

Используемый скрипт: turn.spt.

Определение параметров выбранной спирали: шага, числа элементов на виток и радиуса.

Используемый скрипт: helixpar.spt.

Спираль правозакрученная (вкручивается вдаль по часовой стрелке), перед n ставим "+". Среднее число аминокислотных остатков на виток спирали равно 3,6. Параметр спирали +n=3,6.

Выше изображена пара атомов Cα с номерами остатков 76 и 80, находящихся друг под другом. Расстояние между ними равно 6.12 Å.

Шаг спирали p (минимальное расстояние между двумя эквивалентными точками) находим, измеряя расстояния между парами атомов Cα,находящихся друг под другом. В данном случае p=5.35 Å.

Смещение по оси на один остаток (параметр d спирали) можно вычислить по формуле d=p/n:

d=5,35/3,6=1,49 Å.

Для нахождения радиуса r спирали расположим ее часть перпендикулярно плоскости экрана, для удобства возьмем 72-85 остатки. 5.35 Å.

Рассмотрим 73-77 остатки. В таблице приведены названия атомов Cα данных остатков и их PDB-координаты по осям x и y.

HIS73:A.CA 56.712 -6.190
ALA74:A.CA 56.440 -9.376
LEU75:A.CA 58.159 -7.661
TYR76:A.CA 56.053 -4.496
HIS77:A.CA 52.807 -6.379

Находим расстояние между атомами в декартовой системе координат. Такое измерение позволит получить более точные значения, чем с помощью RasMol, т.к. он ведет измеререния в трехмерной системе координат, что неудобно в данном случае.

По формуле s=((x2-x1)2+(y2-y1)2)1/2 получаем:

s(73;75)=2,063 Å, s(74;76)=4,8953 Å, s(75;77)=5,503 Å.


Среднее значение 4,154 Å. Это диаметр спирали, следовательно, радиус r=2,077 Å, допуская же оценку "сверху", r<2,077 Å.

Описание бета-тяжа как спираль, оценка параметров такой спирали.

Используемый скрипт: sheethelix.spt.

Выше изображена пара атомов Cα с номерами остатков 106 и 108, находящихся друг под другом. Расстояние между ними равно 6.75 Å. Хотя, если совсем уж описывать бета-тяж как спираль, нужно было бы выбирать остатки, между которыми имеются три других.

Вроде бы, тяж вкручивается вдаль по часовой стрелке, т.е. - "правозарученный", перед n ставим "+". Среднее число аминокислотных остатков на такой виток спирали приблизительно равно 2. Параметр +n=2.

Шаг спирали p (минимальное расстояние между двумя эквивалентными точками) находим, измеряя расстояния между парами атомов Cα,находящихся друг под другом. В данном случае p=6.02 Å.

Смещение по оси на один остаток (параметр d) можно вычислить по формуле d=p/n:

d=6,02/2=3,01 Å.


© Eugenia Prokhorova 2011