Mechanism of Degradation of Rice Starch Amylopectin by Oryzenin Using ONIOM Quantum Calculations [DFT/B3LYP/6-31+G(D, P): AM1]


Understanding the molecular factors of rice degradation during its aging concerns our research team. This article emphasizes oryzenin-amylopectin. It aims to reveal the mechanism of amylopectin deterioration during rice aging. The research exploits the Natural Bond Analysis and ONION method at theory level DFT/B3LYP/6-31+G(d, p) and AM1. This methodological approach allows highlighting amylopectin transformation; oryzenin converts amylopectin into amyloidosis in continuous. This led to monosaccharides and disaccharides.

Share and Cite:

Sawaliho, B. and Robert, N. (2022) Mechanism of Degradation of Rice Starch Amylopectin by Oryzenin Using ONIOM Quantum Calculations [DFT/B3LYP/6-31+G(D, P): AM1]. Computational Chemistry, 10, 139-156. doi: 10.4236/cc.2022.104007.

1. Introduction

Rice is the staple food of the Ivorian people. However, Côte d’Ivoire is struggling to fulfill food self-sufficiency in relation to this product. Post-harvest losses are the main cause. They’re explained by the inefficiency of storage techniques or its alteration during storage. Understanding the molecular factors of this modification concerns our research team. This work follows up on two others on the subject. The first focuses on the roles of phospholipids in this process [1]. The second analyzes the oryzenin-amyloidosis interaction [2]. This article emphasizes oryzenin-amylopectin. It aims to elucidate the role of this latter in the degradation of rice.

Previous work on amyloidosis establishes that its degradation occurs due to its transformation into polysaccharides and monosaccharides [2]. On the contrary, authors [3] explain its alteration by the slow transformation of amylose into amylopectin. They conjecture that this transformation results in a compensatory increase in the second. This debate suggests the following research question:

How does oryzenin interact with amylopectin?

This work aims to understand the mechanism that degrades the components of starch; literature ignores it. Moreover, theoretical studies on amylopectin-oryzenin remain unknown. The research targets to contribute to filling this gap. It assumes that the hydrogen bond (HB) between oryzenin and amylopectin decomposes this latter; its formation impacts the stability of many molecules in the solid state and in solution. It justifies those of molecular structures such as DNA, water [4] [5] and carboxylic acids [6] [7].

HB influences some physical constants such as the melting point temperature of chemical compounds. To answer the research question, this work utilizes the resources of theoretical chemistry. It deploys NBO (Natural Bond Orbital) analysis [8]; it uses charge transfer (CT) [9] [10] [11] to describe the formation of HB as in [2]. Additionally, it exploits the branching model to guide calculations of oryzenin-amylopectin interactions. It relies on an amylopectin subunit; this contains three (AMP3G) or four (AMP4G) glucose molecules linked by sugar bridges. Also, the research employs the resources of quantum chemistry.

This work applies the ONION method at precision [DFT/B3LYP/6-31+G(d, p): AM1] and AM1. Furthermore, amylopectin molecules schematize with three (AMP3G) and four (AMP4G) molecules of α-D glucose or synthons connected by osidic bridges. This methodology highlights the mechanism shaping the degradation of amylopectin under the action of oryzenin. This knowledge, associated with that of the amyloidosis degradation, helps to clarify the mechanism underlying the degradation of rice starch during its prolonged storage. This article consists of three parts.

The first details the results of the Natural Bond Orbital (NBO) analysis. The second focuses on the AMP3G-oryzenin interaction. The third examines the one between the latter and AMP4G. These last two parts present the geometrical, energetic, and spectroscopic parameters of these complexes. Also, it discusses the convergent or divergent aspects between the results obtained with amyloidosis. Before, the article summarizes the hardware and the materials and the method of research.

2. Materials and Method of Research

This work retains the strategy used to study the amylose-oryzenin interaction. The amylopectin subunits comprise three (AMP3G) or four (AMP4G) glucose molecules linked by osidic bridges. The arrangement of glucose molecules considers the branching of amylopectin. The amylopectin-oryzenin complex is divided into two parts. The ONIOM/[DFT/B3LYP/6-31+G(d, p): AM1] offers the possibility of calculating its parameters with different levels of precision [12] [13]. The active part concerns the interaction site between the hydrogen of oryzenin and the oxygen of amylopectin. This part undergoes a calculation with high precision DFT/(B3LYP/6-31+G(d, p)). On the contrary, the remaining part of the complex is treated with the low precision semi-empirical AM1 method. Moreover, the research uses Gaussian 09 software to optimize all the parameters of the complexes [14].

Vibration frequency calculations help validate their geometric parameters. From an ideal structure, a “single point” calculation at the same level of theory offers the opportunity to perform the NBO calculation [15]. As in the amylose-oryzenin complexes, the interactions considered are those established between osidic oxygen of amylopectin (by its subunits) and oryzenin-hydrogen. The numbering of the atoms is generated automatically by the GaussView06 software. The electrostatic potential [16] [17] [18] [19] at the level of the B3LYP/ 6-31+G(d, p) theory makes it possible to identify the HB donor or acceptor sites of oryzenin [2]. More, this article presents and discusses its results.

3. Results and Discussion

This section presents the results related to the interactions between the two polysaccharides of the study and oryzenin. It specifies those linking to the NBO analysis. Table 1 and Table 2 present the stabilization energies E(2), the differences

Table 1. Stabilization energy of second-order E(2) perturbations in the AMP3G-oryzenin complex.

Table 2. Stabilization energy of second-order E(2) perturbations in the AMP4G-oryzenin complex.

between the energies of NBO electron donors (i) and NBO electron acceptor (j). They also show the elements F(i, j) of the Fock matrix and the CT. The model of the latter, constructed by Reed and al [11] , makes it possible to describe the interactions in the AMP3G-oryzenin and AMP4G-oryzenin complexes.

3.1. NBO Analysis of the Amylopectin (AMP3G/AMP4G)-Oryzenin Complexes

The NBO analysis presents the results for which the stabilization energies are greater than 1.5 kcal·mol−1. CT are correlated with these quantities. The relative positions of the interacting elements influence them. In the AMP3G-oryzenin complex, for the oryzenin-H30…O49 and oryzenin-H30…O99 interactions, the most important contributions to equilibrium result from the transfer of free chlorine pairs from oryzenin. They correspond to the transitions:

- For the interaction oryzenin-H30…O49:

n Cl 36 ( 4 ) σ C 69 -H 77 * with E(2) = 51.82 kcal·mol−1 and CT = 117.37 me

- For oryzenin-H30…O99:

n Cl 36 ( 4 ) σ O 91 -H 102 * with E(2) = 38.5 kcal·mol−1 and CT = 76.10 me.

A secondary stabilization exists; it’s based on the following transitions:

- For oryzenin-H30…O49,

n O 49 ( 2 ) σ O 29 -H 30 * with E(2) = 8.46 kcal·mol−1 and CT = 16.73 me

- For oryzenin-H30…O99,

n O 99 ( 2 ) σ O 29 -H 30 * with E (2) = 10.59 kcal·mol−1 and CT = 20.98 me

The CT of 20.98 suggests to me that the strong HB in AMP3G-oryzenin comes from the oryzenin-H30…O99 interaction. For those with the N9H10 of oryzenin, no CT is observed with chlorine. Its position prohibits it within the complexes. The most important transfers are those concerning O49 and N9-H10 bonds. They’re

- For the interaction N9-H10…O49,

n O 49 ( 1 ) σ N 9 -H 10 * with E(2) = 5.96 kcal·mol−1 and CT = 9.24 me

- For the interaction N9-H10…O99,

n O 99 ( 1 ) σ N 9 -H 10 * with E(2) = 8.59 kcal·mol−1 and CT = 17.55 me.

In the case of the AMP4G-oryzenin complex, the strongest stabilization comes from free chlorine pairs:

- For the interaction oryzenin-H30…O49,

n Cl 36 ( 4 ) σ C 69 -H 77 * with E(2) = 39.47 kcal·mol−1 and CT = 79.01 me

- For the interaction oryzenin-H30…O91,

n Cl 36 ( 4 ) σ O 119 -H 120 * transition with E(2) = 33.79 kcal·mol−1 CT = 71.90 me

- For oryzenin-H30…O99, two transitions of chlorine bring the strongest stabilization. It’s about:

n Cl 36 ( 4 ) σ C 93 -H 94 * with E(2) = 20.82 kcal·mol−1 and CT = 41.73 me

n Cl 36 ( 4 ) σ O 119 -H 120 * with E(2) = 21.81 kcal·mol−1 and CT = 44.62 me

Besides, the stabilization brought mainly by chlorine, there’s a secondary one.

This result from transitions:

- For oryzenin-H30…O49,

n O 49 ( 1 ) σ O 29 -H 30 * with E(2) = 8.61 kcal·mol−1 and CT = 13.30 me

n O 49 ( 2 ) σ O 29 -H 30 * with E(2) = 6.32 kcal·mol−1 and CT = 12.88 me

- For oryzenin-H30…O91,

n O 91 ( 1 ) σ O 29 -H 30 * with E(2) = 9.99 kcal·mol−1 and CT = 14.83 me.

- For oryzenin-H30…O99,

n O 91 ( 1 ) σ O 29 -H 30 * with E(2) = 10.14 kcal·mol−1 and CT = 16.53 me.

Oxygen O99 participates more in CT. The corresponding value is 16. This proves that the stabilization is maximum for the transition which concerns it. Under these conditions, the oryzenin H30…O99 interaction becomes the most probable. Regarding N9H10…O49, no stabilization comes from chlorine. CT results from transition:

n O 49 ( 2 ) σ N 9 -H 10 * with E(2) = 5.33 kcal·mol−1 and CT = 11.25 me.

- For N9H10…O91, two transitions from the chlorine atom provide the solidest stabilization. They are

n Cl 36 ( 4 ) σ O 54 -H 123 * with E(2) = 28.98 kcal·mol−1 and CT = 59.91 me

n Cl 36 ( 4 ) σ O 51 -H 52 * with E(2) = 17.94 kcal·mol−1 and CT = 36.22 me

- For the interaction N9-H10…O99:

n O 99 ( 1 ) σ N 9 -H 10 * with E(2) = 7.24 kcal·mol−1 and CT = 11.40 me.

HB is strong with the oryzenin-H30 probe associated with O99 in the AMP3G and in AMP4G complexes. This interaction results from the lone chlorine pairs of oryzenin. It initiates the transformation of amylopectin at its branching level. The research exploits this result to construct the ideal complexes. It uses it to model the interactions between oryzenin and AMP4G or AMP3G. In previous research [2] , the electrostatic potential of molecular interaction sites indicated that the oryzenin H30 and H10 constitute the highest potential VS,max. This work appropriates this result to study AMP3G.

3.2. Interaction AMP3G-Oryzenine

These atoms illustrate the hydrogen donor potential of the oryzenin. The geometrical, energetic, spectroscopic parameters and the NBO analysis relate to the interactions involving this hydrogen. Figure 1 shows the geometries of the optimized complexes. A red ball represents an oxygen atom. A big white corresponds to a carbon. A small one schematizes hydrogen.


Figure 1. 3D structures of optimized AMP4G-oryzenin. (a) Interaction N9H10…O49; (b) Interaction N9H10…O99; (c) Interaction N9H10…O99; (d) Interaction oryzenin-H30…O99.

Table 3 shows the geometric parameters. It presents the distance d, the angles of linearity α and of direction β associated with HB. The distance d varies from 1.92 Å to 2.17 Å [20]. It characterizes HB. They’re less than 2.62 Å [21] [22]; this value defines the Vander Waals sum radii of the hydrogen and oxygen atoms [23] [24]. Oryzenin-H30…O49, oryzenin-H30…O99 and N9H10…O99 respect the condition for HB [24]. On the other hand, the α and β associated with the N9H10…O49 interaction deviate from the ideal values. These geometrical parameters establish that oryzenin-H30…O99 is the strongest HB.

Table 4 shows the variations of the enthalpy, the entropy and the free enthalpy of reaction linked to the interactions between AMP3G and oryzenin. Those relating to the enthalpies of reaction are all negative. They fluctuate between −2.730 kcal·mol−1 and −20.175 kcal·mol−1. These values suggest that all reactions are exothermic. Those associated with the oryzenin-H30…O99 interaction represent the weakest.

This interaction corresponds to the most stable. Its free enthalpy of reaction (∆rG = −5.30 kcal·mol−1) indicates that its formation is spontaneous. The other interactions (N9H10…O49 and N9H10…O99) are exothermic but not spontaneous (∆rG > 0). Figure 2 shows free enthalpies of reaction changes associated with AM3G-oryzenin and AMP3G-oryzenin complexes. It illustrates that the others (N9H10…O49 and N9H10…O99) aren’t spontaneous (ArG > 0).

These results suggest the degradation modality of amylopectin; this process begins with the establishment of HB oryzenin-H30…O99. It breaks the osidic bridges of the ramification. It leads to the formation of amyloidosis. Its production linked to this process and its natural disappearance [2] occurs concomitantly. They explain the slow decrease in its rate presence in the starch during the aging of the rice. Furthermore, spectroscopic parameters precise this finding.

Table 3. Geometric parameters of the AMP3G-oryzenine complex at ONIOM (B3LYP/6-31+G(d, p): AM1).

Table 4. Energetic parameters of AMP3G-oryzenine complex at ONIOM (B3LYP/6-31+G(d, p): AM1).

Figure 2. Free enthalpies of reaction changes associated with AM3G-oryzenine and AMP3G-oryzenin complex.

Table 5 outlines the O-H and N-H stretching frequency shift of free and complex oryzenin. It presents their differences. These reach 281 cm−1 for oryzenin-H30…O99 interaction and 236 cm−1 for oryzenin-H30…O49. Those of the N-H bonds remain lower than them. They only attain 126 cm−1 for HB N9-H10…O49. Important variation in vibrational frequencies reflects the strong oxygen attraction of the O99 osidic bridge around H30. It’s greater for the O-H bond; it reaches 281 cm−1 for oryzenin-H30…O99. It reveals that H30 exerts its strongest attraction on oxygen O99 of the amylopectin branching osidic bridge. The latter breaks. This rupture leads to linear chains of α-D-glucose or amyloidosis. This result suggests that the degradation of amylopectin is easier than that of amylose.

Table 5. Elongation frequencies and its variations (cm1) associated with O-H and N-H of free and complex oryzenin in the AMP4G-oryzenin interaction.

The rupture of its osidic bridge conducts to the formation of this starch latter component. This result refutes the conclusions of [3]. Its authors conjecture that the storage of rice favours the fall in the proportion of amylose in the starch; they affirm that an increase in that of amylopectin in the same proportion simultaneously. Thus, degradation of rice starch changes the amylose-amylopectin ratio. This result corroborates the observations of Cao and al [25] regarding the variation over time of this latter. The first degradation of amylopectin also supports the conclusions drawn from the analysis of geometric, energetic, and spectroscopy. Moreover, this research concerns the AMP4G-oryzenin complex.

3.3. Interaction AMP4G-Oryzenin

This section presents and discusses the results related to the geometrical, energetic, spectroscopic parameters and NBO analysis of the AMP4G-oryzenin complexes. The optimized geometries (Figure 3) show three osidic bridges (O49, O91, O99) of the four α-D-glucose “synthons” representing the portion of amylopectin interacting with the two main oryzenin interaction sites (oryzenin-H30 and N9-H10). The distance d materialized by a dotted line in each complex indicates the interaction between the two sites. The geometrical parameters allow identifying it.

Table 6 shows the geometric parameters (d, α, β) of the interactions characterizing HB in the AMP4G-oryzenin complex. Its data agree with those established by Desiraju and al [22]. They indicate that the oryzenin-H30…O49 and oryzenin-H30…O91 interactions lead to strong HB; their distances d are very close (1.93 Å; 1.94 Å).

The environment around the HB influences the angular parameters. The linearity angle α of the oryzenin-H30…O91 interaction equals 170.8˚. It reaches 172.4˚ for oryzenin-H30…O49. The angles α deviate from their ideal values, 180˚. β, associated oryzenin-H30…O49 interaction, is 5.1˚ lower than the expected 109.9˚. The β angle of the oryzenin-H30…O99 interaction exceeds it by 0.6. It’s at the upper limit of the interval of Desiraju and al [22]. These two angles β are different from 109.9˚. Analysis of the three geometric parameters suggests that HB H30…O91 is the most stable of the AMP4G-oryzenin complex.


Figure 3. 3D structures of optimized AMP4G-oryzenin. (a) oryzenin-H30…O49 interaction; (b) oryzenin-N9H10…O49 interaction; (c) oryzenin-H30…O91 interaction; (d) Oryzenin N9H10…O91 interaction; (e) oryzenin-H30…O99 interaction; (f) Interaction N9H10…O99.

Table 6. Geometrical parameters of the AMP4G-oryzenin complex at ONIOM (B3LYP/6-31+G(d, p): AM1).

Its distance d respects the criterion of Desiraju and al [22]. Its angle α and β remain close to their ideal values. These parameters lead to the conclusion that the AMP4G degradation begins in O91 locate on osidic bridge position α (1, 6).

Table 7 compiles the variation of the enthalpies, entropy, and free enthalpies for AMP4G-oryzenin interactions. The variations of the enthalpies are all negative. They appear between −3.525 kcal·mol−1 and −27.968 kcal·mol−1; The HB formation is exothermic. The ∆rH also suggest that the interactions created by oryzenin-H30 and N9H10 are all stable. The lowest ∆rH corresponds to the interaction oryzenin-H30…O99 (∆rH = −27.968 kcal·mol−1).

The negative ∆rG values with oryzenin-H30 show that the interactions oryzenin-H30…O49, oryzenin-H30…O91 and oryzenin-H30…O99 are spontaneous while those established with N9H10 aren’t. ∆rG coincides with the bottommost value (∆rG = −9.719 kcal·mol−1) for the interaction oryzenin-H30…O99. This latter is the strongest HB and the most spontaneous. The HB with O49 and O91 aren’t negligible. These results indicate that the O99 site is the most favourable for HB formation. This suggests that AMP4G degrade primarily from this site.

Table 8 shows the υ(O-H) and υ(N-H) stretching frequencies of free and associated oryzenin. The stretch frequency changes are all positive. The strengths of HB linked with them decrease in the order O99, O91 and O49 for AMP4G. They correspond to the most stable interactions between AMP4G and oryzenin-H30; the most important variation is those related to the oryzenin-H30…O99 interaction of 261 cm−1. It proves that the latter is the strongest of this group.

The O99 site is the most favourable for the formation of HB. This suggests that AMP4G mostly breaks down from there. Its deterioration begins with this osidic bridge in position α (1.6). The rupture of the latter promotes the formation of linear α-D-glucose chains; it leads to amyloidosis. These consist of tetra saccharides. They degrade according to the mechanism described by [2]. Then, this work elucidates the mechanism of starch degradation following the appearance of amyloidosis.

Figure 4 shows the free enthalpy changes associated with interactions between O29H30 or N9H10 and the appropriate proton acceptor sites of AM4G and AMP4G. Its data reveal that all variations linked to AM4G and AMP4G are negative for the O29H30 probe. The modifications related to the interactions between N9H10 and AM4G are also. The formation of HB following oryzenin-O29H30-AMP4G or

Table 7. Energetic parameters of AMP4G-oryzenin complexes at ONIOM (B3LYP/6-31+G(d, p): AM1).

Table 8. Elongation frequencies and changes in O-H and N-H elongation frequencies (cm−1) free and complex oryzenin from AMP4G-oryzenin interactions.

Figure 4. AM4G-oryzenin and AMP4G-oryzenin free enthalpy variations.

oryzenin-N9H10-AM4G interactions correspond to spontaneous processes for the osidic bridges O49, O91, O99 polysaccharides. On the other hand, the positive values associated with the free enthalpy demonstrate that the interactions of oryzenin-N9H10-AMP4G aren’t. This energetic component hardly sheds light on the complete mechanism underlying the degradation of AMP4G.

The interactions between two polysaccharides with O29H30 lead spontaneously to the formation of HB with each osidic bridge O49, O91. O99. In other words, their thermodynamic similarity doesn’t make it possible to detect the second stage of the AMP4G disintegration during its progressive transformation into AM4G. This limit justifies the choice of spectroscopy. This can explain degradation order of the two saccharides in this phase deterioration; the local nature of the phenomenon pleads for this approach.

For the same oryzenin-H30 probe, the two oxygen osidic bridges O(1) and O(2) possible for the oryzenin-H30…O91 interaction, the difference in frequency between the free and associated forms is 257 cm−1. Moreover, that between AM4G and AMP4G is −30 cm−1 for O(1) and +7 cm−1 around O(2) as the receptor site. The frequency variation is greater for amylopectin near O(1) compared to amyloidosis. It reflects a stronger attraction of oxygen from the osidic bridge oxygen O91 towards H30. It breaks the latter more easily.

The degradation of amylopectin into amylose is carried out continuously when the first interact with O91. The same reasoning leads to similar conclusions with O99 and O49. In the latter case, the transformation of AMP4G takes place around O(1) mainly; a smaller proportion comes from its O(2). Furthermore, this research extends to NBO analysis of the AMP4G-oryzenin and AMP3G-oryzenin complexes.

4. Conclusions

This work aims to clarify the chemical processes underlying the degradation of amylopectin with oryzenin. It exploits resources from quantum chemistry as a method. It uses NBO analysis and calculates the geometric, energetic, and spectroscopic parameters of the AMP3G-oryzenin and AMP4G-oryzenin complexes. It employs the ONIOM [DFT/B3LYP/6-31+G(d, p): AM1]. This method makes it possible to suggest a mechanism to explain the degradation of amylopectin by oryzenin. This latter establishes HB through H30 or H10.

This work locates the presence of HB near the oxygen of AMP3G or AMP4G situated at α (1, 6) on the osidic bridge O99. The latter coincides with their branching points. HB results from the strong attraction of oxygen O99 to H30. It causes the rupture of this bridge. It leads to the formation of amyloid. On the other hand, the frequency variation is greater for amylopectin around O(1); it exceeds that associated with the interaction with amyloidosis near O(2). This difference reflects a stronger attraction of oxygen O99, O91 or O49 towards H30. It suggests that the degradation of amylopectin is easier than that of amylose. So, these forces readily break the branching of amylopectin. Its degradation into amyloid becomes continuous. The latter deteriorates gradually according to the mechanism described by [2]; its products are monosaccharides and disaccharides during rice aging or storage. This finding contributes to the debate on rice starch degradation; it refutes the conclusions of [3].

According to [3] , the storage of rice promotes a reduction in the proportion of amylose in the starch; it causes an equivalent increase in that of amylopectin. On the other hand, the observations of Caoet al. [25] confirm the findings of this research. These authors state that amylopectin is degraded before amyloidosis. Moreover, research reinforces the team’s findings [2]. It specifies that the deterioration of rice starch begins with that of its amylopectin. Therefore, it explains the origin of disaccharides and monosaccharides during the aging of rice. Their presence is linked to two types of amyloidosis. Some of it comes from the continual breakdown of amylopectin. The other relates to amyloid, a constituent of starch. Moreover, the article highlights the main difference between amylopectin and amyloid regarding action of oryzenin.

The attraction of oxygen from the saccharide bridge by H30 differs between amylopectin-oryzenin and amyloid-oryzenin. The O-H stretching frequency difference associated with it exceeds that of the second complex. Simply, the attraction of the oryzenin’s H30 for amylopectin’s oxygen on the osidic bridge remains higher than that observed for amyloid. This result constitutes further evidence that amylopectin degrades before amyloid in rice starch.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.


[1] Koffi, K.A., N’Guessan, R.B. and Bamba, E.H.S. (2021) The Phospholipid Degradation in Paddy Rice: A Theoretical Model with DFT/B3LYP 6-311 G. European Journal of Applied Sciences, 9, 162-174.
[2] N’guessan, B.R., Bamba, E.H.S. and Koffi, K.A. (2022) The Oryzenin’s Effect on Di, Tri and Quadri-Saccharide Degradation. An Investigation by a Mixed Method: ONIOM (DFT/B3LYP/6-31 + G(d, p): AM1). Computational Chemistry, 10, 97-119.
[3] Abeysundara, A., Navaratne, S., Wickramasinghe, I. and Ekanayake, D. (2017) Determination of Changes of Amylose and Amylopectin Content of Paddy during Early Storage. International Journal of Science and Research, 6, 2094-2097.
[4] Cipcigan, F., Sokhan, V., Martyna, G. and Crain, J. (2018) Structure and Hydrogen Bonding at the Limits of Liquid Water Stability. Scientific Reports, 8, Article No. 1718.
[5] Liu, J., He, X., Zhang, J.Z.H. and Qi, L.-W. (2018) Hydrogen-Bond Structure Dynamics in Bulk Water: Insights from ab Initio Simulations with Coupled Cluster Theory. Chemical Science, 9, 2065-2073.
[6] Arno, B., Layton, M.L. and Kathleen, K.V. (2000) Hydrogen Bonds in Carboxylic Acid-Carboxylate Systems in Solution. 1. In Anhydrous, Aprotic Media. Organic Letters, 2, 2007-2009.
[7] Zhang, Q. and Du, L. (2016) Hydrogen Bonding in the Carboxylic Acid-Aldehyde Complexes. Computational and Theoretical Chemistry, 1078, 123-128.
[8] Weinhold, F. (2012) Natural Bond Orbital Analysis: A Critical Overview of Relationships to Alternative Bonding Perspectives. Journal of Computational Chemistry, 33, 2363-2379.
[9] Behzadi, H., Esrafili, M.D. and Hadipour, N.L. (2007) A Theoretical Study of 17O, 14N and 2H Nuclear Quadrupole Coupling Tensors in the Real Crystalline Structure of Acetaminophen. Chemical Physics, 333, 97-104.
[10] Esrafili, M.D., Behzadi, H. and Hadipour, N.L. (2008) 14N and 17O Electric Field Gradient Tensors in Benzamide Clusters: Theoretical Evidence for Cooperative and Electronic Delocalization Effects in N-H…O Hydrogen Bonding. Chemical Physics, 348, 175-180.
[11] Reed, A.E., Curtiss, L.A. and Weinhold, F. (1988) Intermolecular Interactions from a Natural Bond Orbital, Donor-Acceptor Viewpoint. Chemical Reviews, 88, 899-926.
[12] Dapprich, S., Komáromi, I., Byun, K.S., Morokuma, K. and Frisch, M.J. (1999) A New ONIOM Implementation in Gaussian98. Part I. The Calculation of Energies, Gradients, Vibrational Frequencies and Electric Field Derivatives. Journal of Molecular Structure: THEOCHEM, 461-462, 1-21.
[13] Vreven, T., Byun, K.S., Komáromi, I., Dapprich, S., Montgomery, J.A., Morokuma, K. and Frisch, M.J. (2006) Combining Quantum Mechanics Methods with Molecular Mechanics Methods in ONIOM. Journal of Chemical Theory and Computation, 2, 815-826.
[14] Frisch, M.J., Trucks, G.W., Schlegel, H.B., et al. (2009) Gaussian 09, Revision A.02. Gaussian, Inc., Wallingford.
[15] Günay, N., Pir, H., Avcı, D. and Atalay, Y. (2013) NLO and NBO Analysis of Sarcosine-Maleic Acid by Using HF and B3LYP Calculations. Journal of Chemistry, 2013, Article ID: 712130.
[16] Murray, J.S. and Sen, K. (1996) Molecular Electrostatic Potentials: Concepts and Applications. Elsevier, Amsterdam, 105-538.
[17] Brinck, T. (1998) The Use of the Electrostatic Potential for Analysis and Prediction of Intermolecular Interactions. Theoretical and Computational Chemistry, 5, 51-93.
[18] Demircioğlu, Z., Albayrak, Ç. and Büyükgüngör, O. (2014) Theoretical and Experimental Investigation of (E)-2-([3,4-dimethylphenyl)imino]methyl)-3-methoxyphenol: Enolketo Tautomerism, Spectroscopic Properties, NLO, NBO and NPA Analysis. Journal of Molecular Structure, 1065-1066, 210-222.
[19] Bader, R.F.W., Carroll, M.T., Cheeseman, J.R. and Chang, C. (1987) Properties of Atoms in Molecules: Atomic Volumes. Journal of the American Chemical Society, 109, 7968-7979.
[20] Arunan, E., Desiraju, G.R., Klein, R.A., Sadlej, J., Scheiner, S., Alkorta, I., Clary, D.C., Crabtree, R.H. and Dannenberg, J.J. (2011) Definition of the Hydrogen Bond (IUPAC Recommendations 2011). Pure and Applied Chemistry, 83, 1637-1641.
[21] Joseph, J. and Jemmis, E.D. (2007) Red-, Blue-, or No-Shift in Hydrogen Bonds: A Unified Explanation. Journal of the American Chemical Society, 129, 4620-4632.
[22] Desiraju, G. and Steiner, T. (2001) The Weak Hydrogen Bond: In Structural Chemistry and Biology. Oxford University Press, Oxford, 480.
[23] Bondi, A. (1964) Van der Waals Volumes and Radii. The Journal of Physical Chemistry, 68, 441-451.
[24] Rowland, R.S. and Taylor, R. (1996) Intermolecular Nonbonded Contact Distances in Organic Crystal Structures: Comparison with Distances Expected from van der Waals Radii. The Journal of Physical Chemistry, 100, 7384-7391.
[25] Cao, Y., Wang, Y., Chen, X. and Ye, J. (2004) Study on Sugar Profile of Rice during Ageing by Capillary Electrophoresis with Electrochemical Detection. Food Chemistry, 86, 131-136.

Copyright © 2022 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.