Efficient Parallel Statistical Model Checking of Biochemical Networks

We consider the problem of verifying stochastic models of biochemical networks against behavioral properties expressed in temporal logic terms. Exact probabilistic verification approaches such as, for example, CSL/PCTL model checking, are undermined by a huge computational demand which rule them out for most real case studies. Less demanding approaches, such as statistical model checking, estimate the likelihood that a property is satisfied by sampling executions out of the stochastic model. We propose a methodology for efficiently estimating the likelihood that a LTL property P holds of a stochastic model of a biochemical network. As with other statistical verification techniques, the methodology we propose uses a stochastic simulation algorithm for generating execution samples, however there are three key aspects that improve the efficiency: first, the sample generation is driven by on-the-fly verification of P which results in optimal overall simulation time. Second, the confidence interval estimation for the probability of P to hold is based on an efficient variant of the Wilson method which ensures a faster convergence. Third, the whole methodology is designed according to a parallel fashion and a prototype software tool has been implemented that performs the sampling/verification process in parallel over an HPC architecture

Analyzing various models of Circadian Clock and Cell Cycle coupling

The daily rhythm can influence the proliferation rate of many cell types. In the mammalian system the transcription of the cell cycle regulatory protein Wee1 is controlled by the circadian clock. Zamborszky et al. (2007) present a computational model coupling the cell cycle and circadian rhythm, showing that this coupling can lead to multimodal cell cycle time distributions. Biological data points to additional couplings, including a link back from the cell cycle to the circadian clock. Proper modelling of this coupling requires a more detailed description of both parts of the model. Hence, we aim at further extending and analysing earlier models using a combination of modelling techniques and computer software, including CoSBI lab, BIOCHAM, and GINsim

Fronto-striatal alterations correlate with apathy severity in behavioral variant frontotemporal dementia

Structural and functional changes in cortical and subcortical regions have been reported in behavioral variant frontotemporal dementia (bvFTD), however, a multimodal approach may provide deeper insights into the neural correlates of neuropsychiatric symptoms. In this multicenter study, we measured cortical thickness (CTh) and subcortical volumes to identify structural abnormalities in 37 bvFTD patients, and 37 age- and sex-matched healthy controls. For seed regions with significant structural changes, whole-brain functional connectivity (FC) was examined in a sub-cohort of N = 22 bvFTD and N = 22 matched control subjects to detect complementary alterations in brain network organization. To explore the functional significance of the observed structural and functional deviations, correlations with clinical and neuropsychological outcomes were tested where available. Significantly decreased CTh was observed in the bvFTD group in caudal middle frontal gyrus, left pars opercularis, bilateral superior frontal and bilateral middle temporal gyrus along with subcortical volume reductions in bilateral basal ganglia, thalamus, hippocampus, and amygdala. Resting-state functional magnetic resonance imaging showed decreased FC in bvFTD between: dorsal striatum and left caudal middle frontal gyrus;putamen and fronto-parietal regions;pallidum and cerebellum. Conversely, bvFTD showed increased FC between: left middle temporal gyrus and paracingulate gyrus;caudate nucleus and insula;amygdala and parahippocampal gyrus. Additionally, cortical thickness in caudal, lateral and superior frontal regions as well as caudate nucleus volume correlated negatively with apathy severity scores of the Neuropsychiatry Inventory Questionnaire. In conclusion, multimodal structural and functional imaging indicates that fronto-striatal regions have a considerable influence on the severity of apathy in bvFTD

Snazer: the simulations and networks analyzer

<p>Abstract</p> <p>Background</p> <p>Networks are widely recognized as key determinants of structure and function in systems that span the biological, physical, and social sciences. They are static pictures of the interactions among the components of complex systems. Often, much effort is required to identify networks as part of particular patterns as well as to visualize and interpret them.</p> <p>From a pure dynamical perspective, simulation represents a relevant <it>way</it>-<it>out</it>. Many simulator tools capitalized on the "noisy" behavior of some systems and used formal models to represent cellular activities as temporal trajectories. Statistical methods have been applied to a fairly large number of replicated trajectories in order to infer knowledge.</p> <p>A tool which both graphically manipulates reactive models and deals with sets of simulation time-course data by aggregation, interpretation and statistical analysis is missing and could add value to simulators.</p> <p>Results</p> <p>We designed and implemented <it>Snazer</it>, the simulations and networks analyzer. Its goal is to aid the processes of visualizing and manipulating reactive models, as well as to share and interpret time-course data produced by stochastic simulators or by any other means.</p> <p>Conclusions</p> <p><it>Snazer </it>is a solid prototype that integrates biological network and simulation time-course data analysis techniques.</p

Improving 3D convolutional neural network comprehensibility via interactive visualization of relevance maps: Evaluation in Alzheimer's disease

Background: Although convolutional neural networks (CNN) achieve high diagnostic accuracy for detecting Alzheimer's disease (AD) dementia based on magnetic resonance imaging (MRI) scans, they are not yet applied in clinical routine. One important reason for this is a lack of model comprehensibility. Recently developed visualization methods for deriving CNN relevance maps may help to fill this gap. We investigated whether models with higher accuracy also rely more on discriminative brain regions predefined by prior knowledge. Methods: We trained a CNN for the detection of AD in N=663 T1-weighted MRI scans of patients with dementia and amnestic mild cognitive impairment (MCI) and verified the accuracy of the models via cross-validation and in three independent samples including N=1655 cases. We evaluated the association of relevance scores and hippocampus volume to validate the clinical utility of this approach. To improve model comprehensibility, we implemented an interactive visualization of 3D CNN relevance maps. Results: Across three independent datasets, group separation showed high accuracy for AD dementia vs. controls (AUC$\geq$0.92) and moderate accuracy for MCI vs. controls (AUC$\approx$0.75). Relevance maps indicated that hippocampal atrophy was considered as the most informative factor for AD detection, with additional contributions from atrophy in other cortical and subcortical regions. Relevance scores within the hippocampus were highly correlated with hippocampal volumes (Pearson's r$\approx$-0.86, p<0.001). Conclusion: The relevance maps highlighted atrophy in regions that we had hypothesized a priori. This strengthens the comprehensibility of the CNN models, which were trained in a purely data-driven manner based on the scans and diagnosis labels.Comment: 24 pages, 9 figures/tables, supplementary material, source code available on GitHu

An analysis of irreversible transitions in a model of the buddying yeast cell cycle

Cells life follows a cycling behaviour which starts at cell birth and leads to cell division through a number of distinct phases. The transitions through the various cell cycle phases are controlled by a complex network of signalling pathways. Many cell cycle transitions are irreversible: once they are started they must reach completion. In this study we investigate the existence of conditions which lead to cases when irreversibility may be broken. Specifically, we characterise the elements of the cell cycle signalling network that are responsible for the irreversibility and we determine conditions for which the irreversible transitions may become reversible. We illustrate our results through a formal approach in which stochastic simulation analysis and model checking verification are combined. Through probabilistic model checking we provide a quantitative measure for the probability of irreversibility in the ``Start" transition of the cell cycle. This is the preliminary version of a paper that was published in Proceedings of the third International Workshop on Practical Applications of Stochastic Modelling (PASM 2008), ENTCS 232, pp. 39-53, 2009. (DOI:10.1016/j.entcs.2009.02.049