RESUMEN
The Python-based program, XMECP, is developed for realizing robust, efficient, and state-of-the-art minimum energy crossing point (MECP) optimization in multiscale complex systems. This article introduces the basic capabilities of the XMECP program by theoretically investigating the MECP mechanism of several example systems including (1) the photosensitization mechanism of benzophenone, (2) photoinduced proton-coupled electron transfer in the cytosine-guanine base pair in DNA, (3) the spin-flip process in oxygen activation catalyzed by an iron-containing 2-oxoglutarate-dependent oxygenase (Fe/2OGX), and (4) the photochemical pathway of flavoprotein adjusted by the intensity of an external electric field. MECPs related to multistate reaction and multistate reactivity in large-scale complex biochemical systems can be well-treated by workflows suggested by the XMECP program. The branching plane updating the MECP optimization algorithm is strongly recommended as it provides derivative coupling vector (DCV) with explicit calculation and can equivalently evaluate contributions from non-QM residues to DCV, which can be nonadiabatic coupling or spin-orbit coupling in different cases. In the discussed QM/MM examples, we also found that the influence on the QM region by DCV can occur through noncovalent interactions and decay with distance. In the example of DNA base pairs, the nonadiabatic coupling occurs across the π-π stacking structure formed in the double-helix system. In contrast to general intuition, in the example of Fe/2OGX, the central ferrous and oxygen part contribute little to the spin-orbit coupling; however, a nearby arginine residue, which is treated by molecular mechanics in the QM/MM method, contributes significantly via two hydrogen bonds formed with α-ketoglutarate (α-KG). This indicates that the arginine residue plays a significant role in oxygen activation, driving the initial triplet state toward the productive quintet state, which is more than the previous knowledge that the arginine residue can bind α-KG at the reaction site by hydrogen bonds.
RESUMEN
The present study investigates the standard model of [1,2]-fluorine migration and that triggered by the rearrangement of cyclopropyl-substituted fluoroepoxides. The [1,2]-fluorine migration reaction proceeds via a synchronous concerted, tight-ion-pair mechanism. When coupled with other reaction coordinates, the whole reaction follows an asynchronous mechanism, while the [1,2]-fluorince migration unit still retains its tight-ion-pair feature and the reaction coordinates of two C-F distances vary synchronously. A general reaction program for α-electron-deficiency-induced [1,2]-fluorine migration is proposed through an analysis of the intermediates generated from nucleophilic addition. The reaction mechanisms associated with α-electron deficiency and the rearrangement are scrutinized using computational chemistry. Two additional reaction programs for [1,2]-fluorine migration are identified. The Gibbs free energy change of [1,2]-fluorine migration exhibits a linear dependence on the value of the Fukui function of the substrate, which could lead to the production of the desired α-monofluoroketone and enhance the utilization of fluorine atoms.