Census of exposed aggregation-prone regions in proteomes.

Loss of solubility usually leads to the detrimental elimination of protein function. In some cases, the protein aggregation is also required for beneficial functions. Given the duality of this phenomenon, it remains a fundamental question how natural selection controls the aggregation. The exponential growth of genomic sequence data and recent progress with in silico predictors of the aggregation allows approaching this problem by a large-scale bioinformatics analysis. Most of the aggregation-prone regions are hidden within the 3D structure, rendering them inaccessible for the intermolecular interactions responsible for aggregation. Thus, the most realistic census of the aggregation-prone regions requires crossing aggregation prediction with information about the location of the natively unfolded regions. This allows us to detect so-called 'exposed aggregation-prone regions' (EARs). Here, we analyzed the occurrence and distribution of the EARs in 76 reference proteomes from the three kingdoms of life. For this purpose, we used a bioinformatics pipeline, which provides a consensual result based on several predictors of aggregation. Our analysis revealed a number of new statistically significant correlations about the presence of EARs in different organisms, their dependence on protein length, cellular localizations, co-occurrence with short linear motifs and the level of protein expression. We also obtained a list of proteins with the conserved aggregation-prone sequences for further experimental tests. Insights gained from this work led to a deeper understanding of the relationship between protein evolution and aggregation.

Pathogenic TRIO variants associated with neurodevelopmental disorders perturb the molecular regulation of TRIO and axon pathfinding in vivo.

The RhoGEF TRIO is known to play a major role in neuronal development by controlling actin cytoskeleton remodeling, primarily through the activation of the RAC1 GTPase. Numerous de novo mutations in the TRIO gene have been identified in individuals with neurodevelopmental disorders (NDDs). We have previously established the first phenotype/genotype correlation in TRIO-associated diseases, with striking correlation between the clinical features of the individuals and the opposite modulation of RAC1 activity by TRIO variants targeting different domains. The mutations hyperactivating RAC1 are of particular interest, as they are recurrently found in patients and are associated with a severe form of NDD and macrocephaly, indicating their importance in the etiology of the disease. Yet, it remains unknown how these pathogenic TRIO variants disrupt TRIO activity at a molecular level and how they affect neurodevelopmental processes such as axon outgrowth or guidance. Here we report an additional cohort of individuals carrying a pathogenic TRIO variant that reinforces our initial phenotype/genotype correlation. More importantly, by performing conformation predictions coupled to biochemical validation, we propose a model whereby TRIO is inhibited by an intramolecular fold and NDD-associated variants relieve this inhibition, leading to RAC1 hyperactivation. Moreover, we show that in cultured primary neurons and in the zebrafish developmental model, these gain-of-function variants differentially affect axon outgrowth and branching in vitro and in vivo, as compared to loss-of-function TRIO variants. In summary, by combining clinical, molecular, cellular and in vivo data, we provide compelling new evidence for the pathogenicity of novel genetic variants targeting the TRIO gene in NDDs. We report a novel mechanism whereby the fine-tuned regulation of TRIO activity is critical for proper neuronal development and is disrupted by pathogenic mutations.

Aspartate-phobia of thermophiles as a reaction to deleterious chemical transformations.

Prokaryotes growing at high temperatures have a high proportion of charged residues in their proteins to stabilize their 3D structure. By mining 175 disparate bacterial and archaeal proteomes we found that, against the general trend for charged residues, the frequency of aspartic acid residues decreases strongly as natural growth temperature increases. In search of the explanation, we hypothesized that the reason for such unusual correlation is the deleterious consequences of spontaneous chemical transformations of aspartate at high temperatures. Our subsequent statistical analysis supported this hypothesis. This finding reveals that organisms have likely adapted to high temperatures by minimizing the harmful consequences of spontaneous chemical transformations.

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.

The repetitive structure of DNA clamps: An overlooked protein tandem repeat.

Structured tandem repeats proteins (STRPs) are a specific kind of tandem repeat proteins characterized by a modular and repetitive three-dimensional structure arrangement. The majority of STRPs adopt solenoid structures, but with the increasing availability of experimental structures and high-quality predicted structural models, more STRP folds can be characterized. Here, we describe "Box repeats", an overlooked STRP fold present in the DNA sliding clamp processivity factors, which has eluded classification although structural data has been available since the late 1990s. Each Box repeat is a ß⍺ßßß module of about 60 residues, which forms a class V "beads-on-a-string" type STRP. The number of repeats present in processivity factors is organism dependent. Monomers of PCNA proteins in both Archaea and Eukarya have 4 repeats, while the monomers of bacterial beta-sliding clamps have 6 repeats. This new repeat fold has been added to the RepeatsDB database, which now provides structural annotation for 66 Box repeat proteins belonging to different organisms, including viruses.

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.

Opposite Modulation of RAC1 by Mutations in TRIO Is Associated with Distinct, Domain-Specific Neurodevelopmental Disorders.

The Rho-guanine nucleotide exchange factor (RhoGEF) TRIO acts as a key regulator of neuronal migration, axonal outgrowth, axon guidance, and synaptogenesis by activating the GTPase RAC1 and modulating actin cytoskeleton remodeling. Pathogenic variants in TRIO are associated with neurodevelopmental diseases, including intellectual disability (ID) and autism spectrum disorders (ASD). Here, we report the largest international cohort of 24 individuals with confirmed pathogenic missense or nonsense variants in TRIO. The nonsense mutations are spread along the TRIO sequence, and affected individuals show variable neurodevelopmental phenotypes. In contrast, missense variants cluster into two mutational hotspots in the TRIO sequence, one in the seventh spectrin repeat and one in the RAC1-activating GEFD1. Although all individuals in this cohort present with developmental delay and a neuro-behavioral phenotype, individuals with a pathogenic variant in the seventh spectrin repeat have a more severe ID associated with macrocephaly than do most individuals with GEFD1 variants, who display milder ID and microcephaly. Functional studies show that the spectrin and GEFD1 variants cause a TRIO-mediated hyper- or hypo-activation of RAC1, respectively, and we observe a striking correlation between RAC1 activation levels and the head size of the affected individuals. In addition, truncations in TRIO GEFD1 in the vertebrate model X. tropicalis induce defects that are concordant with the human phenotype. This work demonstrates distinct clinical and molecular disorders clustering in the GEFD1 and seventh spectrin repeat domains and highlights the importance of tight control of TRIO-RAC1 signaling in neuronal development.

RepeatsDB in 2021: improved data and extended classification for protein tandem repeat structures.

The RepeatsDB database (URL: https://repeatsdb.org/) provides annotations and classification for protein tandem repeat structures from the Protein Data Bank (PDB). Protein tandem repeats are ubiquitous in all branches of the tree of life. The accumulation of solved repeat structures provides new possibilities for classification and detection, but also increasing the need for annotation. Here we present RepeatsDB 3.0, which addresses these challenges and presents an extended classification scheme. The major conceptual change compared to the previous version is the hierarchical classification combining top levels based solely on structural similarity (Class > Topology > Fold) with two new levels (Clan > Family) requiring sequence similarity and describing repeat motifs in collaboration with Pfam. Data growth has been addressed with improved mechanisms for browsing the classification hierarchy. A new UniProt-centric view unifies the increasingly frequent annotation of structures from identical or similar sequences. This update of RepeatsDB aligns with our commitment to develop a resource that extracts, organizes and distributes specialized information on tandem repeat protein structures.

TAPASS: Tool for annotation of protein amyloidogenicity in the context of other structural states.

Numerous studies have demonstrated that the propensity of a protein to form amyloids or amorphous aggregates is encoded by its amino acid sequence. This led to the emergence of several computational programs to predict amyloidogenicity from amino acid sequences. However, a growing number of studies indicate that an accurate prediction of the protein aggregation can only be achieved when also accounting for the overall structural context of the protein, and the likelihood of transition between the initial state and the aggregate. Here, we describe a computational pipeline called TAPASS, which was designed to do just that. The pipeline assigns each residue of a protein as belonging to a structured region or an intrinsically disordered region (IDR). For this purpose, TAPASS uses either several state-of-the-art programs for prediction of IDRs, of transmembrane regions and of structured domains or the artificial intelligence program AlphaFold. In the next step, this assignment is crossed with amyloidogenicity prediction. As a result, TAPASS allows the detection of Exposed Amyloidogenic Regions (EARs) located within intrinsically disordered regions (IDRs) and carrying high amyloidogenic potential. TAPASS can substantially improve the prediction of amyloids and be used in proteome-wide analysis to discover new amyloid-forming proteins. Its results, combined with clinical data, can create individual risk profiles for different amyloidoses, opening up new opportunities for personalised medicine. The architecture of the pipeline is designed so that it makes it easy to add new individual predictors as they become available. TAPASS can be used through the web interface (https://bioinfo.crbm.cnrs.fr/index.php?route=tools&tool=32).

Disentangling the complexity of low complexity proteins.

There are multiple definitions for low complexity regions (LCRs) in protein sequences, with all of them broadly considering LCRs as regions with fewer amino acid types compared to an average composition. Following this view, LCRs can also be defined as regions showing composition bias. In this critical review, we focus on the definition of sequence complexity of LCRs and their connection with structure. We present statistics and methodological approaches that measure low complexity (LC) and related sequence properties. Composition bias is often associated with LC and disorder, but repeats, while compositionally biased, might also induce ordered structures. We illustrate this dichotomy, and more generally the overlaps between different properties related to LCRs, using examples. We argue that statistical measures alone cannot capture all structural aspects of LCRs and recommend the combined usage of a variety of predictive tools and measurements. While the methodologies available to study LCRs are already very advanced, we foresee that a more comprehensive annotation of sequences in the databases will enable the improvement of predictions and a better understanding of the evolution and the connection between structure and function of LCRs. This will require the use of standards for the generation and exchange of data describing all aspects of LCRs. SHORT ABSTRACT: There are multiple definitions for low complexity regions (LCRs) in protein sequences. In this critical review, we focus on the definition of sequence complexity of LCRs and their connection with structure. We present statistics and methodological approaches that measure low complexity (LC) and related sequence properties. Composition bias is often associated with LC and disorder, but repeats, while compositionally biased, might also induce ordered structures. We illustrate this dichotomy, plus overlaps between different properties related to LCRs, using examples.

The C-terminal segment of Leishmania major HslU: Toward potential inhibitors of LmHslVU activity.

It is urgent to develop less toxic and more efficient treatments for leishmaniases and trypanosomiases. We explore the possibility to target the parasite mitochondrial HslVU protease, which is essential for growth and has no analogue in the human host. For this, we develop compounds potentially inhibiting the complex assembly by mimicking the C-terminal (C-ter) segment of the ATPase HslU. We previously showed that a dodecapeptide derived from Leishmania major HslU C-ter segment (LmC12-U2, Cpd 1) was able to bind to and activate the digestion of a fluorogenic substrate by LmHslV. Here, we present the study of its structure-activity relationships. By replacing each essential residue with related non-proteinogenic residues, we obtained more potent analogues. In particular, a cyclohexylglycine residue at position 11 (cpd 24) allowed a more than three-fold gain in potency while reducing the size of compound 24 from twelve to six residues (cpd 50) without significant loss of potency, opening the way toward short HslU C-ter peptidomimetics as potential inhibitors of HslV proteolytic function. Finally, conjugates constituted of LmC6-U2 analogues and a mitochondrial penetrating peptide were found to penetrate into the promastigote form of L. infantum and to inhibit the parasite growth without showing toxicity toward human THP-1 cells at the same concentration (i.e. 30 µM).

DisProt: intrinsic protein disorder annotation in 2020.

The Database of Protein Disorder (DisProt, URL: https://disprot.org) provides manually curated annotations of intrinsically disordered proteins from the literature. Here we report recent developments with DisProt (version 8), including the doubling of protein entries, a new disorder ontology, improvements of the annotation format and a completely new website. The website includes a redesigned graphical interface, a better search engine, a clearer API for programmatic access and a new annotation interface that integrates text mining technologies. The new entry format provides a greater flexibility, simplifies maintenance and allows the capture of more information from the literature. The new disorder ontology has been formalized and made interoperable by adopting the OWL format, as well as its structure and term definitions have been improved. The new annotation interface has made the curation process faster and more effective. We recently showed that new DisProt annotations can be effectively used to train and validate disorder predictors. We believe the growth of DisProt will accelerate, contributing to the improvement of function and disorder predictors and therefore to illuminate the 'dark' proteome.

Molecular Determinants of Fibrillation in a Viral Amyloidogenic Domain from Combined Biochemical and Biophysical Studies.

The Nipah and Hendra viruses (NiV and HeV) are biosafety level 4 human pathogens classified within the Henipavirus genus of the Paramyxoviridae family. In both NiV and HeV, the gene encoding the Phosphoprotein (P protein), an essential polymerase cofactor, also encodes the V and W proteins. These three proteins, which share an intrinsically disordered N-terminal domain (NTD) and have unique C-terminal domains (CTD), are all known to counteract the host innate immune response, with V and W acting by either counteracting or inhibiting Interferon (IFN) signaling. Recently, the ability of a short region within the shared NTD (i.e., PNT3) to form amyloid-like structures was reported. Here, we evaluated the relevance of each of three contiguous tyrosine residues located in a previously identified amyloidogenic motif (EYYY) within HeV PNT3 to the fibrillation process. Our results indicate that removal of a single tyrosine in this motif significantly decreases the ability to form fibrils independently of position, mainly affecting the elongation phase. In addition, we show that the C-terminal half of PNT3 has an inhibitory effect on fibril formation that may act as a molecular shield and could thus be a key domain in the regulation of PNT3 fibrillation. Finally, the kinetics of fibril formation for the two PNT3 variants with the highest and the lowest fibrillation propensity were studied by Taylor Dispersion Analysis (TDA). The results herein presented shed light onto the molecular mechanisms involved in fibril formation.

Census of halide-binding sites in protein structures.

MOTIVATION: Halides are negatively charged ions of halogens, forming fluorides (F-), chlorides (Cl-), bromides (Br-) and iodides (I-). These anions are quite reactive and interact both specifically and non-specifically with proteins. Despite their ubiquitous presence and important roles in protein function, little is known about the preferences of halides binding to proteins. To address this problem, we performed the analysis of halide-protein interactions, based on the entries in the Protein Data Bank. RESULTS: We have compiled a pipeline for the quick analysis of halide-binding sites in proteins using the available software. Our analysis revealed that all of halides are strongly attracted by the guanidinium moiety of arginine side chains, however, there are also certain preferences among halides for other partners. Furthermore, there is a certain preference for coordination numbers in the binding sites, with a correlation between coordination numbers and amino acid composition. This pipeline can be used as a tool for the analysis of specific halide-protein interactions and assist phasing experiments relying on halides as anomalous scatters. AVAILABILITY AND IMPLEMENTATION: All data described in this article can be reproduced via complied pipeline published at https://github.com/rostkick/Halide_sites/blob/master/README.md. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Tally-2.0: upgraded validator of tandem repeat detection in protein sequences.

MOTIVATION: Proteins containing tandem repeats (TRs) are abundant, frequently fold in elongated non-globular structures and perform vital functions. A number of computational tools have been developed to detect TRs in protein sequences. A blurred boundary between imperfect TR motifs and non-repetitive sequences gave rise to necessity to validate the detected TRs. RESULTS: Tally-2.0 is a scoring tool based on a machine learning (ML) approach, which allows to validate the results of TR detection. It was upgraded by using improved training datasets and additional ML features. Tally-2.0 performs at a level of 93% sensitivity, 83% specificity and an area under the receiver operating characteristic curve of 95%. AVAILABILITY AND IMPLEMENTATION: Tally-2.0 is available, as a web tool and as a standalone application published under Apache License 2.0, on the URL https://bioinfo.crbm.cnrs.fr/index.php? route=tools&tool=27. It is supported on Linux. Source code is available upon request. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Tandem repeats lead to sequence assembly errors and impose multi-level challenges for genome and protein databases.

The widespread occurrence of repetitive stretches of DNA in genomes of organisms across the tree of life imposes fundamental challenges for sequencing, genome assembly, and automated annotation of genes and proteins. This multi-level problem can lead to errors in genome and protein databases that are often not recognized or acknowledged. As a consequence, end users working with sequences with repetitive regions are faced with 'ready-to-use' deposited data whose trustworthiness is difficult to determine, let alone to quantify. Here, we provide a review of the problems associated with tandem repeat sequences that originate from different stages during the sequencing-assembly-annotation-deposition workflow, and that may proliferate in public database repositories affecting all downstream analyses. As a case study, we provide examples of the Atlantic cod genome, whose sequencing and assembly were hindered by a particularly high prevalence of tandem repeats. We complement this case study with examples from other species, where mis-annotations and sequencing errors have propagated into protein databases. With this review, we aim to raise the awareness level within the community of database users, and alert scientists working in the underlying workflow of database creation that the data they omit or improperly assemble may well contain important biological information valuable to others.

Amyloidogenicity as a driving force for the formation of functional oligomers.

Insoluble amyloid fibrils formed by self-assembly of amyloidogenic regions of proteins have a cross-ß-structure. In this work, by using targeted molecular dynamics and rigid body simulation, we demonstrate that if a protein consists of an amyloidogenic region and a globular domain(s) and if the linker between them is short enough, such molecules cannot assemble into amyloid fibrils, instead, they form oligomers with a defined and limited number of ß-strands in the cross-ß core. We show that this blockage of the amyloid growth is due to the steric repulsion of the globular structures linked to amyloidogenic regions. Furthermore, we establish a relationship between the linker length and the number of monomers in such nanoparticles. We hypothesise that such oligomerisation can be a yet unrecognised way to form natural protein complexes involved in biological processes. Our results can also be used in protein engineering for designing soluble nanoparticles carrying different functional domains.

Multifunctional Amyloid Oligomeric Nanoparticles for Specific Cell Targeting and Drug Delivery.

Natural selection has endorsed proteins with amazing structures and functionalities that cannot be matched by synthetic means, explaining the exponential interest in developing protein-based materials. Protein self-assembly allows fabricating complex supramolecular structures from relatively simple building blocks, a bottom-up strategy naturally employed by amyloid fibrils. However, the design of amyloid-inspired materials with biological activity is inherently difficult. Here, we exploit a modular procedure to generate functional amyloid nanostructures with tight control of their mesoscopic properties. The soft amyloid core of a yeast prion was fused to dihydrofolate reductase through flexible linkers of different sizes. This enabled us to produce, for the first time, biocompatible protein-only amyloid-like oligomeric nanoparticles with defined dimensions in which the embedded enzyme remained highly active, as assessed by biophysical and enzymatic assays. The modular design allowed one to obtain multifunctional nanoparticles by incorporating the antibody-binding Z-domain to the protein fusion. We show how these assemblies can be exploited for antibody-directed targeting of specific cell types and the localized delivery of methotrexate, resulting in the intracellular uptake of the drug by cancer cells and their death. Overall, the novel protein particles we describe in this work might find applications in areas as diverse as biocatalysis, bioimaging, or targeted therapies.

Porins and Amyloids are Coded by Similar Sequence Motifs.

The relationship between amino acid sequences of the ß-hairpin structures and amyloidogenic ß-arcade-forming motifs are of special interest because, similar to amyloid fibrils, most of the ß-hairpin repeat (BHR) structures have the so-called cross-ß arrangement. Moreover, ß-hairpin is considered as a probable intermediate structure in amyloidogenesis. In this work, a bioinformatics sequence analysis of the known BHR structures is performed in search of amylodogenic motifs able to form ß-arcade fibrils. The analysis shows that the occurrence of the predicted ß-arcade motifs in the BHR regions is very different depending on the BHR structural fold, cellular localization, and phylogeny. One of the most striking observations is the high level of sequence similarity between the BHRs of membranous porins and ß-arcade motifs. This sequence similarity provides additional evidence that the structure of the membranous porins and annular amyloid oligomers may bear a resemblance. Moreover, these results explain how some amyloidogenic sequence can fold in either the ring-like shape oligomers or elongated amyloid fibrils. It has been also found that potentially lethal amyloidogenic ß-arcade motifs are absent in the elongated BHR structures of intracellular eukaryotic proteins. It allows to hypothesize that, in this case, the selective evolutionary pressure acts against aggregation.

BetaSerpentine: a bioinformatics tool for reconstruction of amyloid structures.

Motivation: Numerous experimental studies have suggested that polypeptide chains of large amyloidogenic regions zig-zag in ß-serpentine arrangements. These ß-serpentines are stacked axially and form the superpleated ß-structure. Despite this progress in the understanding of amyloid folds, the determination of their 3D structure at the atomic level is still a problem due to the polymorphism of these fibrils and incompleteness of experimental structural data. Today, the way to get insight into the atomic structure of amyloids is a combination of experimental studies with bioinformatics. Results: We developed a computer program BetaSerpentine that reconstructs ß-serpentine arrangements from individual ß-arches predicted by ArchCandy program and ranks them in order of preference. It was shown that the BetaSerpentine program in combination with the experimental data can be used to gain insight into the detailed 3D structure of amyloids. It opens avenues to the structure-based interpretation and design of the experiments. Availability and implementation: BetaSerpentine webserver can be accessed through website: http://bioinfo.montp.cnrs.fr/b-serpentine. Source code is available in git.hub repository (github.com/stanislavspbgu/BetaSerpentine). Contact: stanislavspbgu@gmail.com or andrey.kajava@crbm.cnrs.fr. Supplementary information: Supplementary data are available at Bioinformatics online.