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.     

Molecular characterization of β1,4-galactosyltransferase 7 genetic mutations linked to the progeroid form of Ehlers-Danlos syndrome (EDS)

β1,4-Galactosyltransferase 7 (β4GalT7) is a key enzyme initiating glycosaminoglycan (GAG) synthesis. Based on in vitro and ex vivo kinetics studies and structure-based modelling, we molecularly characterized β4GalT7 mutants linked to the progeroid form of Ehlers-Danlos syndrome (EDS), a severe connective tissue disorder. Our results revealed that loss of activity upon L206P substitution due to altered protein folding is the primary cause for the GAG synthesis defect in patients carrying the compound A186D and L206P mutations. We showed that R270C substitution strongly reduced β4GalT7 affinity towards xyloside acceptor, thus affecting GAG chains formation. This study establishes the molecular basis for β4GalT7 defects associated with altered GAG synthesis in EDS.     

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.     

The N-terminal rhodanese domain from Azotobacter vinelandii has a stable and folded structure independently of the C-terminal domain

Sulfurtransferase are enzymes involved in the formation, conversion and transport of compounds containing sulfane-sulfur atoms. Although the three-dimensional structure of the rhodanese from the nitrogen-fixing bacterium Azotobacter vinelandii is known, the role of its two domains in the protein conformational stability is still obscure. We have evaluated the susceptibility to proteolytic degradation of the two domains of the enzyme. The two domains show different resistance to the endoproteinases and, in particular, the N-terminal domain shows to be more stable to digestion during time than the C-terminal one. Cloning and overexpression of the N-terminal domain of the protein was performed to better understand its functional and structural role. The recombinant N-terminal domain of rhodanese A. vinelandii is soluble in water solution and the spectroscopic studies by circular dichroism and heteronuclear NMR spectroscopy indicate a stable fold of the protein with the expected alpha/beta topology. The results indicate that this N-terminal domain has already got all the elements necessary for an C-terminal domain independent folding. Its solution structure by NMR, actually under course, will be a valid contribution to understand the role of this domain in the folding process of the sulfurtransferase.     

Self-organization at the origin of life

The concept of an (M,R) system with organizational invariance allows one to understand how a system may be able to maintain itself indefinitely if it is coupled to an external source of energy and materials. However, although this constitutes an important step towards understanding the difference between a living and a non-living system, it is not clear that an (M,R) system with organizational invariance is sufficient to define a living system. To take a further step towards defining what it means to be alive it is necessary to add to a simple (M,R) system some property that represents its identity, and which can be maintained and modified in subsequent generations.     

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.     

Conformational properties of the aggregation precursor state of HypF-N

The conversion of specific proteins or protein fragments into insoluble, ordered fibrillar aggregates is a fundamental process in protein chemistry, biology, medicine and biotechnology. As this structural conversion seems to be a property shared by many proteins, understanding the mechanism of this process will be of extreme importance. Here we present a structural characterisation of a conformational state populated at low pH by the N-terminal domain of Escherichia coli HypF. Combining different biophysical and biochemical techniques, including near- and far-UV circular dichroism, intrinsic and 8-anilinonaphthalene-1-sulfonate-derived fluorescence, dynamic light scattering and limited proteolysis, we will show that this state is largely unfolded but contains significant secondary structure and hydrophobic clusters. It also appears to be more compact than a random coil-like state but less organised than a molten globule state. Increase of the total ionic strength of the solution induces aggregation of such a pre-molten globule state into amyloid-like protofibrils, as revealed by thioflavin T fluorescence and atomic force microscopy. These results show that a pre-molten globule state can be, among other possible conformational states, one of the precursor states of amyloid formation. In addition, the possibility of triggering aggregation by modulating the ionic strength of the solution provides one a unique opportunity to study both the initial precursor state and the aggregation process.     

Plasticity within the obligatory folding nucleus of an immunoglobulin-like domain

A number of β-sandwich immunoglobulin-like domains have been shown to fold using a set of structurally equivalent residues that form a folding nucleus deep within the core of the protein. Formation of this nucleus is sufficient to establish the complex Greek key topology of the native state. These nucleating residues are highly conserved within the immunoglobulin superfamily, but are less well conserved in the fibronectin type III (fnIII) superfamily, where the requirement is simply to have four interacting hydrophobic residues. However, there are rare examples where this nucleation pattern is absent. In this study, we have investigated the folding of a novel member of the fnIII superfamily whose nucleus appears to lack one of the four buried hydrophobic residues. We show that the folding mechanism is unaltered, but the folding nucleus has moved within the hydrophobic core.     

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.     

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).     

Genetic diversity of microsatellite loci in hierarchically structured populations

Microsatellite loci are widely used for investigating patterns of genetic variation within and among populations. Those patterns are in turn determined by population sizes, migration rates, and mutation rates. We provide exact expressions for the first two moments of the allele frequency distribution in a stochastic model appropriate for studying microsatellite evolution with migration, mutation, and drift under the assumption that the range of allele sizes is bounded. Using these results, we study the behavior of several measures related to Wright’s F_ST, including Slatkin’s R_ST. Our analytical approximations for F_ST and R_ST show that familiar relationships between N_em and F_ST or R_ST hold when the migration and mutation rates are small. Using the exact expressions for F_ST and R_ST, our numerical results show that, when the migration and mutation rates are large, these relationships no longer hold. Our numerical results also show that the diversity measures most closely related to F_ST depend on mutation rates, mutational models (stepwise versus two-phase), migration rates, and population sizes. Surprisingly, R_ST is relatively insensitive to the mutation rates and mutational models. The differing behaviors of R_ST and F_ST suggest that properties of the among-population distribution of allele frequencies may allow the roles of mutation and migration in producing patterns of diversity to be distinguished, a topic of continuing investigation.     

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.     

The crystal structure of BamB suggests interactions with BamA and its role within the BAM complex

Escherichia coli BamB is the largest of four lipoproteins in the β-barrel assembly machinery (BAM) complex. It interacts with the periplasmic domain of BamA, an integral outer membrane protein (OMP) essential for OMP biogenesis. Although BamB is not essential, it serves an important function in the BAM complex, significantly increasing the folding efficiency of some OMPs in vivo and in vitro. To learn more about the BAM complex, we solved structures of BamB in three different crystal forms. BamB crystallized in space groups P2₁3, I222, and P2₁2₁2₁, with one molecule per asymmetric unit in each case. Crystals from the space group I222 diffracted to 1. 65-Å resolution. BamB forms an eight-bladed β-propeller with a central pore and is shaped like a doughnut. A DALI search revealed that BamB shares structural homology to several eukaryotic proteins containing WD40 repeat domains, which commonly have β-propeller folds and often serve as scaffolding proteins within larger multi-protein complexes that carry out signal transduction, cell division, and chemotaxis. Using mutagenesis data from previous studies, we docked BamB onto a BamA structural model and assessed known and possible interactions between these two proteins. Our data suggest that BamB serves as a scaffolding protein within the BAM complex by optimally orienting the flexible periplasmic domain of BamA for interaction with other BAM components and chaperones. This may facilitate integration of newly synthesized OMPs into the outer membrane.     

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.     

Sterol chemical configuration influences the thermotropic phase behaviour of dipalmitoylphosphatidylcholine bilayers containing 5α-cholestan-3β- and 3α-ol

It is commonly believed that all membrane sterols are rigid all-trans ring systems with a fully extended alkyl side-chain and that they similarly influence phospholipid bilayer physical properties. Here, we report the sterol concentration-dependent, thermotropic phase behaviour of binary dipalmitoylphosphatidylcholine (DPPC)/sterol mixtures containing two similar 5α-H sterols with different functional group orientations (3α-OH or 3β-OH), which adopt an ideal all-trans planar ring conformation but lack the deformed ring B conformation of cholesterol (Chol) and epicholesterol (Echol), using differential scanning calorimetry (DSC). Our deconvolution of the DSC main phase transition endotherms show differences in the proportions of sterol-poor (sharp) and sterol-rich (broad) domains in the DPPC bilayer with increasing sterol concentration, which delineate gel/liquid-crystalline (P_β′/L_α) and disordered gel (L_β)/liquid-ordered (l_o) phase regions. There are similarities in the DPPC main phase transition temperature, cooperativity and enthalpy for each 3β-ol and 3α-ol pair with increasing sterol concentration and differences in the parameters obtained for both the sterol-poor and sterol-rich regions. The sterol-poor domain persists over a greater concentration range in both 3α-ol/DPPC mixtures, suggesting that either those domains are more stable in the 3α-ols or that those sterols are less miscible in the sterol-rich domain. Corresponding parameters for the sterol-rich domain show that at sterol concentrations up to 20 mol%, the 5α-H,3β-ol is more effective at reducing the phase transition enthalpy of the broad component () than Chol, but is less effective at higher concentrations. Although mixtures containing Echol and 5α-cholestan-3α-ol have similar positive slopes below 7 mol% sterol, suggesting that they abolish the L_β/l_o phase transition equally effectively at low concentrations, Echol is more effective than the saturated 3α-ol at higher sterol concentrations. A comparison of obtained for the saturated and unsaturated pairs suggests that the latter sterols stabilize the l_o phase and broaden and abolish the DPPC main phase transition more effectively than the saturated sterols at physiologically relevant concentrations, supporting the idea that the double bond of Chol and Echol promotes greater sterol miscibility and the formation of l_o phase lipid bilayers relative to corresponding saturated sterols in biological membranes.     

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.     

Mesophile versus thermophile: insights into the structural mechanisms of kinetic stability

Obtaining detailed knowledge of folding intermediate and transition state (TS) structures is critical for understanding protein folding mechanisms. Comparisons between proteins adapted to survive extreme temperatures with their mesophilic homologs are likely to provide valuable information on the interactions relevant to the unfolding transition. For kinetically stable proteins such as alpha-lytic protease (alphaLP) and its family members, their large free energy barrier to unfolding is central to their biological function. To gain new insights into the mechanisms that underlie kinetic stability, we have determined the structure and high temperature unfolding kinetics of a thermophilic homolog, Thermobifida fusca protease A (TFPA). These studies led to the identification of a specific structural element bridging the N and C-terminal domains of the protease (the "domain bridge") proposed to be associated with the enhanced high temperature kinetic stability in TFPA. Mutagenesis experiments exchanging the TFPA domain bridge into alphaLP validate this hypothesis and illustrate key structural details that contribute to TFPA's increased kinetic thermostability. These results lead to an updated model for the unfolding transition state structure for this important class of proteases in which domain bridge undocking and unfolding occurs at or before the TS. The domain bridge appears to be a structural element that can modulate the degree of kinetic stability of the different members of this class of proteases.