MobiDB: 10 years of intrinsically disordered proteins.

The MobiDB database (URL: https://mobidb.org/) is a knowledge base of intrinsically disordered proteins. MobiDB aggregates disorder annotations derived from the literature and from experimental evidence along with predictions for all known protein sequences. MobiDB generates new knowledge and captures the functional significance of disordered regions by processing and combining complementary sources of information. Since its first release 10 years ago, the MobiDB database has evolved in order to improve the quality and coverage of protein disorder annotations and its accessibility. MobiDB has now reached its maturity in terms of data standardization and visualization. Here, we present a new release which focuses on the optimization of user experience and database content. The major advances compared to the previous version are the integration of AlphaFoldDB predictions and the re-implementation of the homology transfer pipeline, which expands manually curated annotations by two orders of magnitude. Finally, the entry page has been restyled in order to provide an overview of the available annotations along with two separate views that highlight structural disorder evidence and functions associated with different binding modes.

DisPhaseDB: An integrative database of diseases related variations in liquid-liquid phase separation proteins.

Motivation: Proteins involved in liquid-liquid phase separation (LLPS) and membraneless organelles (MLOs) are recognized to be decisive for many biological processes and also responsible for several diseases. The recent explosion of research in the area still lacks tools for the analysis and data integration among different repositories. Currently, there is not a comprehensive and dedicated database that collects all disease-related variations in combination with the protein location, biological role in the MLO, and all the metadata available for each protein and disease. Disease-related protein variants and additional features are dispersed and the user has to navigate many databases, with a different focus, formats, and often not user friendly. Results: We present DisPhaseDB, a database dedicated to disease-related variants of liquid-liquid phase separation proteins. It integrates 10 databases, contains 5,741 proteins, 1,660,059 variants, and 4,051 disease terms. It also offers intuitive navigation and an informative display. It constitutes a pivotal starting point for further analysis, encouraging the development of new computational tools.The database is freely available at http://disphasedb.leloir.org.ar.

DisProt in 2022: improved quality and accessibility of protein intrinsic disorder annotation.

The Database of Intrinsically Disordered Proteins (DisProt, URL: https://disprot.org) is the major repository of manually curated annotations of intrinsically disordered proteins and regions from the literature. We report here recent updates of DisProt version 9, including a restyled web interface, refactored Intrinsically Disordered Proteins Ontology (IDPO), improvements in the curation process and significant content growth of around 30%. Higher quality and consistency of annotations is provided by a newly implemented reviewing process and training of curators. The increased curation capacity is fostered by the integration of DisProt with APICURON, a dedicated resource for the proper attribution and recognition of biocuration efforts. Better interoperability is provided through the adoption of the Minimum Information About Disorder (MIADE) standard, an active collaboration with the Gene Ontology (GO) and Evidence and Conclusion Ontology (ECO) consortia and the support of the ELIXIR infrastructure.

ASpli: integrative analysis of splicing landscapes through RNA-Seq assays.

MOTIVATION: Genome-wide analysis of alternative splicing has been a very active field of research since the early days of next generation sequencing technologies. Since then, ever-growing data availability and the development of increasingly sophisticated analysis methods have uncovered the complexity of the general splicing repertoire. A large number of splicing analysis methodologies exist, each of them presenting its own strengths and weaknesses. For instance, methods exclusively relying on junction information do not take advantage of the large majority of reads produced in an RNA-seq assay, isoform reconstruction methods might not detect novel intron retention events, some solutions can only handle canonical splicing events, and many existing methods can only perform pairwise comparisons. RESULTS: In this contribution, we present ASpli, a computational suite implemented in R statistical language, that allows the identification of changes in both, annotated and novel alternative-splicing events and can deal with simple, multi-factor or paired experimental designs. Our integrative computational workflow, that considers the same GLM model applied to different sets of reads and junctions, allows computation of complementary splicing signals. Analyzing simulated and real data, we found that the consolidation of these signals resulted in a robust proxy of the occurrence of splicing alterations. While the analysis of junctions allowed us to uncover annotated as well as non-annotated events, read coverage signals notably increased recall capabilities at a very competitive performance when compared against other state-of-the-art splicing analysis algorithms. AVAILABILITY AND IMPLEMENTATION: ASpli is freely available from the Bioconductor project site https://doi.org/doi:10.18129/B9.bioc.ASpli. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PED in 2021: a major update of the protein ensemble database for intrinsically disordered proteins.

The Protein Ensemble Database (PED) (https://proteinensemble.org), which holds structural ensembles of intrinsically disordered proteins (IDPs), has been significantly updated and upgraded since its last release in 2016. The new version, PED 4.0, has been completely redesigned and reimplemented with cutting-edge technology and now holds about six times more data (162 versus 24 entries and 242 versus 60 structural ensembles) and a broader representation of state of the art ensemble generation methods than the previous version. The database has a completely renewed graphical interface with an interactive feature viewer for region-based annotations, and provides a series of descriptors of the qualitative and quantitative properties of the ensembles. High quality of the data is guaranteed by a new submission process, which combines both automatic and manual evaluation steps. A team of biocurators integrate structured metadata describing the ensemble generation methodology, experimental constraints and conditions. A new search engine allows the user to build advanced queries and search all entry fields including cross-references to IDP-related resources such as DisProt, MobiDB, BMRB and SASBDB. We expect that the renewed PED will be useful for researchers interested in the atomic-level understanding of IDP function, and promote the rational, structure-based design of IDP-targeting drugs.

Bovine Interferon Lambda Is a Potent Antiviral Against SARS-CoV-2 Infection in vitro.

Interferon lambda (IFN-λ) is an antiviral naturally produced in response to viral infections, with activity on cells of epithelial origin and located in the mucosal surfaces. This localized activity results in reduced toxicity compared to type I IFNs, whose receptors are ubiquitously expressed. IFN-λ has been effective in the therapy of respiratory viral infections, playing a crucial role in potentiating adaptive immune responses that initiate at mucosal surfaces. Human IFN-λ has polymorphisms that may cause differences in the interaction with the specific receptor in the human population. Interestingly, bovine IFN-λ3 has an in silico-predicted higher affinity for the human receptor than its human counterparts, with high identity with different human IFN-λ variants, making it a suitable antiviral therapeutic candidate for human health. Here, we demonstrate that a recombinant bovine IFN-λ (rbIFN-λ) produced in HEK-293 cells is effective in preventing SARS-CoV-2 infection of VERO cells, with an inhibitory concentration 50% (IC50) between 30 and 50 times lower than that of human type I IFN tested here (α2b and ß1a). We also demonstrated the absence of toxicity of rbIFN-λ in human PBMCs and the lack of proinflammatory activity on these cells. Altogether, our results show that rbIFN-λ is as an effective antiviral potentially suitable for COVID-19 therapy. Among other potential applications, rbIFN-λ could be useful to preclude virus dispersion to the lungs and/or to reduce transmission from infected people. Moreover, and due to the non-specific activity of this IFN, it can be potentially effective against other respiratory viruses that may be circulating together with SARS-CoV-2.

Chasing coevolutionary signals in intrinsically disordered proteins complexes.

Intrinsically disordered proteins/regions (IDPs/IDRs) are crucial components of the cell, they are highly abundant and participate ubiquitously in a wide range of biological functions, such as regulatory processes and cell signaling. Many of their important functions rely on protein interactions, by which they trigger or modulate different pathways. Sequence covariation, a powerful tool for protein contact prediction, has been applied successfully to predict protein structure and to identify protein-protein interactions mostly of globular proteins. IDPs/IDRs also mediate a plethora of protein-protein interactions, highlighting the importance of addressing sequence covariation-based inter-protein contact prediction of this class of proteins. Despite their importance, a systematic approach to analyze the covariation phenomena of intrinsically disordered proteins and their complexes is still missing. Here we carry out a comprehensive critical assessment of coevolution-based contact prediction in IDP/IDR complexes and detail the challenges and possible limitations that emerge from their analysis. We found that the coevolutionary signal is faint in most of the complexes of disordered proteins but positively correlates with the interface size and binding affinity between partners. In addition, we discuss the state-of-art methodology by biological interpretation of the results, formulate evaluation guidelines and suggest future directions of development to the field.

Prefoldins contribute to maintaining the levels of the spliceosome LSM2-8 complex through Hsp90 in Arabidopsis.

Although originally identified as the components of the complex aiding the cytosolic chaperonin CCT in the folding of actins and tubulins in the cytosol, prefoldins (PFDs) are emerging as novel regulators influencing gene expression in the nucleus. Work conducted mainly in yeast and animals showed that PFDs act as transcriptional regulators and participate in the nuclear proteostasis. To investigate new functions of PFDs, we performed a co-expression analysis in Arabidopsis thaliana. Results revealed co-expression between PFD and the Sm-like (LSM) genes, which encode the LSM2-8 spliceosome core complex, in this model organism. Here, we show that PFDs interact with and are required to maintain adequate levels of the LSM2-8 complex. Our data indicate that levels of the LSM8 protein, which defines and confers the functional specificity of the complex, are reduced in pfd mutants and in response to the Hsp90 inhibitor geldanamycin. We provide biochemical evidence showing that LSM8 is a client of Hsp90 and that PFD4 mediates the interaction between both proteins. Consistent with our results and with the role of the LSM2-8 complex in splicing through the stabilization of the U6 snRNA, pfd mutants showed reduced levels of this snRNA and altered pre-mRNA splicing patterns.

A Role for Pre-mRNA-PROCESSING PROTEIN 40C in the Control of Growth, Development, and Stress Tolerance in Arabidopsis thaliana.

Because of their sessile nature, plants have adopted varied strategies for growing and reproducing in an ever-changing environment. Control of mRNA levels and pre-mRNA alternative splicing are key regulatory layers that contribute to adjust and synchronize plant growth and development with environmental changes. Transcription and alternative splicing are thought to be tightly linked and coordinated, at least in part, through a network of transcriptional and splicing regulatory factors that interact with the carboxyl-terminal domain (CTD) of the largest subunit of RNA polymerase II. One of the proteins that has been shown to play such a role in yeast and mammals is pre-mRNA-PROCESSING PROTEIN 40 (PRP40, also known as CA150, or TCERG1). In plants, members of the PRP40 family have been identified and shown to interact with the CTD of RNA Pol II, but their biological functions remain unknown. Here, we studied the role of AtPRP40C, in Arabidopsis thaliana growth, development and stress tolerance, as well as its impact on the global regulation of gene expression programs. We found that the prp40c knockout mutants display a late-flowering phenotype under long day conditions, associated with minor alterations in red light signaling. An RNA-seq based transcriptome analysis revealed differentially expressed genes related to biotic stress responses and also differentially expressed as well as differentially spliced genes associated with abiotic stress responses. Indeed, the characterization of stress responses in prp40c mutants revealed an increased sensitivity to salt stress and an enhanced tolerance to Pseudomonas syringae pv. maculicola (Psm) infections. This constitutes the most thorough analysis of the transcriptome of a prp40 mutant in any organism, as well as the first characterization of the molecular and physiological roles of a member of the PRP40 protein family in plants. Our results suggest that PRP40C is an important factor linking the regulation of gene expression programs to the modulation of plant growth, development, and stress responses.

Shade delays flowering in Medicago sativa.

Shade-intolerant plants respond to the decrease in the red (R) to far-red (FR) light ratio (R:FR) occurring under shade by elongating stems and petioles and by re-positioning leaves, in a race to outcompete neighbors for the sunlight resource. In some annual species, the shade avoidance syndrome (SAS) is accompanied by the early induction of flowering. Anticipated flowering is viewed as a strategy to set seeds before the resources become severely limiting. Little is known about the molecular mechanisms of SAS in perennial forage crops like alfalfa (Medicago sativa). To study SAS in alfalfa, we exposed alfalfa plants to simulated shade by supplementing with FR light. Low R:FR light produced a classical SAS, with increased internode and petiole lengths, but unexpectedly also with delayed flowering. To understand the molecular mechanisms involved in uncoupling SAS from early flowering, we used a transcriptomic approach. The SAS is likely to be mediated by increased expression of msPIF3 and msHB2 in low R:FR light. Constitutive expression of these genes in Arabidopsis led to SAS, including early flowering, strongly suggesting that their roles are conserved. Delayed flowering was likely to be mediated by the downregulation of msSPL3, which promotes flowering in both Arabidopsis and alfalfa. Shade-delayed flowering in alfalfa may be important to extend the vegetative phase under suboptimal light conditions, and thus assure the accumulation of reserves necessary to resume growth after the next season.

MISTIC2: comprehensive server to study coevolution in protein families.

Correlated mutations between residue pairs in evolutionarily related proteins arise from constraints needed to maintain a functional and stable protein. Identifying these inter-related positions narrows down the search for structurally or functionally important sites. MISTIC is a server designed to assist users to calculate covariation in protein families and provide them with an interactive tool to visualize the results. Here, we present MISTIC2, an update to the previous server, that allows to calculate four covariation methods (MIp, mfDCA, plmDCA and gaussianDCA). The results visualization framework has been reworked for improved performance, compatibility and user experience. It includes a circos representation of the information contained in the alignment, an interactive covariation network, a 3D structure viewer and a sequence logo. Others components provide additional information such as residue annotations, a roc curve for assessing contact prediction, data tables and different ways of filtering the data and exporting figures. Comparison of different methods is easily done and scores combination is also possible. A newly implemented web service allows users to access MISTIC2 programmatically using an API to calculate covariation and retrieve results. MISTIC2 is available at: https://mistic2.leloir.org.ar.

Exploring Genomic, Geographic and Virulence Interactions among Epidemic and Non-Epidemic St. Louis Encephalitis Virus (Flavivirus) Strains.

St. Louis encephalitis virus (SLEV) is a re-emerging arbovirus in South America. In 2005, an encephalitis outbreak caused by SLEV was reported in Argentina. The reason for the outbreak remains unknown, but may have been related to virological factors, changes in vectors populations, avian amplifying hosts, and/or environmental conditions. The main goal of this study was to characterize the complete genome of epidemic and non-epidemic SLEV strains from Argentina. Seventeen amino acid changes were detected; ten were non-conservative and located in proteins E, NS1, NS3 and NS5. Phylogenetic analysis showed two major clades based on geography: the North America and northern Central America (NAnCA) clade and the South America and southern Central America (SAsCA) clade. Interestingly, the presence of SAsCA genotype V SLEV strains in the NAnCA clade was reported in California, Florida and Texas, overlapping with known bird migration flyways. This work represents the first step in understanding the molecular mechanisms underlying virulence and biological variation among SLEV strains.

I-COMS: Interprotein-COrrelated Mutations Server.

Interprotein contact prediction using multiple sequence alignments (MSAs) is a useful approach to help detect protein-protein interfaces. Different computational methods have been developed in recent years as an approximation to solve this problem. However, as there are discrepancies in the results provided by them, there is still no consensus on which is the best performing methodology. To address this problem, I-COMS (interprotein COrrelated Mutations Server) is presented. I-COMS allows to estimate covariation between residues of different proteins by four different covariation methods. It provides a graphical and interactive output that helps compare results obtained using different methods. I-COMS automatically builds the required MSA for the calculation and produces a rich visualization of either intraprotein and/or interprotein covariating positions in a circos representation. Furthermore, comparison between any two methods is available as well as the overlap between any or all four methodologies. In addition, as a complementary source of information, a matrix visualization of the corresponding scores is made available and the density plot distribution of the inter, intra and inter+intra scores are calculated. Finally, all the results can be downloaded (including MSAs, scores and graphics) for comparison and visualization and/or for further analysis.

Family-specific degenerate primer design: a tool to design consensus degenerated oligonucleotides.

Designing degenerate PCR primers for templates of unknown nucleotide sequence may be a very difficult task. In this paper, we present a new method to design degenerate primers, implemented in family-specific degenerate primer design (FAS-DPD) computer software, for which the starting point is a multiple alignment of related amino acids or nucleotide sequences. To assess their efficiency, four different genome collections were used, covering a wide range of genomic lengths: Arenavirus (10 × 10(4) nucleotides), Baculovirus (0.9 × 10(5) to 1.8 × 10(5) bp), Lactobacillus sp. (1 × 10(6) to 2 × 10(6) bp), and Pseudomonas sp. (4 × 10(6) to 7 × 10(6) bp). In each case, FAS-DPD designed primers were tested computationally to measure specificity. Designed primers for Arenavirus and Baculovirus were tested experimentally. The method presented here is useful for designing degenerate primers on collections of related protein sequences, allowing detection of new family members.

Antigen vehiculization particles based on the Z protein of Junin virus.

BACKGROUND: Arenavirus matrix protein Z plays an important role in virus budding and is able to generate enveloped virus-like-particles (VLPs) in absence of any other viral proteins. In these VLPs, Z protein is associated to the plasma membrane inner surface by its myristoyl residue. Budding induction and vesicle formation properties can be exploited to generate enveloped VLPs platform. These structures can be designed to carry specific antigen in the inner side or on the surface of VLPs.Vaccines based on VLPs are a highly effective type of subunit vaccines that mimic the overall structure of virus particles in absence of viral nucleic acid, being noninfectious.In this work we assayed the capacity of Junin Z protein to produce VLPs carrying the green fluorescent protein (eGFP), as a model antigen. RESULTS: In this report the Junin Z protein ability to produce VLPs from 293T cells and its capacity to deliver a specific antigen (eGFP) fused to Z was evaluated. Confocal microscopy showed a particular membrane bending in cells expressing Z and a spot welded distribution in the cytoplasm. VLPs were detected by TEM (transmission electron microscopy) and were purified from cell supernatant. The proteinase protection assay demonstrated the VLPs integrity and the absence of degradation of the fused antigen, thus indicating its internal localization. Finally, immunization of mice with purified VLPs produced high titres of anti-eGFP antibodies compared to the controls. CONCLUSIONS: It was proved that VLPs can be generated from cells transfected with a fusion Junin virus Z-eGFP protein in absence of any other viral protein, and the capacity of Z protein to support fusions at the C-terminal, without impairing its budding activity, allowing vehiculization of specific antigens into VLPs.

The ac53, ac78, ac101, and ac103 genes are newly discovered core genes in the family Baculoviridae.

The family Baculoviridae is a large group of insect viruses containing circular double-stranded DNA genomes of 80 to 180 kbp, which have broad biotechnological applications. A key feature to understand and manipulate them is the recognition of orthology. However, the differences in gene contents and evolutionary distances among the known members of this family make it difficult to assign sequence orthology. In this study, the genome sequences of 58 baculoviruses were analyzed, with the aim to detect previously undescribed core genes because of their remote homology. A routine based on Multi PSI-Blast/tBlastN and Multi HaMStR allowed us to detect 31 of 33 accepted core genes and 4 orthologous sequences in the Baculoviridae which were not described previously. Our results show that the ac53, ac78, ac101 (p40), and ac103 (p48) genes have orthologs in all genomes and should be considered core genes. Accordingly, there are 37 orthologous genes in the family Baculoviridae.

Viral diversity of Junín virus field strains.

The Argentine Hemorrhagic Fever, an endemic disease present in a much of Argentina, is caused by the Junín virus (JUNV). Currently, there are sequences available from several strains of this virus, like those belonging to the vaccine lineage (XJ13, XJ#44 and Candid#1), as well as MC2 (rodent isolate) and IV4454 (human isolate). In this article, we report sequence information on two fragments of genomic segment S of viral isolates from the endemic area. A Nested-RT-PCR was used to amplify discrete genomic regions of 13 isolates of rodent and human origin. The bioinformatics studies revealed a great homogeneity of sequences among the JUNV isolates. The phylogenetic classification showed greater evolutionary distance between the old world arenaviruses (Lassa and LCM virus) than between the new world arenaviruses (JUNV and Machupo virus).

Molecular analysis of the virulence attenuation process in Junín virus vaccine genealogy.

The Junín virus strain Candid#1 was developed as a live attenuated vaccine for Argentine hemorrhagic fever. In this article, we report sequence information of the L and S RNAs of Junín virus Candid#1 and XJ#44 strains, and show the comparisons with the XJ13 wild-type strain and with other Junín virus strains, like Romero, IV4454 and MC2 strains, and other closely and distantly related arenaviruses. Comparisons of the nucleotide and amino acid sequences of all genes of three strains from the same vaccine genealogy, revealed different point mutations that could be associated with the attenuated phenotype. A 91% of the mutations found are consistent with a hypothesis of progressive attenuation of virulence from XJ13 to XJ#44 and to Candid#1; 39% of mutations were observed in XJ#44 and conserved in Candid#1, while another 52% of the mutations appeared only in Candid#1 strain. The remaining 9% corresponded to reverse mutations in the L gene. In summary, the present work shows a set of mutations that could be related to the virulence attenuation phenomenon. This information will serve as a starting point to study this biological phenomenon, provided that a reverse genetics system for Junín virus is developed to allow the generation of infectious virions with specific mutations.

Genomic features of attenuated Junín virus vaccine strain candidate.

Junin virus strain Candid #1 was developed as a live attenuated vaccine for Argentine haemorrhagic fever. In this paper, we report the nucleotide sequences of L RNA of Candid #1 and examine the relationship to its more virulent ancestors Junin virus XJ#44 and XJ 13 (prototype) and other closely and distantly related arenaviruses. Comparisons of the nucleotide and amino acid sequences of L and Z genes of Candid #1 and its progenitor strains revealed twelve point mutations in the L polypeptide that are unique to the vaccine strain. These changes could be provisionally associated with the attenuated phenotype. In contrast, Z ORF was completely conserved among all strains.