A Mathematical Model for Nutrient Metabolic Chemistry

A model based on chemical kinetics for the rate of utilization and/or storage of carbohydrates, fats and proteins is derived and analyzed. This system is studied under different conditions of supply and usage and for short term dynamics and long term dynamics. Both the short term and long term models indicate that starting above an equilibrium threshold leads to growth of the stored species. Results from the short-term and long-term submodels show that the qualitative behavior depends on the levels of certain enzymes. The analysis of a model for enzyme dynamics indicates that the steady-state level of an enzyme should depend on the rate of supply of the substrate.


Introduction
When food is consumed, digestive processes break up the macroscopic chunks of meat, vegetables, bread, beverages and sweets into small organic molecules, minerals, and salts.Carbohydrates, fats and proteins are converted to simple sugars, fatty acids, and amino acids for absorption and utilization by the body for manufacturing tissue and for energy to maintain temperature, activate muscles, and run the signaling and neural systems.If the organism consumes more than it uses, the excess can be stored for use in the future.
Under normal circumstances, simple sugars, fatty acids and amino acids pass through the gastric-blood boundary and enter cells and organs where they can participate in the metabolic processes.The dominant simple sugar is glucose, which can enter the pathway to the Krebs cycle by being converted to pyruvate.
In addition, pyruvate can be converted back to glucose by a different reaction.
vert simple sugars into glycogen or pyruvate, that convert fatty acids to fat or acetyl-CoA, that convert amino acids to pyruvate or acetyl-CoA or proteins, and that convert acetyl-CoA into energy (i.e., ATP) through the Krebs cycle.This model should be regarded as somewhat more than qualitative, but not quite quantitative.First, the actual metabolic reactions are more complex than represented in Figure 1 (see Berg, et al. [8]), with several subspecies produced and further reduced in each reaction represented there.Further, most, if not all, of the reactions are catalyzed by more than one enzyme.And finally, the "reverse" reactions are actually more complicated reactions that are not simply splitting or reuniting of the species displayed.In addition, three of the species Figure 1.A schematic of the reactions.Applied Mathematics shown are actually molecules that are stored in adipose tissue, in liver and in muscle.Nonetheless, the general flow of energy-rich species to energy or storage is represented here.

Enzyme-Controlled Mass Action Model
We posit a mass-action chemistry model (see Figure 1) for the evolution of fats, carbohydrates and proteins.Let SS be the abundance of simple sugars, FA be the abundance of fatty acids, AA be the abundance of amino acids, G be the abundance of glycogen, F be the abundance of fat, and P be the abundance of proteins.The pathway to energy utilization involves two molecules, pyruvate (symbol Pyr), which is an entry molecule, and acetyl-CoA (symbol AcA), which feeds into the Krebs cycle.
These substances undergo enzymatic reactions in the body to extract energy or be processed into the organism.The reaction structure is shown in Figure 1.
Simple sugars are supplied at rate S 1 , stored at rate r 2 , released from storage at rate r 1 , and converted to pyruvate at rate r 4 .Fatty acids are supplied at rate S 2 , converted to fat at rate r 8 , released from storage at rate r 7 , converted to acetyl-CoA at rate r 6 , and synthesized from acetyl-CoA at rate r 9 .Amino acids are supplied at rate S 3 converted to proteins at rate r 12 , and body proteins are broken down into amino acids at rate r 13 .Pyruvate is converted to simple sugar at rate r 3 , and converted to acetyl-CoA at rate r 5 .Acetyl-CoA is converted to fatty acid at rate r 9 and enters the Krebs cycle at rate r 15 .Human biochemistry machinery can synthesize certain of the amino acids from components from the Krebs cycle and cannot synthesize others [8].The amino acids that can be synthesized by the body are called "non-essential", and those that cannot be manufactured by the body are called "essential", and must be supplied by diet.We further note that amino acids are also classified as "ketogenic" and "glucogenic", depending on where they enter the network in terms of energy production.The reaction structure for amino acids accounts for these two different classes of amino acids.We model this by assuming that the meta-class AA can be manufactured from pyruvate (the non-essential amino acids), and can be metabolized to either pyruvate or acetyl-CoA.We note that most of the essential amino acids are ketogenic; consequently, we do not include a reaction that converts acetyl-CoA to (ketogenic) amino acids.For the amino acids that can be converted to pyruvate, we denote this rate as r 11 .Pyruvate can be reacted to form certain amino acids at rate r 10 .The remaining amino acids can be converted to acetyl-CoA, at rate r 14 .
The fatty acid cholesterol can be converted to bile, which can enter the intestines and be reabsorbed or excreted.We model this as a rate constant for excretion of fatty acids.This is a natural, healthy excretion.We notice that proteins and sugars can also be excreted, but the amounts excreted in normal healthy individuals are small, and neglected in this model.
The metabolic reactions are discussed in entry-level biochemistry texts (see, for example, Berg, et al. [8]).These reactions occur in most of the body's cells, Glycogen phosphorylase is regulated by epinephrine, insulin and glucagon.Fatty acids are transported across the outer mitochondrial membrane by carnitine palmitoyl transferase I (CPT1), where they participate in the reactions to produce energy in the form of ATP.The step involving CPT1 is believed to be the rate-limiting step in fatty acid oxidation.CPT1 is inhibited by malonyl-CoA, which is formed by carboxylating acetyl-CoA.The conversion of acetyl-CoA to fatty acid is catalyzed by acetyl-CoA:ACP transacylase, malonyl-CoA:ACP transacylase, 3-ketoacyl-ACP synthase, 3-ketoacyl-ACP reductase, 3-hydroxacyl ACP dehydrase and enoyl-ACP reductase.
Fatty acid synthase (FAS) catalyzes fatty acid synthesis.Its main function is to catalyze the synthesis of palmitate from acetyl-CoA and malonyl-CoA, in the presence of NADPH, into long-chain saturated fatty acids.
It is clear that a model that attempts to account for all these reactions would be cumbersome.Our goal is to add understanding to what determines the reaction rates and fates for classes of nutrients.Consequently, we assume one "effective enzyme" for each reaction in Equation (2.1).

Reaction Models
We assume that all catalyzed reactions are described by a Hill's equation, that is, the reaction rate is of the form . 1 Here i χ is the abundance of species i, which is being converted by the j th reaction.j E is the abundance of enzyme j, j K is the dissociation constant and j k′ is its conversion rate.For simplicity in calculations, we write , and refer to j k as the effective conversion rate.In the context of this model, we expect that eating a meal causes spikes in amino acids, carbohydrates and fatty acids.Their levels rise far above the levels of the the enzymes that convert them to energy or to storable species.In this case, the Hill's model saturates with almost all of the enzymes bound to substrate molecules, so that the rate of conversion is constant.As the reactions proceed, the species levels fall so that not all enzymes are bound, and the conversion rate depends on abundance.
The rates of storage of fats and glycogen involve a reaction which adds the simple molecule to accumulations of the stored molecules.Thus the reaction that converts carbohydrates to glycogen is of the form and the addition of fatty acids to adipose tissue is wall.We assume the excretion of fatty acids in the intestines can be modeled as , where E k is a constant.Finally, we shall assume that the rate of conversion of acetyl-CoA to energy is given, presumably by the exercise routine.
Hence, we assume 15 r U = , the rate of energy usage.
Protein synthesis is accomplished by the mechanism whereby an amino acid, bound to its transfer RNA counterpart, is in turn bound to a peptide chain by the action of a ribosome.A rate of synthesis can be determined in terms of reaction rates for sub-processes [9] [10].There appear to be several controls in this process [8].Protein degradation is accomplished by two main pathways, proteolysis, where all cytosolic proteins are broken down in lysosomes into amino acids, and the more specific process of ubiquitination followed by degradation in a protease complex.We model this set of complex processes by 1) synthesis from amino acids in the form (2.2), and 2) a rate constant k 13 at which proteins are degraded into amino acids.

The CFP Model
With these assumptions, the equations for the evolution of nutrients and intermediaries become (2.5g) This model for the evolution of carbohydrates, fats and proteins (CFP model) derived in this section has eight dependent variables, 27 parameters, three supply rates and the rate of energy usage.Given initial values for each species and bounded smooth supply and usage rates, a unique solution will exist in the eight dimensional phase space for at least a short time, and most such solutions will persist for an infinite time.A fundamental question about such systems is whether solutions grow, decay, or oscillate in phase space.In the next two sections, we show some of the behaviors for this model (or sensible approximations thereof).In Section 3 we examine what we consider to be a model valid for short times by considering the storage and release of an energy molecule, such as occurs in the simple sugar-glycogen subsystem.In Section 4 we examine the system for behavior on a long time scale, with constant supply rates and constant energy usage both representing the long-time average supply and usage.

Analysis of Storage-Activation Model
In this section, we focus on the short-time behavior of a subsystem of the model that gives some qualitative understanding of the dynamics of Equation (2.5).We consider the competition between energy storage and usage in the short term, specifically on a scale of a several minutes to a few hours when energy molecules are used or stored.The sub-models for carbohydrates and for fats are of similar storage-activation form.We assume that the rates at which energy molecules enter the Krebs cycle are well controlled.Specifically, we assume that 0 Ex = and 4 3 r r − and 6 9 r r − are known.These latter rates are related to exercise by the individual.This assumption is further justified by the fact that ATP feeds back to control the activities of several enzymes in the carbohydrate and fatty acid metabolism reactions.

Basic Short-Term Model
Both the glucose-glycogen submodel and the fatty acid-fat submodel are of the same form, with the equations for simple sugars and fatty acids being essentially the same.We treat them as a system of equations for an energetic molecule, ( ) t , being either SS or FA, and a storage molecule, , being either G or F.
The system has the form ( ) f t is the supply rate of the energy molecule, and we denote the rate at which the stored molecule is released as k r and the rate at which the energy molecule is stored as k s .Also, K r is the dissociation constant for the release reaction, and K s is the dissociation constant for storage.
Adding the equations gives the equation for the total molecules.We shall ex- The second represents rapid ingestion and passing the gastric-blood barrier, followed by hard exercise.
This leads us to consider two feeding and exercise patterns: A: constant feeding, constant exercise; and B: concentrated feeding, concentrated exercise.
We have where ( ) ( ) Note that this implies that if intake exceeds usage, glucose and glycogen levels will increase, independent of feeding or exercise patterns.Conversely, if usage exceeds intake, both glucose and glycogen will decrease, also independent of the patterns.

Balanced Constant Feeding, Constant Exercise
Assume that ( ) 0 f t = .In this case, the system is The ε-axis is an equilibrium manifold, and there is a (non-trivial) equilibrium manifold given by Note that the equilibrium manifold will cross the ε axis if Figure 2 shows the progress of the ε-σ dynamics in the phase plane.Several trajectories are shown.Note that the equilibrium is stable for eq r s r k k K ε < and unstable for eq r s r k k K ε ≥ .The interpretation of these results is that if the system starts from an initial point that is out of equilibrium, it will progress toward equilibrium along the line ( ) ( ) the system evolves to the point ( ) 0 ,0 m , with the stored molecule decreasing while the energy molecule increases.If and the initial point lies above the equilibrium manifold, i.e., if ( ) ( ) the system evolves to the point ( ) on the equilibrium manifold given by This result is interesting from the point of view of dieting.If the diet consists of adjusting intake and exercise to balance, there is a possibility that the level of the storage molecule will increase to the stable equilibrium value.This suggests that the diet should start with exercise exceeding intake so that the initial amount of the energy molecules is lessened to a value below the equilibrium.In this case, the energy molecule level increases while the storage molecule level decreases.

Balanced Concentrated Feeding, Concentrated Exercise
Assume that . In this case, the system is identical to equations () except at the feeding and exercise times 2 p t and 2 1 p t + , respectively, where the concentration of the energy molecule jumps by f + units to the right (feeding) or f f The behavior of the model in this scenario is illustrated in Figure 3.It is quite Figure 3. Energy-Storage Evolution, Balanced concentrated intake/output.Consistent with the direction field shown in Figure 2, the two solutions that start above the equilibrium manifold proceed up and to the left, and enter orbits around the equilibrium manifold.The two solutions that start below the equilibrium manifold approach orbits on the ε axis.
similar to the behavior of the constant feeding, constant exercise model, except that the solutions show jumps in the energy molecule at feeding and exercise times, with a drift toward either no storage, or to a nearly periodic behavior around the stable equilibrium.
It is clear that the equilibrium manifold (3.10) is where short-term behavior changes from increasing storage to decreasing storage.In terms of enzyme levels, that manifold is expressed as

FP Model
In this section we analyze the enzyme-controlled mass action model in the long term.It seems plausible that reactions that process the nutrients and intermediaries occur in time periods of several minutes to a few hours, while reactions to store fat and incorporate proteins occur over several days.The basic assumption used here to reduce the equations to equations for fat and protein is that all other reactions are in equilibrium and the rates of supply and usage are constant (representing their average values).The equations could be derived by introducing a small parameter in the reaction rates.Here we simply set the appropri- ate derivatives to zero, forcing the fast reaction equations to equilibrium [11].
The result of this approximation is a model involving differential equations for fat and proteins, similar to the structure of the Alpert-Chow-Hall model, but based on enzyme chemistry.

Equilibrium Solutions
We examine the situation in which the supplies, S 1 , S 2 , S 3 , and the usage, U, are non-negative constants.In the eight-dimensional phase space, we expect an equilibrium solution, at least for some set of supplies and usage values.Note that where the equilibrium is denoted by a subscript 0. From (4.2), we see that in order for 0 FA to be non-negative, we must have 1 2

. U S S S ≤ + + (4.3)
We also see from Equation (4.1d) that the equilibrium fat level F 0 is either Note that the solution given by (4.5) predicts a negative (and therefore non-physical) value of F 0 for ( ) k K k > , two equilibrium solutions exist, while if ( ) , the only non-negative solution is given by (4.4). Figure 4 shows the two possibilities.
We also see from Equation (4.1f) that the equilibrium for the protein level is given by In order to analyze the remainder of the equilibrium values, we note that several of Equations (4.1) contain more than one term that represent reactions involving the same species.We propose to "model" these terms by combining them in a single Hill's rate model.In order to do this, we examine two such generic terms, We wish to combine these terms into one, One way to determine * , k α β and * , K α β is to minimize the L 2 error over some range of χ.We do not pursue this method since there is no natural upper limit for the concentration, and an integral from 0 to ∞ does not converge.Instead, we require the model to be correct for small χ and for very large χ.For small χ, For χ large, we have ( ) Applying the modeling assumption to Equation (4.1) gives Equations (4.13a) and (4.13c) determine the equilibrium for the Pyr-AA subsystem.For given values of S 1 and S 3 , each of these equations determines a manifold on which the equilibrium can exist.A point ( ) AA , Pyr lying on an intersection of two such manifolds determines an equilibrium.Figure 5 shows Note that not all values of S 1 and S 3 result in a steady solution.
From substitution of Equation (4.2) into Equation (4.13a), we see that a further condition for an equilibrium solution is It is clear that equilibrium solutions fail to exist for certain values of S 1 , S 2 , S 3 , and U.This is a result of the model for enzymatic reactions, specifically, Hill's equations lead to reaction rates that are bounded above, so that sufficiently high rates of supply cannot be accommodated.We summarize the situation here: If the constraint (4.3) is satisfied, Equation (4.2) gives the equilibrium solution for FA 0 , and Equations (4.4) and (4.5) give the equilibrium solutions for F 0 .If constraints (4.16) and (4.17) are satisfied, solutions for AA 0 and Pyr 0 can be found from Equations (4.13b) and (4.13c).If constraint (4.18) is satisfied, then Equation (4.13a) can be solved for AcA 0 .With these equilibrium solutions, the solution for P 0 can be found from Equation (4.6), and the solutions for SS 0 and G 0 can be found from Equations (4.1a) and (4.1b).

Stability of Equilibria
For small disturbances ( ) δ from the steady state ( ) 0 , we have Thus small perturbations to the protein equilibrium decay exponentially, so that the equilibrium solution for P is stable.
This exponential behavior in the protein subsystem also generates small disturbance behavior with the same exponent in the fat equation.In addition, the fat equation has a linearly independent solution that can be found by noting that the terms AA δ , Thus we see that the equilibrium solution Hence the equilibrium (4.31) is unstable for all 0 FA 0 > , and the solution 0 F 0 = is either 1) stable for all FA 0 , if 7 8 8 1 k k K > or 2) stable for ( ) and unstable for ( ) 1 k k K < .This analysis, coupled with the lack of critical points elsewhere in the FA-F plane implies that solutions to the differential system with initial conditions starting between the maximal equilibrium solution ( ) ( ) and 0 F 0 = approach the solution 0 F 0 = , while solutions starting above the maximal equilibrium solution must have F → ∞ .
In terms of the enzyme levels, we see that if the initial values of FA and F are such that if F max 0, 1 , FA the system predicts that F → ∞ .

Enzyme Dynamics
We note from Equations (3.10) and (4.33) that enzyme levels are involved in the In this subsection, we examine a model for the level of an enzyme involved in a Michaelis-Menten reaction [12] (see also Lin & Segel [11]).We assume that enzyme molecules are synthesized by ribosomes from mRNA templates, are activated by a process that we do not specify here, and degraded by proteolysis.We further assume that only an enzyme which is not bound to substrate is targeted by proteolysis.
In the Michaelis-Menten model, a substrate molecule A is catabolized by an enzyme E. We shall assume that the possible enzymatic reactions are as shown in Figure 6.An inactivated enzyme E is activated to E + .It is then available to bind reversibly with the substrate molecule to form the complex, E C .The complex can then release the substrate without finishing the reaction, or it can go irreversibly to completion, making a product molecule B and releasing the inactive enzyme.
We shall assume that these reactions change the amount of each of the several substances according to the following mass action reaction equations.
Here S A is the rate of supply of the energy molecule (substrate), and the rate of removal of the storage molecule (product) is S B .Also, S E is the rate at which the enzyme is manufactured by the body.In addition, we assume that inactive and activated enzyme molecules are degraded at rates 1 δ and 2 δ , respectively.We note that the rate of production of the enzyme should depend on the availability of appropriate amino acids for protein manufacture, and may be controlled by the presence or absence of other molecules.
In order to obtain the long-term behavior of this system, let us average the system over a long time interval 0 0 t t t T ≤ ≤ + .For any of the time derivative terms of the enzyme in its various forms, we have as T → ∞ .Then we have ( ) The total amount of the enzyme in its various forms is ( ) ( ) (5.9)According to this model, the steady-state level of enzyme depends on how much it is used, that is, how rapidly substrate is supplied.

Analysis
We have posited and analyzed a model for a set of metabolic reactions involving carbohydrates, fats, and proteins.The model assumes that molecules in each class are stored or metabolized into molecules pyruvate acetyl-CoA, the pathway upstream of the Krebs cycle, where the body's energy is produced.In order to gain an understanding of the behavior of this model, it is analyzed at two different time scales, a short term scale, where the molecule is stored or used, and a long term model that accounts for the storage of fats and proteins.
The CFP model that we propose and analyze here is based on the assumption that each simple sugar, fatty acid or amino acid molecule is reacted in order to be used for energy or to be stored as glycogen, fat, or protein.Each reaction is assumed to be effectively catalyzed by one enzyme.This is a simplification, since the known reactions involve a half-dozen steps, each catalyzed by a different enzyme.There is no explicit justification of this assumption, it simplifies any analysis, while retaining a reasonable reaction structure.
In Section 3, we study the behavior of a simple submodel for energy molecules and their storage.This model shows that there is a one-dimensional manifold relating the energy molecule level and the storage level that divides the behavior into two quite different types.If the level of storage molecule is above the manifold, it will evolve to even higher levels, even if intake and usage are balanced.
Below this manifold, the storage molecule level will decrease.The manifold depends on enzyme levels and the rates at which energy molecules are stored, and at which stored molecules are released.The switch from increasing to decreasing at the manifold persists when a model for feeding and exercise is considered.
The FP model analyzed in Section 4 represents the behavior of the model on long time scales, the time scales where the state of the individual is changing by gaining or losing muscle or fat mass.The model assumes that all reactions except the release-from-storage reactions are fast.This effectively assumes that the carbohydrate model is in quasi-equilibrium.For simplicity, we have assumed that the usage U does not depend on any of the nutrient levels.Energy usage is a complex issue, and as suggested by Pontzer, et al. [13], may "plateau" as a function of activity.
When this model is analyzed by finding equilibria and examining their stability, it separates into carbohydrate, protein, and fat submodels.The protein submodel predicts that the protein equilibrium is stable, while the fat submodel shows equilibrium solutions for fatty acid and fat abundance are such that if the fat level exceeds the equilibrium, the fat level increases without bound, while solutions starting in a region with fat abundance below the equilibrium proceed to a fat level of zero.This dichotomous behavior is similar to that derived from the simplified model of Section 3, and stems from the model for the enzymatic reactions for converting fatty acids to fat, and vice versa.
FP model also predicts that the equilibrium fatty acid level determines whether the system is in a mode where the abundance of fat increases or decreases.The equilibrium fatty acid level is determined by a balance between fatty acid excretion and the excess between total nutrient molecule intake and energy usage (see Equation (4.2)).The equilibrium fatty acid and fat levels then determine whether the fat/fatty acid model dynamics result in growth or decrease of the fat level.Qualitatively, a sufficiently small fatty acid level puts the fat subsystem into the state where the fat level decreases.This seems to agree with the often-quoted weight loss advice, "eat less, move more," since either of those activities would reduce the excess between total nutrient intake and usage, thereby reducing the fatty acid level.It is also interesting that the fatty acid level is a common diagnostic datum.Note that if excretions of proteins or carbohydrates were included, these long-term results would be more complex.
The results from the FP model depend on the reaction rates, i k′ , the dissocia- tion constants i K , and the enzyme levels i E .We expect that the reaction rates i k′ do not change much from place to place in the individual, nor from indi- vidual to individual.However, enzymes are in essence proteins, and the levels will not only be different in different organs and different individuals, but, as the analysis in Section 5 indicates, they could dynamically adapt to reflect how much they are used.One possible consequence of this is that the enzyme levels that have a role in determining the steady state may actually evolve.One scenario that can be imagined in the context of this approach is that instead of the fat level growing without bound, the fat level and the associated enzyme levels may evolve together to a new equilibrium.
The short term and long term qualitative behaviors depend on the levels of certain enzymes.This prompted the development of a model for enzyme production in the style of Michaelis and Menten, which assumes that enzymes which are in use are not degraded, while unused enzymes can be attacked by the protein degradation mechanisms.This model shows that enzyme levels depend on the rate of production and utilization of the substrate, with more processing leading to higher enzyme levels.
The CFP model has features that are reminiscent of the CFP model of Hall [3] [4] [5].In particular, it has equations for the rates of change of glycogen, fats and A shortcoming of the CFP model presented here is that it not consider the important part that insulin plays in the transport of simple sugars across cell boundaries, where the sugars can be used or processed for storage.Pearson, et al.
[6] [7] model this process, assuming that glucose that is released from the breakdown of glycogen is not used in fatty acid synthesis.They argue that the conversion of glycogen to glucose occurs when the system is in the fasted state which has small amounts of insulin.

Conclusions
The model of evolution of enzyme levels also suggests two side issues from this modeling approach.First, it seems quite reasonable that the mechanism that triggers hunger is able to sense the level of active enzymes that are not in complex with the metabolite (i.e., E and E + from Section 5).This would result in the metabolite being processed more quickly, and the individual's appetite returning sooner, and perhaps with stronger signaling, since there would be more non-complexed enzymes.
The second issue is that alcohol is metabolized by a sequence of enzymatic reactions.The speed of clearing alcohol from the blood depends on the enzyme level.If, indeed, the enzyme level depends on the long-time level of alcohol intake, then it seems possible that a higher long-term intake level would imply faster processing, requiring a higher consumption rate to achieve inebriation.
The model proposed and analyzed herein has several shortcomings.Perhaps the biggest is that almost all of the metabolic reactions discussed by Berg, et al. [8] are multi-step, multi-enzyme.Much of the analysis that is done here depends on the simplicity of Hill's model for enzyme reactions.It should also be noted that the analysis in Section 5 suggests that enzymes evolve on a similar time scale as the long-term metabolic dynamics.
Other shortcomings of this model are as follows: 1) Values of parameters and reaction rates are not known.2) We do not consider enzyme dynamics simultaneously with the evolution of metabolites.
3) The body has 23 amino acids, which must be incorporated into each protein in the correct order.(4) Amino acids metabolism has two separate paths for different types of amino acid, termed glucogenic and ketogenic, indicating where they enter the pathway for production of energy, as being converted to pyruvate or acetyl-CoA, respectively.A seemingly related concept is the classification of amino acids as essential or non-essential, indicating whether the body can manufacture them or not.In this model, we lump all amino acids into one species and show them as able to enter either pathway.We also show a reaction to convert pyruvate to amino acid, but no reaction to form amino acids from acetyl-CoA.See Berg, et al. [8].
In spite of the shortcomings, the model has some interesting features, some of which may suggest experiments to confirm or refute.Enzyme levels may be con- trollable by varying the copy number for crucial enzymes and may also be implemented in animal obesity models [14].Also, the presence of an unstable manifold dividing growing storage from declining storage occurs in the short term models of Section 3 and in the long-term model for fatty acid/fat metabolism results from the form of the storage/release rate models of Equations (3.1) and (4.1d), which is dependent on enzyme levels.Simultaneous observations of serum lipids (fatty acids), fat levels, and enzymes may correlate with weight gain or loss.
A mathematical model uses mathematical concepts and language to explain a phenomenon, to study the effects of different parameters, and to make predictions about its behavior.The carbohydrate-fat-protein model presented herein formulates the chemistry of metabolism as a set of differential equations and predicts a threshold for weight gain that depends on enzyme levels.To the extent that this model explains metabolism, it may offer insight into the causes of obesity.

Figure 2 .
Figure 2. Energy-storage evolution, Balanced intake/output.The curve shows the equilibrium manifold.The arrows show the direction field for the system (3.1a)-(3.1b).The solutions from numerical integration are the solid diagonal lines, with directions consistent with the direction field.

Figure 5 .
Figure 5. Intersecting manifolds show the steady solutions (points shown by circles) for certain values of the sources S 1 and S 2 .Equation (4.13b) results in a manifold indicated by a dashed curve, and Equation (4.13c) results in a manifold indicated by a solid curve.

Figure 6 .
Figure 6.Equilibrium manifolds for FA-dynamics.The shaded subregions indicate where the stability analysis implies that F 0 → .
proteins.In contrast to that model, the chemistry of the input molecules is considered in the present model.This level of detail bypasses the de novo liposynthesis (DNL) model, replacing it with a dynamical model for conversion of simple sugars to pyruvate, pyruvate to acetyl-CoA, and acetyl-CoA to fatty acids.The Chow-Hall CFP model has enjoyed significant success in understanding and interpreting the large quantity of weight and weight change data, something that the present model does not purport to do.Thus, we do not offer the present D. A. Drew DOI: 10.4236/am.2018.96045669 Applied Mathematics model as an alternative to the Chow-Hall CFP model, but consider its insights as sufficiently interesting to merit consideration.
mitochondria, and some in specific types of cells such as liver or muscle.However, no hard numbers exist for reaction rates.Moreover, most of the reactions we discuss here are catalyzed by enzymes, so that the reaction rate depends on the level of the enzyme.
D. A. Drew DOI: 10.4236/am.2018.96045651 Applied Mathematics some in the (2.1h)Most of the reactions represented in Equations (2.1) are facilitated by enzymes.Many have been analyzed in some detail, and show that several steps occur in what we write as a one-step reaction, and involve more than one enzyme.Specifically, the conversion of pyruvate to glucose is catalyzed by pyruvate kinase.The conversion of glucose to pyruvate involves several intermediates and enzymes, including hexokinase, phosphoglucose isomerase, phosphofructokinase, fructose biphosphate aldolase, triosephosphate isomerase, glyceraldehyde phosphate dehydrogenase, phosphoglycerate kinase, phosphoglycerate mutase, enolase, and pyruvate kinase.The reaction that produces glycogen from glucose is catalyzed by UDP-glucose pyrophosphorylase.Glycogen is cleaved from its chains by glycogen phosphorylase.Pyruvate is converted to acetyl-CoA by pyruvate dehydrogenase complex.Fatty acids are released from triglycerides by the action of pancreatic lipase.Triglycerides are converted from acetyl-CoA by diglyceride acyltransferase.Citrate synthase catalyzes the reaction of the acetate residue from-CoA and oxaloacetate to form carbon citrate.AMP, by acetyl-CoA and NADH.Citrate synthase is the enzyme that controls the entrance of acetyl-CoA into the Krebs cycle.The action of citrate synthase is regulated by ATP, acetyl-CoA and NADH.Applied Mathematics determination of the manifolds that separate regions where levels of a stored molecule approach zero from regions where the levels increase without bound.Consequently, enzyme levels, and what determines them, are of specific interest to the modeling of nutrition dynamics.More generally, enzymes are ubiquitous and assist most of the known digestive reactions.They add a point of control and are able to speed up reactions by enormous factors or initiate reactions that could not occur without them.The level of enzymes in reaction system regulates the velocity of the reaction.