Enhanced multiresolution wavelet analysis of complex dynamics in nonlinear systems.

Multiresolution wavelet analysis (MWA) is a powerful data processing tool that provides a characterization of complex signals over multiple time scales. Typically, the standard deviations of wavelet coefficients are computed depending on the resolution level and such quantities are used as measures for diagnosing different types of system behavior. To enhance the capabilities of this tool, we propose a combination of MWA with detrended fluctuation analysis (DFA) of detail wavelet coefficients. We find that such an MWA&DFA approach is capable of revealing the correlation features of wavelet coefficients in independent ranges of scales, which provide more information about the complex organization of datasets compared to variances or similar statistical measures of the standard MWA. Using this approach, we consider changes in the dynamics of coupled chaotic systems caused by transitions between different types of complex oscillations. We also demonstrate the potential of the MWA&DFA method for characterizing different physiological conditions by analyzing the electrical brain activity in mice.

Extended detrended fluctuation analysis of electroencephalograms signals during sleep and the opening of the blood-brain barrier.

Detrended fluctuation analysis (DFA) is widely used to characterize long-range power-law correlations in complex signals. However, it has restrictions when nonstationarity is not limited only to slow variations in the mean value. To improve the characterization of inhomogeneous datasets, we have proposed the extended DFA (EDFA), which is a modification of the conventional method that evaluates an additional scaling exponent to take into account the features of time-varying nonstationary behavior. Based on EDFA, here, we analyze rat electroencephalograms to identify specific changes in the slow-wave dynamics of brain electrical activity associated with two different conditions, such as the opening of the blood-brain barrier and sleep, which are both characterized by the activation of the brain drainage function. We show that these conditions cause a similar reduction in the scaling exponents of EDFA. Such a similarity may represent an informative marker of fluid homeostasis of the central nervous system.

Characterizing scaling properties of complex signals with missed data segments using the multifractal analysis.

The scaling properties of complex processes may be highly influenced by the presence of various artifacts in experimental recordings. Their removal produces changes in the singularity spectra and the Hölder exponents as compared with the original artifacts-free data, and these changes are significantly different for positively correlated and anti-correlated signals. While signals with power-law correlations are nearly insensitive to the loss of significant parts of data, the removal of fragments of anti-correlated signals is more crucial for further data analysis. In this work, we study the ability of characterizing scaling features of chaotic and stochastic processes with distinct correlation properties using a wavelet-based multifractal analysis, and discuss differences between the effect of missed data for synchronous and asynchronous oscillatory regimes. We show that even an extreme data loss allows characterizing physiological processes such as the cerebral blood flow dynamics.

Cardiovascular System's States under Stress: Reviews of Methods.

A comparative analysis is made of various methods for processing electrocardiograms and RR-interval sequences. This analysis was carried out by using standard nonlinear-dynamics algorithms and methods. Apart from that, we assessed the expediency of using a number of characteristics to classify the cardiovascular system's state under stress.

Characterization of the chaos-hyperchaos transition based on return times.

We discuss the problem of the detection of hyperchaotic oscillations in coupled nonlinear systems when the available information about this complex dynamical regime is very limited. We demonstrate the ability of diagnosing the chaos-hyperchaos transition from return times into a Poincaré section and show that an appropriate selection of the secant plane allows a correct estimation of two positive Lyapunov exponents (LEs) from even a single sequence of return times. We propose a generalized approach for extracting dynamics from point processes that allows avoiding spurious identification of the dynamical regime caused by artifacts. The estimated LEs are nearly close to their expected values if the second positive LE is essentially different from the largest one. If both exponents become nearly close, an underestimation of the second LE may be obtained. Nevertheless, distinctions between chaotic and hyperchaotic regimes are clearly possible.

Quantifying chaotic dynamics from integrate-and-fire processes.

Characterizing chaotic dynamics from integrate-and-fire (IF) interspike intervals (ISIs) is relatively easy performed at high firing rates. When the firing rate is low, a correct estimation of Lyapunov exponents (LEs) describing dynamical features of complex oscillations reflected in the IF ISI sequences becomes more complicated. In this work we discuss peculiarities and limitations of quantifying chaotic dynamics from IF point processes. We consider main factors leading to underestimated LEs and demonstrate a way of improving numerical determining of LEs from IF ISI sequences. We show that estimations of the two largest LEs can be performed using around 400 mean periods of chaotic oscillations in the regime of phase-coherent chaos. Application to real data is discussed.

Multiresolution analysis of pathological changes in cerebral venous dynamics in newborn mice with intracranial hemorrhage: adrenorelated vasorelaxation.

Intracranial hemorrhage (ICH) is the major problem of modern neonatal intensive care. Abnormalities of cerebral venous blood flow (CVBF) can play a crucial role in the development of ICH in infants. The mechanisms underlying these pathological processes remain unclear; however it has been established that the activation of the adrenorelated vasorelaxation can be an important reason. Aiming to reach a better understanding of how the adrenodependent relaxation of cerebral veins contributes to the development of ICH in newborns, we study here the effects of pharmacological stimulation of adrenorelated dilation of the sagittal sinus by isoproterenol on the cerebral venous hemodynamics. Our study is performed in newborn mice at different stages of ICH using the laser speckle contrast imaging and wavelet analysis of the vascular dynamics of CVBF. We show that the dilation of the sagittal sinus with the decreased velocity of blood flow presides to the stress-induced ICH in newborn mice. These morphofunctional vascular changes are accompanied by an increased variance of the wavelet-coefficients in the areas of endothelial and non-endothelial (KATP-channels activity of vascular muscle) sympathetic components of the CVBF variability. Changes in the cerebral venous hemodynamics at the latent stage of ICH are associated with a high responsiveness of the sagittal sinus to isoproterenol quantifying by wavelet-coefficients related to a very slow region of the frequency domain. The obtained results certify that a high activation of the adrenergic-related vasodilatory responses to severe stress in newborn mice can be one of the important mechanisms underlying the development of ICH. Thus, the venous insufficiency with the decreased blood outflow from the brain associated with changes in the endothelial and the sympathetic components of CVBF-variability can be treated as prognostic criteria for the risk of ICH during the first days after birth.

Analysis of blood pressure dynamics in male and female rats using the continuous wavelet transform.

We study gender-related particularities in cardiovascular responses to stress and nitric oxide (NO) deficiency in rats using HR, mean arterial pressure (MAP) and a proposed wavelet-based approach. Blood pressure dynamics is analyzed: (1) under control conditions, (2) during immobilization stress and recovery and (3) during nitric oxide blockade by N(G)-nitro-L-arginine-methyl ester (L-NAME). We show that cardiovascular sensitivity to stress and NO deficiency depends upon gender. Actually, in females the chronotropic effect of stress is more pronounced, while the pressor effect is weakened compared with males. We conclude that females demonstrate more favorable patterns of cardiovascular responses to stress and more effective NO control of cardiovascular activity than males.

The effect of L-NAME on intra- and inter-nephron synchronization.

Kidney autoregulation involves complicated intra- and inter-nephron synchronization phenomena among oscillatory modes produced, respectively, by the tubuloglomerular feedback (TGF) mechanism and by the myogenic regulation of the afferent arteriolar blood flow. The present study aims at examining to what extent these phenomena are reflected in the overall blood flow to the kidney and how they are affected by intravenous administration of nitro-l-arginine-methyl-ester (L-NAME), a potent NO synthesis inhibitor. Wavelet analysis is applied to detect rhythmic activity at the level of the renal artery and compare the observed fluctuations with blood flow variations recorded from efferent arterioles of individual nephrons. We show that administration of L-NAME increases the gain in both the TGF and the myogenic oscillations, and that both normotensive and hypertensive rats demonstrate reduced stability of the various rhythms. This implies that L-NAME, besides strengthening the gain in the individual feedback mechanisms, also causes more frequent transitions among the various synchronization states. In a broader perspective the purpose of the study is to demonstrate the significance of complex dynamic phenomena in the normal regulation of physiological systems as well as in their response to drugs.

Characterizing multimode interaction in renal autoregulation.

The purpose of this paper is to demonstrate how modern statistical techniques of non-stationary time-series analysis can be used to characterize the mutual interaction among three coexisting rhythms in nephron pressure and flow regulation. Besides a relatively fast vasomotoric rhythm with a period of 5-8 s and a somewhat slower mode arising from an instability in the tubuloglomerular feedback mechanism, we also observe a very slow mode with a period of 100-200 s. Double-wavelet techniques are used to study how the very slow rhythm influences the two faster modes. In a broader perspective, the paper emphasizes the significance of complex dynamic phenomena in the normal and pathological function of physiological systems and discusses how simulation methods can help to understand the underlying biological mechanisms. At the present there is no causal explanation of the very slow mode. However, vascular oscillations with similar frequencies have been observed in other tissues.

Using wavelet analysis to detect the influence of low frequency magnetic fields on human physiological tremor.

The influence of extremely low frequency magnetic fields (ELF-MFs) on human physiological processes and, in particular, on motor activity is still not established with certainty. Using the wavelet-transform approach, changes in the characteristics of human finger micromovement are studied in the presence of a low intensity MF centred at the level of the head. Different approaches to nonstationary signal analysis involving real as well as complex wavelet functions are considered. We find evidence that ELF-MFs lead to more regular postural tremor and more homogeneous energy distribution.

Unraveling cell processes: interference imaging interwoven with data analysis.

The paper presents results on the application of interference microscopy and wavelet-analysis for cell visualization and studies of cell dynamics. We demonstrate that interference imaging of erythrocytes can reveal reorganization of the cytoskeleton and inhomogenity in the distribution of hemoglobin, and that interference imaging of neurons can show intracellular compartmentalization and submembrane structures. We investigate temporal and spatial variations of the refractive index for different cell types: isolated neurons, mast cells and erythrocytes. We show that the refractive dynamical properties differ from cell type to cell type and depend on the cellular compartment. Our results suggest that low frequency variations (0.1-0.6 Hz) result from plasma membrane processes and that higher frequency variations (20-26 Hz) are related to the movement of vesicles. Using double-wavelet analysis, we study the modulation of the 1 Hz rhythm in neurons and reveal its changes under depolarization and hyperpolarization of the plasma membrane. We conclude that interference microscopy combined with wavelet analysis is a useful technique for non-invasive cell studies, cell visualization, and investigation of plasma membrane properties.

Interference microscopy under double-wavelet analysis: a new approach to studying cell dynamics.

This Letter combines a novel experimental approach to the study of intracellular processes with a newly developed technique for multimode time-series analysis. Experiments are performed on isolated pond snail (Lymnaea stagnalis) neurons. Local variations in the cellular refractive index as detected by laser interference microscopy are related to the processes in the cell. A wavelet analysis shows the presence of several identifiable modes in the membrane and intracellular dynamics, and a double-wavelet analysis reveals nonlinear interactions between the regulatory processes in the form of mutual frequency and amplitude modulations.

Double-wavelet approach to studying the modulation properties of nonstationary multimode dynamics.

On the basis of double-wavelet analysis, the paper proposes a method to study interactions in the form of frequency and amplitude modulation in nonstationary multimode data series. Special emphasis is given to the problem of quantifying the strength of modulation for a fast signal by a coexisting slower dynamics and to its physiological interpretation. Application of the approach is demonstrated for a number of model systems, including a model that generates chaotic dynamics. The approach is then applied to proximal tubular pressure data from rat nephrons in order to estimate the degree to which the myogenic dynamics of the afferent arteriole is modulated by the slower tubulo-glomerular dynamics. Our analysis reveals a significantly stronger interaction between the two mechanisms in spontaneously hypertensive rats than in normotensive rats.

Double-wavelet approach to study frequency and amplitude modulation in renal autoregulation.

Biological time series often display complex oscillations with several interacting rhythmic components. Renal autoregulation, for instance, involves at least two separate mechanisms both of which can produce oscillatory variations in the pressures and flows of the individual nephrons. Using double-wavelet analysis we propose a method to examine how the instantaneous frequency and amplitude of a fast mode is modulated by the presence of a slower mode. Our method is applied both to experimental data from normotensive and hypertensive rats showing different oscillatory patterns and to simulation results obtained from a physiologically based model of the nephron pressure and flow control. We reveal a nonlinear interaction between the two mechanisms that regulate the renal blood flow in the form of frequency and amplitude modulation of the myogenic oscillations.

Bimodal oscillations in nephron autoregulation.

The individual functional unit of the kidney (the nephron) displays oscillations in its pressure and flow regulation at two different time scales: fast oscillations associated with a myogenic dynamics of the afferent arteriole, and slower oscillations arising from a delay in the tubuloglomerular feedback. We investigate the intra- and internephron entrainment of the two time scales. In addition to full synchronization, both wavelet analyses of experimental data and numerical simulations reveal a partial entrainment in which neighboring nephrons attain a state of chaotic synchronization with respect to their slow dynamics, but the fast dynamics remain desynchronized.

Comparative analysis of methods for classifying the cardiovascular system's states under stress.

A comparative analysis is made of various methods for processing electrocardiograms and RR-interval sequences. This analysis was carried out by using standard nonlinear-dynamics algorithms and methods. Apart from that, we assessed the expediency of using a number of characteristics to classify the cardiovascular system's state under stress.

Chaotic dynamics from interspike intervals.

Considering two different mathematical models describing chaotic spiking phenomena, namely, an integrate-and-fire and a threshold-crossing model, we discuss the problem of extracting dynamics from interspike intervals (ISIs) and show that the possibilities of computing the largest Lyapunov exponent (LE) from point processes differ between the two models. We also consider the problem of estimating the second LE and the possibility to diagnose hyperchaotic behavior by processing spike trains. Since the second exponent is quite sensitive to the structure of the ISI series, we investigate the problem of its computation.

Extracting dynamics from threshold-crossing interspike intervals: possibilities and limitations.

In this paper we estimate dynamical characteristics of chaotic attractors from sequences of threshold-crossing interspike intervals, and study how the choice of the threshold level (which sets the equation of a secant plane) influences the results of the numerical computations. Under quite general conditions we show that the largest Lyapunov exponent can be estimated from a series of return times to the secant plane, even in the case when some of the loops of the phase space trajectory fail to cross this plane.