Топологические солитоны в двухцепочечной модели ДНК Текст научной статьи по специальности «Физика»

УДК 577.3.01

ТОПОЛОГИЧЕСКИЕ СОЛИТОНЫ В ДВУХЦЕПОЧЕЧНОЙ МОДЕЛИ ДНК

Ф.К. Закирьянов, С.Р. Хамзин, К.Р. Юлмухаметов

Интенсивные исследования последних десятилетий показали, что в нелинейных системах различной физической природы возможны коллективные возбуждения в виде нелинейных уединенных волн - солитонов. Уединенные волны обладают рядом замечательных свойств, сближающих их с частицами - сохранение формы и скорости движения после взаимодействия. Они активно применяются для описания различных кооперативных эффектов в нелинейных динамических системах, в том числе молекулярных цепочках, например, таких как ДНК [1, 2]. Особый интерес представляют квазилинейные молекулярные цепочки, характеризуемые как геометрическим ангармонизмом, так и ангармонизмом потенциала междоузельного взаимодействия. Учет этих факторов приводит к моделям, описываемым системами нелинейных дифференциальных уравнений, решение которых крайне сложно получить аналитическими методами. Это приводит к необходимости разработки методов компьютерного эксперимента для исследования динамики подобного рода систем.

В работах [3, 4] изучена динамика топологических солитонов, описывающих образование открытого состояния двойной спирали ДНК. Показано, что в двойной спирали ДНК могут появляться четыре типа топологических солитонов (см. рис. 4 из [3]). Изучены взаимодействия между солитонами, их взаимодействие с неоднородностями цепи и устойчивость солитонов по отношению к тепловым колебаниям. Основным методом работ [3, 4] был вычислительный эксперимент с использованием дискретной модели ДНК, основанный на разработанной авторами вариационной технике для нахождения солитонов. Вместе с тем такие же решения можно получить и аналитически для модели цепочек в континуальном приближении.

В данной работе использовали дискретную модель ДНК [3, 4], основанную на модели Инг-лэндера [2]. На рис. 1 схематически представлена часть данной модели. Белые и черные прямоугольники соответствуют парам оснований АТ или ОС, а линии - скелету двойной спирали. Гамильтониан для этой модели имеет следующий вид:

Здесь (рп [ и <рп 2 - углы поворота азотистых оснований в я-х узлах первой и второй цепочки соответственно, /, и /2 - моменты инерции оснований первой и второй цепочки соответственно, ех и е2 - коэффициенты упругой связи, кх и к2 - жесткость взаимодействия соседних оснований в одной и разных цепочках соответственно.

В основе численного метода поиска солитонного решения лежит нахождение условного минимума для лагранжиана системы с соответствующими граничными условиями [3, 4]:

Найденные точки минимума служат начальными точками для решения системы уравнений движения. Решалась система дифференциальных уравнений второго порядка (задача Коши)

Рассмотрена математическая модель нелинейной динамики ДНК. Получены солитонные решения, описывающие конформационные изменения в процессе функционирования ДНК. Проведено сравнение аналитических решений с результатами численного эксперимента.

Ключевые слова: ДНК, нелинейная динамика, топологические солитоны.

+ <р„л) ,

где V((рп Х, (рп1) = кх (1 - соэ (рп Х) + кх (1 - сое <рп 2) +к2(\- сое(српХ - <рп2)).

Ф\,\ — Р-00,1 ’ 9\,2 “^-00,2

_ дН{<Рп,1<Рпл) ЯЛ л2 3<РпХ ’

СІГ

йг(Рп,г дн(<Рп,\<Рпл) "’2 Л2 д<р 2

(п = 0, ±1, ± 2...).

Интегрирование уравнений движения проводилось с помощью неявного метода Розенброка 2-го порядка точности.

Кроме численного анализа дискретной модели нами был проведен и теоретический анализ возможности получения решений с использованием теории возмущений. В континуальном приближении с учетом дисперсии длинноволновых колебаний динамика модели описывается системой двух нелинейных дифференциальных уравнений в частных производных:

Рис. 1. Фрагмент двойной спирали ДНК, состоящей из трех пар оснований вида <ЗС или (АТ). Расстояние между основаниями вдоль цепи а = 3.4 А, расстояние между цепями /1 = 16.15 А [3,4]

ЬФхтт ~е\

( 4

а

а4 2

^Фхуууу “у ЧЬ/Ьуу + а <Ріуу '

4^2

тт є2

\ V 12

+ кх + к2 вІпС^і -<р2) = 0,

<Рг

УУУУ

4 >

~(РІу(Р2уу + а2(р2уу

+ кх БІП (р2 - к2 8Іп(<£>| ~(р2) = 0 .

Коэффициенты перед членом с высшей дисперсией <Р\уууу и членом, содержащим нелинейность по пространственным производным, (р2у(р2уу при малых скоростях являются малыми. С увеличением скорости характерный размер кинков, определяемых уравнением вт-Гордон, уменьшается. В этом случае роль указанных слагаемых будет возрастать и особенно сильно проявляться именно вблизи предельной скорости солитона.

Для упрощения теоретических выкладок нами рассмотрена модель однородной молекулы ДНК, состоящей из двух одинаковых цепочек, где для простоты полагается 1Х = /2 = I, ех = е2 = е. Система классических уравнений движения после перехода к новым временным и пространственным переменным путем замены т = 1^1/кх и у = х^еа/кх принимает вид:

<Ры -<Р\хх-Т^-ФХхххх + ^-<Ри<Р\хх + 8ІП<3, + -<р2) = 0 ,

К",

12є К\

Фіи Фіхх у2.є ^хгаг

2є

+ ^є(Ріх<Р2хх + 8ІП9г-а-<р2) = 0’

где безразмерный параметр а = к2/кх.

Покажем теперь, что рассматриваемая система уравнений допускает аналитические решения, подобные полученным в работе [3] в ходе численного эксперимента.

Решения будем искать в виде ряда

<Рх = (р0 + Щ + ¥2 + •••> <Р2 = <Ро + П + ^2 + где % - решение системы, не учитывающей взаимодействие цепочек. Без учета взаимодействия цепочек (при а = О) рассматриваемая система приводит к известному уравнению вш-Гордон, решение которого в виде 2к -кинка имеет вид

% = 4акЛ§ехр00,

где г = (1-£2)_1/2£, л- - скорость солитона, £,=х-$1. Подставляя решение в виде ряда в рассматриваемую систему, получим в первом приближении уравнения для возмущений у/х и V,:

1С

Щ - Ух + ¥х С08 П - ~ *Рйхххх

^-Уі+^СОв^Ц)-

кх

\2е

. +-

1 2

ІХХХХ ' 2Е ФЪхФЪхх

О,

где соэ^о = 1 - 2/сЬ2 2. Полученная система уравнений имеет следующее решение:

¥\

к.

24е(1-я2)2

\

\hz-ik №

сЬ* 2Ае{\-*2)2 _1{] + кг)сЪ

/'

(1 + ^2)сЬ^у

Й1 2 - гк На

с!к + -

15лг,

\2ei\-s ) сЬг

кп

(1к + -

15*:,

12е(1-.у ) сЪг

5 = 0.9

где С] и С2 определяются из условия разрешимости уравнений для возмущений второго порядка малости (в нашем случае Сх - 0 и С2 = 0). Эти константы определяют сдвиг кинков друг относительно друга. В нашем приближении этот сдвиг не удается получить аналитически. Таким образом, взаимодействие цепочек изменяет только форму кинка.

Общий ВИД полученного решения представлен 1200 на рис. 2 и качественно соответствует случаю (е) на рис. 4 [3]. Заметим, что при приближении скорости солитона к предельной ( з —»1) теория возмущений дает большие погрешности.

Солитонное решение задачи можно охарактеризовать топологическим зарядом ц = {д\, д2), где 9/ = (<Роо,, ~ Фоо,1)/2л- 0' = 1,2) является целым числом

(<7, =0, ±1, ±2,...). Для нахождения солитонного решения с топологическим зарядом ч нужно решить задачу на минимум с граничными условиями:

Р-аоД = ^-оо,2 = <Роо,\ = 2^1, (Рю 2 = 2щ2 .

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

8.00

4.00

0.00

-4.00

-8.00

-10.00

-5.00

0.00

5.00

10.00

Рис. 2. Вид стационарных солитонов с топологическим зарядом д = (1,1) при различных значениях безразмерной скорости солитона

1“..йк".ки

Рис. 3. Вид стационарного солитона с топологическим зарядом я = (1,1) при скорости в = 0,592. К = 0,234 Н/м. Ширина 27 звеньев в каждой цепочке. Тип оснований вС

40 80

Ll.sjjI.iLi

Рис. 4. Вид стационарного солитона с топологическим зарядом я = (1,1) при скорости в = 0,804. К = 0,234 Н/м. Ширина 25 звеньев в каждой цепочке. Тип оснований вС

Видно, что происходит сдвиг кинков между двумя цепочками, то есть в одной из цепи кинк уходит немного вперед относительно другого кинка. Анализ показывает, что такой сдвиг происходит не при всех параметрах цепочки и скоростях уединенной волны, в частности, можно задать параметры таким образом и выбрать скорость солитона так, что сдвиг будет отсутствовать на всем протяжении движения солитона по цепи (рис. 4). Данный сдвиг объясняется неоднородностью цепи - у разных пар азотистых оснований различные параметры взаимодействия, масса и инерция, которые влияют на повороты углов в цепи. В случае теоретического исследования сдвиг отсутствует, так как рассматривалась только однородная цепочка (рис. 2). Представленная на рис. 4 динамика кинков без сдвига, имеет хорошее соответствие с аналитическими результатами.

Зависимость энергии топологического солитона я = (1, 1) от его скорости, полученная аналитически, представлена на рис. 5 (Е\). В ходе численного эксперимента была получена зависимость безразмерной энергии топологического солитона (Е) от безразмерной скорости (л) для всех трех интересующих нас типов топологических зарядов я = (1, 1) , я = (1,0) и я = (0, 1), которая представлена на рис. 6.

800-

600-

400-

200

кДж/моль

“Г

"Г

0.00 0.20 0.40 0.60 0.80 1.00

0,1

0.2

0,3

0,4

0,5 0,6

0,7

0,8

- (1.0)-*-(0.1)-

чГ1)1

Рис. 5. Зависимости энергий стационарных солитонов с топологическими зарядами ч = (1, 1) - д = (1, -1) - Е2 и д = (1, 0) - Е3 от безразмерной скорости я (/(1/С = 0,1,0= 0,1)

Рис. 6. Графики зависимости безразмерной энергии топологического солитона (Е) от безразмерной скорости (5), для оснований вс

Видно, что графики зависимости безразмерной энергии от безразмерной скорости для топологических солитонов я = (1, 0) и я = (0, 1) практически совпадают, так же как и при аналитическом решении.

Для получения аналитических решений, описывающих солитоны с топологическим зарядом Я = (1, -1), используем разложение

<Р\ = 9о + ¥\ + ¥г + •••» 9г = ~9о + П + ^ + при этом в первом приближении получим для возмущений следующую систему уравнений:

У>\-У\+¥\ -—^Охххх

+ ^-9ох9охх +«8ш(2р0) = О,

2е

у1-у![+у1со5<р0--^(р0 12 е

'*** +^£(р^(рЬхх-а™(29о) = 0 •

Её решение в виде стационарных волн имеет вид

¥\ =

V, =

15ят,

12^(1-я2)2

-2а

15лг,

12е(1-52)2

+ 2 а

сЪг

а-

лг,

24е(1 -*2)2

а +

К

24е(1-я2)2

г ИЬг-гк

/1 ,2ч 1 ^

-да (1 + к )СЬ —

Ш г - гк ^{\ + к2)сЪ~-

е*ьс1к,

?,ьс!к.

В итоге <рх и <рг имеет солитонные решения с топологическим зарядом ч = (1, -1), представленные на рис. 7 и качественно соответствующие случаю (ж) на рис. 4 [3]. Зависимость энергии топологического солитона д = (1, -1) от скорости солитона представлена на рис. 5 (Е2).

Рис. 7. Вид стационарного солитона с топологическим зарядом ч = (1,-1) для значения безразмерной скорости в = 0,9

Решения, полученные численно, и представленные на рис. 8, имеют качественное сходство с аналитическими решениями (рис. 7). Отличие их связано с упомянутым выше влиянием неоднородности цепочек, приводящим к расхождению кинков с течением времени.

Рис. 8. Вид стационарного солитона с топологическим зарядом я = (1,-1), в разное время численного эксперимента (< - время в программе). Скорости 5 = 0,898. К = 0,234 Н/м. Ширина 32 звена в каждой цепочке. Тип основания вС. Солитон движется справа налево

Солитоны с топологическими зарядами ц = (1, 0) и ц = (0, 1) будем считать практически эквивалентными, поскольку они соответствуют одной физической ситуации распространения кин-ка по одной цепочке, сопровождаемого малым возмущением по другой цепочке. Отличие будет только в значениях энергий рассматриваемых солитонов. Решение рассматриваемой модели в данном случае ищется в виде разложения:

<Р\=<Ро+¥\+¥2+-> Ф2=у\+у2+-’

В первом приближении для возмущений щ и V] имеем систему уравнений:

1С /С о

¥\~¥\+¥\ СО*<Р0-Т^-Рохххх + ~<РОх<РОхх+а*'т(<Ро) = 0 ,

\2є 2є

1С ^

VI -і£+V, -^<р0хххх +-^<Р0Х(Р0ХХ -авіп^о) = 0 .

Полученная система уравнений имеет следующее решение:

¥і =

6е(1-52)2

&2-1к е,ьЛ'

(1 + *2)сЬ^~

(зЬ2іп(2сЬг)-гсЬг) .

6г(1-52)2 сЬг

Таким образом, взаимодействие цепочек изменяет только форму кинка, вид которого представлен на рис. 9, качественно соответствующем случаям (а-г) на рис. 4 [3]. Численное моделирование (рис. 10) указывает на хорошее качественное соответствие с полученными аналитически результатами. Зависимости энергии топологического солитона <\ = (1,0) от его скорости, полученные аналитически и численно представлены на рис. 5 (£з) и рис. 6.

Ф

6

Рис. 9. Вид стационарного солитона с топологическим зарядом ч = (1,0) для значения безразмерной скорости 8 = 0,9

О

1 і - ; і ;

г

і 1

і..

I

'*■ 1

п

40

!ч>1

80 ' >р2

100

Рис. 10. Вид стационарного солитона с топологическим зарядом я = (1,0) при скорости 5 = 0,901. К = 0,8714 Н/м. Ширина 19 звеньев цепи. Тип основания вС

Таким образом, теоретический анализ рассматриваемой модели однородной ДНК показал, что в системе могут быть решения в виде солитонов со следующими топологическими зарядами:

1) Я = (1, 0) (или ч = (0, 1)) (по одной цепочке распространяется возмущенный кинк, а по второй цепочке - малое возмущение), соответствующие локальным движениям, возникающим при транскрипции ДНК;

2) я = (1, 1) (по двум цепочкам распространяются в одном направлении и с одинаковыми скоростями два одинаковых возмущенных кинка), который соответствует переходу ДНК из одного конформационного состояния в другое;

3) я = (1, -1) (по двум цепочкам распространяются кинк и антикинк), который описывает процесс репликации ДНК.

Получены зависимости энергии топологических солитонов от безразмерной скорости 5. При этом энергетически наиболее выгоден солитон с топологическим зарядом я = (1, 0) (рис. 5, 6). Проведено сравнение результатов численного и аналитического решений. Результаты показывают, что решения, полученные аналитическим методом и численным экспериментом, имеют качественное сходство.

Авторы выражают искреннюю благодарность профессору М.А. Шамсутдинову за полезные обсуждения в ходе выполнения работы.

Литература

1. Yakushevich, L.A. Nonlinear Physics of DNA / L.A. Yakushevich. - Wiley, Chichester, 2004.

2. Englander, S.W. / S.W. Englander et al. II Proc. Natl. Acad. Sci. U.S.A. - 1980. - V. 77. -P. 7222-7226.

3. Ковалева, H.A. / H.A. Ковалева, A.B. Савин, Л.И. Маневич и др.// ВМС. - 2006. - Т. 48, № 3. - С. 454—471.

4. Yakushevich, L.V. / L.V. Yakushevich, A.V. Savin, L.I. Manevitch // Phys. Rev. E. - 2002. -V. 66.-P. 016614.

5. Якушевич, JI.B. Методы теоретической физики в исследовании свойств биополимеров / Л.В. Якушевич. - Пущино: ОНТИ НЦБИ АН СССР, 1990.

Поступила в редакцию 12 января 2008 г.

TOPOLOGICAL SOLITONS IN A DNA DOUBLE STRAND MODEL

A mathematical model of nonlinear DNA dynamics is presented. Solitonic solutions describing conformational transitions during DNA functioning are obtained. Analytical solutions are compared with numerical simulations results.

Keywords: DNA, nonlinear dynamics, topological solitons.

Zakir’yanov Farit Kabirovich - Associate Professor of the Theoretical Physics Department of the Bashkir State University.

Закирьянов Фарит Кабирович - доцент кафедры теоретической физики Башкирского государственного университета.

e-mail: [email protected]

Khamzin Salavat Rifovich - Graduate Student of the Theoretical Physics Department of the Bashkir State university.

Хамзин Салават Рифович - аспирант кафедры теоретической физики Башкирского государственного университета.

Yulmukhametov Konstantin Raisovich - junior research fellow of the Institute of computer research at the scientific research of the Ufa State Aviation Engineering University.

Юлмухаметов Константин Раисович - младший научный сотрудник ИКИ НИЧ УГАТУ (Институт компьютерных исследований при научно-исследовательской части Уфимского государственного авиационного технического университета).