Smooth muscle dynamics and maximal expiratory flow in asthma.

A computational model for maximal expiratory flow in constricted lungs is presented. The model was constructed by combining a previous computational model for maximal expiratory flow in normal lungs and a previous mathematical model for smooth muscle dynamics. Maximal expiratory flow-volume curves were computed for different levels of smooth muscle activation. The computed maximal expiratory flow-volume curves agree with data in the literature on flow in constricted nonasthmatic subjects. In the model, muscle force during expiration depends on the balance between the decrease in force that accompanies muscle shortening and the recovery of force that occurs during the time course of expiration, and the computed increase in residual volume (RV) depends on the magnitude of force recovery. The model was also used to calculate RV for a vital capacity maneuver with a slow rate of expiration, and RV was found to be further increased for this maneuver. We propose that the measurement of RV for a vital capacity maneuver with a slow rate of expiration would provide a more sensitive test of smooth muscle activation than the measurement of maximal expiratory flow.     

Slope-Based Stochastic Resonance: How Noise Enables Phasic Neurons to Encode Slow Signals

Fundamental properties of phasic firing neurons are usually characterized in a noise-free condition. In the absence of noise, phasic neurons exhibit Class 3 excitability, which is a lack of repetitive firing to steady current injections. For time-varying inputs, phasic neurons are band-pass filters or slope detectors, because they do not respond to inputs containing exclusively low frequencies or shallow slopes. However, we show that in noisy conditions, response properties of phasic neuron models are distinctly altered. Noise enables a phasic model to encode low-frequency inputs that are outside of the response range of the associated deterministic model. Interestingly, this seemingly stochastic-resonance (SR) like effect differs significantly from the classical SR behavior of spiking systems in both the signal-to-noise ratio and the temporal response pattern. Instead of being most sensitive to the peak of a subthreshold signal, as is typical in a classical SR system, phasic models are most sensitive to the signal''s rising and falling phases where the slopes are steep. This finding is consistent with the fact that there is not an absolute input threshold in terms of amplitude; rather, a response threshold is more properly defined as a stimulus slope/frequency. We call the encoding of low-frequency signals with noise by phasic models a slope-based SR, because noise can lower or diminish the slope threshold for ramp stimuli. We demonstrate here similar behaviors in three mechanistic models with Class 3 excitability in the presence of slow-varying noise and we suggest that the slope-based SR is a fundamental behavior associated with general phasic properties rather than with a particular biological mechanism.     

Adaptation and selective information transmission in the cricket auditory neuron AN2

Sensory systems adapt their neural code to changes in the sensory environment, often on multiple time scales. Here, we report a new form of adaptation in a first-order auditory interneuron (AN2) of crickets. We characterize the response of the AN2 neuron to amplitude-modulated sound stimuli and find that adaptation shifts the stimulus-response curves toward higher stimulus intensities, with a time constant of 1.5 s for adaptation and recovery. The spike responses were thus reduced for low-intensity sounds. We then address the question whether adaptation leads to an improvement of the signal's representation and compare the experimental results with the predictions of two competing hypotheses: infomax, which predicts that information conveyed about the entire signal range should be maximized, and selective coding, which predicts that "foreground" signals should be enhanced while "background" signals should be selectively suppressed. We test how adaptation changes the input-response curve when presenting signals with two or three peaks in their amplitude distributions, for which selective coding and infomax predict conflicting changes. By means of Bayesian data analysis, we quantify the shifts of the measured response curves and also find a slight reduction of their slopes. These decreases in slopes are smaller, and the absolute response thresholds are higher than those predicted by infomax. Most remarkably, and in contrast to the infomax principle, adaptation actually reduces the amount of encoded information when considering the whole range of input signals. The response curve changes are also not consistent with the selective coding hypothesis, because the amount of information conveyed about the loudest part of the signal does not increase as predicted but remains nearly constant. Less information is transmitted about signals with lower intensity.     

An integrated microprocessor-based system enhanced with graphics for the evaluation of slope-ratios and other flow-volume curve parameters

A microprocessor-based system is developed to measure and analyse flow-volume curves produced from forced expiratory manoeuvres. The system evaluates the slope-ratio parameter and displays it as a scaled plot against fractions of vital capacity together with a scaled plot of the flow-volume curves on a graphic display. Other conventional parameters and their predicted values are also calculated. Freedom is allowed for the operator to select or reject manoeuvres from several runs. The system is easily run and guided even by an operator with no programming knowledge.     

Self-similarity and the relationship between abundance and range size

Patterns in the relationships among the range, abundance, and distribution of species within a biome are of fundamental interest in ecology. A self-similarity condition, imposed at the community level and previously demonstrated to lead to the power-law form of the species-area relationship, is extended to the species level and shown to predict testable power-law relationships between range size and both species abundance and area of census cell across scales of spatial resolution. The predicted slopes of plots of log(range size) versus log(abundance) are shown to be in good agreement with data from British breeding bird and mammal censuses and with data on the distribution of fern species in old-growth forest. The predicted slopes of plots of log(range size) versus log (area of census cell) are consistent with the limited available data for British plant species. Self-similarity provides a testable theoretical framework for a unified understanding of patterns among the range, abundance, and distribution of species.     

Sensitivity and specificity of the computational model for maximal expiratory flow

The computational model for forced expiratory flow from human lungs of Lambert and associates (J. Appl. Physiol.: Respirat. Environ. Exercise Physiol. 52: 44-56, 1982) was used to investigate the sensitivity of maximal expiratory flow to lung properties. It was found that maximal flow is very sensitive to recoil pressure and airway areas but not very sensitive to lung volume, airway compliance, and airway length. Linear programming was used to show that a given air flow-pressure curves was compatible with a fairly wide range of airway properties. Additional data for maximal flow with a He-O2 mixture narrowed the range somewhat. It was shown that the flow-pressure curve contains more information about central than peripheral airways and that information about the latter is obtainable only from flows at recoils less than 2 cmH2O. Parameter ranges compatible with individual flow-pressure curves showed differences that demonstrated that such curves give some indication of individual central airway properties.     

Lung inflation does not increase maximal expiratory flow during induced obstruction in the dog

A deep inflation (DI) reverses induced bronchoconstriction in normal human subjects whether assessed by airway resistance before and after a DI or by isovolumic maximal expiratory flows (Vmax) from partial expiratory flow-volume (PEFV) vs. maximum expiratory flow-volume (MEFV) maneuvers. These observations suggest that with induced constriction the hysteresis of airways exceeds that of the parenchyma. In contrast with humans, a previous study of ours on dogs indicated that induced increases in airway resistance were unaffected by DI, suggesting that hysteresis of airways and parenchyma were equal. We hypothesized therefore that in constricted dog lungs, any differences that might arise in isovolumic Vmax between PEFV and MEFV maneuvers would not be due to changes in airway caliber but rather would be wholly determined by isovolumic differences in deflational recoil pressures. Recoil pressures were dynamically measured using six separate alveolar capsules in each of six dogs. At base line there were no significant differences between isovolumic recoil pressures or maximal flows with volume history, suggesting equal degrees of airway and parenchymal hysteresis. After histamine-induced constriction there were also no isovolumic differences in flows, but due to striking nonhomogeneities in dynamic recoil pressure among alveolar capsules, it was not possible to express a single meaningful recoil pressure pertinent to the lungs as a whole. These findings are consistent with the idea that isovolumic comparisons of Vmax serve as a reasonable indicator of changes in the relative degree of airway and parenchymal hysteresis.     

Heterogeneous lung emptying during forced expiration

Several lines of evidence suggest that the healthy mammalian lung empties homogeneously during a maximally forced deflation. Nonetheless, such behavior would appear to be implausible if for no other reason than that airway structure is known to be substantially heterogeneous among parallel pathways of gas conduction. To resolve this paradox we reexamined the degree to which lung emptying is homogeneous, and considered mechanisms that might control differential regional emptying. Twelve excised canine lungs were studied. Regional alveolar pressure relative to pleural pressure was used as an index of regional lung volume. By use of a capsule technique, alveolar pressure was measured simultaneously in each of six regions during flow-limited deflations; flow from the lung was measured plethysmographically. The standard deviation of interregional pressure differences, which was taken as an index of nonuniformity, was 0.0, 0.74, 0.64, and 0.90 cmH2O at mean recoil pressures of 30, 8.4, 4.5, and 2.1 cmH2O (0, 25, 50, and 75% expired vital capacity), indicating that interregional pressure differences increased more rapidly earlier in the deflation. When we examined the time rate of change of regional alveolar pressure as an index of regional flow, we observed an intricate pattern of differential regional behavior that was inapparent in the maximum expiratory flow-volume (MEFV) curve. The most plausible interpretation of these findings is that regions of the healthy excised canine lung empty heterogeneously to a small degree, but in an interdependent compensatory pattern that is inapparent in the configuration of the maximum expiratory flow-volume curve.     

Ventilatory functions of normal rats of different ages

1. Lung volume, respiratory mechanics and flow-volume curve were measured in 96 male Wistar rats weighing from 102 to 521 g. Their ages ranged from 31 to 137 days and therefore, most of them were in the rapidly growing phase. 2. Data except for functional residual capacity were obtained under conditions of artificial ventilation. 3. Total lung capacity and the fractions of lung volume of these animals were found to relate linearly with their body weights. 4. Most of the parameters concerning respiratory mechanics and flow-volume curve also changed generally in proportion to their body weights, but those standardized for lung volume were not constant but varied depending on the weights. 5. These results suggest that the mechanical properties of lungs of young rats are age-related and comparisons of ventilatory functions among those animals should be carefully done.     

Effects of azumolene on doxorubicin-induced Ca2+ release from skeletal and cardiac muscle sarcoplasmic reticulum

The mechanism of doxorubicin-induced Ca2+ release from skeletal and cardiac muscle sarcoplasmic reticulum (SR) was studied by examining the effects of azumolene (a water soluble dantrolene analog) on doxorubicin-mediated Ca2+ release and ryanodine binding. Doxorubicin induced a rapid Ca2+ release from both skeletal and cardiac SR in a similar concentration range (EC50 = 5-10 microM). Maximal doxorubicin-induced Ca2+ release was seen at 2 and 0.2 microM Ca2+ for skeletal and cardiac SR, respectively. Addition of 400 microM azumolene caused approx. 30% inhibition of doxorubicin-induced Ca2+ release from both skeletal and cardiac SR; skeletal SR had significantly higher sensitivity to azumolene than cardiac SR. In the presence of Ca2+, doxorubicin increased [3H]ryanodine binding to both skeletal and cardiac SR; whereas in the absence of Ca2+, doxorubicin led to significant ryanodine binding to skeletal SR, but not to cardiac SR. In both types of SR, doxorubicin-activated, but not Ca2+ activated ryanodine binding was inhibited by azumolene. Azumolene sensitivity for inhibition of doxorubicin-activated ryanodine binding was much higher in skeletal SR than cardiac SR, consistent with the results for effects of azumolene on Ca2+ release. Our results are consistent with the possibility that azumolene inhibits doxorubicin binding by direct competition for the drug receptor(s).     

Preparation of reconstituted acetylcholine receptor membranes suitable for AFM imaging of lipid–protein interactions

The nicotinic acetylcholine receptor (nAChR) has been reconstituted in POPC vesicles at high lipid–protein (L/P) ratios for the preparation of supported lipid bilayers with a low protein density for studies of protein–lipid interactions using atomic force microscopy (AFM). Initial reconstitutions using a standard dialysis method with bulk L/P ratios ranging from 20:1 to 100:1 (w/w) gave heterogeneous samples that contained both empty vesicles and proteoliposomes with a range of L/P ratios. This is problematic because empty vesicles adsorb and rupture to form bilayer patches more rapidly than do protein-rich vesicles, resulting in the loss of protein during sample washing. Although it was not possible to find reconstitution conditions that gave homogeneous populations of vesicles with high L/P ratios, an additional freeze–thaw cycle immediately after dialysis did reproducibly yield a fraction of proteoliposomes with L/P ratios above 100:1. These proteoliposomes were separated by sucrose gradient centrifugation and used to prepare supported bilayers with well-separated individual receptors and minimal adsorbed proteoliposomes. AFM images of such samples showed many small features protruding from the bilayer surface. These features range in height from 1 to 5 nm, consistent with the smaller intracellular domain of the protein exposed, and have lateral dimensions consistent with an individual receptor. Some bilayers with reconstituted protein also had a small fraction of higher features that are assigned to nAChR with the larger extracellular domain exposed and showed evidence for aggregation to give dimers or small oligomers. This work demonstrates the importance of using highly purified reconstituted membranes with uniform lipid–protein ratios for AFM studies of integral membrane protein–lipid interactions.     

CD4+CD25(int) T cells in inflammatory diseases refractory to treatment with glucocorticoids

Up to 30% of patients with autoimmune, allergic, and lymphoproliferative diseases are refractory to glucocorticoid therapy. The present study was undertaken to investigate whether such steroid resistance (SR) is limited to a subpopulation of CD4(+) T cells and, as IL-2 is a putative driver of SR, whether T cell SR is associated with CD25 expression. We show that SR patients have a characteristic subgroup of activated CD4(+) T cells that continue to proliferate despite exposure to high-dose Dexamethasone (Dex), demonstrate that CD4(+)CD25(-) cells are exquisitely sensitive to Dex whereas CD4(+)CD25(int) cells are highly SR, and further find that the combination of an anti-CD25 mAb with Dex enhances suppression of T cell proliferation compared with each agent alone. We therefore conclude that SR is not a general property of all lymphocytes but resides in T cell subpopulations, which are prevalent in SR patients and express intermediary levels of CD25. As a result, we propose a new paradigm for SR disease in which glucocorticoid therapy positively selects SR cells, generating a population of drug-resistant lymphocytes that perpetuate on-going inflammation.     

Effects of thoracic gas compression on maximal and partial flow-volume maneuvers

Airway hysteresis can be evaluated by comparing maximal (MEFV) and partial (PEFV) expiratory flow-volume curves. The maneuvers are often obtained from pulmonary function systems that are subject to gas-compression artifacts. Because gas-compression artifacts might differentially affect PEFV vs. MEFV curves, we simultaneously obtained MEFV and PEFV curves by use of a spirometer and a volume-displacement plethysmograph (a method not subject to gas-compression artifacts) in normal and asthmatic subjects. Plethysmographic flow rates exceeded spirometric flow rates on all MEFV and PEFV maneuvers. When maximal flow exceeded partial flow (or vice versa) in the plethysmograph, the same result was virtually always observed for spirometric measurements. Alveolar pressure (PA) was higher on MEFV than on PEFV maneuvers in asthmatic subjects; comparisons between PA (on PEFV and MEFV maneuvers) in normal subjects varied at different lung volumes. Ratios of Vmax on PEFV maneuvers to Vmax on MEFV maneuvers (Vmax-p/Vmax-c) obtained from a volume-displacement plethysmograph differ quantitatively from ratios determined in systems subject to gas-compression artifacts; qualitatively, however, failure to account for thoracic gas compression ordinarily will not influence the ability to identify airway hysteresis (or lack thereof) by use of Vmax-p-to-Vmax-c ratios.     

Effects of airway tone and volume history on maximal expiratory flow in asthma

We assessed the difference between isovolumic maximal expiratory flows (Vmax) using maneuvers begun at mid-lung volumes, so-called partial expiratory flow-volume curves (P), vs. those begun at full inflation, so-called maximal expiratory flow-volume curves (M), in 10 asthmatic subjects before and following obstruction induced by isocapnic hyperpnea with cold air and before and after bronchodilation with a beta-agonist or antimuscarinic agent. Volume history effects were quantitated as an M-to-P ratio of Vmax at 30% vital capacity (M/P V30). Although M/P V30 was variable among patients at base line, there was a uniform increase in M/P V30 during constriction and a consistent decrease below base line after dilation. Blunting of induced obstruction with beta-agonists also diminished the increase in M/P V30. Antimuscarinics, despite equivalent bronchodilation, failed to alter the degree of obstruction induced by cold air or the increase in M/P V30 seen during obstruction. The level of airway tone, as indicated by specific resistance, related directly to the M/P V30. We conclude that the response of the asthmatic lung to a deep inhalation is relatively predictable when acute changes in airway tone are produced.     

Differential responses of carbon‐degrading enzyme activities to warming: Implications for soil respiration

Extracellular enzymes catalyze rate‐limiting steps in soil organic matter decomposition, and their activities (EEAs) play a key role in determining soil respiration (SR). Both EEAs and SR are highly sensitive to temperature, but their responses to climate warming remain poorly understood. Here, we present a meta‐analysis on the response of soil cellulase and ligninase activities and SR to warming, synthesizing data from 56 studies. We found that warming significantly enhanced ligninase activity by 21.4% but had no effect on cellulase activity. Increases in ligninase activity were positively correlated with changes in SR, while no such relationship was found for cellulase. The warming response of ligninase activity was more closely related to the responses of SR than a wide range of environmental and experimental methodological factors. Furthermore, warming effects on ligninase activity increased with experiment duration. These results suggest that soil microorganisms sustain long‐term increases in SR with warming by gradually increasing the degradation of the recalcitrant carbon pool.     

A model of propagating calcium-induced calcium release mediated by calcium diffusion

The effect of sudden local fluctuations of the free sarcoplasmic [Ca++]i in cardiac cells on calcium release and calcium uptake by the sarcoplasmic reticulum (SR) was calculated with the aid of a simplified model of SR calcium handling. The model was used to evaluate whether propagation of calcium transients and the range of propagation velocities observed experimentally (0.05-15 mm s(-1)) could be predicted. Calcium fluctuations propagate by virtue of focal calcium release from the SR, diffusion through the cytosol (which is modulated by binding to troponin and calmodulin and sequestration by the SR), and subsequently induce calcium release from adjacent release sites of the SR. The minimal and maximal velocities derived from the simulation were 0.09 and 15 mm s(-1) respectively. The method of solution involved writing the diffusion equation as a difference equation in the spatial coordinates. Thus, coupled ordinary differential equations in time with banded coefficients were generated. The coupled equations were solved using Gear's sixth order predictor-corrector algorithm for stiff equations with reflective boundaries. The most important determinants of the velocity of propagation of the calcium waves were the diastolic [Ca++]i, the rate of rise of the release, and the amount of calcium released from the SR. The results are consistent with the assumptions that calcium loading causes an increase in intracellular calcium and calcium in the SR, and an increase in the amount and rate of calcium released. These two effects combine to increase the propagation velocity at higher levels of calcium loading.     

Computerized method for analyzing maximum and partial expiratory flow-volume curves.

Computerized instrumentation and software have been developed to obtain maximum expiratory flow-volume (MEFV) and partial expiratory flow-volume (PEFV) curves. The computerized system calculates and prints out the flow at 25% and 40% of control vital capacity (VC), the expiratory volume, peak expiratory flow rate and expiratory volume at one second (FEV1) divided by VC, the latter expressed as a percent. The flow-volume curves can be displayed on an oscilloscope or plotter and stored on magnetic tape. A pilot study was completed to demonstrate the reliability and validity of the data obtained.     

Model of intracellular calcium cycling in ventricular myocytes

We present a mathematical model of calcium cycling that takes into account the spatially localized nature of release events that correspond to experimentally observed calcium sparks. This model naturally incorporates graded release by making the rate at which calcium sparks are recruited proportional to the whole cell L-type calcium current, with the total release of calcium from the sarcoplasmic reticulum (SR) being just the sum of local releases. The dynamics of calcium cycling is studied by pacing the model with a clamped action potential waveform. Experimentally observed calcium alternans are obtained at high pacing rates. The results show that the underlying mechanism for this phenomenon is a steep nonlinear dependence of the calcium released from the SR on the diastolic SR calcium concentration (SR load) and/or the diastolic calcium level in the cytosol, where the dependence on diastolic calcium is due to calcium-induced inactivation of the L-type calcium current. In addition, the results reveal that the calcium dynamics can become chaotic even though the voltage pacing is periodic. We reduce the equations of the model to a two-dimensional discrete map that relates the SR and cytosolic concentrations at one beat and the previous beat. From this map, we obtain a condition for the onset of calcium alternans in terms of the slopes of the release-versus-SR load and release-versus-diastolic-calcium curves. From an analysis of this map, we also obtain an understanding of the origin of chaotic dynamics.