Nonlinear Enzyme Kinetics Can Lead to High Metabolic Flux Control Coefficients: Implications for the Evolution of Dominance

In a classic study, Kacser & Burns (1981, Genetics 97, 639-666) demonstrated that given certain plausible assumptions, the flux in a metabolic pathway was more or less indifferent to the activity of any of the enzymes in the pathway taken singly. It was inferred from this that the observed dominance of most wild-type alleles with respect to loss-of-function mutations did not require an adaptive, meaning selectionist, explanation. Cornish-Bowden (1987, J. theor. Biol. 125, 333-338) showed that the Kacser-Burns inference was not valid when substrate concentrations were large relative to the relevant Michaelis constants. We find that in a randomly constructed functional pathway, even when substrate levels are small, one can expect high values of control coefficients for metabolic flux in the presence of significant nonlinearities as exemplified by enzymes with Hill coefficients ranging from two to six, or by the existence of oscillatory loops. Under these conditions the flux can be quite sensitive to changes in enzyme activity as might be caused by inactivating one of the two alleles in a diploid. Therefore, the phenomenon of dominance cannot be a trivial ''default'' consequence of physiology but must be intimately linked to the manner in which metabolic networks have been moulded by natural selection.

Efficient and accurate experimental design for enzyme kinetics: Bayesian studies reveal a systematic approach

In areas such as drug development, clinical diagnosis and biotechnology research, acquiring details about the kinetic parameters of enzymes is crucial. The correct design of an experiment is critical to collecting data suitable for analysis, modelling and deriving the correct information. As classical design methods are not targeted to the more complex kinetics being frequently studied, attention is needed to estimate parameters of such models with low variance. We demonstrate that a Bayesian approach (the use of prior knowledge) can produce major gains quantifiable in terms of information, productivity and accuracy of each experiment. Developing the use of Bayesian Utility functions, we have used a systematic method to identify the optimum experimental designs for a number of kinetic model data sets. This has enabled the identification of trends between kinetic model types, sets of design rules and the key conclusion that such designs should be based on some prior knowledge of K-M and/or the kinetic model. We suggest an optimal and iterative method for selecting features of the design such as the substrate range, number of measurements and choice of intermediate points. The final design collects data suitable for accurate modelling and analysis and minimises the error in the parameters estimated. (C) 2003 Elsevier Science B.V. All rights reserved.

Donut-Shaped Fingerprint In Homologous Polypeptide Relationships-A Topological Feature Related To Pathogenic Structural Changes In Conformational Disease

Features of homologous relationship of proteins can provide us a general picture of protein universe, assist protein design and analysis, and further our comprehension of the evolution of organisms. Here we carried Out a Study of the evolution Of protein molecules by investigating homologous relationships among residue segments. The motive was to identify detailed topological features of homologous relationships for short residue segments in the whole protein universe. Based on the data of a large number of non-redundant Proteins, the universe of non-membrane polypeptide was analyzed by considering both residue mutations and structural conservation. By connecting homologous segments with edges, we obtained a homologous relationship network of the whole universe of short residue segments, which we named the graph of polypeptide relationships (GPR). Since the network is extremely complicated for topological transitions, to obtain an in-depth understanding, only subgraphs composed of vital nodes of the GPR were analyzed. Such analysis of vital subgraphs of the GPR revealed a donut-shaped fingerprint. Utilization of this topological feature revealed the switch sites (where the beginning of exposure Of previously hidden "hot spots" of fibril-forming happens, in consequence a further opportunity for protein aggregation is Provided; 188-202) of the conformational conversion of the normal alpha-helix-rich prion protein PrPC to the beta-sheet-rich PrPSc that is thought to be responsible for a group of fatal neurodegenerative diseases, transmissible spongiform encephalopathies. Efforts in analyzing other proteins related to various conformational diseases are also introduced. (C) 2009 Elsevier Ltd. All rights reserved.

Stochastic chemical kinetics and the quasi-steady-state assumption : application to the stochastic simulation algorithm and chemical master equation

Recently the application of the quasi-steady-state approximation (QSSA) to the stochastic simulation algorithm (SSA) was suggested for the purpose of speeding up stochastic simulations of chemical systems that involve both relatively fast and slow chemical reactions [Rao and Arkin, J. Chem. Phys. 118, 4999 (2003)] and further work has led to the nested and slow-scale SSA. Improved numerical efficiency is obtained by respecting the vastly different time scales characterizing the system and then by advancing only the slow reactions exactly, based on a suitable approximation to the fast reactions. We considerably extend these works by applying the QSSA to numerical methods for the direct solution of the chemical master equation (CME) and, in particular, to the finite state projection algorithm [Munsky and Khammash, J. Chem. Phys. 124, 044104 (2006)], in conjunction with Krylov methods. In addition, we point out some important connections to the literature on the (deterministic) total QSSA (tQSSA) and place the stochastic analogue of the QSSA within the more general framework of aggregation of Markov processes. We demonstrate the new methods on four examples: Michaelis–Menten enzyme kinetics, double phosphorylation, the Goldbeter–Koshland switch, and the mitogen activated protein kinase cascade. Overall, we report dramatic improvements by applying the tQSSA to the CME solver.

Catalysis of tRNA Aminoacylation: Single Turnover to Steady-State Kinetics of tRNA Synthetases

Aminoacyl-tRNA synthetases (aaRS) catalyze the bimolecular association reaction between amino acid and tRNA by specifically and unerringly choosing the cognate amino acid and tRNA. There are two classes of such synthetases that perform tRNA-aminoacylation reaction. Interestingly, these two classes of aminoacyl-tRNA synthetases differ not only in their structures but they also exhibit remarkably distinct kinetics under pre-steady-state condition. The class I synthetases show initial burst of product formation followed by a slower steady-state rate. This has been argued to represent the influence of slow product release. In contrast, there is no burst in the case of class H enzymes. The tight binding of product with enzyme for class I enzymes is correlated with the enhancement of rate in presence of elongation factor. EF-TU. In spite of extensive experimental studies, there is no detailed theoretical analysis that can provide a quantitative understanding of this important problem. In this article, we present a theoretical investigation of enzyme kinetics for both classes of aminoacyl-tRNA synthetases. We present an augmented kinetic scheme and then employ the methods of time-dependent probability statistics to obtain expressions for the first passage time distribution that gives both the time-dependent and the steady-state rates. The present study quantitatively explains all the above experimental observations. We propose an alternative path way in the case of class II enzymes showing the tRNA-dependent amino acid activation and the discrepancy between the single-turnover and steady-state rate.

A windows program for the derivation of steady-state equations in enzyme systems

Despite the number of computer-assisted methods described for the derivation of steady-state equations of enzyme systems, most of them are focused on strict steady-state conditions or are not able to solve complex reaction mechanisms. Moreover, many of them are based on computer programs that are either not readily available or have limitations. We present here a computer program called WinStes, which derives equations for both strict steady-state systems and those with the assumption of rapid equilibrium, for branched or unbranched mechanisms, containing both reversible and irreversible conversion steps. It solves reaction mechanisms involving up to 255 enzyme species, connected by up to 255 conversion steps. The program provides all the advantages of the Windows programs, such as a user-friendly graphical interface, and has a short computation time. WinStes is available free of charge on request from the authors. (c) 2006 Elsevier Inc. All rights reserved.

Prediction of Intrinsic kinetics in diffusion disguised immobilized enzyme systems

When immobilized enzyme kinetics are disguised by inter- and intraparticle diffusion effects, an approximate mathematical procedure is indicated whereby experimental data obtained in the limiting ranges of first- and zeroth-order Michaelis-Menten kinetics could be used for the prediction of the kinetic constants.

The Mode of Action of Recombinant Mycobacterium tuberculosis Shikimate Kinase: Kinetics and Thermodynamics Analyses

Tuberculosis remains as one of the main cause of mortality worldwide due to a single infectious agent, Mycobacterium tuberculosis. The aroK-encoded M. tuberculosis Shikimate Kinase (MtSK), shown to be essential for survival of bacilli, catalyzes the phosphoryl transfer from ATP to the carbon-3 hydroxyl group of shikimate (SKH), yielding shikimate-3-phosphate and ADP. Here we present purification to homogeneity, and oligomeric state determination of recombinant MtSK. Biochemical and biophysical data suggest that the chemical reaction catalyzed by monomeric MtSK follows a rapid-equilibrium random order of substrate binding, and ordered product release. Isothermal titration calorimetry (ITC) for binding of ligands to MtSK provided thermodynamic signatures of non-covalent interactions to each process. A comparison of steady-state kinetics parameters and equilibrium dissociation constant value determined by ITC showed that ATP binding does not increase the affinity of MtSK for SKH. We suggest that MtSK would more appropriately be described as an aroL-encoded type II shikimate kinase. Our manuscript also gives thermodynamic description of SKH binding to MtSK and data for the number of protons exchanged during this bimolecular interaction. The negative value for the change in constant pressure heat capacity (ΔCp) and molecular homology model building suggest a pronounced contribution of desolvation of non-polar groups upon binary complex formation. Thermodynamic parameters were deconvoluted into hydrophobic and vibrational contributions upon MtSK:SKH binary complex formation. Data for the number of protons exchanged during this bimolecular interaction are interpreted in light of a structural model to try to propose the likely amino acid side chains that are the proton donors to bulk solvent following MtSK:SKH complex formation. © 2013 Rosado et al.

Studies on active site mutants of P-falciparum adenylosuccinate synthetase: Insights into enzyme catalysis and activation

Adenylosuccinate synthetase catalyzes a reversible reaction utilizing IMP, GTP and aspartate in the presence of Mg2+ to form adenylosuccinate, GDP and inorganic phosphate. Comparison of similarly liganded complexes of Plasmodium falciparum, mouse and Escherichia coil AdSS reveals H-bonding interactions involving nonconserved catalytic loop residues (Asn429, Lys62 and Thr307) that are unique to the parasite enzyme. Site-directed mutagenesis has been used to examine the role of these interactions in catalysis and structural organization of P. falciparum adenylosuccinate synthetase (PfAdSS). Mutation of Asn429 to Val, Lys62 to Leu and Thr307 to Val resulted in an increase in K-m values for IMP, GTP and aspartate, respectively along with a 5 fold drop in the k(cat) value for N429V mutant suggesting the role of these residues in ligand binding and/or catalysis. We have earlier shown that the glycolytic intermediate, fructose 1,6 bisphosphate, which is an inhibitor of mammalian AdSS is an activator of the parasite enzyme. Enzyme kinetics along with molecular docking suggests a mechanism for activation wherein F16BP seems to be binding to the Asp loop and inducing a conformation that facilitates aspartate binding to the enzyme active site. Like in other AdSS, a conserved arginine residue (Arg155) is involved in dimer crosstalk and interacts with IMP in the active site of the symmetry related subunit of PfAdSS. We also report on the iochemical characterization of the arginine mutants (R155L, R155K and R155A) which suggests that unlike in E. coil AdSS, Arg155 in PfAdSS influences both ligand binding and catalysis. (C) 2010 Elsevier B.V. All rights reserved.

From anomalies in neat liquid to structure, dynamics and function in the biological world

Water brings its remarkable thermodynamic and dynamic anomalies in the pure liquid state to biological world where water molecules face a multitude of additional interactions that frustrate its hydrogen bond network. Yet the water molecules participate and control enormous number of biological processes in manners which are yet to be understood at a molecular level. We discuss thermodynamics, structure, dynamics and properties of water around proteins and DNA, along with those in reverse micelles. We discuss the roles of water in enzyme kinetics, in drug-DNA intercalation and in kinetic-proof reading ( the theory of lack of errors in biosynthesis). We also discuss how water may play an important role in the natural selection of biomolecules. (C) 2011 Elsevier B. V. All rights reserved.

On the Presumed Kinetic Consequences of Pre-equilibrium. Implications for the Michaelis-Menten Equation

The overall rate equation for a reaction sequence consisting of a pre-equilibrium and rate-determining steps should not be derived on the basis of the concentration of the intermediate product (X). This is apparently indicated by transition state theory (as the path followed to reach the highest energy transition state is irrelevant), but also proved by a straight-forward mathematical approach. The thesis is further supported by the equations of concurrent reactions as applied to the partitioning of X between the two competing routes (reversal of the pre-equilibrium and formation of product). The rate equation may only be derived rigorously on the basis of the law of mass action. It is proposed that the reactants acquire the overall activation energy prior to the pre-equilibrium, thus forming X in a high-energy state en route to the rate-determining transition state. (It is argued that conventional energy profile diagrams are misleading and need to be reinterpreted.) Also, these arguments invalidate the Michaelis-Menten equation of enzyme kinetics, and necessitate a fundamental revision of our present understanding of enzyme catalysis. (The observed ``saturation kinetics'' possibly arises from weak binding of a second molecule of substrate at the active site; analogous conclusions apply to reactions at surfaces).

Kinetic proofreading at single molecular level: aminoacylation of tRNA(Ile) and the role of water as an editor

Proofreading/editing in protein synthesis is essential for accurate translation of information from the genetic code. In this article we present a theoretical investigation of efficiency of a kinetic proofreading mechanism that employs hydrolysis of the wrong substrate as the discriminatory step in enzyme catalytic reactions. We consider aminoacylation of tRNA(Ile) which is a crucial step in protein synthesis and for which experimental results are now available. We present an augmented kinetic scheme and then employ methods of stochastic simulation algorithm to obtain time dependent concentrations of different substances involved in the reaction and their rates of formation. We obtain the rates of product formation and ATP hydrolysis for both correct and wrong substrates (isoleucine and valine in our case, respectively), in single molecular enzyme as well as ensemble enzyme kinetics. The present theoretical scheme correctly reproduces (i) the amplitude of the discrimination factor in the overall rates between isoleucine and valine which is obtained as (1.8x10(2)).(4.33x10(2)) = 7.8x10(4), (ii) the rates of ATP hydrolysis for both Ile and Val at different substrate concentrations in the aminoacylation of tRNA(Ile). The present study shows a non-michaelis type dependence of rate of reaction on tRNA(Ile) concentration in case of valine. The overall editing in steady state is found to be independent of amino acid concentration. Interestingly, the computed ATP hydrolysis rate for valine at high substrate concentration is same as the rate of formation of Ile-tRNA(Ile) whereas at intermediate substrate concentration the ATP hydrolysis rate is relatively low. We find that the presence of additional editing domain in class I editing enzyme makes the kinetic proofreading more efficient through enhanced hydrolysis of wrong product at the editing CP1 domain.

Electron flow through cytochrome P450

The cytochromes P450 (P450s) are a remarkable class of heme enzymes that catalyze the metabolism of xenobiotics and the biosynthesis of signaling molecules. Controlled electron flow into the thiolate-ligated heme active site allows P450s to activate molecular oxygen and hydroxylate aliphatic C–H bonds via the formation of high-valent metal-oxo intermediates (compounds I and II). Due to the reactive nature and short lifetimes of these intermediates, many of the fundamental steps in catalysis have not been observed directly. The Gray group and others have developed photochemical methods, known as “flash-quench,” for triggering electron transfer (ET) and generating redox intermediates in proteins in the absence of native ET partners. Photo-triggering affords a high degree of temporal precision for the gating of an ET event; the initial ET and subsequent reactions can be monitored on the nanosecond-to-second timescale using transient absorption (TA) spectroscopies. Chapter 1 catalogues critical aspects of P450 structure and mechanism, including the native pathway for formation of compound I, and outlines the development of photochemical processes that can be used to artificially trigger ET in proteins. Chapters 2 and 3 describe the development of these photochemical methods to establish electronic communication between a photosensitizer and the buried P450 heme. Chapter 2 describes the design and characterization of a Ru-P450-BM3 conjugate containing a ruthenium photosensitizer covalently tethered to the P450 surface, and nanosecond-to-second kinetics of the photo-triggered ET event are presented. By analyzing data at multiple wavelengths, we have identified the formation of multiple ET intermediates, including the catalytically relevant compound II; this intermediate is generated by oxidation of a bound water molecule in the ferric resting state enzyme. The work in Chapter 3 probes the role of a tryptophan residue situated between the photosensitizer and heme in the aforementioned Ru-P450 BM3 conjugate. Replacement of this tryptophan with histidine does not perturb the P450 structure, yet it completely eliminates the ET reactivity described in Chapter 2. The presence of an analogous tryptophan in Ru-P450 CYP119 conjugates also is necessary for observing oxidative ET, but the yield of heme oxidation is lower. Chapter 4 offers a basic description of the theoretical underpinnings required to analyze ET. Single-step ET theory is first presented, followed by extensions to multistep ET: electron “hopping.” The generation of “hopping maps” and use of a hopping map program to analyze the rate advantage of hopping over single-step ET is described, beginning with an established rhenium-tryptophan-azurin hopping system. This ET analysis is then applied to the Ru-tryptophan-P450 systems described in Chapter 2; this strongly supports the presence of hopping in Ru-P450 conjugates. Chapter 5 explores the implementation of flash-quench and other phototriggered methods to examine the native reductive ET and gas binding events that activate molecular oxygen. In particular, TA kinetics that demonstrate heme reduction on the microsecond timescale for four Ru-P450 conjugates are presented. In addition, we implement laser flash-photolysis of P450 ferrous–CO to study the rates of CO rebinding in the thermophilic P450 CYP119 at variable temperature. Chapter 6 describes the development and implementation of air-sensitive potentiometric redox titrations to determine the solution reduction potentials of a series of P450 BM3 mutants, which were designed for non-native cyclopropanation of styrene in vivo. An important conclusion from this work is that substitution of the axial cysteine for serine shifts the wild type reduction potential positive by 130 mV, facilitating reduction by biological redox cofactors in the presence of poorly-bound substrates. While this mutation abolishes oxygenation activity, these mutants are capable of catalyzing the cyclopropanation of styrene, even within the confines of an E. coli cell. Four appendices are also provided, including photochemical heme oxidation in ruthenium-modified nitric oxide synthase (Appendix A), general protocols (Appendix B), Chapter-specific notes (Appendix C) and Matlab scripts used for data analysis (Appendix D).

Design, fabrication and characterisation of components for microfluidic enzymatic biofuel cells

The concept of a biofuel cell takes inspiration from the natural capability of biological systems to catalyse the conversion of organic matter with a subsequent release of electrical energy. Enzymatic biofuel cells are intended to mimic the processes occurring in nature in a more controlled and efficient manner. Traditional fuel cells rely on the use of toxic catalysts and are often not easily miniaturizable making them unsuitable as implantable power sources. Biofuel cells however use highly selective protein catalysts and renewable fuels. As energy consumption becomes a global issue, they emerge as important tools for energy generation. The microfluidic platforms developed are intended to maximize the amount of electrical energy extracted from renewable fuels which are naturally abundant in the environment and in biological fluids. Combining microfabrication processes, chemical modification and biological surface patterning these devices are promising candidates for micro-power sources for future life science and electronic applications. This thesis considered four main aspects of a biofuel cell research. Firstly, concept of a miniature compartmentalized enzymatic biofuel cell utilizing simple fuels and operating in static conditions is verified and proves the feasibility of enzyme catalysis in energy conversion processes. Secondly, electrode and microfluidic channel study was performed through theoretical investigations of the flow and catalytic reactions which also improved understanding of the enzyme kinetics in the cell. Next, microfluidic devices were fabricated from cost-effective and disposable polymer materials, using the state-of-the-art micro-processing technologies. Integration of the individual components is difficult and multiple techniques to overcome these problems have been investigated. Electrochemical characterization of gold electrodes modified with Nanoporous Gold Structures is also performed. Finally, two strategies for enzyme patterning and encapsulation are discussed. Several protein catalysts have been effectively immobilized on the surface of commercial and microfabricated electrodes by electrochemically assisted deposition in sol-gel and poly-(o-phenylenediamine) polymer matrices and characterised with confirmed catalytic activity.