- •ПРЕДИСЛОВИЕ
- •ОПРЕДЕЛЕНИЕ И НАЗНАЧЕНИЕ МОДЕЛИРОВАНИЯ
- •1.1. ЧТО ТАКОЕ МОДЕЛЬ?
- •Место моделирования среди методов познания
- •Определение модели
- •Определение модели
- •Цели моделирования
- •1.2. КЛАССИФИКАЦИЯ МОДЕЛЕЙ
- •Материальное моделирование
- •Идеальное моделирование
- •Когнитивные, концептуальные и формальные модели
- •1.3. Классификация математических моделей
- •Классификационные признаки
- •Классификация математических моделей в зависимости от оператора модели
- •Классификация математических моделей в зависимости от целей моделирования (рис. 1.11)
- •Классификация математических моделей в зависимости от методов реализации (рис. 1.12)
- •Задания для самостоятельного выполнения
- •Глава 2
- •ЭТАПЫ ПОСТРОЕНИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ
- •2.1. ОБСЛЕДОВАНИЕ ОБЪЕКТА МОДЕЛИРОВАНИЯ
- •2.2. КОНЦЕПТУАЛЬНАЯ ПОСТАНОВКА ЗАДАЧИ МОДЕЛИРОВАНИЯ
- •2.3. МАТЕМАТИЧЕСКАЯ ПОСТАНОВКА ЗАДАЧИ МОДЕЛИРОВАНИЯ
- •2.4. ВЫБОР И ОБОСНОВАНИЕ ВЫБОРА МЕТОДА РЕШЕНИЯ ЗАДАЧИ
- •stop
- •2.5. РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ В ВИДЕ ПРОГРАММЫ ДЛЯ ЭВМ
- •2.6. ПРОВЕРКА АДЕКВАТНОСТИ МОДЕЛИ
- •2.7. ПРАКТИЧЕСКОЕ ИСПОЛЬЗОВАНИЕ ПОСТРОЕННОЙ МОДЕЛИ И АНАЛИЗ РЕЗУЛЬТАТОВ МОДЕЛИРОВАНИЯ
- •Задания для самостоятельного выполнения
- •ПРИМЕРЫ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
- •3.1. СТАТИЧЕСКИЙ АНАЛИЗ КОНСТРУКЦИЙ
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Математическая постановка задачи
- •Решение задачи
- •Анализ результатов
- •Концептуальная постановка
- •Математическая постановка задачи
- •Методика решения задачи
- •Анализ результатов
- •3.3. ДИНАМИКА ПОПУЛЯЦИИ
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Решение задачи
- •Анализ результатов
- •Математическая постановка задачи для модели Ферхюльста
- •Решение задачи
- •Анализ результатов
- •Численное исследование модели Ферхюльста
- •3.4. МОДЕЛЬ КОНКУРЕНЦИИ ДВУХ ПОПУЛЯЦИЙ
- •Математическая постановка задачи
- •Качественный анализ задачи
- •Численное исследование модели конкуренции популяций
- •3.5. ГАРМОНИЧЕСКИЙ ОСЦИЛЛЯТОР
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Математическая постановка задачи
- •Решение задачи
- •Качественный анализ задачи
- •Численное исследование модели
- •Качественный анализ задачи
- •Решение задачи
- •Задания для самостоятельного выполнения
- •Глава 4
- •СТРУКТУРНЫЕ МОДЕЛИ
- •4.1. ЧТО ТАКОЕ СТРУКТУРНАЯ МОДЕЛЬ?
- •4.2. СПОСОБЫ ПОСТРОЕНИЯ СТРУКТУРНЫХ МОДЕЛЕЙ
- •4.3. ПРИМЕРЫ СТРУКТУРНЫХ МОДЕЛЕЙ
- •Задания для самостоятельного выполнения
- •МОДЕЛИРОВАНИЕ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ
- •5.1. ПРИЧИНЫ ПОЯВЛЕНИЯ НЕОПРЕДЕЛЕННОСТЕЙ И ИХ ВИДЫ
- •5.2. МОДЕЛИРОВАНИЕ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ, ОПИСЫВАЕМОЙ С ПОЗИЦИЙ ТЕОРИИ НЕЧЕТКИХ МНОЖЕСТВ
- •5.3. МОДЕЛИРОВАНИЕ В УСЛОВИЯХ СТОХАСТИЧЕСКОЙ НЕОПРЕДЕЛЕННОСТИ
- •5.4. МОДЕЛИРОВАНИЕ МАРКОВСКИХ СЛУЧАЙНЫХ ПРОЦЕССОВ
- •Задания для самостоятельного выполнения
- •ЛИНЕЙНЫЕ И НЕЛИНЕЙНЫЕ ЯВЛЕНИЯ
- •6.1. О ЗАКОНЕ ГУКА И ГРАНИЦАХ ЛИНЕЙНОСТИ
- •6.3. О ПОСТРОЕНИИ СПЛОШНЫХ МОДЕЛЕЙ. ВЫВОД ВОЛНОВОГО УРАВНЕНИЯ ИЗ ЗАКОНОВ МЕХАНИКИ
- •6.4. РЕШЕНИЕ ВОЛНОВОГО УРАВНЕНИЯ МЕТОДОМ ФУРЬЕ
- •6.6. УРАВНЕНИЯ МАКСВЕЛЛА
- •6.7. О КЛАССИФИКАЦИИ КВАЗИЛИНЕЙНЫХ СИСТЕМ
- •6.8. СВЯЗЬ НЕПРЕРЫВНОГО И ДИСКРЕТНОГО НА ПРИМЕРАХ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ И УРАВНЕНИЯ ШРЕДИНГЕРА
- •6.9. О ПОЛЬЗЕ ФЕНОМЕНОЛОГИИ ПРИ ПОСТРОЕНИИ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
- •6.10. АНАЛИЗ ПОДОБИЯ И РАЗМЕРНОСТИ
- •6.11. АВТОМОДЕЛЬНОСТЬ
- •6.12. САМООРГАНИЗАЦИЯ И СТРУКТУРЫ В НЕЛИНЕЙНЫХ СРЕДАХ
- •6.13. О НЕЛИНЕЙНЫХ ВОЛНАХ В СПЛОШНЫХ СРЕДАХ
- •6.14. ИЕРАРХИЧЕСКИЕ МОДЕЛИ ТУРБУЛЕНТНОСТИ И МНОГОМАСШТАБНЫЕ ФУНКЦИОНАЛЬНЫЕ БАЗИСЫ
- •Иерархический базис для турбулентных полей
- •Одномерный иерархический базис
- •Двумерный базис
- •6.15. ВЕЙВЛЕТЫ
- •Непрерывное вейвлет-преобразование
- •6.16. ВЕЙВЛЕТ-АНАЛИЗ ВРЕМЕННЫХ КОЛЕБАНИЙ
- •6.17. О ФРАКТАЛАХ И ИХ ПРИМЕНЕНИИ
- •Примеры фракталов
- •Подобие и скейлинг
- •Множества Мандельброта и Жюлиа
- •Фрактальная размерность кластеров
- •Экспериментальные методы определения фрактальной размерности
- •Модель случайных фракталов для описания растущих дендритных структур
- •Результаты применения модели случайных фракталов
- •6.18. НЕЛИНЕЙНЫЕ МОДЕЛИ ДНК
- •Структура и физические свойства ДНК
- •Модель Пейрара-Бишопа
- •Задания для самостоятельного выполнения
- •Глава 7
- •7.2. ИМИТАТОР СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ
- •7.3. КЛЕТОЧНЫЕ АВТОМАТЫ
- •7.4. МОДЕЛИРОВАНИЕ ДИСЛОКАЦИЙ В МЕТАЛЛЕ
- •Самоорганизация дислокаций в модели клеточных автоматов
- •Задания для самостоятельного выполнения
- •ЯЗЫК ФОРМАЛЬНОГО ОПИСАНИЯ АЛГОРИТМОВ
- •====== Приложение 2
- •П2.1. Решение уравнений высоких степеней и трансцендентных уравнений с одним неизвестным
- •П2.2. Решение систем линейных уравнений
- •БИБЛИОГРАФИЧЕСКИЙ СПИСОК
- •ПРЕДМЕТНЫЙ УКАЗАТЕЛЬ
- •ОГЛАВЛЕНИЕ
- •Введение в математическое моделирование
6.18. НЕЛИНЕЙНЫЕ МОДЕЛИ ДНК
Малые группы атомов, слишком малые, что бы проявлять точные статистические законы, иг рают главенствующуюроль в весьма упорядочен ных и закономерных процессах внутри организма. Они управляют видимыми признаками большого масштаба, которые организм приобретает в те чение своего развития; они определяют важные особенности его функционирования, и во всем этом проявляются весьма отчетливые и строгие биологические законы.
Э-Шредингер
Структура и физические свойства ДНК
Поскольку реальные механизмы многих биологических про цессов, связанных с ДНК, до настоящего времени не поняты, биологи обращают все большее внимание на структурные аспек ты с целью объяснения основных механизмов функционирования ДНК —транскрипции, репликации и денатурации. Принимая во внимание тот факт, что симметрия ДНК при ее формировании может быть связана с дефектами, эффекты дисторсии в структуре двойной спирали ДНК могут также рассматриваться как зарожде ние и эволюция новых структурных дефектов. Двойная спираль ДНК представляет собой уникальный полимер, в котором моле кулярные цепи образуют так называемую В-конформацию, впер вые установленную по рентгеновским диффракционным картинам. В соответствие с результатами Уотсона и Крика В-структура ДНК имеет следующие основные черты [117]:
1. В состав входят две полинуклеотидные цепи (рис. 6.43), об вивающие с правосторонним поворотом общую ось, образуя двой ную спираль с диаметром ~ А. Двойные цепи являются антипарал-
Рис. 6.43. Структура двойной спирали
лельными и обернуты таким образом друг относительно друга, что не могут быть разделены без взаимного расплетения (эффект плектонемической закрутки). Основания покрывают поверхность спирали, которая обвита сахаро-фосфатными цепями на перифе рии, причем так, что достигается минимизация сил отталкивания между заряженными фосфатными группами.
2. Плоскости оснований почти перпендикулярны оси спира ли; каждая пара расположенных друг против друга оснований имеет водородную связь и образует при этом плоскую пару осно ваний (base pair). Взаимодействие оснований посредством водород ных связей (complementary base pairing phenomenon) и является причиной специфического объединения двух цепей в двойную спираль. Идеальная В-ДНК спираль содержит 10 пар оснований (base pairs, pb) на один поворот (поворот спирали на 1 pb состав ляет 36°), что обеспечивает шаг спирали (rise per turn) и простран ственную периодичность угла закрутки соответственно А = 34 А и
сор = 2я/А = 1.85 нм-1 Ароматические основания имеют ван-дер-
ваальсовский масштаб 3,4 А, и, как следствие этого, наблюдаются эффекты упаковки (base stacking).
3. Персистентные длины поворотных и изгибных мод пред ставляют естественные масштабы для термических флуктуаций (приблизительно 15—20 А). ДНК является достаточно жестким по лимером вследствие растяжения соответствующих пар аминокис лот химическими и водородными связями. Термические персис тентные длины ДНК на изгиб и вращение составляют соответ ственно 150 Ьр (50 нм) и 220 Ьр (75 нм).
Замечательной чертой структуры ДНК, по Уотсону-Крику, яв ляется аккомодация двух типов пар оснований: остаток аденина объединяется с остатком тимина и каждый остаток гуанина —с остатком цитозина. Разнообразие протеинов обеспечивает специ фические микроскопические изменения в двойной спирали.
4. При пространственных масштабах в тысячи пар оснований (~ 0,5 мкм) ДНК in vivo организованы в топологически незави симые петли. Молекулы ДНК с размерами 2—10 kbp (0,6—3,3 мкм) представляют собой «связанные» множества. Основная часть генов (prokariotic cells) организована как единичные замкнутые ДНК с размерами ~2 МЬр.
Статистические модели Д Н К . Структурные переходы в Д Н К
Имеет место формальная аналогия между механизмами неус тойчивости при фазовых переходах (melting instability) в однород ной структуре ДНК и фазовыми переходами в однородных систе мах. Эффекты плавления (денатурации) в растягиваемой двойной цепи ДНК (рис. 6.44) —одна из фундаментальных проблем в ста тистической физике биологических систем —обнаруживают раз личные типы критического поведения, наблюдаемые и других тер модинамических системах [119,129].
Рис. 6.44. Распад ДНК
Ранние модели термического распада ДНК, т.е. разделения двух цепей при нагреве [121], были основаны на моделях типа Изинга, в которых параметр порядка, характеризующий взаимо действие оснований, сводился к расстоянию между последними. Современные исследования связывают эти переходы с нелинейны ми механизмами зарождения конечно-амплитудных флуктуаций [120]. Анализ закономерностей скейлинга указанных переходов был проведен в [128] с использованием асимметричного потенциала, содержащего пологий участок для больших значений перемеще ний, что соответствует распаду двойной спирали при повышенных температурах. Этот тип потенциала включал также взаимодействие между последовательными основаниями в длинноволновом при ближении. В соответствии с этим описанием локальная неустой
чивость плавления возникает в «мягком» режиме, что не подтвер ждается экспериментальными наблюдениями динамики переходов, которые носят резкий характер.
Обобщение данной модели, предложенное в [127], связано с введением параметра жесткости упаковки (stacking parameters), ко торый, учитывая взаимодействие последовательных пар оснований, позволил описать резкий характер перехода.
Исследования однородных ДНК также показали, что указан ная модель дает возможность описать последовательный распад (денатурацию), наблюдаемый в эксперименте. В дополнение уста новлено, что в рамках модели возможна генерация «энтропийно го барьера», с которым связывается «сужение» перехода. При не которых значениях параметра жесткости модель демонстрирует признаки переходов первого рода с конечным значением энтропии плавления, нарушением непрерывности для части связей (что ти пично для ДНК), расходимость продольных и поперечных кор реляционных масштабов.
Степени свободы молекулы Д Н К
Степени свободы, описывающие дисторсии в молекуле ДНК, могут быть разделены на переменные, характеризующие конфор мации центральной оси, и переменные, определяющие внутрен ний поворот двойной спирали [124]. Положение оси ДНК может быть описано пространственной кривой, заданной параметром длины вдоль оси ДНК и радиусом-вектором точки г (s). Касатель
ная t s Э / задается единичным вектором t, положение которого в
пространстве определяется двумя углами или локальной кривизной и дисторсией.
Поворот измеряется как пространственная интенсивность по ворота пар оснований относительно центральной оси, которая для неискаженной структуры ДНК задается как со0. Отклонение в ин тенсивности поворота от to0 может быть описано с использовани ем скалярного поля Q(s). Это дает в целом три степени свободы каждой координате, что необходимо для описания бесконечно-ма лых поворотов координатных осей, связанных с парами оснований при перемещении последних вдоль молекулы.