GRAVY-Hydrophobizitaetsrechner

Ein einfaches arithmetisches Mittel hat eine bequeme und eine gefährliche Folge. Die bequeme Folge ist, dass eine lange Kette nicht allein deshalb hydrophober ist als eine kurze, weil sie länger ist. Jeder Rest trägt dasselbe Gewicht, unabhängig von Position und Umgebung. Die gefährliche Folge ist, dass ein Mittel kein Gedächtnis für die Reihenfolge hat. Eine Kette, deren erstes Drittel stark hydrophob und deren übriger Teil stark hydrophil ist, liefert denselben Wert wie eine, deren Reste abwechseln, und die beiden Moleküle bieten Wasser sehr unterschiedliche Oberflächen.

Die Kyte-Doolittle-Werte sind Rangurteile und keine gemessenen Energien, und das macht ein Mittel über sie überhaupt erst sinnvoll. Sie wurden so vergeben, dass die Skala an einem Ende von stark hydrophil bis am anderen von stark hydrophob läuft, und sie sind konstruktionsbedingt um null symmetrisch. Sie als eine Ordnung des Restverhaltens zu behandeln ist vertretbar. Das daraus entstehende Mittel als freie Energie, als Löslichkeit oder als Verteilungskoeffizienten zu behandeln ist es nicht, und der Leitfaden zum GRAVY-Wert sagt ausdrücklich, wo der Index aufhört sinnvoll zu sein.

Wie weit eine einzelne Substitution den Wert verschieben kann, rechnet man besser, als man es rät. Die Kyte-Doolittle-Skala umfasst neun Punkte von einem Extrem zum anderen, das Ersetzen eines Rests kann die Summe also höchstens um diesen Betrag verschieben. Teilt man durch die Länge, ist der größte Effekt, den ein Rest haben kann, neun geteilt durch n: etwas mehr als zwei Punkte bei einer Kette aus vier Resten und ein Viertelpunkt bei einer Kette aus vierzig. Eine Substitution auf einer langen Kette ist damit eine kleine Änderung des Werts, während der Unterschied zwischen zwei bewusst unterschiedlich gewählten Sequenzen typischerweise um ein Vielfaches größer ist.

Diese Asymmetrie ist das praktische Argument dafür, GRAVY als eine Art zu verstehen, Sequenzen zu sortieren, und nicht als eine Art, eine zu messen. Eine berichtete Änderung von 0,1 zwischen zwei Fassungen eines langen Peptids ist fast nie der interessante Unterschied. Ein Vorzeichenwechsel oder eine Verschiebung um 1 oder mehr ist es. Der Boman-Index-Rechner stellt denselben Punkt aus einer anderen Koeffiziententabelle auf, und der Rechner für den aliphatischen Index sagt ihn noch einmal mit vier statt zwanzig Resten.

Das deutlichste Versagen des Mittelwerts ist die amphipathische Helix. Eine Alpha-Helix hat etwa 3,6 Reste pro Windung, um drei oder vier Positionen voneinander abgetrennte Reste liegen also ungefähr auf derselben Seite. Eine Helix, die an diesen Abständen hydrophobe Reste und an den übrigen polare Reste trägt, ist tatsächlich und stark amphipathisch: Eine Seite wechselt mit Wasser Wechselwirkungen ein, die andere nicht. Über die Kette gemittelt heben sich die beiden Gruppen von Resten auf, und der Wert eines durchgehend amphipathischen Peptids kann sehr nahe an null liegen.

An der Rechnung ist nichts falsch. Dem Mittel wird eine Frage gestellt, die es nicht beantworten kann, nämlich wie die Reste räumlich angeordnet sind, und nicht, was das durchschnittliche Verhalten eines Rests ist. Für ein gefaltetes globuläres Protein ist das Mittel eine vernünftige Abkürzung, weil die Vergrabung durch die ganze Struktur verteilt ist. Für ein kurzes lineares Peptid ist es eine Heuristik. Lautet die Frage, wie viel der Kette chemisch unpolar ist, ist eine Anzahl das bessere Werkzeug, und der Rechner für den hydrophoben Restanteil beantwortet sie, ohne etwas wegzumitteln.

Wörtlich gelesen sagt der GRAVY-Wert genau eines: den Mittelwert der tabellierten Hydopathiewerte der vorhandenen Reste. Oberhalb von null bevorzugt der mittlere Rest eine unpolare Umgebung, unterhalb von null bevorzugt er Wasser. Das ist eine vollständige und zutreffende Aussage der Ausgabe, und sie sollte auch so klein bleiben, denn jede größere Auslegung, die darüber gelegt wird, ist eine Folgerung, die die Rechnung nicht zieht.

Insbesondere sagt er nichts über die Löslichkeit als Konzentration, nichts über Aggregation, nichts über Bindung und nichts über Struktur. Zwei Peptide mit identischem Wert können sich in jedem dieser Punkte unterscheiden. Wenn eine Sequenzeigenschaft gegen die anderen Zusammensetzungskennzahlen gestellt werden muss, ist der geladene Restanteil meist diejenige Größe, die die Auslegung verändert, denn bei mittlerem pH-Wert bestimmt die Ladung das Verhalten in Wasser und nicht die mittlere Hydrophobie. Die Länge ist der Teiler in allem, und der Peptidlängenrechner ist der Ort, an dem diese Zahl feststeht.