All the simulations were performed using the Gromacs 2018.8 software in NPT ensemble with T = 300 K maintained through V-rescale algorithm 59 with τ T = 1.0 ps and P = 1 bar using Berendsen isotropic barostat with τ P = 0.5 ps. A cut-off radius of 1.1 nm was used for