MODELOS MATEMTICOS DE LA FISIOLOGA RESPIRATORIA.

Autores:

Dr. Hctor Quiones Acosta

Centro de Ciberntica Aplicada a la Medicina (CECAM)
Calle: Ave. 146 esq. 31 Cubanacn, Playa.
Tel. :271-1354

Dr. Antonio Barber Gutirrez

Profesor Titular
Instituto Superior de Ciencias Mdicas de La Habana
(ISCM-H)

Dr. Roberto Douglas Pedroso

Profesor Titular

RESUMEN:

Se hace un anlisis de los principales modelos matemticos que aparecen en la literatura, en el

rea de la fisiologa respiratoria, con vistas a seleccionar algunos que sean ms factibles de ser
aplicados en la prctica clnica, no requieran de mucho equipamiento para tomar datos de
pacientes y permitan incorporar patrones de respiracin Yoga como una posible opcin en el
tratamiento. Especficamente para evaluar y tratar a personas que padecen de enfermedades
respiratorias obstructivas, los modelos de dos compartimentos alveolares y corto circuito son
nuestra propuesta. Estos modelos parecen ser ms adecuados a los objetivos mencionados, es
decir, para la evaluacin del paciente y la posible terapia natural con tcnicas respiratorias.

Palabras claves: Modelos matemticos, intercambio de gases, sistema respiratorio humano.

Respiracin Yoga.

Mathematical models of the respiratory physiology. Commented revision.

Abstract

We analyzed the main mathematical models that appeared in the literature in the area of the
respiratory physiology, in order to select some of them for applying into clinical practice, with no
much requirement to obtain patient data and at the same time give the possibility to add pattern
of Yoga breathing as a possible option for the treatment. Specifically to evaluate and treat
patients with obstructive pulmonary disease, models of two alveolar compartments and shunt
are our proposal. These models seem to be more adequate for the purpose mentioned above,
that is, for the evaluation of patients and the possible natural therapy with respiratory
techniques.

Keywords: Mathematical models, gas exchange, human respiratory system, Yoga breathing.
DESARROLLO

La modelacin matemtica penetra en el campo de la Fisiologa Respiratoria con el modelo de

Milhorn, Benton, Ross y Guyton [1]. Este es uno de los primeros modelos pulmonares, con
deducciones matemticas completas. Aqu se consideran tres reservorios: cerebro, pulmones y
tejidos corporales, en cada uno de los cuales se definen los volmenes de oxgeno (O2) y
dixido de carbono (CO2). Por su parte, los trabajos de West [2-5] tratan el complejo y todava
no completamente explorado camino en la modelacin de la ventilacin alveolar y del flujo
sanguneo capilar, as como del intercambio de gases. Especficamente, el anlisis cuantitativo
en las desigualdades de la razn ventilacin/perfusin (VA/Q), puede ayudar a predecir ndices
en las enfermedades obstructivas crnicas del pulmn.

Muchos aos transcurren con publicaciones sobre esta temtica, abordada por West, de la
distribucin (VA/Q). Citemos, por ejemplo, los trabajos con referencia [6-19]. Ya, Vidal Melo et
al. Ver [20,21], en 1990 y 1993, comienzan una modelacin matemtica ms avanzada y
complicada. Modelos de dos compartimentos de intercambio de gases en el pulmn, teniendo
en cuenta corto circuito (shunt), heterogeneidad de (VA/Q) y capacidad de difusin daada,
consideradas simultneamente. Adems, el modelo contiene el sistema dinmico que describe
la variacin del contenido de los gases en sangre arterial (para cada compartimiento) y las
presiones alveolares que pueden obtenerse, las cuales dependen de ciertos parmetros
(fraccin de ventilacin y perfusin). Estos parmetros se estiman usando tcnicas
matemticas de optimizacin (GUESS). Concretamente en el trabajo [21], se presentan dos
aplicaciones, una es un mtodo que considera los gradientes arterial-alveolar de CO2 y O2 en
dos componentes D y (VA/Q). Esta tcnica es evaluada en datos de pacientes con enfermedad
pulmonar obstructiva crnica y permite un mtodo general de diagnstico del estado del
sistema respiratorio. La posibilidad de integrar a esta idea diferentes patrones de respiracin
(inhalacin y exhalacin), con frecuencias respiratorias variables, ayudara en el tratamiento y
prevencin en caso de que el problema sea reversible.

Modelos de un compartimiento

1.Modelo del Sistema de Control Respiratorio.

En el modelo de Milhorn et al. Ver [1] las ecuaciones diferenciales que expresan las variaciones
en la dinmica del pulmn son:

dVAO2/dt = Q( CvO2 - CaO2 ) + VA ( CIO2 -CAO2 ) (1)

dVACO2/dt = Q( CvCO2 - CaCO2 ) + VA ( CICO2 -CACO2 ) (2)

Estas ecuaciones describen los cambios de los volmenes de O2 y CO2 (VAO2 y VACO2) en
el pulmn. Aqu Q es el gasto cardiaco. CvO2, CvCO2, CaO2, CaCO2, CAO2 y CACO2 son las
concentraciones venosas, arteriales y alveolares de O2 y CO2 respectivamente. VA es la
ventilacin alveolar. Las concentraciones inspiradas de O2 y CO2 se denotan por CIO2 y
CICO2. Los dos sumandos de las partes derechas de las ecuaciones (1) y (2) representan la
velocidad de cambio de los volmenes de los gases en la sangre y en el aire, respectivamente,
algunas funciones vienen dadas por relaciones fisiolgicas no lineales (ver [1]). En forma
anloga se escriben las ecuaciones para los dems reservorios.

A continuacin, todas las ecuaciones se expresan en trminos de las concentraciones de los

gases. Las soluciones correspondientes a los estados transitorio y estacionario fueron
obtenidas para diferentes perturbaciones de la entrada de O2 y CO2 utilizando una
computadora, y fueron comparadas con resultados experimentales. Esta idea es la base para la
estimacin de parmetros, en modelos pulmonares ms complejos, que se realiza en la
actualidad.

En este trabajo se investiga adems, la efectividad del sistema respiratorio como un regulador
biolgico. Se asume desde el punto de vista matemtico que el sistema de control respiratorio
humano es normal (lo cual simplifica el modelo), pero entonces para el anlisis o estudio de
personas enfermas o con problemas respiratorios el modelo no es adecuado.

2. Modelo de Mecnica Pulmonar.

En el trabajo [22], Marini presenta un enfoque de la dinmica pulmonar, en el cual describe la

variacin del volumen circulante en el intervalo de inspiracin y expiracin. Se da la posibilidad
de introducir las formas de ondas en cada fase, as como la frecuencia respiratoria. Tomando
esta idea como base, se podra analizar la influencia de la respiracin lenta y profunda de una
persona en estudio para mejorar su funcin pulmonar.

Veamos una pequea descripcin. En el modelo se plantea un balance de presiones, con las
ecuaciones diferenciales:

Ri dVi/dt +Vi/C +Pex =P(t) 0< t < ti (3)

Rx dVx/dt +Vx/C +Pex = 0 ti < t < ttot (4)

Se consideran dos fases. Inspiratoria y Espiratoria. Ri denota la resistencia inspiratoria, Vi el

volumen inspiratorio, C la complianza, Pex la presin alveolar espiratoria final, Rx la resistencia
espiratoria, Vx el volumen espiratorio, ti el tiempo de inspiracin y ttot el tiempo total.

P(t) es la funcin excitatoria (funcin de entrada), la cual se supone en el trabajo que es de tipo
lineal o de forma sinusoidal; pero tericamente puede ser cualquier funcin integrable. Con este
enfoque se pueden introducir diferentes frecuencias y patrones respiratorios.

En el modelo se considera que el sistema respiratorio se infla y se desinfla como un

compartimiento simple. Dentro de ciertos lmites esta consideracin es aproximadamente
correcta para sujetos saludables normales. Sin embargo, como reconoce el propio autor en la
crtica del modelo, esta premisa puede no ser vlida para una variedad de implementaciones
clnicas. Un ejemplo especifico es que son necesarias dos funciones exponenciales para
caracterizar la fase de espiracin en pacientes con obstrucciones del flujo areo severa. Si
nuestro objetivo es estudiar personas con posibles problemas respiratorios obstructivos, este
modelo de mecnica pulmonar no es adecuado. Siguen en esta lnea de investigacin los
artculos con referencia[23-26].

Modelos de dos compartimentos

1. Aplicaciones clnicas.

En el trabajo de Melo [27] se demuestra que la disparidad ventilacin-perfusin y la

capacidadde difusin pulmonar limitada producen un efecto del gasto cardaco (Q) sobre el
intercambio de gases. Concretamente, un mximo en la curva PaO2 vs Q puede ser alcanzado.
Es decir, ambos problemas respiratorios unidos pueden llevar al corazn a trabajar en hipoxia.
De aqu se desprende lo importante que podra ser la prctica de una respiracin correcta.
Mostremos algunos aspectos del contenido del modelo.

La ley de Fick es la que permite expresar las velocidades de cambio de los contenidos
arteriales de los gases; en este caso queda como sigue.

dCa(t)/dt = 0.1 (D/Vc) ( PA - Pa(t) )

donde Ca(t) denota el contenido en sangre arterial del gas correspondiente, D la difusin, PA,
Pa(t) presiones parciales alveolares y arteriales del gas.Para el O2 y CO2 tenemos las
ecuaciones siguientes:

dCaO2/dt = 0.1 (D/Vc) (PAO2 - PaO2) (5)

dCaCO2/dt = 0.1 (20 D/Vc) ( PACO2 - PaCO2 ) (6)

Las ecuaciones (5) y (6) se plantean para cada uno de los dos compartimentos, de acuerdo a
las fracciones de ventilacin y perfusin respectivamente. Las respuestas a las desigualdades
ventilacin perfusin pueden predecir los efectos de problemas pulmonares de carcter
obstructivo. El empleo de modelos de dos compartimentos y shunt se sigue con las
aplicaciones de Whiteley JP [28-30]. Adems, considerando el intercambio de gases, los
trabajos [31,32] son de inters.

2. Terapia con respiracin Yoga.

Las tcnicas respiratorias del Yoga como va para mejorar la salud son milenarias y continan
aplicndose en el mundo moderno. En el libro de Hermgenes [33] se explican los beneficios
de esta terapia.

Ya en el ao 1998, referencia [34] se relata un experimento realizado a un grupo de pacientes

con insuficiencia cardiaca crnica, basado en patrones de respiracin Yoga. Especficamente
muestran beneficios en la reduccin de la frecuencia respiratoria, para mejorar ciertos ndices.
ltimamente Spicuzza et al. [35] concluyen en su trabajo sobre esta temtica que la prctica de
la respiracin Yoga puede tener relevancia clnica en condiciones tales como la enfermedad
cardiaca crnica.

Este enfoque descrito anteriormente podra complementarse con el modelo de Marini [22], ya
expuesto.

CONCLUSIONES

Los modelos de dos compartimentos incluyendo la mecnica pulmonar, el intercambio de gases

y el corto circuito, parece que puedan reflejar con una buena aproximacin el estado del
sistema respiratorio de las personas que sufren de problemas de tipo obstructivo. Adems, se
hace posible la integracin del patrn respiratorio deseado al modelo matemtico. Esto
ayudara en la valoracin del estado del paciente y su tratamiento.

BIBLIOGRAFIA

1. Milhorn HT, Benton R, Ross R, Guyton AC. A Mathematical Model of the Human Respiratory Control System. Biophysical
Journal 1965 Volume 5, pag. 27-44
2. West, JB ventilation-perfusion inequality and overall gas exchange in computer models of the lung
3. Resp. Physiol. 1969. 7:88-110
4. West, J. B. Respiration. Annual Review of Physiology, 1972; vol 34, pag.91
5. West, J. B. Blood flow to the lung and gas exchange. Anesthesiology, 1974; 41, 124-138.
6. West JB New advances in pulmonary gas exchange. Anesth Analg 1975 Jul-Aug;54(4):409-18
7. Zwart A, van Dieren A A simple and non-invasive method to determine the ventilation-perfusion ratio of the lung and the
effective lung perfusion. Acta Anaesthesiol Belg 1975;23 Suppl:53-63
8. West JB Wagner PD Pulmonary Gas Exchange in Bioengineering Aspects of the Lung. ed J.B. West; Lung Biology in Health
and Disease,vol 3, Marcel Dekker Inc. 1977 pp. 361-454
9. Scrimshire DA Theoretical analysis of independent VA and Q inequalities upon pulmonary gas exchange. Respir Physiol 1977
Apr;29(2):163-78
10. Wagner PD, Dantzker DR, Iacovoni VE, Tomlin WC, West JB Ventilation-perfusion inequality in asymptomatic asthma. Am
Rev Respir Dis 1978 Sep;118(3):511-24
11. Beroff MJ, Lelong F, Cherruault Y [Pulmonary gas exchange model: influence of the heterogeneity of distribution on the
ventilation-perfusion and diffusion-perfusion ratios of oxygen transfer]. Int J Biomed Comput 1978 Jul;9(4):287-300
12. Fortune JB, Wagner PD Effects of common dead space on inert gas exchange in mathematical models of the lung. J Appl
Physiol 1979 Oct;47(4):896-906
13. Teplick R, Snider MT, Gilbert JP A comparison of continuous and discrete foreign gas VA/Q distributions. J Appl Physiol 1980
Oct;49(4):684-92
14. Haab P A model for the study of diffusion and perfusion limitation. Fed Proc 1982 Apr;41(6):2119-21
15. Stewart WE, Mastenbrook SM Jr Parametric estimation of ventilation-perfusion ratio distributions. J Appl Physiol 1983
Jul;55(1 Pt 1):37-51
16. Mertens P. A simple model of VA/Q distribution for analysis of inert gas elimination data. J Appl Physiol 1983 Aug;55(2):562-
8
17. Kapitan KS, Wagner PD Linear programming analysis of VA/Q distributions: limits on central moments. J Appl Physiol 1986
May;60(5):1772-81
18. Kapitan KS, Wagner PD Linear programming analysis of VA/Q distributions: average distribution. J Appl Physiol 1987
Apr;62(4):1356-62
19. Jenkins JS, Valcke CP, Ward DS An extended soluble gas exchange model for estimating pulmonary perfusion--II:
Simulation results. IEEE Trans Biomed Eng 1989 Nov;36(11):1105-12
20. Lim LL A statistical model of the VA/Q distribution. J Appl Physiol 1990 Jul;69(1):281-92
21. Vidal Melo MF, Caprihan , A., Luft, UC, Loeppky JA A Distribution of ventilation and diffusion with perfusion in a two-
compartment model of gas exchange.. In "Oxygen Transport to Tissue" (J. Piiper, M. Meyer, and T.K. Gold. Eds.). Vol.XII. pp.
653-664. Plenum, Ney York, 1990.
22. Vidal Melo MF, Loeppky JA, Caprihan A, Luft UC Alveolar ventilation to perfusion heterogeneity and diffusion impairment in a
mathematical model of gas exchange. Comput Biomed Res 1993 Apr;26(2):103-120
23. Marini JJ, Crooke PS A General mathematical model for Respiratory Dynamics Relevant to the Clinical Setting. Am Rev
Repir Dis 1993; Vol 147: 14-24
24. Gavaghan DJ, Hahn CE A tidal breathing model of the forced inspired inert gas sinewave technique. Respir Physiol 1996
Nov;106(2):209-21
25. Sainsbury MC, Lorenzi A, Williams EM, Hahn CE A reconciliation of continuous and tidal ventilation gas exchange models.
Respir Physiol 1997 Aug;109(2):195
26. Busso T, Robbins PA Evaluation of estimates of alveolar gas exchange by using a tidally ventilated nonhomogenous lung
model. J Appl Physiol 1997 Jun;82(6):1963-71
27. Liu CH, Niranjan SC, Clark JW Jr, San KY, Zwischenberger JB, Bidani A Airway mechanics, gas exchange, and blood flow in
a nonlinear model of the normal human lung. J Appl Physiol 1998 Apr;84(4):1447-69
28. Vidal Melo MF ffect of cardiac output on pulmonary gas exchange: role of diffusion limitation with VA/Q mismatch. espir
Physiol 1998 Jul;113(1):23-32
29. Whiteley JP, Gavaghan DJ, Hahn CEW A mathematical evaluation of the multiple breath nitrogen washout(MBNW)
technique and the multiple inert gas elimination technique(MIGET). J Theor Biol 1998 Oct 21;194(4):517-39
30. Whiteley JP, Gavaghan DJ, Hahn CE. A tidal breathing model for the multiple inert gas elimination technique. J Appl Physiol
1999 Jul;87(1):161-9
31. Whiteley JP, Gavaghan DJ, Hahn CE The effect of inspired oxygen concentration on the ventilation-perfusion distribution in
inhomogeneous lungs. J Theor Biol 2000 Jun 21;204(4):575-85
32. Crooke, P.S Hota, S., Marini, J.J. Hotchkiss, J.R. A Mathematical Model for Carbon Disoxide Exchange during Mechanical
Ventilation with Tracheal Gas Insufflation (TGI) by Mathematical and Computer Modelling, 29(1999), pp. 45-61.
33. Niranjan SC, Bidani A, Ghorbel F, Zwischenberger JB, Clark JW Jr. Theoretical Study of inspiratory flow waveforms during
mechanical ventilation on pulmonary blood flow and gas exchange. Comput Biomed Res 1999 Aug; 32(4):355-90
34. Hermgenes de Andrade Jos. Auto perfeccin con Hatha Yoga Buenos Aires, Editorial Kier 1992
35. Andrews. Deep Breathing Can Improve Fitness. The Lancet 1998 May,Volumen 351, number 9112
36. Lucia Spicuzza, Alessandra Gabutti, Cesare Porta, Nicola Montano, Luciano Bernardi. Yoga and chemoreflex response to
hypoxia and hypercapnia The Lancet 2000 October 28, Volumen 356.