Analyzing the Functional Properties of the Creatine Kinase System with Multiscale 'Sloppy' Modeling

In this study the function of the two isoforms of creatine kinase (CK; EC 2.7.3.2) in myocardium is investigated. The ‘phosphocreatine shuttle’ hypothesis states that mitochondrial and cytosolic CK plays a pivotal role in the transport of high-energy phosphate (HEP) groups from mitochondria to myofibrils in contracting muscle. Temporal buffering of changes in ATP and ADP is another potential role of CK. With a mathematical model, we analyzed energy transport and damping of high peaks of ATP hydrolysis during the cardiac cycle. The analysis was based on multiscale data measured at the level of isolated enzymes, isolated mitochondria and on dynamic response times of oxidative phosphorylation measured at the whole heart level. Using ‘sloppy modeling’ ensemble simulations, we derived confidence intervals for predictions of the contributions by phosphocreatine (PCr) and ATP to the transfer of HEP from mitochondria to sites of ATP hydrolysis. Our calculations indicate that only 15±8% (mean±SD) of transcytosolic energy transport is carried by PCr, contradicting the PCr shuttle hypothesis. We also predicted temporal buffering capabilities of the CK isoforms protecting against high peaks of ATP hydrolysis (3750 µM*s(−1)) in myofibrils. CK inhibition by 98% in silico leads to an increase in amplitude of mitochondrial ATP synthesis pulsation from 215±23 to 566±31 µM*s(−1), while amplitudes of oscillations in cytosolic ADP concentration double from 77±11 to 146±1 µM. Our findings indicate that CK acts as a large bandwidth high-capacity temporal energy buffer maintaining cellular ATP homeostasis and reducing oscillations in mitochondrial metabolism. However, the contribution of CK to the transport of high-energy phosphate groups appears limited. Mitochondrial CK activity lowers cytosolic inorganic phosphate levels while cytosolic CK has the opposite effect

Norovirus Genetic Diversity – from within patient viral evolution to global distribution

Noroviruses belong to the family of Caliciviridae and cause acute gastroenteritis. The genetic diversity within the genus Norovirus is extremely large and novel genotypes, recombinants within and between genotypes, and antigenic drift variants are regularly discovered. The distribution and incidence of norovirus genotypes changes over time and by geographic region, with GII.4 as the predominant genotype responsible for approximately 70% of the outbreaks detected worldwide. This thesis contributes to a better understanding of the global norovirus diversity, and the role of chronic norovirus infection among immunocompromised individuals on virus diversity, antigenic variation and evolution. We describe the emergence of a novel GII.4 drift variant and novel GII.17 in chapter 2.1 and 2.2, respectively. Chapter 3 provides an overview of the global norovirus molecular epidemiology using data obtained by the NoroNet surveillance network in the last ten years. Chapter 4 shows that GII.4 was not predominant before the wide spread use of norovirus molecular typing methods in the mid-1990s. Chapter 5 shows that chronic norovirus infection is common among solid organ transplant patients and more than doubled the number of described patients in literature. Finally, chapter 6 shows that virus populations in the immunocompromised host are genetically distinct from viruses circulating in the general population. These hosts therefore may contain a reservoir of newly emerging strains and future studies need to address whether these new strains can be transmitted to other immunocompromised patients or individuals in the general population. Norovirus vaccines are currently tested in clinical trials. The observed global changes in the norovirus molecular epidemiology require a norovirus vaccine with broad immune protection and an antigenic component that can be easily updated within a limited time span. Immunocompromised patients with chronic norovirus infection may suffer of persistent symptoms of diarrhoea with a consequent severe loss of quality of life and future studies need to address how these patients can be treated to eliminate infection

FluxSimulator: An R package to simulate isotopomer distributions in metabolic networks

The representation of biochemical knowledge in terms of fluxes (transformation rates) in a metabolic network is often a crucial step in the development of new drugs and efficient bioreactors. Mass spectroscopy (MS) and nuclear magnetic resonance spectroscopy (NMRS) in combination with 13C labeled substrates are experimental techniques resulting in data that may be used to quantify fluxes in the metabolic network underlying a process. The massive amount of data generated by spectroscopic experiments increasingly requires software which models the dynamics of the underlying biological system. In this work we present an approach to handle isotopomer distributions in metabolic networks using an object-oriented programming approach, implemented using S4 classes in R. The developed package is called FluxSimulator and provides a user friendly interface to specify the topological information of the metabolic network as well as carbon atom transitions in plain text files. The package automatically derives the mathematical representation of the formulated network, and assembles a set of ordinary differential equations (ODEs) describing the change of each isotopomer pool over time. These ODEs are subsequently solved numerically. In a case study FluxSimulator was applied to an example network. Our results indicate that the package is able to reproduce exact changes in isotopomer compositions of the metabolite pools over time at given flux rates

Coordinated behavior of mitochondria in both space and time: a reactive species-activated wave of mitochondrial depolarization.

Reactive oxygen species (ROS) can trigger a transient burst of mitochondrial ROS production via ROS activation of the mitochondrial permeability transition pore (MPTP), a phenomenon termed ROS-induced ROS release (RIRR). The goal of this study was to investigate if the generation of ROS in a discrete region of a cardiomyocyte could serve to propagate RIRR-mediated mitochondrial depolarizations throughout a cell. Our experiments revealed that localized RIRR activated either RIRR-mediated fluctuations in mitochondrial membrane potential (time period: 3–10 min) or a traveling wave of depolarization of the cell's mitochondria (velocity: ∼5 μm/min). Both phenomena appeared to be mediated by the mitochondrial permeability transition pore and eventually encompassed the majority of the mitochondrial population of both isolated rat and rabbit cardiomyocytes. Furthermore, depolarization was often reversible; the waves of depolarization were then followed by a rapid (∼40 μm/min) repolarization wave of the mitochondria. We show that the RIRR can function to communicate the mitochondrial permeability transition from one mitochondrion to another in the isolated adult cardiomyocyte

Computational estimation of tricarboxylic acid cycle fluxes using noisy NMR data from cardiac biopsies

Background: The aerobic energy metabolism of cardiac muscle cells is of major importance for the contractile function of the heart. Because energy metabolism is very heterogeneously distributed in heart tissue, especially during coronary disease, a method to quantify metabolic fluxes in small tissue samples is desirable. Taking tissue biopsies after infusion of substrates labeled with stable carbon isotopes makes this possible in animal experiments. However, the appreciable noise level in NMR spectra of extracted tissue samples makes computational estimation of metabolic fluxes challenging and a good method to define confidence regions was not yet available.Results: Here we present a computational analysis method for nuclear magnetic resonance (NMR) measurements of tricarboxylic acid (TCA) cycle metabolites. The method was validated using measurements on extracts of single tissue biopsies taken from porcine heart in vivo. Isotopic enrichment of glutamate was measured by NMR spectroscopy in tissue samples taken at a single time point after the timed infusion of 13C labeled substrates for the TCA cycle. The NMR intensities for glutamate were analyzed with a computational model describing carbon transitions in the TCA cycle and carbon exchange with amino acids. The model dynamics depended on five flux parameters, which were optimized to fit the NMR measurements. To determine confidence regions for the estimated fluxes, we used the Metropolis-Hastings algorithm for Markov chain Monte Carlo (MCMC) sampling to generate extensive ensembles of feasible flux combinations that describe the data within meas

SARS-CoV-2-Specific Antibody Detection for Seroepidemiology: A Multiplex Analysis Approach Accounting for Accurate Seroprevalence

Background. The COVID-19 pandemic necessitates better understanding of the kinetics of antibody production induced by infection with SARS-CoV-2. We aimed to develop a high-throughput multiplex assay to detect antibodies to SARS-CoV-2 to assess immunity to the virus in the general population. Methods. Spike protein subunits S1 and receptor binding domain, and nucleoprotein were coupled to microspheres. Sera collected before emergence of SARS-CoV-2 (n = 224) and of non-SARS-CoV-2 influenza-like illness (n = 184), and laboratory-confirmed cases of SARS-CoV-2 infection (n = 115) with various severities of COVID-19 were tested for SARS-CoV-2–specific IgG concentrations. Results. Our assay discriminated SARS-CoV-2–induced antibodies and those induced by other viruses. The assay specificity was 95.1%–99.0% with sensitivity 83.6%–95.7%. By merging the test results for all 3 antigens a specificity of 100% was achieved with a sensitivity of at least 90%. Hospitalized COVID-19 patients developed higher IgG concentrations and the rate of IgG production

Comparison of norovirus genogroup I, II and IV seroprevalence among children in the Netherlands, 1963, 1983 and 2006

Noroviruses are a major cause of acute gastroenteritis worldwide and are a genetically diverse group of viruses. Since 2002, an increasing number of norovirus outbreaks have been reported globally, but it is not clear whether this increase has been caused by a higher awareness or reflects the emergence of new genogroup II genotype 4 (GII.4) variants. The hypothesis that norovirus prevalence has increased post-2002 and is related to the emergence of GII.4 is tested in this study. Sera collected from children aged <5 years of three Dutch cross-sectional population based cohorts in 1963, 1983 and 2006/2007 (n=143, n=130 and n=376, respectively) were tested for specific serum IgG by protein array using antigens to GII.4 and a range of other antigens representing norovirus GI, GII and GIV genotypes. The protein array was validated by paired sera of norovirus infected patients and supernatants of B-cell cultures with single epitope specificity. Evidence for norovirus infection was found to be common among Dutch children in each cohort, but the prevalence towards different genotypes changed over time. At the genogroup level, GI seroprevalence decreased significantly between 1963 and 2006/2007, while a significant increase of GII and, in particular, specific antibodies of the genotype GII.4 was detected in the 2006/2007 cohort. There were no children with only GII.4 antibodies in the 1963 cohort. This study shows that the high GII.4 norovirus incidence in very young children is a recent phenomenon. These findings are of importance for vaccine development and trials that are currently focusing mostly on GII.4 viruses

Emergence of a novel GII.17 norovirus – end of the GII.4 era?

In the winter of 2014/15 a novel GII.P17-GII.17 norovirus strain (GII.17 Kawasaki 2014) emerged, as a major cause of gastroenteritis outbreaks in China and Japan. Since their emergence these novel GII.P17-GII.17 viruses have replaced the previously dominant GII.4 genotype Sydney 2012 variant in some areas in Asia but were only detected in a limited number of cases on other continents. This perspective provides an overview of the available information on GII.17 viruses in order to gain insight in the viral and host characteristics of this norovirus genotype. We further discuss the emergence of this novel GII.P17-GII.17 norovirus in context of current knowledge on the epidemiology of noroviruses. It remains to be seen if the currently dominant norovirus strain GII.4 Sydney 2012 will be replaced in other parts of the world. Nevertheless, the public health community and surveillance systems need to be prepared in case of a potential increase of norovirus activity in the next seasons caused by this novel GII.P17-GII.17 norovirus

Serological evidence of influenza a viruses in frugivorous bats from Africa

Bats are likely natural hosts for a range of zoonotic viruses such as Marburg, Ebola, Rabies, as well as for various Corona- and Paramyxoviruses. In 2009/10, researchers discovered RNA of two novel influenza virus subtypes - H17N10 and H18N11 - in Central and South American fruit bats. The identification of bats as possible additional reservoir for influenza A viruses raises questions about the role of this mammalian taxon in influenza A virus ecology and possible public health relevance. As molecular testing can be limited by a short time window in which the virus is present, serological testing provides information about past infections and virus spread in populations after the virus has been cleared. This study aimed at screening available sera from 100 free-ranging, frugivorous bats (Eidolon helvum) sampled in 2009/10 in Ghana, for the presence of antibodies against the complete panel of influenza A haemagglutinin (HA) types ranging from H1 to H18 by means of a protein microarray platform. This technique enables simultaneous serological testing against multiple recombinant HA-types in 5μl of serum. Preliminary results indicate serological evidence against avian influenza subtype H9 in about 30% of the animals screened, with low-level cross-reactivity to phylogenetically closely related subtypes H8 and H12. To our knowledge, this is the first report of serological evidence of influenza A viruses other than H17 and H18 in bats. As avian influenza subtype H9 is associated with human infections, the implications of our findings from

New viruses in idiopathic human diarrhea cases, the Netherlands

Emerging viral infections can be identified by using a viral metagenomics approach for clinical human material. Diarrhea samples of patients with unexplained gastroenteritis from the Netherlands were analyzed by using viral metagenomics. Novel circular DNA viruses, bufaviruses, and genogroup III picobirnaviruses were identified. These data expand our knowledge of the human virome