The PRALINE database

Here we introduce PRALINE - Protein and Rna humAn singLe-nucleotIde-variaNts in condEnsates - a database to interrogate proteins and RNAs contained in human condensates such as stress-granules (SGs), processing bodies (PBs), droplets, amyloids and many more.. For proteins, PRALINE reports the liquid-liquid phase separation (LLPS) and liquid-solid phase separation (LSPT) propensities as well as the predicted and experimentally validated protein and RNA interactions. For RNAs, PRALINE reports the secondary structure content, as well as the predicted and experimentally validated protein and RNA interactions. All the information is also provided considering the presence of human disease-related variants (SNVs). Where available, experimental information on the physical state (LLPS and LSPT) is given with links to relevant literature. All the predictions reported in PRALINE are based on several tools developed in our group in >15 years of activity (e.g., 18514226, 21623348, 27899588, 27320918). The algorithms are widely employed by the scientific community.

RNA molecules

All RNAs sequences were retrieved from Ensembl version 105 (34791404). 1841 Stress Granule (SG) enriched RNAs were taken from (29129640) selecting the ones with log2(fold-change) >=1 and p-value <=0.01. 4852 Processing Body (PB) enriched RNAs are taken from (28965817), retrieving the ones with log2(fold-change) >=1 and q-value <= 0.01) for a total of 5614 unique genes upon removal of obsolete gene ids in the ensembl version 105. For each gene, the longest isoform was taken. The CROSS algorithm for secondary structure prediction of RNAs was launched on the 5614 RNA sequences with Globalscore parameter.

Protein molecules

997 UniprotKB IDs were retrieved from different sources. Stress Granule (SG) proteins were retrieved from Markmiller et al. (29373831), Gingras et al. (29395067) and Jain et al. (26777405) studies. From Jain work we retrieved 411 proteins that were either already known SG proteins or newly discovered through mass spec and IF. From Markmiller and colleagues, we retrieved 397 proteins either already known from previous experiments or found with APEX technique in hek293, NPC and IPSC cells, as well as proteins found to be stress / cell specific or independent. 60 SG proteins were retrieved from Gingras work. We retrieved proteins with Non-negative matrix factorization (NMF) values = 9 or in case of a different NMF value, we collected proteins found to co-localise with G3BP1.

Processing Body (PB) proteins were retrieved from Hustemberger et al. (28965817) and Gingras et al. studies (29395067). From the hustemberger study, we collected 125 proteins found to be significantly enriched in PBs with p-value < 0.025 as reported in the paper. From Gingras and colleagues' work, we retrieved 42 proteins either with NMF value =8 or that co-localize with DCP1A.

From Vendruscolo and Fuxreiter work (34391803) we retrieved 280 droplet forming proteins and 68 amyloid forming ones.

SNVs

DisGeNet curated variant-disease associations (31680165) and ClinVar variants (24234437) with a review status higher than one (“practice guideline”, “reviewed by expert panel” and “criteria provided, multiple submitters, no conflicts”) were used to retrieve disease related single nucleotide changes (SNVs). 13857 SNVs from DisGeNet and 48671 SNVs from ClinVar fell in the coding region of RNAs enriched in SG and PB and their relative position in the transcripts was retrieved from the Ensembl Variation 105 (34791404) in Human Short Variants dataset, excluding insertions and deletions. CROSS algorithm with Globalscore parameter was launched on RNA fragments of 51 nt with the SNV in the center to calculate the difference in secondary structure content between the reference and alternative RNA sequences. To avoid smaller-size fragments, SNVs falling at the beginning and at the end of the RNAs were removed. The reference and alternative CROSS secondary structure propensity profiles of the 51 nt RNA fragments was represented against each other in plots. The mean difference in secondary structure between reference and alternative sequence was computed on a window size of 11 nt upstream and downstream the SNV. Information on the expression quantitative trait loci (eQTL) and splicing quantitative trait loci (sQTL) of the SNVs was obtained from GTEx V8 (32913098).

RNA-RNA interactions

We retrieved human RNA-RNA interactions from RISE database focusing on the experimental interactions obtained with PARIS technique. For each RNA, binding site location were mapped to the longest transcript in Ensembl 105, using blastn algorithm. We retrieved a total of 25.232 RNA-RNA interactions with at least one interactor being an enriched RNA in SG or PB. In 934 of those interactions we found at least a SNV located inside a binding site.

Protein-RNA interactions

Human protein–RNA interactions were collected from eCLIP experiments (32728246) considering the ones found in either of the two cell lines, K562 and HepG230. The dataset contains interactions for 151 unique RBPs: 120 proteins in the K562 cell line and 103 in the HepG2 cell line. We used input-normalized eCLIP data (SM input normalization) with stringent cut-offs [−log10(p value) > 3 and −log2(fold_enrichment) >3]. The data was processed using bioconductor packages for R (https://www.bioconductor.org/) to map genomic locations to Ensembl transcript identifiers. Blastn algorithm was used to map the binding site location location in ensembl 105. We found 114.322 interactions involving SG-PB transcripts and in 5532 of those we found at least a SNV located in the binding site.

Protein-RNA predictions

The RNAct database covers around 6 millions interactions taken from RNAct. RNAct integrates experimentally available interactions (ENCODE, see above) with ab initio predictions (catRAPID), providing a full coverage of the known RNA-binding proteome.

Protein-Protein interactions

Human curated protein-protein interactions were collected from BioGRID database (version 4.4 - 33070389).

LLPS and LSPT Propensities

We computed the propensity to undergo liquid-like and solid-like condensation for the set of 997 proteins detailed previously and for their 6152 natural variants (involving 632 of them), retrieved from UniProtKB. We considered 5949 single point mutations and 203 deletions. To quantify the extent to which each sequence is prone to undergo LSPT and LLPS we used the Zyggregator and catGRANULE algorithms, which compute the liquid-solid and liquid-liquid propensities, respectively.

In the figure we show our LLPS and LSPT predictions considering a reference data set composed of amyloid prone regions APR (29040693) and droplet prone regions DPR (35025997 and 34391803)

For the proteins that do not have natural variants present in UniProtKB we only report the WT scores of Zyggregator and catGRANULE.

Figure 1 Cumulative Distribution Functions of LSPT and LLSP propensitives for experimentally-derived Amyloid-prone (blue) and Droplet-prine (oraange) protein regions. Statistics show that Amyloid and Droplet prone regions can be correctly distinguished with a p-value of 1.4 x 10^-4 and 7.2 x 10^-, respectively. We note that droplets can evolve towards the amyloid state, while the opposite does not occur.