Derivation of rate expressions for nonadiabatic proton-coupled electron transfer reactions in solution

Controlling Electrons and Protons through Theory: Molecular Electrocatalysts to Nanoparticles

The development of renewable energy sources that are environmentally friendly and economical is of critical importance. The effective utilization of such energy sources relies on catalysts to facilitate the interconversion between electrical and chemical energy through multielectron, multiproton reactions. The design of effective catalysts for these types of energy conversion processes requires the ability to control the localization and movement of electrons and protons, as well as the coupling between them. Theoretical calculations, in conjunction with experimental validation and feedback, are playing a key role in these catalyst design efforts. A general theory has been developed for describing proton-coupled electron transfer (PCET) reactions, which encompass all reactions involving the coupled transfer of electrons and protons, including sequential and concerted mechanisms for multielectron, multiproton processes. In addition, computational methods have been devised to compute the input quantities for the PCET rate constant expressions and to generate free energy pathways for molecular electrocatalysts. These methods have been extended to heterogeneous PCET reactions to enable the modeling of PCET processes at electrode and nanoparticle surfaces. Three distinct theoretical studies of PCET reactions relevant to catalyst design for energy conversion processes are discussed. In the first application, theoretical calculations of hydrogen production catalyzed by hangman metalloporphyrins predicted that the porphyrin ligand is reduced, leading to dearomatization and proton transfer from the carboxylic acid hanging group to the meso carbon of the porphyrin rather than the metal center, producing a phlorin intermediate. Subsequent experiments isolated and characterized the phlorin intermediate, validating this theoretical prediction. These molecular electrocatalysts exemplify the potential use of noninnocent ligands to localize electrons and protons on different parts of the catalyst and to direct their motions accordingly. In the second application, theoretical calculations on substituted benzimidazole phenol molecules predicted that certain substituents would lead to multiple intramolecular proton transfer reactions upon oxidation. Subsequent experiments verified these multiproton reactions, as well as the predicted shifts in the redox potentials and kinetic isotope effects. These bioinspired molecular systems demonstrate the potential use of multiproton relays to enable the transport of protons over longer distances along specified pathways, as well as the tuning of redox potentials through this movement of positive charge. In the third application, theoretical studies of heterogeneous PCET in photoreduced ZnO nanoparticles illustrated the significance of proton diffusion through the bulk of the nanoparticle as well as interfacial PCET to an organic radical in solution at its surface. These theoretical calculations were consistent with prior experimental studies of this system, although theoretical methods for heterogeneous PCET have not yet attained the level of predictive capability highlighted for the molecular electrocatalysts. These examples suggest that theory will play a significant role in the design of both molecular and heterogeneous catalysts to control the movement and coupling of electrons and protons. The resulting catalysts will be essential for the development of renewable energy sources to address current energy challenges.

Diabatic model for electrochemical hydrogen evolution based on constrained DFT configuration interaction

The accuracy of density functional theory (DFT) based kinetic models for electrocatalysis is diminished by spurious electron delocalization effects, which manifest as uncertainties in the predicted values of reaction and activation energies. In this work, we present a constrained DFT (CDFT) approach to alleviate overdelocalization effects in the Volmer-Heyrovsky mechanism of the hydrogen evolution reaction (HER). This method is applied a posteriori to configurations sampled along a reaction path to correct their relative stabilities. Concretely, the first step of this approach involves describing the reaction in terms of a set of diabatic states that are constructed by imposing suitable density constraints on the system. Refined reaction energy profiles are then recovered by performing a configuration interaction (CDFT-CI) calculation within the basis spanned by the diabatic states. After a careful validation of the proposed method, we examined HER catalysis on open-ended carbon nanotubes and discovered that CDFT-CI increased activation energies and decreased reaction energies relative to DFT predictions. We believe that a similar approach could also be adopted to treat overdelocalization effects in other electrocatalytic proton-coupled electron transfer reactions, e.g., in the oxygen reduction reaction.

Proton-coupled electron transfer reactions: analytical rate constants and case study of kinetic isotope effects in lipoxygenase

A general theory has been developed for proton-coupled electron transfer (PCET), which is vital to a wide range of chemical and biological processes. This theory describes PCET reactions in terms of nonadiabatic transitions between reactant and product electron-proton vibronic states and includes the effects of thermal fluctuations of the solvent or protein environment, as well as the proton donor-acceptor motion. Within the framework of this general PCET theory, a series of analytical rate constant expressions has been derived for PCET reactions in well-defined regimes. Herein, the application of this theory to PCET in the enzyme soybean lipoxygenase illustrates the regimes of validity for the various rate constant expressions and elucidates the fundamental physical principles dictating PCET reactions. Such theoretical studies provide significant physical insights that guide the interpretation of experimental data and lead to experimentally testable predictions. A combination of theoretical treatments with atomic-level simulations is essential to understanding PCET.

Coupled-perturbed DFTB-QM/MM metadynamics: Application to proton-coupled electron transfer

We present a new concept of free energy calculations of chemical reactions by means of extended sampling molecular dynamics simulations. Biasing potentials are applied on partial atomic charges, which may be combined with atomic coordinates either in a single collective variable or in multi-dimensional biasing simulations. The necessary additional gradients are obtained by solving coupled-perturbed equations within the approximative density-functional tight-binding method. The new computational scheme was implemented in a combination of Gromacs and Plumed. As a prospective application, proton-coupled electron transfer in a model molecular system is studied. Two collective variables are introduced naturally, one for the proton transfer and the other for the electron transfer. The results are in qualitative agreement with the extended free simulations performed for reference. Free energy minima as well as the mechanism of the process are identified correctly, while the topology of the transition region and the height of the energy barrier are only reproduced qualitatively. The application also illustrates possible difficulties with the new methodology. These may be inefficient sampling of spatial coordinates when atomic charges are biased exclusively and a decreased stability of the simulations. Still, the new approach represents a viable alternative for free energy calculations of a certain class of chemical reactions, for instance a proton-coupled electron transfer in proteins.

Theoretical-computational modeling of charge transfer and intersystem crossing reactions in complex chemical systems

In this paper we present a theoretical-computational methodology specifically aimed at describing processes involving internal conversion or intersystem crossing, from atomistic (semiclassical) simulations and, hence, very suitable for treating complex atomic-molecular systems. The core of the presented approach is the evaluation of the diabatic perturbed energy surfaces of a portion of the whole system, treated at the quantum level and therefore preventively selected, in semi-classical interaction with the atomic-molecular environment. Subsequently, the estimation of the coupling between the diabatic surfaces and the inclusion of the obtained observables within a properly designed kinetic model allows the reconstruction of the whole phenomenology directly comparable to the experimental (typically kinetic) data. Application to two systems has demonstrated that the proposed approach can represent a valuable tool, somewhat complementary to other methods based on explicit quantum-dynamical approaches, for the theoretical-computational investigations of large and complex atomic-molecular systems.

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.

A qualitative quantum rate model for hydrogen transfer in soybean lipoxygenase

The hydrogen transfer reaction catalysed by soybean lipoxygenase (SLO) has been the focus of intense study following observations of a high kinetic isotope effect (KIE). Today high KIEs are generally thought to indicate departure from classical rate theory and are seen as a strong signature of tunnelling of the transferring particle, hydrogen or one of its isotopes, through the reaction energy barrier. In this paper, we build a qualitative quantum rate model with few free parameters that describes the dynamics of the transferring particle when it is exposed to energetic potentials exerted by the donor and the acceptor. The enzyme's impact on the dynamics is modelled by an additional energetic term, an oscillatory contribution known as "gating." By varying two key parameters, the gating frequency and the mean donor-acceptor separation, the model is able to reproduce well the KIE data for SLO wild-type and a variety of SLO mutants over the experimentally accessible temperature range. While SLO-specific constants have been considered here, it is possible to adapt these for other enzymes.

Theoretical Insights into Proton-Coupled Electron Transfer from a Photoreduced ZnO Nanocrystal to an Organic Radical

Proton-coupled electron transfer (PCET) at metal-oxide nanoparticle interfaces plays a critical role in many photocatalytic reactions and energy conversion processes. Recent experimental studies have shown that photoreduced ZnO nanocrystals react by PCET with organic hydrogen atom acceptors such as the nitroxyl radical TEMPO. Herein, the interfacial PCET rate constant is calculated in the framework of vibronically nonadiabatic PCET theory, which treats the electrons and transferring proton quantum mechanically. The input quantities to the PCET rate constant, including the electronic couplings, are calculated with density functional theory. The computed interfacial PCET rate constant is consistent with the experimentally measured value for this system, providing validation for this PCET theory. In this model, the electron transfers from the conduction band of the ZnO nanocrystal to TEMPO concertedly with proton transfer from a surface oxygen of the ZnO nanocrystal to the oxygen of TEMPO. Moreover, the proton tunneling at the interface is gated by the relatively low-frequency proton donor-acceptor motion between the TEMPO radical and the ZnO nanocrystal. The ZnO nanocrystal and TEMPO are found to contribute similar amounts to the inner-sphere reorganization energy, implicating structural reorganization at the nanocrystal surface. These fundamental mechanistic insights may guide the design of metal-oxide nanocatalysts for a wide range of energy conversion processes.

Proton transfer dynamics control the mechanism of O2 reduction by a non-precious metal electrocatalyst

Many chemical and biological processes involve the transfer of both protons and electrons. The complex mechanistic details of these proton-coupled electron transfer (PCET) reactions require independent control of both electron and proton transfer. In this report, we make use of lipid-modified electrodes to modulate proton transport to a Cu-based catalyst that facilitates the O2 reduction reaction (ORR), a PCET process important in fuel cells and O2 reduction enzymes. By quantitatively controlling the kinetics of proton transport to the catalyst, we demonstrate that undesired side products such as H2O2 and O2(-) arise from a mismatch between proton and electron transfer rates. Whereas fast proton kinetics induce H2O2 formation and sluggish proton flux produces O2(-), proton transfer rates commensurate with O-O bond breaking rates ensure that only the desired H2O product forms. This fundamental insight aids in the development of a comprehensive framework for understanding the ORR and PCET processes in general.

Dependence of Vibronic Coupling on Molecular Geometry and Environment: Bridging Hydrogen Atom Transfer and Electron-Proton Transfer

The rate constants for typical concerted proton-coupled electron transfer (PCET) reactions depend on the vibronic coupling between the diabatic reactant and product states. The form of the vibronic coupling is different for electronically adiabatic and nonadiabatic reactions, which are associated with hydrogen atom transfer (HAT) and electron-proton transfer (EPT) mechanisms, respectively. Most PCET rate constant expressions rely on the Condon approximation, which assumes that the vibronic coupling is independent of the nuclear coordinates of the solute and the solvent or protein. Herein we test the Condon approximation for PCET vibronic couplings. The dependence of the vibronic coupling on molecular geometry is investigated for an open and a stacked transition state geometry of the phenoxyl-phenol self-exchange reaction. The calculations indicate that the open geometry is electronically nonadiabatic, corresponding to an EPT mechanism that involves significant electronic charge redistribution, while the stacked geometry is predominantly electronically adiabatic, corresponding primarily to an HAT mechanism. Consequently, a single molecular system can exhibit both HAT and EPT character. The dependence of the vibronic coupling on the solvent or protein configuration is examined for the soybean lipoxygenase enzyme. The calculations indicate that this PCET reaction is electronically nonadiabatic with a vibronic coupling that does not depend significantly on the protein environment. Thus, the Condon approximation is shown to be valid for the solvent and protein nuclear coordinates but invalid for the solute nuclear coordinates in certain PCET systems. These results have significant implications for the calculation of rate constants, as well as mechanistic interpretations, of PCET reactions.

Proton-Coupled Electron Transfer: Moving Together and Charging Forward

Proton-coupled electron transfer (PCET) is ubiquitous throughout chemistry and biology. This Perspective discusses recent advances and current challenges in the field of PCET, with an emphasis on the role of theory and computation. The fundamental theoretical concepts are summarized, and expressions for rate constants and kinetic isotope effects are provided. Computational methods for calculating reduction potentials and pKa's for molecular electrocatalysts, as well as insights into linear correlations and non-innocent ligands, are also described. In addition, computational methods for simulating the nonadiabatic dynamics of photoexcited PCET are discussed. Representative applications to PCET in solution, proteins, electrochemistry, and photoinduced processes are presented, highlighting the interplay between theoretical and experimental studies. The current challenges and suggested future directions are outlined for each type of application, concluding with an overall view to the future.

Proton-coupled electron transfer in molecular electrocatalysis: theoretical methods and design principles

Molecular electrocatalysts play an essential role in a wide range of energy conversion processes. The objective of electrocatalyst design is to maximize the turnover frequency and minimize the overpotential for the overall catalytic cycle. Typically, the catalytic cycle is dominated by key proton-coupled electron transfer (PCET) processes comprised of sequential or concerted electron and proton transfer steps. Theoretical methods have been developed to investigate the mechanisms, thermodynamics, and kinetics of PCET processes in electrocatalytic cycles. Electronic structure methods can be used to calculate the reduction potentials and pKa's and to generate thermodynamic schemes, free energy reaction pathways, and Pourbaix diagrams, which indicate the most stable species under certain conditions. These types of calculations have assisted in identifying the thermodynamically favorable mechanisms under specified experimental conditions, such as acid strength and overpotential. Such calculations have also revealed linear correlations among the thermodynamic properties, which can be used to predict the impact of modifying the ligands, substituents, or metal centers. The thermodynamic properties can be tuned with electron-withdrawing or electron-donating substituents. Ligand modification can exploit the role of noninnocent ligands. For example, ligand protonation typically decreases the overpotential. Calculations of rate constants for electron and proton transfer, as well as concerted PCET, have assisted in identifying the kinetically favorable mechanisms under specified conditions. The concerted PCET mechanism is thought to lower the overpotential required for catalysis by avoiding high-energy intermediates. Rate constant calculations have revealed that the concerted mechanism involving intramolecular proton transfer will be favored by designing more flexible ligands that facilitate the proton donor-acceptor motion while also maintaining a sufficiently short equilibrium proton donor-acceptor distance. Overall, theoretical methods have assisted in the interpretation of experimental data and the design of more effective molecular electrocatalysts.

Concerted proton-electron transfers: fundamentals and recent developments

Proton-coupled electron transfers (PCET) are ubiquitous in natural and synthetic processes. This review focuses on reactions where the two events are concerted. Semiclassical models of such reactions allow their kinetic characterization through activation versus driving force relationships, estimates of reorganization energies, effects of the nature of the proton acceptor, and H/D kinetic isotope effect as well as their discrimination from stepwise pathways. Several homogeneous reactions (through stopped-flow and laser flash-quench techniques) and electrochemical processes are discussed in this framework. Once the way has been rid of the improper notion of pH-dependent driving force, water appears as a remarkable proton acceptor in terms of reorganization energy and pre-exponential factor, thanks to its H-bonded and H-bonding properties, similarly to purposely synthesized "H-bond train" molecules. The most recent developments are in modeling and description of emblematic concerted proton-electron transfer (CPET) reactions associated with the breaking of a heavy-atom bond in an all-concerted process.

Theory of the transition from sequential to concerted electrochemical proton-electron transfer

A theory for the calculation of potential energy surfaces of electrochemical proton-coupled electron transfer is considered and parameterized on the basis of thermodynamic relations. The paper discusses the qualitatively different potential energy surfaces predicted by the theory, and their relation to the existence of sequential and concerted proton-electron transfer pathways. The concomitant activation energies for sequential and concerted PET are calculated. The applied overpotential may change the qualitative shape of the PES and therefore the mechanism of the proton-coupled electron transfer reaction.

Gaining insight into the chemistry of lipoxygenases: a computational investigation into the catalytic mechanism of (8R)-lipoxygenase

Lipoxygenases (LOXs) are ubiquitous in nature and catalyze a range of life-essential reactions within organisms. In particular they are critical to the formation of eicosanoids, which are critical for normal cell function. However, a number of important questions about the reactivity and mechanism of these enzymes still remain. Specifically, although the initial step in the mechanism of LOXs has been well studied, little is known of subsequent steps. Thus, with use of a quantum mechanical/molecular mechanical approach, the complete catalytic mechanism of (8R)-LOX was investigated. The results have provided a better understanding of the general chemistry of LOXs as a whole. In particular, from comparisons with soybean LOX-1, it appears that the initial proton-coupled electron transfer may be very similar among all LOXs. Furthermore, LOXs appear to undergo multistate reactivity where potential spin inversion of an electron may occur either in the attack of O(2) or in the regeneration of the active site. Lastly, it is shown that with the explicit modeling of the environment, the regeneration of the active center likely occurs via the rotation of the intermediate followed by an outer-sphere [Formula: see text] transfer as opposed to the formation of a "purple intermediate" complex.

Impaired protein conformational landscapes as revealed in anomalous Arrhenius prefactors

A growing body of data supports a role for protein motion in enzyme catalysis. In particular, the ability of enzymes to sample catalytically relevant conformational substates has been invoked to model kinetic and spectroscopic data. However, direct experimental links between rapidly interconverting conformations and the chemical steps of catalysis remain rare. We report here on the kinetic analysis and characterization of the hydride transfer step catalyzed by a series of mutant thermophilic alcohol dehydrogenases (ht-ADH), presenting evidence for Arrhenius prefactor values that become enormously elevated above an expected value of approximately 10(13) s(-1) when the enzyme operates below its optimal temperature range. Restoration of normal Arrhenius behavior in the ht-ADH reaction occurs at elevated temperatures. A simple model, in which reduced temperature alters the ability of the ht-ADH variants to sample the catalytically relevant region of conformational space, can reproduce the available data. These findings indicate an impaired landscape that has been generated by the combined condition of reduced temperature and mutation at a single, active-site hydrophobic side chain. The broader implication is that optimal enzyme function requires the maintenance of a relatively smooth landscape that minimizes low energy traps.

The diabatic picture of electron transfer, reaction barriers, and molecular dynamics

Diabatic states have a long history in chemistry, beginning with early valence bond pictures of molecular bonding and extending through the construction of model potential energy surfaces to the modern proliferation of methods for computing these elusive states. In this review, we summarize the basic principles that define the diabatic basis and demonstrate how they can be applied in the specific context of constrained density functional theory. Using illustrative examples from electron transfer and chemical reactions, we show how the diabatic picture can be used to extract qualitative insight and quantitative predictions about energy landscapes. The review closes with a brief summary of the challenges and prospects for the further application of diabatic states in chemistry.

Impact of distal mutation on hydrogen transfer interface and substrate conformation in soybean lipoxygenase

The impact of distal mutation on the hydrogen transfer interface properties and on the substrate mobility, conformation, and orientation in soybean lipoxygenase-1 (SLO) is examined. SLO catalyzes a hydrogen abstraction reaction that occurs by a proton-coupled electron transfer mechanism. Mutation of isoleucine 553 to less bulky residues has been found experimentally to increase the magnitude and temperature dependence of the kinetic isotope effect for this reaction. This residue borders the linoleic acid substrate but is approximately 15 A from the active site iron. In the present study, we model these experimental data with a vibronically nonadiabatic theory and perform all-atom molecular dynamics simulations on the complete solvated wild-type and mutant enzymes. Our calculations indicate that the proton transfer equilibrium distance increases and the associated frequency decreases as residue 553 becomes less bulky. The molecular dynamics simulations illustrate that this mutation impacts the mobility, geometrical conformation, and orientation of the linoleic acid within the active site. In turn, these effects alter the proton donor-acceptor equilibrium distance and frequency, leading to the experimentally observed changes in the magnitude and temperature dependence of the kinetic isotope effect. This study provides insight into how the effects of distal mutations may be transmitted in enzymes to ultimately impact the catalytic rates.

Theory of proton-coupled electron transfer in energy conversion processes

Proton-coupled electron transfer (PCET) reactions play an essential role in a broad range of energy conversion processes, including photosynthesis and respiration. These reactions also form the basis of many types of solar fuel cells and electrochemical devices. Recent advances in the theory of PCET enable the prediction of the impact of system properties on the reaction rates. These predictions may guide the design of more efficient catalysts for energy production, including those based on artificial photosynthesis and solar energy conversion. This Account summarizes the theoretically predicted dependence of PCET rates on system properties and illustrates potential approaches for tuning the reaction rates in chemical systems. A general theoretical formulation for PCET reactions has been developed over the past decade. In this theory, PCET reactions are described in terms of nonadiabatic transitions between the reactant and product electron-proton vibronic states. A series of nonadiabatic rate constant expressions for both homogeneous and electrochemical PCET reactions have been derived in various well-defined limits. Recently this theory has been extended to include the effects of solvent dynamics and to describe ultrafast interfacial PCET. Analysis of the rate constant expressions provides insight into the underlying physical principles of PCET and enables the prediction of the dependence of the rates on the physical properties of the system. Moreover, the kinetic isotope effect, which is the ratio of the rates for hydrogen and deuterium, provides a useful mechanistic probe. Typically the PCET rate will increase as the electronic coupling and temperature increase and as the total reorganization energy and equilibrium proton donor-acceptor distance decrease. The rate constant is predicted to increase as the driving force becomes more negative, rather than exhibit turnover behavior in the inverted region, because excited vibronic product states associated with low free energy barriers and relatively large vibronic couplings become accessible. The physical basis for the experimentally observed pH dependence of PCET reactions has been debated in the literature. When the proton acceptor is a buffer species, the pH dependence may arise from the protonation equilibrium of the buffer. It could also arise from kinetic complexity of competing concerted and sequential PCET reaction pathways. In electrochemical PCET, the heterogeneous rate constants and current densities depend strongly on the overpotential. The change in equilibrium proton donor-acceptor distance upon electron transfer may lead to asymmetries in the Tafel plots and deviations of the transfer coefficient from the standard value of one-half at zero overpotential. Applications of this theory to experimentally studied systems illustrate approaches that can be utilized to tune the PCET rate. For example, the rate can be tuned by changing the pH or using different buffer species as proton acceptors. The rate can also be tuned with site-specific mutagenesis in biological systems or chemical modifications that vary the substituents on the redox species in chemical systems. Understanding the impact of these changes on the PCET rate may assist experimental efforts to enhance energy conversion processes.

Driving force dependence of rates for nonadiabatic proton and proton-coupled electron transfer: conditions for inverted region behavior

The driving force dependence of the rate constants for nonadiabatic electron transfer (ET), proton transfer (PT), and proton-coupled electron transfer (PCET) reactions is examined. Inverted region behavior, where the rate constant decreases as the reaction becomes more exoergic (i.e., as DeltaG(0) becomes more negative), has been observed experimentally for ET and PT. This behavior was predicted theoretically for ET but is not well understood for PT and PCET. The objective of this Letter is to predict the experimental conditions that could lead to observation of inverted region behavior for PT and PCET. The driving force dependence of the rate constant is qualitatively different for PT and PCET than for ET because of the high proton vibrational frequency and substantial shift between the reactant and product proton vibrational wave functions. As a result, inverted region behavior is predicted to be experimentally inaccessible for PT and PCET if only the driving force is varied. This behavior may be observed for PT over a limited range of rates and driving forces if the solvent reorganization energy is low enough to cause observable oscillations. Moreover, this behavior may be observed for PT or PCET if the proton donor-acceptor distance increases as DeltaG(0) becomes more negative. Thus, a plausible explanation for experimentally observed inverted region behavior for PT or PCET is that varying the driving force also impacts other properties of the system, such as the proton donor-acceptor distance.