Exibir código-fonte para Simulação de Enovelamento de Proteína

'''Carolina Lenzi'''

Este trabalho tem como objetivo utilizar o método de Monte Carlo para simular o enovelamento de proteínas.


== Proteínas ==

As proteínas são macromoléculas com papel fundamental em diversos processos biológicos, como o transporte de oxigênio, a catálise de reações químicas, a defesa do organismo e a formação dos hormônios. Elas são formadas por centenas de aminoácidos ligados em cadeia por ligações covalentes. Existem 20 aminoácidos proteinogênicos conhecidos que dão origem às proteínas conhecidas na natureza. Cada sequência de aminoácidos ligados covalentemente gera uma proteína diferente, e essa sequência, bem determinada, é chamada de estrutura primária da proteína. A cadeia possui uma certa flexibilidade em suas ligações, que permite a movimentação dos aminoácidos no espaço. Devido à interações não covalentes (como ligações de hidrogênio, ponte salina, interações hidrofóbicas, interações iônicas…) entre os aminoácidos e com o meio, as proteínas tendem a se enovelar, gerando regiões com padrões regulares na macromolécula, chamadas de estrutura secundária da proteína. Essas regiões regulares também tendem a se enovelar, dando origem a estrutura terciária da proteína. As estruturas terciárias são, geralmente, a conformação funcional da proteína, ou seja, sua configuração em estado biologicamente ativo, que é conhecida como conformação nativa.

O enovelamento é uma estratégia de minimização da energia e leva sempre à mesma estrutura: a forma nativa da proteína. Existem muitas estruturas terciárias possíveis para a mesma cadeia de aminoácidos, pois as ligações covalentes permitem a rotação da molécula. Se considerarmos, por exemplo, que os aminoácidos tem 3 conformações possíveis, o que é uma estimativa baixa, para uma cadeia com 100 aminoácidos, existem <math>3^{100} \approx 5\times10^{47}</math> possibilidades de estrutura terciária. Se levasse 10<sup>-13</sup>s para testar cada configuração, o tempo necessário para testar todas seria 10<sup>27</sup> anos (quase 3 vezes a idade do universo!). Entretanto, as proteínas atingem a forma nativa em fração de segundo… essa divergência é conhecida como paradoxo de Levinthal.

== Simulação de Monte Carlo ==

A simulação do enovelamento foi feita com base no livro do Giordano e envolve um modelo muito simplificado do que acontece na realidade. O código foi implementado em C e os gráficos foram gerados no gnuplot.

Consideramos uma cadeia de N aminoácidos, sorteados dentre os 20 possíveis, para montar a estrutura primária da proteína, representada por um vetor de tamanho N. Aminoácidos em posições adjacentes do vetor são considerados ligados covalentemente. A cadeia é colocada em uma rede quadrada de tamanho NxN, para permitir que a proteína esteja completamente desenovelada (esticada). Em cada posição da rede, um aminoácido tem no máximo 4 aminoácidos vizinhos mais próximos com os quais pode interagir ou estar ligado, sendo eles acima, abaixo, à direita e à esquerda. Para as interações de estabilização da proteína, consideramos as seguintes forças atuantes:

*Forças de van der Waals entre aminoácidos não ligados covalentemente: força atrativa para aminoácidos próximos, que perde o efeito com o aumento da distancia entre os aminoacidos;
*Ligações de hidrogênio: ligações entre aminoácidos próximos na rede;
*Interação com a água (hidrofilicidade e hidrofobia): aminoácidos hidrofílicos são atraídos pela água presente no meio e, por isso, tendem a manter a cadeia não enovelada. Já os hidrofóbicos repelem a água do meio e preferem uma estrutura enovelada.

Todas essas forças estão competindo no processo de enovelamento da proteína. Para o modelo, agrupamos essas interações em uma energia <math>J_{ij}</math>, associada a um par de aminoácidos vizinhos na rede, A(m) e A(n) nas posições i e j, e não ligados covalentemente, ou seja, não adjacentes na cadeia. A energia da estrutura é dada pela soma sobre todos os pares de aminoácidos da proteína:

:<math> E = \sum_{<m,n>}\delta_{m,n}J_{A(m), A(n)}</math>

onde <math>\delta_{m, n} = 1</math> se os aminoácidos m e n são vizinhos na rede e não estão ligados covalentemente, e zero caso contrário.

No livro, os autores sugerem pensar as energias <math>J[A(m), A(n)]</math> como uma matriz 20x20 contendo todas as interações entre os 20 aminoácidos considerados. Porém, como a interação depende apenas de quais são os aminoácidos do par, a matriz é simétrica <math>\left(J[A(m), A(n)] = J[A(n), A(m)]\right)</math>, e podemos guardar apenas 210 valores correspondentes aos pares distintos de aminoácidos. 

Com essas definições, podemos partir para a simulação de Monte Carlo. Inicialmente, sorteamos um aminoácido e encontramos a posição dele na rede. Cada posição tem no máximo 8 vizinhos para onde o aminoácido pode (ou não) ser movido, que são os 4 vizinhos mais próximos e as diagonais para cima e para baixo. Calculamos a energia inicial da estrutura e partimos para a tentativa de movimento: sorteamos um de seus vizinhos e analisamos se é possível mover o aminoácido. Primeiro, verificamos se a rede está vazia na posição do vizinho sorteado e, em caso positivo, verificamos se nenhuma ligação entre aminoácidos será comprimida ou esticada com o movimento. Se essas duas condições forem satisfeitas, realizamos o movimento e calculamos a energia final da estrutura. Se a energia final for menor que a inicial, aceitamos o movimento e voltamos ao primeiro passo. Se a energia final for maior que a inicial, aceitamos o movimento apenas se o fator de Boltzmann <math>e^{-\frac{\Delta E}{k_BT}}</math> for maior que um valor aleatório entre 0 e 1. Caso o movimento não seja aceito, recuperamos a estrutura inicial da proteína. Repetimos esses passos diversas vezes e, a cada tentativa de movimento, contamos um tempo de Monte Carlo.

== Resultados ==

[[Arquivo:Redeinicial 15.png|200px|thumb|right|Estrutura inicial da proteína]]
Em todas as simulações, a posição inicial da cadeia foi esticada horizontalmente no meio da matriz da rede. Os valores de <math>J_{ij}</math> foram sorteados dentro do intervalo [-4, -2], como utilizado no livro. As energias foram medidas em unidades de kB por simplicidade.

[[Arquivo:Redefinal 15 t10.png|200px|thumb|right|Estrutura da proteína após 5e5 tempo de monte carlo, T = 10]]