Computational Modeling for the Activation Cycle of G-proteins by G-protein-coupled Receptors

In this paper, we survey five different computational modeling methods. For comparison, we use the activation cycle of G-proteins that regulate cellular signaling events downstream of G-protein-coupled receptors (GPCRs) as a driving example. Starting from an existing Ordinary Differential Equations (ODEs) model, we implement the G-protein cycle in the stochastic Pi-calculus using SPiM, as Petri-nets using Cell Illustrator, in the Kappa Language using Cellucidate, and in Bio-PEPA using the Bio-PEPA eclipse plug in. We also provide a high-level notation to abstract away from communication primitives that may be unfamiliar to the average biologist, and we show how to translate high-level programs into stochastic Pi-calculus processes and chemical reactions.Comment: In Proceedings MeCBIC 2010, arXiv:1011.005

Modeling Robustness Tradeoffs in Yeast Cell Polarization Induced by Spatial Gradients

Cells localize (polarize) internal components to specific locations in response to external signals such as spatial gradients. For example, yeast cells form a mating projection toward the source of mating pheromone. There are specific challenges associated with cell polarization including amplification of shallow external gradients of ligand to produce steep internal gradients of protein components (e.g. localized distribution), response over a broad range of ligand concentrations, and tracking of moving signal sources. In this work, we investigated the tradeoffs among these performance objectives using a generic model that captures the basic spatial dynamics of polarization in yeast cells, which are small. We varied the positive feedback, cooperativity, and diffusion coefficients in the model to explore the nature of this tradeoff. Increasing the positive feedback gain resulted in better amplification, but also produced multiple steady-states and hysteresis that prevented the tracking of directional changes of the gradient. Feedforward/feedback coincidence detection in the positive feedback loop and multi-stage amplification both improved tracking with only a modest loss of amplification. Surprisingly, we found that introducing lateral surface diffusion increased the robustness of polarization and collapsed the multiple steady-states to a single steady-state at the cost of a reduction in polarization. Finally, in a more mechanistic model of yeast cell polarization, a surface diffusion coefficient between 0.01 and 0.001 µm2/s produced the best polarization performance, and this range is close to the measured value. The model also showed good gradient-sensitivity and dynamic range. This research is significant because it provides an in-depth analysis of the performance tradeoffs that confront biological systems that sense and respond to chemical spatial gradients, proposes strategies for balancing this tradeoff, highlights the critical role of lateral diffusion of proteins in the membrane on the robustness of polarization, and furnishes a framework for future spatial models of yeast cell polarization

Reverse Engineering a Signaling Network Using Alternative Inputs

One of the goals of systems biology is to reverse engineer in a comprehensive fashion the arrow diagrams of signal transduction systems. An important tool for ordering pathway components is genetic epistasis analysis, and here we present a strategy termed Alternative Inputs (AIs) to perform systematic epistasis analysis. An alternative input is defined as any genetic manipulation that can activate the signaling pathway instead of the natural input. We introduced the concept of an “AIs-Deletions matrix” that summarizes the outputs of all combinations of alternative inputs and deletions. We developed the theory and algorithms to construct a pairwise relationship graph from the AIs-Deletions matrix capturing both functional ordering (upstream, downstream) and logical relationships (AND, OR), and then interpreting these relationships into a standard arrow diagram. As a proof-of-principle, we applied this methodology to a subset of genes involved in yeast mating signaling. This experimental pilot study highlights the robustness of the approach and important technical challenges. In summary, this research formalizes and extends classical epistasis analysis from linear pathways to more complex networks, facilitating computational analysis and reconstruction of signaling arrow diagrams

Synthetic Morphology Using Alternative Inputs

Designing the shape and size of a cell is an interesting challenge for synthetic biology. Prolonged exposure to the mating pheromone α-factor induces an unusual morphology in yeast cells: multiple mating projections. The goal of this work was to reproduce the multiple projections phenotype in the absence of α-factor using a gain-of-function approach termed “Alternative Inputs (AIs)”. An alternative input is defined as any genetic manipulation that can activate the signaling pathway instead of the natural input. Interestingly, none of the alternative inputs were sufficient to produce multiple projections although some produced a single projection. Then, we extended our search by creating all combinations of alternative inputs and deletions that were summarized in an AIs-Deletions matrix. We found a genetic manipulation (AI-Ste5p ste2Δ) that enhanced the formation of multiple projections. Following up this lead, we demonstrated that AI-Ste4p and AI-Ste5p were sufficient to produce multiple projections when combined. Further, we showed that overexpression of a membrane-targeted form of Ste5p alone could also induce multiple projections. Thus, we successfully re-engineered the multiple projections mating morphology using alternative inputs without α-factor

Risk and prognostic significance of tuberculosis in patients from The TREAT Asia HIV Observational Database

10.1186/1471-2334-9-46BMC Infectious Diseases94

A framework for discrete stochastic simulation on 3D moving boundary domains

We have developed a method for modeling spatial stochastic biochemical reactions in complex, three-dimensional, and time-dependent domains using the reaction-diffusion master equation formalism. In particular, we look to address the fully coupled problems that arise in systems biology where the shape and mechanical properties of a cell are determined by the state of the biochemistry and vice versa. To validate our method and characterize the error involved, we compare our results for a carefully constructed test problem to those of a microscale implementation. We demonstrate the effectiveness of our method by simulating a model of polarization and shmoo formation during the mating of yeast. The method is generally applicable to problems in systems biology where biochemistry and mechanics are coupled, and spatial stochastic effects are critical

Modelling of Yeast Mating Reveals Robustness Strategies for Cell-Cell Interactions

Mating of budding yeast cells is a model system for studying cell-cell interactions. Haploid yeast cells secrete mating pheromones that are sensed by the partner which responds by growing a mating projection toward the source. The two projections meet and fuse to form the diploid. Successful mating relies on precise coordination of dynamic extracellular signals, signaling pathways, and cell shape changes in a noisy background. It remains elusive how cells mate accurately and efficiently in a natural multi-cell environment. Here we present the first stochastic model of multiple mating cells whose morphologies are driven by pheromone gradients and intracellular signals. Our novel computational framework encompassed a moving boundary method for modeling both a-cells and α-cells and their cell shape changes, the extracellular diffusion of mating pheromones dynamically coupled with cell polarization, and both external and internal noise. Quantification of mating efficiency was developed and tested for different model parameters. Computer simulations revealed important robustness strategies for mating in the presence of noise. These strategies included the polarized secretion of pheromone, the presence of the α-factor protease Bar1, and the regulation of sensing sensitivity; all were consistent with data in the literature. In addition, we investigated mating discrimination, the ability of an a-cell to distinguish between α-cells either making or not making α-factor, and mating competition, in which multiple a-cells compete to mate with one α-cell. Our simulations were consistent with previous experimental results. Moreover, we performed a combination of simulations and experiments to estimate the diffusion rate of the pheromone a-factor. In summary, we constructed a framework for simulating yeast mating with multiple cells in a noisy environment, and used this framework to reproduce mating behaviors and to identify strategies for robust cell-cell interactions