Computer simulation of metabolism in palmitate-perfused rat heart. I. Palmitate oxidation

Personalised modelling of clinical heterogeneity between medium-chain acyl-CoA dehydrogenase patients

BACKGROUND

Monogenetic inborn errors of metabolism cause a wide phenotypic heterogeneity that may even differ between family members carrying the same genetic variant. Computational modelling of metabolic networks may identify putative sources of this inter-patient heterogeneity. Here, we mainly focus on medium-chain acyl-CoA dehydrogenase deficiency (MCADD), the most common inborn error of the mitochondrial fatty acid oxidation (mFAO). It is an enigma why some MCADD patients-if untreated-are at risk to develop severe metabolic decompensations, whereas others remain asymptomatic throughout life. We hypothesised that an ability to maintain an increased free mitochondrial CoA (CoASH) and pathway flux might distinguish asymptomatic from symptomatic patients.

RESULTS

We built and experimentally validated, for the first time, a kinetic model of the human liver mFAO. Metabolites were partitioned according to their water solubility between the bulk aqueous matrix and the inner membrane. Enzymes are also either membrane-bound or in the matrix. This metabolite partitioning is a novel model attribute and improved predictions. MCADD substantially reduced pathway flux and CoASH, the latter due to the sequestration of CoA as medium-chain acyl-CoA esters. Analysis of urine from MCADD patients obtained during a metabolic decompensation showed an accumulation of medium- and short-chain acylcarnitines, just like the acyl-CoA pool in the MCADD model. The model suggested some rescues that increased flux and CoASH, notably increasing short-chain acyl-CoA dehydrogenase (SCAD) levels. Proteome analysis of MCADD patient-derived fibroblasts indeed revealed elevated levels of SCAD in a patient with a clinically asymptomatic state. This is a rescue for MCADD that has not been explored before. Personalised models based on these proteomics data confirmed an increased pathway flux and CoASH in the model of an asymptomatic patient compared to those of symptomatic MCADD patients.

CONCLUSIONS

We present a detailed, validated kinetic model of mFAO in human liver, with solubility-dependent metabolite partitioning. Personalised modelling of individual patients provides a novel explanation for phenotypic heterogeneity among MCADD patients. Further development of personalised metabolic models is a promising direction to improve individualised risk assessment, management and monitoring for inborn errors of metabolism.

Identifying Site-Specific Superoxide and Hydrogen Peroxide Production Rates From the Mitochondrial Electron Transport System Using a Computational Strategy

Mitochondrial reactive oxygen species (ROS) play important roles in cellular signaling; however, certain pathological conditions such as ischemia/reperfusion (I/R) injury disrupt ROS homeostasis and contribute to cell death. A major impediment to developing therapeutic measures against oxidative stress-induced cellular damage is the lack of a quantitative framework to identify the specific sources and regulatory mechanisms of mitochondrial ROS production. We developed a thermodynamically consistent, mass-and-charge balanced, kinetic model of mitochondrial ROS homeostasis focused on redox sites of electron transport chain complexes I, II, and III. The model was calibrated and corroborated using comprehensive data sets relevant to ROS homeostasis. The model predicts that complex I ROS production dominates other sources under conditions favoring a high membrane potential with elevated nicotinamide adenine dinucleotide (NADH) and ubiquinol (QH₂) levels. In general, complex I contributes to significant levels of ROS production under pathological conditions, while complexes II and III are responsible for basal levels of ROS production, especially when QH₂ levels are elevated. The model also reveals that hydrogen peroxide production by complex I underlies the non-linear relationship between ROS emission and O₂ at low O₂ concentrations. Lastly, the model highlights the need to quantify scavenging system activity under different conditions to establish a complete picture of mitochondrial ROS homeostasis. In summary, we describe the individual contributions of the electron transport system complex redox sites to total ROS emission in mitochondria respiring under various combinations of NADH- and Q-linked respiratory fuels under varying workloads.

Bistability in fatty-acid oxidation resulting from substrate inhibition

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.

Dynamic modelling of an ACADS genotype in fatty acid oxidation - Application of cellular models for the analysis of common genetic variants

Background

Genome-wide association studies of common diseases or metabolite quantitative traits often identify common variants of small effect size, which may contribute to phenotypes by modulation of gene expression. Thus, there is growing demand for cellular models enabling to assess the impact of gene regulatory variants with moderate effects on gene expression. Mitochondrial fatty acid oxidation is an important energy metabolism pathway. Common noncoding acyl-CoA dehydrogenase short chain (ACADS) gene variants are associated with plasma C4-acylcarnitine levels and allele-specific modulation of ACADS expression may contribute to the observed phenotype.

Methods and findings

We assessed ACADS expression and intracellular acylcarnitine levels in human lymphoblastoid cell lines (LCL) genotyped for a common ACADS variant associated with plasma C4-acylcarnitine and found a significant genotype-dependent decrease of ACADS mRNA and protein. Next, we modelled gradual decrease of ACADS expression using a tetracycline-regulated shRNA-knockdown of ACADS in Huh7 hepatocytes, a cell line with high fatty acid oxidation-(FAO)-capacity. Assessing acylcarnitine flux in both models, we found increased C4-acylcarnitine levels with decreased ACADS expression levels. Moreover, assessing time-dependent changes of acylcarnitine levels in shRNA-hepatocytes with altered ACADS expression levels revealed an unexpected effect on long- and medium-chain fatty acid intermediates.

Conclusions

Both, genotyped LCL and regulated shRNA-knockdown are valuable tools to model moderate, gradual gene-regulatory effects of common variants on cellular phenotypes. Decreasing ACADS expression levels modulate short and surprisingly also long/medium chain acylcarnitines, and may contribute to increased plasma acylcarnitine levels.

Measurement of β-oxidation capacity of biological samples by respirometry: a review of principles and substrates

Oxidation of fatty acids is a major source of energy in the heart, liver, and skeletal muscle. It can be measured accurately using respirometry in isolated mitochondria, intact cells, and permeabilized cells or tissues. This technique directly measures the rate of oxygen consumption or flux at various respiratory states when appropriate substrates, uncouplers, and inhibitors are used. Acylcarnitines such as palmitoylcarnitine or octanoylcarnitine are the commonly used substrates. The β-oxidation pathway is prone to feedforward inhibition resulting from accumulation of short-chain acyl-CoA and depletion of CoA, but inclusion of malate or carnitine prevents accumulation of these intermediaries and CoA depletion.

Biochemical competition makes fatty-acid β-oxidation vulnerable to substrate overload

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.

Lipid metabolism plays an important role in the development of metabolic syndrome, a major risk factor for cardiovascular disease and diabetes. Furthermore, inborn errors in lipid oxidation cause rare, but severe diseases in children. To obtain more insight into the response of lipid oxidation to dietary and medical interventions, we constructed a computational model. The model correctly simulated the rate of lipid oxidation and the time courses of most acyl carnitines. The latter are used as diagnostic markers in blood. Subsequently, we subjected the model to an increased supply of lipids, as often happens in obese people. We discovered that the lipid-oxidation machinery easily becomes overloaded, very much like a highway during rush hours: the more cars enter the road, the slower they proceed and the more they clog the road. Analogously, an overload of lipids slowed down the lipid oxidation and led to an accumulation of intermediate metabolites in the pathway. Potential protection mechanisms of cells consist of restricted entry of lipids into the oxidation pathway or efficient downstream processing of reaction products. In future research we will use the model to test dietary or medical interventions in silico and thereby guide the development of new treatment and prevention strategies.

Modeling to link regional myocardial work, metabolism and blood flows

Given the mono-functional, highly coordinated processes of cardiac excitation and contraction, the observations that regional myocardial blood flows, rMBF, are broadly heterogeneous has provoked much attention, but a clear explanation has not emerged. In isolated and in vivo heart studies the total coronary flow is found to be proportional to the rate-pressure product (systolic mean blood pressure times heart rate), a measure of external cardiac work. The same relationship might be expected on a local basis: more work requires more flow. The validity of this expectation has never been demonstrated experimentally. In this article we review the concepts linking cellular excitation and contractile work to cellular energetics and ATP demand, substrate utilization, oxygen demand, vasoregulation, and local blood flow. Mathematical models of these processes are now rather well developed. We propose that the construction of an integrated model encompassing the biophysics, biochemistry and physiology of cardiomyocyte contraction, then combined with a detailed three-dimensional structuring of the fiber bundle and sheet arrangements of the heart as a whole will frame an hypothesis that can be quantitatively evaluated to settle the prime issue: Does local work drive local flow in a predictable fashion that explains the heterogeneity? While in one sense one can feel content that work drives flow is irrefutable, the are no cardiac contractile models that demonstrate the required heterogeneity in local strain-stress-work; quite the contrary, cardiac contraction models have tended toward trying to show that work should be uniform. The object of this review is to argue that uniformity of work does not occur, and is impossible in any case, and that further experimentation and analysis are necessary to test the hypothesis.

Integrative modeling of the cardiac ventricular myocyte

Cardiac electrophysiology is a discipline with a rich 50-year history of experimental research coupled with integrative modeling which has enabled us to achieve a quantitative understanding of the relationships between molecular function and the integrated behavior of the cardiac myocyte in health and disease. In this paper, we review the development of integrative computational models of the cardiac myocyte. We begin with a historical overview of key cardiac cell models that helped shape the field. We then narrow our focus to models of the cardiac ventricular myocyte and describe these models in the context of their subcellular functional systems including dynamic models of voltage-gated ion channels, mitochondrial energy production, ATP-dependent and electrogenic membrane transporters, intracellular Ca dynamics, mechanical contraction, and regulatory signal transduction pathways. We describe key advances and limitations of the models as well as point to new directions for future modeling research. WIREs Syst Biol Med 2011 3 392-413 DOI: 10.1002/wsbm.122

Computer modeling of mitochondrial tricarboxylic acid cycle, oxidative phosphorylation, metabolite transport, and electrophysiology

A computational model of mitochondrial metabolism and electrophysiology is introduced and applied to analysis of data from isolated cardiac mitochondria and data on phosphate metabolites in striated muscle in vivo. This model is constructed based on detailed kinetics and thermodynamically balanced reaction mechanisms and a strict accounting of rapidly equilibrating biochemical species. Since building such a model requires introducing a large number of adjustable kinetic parameters, a correspondingly large amount of independent data from isolated mitochondria respiring on different substrates and subject to a variety of protocols is used to parameterize the model and ensure that it is challenged by a wide range of data corresponding to diverse conditions. The developed model is further validated by both in vitro data on isolated cardiac mitochondria and in vivo experimental measurements on human skeletal muscle. The validated model is used to predict the roles of NAD and ADP in regulating the tricarboxylic acid cycle dehydrogenase fluxes, demonstrating that NAD is the more important regulator. Further model predictions reveal that a decrease of cytosolic pH value results in decreases in mitochondrial membrane potential and a corresponding drop in the ability of the mitochondria to synthesize ATP at the hydrolysis potential required for cellular function.

Modeling the isolated cardiac myocyte

Computer modeling of cardiac myocytes has flourished in recent years. Models have evolved from mathematical descriptions of ionic channels alone to more sophisticated formulations that include calcium transport mechanisms, ATP production and metabolic pathways. The increased complexity is fueled by the new data available in the field. The continuous production of experimental data has led to the evolution of increasingly refined descriptions of the phenomena by modelers. Integrating the numerous systems involved in cardiac myocyte homeostasis makes the use of computer models necessary due to the unreliability of intuitive approaches. However the complexity of the model should not imply a cumbersome operation of the program. As with any tool, computer models have to be easy to operate or their strength will be diminished and potential users will not benefit fully from them. The contribution of the computer modeler to their respective biological fields will be more successful and enduring if modelers devote sufficient time to implement their equations into a model with user-friendly characteristics.

The computational integrated myocyte: a view into the virtual heart

The presentation outlines an integrative approach for developing a computational model of cardiomyocytes. A modular approach is proposed, and strategies of linking the modules (intermediary metabolism, electrophysiology, and mechanics) of the model are presented. A strong recommendation is given toward an integrated system approach backed by experimental validation.

Regulation of myocardial triacylglycerol synthesis and metabolism

Studies showing a correlation of excess myocardial triacylglycerol stores with apoptosis, fibrosis, and contractile dysfunction indicate that dysregulation of triacylglycerol metabolism may contribute to cardiac disease. This review covers the regulation of heart triacylglycerol accumulation at the critical control points of fatty acid uptake, enzymes of triacylglycerol synthesis, lipolysis, and lipoprotein secretion. These pathways are discussed in the context of the central role myocardial triacylglycerol plays in cardiac energy metabolism and heart disease.

An integrated model of cardiac mitochondrial energy metabolism and calcium dynamics

We present an integrated thermokinetic model describing control of cardiac mitochondrial bioenergetics. The model describes the tricarboxylic acid (TCA) cycle, oxidative phosphorylation, and mitochondrial Ca(2+) handling. The kinetic component of the model includes effectors of the TCA cycle enzymes regulating production of NADH and FADH(2), which in turn are used by the electron transport chain to establish a proton motive force (Delta mu(H)), driving the F(1)F(0)-ATPase. In addition, mitochondrial matrix Ca(2+), determined by Ca(2+) uniporter and Na(+)/Ca(2+) exchanger activities, regulates activity of the TCA cycle enzymes isocitrate dehydrogenase and alpha-ketoglutarate dehydrogenase. The model is described by twelve ordinary differential equations for the time rate of change of mitochondrial membrane potential (Delta Psi(m)), and matrix concentrations of Ca(2+), NADH, ADP, and TCA cycle intermediates. The model is used to predict the response of mitochondria to changes in substrate delivery, metabolic inhibition, the rate of adenine nucleotide exchange, and Ca(2+). The model is able to reproduce, qualitatively and semiquantitatively, experimental data concerning mitochondrial bioenergetics, Ca(2+) dynamics, and respiratory control. Significant increases in oxygen consumption (V(O(2))), proton efflux, NADH, and ATP synthesis, in response to an increase in cytoplasmic Ca(2+), are obtained when the Ca(2+)-sensitive dehydrogenases are the main rate-controlling steps of respiratory flux. These responses diminished when control is shifted downstream (e.g., the respiratory chain or adenine nucleotide translocator). The time-dependent behavior of the model, under conditions simulating an increase in workload, closely reproduces experimentally observed mitochondrial NADH dynamics in heart trabeculae subjected to changes in pacing frequency. The steady-state and time-dependent behavior of the model support the hypothesis that mitochondrial matrix Ca(2+) plays an important role in matching energy supply with demand in cardiac myocytes.

Control of mitochondrial beta-oxidation flux

The control of mitochondrial beta-oxidation, including the delivery of acyl moieties from the plasma membrane to the mitochondrion, is reviewed. Control of beta-oxidation flux appears to be largely at the level of entry of acyl groups to mitochondria, but is also dependent on substrate supply. CPTI has much of the control of hepatic beta-oxidation flux, and probably exerts high control in intact muscle because of the high concentration of malonyl-CoA in vivo. beta-Oxidation flux can also be controlled by the redox state of NAD/NADH and ETF/ETFH(2). Control by [acetyl-CoA]/[CoASH] may also be significant, but it is probably via export of acyl groups by carnitine acylcarnitine translocase and CPT II rather than via accumulation of 3-ketoacyl-CoA esters. The sharing of control between CPTI and other enzymes allows for flexible regulation of metabolism and the ability to rapidly adapt beta-oxidation flux to differing requirements in different tissues.

Subcellular metabolic transients and mitochondrial redox waves in heart cells

Precise matching of energy supply with demand requires delicately balanced control of the enzymes involved in substrate metabolism. In response to a change in substrate supply, the nonlinear properties of metabolic control may induce complex dynamic behavior. Using confocal imaging of flavoprotein redox potential and mitochondrial membrane potential, we show that substrate deprivation leads to subcellular heterogeneity of mitochondrial energization in intact cells. The complex spatiotemporal patterns of redox and matrix potential included local metabolic transients, cell-wide coordinated redox transitions, and propagated metabolic waves both within and between coupled cells. Loss of metabolic synchrony during mild metabolic stress reveals that intra- and intercellular control of mitochondrial function involves diffusible cytoplasmic messengers.

Toward modeling the human physionome

The physionome is the description of the physiological dynamics of the normal intact organism. The march of science brings us now into the era where integration of the various facets of the knowledge of biology and medicine has become a major issue. Modeling is a vehicle for the combining of information from molecular biology, biophysics, and medical biology, but must be combined with strategies for databasing the raw data with greater efficiency than is currently possible. The lessons from the genome project can be applied to the next level major projects, the morphonome and the physionome, the objective being to put integrated forms of the data into the hands of physicians and medical scientists.

A computer model of pancreatic islet glycolysis

We have modeled an experiment with perifused pancreatic islet cells using our BIOSSIM language. The experiment and the resulting model are concerned with glucose uptake and glycolysis by the beta-cells of pancreatic islets. Although glycolysis appears to be involved in insulin release, we do not have enough information to represent insulin release in detail. The rapid entry of glucose into the beta-cell is promoted by a carrier having a very high tissue capacity. Phosphorylation of glucose by the low affinity enzyme glucokinase appears to be limiting for glycolysis. The effects of several hexose diphosphate activators of phosphofructokinase are modeled. Model behavior is described. The kinetic parameters of the enzyme submodels are given. Because of the difficulties of preparing large amounts of experimental material, information on pancreatic islet metabolism is limited. This model is a plausible explanation of the experimental results. Recent work on the genetically engineered glucose transporter and glucokinase is discussed.

Identification of regulatory properties of metabolic networks by graph theoretical modeling

An earlier graph theoretical model of metabolic and gene-expression networks has been modified and extended to include the effect of electrical potentials on binding constants, representation of uncatalyzed processes, and treatment of parallel reactions catalyzed by a single enzyme. Formal operations on the graph, which are facilitated by a set of standardized guidelines, identify the feedback signals in the network and rank them according to their influence. The technique was applied to a model of glycolysis in ascites tumor cells in the absence and presence of 12.5 mM exogenous glucose. Feedback regulation was widely distributed and mostly due to binding of adenine nucleotide cofactors to the enzymes of the network. The major changes in feedback regulation on adding glucose is the relief of inhibition of hexokinase and phosphofructokinase and the activation of pyruvate kinase. We conclude that regulation of tumor cell glycolysis is not restricted to hexokinase or to (Na+,K+)-ATPase as was previously suggested by others.

Computer simulation of metabolism in palmitate-perfused rat heart. III. Sensitivity analysis

The behavior of a computer model of metabolism in glucose- and palmitate-perfused rat hearts was interpreted by sensitivity analysis to explain why the heart preferentially utilizes fatty acids as fuel even in the presence of substantial exogenous glucose. The sensitivity functions identified those metabolites and enzymes which were most important in regulating the metabolic rate and determined which enzymes set the levels of the critical metabolites. Control of the mitochondrial redox potential and the distribution of coenzyme A thioesters regulated the rate of fatty acid utilization while strong inhibition of citrate synthetase resulted in accumulation of acetyl CoA and suppression of pyruvate oxidation. Glycolysis was limited by the cytosolic ATP/ADP ratio set largely by the creatine shuttle. Metabolic control appears to be widely distributed rather than localized at "key" enzymes. Metabolite levels are usually set by enzymes controlled by modifiers whereas metabolic flux is regulated by the enzymes that produce ligands for the modifier-controlled enzymes.

Computer simulation of metabolism in palmitate-perfused rat heart. II. Behavior of complete model

Intermediary metabolism in rat hearts perfused with 11 mM glucose plus 1 mM palmitate was simulated by a computer model. Several enzyme submodels in a previous version of the isolated rat heart computer model were improved, and a new fatty acid oxidation pathway model was added. Compartmentation of metabolites in a pseudo-stationary state was calculated, and its implications are discussed, e.g., citrate level may not regulate glycolysis because it is mostly mitochondrial. Citrate synthetase, controlled largely by its inhibitors, is of key importance in regulating fatty acid metabolism. The response of aconitase to the mitochondrial Mg2+ level is of major importance in setting both the mitochondrial citrate and isocitrate levels. Pyruvate dehydrogenase is about 96% in the inactive phosphorylated form, and the active form is also 15% inhibited by products, severely limiting pyruvate oxidation and causing preferential utilization of palmitate as the metabolic fuel. The simulation is consistent with a creatine phosphate shuttle which delivers high energy phosphate to the site of its utilization for mechanical work.

Computer simulation of metabolism in palmitate-perfused rat heart. I. Palmitate oxidation

A computer model of the fatty acid oxidation pathway in perfused rat heart was constructed. It includes uptake, activation, and beta-oxidation of fatty acids, triglyceride synthesis and hydrolysis, and carnitine-dependent transport of acyl groups across the mitochondrial membrane under pseudosteady state conditions. Fatty acid utilization may be limited by beta-oxidation in hypoxia or ischemia but probably not in aerobic conditions. Nonesterified fatty acids bound to proteins are found to be metabolically available. The model predicts that stearate, but not palmitate, can support the highest observed respiration rate for perfused rat heart without supplementation by other substrates. Fatty acids are preferentially oxidized rather than being stored as triglycerides because the cystosolic acyl CoA level is lower than the Km for triglyceride synthesis. It is suggested that feedback inhibition of triglyceride lipase regulates utilization of triglycerides as fuel in aerobic hearts.