Soliton model of competitive neural dynamics during binocular rivalry

Binocular rivalry is investigated in a continuum model of the primary visual cortex that includes neural excitation and inhibition, stimulus orientation preference, and spike-rate adaptation. Visual stimuli consisting of bars or edges result in localized states of neural activity described by solitons. Stability analysis shows binocular fusion gives way to binocular rivalry when the orientation difference between left-eye and right-eye stimuli destabilizes one or more solitons. The model yields conditions for binocular rivalry, and two types of competitive dynamics are found: either one soliton oscillates between two stimulus regions or two solitons fixed in position at the stimulus regions oscillate out of phase with each other.

Markov analysis of sleep dynamics

A new approach, based on a Markov transition matrix, is proposed to explain frequent sleep and wake transitions during sleep. The matrix is determined by analyzing hypnograms of 113 obstructive sleep apnea patients. Our approach shows that the statistics of sleep can be constructed via a single Markov process and that durations of all states have modified exponential distributions, in contrast to recent reports of a scale-free form for the wake stage and an exponential form for the sleep stage. Hypnograms of the same subjects, but treated with Continuous Positive Airway Pressure, are analyzed and compared quantitatively with the pretreatment ones, suggesting potential clinical applications.

Potential formulation of sleep dynamics

A physiologically based model of the mechanisms that control the human sleep-wake cycle is formulated in terms of an equivalent nonconservative mechanical potential. The potential is analytically simplified and reduced to a quartic two-well potential, matching the bifurcation structure of the original model. This yields a dynamics-based model that is analytically simpler and has fewer parameters than the original model, allowing easier fitting to experimental data. This model is first demonstrated to semiquantitatively match the dynamics of the physiologically based model from which it is derived, and is then fitted directly to a set of experimentally derived criteria. These criteria place rigorous constraints on the parameter values, and within these constraints the model is shown to reproduce normal sleep-wake dynamics and recovery from sleep deprivation. Furthermore, this approach enables insights into the dynamics by direct analogies to phenomena in well studied mechanical systems. These include the relation between friction in the mechanical system and the timecourse of neurotransmitter action, and the possible relation between stochastic resonance and napping behavior. The model derived here also serves as a platform for future investigations of sleep-wake phenomena from a dynamical perspective.

Cerebral nocardia masquerading as metastatic CNS disease in an endometrial cancer patient

Nocardia is a bacterial infection primarily originating from organic rich soil, endemic to several international geographic locations. We present the case of a 61-year-old woman previously treated for endometrial carcinoma, who three years later developed metastatic pulmonary disease and received systemic chemotherapy. After five months, she developed a large right posterior lobe lesion, suspicious for metastatic CNS disease. However, following neurosurgical resection of the lesion and infectious disease consultation, a diagnosis of nocardia was made.

Mean-field modeling of the basal ganglia-thalamocortical system. I Firing rates in healthy and parkinsonian states

Parkinsonism leads to various electrophysiological changes in the basal ganglia-thalamocortical system (BGTCS), often including elevated discharge rates of the subthalamic nucleus (STN) and the output nuclei, and reduced activity of the globus pallidus external (GPe) segment. These rate changes have been explained qualitatively in terms of the direct/indirect pathway model, involving projections of distinct striatal populations to the output nuclei and GPe. Although these populations partly overlap, evidence suggests dopamine depletion differentially affects cortico-striato-pallidal connection strengths to the two pallidal segments. Dopamine loss may also decrease the striatal signal-to-noise ratio, reducing both corticostriatal coupling and striatal firing thresholds. Additionally, nigrostriatal degeneration may cause secondary changes including weakened lateral inhibition in the GPe, and mesocortical dopamine loss may decrease intracortical excitation and especially inhibition. Here a mean-field model of the BGTCS is presented with structure and parameter estimates closely based on physiology and anatomy. Changes in model rates due to the possible effects of dopamine loss listed above are compared with experiment. Our results suggest that a stronger indirect pathway, possibly combined with a weakened direct pathway, is compatible with empirical evidence. However, altered corticostriatal connection strengths are probably not solely responsible for substantially increased STN activity often found. A lower STN firing threshold, weaker intracortical inhibition, and stronger striato-GPe inhibition help explain the relatively large increase in STN rate. Reduced GPe-GPe inhibition and a lower GPe firing threshold can account for the comparatively small decrease in GPe rate frequently observed. Changes in cortex, GPe, and STN help normalize the cortical rate, also in accord with experiments. The model integrates the basal ganglia into a unified framework along with an existing thalamocortical model that already accounts for a wide range of electrophysiological phenomena. A companion paper discusses the dynamics and oscillations of this combined system.

Modeling distributed axonal delays in mean-field brain dynamics

The range of conduction delays between connected neuronal populations is often modeled as a single discrete delay, assumed to be an effective value averaging over all fiber velocities. This paper shows the effects of distributed delays on signal propagation. A distribution acts as a linear filter, imposing an upper frequency cutoff that is inversely proportional to the delay width. Distributed thalamocortical and corticothalamic delays are incorporated into a physiologically based mean-field model of the cortex and thalamus to illustrate their effects on the electroencephalogram (EEG). The power spectrum is acutely sensitive to the width of the thalamocortical delay distribution, and more so than the corticothalamic distribution, because all input signals must travel along the thalamocortical pathway. This imposes a cutoff frequency above which the spectrum is overly damped. The positions of spectral peaks in the resting EEG depend primarily on the distribution mean, with only weak dependences on distribution width. Increasing distribution width increases the stability of fixed point solutions. A single discrete delay successfully approximates a distribution for frequencies below a cutoff that is inversely proportional to the delay width, provided that other model parameters are moderately adjusted. A pair of discrete delays together having the same mean, variance, and skewness as the distribution approximates the distribution over the same frequency range without needing parameter adjustment. Delay distributions with large fractional widths are well approximated by low-order differential equations.

Modeling the impact of impulsive stimuli on sleep-wake dynamics

A neuronal population model of the sleep-wake switch is extended to incorporate impulsive external stimuli. The model includes the mutual inhibition of the sleep-active neurons in the hypothalamic ventrolateral preoptic area (VLPO) and the wake-active monoaminergic brainstem populations (MA), as well as circadian and homeostatic drives. Arbitrary stimuli are described in terms of their relative effects on the VLPO and MA nuclei and represent perturbations on the normal sleep-wake dynamics. By separating the model's intrinsic time scales, an analytic characterization of the dynamics in a reduced model space is developed. Using this representation, the model's response to stimuli is studied, including the latency to return to wake or sleep, or to elicit a transition between the two states. Since sensory stimuli are known to excite the MA, we correspondingly investigate the model's response to auditory tones during sleep, as in clinical sleep fragmentation studies. The arousal threshold is found to vary approximately linearly with the model's total sleep drive, which includes circadian and homeostatic components. This relationship is used to reproduce the clinically observed variation of the arousal threshold across the night, which rises to a maximum near the middle of the night and decreases thereafter. In a further application of the model, time-of-night arousal threshold and body temperature variations in an experimental sleep fragmentation study are replicated. It is proposed that the shift of the extrema of these curves to a greater magnitude later in the night is due to the homeostatic impact of the frequent nocturnal disturbances. By modeling the underlying neuronal interactions, the methods presented here allow the prediction of arousal state responses to external stimuli. This methodology is fundamentally different to previous approaches that model the clinical data within a phenomenological framework. As a result, a broader understanding of how impulsive external stimuli modulate arousal is gained.

Sleep deprivation in a quantitative physiologically based model of the ascending arousal system

A physiologically based quantitative model of the human ascending arousal system is used to study sleep deprivation after being calibrated on a small set of experimentally based criteria. The model includes the sleep-wake switch of mutual inhibition between nuclei which use monoaminergic neuromodulators, and the ventrolateral preoptic area. The system is driven by the circadian rhythm and sleep homeostasis. We use a small number of experimentally derived criteria to calibrate the model for sleep deprivation, then investigate model predictions for other experiments, demonstrating the scope of application. Calibration gives an improved parameter set, in which the form of the homeostatic drive is better constrained, and its weighting relative to the circadian drive is increased. Within the newly constrained parameter ranges, the model predicts repayment of sleep debt consistent with experiment in both quantity and distribution, asymptoting to a maximum repayment for very long deprivations. Recovery is found to depend on circadian phase, and the model predicts that it is most efficient to recover during normal sleeping phases of the circadian cycle, in terms of the amount of recovery sleep required. The form of the homeostatic drive suggests that periods of wake during recovery from sleep deprivation are phases of relative recovery, in the sense that the homeostatic drive continues to converge toward baseline levels. This undermines the concept of sleep debt, and is in agreement with experimentally restricted recovery protocols. Finally, we compare our model to the two-process model, and demonstrate the power of physiologically based modeling by correctly predicting sleep latency times following deprivation from experimental data.

Mean field model of acetylcholine mediated dynamics in the thalamocortical system

A recent continuum model of the large scale electrical activity of the thalamocortical system is generalized to include cholinergic modulation. The model is examined analytically and numerically to determine the effect of acetylcholine (ACh) on its steady states, linear stability, spectrum, and temporal responses. Changing the ACh concentration moves the system between zones of one, three, and five steady states, showing that neuromodulation of synaptic strength is a possible mechanism by which multiple steady states emerge in the brain. The lowest firing rate steady state is always stable, and subsequent fixed points alternate between stable and unstable. Increasing ACh concentration changes the form of the spectrum. Increasing the tonic level of ACh concentration increases the magnitudes of the N100 and P200 in the evoked response potential (ERP), without changing the timing of these peaks. Driving the system with a pulse of cholinergic activity results in a transient increase in the firing rate of cortical neurons that lasts over 10s. Step-like increases in cortical ACh concentration cause increases in the firing rate of cortical neurons, with rapid responses due to fast acting nicotinic receptors and slower responses due to muscarinic receptor suppression of intracortical connections.

Eigenmode structure in solar-wind Langmuir waves

We show that observed spatial- and frequency-domain signatures of intense solar-wind Langmuir waves can be described as eigenmodes trapped in a parabolic density well. Measured solar-wind electric field spectra and waveforms are compared with 1D linear solutions and, in many cases, can be represented by 1-3 low-order eigenstates. To our knowledge, this report is the first observational confirmation of Langmuir eigenmodes in space. These results suggest that linear eigenmodes may be the starting point of the nonlinear evolution, critical for producing solar type II and type III radio bursts.

Controlling limit-cycle behaviors of brain activity

The limit cycles of brain activity are studied using a compact continuum model that reproduces the main features of electroencephalographic signals, including bifurcations of fixed points and limit cycles in seizures. Frequencies and amplitudes are predicted analytically and related to physiology. Gaussian stimuli yield two distinct evoked responses in the linearly stable zone, consistent with experiment. Limit cycles can be initiated or suppressed by control signals or stimuli.

Physiology-based modeling of cortical auditory evoked potentials

Evoked potentials are the transient electrical responses caused by changes in the brain following stimuli. This work uses a physiology-based continuum model of neuronal activity in the human brain to calculate theoretical cortical auditory evoked potentials (CAEPs) from the model's linearized response. These are fitted to experimental data, allowing the fitted parameters to be related to brain physiology. This approach yields excellent fits to CAEP data, which can then be compared to fits of EEG spectra. It is shown that the differences between resting eyes-open EEG and standard CAEPs can be explained by changes in the physiology of populations of neurons in corticothalamic pathways, with notable similarities to certain aspects of slow-wave sleep. This pilot study demonstrates the ability of our model-based fitting method to provide information on the underlying physiology of the brain that is not available using standard methods.

Physiologically based calculation of steady-state evoked potentials and cortical wave velocities

Steady-state evoked potentials (SSEPs) elicited by sinusoidal stimuli are predicted from a physiologically-based model, including bielectrode and volume conduction effects. Comparison with visual SSEPs yields constraints on phase and latency of the retinothalamic transfer function that are consistent with experiment. Predictions of phase velocities measured as SSEPs cross the cortex are consistent with low values measured for slow waves in sleep, while resonant behavior induced by corticothalamic loops, especially near the alpha peak, contributes to wide scatter in waking-state phase velocity measurements comparable to effects from volume conduction. The common use of bielectrode derivations to compensate for volume conduction effects is examined and shown to be incomplete, tending to lead to underestimates of phase velocity, especially at low frequencies and near the alpha peak, due to incorrect elimination of true long-wavelength contributions to the SSEP.

Mean field model of acetylcholine mediated dynamics in the cerebral cortex

A recent continuum model of the large scale electrical activity of the cerebral cortex is generalized to include cholinergic modulation. In this model, dynamic modulation of synaptic strength acts over the time scales of nicotinic and muscarinic receptor action. The cortical model is analyzed to determine the effect of acetylcholine (ACh) on its steady states, linear stability, spectrum, and temporal responses to changes in subcortical input. ACh increases the firing rate in steady states of the system. Changing ACh concentration does not introduce oscillatory behavior into the system, but increases the overall spectral power. Model responses to pulses in subcortical input are affected by the tonic level of ACh concentration, with higher levels of ACh increasing the magnitude firing rate response of excitatory cortical neurons to pulses of subcortical input. Numerical simulations are used to explore the temporal dynamics of the model in response to changes in ACh concentration. Evidence is seen of a transition from a state in which intracortical inputs are emphasized to a state where thalamic afferents have enhanced influence. Perturbations in ACh concentration cause changes in the firing rate of cortical neurons, with rapid responses due to fast acting facilitatory effects of nicotinic receptors on subcortical afferents, and slower responses due to muscarinic suppression of intracortical connections. Together, these numerical simulations demonstrate that the actions of ACh could be a significant factor modulating early components of evoked response potentials.

Laboratory evidence for stochastic plasma-wave growth

The first laboratory confirmation of stochastic growth theory is reported. Floating potential fluctuations are measured in a vacuum arc centrifuge using a Langmuir probe. Statistical analysis of the energy density reveals a lognormal distribution over roughly 2 orders of magnitude, with a high-field nonlinear cutoff whose spatial dependence is consistent with the predicted eigenmode profile. These results are consistent with stochastic growth and nonlinear saturation of a spatially extended eigenmode, the first evidence for stochastic growth of an extended structure.

Visual gamma oscillations: waves, correlations, and other phenomena, including comparison with experimental data

Mean-field theory of brain dynamics is applied to explain the properties of gamma (> or approximately 30 Hz) oscillations of cortical activity often seen during vision experiments. It is shown that mm-scale patchy connections in the primary visual cortex can support collective gamma oscillations with the correct frequencies and spatial structure, even when driven by uncorrelated inputs. This occurs via resonances associated with the the periodic modulation of the network connections, rather than being due to single-cell properties alone. Near-resonant gamma waves are shown to obey the Schrödinger equation, which enables techniques and insights from quantum theory to be used in exploring these classical oscillations. Resulting predictions for gamma responses to stimuli account in a unified way for a wide range of experimental results, including why oscillations and zero-lag synchrony are associated, and variations in correlation functions with time delay, intercellular distance, and stimulus features. They also imply that gamma oscillations may enable a form of frequency multiplexing of neural signals. Most importantly, it is shown that correlations reproduce experimental results that show maximal correlations between cells that respond to related features, but little correlation with other cells, an effect that has been argued to be associated with segmentation of a scene into separate objects. Consistency with infill of missing contours and increase in response with length of bar-shaped stimuli are discussed. Background correlations expected in the absence of stimulation are also calculated and shown to be consistent in form with experimental measurements and similar to stimulus-induced correlations in structure. Finally, possible links of gamma instabilities to certain classes of photically induced seizures and visual hallucinations are discussed.

Energy approach to rivalry dynamics, soliton stability, and pattern formation in neuronal networks

Hopfield's Lyapunov function is used to view the stability and topology of equilibria in neuronal networks for visual rivalry and pattern formation. For two neural populations with reciprocal inhibition and slow adaptation, the dynamics of neural activity is found to include a pair of limit cycles: one for oscillations between states where one population has high activity and the other has low activity, as in rivalry, and one for oscillations between states where both populations have the same activity. Hopfield's Lyapunov function is used to find the dynamical mechanism for oscillations and the basin of attraction of each limit cycle. For a spatially continuous population with lateral inhibition, stable equilibria are found for local regions of high activity (solitons) and for bound states of two or more solitons. Bound states become stable when moving two solitons together minimizes the Lyapunov function, a result of decreasing activity in regions between peaks of high activity when the firing rate is described by a sigmoid function. Lowering the barrier to soliton formation leads to a pattern-forming instability, and a nonlinear solution to the dynamical equations is found to be given by a soliton lattice, which is completely characterized by the soliton width and the spacing between neighboring solitons. Fluctuations due to noise create lattice vacancies analogous to point defects in crystals, leading to activity which is spatially inhomogeneous.

Spike-rate adaptation and neuronal bursting in a mean-field model of brain activity

Spike-rate adaptation is investigated within a mean-field model of brain activity. Two different mechanisms of negative feedback are considered; one involving modulation of the mean firing threshold, and the other, modulation of the mean synaptic strength. Adaptation to a constant stimulus is shown to take place for both mechanisms, and limit-cycle oscillations in the firing rate corresponding to bursts of neuronal activity are investigated. These oscillations are found to result from a Hopf bifurcation when the equilibrium lies between the local maximum and local minimum of a given nullcline. Oscillations with amplitudes significantly below the maximum firing rate are found over a narrow range of possible equilibriums.

Compact dynamical model of brain activity

A compact physiologically based mean-field formulation of brain dynamics is proposed to model observed brain activity and electroencephalographic (EEG) signals. In contrast to existing formulations, which are more detailed and complicated, our model is described by a single second-order delay differential equation that encapsulates salient aspects of the physiology. The model captures essential features of activity mediated by fast corticocortical connections and delayed feedbacks via extracortical pathways and external stimuli. In the linear regime, these features can be simply expressed by three coefficients derived from the properties of these physiological pathways and explicit nonlinear approximations are also derived. This compact model successfully reproduces the main features of experimental EEG's and the predictions of previous models, including resonance peaks in EEG spectra and nonlinear dynamics. As an illustration, key features of the dynamics of epileptic seizures are shown to be reproduced by the model. Due to its compact form, the model will facilitate insight into nonlinear brain dynamics via standard nonlinear techniques and will guide analysis and investigation of more complex models. It is thus a useful tool for analyzing complex brain activity, especially when it exhibits low-dimensional dynamics.

Transformation of arbitrary distributions to the normal distribution with application to EEG test-retest reliability

Many variables in the social, physical, and biosciences, including neuroscience, are non-normally distributed. To improve the statistical properties of such data, or to allow parametric testing, logarithmic or logit transformations are often used. Box-Cox transformations or ad hoc methods are sometimes used for parameters for which no transformation is known to approximate normality. However, these methods do not always give good agreement with the Gaussian. A transformation is discussed that maps probability distributions as closely as possible to the normal distribution, with exact agreement for continuous distributions. To illustrate, the transformation is applied to a theoretical distribution, and to quantitative electroencephalographic (qEEG) measures from repeat recordings of 32 subjects which are highly non-normal. Agreement with the Gaussian was better than using logarithmic, logit, or Box-Cox transformations. Since normal data have previously been shown to have better test-retest reliability than non-normal data under fairly general circumstances, the implications of our transformation for the test-retest reliability of parameters were investigated. Reliability was shown to improve with the transformation, where the improvement was comparable to that using Box-Cox. An advantage of the general transformation is that it does not require laborious optimization over a range of parameters or a case-specific choice of form.

The aggravating role of the ubiquitin–proteasome system in neurodegenerative disease

Intraneuronal inclusion bodies are key pathological features of most age-related neurodegenerative disorders including Parkinson's disease and Alzheimer's disease. These inclusions are commonly characterized both by the presence of ubiquitinated proteins and the sequestration of components of the UPS (ubiquitin–proteasome system). Unfortunately, as we age, the efficiency of the UPS declines, suggesting that the presence of ubiquitinated proteins and UPS components in inclusions may reflect unsuccessful attempts by the (failing) UPS to remove the aggregating proteins. Whether the physical presence of inclusions causes cell death or, conversely, whether they are non-toxic and their presence reflects a cellular protective mechanism remains highly controversial. Animal and in vitro model systems that allow detailed characterization of the inclusions and their effects on the cell have been developed by us and others. Identification of the mechanisms involved in inclusion formation is already aiding the development of novel therapeutic strategies to prevent or alleviate aggregate-associated neurodegenerative diseases.

Hypotensive effects of hawthorn for patients with diabetes taking prescription drugs: a randomised controlled trial

BACKGROUND

Hawthorn (Crataegus laevigata) leaves, flowers and berries are used by herbal practitioners in the UK to treat hypertension in conjunction with prescribed drugs. Small-scale human studies support this approach.

AIM

To investigate the effects of hawthorn for hypertension in patients with type 2 diabetes taking prescribed drugs.

DESIGN OF STUDY

Randomised controlled trial.

SETTING

General practices in Reading, UK.

METHOD

Patients with type 2 diabetes (n = 79) were randomised to daily 1200 mg hawthorn extract (n = 39) or placebo (n = 40) for 16 weeks. At baseline and outcome a wellbeing questionnaire was completed and blood pressure and fasting blood samples taken. A food frequency questionnaire estimated nutrient intake.

RESULTS

Hypotensive drugs were used by 71% of the study population with a mean intake of 4.4 hypoglycaemic and/or hypotensive drugs. Fat intake was lower and sugar intake higher than recommendations, and low micronutrient intake was prevalent. There was a significant group difference in mean diastolic blood pressure reductions (P = 0.035): the hawthorn group showed greater reductions (baseline: 85.6 mmHg, 95% confidence interval [CI] = 83.3 to 87.8; outcome: 83.0 mmHg, 95% CI = 80.5 to 85.7) than the placebo group (baseline: 84.5 mmHg, 95% CI = 82 to 87; outcome: 85.0 mmHg, 95% CI = 82.2 to 87.8). There was no group difference in systolic blood pressure reduction from baseline (3.6 and 0.8 mmHg for hawthorn and placebo groups, respectively; P = 0.329). Although mean fat intake met current recommendations, mean sugar intake was higher and there were indications of potential multiple micronutrient deficiencies. No herb-drug interaction was found and minor health complaints were reduced from baseline in both groups.

CONCLUSIONS

This is the first randomised controlled trial to demonstrate a hypotensive effect of hawthorn in patients with diabetes taking medication.

Multielectrode electroencephalogram power spectra: theory and application to approximate correction of volume conduction effects

Using a physiologically based model of brain activity, electroencephalogram (EEG) power spectra are calculated for signals derived from general linear combinations of voltages from multiple electrodes, with and without filtering by volume conduction. Two simple methods of combining scalp measurements to estimate unfiltered EEG power spectra are then proposed and their accuracy and robustness are explored, using the model predictions as an illustration. It is found that these methods, including a case that uses just three electrodes, enable improved estimation of the underlying spectrum relative to each of several widely used combinations alone.