Figure ci-contre : représentation schématique du processus de transcription puis de traduction et de la structure des molécules issues de ces processus.

Cette étude bioinformatique des toxines de serpents parcourt ces protéines jusqu'à la structure 3D en focalisant sur les cystéines impliquées dans des ponts disulfure.

ExPASy : ensemble d'applications pour l'analyse de séquences peptidiques => sélectionner "Proteins & Proteomes" dans le menu de gauche.

EMBOSS seqret : application pour la conversion de formats de fichiers.

Aller au NCBI.
Rechercher les données pour "multicinctus mRNA kappa-bungarotoxin" (fenêtre du haut) en sélectionnant "All databases" dans le menu déroulant.

A quoi correspond "Y11768".

Quelle est la longueur en nucléotides de la séquence codante ("CDS") ?

Quel codon stop est employé ?

Quelle est la vraie séquence "TATA" ?

Y11768 : numéro d'accession du fichier Genbank du gène de cette toxine

Somme des longueurs des 3 exons constitutifs de la phase de lecture ouverte, soit : (332 - 277) + (1496 - 1397) + (2143 - 2033) = 264 nucléotides

Cliquer sur le lien interne "CDS" pour surligner les exons dans la séquence globale (en bas de la page) : dernier codon du 3è exon => TAA

Cliquer sur le lien interne "regulatory" pour surligner : TATAAA

Pour combien d’acides aminés code la séquence codante ?

Pour combien d’acides aminés code la séquence du peptide signal ?

Le dernier codon est le codon STOP => 261 nucléotides = 87 codons => 87 acides aminés

Peptide signal (lien interne "sig_peptide") : 54 nucléotides = 18 codons => 18 acides aminés
Voir un cours sur les peptides signaux

Faire un schéma du gène ("TATA box", peptide signal, introns et exons avec leurs positions) sur la base des données de ce fichier.

Voir le schéma

Quels sont les mots clés de l'ontologie associés à cette molécule ?

Voir un rappel sur l'ontologie.

Dans la partie "CDS", cliquer sur le lien interne "/db_xref="GOA:P01398"
Termes de l'ontologie en bas de la page : extracellular region - modulation of process of other organisms - pathogenesis - toxin activity
Voir QuickGO

De quoi s'agit-il ?
A quoi d'autre correspond "FASTA" ?

FASTA : format de fichier de séquences
FASTA : programme d'alignement

Aller au site "ExPASy". Trouver le programme "Translate". Cliquer sur la fenêtre "Browse the resource website".

Coller la séquence FASTA dans la fenêtre et lancer l'application (bouton "TRANSLATE!").

Attention : si les traductions semblent peu plausibles, revenir à la fenêtre de soumission, supprimer tout le texte qui n'est pas la séquence proprement dite (">gi|1620372|emb|...") et refaire la traduction.

Que fait ce programme et que signifient "5'3' Frame 1", "5'3' Frame 2", ... ?

Quelle traduction vous semble correcte ?

Traduction d'une séquence nucléotidique sur 6 phases de lecture : 3 sur le brin transcrit.

La traduction a priori correcte (tout du moins plausible) est celle qui contient 87 acides aminés (voir ci-dessus).

Refaire la même opération de traduction in silico avec la séquence FASTA du gène codant la bungarotoxine (fichier Y11768).

Obtient-on un résultat équivalent ? Si non, pourquoi ?

La séquence du gène contient presque 2300 nucléotides qui englobent les séquences signal, les introns et les exons (séquence codante) : le résultat de la traduction n'a donc aucun sens.

On remarque cependant que 3 segments de traduction correspondent aux 3 exons.

Dans le cadre sous "Molecular Components in 2NBT", cliquer sur la barre orange intitulée "snake_toxin superfamily".

A quelle sous-partie du NCBI accède-t-on ?

CDD ("Conserved Domain Database")

Page dédiée à la superfamille des toxines de serpents ("cl11586: snake_toxin Superfamily")

Décrire avec quelques mots-clés le mode d'action des toxines de serpents.

Décrire les particularités du domaine structural caractéristique de ces toxines.

"binding to the nicotinic acetylcholine receptors in the postsynaptic membrane" - "preventing the binding of acetylcholine" - "blocking the excitation of muscles"

Domaine : 60 - 75 acides aminés / structure stabilisée par 4 - 5 ponts disulfure / presque complètement sous forme de feuillets β / soit monomère, soit dimère

Cliquer sur le cadre bleu "cd00206: snake_toxin". Un alignement de séquences peptidiques apparaît en bas de la nouvelle page :

• Cet alignement utilise le code à 1 lettre des acides aminés.
• L'apparence de cet alignement peut-être modifié avec les différents menus ("Row Display", "Color Bits", ..., valider avec "Reformat").

• Interpréter le code couleur des lettres et minuscules/majuscules dans cet alignement.
• Repérer les acides aminés conservés dans cet alignement.
• Lire la syntaxe PROSITE des expressions régulières pour PHI-BLAST.

Ecrire un motif consensus simple au format PROSITE qui mette en évidence certains acides aminés conservés, notamment les cystéines.

Exemple de motif :

C-x(9,15)-C-x(0,2)-G-x(3)-C-x(14,18)-G-C-x(1,3)-C-P-x(9,10)-C-C-x(4,5)-[DE]-[KLNV]-C-N

Revenir à la page "Neuronal Bungarotoxin, NMR, 10 Structures". Cliquer sur le lien "2NBT" en vis-à-vis de "PDB ID:" en haut de la page à droite.

A quelle base de données accède-t-on ?
Par quelle méthode cette structure de toxine a-t-elle été résolue ?
Dans quel journal scientifique et en quelle année ce résultat a-t-il été publié ?

PDB : "Protein DataBank"

Résonance magnétique nucléaire

Publié dans le journal scientifique "Biochemistry" en 1992

• Rechercher les lignes qui décrivent les ponts disulfures ("SSBOND").

• Interpréter la syntaxe de ces lignes : nombre de ponts disulfure, positions des cystéines impliquées dans ces ponts disulfure et chaînes polypeptidiques auxquelles ces cystéines appartiennent.

1 CYS A    3    CYS A   21
2 CYS A   14    CYS A   42
3 CYS A   27    CYS A   31
4 CYS A   46    CYS A   58
5 CYS A   59    CYS A   64

5 ponts pour la chaîne A => Idem pour la chaîne B

S'agit-il d'un monomère, d'un homodimère ou d'un hétrodimère ?
Pourquoi ?

La cystéine est l'acide aminé le moins utilisé dans les chaînes polypeptidiques.
La probabilité que 2 chaînes polypeptidiques contenant 10 cystéines (sur 87 acides aminés) aux mêmes positions soient différentes est quasiment nulle.
Il s'agit d'un homodimère.

Revenir au NCBI. Récupérer la séquence au format FASTA de la protéine dont le N° d'accession est P01398.
Quelle est cette protéine ?

• C'est le précurseur de la bungarotoxine.
• La séquence en acides aminés contient le peptide signal N-terminal.

Aller à ScanProsite (Expasy).

Quelle option choisir ?
Etape 1 : coller la séquence (au format FASTA) du précurseur de la toxine dans la fenêtre.
Suivre les étapes 2 et 3 (notamment le format de sortie dans le menu déroulant).
Lancer le programme. Que recherche-t-il ?

Option 1 : soumission d'une séquence pour rechercher dans la collection de motifs.

ScanProsite recherche, dans la séquence soumise, les correspondances avec les motifs signatures des séquences de protéines de la base de données InterPro.
Cette recherche s'effectue avec le programme InterProScan.

Voir un cours sur les motifs et InterPro

Examiner la page de résultats.

• Si "Graphical view" a été choisi dans le menu déroulant de l'étape 3, cliquer sur le carré vert ou sur le lien "PS00272" (SNAKE_TOXIN Snake toxins signature).
• Autre possibilité : aller à "PDOC00245".

Quel est le motif signature des toxines de serpents ?

Enregistrer le motif obtenu.

Quel acide aminé particulier en fait partie ? Pourquoi ?

SNAKE_TOXIN, PS00272; Snake toxins signature
Consensus pattern:
G-C-x(1,3)-C-P-x(8,10)-C-C-x(2)-[PDEN]

La cystéine puisque les ponts disulfures, et la structure 3D qui en découle, sont la signature de ce type de toxine.

On compare la structure 1LSI de la neurotoxine LSIII du venin du serpent de mer Laticauda semifasciata avec la structure 2NBT (Bungarotoxine).

Aller à la page du calcul du RMSD ("Root Mean Square Deviation") entre 2 structures de protéines.

• Entrer les codes PDB : 2NBT et 1LSI.
• Entrer les positions : 1 et 67 (premier et dernier acides aminés superposés).
• Remarque : la bungarotoxine mature (sans le peptide signal) est de 66 acides aminés.

Une seule chaîne de chaque protéine est représentée.

Que traduit le RMSD et que conclure de la valeur calculée ?

Superposition des carbones α des 2 structures.

RMSD = 7.64 Å. Plus la valeur RMSD est faible, meilleure est la superposition des structures, plus elles sont similaires.

Quelles sont les parties des 2 chaînes polypeptidiques qui ne se superposent le moins bien ?

Pourquoi ?

Les boucles qui relient les brins des feuillets β.

Les boucles sont très flexibles - mobiles : elles sont donc mal résolues dans la structure obtenue par RMN en solution. En d'autres termes, il existe plusieurs positions possibles pour les acides aminés de ces boucles : la structure finale tient compte d'une position moyenne.

Le pourcentage d'identité des acides aminés des 2 protéines est seulement de 8,8 % (valeur faible). Les 2 structures se superposent cependant remarquablement : qu'en conclure ?

Les structures des protéines sont plus conservées que les séquences en acides aminés.

Aller au "SignalP 6.0 Server". Coller la séquence FASTA du précurseur de la bungarotoxine (P01398).
Choisir les paramètres corrects, lancer le programme.

Que recherche ce programme ?

Interpréter les informations dans la figure de prédiction du site de coupure (page de résultats).
Quelles sont la séquence et la position les plus probables du site de coupure ?

• Recherche de site potentiel de coupure d'un petide signal. Il s'agit d'une prédiction basée sur un algorithme et une base de données de séquences de références.

• Site de coupure prédit : GYT₂₁-R₂₂T

Refaire cette recherche avec le programme "Signal-BLAST".

Le résultat est-il identique ? Pourquoi ?

Le résultat est identique.

• Les deux algorithmes doivent suivre une logique proche.
• Les jeux de données de séquences de référence pour lesquelles on a une preuve expérimentale du site de clivage sont suffisamment vastes pour qu'ils soient probablement comparables.