TY - JOUR
T1 - AMADAR: a python-based package for large scale prediction of Diels–Alder transition state geometries and IRC path analysis
AU - Isamura, Bienfait
AU - Lobb, Kevin Alan
N1 - Funding Information:
Bienfait K. Isamura is grateful to the Department of Chemistry at Rhodes University and the Center for High Performance Computing (CHPC) for having provided the resources used to carry out this study (project CHEM0802).
Funding Information:
This work was supported by the BEBUC Scholarship program through the funding granted by the Else-Kroener-Fresenius Stiftung.
Publisher Copyright:
© 2022, The Author(s).
PY - 2022/6/15
Y1 - 2022/6/15
N2 - Predicting transition state geometries is one of the most challenging tasks in computational chemistry, which often requires expert-based knowledge and permanent human intervention. This short communication reports technical details and preliminary results of a python-based tool (AMADAR) designed to generate any Diels–Alder (DA) transition state geometry (TS) and analyze determined IRC paths in a (quasi-)automated fashion, given the product SMILES. Two modules of the package are devoted to performing, from IRC paths, reaction force analyses (RFA) and atomic (fragment) decompositions of the reaction force F and reaction force constant κ. The performance of the protocol has been assessed using a dataset of 2000 DA cycloadducts retrieved from the ZINC database. The sequential location of the corresponding TSs was achieved with a success rate of 95%. RFA plots confirmed the reaction force constant κ to be a good indicator of the (non)synchronicity of the associated DA reactions. Moreover, the atomic decomposition of κ allows for the rationalization of the (a)synchronicity of each DA reaction in terms of contributions stemming from pairs of interacting atoms. The source code of the AMADAR tool is available on GitHub [CMCDD/AMADAR(github.com)] and can be used directly with minor customizations, mostly regarding the local working environment of the user
AB - Predicting transition state geometries is one of the most challenging tasks in computational chemistry, which often requires expert-based knowledge and permanent human intervention. This short communication reports technical details and preliminary results of a python-based tool (AMADAR) designed to generate any Diels–Alder (DA) transition state geometry (TS) and analyze determined IRC paths in a (quasi-)automated fashion, given the product SMILES. Two modules of the package are devoted to performing, from IRC paths, reaction force analyses (RFA) and atomic (fragment) decompositions of the reaction force F and reaction force constant κ. The performance of the protocol has been assessed using a dataset of 2000 DA cycloadducts retrieved from the ZINC database. The sequential location of the corresponding TSs was achieved with a success rate of 95%. RFA plots confirmed the reaction force constant κ to be a good indicator of the (non)synchronicity of the associated DA reactions. Moreover, the atomic decomposition of κ allows for the rationalization of the (a)synchronicity of each DA reaction in terms of contributions stemming from pairs of interacting atoms. The source code of the AMADAR tool is available on GitHub [CMCDD/AMADAR(github.com)] and can be used directly with minor customizations, mostly regarding the local working environment of the user
KW - Diels–Alder reaction
KW - Finite difference approach
KW - Hellman–Feynman force
KW - RDKit
KW - Reaction force analysis
U2 - 10.1186/s13321-022-00618-3
DO - 10.1186/s13321-022-00618-3
M3 - Article
C2 - 35706060
SN - 1758-2946
VL - 14
SP - 39
JO - Journal of Cheminformatics
JF - Journal of Cheminformatics
IS - 1
M1 - 39
ER -