Pinetree: a step-wise gene expression simulator with codon-specific translation rates

Motivation Stochastic gene expression simulations often assume steady-state transcript levels, or they model transcription in more detail than translation. Moreover, they lack accessible programming interfaces, which limits their utility. Results We present Pinetree, a step-wise gene expression simulator with codon-specific translation rates. Pinetree models both transcription and translation in a stochastic framework with individual polymerase and ribosome-level detail. Written in C++ with a Python front-end, any user familiar with Python can specify a genome and simulate gene expression. Pinetree was designed to be efficient and scale to simulate large plasmids or viral genomes. Availability Pinetree is available on GitHub (https://github.com/benjaminjack/pinetree) and the Python Package Index (https://pypi.org/project/pinetree/).

[1]  Michael Margaliot,et al.  On the Ribosomal Density that Maximizes Protein Translation Rate , 2016, PloS one.

[2]  John Yin,et al.  Kinetic Modeling of Virus Growth in Cells , 2018, Microbiology and Molecular Biology Reviews.

[3]  D. Endy,et al.  Intracellular kinetics of a growing virus: a genetically structured simulation for bacteriophage T7. , 1997, Biotechnology and bioengineering.

[4]  Isaac Meilijson,et al.  Genome-Scale Analysis of Translation Elongation with a Ribosome Flow Model , 2011, PLoS Comput. Biol..

[5]  S. J. Plimpton,et al.  A partial-propensity variant of the composition-rejection stochastic simulation algorithm for chemical reaction networks , 2010 .

[6]  Joseph H. Solomon,et al.  The Morphology of the Rat Vibrissal Array: A Model for Quantifying Spatiotemporal Patterns of Whisker-Object Contact , 2011, PLoS Comput. Biol..

[7]  Michael Margaliot,et al.  A model for competition for ribosomes in the cell , 2015, Journal of The Royal Society Interface.

[8]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[9]  J. Plotkin,et al.  Rate-Limiting Steps in Yeast Protein Translation , 2013, Cell.

[10]  Drew Endy,et al.  TABASCO: A single molecule, base-pair resolved gene expression simulator , 2007, BMC Bioinformatics.

[11]  D. Endy,et al.  Computation, prediction, and experimental tests of fitness for bacteriophage T7 mutants with permuted genomes. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[12]  A. Subramaniam,et al.  Kinetic modeling predicts a stimulatory role for ribosome collisions at elongation stall sites in bacteria , 2017, eLife.

[13]  D. Gillespie.  Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .