6) Траектории аттрактора, вычисленные Хао Байлинем для «брюсселятора» при периодическом подводе извне компоненты Х (личное сообщение).
На рис. 8 представлены результаты численных расчетов Хао Байлиня, дающие общее представление об очень сложной структуре такого странного аттрактора для модели, обобщающей «брюсселятор» на случай периодического подвода извне вещества X. Замечательно, что большинство описанных нами типов поведения реально наблюдалось в неорганической химии и в некоторых биологических системах.
Рис. 9. Схема химического реактора, используемого при исследовании колебаний в реакции Белоусова—Жаботинского (однородность реакционной смеси обеспечивает перемешивающее устройство). В реакции участвуют более тридцати продуктов и промежуточных соединений. Эволюция различных путей реакции зависит (помимо других факторов) от концентраций исходных веществ, регулируемых насосами на входе в реактор.
В неорганической химии наиболее известным примером колебательной системы является реакция Белоусова—Жаботинского, открытая в начале 50-х гг. нашего века. Соответствующая схема реакций, получившая название орегонатор, была предложена Нойесом и сотрудниками. По существу, она аналогична «брюсселятору», но отличается большей сложностью. Реакция Белоусова—Жаботинского состоит в окислении органической (малоновой) кислоты броматом калия в присутствии соответствующего катализатора — церия, марганца или ферроина.
В различных экспериментальных условиях у одной и той же системы могут наблюдаться различные формы самоорганизации — химические часы, устойчивая пространственная дифференциация или образование волн химической активности на макроскопических расстояниях[147].
Обратимся теперь к самому интересному вопросу: что дают все эти результаты для понимания функционирования живых систем?
5. Первое знакомство с молекулярной биологией
Ранее в этой главе мы уже показали, что в сильно неравновесных условиях протекают процессы самоорганизации различных типов. Одни из них приводят к установлению химических колебаний, другие — к появлению пространственных структур. Мы видели, что основным условием возникновения явлений самоорганизации является существование каталитических эффектов.
В то время как в неорганическом мире обратная связь между «следствиями» (конечными продуктами) нелинейных реакций и породившими их «причинами» встречается сравнительно редко, в живых системах обратная связь (как установлено молекулярной биологией), напротив, является скорее правилом, чем исключением. Автокатализ (присутствие вещества Х ускоряет процесс образования его в результате реакции), автоингибиция (присутствие вещества Х блокирует катализ, необходимый для производства X) и кросс-катализ (каждое из двух веществ, принадлежащих различным цепям реакций, является катализатором для синтеза другого) лежат в основе классического механизма регуляции, обеспечивающего согласованность метаболической функции.
Нам бы хотелось подчеркнуть одно любопытное различие. В примерах самоорганизации, известных из неорганической химии, молекулы, участвующие в реакциях, просты, тогда как механизмы реакций сложны (например, в реакции Белоусова—Жаботинского удалось установить около тридцати различных промежуточных соединений). С другой стороны, во многих примерах самоорганизации, известных из биологии, схема реакции проста, тогда как молекулы, участвующие в реакции веществ (протеинов нуклеиновых кислот и т. д.), весьма сложны и специфичны. Отмеченное нами различие вряд ли носит случайный характер. В нем проявляется некий первичный элемент, присущий различию между физикой и биологией. У биологических систем есть прошлое. Образующие их молекулы — итог предшествующей эволюции; они были отобраны для участия в автокаталитических механизмах, призванных породить весьма специфические формы процессов организации.
Описание сложной сети метаболической активности и торможения является существенным шагом в понимании функциональной логики биологических систем. К последней мы относим включение в нужный момент синтеза необходимых веществ и блокирование тех химических реакций, неиспользованные продукты которых могли бы угрожать клетке переполнением.
Основной механизм, с помощью которого молекулярная биология объясняет передачу и переработку генетической информации, по существу, является петлей обратной связи, т. е. нелинейным механизмом. Дезоксирибонуклеиновая кислота (ДНК), содержащая в линейно упорядоченном виде всю информацию, необходимую для синтеза различных основных протеинов (без которых невозможно строительство и функционирование клетки), участвует в последовательности реакций, в ходе которых вся информация кодируется в виде определенной последовательности различных протеинов. Некоторые ферменты осуществляют обратную связь среди синтезированных протеинов, активируя и регулируя не только различные стадии превращений, но и автокаталитический механизм репликации ДНК, позволяющий копировать генетическую информацию с такой же скоростью, с какой размножаются клетки.
Молекулярная биология — один из наиболее ярких примеров конвергенции двух наук. Понимание процессов, происходящих на молекулярном уровне в биологических системах, требует взаимно дополняющего развития физики и биологии, первой — в направлении сложного, второй — простого.
Фактически уже сейчас физика имеет дело с исследованием сложных ситуаций, далеких от идеализации, описываемых равновесной термодинамикой, а молекулярная биология добилась больших успехов в установлении связи живых структур с относительно небольшим числом основных биомолекул. Исследуя множество самых различных химических механизмов, молекулярная биология установила мельчайшие детали цепей метаболических реакций, выяснила тонкую, сложную логику регулирования, ингибирования и активации каталитической функции ферментов, связанных с критическими стадиями каждой из метаболических цепей. Тем самым молекулярная биология установила на микроскопическом уровне основы тех неустойчивостей, которые могут происходить в сильно неравновесных условиях.
В некотором смысле живые системы можно сравнить с хорошо налаженным фабричным производством: с одной стороны, они являются вместилищем многочисленных химических превращений, с другой — демонстрируют великолепную пространственно-временную организацию с весьма неравномерным распределением биохимического материала. Ныне перед нами открывается возможность связать воедино функцию и структуру. Рассмотрим кратко два примера, интенсивно исследовавшиеся в последние годы.
Начнем с гликолиза: цепи метаболических реакций, приводящих к расщеплению глюкозы и синтезу аденозинтрифосфата (АТФ) — универсального аккумулятора энергии, общего для всех живых клеток. При расщеплении каждой молекулы глюкозы две молекулы АДФ (аденозиндифосфата) превращаются в две молекулы АТФ. Гликолиз может служить наглядным примером взаимной дополнительности аналитического подхода биологии и физического исследования устойчивости в сильно неравновесной области[148].
В ходе биохимических экспериментов были обнаружены колебания во времени концентраций, связанных с гликолитическим циклом[149]. Было показано, что эти колебания определяются ключевой стадией в цепи реакций — стадией, активируемой АДФ и ингибируемой АТФ. Это — типично нелинейное явление, хорошо приспособленное к регулированию метаболизма. Всякий раз, когда клетка черпает энергию из своих энергетических резервов, она использует фосфатные связи, и АТФ превращается в АДФ. Таким образом, накопление АДФ внутри клетки свидетельствует об интенсивном потреблении энергии и необходимости пополнить энергетические запасы, в то время как накопление АТФ означает, что расщепление глюкозы может происходить в более медленном темпе.
Теоретическое исследование гликолиза показало, что предложенный механизм действительно может порождать концентрационные колебания, т. е. обеспечивать работу химических часов. Вычисленные из теоретических соображений значения концентраций, необходимые для возникновения колебаний, и величина периода цикла согласуются с экспериментальными данными. Гликолитические колебания вызывают модуляцию всех энергетических процессов в клетке, зависящих от концентрации АТФ, и, следовательно, косвенно влияют на другие метаболические цепи.
Можно пойти еще дальше и показать, что в гликолитическом цикле ход реакций регулируется некоторыми ключевыми ферментами, причем сами реакции протекают в сильно неравновесных условиях. Такие расчеты были выполнены Бенно Хессом[150], а полученные результаты обобщены и на другие системы. При обычных условиях; гликолитический цикл соответствует химическим часам, но изменение этих условий может привести к образованию пространственных структур в полном соответствии с предсказаниями на основе существующих теоретических моделей.
С точки зрения термодинамики живая система отличается необычайной сложностью. Одни реакции протекают в слабо неравновесных условиях, другие — в сильно неравновесных условиях. Не все в живой системе «живо». Проходящий через живую систему поток энергии несколько напоминает течение реки — то спокойной и плавной, то низвергающейся водопадом и высвобождающей часть накопленной в ней энергии.
Рассмотрим еще один биологический процесс, также исследованный «на устойчивость»: образование колоний у коллективных амеб Dictyostelium discoideum. Этот процесс[151] интересен как пример явления, пограничного между одноклеточной и многоклеточной биологией.
Когда запас питательных веществ в той среде, в которой живут и размножаются коллективные амебы, иссякает, происходит удивительная перестройка (рис. А): отдельные клетки начинают соединяться в колонию, насчитывающую несколько десятков тысяч клеток. Образовавшийся «псевдоплазмодий» претерпевает дифференциацию, причем очертания его непрерывно изменяются. Образуется «ножка», состоящая примерно из трети всех клеток, с избыточным содержанием целлюлозы. Эта «ножка» несет на себе круглую «головку», наполненную спорами, которые отделяются и распространяются. Как только споры приходят в соприкосновение с достаточно питательной средой, они начинают размножаться и образуют новую колонию коллективных амеб. Перед нами наглядный пример приспособления к окружающей среде. Популяция обитает в некоторой области до тех пор, пока не исчерпывает имеющиеся там ресурсы. Затем она претерпевает метаморфозу, в результате которой обретает способность передвигаться и осваивать другие области.
Исследование первой стадии образования колонии показало, что она начинается с волн перемещения отдельных амеб, распространяющихся по их популяции к спонтанно возникающему «центру притяжения». Экспериментальные исследования и анализ теоретических моделей установили, что миграция является откликом клеток на существование в среде градиента концентрации ключевого вещества — циклической АМФ, периодически испускаемого сначала амебой, ставшей центром притяжения, а затем — после срабатывания механизма задержки — и другими амебами. И в этом случае мы видим, какую важную роль играют химические часы. Как уже неоднократно подчеркивалось, они, по существу, являются новым средством связи. В случае коллективных амеб механизм самоорганизации приводит к установлению связи между клетками.
Мы хотели бы подчеркнуть еще один аспект. Образование колоний коллективных амеб — типичный пример того, что можно было бы назвать «порядком через флуктуации»: возникновение «центра притяжения», испускающего циклическую АМФ, сигнализирует о потере устойчивости нормальной питательной среды, т. е. об исчерпании запаса питательных веществ. То, что при нехватке пищевого ресурса любая амеба может начать испускание химических сигналов — циклической АМФ — и, таким образом, стать «центром притяжения» для остальных амеб, соответствует случайному характеру флуктуации. В данном случае флуктуация усиливается и организует среду.
6. Бифуркации и нарушение симметрии
Рассмотрим теперь более подробно, как возникает самоорганизация и какие процессы начинают происходить, когда ее порог оказывается превзойденным. В равновесном или слабо неравновесном состоянии существует только одно стационарное состояние, зависящее от значений управляющих параметров. Обозначим управляющий параметр через ППП (им может быть, например, концентрация вещества В в «брюсселяторе», описание которого приведено в разд. «За порогом химической неустойчивости»). Проследим за тем, как изменяется состояние системы с возрастанием значения В. Увеличивая концентрацию В, мы как бы уводим систему все дальше и дальше от равновесия. При некотором значении В мы достигаем порога устойчивости термодинамической ветви. Обычно это критическое значение называется точкой бифуркации. [На особую роль этих точек обратил внимание Максвелл, размышляя над отношением между детерминизмом и свободой выбора (см. гл. 2 разд. «Язык динамики»).]
Рис. 10. Бифуркационная диаграмма. Стационарные значения переменной Х представлены на диаграмме как функции параметра бифуркации l.. Сплошные линии соответствуют устойчивым, штриховые — неустойчивым стационарным состояниям. Чтобы достичь ветви D, необходимо выбрать начальную концентрацию Х0 выше значений X, соответствующую ветви Е.
Рассмотрим некоторые типичные бифуркационные диаграммы. В точке бифуркации В термодинамическая ветвь становится неустойчивой относительно флуктуации (см. рис. 10). При критическом значении lс управляющего параметра l система может находиться в трех различных стационарных состояниях: С, Е и D. Два из них устойчивы, третье неустойчиво. Очень важно подчеркнуть, что поведение таких систем зависит от их предыстории. Начав с малых значений управляющего параметра l и медленно увеличивая их, мы с большой вероятностью опишем траекторию АВС. Наоборот, начав с больших значений концентрации Х и поддерживая постоянным значение управляющего параметра l, мы с высокой вероятностью придем в точку D. Таким образом, конечное состояние зависит от предыстории системы. До сих пор история использовалась при интерпретации биологических и социальных явлений. Совершенно неожиданно выяснилось, что предыстория может играть роль и в простых химических процессах.
[147] Более подробно о химических колебательных системах см. в работе: Winfree A. Rotating Chemical Reactions Scientific American, 1974, v. 230, p. 82—95.
[148] Goldbeter A., Nicolis G. An Alosteric Model with Positive Feedback Applied to Glycolitic Oscillations. Progress in Theoretical Biology, 1976, vol. 4, p. 65—160; Goldbeter A., Caplan S. R. Oscillatory Enzymes. Annual Review of Biophysics and Bioengineering, 1976, vol. 5, p. 449—473.
[149] Hess В., Boiteux A., Kruger J. Cooperation of Glycolitic Enzymes. Advances in Enzyme Regulation, 1969, vol 7, p. 149—167; см. также: Hess В., Goldbeter A., Lefever R. Temporal, Spatial and Functional Order in Regulated Biochemical Cellular Systems. Advances in Chemical Physics, 1978. vol. XXXVIII, p. 363—413.
[150] Hess B. Cell Foundation Symposium, 1975, vol. 31, p. 369.
[151] Geresch G. Cell Aggregation and Differentiation in Dictyostelium Discoideum.— In: Developmental Biology, 1968, vol. 3, p. 157—197.