Г.Ю.Ризниченко
Биология математическая

Здесь S - количество сеголеток (штуки) на каждого нерестившегося производителя; x - количество зашедших весной из Черного моря в Азовское производителей хамсы (млрд штук); - среднеквадратичное отклонение.

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

Здесь Y поглощение кислорода, измеренное в мкл(0,25 г)-1ч-1.; D - число дней, в течение которых выдерживались образцы; B - процентное содержание влаги в образцах; Т - температура, измеренная в градусах С.

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

(Из книги: Д.Джефферс "Введение в системный анализ: применение в экологии", М., 1981)

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

Здесь i - номер точки измерения; xe - 'экспериментальные значения переменных; хt - теоретические значения переменных; a1, a2... - параметры, подлежащие оценке; wi - "вес" i-го измерения; N - число точек измерения.

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

Имитационные модели созданы для описания физиологических процессов. происходящих в жизненно важных органах: нервном волокне, сердце, мозге, желудочно-кишечном тракте, кровеносном русле. На них проигрываются "сценарии" процессов, протекающих в норме и при различных патологиях, исследуется влияние на процессы различных внешних воздействий, в том числе лекарственных препаратов. Имитационные модели широко используются для описания продукционного процесса растений и применяются для разработки оптимального режима выращивания растений с целью получения максимального урожая, или получения наиболее равномерно распределенного во времени созревания плодов. Особенно важны такие разработки для дорогостоящего и энергоемкого тепличного хозяйства.

Качественные (базовые) модели

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

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

Неограниченный рост. Экспоненциальный рост. Автокатализ.

В основе любых моделей лежат некоторые предположения. Модель, построенная на основе этих предположений, становится самостоятельным математическим объектом, который можно изучать с помощью арсенала математических методов. Ценность модели определяется тем, насколько характеристики модели соответствуют свойствам моделируемого объекта. Одно из фундаментальных предположений, лежащих в основе всех моделей роста - пропорциональность скорости роста численности популяции, будь то популяция зайцев или популяция клеток. В основе этого предположения лежит тот общеизвестный факт, что важнейшей характеристикой живых систем является их способность к размножению. Для многих одноклеточных организмов или клеток, входящих в состав клеточных тканей - это просто деление, то есть удвоение числа клеток через определенный интервал времени, называемый характерным временем деления. Для сложно организованных растений и животных размножение происходит по более сложному закону, но в простейшей модели можно предположить, что скорость размножения вида пропорциональна численности этого вида.

Математически это записывается с помощью дифференциального уравнения, линейного относительно переменной x, характеризующей численность (концентрацию) особей в популяции:

(1)

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

Предположение о пропорциональности скорости роста популяции ее численности было высказано еще в 18 веке Томасом Робертом Мальтусом (1766-1834) в книге "О росте народонаселения" (1798). Согласно закону (1), если коэффициент пропорциональности R=r=сonst (как это предполагал Мальтус), численность будет расти неограниченно по экспоненте.

В своих работах Мальтус обсуждает последствия этого закона в свете того обстоятельства, что производство продовольствия и других товаров растет линейно, и следовательно, популяция, растущая экспоненциально, обречена на голод.

Для большинства популяций существуют ограничивающие факторы, и по тем или иным причинам рост популяции прекращается. Единственное исключение представляет человеческая популяция, которая на протяжении всего исторического времени растет даже быстрее, чем по экспопненте. (См. Экология математическая, раздел Рост численности человечества). Исследования Мальтуса оказали большое влияние как на экономистов, так и на биологов. В частности, Чарльз Дарвин пишет в своих дневниках, что положенные в основу модели Мальтуса предположения и пропорциональности скорости роста популяции ее численности представляются весьма убедительными, и из этого следует неограниченный экспоненциальный рост численности. В то же время, ни одна из популяций в природе не растет до бесконечности. Следовательно, существуют причины, препятствующие такому росту. Одну из таких причин Дарвин видит в борьбе видов за существование.

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

Ограниченный рост. Уравнение Ферхюльста

Базовой моделью, описывающей ограниченный рост, является модель Ферхюльста (1848):

(2)

Здесь параметр K носит название "емкости популяции" и выражается в единицах численности (или концентрации). Он не имеет какого-либо простого физического или биологического смысла и носит, системный характер, то есть определяется целым рядом различных обстоятельств, среди них ограничения на количество субстрата для микроорганизмов, доступного объема для популяции клеток ткани, пищевой базы или убежищ для высших животных.

График зависимости правой части уравнения (2) от численности x и численности популяции от времени представлены на рис. 1 (а и б).

Рис.1. Ограниченный рост. Зависимость величины скорости роста от численности (а) и численности от времени (б) для логистического уравнения.

В последние десятилетия уравнение Ферхюльста переживает вторую молодость. Изучение дискретного аналога уравнения (2) выявило совершенно новые и замечательные его свойства [3,5]. Рассмотрим численность популяции в последовательные моменты времени. Это соответствует реальной процедуре пересчета особей (или клеток) в популяции. В самом простом виде зависимость численности на временном шаге номер n+1 от численности предыдущем шаге n можно записать в виде:

(3)

Поведение во времени переменной x_n может носить характер не только ограниченного роста, как было для непрерывной модели (2), но также быть колебательным или квазистохастическим (рис.2).

Рис.2. Вид функции зависимости численности на последующем шаге от численности на предыдущем шаге (а) и поведение численности во времени (б) при разных значениях параметра r: 1 - ограниченный рост; 2 - колебания, 3 - хаос.

Тип поведения зависит от величины константы собственной скорости роста r. Кривые, представляющие вид зависимости значения численности в данный момент времени (t+1) от значений численности в предыдущий момент времени t представлены на рис. 2 слева. Справа представлены кривые динамики численности - зависимости числа особей в популяции от времени. Сверху вниз значение параметра собственной скорости роста r увеличивается.

Характер динамики численности определяется видом кривой зависимости F(t+1) от F(t). Эта кривая отражает изменение скорости прироста численности от самой численности. Для всех представленных на рис. 2 слева кривых эта скорость нарастает при малых численностях, и убывает, а затем обращается в нуль при больших численностях. Динамический тип кривой роста популяции зависит от того, насколько быстро происходит рост при малых численностях, т.е. определяется производной (тангенсом угла наклона этой кривой) в нуле, который определяется коэффициентом r - величиной собственной скорости роста. Для небольших r (r<3) численность популяции стремится к устойчивому равновесию. Когда график слева становится более крутым, устойчивое равновесие переходит в устойчивые циклы. По мере увеличения численности длина цикла растет, и значения численности повторяются через 2, 4, 8,..., 2ⁿпоколений. При величине параметра r>5,370 происходит хаотизация решений. При достаточно больших r динамика численности демонстрирует хаотические всплески (вспышки численности насекомых).

Уравнения такого типа неплохо описывают динамику численности сезонно размножающихся насекомых с неперекрывающимися поколениями. При этом некоторые достаточно просто измеряемые характеристики популяций, демонстрирующих квазистохастическое поведение, имеют регулярный характер. В некотором смысле, чем хаотичнее поведение популяции, тем оно предсказуемее. Например, при больших x амплитуда вспышки может быть прямо пропорциональна времени между вспышками.

Дискретное описание оказалось продуктивным для систем самой различной природы. Аппарат представления динамического поведения системы на плоскости в координатах [x_t, x_t+T] позволяет определить, является наблюдаемая система колебательной или квазистохастической. Например, представление данных электрокардиограммы позволило установить, что сокращения человеческого сердца в норме носят нерегулярный характер, а в период приступов стенокардии или в прединфарктном состоянии ритм сокращения сердца становится строго регулярным. Такое "ужесточение" режима является защитной реакцией организма в стрессовой ситуации и свидетельствует об угрозе жизни системы.

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

Еще более интересные математические объекты получаются, если переписать уравнение (3) в виде:

(4)

Рис.3. Изображение множества Мандельброта на комплексной плоскости для процесса xx²+c. Ось абсцисс - Re c, ось ординат - Im c. Изображенная фигура имеет фрактальную структуру, характерной чертой которой является самоподобие.

и рассматривать константу с в комплексной области. При этом получаются объекты, называемые множествами Мандельброта (пример приведен на рис. 3). Подробнее об этих множествах можно прочитать в книге "Красота фракталов" (Образы комплексных динамических систем)[2], там же приведены их многочисленные красочные изображения. Имеют ли эти объекты биологическую интерпретацию, имеющую под собой глубокий смысл, или это просто красивый "сюрприз", который нам преподносит базовая система? Пока на этот вопрос нет окончательного ответа.

Ограничения по субстрату. Модели Моно и Михаэлиса-Ментен.

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

(5)

Здесь К_S - константа, равная концентрации субстрата, при которой скорость роста равна половине максимальной. ₀ - максимальная скорость роста, равная величине r в формуле (2). Это уравнение было впервые написано крупнейшим французким биохимиком. Жаком Моно (1912-1976). Совместно с Франсуа Жакобом им были разработаны представления о роли транспортной рибонуклеиновой кислоты (mRNA) в аппарате размножения клетки. В развитие представлений о генных комплексах, которые были ими названы оперонами, Жакоб и Моно постулировали существование класса генов, которые регулируют функционирование других генов путем воздействия на синтез транспортной РНК. Такой механизм генной регуляции впоследствии полностью подтвердился для бактерий, за что обоим ученым (а также Андре Львову) была присуждена Нобелевская премия 1965 г. Ниже рассмотрена знаменитая модель генной регуляции синтеза двух ферментов, названная триггерной моделью Жакоба и Моно.

Жак Моно был также философом науки и незаурядным писателем. В своей знаменитой книге "Случайность и необходимость", 1971 Моно высказывает мысли о случайности возникновения жизни и эволюции, а также о роли человека и его ответственности за происходящие на Земле процессы.

Любопытно, что модель Моно (5) по форме совпадает с уравнением Михаэлиса-Ментен (1913), которое описывает зависимость скорости ферментативной реакции от концентрации субстрата при условии, когда общее количество молекул фермента постоянно и значительно меньше количества молекул субстрата:

(6)

Здесь К_М - константа Михаэлиса, одна из важнейших для ферментативных реакций величина, определяемая экспериментально, имеющая смысл и размерность концентрации субстрата, при которой скорость реакции равна половине максимальной.

Закон Михаэлиса-Ментен выводится на основании уравнений химической кинетики и описывает скорость образования продукта в соответствии со схемой:

E+S [ES] E+P.

Сходство уравнений (5) и (6) не случайно. Формула Михаэлиса-Ментен (5) отражает более глубокие закономерности кинетики ферментативных реакций, которые в свою очередь определяют жизнедеятельность и рост микроорганизмов, описываемые эмпирической формулой (5).

Базовая модель взаимодействия. Конкуренция. Отбор.

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

В случае, когда взаимодействие происходит случайно, его интенсивность определяется концентрацией взаимодействующих компонентов и их подвижностью обобщенной диффузией. Именно такие представления приняты в базовых моделях взаимодействия видов. Классической книгой, в которой рассматриваются математические модели взаимодействия видов стала книга Вито Вольтерра "Математическая теория борьбы за существование" (1931) [1]. Книга, построена как математический трактат, в ней постулированы в математической форме свойства биологических объектов и их взаимодействий, а затем эти взаимодействия исследуются как математические объекты, Именно с этой работы В.Вольтерра начались современная математическая биология и математическая экология.

Вито Вольтерра (1860-1940) завоевал мировую известность своими работами в области интегральных уравнений и функционального анализа. Кроме чистой математики его всегда интересовали вопросы применения математических методов в биологии, физике, социальных науках. В годы службы в ВВС в Италии, он много работал над вопросами военной техники и технологии (задачи баллистики, бомбометания, эхолокации). В этом человеке сочетался талант ученого и темперамент активного политика, принципиального противника фашизма. Он был единственным итальянским сенатором, проголосовавшим против передачи власти Муссолини. Когда в годы фашистской диктатуры в Италии Вольтерра работал во Франции, Муссолини, желая привлечь на свою сторону всемирно известного ученого, предлагал ему различные высокие посты в фашистской Италии, но всегда получал решительный отказ. Антифашистская позиция привела Вольтерра к отказу от кафедры в Римском университете и от членства в итальянских научных обществах.

Серьезно вопросами динамики популяций В.Вольтерра стал интересоваться с 1925 г. после бесед с молодым зоологом Умберто Д'Анкона, будущим мужем его дочери, Луизы. Д'Анкона, изучая статистику рабных рынков на Адриатике, установил любопытный факт: когда в годы первой мировой войны (и сразу вслед за ней) интенсивность промысла резко сократилась, то в улове увеличилась относительная доля хищных рыб. Такой эффект предсказывался моделью "хищник-жертва", предложенной Вольтерра. Эту модель мы рассмотрим ниже. По сути дела это был первый успех математической биологии.

Вольтерра предположил по аналогии со статистической физикой, что интенсивность взаимодействия пропорциональна вероятности встречи (вероятности столкновения молекул), то есть произведению концентраций. Это и некоторые другие предположения (См. Популяционная динамика) позволили построить математическую теорию взаимодействия популяций одного трофического уровня (конкуренция) или разных трофическиъх уровней (хищник-жертва).

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

1. Пища либо имеется в неограниченном количестве, либо ее поступление с течением времени жестко регламентировано.

2. Особи каждого вида отмирают так, что в единицу времени погибает постоянная доля существующих особей.

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

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

5. Если вид питается пищей, имеющейся в неограниченном количестве, прирост численности вида за единицу времени пропорционален численности вида.

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

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

Действительно, согласно гипотезам Вольтерра, скорость процесса отмирания каждого вида пропорциональна численности вида. В химической кинетике это соответствует мономолекулярной реакции распада некоторого вещества, а в математической модели - отрицательным линейным членам в правых частях уравнений. Согласно представлениям химической кинетики, скорость бимолекулярной реакции взаимодействия двух веществ пропорциональна вероятности столкновения этих веществ, т.е. произведению их концентрации. Точно так же, согласно гипотезам Вольтерра, скорость размножения хищников (гибели жертв) пропорциональна вероятности встреч особей хищника и жертвы, т.е. произведению их численностей. И в том и в другом случае в модельной системе появляются билинейные члены в правых частях соответствующих уравнений. Наконец, линейные положительные члены в правых частях уравнений Вольтерра, отвечающие росту популяций в неограниченных условиях, соответствуют автокаталитическим членам химических реакций. Такое сходство уравнений в химических и экологических моделях позволяет применить для математического моделирования кинетики популяций те же методы исследований, что и для систем химических реакций. Можно покащзать, что вольтерровские уравнения могут быть получены не только из локального "принципа встреч", ведущего свое происхождение из статистической физики, но и исходя из баланса масс каждого из компонентов ценоза и энергетических потоков между этими компонентами.

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

Пусть имеется два совершенно одинаковых вида с одинаковой скоростью размножения, которые являются антагонистами, то есть при встрече они угнетают друг друга. Модель их взаимодействия может быть записана в виде:

(7)

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

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

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

Рис.4. Схема синтеза двух ферментов Жакоба и Моно (а) и фазовый портрет триггерной систем (б) [6].

Еще одной классической триггерной системой является модель альтернативного синтеза двух ферментов Жакоба и Моно. Схема синтеза приведена на рис. 4а. Ген-регулятор каждой системы синтезирует неактивный репрессор. Этот репрессор, соединяясь с продуктом противоположной системы синтеза ферментов, образует активный комплекс. Активный комплекс, обратимо реагируя с участком структурного гена опероном, блокирует синтез mРНК. Таким образом, продукт второй системы Р₂ является корепрессором первой системы, а Р₁ - корепрессором второй. При этом в процессе корепрессии могут участвовать одна, две и более молекул продукта. Очевидно, что при таком характере взаимодействий при интенсивной работе первой системы вторая будет заблокирована, и наоборот. Модель такой системы предложена и подробно изучены в школе проф. Д.С.Чернавского [4] После соответствующих упрощений, уравнения, описывающие синтез продуктов Р₁ и Р₂имеют вид:

(8)

Здесь P₁, P₂ - концентрации продуктов, величины A₁, A₂, B₁, B₂, выражаются через параметры своих систем. Показатель степени m показывает, сколько молекул активного репрессора (соединений молекул продукта с молекулами неактивного репрессора, который предполагается в избытке) соединяются с опероном для блокировки синтеза mRNK.

Фазовый портрет системы, (изображение траекторий системы при разных начальных условиях на координатной плоскости, по осям которой отложены величины переменных системы), для m=2 изображен на рис.4б. Он имеет тот же вид, что и фазовый портрет системы двух конкурирующих видов. Сходство свидетельствует о том, что в основе способности системы к переключениям лежит конкуренция - видов, ферментов, состояний.

Рис.5. Модель химических реакций Лотки. Фазовый портрет системы при значениях параметров, соответствующих затухающим колебаниям.

Классические модели Лотки и Вольтерра

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

Уравнение Лотки рассмотрено им в 1926 г. в книге и описывает систему следующих химических реакций

В некотором объеме находится в избытке вещество А. Молекулы А с некоторой постоянной скоростью превращаются в молекулы вещества X (реакция нулевого порядка). Вещество X может превращаться в вещество Y, причем скорость этой реакции тем больше, чем больше конценрация вещества Y - реакция второго порядка. В схеме это отражено обратной стрелкой над символом y. Молекулы Y в свою очередь необратимо распадаются, в результате образуется вещество B (реакция первого порядка).

Запишем систему уравнений, описывающих реакцию:

(9)

Здесь X, Y, B - концентрации химических компонентов. Первые два уравнения этой системы не зависят от B, поэтому их можно рассматривать отдельно. При определенных значениях параметров в системе возможны затухающие колебания.

Базовой моделью незатухающих колебаний служит классическое уравнение Вольтерра, описывающее взаимодействие видов типа хищник-жертва. Как и в моделях конкуренции (8), взаимодействие видов описывается в соответствии с принципами химической кинетики: скорость убыли количества жертв (x) и скорость прибыли количества хищников (y) считается пропорциональными их произведению

(10)

Рис.6. Модель хищник-жертва Вольтерра, описывающая незатухающие колебания численности. А. Фазовый портрет. Б. Зависимость численности жертвы и хищника от времени.

На рис. 6 представлены фазовый портрет системы, по осям которого отложены численности жертв и хищников - (а) и кинетика численности обоих видов - зависимость численности от времени-(б). Видно, что численности хищников и жертв колеблются в противофазе.

Модель Вольтерра имеет один существенный недостаток. Параметры колебаний ее переменных меняются при флуктуациях параметров и переменных системы. Такую систему называют негрубой.

Этот недостаток устранен в более реалистичных моделях. Модификация модели Вольтерра с учетом ограниченности субстрата в форме Моно (уравнение 5) и учет самоограничения численности (как в уравнении 2) приводит к модели, подробно изученной А.Д.Базыкиным в книге "Биофизика взаимодействующих популяций" (1985).

, (11)

Система (11 представляет собой некий кентавр, составленный из базовых уравнений (1, 2, 5, 10)и объединяющий их свойства. Действительно, при малых численностях и в отсутствие хищника жертва (x) будет размножаться по экспоненциальному закону(1). Хищник (y) в отсутствие жертв будут вымирать также по экспоненте. Если особей того или иного вида много, в соответствии с базовой моделью (2) срабатывает системный ферхюльстовский фактор (член -Ex²в первом уравнении, и -My² - во втором). Интенсивность взаимодействия видов считается пропорциональной произведению их численностей (как в модели (10)) и описывается в форме Моно (модель 5). Здесь роль субстрата играет вид-жертва , а роль микроорганизмов - вид-хищник. Таким образом, модель (11) брала в себя свойства базовых моделей (1), (2), (5), (10).

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

Колебания и ритмы в биологических системах

Для биологических систем характерно периодическое изменение различных характеристик. Период этих колебаний может быть связан с периодическими изменениями условий жизни на Земле - смена времен года, смена дня и ночи. Но многие периодические процессы имеют частоту изменения, не связанную очевидным образом с внешними гео-космическими циклами. Это так называемые "биологические часы" различной природы, начиная от колебаний биомакромолекул, биохимических колебаний, вплоть до популяционных волн.

Рис. 7. Биологические и геофизические ритмы в природе.

С некоторыми из типов периодических движений мы уже имели дело при рассмотрении моделей Лотки и Вольтерра. Однако "биологические часы" имеют свойство, отличающее их от рассмотренных типов колебаний - неизменность во времени периода и амплитуды таких колебаний, означающую стационарность и устойчивость колебательного режима.

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

Рис. 8. Схема реакций гликолиза.

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

Классическим примером колебательной биохимической реакции является гликолиз. В процессе гликолиза осуществляется распад глюкозы и других сахаров , при этом соединения, содержащие шесть молекул углерода, превращаются в трикарбоновые кислоты, включающие три молекулы углерода. За счет избытка свободной энергии в процессе гликолиза на одну молекулу шестиуглеродного сахара образуются две молекулы АТФ. Основную роль в генерации наблюдаемых колебаний концентраций компонентов реакции: фруктозо-6-фосфата, фруктозо-1,6-фосфата и восстановленного НАД играет ключевой фермент гликолитического пути - фосфофруктокиназа (ФФК). Полная схема гликолитических реакций изображена на рис.8

Упрощенная схема реакций может быть представлена в виде:

На схеме [Гл] - глюкоза, Ф6Ф - фруктозо-6-фосфат субстрат ключевой реакции, ФДФ - продукт этой реакции, который является субстратом в следующей стадии. Обе реакции катализируются ферментами. В безразмерных переменных система описывающих реакции уравнений может быть записана в виде:

(12)

Здесь зависимости скоростей реакций от переменных записаны в форме Михаэлиса-Ментен (Моно), как это было представлено в уравнении (6). Кинетика изменений переменных и фазовые портреты системы при разных значениях параметров представлены на рис. 9.

Рис. 9. Модель гликолиза. Кинетика изменений концентраций фруктозо-6-фосфата (x) и фруктозодифосфата (y) (слева) и фазовый портрет системы (справа) при разных значениях параметров системы. а. бесколебательный процесс, б - затухающие колеьания, в-квазигармонические колебания, г - релаксационные колебания. (Higgins, 1967)

Интересно, что колебательные рекции в системе гликолиза были сначала предсказаны на математической модели Хиггинсом (1964), и лишь после этого зарегистрированы экспериментально с помощью метода дифференциальной спектрофотометрии в лаборатории Б.Чанса (1966).

Рис. 10. Схема процессов, приводящих к внутриклеточным колебаниям кальция, (Dupont, Goldbeter, 1983). IP₃ - рецептор, стимулирующий колебания.

Во многих типах живых клеток наблюдаются колебания внутриклеточной концентрации кальция, период которых может варьировать от 0,5 до 10 мин. Простейшая схема процессов, приводящих к гармонально обусловленным колебаниям кальция, основой которых служит процессе кальций индуцированного выхода кальция из клетки, приведена на рис.10. Такие колебания впервые наблюдались Эндо с соавторами (1970) на клетках скелетных мышц, Фабиато (1975) на клетках саркоплазматического ретикулума сердца быка, и позднее - многими другими исследователями

Схема процессов, приводящим к колебаниям Ca, предложенная и описанная на модели Дюпоном и Гольдбетером (1989), представлена на рис. 10. Рассматриваются приток и отток кальция в клетку через плазматическую мембрану (константы скоростей v₁ и v ₂, соответственно); гармонально активируемое освобождение кальция из пула (скорость v₃); активный транспорт цитозольного кальция в пул, (v₄), освобождение кальция из пула , активируемое цитозольным кальцием (v₅); свободный отток кальция из пула в цитозоль (v₆). Модель состоит из двух дифференциальных уравнений

(13)

Здесь S₁- концентрация кальция в цитозоле, S₁- концентрация кальция в гармонально чувствительном пуле.

Рис.11. Модель внутриклеточных колебаний кальция. Кинетика концентрации Са при разных значениях параметров. (Dupont, Goldbeter, 1983).

Выражения для величин скоростей были предложены в Simogyi, Stuckin (1991):

Модель предсказывает колебания концентрации кальция во времени, по форме близкие к экспериментальным.(рис. 11).

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

Внутриклеточные колебания задают эндогенные биологические ритмы (биологические часы), которые свойственны всем живым системам [6]. Именно они определяют периодичность деления клеток, отмеряют время рождения и смерти живых организмов. Модели колебательных систем типа (12.,13) используются в ферментативном катализе, теории иммунитета, в теории трансмембранного ионного переноса, микробиологии и биотехнологии. Колебательные модели популяционной динамики нашли успешное применение для разработки методов борьбы с насекомыми-вредителями полей, садов и лесов.

Волны жизни.

До сих пор мы говорили о базовых моделях поведения живых систем во времени. Стремление к росту и размножению ведет к распространению в пространстве, занятию нового ареала, экспансии живых организмов. Жизнь распространяется так же как пламя по степи во время степного пожара. Эта метафора отражает тот факт, что пожар (в одномерном случае - распространение пламени по бикфордову шнуру) описывается с помощью той же базовой модели, что и распространение вида. Знаменитая в теории горения модель ПКП (Петровского - Колмогорова - Пискунова) впервые была предложена ими в 1937 г. именно в биологической постановке как модель распространения доминирующего вида в пространстве. Все три автора этой работы являются крупнейшими российскими математиками. Академик Иван Георгиевич Петровский (1901-1973) - автор фундаментальных трудов по теории дифференциальных уравнений, алгебре, геометрии, математической физике, в течение более 20 лет был ректором Московского Государственного университета им. М.В.Ломоносова. (1951-1973). Андрей Николаевич Колмогоров (1903-- ) глава российской математической школы по теории вероятностей и теории функций, автор фундаментальных трудов по математической логике, топологии, теории дифференциальных уравнений, теории информации, организатор школьного и университетского математического образования, написал несколько работ, в основу которых положены биологические постановки. В частности в 1936 г. он предложил и подробно исследовал обобщенную модель взаимодействия видов типа хищник-жертва (исправленный и дополненный вариант 1972).(См. Популяционная динамика)

Рассмотрим постановку задачи о распространении вида в активной - богатой энергией (пищей) среде. Пусть в любой точке прямой r>0 размножение вида описывается функцией f(x) = x(1-x). В начальный момент времени вся область слева от нуля занята видом x, концентрация которого близка к единице.. Справа от нуля - пустая территория. В момент времени t=0 вид начинает распространяться (диффундировать) вправо с константой диффузии D. Процесc описывается уравнением:

(14)

При t>0 в такой системе начинает распространяться волна концентраций в область r>0, которая является результатом двух процессов: случайного перемещения особей (диффузии частиц) и размножения, описываемого функцией f(x). С течением времени фронт волны перемещается вправо, причем его форма приближается к определенной предельной форме. Скорость перемещения волны определяется коэффициентом диффузии и формой функции f(x), и для функции f(x), равной нулю при x=0 и x=1 и положительной в промежуточных точках, выражается простой формулой: =2Df'(0).

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

Автоволны и диссипативные структуры. Базовая модель "брюсселятор"

На рассмотренной выше одномерной модели (14) видно, что взаимодействие нелинейной химической реакции и диффузии приводит к нетривиальным режимам. Еще более сложного поведения следует ожидать в двумерных моделях, описывающих взаимодействие компонентов системы. Первая такая модель была изучена Тьюрингом в работе под названием "Химические основы морфогенеза". Алан М.Тьюринг (1912-1954) английский математик и логик, прославился своими работами по компьютерной логике и терии автоматов. В 1952 г. он опубликовал первую часть исследования, посвященного математической теории образования структур в первоначально однородной системе, где одновременно проходят химические реакции, в том числе автокаталитические процессы, сопровождаемые потреблением энергии, и пассивные процессы переноса - диффузия. Это исследование осталось незаконченным, так как он покончил жизнь самоубийством, находясь под действием депрессантов, которыми его принудительно лечили в тюрьме, где он отбывал срок по обвинению в гомосексуализме.

Работа Тьюринга стала классической, ее идеи легли в основу современной теории нелинейных систем, теории самоорганизации и синергетики. Рассматривается система уравнений:

(15)

Уравнения такого типа называются уравнениями "реакция-диффузия". В линейных системах диффузия процесс, который приводит к выравниванию концентраций во всем реакционном объеме. Однако в случае нелинейного взаимодействия переменных x и y, в системе может возникать неустойчивость гомогенного стационарного состояния и образуются сложные пространственно-временные режимы типа автоволн или диссипативных структур - стационарных во времени и неоднородных по пространству распределений концентраций, существование которых поддерживается в активных средах за счет потребления энергии системы в процессах диссипации. Условием возникновения структур в таких системах является различие коэффициентов диффузии реагентов, а именно, наличие близкодействующего "активатора" с малым коэффициентом диффузии и дальнодействующего "ингибитора" с большим коэффициентом диффузии.

Такие режимы в двухкомпонентной системе были изучены в деталях на базовой модели под названием "брюсселятор" (Пригожин и Лефевр, 1968), названной в честь брюссельской научной школы под руководством И.Р.Пригожина, в которой наиболее интенсивно проводились эти исследования. Илья Романович Пригожин (род 1917 г. в Москве) - всю жизнь работал в Бельгии. С 1962 г. он - директор Международного Сольвеевского института физической химии в Брюсселе, а с 1967 г. - директор Центра статистической механики и термодинамики Техасского университета (США). В 1977 г. он получил Нобелевскую премию за работы по нелинейной термодинамике, в частности по теории диссипативных структур - устойчивых во времени неоднородных в пространстве структур. Пригожин является автором и соавтором целого ряда книг ["Термодинамическая теория структуры, устойчивости и флуктуаций", "Порядок из хаоса", "Стрела времени", и др.], в которых он развивает математические, физико-химические, биологические и философские идеи теории самоорганизации в нелинейных системах, исследует причины и закономерности рождения "порядка из хаоса" в богатых энергией открытых для потоков вещества и энергии системах, далеких от термодинамического равновесия, под действием случайных флуктуаций.

Классическая модель "брюсселятор" имеет вид

(16)

и описывает гипотетическую схему химических реакций:

Ключевой является стадия превращения двух молекул x и одной молекулы y в x так называемая тримолекулярная реакция. Такая реакция возможна в процессах с участием ферментов с двумя каталитическими центрами. Нелинейность этой реакции в сочетании с процессами диффузии вещества и обеспечивает возможность пространственно-временных режимов, в том числе образование пространственных структур в первоначально однородной системе морфогенез.

Реакция Белоусова-Жаботинского

Пространственно-временные режимы, предсказанные в моделях типа реакция-диффузия можно наблюдать на химических моделях. Самой знаменитой из них является реакция, которую в 1958 г. описал русский химик Белоусов - окисление лимонной кислоты броматом калия, катализируемое ионной парой Се⁴⁺/Се³⁺. Изучение этой реакции было продолжено Жаботинским (1964), который показал, что в качестве катализатора вместо церия можно применять марганец и железо, а в качестве восстановителя вместо лимонной кислоты можно использовать ряд органических соединений, либо имеющих метиленовую группу, либо образующих ее при окислении. К таким соединениям относятся малоновая и броммалоновая кислоты. Обычно реакция проводят при 25⁰ С в сернокислом растворе смеси бромата калия, малоновой и броммалоновой кислот и сернокислого церия. Экспериментальному и модельному изучению реакции Белоусова-Жаботинского посвящены тысячи работ, так как она дает возможность на простой химической системе наблюдать особенности сложных процессов самиоорганизации и допускает различного типа управление, в том числе с помощью различных режимов освещения (С. Мюллер, В. Зыков, 1993)

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

Рис.12. Типы пространственных режимов в реакции Белоусова-Жаботинского. (Жаботинский, 1974)

Рис.13. Схема реакции Белоусова-Жаботинского.

Механизм реакции весьма сложен с химической точки зрения и содержит десятки промежуточных стадий. Упрощенная схема реакции представлена на рис. 13.

Можно выделить несколько основных стадий, которые определяют колебательный характер реакции. Эти: 1) окисление трехвалентного церия броматом:

2) восстановление четырехвалентного церия малоновой кислотой:

Продукты восстановления бромата, образующиеся на стадии 1, бромируют МК. Получающиеся бромпроизводные МК разрушаются с выделение Br^-. Бромид является сильным ингибитором реакции. Процессы описываются моделью, предложенной А.М. Жаботинским и описанной им в книге "Концентрационные автоколебания", М, Наука, 1974:

в модели А.М.Жаботинского x - концентрация ионов церия , y - концентрация автокатализатора; z - концентрация бромида. Учитывая иерархию констант скоростей реакций можно заменитьб дифференциальное уравнение для z алгебраическим и после введения безразмерных переменных прийти к системе двух уравнений, описывающие колебания концентраций реагентов

Модель хорошо описывает релаксацинные колебания Се⁴⁺, представленные на рис. 14.

Рис.14. Колебания в модели Жаботинского. (Жаботинский, 1974)

Период четко делится на две части - Т1- фаза нарастания и Т2 - фаза спада.

В мировой литературе наибольшее распространение получила модель "орегонатор", предложенная Филдом и Нойесом (1974). В совпременных исследованиях в качестве локального элемента наиболее часто используется модель

Здесь малые параметры e, d, отражают соответствующую иерархию времен процессов, x соответствует безразмерной концентрации HBrO₂, y - Br^-, z - Сe ⁴⁺

Для изучения пространственно-временных структур часто используется модель, описывающая пространственно-временную динамику HBrO₂ (переменная u) и катализатора Сe ⁴⁺ (переменная v).

(17)