Protein repeats: structures, functions, and evolution

Unique ankyrin repeat proteins in the genome of poxviruses-Boon or Wane, a critical review

Ankyrin repeat is a 33-amino acid motif commonly observed in eukaryotes and, to a lesser extent, in prokaryotes and archaea and rarely in viruses. This motif plays a crucial role in regulating various cellular processes like the cell cycle, transcription, cell signaling, and inflammatory responses through interactions between proteins. Poxviruses exhibit a distinctive feature of containing multiple ankyrin repeat proteins within their genomes. All the genera of poxviruses possess these proteins except molluscipox virus, crocodylidpox virus, and red squirrel poxvirus. An intriguing characteristic has generated notable interest in studying the functions of these proteins within poxvirus biology. Within poxviruses, ankyrin repeat proteins exhibit a distinct configuration, featuring ankyrin repeats in the N-terminal region and a cellular F-box homolog in the C-terminal region, which enables interactions with the cellular Skp, Cullin, F-box containing ubiquitin ligase complex. Through the examination of experimental evidences and discussions from current literature, this review elucidates the organization and role of ankyrin repeat proteins in poxviruses. Various research studies have highlighted the significant importance of these proteins in poxviral pathogenesis and, acting as factors that enhance virulence. Consequently, they represent viable targets for developing genetically altered viruses with decreased virulence, thus displaying potential as candidates for vaccines and antiviral therapeutic development contributing to safer and more effective strategies against poxviral infections.

Diversity and structural-functional insights of alpha-solenoid proteins

Alpha-solenoids are a significant and diverse subset of structured tandem repeat proteins (STRPs) that are important in various domains of life. This review examines their structural and functional diversity and highlights their role in critical cellular processes such as signaling, apoptosis, and transcriptional regulation. Alpha-solenoids can be classified into three geometric folds: low curvature, high curvature, and corkscrew, as well as eight subfolds: ankyrin repeats; Huntingtin, elongation factor 3, protein phosphatase 2A, and target of rapamycin; armadillo repeats; tetratricopeptide repeats; pentatricopeptide repeats; Pumilio repeats; transcription activator-like; and Sel-1 and Sel-1-like repeats. These subfolds represent distinct protein families with unique structural properties and functions, highlighting the versatility of alpha-solenoids. The review also discusses their association with disease, highlighting their potential as therapeutic targets and their role in protein design. Advances in state-of-the-art structure prediction methods provide new opportunities and challenges in the functional characterization and classification of this kind of fold, emphasizing the need for continued development of methods for their identification and proper data curation and deposition in the main databases.

De novo gene synthesis by an antiviral reverse transcriptase

Defense-associated reverse transcriptase (DRT) systems perform DNA synthesis to protect bacteria against viral infection, but the identities and functions of their DNA products remain largely unknown. We show that DRT2 systems encode an unprecedented immune pathway that involves de novo gene synthesis through rolling circle reverse transcription of a noncoding RNA (ncRNA). Programmed template jumping on the ncRNA generates a concatemeric cDNA, which becomes double-stranded upon viral infection. This DNA product constitutes a protein-coding, nearly endless open reading frame (neo) gene whose expression leads to potent cell growth arrest, restricting the viral infection. Our work highlights an elegant expansion of genome coding potential through RNA-templated gene creation and challenges conventional paradigms of genetic information encoded along the one-dimensional axis of genomic DNA.

Systematic discovery of DNA-binding tandem repeat proteins

Tandem repeat proteins (TRPs) are widely distributed and bind to a wide variety of ligands. DNA-binding TRPs such as zinc finger (ZNF) and transcription activator-like effector (TALE) play important roles in biology and biotechnology. In this study, we first conducted an extensive analysis of TRPs in public databases, and found that the enormous diversity of TRPs is largely unexplored. We then focused our efforts on identifying novel TRPs possessing DNA-binding capabilities. We established a protein language model for DNA-binding protein prediction (PLM-DBPPred), and predicted a large number of DNA-binding TRPs. A subset was then selected for experimental screening, leading to the identification of 11 novel DNA-binding TRPs, with six showing sequence specificity. Notably, members of the STAR (Short TALE-like Repeat proteins) family can be programmed to target specific 9 bp DNA sequences with high affinity. Leveraging this property, we generated artificial transcription factors using reprogrammed STAR proteins and achieved targeted activation of endogenous gene sets. Furthermore, the members of novel families such as MOON (Marine Organism-Originated DNA binding protein) and pTERF (prokaryotic mTERF-like protein) exhibit unique features and distinct DNA-binding characteristics, revealing interesting biological clues. Our study expands the diversity of DNA-binding TRPs, and demonstrates that a systematic approach greatly enhances the discovery of new biological insights and tools.

Exon shuffling and alternative splicing of ROCO genes in brown algae enables a diverse repertoire of candidate immune receptors

The ROCO family is a family of GTPases characterized by a central ROC-COR tandem domain. Interest in the structure and function of ROCO proteins has increased with the identification of their important roles in human disease. Nevertheless, the functions of most ROCO proteins are still unknown. In the present study, we characterized the structure, evolution, and expression of ROCOs in four species of brown algae. Brown algae have a larger number of ROCO proteins than other organisms reported to date. Phylogenetic analyses showed that ROCOs have an ancient origin, likely originated in prokaryotes. ROCOs in brown algae clustered into four groups and showed no strong relationship with red algae or green algae. Brown algal ROCOs retain the ancestral LRR-ROC-COR domain arrangement, which is found in prokaryotes, plants and some basal metazoans. Remarkably, individual LRR motifs in ROCO genes are each encoded by separate exons and exhibit intense exon shuffling and diversifying selection. Furthermore, the tandem LRR exons exhibit alternative splicing to generate multiple transcripts. Both exon shuffling and alternative splicing of LRR repeats may be important mechanisms for generating diverse ligand-binding specificities as immune receptors. Besides their potential immune role, expression analysis shows that many ROCO genes are responsive to other stress conditions, suggesting they could participate in multiple signal pathways, not limited to the immune response. Our results substantially enhance our understanding of the structure and function of this mysterious gene family.

Analysis of LRBA pathogenic variants and the association with functional protein domains and clinical presentation

LRBA is a cytoplasmic protein that is ubiquitously distributed. Almost all LRBA domains have a scaffolding function. In 2012, it was reported that homozygous variants in LRBA are associated with early-onset hypogammaglobulinemia. Since its discovery, more than 100 pathogenic variants have been reported. This review focuses on the variants reported in LRBA and their possible associations with clinical phenotypes. In this work LRBA deficiency cases reported more than 11 years ago have been revised. A database was constructed to analyze the type of variants, age at onset, clinical diagnosis, infections, autoimmune diseases, and cellular and immunoglobulin levels. The review of cases from 2012 to 2023 showed that LRBA deficiency was commonly diagnosed in patients with a clinical diagnosis of Common Variable Immunodeficiency, followed by enteropathy, neonatal diabetes mellitus, ALPS, and X-linked-like syndrome. Most cases show early onset of presentation at <6 years of age. Most cases lack protein expression, whereas hypogammaglobulinemia is observed in half of the cases, and IgG and IgA levels are isotypes reported at low levels. Patients with elevated IgG levels exhibited more than one autoimmune manifestation. Patients carrying pathogenic variants leading to a premature stop codon show a severe phenotype as they have an earlier onset of disease presentation, severe autoimmune manifestations, premature death, and low B cells and regulatory T cell levels. Missense variants were more common in patients with low IgG levels and cytopenia. This work lead to the conclusion that the type of variant in LRBA has association with disease severity, which leads to a premature stop codon being the ones that correlates with severe disease.

Molecular characterization and functional analysis of Eimeria tenella ankyrin repeat-containing protein

Chicken coccidiosis causes disastrous losses to the poultry industry all over the world. Eimeria tenella is the most prevalent of these disease-causing species. Our former RNA-seq indicated that E. tenella ankyrin repeat-containing protein (EtANK) was expressed differently between drug-sensitive (DS) and drug-resistant strains. In this study, we cloned EtANK and analyzed its translational and transcriptional levels using quantitative real-time PCR (qPCR) and western blotting. The data showed that EtANK was significantly upregulated in diclazuril-resistant (DZR) strain and maduramicin-resistant (MRR) strain compared with the drug-sensitive (DS) strain. In addition, the transcription levels in the DZR strains isolated from the field were higher than in the DS strain. The translation levels of EtANK were higher in unsporulated oocysts (UO) than in sporozoites (SZ), sporulated oocysts (SO), or second-generation merozoites (SM), and the protein levels in SM were significantly higher than in UO, SO, and SZ. The results of the indirect immunofluorescence localization showed that the protein was distributed mainly at the anterior region of SZ and on the surface and in the cytoplasm of SM. The fluorescence intensity increased further with its development in vitro. An anti-rEtANK polyclonal antibody inhibited the invasive ability of E. tenella in DF-1 cells. These results showed that EtANK may be related to host cell invasion, required for the parasite's growth in the host, and may be involved in the development of E. tenella resistance to some drugs.

Comparative genomics of two Asian medicinal leeches Hirudo nipponia and Hirudo tianjinensis: With emphasis on antithrombotic genes and their corresponding proteins

Leeches secrete various biologically active substances which have important medical and pharmaceutical values in antithrombotic treatments. Here, we provide a high quality genome of two Asian medicinal leeches Hirudo nipponia and Hirudo tianjinensis, based on which, we identified 22 antithrombotic gene families, including fourteen coagulation inhibitors, four platelet aggregation inhibitors, three fibrinolysis enhancers, and one tissue penetration enhancer. The total numbers of antithrombotic genes were similar between H. nipponia (N = 86) and H. tianjinensis (N = 83). Molecular evolution analysis showed that no significant differences were detected between the two species in any of the three selection indices (dN, dS, and dN/dS), nor in the number of sites under positive/purifying selection. RNA-Seq based gene expression analysis showed that the overall expression patterns of the antithrombotic gene families were not significantly deviated between the two species. Our results indicated that there were rather close similarities between the two leeches on genomic characteristics, especially for the molecular evolution and expression of antithrombotic genes. Our study provides the most comprehensive collection of antithrombotic biomacromolecules from the two Asian medicinal leeches to date. These results will greatly facilitate the research and application of leech derivatives for medical and pharmaceutical purposes of thrombosis.

Terminal regions of a protein are a hotspot for low complexity regions and selection

Volatile low complexity regions (LCRs) are a novel source of adaptive variation, functional diversification and evolutionary novelty. An interplay of selection and mutation governs the composition and length of low complexity regions. High %GC and mutations provide length variability because of mechanisms like replication slippage. Owing to the complex dynamics between selection and mutation, we need a better understanding of their coexistence. Our findings underscore that positively selected sites (PSS) and low complexity regions prefer the terminal regions of genes, co-occurring in most Tetrapoda clades. We observed that positively selected sites within a gene have position-specific roles. Central-positively selected site genes primarily participate in defence responses, whereas terminal-positively selected site genes exhibit non-specific functions. Low complexity region-containing genes in the Tetrapoda clade exhibit a significantly higher %GC and lower ω (dN/dS: non-synonymous substitution rate/synonymous substitution rate) compared with genes without low complexity regions. This lower ω implies that despite providing rapid functional diversity, low complexity region-containing genes are subjected to intense purifying selection. Furthermore, we observe that low complexity regions consistently display ubiquitous prevalence at lower purity levels, but exhibit a preference for specific positions within a gene as the purity of the low complexity region stretch increases, implying a composition-dependent evolutionary role. Our findings collectively contribute to the understanding of how genetic diversity and adaptation are shaped by the interplay of selection and low complexity regions in the Tetrapoda clade.

MdTPR16, an apple tetratricopeptide repeat (TPR)-like superfamily gene, positively regulates drought stress in apple

The Tetratricopeptide repeat (TPR)-like superfamily with TPR conserved domains is widely involved in the growth and abiotic stress in many plants. In this report, the gene MdTPR16 belongs to the TPR family in apple (Malus domestica). Promoter analysis reveal that MdTPR16 incorporated various stress response elements, including the drought stress response elements. And different abiotic stress treatments, drought especially, significantly induce the response of MdTPR16. Overexpression of MdTPR16 result in better drought tolerance in apple and Arabidopsis by up-regulating the expression levels of drought stress-related genes, achieving a higher chlorophyll content level, more material accumulation, and overall better growth compared to WT in the drought. Under drought stress, the overexpressed MdTPR16 also mitigate the oxidative damage in cells by reducing the electrolyte leakage, malondialdehyde content, and the H2O2 and O2- accumulation in apples and Arabidopsis. In conclusion, MdTPR16 act as a beneficial regulator of drought stress response by regulating the expression of related genes and the cumulation of reactive oxygen species (ROS).

Variation in base composition, structure-function relationships, and origins of structural repetition in bacterial rpsA gene

Protein domain repeats are known to arise due to tandem duplications of internal genes. However, the understanding of the underlying mechanisms of this process is incomplete. The goal of this work was to investigate the mechanism of occurrence of repeat expansion based on studying the sequences of 1324 rpsA genes of bacterial S1 ribosomal proteins containing different numbers of S1 structural domains. The rpsA gene encodes ribosomal S1 protein, which is essential for cell viability as it interacts with both mRNA and proteins. Gene ontology (GO) analysis of S1 domains in ribosomal S1 proteins revealed that bacterial protein sequences in S1 mainly have 3 types of molecular functions: RNA binding activity, nucleic acid activity, and ribosome structural component. Our results show that the maximum value of rpsA gene identity for full-length proteins was found for S1 proteins containing six structural domains (58%). Analysis of consensus sequences showed that parts of the rpsA gene encoding separate S1 domains have no a strictly repetitive structure between groups containing different numbers of S1 domains. At the same time, gene regions encoding some conserved residues that form the RNA-binding site remain conserved. The detected phylogenetic similarity suggests that the proposed fold of the rpsA translation initiation region of Escherichia coli has functional value and is important for translational control of rpsA gene expression in other bacterial phyla, but not only in gamma Proteobacteria.

A novel tetratricopeptide-repeat protein, TTP1, forms complexes with glutamyl-tRNA reductase and protochlorophyllide oxidoreductase during tetrapyrrole biosynthesis

The biosynthesis of the tetrapyrrole end-products chlorophyll and heme depends on a multifaceted control mechanism that acts primarily at the post-translational level upon the rate-limiting step of 5-aminolevulinic acid synthesis and upon light-dependent protochlorophyllide oxidoreductase (POR). These regulatory processes require auxiliary factors that modulate the activity, stability, complex formation, and subplastidal localization of the relevant proteins. Together, they ensure optimal metabolic flow during the day and at night. As an Arabidopsis homolog of the POR-interacting tetratricopeptide-repeat protein (Pitt) first reported in Synechocystis, we characterize tetrapyrrole biosynthesis-regulating tetratricopeptide-repeat protein1 (TTP1). TTP1 is a plastid-localized, membrane-bound factor that interacts with POR, the Mg protoporphyrin monomethylester cyclase CHL27, glutamyl-tRNA reductase (GluTR), GluTR-binding protein, and FLUORESCENCE IN BLUE LIGHT. Lack of TTP1 leads to accumulation of GluTR, enhanced 5-aminolevulinic acid synthesis and lower levels of POR. Knockout mutants show enhanced sensitivity to reactive oxygen species and a slower greening of etiolated seedlings. Based on our studies, the interaction of TTP1 with GluTR and POR does not directly inhibit their enzymatic activity and contribute to the control of 5-aminolevulinic acid synthesis. Instead, we propose that TTP1 sequesters a fraction of these proteins on the thylakoid membrane, and contributes to their stability.

Structured Tandem Repeats in Protein Interactions

Tandem repeats (TRs) in protein sequences are consecutive, highly similar sequence motifs. Some types of TRs fold into structural units that pack together in ensembles, forming either an (open) elongated domain or a (closed) propeller, where the last unit of the ensemble packs against the first one. Here, we examine TR proteins (TRPs) to see how their sequence, structure, and evolutionary properties favor them for a function as mediators of protein interactions. Our observations suggest that TRPs bind other proteins using large, structured surfaces like globular domains; in particular, open-structured TR ensembles are favored by flexible termini and the possibility to tightly coil against their targets. While, intuitively, open ensembles of TRs seem prone to evolve due to their potential to accommodate insertions and deletions of units, these evolutionary events are unexpectedly rare, suggesting that they are advantageous for the emergence of the ancestral sequence but are early fixed. We hypothesize that their flexibility makes it easier for further proteins to adapt to interact with them, which would explain their large number of protein interactions. We provide insight into the properties of open TR ensembles, which make them scaffolds for alternative protein complexes to organize genes, RNA and proteins.

Genome-wide identification and analysis of WD40 proteins reveal that NtTTG1 enhances drought tolerance in tobacco (Nicotiana tabacum)

BACKGROUND

WD40 proteins, which are highly prevalent in eukaryotes, play important roles in plant development and stress responses. However, systematic identification and exploration of WD40 proteins in tobacco have not yet been conducted.

RESULTS

In this study, a total of 399 WD40 regulatory genes were identified in common tobacco (Nicotiana tabacum). Gene structure and motif analysis revealed structural and functional diversity among different clades of tobacco WD40 regulatory genes. The expansion of tobacco WD40 regulatory genes was mainly driven by segmental duplication and purifying selection. A potential regulatory network of NtWD40s suggested that NtWD40s might be regulated by miRNAs and transcription factors in various biological processes. Expression pattern analysis via transcriptome analysis and qRT-PCR revealed that many NtWD40s exhibited tissue-specific expression patterns and might be involved in various biotic and abiotic stresses. Furthermore, we have validated the critical role of NtTTG1, which was located in the nuclei of trichome cells, in enhancing the drought tolerance of tobacco plants.

CONCLUSIONS

Our study provides comprehensive information to better understand the evolution of WD40 regulatory genes and their roles in different stress responses in tobacco.

Apoptosis Inhibitor 5: A Multifaceted Regulator of Cell Fate

Apoptosis, or programmed cell death, is a fundamental process that maintains tissue homeostasis, eliminates damaged or infected cells, and plays a crucial role in various biological phenomena. The deregulation of apoptosis is involved in many human diseases, including cancer. One of the emerging players in the intricate regulatory network of apoptosis is apoptosis inhibitor 5 (API5), also called AAC-11 (anti-apoptosis clone 11) or FIF (fibroblast growth factor-2 interacting factor). While it may not have yet the same level of notoriety as some other cancer-associated proteins, API5 has garnered increasing attention in the cancer field in recent years, as elevated API5 levels are often associated with aggressive tumor behavior, resistance to therapy, and poor patient prognosis. This review aims to shed light on the multifaceted functions and regulatory mechanisms of API5 in cell fate decisions as well as its interest as therapeutic target in cancer.

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.

A STRP-ed definition of Structured Tandem Repeats in Proteins

Tandem Repeat Proteins (TRPs) are a class of proteins with repetitive amino acid sequences that have been studied extensively for over two decades. Different features at the level of sequence, structure, function and evolution have been attributed to them by various authors. And yet many of its salient features appear only when looking at specific subclasses of protein tandem repeats. Here, we attempt to rationalize the existing knowledge on Tandem Repeat Proteins (TRPs) by pointing out several dichotomies. The emerging picture is more nuanced than generally assumed and allows us to draw some boundaries of what is not a "proper" TRP. We conclude with an operational definition of a specific subset, which we have denominated STRPs (Structural Tandem Repeat Proteins), which separates a subclass of tandem repeats with distinctive features from several other less well-defined types of repeats. We believe that this definition will help researchers in the field to better characterize the biological meaning of this large yet largely understudied group of proteins.

Revisiting the Asian Buffalo Leech (Hirudinaria manillensis) Genome: Focus on Antithrombotic Genes and Their Corresponding Proteins

Leeches are well-known annelids due to their obligate blood-feeding habits. Some leech species secrete various biologically active substances which have important medical and pharmaceutical value in antithrombotic treatments. In this study, we provided a high-quality genome of the Asian buffalo leech (Hirudinaria manillensis), based on which we performed a systematic identification of potential antithrombotic genes and their corresponding proteins. Combining automatic and manual prediction, we identified 21 antithrombotic gene families including fourteen coagulation inhibitors, three platelet aggregation inhibitors, three fibrinolysis enhancers, and one tissue penetration enhancer. A total of 72 antithrombotic genes, including two pseudogenes, were identified, including most of their corresponding proteins forming three or more disulfide bonds. Three protein families (LDTI, antistasin, and granulin) had internal tandem repeats containing 6, 10, and 12 conserved cysteines, respectively. We also measured the anticoagulant activities of the five identified hirudins (hirudin_Hman1 ~ hirudin_Hman5). The results showed that three (hirudin_Hman1, hirudin_Hman2, and hirudin_Hman5), but not the remaining two, exhibited anticoagulant activities. Our study provides the most comprehensive collection of antithrombotic biomacromolecules from a leech to date. These results will greatly facilitate the research and application of leech derivatives for medical and pharmaceutical purposes in the treatment of thrombotic diseases.

Genome-Wide Identification, Evolution, and Expression Analysis of the WD40 Subfamily in Oryza Genus

The WD40 superfamily is widely found in eukaryotes and has essential subunits that serve as scaffolds for protein complexes. WD40 proteins play important regulatory roles in plant development and physiological processes, such as transcription regulation and signal transduction; it is also involved in anthocyanin biosynthesis. In rice, only OsTTG1 was found to be associated with anthocyanin biosynthesis, and evolutionary analysis of the WD40 gene family in multiple species is less studied. Here, a genome-wide analysis of the subfamily belonging to WD40-TTG1 was performed in nine AA genome species: Oryza sativa ssp. japonica, Oryza sativa ssp. indica, Oryza rufipogon, Oryza glaberrima, Oryza meridionalis, Oryza barthii, Oryza glumaepatula, Oryza nivara, and Oryza longistaminata. In this study, 383 WD40 genes in the Oryza genus were identified, and they were classified into four groups by phylogenetic analysis, with most members in group C and group D. They were found to be unevenly distributed across 12 chromosomes. A total of 39 collinear gene pairs were identified in the Oryza genus, and all were segmental duplications. WD40s had similar expansion patterns in the Oryza genus. Ka/Ks analyses indicated that they had undergone mainly purifying selection during evolution. Furthermore, WD40s in the Oryza genus have similar evolutionary patterns, so Oryza sativa ssp. indica was used as a model species for further analysis. The cis-acting elements analysis showed that many genes were related to jasmonic acid and light response. Among them, OsiWD40-26/37/42 contained elements of flavonoid synthesis, and OsiWD40-15 had MYB binding sites, indicating that they might be related to anthocyanin synthesis. The expression profile analysis at different stages revealed that most OsiWD40s were expressed in leaves, roots, and panicles. The expression of OsiWD40s was further analyzed by qRT-PCR in 9311 (indica) under various hormone treatments and abiotic stresses. OsiWD40-24 was found to be responsive to both phytohormones and abiotic stresses, suggesting that it might play an important role in plant stress resistance. And many OsiWD40s might be more involved in cold stress tolerance. These findings contribute to a better understanding of the evolution of the WD40 subfamily. The analyzed candidate genes can be used for the exploration of practical applications in rice, such as cultivar culture for colored rice, stress tolerance varieties, and morphological marker development.

Exploring the roles of ZmARM gene family in maize development and abiotic stress response

Armadillo (ARM) was a gene family important to plants, with crucial roles in regulating plant growth, development, and stress responses. However, the properties and functions of ARM family members in maize had received limited attention. Therefore, this study employed bioinformatics methods to analyze the structure and evolution of ARM-repeat protein family members in maize. The maize (Zea mays L.) genome contains 56 ARM genes distributed over 10 chromosomes, and collinearity analysis indicated 12 pairs of linkage between them. Analysis of the physicochemical properties of ARM proteins showed that most of these proteins were acidic and hydrophilic. According to the number and evolutionary analysis of the ARM genes, the ARM genes in maize can be divided into eight subgroups, and the gene structure and conserved motifs showed similar compositions in each group. The findings shed light on the significant roles of 56 ZmARM domain genes in development and abiotic stress, particularly drought stress. RNA-Seq and qRT-PCR analysis revealed that drought stress exerts an influence on specific members of the ZmARM family, such as ZmARM4, ZmARM12, ZmARM34 and ZmARM36. The comprehensive profiling of these genes in the whole genome, combined with expression analysis, establishes a foundation for further exploration of plant gene function in the context of abiotic stress and reproductive development.

Diversity and features of proteins with structural repeats

The review provides information on proteins with structural repeats, including their classification, characteristics, functions, and relevance in disease development. It explores methods for identifying structural repeats and specialized databases. The review also highlights the potential use of repeat proteins as drug design scaffolds and discusses their evolutionary mechanisms.

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.

Engineering of brick and staple components for ordered assembly of synthetic repeat proteins

Synthetic ɑRep repeat proteins are engineered as Brick and Staple protein pairs that together self-assemble into helical filaments. In most cases, the filaments spontaneously form supercrystals. Here, we describe an expanded series of ɑRep Bricks designed to stabilize the interaction between consecutive Bricks, to control the length of the assembled multimers, or to alter the spatial distribution of the Staple on the filaments. The effects of these Brick modifications on the assembly, on the final filament structure and on the crystal symmetry are analyzed by biochemical methods, electron microscopy and small angle X-ray scattering. We further extend the concept of Brick/Staple protein origami by designing a new type of "Janus"-like Brick protein that is equally assembled by orthogonal staples binding its inner or outer surfaces and thus ending inside or outside the filaments. The relative roles of longitudinal and lateral associations in the assembly process are discussed. This set of results demonstrates important proofs-of-principle for engineering these remarkably versatile proteins toward nanometer-to-micron scale constructions.

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.

Why does the first protein repeat often become the only one?

Proteins with two similar motifs in tandem are one of the most common cases of tandem repeat proteins. The question arises: why is the first emerged repeat frequently fixed in the process of evolution, despite the ample opportunities to continue its multiplication at the DNA level? To answer this question, we systematically analyzed the structure and function of these proteins. Our analysis showed that, in the vast majority of cases, the structural repetitive units have a two-fold (C2) internal symmetry. These closed structures provide an internal structural limitation for the subsequent growth of the repeat number. Frequently, the units "swap" their secondary structure elements with each other. Moreover, the duplicated domains, in contrast to other tandem repeat proteins, form binding sites for small molecules around the axis of C2 symmetry. Thus, the closure of the C2 structures and the emergence of new functional sites around the axis of C2 symmetry provide plausible explanations for why a repeat, once appeared, becomes fixed in the evolutionary process. We have placed these structures within the general structural classification of tandem repeat proteins, classifying them as either Class IV or V depending on the size of the repetitive unit.

Evolutionary Study of Protein Short Tandem Repeats in Protein Families

Tandem repeats in proteins are patterns of residues repeated directly adjacent to each other. The evolution of these repeats can be assessed by using groups of homologous sequences, which can help pointing to events of unit duplication or deletion. High pressure in a protein family for variation of a given type of repeat might point to their function. Here, we propose the analysis of protein families to calculate protein short tandem repeats (pSTRs) in each protein sequence and assess their variability within the family in terms of number of units. To facilitate this analysis, we developed the pSTR tool, a method to analyze the evolution of protein short tandem repeats in a given protein family by pairwise comparisons between evolutionarily related protein sequences. We evaluated pSTR unit number variation in protein families of 12 complete metazoan proteomes. We hypothesize that families with more dynamic ensembles of repeats could reflect particular roles of these repeats in processes that require more adaptability.

Structural and biophysical characterization of the secreted, β-helical adhesin EtpA of Enterotoxigenic Escherichia coli

Enterotoxigenic Escherichia coli (ETEC) is a diarrhoeal pathogen associated with high morbidity and mortality especially among young children in developing countries. At present, there is no vaccine for ETEC. One candidate vaccine antigen, EtpA, is a conserved secreted adhesin that binds to the tips of flagellae to bridge ETEC to host intestinal glycans. EtpA is exported through a Gram-negative, two-partner secretion system (TPSS, type Vb) comprised of the secreted EtpA passenger (TpsA) protein and EtpB (TpsB) transporter that is integrated into the outer bacterial membrane. TpsA proteins share a conserved, N-terminal TPS domain followed by an extensive C-terminal domain with divergent sequence repeats. Two soluble, N-terminal constructs of EtpA were prepared and analysed respectively including residues 67 to 447 (EtpA67-447) and 1 to 606 (EtpA1-606). The crystal structure of EtpA67-447 solved at 1.76 Å resolution revealed a right-handed parallel β-helix with two extra-helical hairpins and an N-terminal β-strand cap. Analyses by circular dichroism spectroscopy confirmed the β-helical fold and indicated high resistance to chemical and thermal denaturation as well as rapid refolding. A theoretical AlphaFold model of full-length EtpA largely concurs with the crystal structure adding an extended β-helical C-terminal domain after an interdomain kink. We propose that robust folding of the TPS domain upon secretion provides a template to extend the N-terminal β-helix into the C-terminal domains of TpsA proteins.

Functional annotation of uncharacterized proteins from Fusobacterium nucleatum: identification of virulence factors

Fusobacterium nucleatum is a gram-negative bacteria associated with diverse infections like appendicitis and colorectal cancer. It mainly attacks the epithelial cells in the oral cavity and throat of the infected individual. It has a single circular genome of 2.7 Mb. Many proteins in F. nucleatum genome are listed as "Uncharacterized." Annotation of these proteins is crucial for obtaining new facts about the pathogen and deciphering the gene regulation, functions, and pathways along with discovery of novel target proteins. In the light of new genomic information, an armoury of bioinformatic tools were used for predicting the physicochemical parameters, domain and motif search, pattern search, and localization of the uncharacterized proteins. The programs such as receiver operating characteristics determine the efficacy of the databases that have been employed for prediction of different parameters at 83.6%. Functions were successfully assigned to 46 uncharacterized proteins which included enzymes, transporter proteins, membrane proteins, binding proteins, etc. Apart from the function prediction, the proteins were also subjected to string analysis to reveal the interacting partners. The annotated proteins were also put through homology-based structure prediction and modeling using Swiss PDB and Phyre2 servers. Two probable virulent factors were also identified which could be investigated further for potential drug-related studies. The assigning of functions to uncharacterized proteins has shown that some of these proteins are important for cell survival inside the host and can act as effective drug targets.

Protein repeats evolve and emerge in giant viruses

Nucleocytoplasmatic large DNA viruses (NCLDVs or giant viruses) stand out because of their relatively large genomes encoding hundreds of proteins. These species give us an unprecedented opportunity to study the emergence and evolution of repeats in protein sequences. On the one hand, as viruses, these species have a restricted set of functions, which can help us better define the functional landscape of repeats. On the other hand, given the particular use of the genetic machinery of the host, it is worth asking whether this allows the variations of genetic material that lead to repeats in non-viral species. To support research in the characterization of repeat protein evolution and function, we present here an analysis focused on the repeat proteins of giant viruses, namely tandem repeats (TRs), short repeats (SRs), and homorepeats (polyX). Proteins with large and short repeats are not very frequent in non-eukaryotic organisms because of the difficulties that their folding may entail; however, their presence in giant viruses remarks their advantage for performance in the protein environment of the eukaryotic host. The heterogeneous content of these TRs, SRs and polyX in some viruses hints at diverse needs. Comparisons to homologs suggest that the mechanisms that generate these repeats are extensively used by some of these viruses, but also their capacity to adopt genes with repeats. Giant viruses could be very good models for the study of the emergence and evolution of protein repeats.

Exon shuffling potentiates a diverse repertoire of brown algal NB-ARC-TPR candidate immune receptor proteins via alternative splicing

Like other organisms, brown algae are subject to diseases caused by bacteria, fungi, and viruses. Brown algal immunity mechanisms are not well characterized; however, there is evidence suggesting that pathogen receptors exist in brown algae. One key protein family likely associated with brown algal innate immunity possesses an NB-ARC domain analogous to innate immune proteins in plants and animals. In this study, we conducted an extensive survey of NB-ARC genes in brown algae and obtained insights into the domain organization and evolutionary history of the encoded proteins. Our data show that brown algae possess an ancient NB-ARC-tetratricopeptide repeat (NB-TPR) domain architecture. We identified an N-terminal effector domain, the four-helix bundle, which was not previously found associated with NB-ARC domains. The phylogenetic tree including NB-ARC domains from all kingdoms of life suggests the three clades of brown algal NB-TPRs are likely monophyletic, whereas their TPRs seem to have distinct origins. One group of TPRs exhibit intense exon shuffling, with various alternative splicing and diversifying selection acting on them, suggesting exon shuffling is an important mechanism for evolving ligand-binding specificities. The reconciliation of gene duplication and loss events of the NB-ARC genes reveals that more independent gene gains than losses have occurred during brown algal evolution, and that tandem duplication has played a major role in the expansion of NB-ARC genes. Our results substantially enhance our understanding of the evolutionary history and exon shuffling mechanisms of the candidate innate immune repertoire of brown algae.

Native functions of short tandem repeats

Over a third of the human genome is comprised of repetitive sequences, including more than a million short tandem repeats (STRs). While studies of the pathologic consequences of repeat expansions that cause syndromic human diseases are extensive, the potential native functions of STRs are often ignored. Here, we summarize a growing body of research into the normal biological functions for repetitive elements across the genome, with a particular focus on the roles of STRs in regulating gene expression. We propose reconceptualizing the pathogenic consequences of repeat expansions as aberrancies in normal gene regulation. From this altered viewpoint, we predict that future work will reveal broader roles for STRs in neuronal function and as risk alleles for more common human neurological diseases.

Studying the connection between SF3B1 and four types of cancer by analyzing networks constructed based on published research

Splicing factor 3B subunit 1 (SF3B1) is the largest component of SF3b protein complex which is involved in the pre-mRNA splicing mechanism. Somatic mutations of SF3B1 were shown to be associated with aberrant splicing, producing abnormal transcripts that drive cancer development and/or prognosis. In this study, we focus on the relationship between SF3B1 and four types of cancer, namely myelodysplastic syndrome (MDS), acute myeloid leukemia (AML), and chronic lymphocytic leukemia (CLL) and breast cancer (BC). For this purpose, we identified from the Pubmed library only articles which mentioned SF3B1 in connection with the investigated types of cancer for the period 2007 to 2018 to reveal how the connection has developed over time. We left out all published articles which mentioned SF3B1 in other contexts. We retrieved the target articles and investigated the association between SF3B1 and the mentioned four types of cancer. For this we utilized some of the publicly available databases to retrieve gene/variant/disease information related to SF3B1. We used the outcome to derive and analyze a variety of complex networks that reflect the correlation between the considered diseases and variants associated with SF3B1. The results achieved based on the analyzed articles and reported in this article illustrated that SF3B1 is associated with hematologic malignancies, such as MDS, AML, and CLL more than BC. We found that different gene networks may be required for investigating the impact of mutant splicing factors on cancer development based on the target cancer type. Additionally, based on the literature analyzed in this study, we highlighted and summarized what other researchers have reported as the set of genes and cellular pathways that are affected by aberrant splicing in cancerous cells.

Crystal structure of a polyglycine hydrolase determined using a RoseTTAFold model

Polyglycine hydrolases (PGHs) are secreted fungal proteases that cleave the polyglycine linker of Zea mays ChitA, a defensive chitinase, thus overcoming one mechanism of plant resistance to infection. Despite their importance in agriculture, there has been no previous structural characterization of this family of proteases. The objective of this research was to investigate the proteolytic mechanism and other characteristics by structural and biochemical means. Here, the first atomic structure of a polyglycine hydrolase was identified. It was solved by X-ray crystallography using a RoseTTAFold model, taking advantage of recent technical advances in structure prediction. PGHs are composed of two domains: the N- and C-domains. The N-domain is a novel tertiary fold with an as-yet unknown function that is found across all kingdoms of life. The C-domain shares structural similarities with class C β-lactamases, including a common catalytic nucleophilic serine. In addition to insights into the PGH family and its relationship to β-lactamases, the results demonstrate the power of complementing experimental structure determination with new computational techniques.

Identification, Classification and Characterization Analysis of FBXL Gene in Cotton

F-box/LR (FBXL), Leucine-rich repeats in F-box proteins, belongs to the Skp1-Cullin1-F-box protein (SCF) E3 ligase family. FBXL genes play important roles in plant growth, such as plant hormones, responses to environmental stress, and floral organ development. Here, a total of 518 FBXL genes were identified and analyzed in six plant species. Phylogenetic analysis showed that AtFBXLs, VvFBXLs, and GrFBXLs were clustered into three subfamilies (Ⅰ-Ⅲ). Based on the composition of the F-box domain and carboxyl-terminal amino acid sequence, FBXL proteins were classified into three types (Type-A/-B/-C). Whole-genome duplication (WGD) along with tandem duplications and segmental contributed to the expansion of this gene family. The result indicates that four cotton species are also divided into three subfamilies. FBXLs in cotton were classified into three clades by phylogenetic and structural analyses. Furthermore, expression analyses indicated that the expression patterns of GhFBXLs in different cotton tissues were different. The highly expressed of GH_A07G2363 in 5-8 mm anthers, indicates that this gene might play a role in the reproductive process, providing candidate genes for future studies on cotton fertility materials. This study provides an original functional opinion and a useful interpretation of the FBXL protein family in cotton.

The Difference in Structural States between Canonical Proteins and Their Isoforms Established by Proteome-Wide Bioinformatics Analysis

Alternative splicing is an important means of generating the protein diversity necessary for cellular functions. Hence, there is a growing interest in assessing the structural and functional impact of alternative protein isoforms. Typically, experimental studies are used to determine the structures of the canonical proteins ignoring the other isoforms. Therefore, there is still a large gap between abundant sequence information and meager structural data on these isoforms. During the last decade, significant progress has been achieved in the development of bioinformatics tools for structural and functional annotations of proteins. Moreover, the appearance of the AlphaFold program opened up the possibility to model a large number of high-confidence structures of the isoforms. In this study, using state-of-the-art tools, we performed in silico analysis of 58 eukaryotic proteomes. The evaluated structural states included structured domains, intrinsically disordered regions, aggregation-prone regions, and tandem repeats. Among other things, we found that the isoforms have fewer signal peptides, transmembrane regions, or tandem repeat regions in comparison with their canonical counterparts. This could change protein function and/or cellular localization. The AlphaFold modeling demonstrated that frequently isoforms, having differences with the canonical sequences, still can fold in similar structures though with significant structural rearrangements which can lead to changes of their functions. Based on the modeling, we suggested classification of the structural differences between canonical proteins and isoforms. Altogether, we can conclude that a majority of isoforms, similarly to the canonical proteins are under selective pressure for the functional roles.

Lineage-specific protein repeat expansions and contractions reveal malleable regions of immune genes

Functional diversification, a higher evolutionary rate, and intense positive selection help a limited number of immune genes interact with many pathogens. Repeats in protein-coding regions are a well-known source of functional diversification, adaptive variation, and evolutionary novelty in a short time. Repeats play a crucial role in biochemical functions like functional diversification of transcription regulation, protein kinases, cell adhesion, signaling pathways, morphogenesis, DNA repair, recombination, and RNA processing. Repeat length variation can change the associated protein's interaction, efficacy, and overall protein network. Repeats have an intrinsic unstable nature and can potentially evolve rapidly and expedite the acquisition of complex phenotypic traits and functions. Because of their ability to generate rapid, adaptive variations over short evolutionary distances, repeats are considered "tuning knobs." Repeat length variation in specific genes, like RUNX2 and ALX4, is associated with morphological and physiological changes across vertebrates. Here we study repeat length variation as a potent source of species-specific immune diversification across several clades of tetrapods. Moreover, we provide a clade-wise comprehensive list of immune genes with repeat types for future studies of morphological/evolutionary changes within species groups. We observe significant repeat length variation of FASLG and C1QC in Rodentia and Primates' contrasting species groups, respectively.

The sequence context in poly-alanine regions: structure, function and conservation

MOTIVATION

Poly-alanine (polyA) regions are protein stretches mostly composed of alanines. Despite their abundance in eukaryotic proteomes and their association to nine inherited human diseases, the structural and functional roles exerted by polyA stretches remain poorly understood. In this work we study how the amino acid context in which polyA regions are settled in proteins influences their structure and function.

RESULTS

We identified glycine and proline as the most abundant amino acids within polyA and in the flanking regions of polyA tracts, in human proteins as well as in 17 additional eukaryotic species. Our analyses indicate that the non-structuring nature of these two amino acids influences the α-helical conformations predicted for polyA, suggesting a relevant role in reducing the inherent aggregation propensity of long polyA. Then, we show how polyA position in protein N-termini relates with their function as transit peptides. PolyA placed just after the initial methionine is often predicted as part of mitochondrial transit peptides, whereas when placed in downstream positions, polyA are part of signal peptides. A few examples from known structures suggest that short polyA can emerge by alanine substitutions in α-helices; but evolution by insertion is observed for longer polyA. Our results showcase the importance of studying the sequence context of homorepeats as a mechanism to shape their structure-function relationships.

AVAILABILITY AND IMPLEMENTATION

The datasets used and/or analyzed during the current study are available from the corresponding author onreasonable request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Functional Tuning of Intrinsically Disordered Regions in Human Proteins by Composition Bias

Intrinsically disordered regions (IDRs) in protein sequences are flexible, have low structural constraints and as a result have faster rates of evolution. This lack of evolutionary conservation greatly limits the use of sequence homology for the classification and functional assessment of IDRs, as opposed to globular domains. The study of IDRs requires other properties for their classification and functional prediction. While composition bias is not a necessary property of IDRs, compositionally biased regions (CBRs) have been noted as frequent part of IDRs. We hypothesized that to characterize IDRs, it could be helpful to study their overlap with particular types of CBRs. Here, we evaluate this overlap in the human proteome. A total of 2/3 of residues in IDRs overlap CBRs. Considering CBRs enriched in one type of amino acid, we can distinguish CBRs that tend to be fully included within long IDRs (R, H, N, D, P, G), from those that partially overlap shorter IDRs (S, E, K, T), and others that tend to overlap IDR terminals (Q, A). CBRs overlap more often IDRs in nuclear proteins and in proteins involved in liquid-liquid phase separation (LLPS). Study of protein interaction networks reveals the enrichment of CBRs in IDRs by tandem repetition of short linear motifs (rich in S or P), and the existence of E-rich polar regions that could support specific protein interactions with non-specific interactions. Our results open ways to pin down the function of IDRs from their partial compositional biases.

The PCI domains are “winged” HEAT domains

The HEAT domains are a family of helical hairpin repeat domains, composed of four or more hairpins. HEAT is derived from the names of four family members: huntingtin, eukaryotic translation elongation factor 3 (eEF3), protein phosphatase 2 regulatory A subunit (PP2A), and mechanistic target of rapamycin (mTOR). HEAT domain-containing proteins play roles in a wide range of cellular processes, such as protein synthesis, nuclear transport and metabolism, and cell signaling. The PCI domains are a related group of helical hairpin domains, with a “winged-helix” (WH) subdomain at their C-terminus, which is responsible for multi-subunit complex formation with other PCI domains. The name is derived from the complexes, where these domains are found: the 26S Proteasome “lid” regulatory subcomplex, the COP9 signalosome (CSN), and eukaryotic translation initiation factor 3 (eIF3). We noted that in structure similarity searches using HEAT domains, sometimes PCI domains appeared in the search results ahead of other HEAT domains, which indicated that the PCI domains could be members of the HEAT domain family, and not a related but separate group, as currently thought. Here, we report extensive structure similarity analysis of HEAT and PCI domains, both within and between the two groups of proteins. We present evidence that the PCI domains as a group have greater structural similarity with individual groups of HEAT domains than some of the HEAT domain groups have among each other. Therefore, our results indicate that the PCI domains have evolved from a HEAT domain that acquired a WH subdomain. The WH subdomain in turn mediated self-association into a multi-subunit complex, which eventually evolved into the common ancestor of the Proteasome lid/CSN/eIF3.

Identification of Feronia-interacting proteins in Arabidopsis thaliana

BACKGROUND

Plant growth and development are complex processes modulated by numerous genes, transcription factors, hormones, and peptides. Several reports implicate the membrane-localized Catharanthus roseus receptor-like kinase1 (CrRLK1L) protein, FERONIA (FER), involved in plant development. However, protein targets of FER remain poorly characterized.

OBJECTIVE

FER recombinant proteins were analyzed, and FER-interacting proteins were identified, to better understand the function of the Arabidopsis thaliana FER (AtFER) gene in plant development.

METHODS

AtFER-interacting proteins were identified through Yeast-Two Hybrid (Y2H) and validated by bimolecular fluorescence complementation (BiFC). Autophosphorylation activity was evaluated in AtFER site-directed and deletion mutants.

RESULTS

AtFER cytoplasmic kinase domain (Flag-FER-CD) is autophosphorylated at the Thr residue (s), with T559 and T664 as important sites for AtFER kinase activity. In addition, the carboxy terminal region is essential for AtFER kinase activity. Y2H identified an Armadillo (ARM)-repeat protein (At4g16490) with tandem copies of a degenerate protein sequence motif, a U-BOX 9 (PUB9, At3g07360), IQ-DOMAIN 7 (IQD7, At1g17480), and heteroglycan glucosidase 1 (HGL1, At3g23640) as AtFER-interacting proteins. BiFC confirmed the in vivo interactions between these four proteins and AtFER in tobacco (Nicotiana benthamiana) leaf transient expression assays. The RAPID ALKALINIZATION FACTOR1 (RALF1) peptide, which is a FER ligand, induced the expression of genes encoding the four AtFER-interacting proteins.

CONCLUSION

The AtFER-interacting proteins identified in this study are likely involved in FER-mediated intracellular signaling pathways that are essential in plant growth and development, and possibly plant immunity.

Role of an FNIP Repeat Domain-Containing Protein Encoded by Megavirus Baoshan during Viral Infection

FNIP repeat domain-containing protein (FNIP protein) is a little-studied atypical leucine-rich repeat domain-containing protein found in social amoebae and mimiviruses. Here, a recently reported mimivirus of lineage C, Megavirus baoshan, was analyzed for FNIP protein genes. A total of 82 FNIP protein genes were identified, each containing up to 26 copies of the FNIP repeat, and mostly having an F-box domain at the N terminus. Both nucleotide and amino acid sequences of FNIP repeat were highly conserved. Most of the FNIP protein genes clustered together tandemly in groups of two to 14 genes. Nearly all FNIP protein genes shared similar expression patterns and were expressed 4 to 9 h postinfection. A typical viral FNIP protein, Mb0983, was selected for functional analysis. Protein interactome analysis identified two small GTPases, Rap1B and Rab7A, that interacted with Mb0983 in cytoplasm. The overexpression of Mb0983 in Acanthamoeba castellanii accelerated the degradation of Rap1B and Rab7A during viral infection. Mb0983 also interacted with host SKP1 and cullin-1, which were conserved components of the SKP1-cullin-1-F-box protein (SCF)-type ubiquitin E3 ligase complex. Deletion of the F-box domain of Mb0983 not only abolished its interaction with SKP1 and cullin-1 but also returned the speed of Rap1B and Rab7A degradation to normal in infected A. castellanii. These results suggested that Mb0983 is a part of the SCF-type ubiquitin E3 ligase complex and plays a role in the degradation of Rap1B and Rab7A. They also implied that other viral F-box-containing FNIP proteins might have similar effects on various host proteins. IMPORTANCE Megavirus baoshan encodes 82 FNIP proteins, more than any other reported mimiviruses. Their genetic and transcriptional features suggest that they are important for virus infection and adaption. Since most mimiviral FNIP proteins have the F-box domain, they were predicted to be involved in protein ubiquitylation. FNIP protein Mb0983 interacted with host SKP1 and cullin-1 through the F-box domain, supporting the idea that it is a part of the SCF-type ubiquitin E3 ligase complex. The substrates of Mb0983 for degradation were identified as the host small GTPases Rap1B and Rab7A. Combining the facts of the presence of a large number of FNIP genes in megavirus genomes, the extremely high expression level of the viral ubiquitin gene, and the reported observation that 35% of megavirus-infected amoeba cells died without productive infection, it is likely that megavirus actively explores the host ubiquitin-proteasome pathway in infection and that viral FNIP proteins play roles in the process.

Genome-Wide Analysis of Anthocyanin Biosynthesis Regulatory WD40 Gene FcTTG1 and Related Family in Ficus carica L

WD40 proteins serve as crucial regulators in a broad spectrum of plant developmental and physiological processes, including anthocyanin biosynthesis. However, in fig (Ficus carica L.), neither the WD40 family nor any member involved in anthocyanin biosynthesis has been elucidated. In the present study, 204 WD40 genes were identified from the fig genome and phylogenetically classified into 5 clusters and 12 subfamilies. Bioinformatics analysis prediction localized 109, 69, and 26 FcWD40 proteins to the cytoplasm, nucleus and other cellular compartments, respectively. RNA-seq data mining revealed 127 FcWD40s expressed at FPKM > 10 in fig fruit. Most of these genes demonstrated higher expression in the early stages of fruit development. FcWD40-97 was recruited according to three criteria: high expression in fig fruit, predicted nuclear localization, and closest clustering with TTG1s identified in other plants. FcWD40-97, encoding 339 amino acids including 5 WD-repeat motifs, showed 88.01 and 87.94% amino acid sequence similarity to apple and peach TTG1, respectively. The gene is located on fig chromosome 4, and is composed of 1 intron and 2 exons. Promoter analysis revealed multiple light-responsive elements, one salicylic acid-responsive element, three methyl jasmonate-responsive elements, and one MYB-binding site involved in flavonoid biosynthesis gene regulation. FcWD40-97 was in the FPKM > 100 expression level group in fig fruit, and higher expression was consistently found in the peel compared to the flesh at the same development stages. Expression level did not change significantly under light deprivation, whereas in leaves and roots, its expression was relatively low. Transient expression verified FcWD40-97's localization to the nucleus. Yeast two-hybrid (Y2H) and biomolecular fluorescence complementation (BiFC) assays revealed that FcWD40-97 interacts with FcMYB114, FcMYB123, and FcbHLH42 proteins in vitro and in vivo, showing that FcWD40-97 functions as a member of the MYB-bHLH-WD40 (MBW) complex in anthocyanin-biosynthesis regulation in fig. We therefore renamed FcWD40-97 as FcTTG1. Our results provide the first systematic analysis of the FcWD40 family and identification of FcTTG1 in fig pigmentation.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

New β-Propellers Are Continuously Amplified From Single Blades in all Major Lineages of the β-Propeller Superfamily

β-Propellers are toroidal folds, in which consecutive supersecondary structure units of four anti-parallel β-strands-called blades-are arranged radially around a central axis. Uniquely among toroidal folds, blades span the full range of sequence symmetry, from near identity to complete divergence, indicating an ongoing process of amplification and differentiation. We have proposed that the major lineages of β-propellers arose through this mechanism and that therefore their last common ancestor was a single blade, not a fully formed β-propeller. Here we show that this process of amplification and differentiation is also widespread within individual lineages, yielding β-propellers with blades of more than 60% pairwise sequence identity in most major β-propeller families. In some cases, the blades are nearly identical, indicating a very recent amplification event, but even in cases where such recently amplified β-propellers have more than 80% overall sequence identity to each other, comparison of their DNA sequence shows that the amplification occurred independently.

Structure and Protein-Protein Interactions of Ice Nucleation Proteins Drive Their Activity

Microbially-produced ice nucleating proteins (INpro) are unique molecular structures with the highest known catalytic efficiency for ice formation. Airborne microorganisms utilize these proteins to enhance their survival by reducing their atmospheric residence times. INpro also have critical environmental effects including impacts on the atmospheric water cycle, through their role in cloud and precipitation formation, as well as frost damage on crops. INpro are ubiquitously present in the atmosphere where they are emitted from diverse terrestrial and marine environments. Even though bacterial genes encoding INpro have been discovered and sequenced decades ago, the details of how the INpro molecular structure and oligomerization foster their unique ice-nucleation activity remain elusive. Using machine-learning based software AlphaFold 2 and trRosetta, we obtained and analysed the first ab initio structural models of full length and truncated versions of bacterial INpro. The modeling revealed a novel beta-helix structure of the INpro central repeat domain responsible for ice nucleation activity. This domain consists of repeated stacks of two beta strands connected by two sharp turns. One beta-strand is decorated with a TxT amino acid sequence motif and the other strand has an SxL[T/I] motif. The core formed between the stacked beta helix-pairs is unusually polar and very distinct from previous INpro models. Using synchrotron radiation circular dichroism, we validated the β-strand content of the central repeat domain in the model. Combining the structural model with functional studies of purified recombinant INpro, electron microscopy and modeling, we further demonstrate that the formation of dimers and higher-order oligomers is key to INpro activity. Using computational docking of the new INpro model based on rigid-body algorithms we could reproduce a previously proposed homodimer structure of the INpro CRD with an interface along a highly conserved tyrosine ladder and show that the dimer model agrees with our functional data. The parallel dimer structure creates a surface where the TxT motif of one monomer aligns with the SxL[T/I] motif of the other monomer widening the surface that interacts with water molecules and therefore enhancing the ice nucleation activity. This work presents a major advance in understanding the molecular foundation for bacterial ice-nucleation activity.

The Repeating, Modular Architecture of the HtrA Proteases

A conserved, 26-residue sequence [AA(X₂)[A/G][G/L](X₂)GDV[I/L](X₂)[V/L]NGE(X₁)V(X₆)] and corresponding structure repeating module were identified within the HtrA protease family using a non-redundant set (N = 20) of publicly available structures. While the repeats themselves were far from sequence perfect, they had notable conservation to a statistically significant level. Three or more repetitions were identified within each protein despite being statistically expected to randomly occur only once per 1031 residues. This sequence repeat was associated with a six stranded antiparallel β-barrel module, two of which are present in the core of the structures of the PA clan of serine proteases, while a modified version of this module could be identified in the PDZ-like domains. Automated structural alignment methods had difficulties in superimposing these β-barrels, but the use of a target human HtrA2 structure showed that these modules had an average RMSD across the set of structures of less than 2 Å (mean and median). Our findings support Dayhoff’s hypothesis that complex proteins arose through duplication of simpler peptide motifs and domains.

Genome-wide identification and characterization of F-box family proteins in sweet potato and its expression analysis under abiotic stress

In this study, genome-wide characterization of F-box proteins in sweet potato yielded 243 IbFBX genes, unevenly distributed on the 15 chromosomes of sweet potato. Gene duplication analysis suggested segmental duplication as the principal factor influencing the expansive evolution of IbFBX genes in sweet potato. Phylogenetic analysis clustered F-box proteins in sweet potato, Arabidopsis, and rice into six clades (I-VI). Gene structure analysis of the IbFBX genes revealed that most of the genes within the same clade were highly conserved. Expression profiles of IbFBX family genes in 9 different tissues and under stress conditions revealed that the IbFBXs were highly upregulated or downregulated in response to salt and drought stress, suggesting their significant roles in abiotic stress response and adaptation. Knowledge of the diverse functions and expression patterns of IbFBXs presents a solid theoretical basis for annotating the functions of IbFBXs and further facilitate the molecular breeding of sweet potato.

Numerous variants of leucine rich repeats in proteins from nucleo-cytoplasmic large DNA viruses

Leucine rich repeats (LRRs) occurring in tandem are 20-29 amino acids long. Eleven LRR types have been recognized. Sequence features of LRRs from viruses were investigated using over 600 LRR proteins from 89 species. Directly before, metagenome data of nucleo-cytoplasmic large dsDNA viruses (NCLDVs) have been published; the 2,074 NCLDVs encode 199,021 proteins. From the NCLDVs 547 LRR proteins were identified and 502 were used for analysis. Various variants of known LRR types were identified in viral LRRs. A comprehensive analysis of TpLRR and FNIP that belong to an LRR type was first performed. The repeating unit lengths (RULs) in five types are 19 residues which is the shortest among all LRRs. The RULs of eight LRR types including FNIP are one to five residues shorter than those of the known, corresponding LRR types. The conserved hydrophobic residues such as Leu, Val or Ile in the consensus sequences are frequently substituted by cysteine at one or two positions. Four unique LRR motifs that are different from those identified previously are observed. The present study enhances the previous result. An evolutionary scenario of short or unique LRR was discussed.

Conservation and diversity in allosteric fingerprints of proteins for evolutionary-inspired engineering and design

Hand-in-hand work of physics and evolution delivered protein universe with diversity of forms, sizes, and functions. Pervasiveness and advantageous traits of allostery made it an important component of the protein function regulation, calling for thorough investigation of its structural determinants and evolution. Learning directly from nature, we explored here allosteric communication in several major folds and repeat proteins, including α/β and β-barrels, β-propellers, Ig-like fold, ankyrin and α/β leucine-rich repeat proteins, which provide structural platforms for many different enzymatic and signalling functions. We obtained a picture of conserved allosteric communication characteristic in different fold types, modifications of the structure-driven signalling patterns via sequence-determined divergence to specific functions, as well as emergence and potential diversification of allosteric regulation in multi-domain proteins and oligomeric assemblies. Our observations will be instrumental in facilitating the engineering and de novo design of proteins with allosterically regulated functions, including development of therapeutic biologics. In particular, results described here may guide the identification of the optimal structural platforms (e.g. fold type, size, and oligomerization states) and the types of diversifications/perturbations, such as mutations, effector binding, and order-disorder transition. The tunable allosteric linkage across distant regions can be used as a pivotal component in the design/engineering of modular biological systems beyond the traditional scaffolding function.