Effects of neural drives on breathing in the awake state in humans

Effects of ion channel noise on neural circuits: an application to the respiratory pattern generator to investigate breathing variability

Neural activity generally displays irregular firing patterns even in circuits with apparently regular outputs, such as motor pattern generators, in which the output frequency fluctuates randomly around a mean value. This "circuit noise" is inherited from the random firing of single neurons, which emerges from stochastic ion channel gating (channel noise), spontaneous neurotransmitter release, and its diffusion and binding to synaptic receptors. Here we demonstrate how to expand conductance-based network models that are originally deterministic to include realistic, physiological noise, focusing on stochastic ion channel gating. We illustrate this procedure with a well-established conductance-based model of the respiratory pattern generator, which allows us to investigate how channel noise affects neural dynamics at the circuit level and, in particular, to understand the relationship between the respiratory pattern and its breath-to-breath variability. We show that as the channel number increases, the duration of inspiration and expiration varies, and so does the coefficient of variation of the breath-to-breath interval, which attains a minimum when the mean duration of expiration slightly exceeds that of inspiration. For small channel numbers, the variability of the expiratory phase dominates over that of the inspiratory phase, and vice versa for large channel numbers. Among the four different cell types in the respiratory pattern generator, pacemaker cells exhibit the highest sensitivity to channel noise. The model shows that suppressing input from the pons leads to longer inspiratory phases, a reduction in breathing frequency, and larger breath-to-breath variability, whereas enhanced input from the raphe nucleus increases breathing frequency without changing its pattern.

NEW & NOTEWORTHY

A major source of noise in neuronal circuits is the "flickering" of ion currents passing through the neurons' membranes (channel noise), which cannot be suppressed experimentally. Computational simulations are therefore the best way to investigate the effects of this physiological noise by manipulating its level at will. We investigate the role of noise in the respiratory pattern generator and show that endogenous, breath-to-breath variability is tightly linked to the respiratory pattern.

Manipulation of central blood volume and implications for respiratory control function

The respiratory operating point (ventilatory or arterial Pco2 response) is determined by the intersection point between the controller and plant subsystem elements within the respiratory control system. However, to what extent changes in central blood volume (CBV) influence these two elements and the corresponding implications for the respiratory operating point remain unclear. To examine this, 17 apparently healthy male participants were exposed to water immersion (WI) or lower body negative pressure (LBNP) challenges to manipulate CBV and determine the corresponding changes. The respiratory controller was characterized by determining the linear relationship between end-tidal Pco2 (PetCO2) and minute ventilation (V̇e) [V̇e = S × (PetCO2 − B)], whereas the plant was determined by the hyperbolic relationship between V̇e and PetCO2 (PetCO2 = A/V̇e + C). Changes in V̇e at the operating point were not observed under either WI or LBNP conditions despite altered PetCO2 ( P < 0.01), indicating a moving respiratory operating point. An increase (WI) and a decrease (LBNP) in CBV were shown to reset the controller element (PetCO2 intercept B) rightward and leftward, respectively ( P < 0.05), without any change in S, whereas the plant curve remained unaltered at the operating point. Collectively, these findings indicate that modification of the controller element rather than the plant element is the major factor that contributes toward an alteration of the respiratory operating point during CBV shifts.

Sleep disordered breathing in chronic spinal cord injury

STUDY OBJECTIVES

Spinal cord injury (SCI) is associated with 2-5 times greater prevalence of sleep disordered breathing (SDB) than the general population. The contribution of SCI on sleep and breathing at different levels of injury using two scoring methods has not been assessed. The objectives of this study were to characterize the sleep disturbances in the SCI population and the associated physiological abnormalities using quantitative polysomnography and to determine the contribution of SCI level on the SDB mechanism.

METHODS

We studied 26 consecutive patients with SCI (8 females; age 42.5 ± 15.5 years; BMI 25.9 ± 4.9 kg/m2; 15 cervical and 11 thoracic levels) by spirometry, a battery of questionnaires and by attended polysomnography with flow and pharyngeal pressure measurements. Inclusion criteria for SCI: chronic SCI (> 6 months post injury), level T6 and above and not on mechanical ventilation. Ventilation, end-tidal CO2 (PETCO2), variability in minute ventilation (VI-CV) and upper airway resistance (RUA) were monitored during wakefulness and NREM sleep in all subjects. Each subject completed brief history and exam, Epworth Sleepiness Scale (ESS), Pittsburgh Sleep Quality Index (PSQI), Berlin questionnaire (BQ) and fatigue severity scale (FSS). Sleep studies were scored twice, first using standard 2007 American Academy of Sleep Medicine (AASM) criteria and second using new 2012 recommended AASM criteria.

RESULTS

Mean PSQI was increased to 10.3 ± 3.7 in SCI patients and 92% had poor sleep quality. Mean ESS was increased 10.4 ± 4.4 in SCI patients and excessive daytime sleepiness (ESS ≥ 10) was present in 59% of the patients. Daytime fatigue (FSS > 20) was reported in 96% of SCI, while only 46% had high-risk score of SDB on BQ. Forced vital capacity (FVC) in SCI was reduced to 70.5% predicted in supine compared to 78.5% predicted in upright positions (p < 0.05). Likewise forced expiratory volume in first second (FEV1) was 64.9% predicted in supine compared to 74.7% predicted in upright positions (p < 0.05). Mean AHI in SCI patients was 29.3 ± 25.0 vs. 20.0 ± 22.8 events/h using the new and conventional AASM scoring criteria, respectively (p < 0.001). SCI patients had SDB (AHI > 5 events/h) in 77% of the cases using the new AASM scoring criteria compared to 65% using standard conventional criteria (p < 0.05). In cervical SCI, VI decreased from 7.2 ± 1.6 to 5.5 ± 1.3 L/min, whereas PETCO2 and VI-CV, increased during sleep compared to thoracic SCI.

CONCLUSION

The majority of SCI survivors have symptomatic SDB and poor sleep that may be missed if not carefully assessed. Decreased VI and increased PETCO2 during sleep in patients with cervical SCI relative to thoracic SCI suggests that sleep related hypoventilation may contribute to the pathogenesis SDB in patients with chronic cervical SCI.

Differences in the control of breathing between Himalayan and sea-level residents

We compared the control of breathing of 12 male Himalayan highlanders with that of 21 male sea-level Caucasian lowlanders using isoxic hyperoxic ( = 150 mmHg) and hypoxic ( = 50 mmHg) Duffin's rebreathing tests. Highlanders had lower mean +/- s.e.m. ventilatory sensitivities to CO(2) than lowlanders at both isoxic tensions (hyperoxic: 2.3 +/- 0.3 vs. 4.2 +/- 0.3 l min(1) mmHg(1), P = 0.021; hypoxic: 2.8 +/- 0.3 vs. 7.1 +/- 0.6 l min(1) mmHg(1), P < 0.001), and the usual increase in ventilatory sensitivity to CO(2) induced by hypoxia in lowlanders was absent in highlanders (P = 0.361). Furthermore, the ventilatory recruitment threshold (VRT) CO(2) tensions in highlanders were lower than in lowlanders (hyperoxic: 33.8 +/- 0.9 vs. 48.9 +/- 0.7 mmHg, P < 0.001; hypoxic: 31.2 +/- 1.1 vs. 44.7 +/- 0.7 mmHg, P < 0.001). Both groups had reduced ventilatory recruitment thresholds with hypoxia (P < 0.001) and there were no differences in the sub-threshold ventilations (non-chemoreflex drives to breathe) between lowlanders and highlanders at both isoxic tensions (P = 0.982), with a trend for higher basal ventilation during hypoxia (P = 0.052). We conclude that control of breathing in Himalayan highlanders is distinctly different from that of sea-level lowlanders. Specifically, Himalayan highlanders have decreased central and absent peripheral sensitivities to CO(2). Their response to hypoxia was heterogeneous, with the majority decreasing their VRT indicating either a CO(2)-independent increase in activity of peripheral chemoreceptor or hypoxia-induced increase in [H(+)] at the central chemoreceptor. In some highlanders, the decrease in VRT was accompanied by an increase in sensitivity to CO(2), while in others VRT remained unchanged and their sub-threshold ventilations increased, although these were not statistically significant.

Influence of arousal threshold and depth of sleep on respiratory stability in man: analysis using a mathematical model

We examined the effect of arousals (shifts from sleep to wakefulness) on breathing during sleep using a mathematical model. The model consisted of a description of the fluid dynamics and mechanical properties of the upper airways and lungs, as well as a controller sensitive to arterial and brain changes in CO(2), changes in arterial oxygen, and a neural input, alertness. The body was divided into multiple gas store compartments connected by the circulation. Cardiac output was constant, and cerebral blood flows were sensitive to changes in O(2) and CO(2) levels. Arousal was considered to occur instantaneously when afferent respiratory chemical and neural stimulation reached a threshold value, while sleep occurred when stimulation fell below that value. In the case of rigid and nearly incompressible upper airways, lowering arousal threshold decreased the stability of breathing and led to the occurrence of repeated apnoeas. In more compressible upper airways, to maintain stability, increasing arousal thresholds and decreasing elasticity were linked approximately linearly, until at low elastances arousal thresholds had no effect on stability. Increased controller gain promoted instability. The architecture of apnoeas during unstable sleep changed with the arousal threshold and decreases in elasticity. With rigid airways, apnoeas were central. With lower elastances, apnoeas were mixed even with higher arousal thresholds. With very low elastances and still higher arousal thresholds, sleep consisted totally of obstructed apnoeas. Cycle lengths shortened as the sleep architecture changed from mixed apnoeas to total obstruction. Deeper sleep also tended to promote instability by increasing plant gain. These instabilities could be countered by arousal threshold increases which were tied to deeper sleep or accumulated aroused time, or by decreased controller gains.

A negative interaction between brainstem and peripheral respiratory chemoreceptors modulates peripheral chemoreflex magnitude

Interaction between central (brainstem) and peripheral (carotid body) respiratory chemosensitivity is vital to protect blood gases against potentially deleterious fluctuations, especially during sleep. Previously, using an in situ arterially perfused, vagotomized, decerebrate preparation in which brainstem and peripheral chemoreceptors are perfused separately (i.e. dual perfused preparation; DPP), we observed that the phrenic response to specific carotid body hypoxia was larger when the brainstem was held at 25 Torr P(CO(2)) compared to 50 Torr P(CO(2)). This suggests a negative (i.e. hypo-additive) interaction between chemoreceptors. The current study was designed to (a) determine whether this observation could be generalized to all carotid body stimuli, and (b) exclude the possibility that the hypo-additive response was the simple consequence of ventilatory saturation at high brainstem P(CO(2)). Specifically, we tested how steady-state brainstem P(CO(2)) modulates peripheral chemoreflex magnitude in response to carotid body P(CO(2)) and P(O(2)) perturbations, both above and below eupnoeic levels. We found that the peripheral chemoreflex was more responsive the lower the brainstem P(CO(2)) regardless of whether the peripheral chemoreceptors received stimuli which increased or decreased activation. These findings demonstrate a negative interaction between brainstem and peripheral chemosensitivity in the rat in the absence of ventilatory saturation. We suggest that a negative interaction in humans may contribute to increased controller gain associated with sleep-related breathing disorders and propose that the assumption of simple addition between chemoreceptor inputs used in current models of the respiratory control system be reconsidered.

Role of respiratory control mechanisms in the pathogenesis of obstructive sleep disorders

Obstructive sleep disorders develop when the normal reduction in pharyngeal dilator activity at sleep onset occurs in an individual whose pharynx requires a relatively high level of dilator activity to remain sufficiently open. They range from steady snoring, to slowly evolving hypopneas, to fast-recurring obstructive hypopneas and apneas. A fundamental observation is that the polysomnographic picture differs substantially among subjects with the same pharyngeal collapsibility, and even in the same subject at different times, indicating that the type and severity of the disorder is determined to a large extent by the individual's response to the obstruction. The present report reviews the various mechanisms involved in the response to sleep-induced obstructive events. When the obstructive event takes the form of mild-moderate flow limitation, compensation can take place through an increase in the fraction of time spent in inspiration (Ti/Ttot) without any increase in maximum flow (V(MAX)). With more severe obstructions, V(MAX) must increase. Recent data indicate that the obstructed upper airway can reopen reflexly, without arousal, if chemical drive is allowed to reach a threshold (T(ER)) but that this is often preempted by a low arousal threshold. The relation between T(ER) and arousal threshold, as well as the lung-to-carotid circulation time and the rate of rise of chemical drive during the obstructive event, determine the magnitude of ventilatory overshoot at the end of an event and, by extension, whether initial obstructive events will be followed by stable breathing, slow evolving hypopneas with occasional arousals, or repetitive events.

Analysis of the interplay between neurochemical control of respiration and upper airway mechanics producing upper airway obstruction during sleep in humans

Increased loop gain (a function of both controller gain and plant gain), which results in instability in feedback control, is of major importance in producing recurrent central apnoeas during sleep but its role in causing obstructive apnoeas is not clear. The purpose of this study was to investigate the role of loop gain in producing obstructive sleep apnoeas. Owing to the complexity of factors that may operate to produce obstruction during sleep, we used a mathematical model to sort them out. The model used was based on our previous model of neurochemical control of breathing, which included the effects of chemical stimuli and changes in alertness on respiratory pattern generator activity. To this we added a model of the upper airways that contained a narrowed section which behaved as a compressible elastic tube and was tethered during inspiration by the contraction of the upper airway dilator muscles. These muscles in the model, as in life, responded to changes in hypoxia, hypercapnia and alertness in a manner similar to the action of the chest wall muscles, opposing the compressive action caused by the negative intraluminal pressure generated during inspiration which was magnified by the Bernoulli Effect. As the velocity of inspiratory airflow increased, with sufficiently large increase in airflow velocity, obstruction occurred. Changes in breathing after sleep onset were simulated. The simulations showed that increases in controller gain caused the more rapid onset of obstructive apnoeas. Apnoea episodes were terminated by arousal. With a constant controller gain, as stiffness decreased, obstructed breaths appeared and periods of obstruction recurred longer after sleep onset before disappearing. Decreased controller gain produced, for example, by breathing oxygen eliminated the obstructive apnoeas resulting from moderate reductions in constricted segment stiffness. This became less effective as stiffness was reduced more. Contraction of the upper airway muscles with hypercapnia and hypoxia could prevent obstructed apnoeas with moderate but not with severe reductions in stiffness. Increases in controller gain, as might occur with hypoxia, converted obstructive to central apnoeas. Breathing CO2 eliminated apnoeas when the activity of the upper airway muscles was considered to change as a function of CO2 to some exponent. Low arousal thresholds and increased upper airway resistance are two factors that promoted the occurrence and persistence of obstructive sleep apnoeas.

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.

Toward understanding respiratory sinus arrhythmia: relations to cardiac vagal tone, evolution and biobehavioral functions

Respiratory sinus arrhythmia (RSA, or high-frequency heart-rate variability) is frequently employed as an index of cardiac vagal tone or even believed to be a direct measure of vagal tone. However, there are many significant caveats regarding vagal tone interpretation: 1. Respiratory parameters can confound relations between RSA and cardiac vagal tone.2. Although intraindividual relations between RSA and cardiac vagal control are often strong, interindividual associations may be modest.3. RSA measurement is profoundly influenced by concurrent levels of momentary physical activity, which can bias estimation of individual differences in vagal tone.4. RSA magnitude is affected by beta-adrenergic tone.5. RSA and cardiac vagal tone can dissociate under certain circumstances.6. The polyvagal theory contains evolution-based speculations that relate RSA, vagal tone and behavioral phenomena. We present evidence that the polyvagal theory does not accurately depict evolution of vagal control of heart-rate variability, and that it ignores the phenomenon of cardiac aliasing and disregards the evolution of a functional role for vagal control of the heart, from cardiorespiratory synchrony in fish to RSA in mammals. Unawareness of these issues can lead to misinterpretation of cardiovascular autonomic mechanisms. On the other hand, RSA has been shown to often provide a reasonable reflection of cardiac vagal tone when the above-mentioned complexities are considered. Finally, a recent hypothesis is expanded upon, in which RSA plays a primary role in regulation of energy exchange by means of synchronizing respiratory and cardiovascular processes during metabolic and behavioral change.

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.

Causes of Cheyne-Stokes respiration

Cheyne-Stokes respiration (CSR) is one of several types of unusual breathing with recurrent apneas (dysrhythmias). Reported initially in patients with heart failure or stroke, it was then recognized both in other diseases and as a component of the sleep apnea syndrome. CSR is potentiated and perpetuated by changing states of arousal that occur during sleep. The recurrent hypoxia and surges of sympathetic activity that often occur during the apneas may have serious health consequences. Heart failure and stroke are risk factors for sleep apnea. The recurrent apneas and intermittent hypoxia occurring with sleep apnea further damage the heart and brain. Although all breathing dysrhythmias do not have the same cause, instability in the feedback control involved in the chemical regulation of breathing is the leading cause of CSR. Mathematical models have helped greatly in the understanding of the causes of recurrent apneas.

Mathematical models of periodic breathing and their usefulness in understanding cardiovascular and respiratory disorders

Periodic breathing is an unusual form of breathing with oscillations in minute ventilations and with repetitive apnoeas or near apnoeas. Reported initially in patients with heart failure or stroke, it was later recognized to occur especially during sleep. The recurrent hypoxia and surges of sympathetic activity that often occur during the apnoeas have serious health consequences. Mathematical models have helped greatly in the understanding of the causes of recurrent apnoeas. It is unlikely that every instance of periodic breathing has the same cause, but many result from instability in the feedback control involved in the chemical regulation of breathing caused by increased controller and plant gains and delays in information transfer. Even when it is not the main cause of the periodic breathing, unstable control modifies the ventilatory pattern and sometimes intensifies the recurrent apnoeas. The characteristics of disturbances to breathing and their interaction with the control system can be critical in determining ventilation responses and the occurrence of periodic breathing. Large abrupt changes in ventilation produced, for example, in the transition from waking to sleep and vice versa, or in the transition from breathing to apnoea, are potent factors causing periodic breathing. Mathematical models show that periodic breathing is a 'systems disorder' produced by the interplay of multiple factors. Multiple factors contribute to the occurrence of periodic breathing in congestive heart failure and cerebrovascular disease, increasing treatment options.

Role of acid-base balance in the chemoreflex control of breathing

This paper uses a steady-state modeling approach to describe the effects of changes in acid-base balance on the chemoreflex control of breathing. First, a mathematical model is presented, which describes the control of breathing by the respiratory chemoreflexes; equations express the dependence of pulmonary ventilation on Pco(2) and Po(2) at the central and peripheral chemoreceptors. These equations, with Pco(2) values as inputs to the chemoreceptors, are transformed to equations with hydrogen ion concentrations [H(+)] in brain interstitial fluid and arterial blood as inputs, using the Stewart approach to acid-base balance. Examples illustrate the use of the model to explain the regulation of breathing during acid-base disturbances. They include diet-induced changes in sodium and chloride, altitude acclimatization, and respiratory disturbances of acid-base balance due to chronic hyperventilation and carbon dioxide retention. The examples demonstrate that the relationship between Pco(2) and [H(+)] should not be neglected when modeling the chemoreflex control of breathing. Because pulmonary ventilation controls Pco(2) rather than the actual stimulus to the chemoreceptors, [H(+)], changes in their relationship will alter the ventilatory recruitment threshold Pco(2), and thereby the steady-state resting ventilation and Pco(2).

Physiological markers for anxiety: Panic disorder and phobias

Physiological activation is a cardinal symptom of anxiety, although physiological measurement is still not used for psychiatric diagnosis. An ambulatory study of phobics who were afraid of highway driving showed a concordance between self-reported anxiety during driving, autonomic activation, hypocapnia, and sighing respiration. Patients with panic attacks do not exhibit autonomic activation when they are quietly sitting and not having panic attacks, but do have the same respiratory abnormalities as driving phobics, suggesting that these abnormalities could be a marker for panic disorder. Such abnormalities are compatible with both the false suffocation alarm (D. Klein) and hyperventilation (R. Ley) theories of panic. Hypocapnia, however, is often absent during full-blown panic attacks. Since activation functions as preparation for physical activity, it may not occur when a patient has learned that avoidance of fear by flight or fight is futile. We developed a capnometry feedback assisted breathing training therapy for panic disorder designed to reduce hyperventilation and making breathing regular. Without feedback, conventional therapeutic breathing instructions may actually increase hyperventilation by increasing dyspnea. Five weekly therapy sessions accompanied by daily home practice with a capnometer produced marked clinical improvement compared to changes in an untreated group. Improvement was sustained over a 12-month follow-up period. The therapist avoided any statements or procedures designed to alter cognitions. Improvement occurred regardless of whether patients initially reported mostly respiratory or non-respiratory symptoms during their attacks. There is evidence that modifying any of the three systems comprising a fear network can be therapeutic, as exemplified by cognitive therapy modifying thoughts, exposure therapy modifying avoidance, and breathing training procedures modifying pCO(2).

Introduction of respiratory pattern generators into models of respiratory control

We have adapted two models previously proposed as respiratory pattern generators (RPGs) into a neurochemical feed back control model of ventilation. The RPG models, non-dimensional as originally presented, consisted of oscillating circuits of either two or five interconnected neurons [Matsugu, M., Duffin, J., Poon, C.-S., 1998. Entrainment, instability, quasi-periodicity, and chaos in a compound neural oscillator. J. Comput. Neurosci. 5, 35-51; Botros, S.M., Bruce, E.N., 1990. Neural network implementation of a three-phase model of respiratory rhythm generation. Biol. Cybern. 63, 143-153]. The neurochemical model into which they were integrated [Longobardo, G., Evangelisti, C.J., Cherniack, N.S., 2002. Effects of neural drives on breathing in the awake state in humans. Respir. Physiol. 129, 317-333] included the effects of cerebral blood flow variation with CO2, vagal stretch receptors input and a multicompartment model of carbon dioxide stores. The methodology is described whereby these neuronal oscillator networks were quantified, a necessary step for their inclusion as RPGs in broader models of the overall control of respiration. Subsequent simulations of the ventilation response to carbon dioxide with either respiratory pattern generator model exhibited only a limited range in which tidal volume and frequency increased with increasing respiratory drive. With both models, frequency peaked and then declined, as did ventilation when P CO2 was greater than normal. The range of the models was extended if the respiratory pattern generators were considered to be composed of multiple neuronal oscillators or a single oscillator in which there was increasing phasic input that was gated or pacemaker driven.

Physiological control of diving behaviour in the Weddell sealLeptonychotes weddelli: a model based on cardiorespiratory control theory

Despite being obligate air breathers, many species of marine mammal are capable of spending most of their lives submerged in water. How they do this has been a subject of intense interest to physiologists for over a century,yet we still do not have a detailed understanding of the physiological mechanisms underlying this behaviour. What are the proximate mechanisms that trigger the 'decisions' to submerge and return to the surface? The present study proposes a model intended to address this question, based on fundamental concepts of cardiorespiratory control. Two basic hypotheses are examined by computer simulation, using a mathematical model of the mammalian cardiorespiratory control system with parameter values for an adult Weddell seal: (1) that the control of diving can be considered to be a respiratory control problem, and (2) that dives are initiated and maintained by disfacilitation of respiratory drive, not inhibition. Computer simulations confirmed the plausibility of these hypotheses. Simulated diving behaviour and physiological responses (ventilation, cardiac output, blood and tissue gas tensions) were consistent with published data from freely diving Weddell seals. Dives up to the estimated aerobic dive limit (ADL, 18-25 min) could be simulated without the need for active inhibition of breathing in this model. This theoretical analysis suggests that the most important physiological adjustments occur during the surface interval phase of the dive cycle and include hyperventilation accompanied by high cardiac output, appropriate regulation of cerebral blood flow and central chemoreceptor threshold shifts. During dives, cardiac output, distribution of peripheral blood flow, splenic contraction and peripheral chemoreflex drives were found to modulate physiological and behavioural responses, but were not essential for simulated dives to occur. The main conclusion from this study is that the central chemoreceptor may be an important mechanism involved in the regulation of diving behaviour, implying that CO2, not O2, is the key regulatory variable in this model. This model includes and extends the ADL concept and suggests an explicit mechanism by which the respiratory control system may play a central role in the regulation of diving behaviour. It is likely that respiratory mechanisms are an important component of a hierarchical behavioural control system and further studies are required to test the qualitative and quantitative validity of the model.

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.

Orthostatic modification of ventilatory dynamic response to carbon dioxide perturbations

In order to determine whether changes in ventilatory control contribute to the observed decrease in arterial partial pressure of carbon dioxide (PaCO(2)) during head up tilt, we assessed ventilatory dynamic sensitivity to changes in PaCO(2) during supine and 70 degrees passive head up tilt. In 24 adult normals, we stimulated the ventilatory control system by switching inspired CO(2) between room air and room air+5% CO(2) in a pseudo random binary sequence. A Box-Jenkins model was used to compute ventilatory response to CO(2). Airflow, CO(2), non-invasive beat by beat blood pressure, ECG and cerebral blood flow velocity (Doppler) were recorded. During tilt, sensitivity of the ventilatory controller to CO(2) disturbance increased (from 0.45 to 0.72 L/min/mm Hg, p<0.005); minute ventilation increased (7.63 to 8.47 L/min, p<0.01), end tidal CO(2) (ETCO(2)), cerebral blood flow velocity (CBF) and baroreflex sensitivity decreased (46.9 to 42.9 mm Hg, p<0.001; 84.9 to 72.9 cm/s, p<0.001; and 17.6 to 5.5 ms/mm Hg, p<0.001). The primary observation from our study was that the sensitivity of ventilatory control system to perturbations in ETCO(2) increased during tilt. Taken together with decrease in mean levels of ETCO(2) and an increase in minute ventilation, these results suggest that during tilt, a change in the regulated level or 'set point' of PaCO(2) may occur.

The ventilatory responsiveness to CO(2) below eupnoea as a determinant of ventilatory stability in sleep

Sleep unmasks a highly sensitive hypocapnia-induced apnoeic threshold, whereby apnoea is initiated by small transient reductions in arterial CO(2) pressure (P(aCO(2))) below eupnoea and respiratory rhythm is not restored until P(aCO(2)) has risen significantly above eupnoeic levels. We propose that the 'CO(2) reserve' (i.e. the difference in P(aCO(2)) between eupnoea and the apnoeic threshold (AT)), when combined with 'plant gain' (or the ventilatory increase required for a given reduction in P(aCO(2))) and 'controller gain' (ventilatory responsiveness to CO(2) above eupnoea) are the key determinants of breathing instability in sleep. The CO(2) reserve varies inversely with both plant gain and the slope of the ventilatory response to reduced CO(2) below eupnoea; it is highly labile in non-random eye movement (NREM) sleep. With many types of increases or decreases in background ventilatory drive and P(aCO(2)), the slope of the ventilatory response to reduced P(aCO(2)) below eupnoea remains unchanged from control. Thus, the CO(2) reserve varies inversely with plant gain, i.e. it is widened with hyperventilation and narrowed with hypoventilation, regardless of the stimulus and whether it acts primarily at the peripheral or central chemoreceptors. However, there are notable exceptions, such as hypoxia, heart failure, or increased pulmonary vascular pressures, which all increase the slope of the CO(2) response below eupnoea and narrow the CO(2) reserve despite an accompanying hyperventilation and reduced plant gain. Finally, we review growing evidence that chemoreceptor-induced instability in respiratory motor output during sleep contributes significantly to the major clinical problem of cyclical obstructive sleep apnoea.

A theoretical study of the effect of circadian rhythms on sleep-induced periodic breathing and apnoea

This study employed a mathematical model of the respiratory control system to test the plausibility of the hypothesis that circadian rhythms in respiratory control can significantly influence respiratory stability at sleep onset. Computer simulations utilized a standardized "normal" sleep onset effect, superimposed upon systematic changes in chemoreflex parameters that mimicked the peaks and troughs of normal and high amplitude circadian rhythms. The analysis predicted that circadian influences may augment sleep-induced periodic breathing in nocturnal sleep compared with daytime naps. Furthermore, increased circadian amplitude of chemoreflex threshold, or absence of a circadian rhythm in peripheral chemosensitivity, each acted to stabilize respiration during daytime sleep onset and promote periodic breathing during nocturnal sleep onset. High amplitude circadian rhythms in respiratory control were predicted to cause an increasing number and duration of obstructive apnoeas from early to late night. It is suggested that the circadian timing system creates a nocturnal window of respiratory vulnerability and that abnormal circadian rhythms could potentially induce nocturnal sleep apnoea, even in individuals with normal sleep mechanisms.

Oxygen sensing: applications in humans

Our concepts of oxygen sensing have been transformed over the years. We now appreciate that oxygen sensing is not a unique property limited to “chemoreceptors” but is a common property of tissues and that responses to changes in oxygen levels are not static but can change over time. Respiratory responses initiated at the carotid body are modified by the excitatory and depressant effects of hypoxia at the brain and on the pathways connecting the carotid body to the brain. Equally important is that we are beginning to use our understanding of the cellular and molecular pathways triggered by hypoxia and hyperoxia to identify therapeutic targets to treat diseases such as cancer. We also have a better understanding of the complexities of the human respiratory responses to hypoxia; however, major deficiencies remain in our ability to alter or even measure human ventilatory responses to oxygen deficiency.

Chemoreflex control of ventilation is altered during wakefulness in humans with OSA

We hypothesized that patients with obstructive sleep apnea (OSA) have a different awake ventilatory response to carbon dioxide above and below eupnea compared with normal. Eight male subjects with OSA and control subjects matched for gender, race, age, height and weight voluntarily hyperventilated during wakefulness to reduce the partial pressure of carbon dioxide (PET(CO2)) below 25 mmHg. Subjects were then switched into a rebreathing bag containing a normocapnic (42 mmHg) hypoxic [partial pressure of end tidal oxygen (PET(O2))=50 mmHg (H50)] or hyperoxic [PET(O2)=140 mmHg (H140)] gas mixture. During the trial PET(CO2) increased while PET(O2) was maintained at a constant level. The point at which ventilation and PET(CO2) increased linearly was considered to be the carbon dioxide ventilatory recruitment threshold (VRT(CO2)). Measurements of ventilation and its components (i.e. tidal volume and breathing frequency) were made below this threshold and the slope of the minute ventilation; tidal volume or breathing frequency response above the threshold was determined. Four trials for a given oxygen level were completed. The PET(CO2) that demarcated the VRT(CO2) was increased (H(50)=43.43+/-0.92 vs. 41.05+/-0.67; H(140)=47.65+/-0.80 vs. 45.28+/-0.75), as were measures of ventilation below the threshold (H(50)=18.50+/-2.11 vs. 13.44+/-1.43; H(140)=19.66+/-2.71 vs. 10.83+/-1.24) in the OSA subjects compared with control. In contrast the OSA and control subjects did not respond differently to changes in PET(CO2) above the threshold. We conclude that the PET(CO2) that delineates the VRT(CO2) and ventilation below this threshold is elevated in subjects with OSA.

Piritramide and alfentanil display similar respiratory depressant potency

BACKGROUND

The question whether some opioids exert less respiratory depression than others has not been answered conclusively. We applied pharmacokinetic/pharmacodynamic (PKPD) modeling to obtain an estimate of the C50 for the depression of CO2 elimination as a measure of the respiratory depressant potency of alfentanil and piritramide, two opioids with vastly different pharmacokinetics and apparent respiratory depressant action.

METHODS

Twenty-three patients received either alfentanil (2.3 microg x kg(-1) x min-1, 14 patients, as published previously) or piritramide (17.9 microg x kg(-1) x min(-1), nine patients) until significant respiratory depression occurred. Opioid pharmacokinetics and the arterial PCO2 (PaCO2) were determined from frequent arterial blood samples. An indirect response model accounting for the respiratory stimulation due to increasing PaCO2 was used to describe the PaCO2 data.

RESULTS

The following pharmacodynamic parameters were estimated with NONMEM [population means and interindividual variability (CV)]: k(elCO2) (elimination rate constant of CO2) 0.144 (-) min(-1), F (gain of the CO2 response) 4.0 (fixed according to literature values) (28%), C50 (both drugs) 61.3 microg l-1 (41%), k(eo alfentanil) 0.654 (-) min(-1) and k(eo piritramide) 0.023 (-) min(-1). Assigning separate C50 values for alfentanil and piritramide did not improve the fit compared with a model with the same C50.

CONCLUSION

Since the C50 values did not differ, both drugs are equally potent respiratory depressants. The apparently lower respiratory depressant effect of piritramide when compared with alfentanil is caused by slower equilibration between the plasma and the effect site. Generalizing our results and based on simulations we conclude that slowly equilibrating opioids like piritramide are intrinsically safer with regard to respiratory depression than rapidly equilibrating opioids like alfentanil.

A mathematical model of CO2 variation in the ventilated neonate

A mathematical model of the variation of partial pressure of carbon dioxide in the arterial blood of a ventilated neonate is developed. The model comprises alveolar, arterial, pulmonary, venous and tissue compartments, with gas exchange in the lung determined by inspiration and expiration terms. Gas exchange is modelled through diffusion and convective transfer. Carbon dioxide is produced in the tissue by a metabolic term. Shunting is modelled by allowing blood flow to bypass the pulmonary compartment in which diffusion takes place. The model predicts changes in the carbon dioxide partial pressures that occur following abrupt changes in the ventilation settings, and show broad agreement with actual data obtained from novel sensing technology.