Local Molecular Reactivity of the Colored Dansylglycine in Water and Dioxane Studied through Conceptual DFT

*is study evaluated a fixed long-range corrected range-separated hybrid (RSH) density functional associated with the Def2TZVP basis set alongside the SMD solvation model for the computation of the structure, molecular properties, and chemical reactivity of the M8 intermediate melanoidin pigment in the presence of water and dioxane. *e preference of the active sites pertinent to radical, nucleophilic, and electrophilic attacks is made through linking them with the electrophilic and nucleophilic Parr functions, Fukui function indices, and condensed dual descriptor which are chemical reactivity descriptors that arise from conceptual density functional theory. *e study confirmed the results from previous works showing that the MN12SX density functional is the most appropriate in predicting the chemical reactivity of this molecule in both solvents.


Introduction
Melanoidins are the final product of the Maillard reaction that starts with formation of a Schiff base between a reducing sugar and the amino group of a peptide or protein.Some of these molecules can be isolated, and it has been observed that they possess interesting properties that can be useful not only as dyes for coloring foods and for dye-sensitized solar cells (DSSC) but also for the design and development of drugs for the pharmaceutical industry.
One of these isolated molecules is called dansylglycine (or 2-[ [5-(dimethylamino)naphthalen-1-yl]sulfonylamino] acetic acid) which has interesting properties as a colored molecule and as a fluorescent dye in water and dioxane, and we believe that it could be of interest to study their molecular reactivity by using the ideas of conceptual DFT [1][2][3][4][5], in the same way of our previous works [6][7][8][9][10].
us, in this computational study, we will assess a powerful density functional in calculating the molecular properties and structure of the dansylglycine pigment in the presence of water and dioxane as solvents.Following the same ideas of previous works, we will consider a fixed RSH functional instead of the optimally tuned RSH density functionals that have attained great success [11][12][13][14][15][16][17][18].

Settings and Computational Methods
Following the lines of our previous work [6][7][8][9][10], the computational studies were performed with the Gaussian 09 [34] series of programs with density functional methods as implemented in the computational package.e basis set used in this work was Def2SVP for geometry optimization and frequencies, while Def2TZVP was considered for the calculation of the electronic properties [35,36].All the calculations were performed in the presence of water and dioxane as solvents by doing integral equation formalismpolarized continuum model (IEF-PCM) computations according to the solvation model density (SMD) solvation model [37].

Results and Discussion
In the same way, as in our recent works on melanoidins [6][7][8][9][10], it has been found that the model chemistry formed by the connection between the MN12SX density functional and the Def2TZVP basis set is the best for justifying the fulfilling of the KID procedure for the dansylglycine molecule, both in the presence of water and dioxane.is, instead of presenting the comparison of the values of the orbital energies with the ionization potential I and the electron affinity A for different density functionals, we are showing the results for the global descriptors calculated with the MN12SX density functional in Table 1 comparing the values that arise from both solvents.
It is interesting to observe that the SMD solvation model is powerful enough to allow discriminating between the values for the global descriptors in both solvents.
In the past, various TDDFT studies of molecules of different sizes have used optimally tuned RSH density functionals with great success [11][12][13][14][15][16][17][18].Table 2 provides a comparison between the results involved in the ground-state approximation derived from the HOMO-LUMO gap together with TDDFT results and the experimental value of 340 nm in water and 339 nm in dioxane.It is clearly observed from the results presented in tables that the MN12SX density functional in connection with the Def2TZVP basis set is capable of reproducing accurately the excitation energies of the dansylglycine molecule starting from the HOMO-LUMO gap in both solvents.
Having verified that the MN12SX/Def2TZVP model chemistry is a good choice for the calculation of the global reactivity descriptors and the prediction of the maximum absorption wavelength from the ground state calculation of the HOMO and LUMO, we now present the molecular structures of dansylglycine in water and in dioxane in Figure 1.Meanwhile, the calculated bond lengths and bond angles for both cases are shown in Tables S1-S4 of the Supplementary Information (SI).
e calculations of the condensed Fukui functions and dual descriptor are done by using the Chemcraft molecular analysis program to extract the Mulliken population analysis (MPA), natural population analysis (NPA), and Hirshfeld population analysis (HPA) atomic charges [38] beginning with single-point energy calculations involving the MN12SX density functional that uses the Def2TZVP basis set in line with the SMD solvation model and water and dioxane utilized as solvents.
Considering the potential application the dansylglycine molecule as an antioxidant, it is of interest to get insights into the active sites for radical attack.A graphical representation of the radical Fukui function f 0 calculated with the MN12SX/Def2TZVP model chemistry in water and dioxane is presented in Figure 2.
It can be concluded that the most favorable site for the radical attack will be N7 in both cases, being the chemical reactivity in water a bit larger than in dioxane.
A graphical representation of the dual descriptor Δf k of the dansylglycine molecule calculated with the MN12SX/Def2TZVP model chemistry in water and dioxane is presented in Figure 3.
It can be easily seen from the results for the dual descriptor Δf(r) in Figure 3 that C13 will be the preferred site for a nucleophilic attack and that this atom will act as an electrophilic species in a chemical reaction.In turn, it can be appreciated that N6 will be more prone to electrophilic attacks and that these atomic sites will act as nucleophilic species in chemical reactions that involve the dansylglycine molecule in either solvent.As for the case of radical attacks, in these cases, the chemical reactivity in water will be larger 2 Journal of Chemistry when considering dioxane as a solvent.e results for the condensed descriptor Δf k show the maximum positive value over C13 and the minimum negative value over N6 in both solvents using either population analysis.Finally, the condensed electrophilic and nucleophilic Parr functions P + k and P − k over the atoms of the dansylglycine molecule in water and in dioxane have been calculated by extracting the Mulliken and Hirshfeld (or CM5) atomic charges using the Chemcraft molecular analysis program [38] starting from single-point energy calculations of the ionic species with the MN12SX density functional using the Def2TZVP basis set in the presence of the solvents according to the SMD solvation model.e maximum value of P + k is located over C13 (0.3616 for MPA and 0.2143 for HPA), and the maximum value of P − k is located over N6 (0.4804 for MPA and 0.3719 for HPA) for the calculation in the presence of water, while in dioxane the results are similar with P + k � 0.519 for MPA and 0.2081 for HPA and P − k � 0.3814 for MPA and 0.2962 for HPA.us, there is a nice agreement between the results describing the local reactivity of the dansylglycine molecule provided by Fukui functions, dual descriptor Δf (or its condensed counterpart), and the Parr functions.

Conclusions
A fixed RSH density functional, namely, MN12SX, was examined to determine whether it fulfills the empirical KID procedure.e assessment was conducted by comparing the values from HOMO and LUMO calculations to those generated by the ΔSCF technique for the dansylglycine molecule in water and in dioxane.
is is a compound which is of academic as well as industrial interest.e study has confirmed that the range-separated hybrid meta-NGA density functional MN12SX is the most suited in meeting this goal as it has been observed in previous works.us, this density functional can be a suitable alternative to density functionals where their behavior is optimally tuned using a gap-fitting procedure.It also exhibits the desirable prospect of benefiting future studies aimed at understanding the chemical reactivity of the colored products with larger molecular weights that are obtained when reducing sugars react with proteins and peptides.
From the results of this study, it can be concluded that the model chemistry formed by the connection of the MN12SX density functional and the Def2TZVP basis set is powerful enough to describe adequately the chemical reactivity of the dansylglycine molecule in both solvents through the calculation of the global and local DFT-based reactivity descriptors, including Fukui functions, Parr functions, and the dual descriptor calculations.Furthermore, it is also possible to predict the maximum absorption wavelength for the dansylglycine in water an in dioxane with considerable accuracy.
e prediction involves the MN12SX density functional beginning with the HOMO-LUMO gap instead of TDDFT calculations.Such a finding is particularly crucial considering the likelihood of it being used to inform the alternative determination method on the color that larger systems have such as prosthetic chromophore groups.is becomes necessary in circumstance when it is not possible to carry out TDDFT calculations as in the case of molecules of larger size like peptides and proteins.Journal of Chemistry calculated bond lengths (in Å) of the dansylglycine molecule with the MN12SX density functional using dioxane as a solvent simulated with the SMD solvation model.Table S4: calculated bond angles (in °) of the dansylglycine molecule with the MN12SX density functional using dioxane as a solvent simulated with the SMD solvation model.(Supplementary Materials)

Figure 1 :
Figure 1: A schematic representation of the optimized structure of the dansylglycine molecule in water (a) and in dioxane (b) calculated with the MN12SX density functional showing the numbering of the atoms.

Figure 2 :
Figure 2: A graphical representation of the radical Fukui function f 0 of the dansylglycine molecule calculated with the MN12SX/Def2TZVP model chemistry in water (a) and dioxane (b).

Table 1 :
Global reactivity descriptors for the dansylglycine molecule calculated with the MN12SX/Def2TZVP model chemistry in water and in dioxane.

Table 2 :
Maximum wavelength absorption (λ max ) of the dansylglycine molecule in water and dioxane calculated from the HOMO-LUMO gap and from TDDFT results in comparison with experimental value.