A General Framework for Thermodynamically Consistent Parameterization and Efficient Sampling of Enzymatic Reactions

Kinetic models provide the means to understand and predict the dynamic behaviour of enzymes upon different perturbations. Despite their obvious advantages, classical parameterizations require large amounts of data to fit their parameters. Particularly, enzymes displaying complex reaction and regulatory (allosteric) mechanisms require a great number of parameters and are therefore often represented by approximate formulae, thereby facilitating the fitting but ignoring many real kinetic behaviours. Here, we show that full exploration of the plausible kinetic space for any enzyme can be achieved using sampling strategies provided a thermodynamically feasible parameterization is used. To this end, we developed a General Reaction Assembly and Sampling Platform (GRASP) capable of consistently parameterizing and sampling accurate kinetic models using minimal reference data. The former integrates the generalized MWC model and the elementary reaction formalism. By formulating the appropriate thermodynamic constraints, our framework enables parameterization of any oligomeric enzyme kinetics without sacrificing complexity or using simplifying assumptions. This thermodynamically safe parameterization relies on the definition of a reference state upon which feasible parameter sets can be efficiently sampled. Uniform sampling of the kinetics space enabled dissecting enzyme catalysis and revealing the impact of thermodynamics on reaction kinetics. Our analysis distinguished three reaction elasticity regions for common biochemical reactions: a steep linear region (0> ΔG_r >-2 kJ/mol), a transition region (-2> ΔG_r >-20 kJ/mol) and a constant elasticity region (ΔG_r <-20 kJ/mol). We also applied this framework to model more complex kinetic behaviours such as the monomeric cooperativity of the mammalian glucokinase and the ultrasensitive response of the phosphoenolpyruvate carboxylase of Escherichia coli. In both cases, our approach described appropriately not only the kinetic behaviour of these enzymes, but it also provided insights about the particular features underpinning the observed kinetics. Overall, this framework will enable systematic parameterization and sampling of enzymatic reactions.     

Thermodynamic and kinetic characterisation of individual haems in multicentre cytochromes c3

The characterisation of individual centres in multihaem proteins is difficult due to the similarities in the redox and spectroscopic properties of the centres. NMR has been used successfully to distinguish redox centres and allow the determination of the microscopic thermodynamic parameters in several multihaem cytochromes c(3) isolated from different sulphate-reducing bacteria. In this article we show that it is also possible to discriminate the kinetic properties of individual centres in multihaem proteins, if the complete microscopic thermodynamic characterisation is available and the system displays fast intramolecular equilibration in the time scale of the kinetic experiment. The deconvolution of the kinetic traces using a model of thermodynamic control provides a reference rate constant for each haem that does not depend on driving force and can be related to structural factors. The thermodynamic characterisation of three tetrahaem cytochromes and their kinetics of reduction by sodium dithionite are reported in this paper. Thermodynamic and kinetic data were fitted simultaneously to a model to obtain microscopic reduction potentials, haem-haem and haem-proton interacting potentials, and reference rate constants for the haems. The kinetic information obtained for these cytochromes and recently published data for other multihaem cytochromes is discussed with respect to the structural factors that determine the reference rates. The accessibility for the reducing agent seems to play an important role in controlling the kinetic rates, although is clearly not the only factor.     

Thermodynamic and kinetic characterisation of individual haems in multicentre cytochromes c3

The characterisation of individual centres in multihaem proteins is difficult due to the similarities in the redox and spectroscopic properties of the centres. NMR has been used successfully to distinguish redox centres and allow the determination of the microscopic thermodynamic parameters in several multihaem cytochromes c₃ isolated from different sulphate-reducing bacteria. In this article we show that it is also possible to discriminate the kinetic properties of individual centres in multihaem proteins, if the complete microscopic thermodynamic characterisation is available and the system displays fast intramolecular equilibration in the time scale of the kinetic experiment. The deconvolution of the kinetic traces using a model of thermodynamic control provides a reference rate constant for each haem that does not depend on driving force and can be related to structural factors. The thermodynamic characterisation of three tetrahaem cytochromes and their kinetics of reduction by sodium dithionite are reported in this paper. Thermodynamic and kinetic data were fitted simultaneously to a model to obtain microscopic reduction potentials, haem-haem and haem-proton interacting potentials, and reference rate constants for the haems. The kinetic information obtained for these cytochromes and recently published data for other multihaem cytochromes is discussed with respect to the structural factors that determine the reference rates. The accessibility for the reducing agent seems to play an important role in controlling the kinetic rates, although is clearly not the only factor.     

Chromatophore heterogeneity explains phenomena seen in Rhodobacter sphaeroides previously attributed to supercomplexes

It is generally considered that metabolic reactions are well described by homogeneous kinetic models in which the reaction phase is statistically uniform. In membranes, especially in photosynthetic systems where the protein complement is high, it has recently been recognized that effects of local heterogeneity might contribute additional factors that perturb the kinetic behavior, and require more extensive treatment. We show in this paper that statistical heterogeneity in vesicular systems can also contribute to quite marked discrepancies from the behavior expected from standard kinetic and thermodynamic models, for reactions involving free diffusion in the aqueous phase. We explain the kinetic and thermodynamic effects observed in studies of photosynthetic electron transfer in cells and chromatophores from Rhodobacter sphaeroides previously attributed to supercomplexes, in terms of a model based on heterogeneity in distribution of electron transfer components among the chromatophore population. We discuss examples of data inconsistent with the supercomplex model, but well explained by the heterogeneity model.     

Engineering proteins with tunable thermodynamic and kinetic stabilities

It is widely recognized that enhancement of protein stability is an important biotechnological goal. However, some applications at least, could actually benefit from stability being strongly dependent on a suitable environment variable, in such a way that enhanced stability or decreased stability could be realized as required. In therapeutic applications, for instance, a long shelf-life under storage conditions may be convenient, but a sufficiently fast degradation of the protein after it has performed the planned molecular task in vivo may avoid side effects and toxicity. Undesirable effects associated to high stability are also likely to occur in food-industry applications. Clearly, one fundamental factor involved here is the kinetic stability of the protein, which relates to the time-scale of the irreversible denaturation processes and which is determined to some significant extent by the free-energy barrier for unfolding (the barrier that "separates" the native state from the highly-susceptible-to-irreversible-alterations nonnative states). With an appropriate experimental model, we show that strong environment-dependencies of the thermodynamic and kinetic stabilities can be achieved using robust protein engineering. We use sequence-alignment analysis and simple computational electrostatics to design stabilizing and destabilizing mutations, the latter introducing interactions between like charges which are screened out at high salt. Our design procedures lead naturally to mutating regions which are mostly unstructured in the transition state for unfolding. As a result, the large salt effect on the thermodynamic stability of our consensus plus charge-reversal variant translates into dramatic changes in the time-scale associated to the unfolding barrier: from the order of years at high salt to the order of days at low salt. Certainly, large changes in salt concentration are not expected to occur in biological systems in vivo. Hence, proteins with strong salt-dependencies of the thermodynamic and kinetic stabilities are more likely to be of use in those cases in which high-stability is required only under storage conditions. A plausible scenario is that inclusion of high salt in liquid formulations will contribute to a long protein shelf-life, while the lower salt concentration under the conditions of the application will help prevent the side effects associated with high-stability which may potentially arise in some therapeutic and food-industry applications. From a more general viewpoint, this work shows that consensus engineering and electrostatic engineering can be readily combined and clarifies relevant aspects of the relation between thermodynamic stability and kinetic stability in proteins.     

Experimental concerns in tracer kinetic investigations of apolipoprotein metabolism

The complexity of the metabolism of the plasma lipoproteins makes it impossible to integrate the details of the reactions of specific apolipoproteins and their associated lipids without the use of computerized modeling methods. Because apolipoproteins impart specificity in the transport and chemical processing of plasma lipids, they have been the focus of many in vivo kinetic tracer investigations. The analysis of such kinetic data by modeling techniques has provided important advances in understanding lipoprotein metabolism. An example is the Delipidation Chain, an hypothesis explaining VLDL metabolism in terms of a sequential delipidation process. As a consequence of the advance in knowledge of apolipoprotein structure and metabolism, coupled with progress in computerized modeling of large systems, it has become important to refine the design of in vivo tracer kinetic investigations of the apolipoproteins. Considerations of particular importance include the selection of apolipoprotein tracers which can be shown to undergo the same reactions as the apolipoproteins whose metabolism they trace. If the physical and chemical processes which convert apolipoproteins from one metabolic pool to another are to be analyzed correctly, it is necessary to describe precisely and to measure accurately these pools. Current methods for delineating metabolic pools of apolipoproteins in vivo need to be refined. When accomplished, this will provide new opportunities to investigate the metabolic pathways of the apolipoproteins and their associated lipids. A very important challenge is to design experiments which will differentiate transfer processes, which result in net transport of a reactant, from exchange processes, whereby a tracer and a tracee are exchanged between pools without a net transport event occuring. Since both types of processes occur readily with apolipoproteins, it is important to develop methods to examine them separately. Computerized kinetic modeling provides a means for describing and understanding the complexities of lipoprotein metabolism. A major challenge is for the experimentalist to acquire data which accurately reflect the physiological processes involved in lipoprotein metabolism.     

Unifying temperature effects on the growth rate of bacteria and the stability of globular proteins

The specific growth rate constant for bacterial growth does not obey the Arrhenius-type kinetics displayed by simple chemical reactions. Instead, for bacteria, steep convex curves are observed on an Arrhenius plot at the low- and high-temperature ends of the biokinetic range, with a region towards the middle of the growth range loosely approximating linearity. This central region has been considered by microbiologists to be the "normal physiological range" for bacterial growth, a concept whose meaningfulness we now question. We employ a kinetic model incorporating thermodynamic terms for temperature-induced enzyme denaturation, central to which is a term to account for the large positive heat capacity change during unfolding of the proteins within the bacteria. It is now widely believed by biophysicists that denaturation of complex proteins and/or other macromolecules is due to hydrophobic hydration of non-polar compounds. Denaturation is seen as the process by which enthalpic and entropic forces becomes imbalanced both at high and at low temperatures resulting in conformational changes in the enzyme structure that expose hydrophobic amino acid groups to the surrounding water molecules. The "thermodynamic" rate model, incorporating the heat capacity change and its effect on the enthalpy and entropy of the system, fitted 35 sets of data for psychrophilic, psychrotrophic, mesophilic and thermophilic bacteria well, resulting in biologically meaningful estimates for the important thermodynamic parameters. As these results mirror those obtained by biophysicists for globular proteins, it appears that the same or a similar mechanism applies to bacteria as applies to proteins.     

Kinetic and thermodynamic characterization of the RNA guanylyltransferase reaction

An RNA guanylyltransferase activity is involved in the synthesis of the cap structure found at the 5' end of eukaryotic mRNAs. The RNA guanylyltransferase activity is a two-step ping-pong reaction in which the enzyme first reacts with GTP to produce the enzyme-GMP covalent intermediate with the concomitant release of pyrophosphate. In the second step of the reaction, the GMP moiety is then transferred to a diphosphorylated RNA. Both reactions were previously shown to be reversible. In this study, we report a biochemical and thermodynamic characterization of both steps of the reaction of the RNA guanylyltransferase from Paramecium bursaria Chlorella virus 1, the prototype of a family of viruses infecting green algae. Using a combination of real-time fluorescence spectroscopy, radioactive kinetic assays, and inhibition assays, the complete kinetic parameters of the RNA guanylyltransferase were determined. We produced a thermodynamic scheme for the progress of the reaction as a function of the energies involved in each step. We were able to demonstrate that the second step comprises the limiting steps for both the direct and reverse overall reactions. In both cases, the binding to the RNA substrates is the step requiring the highest energy and generating unstable intermediates that will promote the catalytic activites of the enzyme. This study reports the first thorough kinetic and thermodynamic characterization of the reaction catalyzed by an RNA capping enzyme.     

A probabilistic framework for the exploration of enzymatic capabilities based on feasible kinetics and control analysis

BackgroundAnalysis of limiting steps within enzyme-catalyzed reactions is fundamental to understand their behavior and regulation. Methods capable of unravelling control properties and exploring kinetic capabilities of enzymatic reactions would be particularly useful for protein and metabolic engineering. While single-enzyme control analysis formalism has previously been applied to well-studied enzymatic mechanisms, broader application of this formalism is limited in practice by the limited amount of kinetic data and the difficulty of describing complex allosteric mechanisms.MethodsTo overcome these limitations, we present here a probabilistic framework enabling control analysis of previously unexplored mechanisms under uncertainty. By combining a thermodynamically consistent parameterization with an efficient Sequential Monte Carlo sampler embedded in a Bayesian setting, this framework yields insights into the capabilities of enzyme-catalyzed reactions with modest kinetic information, provided that the catalytic mechanism and a thermodynamic reference point are defined.ResultsThe framework was used to unravel the impact of thermodynamic affinity, substrate saturation levels and effector concentrations on the flux control and response coefficients of a diverse set of enzymatic reactions.ConclusionsOur results highlight the importance of the metabolic context in the control analysis of isolated enzymes as well as the use of statistically sound methods for their interpretation.General SignificanceThis framework significantly expands our current capabilities for unravelling the control properties of general reaction kinetics with limited amount of information. This framework will be useful for both theoreticians and experimentalists in the field.     

Mathematical modelling of glycolysis and adenine nucleotide metabolism of human erythrocytes. I. Reaction-kinetic statements, analysis of in vivo state and determination of starting conditions for in vitro experiments

A mathematical model of energy metabolism of human red cells is presented, which includes besides the glycolytic reactions the adenine nucleotide metabolism. The model is based on the network of chemical reactions, the thermodynamic equilibrium constants of fast reversible reactions and on the kinetic equations for irreversible enzyme reactions. The model consists of a system of 16 differential equations and allows the mathematical evaluation of metabolic levels in the steady state of energy metabolism corresponding to the in vivo state erythrocytes with the kinetic data for the enzymes derived from in vitro experiments. The dependence of the levels of metabolites in the steady state on the activity of some enzymes is analysed to characterize the regulatory properties of the system. The comparison of the steady state levels of the model with experimental data makes it possible to estimate values of some controversial enzyme parameters. Estimates of the kinetic parameters of the following intracellular processes are presented: 1) rate constant of AMP-phosphatase, 2) maximum rate of adenylate deaminase, 3) activity of adenine phosphoribosylpyrophosphate transferase and 4) adenosine transport through the cell membrane. The simulation of the preparatory phase before incubation of erythrocytes indicates, that the model also permits to compute the time course of changes of levels of metabolites. To solve the initial problem the stiff differential equation system is integrated numerically by an efficient program without the application of the quasi-steady-state approximation.     

In vivo kinetics of primary metabolism in Saccharomyces cerevisiae studied through prolonged chemostat cultivation

In this study, prolonged chemostat cultivation is applied to investigate in vivo enzyme kinetics of Saccharomyces cerevisiae. S. cerevisiae was grown in carbon-limited aerobic chemostats for 70-95 generations, during which multiple steady states were observed, characterized by constant intracellular fluxes but significant changes in intracellular metabolite concentrations and enzyme capacities. We provide evidence for two relevant kinetic mechanisms for sustaining constant fluxes: in vivo near-equilibrium of reversible reactions and tight regulation of irreversible reactions by coordinated changes of metabolic effectors. Using linear-logarithmic kinetics, we illustrate that these multiple steady-state measurements provide linear constraints between elasticity parameters instead of their absolute values. Upon perturbation by a glucose pulse, glucose uptake and ethanol excretion in prolonged cultures were remarkably lower, compared to a reference culture perturbed at 10 generations. Metabolome measurements during the transient indicate that the differences might be due to a reduced ATP regeneration capacity in prolonged cultures.     

Simultaneous determination of thermodynamic and kinetic data by isothermal titration calorimetry

BackgroundThermodynamic and binding kinetic data increasingly support and guide the drug optimization process.MethodsBecause ITC thermograms contain binding thermodynamic and kinetic information, an efficient protocol for the simultaneous extraction of thermodynamic and kinetic data for 1:1 protein ligand reactions from AFFINImeter kinITC in one single experiment are presented.ResultsThe effort to apply this protocol requires the same time as for the standard protocol but increases the precision of both thermodynamic and kinetic data.ConclusionsThe protocol enables reliable extraction of both thermodynamic and kinetic data for 1:1 protein-ligand binding reactions with improved precision compared to the ‘standard protocol’.General significanceThermodynamic and kinetic data are recorded under exactly the same conditions in solution without any labeling or immobilization from a protein sample that is not 100% active and would otherwise render the extraction of kinetic parameters impossible.     

Natural selection for kinetic stability is a likely origin of correlations between mutational effects on protein energetics and frequencies of amino acid occurrences in sequence alignments

It appears plausible that natural selection constrains, to some extent at least, the stability in many natural proteins. If, during protein evolution, stability fluctuates within a comparatively narrow range, then mutations are expected to be fixed with frequencies that reflect mutational effects on stability. Indeed, we recently reported a robust correlation between the effect of 27 conservative mutations on the thermodynamic stability (unfolding free energy) of Escherichia coli thioredoxin and the frequencies of residues occurrences in sequence alignments. We show here that this correlation likely implies a lower limit to thermodynamic stability of only a few kJ/mol below the unfolding free energy of the wild-type (WT) protein. We suggest, therefore, that the correlation does not reflect natural selection of thermodynamic stability by itself, but of some other factor which is linked to thermodynamic stability for the mutations under study. We propose that this other factor is the kinetic stability of thioredoxin in vivo, since( i) kinetic stability relates to irreversible denaturation, (ii) the rate of irreversible denaturation in a crowded cellular environment (or in a harsh extracellular environment) is probably determined by the rate of unfolding, and (iii) the half-life for unfolding changes in an exponential manner with activation free energy and, consequently, comparatively small free energy effects can have deleterious consequences for kinetic stability. This proposal is supported by the results of a kinetic study of the WT form and the 27 single-mutant variants of E. coli thioredoxin based on the global analyses of chevron plots and equilibrium unfolding profiles determined from double-jump unfolding assays. This kinetic study suggests, furthermore, one of the factors that may contribute to the high activation free energy for unfolding in thioredoxin (required for kinetic stability), namely the energetic optimization of native-state residue environments in regions, which become disrupted in the transition state for unfolding.     

Thermodynamic and kinetic characterization of Co-C bond homolysis catalyzed by coenzyme B(12)-dependent methylmalonyl-CoA mutase

Methylmalonyl-CoA mutase is a member of the family of coenzyme B(12)-dependent isomerases and catalyzes the 1,2-rearrangement of methylmalonyl-CoA to succinyl-CoA. A common first step in the reactions catalyzed by coenzyme B(12)-dependent enzymes is cleavage of the cobalt-carbon bond of the cofactor, leading to radical-based rearrangement reactions. Comparison of the homolysis rate for the free and enzyme-bound cofactors reveals an enormous rate enhancement which is on the order of a trillion-fold. To address how this large rate acceleration is achieved, we have examined the kinetic and thermodynamic parameters associated with the homolysis reaction catalyzed by methylmalonyl-CoA mutase. Both the rate and the amount of cob(II)alamin formation have been analyzed as a function of temperature with the protiated substrate. These studies yield the following activation parameters for the homolytic reaction at 37 degrees C: DeltaH(f)() = 18.8 +/- 0.8 kcal/mol, DeltaS(f)() = 18.2 +/- 0.8 cal/(mol.K), and DeltaG(f)() = 13.1 +/- 0.6 kcal/mol. Our results reveal that the enzyme lowers the transition state barrier by 17 kcal/mol, corresponding to a rate acceleration of 0.9 x 10(12)-fold. Both entropic and enthalpic factors contribute to the observed rate acceleration, with the latter predominating. The substrate binding step is exothermic, with a DeltaG of -5.2 kcal/mol at 37 degrees C, and is favored by both entropic and enthalpic factors. We have employed the available kinetic and spectroscopic data to construct a qualitative free energy profile for the methylmalonyl-CoA mutase-catalyzed reaction.     

Modeling of uncertainties in biochemical reactions

Mathematical modeling is an indispensable tool for research and development in biotechnology and bioengineering. The formulation of kinetic models of biochemical networks depends on knowledge of the kinetic properties of the enzymes of the individual reactions. However, kinetic data acquired from experimental observations bring along uncertainties due to various experimental conditions and measurement methods. In this contribution, we propose a novel way to model the uncertainty in the enzyme kinetics and to predict quantitatively the responses of metabolic reactions to the changes in enzyme activities under uncertainty. The proposed methodology accounts explicitly for mechanistic properties of enzymes and physico‐chemical and thermodynamic constraints, and is based on formalism from systems theory and metabolic control analysis. We achieve this by observing that kinetic responses of metabolic reactions depend: (i) on the distribution of the enzymes among their free form and all reactive states; (ii) on the equilibrium displacements of the overall reaction and that of the individual enzymatic steps; and (iii) on the net fluxes through the enzyme. Relying on this observation, we develop a novel, efficient Monte Carlo sampling procedure to generate all states within a metabolic reaction that satisfy imposed constrains. Thus, we derive the statistics of the expected responses of the metabolic reactions to changes in enzyme levels and activities, in the levels of metabolites, and in the values of the kinetic parameters. We present aspects of the proposed framework through an example of the fundamental three‐step reversible enzymatic reaction mechanism. We demonstrate that the equilibrium displacements of the individual enzymatic steps have an important influence on kinetic responses of the enzyme. Furthermore, we derive the conditions that must be satisfied by a reversible three‐step enzymatic reaction operating far away from the equilibrium in order to respond to changes in metabolite levels according to the irreversible Michelis–Menten kinetics. The efficient sampling procedure allows easy, scalable, implementation of this methodology to modeling of large‐scale biochemical networks. Biotechnol. Bioeng. 2011;108: 413–423. © 2010 Wiley Periodicals, Inc.     

Mathematical modeling of precipitation and dissolution reactions in microbiological systems

We expand the biogeochemical model CCBATCH to include a precipitation/dissolution sub-model that contains kinetic and equilibrium options. This advancement extends CCBATCH's usefulness to situations in which microbial reactions cause or are affected by formation or dissolution of a solid phase. The kinetic option employs a rate expression that explicitly includes the intrinsic kinetics for reaction ormass-transport control, the differencefrom thermodynamic equilibrium, and the aqueous concentration of the rate-limiting metal or ligand. The equilibrium feature can be used alone, and it also serves as check that the kinetic rate never is too fast and ``overshoots' equilibrium. The features of the expanded CCBATCH are illustrated by an example in which the precipitation of Fe(OH)_{3 (s)} allows the biodegradation of citric acid, even though complexes are strong and not bioavailable. Precipitation releases citrate ligand, and biodegradation of the citrate increases the pH.     

Application of atomic force microscopy for characteristics of single intermolecular interactions

Atomic force microscopy (AFM) can be used to make measurements in vacuum, air, and water. The method is able to gather information about intermolecular interaction forces at the level of single molecules. This review encompasses experimental and theoretical data on the characterization of ligand-receptor interactions by AFM. The advantage of AFM in comparison with other methods developed for the characterization of single molecular interactions is its ability to estimate not only rupture forces, but also thermodynamic and kinetic parameters of the rupture of a complex. The specific features of force spectroscopy applied to ligand-receptor interactions are examined in this review from the stage of the modification of the substrate and the cantilever up to the processing and interpretation of the data. We show the specificities of the statistical analysis of the array of data based on the results of AFM measurements, and we discuss transformation of data into thermodynamic and kinetic parameters (kinetic dissociation constant, Gibbs free energy, enthalpy, and entropy). Particular attention is paid to the study of polyvalent interactions, where the definition of the constants is hampered due to the complex stoichiometry of the reactions.     

The H2 dissociation on the BN, AlN, BP and AlP nanotubes: a comparative study

The thermodynamic and kinetic feasibility of H(2) dissociation on the BN, AlN, BP and AlP zigzag nanotubes has been investigated theoretically by calculating the dissociation and activation energies. We determined the BN and AlP tubes to be inert toward H(2) dissociation, both thermodynamically and kinetically. The reactions are endothermic by 5.8 and 3 kcal mol(-1), exhibiting high activation energies of 38.8 and 30.6 kcal mol(-1), respectively. Our results indicated that H(2) dissociation is thermodynamically favorable on both PB and AlN nanotubes. However, in spite of the thermodynamic feasibility of H(2) dissociation on PB types, this process is kinetically unfavorable due to partly high activation energy. Generally, we concluded that among the four studied tubes, the AlN nanotube may be an appropriate model for H(2) dissociation process, from a thermodynamic and kinetic stand point. We also indicated that H(2) dissociation is not homolytic, rather it takes place via a heterolytic bond cleavage. In addition, a comparative study has been performed on the electrical and geometrical properties of the tubes. Our analysis showed that the electrical conductivity of tubes is as follows: BP>AlP>BN>AlN depending on how to combine the electron rich and electron poor atoms.     

Determination of the kinetics of guanine nucleotide exchange on EF-Tu and EF-Ts: continuing uncertainties

An analysis is made of the rate constants for the reactions involving the interactions of EF-Tu, EF-Ts, GDP, and GTP recently derived by Gromadski et al. [Biochemistry 41 (2002) 162]. Though their measured values appear to allow a reasonable rate of nucleotide exchange sufficient to support rates of protein synthesis in vivo, their data underestimate the thermodynamic barrier involved in nucleotide exchange and therefore cannot be considered definitive. A kinetic scheme consistent with the thermodynamic barrier can be achieved by modification of various rate constants, particularly of those involving the release of EF-Ts from EF-Tu.GTP.EF-Ts, but such constants are markedly different from what are experimentally observed. It thus remains impossible at present satisfactorily to model guanine nucleotide exchange on EF-Tu, catalysed by EF-Ts by a double displacement mechanism, with experimentally derived rate constants. Metabolic control analysis has been applied to determine the degree of flux control of the different steps in the pathway.