Growth states of catalytic reaction networks exhibiting energy metabolism

Origin of exponential growth in nonlinear reaction networks

Exponentially growing systems are prevalent in nature, spanning all scales from biochemical reaction networks in single cells to food webs of ecosystems. How exponential growth emerges in nonlinear systems is mathematically unclear. Here, we describe a general theoretical framework that reveals underlying principles of long-term growth: scalability of flux functions and ergodicity of the rescaled systems. Our theory shows that nonlinear fluxes can generate not only balanced growth but also oscillatory or chaotic growth modalities, explaining nonequilibrium dynamics observed in cell cycles and ecosystems. Our mathematical framework is broadly useful in predicting long-term growth rates from natural and synthetic networks, analyzing the effects of system noise and perturbations, validating empirical and phenomenological laws on growth rate, and studying autocatalysis and network evolution.

Necessary and sufficient conditions for protocell growth

We consider a generic protocell model consisting of any conservative chemical reaction network embedded within a membrane. The membrane results from the self-assembly of a membrane precursor and is semi-permeable to some nutrients. Nutrients are metabolized into all other species including the membrane precursor, and the membrane grows in area and the protocell in volume. Faithful replication through cell growth and division requires a doubling of both cell volume and surface area every division time (thus leading to a periodic surface area-to-volume ratio) and also requires periodic concentrations of the cell constituents. Building upon these basic considerations, we prove necessary and sufficient conditions pertaining to the chemical reaction network for such a regime to be met. A simple necessary condition is that every moiety must be fed. A stronger necessary condition implies that every siphon must be either fed, or connected to species outside the siphon through a pass reaction capable of transferring net positive mass into the siphon. And in the case of nutrient uptake through passive diffusion and of constant surface area-to-volume ratio, a sufficient condition for the existence of a fixed point is that every siphon be fed. These necessary and sufficient conditions hold for any chemical reaction kinetics, membrane parameters or nutrient flux diffusion constants.

Enzyme oscillation can enhance the thermodynamic efficiency of cellular metabolism: consequence of anti-phase coupling between reaction flux and affinity

Cells generally convert nutrient resources to products via energy transduction. Accordingly, the thermodynamic efficiency of this conversion process is one of the most essential characteristics of living organisms. However, although these processes occur under conditions of dynamic metabolism, most studies of cellular thermodynamic efficiency have been restricted to examining steady states; thus, the relevance of dynamics to this efficiency has not yet been elucidated. Here, we develop a simple model of metabolic reactions with anabolism-catabolism coupling catalyzed by enzymes. Through application of external oscillation in the enzyme abundances, the thermodynamic efficiency of metabolism was found to be improved. This result is in strong contrast with that observed in the oscillatory input, in which the efficiency always decreased with oscillation. This improvement was effectively achieved by separating the anabolic and catabolic reactions, which tend to disequilibrate each other, and taking advantage of the temporal oscillations so that each of the antagonistic reactions could progress near equilibrium. In this case, anti-phase oscillation between the reaction flux and chemical affinity through oscillation of enzyme abundances is essential. This improvement was also confirmed in a model capable of generating autonomous oscillations in enzyme abundances. Finally, the possible relevance of the improvement in thermodynamic efficiency is discussed with respect to the potential for manipulation of metabolic oscillations in microorganisms.

Filamentation as primitive growth mode?

Osmotic pressure influences cellular shape. In a growing cell, chemical reactions and dilution induce changes in osmolarity, which in turn influence the cellular shape. Using a protocell model relying upon random conservative chemical reaction networks with arbitrary stoichiometry, we find that when the membrane is so flexible that its shape adjusts itself quasi-instantaneously to balance the osmotic pressure, the protocell either grows filamentous or fails to grow. This behavior is consistent with a mathematical proof. This suggests that filamentation may be a primitive growth mode resulting from the simple physical property of balanced osmotic pressure. We also find that growth is favored if some chemical species are only present inside the protocell, but not in the outside growth medium. Such an insulation requires specific chemical schemes. Modern evolved cells such as E. coli meet these requirements through active transport mechanisms such as the phosphotransferase system.

Minimal conditions for protocell stationary growth

We show that self-replication of a chemical system encapsulated within a membrane growing from within is possible without any explicit feature such as autocatalysis or metabolic closure, and without the need for their emergence through complexity. We use a protocell model relying upon random conservative chemical reaction networks with arbitrary stoichiometry, and we investigate the protocell's capability for self-replication, for various numbers of reactions in the network. We elucidate the underlying mechanisms in terms of simple minimal conditions pertaining only to the topology of the embedded chemical reaction network. A necessary condition is that each moiety must be fed, and a sufficient condition is that each siphon is fed. Although these minimal conditions are purely topological, by further endowing conservative chemical reaction networks with thermodynamically consistent kinetics, we show that the growth rate tends to increase on increasing the Gibbs energy per unit molecular weight of the nutrient and on decreasing that of the membrane precursor.

Entropy production of a steady-growth cell with catalytic reactions

Cells generally convert external nutrient resources to support metabolism and growth. Understanding the thermodynamic efficiency of this conversion is essential to determine the general characteristics of cellular growth. Using a simple protocell model with catalytic reaction dynamics to synthesize the necessary enzyme and membrane components from nutrients, the entropy production per unit-cell-volume growth is calculated analytically and numerically based on the rate equation for chemical kinetics and linear nonequilibrium thermodynamics. The minimal entropy production per unit-cell growth is found to be achieved at a nonzero nutrient uptake rate rather than at a quasistatic limit as in the standard Carnot engine. This difference appears because the equilibration mediated by the enzyme exists only within cells that grow through enzyme and membrane synthesis. Optimal nutrient uptake is also confirmed by protocell models with many chemical components synthesized through a catalytic reaction network. The possible relevance of the identified optimal uptake to optimal yield for cellular growth is also discussed.

Adaptation to optimal cell growth through self-organized criticality

A simple cell model consisting of a catalytic reaction network is studied to show that cellular states are self-organized in a critical state for achieving optimal growth; we consider the catalytic network dynamics over a wide range of environmental conditions, through the spontaneous regulation of nutrient transport into the cell. Furthermore, we find that the adaptability of cellular growth to reach a critical state depends only on the extent of environmental changes, while all chemical species in the cell exhibit correlated partial adaptation. These results are in remarkable agreement with the recent experimental observations of the present cells.