AsiaChem | Chemistry in Japan | December 2021 Volume 2 Issue 1

www.asiachem.news December 2021 | 63 References 1. F. Jensen, Introduction to Computational Chemistry, 3rd ed.; Wiley, Chichester, U.K., 2017. 2. M. Paranjothy, R. Sun, Y. Zhuang, W. L. Hase, WIREs Comput. Mol. Sci. 3, 296 (2013). 3. N. Koga, K. Morokuma, Chem. Rev. 91, 823 (1991). 4. K. N. Houk, P. H.-Y. Cheong, Nature 455, 309 (2008). 5. K. N. Houk, F. Liu, Acc. Chem. Res. 50, 539 (2017). 6. H. B. Schlegel, WIREs Comput. Mol. Sci. 1, 790 (2011). 7. K. Fukui, Acc. Chem. Res. 14, 363 (1981). 8. S. Maeda, Y. Harabuchi, Y. Ono, T. Taketsugu, K. Morokuma, Int. J. Quant. Chem. 115, 258 (2015). 9. S. Maeda, K. Ohno, K. Morokuma, Phys. Chem. Chem. Phys. 15, 3683 (2013). 10. A. L. Dewyer, A. J. Argüelles, P. M. Zimmerman, WIREs Comput. Mol. Sci. 8, e1354 (2018). 11. G. N. Simm, A. C. Vaucher, M. Reiher, J. Phys. Chem. A 123, 385 (2019). 12. S. Maeda, K. Morokuma, J. Chem. Phys. 132, 241102 (2010). 13. S. Maeda, Y. Harabuchi, M. Takagi, T. Taketsugu, K. Morokuma, Chem. Rec. 16, 2232 (2016). 14. M. Hatanaka, T. Yoshimura, S. Maeda, Top. Organomet. Chem. 67, 57 (2020). 15. S. Maeda, Y. Harabuchi, WIREs Comput. Mol. Sci. 11, e1538 (2021). 16. Y. Sumiya, Y. Nagahata, T. Komatsuzaki, T. Taketsugu, S. Maeda, J. Phys. Chem. A 119, 11641 (2015). 17. Y. Sumiya, S. Maeda, Chem. Lett. 49, 553 (2020). 18. Y. Sumiya, S. Maeda, Chem. Lett. 48, 47 (2019). 19. T. Mita, Y. Harabuchi, S. Maeda, Chem. Sci. 11, 7569 (2020). 20. Y. Sumiya, Y. Harabuchi, Y. Nagata, and S. Maeda, ChemRxiv (2021) [DOI: 10.33774/chemrxiv-2021-3xfds]. 21. M. Hada, T. Hasegawa, H. Inoue, M. Takagi, K. Omoto, D. Chujo, S. Iemoto, T. Kuroda, T. Morimoto, T. Hayashi, T. Iijima, T. Tokunaga, N. Ikeda, K. Hujimori, C. Itoh, T. Nishikawa, Y. Yamashita, T. kiwa, S. Koshihara, S. Maeda, Y. Hayashi, ACS Appl. Energy Mater. 2, 7700 (2019). 22. A. Warshel, M. Levitt, J. Mol. Biol. 103, 227 (1976). 23. H. M. Senn, W. Thiel, Angew. Chem. Int. Ed. 48, 1198 (2009). 24. L. W. Chung, W. M. C. Sameera, R. Ramozzi, A. J. Page, M. Hatanaka, G. P. Petrova, T. V. Harris, X. Li, Z. Ke, F. Liu, H. B. Li, L. Ding, K. Morokuma, Chem. Rev. 115, 5678 (2015). 25. T. Vreven, K. Morokuma, Ö. Farkas, H. B. Schlegel, M. J. Frisch, J. Comput. Chem. 24, 760 (2003). 26. T. Vreven, M. J. Frisch, K. N. Kudin, H. B. Schlegel, K. Morokuma, Mol. Phys. 104, 701 (2006). 27. K. Suzuki, K. Morokuma, S. Maeda, J. Comput. Chem. 38, 2213 (2017). 28. K. Suzuki, S. Maeda, K. Morokuma, ACS Omega 4, 1178 (2019). 29. H. Hayashi, H.Takano, H.Katsuyama, Y. Harabuchi, S. Maeda, T. Mita, Chem. Eur. J. 27, 10040 (2021). 30. AFIR-web, see https://afir.sci.hokudai.ac.jp (accessed date: November 16, 2021). of a synthetic method for α,α-difluoroglycine derivatives, which are considered to be bioisosteres of natural glycine and good candidates as drug discovery resources. The flow of this discovery is illustrated in Fig. 5(a). First, the application of QCaRA/AFIR to α,α-difluoroglycine gave a reaction path network, where over 30 reactant candidates were found. Among these candidates, a set of reactants consisting of CF2+NH3+CO2 was selected in consideration of the availability of the species involved. Then, assuming CF3 − and CF 2Br − as sources of CF 2 formation in situ, the reaction path networks for CF3 −+NH 3+CO2 and CF2Br −+NH 3+CO2 were computed by the AFIR method. The results showed that whi le CF3 −+NH 3+CO2 gave CF3CO2 −, CF 2Br −+NH 3+CO2 gave a mixture of the target product and by-products. It was also found that the protons of NH3 promoted the formation of by-products. Then, the reaction path network for CF2Br −+NMe 3+CO2 with the replacement of the proton with the methyl group was obtained by the AFIR method. It was predicted that the desired α,α-difluoroglycine derivatives could be obtained from CF2Br −+NMe 3+CO2 in >99% yield. This article described the reactivity prediction by the AFIR method, which has been developed by the authors, including the latest application results. By applying virtual forces between fragments in a system and inducing chemical changes, the AFIR method gives the reaction paths based on the resulting structural changes. Since the calculated yield of >99% was predicted, the experiment was conducted to confirm the inference. Fig. 5(b) shows the reaction scheme finally discovered experimentally. In the experiment, CF2Br− was generated in situ by mixing Me3SiCF2Br and Ph3SiF2·NBu4. The synthetic experiments afforded α,α-difluoroglycine derivatives from CF2Br −+NMe 3+CO2 in 96% yield. Lastly, we succeeded in isolating the resulting α,α-difluoroglycine derivative as a stable solid after methyl esterification by the treatment with Meerwein reagent. Its three-dimensional structure was confirmed by X-ray crystallography. Furthermore, this three-component reaction proceeded well with various tertiary amines and nitrogen-containing heteroaromatics; hence, it became feasible to synthesize various α,α-difluoroglycine derivatives shown in Fig. 5(c).29 Conclusions This article described the reactivity prediction by the AFIR method, which has been developed by the authors, including the latest application results. By applying virtual forces between fragments in a system and inducing chemical changes, the AFIR method gives the reaction paths based on the resulting structural changes. By systematically testing combinations of various fragment pairs, a systematic automated search for reaction paths is possible. The resulting reaction path network is valuable for elucidating the mechanisms of complicated chemical reactions. The complex reaction path network obtained by the AFIR method can be easily analyzed by the RCMC method based on the kinetics, which is also used to limit the scope of the AFIR search to kinetically accessible paths and to suppress combinatorial explosions. The AFIR method can predict chemical reactions in a system consisting of 30 atoms or less with no previous knowledge. The application to the Strecker synthesis was shown as an example. Even for the system with more than 30 atoms, the reaction mechanism can be analyzed systematically by combining with calculation methods for macromolecular systems such as the semi-empirical quantum chemical calculation method and the QM/ MM method. As examples, this ar ticle presented the structural transition of interfacial amorphous carbon and the conversion of pyruvate to L-lactate by LDH. Lastly, this article introduced the discovery of a synthetic method of difluoroglycine derivatives by QCaRA using the AFIR method as a reaction path search engine. The AFIR method is now widely used in reaction mechanism analysis. This method is implemented in the GRRM20 program and is used by many users.30 The new reaction prediction method combined with QCaRA suggests even more exciting possibilities of the AFIR method. We hope that the new AFIR features in GRRM20 will contribute to many chemical studies in the future. ◆ Acknowledgement This work was partly supported by a grant from JST-CREST (No. JPMJCR14L5), JSTERATO (No. JPMJER1903), and JSPS-WPI. We thank to Ms. Takako Homma for editing a draft of this manuscript.

RkJQdWJsaXNoZXIy NDU2MA==