Hydrogen tunneling links protein dynamics to enzyme catalysis

Bottom-up parameterization of enzyme rate constants: Reconciling inconsistent data

Kinetic models of metabolism are promising platforms for studying complex metabolic systems and designing production strains. Given the availability of enzyme kinetic data from historical experiments and machine learning estimation tools, a straightforward modeling approach is to assemble kinetic data enzyme by enzyme until a desired scale is reached. However, this type of 'bottom up' parameterization of kinetic models has been difficult due to a number of issues including gaps in kinetic parameters, the complexity of enzyme mechanisms, inconsistencies between parameters obtained from different sources, and in vitro-in vivo differences. Here, we present a computational workflow for the robust estimation of kinetic parameters for detailed mass action enzyme models while taking into account parameter uncertainty. The resulting software package, termed MASSef (the Mass Action Stoichiometry Simulation Enzyme Fitting package), can handle standard 'macroscopic' kinetic parameters, including Km, kcat, Ki, Keq, and nh, as well as diverse reaction mechanisms defined in terms of mass action reactions and 'microscopic' rate constants. We provide three enzyme case studies demonstrating that this approach can identify and reconcile inconsistent data either within in vitro experiments or between in vitro and in vivo enzyme function. We further demonstrate how parameterized enzyme modules can be used to assemble pathway-scale kinetic models consistent with in vivo behavior. This work builds on the legacy of knowledge on kinetic behavior of enzymes by enabling robust parameterization of enzyme kinetic models at scale utilizing the abundance of historical literature data and machine learning parameter estimates.

Hydrogen Tunneling Exhibiting Unexpectedly Small Primary Kinetic Isotope Effects

Probing quantum mechanical tunneling (QMT) in chemical reactions is crucial to understanding and developing new transformations. Primary H/D kinetic isotopic effects (KIEs) beyond the semiclassical maximum values of 7-10 (room temperature) are commonly used to assess substantial QMT contributions in one-step hydrogen transfer reactions, because of the much greater QMT probability of protium vs. deuterium. Nevertheless, we report here the discovery of a reaction model occurring exclusively by H-atom QMT with residual primary H/D KIEs. 2-Hydroxyphenylnitrene, generated in N2 matrix, was found to isomerize to an imino-ketone via sequential (domino) QMT involving anti to syn OH-rotamerization (rate determining step) and [1,4]-H shift reactions. These sequential QMT transformations were also observed in the OD-deuterated sample, and unexpected primary H/D KIEs between 3 and 4 were measured at 3 to 20 K. Analogous residual primary H/D KIEs were found in the anti to syn OH-rotamerization QMT of 2-cyanophenol in a N2 matrix. Evidence strongly indicates that these intriguing isotope-insensitive QMT reactivities arise due to the solvation effects of the N2 matrix medium, putatively through coupling with the moving H/D tunneling particle. Should a similar scenario be extrapolated to conventional solution conditions, then QMT may have been overlooked in many chemical reactions.

Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine

We introduce a quantum mechanics/molecular mechanics semiclassical method for studying the solvation process of molecules in water at the nuclear quantum mechanical level with atomistic detail. We employ it in vibrational spectroscopy calculations because this is a tool that is very sensitive to the molecular environment. Specifically, we look at the vibrational spectroscopy of thymidine in liquid water. We find that the C═O frequency red shift and the C═C frequency blue shift, experienced by thymidyne upon solvation, are mainly due to reciprocal polarization effects, that the molecule and the water solvent exert on each other, and nuclear zero-point energy effects. In general, this work provides an accurate and practical tool to study quantum vibrational spectroscopy in solution and condensed phase, incorporating high-level and computationally affordable descriptions of both electronic and nuclear problems.

Ultralow-Friction at Cryogenic Temperature Induced by Hydrogen Correlated Quantum Effect

Temperature is one of the governing factors affecting friction of solids. Undesired high friction state has been generally reported at cryogenic temperatures due to the prohibition of thermally activated processes, following conventional Arrhenius equation. This has brought huge difficulties to lubrication at extremely low temperatures in industry. Here, the study uncovers a hydrogen-correlated sub-Arrhenius friction behavior in hydrogenated amorphous carbon (a-C:H) film at cryogenic temperatures, and a stable ultralow-friction over a wide temperature range (103-348 K) is achieved. This is attributed to hydrogen-transfer-induced mild structural ordering transformation, confirmed by machine-learning-based molecular dynamics simulations. The anomalous sub-Arrhenius temperature dependence of structural ordering transformation rate is well-described by a quantum mechanical tunneling (QMT) modified Arrhenius model, which is correlated with quantum delocalization of hydrogen in tribochemical reactions. This work reveals a hydrogen-correlated friction mechanism overcoming the Arrhenius temperature dependence and provides a new pathway for achieving ultralow friction under cryogenic conditions.

Structural Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects of the NADH/NAD+ Model Reactions in Acetonitrile: Charge-Transfer Complex Tightness Is a Key

It has recently frequently been found that the kinetic isotope effect (KIE) is independent of temperature (T) in H-tunneling reactions in enzymes but becomes dependent on T in their mutants. Many enzymologists found that the trend is related to different donor-acceptor distances (DADs) at tunneling-ready states (TRSs), which could be sampled by protein dynamics. That is, a more rigid system of densely populated short DADs gives rise to a weaker T dependence of KIEs. Theoreticians have attempted to develop H-tunneling theories to explain the observations, but none have been universally accepted. It is reasonable to assume that the DAD sampling concept, if it exists, applies to the H-transfer reactions in solution, as well. In this work, we designed NADH/NAD+ model reactions to investigate their structural effects on the T dependence of hydride KIEs in acetonitrile. Hammett correlations together with N-CH3/CD3 secondary KIEs were used to provide the electronic structure of the TRSs and thus the rigidity of their charge-transfer complexation vibrations. In all three pairs of reactions, a weaker T dependence of KIEs always corresponds to a steeper Hammett slope on the substituted hydride acceptors. It was found that a tighter/rigid charge-transfer complexation system corresponds with a weaker T dependence of KIEs, consistent with the observations in enzymes.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

New classifications for quantum bioinformatics: Q-bioinformatics, QCt-bioinformatics, QCg-bioinformatics, and QCr-bioinformatics

Bioinformatics has revolutionized biology and medicine by using computational methods to analyze and interpret biological data. Quantum mechanics has recently emerged as a promising tool for the analysis of biological systems, leading to the development of quantum bioinformatics. This new field employs the principles of quantum mechanics, quantum algorithms, and quantum computing to solve complex problems in molecular biology, drug design, and protein folding. However, the intersection of bioinformatics, biology, and quantum mechanics presents unique challenges. One significant challenge is the possibility of confusion among scientists between quantum bioinformatics and quantum biology, which have similar goals and concepts. Additionally, the diverse calculations in each field make it difficult to establish boundaries and identify purely quantum effects from other factors that may affect biological processes. This review provides an overview of the concepts of quantum biology and quantum mechanics and their intersection in quantum bioinformatics. We examine the challenges and unique features of this field and propose a classification of quantum bioinformatics to promote interdisciplinary collaboration and accelerate progress. By unlocking the full potential of quantum bioinformatics, this review aims to contribute to our understanding of quantum mechanics in biological systems.

Isotope effects on radical pair performance in cryptochrome: A new hypothesis for the evolution of animal migration: The quantum biology of migration

Mechanisms occurring at the atomic level are now known to drive processes essential for life, as revealed by quantum effects on biochemical reactions. Some macroscopic characteristics of organisms may thus show an atomic imprint, which may be transferred across organisms and affect their evolution. This possibility is considered here for the first time, with the aim of elucidating the appearance of an animal innovation with an unclear evolutionary origin: migratory behaviour. This trait may be mediated by a radical pair (RP) mechanism in the retinal flavoprotein cryptochrome, providing essential magnetic orientation for migration. Isotopes may affect the performance of quantum processes through their nuclear spin. Here, we consider a simple model and then apply the standard open quantum system approach to the spin dynamics of cryptochrome RP. We changed the spin quantum number (I) and g-factor of hydrogen and nitrogen isotopes to investigate their effect on RP's yield and magnetic sensitivity. Strong differences arose between isotopes with I = 1 and I = 1/2 in their contribution to cryptochrome magnetic sensitivity, particularly regarding Earth's magnetic field strengths (25-65 µT). In most cases, isotopic substitution improved RP's magnetic sensitivity. Migratory behaviour may thus have been favoured in animals with certain isotopic compositions of cryptochrome.

Koopmans' theorem for acidic protons

The famous Brønsted acidity, which is relevant in many areas of experimental and synthetic chemistry, but also in biochemistry and other areas, is investigated from a new perspective. Nuclear electronic orbital methods, which explicitly account for the quantum character of selected protons, are applied. The resulting orbital energies of the proton wavefunction are interpreted and related to enthalpies of deprotonation and acid strength in analogy to the Koopmans' theorem for electrons. For a set of organic acids, we observe a correlation which indicates the validity of such a NEO-Koopmans' approach and opens up new opportunities for the computational investigation of more complex acidic systems.

Quantum Biology and the Potential Role of Entanglement and Tunneling in Non-Targeted Effects of Ionizing Radiation: A Review and Proposed Model

It is well established that cells, tissues, and organisms exposed to low doses of ionizing radiation can induce effects in non-irradiated neighbors (non-targeted effects or NTE), but the mechanisms remain unclear. This is especially true of the initial steps leading to the release of signaling molecules contained in exosomes. Voltage-gated ion channels, photon emissions, and calcium fluxes are all involved but the precise sequence of events is not yet known. We identified what may be a quantum entanglement type of effect and this prompted us to consider whether aspects of quantum biology such as tunneling and entanglement may underlie the initial events leading to NTE. We review the field where it may be relevant to ionizing radiation processes. These include NTE, low-dose hyper-radiosensitivity, hormesis, and the adaptive response. Finally, we present a possible quantum biological-based model for NTE.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

Heavy Water Reduces the Electronic Conductance of Protein Wires via Deuteron Interactions with Aromatic Residues

Proteins are versatile, self-assembling nanoelectronic components, but their hopping conductivity is expected to be influenced by solvent fluctuations. The role of the solvent was investigated by measuring the single molecule conductance of several proteins in both H2O and D2O. The conductance of a homologous series of protein wires decreases more rapidly with length in D2O, indicating a 6-fold decrease in carrier diffusion constant relative to the same protein in H2O. The effect was found to depend on the specific aromatic amino acid composition. A tryptophan zipper protein showed a decrease in conductance similar to that of the protein wires, whereas a phenylalanine zipper protein was insensitive to solvent changes. Tryptophan contains an indole amine, whereas the phenylalanine aromatic ring has no exchangeable protons, so the effect of heavy water on conductance is a consequence of specific D- or H-interactions with the aromatic residues.

Effect of solvent viscosity on the activation barrier of hydrogen tunneling in the lipoxygenase reaction

Hydrogen tunneling in enzyme reactions has played an important role in linking protein thermal motions to the chemical steps of catalysis. Lipoxygenases (LOXs) have served as model systems for such reactions, showcasing deep hydrogen tunneling mechanisms associated with enzymatic C-H bond cleavage from polyunsaturated fatty acids. Here, we examined the effect of solvent viscosity on the protein thermal motions associated with LOX catalysis using trehalose and glucose as viscogens. Kinetic analysis of the reaction of the paradigm plant orthologue, soybean lipoxygenase (SLO), with linoleic acid revealed no effect on the first-order rate constants, kcat, or activation energy, Ea. Further studies of SLO active site mutants displaying varying Eas, which have been used to probe catalytically relevant motions, likewise provided no evidence for viscogen-dependent motions. Kinetic analyses were extended to a representative fungal LOX from M. oryzae, MoLOX, and a human LOX, 15-LOX-2. While MoLOX behaved similarly to SLO, we show that viscogens inhibit 15-LOX-2 activity. The latter implicates viscogen sensitive, conformational motions in animal LOX reactions. The data provide insight into the role of water hydration layers in facilitating hydrogen (quantum) tunneling in LOX.

Ultrafast Proton Transfer and Contact Ion-Pair Formation in Formic Acid Clusters

The ultrafast proton transfer dynamics of homogeneous formic acid clusters (FA)n, n < 10, are investigated with femtosecond time-resolved mass spectrometry. We monitor the proton transfer pathway following Rydberg state electronic relaxation and find that successful ion pair formation increases logarithmically with cluster size. Ab initio calculations demonstrate similar excitation/relaxation behavior for each cluster, revealing a contact ion pair forms between two molecules composing the cluster before finally a formate anion (HCOO-) is dissociated by the probe pulse. The sub-ps time scale for rearrangement and proton transfer increases almost linearly with cluster size, requiring ∼67 fs per additional formic acid molecule and ranging from 213 ± 51 fs for the trimer to 667 ± 116 fs for FA9. The near-linear trends measured for both rearrangement lifetime and ion pair formation suggest that proton transfer is unlikely in the formic acid dimer but becomes prominent in small clusters.

Dynamic Domain Links Substrate Binding and Catalysis in the Factor-Inhibiting-HIF-1

The interplay between active-site chemistry and functionally relevant enzyme motions can provide useful insights into selective enzyme modulation. Modulation of the hypoxia-sensing function of factor-inhibiting-HIF-1 (FIH) enzyme is a potential therapeutic strategy in disease states such as ischemia and cancer. The hypoxia-sensing function of FIH relies in major part on the tight coupling of the first half of the catalytic mechanism which involves O2 activation and eventual succinate production to the second half which involves HIF-1α/CTAD substrate hydroxylation. In this study, we demonstrate the role of a loop hinge domain in FIH (FIH102-118) called the 100s loop in maintaining this particular tight coupling. Molecular dynamics patterns from Gaussian Network Model (iGNM) database analysis of FIH identified the 100s loop as one dynamic domain containing a hinge residue (Tyr102) with a potential substrate positioning role. Enzymological and biophysical studies of the 100s loop point mutants revealed altered enzyme kinetics with the exception of the conservative FIH mutant Y102F, which suggests a sterics-related role for this residue. Removal of the bulk of Tyr102 (Y102A) resulted in succinate production, autohydroxylation, and an O2 binding environment comparable to wild-type FIH. However, the HIF-1α/CTAD substrate hydroxylation of this mutant was significantly reduced which implies that (1) the FIH loop hinge residue Tyr102 does not affect O2 activation, (2) the stacking steric interaction of Tyr102 is important in substrate positioning for productive hydroxylation, and (3) Tyr102 is important for the synchronization of O2 activation and substrate hydroxylation.

Nonadiabatic Nuclear-Electron Dynamics: A Quantum Computing Approach

Coupled quantum electron-nuclear dynamics is often associated with the Born-Huang expansion of the molecular wave function and the appearance of nonadiabatic effects as a perturbation. On the other hand, native multicomponent representations of electrons and nuclei also exist, which do not rely on any a priori approximation. However, their implementation is hampered by prohibitive scaling. Consequently, quantum computers offer a unique opportunity for extending their use to larger systems. Here, we propose a quantum algorithm for simulating the time-evolution of molecular systems and apply it to proton transfer dynamics in malonaldehyde, described as a rigid scaffold. The proposed quantum algorithm can be easily generalized to include the explicit dynamics of the classically described molecular scaffold. We show how entanglement between electronic and nuclear degrees of freedom can persist over long times if electrons do not follow the nuclear displacement adiabatically. The proposed quantum algorithm may become a valid candidate for the study of such phenomena when sufficiently powerful quantum computers become available.

Electron Tunneling in Biology: When Does it Matter?

Electrons can tunnel between cofactor molecules positioned along biological electron transport chains up to a distance of ≃ 20 Å on the millisecond time scale of enzymatic turnover. This tunneling range determines the design of biological energy chains facilitating the cross-membrane transport of electrons. Tunneling distance and cofactors' redox potentials become the main physical parameters affecting the rate of electron transport. In addition, universal charge-transport properties are assigned to all proteins, making protein identity, flexibility, and dynamics insignificant. This paradigm is challenged by dynamical models of electron transfer, showing that the electron hopping rate is constant within the crossover distance R* ≃ 12 Å, followed with an exponential falloff at longer distances. If this hypothesis is fully confirmed, natural and man-made energy chains for electron transport should be best designed by placing redox cofactors near the crossover distance R*. Protein flexibility and dynamics affect the magnitude of the maximum hopping rate within the crossover distance. Changes in protein flexibility between forward and backward transitions contribute to vectorial charge transport. For biological energy chains, charge transport through proteins is not defined by universal parameters, and protein identity matters.

Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme's conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme's conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

Water Networks in Photosystem II Using Crystalline Molecular Dynamics Simulations and Room-Temperature XFEL Serial Crystallography

Structural dynamics of water and its hydrogen-bonding networks play an important role in enzyme function via the transport of protons, ions, and substrates. To gain insights into these mechanisms in the water oxidation reaction in Photosystem II (PS II), we have performed crystalline molecular dynamics (MD) simulations of the dark-stable S1 state. Our MD model consists of a full unit cell with 8 PS II monomers in explicit solvent (861 894 atoms), enabling us to compute the simulated crystalline electron density and to compare it directly with the experimental density from serial femtosecond X-ray crystallography under physiological temperature collected at X-ray free electron lasers (XFELs). The MD density reproduced the experimental density and water positions with high fidelity. The detailed dynamics in the simulations provided insights into the mobility of water molecules in the channels beyond what can be interpreted from experimental B-factors and electron densities alone. In particular, the simulations revealed fast, coordinated exchange of waters at sites where the density is strong, and water transport across the bottleneck region of the channels where the density is weak. By computing MD hydrogen and oxygen maps separately, we developed a novel Map-based Acceptor-Donor Identification (MADI) technique that yields information which helps to infer hydrogen-bond directionality and strength. The MADI analysis revealed a series of hydrogen-bond wires emanating from the Mn cluster through the Cl1 and O4 channels; such wires might provide pathways for proton transfer during the reaction cycle of PS II. Our simulations provide an atomistic picture of the dynamics of water and hydrogen-bonding networks in PS II, with implications for the specific role of each channel in the water oxidation reaction.

Atomistic description of the relationship between protein dynamics and catalysis with transition path sampling

Despite initial resistance, it has been increasingly accepted that protein dynamics plays a role in enzymatic catalysis. There have been two lines of research. Some works study slow conformational motions that are not coupled to the reaction coordinate, but guide the system towards catalytically competent conformations. Understanding at the atomistic level how this is accomplished has remained elusive except for a few systems. In this review we focus on fast sub-picosecond motions that are coupled to the reaction coordinate. The use of Transition Path Sampling has allowed us an atomistic description of how these rate-promoting vibrational motions are incorporated in the reaction mechanism. We will also show how we used insights from rate-promoting motions in protein design.

Phonon-assisted electron-proton transfer in [FeFe] hydrogenases: Topological role of clusters

[FeFe] hydrogenases are enzymes that have acquired a unique capacity to synthesize or consume molecular hydrogen (H2). This function relies on a complex catalytic mechanism involving the active site and two distinct electron and proton transfer networks working in concert. By an analysis based on terahertz vibrations of [FeFe] hydrogenase structure, we are able to predict and identify the existence of rate-promoting vibrations at the catalytic site and the coupling with functional residues involved in reported electron and proton transfer networks. Our findings suggest that the positioning of the cluster is influenced by the response of the scaffold to thermal fluctuations, which in turn drives the formation of networks for electron transfer through phonon-assisted mechanisms. Thus, we address the problem of linking the molecular structure to the catalytic function through picosecond dynamics, while raising the functional gain brought by the cofactors or clusters, using the concept of fold-encoded localized vibrations.

Cardiac Contractility Modulation Therapy in Patients with Amyloid Cardiomyopathy and Heart Failure, Case Report, Review of the Biophysics of CCM Function, and AMY-CCM Registry Presentation

Cardiac amyloidosis may result in an aggressive form of heart failure (HF). Cardiac contractility modulation (CCM) has been shown to be a concrete therapeutic option in patients with symptomatic HF, but there is no evidence of its application in patients with cardiac amyloidosis. We present the case of TTR amyloidosis, where CCM therapy proved to be effective. The patient had a history of multiple HF hospitalizations due to an established diagnosis of wild type TTR-Amyloidosis with significant cardiac involvement. Since he was highly symptomatic, except during continuous dobutamine and diuretic infusion, it was opted to pursue CCM therapy device implantation. At follow up, a significant improvement in clinical status was reported with an increase of EF, functional status (6 min walk test improved from zero meters at baseline, to 270 m at 1 month and to 460 m at 12 months), and a reduction in pulmonary pressures. One year after device implantation, no other HF hospital admission was needed. CCM therapy may be effective in this difficult clinical setting. The AMY-CCM Registry, which has just begun, will evaluate the efficacy of CCM in patients with HF and diagnosed TTR amyloidosis to bring new evidence on its potential impact as a therapeutic option.

Effect of Water Deuteration on Protein Electron Transfer

Traditional theories of long-range protein electron transfer describe the reaction rate in terms of the tunneling distance and the reaction free energy. They do not recognize two physical effects: (i) local wetting of the active site by hydration water and (ii) protein identity affecting the rate through dynamics and flexibility. We find, by molecular dynamics simulations, a significant, ∼25 times, slowing down of the rate of protein electron transfer upon deuteration. H/D substitution changes the rate constant pre-exponential factor in the regime of electron transfer controlled by medium dynamics. Switching from light to heavy water increases the effective medium relaxation time. The effect is caused by both a global change in the flexibility of the protein backbone and locally stronger hydrogen bonds to charged residues.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Solvent Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects: Correlation to the Donor-Acceptor Distances

Protein structural effects on the temperature (T) dependence of kinetic isotope effects (KIEs) in H-tunneling reactions have recently been used to discuss about the role of enzyme thermal motions in catalysis. Frequently observed nearly T-independent KIEs in the wild-type enzymes and T-dependent KIEs in variants suggest that H-tunneling in the former is assisted by the naturally evolved protein constructive vibrations that help sample short donor-acceptor distances (DADs) needed. This explanation that correlates the T-dependence of KIEs with DAD sampling has been highly debated as simulations following other H-tunneling models sometimes gave alternative explanations. In this paper, solvent effects on the T-dependence of KIEs of two hydride tunneling reactions of NADH/NAD⁺ analogues (represented by ΔE_a = E_aD - E_aH) were determined in attempts to replicate the observations in enzymes and test the protein vibration-assisted DAD sampling concept. Effects of selected aprotic solvents on the DAD_PRC's of the productive reactant complexes (PRCs) and the DAD_TRS's of the activated tunneling ready states (TRSs) were obtained through computations and analyses of the kinetic data, including 2° KIEs, respectively. A weaker T-dependence of KIEs (i.e., smaller ΔE_a) was found in a more polar aprotic solvent in which the system has a shorter average DAD_PRC and DAD_TRS. Further results show that a charge-transfer (CT) complexation made of a stronger donor/acceptor gives rise to a smaller ΔE_a. Overall, the shorter and less broadly distributed DADs resulting from the stronger CT complexation vibrations give rise to a smaller ΔE_a. Our results appear to support the explanation that links the T-dependence of KIEs to the donor-acceptor rigidity in enzymes.

Dependence of crystallographic atomic displacement parameters on temperature (25-150 K) for complexes of horse liver alcohol dehydrogenase

Enzymes catalyze reactions by binding and orienting substrates with dynamic interactions. Horse liver alcohol dehydrogenase catalyzes hydrogen transfer with quantum-mechanical tunneling that involves fast motions in the active site. The structures and B factors of ternary complexes of the enzyme with NAD+ and 2,3,4,5,6-pentafluorobenzyl alcohol or NAD+ and 2,2,2-trifluoroethanol were determined to 1.1-1.3 Å resolution below the `glassy transition' in order to extract information about the temperature-dependent harmonic motions, which are reflected in the crystallographic B factors. The refinement statistics and structures are essentially the same for each structure at all temperatures. The B factors were corrected for a small amount of radiation decay. The overall B factors for the complexes are similar (13-16 Å2) over the range 25-100 K, but increase somewhat at 150 K. Applying TLS refinement to remove the contribution of pseudo-rigid-body displacements of coenzyme binding and catalytic domains provided residual B factors of 7-10 Å2 for the overall complexes and of 5-10 Å2 for C4N of NAD+ and the methylene carbon of the alcohols. These residual B factors have a very small dependence on temperature and include local harmonic motions and apparently contributions from other sources. Structures at 100 K show complexes that are poised for hydrogen transfer, which involves atomic displacements of ∼0.3 Å and is compatible with the motions estimated from the residual B factors and molecular-dynamics simulations. At 298 K local conformational changes are also involved in catalysis, as enzymes with substitutions of amino acids in the substrate-binding site have similar positions of NAD+ and pentafluorobenzyl alcohol and similar residual B factors, but differ by tenfold in the rate constants for hydride transfer.

Measurement of kinetic isotope effects on peptide hydroxylation using MALDI-MS

Primary kinetic isotope effects (KIEs) provide unique insight into enzymatic reactions, as they can reveal rate-limiting steps and detailed chemical mechanisms. HIF hydroxylases, part of a family of 2-oxoglutarate (2OG) oxygenases are central to the regulation of many crucial biological processes through O2-sensing, but present a challenge to monitor due to the large size of the protein substrate and the similarity between native and hydroxylated substrate. MALDI-TOF MS is a convenient tool to measure peptide masses, which can also be used to measure the discontinuous kinetics of peptide hydroxylation for Factor Inhibiting HIF (FIH). Using this technique, rate data can be observed from the mole-fraction of CTAD and CTAD-OH in small volumes, allowing noncompetitive H/D KIEs to be measured. Slow dCTAD substrate leads to extensive uncoupling of O2 consumption from peptide hydroxylation, leading to enzyme autohydroxylation, which is observed using UV-vis spectroscopy. Simultaneously measuring both the normal product, CTAD-OH, and the uncoupled product, autohydroxylated enzyme, the KIE on the microscopic step of hydrogen atom transfer (HAT) can be estimated. MALDI-MS analysis is a strong method for monitoring reactions that hydroxylate peptides, and can be generalized to other similar reactions, and simultaneous kinetic detection of branched products can provide valuable insight on microscopic KIEs at intermediate mechanistic steps.

Ammonia quantum tunneling in cold rare-gas He and Ar clusters and factorial design approach for methodology evaluation

Quantum tunneling of the ammonia inversion motion and energy level splittings in He and Ar clusters were investigated. It was found that the double well potential (DWP) in He clusters is symmetrical and that the first layer of He atoms is able to model the system. The calculated tunneling splitting was in good agreement with the experimental, 36.4 and 24.6 cm[Formula: see text] respectively. For NH[Formula: see text] in Ar clusters, the DWP becomes slightly asymmetric, which is enough to decrease the resonance and make the symmetric DWP unable to model the system. An asymmetric potential was used and the result was in excellent agreement with the experimental splitting, of 9.0 and 10.6 cm[Formula: see text] respectively. Non-covalent interactions revealed that the asymmetry is caused by dissimilar interactions in each minimum of the double well potential. The effects of different methodologies were analyzed via a design of experiments approach. For the gas-phase NH[Formula: see text] molecule, only diffuse functions were statistically significant while for the NH[Formula: see text] embedded in He cluster both the MP2 method and polarization functions were significant. This tendency suggests higher order polarization functions may be essential to generate accurate barrier heights.

Conformational dynamics modulating electron transfer

Diffusional dynamics of the donor-acceptor distance are responsible for the appearance of a new time scale of diffusion over the distance of electronic tunneling in electron-transfer reactions. The distance dynamics compete with the medium polarization dynamics in the dynamics-controlled electron-transfer kinetics. The pre-exponential factor of the electron-transfer rate constant switches, at the crossover distance, between a distance-independent, dynamics-controlled plateau and exponential distance decay. The crossover between two regimes is controlled by an effective relaxation time slowed down by a factor exponentially depending on the variance of the donor-acceptor displacement. Flexible donor-acceptor complexes must show a greater tendency for dynamics-controlled electron transfer. Energy chains based on electron transport are best designed by placing the redox cofactors near the crossover distance.

Temperature-dependent hydrogen deuterium exchange shows impact of analog binding on adenosine deaminase flexibility but not embedded thermal networks

The analysis of hydrogen deuterium exchange by mass spectrometry as a function of temperature and mutation (TDHDX-MS) has emerged as a generic and efficient tool for the spatial resolution of protein networks that are proposed to function in the thermal activation of catalysis. In this work, we extend TDHDX from apo-enzyme structures to protein-ligand complexes. Using adenosine deaminase as a prototype, we compared the impacts of a substrate analog (1-deaza-adenosine or DAA) and a very tight-binding inhibitor/transition state analog (pentostatin) at single and multiple temperatures. At a single temperature, we observed different HDX-MS properties for the two ligands, as expected from their 10⁶-fold differences in strength of binding. By contrast, analogous patterns for TDHDX-MS emerge in the presence of both DAA and pentostatin, indicating similar impacts of either ligand on the enthalpic barriers for local protein unfolding. We extended TDHDX to a function-altering mutant of adenosine deaminase in the presence of pentostatin and revealed a protein thermal network that is highly similar to that previously reported for the apo-enzyme (Gao et al., 2020, JACS 142, 19936-19949). Finally, we discuss the differential impacts of pentostatin binding on overall protein flexibility vs. site-specific thermal transfer pathways in the context of models for substrate-induced changes to a distributed protein conformational landscape that act in synergy with embedded protein thermal networks to achieve efficient catalysis.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

Controlling Tunneling Oscillation and Quantum Localization in an Asymmetric Double-Well Potential: A Bohmian Perspective

The roles of spatial symmetry and strength of external time-dependent perturbation on the dynamics of a quantum particle, initially localized in one of the wells of an asymmetric double-well potential are studied using the recently developed techniques incorporating quantum theory of motion and time-dependent Fourier grid Hamiltonian methods. The model used here includes a mimic of the related experimental situations which is considered as a perturbation to the static double-well potential. Analysis of localized and delocalized phase space structures and corresponding time-profile of tunneling probability reveal the recipe toward controlling the tunneling oscillations by modulating the parameters of applied perturbation. A study on a stochastic pulsating potential also reveals the root to the quantum localization, even in moderate field strength.

Chronobiology Meets Quantum Biology: A New Paradigm Overlooking the Horizon?

Biological processes and physiological functions in living beings are featured by oscillations with a period of about 24 h (circadian) or cycle at the second and third harmonic (ultradian) of the basic frequency, driven by the biological clock. This molecular mechanism, common to all kingdoms of life, comprising animals, plants, fungi, bacteria, and protists, represents an undoubted adaptive advantage allowing anticipation of predictable changes in the environmental niche or of the interior milieu. Biological rhythms are the field of study of Chronobiology. In the last decade, growing evidence hints that molecular platforms holding up non-trivial quantum phenomena, including entanglement, coherence, superposition and tunnelling, bona fide evolved in biosystems. Quantum effects have been mainly implicated in processes related to electromagnetic radiation in the spectrum of visible light and ultraviolet rays, such as photosynthesis, photoreception, magnetoreception, DNA mutation, and not light related such as mitochondrial respiration and enzymatic activity. Quantum effects in biological systems are the field of study of Quantum Biology. Rhythmic changes at the level of gene expression, as well as protein quantity and subcellular distribution, confer temporal features to the molecular platform hosting electrochemical processes and non-trivial quantum phenomena. Precisely, a huge amount of molecules plying scaffold to quantum effects show rhythmic level fluctuations and this biophysical model implies that timescales of biomolecular dynamics could impinge on quantum mechanics biofunctional role. The study of quantum phenomena in biological cycles proposes a profitable “entanglement” between the areas of interest of these seemingly distant scientific disciplines to enlighten functional roles for quantum effects in rhythmic biosystems.

Multidimensional insights into the repeated electromagnetic field stimulation and biosystems interaction in aging and age-related diseases

We provide a multidimensional sequence of events that describe the electromagnetic field (EMF) stimulation and biological system interaction. We describe this process from the quantum to the molecular, cellular, and organismal levels. We hypothesized that the sequence of events of these interactions starts with the oscillatory effect of the repeated electromagnetic stimulation (REMFS). These oscillations affect the interfacial water of an RNA causing changes at the quantum and molecular levels that release protons by quantum tunneling. Then protonation of RNA produces conformational changes that allow it to bind and activate Heat Shock Transcription Factor 1 (HSF1). Activated HSF1 binds to the DNA expressing chaperones that help regulate autophagy and degradation of abnormal proteins. This action helps to prevent and treat diseases such as Alzheimer’s and Parkinson’s disease (PD) by increasing clearance of pathologic proteins. This framework is based on multiple mathematical models, computer simulations, biophysical experiments, and cellular and animal studies. Results of the literature review and our research point towards the capacity of REMFS to manipulate various networks altered in aging (Reale et al. PloS one 9, e104973, 2014), including delay of cellular senescence (Perez et al. 2008, Exp Gerontol 43, 307-316) and reduction in levels of amyloid-β peptides (Aβ) (Perez et al. 2021, Sci Rep 11, 621). Results of these experiments using REMFS at low frequencies can be applied to the treatment of patients with age-related diseases. The use of EMF as a non-invasive therapeutic modality for Alzheimer’s disease, specifically, holds promise. It is also necessary to consider the complicated and interconnected genetic and epigenetic effects of the REMFS-biological system’s interaction while avoiding any possible adverse effects.

Leveraging intrinsic flexibility to engineer enhanced enzyme catalytic activity

Dynamic motions of enzymes occurring on a broad range of timescales play a pivotal role in all steps of the reaction pathway, including substrate binding, catalysis, and product release. However, it is unknown whether structural information related to conformational flexibility can be exploited for the directed evolution of enzymes with higher catalytic activity. Here, we show that mutagenesis of residues exclusively located at flexible regions distal to the active site of Homo sapiens kynureninase (HsKYNase) resulted in the isolation of a variant (BF-HsKYNase) in which the rate of the chemical step toward kynurenine was increased by 45-fold. Mechanistic pre–steady-state kinetic analysis of the wild type and the evolved enzyme shed light on the underlying effects of distal mutations (>10 Å from the active site) on the rate-limiting step of the catalytic cycle. Hydrogen-deuterium exchange coupled to mass spectrometry and molecular dynamics simulations revealed that the amino acid substitutions in BF-HsKYNase allosterically affect the flexibility of the pyridoxal-5′-phosphate (PLP) binding pocket, thereby impacting the rate of chemistry, presumably by altering the conformational ensemble and sampling states more favorable to the catalyzed reaction.

Hydride Transfer Mechanism of Enzymatic Sugar Nucleotide C2 Epimerization Probed with a Loose-Fit CDP-Glucose Substrate

Transient oxidation–reduction through hydride transfer with tightly bound NAD coenzyme is used by a large class of sugar nucleotide epimerases to promote configurational inversion of carbon stereocenters in carbohydrate substrates. A requirement for the epimerases to coordinate hydride abstraction and re-addition with substrate rotation in the binding pocket poses a challenge for dynamical protein conformational selection linked to enzyme catalysis. Here, we studied the thermophilic C2 epimerase from Thermodesulfatator atlanticus (TaCPa2E) in combination with a slow CDP-glucose substrate (k_cat ≈ 1.0 min^–1; 60 °C) to explore the sensitivity of the enzymatic hydride transfer toward environmental fluctuations affected by temperature (20–80 °C). We determined noncompetitive primary kinetic isotope effects (KIE) due to ²H at the glucose C2 and showed that a normal KIE on the k_cat (^Dk_cat) reflects isotope sensitivity of the hydrogen abstraction to enzyme-NAD⁺ in a rate-limiting transient oxidation. The ^Dk_cat peaked at 40 °C was 6.1 and decreased to 2.1 at low (20 °C) and 3.3 at high temperature (80 °C). The temperature profiles for k_cat with the ¹H and ²H substrate showed a decrease in the rate below a dynamically important breakpoint (∼40 °C), suggesting an equilibrium shift to an impaired conformational landscape relevant for catalysis in the low-temperature region. Full Marcus-like model fits of the rate and KIE profiles provided evidence for a high-temperature reaction via low-frequency conformational sampling associated with a broad distribution of hydride donor–acceptor distances (long-distance population centered at 3.31 ± 0.02 Å), only poorly suitable for quantum mechanical tunneling. Collectively, dynamical characteristics of TaCPa2E-catalyzed hydride transfer during transient oxidation of CDP-glucose reveal important analogies to mechanistically simpler enzymes such as alcohol dehydrogenase and dihydrofolate reductase. A loose-fit substrate (in TaCPa2E) resembles structural variants of these enzymes by extensive dynamical sampling to balance conformational flexibility and catalytic efficiency.

Modeling the Ligand Effect on the Structure of CYP 450 Within the Density Functional Theory

An improved understanding of the P450 structure is relevant to the development of biomimetic catalysts and inhibitors for controlled CH-bond activation, an outstanding challenge of synthetic chemistry. Motivated by the experimental findings of an unusually short Fe-S bond of 2.18 Å for the wild-type (WT) OleT P450 decarboxylase relative to a cysteine pocket mutant form (A369P), a computational model that captures the effect of the thiolate axial ligand on the iron-sulfur distance is presented. With the computational efficiency and streamlined analysis in mind, this model combines a cluster representation of the enzyme─40-110 atoms, depending on the heme and ligand truncation level─with a density functional theory (DFT) description of the electronic structure (ES) and is calibrated against the experimental data. The optimized Fe-S distances show a difference of 0.25 Å between the low and high spin states, in agreement with the crystallographic structures of the OleT WT and mutant forms. We speculate that this difference is attributable to the packing of the ligand; the mutant is bulkier due to an alanine-to-proline replacement, meaning that it is excluded from the energetically favored low-spin minimum because of steric constraints. The presence of pure spin-state pairs and the intersection of the low/high spin states for the enzyme model is indicative of the limitations of single-reference ES methods in such systems and emphasizes the significance of using the proper state when modeling the hydrogen atom transfer (HAT) reaction catalyzed by OleT. At the same time, the correct characterization of both the short and long Fe-S bonds within a small DFT-based model of 42 atoms paves the way for quantum dynamics modeling of the HAT step, which initiates the OleT decarboxylation reaction.

Experimental and Theoretical Examination of the Kinetic Isotope Effect in Cytochrome P450 Decarboxylase OleT

Using a combination of experimental studies, theory, simulation, and modeling, we investigate the hydrogen atom transfer (HAT) reaction by the high-valent ferryl cytochrome P450 (CYP) intermediate known as Compound I, a species that is central to innumerable and important detoxification and biosynthetic reactions. The P450 decarboxylase known as OleT converts fatty acids, a sustainable biological feedstock, into terminal alkenes and thus is of high interest as a potential means to produce fungible biofuels. Previous experimental work has established the intermediacy of Compound I in the C─C scission reaction catalyzed by OleT and an unprecedented ability to monitor the HAT process in the presence of bound fatty acid substrates. Here, we leverage the kinetic simplicity of the OleT system to measure the activation barriers for CYP HAT and the temperature dependence of the substrate ²H kinetic isotope effect. Notably, neither measurement has been previously accessible for a CYP to date. Theoretical analysis alludes to the significance of substrate fatty acid coordination for generating the hydrogen donor/acceptor configurations that are most conducive for HAT to occur. The analysis of the two-dimensional potential energy surface, based on multireference electronic wave functions, illustrates the uncoupled character of the hydrogen motion. Quantum dynamics calculations along the hydrogen reaction path demonstrate that hydrogen tunneling is essential to qualitatively capture the experimental isotope effect, its temperature dependence, and appropriate activation energies. Overall, a more fundamental understanding of the OleT reaction coordinate contributes to the development of biomimetic catalysts for controlled C─H bond activation, an outstanding current challenge for (bio)synthetic chemistry.

L-tyrosine-bound ThiH structure reveals C-C bond break differences within radical SAM aromatic amino acid lyases

2-iminoacetate synthase ThiH is a radical S-adenosyl-L-methionine (SAM) L-tyrosine lyase and catalyzes the L-tyrosine Cα-Cβ bond break to produce dehydroglycine and p-cresol while the radical SAM L-tryptophan lyase NosL cleaves the L-tryptophan Cα-C bond to produce 3-methylindole-2-carboxylic acid. It has been difficult to understand the features that condition one C-C bond break over the other one because the two enzymes display significant primary structure similarities and presumably similar substrate-binding modes. Here, we report the crystal structure of L-tyrosine bound ThiH from Thermosinus carboxydivorans revealing an unusual protonation state of L-tyrosine upon binding. Structural comparison of ThiH with NosL and computational studies of the respective reactions they catalyze show that substrate activation is eased by tunneling effect and that subtle structural changes between the two enzymes affect, in particular, the hydrogen-atom abstraction by the 5´-deoxyadenosyl radical species, driving the difference in reaction specificity.

IR-induced and tunneling reactions in cryogenic matrices: the (incomplete) story of a successful endeavor

In this article, IR-induced and tunneling-driven reactions observed in cryogenic matrices are described in a historical perspective, the entangling of the two types of processes being highlighted. The story of this still ongoing fascinating scientific endeavor is presented here following closely our own involvement in the field for more than 30 years, and thus focuses mostly on our work. It is, because of this reason, also an incomplete story. Nevertheless, it considers a large range of examples, from very selective IR-induced conformational isomerizations to IR-induced bond-breaking/bond-forming reactions and successful observations of rare heavy atom tunneling processes. As a whole, this article provides a rather general overview of the major progress achieved in the field.

Structural and Thermodynamic Effects on the Kinetics of C-H Oxidation by Multisite Proton-Coupled Electron Transfer in Fluorenyl Benzoates

Our recent experimental and theoretical investigations have shown that fluorene C-H bonds can be activated through a mechanism in which the proton and electron are transferred from the C-H bond to a separate base and oxidant in a concerted, elementary step. This multisite proton-coupled electron transfer (MS-PCET) mechanism for C-H bond activation was shown to be facilitated by shorter proton donor-acceptor distances. With the goal of intentionally modulating this donor-acceptor distance, we have now studied C-H MS-PCET in the 3-methyl-substituted fluorenyl benzoate (2-Flr-3-Me-BzO^-). This derivative was readily oxidized by ferrocenium oxidants by initial C-H MS-PCET, with rate constants that were 6-21 times larger than those for 2-Flr-BzO^- with the same oxidants. Structural comparisons by X-ray crystallography and by computations showed that addition of the 3-methyl group caused the expected steric compression; however, the relevant C···O^- proton donor-acceptor distance was longer, due to a twist of the carboxylate group. The structural changes induced by the 3-Me group increased the basicity of the carboxylate, weakened the C-H bond, and reduced the reorganization energy for C-H bond cleavage. Thus, the rate enhancement for 2-Flr-3-Me-BzO^- was due to effects on the thermochemistry and kinetic barrier, rather than from compression of the C···O^- proton donor-acceptor distance. These results highlight both the challenges of controlling molecules on the 0.1 Å length scale and the variety of parameters that affect PCET rate constants.

Ferritin-Like Proteins: A Conserved Core for a Myriad of Enzyme Complexes

Ferritin-like proteins share a common fold, a four α-helix bundle core, often coordinating a pair of metal ions. Although conserved, the ferritin fold permits a diverse set of reactions, and is central in a multitude of macromolecular enzyme complexes. Here, we emphasize this diversity through three members of the ferritin-like superfamily: the soluble methane monooxygenase, the class I ribonucleotide reductase and the aldehyde deformylating oxygenase. They all rely on dinuclear metal cofactors to catalyze different challenging oxygen-dependent reactions through the formation of multi-protein complexes. Recent studies using cryo-electron microscopy, serial femtosecond crystallography at an X-ray free electron laser source, or single-crystal X-ray diffraction, have reported the structures of the active protein complexes, and revealed unprecedented insights into the molecular mechanisms of these three enzymes.

Condensed-phase isomerization through tunnelling gateways

Quantum mechanical tunnelling describes transmission of matter waves through a barrier with height larger than the energy of the wave¹. Tunnelling becomes important when the de Broglie wavelength of the particle exceeds the barrier thickness; because wavelength increases with decreasing mass, lighter particles tunnel more efficiently than heavier ones. However, there exist examples in condensed-phase chemistry where increasing mass leads to increased tunnelling rates². In contrast to the textbook approach, which considers transitions between continuum states, condensed-phase reactions involve transitions between bound states of reactants and products. Here this conceptual distinction is highlighted by experimental measurements of isotopologue-specific tunnelling rates for CO rotational isomerization at an NaCl surface^3,4, showing nonmonotonic mass dependence. A quantum rate theory of isomerization is developed wherein transitions between sub-barrier reactant and product states occur through interaction with the environment. Tunnelling is fastest for specific pairs of states (gateways), the quantum mechanical details of which lead to enhanced cross-barrier coupling; the energies of these gateways arise nonsystematically, giving an erratic mass dependence. Gateways also accelerate ground-state isomerization, acting as leaky holes through the reaction barrier. This simple model provides a way to account for tunnelling in condensed-phase chemistry, and indicates that heavy-atom tunnelling may be more important than typically assumed.

Experimental mapping of short-wavelength phonons in proteins

Phonons are quasi-particles, observed as lattice vibrations in periodic materials, that often dampen in the presence of structural perturbations. Nevertheless, phonon-like collective excitations exist in highly complex systems, such as proteins, although the origin of such collective motions has remained elusive. Here we present a picture of temperature and hydration dependence of collective excitations in green fluorescent protein (GFP) obtained by inelastic neutron scattering. Our results provide evidence that such excitations can be used as a measure of flexibility/softness and are possibly associated with the protein’s activity. Moreover, we show that the hydration water in GFP interferes with the phonon propagation pathway, enhancing the structural rigidity and stability of GFP.

•

Quantum phenomena in biology have long fascinated people around the world

•

This work presents a direct experimental observation of phonons, the quantum vibrations in a protein

•

The collective excitations or phonons in proteins were detected by utilizing inelastic neutron scattering technique at Oak Ridge National Laboratory

•

Our results illustrate the flexibility-activity relationship in proteins by mapping the temperature and hydration dependence of these collective excitations

The role of conformational dynamics on the activity of polymer-conjugated CalB in organic solvents

A perennial interest in enzyme catalysis has been expanding its applicability from aqueous phase where enzymes are naturally evolved to organic solvents in which the majority of industrial chemical synthesis...

Evolution of Optimized Hydride Transfer Reaction and Overall Enzyme Turnover in Human Dihydrofolate Reductase

Evolution of dihydrofolate reductase (DHFR) has been studied using the enzyme from Escherichia coli DHFR (ecDHFR) as a model, but less studies have used the enzyme from Homo sapiens DHFR (hsDHFR). Each enzyme maintains a short and narrow distribution of hydride donor-acceptor distances (DAD) at the tunneling ready state (TRS). Evolution of the enzyme was previously studied in ecDHFR where three key sites were identified as important to the catalyzed reaction. The corresponding sites in hsDHFR are F28, 62-PEKN, and 26-PPLR. Each of these sites was studied here through the creation of mutant variants of the enzyme and measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs) on the reaction. F28 is mutated first to M (F28M) and then to the L of the bacterial enzyme (F28L). The KIEs of the F28M variant are larger and more temperature-dependent than wild-type (WT), suggesting a broader and longer average DAD at the TRS. To more fully mimic ecDHFR, we also study a triple mutant of the human enzyme (F32L-PP26N-PEKN62G). Remarkably, the intrinsic KIEs, while larger in magnitude, are temperature-independent like the WT enzymes. We also construct deletion mutations of hsDHFR removing both the 62-PEKN and 26-PPLR sequences. The results mirror those described previously for insertion mutants of ecDHFR. Taken together, these results suggest a balancing act during DHFR evolution between achieving an optimal TRS for hydride transfer and preventing product inhibition arising from the different intercellular pools of NADPH and NADP⁺ in prokaryotic and eukaryotic cells.

Free Energies of Proton-Coupled Electron Transfer Reagents and Their Applications

We present an update and revision to our 2010 review on the topic of proton-coupled electron transfer (PCET) reagent thermochemistry. Over the past decade, the data and thermochemical formalisms presented in that review have been of value to multiple fields. Concurrently, there have been advances in the thermochemical cycles and experimental methods used to measure these values. This Review (i) summarizes those advancements, (ii) corrects systematic errors in our prior review that shifted many of the absolute values in the tabulated data, (iii) provides updated tables of thermochemical values, and (iv) discusses new conclusions and opportunities from the assembled data and associated techniques. We advocate for updated thermochemical cycles that provide greater clarity and reduce experimental barriers to the calculation and measurement of Gibbs free energies for the conversion of X to XH_n in PCET reactions. In particular, we demonstrate the utility and generality of reporting potentials of hydrogenation, E°(V vs H₂), in almost any solvent and how these values are connected to more widely reported bond dissociation free energies (BDFEs). The tabulated data demonstrate that E°(V vs H₂) and BDFEs are generally insensitive to the nature of the solvent and, in some cases, even to the phase (gas versus solution). This Review also presents introductions to several emerging fields in PCET thermochemistry to give readers windows into the diversity of research being performed. Some of the next frontiers in this rapidly growing field are coordination-induced bond weakening, PCET in novel solvent environments, and reactions at material interfaces.

Direct Observation of Ultrafast Proton Rocking in the BLUF Domain

We present direct observation of ultrafast proton rocking in the central motif of a BLUF domain protein scaffold. The mutant design has taken considerations of modulating the proton-coupled electron transfer (PCET) driving forces by replacing Tyr in the original motif with Trp, as well as of removing the interference of a competing electron transfer pathway. Using femtosecond pump-probe spectroscopy and detailed kinetics analysis, we resolved an electron-transfer-coupled Grotthuss-type forward and reversed proton rocking along the FMN-Gln-Trp proton relay chain. The rates of forward and reversed proton transfer are determined to be very close, namely 51 ps vs 52 ps. The kinetic isotope effect (KIE) constants associated with the forward and reversed proton transfer are 3.9 and 5.3, respectively. The observation of ultrafast proton rocking is not only a crucial step towards revealing the nature of proton relay in BLUF domain, but also provides a new paradigm of proton transfer in proteins for theoretical investigations.

The Proton in Biochemistry: Impacts on Bioenergetics, Biophysical Chemistry, and Bioorganic Chemistry

The proton is the smallest atomic particle, and in aqueous solution it is the smallest hydrated ion, having only two waters in its first hydration shell. In this article we survey key aspects of the proton in chemistry and biochemistry, starting with the definitions of pH and pK _a and their application inside biological cells. This includes an exploration of pH in nanoscale spaces, distinguishing between bulk and interfacial phases. We survey the Eigen and Zundel models of the structure of the hydrated proton, and how these can be used to explain: a) the behavior of protons at the water-hydrophobic interface, and b) the extraordinarily high mobility of protons in bulk water via Grotthuss hopping, and inside proteins via proton wires. Lastly, we survey key aspects of the effect of proton concentration and proton transfer on biochemical reactions including ligand binding and enzyme catalysis, as well as pH effects on biochemical thermodynamics, including the Chemiosmotic Theory. We find, for example, that the spontaneity of ATP hydrolysis at pH ≥ 7 is not due to any inherent property of ATP (or ADP or phosphate), but rather to the low concentration of H⁺. Additionally, we show that acidification due to fermentation does not derive from the organic acid waste products, but rather from the proton produced by ATP hydrolysis.