Rapid automatic detection and alignment of repeats in protein sequences

Reconstitution of a biofilm adhesin system from a sulfate-reducing bacterium in Pseudomonas fluorescens

Biofilms of sulfate-reducing bacterium (SRB) like Desulfovibrio vulgaris Hildenborough (DvH) can facilitate metal corrosion in various industrial and environmental settings leading to substantial economic losses. Although the mechanisms of biofilm formation by DvH are not yet well understood, recent studies indicate the large adhesin, DvhA, is a key determinant of biofilm formation. The dvhA gene neighborhood resembles the biofilm-regulating Lap system of Pseudomonas fluorescens but is curiously missing the c-di-GMP-binding regulator LapD. Instead, DvH encodes an evolutionarily unrelated c-di-GMP-binding protein (DVU1020) that we hypothesized is functionally analogous to LapD. To study this unusual Lap system and overcome experimental limitations with the slow-growing anaerobe DvH, we reconstituted its predicted SRB Lap system in a P. fluorescens strain lacking its native Lap regulatory components (ΔlapGΔlapD). Our data support the model that DvhA is a cell surface-associated LapA-like adhesin with a N-terminal "retention module" and that DvhA is released from the cell surface upon cleavage by the LapG-like protease DvhG. Further, we demonstrate DVU1020 (named here DvhD) represents a distinct class of c-di-GMP-binding, biofilm-regulating proteins that regulates DvhG activity in response to intracellular levels of this second messenger. This study provides insight into the key players responsible for biofilm formation by DvH, thereby expanding our understanding of Lap-like systems.

Target-agnostic identification of human antibodies to Plasmodium falciparum sexual forms reveals cross stage recognition of glutamate-rich repeats

Circulating sexual stages of Plasmodium falciparum (Pf) can be transmitted from humans to mosquitoes, thereby furthering the spread of malaria in the population. It is well established that antibodies (Abs) can efficiently block parasite transmission. In search for naturally acquired Ab targets on sexual stages, we established an efficient method for target-agnostic single B cell activation followed by high-throughput selection of human monoclonal antibodies (mAbs) reactive to sexual stages of Pf in the form of gamete and gametocyte extract. We isolated mAbs reactive against a range of Pf proteins including well-established targets Pfs48/45 and Pfs230. One mAb, B1E11K, was cross-reactive to various proteins containing glutamate-rich repetitive elements expressed at different stages of the parasite life cycle. A crystal structure of two B1E11K Fab domains in complex with its main antigen, RESA, expressed on asexual blood stages, showed binding of B1E11K to a repeating epitope motif in a head-to-head conformation engaging in affinity-matured homotypic interactions. Thus, this mode of recognition of Pf proteins, previously described only for PfCSP, extends to other repeats expressed across various stages. The findings augment our understanding of immune-pathogen interactions to repeating elements of the Plasmodium parasite proteome and underscore the potential of the novel mAb identification method used to provide new insights into the natural humoral immune response against Pf .

Impact Statement

A naturally acquired human monoclonal antibody recognizes proteins expressed at different stages of the Plasmodium falciparum lifecycle through affinity-matured homotypic interactions with glutamate-rich repeats.

Intrinsically disordered Prosystemin discloses biologically active repeat motifs

The in-depth studies over the years on the defence barriers by tomato plants have shown that the Systemin peptide controls the response to a wealth of environmental stress agents. This multifaceted stress reaction seems to be related to the intrinsic disorder of its precursor protein, Prosystemin (ProSys). Since latest findings show that ProSys has biological functions besides Systemin sequence, here we wanted to assess if this precursor includes peptide motifs able to trigger stress-related pathways. Candidate peptides were identified in silico and synthesized to test their capacity to trigger defence responses in tomato plants against different biotic stressors. Our results demonstrated that ProSys harbours several repeat motifs which triggered plant immune reactions against pathogens and pest insects. Three of these peptides were detected by mass spectrometry in plants expressing ProSys, demonstrating their effective presence in vivo. These experimental data shed light on unrecognized functions of ProSys, mediated by multiple biologically active sequences which may partly account for the capacity of ProSys to induce defense responses to different stress agents.

In silico and experimental analysis of the repeated domains in BvaP, a protein important for GBS vaginal colonization

Streptococcus agalactiae (group B strep, GBS) infections in neonates are often fatal and strongly associated with maternal GBS vaginal colonization. Previously, we highlighted the importance of a formerly uncharacterized protein, BvaP, in GBS vaginal colonization. BvaP is highly conserved across GBS and is made up of repeated domains, with a variable number of repeats between strains. Here, we evaluate the prevalence of BvaP repeated domains and their relevance in phenotypes previously associated with vaginal colonization. Using in silico analysis, we found that the number of repeats in the BvaP protein does not generally appear to be associated with serotype, isolation site, or host. Using BvaP truncations in GBS strain A909, we determined that a smaller number of repeats was correlated with decreased bacterial chain length, but adherence to vaginal epithelial cells was complemented using BvaP containing one, two, three, or five repeats. Future research will be geared toward understanding the host immune response to BvaP in vivo and whether vaginal carriage or host response is dependent on the BvaP repeated domains.

Invasive mussels fashion silk-like byssus via mechanical processing of massive horizontally acquired coiled coils

Zebra and quagga mussels (Dreissena spp.) are invasive freshwater biofoulers that perpetrate devastating economic and ecological impact. Their success depends on their ability to anchor onto substrates with protein-based fibers known as byssal threads. Yet, compared to other mussel lineages, little is understood about the proteins comprising their fibers or their evolutionary history. Here, we investigated the hierarchical protein structure of Dreissenid byssal threads and the process by which they are fabricated. Unique among bivalves, we found that threads possess a predominantly β-sheet crystalline structure reminiscent of spider silk. Further analysis revealed unexpectedly that the Dreissenid thread protein precursors are mechanoresponsive α-helical proteins that are mechanically processed into β-crystallites during thread formation. Proteomic analysis of the byssus secretory organ and byssus fibers revealed a family of ultrahigh molecular weight (354 to 467 kDa) asparagine-rich (19 to 20%) protein precursors predicted to form α-helical coiled coils. Moreover, several independent lines of evidence indicate that the ancestral predecessor of these proteins was likely acquired via horizontal gene transfer. This chance evolutionary event that transpired at least 12 Mya has endowed Dreissenids with a distinctive and effective fiber formation mechanism, contributing significantly to their success as invasive species and possibly, inspiring new materials design.

Genetic diversity and in silico analysis of Plasmodium knowlesi Serine Repeat Antigen (SERA) 3 antigen 2 in Malaysia

Plasmodium knowlesi is the leading cause of malaria in Malaysia. Serine Repeat Antigens (SERAs) have an essential role in the parasite life cycle. However, genetic characterization on P. knowlesi SERA3 Ag2 (PkSERA3 Ag2) is lacking. In the present study, nucleotide diversity, natural selection, and haplotypes of PkSERA3 Ag2 in clinical samples from Peninsular Malaysia and Malaysian Borneo were investigated. A total of 50 P. knowlesi clinical samples were collected from Peninsular Malaysia and Malaysian Borneo. The PkSERA3 Ag2 gene was amplified using PCR, and subsequently cloned and sequenced. Genetic diversity, haplotype, natural selection as well as genetic structure and differentiation of PkSERA3 Ag2 were analysed. In addition, in silico analyses were performed to identify repeat motifs, B-cell epitopes, and antigenicity indices of the protein. Analysis of 114 PkSERA3 Ag2 sequences revealed high nucleotide diversity of the gene in Malaysia. A codon-based Z-test indicated that the gene underwent purifying selection. Haplotype and population structure analyses identified two distinct PkSERA3 Ag2 clusters (K = 2, ΔK = 721.14) but no clear genetic distinction between PkSERA3 Ag2 from Peninsular Malaysia and Malaysian Borneo. FST index indicated moderate differentiation of the gene. In silico analyses revealed unique repeat motifs among PkSERA3 Ag2 isolates. Moreover, the amino acid sequence of PkSERA3 Ag2 exhibited potential B-cell epitopes and possessed high antigenicity indices. These findings enhance the understanding of PkSERA3 Ag2 gene as well as its antigenic properties. Further validation is necessary to ascertain the utility of PkSERA3 Ag2 as a serological marker for P. knowlesi infection.

Molecular characterization of the N-terminal half of TasA during amyloid-like assembly and its contribution to Bacillus subtilis biofilm formation

Biofilms are bacterial communities that result from a cell differentiation process leading to the secretion of an extracellular matrix (ECM) by part of the population. In Bacillus subtilis, the main protein component of the ECM is TasA, which forms a fiber-based scaffold that confers structure to the ECM. The N-terminal half of TasA is strongly conserved among Bacillus species and contains a protein domain, the rigid core (RcTasA), which is critical for the structural and functional properties of the recombinant protein. In this study, we demonstrate that recombinantly purified RcTasA in vitro retains biochemical properties previously observed for the entire protein. Further analysis of the RcTasA amino acid sequence revealed two aggregation-prone stretches and a region of imperfect amino acid repeats, which are known to contribute to functional amyloid assembly. Biochemical characterization of these stretches found in RcTasA revealed their amyloid-like capacity in vitro, contributing to the amyloid nature of RcTasA. Moreover, the study of the imperfect amino acid repeats revealed the critical role of residues D64, K68 and D69 in the structural function of TasA. Experiments with versions of TasA carrying the substitutions D64A and K68AD69A demonstrated a partial loss of function of the protein either in the assembly of the ECM or in the stability of the core and amyloid-like properties. Taken together, our findings allow us to better understand the polymerization process of TasA during biofilm formation and provide knowledge into the sequence determinants that promote the molecular behavior of protein filaments in bacteria.

Building alternative splicing and evolution-aware sequence-structure maps for protein repeats

Alternative splicing of repeats in proteins provides a mechanism for rewiring and fine-tuning protein interaction networks. In this work, we developed a robust and versatile method, ASPRING, to identify alternatively spliced protein repeats from gene annotations. ASPRING leverages evolutionary meaningful alternative splicing-aware hierarchical graphs to provide maps between protein repeats sequences and 3D structures. We re-think the definition of repeats by explicitly accounting for transcript diversity across several genes/species. Using a stringent sequence-based similarity criterion, we detected over 5,000 evolutionary conserved repeats by screening virtually all human protein-coding genes and their orthologs across a dozen species. Through a joint analysis of their sequences and structures, we extracted specificity-determining sequence signatures and assessed their implication in experimentally resolved and modelled protein interactions. Our findings demonstrate the widespread alternative usage of protein repeats in modulating protein interactions and open avenues for targeting repeat-mediated interactions.

Deep Learning-based structure modelling illuminates structure and function in uncharted regions of β-solenoid fold space

Repeat proteins are common in all domains of life and exhibit a wide range of functions. One class of repeat protein contains solenoid folds where the repeating unit consists of β-strands separated by tight turns. β-solenoids have distinguishing structural features such as handedness, twist, oligomerisation state, coil shape and size which give rise to their diversity. Characterised β-solenoid repeat proteins are known to form regions in bacterial and viral virulence factors, antifreeze proteins and functional amyloids. For many of these proteins, the experimental structure has not been solved, as they are difficult to crystallise or model. Here we use various deep learning-based structure-modelling methods to discover novel predicted β-solenoids, perform structural database searches to mine further structural neighbours and relate their predicted structure to possible functions. We find both eukaryotic and prokaryotic adhesins, confirming a known functional linkage between adhesin function and the β-solenoid fold. We further identify exceptionally long, flat β-solenoid folds as possible structures of mucin tandem repeat regions and unprecedentedly small β-solenoid structures. Additionally, we characterise a novel β-solenoid coil shape, the FapC Greek key β-solenoid as well as plausible complexes between it and other proteins involved in Pseudomonas functional amyloid fibres.

The role of tandem repeats in bacterial functional amyloids

Repetitivity and modularity of proteins are two related notions incorporated into multiple evolutionary concepts. We discuss whether they may also be essential for functional amyloids. Amyloids are proteins that create very regular and usually highly insoluble fibrils, which are often associated with neurodegeneration. However, recent discoveries showed that amyloid structure of a protein could also be beneficial and desired, e.g., to promote cell adhesion. Functional amyloids are proteins which differ in their characteristics from pathological amyloids, so that the fibril formation could be more under control of an organism. We propose that repeats in the sequence could regulate the aggregation propensity of these proteins. The inclusion of multiple symmetric interactions, due to the presence of the repeats, could be supporting and strengthening the desirable structural properties of functional amyloids. Our results show that tandem repeats in bacterial functional amyloids have a distinct characteristic. The pattern of repeats supports the appropriate level of fibril formation and better controllability of fibril stability. The repeats tend to be more imperfect, which attenuates excessive aggregation propensity. Their desired structure and function are also reinforced by their amino acid profile. Although in the study we focused on bacterial functional amyloids, due to their importance in biofilm formation, we propose that similar mechanisms could be employed in other functional amyloids which are designed by evolution to aggregate in a desirable manner, but not necessarily in pathological amyloids.

Fluorescence Microscopy Study of the Intracellular Sulfur Globule Protein SgpD in the Purple Sulfur Bacterium Allochromatium vinosum

When oxidizing reduced sulfur compounds, the phototrophic sulfur bacterium Allochromatium vinosum forms spectacular sulfur globules as obligatory intracellular-but extracytoplasmic-intermediates. The globule envelope consists of three extremely hydrophobic proteins: SgpA and SgpB, which are very similar and can functionally replace each other, and SgpC which is involved in the expansion of the sulfur globules. The presence of a fourth protein, SgpD, was suggested by comparative transcriptomics and proteomics of purified sulfur globules. Here, we investigated the in vivo function of SgpD by coupling its carboxy-terminus to mCherry. This fluorescent protein requires oxygen for chromophore maturation, but we were able to use it in anaerobically growing A. vinosum provided the cells were exposed to oxygen for one hour prior to imaging. While mCherry lacking a signal peptide resulted in low fluorescence evenly distributed throughout the cell, fusion with SgpD carrying its original Sec-dependent signal peptide targeted mCherry to the periplasm and co-localized it exactly with the highly light-refractive sulfur deposits seen in sulfide-fed A. vinosum cells. Insertional inactivation of the sgpD gene showed that the protein is not essential for the formation and degradation of sulfur globules.

Characterization of two bacterial multi-flavinylated proteins harboring multiple covalent flavin cofactors

In recent years, studies have shown that a large number of bacteria secrete multi-flavinylated proteins. The exact roles and properties, of these extracellular flavoproteins that contain multiple covalently anchored FMN cofactors, are still largely unknown. Herein, we describe the biochemical and structural characterization of two multi-FMN-containing covalent flavoproteins, SaFMN3 from Streptomyces azureus and CbFMN4 from Clostridiaceae bacterium. Based on their primary structure, these proteins were predicted to contain three and four covalently tethered FMN cofactors, respectively. The genes encoding SaFMN3 and CbFMN4 were heterologously coexpressed with a flavin transferase (ApbE) in Escherichia coli, and could be purified by affinity chromatography in good yields. Both proteins were found to be soluble and to contain covalently bound FMN molecules. The SaFMN3 protein was studied in more detail and found to display a single redox potential (-184 mV) while harboring three covalently attached flavins. This is in line with the high sequence similarity when the domains of each flavoprotein are compared. The fully reduced form of SaFMN3 is able to use dioxygen as electron acceptor. Single domains from both proteins were expressed, purified and crystallized. The crystal structures were elucidated, which confirmed that the flavin cofactor is covalently attached to a threonine. Comparison of both crystal structures revealed a high similarity, even in the flavin binding pocket. Based on the crystal structure, mutants of the SaFMN3-D2 domain were designed to improve its fluorescence quantum yield by changing the microenvironment of the isoalloxazine moiety of the flavin cofactor. Residues that quench the flavin fluorescence were successfully identified. Our study reveals biochemical details of multi-FMN-containing proteins, contributing to a better understanding of their role in bacteria and providing leads to future utilization of these flavoprotein in biotechnology.

The saclayvirus Aci01-1 very long and complex fiber and its receptor at the Acinetobacter baumannii surface

The Acinetobacter baumannii bacteriophage Aci01-1, which belongs to the genus Saclayvirus of the order Caudoviricetes, has an icosahedral head and a contractile rigid tail. We report that Aci01-1 has, attached to the tail conical tip, a remarkable 146-nm-long flexible fiber with seven beads and a terminal knot. Its putative gene coding for a 241.36-kDa tail fiber protein is homologous to genes in Aci01-1-related and unrelated phages. Analysis of its 3D structure using AlphaFold provides a structural model for the fiber observed by electron microscopy. We also identified a putative receptor of the phage on the bacterial capsule that is hypothesized to interact with the Aci01-1 long fiber.

Isolation and infection cycle of a polinton-like virus virophage in an abundant marine alga

Virophages are small double stranded DNA (dsDNA) viruses that can only replicate in a host by co-infecting with another virus. Marine algae are commonly associated with virophage-like elements such as Polinton-like viruses (PLVs) that remain largely uncharacterized. Here we isolated a PLV that co-infects the alga Phaeocystis globosa with the Phaeocystis globosa virus-14T (PgV-14T), a close relative of the "Phaeocystis globosa virus-virophage" genomic sequence. We name this PLV 'Gezel-14T. Gezel is phylogenetically distinct from the Lavidaviridae family where all known virophages belong. Gezel-14T co-infection decreases the fitness of its viral host by reducing burst sizes of PgV-14T, yet insufficiently to spare the cellular host population. Genomic screens show Gezel-14T-like PLVs integrated into Phaeocystis genomes, suggesting that these widespread viruses are capable of integration into cellular host genomes. This system presents an opportunity to better understand the evolution of eukaryotic dsDNA viruses as well as the complex dynamics and implications of viral parasitism.

The flax genome reveals orbitide diversity

Background

Ribosomally-synthesized cyclic peptides are widely found in plants and exhibit useful bioactivities for humans. The identification of cyclic peptide sequences and their precursor proteins is facilitated by the growing number of sequenced genomes. While previous research largely focused on the chemical diversity of these peptides across various species, there is little attention to a broader range of potential peptides that are not chemically identified.

Results

A pioneering study was initiated to explore the genetic diversity of linusorbs, a group of cyclic peptides uniquely occurring in cultivated flax (Linum usitatissimum). Phylogenetic analysis clustered the 5 known linusorb precursor proteins into two clades and one singleton. Preliminary tBLASTn search of the published flax genome using the whole protein sequence as query could only retrieve its homologues within the same clade. This limitation was overcome using a profile-based mining strategy. After genome reannotation, a hidden Markov Model (HMM)-based approach identified 58 repeats homologous to the linusorb-embedded repeats in 8 novel proteins, implying that they share common ancestry with the linusorb-embedded repeats. Subsequently, we developed a customized profile composed of a random linusorb-like domain (LLD) flanked by 5 conserved sites and used it for string search of the proteome, which extracted 281 LLD-containing repeats (LLDRs) in 25 proteins. Comparative analysis of different repeat categories suggested that the 5 conserved flanking sites among the non-homologous repeats have undergone convergent evolution driven by functional selection.

Conclusions

The profile-based mining approach is suitable for analyzing repetitive sequences. The 25 LLDR proteins identified herein represent the potential diversity of cyclic peptides within the flax genome and lay a foundation for further studies on the functions and evolution of these protein tandem repeats.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08735-x.

Phase separation modulates the assembly and dynamics of a polarity-related scaffold-signaling hub

Asymmetric cell division (ACD) produces morphologically and behaviorally distinct cells and is the primary way to generate cell diversity. In the model bacterium Caulobacter crescentus, the polarization of distinct scaffold-signaling hubs at the swarmer and stalked cell poles constitutes the basis of ACD. However, mechanisms involved in the formation of these hubs remain elusive. Here, we show that a swarmer-cell-pole scaffold, PodJ, forms biomolecular condensates both in vitro and in living cells via phase separation. The coiled-coil 4-6 and the intrinsically disordered regions are the primary domains that contribute to biomolecular condensate generation and signaling protein recruitment in PodJ. Moreover, a negative regulation of PodJ phase separation by the stalked-cell-pole scaffold protein SpmX is revealed. SpmX impedes PodJ cell-pole accumulation and affects its recruitment ability. Together, by modulating the assembly and dynamics of scaffold-signaling hubs, phase separation may serve as a general biophysical mechanism that underlies the regulation of ACD in bacteria and other organisms.

Thermostable Proteins from HaCaT Keratinocytes Identify a Wide Breadth of Intrinsically Disordered Proteins and Candidates for Liquid-Liquid Phase Separation

Intrinsically disordered proteins (IDPs) move through an ensemble of conformations which allows multitudinous roles within a cell. Keratinocytes, the predominant cell type in mammalian epidermis, have had only a few individual proteins assessed for intrinsic disorder and its possible contribution to liquid-liquid phase separation (LLPS), especially in regard to what functions or structures these proteins provide. We took a holistic approach to keratinocyte IDPs starting with enrichment via the isolation of thermostable proteins. The keratinocyte protein involucrin, known for its resistance to heat denaturation, served as a marker. It and other thermostable proteins were identified by liquid chromatography tandem mass spectrometry and subjected to extensive bioinformatic analysis covering gene ontology, intrinsic disorder, and potential for LLPS. Numerous proteins unique to keratinocytes and other proteins with shared expression in multiple cell types were identified to have IDP traits (e.g., compositional bias, nucleic acid binding, and repeat motifs). Among keratinocyte-specific proteins, many that co-assemble with involucrin into the cell-specific structure known as the cornified envelope scored highly for intrinsic disorder and potential for LLPS. This suggests intrinsic disorder and LLPS are previously unrecognized traits for assembly of the cornified envelope, echoing the contribution of intrinsic disorder and LLPS to more widely encountered features such as stress granules and PML bodies.

Multiple Viral Protein Genome-Linked Proteins Compensate for Viral Translation in a Positive-Sense Single-Stranded RNA Virus Infection

Viral protein genome-linked (VPg) protein plays an essential role in protein-primed replication of plus-stranded RNA viruses. VPg is covalently linked to the 5' end of the viral RNA genome via a phosphodiester bond typically at a conserved amino acid. Whereas most viruses have a single VPg, some viruses have multiple VPgs that are proposed to have redundant yet undefined roles in viral replication. Here, we use cricket paralysis virus (CrPV), a dicistrovirus that has four nonidentical copies of VPg, as a model to characterize the role of VPg copies in infection. Dicistroviruses contain two main open reading frames (ORFs) that are driven by distinct internal ribosome entry sites (IRESs). We systematically generated single and combinatorial deletions and mutations of VPg1 to VPg4 within the CrPV infectious clone and monitored viral yield in Drosophila S2 cells. Deletion of one to three VPg copies progressively decreased viral yield and delayed viral replication, suggesting a threshold number of VPgs for productive infection. Mass spectrometry analysis of CrPV VPg-linked RNAs revealed viral RNA linkage to either a serine or threonine in VPg, mutations of which in all VPgs attenuated infection. Mutating serine 4 in a single VPg abolished viral infection, indicating a dominant negative effect. Using viral minigenome reporters that monitor dicistrovirus 5' untranslated (UTR) and IRES translation revealed a relationship between VPg copy number and the ratio of distinct IRES translation activities. We uncovered a novel viral strategy whereby VPg copies in dicistrovirus genomes compensate for the relative IRES translation efficiencies to promote infection. IMPORTANCE Genetic duplication is exceedingly rare in small RNA viral genomes, as there is selective pressure to prevent RNA genomes from expanding. However, some small RNA viruses encode multiple copies of a viral protein, most notably an unusual viral protein that is linked to the viral RNA genome. Here, we investigate a family of viruses that contains multiple viral protein genome-linked proteins and reveal a novel viral strategy whereby viral protein copy number counterbalances differences in viral protein synthesis mechanisms.

Mutated KLF4(K409Q) in meningioma binds STRs and activates FGF3 gene expression

Krüppel-like factor 4 (KLF4) is a transcription factor that has been proven necessary for both induction and maintenance of pluripotency and self-renewal. Whole-genome sequencing defined a unique mutation in KLF4 (KLF4K409Q) in human meningiomas. However, the molecular mechanism of this tumor-specific KLF4 mutation is unknown. Using genome-wide high-throughput and focused quantitative transcriptional approaches in human cell lines, primary meningeal cells, and meningioma tumor tissue, we found that a change in the evolutionarily conserved DNA-binding domain of KLF4 alters its DNA recognition preference, resulting in a shift in downstream transcriptional activity. In the KLF4K409Q-specific targets, the normally silent fibroblast growth factor 3 (FGF3) is activated. We demonstrated a neomorphic function of KLF4K409Q in stimulating FGF3 transcription through binding to its promoter and in using short tandem repeats (STRs) located within the locus as enhancers.

F/YGG-motif is an intrinsically disordered nucleic-acid binding motif

Heterogeneous nuclear ribonucleoproteins (hnRNP) function in RNA processing, have RNA-recognition motifs (RRMs) and intrinsically disordered, low-complexity domains (LCDs). While RRMs are drivers of RNA binding, there is only limited knowledge about the RNA interaction by the LCD of some hnRNPs. Here, we show that the LCD of hnRNPA2 interacts with RNA via an embedded Tyr/Gly-rich region which is a disordered RNA-binding motif. RNA binding is maintained upon mutating tyrosine residues to phenylalanines, but abrogated by mutating to alanines, thus we term the RNA-binding region ‘F/YGG motif’. The F/YGG motif can bind a broad range of structured (e.g. tRNA) and disordered (e.g. polyA) RNAs, but not rRNA. As the F/YGG otif can also interact with DNA, we consider it a general nucleic acid-binding motif. hnRNPA2 LCD can form dense droplets, by liquid–liquid phase separation (LLPS). Their formation is inhibited by RNA binding, which is mitigated by salt and 1,6-hexanediol, suggesting that both electrostatic and hydrophobic interactions feature in the F/YGG motif. The D290V mutant also binds RNA, which interferes with both LLPS and aggregation thereof. We found homologous regions in a broad range of RNA- and DNA-binding proteins in the human proteome, suggesting that the F/YGG motif is a general nucleic acid-interaction motif.

Identification, In Silico Characterization, and Differential Expression Profiles of Carotenoid, Xanthophyll, Apocarotenoid Biosynthetic Pathways Genes, and Analysis of Carotenoid and Xanthophyll Accumulation in Heracleum moellendorffii Hance

Heracleum moellendorffii Hance is a non-woody forest plant widely used in China, Korea, and Japan because of its various therapeutic properties. However, the genetic details of the carotenoid pathway (CP), xanthophyll pathway (XP), and apocarotenoid pathway (AP) genes have not been studied. Thus, the CP, XP, and AP genes of H. moellendorffii were detected and analyzed. A total of fifteen genes were identified, of which eight, four, and three belonged to CP, XP, and AP, respectively. All identified genes possessed full open reading frames. Phylogenetic characterization of the identified gene sequences showed the highest similarity with other higher plants. Multiple alignments and 3D dimensional structures showed several diverse conserved motifs, such as the carotene-binding motif, dinucleotide-binding motif, and aspartate or glutamate residues. The results of real-time PCR showed that the CP, XP, and AP genes were highly expressed in leaves, followed by the stems and roots. In total, eight different individual carotenoids were identified using HPLC analysis. The highest individual and total carotenoid content were achieved in the leaves, followed by the stems and roots. This study will provide more information on the gene structure of the CP, XP, and AP genes, which may help to increase the accumulation of carotenoids in H. moellendorffii through genetic engineering. These results could be helpful for further molecular and functional studies of CP, XP, and AP genes.

Gene-Guided Discovery and Ribosomal Biosynthesis of Moroidin Peptides

Moroidin is a bicyclic plant octapeptide with tryptophan side-chain cross-links, originally isolated as a pain-causing agent from the Australian stinging tree Dendrocnide moroides. Moroidin and its analog celogentin C, derived from Celosia argentea, are inhibitors of tubulin polymerization and, thus, lead structures for cancer therapy. However, low isolation yields from source plants and challenging organic synthesis hinder moroidin-based drug development. Here, we present biosynthesis as an alternative route to moroidin-type bicyclic peptides and report that they are ribosomally synthesized and posttranslationally modified peptides (RiPPs) derived from BURP-domain peptide cyclases in plants. By mining 793 plant transcriptomes for moroidin core peptide motifs within BURP-domain precursor peptides, we identified a moroidin cyclase in Japanese kerria, which catalyzes the installation of the tryptophan-indole-centered macrocyclic bonds of the moroidin bicyclic motif in the presence of cupric ions. Based on the kerria moroidin cyclase, we demonstrate the feasibility of producing diverse moroidins including celogentin C in transgenic tobacco plants and report specific cytotoxicity of celogentin C against a lung adenocarcinoma cancer cell line. Our study sets the stage for future biosynthetic development of moroidin-based therapeutics and highlights that mining plant transcriptomes can reveal bioactive cyclic peptides and their underlying cyclases from new source plants.

HPVMD-C: a disease-based mutation database of human papillomavirus in China

Human papillomavirus (HPV) can cause condyloma acuminatum and cervical cancer. Some mutations of these viruses are closely related to the persistent infection of cervical cancer and are ideal cancer vaccine targets. Several databases have been developed to collect HPV sequences, but no HPV mutation database has been published. This paper reports a Chinese HPV mutation database (HPVMD-C), which contains 149 HPV genotypes, 468 HPV mutations, 3409 protein sequences, 4727 domains and 236 epitopes. We analyzed the mutation distribution among HPV genotypes, domains and epitopes. We designed a visualization tool to display these mutations, domains and epitopes and provided more detailed information about the disease, region and related literature. We also proposed an HPV genotype prediction tool, which can predict HPV carcinogenic or non-carcinogenic risk genotypes. We expect that HPVMD-C will complement the existing database and provide valuable resources for HPV vaccine research and cervical cancer treatment. HPVMD-C is freely available at Database URL: http://bioinfo.zstu.edu.cn/hpv.

OUP accepted manuscript

Homology of highly divergent genes often cannot be determined from sequence similarity alone. For example, we recently identified in the aphid Hormaphis cornu a family of rapidly evolving bicycle genes, which encode novel proteins implicated as plant gall effectors, and sequence similarity search methods yielded few putative bicycle homologs in other species. Coding sequence-independent features of genes, such as intron-exon boundaries, often evolve more slowly than coding sequences, however, and can provide complementary evidence for homology. We found that a linear logistic regression classifier using only structural features of bicycle genes identified many putative bicycle homologs in other species. Independent evidence from sequence features and intron locations supported homology assignments. To test the potential roles of bicycle genes in other aphids, we sequenced the genome of a second gall-forming aphid, Tetraneura nigriabdominalis and found that many bicycle genes are strongly expressed in the salivary glands of the gall forming foundress. In addition, bicycle genes are strongly overexpressed in the salivary glands of a non-gall forming aphid, Acyrthosiphon pisum, and in the non-gall forming generations of H. cornu. These observations suggest that Bicycle proteins may be used by multiple aphid species to manipulate plants in diverse ways. Incorporation of gene structural features into sequence search algorithms may aid identification of deeply divergent homologs, especially of rapidly evolving genes involved in host-parasite interactions.

Conformational heterogeneity coupled with β-fibril formation of a scaffold protein involved in chronic mental illnesses

Chronic mental illnesses (CMIs) pose a significant challenge to global health due to their complex and poorly understood etiologies and hence, absence of causal therapies. Research of the past two decades has revealed dysfunction of the disrupted in schizophrenia 1 (DISC1) protein as a predisposing factor involved in several psychiatric disorders. DISC1 is a multifaceted protein that serves myriads of functions in mammalian cells, for instance, influencing neuronal development and synapse maintenance. It serves as a scaffold hub forming complexes with a variety (~300) of partners that constitute its interactome. Herein, using combinations of structural and biophysical tools, we demonstrate that the C-region of the DISC1 protein is highly polymorphic, with important consequences for its physiological role. Results from solid-state NMR spectroscopy and electron microscopy indicate that the protein not only forms symmetric oligomers but also gives rise to fibrils closely resembling those found in certain established amyloid proteinopathies. Furthermore, its aggregation as studied by isothermal titration calorimetry (ITC) is an exergonic process, involving a negative enthalpy change that drives the formation of oligomeric (presumably tetrameric) species as well as β-fibrils. We have been able to narrow down the β-core region participating in fibrillization to residues 716-761 of full-length human DISC1. This region is absent in the DISC1^Δ22aa splice variant, resulting in reduced association with proteins from the dynein motor complex, viz., NDE-like 1 (NDEL1) and lissencephaly 1 (LIS1), which are crucial during mitosis. By employing surface plasmon resonance, we show that the oligomeric DISC1 C-region has an increased affinity and shows cooperativity in binding to LIS1 and NDEL1, in contrast to the noncooperative binding mode exhibited by the monomeric version. Based on the derived structural models, we propose that the association between the binding partners involves two neighboring subunits of DISC1 C-region oligomers. Altogether, our findings highlight the significance of the DISC1 C-region as a crucial factor governing the balance between its physiological role as a multifunctional scaffold protein and aggregation-related aberrations with potential significance for disease.

Bioinformatics-Guided Expansion and Discovery of Graspetides

Graspetides are a class of ribosomally synthesized and post-translationally modified peptide natural products featuring ATP-grasp ligase-dependent formation of macrolactones/macrolactams. These modifications arise from serine, threonine, or lysine donor residues linked to aspartate or glutamate acceptor residues. Characterized graspetides include serine protease inhibitors such as the microviridins and plesiocin. Here, we report an update to Rapid ORF Description and Evaluation Online (RODEO) for the automated detection of graspetides, which identified 3,923 high-confidence graspetide biosynthetic gene clusters. Sequence and co-occurrence analyses doubled the number of graspetide groups from 12 to 24, defined based on core consensus sequence and putative secondary modification. Bioinformatic analyses of the ATP-grasp ligase superfamily suggest that extant graspetide synthetases diverged once from an ancestral ATP-grasp ligase and later evolved to introduce a variety of ring connectivities. Furthermore, we characterized thatisin and iso-thatisin, two graspetides related by conformational stereoisomerism from Lysobacter antibioticus. Derived from a newly identified graspetide group, thatisin and iso-thatisin feature two interlocking macrolactones with identical ring connectivity, as determined by a combination of tandem mass spectrometry (MS/MS), methanolytic, and mutational analyses. NMR spectroscopy of thatisin revealed a cis conformation for a key proline residue, while molecular dynamics simulations, solvent-accessible surface area calculations, and partial methanolytic analysis coupled with MS/MS support a trans conformation for iso-thatisin at the same position. Overall, this work provides a comprehensive overview of the graspetide landscape, and the improved RODEO algorithm will accelerate future graspetide discoveries by enabling open-access analysis of existing and emerging genomes.

MedProDB: A database of Mediator proteins

•

Mediator complex is a key component of transcriptional regulation in eukaryotes.

•

Identification of Mediator subunits was done by using computational approaches.

•

Different physicochemical properties, and functions of Mediators were discussed.

•

We have developed first database of Mediator proteins e.g. MedProDB.

•

MedProDB contains different types of search and browse options, and various tools.

In the last three decades, the multi-subunit Mediator complex has emerged as the key component of transcriptional regulation of eukaryotic gene expression. Although there were initial hiccups, recent advancements in bioinformatics tools contributed significantly to in-silico prediction and characterization of Mediator subunits from several organisms belonging to different eukaryotic kingdoms. In this study, we have developed the first database of Mediator proteins named MedProDB with 33,971 Mediator protein entries. Out of those, 12531, 11545, and 9895 sequences belong to metazoans, plants, and fungi, respectively. Apart from the core information consisting of sequence, length, position, organism, molecular weight, and taxonomic lineage, additional information of each Mediator sequence like aromaticity, hydropathy, instability index, isoelectric point, functions, interactions, repeat regions, diseases, sequence alignment to Mediator subunit family, Intrinsically Disordered Regions (IDRs), Post-translation modifications (PTMs), and Molecular Recognition Features (MoRFs) may be of high utility to the users. Furthermore, different types of search and browse options with four different tools namely BLAST, Smith-Waterman Align, IUPred, and MoRF-Chibi_Light are provided at MedProDB to perform different types of analysis. Being a critical component of the transcriptional machinery and regulating almost all the aspects of transcription, it generated lots of interest in structural and functional studies of Mediator functioning. So, we think that the MedProDB database will be very useful for researchers studying the process of transcription. This database is freely available at www.nipgr.ac.in/MedProDB.

Modulation of the enzymatic activity of the flagellar lytic transglycosylase SltF by rod components, and the scaffolding protein FlgJ in Rhodobacter sphaeroides

Macromolecular cell-envelope-spanning structures such as the bacterial flagellum must traverse the cell wall. Lytic transglycosylases enzymes are capable of enlarging gaps in the peptidoglycan meshwork to allow the efficient assembly of supramolecular complexes. In the periplasmic space, the assembly of the flagellar rod requires the scaffold protein FlgJ, which includes a muramidase domain in the canonical models Salmonella enterica and Escherichia coli. In contrast, in Rhodobacter sphaeroides, FlgJ and the dedicated flagellar lytic transglycosylase SltF are separate entities that interact in the periplasm. In this study we show that sltF is expressed along with the genes encoding the early components of the flagellar hierarchy that include the hook-basal body proteins, making SltF available during the rod assembly. Protein-protein interaction experiments demonstrated that SltF interacts with the rod proteins FliE, FlgB, FlgC, FlgF and FlgG through its C-terminal region. A deletion analysis that divides the C-terminus in two halves revealed that the interacting regions for most of the rod proteins are not redundant. Our results also show that the presence of the rod proteins FliE, FlgB, FlgC, and FlgF displace the previously reported SltF-FlgJ interaction. In addition, we observed modulation of the transglycosylase activity of SltF mediated by FlgB and FlgJ that could be relevant to coordinate rod assembly with cell wall remodeling. In summary, different mechanisms regulate the flagellar lytic transglycosylase, SltF ensuring a timely transcription, a proper localization and a controlled enzymatic activity. Importance Several mechanisms participate in the assembly of cell-envelope-spanning macromolecular structures. The sequential expression of substrates to be exported, selective export, and a specific order of incorporation are some of the mechanisms that stand out to drive an efficient assembly process. In this work we analyze how the structural rod proteins, the scaffold protein FlgJ and the flagellar lytic enzyme SltF, interact in an orderly fashion to assemble the flagellar rod into the periplasmic space. A complex arrangement of transient interactions directs a dedicated flagellar muramidase towards the flagellar rod. All these interactions bring this protein to the proximity of the peptidoglycan wall while also modulating its enzymatic activity. This study suggests how a dynamic network of interactions participates in controlling SltF, a prominent component for flagellar formation.

A Novel In Silico Method for Molecular Mimicry Detection Finds a Formin with the Potential to Manipulate the Maize Cell Cytoskeleton

Molecular mimicry is one of the evolutionary strategies that parasites use to manipulate the host metabolism and perform an effective infection. This phenomenon has been observed in several animal and plant pathosystems. Despite the relevance of this mechanism in pathogenesis, little is known about it in fungus-plant interactions. For that reason, we performed an in silico method to select plausible mimicry candidates for the Ustilago maydis-maize interaction. Our methodology used a tripartite sequence comparison between the parasite, the host, and nonparasitic organisms' genomes. Furthermore, we used RNA sequencing information to identify gene coexpression, and we determined subcellular localization to detect potential cases of colocalization in the imitator-imitated pairs. With these approximations, we found a putative extracellular formin in U. maydis with the potential to rearrange the host cell cytoskeleton. In parallel, we detected at least two maize genes involved in the cytoskeleton rearrangement differentially expressed under U. maydis infection; thus, this find increases the expectation for the potential mimicry role of the fungal protein. The use of several sources of data led us to develop a strict and replicable in silico methodology to detect molecular mimicry in pathosystems with enough information available. Furthermore, this is the first time that a genomewide search has been performed to detect molecular mimicry in a U. maydis-maize system. Additionally, to allow the reproducibility of this experiment and the use of this pipeline, we created a Web server called Molecular Mimicry Finder.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Identification of amyloidogenic proteins in the microbiomes of a rat Parkinson's disease model and wild-type rats

Cross seeding between amyloidogenic proteins in the gut is receiving increasing attention as a possible mechanism for initiation or acceleration of amyloid formation by aggregation-prone proteins such as αSN, which is central in the development of Parkinson's disease (PD). This is particularly pertinent in view of the growing number of functional (i.e., benign and useful) amyloid proteins discovered in bacteria. Here we identify two amyloidogenic proteins, Pr12 and Pr17, in fecal matter from PD transgenic rats and their wild type counterparts, based on their stability against dissolution by formic acid (FA). Both proteins show robust aggregation into ThT-positive aggregates that contain higher-order β-sheets and have a fibrillar morphology, indicative of amyloid proteins. In addition, Pr17 aggregates formed in vitro showed significant resistance against FA, suggesting an ability to form highly stable amyloid. Treatment with proteinase K revealed a protected core of approx. 9 kDa. Neither Pr12 nor Pr17, however, affected αSN aggregation in vitro. Thus, amyloidogenicity does not per se lead to an ability to cross-seed fibrillation of αSN. Our results support the use of proteomics and FA to identify amyloidogenic protein in complex mixtures and suggests that there may be numerous functional amyloid proteins in microbiomes.

Functional Domains and Evolutionary History of the PMEL and GPNMB Family Proteins

The ancient paralogs premelanosome protein (PMEL) and glycoprotein nonmetastatic melanoma protein B (GPNMB) have independently emerged as intriguing disease loci in recent years. Both proteins possess common functional domains and variants that cause a shared spectrum of overlapping phenotypes and disease associations: melanin-based pigmentation, cancer, neurodegenerative disease and glaucoma. Surprisingly, these proteins have yet to be shown to physically or genetically interact within the same cellular pathway. This juxtaposition inspired us to compare and contrast this family across a breadth of species to better understand the divergent evolutionary trajectories of two related, but distinct, genes. In this study, we investigated the evolutionary history of PMEL and GPNMB in clade-representative species and identified TMEM130 as the most ancient paralog of the family. By curating the functional domains in each paralog, we identified many commonalities dating back to the emergence of the gene family in basal metazoans. PMEL and GPNMB have gained functional domains since their divergence from TMEM130, including the core amyloid fragment (CAF) that is critical for the amyloid potential of PMEL. Additionally, the PMEL gene has acquired the enigmatic repeat domain (RPT), composed of a variable number of imperfect tandem repeats; this domain acts in an accessory role to control amyloid formation. Our analyses revealed the vast variability in sequence, length and repeat number in homologous RPT domains between craniates, even within the same taxonomic class. We hope that these analyses inspire further investigation into a gene family that is remarkable from the evolutionary, pathological and cell biology perspectives.

Immunoinformatics and Pepscan strategies on the path of a peptide-based serological diagnosis of COVID19

Several diagnostic tools have been developed for clinical and epidemiological assays. RT-PCR and antigen detection tests are more useful for diagnosis of acute disease, while antibody tests allow the estimation of exposure in the population. Currently, there is an urgent need for the development of diagnostic tests for COVID-19 that can be used for large-scale epidemiological sampling. Through a comprehensive strategy, potential 16 mer antigenic peptides suited for antibody-based SARS-CoV-2 diagnosis were identified. A systematic scan of the three structural proteins (S,N and M) and the non-structural proteins (ORFs) present in the SARS-CoV-2 virus was conducted through the combination of immunoinformatic methods, peptide SPOT synthesis and an immunoassay with cellulose-bound peptides (Pepscan). The Pepscan filter paper sheets with synthetic peptides were tested against pools of sera of COVID-19 patients. Antibody recognition showed a strong signal for peptides corresponding to the S, N and M proteins of SARS-CoV-2 virus, but not for the ORFs proteins. The peptides exhibiting higher signal intensity were found in the C-terminal region of the N protein. Several peptides of this region showed strong recognition with all three immunoglobulins in the pools of sera. The differential reactivity observed between the different immunoglobulin isotypes (IgA, IgM and IgG) within different regions of the S and N proteins, can be advantageous for ensuring accurate diagnosis of all infected patients, with different times of exposure to infection. Few peptides of the M protein showed antibody recognition and no recognition was observed for peptides of the ORFs proteins.

Building blocks and blueprints for bacterial autolysins

Bacteria utilize a wide variety of endogenous cell wall hydrolases, or autolysins, to remodel their cell walls during processes including cell division, biofilm formation, and programmed death. We here systematically investigate the composition of these enzymes in order to gain insights into their associated biological processes, potential ways to disrupt them via chemotherapeutics, and strategies by which they might be leveraged as recombinant antibacterial biotherapies. To do so, we developed LEDGOs (lytic enzyme domains grouped by organism), a pipeline to create and analyze databases of autolytic enzyme sequences, constituent domain annotations, and architectural patterns of multi-domain enzymes that integrate peptidoglycan binding and degrading functions. We applied LEDGOs to eight pathogenic bacteria, gram negatives Acinetobacter baumannii, Klebsiella pneumoniae, Neisseria gonorrhoeae, and Pseudomonas aeruginosa; and gram positives Clostridioides difficile, Enterococcus faecium, Staphylococcus aureus, and Streptococcus pneumoniae. Our analysis of the autolytic enzyme repertoires of these pathogens reveals commonalities and differences in their key domain building blocks and architectures, including correlations and preferred orders among domains in multi-domain enzymes, repetitions of homologous binding domains with potentially complementarity recognition modalities, and sequence similarity patterns indicative of potential divergence of functional specificity among related domains. We have further identified a variety of unannotated sequence regions within the lytic enzymes that may themselves contain new domains with important functions.

REP2: A Web Server to Detect Common Tandem Repeats in Protein Sequences

Ensembles of tandem repeats (TRs) in protein sequences expand rapidly to form domains well suited for interactions with proteins. For this reason, they are relatively frequent. Some TRs have known structures and therefore it is advantageous to predict their presence in a protein sequence. However, since most TRs diverge quickly, their detection by classical sequence comparison algorithms is not very accurate. Previously, we developed a method and a web server that used curated profiles and thresholds for the detection of 11 common TRs. Here we present a new web server (REP2) that allows the analysis of TRs in both individual and aligned sequences. We provide currently precomputed analyses for a selection of 78 UniProt reference proteomes. We illustrate how these data can be used to study the evolution of TRs using comparative genomics. REP2 can be accessed at http://cbdm-01.zdv.uni-mainz.de/~munoz/rep/.

Revealing biophysical properties of KfrA-type proteins as a novel class of cytoskeletal, coiled-coil plasmid-encoded proteins

Background

DNA binding KfrA-type proteins of broad-host-range bacterial plasmids belonging to IncP-1 and IncU incompatibility groups are characterized by globular N-terminal head domains and long alpha-helical coiled-coil tails. They have been shown to act as transcriptional auto-regulators.

Results

This study was focused on two members of the growing family of KfrA-type proteins encoded by the broad-host-range plasmids, R751 of IncP-1β and RA3 of IncU groups. Comparative in vitro and in silico studies on KfrA_R751 and KfrA_RA3 confirmed their similar biophysical properties despite low conservation of the amino acid sequences. They form a wide range of oligomeric forms in vitro and, in the presence of their cognate DNA binding sites, they polymerize into the higher order filaments visualized as “threads” by negative staining electron microscopy. The studies revealed also temperature-dependent changes in the coiled-coil segment of KfrA proteins that is involved in the stabilization of dimers required for DNA interactions.

Conclusion

KfrA_R751 and KfrA_RA3 are structural homologues. We postulate that KfrA type proteins have moonlighting activity. They not only act as transcriptional auto-regulators but form cytoskeletal structures, which might facilitate plasmid DNA delivery and positioning in the cells before cell division, involving thermal energy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-020-02079-w.

Leaping into the Unknown World of Sporisorium scitamineum Candidate Effectors

Sporisorium scitamineum is a biotrophic fungus causing sugarcane smut disease. In this study, we set up a pipeline and used genomic and dual transcriptomic data previously obtained by our group to identify candidate effectors of S. scitamineum and their expression profiles in infected smut-resistant and susceptible sugarcane plants. The expression profile of different genes after infection in contrasting sugarcane genotypes assessed by RT-qPCR depended on the plant genotypes and disease progression. Three candidate effector genes expressed earlier only in resistant plants, four expressed in both genotypes, and three later in susceptible plants. Ten genes were cloned and transiently expressed in N. benthamiana leaves to determine their subcellular location, while four localized in more than one compartment. Two candidates, g3890 having a nucleoplasmic and mitochondrial location and g5159 targeting the plant cell wall, were selected to obtain their possible corresponding host targets using co-immunoprecipitation (CoIP) experiments and mass spectrometry. Various potential interactors were identified, including subunits of the protein phosphatase 2A and an endochitinase. We investigated the presence of orthologs in sugarcane and using transcriptome data present their expression profiles. Orthologs of sugarcane shared around 70% similarity. Identifying a set of putative fungal effectors and their plant targets provides a valuable resource for functional characterization of the molecular events leading to smut resistance in sugarcane plants and uncovers further opportunities for investigation.

HPREP: a comprehensive database for human proteome repeats

Amino acid repeats are found to play important roles in both structures and functions of the proteins. These are commonly found in all kingdoms of life, especially in eukaryotes and a larger fraction of human proteins composed of repeats. Further, the abnormal expansions of shorter repeats cause various diseases to humans. Therefore, the analysis of repeats of the entire human proteome along with functional, mutational and disease information would help to better understand their roles in proteins. To fulfill this need, we developed a web database HPREP (http://bioinfo.bdu.ac.in/hprep) for human proteome repeats using Perl and HTML programming. We identified different categories of well-characterized repeats and domain repeats that are present in the human proteome of UniProtKB/Swiss-Prot by using in-house Perl programming and novel repeats by using the repeat detection T-REKS tool as well as XSTREAM web server. Further, these proteins are annotated with functional, mutational and disease information and grouped according to specific repeat types. The developed database enables the users to search by specific repeat type in order to understand their involvement in proteins. Thus, the HPREP database is expected to be a useful resource to gain better insight regarding the different repeats in human proteome and their biological roles.

Computational Identification of Repeat-Containing Proteins and Systems

Repetitive sequence elements in proteins and nucleic acids are often signatures of adaptive or reprogrammable systems in nature. Known examples of these systems, such as transcriptional activator-like effectors (TALE) and CRISPR, have been harnessed as powerful molecular tools with a wide range of applications including genome editing. The continued expansion of genomic sequence databases raises the possibility of prospectively identifying new such systems by computational mining. By leveraging sequence repeats as an organizing principle, here we develop a systematic genome mining approach to explore new types of naturally adaptive systems, five of which are discussed in greater detail. These results highlight the existence of a diverse range of intriguing systems in nature that remain to be explored and also provide a framework for future discovery efforts.

Identification of a Cytopathogenic Toxin from Sneathia amnii

Sneathia amnii is a poorly characterized emerging pathogen that has been implicated in amnionitis and urethritis. We found that S. amnii damages fetal membranes, and we identified and purified a cytotoxic exotoxin that lyses human red blood cells and damages cells from fetal membranes. The gene appears to be cotranscribed with a second gene that encodes a protein with identity to two-partner system transporters, suggesting that it is the "A," or secreted component of a type Vb system. The toxin is 1,881 amino acids with a molecular weight of approximately 200 kDa. It binds to red blood cell membranes and forms pores with a diameter of 2.0 to 3.0 nm, resulting in osmolysis. Because it appears to be the "A" or passenger component of a two-partner system, we propose to name this novel cytotoxin/hemolysin CptA for cytopathogenic toxin component A.IMPORTANCE Sneathia amnii is a very poorly characterized emerging pathogen that can affect pregnancy outcome and cause urethritis and other infections. To date, nothing is known about its virulence factors or pathogenesis. We have identified and isolated a cytotoxin, named CptA for cytopathogenic toxin, component A, that is produced by S. amnii CptA is capable of permeabilizing chorionic trophoblasts and lysing human red blood cells and, thus, may play a role in virulence. Except for small domains conserved among two-partner secretion system passenger proteins, the cytotoxin exhibits little amino acid sequence homology to known toxins. In this study, we demonstrate the pore-forming activity of this novel toxin.

Self-analysis of repeat proteins reveals evolutionarily conserved patterns

BACKGROUND

Protein repeats can confound sequence analyses because the repetitiveness of their amino acid sequences lead to difficulties in identifying whether similar repeats are due to convergent or divergent evolution. We noted that the patterns derived from traditional "dot plot" protein sequence self-similarity analysis tended to be conserved in sets of related repeat proteins and this conservation could be quantitated using a Jaccard metric.

RESULTS

Comparison of these dot plots obviated the issues due to sequence similarity for analysis of repeat proteins. A high Jaccard similarity score was suggestive of a conserved relationship between closely related repeat proteins. The dot plot patterns decayed quickly in the absence of selective pressure with an expected loss of 50% of Jaccard similarity due to a loss of 8.2% sequence identity. To perform method testing, we assembled a standard set of 79 repeat proteins representing all the subgroups in RepeatsDB. Comparison of known repeat and non-repeat proteins from the PDB suggested that the information content in dot plots could be used to identify repeat proteins from pure sequence with no requirement for structural information. Analysis of the UniRef90 database suggested that 16.9% of all known proteins could be classified as repeat proteins. These 13.3 million putative repeat protein chains were clustered and a significant amount (82.9%) of clusters containing between 5 and 200 members were of a single functional type.

CONCLUSIONS

Dot plot analysis of repeat proteins attempts to obviate issues that arise due to the sequence degeneracy of repeat proteins. These results show that this kind of analysis can efficiently be applied to analyze repeat proteins on a large scale.

The origin of domestication genes in goats

Goat domestication was critical for agriculture and civilization, but its underlying genetic changes and selection regimes remain unclear. Here, we analyze the genomes of worldwide domestic goats, wild caprid species, and historical remains, providing evidence of an ancient introgression event from a West Caucasian tur-like species to the ancestor of domestic goats. One introgressed locus with a strong signature of selection harbors the MUC6 gene, which encodes a gastrointestinally secreted mucin. Experiments revealed that the nearly fixed introgressed haplotype confers enhanced immune resistance to gastrointestinal pathogens. Another locus with a strong signal of selection may be related to behavior. The selected alleles at these two loci emerged in domestic goats at least 7200 and 8100 years ago, respectively, and increased to high frequencies concurrent with the expansion of the ubiquitous modern mitochondrial haplogroup A. Tracking these archaeologically cryptic evolutionary transformations provides new insights into the mechanisms of animal domestication.

FiRES: A computational method for the de novo identification of internal structure similarity in proteins

Internal structure similarity in proteins can be observed at the domain and subdomain levels. From an evolutionary perspective, structurally similar elements may arise divergently by gene duplication and fusion events but may also be the product of convergent evolution under physicochemical constraints. The characterization of proteins that contain repeated structural elements has implications for many fields of protein science including protein domain evolution, structure classification, structure prediction, and protein engineering. FiRES (Find Repeated Elements in Structure) is an algorithm that relies on a topology-independent structure alignment method to identify repeating elements in protein structure. FiRES was tested against two hand curated databases of protein repeats: MALIDUP, for very divergent duplicated domains; and RepeatsDB for short tandem repeats. The performance of FiRES was compared to that of lalign, RADAR, HHrepID, CE-symm, ReUPred, and Swelfe. FiRES was the method that most accurately detected proteins either with duplicated domains (accuracy = 0.86) or with multiple repeated units (accuracy = 0.92). FiRES is a new methodology for the discovery of proteins containing structurally similar elements. The FiRES web server is publicly available at http://fires.ifc.unam.mx. The scripts, results, and benchmarks from this study can be downloaded from https://github.com/Claualvarez/fires.

Identification and Analysis of Long Repeats of Proteins at the Domain Level

Amino acid repeats play an important role in the structure and function of proteins. Analysis of long repeats in protein sequences enables one to understand their abundance, structure and function in the protein universe. In the present study, amino acid repeats of length >50 (long repeats) were identified in a non-redundant set of UniProt sequences using the RADAR program. The underlying structures and functions of these long repeats were carried out using the Gene3D for structural domains, Pfam for functional domains and enzyme and non-enzyme functional classification for catalytic and binding of the proteins. From a structural perspective, these long repeats seem to predominantly occur in certain architectures such as sandwich, bundle, barrel, and roll and within these architectures abundant in the superfolds. The lengths of the repeats within each fold are not uniform exhibiting different structures for different functions. We also observed that long repeats are in the domain regions of the family and are involved in the function of the proteins. After grouping based on enzyme and non-enzyme classes, we observed the abundant occurrence of long repeats in specific catalytic and binding of the proteins. In this study, we have analyzed the occurrence of long repeats in the protein sequence universe apart from well-characterized short tandem repeats in sequences and their structures and functions of the proteins at the domain level. The present study suggests that long repeats may play an important role in the structure and function of domains of the proteins.

Repeatability in protein sequences

Low complexity regions (LCRs) in protein sequences have special properties that are very different from those of globular proteins. The rules that define secondary structure elements do not apply when the distribution of amino acids becomes biased. While there is a tendency towards structural disorder in LCRs, various examples, and particularly homorepeats of single amino acids, suggest that very short repeats could adopt structures very difficult to predict. These structures are possibly variable and dependant on the context of intra- or inter-molecular interactions. In general, short repeats in LCRs can induce structure. This could explain the observation that very short (non-perfect) repeats are widespread and many define regions with a function in protein interactions. For these reasons, we have developed an algorithm to quickly analyze local repeatability along protein sequences, that is, how close a protein fragment is from a perfect repeat. Using this algorithm we identified that the proteins of the yeast Saccharomyces cerevisiae are depleted in short repeats (approximate or not) of odd-length, while the human proteins are not, that the fish Danio rerio has many proteins with repeats of length two and that the plant Arabidopsis thaliana has an unusually large amount of repeats of length seven. Our method (REpeatability Scanner, RES, accessible at http://cbdm-01.zdv.uni-mainz.de/~munoz/res/) allows to find regions with approximate short repeats in protein sequences, and helps to characterize the variable use of LCRs and compositional bias in different organisms.

Allelic Diversity of Acetyl Coenzyme A Carboxylase accD/bccp Genes Implicated in Nuclear-Cytoplasmic Conflict in the Wild and Domesticated Pea (Pisum sp.)

Reproductive isolation is an important component of species differentiation. The plastid accD gene coding for the acetyl-CoA carboxylase subunit and the nuclear bccp gene coding for the biotin carboxyl carrier protein were identified as candidate genes governing nuclear-cytoplasmic incompatibility in peas. We examined the allelic diversity in a set of 195 geographically diverse samples of both cultivated (Pisum sativum, P. abyssinicum) and wild (P. fulvum and P. elatius) peas. Based on deduced protein sequences, we identified 34 accD and 31 bccp alleles that are partially geographically and genetically structured. The accD is highly variable due to insertions of tandem repeats. P. fulvum and P. abyssinicum have unique alleles and combinations of both genes. On the other hand, partial overlap was observed between P. sativum and P. elatius. Mapping of protein sequence polymorphisms to 3D structures revealed that most of the repeat and indel polymorphisms map to sequence regions that could not be modeled, consistent with this part of the protein being less constrained by requirements for precise folding than the enzymatically active domains. The results of this study are important not only from an evolutionary point of view but are also relevant for pea breeding when using more distant wild relatives.

A systematic simulation-based meta-analytical framework for prediction of physiological biomarkers in alopecia

Background

Alopecia or hair loss is a complex polygenetic and psychologically devastating disease affecting millions of men and women globally. Since the gene annotation and environmental knowledge is limited for alopecia, a systematic analysis for the identification of candidate biomarkers is required that could provide potential therapeutic targets for hair loss therapy.

Results

We designed an interactive framework to perform a meta-analytical study based on differential expression analysis, systems biology, and functional proteomic investigations. We analyzed eight publicly available microarray datasets and found 12 potential candidate biomarkers including three extracellular proteins from the list of differentially expressed genes with a p-value < 0.05. After expression profiling and functional analysis, we studied protein–protein interactions and observed functional associations of source proteins including WIF1, SPON1, LYZ, GPRC5B, PTPRE, ZFP36L2, HBB, PHF15, LMCD1, KRT35 and VAV3 with target proteins including APCDD1, WNT1, WNT3A, SHH, ESRI, TGFB1, and APP. Pathway analysis of these molecules revealed their role in major physiological reactions including protein metabolism, signal transduction, WNT, BMP, EDA, NOTCH and SHH pathways. These pathways regulate hair growth, hair follicle differentiation, pigmentation, and morphogenesis. We studied the regulatory role of β-catenin, Nf-kappa B, cytokines and retinoic acid in the development of hair growth. Therefore, the differential expression of these significant proteins would affect the normal level and could cause aberrations in hair growth.

Conclusion

Our integrative approach helps to prioritize the biomarkers that ultimately lessen the economic burden of experimental studies. It will also be valuable to discover mutants in genomic data in order to increase the identification of new biomarkers for similar problems.

Electronic supplementary material

The online version of this article (10.1186/s40709-019-0094-x) contains supplementary material, which is available to authorized users.