Centralizing neurofibromatosis experimental tool knowledge with the NF Research Tools Database.

Experimental tools and resources, such as animal models, cell lines, antibodies, genetic reagents and biobanks, are key ingredients in biomedical research. Investigators face multiple challenges when trying to understand the availability, applicability and accessibility of these tools. A major challenge is keeping up with current information about the numerous tools available for a particular research problem. A variety of disease-agnostic projects such as the Mouse Genome Informatics database and the Resource Identification Initiative curate a number of types of research tools. Here, we describe our efforts to build upon these resources to develop a disease-specific research tool resource for the neurofibromatosis (NF) research community. This resource, the NF Research Tools Database, is an open-access database that enables the exploration and discovery of information about NF type 1-relevant animal models, cell lines, antibodies, genetic reagents and biobanks. Users can search and explore tools, obtain detailed information about each tool as well as read and contribute their observations about the performance, reliability and characteristics of tools in the database. NF researchers will be able to use the NF Research Tools Database to promote, discover, share, reuse and characterize research tools, with the goal of advancing NF research. Database URL: https://tools.nf.synapse.org/.

Predicting changes in renal metabolism after compound exposure with a genome-scale metabolic model.

The kidneys are metabolically active organs with importance in several physiological tasks such as the secretion of soluble wastes into the urine and synthesizing glucose and oxidizing fatty acids for energy in fasting (non-fed) conditions. Once damaged, the metabolic capability of the kidneys becomes altered. Here, we define metabolic tasks in a computational modeling framework to capture kidney function in an update to the iRno network reconstruction of rat metabolism using literature-based evidence. To demonstrate the utility of iRno for predicting kidney function, we exposed primary rat renal proximal tubule epithelial cells to four compounds with varying levels of nephrotoxicity (acetaminophen, gentamicin, 2,3,7,8-tetrachlorodibenzodioxin, and trichloroethylene) for six and twenty-four hours, and collected transcriptomics and metabolomics data to measure the metabolic effects of compound exposure. For the transcriptomics data, we observed changes in fatty acid metabolism and amino acid metabolism, as well as changes in existing markers of kidney function such as Clu (clusterin). The iRno metabolic network reconstruction was used to predict alterations in these same pathways after integrating transcriptomics data and was able to distinguish between select compound-specific effects on the proximal tubule epithelial cells. Genome-scale metabolic network reconstructions with coupled omics data can be used to predict changes in metabolism as a step towards identifying novel metabolic biomarkers of kidney function and dysfunction.

Systems-level computational modeling demonstrates fuel selection switching in high capacity running and low capacity running rats.

High capacity and low capacity running rats, HCR and LCR respectively, have been bred to represent two extremes of running endurance and have recently demonstrated disparities in fuel usage during transient aerobic exercise. HCR rats can maintain fatty acid (FA) utilization throughout the course of transient aerobic exercise whereas LCR rats rely predominantly on glucose utilization. We hypothesized that the difference between HCR and LCR fuel utilization could be explained by a difference in mitochondrial density. To test this hypothesis and to investigate mechanisms of fuel selection, we used a constraint-based kinetic analysis of whole-body metabolism to analyze transient exercise data from these rats. Our model analysis used a thermodynamically constrained kinetic framework that accounts for glycolysis, the TCA cycle, and mitochondrial FA transport and oxidation. The model can effectively match the observed relative rates of oxidation of glucose versus FA, as a function of ATP demand. In searching for the minimal differences required to explain metabolic function in HCR versus LCR rats, it was determined that the whole-body metabolic phenotype of LCR, compared to the HCR, could be explained by a ~50% reduction in total mitochondrial activity with an additional 5-fold reduction in mitochondrial FA transport activity. Finally, we postulate that over sustained periods of exercise that LCR can partly overcome the initial deficit in FA catabolic activity by upregulating FA transport and/or oxidation processes.

Estrogen maintains mitochondrial content and function in the right ventricle of rats with pulmonary hypertension.

The typical cause of death in pulmonary hypertension (PH) is right ventricular (RV) failure, with females showing better survival rates than males. Recently, metabolic shift and mitochondrial dysfunction have been demonstrated in RV failure secondary to PH In light of evidence showing that estrogen protects mitochondrial function and biogenesis in noncardiovascular systems, we hypothesized that the mechanism by which estrogen preserves RV function is via protection of mitochondrial content and oxidative capacity in PH We used a well-established model of PH (Sugen+Hypoxia) in ovariectomized female rats with/without estrogen treatment. RV functional measures were derived from pressure-volume relationships measured via RV catheterization in live rats. Citrate synthase activity, a marker of mitochondrial density, was measured in both RV and LV tissues. Respiratory capacity of mitochondria isolated from RV was measured using oxygraphy. We found that RV ventricular-vascular coupling efficiency decreased in the placebo-treated SuHx rats (0.78 ± 0.10 vs. 1.50 ± 0.13 in control, P < 0.05), whereas estrogen restored it. Mitochondrial density decreased in placebo-treated SuHx rats (0.12 ± 0.01 vs. 0.15 ± 0.01 U citrate synthase/mg in control, P < 0.05), and estrogen attenuated the decrease. Mitochondrial quality and oxidative capacity tended to be lower in placebo-treated SuHx rats only. The changes in mitochondrial biogenesis and function paralleled the expression levels of PGC-1α in RV Our results suggest that estrogen protects RV function by preserving mitochondrial content and oxidative capacity. This provides a mechanism by which estrogen provides protection in female PH patients and paves the way to develop estrogen and its targets as a novel RV-specific therapy for PH.

Influence of metabolic dysfunction on cardiac mechanics in decompensated hypertrophy and heart failure.

Alterations in energetic state of the myocardium are associated with decompensated heart failure in humans and in animal models. However, the functional consequences of the observed changes in energetic state on mechanical function are not known. The primary aim of the study was to quantify mechanical/energetic coupling in the heart and to determine if energetic dysfunction can contribute to mechanical failure. A secondary aim was to apply a quantitative systems pharmacology analysis to investigate the effects of drugs that target cross-bridge cycling kinetics in heart failure-associated energetic dysfunction. Herein, a model of metabolite- and calcium-dependent myocardial mechanics was developed from calcium concentration and tension time courses in rat cardiac muscle obtained at different lengths and stimulation frequencies. The muscle dynamics model accounting for the effect of metabolites was integrated into a model of the cardiac ventricles to simulate pressure-volume dynamics in the heart. This cardiac model was integrated into a simple model of the circulation to investigate the effects of metabolic state on whole-body function. Simulations predict that reductions in metabolite pools observed in canine models of heart failure can cause systolic dysfunction, blood volume expansion, venous congestion, and ventricular dilation. Simulations also predict that myosin-activating drugs may partially counteract the effects of energetic state on cross-bridge mechanics in heart failure while increasing myocardial oxygen consumption. Our model analysis demonstrates how metabolic changes observed in heart failure are alone sufficient to cause systolic dysfunction and whole-body heart failure symptoms.

Catalytic Coupling of Oxidative Phosphorylation, ATP Demand, and Reactive Oxygen Species Generation.

Competing models of mitochondrial energy metabolism in the heart are highly disputed. In addition, the mechanisms of reactive oxygen species (ROS) production and scavenging are not well understood. To deepen our understanding of these processes, a computer model was developed to integrate the biophysical processes of oxidative phosphorylation and ROS generation. The model was calibrated with experimental data obtained from isolated rat heart mitochondria subjected to physiological conditions and workloads. Model simulations show that changes in the quinone pool redox state are responsible for the apparent inorganic phosphate activation of complex III. Model simulations predict that complex III is responsible for more ROS production during physiological working conditions relative to complex I. However, this relationship is reversed under pathological conditions. Finally, model analysis reveals how a highly reduced quinone pool caused by elevated levels of succinate is likely responsible for the burst of ROS seen during reperfusion after ischemia.

Feedback Regulation and Time Hierarchy of Oxidative Phosphorylation in Cardiac Mitochondria.

To determine how oxidative ATP synthesis is regulated in the heart, the responses of cardiac mitochondria oxidizing pyruvate to alterations in [ATP], [ADP], and inorganic phosphate ([Pi]) were characterized over a range of steady-state levels of extramitochondrial [ATP], [ADP], and [Pi]. Evolution of the steady states of the measured variables with the flux of respiration shows that: (1) a higher phosphorylation potential is achieved by mitochondria at higher [Pi] for a given flux of respiration; (2) the time hierarchy of oxidative phosphorylation is given by phosphorylation subsystem, electron transport chain, and substrate dehydrogenation subsystems listed in increasing order of their response times; (3) the matrix ATP hydrolysis mass action ratio [ADP] × [Pi]/[ATP] provides feedback to the substrate dehydrogenation flux over the entire range of respiratory flux examined in this study; and finally, (4) contrary to previous models of regulation of oxidative phosphorylation, [Pi] does not modulate the activity of complex III.

Stimulatory effects of calcium on respiration and NAD(P)H synthesis in intact rat heart mitochondria utilizing physiological substrates cannot explain respiratory control in vivo.

Mitochondrial TCA cycle dehydrogenase enzymes have been shown to be stimulated by Ca(2+) under various substrate and ADP incubation conditions in an attempt to determine and understand the role of Ca(2+) in maintaining energy homeostasis in working hearts. In this study, we tested the hypothesis that, at physiological temperature and 1 mM extramitochondrial free magnesium, Ca(2+) can stimulate the overall mitochondrial NAD(P)H generation flux in rat heart mitochondria utilizing pyruvate and malate as substrates at both subsaturating and saturating concentrations. In both cases, we found that, in the physiological regime of mitochondrial oxygen consumption observed in the intact animal and in the physiological range of cytosolic Ca(2+) concentration averaged per beat, Ca(2+) had no observable stimulatory effect. A modest apparent stimulatory effect (22-27%) was observable at supraphysiological maximal ADP-stimulated respiration at 2.5 mM initial phosphate. The stimulatory effects observed over the physiological Ca(2+) range are not sufficient to make a significant contribution to the control of oxidative phosphorylation in the heart in vivo.

Common phenotype of resting mouse extensor digitorum longus and soleus muscles: equal ATPase and glycolytic flux during transient anoxia.

Rates of ATPase and glycolysis are several times faster in actively contracting mouse extensor digitorum longus muscle (EDL) than soleus (SOL), but we find these rates are not distinguishable at rest. We used a transient anoxic perturbation of steady state energy balance to decrease phosphocreatine (PCr) reversibly and to measure the rates of ATPase and of lactate production without muscle activation or contraction. The rate of glycolytic ATP synthesis is less than the ATPase rate, accounting for the continual PCr decrease during anoxia in both muscles. We fitted a mathematical model validated with properties of enzymes and solutes measured in vitro and appropriate for the transient perturbation of these muscles to experimental data to test whether the model accounts for the results. Simulations showed equal rates of ATPase and lactate production in both muscles. ATPase controls glycolytic flux by feedback from its products. Adenylate kinase function is critical because a rise in [AMP] is necessary to activate glycogen phosphorylase. ATPase is the primary source of H+ production. The sum of contributions of the 13 reactions of the glycogenolytic and glycolytic network to total proton load is negligible. The stoichiometry of lactate and H+ production is near unity. These results identify a default state of energy metabolism for resting muscle in which there is no difference in the metabolic phenotype of EDL and SOL. Therefore, additional control mechanisms, involving higher ATPase flux and [Ca2+], must exist to explain the well-known difference in glycolytic rates in fast-twitch and slow-twitch muscles in actively contracting muscle.

Computer modeling of mitochondrial tricarboxylic acid cycle, oxidative phosphorylation, metabolite transport, and electrophysiology.

A computational model of mitochondrial metabolism and electrophysiology is introduced and applied to analysis of data from isolated cardiac mitochondria and data on phosphate metabolites in striated muscle in vivo. This model is constructed based on detailed kinetics and thermodynamically balanced reaction mechanisms and a strict accounting of rapidly equilibrating biochemical species. Since building such a model requires introducing a large number of adjustable kinetic parameters, a correspondingly large amount of independent data from isolated mitochondria respiring on different substrates and subject to a variety of protocols is used to parameterize the model and ensure that it is challenged by a wide range of data corresponding to diverse conditions. The developed model is further validated by both in vitro data on isolated cardiac mitochondria and in vivo experimental measurements on human skeletal muscle. The validated model is used to predict the roles of NAD and ADP in regulating the tricarboxylic acid cycle dehydrogenase fluxes, demonstrating that NAD is the more important regulator. Further model predictions reveal that a decrease of cytosolic pH value results in decreases in mitochondrial membrane potential and a corresponding drop in the ability of the mitochondria to synthesize ATP at the hydrolysis potential required for cellular function.