Three-dimensional architecture of the rod sensory cilium and its disruption in retinal neurodegeneration.

Defects in primary cilia lead to devastating disease because of their roles in sensation and developmental signaling but much is unknown about ciliary structure and mechanisms of their formation and maintenance. We used cryo-electron tomography to obtain 3D maps of the connecting cilium and adjacent cellular structures of a modified primary cilium, the rod outer segment, from wild-type and genetically defective mice. The results reveal the molecular architecture of the cilium and provide insights into protein functions. They suggest that the ciliary rootlet is involved in cellular transport and stabilizes the axoneme. A defect in the BBSome membrane coat caused defects in vesicle targeting near the base of the cilium. Loss of the proteins encoded by the Cngb1 gene disrupted links between the disk and plasma membranes. The structures of the outer segment membranes support a model for disk morphogenesis in which basal disks are enveloped by the plasma membrane.

Coevolutionary signals in metabotropic glutamate receptors capture residue contacts and long-range functional interactions.

Upon ligand binding to a G protein-coupled receptor, extracellular signals are transmitted into a cell through sets of residue interactions that translate ligand binding into structural rearrangements. These interactions needed for functions impose evolutionary constraints so that, on occasion, mutations in one position may be compensated by other mutations at functionally coupled positions. To quantify the impact of amino acid substitutions in the context of major evolutionary divergence in the G protein-coupled receptor subfamily of metabotropic glutamate receptors (mGluRs), we combined two phylogenetic-based algorithms, Evolutionary Trace and covariation Evolutionary Trace, to infer potential structure-function couplings and roles in mGluRs. We found a subset of evolutionarily important residues at known functional sites and evidence of coupling among distinct structural clusters in mGluR. In addition, experimental mutagenesis and functional assays confirmed that some highly covariant residues are coupled, revealing their synergy. Collectively, these findings inform a critical step toward understanding the molecular and structural basis of amino acid variation patterns within mGluRs and provide insight for drug development, protein engineering, and analysis of naturally occurring variants.

CovET: A covariation-evolutionary trace method that identifies protein structure-function modules.

Measuring the relative effect that any two sequence positions have on each other may improve protein design or help better interpret coding variants. Current approaches use statistics and machine learning but rarely consider phylogenetic divergences which, as shown by Evolutionary Trace studies, provide insight into the functional impact of sequence perturbations. Here, we reframe covariation analyses in the Evolutionary Trace framework to measure the relative tolerance to perturbation of each residue pair during evolution. This approach (CovET) systematically accounts for phylogenetic divergences: at each divergence event, we penalize covariation patterns that belie evolutionary coupling. We find that while CovET approximates the performance of existing methods to predict individual structural contacts, it performs significantly better at finding structural clusters of coupled residues and ligand binding sites. For example, CovET found more functionally critical residues when we examined the RNA recognition motif and WW domains. It correlates better with large-scale epistasis screen data. In the dopamine D2 receptor, top CovET residue pairs recovered accurately the allosteric activation pathway characterized for Class A G protein-coupled receptors. These data suggest that CovET ranks highest the sequence position pairs that play critical functional roles through epistatic and allosteric interactions in evolutionarily relevant structure-function motifs. CovET complements current methods and may shed light on fundamental molecular mechanisms of protein structure and function.

A synaptic vesicle-associated Ca2+ channel promotes endocytosis and couples exocytosis to endocytosis.

Synaptic vesicle (SV) exo- and endocytosis are tightly coupled to sustain neurotransmission in presynaptic terminals, and both are regulated by Ca(2+). Ca(2+) influx triggered by voltage-gated Ca(2+) channels is necessary for SV fusion. However, extracellular Ca(2+) has also been shown to be required for endocytosis. The intracellular Ca(2+) levels (<1 microM) that trigger endocytosis are typically much lower than those (>10 microM) needed to induce exocytosis, and endocytosis is inhibited when the Ca(2+) level exceeds 1 microM. Here, we identify and characterize a transmembrane protein associated with SVs that, upon SV fusion, localizes at periactive zones. Loss of Flower results in impaired intracellular resting Ca(2+) levels and impaired endocytosis. Flower multimerizes and is able to form a channel to control Ca(2+) influx. We propose that Flower functions as a Ca(2+) channel to regulate synaptic endocytosis and hence couples exo- with endocytosis.

Recurrent high-impact mutations at cognate structural positions in class A G protein-coupled receptors expressed in tumors.

G protein-coupled receptors (GPCRs) are the largest family of human proteins. They have a common structure and, signaling through a much smaller set of G proteins, arrestins, and effectors, activate downstream pathways that often modulate hallmark mechanisms of cancer. Because there are many more GPCRs than effectors, mutations in different receptors could perturb signaling similarly so as to favor a tumor. We hypothesized that somatic mutations in tumor samples may not be enriched within a single gene but rather that cognate mutations with similar effects on GPCR function are distributed across many receptors. To test this possibility, we systematically aggregated somatic cancer mutations across class A GPCRs and found a nonrandom distribution of positions with variant amino acid residues. Individual cancer types were enriched for highly impactful, recurrent mutations at selected cognate positions of known functional motifs. We also discovered that no single receptor drives this pattern, but rather multiple receptors contain amino acid substitutions at a few cognate positions. Phenotypic characterization suggests these mutations induce perturbation of G protein activation and/or ß-arrestin recruitment. These data suggest that recurrent impactful oncogenic mutations perturb different GPCRs to subvert signaling and promote tumor growth or survival. The possibility that multiple different GPCRs could moonlight as drivers or enablers of a given cancer through mutations located at cognate positions across GPCR paralogs opens a window into cancer mechanisms and potential approaches to therapeutics.

The mGluR6 ligand-binding domain, but not the C-terminal domain, is required for synaptic localization in retinal ON-bipolar cells.

Signals from retinal photoreceptors are processed in two parallel channels-the ON channel responds to light increments, while the OFF channel responds to light decrements. The ON pathway is mediated by ON type bipolar cells (BCs), which receive glutamatergic synaptic input from photoreceptors via a G-protein-coupled receptor signaling cascade. The metabotropic glutamate receptor mGluR6 is located at the dendritic tips of all ON-BCs and is required for synaptic transmission. Thus, it is critically important for delivery of information from photoreceptors into the ON pathway. In addition to detecting glutamate, mGluR6 participates in interactions with other postsynaptic proteins, as well as trans-synaptic interactions with presynaptic ELFN proteins. Mechanisms of mGluR6 synaptic targeting and functional interaction with other synaptic proteins are unknown. Here, we show that multiple regions in the mGluR6 ligand-binding domain are necessary for both synaptic localization in BCs and ELFN1 binding in vitro. However, these regions were not required for plasma membrane localization in heterologous cells, indicating that secretory trafficking and synaptic localization are controlled by different mechanisms. In contrast, the mGluR6 C-terminus was dispensable for synaptic localization. In mGluR6 null mice, localization of the postsynaptic channel protein TRPM1 was compromised. Introducing WT mGluR6 rescued TRPM1 localization, while a C-terminal deletion mutant had significantly reduced rescue ability. We propose a model in which trans-synaptic ELFN1 binding is necessary for mGluR6 postsynaptic localization, whereas the C-terminus has a role in mediating TRPM1 trafficking. These findings reveal different sequence determinants of the multifunctional roles of mGluR6 in ON-BCs.

Defining the layers of a sensory cilium with STORM and cryoelectron nanoscopy.

Primary cilia carry out numerous signaling and sensory functions, and defects in them, "ciliopathies," cause a range of symptoms, including blindness. Understanding of their nanometer-scale ciliary substructures and their disruptions in ciliopathies has been hindered by limitations of conventional microscopic techniques. We have combined cryoelectron tomography, enhanced by subtomogram averaging, with superresolution stochastic optical reconstruction microscopy (STORM) to define subdomains within the light-sensing rod sensory cilium of mouse retinas and reveal previously unknown substructures formed by resident proteins. Domains are demarcated by structural features such as the axoneme and its connections to the ciliary membrane, and are correlated with molecular markers of subcompartments, including the lumen and walls of the axoneme, the membrane glycocalyx, and the intervening cytoplasm. Within this framework, we report spatial distributions of key proteins in wild-type (WT) mice and the effects on them of genetic deficiencies in 3 models of Bardet-Biedl syndrome.

Structure and dynamics of photoreceptor sensory cilia.

The rod and cone photoreceptor cells of the vertebrate retina have highly specialized structures that enable them to carry out their function of light detection over a broad range of illumination intensities with optimized spatial and temporal resolution. Most prominent are their unusually large sensory cilia, consisting of outer segments packed with photosensitive disc membranes, a connecting cilium with many features reminiscent of the primary cilium transition zone, and a pair of centrioles forming a basal body which serves as the platform upon which the ciliary axoneme is assembled. These structures form a highway through which an enormous flux of material moves on a daily basis to sustain the continual turnover of outer segment discs and the energetic demands of phototransduction. After decades of study, the details of the fine structure and distribution of molecular components of these structures are still incompletely understood, but recent advances in cellular imaging techniques and animal models of inherited ciliary defects are yielding important new insights. This knowledge informs our understanding both of the mechanisms of trafficking and assembly and of the pathophysiological mechanisms of human blinding ciliopathies.

MTOR-initiated metabolic switch and degeneration in the retinal pigment epithelium.

The retinal pigment epithelium (RPE) is a particularly vulnerable tissue to age-dependent degeneration. Over the life span, the RPE develops an expanded endo-lysosomal compartment to maintain the high efficiency of phagocytosis and degradation of photoreceptor outer segments (POS) necessary for photoreceptor survival. As the assembly and activation of the mechanistic target of rapamycin complex 1 (mTORC1) occur on the lysosome surface, increased lysosome mass with aging leads to higher mTORC1 activity. The functional consequences of hyperactive mTORC1 in the RPE are unclear. In the current study, we used integrated high-resolution metabolomic and genomic approaches to examine mice with RPE-specific deletion of the tuberous sclerosis 1 (Tsc1) gene which encodes an upstream suppressor of mTORC1. Our data show that RPE cells with constitutively high mTORC1 activity were reprogramed to be hyperactive in glucose and lipid metabolism. Lipolysis was suppressed, mitochondrial carnitine shuttle was inhibited, while genes involved in fatty acid (FA) biosynthesis were upregulated. The metabolic changes occurred prior to structural changes of RPE and retinal degeneration. These findings have revealed cellular events and intrinsic mechanisms that contribute to lipid accumulation in the RPE cells during aging and age-related degeneration.

Residues and residue pairs of evolutionary importance differentially direct signaling bias of D2 dopamine receptors.

The D2 dopamine receptor and the serotonin 5-hydroxytryptamine 2A receptor (5-HT2A) are closely-related G-protein-coupled receptors (GPCRs) from the class A bioamine subfamily. Despite structural similarity, they respond to distinct ligands through distinct downstream pathways, whose dysregulation is linked to depression, bipolar disorder, addiction, and psychosis. They are important drug targets, and it is important to understand how their bias toward G-protein versus ß-arrestin signaling pathways is regulated. Previously, evolution-based computational approaches, difference Evolutionary Trace and Evolutionary Trace-Mutual information (ET-Mip), revealed residues and residue pairs that, when switched in the D2 receptor to the corresponding residues from 5-HT2A, altered ligand potency and G-protein activation efficiency. We have tested these residue swaps for their ability to trigger recruitment of ß-arrestin2 in response to dopamine or serotonin. The results reveal that the selected residues modulate agonist potency, maximal efficacy, and constitutive activity of ß-arrestin2 recruitment. Whereas dopamine potency for most variants was similar to that for WT and lower than for G-protein activation, potency in ß-arrestin2 recruitment for N124H3.42 was more than 5-fold higher. T205M5.54 displayed high constitutive activity, enhanced dopamine potency, and enhanced efficacy in ß-arrestin2 recruitment relative to WT, and L379F6.41 was virtually inactive. These striking differences from WT activity were largely reversed by a compensating mutation (T205M5.54/L379F6.41) at residues previously identified by ET-Mip as functionally coupled. The observation that the signs and relative magnitudes of the effects of mutations in several cases are at odds with their effects on G-protein activation suggests that they also modulate signaling bias.

Single-Atom Fluorescence Switch: A General Approach toward Visible-Light-Activated Dyes for Biological Imaging.

Photoactivatable fluorophores afford powerful molecular tools to improve the spatial and temporal resolution of subcellular structures and dynamics. By performing a single sulfur-for-oxygen atom replacement within common fluorophores, we have developed a facile and general strategy to obtain photoactivatable fluorogenic dyes across a broad spectral range. Thiocarbonyl substitution within fluorophores results in significant loss of fluorescence via a photoinduced electron transfer-quenching mechanism as suggested by theoretical calculations. Significantly, upon exposure to air and visible light residing in their absorption regime (365-630 nm), thio-caged fluorophores can be efficiently desulfurized to their oxo derivatives, thus restoring strong emission of the fluorophores. The effective photoactivation makes thio-caged fluorophores promising candidates for super-resolution imaging, which was realized by photoactivated localization microscopy (PALM) with low-power activation light under physiological conditions in the absence of cytotoxic additives (e.g., thiols, oxygen scavengers), a feature superior to traditional PALM probes. The versatility of this thio-caging strategy was further demonstrated by multicolor super-resolution imaging of lipid droplets and proteins of interest.

Integrative subcellular proteomic analysis allows accurate prediction of human disease-causing genes.

Proteomic profiling on subcellular fractions provides invaluable information regarding both protein abundance and subcellular localization. When integrated with other data sets, it can greatly enhance our ability to predict gene function genome-wide. In this study, we performed a comprehensive proteomic analysis on the light-sensing compartment of photoreceptors called the outer segment (OS). By comparing with the protein profile obtained from the retina tissue depleted of OS, an enrichment score for each protein is calculated to quantify protein subcellular localization, and 84% accuracy is achieved compared with experimental data. By integrating the protein OS enrichment score, the protein abundance, and the retina transcriptome, the probability of a gene playing an essential function in photoreceptor cells is derived with high specificity and sensitivity. As a result, a list of genes that will likely result in human retinal disease when mutated was identified and validated by previous literature and/or animal model studies. Therefore, this new methodology demonstrates the synergy of combining subcellular fractionation proteomics with other omics data sets and is generally applicable to other tissues and diseases.

Intramolecular allosteric communication in dopamine D2 receptor revealed by evolutionary amino acid covariation.

The structural basis of allosteric signaling in G protein-coupled receptors (GPCRs) is important in guiding design of therapeutics and understanding phenotypic consequences of genetic variation. The Evolutionary Trace (ET) algorithm previously proved effective in redesigning receptors to mimic the ligand specificities of functionally distinct homologs. We now expand ET to consider mutual information, with validation in GPCR structure and dopamine D2 receptor (D2R) function. The new algorithm, called ET-MIp, identifies evolutionarily relevant patterns of amino acid covariations. The improved predictions of structural proximity and D2R mutagenesis demonstrate that ET-MIp predicts functional interactions between residue pairs, particularly potency and efficacy of activation by dopamine. Remarkably, although most of the residue pairs chosen for mutagenesis are neither in the binding pocket nor in contact with each other, many exhibited functional interactions, implying at-a-distance coupling. The functional interaction between the coupled pairs correlated best with the evolutionary coupling potential derived from dopamine receptor sequences rather than with broader sets of GPCR sequences. These data suggest that the allosteric communication responsible for dopamine responses is resolved by ET-MIp and best discerned within a short evolutionary distance. Most double mutants restored dopamine response to wild-type levels, also suggesting that tight regulation of the response to dopamine drove the coevolution and intramolecular communications between coupled residues. Our approach provides a general tool to identify evolutionary covariation patterns in small sets of close sequence homologs and to translate them into functional linkages between residues.

ß2-Adrenergic receptor activation mobilizes intracellular calcium via a non-canonical cAMP-independent signaling pathway.

Beta adrenergic receptors (ßARs) are G-protein-coupled receptors essential for physiological responses to the hormones/neurotransmitters epinephrine and norepinephrine which are found in the nervous system and throughout the body. They are the targets of numerous widely used drugs, especially in the case of the most extensively studied ßAR, ß2AR, whose ligands are used for asthma and cardiovascular disease. ßARs signal through Gαs G-proteins and via activation of adenylyl cyclase and cAMP-dependent protein kinase, but some alternative downstream pathways have also been proposed that could be important for understanding normal physiological functioning of ßAR signaling and its disruption in disease. Using fluorescence-based Ca2+ flux assays combined with pharmacology and gene knock-out methods, we discovered a previously unrecognized endogenous pathway in HEK-293 cells whereby ß2AR activation leads to robust Ca2+ mobilization from intracellular stores via activation of phospholipase C and opening of inositol trisphosphate (InsP3) receptors. This pathway did not involve cAMP, Gαs, or Gαi or the participation of the other members of the canonical ß2AR signaling cascade and, therefore, constitutes a novel signaling mechanism for this receptor. This newly uncovered mechanism for Ca2+ mobilization by ß2AR has broad implications for adrenergic signaling, cross-talk with other signaling pathways, and the effects of ßAR-directed drugs.

Correction: the retromer complex is required for rhodopsin recycling and its loss leads to photoreceptor degeneration.

[This corrects the article DOI: 10.1371/journal.pbio.1001847.].

Differential epitope masking reveals synapse-specific complexes of TRPM1.

The transient receptor potential channel TRPM1 is required for synaptic transmission between photoreceptors and the ON subtype of bipolar cells (ON-BPC), mediating depolarization in response to light. TRPM1 is present in the somas and postsynaptic dendritic tips of ON-BPCs. Monoclonal antibodies generated against full-length TRPM1 were found to have differential labeling patterns when used to immunostain the mouse retina, with some yielding reduced labeling of dendritic tips relative to the labeling of cell bodies. Epitope mapping revealed that those antibodies that poorly label the dendritic tips share a binding site (N2d) in the N-terminal arm near the transmembrane domain. A major splice variant of TRPM1 lacking exon 19 does not contain the N2d binding site, but quantitative immunoblotting revealed no enrichment of this variant in synaptsomes. One explanation of the differential labeling is masking of the N2d epitope by formation of a synapse-specific multiprotein complex. Identifying the binding partners that are specific for the fraction of TRPM1 present at the synapses is an ongoing challenge for understanding TRPM1 function.

The retromer complex is required for rhodopsin recycling and its loss leads to photoreceptor degeneration.

Rhodopsin mistrafficking can cause photoreceptor (PR) degeneration. Upon light exposure, activated rhodopsin 1 (Rh1) in Drosophila PRs is internalized via endocytosis and degraded in lysosomes. Whether internalized Rh1 can be recycled is unknown. Here, we show that the retromer complex is expressed in PRs where it is required for recycling endocytosed Rh1 upon light stimulation. In the absence of subunits of the retromer, Rh1 is processed in the endolysosomal pathway, leading to a dramatic increase in late endosomes, lysosomes, and light-dependent PR degeneration. Reducing Rh1 endocytosis or Rh1 levels in retromer mutants alleviates PR degeneration. In addition, increasing retromer abundance suppresses degenerative phenotypes of mutations that affect the endolysosomal system. Finally, expressing human Vps26 suppresses PR degeneration in Vps26 mutant PRs. We propose that the retromer plays a conserved role in recycling rhodopsins to maintain PR function and integrity.

Determinants of endogenous ligand specificity divergence among metabotropic glutamate receptors.

To determine the structural origins of diverse ligand response specificities among metabotropic glutamate receptors (mGluRs), we combined computational approaches with mutagenesis and ligand response assays to identify specificity-determining residues in the group I receptor, mGluR1, and the group III receptors, mGluR4 and mGluR7. Among these, mGluR1 responds to L-glutamate effectively, whereas it binds weakly to another endogenous ligand, L-serine-O-phosphate (L-SOP), which antagonizes the effects of L-glutamate. In contrast, mGluR4 has in common with other group III mGluR that it is activated with higher potency and efficacy by L-SOP. mGluR7 differs from mGluR4 and other group III mGluR in that L-glutamate and L-SOP activate it with low potency and efficacy. Enhanced versions of the evolutionary trace (ET) algorithm were used to identify residues that when swapped between mGluR1 and mGluR4 increased the potency of L-SOP inhibition relative to the potency of L-glutamate activation in mGluR1 mutants and others that diminished the potency/efficacy of L-SOP for mGluR4 mutants. In addition, combining ET identified swaps from mGluR4 with one identified by computational docking produced mGluR7 mutants that respond with dramatically enhanced potency/efficacy to L-SOP. These results reveal that an early functional divergence between group I/II and group III involved variation at positions primarily at allosteric sites located outside of binding pockets, whereas a later divergence within group III occurred through sequence variation both at the ligand-binding pocket and at loops near the dimerization interface and interlobe hinge region. They also demonstrate the power of ET for identifying allosteric determinants of evolutionary importance.

Domain organization and conformational plasticity of the G protein effector, PDE6.

The cGMP phosphodiesterase of rod photoreceptor cells, PDE6, is the key effector enzyme in phototransduction. Two large catalytic subunits, PDE6α and -ß, each contain one catalytic domain and two non-catalytic GAF domains, whereas two small inhibitory PDE6γ subunits allow tight regulation by the G protein transducin. The structure of holo-PDE6 in complex with the ROS-1 antibody Fab fragment was determined by cryo-electron microscopy. The ∼11 Å map revealed previously unseen features of PDE6, and each domain was readily fit with high resolution structures. A structure of PDE6 in complex with prenyl-binding protein (PrBP/δ) indicated the location of the PDE6 C-terminal prenylations. Reconstructions of complexes with Fab fragments bound to N or C termini of PDE6γ revealed that PDE6γ stretches from the catalytic domain at one end of the holoenzyme to the GAF-A domain at the other. Removal of PDE6γ caused dramatic structural rearrangements, which were reversed upon its restoration.

Nonsense mutations in the rhodopsin gene that give rise to mild phenotypes trigger mRNA degradation in human cells by nonsense-mediated decay.

Eight different nonsense mutations in the human rhodopsin gene cause retinitis pigmentosa (RP), an inherited degenerative disease of the retina that can lead to complete blindness. Although all these nonsense mutations lead to premature termination codons (PTCs) in rhodopsin mRNA, some display dominant inheritance, while others are recessive. Because nonsense-mediated decay (NMD) can degrade mRNAs containing PTCs and modulate the inheritance patterns of genetic diseases, we asked whether any of the nonsense mutations in the rhodopsin gene generated mRNAs that were susceptible to degradation by NMD. We hypothesized that nonsense mutations that caused mild RP phenotypes would trigger NMD, whereas those that did not engage NMD would cause more severe RP phenotypes-presumably due to the toxicity of the truncated protein. To test our hypothesis, we transfected human rhodopsin nonsense mutants into HEK293 and HT1080 human cells and measured transcript levels by qRT-PCR. In both cell lines, rhodopsin mutations Q64X and Q344X, which cause severe phenotypes that are dominantly inherited, yielded the same levels of rhodopsin mRNA as wild type. By contrast, rhodopsin mutations W161X and E249X, which cause recessive RP, showed decreased rhodopsin mRNA levels, consistent with NMD. Rhodopsin mutant Y136X, a dominant mutation that causes late-onset RP with a very mild pathology, also gave lower mRNA levels. Treatment of cells with Wortmannin, an inhibitor of NMD, eliminated the degradation of Y136X, W161X, and E249X rhodopsin mRNAs. These results suggest that NMD modulates the severity of RP in patients with nonsense mutations in the rhodopsin gene.