Effect of substrate competition in kinetic models of metabolic networks

Cardiac myosin binding protein-C phosphorylation as a function of multiple protein kinase and phosphatase activities

Phosphorylation of cardiac myosin binding protein-C (cMyBP-C) is a determinant of cardiac myofilament function. Although cMyBP-C phosphorylation by various protein kinases has been extensively studied, the influence of protein phosphatases on cMyBP-C's multiple phosphorylation sites has remained largely obscure. Here we provide a detailed biochemical characterization of cMyBP-C dephosphorylation by protein phosphatases 1 and 2 A (PP1 and PP2A), and develop an integrated kinetic model for cMyBP-C phosphorylation using data for both PP1, PP2A and various protein kinases known to phosphorylate cMyBP-C. We find strong site-specificity and a hierarchical mechanism for both phosphatases, proceeding in the opposite direction of sequential phosphorylation by potein kinase A. The model is consistent with published data from human patients and predicts complex non-linear cMyBP-C phosphorylation patterns that are validated experimentally. Our results suggest non-redundant roles for PP1 and PP2A under both physiological and heart failure conditions, and emphasize the importance of phosphatases for cMyBP-C regulation.

Selection pressure by specialist and generalist insect herbivores leads to optimal constitutive plant defense. A mathematical model

Brassicaceae plants have the glucosinolate-myrosinase defense system, jointly active against herbivory. However, constitutive glucosinolate (GLS) defense is observed to occur at levels that do not deter all insects from feeding. That prompts the question of why Brassicaceae plants have not evolved a higher constitutive defense. The answer may lie in the contrasting relationship between plant defense and host plant preference of specialist and generalist herbivores. GLS content increases a plant's susceptibility to specialist insects. In contrast, generalists are deterred by the plant GLSs. Although GLSs can attract the natural enemies (predators and parasitoids) of these herbivores, enemies can reduce herbivore pressure to some extent only. So, plants can be overrun by specialists if GLS content is too high, whereas generalists can invade the plants if it is too low. Therefore, an optimal constitutive plant defense can minimize the overall herbivore pressure. To explain the optimal defense theoretically, we model the contrasting host selection behavior of insect herbivores and the emergence of their natural enemies by non-autonomous ordinary differential equations, where the independent variable is the plant GLS concentration. From the model, we quantify the optimal amount of GLSs, which minimizes total herbivore (specialists and generalists) pressure. That quite successfully explains the evolution of constitutive defense in plants from the perspective of optimality theory.

Immobilization of Lipases on Chitosan Hydrogels Improves Their Stability in the Presence of the Products of Triglyceride Oxidation

A significant bottleneck for the industrial application of lipases stems from their poor stability in the presence of commercial triglycerides. This is mainly due to the inactivating effect of the products of triglyceride oxidation (PTO), which are usually produced when oils and fats, being imported from far countries, are stored for long periods. In this study, the immobilization of a lipase from Candida rugosa on chitosan hydrogels has been carried out following two alternative approaches based on the enzyme adsorption and entrapment to increase the lipase stability under the operating conditions that are typical of oleochemical transformations. The effect of model compounds representing different classes of PTO on a lipase has been studied to optimize the enzyme immobilization method. Particular attention has been devoted to the characterization of the inactivating effect of PTO in nonaqueous media, which are adopted for most industrial applications of lipases.

Kinetic analyses of two-steps oxidation from l-tyrosine to l-dopaquinone with tyrosinase by capillary electrophoresis/dynamic frontal analysis

Tyrosinase catalyzes the oxidation of l-tyrosine in two stages to produce l-dopa and l-dopaquinone stepwise, and l-dopaquinone is subsequently converted to dopachrome. Most of the conventional analyses subjected only one-step reaction from l-tyrosine to l-dopa or from l-dopa to l-dopaquinone. In this study, kinetic analyses of two-steps oxidation of l-tyrosine with tyrosinase were made by capillary electrophoresis/dynamic frontal analysis (CE/DFA). When l-dopa was introduced into a capillary as a sample plug in a CE/DFA format, the enzymatic oxidation continuously occurred during the electrophoresis, and the product l-dopaquinone was subsequently converted to dopachrome which was detected as a plateau signal. A Michaelis-Menten constant of the second-step kinetic reaction, K_m,Do, was determined as 0.45 ± 0.03 mmol L^-1. In the analysis of the first-step kinetic reaction from l-tyrosine to l-dopa, l-dopa was not resolved by CE/DFA because both l-tyrosine and l-dopa are electrically neutral. The l-dopa formed and co-migrated at the l-tyrosine zone was calibrated beforehand with the final product of dopachrome detected as a plateau signal. Constantly formed l-dopa was successfully detected as a plateau signal of dopachrome, and a Michaelis-Menten constant of K_m,Ty was also determined as 0.061 ± 0.009 mmol L^-1 by the CE/DFA. CE/DFA is applicable to two-steps enzymatic reactions.

Simulations of radioiodine exposure and protective thyroid blocking in a new biokinetic model of the mother-fetus unit at different pregnancy ages

In the case of nuclear incidents, radioiodine may be released. After incorporation, it accumulates in the thyroid and enhances the risk of thyroidal dysfunctions and cancer occurrence by internal irradiation. Pregnant women and children are particularly vulnerable. Therefore, thyroidal protection by administering a large dose of stable (non-radioactive) iodine, blocking radioiodide uptake into the gland, is essential in these subpopulations. However, a quantitative estimation of the protection conferred to the maternal and fetal thyroids in the different stages of pregnancy is difficult. We departed from an established biokinetic model for radioiodine in pregnancy using first-order kinetics. As the uptake of iodide into the thyroid and several other tissues is mediated by a saturable active transport, we integrated an uptake mechanism described by a Michaelis–Menten kinetic. This permits simulating the competition between stable and radioactive iodide at the membrane carrier site, one of the protective mechanisms. The Wollf–Chaikoff effect, as the other protective mechanism, was simulated by adding a total net uptake block for iodide into the thyroid, becoming active when the gland is saturated with iodine. The model’s validity was confirmed by comparing predicted values with results from other models and sparse empirical data. According to our model, in the case of radioiodine exposure without thyroid blocking, the thyroid equivalent dose in the maternal gland increases about 45% within the first weeks of pregnancy to remain in the same range until term. Beginning in the 12th pregnancy week, the equivalent dose in the fetal thyroid disproportionately increases over time and amounts to three times the dose of the maternal gland at term. The maternal and fetal glands’ protection increases concomitantly with the amount of stable iodine administered to the mother simultaneously with acute radioiodine exposure. The dose–effect curves reflecting the combined thyroidal protection by the competition at the membrane carrier site and the Wolff–Chaikoff effect in the mother are characterized by a mean effective dose (ED₅₀) of roughly 1.5 mg all over pregnancy. In the case of the fetal thyroid, the mean effective doses for thyroid blocking, taking into account only the competition at the carrier site are numerically lower than in the mother. Taking into account additionally the Wolff–Chaikoff effect, the dose–effect curves for thyroidal protection in the fetus show a shift to the left over time, with a mean effective dose of 12.9 mg in the 12th week of pregnancy decreasing to 0.5 mg at term. In any case, according to our model, the usually recommended dose of 100 mg stable iodine given at the time of acute radioiodine exposure confers a very high level of thyroidal protection to the maternal and fetal glands over pregnancy. For ethical reasons, the possibilities of experimental studies on thyroid blocking in pregnant women are extremely limited. Furthermore, results from animal studies are associated with the uncertainties related to the translation of the data to humans. Thus model-based simulations may be a valuable tool for better insight into the efficacy of thyroidal protection and improve preparedness planning for uncommon nuclear or radiological emergencies.

Mechanistically Coupled PK (MCPK) Model to Describe Enzyme Induction and Occupancy Dependent DDI of Dabrafenib Metabolism

Dabrafenib inhibits the cell proliferation of metastatic melanoma with the oncogenic BRAF(V600)-mutation. However, dabrafenib monotherapy is associated with pERK reactivation, drug resistance, and consequential relapse. A clinical drug-dose determination study shows increased pERK levels upon daily administration of more than 300 mg dabrafenib. To clarify whether such elevated drug concentrations could be reached by long-term drug accumulation, we mechanistically coupled the pharmacokinetics (MCPK) of dabrafenib and its metabolites. The MCPK model is qualitatively based on in vitro and quantitatively on clinical data to describe occupancy-dependent CYP3A4 enzyme induction, accumulation, and drug–drug interaction mechanisms. The prediction suggests an eight-fold increase in the steady-state concentration of potent desmethyl-dabrafenib and its inactive precursor carboxy-dabrafenib within four weeks upon 150 mg b.d. dabrafenib. While it is generally assumed that a higher dose is not critical, we found experimentally that a high physiological dabrafenib concentration fails to induce cell death in embedded 451LU melanoma spheroids.

Modeling principles of protective thyroid blocking

PURPOSE

In the case of a nuclear incident, the release of radioiodine must be expected. Radioiodine accumulates in the thyroid and by irradiation enhances the risk of cancer. Large doses of stable (non-radioactive) iodine may inhibit radioiodine accumulation and protect the thyroid ('thyroid blocking'). Protection is based on a competition at the active carrier site in the cellular membrane and an additional temporary inhibition of the organification of iodide (Wolff-Chaikoff effect). Alternatively, other agents like e.g. perchlorate that compete with iodide for the uptake into the thyrocytes may also confer thyroidal protection against radioiodine exposure.Biokinetic models for radioiodine mostly describe exchanges between compartments by first order kinetics. This leads to correct predictions only for low (radio)iodide concentrations. These models are not suited to describe the kinetics of iodine if administered at the dosages recommended for thyroid blocking and moreover does not permit to simulate either the protective competition mechanism at the membrane or the Wolff-Chaikoff effect. Models adapted for this purpose must be used. Such models may use a mathematical relation between the serum iodide concentration and a relative uptake suppression or a dependent rate constant determining total thyroidal radioiodine accumulation. Alternatively, the thyroidal uptake rate constant may be modeled as a function of the total iodine content of the gland relative to a saturation amount. Newer models integrate a carrier-mechanism described by Michalis-Menten kinetics in the membrane and in analogy to enzyme kinetics apply the rate law for monomolecular irreversible enzyme reactions with competing substrates to model the competition mechanism. An additional total iodide uptake block, independent on competition but limited in time, is used to simulate the Wolff-Chaikoff effect.

CONCLUSION

The selection of the best model depends on the issue to be studied. Most models cannot quantify the relative contributions of the competition mechanism at the membrane and the Wolff-Chaikoff effect. This makes it impossible or exceedingly difficult to simulate prolonged radioiodine exposure and the effect of repetitive administrations of stable iodine. The newer thyroid blocking models with a separate modeling of competition and Wolff-Chaikoff effect allow better quantitative mechanistic insights and offer the possibility to simulate complex radioiodine exposure scenarios and various protective dosage schemes of stable iodine relatively easily. Moreover, they permit to study the protective effects of other competitors at the membrane carrier site, like e.g. perchlorate, and to draw conclusions on their protective efficacy in comparison to stable iodine.

Finding Ways to Relax: A Revisionistic Analysis of the Chemistry of E. coli GTP Cyclohydrolase II

Guanosine triphosphate (GTP) cyclohydrolase II (RibA) is one of three enzymes that hydrolytically cleave the C8-N9 bond of the GTP guanine. RibA also catalyzes a subsequent hydrolytic attack at the base liberating formate and in addition cleaves the α-β phosphodiester bond of the triphosphate to form pyrophosphate (PPi). These hydrolytic reactions are promoted by tandem active-site metal ions, zinc and magnesium, that respectively function at the GTP guanine and triphosphate moieties. The RibA reaction is part of riboflavin biosynthesis and forms 2,5-diamino-6-β-pyrimidinone 5'-phosphate, an exocyclic pyrimidine nucleotide that ultimately forms the pyrimidine ring of the isoalloxazine of riboflavin. The stoichiometry of the RibA reaction was defined in the study that first identified this activity in Escherichia coli (Foor, F., Brown, G. M. J. Biol. Chem., 1975, 250, 9, 3545-3551) and has not been quantitatively evaluated in subsequent works. Using primarily transient state approaches we examined the interaction of RibA from E. coli with the GTP, inosine triphosphate, and PPi. Our data indicate that PPi is a slow substrate for RibA that is cleaved to form two phosphate ions (Pi). A combination of real-time enzymatically coupled Pi reporter assays and end-point ³¹P NMR revealed that Pi is formed at a catalytically relevant rate in the native reaction of RibA with GTP, redefining the reaction stoichiometry. Furthermore, our data indicate that both PPi and GTP stimulate conformational changes prior to hydrolytic chemistry, and we conclude that the cleavage of PPi bound as a substrate or an intermediate state results in conformational relaxation.

Molecular noise filtering in the β-adrenergic signaling network by phospholamban pentamers

Phospholamban (PLN) is an important regulator of cardiac calcium handling due to its ability to inhibit the calcium ATPase SERCA. β-Adrenergic stimulation reverses SERCA inhibition via PLN phosphorylation and facilitates fast calcium reuptake. PLN also forms pentamers whose physiological significance has remained elusive. Using mathematical modeling combined with biochemical and cell biological experiments, we show that pentamers regulate both the dynamics and steady-state levels of monomer phosphorylation. Substrate competition by pentamers and a feed-forward loop involving inhibitor-1 can delay monomer phosphorylation by protein kinase A (PKA), whereas cooperative pentamer dephosphorylation enables bistable PLN steady-state phosphorylation. Simulations show that phosphorylation delay and bistability act as complementary filters that reduce the effect of random fluctuations in PKA activity, thereby ensuring consistent monomer phosphorylation and SERCA activity despite noisy upstream signals. Preliminary analyses suggest that the PLN mutation R14del could impair noise filtering, offering a new perspective on how this mutation causes cardiac arrhythmias.

A comparison of thyroidal protection by iodine and perchlorate against radioiodine exposure in Caucasians and Japanese

Radioactive iodine released in nuclear accidents may accumulate in the thyroid and by irradiation enhances the risk of cancer. Radioiodine uptake into the gland can be inhibited by large doses of stable iodine or perchlorate. Nutritional iodine daily intake may impact thyroid physiology, so that radiological doses absorbed by the thyroid as well as thyroid blocking efficacy may differ in Japanese with a very rich iodine diet compared to Caucasians. Based on established biokinetic-dosimetric models for the thyroid, we derived the parameters for Caucasians and Japanese to quantitatively compare the effects of radioiodine exposure and the protective efficacy of thyroid blocking by stable iodine at the officially recommended dosages (100 mg in Germany, 76 mg in Japan) or perchlorate. The maximum transport capacity for iodine uptake into the thyroid is lower in Japanese compared to Caucasians. For the same radioiodine exposure pattern, the radiological equivalent thyroid dose is substantially lower in Japanese in the absence of thyroid blocking treatments. In the case of acute radioiodine exposure, stable iodine is less potent in Japanese (ED50 = 41.6 mg) than in Caucasians (ED50 = 2.7 mg) and confers less thyroid protection at the recommended dosages because of a delayed responsiveness to iodine saturation of the gland (Wolff-Chaikoff effect). Perchlorate (ED50 = 10 mg in Caucasians) at a dose of 1000 mg has roughly the same thyroid blocking effect as 100 mg iodine in Caucasians, whereas it confers a much better protection than 76 mg iodine in Japanese. For prolonged exposures, a single dose of iodine offer substantially lower protection than after acute radioiodine exposure in both groups. Repetitive daily iodine administrations improve efficacy without reaching levels after acute radioiodine exposure and achieve only slightly better protection in Japanese than in Caucasians. However, in the case of continuous radioiodine exposure, daily doses of 1000 mg perchlorate achieve a high protective efficacy in Caucasians as well as Japanese (> 0.98). In Caucasians, iodine (100 mg) and perchlorate (1000 mg) at the recommended dosages seem alternatives in case of acute radioiodine exposure, whereas perchlorate has a higher protective efficacy in the case of longer lasting radioiodine exposures. In Japanese, considering protective efficacy, preference should be given to perchlorate in acute as well as prolonged radioiodine exposure scenarios.

Resource uptake and the evolution of moderately efficient enzymes

Enzymes speed up reactions that would otherwise be too slow to sustain the metabolism of selfreplicators. Yet, most enzymes seem only moderately efficient, exhibiting kinetic parameters orders of magnitude lower than their expected physically achievable maxima and spanning over surprisingly large ranges of values. Here, we question how these parameters evolve using a mechanistic model where enzyme efficiency is a key component of individual competition for resources. We show that kinetic parameters are under strong directional selection only up to a point, above which enzymes appear to evolve under near-neutrality, thereby confirming the qualitative observation of other modeling approaches. While the existence of a large fitness plateau could potentially explain the extensive variation in enzyme features reported, we show using a population genetics model that such a widespread distribution is an unlikely outcome of evolution on a common landscape, as mutation–selection–drift balance occupy a narrow area even when very moderate biases towards lower efficiency are considered. Instead, differences in the evolutionary context encountered by each enzyme should be involved, such that each evolves on an individual, unique landscape. Our results point to drift and effective population size playing an important role, along with the kinetics of nutrient transporters, the tolerance to high concentrations of intermediate metabolites, and the reversibility of reactions. Enzyme concentration also shapes selection on kinetic parameters, but we show that the joint evolution of concentration and efficiency does not yield extensive variance in evolutionary outcomes when documented costs to protein expression are applied.

Systems Biology Modeling of the Complement System Under Immune Susceptible Pathogens

The complement system is assembled from a network of proteins that function to bring about the first line of defense of the body against invading pathogens. However, complement deficiencies or invasive pathogens can hijack complement to subsequently increase susceptibility of the body to infections. Moreover, invasive pathogens are increasingly becoming resistant to the currently available therapies. Hence, it is important to gain insights into the highly dynamic interaction between complement and invading microbes in the frontlines of immunity. Here, we developed a mathematical model of the complement system composed of 670 ordinary differential equations with 328 kinetic parameters, which describes all three complement pathways (alternative, classical, and lectin) and includes description of mannose-binding lectin, collectins, ficolins, factor H-related proteins, immunoglobulin M, and pentraxins. Additionally, we incorporate two pathogens: (type 1) complement susceptible pathogen and (type 2) Neisseria meningitidis located in either nasopharynx or bloodstream. In both cases, we generate time profiles of the pathogen surface occupied by complement components and the membrane attack complex (MAC). Our model shows both pathogen types in bloodstream are saturated by complement proteins, whereas MACs occupy <<1.0% of the pathogen surface. Conversely, the MAC production in nasopharynx occupies about 1.5-10% of the total N. meningitidis surface, thus making nasal MAC levels at least about eight orders of magnitude higher. Altogether, we predict complement-imbalance, favoring overactivation, is associated with nasopharynx homeostasis. Conversely, orientating toward complement-balance may cause disruption to the nasopharynx homeostasis. Thus, for sporadic meningococcal disease, our model predicts rising nasal levels of complement regulators as early infection biomarkers.

A comparison of thyroidal protection by stable iodine or perchlorate in the case of acute or prolonged radioiodine exposure

In the case of a nuclear power plant accident, repetitive/prolonged radioiodine release may occur. Radioiodine accumulates in the thyroid and by irradiation enhances the risk of cancer. Large doses of non-radioactive iodine may protect the thyroid by inhibiting radioiodine uptake into the gland (iodine blockade). Protection is based on a competition at the active carrier site in the cellular membrane and the Wolff-Chaikoff effect, the latter being, however, only transient (24-48 h). Perchlorate may alternatively provide protection by a carrier competition mechanism only. Perchlorate has, however, a stronger affinity to the carrier than iodide. Based on an established biokinetic-dosimetric model developed to study iodine blockade, and after its extension to describe perchlorate pharmacokinetics and the inhibition of iodine transport through the carrier, we computed the protective efficacies that can be achieved by stable iodine or perchlorate in the case of an acute or prolonged radioiodine exposure. In the case of acute radioiodine exposure, perchlorate is less potent than stable iodine considering its ED50. A dose of 100 mg stable iodine has roughly the same protective efficacy as 1000 mg perchlorate. For prolonged exposures, single doses of protective agents, whether stable iodine or perchlorate, offer substantially lower protection than after acute radioiodine exposure, and thus repetitive administrations seem necessary. In case of prolonged exposure, the higher affinity of perchlorate for the carrier in combination with the fading Wolff-Chaikoff effect of iodine confers perchlorate a higher protective efficacy compared to stable iodine. Taking into account the frequency and seriousness of adverse effects, iodine and perchlorate at equieffective dosages seem to be alternatives in case of short-term acute radioiodine exposure, whereas preference should be given to perchlorate in view of its higher protective efficacy in the case of longer lasting radioiodine exposures.

Homo-Oligomerisation in Signal Transduction: Dynamics, Homeostasis, Ultrasensitivity, Bistability

Homo-oligomerisation of proteins is a ubiquitous phenomenon whose exact role remains unclear in many cases. To identify novel functions, this paper provides an exploration of general dynamical mathematical models of homo-oligomerisation. Simulation and analysis of these models show that homo-oligomerisation on its own allows for a remarkable variety of complex dynamic and steady state regulatory behaviour such as transient overshoots or homeostatic control of monomer concentration. If post-translational modifications are considered, however, conventional mass action kinetics lead to thermodynamic inconsistencies due to asymmetric combinatorial expansion of reaction routes. Introducing a conservation principle to balance rate equations re-establishes thermodynamic consistency. Using such balanced models it is shown that oligomerisation can lead to bistability by enabling pseudo-multisite modification and kinetic pseudo-cooperativity via multi-enzyme regulation, thereby constituting a novel motif for bistable modification reactions. Due to these potential signal processing capabilities, homo-oligomerisation could play far more versatile roles in signal transduction than previously appreciated.

ENaC regulation by phospholipids and DGK explained through mathematical modeling

Cystic fibrosis is a condition caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR). It is also thought to increase the activity of epithelial sodium channels (ENaC). The altered function of these ion channels is one of the causes of the thick dehydrated mucus that characterizes the disease and is partially responsible for recurrent pulmonary infections and inflammation events that ultimately destroy the lungs of affected subjects. Phosphoinositides are signaling lipids that regulate numerous cellular processes and membrane proteins, including ENaC. Inhibition of diacylglycerol kinase (DGK), an enzyme of the phosphoinositide pathway, reduces ENaC function. We propose a computational analysis that is based on the combination of two existing mathematical models: one representing the dynamics of phosphoinositides and the other explaining how phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) influences ENaC activity and, consequently, airway surface liquid. This integrated model permits, for the first time, a detailed assessment of the intricate interactions between DGK and ENaC and is consistent with available literature data. In particular, the computational approach allows comparisons of two competing hypotheses regarding the regulation of ENaC. The results strongly suggest that the regulation of ENaC is primarily exerted through the control of PI(4,5)P2 production by type-I phosphatidylinositol-4-phosphate 5-kinase (PIP5KI), which in turn is controlled by phosphatidic acid (PA), the product of the DGK reaction.

Kinetic analysis of substrate competition in enzymatic reactions with β-D-galactosidase by capillary electrophoresis / dynamic frontal analysis

Competitive inhibition between two substrates with an enzyme is investigated by capillary electrophoresis/dynamic frontal analysis (CE/DFA). Enzymatic hydrolyses of o-nitrophenyl β-D-galactopyranoside and p-nitrophenyl β-D-galactopyranoside with β-D-galactosidase were examined as a model competitive reaction. A sample solution containing the two substrates was injected into a capillary filled with a separation buffer containing an enzyme. Enzymatic hydrolysis occurred during the electrophoresis, and the products of o-nitrophenol and p-nitrophenol were continuously formed and resolved from the sample zone. Two-steps plateau signal was detected with the two-substrate solutions based on the difference in the effective electrophoretic mobility of o-nitrophenol and p-nitrophenol. Michaelis-Menten constants and inhibition constants were determined with the plateau heights. Usefulness of CE/DFA on competitive inhibition analysis is demonstrated in this study.

Computational Modeling of Anthocyanin Pathway Evolution: Biases, Hotspots, and Trade-offs

The alteration of metabolic pathways is a common mechanism underlying the evolution of new phenotypes. Flower color is a striking example of the importance of metabolic evolution in a complex phenotype, wherein shifts in the activity of the underlying pathway lead to a wide range of pigments. Although experimental work has identified common classes of mutations responsible for transitions among colors, we lack a unifying model that relates pathway function and activity to the evolution of distinct pigment phenotypes. One challenge in creating such a model is the branching structure of pigment pathways, which may lead to evolutionary trade-offs due to competition for shared substrates. In order to predict the effects of shifts in enzyme function and activity on pigment production, we created a simple kinetic model of a major plant pigmentation pathway: the anthocyanin pathway. This model describes the production of the three classes of blue, purple, and red anthocyanin pigments, and accordingly, includes multiple branches and substrate competition. We first studied the general behavior of this model using a naïve set of parameters. We then stochastically evolved the pathway toward a defined optimum and analyzed the patterns of fixed mutations. This approach allowed us to quantify the probability density of trajectories through pathway state space and identify the types and number of changes. Finally, we examined whether our simulated results qualitatively align with experimental observations, i.e., the predominance of mutations which change color by altering the function of branching genes in the pathway. These analyses provide a theoretical framework that can be used to predict the consequences of new mutations in terms of both pigment phenotypes and pleiotropic effects.

Identification of evolutionary and kinetic drivers of NAD-dependent signaling

Nicotinamide adenine dinucleotide (NAD) provides an important link between metabolism and signal transduction and has emerged as central hub between bioenergetics and all major cellular events. NAD-dependent signaling (e.g., by sirtuins and poly-adenosine diphosphate [ADP] ribose polymerases [PARPs]) consumes considerable amounts of NAD. To maintain physiological functions, NAD consumption and biosynthesis need to be carefully balanced. Using extensive phylogenetic analyses, mathematical modeling of NAD metabolism, and experimental verification, we show that the diversification of NAD-dependent signaling in vertebrates depended on 3 critical evolutionary events: 1) the transition of NAD biosynthesis to exclusive usage of nicotinamide phosphoribosyltransferase (NamPT); 2) the occurrence of nicotinamide N-methyltransferase (NNMT), which diverts nicotinamide (Nam) from recycling into NAD, preventing Nam accumulation and inhibition of NAD-dependent signaling reactions; and 3) structural adaptation of NamPT, providing an unusually high affinity toward Nam, necessary to maintain NAD levels. Our results reveal an unexpected coevolution and kinetic interplay between NNMT and NamPT that enables extensive NAD signaling. This has implications for therapeutic strategies of NAD supplementation and the use of NNMT or NamPT inhibitors in disease treatment.

Substrate complex competition is a regulatory motif that allows NFκB RelA to license but not amplify NFκB RelB

Signaling pathways often share molecular components, tying the activity of one pathway to the functioning of another. In the NFκB signaling system, distinct kinases mediate inflammatory and developmental signaling via RelA and RelB, respectively. Although the substrates of the developmental, so-called noncanonical, pathway are induced by inflammatory/canonical signaling, crosstalk is limited. Through dynamical systems modeling, we identified the underlying regulatory mechanism. We found that as the substrate of the noncanonical kinase NIK, the nfkb2 gene product p100, transitions from a monomer to a multimeric complex, it may compete with and inhibit p100 processing to the active p52. Although multimeric complexes of p100 (IκBδ) are known to inhibit preexisting RelA:p50 through sequestration, here we report that p100 complexes can inhibit the enzymatic formation of RelB:p52. We show that the dose-response systems properties of this complex substrate competition motif are poorly accounted for by standard Michaelis-Menten kinetics, but require more detailed mass action formulations. In sum, although tonic inflammatory signaling is required for adequate expression of the noncanonical pathway precursors, the substrate complex competition motif identified here can prevent amplification of the active RelB:p52 dimer in elevated inflammatory conditions to ensure reliable RelB-dependent developmental signaling independent of inflammatory context.

A mechanism for bistability in glycosylation

Glycosyltransferases are a class of enzymes that catalyse the posttranslational modification of proteins to produce a large number of glycoconjugate acceptors from a limited number of nucleotide-sugar donors. The products of one glycosyltransferase can be the substrates of several other enzymes, causing a combinatorial explosion in the number of possible glycan products. The kinetic behaviour of systems where multiple acceptor substrates compete for a single enzyme is presented, and the case in which high concentrations of an acceptor substrate are inhibitory as a result of abortive complex formation, is shown to result in non-Michaelian kinetics that can lead to bistability in an open system. A kinetic mechanism is proposed that is consistent with the available experimental evidence and provides a possible explanation for conflicting observations on the β-1,4-galactosyltransferases. Abrupt switching between steady states in networks of glycosyltransferase-catalysed reactions may account for the observed changes in glycosyl-epitopes in cancer cells.

Simple, selective and fast detection of acrylamide based on glutathione S-transferase

Acrylamide (AA) is a toxic compound formed in thermally prepared foods by Maillard reaction. Besides foods, AA may be found in cosmetic products as an impurity of the widely-used non-toxic polyacrylamide. We present a novel, fast and selective detection method based on the amperometric monitoring of the coupling reaction between reduced glutathione (GSH) and AA catalyzed by glutathione S-transferase (GST) to produce an electrochemically inactive compound. We have used electrodes modified with cobalt-phthalocyanine to monitor the decrease of GHS concentration at +300 mV. Our system is simple, does not require supplementary substrates such as 1-chloro-2,4-dinitrobenzene (CDNB) nor have disadvantageous competitive kinetics characteristic to inhibition like signals. Using the optimum concentration of 100 μM GSH we have obtained a linear calibration graph from 7 to 50 μM AA and a limit of detection of 5 μM AA. The method is not affected by interfering compounds usually found in foods and was applied for real sample analysis.

Data Reduction Approaches for Dissecting Transcriptional Effects on Metabolism

The availability of high-throughput data from transcriptomics and metabolomics technologies provides the opportunity to characterize the transcriptional effects on metabolism. Here we propose and evaluate two computational approaches rooted in data reduction techniques to identify and categorize transcriptional effects on metabolism by combining data on gene expression and metabolite levels. The approaches determine the partial correlation between two metabolite data profiles upon control of given principal components extracted from transcriptomics data profiles. Therefore, they allow us to investigate both data types with all features simultaneously without doing preselection of genes. The proposed approaches allow us to categorize the relation between pairs of metabolites as being under transcriptional or post-transcriptional regulation. The resulting classification is compared to existing literature and accumulated evidence about regulatory mechanism of reactions and pathways in the cases of Escherichia coli, Saccharomycies cerevisiae, and Arabidopsis thaliana.

PqsL uses reduced flavin to produce 2-hydroxylaminobenzoylacetate, a preferred PqsBC substrate in alkyl quinolone biosynthesis in Pseudomonas aeruginosa

Alkyl hydroxyquinoline N-oxides (AQNOs) are antibiotic compounds produced by the opportunistic bacterial pathogen Pseudomonas aeruginosa They are products of the alkyl quinolone (AQ) biosynthetic pathway, which also generates the quorum-sensing molecules 2-heptyl-4(1H)-quinolone (HHQ) and 2-heptyl-3-hydroxy-4(1H)-quinolone (PQS). Although the enzymatic synthesis of HHQ and PQS had been elucidated, the route by which AQNOs are synthesized remained elusive. Here, we report on PqsL, the key enzyme for AQNO production, which structurally resembles class A flavoprotein monooxygenases such as p-hydroxybenzoate 3-hydroxylase (pHBH) and 3-hydroxybenzoate 6-hydroxylase. However, we found that unlike related enzymes, PqsL hydroxylates a primary aromatic amine group, and it does not use NAD(P)H as cosubstrate, but unexpectedly required reduced flavin as electron donor. We also observed that PqsL is active toward 2-aminobenzoylacetate (2-ABA), the central intermediate of the AQ pathway, and forms the unstable compound 2-hydroxylaminobenzoylacetate, which was preferred over 2-ABA as substrate of the downstream enzyme PqsBC. In vitro reconstitution of the PqsL/PqsBC reaction was feasible by using the FAD reductase HpaC, and we noted that the AQ:AQNO ratio is increased in an hpaC-deletion mutant of P. aeruginosa PAO1 compared with the ratio in the WT strain. A structural comparison with pHBH, the model enzyme of class A flavoprotein monooxygenases, revealed that structural features associated with NAD(P)H binding are missing in PqsL. Our study completes the AQNO biosynthetic pathway in P. aeruginosa, indicating that PqsL produces the unstable product 2-hydroxylaminobenzoylacetate from 2-ABA and depends on free reduced flavin as electron donor instead of NAD(P)H.

Developmental Plasticity of the Major Alkyl Cannabinoid Chemotypes in a Diverse Cannabis Genetic Resource Collection

Cannabis is a chemically diverse domesticated plant genus which produces a unique class of biologically active secondary metabolites referred to as cannabinoids. The affinity and selectivity of cannabinoids to targets of the human endocannabinoid system depend on alkyl side chain length, and these structural-activity relationships can be utilized for the development of novel therapeutics. Accurate early screening of germplasm has the potential to accelerate selection of chemical phenotypes (chemotypes) for pharmacological exploitation. However, limited attempts have been made to characterize the plasticity of alkyl cannabinoid composition in different plant tissues and throughout development. A chemotypic diversity panel comprised of 99 individuals from 20 Cannabis populations sourced from the Ecofibre Global Germplasm Collection (ecofibre.com.au and anandahemp.com) was used to examine alkyl cannabinoid variation across vegetative, flowering and maturation stages. A wide range of di-/tri-cyclic as well as C₃-/C₅-alkyl cannabinoid composition was observed between plants. Chemotype at the vegetative and flowering stages was found to be predictive of chemotype at maturation, indicating a low level of plasticity in cannabinoid composition. Chemometric cluster analysis based on composition data from all three developmental stages categorized alkyl cannabinoid chemotypes into three classes. Our results suggest that more extensive chemical and genetic characterization of the Cannabis genepool could facilitate the metabolic engineering of alkyl cannabinoid chemotypes.

Stoichiometry and kinetics of single and mixed substrate uptake in Aspergillus niger

In its natural environment, the filamentous fungus Aspergillus niger grows on decaying fruits and plant material, thereby enzymatically degrading the lignocellulosic constituents (lignin, cellulose, hemicellulose, and pectin) into a mixture of mono- and oligosaccharides. To investigate the kinetics and stoichiometry of growth of this fungus on lignocellulosic sugars, we carried out batch cultivations on six representative monosaccharides (glucose, xylose, mannose, rhamnose, arabinose, and galacturonic acid) and a mixture of these. Growth on these substrates was characterized in terms of biomass yields, oxygen/biomass ratios, and specific conversion rates. Interestingly, in combination, some of the carbon sources were consumed simultaneously and some sequentially. With a previously developed protocol, a sequential chemostat cultivation experiment was performed on a feed mixture of the six substrates. We found that the uptake of glucose, xylose, and mannose could be described with a Michaelis–Menten-type kinetics; however, these carbon sources seem to be competing for the same transport systems, while the uptake of arabinose, galacturonic acid, and rhamnose appeared to be repressed by the presence of other substrates.

Two faces of competition: target-mediated reverse signalling in microRNA and mitogen-activated protein kinase regulatory networks

Biomolecular regulatory networks are organised around hubs, which can interact with a large number of targets. These targets compete with each other for access to their common hubs, but whether the effect of this competition would be significant in magnitude and in function is not clear. In this review, the authors discuss recent in vivo studies that analysed the system level retroactive effects induced by target competition in microRNA and mitogen-activated protein kinase regulatory networks. The results of these studies suggest that downstream targets can regulate the overall state of their upstream regulators, and thus cannot be ignored in analysing biomolecular networks.

Active-Site Flexibility and Substrate Specificity in a Bacterial Virulence Factor: Crystallographic Snapshots of an Epoxide Hydrolase

Pseudomonas aeruginosa secretes an epoxide hydrolase with catalytic activity that triggers degradation of the cystic fibrosis transmembrane conductance regulator (CFTR) and perturbs other host defense networks. Targets of this CFTR inhibitory factor (Cif) are largely unknown, but include an epoxy-fatty acid. In this class of signaling molecules, chirality can be an important determinant of physiological output and potency. Here we explore the active-site chemistry of this two-step α/β-hydrolase and its implications for an emerging class of virulence enzymes. In combination with hydrolysis data, crystal structures of 15 trapped hydroxyalkyl-enzyme intermediates reveal the stereochemical basis of Cif's substrate specificity, as well as its regioisomeric and enantiomeric preferences. The structures also reveal distinct sets of conformational changes that enable the active site to expand dramatically in two directions, accommodating a surprising array of potential physiological epoxide targets. These new substrates may contribute to Cif's diverse effects in vivo, and thus to the success of P. aeruginosa and other pathogens during infection.

Glucosylceramide transferase in Giardia preferentially catalyzes the synthesis of galactosylceramide during encystation

The stage differentiation from trophozoite to cyst (i.e., encystation) is an essential step for Giardia to survive outside its human host and spread the infection via the fecal-oral route. We have previously shown that Giardia expresses glucosylceramide transferase 1 (GlcT1) enzyme, the activity of which is elevated during encystation. We have also reported that blocking the activity of gGlcT1 interferes with the biogenesis of encystation-specific vesicles (ESVs) and cyst viability in Giardia. To further understand the role of this enzyme and how it regulates encystation, we overexpressed, knocked down, and rescued the giardial GlcT1 (gGlcT1) gene and measured its enzymatic activity in live parasites as well as in isolated membrane fractions using NBD-ceramide and UDP-glucose or UDP-galactose. We observed that gGlcT1 is able to catalyze the synthesis of both glucosylceramide (GlcCer) and galactosylceramide (GalCer), however the synthesis of GalCer is 2-3 fold higher than of GlcCer. Although both activities follow Michaelis-Menten kinetics, the bindings of UDP-glucose and UDP-galactose with the enzyme appear to be non-competitive and independent of each other. The modulation of gGlcT1 synthesis concomitantly influenced the expression cyst-wall protein (CWP) and overall encystation. We propose that gGlcT1 is a unique enzyme and that Giardia uses this enzyme to synthesize both GlcCer and GalCer to facilitate the process of encystation/cyst production.

Stochastic modeling of biochemical systems with multistep reactions using state-dependent time delay

To deal with the growing scale of molecular systems, sophisticated modelling techniques have been designed in recent years to reduce the complexity of mathematical models. Among them, a widely used approach is delayed reaction for simplifying multistep reactions. However, recent research results suggest that a delayed reaction with constant time delay is unable to describe multistep reactions accurately. To address this issue, we propose a novel approach using state-dependent time delay to approximate multistep reactions. We first use stochastic simulations to calculate time delay arising from multistep reactions exactly. Then we design algorithms to calculate time delay based on system dynamics precisely. To demonstrate the power of proposed method, two processes of mRNA degradation are used to investigate the function of time delay in determining system dynamics. In addition, a multistep pathway of metabolic synthesis is used to explore the potential of the proposed method to simplify multistep reactions with nonlinear reaction rates. Simulation results suggest that the state-dependent time delay is a promising and accurate approach to reduce model complexity and decrease the number of unknown parameters in the models.

Sjögren's syndrome-associated microRNAs in CD14(+) monocytes unveils targeted TGFβ signaling

Background

Sjögren’s syndrome (SjS) monocytes have a pro-inflammatory phenotype, which may influence SjS pathogenesis. MicroRNAs (miRNAs) are small endogenously expressed molecules that can inhibit protein expression of their targeted genes and have important functions in regulating cell signaling responses. We profiled miRNAs in SjS monocytes to identify a SjS-specific miRNA profile and determine the potential roles of miRNAs in SjS pathogenesis.

Methods

Total RNA was extracted from healthy control (HC, n = 10), SjS (n = 18), systemic lupus erythematosus (SLE, n = 10), and rheumatoid arthritis (RA, n = 10) peripheral blood CD14⁺ monocytes for miRNA microarray analysis. To validate select miRNAs from the microarray analysis, the original cohort and a new cohort of monocyte RNA samples from HC (n = 9), SjS (n = 12), SLE (n = 8), and RA (n = 9) patients were evaluated by quantitative reverse transcription (RT)-PCR. Functional predictions of differentially expressed miRNAs were determined through miRNA target prediction database analyses. Statistical analyses performed included one-way analysis of variance with Bonferroni post tests, linear regression, and receiver operating characteristic curve analyses.

Results

MiRNAs were predominantly upregulated in SjS monocytes in comparison with controls. Quantitative RT-PCR confirmations supported co-regulation of miR-34b-3p, miR-4701-5p, miR-609, miR-300, miR-3162-3p, and miR-877-3p in SjS monocytes (13/30, 43.3 %) in comparison with SLE (1/17, 5.8 %) and RA (1/18, 5.6 %). MiRNA-target pathway predictions identified SjS-associated miRNAs appear to preferentially target the canonical TGFβ signaling pathway as opposed to pro-inflammatory interleukin-12 and Toll-like receptor/NFkB pathways.

Conclusions

Our results underscore a novel underlying molecular mechanism where SjS-associated miRNAs may collectively suppress TGFβ signaling as opposed to pro-inflammatory interleukin-12 and Toll-like receptor/NFκB pathways in SjS pathogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-016-0987-0) contains supplementary material, which is available to authorized users.

Quantifying NAD(P)H production in the upper Entner-Doudoroff pathway from Pseudomonas putida KT2440

•

The first kinetic characterization of PputG6PDH-1 is presented.

•

The relative production of NADH and NADPH by PputG6PDH-1 is quantified.

•

The stoichiometric matrix of in silico metabolic models for Pseudomonas putida must be modified.

Despite the lack of biochemical information, all available in silico metabolic models of Pseudomonas putida KT2440 consider NADP as the only cofactor accepted by the glucose-6-phosphate dehydrogenases. Because the Entner–Doudoroff pathway is the main glycolytic route in this bacterium, determining how much NADH and NADPH are produced in the reaction catalyzed by these enzymes is very important for the correct interpretation of metabolic flux distributions. To determine the actual cofactor preference of the glucose-6-phosphate dehydrogenase encoded by the zwf-1 gene (PputG6PDH-1), the major isoform during growth on glucose, we purified this protein and studied its kinetic properties. Based on simple kinetic principles, we estimated the in vivo relative production of NADH and NADPH during the oxidation of glucose-6-phosphate (G6P). Contrary to the general assumption, our calculations showed that the reaction catalyzed by PputG6PDH-1 yields around 1/3 mol of NADPH and 2/3 mol of NADH per mol of oxidized G6P. Additionally, we obtained data suggesting that the reaction catalyzed by the 6-phosphogluconate dehydrogenase is active during growth on glucose, and it also produces NADH. These results indicate that the stoichiometric matrix of in silico models of P. putida KT2440 must be corrected and highlight the importance of considering the physiological concentrations of the involved metabolites to estimate the actual proportion of NADH and NADPH produced by a dehydrogenase.

Model of tryptophan metabolism, readily scalable using tissue-specific gene expression data

Tryptophan is utilized in various metabolic routes including protein synthesis, serotonin, and melatonin synthesis and the kynurenine pathway. Perturbations in these pathways have been associated with neurodegenerative diseases and cancer. Here we present a comprehensive kinetic model of the complex network of human tryptophan metabolism based upon existing kinetic data for all enzymatic conversions and transporters. By integrating tissue-specific expression data, modeling tryptophan metabolism in liver and brain returned intermediate metabolite concentrations in the physiological range. Sensitivity and metabolic control analyses identified expected key enzymes to govern fluxes in the branches of the network. Combining tissue-specific models revealed a considerable impact of the kynurenine pathway in liver on the concentrations of neuroactive derivatives in the brain. Moreover, using expression data from a cancer study predicted metabolite changes that resembled the experimental observations. We conclude that the combination of the kinetic model with expression data represents a powerful diagnostic tool to predict alterations in tryptophan metabolism. The model is readily scalable to include more tissues, thereby enabling assessment of organismal tryptophan metabolism in health and disease.