Lactate's behavioral switch in the brain: An in-silico model.

Emerging evidence emphasizes lactate's involvement in both physiological processes (energy metabolism, memory, etc.) and disease (traumatic brain injury, epilepsy, etc.). Furthermore, the usefulness of mathematical modeling in deciphering underlying dynamics of the brain to investigate lactate roles and mechanisms of action has been well established. Here, we analyze a novel mathematical model of brain lactate exchanges between four compartments: neurons, astrocytes, capillaries, and extracellular space. A system of four ordinary differential equations is proposed to explain interactions between these compartments. We first optimize and analyze the model's parameters under normal, resting state conditions, and then use it to simulate changes linked to elevated arterial lactate. Our results show that even though increased arterial lactate results in increased uptake by astrocytes and release to the extracellular space, it cannot strongly recover the initial drop in neuronal lactate concentration. Also, we show that the direction of lactate transport between the compartments is influenced by the maximum astrocyte production rate and the transport rate between astrocytes and extracellular space.

A New Computational Model for Neuro-Glio-Vascular Coupling: Astrocyte Activation Can Explain Cerebral Blood Flow Nonlinear Response to Interictal Events.

Developing a clear understanding of the relationship between cerebral blood flow (CBF) response and neuronal activity is of significant importance because CBF increase is essential to the health of neurons, for instance through oxygen supply. This relationship can be investigated by analyzing multimodal (fMRI, PET, laser Doppler…) recordings. However, the important number of intermediate (non-observable) variables involved in the underlying neurovascular coupling makes the discovery of mechanisms all the more difficult from the sole multimodal data. We present a new computational model developed at the population scale (voxel) with physiologically relevant but simple equations to facilitate the interpretation of regional multimodal recordings. This model links neuronal activity to regional CBF dynamics through neuro-glio-vascular coupling. This coupling involves a population of glial cells called astrocytes via their role in neurotransmitter (glutamate and GABA) recycling and their impact on neighboring vessels. In epilepsy, neuronal networks generate epileptiform discharges, leading to variations in astrocytic and CBF dynamics. In this study, we took advantage of these large variations in neuronal activity magnitude to test the capacity of our model to reproduce experimental data. We compared simulations from our model with isolated epileptiform events, which were obtained in vivo by simultaneous local field potential and laser Doppler recordings in rats after local bicuculline injection. We showed a predominant neuronal contribution for low level discharges and a significant astrocytic contribution for higher level discharges. Besides, neuronal contribution to CBF was linear while astrocytic contribution was nonlinear. Results thus indicate that the relationship between neuronal activity and CBF magnitudes can be nonlinear for isolated events and that this nonlinearity is due to astrocytic activity, highlighting the importance of astrocytes in the interpretation of regional recordings.

Modeling of the neurovascular coupling in epileptic discharges.

Despite the interest in simultaneous EEG-fMRI studies of epileptic spikes, the link between epileptic discharges and their corresponding hemodynamic responses is poorly understood. In this context, biophysical models are promising tools for investigating the mechanisms underlying observed signals. Here, we apply a metabolic-hemodynamic model to simulated epileptic discharges, in part generated by a neural mass model. We analyze the effect of features specific to epileptic neuronal activity on the blood oxygen level dependent (BOLD) response, focusing on the issues of linearity in neurovascular coupling and on the origin of negative BOLD signals. We found both sub- and supra-linearity in simulated BOLD signals, depending on whether one observes the early or the late part of the BOLD response. The size of these non-linear effects is determined by the spike frequency, as well as by the amplitude of the excitatory activity. Our results additionally indicate a minor deviation from linearity at the neuronal level. According to a phase space analysis, the possibility to obtain a negative BOLD response to an epileptic spike depends on the existence of a long and strong excitatory undershoot. Moreover, we strongly suggest that a combined EEG-fMRI modeling approach should include spatial assumptions. The present study is a step towards an increased understanding of the link between epileptic spikes and their BOLD responses, aiming to improve the interpretation of simultaneous EEG-fMRI recordings in epilepsy.

Relationship between flow and metabolism in BOLD signals: insights from biophysical models.

In many physiological or pathological situations, the interpretation of BOLD signals remains elusive as the intimate link between neuronal activity and subsequent flow/metabolic changes is not fully understood. During the past decades, a number of biophysical models of the neurovascular coupling have been proposed. It is now well-admitted that these models may bridge between observations (fMRI data) and underlying biophysical and (patho-)physiological mechanisms (related to flow and metabolism) by providing mechanistic explanations. In this study, three well-established models (Buxton's, Friston's and Sotero's) are investigated. An exhaustive parameter sensitivity analysis (PSA) was conducted to study the marginal and joint influences of model parameters on the three main features of the BOLD response (namely the principal peak, the post-stimulus undershoot and the initial dip). In each model, parameters that have the greatest (and least) influence on the BOLD features as well as on the direction of variation of these features were identified. Among the three studied models, parameters were shown to affect the output features in different manners. Indeed, the main parameters revealed by the PSA were found to strongly depend on the way the flow(CBF)-metabolism(CMRO(2)) relationship is implemented (serial vs. parallel). This study confirmed that the model structure which accounts for the representation of the CBF-CMRO(2) relationship (oxygen supply by the flow vs. oxygen demand from neurons) plays a key role. More generally, this work provides substantial information about the tuning of parameters in the three considered models and about the subsequent interpretation of BOLD signals based on these models.

Optimizing the speckle noise for maximum efficacy of data acquisition in coherent imaging.

The impact of multiplicative speckle noise on data acquisition in coherent imaging is studied. This demonstrates the possibility to optimally adjust the level of the speckle noise in order to deliberately exploit, with maximum efficacy, the saturation naturally limiting linear image sensors such as CCD cameras, for instance. This constructive action of speckle noise cooperating with saturation can be interpreted as a novel instance of stochastic resonance or a useful-noise effect.

Constructive action of the speckle noise in a coherent imaging system.

A coherent imaging system with speckle noise is devised and analyzed. This demonstrates the possibility of improving the nonlinear transmission of a coherent image by increasing the level of the multiplicative speckle noise. This noise-assisted image transmission is a novel instance of stochastic resonance phenomena by which nonlinear signal processing benefits from a constructive action of noise.

Noise-enhanced Fisher information in parallel arrays of sensors with saturation.

This paper investigates stochastic resonance in parallel arrays of uncoupled saturating devices. The Fisher information is used to demonstrate the possibility of noise improved parameter estimation for arbitrary parametric signals. Especially, it is shown that improvement by noise always occurs in these arrays, for any configuration of the input signal, even in optimal configuration. The results contribute to establish stochastic resonance in parallel uncoupled arrays as a general mechanism of enhancement by noise, which can occur in wide classes of nonlinearities and for various information processing tasks. It can supplement other mechanisms of stochastic resonance that take place in isolated nonlinearities but generally in restricted configurations of the input signal.