PoincaréMSA is a tool for protein family vizualisation starting from a multiple sequence alignment (either provided by the user or built by homologous search for a target sequence). It is available in the form of an interactive Google Colab notebooks and the underlying algorithm is described in Susmelj et al. [1]. PoincaréMSA takes as input a multiple sequence alignment (MSA) and builds its projection on a Poincaré disk using the method developed by Klimovskaia et al. in [2]. For the detailed tutorial and contacts please see: https://www.dsimb.inserm.fr/POINCARE_MSA
PoincareMSA builds an interactive projection of an input protein multiple sequence alignemnt (MSA) using a method developed by Susmelj et al. [1] based on Poincaré maps [2]. It reproduces both local proximities of protein sequences and hierarchy contained in the given data. Thus, sequences located closer to the center of projection correspond to the proteins sharing the most general functional properites and/or appearing at the earlier stages of evolution.
We provide three different Google Colab notebooks for interactive visualization of multiple sequence alignments:
- PoincareMSA_colab.ipynb takes as input a MSA in
.mfastaformat provided by a user. The user can also provide an annotation in.csvformat which will be used for coloring, as well as an UniProt IDs list used to automatically fetch taxonomy informations for coloring. - PoincareMSA_colab_examples.ipynb builds PoincareMSA projections from the example alignments available in
examplesdirectory. - PoincareMSA_colab_MMseqs2.ipynb performs a homologous sequence search for a target sequence and filtering of the resulting alignment with further projection by PoincaréMSA.
To get a local copy of the software run:
git clone git@github.com:DSIMB/PoincareMSA.git
cd PoincareMSA
The program is implemented in python3.7 using pytorch library for Poincaré disk construction and plotly for interactive visualisation of the resulting projections.
If you are working in Linux, you can use a conda environment to access all the necessary libraries:
conda env create -f env_poincare.yml
conda activate env_poincare
Otherwise here is a list of necessary dependencies to install:
pytorch 1.7.1
sklearn 0.24.1
numpy 1.19.2
pandas 1.2.3
scipy 1.6.0
seaborn 0.11.1
plotly 5.8.0
jax / jaxlib 0.3.25
The best way to try PoincaréMSA is by launching python notebooks with provided examples. To launch a particular example one needs to put the corresponding jupyter notebook to the repository root.
For example, to run PoincaréMSA on kinase dataset on should execute:
cp examples/kinases/PoincareMSA_kinases.ipynb ./
jupyter-notebook PoincareMSA_kinases.ipynb
The notebook can be then easily modified to work with any user-provided dataset.
Otherwise, the used can also launche the projection generation step by step as described below.
Data preparation
The user is invited to provide its MSA in the classical .mfasta format. Each sequence of the alignment is translated to a profile using position-specific scoring matrix (PSSM) according to the pseudo-count algortihm of Henikoff & Henikoff. The related scripts are located in scripts/prepare_data/ directory and driver scripts are provided for every example as create_projection.sh.
The resulting PSSM profiles representing each protein of MSA are stored in the directory fastas0.9, where 0.9 indicate the threshold percentage of gaps per position used to filter initial alignment. To build a Poincaré disk from this data, one needs to run a command from scripts/build_poincare_map/ directory:
python main.py --input_path path_to_files/fastas0.9/ --output_path output_dir --knn 5 --gamma 2 --batchsize 4 --epochs 1000
which will create an output .csv file with protein coordinates in the final projection and .png images reflecting the learning process. The .csv file can be further used to build interactive visualisation.
We provide several examples of PoincareMSA usage for different protien families in the examples folder. Each example comes with a bash script alloqing to reproduce the results starting from MSA and labels contained in data.
When using PoincaréMSA, please cite the following research:
[1] A. K. Susmelj, Y. Ren, Y. Vander Meersche, J.-C. Gelly, T. Galochkina. Poincaré maps for visualization of large protein families, Briefings in Bioinformatics, bbad103 (2023). https://doi.org/10.1093/bib/bbad103
The projection construction is adapted from the original code: https://github.com/facebookresearch/PoincareMaps developed for RNA sequence data visualization as described in the following paper:
[2] A. Klimovskaia, D. Lopez-Paz, L. Bottou et al. Poincaré maps for analyzing complex hierarchies in single-cell data. Nat Commun 11, 2966 (2020). https://doi.org/10.1038/s41467-020-16822-4
For scientific collaboration please contact Dr. Tatiana Galochkina at tatiana.galochkina@u-paris.fr and Dr. Jean-Christophe Gelly at jean-christophe.gelly@u-paris.fr.