Temperature dependence of the heat diffusivity of proteins

Time-resolved spectroscopic mapping of vibrational energy flow in proteins: Understanding thermal diffusion at the nanoscale

Vibrational energy exchange between various degrees of freedom is critical to barrier-crossing processes in proteins. Hemeproteins are well suited for studying vibrational energy exchange in proteins because the heme group is an efficient photothermal converter. The released energy by heme following photoexcitation shows migration in a protein moiety on a picosecond timescale, which is observed using time-resolved ultraviolet resonance Raman spectroscopy. The anti-Stokes ultraviolet resonance Raman intensity of a tryptophan residue is an excellent probe for the vibrational energy in proteins, allowing the mapping of energy flow with the spatial resolution of a single amino acid residue. This Perspective provides an overview of studies on vibrational energy flow in proteins, including future perspectives for both methodologies and applications.

Learning the Hydrophobic, Hydrophilic, and Aromatic Character of Amino Acids from Thermal Relaxation and Interfacial Thermal Conductance

In this study, the thermal relaxation of the 20 naturally occurring amino acids in water and in the protein lysozyme is investigated using transient nonequilibrium molecular dynamics simulations. By modeling the thermal relaxation process, the relaxation times of the amino acids in water occurs over a time scale covering 2-5 ps. For the hydrophobic amino acids, the relaxation time is controlled by the size of the hydrocarbon side chain, while for hydrophilic amino acids, the number of hydrogen bonds does not significantly affect the time scales of the heat dissipation. Our results show that the interfacial thermal conductance at the amino acid-water interface is in the range of 40-80 MW m^-2 K^-1. Hydrophobic and aromatic amino acids tend to have a lower interfacial thermal conductance. Notably, we show that amino acids can be correlated with their thermal relaxation times and molar masses, into simply connected phases with the same hydrophilicity, hydrophobicity, and aromaticity. The thermal relaxation slows down by a factor of up to five in the protein relative to that in water. In the case of the hydrophobic amino acids in the protein lysozyme, the slow down in the thermal relaxation relative to that in water appears to be controlled primarily by the size of the side chain.

Site-selective dynamics of azidolysozyme

The spectroscopic response of and structural dynamics around all azido-modified alanine residues (AlaN₃) in lysozyme are characterized. It is found that AlaN₃ is a positionally sensitive probe for the local dynamics, covering a frequency range of ∼15 cm^-1 for the center frequency of the line shape. This is consistent with findings from selective replacements of amino acids in PDZ2, which reported a frequency span of ∼10 cm^-1 for replacements of Val, Ala, or Glu by azidohomoalanine. For the frequency fluctuation correlation functions, the long-time decay constants τ₂ range from ∼1 to ∼10 ps, which compares with experimentally measured correlation times of 3 ps. Attaching azide to alanine residues can yield dynamics that decays to zero on the few ps time scale (i.e., static component Δ₀ ∼ 0 ps^-1) or to a remaining, static contribution of ∼0.5 ps^-1 (corresponding to 2.5 cm^-1), depending on the local environment on the 10 ps time scale. The magnitude of the static component correlates qualitatively with the degree of hydration of the spectroscopic probe. Although attaching azide to alanine residues is found to be structurally minimally invasive with respect to the overall protein structure, analysis of the local hydrophobicity indicates that the hydration around the modification site differs for modified and unmodified alanine residues, respectively.

Topology, landscapes, and biomolecular energy transport

While ubiquitous, energy redistribution remains a poorly understood facet of the nonequilibrium thermodynamics of biomolecules. At the molecular level, finite-size effects, pronounced nonlinearities, and ballistic processes produce behavior that diverges from the macroscale. Here, we show that transient thermal transport reflects macromolecular energy landscape architecture through the topological characteristics of molecular contacts and the nonlinear processes that mediate dynamics. While the former determines transport pathways via pairwise interactions, the latter reflects frustration within the landscape for local conformational rearrangements. Unlike transport through small-molecule systems, such as alkanes, nonlinearity dominates over coherent processes at even quite short time- and length-scales. Our exhaustive all-atom simulations and novel local-in-time and space analysis, applicable to both theory and experiment, permit dissection of energy migration in biomolecules. The approach demonstrates that vibrational energy transport can probe otherwise inaccessible aspects of macromolecular dynamics and interactions that underly biological function.

Understanding vibrational energy transfer in macromolecules has been challenging to both theory and experiment. Here the authors use non-equilibrium molecular dynamics to reveal the relationship between heat transport in a model peptide, emergent nonlinearity, and the underlying free energy landscape.

Energy Transport Pathways in Proteins: A Non-equilibrium Molecular Dynamics Simulation Study

To facilitate the observation of biomolecular energy transport in real time and with single-residue resolution, recent experiments by Baumann et al. ( Angew. Chem. Int. Ed. 2019 , 58 , 2899 , DOI: 10.1002/anie.201812995 ) have used unnatural amino acids β-(1-azulenyl)alanine (Azu) and azidohomoalanine (Aha) to site-specifically inject and probe vibrational energy in proteins. To aid the interpretation of such experiments, non-equilibrium molecular dynamics simulations of the anisotropic energy flow in proteins TrpZip2 and PDZ3 domains are presented. On this account, an efficient simulation protocol is established that accurately mimics the excitation and probing steps of Azu and Aha. The simulations quantitatively reproduce the experimentally found cooling times of the solvated proteins at room temperature and predict that the cooling slows by a factor 2 below the glass temperature of water. In PDZ3, vibrational energy is shown to travel from the initially excited peptide ligand via a complex network of inter-residue contacts and backbone transport to distal regions of the protein. The supposed connection of these energy transport pathways with pathways of allosteric communication is discussed.

Ballistic and diffusive vibrational energy transport in molecules

Energy transport in molecules is essential for many areas of science and technology. Strong covalent bonds of a molecular backbone can facilitate the involvement of the molecule's high-frequency modes in energy transport, which, under certain conditions, makes the transport fast and efficient. We discuss such conditions and describe various transport regimes in molecules, including ballistic, diffusive, directed diffusion, and intermediate regime cases, in light of recently developed experimental and theoretical approaches.

Relevance of the Resonance Junctions on the Arnold Web to Dynamical Tunneling and Eigenstate Delocalization

We study the competition and correspondence between the classical and quantum routes to intramolecular vibrational energy redistribution (IVR) in a three degrees of freedom model effective Hamiltonian. Specifically, we focus on the classical and the quantum dynamics near the resonance junctions on the Arnold web that are formed by an intersection of independent resonances. The regime of interest models the IVR dynamics from highly excited initial states near dissociation thresholds of molecular systems wherein both classical and purely quantum, involving dynamical tunneling, routes to IVR coexist. In the vicinity of a resonance junction, classical chaos is inevitably present, and hence one expects the quantum IVR pathways to have a strong classical component as well. We show that with increasing resonant coupling strengths the classical component of IVR leads to a transition from coherent dynamical tunneling to incoherent dynamical tunneling. Furthermore, we establish that the quantum IVR dynamics can be predicted based on the structures on the classical Arnold web. In addition, we investigate the nature of the highly excited eigenstates to identify the quantum signatures of the multiplicity-2 junctions. For the parameter regimes studies herein, by projecting the eigenstates onto the Arnold web, we find that eigenstates in the vicinity of the junctions are primarily delocalized due to dynamical tunneling.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

Role of Gasotransmitters in Oxidative Stresses, Neuroinflammation, and Neuronal Repair

To date, three main gasotransmitters, that is, hydrogen sulfide (H₂S), carbon monoxide (CO), and nitric oxide (NO), have been discovered to play major bodily physiological roles. These gasotransmitters have multiple functional roles in the body including physiologic and pathologic functions with respect to the cellular or tissue quantities of these gases. Gasotransmitters were originally known to have only detrimental and noxious effects in the body but that notion has much changed with years; vast studies demonstrated that these gasotransmitters are precisely involved in the normal physiological functioning of the body. From neuromodulation, oxidative stress subjugation, and cardiovascular tone regulation to immunomodulation, these gases perform critical roles, which, should they deviate from the norm, can trigger the genesis of a number of neurodegenerative diseases such as Alzheimer's disease (AD) and Parkinson's disease (PD). The purpose of this review is to discuss at great length physical and chemical properties and physiological actions of H₂S, NO, and CO as well as shedding light on recently researched molecular targets. We particularly put emphasis on the roles in neuronal inflammation and neurodegeneration and neuronal repair.

Infrared Laser Activation of Soluble and Membrane Protein Assemblies in the Gas Phase

Collision-induced dissociation (CID) is the dominant method for probing intact macromolecular complexes in the gas phase by means of mass spectrometry (MS). The energy obtained from collisional activation is dependent on the charge state of the ion and the pressures and potentials within the instrument: these factors limit CID capability. Activation by infrared (IR) laser radiation offers an attractive alternative as the radiation energy absorbed by the ions is charge-state-independent and the intensity and time scale of activation is controlled by a laser source external to the mass spectrometer. Here we implement and apply IR activation, in different irradiation regimes, to study both soluble and membrane protein assemblies. We show that IR activation using high-intensity pulsed lasers is faster than collisional and radiative cooling and requires much lower energy than continuous IR irradiation. We demonstrate that IR activation is an effective means for studying membrane protein assemblies, and liberate an intact V-type ATPase complex from detergent micelles, a result that cannot be achieved by means of CID using standard collision energies. Notably, we find that IR activation can be sufficiently soft to retain specific lipids bound to the complex. We further demonstrate that, by applying a combination of collisional activation, mass selection, and IR activation of the liberated complex, we can elucidate subunit stoichiometry and the masses of specifically bound lipids in a single MS experiment.

Scaling Rules for Vibrational Energy Transport in Globular Proteins

Computational studies of vibrational energy flow in biomolecules have to date mapped out transport pathways on a case-by-case basis. To provide a more general approach, we derive scaling rules for vibrational energy transport in a globular protein, which are identified from extensive nonequilibrium molecular dynamics simulations of vibrational energy flow in the villin headpiece subdomain HP36. We parametrize a master equation based on inter-residue, residue-solvent, and heater-residue energy-transfer rates, which closely reproduces the results of the all-atom simulations. From that fit, two scaling rules emerge, one for energy transport along the protein backbone which relies on a diffusion model and another for energy transport between tertiary contacts, which is based on a harmonic model. Requiring only the calculation of mean and variance of relatively few atomic distances, the approach holds the potential to predict the pathways and time scales of vibrational energy flow in large proteins.

Anisotropic energy flow and allosteric ligand binding in albumin

Allosteric interactions in proteins generally involve propagation of local structural changes through the protein to a remote site. Anisotropic energy transport is thought to couple the remote sites, but the nature of this process is poorly understood. Here, we report the relationship between energy flow through the structure of bovine serum albumin and allosteric interactions between remote ligand binding sites of the protein. Ultrafast infrared spectroscopy is used to probe the flow of energy through the protein backbone following excitation of a heater dye, a metalloporphyrin or malachite green, bound to different binding sites in the protein. We observe ballistic and anisotropic energy flow through the protein structure following input of thermal energy into the flexible ligand binding sites, without local heating of the rigid helix bundles that connect these sites. This efficient energy transport mechanism enables the allosteric propagation of binding energy through the connecting helix structures.

Observing Vibrational Energy Flow in a Protein with the Spatial Resolution of a Single Amino Acid Residue

One of the challenges in physical chemistry has been understanding how energy flows in a condensed phase from the microscopic viewpoint. To address this, space-resolved information at the molecular scale is required but has been lacking due to experimental difficulties. We succeeded in the real-time mapping of the vibrational energy flow in a protein with the spatial resolution of a single amino acid residue by combining time-resolved resonance Raman spectroscopy and site-directed single-Trp mutagenesis. Anti-Stokes Raman intensities of the Trp residues at different sites exhibited different temporal evolutions, reflecting propagation of the energy released by the heme group. A classical heat transport model was not able to reproduce the entire experimental data set, showing that we need a molecular-level description to explain the energy flow in a protein. The systematic application of our general methodology to proteins with different structural motifs may provide a greatly increased understanding of the energy flow in proteins.

Response of villin headpiece-capped gold nanoparticles to ultrafast laser heating

The integrity of a small model protein, the 36-residue villin headpiece HP36, attached to gold nanoparticles (AuNP) is examined, and its response to laser excitation of the AuNPs is investigated. To that end, it is first verified by stationary IR and CD spectroscopy, together with denaturation experiments, that the folded structure of the protein is fully preserved when attached to the AuNP surface. It is then shown by time-resolved IR spectroscopy that the protein does not unfold, even upon the highest pump fluences that lead to local temperature jumps on the order of 1000 K of the phonon system of the AuNPs, since that temperature jump persists for too short a time of a few nanoseconds only to be destructive. Judged from a blue shift of the amide I band, indicating destabilized or a few broken hydrogen bonds, the protein either swells, becomes more unstructured from the termini, or changes its degree of solvation. In any case, it recovers immediately after the excess energy dissipates into the bulk solvent. The process is entirely reversible for millions of laser shots without any indication of aggregation of the protein or the AuNPs and with only a minor fraction of broken protein-AuNP thiol bonds. The work provides important cornerstones in designing laser pulse parameters for maximal heating with protein-capped AuNPs without destroying the capping layer.