Effectively desynchronizing deep brain stimulation based on a coordinated delayed feedback stimulation via several sites: a computational study

In detailed simulations we present a coordinated delayed feedback stimulation as a particularly robust and mild technique for desynchronization. We feed back the measured and band-pass filtered local filed potential via several or multiple sites with different delays, respectively. This yields a resounding desynchronization in a naturally demand-controlled way. Our novel approach is superior to previously developed techniques: It is robust against variations of system parameters, e.g., the mean firing rate. It does not require time-consuming calibration. It also prevents intermittent resynchronization typically caused by all methods employing repetitive administration of shocks. We suggest our novel technique to be used for deep brain stimulation in patients suffering from neurological diseases with pathological synchronization, such as Parkinsonian tremor, essential tremor or epilepsy.     

Cortical activation changes during repeated laser stimulation: a magnetoencephalographic study

Repeated warm laser stimuli produce a progressive increase of the sensation of warmth and heat and eventually that of a burning pain. The pain resulting from repetitive warm stimuli is mediated by summated C fibre responses. To shed more light on the cortical changes associated with pain during repeated subnoxious warm stimulation, we analysed magnetoencephalographic (MEG) evoked fields in eleven subjects during application of repetitive warm laser stimuli to the dorsum of the right hand. One set of stimuli encompassed 10 laser pulses occurring at 2.5 s intervals. Parameters of laser stimulation were optimised to elicit a pleasant warm sensation upon a single stimulus with a rise of skin temperature after repeated stimulation not exceeding the threshold of C mechano-heat fibres. Subjects reported a progressive increase of the intensity of heat and burning pain during repeated laser stimulation in spite of only mild (4.8°C) increase of skin temperature from the first stimulus to the tenth stimulus. The mean reaction time, evaluated in six subjects, was 1.33 s, confirming involvement of C fibres. The neuromagnetic fields were modelled by five equivalent source dipoles located in the occipital cortex, cerebellum, posterior cingulate cortex, and left and right operculo-insular cortex. The only component showing statistically significant changes during repetitive laser stimulation was the late component of the contralateral operculo-insular source peaking at 1.05 s after stimulus onset. The amplitude increases of the late component of the contralateral operculo-insular source dipole correlated with the subjects' numerical ratings of warmth and pain. Results point to a pivotal role of the contralateral operculo-insular region in processing of C-fibre mediated pain during repeated subnoxious laser stimulation.     

Subthalamic local field beta oscillations during ongoing deep brain stimulation in Parkinson's disease in hyperacute and chronic phases

In the past years, local field potential (LFP) signals recorded from the subthalamic nucleus (STN) in patients undergoing deep brain stimulation (DBS) for Parkinson's disease (PD) disclosed that DBS has a controversial effect on STN beta oscillations recorded 2-7 days after surgery for macroelectrode implantation. Nothing is known about these DBS-induced oscillatory changes 30 days after surgery. We recorded STN LFPs during ongoing DBS in 7 patients with PD, immediately (hyperacute phase) and 30 days (chronic phase) after surgery. STN LFP recordings showed stationary intranuclear STN beta LFP activity in hyperacute and chronic phases, confirming that beta peaks were also present in chronic recordings. Power spectra of nuclei with significant beta activity (54% of the sample) showed that it decreased significantly during DBS (p=0.021) under both recording conditions. The time course of beta activity showed more evident DBS-induced changes in the chronic than in the hyperacute phase (p=0.014). DBS-induced changes in STN beta LFPs in patients undergoing DBS in chronic phase provide useful information for developing a new neurosignal-controlled adaptive DBS system.     

Moving toward a generalizable application of central thalamic deep brain stimulation for support of forebrain arousal regulation in the severely injured brain

This review considers the challenges ahead for developing a generalizable strategy for the use of central thalamic deep brain stimulation (CT/DBS) to support arousal regulation mechanisms in the severely injured brain. Historical efforts to apply CT/DBS to patients with severe brain injuries and a proof-of-concept result from a single-subject study are discussed. Circuit and cellular mechanisms underlying the recovery of consciousness are considered for their relevance to the application of CT/DBS, to improve consciousness and cognition in nonprogressive brain injuries. Finally, directions for development, and testing of generalizable criteria for CT/DBS are suggested, which aim to identify neuronal substrates and behavioral profiles that may optimally benefit from support of arousal regulation mechanisms.     

Uniqueness of local myocardial strain patterns with respect to activation time and contractility of the failing heart: a computational study

Background

Myocardial deformation measured by strain is used to detect electro-mechanical abnormalities in cardiac tissue. Estimation of myocardial properties from regional strain patterns when multiple pathologies are present is therefore a promising application of computer modelling. However, if different tissue properties lead to indistinguishable strain patterns (‘degeneracy’), the applicability of any such method will be limited. We investigated whether estimation of local activation time (AT) and contractility from myocardial strain patterns is theoretically possible.

Methods

For four different global cardiac pathologies local myocardial strain patterns for 1025 combinations of AT and contractility were simulated with a computational model (CircAdapt). For each strain pattern, a cohort of similar patterns was found within estimated measurement error using the sum of least-squared differences. Cohort members came from (1) the same pathology only, and (2) all four pathologies. Uncertainty was calculated as accuracy and precision of cohort members in parameter space. Connectedness within the cohorts was also studied.

Results

We found that cohorts drawn from one pathology had parameters with adjacent values although their distribution was neither constant nor symmetrical. In comparison cohorts drawn from four pathologies had disconnected components with drastically different parameter values and accuracy and precision values up to three times higher.

Conclusions

Global pathology must be known when extracting AT and contractility from strain patterns, otherwise degeneracy occurs causing unacceptable uncertainty in derived parameters.

     

Challenges to deep brain stimulation: a pragmatic response to ethical, fiscal, and regulatory concerns

In response to the early success of deep brain stimulation, we offer some common-sense strategies to sustain the work, addressing the need to do so in a fiscally workable, ethically transparent, and scientifically informed manner. After delineating major threats, we will suggest reforms in both the legislative and regulatory spheres that might remediate these challenges. We will recommend (1) revisions to the Bayh-Dole Act of 1980, which governs intellectual property exchange resulting from federally funded research; (2) revisions to the Association of American Medical Colleges recommendations concerning the management of conflicts of interest when scientists with an intellectual property interest participate in clinical research in tandem; (3) revisions to the Food and Drug Administration's pre-market approval process for new devices, including a proposal for a mini-investigational device exemption; and (4) the establishment of a public-private partnership to build ethical and sustainable synergies between the scientific community, industry, and government that would foster discovery and innovation.     

B cell activation triggered by the formation of the small receptor cluster: a computational study

We proposed a spatially extended model of early events of B cell receptors (BCR) activation, which is based on mutual kinase-receptor interactions that are characteristic for the immune receptors and the Src family kinases. These interactions lead to the positive feedback which, together with two nonlinearities resulting from the double phosphorylation of receptors and Michaelis-Menten dephosphorylation kinetics, are responsible for the system bistability. We demonstrated that B cell can be activated by a formation of a tiny cluster of receptors or displacement of the nucleus. The receptors and Src kinases are activated, first locally, in the locus of the receptor cluster or the region where the cytoplasm is the thinnest. Then the traveling wave of activation propagates until activity spreads over the whole cell membrane. In the models in which we assume that the kinases are free to diffuse in the cytoplasm, we found that the fraction of aggregated receptors, capable to initiate B cell activation decreases with the decreasing thickness of cytoplasm and decreasing kinase diffusion. When kinases are restricted to the cell membrane - which is the case for most of the Src family kinases - even a cluster consisting of a tiny fraction of total receptors becomes activatory. Interestingly, the system remains insensitive to the modest changes of total receptor level. The model provides a plausible mechanism of B cells activation due to the formation of small receptors clusters collocalized by binding of polyvalent antigens or arising during the immune synapse formation.     

Calcium dynamics in a bergmann glial cell during metabotropic activation: a simulation study

A mathematical model of calcium dynamics in a Bergmann cell of the cerebellum is proposed. The model adequately describes the experimentally observed behavior of the prototype, including the shape and time-scale of Ca²⁺ responses to single and repetitive metabotropic stimuli and the changes of Ca²⁺ transients caused by inhibition of Ca²⁺ uptake into the store. By means of the model, the role of calcium pumps in regulating the cytoplasmic Ca²⁺ concentration is studied. It is found that the store dimension evaluated by stimulation has the order of magnitude of tens of nanometers, and the Ca²⁺ concentration in the store is about 10 μM.     

Modeling shifts in the rate and pattern of subthalamopallidal network activity during deep brain stimulation

Deep brain stimulation (DBS) of the subthlamic nucleus (STN) represents an effective treatment for medically refractory Parkinson’s disease; however, understanding of its effects on basal ganglia network activity remains limited. We constructed a computational model of the subthalamopallidal network, trained it to fit in vivo recordings from parkinsonian monkeys, and evaluated its response to STN DBS. The network model was created with synaptically connected single compartment biophysical models of STN and pallidal neurons, and stochastically defined inputs driven by cortical beta rhythms. A least mean square error training algorithm was developed to parameterize network connections and minimize error when compared to experimental spike and burst rates in the parkinsonian condition. The output of the trained network was then compared to experimental data not used in the training process. We found that reducing the influence of the cortical beta input on the model generated activity that agreed well with recordings from normal monkeys. Further, during STN DBS in the parkinsonian condition the simulations reproduced the reduction in GPi bursting found in existing experimental data. The model also provided the opportunity to greatly expand analysis of GPi bursting activity, generating three major predictions. First, its reduction was proportional to the volume of STN activated by DBS. Second, GPi bursting decreased in a stimulation frequency dependent manner, saturating at values consistent with clinically therapeutic DBS. And third, ablating STN neurons, reported to generate similar therapeutic outcomes as STN DBS, also reduced GPi bursting. Our theoretical analysis of stimulation induced network activity suggests that regularization of GPi firing is dependent on the volume of STN tissue activated and a threshold level of burst reduction may be necessary for therapeutic effect.     

Rhythmic after-discharge in the cat cerebral cortex and deep brain structures during stimulation of the claustrum

Effects of stimulation of the claustrum and caudate nucleus in the neocortex and various deep brain structures were studied in acute experiments on unanesthetized cats immobilized with tubocurarine. A rhythmic after discharge appeared in neocortical areas 4–7 and 18 (according to Reinoso-Suarez' atlas), and also in the caudate nucleus and various parts of the thalamus. A similar discharge also was observed in the claustrum itself. Diencephalic brain section at the level of the ventral anterior nucleus weakened but did not completely abolish the cortical rhythmic after-discharge in the anterior regions of the neocortex evoked by stimulation of the claustrum. This discharge was completely blocked after sagittal brain section between the claustrum and the rest of the thalamus.I. S. Beritashvili Institute of Physiology, Academy of Sciences of the Georgian SSR, Tbilisi. Translated from Neirofiziologiya, Vol. 15, No. 2, pp. 121–127, March–April, 1983.     

Wide confocal cytometry: a new approach to study proteomic and structural changes in the cell nucleus during the cell cycle

Wide-confocal-cytometry (WCC) is a new method developed in this paper that uses a standard confocal system to gather quantitative information on contents and fine structural details of cells. The system is operated under conditions of non-confocality, in order to capture the maximum amount of light emitted by the specimen (comparable to LSC). After analysis of macromolecule content (DNA, RNA, specific proteins, lipids, etc.), cells can be sampled using conventional confocal microscopy. We analyzed the illumination and acquiring capabilities of WCC. The quantitative power of WCC was validated by analysis of cell cycle stage in Hela cells, looking at DNA content and markers for S phase and mitosis. As an example of the potential of this methodology we have documented changes in cell nucleus during the cell cycle. After mitosis the cell nucleus changes its shape from elongated to ellipsoid and remains constant until G2. This change is associated with nuclear volume increase. As nuclear volume increases, chromatin becomes decondensed in an isometric manner, probably due to the increase in gene expression and factors necessary for RNA metabolism.     

Enthalpic and entropic contributions in the transesterification of sucrose: computational study of lipases and subtilisin

Transesterification of sucrose with fatty acids catalyzed by subtilisin Carlsberg occurs with regioselectivity that is different from that in lipases. Thermomyces lanuginosus lipase (TlL) and Candida antarctica lipase B (CALB) catalyze synthesis at positions 6 and 6', with differing abilities, while subtilisin catalysis leads to the 1'-acylated sucrose. The catalytic machinery in lipases is approximately mirrored in subtilisins but different pocket morphologies including size, shape, and rearrangement of the catalytic elements underlies the differing regioselectivities. The thermodynamic consequences of these differences on the above reactions have been explored systematically using computational methods, determining the free energies of interaction of the putative transition-state adducts. Analysis of the conformers with the lowest transition state energies (protein-ligand interactions and vibrational entropy contributions) indicates that enthalpic factors control specificities in lipases while entropic factors are more important in subtilisin.     

Relative contributions of local and regional factors to species richness and total density of butterflies and moths in semi-natural grasslands

Metapopulation theory predicts that species richness and total population density of habitat specialists increase with increasing area and regional connectivity of the habitat. To test these predictions, we examined the relative contributions of habitat patch area, connectivity of the regional habitat network and local habitat quality to species richness and total density of butterflies and day-active moths inhabiting semi-natural grasslands. We studied butterflies and moths in 48 replicate landscapes situated in southwest Finland, including a focal patch and the surrounding network of other semi-natural grasslands within a radius of 1.5 km from the focal patch. By applying the method of hierarchical partitioning, which can distinguish between independent and joint contributions of individual explanatory variables, we observed that variables of the local habitat quality (e.g. mean vegetation height and nectar plant abundance) generally showed the highest independent effect on species richness and total density of butterflies and moths. Habitat area did not show a significant independent contribution to species richness and total density of butterflies and moths. The effect of habitat connectivity was observed only for total density of the declining butterflies and moths. These observations indicate that the local habitat quality is of foremost importance in explaining variation in species richness and total density of butterflies and moths. In addition, declining butterflies and moths have larger populations in well-connected networks of semi-natural grasslands. Our results suggest that, while it is crucial to maintain high-quality habitats by management, with limited resources it would be appropriate to concentrate grassland management and restoration to areas with well-connected grassland networks in which the declining species currently have their strongest populations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Inhibition of cell proliferation by interferons. Relative contributions of changes in protein synthesis and breakdown to growth control of human lymphoblastoid cells.

Treatment of the Daudi line of human lymphoblastoid cells with concentrations of human interferons within the physiological range progressively inhibits cell proliferation over 1-4 days. Rigorous measurement of the overall rate of protein synthesis during this period, using a concentration of [3H]phenylalanine sufficient to equalize the specific radioactivity of intracellular and extracellular precursor pools, shows that protein synthesis becomes progressively inhibited as the growth inhibition develops. There is a strong correlation between inhibition of amino acid incorporation and inhibition of cell proliferation. In contrast, we find no evidence for any increase in protein degradation rate under these conditions. These results suggest that interferon treatment of susceptible cells can inhibit protein synthesis even in the absence of virus infection and that this inhibition is of a sufficient magnitude to account for the anti-proliferative effect.     

Thalamic deep brain stimulation paradigm to reduce consciousness: Cortico-striatal dynamics implicated in mechanisms of consciousness

Anesthetic manipulations provide much-needed causal evidence for neural correlates of consciousness, but non-specific drug effects complicate their interpretation. Evidence suggests that thalamic deep brain stimulation (DBS) can either increase or decrease consciousness, depending on the stimulation target and parameters. The putative role of the central lateral thalamus (CL) in consciousness makes it an ideal DBS target to manipulate circuit-level mechanisms in cortico-striato-thalamic (CST) systems, thereby influencing consciousness and related processes. We used multi-microelectrode DBS targeted to CL in macaques while recording from frontal, parietal, and striatal regions. DBS induced episodes of abnormally long, vacant staring with low-frequency oscillations here termed vacant, perturbed consciousness (VPC). DBS modulated VPC likelihood in a frequency-specific manner. VPC events corresponded to decreases in measures of neural complexity (entropy) and integration (Φ*), proposed indices of consciousness, and substantial changes to communication in CST circuits. During VPC, power spectral density and coherence at low frequencies increased across CST circuits, especially in thalamo-parietal and cortico-striatal pathways. Decreased consciousness and neural integration corresponded to shifts in cortico-striatal network configurations that dissociated parietal and subcortical structures. Overall, the features of VPC and implicated networks were similar to those of absence epilepsy. As this same multi-microelectrode DBS method–but at different stimulation frequencies–can also increase consciousness in anesthetized macaques, it can be used to flexibly address questions of consciousness with limited confounds, as well as inform clinical investigations of other consciousness disorders.     

Modulatory effects of parallel fiber and molecular layer interneuron synaptic activity on purkinje cell responses to ascending segment input: a modeling study

Based on anatomical, physiological, and model-based studies, it has been proposed that synapses associated with the ascending segment of granule cell axons provide the principle excitatory drive on Purkinje cells which is then modulated by the more numerous parallel fiber synapses. In this study we have evaluated this idea using a detailed compartmental model of a cerebellar Purkinje cell by providing identical ascending segment synaptic inputs during different levels of random parallel fiber and molecular interneuron input. Results suggest that background inputs from parallel fibers and molecular layer interneurons can have a substantial effect on the response of Purkinje cells to ascending segment inputs. Interestingly, these effects are not reflected in the average firing rate of the Purkinje cell and are thus entirely dendritic in effect. These results are considered in the context of the known segregated spatial distribution of the parallel fibers and ascending segment synapses and a new hypothesis concerning the functional organization of cerebellar cortical circuitry.     

Sensitivity of vertical jumping performance to changes in muscle stimulation onset times: a simulation study

The effect of muscle stimulation dynamics on the sensitivity of jumping achievement to variations in timing of muscle stimulation onsets was investigated. Vertical squat jumps were simulated using a forward dynamic model of the human musculoskeletal system. The model calculates the motion of body segments corresponding to STIM(t) of six major muscle groups of the lower extremity, where STIM is muscle stimulation level. For each muscle, STIM was allowed to switch “on” only once. The subsequent rise of STIM to its maximum was described using a sigmoidal curve, the dynamics of which was expressed as rise time (RT). For different values of stimulation RT, the optimal set of onset times was determined using dynamic optimization with height reached by the center of mass as performance criterion. Subsequently, 200 jumps were simulated in which the optimal muscle stimulation onset times were perturbed by adding to each a small number taken from a Gaussian-distributed set of pseudo-random numbers. The distribution of heights achieved in these perturbed jumps was used to quantify the sensitivity of jump height to variations in timing of muscle stimulation onsets. It was found that with increasing RT, the sensitivity of jump height to timing of stimulation onset times decreased. To try and understand this finding, a post-hoc analysis was performed on the perturbed jumps. Jump height was most sensitive to errors in the delay between stimulation onset times of proximal muscles and stimulation onset times of plantar flexors. It is explained how errors in this delay cause aberrations in the configuration of the system, which are regenerative and lead to relatively large jump height deficits. With increasing RT, the initial aberrations due to erroneous timing of muscle stimulation are smaller, and the regeneration is less pronounced. Finally, it is speculated that human subjects decrease or increase RT depending on the relative importance of different performance criteria. Received: 16 February 1998 / Accepted in revised form: 1 March 1999     

Insights into mechanisms of intestinal segmentation in guinea pigs: a combined computational modeling and in vitro study

Segmentation in the guinea pig small intestine consists of a number of discrete motor patterns including rhythmic stationary contractions that occur episodically at specific locations along the intestine. The enteric nervous system regulates segmentation, but the exact circuit is unknown. Using simple computer models, we investigated possible circuits. Our computational model simulated the mean neuron firing rate in the feedforward ascending and descending reflex pathways. A stimulus-evoked pacemaker was located in the afferent pathway or in a feedforward pathway. Output of the feedforward pathways was fed into a simple model to determine the response of the muscle. Predictions were verified in vitro by using guinea pig jejunum, in which segmentation was induced with luminal fatty acid. In the computational model, local stimuli produced an oral contraction and anal dilation, similar to in vitro responses to local distension, but did not produce segmentation. When the stimulus was distributed, representing a nutrient load, the result was either a tonic response or globally synchronized oscillations. However, when we introduced local variations in synaptic coupling, stationary contractions occurred around these locations. This predicts that severing the ascending and descending pathways will induce stationary contractions. An acute lesion in our in vitro model significantly increased the number of stationary contractions immediately oral and anal to the lesion. Our results suggest that spatially localized rhythmic contractions arise from a local imbalance between ascending excitatory and descending inhibitory muscle inputs and require a distributed stimulus and a rhythm generator in the afferent pathway.     

A second-generation computational modeling of cardiac electrophysiology: response of action potential to ionic concentration changes and metabolic inhibition

Background

Cardiac arrhythmias are becoming one of the major health care problem in the world, causing numerous serious disease conditions including stroke and sudden cardiac death. Furthermore, cardiac arrhythmias are intimately related to the signaling ability of cardiac cells, and are caused by signaling defects. Consequently, modeling the electrical activity of the heart, and the complex signaling models that subtend dangerous arrhythmias such as tachycardia and fibrillation, necessitates a quantitative model of action potential (AP) propagation. Yet, many electrophysiological models, which accurately reproduce dynamical characteristic of the action potential in cells, have been introduced. However, these models are very complex and are very time consuming computationally. Consequently, a large amount of research is consecrated to design models with less computational complexity.

Results

This paper is presenting a new model for analyzing the propagation of ionic concentrations and electrical potential in space and time. In this model, the transport of ions is governed by Nernst-Planck flux equation (NP), and the electrical interaction of the species is described by a new cable equation. These set of equations form a system of coupled partial nonlinear differential equations that is solved numerically. In the first we describe the mathematical model. To realize the numerical simulation of our model, we proceed by a finite element discretization and then we choose an appropriate resolution algorithm.

Conclusions

We give numerical simulations obtained for different input scenarios in the case of suicide substrate reaction which were compared to those obtained in literature. These input scenarios have been chosen so as to provide an intuitive understanding of dynamics of the model. By accessing time and space domains, it is shown that interpreting the electrical potential of cell membrane at steady state is incorrect. This model is general and applies to ions of any charge in space and time domains. The results obtained show a complete agreement with literature findings and also with the physical interpretation of the phenomenon. Furthermore, various numerical experiments are presented to confirm the accuracy, efficiency and stability of the proposed method. In particular, we show that the scheme is second-order accurate in space.