Evaluation of steady-state kinetic parameters for enzymes solubilized in water-in-oil microemulsion systems.
Ontology highlight
ABSTRACT: 1. Equations are derived for the steady-state kinetics of substrate conversion by enzymes confined within the water-droplets of water-in-oil microemulsion systems. 2. Water-soluble substrates initially confined within droplets that do not contain enzyme are assumed to be converted into product only after they enter enzyme-containing droplets via the inter-droplet exchange process. 3. Hyperbolic (Michaelis-Menten) kinetics are predicted when the substrate concentration is varied in microemulsions of fixed composition. Both kcat. and Km are predicted to be dependent on the size and concentration of the water-droplets in the microemulsion. 4. The predicted behaviour is shown to be supported by published experimental data. A physical interpretation of the form of the rate equation is presented. 5. The rate equation for an oil-soluble substrate was derived assuming a pseudo-two-phase (oil & water) model for the microemulsion. Both kcat. and Km are shown to be independent of phi aq. Km is larger than the aqueous solution value by a factor approximately equal to the oil/water partition coefficient of the substrate. The validity of the rate equation is confirmed by published data.
Project description:Nicotinamidases are metabolic enzymes that hydrolyze nicotinamide to nicotinic acid. These enzymes are widely distributed across biology, with examples found encoded in the genomes of Mycobacteria, Archaea, Eubacteria, Protozoa, yeast, and invertebrates, but there are none found in mammals. Although recent structural work has improved our understanding of these enzymes, their catalytic mechanism is still not well understood. Recent data show that nicotinamidases are required for the growth and virulence of several pathogenic microbes. The enzymes of Saccharomyces cerevisiae, Drosophila melanogaster, and Caenorhabditis elegans regulate life span in their respective organisms, consistent with proposed roles in the regulation of NAD(+) metabolism and organismal aging. In this work, the steady state kinetic parameters of nicotinamidase enzymes from C. elegans, Sa. cerevisiae, Streptococcus pneumoniae (a pathogen responsible for human pneumonia), Borrelia burgdorferi (the pathogen that causes Lyme disease), and Plasmodium falciparum (responsible for most human malaria) are reported. Nicotinamidases are generally efficient catalysts with steady state k(cat) values typically exceeding 1 s(-1). The K(m) values for nicotinamide are low and in the range of 2 -110 ?M. Nicotinaldehyde was determined to be a potent competitive inhibitor of these enzymes, binding in the low micromolar to low nanomolar range for all nicotinamidases tested. A variety of nicotinaldehyde derivatives were synthesized and evaluated as inhibitors in kinetic assays. Inhibitions are consistent with reaction of the universally conserved catalytic Cys on each enzyme with the aldehyde carbonyl carbon to form a thiohemiacetal complex that is stabilized by a conserved oxyanion hole. The S. pneumoniae nicotinamidase can catalyze exchange of (18)O into the carboxy oxygens of nicotinic acid with H(2)(18)O. The collected data, along with kinetic analysis of several mutants, allowed us to propose a catalytic mechanism that explains nicotinamidase and nicotinic acid (18)O exchange chemistry for the S. pneumoniae enzyme involving key catalytic residues, a catalytic transition metal ion, and the intermediacy of a thioester intermediate.
Project description:Mycobacterium tuberculosis and many other members of the Actinomycetes family produce mycothiol, i.e., 1-d-myo-inosityl-2-(N-acetyl-l-cysteinyl)amido-2-deoxy-alpha-d-glucopyranoside (MSH or AcCys-GlcN-Ins), to act against oxidative and antibiotic stress. The biosynthesis of MSH is essential for cell growth and has been proposed to proceed via a biosynthetic pathway involving four key enzymes, MshA-MshD. The MSH biosynthetic enzymes present potential targets for inhibitor design. With this as a long-term goal, we have carried out a kinetic and mechanistic characterization, using steady-state and pre-steady-state approaches, of the recombinant Mycobacterium smegmatis MshC. MshC catalyzes the ATP-dependent condensation of GlcN-Ins and cysteine to form Cys-GlcN-Ins. Initial velocity and inhibition studies show that the steady-state kinetic mechanism of MshC is a Bi Uni Uni Bi Ping Pong mechanism, with ATP binding followed by cysteine binding, release of PPi, binding of GlcN-Ins, followed by the release of Cys-GlcN-Ins and AMP. The steady-state kinetic parameters were determined to be kcat equal to 3.15 s-1, and Km values of 1.8, 0.1, and 0.16 mM for ATP, cysteine, and GlcN-Ins, respectively. A stable bisubstrate analogue, 5'-O-[N-(l-cysteinyl)sulfamonyl]adenosine, exhibits competitive inhibition versus ATP and noncompetitive inhibition versus cysteine, with an inhibition constant of approximately 306 nM versus ATP. Single-turnover reactions of the first and second half reactions were determined using rapid-quench techniques, giving rates of approximately 9.4 and approximately 5.2 s-1, respectively, consistent with the cysteinyl adenylate being a kinetically competent intermediate in the reaction by MshC.
Project description:The quantitative analysis of biochemical reactions and metabolites is at frontier of biological sciences. The recent availability of high-throughput technology data sets in biology has paved the way for new modelling approaches at various levels of complexity including the metabolome of a cell or an organism. Understanding the metabolism of a single cell and multi-cell organism will provide the knowledge for the rational design of growth conditions to produce commercially valuable reagents in biotechnology. Here, we demonstrate how equations representing steady state mass conservation, energy conservation, the second law of thermodynamics, and reversible enzyme kinetics can be formulated as a single system of linear equalities and inequalities, in addition to linear equalities on exponential variables. Even though the feasible set is non-convex, the reformulation is exact and amenable to large-scale numerical analysis, a prerequisite for computationally feasible genome scale modelling. Integrating flux, concentration and kinetic variables in a unified constraint-based formulation is aimed at increasing the quantitative predictive capacity of flux balance analysis. Incorporation of experimental and theoretical bounds on thermodynamic and kinetic variables ensures that the predicted steady state fluxes are both thermodynamically and biochemically feasible. The resulting in silico predictions are tested against fluxomic data for central metabolism in Escherichia coli and compare favourably with in silico prediction by flux balance analysis.
Project description:In pharmacology and systems biology, it is a fundamental problem to determine how biological systems change their dose-response behavior upon perturbations. In particular, it is unclear how topologies, reactions, and parameters (differentially) affect the dose response. Because parameters are often unknown, systematic approaches should directly relate network structure and function. Here, we outline a procedure to compare general non-monotone dose-response curves and subsequently develop a comprehensive theory for differential dose responses of biochemical networks captured by non-equilibrium steady-state linear framework models. Although these models are amenable to analytical derivations of non-equilibrium steady states in principle, their size frequently increases (super) exponentially with model size. We extract general principles of differential responses based on a model's graph structure and thereby alleviate the combinatorial explosion. This allows us, for example, to determine reactions that affect differential responses, to identify classes of networks with equivalent differential, and to reject hypothetical models reliably without needing to know parameter values. We exemplify such applications for models of insulin signaling.
Project description:It is shown by analysis of a numerical example that kinetic barrier diagrams [Südi (1991) Biochem. J. 276; 265-268] are also useful for displaying the phenomenological resistance of an enzymic turnover to net chemical fluxes observed under steady-state conditions. Most importantly, a new additivity rule is revealed by net flux profiles, which refers to limiting ('ideal') conditions. The rule states that the net fluxes that one obtains under initial steady-state conditions are strictly identical with the overall flux in the corresponding equilibrium system. This finding amounts to defining the relation between unidirectional fluxes and net fluxes, and explains why the total 'one-way flux resistance' to the overall reaction at chemical equilibrium on the one hand, and the 'net flux resistance' to both initial steady-state reactions on the other hand, have to be equal. The conclusions are claimed to be generally valid for consecutive chemical reactions.
Project description:Understanding how cells change their identity and behaviour in living systems is an important question in many fields of biology. The problem of inferring cell trajectories from single-cell measurements has been a major topic in the single-cell analysis community, with different methods developed for equilibrium and non-equilibrium systems (e.g. haematopoeisis vs. embryonic development). We show that optimal transport analysis, a technique originally designed for analysing time-courses, may also be applied to infer cellular trajectories from a single snapshot of a population in equilibrium. Therefore, optimal transport provides a unified approach to inferring trajectories that is applicable to both stationary and non-stationary systems. Our method, StationaryOT, is mathematically motivated in a natural way from the hypothesis of a Waddington's epigenetic landscape. We implement StationaryOT as a software package and demonstrate its efficacy in applications to simulated data as well as single-cell data from Arabidopsis thaliana root development.
Project description:Gap junction (GJ) channels, formed of connexin (Cx) proteins, provide a direct pathway for metabolic and electrical cell-to-cell communication. These specialized channels are not just passive conduits for the passage of ions and metabolites but have been shown to gate robustly in response to transjunctional voltage, Vj, the voltage difference between two coupled cells. Voltage gating of GJs could play a physiological role, particularly in excitable cells, which can generate large transients in membrane potential during the propagation of action potentials. We present a mathematical/computational model of GJ channel voltage gating to assess properties of GJ channels that takes into account contingent gating of two series hemichannels and the distribution of Vj across each hemichannel. From electrophysiological recordings in cell cultures expressing Cx43 or Cx45, the principal isoforms expressed in cardiac tissue, various data sets were fitted simultaneously using global optimization. The results showed that the model is capable of describing both steady-state and kinetic properties of homotypic and heterotypic GJ channels composed of these Cxs. Moreover, mathematical analyses showed that the model can be simplified to a reversible two-state system and solved analytically using a rapid equilibrium assumption. Given that excitable cells are arranged in interconnected networks, the equilibrium assumption allows for a substantial reduction in computation time, which is useful in simulations of large clusters of coupled cells. Overall, this model can serve as a tool for the studying of GJ channel gating and its effects on the spread of excitation in networks of electrically coupled cells.
Project description:Cruzain, an important drug target for Chagas disease, is a member of clan CA of the cysteine proteases. Understanding the catalytic mechanism of cruzain is vital to the design of new inhibitors. To this end, we have determined pH-rate profiles for substrates and affinity agents and solvent kinetic isotope effects in pre-steady-state and steady-state modes using three substrates: Cbz-Phe-Arg-AMC, Cbz-Arg-Arg-AMC, and Cbz-Arg-Ala-AMC. The pH-rate profile of kcat/ Km for Cbz-Arg-Arg-AMC indicated p K1 = 6.6 (unprotonated) and p K2 ∼ 9.6 (protonated) groups were required for catalysis. The temperature dependence of the p K = 6.2-6.6 group exhibited a Δ Hion value of 8.4 kcal/mol, typical of histidine. The pH-rate profile of inactivation by iodoacetamide confirmed that the catalytic cysteine possesses a p Ka of 9.8. Normal solvent kinetic isotope effects were observed for both D2O kcat (1.6-2.1) and D2O kcat/ Km (1.1-1.4) for all three substrates. Pre-steady-state kinetics revealed exponential bursts of AMC production for Cbz-Phe-Arg-AMC and Cbz-Arg-Arg-AMC, but not for Cbz-Arg-Ala-AMC. The overall solvent isotope effect on kcat can be attributed to the solvent isotope effect on the deacylation step. Our results suggest that cruzain is unique among papain-like cysteine proteases in that the catalytic cysteine and histidine have neutral charges in the free enzyme. The generation of the active thiolate of the catalytic cysteine is likely preceded (and possibly triggered) by a ligand-induced conformational change, which could bring the catalytic dyad into the proximity to effect proton transfer.
Project description:Many DNA modification and repair enzymes require access to DNA bases and therefore flip nucleotides. Restriction endonucleases (REases) hydrolyze the phosphodiester backbone within or in the vicinity of the target recognition site and do not require base extrusion for the sequence readout and catalysis. Therefore, the observation of extrahelical nucleotides in a co-crystal of REase Ecl18kI with the cognate sequence, CCNGG, was unexpected. It turned out that Ecl18kI reads directly only the CCGG sequence and skips the unspecified N nucleotides, flipping them out from the helix. Sequence and structure conservation predict nucleotide flipping also for the complexes of PspGI and EcoRII with their target DNAs (/CCWGG), but data in solution are limited and indirect. Here, we demonstrate that Ecl18kI, the C-terminal domain of EcoRII (EcoRII-C) and PspGI enhance the fluorescence of 2-aminopurines (2-AP) placed at the centers of their recognition sequences. The fluorescence increase is largest for PspGI, intermediate for EcoRII-C and smallest for Ecl18kI, probably reflecting the differences in the hydrophobicity of the binding pockets within the protein. Omitting divalent metal cations and mutation of the binding pocket tryptophan to alanine strongly increase the 2-AP signal in the Ecl18kI-DNA complex. Together, our data provide the first direct evidence that Ecl18kI, EcoRII-C and PspGI flip nucleotides in solution.
Project description:A cytosolic cell-free system from rat liver containing the last three enzymes of the urea cycle, a number of cofactors and the substrates aspartate and citrulline was shown to synthesize urea at near-physiological rates ranging between 0.40 and 1.25 mumol/min per g of liver. This system was used to determine the kinetic parameters for arginase. With saturating amounts of Mn2+ (30 microM), arginine remained at a steady-state concentration of 5-35 microM depending on the aspartate and citrulline supply. Vmax. at micromolar arginine concentrations was between 1.10 and 1.25 mumol/min per g of liver, the K0.5 (arginine) between 6.0 and 6.5 microM and positive co-operativity was observed (Hill coefficient 2). Omission of Mn2+ caused a significant accumulation of arginine during the incubation, suggesting a regulatory effect of arginase. Under these conditions, Vmax. was 1.10-1.65 mumol/min per g of liver and the Km (arginine) increased up to 14.4-21.1 microM. The apparent Ka for Mn2+ in the presence of physiological concentrations of ATP, Mg2+ and arginine was calculated to be maximally 8 microM. Initial-velocity experiments with millimolar arginine concentrations as the direct substrate gave the following results, which are in good agreement with literature data. In the absence of Mn2+, Vmax. was 71.3 mumol/min per g of liver and the Km (arginine) 1.58 mM. With 30 microM-Mn2+, Vmax. was 69.4 mumol/min per g of liver and the Km (arginine) decreased to 0.94 mM. On the basis of our results, we propose the presence of high-affinity and low-affinity sites for arginine on rat liver arginase and postulate that alterations in arginase activity arising from changes in the concentration of arginine and of the cofactor Mn2+ may contribute to the regulation of ureagenesis in vivo.