Optimal Bayesian experimental design for models with intractable likelihoods using indirect inference applied to biological process models

This paper addresses the problem of determining optimal designs for biological process models with intractable likelihoods, with the goal of parameter inference. The Bayesian approach is to choose a design that maximises the mean of a utility, and the utility is a function of the posterior distribution. Therefore, its estimation requires likelihood evaluations. However, many problems in experimental design involve models with intractable likelihoods, that is, likelihoods that are neither analytic nor can be computed in a reasonable amount of time. We propose a novel solution using indirect inference (II), a well established method in the literature, and the Markov chain Monte Carlo (MCMC) algorithm of Müller et al. (2004). Indirect inference employs an auxiliary model with a tractable likelihood in conjunction with the generative model, the assumed true model of interest, which has an intractable likelihood. Our approach is to estimate a map between the parameters of the generative and auxiliary models, using simulations from the generative model. An II posterior distribution is formed to expedite utility estimation. We also present a modification to the utility that allows the Müller algorithm to sample from a substantially sharpened utility surface, with little computational effort. Unlike competing methods, the II approach can handle complex design problems for models with intractable likelihoods on a continuous design space, with possible extension to many observations. The methodology is demonstrated using two stochastic models; a simple tractable death process used to validate the approach, and a motivating stochastic model for the population evolution of macroparasites.

Towards the prediction of essential genes by integration of network topology, cellular localization and biological process information

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Polishing of water treated by a biological process with heterogeneous photocatalysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Moving beyond a descriptive aquatic toxicology: the value of biological process and trait information

In order to improve the ability to link chemical exposure to toxicological and ecological effects, aquatic toxicology will have to move from observing what chemical concentrations induce adverse effects to more explanatory approaches, that are concepts which build on knowledge of biological processes and pathways leading from exposure to adverse effects, as well as on knowledge on stressor vulnerability as given by the genetic, physiological and ecological (e.g., life history) traits of biota. Developing aquatic toxicology in this direction faces a number of challenges, including (i) taking into account species differences in toxicant responses on the basis of the evolutionarily developed diversity of phenotypic vulnerability to environmental stressors, (ii) utilizing diversified biological response profiles to serve as biological read across for prioritizing chemicals, categorizing them according to modes of action, and for guiding targeted toxicity evaluation; (iii) prediction of ecological consequences of toxic exposure from knowledge of how biological processes and phenotypic traits lead to effect propagation across the levels of biological hierarchy; and (iv) the search for concepts to assess the cumulative impact of multiple stressors. An underlying theme in these challenges is that, in addition to the question of what the chemical does to the biological receptor, we should give increasing emphasis to the question how the biological receptor handles the chemicals, i.e., through which pathways the initial chemical-biological interaction extends to the adverse effects, how this extension is modulated by adaptive or compensatory processes as well as by phenotypic traits of the biological receptor.

An Architecture forRepresenting Biological Processes based on Networks of Bio-inspired Processors

n this paper we propose the use of Networks of Bio-inspired Processors (NBP) to model some biological phenomena within a computational framework. In particular, we propose the use of an extension of NBP named Network Evolutionary Processors Transducers to simulate chemical transformations of substances. Within a biological process, chemical transformations of substances are basic operations in the change of the state of the cell. Previously, it has been proved that NBP are computationally complete, that is, they are able to solve NP complete problems in linear time, using massively parallel computations. In addition, we propose a multilayer architecture that will allow us to design models of biological processes related to cellular communication as well as their implications in the metabolic pathways. Subsequently, these models can be applied not only to biological-cellular instances but, possibly, also to configure instances of interactive processes in many other fields like population interactions, ecological trophic networks, in dustrial ecosystems, etc.

Acyl glucuronide reactivity in perspective: biological consequences

The metabolic conjugation of exogenous and endogenous carboxylic acid substrates with endogenous glucuronic acid, mediated by the uridine diphosphoglucuronosyl transferase (UGT) superfamily of enzymes, leads to the formation of acyl glucuronide metabolites. Since the late 1970s, acyl glucuronides have been increasingly identified as reactive electrophilic metabolites, capable of undergoing three reactions: intramolecular rearrangement, hydrolysis, and intermolecular reactions with proteins leading to covalent drug-protein adducts. This essential dogma has been accepted for over a decade. The key question proposed by researchers, and now the pharmaceutical industry, is: does or can the covalent modification of endogenous proteins, mediated by reactive acyl glucuronide metabolites, lead to adverse drug reactions, perhaps idiosyncratic in nature? This review evaluates the evidence for acyl glucuronide-derived perturbation of homeostasis, particularly that which might result from the covalent modification of endogenous proteins and other macromolecules. Because of the availability of acyl glucuronides for test tube/in vitro experiments, there is now a substantial literature documenting their rearrangement, hydrolysis and covalent modification of proteins in vitro. It is certain from in vitro experiments that serum albumin, dipeptidyl peptidase IV, tubulin and UGTs are covalently modified by acyl glucuronides. However, these in vitro experiments have been specifically designed to amplify any interference with a biological process in order to find biological effects. The in vivo situation is not at all clear. Certainly it must be concluded that all humans taking carboxylate drugs that form reactive acyl glucuronides will form covalent drug-protein adducts, and it must also be concluded that this in itself is normally benign. However, there is enough in vivo evidence implicating acyl glucuronides, which, when backed up by in vivo circumstantial and documented in vitro evidence, supports the view that reactive acyl glucuronides may initiate toxicity/immune responses. In summary, though acyl glucuronide-derived covalent modification of endogenous macromolecules is well-defined, the work ahead needs to provide detailed links between such modification and its possible biological consequences. (C) 2003 Elsevier Science Ireland Ltd. All rights reserved.

Optimization of integrated chemical-biological degradation of a reactive azo dye using response surface methodology

The integrated chemical-biological degradation combining advanced oxidation by UV/H2O2 followed by aerobic biodegradation was used to degrade C.I. Reactive Azo Red 195A, commonly used in the textile industry in Australia. An experimental design based on the response surface method was applied to evaluate the interactive effects of influencing factors (UV irradiation time, initial hydrogen peroxide dosage and recirculation ratio of the system) on decolourisation efficiency and optimizing the operating conditions of the treatment process. The effects were determined by the measurement of dye concentration and soluble chemical oxygen demand (S-COD). The results showed that the dye and S-COD removal were affected by all factors individually and interactively. Maximal colour degradation performance was predicted, and experimentally validated, with no recirculation, 30 min UV irradiation and 500 mg H2O2/L. The model predictions for colour removal, based on a three-factor/five-level Box-Wilson central composite design and the response surface method analysis, were found to be very close to additional experimental results obtained under near optimal conditions. This demonstrates the benefits of this approach in achieving good predictions while minimising the number of experiments required. (c) 2006 Elsevier B.V. All rights reserved.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

Wastewater treatment by novel hybrid biological – ion exchange processes

This study presents two novel methods for treating important environmental contaminants from two different wastewater streams. One process utilizes the kinetic advantages and reliability of ion exchanging clinoptilolite in combination with biological treatment to remove ammonium from municipal sewage. A second process, HAMBgR (Hybrid Adsorption Membrane Biological Reactor), combines both ion exchange resin and bacteria into a single reactor to treat perchlorate contaminated waters. Combining physicochemical adsorptive treatment with biological treatment can provide synergistic benefits to the overall removal processes. Ion exchange removal solves some of the common operational reliability limitations of biological treatment, like slow response to environmental changes and leaching. Biological activity can in turn help reduce the economic and environmental challenges of ion exchange processes, like regenerant cost and brine disposal. The second section of this study presents continuous flow column experiments, used to demonstrate the ability of clinoptilolite to remove wastewater ammonium, as well as the effectiveness of salt regeneration using highly concentrated sea salt solutions. The working capacity of clinoptilolite more than doubled over the first few loading cycles, while regeneration recovered more than 98% of ammonium. Using the regenerant brine for subsequent halotolerant algae growth allowed for its repeated use, which could lead to cost savings and production of valuable algal biomass. The algae were able to uptake all ammonium in solution, and the brine was able to be used again with no loss in regeneration efficiency. This process has significant advantages over conventional biological nitrification; shorter retention times, wider range of operational conditions, and higher quality effluent free of nitrate. Also, since the clinoptilolite is continually regenerated and the regenerant is rejuvenated by algae, overall input costs are expected to be low. The third section of this study introduces the HAMBgR process for the elimination of perchlorate and presents batch isotherm experiments and pilot reactor tests. Results showed that a variety of ion-exchange resins can be effectively and repeatedly regenerated biologically, and maintain an acceptable working capacity. The presence of an adsorbent in the HAMBgR process improved bioreactor performance during operational fluctuations by providing a physicochemical backup to the biological process. Pilot reactor tests showed that the HAMBgR process reduced effluent perchlorate spikes by up to 97% in comparison to a conventional membrane bio-reactor (MBR) that was subject to sudden changes in influent conditions. Also, the HAMBgR process stimulated biological activity and lead to higher biomass concentrations during increased contaminant loading conditions. Conventional MBR systems can be converted into HAMBgR’s at a low cost, easily justifiable by the realized benefits. The concepts employed in the HAMBgR process can be adapted to treat other target contaminants, not just perchlorate.

Modelling the interaction of keratinocytes and fibroblasts during normal and abnormal wound healing processes

The crosstalk between fibroblasts and keratinocytes is a vital component of the wound healing process, and involves the activity of a number of growth factors and cytokines. In this work, we develop a mathematical model of this crosstalk in order to elucidate the effects of these interactions on the regeneration of collagen in a wound that heals by second intention. We consider the role of four components that strongly affect this process: transforming growth factor-beta, platelet-derived growth factor, interleukin-1 and keratinocyte growth factor. The impact of this network of interactions on the degradation of an initial fibrin clot, as well as its subsequent replacement by a matrix that is mainly comprised of collagen, is described through an eight-component system of nonlinear partial differential equations. Numerical results, obtained in a two-dimensional domain, highlight key aspects of this multifarious process such as reepithelialisation. The model is shown to reproduce many of the important features of normal wound healing. In addition, we use the model to simulate the treatment of two pathological cases: chronic hypoxia, which can lead to chronic wounds; and prolonged inflammation, which has been shown to lead to hypertrophic scarring. We find that our model predictions are qualitatively in agreement with previously reported observations, and provide an alternative pathway for gaining insight into this complex biological process.

An introduction to biologics

An introduction to biologics 23.1 Introduction: Principles of biologics and their use as medicines 23.2 Protein biologics used as drugs 23.2.1 Proteins that function through enzymatic or regulatory activity. 23.2.1.1 Biologics as replacement of a deficient or abnormal protein. 23.2.1.2 Proteins that augment an existing biological process. 23.2.1.3 Proteins that provide a novel function or activity. 23.2.2. Proteins that function through specific targeting activity. 23.2.2.1. Monoclonal antibody nomenclature. 23.2.2.2. Naked monoclonal antibodies. 23.2.2.3. Conjugated monoclonal antibodies. 23.2.3. Recombinant protein vaccines.

Rapid host adaptation by extensive recombination

Experimental investigations into virus recombination can provide valuable insights into the biochemical mechanisms and the evolutionary value of this fundamental biological process. Here, we describe an experimental scheme for studying recombination that should be applicable to any recombinogenic viruses amenable to the production of synthetic infectious genomes. Our approach is based on differences in fitness that generally exist between synthetic chimaeric genomes and the wild-type viruses from which they are constructed. In mixed infections of defective reciprocal chimaeras, selection strongly favours recombinant progeny genomes that recover a portion of wild-type fitness. Characterizing these evolved progeny viruses can highlight both important genetic fitness determinants and the contribution that recombination makes to the evolution of their natural relatives. Moreover, these experiments supply precise information about the frequency and distribution of recombination breakpoints, which can shed light on the mechanistic processes underlying recombination. We demonstrate the value of this approach using the small single-stranded DNA geminivirus, maize streak virus (MSV). Our results show that adaptive recombination in this virus is extremely efficient and can yield complex progeny genomes comprising up to 18 recombination breakpoints. The patterns of recombination that we observe strongly imply that the mechanistic processes underlying rolling circle replication are the prime determinants of recombination breakpoint distributions found in MSV genomes sampled from nature. © 2009 SGM.

A comparison of techniques for the reconstruction of leaf surfaces from scanned data

The foliage of a plant performs vital functions. As such, leaf models are required to be developed for modelling the plant architecture from a set of scattered data captured using a scanning device. The leaf model can be used for purely visual purposes or as part of a further model, such as a fluid movement model or biological process. For these reasons, an accurate mathematical representation of the surface and boundary is required. This paper compares three approaches for fitting a continuously differentiable surface through a set of scanned data points from a leaf surface, with a technique already used for reconstructing leaf surfaces. The techniques which will be considered are discrete smoothing D2-splines [R. Arcangeli, M. C. Lopez de Silanes, and J. J. Torrens, Multidimensional Minimising Splines, Springer, 2004.], the thin plate spline finite element smoother [S. Roberts, M. Hegland, and I. Altas, Approximation of a Thin Plate Spline Smoother using Continuous Piecewise Polynomial Functions, SIAM, 1 (2003), pp. 208--234] and the radial basis function Clough-Tocher method [M. Oqielat, I. Turner, and J. Belward, A hybrid Clough-Tocher method for surface fitting with application to leaf data., Appl. Math. Modelling, 33 (2009), pp. 2582-2595]. Numerical results show that discrete smoothing D2-splines produce reconstructed leaf surfaces which better represent the original physical leaf.

Positional cloning and pathway analysis of the asthma susceptibility gene, NPSR1

In the present study, we identified a novel asthma susceptibility gene, NPSR1 (neuropeptide S receptor 1) on chromosome 7p14.3 by the positional cloning strategy. An earlier significant linkage mapping result among Finnish Kainuu asthma families was confirmed in two independent cohorts: in asthma families from Quebec, Canada and in allergy families from North Karelia, Finland. The linkage region was narrowed down to a 133-kb segment by a hierarchial genotyping method. The observed 77-kb haplotype block showed 7 haplotypes and a similar risk and nonrisk pattern in all three populations studied. All seven haplotypes occur in all three populations at frequences > 2%. Significant elevated relative risks were detected for elevated total IgE (immunoglobulin E) or asthma. Risk effects of the gene variants varied from 1.4 to 2.5. NPSR1 belongs to the G protein-coupled receptor (GPCR) family with a topology of seven transmembrane domains. NPSR1 has 9 exons, with the two main transcripts, A and B, encoding proteins of 371 and 377 amino acids, respectively. We detected a low but ubiquitous expression level of NPSR1-B in various tissues and endogenous cell lines while NPSR1-A has a more restricted expression pattern. Both isoforms were expressed in the lung epithelium. We observed aberrant expression levels of NPSR1-B in smooth muscle in asthmatic bronchi as compared to healthy. In an experimental mouse model, the induced lung inflammation resulted in elevated Npsr1 levels. Furthermore, we demonstrated that the activation of NPSR1 with its endogenous agonist, neuropeptide S (NPS), resulted in a significant inhibition of the growth of NPSR1-A overexpressing stable cell lines (NPSR1-A cells). To determine which target genes were regulated by the NPS-NPSR1 pathway, NPSR1-A cells were stimulated with NPS, and differentially expressed genes were identified using the Affymetrix HGU133Plus2 GeneChip. A total of 104 genes were found significantly up-regulated and 42 down-regulated 6 h after NPS administration. The up-regulated genes included many neuronal genes and some putative susceptibility genes for respiratory disorders. By Gene Ontology enrichment analysis, the biological process terms, cell proliferation, morphogenesis and immune response were among the most altered. The expression of four up-regulated genes, matrix metallopeptidase 10 (MMP10), INHBA (activin A), interleukin 8 (IL8) and EPH receptor A2 (EPHA2), were verified and confirmed by quantitative reverse-transcriptase-PCR. In conclusion, we identified a novel asthma susceptibility gene, NPSR1, on chromosome 7p14.3. NPS-NPSR1 represents a novel pathway that regulates cell proliferation and immune responses, and thus may have functional relevance in the pathogenesis of asthma.