400-ый

Рис.1

Влияние нуклеотидной последовательности ДНК на динамику конфор-мационного возмущения уединенной (солитоноподобной ) волны. Последо-вательность нуклеотидов - вирус саркомы птиц (первые 600 пар оснований). Центр возмущения - 400-ый нуклеотид.

450-ый

Рис.2

То же, что на рис.1, но центр возмущения цепочки ДНК на 450-ом нуклеотиде.

Так работает компьютерная модель динамики солитонов, в определенной мере развитая Салерно после ее выдвижения Инглендером. Салерно дал формализм, описывающий вращательные колебания нуклеотидов молекулы ДНК, для того чтобы объяснить экспериментальные данные по водородно-тритиевому обмену в ДНК. Согласно этой модели по Инглендеру, в цепи ДНК могут возникать (под воздействием теплового шума) и распространяться открытые состояния (“плавление” двойной спи-рали ДНК на коротких участках, обогащенных АТ-парами ) в виде локализованных дислокаций ( уединенных волн). Марио Салерно, про- должая работу Инглендера, в упрощенном варианте выявил влияние последовательности нуклеотидов на нелинейную динамику вращательных колебаний нуклеотидов на однотяжных участках ДНК, образующих такие открытые ("open state") области. Позднее Якушевич, Федянин, Хомма и др. рассмотрели различные обобщения модели Инглендера, с оценкой особенностей строения ДНК, учитывая обрыв водородной связи при открытии оснований, парность цепи ДНК и другие степени свободы, отличные от вращательных. Однако, в указанных работах недостаточно сказано о причинах возникновения дислокаций в ДНК. Мы предлагаем возможный механизм этого процесса в ДНК, альтернативный гипотезе Инглендера о воздействии теплового шума как причины раскрытия пар оснований. Мы считаем, что дислокации на ДНК могут возникать при изменении периода спирали ДНК (основная часть идеи принадлежит М.Ю.Маслову).

В нашей модели нуклеотиды ДНК рассматриваются как осцилляторы, подвешенные на невесомом нерастяжимом стержне; сахаро-фосфатная связь между соседними нуклеотидами в цепи моделируется линейными пружинами; спирализация вдоль цепи не учитывается; водородные связи между комплементарными основаниями моделируется “гравитационным” потенциалом. Гамильтониан по М. Салерно выглядит следующим образом: (1)

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

Коэффициенты в уравнении (1) определяются в соответствии с правилом: в случае АТ и ТА пар, в случае ГЦ и ЦГ пар; - параметр, определенный Федяниным и Якушевич и полученный на основе модели синус-Гордона и экспериментальных данных. Далее для упрощения модели считается, что

Уравнения движения для разности , полученные из (1), имеют по М. Салерно вид:

(2)

где произведена замена .

В случае , в системе (2) можно перейти к безразмерному дифференциальному уравнению синус-Гордона:

, (3)

”непрерывный аналог” системы (2). Это уравнение имеет солитонные решения, в частности, односолитонное решение, или кинк, соответствует дислокации в цепи.