Quantitative assessment of cell fate decision between autophagy and apoptosis

Abstract Autophagy and apoptosis are cellular processes that regulate cell survival and death, the former by eliminating dysfunctional components in the cell, the latter by programmed cell death. Stress signals can induce either process, and it is unclear how cells ‘assess’ cellular damage and make a ‘life’ or ‘death’ decision upon activating autophagy or apoptosis. A computational model of coupled apoptosis and autophagy is built here to analyze the underlying signaling and regulatory network dynamics. The model explains the experimentally observed differential deployment of autophagy and apoptosis in response to various stress signals. Autophagic response dominates at low-to-moderate stress; whereas the response shifts from autophagy (graded activation) to apoptosis (switch-like activation) with increasing stress intensity. The model reveals that cytoplasmic Ca2+ acts as a rheostat that fine-tunes autophagic and apoptotic responses. A G-protein signaling-mediated feedback loop maintains cytoplasmic Ca2+ level, which in turn governs autophagic response through an AMP-activated protein kinase (AMPK)-mediated feedforward loop. Ca2+/calmodulin-dependent kinase kinase β (CaMKKβ) emerges as a determinant of the competing roles of cytoplasmic Ca2+ in autophagy regulation. The study demonstrates that the proposed model can be advantageously used for interrogating cell regulation events and developing pharmacological strategies for modulating cell decisions

Mechanisms of action of autophagy modulators dissected by quantitative systems pharmacology analysis

Autophagy plays an essential role in cell survival/death and functioning. Modulation of autophagy has been recognized as a promising therapeutic strategy against diseases/disorders associated with uncontrolled growth or accumulation of biomolecular aggregates, organelles, or cells including those caused by cancer, aging, neurodegeneration, and liver diseases such as α1-antitrypsin deficiency. Numerous pharmacological agents that enhance or suppress autophagy have been discovered. However, their molecular mechanisms of action are far from clear. Here, we collected a set of 225 autophagy modulators and carried out a comprehensive quantitative systems pharmacology (QSP) analysis of their targets using both existing databases and predictions made by our machine learning algorithm. Autophagy modulators include several highly promiscuous drugs (e.g., artenimol and olanzapine acting as activators, fostamatinib as an inhibitor, or melatonin as a dual-modulator) as well as selected drugs that uniquely target specific proteins (~30% of modulators). They are mediated by three layers of regulation: (i) pathways involving core autophagy-related (ATG) proteins such as mTOR, AKT, and AMPK; (ii) upstream signaling events that regulate the activity of ATG pathways such as calcium-, cAMP-, and MAPK-signaling pathways; and (iii) transcription factors regulating the expression of ATG proteins such as TFEB, TFE3, HIF-1, FoxO, and NF-κB. Our results suggest that PKA serves as a linker, bridging various signal transduction events and autophagy. These new insights contribute to a better assessment of the mechanism of action of autophagy modulators as well as their side effects, development of novel polypharmacological strategies, and identification of drug repurposing opportunities

Simulating Radiating and Magnetized Flows in Multi-Dimensions with ZEUS-MP

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.Comment: Accepted for publication in the ApJ Supplement. 42 pages with 29 inlined figures; uses emulateapj.sty. Discussions in sections 2 - 4 improved per referee comments; several figures modified to illustrate grid resolution. ZEUS-MP source code and documentation available from the Laboratory for Computational Astrophysics at http://lca.ucsd.edu/codes/currentcodes/zeusmp2

An analog of glibenclamide selectively enhances autophagic degradation of misfolded α1-antitrypsin Z

The classical form of α1-antitrypsin deficiency (ATD) is characterized by intracellular accumulation of the misfolded variant α1-antitrypsin Z (ATZ) and severe liver disease in some of the affected individuals. In this study, we investigated the possibility of discovering novel therapeutic agents that would reduce ATZ accumulation by interrogating a C. elegans model of ATD with high-content genome-wide RNAi screening and computational systems pharmacology strategies. The RNAi screening was utilized to identify genes that modify the intracellular accumulation of ATZ and a novel computational pipeline was developed to make high confidence predictions on repurposable drugs. This approach identified glibenclamide (GLB), a sulfonylurea drug that has been used broadly in clinical medicine as an oral hypoglycemic agent. Here we show that GLB promotes autophagic degradation of misfolded ATZ in mammalian cell line models of ATD. Furthermore, an analog of GLB reduces hepatic ATZ accumulation and hepatic fibrosis in a mouse model in vivo without affecting blood glucose or insulin levels. These results provide support for a drug discovery strategy using simple organisms as human disease models combined with genetic and computational screening methods. They also show that GLB and/or at least one of its analogs can be immediately tested to arrest the progression of human ATD liver disease.</div

CemOrange2 fusions facilitate multifluorophore subcellular imaging in C. elegans

Due to its ease of genetic manipulation and transparency, Caenorhabditis elegans (C. elegans) has become a preferred model system to study gene function by microscopy. The use of Aequorea victoria green fluorescent protein (GFP) fused to proteins or targeting sequences of interest, further expanded upon the utility of C. elegans by labeling subcellular structures, which enables following their disposition during development or in the presence of genetic mutations. Fluorescent proteins with excitation and emission spectra different from that of GFP accelerated the use of multifluorophore imaging in real time. We have expanded the repertoire of fluorescent proteins for use in C. elegans by developing a codon-optimized version of Orange2 (CemOrange2). Proteins or targeting motifs fused to CemOrange2 were distinguishable from the more common fluorophores used in the nematode; such as GFP, YFP, and mKate2. We generated a panel of CemOrange2 fusion constructs, and confirmed they were targeted to their correct subcellular addresses by colocalization with independent markers. To demonstrate the potential usefulness of this new panel of fluorescent protein markers, we showed that CemOrange2 fusion proteins could be used to: 1) monitor biological pathways, 2) multiplex with other fluorescent proteins to determine colocalization and 3) gain phenotypic knowledge of a human ABCA3 orthologue, ABT-4, trafficking variant in the C. elegans model organism

Dominant negative variants in KIF5B cause osteogenesis imperfecta via down regulation of mTOR signaling

BACKGROUND: Kinesin motor proteins transport intracellular cargo, including mRNA, proteins, and organelles. Pathogenic variants in kinesin-related genes have been implicated in neurodevelopmental disorders and skeletal dysplasias. We identified de novo, heterozygous variants in KIF5B, encoding a kinesin-1 subunit, in four individuals with osteogenesis imperfecta. The variants cluster within the highly conserved kinesin motor domain and are predicted to interfere with nucleotide binding, although the mechanistic consequences on cell signaling and function are unknown. METHODS: To understand the in vivo genetic mechanism of KIF5B variants, we modeled the p.Thr87Ile variant that was found in two patients in the C. elegans ortholog, unc-116, at the corresponding position (Thr90Ile) by CRISPR/Cas9 editing and performed functional analysis. Next, we studied the cellular and molecular consequences of the recurrent p.Thr87Ile variant by microscopy, RNA and protein analysis in NIH3T3 cells, primary human fibroblasts and bone biopsy. RESULTS: C. elegans heterozygous for the unc-116 Thr90Ile variant displayed abnormal body length and motility phenotypes that were suppressed by additional copies of the wild type allele, consistent with a dominant negative mechanism. Time-lapse imaging of GFP-tagged mitochondria showed defective mitochondria transport in unc-116 Thr90Ile neurons providing strong evidence for disrupted kinesin motor function. Microscopy studies in human cells showed dilated endoplasmic reticulum, multiple intracellular vacuoles, and abnormal distribution of the Golgi complex, supporting an intracellular trafficking defect. RNA sequencing, proteomic analysis, and bone immunohistochemistry demonstrated down regulation of the mTOR signaling pathway that was partially rescued with leucine supplementation in patient cells. CONCLUSION: We report dominant negative variants in the KIF5B kinesin motor domain in individuals with osteogenesis imperfecta. This study expands the spectrum of kinesin-related disorders and identifies dysregulated signaling targets for KIF5B in skeletal development

: A Model System for Anti-Cancer Drug Discovery and Therapeutic Target Identification

The nematode Caenorhabditis elegans (C. elegans) offers a unique opportunity for biological and basic medical researches due to its genetic tractability and well-defined developmental lineage. It also provides an exceptional model for genetic, molecular, and cellular analysis of human disease-related genes. Recently, C. elegans has been used as an ideal model for the identification and functional analysis of drugs (or small-molecules) in vivo. In this review, we describe conserved oncogenic signaling pathways (Wnt, Notch, and Ras) and their potential roles in the development of cancer stem cells. During C. elegans germline development, these signaling pathways regulate multiple cellular processes such as germline stem cell niche specification, germline stem cell maintenance, and germ cell fate specification. Therefore, the aberrant regulations of these signaling pathways can cause either loss of germline stem cells or overproliferation of a specific cell type, resulting in sterility. This sterility phenotype allows us to identify drugs that can modulate the oncogenic signaling pathways directly or indirectly through a high-throughput screening. Current in vivo or in vitro screening methods are largely focused on the specific core signaling components. However, this phenotype-based screening will identify drugs that possibly target upstream or downstream of core signaling pathways as well as exclude toxic effects. Although phenotype-based drug screening is ideal, the identification of drug targets is a major challenge. We here introduce a new technique, called Drug Affinity Responsive Target Stability (DARTS). This innovative method is able to identify the target of the identified drug. Importantly, signaling pathways and their regulators in C. elegans are highly conserved in most vertebrates, including humans. Therefore, C. elegans will provide a great opportunity to identify therapeutic drugs and their targets, as well as to understand mechanisms underlying the formation of cancer

Macrocephaly and developmental delay caused by missense variants in RAB5C

Rab GTPases are important regulators of intracellular vesicular trafficking. RAB5C is a member of the Rab GTPase family that plays an important role in the endocytic pathway, membrane protein recycling and signaling. Here we report on 12 individuals with nine different heterozygous de novo variants in RAB5C. All but one patient with missense variants (n = 9) exhibited macrocephaly, combined with mild-to-moderate developmental delay. Patients with loss of function variants (n = 2) had an apparently more severe clinical phenotype with refractory epilepsy and intellectual disability but a normal head circumference. Four missense variants were investigated experimentally. In vitro biochemical studies revealed that all four variants were damaging, resulting in increased nucleotide exchange rate, attenuated responsivity to guanine exchange factors and heterogeneous effects on interactions with effector proteins. Studies in C. elegans confirmed that all four variants were damaging in vivo and showed defects in endocytic pathway function. The variant heterozygotes displayed phenotypes that were not observed in null heterozygotes, with two shown to be through a dominant negative mechanism. Expression of the human RAB5C variants in zebrafish embryos resulted in defective development, further underscoring the damaging effects of the RAB5C variants. Our combined bioinformatic, in vitro and in vivo experimental studies and clinical data support the association of RAB5C missense variants with a neurodevelopmental disorder characterized by macrocephaly and mild-to-moderate developmental delay through disruption of the endocytic pathway