Abstract

Abstract

A model describing testosterone regulation by distributed input feedback control is considered. It is aimed to hold testosterone concentration above a corresponding level. The distributed feedback control with integral term is proposed. This leads us to analysis of an integro-differential system which describes the connection between the concentration of hormone GnRH, the concentration of the hormone LH and the concentration of testosterone Te. We propose controls in the form of integral operators, i.e. distributed feedback controls, which allows us to hold a suitable testosterone level without loss stability of the biological system.

Full Text

ВВЕДЕНИЕ Математическое моделирование физиологических систем требует понимания биологических процессов. Различные примеры таких моделей можно найти в известной литературе [21]. Для объяснения таких процессов мы выбрали путь трансдукции сигнала, регулируемого осью гипоталамус - гипофиз - гонады (HPG), он отвечает за метаболизм гонадных гормонов, включая тестостерон, стероидный гормон, который вырабатывается в интерстициальных клетках Лейдига в яичко. Колебания активности оси HPG вызывают изменения в регуляции и синтезе гонадных гормонов и могут влиять на гомеостаз организма. Этот сложный путь охватывает ряд химических и биологических явлений и может быть разделен на три этапа (рисунок). Первый этап - на уровне мозга. Гипоталамус вырабатывает гормон, высвобождающий гонадотропин (GnRH), который активирует лютеинизирующий гормон (LH) в гипофизе. Второй этап отражает высвобождение LH в кровоток. Кровоток переносит гормон LH в клетки Лейдига. Третий, последний, этап начинается с активации каскада биохимических событий в клетках Лейдига, что приводит к выработке и последующему выделению тестостерона. Этот путь, как и многие другие в нашем организме, является циклическим и включает в себя механизм контроля обратной связи. Таким образом, когда уровень тестостерона повышается, гипоталамус получает сигнал о достаточном уровне гормона и прекращает продуцировать GnRH, который впоследствии ингибирует высвобождение LH и в результате приводит к снижению концентрации тестостерона. После снижения уровня тестостерона в мозг поступает сигнал о возобновлении процесса. Модели регуляции тестостерона изучались в [3, 5, 6, 17, 20, 21]. Рис. Путь трансдукции сигнала, регулируемый осью гипоталамус - гипофиз - гонады (HPG) В случае, если сигнал нестабилен или не действует, мы добавляем контроль обратной связи для достижения повышения уровня тестостерона. Проблема состоит в том, чтобы сделать это, не разрушая стабильность системы. Это искажение очень опасно и может стать причиной последствий и даже смерти. В нашей модели следует повысить уровень тестостерона до нормального. Контроль может разрушить стабильность модели, что приведет к непредсказуемым реакциям организма, и это представляет опасность для пациентов. Необходимо выбрать параметры нашего распределенного контроля, а именно дозу или интенсивность усвоения лекарства в организме человека в такой форме, чтобы стабильность нашей системы сохранялась. Цель этой статьи - предложить простые неравенства между коэффициентами системы (1), чтобы избежать опасности передозировки. Статья построена следующим образом. Сначала описывается модель регулирования тестостерона. Далее представлена идея подхода и предложен метод, сводящий систему интегро-дифференциальных уравнений к соответствующей системе обыкновенных дифференциальных уравнений. Описание модели Физиологические системы обычно очень сложны. Упрощая модель регуляции тестостерона, мы изобразили путь трансдукции сигнала, инициируемый в мозге, который приводит к выработке тестостерона в клетках Лейдига. Этот биологический процесс регуляции уровня тестостерона может быть представлен в виде контура. Рассмотрим модель регуляции тестостерона с использованием моделирования в [5] наследственной импульсной системы обратной связи и на основе концепций, представленных в [3, 6, 20, 21]. Модель описывает взаимодействие: (I) концентрации гормона (GnRH), который будет обозначаться как с (II) концентрацией гормона (LH) - и (III) концентрацией тестостерона (Te) - ее можно записать в виде (1) (2) где Φ и F - неубывающие и невозрастающие функции соответственно; - соответствующие коэффициенты, а α - скорость усвоения лекарственного средства. Значения соответствуют периодам полураспада GnRH, LH и Te. В здоровом мужском организме все элементы, участвующие в гомеостазе, работают по согласованию. Различные негативные факторы могут разрушить отношения в системе. Предположим, что связь и т.е. импульс, определенный в (2), не действует. В этом случае мы должны добавить дополнительный контроль, благодаря которому возникнет гомеостаз системы. Мы хотели бы предложить механизм, позволяющий удерживать тестостерон на нормальном уровне, хотя в таком случае в биологической системе импульсы не поступают в мозг вовремя. Добавляем элемент управления в форму (3) в правой части второго или третьего уравнения. Функция T(t) в разности под интегралом представляет концентрацию тестостерона, которую мы хотим сохранить. Ясно, что интегральный член увеличивает концентрацию гормона LH, если концентрация Te меньше, чем T(t). Контроль (3) представляется разумным и может быть реализован в соответствии с мнениями клинических экспертов. Подставляя это управление во второе или третье уравнение, мы получаем систему (4) где (5) В случае мы получаем систему, которая описывает контроль уровня тестостерона путем инфузии препарата через систему кровообращения. Во втором случае мы находим систему, которая описывает контроль уровня тестостерона путем введения препарата в непосредственной близости от железы. Мы должны продемонстрировать, что добавление элементов управления не разрушает стабильность. МЕТОД Сложность метода Неавтономные интегродифференциальные уравнения (IDE) (6) где A(t) и k(t, s) - матрицы n × n с непрерывными коэффициентами, которые моделируют многие процессы в медицине и технике. В работах [14, 17] было продемонстрировано, что интегродифференциальные уравнения могут использоваться для моделирования эндокринной регуляции в соответствии с задержкой, возникающей при транспортировке гормона от сайта секреции к рецепторам. Это обычно мотивируется (см., например, [12]) тем, что передача сигналов рецепторами чувствительна к среднему значению концентрации гормонов в течение определенного периода времени, а не к мгновенному значению. Хотя управление с распределенным вводом является часто возникающей сложной проблемой, ей посвящено лишь несколько работ (см., например, работы [2, 12, 19]). Шум в управлении с обратной связью по задержке является основным препятствием, возникающим в математических моделях с распределенными входами: наше управление невозможно основать, используя параметры процесса X(t) только в момент и мы должны использовать среднее значение процесса в соответствующей окрестности Другой причиной появления интегральных слагаемых в такого рода медицинских моделях является то, что для усвоения лекарства организмом требуется время, это некий отдельный процесс со своей внутренней динамикой. Интегральный член с ядром, определяющим вес каждого значения, выполняет эту функцию. Из медицинской практики также ясно, что длительное наблюдение помогает эффективно контролировать лечение пациентов, если принять во внимание множество значений соответствующих параметров, описывающих процесс. Распределенные управления в виде интегральных терминов становятся удобным математическим объектом в этом естественном случае. В [14] и [24] отмечается, что такие модели с распределенными входами могут использоваться при исследовании динамики заселенности, двигателей ракетного топлива и сетевых систем управлений, в [11, 13] - вязкоупругих свойств полимеров. Достаточные условия устойчивости для интегродифференциальных уравнений (6) были получены в нескольких книгах (см., например, [1, 10, 21, 23]). В этой статье мы объединяем анализ устойчивости на основе положительности и специальное приведение интегродифференциальных уравнений к системе обыкновенных дифференциальных уравнений. Отметим, что положительный анализ устойчивости систем с запаздыванием использовался, например, в [8-10, 15, 16, 22, 23]. Прелиминарии Отметим несколько важных для нас результатов позитивного подхода к устойчивости (см., например, [1, 8-10, 15, 16, 22, 23]). Положительность на для системы обыкновенных дифференциальных уравнений (7) означает, что из неравенств (8) Это следует из того, что (9) Если неравенства (8), (9) выполняются для (т.е. ), мы называем систему положительной. Используя нижнюю и верхнюю вектор-функцию и можно гораздо проще найти решение системы неравенств, чем решение системы уравнений. Монотонный интегративный процесс для приближенного решения, который начинался с вектор-функций z(t) и y(t), был предложен в [4]. Заметим, что свойство (8) ⇒ (9) может быть описано с использованием альтернативной терминологии: полупоток сохраняет порядок решений (см. [10]) или положительных систем (см. [23]). Определение Матрица называется Метцлеровой (см., например, [9, 22]), если все ее недиагональные элементы неотрицательны, т.е. (10) T. Wazewski [25] доказано, что для системы (7) обыкновенных дифференциальных уравнений условие (10) необходимо и достаточно для ее положительности (7), т.е. для свойства (8) (9). Автономная система (11) где - постоянные для асимптотически устойчивы (а также экспоненциально устойчивы) тогда и только тогда, когда матрица коэффициента - Гурвицева (все собственные значения отрицательны). Следующее утверждение важно для устойчивости (см., например, [22] или теорему 16.1 в [1]). Лемма 1 Если матрица постоянна и метцлерова, то следующие два утверждения эквивалентны: 1. Существует постоянный вектор с положительными компонентами, так что все компоненты постоянного вектора являются отрицательными. 2. Система обыкновенных дифференциальных уравнений (13) экспоненциально устойчива. Рассмотрим систему (12) где Следующий факт можно получить на основе известного результата (см. [18]). Лемма 2 Если система (12) экспоненциально устойчива, то и система (7) экспоненциально устойчива. Метод сокращения Наш подход по существу основан на следующей идее, сводящей систему интегродифференциальных уравнений к системам обыкновенных дифференциальных уравнений. Рассмотрим систему (13) где неотрицательные константы, а система (14) с начальными условиями Теорема 1 Решение системы (13) и первые три компонента вектора решения системы (14) совпадают. РЕГУЛЯЦИЯ УРОВНЯ ТЕСТОСТЕРОНА Регулирование путем введения препарата в непосредственной близости от железы Подставляя управление (3) в третье уравнение, т.е. в случае получаем систему (15) где - положительные константы и (16) Понятно, что анализ устойчивости можно сделать для системы двух уравнений (17) Следующее утверждение об устойчивости системы (15) будет доказано. Теорема 2 Если тогда система (15) экспоненциально устойчива. Доказательство Используя теорему 1, систему (15) можно свести к системе обыкновенных дифференциальных уравнений (18) с начальными условиями (19) Вектор решения системы (17) и компоненты вектора решения системы (18) совпадают (см. теорему 1). Принимая во внимание теорему Боля - Перрона (см., например, теорему B.20 [1]), можно сделать вывод, что анализ экспоненциальной устойчивости может быть сведен к исследованию системы (20) Чтобы использовать подход, основанный на положительных системах, мы рассмотрим (21) Для экспоненциальной устойчивости системы (21) необходимо и достаточно, согласно лемме 1, чтобы все компоненты вектора-решения алгебраической системы (22) были положительны. Используя классическую теорему Крамера, получим (23) (24) (25) (26) (27) Ясно, что из неравенства в формулировке теоремы 2 следует положительность всех компонент вектора решения. Ссылка на лемму 1 завершает доказательство. Следствие 1 Система (15) всегда может быть стабилизирована с помощью управления (3). Чтобы объяснить это, мы выбираем Условие теоремы 2 будет выполнено. Замечание 1 Значения соответствуют периодам полураспада GnRH, LH и Te и могут быть рассчитаны по биомедицинским данным для типичных уровней гормонов, как указано в литературе [8]: Мы видим, что система (15) устойчива, если Последнее неравенство показывает, что уровень усвоения препарата должен быть достаточно высоким. Регулирование путем контроля концентрации лютеинизирующего гормона Следующее утверждение об устойчивости системы (4) будет доказано при условии, что С медицинской точки зрения это означает, что мы контролируем уровень тестостерона, вводя лекарство через систему кровообращения (28) где являются положительными константами. Аналогично теореме 1 мы можем получить следующее утверждение. Теорема 3 Если тогда система (28) экспоненциально устойчива. Следствие 2 Если (28) выполнено, система (28) всегда может быть стабилизирована управлением (3). Чтобы доказать это, мы выбираем Условие теоремы 3 будет выполнено. Заключение В этой работе определены параметры распределенного контроля, а именно доза и интенсивность всасывания лекарства в организме человека в такой форме, что стабильность системы будет сохранена. Другими словами, предложены простые неравенства между коэффициентами системы (4) во избежание ее неустойчивости.