Combinatorial complexity and dynamical restriction of network flows in signal transduction

The activities and interactions of proteins that govern the cellular response to a signal generate a multitude of protein phosphorylation states and heterogeneous protein complexes. Here, using a computational model that accounts for 307 molecular species implied by specified interactions of four proteins involved in signalling by the immunoreceptor FcepsilonRI, we determine the relative importance of molecular species that can be generated during signalling, chemical transitions among these species, and reaction paths that lead to activation of the protein tyrosine kinase (PTK) Syk. By all of these measures and over two- and ten-fold ranges of model parameters--rate constants and initial concentrations--only a small portion of the biochemical network is active. The spectrum of active complexes, however, can be shifted dramatically, even by a change in the concentration of a single protein, which suggests that the network can produce qualitatively different responses under different cellular conditions and in response to different inputs. Reduced models that reproduce predictions of the full model for a particular set of parameters lose their predictive capacity when parameters are varied over two-fold ranges.     

Quantifying aggregation of IgE-FcepsilonRI by multivalent antigen

Aggregation of cell surface receptors by multivalent ligand can trigger a variety of cellular responses. A well-studied receptor that responds to aggregation is the high affinity receptor for IgE (FcepsilonRI), which is responsible for initiating allergic reactions. To quantify antigen-induced aggregation of IgE-FcepsilonRI complexes, we have developed a method based on multiparameter flow cytometry to monitor both occupancy of surface IgE combining sites and association of antigen with the cell surface. The number of bound IgE combining sites in excess of the number of bound antigens, the number of bridges between receptors, provides a quantitative measure of IgE-FcepsilonRI aggregation. We demonstrate our method by using it to study the equilibrium binding of a haptenated fluorescent protein, 2,4-dinitrophenol-coupled B-phycoerythrin (DNP25-PE), to fluorescein isothiocyanate-labeled anti-DNP IgE on the surface of rat basophilic leukemia cells. The results, which we analyze with the aid of a mathematical model, indicate how IgE-FcepsilonRI aggregation depends on the total concentrations of DNP25-PE and surface IgE. As expected, we find that maximal aggregation occurs at an optimal antigen concentration. We also find that aggregation varies qualitatively with the total concentration of surface IgE as predicted by an earlier theoretical analysis.     

Concentration-dependent diffusion coefficient and dissipative structures

A theoretical study of the Brusselator model with non-uniform distribution of component A and a concentration-dependent diffusion coefficient has been performed. Numerical simulation reveals that a variable diffusion coefficient alters the bifurcation pattern and the stability properties of the steady-state as well as periodic solutions. A simple approximate method, based on one-point collocation, has been proposed to analyze the bifurcation phenomena for the case of fixed boundary conditions and low system size.     

Scaffold-mediated nucleation of protein signaling complexes: elementary principles

Proteins with multiple binding sites play important roles in cell signaling systems by nucleating protein complexes in which, for example, enzymes and substrates are co-localized. Proteins that specialize in this function are called by a variety names, including adapter, linker and scaffold. Scaffold-mediated nucleation of protein complexes can be either constitutive or induced. Induced nucleation is commonly mediated by a docking site on a scaffold that is activated by phosphorylation. Here, by considering minimalist mathematical models, which recapitulate scaffold effects seen in more mechanistically detailed models, we obtain analytical and numerical results that provide insights into scaffold function. These results elucidate how recruitment of a pair of ligands to a scaffold depends on the concentrations of the ligands, on the binding constants for ligand-scaffold interactions, on binding cooperativity, and on the milieu of the scaffold, as ligand recruitment is affected by competitive ligands and decoy receptors. For the case of a bivalent scaffold, we obtain an expression for the unique scaffold concentration that maximally recruits a pair of monovalent ligands. Through simulations, we demonstrate that a bivalent scaffold can nucleate distinct sets of ligands to equivalent extents when the scaffold is present at different concentrations. Thus, the function of a scaffold can potentially change qualitatively with a change in copy number. We also demonstrate how a scaffold can change the catalytic efficiency of an enzyme and the sensitivity of the rate of reaction to substrate concentration. The results presented here should be useful for understanding scaffold function and for engineering scaffolds to have desired properties.     

Rmd9p controls the processing/stability of mitochondrial mRNAs and its overexpression compensates for a partial deficiency of oxa1p in Saccharomyces cerevisiae

Oxa1p is a key component of the general membrane insertion machinery of eukaryotic respiratory complex subunits encoded by the mitochondrial genome. In this study, we have generated a respiratory-deficient mutant, oxa1-E65G-F229S, that contains two substitutions in the predicted intermembrane space domain of Oxa1p. The respiratory deficiency due to this mutation is compensated for by overexpressing RMD9. We show that Rmd9p is an extrinsic membrane protein facing the matrix side of the mitochondrial inner membrane. Its deletion leads to a pleiotropic effect on respiratory complex biogenesis. The steady-state level of all the mitochondrial mRNAs encoding respiratory complex subunits is strongly reduced in the Deltarmd9 mutant, and there is a slight decrease in the accumulation of two RNAs encoding components of the small subunit of the mitochondrial ribosome. Overexpressing RMD9 leads to an increase in the steady-state level of mitochondrial RNAs, and we discuss how this increase could suppress the oxa1 mutations and compensate for the membrane insertion defect of the subunits encoded by these mRNAs.     

Effects of conformational constraint in 2- and 8-cycloleucine analogues of oxytocin and [1-penicillamine] oxytocin examined by circular dichroism and biosassay

The analogues of oxytocin and [1-penicillamine]oxytocin, containing a cycloleucine (Cle) residue in position 2 or 8, were investigated by means of circular dichroism measurements in different solvents, and the results examined in terms of their biological activities. A cycloleucine residue in position 2 substantially reduces the free conformational space of the hormone 20-membered ring moiety (including the disulfide group), and stabilizes a conformation which is close to one of the possible conformations of oxytocin and involves a -turn. In position 8, the Cle residue affects the conformation of the Tyr² side chain, apparently forcing it away from the space above the 20-membered disulfide ring. However, it does not appear that the Cle residue has any significant effect on the overall backbone conformation of the hormone. The steric effect of the penicillamine residue in position 1 on the conformation of the disulfide group and Tyr² side chain from previous investigations is further confirmed. The synthesis and biological potency of [1-penicillamine, 8-cycloleucine]oxytocin is described. This analogue exhibits a strong inhibitory effect on the uterotonic activity of oxytocinin vitro. It also inhibited the vasopressor response to vasopressin.     

Stochastic effects and bistability in T cell receptor signaling

The stochastic dynamics of T cell receptor (TCR) signaling are studied using a mathematical model intended to capture kinetic proofreading (sensitivity to ligand-receptor binding kinetics) and negative and positive feedback regulation mediated, respectively, by the phosphatase SHP1 and the MAP kinase ERK. The model incorporates protein-protein interactions involved in initiating TCR-mediated cellular responses and reproduces several experimental observations about the behavior of TCR signaling, including robust responses to as few as a handful of ligands (agonist peptide-MHC complexes on an antigen-presenting cell), distinct responses to ligands that bind TCR with different lifetimes, and antagonism. Analysis of the model indicates that TCR signaling dynamics are marked by significant stochastic fluctuations and bistability, which is caused by the competition between the positive and negative feedbacks. Stochastic fluctuations are such that single-cell trajectories differ qualitatively from the trajectory predicted in the deterministic approximation of the dynamics. Because of bistability, the average of single-cell trajectories differs markedly from the deterministic trajectory. Bistability combined with stochastic fluctuations allows for switch-like responses to signals, which may aid T cells in making committed cell-fate decisions.