CC BY
54
УДК 531/534: [57+61]
Российский Журнал
www.biomech.ru
ЛОКАЛИЗОВАННЫЕ МОДЫ ДИСТОРСИИ В СТРУКТУРЕ ДВОЙНОЙ СПИРАЛИ ДНК
О.Б. Наймарк
Институт механики сплошных сред Уральского отделения Российской академии наук, Россия, 614013, Пермь, ул. Академика Королева, 1, e-mail: naimark@icmm.ru.
Аннотация. Так как реальные механизмы некоторых биологических процессов, включающие дезоксирибонуклеиновую кислоту (ДНК), до сих пор неизвестны, биологи уделяют все большее внимание структурным аспектам с целью выяснить действительные механизмы репликации (повторение процесса), транскрипции (первого этапа реализации генетической информации в клетке, результатом которого является биосинтез молекул рибонуклеиновой кислоты (РНК) на матрице ДНК) и денатурации (процесса расщепления белковых молекул). В настоящей работе сделаны попытки на основе решения статистической задачи о поведении ансамбля дефектов предложить уравнение движения в виде специальной формы уравнения Гинзбурга-Ландау для параметров порядка, отражающего поведение локализованных дисторсий в структуре ДНК, и исследовать динамические свойства данной системы. Показано, что динамика локализованных мод дисторсии может быть связана с новым типом критических явлений - структурно-скейлинговыми переходами, которые характерны для неравновесных систем с «медленной динамикой». Установлен вид дополнительного параметра порядка для данных систем - структурный параметр скейлинга, который играет роль «эффективных температур» неравновесной системы. Проведенный анализ автомодельных решений выявил существование трех типов коллективных мод в ансамблях дефектов в соответствующих интервалах структурного параметра скейлинга. Показано соответствие характерных коллективных мод трем типам универсальных статистических распределений в суперстатистике неравновновесных систем с медленной динамикой. Обсуждается возможная связь динамики коллективных мод с механизмами транскрипции и репликации.
Ключевые слова: дезоксирибонуклеиновая кислота (ДНК), двойная спираль, структура, ансамбль дефектов, статистический анализ, параметр порядка, скейлинг, неравновесные системы.
Малые группы атомов, слишком малые, чтобы проявлять точные статистические законы, играют главенствующую роль в весьма упорядоченных и закономерных процессах внутри организма. Они управляют видимыми признаками большого масштаба, которые организм приобретает в течение своего развития; они определяют важные особенности его функционирования, и во всем этом проявляются весьма отчетливые и строгие биологические законы.
Э. Шредингер. Что такое жизнь?
© Наймарк О.Б., 2006
09806267
Введение
Знаменитая книга Шредингера [1], цитата из которой вынесена в эпиграф настоящей статьи, явилась одной из первых попыток провести физический анализ основных проявлений «жизни» и, в частности, свойств дезоксирибонуклеиновой кислоты (ДНК) - молекулы, которая кодирует информацию, необходимую для организма в течение всей его жизни для воспроизведения себя. Спустя 50 лет после открытия структуры двойной спирали Уотсоном и Криком [2], ДНК по-прежнему завораживает как физиков, так и биологов, которые пытаются разгадать ее замечательные свойства. Поскольку реальные механизмы многих биологических процессов, связанных с ДНК, до настоящего времени непонятны, биологи обращают все большее внимание на структурные аспекты с целью объяснения основных механизмов функционирования ДНК - транскрипции и репликации. По сравнению с другими биологическими объектами, например, протеинами, ДНК может рассматриваться как относительно простая «сложная система», так как обладает упорядоченной структурой, которую Шредингер назвал «биологическим кристаллом».
В настоящее время широко обсуждается вопрос о том, что исследования статической структуры биологических молекул недостаточны для познания их функций. Это особенно справедливо для ДНК, которая претерпевает значительные конформационные изменения в процессе транскрипции, то есть считывании генетического кода, или при репликации. Именно поэтому представляется важным исследование динамических свойств молекулы ДНК. Однако, чисто динамические исследования недостаточны, так как термические флуктуации играют исключительную роль в функционировании ДНК. Таким образом, представляется естественным развитие статистической механики ДНК.
Принимая во внимание тот факт, что симметрия ДНК при её формировании может быть связана с дефектами, эффекты дисторсии, обусловленные дефектами в структуре двойной спирали ДНК, могут играть существенную роль. В развиваемых в настоящее время моделях микроскопические параметры порядка структуры ДНК вводятся как характеристики так называемых открытых комплексов, представляющих собой локализованные дисторсии, распределенные на малых сегментах 2ё (приблизительно 15-20 оснований). В настоящей работе сделана попытка на основе решения статистической задачи о поведении ансамбля дефектов предложить уравнение движения в виде специальной формы уравнения Гинзбурга-Ландау для параметров порядка, отражающего поведение локализованных дисторсий в структуре ДНК, и исследовать динамические свойства данной системы. Показано, что динамика локализованных мод дисторсии может быть связана с новым типом критических явлений - структурно-скейлинговыми переходами, которые характерны для неравновесных систем с «медленной динамикой». Установлен вид дополнительного параметра порядка для данных систем - структурный параметр скейлинга, который играет роль «эффективных температур» неравновесной системы. Проведенный анализ автомодельных решений выявил существование трех типов коллективных мод в ансамблях дефектов в соответствующих интервалах структурного параметра скейлинга. Показано соответствие характерных коллективных мод трем типам универсальных статистических распределений в суперстатистике неравновновесных систем с медленной динамикой. Обсуждается возможная связь динамики коллективных мод с механизмами транскрипции и репликации.
1. Структура и физические свойства двойной спирали ДНК
Ориентируясь на цели, поставленные в данной статье, приведем некоторые сведения о ДНК как функциональном биологическом объекте.
ДНК в биологии
Существуют два основных класса биологических молекул - протеины и нуклеиновые кислоты. Протеины реализуют основные функции в живом организме, являясь активными механизмами и катализаторами химических реакций, инициируют перенос других молекул, в том числе кислорода (перенос гемоглобина), контролируют проникание ионов через мембраны, выстраивают некоторые жесткие структуры в клетках, активируют деятельность молекулярных «моторов», например, растяжение мышц.
Сценарий строительства протеинов дается генетическим кодом, который сохраняется в структуре ДНК. ДНК выстраивает геном организма и, например, человеческий геном состоит из 46 фрагментов, размещенных в каждой клетке и хромосомах. Эксперименты показывают, что каждая хромосома формируется одной молекулой ДНК, имеющей длину 4-5 см. Геном же, сохраненный в каждой нашей клетке, имеет длину около 2 м.
Структура ДНК [3]
ДНК является полимером, или, более точно, представляет собой два сплетенных полимера. Каждый из мономеров - нуклеотидов, составляющих данный полимер, включает три элемента: фосфатную группу, пятиатомную атомную группу (sugar ring) и основание, которое представляет собой сложную органическую группу.
Важнейшее наблюдение, которое привело Уотсона и Крика к знаменитому открытию структуры двойной спирали ДНК, было связано с пониманием, что основания объединяются в пары с помощью водородных связей и что пары, образуемые purine и pyrimidine, имеющие всегда один размер, дают регулярную структуру двойной цепи нуклеотидов, скрепляемой водородными связями.
Понимание того, почему две связанные цепи нуклеотидов формируют двойную спираль, связано с рассмотрением взаимодействия между группами, формирующими нуклеотиды. Пары оснований представляют собой вытянутые плоские группы, имеющие форму достаточно больших плато, скрепленных друг с другом жесткими сахаро-фосфатными связями, допускающими, однако, ротационную подвижность относительно связей.
Таким образом, основная топология группы и их взаимодействие определяют структуру ДНК как геликоидальную. Притом традиционно отмечается, что длина молекул всегда существенно превышает диаметр, что является основанием для модельных представлений, в которых молекула описывается как деформируемая струна, когда ее поведение не ассоциируется с изменением ее внутренней структуры.
Существенны также два типа сил, которые стабилизируют структуру двойной спирали:
• Водородные связи между дополняющими основаниями, которые, следуя Watson-Crick, объединяются в пары ^(аденин)-Г(тимин) и С(гуанин)-С(цитозин), скрепленные, соответственно, двумя и тремя водородными связями.
• Ковалентные связи, которыми основания соединены вдоль цепи ДНК. При этом существенно, что водородные связи являются значительно более слабыми по сравнению со связями вдоль цепи.
Рис. 2. Формирование открытого комплекса (transcription bubble) и механизм
считывания
В соответствии с результатами Watson-Crick структура ДНК, имеющая так называемую В-конформацию, обладает следующими признаками [3].
• ДНК образуется двумя полинуклеотидными цепями (рис. 1), обвивающими с правосторонним поворотом общую ось, образуя двойную спираль с диаметром ~ 20 А. Двойные цепи являются антипараллельными и обернуты друг относительно друга так, что не могут быть разделены без взаимного расплетения. Основания покрывают поверхность спирали, которая обвита сахаро-фосфатными цепями на периферии, причем таким образом, что достигается минимизация сил отталкивания между заряженными фосфатными группами.
• Плоскости оснований почти перпендикулярны оси спирали и каждая пара расположенных друг против друга оснований имеет водородную связь и образует при этом плоскую пару оснований. Идеальная В-ДНК спираль содержит 10 пар оснований (base pairs, pb) на один поворот (поворот спирали на 1 pb составляет 36°), что обеспечивает шаг спирали (rise per turn) и пространственную периодичность угла закрутки соответственно h = 34 А и ш0 = 2тс/h = 1,85нм_1. Ароматические основания имеют ван-дер-ваальсовский масштаб 3,4 А и, как следствие этого, наблюдается эффект частичного перекрытия (base stacking).
• Персистентные длины поворотных и изгибных мод представляют естественные масштабы для термических флуктуаций (приблизительно 15-20 h). ДНК является достаточно жестким полимером вследствие растяжения соответствующих пар аминокислот химическими и водородными связями. Термические персистентные длины ДНК на изгиб и вращение составляют, соответственно, 150 bp (50 нм) и 220 bp (75 нм).
• При пространственных масштабах в тысячи пар оснований (~ 0,5 мкм) ДНК in vivo организованы в топологически независимые петли. Молекулы ДНК с размерами 2-10 kbp представляют собой «связанные» множества. Основная часть генов (prokariotic cells) организована как единичные замкнутые ДНК с размерами ~2 Mbp.
Чтение генетического кода
Последовательность пар оснований в ДНК кодирует информацию для синтеза протеинов. Чтение генетического кода является достаточно сложной процедурой, так как основания, содержащие информацию, обнаруживают эффект перекрытия («укладки»). Как следствие, считывание кода предполагает наличие дисторсий, «проявляющих» основания, скрытые до этого за счет «укладки». Считывание реализуется молекулой рибонуклеиновой кислоты РНК (ribonucleic acid polymerase),
Рис. 3. Распад ДНК
которая сканирует длину молекулы ДНК, определяя место ее инициирования, то есть зарождения гена. После этого происходит «расплетание» двойной спирали на участке ориентировочно в 20 пар оснований и формирование дисторсии (transcription bubble) (рис. 2). Основания, которые себя проявили при образовании дисторсии, копируются молекулой РНК, которая выполняет функцию «template» для синтеза протеинов. Образовавшаяся «считывающая дисторсия» движется вдоль ДНК, замыкая основания после процедуры считывания и размыкая новые до завершения движения по двойной цепи.
Процесс транскрипции включает дисторсию молекулы, которая имеет настолько большую длину, что ей естественно проявить нелинейную природу силового взаимодействия различных молекулярных групп.
ДНК как динамический объект [4]
Открытие структуры двойной спирали подтвердило определение, что «форма есть функция», но очевидным является также тот факт, что чтение генетического кода предполагает большие конформационные изменения в структуре ДНК. Поэтому ДНК является в высшей степени динамическим объектом, в котором наблюдаются термически активируемые раскрытия участков ДНК, сопровождающие процесс транскрипции. Попытки описания динамики этих процессов предпринимаются последнее время с физических позиций и могут рассматриваться как первые шаги в этом направлении. На «дыхание» (breathing) ДНК биологи обращают пристальное внимание последнее десятилетие, отмечая, что это явление сопровождается временным раскрытием пар оснований. Эксперименты показывают [5], что средняя продолжительность «жизни» связанной пары оснований примерно 10 мс, что свидетельствует о существенных конформационных изменениях, которые приводят к локализации раскрытий пар оснований, сопровождающих механизм транскрипции.
2. Структурные модели ДНК
Структурные переходы в ДНК
Имеет место формальная аналогия между механизмами неустойчивости (melting instability) в структуре ДНК и фазовыми переходами в однородных системах. Эффекты плавления-распада (denaturation) в растягиваемой двойной цепи ДНК (рис. 3) являются одной из фундаментальных проблем в статистической физике биологических систем, которые обнаруживают различные типы критического поведения, наблюдаемые в других нелинейных системах [б].
Ранние модели термического распада ДНК, то есть разделения двух цепей при нагреве, были основаны на моделях типа Изинга, в которых параметр порядка, характеризующий взаимодействие оснований, сводился к расстоянию между последними. Современные исследования связывают эти переходы с нелинейными механизмами зарождения конечно-амплитудных флуктуаций [7-10]. Анализ закономерностей скейлинга указанных переходов был проведен в [7] с использованием асимметричного потенциала, содержащего пологий участок для больших значений перемещений, что соответствует распаду двойной спирали при повышенных температурах. Этот тип потенциала включал также взаимодействие между последовательными основаниями в длинноволновом приближении. В соответствии с этим описанием локальная неустойчивость плавления возникает в «мягком» режиме, что не подтверждается экспериментальными наблюдениями динамики переходов, которые носят резкий характер.
Обобщение данной модели было предложено в [8-10] введением параметра жесткости (“stacking parameters'’), который, учитывая взаимодействие последовательных пар оснований, позволил описать резкий характер перехода. Исследования однородных ДНК также показали, что данная модель дает возможность описания последовательного распада (плавления), наблюдаемого в эксперименте [11]. В дополнение показано, что в рамках модели возможна генерация «энтропийного барьера», с которым связывается «сужение» перехода. При некоторых значениях параметра жесткости модель демонстрирует признаки переходов первого рода с конечным значением энтропии плавления, нарушение непрерывности для части связей (что типично для ДНК), расходимость продольных и поперечных корреляционных масштабов [8].
Степени свободы
Степени свободы, описывающие дисторсии в молекуле ДНК, могут быть разделены на переменные, определяющие конформации центральной оси, и переменные, определяющие внутренний поворот двойной спирали. Положение оси ДНК может быть описано пространственной кривой, заданной параметром длины s вдоль оси ДНК и радиусом-вектором точки r (s). Касательная t = дsr задается единичным вектором, положение которого в пространстве определяется двумя углами или локальной кривизной и дисторсией [12].
Поворот измеряется как пространственная интенсивность поворота пар оснований относительно центральной оси, которая для неискаженной структуры ДНК задается как ш 0. Отклонение в интенсивности поворота от о 0 может быть описано с использованием скалярного поля Q(s). Это дает в целом три степени свободы для каждой координаты s , что необходимо для описания поворотов координатных осей, связанных с парами оснований при перемещении последних вдоль молекулы.
Peyrard-Bishop модель
Проблема распада (denaturation) ДНК вызывает большой интерес в контексте процессов, лежащих в основе механизмов транскрипции - важнейшего свойства молекулярных биологических систем. Эксперименты на разбавленных растворах ДНК демонстрируют существование термически инициируемых процессов денатурации, сопровождающихся разрывом связей пар оснований (opening of base pairs) и образованием локализованных мод дисторсии при некоторой критической температуре Tm [13]. Показано, что процесс денатурации может быть описан в рамках переходов первого рода при контролируемых внешних параметрах - температуре и внешнего
закручивающего момента. Для описания поведения однородной ДНК в [8] предложен потенциал
Н = 1
р2
+ Г (Уп, Уп-1) + V (у„) + ВЪауп 2т
(2.1)
где т - приведенная масса пар оснований; уп определяет трансляционную моду, обусловленную растяжением водородных связей, соединяющих два основания п -й пары, рп = т(ёуп/&). Связь между последовательными парами оснований описывается соотношением
W(Уп, Уп-1) = К I1 + Р ехр(- б(Уп + Уп-1))](Уп - Уп-1)2 .
Этот выбор потенциала взаимодействия обусловлен выводом, что энергия упаковки в основном определяется свойством пар оснований, а не свойствами одиночных оснований. Эффективная константа связи изменяется от значения (1 + р)К до К, когда одна из двух взаимодействующих пар оснований раскрывается. Это предположение находится в качественном соответствии с реальными признаками взаимодействия оснований в ДНК.
Третий член представляет собой одноточечный потенциал, который описывает взаимодействие двух оснований в паре:
V(Уп) = п(в~аУп -1)2.
I = щв
Последний член описывает влияние продольных внешних усилий И. Относительная простота постановки позволяет моделирование в широком временном интервале, что необходимо для оценки экспериментальных условий.
Оценка основных параметров модели принимается для единиц длины в ангстремах А и энергии в электрон-вольтах вV. За единицы измерений времени и силы приняты 1,021-10-14 с и 1602 пН, в соответствии с которыми параметры модели оценивались, как
Щ = 0,04вV, К = 0,04 е^А2, а = 4,5 А, а = 0,35 А-1, т = 300 а.е.м., р = 0,5, 5 = 0,35,
и давали значение температуры денатурации 360 К.
Вращательная мода в ДНК
Реакция молекулы ДНК на внешнее закручивающее напряжение интенсивно исследуется последние десятилетия [12, 14, 15] с целью более глубокого понимания механизмов денатурации. С биологической точки зрения появление закручивающих напряжений в живой клетке не является необычным и может существенно влиять на функционирование ДНК. Было показано, что денатурация может быть описана в рамках фазового перехода первого рода, когда контролирующими параметрами являются температура и приложенная величина внешнего вращения. Эта теория дает объяснение БВЫЛ-ёВЫЛ фазового перехода в экспериментах с единичными молекулами [16, 17].
Эта модель воспроизводит спираль Watson-Crick (ББЫЛ) и схематически представлена на рис. 4. В плоскости, перпендикулярной к оси спирали, для каждой пары оснований (п=1,..., Ы) вводится полярная система координат гп и фп. Сахаро-
Рис. 4. Модель спирали Уотсона-Крика
фосфатное основание моделируется жесткими стержнями с фиксированной длиной Ь, в то время как расстояние кп между плоскостями оснований (п -1) и п может флуктуировать.
Потенциальная энергия V в переменных гп, фп записывается в форме
V = Z [Vm (r„ ) + V, (rn, rn-1) + Vb (rn , rn-1, 0n ) + Vr (0„ )],
(2.2)
где 0n = фn -ф„_! - угол закручивания спирали.
Первый член Vm (rn) = D(e а(r R) -1)2 представляет энергию водородной связи для заданной пары оснований. В данной модели пара оснований с диаметром r > rd = R + 6/а рассматривается как открытая. Глубина потенциала D, обычно порядка 0,1 eV, зависит от типа пары оснований (АТ - аденин-тимин, или GC - гуанин-цитозин), а также от прочности ионной связи.
Сдвиговые усилия, которые оказывают сопротивление скольжению одного основания относительно другого, характеризуются потенциалом упаковки (stacking potential)
Vs (rn, rn-1) = Ee~Ъ{ГпГп- 12R) (rn - rn-1) 2 •
Вследствие уменьшения молекулярной упаковки с появлением открытых пар, сдвиговой фактор экспоненциально уменьшается и становится пренебрежимо малым на расстояниях, превышающих «5Ъ_1 = 10 А, которые близки к диаметру пар оснований.
Флуктуации аксиального расстояния hn между соседними парами оснований описывается упругим потенциалом
Vb (rn, rn-^ б n ) = K[hn-H ]2 = K yl L -rn2 -rn- 1 + 2rnrn-1 COS 0 n -H ^ •
Геликоидальные структуры возникают при H < L : в 5-конформации r < rd минимум потенциала Vb (= 0) достигается при ненулевом угле поворота (рис. 3).
Нелинейные локализованные моды дисторсии в структуре ДНК
Peyrard-Bishop модель позволила провести разнообразные исследования локализованных мод в структуре ДНК. Приближенное решение, полученное на основе данной структуры гамильтониана для одномерной ДНК решетки, показало
существование периодических и пространственно-локализованных решений, получивших название Peyrar-Bishop (РВ) бризеры. Эти коллективные моды исследуются применительно к механизмам термической денатурации ДНК и проявляются на кривой «плавления» в виде локализованных мод тонкой структуры, с динамикой которых связываются основные механизмы транскрипции. Эти локализованные моды, являясь устойчивыми в малоамплитудном пределе, могут распространяться вдоль молекулярной цепи в течение конечного интервала времени с незначительным изменением их энергии и без изменения формы. Имея малую амплитуду, движущиеся РВ бризеры играют роль предвестников, формирующих считывающие домены (transcription bubble) в термолизованной цепи со спонтанно локализованной энергией, и обеспечивают переход к активирующим дисторсиям и связанным с ними механизмами транскрипции и репликации. Бризеры обладают способностью к росту, аккумулируя энергию малых возмущений до возникновения больших доменов денатурации, и в завершение - к полному разделению цепей. Природа этих флуктуаций связана с широкой вариацией скоростей, с которыми растет число денатурированных пар оснований. Такие большие флуктуации, соответствующие «дыханию ДНК», наблюдаются экспериментально [18]. Вследствие экспериментов установлено, что изолированный бризер останавливается, если его фронт достигает «холодной» области. Таким образом, настоящий результат показывает, что разрывы водородных связей не являются чисто механическим процессом и важную роль играют термические флуктуации.
PB-модель служит удачным базовым приближением, позволяющим
интерпретировать ряд механизмов, связанных с динамикой энергетических флуктуаций в ДНК. Однако наблюдаемые при растяжении ДНК локальные вращения спиральной молекулы играют также важную роль в различных биологических функциях ДНК. Вращательная мода «constraining local twist» может быть введена как локальная степень свободы в рамках рассматриваемой «одномерной модели» ДНК.
Численный эксперимент, основанный на использовании PB-модели, позволил установить связь микроскопических механизмов денатурации с конечноамплитудными флуктуациями растягивающих усилий. Аналогичные сценарии нелинейной динамики были исследованы экспериментально и теоретически в [19] при возбуждении колебаний в струнах никелида титана при температурах и уровнях напряжений, соответствующих переходам в режим сверхупругости.
3. Эволюция дефектов и механизмы дисторсии в структуре двойной спирали ДНК
Последовательность переходов в структуре ДНК может быть рассмотрена как свойство эволюции дефектов при структурно-скейлинговых переходах в полимерном кристалле двойной спирали в соответствующем температурном интервале. Роль дефектов представляется важной при формировании полинуклеотидной цепи, если последняя рассматривается как непериодический одномерный кристалл, для которого наиболее устойчивое состояние может быть достигнуто вследствие поворотов, обусловленных дисклинациями, расположенными между соседними парами оснований. Именно с дисклинационной частью дисторсии связан эффект закручивания спирали (рис. 3) на определенный угол, приходящийся на шаг спирали, и, возможно, напряженное (stretching) состояние молекулы ДНК. Данный факт подчеркнул Watson при обсуждении устойчивости структуры двойной спирали BDNA и в [7-10] при обсуждении роли структурных деформаций в реализации биологических механизмов транскрипции, длиннокорреляционных активационных эффектах и, в том числе, при
формировании открытых комплексов. Например, имеются многочисленные данные по усилению эффектов в последовательностях геликоидальных мод, то есть в областях, где повороты между соседними парами оснований, или углы поворота, изменяются с учетом величин, характеризующих условия равновесия. Деформации, обусловленные вращениями, жестко связаны геометрическими соотношениями с условиями формирования открытых комплексов.
Локализованные моды дисторсии в цепи двойной спирали ДНК могут рассматриваться как дислокационные субструктуры. Увеличение дислокационной плотности (в случае «кристалла» двойной спирали ДНК - концентрации разорванных пар оснований) сопровождается уменьшением расстояния между дислокациями и увеличением взаимодействия (за счет роста структурных напряжений) в соответствующих дислокационных субструктурах.
Таким образом, коллективные свойства в дислокационных ансамблях начинают играть определяющую роль при формировании дислокационных субструктур новых структурных масштабов и их эволюции. «Движущей» силой переходов в дислокационных ансамблях является стремление системы достичь относительного минимума полной энергии при формировании дислокационных субструктур [20, 21]. Энергия дислокационных субструктур включает две части: собственную энергию -энергию дислокационной субструктуры и энергию взаимодействия последних. Изменение типа дислокационных субструктур ведет к изменению обеих частей и сопровождается, в целом, уменьшением энергии системы, содержащей новые дислокационные субструктуры.
Основная часть энергии дислокационных субструктур принадлежит собственной энергии дислокаций [22]
Ди = Р^21,А (3.1)
2л г0
где р - плотность дислокаций в дислокационной субструктуре, Ь - вектор Бюргерса,
G - модуль сдвига, г0 - радиус ядра дислокации, Ь - масштаб возмущения упругого
поля, создаваемого дислокацией. Последний масштаб играет важную роль в эволюции дислокационных субструктур: рост дислокационной плотности приводит к
уменьшению Ь при последовательном изменении типов дислокационных субструктур.
Некоторые результаты теории калибровочных полей
Теория калибровочных полей широко используется в настоящее время для анализа структурных и физических свойств сред с дефектами, и ее результаты могут быть применены к системам, содержащим дефекты, но не являющимися кристаллическими. Универсальность приложений теории калибровочных полей обусловлена фундаментальным признаком сред с дефектами - нарушением локальной симметрии среды вследствие изменения дефектами так называемой диффеоморфной структуры поля смещений. Используя формализм теории калибровочных полей (подход Янга-Миллса [22]), указанные изменения могут быть учтены как локализация соответствующей группы симметрии тензора дисторсии и рассматриваться как дополнительные внутренние переменные, кинематически согласованные с полем смещений среды.
Структура так называемых калибровочных полей, сопряженных переменным, связанных с дефектами, определяется типом описываемых дефектов. Основой для приложений теории калибровочных полей является установление внутренней группы симметрии для среды с дефектами. Так, основные механические модели
деформируемого континуума инвариантны относительно однородной группы трансляций T(3) и вращений SO(3), что означает инвариантность лагранжиана системы относительно преобразований
x ^ x' = g&t) + t^
где x и Е, - текущие и начальные координаты, g (^, t) и т(^, t) - операторы вращения и трансляции, соответственно.
Дифференциальные операторы должны быть также инвариантными по отношению к преобразованиям T(3) и SO(3), что учитывается видом ковариантной производной.
DХ = дцx + ГЦx + рц ,
где Гц и - дифференциальные операторы вращения и трансляции. Для среды с дефектами однородность преобразований T(3) и SO(3) нарушается, так как в этом случае локальные вращения и трансляции привносятся дефектами. Внутренняя группа симметрии такой среды определяется так называемым полупрямым произведением групп SO(3) T(3), которые, как функции координат и времени, задают структуру калибровочного поля.
Другим важным результатом теории калибровочных полей является установленная структура лагранжиана среды с дефектами. С учетом введенных дополнительных переменных Гц и лагранжиан может быть записан в виде
L = -2Qgggg-2С2g$g$, (3.2)
где
Q^v = - ^УГЦ + ГЦГУ - ГУГЦ ,
ацу = - ^vPv + ГцА - Ц/Рц + Q^vx ,
так называемые интенсивности калибровочного поля, относящиеся к локальным ротациям и трансляциям, привносимым дефектами; gik - компоненты метрического
тензора. Феноменологические параметры C1 и С2 в лагранжиане L являются константами дислокационного (дисклинационного) взаимодействия. Представление лагранжиана в форме (3.2) будет использовано в дальнейшем при статистическом анализе поведения ансамбля дефектов.
Параметры порядка структуры ДНК с локализованными дисторсиями
Структурные параметры, ассоциированные с трансляционными дефектами, были введены [21] в виде, аналогичном тензору дислокационной плотности. Микроскопический параметр порядка sik характеризует открытые комплексы (open complex), локальная дисторсия которых охватывает сегмент 2d (приблизительно 15-20 пар оснований) с нормальной трансляционной модой B = Bv к площади сегмента
Sd = Sd v (Sd = rod2)
sik = svtvk, (3.3)
где s = SdB.
Для сдвиговой моды микроскопический параметр порядка имеет вид
э1к = 1/25(уг4 + 1^к ), (3.4)
где I - единичный вектор в направлении сдвига; ^ = 8аБ - интенсивность сдвига с
величиной вектора Бюргерса Б = Б1 . Оба параметра порядка характеризуют трансляционные степени свободы для выделенной ячейки двойной спирали ДНК.
Дисклинационные моды ДНК могут быть представлены как антисимметричная мода дисторсии, индуцированная дефектами
*гк = 1/2^А - Ь^к ). (3.5)
Вышеприведенные степени свободы не являются независимыми, и поэтому
могут быть приведены физические аргументы для рассмотрения в первую очередь
поведения симметричных мод дисторсии в структуре ДНК.
Статистическая модель структуры ДНК с локализованными дисторсиями. Метод эффективного поля
Метод эффективного поля часто используется в теоретических моделях в виде вспомогательного (реального или предполагаемого) воздействия для описания сложных межчастичных взаимодействий, которые невозможно (или сложно) описать детально, или же когда характер этих взаимодействий неизвестен. В нашем рассмотрении определение «эффективного поля» предполагает введение вспомогательного многокомпонентного поля, сконструированного таким образом, что гамильтониан системы с учетом структуры эффективного поля соответствует условиям равновесия для данного текущего состояния. Впервые эта идея была предложена Леонтовичем [23] для определения термодинамических потенциалов неравновесных состояний и развита в [24] применительно к некоторым классам конденсированных сред. В соответствии с этим подходом произвольное неравновесное состояние термически однородной системы, которое описывается некоторыми параметрами порядка, может соответствовать этим же параметрам порядка в случае равновесия системы в присутствии некоторого эффективного силового поля. Следуя данному определению, энтропия этого неравновесного состояния будет соответствовать энтропии равновесного состояния в присутствии дополнительного силового поля и зависеть от тех же значений материальных параметров.
Микроскопическая кинетика параметров s. к определяется уравнением Ланжевена
% = Кгк (*1т ) -Ргк , (3.6)
где Кк = Е - энергия дефектов; ^к - случайная 5 -коррелированная часть
/ а^1к,
случайного усилия, удовлетворяющая соотношениям: ^к (1)) = 0 и
(^к (1 ')Ек (1)) = Q5(t -1') . Параметр Q характеризует дисперсионные свойства системы и
в нашем случае соответствует средней энергии зародышей дефектов. Принимая во внимание данные о роли исходной «скрученности» (зпретсоШп^) в структурных превращениях ДНК, величина Q может быть связана со значениями начальной скрученности цепей двойной спирали, что, в свою очередь, может соответствовать
предположению о существовании исходных дефектов дисклинационного типа между соседними парами оснований (stacking nucleotides) вдоль двойной спирали.
Статистическая модель поведения ансамбля дефектов была развита в [21]. В соответствии с методом эффективного поля функция распределения для неравновесного состояния, определяемая решением уравнения Фоккера-Планка
множитель. Как следует из (3.7), статистические свойства ансамбля дефектов могут быть описаны после определения энергии дефектов Е и дисперсионных свойств системы, определяемых величиной Q. С использованием микроскопических и макроскопических переменных и в соответствии с представлением открытых комплексов как дислокационных субструктур, энергия которых может быть записана в виде:
где квадратичный член представляет собственную энергию дефекта (2.2), член Hjksjk описывает взаимодействие между дефектами, а также внешним тензорным полем усилий ojk с использованием приближения эффективного поля, вид которого предполагается следующим
где а, X - константы взаимодействия, определяемые структурой двойной спирали. Среднее значение тензора дисторсии для открытого комплекса дается процедурой усреднения
где G - эффективная характеристика упругости двойной спирали при образовании открытого комплекса; У0 ~ г03 - характерный объем зародыша данного дефекта; Я -расстояние между открытыми комплексами. В результате мы получаем для 5 значение
W,
(З.7)
может быть представлена в виде
где Z - нормирующий
E = E0 - Hiksik +asi2k ,
(З8)
(З.9)
Pik = \SikW(^ V,/ )dsik .
(З.10)
Переход к безразмерным переменным
p lk
Q pik, Sik
VQa
позволяет представить уравнение самосогласования в виде
(З.11)
8 ~ I 1 , что отражает так называемую «статистическую автомодельность» в
V/ го )
поведении дефектов различных структурных уровней. Решение уравнения самосогласования (3.11) приведено на рис. 5, 6 [20].
Решение установило существование трех характерных нелинейных режимов, соответствующих различным областям значений параметра структурного скейлинга 8 (8 > 8* «1,3, 8с < 8 < 8*, 8 < 8с «1), разделенных точками бифуркации решений 8с и
8*.
Статистическое описание позволило предложить вид мезоскопического неравновесного потенциала, описывающего эволюцию открытых комплексов в зависимости от эффективных характеристик жесткости двойной спирали ДНК.
Представленные на рис. 5 зависимости соответствуют решениям уравнения = 0,
где ¥ - мезоскопический потенциал системы. Для одномерного случая
(р = р22, а = а^, 8 = 8^) «минимальное разложение» дается полиномом 6-го порядка по переменной р и имеет вид аналогичный разложению Гинзбурга-Ландау [20].
¥=12 А[%*) р2 _ % вр4+%с ("% ^р6 _ °ар+х^у 1р ^2. (312)
с
Точки бифуркации 8*, 8с играют роль, аналогичную критическим
температурам в теории фазовых переходов Ландау. Градиентный член в (3.12) описывает эффекты нелокальности при взаимодействии дефектов («диффузию» мод дисторсии - открытых комплексов); А, В, С, Б их - параметры, характеризующие структуру ДНК. Отметим, что с учетом вида энергии дислокационного дефекта (3.8) разложение (3.12) отражает вклад основных слагаемых в выражении для гамильтониана (2.1) при формулировке модели Peyrard-Bishop .
Вид разложения по 8 отражает роль данного параметра, как второго параметра порядка, величина которого связана, что будет показано в дальнейшем, с новыми условиями термолизации мезоскопической системы, которые не определяются обычной температурой. Роль последней заключается в формировании критических зародышей (начального значения 8 ), приводящих при соответствующих условиях к формированию открытых комплексов, обладающих различной динамикой в структуре
двойной спирали ДНК, с которыми могут быть связаны различные механизмы функционирования последней. Появление структурного параметра, который выполняет роль эффективной температуры, представляется важным, так как это обеспечивает независимость динамического поведения системы от обычных условий термолизации, которые инициировали процессы «медленной динамики» в рассматриваемых системах, например, формирование «бризеров».
Кинетика введенных параметров порядка соответствует эволюционному неравенству [20]
SF дF dp дF dS _
— =-------------£- +-------------------< 0.
St дp dt дS dt
(3.13)
из которого следуют уравнения движения для параметра дисторсии, характеризующего открытые комплексы, и параметра структурного скейлинга, характеризующего текущую восприимчивость структуры ДНК к формированию открытых комплексов, и, как это будет показано ниже, их динамику
dp
dt
= -г
%
-5 - Da-A (x^)
Л
Sxl дx
l У
(3.14)
(3.15)
*
Локализованные моды дисторсии, индуцированные дефектами
Как следует из решения уравнения (3.11), представленного на рис. 5, переходы через точки бифуркации 5С и 5* приводят к резким изменениям в функции
распределения и формированию коллективных мод тензора дисторсии pik. Тип перехода через критические точки определяется видом бифуркации - групповыми свойствами уравнений эволюции (4.9)-(4.10) в различных областях значений параметра структурного скейлинга 5 (5 >5*, 5С <5 <5*, 5 <5с). Качественные изменения в поведении системы представлены на рис. 7 в виде семейства гетероклиник, соответствующих решениям правой части уравнения (3.16) для различных 5
А( ^5 ]p-bp3 + C( 5/5 Ip5-Da-^ = °. (31б)
(
дxl
дp
В области 5 > 5* это уравнение имеет эллиптический тип, которому соответствуют собственные формы в виде пространственно-периодических мод дисторсии на некоторых масштабах Л со слабой анизотропией формы, определяемой величиной усилия. При 5^5* и пересечении сепаратрисы S2 собственные формы уравнения (3.16) претерпевают качественные изменения, связанные с вырождением пространственного периода Л при его стремлении к значениям Л«-ln(5-5*). Решение в этом случае трансформируется к волновому виду бризеров (breather) при 5 = 5* и к автосолитонным волнам p(Q = p(x - Vt) при появлении метастабильности при 5 < 5*. Амплитуда волны, ее скорость и ширина волнового фронта определяются в этом случае условиями неравновесного перехода
1
p = 2 pa
1 -tahn(CLB-1)], Lb = j-(2AJ 2 (3.17)
*
c
Скорость волновых фронтов определяется соотношением V = %A(pa - pm )/Г2, где (pa - pm) - скачок в значении p в области метастабильности. Переход через точку бифуркации 5 c (сепаратриса S3) сопровождается появлением пространственновременных структур качественно нового типа, характеризующихся резкой кинетикой роста дефекта («blow-up» режим) при t ^ tc на спектре некоторых пространственных масштабов [20, 21].
Предполагая степенной вид зависимостей интенсивности источника и коэффициента нелокальности при 5<5C и p > pC, кинетическое уравнение (3.14) может быть записано в форме
* S (pc ) p “ + Jx o(pc ) p^^Xx) ’ ^3'18^
где S, х0, ю, Р - параметры разложения.
В [20] показано, что для данного вида уравнений развитая стадия кинетики p в пределе некоторых характерных времен t ^ tc соответствует автомодельному решению вида
( , \-т
x
p(x, t) = Ф(t )f (С), С = —, Ф(t) = Ф 0 1 “ , (3.19)
L
t
tc у
где т > 0, Ф0 > 0 - параметры, определяемые типом нелинейности уравнения (4.9); Ьс и їс - характерные пространственные и временные масштабы автомодельного решения. Вид функции /(£) определяется решением соответствующей проблемы на собственные формы и собственные значения данной задачи. Например, для случая ш = Р +1 автомодельное решение уравнения (3.14) имеет вид
і
p(x, t) = [S(t - tc)] p ~)w ' v4 sin2 — + лЄ
P(P + 2)
1 , 2(|3 + 1) „ .n2 ^ ^
V^c у у
L
(3.20)
где 9 - случайное значение в интервале (0,1). Масштаб Ьс, так называемая фундаментальная длина [25], имеет смысл пространственного периода решения (3.20)
Ь„ = 2% + 1)Х о X4 )2.
Автомодельное решение (3.20) описывает кинетику режима роста величины дисторсии в так называемом режиме с обострением при I ^ 1с на спектре
пространственных масштабов Ьн = кЬс, к = 1,2,...,К . В этом случае на масштабах Ьн = кЬс развиваются сложные диссипативные структуры обострения, представляющие каскады простых диссипативных структур этого типа (рис. 7).
Рис. 7. Типы гетероклиник и соответствующие собственные формы
Возможные сценарии эволюции локализованных мод дисторсии, зарождающихся первоначально в форме бризеров при некоторых характерных температурах Тт, когда 5^5*, представлены на рис. 7. Появление новых структурных масштабов (соответствующих бризерам) приводит к резкому изменению значений структурного параметра скейлинга 5, переходу системы в область метастабильности и зарождению автосолитонных мод, способных к перемещению вдоль двойной спирали ДНК. Данная способность локализованной моды дисторсии - открытого комплекса -может лежать в основе важнейших функций ДНК - транскрипции и репликации. Механизм транскрипции в этом случае связывается с движением уединенной волны в выделенном направлении, механизм репликации может быть обусловлен движением волновых фронтов в противоположные направления вдоль двойной спирали, обеспечивая разделение ее цепей. Устойчивость данных механизмов связана с особенностями кинетики для структурного параметра скейлинга, играющего роль структурной восприимчивости двойной спирали к росту локализованной дисторсии. Достижение структурным параметром скейлинга критического значения 5С должно
приводить к неограниченному росту локализованных мод дисторсии и нарушению функций ДНК.
Описанный на основе уравнений (3.14) сценарий эволюции коллективных мод дисторсии в структуре двойной спирали ДНК носит характер «самоорганизованной критичности» в зависимости от значений текущей структурной восприимчивости -структурного параметра скейлинга и позволяет проведение анализа поведения биологической системы на основе спектрального анализа типов коллективных мод, контролирующих поведение последней.
Учитывая роль, которую может играть структурный параметр скейлинга в механизмах, происходящих в двойной спирали ДНК, представляется важным более
детальное исследование влияния данного параметра на статистическое поведение реальных биологических объектов, характеризующихся различной «восприимчивостью» к зарождению и эволюции коллективных мод дисторсии.
4. Динамика ДНК и суперстатистика систем с «медленной динамикой»
Динамические процессы, происходящие в структуре двойной спирали ДНК, учитывая существенную неравновесность последних, можно отнести к процессам с так называемой «медленной динамикой» [26]. Идеи, впервые высказанные Леонтовичем об обобщении статистики Больцмана-Гиббса на случай систем, далеких от равновесия и обнаруживающих черты медленной динамики, получили развитие в последние годы в [27, 28] и известны как неэкстенсивная статистическая механика и термодинамика (nonextensive statistical mechanics and thermodynamics), а также как суперстатистика (superstatistics). Как оказалось, областью приложений этих подходов является описание широкого класса сложных систем, поведение и свойства которых характеризуются как аномальные с точки зрения термодинамики и статистики равновесных или близких к равновесию систем. Основной проблемой при изучении таких сложных систем, к которым относятся и биологические системы, является описание связи между статистическими и динамическими (флуктуационными) закономерностями поведения в ситуациях, когда предположение об эргодичности систем в общем случае не выполняется при развитии неустойчивостей, приводящих к стохастической динамике с высокой чувствительностью макроскопического поведения систем к флуктуациям и, как следствие, развитию пространственно-локализованных коллективных мод с большими временами динамики. Сложная пространственновременная структура таких систем приводит к тому, что нарушается принцип аддитивности для такой важнейшей термодинамической величины, как энтропия.
Обобщение статистики Больцмана-Гиббса для систем с медленной динамикой получило название Tsallis статистики (T-statistics) [27]. Формализм T-statistics аналогичен статистике Больцмана-Гиббса за исключением определения энтропии и специальных преобразований Лежандра, необходимых для формулировки соответствующей термодинамики. На уровне описания динамики систем с «медленной динамикой» предположение об эргодичности заменяется неэргодичностью T-statistics. Отметим, что теория Больцмана-Гиббса охватывает системы с экспоненциальным видом вероятности распределений интенсивных величин, T-statistics описывает степенные законы распределений, которые наблюдаются для большинства неравновесных состояний в системах с медленной динамикой.
Формализм T-statistics, а также предложенной в [28] суперстатистики, можно рассматривать как попытку суперпозиции двух различных статистик, относящихся к неравновесным системам под действием медленно меняющихся или постоянных внешних сил и обнаруживающих макроскопически состояния, близкие к стационарным или с медленной динамикой, но сопровождающихся флуктуациями интенсивных параметров. Таким образом, формализм суперстатистики можно определить как статистику статистик, суть которой заключается в следующем. Рассматриваются системы с медленной динамикой как совокупность доменов, статистика величин для которых может быть охарактеризована известным больцмановским фактором ехр(-р£), где E - эффективная энергия данного домена, Р - интенсивный флуктуирующий параметр - обратная эффективная температура, значения которого принимаются постоянными по результатам усреднения для каждого домена на характерном масштабе времени. Предполагается существование независимой
статистики для интенсивного параметра Р и вводится в рассмотрение усредненный больцмановский фактор
ТО
В(Е )=| / (Р)ехр(-РЕ)ЙР, (4.1)
0
где /(р) - функция распределения для Р. Условие нормировки эффективного больцмановского фактора для стационарной (длинновременной) функции распределения имеет вид:
р(Е) = |в(Е), (4.2)
ТО
где 7 = |В(Е)оЕ .
о
Частным случаем суперстатистики является вариант с определением процедуры нормализации в виде
p(E )=| f (Р) щ e-Edp.
(4.3)
где 7(р) - нормализующий коэффициент для е рЕ при заданном р .
Сопоставляя результаты развитого статистического описания и динамики коллективных мод дисторсии в структуре ДНК, представляется естественным обобщение соотношения (4.3) на случай конкретных типов распределений начальных значений структурного параметра скейлинга. Данное обобщение может быть предложено на основе введения функции распределения, имеющей вид
где f (5) - функция распределения начальных значений восприимчивости; E = Eg .
Предложенное обобщение дает возможность более реалистического описания возможных сценариев развития локализованных мод дисторсии в структуре ДНК, для которых не исключается присутствие всех обсуждаемых выше коллективных мод в представительном объеме, содержащих ансамбли ДНК. Эволюция данных коллективных мод и их взаимодействие формируют в течение заданного интервала времени «флуктуационный рельеф», характеризующий состояние системы в присутствии всех коллективных мод, являющихся степенями свободы системы, по отношению к которым может быть определено понятие «эффективных температур» неравновесной системы.
Суперстатистика и классы универсальности статистического поведения неравновесных систем с медленной динамикой
Стационарные функции распределения систем, рассматриваемых в суперстатистике (в дальнейшем S-систем), обычно обнаруживают негауссово поведение с широкими «хвостами» (fat tails), которые могут затухать по степенному закону, или как «растянутые» (stretched) распределения, или же более сложные. Важный признак S-моделей - существование интенсивного параметра Р,
флуктуирующего на больших пространственно-временных масштабах. В [29] показано, что многочисленные экспериментальные данные достаточно хорошо описываются тремя универсальными типами Б-распределений: х2, обратное х2 и Log-нормальное распределения. Указанные типы Б-распределений, по-видимому, представляют универсальный статистический предел для широкого класса динамических систем. Рассмотрим эту особенность данных распределений более детально применительно к обсуждаемым классам неравновесных систем с медленной динамикой.
Данные неравновесные неоднородные системы содержат много пространственных доменов (ячеек), характеризующихся различными значениями интенсивного параметра р (обратной температурой). Размер домена определяется его корреляционной длиной в непрерывно изменяющемся поле р. Предполагается, что в каждой ячейке достигается состояние локального равновесия, то есть времена соответствующих релаксационных процессов предполагаются малыми. Параметр р предполагается постоянным в каждой ячейке на временном масштабе Т, а затем он принимает другие значения. Для больших времен т > Т стационарные распределения таких неоднородных систем возникают как суперпозиция распределений Больцмана ех Р (-р Е ), «взвешенного» с плотностью вероятности /(р) наблюдаемой переменной р в случайно выбранной ячейке
ТО 1
Р(Е )=/ / (р) 7(р)р(Е УрЕс!р. (4.4)
Здесь Е - эффективная энергия каждой ячейки; р(Е) - плотность состояний; 7 (р) -нормирующий коэффициент для р(Е)е_рЕ для данного р. Строго говоря, распределение / (р) определяется пространственно-временной динамикой распределенной системы. Однако ряд содержательных выводов может быть получен при следующем рассмотрении. Так как по смыслу р - величина положительная,
распределение /(р) не может быть гауссовым. Типичные распределения /(р) для распространенных экспериментальных ситуаций следующие:
а) В случае множественных независимых микроскопических случайных переменных
5, 1 = 1,..., I, определяющих р в «аддитивном» смысле, нормированная сумма
1 1
* = 4Г §51
приближается в соответствии с центральной предельной теоремой к гауссовому распределению Х1 для больших значений I. Так как в общем случае могут существовать много различных случайных переменных, содержащих микроскопические случайные переменные, то имеет место совокупное гауссово распределение переменных Х1,...,Хп, соответствующее макроскопическим степеням свободы системы. Так как р является величиной положительной, то ее значение получается возведением в квадрат этих случайных по Гауссу переменных. В результате
р = §х,2
1=1
дается х2 - распределением п степеней свободы
, _«Р/ -1e /2ро
f (Р)=
где Р0 есть среднее значение Р . В [28, 29] показано, что данный вариант ^-статистики соответствует Т^айзИез и обнаруживает степенные «хвосты» для больших значений переменных. Таким образом, Т^айзИез возникает как универсальный динамический предел, то есть когда детали поведения микроскопических случайных переменных Е, 1 (а также их распределений вероятности) несущественны.
б) Аналогичное вышеприведенному рассмотрение может быть проведено, если «температура» Р-1 более, чем сама величина Р, представляет сумму множества квадратов гауссовых случайных переменных, возникающих в ансамбле микроскопических степеней свободы Е,1 . В результате /(р) равно обратному х2_ распределению вида
/ (р)={^р°'] % р-^2-2 е "”%.
Последнее производит ^-распределение, которое характеризуется экспоненциальным затуханием микроскопических переменных и обладает признаками универсальности, так как детали поведения Е, 1 также не существенны.
в) В отличие от ситуации, когда величина Р дается суммой многочисленных вкладов, для других систем, в частности для турбулентности, случайная переменная Р может быть генерирована мультипликативными случайными процессами. В этом случае мы
I
имеем дело с локальными каскадами случайных переменных Х1 =^^г-, где I - число
1=1
«участников» каскада и Е, 1 - положительные микроскопические случайные
переменные. В соответствии с центральной предельной теоремой для больших I случайные переменные
1 1 1 4ГХ1 = 47 £108
имеют гауссово распределение. Как следствие, X1 имеет Log-нормальное распределение. В общем случае в Р имеют вклад n произведений переменных X
р = Пх.
i=1
n
и, как следствие, log Р = ^ Xi представляет собой сумму гауссовых случайных
i=1
распределений. Таким образом, Р имеет вид Log-нормального распределения
где ц из2 - соответствующие величины среднего значения и дисторсии. Очевидно, что этот Ьо§-нормальный результат для Р также независим от деталей случайного каскада микроскопических переменных £, г- и, таким образом, Ьо§-нормальная
статистика является также универсальной.
Несмотря на то, что могут быть рассмотрены и более сложные случаи, в [29] утверждается, что большинство экспериментальных данных соответствуют трем рассмотренным типам универсальности или простым их комбинациям. В [30] отмечается, что формализм суперстатистики связан с необходимостью определения функций распределения для интенсивного параметра Р , что является следствием недостаточного понимания поведения неравновесных систем с «медленной динамикой». Поскольку системы с дефектами обнаруживают яркие черты «медленной динамики» вследствие зарождения коллективных мод различного типа, то естественно провести сопоставление сценариев эволюции динамических систем этого класса с обсуждаемыми выше типами универсальных статистических распределений.
Периодические возмущения при значениях интенсивного параметра 5 >5*, переходящие в бризерные моды при 5^-5*, инициируются термическими флуктуациями. Тип собственных форм, определяемый нелинейностью неравновесного потенциала и эффектами нелокальности, формирует в системе пространственно-распределенные возмущения, соответствующие признакам аддитивности случая а) и X2 -статистики, когда в приближении больших времен микроскопические переменные подчиняются коллективным модам установленного для 5 > 5* типа.
Переход системы в интервал «эффективных температур» 5* > 5 > 5с, связанный с увеличением энергии дефекта при появлении зародыша с размерами, соответствующими масштабам бризера при 5 *, приводит, как это показано выше, к формированию локализованных автосолитонных конечно-амплитудных мод. При этом наблюдается вырождение пространственных и временных масштабов возмущений, существовавших в интервале 5 > 5*. Данный переход сопровождается резким изменением симметрии неравновесной системы в смысле числа значимых динамических координат системы, поведение которой, однако, также контролируется системой макроскопических переменных, представляющих собой коллективные (автосолитонные) моды. По отношению к данной системе автомодельных переменных, которые допускают аддитивность вышерассмотренного типа, естественной является аналогия с обратной х2 -статистикой.
Зарождение каскадов диссипативных обостряющихся структур в интервале 5 < 5с естественно связать с мультипликативностью случайного процесса при подчинении динамики системы спектру коллективных мод (автомодельных решений) данного типа и Ьо§-нормальной статистикой, характеризующей поведение неравновесной системы.
Установление соответствия динамических свойств биологических систем со статистическими закономерностями, которые могут быть исследованы экспериментально, представляется важным для возможных приложений.
Обсуждение результатов
В литературе обсуждается, что важным для понимания динамики ДНК представляется наблюдение локализованных мод, аналогичных тем, которые имеют место в кристаллических решетках и которые возникают при термической активации. Имея в виду эту аналогию и поведение дефектов дислокационного типа в реальных кристаллических решетках [21], можно высказать предположение, что локализованные моды в ДНК существуют, как и бризеры, механизм зарождения и динамика которых объясняется в рамках Peyrard-Bishop модели.
В реальных ДНК, содержащих генетический код, неоднородность и нелинейность является естественным источником для неоднородности флуктуаций и их развития в коллективные моды различной природы, с которыми могут быть связаны уникальные свойства ДНК. Значимость данных эффектов в настоящее время окончательно не оценена для биологии, однако ряд экспериментов определенно указывают на то, что флуктуационно-индуцированные эффекты локального раскрытия в структуре двойной спирали ДНК имеют место.
Другим важным направлением исследований, которое может ярко проявиться в ближайшее время, является изучение свойств и применение ДНК для чисто физических приложений, имея в виду замечательное свойство ДНК демонстрировать высокоселективную самосборку (self-assembly). Основанием для такого оптимизма служит уровень экспериментальных исследований, реализуемых с единичными ДНК, и позволяющий контроль за их конформационной способностью. Понимание закономерностей денатурации и самосборки ДНК может служить основой для создания наномеханизмов, а также для построения объектов с нанометрической точностью.
Благодарности
Автор благодарит Michel Peyrard и Alain Arneodo за содержательные и стимулирующие дискуссии по обсуждаемой проблеме.
Список литературы
1. Schredinger, E. What is Life? / E. Schredinger. - Cambridge: Cambridge University Press, 1944.
2. Watson, J.D. Molecular structure of nucleic acids / J.D. Watson and F.Y.C.Crick // Nature. - 1953. - Vol. 171. - P. 737-745.
3. Calladine, C.R. Understanding DNA / C.R. Calladine and H.R. Drew. - London: Academic, 1992.
4. Peyrard, M. Nonlinear dynamics and statistical physics of DNA / M. Peyrard // Nonlinearity. - 2004. - Vol.
17. - P. 1-40.
5. Frank-Kamenetskii, M. DNA chemistry. How the double helix breathes / M. Frank-Kamenetskii // Nature. -1987. - Vol. 328. - P. 17-25.
6. Wartell, R.M. Thermal denaturation of DNA molecules: a comparison of theory and experiments / R.M. Wartell, A.S. Benight // Phys. Rev. - 1985. - Vol. 126. - P. 67-98.
7. Peyrard, M. Statistical mechanics of a nonlinear model for DNA denaturation / M. Peyrard, A.R. Bishop // Phys. Rev. Lett. - 1989. - Vol. 62. - P. 2755-2758.
8. Dauxois, T. Entropy-driven transition in one dimension system / T. Dauxois and M. Peyrard // Phys. Rev. E.
- 1995. - Vol. 51. - P. 4027-4040.
9. Peyrard, M. Using DNA to probe nonlinear localized excitations / M. Peyrard // Europhys. Lett. - 1998. -Vol. 44. - P. 271-277.
10. Peyrard, M. The pathway to energy localization in nonlinear lattices / M. Peyrard // Physica D. - 1998. -Vol. 119. - P. 184-199.
11. Cule, D. Denaturation of heterogeneous DNA / D. Cule, T. Hwa // Phys. Rev. E. - 1977. - Vol. 79. - P. 2375-2378.
12. Marko, J.F. Statistical mechanics of supercoiled DNA / J.F. Marko, E.D. Siggia // Phys. Rev. E. - 1995. -Vol. 52. - P. 2912-2937.
13. Saenger, W. Principles of Nucleic Acid Structure / W. Saenger. - New York: Springer-Verlag, 1984.
14. Barbi, M. A twist opening model for DNA / M. Barbi, S. Cocco, M. Peyrard, et al // J. of Biological Physics.
- 1999. - Vol. 24. - P. 358-369.
15. Cocco, S. Statistical mechanics of torque induced denaturation of DNA / S. Cocco, R. Monasson // Phys. Rev. Lett. - 1999. - Vol. 83. - P. 5178-5181.
16. Austin, R.H. Stretch genes / R.H. Austin, J.R. Brody, E.C. Cox, et al // Physics Today. - 1997. - Vol. 50. -P. - 32-38.
17. Strick, T.R. Phase coexistence in single DNA molecule / T.R. Strick, J.F. Allemand, D. Benson, et al // Physica A. - 1999. - Vol. 263. - P. 392-404.
18. Wartel, R.M. Thermal denaturation of DNA molecules: a comparison of theory with experiments / R.M. Wartel, A.S. Baldwin // Phys. Rep. - 1985. - Vol. 126. - P. 67.
19. Наймарк, О.Б. Нелинейная динамика обратимой пластичности в сплавах с эффектом памяти формы / О.Б. Наймарк, Л.В. Филимонова, В. А. Баранников и др. // Физическая мезомеханика. - 2001. - Т. 4, № 5. - С. 13-27.
20. Наймарк, О.Б. Неустойчивости в конденсированных средах, обусловленные дефектами / О.Б. Наймарк // Письма в ЖЭТФ. - 1998. - T. 67, № 9. - C. 751-757.
21. Naimark, O.B. Defect Induced Transitions as Mechanisms of Plasticity and Failure in Multifield Continua / O.B. Naimark // Advances in Multifield Theories of Continua with Substructure / Ed. G. Capriz and P. Mariano, Birkhauser, Boston. - 2004. - P. 75-114.
22. Наймарк, О.Б. Коллективные свойства ансамблей дефектов и некоторые нелинейные проблемы пластичности и разрушения / О.Б. Наймарк // Физическая мезомеханика. - 2003. - Т. 6, № 4. -С. 45-72.
23. Леонтович, М.А. О свободной энергии неравновесного состояния / М.А. Леонтович // ЖЭТФ. - 1938.
- Т. 8, № 7. - С. 844-854.
24. Raikher, Yu.L. The effective field method in the orientational kinetics of magnetic fluids and liquid crystals / Yu.L. Raikher, M.I. Shliomis // Relaxation Phenomena in Condensed Matter, W.Coffey (Ed.), Advances in Chemical Physics Series, Vol. LXXXVI, John Willey&Sons. - 1994. - P. 595-751.
25. Kurdjumov, S.P. Evolution and self-organization laws of complex systems / S.P. Kurdjumov // International Journal of Modern Physics. - 1988. - Vol. 1, No. 4. - P. 299-327.
26. Cugliandolo, L.F. Energy flow, partial equilibration and effective temperatures in systems with slow dynamics / L.F. Cugliandolo, J. Kurchan, L. Peliti // Phys. Rev. E. - 1997. - Vol. 55, No. 4. - P. 3898-3914.
27. Tsallis, C. Possible generalization of Boltzmann-Gibbs statistics / C. Tsallis // J. Stat. Phys. - 1988. - Vol. 52. - P. 479-487.
28. Beck, C. Superstatistics / C. Beck, E.C.D. Cohen // Physica A. - 2003. - Vol. 322. - P. 267-275.
29. Beck, C. From time series to superstatistics / C. Beck, E.C.D. Cohen, H.I. Swinney // Phys. Rev. E. - Vol. 72. - P. 56133.
30. Tsallis, C. What should a statistical mechanics satisfy to reflect nature / C. Tsallis // Physica D. - 2004. -Vol. 193. - P. 3-34.
THE LOCALIZED DISTORTION MODES OF THE DNA DOUBLE HELIX STRUCTURE
O.B. Naimark (Perm, Russia)
Up to the present, real mechanisms of some biological process including deoxyribonucleic acid (DNA) are unknown. Because of this, the biologists concentrate detailed attention to structural aspects to elucidate real mechanisms of replication (processes repetition), transcription (the first stage of genetic cellular information realization which results in biosynthesis of ribonucleic acid (RNA) molecules on the matrix of DNA), and denaturation (process of protein molecules decomposition). In the present study, the attempt is made to suggest the equation of motion as the special form of Ginsburg-Landau equation for
the order parameters, and to investigate the dynamic properties of the system under consideration. The equation is based on the solution of statistical problem of the defects ensemble behavior and describes the localized distortions in the DNA structure. It is shown that the dynamics of localized distortion modes is connected with the new type of critical phenomena - structural-scaling transitions, which are typical for the nonequilibrium systems with “slow dynamics”. The form of the new order parameter - the structural scaling parameter
- is established, which plays the role of “effective temperature” for the nonequilibrium system. The analysis of the automodel equations reveals the existence of three types of collective modes in the defects ensembles in the corresponding intervals of structural scaling parameter. The correspondence between the typical collective modes and three types of universal statistical distributions is shown in the superstatistics of nonequilibrium systems with slow dynamics. The presumable connection between collective modes dynamics and transcription and replication mechanisms is discussed.
Key words: deoxyribonucleic acid (DNA), double helix, structure, defects ensemble, statistical analysis, structural scaling parameter, scaling, nonequilibrium systems.
Получено 24 апреля 2006
