Modeling mitochondrial function

AW, Dash RK. Computational Modeling of Substrate-Dependent Mitochondrial Respiration and Bioenergetics in the Heart and Kidney Cortex and Outer Medulla

Integrated computational modeling provides a mechanistic and quantitative framework to characterize alterations in mitochondrial respiration and bioenergetics in response to different metabolic substrates in-silico. These alterations play critical roles in the pathogenesis of diseases affecting metabolically active organs such as heart and kidney. Therefore, the present study aimed to develop and validate thermodynamically constrained integrated computational models of mitochondrial respiration and bioenergetics in the heart and kidney cortex and outer medulla (OM). The models incorporated the kinetics of major biochemical reactions and transport processes as well as regulatory mechanisms in the mitochondria of these tissues. Intrinsic model parameters such as Michaelis-Menten constants were fixed at previously estimated values, while extrinsic model parameters such as maximal reaction and transport velocities were estimated separately for each tissue. This was achieved by fitting the model solutions to our recently published respirometry data measured in isolated rat heart and kidney cortex and OM mitochondria utilizing various NADH- and FADH2-linked metabolic substrates. The models were validated by predicting additional respirometry and bioenergetics data, which were not used for estimating the extrinsic model parameters. The models were able to predict tissue-specific and substrate-dependent mitochondrial emergent metabolic system properties such as redox states, enzyme and transporter fluxes, metabolite concentrations, membrane potential, and respiratory control index under diverse physiological and pathological conditions. The models were also able to quantitatively characterize differential regulations of NADH- and FADH2-linked metabolic pathways, which contribute differently toward regulations of oxidative phosphorylation and ATP synthesis in the heart and kidney cortex and OM mitochondria.

Hypoxia modeling techniques: A review

Hypoxia is the main cause and effect of a large number of diseases, including the most recent one facing the world, the coronavirus disease (COVID-19). Hypoxia is divided into short-term, long-term, and periodic, it can be the result of diseases, climate change, or living and traveling in the high mountain regions of the world. Since each type of hypoxia can be a cause and a consequence of various physiological changes, the methods for modeling these hypoxias are also different. There are many techniques for modeling hypoxia under experimental conditions. The most common animal for modeling hypoxia is a rat. Hypoxia models (hypoxia simulations) in rats are a tool to study the effect of various conditions on the oxygen supply of the body. These models can provide a necessary information to understand hypoxia and also provide effective treatment, highlighting the importance of various reactions of the body to hypoxia. The main parameters when choosing a model should be reproducibility and the goal that the scientist wants to achieve. Hypoxia in rats can be reproduced both ways exogenously and endogenously. The reason for writing this review was the aim to systematize the models of rats available in the literature in order to facilitate their selection by scientists. The relative strengths and limitations of each model need to be identified and understood in order to evaluate the information obtained from these models and extrapolate these results to humans to develop the necessary generalizations. Despite these problems, animal models have been and remain vital to understanding the mechanisms involved in the development and progression of hypoxia. The eligibility criteria for the selected studies was a comprehensive review of the methods and results obtained from the studies. This made it possible to make generalizations and give recommendations on the application of these methods. The review will assist scientists in choosing an appropriate hypoxia simulation method, as well as assist in interpreting the results obtained with these methods.

Mitochondrial ion channels in cardiac function

Mitochondria have been recognized as key organelles in cardiac physiology and are potential targets for clinical interventions to improve cardiac function. Mitochondrial dysfunction has been accepted as a major contributor to the development of heart failure. The main function of mitochondria is to meet the high energy demands of the heart by oxidative metabolism. Ionic homeostasis in mitochondria directly regulates oxidative metabolism, and any disruption in ionic homeostasis causes mitochondrial dysfunction and eventually contractile failure. The mitochondrial ionic homeostasis is closely coupled with inner mitochondrial membrane potential. To regulate and maintain ionic homeostasis, mitochondrial membranes are equipped with ion transporting proteins. Ion transport mechanisms involving several different ion channels and transporters are highly efficient and dynamic, thus helping to maintain the ionic homeostasis of ions as well as their salts present in the mitochondrial matrix. In recent years, several novel proteins have been identified on the mitochondrial membranes and these proteins are actively being pursued in research for roles in the organ as well as organelle physiology. In this article, the role of mitochondrial ion channels in cardiac function is reviewed. In recent times, the major focus of the mitochondrial ion channel field is to establish molecular identities as well as assigning specific functions to them. Given the diversity of mitochondrial ion channels and their unique roles in cardiac function, they present novel and viable therapeutic targets for cardiac diseases.

Energy metabolism design of the striated muscle cell

The design of the energy metabolism system in striated muscle remains a major area of investigation. Here, we review our current understanding and emerging hypotheses regarding the metabolic support of muscle contraction. Maintenance of ATP free energy, so called energy homeostasis, via mitochondrial oxidative phosphorylation is critical to sustained contractile activity, and this major design criterion is the focus of this review. Cell volume invested in mitochondria reduces the space available for generating contractile force, and this spatial balance between mitochondria acontractile elements to meet the varying sustained power demands across muscle types is another important design criterion. This is accomplished with remarkably similar mass-specific mitochondrial protein composition across muscle types, implying that it is the organization of mitochondria within the muscle cell that is critical to supporting sustained muscle function. Beyond the production of ATP, ubiquitous distribution of ATPases throughout the muscle requires rapid distribution of potential energy across these large cells. Distribution of potential energy has long been thought to occur primarily through facilitated metabolite diffusion, but recent analysis has questioned the importance of this process under normal physiological conditions. Recent structural and functional studies have supported the hypothesis that the mitochondrial reticulum provides a rapid energy distribution system via the conduction of the mitochondrial membrane potential to maintain metabolic homeostasis during contractile activity. We extensively review this aspect of the energy metabolism design contrasting it with metabolite diffusion models and how mitochondrial structure can play a role in the delivery of energy in the striated muscle.

Intracellular Chloride Channels: Novel Biomarkers in Diseases

Ion channels are integral membrane proteins present on the plasma membrane as well as intracellular membranes. In the human genome, there are more than 400 known genes encoding ion channel proteins. Ion channels are known to regulate several cellular, organellar, and physiological processes. Any mutation or disruption in their function can result in pathological disorders, both common or rare. Ion channels present on the plasma membrane are widely acknowledged for their role in various biological processes, but in recent years, several studies have pointed out the importance of ion channels located in intracellular organelles. However, ion channels located in intracellular organelles are not well-understood in the context of physiological conditions, such as the generation of cellular excitability and ionic homeostasis. Due to the lack of information regarding their molecular identity and technical limitations of studying them, intracellular organelle ion channels have thus far been overlooked as potential therapeutic targets. In this review, we focus on a novel class of intracellular organelle ion channels, Chloride Intracellular Ion Channels (CLICs), mainly documented for their role in cardiovascular, neurophysiology, and tumor biology. CLICs have a single transmembrane domain, and in cells, they exist in cytosolic as well as membranous forms. They are predominantly present in intracellular organelles and have recently been shown to be localized to cardiomyocyte mitochondria as well as exosomes. In fact, a member of this family, CLIC5, is the first mitochondrial chloride channel to be identified on the molecular level in the inner mitochondrial membrane, while another member, CLIC4, is located predominantly in the outer mitochondrial membrane. In this review, we discuss this unique class of intracellular chloride channels, their role in pathologies, such as cardiovascular, cancer, and neurodegenerative diseases, and the recent developments concerning their usage as theraputic targets.

BKCa (Slo) Channel Regulates Mitochondrial Function and Lifespan in Drosophila melanogaster

BK_Ca channels, originally discovered in Drosophila melanogaster as slowpoke (slo), are recognized for their roles in cellular and organ physiology. Pharmacological approaches implicated BK_Ca channels in cellular and organ protection possibly for their ability to modulate mitochondrial function. However, the direct role of BK_Ca channels in regulating mitochondrial structure and function is not deciphered. Here, we demonstrate that BK_Ca channels are present in fly mitochondria, and slo mutants show structural and functional defects in mitochondria. slo mutants display an increase in reactive oxygen species and the modulation of ROS affected their survival. We also found that the absence of BK_Ca channels reduced the lifespan of Drosophila, and overexpression of human BK_Ca channels in flies extends life span in males. Our study establishes the presence of BK_Ca channels in mitochondria of Drosophila and ascertains its novel physiological role in regulating mitochondrial structural and functional integrity, and lifespan.

Hesperidin Induces ROS-Mediated Apoptosis along with Cell Cycle Arrest at G2/M Phase in Human Gall Bladder Carcinoma

A natural predominant flavonoid hesperidin rich in citrus fruits exhibits multifunctional medicinal properties. The anticancerous potential of hesperidin has been widely explored; however, the gall bladder carcinoma (GBC) still remains untouched due to the unavailability of efficient experimental model. The aim of our study was to identify the apoptotic and antiproliferative potential of hesperidin in GBC. The promising efficacy of hesperidin was assessed through the generation of reactive oxygen species (ROS), cellular apoptosis, and loss of mitochondrial membrane potential (MMP) in the primary cells generated from surgically removed cancerous gall bladder tissues. Moreover, cell cycle analysis and caspases-3 activity were performed to confirm the apoptosis inducing potential of hesperidin. Results revealed that hesperidin exposure for 24 h at a dose of 200 µM reduced the cell proliferation of GBC cells significantly. In addition, hesperidin treatment further resulted in an increased ROS generation and nuclear condensation at the same dose. Caspase-3 activation and cell cycle arrest at G2/M phase were also accelerated in a dose-dependent manner. Together, these results suggest that hesperidin can be considered as a potential anticancerous compound for the treatment of GBC. Furthermore, evaluation of the pharmacological aspects of hesperidin is desirable for drug development.

Human Cardiac 31P-MR Spectroscopy at 3 Tesla Cannot Detect Failing Myocardial Energy Homeostasis during Exercise

Phosphorus-31 magnetic resonance spectroscopy (³¹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 ³¹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 ³¹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 ³¹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 ³¹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 ³¹P-MRS for diagnostic testing and risk stratification, and warrants developments in cardiac ³¹P-MRS exercise stress testing methodology.

No Muscle Is an Island: Integrative Perspectives on Muscle Fatigue

Muscle fatigue has been studied with a variety approaches, tools and technologies. The foci of these studies have ranged tremendously, from molecules to the entire organism. Single cell and animal models have been used to gain mechanistic insight into the fatigue process. The theme of this review is the concept that the mechanisms of muscle fatigue do not occur in isolation in vivo: muscular work is supported by many complex physiological systems, any of which could fail during exercise and thus contribute to fatigue. To advance our overall understanding of fatigue, a combination of models and approaches is necessary. In this review, we examine the roles that neuromuscular properties, intracellular glycogen, oxygen metabolism, and blood flow play in the fatigue process during exercise and pathological conditions.

Identification of active elementary flux modes in mitochondria using selectively permeabilized CHO cells

Metabolic compartmentation is a key feature of mammalian cells. Mitochondria are the powerhouse of eukaryotic cells, responsible for respiration and the TCA cycle. We accessed the mitochondrial metabolism of the economically important Chinese hamster ovary (CHO) cells using selective permeabilization. We tested key substrates without and with addition of ADP. Based on quantified uptake and production rates, we could determine the contribution of different elementary flux modes to the metabolism of a substrate or substrate combination. ADP stimulated the uptake of most metabolites, directly by serving as substrate for the respiratory chain, thus removing the inhibitory effect of NADH, or as allosteric effector. Addition of ADP favored substrate metabolization to CO2 and did not enhance the production of other metabolites. The controlling effect of ADP was more pronounced when we supplied metabolites to the first part of the TCA cycle: pyruvate, citrate, α-ketoglutarate and glutamine. In the second part of the TCA cycle, the rates were primarily controlled by the concentrations of C4-dicarboxylates. Without ADP addition, the activity of the pyruvate carboxylase-malate dehydrogenase-malic enzyme cycle consumed the ATP produced by oxidative phosphorylation, preventing its accumulation and maintaining metabolic steady state conditions. Aspartate was taken up only in combination with pyruvate, whose uptake also increased, a fact explained by complex regulatory effects. Isocitrate dehydrogenase and α-ketoglutarate dehydrogenase were identified as the key regulators of the TCA cycle, confirming existent knowledge from other cells. We have shown that selectively permeabilized cells combined with elementary mode analysis allow in-depth studying of the mitochondrial metabolism and regulation.

Optical spectroscopy in turbid media using an integrating sphere: mitochondrial chromophore analysis during metabolic transitions

Recent evidence suggests that the activity of mitochondrial oxidative phosphorylation complexes (MOPCs) is modulated at multiple sites. Here, a method of optically monitoring electron distribution within and between MOPCs is described using a center-mounted sample in an integrating sphere (to minimize scattering effects) with a rapid-scanning spectrometer. The redox-sensitive MOPC absorbances (∼465-630 nm) were modeled using linear least squares analysis with individual chromophore spectra. Classical mitochondrial activity transitions (e.g., ADP-induced increase in oxygen consumption) were used to characterize this approach. Most notable in these studies was the observation that intermediates of the catalytic cycle of cytochrome oxidase are dynamically modulated with metabolic state. The MOPC redox state, along with measurements of oxygen consumption and mitochondrial membrane potential, was used to evaluate the conductances of different sections of the electron transport chain. This analysis then was applied to mitochondria isolated from rabbit hearts subjected to ischemia/reperfusion (I/R). Surprisingly, I/R resulted in an inhibition of all measured MOPC conductances, suggesting a coordinated down-regulation of mitochondrial activity with this well-established cardiac perturbation.

Computational modeling of mitochondrial energy transduction

Mitochondria are the power plant of the heart, burning fat and sugars to supply the muscle with the adenosine triphosphate (ATP) free energy that drives contraction and relaxation during each heart beat. This function was first captured in a mathematical model in 1967. Today, interest in such a model has been rekindled by ongoing in silico integrative physiology efforts such as the Cardiac Physiome project. Here, the status of the field of computational modeling of mitochondrial ATP synthetic function is reviewed.

Prediction of muscle energy states at low metabolic rates requires feedback control of mitochondrial respiratory chain activity by inorganic phosphate

The regulation of the 100-fold dynamic range of mitochondrial ATP synthesis flux in skeletal muscle was investigated. Hypotheses of key control mechanisms were included in a biophysical model of oxidative phosphorylation and tested against metabolite dynamics recorded by (31)P nuclear magnetic resonance spectroscopy ((31)P MRS). Simulations of the initial model featuring only ADP and Pi feedback control of flux failed in reproducing the experimentally sampled relation between myoplasmic free energy of ATP hydrolysis (ΔG(p) = ΔG(p)(o')+RT ln ([ADP][Pi]/[ATP]) and the rate of mitochondrial ATP synthesis at low fluxes (<0.2 mM/s). Model analyses including Monte Carlo simulation approaches and metabolic control analysis (MCA) showed that this problem could not be amended by model re-parameterization, but instead required reformulation of ADP and Pi feedback control or introduction of additional control mechanisms (feed forward activation), specifically at respiratory Complex III. Both hypotheses were implemented and tested against time course data of phosphocreatine (PCr), Pi and ATP dynamics during post-exercise recovery and validation data obtained by (31)P MRS of sedentary subjects and track athletes. The results rejected the hypothesis of regulation by feed forward activation. Instead, it was concluded that feedback control of respiratory chain complexes by inorganic phosphate is essential to explain the regulation of mitochondrial ATP synthesis flux in skeletal muscle throughout its full dynamic range.

The effect of isolancifolide on the apoptosis in HL-60 cells through caspase-8-dependent and -independent pathways

Isolancifolide is a compound extracted and isolated from Actinodaphne lancifolia, which is a traditional oriental medicine. To determine whether isolancifolide has therapeutic potential as an anticancer molecule, we assessed its apoptotic effects on HL-60 cells, a human leukemia cell line. Apoptotic activities were investigated using DNA fragmentation assay, immunoblotting, and flow cytometry. We found that the inhibitory concentration 50% of isolancifolide was approximately 20 M. The time- and dose-dependent effects of isolancifolide on apoptosis were determined by DNA fragmentation and propidium iodide staining, and the involvement of caspases and the Bcl-2 family in isolancifolide-induced apoptosis was assessed by Western blotting. During exposure to isolancifolide, the pro-forms or full length of caspases-8, -3, and Bid were decreased, as assessed by Western blotting, while the levels of cleaved forms of caspases-8, -3, and PARP were increased. We observed that the release of cytochrome c and Smac/DIABLO from the mitochondria to the cytosol was accompanied by the loss of mitochondrial membrane potential. The caspase specific inhibitors, z-IETD-fmk and z-LEHD-fmk, blocked the accumulation of sub-G1 cells and the release of cytochrome c, but not that of Smac/DIABLO. These results indicate that isolancifolide induces apoptosis of HL-60 cells through both death receptor and mitochondria pathways, in caspase-8-dependent and -independent manners, suggesting that isolancifolide may be useful in anticancer strategies.

Mitochondrial functional specialization in glycolytic and oxidative muscle fibers: tailoring the organelle for optimal function

In skeletal muscle, two major types of muscle fibers exist: slow-twitch oxidative (type I) fibers designed for low-intensity long-lasting contractions, and fast-twitch glycolytic (type II) fibers designed for high-intensity short-duration contractions. Such a wide range of capabilities has emerged through the selection across fiber types of a narrow set of molecular characteristics suitable to achieve a specific contractile phenotype. In this article we review evidence supporting the existence of distinct functional phenotypes in mitochondria from slow and fast fibers that may be required to ensure optimal muscle function. This includes differences with respect to energy substrate preferences, regulation of oxidative phosphorylation, dynamics of reactive oxygen species, handling of Ca2+, and regulation of cell death. The potential physiological implications on muscle function and the putative mechanisms responsible for establishing and maintaining distinct mitochondrial phenotype across fiber types are also discussed.

Mechanisms of glucose sensing in the pancreatic β-cell: A computational systems-based analysis

Pancreatic β-cells respond to rising blood glucose by increasing oxidative metabolism, leading to an increased ATP/ADP ratio in the cytoplasm with a subsequent influx of calcium and the eventual secretion of insulin. The mechanisms of glucose sensing in the pancreatic β-cell involve the coupling of cytoplasmic and mitochondrial processes. Our analysis, based on mathematical models of data from multiple sources has implications for β-cell function and the treatment of type 2 diabetes (Fridlyand and Philipson, 2010). This β-cell glucose response model correctly predicts changes in the ATP/ADP ratio, cytoplasmic and mitochondrial calcium levels, and other metabolic parameters in response to alterations in substrate delivery at steady-state and during cytoplasmic calcium oscillations. Here we consider how peculiarities of β-cell pathways that result in dysfunction can be a consequence of specific mechanisms of glucose sensitivity, using our computational systems-based analysis. We found that the mitochondrial membrane potential must be relatively low in β-cells compared with other cell types to permit precise mitochondrial regulation of the cytoplasmic ATP/ADP ratio. This key difference may follow from a relative reduction in cellular respiratory activity. Our analysis additionally demonstrates how activity of lactate dehydrogenase, uncoupling proteins, and the redox shuttles all working in concert can regulate β-cell function. We further show that a decreased mitochondrial membrane potential may lead to a low rate of production of reactive oxygen species in β-cells under physiological conditions. This computational systems analysis aids in providing a more complete understanding of the complex process of β-cell glucose sensing.

Glucose sensing in the pancreatic beta cell: a computational systems analysis

Background

Pancreatic beta-cells respond to rising blood glucose by increasing oxidative metabolism, leading to an increased ATP/ADP ratio in the cytoplasm. This leads to a closure of K_ATP channels, depolarization of the plasma membrane, influx of calcium and the eventual secretion of insulin. Such mechanism suggests that beta-cell metabolism should have a functional regulation specific to secretion, as opposed to coupling to contraction. The goal of this work is to uncover contributions of the cytoplasmic and mitochondrial processes in this secretory coupling mechanism using mathematical modeling in a systems biology approach.

Methods

We describe a mathematical model of beta-cell sensitivity to glucose. The cytoplasmic part of the model includes equations describing glucokinase, glycolysis, pyruvate reduction, NADH and ATP production and consumption. The mitochondrial part begins with production of NADH, which is regulated by pyruvate dehydrogenase. NADH is used in the electron transport chain to establish a proton motive force, driving the F₁F₀ ATPase. Redox shuttles and mitochondrial Ca²⁺ handling were also modeled.

Results

The model correctly predicts changes in the ATP/ADP ratio, Ca²⁺ and other metabolic parameters in response to changes in substrate delivery at steady-state and during cytoplasmic Ca²⁺ oscillations. Our analysis of the model simulations suggests that the mitochondrial membrane potential should be relatively lower in beta cells compared with other cell types to permit precise mitochondrial regulation of the cytoplasmic ATP/ADP ratio. This key difference may follow from a relative reduction in respiratory activity. The model demonstrates how activity of lactate dehydrogenase, uncoupling proteins and the redox shuttles can regulate beta-cell function in concert; that independent oscillations of cytoplasmic Ca²⁺ can lead to slow coupled metabolic oscillations; and that the relatively low production rate of reactive oxygen species in beta-cells under physiological conditions is a consequence of the relatively decreased mitochondrial membrane potential.

Conclusion

This comprehensive model predicts a special role for mitochondrial control mechanisms in insulin secretion and ROS generation in the beta cell. The model can be used for testing and generating control hypotheses and will help to provide a more complete understanding of beta-cell glucose-sensing central to the physiology and pathology of pancreatic β-cells.

Mitochondria in cardiomyocyte Ca2+ signaling

Ca(2+) signaling is of vital importance to cardiac cell function and plays an important role in heart failure. It is based on sarcolemmal, sarcoplasmic reticulum and mitochondrial Ca(2+) cycling. While the first two are well characterized, the latter remains unclear, controversial and technically challenging. In mammalian cardiac myocytes, Ca(2+) influx through L-type calcium channels in the sarcolemmal membrane triggers Ca(2+) release from the nearby junctional sarcoplasmic reticulum to produce Ca(2+) sparks. When this triggering is synchronized by the cardiac action potential, a global [Ca(2+)](i) transient arises from coordinated Ca(2+) release events. The ends of intermyofibrillar mitochondria are located within 20 nm of the junctional sarcoplasmic reticulum and thereby experience a high local [Ca(2+)] during the Ca(2+) release process. Both local and global Ca(2+) signals may thus influence calcium signaling in mitochondria and, reciprocally, mitochondria may contribute to the local control of calcium signaling. In addition to the intermyofibrillar mitochondria, morphologically distinct mitochondria are also located in the perinuclear and subsarcolemmal regions of the cardiomyocyte and thus experience a different local [Ca(2+)]. Here we review the literature in regard to several issues of broad interest: (1) the ultrastructural basis for mitochondrion - sarcoplasmic reticulum cross-signaling; (2) mechanisms of sarcoplasmic reticulum signaling; (3) mitochondrial calcium signaling; and (4) the possible interplay of calcium signaling between the sarcoplasmic reticulum and adjacent mitochondria. Finally, this review discusses experimental findings and mathematical models of cardiac calcium signaling between the sarcoplasmic reticulum and mitochondria, identifies weaknesses in these models, and suggests strategies and approaches for future investigations.

Mitochondria as ATP consumers in cellular pathology

ATP provided by oxidative phosphorylation supports highly complex and energetically expensive cellular processes. Yet, in several pathological settings, mitochondria could revert to ATP consumption, aggravating an existing cellular pathology. Here we review (i) the pathological conditions leading to ATP hydrolysis by the reverse operation of the mitochondrial F(o)F(1)-ATPase, (ii) molecular and thermodynamic factors influencing the directionality of the F(o)F(1)-ATPase, (iii) the role of the adenine nucleotide translocase as the intermediary adenine nucleotide flux pathway between the cytosol and the mitochondrial matrix when mitochondria become ATP consumers, (iv) the role of the permeability transition pore in bypassing the ANT, thereby allowing the flux of ATP directly to the hydrolyzing F(o)F(1)-ATPase, (v) the impact of the permeability transition pore on glycolytic ATP production, and (vi) endogenous and exogenous interventions for limiting ATP hydrolysis by the mitochondrial F(o)F(1)-ATPase.

Gallic acid ester derivatives induce apoptosis and cell adhesion inhibition in melanoma cells: The relationship between free radical generation, glutathione depletion and cell death

Malignant melanoma is a lethal disease, and the incidence and mortality associated with it are increasing worldwide. It has a significant tendency to develop both metastasis and resistance to chemotherapy. The tumor cells show abnormal redox regulation, and although the molecular mechanisms involved are not well characterized, they seem to be related to oxidative stress. In a previous study, we showed the antitumoral properties of gallic acid ester derivatives in leukemia cells. Here, we show the effect of octyl, decyl, dodecyl and tetradecyl gallates on B16F10 cells, a melanoma cell line. All compounds induced cytotoxic effects, and the IC(50) values obtained were between 7microM and 17microM after 48h of incubation. Cell death occurred through apoptosis, as demonstrated by the genomic DNA fragmentation pattern. The gallates were able to induce significant production of free radicals, deplete both glutathione and ATP, activate NF-kappaB and promote the inhibition of cell adhesion under the experimental conditions. The glutathione depletion induced by these compounds was related to the inhibition of gamma-glutamylcysteine synthase activity. These results suggest that gallates induce tumoral cell death through apoptosis as a consequence of oxidative stress, though they use different mechanisms to do so. These findings are important since melanoma cells are resistant to death because of their high level of antioxidant defense, adhesion capability and propensity to metastasize.

Alterations of mitochondrial enzymes contribute to cardiac hypertrophy before hypertension development in spontaneously hypertensive rats

Mitochondrial dysfunction is recently thought to be tightly associated with the development of cardiac hypertrophy as well as hypertension. However, the detailed molecular events in mitochondria at early stages of hypertrophic pathogenesis are still unclear. Applying two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) combined with MALDI-TOF/TOF tandem mass spectrometry, here we identified the changed mitochondrial proteins of left ventricular mitochondria in prehypertensive/hypertensive stages of cardiac hypertrophy through comparing spontaneously hypertensive rats (SHR) and the age-matched normotensive Wistar Kyoto (WKY) rats. The results revealed that in the hypertrophic left ventricle of SHR as early as 4 weeks old with normal blood pressure, 33 mitochondrial protein spots presented significant alterations, with 17 down-regulated and 16 up-regulated. Such alterations were much greater than those in 20-week-old SHR with elevated blood pressure. Of the total alterations, the expression of two mitochondrial enzymes, trifunctional enzyme alpha subunit (Hadha) and NADH dehydrogenase 1 alpha subcomplex 10 (Ndufa10), were found to have special expression modification patterns in SHR strain. These data would provide new clues to investigate the potential contribution of mitochondrial dysfunction to the development of cardiac hypertrophy.

Roles of the creatine kinase system and myoglobin in maintaining energetic state in the working heart

BACKGROUND

The heart is capable of maintaining contractile function despite a transient decrease in blood flow and increase in cardiac ATP demand during systole. This study analyzes a previously developed model of cardiac energetics and oxygen transport to understand the roles of the creatine kinase system and myoglobin in maintaining the ATP hydrolysis potential during beat-to-beat transient changes in blood flow and ATP hydrolysis rate.

RESULTS

The theoretical investigation demonstrates that elimination of myoglobin only slightly increases the predicted range of oscillation of cardiac oxygenation level during beat-to-beat transients in blood flow and ATP utilization. In silico elimination of myoglobin has almost no impact on the cytoplasmic ATP hydrolysis potential (DeltaGATPase). In contrast, disabling the creatine kinase system results in considerable oscillations of cytoplasmic ADP and ATP levels and seriously deteriorates the stability of DeltaGATPase in the beating heart.

CONCLUSION

The CK system stabilizes DeltaGATPase by both buffering ATP and ADP concentrations and enhancing the feedback signal of inorganic phosphate in regulating mitochondrial oxidative phosphorylation.

3D imaging of mammalian cells with ion-abrasion scanning electron microscopy

Understanding the hierarchical organization of molecules and organelles within the interior of large eukaryotic cells is a challenge of fundamental interest in cell biology. We are using ion-abrasion scanning electron microscopy (IA-SEM) to visualize this hierarchical organization in an approach that combines focused ion-beam milling with scanning electron microscopy. Here, we extend our previous studies on imaging yeast cells to image subcellular architecture in human melanoma cells and melanocytes at resolutions as high as approximately 6 and approximately 20 nm in the directions parallel and perpendicular, respectively, to the direction of ion-beam milling. The 3D images demonstrate the striking spatial relationships between specific organelles such as mitochondria and membranes of the endoplasmic reticulum, and the distribution of unique cellular components such as melanosomes. We also show that 10nm-sized gold particles and quantum dot particles with 7 nm-sized cores can be detected in single cross-sectional images. IA-SEM is thus a useful tool for imaging large mammalian cells in their entirety at resolutions in the nanometer range.

Nanovehicular intracellular delivery systems

This article provides an overview of principles and barriers relevant to intracellular drug and gene transport, accumulation and retention (collectively called as drug delivery) by means of nanovehicles (NV). The aim is to deliver a cargo to a particular intracellular site, if possible, to exert a local action. Some of the principles discussed in this article apply to noncolloidal drugs that are not permeable to the plasma membrane or to the blood-brain barrier. NV are defined as a wide range of nanosized particles leading to colloidal objects which are capable of entering cells and tissues and delivering a cargo intracelullarly. Different localization and targeting means are discussed. Limited discussion on pharmacokinetics and pharmacodynamics is also presented. NVs are contrasted to micro-delivery and current nanotechnologies which are already in commercial use. Newer developments in NV technologies are outlined and future applications are stressed. We also briefly review the existing modeling tools and approaches to quantitatively describe the behavior of targeted NV within the vascular and tumor compartments, an area of particular importance. While we list "elementary" phenomena related to different level of complexity of delivery to cancer, we also stress importance of multi-scale modeling and bottom-up systems biology approach.

Modelling of cells bioenergetics

This paper presents an integrated model describing the control of Saccharomyces cerevisiae yeast cells bioenergetics. This model describes the oxidative and respirofermentative metabolism. The model assumes that the mitochondria of the Saccharomyces cerevisiae cells are charged with NADH during the tricarboxylic acid cycle, and NADH is discharged from mitochondria later in the electron transport system. Selected effects observed in the Saccharomyces cerevisiae eucaryotic cells, including the Pasteur's and Crabtree effects, are also modeled.

Modeling cellular metabolism and energetics in skeletal muscle: large-scale parameter estimation and sensitivity analysis

Skeletal muscle plays a major role in the regulation of whole-body energy metabolism during physiological stresses such as ischemia, hypoxia, and exercise. Current experimental techniques provide relatively little in vivo data on dynamic responses of metabolite concentrations and metabolic fluxes in skeletal muscle to such physiological stimuli. As a complementary approach to experimental measurements and as a framework for quantitatively analyzing available in vivo data, a physiologically based model of skeletal muscle cellular metabolism and energetics is developed. This model, which incorporates key transport and reaction processes, is based on dynamic mass balances of 30 chemical species in capillary (blood) and tissue (cell) domains. The reaction fluxes in the cellular domain are expressed in terms of a generalized Michaelis?Menten equation involving energy controller ratios ATP/ADP and ATP/ADP and NADH/NAD+ . This formalism introduces a large number of unknown parameters ( approximately 90). Estimating these parameters from in vivo sparse data and evaluating dynamic sensitivities of the model outputs with respect to these parameters is a challenging problem. Parameter estimation is accomplished using an efficient, nonlinear, constraint-based, optimization algorithm that minimizes differences between available experimental data and corresponding model outputs by explicitly utilizing equality constraints on resting fluxes and concentrations. With the estimated parameter values, the model is able to simulate dynamic responses to reduced blood flow (ischemia) of key metabolite concentrations and metabolic fluxes, both measured and nonmeasured. A general parameter sensitivity analysis is carried out to determine and characterize the parameters having the most and least effects on the measured outputs.

A dynamic model of nitric oxide inhibition of mitochondrial cytochrome c oxidase

Nitric oxide can inhibit mitochondrial cytochrome oxidase in both oxygen competitive and uncompetitive modes. A previous model described these interactions assuming equilibrium binding to the reduced and oxidised enzyme respectively (Mason, et al. Proc. Natl. Acad. Sci. U S A 103 (2006) 708-713). Here we demonstrate that the equilibrium assumption is inappropriate as it requires unfeasibly high association constants for NO to the oxidised enzyme. Instead we develop a model which explicitly includes NO binding and its enzyme-bound conversion to nitrite. Removal of the nitrite complex requires electron transfer to the binuclear centre from haem a. This revised model fits the inhibition constants at any value of substrate concentration (ferrocytochrome c or oxygen). It predicts that the inhibited steady state should be a mixture of the reduced haem nitrosyl complex and the oxidized-nitrite complex. Unlike the previous model, binding to the oxidase is always proportional to the degree of inhibition of oxygen consumption. The model is consistent with data and models from a recent paper suggesting that the primary effect of NO binding to the oxidised enzyme is to convert NO to nitrite, rather than to inhibit enzyme activity (Antunes et al. Antioxid. Redox Signal. 9 (2007) 1569-1579).

Synthesis and biological evaluation of thiobenzanilides as anticancer agents

A series of novel thiobenzanilides is described. These compounds have been previously found to show strong biological activity such as antimycotic and antifungal actions. This is the first demonstration on the mechanism of the anticancer effect of thiobenzanilide agents (4a-c) on human melanoma A375 cells. The cytotoxic studies of compounds 4a-c on human melanoma A375 cells indicate thiobenzanilides induced higher cytotoxicity than nitrobenzanilides (3a-c). In addition, DNA flow cytometric analysis shows that 4a-c displays a significant G2/M phase arrest, which progresses to early apoptosis as detected by flow cytometry after double-staining with annexin V and propidium iodide (PI). Because cellular apoptosis is often preceded by the disruption of mitochondrial function, the assessment of mitochondrial function in 4a-c-treated cells is worthy of investigation. Our data revealed that treatment of A375 cells with 4a-c resulted in the loss of mitochondrial membrane potential (DeltaPsi(mt)), a reduction of ATP synthesis, increased reactive oxygen species (ROS) generation, and activation of caspase-3. Thus, we suggest that 4a-c agents are potent inducers of cell apoptosis in A375 cells.

Model of oxygen transport and metabolism predicts effect of hyperoxia on canine muscle oxygen uptake dynamics

Previous studies have shown that increased oxygen delivery, via increased convection or arterial oxygen content, does not speed the dynamics of oxygen uptake, V̇o2m, in dog muscle electrically stimulated at a submaximal metabolic rate. However, the dynamics of transport and metabolic processes that occur within working muscle in situ is typically unavailable in this experimental setting. To investigate factors affecting V̇o2m dynamics at contraction onset, we combined dynamic experimental data across working muscle with a mechanistic model of oxygen transport and metabolism in muscle. The model is based on dynamic mass balances for O2, ATP, and PCr. Model equations account for changes in cellular ATPase, oxidative phosphorylation, and creatine kinase fluxes in skeletal muscle during exercise, and cellular respiration depends on [ADP] and [O2]. Model simulations were conducted at different levels of arterial oxygen content and blood flow to quantify the effects of convection and diffusion of oxygen on the regulation of cellular respiration during step transitions from rest to isometric contraction in dog gastrocnemius muscle. Simulations of arteriovenous O2 differences and V̇o2m dynamics were successfully compared with experimental data (Grassi B, Gladden LB, Samaja M, Stary CM, Hogan MC. J Appl Physiol 85: 1394–1403, 1998; and Grassi B, Gladden LB, Stary CM, Wagner PD, Hogan MC. J Appl Physiol 85: 1404–1412, 1998), thus demonstrating the validity of the model, as well as its predictive capability. The main findings of this study are: 1) the estimated dynamic response of oxygen utilization at contraction onset in muscle is faster than that of oxygen uptake; and 2) hyperoxia does not accelerate the dynamics of diffusion and consequently muscle oxygen uptake at contraction onset due to the hyperoxia-induced increase in oxygen stores. These in silico derived results cannot be obtained from experimental observations alone.