# E1. Introduction to Metabolic Control Analysis - Biology We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

Enzyme kinetics may seem difficult given the complicated mathematical derivations, the number of chemical species involved (an enzyme and all its substrates and products), the number of steps in the mechanism, and the large number of rate, kinetic, and dissociation constants. An example of such a “complicated” reaction explored earlier is shown below:  But single enzymes rarely act in isolation. They are components of complex pathways which have a multitude of steps, many of which are regulated. To fully understand a reaction, it is important to study the concentrations of all species in the entire pathway as a function of time. Imagine deriving the equations and determining all the relevant concentrations and constants of a pathway such as glycolysis!

To study enzyme kinetics in the lab, you have to spend much time in developing assays to measure how the concentration of species changes as a function of time to be able to measure the initial velocities of an enzyme-catalyzed reaction. However, in networks of connected metabolic reactions, the concentration of some species in the system may not change. How can this happen? Two simple examples might help explain how.

Example 1: There is zero input or output from a given reaction. This would occur in a closed system for a reversible reaction at equilibrium.

For a reversible reaction of reactant R going to product P with forward and reverse first order rate constants, the following equation can be written at equilibrium: At equilibrium, R and P don’t change.

Example 2: Consider the reaction as part of a pathway of reactions (like an open system). Now imagine a nonzero input to form reactant R and a nonzero output that consumes product P. If the input and output rates are the same, the concentrations of R and P would not change with time. That is the rate of formation of a reactant R for a given reaction is equal to the rate at which the product P of the give reaction is used. This would lead to a steady state but not equilibrium concentration of the species. Understanding a pure enzyme in vitro and in vivo requires different approaches. Biochemists like to isolate and purify to homogeneity an enzyme found in some tissue and study its mechanism of action. In doing thermodynamic measurements to measure equilibrium constants (Keq) or dissociation constants (KD), from which ΔG0 can be calculated, a protein concentration is usually held constant as the binding ligand concentration is varied (independent variable). A dependent variable signal (often spectroscopic) is measured. Measurements are made when equilibrium is reached.

For enzyme kinetic measurements in vitro, the enzyme concentration is usually held constant while substrate and modifiers are varied (independent variables) to determine how velocity (dependent variable) changes. The velocity is determined by the substrate concentration. When inhibition is studied, the substrate is varied while the inhibitor is held constant at several different fixed concentrations.

In vivo, the substrate concentration and even the enzyme concentration are determined by the velocity. Again compare this to in vitro kinetics when concentrations determine the velocity
For sets of reactions in pathways, it is better to use the term flux, J. In the steady state, the in and out fluxes for a given reaction are identical. Flux J is used to describe the rate of the system whereas rate or velocity v to used describe the rate of an individual enzyme in a system.

Computer programs can find steady states concentrations by finding the roots of the ordinary differential equations (ODE) when set to zero (vf = vr). To model a process at very low concentrations, programs can also use probabilistic or stochastic simulations to model probability distributions for species and their change with time for a finite number of particles. In such simulations, concentrations (mM) are placed with number of particles. ODEs don’t work well to describe these conditions since changes in concentrations are not continuous.

Now back to our earlier rhetorical question of deriving the equations and determining all the relevant concentrations and constants of a pathway such as glycolysis! It actually has been done by Teusink et al for glycolysis in yeast. In fact, many such complicated metabolic and signal transduction pathways have been mathematically modeled in the hopes of better understanding cellular and organismal responses. Quantitatively modeling and predicting input, outputs, and concentrations of all species in complex pathways is the basis of systems biology.

Contributors

• Prof. Henry Jakubowski (College of St. Benedict/St. John's University)

## Metabolic control analysis

Metabolic control analysis (MCA) is a mathematical framework for describing metabolic, signaling, and genetic pathways. MCA quantifies how variables, such as fluxes and species concentrations, depend on network parameters. In particular it is able to describe how network dependent properties, called control coefficients, depend on local properties called elasticities.   

MCA was originally developed to describe the control in metabolic pathways but was subsequently extended to describe signaling and genetic networks. MCA has sometimes also been referred to as Metabolic Control Theory but this terminology was rather strongly opposed by Henrik Kacser, one of the founders [ citation needed ] .

More recent work  has shown that MCA can be mapped directly on to classical control theory and are as such equivalent.

Biochemical systems theory  is a similar formalism, though with a rather different objectives. Both are evolutions of an earlier theoretical analysis by Joseph Higgins. “This book is an introduction to control in biochemical pathways. It introduces students to some of the most important concepts in modern metabolic control. It covers the basics of metabolic control analysis that helps us think about how biochemical networks operate. The book should be suitable for undergraduates in their early to mid years at college.”

Available at Amazon for \$49.95 or directly from me for only \$29.95 at books.analogmachine.org

### Content:

1, Traditional Concepts in Metabolic Regulation
2. Elasticities
3. Introduction to Biochemical Control
4. Linking the Parts to the Whole
5. Experimental Methods
6. Linear Pathways
7. Branched and Cyclic Systems
8. Negative Feedback
9. Stability
10. Stability of Negative Feedback Systems
11. Moiety Conserved Cycles
12. Moiety Conserved Cycles

Appendix A: List of Symbols and Abbreviations
Appendix B: Control Equations

A week ago I wrote about some misconceptions in metabolism that arose from a trip to a systems biology workshop in Holland. This week I discovered that Nature is perpetrating similar misconceptions in their News and Views section. I came &hellip Continue reading &rarr

I’m currently at the ISGSB meeting (International Study group for Systems Biology – formally BTK) in Holland. This is a group that values both theoretical analysis and experimentation in biochemical network analysis. Its an interesting group of people that I have literally &hellip Continue reading &rarr

## II. Analytical Procedures

In the framework of MCA, Reder  developed a generalized linear algebraic method that enables analyzing the sensitivity of metabolic systems to perturbations triggered by either a change in the internal state of the system or by the environment. The departure point of the analysis is the stoichiometric matrix, obtained from the set of differential equations of the model. The stoichiometric matrix defines the structural relationships between the processes and the intermediates participating in the metabolic network under consideration. The information in the stoichiometric matrix is independent of both the enzyme kinetics and the parameters that rule the dynamic behavior of the network. The second piece of information required to perform control analysis is the elasticity matrix defined by the dependence of each process in the metabolic network on the intermediates (e.g., ions or metabolites) included in the model. The elasticity matrix is quantified through the derivatives of the rates of individual processes with respect to each possible effector

By applying matrix algebra, the corresponding control and response coefficients matrices are obtained. The regulation and control in the network is quantified by both kinds of matrices. The regulation exerted by internal or external effectors to a network can be quantified by the response coefficient .

The following matrix relationships were used in the computation of flux and metabolite concentration control coefficients

with C and Γ referring to the flux- and metabolite concentration control coefficients, respectively Idr, the identity matrix of dimension r, or the number of processes in the network under study. Dxv the elasticity matrix Nr the reduced stoichiometric matrix and L, the link matrix that relates the reduced- to the full-stoichiometric matrix of the system (for further details see ).

## INTRODUCTION TO METABOLIC ENGINEERING Chapter 1 of textbook - PowerPoint PPT Presentation

### CE508 LECTURE ONE INTRODUCTION TO METABOLIC ENGINEERING Chapter 1 of textbook CE508 Metabolic Engineering Instructor Mattheos Koffas Course Information . &ndash PowerPoint PPT presentation

PowerShow.com is a leading presentation/slideshow sharing website. Whether your application is business, how-to, education, medicine, school, church, sales, marketing, online training or just for fun, PowerShow.com is a great resource. And, best of all, most of its cool features are free and easy to use.

## Comparing methods for metabolic network analysis and an application to metabolic engineering

Bioinformatics tools have facilitated the reconstruction and analysis of cellular metabolism of various organisms based on information encoded in their genomes. Characterization of cellular metabolism is useful to understand the phenotypic capabilities of these organisms. It has been done quantitatively through the analysis of pathway operations. There are several in silico approaches for analyzing metabolic networks, including structural and stoichiometric analysis, metabolic flux analysis, metabolic control analysis, and several kinetic modeling based analyses. They can serve as a virtual laboratory to give insights into basic principles of cellular functions. This article summarizes the progress and advances in software and algorithm development for metabolic network analysis, along with their applications relevant to cellular physiology, and metabolic engineering with an emphasis on microbial strain optimization. Moreover, it provides a detailed comparative analysis of existing approaches under different categories.

## Metabolic Pathway Analysis as Part of Systems Biology - PowerPoint PPT Presentation

### Friedrich Schiller University Jena, Germany. Friedrich Schiller (1759-1805) . Systems theoretical approaches used in this field for a long time, for example: . &ndash PowerPoint PPT presentation

PowerShow.com is a leading presentation/slideshow sharing website. Whether your application is business, how-to, education, medicine, school, church, sales, marketing, online training or just for fun, PowerShow.com is a great resource. And, best of all, most of its cool features are free and easy to use.

## Introduction: Kleiber and metabolic scaling

In 1932, Kleiber published a paper in an obscure journal  showing that standard metabolic rates among mammals varied with the three-quarters power of body mass: the so-called "elephant to mouse curve", termed "Kleiber's law" in this review. Since that date, this and similar allometric scaling phenomena have been widely and often intensively investigated. These investigations have generated continuing debates. At least three broad issues remain contentious, each compounded on the one hand by the problem of obtaining valid data (in particular, finding procedures by which reliable and reproducible measures of standard metabolic rate can be obtained, especially in poikilotherms) and on the other by statistical considerations (in particular, the validity of fitting scattered points to a straight line on a semi-logarithmic plot).

The first issue is disagreement as to whether any consistent relationship obtains between standard metabolic rate and body mass. Moreover, those who acknowledge such a relationship hold divergent opinions about its range of application. Is it valid only for limited numbers of taxa, or is it universal? Since the 1960s there has been a measure of consensus: a consistent allometric scaling relationship does exist, at least among homoiotherms. Nevertheless, not all biologists agree, and scepticism is widespread, particularly about the alleged universality of Kleiber's law.

Second, assuming that some version of Kleiber's law (a consistent metabolic scaling relationship) applies to at least some taxa, there are disagreements about the gradient of the semi-log plot. That is, if B = aM b , where B = standard metabolic rate, M = body mass, and a and b are constants, what is the value of b? Kleiber  and many subsequent investigators claimed that b = 0.75, and on this matter too a measure of consensus has obtained since the 1960s. Once again, however, not all biologists agree. A significant minority of investigators hold that b = 0.67 and other values have been suggested, at least for some organisms.

Third, assuming a consistent scaling relationship and an agreed value of b, how is Kleiber's law to be interpreted mechanistically? What is its physical or biological basis? For those who claim that b = 0.67, this issue is simple: standard metabolic rate depends on the organism's surface to volume ratio. But for proponents of the majority view, that b = 0.75, the issue is not simple at all. Many interpretations have been proposed, and since several of these are of recent coinage and seem to be mutually incompatible, a critical comparative review seems timely.

Kleiber's initial paper  found support within a decade. The allometric scaling relationship B = aM b (B = standard metabolic rate, M = body mass, a and b are constants and b is taken to be approximately 0.75), was inferred by other investigators during the 1930s [2, 3]. Relevant data have been reviewed periodically since then (e.g. [4–15]) and recent developments have rekindled interest in the field.

Many biological variables other than standard metabolic rate also reportedly fit quarter-power scalings (relationships of the kind V = kM b , where V is the variable in question, k is a constant and b = n/4 n = 3 for metabolic rate). Examples include lifespans, growth rates, densities of trees in forests, and numbers of species in ecosystems (see e.g. ). Some commentators infer that Kleiber's law is, or points to, a universal biological principle, which they have sought to uncover. Others doubt this, not least because it is unclear how (for example) tree densities can be consequences of metabolic scaling or can have the same mechanistic basis. This article focuses on the metabolic rate literature, mentioning other variables only in passing, because most debates in the field have arisen from metabolic rate measurements.

### Variations in the value of b

Most debates about the value of b assume some version of Kleiber's law: i.e. that a single allometric scaling relationship fits metabolic rates over a wide range of organisms. However, as noted in the introduction, there are dissenters. Everyone acknowledges considerable variation both within and among taxa, no matter whether b = 0.75, 0.67, or some other number. The question is whether these variations are deviations from a general law, or whether there is no such law. Conflicting opinions on this fundamental point recall the traditional philosophical difference between physicists and biologists: the former are inclined to see abstract mathematical generalities in any set of numerical data, the latter to see concrete particulars. All recent attempts to explain Kleiber's law by "universal" models have involved physicists and mathematicians the sceptics are predominantly biologists.

Dodds et al.  re-examined published scaling data from Kleiber's original paper onwards and concluded that the consensus (b = 0.75) was not statistically supported. Feldman  found no evidence for any wide-ranging allometric power law in biology and dismissed all attempts to explain scaling relationships by physical or mathematical principles. Atanasov and Dimitrov  found evidence that b ranges from around 0.67 to more than 0.9 over all major animal groups, the values perhaps reflecting complexity of organisation single values such as 0.75 emerge only as averages over each group. Other investigators have been less sceptical publications by Enquist and Niklas [19, 20] give particularly impressive support to the generality of Kleiber's law because Niklas was previously among the doubters.

Whatever one's position, it is indisputable that the Kleiber relationship has many exceptions, even among mammals. Bartels  showed that some mammals, such as shrews, have B values well above those expected from the Kleiber curve. Andersen  discussed the high B values for whales and seals and attributed them to the cold-water habitat. Nevertheless, Kleiber's law has been extended beyond placental mammals to birds and marsupials. Birds have generally higher a values than placental mammals and marsupials have lower ones, but the 0.75-power relationship is still inferred by many investigators (e.g. ). McNab  accepted Kleiber's law as a general approximation but emphasized species variations, which he attributed to differences in diet, habitat and physiological adaptation. Elgar and Harvey  also found variability among groups of species but reasoned that standard metabolic rates vary taxonomically rather than with temperature regulation, food intake or activity. Economos  was also critical of McNab, at least in respect of mammals.

It is difficult to define "standard metabolic rate" in poikilotherms ambient temperature, time since last meal and other variables markedly affect measurements [9, 13, 25]. A heterogeneous array of poikilotherm data  revealed an "average" b value of roughly 0.75. There were wide divergences in some taxa notwithstanding these, Hemmingsen [4, 5] argued that over all animals, plants and protists, metabolic rate scales as the 0.75-power of body mass. More recently published data [26, 27] support this conclusion for a wide range of organisms and body masses. However, a careful re-evaluation of Hemmingsen's data by Prothero  cast further doubt on the applicability of Kleiber's law to unicellular organisms. Scepticism persists, mostly on the grounds of the intrinsic variability of the data, which is too often underestimated because it is disguised in the customary logarithmic plots and is seldom subjected to adequate statistical analysis [11, 29]. However, this too has been debated a suitable choice of procedures for estimating parameters might eliminate inconsistencies and discrepancies from the data, giving more credence to the belief that b = 0.75 over a wide range of taxa . In the following section we shall examine some of the more divergent data in more detail.

In short, there is a clear but by no means total consensus that (i) Kleiber's law is widely (even universally) applicable in biology, (ii) b is approximately 0.75. Variability in the data is generally admitted, so the consensus – and the claim that Kleiber's law manifests a general biological principle – can legitimately be doubted.

### The mass transfer model 

Some of the doubts about the consensus are powerfully supported by studies on small aquatic organisms. Reviewing a large literature on metabolic rates in aquatic invertebrates and algae, Patterson  deployed chemical engineering principles to explain why the b values ranged from about 0.3 to 1.2 in these taxa (his Table 1 provides an excellent summary). Assuming that the delivery of nutrients to each organism entails diffusion through a boundary layer, Patterson showed how water movements and organism size might affect such delivery and hence determine metabolic rate. Using simple geometrical models of organisms (plates, cylinders and spheres), he derived b values ranging from 0.31 to 1.25, more or less consistent with the experimental values.

Patterson plotted two dimensionless numbers against each other, viz. Sherwood number, Sh = hmW/D, where hm = mass transfer coefficient, W = characteristic dimension of organism and D = diffusivity and Reynolds number (a function of organism size), Re = ρUW/μ, where ρ = density, U = water flow speed and μ = coefficient of viscosity. The graphs, which had the form Sh = c.Re d , where d = 0.5 for ideal laminar flow and 0.8 for turbulent flow (c is a constant of proportionality), revealed the relative importance of diffusion and mass transfer (convective movement) in the supply of materials. Patterson was able to derive an expression for hm, and was thus able to relate the supply of materials to body mass.

The two main attractions of this model are (1) good agreement with a wide range of data and (2) derivation from basic physical principles without ad hoc biological or other assumptions. Patterson's approach has implicit support in the literature: Coulson  used chemical engineering principles to argue that mammalian metabolic rates are supply-limited, but he did not develop the argument in mathematical detail. However, Patterson's model has drawbacks. First, it is hard to see how his reasoning can be generalised to other taxa, notwithstanding Coulson's proposal (discussed in a later section). Second, by focusing on diffusion and convective mass transfer, he ignored active processes in the uptake of materials, which are likely to dominate in many organisms. Third, he assumed that metabolism in general is supply-limited in homoiotherms at least, it is more nearly demand-limited under resting conditions, though even this is an oversimplification.

The Patterson model has not been given much attention by other investigators in the field and perhaps it deserves more consideration. Despite its inherent limitations (it is exclusively concerned with small aquatic eukaryotes) it is a potentially fruitful contribution to biophysics.

### Scaling of metabolic rate with surface-to-mass ratio

Several workers accept the reality of allometric scaling but question the value b = 0.75, which a consensus of physiologists has accepted since the 1960s. Many of these sceptics claim that the "true" value of b is 0.66 or 0.67 because the principal determinant of metabolic scaling is the surface-to-volume ratio of the organism hence, assuming constant body density, the surface-to-mass ratio. The first study to suggest this explanation for the mass dependence of B is attributed to Rubner , who studied metabolic rates in various breeds of dog. Heusner  reported that b is approximately 0.67 for any single mammalian species and suggested that the interspecies value of 0.75 is a statistical artefact. Feldman and McMahon  disagreed, but Heusner sustained his position in subsequent articles. For instance, reviewing a substantial body of published data , he argued that metabolic rate data for small and large mammals lie on parallel regression lines, each with a gradient of approximately 0.67 but with different intercepts (i.e. values of a, termed the "specific mass coefficients"). Hayssen and Lacy  found b = 0.65 for small mammals and b = 0.86 for large ones, again suggesting that b = 0.75 is a cross-species "average" with no biological significance but it is questionable whether their data were measurements of standard metabolic rate in all cases. McNab  reported lower values: 0.60 and 0.75, respectively. Heusner  reasoned that if a few large mammals are added to a sample of predominantly small ones, a single regression line for all the data might have a gradient around 0.75. This, however, is misleading, as the following paragraphs will argue.

According to Heusner, the ratio B/M 0.67 is a mass-independent measure of standard metabolism. Variations indicate the effects of factors other than body mass. Other workers broadly share Heusner's opinion (see e.g.  for review and  for a good recent exemplar). Bartels  found a value of 0.66 for mammals Bennett and Harvey  reported 0.67 for birds. Of course, if B varies as M 0.67 , the interesting problem is not the index (b) in the Kleiber equation but the allegedly constant relationship between specific mass coefficient (a) and body size. This point was developed by Wieser , who distinguished the ontogeny of metabolism, which comprises several phases but follows the surface rule (M 0.67 ) overall, from the phylogeny of metabolism, which concerns the mass coefficients (a). Following Heusner's argument, Wieser  wrote the allometric power law in the form B = anM 0.66 and deduced that the specific mass coefficient an = aM 0.09 . Here, a is an interspecific mass coefficient (3.34 w in mammals if M is in kg). Another difficulty with this type of explanation lies in the calculation of body surface area the Meeh coefficient, k, where surface area = kM 0.67 , is difficult to measure unequivocally but is generally taken as

10 (see ). Yet another possible difficulty was identified by Butler et al. , who questioned Heusner's dimensional analysis argument and concluded that no version of Kleiber's law (i.e. no value of b that is constant over a range of species) could be substantiated by his approach.

The claim that b = 0.67 remains a minority view. Those who accept it are faced with the twin difficulties of (i) establishing that their estimates of surface area are correct and (ii) explaining why, in Wieser's notation, an = aM 0.09 . Moreover, even if such arguments as Heusner's are valid for homoiotherms, it is hard to justify their extrapolation to poikilothermic animals, plants and unicellular organisms, all of which are held by consensus to fit Kleiber's law (but see the two preceding sections). Why should temperature fluxes across the body surface be the main determinants of metabolic rate in poikilotherms, particularly microorganisms? Even in mammals, maintenance of body temperature might not be the main contributor to energy turnover at rest (see later). Contrary to the view of Dodds et al. , therefore, b = 0.67 cannot be treated as a "null hypothesis".

Throughout the remainder of this article, the consensus position will be assumed: Kleiber's law is valid for a wide range of organisms, and b = 0.75. This assumption is made tacitly and provisionally and does not imply dismissal of the foregoing sceptical arguments but a field can only be reviewed coherently from the consensus point of view.

### McMahon's model 

A vertical column displaced by a sufficiently large lateral force buckles elastically. The critical length of column, lcr, = k(E/ρ) 1/3 d 2/3 , where d = column diameter, E = Young's modulus and ρ = density. If E and ρ are constant then lcr 3 = cd 2 , where c is a constant of proportionality. McMahon  applied this reasoning to bone dimensions for stationary quadrupeds. In a running quadruped the limbs support bending rather than buckling loads but the vertebral column receives an end thrust that generates a buckling load. It follows that all bone proportions change in the same way with animal size. The mass of a limb, wl, = αld 2 , where α is a constant. If wl is proportional to M, as it generally must be, then M = βld 2 , where β is another constant. Hence (given the above relationship between l and d) M is proportional to l 4 , implying that l is proportional to wl 1/4 hence d is proportional to wl 3/8 , or M 3/8 . Empirical support for this relationship appeared in .

McMahon  also applied this argument to muscles. The work done by a contracting muscle, W, is proportional to σAΔl, where σ is tensile strength, A is the cross-sectional area and Δl is the length change during contraction. The power developed, W/t (t = time), is therefore σAΔl/Δt. Since σ and Δl/Δt are roughly constant and independent of species, W/t varies with A and since A is proportional to d 2 , W/t it is proportional to d 2 , and therefore to (M 3/8 ) 2 = M 3/4 . If this deduction applies to any skeletal muscle (as seems plausible), then it applies to the entire set of metabolic variables supplying the muscular system with nutrients and oxygen. Hence, B varies as M 3/4 . A broadly comparable but simpler argument was advanced by Nevill  large mammals have proportionately more muscle mass than smaller ones. If the contribution of the muscle to B (which Nevill assumes is proportional to M) is partialled out, then the residual B is proportional to M 2/3 . Nevill's paper is seldom cited.

One difficulty with McMahon's model is that little of the energy turnover under conditions of standard metabolic rate measurement entails muscle contraction. The model might still be valid if maximum metabolic rate followed the same allometric scaling law as B this has been widely believed, and Taylor et al.  adduced evidence for it. However, recent detailed studies [46–48] indicate that maximum metabolic rate in birds and mammals scales as M 0.88 , not M 0.75 , although there are disagreements about whether aerobic capacity determines the allometry of maximum metabolic rate [48, 49]. Weibel  presented a large set of data to this effect. (On the other hand, there are reports that in birds the index decreases rather than increases with increasing metabolic output, e.g. .) Another drawback of the McMahon model is that it cannot apply to organisms without muscles, such as protists. This perhaps explains why McMahon's elegant deduction has been largely ignored in recent debates about Kleiber's law.

### The Economos model 

An increased gravitational field increases energy metabolism in animals [52, 53]. Work against gravity is proportional to M 1.0 . If maintenance metabolism were related to surface area (proportional to M 0.67 ) then a combination of the two effects, surface-to-mass ratio and work against gravity, might explain the observed M 3/4 relationship. This model  is difficult to assess: it is not clear why the two proposed factors, surface area dependence and gravitational loading, should combine for all animals (and other taxa) in just the right proportions to generate a 0.75-power dependence on body mass. To take just one example, aquatic microbes are more affected by Brownian motion than by gravity, so why should they show the same balance between surface-to-mass ratio and gravitational effects as mice or elephants? Pace et al.  suggested that the Economos model could be critically tested under conditions of weightlessness in space. No corroboration (or refutation) by studies on astronauts has been reported.

### Allometric scaling in cells and tissues

Before more recent models purporting to explain Kleiber's law are discussed, some comments are needed on scaling of metabolism at the organ, tissue and cell levels. Belief that the Kleiber relationship can be explained in terms of the inherent properties of the cells dates from the 1930s [3, 55] and persists (e.g. [56, 57].

Standard metabolic rate (B) is usually measured as oxygen consumption rate, which correlates with nutrient utilization [9, 15] and rates of excretion of nitrogenous and other wastes  so research in the field has been dominated by respiratory studies. Lung volume, trachaeal volume, vital capacity and tidal volume all scale as M but respiratory frequency varies as M -0.31 , ventilation rate as M 0.77 and oxygen consumption rate as M 0.72 [58–60]. All mammals extract a similar percentage of oxygen (

3%) from respired air . The significance of "pulmonary diffusion capacity" has been debated it scales as M 1.0 so it is disproportionate in bigger animals [17, 61–65].

Stahl  described the scaling of cardiovascular and haematological data. Blood haemoglobin concentration is the same for all mammals except those adapted to high altitudes. Blood volume is

6–7% of body volume for all mammals except aquatic ones. Erythrocyte volume varies with species but bears no obvious relationship to M. The oxygen affinity of haemoglobin varies with body size, being lower in smaller mammals, which unload oxygen to their tissues more rapidly. Capillary density is more or less constant in mammals with bodies larger than a rat's, though it is greater in the smallest mammals . The heart accounts for

0.6% of body mass in all mammals . Heart rate scales as M -0.25 , cardiac output as M 0.81 (60) and circulation time as M 0.25 . The energy cost of supplying the body with 1 ml of oxygen is similar for all mammals .

Standard metabolic rate has two main components: service functions, e.g. the operation of heart and lungs and cellular maintenance functions, e.g. protein and nucleic acid turnover (e.g. ). Krebs  elucidated this second component by studying tissue slices his investigation has since been extended. Oxygen consumption per kg decreases with increasing M in all tissues, but tissues do not all scale identically. Horse brain and kidney have half the oxygen consumption rates of mouse brain and kidney but the difference between these species in respect of liver, lung and spleen is 4-fold [68–70]. Metabolic rate in liver scales as M 0.63 for some organs the exponent is closer to 1.0 the sum of oxygen consumption rates over all tissues gives – approximately – the expected 0.75 index . The difficulty of recalculating B from tissue-slice data is considerable, so the Martin and Fuhrman calculation  has wide confidence limits. Spaargen  suggested that tissues that use little oxygen constitute different percentages of body mass in large and small mammals, leading to a distortion of the surface law (B = M 2/3 ), which would otherwise be valid. More recently, however, Wang et al.  repeated the Martin and Fuhrman calculation using improved data, and found impressive support for the consensus B = M 3/4 .

Cells of any one histological type are size-invariant among mammals but allometric scaling is reported at the cellular level e.g. the metabolic rate of isolated hepatocytes scales as M -0.18 . Numbers of mitochondria per gram of liver (or per hepatocyte), however, scale as M -0.1 [75, 76]. The apparent discrepancy between these values might be illusory (cf. ), or it might indicate a greater proton leak in mitochondria from livers of smaller animals  or allometry in redox slip . Also, larger animals have smaller inner mitochondrial membrane surface areas (the scaling is M -0.1 ) and different fatty acid compositions . The discrepancy between the scalings of hepatocyte and whole-body metabolism is probably explained by the decrease in liver mass, which scales as M 0.82 [75, 80]. Combining liver mass with hepatocyte oxygen consumption, the derived scaling for liver metabolism is M 0.82 .M -0.18 = M 0.64 , consistent with the experimental tissue-slice data (M 0.63 see above). Combining liver mass with mitochondrial number per hepatocyte gives a similar value . Cytochrome c and cytochrome oxidase contents scale roughly as M 0.75 [81–85]. The allometric scaling of mitochondrial inner membrane area, and the body-size-related differences in unsaturated fatty acid content, remain unexplained.

Isolated mammalian cells reportedly attain the same mitochondrial numbers and activities after several generations in culture, irrespective of the tissue of origin or the organism's body mass [86–88]. If allometric scaling is lost at the cellular level after several generations in vitro, then presumably mitochondrial densities, inner membrane areas and cytochrome levels somehow become "normalized". This is a readily testable prediction [see ], but it does not appear to have been subjected to critical experiments. If it is corroborated there will be interesting mechanisms to investigate.

The main conclusions from this section are: (a) different organs make different contributions to the scaling of whole-organism metabolic rates (b) differences at the cellular level make relatively small contributions to scaling at the organ level (c) these differences at cellular level might disappear altogether after several generations in culture. The most striking conclusion is (b). It implies that allometric scaling of metabolic rate does not after all, for the most part, reside in cellular function but at higher levels of physiological organisation. If this is the case, then the alleged applicability of Kleiber's law to unicellular organisms is called into question.

### Resource-flow models

Coulson's flow model  was mentioned earlier. It relates tissue or organ oxygen consumption rates to circulation times, i.e. to the rate of supply of oxygen and nutrients, and these scale as M 0.25 (see previous section). Coulson's approach contrasts with traditional biochemical measurements: the principal variable is not the concentration of a resource but the supply rate metabolic activity depends on encounter frequency not concentration. This perspective merits further development, particularly by extension to the cell internum [89–93]. Obviously, it is within the cell that the reactant molecules are passed over the catalysts and the flow rate increases with the cell's metabolic activity, as Hochachka  cogently described.

However, flow theories advanced to explain Kleiber's law have not followed this line of argument. Banavar et al. [94, 95] and Dreyer and co-workers [27, 96] have shown that the Kleiber relationship can be deduced from the geometries of transport networks, without reference to fluid dynamics. Broadly, these authors argue that as a supply network with local connectivity branches from a single source (in a mammalian circulatory system, the heart is the source), the number of sites supplied by the network increases. Natural selection has optimized the efficiency of supply. A general relationship can be derived between body size and flow rate in the network: delivery rates per unit mass of tissue vary with the quarter-power of body size (M), implying the validity of Kleiber's law.

The most detailed account of this argument  begins with the reasonable assumption that M scales with L D , where L is the physical length of the organism and D is its dimensionality. It proceeds with a theorem: the sum of flows through all parts of the network, F, is proportional to the (dimensionless) length multiplied by the metabolic rate. A quantity measuring the total flow of metabolites per unit mass of organism is then defined: r1 = F/M. r1 (which has units of inverse time) measures the dependence of the network's geometry on body mass, so it indicates the energy cost of metabolite delivery. Another parameter, r2, measures the metabolite demand by the tissues: r2 = the dimensionless length of the "service volume" (the amount of tissue that consumes one unit of metabolite per unit time). It is then deduced that B is proportional to (Mr1/r2) D/(D+1) . Provided that r1 and r2 change proportionately – i.e. supply always matches demand – then for a three-dimensional organism, Kleiber's law follows. According to Banavar et al. , deviations from Kleiber's law indicate inefficiency or some physiological compensation process.

This model has been criticized  because the assumed network does not resemble (e.g.) the mammalian circulatory system, where only terminal nodes (capillaries), not all nodes (as the model implies), are metabolite exchange sites. Also, the model seems to predict that r1/r2 will decrease as B rises from standard to maximal but the best data suggest the opposite trend (see earlier discussion: [46–48]). Banavar et al. do not explicitly allow for differences among tissue types, which are considerable (see above), except perhaps in terms of rather implausible variations among r1/r2 ratios. On the other hand, the model is simple and flexible and it reflects recent developments in the physics of networks. If it could be applied to flow at the cellular level, it might accord with the requirements discussed at the beginning of this section though it is difficult to see how this can be achieved.

Rau  also advanced a fluid-flow model, but his conception is physical not geometrical. Assuming Pouseille flow through an array of similar tubes, such as capillaries, and a roughly constant flow speed, Rau used scaling arguments to derive the relationship t = kM 1/4 , where t is the transport time and k is a constant. If the fluid transport rate (essentially the reciprocal of t) is proportional to B/M, Kleiber's law follows. However, Rau's model appears to assume that because metabolic rate is energy per unit time, it can be equated with the product of fluid volume flow rate and pressure (since energy is equal to pressure times volume). This assumption, which appears to be based exclusively on dimensional analysis, is fallacious.

### Four-dimensional models

Blum  observed that the "volume" of an n-dimensional sphere of radius r is V = π n/2 r n /Γ(n/2 + 1), and that A = dV/dr = nπ n/2 r n-1 /Γ(n/2 + 1). Here, Γ(n) is the gamma-function such that Γ(n + 1) = nn, Γ(2) = 1 and Γ(3/2) = π 1/2 /2. Suppose two objects have "volumes" V1 and V2 and "areas" A1 and A2. From the foregoing, A1/A2 = (V1/V2) (n-1)/n so if n = 4, a 3/4-power relationship between "volumes" (hence, masses?) emerges from a familiar mathematical principle. Might Kleiber's law therefore follow from a four-dimensional description of organisms? Speakman  pointed out that if n = 4, then A is volume (it has three dimensions) and V is hypervolume, the biological significance of which is obscure. However, West et al. [88, 101, 102] have indeed proposed a four-dimensional model to explain the Kleiber relationship, and considerable claims have been made for their account.

This model addresses the supply of materials (particularly oxygen) through space-filling fractal networks of branching tubes. It assumes that as a result of natural selection, organisms maximize their use of resources. The initial account  assumed that energy dissipation is minimised at all branch-points in the network and that the terminal branches are size-invariant (for instance, blood capillaries are the same lengths and diameters in mice and elephants). Kleiber's law and analogous scalings were deduced from these assumptions. In particular, the three-quarters-power exponent was shown to be inherent in the geometry of a branching network that preserves total cross-sectional area at each branch point. The circulatory systems of large animals such as mammals are not exactly area-preserving, but West et al.  reasoned that this objection could be circumvented by considering the pulsatile flow generated in the larger arteries by the action of the heart.

A second, simpler account  developed the model from a geometrical basis. The crucial feature of the branching network is the size-invariance of the terminal units. The effective exchange area, a, is a function of the element lengths at each level of the hierarchy, but one of these, the terminal one (l0), is invariant. Writing Φ as a dimensionless function of the (dimensionless) ratio l1/l2 leads to

Introducing a scaling factor, λ, leads to

which is not proportional to λ 2 because l0 is fixed. The dependence of Φ on λ is not known a priori, but it can be parameterized as Φ(l0/λ l1, l2/l1. ) = λ ε Φ(l0/l1, l2/l1. ), where ε is between 0 and 1. This power law reflects the fractal character of the network's hierarchical organization. Similar reasoning is applied to body volume, hence body mass, and the following expression for the exchange surface area is derived:-

a = kM r , r = (2 + ε)/(3 + ε + ζ),

where k is a constant and ζ (0 < ζ < 1) is an arbitrary exponent of length, just as ε is an arbitrary exponent of area. If natural selection has acted to maximize the scaling of a, then ε must tend to 1 and ζ to 0. This gives r = 0.75. If a limits the supply of oxygen and nutrients, and hence determines standard metabolic rate, then B is proportional to a and Kleiber's law follows.

The model has several attractions: it derives from well-established physical principles, invokes natural selection and is mathematically impeccable. It implies that cells and organelles transport materials internally along space-filling fractal networks rather than by "diffusion", which seems correct [83, 85, 86, 103]. The self-similarity of these transport networks is emphasized particularly in . The dimensionalities of effective exchange surfaces, a, are predicted to be closer to 3 than 2 empirically, the microscopic convolutions of surfaces such as the mammalian intestinal mucosa are well known. The mass of the smallest possible mammal is deduced and shown to be close to the mass of the shrew. Other approaches to exchange networks, assuming minimum energy expenditure and scale-invariance, have led to similar models . The model can be adapted, with no loss of rigour, to new data: Gillooly et al.  showed that the fractal supply network principle can be combined with simple Boltzmann kinetics to explain the effects of both body mass and temperature on metabolic rates. Since mass and temperature are the primary determinants of many physiological and ecological parameters, this work suggests that the model  could revolutionize biology.

This is an impressive range of successes. However, West and his co-workers make claims that are less compelling. The observation that cytochrome oxidase catalytic rates fit the same allometric curve as whole-organism metabolic rates is claimed as corroboration. However, cytochrome oxidase is not an organism, or a cell: it does not have a metabolic rate. It is also debatable whether mitochondria can be said to have "metabolic rates". (In contrast, Hochachka and Somero  noted that oxygen turnover in the whole biosphere can be fitted to the same curve but they recognized this as "a contingent fact with no biological significance".) Also, the explanation derived by West and his colleagues for the alleged body-mass-invariance of the metabolic rates of cultured cells (see earlier) is mathematically neat, but it leads to no experimentally testable predictions, and the heterogeneous data sources cited in this context make the explicandum itself unconvincing. Finally, the model is said to explain the quarter-power scalings of a wide range of biological variables other than metabolic rate, including population densities of trees  and carnivorous animals , plant growth rates, vascular network structure and maturation times [18, 108], and life-spans . It is not clear why any of these variables should depend on the fractal geometries of space-filling supply networks, still less on metabolic rates though there is widespread interest in the application of scaling laws in ecology, for instance in modelling biodiversity  and food webs .

Moreover, there are definite flaws in the model:-

(1) If West et al. were correct, maximal and standard metabolic rates should both scale as M 0.75 . The weight of evidence suggests that maximum rate in homoiotherms scales as M 0.88 (see earlier discussion [46–49] and following section).

(2) During maximal energy output by an organism, the supply of material is likely to be limiting. For example, in mammals, muscle contraction is responsible for most of the energy turnover at maximum output and it is generally believed that the rate is limited by oxygen supply (if anaerobic capacity is ignored). However, under standard metabolic rate conditions, energy demand is generally more significant, i.e. for the service and cellular maintenance functions mentioned previously. Therefore, it is not clear why the geometry and physics of the supply system should predict the allometric scaling of standard rather than maximal metabolic rate. ("Supply" and "demand" under conditions of maximal aerobic metabolism are complex terms because many physiological steps are involved. The extent to which each step limits the maximum metabolic rate might be quantifiable by a suitable extension of metabolic control analysis  this remains an active research area to which West et al. scarcely refer.)

(3) The mathematical derivations given in  are idealisations, but they do not seem to allow for large deviations from b = 0.75. However, there are often wide differences among empirical b values, as discussed earlier these were addressed in, for example,  and . Also, the model does not account, or allow, for the differences in allometric scaling among mammalian tissues and organs [66, 73, 80].

(4) West et al. accept that some of their proposed hierarchical supply networks might be "virtual" (as in mitochondria) rather than explicit (as in mammalian blood circulation), but it is not clear why such networks must always have the same geometry. For instance, why should the intracellular network discussed by Hochachka  show area-preserving branching? There is no evidence that it does. Moreover, the "flow" of reductants through mitochondria presumably takes place in the plane of the inner membrane, which has one dimension fewer than (say) the mammalian circulatory system, so even if mitochondria can be said to have "metabolic rates", the 0.75-power law cannot apply here yet, allegedly, it does apply.

These difficulties show that the West et al. model, despite its impressive economy, elegance, consistency and range, cannot be accepted unreservedly in its present form. The very generality, or "universality", of this model has made it suspect for some biologists . The implication that it reveals a long-suspected universal biological principle implicit in Kleiber's law has ensured its attraction for others .

### The model of Darveau and co-workers 

This group elaborated a multi-cause rather than a single-cause account of allometric scaling. Their "allometric cascade" model holds that each step in the physiological and biochemical pathways involved in ATP biosynthesis and utilization has its own scaling behaviour and makes its own contribution (defined by a control coefficient between 0 and 1) to the whole-organism metabolic rate. Thus, many linked steps rather than a single overarching principle account for Kleiber's law.

This idea is inherently plausible, and the model is attractive because it draws upon recent advances in metabolic control analysis in biochemistry  and physiology . It emphasises that standard metabolic rate is determined by energy demand, not supply and it predicts an exponent for maximal metabolic rate in mammals between 0.8 and 0.9, rather than 0.75, which agrees with experimental findings [46–49] and the data cited by Weibel . Implicitly – though the authors do not emphasize this – it seems capable of explaining b values that are far from 0.75 (cf ). It is hardly surprising, therefore, that many responses to the Darveau et al. model have been positive [e.g. ].

However, Darveau et al. made no attempt to explain why the values of b are typically around 0.75, as West et al. and others have done. The model is phenomenological, not physical and mathematical their equations are not derived from any fundamental principle(s). Moreover, their data cover only some three orders of magnitude of body mass, whereas many studies have involved much wider ranges. This might make their overall b values misleading  or, alternatively, more credible . When their equations are applied to a mass range of eight orders of magnitude, different b values are obtained, not necessarily consistent with published data but on the other hand, the published data might not be correct.

In the first published account of this model  the mathematical argument was flawed. The basic equation was given in the form B = aΣciM b(i) , where a is a constant coefficient, cI is the control coefficient of the i th step in the cascade and b(i) is the exponent of the i th step. By definition, the sum of all the cI values is unity. Darveau et al. did not derive this equation they stated it. They also stated that the overall exponent, the b term in the Kleiber equation, is a weighted average of all the individual b(i) values, the weighting being determined by the relevant control coefficients. It has been suggested that this leads to untenable inferences. For example, since the units of B and a are fixed, the units of cI must depend on those of b(i) but by definition, both b(i) and cI must be dimensionless. Also, according to the basic equation, the contribution made by each step to the overall metabolic rate depends on the units in which body mass is measured. If this criticism is valid then it is impossible to evaluate the model as it stands, because any attempt to align its predictions with experimental data would be meaningless. Another reservation about this model is that it does not purport to apply to all taxa, as the West et al. model does it relates only to metazoa, and in particular to homoiotherms. However, most of the relevant data in the literature concern homoiotherms.

A subsequent publication from this group  re-stated the basic equation in the form B = aΣcI(M/m) b(i) , where the constant a is described as the "characteristic metabolic rate" of an animal with characteristic body mass m. This eliminates the problem of mass units, because the mass term has been rendered dimensionless and it is mathematically simple to express control coefficients in dimensionless form. The revised equation might therefore be immune to some of the criticisms levelled at its predecessor. However, some of the earlier reservations remain: the equation remains phenomenological, not physical or geometrical and the restriction in its range of application is explicit. Nevertheless, these considerations by no means invalidate the model. Indeed, it is supported by data from experiments in exercise physiology .

The models of Darveau et al. [112, 115], Banavar et al. [94, 95] and West et al. [88, 102] all have attractive features but they all have flaws, and they cannot be reconciled with one another. If the positive contributions to biology that these models represent could be further developed, and their defects eliminated, could they be harmonized? If so, the advancement of our understanding would be considerable.

To order the book in print, please select one of the following:

This is an introductory book on metabolic control analysis available as of August 30th, 2018. Page content: 275 pages, 118 Illustrations, 71 Exercises. Retail Price \$49.95. A full-color version of the book in print and pdf form is also available from Gumroad:

A higher-quality printed copy is also available from:

This book is an introduction to control in biochemical pathways. It introduces students to some of the most important concepts in modern metabolic control. It covers the basics of metabolic control analysis that helps us think about how biochemical networks operate. The book should be suitable for undergraduates in their early to mid-years at college.

1. Traditional Concepts in Metabolic Regulation
2. Elasticities
3. Introduction to Biochemical Control
4. Linking the Parts to the Whole
5. Experimental Methods
6. Linear Pathways
7. Branched and Cyclic Systems
8. Negative Feedback
9. Stability
10. Stability of Negative Feedback Systems
11. Moiety Conservation Laws
12. Moiety Conserved Cycles

Appendix A: List of Symbols and Abbreviations
Appendix B: Control Equations