Cofactors revisited – Predicting the impact of flavoprotein-related diseases on a genome scale

Flavin adenine dinucleotide (FAD) and its precursor flavin mononucleotide (FMN) are redox cofactors that are required for the activity of more than hundred human enzymes. Mutations in the genes encoding these proteins cause severe phenotypes, including a lack of energy supply and accumulation of toxic intermediates. Ideally, patients should be diagnosed before they show symptoms so that treatment and/or preventive care can start immediately. This can be achieved by standardized newborn screening tests. However, many of the flavin-related diseases lack appropriate biomarker profiles. Genome-scale metabolic models can aid in biomarker research by predicting altered profiles of potential biomarkers. Unfortunately, current models, including the most recent human metabolic reconstructions Recon and HMR, typically treat enzyme-bound flavins incorrectly as free metabolites. This in turn leads to artificial degrees of freedom in pathways that are strictly coupled. Here, we present a reconstruction of human metabolism with a curated and extended flavoproteome. To illustrate the functional consequences, we show that simulations with the curated model – unlike simulations with earlier Recon versions - correctly predict the metabolic impact of multiple-acyl-CoA-dehydrogenase deficiency as well as of systemic flavin-depletion. Moreover, simulations with the new model allowed us to identify a larger number of biomarkers in flavoproteome-related diseases, without loss of accuracy. We conclude that adequate inclusion of cofactors in constraint-based modelling contributes to higher precision in computational predictions.