Physiological Indirect Response Model to Omics-Powered Quantitative Systems Pharmacology Model

Over the past several decades, mathematical modeling has been applied to increasingly wider scopes of questions in drug development. Accordingly, the range of modeling tools has also been evolving, as showcased by contributions of Jusko and colleagues: from basic pharmacokinetics/pharmacodynamics (PK/PD) modeling to today's platform-based approach of quantitative systems pharmacology (QSP) modeling. Aimed at understanding the mechanism of action of investigational drugs, QSP models characterize systemic effects by incorporating information about cellular signaling networks, which is often represented by omics data. In this perspective, we share a few examples illustrating approaches for the integration of omics into mechanistic QSP modeling. We briefly overview how the evolution of PK/PD modeling into QSP has been accompanied by an increase in available data and the complexity of mathematical methods that integrate it. We discuss current gaps and challenges of integrating omics data into QSP models and propose several potential areas where integrated QSP and omics modeling may benefit drug development.

Predictive analytics identifies key factors driving hyperalgesic priming of muscle sensory neurons

Hyperalgesic priming, a form of neuroplasticity induced by inflammatory mediators, in peripheral nociceptors enhances the magnitude and duration of action potential (AP) firing to future inflammatory events and can potentially lead to pain chronification. The mechanisms underlying the development of hyperalgesic priming are not well understood, limiting the identification of novel therapeutic strategies to combat chronic pain. In this study, we used a computational model to identify key proteins whose modifications caused priming of muscle nociceptors and made them hyperexcitable to a subsequent inflammatory event. First, we extended a previously validated model of mouse muscle nociceptor sensitization to incorporate Epac-mediated interaction between two G protein-coupled receptor signaling pathways commonly activated by inflammatory mediators. Next, we calibrated and validated the model simulations of the nociceptor's AP response to both innocuous and noxious levels of mechanical force after two subsequent inflammatory events using literature data. Then, by performing global sensitivity analyses that simulated thousands of nociceptor-priming scenarios, we identified five ion channels and two molecular processes (from the 18 modeled transmembrane proteins and 29 intracellular signaling components) as potential regulators of the increase in AP firing in response to mechanical forces. Finally, when we simulated specific neuroplastic modifications in Kv1.1 and Nav1.7 alone as well as with simultaneous modifications in Nav1.7, Nav1.8, TRPA1, and Kv7.2, we observed a considerable increase in the fold change in the number of triggered APs in primed nociceptors. These results suggest that altering the expression of Kv1.1 and Nav1.7 might regulate the neuronal hyperexcitability in primed mechanosensitive muscle nociceptors.

Identification of key factors driving inflammation-induced sensitization of muscle sensory neurons

Sensory neurons embedded in muscle tissue that initiate pain sensations, i.e., nociceptors, are temporarily sensitized by inflammatory mediators during musculoskeletal trauma. These neurons transduce peripheral noxious stimuli into an electrical signal [i.e., an action potential (AP)] and, when sensitized, demonstrate lower activation thresholds and a heightened AP response. We still do not understand the relative contributions of the various transmembrane proteins and intracellular signaling processes that drive the inflammation-induced hyperexcitability of nociceptors. In this study, we used computational analysis to identify key proteins that could regulate the inflammation-induced increase in the magnitude of AP firing in mechanosensitive muscle nociceptors. First, we extended a previously validated model of a mechanosensitive mouse muscle nociceptor to incorporate two inflammation-activated G protein-coupled receptor (GPCR) signaling pathways and validated the model simulations of inflammation-induced nociceptor sensitization using literature data. Then, by performing global sensitivity analyses that simulated thousands of inflammation-induced nociceptor sensitization scenarios, we identified three ion channels and four molecular processes (from the 17 modeled transmembrane proteins and 28 intracellular signaling components) as potential regulators of the inflammation-induced increase in AP firing in response to mechanical forces. Moreover, we found that simulating single knockouts of transient receptor potential ankyrin 1 (TRPA1) and reducing the rates of G_αq-coupled receptor phosphorylation and G_αq subunit activation considerably altered the excitability of nociceptors (i.e., each modification increased or decreased the inflammation-induced fold change in the number of triggered APs compared to when all channels were present). These results suggest that altering the expression of TRPA1 or the concentration of intracellular G_αq might regulate the inflammation-induced increase in AP response of mechanosensitive muscle nociceptors.

Metabolic responses in blood-stage malaria parasites associated with increased and decreased sensitivity to PfATP4 inhibitors

BACKGROUND

Spiroindolone and pyrazoleamide antimalarial compounds target Plasmodium falciparum P-type ATPase (PfATP4) and induce disruption of intracellular Na⁺ homeostasis. Recently, a PfATP4 mutation was discovered that confers resistance to a pyrazoleamide while increasing sensitivity to a spiroindolone. Transcriptomic and metabolic adaptations that underlie this seemingly contradictory response of P. falciparum to sublethal concentrations of each compound were examined to understand the different cellular accommodation to PfATP4 disruptions.

METHODS

A genetically engineered P. falciparum Dd2 strain (Dd2^A211V) carrying an Ala211Val (A211V) mutation in PfATP4 was used to identify metabolic adaptations associated with the mutation that results in decreased sensitivity to PA21A092 (a pyrazoleamide) and increased sensitivity to KAE609 (a spiroindolone). First, sublethal doses of PA21A092 and KAE609 causing substantial reduction (30-70%) in Dd2^A211V parasite replication were identified. Then, at this sublethal dose of PA21A092 (or KAE609), metabolomic and transcriptomic data were collected during the first intraerythrocytic developmental cycle. Finally, the time-resolved data were integrated with a whole-genome metabolic network model of P. falciparum to characterize antimalarial-induced physiological adaptations.

RESULTS

Sublethal treatment with PA21A092 caused significant (p < 0.001) alterations in the abundances of 91 Plasmodium gene transcripts, whereas only 21 transcripts were significantly altered due to sublethal treatment with KAE609. In the metabolomic data, a substantial alteration (≥ fourfold) in the abundances of carbohydrate metabolites in the presence of either compound was found. The estimated rates of macromolecule syntheses between the two antimalarial-treated conditions were also comparable, except for the rate of lipid synthesis. A closer examination of parasite metabolism in the presence of either compound indicated statistically significant differences in enzymatic activities associated with synthesis of phosphatidylcholine, phosphatidylserine, and phosphatidylinositol.

CONCLUSION

The results of this study suggest that malaria parasites activate protein kinases via phospholipid-dependent signalling in response to the ionic perturbation induced by the Na⁺ homeostasis disruptor PA21A092. Therefore, targeted disruption of phospholipid signalling in PA21A092-resistant parasites could be a means to block the emergence of resistance to PA21A092.

Metabolic changes accompanying the loss of fumarate hydratase and malate-quinone oxidoreductase in the asexual blood stage of Plasmodium falciparum

In the glucose-rich milieu of red blood cells, asexually replicating malarial parasites mainly rely on glycolysis for ATP production, with limited carbon flux through the mitochondrial tricarboxylic acid (TCA) cycle. By contrast, gametocytes and mosquito-stage parasites exhibit an increased dependence on the TCA cycle and oxidative phosphorylation for more economical energy generation. Prior genetic studies supported these stage-specific metabolic preferences by revealing that six of eight TCA cycle enzymes are completely dispensable during the asexual blood stages of Plasmodium falciparum, with only fumarate hydratase (FH) and malate-quinone oxidoreductase (MQO) being refractory to deletion. Several hypotheses have been put forth to explain the possible essentiality of FH and MQO, including their participation in a malate shuttle between the mitochondrial matrix and the cytosol. However, using newer genetic techniques like CRISPR and dimerizable Cre, we were able to generate deletion strains of FH and MQO in P. falciparum. We employed metabolomic analyses to characterize a double knockout mutant of FH and MQO (ΔFM) and identified changes in purine salvage and urea cycle metabolism that may help to limit fumarate accumulation. Correspondingly, we found that the ΔFM mutant was more sensitive to exogenous fumarate, which is known to cause toxicity by modifying and inactivating proteins and metabolites. Overall, our data indicate that P. falciparum is able to adequately compensate for the loss of FH and MQO, rendering them unsuitable targets for drug development.

Metabolic adjustments of blood-stage Plasmodium falciparum in response to sublethal pyrazoleamide exposure

Due to the recurring loss of antimalarial drugs to resistance, there is a need for novel targets, drugs, and combination therapies to ensure the availability of current and future countermeasures. Pyrazoleamides belong to a novel class of antimalarial drugs that disrupt sodium ion homeostasis, although the exact consequences of this disruption in Plasmodium falciparum remain under investigation. In vitro experiments demonstrated that parasites carrying mutations in the metabolic enzyme PfATP4 develop resistance to pyrazoleamide compounds. However, the underlying mechanisms that allow mutant parasites to evade pyrazoleamide treatment are unclear. Here, we first performed experiments to identify the sublethal dose of a pyrazoleamide compound (PA21A092) that caused a significant reduction in growth over one intraerythrocytic developmental cycle (IDC). At this drug concentration, we collected transcriptomic and metabolomic data at multiple time points during the IDC to quantify gene- and metabolite-level alterations in the treated parasites. To probe the effects of pyrazoleamide treatment on parasite metabolism, we coupled the time-resolved omics data with a metabolic network model of P. falciparum. We found that the drug-treated parasites adjusted carbohydrate metabolism to enhance synthesis of myoinositol-a precursor for phosphatidylinositol biosynthesis. This metabolic adaptation caused a decrease in metabolite flux through the pentose phosphate pathway, causing a decreased rate of RNA synthesis and an increase in oxidative stress. Our model analyses suggest that downstream consequences of enhanced myoinositol synthesis may underlie adjustments that could lead to resistance emergence in P. falciparum exposed to a sublethal dose of a pyrazoleamide drug.

In silico Identification of Key Factors Driving the Response of Muscle Sensory Neurons to Noxious Stimuli

Nociceptive nerve endings embedded in muscle tissue transduce peripheral noxious stimuli into an electrical signal [i.e., an action potential (AP)] to initiate pain sensations. A major contributor to nociception from the muscles is mechanosensation. However, due to the heterogeneity in the expression of proteins, such as ion channels, pumps, and exchangers, on muscle nociceptors, we currently do not know the relative contributions of different proteins and signaling molecules to the neuronal response due to mechanical stimuli. In this study, we employed an integrated approach combining a customized experimental study in mice with a computational model to identify key proteins that regulate mechanical nociception in muscles. First, using newly collected data from somatosensory recordings in mouse hindpaw muscles, we developed and then validated a computational model of a mechanosensitive mouse muscle nociceptor. Next, by performing global sensitivity analyses that simulated thousands of nociceptors, we identified three ion channels (among the 17 modeled transmembrane proteins and four endoplasmic reticulum proteins) as potential regulators of the nociceptor response to mechanical forces in both the innocuous and noxious range. Moreover, we found that simulating single knockouts of any of the three ion channels, delayed rectifier voltage-gated K⁺ channel (Kv1.1) or mechanosensitive channels Piezo2 or TRPA1, considerably altered the excitability of the nociceptor (i.e., each knockout increased or decreased the number of triggered APs compared to when all channels were present). These results suggest that altering expression of the gene encoding Kv1.1, Piezo2, or TRPA1 might regulate the response of mechanosensitive muscle nociceptors.

Inter-study and time-dependent variability of metabolite abundance in cultured red blood cells

BACKGROUND

Cultured human red blood cells (RBCs) provide a powerful ex vivo assay platform to study blood-stage malaria infection and propagation. In recent years, high-resolution metabolomic methods have quantified hundreds of metabolites from parasite-infected RBC cultures under a variety of perturbations. In this context, the corresponding control samples of the uninfected culture systems can also be used to examine the effects of these perturbations on RBC metabolism itself and their dependence on blood donors (inter-study variations).

METHODS

Time-course datasets from five independent studies were generated and analysed, maintaining uninfected RBCs (uRBC) at 2% haematocrit for 48 h under conditions originally designed for parasite cultures. Using identical experimental protocols, quadruplicate samples were collected at six time points, and global metabolomics were employed on the pellet fraction of the uRBC cultures. In total, ~ 500 metabolites were examined across each dataset to quantify inter-study variability in RBC metabolism, and metabolic network modelling augmented the analyses to characterize the metabolic state and fluxes of the RBCs.

RESULTS

To minimize inter-study variations unrelated to RBC metabolism, an internal standard metabolite (phosphatidylethanolamine C18:0/20:4) was identified with minimal variation in abundance over time and across all the samples of each dataset to normalize the data. Although the bulk of the normalized data showed a high degree of inter-study consistency, changes and variations in metabolite levels from individual donors were noted. Thus, a total of 24 metabolites were associated with significant variation in the 48-h culture time window, with the largest variations involving metabolites in glycolysis and synthesis of glutathione. Metabolic network analysis was used to identify the production of superoxide radicals in cultured RBCs as countered by the activity of glutathione oxidoreductase and synthesis of reducing equivalents via the pentose phosphate pathway. Peptide degradation occurred at a rate that is comparable with central carbon fluxes, consistent with active degradation of methaemoglobin, processes also commonly associated with storage lesions in RBCs.

CONCLUSIONS

The bulk of the data showed high inter-study consistency. The collected data, quantification of an expected abundance variation of RBC metabolites, and characterization of a subset of highly variable metabolites in the RBCs will help in identifying non-specific changes in metabolic abundances that may obscure accurate metabolomic profiling of Plasmodium falciparum and other blood-borne pathogens.

Metabolic alterations in the erythrocyte during blood-stage development of the malaria parasite

Background

Human blood cells (erythrocytes) serve as hosts for the malaria parasite Plasmodium falciparum during its 48-h intraerythrocytic developmental cycle (IDC). Established in vitro protocols allow for the study of host–parasite interactions during this phase and, in particular, high-resolution metabolomics can provide a window into host–parasite interactions that support parasite development.

Methods

Uninfected and parasite-infected erythrocyte cultures were maintained at 2% haematocrit for the duration of the IDC, while parasitaemia was maintained at 7% in the infected cultures. The parasite-infected cultures were synchronized to obtain stage-dependent information of parasite development during the IDC. Samples were collected in quadruplicate at six time points from the uninfected and parasite-infected cultures and global metabolomics was used to analyse cell fractions of these cultures.

Results

In uninfected and parasite-infected cultures during the IDC, 501 intracellular metabolites, including 223 lipid metabolites, were successfully quantified. Of these, 19 distinct metabolites were present only in the parasite-infected culture, 10 of which increased to twofold in abundance during the IDC. This work quantified approximately five times the metabolites measured in previous studies of similar research scope, which allowed for more detailed analyses. Enrichment in lipid metabolism pathways exhibited a time-dependent association with different classes of lipids during the IDC. Specifically, enrichment occurred in sphingolipids at the earlier stages, and subsequently in lysophospholipid and phospholipid metabolites at the intermediate and end stages of the IDC, respectively. In addition, there was an accumulation of 18-, 20-, and 22-carbon polyunsaturated fatty acids, which produce eicosanoids and promote gametocytogenesis in infected erythrocyte cultures.

Conclusions

The current study revealed a number of heretofore unidentified metabolic components of the host–parasite system, which the parasite may exploit in a time-dependent manner to grow over the course of its development in the blood stage. Notably, the analyses identified components, such as precursors of immunomodulatory molecules, stage-dependent lipid dynamics, and metabolites, unique to parasite-infected cultures. These conclusions are reinforced by the metabolic alterations that were characterized during the IDC, which were in close agreement with those known from previous studies of blood-stage infection.

A mevalonate bypass system facilitates elucidation of plastid biology in malaria parasites

Malaria parasites rely on a plastid organelle for survival during the blood stages of infection. However, the entire organelle is dispensable as long as the isoprenoid precursor, isopentenyl pyrophosphate (IPP), is supplemented in the culture medium. We engineered parasites to produce isoprenoid precursors from a mevalonate-dependent pathway, creating a parasite line that replicates normally after the loss of the apicoplast organelle. We show that carbon-labeled mevalonate is specifically incorporated into isoprenoid products, opening new avenues for researching this essential class of metabolites in malaria parasites. We also show that essential apicoplast proteins, such as the enzyme target of the drug fosmidomycin, can be deleted in this mevalonate bypass parasite line, providing a new method to determine the roles of other important apicoplast-resident proteins. Several antibacterial drugs kill malaria parasites by targeting basic processes, such as transcription, in the organelle. We used metabolomic and transcriptomic methods to characterize parasite metabolism after azithromycin treatment triggered loss of the apicoplast and found that parasite metabolism and the production of apicoplast proteins is largely unaltered. These results provide insight into the effects of apicoplast-disrupting drugs, several of which have been used to treat malaria infections in humans. Overall, the mevalonate bypass system provides a way to probe essential aspects of apicoplast biology and study the effects of drugs that target apicoplast processes.

Malaria parasites rely on an organelle called the apicoplast for growth and survival. Antimalarial drugs such as azithromycin inhibit basic processes in the apicoplast and result in the disruption of the organelle. Surprisingly, addition of a single metabolite, isopentenyl pyrophosphate (IPP), allows the parasites to survive in culture after disruption of the apicoplast. Unfortunately, using IPP to study this phenomenon has several limitations: IPP is prohibitively expensive, has to be used at high concentrations, and has a half-life less than 5 hours. To address these problems, we engineered parasites to express four enzymes from an alternative pathway capable of producing IPP in the parasites. We validated this new system and used it to metabolically label essential metabolites, to delete an essential apicoplast protein, and to characterize the state of apicoplast-disrupted parasites. A key finding from these studies comes from transcriptomic and metabolomic analysis of parasites treated with the drug azithromycin. We found that apicoplast disruption results in few changes in parasite metabolism. In particular, the expression of hundreds of nuclear-encoded apicoplast proteins are not affected by disruption of the apicoplast organelle, making it likely that apicoplast metabolic pathways and processes are still functional in apicoplast-disrupted parasites.

Multiscale Computational Analysis of Right Ventricular Mechanoenergetics

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.

Short-term metabolic adjustments in Plasmodium falciparum counter hypoxanthine deprivation at the expense of long-term viability

Background

The malarial parasite Plasmodium falciparum is an auxotroph for purines, which are required for nucleic acid synthesis during the intra-erythrocytic developmental cycle (IDC) of the parasite. The capabilities of the parasite and extent to which it can use compensatory mechanisms to adapt to purine deprivation were studied by examining changes in its metabolism under sub-optimal concentrations of hypoxanthine, the primary precursor utilized by the parasite for purine-based nucleic acid synthesis.

Methods

The concentration of hypoxanthine that caused a moderate growth defect over the course of one IDC was determined. At this concentration of hypoxanthine (0.5 μM), transcriptomic and metabolomic data were collected during one IDC at multiple time points. These data were integrated with a metabolic network model of the parasite embedded in a red blood cell (RBC) to interpret the metabolic adaptation of P. falciparum to hypoxanthine deprivation.

Results

At a hypoxanthine concentration of 0.5 μM, vacuole-like structures in the cytosol of many P. falciparum parasites were observed after the 24-h midpoint of the IDC. Parasites grown under these conditions experienced a slowdown in the progression of the IDC. After 72 h of deprivation, the parasite growth could not be recovered despite supplementation with 90 µM hypoxanthine. Simulations of P. falciparum metabolism suggested that alterations in ubiquinone, isoprenoid, shikimate, and mitochondrial metabolism occurred before the appearance of these vacuole-like structures. Alterations were found in metabolic reactions associated with fatty acid synthesis, the pentose phosphate pathway, methionine metabolism, and coenzyme A synthesis in the latter half of the IDC. Furthermore, gene set enrichment analysis revealed that P. falciparum activated genes associated with rosette formation, Maurer’s cleft and protein export under two different nutrient-deprivation conditions (hypoxanthine and isoleucine).

Conclusions

The metabolic network analysis presented here suggests that P. falciparum invokes specific purine-recycling pathways to compensate for hypoxanthine deprivation and maintains a hypoxanthine pool for purine-based nucleic acid synthesis. However, this compensatory mechanism is not sufficient to maintain long-term viability of the parasite. Although P. falciparum can complete a full IDC in low hypoxanthine conditions, subsequent cycles are disrupted.

Electronic supplementary material

The online version of this article (10.1186/s12936-019-2720-3) contains supplementary material, which is available to authorized users.

Metabolic host responses to malarial infection during the intraerythrocytic developmental cycle

BACKGROUND

The malarial parasite Plasmodium falciparum undergoes a complex life cycle, including an intraerythrocytic developmental cycle, during which it is metabolically dependent on the infected human red blood cell (RBC). To describe whole cell metabolic activity within both P. falciparum and RBCs during the asexual reproduction phase of the intraerythrocytic developmental cycle, we developed an integrated host-parasite metabolic modeling framework driven by time-dependent gene expression data.

RESULTS

We validated the model by reproducing the experimentally determined 1) stage-specific production of biomass components and their precursors in the parasite and 2) metabolite concentration changes in the medium of P. falciparum-infected RBC cultures. The model allowed us to explore time- and strain-dependent P. falciparum metabolism and hypothesize how host cell metabolism alters in response to malarial infection. Specifically, the metabolic analysis showed that uninfected RBCs that coexist with infected cells in the same culture decrease their production of 2,3-bisphosphoglycerate, an oxygen-carrying regulator, reducing the ability of hemoglobin in these cells to release oxygen. Furthermore, in response to parasite-induced oxidative stress, infected RBCs downgraded their glycolytic flux by using the pentose phosphate pathway and secreting ribulose-5-phosphate. This mechanism links individually observed experimental phenomena, such as glycolytic inhibition and ribulose-5-phosphate secretion, to the oxidative stress response.

CONCLUSIONS

Although the metabolic model does not incorporate regulatory mechanisms per se, alterations in gene expression levels caused by regulatory mechanisms are manifested in the model as altered metabolic states. This provides the model the capability to capture complex multicellular host-pathogen metabolic interactions of the infected RBC culture. The system-level analysis revealed complex relationships such as how the parasite can reduce oxygen release in uninfected cells in the presence of infected RBCs as well as the role of different metabolic pathways involved in the oxidative stress response of infected RBCs.

Dynamics of cross-bridge cycling, ATP hydrolysis, force generation, and deformation in cardiac muscle

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.

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

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

Markov chain Monte Carlo based analysis of post-translationally modified VDAC gating kinetics

The voltage-dependent anion channel (VDAC) is the main conduit for permeation of solutes (including nucleotides and metabolites) of up to 5 kDa across the mitochondrial outer membrane (MOM). Recent studies suggest that VDAC activity is regulated via post-translational modifications (PTMs). Yet the nature and effect of these modifications is not understood. Herein, single channel currents of wild-type, nitrosated, and phosphorylated VDAC are analyzed using a generalized continuous-time Markov chain Monte Carlo (MCMC) method. This developed method describes three distinct conducting states (open, half-open, and closed) of VDAC activity. Lipid bilayer experiments are also performed to record single VDAC activity under un-phosphorylated and phosphorylated conditions, and are analyzed using the developed stochastic search method. Experimental data show significant alteration in VDAC gating kinetics and conductance as a result of PTMs. The effect of PTMs on VDAC kinetics is captured in the parameters associated with the identified Markov model. Stationary distributions of the Markov model suggest that nitrosation of VDAC not only decreased its conductance but also significantly locked VDAC in a closed state. On the other hand, stationary distributions of the model associated with un-phosphorylated and phosphorylated VDAC suggest a reversal in channel conformation from relatively closed state to an open state. Model analyses of the nitrosated data suggest that faster reaction of nitric oxide with Cys-127 thiol group might be responsible for the biphasic effect of nitric oxide on basal VDAC conductance.

Computational analysis of Ca2+ dynamics in isolated cardiac mitochondria predicts two distinct modes of Ca2+ uptake

Cardiac mitochondria can act as a significant Ca(2+) sink and shape cytosolic Ca(2+) signals affecting various cellular processes, such as energy metabolism and excitation-contraction coupling. However, different mitochondrial Ca(2+) uptake mechanisms are still not well understood. In this study, we analysed recently published Ca(2+) uptake experiments performed on isolated guinea pig cardiac mitochondria using a computer model of mitochondrial bioenergetics and cation handling. The model analyses of the data suggest that the majority of mitochondrial Ca(2+) uptake, at physiological levels of cytosolic Ca(2+) and Mg(2+), occurs through a fast Ca(2+) uptake pathway, which is neither the Ca(2+) uniporter nor the rapid mode of Ca(2+) uptake. This fast Ca(2+) uptake component was explained by including a biophysical model of the ryanodine receptor (RyR) in the computer model. However, the Mg(2+)-dependent enhancement of the RyR adaptation was not evident in this RyR-type channel, in contrast to that of cardiac sarcoplasmic reticulum RyR. The extended computer model is corroborated by simulating an independent experimental dataset, featuring mitochondrial Ca(2+) uptake, egress and sequestration. The model analyses of the two datasets validate the existence of two classes of Ca(2+) buffers that comprise the mitochondrial Ca(2+) sequestration system. The modelling study further indicates that the Ca(2+) buffers respond differentially depending on the source of Ca(2+) uptake. In particular, it suggests that the Class 1 Ca(2+) buffering capacity is auto-regulated by the rate at which Ca(2+) is taken up by mitochondria.

Analysis of cardiovascular dynamics in pulmonary hypertensive C57BL6/J mice

A computer model was used to analyze data on cardiac and vascular mechanics from C57BL6/J mice exposed to 0 (n = 4), 14 (n = 6), 21 (n = 8) and 28 (n = 7) days of chronic hypoxia and treatment with the VEGF receptor inhibitor SUGEN (HySu) to induce pulmonary hypertension. Data on right ventricular pressure and volume, and systemic arterial pressure obtained before, during, and after inferior vena cava occlusion were analyzed using a mathematical model of realistic ventricular mechanics coupled with a simple model of the pulmonary and systemic vascular systems. The model invokes a total of 26 adjustable parameters, which were estimated based on least-squares fitting of the data. Of the 26 adjustable parameters, 14 were set to globally constant values for the entire data set. It was necessary to adjust the remaining 12 parameters to match data from all experimental groups. Of these 12 individually adjusted parameters, three parameters representing pulmonary vascular resistance, pulmonary arterial elastance, and pulmonary arterial narrowing were found to significantly change in HySu-induced remodeling. Model analysis shows a monotonic change in these parameters as disease progressed, with approximately 130% increase in pulmonary resistance, 70% decrease in unstressed pulmonary arterial volume, and 110% increase in pulmonary arterial elastance in the 28-day group compared to the control group. These changes are consistent with prior experimental measurements. Furthermore, the 28-day data could be explained only after increasing the passive elastance of the right free wall compared to the value used for the other data sets, which is likely a consequence of the increased RV collagen accumulation found experimentally. These findings may indicate a compensatory remodeling followed by pathological remodeling of the right ventricle in HySu-induced pulmonary hypertension.

A biophysical model of the mitochondrial ATP-Mg/P(i) carrier

Mitochondrial adenine nucleotide (AdN) content is regulated through the Ca(2+)-activated, electroneutral ATP-Mg/P(i) carrier (APC). The APC is a protein in the mitochondrial carrier super family that localizes to the inner mitochondrial membrane (IMM). It is known to modulate a number of processes that depend on mitochondrial AdN content, such as gluconeogenesis, protein synthesis, and citrulline synthesis. Despite this critical role, a kinetic model of the underlying mechanism has not been developed and validated. Here, a biophysical model of the APC is developed that is thermodynamically balanced and accurately reproduces a number of reported data sets from isolated rat liver and rat kidney mitochondria. The model is based on an ordered bi-bi mechanism for heteroexchange of ATP and P(i) and includes homoexchanges of ATP and P(i) to explain both the initial rate and time course data on ATP and P(i) transport via the APC. The model invokes seven kinetic parameters regarding the APC mechanism and three parameters related to matrix pH regulation by external P(i). These parameters are estimated based on 19 independent data curves; the estimated parameters are validated using six additional data curves. The model takes into account the effects of pH, Mg(2+), and Ca(2+) on ATP and P(i) transport via the APC, and supports the conclusion that the pH gradient across the IMM serves as the primary driving force for AdN uptake or efflux. Moreover, computer simulations demonstrate that extramatrix Ca(2+) modulates the turnover rate of the APC and not the binding affinity of ATP, as previously suggested.

A mathematical model of the tripartite synapse: astrocyte-induced synaptic plasticity

In this paper, we present a biologically detailed mathematical model of tripartite synapses, where astrocytes modulate short-term synaptic plasticity. The model consists of a pre-synaptic bouton, a post-synaptic dendritic spine-head, a synaptic cleft and a peri-synaptic astrocyte controlling Ca(2 + ) dynamics inside the synaptic bouton. This in turn controls glutamate release dynamics in the cleft. As a consequence of this, glutamate concentration in the cleft has been modeled, in which glutamate reuptake by astrocytes has also been incorporated. Finally, dendritic spine-head dynamics has been modeled. As an application, this model clearly shows synaptic potentiation in the hippocampal region, i.e., astrocyte Ca(2 + ) mediates synaptic plasticity, which is in conformity with the majority of the recent findings (Perea and Araque (Science 317, 1083-1086, 2007); Henneberger et al. (Nature 463, 232-236, 2010); Navarrete et al. (PLoS Biol. 10, e1001259, 2012)).

Bisalbuminemia in nephrotic syndrome (a case report)

A rare case of bisalbuminemia in a patient with minimal lesion glomerulonephritis (MLGN) in which the course of the disease and therapeutic response to steroid were typical of minimal lesion glomerulonephritis, is presented and literature reviewed.

Serum amylase and its isoenzymes in cases of viral hepatitis

Serum amylase (SA) and its isoenzymes were studied in 9 control cases and 26 cases of Viral Hepatitis (VH). A significant rise in SA level was found in VH cases (Mean value--287.15 + 100.64) as compared to control cases (Mean value--155.56 + 28.57). The value was markedly increased in 8 cases with hepatic encephalopathy (Mean value--388 + 57.21). Isoenzymes of SA were estimated in all the cases by means of polyacrylamide gel electrophoresis and a direct saccharogenic assay of amylase activity. Besides two regular peaks (Pancreatic and Salivary fractions) a third peak (Slowest of all) was noticed in 4 out of 9 control cases (44.5%) and 19 out of 26 cases of VH (73.1%). Out of 8 cases who developed hepatic encephalopathy during the course, 7 had this peak (87.5%). SA levels and its isoenzymes level correlated well with the extent of hepatic damage (SGPT and Serum bilirubin levels). The possibility of hepatic origin of the third peak has been discussed.