Surface reaction-diffusion kinetics on lattice at the microscopic scale

Engineering a Novel Probiotic Toolkit in Escherichia coli Nissle 1917 for Sensing and Mitigating Gut Inflammatory Diseases

Inflammatory bowel disease (IBD) is characterized by chronic intestinal inflammation with no cure and limited treatment options that often have systemic side effects. In this study, we developed a target-specific system to potentially treat IBD by engineering the probiotic bacterium Escherichia coli Nissle 1917 (EcN). Our modular system comprises three components: a transcription factor-based sensor (NorR) capable of detecting the inflammation biomarker nitric oxide (NO), a type 1 hemolysin secretion system, and a therapeutic cargo consisting of a library of humanized anti-TNFα nanobodies. Despite a reduction in sensitivity, our system demonstrated a concentration-dependent response to NO, successfully secreting functional nanobodies with binding affinities comparable to the commonly used drug Adalimumab, as confirmed by enzyme-linked immunosorbent assay and in vitro assays. This newly validated nanobody library expands EcN therapeutic capabilities. The adopted secretion system, also characterized for the first time in EcN, can be further adapted as a platform for screening and purifying proteins of interest. Additionally, we provided a mathematical framework to assess critical parameters in engineering probiotic systems, including the production and diffusion of relevant molecules, bacterial colonization rates, and particle interactions. This integrated approach expands the synthetic biology toolbox for EcN-based therapies, providing novel parts, circuits, and a model for tunable responses at inflammatory hotspots.

Stochastic simulation algorithms for Interacting Particle Systems

Interacting Particle Systems (IPSs) are used to model spatio-temporal stochastic systems in many disparate areas of science. We design an algorithmic framework that reduces IPS simulation to simulation of well-mixed Chemical Reaction Networks (CRNs). This framework minimizes the number of associated reaction channels and decouples the computational cost of the simulations from the size of the lattice. Decoupling allows our software to make use of a wide class of techniques typically reserved for well-mixed CRNs. We implement the direct stochastic simulation algorithm in the open source programming language Julia. We also apply our algorithms to several complex spatial stochastic phenomena. including a rock-paper-scissors game, cancer growth in response to immunotherapy, and lipid oxidation dynamics. Our approach aids in standardizing mathematical models and in generating hypotheses based on concrete mechanistic behavior across a wide range of observed spatial phenomena.

Quantifying the roles of space and stochasticity in computer simulations for cell biology and cellular biochemistry

Most of the fascinating phenomena studied in cell biology emerge from interactions among highly organized multimolecular structures embedded into complex and frequently dynamic cellular morphologies. For the exploration of such systems, computer simulation has proved to be an invaluable tool, and many researchers in this field have developed sophisticated computational models for application to specific cell biological questions. However, it is often difficult to reconcile conflicting computational results that use different approaches to describe the same phenomenon. To address this issue systematically, we have defined a series of computational test cases ranging from very simple to moderately complex, varying key features of dimensionality, reaction type, reaction speed, crowding, and cell size. We then quantified how explicit spatial and/or stochastic implementations alter outcomes, even when all methods use the same reaction network, rates, and concentrations. For simple cases, we generally find minor differences in solutions of the same problem. However, we observe increasing discordance as the effects of localization, dimensionality reduction, and irreversible enzymatic reactions are combined. We discuss the strengths and limitations of commonly used computational approaches for exploring cell biological questions and provide a framework for decision making by researchers developing new models. As computational power and speed continue to increase at a remarkable rate, the dream of a fully comprehensive computational model of a living cell may be drawing closer to reality, but our analysis demonstrates that it will be crucial to evaluate the accuracy of such models critically and systematically.

Signaling activations through G-protein-coupled-receptor aggregations

Eukaryotic cells transmit extracellular signal information to cellular interiors through the formation of a ternary complex made up of a ligand (or agonist), G-protein, and G-protein-coupled receptor (GPCR). Previously formalized theories of ternary complex formation have mainly assumed that observable states of receptors can only take the form of monomers. Here, we propose a multiary complex model of GPCR signaling activations via the vector representation of various unobserved aggregated receptor states. Our results from model simulations imply that receptor aggregation processes can govern cooperative effects in a regime inaccessible by previous theories. In particular, we show how the affinity of ligand-receptor binding can be largely varied by various oligomer formations in the low concentration range of G-protein stimulus.

Turing and Benjamin–Feir instability mechanisms in non-autonomous systems

The Turing and Benjamin–Feir instabilities are two of the primary instability mechanisms useful for studying the transition from homogeneous states to heterogeneous spatial or spatio-temporal states in reaction–diffusion systems. We consider the case when the underlying reaction–diffusion system is non-autonomous or has a base state which varies in time, as in this case standard approaches, which rely on temporal eigenvalues, break down. We are able to establish respective criteria for the onset of each instability using comparison principles, obtaining inequalities which involve the in general time-dependent model parameters and their time derivatives. In the autonomous limit where the base state is constant in time, our results exactly recover the respective Turing and Benjamin–Feir conditions known in the literature. Our results make the Turing and Benjamin–Feir analysis amenable for a wide collection of applications, and allow one to better understand instabilities emergent due to a variety of non-autonomous mechanisms, including time-varying diffusion coefficients, time-varying reaction rates, time-dependent transitions between reaction kinetics and base states which change in time (such as heteroclinic connections between unique steady states, or limit cycles), to name a few examples.

Self-Avoiding Random Walks as a Model to Study Athermal Linear Polymers under Extreme Plate Confinement

Monte Carlo (MC) simulations, built around chain-connectivity-altering moves and a wall-displacement algorithm, allow us to simulate freely-jointed chains of tangent hard spheres of uniform size under extreme confinement. The latter is realized through the presence of two impenetrable, flat, and parallel plates. Extreme conditions correspond to the case where the distance between the plates approaches the monomer size. An analysis of the local structure, based on the characteristic crystallographic element (CCE) norm, detects crystal nucleation and growth at packing densities well below the ones observed in bulk analogs. In a second step, we map the confined polymer chains into self-avoiding random walks (SAWs) on restricted lattices. We study all realizations of the cubic crystal system: simple, body centered, and face centered cubic crystals. For a given chain size (SAW length), lattice type, origin of SAW, and level of confinement, we enumerate all possible SAWs (equivalently all chain conformations) and calculate the size distribution. Results for intermediate SAW lengths are used to predict the behavior of long, fully entangled chains through growth formulas. The SAW analysis will allow us to determine the corresponding configurational entropy, as it is the driving force for the observed phase transition and the determining factor for the thermodynamic stability of the corresponding crystal morphologies.

pSpatiocyte: a high-performance simulator for intracellular reaction-diffusion systems

Background

Studies using quantitative experimental methods have shown that intracellular spatial distribution of molecules plays a central role in many cellular systems. Spatially resolved computer simulations can integrate quantitative data from these experiments to construct physically accurate models of the systems. Although computationally expensive, microscopic resolution reaction-diffusion simulators, such as Spatiocyte can directly capture intracellular effects comprising diffusion-limited reactions and volume exclusion from crowded molecules by explicitly representing individual diffusing molecules in space. To alleviate the steep computational cost typically associated with the simulation of large or crowded intracellular compartments, we present a parallelized Spatiocyte method called pSpatiocyte.

Results

The new high-performance method employs unique parallelization schemes on hexagonal close-packed (HCP) lattice to efficiently exploit the resources of common workstations and large distributed memory parallel computers. We introduce a coordinate system for fast accesses to HCP lattice voxels, a parallelized event scheduler, a parallelized Gillespie’s direct-method for unimolecular reactions, and a parallelized event for diffusion and bimolecular reaction processes. We verified the correctness of pSpatiocyte reaction and diffusion processes by comparison to theory. To evaluate the performance of pSpatiocyte, we performed a series of parallelized diffusion runs on the RIKEN K computer. In the case of fine lattice discretization with low voxel occupancy, pSpatiocyte exhibited 74% parallel efficiency and achieved a speedup of 7686 times with 663552 cores compared to the runtime with 64 cores. In the weak scaling performance, pSpatiocyte obtained efficiencies of at least 60% with up to 663552 cores. When executing the Michaelis-Menten benchmark model on an eight-core workstation, pSpatiocyte required 45- and 55-fold shorter runtimes than Smoldyn and the parallel version of ReaDDy, respectively. As a high-performance application example, we study the dual phosphorylation-dephosphorylation cycle of the MAPK system, a typical reaction network motif in cell signaling pathways.

Conclusions

pSpatiocyte demonstrates good accuracies, fast runtimes and a significant performance advantage over well-known microscopic particle methods in large-scale simulations of intracellular reaction-diffusion systems. The source code of pSpatiocyte is available at https://spatiocyte.org.

An implicit lipid model for efficient reaction-diffusion simulations of protein binding to surfaces of arbitrary topology

Localization of proteins to a membrane is an essential step in a broad range of biological processes such as signaling, virion formation, and clathrin-mediated endocytosis. The strength and specificity of proteins binding to a membrane depend on the lipid composition. Single-particle reaction-diffusion methods offer a powerful tool for capturing lipid-specific binding to membrane surfaces by treating lipids explicitly as individual diffusible binding sites. However, modeling lipid particle populations is expensive. Here, we present an algorithm for reversible binding of proteins to continuum surfaces with implicit lipids, providing dramatic speed-ups to many body simulations. Our algorithm can be readily integrated into most reaction-diffusion software packages. We characterize changes to kinetics that emerge from explicit vs implicit lipids as well as surface adsorption models, showing excellent agreement between our method and the full explicit lipid model. Compared to models of surface adsorption, which couple together binding affinity and lipid concentration, our implicit lipid model decouples them to provide more flexibility for controlling surface binding properties and lipid inhomogeneity, thus reproducing binding kinetics and equilibria. Crucially, we demonstrate our method's application to membranes of arbitrary curvature and topology, modeled via a subdivision limit surface, again showing excellent agreement with explicit lipid simulations. Unlike adsorption models, our method retains the ability to bind lipids after proteins are localized to the surface (through, e.g., a protein-protein interaction), which can greatly increase the stability of multiprotein complexes on the surface. Our method will enable efficient cell-scale simulations involving proteins localizing to realistic membrane models, which is a critical step for predictive modeling and quantification of in vitro and in vivo dynamics.