Short-Term Rhythms of the Iorespiratory Sy Significance in Neonatology

Latent disturbances in the control of respiration and heart rate (HR) may be important factors in the pathogenesis of life-threatening events during infancy. A method of determining the control of the autonomic nervous system functions involves the analysis of time-dependent ultradian changes of its parameters. The breathing signal and HR variability contain rhythmic components that are generated within the cardiorespiratory network of the brain stem, through reflexes, and by feedback mechanisms. The analysis of these components may provide insights into the functioning of the cardiorespiratory control system. The prominence and precision of the rhythms are correlated with states of vigilance and underlie distinct development during the first months of life. The results of studies on infants at risk (for example, for sudden infant death), with the help of statistical and spectral analysis of time series to obtain new indices, have proved to be inconsistent in their prognostic value of thus studied parameters. Recently, the importance of qualitative and quantitative assessment of the dynamic and complex behavior of time series, based on nonlinear characteristics of the control system, has been emphasized. To what extent, however, the analysis of the dynamic behavior can be utilized for clinical purposes, such as judging the prognosis of deficiencies in control, requires further study regarding physiological baselines and the possible changes resulting from pathological states.     

A formalized criterion of cardiorespiratory synchronization for assessment of dynamic changes in autonomic homeostasis

A comparative study was performed to analyze the diagnostic efficiency of a proposed cardiorespiratory coefficient (CRC) for assessment of dynamic changes in autonomic homeostasis. In the two series of the study, 45 healthy male volunteers were examined using models of emotional tension and acute graded hypoxia. Under stress-inducing conditions, both cardiorespiratory synchronization reflecting the degree of autonomic balance and the changes in the cardiorespiratory system substantially depended on the individual parameters recorded at baseline, specifically, the autonomic tone and the degree of hypoxic resistance. It was found that cardiorespiratory changes could be detected rather quickly during the first minutes of testing if assessed by the CRC changes during test exercises. This allows the accuracy and selectivity of analysis of changes in the integral parameters of autonomic regulation to be increased and the CRC to be recommended as a forecasting criterion for analyzing autonomic homeostasis and autonomic reactivity in humans.     

Dynamic regulation of heart rate during acute hypotension: new insight into baroreflex function

To examine the dynamic properties of baroreflex function, we measured beat-to-beat changes in arterial blood pressure (ABP) and heart rate (HR) during acute hypotension induced by thigh cuff deflation in 10 healthy subjects under supine resting conditions and during progressive lower body negative pressure (LBNP). The quantitative, temporal relationship between ABP and HR was fitted by a second-order autoregressive (AR) model. The frequency response was evaluated by transfer function analysis. Results: HR changes during acute hypotension appear to be controlled by an ABP error signal between baseline and induced hypotension. The quantitative relationship between changes in ABP and HR is characterized by a second-order AR model with a pure time delay of 0.75 s containing low-pass filter properties. During LBNP, the change in HR/change in ABP during induced hypotension significantly decreased, as did the numerator coefficients of the AR model and transfer function gain. Conclusions: 1) Beat-to-beat HR responses to dynamic changes in ABP may be controlled by an error signal rather than directional changes in pressure, suggesting a "set point" mechanism in short-term ABP control. 2) The quantitative relationship between dynamic changes in ABP and HR can be described by a second-order AR model with a pure time delay. 3) The ability of the baroreflex to evoke a HR response to transient changes in pressure was reduced during LBNP, which was due primarily to a reduction of the static gain of the baroreflex.     

Cardiorespiratory Information Dynamics during Mental Arithmetic and Sustained Attention

An analysis of cardiorespiratory dynamics during mental arithmetic, which induces stress, and sustained attention was conducted using information theory. The information storage and internal information of heart rate variability (HRV) were determined respectively as the self-entropy of the tachogram, and the self-entropy of the tachogram conditioned to the knowledge of respiration. The information transfer and cross information from respiration to HRV were assessed as the transfer and cross-entropy, both measures of cardiorespiratory coupling. These information-theoretic measures identified significant nonlinearities in the cardiorespiratory time series. Additionally, it was shown that, although mental stress is related to a reduction in vagal activity, no difference in cardiorespiratory coupling was found when several mental states (rest, mental stress, sustained attention) are compared. However, the self-entropy of HRV conditioned to respiration was very informative to study the predictability of RR interval series during mental tasks, and showed higher predictability during mental arithmetic compared to sustained attention or rest.     

Co-operative components, spatial localization and oscillatory cellular dynamics

The consequences of spatial localization on the dynamic behavior of catalytic components whose outputs control one another's activities but are not the substrates of those activities are explored. For some simple model systems we naturally find limit cycle oscillations, multiple asymptotic states, dependence of the dynamic states upon spatial configuration of the components, and complex switching between dynamic states as a function of external perturbation to the system.     

Heart Rate Dynamics after Combined Strength and Endurance Training in Middle-Aged Women: Heterogeneity of Responses

The loss of complexity in physiological systems may be a dynamical biomarker of aging and disease. In this study the effects of combined strength and endurance training compared with those of endurance training or strength training alone on heart rate (HR) complexity and traditional HR variability indices were examined in middle-aged women. 90 previously untrained female volunteers between the age of 40 and 65 years completed a 21 week progressive training period of either strength training, endurance training or their combination, or served as controls. Continuous HR time series were obtained during supine rest and submaximal steady state exercise. The complexity of HR dynamics was assessed using multiscale entropy analysis. In addition, standard time and frequency domain measures were also computed. Endurance training led to increases in HR complexity and selected time and frequency domain measures of HR variability (P<0.01) when measured during exercise. Combined strength and endurance training or strength training alone did not produce significant changes in HR dynamics. Inter-subject heterogeneity of responses was particularly noticeable in the combined training group. At supine rest, no training-induced changes in HR parameters were observed in any of the groups. The present findings emphasize the potential utility of endurance training in increasing the complex variability of HR in middle-aged women. Further studies are needed to explore the combined endurance and strength training adaptations and possible gender and age related factors, as well as other mechanisms, that may mediate the effects of different training regimens on HR dynamics.     

Personality development in psychotherapy: a synergetic model of state-trait dynamics

Theoretical models of psychotherapy not only try to predict outcome but also intend to explain patterns of change. Studies showed that psychotherapeutic change processes are characterized by nonlinearity, complexity, and discontinuous transitions. By this, theoretical models of psychotherapy should be able to reproduce these dynamic features. Using time series derived from daily measures through internet-based real-time monitoring as empirical reference, we earlier presented a model of psychotherapy which includes five state variables and four trait variables. In mathematical terms, the traits modulate the shape of the functions which define the nonlinear interactions between the variables (states) of the model. The functions are integrated into five coupled nonlinear difference equations. In the present paper, we model how traits (dispositions or competencies of a person) can continuously be altered by new experiences and states (cognition, emotion, behavior). Adding equations that link states to traits, this model not only describes how therapeutic interventions modulate short-term change and fluctuations of psychological states, but also how these can influence traits. Speaking in terms of Synergetics (theory of self-organization in complex systems), the states correspond to the order parameters and the traits to the control parameters of the system. In terms of psychology, trait dynamics is driven by the states—i.e., by the concrete experiences of a client—and creates a process of personality development at a slower time scale than that of the state dynamics (separation of time scales between control and order parameter dynamics).     

Adjustment of the human respiratory system to increased upper airway resistance during sleep

A cardiorespiratory model incorporating control of the human upper airway during sleep is described. Most previous models have not considered the possibility that the upper airway could be a limiting factor for gas exchange. Our model was developed to also predict certain pathophysiological phenomena in the cardiorespiratory system that characterize heavy snoring or sleep apnea. We started by adapting our collapsible upper airway model to include the impact of nasal passage and larynx, and extended the model with equations for gas exchange in the lungs. A feedback loop both to the respiratory pump and the upper airway dilator muscles was included. The model enabled successful breath-by-breath simulations of obstructive events of the upper airway. Although the model incorporates several physiologically relevant components of the system, the simulation results suggest that only few parameters suffice to predict the key adjustments that the cardiorespiratory system is known to make in patients with heavy snoring.     

Metabolic modeling of Saccharomyces cerevisiae using the optimal control of homeostasis: a cybernetic model definition

A model is presented to describe the observed behavior of microorganisms that aim at metabolic homeostasis while growing and adapting to their environment in an optimal way. The cellular metabolism is seen as a network with a multiple controller system with both feedback and feedforward control, i.e., a model based on a dynamic optimal metabolic control. The dynamic network consists of aggregated pathways, each having a control setpoint for the metabolic states at a given growth rate. This set of strategies of the cell forms a true cybernetic model with a minimal number of assumptions. The cellular strategies and constraints were derived from metabolic flux analysis using an identified, biochemically relevant, stoichiometry matrix derived from experimental data on the cellular composition of continuous cultures of Saccharomyces cerevisiae. Based on these data a cybernetic model was developed to study its dynamic behavior. The growth rate of the cell is determined by the structural compounds and fluxes of compounds related to central metabolism. In contrast to many other cybernetic models, the minimal model does not consist of any assumed internal kinetic parameters or interactions. This necessitates the use of a stepwise integration with an optimization of the fluxes at every time interval. Some examples of the behavior of this model are given with respect to steady states and pulse responses. This model is very suitable for describing semiquantitatively dynamics of global cellular metabolism and may form a useful framework for including structured and more detailed kinetic models.     

On the modeling and interpretation of oxygen uptake kinetics from ramp work rate tests

Ramp work rate tests have been used to estimate aerobic parameters in exercise stress testing. Previous studies have suggested an assumption of a linear dynamic system for O2 uptake kinetics. The implication is that model parameters estimated from ramp tests should be similar to those estimated from other dynamic tests. In nine healthy subjects, we found that model parameters used to characterize O2 consumption ramp data were not consistent with those used to characterize step data, when the comparison was made on a subject-to-subject basis. Furthermore the ramp data model parameter values were highly dependent (P less than 0.0001) on the ramp slope. A linear dynamic system interpretation of the ramp data model does not appear to be appropriate, suggesting that caution is needed in the interpretation of ramp data aerobic parameters. The data may be better described by nonlinear or higher order function. Ramp exercise testing is not suitable for assessing dynamic control properties of the cardiorespiratory response to exercise.     

The respiratory consequences of feeding in amphibians and reptiles

Many ectothermic vertebrates ingest very large meals at infrequent intervals. The digestive processes associated with these meals, often coupled with an extensive hypertrophy of the gastrointestinal organs, are energetically expensive and metabolic rate, therefore, increases substantially after feeding (specific dynamic action, SDA). Here, we review the cardio-respiratory consequences of SDA in amphibians and reptiles. For some snakes, the increased oxygen uptake during SDA is of similar magnitude to that of muscular exercise, and the two physiological states, therefore, exert similar and profound demands on oxygen transport by the cardiorespiratory systems. In several species, SDA is attended by increases in heart rate and overall systemic blood flows, but changes in blood flow distribution remain to be investigated. In snakes, the regulation of heart rate appears to involve a non-adrenergic-non-cholinergic mechanism, which may be a regulatory peptide released from the gastrointestinal system during digestion. Digestion is also associated with a net acid secretion to the stomach that causes an increase in plasma HCO3- concentration (the 'alkaline tide'). Experiments on chronically cannulated amphibians and reptiles, show that this metabolic alkalosis is countered by an increased P(CO2), so that the change in arterial pH is reduced. This respiratory compensation of arterial pH is accomplished through a reduction in ventilation relative to metabolism, but the estimated reductions in lung P(O2) are relatively small. The SDA response is also associated with haematological changes, but large interspecific differences exist. The studies on cardiorespiratory responses to digestion may allow for a further understanding of the physiological and structural constraints that limits the ability of reptiles and amphibians to sustain high metabolic rates.     

Inhibition of heart rate variability during sleep in humans by 6700 K pre-sleep light exposure

Two different spectral analyses of heart rate (HR) variability (HRV) were performed on seven young male subjects to evaluate the effects of different color temperatures of light exposure (6700 K, 5000 K, 3000 K) before sleep on cardiac vagal activity. In investigating HRV, we used an ordinary fast Fourier transform (FFT) and coarse graining spectral analysis (CGSA), which selectively extracts random fractal components from a given time series. The results showed that suppressions of HR during sleep after 6700 K light exposure were more inhibited than the other two lighting conditions. Increases in high-frequency (HF) components of HRV during sleep were also inhibited by 6700 K pre-sleep lighting. These results indicate that pre-sleep exposure to light of a higher color temperature may inhibit the enhancement of cardiac vagal activity during sleep. Moreover, significant HF alterations were shown in fractal-free HF (not in ordinary HF) components by CGSA. Because the HF component originates from respiratory sinus arrhythmia with periodical fluctuations, CGSA may be an appropriate approach for HRV evaluation during sleep.     

Chaotic nonlinear dynamics of prestimulus EEG and the structure of auditory evoked potentials

The dynamic pattern of a prestimulus EEG was studied by the fractal analysis technique. The character of this pattern was shown to affect the structure of cortical auditory evoked potentials (EP). A monoperiodic pattern in the control frequency band was accompanied by a formation of several new attractors with low dimensionality attributed to simultaneous functioning of several weakly connected dynamic systems with an indefinite trend of a leading process. During development of multiperiodic processes in EEG segments, the fractal analysis revealed a tendency for a formation of a single complex dynamic system with high dimensionality of the attractor. Characteristic changes in parameters of the primary and middle-latency EP components, their correlation, and factor models are related with the character of nonlinear patterns. Analysis of variance revealed the most effective role of pattern changes in the EEG alpha control band. The structure of combination of parameters of the EP primary components into connected complexes depends on the nonlinear prestimulus EEG patterns. The above predictors determine differently directed and differently pronounced changes in parameters of the middle-latency positive waves.     

A skill-based conditioning games approach to metabolic conditioning for elite rugby football players

The purpose of this study was to evaluate changes in endurance fitness of elite-level rugby union players (n = 35) undertaking skill-based conditioning games for a 9-week preseason training period. Metabolic conditioning was conducted exclusively in the form of skill-based conditioning games in conjunction with heart rate (HR) telemetry. Two markers of cardiorespiratory fitness were assessed at weekly intervals via the recording of HR responses to an intermittent multistage shuttle test. Significant differences post-training were observed for the percentage of maximal HR (% HRmax) reached during the final test stage and the percentage of HR recovery (% HR recovery) from the end of the final stage to the end of the final 1-minute rest period. Significant improvements were demonstrated for % HR recovery at week 7 (p < 0.05) and week 9 (p < 0.01), and % HRmax in the final test stage was significantly lower at weeks 4, 5, and 7 (p < 0.05) and week 9 (p < 0.01). Further improvements from mid-preseason to the end of the preseason training period were observed for % HR recovery scores in week 8 (p < 0.01) and week 9 (p = 0.012) and for % HRmax reached in the final test stage at week 9 (p < 0.05). These results indicate skill-based conditioning games were successful at improving markers of cardiorespiratory endurance for the duration of a 9-week training period in the elite-level professional rugby union players studied. The HR monitoring was demonstrated to be an effective and practical means of quantifying intensity in the conditioning games format and of tracking changes in cardiorespiratory fitness.     

Dynamic simulation and steady-state analysis of a bentho-pelagic coupled ecosystem under different simulation scenarios

Trophic dynamic studies of aquatic ecosystem are usually based on either dynamic or state-state models of energy and material transfer, comprising of a few functional tropic levels such as phytoplankton, primary consumers and secondary consumers. However, though benthic components (and production) are important parts of an aquatic ecosystem, they are often omitted from such approaches.Considering this inadequacy of previous simulation approaches, the aim of the present study, is to investigate the effect of benthic components on overall system dynamics. The objectives of this study are (a) to develop and study a dynamic model for Kakinada bay to understand how benthic components interact with others under a fixed set of conditions, and to find the sensitive parameters for such a benthic-pelagic coupled system (b) to study the effects of different biomasses of benthic components via scenario analysis and (c) to identify the different equilibrium states that a benthic-pelagic coupled system such as Kakinada Bay may reach.On running the simulation for a period of 5000 days (∼10+ years), the system was not seen to reach any equilibrium state; however, all state variables considered in the model were seen to coexist even at the end of this period. While zooplankton was found to be the most sensitive state variable, parameters directly associated with benthic components were indicated to be the most sensitive. Microphytobenthos was negatively correlated with phytoplankton which was corroborated by results of the perturbation scenario analyses as well. Only one equilibrium state – microphytobenthos dominated steady state was found for Kakinada Bay system when certain parameters were slightly changed.     

A Mathematical Model of Uterine Dynamics and its Application to Human Parturition

We have developed a simple mathematical model with three physiologically significant states to describe the changes in intrauterine pressure associated with a contraction during human parturition. The myometrium is modelled as a set of smooth muscle cells, each of which is in one of three states (quiescent, contracted, refractory) at a given time. These states are occupied according to a cycle governed by three temporal parameters. The solutions of the equations describing the model show an oscillatory behavior for particular values of these parameters, which is very similar to the time dependant development of intrauterine pressure during labor. Due to its non-linear terms, our model could lead to chaotic oscillations (in the mathematical sense), whose clinical counterpart may occur in cases of dystocia. Despite its simplicity, this model appears to be a useful guide to further investigations of the oscillatory behavior of the myometrium, or other smooth muscles, in normal and pathological situations.     

Structural time series models with feedback mechanisms

Structural time series models have applications in many different fields such as biology, economics, and meteorology. A structural times series model can be represented as a state-space model where the states of the system represent the unobserved components and the structural parameters have clear interpretations. This paper introduces a class of structural time series models that incorporate feedback from the latent components of the history. An iterative procedure is proposed for estimation. These models allow flexible and robust feedback mechanisms, have clear interpretations, and have a computationally efficient estimation procedure. They are applied to hormone data to characterize hormone secretion and to explore a potential feedback mechanism.     

Cardiorespiratory interactions during resistive load breathing

The addition to the respiratory system of a resistive load results in breathing pattern changes and in negative intrathoracic pressure increases. The aim of this study was to use resistive load breathing as a stimulus to the cardiorespiratory interaction and to examine the extent of the changes in heart rate variability (HRV) and respiratory sinus arrhythmia (RSA) in relation to the breathing pattern changes. HRV and RSA were studied in seven healthy subjects where four resistive loads were applied in a random order during the breath and 8-min recording made in each condition. The HRV spectral power components were computed from the R-R interval sequences, and the RSA amplitude and phase were computed from the sinusoid fitting the instantaneous heart rate within each breath. Adding resistive loads resulted in 1) increasing respiratory period, 2) unchanging heart rate, and 3) increasing HRV and changing RSA characteristics. HRV and RSA characteristics are linearly correlated to the respiratory period. These modifications appear to be linked to load-induced changes in the respiratory period in each individual, because HRV and RSA characteristics are similar at a respiratory period obtained either by loading or by imposed frequency breathing. The present results are discussed with regard to the importance of the breathing cycle duration in these cardiorespiratory interactions, suggesting that these interactions may depend on the time necessary for activation and dissipation of neurotransmitters involved in RSA.