Search

Abdominal aortic aneurysms (AAAs) are a degenerative aortic disease and associated with hallmarks of aging, such as mitophagy. Despite this, the exact associations among mitophagy, aging, and AAA progression remain unknown. In our study, gene expression analysis of human AAA tissue revealed downregulation of mitophagy pathways, mitochondrial structure, and function-related proteins. Human proteomic analyses identified decreased levels of mitophagy mediators PINK1 and Parkin. Aged mice and, separately, a murine AAA model showed reduced mitophagy in aortic vascular smooth muscle cells (VSMCs) and PINK1 and Parkin expression. Parkin knockdown in VSMCs aggravated AAA dilation in murine models, with elevated mitochondrial ROS and impaired mitochondrial function. Importantly, inhibiting USP30, an antagonist of the PINK1/Parkin pathway, increased mitophagy in VSMCs, improved mitochondrial function, and reduced AAA incidence and growth. Our study elucidates a critical mechanism that proposes AAAs as an age-associated disease with altered mitophagy, introducing new potential therapeutic approaches.

2'-deoxy-ATP (dATP) improves cardiac function by increasing the rate of crossbridge cycling and Ca[Formula: see text] transient decay. However, the mechanisms of these effects and how therapeutic responses to dATP are achieved when dATP is only a small fraction of the total ATP pool remain poorly understood. Here, we used a multiscale computational modeling approach to analyze the mechanisms by which dATP improves ventricular function. We integrated atomistic simulations of prepowerstroke myosin and actomyosin association, filament-scale Markov state modeling of sarcomere mechanics, cell-scale analysis of myocyte Ca[Formula: see text] dynamics and contraction, organ-scale modeling of biventricular mechanoenergetics, and systems level modeling of circulatory dynamics. Molecular and Brownian dynamics simulations showed that dATP increases the actomyosin association rate by 1.9 fold. Markov state models predicted that dATP increases the pool of myosin heads available for crossbridge cycling, increasing steady-state force development at low dATP fractions by 1.3 fold due to mechanosensing and nearest-neighbor cooperativity. This was found to be the dominant mechanism by which small amounts of dATP can improve contractile function at myofilament to organ scales. Together with faster myocyte Ca[Formula: see text] handling, this led to improved ventricular contractility, especially in a failing heart model in which dATP increased ejection fraction by 16% and the energy efficiency of cardiac contraction by 1%. This work represents a complete multiscale model analysis of a small molecule myosin modulator from single molecule to organ system biophysics and elucidates how the molecular mechanisms of dATP may improve cardiovascular function in heart failure with reduced ejection fraction., Competing Interests: Competing interests statement:A.D.M. is a cofounder of and has equity interests in Insilicomed Inc. and Vektor Medical, Inc. He serves as scientific advisor to both companies. Some of A.D.M.’s research grants have been identified for conflict of interest management based on the overall scope of the project and its potential benefit to Insilicomed Inc. and Vektor Medical, Inc. The author is required to disclose this relationship in publications acknowledging the grant support; however, the research subject and findings reported in this study did not involve the companies in any way and have no relationship with the business activities or scientific interests of either company. The terms of this arrangement have been reviewed and approved by the University of California San Diego in accordance with its conflict of interest policies.

Computational, or in silico, models are an effective, noninvasive tool for investigating cardiovascular function. These models can be used in the analysis of experimental and clinical data to identify possible mechanisms of (ab)normal cardiovascular physiology. Recent advances in computing power and data management have led to innovative and complex modeling frameworks that simulate cardiovascular function across multiple scales. While commonly used in multiple disciplines, there is a lack of concise guidelines for the implementation of computer models in cardiovascular research. In line with recent calls for more reproducible research, it is imperative that scientists adhere to credible practices when developing and applying computational models to their research. The goal of this manuscript is to provide a consensus document that identifies best practices for in silico computational modeling in cardiovascular research. These guidelines provide the necessary methods for mechanistic model development, model analysis, and formal model calibration using fundamentals from statistics. We outline rigorous practices for computational, mechanistic modeling in cardiovascular research and discuss its synergistic value to experimental and clinical data.

While exercise-mediated vasoregulation in the myocardium is understood to be governed by autonomic, myogenic, and metabolic-mediated mechanisms, we do not yet understand the spatial heterogeneity of vasodilation or its effects on microvascular flow patterns and oxygen delivery. This study uses a simulation and modeling approach to explore the mechanisms underlying the recruitment of myocardial perfusion and oxygen delivery in exercise. The simulation approach integrates model components representing: whole-body cardiovascular hemodynamics, cardiac mechanics and myocardial work; myocardial perfusion; and myocardial oxygen transport. Integrating these systems together, model simulations reveal: (1.) To match expected flow and transmural flow ratios at increasing levels of exercise, a greater degree of vasodilation must occur in the subendocardium compared to the subepicardium. (2.) Oxygen extraction and venous oxygenation are predicted to substantially decrease with increasing exercise level preferentially in the subendocardium, suggesting that an oxygen-dependent error signal driving metabolic mediated recruitment of flow would be operative only in the subendocardium. (3.) Under baseline physiological conditions approximately 4% of the oxygen delivered to the subendocardium may be supplied via retrograde flow from coronary veins., Competing Interests: Declaration of competing interest None., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

The adult mammalian heart has limited regenerative capacity following injury, leading to progressive heart failure and mortality. Recent studies have identified the spiny mouse ( Acomys ) as a unique model for mammalian cardiac isch3emic resilience, exhibiting enhanced recovery after myocardial infarction (MI) compared to commonly used laboratory mouse strains. However, the underlying cellular and molecular mechanisms behind this unique response remain poorly understood. In this study, we comprehensively characterized the metabolic characteristics of cardiomyocytes in Acomys compared to the non-regenerative Mus musculus . We utilized single-nucleus RNA sequencing (snRNA-seq) in sham-operated animals and 1, 3, and 7 days post-myocardial infarction to investigate cardiomyocytes' transcriptomic and metabolomic profiles in response to myocardial infarction. Complementary targeted metabolomics, stable isotope-resolved metabolomics, and functional mitochondrial assays were performed on heart tissues from both species to validate the transcriptomic findings and elucidate the metabolic adaptations in cardiomyocytes following ischemic injury. Transcriptomic analysis revealed that Acomys cardiomyocytes inherently upregulate genes associated with glycolysis, the pentose phosphate pathway, and glutathione metabolism while downregulating genes involved in oxidative phosphorylation (OXPHOS). These metabolic characteristics are linked to decreased reactive oxygen species (ROS) production and increased antioxidant capacity. Our targeted metabolomic studies in heart tissue corroborated these findings, showing a shift from fatty acid oxidation to glycolysis and ancillary biosynthetic pathways in Acomys at baseline with adaptive changes post-MI. Functional mitochondrial studies indicated a higher reliance on glycolysis in Acomys compared to Mus , underscoring the unique metabolic phenotype of Acomys hearts. Stable isotope tracing experiments confirmed a shift in glucose utilization from oxidative phosphorylation in Acomys . In conclusion, our study identifies unique metabolic characteristics of Acomys cardiomyocytes that contribute to their enhanced ischemic resilience following myocardial infarction. These findings provide novel insights into the role of metabolism in regulating cardiac repair in adult mammals. Our work highlights the importance of inherent and adaptive metabolic flexibility in determining cardiomyocyte ischemic responses and establishes Acomys as a valuable model for studying cardiac ischemic resilience in adult mammals.

The coronary circulation has the inherent ability to maintain myocardial perfusion constant over a wide range of perfusion pressures. The phenomenon of pressure-flow autoregulation is crucial in response to flow-limiting atherosclerotic lesions which diminish coronary driving pressure and increase risk of myocardial ischemia and infarction. Despite well over half a century of devoted research, understanding of the mechanisms responsible for autoregulation remains one of the most fundamental and contested questions in the field today. The purpose of this review is to highlight current knowledge regarding the complex interrelationship between the pathways and mechanisms proposed to dictate the degree of coronary pressure-flow autoregulation. Our group recently likened the intertwined nature of the essential determinants of coronary flow control to the symbolically unsolvable "Gordian knot". To further efforts to unravel the autoregulatory "knot", we consider recent challenges to the local metabolic and myogenic hypotheses and the complicated dynamic structural and functional heterogeneity unique to the heart and coronary circulation. Additional consideration is given to interrogation of putative mediators, role of K + and Ca 2+ channels, and recent insights from computational modeling studies. Improved understanding of how specific vasoactive mediators, pathways, and underlying disease states influence coronary pressure-flow relations stands to significantly reduce morbidity and mortality for what remains the leading cause of death worldwide., Competing Interests: Declaration of competing interest The authors have no conflicts of interest to disclose., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Purine nucleotides play central roles in energy metabolism in the heart. Most fundamentally, the free energy of hydrolysis of the adenine nucleotide adenosine triphosphate (ATP) provides the thermodynamic driving force for numerous cellular processes including the actin-myosin crossbridge cycle. Perturbations to ATP supply and/or demand in the myocardium lead to changes in the homeostatic balance between purine nucleotide synthesis, degradation, and salvage, potentially affecting myocardial energetics and, consequently, myocardial mechanics. Indeed, both acute myocardial ischemia and decompensatory remodeling of the myocardium in heart failure are associated with depletion of myocardial adenine nucleotides and with impaired myocardial mechanical function. Yet there remain gaps in the understanding of mechanistic links between adenine nucleotide degradation and contractile dysfunction in heart disease. The scope of this article is to: (i) review current knowledge of the pathways of purine nucleotide depletion and salvage in acute ischemia and in chronic heart disease; (ii) review hypothesized mechanisms linking myocardial mechanics and energetics with myocardial adenine nucleotide regulation; and (iii) highlight potential targets for treating myocardial metabolic and mechanical dysfunction associated with these pathways. It is hypothesized that an imbalance in the degradation, salvage, and synthesis of adenine nucleotides leads to a net loss of adenine nucleotides in both acute ischemia and under chronic high-demand conditions associated with the development of heart failure. This reduction in adenine nucleotide levels results in reduced myocardial ATP and increased myocardial inorganic phosphate. Both of these changes have the potential to directly impact tension development and mechanical work at the cellular level. © 2024 American Physiological Society. Compr Physiol 14:5345-5369, 2024., (Copyright © 2024 American Physiological Society. All rights reserved.)

Objective: Cardiorespiratory fitness (CRF) is tightly linked with health and longevity and is implicated in metabolic flexibility and substrate metabolism. The high capacity runner (HCR) and low capacity runner (LCR) rat lines are a genetically heterogeneous rat model selected and bred for CRF that reflect CRF in humans by exhibiting differences in nutrient handling. This study aims to differentiate the intrinsic substrate preference of the HCR compared to LCR rats to better understand the intersection of mitochondrial respiration and intrinsic CRF., Methods: We performed bulk skeletal muscle RNA-Sequencing on male and female HCR and LCR rats and assessed the effect of rat line on mitochondrial gene expression pathways using the MitoCarta3.0 database. In a separate cohort of rats, mitochondria were isolated from skeletal and cardiac muscle and maximal oxidation rates were measured using an Oroboros O2k when provided either pyruvate or fatty acid substrates., Results: The expression of mitochondrial genes are significantly upregulated in HCR skeletal muscle in both male and female rats. In respirometry experiments, fatty acid oxidative capacities were greater in HCR compared to LCR, and male compared to female rats, as a function of both mitochondrial quality and mitochondrial density. This effect was greater in the skeletal muscle than in the heart. Pyruvate oxidation did not differ significantly between lines., Conclusions: The capacity for increased fatty acid oxidation in the HCR rat is a result of selection for running capacity and is likely a key contributor to the healthy metabolic phenotype of individuals with high CRF., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier GmbH.. All rights reserved.)

Introduction: The left (LV) and right (RV) ventricles are linked biologically, hemodynamically, and mechanically, a phenomenon known as ventricular interdependence. While LV function has long been known to impact RV function, the reverse is increasingly being realized to have clinical importance. Investigating ventricular interdependence clinically is challenging given the invasive measurements required, including biventricular catheterization, and confounding factors such as comorbidities, volume status, and other aspects of subject variability. Methods: Computational modeling allows investigation of mechanical and hemodynamic interactions in the absence of these confounding factors. Here, we use a threesegment biventricular heart model and simple circulatory system to investigate ventricular interdependence under conditions of systolic and diastolic dysfunction of the LV and RV in the presence of compensatory volume loading. We use the end-diastolic pressure-volume relationship, end-systolic pressure-volume relationship, Frank Starling curves, and cardiac power output as metrics. Results: The results demonstrate that LV systolic and diastolic dysfunction lead to RV compensation as indicated by increases in RV power. Additionally, RV systolic and diastolic dysfunction lead to impaired LV filling, interpretable as LV stiffening especially with volume loading to maintain systemic pressure. Discussion: These results suggest that a subset of patients with intact LV systolic function and diagnosed to have impaired LV diastolic function, categorized as heart failure with preserved ejection fraction (HFpEF), may in fact have primary RV failure. Application of this computational approach to clinical data sets, especially for HFpEF, may lead to improved diagnosis and treatment strategies and consequently improved outcomes., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Kim, Randall, Jezek, Beard and Chesler.)

Introduction: As liver disease progresses, scarring results in worsening hemodynamics ultimately culminating in portal hypertension. This process has classically been quantified through the portosystemic pressure gradient (PSG), which is clinically estimated by hepatic venous pressure gradient (HVPG); however, PSG alone does not predict a given patient's clinical trajectory regarding the Baveno stage of cirrhosis. We hypothesize that a patient's PSG sensitivity to venous remodeling could explain disparate disease trajectories., Methods: We created a computational model of the portal system in the context of worsening liver disease informed by physiologic measurements from the field of portal hypertension. We simulated progression of clinical complications, HVPG, and transjugular intrahepatic portosystemic shunt placement while only varying a patient's likelihood of portal venous remodeling., Results: Our results unify hemodynamics, venous remodeling, and the clinical progression of liver disease into a mathematically consistent model of portal hypertension. We find that by varying how sensitive patients are to create venous collaterals with rising PSG we can explain variation in patterns of decompensation for patients with liver disease. Specifically, we find that patients who have higher proportions of portosystemic shunting earlier in disease have an attenuated rise in HVPG, delayed onset of ascites, and less hemodynamic shifting after transjugular intrahepatic portosystemic shunt placement., Discussion: This article builds a computational model of portal hypertension which supports that patient-level differences in venous remodeling may explain disparate clinical trajectories of disease., (Copyright © 2023 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of The American College of Gastroenterology.)

The coronary circulation has an innate ability to maintain constant blood flow over a wide range of perfusion pressures. However, the mechanisms responsible for coronary autoregulation remain a fundamental and highly contested question. This study interrogated the local metabolic hypothesis of autoregulation by testing the hypothesis that hypoxemia-induced exaggeration of the metabolic error signal improves the autoregulatory response. Experiments were performed on open-chest anesthetized swine during stepwise changes in coronary perfusion pressure (CPP) from 140 to 40 mmHg under normoxic (n = 15) and hypoxemic (n = 8) conditions, in the absence and presence of dobutamine-induced increases in myocardial oxygen consumption (MVO 2 ) (n = 5-7). Hypoxemia (PaO 2  < 40 mmHg) decreased coronary venous PO 2 (CvPO 2 ) ~ 30% (P < 0.001) and increased coronary blood flow ~ 100% (P < 0.001), sufficient to maintain myocardial oxygen delivery (P = 0.14) over a wide range of CPPs. Autoregulatory responsiveness during hypoxemia-induced reductions in CvPO 2 were associated with increases of autoregulatory gain (Gc; P = 0.033) but not slope (P = 0.585) over a CPP range of 120 to 60 mmHg. Preservation of autoregulatory Gc (P = 0.069) and slope (P = 0.264) was observed during dobutamine administration ± hypoxemia. Reductions in coronary resistance in response to decreases in CPP predominantly occurred below CvPO 2 values of ~ 25 mmHg, irrespective of underlying vasomotor reserve. These findings support the presence of an autoregulatory threshold under which oxygen-sensing pathway(s) act to preserve sufficient myocardial oxygen delivery as CPP is reduced during increases in MVO 2 and/or reductions in arterial oxygen content., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany.)

Flavin containing monooxygenases (FMOs) are promiscuous enzymes known for metabolizing a wide range of exogenous compounds. In C. elegans, fmo-2 expression increases lifespan and healthspan downstream of multiple longevity-promoting pathways through an unknown mechanism. Here, we report that, beyond its classification as a xenobiotic enzyme, fmo-2 expression leads to rewiring of endogenous metabolism principally through changes in one carbon metabolism (OCM). These changes are likely relevant, as we find that genetically modifying OCM enzyme expression leads to alterations in longevity that interact with fmo-2 expression. Using computer modeling, we identify decreased methylation as the major OCM flux modified by FMO-2 that is sufficient to recapitulate its longevity benefits. We further find that tryptophan is decreased in multiple mammalian FMO overexpression models and is a validated substrate for FMO-2. Our resulting model connects a single enzyme to two previously unconnected key metabolic pathways and provides a framework for the metabolic interconnectivity of longevity-promoting pathways such as dietary restriction. FMOs are well-conserved enzymes that are also induced by lifespan-extending interventions in mice, supporting a conserved and important role in promoting health and longevity through metabolic remodeling., (© 2023. The Author(s).)

The pancreatic tumor microenvironment drives deregulated nutrient availability. Accordingly, pancreatic cancer cells require metabolic adaptations to survive and proliferate. Pancreatic cancer subtypes have been characterized by transcriptional and functional differences, with subtypes reported to exist within the same tumor. However, it remains unclear if this diversity extends to metabolic programming. Here, using metabolomic profiling and functional interrogation of metabolic dependencies, we identify two distinct metabolic subclasses among neoplastic populations within individual human and mouse tumors. Furthermore, these populations are poised for metabolic cross-talk, and in examining this, we find an unexpected role for asparagine supporting proliferation during limited respiration. Constitutive GCN2 activation permits ATF4 signaling in one subtype, driving excess asparagine production. Asparagine release provides resistance during impaired respiration, enabling symbiosis. Functionally, availability of exogenous asparagine during limited respiration indirectly supports maintenance of aspartate pools, a rate-limiting biosynthetic precursor. Conversely, depletion of extracellular asparagine with PEG-asparaginase sensitizes tumors to mitochondrial targeting with phenformin., (© 2022. The Author(s).)

For patients with heart failure, myocardial ATP level can be reduced to one-half of that observed in healthy controls. This marked reduction (from ≈8 mM in healthy controls to as low as 3-4 mM in heart failure) has been suggested to contribute to impaired myocardial contraction and to the decreased pump function characteristic of heart failure. However, in vitro measures of maximum myofilament force generation, maximum shortening velocity, and the actomyosin ATPase activity show effective K M values for MgATP ranging from ≈10 μM to 150 μM, well below the intracellular ATP level in heart failure. Thus, it is not clear that the fall of myocardial ATP observed in heart failure is sufficient to impair the function of the contractile proteins. Therefore, we tested the effect of low MgATP levels on myocardial contraction using demembranated cardiac muscle preparations that were exposed to MgATP levels typical of the range found in non-failing and failing hearts. Consistent with previous studies, we found that a 50% reduction in MgATP level (from 8 mM to 4 mM) did not reduce maximum force generation or maximum velocity of shortening. However, we found that a 50% reduction in MgATP level caused a 20%-25% reduction in maximal power generation (measured during muscle shortening against a load) and a 20% slowing of cross-bridge cycling kinetics. These results suggest that the decreased cellular ATP level occurring in heart failure contributes to the impaired pump function of the failing heart. Since the ATP-myosin ATPase dissociation constant is estimated to be submillimolar, these findings also suggest that MgATP concentration affects cross-bridge dynamics through a mechanism that is more complex than through the direct dependence of MgATP concentration on myosin ATPase activity. Finally, these studies suggest that therapies targeted to increase adenine nucleotide pool levels in cardiomyocytes might be beneficial for treating heart failure., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

Coronary blood flow is tightly regulated to ensure that myocardial oxygen delivery meets local metabolic demand via the concurrent action of myogenic, neural and metabolic mechanisms. Although several competing hypotheses exist, the specific nature of the local metabolic mechanism(s) remains poorly defined. To gain insights into the viability of putative metabolic feedback mechanisms and into the co-ordinated action of parallel regulatory mechanisms, we applied a multiscale modelling framework to analyse experimental data on coronary pressure, flow and myocardial oxygen delivery in the porcine heart in vivo. The modelling framework integrates a previously established lumped-parameter model of myocardial perfusion used to account for transmural haemodynamic variations and a simple vessel mechanics model used to simulate the vascular tone in each of three myocardial layers. Vascular tone in the resistance vessel mechanics model is governed by input stimuli from the myogenic, metabolic and autonomic control mechanisms. Seven competing formulations of the metabolic feedback mechanism are implemented in the modelling framework, and associated model simulations are compared with experimental data on coronary pressures and flows under a range of experimental conditions designed to interrogate the governing control mechanisms. Analysis identifies a maximally probable metabolic mechanism among the seven tested models, in which production of a metabolic signalling factor is proportional to myocardial oxygen consumption and delivery is proportional to flow. Finally, the identified model is validated based on comparisons of simulations with data on the myocardial perfusion response to conscious exercise that were not used for model identification. KEY POINTS: Although several competing hypotheses exist, we lack knowledge of specific nature of the metabolic mechanism(s) governing regional myocardial perfusion. Moreover, we lack an understanding of how parallel myogenic, adrenergic/autonomic and metabolic mechanisms work together to regulatory oxygen delivery in the beating heart. We have developed a multiscale modelling framework to test competing hypotheses against experimental data on coronary pressure, flow and myocardial oxygen delivery in the porcine heart in vivo. The analysis identifies a maximally probable metabolic mechanism among seven tested models, in which the production of a metabolic signalling factor is proportional to myocardial oxygen consumption and delivery is proportional to flow., (© 2022 The Authors. The Journal of Physiology © 2022 The Physiological Society.)

Blood flows and pressures throughout the human cardiovascular system are regulated in response to various dynamic perturbations, such as changes to peripheral demands in exercise, rapid changes in posture, or loss of blood from hemorrhage, via the coordinated action of the heart, the vasculature, and autonomic reflexes. To assess how the systemic and pulmonary arterial and venous circulation, the heart, and the baroreflex work together to effect the whole-body responses to these perturbations, we integrated an anatomically-based large-vessel arterial tree model with the TriSeg heart model, models capturing nonlinear characteristics of the large and small veins, and baroreflex-mediated regulation of vascular tone and cardiac chronotropy and inotropy. The model was identified by matching data from the Valsalva maneuver (VM), exercise, and head-up tilt (HUT). Thirty-one parameters were optimized using a custom parameter-fitting tool chain, resulting in an unique, high-fidelity whole-body human cardiovascular systems model. Because the model captures the effects of exercise and posture changes, it can be used to simulate numerous clinical assessments, such as HUT, the VM, and cardiopulmonary exercise stress testing. The model can also be applied as a framework for representing and simulating individual patients and pathologies. Moreover, it can serve as a framework for integrating multi-scale organ-level models, such as for the heart or the kidneys, into a whole-body model. Here, the model is used to analyze the relative importance of chronotropic, inotropic, and peripheral vascular contributions to the whole-body cardiovascular response to exercise. It is predicted that in normal physiological conditions chronotropy and inotropy make roughly equal contributions to increasing cardiac output and cardiac power output during exercise. Under upright exercise conditions, the nonlinear pressure-volume relationship of the large veins and sympathetic-mediated venous vasoconstriction are both required to maintain preload to achieve physiological exercise levels. The developed modeling framework is built using the open Modelica modeling language and is freely distributed., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

To phenotype mechanistic differences between heart failure with reduced (HFrEF) and preserved (HFpEF) ejection fraction, a closed-loop model of the cardiovascular system coupled with patient-specific transthoracic echocardiography (TTE) and right heart catheterization (RHC) data was used to identify key parameters representing haemodynamics. Thirty-one patient records (10 HFrEF, 21 HFpEF) were obtained from the Cardiovascular Health Improvement Project database at the University of Michigan. Model simulations were tuned to match RHC and TTE pressure, volume, and cardiac output measurements in each patient. The underlying physiological model parameters were plotted against model-based norms and compared between HFrEF and HFpEF. Our results confirm the main mechanistic parameter driving HFrEF is reduced left ventricular (LV) contractility, whereas HFpEF exhibits a heterogeneous phenotype. Conducting principal component analysis, k -means clustering, and hierarchical clustering on the optimized parameters reveal (i) a group of HFrEF-like HFpEF patients (HFpEF1), (ii) a classic HFpEF group (HFpEF2), and (iii) a group of HFpEF patients that do not consistently cluster (NCC). These subgroups cannot be distinguished from the clinical data alone. Increased LV active contractility ( p < 0.001 ) and LV passive stiffness ( p < 0.001 ) at rest are observed when comparing HFpEF2 to HFpEF1. Analysing the clinical data of each subgroup reveals that elevated systolic and diastolic LV volumes seen in both HFrEF and HFpEF1 may be used as a biomarker to identify HFrEF-like HFpEF patients. These results suggest that modelling of the cardiovascular system and optimizing to standard clinical data can designate subgroups of HFpEF as separate phenotypes, possibly elucidating patient-specific treatment strategies. KEY POINTS: Analysis of data from right heart catheterization (RHC) and transthoracic echocardiography (TTE) of heart failure (HF) patients using a closed-loop model of the cardiovascular system identifies key parameters representing haemodynamic cardiovascular function in patients with heart failure with reduced and preserved ejection fraction (HFrEF and HFpEF). Analysing optimized parameters representing cardiovascular function using machine learning shows mechanistic differences between HFpEF groups that are not seen analysing clinical data alone. HFpEF groups presented here can be subdivided into three subgroups: HFpEF1 described as 'HFrEF-like HFpEF', HFpEF2 as 'classic HFpEF', and a third group of HFpEF patients that do not consistently cluster. Focusing purely on cardiac function consistently captures the underlying dysfunction in HFrEF, whereas HFpEF is better characterized by dysfunction in the entire cardiovascular system. Our methodology reveals that elevated left ventricular systolic and diastolic volumes are potential biomarkers for identifying HFrEF-like HFpEF patients., (© 2021 The Authors. The Journal of Physiology © 2021 The Physiological Society.)

We present a computational framework for analyzing and simulating mitochondrial ATP synthesis using basic thermodynamic and kinetic principles. The framework invokes detailed descriptions of the thermodynamic driving forces associated with the processes of the electron transport chain, mitochondrial ATP synthetase, and phosphate and adenine nucleotide transporters. Assembling models of these discrete processes into an integrated model of mitochondrial ATP synthesis, we illustrate how to analyze and simulate in vitro respirometry experiments and how models identified from in vitro experimental data effectively explain cardiac respiratory control in vivo. Computer codes for these analyses are embedded as Python scripts in a Jupyter Book to facilitate easy adoption and modification of the concepts developed here. This accessible framework may also prove useful in supporting educational applications. All source codes are available on at https://beards-lab.github.io/QAMAS&#95;book/., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Background: Proton magnetic resonance spectroscopy ( 1 H-MRS) of the human heart is deemed to be a quantitative method to investigate myocardial metabolite content, but thorough validations of in vivo measurements against invasive techniques are lacking., Purpose: To determine measurement precision and accuracy for quantifications of myocardial total creatine and triglyceride content with localized 1 H-MRS., Study Type: Test-retest repeatability and measurement validation study., Subjects: Sixteen volunteers and 22 patients scheduled for open-heart aortic valve replacement or septal myectomy., Field Strength/sequence: Prospectively ECG-triggered respiratory-gated free-breathing single-voxel point-resolved spectroscopy (PRESS) sequence at 3 T., Assessment: Myocardial total creatine and triglyceride content were quantified relative to the total water content by fitting the 1 H-MR spectra. Precision was assessed with measurement repeatability. Accuracy was assessed by validating in vivo 1 H-MRS measurements against biochemical assays in myocardial tissue from the same subjects., Statistical Tests: Intrasession and intersession repeatability was assessed using Bland-Altman analyses. Agreement between 1 H-MRS measurements and biochemical assay was tested with regression analyses., Results: The intersession repeatability coefficient for myocardial total creatine content was 41.8% with a mean value of 0.083% ± 0.020% of the total water signal, and 36.7% for myocardial triglyceride content with a mean value of 0.35% ± 0.13% of the total water signal. Ex vivo myocardial total creatine concentrations in tissue samples correlated with the in vivo myocardial total creatine content measured with 1 H-MRS: n = 22, r = 0.44; P < 0.05. Likewise, ex vivo myocardial triglyceride concentrations correlated with the in vivo myocardial triglyceride content: n = 20, r = 0.50; P < 0.05., Data Conclusion: We validated the use of localized 1 H-MRS of the human heart at 3 T for quantitative assessments of in vivo myocardial tissue metabolite content by estimating the measurement precision and accuracy., Level of Evidence: 2 TECHNICAL EFFICACY STAGE: 2., (© 2021 The Authors. Journal of Magnetic Resonance Imaging published by Wiley Periodicals LLC. on behalf of International Society for Magnetic Resonance in Medicine.)

The relationship between regional variabilities in airflow (ventilation) and blood flow (perfusion) is a critical determinant of gas exchange efficiency in the lungs. Hypoxic pulmonary vasoconstriction is understood to be the primary active regulator of ventilation-perfusion matching, where upstream arterioles constrict to direct blood flow away from areas that have low oxygen supply. However, it is not understood how the integrated action of hypoxic pulmonary vasoconstriction affects oxygen transport at the system level. In this study we develop, and make functional predictions with a multi-scale multi-physics model of ventilation-perfusion matching governed by the mechanism of hypoxic pulmonary vasoconstriction. Our model consists of (a) morphometrically realistic 2D pulmonary vascular networks to the level of large arterioles and venules; (b) a tileable lumped-parameter model of vascular fluid and wall mechanics that accounts for the influence of alveolar pressure; (c) oxygen transport accounting for oxygen bound to hemoglobin and dissolved in plasma; and (d) a novel empirical model of hypoxic pulmonary vasoconstriction. Our model simulations predict that under the artificial test condition of a uniform ventilation distribution (1) hypoxic pulmonary vasoconstriction matches perfusion to ventilation; (2) hypoxic pulmonary vasoconstriction homogenizes regional alveolar-capillary oxygen flux; and (3) hypoxic pulmonary vasoconstriction increases whole-lobe oxygen uptake by improving ventilation-perfusion matching., Competing Interests: The authors have declared that no competing interests exist.

The spontaneously hypertensive rat (SHR) is a genetic model of primary hypertension with an etiology that includes sympathetic overdrive. To elucidate the neurogenic mechanisms underlying the pathophysiology of this model, we analyzed the dynamic baroreflex response to spontaneous fluctuations in arterial pressure in conscious SHRs, as well as in the Wistar-Kyoto (WKY), the Dahl salt-sensitive, the Dahl salt-resistant, and the Sprague-Dawley rat. Observations revealed the existence of long intermittent periods (lasting up to several minutes) of engagement and disengagement of baroreflex control of heart rate. Analysis of these intermittent periods revealed a predictive relationship between increased mean arterial pressure and progressive baroreflex disengagement that was present in the SHR and WKY strains but absent in others. This relationship yielded the hypothesis that a lower proportion of engagement versus disengagement of the baroreflex in SHR compared with WKY contributes to the hypertension (or increased blood pressure) in SHR compared with WKY. Results of experiments using sinoaortic baroreceptor denervation were consistent with the hypothesis that dysfunction of the baroreflex contributes to the etiology of hypertension in the SHR. Thus, this study provides experimental evidence for the roles of the baroreflex in long-term arterial pressure regulation and in the etiology of primary hypertension in this animal model.

The endoplasmic reticulum (ER) engages mitochondria at specialized ER domains known as mitochondria-associated membranes (MAMs). Here, we used three-dimensional high-resolution imaging to investigate the formation of pleomorphic "megamitochondria" with altered MAMs in brown adipocytes lacking the Sel1L-Hrd1 protein complex of ER-associated protein degradation (ERAD). Mice with ERAD deficiency in brown adipocytes were cold sensitive and exhibited mitochondrial dysfunction. ERAD deficiency affected ER-mitochondria contacts and mitochondrial dynamics, at least in part, by regulating the turnover of the MAM protein, sigma receptor 1 (SigmaR1). Thus, our study provides molecular insights into ER-mitochondrial cross-talk and expands our understanding of the physiological importance of Sel1L-Hrd1 ERAD., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Rationale: Aging is one of the strongest risk factors for atherosclerosis. Yet whether aging increases the risk of atherosclerosis independently of chronic hyperlipidemia is not known. Objective: To determine if vascular aging before the induction of hyperlipidemia enhances atherogenesis. Methods and Results: We analyzed the aortas of young and aged normolipidemic wild type, disease-free mice and found that aging led to elevated IL (interleukin)-6 levels and mitochondrial dysfunction, associated with increased mitophagy and the associated protein Parkin. In aortic tissue culture, we found evidence that with aging mitochondrial dysfunction and IL-6 exist in a positive feedback loop. We triggered acute hyperlipidemia in aged and young mice by inducing liver-specific degradation of the LDL (low-density lipoprotein) receptor combined with a 10-week western diet and found that atherogenesis was enhanced in aged wild-type mice. Hyperlipidemia further reduced mitochondrial function and increased the levels of Parkin in the aortas of aged mice but not young mice. Genetic disruption of autophagy in smooth muscle cells of young mice exposed to hyperlipidemia led to increased aortic Parkin and IL-6 levels, impaired mitochondrial function, and enhanced atherogenesis. Importantly, enhancing mitophagy in aged, hyperlipidemic mice via oral administration of spermidine prevented the increase in aortic IL-6 and Parkin, attenuated mitochondrial dysfunction, and reduced atherogenesis. Conclusions: Before hyperlipidemia, aging elevates IL-6 and impairs mitochondrial function within the aorta, associated with enhanced mitophagy and increased Parkin levels. These age-associated changes prime the vasculature to exacerbate atherogenesis upon acute hyperlipidemia. Our work implies that novel therapeutics aimed at improving vascular mitochondrial bioenergetics or reducing inflammation before hyperlipidemia may reduce age-related atherosclerosis.

Cardiac mechanical function is supported by ATP hydrolysis, which provides the chemical-free energy to drive the molecular processes underlying cardiac pumping. Physiological rates of myocardial ATP consumption require the heart to resynthesize its entire ATP pool several times per minute. In the failing heart, cardiomyocyte metabolic dysfunction leads to a reduction in the capacity for ATP synthesis and associated free energy to drive cellular processes. Yet it remains unclear if and how metabolic/energetic dysfunction that occurs during heart failure affects mechanical function of the heart. We hypothesize that changes in phosphate metabolite concentrations (ATP, ADP, inorganic phosphate) that are associated with decompensation and failure have direct roles in impeding contractile function of the myocardium in heart failure, contributing to the whole-body phenotype. To test this hypothesis, a transverse aortic constriction (TAC) rat model of pressure overload, hypertrophy, and decompensation was used to assess relationships between metrics of whole-organ pump function and myocardial energetic state. A multiscale computational model of cardiac mechanoenergetic coupling was used to identify and quantify the contribution of metabolic dysfunction to observed mechanical dysfunction. Results show an overall reduction in capacity for oxidative ATP synthesis fueled by either fatty acid or carbohydrate substrates as well as a reduction in total levels of adenine nucleotides and creatine in myocardium from TAC animals compared to sham-operated controls. Changes in phosphate metabolite levels in the TAC rats are correlated with impaired mechanical function, consistent with the overall hypothesis. Furthermore, computational analysis of myocardial metabolism and contractile dynamics predicts that increased levels of inorganic phosphate in TAC compared to control animals kinetically impair the myosin ATPase crossbridge cycle in decompensated hypertrophy/heart failure., (© The Author(s) 2020. Published by Oxford University Press on behalf of American Physiological Society.)

Changes in the myocardial energetics associated with aging-reductions in creatine phosphate/ATP ratio, total creatine, and ATP-mirror changes observed in failing hearts compared to healthy controls. Similarly, both aging and heart failure are associated with significant reductions in cardiac performance and maximal left ventricular cardiac power output compared with young healthy individuals. Based on these observations, we hypothesize that reductions in the concentrations of cytoplasmic adenine nucleotide, creatine, and phosphate pools that occur with aging impair the myocardial capacity to synthesize ATP at physiological free energy levels and that the resulting changes to myocardial energetic status impair the mechanical pumping ability of the heart. The purpose of this study is to test these hypotheses using an age-structured population model for myocardial metabolism in the adult female population and to determine the potential impact of reductions in key myocardial metabolite pools in causing metabolic/energetic and cardiac mechanical dysfunction associated with aging. To test these hypotheses, we developed a population model for myocardial energetics to predict myocardial ATP, ADP, creatine phosphate, creatine, and inorganic phosphate concentrations as functions of cardiac work and age in the adult female population. Model predictions support our hypotheses and are consistent with previous experimental observations. The major findings provide a novel, to our knowledge, theoretical and computational framework for further probing complex relationships between the energetics and performance of the heart with aging., (Copyright © 2019 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Objectives: This study sought to determine whether salt-induced ANG II suppression contributes to impaired CBF autoregulation., Methods: Cerebral autoregulation was evaluated with LDF during graded reductions of blood pressure. Autoregulatory responses in rats fed HS (4% NaCl) diet vs LS (0.4% NaCl) diet were analyzed using linear regression analysis, model-free analysis, and a mechanistic theoretical model of blood flow through cerebral arterioles., Results: Autoregulation was intact in LS-fed animals as MAP was reduced via graded hemorrhage to approximately 50 mm Hg. Short-term (3 days) and chronic (4 weeks) HS diet impaired CBF autoregulation, as evidenced by progressive reductions of laser Doppler flux with arterial pressure reduction. Chronic low dose ANG II infusion (5 mg/kg/min, i.v.) restored CBF autoregulation between the pre-hemorrhage MAP and 50 mm Hg in rats fed short-term HS diet. Mechanistic-based model analysis showed a reduced myogenic response and reduced baseline VSM tone with short-term HS diet, which was restored by ANG II infusion., Conclusions: Short-term and chronic HS diet lead to impaired autoregulation in the cerebral circulation, with salt-induced ANG II suppression as a major factor in the initiation of impaired CBF regulation., (© 2018 John Wiley & Sons Ltd.)

During aerobic exercise (>65% of maximum oxygen consumption), the primary source of acetyl-CoA to fuel oxidative ATP synthesis in muscle is the pyruvate dehydrogenase (PDH) reaction. This study investigated how regulation of PDH activity affects muscle energetics by determining whether activation of PDH with dichloroacetate (DCA) alters the dynamics of the phosphate potential of rat gastrocnemius muscle during contraction. Twitch contractions were induced in vivo over a broad range of intensities to sample submaximal and maximal aerobic workloads. Muscle phosphorus metabolites were measured in vivo before and after DCA treatment by phosphorus nuclear magnetic resonance spectroscopy. At rest, DCA increased PDH activation compared with control (90 ± 12% vs. 23 ± 3%, P < 0.05), with parallel decreases in inorganic phosphate (P i ) of 17% (1.4 ± 0.2 vs. 1.7 ± 0.1 mM, P < 0.05) and an increase in the free energy of ATP hydrolysis (ΔG ATP ) (-66.2 ± 0.3 vs. -65.6 ± 0.2 kJ/mol, P < 0.05). During stimulation DCA increased steady-state phosphocreatine (PCr) and the magnitude of ΔG ATP , with concomitant reduction in P i and ADP concentrations. These effects were not due to kinetic alterations in PCr hydrolysis, resynthesis, or glycolytic ATP production and altered the flow-force relationship between mitochondrial ATP synthesis rate and ΔG ATP . DCA had no significant effect at 1.0- to 2.0-Hz stimulation because physiological mechanisms at these high stimulation levels cause maximal activation of PDH. These data support a role of PDH activation in the regulation of the energetic steady state by altering the phosphate potential (ΔG ATP ) at rest and during contraction.

Functional metrics of autonomic control of heart rate, including baroreflex sensitivity, have been shown to be strongly associated with cardiovascular risk. A decrease in baroreflex sensitivity with aging is hypothesized to represent a contributing causal factor in the etiology of primary hypertension. To assess baroreflex function in human subjects, two complementary methods to simulate the response in heart rate elicited by the Valsalva maneuver were developed and applied to data obtained from a cohort of healthy normal volunteers. The first method is based on representing the baroreflex arc as a simple linear filter, transforming changes in arterial pressure to changes in R-R interval. The second method invokes a physiologically based model for arterial mechanics, afferent baroreceptor strain-dependent firing, and control of heart rate via central autonomic response to changes in afferent inputs from aortic and carotid sensors. Analysis based on the linear filter model reveals that the effective response time of the baroreflex arc tends to increase with age in healthy subjects and that the response time/response rate is a predictor of resting systolic pressure. Similar trends were obtained based on the physiologically based model. Analysis of the Valsalva response using the physiologically based model further reveals that different afferent inputs from the carotid sinus and the aortic arch baroreceptors govern different parts of the heart rate response. The observed relationship between baroreflex sensitivity and systolic pressure is surprising because hypertensive subjects were excluded from the study, and there was no observed relationship between arterial pressure and age. NEW & NOTEWORTHY We introduce two methods to assess baroreflex function from data recorded from human subjects performing the Valsalva maneuver. Results demonstrate that the baroreflex response time tends to increase with age in healthy subjects, that response time represents a predictor of resting systolic pressure, and that the Valsalva response reveals different effects mediated by baroreceptors in the carotid sinus compared with those in the aortic arch.

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., Competing Interests: We have the following interests. R.A.K. acknowledges funding support from a BBSRC Industry Case Partnership with GSK (BB/L502273/2), https://bbsrc.ukri.org, https://www.gsk.com. Andy Harrell is employed by GlaxoSmithKline. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials, as detailed online in the guide for authors.

Background Histologic examination of fixed renal tissue is widely used to assess morphology and the progression of disease. Commonly reported metrics include glomerular number and injury. However, characterization of renal histology is a time-consuming and user-dependent process. To accelerate and improve the process, we have developed a glomerular localization pipeline for trichrome-stained kidney sections using a machine learning image classification algorithm. Methods We prepared 4- μ m slices of kidneys from rats of various genetic backgrounds that were subjected to different experimental protocols and mounted the slices on glass slides. All sections used in this analysis were trichrome stained and imaged in bright field at a minimum resolution of 0.92 μ m per pixel. The training and test datasets for the algorithm comprised 74 and 13 whole renal sections, respectively, totaling over 28,000 glomeruli manually localized. Additionally, because this localizer will be ultimately used for automated assessment of glomerular injury, we assessed bias of the localizer for preferentially identifying healthy or damaged glomeruli. Results Localizer performance achieved an average precision and recall of 96.94% and 96.79%, respectively, on whole kidney sections without evidence of bias for or against glomerular injury or the need for manual preprocessing. Conclusions This study presents a novel and robust application of convolutional neural nets for the localization of glomeruli in healthy and damaged trichrome-stained whole-renal section mounts and lays the groundwork for automated glomerular injury scoring., (Copyright © 2018 by the American Society of Nephrology.)

Right ventricular (RV) failure, which occurs in the setting of pressure overload, is characterized by abnormalities in mechanical and energetic function. The effects of these cell- and tissue-level changes on organ-level RV function are unknown. The primary aim of this study was to investigate the effects of myofiber mechanics and mitochondrial energetics on organ-level RV function in the context of pressure overload using a multiscale model of the cardiovascular system. The model integrates the mitochondria-generated metabolite concentrations that drive intracellular actin-myosin cross-bridging and extracellular myocardial tissue mechanics in a biventricular heart model coupled with simple lumped parameter circulations. Three types of pressure overload were simulated and compared to experimental results. The computational model was able to capture a wide range of cardiovascular physiology and pathophysiology from mild RV dysfunction to RV failure. Our results confirm that, in response to pressure overload alone, the RV is able to maintain cardiac output (CO) and predict that alterations in either RV active myofiber mechanics or RV metabolite concentrations are necessary to decrease CO.

Introduction: Pulmonary hypertension (PH) causes pressure overload leading to right ventricular failure (RVF). Myocardial structure and myocyte mechanics are altered in RVF but the direct impact of these cellular level factors on organ level function remain unclear. A computational model of the cardiovascular system that integrates cellular function into whole organ function has recently been developed. This model is a useful tool for investigating how changes in myocyte structure and mechanics contribute to organ function. We use this model to determine how measured changes in myocyte and myocardial mechanics contribute to RVF at the organ level and predict the impact of myocyte-targeted therapy. Methods: A multiscale computational framework was tuned to model PH due to bleomycin exposure in mice. Pressure overload was modeled by increasing the pulmonary vascular resistance (PVR) and decreasing pulmonary artery compliance (CPA). Myocardial fibrosis and the impairment of myocyte maximum force generation (Fmax) were simulated by increasing the collagen content (↑PVR + ↓CPA + fibrosis) and decreasing Fmax (↑PVR + ↓CPA + fibrosis + ↓Fmax). A61603 (A6), a selective α 1A -subtype adrenergic receptor agonist, shown to improve Fmax was simulated to explore targeting myocyte generated Fmax in PH. Results: Increased afterload (RV systolic pressure and arterial elastance) in simulations matched experimental results for bleomycin exposure. Pressure overload alone (↑PVR + ↓CPA) caused decreased RV ejection fraction (EF) similar to experimental findings but preservation of cardiac output (CO). Myocardial fibrosis in the setting of pressure overload (↑PVR + ↓PAC + fibrosis) had minimal impact compared to pressure overload alone. Including impaired myocyte function (↑PVR + ↓PAC + fibrosis + ↓Fmax) reduced CO, similar to experiment, and impaired EF. Simulations predicted that A6 treatment preserves EF and CO despite maintained RV pressure overload. Conclusion: Multiscale computational modeling enabled prediction of the contribution of cellular level changes to whole organ function. Impaired Fmax is a key feature that directly contributes to RVF. Simulations further demonstrate the therapeutic benefit of targeting Fmax, which warrants additional study. Future work should incorporate growth and remodeling into the computational model to enable prediction of the multiscale drivers of the transition from dysfunction to failure.

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.

Phosphorus-31 magnetic resonance spectroscopy ( 31 P-MRS) is a unique non-invasive imaging modality for probing in vivo high-energy phosphate metabolism in the human heart. We investigated whether current 31 P-MRS methodology would allow for clinical applications to detect exercise-induced changes in (patho-)physiological myocardial energy metabolism. Hereto, measurement variability and repeatability of three commonly used localized 31 P-MRS methods [3D image-selected in vivo spectroscopy (ISIS) and 1D ISIS with 1D chemical shift imaging (CSI) oriented either perpendicular or parallel to the surface coil] to quantify the myocardial phosphocreatine (PCr) to adenosine triphosphate (ATP) ratio in healthy humans ( n = 8) at rest were determined on a clinical 3 Tesla MR system. Numerical simulations of myocardial energy homeostasis in response to increased cardiac work rates were performed using a biophysical model of myocardial oxidative metabolism. Hypertrophic cardiomyopathy was modeled by either inefficient sarcomere ATP utilization or decreased mitochondrial ATP synthesis. The effect of creatine depletion on myocardial energy homeostasis was explored for both conditions. The mean in vivo myocardial PCr/ATP ratio measured with 3D ISIS was 1.57 ± 0.17 with a large repeatability coefficient of 40.4%. For 1D CSI in a 1D ISIS-selected slice perpendicular to the surface coil, the PCr/ATP ratio was 2.78 ± 0.50 (repeatability 42.5%). With 1D CSI in a 1D ISIS-selected slice parallel to the surface coil, the PCr/ATP ratio was 1.70 ± 0.56 (repeatability 43.7%). The model predicted a PCr/ATP ratio reduction of only 10% at the maximal cardiac work rate in normal myocardium. Hypertrophic cardiomyopathy led to lower PCr/ATP ratios for high cardiac work rates, which was exacerbated by creatine depletion. Simulations illustrated that when conducting cardiac 31 P-MRS exercise stress testing with large measurement error margins, results obtained under pathophysiologic conditions may still lie well within the 95% confidence interval of normal myocardial PCr/ATP dynamics. Current measurement precision of localized 31 P-MRS for quantification of the myocardial PCr/ATP ratio precludes the detection of the changes predicted by computational modeling. This hampers clinical employment of 31 P-MRS for diagnostic testing and risk stratification, and warrants developments in cardiac 31 P-MRS exercise stress testing methodology.

Multi-scale models can facilitate whole plant simulations by linking gene networks, protein synthesis, metabolic pathways, physiology, and growth. Whole plant models can be further integrated with ecosystem, weather, and climate models to predict how various interactions respond to environmental perturbations. These models have the potential to fill in missing mechanistic details and generate new hypotheses to prioritize directed engineering efforts. Outcomes will potentially accelerate improvement of crop yield, sustainability, and increase future food security. It is time for a paradigm shift in plant modeling, from largely isolated efforts to a connected community that takes advantage of advances in high performance computing and mechanistic understanding of plant processes. Tools for guiding future crop breeding and engineering, understanding the implications of discoveries at the molecular level for whole plant behavior, and improved prediction of plant and ecosystem responses to the environment are urgently needed. The purpose of this perspective is to introduce Crops in silico (cropsinsilico.org), an integrative and multi-scale modeling platform, as one solution that combines isolated modeling efforts toward the generation of virtual crops, which is open and accessible to the entire plant biology community. The major challenges involved both in the development and deployment of a shared, multi-scale modeling platform, which are summarized in this prospectus, were recently identified during the first Crops in silico Symposium and Workshop.

Systems research spanning fields from biology to finance involves the identification of models to represent the underpinnings of complex systems. Formal approaches for data-driven identification of network interactions include statistical inference-based approaches and methods to identify dynamical systems models that are capable of fitting multivariate data. Availability of large data sets and so-called 'big data' applications in biology present great opportunities as well as major challenges for systems identification/reverse engineering applications. For example, both inverse identification and forward simulations of genome-scale gene regulatory network models pose compute-intensive problems. This issue is addressed here by combining the processing power of Graphics Processing Units (GPUs) and a parallel reverse engineering algorithm for inference of regulatory networks. It is shown that, given an appropriate data set, information on genome-scale networks (systems of 1000 or more state variables) can be inferred using a reverse-engineering algorithm in a matter of days on a small-scale modern GPU cluster.

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., (© 2017 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of The Physiological Society and the American Physiological Society.)

Individualized modeling and simulation of blood flow mechanics find applications in both animal research and patient care. Individual animal or patient models for blood vessel mechanics are based on combining measured vascular geometry with a fluid structure model coupling formulations describing dynamics of the fluid and mechanics of the wall. For example, one-dimensional fluid flow modeling requires a constitutive law relating vessel cross-sectional deformation to pressure in the lumen. To investigate means of identifying appropriate constitutive relationships, an automated segmentation algorithm was applied to micro-computerized tomography images from a mouse lung obtained at four different static pressures to identify the static pressure-radius relationship for four generations of vessels in the pulmonary arterial network. A shape-fitting function was parameterized for each vessel in the network to characterize the nonlinear and heterogeneous nature of vessel distensibility in the pulmonary arteries. These data on morphometric and mechanical properties were used to simulate pressure and flow velocity propagation in the network using one-dimensional representations of fluid and vessel wall mechanics. Moreover, wave intensity analysis was used to study effects of wall mechanics on generation and propagation of pressure wave reflections. Simulations were conducted to investigate the role of linear versus nonlinear formulations of wall elasticity and homogeneous versus heterogeneous treatments of vessel wall properties. Accounting for heterogeneity, by parameterizing the pressure/distention equation of state individually for each vessel segment, was found to have little effect on the predicted pressure profiles and wave propagation compared to a homogeneous parameterization based on average behavior. However, substantially different results were obtained using a linear elastic thin-shell model than were obtained using a nonlinear model that has a more physiologically realistic pressure versus radius relationship.

During human heart failure, the balance of cardiac energy use switches from predominantly fatty acids (FAs) to glucose. We hypothesized that this substrate shift was the result of mitochondrial degeneration; therefore, we examined mitochondrial oxidation and ultrastructure in the failing human heart by using respirometry, transmission electron microscopy, and gene expression studies of demographically matched donor and failing human heart left ventricular (LV) tissues. Surprisingly, respiratory capacities for failing LV isolated mitochondria (n = 9) were not significantly diminished compared with donor LV isolated mitochondria (n = 7) for glycolysis (pyruvate + malate)- or FA (palmitoylcarnitine)-derived substrates, and mitochondrial densities, assessed via citrate synthase activity, were consistent between groups. Transmission electron microscopy images also showed no ultrastructural remodeling for failing vs. donor mitochondria; however, the fraction of lipid droplets (LDs) in direct contact with a mitochondrion was reduced, and the average distance between an LD and its nearest neighboring mitochondrion was increased. Analysis of FA processing gene expression between donor and failing LVs revealed 0.64-fold reduced transcript levels for the mitochondrial-LD tether, perilipin 5, in the failing myocardium (P = 0.003). Thus, reduced FA use in heart failure may result from improper delivery, potentially via decreased perilipin 5 expression and mitochondrial-LD tethering, and not from intrinsic mitochondrial dysfunction.-Holzem, K. M., Vinnakota, K. C., Ravikumar, V. K., Madden, E. J., Ewald, G. A., Dikranian, K., Beard, D. A., Efimov, I. R. Mitochondrial structure and function are not different between nonfailing donor and end-stage failing human hearts., (© FASEB.)

Despite extensive study over the past six decades the coupling of chemical reaction and mechanical processes in muscle dynamics is not well understood. We lack a theoretical description of how chemical processes (metabolite binding, ATP hydrolysis) influence and are influenced by mechanical processes (deformation and force generation). To address this need, a mathematical model of the muscle cross-bridge (XB) cycle based on Huxley's sliding filament theory is developed that explicitly accounts for the chemical transformation events and the influence of strain on state transitions. The model is identified based on elastic and viscous moduli data from mouse and rat myocardial strips over a range of perturbation frequencies, and MgATP and inorganic phosphate (Pi) concentrations. Simulations of the identified model reproduce the observed effects of MgATP and MgADP on the rate of force development. Furthermore, simulations reveal that the rate of force re-development measured in slack-restretch experiments is not directly proportional to the rate of XB cycling. For these experiments, the model predicts that the observed increase in the rate of force generation with increased Pi concentration is due to inhibition of cycle turnover by Pi. Finally, the model captures the observed phenomena of force yielding suggesting that it is a result of rapid detachment of stretched attached myosin heads., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

A control system model was developed to analyze data on in vivo coronary blood flow regulation and to probe how different mechanisms work together to control coronary flow from rest to exercise, and under a variety of experimental conditions, including cardiac pacing and with changes in coronary arterial pressure (autoregulation). In the model coronary flow is determined by the combined action of a feedback pathway signal that is determined by the level of plasma ATP in coronary venous blood, an adrenergic open-loop (feed-forward) signal that increases with exercise, and a contribution of pressure-mediated myogenic control. The model was identified based on data from exercise experiments where myocardial oxygen extraction, coronary flow, cardiac interstitial norepinephrine concentration, and arterial and coronary venous plasma ATP concentrations were measured during control and during adrenergic and purinergic receptor blockade conditions. The identified model was used to quantify the relative contributions of open-loop and feedback pathways and to illustrate the degree of redundancy in the control of coronary flow. The results indicate that the adrenergic open-loop control component is responsible for most of the increase in coronary blood flow that occurs during high levels of exercise. However, the adenine nucleotide-mediated metabolic feedback control component is essential. The model was evaluated by predicting coronary flow in cardiac pacing and autoregulation experiments with reasonable fits to the data. The analysis shows that a model in which coronary venous plasma adenine nucleotides are a signal in local metabolic feedback control of coronary flow is consistent with the available data., (Copyright © 2016 the American Physiological Society.)

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., (Copyright © 2016. Published by Elsevier Ltd.)

The Virtual Physiological Human (VPH) project aims to develop integrative, explanatory and predictive computational models (C-Models) as numerical investigational tools to study disease, identify and design effective therapies and provide an in silico platform for drug screening. Ultimately, these models rely on the analysis and integration of experimental data. As such, the success of VPH depends on the availability of physiologically realistic experimental models (E-Models) of human organ function that can be parametrized to test the numerical models. Here, the current state of suitable E-models, ranging from in vitro non-human cell organelles to in vivo human organ systems, is discussed. Specifically, challenges and recent progress in improving the physiological realism of E-models that may benefit the VPH project are highlighted and discussed using examples from the field of research on cardiovascular disease, musculoskeletal disorders, diabetes and Parkinson's disease.

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., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)