Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Нелинейная физика ДНК..pdf
Скачиваний:
9
Добавлен:
15.11.2022
Размер:
16.39 Mб
Скачать

Решения (5.123)—(5.125) показаны схематически на рис. 5.12. Их мож­ но интерпретировать следующим образом: (1) при прохождении кинка псевдовращательный фазовый угол (Р) меняется, скажем, от С2' эндо до СЗ' эндо; (2) в той же области две цепочки растягиваются, в особен­ ности в окрестности центра кинка; (3) углы между парами оснований уменьшаются внутри кинка.

Рис. 5.12. Изображение солитоноподобных решений модели Крумхансла и Алекзандер. Воспроизведено из ссылки [20] с разрешения авторов

5.4.2. Модель Волкова

Чтобы описать конформационные переходы в ДНК, Волков [54] предложил модель, включающую две группы внутренних движений. Первая состоит из движений, связанных с конформационными измене­ ниями внутри мономерного звена двойной спирали. Вторая включает из­ менения положения мономерного звена как целого. Предполагается, что

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

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

Гамильтониан такой двухкомпонентной цепочки равен

I tl + mrl + до (т„ -

T „ _ I )2 + / 0 (гп - rn-i)~ +

 

(5.126)

+ Ф (r„) + kF (г„) (тп+1 -

тп_!) /2

Здесь суммирование производится по всем мономерам внутри цепочки (индекс п); т п приведенная масса нуклеотидного звена

тп = т ъ р г п ъ ь / ( т п Ьр + т п Ьъ ) ,

(5.127)

где тьр — масса пары оснований; тпъъ — масса фрагмента сахарофос­ фатного остова тпьр + гпьь — масса нуклеотидного звена; I — момент инерции нуклеотидного звена относительно оси двойной спирали; до и /о — силовые константы взаимодействий вдоль цепочки (вдоль коор­ динаты z); к — силовая константа взаимодействия между подсистемами г и т; Ф — двухъямная потенциальная энергия конформационного пере­ хода в мономерном звене (рис. 5.13)

Ф (г) = ео + £1 {Г/а) + £ (1 - г2/а 2)2

(5.128)

a F (г) характеризует структурные переходы между подсистемами и вза­ имную ориентацию движений в обоих компонентах

F (г) = {г%- г2) /а2,

(5.129)

где г определяется условием

Рис. 5.13. Форма потенциальной функции конформационного перехода мономер­ ного звена в модели Волкова

а £, £о и £\ — параметры, которые определяют форму двухъямного по­ тенциала, показанного на рис. 5.13.

В континуальном приближении гамильтониан Н можно записать

как

 

 

Н = J d z^I т2 + s\ {дт/dz)2j

+ та2 \^у? + s\ (du/dz)2j +

 

 

(5.131)

+ Ф (и) + 2IkoF (и) dr/dz } /2 h,

где u(z,t) = r(z,t)/a;

= gob?/I\

s\ = foh?/m; ко = kh/I; h — рас­

стояние между мономерами вдоль цепочки. Соответствующие уравнения движения принимают тогда вид

f = s? (dr2/dz2) + к0(dF/du) (du/dz);

(5.132)

u = s\ (<д2и/дг2) - (ЭФ/ди)/2та2 - k0I (dF/du) (dr/dz) /та2

А солитонные решения этих уравнений, найденные Волковым, равны

т = т„ ctg [q(z - v t - zo)\/ {l - R2ctg2 [q (z - v t -

z0)]} ;

(5.133)

и = Щ —6 —26/ {R2 ctg2 [q (z — vt —zo)] —1} ,

 

(5.134)

rv = k0b (2r0 + ab) /q (SQ - v2) ;

R 2 = (2u0 + 6) /

(2u0 -

b);

6 = [2 (ul -

1)

/

(6 - 1)]1/2

q = ± p 1/27-1/2;

 

p = UQ (26 -

3)

+

1;

7 = та2 (у2 -

s|) /е;

 

6 = kli/e (s2 -

v2) ;

 

 

 

величина UQ определяется условием: du/d(z —vt) = 0 при и = UQ.

Рис. 5.14. Солитонные решения модели Волкова. Воспроизведены из ссылки [54] с разрешения автора

Решение уравнения (5.133) имеет форму кинка, а решением уравне­ ния (5.134) является колоколообразная функция от z —vt (рис. 5.14). Оба решения были найдены Волковым и интерпретированы им как решения, описывающие распространение локальных переходов А —►В типа вдоль двойной цепи ДНК.