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

Buffered memory: a hypothesis for the maintenance of functional, virus-specific CD8(+) T cells during cytomegalovirus infection.

Chronic infections have been a major topic of investigation in recent years, but the mechanisms that dictate whether or not a pathogen is successfully controlled are incompletely understood. Cytomegalovirus (CMV) is a herpesvirus that establishes a persistent infection in the majority of people in the world. Like other herpesviruses, CMV is well controlled by an effective immune response and induces little, if any, pathology in healthy individuals. However, controlling CMV requires continuous immune surveillance, and thus, CMV is a significant cause of morbidity and death in immune-compromised individuals. T cells in particular play an important role in controlling CMV and both CD4(+) and CD8(+) CMV-specific T cells are essential. These virus-specific T cells persist in exceptionally large numbers during the infection, traffic into peripheral tissues and remain functional, making CMV an attractive vaccine vector for driving CMV-like T cell responses against recombinant antigens of choice. However, the mechanisms by which these T cells persist and differentiate while remaining functional are still poorly understood, and we have no means to promote their development in immune-compromised patients at risk for CMV disease. In this review, I will briefly summarize our current knowledge of CMV-specific CD8(+) T cells and propose a mechanism that may explain their maintenance and preservation of function during chronic infection

The role of multiple marks in epigenetic silencing and the emergence of a stable bivalent chromatin state

We introduce and analyze a minimal model of epigenetic silencing in budding yeast, built upon known biomolecular interactions in the system. Doing so, we identify the epigenetic marks essential for the bistability of epigenetic states. The model explicitly incorporates two key chromatin marks, namely H4K16 acetylation and H3K79 methylation, and explores whether the presence of multiple marks lead to a qualitatively different systems behavior. We find that having both modifications is important for the robustness of epigenetic silencing. Besides the silenced and transcriptionally active fate of chromatin, our model leads to a novel state with bivalent (i.e., both active and silencing) marks under certain perturbations (knock-out mutations, inhibition or enhancement of enzymatic activity). The bivalent state appears under several perturbations and is shown to result in patchy silencing. We also show that the titration effect, owing to a limited supply of silencing proteins, can result in counter-intuitive responses. The design principles of the silencing system is systematically investigated and disparate experimental observations are assessed within a single theoretical framework. Specifically, we discuss the behavior of Sir protein recruitment, spreading and stability of silenced regions in commonly-studied mutants (e.g., sas2, dot1) illuminating the controversial role of Dot1 in the systems biology of yeast silencing.Comment: Supplementary Material, 14 page

Connectivity of larval stages of sedentary marine communities between hard substrates and offshore structures in the North Sea

Man-made structures including rigs, pipelines, cables, renewable energy devices, and ship wrecks, offer hard substrate in the largely soft-sediment environment of the North Sea. These structures become colonised by sedentary organisms and non-migratory reef fish, and form local ecosystems that attract larger predators including seals, birds, and fish. It is possible that these structures form a system of interconnected reef environments through the planktonic dispersal of the pelagic stages of organisms by ocean currents. Changes to the overall arrangement of hard substrate areas through removal or addition of individual man-made structures will affect the interconnectivity and could impact on the ecosystem. Here, we assessed the connectivity of sectors with oil and gas structures, wind farms, wrecks, and natural hard substrate, using a model that simulates the drift of planktonic stages of seven organisms with sedentary adult stages associated with hard substrate, applied to the period 2001–2010. Connectivity was assessed using a classification system designed to address the function of sectors in the network. Results showed a relatively stable overall spatial distribution of sector function but with distinct variations between species and years. The results are discussed in the context of decommissioning of oil and gas infrastructure in the North Sea

Azimuthal anisotropy and correlations at large transverse momenta in p+pp+p and Au+Au collisions at sNN\sqrt{s_{_{NN}}}= 200 GeV

Results on high transverse momentum charged particle emission with respect to the reaction plane are presented for Au+Au collisions at $\sqrt{s_{_{NN}}}$= 200 GeV. Two- and four-particle correlations results are presented as well as a comparison of azimuthal correlations in Au+Au collisions to those in $p+p$ at the same energy. Elliptic anisotropy, $v_2$, is found to reach its maximum at $p_t \sim 3$ GeV/c, then decrease slowly and remain significant up to $p_t\approx 7$ -- 10 GeV/c. Stronger suppression is found in the back-to-back high-$p_t$ particle correlations for particles emitted out-of-plane compared to those emitted in-plane. The centrality dependence of $v_2$ at intermediate $p_t$ is compared to simple models based on jet quenching.Comment: 4 figures. Published version as PRL 93, 252301 (2004

Systemic hematogenous maintenance of memory inflation by MCMV infection.

Several low-grade persistent viral infections induce and sustain very large numbers of virus-specific effector T cells. This was first described as a response to cytomegalovirus (CMV), a herpesvirus that establishes a life-long persistent/latent infection, and sustains the largest known effector T cell populations in healthy people. These T cells remain functional and traffic systemically, which has led to the recent exploration of CMV as a persistent vaccine vector. However, the maintenance of this remarkable response is not understood. Current models propose that reservoirs of viral antigen and/or latently infected cells in lymph nodes stimulate T cell proliferation and effector differentiation, followed by migration of progeny to non-lymphoid tissues where they control CMV reactivation. We tested this model using murine CMV (MCMV), a natural mouse pathogen and homologue of human CMV (HCMV). While T cells within draining lymph nodes divided at a higher rate than cells elsewhere, antigen-dependent proliferation of MCMV-specific effector T cells was observed systemically. Strikingly, inhibition of T cell egress from lymph nodes failed to eliminate systemic T cell division, and did not prevent the maintenance of the inflationary populations. In fact, we found that the vast majority of inflationary cells, including most cells undergoing antigen-driven division, had not migrated into the parenchyma of non-lymphoid tissues but were instead exposed to the blood supply. Indeed, the immunodominance and effector phenotype of inflationary cells, both of which are primary hallmarks of memory inflation, were largely confined to blood-localized T cells. Together these results support a new model of MCMV-driven memory inflation in which most immune surveillance occurs in circulation, and in which most inflationary effector T cells are produced in response to viral antigen presented by cells that are accessible to the blood supply

Azimuthal anisotropy in Au+Au collisions at sqrtsNN = 200 GeV

The results from the STAR Collaboration on directed flow (v_1), elliptic flow (v_2), and the fourth harmonic (v_4) in the anisotropic azimuthal distribution of particles from Au+Au collisions at sqrtsNN = 200 GeV are summarized and compared with results from other experiments and theoretical models. Results for identified particles are presented and fit with a Blast Wave model. Different anisotropic flow analysis methods are compared and nonflow effects are extracted from the data. For v_2, scaling with the number of constituent quarks and parton coalescence is discussed. For v_4, scaling with v_2^2 and quark coalescence is discussed.Comment: 26 pages. As accepted by Phys. Rev. C. Text rearranged, figures modified, but data the same. However, in Fig. 35 the hydro calculations are corrected in this version. The data tables are available at http://www.star.bnl.gov/central/publications/ by searching for "flow" and then this pape

Mindful breath awareness meditation facilitates efficiency gains in brain networks: A steady-state visually evoked potentials study

The beneficial effects of mindfulness-based therapeutic interventions have stimulated a rapidly growing body of scientific research into underlying psychological processes. Resulting evidence indicates that engaging with mindfulness meditation is associated with increased performance on a range of cognitive tasks. However, the mechanisms promoting these improvements require further investigation. We studied changes in behavioural performance of 34 participants during a multiple object tracking (MOT) task that taps core cognitive processes, namely sustained selective visual attention and spatial working memory. Concurrently, we recorded the steady-state visually evoked potential (SSVEP), an EEG signal elicited by the continuously flickering moving objects, and indicator of attentional engagement. Participants were tested before and after practicing eight weeks of mindful breath awareness meditation or progressive muscle relaxation as active control condition. The meditation group improved their MOT-performance and exhibited a reduction of SSVEP amplitudes, whereas no such changes were observed in the relaxation group. Neither group changed in self-reported positive affect and mindfulness, while a marginal increase in negative affect was observed in the mindfulness group. This novel way of combining MOT and SSVEP provides the important insight that mindful breath awareness meditation may lead to refinements of attention networks, enabling more efficient use of attentional resources

Self-citations at the meso and individual levels: effects of different calculation methods

This paper focuses on the study of self-citations at the meso and micro (individual) levels, on the basis of an analysis of the production (1994–2004) of individual researchers working at the Spanish CSIC in the areas of Biology and Biomedicine and Material Sciences. Two different types of self-citations are described: author self-citations (citations received from the author him/herself) and co-author self-citations (citations received from the researchers’ co-authors but without his/her participation). Self-citations do not play a decisive role in the high citation scores of documents either at the individual or at the meso level, which are mainly due to external citations. At micro-level, the percentage of self-citations does not change by professional rank or age, but differences in the relative weight of author and co-author self-citations have been found. The percentage of co-author self-citations tends to decrease with age and professional rank while the percentage of author self-citations shows the opposite trend. Suppressing author self-citations from citation counts to prevent overblown self-citation practices may result in a higher reduction of citation numbers of old scientists and, particularly, of those in the highest categories. Author and co-author self-citations provide valuable information on the scientific communication process, but external citations are the most relevant for evaluative purposes. As a final recommendation, studies considering self-citations at the individual level should make clear whether author or total self-citations are used as these can affect researchers differently