Search

In this study we demonstrated through analytic considerations and numerical studies that the mitochondrial fatty-acid β-oxidation can exhibit bistable-hysteresis behavior. In an experimentally validated computational model we identified a specific region in the parameter space in which two distinct stable and one unstable steady state could be attained with different fluxes. The two stable states were referred to as low-flux (disease) and high-flux (healthy) state. By a modular kinetic approach we traced the origin and causes of the bistability back to the distributive kinetics and the conservation of CoA, in particular in the last rounds of the β-oxidation. We then extended the model to investigate various interventions that may confer health benefits by activating the pathway, including (i) activation of the last enzyme MCKAT via its endogenous regulator p46-SHC protein, (ii) addition of a thioesterase (an acyl-CoA hydrolysing enzyme) as a safety valve, and (iii) concomitant activation of a number of upstream and downstream enzymes by short-chain fatty-acids (SCFA), metabolites that are produced from nutritional fibers in the gut. A high concentration of SCFAs, thioesterase activity, and inhibition of the p46Shc protein led to a disappearance of the bistability, leaving only the high-flux state. A better understanding of the switch behavior of the mitochondrial fatty-acid oxidation process between a low- and a high-flux state may lead to dietary and pharmacological intervention in the treatment or prevention of obesity and or non-alcoholic fatty-liver disease. Author summary: Obesity is a complex disease which is still poorly understood at the systemic level. Impaired capacity in fat oxidation may contribute to the development of obesity. In this manuscript using a detailed mitochondrial fatty acid oxidation computational model, we demonstrate that the oxidation of fat can exhibit bistability and hysteresis, implying that there is a risk that the pathway gets trapped in a stable state with low activity. We also identify the cause of bistability in the metabolic network and analyse which clinically relevant factors shift the pathway between the high-flux, healthy and low-flux, diseased state. [ABSTRACT FROM AUTHOR]

Mitochondrial fatty-acid beta-oxidation (mFAO) plays a central role in mammalian energy metabolism. Multiple severe diseases are associated with defects in this pathway. Its kinetic structure is characterized by a complex wiring of which the functional implications have hardly been explored. Repetitive cycles of reversible reactions, each cycle shortening the fatty acid by two carbon atoms, evoke competition between intermediates of different chain lengths for a common set of ‘promiscuous’ enzymes (enzymes with activity towards multiple substrates). In our validated kinetic model of the pathway, substrate overload causes a steep and detrimental flux decline. Here, we unravel the underlying mechanism and the role of enzyme promiscuity in it. Comparison of alternative model versions elucidated the role of promiscuity of individual enzymes. Promiscuity of the last enzyme of the pathway, medium-chain ketoacyl-CoA thiolase (MCKAT), was both necessary and sufficient to elicit the flux decline. Subsequently, Metabolic Control Analysis revealed that MCKAT had insufficient capacity to cope with high substrate influx. Next, we quantified the internal metabolic regulation, revealing a vicious cycle around MCKAT. Upon substrate overload, MCKAT’s ketoacyl-CoA substrates started to accumulate. The unfavourable equilibrium constant of the preceding enzyme, medium/short-chain hydroxyacyl-CoA dehydrogenase, worked as an amplifier, leading to accumulation of upstream CoA esters, including acyl-CoA esters. These acyl-CoA esters are at the same time products of MCKAT and inhibited its already low activity further. Finally, the accumulation of CoA esters led to a sequestration of free CoA. CoA being a cofactor for MCKAT, its sequestration limited the MCKAT activity even further, thus completing the vicious cycle. Since CoA is also a substrate for distant enzymes, it efficiently communicated the ‘traffic jam’ at MCKAT to the entire pathway. This novel mechanism provides a basis to explore the role of mFAO in disease and elucidate similar principles in other pathways of lipid metabolism. [ABSTRACT FROM AUTHOR]

At old age, humans generally have declining muscle mass and increased fat deposition, which can increase the risk of developing cardiometabolic diseases. While regular physical activity postpones these age-related derangements, this is not always possible in the elderly because of disabilities or risk of injury. Whole-body vibration (WBV) training may be considered as an alternative to physical activity particularly in the frail population. To explore this possibility, we characterized whole-body and organ-specific metabolic processes in 6-month and 25-month old mice, over a period of 14 weeks of WBV versus sham training. WBV training tended to increase blood glucose turnover rates and stimulated hepatic glycogen utilization during fasting irrespective of age. WBV was effective in reducing white fat mass and hepatic triglyceride content only in old but not in young mice and these reductions were related to upregulation of hepatic mitochondrial uncoupling of metabolism (assessed by high-resolution respirometry) and increased expression of uncoupling protein 2. Because these changes occurred independent of changes in food intake and whole-body metabolic rate (assessed by indirect calorimetry), the liver-specific effects of WBV may be a primary mechanism to improve metabolic health during aging, rather than that it is a consequence of alterations in energy balance. [ABSTRACT FROM AUTHOR]

The dietary fiber guar gum has beneficial effects on obesity, hyperglycemia and hypercholesterolemia in both humans and rodents. The major products of colonic fermentation of dietary fiber, the short-chain fatty acids (SCFAs), have been suggested to play an important role. Recently, we showed that SCFAs protect against the metabolic syndrome via a signaling cascade that involves peroxisome proliferator-activated receptor (PPAR) γ repression and AMP-activated protein kinase (AMPK) activation. In this study we investigated the molecular mechanism via which the dietary fiber guar gum protects against the metabolic syndrome. C57Bl/6J mice were fed a high-fat diet supplemented with 0% or 10% of the fiber guar gum for 12 weeks and effects on lipid and glucose metabolism were studied. We demonstrate that, like SCFAs, also guar gum protects against high-fat diet-induced metabolic abnormalities by PPARγ repression, subsequently increasing mitochondrial uncoupling protein 2 expression and AMP/ATP ratio, leading to the activation of AMPK and culminating in enhanced oxidative metabolism in both liver and adipose tissue. Moreover, guar gum markedly increased peripheral glucose clearance, possibly mediated by the SCFA-induced colonic hormone glucagon-like peptide-1. Overall, this study provides novel molecular insights into the beneficial effects of guar gum on the metabolic syndrome and strengthens the potential role of guar gum as a dietary-fiber intervention. [ABSTRACT FROM AUTHOR]

Studies with dietary supplementation of various types of fibers have shown beneficial effects on symptoms of the metabolic syndrome. Short-chain fatty acids (SCFAs), the main products of intestinal bacterial fermentation of dietary fiber, have been suggested to play a key role. Whether the concentration of SCFAs or their metabolism drives these beneficial effects is not yet clear. In this study we investigated the SCFA concentrations and in vivo host uptake fluxes in the absence or presence of the dietary fiber guar gum. C57Bl/6J mice were fed a high-fat diet supplemented with 0%, 5%, 7.5% or 10% of the fiber guar gum. To determine the effect on SCFA metabolism, 13C-labeled acetate, propionate or butyrate were infused into the cecum of mice for 6 h and the isotopic enrichment of cecal SCFAs was measured. The in vivo production, uptake and bacterial interconversion of acetate, propionate and butyrate were calculated by combining the data from the three infusion experiments in a single steady-state isotope model. Guar gum treatment decreased markers of the metabolic syndrome (body weight, adipose weight, triglycerides, glucose and insulin levels and HOMA-IR) in a dose-dependent manner. In addition, hepatic mRNA expression of genes involved in gluconeogenesis and fatty acid synthesis decreased dose-dependently by guar gum treatment. Cecal SCFA concentrations were increased compared to the control group, but no differences were observed between the different guar gum doses. Thus, no significant correlation was found between cecal SCFA concentrations and metabolic markers. In contrast, in vivo SCFA uptake fluxes by the host correlated linearly with metabolic markers. We argue that in vivo SCFA fluxes, and not concentrations, govern the protection from the metabolic syndrome by dietary fibers. [ABSTRACT FROM AUTHOR]

Dynamic models of metabolism can be useful in identifying potential drug targets, especially in unicellular organisms. A model of glycolysis in the causative agent of human African trypanosomiasis, Trypanosoma brucei, has already shown the utility of this approach. Here we add the pentose phosphate pathway (PPP) of T. brucei to the glycolytic model. The PPP is localized to both the cytosol and the glycosome and adding it to the glycolytic model without further adjustments leads to a draining of the essential bound-phosphate moiety within the glycosome. This phosphate “leak” must be resolved for the model to be a reasonable representation of parasite physiology. Two main types of theoretical solution to the problem could be identified: (i) including additional enzymatic reactions in the glycosome, or (ii) adding a mechanism to transfer bound phosphates between cytosol and glycosome. One example of the first type of solution would be the presence of a glycosomal ribokinase to regenerate ATP from ribose 5-phosphate and ADP. Experimental characterization of ribokinase in T. brucei showed that very low enzyme levels are sufficient for parasite survival, indicating that other mechanisms are required in controlling the phosphate leak. Examples of the second type would involve the presence of an ATP:ADP exchanger or recently described permeability pores in the glycosomal membrane, although the current absence of identified genes encoding such molecules impedes experimental testing by genetic manipulation. Confronted with this uncertainty, we present a modeling strategy that identifies robust predictions in the context of incomplete system characterization. We illustrate this strategy by exploring the mechanism underlying the essential function of one of the PPP enzymes, and validate it by confirming the model predictions experimentally. [ABSTRACT FROM AUTHOR]

Fatty-acid metabolism plays a key role in acquired and inborn metabolic diseases. To obtain insight into the network dynamics of fatty-acid β-oxidation, we constructed a detailed computational model of the pathway and subjected it to a fat overload condition. The model contains reversible and saturable enzyme-kinetic equations and experimentally determined parameters for rat-liver enzymes. It was validated by adding palmitoyl CoA or palmitoyl carnitine to isolated rat-liver mitochondria: without refitting of measured parameters, the model correctly predicted the β-oxidation flux as well as the time profiles of most acyl-carnitine concentrations. Subsequently, we simulated the condition of obesity by increasing the palmitoyl-CoA concentration. At a high concentration of palmitoyl CoA the β-oxidation became overloaded: the flux dropped and metabolites accumulated. This behavior originated from the competition between acyl CoAs of different chain lengths for a set of acyl-CoA dehydrogenases with overlapping substrate specificity. This effectively induced competitive feedforward inhibition and thereby led to accumulation of CoA-ester intermediates and depletion of free CoA (CoASH). The mitochondrial [NAD+]/[NADH] ratio modulated the sensitivity to substrate overload, revealing a tight interplay between regulation of β-oxidation and mitochondrial respiration. [ABSTRACT FROM AUTHOR]

Background and Methodology: Recently, we reported on a new class of naphthoquinone derivatives showing a promising anti-trypanosomatid profile in cell-based experiments. The lead of this series (B6, 2-phenoxy-1,4-naphthoquinone) showed an ED50 of 80 nM against Trypanosoma brucei rhodesiense, and a selectivity index of 74 with respect to mammalian cells. A multitarget profile for this compound is easily conceivable, because quinones, as natural products, serve plants as potent defense chemicals with an intrinsic multifunctional mechanism of action. To disclose such a multitarget profile of B6, we exploited a chemical proteomics approach. Principal Findings: A functionalized congener of B6 was immobilized on a solid matrix and used to isolate target proteins from Trypanosoma brucei lysates. Mass analysis delivered two enzymes, i.e. glycosomal glycerol kinase and glycosomal glyceraldehyde-3-phosphate dehydrogenase, as potential molecular targets for B6. Both enzymes were recombinantly expressed and purified, and used for chemical validation. Indeed, B6 was able to inhibit both enzymes with IC50 values in the micromolar range. The multifunctional profile was further characterized in experiments using permeabilized Trypanosoma brucei cells and mitochondrial cell fractions. It turned out that B6 was also able to generate oxygen radicals, a mechanism that may additionally contribute to its observed potent trypanocidal activity. Conclusions and Significance: Overall, B6 showed a multitarget mechanism of action, which provides a molecular explanation of its promising anti-trypanosomatid activity. Furthermore, the forward chemical genetics approach here applied may be viable in the molecular characterization of novel multitarget ligands. Author Summary: The multitarget approach can represent a promising strategy for the discovery of innovative drug candidates against neglected tropical diseases. However, multitarget drug discovery can be very demanding, because of the highly time-consuming step related to the fine balancing of the biological activities against selected targets. An innovative workflow for discovering multitarget drugs can be envisioned: i) design and synthesis of natural-like compounds; ii) test them using phenotypic cell-based assays; iii) fishing potential targets by means of chemical proteomics. This workflow might rapidly provide new hit candidates that can be further progressed to the hit-to-lead and lead optimization steps of the drug discovery process. The two latter steps can benefit from information on the molecular target(s), which may be identified by chemical proteomics. Herein, we report on the elucidation of the mode of action of a new series of anti-trypanosomal naphthoquinone compounds, previously tested using cell-based assays, by means of chemical proteomics, classical biochemistry, molecular and system biology. [ABSTRACT FROM AUTHOR]

A decade ago, a team of biochemists including two of us, modeled yeast glycolysis and showed that one of the most studied biochemical pathways could not be quite understood in terms of the kinetic properties of the constituent enzymes as measured in cell extract. Moreover, when the same model was later applied to different experimental steady-state conditions, it often exhibited unrestrained metabolite accumulation. Here we resolve this issue by showing that the results of such ab initio modeling are improved substantially by (i) including appropriate allosteric regulation and (ii) measuring the enzyme kinetic parameters under conditions that resemble the intracellular environment. The following modifications proved crucial: (i) implementation of allosteric regulation of hexokinase and pyruvate kinase, (ii) implementation of Vmax values measured under conditions that resembled the yeast cytosol, and (iii) redetermination of the kinetic parameters of glyceraldehyde-3-phosphate dehydrogenase under physiological conditions. Model predictions and experiments were compared under five different conditions of yeast growth and starvation. When either the original model was used (which lacked important allosteric regulation), or the enzyme parameters were measured under conditions that were, as usual, optimal for high enzyme activity, fructose 1,6-bisphosphate and some other glycolytic intermediates tended to accumulate to unrealistically high concentrations. Combining all adjustments yielded an accurate correspondence between model and experiments for all five steady-state and dynamic conditions. This enhances our understanding of in vivo metabolism in terms of in vitro biochemistry. [ABSTRACT FROM AUTHOR]

Kinetic models of metabolism require detailed knowledge of kinetic parameters. However, due to measurement errors or lack of data this knowledge is often uncertain. The model of glycolysis in the parasitic protozoan Trypanosoma brucei is a particularly well analysed example of a quantitative metabolic model, but so far it has been studied with a fixed set of parameters only. Here we evaluate the effect of parameter uncertainty. In order to define probability distributions for each parameter, information about the experimental sources and confidence intervals for all parameters were collected. We created a wiki-based website dedicated to the detailed documentation of this information: the SilicoTryp wiki (http:// silicotryp.ibls.gla.ac.uk/wiki/Glycolysis). Using information collected in the wiki, we then assigned probability distributions to all parameters of the model. This allowed us to sample sets of alternative models, accurately representing our degree of uncertainty. Some properties of the model, such as the repartition of the glycolytic flux between the glycerol and pyruvate producing branches, are robust to these uncertainties. However, our analysis also allowed us to identify fragilities of the model leading to the accumulation of 3-phosphoglycerate and/or pyruvate. The analysis of the control coefficients revealed the importance of taking into account the uncertainties about the parameters, as the ranking of the reactions can be greatly affected. This work will now form the basis for a comprehensive Bayesian analysis and extension of the model considering alternative topologies. [ABSTRACT FROM AUTHOR]

Often, in living cells different molecular species compete for binding to the same molecular target. Typical examples are the competition of genes for the transcription machinery or the competition of mRNAs for the translation machinery. Here we show that such systems have specific regulatory features and how they can be analysed. We derive a theory for molecular competition in parallel reaction networks. Analytical expressions for the response of network fluxes to changes in the total competitor and common target pools indicate the precise conditions for ultrasensitivity and intuitive rules for competitor strength. The calculations are based on measurable concentrations of the competitor-target complexes. We show that kinetic parameters, which are usually tedious to determine, are not required in the calculations. Given their simplicity, the obtained equations are easily applied to networks of any dimension. The new theory is illustrated for competing sigma factors in bacterial transcription and for a genome-wide network of yeast mRNAs competing for ribosomes. We conclude that molecular competition can drastically influence the network fluxes and lead to negative response coefficients and ultrasensitivity. Competitors that bind a large fraction of the target, like bacterial σ70, tend to influence competing pathways strongly. The less a competitor is saturated by the target, the more sensitive it is to changes in the concentration of the target, as well as to other competitors. As a consequence, most of the mRNAs in yeast turn out to respond ultrasensitively to changes in ribosome concentration. Finally, applying the theory to a genome-wide dataset we observe that high and low response mRNAs exhibit distinct Gene Ontology profiles. [ABSTRACT FROM AUTHOR]

Dynamic models of metabolism can be useful in identifying potential drug targets, especially in unicellular organisms. A model of glycolysis in the causative agent of human African trypanosomiasis, Trypanosoma brucei, has already shown the utility of this approach. Here we add the pentose phosphate pathway (PPP) of T. brucei to the glycolytic model. The PPP is localized to both the cytosol and the glycosome and adding it to the glycolytic model without further adjustments leads to a draining of the essential bound-phosphate moiety within the glycosome. This phosphate “leak” must be resolved for the model to be a reasonable representation of parasite physiology. Two main types of theoretical solution to the problem could be identified: (i) including additional enzymatic reactions in the glycosome, or (ii) adding a mechanism to transfer bound phosphates between cytosol and glycosome. One example of the first type of solution would be the presence of a glycosomal ribokinase to regenerate ATP from ribose 5-phosphate and ADP. Experimental characterization of ribokinase in T. brucei showed that very low enzyme levels are sufficient for parasite survival, indicating that other mechanisms are required in controlling the phosphate leak. Examples of the second type would involve the presence of an ATP:ADP exchanger or recently described permeability pores in the glycosomal membrane, although the current absence of identified genes encoding such molecules impedes experimental testing by genetic manipulation. Confronted with this uncertainty, we present a modeling strategy that identifies robust predictions in the context of incomplete system characterization. We illustrate this strategy by exploring the mechanism underlying the essential function of one of the PPP enzymes, and validate it by confirming the model predictions experimentally. [ABSTRACT FROM AUTHOR]