Theoretical Analysis of Coordination Geometries in Transition Metal-Histidine Complexes Using Quantum Chemical Calculations

A theoretical investigation utilizing density functional theory (DFT) calculations was conducted to explore the coordination complexes formed between histidine (His) ligands and various divalent transition metal ions (Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+). Conformational exploration of the His ligand was initially performed to assess its stability upon coordination. Both 1:1 and 1:2 of metal-to-ligand complexes were scrutinized to elucidate their structural features and the relative stability of the complexes. This study examined the ability of His to act as a bidentate or tridentate coordinating ligand, along with the differences in coordination geometry when solvent effects were incorporated. The reduced density gradient (RDG) analysis and local electron attachment energy (LEAE) analysis were employed to elucidate the interaction planes and the nucleophilic and electrophilic properties. The electronic properties were analyzed through electrostatic potential (ESP) maps and natural population analysis (NPA) of atomic charge distributions. This computational study provides valuable insights into the diverse coordination modes of His and its interactions with divalent transition metal ions, contributing to a better understanding of the role of this amino acid ligand in the formation of transition metal complexes. The findings can aid in the design and construction of self-assembled structures involving His-metal coordination.

Quantum Chemical Investigation into the Structural Analysis and Calculated Raman Spectra of Amylose Modeled with Linked Glucose Molecules

This study presents a quantum chemical investigation into the structural analysis and calculated Raman spectra of modeled amylose with varying units of linked glucose molecules. We systematically examined the rotation of hydroxymethyl groups and intramolecular hydrogen bonds within these amylose models. Our study found that as the number of linked glucose units increases, the linear structure becomes more complex, resulting in curled, cyclic, or helical structures facilitated by establishing various intramolecular interactions. The hydroxymethyl groups were confirmed to form interactions with oxygen atoms and with hydroxymethyl and hydroxyl groups from adjacent rings in the molecular structures. We identified distinct peaks and selected specific bands applicable in various analytical contexts by comparing their calculated Raman spectra. Representative vibrational modes within selected regions were identified across the different lengths of amylose models, serving as characteristic signatures for linear and more coiled structural conformations. Our findings contribute to a deeper understanding of amylose structures and spectroscopic signatures, with implications for theoretical studies and potential applications. This work provides valuable reference points for the detailed assignment of Raman peaks of amylose structure, facilitating their application in broader research on carbohydrate structures and their associated spectroscopic properties.

Exploring the Quantum Chemical Energy Landscape with GNN-Guided Artificial Force

Artificial force has been proven useful to get over energy barriers and quickly search a large portion of the energy landscape. This work proposes a method based on graph neural networks to optimize the choice of transformation patterns to examine and accelerate energy landscape exploration. In open search from glutathione, the search efficiency was largely improved in comparison to random selection. We also applied transfer learning from glutathione to tuftsin, resulting in further efficiency gains.

Dynamically Hidden Reaction Paths in the Reaction of CF3+ + CO

Reaction paths on a potential energy surface are widely used in quantum chemical studies of chemical reactions. The recently developed global reaction route mapping (GRRM) strategy automatically constructs a reaction route map, which provides a complete picture of the reaction. Here, we thoroughly investigate the correspondence between the reaction route map and the actual chemical reaction dynamics for the CF₃⁺ + CO reaction studied by guided ion beam tandem mass spectrometry (GIBMS). In our experiments, FCO⁺, CF₂⁺, and CF⁺ product ions were observed, whereas if the collision partner is N₂, only CF₂⁺ is observed. Interestingly, for reaction with CO, GRRM-predicted reaction paths leading to the CF⁺ + F₂CO product channel are found at a barrier height of about 2.5 eV, whereas the experimentally obtained threshold for CF⁺ formation was 7.48 ± 0.15 eV. In other words, the ion was not obviously observed in the GIBMS experiment, unless a much higher collision energy than the requisite energy threshold was provided. On-the-fly molecular dynamics simulations revealed a mechanism that hides these reaction paths, in which a non-statistical energy distribution at the first collisionally reached transition state prevents the reaction from proceeding along some reaction paths. Our results highlight the existence of dynamically hidden reaction paths that may be inaccessible in experiments at specific energies and hence the importance of reaction dynamics in controlling the destinations of chemical reactions.

Effect of alkali metal cations on network rearrangement in polyisoprene ionomers

The effects of cations, Li⁺, Na⁺, and Cs⁺, on the structure of ionic aggregates and network rearrangement in carboxylated polyisoprene (PI) ionomers were studied. We found that network rearrangement via interaggregate hopping of metal carboxylates is improved with a decrease in cation size, even though density functional theory (DFT) calculation indicated the increase in the attractive interaction between metal carboxylates. At the same time, we also found that as the size of the cation decreases, the inclusion of the PI segment in the ionic aggregate increases. The DFT calculation suggested the cation-π interaction between the cation and double bonds in the PI segment as the cause for the inclusion. The inclusion of the PI segment with a low glass transition temperature (T_g) plasticizes the ionic aggregate and would sterically hinder the attractive interaction between metal carboxylates. In fact, the electron spin resonance measurement revealed a decrease in the T_g of the ionic aggregate with a decrease in cation size. Based on our findings, we proposed that the inclusion of PI segments in the ionic aggregate is the possible cause for the enhancement of network rearrangement in the carboxylated PI ionomers with a decrease in the cation size.

Reaction Space Projector (ReSPer) for Visualizing Dynamic Reaction Routes Based on Reduced-Dimension Space

To analyze chemical reaction dynamics based on a reaction path network, we have developed the "Reaction Space Projector" (ReSPer) method with the aid of the dimensionality reduction method. This program has two functions: the construction of a reduced-dimensionality reaction space from a molecular structure dataset, and the projection of dynamic trajectories into the low-dimensional reaction space. In this paper, we apply ReSPer to isomerization and bifurcation reactions of the Au₅ cluster and succeed in analyzing dynamic reaction routes involved in multiple elementary reaction processes, constructing complicated networks (called "closed islands") of nuclear permutation-inversion (NPI) isomerization reactions, and elucidating dynamic behaviors in bifurcation reactions with reference to bundles of trajectories. Interestingly, in the second application, we find a correspondence between the contribution ratios in the ability to visualize and the symmetry of the morphology of closed islands. In addition, the third application suggests the existence of boundaries that determine the selectivity in bifurcation reactions, which was discussed in the phase space. The ReSPer program is a versatile and robust tool to clarify dynamic reaction mechanisms based on the reduced-dimensionality reaction space without prior knowledge of target reactions.

Stereoselective Transesterification of P-Chirogenic Hydroxybinaphthyl Phosphinates

The substitution reaction of phosphinates with a binaphthyloxy group at the phosphorus atom with lithium alkoxides proceeded with good to high efficiencies to give P-chirogenic phosphinates with a high enantiomeric ratio. As alcohols, primary, secondary, and tertiary alcohols could be used, and the use of tert-butyl alcohol yielded the products with a higher enantiomeric ratio. A substrate with two different alkyl groups on the phosphorus atom could also participate in the substitution reaction to give the corresponding products in good yields with excellent selectivity. The molecular structures of one of the substrates and the corresponding products, determined by X-ray analyses, proved that the substitution reaction at the phosphorus atom proceeded with inversion of the absolute configuration. The usefulness of the reaction was demonstrated by using it to prepare a drug candidate for Duchenne muscular dystrophy. Finally, thionation of the resulting phosphinates was carried out to form P-chirogenic phosphinothioates.

Visualization of reaction route map and dynamical trajectory in reduced dimension

In the quantum chemical approach, chemical reaction mechanisms are investigated based on a potential energy surface (PES). Automated reaction path search methods enable us to construct a global reaction route map containing multiple reaction paths corresponding to a series of elementary reaction processes. The on-the-fly molecular dynamics (MD) method provides a classical trajectory exploring the full-dimensional PES based on electronic structure calculations. We have developed two reaction analysis methods, the on-the-fly trajectory mapping method and the reaction space projector (ReSPer) method, by introducing a structural similarity to a pair of geometric structures and revealed dynamic aspects affecting chemical reaction mechanisms. In this review, we will present the details of these analysis methods and discuss the dynamics effects of reaction path curvature and reaction path bifurcation with applications to the CH₃OH + OH^- collision reaction and the Au₅ cluster branching and isomerization reactions.

Geometric analysis of anharmonic downward distortion following paths

A mathematical aspect of the anharmonic downward distortion following (ADDF) path is discussed. The ADDF method is utilized as an automated reaction path search method, which can explore transition state geometries on a potential energy surface from a potential minimum. We show that the maximum number of the ADD stationary paths intersecting the potential minimum is 2^f + 1  - 2, where f denotes the degree of freedom of the system. We also show that the bifurcation of the ADD stationary path is essential to detect all the transition states connected from a given minimum. The ADDF computation is demonstrated for a H₂ O molecule in which pitchfork bifurcation is observed.

Reductive Activation of Small Molecules by Anionic Coinage Metal Atoms and Clusters in the Gas Phase

Metal atoms and clusters exhibit chemical properties that are significantly different or totally absent in comparison to their bulk counterparts. Such peculiarity makes them potential building units for the generation of novel catalysts. Investigations of the gas-phase reactions between size- and charge-selected atoms/clusters and small molecules have provided fundamental insights into their intrinsic reactivity, thus leading to a guiding principle for the rational design of the single-atom and cluster-based catalysts. Especially, recent gas-phase studies have elucidated that small molecules such as O₂ , CO₂ , and CH₃ I can be catalytically activated by negatively-charged atoms/clusters via donation of a partial electronic charge. This Minireview showcases typical examples of such "reductive activation" processes promoted by anions of metal atoms and clusters. Here, we focus on anionic atoms/clusters of coinage metals (Cu, Ag, and Au) owing to the simplicity of their electronic structures. The determination of a correlation between their activation modes and the electronic structures might be helpful for the future development of innovative coinage metal catalysts.

A Theoretical Study on the Mechanism of the Oxidative Deborylation/C-C Coupling Reaction of Borepin Derivatives

One-electron oxidation of borepin derivatives that consists of a boron-containing seven-membered ring has been reported to cause deborylation/C-C coupling, yielding aromatic compounds. The reaction can be achieved not only by transition metal compounds but also by oxidants without transition metal such as O₂ and other organic compounds. Despite numerous experimental attempts, the mechanism of this peculiar reaction as well as the fate of the BCl part eliminated from borepin remain unclear to date. Based on theoretical approaches using the artificial force induced reaction method, here we address the mechanism of the unusual boron-mediated C-C coupling. For this purpose, two borepin derivatives (1 and 35), bearing ethyl and phenyl groups, respectively, were used as reactants, and FeCl₃/MeNO₂ and O₂ were chosen as oxidants. The calculations revealed reaction pathways that provided an overall picture of the mechanism of the target reaction, which features four key steps, namely, (i) quaternization of the boron atom by the coordination of oxidant, (ii) intersystem crossing, (iii) skeletal rearrangement to form a six-membered ring, and (iv) elimination of a boron moiety. The intrinsic nature of boron, i.e., a strong tendency to accept a coordination ligand even under oxidative conditions, is responsible for the oxidative deborylation/C-C coupling of borepin.

Understanding CO oxidation on the Pt(111) surface based on a reaction route network

Kinetic analysis by the rate constant matrix contraction on the reaction route network of CO oxidation on the Pt(111) surface obtained by the artificial force induced reaction reveals the impact of entropic contributions arising from a variety of local minima and transition states.

Analyses of trajectory on-the-fly based on the global reaction route map

A methodology to analyze a trajectory on-the-fly (TOF) based on a global reaction route map consisting of intrinsic reaction coordinate (IRC) pathways is proposed. By using the distance functions in the configurational space, the location of each point on the trajectories is detected, providing a dynamical picture that the molecular system goes over several minima and transition states in the reaction path network. In its application to structural transformations of an Au₅ cluster, a variety of reaction routes are obtained, and the hopping from one IRC to another IRC (IRC-jump) is analyzed. The branching of trajectories over many minima on the potential energy surface via valley-ridge transition points is also discussed.

Implementation and performance of the artificial force induced reaction method in the GRRM17 program

This article reports implementation and performance of the artificial force induced reaction (AFIR) method in the upcoming 2017 version of GRRM program (GRRM17). The AFIR method, which is one of automated reaction path search methods, induces geometrical deformations in a system by pushing or pulling fragments defined in the system by an artificial force. In GRRM17, three different algorithms, that is, multicomponent algorithm (MC-AFIR), single-component algorithm (SC-AFIR), and double-sphere algorithm (DS-AFIR), are available, where the MC-AFIR was the only algorithm which has been available in the previous 2014 version. The MC-AFIR does automated sampling of reaction pathways between two or more reactant molecules. The SC-AFIR performs automated generation of global or semiglobal reaction path network. The DS-AFIR finds a single path between given two structures. Exploration of minimum energy structures within the hypersurface in which two different electronic states degenerate, and an interface with the quantum mechanics/molecular mechanics method, are also described. A code termed SAFIRE will also be available, as a visualization software for complicated reaction path networks. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Propargyl-Assisted Selective Amidation Applied in C-terminal Glycine Peptide Conjugation

Alkyl esters, such as propargyl esters, typically lack the electron-withdrawing inductive effects needed to participate in nucleophilic acyl substitution reactions. Herein, we report an unusual observation in which glycine propargyl ester derivatives displayed selective, base-independent reactivity towards linear alkylamines under mild, metal-free conditions. Through global reaction route mapping (GRRM) modeling calculations, it is predicted that these observations may be governed by factors related to hydrogen-bonding and intermolecular interactions, rather than electron-withdrawing inductive effects. Based on this concept of propargyl-assisted selective amidation, a direct application was made to develop a novel site-specific C-terminal glycine peptide bioconjugation technique as a proof-of-concept, which relies upon the selective reactivity of glycine propargyl esters over that of aspartate and glutamate side-chain-linked propargyl esters.

Study of Potential Energy Surfaces towards Global Reaction Route Mapping

The potential energy surface (PES) is just a theoretical construct based on the Born-Oppenheimer approximation, but it underlies various phenomena, including molecular vibrations, collisional ionizations, and chemical reactions. This account describes how a new idea for global reaction route mapping (GRRM), which had seemed to be impossible for chemical systems with more than three atoms, was born and has been developed during the course of the study of the PES. GRRM has pioneered new fields of chemistry. Furthermore, techniques for GRRM are still developing, and GRRM is further extending its application to various areas of chemistry and chemical physics.

Exploring the role of a single water molecule in the tropospheric reaction of glycolaldehyde with an OH radical: a mechanistic and kinetics study

This work reveals that though a single-water molecule decelerates the atmospheric reaction between the glycolaldehyde and OH radical, however, it facilitates the cis–​trans interconversion along the hydrogen-abstraction pathways.

The effect of Mg2+ incorporation on the structure of calcium carbonate clusters: investigation by the anharmonic downward distortion following method

A structure similar to calcite appears when only four CaCO₃ units aggregate into the cluster, and the addition of Mg changes its structure.

Exploring the mechanism of isomerisation and water-migration in the water-complexes of amino-acid l-proline: electrostatic potential and vibrational analysis

This work reveals interesting pathways for water-migration and neutral ↔ zwitterionic isomerisation in the water complexes of l-proline.