The recent discoveries of microRNA (miRNA) genes and characterization of the first few target genes regulated by miRNAs in Caenorhabditis elegans and Drosophila melanogaster have set the stage for elucidation of a novel network of regulatory control. We present a computational method for whole-genome prediction of miRNA target genes. The method is validated using known examples. For each miRNA, target genes are selected on the basis of three properties: sequence complementarity using a position-weighted local alignment algorithm, free energies of RNA-RNA duplexes, and conservation of target sites in related genomes. Application to the D. melanogaster, Drosophila pseudoobscura and Anopheles gambiae genomes identifies several hundred target genes potentially regulated by one or more known miRNAs.

The discovery of cis-regulatory modules in metazoan genomes is crucial for understanding the connection between genes and organism diversity. It is important to quantify how comparative genomics can improve computational detection of such modules.

Cyclic nucleotide-gated (CNG) ion channels are nonselective cation channels, essential for visual and olfactory sensory transduction. Although the channels include voltage-sensor domains (VSDs), their conductance is thought to be independent of the membrane potential, and their gating regulated by cytosolic cyclic nucleotide-binding domains. Mutations in these channels result in severe, degenerative retinal diseases, which remain untreatable. The lack of structural information on CNG channels has prevented mechanistic understanding of disease-causing mutations, precluded structure-based drug design, and hampered in silico investigation of the gating mechanism. To address this, we built a 3D model of the cone tetrameric CNG channel, based on homology to two distinct templates with known structures: the transmembrane (TM) domain of a bacterial channel, and the cyclic nucleotide-binding domain of the mouse HCN2 channel. Since the TM-domain template had low sequence-similarity to the TM domains of the CNG channels, and to reconcile conflicts between the two templates, we developed a novel, hybrid approach, combining homology modeling with evolutionary coupling constraints. Next, we used elastic network analysis of the model structure to investigate global motions of the channel and to elucidate its gating mechanism. We found the following: (i) In the main mode of motion, the TM and cytosolic domains counter-rotated around the membrane normal. We related this motion to gating, a proposition that is supported by previous experimental data, and by comparison to the known gating mechanism of the bacterial KirBac channel. (ii) The VSDs could facilitate gating (supplementing the pore gate), explaining their presence in such 'voltage-insensitive' channels. (iii) Our elastic network model analysis of the CNGA3 channel supports a modular model of allosteric gating, according to which protein domains are quasi-independent: they can move independently, but are coupled to each other allosterically.

Cells adjust their behaviour in response to redox events by regulating protein activity through the reversible formation of disulfide bridges between cysteine thiols. However, the spatial and temporal control of these modifications remains poorly understood in multicellular organisms. Here we measured the protein thiol-disulfide balance in live Caenorhabditis elegans using a genetically encoded redox sensor and found that it is specific to tissues and is patterned spatially within a tissue. Insulin signalling regulates the sensor's oxidation at both of these levels. Unexpectedly, we found that isogenic individuals exhibit large differences in the sensor's thiol-disulfide balance. This variation contrasts with the general view that glutathione acts as the main cellular redox buffer. Indeed, our work suggests that glutathione converts small changes in its oxidation level into large changes in its redox potential. We therefore propose that glutathione facilitates the sensitive control of the thiol-disulfide balance of target proteins in response to cellular redox events.

Glia play crucial roles in the development and homeostasis of the nervous system. While the GLIA in the Drosophila embryo have been well characterized, their study in the adult nervous system has been limited. Here, we present a detailed description of the glia in the adult nervous system, based on the analysis of some 500 glial drivers we identified within a collection of synthetic GAL4 lines. We find that glia make up ∼10% of the cells in the nervous system and envelop all compartments of neurons (soma, dendrites, axons) as well as the nervous system as a whole. Our morphological analysis suggests a set of simple rules governing the morphogenesis of glia and their interactions with other cells. All glial subtypes minimize contact with their glial neighbors but maximize their contact with neurons and adapt their macromorphology and micromorphology to the neuronal entities they envelop. Finally, glial cells show no obvious spatial organization or registration with neuronal entities. Our detailed description of all glial subtypes and their regional specializations, together with the powerful genetic toolkit we provide, will facilitate the functional analysis of glia in the mature nervous system. GLIA 2017 GLIA 2017;65:606-638.

The removal of apoptotic cells by phagocytic neighbors is essential for metazoan development but remains poorly characterized. Here we report the discovery of a Drosophila phagocytosis receptor, Six-microns-under (SIMU), which is expressed in highly phagocytic cell types during development and required for efficient apoptotic cell clearance by glia in the nervous system and by macrophages elsewhere. SIMU is part of a conserved family of proteins that includes CED-1 and Draper (DRPR). Phenotypic analysis reveals that simu acts upstream of drpr in the same pathway and affects the recognition and engulfment of apoptotic cells, while drpr affects their subsequent degradation. SIMU strongly binds to apoptotic cells, presumably through its EMILIN-like domain, but requires no membrane anchoring, suggesting that it can function as a bridging molecule. Our study introduces an important factor in tissue-resident apoptotic clearance and underscores the prominent role of glia as "semiprofessional" phagocytes in the nervous system.

The complex patterns of gene expression in metazoans are controlled by selective binding of transcription factors (TFs) to regulatory DNA. To improve the quantitative understanding of this process, we have developed a novel method that uses fluorescence anisotropy measurements in a controlled delivery system to determine TF-DNA binding energies in solution with high sensitivity and throughput. Owing to its large dynamic range, the method, named high performance fluorescence anisotropy (HiP-FA), allows for reliable quantification of both weak and strong binding; binding specificities are calculated on the basis of equilibrium constant measurements for mutational DNA variants. We determine the binding preference landscapes for 26 TFs and measure high absolute affinities, but mostly lower binding specificities than reported by other methods. The revised binding preferences give rise to improved predictions of in vivo TF occupancy and enhancer expression. Our approach provides a powerful new tool for the systems-biological analysis of gene regulation.

Natural evolution encodes rich information about the structure and function of biomolecules in the genetic record. Previously, statistical analysis of co-variation patterns in natural protein families has enabled the accurate computation of 3D structures. Here, we explored generating similar information by experimental evolution, starting from a single gene and performing multiple cycles of in vitro mutagenesis and functional selection in Escherichia coli. We evolved two antibiotic resistance proteins, β-lactamase PSE1 and acetyltransferase AAC6, and obtained hundreds of thousands of diverse functional sequences. Using evolutionary coupling analysis, we inferred residue interaction constraints that were in agreement with contacts in known 3D structures, confirming genetic encoding of structural constraints in the selected sequences. Computational protein folding with interaction constraints then yielded 3D structures with the same fold as natural relatives. This work lays the foundation for a new experimental method (3Dseq) for protein structure determination, combining evolution experiments with inference of residue interactions from sequence information. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

The ability to edit the genome is essential for many state-of-the-art experimental paradigms. Since DNA breaks stimulate repair, they can be exploited to target site-specific integration. The clustered, regularly interspaced, short palindromic repeats (CRISPR)/cas9 system from Streptococcus pyogenes has been harnessed into an efficient and programmable nuclease for eukaryotic cells. We thus combined DNA cleavage by cas9, the generation of homologous recombination donors by polymerase chain reaction (PCR) and transient depletion of the non-homologous end joining factor lig4. Using cultured Drosophila melanogaster S2-cells and the phosphoglycerate kinase gene as a model, we reached targeted integration frequencies of up to 50% in drug-selected cell populations. Homology arms as short as 29 nt appended to the PCR primer resulted in detectable integration, slightly longer extensions are beneficial. We confirmed established rules for S. pyogenes cas9 sgRNA design and demonstrate that the complementarity region allows length variation and 5'-extensions. This enables generation of U6-promoter fusion templates by overlap-extension PCR with a standardized protocol. We present a series of PCR template vectors for C-terminal protein tagging and clonal Drosophila S2 cell lines with stable expression of a myc-tagged cas9 protein. The system can be used for epitope tagging or reporter gene knock-ins in an experimental setup that can in principle be fully automated.

Metabolic energy is stored in cells primarily as triacylglycerols in lipid droplets (LDs), and LD dysregulation leads to metabolic diseases. The formation of monolayer-bound LDs from the endoplasmic reticulum (ER) bilayer is poorly understood, but the ER protein seipin is essential to this process. In this study, we report a cryo-electron microscopy structure and functional characterization of Drosophila melanogaster seipin. The structure reveals a ring-shaped dodecamer with the luminal domain of each monomer resolved at ∼4.0 Å. Each luminal domain monomer exhibits two distinctive features: a hydrophobic helix (HH) positioned toward the ER bilayer and a β-sandwich domain with structural similarity to lipid-binding proteins. This structure and our functional testing in cells suggest a model in which seipin oligomers initially detect forming LDs in the ER via HHs and subsequently act as membrane anchors to enable lipid transfer and LD growth.

During the morphological process of sporulation in Bacillus subtilis two adjacent daughter cells (called the mother cell and forespore) follow different programs of gene expression that are linked to each other by signal transduction pathways. At a late stage in development, a signaling pathway emanating from the forespore triggers the proteolytic activation of the mother cell transcription factor σK. Cleavage of pro-σK to its mature and active form is catalyzed by the intramembrane cleaving metalloprotease SpoIVFB (B), a Site-2 Protease (S2P) family member. B is held inactive by two mother-cell membrane proteins SpoIVFA (A) and BofA. Activation of pro-σK processing requires a site-1 signaling protease SpoIVB (IVB) that is secreted from the forespore into the space between the two cells. IVB cleaves the extracellular domain of A but how this cleavage activates intramembrane proteolysis has remained unclear. Structural studies of the Methanocaldococcus jannaschii S2P homolog identified closed (substrate-occluded) and open (substrate-accessible) conformations of the protease, but the biological relevance of these conformations has not been established. Here, using co-immunoprecipitation and fluorescence microscopy, we show that stable association between the membrane-embedded protease and its substrate requires IVB signaling. We further show that the cytoplasmic cystathionine-β-synthase (CBS) domain of the B protease is not critical for this interaction or for pro-σK processing, suggesting the IVB-dependent interaction site is in the membrane protease domain. Finally, we provide evidence that the B protease domain adopts both open and closed conformations in vivo. Collectively, our data support a substrate-gating model in which IVB-dependent cleavage of A on one side of the membrane triggers a conformational change in the membrane-embedded protease from a closed to an open state allowing pro-σK access to the caged interior of the protease.

Recently it has been shown that the quality of protein contact prediction from evolutionary information can be improved significantly if direct and indirect information is separated. Given sufficiently large protein families, the contact predictions contain sufficient information to predict the structure of many protein families. However, since the first studies contact prediction methods have improved. Here, we ask how much the final models are improved if improved contact predictions are used.

Viral gain-of-function mutations frequently evolve during laboratory experiments. Whether the specific mutations that evolve in the lab also evolve in nature and whether they have the same impact on evolution in the real world is unknown. We studied a model virus, bacteriophage λ, that repeatedly evolves to exploit a new host receptor under typical laboratory conditions. Here, we demonstrate that two residues of λ's J protein are required for the new function. In natural λ variants, these amino acid sites are highly diverse and evolve at high rates. Insertions and deletions at these locations are associated with phylogenetic patterns indicative of ecological diversification. Our results show that viral evolution in the laboratory mirrors that in nature and that laboratory experiments can be coupled with protein sequence analyses to identify the causes of viral evolution in the real world. Furthermore, our results provide evidence for widespread host-shift evolution in lambdoid viruses.

Post-transcriptional regulation by microRNAs and siRNAs depends not only on characteristics of individual binding sites in target mRNA molecules, but also on system-level properties such as overall molecular concentrations. We hypothesize that an intracellular pool of microRNAs/siRNAs faced with a larger number of available predicted target transcripts will downregulate each individual target gene to a lesser extent. To test this hypothesis, we analyzed mRNA expression change from 178 microRNA and siRNA transfection experiments in two cell lines. We find that downregulation of particular genes mediated by microRNAs and siRNAs indeed varies with the total concentration of available target transcripts. We conclude that to interpret and design experiments involving gene regulation by small RNAs, global properties, such as target mRNA abundance, need to be considered in addition to local determinants. We propose that analysis of microRNA/siRNA targeting would benefit from a more quantitative definition, rather than simple categorization of genes as 'target' or 'not a target.' Our results are important for understanding microRNA regulation and may also have implications for siRNA design and small RNA therapeutics.

Vitamin K epoxide reductase (VKOR) drives the vitamin K cycle, activating vitamin K-dependent blood clotting factors. VKOR is also the target of the widely used anticoagulant drug, warfarin. Despite VKOR's pivotal role in coagulation, its structure and active site remain poorly understood. In addition, VKOR variants can cause vitamin K-dependent clotting factor deficiency or alter warfarin response. Here, we used multiplexed, sequencing-based assays to measure the effects of 2,695 VKOR missense variants on abundance and 697 variants on activity in cultured human cells. The large-scale functional data, along with an evolutionary coupling analysis, supports a four transmembrane domain topology, with variants in transmembrane domains exhibiting strongly deleterious effects on abundance and activity. Functionally constrained regions of the protein define the active site, and we find that, of four conserved cysteines putatively critical for function, only three are absolutely required. Finally, 25% of human VKOR missense variants show reduced abundance or activity, possibly conferring warfarin sensitivity or causing disease.

The transmembrane (TM) channel-like 1 (TMC1) and TMC2 proteins play a central role in auditory transduction, forming ion channels that convert sound into electrical signals. However, the molecular mechanism of their gating remains unknown. Here, using predicted structural models as a guide, we probed the effects of 12 mutations on the mechanical gating of the transduction currents in native hair cells of Tmc1/2-null mice expressing virally introduced TMC1 variants. Whole-cell electrophysiological recordings revealed that mutations within the pore-lining TM4 and TM6 helices modified gating, reducing the force sensitivity or shifting the open probability of the channels, or both. For some of the mutants, these changes were accompanied by a change in single-channel conductance. Our observations are in line with a model wherein conformational changes in the TM4 and TM6 helices are involved in the mechanical gating of the transduction channel.

In response to starvation, endospore-forming bacteria differentiate into stress-resistant spores that can remain dormant for years yet rapidly germinate and resume growth in response to nutrients. The small molecule dipicolinic acid (DPA) plays a central role in both the stress resistance of the dormant spore and its exit from dormancy during germination. The spoVA locus is required for DPA import during sporulation and has been implicated in its export during germination, but the molecular bases are unclear. Here, we define the minimal set of proteins encoded in the Bacillus subtilis spoVA operon required for DPA import and demonstrate that these proteins form a membrane complex. Structural modeling of these components combined with mutagenesis and in vivo analysis reveal that the C and Eb subunits form a membrane channel, while the D subunit functions as a cytoplasmic plug. We show that point mutations that impair the interactions between D and the C-Eb membrane complex reduce the efficiency of DPA import during sporulation and reciprocally accelerate DPA release during germination. Our data support a model in which DPA transport into spores involves cycles of unplugging and then replugging the C-Eb membrane channel, while nutrient detection during germination triggers DPA release by unplugging it.

Protocadherins (Pcdhs) are cell adhesion and signaling proteins used by neurons to develop and maintain neuronal networks, relying on trans homophilic interactions between their extracellular cadherin (EC) repeat domains. We present the structure of the antiparallel EC1-4 homodimer of human PcdhγB3, a member of the γ subfamily of clustered Pcdhs. Structure and sequence comparisons of α, β, and γ clustered Pcdh isoforms illustrate that subfamilies encode specificity in distinct ways through diversification of loop region structure and composition in EC2 and EC3, which contains isoform-specific conservation of primarily polar residues. In contrast, the EC1/EC4 interface comprises hydrophobic interactions that provide non-selective dimerization affinity. Using sequence coevolution analysis, we found evidence for a similar antiparallel EC1-4 interaction in non-clustered Pcdh families. We thus deduce that the EC1-4 antiparallel homodimer is a general interaction strategy that evolved before the divergence of these distinct protocadherin families.

Protein flexibility ranges from simple hinge movements to functional disorder. Around half of all human proteins contain apparently disordered regions with little 3D or functional information, and many of these proteins are associated with disease. Building on the evolutionary couplings approach previously successful in predicting 3D states of ordered proteins and RNA, we developed a method to predict the potential for ordered states for all apparently disordered proteins with sufficiently rich evolutionary information. The approach is highly accurate (79%) for residue interactions as tested in more than 60 known disordered regions captured in a bound or specific condition. Assessing the potential for structure of more than 1,000 apparently disordered regions of human proteins reveals a continuum of structural order with at least 50% with clear propensity for three- or two-dimensional states. Co-evolutionary constraints reveal hitherto unseen structures of functional importance in apparently disordered proteins.

The ability to design functional sequences and predict effects of variation is central to protein engineering and biotherapeutics. State-of-art computational methods rely on models that leverage evolutionary information but are inadequate for important applications where multiple sequence alignments are not robust. Such applications include the prediction of variant effects of indels, disordered proteins, and the design of proteins such as antibodies due to the highly variable complementarity determining regions. We introduce a deep generative model adapted from natural language processing for prediction and design of diverse functional sequences without the need for alignments. The model performs state-of-art prediction of missense and indel effects and we successfully design and test a diverse 105-nanobody library that shows better expression than a 1000-fold larger synthetic library. Our results demonstrate the power of the alignment-free autoregressive model in generalizing to regions of sequence space traditionally considered beyond the reach of prediction and design.