The application of Kirkwood-Buff theory to study hydration properties of $\alpha$-amino acids

Título: A aplicação da teoria de Kirkwood-Buff para estudar as propriedades de hidratação de $\alpha$-aminoácidos
Autores: Z. Štefanič e B. Hribar-Lee
Publicação: Condensed Matter Physics, Vol. 29, No. 1, 13801 (2026)

1. Problema e Contexto

A estabilidade conformacional e a função das proteínas dependem criticamente de interações não covalentes (eletrostáticas, hidrofóbicas e ligações de hidrogênio) fortemente influenciadas pelo ambiente solvente. Devido à complexidade das proteínas, aminoácidos são frequentemente utilizados como sistemas modelo para entender essas interações. Embora existam estudos sobre a termodinâmica de aminoácidos em soluções eletrolíticas, há uma lacuna no conhecimento sobre como os diferentes tipos de cadeias laterais influenciam a hidratação de aminoácidos em água pura de forma comparativa e abrangente. O objetivo deste trabalho foi investigar as propriedades volumétricas e de hidratação de sete aminoácidos proteicos variados, utilizando uma combinação de dados experimentais e modelagem teórica baseada na Teoria de Kirkwood-Buff (KB).

2. Metodologia

Dados Experimentais:

• Amostras: Foram estudados sete $\alpha$-aminoácidos: arginina, aspartato, glutamato, glicina, lisina, serina e triptofano.
• Medições: As densidades das soluções aquosas foram medidas em quatro temperaturas (278,15 K a 308,15 K) utilizando um densímetro de alta precisão (DSA 5000 M).
• Cálculos: A partir das densidades, foram calculados os volumes molares aparentes ($V_\phi$) e, por extrapolação linear para diluição infinita, os volumes molares padrão ($V^\circ_\phi$).
• Hidratação: Os números de hidratação ($n_H$) foram calculados utilizando o modelo de Millero, baseado em dados de densidade e volumes molares intrínsecos.

Abordagem Teórica e Modelagem:

• Teoria: Foi aplicada a Teoria de Kirkwood-Buff (KB) combinada com a equação integral de Ornstein-Zernike (OZ) com fechamento da cadeia hiperconectada (HNC).
• Modelo Molecular: Utilizou-se um modelo "coarse-grained" (granulado) onde cada aminoácido e a água são representados por uma única "conta" (bead) isotrópica de Lennard-Jones (LJ).
• Parâmetros: As interações foram definidas por potenciais de LJ, com parâmetros de tamanho ($\sigma$) e energia ($\varepsilon$) ajustados para reproduzir a estrutura de curto alcance da água (modelo SPC/E) e as propriedades dos aminoácidos.
• Cálculos: Foram calculadas as funções de distribuição radial ($g_{ij}(r)$) e, subsequentemente, os integrais de Kirkwood-Buff ($G_{ij}$), que permitem obter os volumes molares padrão teoricamente na diluição infinita.

3. Principais Contribuições e Resultados

Propriedades Volumétricas Experimentais:

• Os volumes molares padrão ($V^\circ_\phi$) aumentaram linearmente com a temperatura para todos os aminoácidos estudados.
• Este aumento é atribuído à redução da eletrostricção da água ao redor dos centros zwitteriônicos dos aminoácidos em temperaturas mais elevadas, indicando um número de hidratação menor.
• As contribuições das cadeias laterais para o volume molar padrão foram determinadas em relação à glicina. Observou-se que, embora a lisina tenha uma massa molar ligeiramente inferior à do glutamato, sua cadeia lateral (grupo amino + cadeia alquílica longa) contribui mais significativamente para o aumento do volume molar do que o grupo carboxila do glutamato.

Comparação Teórico-Experimental:

• O modelo teórico baseado em KB/OZ reproduziu com excelente precisão os volumes molares padrão experimentais para cinco dos sete aminoácidos (Ser, Asp, Glu, Lys, Arg).
• Desvios e Correções:
  • Glicina: O modelo inicial subestimou o volume devido à forte camada de hidratação da glicina, que aumenta o raio de exclusão termodinâmica. Ao substituir o raio de van der Waals por um raio de hidratação efetivo ($2R_h$), o valor teórico concordou perfeitamente com o experimental.
  • Triptofano: O modelo isotrópico subestimou o volume devido à anisotropia do anel indol. O uso do raio de hidratação melhorou a concordância, mas uma discrepância residual permaneceu, refletindo a limitação de modelar uma estrutura anisotrópica complexa com uma única esfera.

Números de Hidratação:

• Volumétricos ($n_H$): Calculados experimentalmente, mostraram aumento com o tamanho do aminoácido e diminuição com o aumento da temperatura.
• Estruturais ($N_{shell}$): Calculados teoricamente a partir da integração da função de distribuição radial na primeira casca de solvatação.
• Correlação: Embora os valores absolutos de $N_{shell}$ (que incluem todas as moléculas na primeira casca geométrica) sejam muito maiores que os de $n_H$ (que é uma grandeza termodinâmica de eletrostricção), existe uma forte correlação linear entre eles ($N_{shell} \approx 5.01 \times n_H$). Isso valida o uso de ambos os indicadores para descrever aspectos complementares da hidratação: estrutura local versus resposta macroscópica.

4. Significado e Conclusões

O estudo demonstra que a Teoria de Kirkwood-Buff, quando acoplada a uma equação integral (OZ-HNC) e a um modelo molecular simplificado (coarse-grained), é uma ferramenta poderosa e transparente para conectar a estrutura microscópica da solução com dados volumétricos macroscópicos.

• Validação do Modelo: Mesmo com uma representação simplificada (uma única conta por aminoácido), o modelo fornece estimativas confiáveis e precisas das propriedades de solvatação.
• Insights Físicos: A abordagem permite decompor as contribuições das cadeias laterais e entender como a temperatura afeta a eletrostricção da água.
• Aplicabilidade: Os resultados reforçam a utilidade de aminoácidos como modelos para prever o comportamento de proteínas em solução, oferecendo uma base sólida para estudos bioquímicos e de química de soluções onde a hidratação desempenha um papel crucial.

Em suma, o trabalho estabelece uma ponte eficaz entre simulações teóricas e dados experimentais para a caracterização termodinâmica e estrutural da hidratação de aminoácidos.