Mecanismos de difusão de energia vibracional em biomoléculas têm sido relacionadas a função, alosterismo e sinalização intramolecular. Neste trabalho nós utilizamos uma metodologia computacional para analisar o fluxo de energia em proteínas. Simulações de Dinâmica Molecular são utilizadas para o estudo de difusão térmica, provendo artifícios que não são possíveis experimentalmente: a proteína é esfriada a baixas temperaturas e apenas um resíduo é aquecido através do acoplamento de um banho térmico. Consequentemente, o calor flui do aminoácido aquecido para a proteína, revelando os caminhos da difusão da energia vibracional. Pelo fato de que proteínas termoestáveis possam ter particulares mecanismos de relaxação, distribuição e dissipação da energia vibracional, elas são sistemas interessantes para serem utilizadas por este método. Um padrão de difusão de calor de uma proteína termofílica pode ser identificado e comparado com outro de uma proteína homóloga mesofílica. Aqui estudamos um conjunto de proteínas em particular, as pertencentes à família 11 de Xilanases. O mapa de difusão térmica obtido da proteína no vácuo mostrou diferenças entre xilanases mesofílica e termofílica, e termofílica e hipertermofílica: qualquer que seja o resíduo aquecido, aminoácidos específicos respondem com alta temperatura. Esta resposta em alta energia de certas regiões é decorrente de processos de relaxação. Simulações adicionais e outras análises, como da mobilidade de cada resíduo, levam à hipótese que estas regiões de superfície possuem grande flexibilidade e uma importante interação estrutural com a água. Mapas de difusão térmica para duas proteínas homólogas, diferindo em apenas 7 mutações, sendo 6 delas no N-terminal, apresentam-se diferentes nesta região. Em especial, a mutação Ser35Glu se destaca tanto no mapa quanto em outras medidas realizadas, apresentando-se na proteína mais estável com um maior nível de mobilidade, solvatação e energia de interação. Simulações em água não resultaram em padrões de difusão diferentes, por não apresentarem processos de relaxação. Entretanto, elas evidenciaram as mesmas regiões frias para cinco xilanases da família 11, especialmente o núcleo e a região de ligação do substrato, sugerindo uma possível característica funcional de difusão de calor. Por fim, evidenciado pelos mapas, observou-se que a região do cordão de xilanases termoestáveis, em especial a hipertermoestável, é maior se comparado à proteína mesofílica. Desta forma, através de comparações entre mapas de difusão e estruturas de proteínas similares, esta metodologia pode sugerir novas abordagens em engenharia racional de proteínas com estabilidade modulada.

Vibrational energy dissipation in biomolecules have been related to function, in particular alosterism and intramolecular signaling. In this work, we use a computational method to analyze the energy flux through proteins. Molecular Dynamics simulations are used to study thermal diffusion in protein structure, in an artificial way which is not accessible experimentally: the protein is cooled down to very low temperatures and a single residue is heated by coupling a thermal bath to it. Heat flows from the heated residue to the rest of the protein, revealing the paths of vibrational energy dissipation. Since thermostable proteins may have specific mechanisms for vibrational energy relaxation, dissipation and distribution, they are interesting subjects for the application of the present methodologies. The heat dissipation patterns of thermophilic proteins can be compared to the ones of less stable structures. Here, we focus in a specif set of proteins known as Xylanases of Family 11. The thermal diffusion maps obtained for hiperthermostable, thermostable and mesphilic xylanases in vacuum were different, some of them displaying apparent thermal responses whatever the heated residue. These high temperature regions appeared because of differential structural relaxation processes in each structure. The analysis of the mobility of the structures in equilibrium simulations revealed that these regions are mobile and belong to the surface of the proteins, thus interacting significantly with water molecules. Thermal diffusion maps for homologous proteins differing in only 7 residues, being 6 of them at the N-terminal region of the proteins, were different in this region. One particular mutation was determined to be more mobile in the less thermostable protein, as well as displaying a higher solvation and stronger interaction energies with remaining protein structure. Thermal diffusion simulations in water were not able to discern any difference between the structures, particularly because the relaxation processes observed in vacuum were suppressed. Nevertheless, these maps reveal that every Xylanase display the same cold regions, which were observed to belong the protein core and the catalytic site, suggesting that thermal diffusion may have some functional role. Finally,we observed that a loop which relaxed systematically in thermostable protein resulting in high temperatures is larger than in non-thermostable structures. The addition of the loop to nonthermostable proteins make the maps equivalent. Therefore, the comparison of the thermal diffusion maps of similar structures highlight important structural differences which may be useful for providing insights into the design of proteins with modulated thermal stability.