Enzyme-enzyme interactions and metabolite channelling: alternative mechanisms and their evolutionary significance

A cell is more than the sum of its (dilute) parts: A brief history of quinary structure

Most knowledge of protein structure and function is derived from experiments performed with purified protein resuspended in dilute, buffered solutions. However, proteins function in the crowded, complex cellular environment. Although the first four levels of protein structure provide important information, a complete understanding requires consideration of quinary structure. Quinary structure comprises the transient interactions between macromolecules that provides organization and compartmentalization inside cells. We review the history of quinary structure in the context of several metabolic pathways, and the technological advances that have yielded recent insight into protein behavior in living cells. The evidence demonstrates that protein behavior in isolated solutions deviates from behavior in the physiological environment.

The structural and functional coordination of glycolytic enzymes in muscle: evidence of a metabolon?

Metabolism sustains life through enzyme-catalyzed chemical reactions within the cells of all organisms. The coupling of catalytic function to the structural organization of enzymes contributes to the kinetic optimization important to tissue-specific and whole-body function. This coupling is of paramount importance in the role that muscle plays in the success of Animalia. The structure and function of glycolytic enzyme complexes in anaerobic metabolism have long been regarded as a major regulatory element necessary for muscle activity and whole-body homeostasis. While the details of this complex remain to be elucidated through in vivo studies, this review will touch on recent studies that suggest the existence of such a complex and its structure. A potential model for glycolytic complexes and related subcomplexes is introduced.

Biochemical Systems Theory: A Review

Biochemical systems theory (BST) is the foundation for a set of analytical andmodeling tools that facilitate the analysis of dynamic biological systems. This paper depicts major developments in BST up to the current state of the art in 2012. It discusses its rationale, describes the typical strategies and methods of designing, diagnosing, analyzing, and utilizing BST models, and reviews areas of application. The paper is intended as a guide for investigators entering the fascinating field of biological systems analysis and as a resource for practitioners and experts.

Physiological uncoupling of mitochondrial oxidative phosphorylation. Studies in different yeast species

Under non-phosphorylating conditions a high proton transmembrane gradient inhibits the rate of oxygen consumption mediated by the mitochondrial respiratory chain (state IV). Slow electron transit leads to production of reactive oxygen species (ROS) capable of participating in deleterious side reactions. In order to avoid overproducing ROS, mitochondria maintain a high rate of O(2) consumption by activating different exquisitely controlled uncoupling pathways. Different yeast species possess one or more uncoupling systems that work through one of two possible mechanisms: i) Proton sinks and ii) Non-pumping redox enzymes. Proton sinks are exemplified by mitochondrial unspecific channels (MUC) and by uncoupling proteins (UCP). Saccharomyces. cerevisiae and Debaryomyces hansenii express highly regulated MUCs. Also, a UCP was described in Yarrowia lipolytica which promotes uncoupled O(2) consumption. Non-pumping alternative oxido-reductases may substitute for a pump, as in S. cerevisiae or may coexist with a complete set of pumps as in the branched respiratory chains from Y. lipolytica or D. hansenii. In addition, pumps may suffer intrinsic uncoupling (slipping). Promising models for study are unicellular parasites which can turn off their aerobic metabolism completely. The variety of energy dissipating systems in eukaryote species is probably designed to control ROS production in the different environments where each species lives.

A flexible state-space approach for the modeling of metabolic networks II: advanced interrogation of hybridoma metabolism

Having previously introduced the mathematical framework of topological metabolic analysis (TMA) - a novel optimization-based technique for modeling metabolic networks of arbitrary size and complexity - we demonstrate how TMA facilitates unique methods of metabolic interrogation. With the aid of several hybridoma metabolic investigations as case-studies (Bonarius et al., 1995, 1996, 2001), we first establish that the TMA framework identifies biologically important aspects of the metabolic network under investigation. We also show that the use of a structured weighting approach within our objective provides a substantial modeling benefit over an unstructured, uniform, weighting approach. We then illustrate the strength of TAM as an advanced interrogation technique, first by using TMA to prove the existence of (and to quantitatively describe) multiple topologically distinct configurations of a metabolic network that each optimally model a given set of experimental observations. We further show that such alternate topologies are indistinguishable using existing stoichiometric modeling techniques, and we explain the biological significance of the topological variables appearing within our model. By leveraging the manner in which TMA implements metabolite inputs and outputs, we also show that metabolites whose possible metabolic fates are inadequately described by a given network reconstruction can be quickly identified. Lastly, we show how the use of the TMA aggregate objective function (AOF) permits the identification of modeling solutions that can simultaneously consider experimental observations, underlying biological motivations, or even purely engineering- or design-based goals.

A possible role of intracellular isoelectric focusing in the evolution of eukaryotic cells and multicellular organisms

A new scenario of the origin of eukaryotic cell and multicellularity is presented. A concentric pH-gradient has been shown to exist in the cytosol of eukaryotic cells. The most probable source of such gradient is its self-formation in gradient of electric field between center and periphery of a cell. Theoretical analysis has shown that, for example, a cell of Saccharomyces cerevisiae has enough energy to continuously sustain such gradient of strength about 1.5 kV/cm, the value sufficient for effective isoelectric focusing of cytoplasmic proteins. Focusing of enzymes could highly increase the effectiveness of an otherwise diffusion-limited metabolism of large cells by concentrating enzymes into small and distinct parts of a cytoplasm. By taking away an important physical constraint to the volume of cytoplasm, the intracellular isoelectric focusing enabled evolution of cells 3-4 order of magnitude larger than typical prokaryotic cells. This opened the way for the origin of phagocytosis and lately for the development of different forms of endosymbiosis, some of them resulting in an endosymbiotic origin of mitochondria and plastids. The large volume of a cell-enabled separation of nuclear and cytoplasmic compartments which was a precondition for separation of transcription and translation processes and therefore also for the origin of various RNA-preprocessing mechanisms. The possibility to regulate gene expression by postprocessing RNA and to regulate metabolism by an electrophoretic translocation enzymes between different parts of cytoplasm by changing their isoelectric points opened the way for cell and tissue differentiation and therefore for the origin of complex multicellular organisms.

Stochastic reaction-diffusion simulation of enzyme compartmentalization reveals improved catalytic efficiency for a synthetic metabolic pathway

We have demonstrated the accuracy of a spatial stochastic model of Escherichia coli central carbon metabolism using the next subvolume method (NSM), an efficient implementation of the Gillespie direct method of stochastic simulation. Using this model, we demonstrate that compartmentalization of the enzymes comprising an engineered pathway for biosynthesis of R-1,2-propanediol leads to improved kinetic properties for the pathway enzymes, especially when substrate diffusivities are low. Our results suggest that enzyme compartmentalization is a powerful approach for improving the catalytic turnover of a channeled carbon substrate and should be particularly useful when applied to synthetic metabolic pathways that suffer from poor translation efficiency, are present in highly variable copy numbers, and have low turnover for new substrates. Furthermore, this approach represents a generic modeling framework for simultaneously analyzing spatial and stochastic events in cellular metabolism and should enable quantitative evaluation of the effect of enzyme compartmentalization on virtually any recombinant pathway.

Current status and challenges in connecting models of erythrocyte metabolism to experimental reality

Detailed kinetic models of human erythrocyte metabolism have served to summarize the vast literature and to predict outcomes from laboratory and "Nature's" experiments on this simple cell. Mathematical methods for handling the large array of nonlinear ordinary differential equations that describe the time dependence of this system are well developed, but experimental methods that can guide the evolution of the models are in short supply. NMR spectroscopy is one method that is non-selective with respect to analyte detection but is highly specific with respect to their identification and quantification. Thus time courses of metabolism are readily recorded for easily changed experimental conditions. While the data can be simulated, the systems of equations are too complex to allow solutions of the inverse problem, namely parameter-value estimation for the large number of enzyme and membrane-transport reactions operating in situ as opposed to in vitro. Other complications with the modelling include the dependence of cell volume on time, and the rates of membrane transport processes are often dependent on the membrane potential. These matters are discussed in the light of new modelling strategies.

Mathematical modelling of the urea cycle. A numerical investigation into substrate channelling

Metabolite channelling, the process in which consecutive enzymes have confined substrate transfer in metabolic pathways, has been proposed as a biochemical mechanism that has evolved because it enhances catalytic rates and protects unstable intermediates. Results from experiments on the synthesis of radioactive urea [Cheung, C., Cohen, N.S. & Raijman, L (1989) J. Biol. Chem.264, 4038-4044] have been interpreted as implying channelling of arginine between argininosuccinate lyase and arginase in permeabilized hepatocytes. To investigate this interpretation further, a mathematical model of the urea cycle was written, using Mathematica it simulates time courses of the reactions. The model includes all relevant intermediates, peripheral metabolites, and subcellular compartmentalization. Analysis of the output from the simulations supports the argument for a high degree of, but not absolute, channelling and offers insights for future experiments that could shed more light on the quantitative aspects of this phenomenon in the urea cycle and other pathways.

Interactions between beta-enolase and creatine kinase in the cytosol of skeletal muscle cells

We studied interactions in vivo between the cytosolic muscle isoform of creatine kinase (M-CK) and the muscle isoform of 2-phospho-D-glycerate hydrolyase (beta-enolase) in muscle sarcoplasm by incubating glycerol-skinned fibres with FITC-labelled beta-enolase in the presence or absence of free CK. A small amount of bound beta-enolase was observed in the presence of large concentrations of CK. The mobility of enolase was measured in cultured satellite cells by modulated-fringe-pattern photobleaching. FITC-labelled beta-enolase was totally mobile in both the presence and the absence of CK but its diffusion coefficient was slightly lower in the presence of CK. This suggests a weak interaction in vivo between enolase and CK.

Effects of feedback inhibition on transit time in a linear pathway of Michaelis-Menten-type reactions

An analysis of the effects of external and internal metabolites on the steady-state behavior of linear pathways comprising a sequence of three Michaelis-Menten-type reactions with and without a simple feedback inhibition (i.e. an interaction of an internal metabolite with the pathway) is performed with respect to the transit time tau by its formulation as rectangular-hyperbolic functions of the flux J, instead of direct expressions in terms of the external metabolite concentrations. For a given concentration of the external metabolite M1 (substrate of the pathway) or M4 (product of the pathway), the flux J has a lower value in the pathway with feedback inhibition than in the pathway without feedback inhibition. With variation in the M1 concentration the transit time tau shows a concave relationship with the flux J which is virtually identical for both pathways, yielding a minimum at a certain value of J. With variation in the M4 concentration the transit time tau monotonously decreases with higher value of J, and for a given value of J the feedback inhibition allows a lower transit time. This effect is enhanced with stronger feedback inhibition, and is in turn greatly reduced with higher values of total concentration and rate constants for the first enzyme in the pathway.

A reciprocal allosteric mechanism for efficient transfer of labile intermediates between active sites in CAD, the mammalian pyrimidine-biosynthetic multienzyme polypeptide

Carbamoyl phosphate is the product of carbamoyl phosphate synthetase (CPS II) activity and the substrate of the aspartate transcarbamoylase (ATCase) activity, each of which is found in CAD, a large 240-kDa multienzyme polypeptide in mammals that catalyses the first three steps in pyrimidine biosynthesis. In our study of the transfer of the labile intermediate between the two active sites, we have used assays that differentiate the synthesis of carbamoyl phosphate from the overall reaction of CPS II and ATCase that produces carbamoyl aspartate. We provided excess exogenous carbamoyl phosphate and monitored its access to the respective active sites through the production of carbamoyl phosphate and carbamoyl aspartate from radiolabelled bicarbonate. Three features indicate interactions between the folded CPS II and ATCase domains causing reciprocal conformational changes. First, even in the presence of approximately 1 mM unlabelled carbamoyl phosphate, when the aspartate concentration is high ATCase uses endogenous carbamoyl phosphate for the synthesis of radiolabelled carbamoyl aspartate. In contrast, the isolated CPS II forward reaction is inhibited by excess unlabelled carbamoyl phosphate. Secondly, the affinity of the ATCase for carbamoyl phosphate and aspartate is modulated when substrates bind to CPS II. Thirdly, the transition-state analogue phosphonacetyl-L-aspartate is a less efficient inhibitor of the ATCase when the substrates for CPS II are present. All these effects operate when CPS II is in the more active P state, which is induced by high concentrations of ATP and magnesium ions and when 5'-phosphoribosyl diphosphate (the allosteric activator) is present with low concentrations of ATP; these are conditions that would be met during active biosynthesis in the cell. We propose a phenomenon of reciprocal allostery that encourages the efficient transfer of the labile intermediate within the multienzyme polypeptide CAD. In this model, binding of aspartate to the active site of ATCase causes a conformational change at the active site of the liganded form of CPS II, which protects it from inhibition by its product, carbamoyl phosphate; reciprocally, the substrates for CPS II affect the active site of ATCase by increasing the affinity for its substrates, endogenous carbamoyl phosphate and aspartate, and thus impede access of exogenous carbamoyl phosphate or the transition-state analogue. Reciprocal allostery justifies the close association of the enzyme activities within the polypeptide and ensures that carbamoyl phosphate is efficiently synthesised and is dedicated to the second step of pyrimidine biosynthesis. These conditions fulfill those required for metabolic channeling in the cell.

Interaction of the two proteins of the methoxylation system involved in cephamycin C biosynthesis. Immunoaffinity, protein cross-linking, and fluorescence spectroscopy studies

Cephamycin C-producing microorganisms contain a two-protein enzyme system that converts cephalosporins to 7-methoxycephalosporins. Interaction between the two component proteins P7 (Mr 27,000) and P8 (Mr 32,000) has been studied by immunoaffinity chromatography using anti-P7 and anti-P8 antibodies, cross-linking with glutaraldehyde, and fluorescence spectroscopy analysis. Co-renaturation of the P7 and P8 polypeptides resulted in the formation of a protein complex with a molecular mass of 59 kDa, which corresponds to a heterodimer of P7 and P8. Glutaraldehyde cross-linking of the polypeptides after assembly of the protein complex showed the presence of a single heterodimer form that reacted with antibodies against P7 and P8. Each separate protein did not associate with itself into multimers. The P7.P8 complex co-purified by immunoaffinity chromatography from extracts of Nocardia lactamdurans and Streptomyces clavuligerus, suggesting that both proteins are present as an aggregate in vivo. Fluorescence spectroscopy studies of 5-methylaminonaphthalene-1-sulfonyl-P7 in response to increasing concentrations of P8 showed a blue shift in the fluorophore emission, indicating a conformational change of P7 in response to the interaction of P8 with an apparent dissociation constant of 47 microM. NADH showed affinity for the P7 component. The P7.P8 complex interacted strongly with the substrates S-adenosylmethionine and cephalosporin C, differently from that occurring with the separate P7 or P8 components, resulting in a strong blue shift in the fluorescence emission spectra of the complex.

Application of short column gel permeation in the study of protein-protein interactions

A simple and rapid procedure based on the gel filtration principle is described together with its applicability to the study of protein-protein interactions including subunit-subunit and enzyme-enzyme interactions. Using this procedure, it is shown that phosphoglycerate kinase (PGK) and glyceraldehyde-3-phosphate dehydrogenase (GPDH) interact with a stoichiometry of one PGK molecule combining with one monomeric subunit of GPDH. This interaction has been observed with both enzymes being from the same, as well as from different, species. The Kd values for rabbit muscle PGK and porcine muscle GPDH complex and that for the rabbit muscle PGK and yeast GPDH complex are found to be (4.5 +/- 2.0) x 10(-7) M and (6.5 +/- 1.7) x 10(-7) M, respectively. The specificity of bienzyme association is stronger when enzymes are from the same species than when they are from different species.

Complex formation between mushroom tyrosinase and Manduca dopachrome isomerase

Melanin biosynthesis in animals is initiated by the ubiquitously present tyrosinase and is aided by dopachrome isomerase. We have characterized a novel dopachrome isomerase (decarboxylating) from the hemolymph of Manduca sexta that generates a new quinone methide intermediate during melanogenesis (Sugumaran, M. and Semensi, V. (1991) J. Biol. Chem. 266, 6073-6078). This enzyme has the ability to form a complex with mushroom tyrosinase as judged by a number of physicochemical studies. The isomerase exhibited a marked inhibitory effect on tyrosinase and tyrosinase reciprocated by inhibiting the isomerase. While the isomerase showed no activity toward preformed dopaminechrome, it readily influenced the stability of dopaminechrome generated in situ by tyrosinase. Moreover, mushroom tyrosinase, which lacked specific binding to Concanavalin A Sepharose column, after complexing with the isomerase exhibited binding to this column. The complex formation also affected the pI value as well as mobility on a size exclusion column of these enzymes. Enzymes executing sequential metabolic transformation are known to form complexes called metabolons. Based on these above studies, it is concluded that both the enzymes involved in insect melanogenic pathway--phenoloxidase and dopachrome isomerase--are able to form a metabolon complex.