Cross-reactivity of rPvs48/45, a recombinant Plasmodium vivax protein, with sera from Plasmodium falciparum endemic areas of Africa

Background

Ps48/45, a Plasmodium gametocyte surface protein, is a promising candidate for malaria transmission-blocking (TB) vaccine. Due to its relevance for a multispecies vaccine, we explored the cross-reactivity and TB activity of a recombinant P. vivax Ps48/45 protein (rPvs48/45) with sera from P. falciparum-exposed African donors.

Methods

rPvs48/45 was produced in Chinese hamster ovary cell lines and tested by ELISA for its cross-reactivity with sera from Burkina Faso, Tanzania, Mali, and Nigeria - In addition, BALB/c mice were immunized with the rPvs48/45 protein formulated in Montanide ISA-51 and inoculated with a crude extract of P. falciparum NF-54 gametocytes to evaluate the parasite-boosting effect on rPvs48/45 antibody titers. Specific anti-rPvs48/45 IgG purified from African sera was used to evaluate the ex vivo TB activity on P. falciparum, using standard mosquito membrane feeding assays (SMFA).

Results

rPvs48/45 protein showed cross-reactivity with sera of individuals from all four African countries, in proportions ranging from 94% (Tanzania) to 40% (Nigeria). Also, the level of cross-reactive antibodies varied significantly between countries (p<0.0001), with a higher antibody level in Mali and the lowest in Nigeria. In addition, antibody levels were higher in adults (≥ 17 years) than young children (≤ 5 years) in both Mali and Tanzania, with a higher proportion of responders in adults (90%) than in children (61%) (p<0.0001) in Mali, where male (75%) and female (80%) displayed similar antibody responses. Furthermore, immunization of mice with P. falciparum gametocytes boosted anti-Pvs48/45 antibody responses, recognizing P. falciparum gametocytes in indirect immunofluorescence antibody test. Notably, rPvs48/45 affinity-purified African IgG exhibited a TB activity of 61% against P. falciparum in SMFA.

Conclusion

African sera (exposed only to P. falciparum) cross-recognized the rPvs48/45 protein. This, together with the functional activity of IgG, warrants further studies for the potential development of a P. vivax and P. falciparum cross-protective TB vaccine.

AmyloComp: A Bioinformatic Tool for Prediction of Amyloid Co-aggregation

Typically, amyloid fibrils consist of multiple copies of the same protein. In these fibrils, each polypeptide chain adopts the same β-arc-containing conformation and these chains are stacked in a parallel and in-register manner. In the last few years, however, a considerable body of data has been accumulated about co-aggregation of different amyloid-forming proteins. Among known examples of the co-aggregation are heteroaggregates of different yeast prions and human proteins Rip1 and Rip3. Since the co-aggregation is linked to such important phenomena as infectivity of amyloids and molecular mechanisms of functional amyloids, we analyzed its structural aspects in more details. An axial stacking of different proteins within the same amyloid fibril is one of the most common type of co-aggregation. By using an approach based on structural similarity of the growing tips of amyloids, we developed a computational method to predict amyloidogenic β-arch structures that are able to interact with each other by the axial stacking. Furthermore, we compiled a dataset consisting of 26 experimentally known pairs of proteins capable or incapable to co-aggregate. We utilized this dataset to test and refine our algorithm. The developed method opens a way for a number of applications, including the identification of microbial proteins capable triggering amyloidosis in humans. AmyloComp is available on the website: https://bioinfo.crbm.cnrs.fr/index.php?route=tools&tool=30.

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.

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.

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.

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.

ATR-FTIR spectrum analysis of plasma samples for rapid identification of recovered COVID-19 individuals

The development of fast, cheap and reliable methods to determine seroconversion against infectious agents is of great practical importance. In the context of the COVID-19 pandemic, an important issue is to study the rate of formation of the immune layer in the population of different regions, as well as the study of the formation of post-vaccination immunity in individuals after vaccination. Currently, the main method for this kind of research is enzyme immunoassay (ELISA, enzyme-linked immunosorbent assay). This technique is sufficiently sensitive and specific, but it requires significant time and material costs. We investigated the applicability of attenuated total reflection (ATR) Fourier transform infrared (FTIR) spectroscopy associated with machine learning in blood plasma to detect seroconversion against SARS-CoV-2. The study included samples of 60 patients. Clear spectral differences in plasma samples from recovered COVID-19 patients and conditionally healthy donors were identified using multivariate and statistical analysis. The results showed that ATR-FTIR spectroscopy, combined with principal components analysis (PCA) and linear discriminant analysis (LDA) or artificial neural network (ANN), made it possible to efficiently identify specimens from recovered COVID-19 patients. We built classification models based on PCA associated with LDA and ANN. Our analysis led to 87% accuracy for PCA-LDA model and 91% accuracy for ANN, respectively. Based on this proof-of-concept study, we believe this method could offer a simple, label-free, cost-effective tool for detecting seroconversion against SARS-CoV-2. This approach could be used as an alternative to ELISA.

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.

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.

A structural biology community assessment of AlphaFold2 applications

Most proteins fold into 3D structures that determine how they function and orchestrate the biological processes of the cell. Recent developments in computational methods for protein structure predictions have reached the accuracy of experimentally determined models. Although this has been independently verified, the implementation of these methods across structural-biology applications remains to be tested. Here, we evaluate the use of AlphaFold2 (AF2) predictions in the study of characteristic structural elements; the impact of missense variants; function and ligand binding site predictions; modeling of interactions; and modeling of experimental structural data. For 11 proteomes, an average of 25% additional residues can be confidently modeled when compared with homology modeling, identifying structural features rarely seen in the Protein Data Bank. AF2-based predictions of protein disorder and complexes surpass dedicated tools, and AF2 models can be used across diverse applications equally well compared with experimentally determined structures, when the confidence metrics are critically considered. In summary, we find that these advances are likely to have a transformative impact in structural biology and broader life-science research.

A SARS-CoV-2 Spike Receptor Binding Motif Peptide Induces Anti-Spike Antibodies in Mice andIs Recognized by COVID-19 Patients

The currently devastating pandemic of severe acute respiratory syndrome known as coronavirus disease 2019 or COVID-19 is caused by the coronavirus SARS-CoV-2. Both the virus and the disease have been extensively studied worldwide. A trimeric spike (S) protein expressed on the virus outer bilayer leaflet has been identified as a ligand that allows the virus to penetrate human host cells and cause infection. Its receptor-binding domain (RBD) interacts with the angiotensin-converting enzyme 2 (ACE2), the host-cell viral receptor, and is, therefore, the subject of intense research for the development of virus control means, particularly vaccines. In this work, we search for smaller fragments of the S protein able to elicit virus-neutralizing antibodies, suitable for production by peptide synthesis technology. Based on the analysis of available data, we selected a 72 aa long receptor binding motif (RBM436-507) of RBD. We used ELISA to study the antibody response to each of the three antigens (S protein, its RBD domain and the RBM436-507 synthetic peptide) in humans exposed to the infection and in immunized mice. The seroreactivity analysis showed that anti-RBM antibodies are produced in COVID-19 patients and immunized mice and may exert neutralizing function, although with a frequency lower than anti-S and -RBD. These results provide a basis for further studies towards the development of vaccines or treatments focused on specific regions of the S virus protein, which can benefit from the absence of folding problems, conformational constraints and other advantages of the peptide synthesis production.

Seroreactivity of the Severe Acute Respiratory Syndrome Coronavirus 2 Recombinant S Protein, Receptor-Binding Domain, and Its Receptor-Binding Motif in COVID-19 Patients and Their Cross-Reactivity With Pre-COVID-19 Samples From Malaria-Endemic Areas

Despite the global interest and the unprecedented number of scientific studies triggered by the COVID-19 pandemic, few data are available from developing and low-income countries. In these regions, communities live under the threat of various transmissible diseases aside from COVID-19, including malaria. This study aims to determine the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) seroreactivity of antibodies from COVID-19 and pre-COVID-19 samples of individuals in Mali (West Africa). Blood samples from COVID-19 patients (n = 266) at Bamako Dermatology Hospital (HDB) and pre-COVID-19 donors (n = 283) from a previous malaria survey conducted in Dangassa village were tested by ELISA to assess IgG antibodies specific to the full-length spike (S) protein, the receptor-binding domain (RBD), and the receptor-binding motif (RBM436-507). Study participants were categorized by age, gender, treatment duration for COVID-19, and comorbidities. In addition, the cross-seroreactivity of samples from pre-COVID-19, malaria-positive patients against the three antigens was assessed. Recognition of the SARS-CoV-2 proteins by sera from COVID-19 patients was 80.5% for S, 71.1% for RBD, and 31.9% for RBM (p < 0.001). While antibody responses to S and RBD tended to be age-dependent, responses to RBM were not. Responses were not gender-dependent for any of the antigens. Higher antibody levels to S, RBD, and RBM at hospital entry were associated with shorter treatment durations, particularly for RBD (p < 0.01). In contrast, higher body weights negatively influenced the anti-S antibody response, and asthma and diabetes weakened the anti-RBM antibody responses. Although lower, a significant cross-reactive antibody response to S (21.9%), RBD (6.7%), and RBM (8.8%) was detected in the pre-COVID-19 and malaria samples. Cross-reactive antibody responses to RBM were mostly associated (p < 0.01) with the absence of current Plasmodium falciparum infection, warranting further study.

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).

Predicting Protein Conformational Disorder and Disordered Binding Sites

In the last two decades it has become increasingly evident that a large number of proteins adopt either a fully or a partially disordered conformation. Intrinsically disordered proteins are ubiquitous proteins that fulfill essential biological functions while lacking a stable 3D structure. Their conformational heterogeneity is encoded by the amino acid sequence, thereby allowing intrinsically disordered proteins or regions to be recognized based on their sequence properties. The identification of disordered regions facilitates the functional annotation of proteins and is instrumental for delineating boundaries of protein domains amenable to crystallization. This chapter focuses on the methods currently employed for predicting protein disorder and identifying intrinsically disordered binding sites.

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).

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.

Identification of a Region in the Common Amino-terminal Domain of Hendra Virus P, V, and W Proteins Responsible for Phase Transition and Amyloid Formation

Henipaviruses are BSL-4 zoonotic pathogens responsible in humans for severe encephalitis. Their V protein is a key player in the evasion of the host innate immune response. We previously showed that the Henipavirus V proteins consist of a long intrinsically disordered N-terminal domain (NTD) and a β-enriched C-terminal domain (CTD). These terminals are critical for V binding to DDB1, which is a cellular protein that is a component of the ubiquitin ligase E3 complex, as well as binding to MDA5 and LGP2, which are two host sensors of viral RNA. Here, we serendipitously discovered that the Hendra virus V protein undergoes a liquid-to-hydrogel phase transition and identified the V region responsible for this phenomenon. This region, referred to as PNT3 and encompassing residues 200-310, was further investigated using a combination of biophysical and structural approaches. Congo red binding assays, together with negative-staining transmisison electron microscopy (TEM) studies, show that PNT3 forms amyloid-like fibrils. Fibrillation abilities are dramatically reduced in a rationally designed PNT3 variant in which a stretch of three contiguous tyrosines, falling within an amyloidogenic motif, were replaced by three alanines. Worthy to note, Congo red staining experiments provided hints that these amyloid-like fibrils form not only in vitro but also in cellula after transfection or infection. The present results set the stage for further investigations aimed at assessing the functional role of phase separation and fibrillation by the Henipavirus V proteins.

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.

Immunoreactivity of Sera From Low to Moderate Malaria-Endemic Areas Against Plasmodium vivax rPvs48/45 Proteins Produced in Escherichia coli and Chinese Hamster Ovary Systems

P48/45 is a conserved gametocyte antigen involved in Plasmodium parasite fertilization. A recombinant Plasmodium vivax P48/45 (Pvs48/45) protein expressed in Escherichia coli (E. coli) was highly antigenic and immunogenic in experimental animals and elicited specific transmission-blocking (TB) antibodies in a previous pilot study. Here, a similar Pvs48/45 gene was expressed in Chinese Hamster Ovary (CHO) cells and we compared its immunoreactivity with the E. coli product. Specific antibody titers were determined using plasma from Colombian individuals (n=227) living in endemic areas where both P. vivax and P. falciparum are prevalent and from Guatemala (n=54) where P. vivax is highly prevalent. In Colombia, plasma seroprevalence to CHO-rPvs48/45 protein was 46.3%, while for E. coli-rPvs48/45 protein was 36.1% (p<0.001). In Guatemala, the sero prevalence was 24.1% and 14.8% (p<0.001), respectively. Reactivity index (RI) against both proteins showed an age-dependent increase. IgG2 was the predominant subclass and the antibody avidity index evaluated by ELISA ranged between 4-6 mol/L. Ex vivo P. vivax mosquito direct membrane feeding assays (DMFA) performed in presence of study plasmas, displayed significant parasite transmission-blocking (TB), however, there was no direct correlation between antibody titers and oocysts transmission reduction activity (%TRA). Nevertheless, DMFA with CHO rPvs48/45 affinity purified IgG showed a dose response; 90.2% TRA at 100 μg/mL and 71.8% inhibition at 10 μg/mL. In conclusion, the CHO-rPvs48/45 protein was more immunoreactive in most of the malaria endemic places studied, and CHO-rPvs48/45 specific IgG showed functional activity, supporting further testing of the protein vaccine potential.

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.

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.

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.

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.

Modeling polymorphic ventricular tachycardia at rest using patient-specific induced pluripotent stem cell-derived cardiomyocytes

Background

While mutations in the cardiac type 2 ryanodine receptor (RyR2) have been linked to exercise-induced or catecholaminergic polymorphic ventricular tachycardia (CPVT), its association with polymorphic ventricular tachycardia (PMVT) occurring at rest is unclear. We aimed at constructing a patient-specific human-induced pluripotent stem cell (hiPSC) model of PMVT occurring at rest linked to a single point mutation in RyR2.

Methods

Blood samples were obtained from a patient with PMVT at rest due to a heterozygous RyR2-H29D mutation. Patient-specific hiPSCs were generated from the blood samples, and the hiPSC-derived cardiomyocytes (CMs) were generated via directed differentiation. Using CRIPSR/Cas9 technology, isogenic controls were generated by correcting the RyR2-H29D mutation. Using patch-clamp, fluorescent confocal microscopy and video-image-based analysis, the molecular and functional properties of RyR2-H29D hiPSC—CMs and control hiPSC—CMs were compared.

Findings

RyR2-H29D hiPSC—CMs exhibit intracellular sarcoplasmic reticulum (SR) Ca²⁺ leak through RyR2 under physiological pacing. RyR2-H29D enhances the contribution of inositol 1,4,5-trisphosphate receptors to excitation-contraction coupling (ECC) that exacerbates abnormal Ca²⁺ release in RyR2-H29D hiPSC—CMs. RyR2-H29D hiPSC—CMs exhibit shorter action potentials, delayed afterdepolarizations, arrhythmias and aberrant contractile properties compared to isogenic controls. The RyR2-H29D mutation causes post-translational remodeling that is fully reversed with isogenic controls.

Interpretation

To conclude, in a model based on a RyR2 point mutation that is associated with short-coupled PMVT at rest, RyR2-H29D hiPSC—CMs exhibited aberrant intracellular Ca²⁺ homeostasis, shortened action potentials, arrhythmias and abnormal contractile properties.

Funding

French Muscular Dystrophy Association (AFM; project 16,073, MNM2 2012 and 20,225), “Fondation de la Recherche Médicale” (FRM; SPF20130526710), “Institut National pour la Santé et la Recherche Médicale” (INSERM), National Institutes of Health (ARM; R01 HL145473) and New York State Department of Health (NYSTEM C029156).

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.

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.

Estimation of amyloid aggregate sizes with semi-denaturing detergent agarose gel electrophoresis and its limitations

Semi-denaturing detergent agarose gel electrophoresis (SDD-AGE) was proposed by Vitaly V. Kushnirov in the Michael D. Ter-Avanesyan’s laboratory as a method to compare sizes of amyloid aggregates. Currently, this method is widely used for amyloid investigation, but mostly as a qualitative approach. In this work, we assessed the possibilities and limitations of the quantitative analysis of amyloid aggregate size distribution using SDD-AGE results. For this purpose, we used aggregates of two well-characterized yeast amyloid-forming proteins, Sup35 and Rnq1, and developed a protocol to standardize image analysis and process the result. A detailed investigation of factors that may affect the results of SDD-AGE revealed that both the cell lysis method and electrophoresis conditions can substantially affect the estimation of aggregate size. Despite this, quantitative analysis of SDD-AGE results is possible when one needs to estimate and compare the size of aggregates on the same gel, or even in different experiments, if the experimental conditions are tightly controlled and additional standards are used.

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.

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.

Design of a New [PSI +]-No-More Mutation in SUP35 With Strong Inhibitory Effect on the [PSI +] Prion Propagation

A number of [PSI⁺]-no-more (PNM) mutations, eliminating [PSI⁺] prion, were previously described in SUP35. In this study, we designed and analyzed a new PNM mutation based on the parallel in-register β-structure of Sup35 prion fibrils suggested by the known experimental data. In such an arrangement, substitution of non-charged residues by charged ones may destabilize the fibril structure. We introduced Q33K/A34K amino acid substitutions into the Sup35 protein, corresponding allele was called sup35-M0. The mutagenized residues were chosen based on ArchCandy in silico prediction of high inhibitory effect on the amyloidogenic potential of Sup35. The experiments confirmed that Sup35-M0 leads to the elimination of [PSI⁺] with high efficiency. Our data suggested that the elimination of the [PSI⁺] prion is associated with the decreased aggregation properties of the protein. The new mutation can induce the prion with very low efficiency and is able to propagate only weak [PSI⁺] prion variants. We also showed that Sup35-M0 protein co-aggregates with the wild-type Sup35 in vivo. Moreover, our data confirmed the utility of the strategy of substitution of non-charged residues by charged ones to design new mutations to inhibit a prion formation.

Prion soft amyloid core driven self-assembly of globular proteins into bioactive nanofibrils

Amyloids have been exploited to build amazing bioactive materials. In most cases, short synthetic peptides constitute the functional components of such materials. The controlled assembly of globular proteins into active amyloid nanofibrils is still challenging, because the formation of amyloids implies a conformational conversion towards a β-sheet-rich structure, with a concomitant loss of the native fold and the inactivation of the protein. There is, however, a remarkable exception to this rule: yeast prions. They are singular proteins able to switch between a soluble and an amyloid state. In both states, the structure of their globular domains remains essentially intact. The transit between these two conformations is encoded in prion domains (PrDs): long and disordered sequences to which the active globular domains are appended. PrDs are much larger than typical self-assembling peptides. This seriously limits their use for nanotechnological applications. We have recently shown that these domains contain soft amyloid cores (SACs) that suffice to nucleate their self-assembly reaction. Here we genetically fused a model SAC with different globular proteins. We demonstrate that this very short sequence acts as a minimalist PrD, driving the selective and slow assembly of the initially soluble fusion proteins into amyloid fibrils in which the globular proteins retain their native structure and display high activity. Overall, we provide here a novel, modular and straightforward strategy to build active protein-based nanomaterials at a preparative scale.

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.

Post-Translational Modifications and Diastolic Calcium Leak Associated to the Novel RyR2-D3638A Mutation Lead to CPVT in Patient-Specific hiPSC-Derived Cardiomyocytes

Background: Sarcoplasmic reticulum Ca²⁺ leak and post-translational modifications under stress have been implicated in catecholaminergic polymorphic ventricular tachycardia (CPVT), a highly lethal inherited arrhythmogenic disorder. Human induced pluripotent stem cells (hiPSCs) offer a unique opportunity for disease modeling. Objective: The aims were to obtain functional hiPSC-derived cardiomyocytes from a CPVT patient harboring a novel ryanodine receptor (RyR2) mutation and model the syndrome, drug responses and investigate the molecular mechanisms associated to the CPVT syndrome. Methods: Patient-specific cardiomyocytes were generated from a young athletic female diagnosed with CPVT. The contractile, intracellular Ca²⁺ handling and electrophysiological properties as well as the RyR2 macromolecular remodeling were studied. Results: Exercise stress electrocardiography revealed polymorphic ventricular tachycardia when treated with metoprolol and marked improvement with flecainide alone. We found abnormal stress-induced contractile and electrophysiological properties associated with sarcoplasmic reticulum Ca²⁺ leak in CPVT hiPSC-derived cardiomyocytes. We found inadequate response to metoprolol and a potent response of flecainide. Stabilizing RyR2 with a Rycal compound prevents those abnormalities specifically in CPVT hiPSC-derived cardiomyocytes. The RyR2-D3638A mutation is located in the conformational change inducing-central core domain and leads to RyR2 macromolecular remodeling including depletion of PP2A and Calstabin2. Conclusion: We identified a novel RyR2-D3638A mutation causing 3D conformational defects and aberrant biophysical properties associated to RyR2 macromolecular complex post-translational remodeling. The molecular remodeling is for the first time revealed using patient-specific hiPSC-derived cardiomyocytes which may explain the CPVT proband’s resistance. Our study promotes hiPSC-derived cardiomyocytes as a suitable model for disease modeling, testing new therapeutic compounds, personalized medicine and deciphering underlying molecular mechanisms.

Protein Co-Aggregation Related to Amyloids: Methods of Investigation, Diversity, and Classification

Amyloids are unbranched protein fibrils with a characteristic spatial structure. Although the amyloids were first described as protein deposits that are associated with the diseases, today it is becoming clear that these protein fibrils play multiple biological roles that are essential for different organisms, from archaea and bacteria to humans. The appearance of amyloid, first of all, causes changes in the intracellular quantity of the corresponding soluble protein(s), and at the same time the aggregate can include other proteins due to different molecular mechanisms. The co-aggregation may have different consequences even though usually this process leads to the depletion of a functional protein that may be associated with different diseases. The protein co-aggregation that is related to functional amyloids may mediate important biological processes and change of protein functions. In this review, we survey the known examples of the amyloid-related co-aggregation of proteins, discuss their pathogenic and functional roles, and analyze methods of their studies from bacteria and yeast to mammals. Such analysis allow for us to propose the following co-aggregation classes: (i) titration: deposition of soluble proteins on the amyloids formed by their functional partners, with such interactions mediated by a specific binding site; (ii) sequestration: interaction of amyloids with certain proteins lacking a specific binding site; (iii) axial co-aggregation of different proteins within the same amyloid fibril; and, (iv) lateral co-aggregation of amyloid fibrils, each formed by different proteins.

Identification of the Autochaperone Domain in the Type Va Secretion System (T5aSS): Prevalent Feature of Autotransporters with a β-Helical Passenger

Autotransporters (ATs) belong to a family of modular proteins secreted by the Type V, subtype a, secretion system (T5aSS) and considered as an important source of virulence factors in lipopolysaccharidic diderm bacteria (archetypical Gram-negative bacteria). While exported by the Sec pathway, the ATs are further secreted across the outer membrane via their own C-terminal translocator forming a β-barrel, through which the rest of the protein, namely the passenger, can pass. In several ATs, an autochaperone domain (AC) present at the C-terminal region of the passenger and upstream of the translocator was demonstrated as strictly required for proper secretion and folding. However, considering it was functionally characterised and identified only in a handful of ATs, wariness recently fells on the commonality and conservation of this structural element in the T5aSS. To circumvent the issue of sequence divergence and taking advantage of the resolved three-dimensional structure of some ACs, identification of this domain was performed following structural alignment among all AT passengers experimentally resolved by crystallography before searching in a dataset of 1523 ATs. While demonstrating that the AC is indeed a conserved structure found in numerous ATs, phylogenetic analysis further revealed a distribution into deeply rooted branches, from which emerge 20 main clusters. Sequence analysis revealed that an AC could be identified in the large majority of SAATs (self-associating ATs) but not in any LEATs (lipase/esterase ATs) nor in some PATs (protease autotransporters) and PHATs (phosphatase/hydrolase ATs). Structural analysis indicated that an AC was present in passengers exhibiting single-stranded right-handed parallel β-helix, whatever the type of β-solenoid, but not with α-helical globular fold. From this investigation, the AC of type 1 appears as a prevalent and conserved structural element exclusively associated to β-helical AT passenger and should promote further studies about the protein secretion and folding via the T5aSS, especially toward α-helical AT passengers.

Classification of β-hairpin repeat proteins

In recent years, a number of new protein structures that possess tandem repeats have emerged. Many of these proteins are comprised of tandem arrays of β-hairpins. Today, the amount and variety of the data on these β-hairpin repeat (BHR) structures have reached a level that requires detailed analysis and further classification. In this paper, we classified the BHR proteins, compared structures, sequences of repeat motifs, functions and distribution across the major taxonomic kingdoms of life and within organisms. As a result, we identified six different BHR folds in tandem repeat proteins of Class III (elongated structures) and one BHR fold (up-and-down β-barrel) in Class IV ("closed" structures). Our survey reveals the high incidence of the BHR proteins among bacteria and viruses and their possible relationship to the structures of amyloid fibrils. It indicates that BHR folds will be an attractive target for future structural studies, especially in the context of age-related amyloidosis and emerging infectious diseases. This work allowed us to update the RepeatsDB database, which contains annotated tandem repeat protein structures and to construct sequence profiles based on BHR structural alignments.

Usage of a dataset of NMR resolved protein structures to test aggregation versus solubility prediction algorithms

There has been an increased interest in computational methods for amyloid and (or) aggregate prediction, due to the prevalence of these aggregates in numerous diseases and their recently discovered functional importance. To evaluate these methods, several datasets have been compiled. Typically, aggregation-prone regions of proteins, which form aggregates or amyloids in vivo, are more than 15 residues long and intrinsically disordered. However, the number of such experimentally established amyloid forming and non-forming sequences are limited, not exceeding one hundred entries in existing databases. In this work, we parsed all available NMR-resolved protein structures from the PDB and assembled a new, sevenfold larger, dataset of unfolded sequences, soluble at high concentrations. We proposed to use these sequences as a negative set for evaluating methods for predicting aggregation in vivo. We also present the results of benchmarking cutting edge tools for the prediction of aggregation versus solubility propensity.

RepeatsDB 2.0: improved annotation, classification, search and visualization of repeat protein structures

RepeatsDB 2.0 (URL: http://repeatsdb.bio.unipd.it/) is an update of the database of annotated tandem repeat protein structures. Repeat proteins are a widespread class of non-globular proteins carrying heterogeneous functions involved in several diseases. Here we provide a new version of RepeatsDB with an improved classification schema including high quality annotations for ∼5400 protein structures. RepeatsDB 2.0 features information on start and end positions for the repeat regions and units for all entries. The extensive growth of repeat unit characterization was possible by applying the novel ReUPred annotation method over the entire Protein Data Bank, with data quality is guaranteed by an extensive manual validation for >60% of the entries. The updated web interface includes a new search engine for complex queries and a fully re-designed entry page for a better overview of structural data. It is now possible to compare unit positions, together with secondary structure, fold information and Pfam domains. Moreover, a new classification level has been introduced on top of the existing scheme as an independent layer for sequence similarity relationships at 40%, 60% and 90% identity.

DisProt 7.0: a major update of the database of disordered proteins

The Database of Protein Disorder (DisProt, URL: www.disprot.org) has been significantly updated and upgraded since its last major renewal in 2007. The current release holds information on more than 800 entries of IDPs/IDRs, i.e. intrinsically disordered proteins or regions that exist and function without a well-defined three-dimensional structure. We have re-curated previous entries to purge DisProt from conflicting cases, and also upgraded the functional classification scheme to reflect continuous advance in the field in the past 10 years or so. We define IDPs as proteins that are disordered along their entire sequence, i.e. entirely lack structural elements, and IDRs as regions that are at least five consecutive residues without well-defined structure. We base our assessment of disorder strictly on experimental evidence, such as X-ray crystallography and nuclear magnetic resonance (primary techniques) and a broad range of other experimental approaches (secondary techniques). Confident and ambiguous annotations are highlighted separately. DisProt 7.0 presents classified knowledge regarding the experimental characterization and functional annotations of IDPs/IDRs, and is intended to provide an invaluable resource for the research community for a better understanding structural disorder and for developing better computational tools for studying disordered proteins.

Tally: a scoring tool for boundary determination between repetitive and non-repetitive protein sequences

MOTIVATION

Tandem Repeats (TRs) are abundant in proteins, having a variety of fundamental functions. In many cases, evolution has blurred their repetitive patterns. This leads to the problem of distinguishing between sequences that contain highly imperfect TRs, and the sequences without TRs. The 3D structure of proteins can be used as a benchmarking criterion for TR detection in sequences, because the vast majority of proteins having TRs in sequences are built of repetitive 3D structural blocks. According to our benchmark, none of the existing scoring methods are able to clearly distinguish, based on the sequence analysis, between structures with and without 3D TRs.

RESULTS

We developed a scoring tool called Tally, which is based on a machine learning approach. Tally is able to achieve a better separation between sequences with structural TRs and sequences of aperiodic structures, than existing scoring procedures. It performs at a level of 81% sensitivity, while achieving a high specificity of 74% and an Area Under the Receiver Operating Characteristic Curve of 86%. Tally can be used to select a set of structurally and functionally meaningful TRs from all TRs detected in proteomes. The generated dataset is available for benchmarking purposes.

AVAILABILITY AND IMPLEMENTATION

Source code is available upon request. Tool and dataset can be accessed through our website: http://bioinfo.montp.cnrs.fr/?r=Tally

CONTACT

andrey.kajava@crbm.cnrs.fr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Global genetic diversity of the Plasmodium vivax transmission-blocking vaccine candidate Pvs48/45

Background

Plasmodium vivax 48/45 protein is expressed on the surface of gametocytes/gametes and plays a key role in gamete fusion during fertilization. This protein was recently expressed in Escherichia coli host as a recombinant product that was highly immunogenic in mice and monkeys and induced antibodies with high transmission-blocking activity, suggesting its potential as a P. vivax transmission-blocking vaccine candidate. To determine sequence polymorphism of natural parasite isolates and its potential influence on the protein structure, all pvs48/45 sequences reported in databases from around the world as well as those from low-transmission settings of Latin America were compared.

Methods

Plasmodium vivax parasite isolates from malaria-endemic regions of Colombia, Brazil and Honduras (n = 60) were used to sequence the Pvs48/45 gene, and compared to those previously reported to GenBank and PlasmoDB (n = 222). Pvs48/45 gene haplotypes were analysed to determine the functional significance of genetic variation in protein structure and vaccine potential.

Results

Nine non-synonymous substitutions (E35K, Y196H, H211N, K250N, D335Y, E353Q, A376T, K390T, K418R) and three synonymous substitutions (I73, T149, C156) that define seven different haplotypes were found among the 282 isolates from nine countries when compared with the Sal I reference sequence. Nucleotide diversity (π) was 0.00173 for worldwide samples (range 0.00033–0.00216), resulting in relatively high diversity in Myanmar and Colombia, and low diversity in Mexico, Peru and South Korea. The two most frequent substitutions (E353Q: 41.9 %, K250N: 39.5 %) were predicted to be located in antigenic regions without affecting putative B cell epitopes or the tertiary protein structure.

Conclusions

There is limited sequence polymorphism in pvs48/45 with noted geographical clustering among Asian and American isolates. The low genetic diversity of the protein does not influence the predicted antigenicity or protein structure and, therefore, supports its further development as transmission-blocking vaccine candidate.

α-Synuclein Amyloid Fibrils with Two Entwined, Asymmetrically Associated Protofibrils

Parkinson disease and other progressive neurodegenerative conditions are characterized by the intracerebral presence of Lewy bodies, containing amyloid fibrils of α-synuclein. We used cryo-electron microscopy and scanning transmission electron microscopy (STEM) to study in vitro-assembled fibrils. These fibrils are highly polymorphic. Focusing on twisting fibrils with an inter-crossover spacing of 77 nm, our reconstructions showed them to consist of paired protofibrils. STEM mass per length data gave one subunit per 0.47 nm axial rise per protofibril, consistent with a superpleated β-structure. The STEM images show two thread-like densities running along each of these fibrils, which we interpret as ladders of metal ions. These threads confirmed the two-protofibril architecture of the 77-nm twisting fibrils and allowed us to identify this morphotype in STEM micrographs. Some other, but not all, fibril morphotypes also exhibit dense threads, implying that they also present a putative metal binding site. We propose a molecular model for the protofibril and suggest that polymorphic variant fibrils have different numbers of protofibrils that are associated differently.

Single substitution in bacteriophage T4 RNase H alters the ratio between its exo- and endonuclease activities

The article describes substitutions in bacteriophage T4 RNase H which provide so called das-effect. Phage T4 DNA arrest suppression (das) mutations have been described to be capable of partially suppressing the phage DNA arrest phenotype caused by a dysfunction in genes 46 and/or 47 (also known as Mre11/Rad50 complex). Genetic mapping of das13 (one of the das mutations) has shown it to be in the region of the rnh gene encoding RNase H. Here we report that Das13 mutant of RNase H has substitutions of valine 43 and leucine 242 with isoleucines. To investigate the influence of these mutations on RNase H nuclease properties we have designed a novel in vitro assay that allows us to separate and quantify exo- or endonuclease activities of flap endonuclease. The nuclease assay in vitro showed that V43I substitution increased the ratio between exonuclease/endonuclease activities of RNase H whereas L242I substitution did not affect the nuclease activity of RNase H in vitro. However, both mutations were necessary for the full das effect in vivo. Molecular modelling of the nuclease structure suggests that V43I substitution may lead to disposition of H4 helix, responsible for the interaction with the first base pairs of 5'end of branched DNA. These structural changes may affect unwinding of the first base pairs of gapped or nicked DNA generating a short flap and therefore may stabilize the DNA-enzyme complex. L242I substitution did not affect the structure of RNase H and its role in providing das-effect remains unclear.

TAPO: A combined method for the identification of tandem repeats in protein structures

In recent years, there has been an emergence of new 3D structures of proteins containing tandem repeats (TRs), as a result of improved expression and crystallization strategies. Databases focused on structure classifications (PDB, SCOP, CATH) do not provide an easy solution for selection of these structures from PDB. Several approaches have been developed, but no best approach exists to identify the whole range of 3D TRs. Here we describe the TAndem PrOtein detector (TAPO) that uses periodicities of atomic coordinates and other types of structural representation, including strings generated by conformational alphabets, residue contact maps, and arrangements of vectors of secondary structure elements. The benchmarking shows the superior performance of TAPO over the existing programs. In accordance with our analysis of PDB using TAPO, 19% of proteins contain 3D TRs. This analysis allowed us to identify new families of 3D TRs, suggesting that TAPO can be used to regularly update the collection and classification of existing repetitive structures.

Restricted mobility of side chains on concave surfaces of solenoid proteins may impart heightened potential for intermolecular interactions

Significant progress has been made in the determination of the protein structures with their number today passing over a hundred thousand structures. The next challenge is the understanding and prediction of protein-protein and protein-ligand interactions. In this work we address this problem by analyzing curved solenoid proteins. Many of these proteins are considered as "hub molecules" for their high potential to interact with many different molecules and to be a scaffold for multisubunit protein machineries. Our analysis of these structures through molecular dynamics simulations reveals that the mobility of the side-chains on the concave surfaces of the solenoids is lower than on the convex ones. This result provides an explanation to the observed preferential binding of the ligands, including small and flexible ligands, to the concave surface of the curved solenoid proteins. The relationship between the landscapes and dynamic properties of the protein surfaces can be further generalized to the other types of protein structures and eventually used in the computer algorithms, allowing prediction of protein-ligand interactions by analysis of protein surfaces.