Using HHsearch to tackle proteins of unknown function: A pilot study with PH domains

A refined picture of the native amine dehydrogenase family revealed by extensive biodiversity screening

Native amine dehydrogenases offer sustainable access to chiral amines, so the search for scaffolds capable of converting more diverse carbonyl compounds is required to reach the full potential of this alternative to conventional synthetic reductive aminations. Here we report a multidisciplinary strategy combining bioinformatics, chemoinformatics and biocatalysis to extensively screen billions of sequences in silico and to efficiently find native amine dehydrogenases features using computational approaches. In this way, we achieve a comprehensive overview of the initial native amine dehydrogenase family, extending it from 2,011 to 17,959 sequences, and identify native amine dehydrogenases with non-reported substrate spectra, including hindered carbonyls and ethyl ketones, and accepting methylamine and cyclopropylamine as amine donor. We also present preliminary model-based structural information to inform the design of potential (R)-selective amine dehydrogenases, as native amine dehydrogenases are mostly (S)-selective. This integrated strategy paves the way for expanding the resource of other enzyme families and in highlighting enzymes with original features.

The GTPase activating protein Gyp7 regulates Rab7/Ypt7 activity on late endosomes

Organelles of the endomembrane system contain Rab GTPases as identity markers. Their localization is determined by guanine nucleotide exchange factors (GEFs) and GTPase activating proteins (GAPs). It remains largely unclear how these regulators are specifically targeted to organelles and how their activity is regulated. Here, we focus on the GAP Gyp7, which acts on the Rab7-like Ypt7 protein in yeast, and surprisingly observe the protein exclusively in puncta proximal to the vacuole. Mistargeting of Gyp7 to the vacuole strongly affects vacuole morphology, suggesting that endosomal localization is needed for function. In agreement, efficient endolysosomal transport requires Gyp7. In vitro assays reveal that Gyp7 requires a distinct lipid environment for membrane binding and activity. Overexpression of Gyp7 concentrates Ypt7 in late endosomes and results in resistance to rapamycin, an inhibitor of the target of rapamycin complex 1 (TORC1), suggesting that these late endosomes are signaling endosomes. We postulate that Gyp7 is part of regulatory machinery involved in late endosome function.

Acidithiobacillia class members originating at sites within the Pacific Ring of Fire and other tectonically active locations and description of the novel genus 'Igneacidithiobacillus'

Recent studies have expanded the genomic contours of the Acidithiobacillia, highlighting important lacunae in our comprehension of the phylogenetic space occupied by certain lineages of the class. One such lineage is 'Igneacidithiobacillus', a novel genus-level taxon, represented by 'Igneacidithiobacillus copahuensis' VAN18-1T as its type species, along with two other uncultivated metagenome-assembled genomes (MAGs) originating from geothermally active sites across the Pacific Ring of Fire. In this study, we investigate the genetic and genomic diversity, and the distribution patterns of several uncharacterized Acidithiobacillia class strains and sequence clones, which are ascribed to the same 16S rRNA gene sequence clade. By digging deeper into this data and contributing to novel MAGs emerging from environmental studies in tectonically active locations, the description of this novel genus has been consolidated. Using state-of-the-art genomic taxonomy methods, we added to already recognized taxa, an additional four novel Candidate (Ca.) species, including 'Ca. Igneacidithiobacillus chanchocoensis' (mCHCt20-1TS), 'Igneacidithiobacillus siniensis' (S30A2T), 'Ca. Igneacidithiobacillus taupoensis' (TVZ-G3 TS), and 'Ca. Igneacidithiobacillus waiarikiensis' (TVZ-G4 TS). Analysis of published data on the isolation, enrichment, cultivation, and preliminary microbiological characterization of several of these unassigned or misassigned strains, along with the type species of the genus, plus the recoverable environmental data from metagenomic studies, allowed us to identify habitat preferences of these taxa. Commonalities and lineage-specific adaptations of the seven species of the genus were derived from pangenome analysis and comparative genomic metabolic reconstruction. The findings emerging from this study lay the groundwork for further research on the ecology, evolution, and biotechnological potential of the novel genus 'Igneacidithiobacillus'.

Evolutionary history of the cytochrome P450s from Colletotrichum species and prediction of their putative functional roles during host-pathogen interactions

The genomes of species belonging to the genus Colletotrichum harbor a substantial number of cytochrome P450 monooxygenases (CYPs) encoded by a broad diversity of gene families. However, the biological role of their CYP complement (CYPome) has not been elucidated. Here, we investigated the putative evolutionary scenarios that occurred during the evolution of the CYPome belonging to the Colletotrichum Graminicola species complex (s.c.) and their biological implications. The study revealed that most of the CYPome gene families belonging to the Graminicola s.c. experienced gene contractions. The reductive evolution resulted in species restricted CYPs are predominant in each CYPome of members from the Graminicola s.c., whereas only 18 families are absolutely conserved among these species. However, members of CYP families displayed a notably different phylogenetic relationship at the tertiary structure level, suggesting a putative convergent evolution scenario. Most of the CYP enzymes of the Graminicola s.c. share redundant functions in secondary metabolite biosynthesis and xenobiotic metabolism. Hence, this current work suggests that the presence of a broad CYPome in the genus Colletotrichum plays a critical role in the optimization of the colonization capability and virulence.

Pangenome-level analysis of nucleoid-associated proteins in the Acidithiobacillia class: insights into their functional roles in mobile genetic elements biology

Mobile genetic elements (MGEs) are relevant agents in bacterial adaptation and evolutionary diversification. Stable appropriation of these DNA elements depends on host factors, among which are the nucleoid-associated proteins (NAPs). NAPs are highly abundant proteins that bind and bend DNA, altering its topology and folding, thus affecting all known cellular DNA processes from replication to expression. Even though NAP coding genes are found in most prokaryotic genomes, their functions in host chromosome biology and xenogeneic silencing are only known for a few NAP families. Less is known about the occurrence, abundance, and roles of MGE-encoded NAPs in foreign elements establishment and mobility. In this study, we used a combination of comparative genomics and phylogenetic strategies to gain insights into the diversity, distribution, and functional roles of NAPs within the class Acidithiobacillia with a special focus on their role in MGE biology. Acidithiobacillia class members are aerobic, chemolithoautotrophic, acidophilic sulfur-oxidizers, encompassing substantial genotypic diversity attributable to MGEs. Our search for NAP protein families (PFs) in more than 90 genomes of the different species that conform the class, revealed the presence of 1,197 proteins pertaining to 12 different NAP families, with differential occurrence and conservation across species. Pangenome-level analysis revealed 6 core NAP PFs that were highly conserved across the class, some of which also existed as variant forms of scattered occurrence, in addition to NAPs of taxa-restricted distribution. Core NAPs identified are reckoned as essential based on the conservation of genomic context and phylogenetic signals. In turn, various highly diversified NAPs pertaining to the flexible gene complement of the class, were found to be encoded in known plasmids or, larger integrated MGEs or, present in genomic loci associated with MGE-hallmark genes, pointing to their role in the stabilization/maintenance of these elements in strains and species with larger genomes. Both core and flexible NAPs identified proved valuable as markers, the former accurately recapitulating the phylogeny of the class, and the later, as seed in the bioinformatic identification of novel episomal and integrated mobile elements.

Endoplasmic Reticulum Membrane Contact Sites, Lipid Transport, and Neurodegeneration

The Endoplasmic Reticulum (ER) is an endomembrane system that plays a multiplicity of roles in cell physiology and populates even the most distal cell compartments, including dendritic tips and axon terminals of neurons. Some of its functions are achieved by a cross talk with other intracellular membranous organelles and with the plasma membrane at membrane contacts sites (MCSs). As the ER synthesizes most membrane lipids, lipid exchanges mediated by lipid transfer proteins at MCSs are a particularly important aspect of this cross talk, which synergizes with the cross talk mediated by vesicular transport. Several mutations of genes that encode proteins localized at ER MCSs result in familial neurodegenerative diseases, emphasizing the importance of the normal lipid traffic within cells for a healthy brain. Here, we provide an overview of such diseases, with a specific focus on proteins that directly or indirectly impact lipid transport.

In silico modeling human VPS13 proteins associated with donor and target membranes suggests lipid transfer mechanisms

The VPS13 protein family constitutes a novel class of bridge-like lipid transferases. Autosomal recessive inheritance of mutations in VPS13 genes is associated with the development of neurodegenerative diseases in humans. Bioinformatic approaches previously recognized the domain architecture of these proteins. In this study, we model the first ever full-length structures of the four human homologs VPS13A, VPS13B, VPS13C, and VPS13D in association with model membranes, to investigate their lipid transfer ability and potential structural association with membrane leaflets. We analyze the evolutionary conservation and physicochemical properties of these proteins, focusing on conserved C-terminal amphipathic helices that disturb organelle surfaces and that, adjoined, resemble a traditional Venetian gondola. The gondola domains share significant structural homology with lipid droplet surface-binding proteins. We introduce in silico protein-membrane models displaying the mode of association of VPS13A, VPS13B, VPS13C, and VPS13D to donor and target membranes, and present potential models of action for protein-mediated lipid transfer.

Global phylogenomic novelty of the Cas1 gene from hot spring microbial communities

The Cas1 protein is essential for the functioning of CRISPR-Cas adaptive systems. However, despite the high prevalence of CRISPR-Cas systems in thermophilic microorganisms, few studies have investigated the occurrence and diversity of Cas1 across hot spring microbial communities. Phylogenomic analysis of 2,150 Cas1 sequences recovered from 48 metagenomes representing hot springs (42–80°C, pH 6–9) from three continents, revealed similar ecological diversity of Cas1 and 16S rRNA associated with geographic location. Furthermore, phylogenetic analysis of the Cas1 sequences exposed a broad taxonomic distribution in thermophilic bacteria, with new clades of Cas1 homologs branching at the root of the tree or at the root of known clades harboring reference Cas1 types. Additionally, a new family of casposases was identified from hot springs, which further completes the evolutionary landscape of the Cas1 superfamily. This ecological study contributes new Cas1 sequences from known and novel locations worldwide, mainly focusing on under-sampled hot spring microbial mat taxa. Results herein show that circumneutral hot springs are environments harboring high diversity and novelty related to adaptive immunity systems.

Computational tools to study RNA-protein complexes

RNA is the key player in many cellular processes such as signal transduction, replication, transport, cell division, transcription, and translation. These diverse functions are accomplished through interactions of RNA with proteins. However, protein–RNA interactions are still poorly derstood in contrast to protein–protein and protein–DNA interactions. This knowledge gap can be attributed to the limited availability of protein-RNA structures along with the experimental difficulties in studying these complexes. Recent progress in computational resources has expanded the number of tools available for studying protein-RNA interactions at various molecular levels. These include tools for predicting interacting residues from primary sequences, modelling of protein-RNA complexes, predicting hotspots in these complexes and insights into derstanding in the dynamics of their interactions. Each of these tools has its strengths and limitations, which makes it significant to select an optimal approach for the question of interest. Here we present a mini review of computational tools to study different aspects of protein-RNA interactions, with focus on overall application, development of the field and the future perspectives.

Identification and functional characterization of ORF19.5274, a novel gene involved in both azoles susceptibility and hypha development in Candida albicans

Azole resistance is becoming increasingly serious due to the frequent recurrence of fungal infections and the need for long-term clinical prevention. In our previous study, we discovered ORF19.5274 with an unknown function by TMT^™ quantitative proteomics technology after fluconazole (FLC) treatment of Candida albicans. In this study, we created the target gene deletion strain using CRISPR-Cas9 editing technology to see if ORF19.5274 regulates azole sensitivity. The data showed that ORF19.5274 was involved in hyphal development and susceptibility to antifungal azoles. Deleting this gene resulted in defective hyphal growth in solid medium, while only a weak lag in the initiation of hyphal development and restoring hyphal growth during the hyphal maintenance phase under liquid conditions. Moreover, intracellular reactive oxygen species (ROS) assay and propidium iodide staining assays showed increased endogenous ROS levels and membrane permeability, but decreased metabolic activity of biofilm in orf19.5274Δ/Δ after treatment with FLC in comparison with either SC5314 or orf19.5274Δ/Δ::ORF19.5274 strains. More importantly, orf19.5274Δ/Δ significantly enhanced the FLC efficacy against C. albicans in infected Galleria mellonella larvae. The above characteristics were fully or partially restored in the complemented strain indicating that the changes caused by ORF19.5274 deletion were specific. In summary, the ORF19.5274 gene is required for hyphal development of C. albicans, and is correlated with the response to antifungal azoles in vitro and in vivo. The identification of ORF19.5275 is promising to expand the potential candidate targets for azoles.

A partnership between the lipid scramblase XK and the lipid transfer protein VPS13A at the plasma membrane

Chorea-acanthocytosis (ChAc) and McLeod syndrome are diseases with shared clinical manifestations caused by mutations in VPS13A and XK, respectively. Key features of these conditions are the degeneration of caudate neurons and the presence of abnormally shaped erythrocytes. XK belongs to a family of plasma membrane (PM) lipid scramblases whose action results in exposure of PtdSer at the cell surface. VPS13A is an endoplasmic reticulum (ER)-anchored lipid transfer protein with a putative role in the transport of lipids at contacts of the ER with other membranes. Recently VPS13A and XK were reported to interact by still unknown mechanisms. So far, however, there is no evidence for a colocalization of the two proteins at contacts of the ER with the PM, where XK resides, as VPS13A was shown to be localized at contacts between the ER and either mitochondria or lipid droplets. Here we show that VPS13A can also localize at ER-PM contacts via the binding of its PH domain to a cytosolic loop of XK, that such interaction is regulated by an intramolecular interaction within XK, and that both VPS13A and XK are highly expressed in the caudate neurons. Binding of the PH domain of VPS13A to XK is competitive with its binding to intracellular membranes that mediate other tethering functions of VPS13A. Our findings support a model according to which VPS13A-dependent lipid transfer between the ER and the PM is coupled to lipid scrambling within the PM. They raise the possibility that defective cell surface exposure of PtdSer may be responsible for neurodegeneration.

Disease relevance of rare VPS13B missense variants for neurodevelopmental Cohen syndrome

Autosomal recessive Cohen syndrome is a neurodevelopmental disorder characterized by postnatal microcephaly, intellectual disability, and a typical facial gestalt. Genetic variants in VPS13B have been found to cause Cohen syndrome, but have also been linked to autism, retinal disease, primary immunodeficiency, and short stature. While it is well established that loss-of-function mutations of VPS13B cause Cohen syndrome, the relevance of missense variants for the pathomechanism remains unexplained. Here, we investigate their pathogenic effect through a systematic re-evaluation of clinical patient information, comprehensive in silico predictions, and in vitro testing of previously published missense variants. In vitro analysis of 10 subcloned VPS13B missense variants resulted in full-length proteins after transient overexpression. 6/10 VPS13B missense variants show reduced accumulation at the Golgi complex in the steady state. The overexpression of these 6/10 VPS13B missense variants did not rescue the Golgi fragmentation after the RNAi-mediated depletion of endogenous VPS13B. These results thus validate 6/10 missense variants as likely pathogenic according to the classification of the American College of Medical Genetics through the integration of clinical, genetic, in silico, and experimental data. In summary, we state that exact variant classification should be the first step towards elucidating the pathomechanisms of genetically inherited neuronal diseases.

Structural and biochemical insights into lipid transport by VPS13 proteins

VPS13 proteins are proposed to function at contact sites between organelles as bridges for lipids to move directionally and in bulk between organellar membranes. VPS13s are anchored between membranes via interactions with receptors, including both peripheral and integral membrane proteins. Here we present the crystal structure of VPS13s adaptor binding domain (VAB) complexed with a Pro-X-Pro peptide recognition motif present in one such receptor, the integral membrane protein Mcp1p, and show biochemically that other Pro-X-Pro motifs bind the VAB in the same site. We further demonstrate that Mcp1p and another integral membrane protein that interacts directly with human VPS13A, XK, are scramblases. This finding supports an emerging paradigm of a partnership between bulk lipid transport proteins and scramblases. Scramblases can re-equilibrate lipids between membrane leaflets as lipids are removed from or inserted into the cytosolic leaflet of donor and acceptor organelles, respectively, in the course of protein-mediated transport.

Sequence Analysis and Structural Predictions of Lipid Transfer Bridges in the Repeating Beta Groove (RBG) Superfamily Reveal Past and Present Domain Variations Affecting Form, Function and Interactions of VPS13, ATG2, SHIP164, Hobbit and Tweek

Lipid transfer between organelles requires proteins that shield the hydrophobic portions of lipids as they cross the cytoplasm. In the last decade a new structural form of lipid transfer protein (LTP) has been found: long hydrophobic grooves made of beta-sheet that bridge between organelles at membrane contact sites. Eukaryotes have five families of bridge-like LTPs: VPS13, ATG2, SHIP164, Hobbit and Tweek. These are unified into a single superfamily through their bridges being composed of just one domain, called the repeating beta groove (RBG) domain, which builds into rod shaped multimers with a hydrophobic-lined groove and hydrophilic exterior. Here, sequences and predicted structures of the RBG superfamily were analyzed in depth. Phylogenetics showed that the last eukaryotic common ancestor contained all five RBG proteins, with duplicated VPS13s. The current set of long RBG protein appears to have arisen in even earlier ancestors from shorter forms with 4 RBG domains. The extreme ends of most RBG proteins have amphipathic helices that might be an adaptation for direct or indirect bilayer interaction, although this has yet to be tested. The one exception to this is the C-terminus of SHIP164, which instead has a coiled-coil. Finally, the exterior surfaces of the RBG bridges are shown to have conserved residues along most of their length, indicating sites for partner interactions almost all of which are unknown. These findings can inform future cell biological and biochemical experiments.

The GTPase Arf1 Is a Determinant of Yeast Vps13 Localization to the Golgi Apparatus

VPS13 proteins are evolutionarily conserved. Mutations in the four human genes (VPS13A-D) encoding VPS13A-D proteins are linked to developmental or neurodegenerative diseases. The relationship between the specific localization of individual VPS13 proteins, their molecular functions, and the pathology of these diseases is unknown. Here we used a yeast model to establish the determinants of Vps13's interaction with the membranes of Golgi apparatus. We analyzed the different phenotypes of the arf1-3 arf2Δ vps13∆ strain, with reduced activity of the Arf1 GTPase, the master regulator of Golgi function and entirely devoid of Vps13. Our analysis led us to propose that Vps13 and Arf1 proteins cooperate at the Golgi apparatus. We showed that Vps13 binds to the Arf1 GTPase through its C-terminal Pleckstrin homology (PH)-like domain. This domain also interacts with phosphoinositol 4,5-bisphosphate as it was bound to liposomes enriched with this lipid. The homologous domain of VPS13A exhibited the same behavior. Furthermore, a fusion of the PH-like domain of Vps13 to green fluorescent protein was localized to Golgi structures in an Arf1-dependent manner. These results suggest that the PH-like domains and Arf1 are determinants of the localization of VPS13 proteins to the Golgi apparatus in yeast and humans.

COSMIC Cancer Gene Census 3D database: understanding the impacts of mutations on cancer targets

Mutations in hallmark genes are believed to be the main drivers of cancer progression. These mutations are reported in the Catalogue of Somatic Mutations in Cancer (COSMIC). Structural appreciation of where these mutations appear, in protein-protein interfaces, active sites or deoxyribonucleic acid (DNA) interfaces, and predicting the impacts of these mutations using a variety of computational tools are crucial for successful drug discovery and development. Currently, there are 723 genes presented in the COSMIC Cancer Gene Census. Due to the complexity of the gene products, structures of only 87 genes have been solved experimentally with structural coverage between 90% and 100%. Here, we present a comprehensive, user-friendly, web interface (https://cancer-3d.com/) of 714 modelled cancer-related genes, including homo-oligomers, hetero-oligomers, transmembrane proteins and complexes with DNA, ribonucleic acid, ligands and co-factors. Using SDM and mCSM software, we have predicted the impacts of reported mutations on protein stability, protein-protein interfaces affinity and protein-nucleic acid complexes affinity. Furthermore, we also predicted intrinsically disordered regions using DISOPRED3.

Mapping the glycosyltransferase fold landscape using interpretable deep learning

Glycosyltransferases (GTs) play fundamental roles in nearly all cellular processes through the biosynthesis of complex carbohydrates and glycosylation of diverse protein and small molecule substrates. The extensive structural and functional diversification of GTs presents a major challenge in mapping the relationships connecting sequence, structure, fold and function using traditional bioinformatics approaches. Here, we present a convolutional neural network with attention (CNN-attention) based deep learning model that leverages simple secondary structure representations generated from primary sequences to provide GT fold prediction with high accuracy. The model learns distinguishing secondary structure features free of primary sequence alignment constraints and is highly interpretable. It delineates sequence and structural features characteristic of individual fold types, while classifying them into distinct clusters that group evolutionarily divergent families based on shared secondary structural features. We further extend our model to classify GT families of unknown folds and variants of known folds. By identifying families that are likely to adopt novel folds such as GT91, GT96 and GT97, our studies expand the GT fold landscape and prioritize targets for future structural studies.

A VPS13D spastic ataxia mutation disrupts the conserved adaptor-binding site in yeast Vps13

Mutations in each of the four human VPS13 (VPS13A-D) proteins are associated with distinct neurological disorders: chorea-acanthocytosis, Cohen syndrome, early-onset Parkinson's disease and spastic ataxia. Recent evidence suggests that the different VPS13 paralogs transport lipids between organelles at different membrane contact sites. How each VPS13 isoform is targeted to organelles is not known. We have shown that the localization of yeast Vps13 protein to membranes requires a conserved six-repeat region, the Vps13 Adaptor Binding (VAB) domain, which binds to organelle-specific adaptors. Here, we use a systematic mutagenesis strategy to determine the role of each repeat in recognizing each known adaptor. Our results show that mutation of invariant asparagines in repeats 1 and 6 strongly impacts the binding of all adaptors and blocks Vps13 membrane recruitment. However, we find that repeats 5-6 are sufficient for localization and interaction with adaptors. This supports a model where a single adaptor-binding site is found in the last two repeats of the VAB domain, while VAB domain repeat 1 may influence domain conformation. Importantly, a disease-causing mutation in VPS13D, which maps to the highly conserved asparagine residue in repeat 6, blocks adaptor binding and Vps13 membrane recruitment when modeled in yeast. Our findings are consistent with a conserved adaptor binding role for the VAB domain and suggest the presence of as-yet-unidentified adaptors in both yeast and humans.

PDZD8-mediated lipid transfer at contacts between the ER and late endosomes/lysosomes is required for neurite outgrowth

Membrane contact sites (MCSs) between the endoplasmic reticulum (ER) and late endosomes/lysosomes (LE/lys) are emerging as critical hubs for diverse cellular events, and changes in their extents are linked to severe neurological diseases. While recent studies show that the synaptotagmin-like mitochondrial-lipid-binding (SMP) domain-containing protein PDZD8 may mediate the formation of ER-LE/lys MCSs, the cellular functions of PDZD8 remain largely elusive. Here, we attempt to investigate the lipid transfer activities of PDZD8 and the extent to which its cellular functions depend on its lipid transfer activities. In accordance with recent studies, we demonstrate that PDZD8 is a protrudin (ZFYVE27)-interacting protein and that PDZD8 acts as a tether at ER-LE/lys MCSs. Furthermore, we discover that the SMP domain of PDZD8 binds glycerophospholipids and ceramides both in vivo and in vitro, and that the SMP domain can transport lipids between membranes in vitro. Functionally, PDZD8 is required for LE/lys positioning and neurite outgrowth, which is dependent on the lipid transfer activity of the SMP domain.

SARS-CoV-2 3D database: understanding the coronavirus proteome and evaluating possible drug targets

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a rapidly growing infectious disease, widely spread with high mortality rates. Since the release of the SARS-CoV-2 genome sequence in March 2020, there has been an international focus on developing target-based drug discovery, which also requires knowledge of the 3D structure of the proteome. Where there are no experimentally solved structures, our group has created 3D models with coverage of 97.5% and characterized them using state-of-the-art computational approaches. Models of protomers and oligomers, together with predictions of substrate and allosteric binding sites, protein-ligand docking, SARS-CoV-2 protein interactions with human proteins, impacts of mutations, and mapped solved experimental structures are freely available for download. These are implemented in SARS CoV-2 3D, a comprehensive and user-friendly database, available at https://sars3d.com/. This provides essential information for drug discovery, both to evaluate targets and design new potential therapeutics.

The Vps13 Family of Lipid Transporters and Its Role at Membrane Contact Sites

The conserved VPS13 proteins constitute a new family of lipid transporters at membrane contact sites. These large proteins are suspected to bridge membranes and form a direct channel for lipid transport between organelles. Mutations in the 4 human homologs (VPS13A–D) are associated with a number of neurological disorders, but little is known about their precise functions or the relevant contact sites affected in disease. In contrast, yeast has a single Vps13 protein which is recruited to multiple organelles and contact sites. The yeast model system has proved useful for studying the function of Vps13 at different organelles and identifying the localization determinants responsible for its membrane targeting. In this review we describe recent advances in our understanding of VPS13 proteins with a focus on yeast research.

Role of VPS13, a protein with similarity to ATG2, in physiology and disease

The evolutionarily conserved VPS13 family proteins have been implicated in several cellular processes. Mutations in each of the four human VPS13s cause neurodevelopmental or neurodegenerative disorders. Until recently, the molecular function of VPS13 remained elusive. Genetic, functional and structural studies have now revealed that VPS13 acts at contact sites between intracellular organelles to transport lipids by a novel mechanism: direct transfer between bilayers via a hydrophobic channel that spans its entire rod-like N-terminal half. Predicted similarities to the autophagy protein ATG2 suggested a similar role for ATG2 that has now been confirmed by structural and functional studies. Here, after a brief review of this evidence, we discuss what is known of human VPS13 proteins in physiology and disease.

Characterisation of plasmodial transketolases and identification of potential inhibitors: an in silico study

BACKGROUND

Plasmodial transketolase (PTKT) enzyme is one of the novel pharmacological targets being explored as potential anti-malarial drug target due to its functional role and low sequence identity to the human enzyme. Despite this, features contributing to such have not been exploited for anti-malarial drug design. Additionally, there are no anti-malarial drugs targeting PTKTs whereas the broad activity of these inhibitors against PTKTs from other Plasmodium spp. is yet to be reported. This study characterises different PTKTs [Plasmodium falciparum (PfTKT), Plasmodium vivax (PvTKT), Plasmodium ovale (PoTKT), Plasmodium malariae (PmTKT) and Plasmodium knowlesi (PkTKT) and the human homolog (HsTKT)] to identify key sequence and structural based differences as well as the identification of selective potential inhibitors against PTKTs.

METHODS

A sequence-based study was carried out using multiple sequence alignment, phylogenetic tree calculations and motif discovery analysis. Additionally, TKT models of PfTKT, PmTKT, PoTKT, PmTKT and PkTKT were modelled using the Saccharomyces cerevisiae TKT structure as template. Based on the modelled structures, molecular docking using 623 South African natural compounds was done. The stability, conformational changes and detailed interactions of selected compounds were accessed viz all-atom molecular dynamics (MD) simulations and binding free energy (BFE) calculations.

RESULTS

Sequence alignment, evolutionary and motif analyses revealed key differences between plasmodial and the human TKTs. High quality homodimeric three-dimensional PTKTs structures were constructed. Molecular docking results identified three compounds (SANC00107, SANC00411 and SANC00620) which selectively bind in the active site of all PTKTs with the lowest (better) binding affinity ≤ - 8.5 kcal/mol. MD simulations of ligand-bound systems showed stable fluctuations upon ligand binding. In all systems, ligands bind stably throughout the simulation and form crucial interactions with key active site residues. Simulations of selected compounds in complex with human TKT showed that ligands exited their binding sites at different time steps. BFE of protein-ligand complexes showed key residues involved in binding.

CONCLUSIONS

This study highlights significant differences between plasmodial and human TKTs and may provide valuable information for the development of novel anti-malarial inhibitors. Identified compounds may provide a starting point in the rational design of PTKT inhibitors and analogues based on these scaffolds.

Structure-Based Deep Mining Reveals First-Time Annotations for 46 Percent of the Dark Annotation Space of the 9,671-Member Superproteome of the Nucleocytoplasmic Large DNA Viruses

We conducted an exhaustive search for three-dimensional structural homologs to the proteins of 20 key phylogenetically distinct nucleocytoplasmic DNA viruses (NCLDV). Structural matches covered 429 known protein domain superfamilies, with the most highly represented being ankyrin repeat, P-loop NTPase, F-box, protein kinase, and membrane occupation and recognition nexus (MORN) repeat. Domain superfamily diversity correlated with genome size, but a diversity of around 200 superfamilies appeared to correlate with an abrupt switch to paralogization. Extensive structural homology was found across the range of eukaryotic RNA polymerase II subunits and their associated basal transcription factors, with the coordinated gain and loss of clusters of subunits on a virus-by-virus basis. The total number of predicted endonucleases across the 20 NCLDV was nearly quadrupled from 36 to 132, covering much of the structural and functional diversity of endonucleases throughout the biosphere in DNA restriction, repair, and homing. Unexpected findings included capsid protein-transcription factor chimeras; endonuclease chimeras; enzymes for detoxification; antimicrobial peptides and toxin-antitoxin systems associated with symbiosis, immunity, and addiction; and novel proteins for membrane abscission and protein turnover.IMPORTANCE We extended the known annotation space for the NCLDV by 46%, revealing high-probability structural matches for fully 45% of the 9,671 query proteins and confirming up to 98% of existing annotations per virus. The most prevalent protein families included ankyrin repeat- and MORN repeat-containing proteins, many of which included an F-box, suggesting extensive host cell modulation among the NCLDV. Regression suggested a minimum requirement for around 36 protein structural superfamilies for a viable NCLDV, and beyond around 200 superfamilies, genome expansion by the acquisition of new functions was abruptly replaced by paralogization. We found homologs to herpesvirus surface glycoprotein gB in cytoplasmic viruses. This study provided the first prediction of an endonuclease in 10 of the 20 viruses examined; the first report in a virus of a phenolic acid decarboxylase, proteasomal subunit, or cysteine knot (defensin) protein; and the first report of a prokaryotic-type ribosomal protein in a eukaryotic virus.

Structural bioinformatics predicts that the Retinitis Pigmentosa-28 protein of unknown function FAM161A is a homologue of the microtubule nucleation factor Tpx2

Background: FAM161A is a microtubule-associated protein conserved widely across eukaryotes, which is mutated in the inherited blinding disease Retinitis Pigmentosa-28. FAM161A is also a centrosomal protein, being a core component of a complex that forms an internal skeleton of centrioles. Despite these observations about the importance of FAM161A, current techniques used to examine its sequence reveal no homologies to other proteins.

Methods: Sequence profiles derived from multiple sequence alignments of FAM161A homologues were constructed by PSI-BLAST and HHblits, and then used by the profile-profile search tool HHsearch, implemented online as HHpred, to identify homologues. These in turn were used to create profiles for reverse searches and pair-wise searches. Multiple sequence alignments were also used to identify amino acid usage in functional elements.

Results: FAM161A has a single homologue: the targeting protein for Xenopus kinesin-like protein-2 (Tpx2), which is a strong hit across more than 200 residues. Tpx2 is also a microtubule-associated protein, and it has been shown previously by a cryo-EM molecular structure to nucleate microtubules through two small elements: an extended loop and a short helix. The homology between FAM161A and Tpx2 includes these elements, as FAM161A has three copies of the loop, and one helix that has many, but not all, properties of the one in Tpx2.

Conclusions: FAM161A and ­its homologues are predicted to be a previously unknown variant of Tpx2, and hence bind microtubules in the same way. This prediction allows precise, testable molecular models to be made of FAM161A-microtubule complexes.

Peroxisome retention involves Inp1-dependent peroxisome-plasma membrane contact sites in yeast

Retention of peroxisomes in yeast mother cells requires Inp1, which is recruited to the organelle by the peroxisomal membrane protein Pex3. Here we show that Hansenula polymorpha Inp1 associates peroxisomes to the plasma membrane. Peroxisome-plasma membrane contact sites disappear upon deletion of INP1 but increase upon INP1 overexpression. Analysis of truncated Inp1 variants showed that the C terminus is important for association to the peroxisome, while a stretch of conserved positive charges and a central pleckstrin homology-like domain are important for plasma membrane binding. In cells of a PEX3 deletion, strain Inp1-GFP localizes to the plasma membrane, concentrated in patches near the bud neck and in the cortex of nascent buds. Upon disruption of the actin cytoskeleton by treatment of the cells with latrunculin A, Inp1-GFP became cytosolic, indicating that Inp1 localization is dependent on the presence of an intact actin cytoskeleton.

IMG/VR v.2.0: an integrated data management and analysis system for cultivated and environmental viral genomes

The Integrated Microbial Genome/Virus (IMG/VR) system v.2.0 (https://img.jgi.doe.gov/vr/) is the largest publicly available data management and analysis platform dedicated to viral genomics. Since the last report published in the 2016, NAR Database Issue, the data has tripled in size and currently contains genomes of 8389 cultivated reference viruses, 12 498 previously published curated prophages derived from cultivated microbial isolates, and 735 112 viral genomic fragments computationally predicted from assembled shotgun metagenomes. Nearly 60% of the viral genomes and genome fragments are clustered into 110 384 viral Operational Taxonomic Units (vOTUs) with two or more members. To improve data quality and predictions of host specificity, IMG/VR v.2.0 now separates prokaryotic and eukaryotic viruses, utilizes known prophage sequences to improve taxonomic assignments, and provides viral genome quality scores based on the estimated genome completeness. New features also include enhanced BLAST search capabilities for external queries. Finally, geographic map visualization to locate user-selected viral genomes or genome fragments has been implemented and download options have been extended. All of these features make IMG/VR v.2.0 a key resource for the study of viruses.

Novel pathogenic VPS13A mutation in Moroccan family with Choreoacanthocytosis: a case report

Background

Choreoacanthocytosis (ChAc), is a rare neurodegenerative disease, characterized by movement disorders and acanthocytosis in the peripheral blood smears, and various neurological, neuropsychiatric and neuromuscular signs. It is caused by mutations in VPS13A gene with autosomal recessive pattern of inheritance.

Case presentation

Here we report two patients belonging to a consanguineous Moroccan family who present with movement disorder pathology. They were suspected to have choreoacanthocytosis according to biological, clinical and radiological finding. Thus, whole-exome sequencing was performed for precise diagnosis and identified a homozygous novel nonsense mutation c.337C > T (p.Gln113*) in exon 5 of VPS13A in the two affected siblings.

Conclusion

Here, we report a novel nonsense p.Gln113* mutation in VPS13A identified by whole-exome sequencing, which caused ChAc in a Moroccan family. This is the first description of ChAc in Morocco with genetic confirmation, that expands the mutation diversity of VPS13A and provide clinical, neuroimaging and deep brain stimulation findings.

Structural characterization of the RH1-LZI tandem of JIP3/4 highlights RH1 domains as a cytoskeletal motor-binding motif

JIP3 and JIP4 (JNK-interacting proteins 3 and 4) are adaptors for cargo recruitment by dynein/dynactin and kinesin1 motors. Both are dimers that are stabilised by two sections of leucine zipper coiled coils. The N-terminal Leucine Zipper I (LZI) belongs to a section that binds dynein-DLIC and kinesin1-KHC, whilst the medial Leucine Zipper II (LZII) binds dynactin-p150glued and kinesin1-KLC. Structural data is available for the LZII, but the LZI section is still uncharacterized. Here we characterize the N-terminal part of JIP3/4 which consists of an RH1 (RILP homology 1) domain followed by the LZI coiled coil using bioinformatical, biophysical and structural approaches. The RH1-LZI tandem of JIP3 associates as a high affinity homodimer exhibiting elongated alpha-helical fold. 3D homology modelling of the RH1-LZI tandem reveals that the kinesin1-KHC binding site mainly overlaps with the RH1 domain. A sequence comparison search indicates that only one other protein family has RH1 domains similar to those of JIP3/4, the RILP (Rab-interacting lysosomal protein) family which consists of adaptor proteins linking Rab GTPases to cytoskeletal motors. RILPL2 is recruited through its RH1 domain by the myosin 5a motor. Here, we showed that the RH1 domain of JIP3 also interacts with myosin 5 A in vitro, highlighting JIP3/4 as possible myosin 5a adaptors. Finally, we propose that JIP3/4 and RILP family members define a unique RH1/RH2-architecture adaptor superfamily linking cytoskeletal motors and Rab GTPases.

Systematic prediction of FFAT motifs across eukaryote proteomes identifies nucleolar and eisosome proteins with the predicted capacity to form bridges to the endoplasmic reticulum

The endoplasmic reticulum (ER), the most pervasive organelle, exchanges information and material with many other organelles, but the extent of its inter-organelle connections and the proteins that form bridges are not well known. The integral ER membrane protein VAMP-associated protein (VAP) is found in multiple bridges, interacting with many proteins that contain a short linear motif consisting of "two phenylalanines in an acidic tract" (FFAT). The VAP-FFAT interaction is the most common mechanism by which cytoplasmic proteins, particularly inter-organelle bridges, target the ER. Therefore, predicting new FFAT motifs may both find new individual peripheral ER proteins and identify new routes of communication involving the ER. Here we searched for FFAT motifs across whole proteomes. The excess of eukaryotic proteins with FFAT motifs over background was ≥0.8%, suggesting this is the minimum number of peripheral ER proteins. In yeast, where VAP was previously known to bind 4 proteins with FFAT motifs, a detailed analysis of a subset of proteins predicted 20 FFAT motifs. Extrapolating these findings to the whole proteome estimated the number of FFAT motifs in yeast at approximately 50-55 (0.9% of proteome). Among these previously unstudied FFAT motifs, most have known functions outside the ER, so could be involved in inter-organelle communication. Many of these can target well-characterised membrane contact sites, however some are in nucleoli and eisosomes, organelles previously unknown to have molecular bridges to the ER. We speculate that the nucleolar and eisosomal proteins with predicted motifs may function while bridging to the ER, indicating novel ER-nucleolus and ER-eisosome routes of inter-organelle communication.

Remote homology searches identify bacterial homologues of eukaryotic lipid transfer proteins, including Chorein-N domains in TamB and AsmA and Mdm31p

Background

All cells rely on lipids for key functions. Lipid transfer proteins allow lipids to exit the hydrophobic environment of bilayers, and cross aqueous spaces. One lipid transfer domain fold present in almost all eukaryotes is the TUbular LIPid binding (TULIP) domain. Three TULIP families have been identified in bacteria (P47, OrfX2 and YceB), but their homology to eukaryotic proteins is too low to specify a common origin. Another recently described eukaryotic lipid transfer domain in VPS13 and ATG2 is Chorein-N, which has no known bacterial homologues. There has been no systematic search for bacterial TULIPs or Chorein-N domains.

Results

Remote homology predictions for bacterial TULIP domains using HHsearch identified four new TULIP domains in three bacterial families. DUF4403 is a full length pseudo-dimeric TULIP with a 6 strand β-meander dimer interface like eukaryotic TULIPs. A similar sheet is also present in YceB, suggesting it homo-dimerizes. TULIP domains were also found in DUF2140 and in the C-terminus DUF2993. Remote homology predictions for bacterial Chorein-N domains identified strong hits in the N-termini of AsmA and TamB in diderm bacteria, which are related to Mdm31p in eukaryotic mitochondria. The N-terminus of DUF2993 has a Chorein-N domain adjacent to its TULIP domain.

Conclusions

TULIP lipid transfer domains are widespread in bacteria. Chorein-N domains are also found in bacteria, at the N-terminus of multiple proteins in the intermembrane space of diderms (AsmA, TamB and their relatives) and in Mdm31p, a protein that is likely to have evolved from an AsmA/TamB-like protein in the endosymbiotic mitochondrial ancestor. This indicates that both TULIP and Chorein-N lipid transfer domains may have originated in bacteria.

Planar Cell Polarity Effector Proteins Inturned and Fuzzy Form a Rab23 GEF Complex

A subset of Rab GTPases have been implicated in cilium formation in cultured mammalian cells [1-6]. Rab11 and Rab8, together with their GDP-GTP exchange factors (GEFs), TRAPP-II and Rabin8, promote recruitment of the ciliary vesicle to the mother centriole and its subsequent maturation, docking, and fusion with the cell surface [2-5]. Rab23 has been linked to cilium formation and membrane trafficking at mature cilia [1, 7, 8]; however, the identity of the GEF pathway activating Rab23, a member of the Rab7 subfamily of Rabs, remains unclear. Longin-domain-containing complexes have been shown to act as GEFs for Rab7 subfamily GTPases [9-12]. Here, we show that Inturned and Fuzzy, proteins previously implicated as planar cell polarity (PCP) effectors and in developmentally regulated cilium formation [13, 14], contain multiple longin domains characteristic of the Mon1-Ccz1 family of Rab7 GEFs and form a specific Rab23 GEF complex. In flies, loss of Rab23 function gave rise to defects in planar-polarized trichome formation consistent with this biochemical relationship. In cultured human and mouse cells, Inturned and Fuzzy localized to the basal body and proximal region of cilia, and cilium formation was compromised by depletion of either Inturned or Fuzzy. Cilium formation arrested after docking of the ciliary vesicle to the mother centriole but prior to axoneme elongation and fusion of the ciliary vesicle and plasma membrane. These findings extend the family of longin domain GEFs and define a molecular activity linking Rab23-regulated membrane traffic to cilia and planar cell polarity.

HH-suite3 for fast remote homology detection and deep protein annotation

Background

HH-suite is a widely used open source software suite for sensitive sequence similarity searches and protein fold recognition. It is based on pairwise alignment of profile Hidden Markov models (HMMs), which represent multiple sequence alignments of homologous proteins.

Results

We developed a single-instruction multiple-data (SIMD) vectorized implementation of the Viterbi algorithm for profile HMM alignment and introduced various other speed-ups. These accelerated the search methods HHsearch by a factor 4 and HHblits by a factor 2 over the previous version 2.0.16. HHblits3 is ∼10× faster than PSI-BLAST and ∼20× faster than HMMER3. Jobs to perform HHsearch and HHblits searches with many query profile HMMs can be parallelized over cores and over cluster servers using OpenMP and message passing interface (MPI). The free, open-source, GPLv3-licensed software is available at https://github.com/soedinglab/hh-suite.

Conclusion

The added functionalities and increased speed of HHsearch and HHblits should facilitate their use in large-scale protein structure and function prediction, e.g. in metagenomics and genomics projects.

Nucleotide Second Messenger-Based Signaling in Extreme Acidophiles of the Acidithiobacillus Species Complex: Partition Between the Core and Variable Gene Complements

Cyclic and linear nucleotides are key elements of the signal transduction networks linking perception of the environment to specific cellular behavior of prokaryotes. These molecular mechanisms are particularly important in bacteria exposed to different, and frequently simultaneous, types of extreme conditions. This is the case in acidithiobacilli, a group of extremophilic bacteria thriving in highly acidic biotopes, that must also cope with significant variations in temperature, osmotic potentials and concentrations of various transition metals and metalloids. Environmental cues sensed by bacteria are transduced into differential levels of nucleotides acting as intracellular second messengers, promoting the activation or inhibition of target components and eliciting different output phenotypes. Cyclic (c) di-GMP, one of the most common bacterial second messengers, plays a key role in lifestyle changes in many bacteria, including acidithiobacilli. The presence of functional c-di-GMP-dependent signal transduction pathways in representative strains of the best-known linages of this species complex has been reported. However, a comprehensive panorama of the c-di-GMP modulated networks, the cognate input signals and output responses, are still missing for this group of extremophiles. Moreover, little fundamental understanding has been gathered for other nucleotides acting as second messengers. Taking advantage of the increasing number of sequenced genomes of the taxon, here we address the challenge of disentangling the nucleotide-driven signal transduction pathways in this group of polyextremophiles using comparative genomic tools and strategies. Results indicate that the acidithiobacilli possess all the genetic elements required to establish functional transduction pathways based in three different nucleotide-second messengers: (p)ppGpp, cyclic AMP (cAMP), and c-di-GMP. The elements related with the metabolism and transduction of (p)ppGpp and cAMP appear highly conserved, integrating signals related with nutrient starvation and polyphosphate metabolism, respectively. In contrast, c-di-GMP networks appear diverse and complex, differing both at the species and strain levels. Molecular elements of c-di-GMP metabolism and transduction were mostly found scattered along the flexible genome of the acidithiobacilli, allowing the identification of probable control modules that could be critical for substrate colonization, biofilm development and intercellular interactions. These may ultimately convey increased endurance to environmental stress and increased potential for gene sharing and adaptation to changing conditions.

VPS13A and VPS13C are lipid transport proteins differentially localized at ER contact sites

Structural and biochemical studies show that VPS13 family proteins are lipid transporters. VPS13A and VPS13C tether the ER to mitochondria and to endosomes/lysosomes, respectively, suggesting lipid dyshomeostasis as the cause of chorea acanthocytosis and Parkinson’s disease resulting from their mutations.

Mutations in the human VPS13 genes are responsible for neurodevelopmental and neurodegenerative disorders including chorea acanthocytosis (VPS13A) and Parkinson’s disease (VPS13C). The mechanisms of these diseases are unknown. Genetic studies in yeast hinted that Vps13 may have a role in lipid exchange between organelles. In this study, we show that the N-terminal portion of VPS13 is tubular, with a hydrophobic cavity that can solubilize and transport glycerolipids between membranes. We also show that human VPS13A and VPS13C bind to the ER, tethering it to mitochondria (VPS13A), to late endosome/lysosomes (VPS13C), and to lipid droplets (both VPS13A and VPS13C). These findings identify VPS13 as a lipid transporter between the ER and other organelles, implicating defects in membrane lipid homeostasis in neurological disorders resulting from their mutations. Sequence and secondary structure similarity between the N-terminal portions of Vps13 and other proteins such as the autophagy protein ATG2 suggest lipid transport roles for these proteins as well.

Competitive organelle-specific adaptors recruit Vps13 to membrane contact sites

Targeting of Vps13 to membranes is highly dynamic. Bean et al. identify Ypt35 and Mcp1 as adaptors for Vps13 at endosomes and mitochondria, respectively, and show all known Vps13 adaptors use a related motif to compete for Vps13 membrane recruitment.

The regulated expansion of membrane contact sites, which mediate the nonvesicular exchange of lipids between organelles, requires the recruitment of additional contact site proteins. Yeast Vps13 dynamically localizes to membrane contacts that connect the ER, mitochondria, endosomes, and vacuoles and is recruited to the prospore membrane in meiosis, but its targeting mechanism is unclear. In this study, we identify the sorting nexin Ypt35 as a novel adaptor that recruits Vps13 to endosomal and vacuolar membranes. We characterize an interaction motif in the Ypt35 N terminus and identify related motifs in the prospore membrane adaptor Spo71 and the mitochondrial membrane protein Mcp1. We find that Mcp1 is a mitochondrial adaptor for Vps13, and the Vps13–Mcp1 interaction, but not Ypt35, is required when ER-mitochondria contacts are lost. All three adaptors compete for binding to a conserved six-repeat region of Vps13 implicated in human disease. Our results support a competition-based model for regulating Vps13 localization at cellular membranes.

Yeast and other lower eukaryotic organisms for studies of Vps13 proteins in health and disease

Human Vps13 proteins are associated with several diseases, including the neurodegenerative disorder Chorea-acanthocytosis (ChAc), yet the biology of these proteins is still poorly understood. Studies in Saccharomyces cerevisiae, Dictyostelium discoideum, Tetrahymena thermophila and Drosophila melanogaster point to the involvement of Vps13 in cytoskeleton organization, vesicular trafficking, autophagy, phagocytosis, endocytosis, proteostasis, sporulation and mitochondrial functioning. Recent findings show that yeast Vps13 binds to phosphatidylinositol lipids via 4 different regions and functions at membrane contact sites, enlarging the list of Vps13 functions. This review describes the great potential of simple eukaryotes to decipher disease mechanisms in higher organisms and highlights novel insights into the pathological role of Vps13 towards ChAc.

Vps13-Mcp1 interact at vacuole-mitochondria interfaces and bypass ER-mitochondria contact sites

Interorganelle membrane contacts work in a networked fashion to allow exchange of metabolites throughout the cell. In yeast, mitochondria–vacuole contacts act redundantly with ER–mitochondria contacts. We show that the yeast mitochondrial protein Mcp1 binds the endosomal/vacuolar protein Vps13 to mediate the physiological function of vacuole–mitochondria contacts.

Membrane contact sites between endoplasmic reticulum (ER) and mitochondria, mediated by the ER–mitochondria encounter structure (ERMES) complex, are critical for mitochondrial homeostasis and cell growth. Defects in ERMES can, however, be bypassed by point mutations in the endosomal protein Vps13 or by overexpression of the mitochondrial protein Mcp1. How this bypass operates remains unclear. Here we show that the mitochondrial outer membrane protein Mcp1 functions in the same pathway as Vps13 by recruiting it to mitochondria and promoting its association to vacuole–mitochondria contacts. Our findings support a model in which Mcp1 and Vps13 work as functional effectors of vacuole–mitochondria contact sites, while tethering is mediated by other factors, including Vps39. Tethered and functionally active vacuole–mitochondria interfaces then compensate for the loss of ERMES-mediated ER–mitochondria contact sites.

Tubular lipid binding proteins (TULIPs) growing everywhere

Tubular lipid binding proteins (TULIPs) have become a focus of interest in the cell biology of lipid signalling, lipid traffic and membrane contact sites. Each tubular domain has an internal pocket with a hydrophobic lining that can bind a hydrophobic molecule such as a lipid. This allows TULIP proteins to carry lipids through the aqueous phase. TULIP domains were first found in a large family of extracellular proteins related to the bacterial permeability-inducing protein (BPI) and cholesterol ester transfer protein (CETP). Since then, the same fold and lipid transfer capacity have been found in SMP domains (so-called for their occurrence in synaptotagmin, mitochondrial and lipid binding proteins), which localise to intracellular membrane contact sites. Here the methods for identifying known TULIPs are described, and used to find previously unreported TULIPs, one in the silk polymer and another in prokaryotes illustrated by the E. coli protein YceB. The bacterial TULIP alters views on the likely evolution of the domain, suggesting its presence in the last universal common ancestor. The major function of TULIPs is to handle lipids, but we still do not know how they work in detail, or how many more remain to be discovered. This article is part of a Special Issue entitled: Membrane Contact Sites edited by Christian Ungermann and Benoit Kornmann.

•

Proteins with the tubular lipid binding fold exist in a wider variety than is usually appreciated.

•

TULIPs are found in prokaryotes, altering views on their evolution.

•

It is not yet known whether TULIPs transfer lipids as tunnels or as shuttles.

•

Tests have not yet been done to say if TULIPs with SMP domains (for example E-syts and ERMES components) tether contact sites.

•

It is likely that more TULIPs remain to be discovered.

Vps13 and Cdc31/centrin: Puzzling partners in membrane traffic

Myers and Payne discuss work by De et al. describing the unique structure of Vps13 and its role in trafficking

Yeast Vps13 is a member of a conserved protein family that includes human homologues associated with neurodegenerative and developmental disorders. In this issue, De et al. (2017. J. Cell Biol.https://doi.org/10.1083/jcb.201606078) establish direct roles for Vps13 and its surprising binding partner, the calcium-binding centrin Cdc31, in trans-Golgi network (TGN) to endosome traffic and TGN homotypic fusion.