Communication maps computed for homodimeric hemoglobin: computational study of water-mediated energy transport in proteins.

Frequency-resolved communication maps provide a coarse-grained picture of energy transport in nanoscale systems. We calculate communication maps for homodimeric hemoglobin from Scapharca inaequivalvis and sample them to elucidate energy transfer pathways between the binding sites and other parts of the protein with focus on the role of the cluster of water molecules at the interface between the globules. We complement analysis of communication maps with molecular simulations of energy flow. Both approaches reveal that excess energy in one heme flows mainly to regions of the interface where early hydrogen bond rearrangements occur in the allosteric transition. In particular, energy is carried disproportionately by the water molecules, consistent with the larger thermal conductivity of water compared to proteins.

Quantum energy flow and the kinetics of water shuttling between hydrogen bonding sites on trans-formanilide.

A potential energy surface for trans-formanilide (TFA)-H2O is calculated and applied to study energy flow in the complex as well as the kinetics of water shuttling between hydrogen bonding sites on TFA. In addition to the previously identified H2O-TFA(C[Double Bond]O) and H2O-TFA(NH) minima, with the water monomer bound to the C[Double Bond]O and NH groups, respectively, the new surface reveals a second local minimum with the water bound to the C[Double Bond]O group, and which lies energetically 310 cm(-1) above the previously identified H2O-TFA(C[Double Bond]O) global minimum. On this surface, the energy barrier for water shuttling from H2O-TFA(C[Double Bond]O) global minimum to H2O-TFA(N-H) is 984 cm(-1), consistent with the experimental bounds of 796 and 988 cm(-1) [J. R. Clarkson et al. Science 307, 1443 (2005)]. The ergodicity threshold of TFA is calculated to be 1450 cm(-1); for the TFA-H2O complex, the coupling to the water molecule is found to lower the ergodicity threshold to below the isomerization barrier. Energy transfer between the activated complex and the vibrational modes of TFA is calculated to be sufficiently rapid that the Rice-Ramsperger-Kassel-Marcus (RRKM) theory does not overestimate the rate of water shuttling. The possibility that the rate constant for water shuttling is higher than the RRKM theory estimate is discussed in light of the relatively high energy of the ergodicity threshold calculated for TFA.

Influence of vibrational energy flow on isomerization of flexible molecules: incorporating non-Rice-Ramsperger-Kassel-Marcus kinetics in the simulation of dipeptide isomerization.

The conformational isomerization of a dipeptide, N-acetyl-tryptophan methyl amide (NATMA), is studied computationally by including important dynamical corrections to Rice-Ramsperger-Kassel-Marcus (RRKM) theory for the transition rate between pairs of isomers. The dynamical corrections arise from incomplete or sluggish vibrational energy flow in the dipeptide, a property suggested by the mode-selective chemistry that has been observed by Dian et al. [J. Chem. Phys. 120, 133 (2004)]. We compute the extent and rate of vibrational energy flow in NATMA quantum mechanically using local random matrix theory, which we then use to correct the RRKM theory rates. The latter rates are then introduced into a master equation to study the population dynamics of the dipeptide. Incomplete or slow vibrational energy flow is found to enhance the conformational selectivity of NATMA over RRKM estimates.