Математическое моделирование солитонов на ДНК

Страница 5

(поли-A-последовательность). Параметр системы . Варьировался параметр A (угол наклона уступа функции ).

Численное интегрирование системы (2) ( рис. 1) показало, что образуются две уединенных волны, движущихся справа налево по цепи с постоянной скоростью. Первая волна имеет форму квазикинка, а вторая волна имеет форму квазибризера, причем скорость первой волны превосходит таковую для второй. Обе волны за счет “квазициклических” граничных условий, доходя до левого конца, появляются на правом конце без изменения своей формы. Квазикинк, проходя по цепи маятников, изменяет координату каждого маятника на угол (маятник делает полный оборот). Поэтому, проходя по замкнутой цепи маятников К раз, он изменяет координату каждого маятника на угол Этим объясняется “уступообразная” форма графика на рис. 1.

На рис. 2 представлены результаты интегрирования системы (4) при тех же условиях. Из рисунка видно, что образуются те же две уединенных волны - квазикинк и квазибризер. Но принципиальное отличие от рассмотренного случая состоит в том, что квазикинк в самом начале движется с отрицательным ускорением, так что в результате его скорость оказывается меньше скорости квазибризера. Заметим, что исследования проводились на однородной поли-A-последовательности; так что изменение скорости квазикинка нельзя объяснить влиянием неоднородности цепочки. Этот эффект объясняется нелинейным взаимодействием между ее мономерами.

Рис. 3 иллюстрирует результаты интегрирования системы (4) при тех же условиях за исключением того, что A=2. В данном случае реализуется только квазикинк и его отрицательное ускорение в начале движения таково, что в результате он движется в направлении, противоположном первоначальному. При интегрировании системы (2) в аналогичных условиях также образуется только квазикинк. Его скорость не меняется по сравнению со случаем рис. 1.

Существенно, что при соответствующих условиях в системе типа ДНК или РНК могут возникнуть перевзбужденные ровибронные состояния. На квантовом языке это было бы адекватно перезаселению высоко лежащих квантовых уровней по сравнению с основным (реализации инверсной заселенности). В этом случае возникает заманчивая мысль, связанная с принципиальной возможностью создания биосолитонного лазера (БСЛ) на молекулах ДНК.

Однако, в теории динамики биополимеров хорошо известно, что конформационные движения реализуются по механизму ограниченной диффузии ввиду сильного влияния диссипативных сил со стороны микроокружения. По этой причине решение проблемы создания БСЛ на ДНК представляется весьма проблематичным, по крайней мере, для подтверждения идеи необходимо выполнение условий: где и - ширина и скорость солитона соответственно, - время диссипации. Положив 5 A и (скорость звука), имеем оценку . Отметим, что характерное время диссипации за счёт водных гидродинамических сил а время затухания, обусловливаемое процессами внутри самой молекулы (см., напр., Шайтан К.В. Биофизика. М.,1994. Т.39. С.949.; Чернавский и др. 1986. № 287. С. 21.).

Существует также и другая сложность в отношении самосогласования биосолитонов и волны электромагнитного переизлучения. Напомним, что математическое моделирование в данном случае проводилось на монотонной поли-A ДНК и поэтому оставалось неясным влияет ли гетерогенная естественная последовательность ДНК на динамику солитонного возбуждения в молекуле. Чтобы проверить это, как и ранее, был взят С-район ДНК на 3’-конце вируса саркомы птиц в качестве полигона для запуска солитонов на разных участках полимера. На этот раз вычисляли производную от функции с тем, чтобы нагляднее показать движения солитонов.

Страницы: 1 2 3 4 5 6