Mechanistic Mathematical Modeling Tests Hypotheses of the Neurovascular Coupling in fMRI.

Functional magnetic resonance imaging (fMRI) measures brain activity by detecting the blood-oxygen-level dependent (BOLD) response to neural activity. The BOLD response depends on the neurovascular coupling, which connects cerebral blood flow, cerebral blood volume, and deoxyhemoglobin level to neuronal activity. The exact mechanisms behind this neurovascular coupling are not yet fully investigated. There are at least three different ways in which these mechanisms are being discussed. Firstly, mathematical models involving the so-called Balloon model describes the relation between oxygen metabolism, cerebral blood volume, and cerebral blood flow. However, the Balloon model does not describe cellular and biochemical mechanisms. Secondly, the metabolic feedback hypothesis, which is based on experimental findings on metabolism associated with brain activation, and thirdly, the neurotransmitter feed-forward hypothesis which describes intracellular pathways leading to vasoactive substance release. Both the metabolic feedback and the neurotransmitter feed-forward hypotheses have been extensively studied, but only experimentally. These two hypotheses have never been implemented as mathematical models. Here we investigate these two hypotheses by mechanistic mathematical modeling using a systems biology approach; these methods have been used in biological research for many years but never been applied to the BOLD response in fMRI. In the current work, model structures describing the metabolic feedback and the neurotransmitter feed-forward hypotheses were applied to measured BOLD responses in the visual cortex of 12 healthy volunteers. Evaluating each hypothesis separately shows that neither hypothesis alone can describe the data in a biologically plausible way. However, by adding metabolism to the neurotransmitter feed-forward model structure, we obtained a new model structure which is able to fit the estimation data and successfully predict new, independent validation data. These results open the door to a new type of fMRI analysis that more accurately reflects the true neuronal activity.

Positive Allosteric Modulator of GABA Lowers BOLD Responses in the Cingulate Cortex.

Knowledge about the neural underpinnings of the negative blood oxygen level dependent (BOLD) responses in functional magnetic resonance imaging (fMRI) is still limited. We hypothesized that pharmacological GABAergic modulation attenuates BOLD responses, and that blood concentrations of a positive allosteric modulator of GABA correlate inversely with BOLD responses in the cingulate cortex. We investigated whether or not pure task-related negative BOLD responses were co-localized with pharmacologically modulated BOLD responses. Twenty healthy adults received either 5 mg diazepam or placebo in a double blind, randomized design. During fMRI the subjects performed a working memory task. Results showed that BOLD responses in the cingulate cortex were inversely correlated with diazepam blood concentrations; that is, the higher the blood diazepam concentration, the lower the BOLD response. This inverse correlation was most pronounced in the pregenual anterior cingulate cortex and the anterior mid-cingulate cortex. For subjects with diazepam plasma concentration > 0.1 mg/L we observed negative BOLD responses with respect to fixation baseline. There was minor overlap between cingulate regions with task-related negative BOLD responses and regions where the BOLD responses were inversely correlated with diazepam concentration. We interpret that the inverse correlation between the BOLD response and diazepam was caused by GABA-related neural inhibition. Thus, this study supports the hypothesis that GABA attenuates BOLD responses in fMRI. The minimal overlap between task-related negative BOLD responses and responses attenuated by diazepam suggests that these responses might be caused by different mechanisms.

Mathematical modeling improves EC50 estimations from classical dose-response curves.

UNLABELLED: The ß-adrenergic response is impaired in failing hearts. When studying ß-adrenergic function in vitro, the half-maximal effective concentration (EC50 ) is an important measure of ligand response. We previously measured the in vitro contraction force response of chicken heart tissue to increasing concentrations of adrenaline, and observed a decreasing response at high concentrations. The classical interpretation of such data is to assume a maximal response before the decrease, and to fit a sigmoid curve to the remaining data to determine EC50 . Instead, we have applied a mathematical modeling approach to interpret the full dose-response curve in a new way. The developed model predicts a non-steady-state caused by a short resting time between increased concentrations of agonist, which affect the dose-response characterization. Therefore, an improved estimate of EC50 may be calculated using steady-state simulations of the model. The model-based estimation of EC50 is further refined using additional time-resolved data to decrease the uncertainty of the prediction. The resulting model-based EC50 (180-525 nm) is higher than the classically interpreted EC50 (46-191 nm). Mathematical modeling thus makes it possible to re-interpret previously obtained datasets, and to make accurate estimates of EC50 even when steady-state measurements are not experimentally feasible. DATABASE: The mathematical models described here have been submitted to the JWS Online Cellular Systems Modelling Database, and may be accessed at http://jjj.bio.vu.nl/database/nyman.

Fetal development of baroreflex sensitivity: the chicken embryo as a case model.

The baroreflex is the main short term compensatory mechanism to buffer blood pressure changes and maintain circulatory homeostasis. Its ontogeny and importance during prenatal life is not fully understood so we used broiler chickens to investigate the maturation of the baroreflex in late incubation using a novel method that measured changes in heart rate during spontaneous fluctuations in blood pressure. Our results suggest that a baroreflex is already functional at d17 with no indication of further maturation in terms of sensitivity (gain at 17 d was 52.9±8.3 and at 20 d 69.5±16.2 ms kPa(-1)). The physiological relevance of these values is shown using data surrogation methods. Although the results contrast with the progressive baroreflex maturation indicated by the pharmacological method, we sustain that both methods provide information on baroreflex regulation. While the spontaneous method evaluates truly physiological (but small) pressure changes, the pharmacological method provides a more consistent and repetitive challenge for the reflex that requires a different recruitment of baroreflex effectors.