Ensemble kinetic modelling links residual enzyme activity to clinical symptoms in mitochondrial β-oxidation defects
Ensemble kinetic modelling links residual enzyme activity to clinical symptoms in mitochondrial β-oxidation defects
Odendaal, C.; Krebs, O.; Bakker, B. M.
AbstractThe mitochondrial fatty acid {beta}-oxidation (mFAO) is an important source of energy when carbohydrate stores are depleted. It is also involved in many diseases, including inherited fatty-acid oxidation deficiencies (mFAODs). Patients with the same genetic variant often present with clinically heterogeneous phenotypes, but the mechanisms contributing to this heterogeneity are poorly understood. To investigate the underlying pathophysiology of different mFAODs, we constructed a computational model of mFAO in human liver, based on experimentally determined enzyme kinetics. A recognised, but seldom addressed challenge in metabolic modelling is the substantial uncertainty about kinetic parameter values. Whereas experimental values of some mFAO parameters are quite reproducible, others vary by up to four orders of magnitude between different reports. To address this, we generated an ensemble of kinetic models, each with the same reaction stoichiometry and rate equations, but different kinetic parameters, sampled from distributions of literature-derived values. We also comprehensively report these values and the arguments based on which they were evaluated. The resulting models were validated against available flux data, yielding a final ensemble of 51 valid models. These models recapitulate recent findings about the accumulation of medium-chain acyl-CoAs and the concomitant depletion of free CoA (CoASH) in medium-chain acyl-CoA dehydrogenase deficiency. We applied the ensemble to a set of known mFAODs, separating them into long-chain (LC-) and short-/medium-chain (S/MC-)mFAODs. The residual activity at which clinical symptoms are known to occur corresponded well with the residual activity in the model at which pathway flux was significantly decreased in LC-mFAODs. Residual activity in S/MC-mFAODs correlated less strongly with pathway flux, but these deficiencies did show a combined flux- and CoASH-reduction effect. This comparison is of importance to researchers and clinicians, as it identifies possible ways in which insights about one mFAOD may be applied to another based on shared biochemical properties.