Prediction of Oscillations in Glycolysis in Ethanol-Consuming Erythrocyte-Bioreactors

A mathematical model of energy metabolism in erythrocyte-bioreactors loaded with alcohol dehydrogenase and acetaldehyde dehydrogenase was constructed and analyzed. Such erythrocytes can convert ethanol to acetate using intracellular NAD and can therefore be used to treat alcohol intoxication. Analysis of the model revealed that the rate of ethanol consumption by the erythrocyte-bioreactors increases proportionally to the activity of incorporated ethanol-consuming enzymes until their activity reaches a specific threshold level. When the ethanol-consuming enzyme activity exceeds this threshold, the steady state in the model becomes unstable and the model switches to an oscillation mode caused by the competition between glyceraldehyde phosphate dehydrogenase and ethanol-consuming enzymes for NAD. The amplitude and period of metabolite oscillations first increase with the increase in the activity of the encapsulated enzymes. A further increase in these activities leads to a loss of the glycolysis steady state, and a permanent accumulation of glycolytic intermediates. The oscillation mode and the loss of the steady state can lead to the osmotic destruction of erythrocyte-bioreactors due to an accumulation of intracellular metabolites. Our results demonstrate that the interaction of enzymes encapsulated in erythrocyte-bioreactors with erythrocyte metabolism should be taken into account in order to achieve the optimal efficacy of these bioreactors.