The T cell receptor triggering apparatus is composed of monovalent or monomeric proteins.

Understanding the component stoichiometry of the T cell antigen receptor (TCR) triggering apparatus is essential for building realistic models of signal initiation. Recent studies suggesting that the TCR and other signaling-associated proteins are preclustered on resting T cells relied on measurements of the behavior of membrane proteins at interfaces with functionalized glass surfaces. Using fluorescence recovery after photobleaching, we show that, compared with the apical surface, the mobility of TCRs is significantly reduced at Jurkat T cell/glass interfaces, in a signaling-sensitive manner. Using two biophysical approaches that mitigate these effects, bioluminescence resonance energy transfer and two-color coincidence detection microscopy, we show that, within the uncertainty of the methods, the membrane components of the TCR triggering apparatus, i.e. the TCR complex, MHC molecules, CD4/Lck and CD45, are exclusively monovalent or monomeric in human T cell lines, implying that TCR triggering depends only on the kinetics of TCR/pMHC interactions. These analyses also showed that constraining proteins to two dimensions at the cell surface greatly enhances random interactions versus those between the membrane and the cytoplasm. Simulations of TCR-pMHC complex formation based on these findings suggest how unclustered TCR triggering-associated proteins might nevertheless be capable of generating complex signaling outputs via the differential recruitment of cytosolic effectors to the cell membrane.

Dynamic modelling of cell death during biofilm development.

Biofilms are currently recognised as the predominant bacterial life-style and it has been suggested that biofilm development is influenced by a number of different processes such as adhesion, detachment, mass transport, quorum sensing, cell death and active dispersal. One of the least understood processes and its effects on biofilm development is cell death. However, experimental studies suggest that bacterial death is an important process during biofilm development and many studies show a relationship between cell death and dispersal in microbial biofilms. We present a model of the process of cell death during biofilm development, with a particular focus on the spatial localisation of cell death or cell damage. Three rules governing cell death or cell damage were evaluated which compared the effects of starvation, damage accumulation, and viability during biofilm development and were also used to design laboratory based experiments to test the model. Results from model simulations show that actively growing biofilms develop steep nutrient gradients within the interior of the biofilm that affect neighbouring microcolonies resulting in cell death and detachment. Two of the rules indicated that high substrate concentrations lead to accelerated cell death, in contrast to the third rule, based on the accumulation of damage, which predicted earlier cell death for biofilms grown with low substrate concentrations. Comparison of the modelling results with experimental results suggests that cell death is favoured under low nutrient conditions and that the accumulation of damage may be the main cause of cell death during biofilm development.

Mathematical modelling of the regulation of spa (protein A) transcription in Staphylococcus aureus.

In the present work a general systems biology approach has been used to study the complex regulatory network controlling the transcription of the spa gene, encoding protein A, a major surface protein and an important virulence factor of Staphylococcus aureus. A valid mathematical model could be formulated using parameter values, which were fitted to quantitative Northern blot data from various S. aureus regulatory mutants using a gradient search method. The model could correctly predict spa expression levels in 4 different regulatory mutants not included in the parameter value search, and in 2 other S. aureus strains, SH1000 and UAMS-1. The mathematical model revealed that sarA and sarS seem to balance each other in a way that when the activating impact of sarS is small, e.g. in the wild-type, the repressive impact of sarA is small, while in an agr-deficient background, when the impact of sarS is maximal, the repressive impact of sarA is close to its maximum. Furthermore, the model revealed that Rot and SarS act synergistically to stimulate spa expression, something that was not obvious from experimental data. We believe that this mathematical model can be used to evaluate the significance of other putative interactions in the regulatory network governing spa transcription.

Invasive Staphylococcus aureus strains are highly variable in PFGE patterns, agr group and exoprotein production.

In the present study, we have investigated 37 invasive Staphylococcus aureus strains (collected between 1997 and 2005) from 33 human episodes of septicaemia causing either endocarditis or vertebral osteomyelitis. All S. aureus strains were typed using pulsed field gel electrophoresis (PFGE), and most strains belonged to any of 4 different PFGE clusters. There was no correlation between any of the PFGE clusters with site of infection. All strains showed highly different expression patterns of extracellular proteins, i.e. we found a vast variation in the number of proteins and amount of individual proteins expressed by the different strains. There was no correlation between any cluster of exoprotein patterns with endocarditis or with vertebral osteomyelitis. We did not find any correlation between agr group and endocarditis, as previously reported. On the other hand, a correlation between some of the PFGE clusters with a certain agr group was found. Known risk factors for S. aureus infections were observed in a majority of the patients.

Cardiomyogenic gene expression profiling of differentiating human embryonic stem cells.

Human embryonic stem cells (hESCs) can differentiate into a variety of specialized cell types. Thus, they provide a model system for embryonic development to investigate the molecular processes of cell differentiation and lineage commitment. The development of the cardiac lineage is easily detected in mixed cultures by the appearance of spontaneously contracting areas of cells. We performed gene expression profiling of undifferentiated and differentiating hESCs and monitored 468 genes expressed during cardiac development and/or in cardiac tissue. Their transcription during early differentiation of hESCs through embryoid bodies (EBs) was investigated and compared with spontaneously differentiating hESCs maintained on feeders in culture without passaging (high-density (HD) protocol). We observed a larger variation in the gene expression between cells from a single cell line that were differentiated using two different protocols than in cells from different cell lines that were cultured according to the same protocol. Notably, the EB protocol resulted in more reproducible transcription profiles than the HD protocol. The results presented here provide new information about gene regulation during early differentiation of hESCs with emphasis on the cardiomyogenic program. In addition, we also identified regulatory elements that could prove critical for the development of the cardiomyocyte lineage.

An individual-based network model to evaluate interventions for controlling pneumococcal transmission.

BACKGROUND: Streptococcus pneumoniae is a major cause of morbidity and mortality worldwide, but also a common colonizer of the upper respiratory tract. The emergence and spread of antibiotic resistant pneumococcal strains has threatened effective therapy. The long-term effects of measures aiming to limit pneumococcal spread are poorly understood. Computational modeling makes it possible to conduct virtual experiments that are impractical to perform in real life and thereby allows a more full understanding of pneumococcal epidemiology and control efforts. METHODS: We have developed a contact network model to evaluate the efficacy of interventions aiming to control pneumococcal transmission. Demographic data from Sweden during the mid-2000s were employed. Analyses of the model's parameters were conducted to elucidate key determinants of pneumococcal spread. Also, scenario simulations were performed to assess candidate control measures. RESULTS: The model made good predictions of previous findings where a correlation has been found between age and pneumococcal carriage. Of the parameters tested, group size in day-care centers was shown to be one of the most important factors for pneumococcal transmission. Consistent results were generated from the scenario simulations. CONCLUSION: We recommend, based on the model predictions, that strategies to control pneumococcal disease and organism transmission should include reducing the group size in day-care centers.

Modeling the regulation of the competence-evoking quorum sensing network in Streptococcus pneumoniae.

Competence for genetic transformation seems to play a fundamental role in the biology of Streptococcus pneumoniae and is believed to account for serotype switching, evolution of virulence factors, and rapid emergence of antibiotic resistance. The initiation of competence is regulated by the quorum sensing system referred as the ComABCDE pathway. Experimental studies reveal that competence is down-regulated a short time after its induction and several hypotheses about the mechanism(s) responsible for this shut-down have been presented. Possibly, a ComX-dependent gene product, such as a repressor or a phosphatase, is involved. To better understand the down-regulation of the competence-evoking system in S. pneumoniae, a mathematical model was set up. By analyzing the model, we suggest that shut-down of competence possibly occurs at the transcriptional level on the comCDE operon. As a result of introducing a putative comX-dependent repressor, which inhibits expression of comCDE and comX, in the mathematical model, competence is demonstrated to appear in waves. This is supported by experimental studies showing the appearance of successive competence cycles in pneumococcal batch cultures.

Modeling the effect of acylated homoserine lactone antagonists in Pseudomonas aeruginosa.

Pseudomonas aeruginosa is a gram-negative bacterium that causes serious illnesses, particularly in immunocompromised individuals, often with a fatal outcome. The finding that the acylated homoserine lactone quorum sensing (QS) system controls the production of virulence factors in P. aeruginosa makes this system a possible target for antimicrobial therapy. It has been suggested that an N-(3-oxododecanoyl)-homoserine lactone (3O-C12-HSL) antagonist, a QS blocker (QSB), would interfere efficiently with the quorum sensing system in P. aeruginosa and thus reduce the virulence of this pathogen. In this work, a mathematical model of the QS system in P. aeruginosa has been developed. The model was used to virtually add 3O-C12-HSL antagonists that differed in their affinity for the receptor protein and for their ability to mediate degradation of the receptor. The model suggests that very small differences in these parameters for different 3O-C12-HSL antagonists can greatly affect the success of QSB based inhibition of the QS system in P. aeruginosa. Most importantly, it is proposed that the ability of the 3O-C12-HSL antagonist to mediate degradation of LasR is the core parameter for successful QSB based inhibition of the QS system in P. aeruginosa. Finally, this study demonstrates that QSBs can shift the system to a low steady state, corresponding to an uninduced state and thus, suggests that the use of 3O-C12-HSL antagonists may constitute a promising therapeutic approach against P. aeruginosa involved infections.

Asynchronous population dynamics of Siberian lemmings across the Palaearctic tundra.

The synchrony of Siberian lemming (Lemmus sibiricus L.) population dynamics was investigated during a ship-borne expedition along the Palaearctic tundra coast in the summer of 1994. On 12 sites along the coast from the Kola Peninsula to Wrangel Island, relative densities of lemmings were recorded using a standardised snap-trapping programme. The phase position of the lemming cycle in each of the studied populations was determined based on current density estimates, signs of previous density and the age profile of each population (ageing based on eye lens mass). In addition, dendrochronological methods were used to determine when the last peak in the density of microtine populations occurred at each site. The examined lemming populations were in different phases of the lemming cycle. Some populations were in the peak phase, as indicated by high current densities, an age profile in which older individuals were well represented, and signs of high previous density (abundant old lemming faeces). Other populations were in the decline phase, as reflected in a moderate current density, a predominance of older individuals and signs of high previous density. Populations in the low phase had an extremely low current density and showed signs of high previous density, while populations in the increase phase had a moderate current density, a predominance of younger individuals and showed signs of low previous density. The results of phase determinations based on dendrochronological methods support the findings based on lemming demography. Recent Russian studies carried out on some of the sites also agreed with our phase determination results. Thus, on a regional scale (across the whole Palaearctic tundra), the population dynamics of Siberian lemmings can be considered asynchronous. However, sites situated adjacent to each other were often phase synchronous, suggesting a more fine-grained pattern of dynamics with synchrony over distances as long as 1000 km or so, e.g. the Yamal and Taymyr Peninsulas.

Toward quantifying the thymic dysfunctional state in mouse models of inflammatory bowel disease.

Inflammatory bowel disease is characterized by a number of immunological alterations, not the least in the T-cell compartment. Numerous animal models of colitis have revealed aberrant thymocyte dynamics associated with skewed thymocyte development. The recent advancements in quantitative methods have proposed critical kinetic alterations in the thymocyte development during the progression of colitis. This review focuses on the aberrant thymocyte dynamics in Gαi2-deficient mice as this mouse model provides most quantitative data of the thymocyte development associated with colitis. Herein, we discuss several dynamic changes during the progression of colitis and propose a hypothesis for the underlying causes for the skewed proportions of the thymocyte populations seen in the Gαi2-deficient mice and in other mouse models of colitis.

3D computation modelling of the influence of cytokine secretion on Th-cell development suggests that negative selection (inhibition of Th1 cells) is more effective than positive selection by IL-4 for Th2 cell dominance.

Th-cell development has been suggested to include selective mechanisms in which certain cytokines select either Th1 or Th2 cells to proliferate and grow. The selective theory is based on the observation that Th2 cells secrete IL-4, a cytokine that promotes Th2 development, whereas Th1 cells secrete interferon-gamma (IFN-gamma) that favours Th1 development, and both positive and negative selective influences have been suggested to operate. In this study, we investigate the role of autocrine secretion and utilization of IL-4 by Th2 cells and address the question of whether an activated Th2 cell can be positively selected by IL-4 secreted from other Th2 cells. We present a spatial three dimensional (3D) modelling approach to simulate the interaction between the IL-4 ligand and its IL-4 receptors expressed on discrete IL-4 secreting cells. The simulations, based on existing experimental data on the IL-4 receptor-ligand system, illustrate how Th-cell development is highly dependent on the distance between cells that are communicating. The model suggests that a single Th2 cell is likely to communicate with possible target cells within a range of approximately 100 microm and that an activated Th2 cell manages to fill most of its own IL-4 receptors, even at a low secretion rate. The predictions made by the model suggest that negative selection against Th1 cells is more effective than positive selection by IL-4 for promoting Th2 dominance.

Differentiating human embryonic stem cells express a unique housekeeping gene signature.

Housekeeping genes (HKGs) are involved in basic functions needed for the sustenance of the cell and are assumed to be constitutively expressed at a constant level. Based on these features, HKGs are frequently used for normalization of gene expression data. In the present study, we used the CodeLink Gene Expression Bioarray system to interrogate changes in gene expression occurring during differentiation of human ESCs (hESCs). Notably, in the three hESC lines used for the study, we observed that the RNA levels of 56 frequently used HKGs varied to a degree that rendered them inappropriate as reference genes. Therefore, we defined a novel set of HKGs specifically for hESCs. Here we present a comprehensive list of 292 genes that are stably expressed (coefficient of variation <20%) in differentiating hESCs. These genes were further grouped into high-, medium-, and low-expressed genes. The expression patterns of these novel HKGs show very little overlap with results obtained from somatic cells and tissues. We further explored the stability of this novel set of HKGs in independent, publicly available gene expression data from hESCs and observed substantial similarities with our results. Gene expression was confirmed by real-time quantitative polymerase chain reaction analysis. Taken together, these results suggest that differentiating hESCs have a unique HKG signature and underscore the necessity to validate the expression profiles of putative HKGs. In addition, this novel set of HKGs can preferentially be used as controls in gene expression analyses of differentiating hESCs.

In silico simulations suggest that Th-cell development is regulated by both selective and instructive mechanisms.

Th-cell differentiation is highly influenced by the local cytokine environment. Although cytokines such as IL-12 and IL-4 are known to polarize the Th-cell response towards Th1 or Th2, respectively, it is not known whether these cytokines instruct the developmental fate of uncommitted Th cells or select cells that have already been committed through a stochastic process. We present an individual based model that accommodates both stochastic and deterministic processes to simulate the dynamic behaviour of selective versus instructive Th-cell development. The predictions made by each model show distinct behaviours, which are compared with experimental observations. The simulations show that the instructive model generates an exclusive Th1 or Th2 response in the absence of an external cytokine source, whereas the selective model favours coexistence of the phenotypes. A hybrid model, including both instructive and selective development, shows behaviour similar to either the selective or the instructive model dependent on the strength of activation. The hybrid model shows the closest qualitative agreement with a number of well-established experimental observations. The predictions by each model suggest that neither pure selective nor instructive Th development is likely to be functional as exclusive mechanisms in Th1/Th2 development.

A theoretical framework for quantitative analysis of the molecular basis of costimulation.

We present a theoretical framework for simulating the synaptic accumulation of the costimulatory molecules CD28, CTLA-4, B7-1, and B7-2, based on a system of mean-field, ordinary differential equations, and rigorous biophysical and expression data. The simulations show that binding affinity, stoichiometric properties, expression levels, and, in particular, competition effects all profoundly influence complex formation at cellular interfaces. B7-2 engages 33-fold more CD28 than CTLA-4 at the synapse in contrast to B7-1, which ligates approximately 7-fold more CTLA-4 than CD28. Although B7-1 completely dominates interactions with CTLA-4, forming linear arrays of 7-18 receptor-ligand pairs, CTLA-4 is fully engaged by B7-2 when B7-1 is absent. Additional simulations reveal the sensitivity of CD28 interactions to modeled transport processes. The results support the concept that B7-2 and B7-1 are the dominant ligands of CD28 and CTLA-4, respectively, and indicate that the inability of B7-2 to recruit CTLA-4 to the synapse cannot be due to the differential binding properties of B7-1 and B7-2 only. We discuss the apparent redundancy of B7-1 in the context of a potentially dynamic synaptic microenvironment, and in light of functions other than the direct enhancement of T cell inhibition by CTLA-4.

Characterizing the dynamics of the quorum-sensing system in Staphylococcus aureus.

The virulence determinants of Staphylococcus aureus are expressed in a growth phase-dependent manner governed by the autoinducible quorum-sensing system agr. Activation of the agr system results in a rapid increase in the regulator RNAIII and occurs in response to accumulation of AIP. In order to activate the agr system, a basal transcription of agr components must be assumed. This basal transcription of agr components seems to be stimulated by sarA. To better understand how SarA would affect activation of the agr system by modulating the basal agr activity, a mathematical model for autoactivation of the agr system was set up. The model predicted that the agr system is hysteretic, meaning that the agr system is activated at a specific concentration of autoinducing peptide (AIP), whereas it is inactivated at a specific lower concentration of AIP. According to the model, changing the basal agr activity only had a marginal effect on steady-state levels of RNAIII but changed the sensitivity of the cells to AIP. This was supported by Northern blot analysis of RNAIII in S. aureus mutants with different levels of SarA expression. Since natural antagonistic AIPs have been demonstrated, the effect of adding inhibitors to the system was analyzed.

Statistical model of the interactions between Helicobacter pylori infection and gastric cancer development.

BACKGROUND: The bacterium Helicobacter pylori is associated with a number of gastrointestinal diseases, such as gastric ulcer, duodenal ulcer and gastric cancer. Several histological changes may be observed during the course of infection; some may influence the progression towards cancer. The aim of this study was to build a statistical model to discover direct interactions between H. pylori and different precancerous changes of the gastric mucosa, and in what order and to what degree those may influence the development of the intestinal type of gastric cancer. METHODS: To find direct and indirect interactions between H. pylori and different histological variables, log-linear analyses were used on a case-control study. To generate mathematically and biologically relevant statistical models, a designed algorithm and observed frequency tables were used. RESULTS: The results show that patients with H. pylori infection need to present with proliferation and intestinal metaplasia to develop gastric cancer of the intestinal type. Proliferation and intestinal metaplasia interacted with the variables atrophy and foveolar hyperplasia. Intestinal metaplasia was the only variable with direct interaction with gastric cancer. Gender had no effect on the variables examined. CONCLUSION: The direct interactions observed in the final statistical model between H. pylori, changes of the mucosa and gastric cancer strengthens and supports previous theories about the progression towards gastric cancer. The results suggest that gastric cancer of the intestinal type may develop from H. pylori infection, proliferation and intestinal metaplasia, while atrophy and foveolar hyperplasia interplay with the other histological variables in the disease process.

The role of regulators in the expression of quorum-sensing signals in Pseudomonas aeruginosa.

Quorum-sensing systems provide Pseudomonas aeruginosa with a sensitive regulatory mechanism that allows for the induction of several phenotypic genes in a cell density fashion. In this work, a mathematical model of the acylated homoserine lactones regulatory network system in P. aeruginosa has been developed. It is the first integrated model to consider both quorum-sensing systems. The model has allowed us to disentangle the complex behavior exhibited by the system as the concentration of extracellular OdDHL is increased. At either low or high levels of extracellular OdDHL, the bacterium remains in an uninduced or induced state, respectively. At moderate levels, the behavior is characterized by several states. Here, the bacteria can switch suddenly from an uninduced to an induced phenotype in response to small changes in the concentration of extracellular OdDHL. Additionally, we have been able to address the roles of RsaL and Vfr as regulators of the quorum-sensing system. An important result from this analysis suggests that RsaL will increase the concentration of extracellular OdDHL required to induce the system, and it is a key regulator of the inhibition of the quorum-sensing system under low cell densities. Most importantly, our results suggest that Vfr has strong regulatory effects on the system as an increased affinity between the LasR/OdDHL complex, and the lasR promoter leads to significant qualitative changes in induction patterns. We also show experimental data that demonstrate that Vfr is required for signal production in the early phase of growth, but that in the latter stages of growth, the vfr mutant is able to synthesize wild-type levels of signal.