Network-level analysis of metabolic regulation in the human red blood cell using random sampling and singular value decomposition
© Barrett et al; licensee BioMed Central Ltd. 2006
Received: 13 August 2005
Accepted: 13 March 2006
Published: 13 March 2006
Extreme pathways (ExPas) have been shown to be valuable for studying the functions and capabilities of metabolic networks through characterization of the null space of the stoichiometric matrix (S). Singular value decomposition (SVD) of the ExPa matrix P has previously been used to characterize the metabolic regulatory problem in the human red blood cell (hRBC) from a network perspective. The calculation of ExPas is NP-hard, and for genome-scale networks the computation of ExPas has proven to be infeasible. Therefore an alternative approach is needed to reveal regulatory properties of steady state solution spaces of genome-scale stoichiometric matrices.
We show that the SVD of a matrix (W) formed of random samples from the steady-state solution space of the hRBC metabolic network gives similar insights into the regulatory properties of the network as was obtained with SVD of P. This new approach has two main advantages. First, it works with a direct representation of the shape of the metabolic solution space without the confounding factor of a non-uniform distribution of the extreme pathways and second, the SVD procedure can be applied to a very large number of samples, such as will be produced from genome-scale networks.
These results show that we are now in a position to study the network aspects of the regulatory problem in genome-scale metabolic networks through the use of random sampling.
In recent years, chemically-detailed metabolic networks capable of simulating growth have been reconstructed at the genome-scale for nearly twenty organisms . Any metabolic network can be represented as a stoichiometric matrix, and by assuming steady-state conditions one can compute the "space" in which all allowable steady-state flux distributions reside. This space can be spanned by a set of extreme pathways (ExPas) [2, 3], and for a given metabolic network this set is unique, invariant and irreducible . ExPas are so named because they form the edges of the steady-state solution space of the network and thus characterize its "extreme" functions. Mathematically, the ExPas are a set of conical basis vectors that forms a convex cone in high-dimensional space. In this manner, the extreme pathways circumscribe all of the possible steady-state flux distributions possible in the metabolic network. Biochemically, each ExPa corresponds to a steady-state flux distribution through the network, and is more than just a linear set of reactions linking substrate to product. Specifically, each ExPa describes the relative fluxes among a set of reactions needed for the conversion of substrate(s) into product(s), and does so by creating all byproducts needed to maintain systemic balances and by simultaneously maintaining all cofactor pools at steady state. The development of such network-based pathway definitions, including the elementary flux modes [5–7], has led to a systems view of network properties . Genome-scale metabolic pathway analysis [5, 9, 10] holds much promise for understanding fundamental aspects of network functionality, robustness, and regulation [11, 12], for metabolic engineering [13, 14], and for elucidating many other systems properties [12, 15–23].
The number of ExPas for a reaction network grows exponentially with network size, limiting their use for the analysis of large networks. Because of this limitation, ExPas have mainly been used to analyze realistic, but small, networks such as the hRBC . In further work [4, 24], singular value decomposition (SVD) was applied to the ExPa matrices (P), resulting in two key insights. First, it was shown that SVD of P can compute the effective dimensionality of the regulatory problem from a network perspective. Second, it was shown that the 'eigenpathways' that correspond to the modes of the SVD identified key network branch points that can represent key control points (reactions) in the network. These points constitute a minimal set of reactions whose control sets the steady state flux state for the entire network, represented as a specific location in the steady-state solution space. SVD rank orders these key branch points in accordance with the shape of the solution space and thus rank orders the effectiveness by which the branch points can be used to move the network's flux distribution around in the solution space. These features were demonstrated by showing that the dominant features of the first five eigenpathways correspond to actual major points of control in the hRBC and can effectively locate the solution in a 5-dimensional space . SVD of P can thus define the magnitude and nature of the regulatory problem associated with controlling a metabolic network.
Given the fundamental understanding of network properties that can be derived from the SVD of P, it would be natural to perform the same analysis on the much larger genome-scale metabolic networks that have already been reconstructed. Unfortunately, this is not currently possible due to the NP-hard nature of calculating ExPas. Even if computing ExPas were possible on a genome-scale, the likely result would be matrices containing extremely large numbers of ExPas, many more than would be needed using uniform random sampling. This number of ExPas presents an additional formidable hurdle for the SVD calculation, which is an O(n2) algorithm .
In this work, we demonstrate an approach that addresses this problem. As mentioned above, the set of ExPas for a metabolic network defines the shape of the steady-state solution space of the network. The information that SVD of the set of ExPas provides is a description of the space. The eigenvectors that are produced by the SVD indicate, in rank order, the directions corresponding to the greatest variance in the high-dimensional space, and the singular values give the relative magnitude of this variation. A potentially better approach is to characterize the shape of the solution spaces by generating uniform random samples within it. Each sample solution is used as a column in a matrix (W). With enough samples, a "picture" of the solution space emerges, which can then be analyzed using SVD in the same manner as was done previously for the ExPas. In this work, we performed uniform random sampling of the steady-state solution space of the hRBC metabolic network and applied SVD to the matrix, W, of these samples. We show that similar, though not identical, results are obtained through SVD of W as were obtained through SVD of P. Thus, this study provides a new method for the calculation of Principle Regulatory Modes (PRMs) that improves upon the previous method by being more tractable for genome-scale networks due to its independence from ExPas. Additionally, the interpretation of the results of the new method are more straightforward because it characterizes the shape of the solution space uniformly, without being affected by the non-uniformity of the distribution of ExPas on the edges of the solution space.
Ascertaining sampling sufficiency
Comparing the modes derived from P and W
Metabolite abbreviations used in Figures 2–6.
Reaction abbreviations used in Figures 2–6. The corresponding reaction stoichiometries are also included.
ADO + H2O → INO + NH3
Adenine phosphoribosyl transferase
PRPP + ADE + H2O → AMP + 2 Pi
ADO + ATP → ADP + AMP
FDP ↔ GA3P + DHAP
Adenosine monophosphate phosphohydrolase
AMP + H2O → ADO + Pi + H
Adenosine monophosphate deaminase
AMP + H2O → IMP + NH3
2 ADP ↔ ATP + AMP
23DPG + H2O → 3PG + Pi
13DPG → 23DPG + H
2PG ↔ PEP + H2O
G6P + NADP → 6PGL + NADPH + H
Glyceraldehyde phosphate dehydrogenase
GA3P + NAD + Pi ↔ 13DPG + NADH + H
Hypoxanthine guanine phosphoryl transferase
PRPP + HX + H2O → IMP + 2 Pi
GLU + ATP → G6P + ADP + H
IMP + H2O → INO + Pi + H
PYR + NADH + H ↔ LAC + NAD
6PGC + NADP → RL5P + NADPH + CO2
F6P + ATP → FDP + ADP + H
G6P ↔ F6P
13DPG + ADP ↔ 3PG + ATP
6PGL + H2O ↔ 6PGC + H
3PG ↔ 2PG
PEP + ADP + H → PYR + ATP
Purine nucleoside phosphorylase
INO + Pi ↔ HX + R1P
R1P ↔ R5P
Phosphoribosyl pyrophosphate synthetase
R5P + ATP → PRPP + AMP
RL5P ↔ R5P
GA3P + S7P ↔ E4P + F6P
X5P + R5P ↔ S7P + E4P
X5P + E4P ↔ F6P + GA3P
Triose phosphate isomerase
DHAP ↔ GA3P
RL5P ↔ X5P
Assessing the correspondence of inferred reaction regulatory importance
One of the primary insights resulting from the application of SVD to P was an understanding of the relative importance amongst all of the reactions in terms of regulatory control. The aim of this work was to show that similar insight can be obtained through uniform random sampling. To do this, we followed a previously developed procedure  and scored each reaction by the l2-norm (length) of its loading vector, which is a vector of the reaction's contribution to each of the eigenvectors resulting from the SVD procedure (see Methods section). The rationale for this scoring method is based on our goal of identifying those reactions most critical for controlling flux through the metabolic network. These control reactions are expected to have larger contributions to the eigenvectors resulting from the SVD procedure. The impact of these larger contributions is reflected in a larger absolute value of the loading vector, signifying a larger distance of the reaction from the origin. So in this context the l2-norm quantifying the length of the loading vector is appropriate for our purposes.
Scores allow us to impart a rank ordering on the reactions. The degree to which the SVD of W and the SVD of P provide similar information about the regulatory importance of the reactions, then, is measured by the similarity in which they rank order the reactions by their scores. We measured this similarity by computing the rank correlation coefficient between the two rank orderings of reactions. Just like Pearson's correlation coefficient, the rank correlation has a range [-1, 1]. We used two related rank correlation measures : Kendall's tau and Spearman's rho. We found the reaction rankings to be highly correlated, resulting in a tau value of 0.88 with P-value 1.1*10-15 and a rho value of 0.97 with P-value 4.0*10-12. Based on these results, we claim that the sampling methodology produces similar insight into identifying the most important reactions for flux control as did the ExPa-based approach.
Reconstruction of physiological flux distributions
In previous work  it was demonstrated that SVD of the extreme pathway matrix for the hRBC metabolic network revealed the essential dimensionality of the regulatory problem. Furthermore, it was shown that the dominant features of the eigenpathways corresponding to the five largest singular values corresponded to control points in the hRBC metabolic network. These results are significant because they demonstrate a potentially powerful method for developing network-based understanding of regulation in genome-scale metabolic networks. Unfortunately, the method relies on the computation of the extreme pathways for these networks – a computational exercise that is currently infeasible. In an effort to circumvent this obstacle, we investigated the utility of randomly sampling the steady-state solution space of the hRBC metabolic network, and through SVD of the matrix of sampled solution points arrived at comparable results to those obtained via SVD of P. We present three main results. First, we found that the sampling procedure did lead to the same estimates of the dimensionality of the regulatory problem for the metabolic network studied. Secondly, we demonstrated that the resulting PRMs reconstruct physiological flux distributions as efficiently as did the eigenpathways. Thirdly, we showed that the flux adjustments represented by the PRMs were similar, and similarly interpretable, as those represented by the eigenpathways. In particular, we found that the ranking of reactions by their regulatory "importance" in the two approaches was highly correlated with high statistical significance.
It is important to note that there is no a priori reason for believing that SVD of P would be superior to SVD of W. Thus, while the focus of this work was to explain the results of SVD of W in the context of what had already been shown for SVD of P, it is not necessary that an effort be made to force SVD of W to reproduce SVD of P. Instead, what was shown herein is that the approaches can be thought of as attacking the same problem and that SVD of W provides a way forward into genome-scale matrices to the line of inquiry that was begun with SVD of P, since sampling of genome-scale metabolic states has been accomplished at the genome-scale for E. coli and H. pylori [29, 30] – whereas the computation of P for these networks has proven infeasible. The results of the two methods will not generally be identical. In the case of SVD of P not only the shape of the space, but, more precisely, the distribution of the pathways along the edges of the space is what is being characterized. Because the samples are generated uniformly in the solution space, it is the shape of the solution space alone that is being characterized in this approach. Thus, the interpretation of SVD of W is more straightforward than is the interpretation of SVD of P.
We used physiological flux distributions computed from a kinetic model and compared how well PRMs and eigenpathways reconstructed them. As can be seen in Figure 7, the PRMs performed just as well, if not better, than the eigenpathways. Indeed, in three instances, the PRMs arrived at essentially complete reconstructions utilizing fewer PRMs. Based on these results, we estimate the regulatory problem for the hRBC to be roughly five-dimensional – in agreement with the estimate based on the eigenpathways.
The fractional projection of eigenpathways and PRMs onto one another.
The flux maps in Figures 2, 3, 4, 5, 6 demonstrate additional satisfying correspondence between PRMs and eigenpathways. First, the principal singular mode (first PRM) is a valid flux distribution, adhering to mass balance and thermodynamic constraints just as is the case with the principal eigenpathway. Second, the regulatory correspondence that was achieved with the eigenpathways is also possible with the PRMs. For example, the dominant feature of the second PRM is also the flux split between glycolysis and the pentose phosphate pathway. The third singular mode describes, similarly to the third eigenpathway, the glycolytic pathway with the production of ATP, but instead of proceeding through to pyruvate with the production of NADH, it proceeds through to lactate with the utilization of NADH to produce NAD. The fourth PRM, as can be seen in Figure 5, is closer to eigenpathway 5, and this is due to the dominant feature of pyruvate uptake and subsequent conversion to lactate. Finally, the fifth PRM shares the same dominant characterization as eigenpathway 4, this being the flux split between lower glycolysis and the Rapoport-Luebering shunt.
The PRMs of W have a slightly different interpretation than the eigenpathways derived from P. Mathematically, the modes of W are the orthogonal directions of maximal variance in the steady-state solution space of the stoichiometric matrix. Biochemically, they correspond to the directions of greatest flexibility in achieving a network-wide flux distribution. The unequal variances in these orthogonal directions imply a hierarchy, wherein control of the more dominant modes gives the cell more control of achieving a desired flux distribution.
From these results we conclude that SVD of a sample of uniformly obtained steady state flux distributions for a metabolic network leads to a similar means of characterizing the regulatory problem as was obtained previously with ExPas. This was a necessary demonstration on a small system, before applying the technique to genome-scale networks. With a number of genome-scale metabolic networks now significantly reconstructed for a number of organisms, the results presented in this work indicate that it should now be possible – given that feasible genome-scale sampling has been demonstrated [29, 30] – to realize the dimensionality and mechanisms of metabolic regulation implicit in the topology of these networks.
Human red blood cell model
For this work we utilized a previously-constructed kinetic model  and a stoichiometric model [22, 24] of the human red blood cell. The stoichiometric model accounts for 39 metabolites and 32 internal metabolic reactions, and also contains 12 primary exchange fluxes and 7 currency exchange fluxes.
Physiologic steady-state flux distributions
Four steady-state flux distributions were calculated using an in silico kinetic model , as detailed in earlier work , and used to represent physiological states of interest in the solution space. The first distribution was for an unstressed state or normal physiological state, corresponding to no metabolic load. The other three corresponded to the primary metabolic demands placed upon the hRBC based on its physiologic role in the human body. These demands were for ATP (energy load), NADPH (oxidative load), and NADH (met-hemoglobin load). In each case, the load was maximized in the kinetic model to the highest value possible where the model could still be solved. Thus, they represented as high a load as the network could handle. All these states were previously computed and used in a previous study . These maximum values were 0.37 mM/h for ATP load (above a mandatory ~1 mM/hr running of the Na+/K+ pump), 2.5 mM/h for the NADPH load, and 1.7 mM/h for the NADH load.
Physiologic flux constraints
The null space of the stoichiometric matrix of the hRBC metabolic network is a (mathematically) unbounded space. To transform this space into one that is physiologically relevant, we applied, as before, constraints by bounding the minimum and maximum fluxes for each reaction. All irreversible reactions were constrained with a minimum flux of 0 mM/h and a maximum flux of 2000 mM/h unless experimentally-determined maximum flux values were available. Such maximum reaction rates were available for 22 reactions, and have been published previously [22, 32]. Additionally, certain fluxes must have a minimum, non-zero value in order to sustain the hRBC. These values were obtained from a previous study on the hRBC [32, 33], and since they are few in number are given here: DPGase, 0.3 mM/h; ATP production, 1 mM/h; NADH production, 0.4 mM/h; and NADPH production, 0.05 mM/h.
Hit-and-Run sampling of the solution space
We use a variant of an artificial centering hit-and-run algorithm  to sample random flux distributions in the steady-state solution space of the hRBC. The algorithm is a random walk algorithm entailing three steps. The first step is to generate a set of "warmup" points in the solution space. This generation of warmup points was performed by slightly reducing the maximum flux and increasing the minimum flux for all reactions and then calculating a candidate flux distribution along a reaction coordinate using linear programming. Ten thousand such warmup points were created. The second step of the algorithm is to select an initial point from one of the "warm up" points generated. The warm up points were stored in a matrix W and were used to calculate an approximate center s of the solution space. The third step calculates the sample points. Beginning with an initial, randomly selected point x m from W, a direction vector is constructed by choosing a randomly selected point y from W and moving along the direction vector d = y-s from x m resulting in a new point xm+1. After each such move xm+1was randomly substituted into W and the approximate center recalculated. For the results presented herein, we saved a sample point every 1000 moves in the space. It should be noted that the red blood cell metabolic network does not have any non-trivial biochemical loops (type III ExPas) , and thus the samples generated do not violate "loop law" thermodynamic constraints [36, 37] that would otherwise need to be explicitly addressed in sampling . Note that in previous work , it was stated that the hRBC contains 16 Type III ExPas. The computation of ExPas in that work entailed the decoupling of reversible reactions into forward and backward reactions, and these resulting trivial loops were the 16 Type III pathways. In this work, reversible reactions were not decoupled, and so trivial Type III pathways were not an issue. One issue not encountered in this work which will is important for sampling genome-scale networks is the issue of non-trivial Type III ExPas (more than 2 reactions). Since Type III ExPas correspond to thermodynamically infeasible flux distributions, the generation of samples from genome-scale systems will need to be cognizant of this issue. This issue has been addressed elsewhere and will not be an impediment to sampling genome-scale metabolic networks .
The extreme pathways were calculated from the stoichiometric matrix of the hRBC using a previously published algorithm . Vmax values for the reactions in the hRBC have been previously published [22, 32]. Each extreme pathway was scaled by the most limiting Vmax amongst all reactions that it contained. All of the extreme pathways had at least one reaction with an experimentally defined Vmax.
Angle between vectors
Angles between modes were computed using standard linear algebra. For two vectors u and v, the angle was computed as Θ = cos-1 (u·v/||u|| ||v||).
Regulatory importance of reactions
The result of applying the SVD to either the matrix of ExPas or the matrix of uniformly sampled flux distributions is a set of eigenvectors. Each reaction has a corresponding loading vector, which is the vector of the reaction's contribution to each of the eigenvectors. Specifically, if vij is the contribution of reaction j to eigenvector i, and there are m eigenvectors, then the loading vector for reaction j is l j = (v1j, v2j, ..., v mj ). The score for reaction j is then given by the l2-norm of the loading vector: | l j | = (v1j2 + v2j2 + ... + v mj 2)1/2.
singular value decomposition
principal regulatory mode
This work was funded in part by NIHGMS grant number GM068837. BOP serves on the Scientific Advisory Board of Genomatica Inc.
- Reed JL, Famili I, Thiele I, Palsson BO: Towards multidimensional genome annotation. Nat Rev Genet 2006, 7: 130–141. 10.1038/nrg1769View ArticlePubMedGoogle Scholar
- Schilling CH, Letscher D, Palsson BO: Theory for the systemic definition of metabolic pathways and their use in interpreting metabolic function from a pathway-oriented perspective. J Theor Biol 2000, 203: 229–248. 10.1006/jtbi.2000.1073View ArticlePubMedGoogle Scholar
- Schilling CH, Palsson BO: The underlying pathway structure of biochemical reaction networks. Proc Natl Acad Sci U S A 1998, 95: 4193–4198. 10.1073/pnas.95.8.4193PubMed CentralView ArticlePubMedGoogle Scholar
- Price ND, Reed JL, Papin JA, Famili I, Palsson BO: Analysis of metabolic capabilities using singular value decomposition of extreme pathway matrices. Biophys J 2003, 84: 794–804.PubMed CentralView ArticlePubMedGoogle Scholar
- Schuster S, Fell DA, Dandekar T: A general definition of metabolic pathways useful for systematic organization and analysis of complex metabolic networks. Nat Biotechnol 2000, 18: 326–332. 10.1038/73786View ArticlePubMedGoogle Scholar
- Schuster S, Dandekar T, Fell DA: Detection of elementary flux modes in biochemical networks: a promising tool for pathway analysis and metabolic engineering. Trends Biotechnol 1999, 17: 53–60. 10.1016/S0167-7799(98)01290-6View ArticlePubMedGoogle Scholar
- Schuster S, Hilgetag C: On elementary flux modes in biochemical reaction systems at steady state. J Biol Syst 1994, 2: 165–182. 10.1142/S0218339094000131View ArticleGoogle Scholar
- Papin JA, Reed JL, Palsson BO: Hierarchical thinking in network biology: the unbiased modularization of biochemical networks. Trends Biochem Sci 2004, 29: 641–647. 10.1016/j.tibs.2004.10.001View ArticlePubMedGoogle Scholar
- Schilling CH, Schuster S, Palsson BO, Heinrich R: Metabolic pathway analysis: basic concepts and scientific applications in the post-genomic era. Biotechnol Prog 1999, 15: 296–303. 10.1021/bp990048kView ArticlePubMedGoogle Scholar
- Papin JA, Price ND, Wiback SJ, Fell DA, Palsson BO: Metabolic pathways in the post-genome era. Trends Biochem Sci 2003, 28: 250–258. 10.1016/S0968-0004(03)00064-1View ArticlePubMedGoogle Scholar
- Stelling J, Klamt S, Bettenbrock K, Schuster S, Gilles ED: Metabolic network structure determines key aspects of functionality and regulation. Nature 2002, 420: 190–193. 10.1038/nature01166View ArticlePubMedGoogle Scholar
- Price ND, Papin JA, Palsson BO: Determination of Redundancy and Systems Properties of Helicobacter pylori's Metabolic Network Using Genome-Scale Extreme Pathway Analysis. Genome Res 2002, 12: 760–769. 10.1101/gr.218002. Article published online before print in April 2002PubMed CentralView ArticlePubMedGoogle Scholar
- Liao JC, Hou SY, Chao YP: Pathway analysis, engineering and physiological considerations for redirecting central metabolism. Biotechnol Bioeng 1996, 52: 129–140. Publisher Full Text 10.1002/(SICI)1097-0290(19961005)52:1%3C129::AID-BIT13%3E3.0.CO;2-JView ArticlePubMedGoogle Scholar
- Carlson R, Fell D, Srienc F: Metabolic pathway analysis of a recombinant yeast for rational strain development. Biotechnol Bioeng 2002, 79: 121–134. 10.1002/bit.10305View ArticlePubMedGoogle Scholar
- Pfeiffer T, I SV, Nuno JC, Montero F, Schuster S: METATOOL: for studying metabolic networks. Bioinformatics 1999, 15: 251–257. 10.1093/bioinformatics/15.3.251View ArticlePubMedGoogle Scholar
- Schuster S, Fell DA, Pfeiffer T, Dandekar T, Bork P: Elementary Mode Analysis Illustrated with Human Red Cell Metabolism. In BioThermoKinetics in the Post Genomic Era. Edited by: Larsson D, Pahlman IL and Gustaffson L. Goeteborg, Chalmers; 1998:332–339.Google Scholar
- Reed JL, Palsson BO: Genome-Scale In Silico Models of E. coli Have Multiple Equivalent Phenotypic States: Assessment of Correlated Reaction Subsets That Comprise Network States. Genome Res 2004, 14: 1797–1805. 10.1101/gr.2546004PubMed CentralView ArticlePubMedGoogle Scholar
- Papin JA, Price ND, Palsson BO: Extreme pathway lengths and reaction participation in genome-scale metabolic networks. Genome Res 2002, 12: 1889–1900. 10.1101/gr.327702PubMed CentralView ArticlePubMedGoogle Scholar
- Schilling CH, Edwards JS, Letscher D, Palsson BO: Combining pathway analysis with flux balance analysis for the comprehensive study of metabolic systems. Biotechnol Bioeng 2000, 71: 286–306. 10.1002/1097-0290(2000)71:4<286::AID-BIT1018>3.0.CO;2-RView ArticlePubMedGoogle Scholar
- Famili I, Palsson BO: Systemic metabolic reactions are obtained by singular value decomposition of genome-scale stoichiometric matrices. Journal of Theoretical Biology 2003, 224: 87–96. 10.1016/S0022-5193(03)00146-2View ArticlePubMedGoogle Scholar
- Schilling CH, Palsson BO: Assessment of the Metabolic Capabilities of Haemophilus influenzae Rd through a Genome-scale Pathway Analysis. J Theor Biol 2000, 203: 249–283. 10.1006/jtbi.2000.1088View ArticlePubMedGoogle Scholar
- Wiback SJ, Palsson BO: Extreme pathway analysis of human red blood cell metabolism. Biophys J 2002, 83: 808–818.PubMed CentralView ArticlePubMedGoogle Scholar
- Schuster S, Pfeiffer T, Moldenhauer F, Koch I, Dandekar T: Exploring the pathway structure of metabolism: decomposition into subnetworks and application to Mycoplasma pneumoniae. Bioinformatics 2002, 18: 351–361. 10.1093/bioinformatics/18.2.351View ArticlePubMedGoogle Scholar
- Price ND, Reed JL, Papin JA, Wiback SJ, Palsson BO: Network-based Analysis of Metabolic Regulation in the Human Red Blood Cell. Journal of Theoretical Biology 2003, 225: 185–194. 10.1016/S0022-5193(03)00237-6View ArticlePubMedGoogle Scholar
- Golub GH, Kahan W: Calculating the Singular Values and Pseudoinverse of a Matrix. J of the Society for Industrial and Mathematics: Series B, Numerical Analysis 1965, 2: 205–224.Google Scholar
- Yoon H, Yang K, C. S: Feature Subset Selection and Feature Ranking for Multivariate Time Series. IEEE Transactions on Knowledge and Data Engineering 2005., 17:Google Scholar
- Kendall M, Gibbons JD: Rank Correlation Methods. 5th edition. 1990.Google Scholar
- Thiele I, Vo TD, Price ND, Palsson B: An Expanded Metabolic Reconstruction of Helicobacter pylori (iIT341 GSM/GPR): An in silico genome-scale characterization of single and double deletion mutants. J Bacteriol 2005, 187: 5818–5830. 10.1128/JB.187.16.5818-5830.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Almaas E, Kovacs B, Vicsek T, Oltvai ZN, Barabasi AL: Global organization of metabolic fluxes in the bacterium Escherichia coli. Nature 2004, 427: 839–843. 10.1038/nature02289View ArticlePubMedGoogle Scholar
- Price ND, Thiele I, B.O. P: Candidate states of Helicobacter pylori's genome-scale metabolic network uponapplication of loop law thermodynamic constraints. Biophsyical Journal 2006., accepted:Google Scholar
- Jamshidi N, Edwards JS, Fahland T, Church GM, Palsson BO: Dynamic simulation of the human red blood cell metabolic network. Bioinformatics 2001, 17: 286–287. 10.1093/bioinformatics/17.3.286View ArticlePubMedGoogle Scholar
- Joshi A, Palsson BO: Metabolic dynamics in the human red cell. Part III--Metabolic reaction rates. J theor biol 1990, 142: 41–68.View ArticlePubMedGoogle Scholar
- Price ND, Schellenberger J, Palsson BO: Uniform Sampling of Steady State Flux Spaces: Means to Design Experiments and to Interpret Enzymopathies. Biophysical Journal 2004, 87: 2172–2186. 10.1529/biophysj.104.043000PubMed CentralView ArticlePubMedGoogle Scholar
- Kaufman DE, Smith RL: Direction Choice for Accelerated Convergence in Hit-and-Run Sampling. Operations Research 1998., Jan-Feb:Google Scholar
- Price ND, Famili I, Beard DA, Palsson BO: Extreme pathways and Kirchhoff's second law. Biophys J 2002, 83: 2879–2882.PubMed CentralView ArticlePubMedGoogle Scholar
- Beard DA, Liang SD, Qian H: Energy balance for analysis of complex metabolic networks. Biophys J 2002, 83: 79–86.PubMed CentralView ArticlePubMedGoogle Scholar
- Qian H, Beard DA: Thermodynamics of stoichiometric biochemical networks in living systems far from equilibrium. Biophys Chem 2005, 114: 213–220. 10.1016/j.bpc.2004.12.001View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.