Abstract:Multistep electrocatalytic cascades are capable of expanding the range and extent of electrochemical reactions, leading to fuel cells with high energy density and sensors with enhanced sensitivity[1]. In such processes, a key limitation is the mass transport of reaction intermediate between spatially organized active sites. Substrate channeling, by which reaction intermediates are directly transported to downstream active sites without equilibrating to the bulk environment, is a key aspect of efficient cascades[2]. However, the application to artificial systems is hindered by the need for molecular-level design and experiments often on nanosecond time scales. Multiscale simulation of multistep cascades bridges the gap between molecular interaction and experiment results. Channeling via surface interactions involves simultaneously occurring reaction, adsorption, desorption, surface diffusion and bulk diffusion. Kinetic Monte Carlo (KMC) method is widely used to characterize the overall kinetics of heterogeneous catalysis with complex and coupled events[3]. Here, we describe the use of KMC simulation, aided with Molecular Dynamics (MD), to study the channeling in an enzymatic complex system with electrostatic bridge, aiming to quantify the overall kinetics, identify the potential limitation step, and bridge the gap between molecular simulation and experiment result (e.g., stop-flow lag time experiment). In this project, we proposed a channeling model with discrete reaction and hopping sites, as shown in Fig. 1a. Events on each individual site is allowed to interact with the external environment (e.g., bulk intermediate concentration). For rate-limiting enzyme kinetics, Fig. 1c shows an example of lag time, an experimentally accessible metric, on channeling efficiency. Specifically, when the desorption rate for intermediate on E2 is low (k2 = 0.1 kcat), the lag time significantly decreases (from 212 s to 19 s) with improved substrate channeling. However, Fig. 1e shows that the bridge has minimal impact when the desorption rates on E2 represents typical enzyme systems (k2 = 10 kcat). This indicates that in an ideal channeling system, the functional bridge should be able to protect the adsorbed intermediate from equilibrating to bulk environment[4,5]. In addition, if the bulk intermediate is allow to re-adsorb onto the bridge, the bridge impact saw a great increase (Fig.1e) in spite of high k2 value. Therefore, the behavior of intermediate molecules around bridge/enzyme boundary is of great significance. Based on these results, we are able to identify the potential limiting factors on overall kinetics, which highlights the primary interactions to be studied by integrated molecular simulation and experiment. Acknowledgement We gratefully acknowledge support from Army Research Office MURI (#W911NF1410263) via The University of Utah References 1. M. J. Moehlenbrock, M. T. Meredith, and S. D. Minteer, ACS Catalysis, 2, 17–25 (2012) http://pubs.acs.org/doi/abs/10.1021/cs200482v. 2. I. Wheeldon et al., Nature Chemistry, 8, 299–309 (2016) http://www.nature.com/doifinder/10.1038/nchem.2459. 3. M. Stamatakis and D. G. Vlachos, ACS Catalysis, 2, 2648–2663 (2012) http://pubs.acs.org/doi/abs/10.1021/cs3005709. 4. A. H. Elcock, G. A. Huber, and J. A. McCammon, Biochemistry, 36, 16049–16058 (1997) http://pubs.acs.org/doi/abs/10.1021/bi971709u. 5. M. Trujillo, R. G. K. Donald, D. S. Roos, P. J. Greene, and D. V. Santi, Biochemistry, 35, 6366–6374 (1996) http://pubs.acs.org/doi/abs/10.1021/bi952923q. Figure 1

Molecular Dynamics Simulation of Substrate Transport on Alpha-Helix Peptide for Multistep Catalysis

Molecular Dynamics Simulations on Apo ADP/ATP Carrier Shed New Lights on the Featured Motif of the Mitochondrial Carriers.

Kinetic Monte Carlo Simulation of Substrate Channeling Via Electrostatic Interactions

Substrate Channeling in an Artificial Metabolon: A Molecular Dynamics Blueprint for an Experimental Peptide Bridge

Intermediate Hopping on Flexible Cascade Surface with Electrostatic Channeling

Markov-State Transition Path Analysis of Electrostatic Channeling

Quantitative Assessment of Channeling Mechanisms in Nanoscale Catalytic Architectures

Energetics of substrate transport in proton-dependent oligopeptide transporters

Structural studies of catalytic peptides using molecular dynamics simulations

Novel Quantification of Cascade Kinetics of Electrostatic Channeling

Molecular Dynamics for Electrocatalysis: Mechanism Explanation and Performance Prediction

Molecular Dynamics Simulation on Quaternized Ammonium Poly (ether Ether Ketone) for Anion-Exchange Membrane Fuel Cells

Resolving Chemical Dynamics in Biological Energy Conversion: Long-Range Proton-Coupled Electron Transfer in Respiratory Complex I

Qm/Mm Study Of The Enzymatic Catalysis

Electrostatic Interactions Determine Entrance/Release order of Substrates in the Catalytic Cycle of Adenylate Kinase.

Simulation of Proton Transport in Proton Exchange Membranes with Reactive Molecular Dynamics

Dynamic Charge Distribution as a Key Driver of Catalytic Reactivity in an Artificial Metalloenzyme

(Invited) First-Principles Molecular Dynamics Investigations of Proton and Hydroxide Transport in Model Anion-Exchange- and Proton-Exchange Membranes in Nanoconfined, Low-Hydration Environments

Brownian dynamic study of an enzyme metabolon in the TCA cycle: Substrate kinetics and channeling

Modulation of Proton-Coupled Electron Transfer Reactions in Lysine-Containing Alpha-Helixes: Alpha-Helixes Promoting Long-Range Electron Transfer.

Global and Kinetic Profiles of Substrate Diffusion in Candida Antarctica Lipase B: Molecular Dynamics with the Markov-State Model