Vibrational energy relaxation in proteins

Unraveling the electronic origin of a special feature in the triplet-minus-singlet spectra of carotenoids in natural photosystems

The influence of carotenoid triplet states on the Qy electronic transitions of chlorophylls has been observed in experiments on light-harvesting complexes over the past three decades, but the interpretation of the resulting spectral feature in the triplet minus singlet (T-S) absorption spectra of photosystems is still debated, as the physical-chemical explanation of this feature has been elusive. Here, we resolve this debate, by explaining the T-S spectra of pigment complexes over the Qy-band spectral region through a comparative study of chlorophyll-carotenoid model dyads and larger pigment complexes from the main light harvesting complex of higher plants (LHCII). This goal is achieved by combining state-of-the-art time-dependent density functional theory with analysis of the relationship between electronic properties and nuclear structure, and by comparison to the experiment. We find that the special signature in the T-S spectra of both model and natural photosystems is determined by singlet-like triplet excitations that can be described as effective singlet excitations on chlorophylls influenced by a stable electronic triplet on the carotenoid. The comparison with earlier experiments on different light-harvesting complexes confirms our theoretical interpretation of the T-S spectra in the Qy spectral region. Our results indicate an important role for the chlorophyll-carotenoid electronic coupling, which is also responsible for the fast triplet-triplet energy transfer, suggesting a fast trapping of the triplet into the relaxed carotenoid structure. The gained understanding of the interplay between the electronic and nuclear structures is potentially informative for future studies of the mechanism of photoprotection by carotenoids.

Electronic-state chaos, intramolecular electronic energy redistribution, and chemical bonding in persisting multidimensional nonadiabatic systems

We study the chaotic, huge fluctuation of electronic state, resultant intramolecular energy redistribution, and strong chemical bonding surviving the fluctuation with exceedingly long lifetimes of highly excited boron clusters. Those excited states constitute densely quasi-degenerate state manifolds. The huge fluctuation is induced by persisting multidimensional nonadiabatic transitions among the states in the manifold. We clarify the mechanism of their coexistence and its physical significance. In doing so, we concentrate on two theoretical aspects. One is quantum chaos and energy randomization, which are to be directly extracted from the properties of the total electronic wavefunctions. The present dynamical chaos takes place through frequent transitions from adiabatic states to others, thereby making it very rare for the system to find dissociation channels. This phenomenon leads to the concept of what we call intramolecular nonadiabatic electronic-energy redistribution, which is an electronic-state generaliztion of the notion of intramolecular vibrational energy redistribution. The other aspect is about the peculiar chemical bonding. We investigate it with the energy natural orbitals (ENOs) to see what kind of theoretical structures lie behind the huge fluctuation. The ENO energy levels representing the highly excited states under study appear to have four robust layers. We show that the energy layers responsible for chaotic dynamics and those for chemical bonding are widely separated from each other, and only when an event of what we call "inter-layer crossing" happens to burst can the destruction of these robust energy layers occur, resulting in molecular dissociation. This crossing event happens only rarely because of the large energy gaps between the ENO layers. It is shown that the layers of high energy composed of complex-valued ENOs induce the turbulent flow of electrons and electronic-energy in the cluster. In addition, the random and fast time-oscillations of those high energy ENOs serve as a random force on the nuclear dynamics, which can work to prevent a concentration of high nuclear kinetic energy in the dissociation channels.

Physical pictures of rotation mechanisms of F1- and V1-ATPases: Leading roles of translational, configurational entropy of water

We aim to develop a theory based on a concept other than the chemo-mechanical coupling (transduction of chemical free energy of ATP to mechanical work) for an ATP-driven protein complex. Experimental results conflicting with the chemo-mechanical coupling have recently emerged. We claim that the system comprises not only the protein complex but also the aqueous solution in which the protein complex is immersed and the system performs essentially no mechanical work. We perform statistical-mechanical analyses on V₁-ATPase (the A₃B₃DF complex) for which crystal structures in more different states are experimentally known than for F₁-ATPase (the α₃β₃γ complex). Molecular and atomistic models are employed for water and the structure of V₁-ATPase, respectively. The entropy originating from the translational displacement of water molecules in the system is treated as a pivotal factor. We find that the packing structure of the catalytic dwell state of V₁-ATPase is constructed by the interplay of ATP bindings to two of the A subunits and incorporation of the DF subunit. The packing structure represents the nonuniformity with respect to the closeness of packing of the atoms in constituent proteins and protein interfaces. The physical picture of rotation mechanism of F₁-ATPase recently constructed by Kinoshita is examined, and common points and differences between F₁- and V₁-ATPases are revealed. An ATP hydrolysis cycle comprises binding of ATP to the protein complex, hydrolysis of ATP into ADP and Pi in it, and dissociation of ADP and Pi from it. During each cycle, the chemical compounds bound to the three A or β subunits and the packing structure of the A₃B₃ or α₃β₃ complex are sequentially changed, which induces the unidirectional rotation of the central shaft for retaining the packing structure of the A₃B₃DF or α₃β₃γ complex stabilized for almost maximizing the water entropy. The torque driving the rotation is generated by water with no input of chemical free energy. The presence of ATP is indispensable as a trigger of the torque generation. The ATP hydrolysis or synthesis reaction is tightly coupled to the rotation of the central shaft in the normal or inverse direction through the water-entropy effect.

Quantum Chaos in the Dynamics of Molecules

Quantum chaos is reviewed from the viewpoint of "what is molecule?", particularly placing emphasis on their dynamics. Molecules are composed of heavy nuclei and light electrons, and thereby the very basic molecular theory due to Born and Oppenheimer gives a view that quantum electronic states provide potential functions working on nuclei, which in turn are often treated classically or semiclassically. Therefore, the classic study of chaos in molecular science began with those nuclear dynamics particularly about the vibrational energy randomization within a molecule. Statistical laws in probabilities and rates of chemical reactions even for small molecules of several atoms are among the chemical phenomena requiring the notion of chaos. Particularly the dynamics behind unimolecular decomposition are referred to as Intra-molecular Vibrational energy Redistribution (IVR). Semiclassical mechanics is also one of the main research fields of quantum chaos. We herein demonstrate chaos that appears only in semiclassical and full quantum dynamics. A fundamental phenomenon possibly giving birth to quantum chaos is "bifurcation and merging" of quantum wavepackets, rather than "stretching and folding" of the baker's transformation and the horseshoe map as a geometrical foundation of classical chaos. Such wavepacket bifurcation and merging are indeed experimentally measurable as we showed before in the series of studies on real-time probing of nonadiabatic chemical reactions. After tracking these aspects of molecular chaos, we will explore quantum chaos found in nonadiabatic electron wavepacket dynamics, which emerges in the realm far beyond the Born-Oppenheimer paradigm. In this class of chaos, we propose a notion of Intra-molecular Nonadiabatic Electronic Energy Redistribution (INEER), which is a consequence of the chaotic fluxes of electrons and energy within a molecule.

Energy Transport and Its Function in Heptahelical Transmembrane Proteins

Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β₂ adrenergic receptor (β₂AR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [J. Chem. Phys.2020, 152, 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure. Despite their overall structural similarity, the three 7TM proteins reveal quite different strategies to redistribute excess energy. While bR quickly removes the energy using the TM7 helix as a "lightning rod", Rho exhibits a rather poor energy dissipation, which might eventually require the hydrolysis of the Schiff base between the protein and the retinal chromophore to prevent overheating. Heating the ligand adrenaline of β₂AR, the resulting energy transport network of the protein is found to change significantly upon switching from the active state to the inactive state. While the energy flow may highlight aspects of the inter-residue couplings of β₂AR, it seems not particularly suited to explain allosteric phenomena.

Origin of Protein Quake: Energy Waves Conducted by a Precise Mechanical Machine

A long-standing challenge in protein biophysics is to understand protein quake in myoglobin─the structural dynamics responsible for redistributing the excess heme energy after photolysis. Despite extensive efforts, the molecular mechanism of this process remains elusive. Using the energy flow theory, we uncovered a fundamental new phenomenon: the heme energy is redistributed by sinusoidal waves with a ubiquitous fundamental frequency and two overtones. The energy waves emanate from the heme into the myoglobin backbone via a conduit of five consecutive dihedrals of the proximal histidine and then travel quickly along the backbone to reach sidechains across the protein. This mechanism is far more effective than the diffusion-based mechanism from previous studies because waves are systematic while diffusion is random. To propagate energy waves, coordinates must cooperate, resulting in collective modes that are singular vectors of the generalized work functional. These modes show task partitioning: a handful of high-energy modes generate large-scale breathing motion, which loosens up the protein matrix to enable hundreds of low-energy vibrational modes for energy transduction.

Versatile Vibrational Energy Sensors for Proteins

Vibrational energy transfer (VET) is emerging as key mechanism for protein functions, possibly playing an important role for energy dissipation, allosteric regulation, and enzyme catalysis. A deep understanding of VET is required to elucidate its role in such processes. Ultrafast VIS-pump/IR-probe spectroscopy can detect pathways of VET in proteins. However, the requirement of having a VET donor and a VET sensor installed simultaneously limits the possible target proteins and sites; to increase their number we compare six IR labels regarding their utility as VET sensors. We compare these labels in terms of their FTIR, and VET signature in VET donor-sensor dipeptides in different solvents. Furthermore, we incorporated four of these labels in PDZ3 to assess their capabilities in more complex systems. Our results show that different IR labels can be used interchangeably, allowing for free choice of the right label depending on the system under investigation and the methods available.

Information Flow and Allosteric Communication in Proteins

Based on Schreiber's work on transfer entropy, a molecular theory of nonlinear information transfer in proteins is developed. The joint distribution function for residue fluctuations is expressed in terms of tensor Hermite polynomials which conveniently separate harmonic and nonlinear contributions to information transfer. The harmonic part of information transfer is expressed as the difference between time dependent and independent mutual information. Third order nonlinearities are discussed in detail. Amount and speed of information transfer between residues, important for understanding allosteric activity in proteins, are discussed. While mutual information shows the maximum amount of information that may be transferred between two residues, it does not explain the actual amount of transfer nor the transfer rate of information. For this, dynamic equations of the system are needed. The solution of the Langevin equation and molecular dynamics trajectories are used in the present work for this purpose. Allosteric communication in Human NAD-dependent isocitrate dehydrogenase is studied as an example. Calculations show that several paths contribute collectively to information transfer. Important residues on these paths are identified. Time resolved information transfer between these residues, their amplitudes and transfer rates, which are in agreement with time resolved ultraviolet resonance Raman measurements in general, are estimated. Estimated transfer rates are in the order of 1-20 megabits per second. Information transfer from third order contributions are one to two orders of magnitude smaller than the harmonic terms, showing that harmonic analysis is a good approximation to information transfer.

Computational Study on the Thermal Conductivity of a Protein

Protein molecules are thermally fluctuating and tightly packed amino acid residues strongly interact with each other. Such interactions are characterized in terms of heat current at the atomic level. We calculated the thermal conductivity of a small globular protein, villin headpiece subdomain, based on the linear response theory using equilibrium molecular dynamics simulation. The value of its thermal conductivity was 0.3 ± 0.01 [W m^-1 K^-1], which is in good agreement with experimental and computational studies on the other proteins in the literature. Heat current along the main chain was dominated by local vibrations in the polypeptide bonds, with amide I, II, III, and A bands on the Fourier transform of the heat current autocorrelation function.

Mechanism for the rare fluctuation that powers protein conformational change

Most functional processes of biomolecules are rare events. Key to a rare event is the rare fluctuation that enables the energy activation process that precedes and powers crossing of the activation barrier. However, the physical nature of this rare fluctuation and how it enables energy activation and subsequently barrier crossing are unknown. We developed a novel metric, the reaction capacity pC, that rigorously defines the beginning and parameterizes the progress of energy activation. This enabled us to identify the rare fluctuation as a special phase-space condition that is necessary and sufficient for initiating systematic energy flow from the non-reaction coordinates into the reaction coordinates. The energy activation of a prototype biomolecular isomerization reaction is dominated by kinetic energy transferring into and accumulating in the reaction coordinates, administered by inertial forces alone. This mechanism for energy activation is fundamentally different from the mechanism suggested by Kramers theory.

Identification of the invariant manifolds of the LiCN molecule using Lagrangian descriptors

In this paper, we apply Lagrangian descriptors to study the invariant manifolds that emerge from the top of two barriers existing in the LiCN⇌LiNC isomerization reaction. We demonstrate that the integration times must be large enough compared with the characteristic stability exponents of the periodic orbit under study. The invariant manifolds manifest as singularities in the Lagrangian descriptors. Furthermore, we develop an equivalent potential energy surface with 2 degrees of freedom, which reproduces with a great accuracy previous results [F. Revuelta, R. M. Benito, and F. Borondo, Phys. Rev. E 99, 032221 (2019)2470-004510.1103/PhysRevE.99.032221]. This surface allows the use of an adiabatic approximation to develop a more simplified potential energy with solely 1 degree of freedom. The reduced dimensional model is still able to qualitatively describe the results observed with the original 2-degrees-of-freedom potential energy landscape. Likewise, it is also used to study in a more simple manner the influence on the Lagrangian descriptors of a bifurcation, where some of the previous invariant manifolds emerge, even before it takes place.

Through bonds or contacts? Mapping protein vibrational energy transfer using non-canonical amino acids

Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the β-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway.

Vibrational Energy Relaxation of Deuterium Fluoride in d-Dichloromethane: Insights from Different Potentials

Vibrationally excited deuterium fluoride (DF) formed by fluorine atom reaction with a solvent was found (Science, 2015, 347, 530) to relax rapidly (less than 10 ps) in acetonitrile-d₃ (CD₃CN) and dichloromethane-d₂ (CD₂Cl₂). However, insights into how CD₂Cl₂ facilitates this energy relaxation have so far been lacking, given the weak interaction between DF and a single CD₂Cl₂. In this work, we report the results of reactive simulations with a two-state reactive empirical valence bond (EVB) potential to study the energy deposited into nascent DF after transition-state passage and of nonequilibrium molecular dynamics simulations using multiple different potential energy functions to model the relaxation dynamics. For these second simulations, we used the standard Merck molecular force field (MMFF) potential, an MMFF-based covalent-ionic empirical valence bond (EVB) potential (EVB^CI), a newly developed potential [referred to as MMFF(rDF)] which extends upon the MMFF potential by making the DF/CD₂Cl₂ interaction depend on the value of the D-F bond stretching coordinate and by taking the anisotropic charge distribution of the solvent molecules into account, the polarizable atomic multipole optimized energetics for biomolecular applications (AMOEBA) potential, and the quantum mechanics/molecular mechanics (QM/MM) potential. The relaxation is revealed to be highly sensitive to the potential used. Neither standard MMFF nor EVB^CI reproduces the experimentally observed rapid relaxation dynamics, and they also fail to provide a good description of the interaction potential between DF and CD₂Cl₂ as calculated using CCSD(T)-F12. This is attributed to the use of a point-charge model for the solute and to failing to model the anisotropic electrostatic properties of CD₂Cl₂. The MMFF(rDF), AMOEBA, and QM/MM potentials all reproduce the CCSD(T)-F12 two-body DF---CD₂Cl₂ interaction potential rather well but only with the QM/MM approach is fast vibrational relaxation obtained (lifetimes of ∼288, ∼186, and ∼8 ps, respectively), which we attribute to differences in the solute-solvent local structure. With QM/MM, a unique "many-body" interaction pattern in which DF is in close contact with two solvent Cl atoms and more than three solvent D atoms is found, but this structure is not seen with other potentials. The QM/MM dynamics also display enhanced solute-solvent interactions with vibrationally excited DF that induce a DF band redshift and hence a resonant overlap with solvent C-D modes, which facilitate the intermolecular energy transfer. Our work also suggests that potentials used to model energy relaxation need to capture the fine structure of solute-solvent interactions and not just the two-body part.

Acceleration of catalysis in dihydrofolate reductase by transient, site-specific photothermal excitation

We have studied the role of protein dynamics in chemical catalysis in the enzyme dihydrofolate reductase (DHFR), using a pump-probe method that employs pulsed-laser photothermal heating of a gold nanoparticle (AuNP) to directly excite a local region of the protein structure and transient absorbance to probe the effect on enzyme activity. Enzyme activity is accelerated by pulsed-laser excitation when the AuNP is attached close to a network of coupled motions in DHFR (on the FG loop, containing residues 116-132, or on a nearby alpha helix). No rate acceleration is observed when the AuNP is attached away from the network (distal mutant and His-tagged mutant) with pulsed excitation, or for any attachment site with continuous wave excitation. We interpret these results within an energy landscape model in which transient, site-specific addition of energy to the enzyme speeds up the search for reactive conformations by activating motions that facilitate this search.

Mimicking Natural Photosynthesis: Designing Ultrafast Photosensitized Electron Transfer into Multiheme Cytochrome Protein Nanowires

Efficient nanomaterials for artificial photosynthesis require fast and robust unidirectional electron transfer (ET) from photosensitizers through charge-separation and accumulation units to redox-active catalytic sites. We explored the ultrafast time-scale limits of photo-induced charge transfer between a Ru(II)tris(bipyridine) derivative photosensitizer and PpcA, a 3-heme c-type cytochrome serving as a nanoscale biological wire. Four covalent attachment sites (K28C, K29C, K52C, and G53C) were engineered in PpcA enabling site-specific covalent labeling with expected donor-acceptor (DA) distances of 4-8 Å. X-ray scattering results demonstrated that mutations and chemical labeling did not disrupt the structure of the proteins. Time-resolved spectroscopy revealed three orders of magnitude difference in charge transfer rates for the systems with otherwise similar DA distances and the same number of covalent bonds separating donors and acceptors. All-atom molecular dynamics simulations provided additional insight into the structure-function requirements for ultrafast charge transfer and the requirement of van der Waals contact between aromatic atoms of photosensitizers and hemes in order to observe sub-nanosecond ET. This work demonstrates opportunities to utilize multi-heme c-cytochromes as frameworks for designing ultrafast light-driven ET into charge-accumulating biohybrid model systems, and ultimately for mimicking the photosynthetic paradigm of efficiently coupling ultrafast, light-driven electron transfer chemistry to multi-step catalysis within small, experimentally versatile photosynthetic biohybrid assemblies.

Ultrafast Dynamics and Vibrational Relaxation in Six-Coordinate Heme Proteins Revealed by Femtosecond Stimulated Raman Spectroscopy

Identifying the structural rearrangements during photoinduced reactions is a fundamental challenge for understanding from a microscopic perspective the dynamics underlying the functional mechanisms of heme proteins. Here, femtosecond stimulated Raman spectroscopy is applied to follow the ultrafast evolution of two different proteins, each bearing a six-coordinate heme with two amino acid axial ligands. By exploiting the sensitivity of Raman spectra to the structural configuration, we investigate the effects of photolysis and the binding of amino acid residues in cytochrome c and neuroglobin. By comparing the system response for different time delays and Raman pump resonances, we show how detailed properties of atomic motions and energy redistribution can be unveiled. In particular, we demonstrate substantially faster energy flow from the dissociated heme to the protein moiety in cytochrome c, which we assign to the presence of covalent heme–protein bonds.

Oxygen K-shell spectroscopy of isolated progressively solvated peptide

X-ray spectroscopy of an isolated controllably hydrated peptide: core excitation of the first solvation shell enhances peptide backbone fragmentation.

Normal mode analysis and beyond

Normal mode analysis provides a powerful tool in biophysical computations. Particularly, we shed light on its application to protein properties because they directly lead to biological functions. As a result of normal mode analysis, the protein motion is represented as a linear combination of mutually independent normal mode vectors. It has been widely accepted that the large amplitude motions throughout the entire protein molecule can be well described with a few low-frequency normal modes. Furthermore, it is possible to represent the effect of external perturbations, e.g., ligand binding, hydrostatic pressure, as the shifts of normal mode variables. Making use of this advantage, we are able to explore mechanical properties of proteins such as Young's modulus and compressibility. Within thermally fluctuating protein molecules under physiological conditions, tightly packed amino acid residues interact with each other through heat and energy exchanges. Since the structure and dynamics of protein molecules are highly anisotropic, the flow of energy and heat should also be anisotropic. Based on the harmonic approximation of the heat current operator, it is possible to analyze the communication map of a protein molecule. By using this method, the energy transfer pathways of photoactive yellow protein were calculated. It turned out that these pathways are similar to those obtained via the Green-Kubo formalism with equilibrium molecular dynamics simulations, indicating that normal mode analysis captures the intrinsic nature of the transport properties of proteins.

Ultrafast Structural Changes Decomposed from Serial Crystallographic Data

Direct visualization of electronic and molecular events during biochemical reactions is essential to mechanistic insights. This Letter presents an in-depth analysis of the serial crystallographic data sets collected by Barends and Schlichting et al. ( Science 2015 , 350 , 445 ) that probe the ligand photodissociation in carbonmonoxy myoglobin. This analysis reveals electron density changes caused by the formation of high-spin 3d atomic orbitals of the heme iron upon photolysis and their dynamic behaviors within the first few picoseconds. The heme iron is found popping out of and recoiling back into the heme plane in succession. These findings provide long-awaited visual validations for previous works using ultrafast spectroscopy and molecular dynamics simulations. Electron density variations are also found largely in the solvent during the first period of a low-frequency oscillation. This work demonstrates the importance of the analytical methods in detecting and isolating weak, transient signals of electronic changes arising from chemical reactions.

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.

Orthogonal beam ballistic backscatter stimulated Raman microscopy

When the axial gain length of a stimulated Raman microscope is less than about 40% of the emission wavelength significant dipole-like ballistic backscatter will occur. Here we analyze a scanning microscope configured with orthogonal water dipping pump and probe objectives that satisfies this criterion. The pump beam focus may be a Gaussian spot or a droplet Bessel beam which minimizes the secondary Bessel beam lobes and provides multiple simultaneous pump focal spot regions. Radial and linearly polarized pump beams enable backscattered polarized signals along both transverse axes of the probe beam. Low level Mie backscatter is the primary photon noise source which should enable rapid sub-wavelength resolution 3-dimensional imaging of label-free Raman contrast for in-vivo pathology, as well as, imaging physiologic concentrations of Raman labelled metabolites and drugs.

Spiers Memorial Lecture. Introductory lecture: the impact of structure on photoinduced processes in nucleic acids and proteins

Light is an important environmental variable and most organisms have evolved means to sense, exploit or avoid it and to repair detrimental effects on their genome. In general, light absorption is the task of specific chromophores, however other biomolecules such as oligonucleotides also do so which can result in undesired outcomes such as mutations and cancer. Given the biological importance of light-induced processes and applications for imaging, optogenetics, photodynamic therapy or photovoltaics, there is a great interest in understanding the detailed molecular mechanisms of photoinduced processes in proteins and nucleic acids. The processes are typically characterized by time-resolved spectroscopic approaches or computation, inferring structural information on transient species from stable ground state structures. Recently, however, structure determination of excited states or other short-lived species has become possible with the advent of X-ray free-electron lasers. This review gives an overview of the impact of structure on the understanding of photoinduced processes in macromolecules, focusing on systems presented at this Faraday Discussion meeting.

Raman, photoluminescence, and a.c. electrical studies of pure and l-serine doped ammonium dihydrogen phosphate single crystals: an understanding of defect chemistry in hydrogen bonding

Ammonium dihydrogen phosphate (ADP) crystals have a wide range of nonlinear optical applications.

Reaction mechanism and reaction coordinates from the viewpoint of energy flow

Reaction coordinates are of central importance for correct understanding of reaction dynamics in complex systems, but their counter-intuitive nature made it a daunting challenge to identify them. Starting from an energetic view of a reaction process as stochastic energy flows biased towards preferred channels, which we deemed the reaction coordinates, we developed a rigorous scheme for decomposing energy changes of a system, both potential and kinetic, into pairwise components. The pairwise energy flows between different coordinates provide a concrete statistical mechanical language for depicting reaction mechanisms. Application of this scheme to the C7eq → C7ax transition of the alanine dipeptide in vacuum revealed novel and intriguing mechanisms that eluded previous investigations of this well studied prototype system for biomolecular conformational dynamics. Using a cost function developed from the energy decomposition components by proper averaging over the transition path ensemble, we were able to identify signatures of the reaction coordinates of this system without requiring any input from human intuition.

Probing the communication of deoxythymidine triphosphate in HIV-1 reverse transcriptase by communication maps and interaction energy studies

We calculate communication maps for HIV-1 Reverse Transcriptase (RT) to elucidate energy transfer pathways between deoxythymidine triphosphate (dTTP) and other parts of the protein.

Resolving the Origins of Crystalline Anharmonicity Using Terahertz Time-Domain Spectroscopy and ab Initio Simulations

Anharmonicity has been shown to be an important piece of the fundamental framework that dictates numerous observable phenomena. In particular, anharmonicity is the driving force of vibrational relaxation processes, mechanisms that are integral to the proper function of numerous chemical processes. However, elucidating its origins has proven difficult due to experimental and theoretical challenges, specifically related to separating the anharmonic contributions from other unrelated effects. While no one technique is particularly suited for providing a complete picture of anharmonicity, by combining multiple complementary methods such a characterization can be made. In this study the role of individual atomic interactions on the anharmonic properties of crystalline purine, the building block of many DNA and RNA nucleobases, is studied by experimental terahertz time-domain spectroscopy and first-principles density functional theory (DFT) and ab initio molecular dynamics simulations (AIMD). In particular, the detailed vibrational information provided by the DFT calculations is used to interpret the atomic origins of anharmonic-related effects as determined by the AIMD calculations, which are in good agreement with the experimental data. The results highlight that anharmonicity is especially pronounced in the intermolecular interactions, particularly along the amine hydrogen bond coordinate, and yields valuable insight into what is similarly observed complex biosystems and crystalline solids.

Vibrational Heat Transport in Molecular Junctions

We review studies of vibrational energy transfer in a molecular junction geometry, consisting of a molecule bridging two heat reservoirs, solids or large chemical compounds. This setup is of interest for applications in molecular electronics, thermoelectrics, and nanophononics, and for addressing basic questions in the theory of classical and quantum transport. Calculations show that system size, disorder, structure, dimensionality, internal anharmonicities, contact interaction, and quantum coherent effects are factors that combine to determine the predominant mechanism (ballistic/diffusive), effectiveness (poor/good), and functionality (linear/nonlinear) of thermal conduction at the nanoscale. We review recent experiments and relevant calculations of quantum heat transfer in molecular junctions. We recount the Landauer approach, appropriate for the study of elastic (harmonic) phononic transport, and outline techniques that incorporate molecular anharmonicities. Theoretical methods are described along with examples illustrating the challenge of reaching control over vibrational heat conduction in molecules.

Kinetic Isotope Effect Provides Insight into the Vibrational Relaxation Mechanism of Aromatic Molecules: Application to Cyano-phenylalanine

Varying the reduced mass of an oscillator via isotopic substitution provides a convenient means to alter its vibrational frequency and hence has found wide applications. Herein, we show that this method can also help delineate the vibrational relaxation mechanism, using four isotopomers of the unnatural amino acid p-cyano-phenylalanine (Phe-CN) as models. In water, the nitrile stretching frequencies of these isotopomers, Phe-(12)C(14)N (1), Phe-(12)C(15)N (2), Phe-(13)C(14)N (3), and Phe-(13)C(15)N (4), are found to be equally separated by ∼27 cm(-1), whereas their vibrational lifetimes are determined to be 4.0 ± 0.2 (1), 2.2 ± 0.1 (2), 3.4 ± 0.2 (3), and 7.9 ± 0.5 ps (4), respectively. We find that an empirical relationship that considers the effective reduced mass of CN can accurately account for the observed frequency gaps, while the vibrational lifetime distribution, which suggests an intramolecular relaxation mechanism, can be rationalized by the order-specific density of states near the CN stretching frequency.

New Insights from Sum Frequency Generation Vibrational Spectroscopy into the Interactions of Islet Amyloid Polypeptides with Lipid Membranes

Studies of amyloid polypeptides on membrane surfaces have gained increasing attention in recent years. Several studies have revealed that membranes can catalyze protein aggregation and that the early products of amyloid aggregation can disrupt membrane integrity, increasing water permeability and inducing ion cytotoxicity. Nonetheless, probing aggregation of amyloid proteins on membrane surfaces is challenging. Surface-specific methods are required to discriminate contributions of aggregates at the membrane interface from those in the bulk phase and to characterize protein secondary structures in situ and in real time without the use of perturbing spectroscopic labels. Here, we review the most recent applications of sum frequency generation (SFG) vibrational spectroscopy applied in conjunction with computational modeling techniques, a joint experimental and computational methodology that has provided valuable insights into the aggregation of islet amyloid polypeptide (IAPP) on membrane surfaces. These applications show that SFG can provide detailed information about structures, kinetics, and orientation of IAPP during interfacial aggregation, relevant to the molecular mechanisms of type II diabetes. These recent advances demonstrate the promise of SFG as a new approach for studying amyloid diseases at the molecular level and for the rational drug design targeting early aggregation products on membrane surfaces.

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.

Direct observation of ultrafast collective motions in CO myoglobin upon ligand dissociation

The hemoprotein myoglobin is a model system for the study of protein dynamics. We used time-resolved serial femtosecond crystallography at an x-ray free-electron laser to resolve the ultrafast structural changes in the carbonmonoxy myoglobin complex upon photolysis of the Fe-CO bond. Structural changes appear throughout the protein within 500 femtoseconds, with the C, F, and H helices moving away from the heme cofactor and the E and A helices moving toward it. These collective movements are predicted by hybrid quantum mechanics/molecular mechanics simulations. Together with the observed oscillations of residues contacting the heme, our calculations support the prediction that an immediate collective response of the protein occurs upon ligand dissociation, as a result of heme vibrational modes coupling to global modes of the protein.

Energy exchange network of inter-residue interactions within a thermally fluctuating protein molecule: A computational study

Protein function is regulated not only by the structure but also by physical dynamics and thermal fluctuations. We have developed the computer program, CURrent calculation for proteins (CURP), for the flow analysis of physical quantities within thermally fluctuating protein media. The CURP program was used to calculate the energy flow within the third PDZ domain of the neuronal protein PSD-95, and the results were used to illustrate the energy exchange network of inter-residue interactions based on atomistic molecular dynamics simulations. The removal of the α3 helix is known to decrease ligand affinity by 21-fold without changing the overall protein structure; nevertheless, we demonstrated that the helix constitutes an essential part of the network graph.

Vibrational control of electron-transfer reactions: a feasibility study for the fast coherent transfer regime

Electron donors are connectedvialeft and right bridges to electron acceptors. Following electron-transfer initiation, the IR excitation of selected bridge vibrational modes can tune the directionality of electron transfer.