Markov chain aggregation and its application to rule-based modelling

Rule-based modelling allows to represent molecular interactions in a compact and natural way. The underlying molecular dynamics, by the laws of stochastic chemical kinetics, behaves as a continuous-time Markov chain. However, this Markov chain enumerates all possible reaction mixtures, rendering the analysis of the chain computationally demanding and often prohibitive in practice. We here describe how it is possible to efficiently find a smaller, aggregate chain, which preserves certain properties of the original one. Formal methods and lumpability notions are used to define algorithms for automated and efficient construction of such smaller chains (without ever constructing the original ones). We here illustrate the method on an example and we discuss the applicability of the method in the context of modelling large signalling pathways

A novel brain receptor is expressed in a distinct population of olfactory sensory neurons

Three novel G-protein-coupled receptor genes related to the previously described RA1c gene have been isolated from the mouse genome. Expression of these genes has been detected in distinct areas of the brain and also in the olfactory epithelium of the nose. Developmental studies revealed a differential onset of expression: in the brain at embryonic stage 17, in the olfactory system at stage E12. In order to determine which cell type in the olfactory epithelium expresses this unique receptor type, a transgenic approach was employed which allowed a coexpression of histological markers together with the receptor and thus visualization of the appropriate cell population. It was found that the receptor-expressing cells were located very close to the basal membrane of the epithelium; however, the cells extended a dendritic process to the epithelial surface and their axons projected into the main olfactory bulb where they converged onto two or three glomeruli in the dorsal and posterior region of the bulb. Thus, these data provide evidence that this unique type of receptor is expressed in mature olfactory neurons and suggests that it may be involved in the detection of special odour molecules

Is fall prevention by vitamin D mediated by a change in postural or dynamic balance?

Introduction: The objectives were:(1) to validate a quantitative balance assessment method for fall risk prediction; (2) to investigate whether the effect of vitamin D and calcium on the risk of falling is mediated through postural or dynamic balance, as assessed by this method. Materials and methods: A secondary analysis of a double blind randomized controlled trial was employed, which included 64 institutionalized elderly women with complete balance assessment (age range: 65-97; mean 25-hydroxyvitamin D levels: 16.4ng/ml (SD ±9.9). Participants received 1,200mg calcium plus 800IU cholecalciferol (n=33) or 1,200mg calcium (n=31) per day over a 3-month treatment period. Using an electronic device attached to the lower back of the participant, balance was assessed as the degree of trunk angular displacement and angular velocity during a postural task (standing on two legs, eyes open, for 20 s) and a dynamic task (get up from a standard height chair with arm rests, sit down and then stand up again and remain standing). Results: It was found that both postural and dynamic balance independently and significantly predicted the rate of falling within the 3-month follow-up. Vitamin D plus calcium reduced the rate of falls by 60% [relative risk (RR)=0.40; 95% CI: 0.17, 0.94] if compared with calcium alone. Once postural and dynamic balance were added to the regression analysis, they both attenuated the effect of vitamin D plus calcium on the rate of falls. For postural balance, the RR changed by 22% from 0.40 to 0.62 if angular displacement was added to the model, and by 9% from 0.40 to 0.49 if angular velocity was added. For dynamic balance, it changed by 1% from 0.40 to 0.41 if angular displacement was added, and by 14% from 0.40 to 0.54 if angular velocity was added. Discussion: Thus, balance assessment using trunk angular displacement is a valid method for the prediction of falls in older women. Of the observed 60% reduction in the rate of falls by vitamin D plus calcium supplementation compared with calcium alone, up to 22% of the treatment effect was explained by a change in postural balance and up to 14% by dynamic balanc

Syntactic Markovian Bisimulation for Chemical Reaction Networks

In chemical reaction networks (CRNs) with stochastic semantics based on continuous-time Markov chains (CTMCs), the typically large populations of species cause combinatorially large state spaces. This makes the analysis very difficult in practice and represents the major bottleneck for the applicability of minimization techniques based, for instance, on lumpability. In this paper we present syntactic Markovian bisimulation (SMB), a notion of bisimulation developed in the Larsen-Skou style of probabilistic bisimulation, defined over the structure of a CRN rather than over its underlying CTMC. SMB identifies a lumpable partition of the CTMC state space a priori, in the sense that it is an equivalence relation over species implying that two CTMC states are lumpable when they are invariant with respect to the total population of species within the same equivalence class. We develop an efficient partition-refinement algorithm which computes the largest SMB of a CRN in polynomial time in the number of species and reactions. We also provide an algorithm for obtaining a quotient network from an SMB that induces the lumped CTMC directly, thus avoiding the generation of the state space of the original CRN altogether. In practice, we show that SMB allows significant reductions in a number of models from the literature. Finally, we study SMB with respect to the deterministic semantics of CRNs based on ordinary differential equations (ODEs), where each equation gives the time-course evolution of the concentration of a species. SMB implies forward CRN bisimulation, a recently developed behavioral notion of equivalence for the ODE semantics, in an analogous sense: it yields a smaller ODE system that keeps track of the sums of the solutions for equivalent species.Comment: Extended version (with proofs), of the corresponding paper published at KimFest 2017 (http://kimfest.cs.aau.dk/

An overview of the mid-infrared spectro-interferometer MATISSE: science, concept, and current status

MATISSE is the second-generation mid-infrared spectrograph and imager for the Very Large Telescope Interferometer (VLTI) at Paranal. This new interferometric instrument will allow significant advances by opening new avenues in various fundamental research fields: studying the planet-forming region of disks around young stellar objects, understanding the surface structures and mass loss phenomena affecting evolved stars, and probing the environments of black holes in active galactic nuclei. As a first breakthrough, MATISSE will enlarge the spectral domain of current optical interferometers by offering the L and M bands in addition to the N band. This will open a wide wavelength domain, ranging from 2.8 to 13 um, exploring angular scales as small as 3 mas (L band) / 10 mas (N band). As a second breakthrough, MATISSE will allow mid-infrared imaging - closure-phase aperture-synthesis imaging - with up to four Unit Telescopes (UT) or Auxiliary Telescopes (AT) of the VLTI. Moreover, MATISSE will offer a spectral resolution range from R ~ 30 to R ~ 5000. Here, we present one of the main science objectives, the study of protoplanetary disks, that has driven the instrument design and motivated several VLTI upgrades (GRA4MAT and NAOMI). We introduce the physical concept of MATISSE including a description of the signal on the detectors and an evaluation of the expected performances. We also discuss the current status of the MATISSE instrument, which is entering its testing phase, and the foreseen schedule for the next two years that will lead to the first light at Paranal.Comment: SPIE Astronomical Telescopes and Instrumentation conference, June 2016, 11 pages, 6 Figure

Computational exploration of molecular receptive fields in the olfactory bulb reveals a glomerulus-centric chemical map

© The Author(s) 2020. This article is licensed under a Creative Commons Attribution 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.Progress in olfactory research is currently hampered by incomplete knowledge about chemical receptive ranges of primary receptors. Moreover, the chemical logic underlying the arrangement of computational units in the olfactory bulb has still not been resolved. We undertook a large-scale approach at characterising molecular receptive ranges (MRRs) of glomeruli in the dorsal olfactory bulb (dOB) innervated by the MOR18-2 olfactory receptor, also known as Olfr78, with human ortholog OR51E2. Guided by an iterative approach that combined biological screening and machine learning, we selected 214 odorants to characterise the response of MOR18-2 and its neighbouring glomeruli. We found that a combination of conventional physico-chemical and vibrational molecular descriptors performed best in predicting glomerular responses using nonlinear Support-Vector Regression. We also discovered several previously unknown odorants activating MOR18-2 glomeruli, and obtained detailed MRRs of MOR18-2 glomeruli and their neighbours. Our results confirm earlier findings that demonstrated tunotopy, that is, glomeruli with similar tuning curves tend to be located in spatial proximity in the dOB. In addition, our results indicate chemotopy, that is, a preference for glomeruli with similar physico-chemical MRR descriptions being located in spatial proximity. Together, these findings suggest the existence of a partial chemical map underlying glomerular arrangement in the dOB. Our methodology that combines machine learning and physiological measurements lights the way towards future high-throughput studies to deorphanise and characterise structure-activity relationships in olfaction.Peer reviewe

The Process-Interaction-Model: a common representation of rule-based and logical models allows studying signal transduction on different levels of detail

BACKGROUND: Signaling systems typically involve large, structured molecules each consisting of a large number of subunits called molecule domains. In modeling such systems these domains can be considered as the main players. In order to handle the resulting combinatorial complexity, rule-based modeling has been established as the tool of choice. In contrast to the detailed quantitative rule-based modeling, qualitative modeling approaches like logical modeling rely solely on the network structure and are particularly useful for analyzing structural and functional properties of signaling systems. RESULTS: We introduce the Process-Interaction-Model (PIM) concept. It defines a common representation (or basis) of rule-based models and site-specific logical models, and, furthermore, includes methods to derive models of both types from a given PIM. A PIM is based on directed graphs with nodes representing processes like post-translational modifications or binding processes and edges representing the interactions among processes. The applicability of the concept has been demonstrated by applying it to a model describing EGF insulin crosstalk. A prototypic implementation of the PIM concept has been integrated in the modeling software ProMoT. CONCLUSIONS: The PIM concept provides a common basis for two modeling formalisms tailored to the study of signaling systems: a quantitative (rule-based) and a qualitative (logical) modeling formalism. Every PIM is a compact specification of a rule-based model and facilitates the systematic set-up of a rule-based model, while at the same time facilitating the automatic generation of a site-specific logical model. Consequently, modifications can be made on the underlying basis and then be propagated into the different model specifications – ensuring consistency of all models, regardless of the modeling formalism. This facilitates the analysis of a system on different levels of detail as it guarantees the application of established simulation and analysis methods to consistent descriptions (rule-based and logical) of a particular signaling system