Model estimation of cerebral hemodynamics between blood flow and volume changes: A data-based modeling approach

It is well known that there is a dynamic relationship between cerebral blood flow (CBF) and cerebral blood volume (CBV). With increasing applications of functional MRI, where the blood oxygen-level-dependent signals are recorded, the understanding and accurate modeling of the hemodynamic relationship between CBF and CBV becomes increasingly important. This study presents an empirical and data-based modeling framework for model identification from CBF and CBV experimental data. It is shown that the relationship between the changes in CBF and CBV can be described using a parsimonious autoregressive with exogenous input model structure. It is observed that neither the ordinary least-squares (LS) method nor the classical total least-squares (TLS) method can produce accurate estimates from the original noisy CBF and CBV data. A regularized total least-squares (RTLS) method is thus introduced and extended to solve such an error-in-the-variables problem. Quantitative results show that the RTLS method works very well on the noisy CBF and CBV data. Finally, a combination of RTLS with a filtering method can lead to a parsimonious but very effective model that can characterize the relationship between the changes in CBF and CBV.

The hemodynamic impulse response to a single neural event

This article investigates the relation between stimulus-evoked neural activity and cerebral hemodynamics. Specifically, the hypothesis is tested that hemodynamic responses can be modeled as a linear convolution of experimentally obtained measures of neural activity with a suitable hemodynamic impulse response function. To obtain a range of neural and hemodynamic responses, rat whisker pad was stimulated using brief (less than or equal to2 seconds) electrical stimuli consisting of single pulses (0.3 millisecond, 1.2 mA) combined both at different frequencies and in a paired-pulse design. Hemodynamic responses were measured using concurrent optical imaging spectroscopy and laser Doppler flowmetry, whereas neural responses were assessed through current source density analysis of multielectrode recordings from a single barrel. General linear modeling was used to deconvolve the hemodynamic impulse response to a single "neural event" from the hemodynamic and neural responses to stimulation. The model provided an excellent fit to the empirical data. The implications of these results for modeling schemes and for physiologic systems coupling neural and hemodynamic activity are discussed.

Therapeutic potential of computer to cerebral cortex implantable devices

In this article, an overview of some of the latest developments in the field of cerebral cortex to computer interfacing (CCCI) is given. This is posed in the more general context of Brain-Computer Interfaces in order to assess advantages and disadvantages. The emphasis is clearly placed on practical studies that have been undertaken and reported on, as opposed to those speculated, simulated or proposed as future projects. Related areas are discussed briefly only in the context of their contribution to the studies being undertaken. The area of focus is notably the use of invasive implant technology, where a connection is made directly with the cerebral cortex and/or nervous system. Tests and experimentation which do not involve human subjects are invariably carried out a priori to indicate the eventual possibilities before human subjects are themselves involved. Some of the more pertinent animal studies from this area are discussed. The paper goes on to describe human experimentation, in which neural implants have linked the human nervous system bidirectionally with technology and the internet. A view is taken as to the prospects for the future for CCCI, in terms of its broad therapeutic role.

Increased protein SUMOylation following focal cerebral ischemia

Stroke is a major cause of death and disability, which involves excessive glutamate receptor activation leading to excitotoxic cell death. We recently reported that SUMOylation can regulate kainate receptor (KAR) function. Here we investigated changes in protein SUMOylation and levels of KAR and AMPA receptor subunits in two different animal stroke models: a rat model of focal ischemia with reperfusion and a mouse model without reperfusion. In rats, transient middle cerebral artery occlusion (MCAO) resulted in a striatal and cortical infarct. A dramatic increase in SUMOylation by both SUMO-1 and SUMO-2/3 was observed at 6h and 24h in the striatal infarct area and by SUMO-2/3 at 24h in the hippocampus, which was not directly subjected to ischemia. In mice, permanent MCAO resulted in a selective cortical infarct. No changes in SUMOylation occurred at 6h but there was increased SUMO-1 conjugation in the cortical infarct and non-ischemic hippocampus at 24h after MCAO. Interestingly, SUMOylation by SUMO-2/3 occurred only outside the infarct area. In both rat and mouse levels of KARs were only decreased in the infarct regions whereas AMPARs were decreased in the infarct and in other brain areas. These results suggest that posttranslational modification by SUMO and down-regulation of AMPARs and KARs may play important roles in the pathophysiological response to ischemia.

Nitric oxide suppresses cerebral vasomotion by sGC-independent effects on ryanodine receptors and voltage-gated calcium channels

Background/Aims: In cerebral arteries, nitric oxide (NO) release plays a key role in suppressing vasomotion. Our aim was to establish the pathways affected by NO in rat middle cerebral arteries. Methods: In isolated segments of artery, isometric tension and simultaneous measurements of either smooth muscle membrane potential or intracellular [Ca 2+ ] ([Ca 2+ ] SMC ) changes were recorded. Results: In the absence of L -NAME, asynchronous propagating Ca 2+ waves were recorded that were sensitive to block with ryanodine, but not nifedipine. L -NAME stimulated pronounced vasomotion and synchronous Ca 2+ oscillations with close temporal coupling between membrane potential, tone and [Ca 2+ ] SMC . If nifedipine was applied together with L -NAME, [Ca 2+ ] SMC decreased and synchronous Ca 2+ oscillations were lost, but asynchronous propagating Ca 2+ waves persisted. Vasomotion was similarly evoked by either iberiotoxin, or by ryanodine, and to a lesser extent by ODQ. Exogenous application of NONOate stimulated endothelium-independent hyperpolarization and relaxation of either L -NAME-induced or spontaneous arterial tone. NO-evoked hyperpolarization involved activation of BK Ca channels via ryanodine receptors (RYRs), with little involvement of sGC. Further, in whole cell mode, NO inhibited current through L-type voltage-gated Ca 2+ channels (VGCC), which was independent of both voltage and sGC. Conclusion: NO exerts sGC-independent actions at RYRs and at VGCC, both of which normally suppress cerebral artery myogenic tone.

Evidence both L-type and non-L-type voltage-dependent calcium channels contribute to cerebral artery vasospasm following loss of NO in the rat

We recently found block of NO synthase in rat middle cerebral artery caused spasm, associated with depolarizing oscillations in membrane potential (Em) similar in form but faster in frequency (circa 1 Hz) to vasomotion. T-type voltage-gated Ca2+ channels contribute to cerebral myogenic tone and vasomotion, so we investigated the significance of T-type and other ion channels for membrane potential oscillations underlying arterial spasm. Smooth muscle cell membrane potential (Em) and tension were measured simultaneously in rat middle cerebral artery. NO synthase blockade caused temporally coupled depolarizing oscillations in cerebrovascular Em with associated vasoconstriction. Both events were accentuated by block of smooth muscle BKCa. Block of T-type channels or inhibition of Na+/K+-ATPase abolished the oscillations in Em and reduced vasoconstriction. Oscillations in Em were either attenuated or accentuated by reducing [Ca2+]o or block of KV, respectively. TRAM-34 attenuated oscillations in both Em and tone, apparently independent of effects against KCa3.1. Thus, rapid depolarizing oscillations in Em and tone observed after endothelial function has been disrupted reflect input from T-type calcium channels in addition to L-type channels, while other depolarizing currents appear to be unimportant. These data suggest that combined block of T and L-type channels may represent an effective approach to reverse cerebral vasospasm.

Thromboxane A2 inhibition of SKCa after NO synthase block in rat middle cerebral artery

NO/prostanoid independent, EDHF-mediated hyperpolarization and dilation in rat middle cerebral arteries is mediated solely by endothelial cell IK(Ca). However, when the NO-pathway is also active, both SK(Ca) and IK(Ca) contribute to EDHF responses. As the SK(Ca) component can be inhibited by stimulation of thromboxane A(2) (TxA(2)) TP receptors and NO has the potential ability to inhibit thromboxane synthesis, we investigated whether TxA(2) might explain loss of functional input from SK(Ca) during NOS inhibition in cerebral arteries. EXPERIMENTAL APPROACH: Rat middle cerebral arteries were mounted in a wire myograph. Endothelium-dependent responses to the PAR2 agonist, SLIGRL were assessed as simultaneous changes in smooth muscle membrane potential and tension. KEY RESULTS: Responses were obtained in the presence of L-NAME as appropriate. Inhibition of TP receptors with either ICI 192,605 or SQ 29,548, did not affect EDHF mediated hyperpolarization and relaxation, but in their presence neither TRAM-34 nor apamin (to block IK(Ca) and SK(Ca) respectively) individually affected the EDHF response. However, in combination they virtually abolished it. Similar effects were obtained in the presence of the thromboxane synthase inhibitor, furegrelate, which additionally revealed an iberiotoxin-sensitive residual EDHF hyperpolarization and relaxation in the combined presence of TRAM-34 and apamin. CONCLUSIONS AND IMPLICATIONS: In the rat middle cerebral artery, inhibition of NOS leads to a loss of the SK(Ca) component of EDHF responses. Either antagonism of TP receptors or block of thromboxane synthase restores an input through SK(Ca). These data indicate that NO normally enables SK(Ca) activity in rat middle cerebral arteries.

Evidence for involvement of both IKCa and SKCa channels in hyperpolarizing responses of the rat middle cerebral artery

Endothelium-derived hyperpolarizing factor responses in the rat middle cerebral artery are blocked by inhibiting IKCa channels alone, contrasting with peripheral vessels where block of both IKCa and SKCa is required. As the contribution of IKCa and SKCa to endothelium-dependent hyperpolarization differs in peripheral arteries, depending on the level of arterial constriction, we investigated the possibility that SKCa might contribute to equivalent hyperpolarization in cerebral arteries under certain conditions. METHODS: Rat middle cerebral arteries (approximately 175 microm) were mounted in a wire myograph. The effect of KCa channel blockers on endothelium-dependent responses to the protease-activated receptor 2 agonist, SLIGRL (20 micromol/L), were then assessed as simultaneous changes in tension and membrane potential. These data were correlated with the distribution of arterial KCa channels revealed with immunohistochemistry. RESULTS: SLIGRL hyperpolarized and relaxed cerebral arteries undergoing variable levels of stretch-induced tone. The relaxation was unaffected by specific inhibitors of IKCa (TRAM-34, 1 micromol/L) or SKCa (apamin, 50 nmol/L) alone or in combination. In contrast, the associated smooth-muscle hyperpolarization was inhibited, but only with these blockers in combination. Blocking nitric oxide synthase (NOS) or guanylyl cyclase evoked smooth-muscle depolarization and constriction, with both hyperpolarization and relaxation to SLIGRL being abolished by TRAM-34 alone, whereas apamin had no effect. Immunolabeling showed SKCa and IKCa within the endothelium. CONCLUSIONS: In the absence of NO, IKCa underpins endothelium-dependent hyperpolarization and relaxation in cerebral arteries. However, when NOS is active SKCa contributes to hyperpolarization, whatever the extent of background contraction. These changes may have relevance in vascular disease states where NO release is compromised and when the levels of SKCa expression may be altered.

Possible role for K+ in endothelium-derived hyperpolarizing factor–linked dilatation in rat middle cerebral artery

Background and Purpose— Endothelium-derived hyperpolarizing factor (EDHF) and K+ are vasodilators in the cerebral circulation. Recently, K+ has been suggested to contribute to EDHF-mediated responses in peripheral vessels. The EDHF response to the protease-activated receptor 2 ligand SLIGRL was characterized in cerebral arteries and used to assess whether K+ contributes as an EDHF. Methods— Rat middle cerebral arteries were mounted in either a wire or pressure myograph. Concentration-response curves to SLIGRL and K+ were constructed in the presence and absence of a variety of blocking agents. In some experiments, changes in tension and smooth muscle cell membrane potential were recorded simultaneously. Results— SLIGRL (0.02 to 20 μmol/L) stimulated concentration and endothelium-dependent relaxation. In the presence of NG-nitro-L-arginine methyl ester, relaxation to SLIGRL was associated with hyperpolarization and sensitivity to a specific inhibitor of IKCa, 1-[(2-chlorophenyl)diphenylmethyl]-1H-pyrazole (1μmol/L), reflecting activation of EDHF. Combined inhibition of KIR with Ba2+ (30μmol/L) and Na+/K+-ATPase with ouabain (1 μmol/L) markedly attenuated the relaxation to EDHF. Raising extracellular [K+] to 15 mmol/L also stimulated smooth muscle relaxation and hyperpolarization, which was also attenuated by combined application of Ba2+ and ouabain. Conclusions— SLIGRL evokes EDHF-mediated relaxation in the rat middle cerebral artery, underpinned by hyperpolarization of the smooth muscle. The profile of blockade of EDHF-mediated hyperpolarization and relaxation supports a pivotal role for IKCa channels. Furthermore, similar inhibition of responses to EDHF and exogenous K+ with Ba2+ and ouabain suggests that K+ may contribute as an EDHF in the middle cerebral artery.

Statins and selective inhibition of Rho kinase protect small conductance calcium-activated potassium channel function (KCa2.3) in cerebral arteries

Background: In rat middle cerebral and mesenteric arteries the KCa2.3 component of endothelium-dependent hyperpolarization (EDH) is lost following stimulation of thromboxane (TP) receptors, an effect that may contribute to the endothelial dysfunction associated with cardiovascular disease. In cerebral arteries, KCa2.3 loss is associated with NO synthase inhibition, but is restored if TP receptors are blocked. The Rho/Rho kinase pathway is central for TP signalling and statins indirectly inhibit this pathway. The possibility that Rho kinase inhibition and statins sustain KCa2.3 hyperpolarization was investigated in rat middle cerebral arteries (MCA). Methods: MCAs were mounted in a wire myograph. The PAR2 agonist, SLIGRL was used to stimulate EDH responses, assessed by simultaneous measurement of smooth muscle membrane potential and tension. TP expression was assessed with rt-PCR and immunofluorescence. Results: Immunofluorescence detected TP in the endothelial cell layer of MCA. Vasoconstriction to the TP agonist, U46619 was reduced by Rho kinase inhibition. TP receptor stimulation lead to loss of KCa2.3 mediated hyperpolarization, an effect that was reversed by Rho kinase inhibitors or simvastatin. KCa2.3 activity was lost in L-NAME-treated arteries, but was restored by Rho kinase inhibition or statin treatment. The restorative effect of simvastatin was blocked after incubation with geranylgeranyl-pyrophosphate to circumvent loss of isoprenylation. Conclusions: Rho/Rho kinase signalling following TP stimulation and L-NAME regulates endothelial cell KCa2.3 function. The ability of statins to prevent isoprenylation and perhaps inhibit of Rho restores/protects the input of KCa2.3 to EDH in the MCA, and represents a beneficial pleiotropic effect of statin treatment.

Early and late stimulus-evoked cortical hemodynamic responses provide insight into the neurogenic nature of neurovascular coupling

Understanding neurovascular coupling is a prerequisite for the interpretation of results obtained from modern neuroimaging techniques. This study investigated the hemodynamic and neural responses in rat somatosensory cortex elicited by 16 seconds electrical whisker stimuli. Hemodynamics were measured by optical imaging spectroscopy and neural activity by multichannel electrophysiology. Previous studies have suggested that the whisker-evoked hemodynamic response contains two mechanisms, a transient ‘backwards’ dilation of the middle cerebral artery, followed by an increase in blood volume localized to the site of neural activity. To distinguish between the mechanisms responsible for these aspects of the response, we presented whisker stimuli during normocapnia (‘control’), and during a high level of hypercapnia. Hypercapnia was used to ‘predilate’ arteries and thus possibly ‘inhibit’ aspects of the response related to the ‘early’ mechanism. Indeed, hemodynamic data suggested that the transient stimulus-evoked response was absent under hypercapnia. However, evoked neural responses were also altered during hypercapnia and convolution of the neural responses from both the normocapnic and hypercapnic conditions with a canonical impulse response function, suggested that neurovascular coupling was similar in both conditions. Although data did not clearly dissociate early and late vascular responses, they suggest that the neurovascular coupling relationship is neurogenic in origin.

A three-compartment model of the hemodynamic response and oxygen delivery to brain

We describe a mathematical model linking changes in cerebral blood flow, blood volume and the blood oxygenation state in response to stimulation. The model has three compartments to take into account the fact that the cerebral blood flow and volume as measured concurrently using laser Doppler flowmetry and optical imaging spectroscopy have contributions from the arterial, capillary as well as the venous compartments of the vasculature. It is an extension to previous one-compartment hemodynamic models which assume that the measured blood volume changes are from the venous compartment only. An important assumption of the model is that the tissue oxygen concentration is a time varying state variable of the system and is driven by the changes in metabolic demand resulting from changes in neural activity. The model takes into account the pre-capillary oxygen diffusion by flexibly allowing the saturation of the arterial compartment to be less than unity. Simulations are used to explore the sensitivity of the model and to optimise the parameters for experimental data. We conclude that the three-compartment model was better than the one-compartment model at capturing the hemodynamics of the response to changes in neural activation following stimulation.

Reduced-order modeling of light transport in tissue for real-time monitoring of brain hemodynamics using diffuse optical tomography

This paper proposes a new reconstruction method for diffuse optical tomography using reduced-order models of light transport in tissue. The models, which directly map optical tissue parameters to optical flux measurements at the detector locations, are derived based on data generated by numerical simulation of a reference model. The reconstruction algorithm based on the reduced-order models is a few orders of magnitude faster than the one based on a finite element approximation on a fine mesh incorporating a priori anatomical information acquired by magnetic resonance imaging. We demonstrate the accuracy and speed of the approach using a phantom experiment and through numerical simulation of brain activation in a rat's head. The applicability of the approach for real-time monitoring of brain hemodynamics is demonstrated through a hypercapnic experiment. We show that our results agree with the expected physiological changes and with results of a similar experimental study. However, by using our approach, a three-dimensional tomographic reconstruction can be performed in ∼3  s per time point instead of the 1 to 2 h it takes when using the conventional finite element modeling approach

Regulation of KCa2.3 andendothelium-dependent hyperpolarization (EDH) in the rat middle cerebral artery: the role of lipoxygenase metabolites and isoprostanes

Background and Purpose. In rat middle cerebral arteries, endothelium-dependent hyperpolarization (EDH) is mediated by activation of calcium-activated potassium(KCa) channels specifically KCa2.3 and KCa3.1. Lipoxygenase (LOX) products function as endothelium-derived hyperpolarizing factors (EDHFs) in rabbit arteries by stimulating KCa2.3. We investigated if LOX products contribute to EDH in rat cerebral arteries. Methods. Arachidonic acid (AA) metabolites produced in middle cerebral arteries were measured using HPLC and LC/MS. Vascular tension and membrane potential responses to SLIGRL were simultaneously recorded using wire myography and intracellular microelectrodes. Results. SLIGRL, an agonist at PAR2 receptors, caused EDH that was inhibited by a combination of KCa2.3 and KCa3.1 blockade. Non-selective LOX-inhibition reduced EDH, whereas inhibition of 12-LOX had no effect. Soluble epoxide hydrolase (sEH) inhibition enhanced the KCa2.3 component of EDH. Following NO synthase (NOS) inhibition, the KCa2.3 component of EDH was absent. Using HPLC, middle cerebral arteries metabolized 14C-AA to 15- and 12-LOX products under control conditions. With NOS inhibition, there was little change in LOX metabolites, but increased F-type isoprostanes. 8-iso-PGF2α inhibited the KCa2.3 component of EDH. Conclusions. LOX metabolites mediate EDH in rat middle cerebral arteries. Inhibition of sEH increases the KCa2.3 component of EDH. Following NOS inhibition,loss of KCa2.3 function is independent of changes in LOX production or sEH inhibition but due to increased isoprostane production and subsequent stimulation of TP receptors. These findings have important implications in diseases associated with loss of NO signaling such as stroke; where inhibition of sEH and/or isoprostane formation may of benefit.

Exploration of the neural correlates of cerebral palsy for sensorimotor BCI control

Cerebral palsy (CP) includes a broad range of disorders, which can result in impairment of posture and movement control. Brain-computer interfaces (BCIs) have been proposed as assistive devices for individuals with CP. Better understanding of the neural processing underlying motor control in affected individuals could lead to more targeted BCI rehabilitation and treatment options. We have explored well-known neural correlates of movement, including event-related desynchronization (ERD), phase synchrony, and a recently-introduced measure of phase dynamics, in participants with CP and healthy control participants. Although present, significantly less ERD and phase locking were found in the group with CP. Additionally, inter-group differences in phase dynamics were also significant. Taken together these findings suggest that users with CP exhibit lower levels of motor cortex activation during motor imagery, as reflected in lower levels of ongoing mu suppression and less functional connectivity. These differences indicate that development of BCIs for individuals with CP may pose additional challenges beyond those faced in providing BCIs to healthy individuals.