Global Stability of Reversible Enzymatic Metabolic Chains

We consider metabolic networks with reversible enzymatic reactions. The model is written as a system of ordinary differential equations, possibly with inputs and outputs. We prove the global stability of the equilibrium (if it exists), using techniques of monotone systems and compartmental matrices. We show that the equilibrium does not always exist. Finally, we consider a metabolic system coupled with a genetic network, and we study the dependence of the metabolic equilibrium (if it exists) with respect to concentrations of enzymes. We give some conclusions concerning the dynamical behavior of coupled genetic/metabolic systems.     

Three-species ecosystems in a solvable model

A detailed discussion of the three-species ecosystems is presented in an exactly solvable model with interactions of the Gompertz form. Three different possibilities, namely, a one-prey-two-predator system, a two-prey-one-predator system and a three-step prey-predator food chain are considered. These systems are studied not only when they include their basic prey-predator interactions, but also when various self-interactions as well as competition between like species, in different possible combinations, are included. It is then inferred, by obtaining and examining the exact solutions, as to when these systems possess stable equilibrium and when not, or when they are purely oscillatory, etc. We also study, within our model, the two-species versus three-species situation. It is seen that there are situations when the three-species system possesses stable equilibrium even under circumstances under which the corresponding two-species system is unstable. We also come across cases when the addition of the third species destroys the possibility of stable equilibrium which the initial two-species system possessed. Some other results also follow. Of particular interest is the one where the initial two-species system is purely oscillatory but the enlarged system, which is a three-step prey-predator chain, has the first and the last populations of the chain rising indefinitely and the middle population remains oscillatory. A comparison of our results with results of other authors, wherever possible, has also been made.     

Effect of hunting cooperation and fear in a predator-prey model

Dynamics of predator-prey systems under the influence of cooperative hunting among predators and the fear thus imposed on the prey population is of great importance from ecological point of view. The role of hunting cooperation and the fear effect in the predator-prey system is gaining considerable attention by the researchers recently. But the study on combined effect of hunting cooperation and fear in the predator-prey system is not yet studied. In the present paper, we investigate the impact of hunting cooperation among predators and predator induced fear in prey population by using the classical predator-prey model. We consider that predator populations cooperate during hunting. We also consider that hunting cooperation induces fear among prey, which has far richer and complex dynamics. We observe that without hunting cooperation, the unique coexistence equilibrium point is globally asymptotically stable. However, an increase in the hunting cooperation induced fear may destabilize the system and produce periodic solution via Hopf-bifurcation. The stability of the Hopf-bifurcating periodic solution is obtained by computing the Lyapunov coefficient. The limit cycles thus obtained may be supercritical or subcritical. We also observe that the system undergoes the Bogdanov-Takens bifurcation in two-parameter space. Further, we observe that the system exhibits backward bifurcation between predator-free equilibrium and coexisting equilibrium. The system also exhibits two different types of bi-stabilities due to subcritical Hopf-bifurcation (between interior equilibrium and stable limit cycle) and backward bifurcation (between predator-free and interior equilibrium points). Further, we observe strong demographic Allee phenomenon in the system. To visualize the dynamical behavior of the system, extensive numerical experiments are performed by using MATLAB and MATCONT softwares.     

The effect of evolution on host-parasitoid systems

It is well known that a simple first-order difference equation can exhibit complex population dynamics, such as sustained oscillations and chaos. An interesting problem is whether such oscillatory dynamics are expected to occur in real populations. This paper assumes that the resident system is composed of 1-host and 1-parasitoid and that only the host is allowed to evolve, but not the parasitoid. Based on the invasibility of a host to host-parasitoid systems, we investigate the dynamics of the host-parasitoid system favored by natural selection. We consider two cases. In the first case, the host's evolution involving both the intrinsic growth rate and the sensitivity to density is considered. In the second case, the host's evolution involving both the intrinsic growth rate and the vulnerability to the parasitoid is considered. In both cases, we see that the dynamics with a stable equilibrium will not be favored by natural selection without the trade-off between the host's traits which are allowed to evolve. The host-parasitoid system with a stable equilibrium will be eventually invaded by a host type that develops an unstable equilibrium with the parasitoid. If there is a trade-off between the host's traits which are allowed to evolve, a host-parasitoid system with a stable equilibrium can be favored by natural selection.     

On the relationship between dissipation and the rate of spontaneous entropy production from linear irreversible thermodynamics

When systems are far from equilibrium, the temperature, the entropy and the thermodynamic entropy production are not defined and the Gibbs entropy does not provide useful information about the physical properties of a system. Furthermore, far from equilibrium, or if the dissipative field changes in time, the spontaneous entropy production of linear irreversible thermodynamics becomes irrelevant. In 2000 we introduced a definition for the dissipation function and showed that for systems of arbitrary size, arbitrarily near or far from equilibrium, the time integral of the ensemble average of this quantity can never decrease. In the low-field limit, its ensemble average becomes equal to the spontaneous entropy production of linear irreversible thermodynamics. We discuss how these quantities are related and why one should use dissipation rather than entropy or entropy production for non-equilibrium systems.     

Linking saturation,stability and sustainability in food webs with observed equilibrium structure

Stability of a dynamic equilibrium in a predator-prey system depends both on the type of functional response and on the point of equilibrium on the response curve. Saturation effects from Holling type II responses are known to destabilise prey populations, while a type III (sigmoid) response curve has been shown to provide stability at lower levels of saturation. These effects have also been shown in multi-trophic model systems. However, stability analyses of observed equilibria in real complex ecosystems have as yet not assumed non-linear functional responses. Here, we evaluate the implications of saturation in observed balanced material-flow structures, for system stability and sustainability. We first make the effects of the non-linear functional responses on the interaction strengths in a food web transparent by expressing the elements of Jacobian ‘community’ matrices for type II and III systems as simple functions of their linear (type I) counterparts. We then determine the stability of the systems and distinguish two critical saturation levels: (1) a level where the system is just as stable as a type I system and (2) a level above which the system cannot be stable unless it is subsidised, separating a stable materially sustainable regime from an unsustainable one. We explain the stabilising and destabilising effects in terms of the feedbacks in the systems. The results shed light on the robustness of observed patterns of interaction strengths in complex food webs and suggest the implausibility of saturation playing a significant role in the equilibrium dynamics of sustainable ecosystems.     

Influence of ligand valency on ligand-influenced monomer-dimer equilibrium systems and biological control mechanisms

Two specific ligand-influenced monomer-dimer equilibrium systems are discussed. Each has a ligand-to-dimer subunit ratio of 0.5. Equilibrium characteristics of the system are described in terms of the effect of a bivalent ligand on both experimental and theoretical analysis. It is shown that Hill expressions need to be modified for multivalent ligand systems, and care is needed in equilibrium parameter determination. Some unique properties of these systems are retained, yet others are altered in the presence of a multivalent ligand. A suggestion as to the minimum amount of information needed to describe completely a ligand-influenced monomerdimer system is given. The estrogen-receptor system is presented as an attractive biological model for both theoretical and experimental study of ligand-influenced polymerizing systems and their role in cellular control requirements.     

Product-Form Stationary Distributions for Deficiency Zero Chemical Reaction Networks

We consider stochastically modeled chemical reaction systems with mass-action kinetics and prove that a product-form stationary distribution exists for each closed, irreducible subset of the state space if an analogous deterministically modeled system with mass-action kinetics admits a complex balanced equilibrium. Feinberg’s deficiency zero theorem then implies that such a distribution exists so long as the corresponding chemical network is weakly reversible and has a deficiency of zero. The main parameter of the stationary distribution for the stochastically modeled system is a complex balanced equilibrium value for the corresponding deterministically modeled system. We also generalize our main result to some non-mass-action kinetics.     

Age structure in predator-prey systems. I. A general model and a specific example

To study the effects of age structure in predator-prey systems, a general, analytically tractable model is formulated and solved. We demonstrate the usefulness of the model in a study of a specific system of two mites. We show that to maintain stable equilibrium between the herbaceous (pest) mite and the predacious mite, the nonintuitive strategy of reducing the growth rate of the predator may be necessary. The modelling technique allows a determination of the magnitude of the effect of age structure on stability.     

The stability of predator-prey systems subject to the Allee effects

In recent years, many theoreticians and experimentalists have concentrated on the processes that affect the stability of predator-prey systems. But few papers have addressed the Allee effect with focus on the their stability. In this paper, we select two classical models describing predator-prey systems and introduce the Allee effects into the dynamics of both the predator and prey populations in these models, respectively. By combining mathematical analysis with numerical simulation, we have shown that the Allee effect may be a destabilizing force in predator-prey systems: the equilibrium point of the system could be changed from stable to unstable or otherwise, the system, even when it is stable, will take much longer time to reach the stable state. We also conclude that the equilibrium of the prey population will be enlarged due to the Allee effect of the predator, but the Allee effects of the prey may decrease the equilibrium value of the predator, or that of both the predator and prey. It should also be pointed out that the impact of the Allee effects of predator and prey due to different mechanisms on different predator-prey systems could also vary.     

The effects of intraspecific competition and stabilizing selection on a polygenic trait

The equilibrium properties of an additive multilocus model of a quantitative trait under frequency- and density-dependent selection are investigated. Two opposing evolutionary forces are assumed to act: (i) stabilizing selection on the trait, which favors genotypes with an intermediate phenotype, and (ii) intraspecific competition mediated by that trait, which favors genotypes whose effect on the trait deviates most from that of the prevailing genotypes. Accordingly, fitnesses of genotypes have a frequency-independent component describing stabilizing selection and a frequency- and density-dependent component modeling competition. We study how the equilibrium structure, in particular, number, degree of polymorphism, and genetic variance of stable equilibria, is affected by the strength of frequency dependence, and what role the number of loci, the amount of recombination, and the demographic parameters play. To this end, we employ a statistical and numerical approach, complemented by analytical results, and explore how the equilibrium properties averaged over a large number of genetic systems with a given number of loci and average amount of recombination depend on the ecological and demographic parameters. We identify two parameter regions with a transitory region in between, in which the equilibrium properties of genetic systems are distinctively different. These regions depend on the strength of frequency dependence relative to pure stabilizing selection and on the demographic parameters, but not on the number of loci or the amount of recombination. We further study the shape of the fitness function observed at equilibrium and the extent to which the dynamics in this model are adaptive, and we present examples of equilibrium distributions of genotypic values under strong frequency dependence. Consequences for the maintenance of genetic variation, the detection of disruptive selection, and models of sympatric speciation are discussed.     

Duration time of a one-dimensional random walk as a function of the energies of the intermediate states: application for dissociation and relaxation processes in DNA hybrids.

Kinetic parameters of macromolecular systems are important for their function in vitro and in vivo. These parameters describe how fast the system dissociates (the characteristic dissociation time), and how fast the system reaches equilibrium (characteristic relaxation time). For many macromolecular systems, the transitions within the systems are described as a random walk through a number of states with various free energies. The rate of transition between two given states within the system is characterized by the average time which passes between starting the movement from one state, and reaching the other state. This time is referred to as the mean first-passage time between two given states. The characteristic dissociation and relaxation times of the system depend on the first-passages times between the states within the system. Here, for a one-dimensional random walk we derived an equation, which connects the mean first-passage time between two states with the free energies of the states within the system. We also derived the general equation, which is not restricted to one-dimensional systems, connecting the relaxation time of the system with the first-passage times between states. The application of these equations to DNA branch migration, DNA structural transitions and other processes is discussed.     

Enthalpy distribution functions for protein-DNA complexes: example of the binding of AT-hooks to target DNA

In this article we use the published heat capacity data of Dragan et al. [A.I. Dragan, et al., The energetics of specific binding of AT-hooks from HMGA1 to target DNA, J. Mol. Biol. 327 (2003) 393-411] on the association of proteins with DNA duplexes to construct enthalpy probability distributions for the protein/DNA complexes formed in these systems. We first analyze the multistep equilibrium that determines the species concentrations in this system to determine whether or not the DNA-peptide complex goes cleanly to DNA single-strands and peptide. Using the heat capacity data for this case we employ the maximum-entropy method to construct enthalpy probability distribution functions for the species involved in this equilibrium. We find that the distribution functions for this system clearly show bimodal behavior indicating a two-state transition from complex to non-complex form.     

Analysis of solvent-free ethyl oleate enzymatic synthesis at equilibrium conditions

This work reports experimental equilibrium data for the esterification of pure oleic acid and a fatty acid mixture with ethanol, using an immobilized Candida antarctica B lipase as catalyst. Reactions are performed in a solvent-free system, containing a mixture of substrates and different amounts of distilled water. According to the initial amount of water and the extent of the reaction, one or two liquid phases are present. Therefore, when the equilibrium is achieved, the liquid–liquid and chemical reaction equilibria have to be simultaneously satisfied.

Several reports dealing with enzymatic reactions performed in two-phase systems have found that the value of the reaction equilibrium constant calculated from overall experimental concentrations varies not only with temperature but also with substrate ratio and water content. Although this approach is a valuable way to explore equilibrium shifts in biphasic systems, it is limited to ideal systems with constant partition coefficients. The aim of this work is to consider the biphasic nature of the reactive mixture through a computational procedure that simultaneously takes into account liquid–liquid and reaction equilibria. This approach enables the determination of a classical temperature-dependent thermodynamic equilibrium constant, which accurately fits experimental equilibrium conversions over a wide range of operating conditions.     

The Critical Roles of Information and Nonequilibrium Thermodynamics in Evolution of Living Systems

Living cells are spatially bounded, low entropy systems that, although far from thermodynamic equilibrium, have persisted for billions of years. Schrödinger, Prigogine, and others explored the physical principles of living systems primarily in terms of the thermodynamics of order, energy, and entropy. This provided valuable insights, but not a comprehensive model. We propose the first principles of living systems must include: (1) Information dynamics, which permits conversion of energy to order through synthesis of specific and reproducible, structurally-ordered components; and (2) Nonequilibrium thermodynamics, which generate Darwinian forces that optimize the system. Living systems are fundamentally unstable because they exist far from thermodynamic equilibrium, but this apparently precarious state allows critical response that includes: (1) Feedback so that loss of order due to environmental perturbations generate information that initiates a corresponding response to restore baseline state. (2) Death due to a return to thermodynamic equilibrium to rapidly eliminate systems that cannot maintain order in local conditions. (3) Mitosis that rewards very successful systems, even when they attain order that is too high to be sustainable by environmental energy, by dividing so that each daughter cell has a much smaller energy requirement. Thus, nonequilibrium thermodynamics are ultimately responsible for Darwinian forces that optimize system dynamics, conferring robustness sufficient to allow continuous existence of living systems over billions of years.     

平衡与非平衡生态学在锡林河流域典型草原放牧系统中的应用

为认识放牧系统的复杂性和稳定性 ,产生了放牧系统的平衡生态学和非平衡生态学原理。放牧系统的平衡生态学原理假定 :一旦干扰在系统中发生 ,系统将偏离平衡态 ;而当干扰解除后 ,系统将自动返回原来的状态或在新的领域实现平衡。在对内蒙古锡林河流域典型草原放牧系统动态的研究中 ,来自平衡生态学的 Clem ents- Duksterhuis演替理论提供了一个基本的研究框架。尽管已经证实对退化不太严重的典型草原放牧系统 ,平衡生态学原理是适用的 ,但是对于这一地区严重退化的放牧系统的动态 ,它显然并不能给予合理的解释。事实上许多放牧系统动态遵循非平衡生态学原理。在非平衡放牧系统中 ,稳定的状态是不会实现的 ,因为在这样的系统中 ,非生物变量对于植被的动态似乎起着决定性的影响 ,从而也决定着草食动物的种群动态。状态与过渡模型基于非平衡生态学原理 ,它能够解释过度放牧下典型草原生态系统的崩溃或灌丛化 ,因此它适于该地区严重退化的典型草原放牧系统的动态。鉴于内蒙古锡林河流域典型草原放牧系统普遍严重退化的现实 ,未来该地退化放牧系统的研究应更多地应用非平衡生态学原理 ,并且严重退化的草原生态系统的恢复试验 ,特别是灌丛化草原的重建也应置于它的指导之下     

On computer simulation methods used to study models of two-component lipid bilayers

We show that the use of a computer simulation method introduced to calculate the equilibrium thermodynamic properties of a model of a two-component lipid bilayer membrane [Freire, E., & Snyder, B. (1980) Biochemistry 19, 88-94] is incorrect. This is done by comparing the method to that of Metropolis, which has been proven to generate equilibrium distribution of that model, and by showing that back-processes have been omitted in the implicit master equation of Freire and Snyder. We have illustrated this explicitly by first generating distributions according to the method of Freire and Snyder and then allowing the system to relax via the Kawasaki method, which uses the technique of Metropolis. We show that relaxation to a different distribution occurs. We also remark that the cluster distributions generated by the Freire-Snyder method are substantially different from those occurring in equilibrium distributions. Thus, conclusions about equilibrium thermodynamic properties such as specific heats and transition enthalpies or about transport properties or cluster properties at equilibrium cannot be drawn from the results obtained by using this method. Finally, we point out that the method of Freire and Snyder is appropriate to so-called aggregation models, which have been used to study irreversible growth; and we suggest biological systems that might be simulated by their method.     

Model of CE enantioseparation systems with a mixture of chiral selectors. Part I. Theory of migration and interconversion

Theory of equilibria, migration and dynamics of interconversion of a chiral analyte in electromigration enantioseparation systems involving a mixture of chiral selectors for the chiral recognition (separation) are proposed. The model assumes that each individual analyte-CS interaction is fast, fully independent on other interactions and the analyte can interact with CS in 1:1 ratio and that the analyte is present in the concentration small enough not to considerably change the concentration of free CSs. Under these presumptions, the system behaves as there was only one chiral selector with a certain overall equilibrium constant, overall mobility of analyte-selector complex (associate) and overall rate constant of interconversion in a chiral environment. We give the mathematical equations of the overall parameters. A special interest is devoted to the dynamics of interconversion. Interconversion in systems with mixture of chiral selectors is governed by two apparent rate constants of interconversion in the same way as in case of singe-selector systems. We propose the experimental design that allows to determine rates of interconversion in both chiral and achiral parts of the enantioseparation system separately. The approach is verified experimentally in the second part of the article.