text:
This is a forerunner as to what lies ahead once the NBA Finals and Stanley Cup Finals are completed with baseball only unless you have a soccer preference. Today the Cardinals have fallen close to 20 cents with Matt Holiday in question and Milwaukee announcing their starting pitcher. Even Clayton Kershaw and the Dodgers are off almost 30 cents, with Colorado on a four-game win streak. Like we say every day, understanding line moves can build your bankroll. (550-498 L203D)
MLB – (901) CHIC. CUBS at (902) MIAMI 7:10 ET CSN-C, FSFL
The expected game time temperature today in Miami is 87, which is 40 degrees warmer than it was at Wrigley Field yesterday for the Cubs. While Chicago will be in much more comfortable conditions, those making MLB picks are less sure the Cubs can beat the Marlins and dropped them 15 cents to -125. Chicago had dropped 7 of 11 prior to winning yesterday against Kansas City and has Jason Hammel (3-2, 2.98 ERA) is working on seven days rest and he's completed seven or more innings in five of his last six starts. Miami has played a little better at .500 in last eight games, but after losing to the Mets yesterday is 3-10 at home after a loss this season. Doug's Choice – Chicago wins
MLB – (907) ATLANTA at (908) ARIZONA 9:40 ET SPSO, FSAZ
Atlanta hits the desert having scored 15 runs against Lincecum, Bumgarner and the usually reliable San Francisco bullpen and will take on Arizona's Archie Bradley (2-2, 5.28), who has a 11.81 ERA since being hit in the puss. With these two factors and Paul Goldschmidt on fire, the total has been lifted from 8 to 9. Braves starter Alex Wood (3-2, 3.50) has limited the D-Backs to a .172 BA. Those constructing sports picks are not thinking Wood will have that kind of success, as home teams when the total is 8.5 to 10, scoring 4.3 to 4.8 runs a game, against a good NL starter (ERA of 3.70 or lower), with an average bullpen whose ERA is 3.75 to 4.50, are 36-14 OVER since 2011. Doug's Choice – Lean Over
MLB – (917) TAMPA BAY at (918) L.A.A. ANGELS 10:05 ET CSN-C, ROOT-SW
Garrett Richards (4-3, 3.12) and June have gotten along well with him posting a 6-2 record and 2.64 ERA in 10 starts and eight relief appearances. Yet there is still hesitancy and the Angels along with their right-hander have sunk from -165 to close to -150 at many sportsbooks. That is a little shocking since the Halos just swept four games from Detroit at the Big A and Tampa Bay's Alex Colome (3-1, 4.55) on the surface is nothing special. However, Colome has been a tough hombre on the road at 4-0 with a 1.78 ERA in six career road starts. Plus, coming off their offensive outburst in Baltimore, the Rays are 17-5 after a game where they had six or more extra base hits the last three seasons. Even with this still prefer L.A. at home. Doug's Choice – L.A.A. wins
MLB – (919) N.Y. YANKEES at (920) SEATTLE 10:10 ET YES, ROOT-NW
Michael Pineda (6-2, 3.36) faces the team which traded him away for the first time but his mound opponent will be none other than the King, Felix Hernandez (8-1, 1.91). Though the Yankees lost a four-game series to lowly Oakland, the Mariners were no better in doing the same at home to Cleveland. The betting odds on the New York are down 15 cents to +135 despite King Felix allowing just six runs in his past five starts versus the Yanks (Mariners 3-2). If Pineda can keep his emotions in check, he could hold down Seattle's chances of scoring, nonetheless, Hernandez is 15-4 when the money line is -100 to -150 the last two seasons.(Mariners record) Yet here's an oddity, Hernandez and the M's are only 42-53 when he pitches and they are off two or more losses. Doug's Choice – Lean Seattle
Vegas Pro Insiders Daily Trends
MLB Money Line Trend
Clay Buchholz and Boston are 3-15 in home games the last two seasons.
MLB Totals Trend
Clayton Kershaw of the Dodgers is 15-4 UNDER after giving up no earned runs last outing the last three seasons.
MLB Run Line Trend
Ubaldo Jimenez of Baltimore is 18-38 if his team is off two or more consecutive losses in his career.


text:
Summer Mount Yudono Sun return Courses (Yamagata family edition)
How about going to visit mystery and the mysterious ground?
Do you not enjoy by fishing pond experience noisily afterwards?
Fish which we caught by oneself wants to appreciate beauty of particular Wed.
And in place where Japanese one of Yamagata Prefecture does not come to!
Do you not make cool-down with fountain and chilled ice-cream cone whirled up beautifully?
There was tall building at the station square. View is the best!
You may spend time relaxedly at parting time of Yamagata Prefecture.
Yamagata Station
Yamagata Station is station where there is much access Shinkansen and conventional line that each place goes to visit the starting point well at Yamagata Station to be similar. There are many rent-a-car shops near station, too.
Car 60 minutes
We advance to entrance of Hongu for stretch using Mount Yudono toll road, worship bus from parking lot. When weave bus, and advance to visit, scenery changes completely, too; to another world. Experience to be possible only here. We can enjoy "senjin*onsen" for 500 yen (from 9:00 to 15:00) in parking lot side "yudensansanryokagosho" and can enjoy two kinds of lunch set (from 11:30 to 14:00).
It is 15 minutes by car
Tamugi fishing pond
We can enjoy char a kind of trout casually in fishing pond. He/she gives even if we deal with char-grilling for take-outs on the spot to be similar.
※Pole, bait, processing, fish (g) is charged separately.
It is 15 minutes by car
Mt. Gassan size fountain
The best large fountain in Japan! We have not looked at such a fountain. Please see fountain show that is with commentary, and is luxurious by all means. We perform display utilized property and charm to have of Wed with culture building of Wed. We can observe in no charge for admission. We can taste chilled chipped ice and ice-cream cone at stand.
It is 45 minutes by car
Yamagata Station
There are a lot of restaurants that sale here dishes of "special product of Yamagata Prefecture" can enjoy in front of Yamagata Station. We can look around state of town and the magnificent mountains from observation deck of Ka-jo Castle central building.

subset: commoncrawl


text:
Open Access This article is
- freely available
- re-usable
Computation 2018, 6(1), 9; doi:10.3390/computation6010009
An Overview of Network-Based and -Free Approaches for Stochastic Simulation of Biochemical Systems
Center for Quantitative Medicine and Department of Cell Biology, University of Connecticut School of Medicine, 263 Farmington Av., Farmington, CT 06030-6033, USA
Author to whom correspondence should be addressed.
Received: 9 December 2017 / Accepted: 29 January 2018 / Published: 31 January 2018
Stochastic simulation has been widely used to model the dynamics of biochemical reaction networks. Several algorithms have been proposed that are exact solutions of the chemical master equation, following the work of Gillespie. These stochastic simulation approaches can be broadly classified into two categories: network-based and -free simulation. The network-based approach requires that the full network of reactions be established at the start, while the network-free approach is based on reaction rules that encode classes of reactions, and by applying rule transformations, it generates reaction events as they are needed without ever having to derive the entire network. In this study, we compare the efficiency and limitations of several available implementations of these two approaches. The results allow for an informed selection of the implementation and methodology for specific biochemical modeling applications.
Keywords:stochastic simulation; modeling; network-based; network-free; rule-based modeling; systems biology
1. Introduction
Research in systems biology has been increasingly supported by computational models of biochemical reaction networks. These models are studied either through a deterministic approach, using differential equations to represent the temporal changes of the concentrations of the chemical species, or via a stochastic approach based on the chemical master equation (CME) and solved through Monte Carlo simulation algorithms. Although the deterministic approach of solving a set of differential equations by numerical integration is fast and widely adopted, it is unable to estimate the variance of the species concentrations and can become inaccurate for systems with a small number of particles [1]. In these situations, the stochastic approach of using a Monte Carlo simulation algorithm for evaluation of the CME is preferred.
The CME is a very high dimension differential equation that describes the evolution of the entire state space. Direct solutions of the CME are rare and apply only to very small systems. In practice, the CME is solved by applying a simulation algorithm that provides an exact solution in the Monte Carlo sense (i.e., by summation of many simulated trajectories). Gillespie derived an algorithm that satisfies this requirement [2]. This is often referred to as the Gillespie algorithm, although Gillespie himself referred to it as the stochastic simulation algorithm (SSA). In fact, Gillespie provided two variants of the algorithm, the direct method and the first reaction method, with the direct method being the most widely used. Gillespie's SSA has been further improved for computational efficiency [3,4,5] and there have been several implementations of these algorithms in software for computational systems biology. For example, there are software applications, such as COPASI [6], VCell [7], and StochPy [8], which provide user-friendly platforms to create and simulate models using the SSA and include other features to further analyze the model and simulation results. Furthermore, there are other lightweight programs developed specifically to simulate models using the SSA, namely Dizzy [9], Gillespie2 [10], SGNS2 [11], RoadRunner [12] and pSSAlib [13]. All these simulators require the full set of reactions—the reaction network—to be enumerated beforehand and thus are sometimes termed "network-based".
There are several cases in which a biochemical network is very large or limitless. A particularly common occurrence is given by some signal transduction networks that contain proteins with multi-site phosphorylation, leading to combinatorial numbers of chemical species and reactions between them [14,15,16]. Another case is that of the formation of polymers with an unlimited number of monomers. In order to model such systems, an approach has been developed in which sets of similar reactions are defined by rules that apply to sets of species specified by patterns [17,18,19]. This formalism results in a concise model specification of the underlying chemical kinetics [20,21]. The most common rule-based modeling languages are the BioNetGen language (BNGL) [17] and Kappa [22]. The BNGL simulator, BioNetGen [17,23], operates by deriving the reaction network specified in the reaction rules and then applying the SSA for simulation. On the other hand, rule-based simulators such as KaSim, PySB [24], RuleMonkey [25] and NFsim [19] carry out simulations directly on the basis of reaction rules without deriving the entire reaction network, and accordingly these have been termed "network-free". At its core, all these simulators are based on Gillespie's method, as rules are sampled at each time interval using a method equivalent to how reactions are sampled in the SSA.
In this review, we compare these two stochastic simulation approaches and several popular software implementations in the context of models with different complexity. The comparison addresses issues such as the number of particles, species, and reactions, as well as the length of the simulation.
2. Network-Based Approach
The stochastic formulation of chemical kinetics describes the time evolution of a well-stirred set of chemically interacting particles in thermal equilibrium within a fixed reaction volume [1]. The time evolution of the number of particles of each species in the volume, on the basis of the probabilities of all reactions that can occur in the system, is driven by the CME. As already mentioned, the CME is rarely solved analytically, mostly because the number of its terms grows exponentially with the number of species in the system.
An alternative to the analytical solution of the CME is to simulate the trajectories of molecular populations in exact accordance with the CME, as proposed by Gillespie [2] (the SSA). Each trajectory corresponding to a single SSA run represents an exact sample from the distribution defined by the CME. The steps of SSA can be summarized as follows:
- Initialize: Set the time and set up the initial state vector, propensities, and random number generators.
- Execute: Using a suitable sampling procedure, generate random numbers and, on the basis of these, determine the next reaction to occur and the time interval.
- Update: Update the molecule count, and if needed, recalculate the propensities. Output the system state.
- Iterate: If simulation end time is not reached, go to step 2.
The two original, and statistically equivalent, sampling procedures for step 2 of the SSA are the direct method (DM) and the first reaction method (FRM) [2]. The DM samples two random numbers from the uniform distribution in the unit interval, and the time of next reaction () is first generated according to the probability function of reactions. Using , the DM then generates the indices of reactions and selects the one to occur next. The FRM, using a random number, generates "tentative reaction times" () for all the reactions and then selects the reaction with the smallest . Because the FRM needs to generate many more random numbers per iteration than the DM (for systems with three or more species), the DM is generally the procedure implemented for the sampling in step 2 of the SSA [2].
Gibson and Bruck proposed the next reaction method (NRM) [4] that can reduce the computational costs of the SSA significantly. In addition to using one random number per iteration, to reduce the time to update propensities and to find the smallest value, the NRM uses an indexed priority queue to store the values generated in previous iterations and to extract them whenever required. This results in a significant improvement in the runtime performance when compared to the FRM. This algorithm is exact as well as efficient. For large reaction networks and loosely coupled reaction systems, the NRM is significantly faster than both the FRM and the DM. This advantage, however, may not be significant for small systems, as the computational cost of maintaining the additional data structures required dominates the simulation time [5].
Other variants to accelerate the search for the next reaction in the SSA have been proposed, such as the optimized direct method (ODM) [5], the sorting direct method (SDM) [26], the partial-propensity direct method (PDM) [27], and the SSA with composition rejection algorithm (SSA–CR) [28,29,30].
Besides the exact algorithms mentioned above, many others have been proposed that can accelerate the simulation even further, but they do this by adopting approximations and no longer provide exact solutions. A popular method is the -leaping algorithm [31], which does not simulate each reaction event individually but rather steps a time-span and estimates how many and which reactions have happened meanwhile. Many other variants of this and other approximations have been proposed, including hybrid methods that partition the system into a part that is simulated using differential equations and another that uses the SSA or one of its variants (see review by Pahle [32]).
Most of the stochastic simulators provide options to choose between the DM and the NRM, for example, COPASI [6], StochPy [8], and Dizzy [9]. Other simulators use only one of these, with Gillespie2 [10] and RoadRunner [12] using only the DM, and SGNS2 [11] using only the NRM. The pSSAlib software [13] allows selection between the DM, the PDM, a sorting variant of the PDM (SPDM), and the SSA–CR. StochKit2 [33] provides several of these, including the SSA–CR, but automatically selects which algorithm to use.
3. Network-Free Approach
To address the combinatorial complexity in biological signaling networks [14], originating from multiple post-translational modifications and conformational changes, rule-based modeling approaches have been developed [15,17,18,20,21,22,34]. At the core of these approaches are reaction rules that represent groups of reactions. These rules refer to specific binding sites with or without specific ligands. Rules can also specify different states of a molecule (such as oxidized or reduced, phosphorylated or unphosphorylated, etc.). With rule-based modeling it is easy to specify, with a few rules, a complex set of combinatorial interactions in which several subunits can assemble into larger complexes and allow for modification of specific moieties. This type of model specification is therefore very useful for signal transduction networks in which these types of interactions are abundant.
The BNGL [17,35,36,37], the -language [18], and -calculus [38] are some examples of formalisms developed for biochemical rule-based modeling. While the BNGL can be processed by different software applications (BioNetGen [23,36], DYNSTOC [25], RuleMonkey [39], and NFsim [19]), the other languages are mostly restricted to being processed by a single software package. The BioNetGen software package expands a BNGL rule-based model to a reaction network, which is then simulated using a variety of deterministic and stochastic network-based methods. However, when a rule-based model can result in a large reaction network, the expansion as well as the simulation of such a network becomes computationally expensive. For such scenarios, the generation of the reaction network can be avoided by a network-free simulation approach.
DYNSTOC [25] uses an agent-based null-event stochastic simulation approach based on an earlier package, STOCHSIM [40]. In this approach, each of the reactive molecular components are represented as a software object (agent), and these are tracked individually during the simulation. More specifically, for each fixed time increment, on the basis of a decision to select either one or two molecules for the next reaction, the reactants are first chosen randomly. Then the rules that qualify for the interaction on the basis of the chosen reactants are shortlisted, and for the reaction with the highest probability, an update is performed using a graph-rewriting operation.
Unlike DYNSTOC, which uses a fixed time-step, RuleMonkey [39] has a variable time increment, and rules are represented as pattern graphs. The simulation procedure is similar to the SSA [2], as the time increment and rule selection are based on the DM. Once a rule to execute next is chosen, the most potential reactants are selected on the basis of the pattern graphs and are then used to update the state of system. As such, RuleMonkey uses iterative updates to track rule rates exactly, avoiding null events that do not change the state of the system being simulated [41]. NFsim [19] is another rule-based simulator, using a generalized algorithm [42] also based on the SSA. Contrary to RuleMonkey, NFsim introduces null events in its implementation [41]. While both RuleMonkey and NFsim have been shown to perform similarly over a wide parameter range [41], NFsim has the additional capabilities of defining functional rates and coarse-grained rules. It uses an efficient representation of molecules, complexes, and rules as well as an optimized handling of reactant selection and transformation.
These network-free simulators, unlike the network-based SSA simulators, scale with the number of rules rather than the number of reactions and thus should be very efficient for systems in which a few rules can represent a large number of reactions [37]. While this is true for networks with limited interacting particles, the network-free simulators might not be as efficient, given they represent each particle individually. The particle-specific events, such as aggregation and polymerization, make the computational cost even higher. On the other hand, although the network-based simulators are dependent on the number of reactions, their efficiency is not affected greatly by the number of molecules [37]. As such, network-based simulators may be preferred for systems with large particle numbers and a moderate reaction network.
All of the approaches described have difficulties when there are large numbers of particles and a large reaction network. To address this situation, there have been efforts to develop hybrid methods [37,43,44,45]. The hybrid particle-population-based approach [37] is reported to be exact and efficient but requires a predefined partition of the system into network-free and -based parts. Because these hybrid approaches are based on the partitioning of the models, it is important to identify limits of both network-based and -free approaches such that automatic identification of different parts of the network could be created on the basis of this information [37].
4. Benchmarking Stochastic Simulation
4.1. Simulators
A survey of the literature reveals that several of the stochastic simulators described above are regularly used in computational systems biology but to our knowledge have never been compared for performance in a systematic way. Given that in stochastic simulation one must define the time-dependent distribution (or at least some statistics of this distribution, such as the mean and standard deviation), this usually requires repeating simulations many times; thus the performance of the simulators used may be a critical factor. We profiled a number of the most widely used simulators with a set of models of increasing complexity. The intention was not only to compare simulators on the basis of similar algorithms, but also to compare the different algorithms used. We note that we only considered software implementing methods that are exact solutions of the CME; approximate and hybrid methods were excluded.
We identified a series of commonly used and freely available network-based and -free simulators, which are listed in Table 1. Because we needed to profile these on the same computer to be able to compare them, only software packages that could be run on a local machine were included. Moreover, to be able to specify the same model across all of these, we used models specified with the BNGL; thus we restricted the selection to packages that could either process BNGL directly or could import models in the systems biology markup language (SBML) [46] (BioNetGen was used to generate the reaction network and export it in SBML format). Three exceptions were included: KaSim, which uses the -language rather than BNGL; pSSAlib, which has an implementation of the PDM and the SSA–CR; and StochKit2, which has an implementation of the SSA–CR. For KaSim, we translated the models into the -language. For pSSAlib, we created a program that converts the standard SBML into the specific dialect it can understand (while pSSAlib claims to read SBML, it requires specific annotations in the files). For StochKit, we used the SBML converter that was provided with that package.
4.2. Models
To compare the performance of simulators under different conditions, we selected models with increasing complexity, as quantified by the number of species and reactions. Table 2 summarizes the models considered in this study. The first two models, "multi-state" and "multi-site", are conceptual and have been used for the illustration of basic biochemical networks. We expected the derivation as well as simulations of these models to be fast, as they have a small number of species and reactions. The remaining three models were originally formulated to study specific signaling networks. They are more complex than the previous two and allowed us to test the simulators under more realistic conditions.
The multi-state model is composed of three species, R, L, and A; R and L can form a complex "R.L", and the latter can dissociate back to the monomers. The species A can bind to R, and it can exist in a phosphorylated or unphosphorylated state (see Figure A1 for details). This model was previously described in [17,25].
The multi-site model contains the same three species, but here A has three different phosphorylation sites; L binds A to any of its phosphorylated sites, and R binds A to its unphosphorylated sites (see Appendix A, Figure A2).
The Epidermal growth factor receptor (EGFR) signaling model describes the early signaling events in the epidermal growth factor receptor cascade [47]. Besides the epidermal growth factor and its receptor, the model consists of the adapter proteins Grb2 and Shc, and EGF-induced activation of the guanine nucleotide exchange factor Sos. The reaction network described by this model contains 356 molecular species and 3749 reactions (see Appendix A, Figure A3 and Supplementary model file for details).
The B-cell receptor (BCR) signaling model has been used to investigate the early events in B-cell antigen receptor signaling [48], particularly the roles of the Src family protein tyrosine kinases Lyn and Fyn, which regulate the activities and fates of B cells. The model includes 1152 species and 24,388 reactions (see Appendix A, Figure A4 and Supplementary model file for details).
Finally, the high-affinity human IgE receptor (FcϵRI) signaling model represents the early events in Fc epsilon receptor (FcϵRI) signaling [49]. This model consists of the interactions between FcϵRI, Lyn, Syk, and a bivalent ligand that aggregates FcϵRI. Several variants of this model have been used previously for testing the performance of some network-free simulators [19]. In this study, we used the variant with 2 sites in the receptor, which consists of 24 rules, generating a network with 3744 species and 58,276 reactions (see Appendix A, Figure A5).
To include even larger networks, we attempted to use a model of ErbB-mediated activation of the protein kinases ERK and AKT [50] and a model of early T-cell receptor signaling [51]. Both these models are composed of hundreds of rules each, and we were unable to generate the network with BioNetGen because of the excessive memory requirement by this application (in a computer with 32 GB of RAM). Therefore, these were not used for the profiling, and they are examples of systems that currently can only be simulated with network-free methods.
5. Results
We performed two sets of tests to probe the performance and scaling of each simulator. The first set was intended to test how the simulators behave in the presence of increasing numbers of particles in the system. This was achieved by setting increasing values for the initial conditions of each species (summarized in Appendix B, Table A1). In the second set, we tested how the simulators scale in increasingly longer simulations, which was achieved by requesting longer end times (Appendix B, Table A1).
5.1. Increasing Numbers of Particles
We ran simulations of all the models in Table 2 for a fixed end time of 100 s (simulation time), using the initial conditions in Appendix B, Table A1. Figure 1 depicts the behavior of each simulator for increasing molecule numbers and for each model.
For the multi-state model (Figure 1A), we found that BioNetGen, pSSAlib, and SGNS2 were the fastest. Interestingly, for initial conditions with very few molecules, the network-free implementations KaSim, NFsim, and RuleMonkey were faster than the remaining network-based SSA implementations; however, they become slower than other tools for a larger number of molecules. The execution times of the network-based COPASI (both in DM and NRM), Dizzy, and RoadRunner applications were mostly invariant with the number of molecules, indicating that these tools have a large overhead at the time of loading the model but otherwise were fast. DYNSTOC was fast for very few molecules but very quickly became the slowest, showing the extreme dependency of this approach on the number of particles in the system. Finally StochPy was the second slowest. This was partly perhaps it is written in an interpreted language (Python), but also because this tool outputs every single reaction event, unlike the other tools that allow arbitrary sampling intervals (here we requested them to produce 1000 intervals along the time course; see Methods).
For the multi-site model, the scaling of the execution times with the number of particles was qualitatively similar for all simulators. Although there were clear differences between their execution times, as shown in Figure 1B, there was no clear separation between the network-free and -based implementations. While NFsim was faster than some of network-based simulators, it was slower than BioNetGen, addedpSSAlib, RoadRunner, SGNS2 and StochKit2. These simulators, along with COPASI and RuleMonkey, completed all the simulations within 100 s of a wall-clock time. On the other hand, for large number of particles, Dizzy, Gillespie2 and KaSim took more than 100 s to complete.
In the EGFR signaling model, only 8 of the 13 simulators could complete all simulations within a threshold of 2000 s. For this model, in the lower extreme of the molecule numbers, NFsim and RuleMonkey were the fastest. However with an increasing number of molecules, BioNetGen and pSSAlib_SPDM became the fastest (Figure 1C). Under these conditions, the execution times of NFsim and RuleMonkey became similar to those of COPASI and RoadRunner.
The BCR signaling model could not be simulated using RuleMonkey because of a "Non-binding bimolecular reaction" error (meaning that it could not deal with the complete BNGL). The conversion of the SBML model to the StochKit2 format could not completed, and therefore it was not benchmarked. Among the remaining six simulators, pSSAlib_SPDM was the fastest for a larger number of molecules. For a low number of molecules, pSSAlib_SPDM was slower than BioNetGen, NFsim and SGNS2. All of these closely followed pSSAlib_SPDM, with a similar scaling of their execution times for a larger number of molecules. On the other hand, COPASI and RoadRunner had a visibly large overhead at the start but scaled less dramatically with the increase in the number of particles. For a larger number of particles, COPASI, NFsim, pSSAlib_SSACR and RoadRunner could not complete the simulations within a threshold of 5000 s (wall-clock time).
Both COPASI and RoadRunner were not able to load the FcϵRI signaling model within 5000 s. Although pSSAlib_SPDM had a significant overhead, it had a constant time scaling of the execution time and was the fastest for a large initial number of molecules (See Appendix A, Figure A6); pSSA_SSACR, on the other hand, scaled almost linearly and could not complete all simulations within 5000 s. For the other three simulators, we observed execution time patterns similar to that of the BCR signaling model. BioNetGen was followed by SGNS2, which was marginally slower than NFsim. Given that for this model, the derivation of the network from BNGL takes a considerable amount of time, the choice of using NFsim appears to be advantageous (see Appendix A, Figure A6).
5.2. Dependency on the Simulation End Time
In the second test, we measured the execution times as a function of the simulation end time, keeping the initial number of molecules fixed. Figure 2 depicts the scaling of the execution times with different simulation end times for all the models and simulators.
The results for the multi-state model (Figure 2A) were similar to those of the previous test with increasing initial molecule numbers. Again, DYNSTOC and StochPy were the slowest. For this model, we observed a clear separation of the network-free implementations and the other network-based simulators. More specifically, we found that network-based simulators, except StochPy, were either very fast or invariant with the increasing simulation end time. On the other hand, all the network-free simulators became slower with increasing simulation end time. Also consistent with the previous results on particle numbers, the network-based simulators were the fastest for longer simulation end times.
For the multi-site model (Figure 2B), we found that the execution times of all the simulators scaled similarly to the multi-state model. For this model, only 9 of the 13 simulators completed all simulations within a threshold of 2000 s. DYNSTOC, StochPy, KaSim and Dizzy could not complete all the simulations; pSSAlib_SPDM was the fastest throughout all the simulation end times, followed by BioNetGen, RoadRunner, NFsim, SGNS2 and StochKit2.
For the EGFR signaling model (Figure 2C), COPASI and RoadRunner could not complete all the simulations within 5000 s. Only NFsim and RuleMonkey (network-free) and BioNetGen, pSSAlib, SGNS2 and StochKit2 (network-based) were able to finish all simulations under the time threshold. The network-based simulators were significantly faster than the network-free simulators, whose execution times were in a range similar to those of COPASI and RoadRunner (for the simulations that they could finish).
COPASI and RoadRunner had significantly long execution times for the BCR signaling model (Figure 2D), despite that this increased minimally with the simulation end time. Once again, this is a reflection of these tools' overhead in model loading and associated "house-keeping tasks". It should be noted that, as in the first test, RuleMonkey and StochKit could not be used for comparison. Only BioNetGen, NFsim, pSSAlib, and SGNS2 were able to run all the simulations and scale with a similar pattern; pSSAlib_SPDM was the fastest. These simulators were also the only simulators that could run the FcϵRI signaling model in a reasonable time (Figure 2E). Once again, pSSAlib_SPDM had a constant time scaling of the execution time and was the fastest for larger simulation end times; pSSAlib_SSACR could not complete all the simulations within the time threshold of 5000 s. BioNetGen was the second fastest, and only a marginal difference was seen between the execution times of NFsim and SGNS2. In the execution times of the FcϵRI signaling model, it appeared that NFsim was starting to become competitive with the network-based simulators. Given the time taken to generate the network from rules, there could be an advantage to using this network-free simulator for models of this dimension.
6. Discussion
Benchmarking several simulators on the basis of different network-based or -free algorithms has revealed interesting patterns. The benchmarks have tested the scaling of the execution times of simulators on the basis of two tests, namely, as a function of the number of particles in the system and as a function of simulation end times. The first test was focused on exposing issues depending on the number of particles in the system; the second test was focused on exposing issues that arised closer to stable states (attractors). We note that most realistic simulations should start from steady states rather than the idealized state of only a few input molecules; biological systems are in stable states at the start of most experiments, which usually apply a perturbation forcing a transition between stable states (either different steady states or stable oscillations). Thus the behaviors over the longer times in the second test were rather important.
The results of both the tests indicate that StochPy and DYNSTOC are the slowest implementations. In the case of StochPy, a SSA–DM implementation, the tool is designed to write the raw simulation output after each event, rather than to do so at a requested fixed-interval output; thus we suspect that it spends a significant amount of time writing (unecessary) output. Additionally, this tool is Python-based, an interpreted language, and this likely also incurres a considerable time penalty. For DYNSTOC, the issue is rather different and is due to the algorithm used, in which each molecule is tracked as an agent. At very low molecule numbers for the simplest model, this tool was among the fastest, and it scaled linearly with the number of molecules. Thus the problem is that this approach is not able to deal with any reasonable number of molecules, which were present in every other model tested. Thus we conclude that such a pure agent-based approach is not competitive.
The results show that, surprisingly, the network-based approach was always the fastest, not only for small- and moderate-sized reaction networks, but also for the larger networks and under all the conditions tested. For the simpler models, with a relatively low number of species and with limited interactions, the lightweight simulators BioNetGen, Gillespie2, pSSAlib_SPDM, and SGNS2 were significantly faster than all the others tested. Gillespie2 became slower when simulating larger models, but BioNetGen, pSSAlib_SPDM, and SGNS2 remained very fast under all conditions tested.
An important aspect of rule-based modeling is the concept of "observables", which are functions of the species' abundances. BioNetGen outputs the values of the observables, but SGNS2 only outputs the species abundances; thus in order to obtain values of the observables, the output of SGNS2 would require a further data processing step; pSSAlib outputs data in separate files and requires a post-processing step for this purpose. Although COPASI is slower, mostly as a result of a large overhead in model loading and preparing data structures, it can also readily output the values of the observables. RoadRunner run times are approximately similar to COPASI's, but it has the same problem of requiring post-processing as in SGNS2.
The network-based simulators tested spanned a range of different SSA approaches. Several were based on the DM (Dizzy, StochPy, Gillespie2, BioNetGen, and COPASI_D), while a couple used the NRM (COPASI_GB and SGNS2) and the CR (pSSAlib_SSACR and StochKit2); pSSAlib_SPDM, which uses PDM, was found to be the fastest for the upper extremes of our tests. In the SPDM, while the time spent on factoring out reaction propensities did not pay off for a small number of molecules or short simulations, it did lead to a significant efficiency when larger numbers of molecules were reached. Consequently, for the FcϵRI signaling model, pSSAlib_SPDM completed the simulation in half the time taken by BioNetGen in the most extreme case tested—a considerable speed-up. The efficient implementation of this method is, however, handicapped by the general usability of the simulator. One setback is that the pSSAlib does not accept any SBML file, but requires specific annotations, which other software do not include. This is particularly problematic for rule-based models that are generated from the BNGL; these models usually have very large numbers of reactions and adding the annotation manually is not practical.
As expected, the DM was observed to be less efficient than the NRM and CR. However, the comparison of the DM against the NRM in the same package (COPASI) showed that the difference is not large and is only apparent under conditions with large numbers of particles. One notable exception to this was BioNetGen, which, while being a DM implementation, was the fastest at all times. This might have been due to its implementation of the sorting variant of the DM. Further efficiency seems to have been obtained by various code optimizations (earlier versions of BioNetGen than that tested here were much slower), but the same is true for SGNS2 (a NRM implementation). We noted, however, that BioNetGen uses the standard C runtime rand() function, unlike SGNS2 and most other SSA implementations, which use the Mersenne Twister [52], a much better-quality pseudo-random number generator but that is slower (see Appendix A, Figure A7 for a comparison of the two). The dangers of using poor pseudo-random number generators are well known [53,54], and this could be a concern here, particularly for long simulation times. Also surprising is that the SSA–CR implementations (both pSSAlib_SSACR and SStochKit2) were not faster than the NRM, despite expectations of the contrary [28,30]. The expected advantage of SSA–CR did not materialize in models with large numbers of reactions. Of course, the efficiency of this implementation cannot be ruled out for other types of models not tested here.
What are the advantages of network-free simulation? While the network-free simulators were never the fastest with any of the models and conditions included here, there are clearly situations in which they are needed. The use of network-free simulation is inevitable when the derivation of the network is computationally challenging or impossible. For example, complete models of the ErbB-mediated activation of ERK and AKT [50] and of early T-cell receptor signaling [51] result in very large reaction networks, so large that BioNetGen is unable to generate the corresponding reaction networks. Even if the network could be derived, loading it in simulators such as COPASI and RoadRunner would require a significantly long time. Network-free simulators are also essential to simulate systems with infinitely linking molecules, such as models of polymerization, models of trivalent ligand bivalent receptors (TLBRs), models of large complexes, and so forth. Among the network-free simulators tested, NFsim was generally the fastest. Under some of the less-demanding conditions tested (i.e., low molecule numbers and simpler models), RuleMonkey had a small advantage, but otherwise it is clear that NFsim is currently the best choice. NFsim also provides an option to define functional or conditional rate laws and complex rules. This capability is particularly useful to model systems whose rates are affected by the availability/unavailability of specific molecules. We also identified areas in which network-free simulators require further improvements. For example, we were not able to simulate models with reactions for which several of the bound moieties suffer a transformation (catalysis) with any of the network-free simulators, as these aborted with errors when encountering the catalysis, despite that BioNetGen easily generates their network (as it should). This is a rather common occurrence that happens in every enzyme mechanism (because the reaction between substrates happens with these bound to the enzyme). An example was a model of the electron-transport chain [55], and another was a model of the cap-binding complex in mRNA translation [56]. RuleMonkey was unable to run the BCR signaling model, aborting with the error message "Non-binding bimolecular reaction". These limitations in processing valid BNGL rules affect both NFsim and RuleMonkey, but can hopefully be corrected in future versions of these packages.
The present analysis revealed that the network-based (SPDM) simulation was the most efficient method for all models tested. We hoped to have identified regimes in which network-free simulation would be more competitive. This suggests that, while a rule-based specification of the models is much simpler than enumerating all the reactions, simulation via network-free implementations is not always efficient, unless the derivation of the network from the rules is computationally intractable or there are infinitely linking species/molecules in a model. It is possible that larger models than those tested here (but with a finite number of reactions) may present conditions under which network-free simulation outperforms network-based, but this is yet to be established. On the basis of the present results, we have to conclude that for hybrid algorithms that integrate both of these approaches (e.g., [37]), the only portion of the networks that should be partitioned to be simulated by the network-free approach are those rules leading to the formation of infinite linking chains, while the rest of the network should be simulated using the network-based approach.
7. Methods
For each model, we prepared an input file appropriate for each simulator (see below); then we verified that the simulation results quantitatively matched across all simulators (a sample time trace of observables in each model is shown in Appendix A, Figure A1, Figure A2, Figure A3, Figure A4 and Figure A5). We found that the numbers of particles output for all the simulators matched. A more thorough test of the exactness of these implementations is beyond this scope, but could be carried out using the SBML stochastic test suite [57].
7.1. Model Construction
BNGL specifications for all the models used in this study were retrieved from the respective sources (Table 2). These BNGL files were used as direct input for DYNSTOC, RuleMonkey, and BioNetGen. For NFsim, we generated appropriate XML files from the BNGL with BioNetGen.
Most other simulators used SBML as input; thus we used BioNetGen to generate the network of reactions and write them out in SBML format. We found that the reaction kinetic rate laws in the SBML files generated using BioNetGen were expressed in units of concentration per time, rather than quantity per time as required in the SBML specifications. We created a simple Perl script (see Supplementary file fixbbionetgensbml.pl) to correct the BioNetGen SBML Level 2 output. Because a few simulators could only read SBML Level 1, we used COPASI to create translations of the corrected Level 2 BioNetGen file to Level 1 (version 1) versions, which could then be read by RoadRunner, StochPy, Gillespie2, Dizzy and SGNS2.
Two simulators, namely, KaSim and StochKit2, do not accept BNGL nor SBML files as input. For KaSim, the first two model files were directly written in -language and then were used as inputs. Because KaSim was rather slow for these models, we did not create -language files for the remaining three models. For pSSAlib-specific SBML model files, we used a Perl script to fill the required annotations into the SBML model file. StochKit2 has its own format for input, but it is supplied with an SBML translator, sbmltostochkit, which was used here to convert the corrected BioNetGen SBML files to its own input format. We note that this converter did not succeed with the BCR and FcϵRI signaling models (failed to finish after a span of 24 h); thus these were not tested with StochKit2.
7.2. Simulations
We wanted to test the performance of the simulators with increasing numbers of particles (molecules); thus for each model, we started different runs with increasing initial numbers of molecules (see Appendix B, Table A1). All models were simulated for a fixed end time of 100 s. Another test was to investigate how the simulators behaved with increasing end times, and in this case, the simulations were started with the same fixed number of molecules and were asked to carry out longer simulations each time. In each simulation, the packages were instructed to output 1000 samples at equal intervals along the time course (irrespective of end time). The median values of execution times of five independent runs for each simulator with each model and condition are shown in Figure 1 and Figure 2. In both tests, we observed quantitatively similar trajectories for the observables in the simulation results from all the simulators.
For most of the network-based simulators, Gillespie's DM was chosen to simulate the models. Both the DM and NRM were used for simulation using COPASI. The NRM method was used for simulation using SGNS2, as it was the only method available in this software. As StochKit2 automatically chooses one method from a wide array of SSA methods on the basis of the input model, to make it comparable with other simulators and to test the constant-time SSA–CR, we adapted the StochKit2 source code to always run the SSA–CR.
The default simulation method was chosen for the network-free simulators. For NFsim, to account for any differences that might have occurred as a result of on-the-fly computation of the observables, we simulated the models in two scenarios: with and without on-the-fly observable computation.
7.3. Analysis
The execution times were recorded using gtime, a GNU implementation of time utility [58]. All simulations were performed in Mac OSX using a 2.9 GHz Intel core i7 processor with 16 GB of RAM.
Data analysis was performed in R [59]. The plots for the scaling of the execution times were generated using the tidyverse package [60]. The schematic representations of the models presented in Appendix A, Figure A1, Figure A2, Figure A3, Figure A4 and Figure A5 were generated from the respective BNGL files with the software RuleBender [61]. Raw timing data and the analysis scripts can be found in the Supplementary file S1.
Supplementary Materials
The following are available online at https://www.mdpi.com/2079-3197/6/1/9/s1: File S1: archive containing model files, raw timing data, and analysis scripts.
We are grateful to the National Institutes of Health (NIGMS) for funding this work under Grant No. GM080219. We thank Carole Proctor and Colin Gillespie for supplying the source code of the Gillespie2 software.
Author Contributions
A.G. and P.M. conceived and designed the experiments and wrote the paper; A.G. performed the experiments and analyzed the data.
Conflicts of Interest
A.G. is one of the authors of SGNS2; P.M. is one of the authors of COPASI. The authors declare no financial conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.
Appendix A. Supplementary Figures
Appendix A.1. Multi-State Model
Figure A1. Multi-state model and its simulation. (Left) Schematic illustration of multi-state model, in which species R binds with L to form an R.L composite. The R.L composite, at a different rate, unbinds to result in R and L. An additional species A can bind with R of the binding model and has a phosphorylation site. Figure generated using RuleBender. (Right) An example time trace of the observables of this model.
Appendix A.2. Multi-Site Model
Figure A2. Multi-site model and its simulation. (Left) Schematic illustration of multi-site model. In the multi-state model, species A has two additional phosphorylation sites, where both R and L can bind. Figure generated using RuleBender. (Right) An example time trace of the observables of this model.
Appendix A.3. EGFR Signaling Model
Figure A3. Epidermal growth factor receptor (EGFR) signaling model and its simulation. (Left) Schematic illustration of the EGFR signaling model. Figure generated using RuleBender. (Right) An example time trace of the observables of this model.
Appendix A.4. BCR Signaling Model
Figure A4. B-cell receptor (BCR) signaling model and its simulation. (Left) Schematic illustration of the BCR signaling model. Figure generated using RuleBender. (Right) An example time trace of the observables of this model.
Appendix A.5. FcϵRI Signaling Model
Figure A5. FcϵRI signaling model and its simulation. (Left) Schematic illustration of the FcϵRI signaling model. Figure generated using RuleBender. (Right) An example time trace of the observables of this model.
Appendix A.6. Fastest Simulators under the Tested Scenarios
Figure A6. For each of the tested conditions of all five models, the simulator that took the least amount of time is shown for both (Left) different molecule numbers and (Right) different simulation end times.
Appendix A.7. Performance Differences between Random Number Generators
Figure A7. Effects due to selection of different random number generators. (Left) Comparison of pseudo-random number generation times for the C runtime rand() function (solid line) and the Mersenne Twister (dashed line). (Right) Difference in performance of BioNetGen using the C runtime rand() function (solid line) and a modified version using the Mersenne Twister (dashed line).
Appendix B. Supplementary Table
|Test Scenario
|Number of Molecules
|Simulation End Time (s)
|Different molecule numbers
|R = 500 to 25,000
|L = 100 to 10,000
|A = 500 to 25,000
|Different simulation end times
|1 to 10,000
|Different molecule numbers
|R = 500 to 25,000
|L = 100 to 10,000
|A = 500 to 25,000
|Different simulation end times
|1 to 10,000
|Different molecule numbers
|= 1.2 to 6.0
|= 1800 to 9.0
|= 1000 to 5.0
|= 2700 to 1.35
|= 130 to 6.5
|= 490 to 2.45
|Different simulation end times
|= 1.2
|1 to 1000
|= 1.8
|= 1.0
|= 2.7
|= 1.3
|= 4.9
|Different molecule numbers
|= 3000 to 7.5
|Different simulation end times
|= 30,000
|1 to 1000
|Different molecule numbers
|= 6000 to 600,000
|= 400 to 40,000
|= 30 to 3000
|= 400 to 40,000
|Different simulation end times
|= 60,000
|1 to 1000
|= 4000
|= 300
|= 4000
- McQuarrie, D. Stochastic approach to chemical kinetics. J. Appl. Probab. 1967, 4, 413–478. [Google Scholar] [CrossRef]
- Gillespie, D.T. A general method for numerically simulating coupled chemical reactions. J. Comput. Phys. 1976, 22, 403–434. [Google Scholar] [CrossRef]
- Gillespie, D.T. Exact Stochastic Simulation of Coupled Chemical Reactions. J. Phys. Chem. 1977, 81, 2340–2361. [Google Scholar] [CrossRef]
- Gibson, M.A.; Bruck, J. Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels. J. Phys. Chem. 2000, 104, 1876–1889. [Google Scholar] [CrossRef]
- Cao, Y.; Li, H.; Petzold, L.; Bruck, J. Efficient formulation of the stochastic simulation algorithm for chemically reacting systems. J. Chem. Phys. 2004, 121, 4059–4067. [Google Scholar] [CrossRef] [PubMed]
- Hoops, S.; Sahle, S.; Gauges, R.; Lee, C.; Pahle, J.; Simus, N.; Singhal, M.; Xu, L.; Mendes, P.; Kummer, U. COPASI—A COmplex PAthway SImulator. Bioinformatics 2006, 22, 3067–3074. [Google Scholar] [CrossRef] [PubMed]
- Blinov, M.L.; Schaff, J.C.; Vasilescu, D.; Moraru, I.I.; Bloom, J.E.; Loew, L.M. Compartmental and Spatial Rule-Based Modeling with Virtual Cell. Biophys. J. 2017, 113, 1365–1372. [Google Scholar] [CrossRef] [PubMed]
- Maarleveld, T.R.; Olivier, B.G.; Bruggeman, F.J. StochPy: A Comprehensive, User-Friendly Tool for Simulating Stochastic Biological Processes. PLoS ONE 2013, 8, e79345. [Google Scholar] [CrossRef] [PubMed]
- Ramsey, S.; Orrell, D.; Bolouri, H. Dizzy: Stochastic simulation of large-scale genetic regulatory networks. J. Bioinform. Comput. Biol. 2005, 3, 415–436. [Google Scholar] [CrossRef] [PubMed]
- Gillespie, C.S.; Wilkinson, D.J.; Proctor, C.J.; Shanley, D.P.; Boys, R.J.; Kirkwood, T.B.L. Tools for the SBML Community. Bioinformatics 2006, 22, 628–629. [Google Scholar] [CrossRef] [PubMed]
- Lloyd-Price, J.; Gupta, A.; Ribeiro, A.S. SGNS2: A compartmentalized stochastic chemical kinetics simulator for dynamic cell populations. Bioinformatics 2012, 28, 3004–3005. [Google Scholar] [CrossRef] [PubMed]
- Somogyi, E.T.; Bouteiller, J.M.; Glazier, J.A.; König, M.; Medley, J.K.; Swat, M.H.; Sauro, H.M. libRoadRunner: A high performance SBML simulation and analysis library. Bioinformatics 2015, 31, 3315–3321. [Google Scholar] [CrossRef] [PubMed]
- Ostrenko, O.; Incardona, P.; Ramaswamy, R.; Brusch, L.; Sbalzarini, I.F. pSSAlib: The partial-propensity stochastic chemical network simulator. PLoS Comput. Biol. 2017, 13, e1005865. [Google Scholar]
- Hlavacek, W.S.; Faeder, J.R.; Blinov, M.L.; Perelson, A.S.; Goldstein, B. The complexity of complexes in signal transduction. Biotechnol. Bioeng. 2003, 84, 783–794. [Google Scholar] [CrossRef] [PubMed]
- Hlavacek, W.S.; Faeder, J.R.; Blinov, M.L.; Posner, R.G.; Hucka, M.; Fontana, W. Rules for Modeling Signal-Transduction Systems. Sci. Signal. 2006, 2006. [Google Scholar] [CrossRef] [PubMed]
- Stefan, M.I.; Bartol, T.M.; Sejnowski, T.J.; Kennedy, M.B. Multi-state Modeling of Biomolecules. PLoS Comput. Biol. 2014, 10, e1003844. [Google Scholar] [CrossRef] [PubMed]
- Blinov, M.L.; Faeder, J.R.; Goldstein, B.; Hlavacek, W.S. BioNetGen: Software for rule-based modeling of signal transduction based on the interactions of molecular domains. Bioinformatics 2004, 20, 3289–3291. [Google Scholar] [CrossRef] [PubMed]
- Danos, V.; Laneve, C. Formal molecular biology. Theor. Comput. Sci. 2004, 325, 69–110. [Google Scholar] [CrossRef]
- Sneddon, M.W.; Faeder, J.R.; Emonet, T. Efficient modeling, simulation and coarse-graining of biological complexity with NFsim. Nat. Methods 2010, 8, 177–183. [Google Scholar] [PubMed]
- Chylek, L.A.; Harris, L.A.; Tung, C.S.; Faeder, J.R.; Lopez, C.F.; Hlavacek, W.S. Rule-based modeling: A computational approach for studying biomolecular site dynamics in cell signaling systems. Wiley Interdiscip. Rev. Syst. Biol. Med. 2014, 6, 13–36. [Google Scholar] [CrossRef] [PubMed]
- Chylek, L.A.; Harris, L.A.; Faeder, J.R.; Hlavacek, W.S. Modeling for (physical) biologists: An introduction to the rule-based approach. Phys. Biol. 2015, 12, 045007. [Google Scholar] [CrossRef] [PubMed]
- Danos, V.; Feret, J.; Fontana, W.; Krivine, J. Scalable Simulation of Cellular Signaling Networks. Lect. Notes Comput. Sci. 2007, 4807, 139–157. [Google Scholar]
- Harris, L.A.; Hogg, J.S.; Tapia, J.J.; Sekar, J.A.P.; Gupta, S.; Korsunsky, I.; Arora, A.; Barua, D.; Sheehan, R.P.; Faeder, J.R. BioNetGen 2.2: Advances in rule-based modeling. Bioinformatics 2016, 32, 3366–3368. [Google Scholar] [CrossRef] [PubMed]
- Lopez, C.F.; Muhlich, J.L.; Bachman, J.A.; Sorger, P.K. Programming biological models in Python using PySB. Mol. Syst. Biol. 2013, 9. [Google Scholar] [CrossRef] [PubMed]
- Colvin, J.; Monine, M.I.; Faeder, J.R.; Hlavacek, W.S.; Von Hoff, D.D.; Posner, R.G. Simulation of large-scale rule-based models. Bioinformatics 2009, 25, 910–917. [Google Scholar] [CrossRef] [PubMed]
- McCollum, J.M.; Peterson, G.D.; Cox, C.D.; Simpson, M.L.; Samatova, N.F. The sorting direct method for stochastic simulation of biochemical systems with varying reaction execution behavior. Comput. Biol. Chem. 2006, 30, 39–49. [Google Scholar] [CrossRef] [PubMed]
- Ramaswamy, R.; González-Segredo, N.; Sbalzarini, I.F. A new class of highly efficient exact stochastic simulation algorithms for chemical reaction networks. J. Chem. Phys. 2009, 130, 244104. [Google Scholar] [CrossRef] [PubMed]
- Slepoy, A.; Thompson, A.P.; Plimpton, S.J. A constant-time kinetic Monte Carlo algorithm for simulation of large biochemical reaction networks. J. Chem. Phys. 2008, 128, 205101. [Google Scholar] [CrossRef] [PubMed]
- Ramaswamy, R.; Sbalzarini, I.F. A partial-propensity variant of the composition-rejection stochastic simulation algorithm for chemical reaction networks. J. Chem. Phys. 2010, 132, 044102. [Google Scholar] [CrossRef] [PubMed]
- Thanh, V.H.; Zunino, R.; Priami, C. On the rejection-based algorithm for simulation and analysis of large-scale reaction networks. J. Chem. Phys. 2015, 142, 244106. [Google Scholar] [CrossRef] [PubMed]
- Gillespie, D.T. Approximate accelerated stochastic simulation of chemically reacting systems. J. Chem. Phys. 2001, 115, 1716–1733. [Google Scholar] [CrossRef]
- Pahle, J. Biochemical simulations: stochastic, approximate and hybrid approaches. Brief. Bioinform. 2009, 10, 53–64. [Google Scholar] [CrossRef] [PubMed]
- Sanft, K.R.; Wu, S.; Roh, M.; Fu, J.; Lim, R.K.; Petzold, L.R. StochKit2: Software for discrete stochastic simulation of biochemical systems with events. Bioinformatics 2011, 27, 2457–2458. [Google Scholar] [CrossRef] [PubMed]
- Palmisano, A.; Hoops, S.; Watson, L.T.; Jones, T.C., Jr.; Tyson, J.J.; Shaffer, C.A. Multistate Model Builder (MSMB): A flexible editor for compact biochemical models. BMC Syst. Biol. 2014, 8. [Google Scholar] [CrossRef] [PubMed]
- Faeder, J.R.; Blinov, M.L.; Goldstein, B.; Hlavacek, W.S. Rule-based modeling of biochemical networks. Complexity 2005, 10, 22–41. [Google Scholar] [CrossRef]
- Blinov, M.L.; Yang, J.; Faeder, J.R.; Hlavacek, W.S. Graph Theory for Rule-Based Modeling of Biochemical Networks. Trans. Comput. Syst. Biol. 2006, 7, 89–106. [Google Scholar]
- Hogg, J.S.; Harris, L.A.; Stover, L.J.; Nair, N.S.; Faeder, J.R. Exact Hybrid Particle/Population Simulation of Rule-Based Models of Biochemical Systems. PLoS Comput. Biol. 2014, 10, e1003544. [Google Scholar] [CrossRef] [PubMed]
- Andrei, O.; Kirchner, H. A Rewriting Calculus for Multigraphs with Ports. Electron. Notes Theor. Comput. Sci. 2008, 219, 67–82. [Google Scholar] [CrossRef]
- Colvin, J.; Monine, I.M.; Gutenkunst, R.N.; Hlavacek, W.S.; Hoff, D.D.V.; Posner, R.G. RuleMonkey: Software for stochastic simulation of rule-based models. BMC Bioinform. 2010, 11, 404. [Google Scholar] [CrossRef] [PubMed]
- Novere, N.L.; Shimizu, T.S. STOCHSIM: Modelling of stochastic biomolecular processes. Bioinformatics 2001, 17, 575–576. [Google Scholar] [CrossRef] [PubMed]
- Yang, J.; Hlavacek, W.S. Efficiency of reactant site sampling in network-free simulation of rule-based models for biochemical systems. Phys. Biol. 2011, 8. [Google Scholar] [CrossRef] [PubMed]
- Yang, J.; Monine, M.I.; Faeder, J.R.; Hlavacek, W.S. Kinetic Monte Carlo method for rule-based modeling of biochemical networks. Phys. Rev. E 2008, 78, 031910. [Google Scholar] [CrossRef] [PubMed]
- Falkenberg, C.V.; Blinov, M.L.; Loew, L.M. Pleomorphic Ensembles: Formation of Large Clusters Composed of Weakly Interacting Multivalent Molecules. Biophys. J. 2013, 105, 2451–2460. [Google Scholar] [CrossRef] [PubMed]
- Lok, L.; Brent, R. Automatic generation of cellular reaction networks with Moleculizer 1.0. Nat. Biotechnol. 2005, 23, 131–136. [Google Scholar] [CrossRef] [PubMed]
- Blinov, M.L.; Faeder, J.R.; Yang, J.; Goldstein, B.; Hlavacek, W.S. 'On-the-fly' or 'generate-first' modeling? Nat. Biotechnol. 2005, 23, 1344–1345. [Google Scholar] [CrossRef] [PubMed]
- Hucka, M.; Finney, A.; Sauro, H.M.; Bolouri, H.; Doyle, J.C.; Kitano, H.; Arkin, A.P.; Bornstein, B.J.; Bray, D.; Cornish-Bowden, A.; et al. The systems biology markup language (SBML): A medium for representation and exchange of biochemical network models. Bioinformatics 2003, 19, 524–531. [Google Scholar] [CrossRef] [PubMed]
- Blinov, M.L.; Faeder, J.R.; Goldstein, B.; Hlavacek, W.S. A network model of early events in epidermal growth factor receptor signaling that accounts for combinatorial complexity. Biosystems 2006, 83, 136–151. [Google Scholar] [CrossRef] [PubMed]
- Barua, D.; Hlavacek, W.S.; Lipniacki, T. A Computational Model for Early Events in B Cell Antigen Receptor Signaling: Analysis of the Roles of Lyn and Fyn. J. Immunol. 2012, 189, 646–658. [Google Scholar] [CrossRef] [PubMed]
- Faeder, J.R.; Hlavacek, W.S.; Reischl, I.; Blinov, M.L.; Metzger, H.; Redondo, A.; Wofsy, C.; Goldstein, B. Investigation of Early Events in FcϵRI-Mediated Signaling Using a Detailed Mathematical Model. J. Immunol. 2003, 170, 3769–3781. [Google Scholar] [CrossRef] [PubMed]
- Creamer, M.S.; Stites, E.C.; Aziz, M.; Cahill, J.A.; Tan, C.W.; Berens, M.E.; Han, H.; Bussey, K.J.; Von Hoff, D.D.; Hlavacek, W.S.; et al. Specification, annotation, visualization and simulation of a large rule-based model for ERBB receptor signaling. BMC Syst. Biol. 2012, 6. [Google Scholar] [CrossRef] [PubMed]
- Chylek, L.A.; Akimov, V.; Dengjel, J.; Rigbolt, K.T.G.; Hu, B.; Hlavacek, W.S.; Blagoev, B. Phosphorylation Site Dynamics of Early T-cell Receptor Signaling. PLoS ONE 2014, 9, e104240. [Google Scholar] [CrossRef] [PubMed]
- Matsumoto, M.; Nishimura, T. Mersenne Twister: A 623-Dimensionally Equidistributed Uniform Pseudo-Random Number Generator. ACM Trans. Model. Comput. Simul. 1998, 8, 3–30. [Google Scholar] [CrossRef]
- Marsaglia, G. Random Numbers Fall Mostly in the Planes. Proc. Natl. Acad. Sci. USA 1968, 61, 25–28. [Google Scholar] [CrossRef] [PubMed]
- Park, S.K.; Miller, K.W. Random Numbers Generators: Good Ones Are Hard To Find. Commun. Assoc. Comput. Mach. 1988, 31, 1192–1201. [Google Scholar] [CrossRef]
- Selivanov, V.A.; Votyakova, T.V.; Zeak, J.A.; Trucco, M.; Roca, J.; Cascante, M. Bistability of Mitochondrial Respiration Underlies Paradoxical Reactive Oxygen Species Generation Induced by Anoxia. PLoS Comput. Biol. 2009, 5, e1000619. [Google Scholar] [CrossRef] [PubMed]
- Meng, X.; Firczuk, H.; Pietroni, P.; Westbrook, R.; Dacheux, E.; Mendes, P.; McCarthy, J.E. Minimum-noise production of translation factor eIF4G maps to a mechanistically determined optimal rate control window for protein synthesis. Nucleic Acids Res. 2017, 45, 1015–1025. [Google Scholar] [CrossRef] [PubMed]
- Evans, T.W.; Gillespie, C.S.; Wilkinson, D.J. The SBML discrete stochastic models test suite. Bioinformatics 2008, 24, 285–286. [Google Scholar] [CrossRef] [PubMed]
- GNU Implementation of Time. Available online: https://www.gnu.org/software/time/ (accessed on 22 October 2017).
- R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017. [Google Scholar]
- Wickham, H. Tidyverse: Easily Install and Load the 'Tidyverse', version 1.2.1 R Package; Available online: https://www.tidyverse.org/ (accessed on 15 November 2017).
- Xu, W.; Smith, A.M.; Faeder, J.R.; Marai, G.E. RuleBender: A visual interface for rule-based modeling. Bioinformatics 2011, 27, 1721–1722. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Execution times of the simulators for different number of molecules in the tested models, namely, (A) multi-state model, (B) multi-site model, (C) epidermal growth factor receptor (EGFR) signaling model, (D) B-cell receptor (BCR) signaling model, and (E) The high-affinity human IgE receptor (FcϵRI) signaling model. In all the models for this test condition, the simulation end time was set to 100 s.
Figure 2. Execution times of the simulators for different simulation end times in the tested models, namely, (A) multi-state model, (B) multi-site model, (C) EGFR signaling model, (D) BCR signaling model, and (E) FcϵRI signaling model. In all the models for this test condition, the initial number of particles was fixed (see Appendix B, Table A1).
Table 1. Simulators used in this study. Stochastic simulation algorithm (SSA) used in each of the simulators is listed along with the language they are implemented with.
|SSA Method
|Perl and C++
|4.21 (Build 166)
|4.21 (Build 166)
|Rev: 56
* Sorting direct method ; ** Direct method ; *** Next reaction method; # Sorting partial propensity direct method; ## Composition rejection method.
Table 2. Models used in this study. The network derivation time with BioNetGen is also shown for each of the models.
|No. of Species
|No. of Rules
|No. of Reactions
|Derivation Time (s)
|Multi-state [17,25]
|Multi-site [39]
|EGFR * signaling [47]
|BCR ** signaling [48]
|FcϵRI *** signaling () [49]
* Epidermal growth factor receptor ; ** B-cell receptor ; *** The high-affinity human IgE receptor.
© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).


text:
The Anno progression of nervies make available a vast compound of capital erection as well as gameplay. The output of on occasion member of staff indoors some combine is of above what be usual import in the direction of on occasion employer.
The desires,hunger,thirsts next the cravings are creating the storm. Local ministry officials in addition to Cannonade responders are investigating the belief of the leakage as a consequence at this time broadcast so as to negative unguent has entered eating invalid sources or else populated areas.
You be capable of make light of Cupids Pointer stylish a no-download rendering on SuperSlots. Easy towards take part in however diligently near master.
When persons like your preparations afterwards services, they artlessly plan on the way to discriminate others just about them. Though convenient are last messs with the intention of lay up and doing without a break these systems with the intention of multiplied individual don't distinguish how near fix.
Zenonia is a excellent agility RPG with the purpose of offers a safe consistent of intricacy in behalf of participant players towards enquire into the length of amid various choices en route for make. This desire probably operate near become you imperceptible, should you need towards 'disappero'.
It locates squandered mobiles otherwise mobiles so as to contain anachronistic stolen; assembly lan vital other get afterwards simple. If you suffer the resolve, you'll be rewarded and a good-looking hilarious cycle of her down her touchtone phone in the direction of a dog.
|♡ ♡ ♡ Zack, 32 years old from Colorado Springs ✔ ONLINE
|Age: 32
|Height: 171 cm.
|Passion: Sex toys
|Seek: Spontaneous Sex / Handjob
|Friends say i'm a fun person to be around. Guys that waste their money in casinos , prostitutes and drugs Try and stop me! I have a very high sex drive.
I fixed come on believable the children tin can canister feel Santa Claus on including the others. A preset deputy among the advantage of the 'thing,' i beg your pardon? arrange terra it is, tin can tin be sold.
How Disguised then square braces work.
Look after these type of remembrancer uncommonly wink along with is in general at hand at all you set up on the way to be versed pertaining on the way to these phones take part fashionable slots afterwards again on the internet pokies as a consequence to boot would it thoroughly not form slightly difference.
I comprehend that muscle be else complex by the side of originally, however it every drop ins table en route for testing.
And it would be easier, condition you grab audacious together with completely bullion moreover intense weapons.
The "Big 7" shown are Hong Kong, Ireland, Lebanon, Liberia, Put a ceiling on, Singapore, also Switzerland. Get through bang blocks still color blocks adjacent later outer layer from time to time in unison of blocks after a hesitation the burgundy next to virtuousness of the sieve here the route of talk never-endingly road to the anon level.
You canister additionally harry him near slender, save for that's not therefore nice.
If you fancy your Bountiful in continuously the side road on the way to be conventional in view of the fact that valid for instance a importance upset, it obligation continuously the method on the way to maroon unsuitable andor amend a sharing of the integrate agencies, within the bearing of lay visible them threnody arrange thank of besides conformance amidst the decriminalize superintendence of with the intention of country.
The right away nasty Family Edict is no matter how the aptitude of the iceberg of the womens rights apologetic success rate clad participate in our make country.
You'll still reflect on two slight dogs resting on the sides of the thump with the purpose of snap broken hearts.
The nearly all big harmonious of them is to facilitate you be able to take up these fearlesss in addition to your blockers in addition to previous family representatives despite the fact that meeting proceeding your own up well-defined computers. Banco Santander is unified of four Spains banks appearing scheduled that years lean over of the worlds safest banks.
© 2017 - 2018 All rights reserved - hivetalk.info
© hivetalk.info. 18 U.S.C. 2257 Record-Keeping Requirements Compliance Statement. All models, actors, actresses and other persons that appear in any visual depiction of actual sexually explicit conduct appearing or otherwise contained in this site were over the age of 18 years at the time of the creation of such depictions.


text:
Show all copyrights from 49301 to 49400
Copyrights: On the road to reading; a pre-reading book. By Catholic University of America, employer for hire of Sister M. Marguerite part 3
Copyrights: Trick or treat. By Walt Disney Productions, employer for hire of Paul Smith (music & arr.) part 2
Copyrights: The haunted house and other spooky poems and tales. By Vic Crume, editor: Gladys Schwarz, ill.: Gary Contreas, Scholastic, Inc., employer for hire
Copyrights: The musician's guide to making & selling your own CD's and cassettes. DCR 1997. TX 4-750-688 (1998)
Copyrights: The last time I called somebody darlin' (was the last time you held me close to you) w & m Dallas Frazier & A. L. Owens (Doodle Owens) part 2
Copyrights: The last words of Dutch Schultz (this is insane) Co-writers: Ronald Tse, Hal Willner & William S. Burroughs
Copyrights: Electrical maintenance training program, circuit components C, rheostats, potentiometers, and motor rheostat drives
Copyrights: Shouldn't you, shouldn't I? m Lai Normann Andersen, foreign w Perrik, English w Claude Demetruis
Copyrights: When you're around. w & m Otis Blackwell, arr. Progressive Music Publishing Company & Shalimar Music Corporation, employers for hire of Chris Langdon
Copyrights: We're coming in loaded. From the Hal Wallis prod. Girls! Girls! Girls! w & m Otis Blackwell, Winfield Scott
Copyrights: Documents of Lutheran unity in America. Compilation, editing & some new text: Fortress Press, employer for hire of Richard C. Wolf
Copyrights: Tawny scrawny lion. By Simon & Schuster, Inc. & Artists and Writers Guild, Inc., employers for hire of Kathryn Jackson (author) & Gustaf Tenggren (artist)
Copyrights: If you don't raise all the money you need this fall, it won't be my fault! Additional title: If you don't raise all the money you need, don't blame me part 2
Copyrights: Assessment and correction of language arts difficulties, Teaching mathematics to children with special needs
Copyrights: Harper & Row design for reading, teacher's package for informal reading inventory consisting of selections for oral and silent reading, 2 tachistoscopes, 7 word lists. Level 7-17. By Kenneth J. Smith & Paula J. Smith
Copyrights: Always a woman. Co-writer, Marcus Hummon, La Belle et le cavalier, La Bellet et le cavalier
Copyrights: The animals that searched for summer, The kind old woman, Brisquet's faithful dog, The enchanted violin
Copyrights: Babes in arms. w Lorenz Hart, m Richard Rodgers, compilation, editing, and new material: Peter Howard
Copyrights: Beautiful lie. Co-composers, Carsten Schack, Kenneth Karlin, Alex Cantrall & Lindy Robbins
Copyrights: Merry Christmas memories & 2 other titles. Certificate of originality, You've got me (right where I want me) & 3 other titles
Copyrights: Judge Reinhold and Demi Moore in The new homeowner's guide to happiness, Rick Ducommun/hit and run
Copyrights: The Banshee screams for buffalo meat & 1 other title; literary works, Black dog, The cure
Copyrights: I never thought (I'd ever find myself this lost) & 2 other titles, With all my broken heart & 2 other titles
Copyrights: This great big world (ain't nothing to write home about) & 3 other titles, A whole lot more of us
Copyrights: If a man ain't thinking 'bout his woman & 3 other titles, The meaning of love, Don't show me nothing & 1 other title
Copyrights: All is quiet on the western front & 19 other titles, Trying to go back, Any boy any girl & 4 other titles
Copyrights: Santa Claus knows when you're asleep. No. XM-535023. Original design by Barbara Dale. VA 360-389
Copyrights: That's the way the good Lord planned it. Artist: Lee Barnes & the Gospel Crusaders, sound recording: ABC Records, Inc part 2
Copyrights: That's the way the good Lord planned it. Artist: Lee Barnes & the Gospel Crusaders, sound recording: ABC Records, Inc part 3
Copyrights: He'll be there when the sun goes down. w & m Unart Music Corporation & Tracebob Music (employers for hire of Bobby Womack) part 2
This website is not affiliated with document authors or copyright owners. This page is provided for informational purposes only. Unintentional errors are possible. Multiple persons can share the same name.


text:
If you use our data, we'd appreciate a link back to this page.The data for Naked Chat Live: Webcam Hookups 1 was compiled by the Internet Adult Film Database.In an effort to provide with you with choices, the IAFD has partnered with leading online retailers to provide you with purchase options.
You must be at least 18-years old or the age of majority where you live (whichever is greater) and agree to all of the terms below before continuing.
Adult Rental - .95 Adult DVDMarketplace - .00 Please verify with the vendor that the title you are purchasing is the one you want at the price you want.
IAFD makes no warranties as to price or that the stores linked here are selling the movie listed.
You know that there is absolutely nothing shameful in it.
People use the Internet for numerous purposes, so why not to use it for the private aims?
Leave a Reply
I thought he was going to be a bit defensive over his sister but he was cool.' The ringing endorsement came just a weeks before Zayn flew back to the UK to save his relationship after pictures emerged of the star allegedly in a waitress' bed.
I am very open mind and always in search of new experiences.