The tricarboxylic acid cycle in Dictyostelium discoideum. III. Analysis of steady state and dynamic behavior.

The examination of model robustness in the previous paper (Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22919-22925 led to the suggestion that the current model for the tricarboxylic acid cycle in Dictyostelium discoideum is ill-determined with respect to one or more of the features reflecting pyruvate metabolism. This conclusion is further supported here by results of steady state and dynamic analyses. The tricarboxylic acid cycle, according to the current model, is poised on a knife's edge with its behavior rigidly determined; any alteration of the system's components leads to nonviable behavior, as exemplified by explosive accumulation of pyruvate and loss of steady state in response to a minute change in the level of malate dehydrogenase. With the additional results in this paper, we are able to refine the diagnosis of the problem and suggest three different areas of the current model that might profitably be re-examined by experiment. These include the kinetics of the reactions at the malate branch point, the turnover times for the alanine, glutamate, and aspartate pools in vivo, and the dynamic mass balances for the cofactor NAD. We also suggest a minimal modification in the current model that could alleviate or circumvent some of these problems.     

The tricarboxylic acid cycle in Dictyostelium discoideum. II. Evaluation of model consistency and robustness.

When kinetic models of complex biochemical systems are reconstructed from knowledge of the component reactions that have been characterized in vitro, or when values must be assumed for some of the parameters, errors are invariably encountered, and, as a consequence, the resulting model is frequently internally inconsistent. The simplest and most basic manifestations of such logical inconsistency are the failure of the model to exhibit a steady state or to yield a steady state that is in agreement with the actual steady state of the integrated system, or to yield a steady state that is dynamically stable. Models that are consistent may nonetheless be lacking in robustness, which is manifested as a pathological sensitivity to small changes in the values of their parameters. In this paper, we examine the current model of the tricarboxylic acid cycle in Dictyostelium discoideum (see Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22912-22918) with regard to these basic indicators of model quality. This may be viewed as a preliminary analysis; the object is to determine whether or not the model is reasonable and worthy of a more refined analysis and, if not, to diagnose the areas in need of modification before further analysis is undertaken. The results demonstrate that the current model of the tricarboxylic acid cycle is self-consistent and possesses a steady state that is in agreement with experimental evidence. However, the results also suggest that this model is not very robust. The high sensitivities of parameters influencing pyruvate metabolism indicate that the experimental characterization of these reactions might be fruitfully re-examined. These high sensitivities lead us to predict that this model of the tricarboxylic acid cycle should be accurate only over a very narrow range in variation of the independent variables. This is verified by the results presented in the following paper (Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22926-22933).     

The tricarboxylic acid cycle in Dictyostelium discoideum. I. Formulation of alternative kinetic representations.

Enzyme systems within living cells have recently been shown to be highly ordered structures that violate classic assumptions of the Michaelis-Menten formalism, which originally was developed for the characterization of isolated reactions in vitro. This evidence suggests that a thorough examination of alternative kinetic formalisms for integrated biochemical systems is in order. The purpose of this series of papers is to assess the utility of an alternative power-law formalism by carrying out a detailed comparative analysis of a relatively large, representative system--the tricarboxylic acid cycle of Dictyostelium discoideum. This system was chosen because considerable experimental information already has been synthesized into a detailed kinetic model of the intact system. In this first paper, we set the stage for subsequent analysis within the framework of the power-law formalism: we review the underlying theory, emphasizing recent developments, formulate the model in terms that are convenient for the analysis to follow, and develop the system representation in both the Michaelis-Menten and power-law forms. In the second paper (Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22919-22925), these alternative representations are shown to be internally consistent and locally equivalent. The third paper (Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22926-22933) provides a complete analysis of the steady state behavior and also treats the dynamic behavior of the model.     

Engineering cytochrome c peroxidase into cytochrome P450: a proximal effect on heme-thiolate ligation.

In an effort to investigate factors required to stabilize heme-thiolate ligation, key structural components necessary to convert cytochrome c peroxidase (CcP) into a thiolate-ligated cytochrome P450-like enzyme have been evaluated and the H175C/D235L CcP double mutant has been engineered. The UV-visible absorption, magnetic circular dichroism (MCD) and electron paramagnetic resonance (EPR) spectra for the double mutant at pH 8.0 are reported herein. The close similarity between the spectra of ferric substrate-bound cytochrome P450cam and those of the exogenous ligand-free ferric state of the double mutant with all three techniques support the conclusion that the latter has a pentacoordinate, high-spin heme with thiolate ligation. Previous efforts to prepare a thiolate-ligated mutant of CcP with the H175C single mutant led to Cys oxidation to cysteic acid [Choudhury et al. (1994) J. Biol. Chem. 267, 25656-25659]. Therefore it is concluded that changing the proximal Asp235 residue to Leu is critical in forming a stable heme-thiolate ligation in the resting state of the enzyme. To further probe the versatility of the CcP double mutant as a ferric P450 model, hexacoordinate low-spin complexes have also been prepared. Addition of the neutral ligand imidazole or of the anionic ligand cyanide results in formation of hexacoordinate adducts that retain thiolate ligation as determined by spectral comparison to the analogous derivatives of ferric P450cam. The stability of these complexes and their similarity to the analogous forms of P450cam illustrates the potential of the H175C/D235L CcP double mutant as a model for ferric P450 enzymes. This study marks the first time a stable cyanoferric complex of a model P450 has been made and demonstrates the importance of the environment around the primary coordination ligands in stabilizing metal-ligand ligation.     

Human trifunctional protein deficiency: a new disorder of mitochondrial fatty acid beta-oxidation.

In this paper we report the identification of a new disorder of mitochondrial fatty acid beta-oxidation in a patient which presented with clear manifestations of a mitochondrial beta-oxidation disorder. Subsequent studies in fibroblasts revealed an impairment in palmitate beta-oxidation and in addition, a combined deficiency of long-chain enoyl-CoA hydratase, long-chain 3-hydroxyacyl-CoA-dehydrogenase and long-chain 3-oxoacyl-CoA thiolase. The recent identification of a multifunctional, membrane-bound beta-oxidation enzyme protein catalyzing all these three enzyme activities (Carpenter et al. (1992) Biochem. Biophys. Res. Commun. 183, 443-448; Uchida et al. (1992) J. Biol. Chem. 267, 1034-1041) suggested an underlying basis for this peculiar combination of three enzyme deficiencies. We show by means of size-exclusion chromatography that there is, indeed, a deficiency of the multifunctional beta-oxidation enzyme protein in this patient.     

Information transfer in metabolic pathways. Effects of irreversible steps in computer models.

Various metabolic models have been studied by computer simulation in an effort to understand why allowing for the reversibility of the reaction catalysed by pyruvate kinase, normally considered as irreversible for all practical purposes, significantly altered the behaviour of the model of glycolysis in Trypanosoma brucei [Eisenthal, R. & Cornish-Bowden, A. (1998) J. Biol. Chem. 273, 5500-5505]. Studies of several much simpler models indicate that the enzymes catalysing early steps in a pathway must receive information about the concentrations of the metabolites at the end of the pathway if a model is to be able to reach a steady state; treating all internal steps as reversible is just one way of ensuring this. Feedback inhibition provides a much better way, and as long as feedback loops are present in a model it makes almost no difference to the behaviour whether the intermediate steps with large equilibrium constants are treated as irreversible. In the absence of feedback loops, ordinary product inhibition of all the enzymes in the chain can also transfer information; this is efficient for regulating fluxes but very inefficient for regulating intermediate concentrations. More complicated patterns of regulation, such as activation of a competing branch or forcing flux through a parallel route, can also serve to some degree as ways of passing information around an irreversible step. However, they normally do so less efficiently than inhibition, because the extent to which an enzyme or a pathway can be activated always has an upper limit (which may be below what is required), whereas most enzymes are inhibited completely at saturating concentrations of inhibitor.     

Lipoamide dehydrogenase from Azotobacter vinelandii: site-directed mutagenesis of the His450-Glu455 diad. Spectral properties of wild type and mutated enzymes.

Three amino acid residues in the active site of lipoamide dehydrogenase from Azotobacter vinelandii were replaced by other residues. His450, the active-site base, was changed into Ser, Tyr and Phe. Pro451, in cis conformation, was changed into Ala. Glu455 was replaced with Asp and Gln. Absorption, fluorescence and CD spectroscopy of the mutated enzymes in their oxidized state (Eox) showed only minor changes with respect to the wild-type enzyme, whereas considerable changes were observed in the spectra of the two-electron-reduced (EH2) species of the enzymes upon reduction by the substrate dihydrolipoamide. Differences in extent of reduction of the flavin by NADH indicate that the redox potential of the flavin is altered by the mutations. Enzyme Pro451----Ala [corrected] showed the greatest deviation from wild type. The enzyme is very easily over-reduced to the four-electron reduced state (EH4) by dihydrolipoamide. This is probably due to a change in the backbone conformation caused by the cis-trans conversion. From studies on the pH dependence of the thiolate charge-transfer absorption and the relative fluorescence of EH2 of the enzymes, it is concluded that mutation of His450 results in a relatively simple and easily interpreted distribution of electronic species at the EH2 level. For all three His450-mutated enzymes an apparent pKa1 near 5.5 is calculated that is assigned to the interchange thiol. A second apparent pKa2 is calculated of 6.9, 7.5 and 7.1 for the His450----Phe, -Ser and -Tyr enzymes, respectively, and signifies the deprotonation of the tautomeric equilibrium between the interchange and charge-transfer thiols. The difference in apparent pKa2 values between the His450-mutated enzymes is explained by changes in micropolarity. At the EH2 level the wild-type enzyme consists of multiple electronic forms as reported for the Escherichia coli enzyme [Wilkinson, K. D. and Williams C. H. Jr (1979) J. Biol. Chem. 254, 852-862]. Based on the results obtained with the His450-mutated enzymes, it is concluded that the lowest pKa is associated with the interchange thiol. A model for the equilibrium species of the wild-type enzyme at the EH2 level is presented which takes three pKa values into account. The results of the pH dependence of the electronic species at the EH2 level of Glu455-mutated enzymes essentially follow the model proposed for the wild-type enzyme. However mutation of Glu455 shifts the tautomeric equilibrium of EH2 in favor of the charge-transfer species.(ABSTRACT TRUNCATED AT 400 WORDS)     

Effect of species differences on stromelysin-1 (MMP-3) inhibitor potency. An explanation of inhibitor selectivity using homology modeling and chimeric proteins.

For an animal model to predict a compound's potential for treating human disease, inhibitor interactions with the cognate enzymes of separate species must be comparable. Rabbit and human isoforms of stromelysin-1 are highly homologous, yet there are clear and significant compound-specific differences in inhibitor potencies between these two enzymes. Using crystal structures of discordant inhibitors complexed with the human enzyme, we generated a rabbit enzyme homology model that was used to identify two unmatched residues near the active site that could explain the observed disparities. To test these observations, we designed and synthesized three chimeric mutants of the human enzyme containing the single (H224N and L226F) and double (H224N/L226F) mutations. A comparison of inhibitor potencies among the mutant and wild-type enzymes shows that the mutation of a single amino acid in the human enzyme, histidine 224 to asparagine, is sufficient to change the selectivity profile of the mutant to that of the rabbit isoform. These studies emphasize the importance of considering species differences, which can result from even minor protein sequence variations, for the critical enzymes in an animal disease model. Homology modeling provides a tool to identify key differences in isoforms that can significantly affect native enzyme activity.     

Systems analysis of the tricarboxylic acid cycle in Dictyostelium discoideum. II. Control analysis.

A steady-state computer model of the tricarboxylic cycle in Dictyostelium discoideum was analyzed using metabolic control theory. The steady state had variations of less than 0.04% over the last half of the simulation for both metabolite concentrations and fluxes. Metabolite and flux control coefficients were determined by varying enzymatic activities within 2% of their initial values and simulating the responses of metabolite concentrations and fluxes to these changes. Under these conditions, summation properties were met for most metabolite and all flux control coefficients. Maximum flux control coefficients were found for succinate dehydrogenase (0.35), malic enzyme (0.24), and malate dehydrogenase (-0.18). Comparable control was found for the reaction supplying pyruvate (0.14) and for the sum of the input amino acids (0.43), which serve as an energy source for D. discoideum. The time-dependent processes by which a new steady state was established were examined after increasing malic enzyme or malate dehydrogenase activities. This provided a method for an analysis of the mechanisms by which the observed control coefficients were generated. In addition, the effects of increasing the stimuli within 5-20% of the original enzyme activity were examined. Under these conditions, more typical of experimental stimuli and measurable responses, the metabolic model failed to return to steady state, and thus summation properties were not met. Whether "true" steady states ever occur or whether metabolic control theory can be applied in vivo is discussed.     

Fluoride-inhibited calcium ATPase of sarcoplasmic reticulum. Magnesium and fluoride stoichiometry.

The sarcoplasmic reticulum (SR) CaATPase is inactivated by fluoride in the presence of magnesium (Murphy, A. J., and Coll, R. J. (1992) J. Biol. Chem. 267, 5229-5235). The inactive complex is very stable and can be isolated free of other components by 48 h of dialysis at 4 degrees C (Murphy, A. J., and Coll, R. J. (1992) J. Biol. Chem. 267, 16990-16994). In this study, we used a fluoride-specific electrode to determine that the amount of tightly bound fluoride in the complex was 9.4 +/- 2 nmol mg-1 SR protein. The rate constant of inactivation was very similar to the rate constant of fluoride incorporation and varied directly as the square of the fluoride concentration. Luminal Ca2+ accelerated reactivation of the inhibited enzyme, and the rate constants of activity regain and fluoride release were very similar. Although required for inhibition, added magnesium did not accelerate reactivation. Analysis for magnesium using antipyrylazo III of the inhibited enzyme showed 4.1 +/- 0.4 nmol mg-1 SR protein. As there is much evidence in the literature supportive of an estimate of calcium pumps equal to approximately 4-5 nmol mg-1 SR protein, our results indicate that each inhibited enzyme contains two tightly bound fluorides and one tightly bound magnesium.     

Steady-state kinetics of models of respiratory chain enzymes with isopotential pools and conformational site enzymes.

The formal kinetics of the respiratory chain enzymes is a very complicated mathematical problem. Our main objective here is to introduce methodology that should be useful in approaching this problem, especially at steady state. Most of the work reported here is on models that involve one or two isopotential pools of enzymes and also possibly a site enzyme between two pools that may or may not undergo a conformation change. Preliminary topics treated are the “cross-over” phenomenon and the “mass action” approximation.     

Role of the 2-amino group of deoxyguanosine in sequence recognition by EcoRI restriction and modification enzymes.

The dG residues within the EcoRI recognition sequence of ColE1 DNA have been selectively replaced with dI. Methylation of the altered sequence by the EcoRI modification enzyme is extremely slow as compared with methyl transfer to the natural recognition site. Since the affinity of the modification enzyme for the dI-containing sequence is considerably less than that for the natural sequence, we have concluded that the 2-amino group of dG has an important role in DNA site recognition by this enzyme. In contrast, the altered site is subject to cleavage by EcoRI endonuclease at rates essentially identical with those observed with the natural sequence. These results strongly suggest that the two enzymes utilize different contacts within the EcoRI site and are consisted with our conclusion (Rubin, R. A., and Modrich, P. (1977) J. Biol. Chem. 252, 7265-7272) that the two proteins interact with their common recognition sequence in different ways.     

Comparison of fungal glucose oxidases. Chemical, physicochemical and immunological studies.

Glucose oxidases (beta-D-glucose:oxygen 1-oxidoreductase, EC 1.1.3.4) from two fungal genera (Aspergillus and Penicillium) were studied chemically, physicochemically and immunologically to elucidate the similarities and dissimilarities between these enzymes. Investigation of circular dichroism spectra revealed that these enzymes proteins possess essentially identical conformations. However, differences found in thermal inactivation parameters, catalytic parameters and quantitative immunological reactivities indicate that these enzymes must have some minor but distinct variations in their structures. Interestingly, it was observed that the Penicillium enzyme cross-reacted with the antiserum against the Aspergillus enzyme with an association constant of two orders of magnitude lower than that of the Aspergillus enzyme, and that the precipitin one of the Penicillium enzyme fused together with that of the Aspergillus enzyme in the immunodouble diffusion test. These results lead to the conclusion that these enzymes are closely related but not completely identical, and suggest that they might have evolved from a common ancestral precursor.     

Identification of two discrete peptide: N-glycanases in Oryzias latipes during embryogenesis.

Two different types of peptide:N-glycanase (PNGase) were identified in developing embryos of medaka fish ( Oryzias latipes ). Because the optimum pH values for their activities were acidic and neutral, they were designated as acid PNGase M and neutral PNGase M, respectively. The acid PNGase M corresponded to the enzyme that had been partially purified from medaka embryos (Seko,A., Kitajima,K., Inoue,Y. and Inoue,S. (1991) J. Biol. Chem., 266, 22110-22114). The apparent molecular weight of this enzyme was 150 K, and the optimal pH was 3.5-4.0, and the K m for L-hyosophorin was 44 microM. L-Hyosophorin is a cortical alveolus-derived glycononapeptide with a large N-linked glycan chain present in the perivitelline space of the developing embryo. Acid PNGase M was competitively inhibited by a free de-N-glycosylated nonapeptide derived from L-hyosophorin. This enzyme was expressed in ovaries and embryos at all developmental stages after gastrulation, but activity was not detected in embryos at developmental stages between fertilization and gastrula. Several independent lines of evidence suggested that acid PNGase M may be responsible for the unusual accumulation of free N-glycans derived from yolk glycoproteins (Iwasaki,M., Seko,A., Kitajima,K., Inoue,Y. and Inoue,S. (1992) J. Biol. Chem., 267, 24287-24296). In contrast, the neutral PNGase M was expressed in blastoderms from the 4-8 cell stage and in cells up to early gastrula. The general significance of these findings is that they show a developmental stage-dependent expression of the two PNGase activities, and that expression of the neutral PNGase M activity occurs concomitantly with the de-N-glycosylation of L-hyosophorin. These data thus support our conclusion that the neutral PNGase M is responsible for the developmental-stage-related de-N-glycosylation of the L-hyosophorin.     

Effect of Carbon Source and the Role of Cyclic Adenosine 3′,5′-Monophosphate on the Caulobacter Cell Cycle

The expression of cell cycle events in Caulobacter crescentus CB13 has been shown to be associated with regulation of carbohydrate utilization. Growth on lactose and galactose depends on induction of specific enzymes. Prior growth on glucose results in a delay in enzyme expression and cell cycle arrest at the nonmotile, predivisional stage. Dibutyryl cyclic adenosine 3',5'-monophosphate (AMP) was shown to stimulate expression of the inducible enzymes and, thus, the initiation of the cell cycle. beta-Galactosidase-constitutive mutants did not exhibit a cell cycle arrest upon transfer of cultures from glucose to lactose. Furthermore, carbon source starvation results in accumulation of the cells at the predivisional stage. The cell cycle arrest therefore results from nutritional deprivation and is analogous to the general control system exhibited by yeast (Hartwell, Bacteriol. Rev. 38:164-198, 1974; Wolfner et al., J. Mol. Biol. 96:273-290, 1975), which coordinates cell cycle initiation with metabolic state. Transfer of C. crescentus CB13 from glucose to mannose did not result in a cell cycle arrest, and it was demonstrated that this carbon source is metabolized by constitutive enzymes. Growth on mannose, however, is stimulated by exogenous dibutyryl cyclic AMP without a concomitant increase in the specific activity of the mannose catabolic enzymes. The effect of cyclic AMP on growth on sugars metabolized by inducible enzymes, as well as on sugars metabolized by constitutive enzymes, may represent a regulatory system common to both types of sugar utilization, since they share features that differ from glucose utilization, namely, temperature-sensitive growth and low intracellular concentrations of cyclic guanosine 3',5'-monophosphate.     

Comparative studies of two cathepsin B isozymes from porcine spleen. Isolation, polypeptide chain arrangements, and enzyme specificity

Our previous studies on carbohydrate structures of purified porcine spleen cathepsin B indicated that there are two cathepsin B isozymes, each containing a different carbohydrate (Takahashi, T., Schmidt, P.G., and Tang, J. (1984) J. Biol. Chem. 259, 6059-6062). We have now isolated these two enzymes and carried out a comparative study on their structures and enzymic properties. The major isozyme (CB-I) is a two-chain enzyme (Mr = 28,000) with a light chain (Mr = 5,000) and a heavy chain (Mr = 23,000), whereas the minor enzyme (CB-II) is a single chain enzyme (Mr = 27,000). The NH2-terminal amino acid residues of CB-I were leucine and valine for the light and heavy chain, respectively. However, the NH2-terminal residue of CB-II was not available for automated Edman degradation. In addition, peptide mapping experiments indicated a difference in the primary structure of these two proteins. Despite such structural differences, they are similar in many enzymic properties. CB-I was more catalytically efficient than CB-II toward synthetic substrates, except for the substrate benzoyl-L-arginine beta-naphthylamide for which the relative catalytic efficiency is reversed. Both isozymes degraded glucagon by a dipeptidyl carboxypeptidase activity. Under the same conditions, CB-I was 4-5 times more efficient than CB-II. The results indicate that the cathepsin B isozymes are two separate gene products, but they are similar in enzymic properties.     

Computational functions in biochemical reaction networks.

In prior work we demonstrated the implementation of logic gates, sequential computers (universal Turing machines), and parallel computers by means of the kinetics of chemical reaction mechanisms. In the present article we develop this subject further by first investigating the computational properties of several enzymatic (single and multiple) reaction mechanisms: we show their steady states are analogous to either Boolean or fuzzy logic gates. Nearly perfect digital function is obtained only in the regime in which the enzymes are saturated with their substrates. With these enzymatic gates, we construct combinational chemical networks that execute a given truth-table. The dynamic range of a network's output is strongly affected by "input/output matching" conditions among the internal gate elements. We find a simple mechanism, similar to the interconversion of fructose-6-phosphate between its two bisphosphate forms (fructose-1,6-bisphosphate and fructose-2,6-bisphosphate), that functions analogously to an AND gate. When the simple model is supplanted with one in which the enzyme rate laws are derived from experimental data, the steady state of the mechanism functions as an asymmetric fuzzy aggregation operator with properties akin to a fuzzy AND gate. The qualitative behavior of the mechanism does not change when situated within a large model of glycolysis/gluconeogenesis and the TCA cycle. The mechanism, in this case, switches the pathway's mode from glycolysis to gluconeogenesis in response to chemical signals of low blood glucose (cAMP) and abundant fuel for the TCA cycle (acetyl coenzyme A).     

Competitive inhibition of lipolytic enzymes. I. A kinetic model applicable to water-insoluble competitive inhibitors

It is now becoming clear from the abundant lipolytic enzyme literature that any meaningful interpretation of inhibition data has to take into account the kinetics of enzyme action at the lipid/water interface. We attempt in the present paper to provide a kinetic model applicable to water-insoluble competitive inhibitors, in order to quantitatively compare the results obtained at several laboratories. We derived kinetic equations applicable to the pre-steady state as well as steady state. By measuring the inhibitory power, as described in the present paper, it is possible to obtain a normalized estimation of the relative efficiency of various potential inhibitors. Furthermore, with the kinetic treatment developed here, it is possible to make quantitative comparisons with the same inhibitor placed under various physico-chemical situations, i.e., micellar or monolayer states.     

Aminoethylcysteine can replace the function of the essential active site lysine of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl coenzyme A reductase.

The biodegradative 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase of Pseudomonas mevalonii catalyzes the NAD(+)-dependent conversion of (S)-HMG-CoA to (R)-mevalonate. Crystallographic analysis of abortive ternary complexes of this enzyme revealed lysine 267 located at a position in the active site, suggesting that it might serve as the general acid/base for catalysis. Site-directed mutagenesis and subsequent chemical derivatization were therefore employed to investigate this active site lysine. Replacement of lysine 267 by alanine, histidine, or arginine resulted in mutant enzymes that lacked detectable activity. Lysine 267 was next replaced by the lysine analogues aminoethylcysteine and carboxyamidomethylcysteine. Using instead of the wild-type enzyme the fully active, cysteine-free mutant enzyme C156A/C296A, lysine 267 was first replaced by cysteine. Cysteine 267 of mutant enzyme C156A/C296A/K267C was then treated with bromoethylamine or iodoacetamide to insert aminoethylcysteine (AEC) or carboxyamidomethylcysteine at position 267. The carboxyamidomethylcysteine derivative was inactive, whereas mutant enzyme C156A/C296A/K267AEC exhibited high catalytic activity. That aminoethylcysteine, but not other basic amino acids, can replace the function of lysine 267 documents both the importance of this residue and the requirement for a precisely positioned positive charge at the active site of the enzyme.     

The dynamics of single-substrate continuous cultures: the role of ribosomes

When a chemostat is perturbed from its steady state, it displays complex dynamics. For instance, if the identity of the growth-limiting substrate is switched abruptly, the substrate concentration and cell density undergo a pronounced excursion from the steady state that can last several days. These dynamics occur because certain physiological variables respond slowly. In the literature, several physiological variables have been postulated as potential sources of the slow response. These include transport enzymes, biosynthetic enzymes, and ribosomes. We have been addressing this problem by systematically exploring the role of these variables. In previous work Shoemaker et al. (J. Theor. Biol., 222 (2003) 307-322), we studied the role of transport enzymes, and we showed that transients starting from low transport enzyme levels could be quantitatively captured by a model taking due account of transport enzyme synthesis. However, there is some experimental data indicating that slow responses occur even if the initial enzyme levels are high. Here, we analyse this data to show that in these cases, the sluggish response is most probably due to slow adjustment of the ribosome levels. To test this hypothesis, we extend our previous model by accounting for the evolution of both the transport enzyme and the ribosomes. Based on a kinetic analysis of the data in the literature, we assume that the specific protein synthesis rate is proportional to the ribosome level, and the specific ribosome synthesis rate is autocatalytic. Simulations of the model show remarkable agreement with experimentally observed steady states and the transients. Specifically, the model predictions are in good agreement with (1) the steady-state profiles of the cell density, substrate concentration, RNA, proteins, and transport enzymes, (2) the instantaneous specific substrate uptake, growth, and respiration rates in response to a continuous-to-batch shift, and (3) the transient profiles of the cell density, substrate concentration, and RNA in response to feed switches and dilution rate shifts. Time-scale analysis of the model reveals that every transient response is a combination of two fundamental (and simpler) dynamics, namely, substrate-sufficient batch dynamics and cell-sufficient fed-batch dynamics. We obtain further insight into the transient response by analysing the equations describing these fundamental dynamics. The analysis reveals that in feed switches or dilution rate shift-ups, the transport enzyme reaches a maximum before RNA achieves its maximum, and in dilution rate shift-downs the cell density reaches a maximum before RNA achieves a minimum.