“…The lowest energy structures were further optimized by undergoing ( Turner et al, 2017 ) (PM7) ( Stewart, 2013 ; Turner et al, 2017 ) in Molecular Orbital PACkage (MOPAC 2016) method. Then we used Schrödinger software to carry out the optimization function applying the Becke 3-parameter (exchange) with correlation by Lee Yang and Parr (B3LYP) ( Becke, 1988 ; Lee et al, 1988 ; Becke, 1993 ; Stephens et al, 1994 ; Adekoya et al, 2022 ; Sakr et al, 2022 ) functional and the 6-31G (**) ( Frisch et al, 1984 ; Adekoya et al, 2022 ; Sakr et al, 2022 ) basis set, density functional theory (DFT) ( Calais, 1993 ; Koch and Holthausen, 2001 ; Adekoya et al, 2022 ; Sakr et al, 2022 ), and the standard Poisson-Boltzmann continuum solvation function (PBF) ( Gilson et al, 1985 ; Sharp and Honig, 1990 ; Tomasi and Persico, 1994 ; Wang et al, 2017 ) for further geometry optimization.…”