BIOPHYSICAL THEORY AND MODELING |
Coarse-grained molecular simulation of diffusion and reaction kinetics in a crowded virtual cytoplasm
Douglas Ridgway 1, Gordon Broderick 1, Ana Lopez-Campistrous 1, Melania Ru'aini 1, Philip Winter 1, Matthew Hamilton 1, Pierre Boulanger 1, Andriy Kovalenko 2 and Michael J Ellison 1*
1 University of Alberta
2 National Institute for Nanotechnology National Research Council of Canada
* To whom correspondence should be addressed. E-mail: mike.ellison{at}ualberta.ca.
Submitted on June 22, 2007
Revised on August 29, 2007
Accepted on 27 November 2007
 |
Abstract |
|---|
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 (time-step, move length, 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 catalogue from E. coli K12. Simulated diffusion constants exhibited anomalous behaviour 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 time scales. These observations illustrate that quantitative realism in biosimulation will depend to some extent on meso-scale phenomena that is not currently well understood.
Key Words:
anomalous diffusion, computer modeling, enzyme kinetics, intracellular environment, proteomics, systems biology
