Выпуск 284
Лаборатория Наномир
Когда реальность открывает тайны,
уходят в тень и
меркнут чудеса ...
Нанотех против Пикотех, что сабля против пулемёта...
Начинается с отметки 7 минут 20 секунд. Длительность ~2 минуты. Смотреть всем обязательно!
Проверяем альтернативные пикотехнологические модели аминокислот.
Victoria: Александр Юрьевич, а Вы могли бы написать скрипт для такого вида шаблона аминокислоты:
Кушелев: Первые четыре атома:
N = Hedra family:1 scalep:100 scaleq:100 scaler:100 radius:0.92 pos:[0,0,0] p:0.4 wirecolor:[0,0,250]; Converttomesh N
rotate N 54.73 [0,1,0]
rotate N 60 [1,0,0]
move N [1.33,0,0]
CA = Hedra family:1 scalep:100 scaleq:100 scaler:100 radius:0.91 pos:[0,0,0]
p:0.4 wirecolor:[200,250,250]; Converttomesh CA
rotate CA 54.73 [0,1,0]
move CA [2.79,0,0]
C = Hedra family:1 scalep:100 scaleq:100 scaler:100 radius:0.91 pos:[0,0,0] p:0.4 wirecolor:[200,250,250]; Converttomesh C
move C [3.28,1.43,0]
rotate C 15 [0,0,1]
C.pivot = [0,0,0]
rotate C -30 [1,0,0]
move C [0,-0.06,0]
O = copy C wirecolor:[250,0,0]
move O [0.33,0.5,0.06]
O.pivot = [3.28,1.43,-0.5]
scale O [0.99,0.99,0.99]
В принципе из этой модели уже можно пытаться собирать альфа-спираль и другие типы структур (пи-спираль, бета-спираль, 310-спираль)...
Victoria: А.Ю., в отличие от остальных атомов аминокислоты (N, CA, C) кислород (О) Вы не прописали в скрипте заново, а скопировали с имеющегося октаэдра для углерода С. Хотелось бы поправить его радиус, уменьшим до 0.88 и установить координаты его центра. Вы многократно использовали в скрипте команды move и rotate поэтому я не пойму, где конечные координаты центра красного октаэдра.
Кушелев: Координаты можно измерить в параллельных проекциях 3DS Max.
Victoria: Модель правой альфа спирали из последней модели шаблона для аминокислот:
Кушелев: Уважаемая Виктория!
К сожалению, в Вашей модели альфа-спирали на один виток приходится не 3.6 аминокислотных остатка, а меньше трёх, т.е. ~2/5
Попытайтесь ещё раз сделать модель альфа-спирали.
Виктория: Не может быть. А Вы не ошиблись? Вертикальные водородные связи образуются в этой модели как раз между кислородом первого а.о. и азотом четвёртого. Это значит, что всё-таки 3,5 - 4 а.о. на виток.
Кушелев: Водородные связи - дело десятое. Важно, на сколько аминокислотных остатков поворачивается спираль, совмещаясь со своим прежним положением. В Вашей модели явно меньше 2.5 аминокислотных остатков на виток, а экспериментальное значение 3.6 ... 3.7
Вы однозначно проделали большую работу, провели ряд модельных экспериментов. Фактически Вы получили доказательство того, что от модели альфа-спирали из современного учебника неизбежен переход к кольцегранной пикотехнологической модели. Строение аминокислот и аминокислотных остатков имеют особенности, характерные для полуацеталей, в которых углы между направлениями на центры атомов не просто отличаются от тетраэдрических, а могут плавно меняться.
Альфа-спираль из многогранников никак не получается, т.е. поворот на каждом аминокислотном остатке получается существенно больше, чем на угол 360/3.6=100 градусов.
В частности, в этой модели поворот на каждом аминокислотном остатке происходит на угол, превышающий 150 градусов:
Можно сделать, как в учебнике, т.е. уменьшить этот угол, но при этом грани многогранников уже не сойдутся. Более того, из учебника органической химии известно, что углы в альфа-спирали не тетраэдрические:
Модель группы R-COOH из учебника.
Все видят, что атомы лежат в одной плоскости?
Модель альфа-спирали из учебника. Угол, на который поворачивает спираль на каждом аминокислотном остатке ~100 градусов.На этой модели альфа-спирали (из учебника) видно, что атомы группы СО, атом Сальфа и атом азота следующего аминокислотного остатка лежат в одной плоскости. Единственная неточность - радикалы перепутаны с атомами водорода. Отсюда и разговор об L-форме аминокислот. Достаточно переставить радикалы, и форма L1 перейдёт в D2
4 типа спиралей (альфа-, 310-, бета- и пи-) можно посмотреть в 205-ом выпуске рассылки "Новости лаборатории Наномир": http://nanoworld88.narod.ru/data/205.htm
В частности, пикотехнологическая
модель альфа-спирали:http://video.yandex.ru/users/nanoworld/view/2352/#hq
Подробнее: http://nanoworld88.narod.ru/data/220.htm
Victoria: Примерно такая альфа-спираль:
Кушелев: На вид очень симпатичная спираль...
Я вижу в Вашей модели разные композиционные углы. Пикантная идея, которая заключается в том, что по коду альфа-спирали, например, ggc ggc ggc ggc ggc ggc ggc ggc
ggc ggc ggc ggc, рибосома ставит аминокислоты в растущую белковую цепь под разными композиционными углами, меняя их таким образом, чтобы в среднем получалось 3.6 ... 3.7 аминокислотных остатков на виток альфа-спирали. Я правильно понял Вашу идею? Если так, то давайте её проверим.
Один из способов проверки - измерить длину
водородных связей. Известно, что в реальной альфа-спирали все водородные связи, соединяющие атомы азота с атомами кислорода, практически одинаковы. Об этом говорит спектр, с помощью которого длину водородной связи можно определить с точностью не хуже 0.1%. Это значит, что расстояние между центрами атомов кислорода и азота, которые соединены водородной связью, не могут отличаться больше, чем на 0.1%. В Вашей модели видно, что это расстояние меняется в разы, а в реальной альфа-спирали все композицзионные углы
одинаковы. Так что модель нужно подправить.
Victoria: Вообще то "композиционные углы" планировались одинаковыми . Но если в ручном режиме я присоединила аминокислоты не совсем точно, то давайте напишем скрипт, по которому аминокислоты будут единообразно присоединяться друг к другу, водородные связи будут одинаковыми и мы сможем посчитать количество а.о. на виток в такой модели альфа-спирали.
Вы могли бы это сделать? У Вас со скриптами лучше получается, Вы практически на ты с ними.
Кушелев: Я уже пробовал. Одинаковые тетраэдрические углы не дают параметров альфа-спирали ни при каких комбинациях. Поэтому и удивился, когда Вы показали последнюю модель.
Кушелев: Мы с Дмитрием Кожевниковым ещё в 1992-ом году поняли, что альфа-спираль из тетраэдрических углов не составить. Я ему сразу сказал, что делать модель аминокислоты - практически безнадёжное дело. Но он парень дотошный и настойчивый. Он достал толстенный учебник по органической химии и нашёл там структурную схему полуацетали, где связи между атомами имеют ту же природу, что и в аминокислоте. Я ему говорил: "Брось ты это безнадёжное дело!", а он мне и отвечает: "Раз модель бензола и модели нуклеотидов получились, хотя ты тоже говорил "безнадёжное дело", то может быть и модель аминокислоты получится". Ну и короче, добился того, что я тоже начал вместе с ним делать модель аминокислоты. С очередной попытки мы получили эту модель полуацетали:
или ацетали. Я уже точно не помню название...
Потом была сделана модель молекулы мочевины
Далее к этой модели была добавлена группа азота:
Так и появилась первая пикотехнологическая модель аминокислоты...
Мы выбрали L-форму, согласно учебнику. Позднее я понял, что радикал нужно переставить. При этом L1-форма превратилась в D2-форму, которую я заменил на зеркально-симметричную, т.е. L2-форму. Из неё получается зеркальная спираль. Теперь это недоразумение устранено, и правая альфа-спираль строится из D2-форм аминокислотных остатков.
Виктория: Мой вариант моделей аминокислот:
Кушелев: Я правильно понял, что осталось одно несущественное отличие, радикал. Переставим?
Кстати, на Вашей иллюстрации упрощённая модель не соответствует кольцегранной. Причём по угловым величинам сильно не соответствует. И по линейным величинам тоже...
А самое большое несоответствие Вашей и моей модели конечно же в положении радикала. Я бы сказал "радикальное несоответствие". Но в случае глицина этого радикального несоответствия нет, поэтому давайте разберёмся, есть ли различия между нашими моделями глицина ?
Модель B и C с Вашей иллюстрации (если не считать радикала) в точности совпадают с моей, т.к. они из неё и сделаны путём перестановки радикала
А раз совпадают кольцегранные, значит должны совпасть и упрощённые, т.е. многогранные.
А раз совпадают упрощённые, значит совпадает и скрипт.
Единственное отличие, которое осталось - это другое расположение радикала. Тут можно проверить на уровне белков, где существенно положение радикалов.
Попробуем проверить на пептидах с дисульфидными мостиками...
фрагмент препроокситоцина, состоящий из 9 аминокислотных остатков. После сборки рибосомой этот фрагмент модифицируется и вырезается из белка-предшественника.
Есть похожий белок, вазопрессин.
Композиционный код фрагмента пропревазопрессина: 1,1,1,1,1,1,4,1,1
dne-файл:
60FT
FT source 1..545
FT CDS <108..134
XX
SQ Sequence 250 BP; 36 A; 105 C; 78 G; 31 T; 0 other;
acagagccac
caagcagtgc tgcatacggg gtccacctgt gtgcaccagg atgcctgaca
ccatgctgcc cgcctgcttc ctcggcctac tggccttctc ctccgcgtgc tacttccaga
actgcccgag gggcggcaag agggccatgt ccgacctgga gctgagacag tgcctcccct
gcggccccgg gggcaaaggc cgctgcttcg ggcccagcat ctgctgcgcg gacgagctgg
gctgcttcgt gggcacggct gaggcgctgc gctgccagga ggagaactac ctgccgtcgc
cctgccagtc cggccagaag gcgtgcggga gcgggggccg ctgcgccgcc ttcggcgttt
gctgcaacga cgagagctgc gtgaccgagc ccgagtgccg cgagggcttt caccgccgcg
cccgcgccag cgaccggagc aacgccacgc agctggacgg gccggccggg gccttgctgc
tgcggctggt gcagctggcc ggggcgcccg agcccttcga gcccgcccag cccgacgcct
actgagcccc gcgctcgccc caccggcgcg ctcttcgcgc ccgcccctgc agcacggaca
ataaacctcc gccaatgcaa a
//
Вазопрессин / vasopressin получается из этого фрагмента путём модификации, т.е. изменения композиционных углов после сборки белка-предшественника рибосомой.
Фрагмент А-цепи инсулина, который должен быть замкнут через дисульфидный мостик.
Однако, для этого нужно изменить один или несколько композиционных углов после сборки А-цепи рибосомой. Интересно разобраться, какие углы нужно изменить, т.е. какие аминокислотные остатки повернуть, чтобы цикл замкнулся через дисульфидный мостик?
Я припоминаю, что приблизительно такая же последовательность композиционных углов у меня получалась в пластмассовой модели цикла окситоцина. Именно модель окситоцина помогла мне понять, что кроме альфа-спирали, 310-спирали и бета-спирали есть ещё одна спираль, которую я назвал кси-спиралью, но оказывается у неё уже есть общепринятое название, пи-спираль.
Сейчас я подберу композиционный код, чтобы замкнуть цикл из 6 аминокислотных остатков. При этом в случае изменения позиции радикалов цикл должен разомкнуться. Это и будет тот эксперимент, который позволяет выбрать один из двух положений радикалов...
Начальная последовательность инсулина в этой зоне была: 333123
После модификации, когда композиционный код стал соответствовать: 333131, дисульфидный мостик (2 атома вверху, показанные жёлтым цветом) замкнулся.
Теперь интересно проверить, возможно ли такое, если переставить радикалы по варианту Виктории Соколик?
Уважаемая Виктория! Я выкладываю файл для 3DS Max, в котором радикалы цистеинов отделены от остальной части модели.
http://nanoworld88.narod.ru/3DS/333131_d2.max
Вы можете переставить их по своему усмотрению. Если получится альтернативная модель, в которой радикалы цистеинов дотянутся друг до друга, тогда Ваш вариант модели не будет противоречить результату модельного эксперимента.
Кушелев: А как быть с кодом gga ?
Код этого фрагмента белка: tgt tgt ctt agc gga tgc
Так что Ваш вариант без модификации этого участка белка после рибосомальной сборки тоже не получается...
Кушелев: Так я об этом сразу написал. И в Вашей точно не получится, т.к. в нуклеотидной последовательности есть код бета-спирали gga.
Но мы с Вами знаем, что после сборки на рибосоме белок-предшественник (преинсулин) модифицируется. Углы меняются. Пикотехнологическая модель показывает, что сразу цикл окситоцина, вазопрессина, инсулина... сделать не получается. Не влезают в цикл АСС-концы тРНК. Поэтому и приходится делать другую структуру (участок альфа-спирали, как в случае окситоцина и вазопрессина или участок пи-спирали, как в случае инсулина), из которой потом путём поворота двух (в случае инсулина) аминокислотных остатков получается цикл из 6 остатков, замкнутых через дисульфидный мостик.
Модель окситоцинового цикла я сделал из пластмассовых колец ещё в 1992-ом году:
структура окситоцина из видеосюжета "Формы и механизмы молекул":http://video.mail.ru/bk/yaroslav2008/22/39.html
Кушелев: Уважаемая Виктория!
Я Вам открою ещё один секрет, который позволяет сделать однозначный выбор между положениями радикалов аминокислот.
Напомню, что метионин гнёт альфа-спираль. Как он это делает? Он вставляет между собственным атомом кислорода группы СО и атомом азота аминокислотного остатка следующего витка дополнительный атом углерода, точнее группу CH3. Если радикал переставить, то он не дотянется...
Идёт создание новой версии программы PikoTech,
которая уже начинает строить упрощённые 3D-структуры белков по композиционному коду...
Структуру белка ещё "плющит", но это - временные проблемы
Ну вот, уже не плющит...
***
[0:26:18] Кушелев Александр Юрьевич: Точка стоит где нужно
[0:26:25]
V: нашёл ошибку )
[0:27:53] Кушелев Александр Юрьевич: Ура, получилось! Первый релиз работает
[0:28:17] V: вопрос на сколько точно он рабтает
)
[0:28:23] Кушелев Александр Юрьевич: Он уже сейчас подходит для анализа структур крупных белков. Он показывает один атом каждого аминокислотного остатка
[0:28:47] Кушелев Александр Юрьевич: Этого достаточно, чтобы видеть контуры крупных белков. Следующий релиз - это прорисовка вместо одного атома - 4 общих для всех аминокислот атомов. Координаты их известны...
[0:29:39] Кушелев Александр Юрьевич: Так что поздравляю с очередной победой!
[0:29:53] Кушелев Александр Юрьевич: Ты можешь скрипт-то прислать?
[0:30:35] *** V sent Colmodel.m ***
[0:30:49] V: это скрипт не для 3дМаши!!!!
[0:30:59] Кушелев Александр Юрьевич: А для кого?
[0:31:10] V: это скрипт для MATLAB
[0:31:17] Кушелев Александр Юрьевич: Тоже интересно. А для 3DS Max сможешь сделать?
[0:32:17] V: Буду разбираться ещё ...
[0:34:48] Кушелев Александр Юрьевич: Кстати, скрипт матлаба тоже неплохо смотрится...
[0:35:51] Кушелев Александр Юрьевич: Возможно, что имеет смысл параллельно работать и со скриптом матлаба, и со скриптом Макса
Создан скрипт для MatLab, который строит пикотехнологическую модель стебля белковой молекулы по композиционному коду.
% p = #()
% angx = #(-24,180,60,-45,-45); angy = #(-10.89167,83,0,15,15); angz = #(92,120,60,110,110)
% dx = #(2,1.5,1.6,2,2); dy = #(0.6,1,0.8,1,0.6); dz = #(-0.45,0,-0.6,-0.15,-0.45)
% --collagen
% kk = #(3,3,1,3,3,1,1,1,1,2,1,1,1,3,3,1,1,3,3,3,3,2,3,1,3,3,1,3,1,3,3,1,1,3,
% 1,1,3,3,1,3,1,1,3,1,3,3,1,4,1,3,1,1,1,3,3,3,1,3,1,1,3,1,1,1,1,1,3,3,
% 3,1,3,3,2,3,1,3,1,1,1,1,3,3,1,1,1,1,1,3,1,3,3,3,3,3,3,3,3,1,3,3,3,1,
% 1,1,3,1,1,1,3,3,1,2,1,3,1,1,1,1,1,3,1,1,1,3,1,1,1,2,1,3,3,1,1,3,3,3,
% 3,3,3,1,1,3,3,1,3,1,1,3,2,1,3,3,3,3,1,1,3,1,1,1,3,1,3,1,3,4,4,1,1,1,
% 1,1,1,1,1,1,1,3,1,1,3,3,1,1,1,3,3,1,1,1,3,1,1,1,4,3,3,1,1,1,1,1,1,1,
% 1,1,1,1,1,4,3,3,4,1,3,1,1,4,1,1,1,1,1,1,1,2,1,1,1,1,1,4,1,1,1,1,1,4,
% 1,4,1,1,1,4,1,1,1,1,3,1,1,1,1,1,1,1,1,3,1,1,1,1,1,4,1,1,1,1,4,1,4,1,
% 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,4,3,1,1,1,1,1,1,1,1,1,4,3,1,
% 4,1,1,1,1,4,3,1,1,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,4,1,1,
% 1,1,1,1,3,1,1,1,3,1,1,1,1,4,3,3,1,1,1,3,4,1,1,4,1,4,1,1,1,1,1,1,1,1,
% 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,4,1,1,3,4,4,1,1,2,1,1,1,1,1,1,1,
% 1,3,2,1,1,1,1,1,3,4,1,1,1,3,1,1,4,3,1,3,2,2,1,1,3,2,4,3,1,1,1,4,4,1,
% 1,4,1,4,1,1,3,1,1,3,3,1,1,3,1,3,4,1,1,3,3,1,3,1,1,1,1,1,1,1,1,1,3,1)
% peptide = box length:1 width:1 height:1 position:[0,0,-.5] wirecolor:[250,250,250]; Converttomesh peptide
% for k = 1 to 473 do(
% element = Hedra family:1 scalep:100 scaleq:100 scaler:100 radius:1 pos:[0,0,0] p:0.4; Converttomesh element
% attach peptide element; peptide.pivot = [0,0,0]; move peptide [dx[kk[k]],dy[kk[k]],dz[kk[k]]]
% rotate peptide angx[kk[k]] [1,0,0]; rotate peptide angy[kk[k]] [0,1,0]; rotate peptide angz[kk[k]] [0,0,1])
% rotate peptide 30 [0,1,0]; move peptide [-21.5,5,0]
% TAtom p:[ 27.85, -19, -26.5] ch1:"N" ch2:"
" ch3:" ",
% TAtom p:[ 9.7, -14.4, 0 ] ch1:"C" ch2:"A" ch3:" " con:#(3,5),
% TAtom p:[ -23.5, -21.0, 0 ] ch1:"C" ch2:" " ch3:" ",
% TAtom p:[ -44.5, -11.6, 0 ] ch1:"O" ch2:" " ch3:" " con:#()
function ColModel()
%global Cnt Peptide;
% задаём углы и смещения
angx = [-24 180 60 -45 -45];
angy = [-10.89167 83 0 15 15];
angz = [ 92 120 60 110 110];
dx = [2 1.5 1.6 2 2];
dy = [0.6 1 0.8 1 0.6];
dz = [-0.45 0 -0.6 -0.15 -0.45];
% --collagen
kk = [3 3 1 3 3 1 1 1 1 2 1 1 1 3 3 1 1 3 3 3 3 2 3 1 3 3 1 3 1 3 3 1 1 3 ...
1 1 3 3 1 3 1 1 3 1 3 3 1 4 1 3 1 1 1 3 3 3 1 3 1 1 3 1 1 1 1 1 3 3 ...
3 1 3 3 2 3 1 3 1 1 1 1 3 3 1 1 1 1 1 3 1 3 3 3 3 3 3 3 3 1 3 3 3 1 ...
1 1 3 1 1 1 3 3 1 2 1 3 1 1 1 1 1 3 1 1 1 3 1 1 1 2 1 3 3 1 1 3 3 3 ...
3 3 3 1 1 3 3 1 3 1 1 3 2 1 3 3 3 3 1 1 3 1 1 1 3 1 3 1 3 4 4 1 1 1 ...
1 1 1 1 1 1 1 3 1 1 3 3 1 1 1 3 3 1 1 1 3 1 1 1 4 3 3 1 1 1 1 1 1 1 ...
1 1 1 1 1 4 3 3 4 1 3 1 1 4 1 1 1 1 1 1 1 2 1 1 1 1 1 4 1 1 1 1 1 4 ...
1 4 1 1 1 4 1 1 1 1 3 1 1 1 1 1 1 1 1 3 1 1 1 1 1 4 1 1 1 1 4 1 4 1 ...
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4 3 1 1 1 1
1 1 1 1 1 4 3 1 ...
4 1 1 1 1 4 3 1 1 3 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4 1 1 ...
1 1 1 1 3 1 1 1 3 1 1 1 1 4 3 3 1 1 1 3 4 1 1 4 1 4 1 1 1 1 1 1 1 1 ...
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4 1 1 3 4 4 1 1 2 1 1 1 1 1 1 1 ...
1 3 2 1 1 1 1 1 3 4 1 1 1 3 1 1 4 3 1 3 2 2 1 1 3 2 4 3 1 1 1 4 4 1 ...
1 4 1 4 1 1 3 1 1 3 3 1 1 3 1 3 4 1 1 3 3 1 3 1 1 1 1 1 1 1 1 1 3 1];
% начальная точка
%Cnt1 = [ 0, 0, 0];
Cnt1 = [ 27.85, -19, -26.5]; Cnt1=Cnt1*0.051;
%Cnt2 = [ 10, 0, 0 ];
Cnt2 = [ 9.7,
-14.4, 0 ];Cnt2=Cnt2*0.051;
Cnt3 = [-23.5, -21.0, 0 ];Cnt3=Cnt3*0.051;
Cnt4 = [-44.5,
-11.6, 0 ];Cnt4=Cnt4*0.051;
for i = 1:1:473
Peptide1(i,1:3) = Cnt1;
Peptide2(i,1:3) = Cnt2;
Peptide3(i,1:3) = Cnt3;
Peptide4(i,1:3) = Cnt4;
% сдвинем
delta = [dx(kk(i)), dy(kk(i)), dz(kk(i)) ];
for j=1:1:i
Peptide1(j,1:3) = Peptide1(j,1:3) + delta;
Peptide2(j,1:3) = Peptide2(j,1:3) + delta;
Peptide3(j,1:3) = Peptide3(j,1:3) + delta;
Peptide4(j,1:3) = Peptide4(j,1:3) + delta;
end
%повернём
q = angle2quat(angx(kk(i))*0.0175,angz(kk(i))*0.0175,angz(kk(i))*0.0175,'XYZ');
M = quat2dcm(q);
for j=1:1:i
Peptide1(j,1:3) = Peptide1(j,1:3)*M;
Peptide2(j,1:3) = Peptide2(j,1:3)*M;
Peptide3(j,1:3) = Peptide3(j,1:3)*M;
Peptide4(j,1:3) = Peptide4(j,1:3)*M;
end
end
plot3(Peptide1(1:473,1),Peptide1(1:473,2),Peptide1(1:473,3),'DisplayName','Peptide');hold on;
plot3(Peptide2(1:473,1),Peptide2(1:473,2),Peptide2(1:473,3),'DisplayName','Peptide');
plot3(Peptide3(1:473,1),Peptide3(1:473,2),Peptide3(1:473,3),'DisplayName','Peptide');
plot3(Peptide4(1:473,1),Peptide4(1:473,2),Peptide4(1:473,3),'DisplayName','Peptide');
%Pause();
fpdb=fopen('collagen.pdb','w');
fprintf(fpdb,'PFRMAT TS\n');
fprintf(fpdb,'TARGET\n');
fprintf(fpdb,'AUTHOR TEST\n');
fprintf(fpdb,'REMARK Nanoworld laboratory\n');
fprintf(fpdb,'METHOD Method description\n');
fprintf(fpdb,'MODEL 1\n');
fprintf(fpdb,'PARENT N/A\n');
%fwrite( fpdb,SH);
for i=1:1:473
str='ATOM
N GLY A 1 1.00 0.50 A N\r\n';
u=num2str((i-1)*4+1,'%d'); uu=size(u); str(12-uu(2):11)
= u;
u=num2str(Peptide1(i,1),'%6.3f'); uu=size(u); str(39-uu(2):38) = u;
u=num2str(Peptide1(i,2),'%6.3f');
uu=size(u); str(47-uu(2):46) = u;
u=num2str(Peptide1(i,3),'%6.3f'); uu=size(u); str(55-uu(2):54) = u;
fprintf(fpdb,str);
str='ATOM CA GLY A
1 1.00 0.50 A N\r\n';
u=num2str((i-1)*4+2,'%d'); uu=size(u); str(12-uu(2):11) = u;
u=num2str(Peptide2(i,1),'%6.3f'); uu=size(u); str(39-uu(2):38) = u;
u=num2str(Peptide2(i,2),'%6.3f'); uu=size(u); str(47-uu(2):46) = u;
u=num2str(Peptide2(i,3),'%6.3f'); uu=size(u); str(55-uu(2):54) = u;
fprintf(fpdb,str);
str='ATOM
C GLY A 1 1.00 0.50 A N\r\n';
u=num2str((i-1)*4+3,'%d'); uu=size(u); str(12-uu(2):11) = u;
u=num2str(Peptide3(i,1),'%6.3f'); uu=size(u); str(39-uu(2):38) = u;
u=num2str(Peptide3(i,2),'%6.3f');
uu=size(u); str(47-uu(2):46) = u;
u=num2str(Peptide3(i,3),'%6.3f'); uu=size(u); str(55-uu(2):54) = u;
fprintf(fpdb,str);
str='ATOM O GLY A 1 1.00 0.50 A N\r\n';
u=num2str((i-1)*4+4,'%d');
uu=size(u); str(12-uu(2):11) = u;
u=num2str(Peptide4(i,1),'%6.3f'); uu=size(u); str(39-uu(2):38) = u;
u=num2str(Peptide4(i,2),'%6.3f'); uu=size(u); str(47-uu(2):46) = u;
u=num2str(Peptide4(i,3),'%6.3f'); uu=size(u); str(55-uu(2):54) = u;
fprintf(fpdb,str);
%fwrite(fpdb,
str);
end
fclose(fpdb);
end
В данном варианте скрипта входными данными является композиционный код белка коллагена, который выводит программа от Prosolver. В результате получаем pdb-файл стебля белка, т.е. без радикалов:
Фрагмент pdb-файла, который выводит этот скрипт:
ATOM 2 CA GLY A 1 -54.250-156.124 84.413 1.00 0.50 A N
ATOM
3 C GLY A 1 -54.520-156.695 86.020 1.00 0.50 A N
ATOM 4 O GLY A 1 -54.192-156.588 87.142 1.00 0.50 A
N
ATOM 5 N GLY A 1 -53.964-159.070 87.129 1.00 0.50 A N
ATOM 6 CA GLY A 1 -53.345-157.584 86.745 1.00 0.50 A N
ATOM 7 C GLY A 1 -53.650-156.200 87.732 1.00
0.50 A N
ATOM 8 O GLY A 1 -54.243-155.191 87.807 1.00 0.50 A N
ATOM 9 N GLY A 1 -52.051-154.151 88.379 1.00 0.50 A N
ATOM 10 CA GLY A 1 -53.088-154.741
87.232 1.00 0.50 A N
ATOM 11 C GLY A 1 -54.800-154.616 87.417 1.00 0.50 A N
ATOM 12 O GLY A 1 -55.783-155.212 87.650 1.00 0.50 A N
ATOM 13 N GLY
A 1 -56.361-154.056 85.520 1.00 0.50 A N
ATOM 14 CA GLY A 1 -56.338-155.183 86.732 1.00 0.50 A N
ATOM 15 C GLY A 1 -56.475-154.723 88.390 1.00 0.50 A N
ATOM
16 O GLY A 1 -55.892-154.603 89.401 1.00 0.50 A N
ATOM 17 N GLY A 1 -57.742-155.925 90.421 1.00 0.50 A
N
ATOM 18 CA GLY A 1 -56.341-155.886 89.541 1.00 0.50 A N
ATOM 19 C GLY A 1 -55.124-154.699 89.841 1.00 0.50 A N
ATOM 20 O GLY A 1 -54.619-153.721
89.433 1.00 0.50 A N
ATOM 21 N GLY A 1 -52.454-154.905 89.780 1.00 0.50 A N
ATOM 22 CA GLY A 1 -53.807-154.659 88.860 1.00 0.50 A N
ATOM 23 C GLY
A 1 -54.574-153.115 88.783 1.00 0.50 A N
ATOM 24 O GLY A 1 -55.508-152.492 89.124 1.00 0.50 A N
ATOM 25 N GLY A 1 -55.513-152.138 86.658 1.00 0.50 A N
Технологии Гермеса...
Гермес: добрый день!Кушелев: pdb-файл содержит координаты атомов в молекуле белка. Нужен тем, кто использует структуру белка в своей работе. Фармацевты, генные инженеры, биотехнологи и все, кто связан с нанотехнологией. Пикотехнология - это развитие нанотехнологии, только в ~1000 раз точнее.
Получают pdb-файлы с помощью рентгеноструктурного анализа или с помощью нейтронной дифракции. Для этого несколько месяцев выращивают кристалл из миллиардов миллиардов одинаковых молекул исследуемого белка. После этого просвечивают белковый кристалл рентгеном или нейтронами. Получают дифракционные картины в разных ракурсах и анализируют на суперкомпьютерах.
Программа "Пикотех" позволяет смоделировать работу рибосомы и узнать структуру белка в ~1000 раз точнее и в миллион раз быстрее. Кроме того рентген и нейтроны пригодны лишь для 3% белков, которые кристаллизуются. Пикотех позволяет определать структуру белка независимо от способности кристаллизоваться.
Ежегодно продаётся более 1000 pdb-файлов по средней цене ~10 000 евро. Все покупатели - наши потенциальные клиенты.
Пригоден ли конкурс CASP для пикотехнологов?
Victoria: На конкурсе вывесили новую задачу:
>R0011 CD27L , , 81 residues
HTIVYDGEVDKISATVVGWGYNDGKILICDIKDYVPGQTQNLYVVGGGACEKISSITKEK
FIMIKGNDRFDTLYKALDFIN
Это кусок рецептора к цитокину CD27. Вот только организм они опять не указали, поэтому нуклеотидную последовательность найти будет почти невозможно.
Кушелев: Нашёл другой цитокин:
http://www.genome.jp/dbget-bin/www_bget?rno:500318
http://www.genscript.com/product_001/ge … 335.1.html
60FT
FT source 1..755
FT CDS <1..755
XX
SQ Sequence 250 BP; 36 A; 105 C; 78 G; 31 T; 0 other;
atggcatggccacccctctactggctctgcatgctggggaccttggtagggctcttagct
accccagccccaaacaactgtccagacagacactactggattggggcaggactctgctgc
cagatgtgtgggccagggacattccttgtgaagcactgtgaccaggacagagccgctgct
cagtgtgatccgtgtattccaggcacgtccttctctccagactaccacacccggccccac
tgcgagagctgcaggcactgtaactctggttttctcatccgcaactgcacggtcactgcc
aatgctgagtgcacctgttccaagggctggcagtgcagggaccaggaatgtacggagtgt
gaccctcctctaaaccctgcactgactagccagccctctgaggcccccagcccccagctg
ccaccacccacccacttaccttatgccacagagaagccatcctggcccccacaaaggcag
cttcccgactcaactgtctatagccgcctgccatcccaaagacctctgtgcagctccgac
tgtatccggatctttgtgaccttctccagcatgcttcttgtcttcgtcctgggtgggatc
ttgttcttccatcaaagaagaaatcatgggccaaatgaagatagccaggcagtacctgaa
gagctttgtccttacagctgccctagggaggaggagggcagtgtcatccctatccaggaa
gactaccggaaacctgagcctgcttcctacccttgatgcggtgctggggggtcggag
//
Композиционный код:
1,2,1,2,1,1,1,1,1,1,1,4,4,1,1,2,4,1,3,3,1,2,1,2,1,1,3,2,1,3,1,1,1,3,4,2,2,1,1,1,1,1,3,4,2,4,2,1,3,4,
1,1,3,1,1,1,3,1,3,3,1,3,3,4,3,3,2,1,4,1,1,3,2,1,1,1,1,4,1,1,1,1,1,1,1,1,3,1,3,3,3,1,1,1,1,1,4,1,3,1,
3,3,1,1,1,3,1,1,1,1,1,1,1,1,1,3,3,4,1,3,1,3,3,2,1,3,2,4,3,1,1,1,3,1,1,1,1,1,1,4,2,2,1,1,1,3,3,3,1,2,
1,1,2,1,1,1,2,3,1,1,3,1,1,2,3,1,3,1,1,4,2,1,3,3,3,4,1,1,1,1,3,1,4,1,3,4,1,1,1,1,1,3,3,1,1,1,4,3,4,1,
1,1,1,3,3,3,3,3,3,4,2,3,3,3,1,1,2,2,3,3,1,3,3,3,1,1,1,3,1,1,1,1,1,3,1,1,3,1,1,3,1,1,4,3,3,1,3,3,1,1,
3
Cytocine (PyMOL, mesh)
Если нужен pdb-файл, моглу выложить...
Готовим встречу с Даниилом Шостиным...
Skype:
2012-02-20
[14:39:07] Alexey Komanyov: приветствую)
[14:41:31] Кушелев Александр Юрьевич: Рад видеть!
[14:41:37] Alexey Komanyov: взаимно)
[14:41:53]
Кушелев Александр Юрьевич: Ну что, организуем встречу с Даниилом Шостиным?
[14:42:05] Alexey Komanyov: в каком формате?
[14:43:33] Кушелев
Александр Юрьевич: Интересно было бы для начала передать в Екатеринбург учебные пособия с пикотехнологическими моделями и ... объяснить, что он фактически создал разновидность пикотехнологического конструктора, который является отличным дополнением к учебным пособиям по пикотехнологии, которые уже утверждены министерством образования и науки России
[14:44:41] Кушелев Александр Юрьевич: А речь можно организовать через прямую интернет-трансляцию, которые организует Ярослав Старухин.
[14:45:02] Кушелев Александр Юрьевич: Вообще-то я хотел лично приехать в Екатернибург, но эту поездку нужно слегка подготовить
[14:45:25] Кушелев Александр Юрьевич: Первый контакт с ***** интересно было бы организовать заранее, чтобы народ уже был готов к приезду делегации из Москвы
[14:45:44] Alexey Komanyov:
если вы сами сюда приедете, я готов составить вам компанию и оказать определенное организационное содействие, в т.ч. освещение в СМИ и на РП
[14:46:04] Кушелев Александр Юрьевич: Это я уже понял ещё при первой нашей встрече
[14:46:33] Кушелев Александр Юрьевич: Два мероприятия вместо одного ...
[14:46:39] Alexey Komanyov: два?
[14:47:09] Кушелев Александр Юрьевич: Ну да. Сначала присылаем учебные плакаты и организумем прямую интернет-трансляцию, а второе мероприятие - едем из Москвы в Екатеринбург
[14:48:34] Кушелев Александр Юрьевич: Кстати, видели материал из Храма Солнца?
[14:48:50] Alexey Komanyov: ... могу обещать промоушен этого события со всех сторон
[14:50:19] Кушелев Александр Юрьевич: http://nanoworld88.narod.ru/data/251.htm
[14:52:27] Кушелев Александр Юрьевич: Понятно. Тогда нам нужно связаться с ... Подробнее: http://nanoworld88.narod.ru/data/274.htm
[14:53:09] Кушелев Александр Юрьевич: Речь
идёт о новой (пико-) технологии, но с точки зрения обывателя научные открытия - тоже аномальные явления
[14:53:18] Кушелев Александр Юрьевич: А тут конкретно привязка к Екатеринбургу
[14:53:33] Кушелев Александр Юрьевич: Ведь Даниил Шостин - ученик 9-ой школы Екатеринбурга ...
[14:54:04] Кушелев Александр Юрьевич: И ему удалось создать и запатентовать новый конструктор, который, как выяснилось, является научным открытием...
[14:54:41] Кушелев Александр Юрьевич: Слабо взять в руки конструктор Даниила Шостина собрать из него что-нибудь?
[14:56:15] Кушелев Александр Юрьевич: Главное - ...
[14:56:34] Кушелев Александр Юрьевич: Ну и очень важно - организовать прямую интернет-трансляцию
[14:56:40] Кушелев Александр Юрьевич: Это
тоже могут сделать ...
[14:57:36] Кушелев Александр Юрьевич: ... на форуме паранормал можно начать обсуждение ...
[14:57:46] Alexey Komanyov:
я вам могу предоставить ...
[14:58:10] Кушелев Александр Юрьевич: ОК.
[14:59:44] Alexey Komanyov: вы готовы им красочно все объяснить?
[15:00:31] Кушелев Александр Юрьевич: Я могу красночно объяснить
[15:02:54] Alexey Komanyov: Тогда я от вас жду следующее: 1) пресс-релиз
грядущего события , его смысл и красоту (для СМИ) 2) письмо для ...
[15:32:40] Кушелев Александр Юрьевич: ОК!
Представим себе уменьшенную копию Стоунхэнджа, в которой все линейные размеры уменьшены в 10 раз. При этом объём каждого камня уменьшится в 10^3=1000 раз. Самый крупный камень в модели будет весить не 55 тонн, а 55 кг. Это - приблизительно средний вес студентки
Таким образом, мы получили игрушечный Стоунхэндж, который меньше настоящего в 10 раз по "росту" и в 1000 раз по весу.
А теперь попробуем сделать наборот, т.е. не уменьшим, а увеличим массу камней Стоунхэнджа в 1000 раз. Масса самого крупного камня (55 тонн) станет в 1000 раз больше, т.е. 55 000 тонн. Теперь вспомним, что масса одного из камней Храма Солнца в Ласпи (под Севастополем) превышает 70 000 тонн. Получается, что по сравнению с Храмом Солнца Стоунхэндж - игрушка
"Летающие тарелки" бороздят ... мониторы
2012-02-23Китайская клавиатура из ... Заповедника под голубым Солнцем
Китайская клавиатура подходит к монитору из фантастического рассказа "Заповедник под голубым Солнцем"
Цитата: "Несколько "человек" сидели за круглыми экранами..."
Это - не Пандора! Это - Испания...
Фантастический лес. Кантабрии, Испания.
Лавандовые поля, Прованс, Франции
Кристаллические мыльные пузыри...
При температуре от -6 до -9 пузыри замерзают примерно в течении двух минут. Если температура ниже, то они замерзают быстрее и бьются как стекло.Снежинка соприкоснется
с пузырем и от нее во все стороны побегут бусинки льда. Причем тип кристаллизации будет такой же, как и у снежинки, то есть можно выбрать его заранее, отлавливая снежинки и рассматривая их на варежке...
Ещё можно заморозить пузырь, если аккуратно положить его на снег.
Из-за своей маленькой толщины пленка пузыря при этом остаётся пластичной, так что если уронить бывший мыльный пузырь на землю, то он не разобьется на множество ледяных осколков. Вместо этого на пузыре образуются вмятины, а если уронить шарик посильнее, он распадется на кусочки, которые тут же скрутятся в трубочки
Приглашение к сотрудничеству
Лаборатория Наномир готова к любому взаимовыгодному сотрудничеству. У нас есть сторонники как явные, которые помогают морально и материально, есть очень много пассивных наблюдателей, есть и ярые противники, которые используют любые методы и средства (аморальные и просто преступные), чтобы уничтожить работу лаборатории и дискредитировать ее.
В одиночку
внедрить технологии, выводящие цивилизацию на новый уровень, невозможно. Благодаря поддержке множества заинтересованных людей проделана огромная работа. Ознакомиться с её результатами можно изучив материал рассылки "Новости лаборатории Наномир". Люди науки могут изучить научные труды.
Вклад каждого не останется незамеченным в случае успеха в реализации научных проектов. Результаты совместной деятельности принадлежат участникам проекта пропорционально коэффициентам творческого и финансового участия.
Основные разработки:
Микроволновый источник энергии
В 2011 году были куплены рубиновые шарики для эксперимента на сумму ~1000 долл. В результате было сделано научное открытие, проверена защита диэлектрических резонаторов от перенапряжения. В этом году, вероятно, можно будет создать микроволновую энергетику, т.к. удалось найти сырьё (рубин #8), из которого сделаны рубиновые шарики для эксперимента в Дубне.
Параллельно началась серия экспериментов на ускорителе "B" в диапазоне 20 ... 140 ГГц.
"Эликсир вечной молодости"
28 сентября 2011г. в институте геронтологии (г.Киев) начался эксперимент по созданию "эликсира вечной молодости". Благодаря первому взносу (в размере 500 долларов) Золдракса и поддержке других соинвесторов. Продолжаются переговоры с потенциальными инвесторами по поводу финансирования этого проекта.
Программа исследований Презентация
Пикотехнология
Созданы первые версии пикотехнологии (выпуски рассылки 212 270 271 272 ), с помощью которой Александр Кушелев и Виктория Соколик сделали более10 научных открытий.
Виктория Соколик: Уважаемые коллеги, Вашему вниманию предоставляется услуга -- моделирование 2D и 3D структуры любого белка без ограничений в его размере и степени изученности с помощью программного обеспечения, базирующемся на принципиально новом подходе декодирования нуклеотидной последовательности, детерминирующей данный белок.
Всё, что необходимо от заказчика, это нуклеотидная последовательность мРНК интересующего его белка (или код этой нуклеотидной последовательности в EMBL, или хотя бы код самого белка в PDB).
В течение 1-3 суток мы готовы предоставить Вам схему вторичной структуры заказанного белка (2D), модель его пространственной структуры (3D) в виртуальном пространстве, а также файл .pdb с координатами каждого атома белка.
Файл .pdb может быть использован по аналогии с файлами закристаллизованных белков из PDB банка для дальнейшего конформационного анализа белка методами молекулярной динамики с учётом физико-химической специфики микроокружения белка или его взаимодействия с лигандами.
Таким образом, Вы сможете максимально быстро удобным для Вас способом (по электронной почте, на сайте либо на электронном носителе) получить информацию о структурном шаблоне Вашего белка.
Первые 10 заказов -- бесплатно . Цена следующих заказов -- договорная.
Сотрудничество может быть различным:
- участие в научных дискуссиях на форуме (конструктивное)
- совместное создание коммерческого продукта
- поиск инвесторов
- выступить менеджером по продаже готовых коммерческих продуктов
- конструктивные предложения по продвижению идей лаборатории Наномир
- содействие в проведении экспериментов
и т.п.
- написание совместных научных статей и т.п.
- материальный вклад (денежный или обеспечение оборудованием и материалами)
Пожалуйста, сообщайте о своем вкладе, чтобы мы зачли Вас как партнера лаборатории Наномир.
+7-926-5101703, +7-903-2003424, +7-916-8265031, Skype: Kushelev2009, mail: kushelev2011@yandex.ru
веб-мани: WM-кошелек R426964799301