A computational model of mitochondrial AZT metabolism.
Ontology highlight
ABSTRACT:
SUBMITTER: Bradshaw PC
PROVIDER: S-EPMC1316272 | biostudies-literature | 2005 Dec
REPOSITORIES: biostudies-literature
Publications
A computational model of mitochondrial AZT metabolism.
Bradshaw Patrick C PC Li Jiaxin J Samuels David C DC
The Biochemical journal 20051201 Pt 2
The mechanisms of the mitochondrial toxicity of AZT (azidothymidine; zidovudine) are not clear. The two main contenders are the incorporation of phosphorylated AZT into the mtDNA (mitochondrial DNA) and the competitive inhibition of natural deoxynucleotide metabolism. We have built a computational model of AZT metabolism in mitochondria in order to better understand these toxicity mechanisms. The model includes the transport of non-phosphorylated and phosphorylated forms of AZT into mitochondria ...[more]
Similar Datasets
Project description:Iron is essential for all known life due to its redox properties; however, these same properties can also lead to its toxicity in overload through the production of reactive oxygen species. Robust systemic and cellular control are required to maintain safe levels of iron, and the liver seems to be where this regulation is mainly located. Iron misregulation is implicated in many diseases, and as our understanding of iron metabolism improves, the list of iron-related disorders grows. Recent developments have resulted in greater knowledge of the fate of iron in the body and have led to a detailed map of its metabolism; however, a quantitative understanding at the systems level of how its components interact to produce tight regulation remains elusive. A mechanistic computational model of human liver iron metabolism, which includes the core regulatory components, is presented here. It was constructed based on known mechanisms of regulation and on their kinetic properties, obtained from several publications. The model was then quantitatively validated by comparing its results with previously published physiological data, and it is able to reproduce multiple experimental findings. A time course simulation following an oral dose of iron was compared to a clinical time course study and the simulation was found to recreate the dynamics and time scale of the systems response to iron challenge. A disease state simulation of haemochromatosis was created by altering a single reaction parameter that mimics a human haemochromatosis gene (HFE) mutation. The simulation provides a quantitative understanding of the liver iron overload that arises in this disease. This model supports and supplements understanding of the role of the liver as an iron sensor and provides a framework for further modelling, including simulations to identify valuable drug targets and design of experiments to improve further our knowledge of this system.
| S-EPMC3820522 | biostudies-literature
Project description:Mitochondrial dysfunction causes oxidative stress and cardiomyopathy. Oxidative stress also is a side effect of dideoxynucleoside antiretrovirals (NRTI) and is observed in NRTI-induced cardiomyopathy. We show here that treatment with the NRTI AZT {1-[(2R,4S,5S)-4-azido-5-(hydroxymethyl)oxolan-2-yl]-5-methylpyrimidine-2,4-dione} modulates cardiac gene expression epigenetically through production of mitochondrially derived reactive oxygen species. Transgenic mice with ubiquitous expression of mitochondrially targeted catalase (MCAT) and C57Bl/6 wild-type mice littermates (WT) were administered AZT (0.22 mg/day po, 35 days), and cardiac DNA and mRNA were isolated. In AZT-treated WT, 95 cardiac genes were differentially expressed compared with vehicle-treated WTs. When MCAT mice were treated with AZT, each of those 95 genes reverted toward the expression of vehicle-treated WTs. In AZT-treated WT hearts, Mthfr [5,10-methylenetetrahydrofolate reductase; a critical enzyme in synthesis of methionine cycle intermediates including S-adenosylmethionine (SAM)], was overexpressed. Steady-state abundance of SAM in cardiac extracts from AZT-treated MCAT mice increased 60% above that of vehicle-treated MCAT. No such change occurred in WT. AZT caused hypermethylation (47%) and hypomethylation (53%) of differentially methylated DNA regions in WT cardiac DNA. AZT-treated MCAT heart DNA exhibited greater hypermethylation (91%) and less hypomethylation (9%) compared with vehicle-treated MCAT controls. The gene encoding protein kinase C-? displayed multifocal epigenetic regulation caused by oxidative stress. Results show that mitochondrially derived oxidative stress in the heart hinders cardiac DNA methylation, alters steady-state abundance of SAM, alters cardiac gene expression, and promotes characteristic pathophysiological changes of cardiomyopathy. This mechanism for NRTI toxicity offers insight into long-term side effects from these commonly used antiviral agents.
| S-EPMC4593827 | biostudies-literature
Project description:BACKGROUND: Dyslipidemia is an important risk factor for cardiovascular disease and type II diabetes. Lipoprotein diagnostics, such as LDL cholesterol and HDL cholesterol, help to diagnose these diseases. Lipoprotein profile measurements could improve lipoprotein diagnostics, but interpretational complexity has limited their clinical application to date. We have previously developed a computational model called Particle Profiler to interpret lipoprotein profiles. In the current study we further developed and calibrated Particle Profiler using subjects with specific genetic conditions. We subsequently performed technical validation and worked at an initial indication of clinical usefulness starting from available data on lipoprotein concentrations and metabolic fluxes. Since the model outcomes cannot be measured directly, the only available technical validation was corroboration. For an initial indication of clinical usefulness, pooled lipoprotein metabolic flux data was available from subjects with various types of dyslipidemia. Therefore we investigated how well lipoprotein metabolic ratios derived from Particle Profiler distinguished reported dyslipidemic from normolipidemic subjects. RESULTS: We found that the model could fit a range of normolipidemic and dyslipidemic subjects from fifteen out of sixteen studies equally well, with an average 8.8% ± 5.0% fit error; only one study showed a larger fit error. As initial indication of clinical usefulness, we showed that one diagnostic marker based on VLDL metabolic ratios better distinguished dyslipidemic from normolipidemic subjects than triglycerides, HDL cholesterol, or LDL cholesterol. The VLDL metabolic ratios outperformed each of the classical diagnostics separately; they also added power of distinction when included in a multivariate logistic regression model on top of the classical diagnostics. CONCLUSIONS: In this study we further developed, calibrated, and corroborated the Particle Profiler computational model using pooled lipoprotein metabolic flux data. From pooled lipoprotein metabolic flux data on dyslipidemic patients, we derived VLDL metabolic ratios that better distinguished normolipidemic from dyslipidemic subjects than standard diagnostics, including HDL cholesterol, triglycerides and LDL cholesterol. Since dyslipidemias are closely linked to cardiovascular disease and diabetes type II development, lipoprotein metabolic ratios are candidate risk markers for these diseases. These ratios can in principle be obtained by applying Particle Profiler to a single lipoprotein profile measurement, which makes clinical application feasible.
| S-EPMC3305892 | biostudies-other
Project description:The liver is to date the best example of a sexually dimorphic non-reproductive organ. Over 1,000 genes are differentially expressed between sexes indicating that female and male livers are two metabolically distinct organs. The spectrum of liver diseases is broad and is usually prevalent in one or the other sex, with different contributing genetic and environmental factors. It is thus difficult to predict individual's disease outcomes and treatment options. Systems approaches including mathematical modeling can aid importantly in understanding the multifactorial liver disease etiology leading toward tailored diagnostics, prognostics and therapy. The currently established computational models of hepatic metabolism that have proven to be essential for understanding of non-alcoholic fatty liver disease (NAFLD) and hepatocellular carcinoma (HCC) are limited to the description of gender-independent response or reflect solely the response of the males. Herein we present LiverSex, the first sex-based multi-tissue and multi-level liver metabolic computational model. The model was constructed based on in silico liver model SteatoNet and the object-oriented modeling. The crucial factor in adaptation of liver metabolism to the sex is the inclusion of estrogen and androgen receptor responses to respective hormones and the link to sex-differences in growth hormone release. The model was extensively validated on literature data and experimental data obtained from wild type C57BL/6 mice fed with regular chow and western diet. These experimental results show extensive sex-dependent changes and could not be reproduced in silico with the uniform model SteatoNet. LiverSex represents the first large-scale liver metabolic model, which allows a detailed insight into the sex-dependent complex liver pathologies, and how the genetic and environmental factors interact with the sex in disease appearance and progression. We used the model to identify the most important sex-dependent metabolic pathways, which are involved in accumulation of triglycerides representing initial steps of NAFLD. We identified PGC1A, PPARα, FXR, and LXR as regulatory factors that could become important in sex-dependent personalized treatment of NAFLD.
| S-EPMC5907313 | biostudies-literature
Project description:Azidothymidine (AZT) is a reverse transcriptase (RT) inhibitor that efficiently blocks the replication of spumaretroviruses or foamy viruses (FVs). To more precisely elucidate the mechanism of action of the FV RT enzyme, we generated an AZT-resistant FV in cell culture. Biologically resistant virus was obtained for simian foamy virus from macaque (SFVmac), which was insensitive to AZT concentrations of 1 mM, but not for FVs derived from chimpanzees. Nucleotide sequencing revealed four non-silent mutations in the pol gene. Introduction of these mutations into an infectious molecular clone identified all changes to be required for the fully AZT-resistant phenotype of SFVmac. The alteration of individual sites showed that AZT resistance in SFVmac was likely acquired by consecutive acquisition of pol mutations in a defined order, because some alterations on their own did not result in an efficiently replicating virus, neither in the presence nor in the absence of AZT. The introduction of the mutations into the RT of the closely related prototypic FV (PFV) did not yield an AZT-resistant virus, instead they significantly impaired the viral fitness.
| S-EPMC2276251 | biostudies-literature
Project description:Reprogramming of energy metabolism is a hallmark of cancer that enables the cancer cells to meet the increased energetic requirements due to uncontrolled proliferation. One prominent example is pancreatic ductal adenocarcinoma, an aggressive form of cancer with an overall 5-year survival rate of 5%. The reprogramming mechanism in pancreatic cancer involves deregulated uptake of glucose and glutamine and other opportunistic modes of satisfying energetic demands in a hypoxic and nutrient-poor environment. In the current study, we apply systems biology approaches to enable a better understanding of the dynamics of the distinct metabolic alterations in KRAS-mediated pancreatic cancer, with the goal of impeding early cell proliferation by identifying the optimal metabolic enzymes to target. We have constructed a kinetic model of metabolism represented as a set of ordinary differential equations that describe time evolution of the metabolite concentrations in glycolysis, glutaminolysis, tricarboxylic acid cycle and the pentose phosphate pathway. The model is comprised of 46 metabolites and 53 reactions. The mathematical model is fit to published enzyme knockdown experimental data. We then applied the model to perform in silico enzyme modulations and evaluate the effects on cell proliferation. Our work identifies potential combinations of enzyme knockdown, metabolite inhibition, and extracellular conditions that impede cell proliferation. Excitingly, the model predicts novel targets that can be tested experimentally. Therefore, the model is a tool to predict the effects of inhibiting specific metabolic reactions within pancreatic cancer cells, which is difficult to measure experimentally, as well as test further hypotheses toward targeted therapies.
| S-EPMC5388762 | biostudies-literature
Project description: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.
| S-EPMC1201507 | biostudies-literature
Project description:BackgroundIn eukaryotes, folate metabolism is compartmentalized and occurs in both the cytosol and the mitochondria. The function of this compartmentalization and the great changes that occur in the mitochondrial compartment during embryonic development and in rapidly growing cancer cells are gradually becoming understood, though many aspects remain puzzling and controversial.ApproachWe explore the properties of cytosolic and mitochondrial folate metabolism by experimenting with a mathematical model of hepatic one-carbon metabolism. The model is based on known biochemical properties of mitochondrial and cytosolic enzymes. We use the model to study questions about the relative roles of the cytosolic and mitochondrial folate cycles posed in the experimental literature. We investigate: the control of the direction of the mitochondrial and cytosolic serine hydroxymethyltransferase (SHMT) reactions, the role of the mitochondrial bifunctional enzyme, the role of the glycine cleavage system, the effects of variations in serine and glycine inputs, and the effects of methionine and protein loading.ConclusionThe model reproduces many experimental findings and gives new insights into the underlying properties of mitochondrial folate metabolism. Particularly interesting is the remarkable stability of formate production in the mitochondria in the face of large changes in serine and glycine input. The model shows that in the presence of the bifunctional enzyme (as in embryonic tissues and cancer cells), the mitochondria primarily support cytosolic purine and pyrimidine synthesis via the export of formate, while in adult tissues the mitochondria produce serine for gluconeogenesis.
| S-EPMC1713227 | biostudies-literature
Project description:We carried out a series of experiment demonstrating the role of mitochondria in the cytosolic and mitochondrial Ca(2+) transients and compared the results with those from computer simulation. In rat ventricular myocytes, increasing the rate of stimulation (1~3 Hz) made both the diastolic and systolic [Ca(2+)] bigger in mitochondria as well as in cytosol. As L-type Ca(2+) channel has key influence on the amplitude of Ca(2+)-induced Ca(2+) release, the relation between stimulus frequency and the amplitude of Ca(2+) transients was examined under the low density (1/10 of control) of L-type Ca(2+) channel in model simulation, where the relation was reversed. In experiment, block of Ca(2+) uniporter on mitochondrial inner membrane significantly reduced the amplitude of mitochondrial Ca(2+) transients, while it failed to affect the cytosolic Ca(2+) transients. In computer simulation, the amplitude of cytosolic Ca(2+) transients was not affected by removal of Ca(2+) uniporter. The application of carbonyl cyanide 4-(trifluoromethoxy) phenylhydrazone (FCCP) known as a protonophore on mitochondrial membrane to rat ventricular myocytes gradually increased the diastolic [Ca(2+)] in cytosol and eventually abolished the Ca(2+) transients, which was similarly reproduced in computer simulation. The model study suggests that the relative contribution of L-type Ca(2+) channel to total transsarcolemmal Ca(2+) flux could determine whether the cytosolic Ca(2+) transients become bigger or smaller with higher stimulus frequency. The present study also suggests that cytosolic Ca(2+) affects mitochondrial Ca(2+) in a beat-to-beat manner, however, removal of Ca(2+) influx mechanism into mitochondria does not affect the amplitude of cytosolic Ca(2+) transients.
| S-EPMC3186923 | biostudies-literature
Project description:In non-alcoholic fatty liver disease (NAFLD), lipid build-up and the resulting damage is known to occur more severely in pericentral cells. Due to the complexity of studying individual regions of the sinusoid, the causes of this zone specificity and its implications on treatment are largely ignored. In this study, a computational model of liver glucose and lipid metabolism is presented which treats the sinusoid as the repeating unit of the liver rather than the single hepatocyte. This allows for inclusion of zonated enzyme expression by splitting the sinusoid into periportal to pericentral compartments. By simulating insulin resistance (IR) and high intake diets leading to the development of steatosis in the model, we identify key differences between periportal and pericentral cells accounting for higher susceptibility to pericentral steatosis. Secondly, variation between individuals is seen in both susceptibility to steatosis and in its development across the sinusoid. Around 25% of obese individuals do not show excess liver fat, whilst 16% of lean individuals develop NAFLD. Furthermore, whilst pericentral cells tend to show higher lipid levels, variation is seen in the predominant location of steatosis from pericentral to pan-sinusoidal or azonal. Sensitivity analysis was used to identify the processes which have the largest effect on both total hepatic triglyceride levels and on the sinusoidal location of steatosis. As is seen in vivo, steatosis occurs when simulating IR in the model, predominantly due to increased uptake, along with an increase in de novo lipogenesis. Additionally, concentrations of glucose intermediates including glycerol-3-phosphate increased when simulating IR due to inhibited glycogen synthesis. Several differences between zones contributed to a higher susceptibility to steatosis in pericentral cells in the model simulations. Firstly, the periportal zonation of both glycogen synthase and the oxidative phosphorylation enzymes meant that the build-up of glucose intermediates was less severe in the periportal hepatocyte compartments. Secondly, the periportal zonation of the enzymes mediating β-oxidation and oxidative phosphorylation resulted in excess fats being metabolised more rapidly in the periportal hepatocyte compartments. Finally, the pericentral expression of de novo lipogenesis contributed to pericentral steatosis when additionally simulating the increase in sterol-regulatory element binding protein 1c (SREBP-1c) seen in NAFLD patients in vivo. The hepatic triglyceride concentration was predicted to be most sensitive to inter-individual variation in the activity of enzymes which, either directly or indirectly, determine the rate of free fatty acid (FFA) oxidation. The concentration was most strongly dependent on the rate constants for β-oxidation and oxidative phosphorylation. It also showed moderate sensitivity to the rate constants for processes which alter the allosteric inhibition of β-oxidation by acetyl-CoA. The predominant sinusoidal location of steatosis meanwhile was most sensitive variations in the zonation of proteins mediating FFA uptake or triglyceride release as very low density lipoproteins (VLDL). Neither the total hepatic concentration nor the location of steatosis showed strong sensitivity to variations in the lipogenic rate constants.
| S-EPMC5025084 | biostudies-literature