A genome-scale metabolic reconstruction of Pseudomonas putida KT2440: iJN746 as a cell factory

Background

The study of biological interaction networks is a central theme of systems biology. Here, we investigate the relationships between two distinct types of interaction networks: the metabolic pathway map and the protein-protein interaction network (PIN). It has long been established that successive enzymatic steps are often catalyzed by physically interacting proteins forming permanent or transient multi-enzymes complexes. Inspecting high-throughput PIN data, it was shown recently that, indeed, enzymes involved in successive reactions are generally more likely to interact than other protein pairs. In our study, we expanded this line of research to include comparisons of the underlying respective network topologies as well as to investigate whether the spatial organization of enzyme interactions correlates with metabolic efficiency.

Results

Analyzing yeast data, we detected long-range correlations between shortest paths between proteins in both network types suggesting a mutual correspondence of both network architectures. We discovered that the organizing principles of physical interactions between metabolic enzymes differ from the general PIN of all proteins. While physical interactions between proteins are generally dissortative, enzyme interactions were observed to be assortative. Thus, enzymes frequently interact with other enzymes of similar rather than different degree. Enzymes carrying high flux loads are more likely to physically interact than enzymes with lower metabolic throughput. In particular, enzymes associated with catabolic pathways as well as enzymes involved in the biosynthesis of complex molecules were found to exhibit high degrees of physical clustering. Single proteins were identified that connect major components of the cellular metabolism and may thus be essential for the structural integrity of several biosynthetic systems.

Conclusion

Our results reveal topological equivalences between the protein interaction network and the metabolic pathway network. Evolved protein interactions may contribute significantly towards increasing the efficiency of metabolic processes by permitting higher metabolic fluxes. Thus, our results shed further light on the unifying principles shaping the evolution of both the functional (metabolic) as well as the physical interaction network.     

Physical interactions between tricarboxylic acid cycle enzymes in Bacillus subtilis: evidence for a metabolon

The majority of all proteins of a living cell is active in complexes rather than in an isolated way. These protein-protein interactions are of high relevance for many biological functions. In addition to many well established protein complexes an increasing number of protein-protein interactions, which form rather transient complexes has recently been discovered. The formation of such complexes seems to be a common feature especially for metabolic pathways. In the Gram-positive model organism Bacillus subtilis, we identified a protein complex of three citric acid cycle enzymes. This complex consists of the citrate synthase, the isocitrate dehydrogenase, and the malate dehydrogenase. Moreover, fumarase and aconitase interact with malate dehydrogenase and with each other. These five enzymes catalyze sequential reaction of the TCA cycle. Thus, this interaction might be important for a direct transfer of intermediates of the TCA cycle and thus for elevated metabolic fluxes via substrate channeling. In addition, we discovered a link between the TCA cycle and gluconeogenesis through a flexible interaction of two proteins: the association between the malate dehydrogenase and phosphoenolpyruvate carboxykinase is directly controlled by the metabolic flux. The phosphoenolpyruvate carboxykinase links the TCA cycle with gluconeogenesis and is essential for B. subtilis growing on gluconeogenic carbon sources. Only under gluconeogenic growth conditions an interaction of these two proteins is detectable and disappears under glycolytic growth conditions.     

Cysteine synthesis in plants: protein-protein interactions of serine acetyltransferase from Arabidopsis thaliana

The biosynthesis of cysteine represents the final step of sulfate assimilation in bacteria and plants. It is catalyzed by the sequential action of serine acetyltransferase (SAT) and O -acetylserine (thiol) lyase (OAS-TL) which form a cysteine synthase (CS) complex in vitro . SAT and OAS-TL from Arabidopsis thaliana have previously been cloned, and now the first evidence is presented for the CS complex and SAT self-interaction in vivo employing the yeast two-hybrid system. Application of this method proved to be an efficient tool for the analysis of protein-protein interactions within a plant metabolic protein complex. Mapping of SAT domain structure revealed two new, independent domains with specific functions in protein-protein interaction. Analysis using truncated proteins proved the C-terminus of SAT to be sufficient for association with OAS-TL and to correlate with the putative transferase activity domain. SAT/SAT interaction was localized in the central region of the protein and occured also between SAT isoforms. Both protein interaction domains coincided with distinct α-helical and β-sheet clusters and together correlated with the minimal protein structure required for SAT catalysis as shown by functional complementation of an Escherichia coli mutant. The homo- and hetero-oligomerization properties are discussed with respect to the assumed function of the CS complex in metabolic channeling and activation of SAT by interaction with OAS-TL.     

On the origin of intracellular compartmentation and organized metabolic systems

The history of the development of the ideas and research of organized metabolic systems during last three decades is shortly reviewed. The cell cytoplasm is crowded with solutes, soluble macromolecules such as enzymes, nucleic acids, structural proteins and membranes. The high protein density within the large compartments of the cells predominantly determines the major characteristics of cellular environment such as viscosity, diffusion and inhomogeneity. The fact that the solvent viscosity of cytoplasm is not substantially different from the water is explained by intracellular structural heterogeneity: the intrinsic macromolecular density is relatively low within the interstitial voids in the cell because many soluble enzymes are apparently integral parts of the insoluble cytomatrix and are not distributed homogeneously. The molecular crowding and sieving restrict the mobility of very large solutes, binding severely restrict the mobility of smaller solutes. One of consequence of molecular crowding and hindered diffusion is the need to compartmentalize metabolic pathway to overcome diffusive barriers. Although the movement of small molecules is slowed down in the cytoplasm, the metabolism can successfully proceed and even be facilitated by metabolite channeling which directly transfers the intermediate from one enzyme to an adjacent enzyme without the need of free aqueous-phase diffusion. The enhanced probability for intermediates to be transferred from one active site to the other by sequential enzymes requires stable or transient interactions of the relevant enzymes, which associate physically in non-dissociable, static multienzyme complexes--metabolones, particles containing enzymes of a part or whole metabolic systems. Therefore, within the living cell the metabolism depends on the structural organization of enzymes forming microcompartments. Since cells contain many compartments and microenvironments, the measurement of the concentration of metabolites in whole cells or tissues gives an average cellular concentration and not that which is actually sensed by the active site of a specific enzyme. Thus, the microcompartmentation could provide a mechanism which can control metabolic pathways. Independently and in parallel to the developments described above, the ideas of compartmentation came into existence from the necessity to explain important physiological phenomena, in particular in heart research and in cardiac electrophysiology. These phenomena demonstrated the physiological importance of the biophysical and biochemical mechanisms described in this review.     

A systematic study of yeast sterol biosynthetic protein-protein interactions using the split-ubiquitin system

Sterol biosynthesis occurs in the ER and most sterol biosynthetic enzymes have transmembrane domains. However, due to difficulties in characterizing membrane protein-protein interactions, the nature of the sterol biosynthetic complex as well as in vivo interactions between various enzymes have not been described. We employed a split-ubiquitin membrane protein yeast two-hybrid system to characterize interactions between sterol biosynthetic proteins. Fourteen bait constructs were co-transformed into a reporter yeast strain with 14 prey constructs representing all sterol enzymatic reactions beginning with the synthesis of squalene. Our results not only confirmed several previous interactions, but also allowed us to identify novel interactions. Based on these results, ergosterol biosynthetic enzymes display specific protein-protein interactions forming a functional complex we designate, the ergosome. In this complex, Erg11p, Erg25p, Erg27p, and Erg28p appear to form a core center that can interact with other enzymes in the pathway. Also Erg24p and Erg2p, two enzymes that are sensitive to morpholine antifungals, appear to interact with one another; however, the profile of protein interaction partners appears to be unique. Erg2p and Erg3p, two enzymes catalyzing sequential reactions also appear to have different interaction partners. Our results provide a working model as to how sterol biosynthetic enzymes are topologically organized not only in yeast but in plant and animal systems that share many of these biosynthetic reactions.     

Diversion of flux toward sesquiterpene production in Saccharomyces cerevisiae by fusion of host and heterologous enzymes

The ability to transfer metabolic pathways from the natural producer organisms to the well-characterized cell factory Saccharomyces cerevisiae is well documented. However, as many secondary metabolites are produced by collaborating enzymes assembled in complexes, metabolite production in yeast may be limited by the inability of the heterologous enzymes to collaborate with the native yeast enzymes. This may cause loss of intermediates by diffusion or degradation or due to conversion of the intermediate through competitive pathways. To bypass this problem, we have pursued a strategy in which key enzymes in the pathway are expressed as a physical fusion. As a model system, we have constructed several fusion protein variants in which farnesyl diphosphate synthase (FPPS) of yeast has been coupled to patchoulol synthase (PTS) of plant origin (Pogostemon cablin). Expression of the fusion proteins in S. cerevisiae increased the production of patchoulol, the main sesquiterpene produced by PTS, up to 2-fold. Moreover, we have demonstrated that the fusion strategy can be used in combination with traditional metabolic engineering to further increase the production of patchoulol. This simple test case of synthetic biology demonstrates that engineering the spatial organization of metabolic enzymes around a branch point has great potential for diverting flux toward a desired product.     

Spatial organization of enzymes for metabolic engineering

As synthetic pathways built from exogenous enzymes become more complicated, the probability of encountering undesired interactions with host organisms increases, thereby lowering product titer. An emerging strategy to combat this problem is to spatially organize pathway enzymes into multi-protein complexes, where high local concentrations of enzymes and metabolites may enhance flux and limit problematic interactions with the cellular milieu. Co-localizing enzymes using synthetic scaffolds has improved titers for multiple pathways. While lacking physical diffusion barriers, scaffolded systems could concentrate intermediates locally through a mechanism analogous to naturally occurring microdomains. A more direct strategy for compartmentalizing pathway components would be to encapsulate them within protein shells. Several classes of shells have been loaded with exogenous proteins and expressed successfully in industrial hosts. A critical challenge for achieving ideal pathway compartmentalization with protein shells will likely be evolving pores to selectively limit intermediate diffusion. Eventually, these tools should enhance our ability to rationally design metabolic pathways.     

The Bacillus subtilis cell division proteins FtsL and DivIC are intrinsically unstable and do not interact with one another in the absence of other septasomal components

The bacterial septum appears to comprise a macromolecular assembly of essential cell division proteins (the 'septasome') that are responsible for physically dividing the cell during cytokinesis. FtsL and DivIC are essential components of this division machinery in Bacillus subtilis. We used yeast two-hybrid analysis as well as a variety of biochemical and biophysical methods to examine the proposed interaction between Bacillus subtilis FtsL and DivIC. We show that FtsL and DivIC are thermodynamically unstable proteins that are likely to be unfolded and therefore targeted for degradation unless stabilized by interactions with other components of the septasome. However, we show that this stabilization does not result from a direct interaction between FtsL and DivIC. We propose that the observed interdependence of DivIC and FtsL stability is a result of indirect interactions that are mediated by other septasomal proteins.     

Production and characterization of bifunctional enzymes. Substrate channeling in the aspartate pathway

The direct channeling of an intermediate between enzymes that catalyze consecutive reactions in a pathway offers the possibility of an efficient, exclusive, and protected means of metabolite delivery. Aspartokinase-homoserine dehydrogenase I (AK-HDH I) from Escherichia coli is an unusual bifunctional enzyme in that it does not catalyze consecutive reactions. The potential channeling of the intermediate beta-aspartyl phosphate between the aspartokinase of this bifunctional enzyme and aspartate semialdehyde dehydrogenase (ASADH), the enzyme that catalyzes the intervening reaction, has been examined. The introduction of increasing levels of inactivated ASADH has been shown to compete against enzyme-enzyme interactions and direct intermediate channeling, leading to a decrease in the overall reaction flux through these consecutive enzymes. These same results are obtained whether these experiments are conducted with aspartokinase III, a naturally occurring monofunctional isozyme, with an artificially produced monofunctional aspartokinase I, or with a fusion construct of AK I-ASADH. These results provide definitive evidence for the channeling of beta-aspartyl phosphate between aspartokinase and aspartate semialdehyde dehydrogenase in E. coli and suggest that ASADH may provide a bridge to channel the intermediates between the non-consecutive reactions of AK-HDH I.     

Starch biosynthetic enzymes from developing maize endosperm associate in multisubunit complexes

Mutations affecting specific starch biosynthetic enzymes commonly have pleiotropic effects on other enzymes in the same metabolic pathway. Such genetic evidence indicates functional relationships between components of the starch biosynthetic system, including starch synthases (SSs), starch branching enzymes (BEs), and starch debranching enzymes; however, the molecular explanation for these functional interactions is not known. One possibility is that specific SSs, BEs, and/or starch debranching enzymes associate physically with each other in multisubunit complexes. To test this hypothesis, this study sought to identify stable associations between three separate SS polypeptides (SSI, SSIIa, and SSIII) and three separate BE polypeptides (BEI, BEIIa, and BEIIb) from maize (Zea mays) amyloplasts. Detection methods included in vivo protein-protein interaction tests in yeast (Saccharomyces cerevisiae) nuclei, immunoprecipitation, and affinity purification using recombinant proteins as the solid phase ligand. Eight different instances were detected of specific pairs of proteins associating either directly or indirectly in the same multisubunit complex, and direct, pairwise interactions were indicated by the in vivo test in yeast. In addition, SSIIa, SSIII, BEIIa, and BEIIb all comigrated in gel permeation chromatography in a high molecular mass form of approximately 600 kD, and SSIIa, BEIIa, and BEIIb also migrated in a second high molecular form, lacking SSIII, of approximately 300 kD. Monomer forms of all four proteins were also detected by gel permeation chromatography. The 600- and 300-kD complexes were stable at high salt concentration, suggesting that hydrophobic effects are involved in the association between subunits.     

Old pathway — new concept: control of glycolysis by metabolite- modulated dynamic enzyme associations

The presence and kinetic advantages of heterologous interactions between sequential enzymes has been established. Dynamic interactions between the sub-units of oligomeric enzymes may also affect enzymatic activity. If the interactions are specific for certain oligomeric forms of heterologous enzyme pairs, metabolic flux can be modulated by the state of association of enzymes controlled by participating alternatively in homologous and heterologous interactions. A molecular model including the regulatory role of specific metabolites is suggested for the control of glycolysis by pairwise dynamic enzyme associations.     

The Carbon Assimilation Network in Escherichia coli Is Densely Connected and Largely Sign-Determined by Directions of Metabolic Fluxes

Gene regulatory networks consist of direct interactions but also include indirect interactions mediated by metabolites and signaling molecules. We describe how these indirect interactions can be derived from a model of the underlying biochemical reaction network, using weak time-scale assumptions in combination with sensitivity criteria from metabolic control analysis. We apply this approach to a model of the carbon assimilation network in Escherichia coli. Our results show that the derived gene regulatory network is densely connected, contrary to what is usually assumed. Moreover, the network is largely sign-determined, meaning that the signs of the indirect interactions are fixed by the flux directions of biochemical reactions, independently of specific parameter values and rate laws. An inversion of the fluxes following a change in growth conditions may affect the signs of the indirect interactions though. This leads to a feedback structure that is at the same time robust to changes in the kinetic properties of enzymes and that has the flexibility to accommodate radical changes in the environment.     

Interactomic study on interaction between lipid droplets and mitochondria

An increasing body of evidence shows that the lipid droplet, a neutral lipid storage organelle, plays a role in lipid metabolism and energy homeostasis through its interaction with mitochondria. However, the cellular functions and molecular mechanisms of the interaction remain ambiguous. Here we present data from transmission electron microscopy, fluorescence imaging, and reconstitution assays, demonstrating that lipid droplets physically contact mitochondria in vivo and in vitro. Using a bimolecular fluorescence complementation assay in Saccharomyces cerevisiae, we generated an interactomic map of protein-protein contacts of lipid droplets with mitochondria and peroxisomes. The lipid droplet proteins Erg6 and Pet10 were found to be involved in 75% of the interactions detected. Interestingly, interactions between 3 pairs of lipid metabolic enzymes were detected. Collectively, these data demonstrate that lipid droplets make physical contacts with mitochondria and peroxisomes, and reveal specific molecular interactions that suggest active participation of lipid droplets in lipid metabolism in yeast.     

Evidence for iron channeling in the Fet3p-Ftr1p high-affinity iron uptake complex in the yeast plasma membrane

In high-affinity iron uptake in the yeast Saccharomyces cerevisiae, Fe(II) is oxidized to Fe(III) by the multicopper oxidase, Fet3p, and the Fe(III) produced is transported into the cell via the iron permease, Ftr1p. These two proteins are likely part of a heterodimeric or higher order complex in the yeast plasma membrane. We provide kinetic evidence that the Fet3p-produced Fe(III) is trafficked to Ftr1p for permeation by a classic metabolite channeling mechanism. We examine the (59)Fe uptake kinetics for a number of complexes containing mutant forms of both Fet3p and Ftr1p and demonstrate that a residue in one protein interacts with one in the other protein along the iron trafficking pathway as would be expected in a channeling process. We show that, as a result of some of these mutations, iron trafficking becomes sensitive to an added Fe(III) chelator that inhibits uptake in a strictly competitive manner. This inhibition is not strongly dependent on the chelator strength, however, suggesting that Fe(III) dissociation from the iron uptake complex, if it occurs, is kinetically slow relative to iron permeation. Metabolite channeling is a common feature of multifunctional enzymes. We constructed the analogous ferroxidase, permease chimera and demonstrate that it supports iron uptake with a kinetic pattern consistent with a channeling mechanism. By analogy to the Fe(III) trafficking that leads to the mineralization of the ferritin core, we propose that ferric iron channeling is a conserved feature of iron homeostasis in aerobic organisms.     

Biochemical and structural insights into bacterial organelle form and biogenesis

Many heterotrophic bacteria have the ability to make polyhedral structures containing metabolic enzymes that are bounded by a unilamellar protein shell (metabolosomes or enterosomes). These bacterial organelles contain enzymes associated with a specific metabolic process (e.g. 1,2-propanediol or ethanolamine utilization). We show that the 21 gene regulon specifying the pdu organelle and propanediol utilization enzymes from Citrobacter freundii is fully functional when cloned in Escherichia coli, both producing metabolosomes and allowing propanediol utilization. Genetic manipulation of the level of specific shell proteins resulted in the formation of aberrantly shaped metabolosomes, providing evidence for their involvement as delimiting entities in the organelle. This is the first demonstration of complete recombinant metabolosome activity transferred in a single step and supports phylogenetic evidence that the pdu genes are readily horizontally transmissible. One of the predicted shell proteins (PduT) was found to have a novel Fe-S center formed between four protein subunits. The recombinant model will facilitate future experiments establishing the structure and assembly of these multiprotein assemblages and their fate when the specific metabolic function is no longer required.     

Relationships among carbohydrate intermediate metabolites and DNA damage and repair in yeast from a systems biology perspective

Glucose and fructose are major dietary carbohydrates that are essential for general metabolism. The elevated consumption of these two monosaccharides by the human population is related to the development of pluri-metabolic syndromes (e.g., diabetes mellitus and obesity). Glucose and fructose are metabolized by specific biochemical pathways to generate energy and metabolites. Many of these metabolites are mono- and bi-phosphorylated compounds, which renders them likely to generate reactive carbonyl species (RCS). Under physiological conditions, RCS react non-enzymatically with macromolecules and small molecules by means of Maillard reactions, forming stable glycated/fructated compounds called advanced glycation end products (AGEs). DNA and dNTPs are prone to react with RCS, forming DNA- and dNTP-AGEs, and many of these compounds are genotoxic and/or mutagenic. Unfortunately, little is understood about the genotoxicity and/or mutagenicity of carbohydrate intermediate metabolites or their interactions with DNA repair and carbohydrate metabolic-associated proteins. To elucidate these associations between carbohydrate metabolic pathways, DNA repair mechanisms, and dNTP-/DNA-AGEs, a systems biology study was performed by employing algorithms to mine literature data and construct physical protein-protein interactions. The results obtained in this work allow us to construct a model suggesting that yeast carbohydrate metabolic-associated enzymes activate different mechanisms for DNA repair and dNTP synthesis and act during DNA replication to protect the genome against the effects of RCS.     

Transmembrane homodimerization of receptor-like protein tyrosine phosphatases

Receptor-like protein tyrosine phosphatases (RPTPs) are type I integral membrane proteins. Together with protein tyrosine kinases, RPTPs regulate the phosphotyrosine levels in the cell. Studies of two RPTPs, CD45 and PTPalpha, have provided strong evidence that dimerization leads to inactivation of the receptors, and that the dimerization of PTPalpha involves interactions in the transmembrane domain (TMD). Using the TOXCAT assay, a genetic approach for analyzing TM interactions in Escherichia coli membranes, we show that the TMD of RPTPs interact in the membrane, albeit to different extents. Using fusion proteins of TMDs, we also observe an equilibrium between monomer and dimer in sodium dodecyl sulfate (SDS) micelles. Through a mutational study of the DEP1 TMD, we demonstrate that these interactions are specific. Taken together, our results define a subset of the RPTP family in which TM homodimerization may act as a mediator of protein function.     

Temporal Regulation of the Metabolome and Proteome in Photosynthetic and Photorespiratory Pathways Contributes to Maize Heterosis

Heterosis or hybrid vigor is widespread in plants and animals. Although the molecular basis for heterosis has been extensively studied, metabolic and proteomic contributions to heterosis remain elusive. Here we report an integrative analysis of time-series metabolome and proteome data in maize (Zea mays) hybrids and their inbred parents. Many maize metabolites and proteins are diurnally regulated, and many of these show nonadditive abundance in the hybrids, including key enzymes and metabolites involved in carbon assimilation. Compared with robust trait heterosis, metabolic heterosis is relatively mild. Interestingly, most amino acids display negative mid-parent heterosis (MPH), i.e., having lower values than the average of the parents, while sugars, alcohols, and nucleoside metabolites show positive MPH. From the network perspective, metabolites in the photosynthetic pathway show positive MPH, whereas metabolites in the photorespiratory pathway show negative MPH, which corresponds to nonadditive protein abundance and enzyme activities of key enzymes in the respective pathways in the hybrids. Moreover, diurnally expressed proteins that are upregulated in the hybrids are enriched in photosynthesis-related gene-ontology terms. Hybrids may more effectively remove toxic metabolites generated during photorespiration, and thus maintain higher photosynthetic efficiency. These metabolic and proteomic resources provide unique insight into heterosis and its utilization for high yielding maize and other crop plants.