Analysis of the Acting Forces in a Theory of Catalysis and Mechanochemistry

Theory and Examples of Catch Bonds

We discuss slip bonds, catch bonds, and the tug-of-war mechanism using mathematical arguments. The aim is to explain the theoretical tool of molecular potential energy surfaces (PESs). For this, we propose simple 2-dimensional surface models to demonstrate how a molecule under an external force behaves. Examples are selectins. Catch bonds, in particular, are explained in more detail, and they are contrasted to slip bonds. We can support special two-dimensional molecular PESs for E- and L-selectin which allow the catch bond property. We demonstrate that Newton trajectories (NT) are powerful tools to describe these phenomena. NTs form the theoretical background of mechanochemistry.

Topology of molecular deformations induces triphasic catch bonding in selectin-ligand bonds

Among the long-standing efforts to elucidate the physical mechanisms of protein-ligand catch bonding, particular attention has been directed at the family of selectin proteins. Selectins exhibit slip, catch-slip, and slip-catch-slip bonding, with minor structural modifications causing major changes in selectins' response to force. How can a single structural mechanism allow interconversion between these various behaviors? We present a unifying theory of selectin-ligand catch bonding, using a structurally motivated free energy landscape to show how the topology of force-induced deformations of the molecular system produces the full range of observed behaviors. We find that the pathway of bond rupture deforms in non-trivial ways, such that unbinding dynamics depend sensitively on force. This implies a severe breakdown of Bell's theory-a paradigmatic theory used widely in catch bond modeling-raising questions about the suitability of Bell's theory in modeling other catch bonds. Our approach can be applied broadly to other protein-ligand systems.

Controlling Chemical Reactivity with Optimally Oriented Electric Fields: A Generalization of the Newton Trajectory Method

The use of oriented external electric fields (OEEF) as a tool to accelerate chemical reactions has recently attracted much interest. A new model to calculate the optimal OEEF of the least intensity to induce a barrierless chemical reaction path is presented. A suitable ansatz is provided by defining an effective potential energy surface (PES), which considers the unperturbed or original PES of the molecular reactive system and the action of a constant OEEF on the overall dipole moment of system. Based on a generalization of the Newton Trajectories (NT) method, it is demonstrated that the optimal OEEF can be determined upon locating a special point of the potential energy surface (PES), the so-called "optimal bond-breaking point" (optimal BBP), for which two different algorithms are proposed. At this point, the gradient of the original or unperturbed PES is an eigenvector of zero eigenvalue of the Hessian matrix of the effective PES. A thorough discussion of the geometrical aspects of the optimal BBP and the optimal OEEF is provided using a two-dimensional model, and numerical calculations of the optimal OEEF for a S_N2 reaction and the 1,3-dipolar retrocycloaddition of isoxazole to fulminic acid plus acetylene reaction serve as a proof of concept. The knowledge of the orientation of optimal OEEF provides a practical way to reduce the effective barrier of a given chemical process.

Dissecting Biological and Synthetic Soft-Hard Interfaces for Tissue-Like Systems

Soft and hard materials at interfaces exhibit mismatched behaviors, such as mismatched chemical or biochemical reactivity, mechanical response, and environmental adaptability. Leveraging or mitigating these differences can yield interfacial processes difficult to achieve, or inapplicable, in pure soft or pure hard phases. Exploration of interfacial mismatches and their associated (bio)chemical, mechanical, or other physical processes may yield numerous opportunities in both fundamental studies and applications, in a manner similar to that of semiconductor heterojunctions and their contribution to solid-state physics and the semiconductor industry over the past few decades. In this review, we explore the fundamental chemical roles and principles involved in designing these interfaces, such as the (bio)chemical evolution of adaptive or buffer zones. We discuss the spectroscopic, microscopic, (bio)chemical, and computational tools required to uncover the chemical processes in these confined or hidden soft-hard interfaces. We propose a soft-hard interaction framework and use it to discuss soft-hard interfacial processes in multiple systems and across several spatiotemporal scales, focusing on tissue-like materials and devices. We end this review by proposing several new scientific and engineering approaches to leveraging the soft-hard interfacial processes involved in biointerfacing composites and exploring new applications for these composites.

Understanding Ball Milling Mechanochemical Processes with DFT Calculations and Microkinetic Modeling

Mechanochemistry is an emerging field with many potential applications in sustainable chemistry. But despite the growing interest in the field, its underlying mechanistic foundations are not fully understood yet. This work presents the application of computational tools, such as DFT calculations in continuum and microkinetic modeling, to the analysis of mechanically activated procedures. Two reactions reported in previous experimental publications were studied: (i) a series of Diels-Alder reactions and (ii) the synthesis of sulfonylguanidines. Calculations succeed in reproducing experimentally reported reaction times. The procedures were mostly standard, coupled with some sensitive choices in terms of starting concentrations and dielectric constant. This means that these particular reactions accelerated by ball milling followed the same mechanism as the equivalent reactions in solution. The implications of this result on the general picture of mechanochemical processes are discussed.

Barnes Update Applied in the Gauss-Newton Method: An Improved Algorithm to Locate Bond Breaking Points

A mechanochemical reaction is a reaction induced by mechanical energy. A general accepted model for this type of reaction consists of a first-order perturbation on the associated potential energy surface (PES) of the unperturbed molecular system due to mechanical stress or pulling force. Within this theoretical framework, the so-called optimal barrier breakdown points or optimal bond breaking points (BBPs) are critical points of the unperturbed PES where the Hessian matrix has a zero eigenvector that coincides with the gradient vector. Optimal BBPs are "catastrophe points" that are particularly important because their associated gradient indicates how to optimally harness tensile forces to induce reactions by transforming a chemical reaction into a barrierless process. Building on a previous method based on a nonlinear least-squares minimization to locate BBPs (Bofill et al., J. Chem. Phys. 2017, 147, 152710-10), we propose a new algorithm to locate BBPs of any molecular system based on the Gauss-Newton method combined with the Barnes update for a nonsymmetric Jacobian matrix, which is shown to be more appropriate than the Broyden update. The efficiency of the new method is demonstrated for a multidimensional model PES and two medium size molecular systems of interest in enzymatic catalysis and mechanochemistry.

Some Mathematical Reasoning on the Artificial Force Induced Reaction Method

There are works of the Maeda-Morokuma group, which propose the artificial force induced reaction (AFIR) method (Maeda et al., J. Comput. Chem. 2014, 35, 166 and 2018, 39, 233). We study this important method from a theoretical point of view. The understanding of the proposers does not use the barrier breakdown point of the AFIR parameter, which usually is half of the reaction path between the minimum and the transition state which is searched for. Based on a comparison with the theory of Newton trajectories, we could better understand the method. It allows us to follow along some reaction pathways from minimum to saddle point, or vice versa. We discuss some well-known two-dimensional test surfaces where we calculate full AFIR pathways. If one has special AFIR curves at hand, one can also study the behavior of the ansatz. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Comment on "Exploring Potential Energy Surface with External Forces"

Recently, a work (Wolinski, K., J. Chem. Theory Comput. 2018, 14, 6306, 10.1021/acs.jctc.8b00885 ) was published in which the SEGO method (standard and enforced geometry optimization) was proposed to find new minimums on potential energy surfaces. We study this important method from a theoretical point of view. Up to now, the understanding of the proposer does not take into account the barrier breakdown point on a SEGO path being usually half of the path, which is searched for. However, a better understanding of the method allows us to follow along the reaction pathway from a minimum to a saddle point or vice versa. We discuss the well-known two-dimensional MB test surface where we calculate full SEGO pathways. If one has special SEGO curves at hand, one can also detect some weaknesses of the ansatz.

Plasmonic Nanocavities Enable Self-Induced Electrostatic Catalysis

The potential of strong interactions between light and matter remains to be further explored within a chemical context. Towards this end herein we study the electromagnetic interaction between molecules and plasmonic nanocavities. By means of electronic structure calculations, we show that self-induced catalysis emerges without any external stimuli through the interaction of the molecular permanent and fluctuating dipole moments with the plasmonic cavity modes. We also exploit this scheme to modify the transition temperature T1/2 of spin-crossover complexes as an example of how strong light-matter interactions can ultimately be used to control a materials responses.

Force-conductance spectroscopy of a single-molecule reaction

We demonstrate how simultaneous measurements of conductance and force can be used to monitor the step-by-step progress of a mechanically-activated cis-to-trans isomerization single-molecule reaction, including events that cannot be distinguished using force or conductance alone. To do so, we simulated the force-conductance profile of cyclopropane oligomers connected to graphene nanoribbon electrodes that undergo a cis-to-trans isomerization during mechanical elongation. This was done using a combination of classical molecular dynamics simulation of the pulling using a reactive force field, and Landauer transport computations of the conductance with nonequilibrium Green's function methods. The isomerization events can be distinguished in both force and conductance profiles. However, the conductance profile during the mechanical elongation distinguishes between reaction intermediates that cannot be resolved using force. In turn, the force signals non-reactive deformations in the molecular backbone which are not visible in the conductance profile. These observations are shown to be robust to the choice of electrode and Hamiltonian model. The computations exemplify the potential of the integration of covalent mechanochemistry with molecular conductance to investigate chemical reactivity at the single-entity limit.

An algorithm to locate optimal bond breaking points on a potential energy surface for applications in mechanochemistry and catalysis

The reaction path of a mechanically induced chemical transformation changes under stress. It is well established that the force-induced structural changes of minima and saddle points, i.e., the movement of the stationary points on the original or stress-free potential energy surface, can be described by a Newton Trajectory (NT). Given a reactive molecular system, a well-fitted pulling direction, and a sufficiently large value of the force, the minimum configuration of the reactant and the saddle point configuration of a transition state collapse at a point on the corresponding NT trajectory. This point is called barrier breakdown point or bond breaking point (BBP). The Hessian matrix at the BBP has a zero eigenvector which coincides with the gradient. It indicates which force (both in magnitude and direction) should be applied to the system to induce the reaction in a barrierless process. Within the manifold of BBPs, there exist optimal BBPs which indicate what is the optimal pulling direction and what is the minimal magnitude of the force to be applied for a given mechanochemical transformation. Since these special points are very important in the context of mechanochemistry and catalysis, it is crucial to develop efficient algorithms for their location. Here, we propose a Gauss-Newton algorithm that is based on the minimization of a positively defined function (the so-called σ-function). The behavior and efficiency of the new algorithm are shown for 2D test functions and for a real chemical example.