Bench and mathematical modeling of the effects of breathing a helium/oxygen mixture on expiratory time constants in the presence of heterogeneous airway obstructions

ABSTRACT: BACKGROUND: Expiratory time constants are used to quantify emptying of the lung as a whole, and emptying of individual lung compartments. Breathing low-density helium/oxygen mixtures may modify regional time constants so as to redistribute ventilation, potentially reducing gas trapping and hyperinflation for patients with obstructive lung disease. In the present work, bench and mathematical models of the lung were used to study the influence of heterogeneous patterns of obstruction on compartmental and whole-lung time constants. METHODS: A two-compartment mechanical test lung was used with the resistance in one compartment held constant, and a series of increasing resistances placed in the opposite compartment. Measurements were made over a range of lung compliances during ventilation with air or with a 8/22% mixture of helium/oxygen. The resistance imposed by the breathing circuit was assessed for both gases. Experimental results were compared with predictions of a mathematical model applied to the test lung and breathing circuit. In addition, compartmental and whole-lung time constants were compared with those reported by the ventilator. RESULTS: Time constants were greater for larger minute ventilation, and were reduced by substituting helium/oxygen in place of air. Notably, where time constants were long due to high lung compliance (i.e. low elasticity), helium/oxygen improved expiratory flow even for a low level of resistance representative of healthy, adult airways. In such circumstances, the resistance imposed by the external breathing circuit was significant. Mathematical predictions were in agreement with experimental results. Time constants reported by the ventilator were wellcorrelated with those determined for the whole-lung and for the low-resistance compartment, but poorly correlated with time constants determined for the high-resistance compartment. CONCLUSIONS: It was concluded that breathing a low-density gas mixture, such as helium/oxygen, can improve expiratory flow from an obstructed lung compartment, but that such improvements will not necessarily affect time constants measured by the ventilator. Further research is required to determine if alternative measurements made at the ventilator level are predictive of regional changes in ventilation. It is anticipated that such efforts will be aided by continued development of mathematical models to include pertinent physiological and pathophysiological phenomena that are difficult to reproduce in mechanical test systems.     

Effect of inspiratory muscle fatigue on breathing pattern

Our aim was to determine whether inspiratory muscle fatigue changes breathing pattern and whether any changes seen occur before mechanical fatigue develops. Nine normal subjects breathed through a variable inspiratory resistance with a predetermined mouth pressure (Pm) during inspiration and a fixed ratio of inspiratory time to total breath duration. Breathing pattern after resistive breathing (recovery breathing pattern) was compared with breathing pattern at rest and during CO2 rebreathing (control breathing pattern) for each subject. Relative rapid shallow breathing was seen after mechanical fatigue and also in experiments with electromyogram evidence of diaphragmatic fatigue where Pm was maintained at the predetermined level during the period of resistive breathing. In contrast there was no significant difference between recovery and control breathing patterns when neither mechanical nor electromyogram fatigue was seen. It is suggested that breathing pattern after inspiratory muscle fatigue changes in order to minimize respiratory sensation.     

Ontogeny of respiratory control and pulmonary mechanics in newborn rabbits

Maturation of the respiratory pattern and the active and passive mechanical properties of the respiratory system were assessed in 19 tracheotomized rabbits (postnatal age range: 1-26 days) placed in a body plethysmograph. With maturation both minute ventilation and tidal volume significantly increased, whereas respiratory frequency decreased. When normalized for body weight (kg) both the passive (Rrs X kg) and active (R'rs X kg) resistances of the respiratory system significantly increased with age, whereas the corresponding passive (Crs X kg-1) and active (C'rs X kg-1) compliances significantly decreased. At any given age R'rs X kg only slightly exceeded Rrs X kg, whereas C'rs X kg-1 was significantly lower than Crs X kg-1. Moreover, the maturational increases in Rrs X kg and R'rs X kg exceeded the corresponding decreases in Crs X kg-1 and C'rs X kg-1, resulting in significant age-related increases in both the passive (tau rs) and active (tau'rs) time constants of the respiratory system. Due to the age-related increases in tau'rs, producing a delayed volume response to any given inspiratory driving pressure, the relative volume loss obtained at any time during inspiration was greater in the maturing rabbit. On the other hand, because of concomitant compensatory changes in respiratory pattern, evidenced by increases in inspiratory duration with age, the end-inspiratory tidal volume loss in the maturing animal was maintained generally less than 10% at all postnatal ages. Thus maturational changes in respiratory pattern appear coupled to changes in the active mechanical properties of the respiratory system. The latter coupling serves to optimize the transduction of inspiratory pressure into volume change in a manner consistent with establishing the minimum inspiratory work of breathing during postnatal development.     

Mechanical impedance as determinant of inspiratory neural drive during exercise in humans

Five healthy males exercised progressively with small 2-min increments in work load. We measured inspiratory drive (occlusion pressure, P0.1), pulmonary resistance (RL), dynamic pulmonary compliance (Cdyn), transdiaphragmatic pressure (Pdi), and diaphragmatic electromyogram (EMGdi). Minute ventilation (VE), mean inspiratory flow rate (VT/TI), and P0.1 all increased exponentially with increased work load, but P0.1 increased at a faster rate than did VT/TI or VE. Thus effective impedance (P0.1/VT/TI) rose throughout exercise. The increasing P0.1 was mostly due to augmented Pdi and coincided with increased EMGdi during this initial portion of inspiration. We found no consistent change in RL or Cdyn throughout exercise. With He breathing (80% He-20% O2), RL was reduced at all work loads; P0.1 fell in comparison with air-breathing values and VE, VT, and VT/TI rose in moderate and heavy work; and P0.1/VT/TI was unchanged with increasing exercise loads. Step reductions in gas density at a constant work load of any intensity showed an immediate reduction in the rate of rise of EMGdi and Pdi followed by increased VT/TI, breathing frequency, and hypocapnia. These changes were maintained during prolonged periods of unloading and were immediately reversible on return to air breathing. These data are consistent with the existence of a reflex effect on the magnitude of inspiratory neural drive during exercise that is sensitive to the load presented by the normal mechanical time constant of the respiratory system. This "load" is a significant determinant of the hyperpneic response and thus of the maintenance of normocapnia during exercise.     

Effects of vagotomy on respiratory mechanics in newborn and adult pigs

We have examined breathing patterns and respiratory mechanics in anesthetized tracheostomized newborn piglets and adult pigs and the changes determined by cervical bilateral vagotomy. Piglets had a respiratory system compliance and resistance, on a per kilogram basis, respectively, higher and smaller than the adults. After vagotomy neither variable changed in the newborn, but resistance dropped in the adult. This may suggest that efferent vagal control of bronchomotor tone is more pronounced in the adult. Respiratory system time constant was longer in newborns both before and after vagotomy. The distortion of the chest wall, examined as the ratio between the volume inhaled spontaneously and the passive volume for the same abdominal motion, was more marked in newborns, reflecting their higher chest wall compliance. The work per minute, computed from the pressure and volume changes, was larger in piglets. After vagotomy the external work per minute was not different; however, the larger tidal volumes were accompanied by a larger chest distortion. This may indicate that vagal control of the breathing pattern, by limiting the depth of inspiration and hence the amount of chest distortion, has implications on the energetics of breathing.     

Pulmonary mechanics during rapid mechanical ventilation in rabbits with saline-lavaged lungs

Infants with respiratory failure are frequently mechanically ventilated at rates exceeding 60 breaths/min. We analyzed the effect of ventilatory rates of 30, 60, and 90 breaths/min (inspiratory times of 0.6, 0.3, and 0.2 s, respectively) on the pressure-flow relationships of the lungs of anesthetized paralyzed rabbits after saline lavage. Tidal volume and functional residual capacity were maintained constant. We computed effective inspiratory and expiratory resistance and compliance of the lungs by dividing changes in transpulmonary pressure into resistive and elastic components with a multiple linear regression. We found that mean pulmonary resistance was lower at higher ventilatory rates, while pulmonary compliance was independent of ventilatory rate. The transpulmonary pressure developed by the ventilator during inspiration approximated a linear ramp. Gas flow became constant and the pressure-volume relationship linear during the last portion of inspiration. Even at a ventilatory rate of 90 breaths/min, 28-56% of the tidal volume was delivered with a constant inspiratory flow. Our findings are consistent with the model of Bates et al. (J. Appl. Physiol. 58: 1840-1848, 1985), wherein the distribution of gas flow within the lungs depends predominantly on resistive factors while inspiratory flow is increasing, and on elastic factors while inspiratory flow is constant. This dynamic behavior of the surfactant-depleted lungs suggests that, even with very short inspiratory times, distribution of gas flow within the lungs is in large part determined by elastic factors. Unless the inspiratory time is further shortened, gas flow may be directed to areas of increased resistance, resulting in hyperinflation and barotrauma.     

Breathing cardiovascular variability and baroreflex in mechanically ventilated patients

Heart rate and blood pressure variations during spontaneous ventilation are related to the negative airway pressure during inspiration. Inspiratory airway pressure is positive during mechanical ventilation, suggesting that reversal of the normal baroreflex-mediated pattern of variability may occur. We investigated heart rate and blood pressure variability and baroreflex sensitivity in 17 mechanically ventilated patients. ECG (RR intervals), invasive systolic blood pressure (SBP), and respiratory flow signals were recorded. High-frequency (HF) amplitude of RR and SBP time series and HF phase differences between RR, SBP, and ventilatory signals were continuously computed by Complex DeModulation (CDM). Cross-spectral analysis was used to assess the coherence and the gain functions between RR and SBP, yielding baroreflex sensitivity indices. The HF phase difference between SBP and ventilatory signals was nearly constant in all patients with inversion of SBP variability during the ventilator cycle compared with cycling with negative inspiratory pressure to replicate spontaneous breathing. In 12 patients (group 1), the phase difference between RR and ventilatory signals changed over time and the HF-RR amplitude varied. In the remaining five patients (group 2), RR-ventilatory signal phase and HF-RR amplitude showed little change; however, only one of these patients exhibited a RR-ventilatory signal phase difference mimicking the normal pattern of respiratory sinus arrhythmia. Spectral coherence between RR and SBP was lower in the group with phase difference changes. Positive pressure ventilation exerts mainly a mechanical effect on SBP, whereas its influence on HR variability seems more complex, suggesting a role for neural influences.     

Changes in mechanics of breathing during experimental cough and sneeze in anaesthetized cats

Pleural pressure, airflow and tidal volume during experimental cough and sneeze elicited by mechanical stimulation of the tracheobronchial and nasal mucous membranes were investigated in fifty anaesthetized cats (pentobarbital, 40 mg/kg i.p.). Pressure-volume, pressure-flow and flow-volume relations were studied during these expulsive processes. In comparison to quiet breathing there was a decrease in dynamic lung compliance in both respiratory tract reflexes (p less than 0.001), especially in their expiratory phases. As compared to quiet breathing, the total work of breathing was significantly increased (p less than 0.001) in cough (20 times) as well as in sneeze (13 times). The total lung resistance increased markedly (p less than 0.001) in both cough and sneeze compared to quiet breathing. In these expulsive processes there was also a high "cough index" (resistance calculated from the peak flow and instantaneous pressure). The flow-volume curve in cough, in contradistinction to sneeze, indicated a significantly reduced airflow of the end of expiration (at 85% of the expired volume), demonstrating a concomitant bronchoconstriction.     

Effects of posture on the respiratory mechanics of the chick embryo

In the chicken embryo, pulmonary ventilation and pulmonary gas exchange begin approximately one day before the completion of hatching. We asked to what extent the posture inside the egg, and the presence of the eggshell and membranes, may alter the mechanical behaviour of the respiratory system. The passive mechanical properties of the respiratory system were studied in chicken embryos during the internal pipping phase (rupture of the air cell) or the external pipping phase (hole in the eggshell). Tracheal pressure and changes in lung volume were recorded during mechanical ventilation, first, with the embryo curled up inside the egg, then again after exteriorization from the eggshell. In the internal pippers, respiratory system compliance increased and expiratory resistance decreased after exteriorization, whereas the mean inspiratory impedance did not change. In the external pippers, exteriorization had no significant effects on respiratory compliance, resistance, or impedance, and the values were similar to those of newly hatched chicks. We conclude that, in the chicken embryo, at a time when pulmonary ventilation becomes an important mechanism for gas exchange, the curled up posture inside the egg does not provide any significant mechanical constraint to breathing.     

Transmural pressure in the right atrium during positive pressure respiration in cats

The influence of continuous positive pressure breathing (cm H2O) on the breathing mechanics, central venous pressure, and transmural pressure in the right atrium, were studied in anaesthetised cats separately during inspiration and expiration. It's shown that hemodynamics effects are directly connected with the influence of increased intrathoracic pressure during whole breathing cycles in contrast with the phase changes in natural expiration and inspiration. The inversion of relation of intrathoracic and central venous pressure due to displacement of the mechanical respiratory characteristics became the factors defining the fall of the right atrium filling pressure.     

Alveolar pressure inhomogeneity during low-frequency oscillation of excised canine lobes

To quantify the inhomogeneity of alveolar pressures (PA) during cyclic changes in lung volume similar to those present during spontaneous breathing, inhomogeneity of PA was measured with an alveolar capsule technique in six excised canine lungs. The lungs were ventilated by a quasi-sinusoidal pump with a constant end-expiratory lung volume and tidal volumes of 10, 20, and 40% of vital capacity at breathing frequencies ranging from 5 to 45 breaths/min. Inhomogeneity of PA was quantified as the sample standard deviation of pressures measured in three capsules. A component of inhomogeneity in phase with flow and a smaller component out of phase with flow were present. The in-phase component increased approximately linearly with flow. The ratio of inhomogeneity to flow was smaller at large tidal volumes and, at the two higher tidal volumes studied, the ratio was greater during inspiration than during expiration. If these data are interpreted in terms of a simple circuit model, this degree of inhomogeneity implies an approximately twofold variation in regional time constants. Despite these considerable differences in time constants, the absolute amount of inhomogeneity as defined by the sample standard deviation of the three PA's was small (maximum 0.57 +/- 0.32 cmH2O at the highest breathing frequency and tidal volume) because airway resistance in the canine lung was small.     

Chest wall and respiratory system mechanics in cats: effects of flow and volume

The effects of inspiratory flow rate and inflation volume on the resistive properties of the chest wall were investigated in six anesthetized paralyzed cats by use of the technique of rapid airway occlusion during constant flow inflation. This allowed measurement of the intrinsic resistance (Rw,min) and overall dynamic inspiratory impedance (Rw,max), which includes the additional pressure losses due to time constant inequalities within the chest wall tissues and/or stress adaptation. These results, together with our previous data pertaining to the lung (Kochi et al., J. Appl. Physiol. 64: 441-450, 1988), allowed us to determine Rmin and Rmax of the total respiratory system (rs). We observed that 1) Rw,max and Rrs,max exhibited marked frequency dependence; 2) Rw,min was independent of flow (V) and inspired volume (delta V), whereas Rrs,min increased linearly with V and decreased with increasing delta V; 3) Rw,max decreased with increasing V, whereas Rrs,max exhibited a minimum value at a flow rate substantially higher than the resting range of V; 4) both Rw,max and Rrs,max increased with increasing delta V. We conclude that during resting breathing, flow resistance of the chest wall and total respiratory system, as conventionally measured, includes a significant component reflecting time constant inequalities and/or stress adaptation phenomena.     

Changes in breathing pattern at loads near perceptual threshold at different work levels

Five subjects were tested to determine the threshold for detection of an added resistance to inspiration in three tests, one at rest and two with exercise (mild = 50 W; moderate = 100 W) on a cycle ergometer. Changes in the breathing pattern were examined at added resistances near the perceptual threshold. Added inspiratory resistances with a 50% probability of detection were very variable at rest; they decreased significantly from rest (250 Pa.l-1.s-1) to moderate exercise (98 Pa.l-1.s-1) in four subjects. It is suggested that physical exercise may cause discomfort even when workers wearing a respirator do not have any abnormal sensation during sedentary work. Breathing patterns were compared between resistance loaded and unloaded breathing during each test. Decreases in inspiratory peak flow and acceleration of flow early in inspiration were found in resistance loaded breathing in almost all tests and a tendency for tidal volume to decrease was found during moderate exercise only. The ratios of resistance loaded to unloaded breathing for inspiratory time (ti) and total time (tt) tended to be greater in the detected than in the undetected responses at rest and during mild exercise but not during moderate exercise. This would imply that further prolongation of ti and tt in the detected responses was attributable to conscious or subconscious aspects of the resistance leading responses: however, these adjustments in breathing, which reduce frequency, would be less likely to occur as the work rate increases.     

Robust Unidirectional Airflow through Avian Lungs: New Insights from a Piecewise Linear Mathematical Model

Avian lungs are remarkably different from mammalian lungs in that air flows unidirectionally through rigid tubes in which gas exchange occurs. Experimental observations have been able to determine the pattern of gas flow in the respiratory system, but understanding how the flow pattern is generated and determining the factors contributing to the observed dynamics remains elusive. It has been hypothesized that the unidirectional flow is due to aerodynamic valving during inspiration and expiration, resulting from the anatomical structure and the fluid dynamics involved, however, theoretical studies to back up this hypothesis are lacking. We have constructed a novel mathematical model of the airflow in the avian respiratory system that can produce unidirectional flow which is robust to changes in model parameters, breathing frequency and breathing amplitude. The model consists of two piecewise linear ordinary differential equations with lumped parameters and discontinuous, flow-dependent resistances that mimic the experimental observations. Using dynamical systems techniques and numerical analysis, we show that unidirectional flow can be produced by either effective inspiratory or effective expiratory valving, but that both inspiratory and expiratory valving are required to produce the high efficiencies of flows observed in avian lungs. We further show that the efficacy of the inspiratory and expiratory valving depends on airsac compliances and airflow resistances that may not be located in the immediate area of the valving. Our model provides additional novel insights; for example, we show that physiologically realistic resistance values lead to efficiencies that are close to maximum, and that when the relative lumped compliances of the caudal and cranial airsacs vary, it affects the timing of the airflow across the gas exchange area. These and other insights obtained by our study significantly enhance our understanding of the operation of the avian respiratory system.     

Termination of inspiration by phase-dependent respiratory vagal feedback in awake normal humans.

Imperceptible levels of proportional assist ventilation applied throughout inspiration reduced inspiratory time (TI) in awake humans. More recently, the reduction in TI was associated with flow assist, but flow assist also reaches a maximum value early during inspiration. To test the separate effects of flow assist and timing of assist, we applied a pseudorandom binary sequence of flow-assisted breaths during early, late, or throughout inspiration in eight normal subjects. We hypothesized that imperceptible flow assist would shorten TI most effectively when applied during early inspiration. Tidal volume, integrated respiratory muscle pressure per breath, TI, and TE were recorded. All stimuli (early, late, or flow assist applied throughout inspiration) resulted in a significant increase in inspiratory flow; however, only when the flow assist was applied during early inspiration was there a significant reduction in TI and the integrated respiratory muscle pressure per breath. These results provide further evidence that vagal feedback modulates breathing on a breath-by-breath basis in conscious humans within a physiological range of breath sizes.     

Expiratory pattern of newborn mammals

The passive mechanical time constant (tau pass) of the respiratory system is relatively similar among newborn mammalian species, approximately 0.15-0.2 s. However, breathing rate (f) is higher in smaller species than larger species in order to accommodate the relatively larger metabolic demands. Since tidal volume per kilogram is an interspecies constant, in the fastest breathing species the short expiratory time should determine a substantial dynamic elevation of the functional residual capacity (FRC). We examined the possibility of a difference in expiratory time constant between dynamic and passive conditions by analyzing the expiratory flow pattern of nine newborn unanesthetized species during resting breathing. In most newborns the late portion of the expiratory flow-volume curve was linear, suggesting muscle relaxation. The slope of the curve, which represents the dynamic expiratory time constant of the respiratory system (tau exp), varied considerably among animals (from 0.1 to 0.7 s), being directly related to the inspiratory time and inversely proportional to f. In relatively slow-breathing newborns, such as infants and piglets, tau exp is longer than tau pass most likely due to an increase in the expiratory laryngeal resistance and FRC is substantially elevated. On the contrary, in the fastest breathing newborns (such as rats and mice) tau exp is similar or even less than tau pass, because at these high rates dynamic lung compliance is lower than its passive value and the dynamic elevation of FRC is small. In dynamic conditions, therefore, the product of tau exp and f is maintained within narrow limits.(ABSTRACT TRUNCATED AT 250 WORDS)     

Laryngeal muscle response to phasic and tonic upper airway pressure and flow.

The hypothesis that respiratory modulation due to upper airway (UA) pressure and flow is dependent on stimulus modality and respiratory phase-specific activation was assessed in anesthetized, tracheotomized, spontaneously breathing piglets. Negative pressure and flow applied to the isolated UA at room or body temperature during inspiration only enhanced posterior cricoarytenoid muscle activity from that present without UA pressure and flow (baseline) by 15--20%. Time shifting the onset of UA flow relative to tracheal flow decreased this enhancement. The same enhancement was observed with oscillatory or constant airflow. UA positive pressure and flow at room or body temperature applied during expiration only enhanced thyroarytenoid muscle activity from baseline by 50--160%. The same enhancement was observed with oscillatory or constant airflow at body temperature. Constant positive pressure and flow enhanced thyroarytenoid muscle activity more than oscillatory pressure and flow at room temperature. We conclude that the respiratory modulation of UA afferents is processed in a phase-specific fashion and is dependent on stimulus modality (tonic vs. phasic).     

Decay of inspiratory muscle pressure during expiration in conscious humans

In eight conscious spontaneously breathing adults we studied the decay of pressure developed by the inspiratory muscles during expiration (PmusI). PmusI was obtained according to the following equation: PmusI(t) = Ers X V(t) - Rrs X V(t), where V is volume and V is flow at any instant t during spontaneous expiration, and Ers and Rrs are, respectively, the passive elastance and resistance of the total respiratory system. Ers was determined with the relaxation method, and resistance with the interrupter method. All subjects showed marked braking of expiratory flow by PmusI. The mean time for PmusI to reduce to 50 and 0% amounted, respectively, to 23 and 79% of expiratory time. During expiration, 24-55% of the elastic energy stored during inspiration was used as resistive work and the remainder (45-76%) as negative work.     

Localization of the sites of pulmonary vasomotion by use of arterial and venous occlusion

In this study, we present a new approach for using the pressure vs. time data obtained after various vascular occlusion maneuvers in pump-perfused lungs to gain insight into the longitudinal distribution of vascular resistance with respect to vascular compliance. Occlusion data were obtained from isolated dog lung lobes under normal control conditions, during hypoxia, and during histamine or serotonin infusion. The data used in the analysis include the slope of the arterial pressure curve and the zero time intercept of the extrapolated venous pressure curve after venous occlusion, the equilibrium pressure after simultaneous occlusion of both the arterial inflow and venous outflow, and the area bounded by equilibrium pressure and the arterial pressure curve after arterial occlusion. We analyzed these data by use of a compartmental model in which the vascular bed is represented by three parallel compliances separated by two series resistances, and each of the three compliances and the two resistances can be identified. To interpret the model parameters, we view the large arteries and veins as mainly compliance vessels and the small arteries and veins as mainly resistance vessels. The capillary bed is viewed as having a high compliance, and any capillary resistance is included in the two series resistances. With this view in mind, the results are consistent with the major response to serotonin infusion being constriction of large and small arteries (a decrease in arterial compliance and an increase in arterial resistance), the major response to histamine infusion being constriction of small and large veins (an increase in venous resistance and a decrease in venous compliance), and the major response to hypoxia being constriction of the small arteries (an increase in arterial resistance). The results suggest that this approach may have utility for evaluation of the sites of action of pulmonary vasomotor stimuli.     

Endothelin involvement in respiratory centre activity.

To evaluate the role of endothelin (ET) in respiratory homeostasis we studied the effects of the ET(A) and ET(B) receptor blocking agent bosentan on respiratory mechanics and control in seven anaesthetised spontaneously breathing pigs, for 180 min after single bolus administration (20 mg/kg i.v.). The results show that the block of ET receptors induced a significant increase in compliance and decrease in resistance of the respiratory system, entailing a significant reduction of diaphragmatic electromyographic activity, without affecting the centroid frequency of the power spectrum. Bosentan administration induced a significant increase in tidal volume (V(T)), accompanied by a significant decrease in respiratory frequency, without any significant change in pulmonary ventilation, CO(2) arterial blood gas pressure or pH. Since the relationship between V(T) and inspiratory time remained substantially constant after bosentan administration, the changes in respiratory pattern appear to be the result of an upward shift in inspiratory off-switch threshold. Both inspiratory and expiratory times during occluded breathing were increased by block of ET receptors, suggesting also a central respiratory neuromodulator effect of ET. In conclusion the present results suggest that the block of ET receptors in spontaneously breathing pigs exerts a role on mechanical properties of the respiratory system as well as on peripheral and central mechanisms of breathing control.