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.     

The shared reward dilemma

One of the most direct human mechanisms of promoting cooperation is rewarding it. We study the effect of sharing a reward among cooperators in the most stringent form of social dilemma, namely the prisoner's dilemma (PD). Specifically, for a group of players that collect payoffs by playing a pairwise PD game with their partners, we consider an external entity that distributes a fixed reward equally among all cooperators. Thus, individuals confront a new dilemma: on the one hand, they may be inclined to choose the shared reward despite the possibility of being exploited by defectors; on the other hand, if too many players do that, cooperators will obtain a poor reward and defectors will outperform them. By appropriately tuning the amount to be shared a vast variety of scenarios arises, including the traditional ones in the study of cooperation as well as more complex situations where unexpected behavior can occur. We provide a complete classification of the equilibria of the n-player game as well as of its evolutionary dynamics.     

Evolutionary stability of first-order-information indirect reciprocity in sizable groups

Indirect reciprocity is considered as a key mechanism for explaining the evolution of cooperation in situations where the same individuals interact only a few times. Under indirect reciprocity, an individual who helps others gets returns indirectly from others who know her good reputation. Recently, many studies have discussed the effect of reputation criteria based only on the former actions of the others (first-order information) and of those based also on the former reputation of opponents of the others (second-order information) on the evolution of indirect reciprocity. In this study, we investigate the evolutionary stability of the indirectly reciprocal strategy (discriminating strategy: DIS), which cooperates only with opponents who have good reputations, in -person games where more than two individuals take part in a single group (interaction). We show that in n-person games, DIS is an evolutionarily stable strategy (ESS) even under the image-scoring reputation criterion, which is based only on first-order information and where cooperations (defections) are judged to be good (bad). This result is in contrast to that of 2-person games where DIS is not an ESS under reputation criteria based only on first-order information.     

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.     

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.     

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.     

Uncovering operational interactions in genetic networks using asynchronous Boolean dynamics

Biological networks of large dimensions, with their diagram of interactions, are often well represented by a Boolean model with a family of logical rules. The state space of a Boolean model is finite, and its asynchronous dynamics are fully described by a transition graph in the state space. In this context, a model reduction method will be developed for identifying the active or operational interactions responsible for a given dynamic behaviour. The first step in this procedure is the decomposition of the asynchronous transition graph into its strongly connected components, to obtain a “reduced” and hierarchically organized graph of transitions. The second step consists of the identification of a partial graph of interactions and a sub-family of logical rules that remain operational in a given region of the state space. This model reduction method and its usefulness are illustrated by an application to a model of programmed cell death. The method identifies two mechanisms used by the cell to respond to death-receptor stimulation and decide between the survival and apoptotic pathways.     

KCNE3 is an inhibitory subunit of the Kv4.3 potassium channel

The mammalian Kv4.3 potassium channel is a fast activating and inactivating K+ channel widely distributed in mammalian tissues. Kv4.3 is the major component of various physiologically important currents ranging from A-type currents in the CNS to the transient outward potassium conductance in the heart (I(to)). Here we show that the KCNE3 beta-subunit has a strong inhibitory effect on current conducted by heterologously expressed Kv4.3 channels. KCNE3 reduces the Kv4.3 current amplitude, and it slows down the channel activation and inactivation as well as the recovery from inactivation. KCNE3 also inhibits currents generated by Kv4.3 in complex with the accessory subunit KChIP2. We find the inhibitory effect of KCNE3 to be specific for Kv4.3 within the Kv4 channel family. Kv4.3 has previously been shown to interact with a number of beta-subunits, but none of the described subunit-interactions exert an inhibitory effect on the Kv4.3 current.     

Melanoblast proliferation dynamics during mouse embryonic development. Modeling and validation

In this paper, we are looking for mathematical modeling of mouse embryonic melanoblast proliferation dynamics, taking into account, the expression level of β‐catenin. This protein plays an important role into the whole signal pathway process. Different assumptions on some unobservable features lead to different candidate models. From real data measured, from biological experiments and from a priori biological knowledge, it was able to validate or invalidate some of the candidate models. Data assimilation and parameter identification allowed us to derive a mathematical model that is in very good agreement with biological data. As a result, the produced model can give tracks for biologists into their biological investigations and experimental evidence. Another interest is the use of this model for robust hidden parameter identification like double times or number of founder melanoblasts.     

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.     

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.     

Application of HN(C)N to rapid estimation of 1J(N-C(alpha)) coupling constants correlated to psi torsion angles in proteins: implication to structural genomics

We recently described a triple resonance experiment, HN(C)N, for sequential correlation of H(N) and 15N atoms in (15N, 13C) labeled proteins [J. Biomol. NMR. 20 (2001) 135]. Here, we describe an approach based on this experiment for estimation of one bond N-C(alpha) J-couplings in medium size labeled proteins, which seem to show good correlations with psi torsion angles along the protein backbone. The approach uses the ratio of the intensities of the sequential and diagonal peaks in the F(2)-F(3) planes of the HN(C)N spectrum. The reliability of the approach has been demonstrated using a short peptide wherein the coupling constants have been measured by the present method and also independently from peak splittings in HSQC spectra. The two results agree within 10%. The applicability of the procedure to proteins has been demonstrated using doubly labeled FK506 binding protein (FKBP, molecular mass approximately 12 kDa). Coupling constant estimates have been obtained for 62 out of 100 non-proline residues and they show a correlation with psi torsion angles, as has been reported before. This semi-quantitative application of HN(C)N extends the significance of the experiment especially, in the context of structural genomics, since the single experiment, not only provides a great enhancement in the speed of resonance assignment, but also provides quantitative structural information.     

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.     

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.     

Crystal structure of oxidized cytochrome c(6A) from Arabidopsis thaliana

Compared with algal and cyanobacterial cytochrome c(6), cytochrome c(6A) from higher plants contains an additional loop of 12 amino acid residues. We have determined the first crystal structure of cytochrome c(6A) from Arabidopsis thaliana at 1.5 Angstrom resolution in order to help elucidate its function. The overall structure of cytochrome c(6A) follows the topology of class I c-type cytochromes in which the heme prosthetic group covalently binds to Cys16 and Cys19, and the iron has octahedral coordination with His20 and Met60 as the axial ligands. Two cysteine residues (Cys67 and Cys73) within the characteristic 12 amino acids loop form a disulfide bond, contributing to the structural stability of cytochrome c(6A). Our model provides a chemical basis for the known low redox potential of cytochrome c(6A) which makes it an unsuitable electron carrier between cytochrome b(6)f and PSI.     

A divergent synthesis of [1-14C]-mono-E isomers of fatty acids

A convenient synthesis of [1-14C]-mono-trans fatty acid using olefin inversion as a key-step is described. This methodology allows for a facile synthesis of [1-14C]-labelled mono-trans analogues of oleic, linoleic and linolenic acids. As an example, only eleven steps were necessary to obtain the [1-14C]-mono-E isomers of linolenic acid from its commercial all-Z form. In the first step, Barton's decarboxylation procedure yielded a bromo intermediate. Epoxidation of this compound resulted in the formation of three monoepoxides, which could be separated by HPLC. After identification by 1H NMR and MS, the pure monoepoxides were then subjected to inversion consisting of a stereospecific deoxygenation followed by a beta-elimination step. Finally, the labelling was introduced by substitution of the bromine by a [14C]-cyano group followed by hydrolysis.     

Modeling species-specific diacylglycerol dynamics in the RAW 264.7 macrophage

A mathematical model of the G protein signaling pathway in RAW 264.7 macrophages downstream of P2Y₆ receptors activated by the ubiquitous signaling nucleotide uridine 5’-diphosphate is developed. The model, which is based on time-course measurements of inositol trisphosphate, cytosolic calcium, and diacylglycerol, focuses particularly on differential dynamics of multiple chemical species of diacylglycerol. When using the canonical pathway representation, the model predicted that key interactions were missing from the current network structure. Indeed, the model suggested that accurate depiction of experimental observations required an additional branch to the signaling pathway. An intracellular pool of diacylglycerol is immediately phosphorylated upon stimulation of an extracellular receptor for uridine 5’-diphosphate and subsequently used to aid replenishment of phosphatidylinositol. As a result of sensitivity analysis of the model parameters, key predictions can be made regarding which of these parameters are the most sensitive to perturbations and are therefore most responsible for output uncertainty.     

Identification of an actin-binding site in p47phox an organizer protein of NADPH oxidase

Actin has been reported to enhance the superoxide-generating activity of neutrophil NADPH oxidase in a cell-free system and to interact with p47phox, a regulatory subunit of the oxidase. In the present study, we searched for an actin-binding site in p47phox by far-western blotting and blot-binding assays using truncated forms of p47phox. The amino-acid sequence 319-337 was identified as an actin-binding site, and a synthetic peptide of this sequence bound to actin. The sequence shows no homology to other known actin-binding motifs. It is located in the autoinhibitory region of p47phox and includes Ser-328, a phosphorylation site essential for unmasking. Although a phosphorylation-mimetic p47phox mutant bound to actin with a lower affinity than the wild type, the same mutant interacted with filamentous actin more efficiently than the wild type. A mutant peptide p47phox (319-337, Ser328Glu) bound to filamentous actin more tightly than to monomer actin. These results suggest that p47phox moves to cortical actin when it becomes unmasked in the cells.