Relating simulation studies by provenance—Developing a family of Wnt signaling models

For many biological systems, a variety of simulation models exist. A new simulation model is rarely developed from scratch, but rather revises and extends an existing one. A key challenge, however, is to decide which model might be an appropriate starting point for a particular problem and why. To answer this question, we need to identify entities and activities that contributed to the development of a simulation model. Therefore, we exploit the provenance data model, PROV-DM, of the World Wide Web Consortium and, building on previous work, continue developing a PROV ontology for simulation studies. Based on a case study of 19 Wnt/β-catenin signaling models, we identify crucial entities and activities as well as useful metadata to both capture the provenance information from individual simulation studies and relate these forming a family of models. The approach is implemented in WebProv, a web application for inserting and querying provenance information. Our specialization of PROV-DM contains the entities Research Question, Assumption, Requirement, Qualitative Model, Simulation Model, Simulation Experiment, Simulation Data, and Wet-lab Data as well as activities referring to building, calibrating, validating, and analyzing a simulation model. We show that most Wnt simulation models are connected to other Wnt models by using (parts of) these models. However, the overlap, especially regarding the Wet-lab Data used for calibration or validation of the models is small. Making these aspects of developing a model explicit and queryable is an important step for assessing and reusing simulation models more effectively. Exposing this information helps to integrate a new simulation model within a family of existing ones and may lead to the development of more robust and valid simulation models. We hope that our approach becomes part of a standardization effort and that modelers adopt the benefits of provenance when considering or creating simulation models.     

α1-AMP-Activated Protein Kinase Regulates Hypoxia-Induced Na,K-ATPase Endocytosis via Direct Phosphorylation of Protein Kinase Cζ

Hypoxia promotes Na,K-ATPase endocytosis via protein kinase Cζ (PKCζ)-mediated phosphorylation of the Na,K-ATPase α subunit. Here, we report that hypoxia leads to the phosphorylation of 5′-AMP-activated protein kinase (AMPK) at Thr172 in rat alveolar epithelial cells. The overexpression of a dominant-negative AMPK α subunit (AMPK-DN) construct prevented the hypoxia-induced endocytosis of Na,K-ATPase. The overexpression of the reactive oxygen species (ROS) scavenger catalase prevented hypoxia-induced AMPK activation. Moreover, hypoxia failed to activate AMPK in mitochondrion-deficient ρ⁰-A549 cells, suggesting that mitochondrial ROS play an essential role in hypoxia-induced AMPK activation. Hypoxia-induced PKCζ translocation to the plasma membrane and phosphorylation at Thr410 were prevented by the pharmacological inhibition of AMPK or by the overexpression of the AMPK-DN construct. We found that AMPK α phosphorylates PKCζ on residue Thr410 within the PKCζ activation loop. Importantly, the activation of AMPK α was necessary for hypoxia-induced AMPK-PKCζ binding in alveolar epithelial cells. The overexpression of T410A mutant PKCζ prevented hypoxia-induced Na,K-ATPase endocytosis, confirming that PKCζ Thr410 phosphorylation is essential for this process. PKCζ activation by AMPK is isoform specific, as small interfering RNA targeting the α1 but not the α2 catalytic subunit prevented PKCζ activation. Accordingly, we provide the first evidence that hypoxia-generated mitochondrial ROS lead to the activation of the AMPK α1 isoform, which binds and directly phosphorylates PKCζ at Thr410, thereby promoting Na,K-ATPase endocytosis.When exposed to low oxygen levels (hypoxia), cells develop adaptative strategies to maintain adequate levels of ATP (21). These strategies include increasing the efficiency of energy-producing pathways, mostly through anaerobic glycolysis, while decreasing energy-consuming processes such as Na,K-ATPase activity (30). Alveolar hypoxia occurs in many respiratory disorders, and it has been shown to decrease epithelial active Na⁺ transport, leading to impaired fluid reabsorption (37, 41, 42). Active Na⁺ transport and, thus, alveolar fluid reabsortion are effected mostly via apical sodium channels and the basolateral Na,K-ATPase (32, 38, 42). We have reported previously that hypoxia inhibits Na,K-ATPase activity by promoting its endocytosis from the plasma membrane by a mechanism that requires the generation of mitochondrial reactive oxygen species (ROS) and the phosphorylation of the Na,K-ATPase α subunit at Ser18 by protein kinase Cζ (PKCζ) (8, 9).The 5′-AMP-activated protein kinase (AMPK) is a heterotrimeric Ser/Thr kinase composed of a catalytic α subunit and regulatory β and γ subunits. Both isoforms of the AMPK catalytic subunit (α1 and α2) form complexes with noncatalytic subunits. The α1 subunit is ubiquitously expressed, whereas the α2 subunit isoform is expressed predominantly in tissues like the liver, heart, and skeletal muscle (36). The α1 and α2 subunit isoforms have ∼90% homology in their N-terminal catalytic domains and ∼60% homology in their C-terminal domains (36), suggesting that they may have distinct downstream targets (31). AMPK activation requires phosphorylation at Thr172 in the activation loop of the α subunit by upstream kinases (12, 19). Findings from recent studies suggest that AMPK is an important signaling intermediary in coupling ion transport and metabolism (15). Indeed, it has been reported that the pharmacological activation of AMPK inhibits amiloride- and ouabain-sensitive epithelial Na⁺ transport (15). Moreover, the activities of the epithelial Na⁺ channel (ENaC) (2, 17), the Na,K-ATPase (40), and the cystic fibrosis transmembrane conductance regulator (17) have been shown to be inhibited by AMPK. Here, we provide evidence that hypoxia, via mitochondrial ROS, leads to AMPK activation and that AMPK binds to and directly phosphorylates PKCζ in an isoform-specific manner, thus promoting Na,K-ATPase endocytosis in alveolar epithelial cells (AEC).     

The mutational structure of metabolism in Caenorhabditis elegans

A properly functioning organism must maintain metabolic homeostasis. Deleterious mutations degrade organismal function, presumably at least in part via effects on metabolic function. Here we present an initial investigation into the mutational structure of the Caenorhabditis elegans metabolome by means of a mutation accumulation experiment. We find that pool sizes of 29 metabolites vary greatly in their vulnerability to mutation, both in terms of the rate of accumulation of genetic variance (the mutational variance, VM) and the rate of change of the trait mean (the mutational bias, ΔM). Strikingly, some metabolites are much more vulnerable to mutation than any other trait previously studied in the same way. Although we cannot statistically assess the strength of mutational correlations between individual metabolites, principal component analysis provides strong evidence that some metabolite pools are genetically correlated, but also that there is substantial scope for independent evolution of different groups of metabolites. Averaged over mutation accumulation lines, PC3 is positively correlated with relative fitness, but a model in which metabolites are uncorrelated with fitness is nearly as good by Akaike's Information Criterion.     

Microtubule-associated Proteins 1 (MAP1) Promote Human Immunodeficiency Virus Type I (HIV-1) Intracytoplasmic Routing to the Nucleus

After cell entry, HIV undergoes rapid transport toward the nucleus using microtubules and microfilaments. Neither the cellular cytoplasmic components nor the viral proteins that interact to mediate transport have yet been identified. Using a yeast two-hybrid screen, we identified four cytoskeletal components as putative interaction partners for HIV-1 p24 capsid protein: MAP1A, MAP1S, CKAP1, and WIRE. Depletion of MAP1A/MAP1S in indicator cell lines and primary human macrophages led to a profound reduction in HIV-1 infectivity as a result of impaired retrograde trafficking, demonstrated by a characteristic accumulation of capsids away from the nuclear membrane, and an overall defect in nuclear import. MAP1A/MAP1S did not impact microtubule network integrity or cell morphology but contributed to microtubule stabilization, which was shown previously to facilitate infection. In addition, we found that MAP1 proteins interact with HIV-1 cores both in vitro and in infected cells and that interaction involves MAP1 light chain LC2. Depletion of MAP1 proteins reduced the association of HIV-1 capsids with both dynamic and stable microtubules, suggesting that MAP1 proteins help tether incoming viral capsids to the microtubular network, thus promoting cytoplasmic trafficking. This work shows for the first time that following entry into target cells, HIV-1 interacts with the cytoskeleton via its p24 capsid protein. Moreover, our results support a role for MAP1 proteins in promoting efficient retrograde trafficking of HIV-1 by stimulating the formation of stable microtubules and mediating the association of HIV-1 cores with microtubules.     

Allelic Variation of Cytochrome P450s Drives Resistance to Bednet Insecticides in a Major Malaria Vector

Scale up of Long Lasting Insecticide Nets (LLINs) has massively contributed to reduce malaria mortality across Africa. However, resistance to pyrethroid insecticides in malaria vectors threatens its continued effectiveness. Deciphering the detailed molecular basis of such resistance and designing diagnostic tools is critical to implement suitable resistance management strategies. Here, we demonstrated that allelic variation in two cytochrome P450 genes is the most important driver of pyrethroid resistance in the major African malaria vector Anopheles funestus and detected key mutations controlling this resistance. An Africa-wide polymorphism analysis of the duplicated genes CYP6P9a and CYP6P9b revealed that both genes are directionally selected with alleles segregating according to resistance phenotypes. Modelling and docking simulations predicted that resistant alleles were better metabolizers of pyrethroids than susceptible alleles. Metabolism assays performed with recombinant enzymes of various alleles confirmed that alleles from resistant mosquitoes had significantly higher activities toward pyrethroids. Additionally, transgenic expression in Drosophila showed that flies expressing resistant alleles of both genes were significantly more resistant to pyrethroids compared with those expressing the susceptible alleles, indicating that allelic variation is the key resistance mechanism. Furthermore, site-directed mutagenesis and functional analyses demonstrated that three amino acid changes (Val¹⁰⁹Ile, Asp³³⁵Glu and Asn³⁸⁴Ser) from the resistant allele of CYP6P9b were key pyrethroid resistance mutations inducing high metabolic efficiency. The detection of these first DNA markers of metabolic resistance to pyrethroids allows the design of DNA-based diagnostic tools to detect and track resistance associated with bednets scale up, which will improve the design of evidence-based resistance management strategies.     

Fast Approximate Quadratic Programming for Graph Matching

Quadratic assignment problems arise in a wide variety of domains, spanning operations research, graph theory, computer vision, and neuroscience, to name a few. The graph matching problem is a special case of the quadratic assignment problem, and graph matching is increasingly important as graph-valued data is becoming more prominent. With the aim of efficiently and accurately matching the large graphs common in big data, we present our graph matching algorithm, the Fast Approximate Quadratic assignment algorithm. We empirically demonstrate that our algorithm is faster and achieves a lower objective value on over 80% of the QAPLIB benchmark library, compared with the previous state-of-the-art. Applying our algorithm to our motivating example, matching C. elegans connectomes (brain-graphs), we find that it efficiently achieves performance.     

Comparison between the effects of diallyl tetrasulfide on human retina pigment epithelial cells (ARPE-19) and HCT116 cells

Background

Diallyl mono- and polysulfanes from garlic are known to induce an adaptive cell response and the formation of antioxidants in cancer cells. In the case of a severe ER stress and a failure in the response, cancer cells eventually go into apoptosis. Only little is known about the response of normal cells upon treatment.

Methods

Normal ARPE-19 cells were treated with diallyl tetrasulfide to study their cellular response and the results were compared with those of HCT116 cancer cells. Cell viability was checked by an MTT assay and cytofluorimetry. The formation of superoxide radicals, H₂O₂ and thiols were determined and proteins involved in the ER stress response were also detected by Western blot analysis.

Results

We found that diallyl tetrasulfide induced reactive oxygen species (ROS) in normal cells similar to cancer cells in a time (0 to 60 min) and dose dependent manner (0 to 50 μM). The level of heme oxigenase-1 (HO-1) was up-regulated in both cell types. Initially, we found a decrease in the total thiol level in both cell types but in contrast to cancer cells, normal cells recovered from the decrease in the total thiol concentration within 60 min of treatment.

Conclusions

The recovery of the thiol concentration in normal cells treated with diallyl tetrasulfide seems to be responsible for the failure to induce the ER stress signalling pathway and finally apoptosis in normal cells.

General Significance

The difference in the recovery of the thiol status might be an explanation for the anti-carcinogenic effects of garlic compounds.