MOLNs: A cloud platform for interactive, reproducible and scalable spatial stochastic computational experiments in systems biology using PyURDME

Computational experiments using spatial stochastic simulations have led to important new biological insights, but they require specialized tools, a complex software stack, as well as large and scalable compute and data analysis resources due to the large computational cost associated with Monte Carlo computational workflows. The complexity of setting up and managing a large-scale distributed computation environment to support productive and reproducible modeling can be prohibitive for practitioners in systems biology. This results in a barrier to the adoption of spatial stochastic simulation tools, effectively limiting the type of biological questions addressed by quantitative modeling. In this paper, we present PyURDME, a new, user-friendly spatial modeling and simulation package, and MOLNs, a cloud computing appliance for distributed simulation of stochastic reaction-diffusion models. MOLNs is based on IPython and provides an interactive programming platform for development of sharable and reproducible distributed parallel computational experiments

Systems Level Modeling of the Cell Cycle Using Budding Yeast

Proteins involved in the regulation of the cell cycle are highly conserved across all eukaryotes, and so a relatively simple eukaryote such as yeast can provide insight into a variety of cell cycle perturbations including those that occur in human cancer. To date, the budding yeast Saccharomyces cerevisiae has provided the largest amount of experimental and modeling data on the progression of the cell cycle, making it a logical choice for in-depth studies of this process. Moreover, the advent of methods for collection of high-throughput genome, transcriptome, and proteome data has provided a means to collect and precisely quantify simultaneous cell cycle gene transcript and protein levels, permitting modeling of the cell cycle on the systems level. With the appropriate mathematical framework and sufficient and accurate data on cell cycle components, it should be possible to create a model of the cell cycle that not only effectively describes its operation, but can also predict responses to perturbations such as variation in protein levels and responses to external stimuli including targeted inhibition by drugs. In this review, we summarize existing data on the yeast cell cycle, proteomics technologies for quantifying cell cycle proteins, and the mathematical frameworks that can integrate this data into representative and effective models. Systems level modeling of the cell cycle will require the integration of high-quality data with the appropriate mathematical framework, which can currently be attained through the combination of dynamic modeling based on proteomics data and using yeast as a model organism

Ribosome traffic on mRNAs maps to gene ontology : genome-wide quantification of translation initiation rates and polysome size regulation

Peer reviewedPublisher PD

A Hybrid Model of Mammalian Cell Cycle Regulation

The timing of DNA synthesis, mitosis and cell division is regulated by a complex network of biochemical reactions that control the activities of a family of cyclin-dependent kinases. The temporal dynamics of this reaction network is typically modeled by nonlinear differential equations describing the rates of the component reactions. This approach provides exquisite details about molecular regulatory processes but is hampered by the need to estimate realistic values for the many kinetic constants that determine the reaction rates. It is difficult to estimate these kinetic constants from available experimental data. To avoid this problem, modelers often resort to ‘qualitative’ modeling strategies, such as Boolean switching networks, but these models describe only the coarsest features of cell cycle regulation. In this paper we describe a hybrid approach that combines the best features of continuous differential equations and discrete Boolean networks. Cyclin abundances are tracked by piecewise linear differential equations for cyclin synthesis and degradation. Cyclin synthesis is regulated by transcription factors whose activities are represented by discrete variables (0 or 1) and likewise for the activities of the ubiquitin-ligating enzyme complexes that govern cyclin degradation. The discrete variables change according to a predetermined sequence, with the times between transitions determined in part by cyclin accumulation and degradation and as well by exponentially distributed random variables. The model is evaluated in terms of flow cytometry measurements of cyclin proteins in asynchronous populations of human cell lines. The few kinetic constants in the model are easily estimated from the experimental data. Using this hybrid approach, modelers can quickly create quantitatively accurate, computational models of protein regulatory networks in cells

Establishment of a robust single axis of cell polarity by coupling multiple positive feedback loops

Establishment of cell polarity-or symmetry breaking-relies on local accumulation of polarity regulators. Although simple positive feedback is sufficient to drive symmetry breaking, it is highly sensitive to stochastic fluctuations typical for living cells. Here, by integrating mathematical modelling with quantitative experimental validations, we show that in the yeast Saccharomyces cerevisiae a combination of actin- and guanine nucleotide dissociation inhibitor-dependent recycling of the central polarity regulator Cdc42 is needed to establish robust cell polarity at a single site during yeast budding. The guanine nucleotide dissociation inhibitor pathway consistently generates a single-polarization site, but requires Cdc42 to cycle rapidly between its active and inactive form, and is therefore sensitive to perturbations of the GTPase cycle. Conversely, actin-mediated recycling of Cdc42 induces robust symmetry breaking but cannot restrict polarization to a single site. Our results demonstrate how cells optimize symmetry breaking through coupling between multiple feedback loops

Deterministic Intracellular Modeling

The United States Air Force is interested in the potential side effects at the cellular level from exposure to mission-essential chemicals. Presently, Air Force toxicology studies are conducted to help shed light in identifying potential hazards to workers. However, it takes a considerable amount of money, resources, and time to obtain and analyze experimental results from toxicology studies. The necessity for innovative methods that enable researchers to more effectively generate and analyze data is apparent