Diffusion and retention are major determinants of protein targeting to the inner nuclear membrane.

Newly synthesized membrane proteins are constantly sorted from the endoplasmic reticulum (ER) to various membranous compartments. How proteins specifically enrich at the inner nuclear membrane (INM) is not well understood. We have established a visual in vitro assay to measure kinetics and investigate requirements of protein targeting to the INM. Using human LBR, SUN2, and LAP2ß as model substrates, we show that INM targeting is energy-dependent but distinct from import of soluble cargo. Accumulation of proteins at the INM relies on both a highly interconnected ER network, which is affected by energy depletion, and an efficient immobilization step at the INM. Nucleoporin depletions suggest that translocation through nuclear pore complexes (NPCs) is rate-limiting and restricted by the central NPC scaffold. Our experimental data combined with mathematical modeling support a diffusion-retention-based mechanism of INM targeting. We experimentally confirmed the sufficiency of diffusion and retention using an artificial reporter lacking natural sorting signals that recapitulates the energy dependence of the process in vivo.

Reaction schemes, escape times and geminate recombinations in particle-based spatial simulations of biochemical reactions.

Modeling the spatiotemporal dynamics of biochemical reaction systems at single-molecule resolution has become feasible with the increase of computing power and is applied especially to cellular signal transduction. For an association reaction the two molecules have to be in contact. Hence, a physically faithful model of the molecular interaction assumes non-overlapping molecules that interact at their surfaces (boundary scheme). For performance reasons, this model can be replaced by particles that can overlap and react when they are closer than a certain distance with a reaction probability (volume scheme). Here we present an analytical approximation for the reaction probability in the volume scheme and compare the volume- with the boundary scheme. A dissociation reaction, in contrast, creates two molecules next to each other. If the reaction is reversible, these two products can directly re-bind again, leading to an overestimation of the dimerized state in the simulation. We show how the correct recombination rate can be achieved if the products of the dissociation are placed at identical positions, but cannot react for a certain timespan. This refractory time corresponds to the completion of the diffusion-controlled dissociation of the two molecules to their contact distance r(i)+r(j) at t = τ ×(r(i)+r(j))²/(D(i)+D(j) with τ = 1/10 for molecules with radii r(i) and r(j) and diffusion coefficients D(i) and D(j), respectively.

Hybrid spatial Gillespie and particle tracking simulation.

MOTIVATION: Cellular signal transduction involves spatial-temporal dynamics and often stochastic effects due to the low particle abundance of some molecular species. Others can, however, be of high abundances. Such a system can be simulated either with the spatial Gillespie/Stochastic Simulation Algorithm (SSA) or Brownian/Smoluchowski dynamics if space and stochasticity are important. To combine the accuracy of particle-based methods with the superior performance of the SSA, we suggest a hybrid simulation. RESULTS: The proposed simulation allows an interactive or automated switching for regions or species of interest in the cell. Especially we see an application if for instance receptor clustering at the membrane is modeled in detail and the transport through the cytoplasm is included as well. The results show the increase in performance of the overall simulation, and the limits of the approach if crowding is included. Future work will include the development of a GUI to improve control of the simulation. AVAILABILITY OF IMPLEMENTATION: www.bison.ethz.ch/research/spatial_simulations. CONTACT: mklann@ee.ethz.ch or koeppl@ethz.ch Supplementary/Information: Supplementary data are available at Bioinformatics online.

Spatial simulations in systems biology: from molecules to cells.

Cells are highly organized objects containing millions of molecules. Each biomolecule has a specific shape in order to interact with others in the complex machinery. Spatial dynamics emerge in this system on length and time scales which can not yet be modeled with full atomic detail. This review gives an overview of methods which can be used to simulate the complete cell at least with molecular detail, especially Brownian dynamics simulations. Such simulations require correct implementation of the diffusion-controlled reaction scheme occurring on this level. Implementations and applications of spatial simulations are presented, and finally it is discussed how the atomic level can be included for instance in multi-scale simulation methods.

From microscopy data to in silico environments for in vivo-oriented simulations.

ABSTRACT: : In our previous study, we introduced a combination methodology of Fluorescence Correlation Spectroscopy (FCS) and Transmission Electron Microscopy (TEM), which is powerful to investigate the effect of intracellular environment to biochemical reaction processes. Now, we developed a reconstruction method of realistic simulation spaces based on our TEM images. Interactive raytracing visualization of this space allows the perception of the overall 3D structure, which is not directly accessible from 2D TEM images. Simulation results show that the diffusion in such generated structures strongly depends on image post-processing. Frayed structures corresponding to noisy images hinder the diffusion much stronger than smooth surfaces from denoised images. This means that the correct identification of noise or structure is significant to reconstruct appropriate reaction environment in silico in order to estimate realistic behaviors of reactants in vivo. Static structures lead to anomalous diffusion due to the partial confinement. In contrast, mobile crowding agents do not lead to anomalous diffusion at moderate crowding levels. By varying the mobility of these non-reactive obstacles (NRO), we estimated the relationship between NRO diffusion coefficient (Dnro) and the anomaly in the tracer diffusion (α). For Dnro=21.96 to 44.49 µm2/s, the simulation results match the anomaly obtained from FCS measurements. This range of the diffusion coefficient from simulations is compatible with the range of the diffusion coefficient of structural proteins in the cytoplasm. In addition, we investigated the relationship between the radius of NRO and anomalous diffusion coefficient of tracers by the comparison between different simulations. The radius of NRO has to be 58 nm when the polymer moves with the same diffusion speed as a reactant, which is close to the radius of functional protein complexes in a cell.

Spatial modeling of vesicle transport and the cytoskeleton: the challenge of hitting the right road.

The membrane trafficking machinery provides a transport and sorting system for many cellular proteins. We propose a mechanistic agent-based computer simulation to integrate and test the hypothesis of vesicle transport embedded into a detailed model cell. The method tracks both the number and location of the vesicles. Thus both the stochastic properties due to the low numbers and the spatial aspects are preserved. The underlying molecular interactions that control the vesicle actions are included in a multi-scale manner based on the model of Heinrich and Rapoport (2005). By adding motor proteins we can improve the recycling process of SNAREs and model cell polarization. Our model also predicts that coat molecules should have a high turnover at the compartment membranes, while the turnover of motor proteins has to be slow. The modular structure of the underlying model keeps it tractable despite the overall complexity of the vesicle system. We apply our model to receptor-mediated endocytosis and show how a polarized cytoskeleton structure leads to polarized distributions in the plasma membrane both of SNAREs and the Ste2p receptor in yeast. In addition, we can couple signal transduction and membrane trafficking steps in one simulation, which enables analyzing the effect of receptor-mediated endocytosis on signaling.

Agent-based simulation of reactions in the crowded and structured intracellular environment: Influence of mobility and location of the reactants.

BACKGROUND: In this paper we apply a novel agent-based simulation method in order to model intracellular reactions in detail. The simulations are performed within a virtual cytoskeleton enriched with further crowding elements, which allows the analysis of molecular crowding effects on intracellular diffusion and reaction rates. The cytoskeleton network leads to a reduction in the mobility of molecules. Molecules can also unspecifically bind to membranes or the cytoskeleton affecting (i) the fraction of unbound molecules in the cytosol and (ii) furthermore reducing the mobility. Binding of molecules to intracellular structures or scaffolds can in turn lead to a microcompartmentalization of the cell. Especially the formation of enzyme complexes promoting metabolic channeling, e.g. in glycolysis, depends on the co-localization of the proteins. RESULTS: While the co-localization of enzymes leads to faster reaction rates, the reduced mobility decreases the collision rate of reactants, hence reducing the reaction rate, as expected. This effect is most prominent in diffusion limited reactions. Furthermore, anomalous diffusion can occur due to molecular crowding in the cell. In the context of diffusion controlled reactions, anomalous diffusion leads to fractal reaction kinetics. The simulation framework is used to quantify and separate the effects originating from molecular crowding or the reduced mobility of the reactants. We were able to define three factors which describe the effective reaction rate, namely f diff for the diffusion effect, f volume for the crowding, and f access for the reduced accessibility of the molecules. CONCLUSIONS: Molecule distributions, reaction rate constants and structural parameters can be adjusted separately in the simulation allowing a comprehensive study of individual effects in the context of a realistic cell environment. As such, the present simulation can help to bridge the gap between in vivo and in vitro kinetics.

Multi-scale spatio-temporal modeling: lifelines of microorganisms in bioreactors and tracking molecules in cells.

Agent-based models are rigorous tools for simulating the interactions of individual entities, such as organisms or molecules within cells and assessing their effects on the dynamic behavior of the system as a whole. In context with bioprocess and biosystems engineering there are several interesting and important applications. This contribution aims at introducing this strategy with the aid of two examples characterized by striking distinctions in the scale of the individual entities and the mode of their interactions. In the first example a structured-segregated model is applied to travel along the lifelines of single cells in the environment of a three-dimensional turbulent field of a stirred bioreactor. The modeling approach is based on an Euler-Lagrange formulation of the system. The strategy permits one to account for the heterogeneity present in real reactors in both the fluid and cellular phases, respectively. The individual response of the cells to local variations in the extracellular concentrations is pictured by a dynamically structured model of the key reactions of the central metabolism. The approach permits analysis of the lifelines of individual cells in space and time.The second application of the individual modeling approach deals with dynamic modeling of signal transduction pathways in individual cells. Usually signal transduction networks are portrayed as being wired together in a spatially defined manner. Living circuitry, however, is placed in highly malleable internal architecture. Creating a homogenous bag of molecules, a well-mixed system, the dynamic behavior of which is modeled with a set of ordinary differential equations is normally not valid. The dynamics of the MAP kinase and a steroid hormone pathway serve as examples to illustrate how single molecule tracking can be linked with the stochasticity of biochemical reactions, where diffusion and reaction occur in a probabilistic manner. The problem of hindered diffusion caused by macromolecular crowding is also taken into account.

Stochastic simulation of signal transduction: impact of the cellular architecture on diffusion.

The transduction of signals depends on the translocation of signaling molecules to specific targets. Undirected diffusion processes play a key role in the bridging of spaces between different cellular compartments. The diffusion of the molecules is, in turn, governed by the intracellular architecture. Molecular crowding and the cytoskeleton decrease macroscopic diffusion. This article shows the use of a stochastic simulation method to study the effects of the cytoskeleton structure on the mobility of macromolecules. Brownian dynamics and single particle tracking were used to simulate the diffusion process of individual molecules through a model cytoskeleton. The resulting average effective diffusion is in line with data obtained in the in vitro and in vivo experiments. It shows that the cytoskeleton structure strongly influences the diffusion of macromolecules. The simulation method used also allows the inclusion of reactions in order to model complete signaling pathways in their spatio-temporal dynamics, taking into account the effects of the cellular architecture.