File size: 4,124 Bytes
35905be
 
5e4e8b3
 
35905be
954c9f3
 
 
 
 
 
5e4e8b3
 
771d93b
5e4e8b3
954c9f3
f77377b
 
 
9f7c1ad
f77377b
 
 
 
 
 
 
 
 
 
 
 
954c9f3
 
 
fe0550a
f77377b
9f7c1ad
 
 
 
f77377b
 
 
9f7c1ad
f77377b
9f7c1ad
 
 
954c9f3
 
 
 
 
 
 
 
 
 
13b5d09
954c9f3
 
13b5d09
 
 
954c9f3
13b5d09
954c9f3
 
13b5d09
 
 
 
 
 
 
954c9f3
 
 
 
 
 
 
 
 
e7cc24d
c91eb3d
954c9f3
 
 
 
13b5d09
954c9f3
 
 
 
 
 
 
 
 
 
 
 
 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
---
license: cc-by-nc-sa-4.0
pipeline_tag: feature-extraction
library_name: transformers
---

# Model Details:
Ankh3 is a protein language model that is jointly optimized on two objectives:
* Masked language modeling with multiple masking probabilities
* Protein sequence completion.

This is the model of the paper [Ankh3: Multi-Task Pretraining with Sequence Denoising and Completion Enhances Protein Representations](https://huggingface.co/papers/2505.20052).

Code: https://github.com/agemagician/Ankh

1. Masked Language Modeling:
  - The idea of this task is to intentionally 'corrupt' an input protein sequence by
    masking a certain percentage (X%) of its individual tokens (amino acids),
    and then train the model to reconstruct the original sequence.
  
  - Example on a protein sequence before and after corruption:

    Original protein sequence: MKAYVLINSRGP

    This sequence will be masked/corrupted using sentinel tokens as shown below:
    Sequence after corruption: M <extra_id_0> A Y <extra_id_1> L I <extra_id_2> S R G <extra_id_3>


    The decoder learns to correspond each sentinel token to the actual amino acid that was masked.
    In this example: <extra_id_0> K means that <extra_id_0> corresponds to the "K" amino acid and so on.

    Decoder output: <extra_id_0> K <extra_id_1> V <extra_id_2> N <extra_id_3> P



2. Protein Sequence Completion:
- The idea of this task is to cut the input sequence into
  two segments, where the first segment is fed to the encoder
  and the decoder is tasked to auto-regressively generate the
  second segment conditioned on the first segment representation
  outputted from the encoder.

- Example on protein sequence completion:

  Original sequence: MKAYVLINSRGP

  We will pass "MKAYVL" of it to the encoder, and the decoder is trained
  that given the representation of the first part provided by the encoder,
  it should output the second part which is: "INSRGP"



# How to use:

## For Embedding Extraction:
```python
from transformers import T5ForConditionalGeneration, T5Tokenizer, T5EncoderModel
import torch

# Random sequence from uniprot, most likely Ankh3 saw it during pre-training.
sequence = "MDTAYPREDTRAPTPSKAGAHTALTLGAPHPPPRDHLIWSVFSTLYLNLCCLGFLALAYSIKARDQKVVGDLEAARRFGSKAKCYNILAAMWTLVPPLLLLGLVVTGALHLARLAKDSAAFFSTKFDDADYD"

ckpt = "ElnaggarLab/ankh3-large"

# Make sure that you must use `T5Tokenizer` not `AutoTokenizer`.
tokenizer = T5Tokenizer.from_pretrained(ckpt)

# To use the encoder representation using the NLU prefix:
encoder_model = T5EncoderModel.from_pretrained(ckpt).eval()


# For extracting embeddings, consider trying the '[S2S]' prefix.
# Since this prefix was specifically used to denote sequence completion
# during the model's pre-training, its use can sometimes
# lead to improved embedding quality.

nlu_sequence = "[NLU]" + sequence
encoded_nlu_sequence = tokenizer(nlu_sequence, add_special_tokens=True, return_tensors="pt", is_split_into_words=False)

with torch.no_grad():
  embedding = encoder_model(**encoded_nlu_sequence)
```

## For Sequence Completion:
```python
from transformers import T5ForConditionalGeneration, T5Tokenizer
from transformers.generation import GenerationConfig
import torch

sequence = "MDTAYPREDTRAPTPSKAGAHTALTLGAPHPPPRDHLIWSVFSTLYLNLCCLGFLALAYSIKARDQKVVGDLEAARRFGSKAKCYNILAAMWTLVPPLLLLGLVVTGALHLARLAKDSAAFFSTKFDDADYD"

ckpt = "ElnaggarLab/ankh3-large"
tokenizer = T5Tokenizer.from_pretrained(ckpt)
# To use the sequence to sequence task using the S2S prefix:
model = T5ForConditionalGeneration.from_pretrained(ckpt).eval()


half_length = int(len(sequence) * 0.5)
s2s_sequence = "[S2S]" + sequence[:half_length]
encoded_s2s_sequence = tokenizer(s2s_sequence, add_special_tokens=True, return_tensors="pt", is_split_into_words=False)
# + 1 to account for the start of sequence token.
gen_config = GenerationConfig(min_length=half_length + 1, max_length=half_length + 1, do_sample=False, num_beams=1)
generated_sequence = model.generate(encoded_s2s_sequence["input_ids"], gen_config, )
predicted_sequence = sequence[:half_length] + tokenizer.batch_decode(generated_sequence)[0]
```