The functional design of the rotary enzyme ATP synthase is consistent with maximum entropy production

Predicting the Rate Structure of an Evolved Metabolic Network

Background: When glucose molecules are metabolized by a biological cell, the molecules are constrained to flow along distinct reaction trajectories, which are defined by the cell's underlying metabolic network. Methods: Using the computational technique of Elementary Mode Analysis, the entire set of all possible trajectories can be enumerated, effectively allowing metabolism to be viewed in a discretized space. Results: With the resulting set of Elementary Flux Modes (EMs), macroscopic fluxes, (of both mass and energy) that cross the cell envelope can be computed by a simple, linear combination of the individual EM trajectories. The challenge in this approach is that the usage probability of each EM is unknown. But, because the analytical framework we have adopted allows metabolism to be viewed in a discrete space, we can use the mathematics of statistical thermodynamics to derive the usage probabilities when the system entropy is maximized. The resulting probabilities, which obey a Boltzmann-type distribution, predict a rate structure for the metabolic network that is in remarkable agreement with experimentally measured rates of adaptively evolved E. coli strains. Conclusions: Thus, in principle, the intracellular dynamic properties of such bacteria can be predicted, using only the knowledge of the DNA sequence, to reconstruct the metabolic reaction network, and the measurement of the specific glucose uptake rate.

Self-Organization of Enzyme-Catalyzed Reactions Studied by the Maximum Entropy Production Principle

The self-organization of open reaction systems is closely related to specific mechanisms that allow the export of internally generated entropy from systems to their environment. According to the second law of thermodynamics, systems with effective entropy export to the environment are better internally organized. Therefore, they are in thermodynamic states with low entropy. In this context, we study how self-organization in enzymatic reactions depends on their kinetic reaction mechanisms. Enzymatic reactions in an open system are considered to operate in a non-equilibrium steady state, which is achieved by satisfying the principle of maximum entropy production (MEPP). The latter is a general theoretical framework for our theoretical analysis. Detailed theoretical studies and comparisons of the linear irreversible kinetic schemes of an enzyme reaction in two and three states are performed. In both cases, in the optimal and statistically most probable thermodynamic steady state, a diffusion-limited flux is predicted by MEPP. Several thermodynamic quantities and enzymatic kinetic parameters, such as the entropy production rate, the Shannon information entropy, reaction stability, sensitivity, and specificity constants, are predicted. Our results show that the optimal enzyme performance may strongly depend on the number of reaction steps when linear reaction mechanisms are considered. Simple reaction mechanisms with a smaller number of intermediate reaction steps could be better organized internally and could allow fast and stable catalysis. These could be features of the evolutionary mechanisms of highly specialized enzymes.

Flexibility of enzymatic transitions as a hallmark of optimized enzyme steady-state kinetics and thermodynamics

We investigate the relations between the enzyme kinetic flexibility, the rate of entropy production, and the Shannon information entropy in a steady-state enzyme reaction. All these quantities are maximized with respect to enzyme rate constants. We show that the steady-state, which is characterized by the most flexible enzymatic transitions between the enzyme conformational states, coincides with the global maxima of the Shannon information entropy and the rate of entropy production. This steady-state of an enzyme is referred to as globally optimal. This theoretical approach is then used for the analysis of the kinetic and the thermodynamic performance of the enzyme triose-phosphate isomerase. The analysis reveals that there exist well-defined maxima of the kinetic flexibility, the rate of entropy production, and the Shannon information entropy with respect to any arbitrarily chosen rate constant of the enzyme and that these maxima, calculated from the measured kinetic rate constants for the triose-phosphate isomerase are lower, however of the same order of magnitude, as the maxima of the globally optimal state of the enzyme. This suggests that the triose-phosphate isomerase could be a well, but not fully evolved enzyme, as it was previously claimed. Herein presented theoretical investigations also provide clear evidence that the flexibility of enzymatic transitions between the enzyme conformational states is a requirement for the maximal Shannon information entropy and the maximal rate of entropy production.

Entropy Production and Its Application to the Coupled Nonequilibrium Processes of ATP Synthesis

Starting from the universal concept of entropy production, a large number of new results are obtained and a wealth of novel thermodynamic, kinetic, and molecular mechanistic insights are provided into the coupling of oxidation and ATP synthesis in the vital process of oxidative phosphorylation (OX PHOS). The total dissipation, Φ, in OX PHOS with succinate as respiratory substrate is quantified from measurements, and the partitioning of Φ into the elementary components of ATP synthesis, leak, slip, and other losses is evaluated for the first time. The thermodynamic efficiency, η, of the coupled process is calculated from the data on Φ and shown to agree well with linear nonequilibrium thermodynamic calculations. Equations for the P/O ratio based on total oxygen consumed and extra oxygen consumed are derived from first principles and the source of basal (state 4) mitochondrial respiration is postulated from molecular mechanistic considerations based on Nath’s two-ion theory of energy coupling within the torsional mechanism of energy transduction and ATP synthesis. The degree of coupling, q, between oxidation and ATP synthesis is determined from the experimental data and the irreversible thermodynamics analysis. The optimality of biological free energy converters is explored in considerable detail based on (i) the standard biothermodynamic approach, and (ii) a new biothermokinetic approach developed in this work, and an effective solution that is shown to arise from consideration of the molecular aspects in Nath’s theory is formulated. New experimental data in state 4 with uncouplers and redox inhibitors of OX PHOS and on respiratory control in the physiological state 3 with ADP and uncouplers are presented. These experimental observations are shown to be incompatible with Mitchell’s chemiosmotic theory. A novel scheme of coupling based on Nath’s two-ion theory of energy coupling within the torsional mechanism is proposed and shown to explain the data and also pass the test of consistency with the thermodynamics, taking us beyond the chemiosmotic theory. It is concluded that, twenty years since its first proposal, Nath’s torsional mechanism of energy transduction and ATP synthesis is now well poised to catalyze the progress of experimental and theoretical research in this interdisciplinary field.

Maximum Entropy Production Theorem for Transitions between Enzyme Functional States and Its Applications

Transitions between enzyme functional states are often connected to conformational changes involving electron or proton transport and directional movements of a group of atoms. These microscopic fluxes, resulting in entropy production, are driven by non-equilibrium concentrations of substrates and products. Maximal entropy production exists for any chosen transition, but such a maximal transitional entropy production (MTEP) requirement does not ensure an increase of total entropy production, nor an increase in catalytic performance. We examine when total entropy production increases, together with an increase in the performance of an enzyme or bioenergetic system. The applications of the MTEP theorem for transitions between functional states are described for the triosephosphate isomerase, ATP synthase, for β-lactamases, and for the photochemical cycle of bacteriorhodopsin. The rate-limiting steps can be easily identified as those which are the most efficient in dissipating free-energy gradients and in performing catalysis. The last step in the catalytic cycle is usually associated with the highest free-energy dissipation involving proton nanocurents. This recovery rate-limiting step can be optimized for higher efficiency by using corresponding MTEP requirements. We conclude that biological evolution, leading to increased optimal catalytic efficiency, also accelerated the thermodynamic evolution, the synergistic relationship we named the evolution-coupling hypothesis.

Integration of demand and supply sides in the ATP energy economics of cells

The central aspects of the energy economics of living cells revolve around the synthesis and utilization of molecules of adenosine triphosphate (ATP). Current descriptions of cell metabolism and its regulation in most textbooks of biochemistry assume that enzymes and transporters behave in the same way in isolation and in a cell. Calculations of the mechanistic or maximal P/O ratios in oxidative phosphorylation by mammalian cells generally consider only the supply side of the problem without linking to ATP-demand processes. The purpose of this article is to calculate the mechanistic P/O ratio by integration of the supply and demand sides of ATP reactions. The mechanistic stoichiometry calculated from an integrated approach is compared with that obtained from the standard model that considers only ATP supply. After accounting for leaks, slips, and other losses, the actual or operative P/O calculated by the integrated method is found to be in good agreement with the experimental values of the P/O ratio determined in mitochondria for both succinate and NADH-linked respiratory substrates. The thermodynamic consequences of these results and the biological implications are discussed. An integrated model of oxidative phosphorylation that goes beyond the chemiosmotic theory is presented, and a solution to the longstanding fundamental problem of respiratory control is found.

The maximum entropy production requirement for proton transfers enhances catalytic efficiency for β-lactamases

Movement of charges during enzyme catalytic cycle may be due to conformational changes, or to fast electron or proton transfer, or to both events. In each case, entropy production can be calculated using Terrel L. Hill's method, if relevant microscopic rate constants are known. When ranked by their evolutionary distance from putative common ancestor, three β-lactamases considered in this study show correspondingly increased catalytic constant, catalytic efficiency, and overall entropy production. The acylation and deacylation steps with concomitant proton shuttles are the most important contributors to overall entropy production. The maximal entropy production requirement for the ES↔EP or EP↔E + P step leads to optimal rate constants, performance parameters, and entropy production values, which are close to those extracted from experiments and also rank in accordance with evolutionary distances. Concurrent maximization of entropy productions for both proton transfer steps revealed that evolvability potential of different β-lactamases is similarly high. These results may have implications in particular for latent potential of β-lactamases to evolve further and in general for selection of optimized enzymes through natural or directed evolution.

Theoretical and Biological Evaluation of the Link between Low Exercise Capacity and Disease Risk

Large-scale epidemiological studies show that low exercise capacity is the highest risk factor for all-cause morbidity and mortality relative to other conditions including diabetes, hypertension, and obesity. This led us to formulate the energy transfer hypothesis (ETH): Variation in capacity for energy transfer is the central mechanistic determinant of the divide between disease and health. As a test of this hypothesis, we predicted that two-way selective breeding of genetically heterogeneous rats for low and high intrinsic treadmill running capacity (a surrogate for energy transfer) would also produce rats that differ for disease risks. The lines are termed low-capacity runners (LCRs) and high-capacity runners (HCRs) and, after 36 generations of selection, they differ by more than eightfold in running capacity. Consistent with the ETH, the LCRs score high for developing disease risks, including metabolic syndrome, neurodegeneration, cognitive impairment, fatty liver disease, susceptibility to cancer, and reduced longevity. The HCRs are resistant to the development of these disease risks. Here we synthesize ideas on nonequilibrium thermodynamics and evolution from Ilya Prigogine, Hans Krebs, and Peter Mitchell to formulate theoretic explanations for the ETH. First, at every moment in time, the atoms and molecules of organisms are reorganizing to pursue avenues for energy transfer. Second, this continuous organization is navigating in a constantly changing environment such that "strategies" are perpetually in flux and do not leave a simple footprint (evolution). Third, as a consequence, human populations demonstrate a wide variation in capacity for energy transfer that mirrors mechanistically the divide between disease and health.

Energy conservation and maximal entropy production in enzyme reactions

A procedure for maximization of the density of entropy production in a single stationary two-step enzyme reaction is developed. Under the constraints of mass conservation, fixed equilibrium constant of a reaction and fixed products of forward and backward enzyme rate constants the existence of maximum in the density of entropy production is demonstrated. In the state with maximal density of entropy production the optimal enzyme rate constants, the stationary concentrations of the substrate and the product, the stationary product yield as well as the stationary reaction flux are calculated. The test, whether these calculated values of the reaction parameters are consistent with their corresponding measured values, is performed for the enzyme Glucose Isomerase. It is found that calculated and measured rate constants agree within an order of magnitude, whereas the calculated reaction flux and the product yield differ from their corresponding measured values for less than 20 % and 5 %, respectively. This indicates that the enzyme Glucose Isomerase, considered in a non-equilibrium stationary state, as found in experiments using the continuous stirred tank reactors, possibly operates close to the state with the maximum in the density of entropy production.

Is the catalytic activity of triosephosphate isomerase fully optimized? An investigation based on maximization of entropy production

Triosephosphate isomerase (TIM) is often described as a fully evolved housekeeping enzyme with near-maximal possible reaction rate. The assumption that an enzyme is perfectly evolved has not been easy to confirm or refute. In this paper, we use maximization of entropy production within known constraints to examine this assumption by calculating steady-state cyclic flux, corresponding entropy production, and catalytic activity in a reversible four-state scheme of TIM functional states. The maximal entropy production (MaxEP) requirement for any of the first three transitions between TIM functional states leads to decreased total entropy production. Only the MaxEP requirement for the product (R-glyceraldehyde-3-phosphate) release step led to a 30% increase in enzyme activity, specificity constant kcat/KM, and overall entropy production. The product release step, due to the TIM molecular machine working in the physiological direction of glycolysis, has not been identified before as the rate-limiting step by using irreversible thermodynamics. Together with structural studies, our results open the possibility for finding amino acid substitutions leading to an increased frequency of loop six opening and product release.

Acclimation of leaves to low light produces large grana: the origin of the predominant attractive force at work

Photosynthetic membrane sacs (thylakoids) of plants form granal stacks interconnected by non-stacked thylakoids, thereby being able to fine-tune (i) photosynthesis, (ii) photoprotection and (iii) acclimation to the environment. Growth in low light leads to the formation of large grana, which sometimes contain as many as 160 thylakoids. The net surface charge of thylakoid membranes is negative, even in low-light-grown plants; so an attractive force is required to overcome the electrostatic repulsion. The theoretical van der Waals attraction is, however, at least 20-fold too small to play the role. We determined the enthalpy change, in the spontaneous stacking of previously unstacked thylakoids in the dark on addition of Mg(2+), to be zero or marginally positive (endothermic). The Gibbs free-energy change for the spontaneous process is necessarily negative, a requirement that can be met only by an increase in entropy for an endothermic process. We conclude that the dominant attractive force in thylakoid stacking is entropy-driven. Several mechanisms for increasing entropy upon stacking of thylakoid membranes in the dark, particularly in low-light plants, are discussed. In the light, which drives the chloroplast far away from equilibrium, granal stacking accelerates non-cyclic photophosphorylation, possibly enhancing the rate at which entropy is produced.

Biological evolution of replicator systems: towards a quantitative approach

The aim of this work is to study the features of a simple replicator chemical model of the relation between kinetic stability and entropy production under the action of external perturbations. We quantitatively explore the different paths leading to evolution in a toy model where two independent replicators compete for the same substrate. To do that, the same scenario described originally by Pross (J Phys Org Chem 17:312-316, 2004) is revised and new criteria to define the kinetic stability are proposed. Our results suggest that fast replicator populations are continually favored by the effects of strong stochastic environmental fluctuations capable to determine the global population, the former assumed to be the only acting evolution force. We demonstrate that the process is continually driven by strong perturbations only, and that population crashes may be useful proxies for these catastrophic environmental fluctuations. As expected, such behavior is particularly enhanced under very large scale perturbations, suggesting a likely dynamical footprint in the recovery patterns of new species after mass extinction events in the Earth's geological past. Furthermore, the hypothesis that natural selection always favors the faster processes may give theoretical support to different studies that claim the applicability of maximum principles like the Maximum Metabolic Flux (MMF) or Maximum Entropy Productions Principle (MEPP), seen as the main goal of biological evolution.

Modeling carbon allocation in trees: a search for principles

We review approaches to predicting carbon and nitrogen allocation in forest models in terms of their underlying assumptions and their resulting strengths and limitations. Empirical and allometric methods are easily developed and computationally efficient, but lack the power of evolution-based approaches to explain and predict multifaceted effects of environmental variability and climate change. In evolution-based methods, allocation is usually determined by maximization of a fitness proxy, either in a fixed environment, which we call optimal response (OR) models, or including the feedback of an individual's strategy on its environment (game-theoretical optimization, GTO). Optimal response models can predict allocation in single trees and stands when there is significant competition only for one resource. Game-theoretical optimization can be used to account for additional dimensions of competition, e.g., when strong root competition boosts root allocation at the expense of wood production. However, we demonstrate that an OR model predicts similar allocation to a GTO model under the root-competitive conditions reported in free-air carbon dioxide enrichment (FACE) experiments. The most evolutionarily realistic approach is adaptive dynamics (AD) where the allocation strategy arises from eco-evolutionary dynamics of populations instead of a fitness proxy. We also discuss emerging entropy-based approaches that offer an alternative thermodynamic perspective on allocation, in which fitness proxies are replaced by entropy or entropy production. To help develop allocation models further, the value of wide-ranging datasets, such as FLUXNET, could be greatly enhanced by ancillary measurements of driving variables, such as water and soil nitrogen availability.

Minimization of a free-energy-like potential for non-equilibrium flow systems at steady state

This study examines a new formulation of non-equilibrium thermodynamics, which gives a conditional derivation of the 'maximum entropy production' (MEP) principle for flow and/or chemical reaction systems at steady state. The analysis uses a dimensionless potential function (st) for non-equilibrium systems, analogous to the free energy concept of equilibrium thermodynamics. Spontaneous reductions in (st) arise from increases in the 'flux entropy' of the system--a measure of the variability of the fluxes--or in the local entropy production; conditionally, depending on the behaviour of the flux entropy, the formulation reduces to the MEP principle. The inferred steady state is also shown to exhibit high variability in its instantaneous fluxes and rates, consistent with the observed behaviour of turbulent fluid flow, heat convection and biological systems; one consequence is the coexistence of energy producers and consumers in ecological systems. The different paths for attaining steady state are also classified.

Bacterial chemotaxis and entropy production

Entropy production is calculated for bacterial chemotaxis in the case of a migrating band of bacteria in a capillary tube. It is found that the speed of the migrating band is a decreasing function of the starting concentration of the metabolizable attractant. The experimentally found dependence of speed on the starting concentration of galactose, glucose and oxygen is fitted with power-law functions. It is found that the corresponding exponents lie within the theoretically predicted interval. The effect of the reproduction of bacteria on band speed is considered, too. The acceleration of the band is predicted due to the reproduction rate of bacteria. The relationship between chemotaxis, the maximum entropy production principle and the formation of self-organizing structure is discussed.

Maximum entropy production and plant optimization theories

Plant ecologists have proposed a variety of optimization theories to explain the adaptive behaviour and evolution of plants from the perspective of natural selection ('survival of the fittest'). Optimization theories identify some objective function--such as shoot or canopy photosynthesis, or growth rate--which is maximized with respect to one or more plant functional traits. However, the link between these objective functions and individual plant fitness is seldom quantified and there remains some uncertainty about the most appropriate choice of objective function to use. Here, plants are viewed from an alternative thermodynamic perspective, as members of a wider class of non-equilibrium systems for which maximum entropy production (MEP) has been proposed as a common theoretical principle. I show how MEP unifies different plant optimization theories that have been proposed previously on the basis of ad hoc measures of individual fitness--the different objective functions of these theories emerge as examples of entropy production on different spatio-temporal scales. The proposed statistical explanation of MEP, that states of MEP are by far the most probable ones, suggests a new and extended paradigm for biological evolution--'survival of the likeliest'--which applies from biomacromolecules to ecosystems, not just to individuals.

Steady state of a dissipative flow-controlled system and the maximum entropy production principle

A theory to predict the steady-state position of a dissipative flow-controlled system, as defined by a control volume, is developed based on the maximum entropy principle of Jaynes, involving minimization of a generalized free-energy-like potential. The analysis provides a theoretical justification of a local, conditional form of the maximum entropy production principle, which successfully predicts the observable properties of many such systems. The analysis reveals a very different manifestation of the second law of thermodynamics in steady-state flow systems, which provides a driving force for the formation of complex systems, including life.