The MPI bioinformatics Toolkit as an integrative platform for advanced protein sequence and structure analysis

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.

The Untapped Australasian Diversity of Astaxanthin-Producing Yeasts with Biotechnological Potential-Phaffia australis sp. nov. and Phaffia tasmanica sp. nov

Phaffia is an orange-colored basidiomycetous yeast genus of the order Cystofilobasidiales that contains a single species, P. rhodozyma. This species is the only fungus known to produce the economically relevant carotenoid astaxanthin. Although Phaffia was originally found in the Northern hemisphere, its diversity in the southern part of the globe has been shown to be much greater. Here we analyze the genomes of two Australasian lineages that are markedly distinct from P. rhodozyma. The two divergent lineages were investigated within a comprehensive phylogenomic study of representatives of the Cystofilobasidiales that supported the recognition of two novel Phaffia species, for which we propose the names of P. australis sp. nov. and P. tasmanica sp. nov. Comparative genomics and other analyses confirmed that the two new species have the typical Phaffia hallmark-the six genes necessary for the biosynthesis of astaxanthin could be retrieved from the draft genome sequences, and this carotenoid was detected in culture extracts. In addition, the organization of the mating-type (MAT) loci is similar to that of P. rhodozyma, with synteny throughout most regions. Moreover, cases of trans-specific polymorphism involving pheromone receptor genes and pheromone precursor proteins in the three Phaffia species, together with their shared homothallism, provide additional support for their classification in a single genus.

The Evolutionary History of Ephs and Ephrins: Toward Multicellular Organisms

Eph receptor (Eph) and ephrin signaling regulate fundamental developmental processes through both forward and reverse signaling triggered upon cell–cell contact. In vertebrates, they are both classified into classes A and B, and some representatives have been identified in many metazoan groups, where their expression and functions have been well studied. We have extended previous phylogenetic analyses and examined the presence of Eph and ephrins in the tree of life to determine their origin and evolution. We have found that 1) premetazoan choanoflagellates may already have rudimental Eph/ephrin signaling as they have an Eph-/ephrin-like pair and homologs of downstream-signaling genes; 2) both forward- and reverse-downstream signaling might already occur in Porifera since sponges have most genes involved in these types of signaling; 3) the nonvertebrate metazoan Eph is a type-B receptor that can bind ephrins regardless of their membrane-anchoring structure, glycosylphosphatidylinositol, or transmembrane; 4) Eph/ephrin cross-class binding is specific to Gnathostomata; and 5) kinase-dead Eph receptors can be traced back to Gnathostomata. We conclude that Eph/ephrin signaling is of older origin than previously believed. We also examined the presence of protein domains associated with functional characteristics and the appearance and conservation of downstream-signaling pathways to understand the original and derived functions of Ephs and ephrins. We find that the evolutionary history of these gene families points to an ancestral function in cell–cell interactions that could contribute to the emergence of multicellularity and, in particular, to the required segregation of cell populations.

Mycobacterium Phage Butters-Encoded Proteins Contribute to Host Defense against Viral Attack

A diverse set of prophage-mediated mechanisms protecting bacterial hosts from infection has been recently uncovered within cluster N mycobacteriophages isolated on the host, Mycobacterium smegmatis mc2155. In that context, we unveil a novel defense mechanism in cluster N prophage Butters. By using bioinformatics analyses, phage plating efficiency experiments, microscopy, and immunoprecipitation assays, we show that Butters genes located in the central region of the genome play a key role in the defense against heterotypic viral attack. Our study suggests that a two-component system, articulated by interactions between protein products of genes 30 and 31, confers defense against heterotypic phage infection by PurpleHaze (cluster A/subcluster A3) or Alma (cluster A/subcluster A9) but is insufficient to confer defense against attack by the heterotypic phage Island3 (cluster I/subcluster I1). Therefore, based on heterotypic phage plating efficiencies on the Butters lysogen, additional prophage genes required for defense are implicated and further show specificity of prophage-encoded defense systems.IMPORTANCE Many sequenced bacterial genomes, including those of pathogenic bacteria, contain prophages. Some prophages encode defense systems that protect their bacterial host against heterotypic viral attack. Understanding the mechanisms undergirding these defense systems is crucial to appreciate the scope of bacterial immunity against viral infections and will be critical for better implementation of phage therapy that would require evasion of these defenses. Furthermore, such knowledge of prophage-encoded defense mechanisms may be useful for developing novel genetic tools for engineering phage-resistant bacteria of industrial importance.

Tailor-made multiple sequence alignments using the PRALINE 2 alignment toolkit

Summary

PRALINE 2 is a toolkit for custom multiple sequence alignment workflows. It can be used to incorporate sequence annotations, such as secondary structure or (DNA) motifs, into the alignment scoring, as well as to customize many other aspects of a progressive multiple alignment workflow.

Availability and implementation

PRALINE 2 is implemented in Python and available as open source software on GitHub: https://github.com/ibivu/PRALINE/.

Characterization of the Rifamycin-Degrading Monooxygenase From Rifamycin Producers Implicating Its Involvement in Saliniketal Biosynthesis

Rifamycin derivatives, such as rifampicin, have potent antibiotic activity and have long been used in the clinic as mainstay components for the treatment of tuberculosis, leprosy, and AIDS-associated mycobacterial infections. However, the extensive usage of these antibiotics has resulted in the rapid development of bacterial resistance. The resistance mechanisms mainly include mutations of the rifamycin target RNA polymerase of bacteria and enzymatic modifications of rifamycin antibiotics. One modification is the recently characterized rifamycin degradation catalyzed by Rox enzymes, which belong to the widely occurring flavin monooxygenases. Intriguingly, our recent sequence analysis revealed the rifamycin producers also encode Rox homologs that are not yet characterized. In this work, we expanded the study of the Rox-catalyzed rifamycin degradation. We first showed that the Rox proteins from rifamycin producers have the enzymatic rifamycin SV-degrading activity. Then we used the structurally diverse rifamycin compounds rifampicin and 16-demethylrifamycin W to probe the substrate scope and found that they each have a slightly different substrate scope. Finally, we demonstrated that Rox proteins can also catalyze the transformation of 16-demethylsalinisporamycin to 16-demethylsaliniketal A. Since 16-demethylsalinisporamycin and 16-demethylsaliniketal A are the counterpart analogs of salinisporamycin and saliniketal A, our biochemical findings not only uncover a previously uncharacterized self-resistance mechanism in the rifamycin producers, but also bridge the gap between the biosynthesis of the potential antitumor compound saliniketal A.

Characterization of the primate TRIM gene family reveals the recent evolution in primates

The tripartite motif (TRIM) gene family encodes diverse distinct proteins that play important roles in many biological processes. However, the molecular evolution and phylogenetic relationships of TRIM genes in primates are still elusive. We performed a genomic approach to identify and characterize TRIM genes in human and other six primate genomes. In total, 537 putative functional TRIM genes were identified and TRIM members varied among primates. A neighbor joining (NJ) tree based on the protein sequences of 82 human TRIM genes indicates seven TRIM groups, which is consistent with the results based on the architectural motifs. Many TRIM gene duplication events were identified, indicating a recent expansion of TRIM family in primate lineages. Interestingly, the chimpanzee genome shows the greatest TRIM gene expansion among the primates; however, its congeneric species, bonobo, has the least number of TRIM genes and no duplication event. Moreover, we identified a ~ 200 kb deletion on chromosome 11 of bonobos that results in a loss of cluster3 TRIM genes. The loss of TRIM genes might have occurred within the last 2 mys. Analysis of positive selection recovered 9 previously reported and 21 newly identified positively selected TRIM genes. In particular, most positive selected sites are located in the B30.2 domains. Our results have provided new insight into the evolution of primate TRIM genes and will broaden our understanding on the functions of the TRIM family.

Complete genome analysis of Pantoea agglomerans-infecting bacteriophage vB_PagM_AAM22

A novel myovirus, vB_PagM_AAM22 (AAM22), was isolated in Lithuania using Pantoea agglomerans as the host for phage propagation. The 49,744-bp genome of AAM22 has a G + C content of 48.4% and contains 96 probable protein-encoding genes and no genes for tRNA. In total, 34 ORFs were given a putative functional annotation, including genes associated with virion morphogenesis, DNA metabolism, and phage-host interactions. Based on comparative phylogenetic analysis, AAM22 cannot be assigned to any genus currently recognized by the ICTV and is a potential candidate to form a new genus within the family Myoviridae.

Molecular Topology of RNA Polymerase I Upstream Activation Factor

Upstream activation factor (UAF) is a multifunctional transcription factor in Saccharomyces cerevisiae that plays dual roles in activating RNA polymerase I (Pol I) transcription and repression of Pol II. For Pol I, UAF binds to a specific upstream element in the ribosomal DNA (rDNA) promoter and interacts with two other Pol I initiation factors, the TATA-binding protein (TBP) and core factor (CF). We used an integrated combination of chemical cross-linking mass spectrometry (CXMS), molecular genetics, protein biochemistry, and structural modeling to understand the topological framework responsible for UAF complex formation. Here, we report the molecular topology of the UAF complex, describe new structural and functional domains that play roles in UAF complex integrity, assembly, and biological function, and provide roles for previously identified UAF domains that include the Rrn5 SANT and histone fold domains. We highlight the role of new domains in Uaf30 that include an N-terminal winged helix domain and a disordered tethering domain as well as a BORCS6-like domain found in Rrn9. Together, our results reveal a unique network of topological features that coalesce around a histone tetramer-like core to form the dual-function UAF complex.

Molecular Identification, Characterization, and Expression Analysis of a Gonadotropin-Releasing Hormone Receptor (GnRH-R) in Pacific Abalone, Haliotis discus hannai

A full-length cDNA sequence encoding a GnRH receptor was cloned from the pleuropedal ganglion of the Pacific abalone, Haliotis discus hannai. The cloned sequence is 1499-bp in length encoding a protein of 460 amino acid residues, with a molecular mass of 52.22 kDa and an isoelectric point (pI) of 9.57. The architecture of HdhGnRH-R gene exhibited key features of G protein-coupled receptors (GPCRs), including seven membrane spanning domains, putative N-linked glycosylation motifs, and phosphorylation sites of serine and threonine residues. It shared 63%, 52%, and 30% sequence identities with Octopus vulgaris, Limulus polyphemus, and Mizuhopecten yessoensis GnRH-R II sequences, respectively. Phylogenetic analysis indicated that HdhGnRH-R gene was clustered with GnRH-R II of O. vulgaris and O. bimaculoides. qPCR assay demonstrated that the mRNA expression level of this receptor was significantly higher in the pleuropedal ganglion than that in any other examined tissue. Transcriptional activities of this gene in gonadal tissues were significantly higher in the ripening stage. The mRNA expression of this gene was significantly higher in pleuropedal ganglion, testis, and ovary at higher effective accumulative temperature (1000 °C). In situ hybridization revealed that HdhGnRH-R mRNA was expressed in neurosecretory cells of pleuropedal ganglion. Our results suggest that HdhGnRH-R gene synthesized in the neural ganglia might be involved in the control of gonadal maturation and gametogenesis of H. discus hannai. This is the first report of GnRH-R in H. discus hannai and the results may contribute to further studies of GPCRs evolution or may useful for the development of aquaculture method of this abalone species.

Identification and Functional Analysis of Temperate Siphoviridae Bacteriophages of Acinetobacter baumannii

Acinetobacter baumannii is an opportunistic pathogen that presents a serious clinical challenge due to its increasing resistance to all available antibiotics. Phage therapy has been introduced recently to treat antibiotic-incurable A. baumannii infections. In search for new A. baumannii specific bacteriophages, 20 clinical A. baumannii strains were used in two pools in an attempt to enrich phages from sewage. The enrichment resulted in induction of resident prophage(s) and three temperate bacteriophages, named vB_AbaS_fEg-Aba01, vB_AbaS_fLi-Aba02 and vB_AbaS_fLi-Aba03, all able to infect only one strain (#6597) of the 20 clinical strains, were isolated. Morphological characteristics obtained by transmission electron microscopy together with the genomic information revealed that the phages belong to the family Siphoviridae. The ca. 35 kb genomic sequences of the phages were >99% identical to each other. The linear ds DNA genomes of the phages contained 10 nt cohesive end termini, 52–54 predicted genes, an attP site and one tRNA gene each. A database search revealed an >99% identical prophage in the genome of A. baumannii strain AbPK1 (acc. no. CP024576.1). Over 99% identical prophages were also identified from two of the original 20 clinical strains (#5707 and #5920) and both were shown to be spontaneously inducible, thus very likely being the origins of the isolated phages. The phage vB_AbaS_fEg-Aba01 was also able to lysogenize the susceptible strain #6597 demonstrating that it was fully functional. The phages showed a very narrow host range infecting only two A. baumannii strains. In conclusion, we have isolated and characterized three novel temperate Siphoviridae phages that infect A.baumannii.

Structural basis for effector protein recognition by the Dot/Icm Type IVB coupling protein complex

The Legionella pneumophila Dot/Icm type IVB secretion system (T4BSS) is extremely versatile, translocating ~300 effector proteins into host cells. This specialized secretion system employs the Dot/Icm type IVB coupling protein (T4CP) complex, which includes IcmS, IcmW and LvgA, that are known to selectively assist the export of a subclass of effectors. Herein, the crystal structure of a four-subunit T4CP subcomplex bound to the effector protein VpdB reveals an interaction between LvgA and a linear motif in the C-terminus of VpdB. The same binding interface of LvgA also interacts with the C-terminal region of three additional effectors, SidH, SetA and PieA. Mutational analyses identified a FxxxLxxxK binding motif that is shared by VpdB and SidH, but not by SetA and PieA, showing that LvgA recognizes more than one type of binding motif. Together, this work provides a structural basis for how the Dot/Icm T4CP complex recognizes effectors, and highlights the multiple substrate-binding specificities of its adaptor subunit.

Intrinsic Disorder in Tetratricopeptide Repeat Proteins

Among the realm of repeat containing proteins that commonly serve as “scaffolds” promoting protein-protein interactions, there is a family of proteins containing between 2 and 20 tetratricopeptide repeats (TPRs), which are functional motifs consisting of 34 amino acids. The most distinguishing feature of TPR domains is their ability to stack continuously one upon the other, with these stacked repeats being able to affect interaction with binding partners either sequentially or in combination. It is known that many repeat-containing proteins are characterized by high levels of intrinsic disorder, and that many protein tandem repeats can be intrinsically disordered. Furthermore, it seems that TPR-containing proteins share many characteristics with hybrid proteins containing ordered domains and intrinsically disordered protein regions. However, there has not been a systematic analysis of the intrinsic disorder status of TPR proteins. To fill this gap, we analyzed 166 human TPR proteins to determine the degree to which proteins containing TPR motifs are affected by intrinsic disorder. Our analysis revealed that these proteins are characterized by different levels of intrinsic disorder and contain functional disordered regions that are utilized for protein-protein interactions and often serve as targets of various posttranslational modifications.

PL-search: a profile-link-based search method for protein remote homology detection

Protein remote homology detection is a fundamental and important task for protein structure and function analysis. Several search methods have been proposed to improve the detection performance of the remote homologues and the accuracy of ranking lists. The position-specific scoring matrix (PSSM) profile and hidden Markov model (HMM) profile can contribute to improving the performance of the state-of-the-art search methods. In this paper, we improved the profile-link (PL) information for constructing PSSM or HMM profiles, and proposed a PL-based search method (PL-search). In PL-search, more robust PLs are constructed through the double-link and iterative extending strategies, and an accurate similarity score of sequence pairs is calculated from the two-level Jaccard distance for remote homologues. We tested our method on two widely used benchmark datasets. Our results show that whether HHblits, JackHMMER or position-specific iterated-BLAST is used, PL-search obviously improves the search performance in terms of ranking quality as well as the number of detected remote homologues. For ease of use of PL-search, both its stand-alone tool and the web server are constructed, which can be accessed at http://bliulab.net/PL-search/.

The PSIPRED Protein Analysis Workbench: 20 years on

The PSIPRED Workbench is a web server offering a range of predictive methods to the bioscience community for 20 years. Here, we present the work we have completed to update the PSIPRED Protein Analysis Workbench and make it ready for the next 20 years. The main focus of our recent website upgrade work has been the acceleration of analyses in the face of increasing protein sequence database size. We additionally discuss any new software, the new hardware infrastructure, our webservices and web site. Lastly we survey updates to some of the key predictive algorithms available through our website.

Pantoeaagglomerans-Infecting Bacteriophage vB_PagS_AAS21: A Cold-Adapted Virus Representing a Novel Genus within the Family Siphoviridae

A novel cold-adapted siphovirus, vB_PagS_AAS21 (AAS21), was isolated in Lithuania using Pantoea agglomerans as the host for phage propagation. AAS21 has an isometric head (~85 nm in diameter) and a non-contractile flexible tail (~174 × 10 nm). With a genome size of 116,649 bp, bacteriophage AAS21 is the largest Pantoea-infecting siphovirus sequenced to date. The genome of AAS21 has a G+C content of 39.0% and contains 213 putative protein-encoding genes and 29 genes for tRNAs. A comparative sequence analysis revealed that 89 AAS21 open reading frames (ORFs) code for unique proteins that have no reliable identity to database entries. In total, 63 AAS21 ORFs were functionally annotated, including those coding for the proteins responsible for virion morphogenesis, phage-host interactions, and DNA metabolism. Proteomic analysis led to the experimental identification of 19 virion proteins, including 11 that were predicted by bioinformatics approaches. Based on comparative phylogenetic analysis, AAS21 cannot be assigned to any genus currently recognized by ICTV and may represents a new branch of viruses within the family Siphoviridae.

Expansion of the Transporter-Opsin-G protein-coupled receptor superfamily with five new protein families

Here we provide bioinformatic evidence that the Organo-Arsenical Exporter (ArsP), Endoplasmic Reticulum Retention Receptor (KDELR), Mitochondrial Pyruvate Carrier (MPC), L-Alanine Exporter (AlaE), and the Lipid-linked Sugar Translocase (LST) protein families are members of the Transporter-Opsin-G Protein-coupled Receptor (TOG) Superfamily. These families share domains homologous to well-established TOG superfamily members, and their topologies of transmembranal segments (TMSs) are compatible with the basic 4-TMS repeat unit characteristic of this Superfamily. These repeat units tend to occur twice in proteins as a result of intragenic duplication events, often with subsequent gain/loss of TMSs in many superfamily members. Transporters within the ArsP family allow microbial pathogens to expel toxic arsenic compounds from the cell. Members of the KDELR family are involved in the selective retrieval of proteins that reside in the endoplasmic reticulum. Proteins of the MPC family are involved in the transport of pyruvate into mitochondria, providing the organelle with a major oxidative fuel. Members of family AlaE excrete L-alanine from the cell. Members of the LST family are involved in the translocation of lipid-linked glucose across the membrane. These five families substantially expand the range of substrates of transport carriers in the superfamily, although KDEL receptors have no known transport function. Clustering of protein sequences reveals the relationships among families, and the resulting tree correlates well with the degrees of sequence similarity documented between families. The analyses and programs developed to detect distant relatedness, provide insights into the structural, functional, and evolutionary relationships that exist between families of the TOG superfamily, and should be of value to many other investigators.

Isolation and Characterization of a Novel Phage for Controlling Multidrug-Resistant Klebsiella pneumoniae

The emergence of multidrug-resistant bacterial pathogens has severely threatened global health. A phage with the ability to efficiently and specifically lyse bacteria is considered an alternative for controlling multidrug-resistant bacterial pathogens. The discovery of novel agents for controlling the infections caused by K. pneumoniae is urgent due to the broad multidrug-resistance of K. pneumoniae. Only a few phage isolates have been reported to infect multidrug-resistant K. pneumoniae. In this study, by using the multidrug-resistant K. pneumoniae strain as an indicator, a novel phage called vB_KleS-HSE3, which maintains high antibacterial activity and high physical stability, was isolated from hospital sewage. This phage infected one of four tested multidrug-resistant K. pneumoniae strains. This phage belongs to the Siphoviridae family and a comparative genomic analysis showed that this phage is part of a novel phage lineage among the Siphoviridae family of phages that infect strains of Klebsiella. Based on its features, the vB_KleS-HSE3 phage has potential for controlling infections caused by multidrug-resistant K. pneumoniae.

Molecular characterization and spatiotemporal expression of prohormone convertase 2 in the Pacific abalone, Haliotis discus hannai

Prohormone convertases (PCs) are subtilisin-like proteases responsible for the intracellular processing of prohormones and proneuropeptides in vertebrates and invertebrates. The full-length PC2 cDNA sequence was cloned from pleuropedal ganglion of Haliotis discus hannai, consisted of 2254-bp with an open reading frame of 1989-bp and encoded a protein of 662 amino acid residues. The architecture of Hdh PC2 displayed key features of PCs, including a signal peptide, a pro-segment domain with sites for autocatalytic activation, a catalytic domain, and a pro-protein domain (P-domain). It shares the highest homology of its amino acid sequence with the PC2 from H. asinina and to lesser extent with that of Homo sapiens and Rana catesbeiana PC2. Sequence alignment analysis indicated that Hdh PC2 was highly conserved in the catalytic domain, including a catalytic triad of serine proteinases of the subtilisin family at positions Asp-195, His-236, and Ser-412. The cloned sequence contained a canonical integrin binding sequence, and four cysteine residues involved in the formation of an intramolecular disulfide link. Phylogenetic analysis revealed that the Hdh PC2 is robustly clustered with the Has PC2. Quantitative PCR assay demonstrated that the Hdh PC2 was predominantly expressed in the pleuropedal ganglion rather than in other examined tissues. Although PC2 mRNA was expressed throughout the gametogenetic cycle of male and female abalone, the expression level was significantly higher in the ripening stage of female abalone. Also, a significantly higher expression was observed in the pleuropedal ganglion and gonadal tissues at a higher effective accumulative temperature (1000°C). In situ hybridization revealed that the PC2 mRNA expressing neurosecretory cells were distributed in the cortex region of the pleuropedal ganglion. According to the results, it can be concluded that pleuropedal ganglion is the highest site of PC2 activity, and this enzyme might be involved in the abalone reproduction process.

The use of local structural similarity of distant homologues for crystallographic model building from a molecular-replacement solution

The performance of automated protein model building usually decreases with resolution, mainly owing to the lower information content of the experimental data. This calls for a more elaborate use of the available structural information about macromolecules. Here, a new method is presented that uses structural homologues to improve the quality of protein models automatically constructed using ARP/wARP. The method uses local structural similarity between deposited models and the model being built, and results in longer main-chain fragments that in turn can be more reliably docked to the protein sequence. The application of the homology-based model extension method to the example of a CFA synthase at 2.7 Å resolution resulted in a more complete model with almost all of the residues correctly built and docked to the sequence. The method was also evaluated on 1493 molecular-replacement solutions at a resolution of 4.0 Å and better that were submitted to the ARP/wARP web service for model building. A significant improvement in the completeness and sequence coverage of the built models has been observed.

Functional Innovation in the Evolution of the Calcium-Dependent System of the Eukaryotic Endoplasmic Reticulum

The origin of eukaryotes was marked by the emergence of several novel subcellular systems. One such is the calcium (Ca²⁺)-stores system of the endoplasmic reticulum, which profoundly influences diverse aspects of cellular function including signal transduction, motility, division, and biomineralization. We use comparative genomics and sensitive sequence and structure analyses to investigate the evolution of this system. Our findings reconstruct the core form of the Ca²⁺-stores system in the last eukaryotic common ancestor as having at least 15 proteins that constituted a basic system for facilitating both Ca²⁺ flux across endomembranes and Ca²⁺-dependent signaling. We present evidence that the key EF-hand Ca²⁺-binding components had their origins in a likely bacterial symbiont other than the mitochondrial progenitor, whereas the protein phosphatase subunit of the ancestral calcineurin complex was likely inherited from the asgard archaeal progenitor of the stem eukaryote. This further points to the potential origin of the eukaryotes in a Ca²⁺-rich biomineralized environment such as stromatolites. We further show that throughout eukaryotic evolution there were several acquisitions from bacteria of key components of the Ca²⁺-stores system, even though no prokaryotic lineage possesses a comparable system. Further, using quantitative measures derived from comparative genomics we show that there were several rounds of lineage-specific gene expansions, innovations of novel gene families, and gene losses correlated with biological innovation such as the biomineralized molluscan shells, coccolithophores, and animal motility. The burst of innovation of new genes in animals included the wolframin protein associated with Wolfram syndrome in humans. We show for the first time that it contains previously unidentified Sel1, EF-hand, and OB-fold domains, which might have key roles in its biochemistry.

Regulation of Phosphoribosyl-Linked Serine Ubiquitination by Deubiquitinases DupA and DupB

The family of bacterial SidE enzymes catalyzes non-canonical phosphoribosyl-linked (PR) serine ubiquitination and promotes infectivity of Legionella pneumophila. Here, we describe identification of two bacterial effectors that reverse PR ubiquitination and are thus named deubiquitinases for PR ubiquitination (DUPs; DupA and DupB). Structural analyses revealed that DupA and SidE ubiquitin ligases harbor a highly homologous catalytic phosphodiesterase (PDE) domain. However, unlike SidE ubiquitin ligases, DupA displays increased affinity to PR-ubiquitinated substrates, which allows DupA to cleave PR ubiquitin from substrates. Interfering with DupA-ubiquitin binding switches its activity toward SidE-type ligase. Given the high affinity of DupA to PR-ubiquitinated substrates, we exploited a catalytically inactive DupA mutant to trap and identify more than 180 PR-ubiquitinated host proteins in Legionella-infected cells. Proteins involved in endoplasmic reticulum (ER) fragmentation and membrane recruitment to Legionella-containing vacuoles (LCV) emerged as major SidE targets. The global map of PR-ubiquitinated substrates provides critical insights into host-pathogen interactions during Legionella infection.

Attempts to develop an enzyme converting DHIV to KIV

Dihydroxy-acid dehydratase (DHAD) catalyzes the dehydration of R-2,3-dihydroxyisovalerate (DHIV) to 2-ketoisovalerate (KIV) using an Fe-S cluster as a cofactor, which is sensitive to oxidation and expensive to synthesize. In contrast, sugar acid dehydratases catalyze the same chemical reactions using a magnesium ion. Here, we attempted to substitute the high-cost DHAD with a cost-efficient engineered sugar acid dehydratase using computational protein design (CPD). First, we tried without success to modify the binding pocket of a sugar acid dehydratase to accommodate the smaller, more hydrophobic DHIV. Then, we used a chemically activated substrate analog to react with sugar acid dehydratases or other enolase superfamily enzymes. Mandelate racemase from Pseudomonas putida (PpManR) and the putative sugar acid dehydratase from Salmonella typhimurium (StPutD) showed beta-elimination activity towards chlorolactate (CLD). CPD combined with medium-throughput selection improved the PpManR kcat/KM for CLD by four-fold. However, these enzyme variants did not show dehydration activity towards DHIV. Lastly, assuming phosphorylation could also be a good activation mechanism, we found that mevalonate-3-kinase (M3K) from Picrophilus torridus (PtM3K) exhibited adenosine triphosphate (ATP) hydrolysis activity when mixed with DHIV, indicating phosphorylation activity towards DHIV. Engineering PpManR or StPutD to accept 3-phospho-DHIV as a substrate was performed, but no variants with the desired activity were obtained.

Variability of mitochondrial ORFans hints at possible differences in the system of doubly uniparental inheritance of mitochondria among families of freshwater mussels (Bivalvia: Unionida)

BACKGROUND

Supernumerary ORFan genes (i.e., open reading frames without obvious homology to other genes) are present in the mitochondrial genomes of gonochoric freshwater mussels (Bivalvia: Unionida) showing doubly uniparental inheritance (DUI) of mitochondria. DUI is a system in which distinct female-transmitted and male-transmitted mitotypes coexist in a single species. In families Unionidae and Margaritiferidae, the transition from dioecy to hermaphroditism and the loss of DUI appear to be linked, and this event seems to affect the integrity of the ORFan genes. These observations led to the hypothesis that the ORFans have a role in DUI and/or sex determination. Complete mitochondrial genome sequences are however scarce for most families of freshwater mussels, therefore hindering a clear localization of DUI in the various lineages and a comprehensive understanding of the influence of the ORFans on DUI and sexual systems. Therefore, we sequenced and characterized eleven new mitogenomes from poorly sampled freshwater mussel families to gather information on the evolution and variability of the ORFan genes and their protein products.

RESULTS

We obtained ten complete plus one almost complete mitogenome sequence from ten representative species (gonochoric and hermaphroditic) of families Margaritiferidae, Hyriidae, Mulleriidae, and Iridinidae. ORFan genes are present only in DUI species from Margaritiferidae and Hyriidae, while non-DUI species from Hyriidae, Iridinidae, and Mulleriidae lack them completely, independently of their sexual system. Comparisons among the proteins translated from the newly characterized ORFans and already known ones provide evidence of conserved structures, as well as family-specific features.

CONCLUSIONS

The ORFan proteins show a comparable organization of secondary structures among different families of freshwater mussels, which supports a conserved physiological role, but also have distinctive family-specific features. Given this latter observation and the fact that the ORFans can be either highly mutated or completely absent in species that secondarily lost DUI depending on their respective family, we hypothesize that some aspects of the connection among ORFans, sexual systems, and DUI may differ in the various lineages of unionids.

Reconstructing the Remote Origins of a Fold Singleton from a Flavodoxin-Like Ancestor

Evolutionary processes that led to the emergence of structured protein domains left footprints in the sequences of modern proteins. We searched for such hints employing state-of-the-art sequence analysis and found evidence that the HemD-like fold emerged from the flavodoxin-like fold through segment swap and gene duplication. To verify this hypothesis, we reverted these evolutionary steps experimentally, constructing a HemD-half that resulted in a protein with the canonical flavodoxin-like architecture. These results of fold reconstruction from the sequence of a different fold strongly support our hypothesis of common ancestry. It further illustrates the plasticity of modern proteins to form new folded proteins.

Genome-wide analysis of Homo sapiens, Arabidopsis thaliana, and Saccharomyces cerevisiae reveals novel attributes of tail-anchored membrane proteins

BACKGROUND

Tail-anchored membrane proteins (TAMPs) differ from other integral membrane proteins, because they contain a single transmembrane domain at the extreme carboxyl-terminus and are therefore obliged to target to membranes post-translationally. Although 3-5% of all transmembrane proteins are predicted to be TAMPs only a small number are well characterized.

RESULTS

To identify novel putative TAMPs across different species, we used TAMPfinder software to identify 859, 657 and 119 putative TAMPs in human (Homo sapiens), plant (Arabidopsis thaliana), and yeast (Saccharomyces cerevisiae), respectively. Bioinformatics analyses of these putative TAMP sequences suggest that the list is highly enriched for authentic TAMPs. To experimentally validate the software predictions several human and plant proteins identified by TAMPfinder that were previously uncharacterized were expressed in cells and visualized at subcellular membranes by fluorescence microscopy and further analyzed by carbonate extraction or by bimolecular fluorescence complementation. With the exception of the pro-apoptotic protein harakiri, which is, peripherally bound to the membrane this subset of novel proteins behave like genuine TAMPs. Comprehensive bioinformatics analysis of the generated TAMP datasets revealed previously unappreciated common and species-specific features such as the unusual size distribution of and the propensity of TAMP proteins to be part of larger complexes. Additionally, novel features of the amino acid sequences that anchor TAMPs to membranes were also revealed.

CONCLUSIONS

The findings in this study more than double the number of predicted annotated TAMPs and provide new insights into the common and species-specific features of TAMPs. Furthermore, the list of TAMPs and annotations provide a resource for further investigation.

Identification, characterization, and expression analysis of a serotonin receptor involved in the reproductive process of the Pacific abalone, Haliotis discus hannai

Serotonin receptor (5-HT) is a biogenic amine acting as a neurotransmitter and neuromodulator that mediates various aspects of reproduction and gametogenesis. The full-length nucleotide sequence of Haliotis discus hannai encodes a protein of 417 amino acids with a predicted molecular mass of 46.54 kDa and isoelectric point of 8.94. The structural profile of 5-HT_Hdh displayed key features of G protein-coupled receptors, including seven hydrophobic transmembrane domains, putative N-linked glycosylation sites, and several phosphorylation consensus motifs. It shares the highest homology of its amino acid sequence with the 5-HT receptor from Haliotis asinina, and to lesser extent of human 5-HT receptor. The cloned sequence possesses two cysteine residues (Cys-115 and Cys-193), which are likely to form a disulfide bond. Phylogenetic comparison with other known 5-HT receptor genes revealed that the 5-HT_Hdh is most closely related to the 5-HT_Ha receptor. The three-dimensional structure of the 5-HT_Hdh showed multiple alpha helices which is separated by a helix-loop-helix (HLH) structure. Quantitative PCR demonstrated that the receptor mRNA was predominantly expressed in the pleuropedal ganglion. Significant differences in the transcriptional activity of the 5-HT_Hdh gene were observed in the ovary at the ripening stage. An exclusive expression was detected in pleuropedal ganglion, testis, and ovary at higher effective accumulative temperature (1000 °C). In situ hybridization showed that the 5-HT_Hdh expressing neurosecretory cells were distributed in the cortex of the pleuropedal ganglion. Our results suggest that 5-HT_Hdh synthesized in the neural ganglia may be involved in oocyte maturation and spawning of H. discus hannai.

The ancestral KH peptide at the root of a domain family with three different folds

Motivation

The direct ancestor of the DNA-protein world of today is considered to have been an RNA-peptide world, in which peptides were co-factors of RNA-mediated catalysis and replication. Evidence for these ancestral peptides, from which folded proteins evolved, can be derived even today from regions of local sequence similarity within globally dissimilar folds. One of these is the 45-residue motif common to both folds of the hnRNP K homology (KH) domain.

Results

In a survey of KH domains, we found a third fold that contains the KH motif at its core. This corresponds to the Small Domain of bacterial Ribonucleases G/E and, like type I and type II KH domains, it cannot be related to the others by a single genetic event, providing further support for the KH motif as an ancestral peptide predating folded proteins.

Supplementary information

Supplementary data are available at Bioinformatics online.

The WW1 Domain Enhances Autoinhibition in Smurf Ubiquitin Ligases

Downregulation of ubiquitin (Ub) ligase activity prevents premature ubiquitination and is critical for cellular homeostasis. Nedd4 Ub ligases share a common domain architecture and yet are regulated in distinct ways through interactions of the catalytic HECT domain with the N-terminal C2 domain or the central WW domain region. Smurf1 and Smurf2 are two highly related Nedd4 ligases with ~70% overall sequence identity. Here, we show that the Smurf1 C2 domain interacts with the HECT domain and inhibits ligase activity in trans. However, in contrast to Smurf2, we find that full-length Smurf1 is a highly active Ub ligase, and we can attribute this striking difference in regulation to the lack of one WW domain (WW1) in Smurf1. Using NMR spectroscopy and biochemical assays, we identified the WW1 region as an additional inhibitory element in Smurf2 that cooperates with the C2 domain to enhance HECT domain binding and Smurf2 inhibition. Our work provides important insights into Smurf regulation and highlights that the activities of highly related proteins can be controlled in distinct ways.

Analysis of Heterodimeric "Mutual Synergistic Folding"-Complexes

Several intrinsically disordered proteins (IDPs) are capable to adopt stable structures without interacting with a folded partner. When the folding of all interacting partners happens at the same time, coupled with the interaction in a synergistic manner, the process is called Mutual Synergistic Folding (MSF). These complexes represent a discrete subset of IDPs. Recently, we collected information on their complexes and created the MFIB (Mutual Folding Induced by Binding) database. In a previous study, we compared homodimeric MSF complexes with homodimeric and monomeric globular proteins with similar amino acid sequence lengths. We concluded that MSF homodimers, compared to globular homodimeric proteins, have a greater solvent accessible main-chain surface area on the contact surface of the subunits, which becomes buried during dimerization. The main driving force of the folding is the mutual shielding of the water-accessible backbones, but the formation of further intermolecular interactions can also be relevant. In this paper, we will report analyses of heterodimeric MSF complexes. Our results indicate that the amino acid composition of the heterodimeric MSF monomer subunits slightly diverges from globular monomer proteins, while after dimerization, the amino acid composition of the overall MSF complexes becomes more similar to overall amino acid compositions of globular complexes. We found that inter-subunit interactions are strengthened, and additionally to the shielding of the solvent accessible backbone, other factors might play an important role in the stabilization of the heterodimeric structures, likewise energy gain resulting from the interaction of the two subunits with different amino acid compositions. We suggest that the shielding of the β-sheet backbones and the formation of a buried structural core along with the general strengthening of inter-subunit interactions together could be the driving forces of MSF protein structural ordering upon dimerization.

Jwalk and MNXL web server: model validation using restraints from crosslinking mass spectrometry

Motivation

Crosslinking Mass Spectrometry generates restraints that can be used to model proteins and protein complexes. Previously, we have developed two methods, to help users achieve better modelling performance from their crosslinking restraints: Jwalk, to estimate solvent accessible distances between crosslinked residues and MNXL, to assess the quality of the models based on these distances.

Results

Here, we present the Jwalk and MNXL webservers, which streamline the process of validating monomeric protein models using restraints from crosslinks. We demonstrate this by using the MNXL server to filter models made of varying quality, selecting the most native-like.

Availability and implementation

The webserver and source code are freely available from jwalk.ismb.lon.ac.uk and mnxl.ismb.lon.ac.uk.

Supplementary information

Supplementary data are available at Bioinformatics online.

Current Progress on Equilibrative Nucleoside Transporter Function and Inhibitor Design

Physiological nucleosides are used for the synthesis of DNA, RNA, and ATP in the cell and serve as universal mammalian signaling molecules that regulate physiological processes such as vasodilation and platelet aggregation by engaging with cell surface receptors. The same pathways that allow uptake of physiological nucleosides mediate the cellular import of synthetic nucleoside analogs used against cancer, HIV, and other viral diseases. Physiological nucleosides and nucleoside drugs are imported by two families of nucleoside transporters: the SLC28 concentrative nucleoside transporters (CNTs) and SLC29 equilibrative nucleoside transporters (ENTs). The four human ENT paralogs are expressed in distinct tissues, localize to different subcellular sites, and transport a variety of different molecules. Here we provide an overview of the known structure-function relationships of the ENT family with a focus on ligand binding and transport in the context of a new hENT1 homology model. We provide a generic residue numbering system for the different ENTs to facilitate the interpretation of mutational data produced using different ENT homologs. The discovery of paralog-selective small-molecule modulators is highly relevant for the design of new therapies and for uncovering the functions of poorly characterized ENT family members. Here, we discuss recent developments in the discovery of new paralog-selective small-molecule ENT inhibitors, including new natural product-inspired compounds. Recent progress in the ability to heterologously produce functional ENTs will allow us to gain insight into the structure and functions of different ENT family members as well as the rational discovery of highly selective inhibitors.

Characterization of intellectual disability and autism comorbidity through gene panel sequencing

Intellectual disability (ID) and autism spectrum disorder (ASD) are clinically and genetically heterogeneous diseases. Recent whole exome sequencing studies indicated that genes associated with different neurological diseases are shared across disorders and converge on common functional pathways. Using the Ion Torrent platform, we developed a low-cost next-generation sequencing gene panel that has been transferred into clinical practice, replacing single disease-gene analyses for the early diagnosis of individuals with ID/ASD. The gene panel was designed using an innovative in silico approach based on disease networks and mining data from public resources to score disease-gene associations. We analyzed 150 unrelated individuals with ID and/or ASD and a confident diagnosis has been reached in 26 cases (17%). Likely pathogenic mutations have been identified in another 15 patients, reaching a total diagnostic yield of 27%. Our data also support the pathogenic role of genes recently proposed to be involved in ASD. Although many of the identified variants need further investigation to be considered disease-causing, our results indicate the efficiency of the targeted gene panel on the identification of novel and rare variants in patients with ID and ASD.

Multiple C2 domains and transmembrane region proteins (MCTPs) tether membranes at plasmodesmata

In eukaryotes, membrane contact sites (MCS) allow direct communication between organelles. Plants have evolved a unique type of MCS, inside intercellular pores, the plasmodesmata, where endoplasmic reticulum (ER)-plasma membrane (PM) contacts coincide with regulation of cell-to-cell signalling. The molecular mechanism and function of membrane tethering within plasmodesmata remain unknown. Here, we show that the multiple C2 domains and transmembrane region protein (MCTP) family, key regulators of cell-to-cell signalling in plants, act as ER-PM tethers specifically at plasmodesmata. We report that MCTPs are plasmodesmata proteins that insert into the ER via their transmembrane region while their C2 domains dock to the PM through interaction with anionic phospholipids. A Atmctp3/Atmctp4 loss of function mutant induces plant developmental defects, impaired plasmodesmata function and composition, while MCTP4 expression in a yeast Δtether mutant partially restores ER-PM tethering. Our data suggest that MCTPs are unique membrane tethers controlling both ER-PM contacts and cell-to-cell signalling.

Rescrutiny of the sansanmycin biosynthetic gene cluster leads to the discovery of a novel sansanmycin analogue with more potency against Mycobacterium tuberculosis

A novel sansanmycin analogue, sansanmycin Q (1), was identified by genome mining from the fermentation broth of Streptomyces sp. SS (CPCC 200442). In comparison with other sansanmycin compounds, sansanmycin Q has an extra glycine residue at the N-terminus of the pseudopeptide backbone. The additional glycine was proved to be assembled to sansanmycin A by SsaB, a tRNA-dependent aminoacyltransferase, based on the results of rescrutiny of sansanmycin biosynthetic gene cluster, and then overexpression and knockout of ssaB in the wild-type strain. The structure of sansanmycin Q was assigned by interpretation of NMR and mass spectral data. The results of the bioassay disclosed that sansanmycin Q exhibited more potency against Mycobacterium tuberculosis H₃₇Rv and a rifampicin- and isoniazid-resistant strain than sansanmycin A.

Cryptic inoviruses revealed as pervasive in bacteria and archaea across Earth's biomes

Bacteriophages from the Inoviridae family (inoviruses) are characterized by their unique morphology, genome content and infection cycle. One of the most striking features of inoviruses is their ability to establish a chronic infection whereby the viral genome resides within the cell in either an exclusively episomal state or integrated into the host chromosome and virions are continuously released without killing the host. To date, a relatively small number of inovirus isolates have been extensively studied, either for biotechnological applications, such as phage display, or because of their effect on the toxicity of known bacterial pathogens including Vibrio cholerae and Neisseria meningitidis. Here, we show that the current 56 members of the Inoviridae family represent a minute fraction of a highly diverse group of inoviruses. Using a machine learning approach leveraging a combination of marker gene and genome features, we identified 10,295 inovirus-like sequences from microbial genomes and metagenomes. Collectively, our results call for reclassification of the current Inoviridae family into a viral order including six distinct proposed families associated with nearly all bacterial phyla across virtually every ecosystem. Putative inoviruses were also detected in several archaeal genomes, suggesting that, collectively, members of this supergroup infect hosts across the domains Bacteria and Archaea. Finally, we identified an expansive diversity of inovirus-encoded toxin–antitoxin and gene expression modulation systems, alongside evidence of both synergistic (CRISPR evasion) and antagonistic (superinfection exclusion) interactions with co-infecting viruses, which we experimentally validated in a Pseudomonas model. Capturing this previously obscured component of the global virosphere may spark new avenues for microbial manipulation approaches and innovative biotechnological applications.

A machine learning approach was used to recover over 10,000 inovirus-like sequences from existing microbial genomes and metagenomes, consequently proposing the reclassification of the Inoviridae family to a viral order, and uncover the previously unrecognized diversity of these viruses across hosts and environments.

Structure and Function of the Branched Receptor-Binding Complex of Bacteriophage CBA120

Bacteriophages recognize their host cells with the help of tail fiber and tailspike proteins that bind, cleave, or modify certain structures on the cell surface. The spectrum of ligands to which the tail fibers and tailspikes can bind is the primary determinant of the host range. Bacteriophages with multiple tailspike/tail fibers are thought to have a wider host range than their less endowed relatives but the function of these proteins remains poorly understood. Here, we describe the structure, function, and substrate specificity of three tailspike proteins of bacteriophage CBA120-TSP2, TSP3 and TSP4 (orf211 through orf213, respectively). We show that tailspikes TSP2, TSP3 and TSP4 are hydrolases that digest the O157, O77, and O78 Escherichia coli O-antigens, respectively. We demonstrate that recognition of the E. coli O157:H7 host by CBA120 involves binding to and digesting the O157 O-antigen by TSP2. We report the crystal structure of TSP2 in complex with a repeating unit of the O157 O-antigen. We demonstrate that according to the specificity of its tailspikes TSP2, TSP3, and TSP4, CBA120 can infect E. coli O157, O77, and O78, respectively. We also show that CBA120 infects Salmonella enterica serovar Minnesota, and this host range expansion is likely due to the function of TSP1. Finally, we describe the assembly pathway and the architecture of the TSP1-TSP2-TSP3-TSP4 branched complex in CBA120 and its related ViI-like phages.

Genomic differences within the phylum Marinimicrobia: From waters to sediments in the Mariana Trench

Marinimicrobia are widespread from the marine surface to the hadal zone. Major clades of Marinimicrobia have evolved to different ecotypes along with energy gradients, but their genomes in deeper waters and sediments have rarely been studied. Here we obtained 11 Marinimicrobia draft genomes from the water column in the full-ocean depth and the hadal sediments in the Mariana Trench. All the predicted genomic capabilities of the metagenome-assembled genomes (MAGs) are indicative of heterotrophic lifestyle. The MAGs from the hadal depths are distinct from those from the mesopelagic and bathypelagic depths by enrichment of the genes involved in amino acids metabolism and mismatch repair. Compared with the MAGs from waters, those from the sediments were dramatically expanded by acquiring the genes responsible for chemotaxis, mobility and the two-component systems. Marinimicrobia were apparently differentiated in the environments with different depths, organic matters and electronic acceptors. Our results also posit a potential evolutionary relationship between the species inhabiting the waters and sediments, indicating the occurrence of allopatric speciation in Marinimicrobia.

Coupling of ATPase activity, microtubule binding, and mechanics in the dynein motor domain

The movement of a molecular motor protein along a cytoskeletal track requires communication between enzymatic, polymer-binding, and mechanical elements. Such communication is particularly complex and not well understood in the dynein motor, an ATPase that is comprised of a ring of six AAA domains, a large mechanical element (linker) spanning over the ring, and a microtubule-binding domain (MTBD) that is separated from the AAA ring by a ~ 135 Å coiled-coil stalk. We identified mutations in the stalk that disrupt directional motion, have microtubule-independent hyperactive ATPase activity, and nucleotide-independent low affinity for microtubules. Cryo-electron microscopy structures of a mutant that uncouples ATPase activity from directional movement reveal that nucleotide-dependent conformational changes occur normally in one-half of the AAA ring, but are disrupted in the other half. The large-scale linker conformational change observed in the wild-type protein is also inhibited, revealing that this conformational change is not required for ATP hydrolysis. These results demonstrate an essential role of the stalk in regulating motor activity and coupling conformational changes across the two halves of the AAA ring.

Biological properties and genomics analysis of vB_KpnS_GH-K3, a Klebsiella phage with a putative depolymerase-like protein

Bacteriophages have been recently revisited as an alternative biocontrol tool due to the limitations of antibiotic treatment. In this study, we reported on the biological characteristics and genomic information of vB_KpnS_GH-K3 (abbreviated as GH-K3), a Klebsiella phage of the Siphoviridae family, which was previously isolated from a hospital sewage system. One-step growth curve analysis indicated that the burst size of GH-K3 was 291 PFU/cell. GH-K3 maintained a stable titer in a broad range of pH values (6-10) and temperature (up to 50 °C). Based on bioinformatics analysis, GH-K3 comprises of 49,427 bp containing a total of 77 open reading frames (ORFs), which share high degree of nucleotide similarity and close evolutionary relationships with at least 12 other Klebsiella phages. Of note, GH-K3 gp32 was identified as a unique ORF. The major segment of gp32 sequence at the C-terminus (residues 351-907) was found highly variable as determined by its mismatch with the nucleotide and protein sequences available at NCBI database. Furthermore, HHpred analysis indicated that GH-K3 gp32 contains three domains (PDB ID: 5W6S_A, 3GQ8_A and 1BHE_A) similar to depolymerase (depoKP36) of Klebsiella phage KP36 suggestive of a potential depolymerase activity during host receptor-binding in the processes of phage infection. Altogether, the current data revealed a novel putative depolymerase-like protein which is most likely to play an important role in phage-host interaction.

Novel mutations in KMT2B offer pathophysiological insights into childhood-onset progressive dystonia

Rapid progress has recently been made in the elucidation of the genetic basis of childhood-onset inherited generalized dystonia (IGD) due to the implementation of genomic sequencing methodologies. We identified four patients with childhood-onset IGD harboring novel disease-causing mutations in lysine-specific histone methyltransferase 2B gene (KMT2B) by whole-exome sequencing. The main focus of this paper is to gain novel pathophysiological insights through understanding the molecular consequences of these mutations. The disease course is mostly progressive, evolving from lower limbs into generalized dystonia, which could be associated with dysarthria, dysphonia, intellectual disability, orofacial dyskinesia, and sometimes distinct dysmorphic facial features. In two patients, motor performances improved after bilateral implantation of deep brain stimulation in the globus pallidus internus (GPi-DBS). Pharmacotherapy with trihexyphenidyl reduced dystonia in two patients. We discovered three novel KMT2B mutations. Our analyses revealed that the mutation in patient 1 (c.7463A > G, p.Y2488C) is localized in the highly conserved FYRC domain of KMT2B. This mutation holds the potential to alter the inter-domain FYR interactions, which could lead to KMT2B instability. The mutations in patients 2 and 3 (c.3596_3697insC, p.M1202Dfs*22; c.4229delA, p.Q1410Rfs*12) lead to predicted unstable transcripts, likely to be subject to degradation by non-sense-mediated decay. Childhood-onset progressive dystonia with orofacial involvement is one of the main clinical manifestations of KMT2B mutations. In all, 26% (18/69) of the reported cases have T2 signal alterations of the globus pallidus internus, mostly at a younger age. Anticholinergic medication and GPi-DBS are promising treatment options and shall be considered early.

A Three-Dimensional Model of Human Lysyl Oxidase, a Cross-Linking Enzyme

Lysyl oxidase (LOX) is a cross-linking enzyme identified 50 years ago, but its 3D structure is still unknown. We have thus built a 3D model of human LOX by homology modeling using the X-ray structure of human lysyl oxidase-like 2 as a template. This model is the first one to recapitulate all known biochemical features of LOX, namely, the coordination of the copper ion and the formation of the lysine tyrosylquinone cofactor and the disulfide bridges. Furthermore, this model is stable during a 1 μs molecular dynamics simulation. The catalytic site is located in a groove surrounded by two loops. The distance between these loops fluctuated during the simulations, which suggests that the groove forms a hinge with a variable opening, which is able to accommodate the various sizes of LOX substrates. This 3D model is a pre-requisite to perform docking experiments with LOX substrates and other partners to identify binding sites and to design new LOX inhibitors specific for therapeutic purpose.

An efficient use of X-ray information, homology modeling, molecular dynamics and knowledge-based docking techniques to predict protein-monosaccharide complexes

Unraveling the structure of lectin-carbohydrate complexes is vital for understanding key biological recognition processes and development of glycomimetic drugs. Molecular Docking application to predict them is challenging due to their low affinity, hydrophilic nature and ligand conformational diversity. In the last decade several strategies, such as the inclusion of glycan conformation specific scoring functions or our developed solvent-site biased method, have improved carbohydrate docking performance but significant challenges remain, in particular, those related to receptor conformational diversity. In the present work we have analyzed conventional and solvent-site biased autodock4 performance concerning receptor conformational diversity as derived from different crystal structures (apo and holo), Molecular Dynamics snapshots and Homology-based models, for 14 different lectin-monosaccharide complexes. Our results show that both conventional and biased docking yield accurate lectin-monosaccharide complexes, starting from either apo or homology-based structures, even when only moderate (45%) sequence identity templates are available. An essential element for success is a proper combination of a middle-sized (10-100 structures) conformational ensemble, derived either from Molecular dynamics or multiple homology model building. Consistent with our previous works, results show that solvent-site biased methods improve overall performance, but that results are still highly system dependent. Finally, our results also show that docking can select the correct receptor structure within the ensemble, underscoring the relevance of joint evaluation of both ligand pose and receptor conformation.

The complete genome of Rachiplusia nu nucleopolyhedrovirus (RanuNPV) and the identification of a baculoviral CPD-photolyase homolog

We described a novel baculovirus isolated from the polyphagous insect pest Rachiplusia nu. The virus presented pyramidal-shaped occlusion bodies (OBs) with singly-embed nucleocapsids and a dose mortality response of 6.9 × 10³ OBs/ml to third-instar larvae of R. nu. The virus genome is 128,587 bp long with a G + C content of 37.9% and 134 predicted ORFs. The virus is an alphabaculovirus closely related to Trichoplusia ni single nucleopolyhedrovirus, Chrysodeixis chalcites nucleopolyhedrovirus, and Chrysodeixis includens single nucleopolyhedrovirus and may constitute a new species. Surprisingly, we found co-evolution among the related viruses and their hosts at species level. Besides, auxiliary genes with homologs in other baculoviruses were found, e.g. a CPD-photolyase. The gene seemed to be result of a single event of horizontal transfer from lepidopterans to alphabaculovirus, followed by a transference from alpha to betabaculovirus. The predicted protein appears to be an active enzyme that ensures likely DNA protection from sunlight.

Exploring lectin-like activity of the S-layer protein of Lactobacillus acidophilus ATCC 4356

The surface layer (S-layer) protein of Lactobacillus acidophilus is a crystalline array of self-assembling, proteinaceous subunits non-covalently bound to the outmost bacterial cell wall envelope and is involved in the adherence of bacteria to host cells. We have previously described that the S-layer protein of L. acidophilus possesses anti-viral and anti-bacterial properties. In this work, we extracted and purified S-layer proteins from L. acidophilus ATCC 4356 cells to study their interaction with cell wall components from prokaryotic (i.e., peptidoglycan and lipoteichoic acids) and eukaryotic origin (i.e., mucin and chitin), as well as with viruses, bacteria, yeast, and blood cells. Using chimeric S-layer fused to green fluorescent protein (GFP) from different parts of the protein, we analyzed their binding capacity. Our results show that the C-terminal part of the S-layer protein presents lectin-like activity, interacting with different glycoepitopes. We further demonstrate that lipoteichoic acid (LTA) serves as an anchor for the S-layer protein. Finally, a structure for the C-terminal part of S-layer and possible binding sites were predicted by a homology-based model.

Binding of IFT22 to the intraflagellar transport complex is essential for flagellum assembly

Intraflagellar transport (IFT) relies on motor proteins and the IFT complex to construct cilia and flagella. The IFT complex subunit IFT22/RabL5 has sequence similarity with small GTPases although the nucleotide specificity is unclear because of non-conserved G4/G5 motifs. We show that IFT22 specifically associates with G-nucleotides and present crystal structures of IFT22 in complex with GDP, GTP, and with IFT74/81. Our structural analysis unravels an unusual GTP/GDP-binding mode of IFT22 bypassing the classical G4 motif. The GTPase switch regions of IFT22 become ordered upon complex formation with IFT74/81 and mediate most of the IFT22-74/81 interactions. Structure-based mutagenesis reveals that association of IFT22 with the IFT complex is essential for flagellum construction in Trypanosoma brucei although IFT22 GTP-loading is not strictly required.

A single ectomycorrhizal plant root system includes a diverse and spatially structured fungal community

Although only a relatively small proportion of plant species form ectomycorrhizae with fungi, it is crucial for growth and survival for a number of widespread woody plant species. Few studies have attempted to investigate the fine scale spatial structure of entire root systems of adult ectomycorrhizal (EcM) plants. Here, we use the herbaceous perennial Bistorta vivipara to map the entire root system of an adult EcM plant and investigate the spatial structure of its root-associated fungi. All EcM root tips were sampled, mapped and identified using a direct PCR approach and Sanger sequencing of the internal transcribed spacer region. A total of 32.1% of all sampled root tips (739 of 2302) were successfully sequenced and clustered into 41 operational taxonomic units (OTUs). We observed a clear spatial structuring of the root-associated fungi within the root system. Clusters of individual OTUs were observed in the younger parts of the root system, consistent with observations of priority effects in previous studies, but were absent from the older parts of the root system. This may suggest a succession and fragmentation of the root-associated fungi even at a very fine scale, where competition likely comes into play at different successional stages within the root system.