Search for: All records

Abstract SummaryMolecular mechanisms of biological functions and disease processes are exceptionally complex, and our ability to interrogate and understand relationships is becoming increasingly dependent on the use of computational modeling. We have developed “BioModME,” a standalone R-based web application package, providing an intuitive and comprehensive graphical user interface to help investigators build, solve, visualize, and analyze computational models of complex biological systems. Some important features of the application package include multi-region system modeling, custom reaction rate laws and equations, unit conversion, model parameter estimation utilizing experimental data, and import and export of model information in the Systems Biology Matkup Language format. The users can also export models to MATLAB, R, and Python languages and the equations to LaTeX and Mathematical Markup Language formats. Other important features include an online model development platform, multi-modality visualization tool, and efficient numerical solvers for differential-algebraic equations and optimization. Availability and implementationAll relevant software information including documentation and tutorials can be found at https://mcw.marquette.edu/biomedical-engineering/computational-systems-biology-lab/biomodme.php. Deployed software can be accessed at https://biomodme.ctsi.mcw.edu/. Source code is freely available for download at https://github.com/MCWComputationalBiologyLab/BioModME. 

ABSTRACT Voltage‐dependent anion channel (VDAC) is the primary conduit for regulated passage of ions and metabolites into and out of a mitochondrion. Calculating the solvation free energy for VDAC is crucial for understanding its stability, function, and interactions within the cellular environment. In this article, numerical schemes for computing the total solvation free energy for VDAC—comprising electrostatic, ideal gas, and excess free energies plus the nonpolar energy—are developed based on a nonuniform size modified Poisson–Boltzmann ion channel (nuSMPBIC) finite element solver along with tetrahedral meshes for VDAC proteins. The current mesh generation package is also updated to improve mesh quality and accelerate mesh generation. A VDAC Solvation Free Energy Calculation (VSFEC) package is then created by integrating these schemes with the updated mesh package, the nuSMPBIC finite element package, the PDB2PQR package, and the OPM database, as well as one uniform SMPBIC finite element package and one Poisson–Boltzmann ion channel (PBIC) finite element package. With the VSFEC package, many numerical experiments are made using six VDAC proteins, eight ionic solutions containing up to four ionic species, including ATP⁴⁻and Ca²⁺, two reference states, different boundary values, and different permittivity constants. The test results underscore the importance of considering nonuniform ionic size effects to explore the varying patterns of the total solvation free energy, and demonstrate the high performance of the VSFEC package for VDAC solvation free energy calculation. 

Abstract 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. 

We developed a computational model of sodium fluorescein (SF) biliary excretion in ex vivo machine perfusion and used this model to assess changes in model parameters associated with the activity of MRP2, a hepatocyte membrane transporter, in response to increasing warm ischemia time. We found a significant decrease in the parameter value describing MRP2 activity, consistent with a role of decreased MRP2 function in ischemia-reperfusion injury leading to decreased secretion of SF into bile. 

We assessed lung tissue mitochondrial bioenergetics in rats with tolerance (H-T) or susceptibility (H-S) to hyperoxia-induced ARDS. Results from studies in isolated mitochondria, tissue homogenate, and isolated perfused lungs show that mitochondrial bioenergetics are differentially altered in H-T and H-S lungs suggesting a potential role for mitochondrial bioenergetics in hyperoxia-induced ARDS. Results are clinically relevant since hyperoxia exposure is a primary therapy for patients with ARDS, and differential sensitivity to hyperoxia surely occurs in humans. 

Mitochondria are major sources of reactive oxygen species (ROS), which play important roles in both physiological and pathological processes. However, the specific contributions of different ROS production and scavenging components in the mitochondria of metabolically active tissues such as heart and kidney cortex and outer medulla (OM) are not well understood. Therefore, the goal of this study was to determine contributions of different ROS production and scavenging components and provide detailed comparisons of mitochondrial respiration, bioenergetics, ROS emission between the heart and kidney cortex and OM using tissues obtained from the same Sprague-Dawley rat under identical conditions and perturbations. Specifically, data were obtained using both NADH-linked substrate pyruvate + malate and FADH2-linked substrate succinate followed by additions of inhibitors of different components of the electron transport chain (ETC) and oxidative phosphorylation (OxPhos) and other ROS production and scavenging systems. Currently, there is limited data available for the mitochondria of kidney cortex and OM, the two major energy-consuming tissues in the body only next to the heart, and scarce quantitative information on the interplay between mitochondrial ROS production and scavenging systems in the three tissues. The findings from this study demonstrate significant differences in mitochondrial respiratory and bioenergetic functions and ROS emission among the three tissues. The results quantify the rates of ROS production from different complexes of the ETC, identify the complexes responsible for variations in mitochondrial membrane depolarization and regulations of ROS production, and quantify the contributions of ROS scavenging enzymes towards overall mitochondrial ROS emission. These findings advance our fundamental knowledge of tissue-specific and substrate-dependent mitochondrial respiratory and bioenergetic functions and ROS emission. This is important given the critical role that excess ROS production, oxidative stress, and mitochondrial dysfunction in the heart and kidney cortex and OM play in the pathogenesis of cardiovascular and renal diseases, including salt-sensitive hypertension. Keywords: Forward and reverse electron transfer; Mitochondrial metabolism; NADPH oxidase; Oxidative stress; ROS emission; ROS production and scavenging; Respiration and bioenergetics. 

Abstract— During liver transplantation, ischemia-reperfusion injury (IRI) is inevitable and decreases the overall success of the surgery. While guidelines exist, there is no reliable way to quantitatively assess the degree of IRI present in the liver. Our recent study has shown a correlation between the bile-to-plasma ratio of FDA-approved sodium fluorescein (SF) and the degree of hepatic IRI, presumably due to IRI-induced decrease in the activity of the hepatic multidrug resistance-associated protein 2 (MRP2); however, the contribution of SF blood clearance via the bile is still convoluted with other factors, such as renal clearance. In this work, we sought to computationally model SF blood clearance via the bile. First, we converted extant SF fluorescence data from rat whole blood, plasma, and bile to concentrations using calibration curves. Next, based on these SF concentration data, we generated a “liver-centric”, physiologically-based pharmacokinetic (PBPK) model of SF liver uptake and clearance via the bile. Model simulations show that SF bile concentration is highly sensitive to a change in the activity of hepatic MPR2. These simulations suggest that SF bile clearance along with the PBPK model can be used to quantify the effect of IRI on the activity of MRP2. Clinical Relevance— This study establishes the theory necessary to generate a model for predicting the degree of IRI during liver transplantation.