The stoichiometric matrix of a simplified metabolic network inBacillus Subtillis was constructed from the flux balance equations,which were used for reconciliation of the measured rates anddetermination of the inner m...The stoichiometric matrix of a simplified metabolic network inBacillus Subtillis was constructed from the flux balance equations,which were used for reconciliation of the measured rates anddetermination of the inner metabolic rates. Thus more reliableresults of the true and empirical maintenance coefficients wereobtained. The true maintenance coefficient is linearly related to thespecific growth rate and changes with the P/O ratio. The measuredbiomass yield of adenosine triphosphate (ATP) is also linearlyrelated to the P/O ratio.展开更多
Based on the gene-protein-reaction (GPR) model of S. cerevisiae_iND750 and the method of constraint-based analysis, we first calculated the metabolic flux distribution of S. cere-visiae_iND750. Then we calculated the ...Based on the gene-protein-reaction (GPR) model of S. cerevisiae_iND750 and the method of constraint-based analysis, we first calculated the metabolic flux distribution of S. cere-visiae_iND750. Then we calculated the deletion impact of 438 calculable genes, one by one, on the metabolic flux redistribution of S. cere-visiae_iND750. Next we analyzed the correlation between v (describing deletion impact of one gene) and d (connection degree of one gene) and the correlation between v and Vgene (flux sum controlled by one gene), and found that both of them were not of linear relation. Furthermore, we sought out 38 important genes that most greatly affected the metabolic flux distribution, and determined their functional subsystems. We also found that many of these key genes were related to many but not several subsystems. Because the in silico model of S. cere-visiae_iND750 has been tested by many ex-periments, thus is credible, we can conclude that the result we obtained has biological sig-nificance.展开更多
Elementary flux mode (EFM) analysis was used in the metabolic analysis of central carbon metabolism in Saccharomyces cerevisiae based on constructed cellular network. Calculated from the metabolic model, the ethanol...Elementary flux mode (EFM) analysis was used in the metabolic analysis of central carbon metabolism in Saccharomyces cerevisiae based on constructed cellular network. Calculated from the metabolic model, the ethanol-producing pathway No. 37 furthest converts the substrate into ethanol among the 78 elementary flux modes. The in silico metabolic phenotypes predicted based on this analysis fit well with the fermentation performance of the engineered strains, KAM3 and KAMll, which confirmed that EFM analysis is valid to direct the construction of Saccharomyces cerevisiae engineered strains, to increase the ethanol yield.展开更多
Metabolic modeling and machine learning(ML)are crucial components of the evolving next-generation tools in systems and synthetic biology,aiming to unravel the intricate relationship between genotype,phenotype,and the ...Metabolic modeling and machine learning(ML)are crucial components of the evolving next-generation tools in systems and synthetic biology,aiming to unravel the intricate relationship between genotype,phenotype,and the environment.Nonetheless,the comprehensive exploration of integrating these two frameworks,and fully harnessing the potential of fluxomic data,remains an unexplored territory.In this study,we present,rigorously evaluate,and compare ML-based techniques for data integration.The hybrid model revealed that the overexpression of six target genes and the knockout of seven target genes contribute to enhanced ethanol production.Specifically,we investigated the influence of succinate dehydrogenase(SDH)on ethanol biosynthesis in Saccharomyces cerevisiae through shake flask experiments.The findings indicate a noticeable increase in ethanol yield,ranging from 6%to 10%,in SDH subunit gene knockout strains compared to the wild-type strain.Moreover,in pursuit of a high-yielding strain for ethanol production,dual-gene deletion experiments were conducted targeting glycerol-3-phosphate dehydrogenase(GPD)and SDH.The results unequivocally demonstrate significant enhancements in ethanol production for the engineered strains Δsdh4Δgpd1,Δsdh5Δgpd1,Δsdh6Δgpd1,Δsdh4Δgpd2,Δsdh5Δgpd2,and Δsdh6Δgpd2,with improvements of 21.6%,27.9%,and 22.7%,respectively.Overall,the results highlighted that integrating mechanistic flux features substantially improves the prediction of gene knockout strains not accounted for in metabolic reconstructions.In addition,the finding in this study delivers valuable tools for comprehending and manipulating intricate phenotypes,thereby enhancing prediction accuracy and facilitating deeper insights into mechanistic aspects within the field of synthetic biology.展开更多
Background:Human immunodeficiency virus type 1(HIV-1)remains a persistent global health challenge.Therefore,a continuous exploration of novel therapeutic strategies is essential.A comprehensive understanding of how HI...Background:Human immunodeficiency virus type 1(HIV-1)remains a persistent global health challenge.Therefore,a continuous exploration of novel therapeutic strategies is essential.A comprehensive understanding of how HIV-1 utilizes the cellular metabolism machinery for replication can provide insights into new therapeutic approaches.Methods:In this study,we performed a flux balance analysis using a genome-scale metabolic model(GEM)integrated with an HIV-1 viral biomass objective function to identify potential targets for anti–HIV-1 interventions.We generated a GEM by integrating an HIV-1 production reaction into CD4+T cells and optimized for both host and virus optimal states as objective functions to depict metabolic profiles of cells in the status for optimal host biomass maintenance or for optimal HIV-1 virion production.Differential analysis was used to predict biochemical reactions altered optimal for HIV-1 production.In addition,we conducted in silico simulations involving gene and reaction knock-outs to identify potential anti–HIV-1 targets,which were subsequently validated by human phytohemagglutinin(PHA)blasts infected with HIV-1.Results:Differential analysis identified several altered biochemical reactions,including increased lysine uptake and oxidative phosphorylation(OXPHOS)activities in the virus optima compared with the host optima.In silico gene and reaction knock-out simulations revealed de novo pyrimidine synthesis,and OXPHOS could serve as potential anti–HIV-1 metabolic targets.In vitro assay confirmed that targeting OXPHOS using metformin could suppress the replication of HIV-1 by 56.6%(385.4±67.5 pg/mL in the metformintreated group vs.888.4±32.3 pg/mL in the control group,P<0.001).Conclusion:Our integrated host-virus genome-scale metabolic study provides insights on potential targets(OXPHOS)for anti-HIV therapies.展开更多
OBJECTIVE One of the long-expected goals of genome-scale metabolic modeling is to evaluate the influence of the perturbed enzymes to the yield of an expected end product.METHDOS Metabolic control analysis(MCA)performs...OBJECTIVE One of the long-expected goals of genome-scale metabolic modeling is to evaluate the influence of the perturbed enzymes to the yield of an expected end product.METHDOS Metabolic control analysis(MCA)performs such role to calculate the sensitivity of flux change upon that of enzymes under the framework of ordinary differential equation(ODE)models,which are restricted in small-scale networks and require explicit kinetic parameters.The constraint-based models,like flux balance analysis(FBA),lack of the room of performing MCA because they are parameters-free.In this study,we developed a hyper-cube shrink algorithm(HCSA)to incorporate the enzymatic properties to the FBA model by introducing a pair of parameters for each reaction.Our algorithm was able to handle not only prediction of knockout strains but also strains with an adjustment of expression level of certain enzymes.RESULTS We first showed the concept by applying HCSA to a simplest three-nodes network.Then we show the HCSA possesses Michaelis-Menten like behaviors characterized by steady state of ODE.We obtained good prediction of a synthetic network in Saccharomyces cerevisiae producing voilacein and analogues.Finally we showed its capability of predicting the flux distribution in genome-scale networks by applying it to sporulation in yeast.CONCLUSION We have developed an algorithm the impact on fluxes when certain enzymes were inhibited or activated.It provides us a powerful tool to evaluate the consequences of enzyme inhibitor or activator.展开更多
Soil heat flux is important for surface energy balance (SEB), and inaccurate estimation of soil heat flux often leads to surface energy imbalance. In this paper, by using observations of surface radiation fluxes and...Soil heat flux is important for surface energy balance (SEB), and inaccurate estimation of soil heat flux often leads to surface energy imbalance. In this paper, by using observations of surface radiation fluxes and soil temperature gradients at a semi-arid grassland in Xilingguole, Inner Mongolia, China from June to September 2008, the characters of the SEB for the semi-arid grassland were analyzed. Firstly, monthly averaged diurnal variations of SEB components were revealed. A 30-min forward phase displacement of soil heat flux (G) observed by a fluxplate at the depth of 5-em below the soil surface was conducted and its effect on the SEB was studied. Secondly, the surface soil heat flux (Gs) was computed by using harmonic analysis and the effect of the soil heat storage between the surface and the fluxplate on the SEB was examined. The results show that with the 30-min forward phase displacement of observed G, the slope of the ordinary linear regression (OLR) of turbulent fluxes (H+LE) against available energy (Rn G) increased from 0.835 to 0.842, i.e., the closure ratio of SEB increased by 0.7%, yet energy imclosure of 15.8% still existed in the SEB. When Gs, instead of G was used in the SEB equation, the slope of corresponding OLR of (H+LE) against (Rn-Gs) reached 0.979, thereby the imelosure ratio of SEB was reduced to only 2.1%.展开更多
Background Several enzymes and cofactors have been identified as contributing to the slow utilization of xylose by xylose-fermenting strains of Saccharomyces cerevisiae.However,there has been no consensus on which of ...Background Several enzymes and cofactors have been identified as contributing to the slow utilization of xylose by xylose-fermenting strains of Saccharomyces cerevisiae.However,there has been no consensus on which of these possible bottle-necks are the most important to address.A previous strain characterization study from our lab suggested that insufficient NAD+limits fermentation and may be the most important bottleneck affecting utilization of xylose for the production of ethanol.The development and validation of a genome scale dynamic flux balance model would help to verify the existence and extent of this and other metabolic bottlenecks and suggest solutions to guide future strain development thereby minimiz-ing bottleneck impact on process economics.Results A dynamic flux balance model was developed to identify bottlenecks in several strains of S.cerevisiae,both with wild-type pentose phosphate pathway expression and with the pathway over expressed.ZWF1 was found to be limiting in the oxidative portion of the pentose phosphate pathway under oxygen replete conditions.This pathway is used to regenerate NADPH.Under oxygen limiting conditions,respiration of xylose was limited by the lack of oxygen as a terminal electron acceptor.Ethanol production was also limited under these conditions due to the inability to balance NAD+/NADH.The model suggests the use of the anaplerotic glyoxylate pathway to improve NAD+/NADH balance,increasing ethanol produc-tion by 50%while producing succinate as a coproduct at upwards of 20 g/l.Conclusion In the production of high value chemicals from biomass,the use of the respiratory metabolism is a waste of feedstock carbon.Bottlenecks previously identified in the oxidative pentose phosphate pathway are currently only relevant under oxygen-replete conditions and cannot impact the partitioning of carbon between the respiratory and fermentative pathways.Focusing future efforts on the non-respiratory balancing of NAD+/NADH,perhaps through the glyoxylate pathway,would improve the economics of ethanol production both directly and through coproduct formation.展开更多
Flux balance analysis, based on the mass conservation law in a cellular organism, has been extensively employed to study the interplay between structures and functions of cellular metabolic networks. Consequently, the...Flux balance analysis, based on the mass conservation law in a cellular organism, has been extensively employed to study the interplay between structures and functions of cellular metabolic networks. Consequently, the phenotypes of the metabolism can be well elucidated. In this paper, we introduce the Expanded Flux Variability Analysis (EFVA) to characterize the intrinsic nature of metabolic reactions, such as flexibility, modularity and essentiality, by exploring the trend of the range, the maximum and the minimum flux of reactions. We took the metabolic network of Escherichia coli as an example and analyzed the variability of reaction fluxes under different growth rate constraints. The average variabil-ity of all reactions decreases dramatically when the growth rate increases. Consider the noise effect on the metabolic system, we thus argue that the microorganism may practically grow under a suboptimal state. Besides, under the EFVA framework, the reactions are easily to be grouped into catabolic and anabolic groups. And the anabolic groups can be further assigned to specific biomass constitute. We also discovered the growth rate dependent essentiality of reactions.展开更多
The metabolic network has become a hot topic in the area of system biology and flux-based analysis plays a very important role in understanding the characteristics of organism metabolic networks. We review mainly the ...The metabolic network has become a hot topic in the area of system biology and flux-based analysis plays a very important role in understanding the characteristics of organism metabolic networks. We review mainly the static methods for analyzing metabolic networks such as flux balance analysis (FBA), minimization of metabolic adjustment (MOMA), regulatory on / off minimization (ROOM), and dynamic flux balance analysis with linear quadratic regulator (DFBA-LQR). Then several kinds of commonly used software for flux analysis are introduced briefly and compared with each other. Finally, we highlight the applications of metabolic network flux analysis, especially its usage combined with other biological characteristics and its usage for drug design. The idea of combining the analysis of metabolic networks and other biochemical data has been gradually promoted and used in several aspects such as the combination of metabolic flux and the regulation of gene expression, the influence of protein evolution caused by metabolic flux, the relationship between metabolic flux and the topological characteristics, the optimization of metabolic engineering. More comprehensive and accurate properties of metabolic networks will be obtained by integrating metabolic flux analysis, network topological characteristics and dynamic modeling.展开更多
基金Supported by the Key Program of National Natural Science Foundation of China (No. 20036010) and the National Science Fund for Distinguished Young Scholars (No. 20028607).
文摘The stoichiometric matrix of a simplified metabolic network inBacillus Subtillis was constructed from the flux balance equations,which were used for reconciliation of the measured rates anddetermination of the inner metabolic rates. Thus more reliableresults of the true and empirical maintenance coefficients wereobtained. The true maintenance coefficient is linearly related to thespecific growth rate and changes with the P/O ratio. The measuredbiomass yield of adenosine triphosphate (ATP) is also linearlyrelated to the P/O ratio.
文摘Based on the gene-protein-reaction (GPR) model of S. cerevisiae_iND750 and the method of constraint-based analysis, we first calculated the metabolic flux distribution of S. cere-visiae_iND750. Then we calculated the deletion impact of 438 calculable genes, one by one, on the metabolic flux redistribution of S. cere-visiae_iND750. Next we analyzed the correlation between v (describing deletion impact of one gene) and d (connection degree of one gene) and the correlation between v and Vgene (flux sum controlled by one gene), and found that both of them were not of linear relation. Furthermore, we sought out 38 important genes that most greatly affected the metabolic flux distribution, and determined their functional subsystems. We also found that many of these key genes were related to many but not several subsystems. Because the in silico model of S. cere-visiae_iND750 has been tested by many ex-periments, thus is credible, we can conclude that the result we obtained has biological sig-nificance.
基金Supported by the National Natural Science Foundation of China (No.2002AA647040)
文摘Elementary flux mode (EFM) analysis was used in the metabolic analysis of central carbon metabolism in Saccharomyces cerevisiae based on constructed cellular network. Calculated from the metabolic model, the ethanol-producing pathway No. 37 furthest converts the substrate into ethanol among the 78 elementary flux modes. The in silico metabolic phenotypes predicted based on this analysis fit well with the fermentation performance of the engineered strains, KAM3 and KAMll, which confirmed that EFM analysis is valid to direct the construction of Saccharomyces cerevisiae engineered strains, to increase the ethanol yield.
基金financially supported by the National Natural Science Foundation of China(Grant NO.32071461)the National Key Research and Development Program of China(Grant NO.2019YFA0904300).
文摘Metabolic modeling and machine learning(ML)are crucial components of the evolving next-generation tools in systems and synthetic biology,aiming to unravel the intricate relationship between genotype,phenotype,and the environment.Nonetheless,the comprehensive exploration of integrating these two frameworks,and fully harnessing the potential of fluxomic data,remains an unexplored territory.In this study,we present,rigorously evaluate,and compare ML-based techniques for data integration.The hybrid model revealed that the overexpression of six target genes and the knockout of seven target genes contribute to enhanced ethanol production.Specifically,we investigated the influence of succinate dehydrogenase(SDH)on ethanol biosynthesis in Saccharomyces cerevisiae through shake flask experiments.The findings indicate a noticeable increase in ethanol yield,ranging from 6%to 10%,in SDH subunit gene knockout strains compared to the wild-type strain.Moreover,in pursuit of a high-yielding strain for ethanol production,dual-gene deletion experiments were conducted targeting glycerol-3-phosphate dehydrogenase(GPD)and SDH.The results unequivocally demonstrate significant enhancements in ethanol production for the engineered strains Δsdh4Δgpd1,Δsdh5Δgpd1,Δsdh6Δgpd1,Δsdh4Δgpd2,Δsdh5Δgpd2,and Δsdh6Δgpd2,with improvements of 21.6%,27.9%,and 22.7%,respectively.Overall,the results highlighted that integrating mechanistic flux features substantially improves the prediction of gene knockout strains not accounted for in metabolic reconstructions.In addition,the finding in this study delivers valuable tools for comprehending and manipulating intricate phenotypes,thereby enhancing prediction accuracy and facilitating deeper insights into mechanistic aspects within the field of synthetic biology.
基金the National Natural Science Foundation of China(82071784)the Fundamental Research Funds for the Central Universities(2042022dx0003 and PTPP2023002)+1 种基金the Key Research and Development Project of Hubei Province(2020BCA069)the Translational Medicine and Interdisciplinary Research Joint Fund of Zhongnan Hospital of Wuhan University(ZNJC202007).
文摘Background:Human immunodeficiency virus type 1(HIV-1)remains a persistent global health challenge.Therefore,a continuous exploration of novel therapeutic strategies is essential.A comprehensive understanding of how HIV-1 utilizes the cellular metabolism machinery for replication can provide insights into new therapeutic approaches.Methods:In this study,we performed a flux balance analysis using a genome-scale metabolic model(GEM)integrated with an HIV-1 viral biomass objective function to identify potential targets for anti–HIV-1 interventions.We generated a GEM by integrating an HIV-1 production reaction into CD4+T cells and optimized for both host and virus optimal states as objective functions to depict metabolic profiles of cells in the status for optimal host biomass maintenance or for optimal HIV-1 virion production.Differential analysis was used to predict biochemical reactions altered optimal for HIV-1 production.In addition,we conducted in silico simulations involving gene and reaction knock-outs to identify potential anti–HIV-1 targets,which were subsequently validated by human phytohemagglutinin(PHA)blasts infected with HIV-1.Results:Differential analysis identified several altered biochemical reactions,including increased lysine uptake and oxidative phosphorylation(OXPHOS)activities in the virus optima compared with the host optima.In silico gene and reaction knock-out simulations revealed de novo pyrimidine synthesis,and OXPHOS could serve as potential anti–HIV-1 metabolic targets.In vitro assay confirmed that targeting OXPHOS using metformin could suppress the replication of HIV-1 by 56.6%(385.4±67.5 pg/mL in the metformintreated group vs.888.4±32.3 pg/mL in the control group,P<0.001).Conclusion:Our integrated host-virus genome-scale metabolic study provides insights on potential targets(OXPHOS)for anti-HIV therapies.
基金The project supported by 985 Startup Funding in PKU
文摘OBJECTIVE One of the long-expected goals of genome-scale metabolic modeling is to evaluate the influence of the perturbed enzymes to the yield of an expected end product.METHDOS Metabolic control analysis(MCA)performs such role to calculate the sensitivity of flux change upon that of enzymes under the framework of ordinary differential equation(ODE)models,which are restricted in small-scale networks and require explicit kinetic parameters.The constraint-based models,like flux balance analysis(FBA),lack of the room of performing MCA because they are parameters-free.In this study,we developed a hyper-cube shrink algorithm(HCSA)to incorporate the enzymatic properties to the FBA model by introducing a pair of parameters for each reaction.Our algorithm was able to handle not only prediction of knockout strains but also strains with an adjustment of expression level of certain enzymes.RESULTS We first showed the concept by applying HCSA to a simplest three-nodes network.Then we show the HCSA possesses Michaelis-Menten like behaviors characterized by steady state of ODE.We obtained good prediction of a synthetic network in Saccharomyces cerevisiae producing voilacein and analogues.Finally we showed its capability of predicting the flux distribution in genome-scale networks by applying it to sporulation in yeast.CONCLUSION We have developed an algorithm the impact on fluxes when certain enzymes were inhibited or activated.It provides us a powerful tool to evaluate the consequences of enzyme inhibitor or activator.
基金Supported by the National Basic Research Program of China(2012CB955304)National Natural Science Foundation of China(40830957and40175008)China Postdoctoral Scientific Research Fund(20110490854)
文摘Soil heat flux is important for surface energy balance (SEB), and inaccurate estimation of soil heat flux often leads to surface energy imbalance. In this paper, by using observations of surface radiation fluxes and soil temperature gradients at a semi-arid grassland in Xilingguole, Inner Mongolia, China from June to September 2008, the characters of the SEB for the semi-arid grassland were analyzed. Firstly, monthly averaged diurnal variations of SEB components were revealed. A 30-min forward phase displacement of soil heat flux (G) observed by a fluxplate at the depth of 5-em below the soil surface was conducted and its effect on the SEB was studied. Secondly, the surface soil heat flux (Gs) was computed by using harmonic analysis and the effect of the soil heat storage between the surface and the fluxplate on the SEB was examined. The results show that with the 30-min forward phase displacement of observed G, the slope of the ordinary linear regression (OLR) of turbulent fluxes (H+LE) against available energy (Rn G) increased from 0.835 to 0.842, i.e., the closure ratio of SEB increased by 0.7%, yet energy imclosure of 15.8% still existed in the SEB. When Gs, instead of G was used in the SEB equation, the slope of corresponding OLR of (H+LE) against (Rn-Gs) reached 0.979, thereby the imelosure ratio of SEB was reduced to only 2.1%.
基金Partial funding for this study was provided through a multistate hatch Grant from Oregon State University Agricultural Experiment Station to the corresponding author.
文摘Background Several enzymes and cofactors have been identified as contributing to the slow utilization of xylose by xylose-fermenting strains of Saccharomyces cerevisiae.However,there has been no consensus on which of these possible bottle-necks are the most important to address.A previous strain characterization study from our lab suggested that insufficient NAD+limits fermentation and may be the most important bottleneck affecting utilization of xylose for the production of ethanol.The development and validation of a genome scale dynamic flux balance model would help to verify the existence and extent of this and other metabolic bottlenecks and suggest solutions to guide future strain development thereby minimiz-ing bottleneck impact on process economics.Results A dynamic flux balance model was developed to identify bottlenecks in several strains of S.cerevisiae,both with wild-type pentose phosphate pathway expression and with the pathway over expressed.ZWF1 was found to be limiting in the oxidative portion of the pentose phosphate pathway under oxygen replete conditions.This pathway is used to regenerate NADPH.Under oxygen limiting conditions,respiration of xylose was limited by the lack of oxygen as a terminal electron acceptor.Ethanol production was also limited under these conditions due to the inability to balance NAD+/NADH.The model suggests the use of the anaplerotic glyoxylate pathway to improve NAD+/NADH balance,increasing ethanol produc-tion by 50%while producing succinate as a coproduct at upwards of 20 g/l.Conclusion In the production of high value chemicals from biomass,the use of the respiratory metabolism is a waste of feedstock carbon.Bottlenecks previously identified in the oxidative pentose phosphate pathway are currently only relevant under oxygen-replete conditions and cannot impact the partitioning of carbon between the respiratory and fermentative pathways.Focusing future efforts on the non-respiratory balancing of NAD+/NADH,perhaps through the glyoxylate pathway,would improve the economics of ethanol production both directly and through coproduct formation.
基金Supported by the National Natural Science Foundation of China (Grant No. 10721403)National Basic Research Program of China (Grant Nos. 2006CB910706, 2007CB814800, 2009CB918500)Chun-Tsung endowment at Peking University and National Fund for Fostering Talents of Basic Science (Grant No. J0630311)
文摘Flux balance analysis, based on the mass conservation law in a cellular organism, has been extensively employed to study the interplay between structures and functions of cellular metabolic networks. Consequently, the phenotypes of the metabolism can be well elucidated. In this paper, we introduce the Expanded Flux Variability Analysis (EFVA) to characterize the intrinsic nature of metabolic reactions, such as flexibility, modularity and essentiality, by exploring the trend of the range, the maximum and the minimum flux of reactions. We took the metabolic network of Escherichia coli as an example and analyzed the variability of reaction fluxes under different growth rate constraints. The average variabil-ity of all reactions decreases dramatically when the growth rate increases. Consider the noise effect on the metabolic system, we thus argue that the microorganism may practically grow under a suboptimal state. Besides, under the EFVA framework, the reactions are easily to be grouped into catabolic and anabolic groups. And the anabolic groups can be further assigned to specific biomass constitute. We also discovered the growth rate dependent essentiality of reactions.
基金supported by the National Natural Science Foundation of China (30800199, 20773085 and 30770502)National High-Tech Research and Development Program of China (2007AA02Z333)
文摘The metabolic network has become a hot topic in the area of system biology and flux-based analysis plays a very important role in understanding the characteristics of organism metabolic networks. We review mainly the static methods for analyzing metabolic networks such as flux balance analysis (FBA), minimization of metabolic adjustment (MOMA), regulatory on / off minimization (ROOM), and dynamic flux balance analysis with linear quadratic regulator (DFBA-LQR). Then several kinds of commonly used software for flux analysis are introduced briefly and compared with each other. Finally, we highlight the applications of metabolic network flux analysis, especially its usage combined with other biological characteristics and its usage for drug design. The idea of combining the analysis of metabolic networks and other biochemical data has been gradually promoted and used in several aspects such as the combination of metabolic flux and the regulation of gene expression, the influence of protein evolution caused by metabolic flux, the relationship between metabolic flux and the topological characteristics, the optimization of metabolic engineering. More comprehensive and accurate properties of metabolic networks will be obtained by integrating metabolic flux analysis, network topological characteristics and dynamic modeling.