Modeling physiological processes in plankton on enzyme kinetic principles

[EN] Many ecologically important chemical transformations in the ocean are controlled by biochemical enzyme reactions in plankton. Nitrogenase regulates the transformation of N2 to ammonium in some cyanobacteria and serves as the entryway for N2 into the ocean biosphere. Nitrate reductase controls the reduction of NO3 to NO2 and hence new production in phytoplankton. The respiratory electron transfer system in all organisms links the carbon oxidation reactions of intermediary metabolism with the reduction of oxygen in respiration. Rubisco controls the fixation of CO2 into organic matter in phytoplankton and thus is the major entry point of carbon into the oceanic biosphere. In addition to these, there are the enzymes that control CO2 production, NH4 excretion and the fluxes of phosphate. Some of these enzymes have been recognized and researched by marine scientists in the last thirty years. However, until recently the kinetic principles of enzyme control have not been exploited to formulate accurate mathematical equations of the controlling physiological expressions. Were such expressions available they would increase our power to predict the rates of chemical transformations in the extracellular environment of microbial populations whether this extracellular environment is culture media or the ocean. Here we formulate from the principles of bisubstrate enzyme kinetics, mathematical expressions for the processes of NO3 reduction, O2 consumption, N2 fixation, total nitrogen uptake.

Efficient and cost-effective experimental determination of kinetic constants and data: the success of a Bayesian systematic approach to drug transport, receptor binding, continuous culture and cell transport kinetics

Details about the parameters of kinetic systems are crucial for progress in both medical and industrial research, including drug development, clinical diagnosis and biotechnology applications. Such details must be collected by a series of kinetic experiments and investigations. The correct design of the experiment is essential to collecting data suitable for analysis, modelling and deriving the correct information. We have developed a systematic and iterative Bayesian method and sets of rules for the design of enzyme kinetic experiments. Our method selects the optimum design to collect data suitable for accurate modelling and analysis and minimises the error in the parameters estimated. The rules select features of the design such as the substrate range and the number of measurements. We show here that this method can be directly applied to the study of other important kinetic systems, including drug transport, receptor binding, microbial culture and cell transport kinetics. It is possible to reduce the errors in the estimated parameters and, most importantly, increase the efficiency and cost-effectiveness by reducing the necessary amount of experiments and data points measured. (C) 2003 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.

Bayesian L-optimal exact design of experiments for biological kinetic models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Molecular and Kinetic Characterization of Two Extracellular Xylanases Isolated from Leucoagaricus gongylophorus

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

MONITORING ENZYME REACTIONS BY FAST ATOM BOMBARDMENT MASS SPECTROMETRY

The potential for the direct analysis of enzyme reactions by fast atom bombardment (FAB) mass spectrometry has been investigated. Conditions are presented for the maintenance of enzymatic activity under FAB conditions along with FAB mass spectrometric data showing that these conditions can be applied to solutions of enzyme and substrate to follow enzymatic reactions inside the mass spectrometer in real-time. In addition, enzyme kinetic behavior under FAB mass spectrometric conditions is characterized using trypsin and its assay substrate, TAME, as an enzyme-substrate reaction model. These results show that two monitoring methods can be utilized to follow reactions by FAB mass spectrometry. The advantages of each method are discussed and illustrated by obtaining kinetic parameters from the direct analysis of enzyme reactions with assay or peptide substrates. ^

A potential role for isothermal calorimetry in studies of the effects of thermodynamic non-ideality in enzyme-catalyzed reactions

Attention is drawn to the feasibility of using isothermal calorimetry for the characterization of enzyme reactions under conditions bearing greater relevance to the crowded biological environment, where kinetic parameters are likely to differ significantly from those obtained by classical enzyme kinetic studies in dilute solution. An outline of the application of isothermal calorimetry to the determination of enzyme kinetic parameters is followed by considerations of the nature and consequences of crowding effects in enzyme catalysis. Some of those effects of thermodynamic non-ideality are then illustrated by means of experimental results from calorimetric studies of the effect of molecular crowding on the kinetics of catalysis by rabbit muscle pyruvate kinase. This review concludes with a discussion of the potential of isothermal calorimetry for the experimental determination of kinetic parameters for enzymes either in biological environments or at least in media that should provide reasonable approximations of the crowded conditions encountered in vivo. Copyright (C) 2004 John Wiley Sons, Ltd.

Kinetic modelling of in vitro data of PI3K, mTOR1, PTEN enzymes and on-target inhibitors Rapamycin, BEZ235, and LY294002

The phosphatidylinositide 3-kinases (PI3K) and mammalian target of rapamycin-1 (mTOR1) are two key targets for anti-cancer therapy. Predicting the response of the PI3K/AKT/mTOR1 signalling pathway to targeted therapy is made difficult because of network complexities. Systems biology models can help explore those complexities but the value of such models is dependent on accurate parameterisation. Motivated by a need to increase accuracy in kinetic parameter estimation, and therefore the predictive power of the model, we present a framework to integrate kinetic data from enzyme assays into a unified enzyme kinetic model. We present exemplar kinetic models of PI3K and mTOR1, calibrated on in vitro enzyme data and founded on Michaelis-Menten (MM) approximation. We describe the effects of an allosteric mTOR1 inhibitor (Rapamycin) and ATP-competitive inhibitors (BEZ2235 and LY294002) that show dual inhibition of mTOR1 and PI3K. We also model the kinetics of phosphatase and tensin homolog (PTEN), which modulates sensitivity of the PI3K/AKT/mTOR1 pathway to these drugs. Model validation with independent data sets allows investigation of enzyme function and drug dose dependencies in a wide range of experimental conditions. Modelling of the mTOR1 kinetics showed that Rapamycin has an IC50 independent of ATP concentration and that it is a selective inhibitor of mTOR1 substrates S6K1 and 4EBP1: it retains 40% of mTOR1 activity relative to 4EBP1 phosphorylation and inhibits completely S6K1 activity. For the dual ATP-competitive inhibitors of mTOR1 and PI3K, LY294002 and BEZ235, we derived the dependence of the IC50 on ATP concentration that allows prediction of the IC50 at different ATP concentrations in enzyme and cellular assays. Comparison of the drug effectiveness in enzyme and cellular assays showed that some features of these drugs arise from signalling modulation beyond the on-target action and MM approximation and require a systems-level consideration of the whole PI3K/PTEN/AKT/mTOR1 network in order to understand mechanisms of drug sensitivity and resistance in different cancer cell lines. We suggest that using these models in systems biology investigation of the PI3K/AKT/mTOR1 signalling in cancer cells can bridge the gap between direct drug target action and the therapeutic response to these drugs and their combinations.

Some novel techniques of parameter estimation for the dynamical models in biological systems

Inverse problems based on using experimental data to estimate unknown parameters of a system often arise in biological and chaotic systems. In this paper, we consider parameter estimation in systems biology involving linear and non-linear complex dynamical models, including the Michaelis–Menten enzyme kinetic system, a dynamical model of competence induction in Bacillus subtilis bacteria and a model of feedback bypass in B. subtilis bacteria. We propose some novel techniques for inverse problems. Firstly, we establish an approximation of a non-linear differential algebraic equation that corresponds to the given biological systems. Secondly, we use the Picard contraction mapping, collage methods and numerical integration techniques to convert the parameter estimation into a minimization problem of the parameters. We propose two optimization techniques: a grid approximation method and a modified hybrid Nelder–Mead simplex search and particle swarm optimization (MH-NMSS-PSO) for non-linear parameter estimation. The two techniques are used for parameter estimation in a model of competence induction in B. subtilis bacteria with noisy data. The MH-NMSS-PSO scheme is applied to a dynamical model of competence induction in B. subtilis bacteria based on experimental data and the model for feedback bypass. Numerical results demonstrate the effectiveness of our approach.

Understanding Human Drug Conjugating Enzymes; Regio- and Stereoselectivity in Sulfotransferase 1A3 and the UDP-Glucuronosyltransferases

Sulfotransferases (SULTs) and UDP-glucuronosyltransferases (UGTs) are important detoxification enzymes and they contribute to bioavailability and elimination of many drugs. SULT1A3 is an extrahepatic enzyme responsible for the sulfonation of dopamine, which is often used as its probe substrate. A new method for analyzing dopamine-3-O-sulfate and dopamine-4-O-sulfate by high-performance liquid chromatography was developed and the enzyme kinetic parameters for their formation were determined using purified recombinant human SULT1A3. The results show that SULT1A3 strongly favors the 3-hydroxy group of dopamine, which indicates that it may be the major enzyme responsible for the difference between the circulating levels of dopamine sulfates in human blood. All 19 known human UGTs were expressed as recombinant enzymes in baculovirus infected insect cells and their activities toward dopamine and estradiol were studied. UGT1A10 was identified as the only UGT capable of dopamine glucuronidation at a substantial level. The results were supported by studies with human intestinal and liver microsomes. The affinity was low indicating that UGT1A10 is not an important enzyme in dopamine metabolism in vivo. Despite the low affinity, dopamine is a potential new probe substrate for UGT1A10 due to its selectivity. Dopamine was used to study the importance of phenylalanines 90 and 93 in UGT1A10. The results revealed distinct effects that are dependent on differences in the size of the side chain and on the differences in their position within the protein. Examination of twelve mutants revealed lower activity in all of them. However, the enzyme kinetic studies of four mutants showed that their affinities were similar to that of UGT1A10 suggesting that F90 and F93 are not directly involved in dopamine binding in the active site. The glucuronidation of β-estradiol and epiestradiol (α-estradiol) was studied to elucidate how the orientation of the 17-OH group affects conjugation at the 3-OH or the 17-OH of either diastereomer. The results show that there are clear differences in the regio- and stereoselectivities of UGTs. The most active isoforms were UGT1A10 and UGT2B7 demonstrating opposite regioselectivity. The stereoselectivities of UGT2Bs were more complex than those of UGT1As. The amino acid sequences of the human UGTs 1A9 and 1A10 are 93% identical, yet there are large differences in their activity and substrate selectivity. Several mutants were constructed to identify the residues responsible for the activity differences. The results revealed that the residues between Leu86 and Tyr176 of UGT1A9 determine the differences between UGT1A9 and UGT1A10. Phe117 of UGT1A9 participated in 1-naphthol binding and the residues at positions 152 and 169 contributed to the higher glucuronidation rates of UGT1A10. In summary, the results emphasize that the substrate selectivities, including regio- and stereoselectivities, of UGTs are complex and they are controlled by many amino acids rather than one critical residue.

A Kunitz-type serine protease inhibitor from Butea monosperma seed and its influence on developmental physiology of Helicoverpa armigera

A protease inhibitor from the seeds of Butea monosperma (BmPI) was purified, characterized and studied for its influence on developmental physiology of Helicover-pa armigera. BmPI on two-dimensional separations indicated the presence of a 14 kDa protein with an isoelectric point in the acidic region (pl 5.6). Multiple Sequence Analysis data suggested that the BmPI contains a sequence motif which is conserved in various trypsin and chymotrypsin inhibitors of Kunitz-type. The inhibitor exhibited trypsin inhibitory activity in a broad range of pH (4-10) and temperature (10-80 degrees C). The enzyme kinetic studies revealed BmPI as a competitive inhibitor with a K-i value of 1.2 x 10(-9) M. In vitro studies with BmPI indicated measurable inhibitory activity on total gut proteolytic enzymes of H. armigera (IC(50)2.0 mu g/ml) and bovine trypsin. BmPI supplemented artificial diet caused dose dependent mortality and reduction in growth and weight. The fertility and fecundity of H. armigera, declined whereas the larval-pupal duration of the insect life cycle extended. These detrimental effects on H. armigera suggest the usefulness of BmPl in insect pest management of food crops. (C) 2014 Elsevier Ltd. All rights reserved.

Arylsulfatase activity in marine polychaetes

Marine polychaetes, collected from the Vellar Estuary exhibited arylsulfatase activity. Lumbriconeries sp. Polydora sp. Monojis sp. and Heteromastus sp. were selected for this study. Of these, Heteromastus sp. showed maximum enzymatic activity and it has been chosen for the enzyme kinetic studies such as pH, optimal temperature, period of incubation and the effect of DDT. Enzyme activity showed single peak at pH 6.2 possibly indicating the presence of one type of arylsulfatase. Maximum activity was attained after 12h of incubation at 29°C. DDT has an inhibiting effect on the arylsulfatase activity even at the concentration of 10 p.p.m. and the activity was completely lost at 100 p.p.m.

Characterization of the highly conserved VFMGD motif in a bacterial polyisoprenyl-phosphate N-acetylaminosugar-1-phosphate transferase

Polyisoprenyl-phosphate N-acetylaminosugar-1-phosphate transferases (PNPTs) constitute a family of eukaryotic and prokaryotic membrane proteins that catalyze the transfer of a sugar-1-phosphate to a phosphoisoprenyl lipid carrier. All PNPT members share a highly conserved 213-Valine-Phenylalanine-Methionine-Glycine-Aspartic acid-217 (VFMGD) motif. Previous studies using the MraY protein suggested that the aspartic acid residue in this motif, D267, is a nucleophile for a proposed double-displacement mechanism involving the cleavage of the phosphoanhydride bond of the nucleoside. Here, we demonstrate that the corresponding residue in the E. coli WecA, D217, is not directly involved in catalysis, as its replacement by asparagine results in a more active enzyme. Kinetic data indicate that the D217N replacement leads to more than twofold increase in V(max) without significant change in the K(m) for the nucleoside sugar substrate. Furthermore, no differences in the binding of the reaction intermediate analog tunicamycin were found in D217N as well as in other replacement mutants at the same position. We also found that alanine substitutions in various residues of the VFMGD motif affect to various degrees the enzymatic activity of WecA in vivo and in vitro. Together, our data suggest that the highly conserved VFMGD motif defines a common region in PNPT proteins that contributes to the active site and is likely involved in the release of the reaction product.

The Saccharomyces cerevisiae quinone oxidoreductase Lot6p: stability, inhibition and cooperativity

Lot6p (EC 1.5.1.39; Ylr011wp) is the sole quinone oxidoreductase in the budding yeast, Saccharomyces cerevisiae. Using hexahistidine tagged, recombinant Lot6p, we determined the steady-state enzyme kinetic parameters with both NADH and NADPH as electron donors; no cooperativity was observed with these substrates. The NQO1 inhibitor curcumin, the NQO2 inhibitor resveratrol, the bacterial nitroreductase inhibitor nicotinamide and the phosphate mimic vanadate all stabilise the enzyme towards thermal denaturation as judged by differential scanning fluorimetry. All except vanadate have no observable effect on the chemical cross-linking of the two subunits of the Lot6p dimer. These compounds all inhibit Lot6p's oxidoreductase activity, and all except nicotinamide exhibit negative cooperativity. Molecular modelling suggests that curcumin, resveratrol and nicotinamide all bind over the isoalloxazine ring of the FMN cofactor in Lot6p. Resveratrol was predicted to contact an α-helix that links the two active sites. Mutation of Gly-142 (which forms part of this helix) to serine does not greatly affect the thermal stability of the enzyme. However, this variant shows less cooperativity towards resveratrol than the wild type. This suggests a plausible hypothesis for the transmission of information between the subunits and, thus, the molecular mechanism of negative cooperativity in Lot6p.

Inhibition réversible et photomarquage de la transglutaminase tissulaire

La transglutaminase tissulaire est une enzyme dépendante du calcium qui catalyse la formation de liens isopeptidiques, entre les chaînes latérales de résidus glutamine et lysine, permettant, par le fait même, la réticulation des protéines dans les systèmes biologiques. Elle joue un rôle, entre autres, dans l’endocytose, la régulation du développement des cellules, et même dans l’apoptose. Néanmoins, une dérégulation de l’activité biologique de cette enzyme peut entrainer différentes pathologies, comme la formation de cataractes, de plaques amyloïdes dans la maladie d’Alzheimer, ou encore peut mener au développement de la maladie céliaque. C’est pourquoi une meilleure connaissance du mécanisme d’action de cette enzyme et la possibilité de réguler son action à l’aide de substrats ou d’inhibiteurs sont nécessaires. Dans cette optique, une méthode d’expression et de purification de la transglutaminase humaine a été développée, permettant de travailler directement avec la cible pharmacologique désirée. De plus, une étude du mode d’inhibition et de liaison d’une classe d’inhibiteurs réversibles précédemment découverte dans le groupe, soit la famille des trans-cinnamoyles, a permis d’identifier que la puissance de ces molécules est influencée par la présence du calcium et qu’une inhibition dépendante du temps est observée, en lien avec un potentiel équilibre conformationnel lent de la transglutaminase. D’un autre côté, la susceptibilité à une attaque nucléophile par des thiols de cette classe de molécule rend leur potentiel pharmacologique grandement diminué, et c’est pourquoi une nouvelle famille de molécules a été identifiée, basée sur un squelette ynone, avec une valeur d’IC50 très prometteuse de 2,6 μM, en faisant un des meilleurs inhibiteurs réversibles de la transglutaminase développés à ce jour. Finalement, une stratégie de photomarquage jumelée à une analyse de spectrométrie de masse en tandem a été développée pour la découverte du site de liaison du substrat dérivé de la lysine, dans le but de mieux comprendre le mécanisme complexe de cette enzyme.