Investigating the robustness of the classical enzyme kinetic equations in small intracellular compartments
- Ramon Grima^{1}Email author
https://doi.org/10.1186/1752-0509-3-101
© Grima; licensee BioMed Central Ltd. 2009
Received: 2 June 2009
Accepted: 8 October 2009
Published: 8 October 2009
Abstract
Background
Classical descriptions of enzyme kinetics ignore the physical nature of the intracellular environment. Main implicit assumptions behind such approaches are that reactions occur in compartment volumes which are large enough so that molecular discreteness can be ignored and that molecular transport occurs via diffusion. Though these conditions are frequently met in laboratory conditions, they are not characteristic of the intracellular environment, which is compartmentalized at the micron and submicron scales and in which active means of transport play a significant role.
Results
Starting from a master equation description of enzyme reaction kinetics and assuming metabolic steady-state conditions, we derive novel mesoscopic rate equations which take into account (i) the intrinsic molecular noise due to the low copy number of molecules in intracellular compartments (ii) the physical nature of the substrate transport process, i.e. diffusion or vesicle-mediated transport. These equations replace the conventional macroscopic and deterministic equations in the context of intracellular kinetics. The latter are recovered in the limit of infinite compartment volumes. We find that deviations from the predictions of classical kinetics are pronounced (hundreds of percent in the estimate for the reaction velocity) for enzyme reactions occurring in compartments which are smaller than approximately 200 nm, for the case of substrate transport to the compartment being mediated principally by vesicle or granule transport and in the presence of competitive enzyme inhibitors.
Conclusion
The derived mesoscopic rate equations describe subcellular enzyme reaction kinetics, taking into account, for the first time, the simultaneous influence of both intrinsic noise and the mode of transport. They clearly show the range of applicability of the conventional deterministic equation models, namely intracellular conditions compatible with diffusive transport and simple enzyme mechanisms in several hundred nanometre-sized compartments. An active transport mechanism coupled with large intrinsic noise in enzyme concentrations is shown to lead to huge deviations from the predictions of deterministic models. This has implications for the common approach of modeling large intracellular reaction networks using ordinary differential equations and also for the calculation of the effective dosage of competitive inhibitor drugs.
Background
The inside of a cell is a highly complex environment. In the past two decades, detailed measurements of the chemical and biophysical properties of the cytoplasm have established that the conditions in which intracellular reactions occur are, by and large, very different than those typically maintained in laboratory conditions. One of the outstanding differences between in vivo and in vitro conditions, is that in the former, biochemical reactions typically occur in minuscule reaction volumes [1]. For example, in eukaryotic cells, many biochemical pathways are sequestered within membrane-bound compartments, ranging from ~50 nm diameter vesicles to the nucleus, which can be several microns in size [2]. It is also found that the total concentration of macromolecules inside both prokaryotic and eukaryotic cells is very large [3, 4], of the order of 50 - 400 mg/ml which implies that between 5% and 40% of the total intracellular volume is physically occupied by these molecules [5]. The concentration of these crowding molecules is highly heterogeneous (see for example [6]), meaning that typically one will find small pockets of intracellular space, characterized by low macromolecular crowding, surrounded by a "sea" of high crowding; such pockets of space may serve as effective compartments where reactions may occur more easily than in the rest of the cytosol. Analysis of experimental data for the dependence of diffusion coefficients with molecular size suggests the length scale of such effective compartments is in the range 35-50 nm [7], a size comparable to that of the smallest vesicles. The significant crowding also suggests that frequently an active means of transport such as vesicle-mediated transport, may be more desirable than simple diffusion as a means of intracellular transport.
The volume of a spherical cavity of space of diameter 50 nm is merely ~6.5 × 10^{-20} liters, an extremely small number compared to the typical macroscopic reaction volumes of in vitro experiments (experimental attolitre biochemistry is still in its infancy - see for example [8]). These very small reaction volumes imply that at physiologically relevant concentrations (nano to millimolar), the copy number of a significant number of intracellular molecules is very small [1] and consequently that intrinsic noise cannot be ignored; for example 255 μM corresponds to an average of just 10 molecules in a 50 nm vesicle and fluctuations about this mean of the order of 3 molecules [9].
The traditional mathematical framework of physical chemistry ignores the basic physical properties of the intracellular environment. Kinetics are described by a set of coupled ordinary differential equations which implicitly assume (i) that the reaction compartment is so large that molecular discreteness can be ignored and that hence integer numbers of molecules per unit volume can be replaced by a continuous variable, the molar concentration. Since the number of molecules is assumed to be very large, stochastic fluctuations are deemed negligible and the equations are hence deterministic; (ii) the reaction compartment is well-stirred so that homogeneous conditions prevail throughout [9]. Both assumptions can be justified for reactions occurring in a constantly stirred reactor of macroscopic dimensions. However if diffusion is the dominant transport process inside the compartment then the homogeneity assumption holds only if the volume is small enough so that in the time between successive reactions, a molecule will diffuse a distance much larger than the size of the compartment. This comes at the expense of the first assumption. It hence appears natural that for intracellular applications, the first assumption, namely that of deterministic kinetics cannot be justified a priori. The second assumption can be justified if reactions are localized in sufficiently small parts of the cell and in particular for reaction-limited processes i.e. those for which the typical time for two molecules to meet each other via diffusion is much less than the typical time for them to react if they are in close proximity. For such conditions, a molecule will come within reaction range several times before participating in a successful reaction, in the process sampling the compartment many times which naturally leads to well-mixed conditions [9–11].
Results
Model I: Michaelis-Menten reaction occurring in a compartment volume of sub-micron dimensions. Substrate input into compartment occurs via a Poisson process
Hence the prediction of the MM equation is only correct, i.e. α = α_{ M }, in the limit of infinitely large compartment volumes, in which case the second term on the left hand side of Eq. (7) will become vanishingly small and can be neglected. For finite compartment volumes, the MM equation is not exact (except in the two limiting cases of α_{ M }→ 0 and α_{ M }→ 1) but is at best an approximation, even though steady-state conditions are imposed; this is at odds with the prediction of the conventional deterministic theory. An inspection of Eqs. (7) and (8) shows that the magnitude of the deviations from the MM equation depends on the two non-dimensional quantities: (i) K_{ M }Ω, a measure of the rate at which enzyme-substrate combination events occur relative to the rate of decay of complex molecules; (ii) [E_{ T }]Ω, the total integer number of enzyme molecules in the compartment.
As shown in the Methods section, the MM equation is found to implicity assume that the noise about the macroscopic substrate and enzyme concentrations is uncorrelated (this assumption has generally been found to be at the heart of many macroscopic models - for example see [16]); properly taking into account these non-zero correlations leads to the corrections encapsulated by Eqs. (7) and (8). These correlations are expected to be small in two particular cases: (i) if K_{ M }is large; in this case when substrate molecules combine with an enzyme to form a complex, the latter dissociates very quickly back into free enzyme and thus successive enzyme-substrate events to the same enzyme molecule are bound to be almost independent of each other. The opposite situation of small K_{ M }would imply that the bottleneck in the catalytic process is the decay of complex rather than enzyme-substrate combination; if a successful combination occurs, the next substrate to arrive to the same enzyme molecule would have to wait until the complex decays, naturally leading to correlations between successive enzyme-substrate combination events. (ii) if the total number of enzyme molecules is large; in such a case, at any one time, the noise about the macroscopic concentrations will be the sum total from a large number of enzymes, each at a different stage in the catalytic process and each independent from all others, which naturally dilutes any temporal correlations.
Maximum Percentage error in reaction velocity from prediction of the MM equation for Model I.
D/nm | K_{ M }= 10 μM | 100 μM | 1000 μM | Copy No. |
---|---|---|---|---|
50 | 11.83 [17.00] | 4.09 [4.33] | 0.59 | 10 |
100 | 4.74 [5.00] | 0.73 [0.74] | 0.08 | 10 |
200 | 0.90 | 0.10 | 0.01 | 10 |
50 | 3.98 [5.33] | 1.88 [2.02] | 0.43 | 100 |
100 | 2.10 [2.23] | 0.52 [0.52] | 0.07 | 100 |
200 | 0.61 | 0.09 | 0.01 | 100 |
The theory is always found to underestimate the actual deviations predicted by simulations; hence the theoretical expressions provide a quick, convenient way by which one can generally estimate a lower bound on the deviations to be expected from the MM equation without the need to perform extensive stochastic simulation.
Model II: Michaelis-Menten reaction occurring in a compartment volume of sub-micron dimensions. Substrate is input into compartment in groups or bursts of M molecules at a time
Model I captures the basics of a general enzyme-catalyzed process occurring in a small intracellular compartment. In this section we build upon this model to incorporate further biological realism. In particular, in the previous model we assumed that substrate input can be well described by a Poisson process, where one molecule at a time is fed into the compartment with some average rate k_{ in }. This is the simplest possible assumption and approximates well the situation in which molecules are brought to the compartment via normal diffusion. However there are many situations where this may not be the case; we now describe two such cases.
The intracellular condition of macromolecular crowding limits the Brownian motion of molecules in the cytoplasm, this being reflected in the relatively small diffusion coefficients measured in vivo compared to those known in vitro for moderately to relatively large molecules. Experiments with inert tracer particles in the cytoplasm of Swiss 3T3 cells show that the in vivo diffusion coefficient is an order of magnitude less than that in vitro for molecules with hydrodynamic radius 14 nm and diffusion becomes negligibly small for molecules larger than approximately 25 nm [7]; similar results have been obtained in Xenopus neurons [18] and skeletal muscle myotubes [19]. If diffusion is considerably hindered, one expects active transport to become a more desirable mode of transport. Indeed there exists ample evidence for the active transport of macromolecules: they are typically packaged in a vesicle or a granule which is then transported along microtubules or by some other means. It is also found that each vesicle or granule typically contains several of these molecules (examples are: mRNA molecules - several estimated per granule [20, 21]; cholesterol molecules which are transported in low-density lipoproteins [2] - approximately 1500 per lipoprotein).
Generally an active means of transport is not exclusively linked with the transport of large substrate molecules. The cell being a highly compartmentalized and dynamic entity requires for its survival the precise transport of certain molecules from one compartment to another and a regulation of this transport depending on its current physiological state. Brownian motion leads to an isotropic movement of molecules down the concentration gradient and to a consequent damping of the substrate concentration with distance. In contrast active transport provides a directed (anisotropic) means of transport with little or no loss of substrate with distance, is independent of the concentration gradient and it is also easily amenable to modulation.
Maximum Percentage error in reaction velocity from prediction of the MM equation for Model II.
D/nm | K_{ M }= 10 μM | 100 μM | 1000 μM | Copy No. |
---|---|---|---|---|
50 | 225.40 | 152.83 [291.56] | 45.43 | 10 |
100 | 161.59 [331.66] | 52.74 [58.39] | 6.82 [6.99] | 10 |
200 | 65.09 | 8.45 [8.50] | 0.88 | 10 |
50 | 32.97 | 30.17 [61.66] | 18.14 | 100 |
100 | 30.78 [66.03] | 19.76 [24.52] | 5.57 [6.06] | 100 |
200 | 21.27 | 6.61 [6.91] | 0.85 | 100 |
Model III: Michaelis-Menten reaction with competitive inhibitor occurring in a compartment volume of sub-micron dimensions. Substrate input as in previous models
Maximum Percentage error in reaction velocity from prediction of the MM equation for Model III.
D/nm | K_{ M }= 10 μM | 100 μM | 1000 μM | M (burst size) |
---|---|---|---|---|
50 | 67.8 [76.8] | 67.8 [76.5] | 67.8 | 1 |
100 | 20.8 [26.4] | 20.6 [26.1] | 20.6 | 1 |
200 | 2.8 | 2.7 | 2.7 | 1 |
50 | 1001.8 | 234.9 [169.4] | 86 | 50 |
100 | 343.7 [345.5] | 73.4 [75.2] | 26.2 [31.5] | 50 |
200 | 71.4 | 11.3 [11.5] | 3.6 | 50 |
Discussion and Conclusion
In this last section we discuss some fine points regarding: (i)the assumptions behind the use of master equations which throws light on the range of use of the derived mesoscopic equations, (ii) the use of the system-size expansion to perturbatively solve the master equation and (iii) the assumption of steady-state metabolic conditions. We conclude by placing our work in the context of previous recent studies of stochastic enzyme kinetics and discuss possible experiments to verify some of the conclusions we have reached.
We have implicitly assumed throughout the article that a single (global) master equation model suffices to capture the deviations from classical kinetics due to fluctuations in chemical concentrations inside a single subcellular compartment. As noted by Baras and Mansour [23], "the global master equation selects the very limited class of exceptionally large fluctuations that appear at the level of the entire system, disregarding important nonequilibrium features originated by local fluctuations." Hence the results presented here necessarily underestimate the possible deviations from classical kinetics, in particular the local fluctuations due to diffusion of molecules inside the compartment. These local fluctuations are typically small for reaction-limited processes (as in this article) but significant for diffusion-limited ones. To capture them effectively, one would be required to spatially discretize the compartment into many small elements and describe the reaction-diffusion processes between these elements by means of a multivariate master equation [12, 23]. The latter is known as a reaction-diffusion master equation; typically it does not allow detailed analytical investigation as for a global master equation and one is limited to stochastic simulation. Use of the global master equation is also restricted for compartments which are not too small: in particular the linear dimensions of the compartment should be larger than the average distance traveled by a molecule before undergoing a successful reaction with another molecule i.e. the length scale is much larger than that inherent in molecular dynamics simulation [23].
We have applied the systematic expansion due to van Kampen to perturbatively solve the master equation. It is sometimes a priori assumed that because this expansion is about the macroscopic concentrations, it cannot give information regarding the stochastic kinetics of few particle/small volume systems. This is true if one restricts oneself to the expansion to order Ω^{0} i.e. the linear-noise approximation; this is commonly the case found in the literature since the algebra becomes tedious if one considers more terms. However we note that as argued and shown by van Kampen himself [12], terms beyond the linear-noise approximation in the system-size expansion add terms to the fluctuations that are of order of a single particle relative to the macroscopic quantities and are essential to understanding how fluctuations are affected by the presence of non-linear terms in the macroscopic equation (substrate-enzyme binding in our case). In our theory we went beyond the linear-noise approximation. We find that the predicted theoretical results are in reasonable agreement, in many cases (comparison of bold and italic values in Tables 1, 2 and 3), with stochastic simulations of just a few tens of enzyme molecules in sub-micron compartments, which justifies our methodology.
We have also imposed metabolic steady-state conditions inside the subcellular compartment. Technically this is convenient since in such a case one does not deal with complex transients. Also since under such conditions the MM equation is exact from a deterministic point of view, it provides a very useful reference point versus which to accurately compute deviations due to intrinsic noise. In reality one may not always have steady-state conditions inside cells, this depending strongly on the rate of substrate input relative to the maximum rate at which the enzyme can process substrate. Another possibility is that one is dealing with a batch reaction i.e. one in which a number of substrate molecules are transported at one go and just once to the subcellular compartment (e.g. via vesicle-mediated transport) and the reaction proceeds thereafter without any further substrate replenishment. This latter scenario is compatible with the presentation of the MM equation typical in standard physical chemistry textbooks. The MM equation is then an approximation (not exact as in steady-state case) to the deterministic kinetics, when substrate is present in much larger concentration than enzyme. This case is currently under investigation using the same perturbative framework used in this article.
We note that this is not the first attempt to study stochastic enzyme kinetics. The bulk of recent studies [24–27] have focused on understanding the kinetics of a Michaelis-Menten type reaction catalyzed by a single enzyme molecule. Deviations from classical kinetics were found to be most pronounced when one takes into account substrate fluctuations [26]. These pioneering studies were restricted to a single-enzyme assisted reaction which reduces complexity thereby making it ideal from a theoretical perspective; since the reaction is dependent on just a single enzyme molecule one also finds maximum deviations from deterministic kinetics. In reality, it is unlikely to find just one enzyme molecule inside a subcellular compartment - as mentioned in the introduction a physiological concentration of just a few hundred micromolar would correspond to few tens inside the typically smallest subcellular compartment. It is also the case that diffusion may not always be the main means of substrate transport to the compartment and that the reaction maybe more complex than the simple Michaelis-Menten type reaction of these previous studies. The present study fills in these gaps by using a systematic method to derive approximate and relatively simple analytic expressions for mesoscopic rate equations describing the kinetics of the general case of N enzyme molecules in a subcellular compartment with or without active transport of substrate and in the presence of enzyme inhibitors. Most importantly our approach shows the effects of intrinsic noise on the kinetics can be captured via effective ordinary differential equations. This enables quick estimation of the magnitude of stochastic effects on reaction kinetics and thus gives insight into whether a model or parts of a model should be designed to be stochastic or deterministic without the need for extensive stochastic simulation. In the present study, this approach enabled us to readily compute, for the first time, the deviations from deterministic kinetics for a broad range of realistic in vivo parameter constants (Tables 1, 2 and 3), a task which would be considerably lengthy if one had to rely solely on data obtained from ensemble-averaged stochastic simulations.
We conclude by briefly discussing possible experiments which can verify the predictions made in this article. It is arguably not an easy task to perform the required experiments in real-time in a living cell. A viable alternative would consist of monitoring reaction kinetics inside single artificially-made vesicles. Pick et al [8] have shown that the addition of cytochalasin to mammalian cells induces them to extrude from their plasma membrane minuscule vesicles of attolitre volume with fully functional cell surface receptors and also retaining cytosolic proteins in their interior. The change in the intra-vesicular calcium ion concentration in response to surface ligand binding was measured using fluorescence confocal microscopy (FCM). Since the vesicle sizes are of typical small sub-cellular compartment dimensions (1 attolitre corresponds to a spherical vesicle of approximate diameter 120 nm) and FCM allows the measurement of the concentration of a fluorescent probe (via a calibration procedure), this experimental technique appears ideal to verify the predictions of Model I and of Model III for the case of diffusive substrate transport. Model II and Model III with vesicle-transport of substrate are probably much more challenging to verify since one then needs to construct the in vitro equivalent of microtubules. This is within the scope of synthetic biology and may be a possibility in the next few years.
Methods
We here provide full details of the calculations reported in the Results section. The system size-expansion which is at the heart of the analysis has to-date not been applied extensively to biological problems and thus we go into some detail in its elucidation in Sub section I, which is dedicated exclusively to Model I. For other recent applications of the general method in the context of reaction kinetics, see for example [28] and [29]. Subsections II and III (treating Model II and Model III, respectively) naturally build on the results of the first subsection and thus we only give the main steps of the calculations in these last two cases. Sub section IV has a brief discussion of the simulation methods used to verify the theoretical results.
Model I: Michaelis-Menten reaction occurring in a compartment volume of sub-micron dimensions. Substrate input into compartment is modeled as a Poisson process
where π = π(n_{ C }, n_{ P }, n_{ S }), n_{ X }is the integer number of molecules of type X (where X = C, P, S), Ω is the compartment volume, and are the step operators defined by their action on a general function g(n_{ X }) as: g(n_{ X }) = g(n_{ X }± 1). Note that the relevant variables are three, not four: the integer number of molecules of free enzyme (n_{ E }) is not an independent variable due to the fact that the total amount of enzyme is conserved. The master equation cannot be solved exactly but it is possible to systematically approximate it by using an expansion in powers of the inverse square root of the volume of the compartments. This is generally called the system-size expansion [12].
Note that in Eq. (18) terms which involve products of first and second-order derivatives, third-order derivatives or higher have been omitted - these do not affect the low-order moment equations which we will be calculating.
Analysis of Ω^{1/2} terms
where v_{ max }= k_{2} [E_{ T }] is the maximum reaction velocity, [E_{ T }] = [E] + [C] is the total enzyme concentration which is a constant at all times and K_{ M }= (k_{1} + k_{2})/k_{0} is the Michaelis-Menten constant.
Analysis of Ω^{0} terms
where α = k_{ in }/v_{ max }is the normalized reaction velocity of the enzyme.
Analysis of Ω^{-1/2} terms
An alternative mesoscopic rate equation replacing the MM equation
We finish this section by noting that Eq. (30) will be found to be valid generally and not only for the simple Michaelis-Menten scheme treated in this section; the details of the reaction network come in through the form of Eq. (31) which is reaction-specific.
Model II: Michaelis-Menten reaction occurring in a compartment volume of sub-micron dimensions. Substrate is input into compartment in groups or bursts of M molecules at a time
This clearly shows that generally larger deviations from the predictions of the MM equation are expected in this case compared to those computed for Model I.
Model III: Michaelis-Menten reaction with competitive inhibitor occurring in a compartment volume of sub-micron dimensions. Substrate input as in two previous models
The change of variables from n_{ X }to ϵ_{ X }is done as before, however note that now the conservation law for enzyme is different than in the two previous models. The total enzyme concentration is now equal to [E_{ T }] = [E] + [C] + [V] from which it follows that n_{ E }= Ω[E] - Ω^{1/2}(ϵ_{ C }+ ϵ_{ V }). The description is chosen to be in terms of numbers of molecules of types C, S and V and thus E being a dependent variable does not show up explicitly in the step operators of the master equation above.
Analysis of Ω^{1/2} terms
where β = [I]/K_{ i }and K_{ i }= k_{3}/ is the dissociation constant of the inhibitor.
Analysis of Ω^{0} and Ω^{-1/2} terms
In the above equations we have defined k' = k_{3} + k_{4}. Note also that the system of equations has been simplified through the application of a few row operations.
Note that = 0 such that at α = 0, there is no correction to the MM equation i.e. α_{ M }= 0 also. The case β = 0 reduces to Model II, i.e. f(α) is given by Eq. (34).
Stochastic simulation
In this section we briefly describe the simulation methods used to verify the theoretical results which are described in detail in the Results section. All simulations were carried out using Gillespie's exact stochastic simulation algorithm, conveniently implemented in the standard simulation platform, Dizzy [32].
The data points in Figure 2 were generated by iterating the following four-step procedure: (i) pick a value for α between 0 and 1. This gives the substrate input rate k_{ in }= α v_{ max }; (ii) run the simulation and measure the ensemble-averaged substrate concentration, ⟨n_{ S }/Ω⟩ = [S*] at long times; (iii) compute α_{ M }using the MM equation, α_{ M }= [S*]/([S*]+ K_{ M }); (iv) compute the absolute percentage error R_{ p }= 100|(1 - α_{ M }/α)|. The solid curves in Figure 2 were obtained by numerically solving the cubic polynomial in α given by Eqs. (7) and (8) in the Results section for given values of α_{ M }and then using the above expression for R_{ p }. Figure 3 is generated in the same manner as Figure 2, except that: in step (i) we fix M and pick a value for α between 0 and 1. Since k_{ in }= M , the required simulation parameter is = α v_{ max }/M; step (iv) is not needed. The solid curves were obtained by numerically solving the cubic polynomial in α given by Eqs. (7) and (9) in the Results section for given values of [S*]. The y-axis for this figure is v/v_{ max }= α_{ M }for the MM equation and v/v_{ max }= α for the stochastic model. Figure 4 is obtained by numerically solving the quintic polynomial in α given by Eqs. (7) and (12) in the Results section together with the coefficients given by Eqs. (52)-(62) in the present section; the inhibitor concentration, [I], is varied while the substrate concentration, [S*], is kept fixed. The substrate concentration is chosen so that at [I] = 0, v/v_{ max }= 0.909 in all cases. Note that for models I and II, α_{ M }= [S*]/([S*] + K_{ M }) while for Model III, α_{ M }= [S*]/([S*] + (1 + β)K_{ M }). Note that the error bars are very small on the scale of the figures and thus are not shown.
Declarations
Acknowledgements
It is a pleasure to thank Arthur Straube and Philipp Thomas for interesting discussions. The author gratefully acknowledges support from SULSA (Scottish Universities Life Sciences Alliance).
Authors’ Affiliations
References
- Luby-Phelps K: Cytoarchitecture and Physical properties of Cytoplasm: Volume, Viscosity, Diffusion, Intracellular Surface Area. Int Rev Cytology. 1999, 192: 189-221. full_text. full_textView ArticleGoogle Scholar
- Alberts B, et al.: Molecular Biology of the Cell. 1994, Garland PublishingGoogle Scholar
- Minton AP: How can biochemical reactions within cells differ from those in test tubes?. J Cell Sci. 2006, 119: 2863-2869. 10.1242/jcs.03063View ArticlePubMedGoogle Scholar
- Trepat X, et al.: Universal physical responses to stretch in the living cell. Nature. 2007, 447: 592-596. 10.1038/nature05824PubMed CentralView ArticlePubMedGoogle Scholar
- Schnell S, Turner TE: Reaction kinetics in intracellular environments with macro-molecular crowding: simulations and rate laws. Prog Biophys Mol Biol. 2004, 85: 235-260. 10.1016/j.pbiomolbio.2004.01.012View ArticlePubMedGoogle Scholar
- Medalia O, et al.: Macromolecular Architecture in Eukaryotic Cells Visualized by Cryoelectron tomography. Science. 2004, 298: 1209-1213. 10.1126/science.1076184.View ArticleGoogle Scholar
- Luby-Phelps K, Castle PE, Taylor DL, Lanni F: Hindered diffusion of inert tracer particles in the cytoplam of mouse 3T3 cells. Proc Natl Acad Sci USA. 1987, 84: 4910-4913. 10.1073/pnas.84.14.4910PubMed CentralView ArticlePubMedGoogle Scholar
- Pick H, et al.: Investigating Cellular Signaling Reactions in Single Attoliter vesicles. J Am Chem Soc. 2005, 127: 2908-2912. 10.1021/ja044605xView ArticlePubMedGoogle Scholar
- Grima R, Schnell S: Modelling reaction kinetics inside cells. Essays in Biochemistry. 2008, 45: 41-56. 10.1042/BSE0450041PubMed CentralView ArticlePubMedGoogle Scholar
- Gillespie DT: Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977, 81: 2340-2361. 10.1021/j100540a008.View ArticleGoogle Scholar
- Ben-Avraham D, Havlin S: Diffusion and Reactions in Fractals and Disordered Systems. 2000, Cambridge University PressView ArticleGoogle Scholar
- van Kampen NG: Stochastic processes in physics and chemistry. 2007, Amsterdam: ElsevierGoogle Scholar
- Bartholomay AF: A Stochastic Approach to Statistical Kinetics with Application to Enzyme Kinetics. Biochemistry. 1962, 1: 223-230. 10.1021/bi00908a005View ArticlePubMedGoogle Scholar
- Bartholomay AF: Enzymatic Reaction-Rate Theory - A Stochastic Approach. Annals of the New York Academy of Sciences. 1962, 96: 897- 10.1111/j.1749-6632.1962.tb54110.xView ArticlePubMedGoogle Scholar
- Jachimowski CJ, McQuarrie DA, Russell ME: A Stochastic Approach to Enzyme-Substrate Reactions. Biochemistry. 1964, 3: 1732-1736. 10.1021/bi00899a025View ArticlePubMedGoogle Scholar
- Grima R: Multiscale modeling of biological pattern formation. Curr Top Dev Biol. 2008, 81: 435-460. full_textView ArticlePubMedGoogle Scholar
- Berg JM, Tymoczko JL, Stryer L: Biochemistry. 2002, New York: Freeman, 5Google Scholar
- Popov S, Poo MM: Diffusional transport of macromolecules in developing nerve processes. J Neurosci. 1992, 12: 77-85.PubMedGoogle Scholar
- Arrio-Dupont M, Cribier S, Foucault G, Devaux PF, d'Albis A: Diffusion of fluorescently labelled macromolecules in cultured muscle cells. Biophys J. 1996, 70: 2327-2332. 10.1016/S0006-3495(96)79798-9PubMed CentralView ArticlePubMedGoogle Scholar
- Ainger K, et al.: Transport and localization of exogenous myelin basic protein mRNA microinjected into Oligodendrocytes. J Cell Biol. 1993, 123: 431-441. 10.1083/jcb.123.2.431View ArticlePubMedGoogle Scholar
- Bassell G, Singer RH: mRNA and cytoskeletal filaments. Curr Opin Cell Biol. 1997, 9: 109-115. 10.1016/S0955-0674(97)80159-7View ArticlePubMedGoogle Scholar
- Fersht A: Structure and mechanism in protein science. 1998, New York: FreemanGoogle Scholar
- Baras F, Mansour MM: Reaction-diffusion master equation: A comparison with microscopic simulations. Phys Rev E. 1996, 54: 6139-6148. 10.1103/PhysRevE.54.6139.View ArticleGoogle Scholar
- English BP, Min W, van Oijen AM, Lee KT, Luo G, et al.: Ever-fluctuating single enzyme molecules: Michaelis-Menten equation revisited. Nat Chem Biol. 2005, 2: 87-94. 10.1038/nchembio759View ArticlePubMedGoogle Scholar
- Kou SC, Cherayil BJ, Min W, English BP, Xie SX: Single-molecule Michaelis-Menten equations. J Phys Chem B. 2005, 109: 19068-19081. 10.1021/jp051490qView ArticlePubMedGoogle Scholar
- Stefanini MO, McKane AJ, Newman TJ: Single enzyme pathways and substrate fluctuations. Nonlinearity. 2005, 18: 1575-1595. 10.1088/0951-7715/18/4/008.View ArticleGoogle Scholar
- Qian H, Elson EL: Single-molecule enzymology: stochastic Michaelis-Menten kinetics. Biophys Chem. 2002, 101-102: 565-576. 10.1016/S0301-4622(02)00145-XView ArticlePubMedGoogle Scholar
- Grima R: Noise-induced breakdown of the Michaelis-Menten equation in steady-state conditions. Phys Rev Letts. 2009, 102: 218103-10.1103/PhysRevLett.102.218103.View ArticleGoogle Scholar
- McKane AJ, Nagy J, Newman TJ, Stefanini M: Amplified biochemical oscillations in cellular systems. J Stat Phys. 2007, 128: 165-191. 10.1007/s10955-006-9221-9.View ArticleGoogle Scholar
- Ngo LG, Roussel MR: A new class of biochemical oscillators based on competitive binding. Eur J Biochem. 1997, 245: 182-190. 10.1111/j.1432-1033.1997.00182.xView ArticlePubMedGoogle Scholar
- Davis KL, Roussel MR: Optimal observability of sustained stochastic competitive inhibition oscillations at organellar volumes. FEBS journal. 2006, 273: 84-95. 10.1111/j.1742-4658.2005.05043.xView ArticlePubMedGoogle Scholar
- Ramsey S, Orrell D, Bolouri H: Dizzy: Stochastic simulation of large-scale genetic regulatory networks. J Bioinform Comput Biol. 2005, 3: 415-436. 10.1142/S0219720005001132View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.