Novos modelos para cálculo de energia livre de solvatação em simulações de dinâmica molecular

Combinando métodos computacionais da eletrostática e dinâmica molecular, este trabalho teve como objetivo descrever processos de solvatação em sistemas quimicamente importantes. Foram obtidas propriedades termodinâmicas necessárias para o entendimento do processo de solvatação. Estão descritos e testados modelos para descrever a interação soluto-solvente, possibilitando, assim, aprimorar a descrição físico-química dos processos de solvatação. Utilizaram-se programas desenvolvidos em nosso grupo e programas comerciais que permitem os cálculos de dinâmica molecular e química quântica. Uma nova abordagem para o cálculo de energia livre de solvatação foi desenvolvida proporcionando a obtenção acurada e eficiente dessa propriedade, dentro do enfoque da dinâmica molecular. Nessa nova abordagem, novas metodologias para a geração de cavidades moleculares foram propostas e avaliadas. As energias livres de solvatação obtidas estão em boa concordância com os valores experimentais.

Defeitos pontuais nos compostos intermetálicos ZrNi e Zr/sub 2/Ni estudados por dinâmica molecular

Empregamos a técnica de Dinâmica Molecular para estudar propriedades de defeitos pontuais nos compostos intermetálicos ZrNi e Zr2Ni. Descrevemos as configurações estáveis de defeitos e mecanismos de migração, assim como as energias envolvidas. Os potenciais interatômicos foram derivados do Embedded Atom Model. No intuito de levar em conta a variação de estequiometria causada pela presença de alguns tipos de defeitos em intermetálicos, apresentamos um método numérico que fornece a energia efetiva de formação de defeitos e aplicamos o método ao ZrNi e Zr2Ni. Os resultados mostraram que vacâncias são mais estáveis na sub—rede do Ni, com energia de formação de 0,83 e-0,61 eV em ZrNi e Zr2Ni, respectivamente. Vacâncias de Zr são instáveis em ambos compostos; elas decaem espontaneamente em pares anti—sítio e vacância de Ni. Configurações e energias de formação de intersticiais também foram calculadas e mostraram comportamentos similares. Em ZrNi, a migração de vacâncias ocorre preferencialmente nas direções [025] e [100], com as respectivas energias de migração 0,67 e 0,73 eV, e é um processo essencialmente bidimensional no plano (001). Em Zr,Ni, a migração de vacâncias é unidimensional, ocorrendo na direção [001], com energia de migração de 0,67 eV. Em ambos compostos a presença de defeitos de anti—sítio de Ni diminui a energia de migração da vacância de Ni em até 3 vezes e facilita a movimentação em três dimensões. Mecanismos de anel não são energeticamente eficientes em comparação com saltos diretos. As configurações estáveis de intersticiais em ambos compostos consistem em um átomo de Ni sobre o plano (001) entre dois vizinhos de Zr fora do plano. Intersticiais de Zr são instáveis e tendem a deslocar um átomo de Ni, ocupando seu sítio. Energias de deslocamento foram estudadas através de simulações de irradiação de ambos compostos. Durante o processo de colisão binária, um potencial universal ZBL foi usado para colisões a curta distância. Para distâncias intermediárias usamos um potencial de união arbitrário. Zr mostrou—se mais difícil de ser arrancado de seu sítio do que Ni. Encontramos valores de energia de deslocamento no intervalo de aproximadamente 29 eV até 546 eV. Alguns resultados experimentais são mostrados e apresentam boa concordância com os cálculos.

Simulação computacional por dinâmica molecular de propriedades de equilíbrio e espectroscópicas de misturas líquidas CS2/C6H6

Propriedades de equilíbrio e espectroscópicas são determinadas para a mistura líquida CS2/C6H6 na temperatura de 298K e nas frações molares de benzeno iguais a 0,25, 0,50 e 0,75, usando simulação computacional por Dinâmica Molecular. As interações intermoleculares são descritas em relaçãao a todos os átomos por um potencial Lennard-Jones (12/6), com as interações eletrostáticas representadas por quadrupolos pontuais, localizados nos centros de massa das mol´eculas. Um segundo potencial também é considerado, onde a distribuição de cargas na molécula é usada para descrever as interações eletrostáticas. Dados termodinâmicos e propriedades estruturais, descritas através de funções de distribuição radial e de correlação angular, são calculados. Os resultados obtidos com os dois potenciais mostram boa concordância com os dados experimentais. A análise das frações molares locais e das energias internas indica que a mistura líquida CS2/C6H6 apresenta comportamento aproximadamente ideal. As funções de distribuição radial e de correlação angular evidenciam uma configuração preferencialmente paralela entre as moléculas de CS2 e C6H6 a distâncias menores em relação a distância do máximo da função de distribução , correspondente a primeira camada de coordenação . Nesta região, a configuração perpendicularé fracamente favorecida em relação as demais orientações. Funções de correlação temporal da polarizabilidade coletiva e da polarizabilidade molecular, correspondentes aos espalhamentos Rayleigh e Raman, respectivamente, são calculadas a partir do modelo de interação dipolo induzido por dipolo de 1a ordem. As funções de correlação são divididas nas contribuições orientacional, induzida por colisão e cruzada, e nas contribuições dos componentes. Propriedades de equilíbrio como anisotropia efetiva e intensidades integradas são determinadas para as misturas e líquidos puros. Os resultados indicam que a correlação orientacional entre as moléculas de CS2 é a maior responsável pelos espectros Rayleigh e espectros Raman deste componente. A maior participação do benzeno ocorre no espalhamento Raman com uma maior correlação na contribuição induzida por colisão. A separação na escala temporal entre a dinâmica reorientacional e induzida por colisão determinada na simulação é menor em relação ao experimento. Tempos de correlação da contribuição induzida por colisão são maiores nas misturas em relação aos líquidos puros no espalhamento Rayleigh, conforme o experimento, e estão relacionados a contribuições significativas das funções de correlação entre os componentes CS2 e benzeno a tempo longo. Funções de correlação temporal do momento dipolar coletivo são determinadas a partir do modelo dipolo induzido por quadrupolo e analisadas em termos de contribuições isotrópicas, anisotrópicas e cruzadas, como também de contribuições dos componentes. A relação entre estes termos é concordante com os parâmetros moleculares de anisotropia e quadrupolo dos componentes. Os resultados da simulação indicam absorbância de excesso para as misturas em uma ordem semelhante a do experimento. O tempo de correlação maior na mistura de fração molar de benzeno igual a 0,25, também encontrado no experimento, está relacionado a elevada contribuição a tempo longo da função de correlação entre dipolos induzidos em moléculas de CS2 por moléculas de C6H6. Os espectros calculados na simulação e os espectros experimentais do espalhamento de luz despolarizado Rayleigh apresentam boa concordância, validando os modelos de potenciais de interação e de indução usados. No entanto, diferenças são encontradas no infravermelho longínquo. Modi cações no modelo simplificado de indução são propostas.

Dinâmica molecular investigada com pulsos de femto-segundos

Técnicas experimentais de óptica não-linear resolvidas no tempo são capazes, em diversas situações, de fornecer mais informações sobre a dinâmica e estrutura molecular do que técnicas sem nenhuma resolução temporal. Nesta tese investigou-se três sistemas com técnicas espectroscópicas resolvidas no tempo, utilizando pulsos ultra-curtos de um sistema LASER amplificado, baseado em cristal de Ti:Safira, e de amplificadores ópticos paramétricos. Para este estudo foram implementadas e analisadas quatro diferentes técnicas espectroscópicas resolvidas no tempo: Espalhamento Raman anti-Stokes Coerente (CARS), Espalhamento hiper-Rayleigh Resolvido no Tempo (TRHRS), Bombeio-Prova e Bombeio- Depleção-Prova. O três sistemas investigados apresentam um grau crescente de complexidade nas suas dinâmicas e interações com o meio. O primeiro dos sistemas estudados foi a dinâmica da molécula de H2 no regime de impacto, na presença de moléculas e átomos perturbadores (N2 e He). Foi possível determinar, pela primeira vez, coeficientes de alargamento para o ramo Q devido à colisões entre H2+H2 utilizando a técnica de CARS resolvido temporalmente, bem como novos coeficientes para o deslocamento das linhas. O emprego da mesma técnica em sistemas binários lançou mais luz sobre a possibilidade de existir alguma não-homogeneidade em tais tipos de colisões. O segundo experimento investigou a dinâmica de relaxação orientacional da acetonitrila em solução utilizando a nova técnica de TRHRS; esta foi desenvolvida e testada com sucesso pela primeira vez neste trabalho. Um modelamento teórico mostra que dentro do modelo de Debye para a difusão rotacional, tempos associados ao terceiro momento do tempo de difusão (τ3) e ao primeiro momento (τ1) devem ser observados. Os resultados experimentais para a acetonitrila apresentam um decaimento consistente com este modelo (τ3~850 fs), concordando também com valores derivados de resultados experimentais e teóricos obtidos por outros métodos. Além disso, também foi observado experimentalmente uma segunda componente mais rápida (~50 fs), que não pôde ser explicada dentro do modelo de difusão rotacional livre. Esta componente foi atribuída a efeitos reorientacionais coletivos em conjunto com efeitos de superposição temporal de pulsos. A interpretação do sinal desta nova técnica foi testada com sucesso realizando o experimento em tetracloreto de carbono.O terceiro sistema investigado foi a dinâmica de relaxação de energia intramolecular de carotenóides. A principal questão, da vii existência ou não de estados eletrônicos adicionais, foi abordada utilizando uma modificação da técnica de Bombeio e Prova, denominada de Bombeio-Depleção-Prova. A técnica de Bombeio-Depleção-Prova foi pela primeira vez utilizada em carotenóides, possibilitando observar dinâmicas nunca vistas por qualquer outra técnica experimental. Dentro da resolução temporal dos experimentos realizados, a presença de novos estados eletrônicos em carotenóides livres em solução é descartada. Além disso, uma nova ferramenta de análise de dados foi desenvolvida. A análise de alvos espectrais baseada em um algoritmo evolutivo trouxe mais argumentos para o modelo clássico de três níveis. A dinâmica de carotenóides livres em solução foi comparada com sua dinâmica quando inseridos no complexo fotossintético de coleta de luz (LH2) da bactéria púrpura, utilizando para isso a técnica de Bombeio e Prova. Os resultados obtidos mostram um sistema de níveis de energia mais complexo do que quando o carotenóide está livre em solução. Esta modificação devido ao ambiente é possivelmente gerada pela quebra de simetria do carotenóide quando inserido no complexo LH2, possibilitando o cruzamento entre sistemas (singlete-triplete).

Espectroscopia resolvida no tempo : caracterização de pulsos curtos, dinâmica molecular em líquidos, modelamento de luz incoerente

Nesta tese realizamos três trabalhos que envolvem diferentes aspectos da espectroscopia resolvida no tempo. No primeiro discutimos brevemente a teoria e as principais técnicas para a caracterização de pulsos curtos. Analisamos detalhamente uma destas técnicas e propusemos modi cações que possibilitaram o barateamento dos custos da montagem e, além disso, introduziram novas características que sanaram alguns problemas que a montagem original apresentava e que também possibilitaram uma melhor caracterização da própria técnica. Descrevemos cuidadosamente as condições que devem ser satisfeitas pela geometria dos feixes e pelos componentes da montagem para obter uma caracterização correta dos pulsos curtos. Também apresentamos o procedimento de calibração do sistema. Pulsos com diferentes tempos e funções de fase foram caracterizados e os resultados foram validados por testes de con abilidade da informação recuperada. O trabalho seguinte foi o estudo da dinâmica molecular em líquidos puros e em misturas através da técnica efeito Kerr óptico resolvido no tempo, que é uma técnica do tipo bombeio e prova não ressonante, usando um sistema laser Ti:Sa ra com pulsos de 170 fs, centrados em 800 nm. As moléculas estudadas foram o dissulfeto de carbono (CS2), benzeno (C6H6), alilbenzeno (C9H10) e o poliestireno (PS). A teoria necessária para descrever os resultados medidos foi desenvolvida no regime temporal. O modelo de Debye para a relaxação difusiva da anisotropia orientacional descreve a componente de tempos longa, acima de um picosegundo. Partindo do Hamiltoniano de interação, desenvolvemos a teoria da função resposta linear, chegando a uma expressão para a relaxação da polarizabilidade anisotrópica, necessária para descrever os tempos curtos (subpicosegundos). Além disso, a passagem para o regime espectral utilizando somente dados experimentais, ou seja, sem a necessidade de levar em conta modelos especí cos, também foi discutida. Os resultados mostram que os tempos difusivos tanto nos líquidos puros quanto nas misturas seguem a equação de Debye-Stokes-Einstein que prevê um aumento deste tempo para viscosidades maiores. Os tempos curtos são analisados em termos da componente não difusiva da resposta espectral associada à dinâmica molecular. As alterações do espectro foram quanti cadas e a explicação para as mudanças observadas foi dada em termos das con gurações estruturais de interação que levam a uma alteração do potencial intermolecular dentro do qual as moléculas executam movimentos libracionais. Por último, investigamos a questão do modelamento de pulsos de luz incoerente. Para isto trabalhamos com um laser de corante banda larga sem elemento de seleção espectral intracavidade. A técnica utilizada foi o espalhamento forçado de luz ao qual foi acoplado um modelador composto por uma grade, uma lente e uma máscara que alterava a função de fase espectral relativa entre os feixes formadores da grade transiente na amostra de interesse. Realizamos uma análise detalhada desta montagem e obtivemos expressões para ajustar os dados medidos. Os resultados mostram que a função de correlação pode ser alterada de forma especí ca através da escolha de determinadas funções de fase espectrais.

Modelo unidimensional do atrito em escala atômica : um estudo por dinâmica molecular

o trabalho de Krim e Windom [Phys. Rev B 38, 12184(1988)], que revelou a natureza viscosa do atrito em escala atômica, gerou uma intensa atividade tanto teórica quanto experimental. Contudo, questões fundamentais ainda se matém em aberto com opor exemplo a relação entre o coeficiente de atrito viscoso e a topologia do substrato assim como a dependência com a temperatura do substrato. Neste trabalho apresentamos os resultados, obtidos usando dinâmica molecular, para um modelo unidimensional de um sistema adsorvato/substrato. Pesquisamos diferentes relações de comensuração entre o adsorvato e o substrato assim como também a dependência do coeficiente de atrito fonônico, Nph, com a temperatura. Para todas as configurações estudadas obtivemos que o coeficiente de atrito fonônico depende quadráticamente da amplitude de corrugação do substrato, mas tem uma dependência não trivial com a razão de comensuração adsorvato/substrato. O resultado mais impressionante é que para relações de comensuração entre 0.65 e 0.9 o atrito fonônico é muito fraco. Nosso resultados podem explicar as diferencias encontradas na literatura em relação à magnitude do atrito de origem eletrônico e o de origem fonônico.

Interação heparina-antitrombina : reconhecimento molecular caracterizado por ferramentas de modelagem molecular

A heparina foi isolada no início do século XX e permanece, até os dias atuais, como um dos mais importantes e eficientes agentes terapêuticos de ação antitrombótica. Sua atividade anticoagulante deve-se à ativação da antitrombina (AT), uma serpina responsável pela inibição fisiológica de serino-proteinases plasmáticas, tais como fIIa e fXa. Os esforços no sentido da elucidação dos aspectos estruturais e dinâmicos associados ao reconhecimento molecular da heparina pela AT, em nível atômico, vêm encontrando diversas dificuldades, principalmente associadas aos compostos sacarídicos. Em decorrência de sua elevada polaridade e flexibilidade, glicosaminoglicanos como a heparina são difíceis de estudar e modelar. Soma-se a isto o fato de que os resíduos de iduronato presentes na heparina (IdoA) apresentam um incomum equilíbrio conformacional entre estados de cadeira (1C4) e bote-torcido (2SO), sendo esta última estrutura postulada como a possível conformação bioativa. Sendo assim, este trabalho apresenta um estudo de modelagem molecular do perfil conformacional da heparina, tanto em solução quanto complexada à AT, utilizando cálculos ab initio e simulações de dinâmica molecular (DM) Em decorrência da ausência de parâmetros capazes de descrever polissacarídeos nos campos de força atualmente disponíveis, cargas atômicas foram geradas, utilizando-se os esquemas de Mulliken, Löwdin e cargas ajustadas ao Potencial Eletrostático, através de cálculos quantum-mecânicos na base 6-31G** e testadas na simulação de DM da heparina. Diversas condições de simulação, tais como modelos de água, concentrações de sais e as conformações 1C4 e 2SO do IdoA, foram avaliadas de forma a identificar as melhores condições para a descrição conformacional da heparina em solução. O protocolo de DM obtido foi então utilizado no estudo do complexo AT-heparina. Os resultados obtidos estão de acordo com os dados experimentais atualmente disponíveis acerca da conformação da heparina e da contribuição energética de cada resíduo de aminoácido da AT para a formação do complexo com a heparina. A partir dos cálculos ab initio realizados, foi proposto um refinamento na estrutura da heparina determinada por RMN, enquanto que as simulações de DM permitiram reinterpretar aspectos da estrutura tridimensional da AT determinada por cristalografia de raios-X. Adicionalmente, os dados obtidos sugerem que não há requerimento conformacional para a interação do IdoA com a AT Globalmente, os dados indicam que simulações de DM podem ser utilizadas para representar adequadamente a conformação da heparina, assim como para caracterizar e quantificar suas interações com a AT. Assim sendo, propomos o uso de simulações de DM do complexo entre a AT e a heparina, ou entre a AT e compostos derivados da heparina, como uma ferramenta útil no processo de desenvolvimento de novos agentes antitrombóticos e anticoagulantes.

Influência de parâmetros moleculares em funções de correlação temporal na dinâmica de solvatação mecânica

No presente trabalho descrevemos nossos resultados relativos à investigação da dinâmica de solvatação mecânica por meio de simulações por dinâmica molecular, respeitando o regime da resposta linear, em sistemas-modelo de argônio líquido com um soluto monoatômico ou diatômico dissolvido. Estudamos sistematicamente a influência dos parâmetros moleculares dos solutos (tamanho, polarizabilidade) e da densidade frente a vários modelos de solvatação. Funções de Correlação Temporal da Energia de Solvatação foram calculadas com relação à correlações de n-corpos (n = 2; 3) distinguindo interações repulsivas e atrativas para ambos os sistemas líquidos. Também obtivemos segundas derivadas temporais dessas funções referindo-se à parcelas translacionais, rotacionais e roto-translacionais na solução do diatômico. Encontramos que funções de correlação temporal coletivas podem ser razoavelmente bem aproximadas por correlações binárias a densidades baixas e, a densidades altas, correlações ternárias tornam-se mais importantes produzindo um descorrelacionamento mais rápido das funções coletivas devido a efeitos de cancelamento parciais. As funções de correlação para interações repulsivas e atrativas exibem comportamentos dinâmicos independentes do modelo de solvatação devido a fatores de escalonamento linear que afetam apenas as amplitudes das dessas funções de correlação temporal. Em geral, os sistemas com grau de liberdade rotacional apresentam tempos de correlação mais curtos para a dinâmica coletiva e tempos de correlação mais longos para as funções binárias e ternárias. Finalmente, esse estudo mostra que os sistemas contendo o diatômico relaxam-se predominantemente por mecanismos translacionais binários em modelos de solvatação envolvendo alterações apenas na polarizabilidade do soluto, e por mecanismos rotacionais atrativos binários em modelos envolvendo alterações no comprimento de ligação.

Estudo computacional de líquidos iônicos do tipo dialquilimidazólio

Neste trabalho foram desenvolvidos parâmetros de campo de força dos ânions tetrafenilborato (BPh4 -) e hexafluorfosfato (PF6 -) dentro da metodologia AMBER para a simulação computacional por dinâmica molecular de Líquidos Iônicos formados por estes ânions e cátions do tipo dialquilimidazólio, o 1-n-butil-3-metilimidazólio (BMI+), o 1-etil-3-metilimidazólio (EMI+) e o 1,3-dimetilimidazólio (MMI+). A validação destes parâmetros foi realizada por comparação entre as freqüências dos modos normais obtido por cálculo ab initio com aquelas obtidas por mecânica molecular, juntamente com uma comparação entre as estruturas moleculares e momentos multipolares obtidos pelas duas metodologias. Seguiu-se então a validação por comparação dos resultados dos cálculos de dinâmica molecular com dados experimentais, como densidades, entalpias de vaporização, condutividade elétrica, estrutura radial e espacial e também dados de difração de nêutrons. Foi atingida uma concordância bastante grande entre dados experimentais e cálculo teórico principalmente no que diz respeito à estrutura dos Líquidos Iônicos e foi possível racionalizar em termos de tamanho do grupo alquila do cátion (n-butil, metil ou etil) e do tamanho do ânion tanto propriedades estruturais deste líquidos quanto características dinâmicas dos mesmos.

Estudo teórico da solvatação de ânions em solventes orgânicos

O presente trabalho relata a parametrização do modelo de contínuo IEF-PCM para o estudo da solvatação de ânions monovalentes em acetonitrila e N,N-dimetilformamida. Nós propomos duas metodologias de parametrização do modelo de contínuo polarizável para este estudo teórico. A primeira proposta é baseada unicamente no cálculo do termo eletrostático da energia livre e na equiparação deste termo com a energia livre de solvatação experimental. Na segunda proposta de metodologia, a parametrização do modelo de contínuo polarizável foi realizada com o auxílio das simulações de dinâmica molecular e dentro desta segunda proposta nós calculamos todos os termos que contribuem para a energia livre de solvatação, a energia de cavitação. a energia eletrostática e a energia de van der Waals. Com a aplicação da primeira metodologia, é obtido um erro médio absoluto de 2,1 kcal/mol nas energias de solvatação dos ânions em acetonitrila enquanto que o erro médio absoluto calculado para energias de solvatação dos ânions em N,N-dimetilformamida é 2,8 kcal/mol. A aplicação da segunda metodologia resulta num erro médio absoluto de 1,4 kcal/mol nas energias de solvatação dos ânions em acetonitrila e um erro médio absoluto de 1,2 kcal/mol nas energias de solvatação calculadas para os ânions em N,N-dimetilformamida. Estes resultados indicam que o modelo de contínuo polarizável IEF-PCM é adequado para o estudo dos efeitos termodinâmicos nas soluções de ânions nos solventes acetonitrila e N,N-dimetilformamida.

Estudo do potencial carcinogênico de poluentes utilizando simulação molecular

Substâncias carcinogênicas são diariamente lançadas no meio ambiente, tanto na atmosfera quanto em corpos hídricos. Estas são capazes de ligar-se a proteínas constituintes de tecidos vivos produzindo um carcinoma, cuja probabilidade de formação depende da afinidade do poluente com os grupos funcionais presentes nos substratos protéicos. Contudo, os mecanismos pelos quais ocorre a formação do carcinoma não estão totalmente esclarecidos. Alguns modelos baseados em propriedades moleculares foram formulados na tentativa de prever quais serão os mecanismos, compostos intermediários e produtos finais de reação. Porém, esses modelos apresentam sérias limitações por não levarem em conta a dinâmica do processo reativo. Para que se possa estimar os mecanismos, é preciso detectar a formação ou ruptura de ligações ao longo do tempo, o que torna necessário utilizar modelos transientes. O presente trabalho apresenta um modelo que resolve a equação de Schrödinger dependente do tempo para verificar qual o mecanismo de reação entre uma substância carcinogênica e um aminoácido. A simulação do cenário transiente proposto requer um baixo tempo de processamento e possibilita uma fácil interpretação dos resultados obtidos.

Efeitos dos ácidos propiônico e metilmalônico sobre o imunoconteúdo da subunidade de alto peso molecular dos neurofilamentos em fatias de córtex cerebral de ratos

As acidemias orgânicas são desordens metabólicas caracterizadas, bioquimicamente, pelo acúmulo de um ou mais ácidos orgânicos e, clinicamente, por disfunções neurológicas graves. Os ácidos propiônico (PA) e metilmalônico (MMA) são metabólitos presentes, em altas concentrações, nos tecidos e fluídos biológicos de pacientes afetados pelas acidemias propiônica e metilmalônica, respectivamente. Os neurofilamentos (NFs) são proteínas do citoesqueleto neuronal formados pela polimerização de três subunidades, denominadas: NF de alto peso molecular (NFH), NF de médio peso molecular (NF-M) e NF de baixo peso molecular (NF-L), com pesos moleculares aparentes de 200, 150 e 68 kDa, respectivamente. Essas proteínas são amplamente fosforiladas, sendo que seu nível de fosforilação está relacionado com a polimerização dos NFs e com a regulação das interações entre os constituintes do citoesqueleto. Neste trabalho fizemos um estudo ontogenético dos efeitos in vitro do PA e do MMA 2,5 mM sobre a imunorreatividade da subunidade NF-H em fatias de córtex cerebral de ratos de 9, 12, 17, 21 e 60 dias de idade. Para tanto utilizamos uma abordagem metodológica baseada na imunodetecção da NF-H com dois anticorpos diferentes: o anticorpo monoclonal anti NF-200 clone NE14 que se liga somente a epítopos fosforilados da NF-H e, portanto, reconhece a NF-H fosforilada; e o anticorpo monoclonal anti NF-200 clone N52 que reconhece a NF-H total, independentemente do nível de fosforilação. Nossos resultados mostraram que os metabólitos induzem alterações na imunorreatividade da NF-H presente na fração citoesquelética de maneira dependente do estágio de desenvolvimento do animal Considerando o conjunto de dados obtidos com os anticorpos NE14 e N52 na fração citoesquelética, podemos dizer que o tratamento das fatias de córtex cerebral de ratos de 9 dias aumentou o nível de fosforilação sem alterar o imunoconteúdo total; em ratos de 12 e 17 dias de vida levou a um acúmulo da subunidade NF-H na fração citoesquelética e que esta proteína está na forma fosforilada; e finalmente, em ratos de 21 e 60 dias de vida os metabólitos não alteraram os parâmetros estudados. Mostramos, também, que o acúmulo da NF-H provocado pelos metabólitos em ratos de 12 e 17 dias de vida não decorreu de um aumento na síntese da proteína, mas, provalvelmente devido a um desequilíbrio na relação NF-H solúvel/NF-H polimerizada, que favoreceu a forma polimerizada. Finalmente, mostramos um possível envolvimento de mecanismos de neurotransmissão GABAérgica e glutamatérgica para os efeitos observados em ratos de 12 e 17 dias de idade, respectivamente. Considerando que a fosforilação é um mecanismo importante para a regulação da dinâmica de polimerização dos NFs e das interações com outros componentes do citoesqueleto, sugerimos que uma desorganização na sua homeostase possa causar uma neurodegeneração, que é característica dos pacientes portadores das acidemias propiônica e metilmalônica.