Modeling the Emergence of Circadian Rhythms in a Clock Neuron Network

Circadian rhythms in pacemaker cells persist for weeks in constant darkness, while in other types of cells the molecular oscillations that underlie circadian rhythms damp rapidly under the same conditions. Although much progress has been made in understanding the biochemical and cellular basis of circadian rhythms, the mechanisms leading to damped or self-sustained oscillations remain largely unknown. There exist many mathematical models that reproduce the circadian rhythms in the case of a single cell of the Drosophila fly. However, not much is known about the mechanisms leading to coherent circadian oscillation in clock neuron networks. In this work we have implemented a model for a network of interacting clock neurons to describe the emergence (or damping) of circadian rhythms in Drosophila fly, in the absence of zeitgebers. Our model consists of an array of pacemakers that interact through the modulation of some parameters by a network feedback. The individual pacemakers are described by a well-known biochemical model for circadian oscillation, to which we have added degradation of PER protein by light and multiplicative noise. The network feedback is the PER protein level averaged over the whole network. In particular, we have investigated the effect of modulation of the parameters associated with (i) the control of net entrance of PER into the nucleus and (ii) the non-photic degradation of PER. Our results indicate that the modulation of PER entrance into the nucleus allows the synchronization of clock neurons, leading to coherent circadian oscillations under constant dark condition. On the other hand, the modulation of non-photic degradation cannot reset the phases of individual clocks subjected to intrinsic biochemical noise.

Controlling chaos in wave-particle interactions

We analyze the behavior of a relativistic particle moving under the influence of a uniform magnetic field and a stationary electrostatic wave. We work with a set of pulsed waves that allows us to obtain an exact map for the system. We also use a method of control for near-integrable Hamiltonians that consists of the addition of a small and simple control term to the system. This control term creates invariant tori in phase space that prevent chaos from spreading to large regions, making the controlled dynamics more regular. We show numerically that the control term just slightly modifies the system but is able to drastically reduce chaos with a low additional cost of energy. Moreover, we discuss how the control of chaos and the consequent recovery of regular trajectories in phase space are useful to improve regular particle acceleration.

A Chaotic Approach of Differential Evolution Optimization Applied to Loudspeaker Design Problem

Over the past few years, the field of global optimization has been very active, producing different kinds of deterministic and stochastic algorithms for optimization in the continuous domain. These days, the use of evolutionary algorithms (EAs) to solve optimization problems is a common practice due to their competitive performance on complex search spaces. EAs are well known for their ability to deal with nonlinear and complex optimization problems. Differential evolution (DE) algorithms are a family of evolutionary optimization techniques that use a rather greedy and less stochastic approach to problem solving, when compared to classical evolutionary algorithms. The main idea is to construct, at each generation, for each element of the population a mutant vector, which is constructed through a specific mutation operation based on adding differences between randomly selected elements of the population to another element. Due to its simple implementation, minimum mathematical processing and good optimization capability, DE has attracted attention. This paper proposes a new approach to solve electromagnetic design problems that combines the DE algorithm with a generator of chaos sequences. This approach is tested on the design of a loudspeaker model with 17 degrees of freedom, for showing its applicability to electromagnetic problems. The results show that the DE algorithm with chaotic sequences presents better, or at least similar, results when compared to the standard DE algorithm and other evolutionary algorithms available in the literature.

Dynamical changes from harmonic vibrations of a limited power supply driving a Duffing oscillator

The harmonic oscillations of a Duffing oscillator driven by a limited power supply are investigated as a function of the alternative strength of the rotor. The semi-trivial and non-trivial solutions are derived. We examine the stability of these solutions and then explore the complex behaviors associated with the bifurcations sequences. Interestingly, a 3D diagram provides a global view of the effects of alternate strength on the appearance of chaos and hyperchaos on the system.

Mathematical models of cognitive processes

The research activity carried out during the PhD course was focused on the development of mathematical models of some cognitive processes and their validation by means of data present in literature, with a double aim: i) to achieve a better interpretation and explanation of the great amount of data obtained on these processes from different methodologies (electrophysiological recordings on animals, neuropsychological, psychophysical and neuroimaging studies in humans), ii) to exploit model predictions and results to guide future research and experiments. In particular, the research activity has been focused on two different projects: 1) the first one concerns the development of neural oscillators networks, in order to investigate the mechanisms of synchronization of the neural oscillatory activity during cognitive processes, such as object recognition, memory, language, attention; 2) the second one concerns the mathematical modelling of multisensory integration processes (e.g. visual-acoustic), which occur in several cortical and subcortical regions (in particular in a subcortical structure named Superior Colliculus (SC)), and which are fundamental for orienting motor and attentive responses to external world stimuli. This activity has been realized in collaboration with the Center for Studies and Researches in Cognitive Neuroscience of the University of Bologna (in Cesena) and the Department of Neurobiology and Anatomy of the Wake Forest University School of Medicine (NC, USA). PART 1. Objects representation in a number of cognitive functions, like perception and recognition, foresees distribute processes in different cortical areas. One of the main neurophysiological question concerns how the correlation between these disparate areas is realized, in order to succeed in grouping together the characteristics of the same object (binding problem) and in maintaining segregated the properties belonging to different objects simultaneously present (segmentation problem). Different theories have been proposed to address these questions (Barlow, 1972). One of the most influential theory is the so called “assembly coding”, postulated by Singer (2003), according to which 1) an object is well described by a few fundamental properties, processing in different and distributed cortical areas; 2) the recognition of the object would be realized by means of the simultaneously activation of the cortical areas representing its different features; 3) groups of properties belonging to different objects would be kept separated in the time domain. In Chapter 1.1 and in Chapter 1.2 we present two neural network models for object recognition, based on the “assembly coding” hypothesis. These models are networks of Wilson-Cowan oscillators which exploit: i) two high-level “Gestalt Rules” (the similarity and previous knowledge rules), to realize the functional link between elements of different cortical areas representing properties of the same object (binding problem); 2) the synchronization of the neural oscillatory activity in the γ-band (30-100Hz), to segregate in time the representations of different objects simultaneously present (segmentation problem). These models are able to recognize and reconstruct multiple simultaneous external objects, even in difficult case (some wrong or lacking features, shared features, superimposed noise). In Chapter 1.3 the previous models are extended to realize a semantic memory, in which sensory-motor representations of objects are linked with words. To this aim, the network, previously developed, devoted to the representation of objects as a collection of sensory-motor features, is reciprocally linked with a second network devoted to the representation of words (lexical network) Synapses linking the two networks are trained via a time-dependent Hebbian rule, during a training period in which individual objects are presented together with the corresponding words. Simulation results demonstrate that, during the retrieval phase, the network can deal with the simultaneous presence of objects (from sensory-motor inputs) and words (from linguistic inputs), can correctly associate objects with words and segment objects even in the presence of incomplete information. Moreover, the network can realize some semantic links among words representing objects with some shared features. These results support the idea that semantic memory can be described as an integrated process, whose content is retrieved by the co-activation of different multimodal regions. In perspective, extended versions of this model may be used to test conceptual theories, and to provide a quantitative assessment of existing data (for instance concerning patients with neural deficits). PART 2. The ability of the brain to integrate information from different sensory channels is fundamental to perception of the external world (Stein et al, 1993). It is well documented that a number of extraprimary areas have neurons capable of such a task; one of the best known of these is the superior colliculus (SC). This midbrain structure receives auditory, visual and somatosensory inputs from different subcortical and cortical areas, and is involved in the control of orientation to external events (Wallace et al, 1993). SC neurons respond to each of these sensory inputs separately, but is also capable of integrating them (Stein et al, 1993) so that the response to the combined multisensory stimuli is greater than that to the individual component stimuli (enhancement). This enhancement is proportionately greater if the modality-specific paired stimuli are weaker (the principle of inverse effectiveness). Several studies have shown that the capability of SC neurons to engage in multisensory integration requires inputs from cortex; primarily the anterior ectosylvian sulcus (AES), but also the rostral lateral suprasylvian sulcus (rLS). If these cortical inputs are deactivated the response of SC neurons to cross-modal stimulation is no different from that evoked by the most effective of its individual component stimuli (Jiang et al 2001). This phenomenon can be better understood through mathematical models. The use of mathematical models and neural networks can place the mass of data that has been accumulated about this phenomenon and its underlying circuitry into a coherent theoretical structure. In Chapter 2.1 a simple neural network model of this structure is presented; this model is able to reproduce a large number of SC behaviours like multisensory enhancement, multisensory and unisensory depression, inverse effectiveness. In Chapter 2.2 this model was improved by incorporating more neurophysiological knowledge about the neural circuitry underlying SC multisensory integration, in order to suggest possible physiological mechanisms through which it is effected. This endeavour was realized in collaboration with Professor B.E. Stein and Doctor B. Rowland during the 6 months-period spent at the Department of Neurobiology and Anatomy of the Wake Forest University School of Medicine (NC, USA), within the Marco Polo Project. The model includes four distinct unisensory areas that are devoted to a topological representation of external stimuli. Two of them represent subregions of the AES (i.e., FAES, an auditory area, and AEV, a visual area) and send descending inputs to the ipsilateral SC; the other two represent subcortical areas (one auditory and one visual) projecting ascending inputs to the same SC. Different competitive mechanisms, realized by means of population of interneurons, are used in the model to reproduce the different behaviour of SC neurons in conditions of cortical activation and deactivation. The model, with a single set of parameters, is able to mimic the behaviour of SC multisensory neurons in response to very different stimulus conditions (multisensory enhancement, inverse effectiveness, within- and cross-modal suppression of spatially disparate stimuli), with cortex functional and cortex deactivated, and with a particular type of membrane receptors (NMDA receptors) active or inhibited. All these results agree with the data reported in Jiang et al. (2001) and in Binns and Salt (1996). The model suggests that non-linearities in neural responses and synaptic (excitatory and inhibitory) connections can explain the fundamental aspects of multisensory integration, and provides a biologically plausible hypothesis about the underlying circuitry.

A 2-dimensional computational model to analyze the effects of cellular heterogeinity on cardiac pacemaking

The mechanical action of the heart is made possible in response to electrical events that involve the cardiac cells, a property that classifies the heart tissue between the excitable tissues. At the cellular level, the electrical event is the signal that triggers the mechanical contraction, inducing a transient increase in intracellular calcium which, in turn, carries the message of contraction to the contractile proteins of the cell. The primary goal of my project was to implement in CUDA (Compute Unified Device Architecture, an hardware architecture for parallel processing created by NVIDIA) a tissue model of the rabbit sinoatrial node to evaluate the heterogeneity of its structure and how that variability influences the behavior of the cells. In particular, each cell has an intrinsic discharge frequency, thus different from that of every other cell of the tissue and it is interesting to study the process of synchronization of the cells and look at the value of the last discharge frequency if they synchronized.

Frontal areas contribute to reduced global coordination of resting-state gamma activities in drug-naïve patients with schizophrenia

Schizophrenia has been postulated to involve impaired neuronal cooperation in large-scale neural networks, including cortico-cortical circuitry. Alterations in gamma band oscillations have attracted a great deal of interest as they appear to represent a pathophysiological process of cortical dysfunction in schizophrenia. Gamma band oscillations reflect local cortical activities, and the synchronization of these activities among spatially distributed cortical areas has been suggested to play a central role in the formation of networks. To assess global coordination across spatially distributed brain regions, Omega complexity (OC) in multichannel EEG was proposed. Using OC, we investigated global coordination of resting-state EEG activities in both gamma (30–50 Hz) and below-gamma (1.5–30 Hz) bands in drug-naïve patients with schizophrenia and investigated the effects of neuroleptic treatment. We found that gamma band OC was significantly higher in drug-naïve patients with schizophrenia compared to control subjects and that a right frontal electrode (F3) contributed significantly to the higher OC. After neuroleptic treatment, reductions in the contribution of frontal electrodes to global OC in both bands correlated with the improvement of schizophrenia symptomatology. The present study suggests that frontal brain processes in schizophrenia were less coordinated with activity in the remaining brain. In addition, beneficial effects of neuroleptic treatment were accompanied by improvement of brain coordination predominantly due to changes in frontal regions. Our study provides new evidence of improper intrinsic brain integration in schizophrenia by investigating the resting-state gamma band activity.

Trends in Migrations of Refugees in Serbia

This project aimed to establish a comprehensive long-term projection of refugee migration trends in Serbia which would allow the different organisations engaged in caring for refugees to plan their work more effectively. Mr. Cvetkovic studied first the official records of numbers and origins of refugees and the various definitions used to describe refugees, exiled persons, etc., considering also the indication of nationality, i.e. Serb or Yugoslav, and the future intentions expressed by refugees. He concluded that more than three quarters of the total number of refugees in Serbia wish to remain in Serbia/Federal Republic of Yugoslavia, and only 20% wish to repatriate. He concludes that the situation in relation to the ethnic-national structure of Serbia is extremely complicated, even to the point of chaos. One certainty is that sooner or later a huge majority of the refugees in Serbia today will be granted full citizenship and will then participate in the choice of the political system and the dominant national values and institutions of the country. The experiences of the refugees, as well as of their fellow nationals in Serbia, make it relatively unlikely that they will make rational choices that could produce a balance between civil democracy and national totalitarianism. Insofar as overall political relations within Serbia and its Federation with Montenegro develop along democratic lines, the civic identity of Serbs (as opposed to their national/ethnic status) can be expected to become stronger, and the civil and national to become an integral part of thinking which does not represent a challenge to the democratic state. This would help ensure that the migratory movements of the Serbs and of other national and ethnic groups in the region are motivated by economic rather than political or ethnic reasons.

Duration of Greenland Stadial 22 and ice-gas Δage from counting of annual layers in Greenland NGRIP ice core

High-resolution measurements of chemical impurities and methane concentrations in Greenland ice core samples from the early glacial period allow the extension of annual-layer counted chronologies and the improvement of gas age-ice age difference (Δage) essential to the synchronization of ice core records. We report high-resolution measurements of a 50 m section of the NorthGRIP ice core and corresponding annual layer thicknesses in order to constrain the duration of the Greenland Stadial 22 (GS-22) between Greenland Interstadials (GIs) 21 and 22, for which inconsistent durations and ages have been reported from Greenland and Antarctic ice core records as well as European speleothems. Depending on the chronology used, GS-22 occurred between approximately 89 (end of GI-22) and 83 kyr b2k (onset of GI-21). From annual layer counting, we find that GS-22 lasted between 2696 and 3092 years and was followed by a GI-21 pre-cursor event lasting between 331 and 369 yr. Our layer-based counting agrees with the duration of stadial 22 as determined from the NALPS speleothem record (3250 ± 526 yr) but not with that of the GICC05modelext chronology (2620 yr) or an alternative chronology based on gas-marker synchronization to EPICA Dronning Maud Land ice core. These results show that GICC05modelext overestimates accumulation and/or underestimates thinning in this early part of the last glacial period. We also revise the possible ranges of NorthGRIP Δdepth (5.49 to 5.85 m) and Δage (498 to 601 yr) at the warming onset of GI-21 as well as the Δage range at the onset of the GI-21 precursor warming (523 to 654 yr), observing that temperature (represented by the δ15N proxy) increases before CH4 concentration by no more than a few decades.

Dosimetric impact of geometric errors due to respiratory motion prediction on dynamic multileaf collimator-based four-dimensional radiation delivery

The synchronization of dynamic multileaf collimator (DMLC) response with respiratory motion is critical to ensure the accuracy of DMLC-based four dimensional (4D) radiation delivery. In practice, however, a finite time delay (response time) between the acquisition of tumor position and multileaf collimator response necessitates predictive models of respiratory tumor motion to synchronize radiation delivery. Predicting a complex process such as respiratory motion introduces geometric errors, which have been reported in several publications. However, the dosimetric effect of such errors on 4D radiation delivery has not yet been investigated. Thus, our aim in this work was to quantify the dosimetric effects of geometric error due to prediction under several different conditions. Conformal and intensity modulated radiation therapy (IMRT) plans for a lung patient were generated for anterior-posterior/posterior-anterior (AP/PA) beam arrangements at 6 and 18 MV energies to provide planned dose distributions. Respiratory motion data was obtained from 60 diaphragm-motion fluoroscopy recordings from five patients. A linear adaptive filter was employed to predict the tumor position. The geometric error of prediction was defined as the absolute difference between predicted and actual positions at each diaphragm position. Distributions of geometric error of prediction were obtained for all of the respiratory motion data. Planned dose distributions were then convolved with distributions for the geometric error of prediction to obtain convolved dose distributions. The dosimetric effect of such geometric errors was determined as a function of several variables: response time (0-0.6 s), beam energy (6/18 MV), treatment delivery (3D/4D), treatment type (conformal/IMRT), beam direction (AP/PA), and breathing training type (free breathing/audio instruction/visual feedback). Dose difference and distance-to-agreement analysis was employed to quantify results. Based on our data, the dosimetric impact of prediction (a) increased with response time, (b) was larger for 3D radiation therapy as compared with 4D radiation therapy, (c) was relatively insensitive to change in beam energy and beam direction, (d) was greater for IMRT distributions as compared with conformal distributions, (e) was smaller than the dosimetric impact of latency, and (f) was greatest for respiration motion with audio instructions, followed by visual feedback and free breathing. Geometric errors of prediction that occur during 4D radiation delivery introduce dosimetric errors that are dependent on several factors, such as response time, treatment-delivery type, and beam energy. Even for relatively small response times of 0.6 s into the future, dosimetric errors due to prediction could approach delivery errors when respiratory motion is not accounted for at all. To reduce the dosimetric impact, better predictive models and/or shorter response times are required.

Theta burst transcranial magnetic stimulation is associated with increased EEG synchronization in the stimulated relative to unstimulated cerebral hemisphere

Theta burst transcranial magnetic stimulation (TBS) may induce behavioural changes that outlast the stimulation period. The neurophysiological basis of these behavioural changes are currently under investigation. Given the evidence that cortical information processing relies on transient synchronization and desynchronization of neuronal assemblies, we set out to test whether TBS is associated with changes of neuronal synchronization as assessed by surface EEG. In four healthy subjects one TBS train of 600 pulses (200 bursts, each burst consisting of 3 pulses at 30 Hz, repeated at intervals of 100 ms) was applied over the right frontal eye field and EEG synchronization was assessed in a time-resolved manner over 60 min by using a non-overlapping moving window. For each time step the linear cross-correlation matrix for six EEG channels of the right and for the six homotopic EEG channels of the left hemisphere were computed and their largest eigenvalues used to assess changes of synchronization. Synchronization was computed for broadband EEG and for the delta, theta, alpha, beta and gamma frequency bands. In all subjects EEG synchronization of the stimulated hemisphere was significantly and persistently increased relative to EEG synchronization of the unstimulated hemisphere. This effect occurred immediately after TBS for the theta, alpha, beta and gamma frequency bands and 10-20 min after TBS for broadband and delta frequency band EEG. Our results demonstrate that TBS is associated with increased neuronal synchronization of the cerebral hemisphere ipsilateral to the stimulation site relative to the unstimulated hemisphere. We speculate that enhanced synchronization interferes with cortical information processing and thus may be a neurophysiological correlate of the impaired behavioural performance detected previously.

Cell cycle-dependent regulation of extra-adrenal glucocorticoid synthesis in murine intestinal epithelial cells

Glucocorticoids are anti-inflammatory steroids with important applications in the treatment of inflammatory diseases. Endogenous glucocorticoids are mainly produced by the adrenal glands, although there is increasing evidence for extra-adrenal sources. Recent findings show that intestinal crypt cells produce glucocorticoids, which contribute to the maintenance of intestinal immune homeostasis. Intestinal glucocorticoid synthesis is critically regulated by the transcription factor liver receptor homologue-1 (LRH-1). As expression of steroidogenic enzymes and LRH-1 is restricted to the proliferating cells of the crypts, we aimed to investigate the role of the cell cycle in the regulation of LRH-1 activity and intestinal glucocorticoid synthesis. We here show that either pharmacological or molecular modulation of cell cycle progression significantly inhibited expression of steroidogenic enzymes and synthesis of glucocorticoids in intestinal epithelial cells. Synchronization of intestinal epithelial cells in the cell cycle revealed that expression of steroidogenic enzymes is preferentially induced at the G(1)/S stage. Differentiation of immature intestinal epithelial cells to mature nonproliferating cells also resulted in reduced expression of steroidogenic enzymes. This cell cycle-related effect on intestinal steroidogenesis was found to be mediated through the regulation of LRH-1 transcriptional activity. This mechanism may restrict intestinal glucocorticoid synthesis to the proliferating cells of the crypts.

Trial-to-trial size variability of motor-evoked potentials. A study using the triple stimulation technique

Motor-evoked potentials (MEPs) vary in size from one stimulus to the next. The objective of this study was to determine the cause and source of trial-to-trial MEP size variability. In two experiments involving 10 and 14 subjects, the variability of MEPs to cortical stimulation (cortical-MEPs) in abductor digiti minimi (ADM) and abductor hallucis (AH) was compared to those responses obtained using the triple stimulation technique (cortical-TST). The TST eliminates the effects of motor neuron (MN) response desynchronization and of repetitive MN discharges. Submaximal stimuli were used in both techniques. In six subjects, cortical-MEP variability was compared to that of brainstem-MEP and brainstem-TST. Variability was greater for MEPs than that for TST responses, by approximately one-third. The variability was the same for cortical- and brainstem-MEPs and was similar in ADM and AH. Variability concerned at least 10-15% of the MN pool innervating the target muscle. With the stimulation parameters used, repetitive MN discharges did not influence variability. For submaximal stimuli, approximately two-third of the observed MEP size variability is caused by the variable number of recruited alpha-MNs and approximately one-third by changing synchronization of MN discharges. The source of variability is most likely localized at the spinal segmental level.

BOLD correlates of EEG alpha phase-locking and the fMRI default mode network

Phase locking or synchronization of brain areas is a key concept of information processing in the brain. Synchronous oscillations have been observed and investigated extensively in EEG during the past decades. EEG oscillations occur over a wide frequency range. In EEG, a prominent type of oscillations is alpha-band activity, present typically when a subject is awake, but at rest with closed eyes. The spectral power of alpha rhythms has recently been investigated in simultaneous EEG/fMRI recordings, establishing a wide-range cortico-thalamic network. However, spectral power and synchronization are different measures and little is known about the correlations between BOLD effects and EEG synchronization. Interestingly, the fMRI BOLD signal also displays synchronous oscillations across different brain regions. These oscillations delineate so-called resting state networks (RSNs) that resemble the correlation patterns of simultaneous EEG/fMRI recordings. However, the nature of these BOLD oscillations and their relations to EEG activity is still poorly understood. One hypothesis is that the subunits constituting a specific RSN may be coordinated by different EEG rhythms. In this study we report on evidence for this hypothesis. The BOLD correlates of global EEG synchronization (GFS) in the alpha frequency band are located in brain areas involved in specific RSNs, e.g. the 'default mode network'. Furthermore, our results confirm the hypothesis that specific RSNs are organized by long-range synchronization at least in the alpha frequency band. Finally, we could localize specific areas where the GFS BOLD correlates and the associated RSN overlap. Thus, we claim that not only the spectral dynamics of EEG are important, but also their spatio-temporal organization.

Cortical local and long-range synchronization interplay in human absence seizure initiation

Brain activity relies on transient, fluctuating interactions between segregated neuronal populations. Synchronization within a single and between distributed neuronal clusters reflects the dynamics of these cooperative patterns. Thus absence epilepsy can be used as a model for integrated, large-scale investigation of the emergence of pathological collective dynamics in the brain. Indeed, spike-wave discharges (SWD) of an absence seizure are thought to reflect abnormal cortical hypersynchronization. In this paper, we address two questions: how and where do SWD arise in the human brain? Therefore, we explored the spatio-temporal dynamics of interactions within and between widely distributed cortical sites using magneto-encephalographic recordings of spontaneous absence seizures. We then extracted, from their time-frequency analysis, local synchronization of cortical sources and long-range synchronization linking distant sites. Our analyses revealed a reproducible sequence of 1) long-range desynchronization, 2) increased local synchronization and 3) increased long-range synchronization. Although both local and long-range synchronization displayed different spatio-temporal profiles, their cortical projection within an initiation time window overlap and reveal a multifocal fronto-central network. These observations contradict the classical view of sudden generalized synchronous activities in absence epilepsy. Furthermore, they suggest that brain states transition may rely on multi-scale processes involving both local and distant interactions.