Role of Microsolvation and Quantum Effects in the Accurate Prediction of Kinetic Isotope Effects: The Case of Hydrogen Atom Abstraction in Ethanol by Atomic Hydrogen in Aqueous Solution

Solvate Suite: A Command-Line Interface for Molecular Simulations and Multiscale Microsolvation Modeling

In this work, we introduce the Solvate Suite, a comprehensive and modular command-line interface designed for molecular simulation and microsolvation modeling. The suite interfaces with widely used scientific software, streamlining computational experiments for liquid systems through the automated creation of simulation boxes and topology with adjustable simulation parameters. Furthermore, it has features for graphical and statistical analysis of simulated properties and extraction of trajectory configurations with various filters. Additionally, it introduces innovative strategies for microsolvation modeling with a multiscale approach, employing equilibrated dynamics to identify favorable solute-solvent interactions and enabling full cluster optimization for free-energy calculations without imaginary frequency contamination.

Reactions of a hydrogen atom with haloacetates in aqueous solutions: Computational evidence for proton-coupled electron transfer and competing mechanisms

A computational study of the mechanisms and kinetics of the aqueous reactions of a hydrogen atom with haloacetates is presented. Several mechanisms in the close competition are observed, such as proton-coupled electron transfer (PCET), hydrogen atom transfer (HAT), and halogen abstraction (XA). Computations predict that dechlorination takes place via PCET mechanisms and not via XA, as stated earlier, while XA is the fastest mechanism forIAc - . The reaction rate constants are reasonably well predicted within the theoretically most reliable canonical variational transition state theory with small curvature tunneling corrections and compared with the experimental ones. To reproduce the experimental rate constants of the debromination process it is necessary to include the PCET and XA cumulative values. Small curvature tunneling corrections to the rate constants are the highest for HAT and PCET mechanisms, up to 70 times larger than the Wigner, while variational effects for XA mechanisms are very small.

Combined QM/MM, Machine Learning Path Integral Approach to Compute Free Energy Profiles and Kinetic Isotope Effects in RNA Cleavage Reactions

We present a fast, accurate, and robust approach for determination of free energy profiles and kinetic isotope effects for RNA 2'-O-transphosphorylation reactions with inclusion of nuclear quantum effects. We apply a deep potential range correction (DPRc) for combined quantum mechanical/molecular mechanical (QM/MM) simulations of reactions in the condensed phase. The method uses the second-order density-functional tight-binding method (DFTB2) as a fast, approximate base QM model. The DPRc model modifies the DFTB2 QM interactions and applies short-range corrections to the QM/MM interactions to reproduce ab initio DFT (PBE0/6-31G*) QM/MM energies and forces. The DPRc thus enables both QM and QM/MM interactions to be tuned to high accuracy, and the QM/MM corrections are designed to smoothly vanish at a specified cutoff boundary (6 Å in the present work). The computational speed-up afforded by the QM/MM+DPRc model enables free energy profiles to be calculated that include rigorous long-range QM/MM interactions under periodic boundary conditions and nuclear quantum effects through a path integral approach using a new interface between the AMBER and i-PI software. The approach is demonstrated through the calculation of free energy profiles of a native RNA cleavage model reaction and reactions involving thio-substitutions, which are important experimental probes of the mechanism. The DFTB2+DPRc QM/MM free energy surfaces agree very closely with the PBE0/6-31G* QM/MM results, and it is vastly superior to the DFTB2 QM/MM surfaces with and without weighted thermodynamic perturbation corrections. 18O and 34S primary kinetic isotope effects are compared, and the influence of nuclear quantum effects on the free energy profiles is examined.

An integrated protocol to study hydrogen abstraction reactions by atomic hydrogen in flexible molecules: application to butanol isomers

This work presents a protocol designed to study hydrogen abstraction reactions by atomic hydrogen in molecules with multiple conformations.

Reply to the 'Comment on "Topography of the Free Energy Landscape on the Claisen-Schmidt Condensation: Solvent and Temperature Effect in the Rate-Controlling Step"' by N. D. Coutinho, H. G. Machado, V. H. Carvalho-Silva and W. A. da Silva, Phys. Chem. Chem. Phys., 2021, , 6738

In the Comment on our paper on the description of the Gibbs Free energy profile of Claisen-Schmidt condensation, it is claimed that our calculations are flawed due to inconsistencies with experimental isotope effects in aqueous acetonitrile. In this Reply, we presented rigorous arguments, ambiguities in the Comment and new calculations confirming the consistency of our results: (i) small differences in the relative energetic barriers in both experimental and theoretical curves make the assignment of the rate-limiting step debatable, making the concept of RCS questionable; (ii) it is shown how the misinterpretation of the elementary steps and of the overall processes rate constants led the Comment to incorrect conclusions about the behavior of the inverse isotopic effect; (iii) neglect in the Comment of the inverse kinetic isotope effect in step R2 due to the hybridization conversion, and of the inverse equilibrium isotopic effect for step R1 to describe an overall iKIE > 1, (iv) an erroneous suggestion in the Comment that the disagreement between experimental kinetic parameters is due to the fact that acetonitrile is not used in previous experimental works, when contradictorily the literature recommends it as being indispensable to allow kinetic accuracy; and (v) new calculations improved by explicit-implicit hybrid treatment again ensure that step R4, and not step R5, can assume the role of RCS in protic solvents. Recognizing that questioning is an excellent path for promoting understanding, we hope that the answers provided here will help to clarify and expand the pertinent topics under discussion.

Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations

Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics-molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.

Kinetics of chain reaction driven by proton-coupled electron transfer: α-hydroxyethyl radical and bromoacetate in buffered aqueous solutions

We measured and computed the rate constants of the reaction between the α-hydroxyethyl radical (˙CH(CH3)OH) and bromoacetate (BrCH2CO2-) in the non-buffered (NB), as well as in the bicarbonate (HCO3-) and hydrogen phosphate (HPO42-) buffered aqueous solutions in the presence of ethanol. These complex multistep reactions are initiated by the proton-coupled electron transfer (PCET) which reduces BrCH2CO2- and incites its debromination. The PCET is followed by the step in which the resulting carboxymethyl radical propagates a radical chain reaction thus recovering ˙CH(CH3)OH and enhancing the debromination yields. It is found that the rate constants for the initial PCET step (k1) are raised by ca. an order of magnitude in the presence of the buffers (k1(NB) = 1.4 × 105 dm3 mol-1 s-1; k1(HCO3-) = 1.4 × 106 dm3 mol-1 s-1; k1(HPO42-) = 1.1 × 106 dm3 mol-1 s-1). To rationalize this, we used density functional theory at the M06-2X-D3/6-311+G(2d,p) level in conjunction with the polarizable continuum model (PCM) for an implicit description of the aqueous environment. To acceptably reproduce the measured rate constants, the minimal solute, consisting of ˙CH(CH3)OH, BrCH2CO2- and the buffer anion, has to be expanded by at least 2-3 explicit molecules of the water solvent. The used kinetic model consisting of a set of coupled differential equations indicates the sigmoid dependence of yields vs. k1 thereby confirming the autocatalytic trait of these reactions. The computations unravel the profound influence of the presence of buffers on these reaction systems. On the one hand, the buffer anions promote the PCET by accelerating the proton transfer; on the other hand, they slow down the propagation step by forming the strong hydrogen bonds with the carboxymethyl radical. The two opposing effects cancel out and cause the Br- yields to remain approximately comparable in the non-buffered and buffered media.

Chemical reactivity from the vibrational ground-state level. The role of the tunneling path in the tautomerization of urea and derivatives

Recent developments of low-temperature techniques are providing valuable knowledge about chemical processes that manifest in the quantum regimen. The tunneling effect from the vibrational ground-state is the main mechanism of these reactions, which usually involves the motion or transfer of hydrogen atoms. Theoretical methods can enrich the information supplied by these experimental methods through an insightful analysis of the tunneling process. In this context, canonical variational transition state theory with multidimensional tunneling corrections (CVT/MT) can handle this type of reaction, and it has been applied to several systems within the small-curvature approximation for tunneling (SCT). This method is of proven reliability for polyatomic reactions occurring at room temperature and above, but no tests have been performed to check its performance when only the lowest energy level is populated. In this work, we compare SCT against the least-action tunneling (LAT) method to study the tautomerization and cis-trans interconversion reactions in the enol forms of urea, thiourea, and selenourea. To the best of our knowledge, this is the first time that the LAT method is applied to a polyatomic reaction occurring in the deep-tunneling region. The theoretical results indicate that the reaction mechanisms are controlled by tunneling. The SCT and LAT tautomerization reaction times are in good agreement with the experimental values; however, LAT seems superior to SCT for reactions (tautomerizations) that involve moderate reaction path curvature, whereas the opposite is true for reactions with small curvature (interconversions). These results led us to introduce and recommend the microcanonically optimized tunneling path that selects the tunneling probability as the maximum between the SCT and LAT tunneling probabilities.

Influence of local microenvironment on the double hydrogen transfer in porphycene

We performed time-resolved transient absorption and fluorescence anisotropy measurements in order to study tautomerization of porphycene in rigid polymer matrices at cryogenic temperatures. Studies were carried out in poly(methyl methacrylate) (PMMA), poly(vinyl butyral) (PVB), and poly(vinyl alcohol) (PVA). The results prove that in all studied media hydrogen tunnelling plays a significant role in the double hydrogen transfer which becomes very sensitive to properties of the environment below approx. 150 K. We also demonstrate that there exist two populations of porphycene molecules in rigid media: "hydrogen-transferring" molecules, in which tautomerization occurs on time scales below 1 ns and "frozen" molecules in which double hydrogen transfer is too slow to be monitored with nanosecond techniques. The number of "frozen" molecules increases when the sample is cooled. We explain this effect by interactions of guest molecules with a rigid host matrix which disturbs symmetry of porphycene and hinders tunnelling. Temperature dependence of the number of hydrogen-transferring molecules suggests that the factor which restores the symmetry of the double-minimum potential well in porphycene are intermolecular vibrations localized in separated regions of the amorphous polymer.