Within the protein interior, where we observe various types of interactions, nonuniform flow of thermal energy occurs along the polypeptide chain and through nonbonded native contacts, leading to inhomogeneous transport efficiencies from one site to another. The folded native protein serves not merely as thermal transfer medium but, more importantly, as sophisticated molecular nanomachines in cells. Therefore, we are particularly interested in what sort of "communication" is mediated through native contacts in the folded proteins and how such features are quantitatively depicted in terms of local transport coefficients of heat currents. To address the issue, we introduced a concept of inter-residue thermal conductivity and characterized the nonuniform thermal transport properties of a small globular protein, HP36, using equilibrium molecular dynamics simulation and the Green-Kubo formula. We observed that the thermal transport of the protein was dominated by that along the polypeptide chain, while the local thermal conductivity of nonbonded native contacts decreased in the order of H-bonding > π-stacking > electrostatic > hydrophobic contacts. Furthermore, we applied machine learning techniques to analyze the molecular mechanism of protein thermal transport. As a result, the contact distance, variance in contact distance, and H-bonding occurrence probability during MD simulations are found to be the top three important determinants for predicting local thermal transport coefficients.
Read full abstract