Redox cycling is an understated mechanism of toxicity associated with a plethora of xenobiotics, responsible for preventing the effective treatment of serious conditions such as malaria and cardiomyopathy. Quinone compounds are notorious redox cyclers, present in drugs such as doxorubicin, which is used to treat a host of human cancers. However, the therapeutic index of doxorubicin is undermined by dose-dependent cardiotoxicity, which may be a function of futile redox cycling. In this study, a doxorubicin-specific in silico quinone redox metabolism model is described. Doxorubicin-GSH adduct formation kinetics are thermodynamically estimated from its reduction potential, while the remainder of the model is parameterised using oxygen consumption rate data, indicative of hydroquinone auto-oxidation. The model is then combined with a comprehensive glutathione metabolism model, facilitating the simulation of quinone redox cycling, and adduct-induced GSH depletion. Simulations suggest that glutathione pools are most sensitive to exposure duration at pharmacologically and supra-pharmacologically relevant doxorubicin concentrations. The model provides an alternative method of investigating and quantifying redox cycling induced oxidative stress, circumventing the experimental difficulties of measuring and tracking radical species. This in silico framework provides a platform from which GSH depletion can be explored as a function of a compound’s physicochemical properties.

Cancer cells depend on glucose metabolism via glycolysis as a primary energy source, despite the presence of oxygen and fully functioning mitochondria, in order to promote growth, proliferation and longevity. Glycolysis relies upon NAD+ to accept electrons in the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) reaction, linking the redox state of the cytosolic NAD+ pool to glycolytic rate. The free cytosolic NAD+/NADH ratio is involved in over 700 oxidoreductive enzymatic reactions and as such, the NAD+/NADH ratio is regarded as a metabolic readout of overall cellular redox state. Many experimental techniques that monitor or measure total NAD+ and NADH are unable to distinguish between protein-bound and unbound forms. Yet total NAD+/NADH measurements yield little information, since it is the free forms of NAD+ and NADH that determine the kinetic and thermodynamic influence of redox potential on glycolytic rate. Indirect estimations of free NAD+/NADH are based on the lactate/pyruvate (L/P) ratio at chemical equilibrium, but these measurements are often undermined by high lability. To elucidate the sensitivity of the free NAD+/NADH ratio to changes in extracellular substrate, an in silico model of hepatocarcinoma glycolysis was constructed and validated against in vitro data. Model simulations reveal that over experimentally relevant concentrations, changes in extracellular glucose and lactate concentration during routine cancer cell culture can lead to significant deviations in the NAD+/NADH ratio. Based on the principles of chemical equilibrium, the model provides a platform from which experimentally challenging situations may be examined, suggesting that extracellular substrates play an important role in cellular redox and bioenergetic homeostasis.