Elementary vectors and autocatalytic sets for resource allocation in next-generation models of cellular growth

Traditional (genome-scale) metabolic models of cellular growth involve an approximate biomass “reaction”, which specifies biomass composition in terms of precursor metabolites (such as amino acids and nucleotides). On the one hand, biomass composition is often not known exactly and may vary drastically between conditions and strains. On the other hand, the predictions of computational models crucially depend on biomass. Also elementary flux modes (EFMs), which generate the flux cone, depend on the biomass reaction. To better understand cellular phenotypes across growth conditions, we introduce and analyze new classes of elementary vectors for comprehensive (next-generation) metabolic models, involving explicit synthesis reactions for all macromolecules. Elementary growth modes (EGMs) are given by stoichiometry and generate the growth cone. Unlike EFMs, they are not support-minimal, in general, but cannot be decomposed “without cancellations”. In models with additional (capacity) constraints, elementary growth vectors (EGVs) generate a growth polyhedron and depend also on growth rate. However, EGMs/EGVs do not depend on the biomass composition. In fact, they cover all possible biomass compositions and can be seen as unbiased versions of elementary flux modes/vectors (EFMs/EFVs) used in traditional models. To relate the new concepts to other branches of theory, we consider autocatalytic sets of reactions. Further, we illustrate our results in a small model of a self-fabricating cell, involving glucose and ammonium uptake, amino acid and lipid synthesis, and the expression of all enzymes and the ribosome itself. In particular, we study the variation of biomass composition as a function of growth rate. In agreement with experimental data, low nitrogen uptake correlates with high carbon (lipid) storage.     

High expression of cloned genes in E. coli and its consequences

Using recombinant DNA technology it is relatively easy to construct strains of bacteria which over-produce foreign proteins such as those from human, animal or viral sources. In the near future the same methodology will be used to increase the production of natural products of microorganism by manipulating metabolic pathways, e.g. those involved in antibiotic or amino acid syntheses. This article reviews the factors which lead to over-expression of proteins in E. coli and the possible consequences when large scale culture is undertaken.     

Metabolic pathway analysis of a recombinant yeast for rational strain development

Elementary mode analysis has been used to study a metabolic pathway model of a recombinant Saccharomyces cerevisiae system that was genetically engineered to produce the bacterial storage compound poly-beta-hydroxybutyrate (PHB). The model includes biochemical reactions from the intermediary metabolism and takes into account cellular compartmentalization as well as the reversibility/irreversibility of the reactions. The reaction network connects the production and/or consumption of eight external metabolites including glucose, acetate, glycerol, ethanol, PHB, CO(2), succinate, and adenosine triphosphate (ATP). Elementary mode analysis of the wild-type S. cerevisiae system reveals 241 unique reaction combinations that balance the eight external metabolites. When the recombinant PHB pathway is included, and when the reaction model is altered to simulate the experimental conditions when PHB accumulates, the analysis reveals 20 unique elementary modes. Of these 20 modes, 7 produce PHB with the optimal mode having a theoretical PHB carbon yield of 0.67. Elementary mode analysis was also used to analyze the possible effects of biochemical network modifications and altered culturing conditions. When the natively absent ATP citrate-lyase activity is added to the recombinant reaction network, the number of unique modes increases from 20 to 496, with 314 of these modes producing PHB. With this topological modification, the maximum theoretical PHB carbon yield increases from 0.67 to 0.83. Adding a transhydrogenase reaction to the model also improves the theoretical conversion of substrate into PHB. The recombinant system with the transhydrogenase reaction but without the ATP citrate-lyase reaction has an increase in PHB carbon yield from 0.67 to 0.71. When the model includes both the ATP citrate-lyase reaction and the transhydrogenase reaction, the maximum theoretical carbon yield increases to 0.84. The reaction model was also used to explore the possibility of producing PHB under anaerobic conditions. In the absence of oxygen, the recombinant reaction network possesses two elementary modes capable of producing PHB. Interestingly, both modes also produce ethanol. Elementary mode analysis provides a means of deconstructing complex metabolic networks into their basic functional units. This information can be used for analyzing existing pathways and for the rational design of further modifications that could improve the system's conversion of substrate into product.     

Metabolic flux analysis for optimizing the specific growth rate of recombinant <Emphasis Type="Italic">Aspergillus niger</Emphasis>

A comprehensive metabolic network comprising three intracellular compartments (cytoplasm, mitochondrion and peroxisome) was developed for Aspergillus niger. The metabolic flux network includes carbohydrate and amino acid metabolism in both anabolic and catabolic reactions. Linear programming was used for the optimization of the specific growth rates in combination with 37 measured input and output fluxes of the key metabolites to evaluate corresponding intracellular flux distributions throughout the batch fermentations. Logarithmic sensitivity analysis revealed that the addition of proline, alanine and glutamate benefited growth in defined media. The experimental observations and flux analysis showed that tyrosine was a potential candidate for biomass production improvement. Model predictions was verified by conducting batch and fed-batch fermentations and it was found that the addition of the four amino acids according to the predetermined schedule resulted in a 44 and 41% improvements in biomass and recombinant protein productions, respectively.     

Phenotypic characterization of <Emphasis Type="Italic">Corynebacterium glutamicum</Emphasis> under osmotic stress conditions using elementary mode analysis

Corynebacterium glutamicum, a soil bacterium, is used to produce amino acids such as lysine and glutamate. C. glutamicum is often exposed to osmolality changes in its medium, and the bacterium has therefore evolved several adaptive response mechanisms to overcome them. In this study we quantify the metabolic response of C. glutamicum under osmotic stress using elementary mode analysis (EMA). Further, we obtain the optimal phenotypic space for the synthesis of lysine and formation of biomass. The analysis demonstrated that with increasing osmotic stress, the flux towards trehalose formation and energy-generating pathways increased, while the flux of anabolic reactions diminished. Nodal analysis indicated that glucose-6-phosphate, phosphoenol pyruvate, and pyruvate nodes were capable of adapting to osmotic stress, whereas the oxaloacetic acid node was relatively unresponsive. Fewer elementary modes were active under stress indicating the rigid behavior of the metabolism in response to high osmolality. Optimal phenotypic space analysis revealed that under normal conditions the organism optimized growth during the initial log phase and lysine and trehalose formation during the stationary phase. However, under osmotic stress, the analysis demonstrated that the organism operates under suboptimal conditions for growth, and lysine and trehalose formation.     

The geometry of the flux cone of a metabolic network

The analysis of metabolic networks has become a major topic in biotechnology in recent years. Applications range from the enhanced production of selected outputs to the prediction of genotype-phenotype relationships. The concepts used are based on the assumption of a pseudo steady-state of the network, so that for each metabolite inputs and outputs are balanced. The stoichiometric network analysis expands the steady state into a combination of nonredundant subnetworks with positive coefficients called extremal currents. Based on the unidirectional representation of the system these subnetworks form a convex cone in the flux-space. A modification of this approach allowing for reversible reactions led to the definition of elementary modes. Extreme pathways are obtained with the same method but splitting up internal reactions into forward and backward rates. In this study, we explore the relationship between these concepts. Due to the combinatorial explosion of the number of elementary modes in large networks, we promote a further set of metabolic routes, which we call the minimal generating set. It is the smallest subset of elementary modes required to describe all steady states of the system. For large-scale networks, the size of this set is of several magnitudes smaller than that of elementary modes and of extreme pathways.     

A metabolic study to decipher amino acid catabolism-directed biofuel synthesis in Acetoanaerobium sticklandii DSM 519

Acetoanaerobium sticklandii DSM 519 is a hyper-ammonia-producing anaerobe. It has the ability to produce organic solvents and acids from protein catabolism through Stickland reactions and specialized pathways. Nevertheless, its protein catabolism-directed biofuel production has not yet been understood. The present study aimed to decipher such growth-associated metabolic potential of this organism at different growth phases using metabolic profiling. A seed culture of this organism was grown separately in metabolic assay media supplemented with gelatin and or a mixture of amino acids. The extracellular metabolites produced by this organism were qualitatively analyzed by gas chromatography–mass spectrometry platform. The residual amino acids after protein degradation and amino acids assimilation were identified and quantitatively measured by high-performance liquid chromatography (HPLC). Organic solvents and acids produced by this organism were detected and the quantity of them determined with HPLC. Metabolic profiling data confirmed the presence of amino acid catabolic products including tyramine, cadaverine, methylamine, and putrescine in fermented broth. It also found products including short-chain fatty acids and organic solvents of the Stickland reactions. It reported that amino acids were more appropriate for its growth yield compared to gelatin. Results of quantitative analysis of amino acids indicated that many amino acids either from gelatin or amino acid mixture were catabolised at a log-growth phase. Glycine and proline were poorly consumed in all growth phases. This study revealed that apart from Stickland reactions, a specialized system was established in A. sticklandii for protein catabolism-directed biofuel production. Acetone–butanol–ethanol (ABE), acetic acid, and butyric acid were the most important biofuel components produced by this organism. The production of these components was achieved much more on gelatin than amino acids. Thus, A. sticklandii is suggested herein as a potential organism to produce butyric acid along with ABE from protein-based wastes (gelatin) in bio-energy sectors.

     

Generation of recombinant antibodies and means for increasing their affinity

Highly specific interaction with foreign molecules is a unique feature of antibodies. Since 1975, when Keller and Milstein proposed the method of hybridoma technology and prepared mouse monoclonal antibodies, many antibodies specific to various antigens have been obtained. Recent development of methods for preparation of recombinant DNA libraries and in silico bioinformatics approaches for protein structure analysis makes possible antibody preparation using gene engineering approaches. The development of gene engineering methods allowed creating recombinant antibodies and improving characteristics of existing antibodies; this significantly extends the applicability of antibodies. By modifying biochemical and immunochemical properties of antibodies by changing their amino acid sequences it is possible to create antibodies with properties optimal for certain tasks. For example, application of recombinant technologies resulted in antibody preparation of high affinity significantly exceeding the initial affinity of natural antibodies. In this review we summarize information about the structure, modes of preparation, and application of recombinant antibodies and their fragments and also consider the main approaches used to increase antibody affinity.     

Conditions affecting production of functional muscle recombinant alpha-tropomyosin in Saccharomyces cerevisiae

Yeasts are attractive hosts for heterologous protein production as they follow the general eukaryotic post-translational modification pattern. The well-known Saccharomyces cerevisiae has been used to produce a large variety of foreign proteins. The proper function of muscle tropomyosin depends on a specific modification at its N-terminus. Although tropomyosin has been produced in different expression systems, only the recombinant protein produced in the yeast Pichia pastoris has native-like functional properties. In this paper we describe the production of functional skeletal muscle tropomyosin in the yeast S. cerevisiae. The recombinant protein was produced in high amounts and production was strongly affected by genetic and environmental factors, including plasmid copy number, promoter strength, and growth media composition.     

Reaction routes in biochemical reaction systems: algebraic properties,validated calculation procedure and example from nucleotide metabolism

 Elementary flux modes (direct reaction routes) are minimal sets of enzymes that can operate at steady state, with all irreversible reactions used in the appropriate direction. They can be interpreted as component pathways of a (bio)chemical reaction network. Here, two different definitions of elementary modes are given and their equivalence is proved. Several algebraic properties of elementary modes are then presented and proved. This concerns, amongst other features, the minimal number of enzymes of the network not used in an elementary mode and the situations where irreversible reactions are replaced by reversible ones. Based on these properties, a refined algorithm is presented, and it is formally proved that this algorithm will exclusively generate all the elementary flux modes of an arbitrary network containing reversible or irreversible reactions or both. The algorithm is illustrated by a biochemical example relevant in nucleotide metabolism. The computer implementation in two different programming languages is discussed. Received: 1 January 2001 / Revised version: 17 December 2001 / Published online: 17 July 2002     

Analysis of growth of Lactobacillus plantarum WCFS1 on a complex medium using a genome-scale metabolic model

A genome-scale metabolic model of the lactic acid bacterium Lactobacillus plantarum WCFS1 was constructed based on genomic content and experimental data. The complete model includes 721 genes, 643 reactions, and 531 metabolites. Different stoichiometric modeling techniques were used for interpretation of complex fermentation data, as L. plantarum is adapted to nutrient-rich environments and only grows in media supplemented with vitamins and amino acids. (i) Based on experimental input and output fluxes, maximal ATP production was estimated and related to growth rate. (ii) Optimization of ATP production further identified amino acid catabolic pathways that were not previously associated with free-energy metabolism. (iii) Genome-scale elementary flux mode analysis identified 28 potential futile cycles. (iv) Flux variability analysis supplemented the elementary mode analysis in identifying parallel pathways, e.g. pathways with identical end products but different co-factor usage. Strongly increased flexibility in the metabolic network was observed when strict coupling between catabolic ATP production and anabolic consumption was relaxed. These results illustrate how a genome-scale metabolic model and associated constraint-based modeling techniques can be used to analyze the physiology of growth on a complex medium rather than a minimal salts medium. However, optimization of biomass formation using the Flux Balance Analysis approach, reported to successfully predict growth rate and by product formation in Escherichia coli and Saccharomyces cerevisiae, predicted too high biomass yields that were incompatible with the observed lactate production. The reason is that this approach assumes optimal efficiency of substrate to biomass conversion, and can therefore not predict the metabolically inefficient lactate formation.     

Production of swine pepsinogen by protein-producing Bacillus brevis carrying swine pepsinogen cDNA

Summary An expression-secretion vector, pNU100, was constructed, utilizing the promoter and coding sequences for the signal peptide and nine amino-terminal amino acids of the middle wall protein, to produce foreign proteins by protein-producing Bacillus brevis. Expression of swine pepsinogen cDNA in B. brevis was examined with pNU100 as a vector. The recombinant swine pepsinogen synthesized by B. brevis was found to accumulate extracellularly in the form of a soluble protein and to have acid protease activity. The acid protease activity was completely inhibited by pepstatin. Furthermore, the recombinant pepsinogen was converted autocatalytically to pepsin under acidic conditions. This indicates that B. brevis produces a pepsinogen with the same conformation as authentic pepsinogen. Efficient production of the enzyme (11 mg/l) was achieved by regulating the pH of the medium. The enzyme produced by B. brevis remained stable on cultivation for a long period, up to 40 h. This is suggested to be due to a unique property of protein-producing B. brevis, i. e. a deficiency in extracellular protease production.     

Bioconjugation of therapeutic proteins and enzymes using the expanded set of genetically encoded amino acids

The last decade has witnessed striking progress in the development of bioorthogonal reactions that are strictly directed towards intended sites in biomolecules while avoiding interference by a number of physical and chemical factors in biological environment. Efforts to exploit bioorthogonal reactions in protein conjugation have led to the evolution of protein translational machineries and the expansion of genetic codes that systematically incorporate a range of non-natural amino acids containing bioorthogonal groups into recombinant proteins in a site-specific manner. Chemoselective conjugation of proteins has begun to find valuable applications to previously inaccessible problems. In this review, we describe bioorthogonal reactions useful for protein conjugation, and biosynthetic methods that produce proteins amenable to those reactions through an expanded genetic code. We then provide key examples in which novel protein conjugates, generated by the genetic incorporation of a non-natural amino acid and the chemoselective reactions, address unmet needs in protein therapeutics and enzyme engineering.     

Amino acid consumption in naïve and recombinant CHO cell cultures: producers of a monoclonal antibody

Most commercial media for mammalian cell culture are designed to satisfy the amino acid requirements for cell growth, but not necessarily those for recombinant protein production. In this study, we analyze the amino acid consumption pattern in naïve and recombinant Chinese hamster ovary (CHO) cell cultures. The recombinant model we chose was a CHO-S cell line engineered to produce a monoclonal antibody. We report the cell concentration, product concentration, and amino acid concentration profiles in naïve and recombinant cell cultures growing in CD OptiCHO™ medium with or without amino acid supplementation with a commercial supplement (CHO CD EfficientFeed™ B). We quantify and discuss the amino acid demands due to cell growth and recombinant protein production during long term fed batch cultivation protocols. We confirmed that a group of five amino acids, constituting the highest mass fraction of the product, shows the highest depletion rates and could become limiting for product expression. In our experiments, alanine, a non-important mass constituent of the product, is in high demand during recombinant protein production. Evaluation of specific amino acid demands could be of great help in the design of feeding/supplementation strategies for recombinant mammalian cell cultures.     

Enhancement of human granulocyte-colony stimulating factor production in recombinant E. coli using batch cultivation

Development of inexpensive and simple culture media is always favorable for recombinant protein over-expression in E. coli. The effects of medium composition on the production of recombinant human granulocyte-colony stimulating factor (rh-GCSF) were investigated in batch culture of E. coli BL21 (DE3) [pET23a-hgcsf]. First, the optimum medium for production of rh-GCSF was determined; and, then it was shown that mixture of amino acid addition at induction time, which was determined on the basis of amino acids frequency in the recombinant protein, increases recombinant protein expression level significantly. Furthermore, the effect of glucose concentration on productivity of rh-GCSF was investigated; 20 g/l of glucose will result in maximum attainable biomass and rh-GCSF in this process. At optimum conditions, a cell dry weight of 10.5 g/l, an expression level of about 35% of total cellular protein, rh-GCSF concentration of 1.75 ± 0.1 g/l, and overall rh-GCSF yield of 165 ± 5 mg/g were obtained.     

Genetic encoding of non‐natural amino acids in Drosophila melanogaster Schneider 2 cells

Insect cells are useful for the high‐yield production of recombinant proteins including chemokines and membrane proteins. In this study, we developed an insect cell‐based system for incorporating non‐natural amino acids into proteins at specific sites. Three types of promoter systems were constructed, and their efficiencies were compared for the expression of the prokaryotic amber suppressor tRNA^Tyr in Drosophila melanogaster Schneider 2 cells. When paired with a variant of Escherichia coli tyrosyl‐tRNA synthetase specific for 3‐iodo‐L ‐tyrosine, the suppressor tRNA transcribed from the U6 promoter most efficiently incorporated the amino acid into proteins in the cells. The transient and stable introductions of these prokaryotic molecules into the insect cells were then compared in terms of the yield of proteins containing non‐natural amino acids, and the “transient” method generated a sevenfold higher yield. By this method, 4‐azido‐L ‐phenylalanine was incorporated into human interleukin‐8 at a specific site. The yield of the azido‐containing IL‐8 was 1 μg/1 mL cell culture, and the recombinant protein was successfully labeled with a fluorescent probe by the Staudinger–Bertozzi reaction.     

Fed-batch high-cell-density fermentation strategies for Pichia pastoris growth and production

Pichia pastoris is extensively used to produce various heterologous proteins. Amounts of biopharmaceutical drugs and industrial enzymes have been successfully produced by fed-batch high-cell-density fermentation (HCDF) of this cell factory. High levels of cell mass in defined media can be easily achieved and therefore large quantities of recombinant proteins with enhanced activities and lower costs can be obtained through HCDF technology. A robust HCDF process makes a successful transition to commercial production. Recently, efforts have been made to increase the heterologous protein production and activity by the HCDF of P. pastoris. However, challenges around selecting a suitable HCDF strategy exist. The high-level expression of a specific protein in P. pastoris is still, at least in part, limited by optimizing the methanol feeding strategy. Here, we review the progress in developments and applications of P. pastoris HCDF strategies for enhanced expression of recombinant proteins. We focus on the methanol induction strategies for efficient fed-batch HCDF in bioreactors, mainly focusing on various stat-induction strategies, co-feeding, and the limited induction strategy. These processes control strategies have opened the door for expressing foreign proteins in P. pastoris and are expected to enhance the production of recombinant proteins.     

Immobilized palladium(II) ion affinity chromatography for recovery of recombinant proteins with peptide tags containing histidine and cysteine

Fusion of peptide‐based tags to recombinant proteins is currently one of the most used tools for protein production. Also, immobilized metal ion affinity chromatography (IMAC) has a huge application in protein purification, especially in research labs. The combination of expression systems of recombinant tagged proteins with this robust chromatographic system has become an efficient and rapid tool to produce milligram‐range amounts of proteins. IMAC‐Ni(II) columns have become the natural partners of 6xHis‐tagged proteins. The Ni(II) ion is considered as the best compromise of selectivity and affinity for purification of a recombinant His‐tagged protein. The palladium(II) ion is also able to bind to side chains of amino acids and form ternary complexes with iminodiacetic acid and free amino acids and other sulfur‐containing molecules. In this work, we evaluated two different cysteine‐ and histidine‐containing six amino acid tags linked to the N‐terminal group of green fluorescent protein (GFP) and studied the adsorption and elution conditions using novel eluents. Both cysteine‐containing tagged GFPs were able to bind to IMAC‐Pd(II) matrices and eluted successfully using a low concentration of thiourea solution. The IMAC‐Ni(II) system reaches less than 20% recovery of the cysteine‐containing tagged GFP from a crude homogenate of recombinant Escherichia coli, meanwhile the IMAC‐Pd(II) yields a recovery of 45% with a purification factor of 13. Copyright © 2014 John Wiley & Sons, Ltd.