Effects of quenched and annealed macromolecular crowding elements on a simple model for signaling in T lymphocytes

Three-dimensional stochastic off-lattice model of binding chemistry in crowded environments

Molecular crowding is one of the characteristic features of the intracellular environment, defined by a dense mixture of varying kinds of proteins and other molecules. Interaction with these molecules significantly alters the rates and equilibria of chemical reactions in the crowded environment. Numerous fundamental activities of a living cell are strongly influenced by the crowding effect, such as protein folding, protein assembly and disassembly, enzyme activity, and signal transduction. Quantitatively predicting how crowding will affect any particular process is, however, a very challenging problem because many physical and chemical parameters act synergistically in ways that defy easy analysis. To build a more realistic model for this problem, we extend a prior stochastic off-lattice model from two-dimensional (2D) to three-dimensional (3D) space and examine how the 3D results compare to those found in 2D. We show that both models exhibit qualitatively similar crowding effects and similar parameter dependence, particularly with respect to a set of parameters previously shown to act linearly on total reaction equilibrium. There are quantitative differences between 2D and 3D models, although with a generally gradual nonlinear interpolation as a system is extended from 2D to 3D. However, the additional freedom of movement allowed to particles as thickness of the simulation box increases can produce significant quantitative change as a system moves from 2D to 3D. Simulation results over broader parameter ranges further show that the impact of molecular crowding is highly dependent on the specific reaction system examined.

Unified regression model of binding equilibria in crowded environments

Molecular crowding is a critical feature distinguishing intracellular environments from idealized solution-based environments and is essential to understanding numerous biochemical reactions, from protein folding to signal transduction. Many biochemical reactions are dramatically altered by crowding, yet it is extremely difficult to predict how crowding will quantitatively affect any particular reaction systems. We previously developed a novel stochastic off-lattice model to efficiently simulate binding reactions across wide parameter ranges in various crowded conditions. We now show that a polynomial regression model can incorporate several interrelated parameters influencing chemistry under crowded conditions. The unified model of binding equilibria accurately reproduces the results of particle simulations over a broad range of variation of six physical parameters that collectively yield a complicated, non-linear crowding effect. The work represents an important step toward the long-term goal of computationally tractable predictive models of reaction chemistry in the cellular environment.

Coarse-grained molecular simulation of diffusion and reaction kinetics in a crowded virtual cytoplasm

We present a general-purpose model for biomolecular simulations at the molecular level that incorporates stochasticity, spatial dependence, and volume exclusion, using diffusing and reacting particles with physical dimensions. To validate the model, we first established the formal relationship between the microscopic model parameters (timestep, move length, and reaction probabilities) and the macroscopic coefficients for diffusion and reaction rate. We then compared simulation results with Smoluchowski theory for diffusion-limited irreversible reactions and the best available approximation for diffusion-influenced reversible reactions. To simulate the volumetric effects of a crowded intracellular environment, we created a virtual cytoplasm composed of a heterogeneous population of particles diffusing at rates appropriate to their size. The particle-size distribution was estimated from the relative abundance, mass, and stoichiometries of protein complexes using an experimentally derived proteome catalog from Escherichia coli K12. Simulated diffusion constants exhibited anomalous behavior as a function of time and crowding. Although significant, the volumetric impact of crowding on diffusion cannot fully account for retarded protein mobility in vivo, suggesting that other biophysical factors are at play. The simulated effect of crowding on barnase-barstar dimerization, an experimentally characterized example of a bimolecular association reaction, reveals a biphasic time course, indicating that crowding exerts different effects over different timescales. These observations illustrate that quantitative realism in biosimulation will depend to some extent on mesoscale phenomena that are not currently well understood.