Awesome

ProtTrans is providing state of the art pre-trained models for proteins. ProtTrans was trained on thousands of GPUs from Summit and hundreds of Google TPUs using various Transformer models.

Have a look at our paper ProtTrans: cracking the language of life’s code through self-supervised deep learning and high performance computing for more information about our work.

This repository will be updated regulary with new pre-trained models for proteins as part of supporting bioinformatics community in general, and Covid-19 research specifically through our Accelerate SARS-CoV-2 research with transfer learning using pre-trained language modeling models project.

Table of Contents

• ⌛️  News
• 🚀  Installation
• 🚀  Quick Start
• ⌛️  Models Availability
• ⌛️  Dataset Availability
• 🚀  Usage
  • 🧬  Feature Extraction (FE)
  • 🚀  Logits extraction
  • 💥  Fine Tuning (FT)
  • 🧠  Prediction
  • ⚗️  Protein Sequences Generation
  • 🧐  Visualization
  • 📈  Benchmark
• 📊  Original downstream Predictions
• 📊  Followup use-cases
• 📊  Comparisons to other tools
• ❤️  Community and Contributions
• 📫  Have a question?
• 🤝  Found a bug?
• ✅  Requirements
• 🤵  Team
• 💰  Sponsors
• 📘  License
• ✏️  Citation

<a name="news"></a>

⌛️  News

• 2023/07/14: FineTuning with LoRA provides a notebooks on how to fine-tune ProtT5 on both, per-residue and per-protein tasks, using Low-Rank Adaptation (LoRA) for efficient finetuning (thanks @0syrys !).
• 2022/11/18: Availability: LambdaPP offers a simple web-service to access ProtT5-based predictions and UniProt now offers to download pre-computed ProtT5 embeddings for a subset of selected organisms.

<a name="install"></a>

🚀  Installation

All our models are available via huggingface/transformers:

pip install torch
pip install transformers
pip install sentencepiece

For more details, please follow the instructions for transformers installations.

A recently introduced change in the T5-tokenizer results in UnboundLocalError: cannot access local variable 'sentencepiece_model_pb2 and can either be fixed by installing this PR or by manually installing:

pip install protobuf

If you are using a transformer version after this PR, you will see this warning. Explicitly setting legacy=True will result in expected behavor and will avoid the warning. You can also safely ignore the warning as legacy=True is the default.

<a name="quick"></a>

🚀  Quick Start

Example for how to derive embeddings from our best-performing protein language model, ProtT5-XL-U50 (aka ProtT5); also available as colab:

from transformers import T5Tokenizer, T5EncoderModel
import torch
import re

device = torch.device('cuda:0' if torch.cuda.is_available() else 'cpu')

# Load the tokenizer
tokenizer = T5Tokenizer.from_pretrained('Rostlab/prot_t5_xl_half_uniref50-enc', do_lower_case=False)

# Load the model
model = T5EncoderModel.from_pretrained("Rostlab/prot_t5_xl_half_uniref50-enc").to(device)

# only GPUs support half-precision currently; if you want to run on CPU use full-precision (not recommended, much slower)
if device == torch.device("cpu"):
    model.to(torch.float32)

# prepare your protein sequences as a list
sequence_examples = ["PRTEINO", "SEQWENCE"]

# replace all rare/ambiguous amino acids by X and introduce white-space between all amino acids
sequence_examples = [" ".join(list(re.sub(r"[UZOB]", "X", sequence))) for sequence in sequence_examples]

# tokenize sequences and pad up to the longest sequence in the batch
ids = tokenizer(sequence_examples, add_special_tokens=True, padding="longest")

input_ids = torch.tensor(ids['input_ids']).to(device)
attention_mask = torch.tensor(ids['attention_mask']).to(device)

# generate embeddings
with torch.no_grad():
    embedding_repr = model(input_ids=input_ids, attention_mask=attention_mask)

# extract residue embeddings for the first ([0,:]) sequence in the batch and remove padded & special tokens ([0,:7]) 
emb_0 = embedding_repr.last_hidden_state[0,:7] # shape (7 x 1024)
# same for the second ([1,:]) sequence but taking into account different sequence lengths ([1,:8])
emb_1 = embedding_repr.last_hidden_state[1,:8] # shape (8 x 1024)

# if you want to derive a single representation (per-protein embedding) for the whole protein
emb_0_per_protein = emb_0.mean(dim=0) # shape (1024)

We also have a script which simplifies deriving per-residue and per-protein embeddings from ProtT5 for a given FASTA file:

python prott5_embedder.py --input sequences/some.fasta --output embeddings/residue_embeddings.h5
python prott5_embedder.py --input sequences/some.fasta --output embeddings/protein_embeddings.h5 --per_protein 1

<a name="models"></a>

⌛️  Models Availability

<a name="datasets"></a>

⌛️  Datasets Availability

<a name="usage"></a>

🚀  Usage

How to use ProtTrans:

<a name="feature-extraction"></a>

• <b>🧬  Feature Extraction (FE):</b><br/> Please check: Embedding Section. Colab example for feature extraction via ProtT5-XL-U50

<a name="logits-extraction"></a>

• <b>🚀  Logits Extraction:</b><br/> For ProtT5-logits extraction, please check: VESPA logits script.

<a name="fine-tuning"></a>

• <b>💥  Fine Tuning (FT):</b><br/> Please check: Fine Tuning Section. More information coming soon.

<a name="prediction"></a>

• <b>🧠  Prediction:</b><br/> Please check: Prediction Section. Colab example for secondary structure prediction via ProtT5-XL-U50 and Colab example for subcellular localization prediction as well as differentiation between membrane-bound and water-soluble proteins via ProtT5-XL-U50.

<a name="protein-generation"></a>

• <b>⚗️  Protein Sequences Generation:</b><br/> Please check: Generate Section. More information coming soon.

<a name="visualization"></a>

• <b>🧐  Visualization:</b><br/> Please check: Visualization Section. More information coming soon.

<a name="benchmark"></a>

• <b>📈  Benchmark:</b><br/> Please check: Benchmark Section. More information coming soon.

<a name="results"></a>

📊  Original downstream Predictions

<a name="q3"></a>

• <b>🧬  Secondary Structure Prediction (Q3):</b><br/>

🆕 Predict your sequence live on predictprotein.org.

<a name="q8"></a>

• <b>🧬  Secondary Structure Prediction (Q8):</b><br/>

🆕 Predict your sequence live on predictprotein.org.

<a name="q2"></a>

• <b>🧬  Membrane-bound vs Water-soluble (Q2):</b><br/>

<a name="q10"></a>

• <b>🧬  Subcellular Localization (Q10):</b><br/>

<a name="inaction"></a>

📊  Use-cases

<a name="comparison"></a>

📊  Comparison to other protein language models (pLMs)

While developing the use-cases, we compared ProtTrans models to other protein language models, for instance the ESM models. To focus on the effect of changing input representaitons, the following comparisons use the same architectures on top on different embedding inputs.

Important note on ProtT5-XL-UniRef50 (dubbed ProtT5-XL-U50): all performances were measured using only embeddings extracted from the encoder-side of the underlying T5 model as described here. Also, experiments were ran in half-precision mode (model.half()), to speed-up embedding generation. No performance degradation could be observed in any of the experiments when running in half-precision.

<a name="community"></a>

❤️  Community and Contributions

The ProtTrans project is a open source project supported by various partner companies and research institutions. We are committed to share all our pre-trained models and knowledge. We are more than happy if you could help us on sharing new ptrained models, fixing bugs, proposing new feature, improving our documentation, spreading the word, or support our project.

<a name="question"></a>

📫  Have a question?

We are happy to hear your question in our issues page ProtTrans! Obviously if you have a private question or want to cooperate with us, you can always reach out to us directly via our RostLab email

<a name="bug"></a>

🤝  Found a bug?

Feel free to file a new issue with a respective title and description on the the ProtTrans repository. If you already found a solution to your problem, we would love to review your pull request!.

<a name="requirements"></a>

✅  Requirements

For protein feature extraction or fine-tuninng our pre-trained models, Pytorch and Transformers library from huggingface is needed. For model visualization, you need to install BertViz library.

<a name="team"></a>

🤵  Team

• <b>Technical University of Munich:</b><br/>

• <b>Med AI Technology:</b><br/>

• <b>Google:</b><br/>

• <b>Nvidia:</b><br/>

• <b>Seoul National University:</b><br/>

• <b>ORNL:</b><br/>

<a name="sponsors"></a>

💰  Sponsors

Nvidia	Google	Google	ORNL	Software Campus

<a name="license"></a>

📘  License

The ProtTrans pretrained models are released under the under terms of the Academic Free License v3.0 License.

<a name="citation"></a>

✏️  Citation

If you use this code or our pretrained models for your publication, please cite the original paper:

@ARTICLE
{9477085,
author={Elnaggar, Ahmed and Heinzinger, Michael and Dallago, Christian and Rehawi, Ghalia and Yu, Wang and Jones, Llion and Gibbs, Tom and Feher, Tamas and Angerer, Christoph and Steinegger, Martin and Bhowmik, Debsindhu and Rost, Burkhard},
journal={IEEE Transactions on Pattern Analysis and Machine Intelligence},
title={ProtTrans: Towards Cracking the Language of Lifes Code Through Self-Supervised Deep Learning and High Performance Computing},
year={2021},
volume={},
number={},
pages={1-1},
doi={10.1109/TPAMI.2021.3095381}}