Abstract
(Macro)molecular crowding and the ability of the ubiquitous cytoskeleton to dynamically polymerize–depolymerize are prevalent cytoplasmic conditions in prokaryotic and eukaryotic cells. Protein interactions, enzymatic or signaling reactions - single, sequential or in complexes - whole metabolic pathways and organelles can be affected by crowding, the type and polymeric status of cytoskeletal proteins (e.g., tubulin, actin), and their imparted organization. The self-organizing capability of the cytoskeleton can orchestrate metabolic fluxes through entire pathways while its fractal organization can frame the scaling of activities in several levels of organization. The intracellular environment dynamics (e.g., biochemical reactions) is dominated by the orderly cytoskeleton and the intrinsic randomness of molecular crowding. Existing evidence underscores the inherent capacity of intracellular organization to generate emergent global behavior. Yet unknown is the relative impact on cell function provided by organelle or functional compartmentation based on transient proteins association driven by weak interactions (quinary structures) under specific environmental challenges or functional conditions (e.g., hypoxia, division, differentiation). We propose a qualitative, integrated structural–functional model of cytoplasmic organization based on a modified version of the Sierspinsky–Menger–Mandelbrot sponge, a 3D representation of a percolation cluster, and examine its capacity to accommodate established experimental facts.
It seems to me that cells leave very little to random processes and that they have evolved the capacity to escape much of the chaos of solutions, (…). It appears unlikely that a messy alphabet soup would be used to spell out the elaborate prose of intermediary metabolism.
James S. Clegg, 1984
Introduction
Cells are very far from random mixtures of molecules. The classical experiments of Kempner and Miller (1968) showed that cells are not bags of freely floating enzymes. Using cells from the unicellular eukaryote Euglena as “centrifuge tubes” (Clegg, 1984a), these authors (Kempner and Miller, 1968) could distinguish several layers within cells (that remained viable) after centrifugation. No macromolecules could be detected in the “soluble phase” and many of the main enzymes considered to exist free in solution were instead associated with layers containing organelles (mitochondria, lysosomes, nucleus) and subcellular structures (ribosomes).
Ideas about cytoplasmic organization have a long history that can be traced back to the notion of protoplasm as the substratum of cellular activity (Welch and Clegg, 2010). Early microscopic techniques (reviewed in Aon and Cortassa, 1997, Chapter 6) and the more recent cryoelectron tomography (Medalia et al., 2002) unveiled the overall crowded nature of the cellular cytoplasm populated by complex macromolecular assemblies besides subcellular organelles (Luby-Phelps, 2000; Minton, 2001; Grunewald et al., 2003; Ovádi and Norris, 2013) (see Figure 2A). Studies starting in the twentieth and well into the twenty-first centuries used advanced molecular–cellular biology methods including fluorescence recovery after photobleaching, in conjunction with confocal microscopy and time-resolved anisotropy methods, to establish a decrease in the diffusion coefficients of proteins in the cytoplasm and in the endoplasmic reticulum lumen compared with water (Luby-Phelps, 2000; Verkman, 2002; Rivas et al., 2004). The consequences of these results are very significant for cell function because diffusion of solutes and macromolecules in cellular compartments mediates many physiological processes, including metabolism and signaling events. On the other hand, active transport via motor proteins leads to a significantly higher mobility compared to diffusion processes (Fakhri et al., 2014).
Obtaining information about the molecular properties of proteins in the living cell is becoming an active field of research (Wirth and Gruebele, 2013). Existent and new molecular techniques and methodologies enable monitoring of native protein activity and folding in cells, offering information on concentration, dynamics, location, interactions and protein proximity (Diekmann and Hoischen, 2014; Fakhri et al., 2014).
Enzymatic reactions in organized crowded media
Supramolecular organization and crowding are two main traits of the intracellular milieu (Aon and Cortassa, 1997; Aon et al., 2001). The cellular cytoplasm contains (macro)molecules: a mixture of molecules of low (e.g., ATP, glutathione, NADH) and higher molecular weight (e.g., proteins, lipids, polysaccharides) and macromolecular arrays (e.g., tubulin and F-actin polymers, glycogen granules), at concentrations such that they occupy a large fraction of its total volume. Such media are “crowded” but no individual (macro)molecular species is present at a high concentration per se (Minton, 1997).
“Background” species concern (macro)molecules that do not interact specifically with either the reactants or products of a particular reaction. Proteins in the crowded cellular environment can stick to each other through non-specific interactions (e.g., electrostatic, hydrophobic). “Background” species can also contribute large steric repulsive forces in crowded environments (Minton, 2000) that may not be observed directly because they do not lead to the formation of complexes (Zhou et al., 2008). The excluded volume effect depends on the size and concentration of molecular crowders and refers to the volume between, e.g., a pair of interacting proteins, that cannot be occupied by a third protein (Wirth and Gruebele, 2013). If the size of a (macro)molecule is comparable to the size of background species, the available volume is considerably smaller (i.e., excluded volume higher) than in the case that the (macro)molecule is relatively tinier. Lower available volume increases the contribution of steric repulsion to reduce entropy and increase free energy (Rivas et al., 2004). If we consider that to maximize the available volume is a way to reduce free energy then (macro)molecular crowding facilitates a decrease in excluded (occupied) volume via, e.g., molecular compaction and association (Minton, 2000, 2001; Ellis, 2001). Another consequence of crowding is that, via decrease of excluded volume, the folded over the unfolded state of a protein or protein complex is favored (Wirth and Gruebele, 2013).
The reaction rate of an enzymatic reaction may be controlled at diffusional (substrate(s), S, access to the enzyme's, E, active site) and/or kinetic (an intrinsic step in the reaction scheme limits the rate) levels. Broadly speaking, diffusional and kinetic control may be assessed through diffusion, percolation or transport of the species S and E involved, and k2, the rate constant of the enzyme–substrate, ES, complex transformation into product according to the Henri–Michaelis–Menten (HMM) formalism (Segel, 1975).
For conceptual purpose, let us consider a simple reaction converting S into P, that involves one ligand, S, one catalytic site, E, and one enzyme–substrate complex, ES (Aon et al., 2004b):
Diffusion, percolation, transport of S and/or Ek2, also known as the catalytic rate constant, kcat, and k−1 are monomolecular rate constants whereas k1 is a bimolecular rate constant. In the derivation of the HMM equation from a quasi–steady-state assumption, the dynamics of the ES complex association–dissociation is considered to be so fast that its concentration can be treated as if it were in steady state. Accordingly, the k1 step is considered not to be limiting the S → P conversion. Importantly, from experiments performed over the past two decades, single-molecule enzymology has provided insights into how specific enzymes—particularly molecular motors and nucleic acid enzymes—work at the molecular level. These studies confirmed that the HMM mechanism expressed in Equation (1) holds at the single molecule level (Xie, 2013).
In heterogeneous, organized, media the rate of encounter between E and S may be subjected to transport restrictions generated by anomalous diffusion. Protein stickiness will increase the apparent viscosity of the cytoplasm thus decreasing the diffusion coefficient, D, since both (viscosity vs. D) are inversely related (Dix and Verkman, 2008). Anomalous diffusion introduces a time dependence in D, essentially due to the medium heterogeneity (Wirth and Gruebele, 2013) (see below Section Fractal Kinetics in Organized Crowded Media). Both translational and rotational diffusion can be influenced by the excluded volume and the shape of the crowding protein more than other factors such as hydrodynamic or direct interactions (Balbo et al., 2013). Thus, molecular crowding can affect k1 decreasing enzymatic rates, essentially because as higher the excluded volume by (macro)molecules the higher the rate limitation of the E–S encounter (Homchaudhuri et al., 2006; Pastor et al., 2014). In support of this interpretation, enzymatic reactions occurring in the presence of increasing dextran concentrations exhibited lower Vmax and higher KM, the Michaelis–Menten constant (Pastor et al., 2014). The volume occupied by dextran, independent of its size, had an important role on the initial velocity of the hydrolysis of N-succinyl-L-phenyl-Ala-p-nitroanilide catalyzed by alpha-chymotrypsin (Pastor et al., 2011). The KM increase could be attributed to a slower diffusion of the protein due to the presence of crowding, whereas the decrease in Vmax could be explained by the effect of mixed inhibition by product, which is enhanced in crowded media (Pastor et al., 2014). These results also underscore the relevant role of enzyme size in the initial velocity of reactions occurring in dextran crowded media. When enzymes are small the reaction's initial velocity mainly depends on the excluded volume. However, for large enzymes, the initial velocity of the reaction is also affected by the size of obstacles present in the environment.
Modeling and experimental work (Kim and Yethiraj, 2009; Pastor et al., 2014) also supports the idea that macromolecular crowding can contribute significantly to changes in enzymatic reactions (Vasilescu et al., 2013). To simulate rapid metabolite transfer between the enzymatic components of the phosphotransferase system (PTS), macromolecular crowding had to be assumed both to increase the association rate constants and to decrease the dissociation rate constants of the PTS complexes (Rohwer et al., 1998). However, crowding was not necessary to simulate yeast glycolysis suggesting that in this eukaryote it does not affect the glycolytic pathway (referred in Rivas et al., 2004). Cortassa and colleagues have also shown that, at least for two glycolytic enzyme couples, the effects of tubulin cytoskeleton proteins was specific and independent from crowding (Cortassa et al., 1994).
Very recent studies using carbon nanotubes for intracellular tracking of kinesin-1 motility highlight stirring dynamics as another important mode of active intracellular transport. Recorded kinesin-1 motility in COS-7 cells over five orders of magnitude in time (Fakhri et al., 2014) enabled the detection of different dynamic regimes ranging between the extremes of random thermal diffusion and kinesin-driven directed transport propelled by stirring dynamics as a non-equilibrium regime between those extremes.
Weak interactions can mediate transient protein–protein interactions collectively known as “quinary structures,” a term introduced by McConkey (1982) to define a fifth level (beyond the quaternary) of inherently transient protein structural organization. Thermodynamically, the transience of a quinary structure is based on the low stability of the molecular interaction as well as the low energetic barrier between molecular states (Wirth and Gruebele, 2013). Quinary structure has been implicated in a number of cellular processes from metabolism [e.g., the metabolon (Srere, 1987), the protein synthesis pathway (Dang et al., 1985)] to cell signaling (Li et al., 2012). Its inherent transience facilitates dynamic spatial organization of macromolecules in the cytoplasm via loose groupings of, e.g., proteins, when they are working together, but not otherwise (Wirth and Gruebele, 2013).
Overall, (macro)molecular crowding can drive molecular associations in the cytoplasm, and via modulation of the available volume or transient quinary structures influence the kinetics of biochemical reactions.
Cytoskeleton organization and metabolic fluxes
The intracellular environment is not only highly crowded but exhibits a high degree of dynamic organization governed by the principles of self-organization, as they apply to thermodynamically open non-equilibrium systems such as cells (Nicolis and Prigogine, 1977; Aon and Cortassa, 1997; de la Fuente, 2013). By exchanging energy, matter or information with their environment cells or tissues can exhibit emergence that is they self-organize their internal structure and dynamics with novel and sometimes surprising macroscopic properties. For example, the ubiquitous cytoskeletal protein network (actin or tubulin) behaves as a non-linear dissipative system, i.e., it consumes adenine nucleotides to polymerize, with the ability to self-organize, e.g., oscillate (Mandelkow et al., 1989; Mandelkow and Mandelkow, 1992; Tabony and Job, 1992) or alternate in a bistable manner (Aon et al., 1996b; Aon and Cortassa, 1997) between polymerized and depolymerized states. Structurally, the cytoskeletal network exhibits fractal properties (Mandelbrot, 1982; Feder, 1988), i.e., spatially organized in a self-similar manner thus exhibiting an alike form when observed at different degrees of magnification (Rabouille et al., 1992; Aon and Cortassa, 1994; Losa and Nonnenmacher, 1996).
The highly dynamic polymer composite of the cellular cytoplasm is dominated by protein polymers, e.g., microtubules, F-actin and intermediate filaments. In axons in the spinal cord the interaction between neurofilaments and F-actin results in a gel with particular viscoelastic properties (Leterrier et al., 1996). The nucleotide triphosphate hydrolysis-driven polymerization–depolymerization dynamics of cytoskeletal proteins is reciprocally influencing and being influenced by the myriad of reactions involved in intracellular transport (e.g., motors), metabolism, cell locomotion, and muscle contraction among other functions. For example, microtubule tracks are embedded in the viscoelastic actin cytoskeleton, which in turn fluctuates as a result of stresses generated by cytoplasmic myosins; myosin locally contracts the actin network with an attachment time of several seconds, followed by sudden release (Fakhri et al., 2014).
Cytoskeleton organization is also influenced by cytoplasmic molecular crowding because it favors protein interactions that may form modular complexes (Spirin and Mirny, 2003); some of these molecular complexes constitute stable or transient multienzyme associations (metabolon) (Srere, 1987) capable of metabolic channeling (Welch, 1977; Ovádi and Srere, 2000). In metabolic channels, reactions are facilitated by product–substrate transfer between closely associated enzymes (Ovádi and Srere, 2000).
Microtubules, actin microfilaments and intermediate filaments represent an enormous protein surface in the cell with an estimated area of 3000 μm2 for a typical mammalian cell in culture (Luby-Phelps, 2000). Thus, the cytoskeleton provides an interface for binding a variety of proteins and enzymes. The binding affinity depends on several factors, including phosphorylation status (Luther and Lee, 1986; Roberts and Somero, 1987; Pedrotti et al., 1996), the presence of other proteins, e.g., Microtubule Associated Proteins (MAPs) (Cortassa et al., 1994; Aon and Cortassa, 1997), or enzymatic activity (Cortassa and Aon, 1994; Vertessy et al., 1997; Cassimeris et al., 2012). Protein binding also affects microtubule and enzyme dynamics (Aon et al., 2001; Ovádi and Norris, 2013; Olah et al., 2015). The enhancement of metabolic flux depends upon several factors: (i) the presence for some enzymatic reactions of MAPs apart from tubulin; (ii) the concentration of microtubular protein (MTP); and (iii) the polymeric status (Cortassa et al., 1994). For example, an increase in flux through pyruvate kinase coupled to lactate dehydrogenase was elicited by MTP in a particular concentration range (Cortassa et al., 1994; Aon et al., 1996a) (Figures 1A–C; see figure legend for details). Paclitaxel and nocodazole, two drugs affecting microtubule organization and dynamics in opposite ways, were able to alter the secretion of proteolytic enzymes associated with invasion and metastasis of tumor cells (Alonso et al., 1999). While paclitaxel promotes microtubules polymerization, nocodazole elicits de-polymerization and as such they enhanced or reduced, respectively, the secretion of the urokinase-type plasminogen activator and the matrix metalloproteinase 9 to the culture medium in F3II mammary-carcinoma cells (Alonso et al., 1999).
Figure 1
Ultrasensitivity is a more sensitive response than the one expected from the classical hyperbola of Michaelis–Menten kinetics (Goldbeter and Koshland, 1982; Koshland et al., 1982). The “normal” hyperbolic response requires an 81-fold change in ligand (e.g., substrate, effector) to increase the reaction rate from 10 to 90% of the maximal velocity. Thus, ultrasensitive systems are those that need less than 81-fold change whereas sub-sensitive ones demand more. Depending on their polymeric status, cytoskeletal protein dynamics (actin, MTP) can modulate the ultrasensitive response of enzymatic systems (Aon et al., 2001). This modulation can be mediated by the cytoskeleton in the cell stress response involving changes in volume due to osmotic regulation (Busch et al., 1994; Haussinger et al., 1994a,b; Aon et al., 2000b).
More recent studies show that the heterodimer tubulin composed of α and β subunits can selectively modulate the mitochondrial outer membrane (MOM) permeability via the voltage-dependent anion channel (VDAC) (Rohwer et al., 1998; Rostovtseva and Bezrukov, 2012) thus having an impact on cellular and mitochondrial energetics (Guzun et al., 2011; Gonzalez-Granillo et al., 2012; Rostovtseva and Bezrukov, 2012). Tubulin in the nM range can influence the voltage sensitivity of VDAC reconstituted into planar phospholipid membranes, and ADP availability to the adenine nucleotide translocator in isolated mitochondria (Rostovtseva et al., 2008). Indeed, VDAC (or porin) is the most abundant protein in the MOM and is primarily involved in the ATP/ADP exchange between the cytoplasm and mitochondria (Rostovtseva and Colombini, 1997; Colombini, 2004).
Fractal kinetics in organized crowded media
The impact of self-organized cytoskeletal proteins in fractal forms upon the dynamics of cellular biochemistry started to be explored more than 20 years ago. Studying chemical reactions in heterogeneous media, Kopelman (1988) made two crucial observations: (i) reactions proceed faster in disconnected (shredded) topologies than in connected ones, and (ii) the rate constants become time-dependent, scaling with time as t−h; the h parameter reflects the characteristics of the medium, e.g., obstacle density, through the fractal dimension, and provides a link with the kinetic rate constant (Kang and Redner, 1984; Dewey, 1995; Aon et al., 2004b). In a medium like the cellular cytoplasm, molecular diffusion becomes anomalous and diffusion coefficients time-dependent due to heterogeneity given not only by (macro)molecular crowding but also by rheological changes impinging on viscosity (Forgacs and Newman, 1994; Aon and Cortassa, 1997).
Aon and Cortassa (1994) and Forgacs (1995) proposed that the cellular cytoplasm is organized as a percolation cluster. It was further conjectured that as a highly shredded object, a percolation cluster may, in principle, bestow more catalytic power to cytoplasmic enzymes. Indeed, medium organization in percolation clusters may enhance reaction rates at short times (Aon et al., 2004b; Hiroi et al., 2011) (Figure 1D). The kinetics of biochemical reactions in the cytoplasm depend on (macro)molecular crowding given by the cytoskeleton organization. Hiroi et al. (2011) showed that changing the reaction rate may be possible when the degree of intracellular macromolecular crowding is modified by experimentally manipulating the structure of the cytoskeleton. Cytoskeleton disruption with cytochalasin B and colchicine changed the anomalous diffusion parameter exhibited by enzymes and substrates/products in the cellular cytoplasm. Since the total protein concentration was maintained this suggested that not merely the concentration of intracellular proteins, but also their physiological organization profoundly affects diffusion of free molecules in a cell (Hiroi et al., 2011).
In cellular biochemistry, the fractal approach has been directed to understanding the organization and behavior of (macro)molecules in cells (Rabouille et al., 1992; Savageau, 1995; Liebovitch and Todorov, 1996; Aon et al., 2004b; Schnell and Turner, 2004; Aon and Cortassa, 2009). The approach to fractal kinetics in cells differs between authors. The dependence of the rate constant upon h has been modeled according to a fractal (Zipf–Mandelbrot) distribution (Schnell and Turner, 2004); assuming HMM kinetics in 2D lattices using Monte Carlo simulations with time-dependent rate constant (Berry, 2002), or in terms of the dependence of the parameter h on the spectral dimension, Ds, for HMM or sigmoidal kinetics (Aon et al., 2004b; Hiroi et al., 2011). Main findings show that: (i) spatial segregation of substrate and products increase with the degree of obstruction in the medium making stronger the deviation of the rate constants at longer times and, consequently, the fractal kinetic description as compared with the classical approach (Berry, 2002; Schnell and Turner, 2004); (ii) at short times the reaction rate becomes much larger, at low substrate levels, in fractal than in Euclidean space (Figure 1D); this behavior depends on the time-dependence of the KM, or an increase in cooperativity and reaction amplification in allosteric kinetics (Aon et al., 2004b). The quickly relaxing molecular mechanisms, when cells are challenged by sudden changes in environmental conditions, would provide fast and precise adaptation. Indeed, fast responses can lead to slow exhaustion processes preventing lack of substrates, effectors for reactions locally, as found in fractal media organized like percolation clusters (Hiroi et al., 2011).
Computational modeling of cytoplasmic structure–function
Interactivity in complex spatiotemporally organized systems like the cellular cytoplasm is fundamental to their counterintuitive behavior and one of the main reasons justifying the need of mathematical modeling for their study. What we seek to understand is how function is coordinated in a cell that exhibits spatially distributed heterogeneous and compartmentalized subsystems with simultaneously unfolding dynamics (Aon, 2013). Factually, (macro)molecular crowding and cytoskeleton organization (i.e., structural) are able to influence the dynamics of biochemical reactions (i.e., functional), yet how does the structural–functional coupling unfolds in physiologically meaningful spatiotemporal patterns is far from clear. One reason is that computational modeling has been mainly concerned with structural or biochemical networks but not with their integrated function. However, forced by computational burden, modeling of cytoplasmic structure–function faces some daunting challenges, beyond the complexity of the task, that demand choices between atomistic molecular-dynamic simulations (McGuffee and Elcock, 2010; Mereghetti and Wade, 2012) and lower resolution coarse-grained models (Moore et al., 2014). Mesoscale models represent a reasonable trade-off between higher simplicity (e.g., treating macromolecules as single interacting centers) while amenable to include finer biopolymer representations to address mustiscale problems of diffusion and interaction, as recently reported for the Escherichia coli cytoplasm (Trovato and Tozzini, 2014). Modeling can help decide quantitative issues such as whether moderate attraction between proteins and crowding molecules, on the order of 1 kJ/mol, can counteract the excluded volume effect (Rosen et al., 2011) or which macromolecule sizes will experience the strongest attraction and anomalous diffusion (Trovato and Tozzini, 2014).
Some insight into how structural dynamics can affect biochemical function comes from computational modeling of MTP dynamics coupled to the glycolytic pathway and its branches to the Krebs cycle, ethanolic fermentation, and the pentose phosphate (PP) pathways. This study showed that MTP dynamics can coordinately increase or decrease the flux through glycolysis, and that depending on the degree of MTP polymerization a negative control may be exerted by the PP pathway on glycolysis (Aon and Cortassa, 2002). These results may be relevant for cancer therapy because the PP pathway is critical for tumor cells to generate intermediates for nucleic acid synthesis and provide NADPH required both for the synthesis of fatty acids and cell survival under oxidative stress (Patra and Hay, 2014).
A key for progress in this complicated research field will be to adopt an experimental–modeling synergy involving iteration of the loop: simulation–validation and prediction–experimentation (Cortassa and Aon, 2013).
Emergence in subcellular organelle networks
Cytoplasmic organization comprises not only biochemical reactions but also organelles representing membrane-bound subcellular compartments. Depending on their specific function and situation, subcellular compartments such as mitochondria can play substantial roles in physiology as well as pathophysiology. In heart muscle, for example, mitochondria appear as a network in the form of a regular lattice, spanning the whole myocardial tissue like a power grid (Slodzinski et al., 2008; Aon et al., 2009). Certainly, other cytoplasmic organelles can be subject to the principles of self-organization such as the nucleus and the Golgi complex (Misteli, 2001) as well as the genome itself (Misteli, 2009), but this possibility needs experimental support.
Self-organized collective dynamics in cardiac myocytes arise from synergistically coupled subcellular networks of, e.g., mitochondria (Aon et al., 2006b) or Ca2+ release units, the latter constituted by four compartments (sarcoplasmic and junctional reticulum, myoplasmic and dyadic space) (Nivala et al., 2012). Emergence in these networks occurs through signaling via second messengers such as reactive oxygen species (ROS) or Ca2+. In addition to the normal excitation–contraction–metabolism coupling, a rich dynamic spectrum results, including oscillations, electrical or chemical waves, action potential duration alternans, early or delayed after depolarizations among others (Aon et al., 2004a; Zhou et al., 2010; Nivala et al., 2012; Qu, 2013; Zorov et al., 2014).
In heart mitochondria it was found that the transition from physiological to pathophysiological behavior happens as an emergent phenomenon of the cardiac mitochondrial network, with all the characteristics of systems at critical state (Aon et al., 2006b). This transition occurs at the percolation threshold as determined by applying percolation theory (Aon et al., 2004a). Percolation describes how local neighbor–neighbor interactions among elements in a lattice can scale to produce a macroscopic response spanning from one to the other end of a mitochondrial array (Stauffer and Aharony, 1994). Such a “spanning cluster” forms when there is a critical density of elements close to the threshold for a transition (the percolation threshold). A mitochondrial percolation cluster attains criticality at a certain threshold level of ROS (Aon et al., 2004a). The transition is self-organized, occurs with all the traits of universality—that is, with similar critical exponents as predicted by percolation theory (Schroeder, 1991; Stauffer and Aharony, 1994; Sornette, 2000), and the mitochondrial cluster exhibits fractal organization (Aon et al., 2004b). The ensuing collective oscillations—which involve at least 60% of the mitochondrial network—are synchronized by ROS via ROS-induced ROS release (Zorov et al., 2000, 2014; Aon et al., 2003; Brady et al., 2006) through a diffusion-based mechanism (Zhou et al., 2010).
In yeast, spontaneous oscillations of Saccharomyces cerevisiae mitochondrial redox states and membrane potential occur within individual yeasts (Aon et al., 2007b), and synchrony of yeast population indicates the operation of an efficient system of cell–cell interaction to produce concerted metabolic multicellular behavior (Murray et al., 2003, 2013; Lloyd and Murray, 2006, 2007; Roussel and Lloyd, 2007).
Toward an integrated structural–functional model of cytoplasmic organization
A model compatible with the 3D visualization of the cellular cytoplasmic organization as a percolation cluster (Aon and Cortassa, 1994; Forgacs, 1995) as suggested by its crowded-organized nature (Medalia et al., 2002) (Figure 2A) is the Sierpinsky–Menger–Mandelbrot sponge (Mandelbrot, 1982; Raicu and Popescu, 2008) (Figure 2B) or a modified version introduced by (Welch and Clegg, 2010) to account also for functional aspects given by “confined regions, ranging from organelles to protein complexes…, responsible for the execution of localized metabolic processes…” (Figure 2C).
Figure 2
The modified Sierpinsky–Menger–Mandelbrot sponge captures several structural–functional features discussed in this review: (1) the cytoplasmic structure represents a surface with the appearance of a 3D object, because the surface area increases as the volume shrinks. It is straightforward to imagine that the surface area will be modulated by the degree of polymerization of the cytoskeleton in turn influenced by (macro)molecular crowding. (2) It embodies the orderly cytoskeletal organization and the intrinsic randomness of (macro)molecular crowding. These two conditions can create zones of heterogeneity via exclusion (occupied) volume based on attracting (e.g., electrostatic, hydrophobic) or repulsing (e.g., steric) forces. Locally, these regions possess levels of free energy that modulate molecular association, compaction and folding/unfolding of proteins (Zhou et al., 2008; Wirth and Gruebele, 2013). (3) Is a good model for the percolation of fluids through the cytoplasm, including “confined regions” or locally separate clusters of biochemical activity that may extend to other cytoplasmic regions depending on the concentration and status of e.g., enzymatic or organelles' physiology. These localized clusters determine a certain distance with respect to the percolation threshold where local activities become global thus providing a principle of coordinated functional organization (Aon et al., 2004a; Nivala et al., 2012). (4) Enzymes or enzymatic complexes through their binding to the cytoskeleton, and substrates/effectors/messengers percolating through the “sponge” would determine reaction rates (Figure 1) defining local clusters of activity that result in different product concentrations and gradients. (5) It accommodates the existence of “quinary” structures (McConkey, 1982) that drive transiently and loosely grouped ensemble of proteins working together in a dynamic and spatially organized way, e.g., a metabolon or supramolecularly organized enzymatic complex (Srere, 1987), thus explaining compartmentalization in cytoplasmic regions that are not bounded by membranes (Wirth and Gruebele, 2013) (Figure 2C). (6) The sponge-like model of the cytoplasm is sound from the status of intracellular water (Clegg, 1984a,b) and rheological standpoints. The ground plan of living cells has been pictured as a reversible, non-covalent gel network (Luby-Phelps et al., 1986; Rabouille et al., 1990) that can be subjected to sol–gel transitions. A colloidal sol state has liquid properties with well-defined viscosity whereas in a gel viscosity becomes practically infinite with the percolation threshold, given by the concentration of polymers in solution, as the critical parameter at which the sol–gel transitions happen (Forgacs and Newman, 1994).
With a model at hand, we can now put it to test and ask questions and/or verify certain predictions. For example, can confined cytoplasmic clusters of biochemical activity be experimentally demonstrated? Could the same enzyme(s), enzymatic complex or organelle behave differently according to local cytoplasmic conditions, e.g., higher presence of F-actin (G-actin) with respect to microtubules (tubulin)? Can clusters of biochemical activity based on quinary structures be demonstrated, and are these clusters sensitive to the dynamic intracellular environment? Which intracellular conditions define percolation thresholds, and can we show local activity become globally spread? Do these conditions change with cellular stages of growth, division (e.g., G1, S phases) or differentiation? Is the status of cellular cytoplasmic organization different in differentiated with respect to pluripotent stem cells? How, and to what extent, does the intracellular dynamic field of interrelating polymeric forces play out for the dynamics of metabolic and signaling pathways interacting or being influenced by those polymers, or forming supramolecular complexes themselves? All of these are fascinating and relevant questions that now can be addressed experimentally with emerging new methodologies and a unified structural–functional theoretical framework for cytoplasmic organization.
Conflict of interest statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Statements
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AlonsoD. F.FarinaH. G.ArreguiC.AonM. A.GomezD. E. (1999). Modulation of urokinase-type plasminogen activator and metalloproteinase activities in cultured mouse mammary-carcinoma cells: enhancement by paclitaxel and inhibition by nocodazole. Int. J. Cancer83, 242–246.
2
AonM. A. (2013). Complex systems biology of networks: the riddle and the challenge, in Systems Biology of Metabolic and Signaling Networks. Energy, Mass and Information Transfer, eds AonM. A.SaksV.SchlattnerU. (Heidelberg; New York; Dordrecht; London: Springer-Verlag Berlin Heidelberg), 19–35.
3
AonM. A.CaceresA.CortassaS. (1996a). Heterogeneous distribution and organization of cytoskeletal proteins drive differential modulation of metabolic fluxes. J. Cell. Biochem. 60, 271–278.
4
AonM. A.CortassaS. (1994). On the fractal nature of cytoplasm. FEBS Lett. 344, 1–4. 10.1016/0014-5793(94)00321-1
5
AonM. A.CortassaS. (1997). Dynamic Biological Organization. Fundamentals as Applied to Cellular Systems. London: Chapman and Hall. 10.1007/978-94-011-5828-2
6
AonM. A.CortassaS. (2002). Coherent and robust modulation of a metabolic network by cytoskeletal organization and dynamics. Biophys. Chem. 97, 213–231. 10.1016/S0301-4622(02)00056-X
7
AonM. A.CortassaS. (2009). Chaotic dynamics, noise and fractal space in biochemistry, in Encyclopedia of Complexity and Systems Science, ed MeyersR. (New York, NY: Springer), 476–489.
8
AonM. A.CortassaS.AkarF. G.BrownD. A.ZhouL.O'RourkeB. (2009). From mitochondrial dynamics to arrhythmias. Int. J. Biochem. Cell Biol. 41, 1940–1948. 10.1016/j.biocel.2009.02.016
9
AonM. A.CortassaS.CaceresA. (1996b). Models of cytoplasmic structure and function, in Computation in Cellular and Molecular Biological Systems, eds CuthbertsonR.HolcombeR.PatonR. (London: World Scientific), 195–207.
10
AonM. A.CortassaS.Gomez CasatiD. F.IglesiasA. A. (2000b). Effects of stress on cellular infrastructure and metabolic organization in plant cells. Int. Rev. Cytol. 194, 239–273. 10.1016/S0074-7696(08)62398-0
11
AonM. A.CortassaS.LemarK. M.HayesA. J.LloydD. (2007b). Single and cell population respiratory oscillations in yeast: a 2-photon scanning laser microscopy study. FEBS Lett. 581, 8–14. 10.1016/j.febslet.2006.11.068
12
AonM. A.CortassaS.MarbanE.O'RourkeB. (2003). Synchronized whole cell oscillations in mitochondrial metabolism triggered by a local release of reactive oxygen species in cardiac myocytes. J. Biol. Chem. 278, 44735–44744. 10.1074/jbc.M302673200
13
AonM. A.CortassaS.O'RourkeB. (2004a). Percolation and criticality in a mitochondrial network. Proc. Natl. Acad. Sci. U.S.A. 101, 4447–4452. 10.1073/pnas.0307156101
14
AonM. A.CortassaS.O'RourkeB. (2006b). The fundamental organization of cardiac mitochondria as a network of coupled oscillators. Biophys. J. 91, 4317–4327. 10.1529/biophysj.106.087817
15
AonM. A.Gomez-CasatiD. F.IglesiasA. A.CortassaS. (2001). Ultrasensitivity in (supra)molecularly organized and crowded environments. Cell Biol. Int. 25, 1091–1099. 10.1006/cbir.2001.0804
16
AonM. A.O'RourkeB.CortassaS. (2004b). The fractal architecture of cytoplasmic organization: scaling, kinetics and emergence in metabolic networks. Mol. Cell. Biochem. 256–257, 169–184. 10.1023/B:MCBI.0000009867.54552.09
17
BalboJ.MereghettiP.HertenD. P.WadeR. C. (2013). The shape of protein crowders is a major determinant of protein diffusion. Biophys. J. 104, 1576–1584. 10.1016/j.bpj.2013.02.041
18
BerryH. (2002). Monte carlo simulations of enzyme reactions in two dimensions: fractal kinetics and spatial segregation. Biophys. J. 83, 1891–1901. 10.1016/S0006-3495(02)73953-2
19
BradyN. R.Hamacher-BradyA.WesterhoffH. V.GottliebR. A. (2006). A wave of reactive oxygen species (ROS)-induced ROS release in a sea of excitable mitochondria. Antioxid. Redox Signal. 8, 1651–1665. 10.1089/ars.2006.8.1651
20
BuschG. L.SchreiberR.DartschP. C.VolklH.Vom DahlS.HaussingerD.et al. (1994). Involvement of microtubules in the link between cell volume and pH of acidic cellular compartments in rat and human hepatocytes. Proc. Natl. Acad. Sci. U.S.A. 91, 9165–9169. 10.1073/pnas.91.19.9165
21
CassimerisL.SilvaV. C.MillerE.TonQ.MolnarC.FongJ. (2012). Fueled by microtubules: does tubulin dimer/polymer partitioning regulate intracellular metabolism?Cytoskeleton (Hoboken). 69, 133–143. 10.1002/cm.21008
22
CleggJ. S. (1984a). Properties and metabolism of the aqueous cytoplasm and its boundaries. Am. J. Physiol. 246, R133–R151.
23
CleggJ. S. (1984b). Intracellular water and the cytomatrix: some methods of study and current views. J. Cell Biol. 99, 167s–171s. 10.1083/jcb.99.1.167s
24
ColombiniM. (2004). VDAC: the channel at the interface between mitochondria and the cytosol. Mol. Cell. Biochem. 256–257, 107–115. 10.1023/B:MCBI.0000009862.17396.8d
25
CortassaS.AonM. A. (1994). Spatio-temporal regulation of glycolysis and oxidative phosphorylation in vivo in tumor and yeast cells. Cell Biol. Int. 18, 687–713. 10.1006/cbir.1994.1099
26
CortassaS.AonM. A. (2013). Dynamics of mitochondrial redox and energy networks: insights from an experimental-computational synergy, in Systems Biology of Metabolic and Signaling Networks. Energy, Mass and Information Transfer, eds AonM. A.SaksV.SchlattnerU. (Heidelberg; New York; Dordrecht; London: Springer-Verlag Berlin Heidelberg), 115–144.
27
CortassaS.CaceresA.AonM. A. (1994). Microtubular protein in its polymerized or nonpolymerized states differentially modulates in vitro and intracellular fluxes catalyzed by enzymes of carbon metabolism. J. Cell. Biochem. 55, 120–132. 10.1002/jcb.240550114
28
DangC. V.FergusonB.BurkeD. J.GarciaV.YangD. C. (1985). Interactions of aminoacyl-tRNA synthetases in high-molecular-weight multienzyme complexes from rat liver. Biochim. Biophys. Acta829, 319–326. 10.1016/0167-4838(85)90239-0
29
de la FuenteI. M. (2013). Metabolic dissipative structures, in Systems Biology of Metabolic and Signaling Networks. Energy, Mass and Information Transfer, eds AonM. A.SaksV.SchlattnerU. (Heidelberg; New York; Dordrecht; London: Springer-Verlag Berlin Heidelberg), 179–211.
30
DeweyT. G. (1995). Chemically-controlled reaction kinetics on fractals: application to hydrogen exchange in lysozyme. Fractals3, 251–267. 10.1142/S0218348X95000217
31
DiekmannS.HoischenC. (2014). Biomolecular dynamics and binding studies in the living cell. Phys. Life Rev. 11, 1–30. 10.1016/j.plrev.2013.11.011
32
DixJ. A.VerkmanA. S. (2008). Crowding effects on diffusion in solutions and cells. Annu. Rev. Biophys. 37, 247–263. 10.1146/annurev.biophys.37.032807.125824
33
EllisR. J. (2001). Macromolecular crowding: an important but neglected aspect of the intracellular environment. Curr. Opin. Struct. Biol. 11, 114–119. 10.1016/S0959-440X(00)00172-X
34
FakhriN.WesselA. D.WillmsC.PasqualiM.KlopfensteinD. R.MacKintoshF. C.et al. (2014). High-resolution mapping of intracellular fluctuations using carbon nanotubes. Science344, 1031–1035. 10.1126/science.1250170
35
FederJ. (1988). Fractals. New York, NY: Plenum Press. 10.1007/978-1-4899-2124-6
36
ForgacsG. (1995). On the possible role of cytoskeletal filamentous networks in intracellular signaling: an approach based on percolation. J. Cell Sci. 108(Pt 6), 2131–2143.
37
ForgacsG.NewmanS. A. (1994). Phase transitions, interfaces, and morphogenesis in a network of protein fibers. Int. Rev. Cytol. 150, 139–148. 10.1016/S0074-7696(08)61540-5
38
GoldbeterA.KoshlandD. E.Jr. (1982). Sensitivity amplification in biochemical systems. Q. Rev. Biophys. 15, 555–591. 10.1017/S0033583500003449
39
Gonzalez-GranilloM.GrichineA.GuzunR.UssonY.TeppK.ChekulayevV.et al. (2012). Studies of the role of tubulin beta II isotype in regulation of mitochondrial respiration in intracellular energetic units in cardiac cells. J. Mol. Cell. Cardiol. 52, 437–447. 10.1016/j.yjmcc.2011.07.027
40
GrunewaldK.MedaliaO.GrossA.StevenA. C.BaumeisterW. (2003). Prospects of electron cryotomography to visualize macromolecular complexes inside cellular compartments: implications of crowding. Biophys. Chem. 100, 577–591. 10.1016/S0301-4622(02)00307-1
41
GuzunR.Karu-VarikmaaM.Gonzalez-GranilloM.KuznetsovA. V.MichelL.Cottet-RousselleC.et al. (2011). Mitochondria-cytoskeleton interaction: distribution of beta-tubulins in cardiomyocytes and HL-1 cells. Biochim. Biophys. Acta1807, 458–469. 10.1016/j.bbabio.2011.01.010
42
HaussingerD.LangF.GerokW. (1994a). Regulation of cell function by the cellular hydration state. Am. J. Physiol. 267, E343–E355.
43
HaussingerD.StollB.vom DahlS.TheodoropoulosP. A.MarkogiannakisE.GravanisA.et al. (1994b). Effect of hepatocyte swelling on microtubule stability and tubulin mRNA levels. Biochem. Cell Biol. 72, 12–19. 10.1139/o94-003
44
HiroiN.LuJ.IbaK.TabiraA.YamashitaS.OkadaY.et al. (2011). Physiological environment induces quick response – slow exhaustion reactions. Front. Physiol. 2:50. 10.3389/fphys.2011.00050
45
HomchaudhuriL.SarmaN.SwaminathanR. (2006). Effect of crowding by dextrans and Ficolls on the rate of alkaline phosphatase-catalyzed hydrolysis: a size-dependent investigation. Biopolymers83, 477–486. 10.1002/bip.20578
46
KangK.RednerS. (1984). Scaling approach for the kinetics of recombination processes. Phys. Rev. Lett. 52, 955–958. 10.1103/PhysRevLett.52.955
47
KempnerE. S.MillerJ. H. (1968). The molecular biology of Euglena gracilis. IV. Cellular stratification by centrifuging. Exp. Cell Res. 51, 141–149. 10.1016/0014-4827(68)90164-X
48
KimJ. S.YethirajA. (2009). Effect of macromolecular crowding on reaction rates: a computational and theoretical study. Biophys. J. 96, 1333–1340. 10.1016/j.bpj.2008.11.030
49
KopelmanR. (1988). Fractal reaction kinetics. Science241, 1620–1626. 10.1126/science.241.4873.1620
50
KoshlandD. E.Jr.GoldbeterA.StockJ. B. (1982). Amplification and adaptation in regulatory and sensory systems. Science217, 220–225. 10.1126/science.7089556
51
LeterrierJ. F.KasJ.HartwigJ.VegnersR.JanmeyP. A. (1996). Mechanical effects of neurofilament cross-bridges. Modulation by phosphorylation, lipids, and interactions with F-actin. J. Biol. Chem. 271, 15687–15694. 10.1074/jbc.271.26.15687
52
LiP.BanjadeS.ChengH. C.KimS.ChenB.GuoL.et al. (2012). Phase transitions in the assembly of multivalent signalling proteins. Nature483, 336–340. 10.1038/nature10879
53
LiebovitchL. S.TodorovA. T. (1996). Using fractals and nonlinear dynamics to determine the physical properties of ion channel proteins. Crit. Rev. Neurobiol. 10, 169–187. 10.1615/CritRevNeurobiol.v10.i2.20
54
LloydD.MurrayD. B. (2006). The temporal architecture of eukaryotic growth. FEBS Lett. 580, 2830–2835. 10.1016/j.febslet.2006.02.066
55
LloydD.MurrayD. B. (2007). Redox rhythmicity: clocks at the core of temporal coherence. Bioessays29, 465–473. 10.1002/bies.20575
56
LosaG. A.NonnenmacherT. F. (1996). Self-similarity and fractal irregularity in pathologic tissues. Mod. Pathol. 9, 174–182.
57
Luby-PhelpsK. (2000). Cytoarchitecture and physical properties of cytoplasm: volume, viscosity, diffusion, intracellular surface area. Int. Rev. Cytol. 192, 189–221. 10.1016/S0074-7696(08)60527-6
58
Luby-PhelpsK.TaylorD. L.LanniF. (1986). Probing the structure of cytoplasm. J. Cell Biol. 102, 2015–2022. 10.1083/jcb.102.6.2015
59
LutherM. A.LeeJ. C. (1986). The role of phosphorylation in the interaction of rabbit muscle phosphofructokinase with F-actin. J. Biol. Chem. 261, 1753–1759.
60
MandelbrotB. B. (1982). The Fractal Geometry of Nature. New York, NY: W.H. Freeman.
61
MandelkowE.MandelkowE. M.HotaniH.HessB.MullerS. C. (1989). Spatial patterns from oscillating microtubules. Science246, 1291–1293. 10.1126/science.2588005
62
MandelkowE. M.MandelkowE. (1992). Microtubule oscillations. Cell Motil. Cytoskeleton22, 235–244. 10.1002/cm.970220403
63
McConkeyE. H. (1982). Molecular evolution, intracellular organization, and the quinary structure of proteins. Proc. Natl. Acad. Sci. U.S.A. 79, 3236–3240. 10.1073/pnas.79.10.3236
64
McGuffeeS. R.ElcockA. H. (2010). Diffusion, crowding & protein stability in a dynamic molecular model of the bacterial cytoplasm. PLoS Comput. Biol. 6:e1000694. 10.1371/journal.pcbi.1000694
65
MedaliaO.WeberI.FrangakisA. S.NicastroD.GerischG.BaumeisterW. (2002). Macromolecular architecture in eukaryotic cells visualized by cryoelectron tomography. Science298, 1209–1213. 10.1126/science.1076184
66
MereghettiP.WadeR. C. (2012). Atomic detail brownian dynamics simulations of concentrated protein solutions with a mean field treatment of hydrodynamic interactions. J. Phys. Chem. B116, 8523–8533. 10.1021/jp212532h
67
MintonA. P. (1997). Influence of excluded volume upon macromolecular structure and associations in “crowded” media. Curr. Opin. Biotechnol. 8, 65–69. 10.1016/S0958-1669(97)80159-0
68
MintonA. P. (2000). Implications of macromolecular crowding for protein assembly. Curr. Opin. Struct. Biol. 10, 34–39. 10.1016/S0959-440X(99)00045-7
69
MintonA. P. (2001). The influence of macromolecular crowding and macromolecular confinement on biochemical reactions in physiological media. J. Biol. Chem. 276, 10577–10580. 10.1074/jbc.R100005200
70
MisteliT. (2001). The concept of self-organization in cellular architecture. J. Cell Biol. 155, 181–185. 10.1083/jcb.200108110
71
MisteliT. (2009). Self-organization in the genome. Proc. Natl. Acad. Sci. U.S.A. 106, 6885–6886. 10.1073/pnas.0902010106
72
MooreT.WuS. K.MichaelM.YapA. S.GomezG. A.NeufeldZ. (2014). Self-organizing actomyosin patterns on the cell cortex at epithelial cell-cell junctions. Biophys. J. 107, 2652–2661. 10.1016/j.bpj.2014.10.045
73
MurrayD. B.AmarielC.SasidharanK.MachneR.AonM. A.LloydD. (2013). Temporal partitioning of the yeast cellular network, in Systems Biology of Metabolic and Signaling Networks. Energy, Mass and Information Transfer, eds AonM. A.SaksV.SchlattnerU. (Heidelberg; New York; Dordrecht; London: Springer-Verlag Berlin Heidelberg), 323–349.
74
MurrayD. B.KleveczR. R.LloydD. (2003). Generation and maintenance of synchrony in Saccharomyces cerevisiae continuous culture. Exp. Cell Res. 287, 10–15. 10.1016/S0014-4827(03)00068-5
75
NicolisG.PrigogineI. (1977). Self-organization in Nonequilibrium Systems: From Dissipative Structures to Order Through Fluctuations. New York, NY: Wiley.
76
NivalaM.KoC. Y.WeissJ. N.QuZ. (2012). Criticality in intracellular calcium signaling in cardiac myocytes. Biophys. J. 102, 2433–2442. 10.1016/j.bpj.2012.05.001
77
OlahJ.NorrisV.OvádiJ. (2015). Modeling of sensing potency of cytoskeletal systems decorated with metabolic enzymes. J. Theor. Biol. 365, 190–196. 10.1016/j.jtbi.2014.10.018
78
OvádiJ.NorrisV. (2013). Moonlighting function of the tubulin cytoskeleton: macromolecular architectures in the cytoplasm, in Systems Biology of Metabolic and Signaling Networks. Energy, Mass and Information Transfer, eds AonM. A.SaksV.SchlattnerU. (Heidelberg; New York; Dordrecht; London: Springer-Verlag Berlin Heidelberg), 165–178.
79
OvádiJ.SrereP. A. (2000). Macromolecular compartmentation and channeling. Int. Rev. Cytol. 192, 255–280. 10.1016/S0074-7696(08)60529-X
80
PastorI.PituliceL.BalcellsC.VilasecaE.MadurgaS.IsvoranA.et al. (2014). Effect of crowding by Dextrans in enzymatic reactions. Biophys. Chem. 185, 8–13. 10.1016/j.bpc.2013.10.006
81
PastorI.VilasecaE.MadurgaS.GarcesJ. L.CascanteM.MasF. (2011). Effect of crowding by dextrans on the hydrolysis of N-Succinyl-L-phenyl-Ala-p-nitroanilide catalyzed by alpha-chymotrypsin. J. Phys. Chem. B115, 1115–1121. 10.1021/jp105296c
82
PatraK. C.HayN. (2014). The pentose phosphate pathway and cancer. Trends Biochem. Sci. 39, 347–354. 10.1016/j.tibs.2014.06.005
83
PedrottiB.UlloaL.AvilaJ.IslamK. (1996). Characterization of microtubule-associated protein MAP1B: phosphorylation state, light chains, and binding to microtubules. Biochemistry35, 3016–3023. 10.1021/bi951314f
84
QuZ. (2013). Network dynamics in cardiac electrophysiology, in Systems Biology of Metabolic and Signaling Networks. Energy, Mass and Information Transfer, eds AonM. A.SaksV.SchlattnerU. (Heidelberg; New York; Dordrecht; London: Springer-Verlag Berlin Heidelberg), 243–260.
85
RabouilleC.AonM. A.MullerG.CartaudJ.ThomasD. (1990). The supramolecular organization of ovomucin. Biophysical and morphological studies. Biochem. J. 266, 697–706.
86
RabouilleC.CortassaS.AonM. A. (1992). Fractal organisation in biological macromolecular lattices. J. Biomol. Struct. Dyn. 9, 1013–1024. 10.1080/07391102.1992.10507973
87
RaicuV.PopescuA. (2008). Integrated Molecular and Cellular Biophysics. New York, NY: Springer. 10.1007/978-1-4020-8268-9
88
RivasG.FerroneF.HerzfeldJ. (2004). Life in a crowded world. EMBO Rep. 5, 23–27. 10.1038/sj.embor.7400056
89
RobertsS. J.SomeroG. N. (1987). Binding of phosphofructokinase to filamentous actin. Biochemistry26, 3437–3442. 10.1021/bi00386a028
90
RohwerJ. M.PostmaP. W.KholodenkoB. N.WesterhoffH. V. (1998). Implications of macromolecular crowding for signal transduction and metabolite channeling. Proc. Natl. Acad. Sci. U.S.A. 95, 10547–10552. 10.1073/pnas.95.18.10547
91
RosenJ.KimY. C.MittalJ. (2011). Modest protein-crowder attractive interactions can counteract enhancement of protein association by intermolecular excluded volume interactions. J. Phys. Chem. B115, 2683–2689. 10.1021/jp200625k
92
RostovtsevaT.ColombiniM. (1997). VDAC channels mediate and gate the flow of ATP: implications for the regulation of mitochondrial function. Biophys. J. 72, 1954–1962. 10.1016/S0006-3495(97)78841-6
93
RostovtsevaT. K.BezrukovS. M. (2012). VDAC inhibition by tubulin and its physiological implications. Biochim. Biophys. Acta1818, 1526–1535. 10.1016/j.bbamem.2011.11.004
94
RostovtsevaT. K.SheldonK. L.HassanzadehE.MongeC.SaksV.BezrukovS. M.et al. (2008). Tubulin binding blocks mitochondrial voltage-dependent anion channel and regulates respiration. Proc. Natl. Acad. Sci. U.S.A. 105, 18746–18751. 10.1073/pnas.0806303105
95
RousselM. R.LloydD. (2007). Observation of a chaotic multioscillatory metabolic attractor by real-time monitoring of a yeast continuous culture. FEBS J. 274, 1011–1018. 10.1111/j.1742-4658.2007.05651.x
96
SavageauM. A. (1995). Michaelis-Menten mechanism reconsidered: implications of fractal kinetics. J. Theor. Biol. 176, 115–124. 10.1006/jtbi.1995.0181
97
SchnellS.TurnerT. E. (2004). Reaction kinetics in intracellular environments with macromolecular crowding: simulations and rate laws. Prog. Biophys. Mol. Biol. 85, 235–260. 10.1016/j.pbiomolbio.2004.01.012
98
SchroederM. (1991). Fractals, Chaos, Power Laws. Minutes from an Infinite Paradise. New York, NY: W.H. Freeman and Company.
99
SegelI. H. (1975). Enzyme Kinetics. Behavior and Analysis of Rapid Equilibrium and Steady State Enzyme Systems. New York, NY: Wiley-Interscience Publications.
100
SlodzinskiM. K.AonM. A.O'RourkeB. (2008). Glutathione oxidation as a trigger of mitochondrial depolarization and oscillation in intact hearts. J. Mol. Cell. Cardiol. 45, 650–660. 10.1016/j.yjmcc.2008.07.017
101
SornetteD. (2000). Critical Phenomena in Natural Sciences. Chaos, Fractals, Selforganization and Disorder: Concepts and Tools. Berlin; Heidelberg: Springer-Verlag. 10.1007/978-3-662-04174-1
102
SpirinV.MirnyL. A. (2003). Protein complexes and functional modules in molecular networks. Proc. Natl. Acad. Sci. U.S.A. 100, 12123–12128. 10.1073/pnas.2032324100
103
SrereP. A. (1987). Complexes of sequential metabolic enzymes. Annu. Rev. Biochem. 56, 89–124. 10.1146/annurev.bi.56.070187.000513
104
StaufferD.AharonyA. (1994). Introduction To Percolation Theory. London: Taylor and Francis.
105
TabonyJ.JobD. (1992). Gravitational symmetry breaking in microtubular dissipative structures. Proc. Natl. Acad. Sci. U.S.A. 89, 6948–6952. 10.1073/pnas.89.15.6948
106
TrovatoF.TozziniV. (2014). Diffusion within the cytoplasm: a mesoscale model of interacting macromolecules. Biophys. J. 107, 2579–2591. 10.1016/j.bpj.2014.09.043
107
VasilescuC.OlteanuM.FlondorP.CalinG. A. (2013). Fractal-like kinetics of intracellular enzymatic reactions: a chemical framework of endotoxin tolerance and a possible non-specific contribution of macromolecular crowding to cross-tolerance. Theor. Biol. Med. Model. 10, 55. 10.1186/1742-4682-10-55
108
VerkmanA. S. (2002). Solute and macromolecule diffusion in cellular aqueous compartments. Trends Biochem. Sci. 27, 27–33. 10.1016/S0968-0004(01)02003-5
109
VertessyB. G.OroszF.KovacsJ.OvádiJ. (1997). Alternative binding of two sequential glycolytic enzymes to microtubules. Molecular studies in the phosphofructokinase/aldolase/microtubule system. J. Biol. Chem. 272, 25542–25546. 10.1074/jbc.272.41.25542
110
WelchG. R. (1977). On the role of organized multienzyme systems in cellular metabolism: a general synthesis. Prog. Biophys. Mol. Biol. 32, 103–191. 10.1016/0079-6107(78)90019-6
111
WelchG. R.CleggJ. S. (2010). From protoplasmic theory to cellular systems biology: a 150-year reflection. Am. J. Physiol. Cell Physiol. 298, C1280–C1290. 10.1152/ajpcell.00016.2010
112
WirthA. J.GruebeleM. (2013). Quinary protein structure and the consequences of crowding in living cells: leaving the test-tube behind. Bioessays35, 984–993. 10.1002/bies.201300080
113
XieX. S. (2013). Biochemistry. Enzyme kinetics, past and present. Science342, 1457–1459. 10.1126/science.1248859
114
ZhouH. X.RivasG.MintonA. P. (2008). Macromolecular crowding and confinement: biochemical, biophysical, and potential physiological consequences. Annu. Rev. Biophys. 37, 375–397. 10.1146/annurev.biophys.37.032807.125817
115
ZhouL.AonM. A.AlmasT.CortassaS.WinslowR. L.O'RourkeB. (2010). A reaction-diffusion model of ROS-induced ROS release in a mitochondrial network. PLoS Comput. Biol. 6:e1000657. 10.1371/journal.pcbi.1000657
116
ZorovD. B.FilburnC. R.KlotzL. O.ZweierJ. L.SollottS. J. (2000). Reactive oxygen species (ROS)-induced ROS release: a new phenomenon accompanying induction of the mitochondrial permeability transition in cardiac myocytes. J. Exp. Med. 192, 1001–1014. 10.1084/jem.192.7.1001
117
ZorovD. B.JuhaszovaM.SollottS. J. (2014). Mitochondrial reactive oxygen species (ROS) and ROS-induced ROS release. Physiol. Rev. 94, 909–950. 10.1152/physrev.00026.2013
Summary
Keywords
enzyme kinetics, metabolism, quinary structures, cytoskeleton, molecular crowding, fractal, Sierpinsky sponge, percolation
Citation
Aon MA and Cortassa S (2015) Function of metabolic and organelle networks in crowded and organized media. Front. Physiol. 5:523. doi: 10.3389/fphys.2014.00523
Received
15 October 2014
Accepted
19 December 2014
Published
21 January 2015
Volume
5 - 2014
Edited by
Noriko Hiroi, Keio University, Japan
Reviewed by
Julio Vera González, University Hospital Erlangen, Germany; Akatsuki Kimura, National Institute of Genetics, Japan
Copyright
© 2015 Aon and Cortassa.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Miguel A. Aon, Department of Medicine, School of Medicine, Johns Hopkins University, 720 Rutland Avenue, Ross Bldg. 1059, Baltimore, MD 21205, USA e-mail: maon1@jhmi.edu
This article was submitted to Systems Biology, a section of the journal Frontiers in Physiology.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.