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 impact assessment for heterologous protein production in Streptomyces lividans based on genome-scale metabolic network modeling

The metabolic impact exerted on a microorganism due to heterologous protein production is still poorly understood in Streptomyces lividans. In this present paper, based on exometabolomic data, a proposed genome-scale metabolic network model is used to assess this metabolic impact in S. lividans. Constraint-based modeling results obtained in this work revealed that the metabolic impact due to heterologous protein production is widely distributed in the genome of S. lividans, causing both slow substrate assimilation and a shift in active pathways. Exchange fluxes that are critical for model performance have been identified for metabolites of mouse tumor necrosis factor, histidine, valine and lysine, as well as biomass. Our results unravel the interaction of heterologous protein production with intracellular metabolism of S. lividans, thus, a possible basis for further studies in relieving the metabolic burden via metabolic or bioprocess engineering.     

Metabolic flux analysis of CHO cells at growth and non-growth phases using isotopic tracers and mass spectrometry

Chinese hamster ovary (CHO) cells are the main platform for production of biotherapeutics in the biopharmaceutical industry. However, relatively little is known about the metabolism of CHO cells in cell culture. In this work, metabolism of CHO cells was studied at the growth phase and early stationary phase using isotopic tracers and mass spectrometry. CHO cells were grown in fed-batch culture over a period of six days. On days 2 and 4, [1,2-¹³C] glucose was introduced and the labeling of intracellular metabolites was measured by gas chromatography-mass spectrometry (GC–MS) at 6, 12 and 24 h following the introduction of tracer. Intracellular metabolic fluxes were quantified from measured extracellular rates and ¹³C-labeling dynamics of intracellular metabolites using non-stationary ¹³C-metabolic flux analysis (¹³C-MFA). The flux results revealed significant rewiring of intracellular metabolic fluxes in the transition from growth to non-growth, including changes in energy metabolism, redox metabolism, oxidative pentose phosphate pathway and anaplerosis. At the exponential phase, CHO cell metabolism was characterized by a high flux of glycolysis from glucose to lactate, anaplerosis from pyruvate to oxaloacetate and from glutamate to α-ketoglutarate, and cataplerosis though malic enzyme. At the stationary phase, the flux map was characterized by a reduced flux of glycolysis, net lactate uptake, oxidative pentose phosphate pathway flux, and reduced rate of anaplerosis. The fluxes of pyruvate dehydrogenase and TCA cycle were similar at the exponential and stationary phase. The results presented here provide a solid foundation for future studies of CHO cell metabolism for applications such as cell line development and medium optimization for high-titer production of recombinant proteins.     

Comparative analysis of high butanol tolerance and production in clostridia

2016, was the 100 years anniversary from launching of the first industrial acetone-butanol-ethanol (ABE) microbial production process. Despite this long period and also revival of scientific interest in this fermentative process over the last 20 years, solventogenic clostridia, mainly Clostridium acetobutylicum, Clostridium beijerinckii, Clostridium saccharoperbutylacetonicum and Clostridium pasteurianum, still have most of their secrets. One such poorly understood mechanism is butanol tolerance, which seems to be one of the most significant bottlenecks obstructing industrial exploitation of the process because the maximum achievable butanol concentration is only about 21 g/L. This review describes all the known cellular responses elicited by butanol, such as modifications of cell membrane and cell wall, formation of stress proteins, extrusion of butanol by efflux pumps, response of regulatory pathways, and also maps both random and targeted mutations resulting in high butanol production phenotypes. As progress in the field is inseparably associated with emerging methods, enabling a deeper understanding of butanol tolerance and production, progress in these methods, including genome mining, RNA sequencing and constructing of genome scale models are also reviewed. In conclusion, a comparative analysis of both phenomena is presented and a theoretical relationship is described between butanol tolerance/high production and common features including efflux pump formation/activity, stress protein production, membrane modifications and biofilm growth.     

Structural Basis for Catalysis of a Tetrameric Class IIa Fructose 1,6-Bisphosphate Aldolase from Mycobacterium tuberculosis

Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), currently infects one-third of the world's population in its latent form. The emergence of multidrug-resistant and extensive drug-resistant strains has highlighted the need for new pharmacological targets within M. tuberculosis. The class IIa fructose 1,6-bisphosphate aldolase (FBA) enzyme from M. tuberculosis (MtFBA) has been proposed as one such target since it is upregulated in latent TB. Since the structure of MtFBA has not been determined and there is little information available on its reaction mechanism, we sought to determine the X-ray structure of MtFBA in complex with its substrates. By lowering the pH of the enzyme in the crystalline state, we were able to determine a series of high-resolution X-ray structures of MtFBA bound to dihydroxyacetone phosphate, glyceraldehyde 3-phosphate, and fructose 1,6-bisphosphate at 1.5, 2.1, and 1.3 Å, respectively. Through these structures, it was discovered that MtFBA belongs to a novel tetrameric class of type IIa FBAs. The molecular details at the interface of the tetramer revealed important information for better predictability of the quaternary structures among the FBAs based on their primary sequences. These X-ray structures also provide interesting and new details on the reaction mechanism of class II FBAs. Substrates and products were observed in geometries poised for catalysis; in addition, unexpectedly, the hydroxyl-enolate intermediate of dihydroxyacetone phosphate was also captured and resolved structurally. These concise new details offer a better understanding of the reaction mechanisms for FBAs in general and provide a structural basis for inhibitor design efforts aimed at this class of enzymes.     

Regulatory networks,genes and glycerophospholipid biosynthesis pathway in schistosomiasis: A systems biology view for pharmacological intervention

Understanding network topology through embracing the global dynamical regulation of genes in an active state space rather than traditional one-gene–one trait approach facilitates the rational drug development process. Schistosomiasis, a neglected tropical disease, has glycerophospholipids as abundant molecules present on its surface. Lack of effective clinical solutions to treat pathogens encourages us to carry out systems-level studies that could contribute to the development of an effective therapy. Development of a strategy for identifying drug targets by combined genome-scale metabolic network and essentiality analyses through in silico approaches provides tantalizing opportunity to investigate the role of protein/substrate metabolism. A genome-scale metabolic network model reconstruction represents choline–phosphate cytidyltransferase as the rate limiting enzyme and regulates the rate of phosphatidylcholine (PC) biosynthesis. The uptake of choline was regulated by choline concentration, promoting the regulation of phosphocholine synthesis. In Schistosoma, the change in developmental stage could result from the availability of choline, hampering its developmental cycle. There are no structural reports for this protein. In order to inhibit the activity of choline–phosphate cytidyltransferase (CCT), it was modeled by homology modeling using 1COZ as the template from Bacillus subtilis. The transition-state stabilization and catalytic residues were mapped as ‘HXGH’ and ‘RTEGISTT’ motif. CCT catalyzes the formation of CDP-choline from phosphocholine in which nucleotidyltransferase adds CTP to phosphocholine. The presence of phosphocholine permits the parasite to survive in an immunologically hostile environment. This feature endeavors development of an inhibitor specific for cytidyltransferase in Schistosoma. Flavonolignans were used to inhibit this activity in which hydnowightin showed the highest affinity as compared to miltefosine.     

Computational analysis of phenotypic space in heterologous polyketide biosynthesis—Applications to Escherichia coli, Bacillus subtilis, and Saccharomyces cerevisiae

Polyketides represent a class of natural product small molecules with an impressive range of medicinal activities. In order to improve access to therapeutic polyketide compounds, heterologous metabolic engineering has been applied to transfer polyketide genetic pathways from often fastidious native hosts to more industrially-amenable heterologous hosts such as Escherichia coli, Saccharomyces cerevisiae, or Streptomyces coelicolor. Efforts thus far have resulted in titers either inferior to the native host and significantly below the theoretical yield, emphasizing the need to computationally investigate and engineer the interaction between native and heterologous metabolism for the improved production of heterologous polyketide compounds. In this work, we applied flux balance analysis on genome-scale models to simulate cellular metabolism and 6-deoxyerythronolide B (the cyclized polyketide precursor to erythromycin) production in three common heterologous hosts (E. coli, Bacillus subtilis, and S. cerevisiae) under a variety of carbon-source and medium compositions. We then undertook minimization of metabolic adjustment optimization to identify single and double gene-knockouts that resulted in increased polyketide production while maintaining cellular growth. For the production of 6-deoxyerythronolide B, the results suggest B. subtilis and E. coli are better heterologous hosts when compared to S. cerevisiae and that several single and multiple gene-knockout mutants are computationally predicted to improve specific production, in some cases, over 25-fold.     

Multivalent II [beta-D-Galp-(1-->4)-beta-D-GlcpNAc] and Talpha [beta-D-Galp-(1-->3)-alpha-D-GalpNAc] specific Moraceae family plant lectin from the seeds of Ficus bengalensis fruits

A galactose specific lectin was isolated from the seeds of Ficus bengalensis (Moraceae) fruits and designated as F. bengalensis agglutinin (FBA). The lectin was purified by affinity repulsion chromatography on fetuin-agarose and was a monomer of molecular mass 33kDa. Like other Moraceae family lectins, carbohydrate-binding activity of FBA was independent of any divalent cation. FBA did not bind with simple saccharides, however sugar ligands with aromatic aglycons showed pronounced binding. The combining site of FBA recognized preferably Galbeta1,4GlcNAcbeta1-(II) followed by Galbeta1,3GalNAcalpha1-(Talpha) containing glycotopes. Interaction with saccharides revealed that the combining site of FBA could well accommodate a tetrasaccharide, asialo GM1 glycan (Galbeta1,3GalNAcbeta1,4Galbeta1,4Glcbeta1-), whereas polyvalent Tn (GalNAcalpha1-Ser/Thr), one of the well-recognized ligands of Moraceae family lectin, did not interact well with FBA.     

The post-illumination chlorophyll fluorescence transient indicates the RuBP regeneration limitation of photosynthesis in low light in Arabidopsis

The mechanism of post-illumination chlorophyll fluorescence transient (PIFT) was investigated in Arabidopsis. PIFT was detected in the wild type after illumination with low light. In the fba3-2 (fructose-1,6-bisphosphate aldolase) mutant, in which PIFT is enhanced, strong light also induced PIFT. PIFT was suppressed not only in the triose phosphate/phosphate translocator (tpt-2) mutant, but also in tpt-2 fba3-2, suggesting that triose phosphates, such as dihydroxyacetone phosphate (DHAP), are involved in the PIFT mechanism. We concluded that PIFT is associated with ribulose-1,5-bisphosphate (RuBP)-regeneration limitation of photosynthesis in low light.