BioSimulators: a central registry of simulation engines and services for recommending specific tools.

Computational models have great potential to accelerate bioscience, bioengineering, and medicine. However, it remains challenging to reproduce and reuse simulations, in part, because the numerous formats and methods for simulating various subsystems and scales remain siloed by different software tools. For example, each tool must be executed through a distinct interface. To help investigators find and use simulation tools, we developed BioSimulators (https://biosimulators.org), a central registry of the capabilities of simulation tools and consistent Python, command-line and containerized interfaces to each version of each tool. The foundation of BioSimulators is standards, such as CellML, SBML, SED-ML and the COMBINE archive format, and validation tools for simulation projects and simulation tools that ensure these standards are used consistently. To help modelers find tools for particular projects, we have also used the registry to develop recommendation services. We anticipate that BioSimulators will help modelers exchange, reproduce, and combine simulations.

Epidemiological modeling in StochSS Live!

SUMMARY: We present StochSS Live!, a web-based service for modeling, simulation and analysis of a wide range of mathematical, biological and biochemical systems. Using an epidemiological model of COVID-19, we demonstrate the power of StochSS Live! to enable researchers to quickly develop a deterministic or a discrete stochastic model, infer its parameters and analyze the results. AVAILABILITY AND IMPLEMENTATION: StochSS Live! is freely available at https://live.stochss.org/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Phylogenomics and historical biogeography of West Indian Rock Iguanas (genus Cyclura).

The genus Cyclura includes nine extant species and six subspecies of West Indian Rock Iguanas and is one of the most imperiled genera of squamate reptiles globally. An understanding of species diversity, evolutionary relationships, diversification, and historical biogeography in this group is crucial for implementing sound long-term conservation strategies. We collected DNA samples from 1 to 10 individuals per taxon from all Cyclura taxa (n = 70 ingroup individuals), focusing where possible on incorporating individuals from different populations of each species. We also collected 1-2 individuals from each of seven outgroup species of iguanas (Iguana delicatissima; five Ctenosaura species) and Anolis sagrei (n = 12 total outgroup individuals). We used targeted genomic sequence capture to isolate and to sequence 1,872 loci comprising of 687,308 base pairs (bp) from each of the 82 individuals from across the nuclear genome. We extracted mitochondrial reads and assembled and annotated mitogenomes for all Cyclura taxa plus outgroup species. We present well-supported phylogenomic gene tree/species tree analyses for all extant species of Cyclura using ASTRAL-III, SVDQuartets, and StarBEAST2 methods, and discuss the taxonomic, biogeographic, and conservation implications of these data. We find a most recent common ancestor of the genus 9.91 million years ago. The earliest divergence within Cyclura separates C. pinguis from a clade comprising all other Cyclura. Within the latter group, a clade comprising C. carinata from the southern Lucayan Islands and C. ricordii from Hispaniola is the sister taxon to a clade comprising the other Cyclura. Among the other Cyclura, the species C. cornuta and C. stejnegeri (from Hispaniola and Isla Mona) form the sister taxon to a clade of species from Jamaica (C. collei), Cuba and Cayman Islands (C. nubila and C. lewisi), and the eastern (C. rileyi) and western (C. cychlura) Lucayan Islands. Cyclura cychlura and C. rileyi form a clade whose sister taxa are C. nubila and C. lewisi. Cyclura collei is the sister taxon to these four species combined.

Coordinating cell polarization and morphogenesis through mechanical feedback.

Many cellular processes require cell polarization to be maintained as the cell changes shape, grows or moves. Without feedback mechanisms relaying information about cell shape to the polarity molecular machinery, the coordination between cell polarization and morphogenesis, movement or growth would not be possible. Here we theoretically and computationally study the role of a genetically-encoded mechanical feedback (in the Cell Wall Integrity pathway) as a potential coordination mechanism between cell morphogenesis and polarity during budding yeast mating projection growth. We developed a coarse-grained continuum description of the coupled dynamics of cell polarization and morphogenesis as well as 3D stochastic simulations of the molecular polarization machinery in the evolving cell shape. Both theoretical approaches show that in the absence of mechanical feedback (or in the presence of weak feedback), cell polarity cannot be maintained at the projection tip during growth, with the polarization cap wandering off the projection tip, arresting morphogenesis. In contrast, for mechanical feedback strengths above a threshold, cells can robustly maintain cell polarization at the tip and simultaneously sustain mating projection growth. These results indicate that the mechanical feedback encoded in the Cell Wall Integrity pathway can provide important positional information to the molecular machinery in the cell, thereby enabling the coordination of cell polarization and morphogenesis.

An arbitrary Lagrangian Eulerian smoothed particle hydrodynamics (ALE-SPH) method with a boundary volume fraction formulation for fluid-structure interaction.

We present a new weakly-compressible smoothed particle hydrodynamics (SPH) method capable of modeling non-slip fixed and moving wall boundary conditions. The formulation combines a boundary volume fraction (BVF) wall approach with the transport-velocity SPH method. The resulting method, named SPH-BVF, offers detection of arbitrarily shaped solid walls on-the-fly, with small computational overhead due to its local formulation. This simple framework is capable of solving problems that are difficult or infeasible for standard SPH, namely flows subject to large shear stresses or at moderate Reynolds numbers, and mass transfer in deformable boundaries. In addition, the method extends the transport-velocity formulation to reaction-diffusion transport of mass in Newtonian fluids and linear elastic solids, which is common in biological structures. Taken together, the SPH-BVF method provides a good balance of simplicity and versatility, while avoiding some of the standard obstacles associated with SPH: particle penetration at the boundaries, tension instabilities and anisotropic particle alignments, that hamper SPH from being applied to complex problems such as fluid-structure interaction in a biological system.

The effect of cell geometry on polarization in budding yeast.

The localization (or polarization) of proteins on the membrane during the mating of budding yeast (Saccharomyces cerevisiae) is an important model system for understanding simple pattern formation within cells. While there are many existing mathematical models of polarization, for both budding and mating, there are still many aspects of this process that are not well understood. In this paper we set out to elucidate the effect that the geometry of the cell can have on the dynamics of certain models of polarization. Specifically, we look at several spatial stochastic models of Cdc42 polarization that have been adapted from published models, on a variety of tip-shaped geometries, to replicate the shape change that occurs during the growth of the mating projection. We show here that there is a complex interplay between the dynamics of polarization and the shape of the cell. Our results show that while models of polarization can generate a stable polarization cap, its localization at the tip of mating projections is unstable, with the polarization cap drifting away from the tip of the projection in a geometry dependent manner. We also compare predictions from our computational results to experiments that observe cells with projections of varying lengths, and track the stability of the polarization cap. Lastly, we examine one model of actin polarization and show that it is unlikely, at least for the models studied here, that actin dynamics and vesicle traffic are able to overcome this effect of geometry.

Macromolecular Crowding Regulates the Gene Expression Profile by Limiting Diffusion.

We seek to elucidate the role of macromolecular crowding in transcription and translation. It is well known that stochasticity in gene expression can lead to differential gene expression and heterogeneity in a cell population. Recent experimental observations by Tan et al. have improved our understanding of the functional role of macromolecular crowding. It can be inferred from their observations that macromolecular crowding can lead to robustness in gene expression, resulting in a more homogeneous cell population. We introduce a spatial stochastic model to provide insight into this process. Our results show that macromolecular crowding reduces noise (as measured by the kurtosis of the mRNA distribution) in a cell population by limiting the diffusion of transcription factors (i.e. removing the unstable intermediate states), and that crowding by large molecules reduces noise more efficiently than crowding by small molecules. Finally, our simulation results provide evidence that the local variation in chromatin density as well as the total volume exclusion of the chromatin in the nucleus can induce a homogenous cell population.

Stochastic Simulation Service: Bridging the Gap between the Computational Expert and the Biologist.

We present StochSS: Stochastic Simulation as a Service, an integrated development environment for modeling and simulation of both deterministic and discrete stochastic biochemical systems in up to three dimensions. An easy to use graphical user interface enables researchers to quickly develop and simulate a biological model on a desktop or laptop, which can then be expanded to incorporate increasing levels of complexity. StochSS features state-of-the-art simulation engines. As the demand for computational power increases, StochSS can seamlessly scale computing resources in the cloud. In addition, StochSS can be deployed as a multi-user software environment where collaborators share computational resources and exchange models via a public model repository. We demonstrate the capabilities and ease of use of StochSS with an example of model development and simulation at increasing levels of complexity.

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.

Spatial stochastic dynamics enable robust cell polarization.

Although cell polarity is an essential feature of living cells, it is far from being well-understood. Using a combination of computational modeling and biological experiments we closely examine an important prototype of cell polarity: the pheromone-induced formation of the yeast polarisome. Focusing on the role of noise and spatial heterogeneity, we develop and investigate two mechanistic spatial models of polarisome formation, one deterministic and the other stochastic, and compare the contrasting predictions of these two models against experimental phenotypes of wild-type and mutant cells. We find that the stochastic model can more robustly reproduce two fundamental characteristics observed in wild-type cells: a highly polarized phenotype via a mechanism that we refer to as spatial stochastic amplification, and the ability of the polarisome to track a moving pheromone input. Moreover, we find that only the stochastic model can simultaneously reproduce these characteristics of the wild-type phenotype and the multi-polarisome phenotype of a deletion mutant of the scaffolding protein Spa2. Significantly, our analysis also demonstrates that higher levels of stochastic noise results in increased robustness of polarization to parameter variation. Furthermore, our work suggests a novel role for a polarisome protein in the stabilization of actin cables. These findings elucidate the intricate role of spatial stochastic effects in cell polarity, giving support to a cellular model where noise and spatial heterogeneity combine to achieve robust biological function.

GillesPy2: A Biochemical Modeling Framework for Simulation Driven Biological Discovery.

Stochastic modeling has become an essential tool for studying biochemical reaction networks. There is a growing need for user-friendly and feature-complete software for model design and simulation. To address this need, we present GillesPy2, an open-source framework for building and simulating mathematical and biochemical models. GillesPy2, a major upgrade from the original GillesPy package, is now a stand-alone Python 3 package. GillesPy2 offers an intuitive interface for robust and reproducible model creation, facilitating rapid and iterative development. In addition to expediting the model creation process, GillesPy2 offers efficient algorithms to simulate stochastic, deterministic, and hybrid stochastic-deterministic models.

A positive feedback loop involving the Spa2 SHD domain contributes to focal polarization.

Focal polarization is necessary for finely arranged cell-cell interactions. The yeast mating projection, with its punctate polarisome, is a good model system for this process. We explored the critical role of the polarisome scaffold protein Spa2 during yeast mating with a hypothesis motivated by mathematical modeling and tested by in vivo experiments. Our simulations predicted that two positive feedback loops generate focal polarization, including a novel feedback pathway involving the N-terminal domain of Spa2. We characterized the latter using loss-of-function and gain-of-function mutants. The N-terminal region contains a Spa2 Homology Domain (SHD) which is conserved from yeast to humans, and when mutated largely reproduced the spa2Δ phenotype. Our work together with published data show that the SHD domain recruits Msb3/4 that stimulates Sec4-mediated transport of Bud6 to the polarisome. There, Bud6 activates Bni1-catalyzed actin cable formation, recruiting more Spa2 and completing the positive feedback loop. We demonstrate that disrupting this loop at any point results in morphological defects. Gain-of-function perturbations partially restored focal polarization in a spa2 loss-of-function mutant without restoring localization of upstream components, thus supporting the pathway order. Thus, we have collected data consistent with a novel positive feedback loop that contributes to focal polarization during pheromone-induced polarization in yeast.

Saving the Devils Is in the Details: Tasmanian Devil Facial Tumor Disease Can Be Eliminated with Interventions.

Tasmanian Devils facial tumor disease (DFTD) is severely impacting the population of this wild animal. We developed a computational model of the population of Tasmanian Devils, and the change induced by DFTD. We use this model to test possible intervention strategies Tasmanian conservationists could do. We investigate bait drop vaccination programs, diseased animal removals programs, and evolution of natural immunity. We conclude that a combination of intervention strategies gives the most favorable outcome. An additional goal of this paper is reproducibility of our results. Our StochSS software platform features the ability to share and reproduce the computational notebooks that created all of the results in the paper. We endeavor that all readers should be able to reproduce our results with minimum effort.

Noise resistant synchronization and collective rhythm switching in a model of animal group locomotion.

Biology is suffused with rhythmic behaviour, and interacting biological oscillators often synchronize their rhythms with one another. Colonies of some ant species are able to synchronize their activity to fall into coherent bursts, but models of this phenomenon have neglected the potential effects of intrinsic noise and interspecific differences in individual-level behaviour. We investigated the individual and collective activity patterns of two Leptothorax ant species. We show that in one species (Leptothorax sp. W), ants converge onto rhythmic cycles of synchronized collective activity with a period of about 20 min. A second species (Leptothorax crassipilis) exhibits more complex collective dynamics, where dominant collective cycle periods range from 16 min to 2.8 h. Recordings that last 35 h reveal that, in both species, the same colony can exhibit multiple oscillation frequencies. We observe that workers of both species can be stimulated by nest-mates to become active after a refractory resting period, but the durations of refractory periods differ between the species and can be highly variable. We model the emergence of synchronized rhythms using an agent-based model informed by our empirical data. This simple model successfully generates synchronized group oscillations despite the addition of noise to ants' refractory periods. We also find that adding noise reduces the likelihood that the model will spontaneously switch between distinct collective cycle frequencies.

Progressive arm muscle weakness in ALS follows the same sequence regardless of onset site: use of TOMS, a novel analytic method to track limb strength.

Objective: Examine sequence of weakness in arm muscles from longitudinal hand-held dynamometry (HHD) data in ALS for congruence with contiguous spread of neurodegeneration along spinal cord segments. Methods: Longitudinal HHD data from the Ceftriaxone clinical trial were examined using nonlinear mixed models, assuming a logistic trajectory from normal to zero strength. Unobserved baseline normal strength of weak muscles was assumed using strength of the best-preserved muscle. A novel metric called "time from onset to midway strength" (TOMS) was estimated for each muscle group, and TOMS ratios were examined to identify sequence of weakness, overall and by onset site. Results: Shoulder flexion (SF), elbow flexion (EF), elbow extension (EE), wrist extension (WE), and first dorsal interosseous (FDI) were measured on each side. Over a median of 36 weeks, 513 subjects provided 2589 sets of HHD measures. TOMS increased sequentially in the following order: FDI, WE, SF, EF, and EE. TOMS ratios estimates with 95% CIs (adjusted for multiple comparisons) were: WE/FDI 1.32 (1.24-1.41), SF/WE 1.06 (1.01-1.10), EF/SF 1.06 (1.02-1.10), and EE/EF 1.18 (1.12-1.23). Elbow and shoulder flexors weakened sooner than did elbow extensors. The sequence of arm muscle weakness progression was similar regardless of onset site. Conclusion: Nonsegmental progression of arm muscle weakness that is similar for different onset sites favors cortical influence/network spread over contiguous spread of neurodegeneration in the spinal cord. Furthermore, this study confirms the "split elbow" pattern. TOMS and other proposed methods may have value as outcome measures in clinical research.

The diffusive finite state projection algorithm for efficient simulation of the stochastic reaction-diffusion master equation.

We have developed a computational framework for accurate and efficient simulation of stochastic spatially inhomogeneous biochemical systems. The new computational method employs a fractional step hybrid strategy. A novel formulation of the finite state projection (FSP) method, called the diffusive FSP method, is introduced for the efficient and accurate simulation of diffusive transport. Reactions are handled by the stochastic simulation algorithm.

Validation data for a hybrid smoothed dissipative particle dynamics (SDPD) spatial stochastic simulation algorithm (sSSA) method.

We present the validation of the hybrid sSSA-SDPD method for advection-diffusion-reaction problems coupled to discrete biochemical systems, as presented in the publication "A hybrid smoothed dissipative particle dynamics (SDPD) spatial stochastic simulation algorithm (sSSA) for advection-diffusion-reaction problems" (Drawert et al., 2019). We validate 1D diffusion, and 2D diffusion cases against their analytical solutions. We present graphs and tables of data showing the error in the simulation method.

Sources of intraspecific variation in the collective tempo and synchrony of ant societies.

Populations of independently oscillating agents can sometimes synchronize. In the context of animal societies, conspicuous synchronization of activity is known in some social insects. However, the causes of variation in synchrony within and between species have received little attention. We repeatedly assessed the short-term activity cycle of ant colonies (Temnothorax rugatulus) and monitored the movements of individual workers and queens within nests. We detected persistent differences between colonies in the waveform properties of their collective activity oscillations, with some colonies consistently oscillating much more erratically than others. We further demonstrate that colony crowding reduces the rhythmicity (i.e., the consistent timing) of oscillations. Workers in both erratic and rhythmic colonies spend less time active than completely isolated workers, but workers in erratic colonies oscillate out of phase with one another. We further show that the queen's absence can impair the ability of colonies to synchronize worker activity and that behavioral differences between queens are linked with the waveform properties of their societies.

Using stochastic epidemiological models to evaluate conservation strategies for endangered amphibians.

Recent outbreaks of chytridiomycosis, the disease of amphibians caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), have contributed to population declines of numerous amphibian species worldwide. The devastating impacts of this disease have led researchers to attempt drastic conservation measures to prevent further extinctions and loss of biodiversity. The conservation measures can be labour-intensive or expensive, and in many cases have been unsuccessful. We developed a mathematical model of Bd outbreaks that includes the effects of demographic stochasticity and within-host fungal load dynamics. We investigated the impacts of one-time treatment conservation strategies during the disease outbreak that occurs following the initial arrival of Bd into a previously uninfected frog population. We found that for all versions of the model, for a large fraction of parameter space, none of the one-time treatment strategies are effective at preventing disease-induced extinction of the amphibian population. Of the strategies considered, treating frogs with antifungal agents to reduce their fungal load had the greatest likelihood of a beneficial outcome and the lowest risk of decreasing the persistence of the frog population, suggesting that this disease mitigation strategy should be prioritized over disinfecting the environment or reducing host density.

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 and 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.