Active Site Dynamical Effects in the Hydrogen Transfer Rate-limiting Step in the Catalysis of Linoleic Acid by Soybean Lipoxygenase-1 (SLO-1): Primary and Secondary Isotope Contributions

Rare Events Sampling Methods for Quantum and Classical Ab Initio Molecular Dynamics

We provide an approach to sample rare events during classical ab initio molecular dynamics and quantum wavepacket dynamics. For classical AIMD, a set of fictitious degrees of freedom are introduced that may harmonically interact with the electronic and nuclear degrees of freedom to steer the dynamics in a conservative fashion toward energetically forbidden regions. A similar approach when introduced for quantum wavepacket dynamics has the effect of biasing the trajectory of the wavepacket centroid toward the regions of the potential surface that are difficult to sample. The approach is demonstrated for a phenol-amine system, which is a prototypical problem for condensed phase-proton transfer, and for model potentials undergoing wavepacket dynamics. In all cases, the approach yields trajectories that conserve energy while sampling rare events.

Weighted-Graph-Theoretic Methods for Many-Body Corrections within ONIOM: Smooth AIMD and the Role of High-Order Many-Body Terms

We present a weighted-graph-theoretic approach to adaptively compute contributions from many-body approximations for smooth and accurate post-Hartree-Fock (pHF) ab initio molecular dynamics (AIMD) of highly fluxional chemical systems. This approach is ONIOM-like, where the full system is treated at a computationally feasible quality of treatment (density functional theory (DFT) for the size of systems considered in this publication), which is then improved through a perturbative correction that captures local many-body interactions up to a certain order within a higher level of theory (post-Hartree-Fock in this publication) described through graph-theoretic techniques. Due to the fluxional and dynamical nature of the systems studied here, these graphical representations evolve during dynamics. As a result, energetic "hops" appear as the graphical representation deforms with the evolution of the chemical and physical properties of the system. In this paper, we introduce dynamically weighted, linear combinations of graphs, where the transition between graphical representations is smoothly achieved by considering a range of neighboring graphical representations at a given instant during dynamics. We compare these trajectories with those obtained from a set of trajectories where the range of local many-body interactions considered is increased, sometimes to the maximum available limit, which yields conservative trajectories as the order of interactions is increased. The weighted-graph approach presents improved dynamics trajectories while only using lower-order many-body interaction terms. The methods are compared by computing dynamical properties through time-correlation functions and structural distribution functions. In all cases, the weighted-graph approach provides accurate results at a lower cost.

Challenges in constructing accurate methods for hydrogen transfer reactions in large biological assemblies: rare events sampling for mechanistic discovery and tensor networks for quantum nuclear effects

We present two methods that address the computational complexities arising in hydrogen transfer reactions in enzyme active sites.

Comparative kinetic isotope effects on first- and second-order rate constants of soybean lipoxygenase variants uncover a substrate-binding network

Lipoxygenases are widespread enzymes found in virtually all eukaryotes, including fungi, and, more recently, in prokaryotes. These enzymes act on long-chain polyunsaturated fatty acid substrates (C18 to C20), raising questions regarding how the substrate threads its way from solvent to the active site. Herein, we report a comparison of the temperature dependence of isotope effects on first- and second-order rate constants among single-site variants of the prototypic plant enzyme soybean lipoxygenase-1 substituted at amino acid residues inferred to impact substrate binding. We created 10 protein variants including four amino acid positions, Val-750, Ile-552, Ile-839, and Trp-500, located within a previously proposed substrate portal. The conversion of these bulky hydrophobic side chains to smaller side chains is concluded to increase the mobility of flanking helices, giving rise to increased off rates for substrate dissociation from the enzyme. In this manner, we identified a specific "binding network" that can regulate movement of the substrate from the solvent to the active site. Taken together with our previous findings on C-H and O2 activation of soybean lipoxygenase-1, these results support the emergence of multiple complementary networks within a single protein scaffold that modulate different steps along the enzymatic reaction coordinate.

Photoelectrons Are Not Always Quite Free

The photoelectron spectra of Sm₂O^- obtained over a range of photon energies exhibit anomalous changes in relative excited-state band intensities. Specifically, the excited-state transition intensities increase relative to the transition to the neutral ground state with decreasing photon energy, the opposite of what is expected from threshold effects. This phenomenon was previously observed in studies on several Sm-rich homo- and heterolanthanide oxides collected with two different harmonic outputs of a Nd:YAG (2.330 and 3.495 eV) [ J. Chem. Phys. 2017, 146, 194310]. We relate these anomalous intensities to populations of ground and excited anionic and neutrals states through the inspection of time-dependent perturbation theory within the adiabatic and sudden limits and for the first time show that transition intensities in photoelectron spectroscopy have a deep significance in gauging participation from excited states. We believe our results will have significance in the study of other electron-rich systems that have especially high density of accessible spin states.

Impact of Mutations on the Binding Pocket of Soybean Lipoxygenase: Implications for Proton-Coupled Electron Transfer

Soybean lipoxygenase catalyzes a proton-coupled electron transfer (PCET) reaction and serves as a prototype for hydrogen tunneling in enzymes due to the unusually high kinetic isotope effect and significant modulation of the rate constant and kinetic isotope effect by mutation. Herein these experimental observations are interpreted in the context of changes to the substrate binding pocket in microsecond molecular dynamics simulations of wild-type and mutant soybean lipoxygenase. The binding pocket exhibits an hourglass shape with residues L546 and L754 bracketing the bottleneck, positioning the linoleic acid substrate for PCET. Mutation of I553 to less bulky residues slightly increases the width of the bottleneck and the volume of the binding pocket. Mutating L546 or L754 to a smaller residue also enlarges this width and volume, and mutating both has an even more dramatic effect. This analysis illustrates how mutation of the substrate binding pocket can be used as a strategy to tune the kinetics.

Fundamental Insights into Proton-Coupled Electron Transfer in Soybean Lipoxygenase from Quantum Mechanical/Molecular Mechanical Free Energy Simulations

The proton-coupled electron transfer (PCET) reaction catalyzed by soybean lipoxygenase has served as a prototype for understanding hydrogen tunneling in enzymes. Herein this PCET reaction is studied with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The free energy surfaces are computed as functions of the proton donor-acceptor (C-O) distance and the proton coordinate, and the potential of mean force is computed as a function of the C-O distance, inherently including anharmonicity. The simulation results are used to calculate the kinetic isotope effects for the wild-type enzyme (WT) and the L546A/L754A double mutant (DM), which have been measured experimentally to be ∼80 and ∼700, respectively. The PCET reaction is found to be exoergic for WT and slightly endoergic for the DM, and the equilibrium C-O distance for the reactant is found to be ∼0.2 Å greater for the DM than for WT. The larger equilibrium distance for the DM, which is due mainly to less optimal substrate binding in the expanded binding cavity, is primarily responsible for its higher kinetic isotope effect. The calculated potentials of mean force are anharmonic and relatively soft at shorter C-O distances, allowing efficient thermal sampling of the shorter distances required for effective hydrogen tunneling. The primarily local electrostatic field at the transferring hydrogen is ∼100 MV/cm in the direction to facilitate proton transfer and increases dramatically as the C-O distance decreases. These simulations suggest that the overall protein environment is important for conformational sampling of active substrate configurations aligned for proton transfer, but the PCET reaction is influenced primarily by local electrostatic effects that facilitate conformational sampling of shorter proton donor-acceptor distances required for effective hydrogen tunneling.