Search

Spatial stochastic molecular simulations in biology are limited by the intense computation required to track molecules in space either in a discrete time or discrete space framework, which has led to the development of parallel methods that can take advantage of the power of modern supercomputers in recent years. We systematically test suggested components of stochastic reaction-diffusion operator splitting in the literature and discuss their effects on accuracy. We introduce an operator splitting implementation for irregular meshes that enhances accuracy with minimal performance cost. We test a range of models in small-scale MPI simulations from simple diffusion models to realistic biological models and find that multi-dimensional geometry partitioning is an important consideration for optimum performance. We demonstrate performance gains of 1-3 orders of magnitude in the parallel implementation, with peak performance strongly dependent on model specification. [ABSTRACT FROM AUTHOR]

Synchronous spiking during cerebellar tasks has been observed across Purkinje cells: however, little is known about the intrinsic cellular mechanisms responsible for its initiation, cessation and stability. The Phase Response Curve (PRC), a simple input-output characterization of single cells, can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically. Recently, the PRC theory applied to cerebellar Purkinje cells revealed that these behave as phase-independent integrators at low firing rates, and switch to a phase-dependent mode at high rates. Given the implications for computation and information processing in the cerebellum and the possible role of synchrony in the communication with its post-synaptic targets, we further explored the firing rate dependency of the PRC in Purkinje cells. We isolated key factors for the experimental estimation of the PRC and developed a closed-loop approach to reliably compute the PRC across diverse firing rates in the same cell. Our results show unambiguously that the PRC of individual Purkinje cells is firing rate dependent and that it smoothly transitions from phase independent integrator to a phase dependent mode. Using computational models we show that neither channel noise nor a realistic cell morphology are responsible for the rate dependent shift in the phase response curve. [ABSTRACT FROM AUTHOR]

Abstract Blood oxygenation level dependent contrast (BOLD) functional MRI (fMRI) responses, in a 7-T magnet, were observed in the cerebellum of alpha-chloralose anaesthetized rats in response to innocuous electrical stimulation of a forepaw or hindpaw. The responses were imaged in both coronal and sagittal slices which allowed for a clear delineation and localization of the observed activations. We demonstrate the validity of our fMRI protocol by imaging the responses in somatosensory cortex to the same stimuli and by showing reproducibility of the cerebellar responses. Widespread bilateral activations were found with mainly a patchy and mediolateral band organization, more pronounced ipsilaterally. Possible parasagittal bands were observed only in contralateral lobule VI. There was no overlap between the cerebellar activations caused by forepaw and hindpaw stimuli. The overall horizontal organization of these responses was quite remarkable. For both stimulation paradigms most of the activation patches were positioned in either a rostral or caudal broad plane running anteroposteriorly through both anterior and posterior cerebellum. The rostral planes were completely separated, with the forepaw activation closer to the surface, while the caudal plane was common to both stimulation protocols. We relate our findings to the known projection patterns of spinocerebellar and cuneocerebellar mossy fibres, and to human fMRI studies. [ABSTRACT FROM AUTHOR]

A letter to the editor is presented in response to the article "What Makes Doctors Use Computers?," in the August 1984 issue.

Acute ionizing radiation (IR) causes severe DNA damage, leading to cell cycle arrest, cell death, and activation of the innate immune system. The role and signaling pathway of stimulator of interferon genes (STING) in IR-induced tissue damage and cell death are not well understood. This study revealed that STING is crucial for promoting apoptosis in response to DNA damage caused by acute IR both in vitro and in vivo. STING binds to poly (ADP‒ribose) (PAR) produced by activated poly (ADP‒ribose) polymerase-1 (PARP1) upon IR. Compared with that in WT cells, apoptosis was suppressed in Sting gt-/gt- cells. Excessive PAR production by PARP1 due to DNA damage enhances STING phosphorylation, and inhibiting PARP1 reduces cell apoptosis after IR. In vivo, IR-induced crypt cell death was significantly lower in Sting gt-/gt- mice or with low-dose PARP1 inhibitor, PJ34, resulting in substantial resistance to abdominal irradiation. STING deficiency or inhibition of PARP1 function can reduce the expression of the proapoptotic gene PUMA, decrease the localization of Bax on the mitochondrial membrane, and thus reduce cell apoptosis. Our findings highlight crucial roles for STING and PAR in the IR-mediated induction of apoptosis, which may have therapeutic implications for controlling radiation-induced apoptosis or acute radiation symptoms. STING responds to acute ionizing radiation-mediated DNA damage by directly binding to poly (ADP-ribose) (PAR) produced by activated poly (ADP-ribose) polymerase-1 (PARP1), and mainly induces cell apoptosis through Puma-Bax interaction. STING deficiency or reduced production of PAR protected mice against Acute Radiation Syndrome., Competing Interests: Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Experimental evidence showed that an increase in intracellular chloride concentration [Formula: see text] caused by gamma-aminobutyric acid (GABA) input can promote epileptic firing activity, but the actual mechanisms remain elusive. Here in this theoretical work, we show that influx of chloride and concomitant bicarbonate ion [Formula: see text] efflux upon GABA receptor activation can induce epileptic firing activity by transition of GABA from inhibition to excitation. We analyzed the intrinsic property of neuron firing states as a function of [Formula: see text] We found that as [Formula: see text] increases, the system exhibits a saddle-node bifurcation, above which the neuron exhibits a spectrum of intensive firing, periodic bursting interrupted by depolarization block (DB) state, and eventually a stable DB through a Hopf bifurcation. We demonstrate that only GABA stimuli together with [Formula: see text] efflux can switch GABA's effect to excitation which leads to a series of seizure-like events (SLEs). Exposure to a low [Formula: see text] can drive neurons with high concentrations of [Formula: see text] downward to lower levels of [Formula: see text], during which it could also trigger SLEs depending on the exchange rate with the bath. Our analysis and simulation results show how the competition between GABA stimuli-induced accumulation of [Formula: see text] and [Formula: see text] application-induced decrease of [Formula: see text] regulates the neuron firing activity, which helps to understand the fundamental ionic dynamics of SLE., Competing Interests: The authors declare no competing financial interests., (Copyright © 2024 Liu et al.)

Most central neurons have intricately branched dendritic trees that integrate massive numbers of synaptic inputs. Intrinsic active mechanisms in dendrites can be heterogeneous and be modulated in a branch-specific way. However, it remains poorly understood how heterogeneous intrinsic properties contribute to processing of synaptic input. We propose the first computational model of the cerebellar Purkinje cell with dendritic heterogeneity, in which each branch is an individual unit and is characterized by its own set of ion channel conductance densities. When simultaneously activating a cluster of parallel fiber synapses, we measure the peak amplitude of a response and observe how changes in P-type calcium channel conductance density shift the dendritic responses from a linear one to a bimodal one including dendritic calcium spikes and vice-versa. These changes relate to the morphology of each branch. We show how dendritic calcium spikes propagate and how Kv4.3 channels block spreading depolarization to nearby branches., Competing Interests: The authors declare no competing interests., (© 2024 The Author(s).)

Vesicles carry out many essential functions within cells through the processes of endocytosis, exocytosis, and passive and active transport. This includes transporting and delivering molecules between different parts of the cell, and storing and releasing neurotransmitters in neurons. To date, computational simulation of these key biological players has been rather limited and has not advanced at the same pace as other aspects of cell modeling, restricting the realism of computational models. We describe a general vesicle modeling tool that has been designed for wide application to a variety of cell models, implemented within our software STochastic Engine for Pathway Simulation (STEPS), a stochastic reaction-diffusion simulator that supports realistic reconstructions of cell tissue in tetrahedral meshes. The implementation is validated in an extensive test suite, parallel performance is demonstrated in a realistic synaptic bouton model, and example models are visualized in a Blender extension module., (© 2024. The Author(s).)

The fetal phenotype of MPDZ-associated congenital hydrocephalus type 2 with or without brain or eye anomalies (HYC2) (OMIM 615219) is not well described in the literature. The present case shows not previously published clinical fetal features that are detected during routine second trimester ultrasound screening at 21 weeks of gestation such as bilateral ventriculomegaly, lean cavum septum pellucidum, suspicion of hypoplastic corpus callosum, and suspicion of gyration disorder with normal fossa posterior. Combination of clinical features and a gene panel for congenital malformation syndromes detected a homozygous, likely pathogenic nonsense variant in the MPDZ gene. HYC2 is a rare autosomal recessive disorder with prenatal onset. Clinical presentation is highly variable, varying from stillbirth and severe neurodevelopmental problems with death in infancy to adult patients. Other reported associated congenital anomalies are mainly heart defects and ophthalmologic abnormalities. The present case so far is the first prenatally well described case of HYC2 in an ongoing pregnancy., (© 2024 John Wiley & Sons Ltd.)

The cerebellum has been a popular topic for theoretical studies because its structure was thought to be simple. Since David Marr and James Albus related its function to motor skill learning and proposed the Marr-Albus cerebellar learning model, this theory has guided and inspired cerebellar research. In this review, we summarize the theoretical progress that has been made within this framework of error-based supervised learning. We discuss the experimental progress that demonstrates more complicated molecular and cellular mechanisms in the cerebellum as well as new cell types and recurrent connections. We also cover its involvement in diverse non-motor functions and evidence of other forms of learning. Finally, we highlight the need to explain these new experimental findings into an integrated cerebellar model that can unify its diverse computational functions., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

We investigate the relationship between primary dendrite selection of Purkinje cells and migration of their presynaptic partner granule cells during early cerebellar development. During postnatal development, each Purkinje cell grows more than three dendritic trees, from which a primary tree is selected for development, whereas the others completely retract. Experimental studies suggest that this selection process is coordinated by physical and synaptic interactions with granule cells, which undergo a massive migration at the same time. However, technical limitations hinder continuous experimental observation of multiple cell populations. To explore possible mechanisms underlying this selection process, we constructed a computational model using a new computational framework, NeuroDevSim. The study presents the first computational model that simultaneously simulates Purkinje cell growth and the dynamics of granule cell migrations during the first two postnatal weeks, allowing exploration of the role of physical and synaptic interactions upon dendritic selection. The model suggests that interaction with parallel fibers is important to establish the distinct planar morphology of Purkinje cell dendrites. Specific rules to select which dendritic trees to keep or retract result in larger winner trees with more synaptic contacts than using random selection. A rule based on afferent synaptic activity was less effective than rules based on dendritic size or numbers of synapses., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Kato, De Schutter. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The Neural Development Simulator, NeuroDevSim, is a Python module that simulates the most important aspects of brain development: morphological growth, migration, and pruning. It uses an agent-based modeling approach inherited from the NeuroMaC software. Each cycle has agents called fronts execute model-specific code. In the case of a growing dendritic or axonal front, this will be a choice between extension, branching, or growth termination. Somatic fronts can migrate to new positions and any front can be retracted to prune parts of neurons. Collision detection prevents new or migrating fronts from overlapping with existing ones. NeuroDevSim is a multi-core program that uses an innovative shared memory approach to achieve parallel processing without messaging. We demonstrate linear strong parallel scaling up to 96 cores for large models and have run these successfully on 128 cores. Most of the shared memory parallelism is achieved without memory locking. Instead, cores have only write privileges to private sections of arrays, while being able to read the entire shared array. Memory conflicts are avoided by a coding rule that allows only active fronts to use methods that need writing access. The exception is collision detection, which is needed to avoid the growth of physically overlapping structures. For collision detection, a memory-locking mechanism was necessary to control access to grid points that register the location of nearby fronts. A custom approach using a serialized lock broker was able to manage both read and write locking. NeuroDevSim allows easy modeling of most aspects of neural development for models simulating a few complex or thousands of simple neurons or a mixture of both., Code Available at: https://github.com/CNS-OIST/NeuroDevSim., Competing Interests: The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 De Schutter.)