Calculateur de point isoelectrique (pI)

Un pI n’appartient pas à un peptide seul. Il appartient au peptide et à la table de valeurs de pKa utilisée pour raisonner sur lui, et ces deux choses sont séparables. Changez l’ensemble de pKa et la même séquence donne une réponse différente, pas une réponse légèrement différente mais un nombre différent, parce que le passage par zéro se déplace avec la courbe sur laquelle il se trouve. Les échelles publiées ne sont pas interchangeables, et aucun effort n’est fait ici pour les réconcilier. Deux outils qui annoncent des valeurs de pI qui se contredisent ne sont généralement pas en désaccord sur la chimie. Ils injectent des constantes différentes dans la même définition.

C’est aussi la raison honnête de traiter un pI comme approximatif, même lorsqu’il est annoncé à plusieurs décimales. La définition elle-même est exacte à l’intérieur du modèle, et le modèle n’est meilleur que ses constantes. Les valeurs de pKa tabulées ont été mesurées dans des conditions particulières, et l’écart entre les valeurs publiées pour un même groupe est de l’ordre d’une fraction d’unité de pH. L’explication des points isoélectriques des peptides couvre ce que la définition fixe et ce qu’elle ne fixe pas, et mérite d’être lue avant de citer un chiffre à plus de précision que la méthode ne le permet.

Il n’existe pas de forme fermée pour le pH auquel une somme de nombreux termes logistiques traverse zéro. Chaque groupe apporte une fonction rationnelle de dix élevé à une puissance, et une somme de plusieurs termes de ce genre ne peut pas être inversée algébriquement. Ce que l’on peut faire est de remarquer que la somme décroît de façon monotone lorsque le pH monte, qu’elle ne traverse donc zéro qu’au plus une fois, et que tout changement de signe encadre donc une racine unique. Encadrer puis diviser l’intervalle par deux suffit donc, et c’est la méthode employée ici.

La bissection est choisie plutôt qu’une méthode fondée sur la dérivée pour une raison simple : elle ne peut pas diverger. Un pas dérivé d’une pente peut dépasser et quitter l’intervalle, après quoi la recherche doit être relancée. La division par deux ne quitte jamais l’intervalle, ne l’élargit jamais, et converge à une vitesse connue à l’avance, si bien que le nombre de passes nécessaire pour atteindre une tolérance peut être annoncé avant même que la recherche commence. Sur une fonction lisse et monotone, la méthode ennuyeuse est la bonne. Le calculateur de charge nette fournit la fonction à résoudre.

Un pI annoncé à quatre décimales est une solution précise d’une fonction approximative, et les deux descriptions ne doivent pas se confondre. La bissection amène la charge résiduelle sous une tolérance bien inférieure à toute grandeur significative, parce que resserrer la tolérance ne coûte rien d’autre que des itérations supplémentaires. Les valeurs de pKa qui l’alimentent sont des constantes tabulées, et le vrai désaccord entre les valeurs publiées pour un même groupe est d’une fraction d’unité de pH. Un résultat cité à quatre décimales décrit donc l’arithmétique exactement et la chimie à une décimale près.

En pratique, cela signifie que les trois derniers chiffres sont un résidu arithmétique plutôt qu’une information. Deux décimales constituent déjà une précision généreuse pour la plupart des usages, et une décimale est généralement la précision sur laquelle agissent une formulation de tampon, un chiffre de fournisseur ou une méthode chromatographique. Rapporter plus n’est pas plus exact, seulement plus précis. Si un calcul en aval doit être fait à partir d’un pI, gardez-le à une décimale et laissez l’arrondi se produire tôt, plutôt que de propager quatre décimales dans le calcul suivant.

La composition de la chaîne décide de la réponse plus que la plupart des lecteurs ne l’attendent. Sur une séquence qui possède peu ou pas de chaînes latérales basiques, la seule charge positive disponible est celle du terminus N, et le passage par zéro est fixé par la région où les groupes carboxyle sont encore partiellement protonés. C’est pourquoi des peptides sans lysine, arginine ni histidine donnent légitimement un pI proche de 3, et pourquoi une valeur aussi basse est un résultat de structure et non une erreur. Le calculateur de charge en fonction du pH montre la même chose sous forme de courbe, avec le passage par zéro situé entre deux paliers rapprochés.

La même sensibilité joue dans l’autre sens. Capuchonner un terminus supprime une charge entière, ce qui déplace le passage par zéro bien plus que toute modification plausible du reste de la chaîne, et c’est le calculateur de modification terminale qu’il faut consulter pour établir ce qu’un capuchon donné fait à la masse avant de raisonner sur cet effet. Un pI ne dit rien non plus de la solubilité ni de la vitesse de migration. C’est le pH auquel la charge s’annule, ce qui est un fait unique et non une description générale du comportement d’un peptide autour de lui. Déterminer les conditions de tampon qui produisent un pH donné relève d’une autre arithmétique, prise en charge par le calculateur de dilution de tampon .