Resolving medulloblastoma cellular architecture by single-cell genomics.

Medulloblastoma is a malignant childhood cerebellar tumour type that comprises distinct molecular subgroups. Whereas genomic characteristics of these subgroups are well defined, the extent to which cellular diversity underlies their divergent biology and clinical behaviour remains largely unexplored. Here we used single-cell transcriptomics to investigate intra- and intertumoral heterogeneity in 25 medulloblastomas spanning all molecular subgroups. WNT, SHH and Group 3 tumours comprised subgroup-specific undifferentiated and differentiated neuronal-like malignant populations, whereas Group 4 tumours consisted exclusively of differentiated neuronal-like neoplastic cells. SHH tumours closely resembled granule neurons of varying differentiation states that correlated with patient age. Group 3 and Group 4 tumours exhibited a developmental trajectory from primitive progenitor-like to more mature neuronal-like cells, the relative proportions of which distinguished these subgroups. Cross-species transcriptomics defined distinct glutamatergic populations as putative cells-of-origin for SHH and Group 4 subtypes. Collectively, these data provide insights into the cellular and developmental states underlying subtype-specific medulloblastoma biology.

SH2 Domains Serve as Lipid-Binding Modules for pTyr-Signaling Proteins.

The Src-homology 2 (SH2) domain is a protein interaction domain that directs myriad phosphotyrosine (pY)-signaling pathways. Genome-wide screening of human SH2 domains reveals that ∼90% of SH2 domains bind plasma membrane lipids and many have high phosphoinositide specificity. They bind lipids using surface cationic patches separate from pY-binding pockets, thus binding lipids and the pY motif independently. The patches form grooves for specific lipid headgroup recognition or flat surfaces for non-specific membrane binding and both types of interaction are important for cellular function and regulation of SH2 domain-containing proteins. Cellular studies with ZAP70 showed that multiple lipids bind its C-terminal SH2 domain in a spatiotemporally specific manner and thereby exert exquisite spatiotemporal control over its protein binding and signaling activities in T cells. Collectively, this study reveals how lipids control SH2 domain-mediated cellular protein-protein interaction networks and suggest a new strategy for therapeutic modulation of pY-signaling pathways.

Genome-wide functional annotation of dual-specificity protein- and lipid-binding modules that regulate protein interactions.

Emerging evidence indicates that membrane lipids regulate protein networking by directly interacting with protein-interaction domains (PIDs). As a pilot study to identify and functionally annodate lipid-binding PIDs on a genomic scale, we performed experimental and computational studies of PDZ domains. Characterization of 70 PDZ domains showed that ~40% had submicromolar membrane affinity. Using a computational model built from these data, we predicted the membrane-binding properties of 2,000 PDZ domains from 20 species. The accuracy of the prediction was experimentally validated for 26 PDZ domains. We also subdivided lipid-binding PDZ domains into three classes based on the interplay between membrane- and protein-binding sites. For different classes of PDZ domains, lipid binding regulates their protein interactions by different mechanisms. Functional studies of a PDZ domain protein, rhophilin 2, suggest that all classes of lipid-binding PDZ domains serve as genuine dual-specificity modules regulating protein interactions at the membrane under physiological conditions.

Lipids Regulate Lck Protein Activity through Their Interactions with the Lck Src Homology 2 Domain.

Lymphocyte-specific protein-tyrosine kinase (Lck) plays an essential role in T cell receptor (TCR) signaling and T cell development, but its activation mechanism is not fully understood. To explore the possibility that plasma membrane (PM) lipids control TCR signaling activities of Lck, we measured the membrane binding properties of its regulatory Src homology 2 (SH2) and Src homology 3 domains. The Lck SH2 domain binds anionic PM lipids with high affinity but with low specificity. Electrostatic potential calculation, NMR analysis, and mutational studies identified the lipid-binding site of the Lck SH2 domain that includes surface-exposed basic, aromatic, and hydrophobic residues but not the phospho-Tyr binding pocket. Mutation of lipid binding residues greatly reduced the interaction of Lck with the ζ chain in the activated TCR signaling complex and its overall TCR signaling activities. These results suggest that PM lipids, including phosphatidylinositol 4,5-bisphosphate and phosphatidylinositol 3,4,5-trisphosphate, modulate interaction of Lck with its binding partners in the TCR signaling complex and its TCR signaling activities in a spatiotemporally specific manner via its SH2 domain.

p21-activated Kinases (PAKs) Mediate the Phosphorylation of PREX2 Protein to Initiate Feedback Inhibition of Rac1 GTPase.

Phosphatidylinositol 3,4,5-trisphosphate (PIP3)-dependent Rac exchanger 2 (PREX2) is a guanine nucleotide exchange factor (GEF) for the Ras-related C3 botulinum toxin substrate 1 (Rac1) GTPase, facilitating the exchange of GDP for GTP on Rac1. GTP-bound Rac1 then activates its downstream effectors, including p21-activated kinases (PAKs). PREX2 and Rac1 are frequently mutated in cancer and have key roles within the insulin-signaling pathway. Rac1 can be inactivated by multiple mechanisms; however, negative regulation by insulin is not well understood. Here, we show that in response to being activated after insulin stimulation, Rac1 initiates its own inactivation by decreasing PREX2 GEF activity. Following PREX2-mediated activation of Rac1 by the second messengers PIP3 or Gßγ, we found that PREX2 was phosphorylated through a PAK-dependent mechanism. PAK-mediated phosphorylation of PREX2 reduced GEF activity toward Rac1 by inhibiting PREX2 binding to PIP3 and Gßγ. Cell fractionation experiments also revealed that phosphorylation prevented PREX2 from localizing to the cellular membrane. Furthermore, the onset of insulin-induced phosphorylation of PREX2 was delayed compared with AKT. Altogether, we propose that second messengers activate the Rac1 signal, which sets in motion a cascade whereby PAKs phosphorylate and negatively regulate PREX2 to decrease Rac1 activation. This type of regulation would allow for transient activation of the PREX2-Rac1 signal and may be relevant in multiple physiological processes, including diseases such as diabetes and cancer when insulin signaling is chronically activated.

Genome-wide structural analysis reveals novel membrane binding properties of AP180 N-terminal homology (ANTH) domains.

An increasing number of cytosolic proteins are shown to interact with membrane lipids during diverse cellular processes, but computational prediction of these proteins and their membrane binding behaviors remains challenging. Here, we introduce a new combinatorial computation protocol for systematic and robust functional prediction of membrane-binding proteins through high throughput homology modeling and in-depth calculation of biophysical properties. The approach was applied to the genomic scale identification of the AP180 N-terminal homology (ANTH) domain, one of the modular lipid binding domains, and prediction of their membrane binding properties. Our analysis yielded comprehensive coverage of the ANTH domain family and allowed classification and functional annotation of proteins based on the differences in local structural and biophysical features. Our analysis also identified a group of plant ANTH domains with unique structural features that may confer novel functionalities. Experimental characterization of a representative member of this subfamily confirmed its unique membrane binding mechanism and unprecedented membrane deforming activity. Collectively, these studies suggest that our new computational approach can be applied to genome-wide functional prediction of other lipid binding domains.

Antigen cross-presentation in young tumor-bearing hosts promotes CD8+ T cell terminal differentiation.

The immune system undergoes a progressive functional remodeling with age. Understanding how age bias shapes antitumor immunity is essential in designing effective immunotherapies, especially for pediatric patients. Here, we explore antitumor CD8+ T cell responses generated in young (prepubescent) and adult (presenescent) mice. Using an MHCI-deficient tumor model, we observed that tumor-reactive CD8+ T cells expanded in young tumor-bearing (TB) mice acquired a terminally differentiated phenotype characterized by overexpression of inhibitory receptors and the transcription factor Tox1. Furthermore, tumor-infiltrating CD8+ T cells from young tumors yielded a poor cytokine response compared with CD8+ T cells infiltrating adult tumors. Young migratory dendritic cells (migDCs) from the draining lymph nodes (dLNs), and mononuclear phagocytic cells (MPCs) infiltrating young tumors, were more competent in capturing and cross-presenting tumor antigen, leading to enhanced priming of CD8+ T cells in dLNs and their subsequent terminal differentiation in the tumors. Single-cell transcriptional profiling of tumor-infiltrating MPCs demonstrated that young MPCs are polarized toward an inflammatory, effector phenotype. Consistent with our observations in young versus adult TB mice, analysis of immune infiltrates from pediatric solid tumors showed a correlation between tumor-infiltrating CD8+ T cells with an exhaustion phenotype and the frequency of PD-L1-expressing monocytes/macrophages. Collectively, these data indicate that a young tissue microenvironment contributes to the generation of an immune response skewed toward a less pliable terminal effector state, thus narrowing the window for immunotherapeutic interventions.

Genomes for Kids: The Scope of Pathogenic Mutations in Pediatric Cancer Revealed by Comprehensive DNA and RNA Sequencing.

Genomic studies of pediatric cancer have primarily focused on specific tumor types or high-risk disease. Here, we used a three-platform sequencing approach, including whole-genome sequencing (WGS), whole-exome sequencing (WES), and RNA sequencing (RNA-seq), to examine tumor and germline genomes from 309 prospectively identified children with newly diagnosed (85%) or relapsed/refractory (15%) cancers, unselected for tumor type. Eighty-six percent of patients harbored diagnostic (53%), prognostic (57%), therapeutically relevant (25%), and/or cancer-predisposing (18%) variants. Inclusion of WGS enabled detection of activating gene fusions and enhancer hijacks (36% and 8% of tumors, respectively), small intragenic deletions (15% of tumors), and mutational signatures revealing of pathogenic variant effects. Evaluation of paired tumor-normal data revealed relevance to tumor development for 55% of pathogenic germline variants. This study demonstrates the power of a three-platform approach that incorporates WGS to interrogate and interpret the full range of genomic variants across newly diagnosed as well as relapsed/refractory pediatric cancers. SIGNIFICANCE: Pediatric cancers are driven by diverse genomic lesions, and sequencing has proven useful in evaluating high-risk and relapsed/refractory cases. We show that combined WGS, WES, and RNA-seq of tumor and paired normal tissues enables identification and characterization of genetic drivers across the full spectrum of pediatric cancers. This article is highlighted in the In This Issue feature, p. 2945.

High-throughput computational structure-based characterization of protein families: START domains and implications for structural genomics.

SkyLine, a high-throughput homology modeling pipeline tool, detects and models true sequence homologs to a given protein structure. Structures and models are stored in SkyBase with links to computational function annotation, as calculated by MarkUs. The SkyLine/SkyBase/MarkUs technology represents a novel structure-based approach that is more objective and versatile than other protein classification resources. This structure-centric strategy provides a multi-dimensional organization and coverage of protein space at the levels of family, function, and genome. The concept of "modelability", the ability to model sequences on related structures, provides a reliable criterion for membership in a protein family ("leverage") and underlies the unique success of this approach. The overall procedure is illustrated by its application to START domains, which comprise a Biomedical Theme for the Northeast Structural Genomics Consortium as part of the Protein Structure Initiative. START domains are typically involved in the non-vesicular transport of lipids. While 19 experimentally determined structures are available, the family, whose evolutionary hierarchy is not well determined, is highly sequence diverse, and the ligand-binding potential of many family members is unknown. The SkyLine/SkyBase/MarkUs approach provides significant insights and predicts: (1) many more family members (approximately 4,000) than any other resource; (2) the function for a large number of unannotated proteins; (3) instances of START domains in genomes from which they were thought to be absent; and (4) the existence of two types of novel proteins, those containing dual START domain and those containing N-terminal START domains.

Clinical genome sequencing uncovers potentially targetable truncations and fusions of MAP3K8 in spitzoid and other melanomas.

Spitzoid melanoma is a specific morphologic variant of melanoma that most commonly affects children and adolescents, and ranges on the spectrum of malignancy from low grade to overtly malignant. These tumors are generally driven by fusions of ALK, RET, NTRK1/3, MET, ROS1 and BRAF1,2. However, in approximately 50% of cases no genetic driver has been established2. Clinical whole-genome and transcriptome sequencing (RNA-Seq) of a spitzoid tumor from an adolescent revealed a novel gene fusion of MAP3K8, encoding a serine-threonine kinase that activates MEK3,4. The patient, who had exhausted all other therapeutic options, was treated with a MEK inhibitor and underwent a transient clinical response. We subsequently analyzed spitzoid tumors from 49 patients by RNA-Seq and found in-frame fusions or C-terminal truncations of MAP3K8 in 33% of cases. The fusion transcripts and truncated genes all contained MAP3K8 exons 1-8 but lacked the autoinhibitory final exon. Data mining of RNA-Seq from the Cancer Genome Atlas (TCGA) uncovered analogous MAP3K8 rearrangements in 1.5% of adult melanomas. Thus, MAP3K8 rearrangements-uncovered by comprehensive clinical sequencing of a single case-are the most common genetic event in spitzoid melanoma, are present in adult melanomas and could be amenable to MEK inhibition.

Structural alignment of protein--DNA interfaces: insights into the determinants of binding specificity.

A new method is introduced to structurally align interfaces observed in protein--DNA complexes. The method is based on a procedure that describes the interfacial geometry in terms of the spatial relationships between individual amino acid--nucleotide pairs. An amino acid--amino acid similarity matrix, S, is defined that provides a quantitative measure of the geometric relationships of amino acids in different interfaces and the entire stretch of "local" DNA within some distance of each amino acid. S is used as a substitution matrix in a dynamic programming algorithm that aligns the interfacial amino acids of the two complexes. The quality of the alignment is determined by an interface alignment score, IAS, that provides a quantitative measure of the similarity in the docking geometry between two protein--DNA complexes. We have clustered a large set of protein--DNA complexes based on their IAS values. In general, proteins within a single family form identifiable clusters. Subgroup clustering is often observed within families offering a fine-grained description of docking geometries. Although proteins with similar folds tend to dock in similar ways, important differences are observed even for structural motifs that almost perfectly align. Relationships are observed between the interfaces formed in cognate and non-cognate complexes involving the same proteins indicating a strong driving force to maintain certain contacts, even if this requires a distortion of the DNA. There are cases where inter-family similarities are greater than intra-family similarities. Our method offers the possibility of comparing different protein--DNA interfaces in a detailed, objective and quantitative fashion. This offers the possibility of new approaches to the description of the determinants of molecular recognition and to the prediction of protein and DNA sequence combinations that are optimal for binding.

Template-based prediction of protein function.

We discuss recent approaches for structure-based protein function annotation. We focus on template-based methods where the function of a query protein is deduced from that of a template for which both the structure and function are known. We describe the different ways of identifying a template. These are typically based on sequence analysis but new methods based on purely structural similarity are also being developed that allow function annotation based on structural relationships that cannot be recognized by sequence. The growing number of available structures of known function, improved homology modeling techniques and new developments in the use of structure allow template-based methods to be applied on a proteome-wide scale and in many different biological contexts. This progress significantly expands the range of applicability of structural information in function annotation to a level that previously was only achievable by sequence comparison.

An ankyrin repeat domain of AKR2 drives chloroplast targeting through coincident binding of two chloroplast lipids.

In organellogenesis of the chloroplast from endosymbiotic cyanobacteria, the establishment of protein-targeting mechanisms to the chloroplast should have been pivotal. However, it is still mysterious how these mechanisms were established and how they work in plant cells. Here we show that AKR2A, the cytosolic targeting factor for chloroplast outer membrane (COM) proteins, evolved from the ankyrin repeat domain (ARD) of the host cell by stepwise extensions of its N-terminal domain and that two lipids, monogalactosyldiacylglycerol (MGDG) and phosphatidylglycerol (PG), of the endosymbiont were selected to function as the AKR2A receptor. Structural analysis, molecular modeling, and mutational analysis of the ARD identified two adjacent sites for coincidental and synergistic binding of MGDG and PG. Based on these findings, we propose that the targeting mechanism of COM proteins was established using components from both the endosymbiont and host cell through a modification of the protein-protein-interacting ARD into a lipid binding domain.

Enhanced apoptosis of B and T lymphocytes in TAFII105 dominant-negative transgenic mice is linked to nuclear factor-kappa B.

The general transcription factor TFIID is composed of the TATA-binding protein (TBP) and 12-14 TBP-associated factors (TAF(II)s). Some TAF(II)s act as bridges between transcription activators and the general transcription machinery through direct interaction with activation domains. Although TAF-mediated transcription activation has been established, there is little genetic evidence connecting it to binding of an activator. TAF(II)105 is a substoichiometric subunit of transcription factor IID highly expressed in B lymphocytes. In this study, we examined the physiological role of TAF(II)105 and its mechanism of action in vivo by expressing two forms of dominant-negative mutant TAF(II)105 in mice. We show that TAF(II)105 has a pro-survival role in B and T lymphocytes, where the native protein is expressed. In addition, TAF(II)105 is important for T cell maturation and for production of certain antibody isotypes. These phenotypic alterations were absent in mice expressing a dominant-negative mutant that lacks one of the domains mediating p65/RelA binding in vitro. These findings provide support to the notion that interaction between the activator and TAF is important for their function in vivo.