skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on October 2, 2026

Title: The Shaping of Enzymatic Free Energy Barriers through the Creation of Rate-Promoting Vibrations via Inter-Residue Cross-Talk on Multiple Time Scales
Years ago, we identified a rapid vibrational motion  termed a rate-promoting vibration (RPV) as central to the reaction coordinate of some enzymes. This study addresses two key questions for one example enzyme we have studied, lactate dehydrogenase (LDH): first, what is a lower bound on the RPV’s contribution to catalytic efficiency, and second, what is the mechanism of RPV formation via allosteric transmission. The goal is to understand how we can artificially create such a system. LDH catalyzes the interconversion of pyruvate and lactate via hydride and proton transfer. Altering the motion range between Val31 and Arg106, central residues in the promoting vibration, with a modest constraint reduces the reaction rate (through a raising of the free energy barrier) by over 3 orders of magnitude. Committor analysis shows that shorter distances in the constrained system shift the transition state toward proton transfer, while natural or longer distances favor a transition state formed in hydride transfer. PCA confirms the anticorrelated motion between Val31 and Arg106, aligning with vibrational modes to optimize the reaction path. Critically, we find that a breathing motion among alpha helices is used to create the necessary distance over which a rapid and short RPV can be effective. Network analysis reveals that Val31, Ala34, and Cys35 have higher eigenvector centrality in reactive trajectories, indicating enhanced inter-residue communication. These findings underscore RPVs as crucial modulators of enzymatic function through dynamic and allosteric mechanisms and suggest an approach to the generation of nonbiologic protein catalysts that include a promoting vibration.  more » « less
Award ID(s):
2244981
PAR ID:
10651537
Author(s) / Creator(s):
; ;
Publisher / Repository:
ACS JPC B
Date Published:
Journal Name:
The Journal of Physical Chemistry B
Volume:
129
Issue:
39
ISSN:
1520-6106
Page Range / eLocation ID:
9973 to 9982
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Conspectus: The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with strong energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate—the pathway from reactants to products—is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic-vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photo-induced electron transfer, singlet-triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into rapid decoherence—loss of phase and amplitude correlations—of vibrational coherences along promoter vibrations during a sub-picosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from Franck-Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics—rather than by the laser pulse—in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine-molybdenum complex hosting a nonreactive dinitrogen substrate. A slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck-Condon active terpyridine vibration to dinitrogen stretching vibration, energizing the dinitrogen bond. 
    more » « less
  2. We present numerically exact quantum dynamics simulations using the hierarchical equation of motion approach to investigate the resonance enhancement of chemical reactions due to the vibrational strong coupling (VSC) in polariton chemistry. The results reveal that the cavity mode acts like a “rate-promoting vibrational mode” that enhances the ground state chemical reaction rate constant when the cavity mode frequency matches the vibrational transition frequency. The exact simulation predicts that the VSC-modified rate constant will change quadratically as the light–matter coupling strength increases. When changing the cavity lifetime from the lossy limit to the lossless limit, the numerically exact results predict that there will be a turnover of the rate constant. Based on the numerical observations, we present an analytic rate theory to explain the observed sharp resonance peak of the rate profile when tuning the cavity frequency to match the quantum transition frequency of the vibrational ground state to excited states. This rate theory further explains the origin of the broadening of the rate profile. The analytic rate theory agrees with the numerical results under the golden rule limit and the short cavity lifetime limit. To the best of our knowledge, this is the first analytic theory that is able to explain the sharp resonance behavior of the VSC-modified rate profile when coupling an adiabatic ground state chemical reaction to the cavity. We envision that both the numerical analysis and the analytic theory will offer invaluable theoretical insights into the fundamental mechanism of the VSC-induced rate constant modifications in polariton chemistry. 
    more » « less
  3. The nonadiabatic states and dynamics are investigated for a linear vibronic coupling Hamiltonian with a static electronic splitting and weak off-diagonal Jahn-Teller coupling through a single vibration with a vibrational-electronic resonance. With a transformation of the electronic basis, this Hamiltonian is also applicable to the anti-correlated vibration in a symmetric homodimer with marginally strong constant off-diagonal coupling, where the non-adiabatic states and dynamics model electronic excitation energy transfer or self-exchange electron transfer. For parameters modeling a free-base naphthalocyanine, the nonadiabatic couplings are deeply quantum mechanical and depend on wavepacket width; scalar couplings are as important as the derivative couplings that are usually interpreted to depend on vibrational velocity in semiclassical curve crossing or surface hopping theories. A colored visualization scheme that fully characterizes the non-adiabatic states using the exact factorization is developed. The nonadiabatic states in this nested funnel have nodeless vibrational factors with strongly avoided zeroes in their vibrational probability densities. Vibronic dynamics are visualized through the vibrational coordinate dependent density of the time-dependent dipole moment in free induction decay. Vibrational motion is amplified by the nonadiabatic couplings, with asymmetric and anisotropic motions that depend upon the excitation polarization in the molecular frame and can be reversed by a change in polarization. This generates a vibrational quantum beat anisotropy in excess of 2/5. The amplitude of vibrational motion can be larger than that on the uncoupled potentials, and the electronic population transfer is maximized within one vibrational period. Most of these dynamics are missed by the adiabatic approximation, and some electronic and vibrational motions are completely suppressed by the Condon approximation of a coordinate-independent transition dipole between adiabatic states. For all initial conditions investigated, the initial nonadiabatic electronic motion is driven towards the lower adiabatic state, and criteria for this directed motion are discussed. 
    more » « less
  4. Unimolecular decay of the formaldehyde oxide (CH2OO) Criegee intermediate proceeds via a 1,3 ring-closure pathway to dioxirane and subsequent rearrangement and/or dissociation to many products including hydroxyl (OH) radicals that are detected. Vibrational activation of jet-cooled CH2OO with two quanta of CH stretch (17-18 kcal mol-1) leads to unimolecular decay at an energy significantly below the transition state barrier of 19.46  0.25 kcal mol-1, refined utilizing a high-level electronic structure method HEAT-345(Q)Λ. The observed unimolecular decay rate of 1.6 +/- 0.4 x 106 s-1 is two orders of magnitude slower than that predicted by statistical unimolecular reaction theory using several different models for quantum mechanical tunneling. The nonstatistical behavior originates from excitation of a CH stretch vibration that is orthogonal to the heavy atom motions along the reaction coordinate and slow intramolecular vibrational energy redistribution due to the sparse density of states. 
    more » « less
  5. Large-amplitude vibrational motion influences the rovibrational structure of molecules that tunnel between multiple wells. Reaction path (RP) Hamiltonians, and curvilinear coordinates more gen- erally, are useful for modelling pure vibrational motion in these systems and provide a practical framework for calculating accurate ab initio anharmonic vibrational energies and tunnelling split- tings with perturbation theory. These computational tools also offer the means to address rotation- vibration coupling associated with large-amplitude motion in rotating molecules. In this paper, we incorporate the reduced axis system (RAS) frameembeddingwithRPHamiltoniansandsecond-order vibrational Møller-Plesset perturbation theory (VMP2). Because the RP-RAS Hamiltonian eliminates rotation-vibration momentumcoupling everywhere along a one-dimensional reaction path, it is well suited for rovibrational VMP2 methods, the convergence of which relies critically on approximate vibration-vibration and vibration-rotation separability. The accuracy of this combined RP-RAS-VMP2 scheme is demonstrated by comparisons with numerically exact variational calculations and VMP2 parameters based on traditional Eckart embeddings for reduced-dimension models of torsional tunnelling in hydrogenperoxideandinversion tunnelling in cyclopropyl radical. Thefavourablecom- putational scaling ofVMP2makes it a promising strategy for calculating accurate tunnelling-rotation parameters for medium-sized and larger molecules in full dimensionality. 
    more » « less