Chaste: an open source C++ library for computational physiology and biology

Chaste - Cancer, Heart And Soft Tissue Environment - is an open source C++ library for the computational simulation of mathematical models developed for physiology and biology. Code development has been driven by two initial applications: cardiac electrophysiology and cancer development. A large number of cardiac electrophysiology studies have been enabled and performed, including high performance computational investigations of defibrillation on realistic human cardiac geometries. New models for the initiation and growth of tumours have been developed. In particular, cell-based simulations have provided novel insight into the role of stem cells in the colorectal crypt. Chaste is constantly evolving and is now being applied to a far wider range of problems. The code provides modules for handling common scientific computing components, such as meshes and solvers for ordinary and partial differential equations (ODEs/PDEs). Re-use of these components avoids the need for researchers to "re-invent the wheel" with each new project, accelerating the rate of progress in new applications. Chaste is developed using industrially-derived techniques, in particular test-driven development, to ensure code quality, re-use and reliability. In this article we provide examples that illustrate the types of problems Chaste can be used to solve, which can be run on a desktop computer. We highlight some scientific studies that have used or are using Chaste, and the insights they have provided. The source code, both for specific releases and the development version, is available to download under an open source Berkeley Software Distribution (BSD) licence at http://www.cs.ox.ac.uk/chaste, together with details of a mailing list and links to documentation and tutorials

Predator-Induced Vertical Behavior of a Ctenophore

Although many studies have focused on Mnemiopsis leidyi predation, little is known about the role of this ctenophore as prey when abundant in native and invaded pelagic systems. We examined the response of the ctenophore M. leidyi to the predatory ctenophore Beroe ovata in an experiment in which the two species could potentially sense each other while being physically separated. On average, M. leidyi responded to the predator’s presence by increasing variability in swimming speeds and by lowering their vertical distribution. Such behavior may help explain field records of vertical migration, as well as stratified and near-bottom distributions of M. leidyi

Ten Simple Rules for Effective Computational Research

<p>Ten Simple Rules for Effective Computational Research</p

A Multicellular Model of Intestinal Crypt Buckling and Fission

Crypt fission is an in vivo tissue deformation process that is involved in both intestinal homeostasis and colorectal tumourigenesis. Despite its importance, the mechanics underlying crypt fission are currently poorly understood. Recent experimental development of organoids, organ-like buds cultured from crypt stem cells in vitro, has shown promise in shedding light on crypt fission. Drawing inspiration from observations of organoid growth and fission in vivo, we develop a computational model of a deformable epithelial tissue layer. Results from in silico experiments show the stiffness of cells and the proportions of cell subpopulations affect the nature of deformation in the epithelial layer. In particular, we find that increasing the proportion of stiffer cells in the layer increases the likelihood of crypt fission occurring. This is in agreement with and helps explain recent experimental work

Comparing individual-based approaches to modelling the self-organization of multicellular tissues.

The coordinated behaviour of populations of cells plays a central role in tissue growth and renewal. Cells react to their microenvironment by modulating processes such as movement, growth and proliferation, and signalling. Alongside experimental studies, computational models offer a useful means by which to investigate these processes. To this end a variety of cell-based modelling approaches have been developed, ranging from lattice-based cellular automata to lattice-free models that treat cells as point-like particles or extended shapes. However, it remains unclear how these approaches compare when applied to the same biological problem, and what differences in behaviour are due to different model assumptions and abstractions. Here, we exploit the availability of an implementation of five popular cell-based modelling approaches within a consistent computational framework, Chaste (http://www.cs.ox.ac.uk/chaste). This framework allows one to easily change constitutive assumptions within these models. In each case we provide full details of all technical aspects of our model implementations. We compare model implementations using four case studies, chosen to reflect the key cellular processes of proliferation, adhesion, and short- and long-range signalling. These case studies demonstrate the applicability of each model and provide a guide for model usage

Underrepresentation of Elderly People in Randomised Controlled Trials. The Example of Trials of 4 Widely Prescribed Drugs

BACKGROUND: We aimed to determine the representation of elderly people in published reports of randomized controlled trials (RCTs). We focused on trials of 4 medications--pioglitazone, rosuvastatin, risedronate, and valsartan-frequently used by elderly patients with chronic medical conditions. METHODS AND FINDINGS: We selected all reports of RCTs indexed in PubMed from 1966 to April 2008 evaluating one of the 4 medications of interest. Estimates of the community-based "on-treatment" population were from a national health insurance database (SNIIR-AM) covering approximately 86% of the population in France. From this database, we evaluated data claims from January 2006 to December 2007 for 1,958,716 patients who received one of the medications of interest for more than 6 months. Of the 155 RCT reports selected, only 3 studies were exclusively of elderly patients (2 assessing valsartan; 1 risedronate). In only 4 of 37 reports (10.8%) for pioglitazone, 4 of 22 (18.2%) for risedronate, 3 of 29 (10.3%) for rosuvastatine and 9 of 67 (13.4%) for valsartan, the proportion of patients aged 65 or older was within or above that treated in clinical practice. In 62.2% of the reports for pioglitazone, 40.9% for risedronate, 37.9% for rosuvastatine, and 70.2% for valsartan, the proportion of patients aged 65 or older was lower than half that in the treated population. The representation of elderly people did not differ by publication date or sample size. CONCLUSIONS: Elderly patients are poorly represented in RCTs of drugs they are likely to receive

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

Search for squarks and gluinos in events with isolated leptons, jets and missing transverse momentum at s√=8 TeV with the ATLAS detector

The results of a search for supersymmetry in final states containing at least one isolated lepton (electron or muon), jets and large missing transverse momentum with the ATLAS detector at the Large Hadron Collider are reported. The search is based on proton-proton collision data at a centre-of-mass energy s√=8 TeV collected in 2012, corresponding to an integrated luminosity of 20 fb−1. No significant excess above the Standard Model expectation is observed. Limits are set on supersymmetric particle masses for various supersymmetric models. Depending on the model, the search excludes gluino masses up to 1.32 TeV and squark masses up to 840 GeV. Limits are also set on the parameters of a minimal universal extra dimension model, excluding a compactification radius of 1/R c = 950 GeV for a cut-off scale times radius (ΛR c) of approximately 30

Search for squarks and gluinos with the ATLAS detector in final states with jets and missing transverse momentum using √s=8 TeV proton-proton collision data

A search for squarks and gluinos in final states containing high-p T jets, missing transverse momentum and no electrons or muons is presented. The data were recorded in 2012 by the ATLAS experiment in s√=8 TeV proton-proton collisions at the Large Hadron Collider, with a total integrated luminosity of 20.3 fb−1. Results are interpreted in a variety of simplified and specific supersymmetry-breaking models assuming that R-parity is conserved and that the lightest neutralino is the lightest supersymmetric particle. An exclusion limit at the 95% confidence level on the mass of the gluino is set at 1330 GeV for a simplified model incorporating only a gluino and the lightest neutralino. For a simplified model involving the strong production of first- and second-generation squarks, squark masses below 850 GeV (440 GeV) are excluded for a massless lightest neutralino, assuming mass degenerate (single light-flavour) squarks. In mSUGRA/CMSSM models with tan β = 30, A 0 = −2m 0 and μ > 0, squarks and gluinos of equal mass are excluded for masses below 1700 GeV. Additional limits are set for non-universal Higgs mass models with gaugino mediation and for simplified models involving the pair production of gluinos, each decaying to a top squark and a top quark, with the top squark decaying to a charm quark and a neutralino. These limits extend the region of supersymmetric parameter space excluded by previous searches with the ATLAS detector

Phthiocerol Dimycocerosates of M. tuberculosis Participate in Macrophage Invasion by Inducing Changes in the Organization of Plasma Membrane Lipids

Phthiocerol dimycocerosates (DIM) are major virulence factors of Mycobacterium tuberculosis (Mtb), in particular during the early step of infection when bacilli encounter their host macrophages. However, their cellular and molecular mechanisms of action remain unknown. Using Mtb mutants deleted for genes involved in DIM biosynthesis, we demonstrated that DIM participate both in the receptor-dependent phagocytosis of Mtb and the prevention of phagosomal acidification. The effects of DIM required a state of the membrane fluidity as demonstrated by experiments conducted with cholesterol-depleting drugs that abolished the differences in phagocytosis efficiency and phagosome acidification observed between wild-type and mutant strains. The insertion of a new cholesterol-pyrene probe in living cells demonstrated that the polarity of the membrane hydrophobic core changed upon contact with Mtb whereas the lateral diffusion of cholesterol was unaffected. This effect was dependent on DIM and was consistent with the effect observed following DIM insertion in model membrane. Therefore, we propose that DIM control the invasion of macrophages by Mtb by targeting lipid organisation in the host membrane, thereby modifying its biophysical properties. The DIM-induced changes in lipid ordering favour the efficiency of receptor-mediated phagocytosis of Mtb and contribute to the control of phagosomal pH driving bacilli in a protective niche