Screen for Chemical Modulators of Autophagy Reveals Novel Therapeutic Inhibitors of mTORC1 Signaling

Background

Mammalian target of rapamycin complex 1 (mTORC1) is a protein kinase that relays nutrient availability signals to control numerous cellular functions including autophagy, a process of cellular self-eating activated by nutrient depletion. Addressing the therapeutic potential of modulating mTORC1 signaling and autophagy in human disease requires active chemicals with pharmacologically desirable properties.

Methodology/Principal Findings

Using an automated cell-based assay, we screened a collection of >3,500 chemicals and identified three approved drugs (perhexiline, niclosamide, amiodarone) and one pharmacological reagent (rottlerin) capable of rapidly increasing autophagosome content. Biochemical assays showed that the four compounds stimulate autophagy and inhibit mTORC1 signaling in cells maintained in nutrient-rich conditions. The compounds did not inhibit mTORC2, which also contains mTOR as a catalytic subunit, suggesting that they do not inhibit mTOR catalytic activity but rather inhibit signaling to mTORC1. mTORC1 inhibition and autophagosome accumulation induced by perhexiline, niclosamide or rottlerin were rapidly reversed upon drug withdrawal whereas amiodarone inhibited mTORC1 essentially irreversibly. TSC2, a negative regulator of mTORC1, was required for inhibition of mTORC1 signaling by rottlerin but not for mTORC1 inhibition by perhexiline, niclosamide and amiodarone. Transient exposure of immortalized mouse embryo fibroblasts to these drugs was not toxic in nutrient-rich conditions but led to rapid cell death by apoptosis in starvation conditions, by a mechanism determined in large part by the tuberous sclerosis complex protein TSC2, an upstream regulator of mTORC1. By contrast, transient exposure to the mTORC1 inhibitor rapamycin caused essentially irreversible mTORC1 inhibition, sustained inhibition of cell growth and no selective cell killing in starvation.

Conclusion/Significance

The observation that drugs already approved for human use can reversibly inhibit mTORC1 and stimulate autophagy should greatly facilitate the preclinical and clinical testing of mTORC1 inhibition for indications such as tuberous sclerosis, diabetes, cardiovascular disease and cancer.     

Comparative Analysis of RNA Families Reveals Distinct Repertoires for Each Domain of Life

The RNA world hypothesis, that RNA genomes and catalysts preceded DNA genomes and genetically-encoded protein catalysts, has been central to models for the early evolution of life on Earth. A key part of such models is continuity between the earliest stages in the evolution of life and the RNA repertoires of extant lineages. Some assessments seem consistent with a diverse RNA world, yet direct continuity between modern RNAs and an RNA world has not been demonstrated for the majority of RNA families, and, anecdotally, many RNA functions appear restricted in their distribution. Despite much discussion of the possible antiquity of RNA families, no systematic analyses of RNA family distribution have been performed. To chart the broad evolutionary history of known RNA families, we performed comparative genomic analysis of over 3 million RNA annotations spanning 1446 families from the Rfam 10 database. We report that 99% of known RNA families are restricted to a single domain of life, revealing discrete repertoires for each domain. For the 1% of RNA families/clans present in more than one domain, over half show evidence of horizontal gene transfer (HGT), and the rest show a vertical trace, indicating the presence of a complex protein synthesis machinery in the Last Universal Common Ancestor (LUCA) and consistent with the evolutionary history of the most ancient protein-coding genes. However, with limited interdomain transfer and few RNA families exhibiting demonstrable antiquity as predicted under RNA world continuity, our results indicate that the majority of modern cellular RNA repertoires have primarily evolved in a domain-specific manner.     

Global Analysis of ATM Polymorphism Reveals Significant Functional Constraint

ATM, the gene that is mutated in ataxia-telangiectasia, is associated with cerebellar degeneration, abnormal proliferation of small blood vessels, and cancer. These clinically important manifestations have stimulated interest in defining the sequence variation in the ATM gene. Therefore, we undertook a comprehensive survey of sequence variation in ATM in diverse human populations. The protein-encoding exons of the gene (9,168 bp) and the adjacent intron and untranslated sequences (14,661 bp) were analyzed in 93 individuals from seven major human populations. In addition, the coding sequence was analyzed in one chimpanzee, one gorilla, one orangutan, and one Old World monkey. In human ATM, 88 variant sites were discovered by denaturing high-performance liquid chromatography, which is 96%-100% sensitive for detection of DNA sequence variation. ATM was compared to 14 other autosomal genes for nucleotide diversity. The noncoding regions of ATM had diversity values comparable to other genes, but the coding regions had very low diversity, especially in the last 29% of the protein sequence. A test of the neutral evolution hypothesis, through use of the Hudson/Kreitman/Aguadé statistic, revealed that this region of the human ATM gene was significantly constrained relative to that of the orangutan, the Old World monkey, and the mouse, but not relative to that of the chimpanzee or the gorilla. ATM displayed extensive linkage disequilibrium, consistent with suppression of meiotic recombination at this locus. Seven haplotypes were defined. Two haplotypes accounted for 82% of all chromosomes analyzed in all major populations; two others carrying the same D126E missense polymorphism accounted for 33% of chromosomes in Africa but were never observed outside of Africa. The high frequency of this polymorphism may be due either to a population expansion within Africa or to selective pressure.     

Engineering of a Polydisperse Small Heat-Shock Protein Reveals Conserved Motifs of Oligomer Plasticity

1. Download high-res image (186KB)
2. Download full-size image

     

Meta-Analysis of Arabidopsis thaliana Phospho-Proteomics Data Reveals Compartmentalization of Phosphorylation Motifs

Protein (de)phosphorylation plays an important role in plants. To provide a robust foundation for subcellular phosphorylation signaling network analysis and kinase-substrate relationships, we performed a meta-analysis of 27 published and unpublished in-house mass spectrometry–based phospho-proteome data sets for Arabidopsis thaliana covering a range of processes, (non)photosynthetic tissue types, and cell cultures. This resulted in an assembly of 60,366 phospho-peptides matching to 8141 nonredundant proteins. Filtering the data for quality and consistency generated a set of medium and a set of high confidence phospho-proteins and their assigned phospho-sites. The relation between single and multiphosphorylated peptides is discussed. The distribution of p-proteins across cellular functions and subcellular compartments was determined and showed overrepresentation of protein kinases. Extensive differences in frequency of pY were found between individual studies due to proteomics and mass spectrometry workflows. Interestingly, pY was underrepresented in peroxisomes but overrepresented in mitochondria. Using motif-finding algorithms motif-x and MMFPh at high stringency, we identified compartmentalization of phosphorylation motifs likely reflecting localized kinase activity. The filtering of the data assembly improved signal/noise ratio for such motifs. Identified motifs were linked to kinases through (bioinformatic) enrichment analysis. This study also provides insight into the challenges/pitfalls of using large-scale phospho-proteomic data sets to nonexperts.     

Network-Based Elucidation of Human Disease Similarities Reveals Common Functional Modules Enriched for Pluripotent Drug Targets

Current work in elucidating relationships between diseases has largely been based on pre-existing knowledge of disease genes. Consequently, these studies are limited in their discovery of new and unknown disease relationships. We present the first quantitative framework to compare and contrast diseases by an integrated analysis of disease-related mRNA expression data and the human protein interaction network. We identified 4,620 functional modules in the human protein network and provided a quantitative metric to record their responses in 54 diseases leading to 138 significant similarities between diseases. Fourteen of the significant disease correlations also shared common drugs, supporting the hypothesis that similar diseases can be treated by the same drugs, allowing us to make predictions for new uses of existing drugs. Finally, we also identified 59 modules that were dysregulated in at least half of the diseases, representing a common disease-state “signature”. These modules were significantly enriched for genes that are known to be drug targets. Interestingly, drugs known to target these genes/proteins are already known to treat significantly more diseases than drugs targeting other genes/proteins, highlighting the importance of these core modules as prime therapeutic opportunities.     

Structural and Inhibition Analysis Reveals the Mechanism of Selectivity of a Series of Aggrecanase Inhibitors

Several inhibitors of a series of cis-1(S)2(R)-amino-2-indanol-based compounds were reported to be selective for the aggrecanases, ADAMTS-4 and -5 over other metalloproteases. To understand the nature of this selectivity for aggrecanases, the inhibitors, along with the broad spectrum metalloprotease inhibitor marimastat, were independently bound to the catalytic domain of ADAMTS-5, and the corresponding crystal structures were determined. By comparing the structures, it was determined that the specificity of the relative inhibitors for ADAMTS-5 was not driven by a specific interaction, such as zinc chelation, hydrogen bonding, or charge interactions, but rather by subtle and indirect factors, such as water bridging, ring rigidity, pocket size, and shape, as well as protein conformation flexibility.Osteoarthritis (OA)³ pathology includes degradation of articular cartilage, along with subchondral bone sclerosis and osteophyte formation, all contributing to impaired joint function. Pain, restricted movement, and joint instability accompany these structural changes and often result in the need for total joint replacement. Current therapies alleviate the mild to moderate pain and inflammation associated with OA, but do not protect the cartilage from further damage and have not demonstrated an effect on disease progression (1). Therefore, therapeutics that prevent or slow the alteration of joint structure and function will address a major unmet medical need.Loss of aggrecan, a macromolecular proteoglycan providing cartilage with its properties of compressibility and resilience, is a major phenotype associated with OA and is believed to be a critical event in driving disease progression (2, 3). Both ex vivo and in vivo proof of concept studies support ADAMTS-4 and ADAMTS-5, commonly referred to as aggrecanase-1 and -2, respectively, as the two major enzymes responsible for the proteolytic breakdown of cartilage aggrecan (reviewed in Ref. 4). Blocking their activity may be an attractive strategy to stop or slow down the progression of the disease, as suggested by studies in knock-out mice (5). Given the chronic nature of the disease, long term treatment will be likely, demanding very safe therapeutic interventions only achievable with ADAMTS-4- and ADAMTS-5-specific inhibitors lacking off-target side effects. Designing selectivity has been very challenging and a major source of difficulty is that at least 57 metalloproteases (MP) divided in three major families, 1) matrix metalloproteases (MMP); 2) a disintegrin and metalloproteinase (ADAM); and 3) a disintegrin and metalloproteinase with thrombospondin motifs (ADAMTS), are present in humans. ADAMTS-4 and -5 belong to the ADAMTS family and share common catalytic and structural features with the other MP members. These features include the highly conserved amino acid sequence, HEXXHXXGXXH, harboring a catalytic zinc cation, required for activation of the peptide bond toward hydrolysis. In addition, many MPs share significant structural topology in the active site, such as a flexible S1′ loop. To complicate matters further, only a handful of MP structures, mostly in the MMP family, have been determined, and the functions of most MPs still remain unknown, a fact that has earned them an orphan status denomination. Lack of structural information has hindered the design for inhibitor specificity and without known substrates for many MPs, assays for screening inhibitors are often not available, making determinations of selectivity difficult (6, 7).Yao et al. (8, 9) reported the discovery of a series of (2R)-N4-hydroxy-2-(3-hydroxybenzyl)-N1-[(1S,2R)-2-hydroxy-,3-dihydro-1H-inden-1-yl]butanediamide derivatives as potent and selective inhibitors of aggrecanase activity. Using a homology model of aggrecanase based on the active site of atrolysin C and adamalysin II and docking compound 8, shown in Fig. 1, the authors concluded that the 3-hydroxyl group of inhibitor 8 achieved selectivity through a specific hydrogen-bonding interaction with Thr⁴⁴⁰ (numerical numbering is based on the human sequence of ADAMTS-5) in the S1′ pocket of aggrecanase.Open in a separate windowFIGURE 1.Aggrecanase inhibitors. Structures of the inhibitors evaluated in these studies. The P1 moiety of each inhibitor is marked in red.Whereas both ADAMTS-4 and -5 have a threonine at this position, MMP-1, -2, -3, -7, -8, -9, -10, -13, -14, -16, and ADAM-17 have a valine, lending credence to the proposed hypothesis around selectivity. Recently, our group has established a protocol for crystallizing the catalytic domain of ADAMTS-5 and determined its three-dimensional structure (10). Thus, experimental validation or invalidation of the hypothesis that inhibitor 8 and related molecules form a specific hydrogen bond with the hydroxyl group of threonine in the S1′ pocket of aggrecanases could now be determined by protease-inhibitor crystallographic analysis.In the current study, we wanted to confirm and extend the selectivity profile of compound 8 and 11 against a wide array of MPs. Moreover, we wanted to elucidate the key molecular interactions responsible for the enhanced selectivity profile of this series of compounds. For this purpose we generated co-crystals and solved the structures of marimastat, compound 8, and compound 11 in complex with the catalytic domain of recombinant human ADAMTS-5.     

HGDP and HapMap Analysis by Ancestry Mapper Reveals Local and Global Population Relationships

Knowledge of human origins, migrations, and expansions is greatly enhanced by the availability of large datasets of genetic information from different populations and by the development of bioinformatic tools used to analyze the data. We present Ancestry Mapper, which we believe improves on existing methods, for the assignment of genetic ancestry to an individual and to study the relationships between local and global populations. The principle function of the method, named Ancestry Mapper, is to give each individual analyzed a genetic identifier, made up of just 51 genetic coordinates, that corresponds to its relationship to the HGDP reference population. As a consequence, the Ancestry Mapper Id (AMid) has intrinsic biological meaning and provides a tool to measure similarity between world populations. We applied Ancestry Mapper to a dataset comprised of the HGDP and HapMap data. The results show distinctions at the continental level, while simultaneously giving details at the population level. We clustered AMids of HGDP/HapMap and observe a recapitulation of human migrations: for a small number of clusters, individuals are grouped according to continental origins; for a larger number of clusters, regional and population distinctions are evident. Calculating distances between AMids allows us to infer ancestry. The number of coordinates is expandable, increasing the power of Ancestry Mapper. An R package called Ancestry Mapper is available to apply this method to any high density genomic data set.     

Global Gene Expression Analysis Reveals a Link between NDRG1 and Vesicle Transport

N-myc downstream-regulated gene 1 (NDRG1) is induced by cellular stress such as hypoxia and DNA damage, and in humans, germ line mutations cause Charcot-Marie-Tooth disease. However, the cellular roles of NDRG1 are not fully understood. Previously, NDRG1 was shown to mediate doxorubicin resistance under hypoxia, suggesting a role for NDRG1 in cell survival under these conditions. We found decreased apoptosis in doxorubicin-treated cells expressing NDRG1 shRNAs under normoxia, demonstrating a requirement for NDRG1 in apoptosis in breast epithelial cells under normal oxygen pressure. Also, different cellular stress regimens, such as hypoxia and doxorubicin treatment, induced NDRG1 through different stress signalling pathways. We further compared expression profiles in human breast epithelial cells ectopically over-expressing NDRG1 with cells expressing NDRG1 shRNAs in order to identify biological pathways where NDRG1 is involved. The results suggest that NDRG1 may have roles connected to vesicle transport.     

Genome Sequencing of Autism-Affected Families Reveals Disruption of Putative Noncoding Regulatory DNA

We performed whole-genome sequencing (WGS) of 208 genomes from 53 families affected by simplex autism. For the majority of these families, no copy-number variant (CNV) or candidate de novo gene-disruptive single-nucleotide variant (SNV) had been detected by microarray or whole-exome sequencing (WES). We integrated multiple CNV and SNV analyses and extensive experimental validation to identify additional candidate mutations in eight families. We report that compared to control individuals, probands showed a significant (p = 0.03) enrichment of de novo and private disruptive mutations within fetal CNS DNase I hypersensitive sites (i.e., putative regulatory regions). This effect was only observed within 50 kb of genes that have been previously associated with autism risk, including genes where dosage sensitivity has already been established by recurrent disruptive de novo protein-coding mutations (ARID1B, SCN2A, NR3C2, PRKCA, and DSCAM). In addition, we provide evidence of gene-disruptive CNVs (in DISC1, WNT7A, RBFOX1, and MBD5), as well as smaller de novo CNVs and exon-specific SNVs missed by exome sequencing in neurodevelopmental genes (e.g., CANX, SAE1, and PIK3CA). Our results suggest that the detection of smaller, often multiple CNVs affecting putative regulatory elements might help explain additional risk of simplex autism.     

The Proteome of the Isolated Chlamydia trachomatis Containing Vacuole Reveals a Complex Trafficking Platform Enriched for Retromer Components

Chlamydia trachomatis is an important human pathogen that replicates inside the infected host cell in a unique vacuole, the inclusion. The formation of this intracellular bacterial niche is essential for productive Chlamydia infections. Despite its importance for Chlamydia biology, a holistic view on the protein composition of the inclusion, including its membrane, is currently missing. Here we describe the host cell-derived proteome of isolated C. trachomatis inclusions by quantitative proteomics. Computational analysis indicated that the inclusion is a complex intracellular trafficking platform that interacts with host cells’ antero- and retrograde trafficking pathways. Furthermore, the inclusion is highly enriched for sorting nexins of the SNX-BAR retromer, a complex essential for retrograde trafficking. Functional studies showed that in particular, SNX5 controls the C. trachomatis infection and that retrograde trafficking is essential for infectious progeny formation. In summary, these findings suggest that C. trachomatis hijacks retrograde pathways for effective infection.