Calcolatore dell'indice alifatico

L’indice somma la percentuale molare di quattro residui con pesi fissi e ignora del tutto gli altri sedici. L’alanina ha un peso di uno, la valina 2.9 e l’isoleucina e la leucina condividono un peso di 3.9. Ogni residuo carico, ogni residuo polare e fenilalanina, triptofano e cisteina gli sono invisibili. I pesi sono coefficienti di una regressione ricavata nel lavoro originale, non costanti fisiche, quindi codificano una relazione osservata in proteine globulari ripiegate e niente di più ampio di così. Quell’origine vale la pena ricordarla, perché fissa la popolazione che il numero descrive ed è il motivo per cui l’interpretazione si ferma dove si ferma.

Una conseguenza immediata è che la scala è limitata soltanto dalla composizione che la alimenta. Una catena di alanina pura segna esattamente 100 qualunque ne sia la lunghezza, perché la sua percentuale molare di alanina è 100 e il peso è 1. Una catena di trenta due residui fatta di otto alanine, otto valine, otto isoleucine e otto leucine segna 25 + 2.9 x 25 + 3.9 x 50, che è 292.5. Nulla in quella cifra è sbagliato; mostra che il consueto intervallo da 60 a 100 citato per le proteine reali è una proprietà del modo in cui quelle proteine tendono a essere composte, e non un tetto inscritto nell’indice.

Su una catena breve il denominatore conta più dei pesi. Una sola alanina in una catena di quattro residui è il 25 per cento, quindi AGDE segna 25, un numero che sembra rassicurantemente simile a un indice proteico e che in realtà è soltanto la sua frazione di alanina. Metti la stessa alanina in una catena di quaranta residui e hai il 2.5 per cento, quindi l’indice diventa 2.5. Il calcolatore della lunghezza del peptide è ciò che dice in quale dei due regimi ci si trova prima di leggere qualsiasi cosa nel risultato, e su una catena breve è lui a decidere quanto il numero significhi.

La formula è definita sulla percentuale molare, che è il conteggio diviso per la lunghezza e moltiplicato per cento, e non sulla frazione molare semplice. Ricalcolandola a mano con 0.25 dove dovrebbe esserci 25 si ottiene un risultato cento volte più piccolo, e l’errore è completamente silenzioso perché nell’aritmetica non fallisce nulla: AGDE esce a 0.25 invece che a 25, che sembra ancora un indice ed è sbagliato di due ordini di grandezza. Il motivo per cui sopravvive a una revisione è che la risposta sbagliata ha la stessa forma di quella giusta: un numero piccolo in un campo che di solito contiene un numero piccolo.

C’è un controllo rapido che la intercetta. Per una catena che contiene alanina e nessuno degli altri tre residui contati, l’indice deve coincidere esattamente con la percentuale molare di alanina, quindi la risposta è sempre fra 0 e 100 in quel caso. Qualsiasi cosa sopra 100 dice che la formula è stata applicata a frazioni anziché a percentuali. La strada alternativa è saltare del tutto le percentuali: moltiplicare ogni conteggio grezzo per il suo peso, sommare i quattro termini e dividere una volta sola per la lunghezza alla fine. Una divisione è molto più difficile da sbagliare di due, e i conteggi provengono direttamente dal calcolatore della composizione degli amminoacidi .

C’è un caso correlato che sembra un fallimento e non lo è. Una catena priva dei quattro residui contati restituisce zero, che è una vera risposta aritmetica e non un dato respinto. È anche la risposta meno informativa che la pagina possa dare, perché su una scala dove le proteine ripiegate reali stanno intorno alla metà delle centinaia uno zero descrive una catena costruita interamente dai residui che la formula è stata definita per ignorare. Riportarlo come valore calcolato e non come errore è il comportamento corretto, e leggerlo come un punteggio basso anziché come un’assenza di residui rilevanti è l’errore.

La correlazione dietro ai pesi è stata ricavata fra composizione e comportamento di fusione di proteine globulari ripiegate. Un peptide lineare breve non ha ripiegamento e quindi non ha comportamento di fusione che l’indice possa prevedere, per cui la lettura di stabilità non ha nulla a cui agganciarsi. Il numero in sé rimane un fedele riassunto aritmetico di quattro percentuali molari. Leggerlo come una cifra di stabilità per una catena non ripiegata è un errore di categoria e non una piccola imprecisione, ed è l’errore più facile da commettere perché la parola stabilità è stampata proprio accanto al risultato.

Oltre il ripiegamento, l’indice tace sulla carica, e la carica è ciò che governa se una catena resta in soluzione. Due peptidi con percentuali identiche di alanina, valina, isoleucina e leucina possono differire in ogni residuo ionizzabile che contengono e segnare esattamente lo stesso valore. La frazione dei residui idrofobici conta i residui non polari senza pesarli e quindi si sovrappone solo in parte, mentre la guida al punteggio GRAVY e all’idropatia di Kyte-Doolittle espone perché una media pesata per residuo è il riassunto che risponde davvero a questo tipo di variazione. La carica è la variabile che più spesso decide come si comporta un peptide in soluzione, e un indice che non la vede sta rispondendo a una domanda più stretta di quella che di solito viene posta.

La lunghezza è invisibile in un modo che prende alla sprovvista. Due catene con composizione identica segnano lo stesso valore per quanto lunghe siano, perché le percentuali sono le stesse nei due casi. È una conseguenza diretta della definizione, e significa che l’indice non è un descrittore di lunghezza. Confrontarlo fra catene di dimensioni molto diverse significa confrontare anche la composizione oltre che la dimensione, altrimenti il confronto non è alla pari, e una differenza nel risultato può non dire assolutamente nulla sulla variazione che volevi introdurre.

Il risultato è adimensionale e non è una temperatura. Non esiste una conversione pubblicata dall’indice alifatico a un punto di fusione in gradi, perché l’indice è un coefficiente di regressione applicato a una composizione e non una grandezza fisica. La sua forma utile è una classifica: fra catene di lunghezza e ripiegamento simili, un valore più alto si trova all’estremo dell’intervallo che nei dati originali era associato a un comportamento di fusione più alto. Usato in qualunque altro modo smette di portare significato, e un numero citato senza la composizione che vi sta dietro non può essere controllato da nessun altro.

Quindi l’uso onesto è comparativo e stretto. Calcolalo su due sequenze, guarda quale delle quattro percentuali è davvero cambiata, e trattala come un motivo per andare a guardare qualcos’altro anziché come un risultato in sé. Se la domanda sotto il numero è quanto sia idrofoba una catena, il calcolatore del punteggio GRAVY la risponde in modo più diretto; se è quanto fortemente una sequenza possa appoggiarsi a una superficie proteica, l’ indice di Boman è stato costruito per questo; e se è semplicemente quanta materia contiene una catena, la guida sul peso molecolare e sui calcoli dei peptidi è l’aritmetica rilevante.