Fasting glucose and insulin are intermediate traits for type 2 diabetes. Here we explore the role of coding variation on these traits by analysis of variants on the HumanExome BeadChip in 60,564 non-diabetic individuals and in 16,491 T2D cases and 81,877 controls. We identify a novel association of a low-frequency nonsynonymous SNV in GLP1R (A316T; rs10305492; MAF=1.4%) with lower FG (β=−0.09±0.01 mmol l⁻¹, P=3.4 × 10⁻¹²), T2D risk (OR[95%CI]=0.86[0.76–0.96], P=0.010), early insulin secretion (β=−0.07±0.035 pmol_insulin mmol_glucose⁻¹, P=0.048), but higher 2-h glucose (β=0.16±0.05 mmol l⁻¹, P=4.3 × 10⁻⁴). We identify a gene-based association with FG at G6PC2 (p_SKAT=6.8 × 10⁻⁶) driven by four rare protein-coding SNVs (H177Y, Y207S, R283X and S324P). We identify rs651007 (MAF=20%) in the first intron of ABO at the putative promoter of an antisense lncRNA, associating with higher FG (β=0.02±0.004 mmol l⁻¹, P=1.3 × 10⁻⁸). Our approach identifies novel coding variant associations and extends the allelic spectrum of variation underlying diabetes-related quantitative traits and T2D susceptibility.

The mesoderm- and epithelial-mesenchymal transition-associated transcription factor FOXC1 is specifically overexpressed in basal-like breast cancer (BLBC), but its biochemical function is not understood. Here, we demonstrate that FOXC1 controls cancer stem cell (CSC) properties enriched in BLBC cells via activation of Smoothened (SMO)-independent Hedgehog (Hh) signaling. This non-canonical activation of Hh is specifically mediated by Gli2. Furthermore, we show that the N-terminal domain of FOXC1 (aa 1-68) binds directly to an internal region (aa 898-1168) of Gli2, enhancing the DNA-binding and transcription-activating capacity of Gli2. FOXC1 expression correlates with that of Gli2 and its targets in human breast cancers. Moreover, FOXC1 overexpression reduces sensitivity to anti-Hedgehog (Hh) inhibitors in BLBC cells and xenograft tumors. Together, these findings reveal FOXC1-mediated non-canonical Hh signaling that determines the BLBC stem-like phenotype and anti-Hh sensitivity, supporting inhibition of FOXC1 pathways as potential approaches for improving BLBC treatment.

Glycosaminoglycan polysaccharides play critical roles in many cellular processes, ranging from viral invasion and angiogenesis to spinal cord injury. Their diverse biological activities are derived from an ability to regulate a remarkable number of proteins. However, few methods exist for the rapid identification of glycosaminoglycan–protein interactions and for studying the potential of glycosaminoglycans to assemble multimeric protein complexes. Here, we report a multidisciplinary approach that combines new carbohydrate microarray and computational modeling methodologies to elucidate glycosaminoglycan–protein interactions. The approach was validated through the study of known protein partners for heparan and chondroitin sulfate, including fibroblast growth factor 2 (FGF2) and its receptor FGFR1, the malarial protein VAR2CSA, and tumor necrosis factor-α (TNF-α). We also applied the approach to identify previously undescribed interactions between a specific sulfated epitope on chondroitin sulfate, CS-E, and the neurotrophins, a critical family of growth factors involved in the development, maintenance, and survival of the vertebrate nervous system. Our studies show for the first time that CS is capable of assembling multimeric signaling complexes and modulating neurotrophin signaling pathways. In addition, we identify a contiguous CS-E-binding site by computational modeling that suggests a potential mechanism to explain how CS may promote neurotrophin-tyrosine receptor kinase (Trk) complex formation and neurotrophin signaling. Together, our combined microarray and computational modeling methodologies provide a general, facile means to identify new glycosaminoglycan–protein–protein interactions, as well as a molecular-level understanding of those complexes.

The understanding of G-protein coupled receptors (GPCRs) is undergoing a revolution due to increased information about their signaling and the experimental determination of structures for more than 25 receptors. The availability of at least one receptor structure for each of the GPCR classes, well separated in sequence space, enables an integrated superfamily-wide analysis to identify signatures involving the role of conserved residues, conserved contacts, and downstream signaling in the context of receptor structures. In this study, we align the transmembrane (TM) domains of all experimental GPCR structures to maximize the conserved inter-helical contacts. The resulting superfamily-wide GpcR Sequence-Structure (GRoSS) alignment of the TM domains for all human GPCR sequences is sufficient to generate a phylogenetic tree that correctly distinguishes all different GPCR classes, suggesting that the class-level differences in the GPCR superfamily are encoded at least partly in the TM domains. The inter-helical contacts conserved across all GPCR classes describe the evolutionarily conserved GPCR structural fold. The corresponding structural alignment of the inactive and active conformations, available for a few GPCRs, identifies activation hot-spot residues in the TM domains that get rewired upon activation. Many GPCR mutations, known to alter receptor signaling and cause disease, are located at these conserved contact and activation hot-spot residue positions. The GRoSS alignment places the chemosensory receptor subfamilies for bitter taste (TAS2R) and pheromones (Vomeronasal, VN1R) in the rhodopsin family, known to contain the chemosensory olfactory receptor subfamily. The GRoSS alignment also enables the quantification of the structural variability in the TM regions of experimental structures, useful for homology modeling and structure prediction of receptors. Furthermore, this alignment identifies structurally and functionally important residues in all human GPCRs. These residues can be used to make testable hypotheses about the structural basis of receptor function and about the molecular basis of disease-associated single nucleotide polymorphisms. G-protein coupled receptors (GPCRs) are a large superfamily of integral membrane proteins that share a characteristic 7 transmembrane helix fold. They detect various molecules outside of the cell and signal their presence to the inside of the cell. At least half of the 800 human GPCRs are potential drug targets, so understanding their structure and function is critical. Experimental structures are now available for at least one receptor from each GPCR class. The structure of the 7 helix fold is highly conserved even for receptors with very low sequence similarity. We analyze the available experimental structures and compare the common inter-helical contacts. Our analysis leads to a unified sequence-structure alignment of the GPCR superfamily that can then be used as the starting point for structure prediction of all other GPCRs. A key result of our analysis is a list of conserved contact residues and activation “hot-spots” residues that are critical for GPCR folding and function. We propose that mutations and natural variants of amino acids at these locations in the GPCRs can dramatically influence their activation state and alter intracellular signaling. This provides hypotheses for the molecular mechanisms underlying disease causing mutants for any GPCR.

A major challenge in the application of structure-based drug design methods to proteins belonging to the superfamily of G protein-coupled receptors (GPCRs) is the paucity of structural information (1). The 19 chemokine receptors, belonging to the Class A family of GPCRs, are important drug targets not only for autoimmune diseases like multiple sclerosis but also for the blockade of human immunodeficiency virus type 1 entry (2). Using the MembStruk computational method (3), we predicted the three-dimensional structure of the human CCR1 receptor. In addition, we predicted the binding site of the small molecule CCR1 antagonist BX 471, which is currently in Phase II clinical trials (4). Based on the predicted antagonist binding site we designed 17 point mutants of CCR1 to validate the predictions. Subsequent competitive ligand binding and chemotaxis experiments with these mutants gave an excellent correlation to these predictions. In particular, we find that Tyr-113 and Tyr-114 on transmembrane domain 3 and Ile-259 on transmembrane 6 contribute significantly to the binding of BX 471. Finally, we used the predicted and validated structure of CCR1 in a virtual screening validation of the Maybridge data base, seeded with selective CCR1 antagonists. The screen identified 63% of CCR1 antagonists in the top 5% of the hits. Our results indicate that rational drug design for GPCR targets is a feasible approach.