Conditions for duality between fluxes and concentrations in biochemical networks

Mathematical and computational modelling of biochemical networks is often done in terms of either the concentrations of molecular species or the fluxes of biochemical reactions. When is mathematical modelling from either perspective equivalent to the other? Mathematical duality translates concepts, theorems or mathematical structures into other concepts, theorems or structures, in a one-to-one manner. We present a novel stoichiometric condition that is necessary and sufficient for duality between unidirectional fluxes and concentrations. Our numerical experiments, with computational models derived from a range of genome-scale biochemical networks, suggest that this flux-concentration duality is a pervasive property of biochemical networks. We also provide a combinatorial characterisation that is sufficient to ensure flux-concentration duality. That is, for every two disjoint sets of molecular species, there is at least one reaction complex that involves species from only one of the two sets. When unidirectional fluxes and molecular species concentrations are dual vectors, this implies that the behaviour of the corresponding biochemical network can be described entirely in terms of either concentrations or unidirectional fluxes

Optimal and Approximate Q-value Functions for Decentralized POMDPs

Decision-theoretic planning is a popular approach to sequential decision making problems, because it treats uncertainty in sensing and acting in a principled way. In single-agent frameworks like MDPs and POMDPs, planning can be carried out by resorting to Q-value functions: an optimal Q-value function Q* is computed in a recursive manner by dynamic programming, and then an optimal policy is extracted from Q*. In this paper we study whether similar Q-value functions can be defined for decentralized POMDP models (Dec-POMDPs), and how policies can be extracted from such value functions. We define two forms of the optimal Q-value function for Dec-POMDPs: one that gives a normative description as the Q-value function of an optimal pure joint policy and another one that is sequentially rational and thus gives a recipe for computation. This computation, however, is infeasible for all but the smallest problems. Therefore, we analyze various approximate Q-value functions that allow for efficient computation. We describe how they relate, and we prove that they all provide an upper bound to the optimal Q-value function Q*. Finally, unifying some previous approaches for solving Dec-POMDPs, we describe a family of algorithms for extracting policies from such Q-value functions, and perform an experimental evaluation on existing test problems, including a new firefighting benchmark problem

Assessment of progressive collapse in multi-storey buildings

Accepted versio

Troppo - A Python framework for the reconstruction of context-specific metabolic models

The surge in high-throughput technology availability for molecular biology has enabled the development of powerful predictive tools for use in many applications, including (but not limited to) the diagnosis and treatment of human diseases such as cancer. Genome-scale metabolic models have shown some promise in clearing a path towards precise and personalized medicine, although some challenges still persist. The integration of omics data and subsequent creation of context-specific models for specific cells/tissues still poses a significant hurdle, and most current tools for this purpose have been implemented using proprietary software. Here, we present a new software tool developed in Python, troppo - Tissue-specific RecOnstruction and Phenotype Prediction using Omics data, implementing a large variety of context-specific reconstruction algorithms. Our framework and workflow are modular, which facilitates the development of newer algorithms or omics data sources.This study was supported by the Portuguese Foundation for Science and Technology (FCT) under the scope of the strategic funding of UID/BIO/04469/2019 unit and BioTecNorte operation (NORTE-01-0145-FEDER-000004) funded by the European Regional Development Fund under the scope of Norte2020 - Programa Operacional Regional do Norte. The authors also thank the PhD scholarships funded by national funds through Fundacao para a Ciencia e Tecnologia, with references: SFRH/BD/133248/2017 (J.F.), SFRH/BD/118657/2016 (V.V.).info:eu-repo/semantics/publishedVersio

Gaussian Mixture Model of Heart Rate Variability

Heart rate variability (HRV) is an important measure of sympathetic and parasympathetic functions of the autonomic nervous system and a key indicator of cardiovascular condition. This paper proposes a novel method to investigate HRV, namely by modelling it as a linear combination of Gaussians. Results show that three Gaussians are enough to describe the stationary statistics of heart variability and to provide a straightforward interpretation of the HRV power spectrum. Comparisons have been made also with synthetic data generated from different physiologically based models showing the plausibility of the Gaussian mixture parameters

Evaluation of global simulations of aerosol particle and cloud condensation nuclei number, with implications for cloud droplet formation

A total of 16 global chemistry transport models and general circulation models have participated in this study; 14 models have been evaluated with regard to their ability to reproduce the near-surface observed number concentration of aerosol particles and cloud condensation nuclei (CCN), as well as derived cloud droplet number concentration (CDNC). Model results for the period 2011-2015 are compared with aerosol measurements (aerosol particle number, CCN and aerosol particle composition in the submicron fraction) from nine surface stations located in Europe and Japan. The evaluation focuses on the ability of models to simulate the average across time state in diverse environments and on the seasonal and short-term variability in the aerosol properties. There is no single model that systematically performs best across all environments represented by the observations. Models tend to underestimate the observed aerosol particle and CCN number concentrations, with average normalized mean bias (NMB) of all models and for all stations, where data are available, of -24 % and -35 % for particles with dry diameters > 50 and > 120 nm, as well as -36 % and -34 % for CCN at supersaturations of 0.2 % and 1.0 %, respectively. However, they seem to behave differently for particles activating at very low supersaturations (<0.1 %) than at higher ones. A total of 15 models have been used to produce ensemble annual median distributions of relevant parameters. The model diversity (defined as the ratio of standard deviation to mean) is up to about 3 for simulated N-3 (number concentration of particles with dry diameters larger than 3 nm) and up to about 1 for simulated CCN in the extra-polar regions. A global mean reduction of a factor of about 2 is found in the model diversity for CCN at a supersaturation of 0.2 % (CCN0.2) compared to that for N-3, maximizing over regions where new particle formation is important. An additional model has been used to investigate potential causes of model diversity in CCN and bias compared to the observations by performing a perturbed parameter ensemble (PPE) accounting for uncertainties in 26 aerosol-related model input parameters. This PPE suggests that biogenic secondary organic aerosol formation and the hygroscopic properties of the organic material are likely to be the major sources of CCN uncertainty in summer, with dry deposition and cloud processing being dominant in winter. Models capture the relative amplitude of the seasonal variability of the aerosol particle number concentration for all studied particle sizes with available observations (dry diameters larger than 50, 80 and 120 nm). The short-term persistence time (on the order of a few days) of CCN concentrations, which is a measure of aerosol dynamic behavior in the models, is underestimated on average by the models by 40 % during winter and 20 % in summer. In contrast to the large spread in simulated aerosol particle and CCN number concentrations, the CDNC derived from simulated CCN spectra is less diverse and in better agreement with CDNC estimates consistently derived from the observations (average NMB -13 % and -22 % for updraft velocities 0.3 and 0.6 m s(-1), respectively). In addition, simulated CDNC is in slightly better agreement with observationally derived values at lower than at higher updraft velocities (index of agreement 0.64 vs. 0.65). The reduced spread of CDNC compared to that of CCN is attributed to the sublinear response of CDNC to aerosol particle number variations and the negative correlation between the sensitivities of CDNC to aerosol particle number concentration (partial derivative N-d/partial derivative N-a) and to updraft velocity (partial derivative N-d/partial derivative w). Overall, we find that while CCN is controlled by both aerosol particle number and composition, CDNC is sensitive to CCN at low and moderate CCN concentrations and to the updraft velocity when CCN levels are high. Discrepancies are found in sensitivities partial derivative N-d/partial derivative N-a and partial derivative N-d/partial derivative w; models may be predisposed to be too "aerosol sensitive" or "aerosol insensitive" in aerosol-cloud-climate interaction studies, even if they may capture average droplet numbers well. This is a subtle but profound finding that only the sensitivities can clearly reveal and may explain intermodel biases on the aerosol indirect effect.Peer reviewe

Evaluation of Global Simulations of Aerosol Particle and Cloud Condensation Nuclei Number, with Implications for Cloud Droplet Formation

A total of 16 global chemistry transport models and general circulation models have participated in this study; 14 models have been evaluated with regard to their ability to reproduce the near-surface observed number concentration of aerosol particles and cloud condensation nuclei (CCN), as well as derived cloud droplet number concentration (CDNC). Model results for the period 2011-2015 are compared with aerosol measurements (aerosol particle number, CCN and aerosol particle composition in the submicron fraction) from nine surface stations located in Europe and Japan. The evaluation focuses on the ability of models to simulate the average across time state in diverse environments and on the seasonal and short-term variability in the aerosol properties. There is no single model that systematically performs best across all environments represented by the observations. Models tend to underestimate the observed aerosol particle and CCN number concentrations, with average normalized mean bias (NMB) of all models and for all stations, where data are available, of -24% and -35% for particles with dry diameters > 50 and > 120nm, as well as -36% and -34% for CCN at supersaturations of 0.2% and 1.0%, respectively. However, they seem to behave differently for particles activating at very low supersaturations (< 0.1%) than at higher ones. A total of 15 models have been used to produce ensemble annual median distributions of relevant parameters. The model diversity (defined as the ratio of standard deviation to mean) is up to about 3 for simulated N3 (number concentration of particles with dry diameters larger than 3 nm) and up to about 1 for simulated CCN in the extra-polar regions. A global mean reduction of a factor of about 2 is found in the model diversity for CCN at a supersaturation of 0.2% (CCN(0.2)) compared to that for N3, maximizing over regions where new particle formation is important. An additional model has been used to investigate potential causes of model diversity in CCN and bias compared to the observations by performing a perturbed parameter ensemble (PPE) accounting for uncertainties in 26 aerosol-related model input parameters. This PPE suggests that biogenic secondary organic aerosol formation and the hygroscopic properties of the organic material are likely to be the major sources of CCN uncertainty in summer, with dry deposition and cloud processing being dominant in winter. Models capture the relative amplitude of the seasonal variability of the aerosol particle number concentration for all studied particle sizes with available observations (dry diameters larger than 50, 80 and 120nm). The short-term persistence time (on the order of a few days) of CCN concentrations, which is a measure of aerosol dynamic behavior in the models, is underestimated on average by the models by 40% during winter and 20% in summer