Protein aggregation, structural disorder and RNA-binding ability: a new approach for physico-chemical and gene ontology classification of multiple datasets

Computational optimization of angiotensin-converting enzyme 2 for SARS-CoV-2 Spike molecular recognition

Since the beginning of the Covid19 pandemic, many efforts have been devoted to identifying approaches to neutralize SARS-CoV-2 replication within the host cell. A promising strategy to block the infection consists of using a mutant of the human receptor angiotensin-converting enzyme 2 (ACE2) as a decoy to compete with endogenous ACE2 for the binding to the SARS-CoV-2 Spike protein, which decreases the ability of the virus to enter the host cell. Here, using a computational framework based on the 2D Zernike formalism we investigate details of the molecular binding and evaluate the changes in ACE2-Spike binding compatibility upon mutations occurring in the ACE2 side of the molecular interface. We demonstrate the efficacy of our method by comparing our results with experimental binding affinities changes upon ACE2 mutations, separating ones that increase or decrease binding affinity with an Area Under the ROC curve ranging from 0.66 to 0.93, depending on the magnitude of the effects analyzed. Importantly, the iteration of our approach leads to the identification of a set of ACE2 mutants characterized by an increased shape complementarity with Spike. We investigated the physico-chemical properties of these ACE2 mutants and propose them as bona fide candidates for Spike recognition.

Structural analysis of SARS-CoV-2 genome and predictions of the human interactome

Specific elements of viral genomes regulate interactions within host cells. Here, we calculated the secondary structure content of >2000 coronaviruses and computed >100 000 human protein interactions with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genomic regions display different degrees of conservation. SARS-CoV-2 domain encompassing nucleotides 22 500-23 000 is conserved both at the sequence and structural level. The regions upstream and downstream, however, vary significantly. This part of the viral sequence codes for the Spike S protein that interacts with the human receptor angiotensin-converting enzyme 2 (ACE2). Thus, variability of Spike S is connected to different levels of viral entry in human cells within the population. Our predictions indicate that the 5' end of SARS-CoV-2 is highly structured and interacts with several human proteins. The binding proteins are involved in viral RNA processing, include double-stranded RNA specific editases and ATP-dependent RNA-helicases and have strong propensity to form stress granules and phase-separated assemblies. We propose that these proteins, also implicated in viral infections such as HIV, are selectively recruited by SARS-CoV-2 genome to alter transcriptional and post-transcriptional regulation of host cells and to promote viral replication.

Structural analysis of SARS-CoV-2 genome and predictions of the human interactome

Specific elements of viral genomes regulate interactions within host cells. Here, we calculated the secondary structure content of >2000 coronaviruses and computed >100 000 human protein interactions with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genomic regions display different degrees of conservation. SARS-CoV-2 domain encompassing nucleotides 22 500-23 000 is conserved both at the sequence and structural level. The regions upstream and downstream, however, vary significantly. This part of the viral sequence codes for the Spike S protein that interacts with the human receptor angiotensin-converting enzyme 2 (ACE2). Thus, variability of Spike S is connected to different levels of viral entry in human cells within the population. Our predictions indicate that the 5' end of SARS-CoV-2 is highly structured and interacts with several human proteins. The binding proteins are involved in viral RNA processing, include double-stranded RNA specific editases and ATP-dependent RNA-helicases and have strong propensity to form stress granules and phase-separated assemblies. We propose that these proteins, also implicated in viral infections such as HIV, are selectively recruited by SARS-CoV-2 genome to alter transcriptional and post-transcriptional regulation of host cells and to promote viral replication.

RNA-binding and prion domains: the Yin and Yang of phase separation

Proteins and RNAs assemble in membrane-less organelles that organize intracellular spaces and regulate biochemical reactions. The ability of proteins and RNAs to form condensates is encoded in their sequences, yet it is unknown which domains drive the phase separation (PS) process and what are their specific roles. Here, we systematically investigated the human and yeast proteomes to find regions promoting condensation. Using advanced computational methods to predict the PS propensity of proteins, we designed a set of experiments to investigate the contributions of Prion-Like Domains (PrLDs) and RNA-binding domains (RBDs). We found that one PrLD is sufficient to drive PS, whereas multiple RBDs are needed to modulate the dynamics of the assemblies. In the case of stress granule protein Pub1 we show that the PrLD promotes sequestration of protein partners and the RBD confers liquid-like behaviour to the condensate. Our work sheds light on the fine interplay between RBDs and PrLD to regulate formation of membrane-less organelles, opening up the avenue for their manipulation.

In silico, in vitro, and in vivo Approaches to Identify Molecular Players in Fragile X Tremor and Ataxia Syndrome

Fragile X-associated tremor/ataxia syndrome (FXTAS) is a late-onset neurodegenerative monogenetic disorder affecting carriers of premutation (PM) forms of the FMR1 gene, resulting in a progressive development of tremors, ataxia, and neuropsychological problems. This highly disabling disease is quite common in the general population with an estimation of about 20 million PM carriers worldwide. The chances of developing FXTAS increase dramatically with age, with about 45% of male carriers over the age of 50 being affected. Both the gene and pathogenic trigger, a mutant expansion of CGG RNA, causing FXTAS are known. This makes it an interesting disease to develop targeted therapeutic interventions for. Yet, no such interventions are available at this moment. Here we discuss in silico, in vitro, and in vivo approaches and how they have been used to identify the molecular determinants of FXTAS pathology. These approaches have yielded substantial information about FXTAS pathology and, consequently, many markers have emerged to play a key role in understanding the disease mechanism. Integration of the different approaches is expected to provide crucial information about the value of these markers as either therapeutic target or biomarker, essential to monitor therapeutic interventions in the future.

An Integrative Study of Protein-RNA Condensates Identifies Scaffolding RNAs and Reveals Players in Fragile X-Associated Tremor/Ataxia Syndrome

Recent evidence indicates that specific RNAs promote the formation of ribonucleoprotein condensates by acting as scaffolds for RNA-binding proteins (RBPs). We systematically investigated RNA-RBP interaction networks to understand ribonucleoprotein assembly. We found that highly contacted RNAs are structured, have long UTRs, and contain nucleotide repeat expansions. Among the RNAs with such properties, we identified the FMR1 3' UTR that harbors CGG expansions implicated in fragile X-associated tremor/ataxia syndrome (FXTAS). We studied FMR1 binding partners in silico and in vitro and prioritized the splicing regulator TRA2A for further characterization. In a FXTAS cellular model, we validated the TRA2A-FMR1 interaction and investigated implications of its sequestration at both transcriptomic and post-transcriptomic levels. We found that TRA2A co-aggregates with FMR1 in a FXTAS mouse model and in post-mortem human samples. Our integrative study identifies key components of ribonucleoprotein aggregates, providing links to neurodegenerative disease and allowing the discovery of therapeutic targets.

CROSSalive: a web server for predicting the in vivo structure of RNA molecules

MOTIVATION

RNA structure is difficult to predict in vivo due to interactions with enzymes and other molecules. Here we introduce CROSSalive, an algorithm to predict the single- and double-stranded regions of RNAs in vivo using predictions of protein interactions.

RESULTS

Trained on icSHAPE data in presence (m6a+) and absence of N6 methyladenosine modification (m6a-), CROSSalive achieves cross-validation accuracies between 0.70 and 0.88 in identifying high-confidence single- and double-stranded regions. The algorithm was applied to the long non-coding RNA Xist (17 900 nt, not present in the training) and shows an Area under the ROC curve of 0.83 in predicting structured regions.

AVAILABILITY AND IMPLEMENTATION

CROSSalive webserver is freely accessible at http://service.tartaglialab.com/new_submission/crossalive.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

RNA structure drives interaction with proteins

The combination of high-throughput sequencing and in vivo crosslinking approaches leads to the progressive uncovering of the complex interdependence between cellular transcriptome and proteome. Yet, the molecular determinants governing interactions in protein-RNA networks are not well understood. Here we investigated the relationship between the structure of an RNA and its ability to interact with proteins. Analysing in silico, in vitro and in vivo experiments, we find that the amount of double-stranded regions in an RNA correlates with the number of protein contacts. This relationship -which we call structure-driven protein interactivity- allows classification of RNA types, plays a role in gene regulation and could have implications for the formation of phase-separated ribonucleoprotein assemblies. We validate our hypothesis by showing that a highly structured RNA can rearrange the composition of a protein aggregate. We report that the tendency of proteins to phase-separate is reduced by interactions with specific RNAs.

Advances in the characterization of RNA-binding proteins

From transcription, to transport, storage, and translation, RNA depends on association with different RNA-binding proteins (RBPs). Methods based on next-generation sequencing and protein mass-spectrometry have started to unveil genome-wide interactions of RBPs but many aspects still remain out of sight. How many of the binding sites identified in high-throughput screenings are functional? A number of computational methods have been developed to analyze experimental data and to obtain insights into the specificity of protein-RNA interactions. How can theoretical models be exploited to identify RBPs? In addition to oligomeric complexes, protein and RNA molecules can associate into granular assemblies whose physical properties are still poorly understood. What protein features promote granule formation and what effects do these assemblies have on cell function? Here, we describe the newest in silico, in vitro, and in vivo advances in the field of protein-RNA interactions. We also present the challenges that experimental and computational approaches will have to face in future studies. WIREs RNA 2016, 7:793-810. doi: 10.1002/wrna.1378 For further resources related to this article, please visit the WIREs website.

Non-random distribution of homo-repeats: links with biological functions and human diseases

The biological function of multiple repetitions of single amino acids, or homo-repeats, is largely unknown, but their occurrence in proteins has been associated with more than 20 hereditary diseases. Analysing 122 bacterial and eukaryotic genomes, we observed that the number of proteins containing homo-repeats is significantly larger than expected from theoretical estimates. Analysis of statistical significance indicates that the minimal size of homo-repeats varies with amino acid type and proteome. In an attempt to characterize proteins harbouring long homo-repeats, we found that those containing polar or small amino acids S, P, H, E, D, K, Q and N are enriched in structural disorder as well as protein- and RNA-interactions. We observed that E, S, Q, G, L, P, D, A and H homo-repeats are strongly linked with occurrence in human diseases. Moreover, S, E, P, A, Q, D and T homo-repeats are significantly enriched in neuronal proteins associated with autism and other disorders. We release a webserver for further exploration of homo-repeats occurrence in human pathology at http://bioinfo.protres.ru/hradis/.