Extent of over-expression of hepatocyte growth factor receptor in colorectal tumours is dependent on the choice of normaliser

For reliable results from quantitative RT-PCR, the starting quantity of total RNA and other parameters need to be controlled. Most studies do this by normalising their results to a single reference gene. This study quantified the mRNA expression of three putative reference genes (ubiquitin C, cyclophilin E, and porphobilinogen deaminase) and the target gene hepatocyte growth factor receptor (HGFR) in matched colorectal tumour and normal mucosa samples. Each of the putative reference genes was found to be significantly over-expressed in the tumour samples compared to the normal samples. When HGFR expression was normalised to each of these reference genes using the 2 (-DeltaDeltaC(T)) method of relative quantification, the number of tumour samples in which HGFR was found to be over-expressed varied from 30% to 63% depending on the reference gene chosen for normalisation. This shows that normalising to a single reference gene without prior validation is inappropriate.     

One-hit stochastic decline in a mechanochemical model of cytoskeleton-induced neuron death II: transition state metastability

This is the second of two papers in which we study a mathematical model of cytoskeleton-induced neuron death. Recent evidence indicates that aggravated assembly or destruction of the cytoskeleton can trigger programmed death in neurons, by mechanisms as yet poorly understood. In our model, assembly control of the neuronal cytoskeleton interacts with both cellular stress levels and cytosolic free radical concentrations to trigger neurodegeneration. This trigger mechanism is further modulated by a diffusible toxic factor released from dying neurons. In the companion report we established that the model relates the observed general patterns of neuron decline to specific scales of cytoskeleton reorganization and cell-cell interaction strength. In this paper we study the transit of neurons through states intermediate between initial viability and cell death in our model. We find that the stochastic flow of neuron fate, from viability to cell death, self-organizes into two distinct temporal phases. There is a rapid relaxation of the initial neuron population to a more disordered phase that is long-lived, or metastable, with respect to the time scales of change in single cells. Strikingly, cellular egress from this metastable phase follows the one-hit kinetic pattern of exponential decline now established as a principal hallmark of cell death in neurodegenerative disorders. Intermediate state metastability may therefore be an important element in the systems biology of one-hit neurodegeneration.     

Three-person game facilitates indirect reciprocity under image scoring

Reputation building plays an important role in the evolution of reciprocal altruism when the same individuals do not interact repeatedly because, by referring to reputation, a reciprocator can know which partners are cooperative and can reciprocate with a cooperator. This reciprocity based on reputation is called indirect reciprocity. Previous studies of indirect reciprocity have focused only on two-person games in which only two individuals participate in a single interaction, and have claimed that indirectly reciprocal cooperation cannot be established under image scoring reputation criterion where the reputation of an individual who has cooperated (defected) becomes good (bad). In this study, we specifically examine three-person games, and reveal that indirectly reciprocal cooperation can be formed and maintained stably, even under image scoring, by a nucleus shield mechanism. In the nucleus shield, reciprocators are a shield that keeps out unconditional defectors, whereas unconditional cooperators are the backbone of cooperation that retains a good reputation among the population.     

A multilevel approach to cancer growth modeling

Cancer growth models may be divided into macroscopic models, which describe the tumor as a single entity, and microscopic ones, which consider the tumor as a complex system whose behavior emerges from the local dynamics of its basic components, the neoplastic cells. Mesoscopic models (e.g. as based on the Local Interaction Simulation Approach [Delsanto, P.P., Mignogna, R., Scalerandi, M., Schechter, R., 1998. In: Delsanto, P.P. Saenz, A.W. (Eds.), New Perspectives on Problems in Classical and Quantum Physics, vol. 2. Gordon & Breach, New Delhi, p. 5174]), which explicitly consider the behavior of cell clusters and their interactions, may be used instead of the microscopic ones, in order to study the properties of cancer biology that strongly depend on the interactions of small groups of cells at intermediate spatial and temporal scales. All these approaches have been developed independently, which limits their usefulness, since they all include relevant features and information that should be cross-correlated for a deeper understanding of the mechanisms involved. In this contribution we consider multicellular tumor spheroids as biological reference systems and propose an intermediate model to bridge the gap between a macroscopic formulation of tumor growth and a mesoscopic one. Thus we are able to establish, as an important result of our formalism, a direct correspondence between parameters characterizing processes occurring at different scales. In particular, we analyze their dependence on an important limiting factor to tumor growth, i.e. the extra-cellular matrix pressure. Since the macro and meso-models stem from totally different roots (energy conservation and clinical observations vs. cell groups dynamics), their consistency may be used to validate both approaches. It may also be interesting to note that the proposed formalism fits well into a recently proposed conjecture of growth laws universality.     

Patterning, prestress, and peeling dynamics of myocytes

As typical anchorage-dependent cells myocytes must balance contractility against adequate adhesion. Skeletal myotubes grown as isolated strips from myoblasts on micropatterned glass exhibited spontaneous peeling after one end of the myotube was mechanically detached. Such results indicate the development of a prestress in the cells. To assess this prestress and study the dynamic adhesion strength of single myocytes, the shear stress of fluid aspirated into a large-bore micropipette was then used to forcibly peel myotubes. The velocity at which cells peeled from the surface, V(peel), was measured as a continuously increasing function of the imposed tension, T(peel), which ranges from approximately 0 to 50 nN/ micro m. For each cell, peeling proved highly heterogeneous, with V(peel) fluctuating between 0 micro m/s ( approximately 80% of time) and approximately 10 micro m/s. Parallel studies of smooth muscle cells expressing GFP-paxillin also exhibited a discontinuous peeling in which focal adhesions fractured above sites of strong attachment (when pressure peeled using a small-bore pipette). The peeling approaches described here lend insight into the contractile-adhesion balance and can be used to study the real-time dynamics of stressed adhesions through both physical detection and the use of GFP markers; the methods should prove useful in comparing normal versus dystrophic muscle cells.     

Perturbation avalanches and criticality in gene regulatory networks

Boolean networks are simplified models of gene regulatory networks. We derive an approximation of the size distribution of perturbation avalanches in Boolean networks based on known results in the theory of branching processes. We show numerically that the approximation works well for different kinds of Boolean networks. It has been suggested that gene regulatory networks may be dynamically critical. To study this, as an application of the presented theory we present a novel method for estimating an order parameter from microarray data. According to the available data and our method, we find that gene regulatory networks appear to be stable and reside near the phase transition between order and chaos.     

Computation of threshold conditions for epidemiological models and global stability of the disease-free equilibrium (DFE)

One goal of this paper is to give an algorithm for computing a threshold condition for epidemiological systems arising from compartmental deterministic modeling. We calculate a threshold condition T(0) of the parameters of the system such that if T(0)<1 the disease-free equilibrium (DFE) is locally asymptotically stable (LAS), and if T(0)>1, the DFE is unstable. The second objective, by adding some reasonable assumptions, is to give, depending on the model, necessary and sufficient conditions for global asymptotic stability (GAS) of the DFE. In many cases, we can prove that a necessary and sufficient condition for the global asymptotic stability of the DFE is R(0)< or =1, where R(0) is the basic reproduction number [O. Diekmann, J.A. Heesterbeek, Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation, Wiley, New York, 2000]. To illustrate our results, we apply our techniques to examples taken from the literature. In these examples we improve the results already obtained for the GAS of the DFE. We show that our algorithm is relevant for high dimensional epidemiological models.     

The evolution of n-player cooperation-threshold games and ESS bifurcations

An evolutionary game of individuals cooperating to obtain a collective benefit is here modelled as an n-player Prisoner's Dilemma game. With reference to biological situations, such as group foraging, we introduce a threshold condition in the number of cooperators required to obtain the collective benefit. In the simplest version, a three-player game, complex behaviour appears as the replicator dynamics exhibits a catastrophic event separating a parameter region allowing for coexistence of cooperators and defectors and a region of pure defection. Cooperation emerges through an ESS bifurcation, and cooperators only thrive beyond a critical point in cost-benefit space. Moreover, a repelling fixed point of the dynamics acts as a barrier to the introduction of cooperation in defecting populations. The results illustrate the qualitative difference between two-player games and multiple player games and thus the limitations to the generality of conclusions from two-player games. We present a procedure to find the evolutionarily stable strategies in any n-player game with cost and benefit depending on the number of cooperators. This was previously done by Motro [1991. Co-operation and defection: playing the field and the ESS. J. Theor. Biol. 151, 145-154] in the special cases of convex and concave benefit functions and constant cost.     

The role of low avidity T cells in the protection against type 1 diabetes: a modeling investigation

Cytotoxic T lymphocytes (CTLs) play a dominant role in the pathogenesis of autoimmune diabetes, commonly denoted Type 1 Diabetes (T1D). These CTLs (notably CD8⁺ T cells) recognize and kill insulin-secreting pancreatic β cells, reducing their number by ∼90%. The resulting reduction of insulin secretion causes the defective regulation of glucose metabolism, leading to the characteristic symptoms of diabetes. Recognition of β cells as targets by CTLs depends on the interactions between MHC-peptide complexes on the surface of β cells and receptors (TCRs) on T cells. Those CTLs with high affinity TCRs (also called high avidity T cells) cause most of the harm, while those with low affinity TCRs (also called low avidity T cells) play a more mysterious role. Recent experimental evidence suggests that low avidity T cells accumulate as memory T cells during the disease and may be protective in NOD mice (a strain prone to developing T1D), delaying disease progression. It has been hypothesized that such low avidity T cells afford disease protection either by crowding the islets of Langerhans, where β cells reside, or by killing antigen presenting cells (APCs).In this paper, we explore the hypothesized mechanisms for this protective effect in the context of a series of models for (1) the interactions of low and high avidity T cells, (2) the effect of APCs and (3) the feedback from β cell killing to autoantigen-induced T cell proliferation. We analyze properties of these models, noting consistency of predictions with observed behaviour. We then use the models to examine the influence of various treatment strategies on the progression of the disease. The model reveals that progressive accumulation of memory low avidity autoreactive T cells during disease progression makes treatments aimed at expanding these protective T cell types more effective close to, or at the onset of clinical disease. It also provides evidence for the hypothesis that low avidity T cells kill APCs (rather than the alternate hypothesis that they crowd the islets).     

Molecular cloning and functional expression of a sodium bicarbonate cotransporter from guinea-pig parotid glands

We recently found that the concentration of HCO3- in guinea-pig saliva is very similar to that of human saliva; however, the entity that regulates HCO3- transport has not yet been fully characterized. In order to investigate the mechanism of HCO3- transport, we identified, cloned, and characterized a sodium bicarbonate (Na(+)/HCO3- cotransporter found in guinea-pig parotid glands (gpNBC1). The gpNBC1 gene encodes a 1079-amino acid protein that has 95% and 96% homology with human and mouse parotid NBC1, respectively. Oocytes expressing gpNBC1 were exposed to HCO3- or Na(+)-free solutions, which resulted in a marked change in membrane potentials (V(m)), suggesting that gpNBC1 is electrogenic. Likewise, a gpNBC1-mediated pH recovery was observed in gpNBC1 transfected human hepatoma cells; however, in the presence of 4, 4-diisothiocyanostilbene-2,2-disulfonic acid, a specific NBC1 inhibitor, such changes in V(m) and pH(i) were not observed. Together, the data show that the cloned guinea-pig gene is a functional, as well as sequence homologue of human NBC1.     

Changes of enzyme activity in lipid signaling pathways related to substrate reordering

The static fluid mosaic model of biological membranes has been progressively complemented by a dynamic membrane model that includes phospholipid reordering in domains that are proposed to extend from nanometers to microns. Kinetic models for lipolytic enzymes have only been developed for homogeneous lipid phases. In this work, we develop a generalization of the well-known surface dilution kinetic theory to cases where, in a same lipid phase, both domain and nondomain phases coexist. Our model also allows understanding the changes in enzymatic activity due to a decrease of free substrate concentration when domains are induced by peptides. This lipid reordering and domain dynamics can affect the activity of lipolytic enzymes, and can provide a simple explanation for how basic peptides, with a strong direct interaction with acidic phospholipids (such as beta-amyloid peptide), may cause a complex modulation of the activities of many important enzymes in lipid signaling pathways.     

Functional expression and characterization of Aedes aegypti dopachrome conversion enzyme.

A full-length mosquito dopachrome conversion enzyme (DCE) and its truncated form lacking the last 54 carboxyl-terminal amino acid residues are expressed using a baculovirus/insect cell expression system. The full-length recombinant DCE displayed multiple bands during native PAGE with substrate staining, but only one active band was detected when the truncated recombinant DCE was analyzed under identical analysis conditions. Our data suggest that the last 50 some carboxyl-terminal residues are involved in the polymerization of the DCE molecules and that the proposed DCE isozymes likely reflect the presence of multimers of the same DCE molecules. The significance of the recombinant DCE in accelerating the melanization pathway is demonstrated by a rapid production of melanin in a dopa and tyrosinase reaction mixture in the presence of recombinant DCE. The DCE sequence data obtained in our previous study, together with results of functional expression and biochemical characterization achieved in this study, provide a necessary reference for the study of other insect DCEs.     

Pressure equilibrium and jump study on unfolding of 23-kDa protein from spinach photosystem II

Pressure-induced unfolding of 23-kDa protein from spinach photosystem II has been systematically investigated at various experimental conditions. Thermodynamic equilibrium studies indicate that the protein is very sensitive to pressure. At 20 degrees C and pH 5.5, 23-kDa protein shows a reversible two-state unfolding transition under pressure with a midpoint near 160 MPa, which is much lower than most natural proteins studied to date. The free energy (DeltaG(u)) and volume change (DeltaV(u)) for the unfolding are 5.9 kcal/mol and -160 ml/mol, respectively. It was found that NaCl and sucrose significantly stabilize the protein from unfolding and the stabilization is associated not only with an increase in DeltaG(u) but also with a decrease in DeltaV(u). The pressure-jump studies of 23-kDa protein reveal a negative activation volume for unfolding (-66.2 ml/mol) and a positive activation volume for refolding (84.1 ml/mol), indicating that, in terms of system volume, the protein transition state lies between the folded and unfolded states. Examination of the temperature effect on the unfolding kinetics indicates that the thermal expansibility of the transition state and the unfolded state of 23-kDa protein are closer to each other and they are larger than that of the native state. The diverse pressure-refolding pathways of 23-kDa protein in some conditions were revealed in pressure-jump kinetics.