Период спирали в системе (5) длина poly(A)-цепи - 300 пар оснований. То есть параметры периода спирали в начальных условиях и в системе (5) различны. Таким образом смоделирован перенос ДНК из кристаллического состояния в увлажненное.
Граничные условия следующие (назовем их “квазициклическими”):
Особенностью данной модели является то, что при переходе из состояния с периодом в 10 пар в состояние с периодом в 10, 5 пар почти вся цепь оказывается денатурированной (“расплавленной”). Приведенные ниже результаты описывают процесс ренатурации такой цепи с возникновением дислокаций.
В этих экспериментах варьировались параметры: 1) диссипация 2) отношение параметров упругости 3) угол обрыва водородных связей .
На рис. 3 и 4 представлены результаты численного интегрирования системы (5). Показана не сама функция , а разница , поскольку область изменения функции (приблизительно от до ) велика по сравнению с характерными изменениями в системе (приблизительно от 0 до 9). Горизонтальная часть графиков соответствует нераспаренному участку цепи с периодом спирали . Наклонная часть графиков на рис. 3(a), 4(а) соответствует дислокации.
Можно сделать следующие выводы:
1) Способность к образованию дислокации в этой модели сильно зависит от . При дислокация возникла во всех рассмот-ренных случаях.
2) Способность к образованию дислокации также сильно зависит от параметра. Во всех случаях, когда параметр велик (
на рис. 1.а, 2.а ), дислокация возникла. В пользу этого утверждения также свидетельствует сравнение рис. 3(а) и 4(г).
Как показывают дополнительные расчеты, влияниена эффект проявляется в меньшей степени. Дислокация образуется или не образуется вне зависимости от значения ( или ). При больших значениях дислокация образуется медленнее, чем при меньших.
3) На рис. 3(а), 4(в,г) видно, что дислокация имеет кинкообразную форму.
Ширина дислокации зависит от параметров (чем больше , тем меньше ширина дислокации) и (чем больше , тем меньше ширина дислокации).
Развивая дальше модели солитонных возбуждений в ДНК (совместно с М.Ю.Масловым и др.) мы использовали условия, при которых цепочки ДНК моделируются набором ровибронных осцилляторов, подвешенных на невесомом нерастяжимом стержне; для простоты спирализация цепи не учитывается, а ровибронные степени свободы одной из цепочек считаются “замороженными”.
В этом случае гамильтониан для “активной” цепочки записывается в следующем виде:
H=H0+H1+H2
(1)
где: - число пар оснований в цепи; - гамильтониан, описывающий собственные осцилляции мономеров ( - углы вращения нуклеотидов в цепочке, - момент инерции оснований); - гамильтониан , характеризующий нелинейно-периодическую связь между осцилляторами (- константа упругости цепочки, ), - гамильтониан,
(а)
(б)
а)x0=200 б)x0=250
Рис.3
в) г)
в) x0=300 г) x0=350
Рис. 4
описывающий нелинейную связь между “активной” и “замороженной” () цепочками ДНК (- константа упругости водородных связей между комплементарными основаниями, коэффициенты в уравнении (1) определяются в соответствии с правилом: в случае АТ и ТА пар, в случае ГЦ и ЦГ пар; - параметр, полученный ранее (см. выше) и определяемый на основе модели синус-Гордона).
При малых гамильтониан, что совпадает с соответствующей частью общего гамильтониана, использованного ранее (см. выше). В этом случае уравнения движения для , полученные из (1),
имеют вид:
(2)
где произведена замена .
В случае в системе (2) можно перейти к безразмерному дифференциальному уравнению синус-Гордона:
, (3)
”непрерывный аналог” системы (2). Это уравнение имеет солитонные решения, в частности, односолитонное решение, или кинк, характеризующий динамику распространения дислокации в цепи.
В соответствии с (1) система нелинейных уравнений движения записывается следующим образом:
(4)
Как видим, системы (2) и (4) существенно различаются. Отметим, однако, что проведенное нами численное моделирование динамики систем (2) и (4) показало следующее: если в качестве начальных условий для численного интегрирования (2) выбрать односолитонное решение его “непрерывного аналога” (3) - кинк (см. выше), то обнаруживается принципиальное сходство в характере решений.
Однако, при задании начальных условий в следующем виде:
(5)
где - ”ступенчатая” функция с высотой ступени и углом наклона уступа A, выявилось различие динамики данных систем (срав. рис.1 и 2,3). Более точно, системы (2) и (4) численно интегрировались методом Рунге-Кутта четвертого порядка с начальными условиями, заданными в виде (7), в интервале с шагом . Граничные условия - “квази-циклические”:
(поли-A-последовательность). Параметр системы . Варьировался параметр A (угол наклона уступа функции ).
Численное интегрирование системы (2) ( рис. 1) показало, что образуются две уединенных волны, движущихся справа налево по цепи с постоянной скоростью. Первая волна имеет форму квазикинка, а вторая волна имеет форму квазибризера, причем скорость первой волны превосходит таковую для второй. Обе волны за счет “квазициклических” граничных условий, доходя до левого конца, появляются на правом конце без изменения своей формы. Квазикинк, проходя по цепи маятников, изменяет координату каждого маятника на угол (маятник делает полный оборот). Поэтому, проходя по замкнутой цепи маятников К раз, он изменяет координату каждого маятника на угол Этим объясняется “уступообразная” форма графика на рис. 1.
На рис. 2 представлены результаты интегрирования системы (4) при тех же условиях. Из рисунка видно, что образуются те же две уединенных волны - квазикинк и квазибризер. Но принципиальное отличие от рассмотренного случая состоит в том, что квазикинк в самом начале движется с отрицательным ускорением, так что в результате его скорость оказывается меньше скорости квазибризера. Заметим, что исследования проводились на однородной поли-A-последовательности; так что изменение скорости квазикинка нельзя объяснить влиянием неоднородности цепочки. Этот эффект объясняется нелинейным взаимодействием между ее мономерами.
Рис. 3 иллюстрирует результаты интегрирования системы (4) при тех же условиях за исключением того, что A=2. В данном случае реализуется только квазикинк и его отрицательное ускорение в начале движения таково, что в результате он движется в направлении, противоположном первоначальному. При интегрировании системы (2) в аналогичных условиях также образуется только квазикинк. Его скорость не меняется по сравнению со случаем рис. 1.
Существенно, что при соответствующих условиях в системе типа ДНК или РНК могут возникнуть перевзбужденные ровибронные состояния. На квантовом языке это было бы адекватно перезаселению высоко лежащих квантовых уровней по сравнению с основным (реализации инверсной заселенности). В этом случае возникает заманчивая мысль, связанная с принципиальной возможностью создания биосолитонного лазера (БСЛ) на молекулах ДНК.
Однако, в теории динамики биополимеров хорошо известно, что конформационные движения реализуются по механизму ограниченной диффузии ввиду сильного влияния диссипативных сил со стороны микроокружения. По этой причине решение проблемы создания БСЛ на ДНК представляется весьма проблематичным, по крайней мере, для подтверждения идеи необходимо выполнение условий: где и - ширина и скорость солитона соответственно, - время диссипации. Положив 5 A и (скорость звука), имеем оценку . Отметим, что характерное время диссипации за счёт водных гидродинамических сил а время затухания, обусловливаемое процессами внутри самой молекулы (см., напр., Шайтан К.В. Биофизика. М.,1994. Т.39. С.949.; Чернавский и др. 1986. № 287. С. 21.).
Существует также и другая сложность в отношении самосогласования биосолитонов и волны электромагнитного переизлучения. Напомним, что математическое моделирование в данном случае проводилось на монотонной поли-A ДНК и поэтому оставалось неясным влияет ли гетерогенная естественная последовательность ДНК на динамику солитонного возбуждения в молекуле. Чтобы проверить это, как и ранее, был взят С-район ДНК на 3’-конце вируса саркомы птиц в качестве полигона для запуска солитонов на разных участках полимера. На этот раз вычисляли производную от функции с тем, чтобы нагляднее показать движения солитонов.
На рис.5,6 (см. ниже) хорошо видно, как при сдвиге области возбуждения солитонной волны от правой нижней части графика налево траектория волны претерпевает существенные изменения, т.е. “словесно-речевое” наполнение ДНК отображается в поведении солитона. Но главное здесь не только и не столько в этом. На этот раз характерно не качание волны около некоторого положения равновесия, а движение ее в левую часть цепочки после определенного временного интервала. В этом видится определенный биологический смысл. Солитон как потенциальный “субъект чтения” ДНК должен “просматривать” протяженные контекстные зоны, а не застревать на одних и тех же “словах”.
а) б)
Рис.5
а) Результаты численного моделирования динамики распространения возмущений в ДНК на основе системы (2) при значении параметра A=1.
б) То же, вид сверху.
а)
б)
Рис.6
а) Результаты численного моделирования динамики распространения возмущений в ДНК на основе системы (4) при значении параметра A=1.
б) То же, вид сверху.
a)
б)
Рис.7
а) Результаты численного моделирования динамики распространения возмущений в ДНК на основе системы (4) при значении параметра A=2.
б) То же, вид сверху.
200-ый
Рис.8
Солитонное возбуждение ДНК, но с учётом нелинейности ковалентных связей в сахаро-фосфатном остове ДНК. Последовательность нуклеотидов - вирус саркомы птиц (первые 600 пар оснований). Центр возмущения - 200-ый нуклеотид.
400-ый
Рис.9
То же, что на рис. 8, но центр возмущения - 400-ый нуклеотид.
500-ый
Рис.10
То же, что на рис. 9, но центр возмущения - 500-ый нуклеотид.
ЭКСПЕРИМЕНТАЛЬНЫЕ ДОКАЗАТЕЛЬСТВА СОЛИТОНООБРАЗОВАНИЯ НА ИНФОРМАЦИОННЫХ БИОПОЛИМЕРАХ “IN VITRO”
Способны ли молекулы ДНК и белков к солитонным возбуждениям, предсказанным в многочисленных теоретических моделях? Нами предприняты попытки фиксации нелинейных волн такого рода in vitro методом спектроскопии корреляции фотонов. Выявлены устойчивые эффекты, которые по ряду признаков соответствуют, в частности, процессу спонтанного солитонообразования в рамках явления возврата Ферми-Паста-Улама. [8,19,25,31,32]. Обнаружилось, что при переходе от разбавленного раствора ДНК к полуразбавленному можно зарегистрировать аномально долго затухающие акустические колебания гелевого континуума ДНК. Слабо затухающие колебания исчезают по мере перехода от полуразбавленного к разбавленному раствору и в результате уменьшения длины фрагментов ДНК. Эти данные подтверждают ранние работы для агарозы и коллагена, где впервые обнаружен феномен аномально слабой затухаемости плотностных колебаний биогелей. Аномальное поведение ДНК зарегистрировали после наших наблюдений и японские авторы методом прямой регистрации броуновской динамики флуоресцентно-меченой ДНК. Причем, в работе японцев выявились и другие необычные особенности нелинейной динамики ДНК, которые не укладываются в хорошо разработанные теоретические модели Цимма и Роуза, но которые хорошо соответствуют нашим наблюдениям и трактовке молекул ДНК как структур, резонирующих на особые внешние волновые регуляторные сигналы [25,6,7,15,16,29] (также см. ниже). Такая самоорганизация волновых процессов в ДНК может происходить и при таких физических условиях, когда существенную роль играют кооперативные процессы на уровне макромолекулярного континуума молекул ДНК, приближающегося к структуре хромосом. Чем более структура растворов ДНК отличается от архитектоники ДНК в хромосомах (в приводимых нами экспериментах это относительно короткие фрагменты полимера), тем менее существенны коллективные дальние (в масштабах макромолекулярных протяженностей полинуклеотида) взаимодействия между цепями ДНК, столь важные для эпигенетических функций генома. Ключевым звеном в данных экспериментах является четкая регистрация поведения ДНК in vitro, которое ранее зафиксировано Бреннером и Носсалом для агарозы и коллагена в аналогичных условиях. Это позволяет рассматривать нелинейную динамику такого рода для ДНК и других информационных биополимеров как проявление солитонных свойств в рамках явления возврата Ферми-Паста-Улама (ФПУ). Нелинейная динамика ДНК, ее гидродинамическое поведение и акустика чрезвычайно чувствительны к внешним физическим воздействиям in vitro - энзиматической рестрикции, разбавлению-концентрированию, нагреву-охлаждению, ультразвуковой обработке, слабым механическим воздействиям, облучению ИК-лазерным полем, излучением ФПУ-генератора с широкополосным электромагнитным спектром. Эти и аналогичные воздействия могут и должны в той или иной мере оказывать влияние на генетический аппарат в условиях in vivo, искажающее нормальные эпигенознаковые функции хромосом, что также подтверждается в наших экспериментах. Нелинейная динамика ДНК обнаруживает и другие “аномальные” свойства. Мы зафиксировали резкое различие коэффициентов диффузии для кольцевых и линеаризованных плазмидных ДНК [33], которое также не укладывается в циммовскую теорию поведения полимеров в водных растворах и в этом плане находит подтверждение в работах группы Роберта Пекоры (США) и упоминавшемся исследовании Матсумото с соавторами. Эти необычные свойства ДНК, вероятно, играют важную роль, например, для понимания механизмов управляемого “пилотирования” и точной “посадки” транспозонов ДНК (аналогов плазмид) в пределах жидкокристаллического сверхвязкого и сверхплотного континуума хромосом. Эта задача находится в области общей и нерешенной проблемы молекулярной биологии - проблемы самоорганизации внутриклеточных, межклеточных и межтканевых структур, их “взаимоузнаваний”. Ясно, что, зная волновые, гидродинамические и иные механизмы точного пилотирования таких немаловажных для человека транспозонов, как онкогены и обратнотранскриптазный геном вируса иммунодефицита человека, мы будем иметь возможность корректировать их в необходимом направлении, исключающем патогенез. Не менее существенным представляется факт обнаружения нелинейной динамики ДНК с признаками поведения солитонов по типу явления возврата ФПУ. Это также дает вклад в осознание принципов макромолекулярных и надмолекулярных взаимоузнаваний в организме по линии солитонно-резонансных дальних взаимодействий и делает более реалистичной попытку дать новую версию работы генома эукариот, обсуждавшуюся выше. Мы обнаружили и другие необычные проявления физических свойств ДНК - ее последействие или следовую память [25]. Этот феномен ставит проблему новых типов геномных функций. Возможно, это явление тесно связано с особой памятью генома высших биосистем, а также, вероятно, и с памятью коры головного мозга. Но если для ассоциативной корковой памяти и памяти генома растений нами и другими даны физико-математические модели в терминах и понятиях голографических и солитонных процессов, то память последействия ДНК - явление далеко не ясное и нуждающееся в более глубоком исследовании и осторожной трактовке. Этот эффект зарегистрирован нами при динамическом лазерном светорассеянии на препаратах высокоочищенных ядер из эритроцитов кур и на высокополимерной чистой ДНК из зобной железы теленка [25].По сути, аналогичное явление наблюдала группа Р.Пекоры (США) и назвала его “MED-effect” (Mimicing Effect of Dust), т. е. эффект, имитирующий пыль. Так же как и в наших работах, этот феномен обнаружен методом корреляционной лазерной спектроскопии на рестриктных фрагментах ДНК строго определенной длины. И в этом случае ДНК вела себя “аномальным”образом: зондирующие фотоны дифрагировали не только на полинуклеотидных цепях, но и на “посторонних частицах”, которых в препарате заведомо не было, что обеспечивалось специальным обеспыливанием. Этот никак не прокомментированный группой Р. Пекоры эффект сильно затруднил ей попытки объяснить поведение ДНК с позиций казалось бы хорошо разработанной теории Цимма и Роуза для динамики полимеров в водных растворах. И это еще раз было подтверждено в Японии Матсумото и др. прямым наблюдением “аномально” броунирующей флуоресцентно-меченой ДНК. Представляется, что в работе группы Пекоры cветорассеяние происходило не только на реальных фрагментах ДНК, но и на волновых следовых структурах ДНК, оставляемых броунирующими молекулами этого суперинформационного биополимера в духе теории физического вакуума, где постулируется идея генерации фантомных торсионных аксионно-кластерных эквивалентов физических тел.
Страницы: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15
