Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

MATHEMATICAL MODELS OF BRONCHIAL PERMEABILITY IN CLINICAL PRACTICE

Karzilov A.I. 1 Tetenev F.F. 1 Bodrova T.N. 1 Mesko P.E. 1 Kalinina O.V. 1 Yakis O.V. 1
1 Siberian state medical university, Tomsk
Health subjects (n = 15), patients with bronchial asthma (n = 30), COPD (n = 16), community-acquired pneumonia (n = 17) have been inspected, – in total was 78 person (40 men, 38 women, 15–64 years). Recorded a pneumotachogram of quiet breathing with discontinuing during 0,5 s an airflow in inspiration and in expiration. Have been calculated regression models of FEV1/FVC % and FEV1 % indexes by stepwise method, including predictors: Sex (S), Age (A), Height (H), pneumotachogram parameters (Spi, Api, Ae): FEV1/FVC % = 5,604∙S – 0,438∙A + 0,462∙H – 355,249∙Spi + 0,056∙Api + 1,671∙Ae (1); FEV1 % = 0,606 H – 665,346∙Spi (2). Parameters of models: р = 0,000 (1, 2), adjusted R2 – 97,4 % (1) and 93,5 % (2). Have been calculated sensitivity (Se), specificity (Sp) and other parameters. Cut of points were chosen at a minimum value dS = |Se – Sp| and there matched prognosticated 72 % of FEV1/FVC % (Se – 80,6 %, Sp – 78,6 %) and 89 % of FEV1 % (Se – 78,8 %, Sp – 73,3 %).
mathematical model
bronchial permeability
pneumotachogram
breathing patterns
forced expiration
bronchial asthma
COPD
1. Bühl A., Zöfel P. SPSS Version 10. Einführung in die moderne Datenanalyse unter Windows. Sankt-Petersburg: DiaSoft, 2005. 608 p.
2. Vostrikova E.A., Bagrova L.O., Osipov A.G. et al. Pulmonology, 2004, no. 4, рр. 45–50.
3. Grippi M.A. Pulmonary Patophysiology. Moscow: Binom, 1997. 344 p.
4. Karzilov A.I., Tetenev F.F., Bodrova T.N. Pulmonology, 2005, no. 3, рр. 77–82.
5. Karzilov A.I. The respiratory system’s biomechanical homeostasis and its maintenance mechanisms in normal conditions and at obstructive pulmonary diseases. Bulletin of Siberian Medicine. 2007, no. 1, pp. 13–38.
6. Karzilov A.I. The adjustment of breath stability at patients with obstructive pulmonary diseases. The dissertation author’s abstract for degree of medicine doctor’s. Barnaul, 2009, 39 p.
7. Klement R.F. Disease of respiratory organs. Moscow: Medicina, 1989, Vol. 1, pp. 49–70.
8. Clinical recommendations. Pulmonology. Moscow: GEOTAR-Media, 2007, 240 p.
9. Patent. Russia. № 200118413/14, 10.07.2000.
10. Tetenev F.F. Bjul. sib. mediciny, 2005, no. 4, pp. 14–26.
11. Fletcher R., Fletcher S., Wagner E. Clinical Epidemiology. The Essentials. 3rd ed. Moscow: Media Sfera, 1998, 352 p.
12. Global Strategy for Asthma Management and Prevention, Global Initiative for Asthma (GINA). 2011, Available from: http://www.ginasthma.org.
13. Global Strategy for t he Diagnosis, Management and Prevention of COPD, Global Initiative for Chronic Obstructive Lung Disease (GOLD). 2011, Available from: http://www.goldcopd.org.
14. Hage R., Aerts J.G.J.V., Verbraak A.F.M. et al. Eur. Respir. J., 1995, no. 8, pp. 1910–1914.

Бронхиальная проходимость является ведущим фактором, определяющим вентиляционные возможности аппарата внешнего дыхания и жизнедеятельность человека в норме, при бронхообструктивных нарушениях [3, 7, 10, 12, 13] и относится к параметрам гомеостатической регуляции в аппарате внешнего дыхания [5].

Референтным методом диагностики НБП у больных бронхиальной астмой и хронической обструктивной болезнью легких, согласно международным рекомендациям Глобальной стратегии лечения и профилактики бронхиальной астмы (GINA, редакция 2011) [12], Глобальной стратегии диагностики, лечения и профилактики хронической обструктивной болезни легких (GOLD, редакция 2011) [13] и российским клиническим рекомендациям по пульмонологии, 2007 [8], являются параметры форсированного выдоха за первую секунду - ОФВ1, соотнесенного в процентах с должной величиной - ОФВ1 % и процентного отношения к форсированной жизненной емкости легких - ОФВ1/ФЖЕЛ %.

Результаты пробы с форсированным выдохом зависят не только от состояния механических свойств бронхолегочной системы, но и целого ряда внелегочных факторов: ригидность грудной клетки; регуляция тонуса и сокращения дыхательной мускулатуры; состояние дыхательного центра; утомление, желание и способность больного выполнить дыхательный маневр и др. [3, 7, 10].

В условиях скрининга респираторной функции показатели форсированного выдоха являются не референтным методом, а способом прогнозирования НБП с низким процентом верификации НБП - так, при найденном снижении ОФВ1 % у работающих людей на промышленном предприятии диагноз был верифицирован у 22,7 % мужчин и 8,9 % женщин, сниженной пиковой скорости выдоха - у 18 % мужчин и 7,4 % женщин [2].

Поэтому необходимо проводить поиск и изучение простых, необременительных способов прогнозирования НБП, которые в последующем могут стать объектом внедрения инновационных технологий. Таким методом может быть пневмотахограмма спонтанного дыхания с прерыванием воздушного потока. Показано, что у здоровых лиц и больных с бронхообструктивным синдромом наблюдаются различные значимые корреляционные связи между паттернами спокойного дыхания [4], на пневмотахограмме спокойного дыхания после фазы прерывания воздушного потока регистрируется «дельта-пик», выраженность которого коррелирует со степенью бронхиальной обструкции [14].

Цель исследования: изучить регрессионные модели, отражающие зависимость показателей бронхиальной проходимости ОФВ1 % и ОФВ1/ФЖЕЛ % от параметров паттернов спокойного дыхания вдоха и выдоха при прерывании воздушного потока, антропометрических данных и оценить их диагностические возможности в прогнозировании нарушения бронхиальной проходимости.

Материал и методы исследования

Исследование было проведено на смешанной группе лиц, включающей здоровых некурящих людей (n = 15), больных бронхиальной астмой (БА) (n = 30), хронической обструктивной болезнью легких (ХОБЛ) (n = 16) и внебольничной пневмонией (ВП) (n = 17), - всего 78 человек, из них мужчин - 40, женщин - 38, возраст 15-64 года.

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

Референтное исследование бронхиальной проходимости проводилось при помощи пневмотахографа (Медфизприбор, Казань). Определялись ФЖЕЛ, л, ОФВ1, л. Рассчитывали отношение ОФВ1/ФЖЕЛ % в процентах; процентное отношение ОФВ1 к должной величине - ОФВ1 %, рассчитанной по формулам Клемента Р.Ф. и соавт. [7].

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

19 

Рис. 1. Схема оригинальной установки проведенного исследования: A - автоматический прерыватель воздушного потока; B - пневмотахографическая трубка Вотчала; C - блок манометров; D - многоканальный регистратор; Pa - ротовое давление, Ṽ - пневмотахограмма

В условиях относительного покоя исследуемого проводили запись при скорости движения ленты самописца 100 мм∙с-1 кривой пневмотахограммы спонтанного дыхания при однократном прерывании воздушного потока на 0,5 с в фазу вдоха и выдоха в 3-4 дыхательных циклах. Исследование продолжалось 2-3 мин и не вызывало каких-либо неприятных ощущений у обследуемых лиц.

Из зарегистрированных циклов выбирали средний и рассчитывали параметры пневмотахограммы в периоды вдоха и выдоха (рис. 2): продолжительность основных пиков в фазу вдоха (Lin) до прерывания воздушного потока Li, после прерывания Lmi, с; аналогично в фазу выдоха (Lex) - Le и Lme, с; дополнительного пика в фазу вдоха Lpi и выдоха Lpe, с; амплитуда этих пиков в фазу вдоха: Hi, Hmi, Hpi и выдоха: He; Hme; Hpe, л∙с-1; площадь пиков на вдохе: Si = 0,5∙Hi∙Li, Spi = 0,5∙Hpi∙Lpi и выдохе: Se = 0,5∙He∙Le, Spe = 0,5∙Hpe∙Lpe, л; скорость отклонения кривой пневмотахограммы от изолинии на вдохе - Ai и выдохе - Ae, скорость приближения к изолинии на вдохе Api, Ami и выдохе Ape, Ame, л∙с-2. При регрессионном анализе также учитывались рост Р, см; вес W, кг; возраст B, годы; пол П, 1 - мужской, 2 - женский.

Математические модели прогнозируемых показателей строились при помощи множественного линейного регрессионного анализа с использованием программного пакета SPSS (Release 14, SPSS Inc., США). Проверку на нормальность распределения значений признака проводили по критерию Колмогорова-Смирнова.

20 

Рис. 2. Анализируемые параметры пневмотахограммы спокойного дыхания во время прерывания воздушного потока в фазу вдоха (Lin) и выдоха (Lex). Объяснение в тексте

Операционные характеристики теста - чувствительность (Se), специфичность (Sp), прогностичность положительного (Pp) и отрицательного результата (Pn) [11] рассчитывали по точкам разделения теста, поиск которых проводили по методике расчета минимальных абсолютных значений dS = |Se - Sp|, предложенной [6].

Результаты исследования и их обсуждение

Построение регрессионных моделей зависимых переменных ОФВ1/ФЖЕЛ % и ОФВ1 % проводилось пошаговым методом [1]. Распределение фактических значений ОФВ1/ФЖЕЛ % и ОФВ1 % соответствовало нормальному закону. В качестве предикторов использовали параметры антропометрии - Р, W, B, П и пневмотахограммы - Hpi, Lpi, Hmi, Hi, Lmi, Li, Si, Smi, Spi, Api, Ai, Ami, Hpe, Lpe, Hme, He, Lme, Le, Se, Sme, Spe, Ape, Ae, Ame. Были получены следующие уравнения регрессии:

ОФВ1/ФЖЕЛ % = 5,604∙П - ‒ 0,438∙В + 0,462∙P - 355,249∙Spi + + 0,056∙Api + 1,671∙Ae; (1)

ОФВ1 % = 0,606∙P - 665,346∙Spi. (2)

При дисперсионном анализе моделей наблюдался высокий уровень статистической значимости (p = 0,000). Процент отражения моделью изменчивости прогнозируемых переменных (скорректированный R-квадрат) составил 97,4 % (1) и 93,5 % (2). Гистограммы остатков соответствовали нормальному закону распределения. Распределение прогнозируемых значений ОФВ1/ФЖЕЛ %, ОФВ1 % соответствовало нормальному закону.

Операционные характеристики теста чувствительность (Se), специфичность (Sp), прогностичность положительного (Pp) и отрицательного результата (Pn), рассчитанные по методике [6], соответствовали прогнозируемым 72 % ОФВ1/ФЖЕЛ % (Se - 80,6 %, Sp - 78,6 %, Pp - 76,3 %, Pn - 82,5 %) и 89 % ОФВ1 % (Se - 78,8 %, Sp - 73,3 %, Pp - 68,4 %, Pn - 82,5 %). Полученные данные свидетельствуют о достаточно хороших диагностических возможностях моделей прогнозирования нарушения бронхиальной проходимости по математическим моделям ОФВ1/ФЖЕЛ % и ОФВ1 %.

На данный способ определения нарушения бронхиальной проходимости был получен патент Российской Федерации на изобретение [9].

Преимущества предлагаемого метода: одновременное прогнозирование показателей ОФВ1/ФЖЕЛ % и ОФВ1 %; не требуется сотрудничество с обследуемым при минимальной функциональной нагрузке на него; непродолжительность исследования; может применяться в клинических условиях у больных, страдающих пневмонией, бронхиальной астмой, ХОБЛ, которые не могут выполнить спирометрические диагностические пробы, а также у тяжелобольных, находящихся в клиностатическом положении; возможность продолжительного мониторирования состояния бронхиальной проходимости у больных в отделениях интенсивной терапии. Данная технология в качестве инновационной может быть использована в специализированных респираторах, предназначенных для лиц, находящихся в экстремальных и чрезвычайных ситуациях, связанных с присутствием бронхоконстрикорных веществ в окружающей среде.

Поскольку проведенное исследование не является популяционным (высокий преваленс нарушения бронхиальной проходимости, недостаточное число обследуемых), необходимо в условиях скрининга респираторной функции изучить диагностическую ценность прогнозирования показателей ОФВ1/ФЖЕЛ %, ОФВ1 % по параметрам пневмотахограммы и антропометрии.

Выводы

1. Регрессионные модели ОФВ1/ФЖЕЛ % (R2 - 97,4 %) и ОФВ1 % - (R2 - 93,5 %), построенные по параметрам роста, возраста, пола и показателям пневмотахограммы вдоха и выдоха при прерывании воздушного потока, при высоких уровнях статистической значимости (p = 0,000), полноценно отражают изменчивость прогнозируемых переменных.

2. Операционные характеристики свидетельствуют о достаточно хороших диагностических возможностях полученных моделей в прогнозировании нарушения бронхиальной проходимости по параметрам ОФВ1/ФЖЕЛ % и ОФВ1 %.

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

Рецензенты:

Черногорюк Г.Э., д.м.н., профессор, зав. кафедрой госпитальной терапии с курсом физической реабилитации и спортивной медицины ГБОУ ВПО «Сибирский государственный медицинский университет», г. Томск;

Букреева Е.Б., д.м.н., профессор кафедры внутренних болезней педиатрического факультета ГБОУ ВПО «Сибирский государственный медицинский университет», г. Томск.

Работа поступила в редакцию 14.12.2012.