Identification of neural activity based on fMRI data: A simulation study

This study proposes a technique for determining effective connectivity among brain regions which operates at the level of neuronal dynamics. We propose an alternative time-variant dynamic causal model (TV-DCM) where neuronal dynamic activity evolves based on generalized dynamic neural networks (GDNNs). The identification of brain architecture connectivity is carried out based on a least squares criterion and on a global search technique. Computer simulations carried out in the paper show that TV-DCM may provide multiple solutions, i.e. a set of different architectures all of which approximate the data equally well. Numerical comparisons between TV-DCM and DCM are also given. In order to determine the unique causal structure of brain regions, we apply an additional criterion, i.e. an estimation of generalization error, known from the theory of neural networks. Computer simulations also confirm the validity of our techniques.     

Dynamic causal models and autopoietic systems

Dynamic Causal Modelling (DCM) and the theory of autopoietic systems are two important conceptual frameworks. In this review, we suggest that they can be combined to answer important questions about self-organising systems like the brain. DCM has been developed recently by the neuroimaging community to explain, using biophysical models, the non-invasive brain imaging data are caused by neural processes. It allows one to ask mechanistic questions about the implementation of cerebral processes. In DCM the parameters of biophysical models are estimated from measured data and the evidence for each model is evaluated. This enables one to test different functional hypotheses (i.e., models) for a given data set. Autopoiesis and related formal theories of biological systems as autonomous machines represent a body of concepts with many successful applications. However, autopoiesis has remained largely theoretical and has not penetrated the empiricism of cognitive neuroscience. In this review, we try to show the connections that exist between DCM and autopoiesis. In particular, we propose a simple modification to standard formulations of DCM that includes autonomous processes. The idea is to exploit the machinery of the system identification of DCMs in neuroimaging to test the face validity of the autopoietic theory applied to neural subsystems. We illustrate the theoretical concepts and their implications for interpreting electroencephalographic signals acquired during amygdala stimulation in an epileptic patient. The results suggest that DCM represents a relevant biophysical approach to brain functional organisation, with a potential that is yet to be fully evaluated.     

A new method of identifying a systems based on the minimal quadratic discrepancy criterion for biophysics problems

A wide range of biophysical systems are described by nonlinear dynamic models mathematically presented as a set of ordinary differential equations in the Cauchy explicit form: [formula: see text] Fij(X1(t),..,XN(t),t), (i = 1,...,N, j = 1,...,M), where Fij (X1(t), ..., XN(t), t) is a set of basis functions satisfying the Lipschitz condition. We investigate the problem of evaluation of model constants aij (the system identification) using experimental data about the time dependence of the dynamic parameters of the system Xi(t). A new method of system identification for the class of similar nonlinear dynamic models is proposed. It is shown that the problem of identifying an initial nonlinear model can be reduced to the solution of a system of linear equations for the matrix of the dynamic model constants [aj]i. It is proposed to determine the set of dynamic model constants aij using the criterion of minimal quadratic discrepancy for the time dependence of the set of dynamic parameters Xi(t). An important special case of the nonlinear model, the quadratic model, is considered. Test problems of identification using this method are presented for two nonlinear systems: the Van der Pol type multiparametric nonlinear oscillator and the strange attractor of Ressler, a widely known example of dynamic systems showing the stochastic behavior.     

The research of constructing dynamic cognition model based on brain network

Estimating the functional interactions and connections between brain regions to corresponding process in cognitive, behavioral and psychiatric domains is a central pursuit for understanding the human connectome. Few studies have examined the effects of dynamic evolution on cognitive processing and brain activation using brain network model in scalp electroencephalography (EEG) data. Aim of this study was to investigate the brain functional connectivity and construct dynamic programing model from EEG data and to evaluate a possible correlation between topological characteristics of the brain connectivity and cognitive evolution processing. Here, functional connectivity between brain regions is defined as the statistical dependence between EEG signals in different brain areas and is typically determined by calculating the relationship between regional time series using wavelet coherence. We present an accelerated dynamic programing algorithm to construct dynamic cognitive model that we found that spatially distributed regions coherence connection difference, the topologic characteristics with which they can transfer information, producing temporary network states. Our findings suggest that brain dynamics give rise to variations in complex network properties over time after variation audio stimulation, dynamic programing model gives the dynamic evolution processing at different time and frequency. In this paper, by applying a new construct approach to understand whole brain network dynamics, firstly, brain network is constructed by wavelet coherence, secondly, different time active brain regions are selected by network topological characteristics and minimum spanning tree. Finally, dynamic evolution model is constructed to understand cognitive process by dynamic programing algorithm, this model is applied to the auditory experiment, results showed that, quantitatively, more correlation was observed after variation audio stimulation, the EEG function connection dynamic evolution model on cognitive processing is feasible with wavelet coherence EEG recording.     

Optimizing experimental design for comparing models of brain function

This article presents the first attempt to formalize the optimization of experimental design with the aim of comparing models of brain function based on neuroimaging data. We demonstrate our approach in the context of Dynamic Causal Modelling (DCM), which relates experimental manipulations to observed network dynamics (via hidden neuronal states) and provides an inference framework for selecting among candidate models. Here, we show how to optimize the sensitivity of model selection by choosing among experimental designs according to their respective model selection accuracy. Using Bayesian decision theory, we (i) derive the Laplace-Chernoff risk for model selection, (ii) disclose its relationship with classical design optimality criteria and (iii) assess its sensitivity to basic modelling assumptions. We then evaluate the approach when identifying brain networks using DCM. Monte-Carlo simulations and empirical analyses of fMRI data from a simple bimanual motor task in humans serve to demonstrate the relationship between network identification and the optimal experimental design. For example, we show that deciding whether there is a feedback connection requires shorter epoch durations, relative to asking whether there is experimentally induced change in a connection that is known to be present. Finally, we discuss limitations and potential extensions of this work.     

Brain function and higher nervous activity

Some actual problems of higher nervous activity have been analysed on the peculiarities of brain states in the course of different types of conditioning and reactivity of the nervous structures which depend on the brain state have been considered. A concept of brain state as a specific mechanism of each temporal connection forming during the learning process has been formulated for the first time. The authors suggest that the brain represents the dynamic system with changeable structure which reveals itself in multitude nervous set constellation during various types of activity. This concept is presented to be the theoretical basis for integral evaluation of functional capacities of central nervous system.     

Modeling disordered regions in proteins using Rosetta

Protein structure prediction methods such as Rosetta search for the lowest energy conformation of the polypeptide chain. However, the experimentally observed native state is at a minimum of the free energy, rather than the energy. The neglect of the missing configurational entropy contribution to the free energy can be partially justified by the assumption that the entropies of alternative folded states, while very much less than unfolded states, are not too different from one another, and hence can be to a first approximation neglected when searching for the lowest free energy state. The shortcomings of current structure prediction methods may be due in part to the breakdown of this assumption. Particularly problematic are proteins with significant disordered regions which do not populate single low energy conformations even in the native state. We describe two approaches within the Rosetta structure modeling methodology for treating such regions. The first does not require advance knowledge of the regions likely to be disordered; instead these are identified by minimizing a simple free energy function used previously to model protein folding landscapes and transition states. In this model, residues can be either completely ordered or completely disordered; they are considered disordered if the gain in entropy outweighs the loss of favorable energetic interactions with the rest of the protein chain. The second approach requires identification in advance of the disordered regions either from sequence alone using for example the DISOPRED server or from experimental data such as NMR chemical shifts. During Rosetta structure prediction calculations the disordered regions make only unfavorable repulsive contributions to the total energy. We find that the second approach has greater practical utility and illustrate this with examples from de novo structure prediction, NMR structure calculation, and comparative modeling.     

Computation of a stabilizing set of feedback matrices of a large-scale nonlinear musculoskeletal dynamic model

The purpose of this study is to present a general mathematical framework to compute a set of feedback matrices which stabilize an unstable nonlinear anthropomorphic musculoskeletal dynamic model. This method is activity specific and involves four fundamental stages. First, from muscle activation data (input) and motion degrees-of-freedom (output) a dynamic experimental model is obtained using system identification schemes. Second, a nonlinear musculoskeletal dynamic model which contains the same number of muscles and degrees-of-freedom and best represents the activity being considered is proposed. Third, the nonlinear musculoskeletal model (anthropomorphic model) is replaced by a family of linear systems, parameterized by the same set of input/output data (nominal points) used in the identification of the experimental model. Finally, a set of stabilizing output feedback matrices, parameterized again by the same set of nominal points, is computed such that when combined with the anthropomorphic model, the combined system resembles the structural form of the experimental model. The method is illustrated in regard to the human squat activity.     

Computation of a Stabilizing Set of Feedback Matrices bf a Large-scale Nonlinear Musculoskeletal Dynamic Model

The purpose of this study is to present a general mathematical framework to compute a set of feedback matrices which stabilize an unstable nonlinear anthropomorphic musculoskeletal dynamic model. This method is activity specific and involves four fundamental stages. First, from muscle activation data (input) and motion degrees-of-freedom (output) a dynamic experimental model is obtained using system identification schemes. Second, a nonlinear musculoskeletal dynamic model which contains the same number of muscles and degrees-of-freedom and best represents the activity being considered is proposed. Third, the nonlinear musculoskeletal model (anthropomorphic model) is replaced by a family of linear systems, parameterized by the same set of input/ output data (nominal points) used in the identification of the experimental model. Finally, a set of stabilizing output feedback matrices, parameterized again by the same set of nominal points, is computed such that when combined with the anthropomorphic model, the combined system resembles the structural form of the experimental model. The method is illustrated in regard to the human squat activity.     

Identification of the Pain Process by Cold Stimulation: Using Dynamic Causal Modeling of Effective Connectivity in Functional Near-Infrared Spectroscopy (fNIRS)

Background

Pain is an unpleasant sensory and emotional experience followed by anxiety, depression, and frustration. Functional Near-Infrared Spectroscopy (fNIRS) as an optical technique identifies the brain functional networks by investigating connectivity between functionally linked of different anatomical regions in response to pain stimulation.

Dynamic causal models and physiological inference: a validation study using isoflurane anaesthesia in rodents

Generative models of neuroimaging and electrophysiological data present new opportunities for accessing hidden or latent brain states. Dynamic causal modeling (DCM) uses Bayesian model inversion and selection to infer the synaptic mechanisms underlying empirically observed brain responses. DCM for electrophysiological data, in particular, aims to estimate the relative strength of synaptic transmission at different cell types and via specific neurotransmitters. Here, we report a DCM validation study concerning inference on excitatory and inhibitory synaptic transmission, using different doses of a volatile anaesthetic agent (isoflurane) to parametrically modify excitatory and inhibitory synaptic processing while recording local field potentials (LFPs) from primary auditory cortex (A1) and the posterior auditory field (PAF) in the auditory belt region in rodents. We test whether DCM can infer, from the LFP measurements, the expected drug-induced changes in synaptic transmission mediated via fast ionotropic receptors; i.e., excitatory (glutamatergic) AMPA and inhibitory GABA(A) receptors. Cross- and auto-spectra from the two regions were used to optimise three DCMs based on biologically plausible neural mass models and specific network architectures. Consistent with known extrinsic connectivity patterns in sensory hierarchies, we found that a model comprising forward connections from A1 to PAF and backward connections from PAF to A1 outperformed a model with forward connections from PAF to A1 and backward connections from A1 to PAF and a model with reciprocal lateral connections. The parameter estimates from the most plausible model indicated that the amplitude of fast glutamatergic excitatory postsynaptic potentials (EPSPs) and inhibitory postsynaptic potentials (IPSPs) behaved as predicted by previous neurophysiological studies. Specifically, with increasing levels of anaesthesia, glutamatergic EPSPs decreased linearly, whereas fast GABAergic IPSPs displayed a nonlinear (saturating) increase. The consistency of our model-based in vivo results with experimental in vitro results lends further validity to the capacity of DCM to infer on synaptic processes using macroscopic neurophysiological data.     

An autonomous system for identifying and governing a cell's state in yeast

We present an approach for an autonomous system that detects a particular state of interest in a living cell and can govern cell fate accordingly. Cell states could be better identified by the expression pattern of several genes than of a single one. Therefore, autonomous identification can be achieved by a system that measures the expression of these several genes and integrates their activities into a single output. We have constructed a system that diagnoses a unique state in yeast, in which two independent pathways, methionine anabolism and galactose catabolism, are active. Our design is based on modifications of the yeast two-hybrid system. We show that cells could autonomously report on their state, identify the state of interest, and inhibit their growth accordingly. The system's sensitivity is adjustable to detect states with limited dynamic range of inputs. The system's output depends only on the activity of input pathways, not on their identity; hence it is straightforward to diagnose any pair of inputs. A simple model is presented that accounts for the data and provides predictive power. We propose that such systems could handle real-life states-of-interest such as identification of aberrant versus normal growth.     

An Entangled Model for Sustainability Indicators

Nowadays the challenge for humanity is to find pathways towards sustainable development. Decision makers require a set of sustainability indicators to know if the sustainability strategies are following those pathways. There are more than one hundred sustainability indicators but they differ on their relative importance according to the size of the locality and change on time. The resources needed to follow these sustainability indicators are scarce and in some instances finite, especially in smaller regions. Therefore strategies to select set of these indicators are useful for decision makers responsible for monitoring sustainability. In this paper we propose a model for the identification and selection of a set of sustainability indicators that adequately represents human systems. In developing this model, we applied evolutionary dynamics in a space where sustainability indicators are fundamental entities interconnected by an interaction matrix. we used a fixed interaction that simulates the current context for the city of Cuernavaca, México as an example. We were able to identify and define relevant sets indicators for the system by using the Pareto principle. In this case we identified a set of sixteen sustainability indicators with more than 80% of the total strength. This set presents resilience to perturbations. For the Tangled Nature framework we provided a manner of treating different contexts (i.e., cities, counties, states, regions, countries, continents or the whole planet), dealing with small dimensions. This model provides decision makers with a valuable tool to select sustainability indicators set for towns, cities, regions, countries, continents or the entire planet according to a coevolutionary framework. The social legitimacy can arise from the fact that each individual indicator must be selected from those that are most important for the subject community.     

Brain Network Adaptability across Task States

Activity in the human brain moves between diverse functional states to meet the demands of our dynamic environment, but fundamental principles guiding these transitions remain poorly understood. Here, we capitalize on recent advances in network science to analyze patterns of functional interactions between brain regions. We use dynamic network representations to probe the landscape of brain reconfigurations that accompany task performance both within and between four cognitive states: a task-free resting state, an attention-demanding state, and two memory-demanding states. Using the formalism of hypergraphs, we identify the presence of groups of functional interactions that fluctuate coherently in strength over time both within (task-specific) and across (task-general) brain states. In contrast to prior emphases on the complexity of many dyadic (region-to-region) relationships, these results demonstrate that brain adaptability can be described by common processes that drive the dynamic integration of cognitive systems. Moreover, our results establish the hypergraph as an effective measure for understanding functional brain dynamics, which may also prove useful in examining cross-task, cross-age, and cross-cohort functional change.     

Distributed parameters deterministic model for treatment of brain tumors using Galerkin finite element method

In this paper, we present a distributed parameters deterministic model for treatment of brain tumors using Galerkin finite element method. The dynamic model comprises system of three coupled reaction-diffusion models, involving the tumor cells, the normal tissues and the drug concentration. An optimal control problem is formulated with the goal of minimizing the tumor cell density and reducing the side effects of the drug. A distributed parameters method based on the application of variational calculus is used on an integral-Hamiltonian, which is then used to obtain an optimal coupled system of forward state equations and backward co-state equations. The Galerkin finite element method is used to realistically represent the brain structure as well as to facilitate computation. Finally a three-dimensional test case is considered and partitioned into a set of spherical finite elements, using tri-linear basis functions, except for the elements affected by singularities of polar and azimuthal angles, as well as the origin.     

Methods to study protein dynamics and folding by mass spectrometry

It is clear that merely knowing the structure of a protein alone is not sufficient to fully understand its behavior: knowledge also of the dynamic events that occur within proteins is vital to elucidate their function and folding. In recent years, mass spectrometry has come to the forefront as a powerful biophysical method, which can shed light both on the structure and dynamics of proteins. Hydrogen exchange monitored by mass spectrometry is a highly sensitive marker of the backbone dynamics in solution that, combined with gas phase dissociation techniques, can provide a high resolution tool to locate the dynamic regions of a protein. Additionally, charge state distributions in electrospray mass spectra yield insight into the nature and population of alternate structural states present at equilibrium. In this paper, we describe several applications of these methodologies to probe the dynamic events key to the structure, folding, and biological functions of proteins.     

A mosaic pattern characterizes the evolution of the avian brain

Diversity in vertebrate brain size and composition is thought to arise from either developmental constraints that cause coordinated changes between brain regions or a mosaic model, whereby changes in individual brain regions are independent of changes in other brain regions. These two mechanisms were tested in birds using multiple regression analyses. Across 13 orders, significant correlations were present between some brain regions, but not all. Most of the correlated changes reflect the connectivity between different brain components, such that regions with the most interconnections are correlated with one another but not other brain regions. Whether mosaic changes are characteristic of brain regions or systems in birds, however, to our knowledge, remains to be investigated.     

Dynamic causal models of neural system dynamics: current state and future extensions

Complex processes resulting from interaction of multiple elements can rarely be understood by analytical scientific approaches alone; additional, mathematical models of system dynamics are required. This insight, which disciplines like physics have embraced for a long time already, is gradually gaining importance in the study of cognitive processes by functional neuroimaging. In this field, causal mechanisms in neural systems are described in terms of effective connectivity. Recently, dynamic causal modelling (DCM) was introduced as a generic method to estimate effective connectivity from neuroimaging data in a Bayesian fashion. One of the key advantages of DCM over previous methods is that it distinguishes between neural state equations and modality-specific forward models that translate neural activity into a measured signal. Another strength is its natural relation to Bayesian model selection (BMS) procedures. In this article, we review the conceptual and mathematical basis of DCM and its implementation for functional magnetic resonance imaging data and event-related potentials. After introducing the application of BMS in the context of DCM, we conclude with an outlook to future extensions of DCM. These extensions are guided by the long-term goal of using dynamic system models for pharmacological and clinical applications, particularly with regard to synaptic plasticity.     

A quantitative model of sleep-wake dynamics based on the physiology of the brainstem ascending arousal system

A quantitative, physiology-based model of the ascending arousal system is developed, using continuum neuronal population modeling, which involves averaging properties such as firing rates across neurons in each population. The model includes the ventrolateral preoptic area (VLPO), where circadian and homeostatic drives enter the system, the monoaminergic and cholinergic nuclei of the ascending arousal system, and their interconnections. The human sleep-wake cycle is governed by the activities of these nuclei, which modulate the behavioral state of the brain via diffuse neuromodulatory projections. The model parameters are not free since they correspond to physiological observables. Approximate parameter bounds are obtained by requiring consistency with physiological and behavioral measures, and the model replicates the human sleep-wake cycle, with physiologically reasonable voltages and firing rates. Mutual inhibition between the wake-promoting monoaminergic group and sleep-promoting VLPO causes ;;flip-flop' behavior, with most time spent in 2 stable steady states corresponding to wake and sleep, with transitions between them on a timescale of a few minutes. The model predicts hysteresis in the sleep-wake cycle, with a region of bistability of the wake and sleep states. Reducing the monoaminergic-VLPO mutual inhibition results in a smaller hysteresis loop. This makes the model more prone to wake-sleep transitions in both directions and makes the states less distinguishable, as in narcolepsy. The model behavior is robust across the constrained parameter ranges, but with sufficient flexibility to describe a wide range of observed phenomena.