Gears in chemical reaction networks for optimizing energy transduction efficiency

Abstract

Abstract Similarly to gear systems in vehicles, most chemical reaction networks (CRNs) involved in energy transduction have at their disposal multiple transduction pathways, each characterized by distinct efficiencies. We conceptualize these pathways as ‘chemical gears’ and demonstrate their role in refining the second law of thermodynamics. This allows us to determine the optimal efficiency of a CRN, and the gear enabling it, solely based on its topology and operating conditions, defined by the chemical potentials of its input and output species. By suitably tuning reaction kinetics, a CRN can be engineered to self-regulate its gear settings, maintaining optimal efficiency under varying external conditions. We demonstrate this principle in a biological context with a CRN where enzymes function as gear shifters, autonomously adapting the system to achieve near-optimal efficiency across changing environments. Additionally, we analyze the gear system of an artificial molecular motor, identifying numerous counterproductive gears and providing insights into its transduction capabilities and optimization.

Locations

  • Nature Communications
  • arXiv (Cornell University)
  • PubMed

Ask a Question About This Paper

Summary

Chemical Reaction Networks (CRNs) that transduce free energy, often termed chemical machines, are fundamental to both biological processes and the design of artificial molecular devices. While conventional thermodynamics, particularly the second law, establishes broad limits on their efficiency, this work introduces a refined framework that significantly deepens the understanding of how these machines can achieve and maintain optimal performance.

The central innovation is the formal concept of “chemical gears.” Analogous to mechanical gears in vehicles, these chemical gears represent distinct energy transduction pathways within a CRN, each characterized by a unique efficiency. These gears are rigorously derived from the well-established mathematical tool of Elementary Flux Modes (EFMs), a concept previously widely used in metabolic network analysis to identify fundamental reaction pathways. By conceptualizing these EFMs as chemical gears, the paper demonstrates that the maximum achievable efficiency of a CRN, along with the specific gear enabling it, can be determined solely from the network’s topological structure and the chemical potentials of its input and output species (its operating conditions). This provides a tighter, more informative upper bound on efficiency than the general limits imposed by the second law of thermodynamics.

A key implication of this framework is the potential for CRNs to self-regulate their gearing. The paper demonstrates that by suitably tuning reaction kinetics—for instance, through enzyme regulation in biological systems—a CRN can autonomously switch between different gears to maintain optimal efficiency even as external conditions vary. This mechanism mirrors metabolic switching observed in living organisms, where they adapt their metabolism to changing environmental demands. The work presents a biologically inspired CRN model where enzymes function as “gear shifters,” allowing the system to adapt and achieve near-optimal efficiency across varying conditions, highlighting a trade-off between power and efficiency.

Beyond biology, this gear-based perspective offers a powerful benchmark for assessing the efficiency of chemical machines and provides a novel avenue for designing more efficient artificial molecular motors. The paper applies its framework to existing models of artificial molecular motors, revealing that many currently operate far from optimally due to the contribution of “futile” or suboptimal gears. This suggests that by implementing sophisticated gear regulation, similar to biological systems, the performance of synthetic nanomachines could be substantially enhanced.

The main prior ingredients underpinning this research include the fundamental theory of open Chemical Reaction Networks, particularly their steady-state properties and the principles of free energy transduction, building upon seminal works in non-equilibrium thermodynamics. Crucially, the mathematical foundation of Elementary Flux Modes (EFMs), originating from constraint-based metabolic modeling, is directly leveraged and re-interpreted as chemical gears. Concepts from enzyme kinetics and regulation are essential for illustrating the self-regulation mechanisms. Finally, the broader field of molecular motors, both biological and artificial, provides the practical context and application domain for these theoretical advancements. The paper also uses the mathematical concept of “conformal vectors” to decompose steady-state fluxes, which aids in proving the upper bound on transduction efficiency.

Abstract Light‐induced processes are the basis of many fundamental natural phenomena as well as of a variety of applications. Since the functions that can arise from the interaction between light … Abstract Light‐induced processes are the basis of many fundamental natural phenomena as well as of a variety of applications. Since the functions that can arise from the interaction between light and matter depend on the degree of complexity and organization of the receiving “matter,” the research on these processes has progressively moved from molecular to supramolecular (multicomponent) systems, thereby originating the field of supramolecular photochemistry. In this context, several research groups have synthesized and investigated multicomponent chemical systems capable of performing specific light‐induced functions, such as elaboration of information in the form of input/output signals and relative mechanical motions of the molecular components. Systems of this type can be viewed as simple examples of molecular devices and machines. These studies are of interest not only for increasing the basic understanding of photoinduced processes but also for the growth of nanoscience and the development of nanotechnology.
One of the most fascinating properties of the biotechnologically important organism Saccharomyces cerevisiae is its ability to perform simultaneous respiration and fermentation at high growth rate even under fully aerobic … One of the most fascinating properties of the biotechnologically important organism Saccharomyces cerevisiae is its ability to perform simultaneous respiration and fermentation at high growth rate even under fully aerobic conditions. In the present work, this Crabtree effect called phenomenon was investigated in detail by comparative 13C metabolic flux analysis of S. cerevisiae growing under purely oxidative, respiro-fermentative and predominantly fermentative conditions.The metabolic shift from oxidative to fermentative growth was accompanied by complex changes of carbon flux throughout the whole central metabolism. This involved a flux redirection from the pentose phosphate pathway (PPP) towards glycolysis, an increased flux through pyruvate carboxylase, the fermentative pathways and malic enzyme, a flux decrease through the TCA cycle, and a partial relocation of alanine biosynthesis from the mitochondrion to the cytosol. S. cerevisiae exhibited a by-pass of pyruvate dehydrogenase in all physiological regimes. During oxidative growth this by-pass was mainly provided via pyruvate decarboxylase, acetaldehyde dehydrogenase, acetyl-CoA synthase and transport of acetyl-CoA into the mitochondrion. During fermentative growth this route, however, was saturated due to limited enzyme capacity. Under these conditions the cells exhibited high carbon flux through a chain of reactions involving pyruvate carboxylase, the oxaloacetate transporter and malic enzyme. During purely oxidative growth the PPP alone was sufficient to completely supply NADPH for anabolism. During fermentation, it provided only 60 % of the required NADPH.We conclude that, in order to overcome the limited capacity of pyruvate dehydrogenase, S. cerevisiae possesses different metabolic by-passes to channel carbon into the mitochondrion. This involves the conversion of cytosolic pyruvate either into acetyl CoA or oxaloacetate followed by intercompartmental transport of these metabolites. During oxidative growth mainly the NAD specific isoforms of acetaldehyde dehydrogenase and isocitrate dehydrogenase catalyze the corresponding reactions in S. cerevisiae, whereas NADPH supply under fermentative conditions involves significant contribution of sources other than the PPP such as e. g. NADPH specific acetaldehyde dehydrogenase or isocitrate dehydrogenase.
Abstract Elementary flux mode (EFM) analysis allows the unbiased decomposition of a metabolic network into minimal functional units, making it a powerful tool for metabolic engineering. While the use of … Abstract Elementary flux mode (EFM) analysis allows the unbiased decomposition of a metabolic network into minimal functional units, making it a powerful tool for metabolic engineering. While the use of EFM analysis (EFMA) is still limited by the size of the models it can handle, EFMA has been successfully applied to solve real‐world metabolic engineering problems. Here we provide a user‐oriented introduction to EFMA, provide examples of recent applications, analyze current research strategies to overcome the computational restrictions and give an overview over current approaches, which aim to identify and calculate only biologically relevant EFMs.
Elementary flux modes (EFMs)--non-decomposable minimal pathways--are commonly accepted tools for metabolic network analysis under steady state conditions. Valid states of the network are linear superpositions of elementary modes shaping a … Elementary flux modes (EFMs)--non-decomposable minimal pathways--are commonly accepted tools for metabolic network analysis under steady state conditions. Valid states of the network are linear superpositions of elementary modes shaping a polyhedral cone (the flux cone), which is a well-studied convex set in computational geometry. Computing EFMs is thus basically equivalent to extreme ray enumeration of polyhedral cones. This is a combinatorial problem with poorly scaling algorithms, preventing the large-scale analysis of metabolic networks so far.Here, we introduce new algorithmic concepts that enable large-scale computation of EFMs. Distinguishing extreme rays from normal (composite) vectors is one critical aspect of the algorithm. We present a new recursive enumeration strategy with bit pattern trees for adjacent rays--the ancestors of extreme rays--that is roughly one order of magnitude faster than previous methods. Additionally, we introduce a rank updating method that is particularly well suited for parallel computation and a residue arithmetic method for matrix rank computations, which circumvents potential numerical instability problems. Multi-core architectures of modern CPUs can be exploited for further performance improvements. The methods are applied to a central metabolism network of Escherichia coli, resulting in approximately 26 Mio. EFMs. Within the top 2% modes considering biomass production, most of the gain in flux variability is achieved. In addition, we compute approximately 5 Mio. EFMs for the production of non-essential amino acids for a genome-scale metabolic network of Helicobacter pylori. Only large-scale EFM analysis reveals the >85% of modes that generate several amino acids simultaneously.An implementation in Java, with integration into MATLAB and support of various input formats, including SBML, is available at http://www.csb.ethz.ch in the tools section; sources are available from the authors upon request.
The authors present general considerations and simple models for the operation of isothermal motors at small scales, in asymmetric environments. Their work is inspired by recent observations on the behavior … The authors present general considerations and simple models for the operation of isothermal motors at small scales, in asymmetric environments. Their work is inspired by recent observations on the behavior of molecular motors in the biological realm, where chemical energy is converted into mechanical energy. A generic Onsager-like description of the linear (close to equilibrium) regime is presented, which exhibits structural differences from the usual Carnot engines. Turning to more explicit models for a single motor, the authors show the importance of the time scales involved and of the spatial dependence of the motor's chemical activity. Considering the situation in which a large collection of such motors operates together. The authors exhibit new features among which are dynamical phase transitions formally similar to paramagnetic-ferromagnetic and liquid-vapor transitions.
Metabolically versatile free-living bacteria have global regulation systems that allow cells to selectively assimilate a preferred compound among a mixture of several potential carbon sources. This process is known as … Metabolically versatile free-living bacteria have global regulation systems that allow cells to selectively assimilate a preferred compound among a mixture of several potential carbon sources. This process is known as carbon catabolite repression (CCR). CCR optimizes metabolism, improving the ability of bacteria to compete in their natural habitats. This review summarizes the regulatory mechanisms responsible for CCR in the bacteria of the genus Pseudomonas, which can live in many different habitats. Although the information available is still limited, the molecular mechanisms responsible for CCR in Pseudomonas are clearly different from those of Enterobacteriaceae or Firmicutes. An understanding of the molecular mechanisms underlying CCR is important to know how metabolism is regulated and how bacteria degrade compounds in the environment. This is particularly relevant for compounds that are degraded slowly and accumulate, creating environmental problems. CCR has a major impact on the genes involved in the transport and metabolism of nonpreferred carbon sources, but also affects the expression of virulence factors in several bacterial species, genes that are frequently directed to allow the bacterium to gain access to new sources of nutrients. Finally, CCR has implications in the optimization of biotechnological processes such as biotransformations or bioremediation strategies.
Motor proteins are nature's solution for directing movement at the molecular level. The field of artificial molecular motors takes inspiration from these tiny but powerful machines. Although directional motion on … Motor proteins are nature's solution for directing movement at the molecular level. The field of artificial molecular motors takes inspiration from these tiny but powerful machines. Although directional motion on the nanoscale performed by synthetic molecular machines is a relatively new development, significant advances have been made. In this review an overview is given of the principal designs of artificial molecular motors and their modes of operation. Although synthetic molecular motors have also found widespread application as (multistate) switches, we focus on the control of directional movement, both at the molecular scale and at larger magnitudes. We identify some key challenges remaining in the field.
Many biomolecular motors catalyze the hydrolysis of chemical fuels, such as adenosine triphosphate, and use the energy released to direct motion through information ratchet mechanisms. Here we describe chemically-driven artificial … Many biomolecular motors catalyze the hydrolysis of chemical fuels, such as adenosine triphosphate, and use the energy released to direct motion through information ratchet mechanisms. Here we describe chemically-driven artificial rotary and linear molecular motors that operate through a fundamentally different type of mechanism. The directional rotation of [2]- and [3]catenane rotary molecular motors and the transport of substrates away from equilibrium by a linear molecular pump are induced by acid-base oscillations. The changes simultaneously switch the binding site affinities and the labilities of barriers on the track, creating an energy ratchet. The linear and rotary molecular motors are driven by aliquots of a chemical fuel, trichloroacetic acid. A single fuel pulse generates 360° unidirectional rotation of up to 87% of crown ethers in a [2]catenane rotary motor.
We formulate a nonequilibrium thermodynamic description for open chemical reaction networks (CRNs) described by a chemical master equation. The topological properties of the CRN and its conservation laws are shown … We formulate a nonequilibrium thermodynamic description for open chemical reaction networks (CRNs) described by a chemical master equation. The topological properties of the CRN and its conservation laws are shown to play a crucial role. They are used to decompose the entropy production into a potential change and two work contributions, the first due to time dependent changes in the externally controlled chemostats concentrations and the second due to flows maintained across the system by nonconservative forces. These two works jointly satisfy a Jarzynski and Crooks fluctuation theorem. In the absence of work, the potential is minimized by the dynamics as the system relaxes to equilibrium and its equilibrium value coincides with the maximum entropy principle. A generalized Landauer’s principle also holds: the minimal work needed to create a nonequilibrium state is the relative entropy of that state to its equilibrium value reached in the absence of any work.
Organisms adapt to changing environments by adjusting their development, metabolism, and behavior to improve their chances of survival and reproduction. To achieve such flexibility, organisms must be able to sense … Organisms adapt to changing environments by adjusting their development, metabolism, and behavior to improve their chances of survival and reproduction. To achieve such flexibility, organisms must be able to sense and respond to changes in external environmental conditions and their internal state. Metabolic adaptation in response to altered nutrient availability is key to maintaining energy homeostasis and sustaining developmental growth. Furthermore, environmental variables exert major influences on growth and final adult body size in animals. This developmental plasticity depends on adaptive responses to internal state and external cues that are essential for developmental processes. Genetic studies have shown that the fruit fly Drosophila, similarly to mammals, regulates its metabolism, growth, and behavior in response to the environment through several key hormones including insulin, peptides with glucagon-like function, and steroid hormones. Here we review emerging evidence showing that various environmental cues and internal conditions are sensed in different organs that, via inter-organ communication, relay information to neuroendocrine centers that control insulin and steroid signaling. This review focuses on endocrine regulation of development, metabolism, and behavior in Drosophila, highlighting recent advances in the role of the neuroendocrine system as a signaling hub that integrates environmental inputs and drives adaptive responses.
Biomolecular machines are protein complexes that convert between different forms of free energy. They are utilized in nature to accomplish many cellular tasks. As isothermal nonequilibrium stochastic objects at low … Biomolecular machines are protein complexes that convert between different forms of free energy. They are utilized in nature to accomplish many cellular tasks. As isothermal nonequilibrium stochastic objects at low Reynolds number, they face a distinct set of challenges compared to more familiar human-engineered macroscopic machines. Here we review central questions in their performance as free energy transducers, outline theoretical and modeling approaches to understand these questions, identify both physical limits on their operational characteristics and design principles for improving performance, and discuss emerging areas of research.
Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic … Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.
Cells are the basic units of all living matter which harness the flow of energy to drive the processes of life. While the biochemical networks involved in energy transduction are … Cells are the basic units of all living matter which harness the flow of energy to drive the processes of life. While the biochemical networks involved in energy transduction are well-characterized, the energetic costs and constraints for specific cellular processes remain largely unknown. In particular, what are the energy budgets of cells? What are the constraints and limits energy flows impose on cellular processes? Do cells operate near these limits, and if so how do energetic constraints impact cellular functions? Physics has provided many tools to study nonequilibrium systems and to define physical limits, but applying these tools to cell biology remains a challenge. Physical bioenergetics, which resides at the interface of nonequilibrium physics, energy metabolism, and cell biology, seeks to understand how much energy cells are using, how they partition this energy between different cellular processes, and the associated energetic constraints. Here we review recent advances and discuss open questions and challenges in physical bioenergetics.
We provide a rigorous definition of free-energy transduction and its efficiency in arbitrary -- linear or nonlinear -- open chemical reaction networks (CRNs) operating at steady state. Our method is … We provide a rigorous definition of free-energy transduction and its efficiency in arbitrary -- linear or nonlinear -- open chemical reaction networks (CRNs) operating at steady state. Our method is based on the knowledge of the stoichiometric matrix and of the chemostatted species (i.e. the species maintained at constant concentration by the environment) to identify the fundamental currents and forces contributing to the entropy production. Transduction occurs when the current of a stoichiometrically balanced process is driven against its spontaneous direction (set by its force) thanks to other processes flowing along their spontaneous direction. In these regimes, open CRNs operate as thermodynamic machines. After exemplifying these general ideas using toy models, we analyze central energy metabolism. We relate the fundamental currents to metabolic pathways and discuss the efficiency with which they are able to transduce free energy.
Macroscopic electric motors continue to have a large impact on almost every aspect of modern society. Consequently, the effort towards developing molecular motors Macroscopic electric motors continue to have a large impact on almost every aspect of modern society. Consequently, the effort towards developing molecular motors
Confronted with thermodynamically adverse output processes, free-energy transducers may shift to lower gears, thereby reducing output per unit input. This option is well known for inanimate machines such as automobiles, … Confronted with thermodynamically adverse output processes, free-energy transducers may shift to lower gears, thereby reducing output per unit input. This option is well known for inanimate machines such as automobiles, but unappreciated in biology. The present study extends existing non-equilibrium thermodynamic principles to underpin biological gear shifting and identify possible mechanisms. It shows that gear shifting differs from altering the degree of coupling and that living systems may use it to optimize their performance: microbial growth is ultimately powered by the Gibbs energy of catabolism, which is partially transformed into Gibbs energy ('output force') in the ATP that is produced. If this output force is high, the cell may turn to a catabolic pathway with a lower ATP stoichiometry. Notwithstanding the reduced stoichiometry, the ATP synthesis flux may then actually increase as compared to that in a system without gear shift, in which growth might come to a halt. A 'variomatic' gear switching strategy should be optimal, explaining why organisms avail themselves of multiple catabolic pathways, as these enable them to shift gears when the growing gets tough.
Metabolic switches are a crucial hallmark of cellular development and regeneration. In response to changes in their environment or physiological state, cells undergo coordinated metabolic switching that is necessary to … Metabolic switches are a crucial hallmark of cellular development and regeneration. In response to changes in their environment or physiological state, cells undergo coordinated metabolic switching that is necessary to execute biosynthetic demands of growth and repair. In this Review, we discuss how metabolic switches represent an evolutionarily conserved mechanism that orchestrates tissue development and regeneration, allowing cells to adapt rapidly to changing conditions during development and postnatally. We further explore the dynamic interplay between metabolism and how it is not only an output, but also a driver of cellular functions, such as cell proliferation and maturation. Finally, we underscore the epigenetic and cellular mechanisms by which metabolic switches mediate biosynthetic needs during development and regeneration, and how understanding these mechanisms is important for advancing our knowledge of tissue development and devising new strategies to promote tissue regeneration.