Using argument notation to engineer biological simulations with increased confidence

The application of computational and mathematical modelling to explore the mechanics of biological systems is becoming prevalent. To significantly impact biological research, notably in developing novel therapeutics, it is critical that the model adequately represents the captured system. Confidence in adopting in silico approaches can be improved by applying a structured argumentation approach, alongside model development and results analysis. We propose an approach based on argumentation from safety-critical systems engineering, where a system is subjected to a stringent analysis of compliance against identified criteria. We show its use in examining the biological information upon which a model is based, identifying model strengths, highlighting areas requiring additional biological experimentation and providing documentation to support model publication. We demonstrate our use of structured argumentation in the development of a model of lymphoid tissue formation, specifically Peyer's Patches. The argumentation structure is captured using Artoo (www.york.ac.uk/ycil/software/artoo), our Web-based tool for constructing fitness-for-purpose arguments, using a notation based on the safety-critical goal structuring notation. We show how argumentation helps in making the design and structured analysis of a model transparent, capturing the reasoning behind the inclusion or exclusion of each biological feature and recording assumptions, as well as pointing to evidence supporting model-derived conclusions

An ecological approach to anomaly detection: the EIA Model.

The presented work proposes a new approach for anomaly detection. This approach is based on changes in a population of evolving agents under stress. If conditions are appropriate, changes in the population (modeled by the bioindicators) are representative of the alterations to the environment. This approach, based on an ecological view, improves functionally traditional approaches to the detection of anomalies. To verify this assertion, experiments based on Network Intrussion Detection Systems are presented. The results are compared with the behaviour of other bioinspired approaches and machine learning techniques

Management of infants with brief resolved unexplained events (Brue) and apparent life-threatening events (alte): A rand/ucla appropriateness approach

Unexpected events of breath, tone, and skin color change in infants are a cause of consider-able distress to the caregiver and there is still debate on their appropriate management. The aim of this study is to survey the trend in prevention, decision-making, and management of brief resolved unexplained events (BRUE)/apparent life-threatening events (ALTE) and to develop a shared proto-col among hospitals and primary care pediatricians regarding hospital admission criteria, work-up and post-discharge monitoring of patients with BRUE/ALTE. For the study purpose, a panel of 54 experts was selected to achieve consensus using the RAND/UCLA appropriateness method. Twelve scenarios were developed: one addressed to primary prevention of ALTE and BRUE, and 11 focused on hospital management of BRUE and ALTE. For each scenario, participants were asked to rank each option from ‘1’ (extremely inappropriate) to ‘9’ (extremely appropriate). Results derived from panel meeting and discussion showed several points of agreement but also disagreement with different opinion emerged and the need of focused education on some areas. However, by combining previous recommendations with expert opinion, the application of the RAND/UCLA appropriateness permit-ted us to drive pediatricians to reasoned and informed decisions in term of evaluation, treatment and follow-up of infants with BRUE/ALTE, reducing inappropriate exams and hospitalisation and highlighting priorities for educational interventions

Statistical Techniques Complement UML When Developing Domain Models of Complex Dynamical Biosystems

Computational modelling and simulation is increasingly being used to complement traditional wet-lab techniques when investigating the mechanistic behaviours of complex biological systems. In order to ensure computational models are fit for purpose, it is essential that the abstracted view of biology captured in the computational model, is clearly and unambiguously defined within a conceptual model of the biological domain (a domain model), that acts to accurately represent the biological system and to document the functional requirements for the resultant computational model. We present a domain model of the IL-1 stimulated NF-κB signalling pathway, which unambiguously defines the spatial, temporal and stochastic requirements for our future computational model. Through the development of this model, we observe that, in isolation, UML is not sufficient for the purpose of creating a domain model, and that a number of descriptive and multivariate statistical techniques provide complementary perspectives, in particular when modelling the heterogeneity of dynamics at the single-cell level. We believe this approach of using UML to define the structure and interactions within a complex system, along with statistics to define the stochastic and dynamic nature of complex systems, is crucial for ensuring that conceptual models of complex dynamical biosystems, which are developed using UML, are fit for purpose, and unambiguously define the functional requirements for the resultant computational model

Unlocking the potential of publicly available microarray data using inSilicoDb and inSilicoMerging R/Bioconductor packages

BACKGROUND: With an abundant amount of microarray gene expression data sets available through public repositories, new possibilities lie in combining multiple existing data sets. In this new context, analysis itself is no longer the problem, but retrieving and consistently integrating all this data before delivering it to the wide variety of existing analysis tools becomes the new bottleneck. RESULTS: We present the newly released inSilicoMerging R/Bioconductor package which, together with the earlier released inSilicoDb R/Bioconductor package, allows consistent retrieval, integration and analysis of publicly available microarray gene expression data sets. Inside the inSilicoMerging package a set of five visual and six quantitative validation measures are available as well. CONCLUSIONS: By providing (i) access to uniformly curated and preprocessed data, (ii) a collection of techniques to remove the batch effects between data sets from different sources, and (iii) several validation tools enabling the inspection of the integration process, these packages enable researchers to fully explore the potential of combining gene expression data for downstream analysis. The power of using both packages is demonstrated by programmatically retrieving and integrating gene expression studies from the InSilico DB repository [https://insilicodb.org/app/]

Validating Antimetastatic Effects of Natural Products in an Engineered Microfluidic Platform Mimicking Tumor Microenvironment

Development of new, antimetastatic drugs from natural products has been substantially constrained by the lack of a reliable in vitro screening system. Such a system should ideally mimic the native, three-dimensional (3D) tumor microenvironment involving different cell types and allow quantitative analysis of cell behavior critical for metastasis. These requirements are largely unmet in the current model systems, leading to poor predictability of the in vitro collected data for in vivo trials, as well as prevailing inconsistency among different in vitro tests. In the present study, we report application of a 3D, microfluidic device for validation of the antimetastatic effects of 12 natural compounds. This system supports co-culture of endothelial and cancer cells in their native 3D morphology as in the tumor microenvironment and provides real-time monitoring of the cells treated with each compound. We found that three compounds, namely sanguinarine, nitidine, and resveratrol, exhibited significant antimetastatic or antiangiogenic effects. Each compound was further examined for its respective activity with separate conventional biological assays, and the outcomes were in agreement with the findings collected from the microfluidic system. In summary, we recommend use of this biomimetic model system as a new engineering tool for high-throughput evaluation of more diverse natural compounds with varying anticancer potentials

Representation in the (Artificial) Immune System

Much of contemporary research in Artificial Immune Systems (AIS) has partitioned into either algorithmic machine learning and optimisation, or, modelling biologically plausible dynamical systems, with little overlap between. We propose that this dichotomy is somewhat to blame for the lack of significant advancement of the field in either direction and demonstrate how a simplistic interpretation of Perelson’s shape-space formalism may have largely contributed to this dichotomy. In this paper, we motivate and derive an alternative representational abstraction. To do so we consider the validity of shape-space from both the biological and machine learning perspectives. We then take steps towards formally integrating these perspectives into a coherent computational model of notions such as life-long learning, degeneracy, constructive representations and contextual recognition—rhetoric that has long inspired work in AIS, while remaining largely devoid of operational definition

Detecting Network Communities: An Application to Phylogenetic Analysis

This paper proposes a new method to identify communities in generally weighted complex networks and apply it to phylogenetic analysis. In this case, weights correspond to the similarity indexes among protein sequences, which can be used for network construction so that the network structure can be analyzed to recover phylogenetically useful information from its properties. The analyses discussed here are mainly based on the modular character of protein similarity networks, explored through the Newman-Girvan algorithm, with the help of the neighborhood matrix . The most relevant networks are found when the network topology changes abruptly revealing distinct modules related to the sets of organisms to which the proteins belong. Sound biological information can be retrieved by the computational routines used in the network approach, without using biological assumptions other than those incorporated by BLAST. Usually, all the main bacterial phyla and, in some cases, also some bacterial classes corresponded totally (100%) or to a great extent (>70%) to the modules. We checked for internal consistency in the obtained results, and we scored close to 84% of matches for community pertinence when comparisons between the results were performed. To illustrate how to use the network-based method, we employed data for enzymes involved in the chitin metabolic pathway that are present in more than 100 organisms from an original data set containing 1,695 organisms, downloaded from GenBank on May 19, 2007. A preliminary comparison between the outcomes of the network-based method and the results of methods based on Bayesian, distance, likelihood, and parsimony criteria suggests that the former is as reliable as these commonly used methods. We conclude that the network-based method can be used as a powerful tool for retrieving modularity information from weighted networks, which is useful for phylogenetic analysis