Evaluation of rate law approximations in bottom-up kinetic models of metabolism.

BackgroundThe mechanistic description of enzyme kinetics in a dynamic model of metabolism requires specifying the numerical values of a large number of kinetic parameters. The parameterization challenge is often addressed through the use of simplifying approximations to form reaction rate laws with reduced numbers of parameters. Whether such simplified models can reproduce dynamic characteristics of the full system is an important question.ResultsIn this work, we compared the local transient response properties of dynamic models constructed using rate laws with varying levels of approximation. These approximate rate laws were: 1) a Michaelis-Menten rate law with measured enzyme parameters, 2) a Michaelis-Menten rate law with approximated parameters, using the convenience kinetics convention, 3) a thermodynamic rate law resulting from a metabolite saturation assumption, and 4) a pure chemical reaction mass action rate law that removes the role of the enzyme from the reaction kinetics. We utilized in vivo data for the human red blood cell to compare the effect of rate law choices against the backdrop of physiological flux and concentration differences. We found that the Michaelis-Menten rate law with measured enzyme parameters yields an excellent approximation of the full system dynamics, while other assumptions cause greater discrepancies in system dynamic behavior. However, iteratively replacing mechanistic rate laws with approximations resulted in a model that retains a high correlation with the true model behavior. Investigating this consistency, we determined that the order of magnitude differences among fluxes and concentrations in the network were greatly influential on the network dynamics. We further identified reaction features such as thermodynamic reversibility, high substrate concentration, and lack of allosteric regulation, which make certain reactions more suitable for rate law approximations.ConclusionsOverall, our work generally supports the use of approximate rate laws when building large scale kinetic models, due to the key role that physiologically meaningful flux and concentration ranges play in determining network dynamics. However, we also showed that detailed mechanistic models show a clear benefit in prediction accuracy when data is available. The work here should help to provide guidance to future kinetic modeling efforts on the choice of rate law and parameterization approaches

Thermal inactivation and conformational lock studies on glucose oxidase

In this study, the dissociative thermal inactivation and conformational lock theories are applied for the homodimeric enzyme glucose oxidase (GOD) in order to analyze its structure. For this purpose, the rate of activity reduction of glucose oxidase is studied at various temperatures using b-D-glucose as the substrate by incubation of enzyme at various temperatures in the wide range between 40 and 70 �C using UV–Vis spectrophotometry. It was observed that in the two ranges of temperatures, the enzyme has two different forms. In relatively low temperatures, the enzyme is in its dimeric state and has normal activity. In high temperatures, the activity almost disappears and it aggregates. The above achievements are confirmed by dynamic light scattering. The experimental parameter ‘‘n’’ as the obvious number of conformational locks at the dimer interface of glucose oxidase is obtained by kinetic data, and the value is near to two. To confirm the above results, the X-ray crystallography structure of the enzyme, GOD (pdb, 1gal), was also studied. The secondary and tertiary structures of the enzyme to track the thermal inactivation were studied by circular dichroism and fluorescence spectroscopy, respectively. We proposed a mechanism model for thermal inactivation of GOD based on the absence of the monomeric form of the enzyme by circular dichroism and fluorescence spectroscopy

The interplay of intrinsic and extrinsic bounded noises in genetic networks

After being considered as a nuisance to be filtered out, it became recently clear that biochemical noise plays a complex role, often fully functional, for a genetic network. The influence of intrinsic and extrinsic noises on genetic networks has intensively been investigated in last ten years, though contributions on the co-presence of both are sparse. Extrinsic noise is usually modeled as an unbounded white or colored gaussian stochastic process, even though realistic stochastic perturbations are clearly bounded. In this paper we consider Gillespie-like stochastic models of nonlinear networks, i.e. the intrinsic noise, where the model jump rates are affected by colored bounded extrinsic noises synthesized by a suitable biochemical state-dependent Langevin system. These systems are described by a master equation, and a simulation algorithm to analyze them is derived. This new modeling paradigm should enlarge the class of systems amenable at modeling. We investigated the influence of both amplitude and autocorrelation time of a extrinsic Sine-Wiener noise on: $(i)$ the Michaelis-Menten approximation of noisy enzymatic reactions, which we show to be applicable also in co-presence of both intrinsic and extrinsic noise, $(ii)$ a model of enzymatic futile cycle and $(iii)$ a genetic toggle switch. In $(ii)$ and $(iii)$ we show that the presence of a bounded extrinsic noise induces qualitative modifications in the probability densities of the involved chemicals, where new modes emerge, thus suggesting the possibile functional role of bounded noises

Dynamic Disorder in Quasi-Equilibrium Enzymatic Systems

Conformations and catalytic rates of enzymes fluctuate over a wide range of timescales. Despite these fluctuations, there exist some limiting cases in which the enzymatic catalytic rate follows the macroscopic rate equation such as the Michaelis-Menten law. In this paper we investigate the applicability of macroscopic rate laws for fluctuating enzyme systems in which catalytic transitions are slower than ligand binding-dissociation reactions. In this quasi-equilibrium limit, for an arbitrary reaction scheme we show that the catalytic rate has the same dependence on ligand concentrations as obtained from mass-action kinetics even in the presence of slow conformational fluctuations. These results indicate that the timescale of conformational dynamics – no matter how slow – will not affect the enzymatic rate in quasi-equilibrium limit. Our numerical results for two enzyme-catalyzed reaction schemes involving multiple substrates and inhibitors further support our general theory

PPARγ Pro12Ala polymorphism and risk of acute coronary syndrome in a prospective study of Danes

<p>Abstract</p> <p>Background</p> <p>Acute coronary syndrome (ACS) is a major cause of morbidity and mortality in the western world. Peroxisome proliferator-activated receptor γ (PPARγ) plays a key role in the regulation of the energy balance, adipocyte differentiation and lipid biosynthesis. The aim was to investigate if the polymorphism <it>PPARγ2 </it>Pro¹²Ala, which encodes a less efficient transcription factor, was associated with risk of acute coronary disease and if there were interactions between this polymorphism and factors that modify PPARγ activity, such as alcohol intake, smoking, and use of non-steroidal anti-inflammatory medicine.</p> <p>Methods</p> <p>A case-cohort study including 1031 ACS cases and a sub-cohort of 1703 persons was nested within the population-based prospective study Diet, Cancer and Health of 57,053 individuals.</p> <p>Results</p> <p>Homozygous male variant allele carriers of <it>PPARγ2 </it>Pro¹²Ala were at higher risk of ACS (HR = 2.12, 95% CI: 1.00–4.48) than homozygous carriers of the Pro-allele. Among men, there was a statistically significant interaction between genotypes and alcohol intake such that homozygous variant allele carriers with a low alcohol intake were at higher risk of ACS (HR = 25.3, CI: 16.5–38.7) compared to homozygous common allele carriers (p for interaction < 0.0001). Overall, the association was only observed among homozygous variant allele carriers. Thus, all the observed associations were obtained in subgroups including small numbers of cases. It is therefore possible that the observed associations were due to chance.</p> <p>Conclusion</p> <p>In the present study, there were no consistent associations between PPARγ Pro¹²Ala and risk of ACS, and no consistent interaction with alcohol, BMI, NSAID or smoking in relation to ACS.</p

Pulmonary Vaccination as a Novel Treatment for Lung Fibrosis

Pulmonary fibrosis is an untreatable, uniformly fatal disease of unclear etiology that is the result of unremitting chronic inflammation. Recent studies have implicated bone marrow derived fibrocytes and M2 macrophages as playing key roles in propagating fibrosis. While the disease process is characterized by the accumulation of lymphocytes in the lung parenchyma and alveolar space, their role remains unclear. In this report we definitively demonstrate the ability of T cells to regulate lung inflammation leading to fibrosis. Specifically we demonstrate the ability of intranasal vaccinia vaccination to inhibit M2 macrophage generation and fibrocyte recruitment and hence the accumulation of collagen and death due to pulmonary failure. Mechanistically, we demonstrate the ability of lung Th1 cells to prevent fibrosis as vaccinia failed to prevent disease in Rag−/− mice or in mice in which the T cells lacked IFN-γ. Furthermore, vaccination 3 months prior to the initiation of fibrosis was able to mitigate the disease. Our findings clearly demonstrate the role of T cells in regulating pulmonary fibrosis as well as suggest that vaccinia-induced immunotherapy in the lung may prove to be a novel treatment approach to this otherwise fatal disease

Monotone and near-monotone biochemical networks

Monotone subsystems have appealing properties as components of larger networks, since they exhibit robust dynamical stability and predictability of responses to perturbations. This suggests that natural biological systems may have evolved to be, if not monotone, at least close to monotone in the sense of being decomposable into a “small” number of monotone components, In addition, recent research has shown that much insight can be attained from decomposing networks into monotone subsystems and the analysis of the resulting interconnections using tools from control theory. This paper provides an expository introduction to monotone systems and their interconnections, describing the basic concepts and some of the main mathematical results in a largely informal fashion

Salivary gland branching morphogenesis: a quantitative systems analysis of the Eda/Edar/NFκB paradigm

<p>Abstract</p> <p>Background</p> <p>Ectodysplasin-A appears to be a critical component of branching morphogenesis. Mutations in mouse <it>Eda </it>or human <it>EDA </it>are associated with absent or hypoplastic sweat glands, sebaceous glands, lacrimal glands, salivary glands (SMGs), mammary glands and/or nipples, and mucous glands of the bronchial, esophageal and colonic mucosa. In this study, we utilized <it>Eda</it>^{<it>Ta </it>}(Tabby) mutant mice to investigate how a marked reduction in functional Eda propagates with time through a defined genetic subcircuit and to test the proposition that canonical NFκB signaling is sufficient to account for the differential expression of developmentally regulated genes in the context of <it>Eda </it>polymorphism.</p> <p>Results</p> <p>The quantitative systems analyses do not support the stated hypothesis. For most NFκB-regulated genes, the observed time course of gene expression is nearly unchanged in Tabby (<it>Eda</it>^<it>Ta</it>) as compared to wildtype mice, as is NFκB itself. Importantly, a subset of genes is dramatically differentially expressed in Tabby (<it>Edar</it>, <it>Fgf8</it>, <it>Shh</it>, <it>Egf</it>, <it>Tgfa</it>, <it>Egfr</it>), strongly suggesting the existence of an alternative Eda-mediated transcriptional pathway pivotal for SMG ontogeny. Experimental and <it>in silico </it>investigations have identified C/EBPα as a promising candidate.</p> <p>Conclusion</p> <p>In Tabby SMGs, upregulation of the Egf/Tgfα/Egfr pathway appears to mitigate the potentially severe abnormal phenotype predicted by the downregulation of Fgf8 and Shh. Others have suggested that the buffering of the phenotypic outcome that is coincident with variant Eda signaling could be a common mechanism that permits viable and diverse phenotypes, normal and abnormal. Our results support this proposition. Further, if branching epithelia use variations of a canonical developmental program, our results are likely applicable to understanding the phenotypes of other branching organs affected by <it>Eda </it>(<it>EDA</it>) mutation.</p