Next Article in Journal
Theoretical Study of Intramolecular Interactions in Peri-Substituted Naphthalenes: Chalcogen and Hydrogen Bonds
Previous Article in Journal
Regulation of T Cell Activation and Differentiation by Extracellular Vesicles and Their Pathogenic Role in Systemic Lupus Erythematosus and Multiple Sclerosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chemical Reactivity Theory Study of Advanced Glycation Endproduct Inhibitors

by
Juan Frau
1 and
Daniel Glossman-Mitnik
1,2,*
1
Departament de Química, Universitat de les Illes Balears, 07122 Palma de Mallorca, Spain
2
Laboratorio Virtual NANOCOSMOS, Centro de Investigación en Materiales Avanzados, Departamento de Medio Ambiente y Energía, Chihuahua, Chih 31136, Mexico
*
Author to whom correspondence should be addressed.
Molecules 2017, 22(2), 226; https://doi.org/10.3390/molecules22020226
Submission received: 14 December 2016 / Accepted: 25 January 2017 / Published: 2 February 2017
(This article belongs to the Section Computational and Theoretical Chemistry)

Abstract

:
Several compounds with the known ability to perform as inhibitors of advanced glycation endproducts (AGE) have been studied with Density Functional Theory (DFT) through the use of a number of density functionals whose accuracy has been tested across a broad spectrum of databases in Chemistry and Physics. The chemical reactivity descriptors for these systems have been calculated through Conceptual DFT in an attempt to relate their intrinsic chemical reactivity with the ability to inhibit the action of glycating carbonyl compounds on amino acids and proteins. This knowledge could be useful in the design and development of new drugs which can be potential medicines for diabetes and Alzheimer’s disease.

Graphical Abstract

1. Introduction

It is well established that advanced glycation endproducts (AGEs) are the final products of a series of chemical reactions initiated by the attachment of reducing sugars to free amino groups in proteins, lipids, and nucleic acids [1,2]. AGEs are a group of complex and heterogeneous compounds that are implicated in diabetic complications and Alzheimer’s disease.
Therefore, there is a large amount of research dedicated to the discovery and investigation of AGE inhibitors that can be used as a therapeutic approach for the treatment of the pathologies associated with them. AGE inhibitors may act by various mechanisms at different steps of advanced glycation endproduct (AGE) formation (depending on oxidative stress and/or carbonyl stress) and AGE-mediated damage: trapping of carbonyl species, antioxidant activity by transition metal chelation, or free radical scavenging.
Chemical Reactivity Theory or Conceptual Density Functional Theory (DFT) [3,4] provides a series of descriptors that can be useful for the analysis of the molecular properties of known AGE inhibitors, and a starting point for the development of new ones. These descriptors are based on the calculation of the energies of the neutral system as well as the corresponding radical cation and anion by means of a ΔSCF technique, that is, the difference in energies between the neutral and the radical ionic species. This can be made simpler by applying the so-called KID procedure (for “Koopmans in DFT”) that consists of the identification of the vertical ionization potential I with ϵ H (the energy of the HOMO) and the vertical electron affinity A with ϵ L (the energy of the LUMO). This KID procedure is an approximation, because it is well known that the Koopmans’s theorem is not valid within DFT [5,6,7,8]. Notwithstanding, it can be useful for faster calculations of Conceptual DFT descriptors for large molecular systems where the determination of the electronic energy of the radical cation and anion could be computationally costly and difficult to converge.
As a follow-up of our recently published work on the chemistry of antioxidant molecules [9], the objective of this work is to perform a comparative study of the performance of several model chemistries based on the latest Minnesota density functionals [10] for the prediction of the chemical reactivity of twelve inhibitors of advanced glycation endproducts presented in the literature [11,12]. The comparison will be performed through the calculation of the global descriptors that arise from Conceptual DFT and the verification of the KID procedure. It must be stressed that it is not our intention to calculate the vertical I and A to compare them to experimental results, but to predict trends of chemical reactivity that can be useful for other researchers in order to develop new AGE inhibitors.

2. Theoretical Background

As this work is part of an ongoing project, the theoretical background is similar to that presented in previous research [9,13,14,15,16,17,18], and will be shown here for the sake of completeness. Within the conceptual framework of DFT [19,20], the chemical potential μ is defined as
μ = E N v ( r ) = χ
where χ is the electronegativity, while the global hardness η is
η = 2 E N 2 v ( r )
Using a finite difference approximation and the KID procedure, the above expressions can be written as:
μ = 1 2 ( I + A ) 1 2 ( ϵ L + ϵ H ) = χ K
η = ( I A ) ( ϵ L ϵ H ) = η K
where ϵ H and ϵ L are the energies of the highest occupied and the lowest unoccupied molecular orbitals (HOMO and LUMO), respectively.
The electrophilicity index ω has been defined as [21]:
ω = μ 2 2 η = ( I + A ) 2 4 ( I A ) ( ϵ L + ϵ H ) 2 4 ( ϵ L ϵ H ) = ω K
The electrodonating ( ω ) and electroaccepting ( ω + ) powers have been defined as [22]:
ω = ( 3 I + A ) 2 16 ( I A ) ( 3 ϵ H + ϵ L ) 2 16 η K = ω K
ω + = ( I + 3 A ) 2 16 ( I A ) ( ϵ H + 3 ϵ L ) 2 16 η K = ω K +
It follows that a larger ω + value corresponds to a better capability of accepting charge, whereas a smaller value of ω makes it a better electron donor. In order to compare ω + with - ω , the following definition of net electrophilicity has been proposed [23]:
Δ ω ± = ω + ( ω ) = ω + + ω ω K + ( ω K ) = ω K + + ω K = Δ ω K ±
that is, the electroaccepting power relative to the electrodonating power.

3. Settings and Computational Methods

Along the lines of our previous work [9,13,14,15,16,17,18], the computational studies were performed with the Gaussian 09 [24] series of programs, with density functional methods as implemented in the computational package. The equilibrium geometries of the molecules were determined by means of the gradient technique. The force constants and vibrational frequencies were determined by computing analytical frequencies on the stationary points obtained after the optimization to check if there were true minima. The basis set used in this work was Def2SVP for geometry optimization and frequencies, while Def2TZVP was considered for the calculation of the electronic properties [25,26].
For the calculation of the molecular structure and properties of the studied systems [9,13,14,15,16,17,18], we chose several density functionals from the latest Minnesota density functionals family, which consistently provide satisfactory results for several structural and thermodynamic properties [10]: M11, which is a is a range-separated hybrid meta-GGA [27]; M11L, which is a dual-range local meta-GGA [28]; MN12L, which is a nonseparable local meta-NGA [29]; MN12SX, which is a range-separated hybrid nonseparable meta-NGA [30]; N12, which is a nonseparable gradient approximation [31]; N12SX, which is a range-separated hybrid nonseparable gradient approximation [30]; SOGGA11, which is a GGA density functional [32]; and SOGGA11X, which is a hybrid GGA density functional [33]. In these functionals, GGA stands for generalized gradient approximation (in which the density functional depends on the up and down spin densities and their reduced gradient), and NGA stands for nonseparable gradient approximation (in which the density functional depends on the up/down spin densities and their reduced gradient, and also adopts a nonseparable form). All calculations were performed in the presence of water as a solvent, by doing IEF-PCM (Integral Equation Formalism for the Polarized Continuum Model) computations according to the Solvation Model Density (SMD) solvation model [34].

4. Results and Discussion

The molecular structures of ALT-711, ALT-946, aminoguanidine, carnosine, GLY-230, LR-9, metformin, OPB-9195, pentoxifylline, pioglitazone, pyridoxamine, and tenilsetam [11,12] were pre-optimized by starting with the readily available structures shown in Figure 1 and finding the most stable conformers by means of the Avogadro 1.2.0 program [35,36] through a random sampling with molecular mechanics techniques and a consideration of all the torsional angles through the general AMBER force field [37]. The structures of the resulting conformers were then reoptimized with the eight density functionals mentioned in the previous section in conjunction with the Def2SVP basis set and the SMD solvation model, using water as a solvent.
It is worth calculating the electronegativity χ, the global hardness η, and the global electrophilicity ω for the studied systems using the ΔSCF and KID in order to verify the quality of the procedures. Additionally, we will include the electrodonating ( ω ) and electroaccepting ( ω + ) powers as well as the net electrophilicity Δ ω ± in the calculations for further verification.
The HOMO and LUMO orbital energies (in eV), ionization potentials I, electron affinities A (in eV), and global electronegativity χ, total hardness η, global electrophilicity ω, electrodonating power, ( ω ), electroaccepting power ( ω + ), and net electrophilicity Δ ω ± of the AGE inhibitors calculated with the eight density functionals and the Def2TZVP basis set using water as solvent simulated with the SMD parametrization of the IEF-PCM model are presented in Tables S1A to S8A of the Electronic Supplementary Materials (ESM). The upper part of the tables shows the results derived assuming the validity of the KID procedure (hence the subscript K), and the lower part shows the results derived from the calculated ΔSCF energies.
We have previously designed several descriptors that relate the results obtained through the HOMO and LUMO calculations with those obtained by means of the vertical I and A with a ΔSCF procedure. However, it must be stressed that it is not our intention to perform a gap-fitting by minimizing a descriptor by choosing optimal range-separation parameter γ, but to check if the density functionals considered in this study—in which, some of them contain a fixed range-separation parameter γ—obey the KID procedure. In fact, there is no range-separation parameter γ in our designed descriptors. Moreover, we have considered A as minus the energy of the LUMO of the neutral system instead of considering A as minus the energy of the HOMO of the N+1 electron system, as it was in References [38,39].
The first three descriptors are related to the simplest fulfillment of the KID procedure by relating ϵ H with -I, ϵ L with -A, and the behavior of them in the description of the HOMO–LUMO gap: J I = | ϵ H + E g s ( N 1 ) E g s ( N ) | , J A = | ϵ L + E g s ( N ) E g s ( N + 1 ) | , and J H L = J I 2 + J A 2 .
Next, we consider four other descriptors that analyze how well the studied density functionals are useful for the prediction of the electronegativity χ, the global hardness η, and the global electrophilicity ω, and for a combination of these Conceptual DFT descriptors, just considering the energies of the HOMO and LUMO or the vertical I and A: J χ = | χ χ K | , J η = | η η K | , J ω = | ω ω K | and J D 1 = J χ 2 + J η 2 + J ω 2 , where D1 stands for the first group of Conceptual DFT descriptors.
Finally, we designed four other descriptors to verify the goodness of the studied density functionals for the prediction of the electroaccepting power ω + , the electrodonating power ω , the net electrophilicity Δ ω ± , and for a combination of these Conceptual DFT descriptors, just considering the energies of the HOMO and LUMO or the vertical I and A: J ω + = | ω + ω K + | , J ω = | ω ω K | , J Δ ω ± = | Δ ω ± Δ ω K ± | and J D 2 = J ω 2 + J ω + 2 + J Δ ω ± 2 , where D2 stands for the first group of Conceptual DFT descriptors.
The results of the calculations of J I , J A , J H L , J χ , J η , J ω , J D 1 , J ω + , J ω , J Δ ω ± , and J D 2 for the AGE inhibitors considered in this work are displayed in Tables S1B to S8B of the Electronic Supplementary Materials.
On the basis of the results for the descriptors presented on Tables S1B–S8B of the ESM, we have compiled the average values for for each density functional on the whole group of AGE inhibitors, and the calculated results are displayed on Table 1.
As can be seen from Table 1, the KID procedure holds with great accuracy for the MN12SX and N12SX density functionals, which are a range-separated hybrid meta-NGA and a range-separated hybrid NGA density functionals, respectively. Indeed, the values of J I , J A , and J H L are not exactly zero. However, their values can be favorably compared with the results presented for these quantities in the work of Lima et al. [39], where the minima have been obtained by choosing a parameter that enforces that behavior.
It is interesting to see that the same density functionals also fulfill the KID procedure for the other descriptors, namely J χ , J η , J ω , and J D 1 , as well as for J ω , J ω + , J Δ ω ± , and J D 2 . These results are very important, because they show that it is not enough to rely only on J I , J A , and J H L . For example, if we consider only J χ , the values are very close to zero for all of the density functionals considered. As for the other descriptors, only the MN12SX and N12SX density functionals show this behavior. This means that the results for J χ are due to a fortuitous cancellation of errors.
The usual GGA (SOGGA11) and hybrid-GGA (SOGGA11X) are not good for the fulfillment of the KID procedure, and the same conclusion is valid for the local functionals M11L, MN12L and N12.
An important fact is that although the range-separated hybrid NGA and range-separated hybrid meta-NGA density functionals can be useful for the calculation of the Conceptual DFT descriptors, it is not the same for the range-separated hybrid GG A (M11) density functional. An inspection of Table S1A of the ESM shows that this is due to the fact that this functional inadequately describes the energy of the LUMO, leading to negative values of A, which are in contradiction with the ΔSCF results.
As the mode of action of the studied AGE inhibitors can be different, involving one, two or three of the mechanisms presented in the Introduction, it is very difficult to say which of them will be more effective as an inhibitor of the glycation process. However, in a recent experimental work on the reactivity, selectivity, and reaction mechanisms of some sequestering agents of reactive carbonyl species [40], the authors attempted to relate their results with several Conceptual DFT descriptors. Although only three of our studied AGEs inhibitors were considered in that work, it is worth seeing if a similar comparison can be performed using the results of our calculations. Notwithstanding, it is not fair to attempt a linear regression analysis with only three points, but potentially good qualitative trends for their ability to perform as inhibitors of AGEs could be probably obtained.
Indeed, the key factors in the study of the reactivity of the considered AGEs inhibitors is their nucleophilic character, because we are only assuming the interaction through the Maillard reaction. Thus, we are going to consider the published results for the interaction of the sequestering agents with methylglyoxal (MGO). There are several definitions of nucleophilicity available in the literature of Conceptual DFT, and the interested reader is referred to the recent work of Domingo and Pérez [41]. For our purposes, the most significant definitions are those that relate the nucleophilicity N with the inverse of the electrophilicity ω [42] or with the inverse of the electrodonating power ω [43]. It is worth observing that the trend for the ability to react with glycating carbonyls that arises from Table 1 of the work of Colzani et al. [40]—that is, aminoguanidine > carnosine > pyridoxamine—correlates well with the previous definitions of nucleophilicity N, as calculated from the results of Tables S4A and S6A of the ESM. It is interesting to see that the same correlations are found if we consider the inverse of the net electrophilicity Δ ω ± results that are shown in Tables S4A and S6A of the ESM. This suggests that the inverse of Δ ω ± could be an alternative definition for the nucleophilicity N.
On the basis of the previous analysis, and extending the correlation to the other molecular systems under study, we can find the approximate qualitative trends:
ALT 946 > Aminoguanidine > Metformin > Carnosine > GLY 230 > Tenilsetam >
> Pyridoxamine > Pentoxifylline Pioglitazone > LR 9 > OPB 9195 > ALT 711
for the MN12SX density functional, and
Aminoguanidine > ALT 946 > Metformin > Carnosine > GLY 230 > Tenilsetam >
> Pyridoxamine > Pentoxifylline Pioglitazone > LR 9 > OPB 9195 > ALT 711
for the N12SX density functional.
These qualitative trends are very similar (except for the different rank of aminoguanidine in each case), and they are representative of the known pharmacological properties of the studied AGE inhibitors [11,12]. Therefore, the calculation of the nucleophilicity N (as the inverse of ω, the inverse of ω , or probably the inverse of Δ ω ± ) within both model chemistries could be useful for the prediction of the AGE inhibition ability of potential molecular systems designed with this purpose in mind.

5. Conclusions

The latest Minnesota family of density functionals (M11, M11L, MN12L, MN12SX, N12, N12SX, SOGGA11, and SOGGA11X) have been tested for the fulfillment of the KID procedure by comparison of the HOMO- and LUMO-derived values with those obtained through a ΔSCF technique. It has been shown that the range-separated hybrid meta-NGA density functional (MN12SX) and the range-separated hybrid NGA density functional (N12SX) are the best for the accomplishment of this objective. As such, they are a good alternative to those density functionals whose behavior have been tuned through a gap-fitting procedure and a good prospect for their usefulness in the description of the chemical reactivity of molecular systems of large size.
From the whole of the results presented in this work, it can be seen that the sites of interaction of the studied AGE inhibitors can be predicted by using DFT-based reactivity descriptors such as the electronegativity, global hardness, global electrophilicity, electrodonating and electroaccepting powers, and net electrophilicity. These descriptors were used in the characterization and successful description of the proposed AGE inhibitors, and provide a firm explanation for the reactivity of those molecules.
Moreover, a qualitative trend has been found between the potential pharmacological activity and the nucleophilicity N of the studied molecules. This is based on calculations performed with the MN12SX and N12SX density functionals in connection with the Def2TZVP basis set and the SMD parametrization of the IEF-PCM model using water as a solvent. It can be concluded that these model chemistries (MN12SX/Def2TZVP/SMD(Water) and N12SX/Def2TZVP/SMD(Water)) are the best for fulfilling the KID procedure and for the prediction of the of the potential AGE inhibition ability of molecular systems.

Supplementary Materials

The following are available online at: https://www.mdpi.com/1420-3049/22/2/226/s1, ESM.pdf : Electronic Supplementary Materials.

Acknowledgments

This work has been partially supported by CIMAV, SC and Consejo Nacional de Ciencia y Tecnología (CONACYT, Mexico) through Grant 219566/2014 for Basic Science Research and Grant 265217/2016 for a Foreign Sabbatical Leave. Daniel Glossman-Mitnik conducted this work while a Sabbatical Fellow at the University of the Balearic Islands from which support is gratefully acknowledged. This work was cofunded by the Ministerio de Economía y Competitividad (MINECO) and the European Fund for Regional Development (FEDER) (CTQ2014-55835-R).

Author Contributions

Daniel Glossman-Mitnik conceived and designed the research and headed, wrote and revised the manuscript, while Juan Frau contributed to the writing and the revision of the article. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DFTDensity Functional Theory
AGEAdvanced Glycation Endproducts
SCFSelf Consistent Field
KIDKoopmans in DFT
HOMOHigher Ocuppied Molecular Orbital
LUMOLower Unocuppied Molecular Orbital
IEF-PCMIntegral Equation Formalism - Polarized Continuum Model
SMDSolvation Model Density

References

  1. Singh, R.; Barden, A.; Mori, T.; Beilin, L. Advanced Glycation End-Products: A Review. Diabetologia 2001, 44, 129–146. [Google Scholar] [CrossRef] [PubMed]
  2. Peng, X.; Ma, J.; Chen, F.; Wang, M. Naturally Occurring Inhibitors Against the Formation of Advanced Glycation End-Products. Food Funct. 2011, 2, 289–301. [Google Scholar] [CrossRef] [PubMed]
  3. Parr, R.; Yang, W. Density-Functional Theory of Atoms and Molecules; Oxford University Press: New York, NY, USA, 1989. [Google Scholar]
  4. Chattaraj, P. (Ed.) Chemical Reactivity Theory—A Density Functional View; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2009.
  5. Lewars, E. Computational Chemistry—Introduction to the Theory and Applications of Molecular and Quantum Mechanics; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2003. [Google Scholar]
  6. Young, D. Computational Chemistry—A Practical Guide for Applying Techniques to Real-World Problems; John Wiley & Sons: New York, NY, USA, 2001. [Google Scholar]
  7. Jensen, F. Introduction to Computational Chemistry, 2nd ed.; John Wiley & Sons: Chichester, UK, 2007. [Google Scholar]
  8. Cramer, C. Essentials of Computational Chemistry—Theories and Models, 2nd ed.; John Wiley & Sons: Chichester, UK, 2004. [Google Scholar]
  9. Frau, J.; Muñoz, F.; Glossman-Mitnik, D. A Molecular Electron Density Theory Study of the Chemical Reactivity of Cis- and Trans-Resveratrol. Molecules 2016, 21, 1650. [Google Scholar] [CrossRef] [PubMed]
  10. Peverati, R.; Truhlar, D.G. Quest for a Universal Density Functional: The Accuracy of Density Functionals Across a Broad Spectrum of Databases in Chemistry and Physics. Philos. Trans. Ser. A Math. Phys. Eng. Sci. 2014, 372. [Google Scholar] [CrossRef] [PubMed]
  11. Rahbar, S.; Figarola, J.L. Novel Inhibitors of Advanced Glycation Endproducts. Arch. Biochem. Biophys. 2003, 419, 63–79. [Google Scholar] [CrossRef] [PubMed]
  12. Peyroux, J.; Sternberg, M. Advanced Glycation Endproducts (AGEs): Pharmacological Inhibition in Diabetes. Pathol. Biol. 2006, 54, 405–419. [Google Scholar] [CrossRef] [PubMed]
  13. Glossman-Mitnik, D. A Comparison of the Chemical Reactivity of Naringenin Calculated with M06 Family of Density Functionals. Chem. Cent. J. 2013, 7, 155–161. [Google Scholar] [CrossRef] [PubMed]
  14. Martínez-Araya, J.I.; Salgado-Morán, G.; Glossman-Mitnik, D. Computational Nanochemistry Report on the Oxicams—Conceptual DFT and Chemical Reactivity. J. Phys. Chem. B 2013, 117, 6639–6651. [Google Scholar] [CrossRef] [PubMed]
  15. Glossman-Mitnik, D. Computational Nanochemistry Study of the Chemical Reactivity Properties of the Rhodamine B Molecule. Proc. Comput. Sci. 2013, 18, 816–825. [Google Scholar] [CrossRef]
  16. Martínez-Araya, J.I.; Salgado-Morán, G.; Glossman-Mitnik, D. Computational Nutraceutics: Chemical Reactivity Properties of the Flavonoid Naringin by Means of Conceptual DFT. J. Chem. 2013, 2013, 850270. [Google Scholar] [CrossRef]
  17. Glossman-Mitnik, D. Chemical Reactivity Theory within DFT Applied to the Study of the Prunin Flavonoid. Eur. Int. J. Sci. Techno. 2014, 3, 195–207. [Google Scholar]
  18. Glossman-Mitnik, D. Computational Chemistry of Natural Products: A Comparison of the Chemical Reactivity of Isonaringin Calculated with the M06 Family of Density Functionals. J. Mol. Model. 2014, 20, 1–7. [Google Scholar] [CrossRef] [PubMed]
  19. Parr, R.; Yang, W. Density Functional Approach to the Frontier-Electron Theory of Chemical Reactivity. J. Am. Chem. Soc. 1984, 106, 4049–4050. [Google Scholar] [CrossRef]
  20. Geerlings, P.; de Proft, F.; Langenaeker, W. Conceptual Density Functional Theory. Chem. Rev. 2003, 103, 1793–1873. [Google Scholar] [CrossRef] [PubMed]
  21. Parr, R.; Szentpaly, L.; Liu, S. Electrophilicity Index. J. Am. Chem. Soc. 1999, 121, 1922–1924. [Google Scholar] [CrossRef]
  22. Gázquez, J.; Cedillo, A.; Vela, A. Electrodonating and Electroaccepting Powers. J. Phys. Chem. A 2007, 111, 1966–1970. [Google Scholar] [CrossRef] [PubMed]
  23. Chattaraj, P.; Chakraborty, A.; Giri, S. Net Electrophilicity. J. Phys. Chem. A 2009, 113, 10068–10074. [Google Scholar] [CrossRef] [PubMed]
  24. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09 Revision D.01; Gaussian Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  25. Weigend, F.; Ahlrichs, R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef] [PubMed]
  26. Weigend, F. Accurate Coulomb-fitting Basis Sets for H to R. Phys. Chem. Chem. Phys. 2006, 8, 1057–1065. [Google Scholar] [CrossRef] [PubMed]
  27. Peverati, R.; Truhlar, D.G. Improving the Accuracy of Hybrid Meta-GGA Density Functionals by Range Separation. J. Phys. Chem. Lett. 2011, 2, 2810–2817. [Google Scholar] [CrossRef]
  28. Peverati, R.; Truhlar, D.G. M11-L: A Local Density Functional That Provides Improved Accuracy for Electronic Structure Calculations in Chemistry and Physics. J. Phys. Chem. Lett. 2012, 3, 117–124. [Google Scholar] [CrossRef]
  29. Peverati, R.; Truhlar, D.G. An Improved and Broadly Accurate Local Approximation to the Exchange-Correlation Density Functional: the MN12-L Functional for Electronic Structure Calculations in Chemistry and Physics. Phys. Chem. Chem. Phys. 2012, 14, 13171–13174. [Google Scholar] [CrossRef] [PubMed]
  30. Peverati, R.; Truhlar, D.G. Screened-Exchange Density Functionals with Broad Accuracy for Chemistry and Solid-State Physics. Phys. Chem. Chem. Phys. 2012, 14, 16187–16191. [Google Scholar] [CrossRef] [PubMed]
  31. Peverati, R.; Truhlar, D.G. Exchange-Correlation Functional with Good Accuracy for Both Structural and Energetic Properties while Depending Only on the Density and Its Gradient. J. Chem. Theory Comput. 2012, 8, 2310–2319. [Google Scholar] [CrossRef] [PubMed]
  32. Peverati, R.; Zhao, Y.; Truhlar, D.G. Generalized Gradient Approximation That Recovers the Second-Order Density-Gradient Expansion with Optimized Across-the-Board Performance. J. Phys. Chem. Lett. 2011, 2, 1991–1997. [Google Scholar] [CrossRef]
  33. Peverati, R.; Truhlar, D.G. Communication: A Global Hybrid Generalized Gradient Approximation to the Exchange-Correlation Functional That Satisfies the Second-Order Density-Gradient Constraint and Has Broad Applicability in Chemistry. J. Chem. Phys. 2011, 135, 191102. [Google Scholar] [CrossRef] [PubMed]
  34. Marenich, A.; Cramer, C.; Truhlar, D. Universal Solvation Model Based on Solute Electron Density and a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378–6396. [Google Scholar] [CrossRef] [PubMed]
  35. Avogadro: An Open-Source Molecular Builder and Visualization Tool—Version 1.2.0. 2016. Available online: http://avogadro.openmolecules.net (accessed on 6 October 2016).
  36. Hanweel, M.; Curtis, D.E.; Lonie, D.; Vandermeersch, T.; Zurek, E.; Hutchison, G. Avogadro: An Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J. Cheminform. 2012, 4, 17. [Google Scholar] [CrossRef] [PubMed]
  37. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and Testing of a General AMBER Force fFeld. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef] [PubMed]
  38. Kronik, L.; Stein, T.; Refaely-Abramson, S.; Baer, R. Excitation Gaps of Finite-Sized Systems from Optimally Tuned Range-Separated Hybrid Functionals. J. Chem. Theory Comput. 2012, 8, 1515–1531. [Google Scholar] [CrossRef] [PubMed]
  39. Lima, I.T.; Prado, A.d.S.; Martins, J.B.L.; de Oliveira Neto, P.H.; Ceschin, A.M.; da Cunha, W.F.; da Silva Filho, D.A. Improving the Description of the Optical Properties of Carotenoids by Tuning the Long-Range Corrected Functionals. J. Phys. Chem. A 2016, 120, 4944–4950. [Google Scholar] [CrossRef] [PubMed]
  40. Colzani, M.; de Maddis, D.; Casali, G.; Carini, M.; Vistoli, G.; Aldini, G. Reactivity, Selectivity, and Reaction Mechanisms of Aminoguanidine, Hydralazine, Pyridoxamine, and Carnosine as Sequestering Agents of Reactive Carbonyl Species: A Comparative Study. ChemMedChem 2016, 11, 1778–1789. [Google Scholar] [CrossRef] [PubMed]
  41. Domingo, L.R.; Perez, P. The Nucleophilicity N Index in Organic Chemistry. Org. Biomol. Chem. 2011, 9, 7168–7175. [Google Scholar] [CrossRef] [PubMed]
  42. Domingo, L.R.; Chamorro, E.; Perez, P. Understanding the Reactivity of Captodative Ethylenes in Polar Cycloaddition Reactions. A Theoretical Study. J. Org. Chem. 2008, 73, 4615–4624. [Google Scholar] [CrossRef] [PubMed]
  43. Pratihar, S.; Roy, S. Nucleophilicity and Site Selectivity of Commonly Used Arenes and Heteroarenes. J. Org. Chem. 2010, 75, 4957–4963. [Google Scholar] [CrossRef] [PubMed]
  • Sample Availability: Not Available.
Figure 1. Molecular Structures of (a) ALT-711; (b) ALT-946; (c) Aminoguanidine; (d) Carnosine; (e) GLY-230; (f) LR-9; (g) Metformin; (h) OPB-9195; (i) Pentoxifylline; (j) Pioglitazone; (k) Pyridoxamine; and (l) Tenilsetam.
Figure 1. Molecular Structures of (a) ALT-711; (b) ALT-946; (c) Aminoguanidine; (d) Carnosine; (e) GLY-230; (f) LR-9; (g) Metformin; (h) OPB-9195; (i) Pentoxifylline; (j) Pioglitazone; (k) Pyridoxamine; and (l) Tenilsetam.
Molecules 22 00226 g001
Table 1. Average descriptors J I , J A , J H L , J χ , J η , J ω , J D 1 , J ω + , J ω , J Δ ω ± , and J D 2 for the ALT-711, ALT-946, aminoguanidine, carnosine, GLY-230, LR-9, metformin, OPB-9195, pentoxifylline, pioglitazone, pyridoxamine, and tenilsetam calculated with the M11, M11L, MN12L, MN12SX, N12, N12SX, SOGGA11, and SOGGA11X density functionals and the Def2TZVP basis set using water as solvent simulated with the solvation model density (SMD) parametrization of the IEF-PCM (Integral Equation Formalism for the Polarized Continuum Model) model.
Table 1. Average descriptors J I , J A , J H L , J χ , J η , J ω , J D 1 , J ω + , J ω , J Δ ω ± , and J D 2 for the ALT-711, ALT-946, aminoguanidine, carnosine, GLY-230, LR-9, metformin, OPB-9195, pentoxifylline, pioglitazone, pyridoxamine, and tenilsetam calculated with the M11, M11L, MN12L, MN12SX, N12, N12SX, SOGGA11, and SOGGA11X density functionals and the Def2TZVP basis set using water as solvent simulated with the solvation model density (SMD) parametrization of the IEF-PCM (Integral Equation Formalism for the Polarized Continuum Model) model.
J I J A J HL J χ J η J ω J D 1 J ω J ω + J Δ ω ± J D 2
M112.4462.3333.3820.0744.7790.6614.8380.9991.0522.0462.513
M11L0.3520.1960.4200.0920.5480.1810.6010.3440.3470.6810.846
MN12L0.3110.1980.3770.0720.4910.1330.5280.2420.2490.4870.604
MN12SX0.0690.1080.1430.0670.1250.0540.1640.1370.0700.2070.260
N120.4590.3550.6010.0880.8110.3060.8960.5680.5671.1351.396
N12SX0.0700.1020.1460.0540.1640.0430.1840.1030.0680.1560.204
SOGGA110.4890.4480.6990.1340.9370.4121.0600.8090.7591.5511.918
SOGGA11X0.9440.9581.3460.0321.9020.3811.9530.6460.6391.2861.575

Share and Cite

MDPI and ACS Style

Frau, J.; Glossman-Mitnik, D. Chemical Reactivity Theory Study of Advanced Glycation Endproduct Inhibitors. Molecules 2017, 22, 226. https://doi.org/10.3390/molecules22020226

AMA Style

Frau J, Glossman-Mitnik D. Chemical Reactivity Theory Study of Advanced Glycation Endproduct Inhibitors. Molecules. 2017; 22(2):226. https://doi.org/10.3390/molecules22020226

Chicago/Turabian Style

Frau, Juan, and Daniel Glossman-Mitnik. 2017. "Chemical Reactivity Theory Study of Advanced Glycation Endproduct Inhibitors" Molecules 22, no. 2: 226. https://doi.org/10.3390/molecules22020226

APA Style

Frau, J., & Glossman-Mitnik, D. (2017). Chemical Reactivity Theory Study of Advanced Glycation Endproduct Inhibitors. Molecules, 22(2), 226. https://doi.org/10.3390/molecules22020226

Article Metrics

Back to TopTop