Matriu de pesos posicionals
Una matriu de pesos posicionals o MPP (en anglés, position weight matrix (PWM), position-specific weight matrix (PSWM) o position-specific scoring matrix (PSSM)) és una forma de representació i predicció de motius o patrons en seqüències biològiques (compostes per nucleòtits o aminoàcits) que té com a objectiu descriure les variacions intrínseques en els seus patrons.
Per este motiu, estes matrius deriven normalment d'un conjunt de seqüències alineades que se sospita que es troben funcionalmente relacionades.
Des de la seua creació, les MPP s'han anat adaptant a distints tipos de seqüències i s'han desenrollat múltiples aproximacions per a determinar els paràmetros necessaris,[1] en lo que han terminat per constituir una part molt important de moltes ferramentes de software per al descobriment de patrons.
Context
[editar | editar còdic]La matriu de pesos posicionals va ser introduïda pel genetista americà Gary Stormo en colaboració en els seus companyers en l'any 1982,[1][2] com a alternativa a les seqüències consens.
Estes seqüències havien segut amprades anteriorment per a representar patrons en les seqüències biològiques, encara que posseïen problemes en la predicció de noves ocurrències d'estos patrons.
Les MPP s'usen comunament per al modelaje de les interaccions específiques entre l'ADN i les proteïnes.[1] No obstant, la primera volta que es varen amprar estes matrius va ser en seqüències d'ARN de E. coli, a on es varen descriure distints llocs que funcionen com a dominis d'iniciació de la traducció.[3]
En la seua creació, el matemàtic Andrzej Ehrenfeucht va propondre l'algoritme perceptrón per a que les matrius foren capaces de diferenciar entre les regions reals d'unió, d'aquelles no funcionals que posseïxen seqüències similars. Es va portar a terme un entrenament del perceptrón en dos conjunts de regions que va resultar en una matriu i un llímit capaç de distinguir entre abdós conjunts.[2][3]
L'us de la matriu per a escanejar noves seqüències que no formen part de l'entrenament ha demostrat que este método resulta més sensitivo i precís que la millor seqüència consens existent.[3]
Les ventages de l'ocupació de les matrius de pesos posicionals sobre les seqüències consens han fet d'elles un método popular per a representar patrons en seqüències biològiques i un component essencial en els algoritmes moderns en el descobriment de motius.[4]
Càlcul de la matriu de pesos posicionals (MPP)
[editar | editar còdic]Les matrius de pesos posicionals expressen els patrons inherents a una alliniació de seqüències múltiples d'un conjunt de seqüències homòlogues.
L'objectiu és trobar coincidències dites seqüències alineades i les seqüències de la base de senyes, assignant un pes major a aquelles posicions que es troben conservades respecte a aquelles que són variables.[5]
Obtenció de matrius de freqüències
[editar | editar còdic]La matriu de pesos posicionals es calcula a partir del conjunt de seqüències homòlogues alineades i té una fila per cada símbol de l'alfabet (4 files per a nucleòtits en seqüències d'ADN o 20 files per a aminoàcits en seqüències de proteïnes) i una columna per a cada posició del patró.
En el següent eixemple, s'observen N=10 seqüències alineades de ADN, en L=9 nucleòtits cada una, per a les quals volem descriure un patró:[6]
|
GAGGTAAAC |
Matriu de freqüències absolutes
[editar | editar còdic]El primer pas per a construir una matriu de pesos posicionals és crear una matriu de freqüències absolutes (F). Esta matriu es construïx contant les ocurrències observades de cada nucleòtit en cada posició concreta.
Com es pot observar en l'eixemple, en la posició 1 s'han trobat 3A, 2C, 1G i 4T. Si es realisa el mateix contaje per a totes les posicions s'obtindrà la següent matriu F d'un tamany de 4x9:
___MATH_0___[6]
___MATH_1___: El número de files és 4, perque solament existixen 4 nucleòtits possibles (A,G,C,T) per a cada posició i el número de columnes és 9 perque la llongitut de les seqüències utilisades és de 9 nucleòtits.
Matriu de freqüències relatives o provabilitats
[editar | editar còdic]Per a poder interpretar els coeficients de la matriu com a provabilitats, es dividix cada freqüència absoluta entre el número total de seqüències (N) per a l'obtenció de la matriu de freqüències relatives o provabilitats(P).
Per eixemple, per a calcular la provabilitat de trobar una A en la posició 1, es dividix la cantitat d'observades en eixa posició (3) entre el número de seqüències totals (10): 3/10 = 0,3 . Este càlcul es realisa per a tots els valors de la matriu F i s'obté la següent matriu P:
___MATH_2___[6]
Esta matriu permet calcular, per a una seqüència determinada, la provabilitat de trobar dita seqüència en la regió funcional descrita per la matriu.
Donada la seqüència S = GAGGTAAAC es realisa el següent càlcul:
___MATH_3___
Es multipliquen els coeficients de la matriu corresponents a cada nucleòtit de cada posició de la seqüència S. El resultat obtingut s'interpreta com que existix una provabilitat d'un 0.07% de que la seqüència S es trobe en un lloc funcional real.
Definició matemàtica
[editar | editar còdic]Cada element de la matriu de provabilitat ___MATH_4___ és igual a la provabilitat de trobar una lletra ___MATH_5___ en la posició ___MATH_6___ del patró. Donat un conjunt de ___MATH_7___ seqüències alineades de llongitut ___MATH_8___, els elements de la matriu ___MATH_9___ es calculen com:
___MATH_10___
a on:
___MATH_11___ és el número total de seqüències alineades
___MATH_12___: ___MATH_13___ fa referència a dites seqüències, i pren valors d'1 fins a ___MATH_14___
___MATH_15___: ___MATH_16___ fa referència a les posicions dins de les seqüències, i pren valors des d'1 fins a ___MATH_17___.
___MATH_18___: ___MATH_19___ és una de les lletres de l'alfabet, en este cas A, C, T o G
___MATH_20___ és la provabilitat de trobar la lletra ___MATH_21___ en la posició ___MATH_22___.
___MATH_23___ és una funció indicatriz en la qual ___MATH_24___ és la lletra corresponent a la posició ___MATH_25___ en la seqüència ___MATH_26___, i seguix la fòrmula:
___MATH_27___, és dir, la funció ___MATH_28___ pren el valor d'1 quan ___MATH_29___ és equivalent al valor de ___MATH_30___, i el valor de 0 quan ___MATH_31___ és diferent de ___MATH_32___.
Obtenció de la matriu de pesos posicionals
[editar | editar còdic]En la majoria dels casos, els valors de la matriu de pesos posicionals M es calculen per mig de logaritmo de la raó de verosimilitut.
Una volta obtinguda la matriu de freqüències relatives o provabilitats, els valors es normalisen segons les freqüències esperades o provabilitat a priori de ___MATH_33___ (___MATH_34___) i s'obté la raó de verosimilitut. Finalment esta raó de verosimilitut es transforma a escala logarítmica.
D'esta manera, els valors que s'obtenen representen cóm de freqüent és cada nucleòtit en cada posició respecte a l'encert. Valors superiors a 0 indiquen una freqüència major a l'esperada per encert, mentres que valors inferiors a 0 indiquen una freqüència inferior a l'esperada per encert.
Per a la conversió del valor d'en la primera posició, es deu dividir l'element d'en la posició 1 de la matriu de freqüències relatives (0,3) entre la provabilitat a priori de (0,25): 0,3/0,25=1,2 i posteriorment es transforma per mig del logaritmo neperiano: ln(1,2)=0.18.
Este càlcul es realisa per a tots els valors de la matriu F i s'obté la següent matriu M:
___MATH_35___[6]
Comentari: els valors entre paréntesis són negatius, encara que és possible que no s'observen per un error de visualisació.
Ya que es tracta d'una seqüència de nucleòtits, s'assumix que la provabilitat esperada a l'encert és de 0,25. No obstant, no sempre existix una distribució uniforme en les unitats de k; per eixemple, quan s'estudien organismes en una alta composició de CG: en este cas, les provabilitats de C i G seran majors que les d'i T.
Els valors resultants d'esta matriu permeten calcular una puntuació per a una seqüència determinada. Esta puntuació serà 0 quan siga igual de provable trobar dita seqüència en la regió funcional descrita per la matriu que en qualsevol regió aleatòria del genoma, major que 0 quan siga més provable trobar la seqüència en la regió funcional que en una regió aleatòria i menor que 0 quan siga més provable trobar la seqüència en una regió aleatòria que en la regió funcional.[6]
Donada la seqüència S = GAGGTAAAC, la seua puntuació es calcularia del següent modo:
___MATH_36___
Es sumen els coeficients de la MPP corresponents a cada nucleòtit de cada posició de la seqüència donada. El resultat obtingut és major que 0, per lo que és més provable trobar la seqüència S en un lloc funcional real que en un lloc aleatori del genoma.
Definició matemàtica
[editar | editar còdic]___MATH_37___
a on ___MATH_38___ és la matriu de pesos posicionals, ___MATH_39___ és la matriu de freqüències relatives o provabilitats i ___MATH_40___ representa la provabilitat a priori o freqüències esperades de les unitats de ___MATH_41___ (0,25 per a la seqüència de nucleòtits i 0,05 per a la seqüència d'aminoàcits si assumim una distribució uniforme).
Independència estadística
[editar | editar còdic]Cal tindre en conte que totes les matrius assumixen una independència estadística entre les posicions del patró, ya que els valors per cada posició s'han calculat independentment del restant de posicions. Açò pot resultar una llimitació en alguns casos[7] i es podria rectificar per mig de la construcció de matrius més complexes afegint dimensions adicionals que tinguen en conte les posicions adjacents.
Referències
[editar | editar còdic]- ↑ 1,0 1,1 1,2 Quantitative Biology.1(2)
- 115–130.ISSN 2095-4697.doi:10.1007/s40484-013-0012-4.Consultat el 2020-12-19.
- ↑ 2,0 2,1 Nucleic Acids Research.10(9)
- 2997–3011.ISSN 0305-1048.doi:10.1093/nar/10.9.2997.Consultat el 2020-12-19.
- ↑ 3,0 3,1 3,2 Bioinformatics.16(1)
- 16–23.ISSN 1367-4803.doi:10.1093/bioinformatics/16.1.16.Consultat el 2020-12-19.
- ↑ Bioinformatics.22(14)
- i454–i463.ISSN 1367-4803.doi:10.1093/bioinformatics/btl227.Consultat el 2020-12-20.
- ↑ Gromiha, M. Michael (2010). Protein Bioinformatics (en en), Elsevier, pp. 29–62. doi:10.1016/b978-8-1312-2297-3.50002-3. ISBN 978-81-312-2297-3.
- ↑ 6,0 6,1 6,2 6,3 6,4 Guigó, Roderic. «An Introduction to Position Specific Scoring Matrius». bioinformaticaupf.crg.eu. Consultat el 2020-12-19.
- ↑ Current Protocols in Bioinformatics.51(1)ISSN 1934-3396.doi:10.1002/0471250953.bi0215s51.Consultat el 2020-12-20.
- Este artícul conté una traducció derivada de «Matriz de pesos posicionales» de Wikipedia en castellà publicada baix la Llicència de documentació lliure de GNU i la Llicència Creative Commons Reconeiximent-CompartirIgual 4.0 Internacional.