Calculateur d'hydrophobicite GRAVY

Une moyenne arithmétique simple a une conséquence commode et une conséquence dangereuse. La conséquence commode est qu’une chaîne longue n’est pas plus hydrophobe qu’une chaîne courte simplement parce qu’elle est plus longue. Chaque résidu porte le même poids quelle que soit sa position ou son contexte. La conséquence dangereuse est qu’une moyenne n’a aucune mémoire de l’ordre. Une chaîne dont le premier tiers est fortement hydrophobe et dont le reste est fortement hydrophile renvoie le même score qu’une chaîne dont les résidus alternent, et les deux molécules présentent à l’eau des surfaces complètement différentes.

Les valeurs de Kyte-Doolittle sont des jugements de rang et non des énergies mesurées, et c’est précisément ce qui rend raisonnable une moyenne prise sur elles. Elles ont été attribuées pour que l’échelle aille d’un bout fortement hydrophile à l’autre bout fortement hydrophobe, et elles sont symétriques par rapport à zéro par construction. Les traiter comme un classement du comportement des résidus est défendable. Traiter la moyenne qui en résulte comme une énergie libre, une solubilité ou un coefficient de partition ne l’est pas, et le guide sur le score GRAVY dit explicitement où l’indice cesse d’avoir un sens.

Ce qu’une substitution unique peut déplacer sur le score vaut mieux d’être calculé que deviné. L’échelle de Kyte-Doolittle s’étend sur neuf points d’un extrême à l’autre, donc remplacer un résidu peut déplacer la somme d’au plus cette quantité. En divisant par la longueur, le plus grand effet qu’un résidu puisse avoir est neuf divisé par n : un peu plus de deux points sur une chaîne de quatre, et un quart de point sur une chaîne de quarante. Une substitution sur une chaîne longue est donc un faible changement du score, tandis que la différence entre deux séquences choisies pour être différentes est typiquement plusieurs fois plus grande.

Cette asymétrie est l’argument pratique pour traiter GRAVY comme un moyen de trier des séquences plutôt que comme un moyen d’en mesurer une. Un écart annoncé de 0.1 entre deux versions d’un peptide long n’est presque jamais l’écart intéressant. Un changement de signe, ou un décalage de 1 ou plus, l’est. Le calculateur d’indice de Boman fait le même constat à partir d’une autre table de coefficients, et le calculateur d’indice aliphatique le refait avec quatre résidus au lieu de vingt.

L’échec le plus net de la moyenne est l’hélice amphipathe. Une hélice alpha compte environ 3.6 résidus par tour, si bien que des résidus séparés de trois ou quatre positions se retrouvent à peu près sur la même face. Une hélice construite avec des résidus hydrophobes à ces espacements et des résidus polaires aux autres est réellement et fortement amphipathe : une face interagit avec l’eau et l’autre non. Moyennée sur la chaîne, les deux ensembles de résidus se compensent, et le score d’un peptide franchement amphipathe peut se situer très près de zéro.

Rien n’est faux dans le calcul. On demande à la moyenne une question qu’elle ne peut pas traiter, à savoir le comportement moyen d’un résidu et non la façon dont ces résidus sont disposés dans l’espace. Pour une protéine globulaire repliée, la moyenne est un raccourci raisonnable parce que l’enfouissement est réparti dans toute la structure. Pour un peptide linéaire court, c’est une heuristique. Si la question est de savoir quelle part de la chaîne est chimiquement non polaire, un compte est un meilleur outil, et le calculateur de fraction de résidus hydrophobes y répond sans moyenner quoi que ce soit.

Lu à la lettre, le score GRAVY dit une seule chose : la moyenne des valeurs d’hydropathie tabulées des résidus présents. Au-dessus de zéro, le résidu moyen préfère un environnement non polaire, et en dessous de zéro il préfère l’eau. C’est un énoncé complet et exact du résultat, et il vaut mieux le laisser aussi petit, parce que toute interprétation plus vaste qui vient s’y superposer est une inférence que le calcul ne fait pas.

En particulier, il ne dit rien de la solubilité en tant que concentration, rien de l’agrégation, rien de la liaison et rien de la structure. Deux peptides de scores identiques peuvent différer sur tous ces points. Lorsqu’une propriété de séquence doit être confrontée aux autres descripteurs de composition, c’est la fraction de résidus chargés qui change d’ordinaire l’interprétation, puisque c’est la charge et non l’hydrophobie moyenne qui gouverne le comportement dans l’eau à pH modéré. La longueur est le diviseur dans tout cela, et c’est le calculateur de longueur de peptide qui établit ce nombre.