An integrated model of the human ventilatory control system: the response to hypoxia

Comparison of intranasal naloxone and intranasal nalmefene in a translational model assessing the impact of synthetic opioid overdose on respiratory depression and cardiac arrest

Introduction

Using a validated translational model that quantitatively predicts opioid-induced respiratory depression and cardiac arrest, we compared cardiac arrest events caused by synthetic opioids (fentanyl, carfentanil) following rescue by intranasal (IN) administration of the μ-opioid receptor antagonists naloxone and nalmefene.

Methods

This translational model was originally developed by Mann et al. (Clin Pharmacol Ther 2022) to evaluate the effectiveness of intramuscular (IM) naloxone. We initially implemented this model using published codes, reproducing the effects reported by Mann et al. on the incidence of cardiac arrest events following intravenous doses of fentanyl and carfentanil as well as the reduction in cardiac arrest events following a standard 2 mg IM dose of naloxone. We then expanded the model in terms of pharmacokinetic and µ-opioid receptor binding parameters to simulate effects of 4 mg naloxone hydrochloride IN and 3 mg nalmefene hydrochloride IN, both FDA-approved for the treatment of opioid overdose. Model simulations were conducted to quantify the percentage of cardiac arrest in 2000 virtual patients in both the presence and absence of IN antagonist treatment.

Results

Following simulated overdoses with both fentanyl and carfentanil in chronic opioid users, IN nalmefene produced a substantially greater reduction in the incidence of cardiac arrest compared to IN naloxone. For example, following a dose of fentanyl (1.63 mg) producing cardiac arrest in 52.1% (95% confidence interval, 47.3-56.8) of simulated patients, IN nalmefene reduced this rate to 2.2% (1.0-3.8) compared to 19.2% (15.5-23.3) for IN naloxone. Nalmefene also produced large and clinically meaningful reductions in the incidence of cardiac arrests in opioid naïve subjects. Across dosing scenarios, simultaneous administration of four doses of IN naloxone were needed to reduce the percentage of cardiac arrest events to levels that approached those produced by a single dose of IN nalmefene.

Conclusion

Simulations using this validated translational model of opioid overdose demonstrate that a single dose of IN nalmefene produces clinically meaningful reductions in the incidence of cardiac arrest compared to IN naloxone following a synthetic opioid overdose. These findings are especially impactful in an era when >90% of all opioid overdose deaths are linked to synthetic opioids such as fentanyl.

Physiologically based modeling reveals different risk of respiratory depression after fentanyl overdose between adults and children

Despite a rapid increase in pediatric mortality rate from prescription and illicit opioids, there is limited research on the dose-dependent impact of opioids on respiratory depression in children, the leading cause of opioid-associated death. In this article, we extend a previously developed translational model to cover pediatric populations by incorporating age-dependent pharmacokinetic, pharmacodynamic, and physiological changes compared to adults. Our model reproduced previous perioperative clinical findings that adults and children have similar risk of respiratory depression at the same plasma fentanyl concentration when specific endpoints (minute ventilation, CO2 tension in the blood) were used. However, our model points to a potential caveat that, in a perioperative setting, routine use of mechanical ventilation and supplemental oxygen maintained the blood and tissue oxygen partial pressures in patients and prevented the use of oxygen-related endpoints to evaluate the consequences of respiratory depression. In a community setting when such oxygenation procedures are not immediately available, our model suggests that the higher oxygen demand and reduced cerebrovascular reactivity could make children more susceptible to severe hypoxemia and brain hypoxia, even with the same plasma fentanyl concentration as adults. Our work indicates that when developing intervention strategies to protect children from opioid overdose in a community setting, these pediatric-specific factors may need to be considered.

Intranasal Naloxone Repeat Dosing Strategies and Fentanyl Overdose: A Simulation-Based Randomized Clinical Trial

Importance

Questions have emerged as to whether standard intranasal naloxone dosing recommendations (ie, 1 dose with readministration every 2-3 minutes if needed) are adequate in the era of illicitly manufactured fentanyl and its derivatives (hereinafter, fentanyl).

Objective

To compare naloxone plasma concentrations between different intranasal naloxone repeat dosing strategies and to estimate their effect on fentanyl overdose.

Design, Setting, and Participants

This unblinded crossover randomized clinical trial was conducted with healthy participants in a clinical pharmacology unit (Spaulding Clinical Research, West Bend, Wisconsin) in March 2021. Inclusion criteria included age 18 to 55 years, nonsmoking status, and negative test results for the presence of alcohol or drugs of abuse. Data analysis was performed from October 2021 to May 2023.

Intervention

Naloxone administered as 1 dose (4 mg/0.1 mL) at 0, 2.5, 5, and 7.5 minutes (test), 2 doses at 0 and 2.5 minutes (test), and 1 dose at 0 and 2.5 minutes (reference).

Main Outcomes and Measures

The primary outcome was the first prespecified time with higher naloxone plasma concentration. The secondary outcome was estimated brain hypoxia time following simulated fentanyl overdoses using a physiologic pharmacokinetic-pharmacodynamic model. Naloxone concentrations were compared using paired tests at 3 prespecified times across the 3 groups, and simulation results were summarized using descriptive statistics.

Results

This study included 21 participants, and 18 (86%) completed the trial. The median participant age was 34 years (IQR, 27-50 years), and slightly more than half of participants were men (11 [52%]). Compared with 1 naloxone dose at 0 and 2.5 minutes, 1 dose at 0, 2.5, 5, and 7.5 minutes significantly increased naloxone plasma concentration at 10 minutes (7.95 vs 4.42 ng/mL; geometric mean ratio, 1.95 [1-sided 97.8% CI, 1.28-∞]), whereas 2 doses at 0 and 2.5 minutes significantly increased the plasma concentration at 4.5 minutes (2.24 vs 1.23 ng/mL; geometric mean ratio, 1.98 [1-sided 97.8% CI, 1.03-∞]). No drug-related serious adverse events were reported. The median brain hypoxia time after a simulated fentanyl 2.97-mg intravenous bolus was 4.5 minutes (IQR, 2.1-∞ minutes) with 1 naloxone dose at 0 and 2.5 minutes, 4.5 minutes (IQR, 2.1-∞ minutes) with 1 naloxone dose at 0, 2.5, 5, and 7.5 minutes, and 3.7 minutes (IQR, 1.5-∞ minutes) with 2 naloxone doses at 0 and 2.5 minutes.

Conclusions and Relevance

In this clinical trial with healthy participants, compared with 1 intranasal naloxone dose administered at 0 and 2.5 minutes, 1 dose at 0, 2.5, 5, and 7.5 minutes significantly increased naloxone plasma concentration at 10 minutes, whereas 2 doses at 0 and 2.5 minutes significantly increased naloxone plasma concentration at 4.5 minutes. Additional research is needed to determine optimal naloxone dosing in the community setting.

Trial Registration

ClinicalTrials.gov Identifier: NCT04764630.

Prediction of exhaled carbon dioxide concentration using a computer-simulated person that included alveolar gas exchange

Accurate prediction of inhaled CO₂ concentration and alveolar gas exchange efficiency would improve the prediction of CO₂ concentrations around the human body, which is essential for advanced ventilation design in buildings. We therefore, developed a computer-simulated person (CSP) that included a computational fluid dynamics approach. The CSP simulates metabolic heat production at the skin surface and carbon dioxide (CO₂ ) gas exchange at the alveoli during the transient breathing cycle. This makes it possible to predict the CO₂ distribution around the human body. The numerical model of the CO₂ gas exchange mechanism includes both the upper and lower airways and makes it possible to calculate the alveolar CO₂ partial pressure; this improves the prediction accuracy. We used the CSP to predict emission rates of metabolically generated CO₂ exhaled by a person and assumed that the tidal volume will be unconsciously reduced as a result of exposure to poor indoor air quality. A reduction in tidal volume resulted in a decrease in CO₂ emission rates of the same magnitude as was observed in our published experimental data. We also observed that the predicted inhaled CO₂ concentration depended on the flow pattern around the human body, as would be expected.

Development of a Translational Model to Assess the Impact of Opioid Overdose and Naloxone Dosing on Respiratory Depression and Cardiac Arrest

In response to a surge of deaths from synthetic opioid overdoses, there have been increased efforts to distribute naloxone products in community settings. Prior research has assessed the effectiveness of naloxone in the hospital setting; however, it is challenging to assess naloxone dosing regimens in the community/first-responder setting, including reversal of respiratory depression effects of fentanyl and its derivatives (fentanyls). Here, we describe the development and validation of a mechanistic model that combines opioid mu receptor binding kinetics, opioid agonist and antagonist pharmacokinetics, and human respiratory and circulatory physiology, to evaluate naloxone dosing to reverse respiratory depression. Validation supports our model, which can quantitatively predict displacement of opioids by naloxone from opioid mu receptors in vitro, hypoxia-induced cardiac arrest in vivo, and opioid-induced respiratory depression in humans from different fentanyls. After validation, overdose simulations were performed with fentanyl and carfentanil followed by administration of different intramuscular naloxone products. Carfentanil induced more cardiac arrest events and was more difficult to reverse than fentanyl. Opioid receptor binding data indicated that carfentanil has substantially slower dissociation kinetics from the opioid receptor compared to 9 other fentanyls tested, which likely contributes to the difficulty in reversing carfentanil. Administration of the same dose of naloxone intramuscularly from 2 different naloxone products with different formulations resulted in differences in the number of virtual patients experiencing cardiac arrest. This work provides a robust framework to evaluate dosing regimens of opioid receptor antagonists to reverse opioid-induced respiratory depression, including those caused by newly emerging synthetic opioids.

MATHEMATICAL MODELS OF HUMAN RESPIRATORY AND BLOOD CIRCULATORY SYSTEMS

Aim. To analyze modern approaches to mathematical modeling of human respiratory and blood circulatory systems. Methods. Comprehensive review of scientific literature sources extracted from domestic and international resources databases. Results. Historical information and modern data concerning mathematical modeling of human functional respiratory and blood circulatory systems were summarized and analyzed in present ¬review; current trends in approaches to the construction of these models were revealed. Conclusions. Currently, two main approaches to the mathematical modeling of respiratory and blood circulatory systems exist. One of them is the construction of models of the mechanics of respiration and blood circulation. They are based on the models of mechanics of solid deformable body, thermomechanics, hydromechanics, and continuum mechanics. This approach uses complex mathematical apparatus, including Navier-Stokes equation, which makes it possible to obtain a number of theoretical results, but it is hardly usable for real problems solutions at present time. The second approach is based on the model of F. Grodins, who represented the process of breathing as a controlled dynamic system, described by ordinary differential equations, in which the process control is carried out according to the feedback principle. There is a significant number of modifications of this model, which made it possible to simulate various disturbing influences, such as physical activity, hypoxia and hyperemia, and to predict parameters characterizing functional respiratory system under these disturbing influences.

Seven Mathematical Models of Hemorrhagic Shock

Although mathematical modelling of pressure-flow dynamics in the cardiocirculatory system has a lengthy history, readily finding the appropriate model for the experimental situation at hand is often a challenge in and of itself. An ideal model would be relatively easy to use and reliable, besides being ethically acceptable. Furthermore, it would address the pathogenic features of the cardiovascular disease that one seeks to investigate. No universally valid model has been identified, even though a host of models have been developed. The object of this review is to describe several of the most relevant mathematical models of the cardiovascular system: the physiological features of circulatory dynamics are explained, and their mathematical formulations are compared. The focus is on the whole-body scale mathematical models that portray the subject's responses to hypovolemic shock. The models contained in this review differ from one another, both in the mathematical methodology adopted and in the physiological or pathological aspects described. Each model, in fact, mimics different aspects of cardiocirculatory physiology and pathophysiology to varying degrees: some of these models are geared to better understand the mechanisms of vascular hemodynamics, whereas others focus more on disease states so as to develop therapeutic standards of care or to test novel approaches. We will elucidate key issues involved in the modeling of cardiovascular system and its control by reviewing seven of these models developed to address these specific purposes.

Ventilatory oscillations at exercise in hypoxia: A mathematical model

We evaluated the mechanisms responsible for the instability of ventilation control system under simultaneous metabolic (exercise) and environmental (hypoxia) stresses, promoting the genesis of periodic breathing. A model following the main concepts of ventilatory control has been tested, including cardiovascular and respiratory parameters, characteristics of peripheral and central chemoreceptors, at mild exercise in hypoxia (FIO₂=0.145). Interaction between O₂ and CO₂ sensing was introduced following three different modalities. A sensitivity and multivariate regression analyses closely matched with physiological data for magnitude and period of oscillations. Low FIO₂ and long circulatory delay from lungs to peripheral chemoreceptors (DeltaTp) lengthen the period of oscillations, while high peripheral and central chemoresponses to O₂ and CO₂, low FIO₂ and high DeltaTp increased their magnitude. Peripheral and central O₂/CO₂ interactions highlight the role of CO₂ on peripheral gain to O₂ and the contribution of peripheral afferences on central gain to CO₂. Our model supports the key role of peripheral chemoreceptors in the genesis of ventilatory oscillations. Differences in the dynamics of central and peripheral components might be determinant for the system stability.

Development of a cardiopulmonary mathematical model incorporating a baro-chemoreceptor reflex control system

This article describes the development of a comprehensive mathematical model of the human cardiopulmonary system that combines the respiratory and cardiovascular systems and their associated autonomous nervous control actions. The model is structured to allow the complex interactions between the two systems and the responses of the combined system to be predicted under different physiological conditions. The cardiovascular system model contains 13 compartments, including the heart chambers operating as a pump and the blood vessels represented as distensible tubes configured in a serial and parallel arrangement. The accurate representation of the hemodynamics in the system and the good fit to published pressure and flow waveforms gave confidence in the modelling approach adopted for the cardiovascular system prior to the incorporation of the baroreflex control and the respiratory models. An improved baroreceptor reflex model is developed in this research, incorporating afferent, central and efferent compartments. A sigmoid function is included in the efferent compartment to produce sympathetic and parasympathetic nerve outflow to the effector sites. The baroreflex action is modelled using physiological data, its interaction with the chemoreflex control is explained and the simulation results presented show the ability of the model to predict the static and dynamic hemodynamic responses to environmental disturbances. A previously published respiratory model that includes the mechanics of breathing, gas exchange process and the regulation of the system is then combined with the cardiovascular model to form the cardiopulmonary model. Through comparison with published data, the cardiopulmonary model with the baro-chemoreflex control is validated during hypoxia and hypercapnia. The percentage difference between the predicted and measured changes in the heart rates and the mean arterial pressures are within 3% in both cases. The total peripheral resistance correlates well for hypoxia but is less good for hypercapnia, where the predicted change from normal condition is around 7% compared with a measured change of 23%. An example showing the application of the proposed model in sport science is also included.

Systemic modelling of human bioenergetics and blood circulation

This work reviews the main aspects of human bioenergetics and the dynamics of the cardiovascular system, with emphasis on modelling their physiological characteristics. The methods used to study human bioenergetics and circulation dynamics, including the use of mathematical models, are summarised. The main characteristics of human bioenergetics, including mitochondrial metabolism and global energy balance, are first described, and the systemic aspects of blood circulation and related physiological issues are introduced. The authors also discuss the present status of studies of human bioenergetics and blood circulation. Then, the limitations of the existing studies are described in an effort to identify directions for future research towards integrated and comprehensive modelling. This review emphasises that a multi-scale and multi-physical approach to bioenergetics and blood circulation that considers multiple scales and physiological factors are necessary for the appropriate clinical application of computational models.

Computational models for the study of heart-lung interactions in mammals

The operation and regulation of the lungs and the heart are closely related. This is evident when examining the anatomy within the thorax cavity, in the brainstem and in the aortic and carotid arteries where chemoreceptors and baroreceptors, which provide feedback affecting the regulation of both organs, are concentrated. This is also evident in phenomena such as respiratory sinus arrhythmia where the heart rate increases during inspiration and decreases during expiration, in other types of synchronization between the heart and the lungs known as cardioventilatory coupling and in the association between heart failure and sleep apnea where breathing is interrupted periodically by periods of no-breathing. The full implication and physiological significance of the cardiorespiratory coupling under normal, pathological, or extreme physiological conditions are still unknown and are subject to ongoing investigation both experimentally and theoretically using mathematical models. This article reviews mathematical models that take heart-lung interactions into account. The main ideas behind low dimensional, phenomenological models for the study of the heart-lung synchronization and sleep apnea are described first. Higher dimensions, physiology-based models are described next. These models can vary widely in detail and scope and are characterized by the way the heart-lung interaction is taken into account: via gas exchange, via the central nervous system, via the mechanical interactions, and via time delays. The article emphasizes the need for the integration of the different sources of heart-lung coupling as well as the different mathematical approaches.

Prediction of the susceptibility to AMS in simulated altitude

Acute mountain sickness (AMS) develops when rapidly ascending to high altitudes. However, some mountaineers will suffer from AMS even at 2,000 m and others not until 5,000 m. The awareness of the individual susceptibility for AMS would be helpful for preventive strategies. Thus, the main purpose of this paper is the comparison of existing studies dealing with the prediction of AMS susceptibility and to draw conclusions on presently most valuable tests.

DATA SOURCE

A PubMed search has been performed, and preliminary observations from our laboratory have been included. The cautious conclusion derived from the reviewed 16 studies is that values of arterial oxygen saturation (SaO(2)), determined 20-30 min after exposure to simulated hypoxia equivalent to 2,300-4,200 m, seem to be the most useful predictors of AMS susceptibility (>80% correct prediction). Because the sympathetic activation during acute exposure to hypoxia may well contribute to the AMS development, parameters like heart rate variability or blood lactate could even enhance this predictability. The ventilatory response to hypoxia is easily trainable by pre-exposures to hypoxia but considers only part of the complex acclimatization process.

A model for control of breathing in mammals: coupling neural dynamics to peripheral gas exchange and transport

A new model for aspects of the control of respiration in mammals has been developed. The model integrates a reduced representation of the brainstem respiratory neural controller together with peripheral gas exchange and transport mechanisms. The neural controller consists of two components. One component represents the inspiratory oscillator in the pre-Bötzinger complex (pre-BötC) incorporating biophysical mechanisms for rhythm generation. The other component represents the ventral respiratory group (VRG), which is driven by the pre-BötC for generation of inspiratory (pre)motor output. The neural model was coupled to simplified models of the lungs incorporating oxygen and carbon dioxide transport. The simplified representation of the brainstem neural circuitry has regulation of both frequency and amplitude of respiration and is done in response to partial pressures of oxygen and carbon dioxide in the blood using proportional (P) and proportional plus integral (PI) controllers. We have studied the coupled system under open and closed loop control. We show that two breathing regimes can exist in the model. In one regime an increase in the inspiratory frequency is accompanied by an increase in amplitude. In the second regime an increase in frequency is accompanied by a decrease in amplitude. The dynamic response of the model to changes in the concentration of inspired O2 or inspired CO2 was compared qualitatively with experimental data reported in the physiological literature. We show that the dynamic response with a PI-controller fits the experimental data better but suggests that when high levels of CO2 are inspired the respiratory system cannot reach steady state. Our model also predicts that there could be two possible mechanisms for apnea appearance when 100% O2 is inspired following a period of 5% inspired O2. This paper represents a novel attempt to link neural control and gas transport mechanisms, highlights important issues in amplitude and frequency control and sets the stage for more complete neurophysiological control models.

A Mathematical Model of Human Respiration at Altitude

We developed a mathematical model of human respiration in the awake state that can be used to predict changes in ventilation, blood gases, and other critical variables during conditions of hypocapnia, hypercapnia and these conditions combined with hypoxia. Hence, the model is capable of describing ventilation changes due to the hypocapnic-hypoxia of high altitude. The basic model is that of Grodins et al. [Grodins, F. S., J. Buell, and A. J. Bart. J. Appl. Physiol. 22:260-276, 1967]. We updated the descriptions of (1) the effects of blood gases on cardiac output and cerebral blood flow, (2) acid-base balance in blood and tissues, (3) O2 and CO2 binding to hemoglobin and most importantly, (4) the respiratory-chemostat controller. The controller consists of central and peripheral sections. The central chemoceptor-induced ventilation response is simply a linear function of brain P(CO2) above a threshold value. The peripheral response has both a linear term similar to that for the central chemoceptors, but dependent upon carotid body P(CO2) and with a different threshold and a complex, nonlinear term that includes multiplication of separate terms involving carotid body P(O2) and P(CO2). Together, these terms produce 'dogleg'-shaped curves of ventilation plotted against P(CO2) which form a fan-like family for different values of P(CO2). With this chemical controller, our model closely describes a wide range of experimental data under conditions of solely changes in P(CO2) and for short-term hypoxia coupled with P(CO2) changes. This model can be used to accurately describe changes in ventilation and respiratory gases during ascent and during short-term residence at altitude. Hence, it has great applicability to studying O2-delivery systems in aircraft.

Multi-organ system model of O2 and CO2 transport during isocapnic and poikilocapnic hypoxia

A multi-organ systems model of O(2) and CO(2) transport is developed to analyze the control of ventilation and blood flow during hypoxia. Among the aspects of the control processes that this model addressed are possible mechanisms responsible for the second phase of the ventilatory hypoxic response to mild hypoxia, i.e., hypoxic ventilatory decline (HVD). Species mass transport processes are described by compartmental mass balances in brain, heart, skeletal muscle, and "other tissues" connected in parallel via the circulation. In pulmonary and systemic capillaries and in the vasculature connecting the systemic tissues, species transport processes are represented by a one-dimensional, convection-dispersion model. The effects of bicarbonate acid-base buffering, hemoglobin, and myoglobin on the transport processes are included. The model incorporates feedback control mechanisms through a cardiorespiratory control system in which peripheral and central chemoreceptors sense O(2) and CO(2) partial pressures. Model simulations of the ventilatory responses to isocapnic and poikilocapnic hypoxia show two phases with distinct dynamics. A fast phase is discernable immediately after switching from normoxic to hypoxic conditions, while a delayed slow phase (HVD) typically becomes manifested after several minutes. Model simulations allow quantitative evaluation of several proposed mechanisms to account for HVD. Under isocapnic hypoxia, simulations indicate that an increase in brain blood flow has no effect on HVD, but that HVD can be entirely described by central ventilatory depression (CVD). Under poikilocapnic hypoxia, the hypocapnia caused by hypoxic hyperventilation has no effect on HVD.

Model-Based Analysis of Mechanisms Responsible for Sleep-Induced Carbon Dioxide Differences

This work describes a comprehensive mathematical model of the human respiratory control system which incorporates the central mechanisms for predicting sleep-induced changes in chemical regulation of ventilation. The model integrates four individual compartments for gas storage and exchange, namely alveolar air, pulmonary blood, tissue capillary blood, body tissues, and gas transport between them. An essential mechanism in the carbon dioxide transport is its dissociation into bicarbonate and acid, where a buffering mechanism through hemoglobin is used to prevent harmfully low pH levels. In the current model, we assume high oxygen levels and consider intracellular hydrogen ion concentration as the principal respiratory control variable. The resulting system of delayed differential equations is solved numerically. With an appropriate choice of key parameters, such as velocity of blood flow and gain of a non-linear controller function, the model provides steady-state results consistent with our experimental observations measured in subjects across sleep onset. Dynamic predictions from the model give new insights into the behaviour of the system in subjects with different buffering capacities and suggest novel hypotheses for future experimental and clinical studies.

Mathematical modeling of cardiovascular system dynamics using a lumped parameter method

This work reviews the main aspects of cardiovascular system dynamics with an emphasis on modeling hemodynamic characteristics by the use of a lumped parameter approach. The methodological and physiological aspects of the circulation dynamics are summarized with the help of existing mathematical models. The main characteristics of the hemodynamic elements, such as the heart and arterial and venous systems, are first described. Distributed models of an arterial network are introduced, and their characteristics are compared with those of lumped parameter models. We also discuss the nonlinear characteristics of the pressure-volume relationship in veins. Then the control pathways that participate in feedback mechanisms (baroreceptors and cardiopulmonary receptors) are described to explain the interaction between hemodynamics and autonomic nerve control in the circulation. Based on a set-point model, the computational aspects of reflex control are explained.

Specific carotid body chemostimulation is sufficient to elicit phrenic poststimulus frequency decline in a novel in situ dual-perfused rat preparation

Time-dependent ventilatory responses to hypoxic and hypercapnic challenges, such as posthypoxic frequency decline (PHxFD) and posthypercapnic frequency decline (PHcFD), could profoundly affect breathing stability. However, little is known about the mechanisms that mediate these phenomena. To determine the contribution of specific carotid body chemostimuli to PHxFD and PHcFD, we developed a novel in situ arterially perfused, vagotomized, decerebrate rat preparation in which central and peripheral chemoreceptors are perfused separately (i.e., a nonanesthetized in situ dual perfused preparation). We confirmed that 1) the perfusion of central and peripheral chemoreceptor compartments was independent by applying specific carotid body hypoxia and hypercapnia before and after carotid sinus nerve transection, 2) the PCO(2) chemoresponse of the dual perfused preparation was similar to other decerebrate preparations, and 3) the phrenic output was stable enough to allow investigation of time-dependent phenomena. We then applied four 5-min bouts (separated by 5 min) of specific carotid body hypoxia (40 Torr PO(2) and 40 Torr PCO(2)) or hypercapnia (100 Torr PO(2) and 60 Torr PCO(2)) while holding the brain stem PO(2) and PCO(2) constant. We report the novel finding that specific carotid body chemostimuli were sufficient to elicit several phrenic time-dependent phenomena in the rat. Hypoxic challenges elicited PHxFD that increased with bout, leading to progressive augmentation of the phrenic response. Conversely, hypercapnia elicited short-term depression and PHcFD, neither of which was bout dependent. These results, placed in the context of previous findings, suggest multiple physiological mechanisms are responsible for PHxFD and PHcFD, a redundancy that may illustrate that these phenomena have significant adaptive advantages.

A mathematical model of ventilation response to inhaled carbon monoxide

A comprehensive mathematical model, describing the respiration, circulation, oxygen metabolism, and ventilatory control, is assembled for the purpose of predicting acute ventilation changes from exposure to carbon monoxide in both humans and animals. This Dynamic Physiological Model is based on previously published work, reformulated, extended, and combined into a single model. Model parameters are determined from literature values, fitted to experimental data, or allometrically scaled between species. The model predictions are compared with ventilation-time history data collected in goats exposed to carbon monoxide, with quantitatively good agreement. The model reaffirms the role of brain hypoxia on hyperventilation during carbon monoxide exposures. Improvement in the estimation of total ventilation, through a more complete knowledge of ventilation control mechanisms and validated by animal data, will increase the accuracy of inhalation toxicology estimates.

Prediction of susceptibility to acute mountain sickness by SaO2 values during short-term exposure to hypoxia

Prediction of the development of acute mountain sickness (AMS) in individuals going to high altitudes is still a matter of debate. Whereas some studies found that subjects with a blunted hypoxic ventilatory response (HVR) are predisposed to AMS, others did not. However, the HVR has often been determined under very acute (5 to 10 min) isocapnic hypoxia without consideration of the subsequent hypoxic ventilatory decline (HVD), and the assessment of AMS susceptibility was based on a single altitude exposure. Therefore, the aim of the present study was to evaluate the relationship between the individual arterial oxygen saturation (Sa(O2)) after a 20- to 30-min exposure to poikilocapnic hypoxia and the AMS susceptibility based on repeated observations. A total of 150 healthy male and female mountaineers (ages: 42 +/- 13 yr), 63 of whom had known susceptibility to AMS and 87 of whom never suffered from AMS, were exposed to various degrees of normobaric and hypobaric hypoxia. Sa(O2) values were taken by finger pulseoximetry after 20 to 30 min of hypoxic exposure. Sa(O2) values after 20 to 30 min of hypoxia were on average 4.9% lower in subjects susceptible to AMS than in those who were not. Logistic regression analysis revealed altitude-dependent Sa(O2) values to be predictive for AMS susceptibility. Based on the derived model, AMS susceptibility was correctly predicted in 86% of the selected individuals exposed to short-term hypoxia. In conclusion, Sa(O2) values after 20 to 30 min of exposure to normobaric or hypobaric hypoxia represent a useful tool to detect subjects highly susceptible to AMS.

Cerebral autoregulation and gas exchange studied using a human cardiopulmonary model

The goal of this work is to study the cerebral autoregulation, brain gas exchange, and their interaction by means of a mathematical model. We have previously developed a model of the human cardiopulmonary (CP) system, which included the whole body circulatory system, lung and peripheral tissue gas exchange, and the central nervous system control of arterial pressure and ventilation. In this study, we added a more detailed description of cerebral circulation, cerebrospinal fluid (CSF) dynamics, brain gas exchange, and cerebral blood flow (CBF) autoregulation. Two CBF regulatory mechanisms are included: autoregulation and CO(2) reactivity. Central chemoreceptor control of ventilation is also included. We first established nominal operating conditions for the cerebral model in an open-loop configuration using data generated by the CP model as inputs. The cerebral model was then integrated into the larger CP model to form a new integrated CP model, which was subsequently used to study cerebral hemodynamic and gas exchange responses to test protocols commonly used in the assessment of CBF autoregulation (e.g., carotid artery compression and the thigh-cuff deflation test). The model can closely mimic the experimental findings and provide biophysically based insights into the dynamics of cerebral autoregulation and brain tissue gas exchange as well as the mechanisms of their interaction during test protocols, which are aimed at assessing the degree of autoregulation. With further refinement, our CP model may be used on measured data associated with the clinical evaluation of the cerebral autoregulation and brain oxygenation in patients.

An integrated model of the human ventilatory control system: the response to hypercapnia

This work presents a mathematical model of the human respiratory control system, based on physiological knowledge. It includes three compartments for gas storage and exchange (lungs, brain tissue and other body tissues), and various kinds of feedback mechanisms. These comprehend peripheral chemoreceptors in the carotid body, central chemoreceptors in the medulla and a central ventilatory depression. The latter acts by reducing the response of the central neural system to the afferent peripheral chemoreceptor activity during prolonged hypoxia of the brain tissue. Furthermore, the model considers local blood flow adjustments in response to O2 and CO2 arterial pressure changes. In this study, the model has been validated by simulating the response to square changes in alveolar PCO2, performed at different constant levels of alveolar PO2. A good agreement with data reported in the literature has been checked. Subsequently, a sensitivity analysis on the role of the main feedback mechanisms on ventilation response to CO2 has been performed. The results suggest that the ventilatory response to CO2 challenges during hyperoxia can be almost completely ascribed to the central chemoreflex, while, during normoxia, the peripheral chemoreceptors provide a modest contribution too. By contrast, the response to hypercapnic stimuli during hypoxia involves a complex superimposition among different factors with disparate dynamics. Hence, results suggest that the ventilatory response to hypercapnia during hypoxia is more complex than that provided by simple empirical models, and that discrimination between the central and peripheral components based on time constants may be misleading.