Calcolatore di idrofobicita GRAVY

Una media aritmetica semplice ha una conseguenza comoda e una pericolosa. La comoda è che una catena lunga non è più idrofobica di una corta solo perché è più lunga. Ogni residuo ha lo stesso peso indipendentemente dalla posizione o dal contesto. La pericolosa è che una media non ha memoria dell’ordine. Una catena il cui primo terzo è fortemente idrofobico e il cui resto è fortemente idrofilo restituisce lo stesso punteggio di una catena con residui alternati, e le due molecole presentano all’acqua superfici completamente diverse.

I valori di Kyte-Doolittle sono giudizi di ordinamento e non energie misurate, ed è questo che rende sensata una media su di essi. Sono stati assegnati in modo che la scala vada da fortemente idrofila a un estremo a fortemente idrofobica all’altro, e per costruzione sono simmetrici rispetto a zero. Trattarli come un ordinamento del comportamento dei residui è difendibile. Trattare la media che ne risulta come un’energia libera, una solubilità o un coefficiente di partizione non lo è, e la guida al punteggio GRAVY è esplicita su dove l’indice smette di avere significato.

Quanto una singola sostituzione possa spostare il punteggio vale la pena calcolarlo invece di indovinarlo. La scala di Kyte-Doolittle copre nove punti da un estremo all’altro, quindi sostituire un residuo può spostare la somma al massimo di quella quantità. Dividendo per la lunghezza, l’effetto massimo che un residuo può avere è nove su n: poco più di due punti su una catena di quattro e un quarto di punto su una catena di quaranta. Una sostituzione su una catena lunga è quindi una piccola variazione del punteggio, mentre la differenza fra due sequenze scelte per essere diverse è tipicamente diverse volte più grande.

Quell’asimmetria è l’argomento pratico per trattare GRAVY come un modo di ordinare le sequenze e non di misurarne una. Una variazione riportata di 0.1 fra due versioni di un peptide lungo non è quasi mai la differenza interessante. Lo è un cambio di segno, o uno spostamento di 1 o più. Il calcolatore dell’indice di Boman dice la stessa cosa partendo da una tabella di coefficienti diversa, e il calcolatore dell’indice alifatico la ripete con quattro residui invece di venti.

Il fallimento più chiaro della media è l’elica amfipatica. Un’elica alfa ha circa 3.6 residui per giro, quindi residui separati da tre o quattro posizioni cadono all’incirca sulla stessa faccia di essa. Un’elica costruita con residui idrofobici a quelle distanze e residui polari nelle altre posizioni è genuinamente e fortemente amfipatica: una faccia interagisce con l’acqua e l’altra no. Mediando sulla catena i due insiemi di residui si annullano, e il punteggio di un peptide fortemente amfipatico può stare molto vicino a zero.

Non c’è nulla di sbagliato nel calcolo. Alla media viene posta una domanda a cui non può rispondere, che è quale sia il comportamento medio di un residuo, e non come quei residui siano disposti nello spazio. Per una proteina globulare piegata la media è una scorciatoia ragionevole, perché il seppellimento è distribuito nella struttura. Per un peptide lineare breve è un’euristica. Se la domanda è quanto della catena sia chimicamente non polare, un conteggio è lo strumento migliore, e il calcolatore della frazione di residui idrofobici risponde senza mediare via nulla.

Letto alla lettera, il punteggio GRAVY dice una cosa sola: la media dei valori di idropatia tabulati dei residui presenti. Sopra zero il residuo medio preferisce un ambiente non polare, sotto zero preferisce l’acqua. È un’affermazione completa e accurata dell’uscita, e vale la pena tenerla così piccola, perché ogni interpretazione più grande che vi si aggiunga sopra è un’inferenza che il calcolo non fa.

In particolare non dice niente sulla solubilità come concentrazione, niente sull’aggregazione, niente sul legame e niente sulla struttura. Due peptidi con punteggi identici possono differire su tutti questi aspetti. Quando una proprietà di sequenza va confrontata con gli altri descrittori di composizione, di solito è la frazione di residui carichi a cambiare l’interpretazione, perché a pH moderato in acqua a governare il comportamento è la carica e non l’idrofobicità media. In tutto questo la lunghezza è il divisore, ed è il calcolatore della lunghezza del peptide a stabilire quel numero.