- Research article
- Open Access
Computational molecular characterization of the flavonoid rutin
Chemistry Central Journalvolume 4, Article number: 12 (2010)
In this work, we make use of a model chemistry within Density Functional Theory (DFT) recently presented, which is called M05-2X, to calculate the molecular structure of the flavonoid Rutin, as well as to predict the infrared (IR) and ultraviolet (UV-Vis) spectra, the dipole moment and polarizability, the free energy of solvation in different solvents as an indication of solubility, the HOMO and LUMO orbitals, and the chemical reactivity parameters that arise from Conceptual DFT. The calculated values are compared with the available experimental data for this molecule as a means of validation of the used model chemistry.
Flavonoids are phenolic substances characterized for a low molecular weight and they are abundant in plant tissues, apple being one of the most important (particularly its skin) [1, 2]. In the human body they show a lot of biological properties as antioxidants, antiallergenic, antibacterial, antifungal, antiviral and anticarcinogenic agents. These characteristics confer to them pharmacological properties useful for the treatment of diseases that go from allergies, bacterial and viral infectious processes, to those of greater risk like the coronary diseases, cancer and HIV [3–5]. The mechanism by which flavonoids carry out their properties, mainly their antioxidant power, is either by inhibiting the formation or activity of reactive oxygen species, or by direct interaction with DNA, enzymes and membrane receptors.
Theoretical investigations of the physical and chemical properties of flavonoids are very important in order to disclose the relationship between the structure, properties and performance, and to help in the design and synthesis of new derivatives with improved properties. We have experimentally found that some natural flavonoids have a strong ability for complexing metal ions, in particular, those related to heavy metals [6–8]. Thus, natural flavonoids could be useful in water treatment, cleaning and purification. The objective of this letter is to report the results of the calculation of the molecular structure and properties of the flavonoid Rutin using a recently developed density functional . The IR and UV-Vis spectra, the dipole moment and polarizability, the free energy of solvation in different solvents as an indication of solubility, the HOMO and LUMO orbitals, and the chemical reactivity parameters that arise from Conceptual DFT [10, 11] are reported. The calculated values are compared with the available experimental data for this molecule as a means of validation of the used model chemistry. The spectra and the calculated values are important in the sense that they are an indication of the chemical stability, the thermochemistry, the color and the region of the solar spectrum where the absorption takes place, the solubility and the chemical reactivity which is useful to predict the possible complexation sites.
Theory and Computational Details
For all the calculations, we have chosen the hybrid meta-GGA M05-2X functional , which consistently provides satisfactory results for several structural and thermodynamic properties. Although there are a new class of functionals, the so called M06 functionals, our own experience indicates that the improvement in the calculated molecular structure and properties of systems of the size that we are considering in this paper is only marginal. The 3-21G(d) basis set was used for the geometry optimizations and evaluations of harmonic frequencies both in the gas phase and in aqueous solution of the flavonoid. It has been found that this basis set has a remarkable ability to predict the molecular structure and properties of large systems when coupled withe B3LYP density functional , and the same has been our own experience when the M05-2X functional is considered. The equilibrium geometry of the studied molecule was determined by means of the gradient technique. The force constants and vibrational frequencies were determined from calculation using the FREQ keyword on the stationary points obtained after the optimization to check if they were true minima. The electronic properties were calculated with the 6-31+G(d) basis set. A suitable description of this basis set is provided in some of the most important Computational Chemistry recent books [13–16]. Solvation energies were computed by the Integral Equation Formalism - Polarizable Continuum Model (IEF-PCM) , including the UAHF model. All the calculations have been performed with the Gaussian 03W series of programs .
The calculation of the ultraviolet (UV-Vis) spectra of the flavonoid and their metallic complexes has been performed by solving the time dependent Kohn-Sham equations according to the method implemented in Gaussian 03W [13, 19–21]. The equations have been solved for 10 excited states.
The infrared (IR) and ultraviolet (UV-Vis) spectra were calculated and visualized using the Swizard program . In all cases the displayed spectra show the calculated frequencies and absorption wavelengths. The vertical ionization potential I and electron affinity A were calculated in two ways: i) as the difference between the total energy of the neutral molecule and the corresponding ions, taken at the geometry of the neutral molecule in order to keep the external potential constant, and ii) considering the approximation given by the Koopmans' theorem [13–16], where the HOMO energy is equal to -I and the LUMO energy is equal to -A.
Results and Discussion
The results for the equilibrium conformation of the neutral molecule of 2-(3,4-dihydroxyphenyl)-5,7-dihydroxy-3-[(2S,3R,4 S,5S,6R)-3,4,5-trihydroxy-6-([(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl) oxan-2 yl]oxy-4H-chromen-4-one (or rutin, for short) calculated with the M05-2X/3-21G(d) model chemistry are presented in Figure 1 through a representation of the molecular structure showing the atomic labeling and numbering as well as the interatomic bond lengths and several selected angles. The results have been compared with those of the quercetin molecule [25, 26]. As the comparison of the computed molecular structure of quercetin with the X-ray results has been presented already [25, 26], we are not repeating it here. However, it can be said that the agreement between the computed quercetin moiety of rutin and the quercetin X-ray results is very good. It should be remarked that there are two H-bonds belonging to O30 with H29, and O68 with H31. This could be an explanation of the increased water, methanol and ethanol solubility (see below).
The infrared spectrum (IR) for the rutin molecule calculated with the M05-2X/3-21G(d) model chemistry is displayed in Figure 2. The vibrational band assignments have been done using the ChemCraft for Windows molecular visualization program. By comparison with the experimental IR spectrum , an average scaling factor of 0.995 could be estimated. In Table 1, we present a comparison of the experimental, computed and scaled frequencies for the rutin molecule as an assessment of the M05-2X functional for calculating vibrational frequencies. However, it must be noted that a strong peak at 2747 cm-1 related to an internal H-bond has been omitted in order not to obscure the rest of the spectrum and for the sake of clarity. Thus, it is expected that the model chemistry used in this work can reproduce the experimental spectrum of the rutin molecule with a certain degree of accuracy.
The ultraviolet spectrum (UV-Vis) of the rutin molecule was calculated with TDDFT using the M05-2X/6-31+G(d,p) model chemistry. The results in Table 2 show the first ten electronic transitions states of rutin, both in nanometers (nm) and electron-volts (eV), the oscillators strengths (f) that can give an idea of the intensity of the transition, and the orbital assignments, indicating the percentage of any H-n → L+n transition involved. The wavelength belonging to the HOMO-LUMO transition will take place at 290 nm. As the HOMO-LUMO transition takes place in the ultraviolet region, close to but out of the visible zone, it can be predicted that this molecule will be colorless or slightly colored.
The molecular dipole moment is perhaps the simplest experimental measure of charge distribution in a molecule. The accuracy of the overall distribution of electrons in a molecule is hard to quantify, since it involves all the multipoles. The polarizability a contributes to the understanding of the response of the system when the external field is changed, while the number of electrons N is kept fixed. The polarizability is calculated as the average of the polarizability tensor . From the present calculations, the total energy, the total dipole moment and the isotropic polarizability of the ground state with the 6-31+G(d,p) model chemistry are -2250.341 au, 7.6877 Debye and 177.57 Bohr3 for the rutin molecule. These results for the dipole moment and the isotropic polarizability could be of interest as an indication of the solubility and chemical reactivity of the studied molecule, not only for it synthesis but for the potential application in complexation of metal cations for water cleaning and purification.
The free energy of solvation ΔG(solv) of the molecule have been calculated for rutin by resorting to the M05-2X/6-31+G(d,p) model chemistry coupled with the Integral-Equation-Formalism of the Polarized Continuum Model (IEF-PCM) for different solvents as implemented in Gaussian 03. The solubility of a molecule will depend on several kinetic and thermodynamic factors. However, it can be said that the magnitude and the sign of ΔG(solv) could be a good approximation as an index of solubility. In this way, a negative sign and a large magnitude will be an indication of increased solubility. The results of these calculations for the studied molecule can be summarized as follows: Acetone = -18.36 Kcal/mol, Acetonitrile = -6.05 Kcal/mol, Aniline = 11.36 Kcal/mol, Benzene = 0.10 Kcal/mol, CCl4 = -0.14 Kcal/mol, Chlorobenzene = -5.42 Kcal/mol, Chloroform = -9.03 Kcal/mol, Cyclohexane = -6.32 kcal/mol, Dichloroethane = -13.37 Kcal/mol, Dichloromethane = -15.59 Kcal/mol, Diethylether= -13.74 Kcal/mol, DMSO = -15.05 Kcal/mol, Ethanol = -52.88 Kcal/mol, Heptane = -7.73 Kcal/mol, Methanol = -56.35 Kcal/mol, Nitromethane = -14.48 Kcal/mol, THF = - 10.58 Kcal/mol, Toluene = -2.36 Kcal/mol, and Water = -49.42 Kcal/mol. These values could be an indication that the studied molecule will be mostly soluble in ethanol, methanol, and water, and this can be related to the results obtained for the dipole moment and polarizability.
The HOMO and LUMO of rutin calculated with the M05-2X/6-31+G(d,p) model chemistry are displayed in Figure 3. The HOMO and LUMO densities are over the flavonoid moiety, but not over the glycoside rest. This can give us an idea of the reactivity of the molecule and means that only the flavonoid moiety will be reactive, for example, in complexation with metal cations.
Within the conceptual framework of DFT , the chemical potential μ, which measures the escaping tendency of an electron from equilibrium is defined as:
where χ is the electronegativity.
The global hardness η can be seen as the resistance to charge transfer:
where ϵ H and ϵ L are the energies of the highest occupied and the lowest unoccupied molecular orbitals, HOMO and LUMO, respectively.
The electrophilicity index ω represents the stabilization energy of the system when it gets saturated by electrons coming from the surrounding:
The validity of the Koopmans' theorem within the DFT approximation is controversial. However, it has been shown  that although the KS orbitals may differ in shape and energy from the HF orbitals, the combination of them produces Conceptual DFT reactivity descriptors that correlate quite well with the reactivity descriptors obtained through Hartree-Fock calculations. Thus, it is worth to calculate the electronegativity, global hardness and global electrophilicity for the rutin molecule using both approximations in order to verify the quality of the procedures.
The results for the vertical I and A of the rutin molecule obtained through energy differences between the ionized and the neutral state, calculated at the geometry of the neutral molecule are I = 7.284 eV and A = 0.067 eV. The HOMO and LUMO energies are -6.836 eV and -0.522 eV, respectively. It can be seen that there is a good qualitative agreement between both results for I, but not for A. The calculated values of the electronegativity, global hardness and global electrophilicity using the I and A are χ = 3.676 eV, = 3.608 eV, and ω = 1.873 eV. Using the HOMO and LUMO energies, within the Koopmans' theorem, the corresponding values are χ = 3.679 eV, η = 3.158 eV, and ω = 2.143 eV. Again, there is a good qualitative agreement for the reactivity parameters calculated through both procedures. It can be concluded that for the particular case of the rutin molecule, the M05-2X/6-31+G(d,p) model chemistry is able to predict the Conceptual DFT reactivity indices calculated through HOMO and LUMO energies as well as from the I and A obtained through energy differences with qualitative similar good accuracy.
The condensed Fukui functions can also be employed to determine the reactivity of each atom in the molecule. The corresponding condensed functions are given by (for nucleophilic attack), (for electrophilic attack), and (for radical attack), where qk is the gross charge of atom k in the molecule.
It is possible to evaluate condensed Fukui functions from single-points calculations directly, without resorting to additional calculations involving the systems with N-1 and N+1 electrons:
with c ai being the LCAO coefficients and S ab the overlap matrix. The condensed Fukui functions are normalized, thus and .
The results from the calculation of the condensed Fukui functions for nucleophilic, electrophilic and radical attack have been obtained by resorting to the AOMix molecular analysis program . The sites for electrophilic attack will be those atoms bearing a negative charge and where the Fukui function is a maximum. These values confirm that the sites for the electrophilic attack are the C12 and C22 atoms. The site for potential nucleophilic attack would depend on the values of on the atoms with a positive charge density. The results indicate that the site for nucleophilic attack will be the C11 and C13 atoms. Finally, the site for radical attack, governed by the values of will be the C12 atom. The glycoside rings are not very reactive in this context.
In this work, the M05-2X/3-21G(d) and M05-2X/6-31+G(d,p) model chemistries have been applied to the study of a molecule which is potentially useful for water cleaning and purification. The molecular structure for the rutin molecule has been determined by using the M05-2X/3-21G(d) model chemistry. A comparison has been made with the results from the experimental X-ray crystallography for the flavonoid quercetin. The agreement is generally good. Two internal H-bonds have been described that could be an explanation for the increased solubility of rutin in water, ethanol and methanol.
The shape of the frontier orbitals of this molecule were displayed as well as some electronic parameters like the total energy, the dipole moment and the polarizability and the infrared (IR) and ultraviolet (UV-Vis) spectra for the rutin molecule have been predicted according to the M05-2X/6-31+G(d,p) model chemistry, and an assignment of the principal peaks has been achieved.
The free energy of solvation ΔG(solv) of the rutin molecule has been calculated by resorting to the M05-2X/6-31+G(d,p) model chemistry coupled with the Integral-Equation-Formalism of the Polarized Continuum Model (IEF-PCM) for different solvents and the results gave an indication of water, ethanol and methanol as the solvents in which this molecule could be potentially soluble.
The ionization potential I and the electron affinity A have been calculated through energy differences between the ionic and the neutral states, all at the geometry of the neutral molecule, and they have been compared well with the results obtained from the HOMO and LUMO energies obtained through the Koopmans' theorem procedure. The results indicate a qualitative good agreement, which can be consider an indication of the goodness of the proposed model chemistry.
The M05-2X density functional in combination with several basis sets appears to be a useful tool for the study of the molecular structure and electronic properties of flavonoids and the possible nanostructures derived from them, and further applications to several molecular systems of this kind are being pursued in our laboratory.
Pietta PG: Flavonoids as antioxidants. J Nat Prod. 2000, 63: 1025-1042. 10.1021/np9904509.
Wolfe K, Wu X, Liu RH: Antioxidant activity of apple peels. J Agric Food Chem. 2003, 51: 609-614. 10.1021/jf020782a.
Scalbert A, Williamson G: Dietary intake and bioavailability of polyphenols. J Nutr. 2000, 130: 2073S-2085S.
Erkoc S, Erkoc F, Keskin N: Theoretical investigation of quercetin and its radical isomers. J Mol Struct - THEOCHEM. 2003, 631: 141-146. 10.1016/S0166-1280(03)00237-9.
Russo N, Toscano M, Uccella N: Semiempirical molecular modeling into quercetin reactive site: Structural, conformational, and electronic features. J Agric Food Chem. 2000, 48: 3232-3237. 10.1021/jf990469h.
Payán-Gómez SA, Flores-Holguín N, Pérez-Hernández A, Pinón-Miramontes M, Glossman-Mitnik D: Computational molecular characterization of the flavonoid morin and its Pt(II), Pd(II) and Zn(II) complexes. J Mol Mod
Lekka CE, Ren J, Meng S, Kaxiras E: Structural, electronic and optical properties of representative Cu-flavonoid complexes. J Phys Chem B. 2009, 113: 6478-6483. 10.1021/jp807948z.
Malesev D, Kuntic V: Investigation of metal-flavonoid chelates and the determination of flavonoids via metal-flavonoid complexing reactions. J Serb Chem Soc. 2007, 72: 921-939. 10.2298/JSC0710921M.
Zhao Y, Truhlar DG: Density functionals with broad applicability in chemistry. Acc Chem Res. 2008, 41: 157-167. 10.1021/ar700111a.
Parr RG, Yang W: Density Functional Theory of Atoms and Molecules. 1989, New York: Oxford University Press
Geerlings P, DeProft F, Langenaeker W: Conceptual density functional theory. Chem Rev. 2003, 103: 1793-1873. 10.1021/cr990029p.
Zandler ME, D'Souza F: The remarkable ability of B3LYP/3-21G(*) calculations to describe geometry, spectral and electrochemical properties of molecular and supramolecular porphyrin-fullerene conjugates. Comptes Rendus Chimie. 2006, 9: 960-981. 10.1016/j.crci.2005.12.008.
Lewars E: Computational Chemistry - Introduction to the Theory and Applications of Molecular and Quantum Mechanics. 2003, Norwell, MA: Kluwer Academic Publishers
Young D: Computational Chemistry: A Practical Guide for Applying Techniques to Real-World Problems. 2001, New York: John Wiley & Sons
Jensen F: Introduction to Computational Chemistry. 2007, Chichester: John Wiley & Sons
Cramer CJ: Essentials of Computational Chemistry - Theories and Models. 2002, Chichester: John Wiley & Sons
Tomasi J, Mennucci B, Cancès E: The IEF version of the PCM solvation method: An overview of a new method addressed to study molecular solutes at the QM ab initio level. J Mol Struct - THEOCHEM. 1999, 464: 211-226. 10.1016/S0166-1280(98)00553-3.
Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Montgomery JA, Vreven T, Kudin KN, Burant JC, Millam JM, Iyengar SS, Tomasi J, Barone V, Mennucci B, Cossi M, Scalmani G, Rega N, Petersson GA, Nakatsuji H, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Klene M, Li X, Knox JE, Hratchian HP, Cross JB, Bakken V, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Ayala PY, Morokuma K, Voth GA, Salvador P, Dannenberg JJ, Zakrzewski VG, Dapprich S, Daniels AD, Strain MC, Farkas O, Malick DK, Rabuck AD, Raghavachari K, Foresman JB, Ortiz JV, Cui Q, Baboul AG, Clifford S, Cioslowski J, Stefanov BB, Liu G, Liashenko A, Piskorz P, Komaromi I, Martin RL, Fox DJ, Keith T, Al-Laham MA, Peng CY, Nanayakkara A, Challacombe M, Gill PMW, Johnson B, Chen W, Wong MW, Gonzalez C, Pople JA: Gaussian 03, Revision E.01. 2004, Gaussian, Inc., Wallingford, CT
Stratmann RE, Scuseria GE, Frisch MJ: An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules. J Chem Phys. 1998, 109: 8218-8224. 10.1063/1.477483.
Bauernschmitt R, Ahlrichs R: Treatment of electronic excitations within the adiabatic approximation of time-dependent density functional theory. Chem Phys Lett. 1996, 256: 454-464. 10.1016/0009-2614(96)00440-X.
Casida ME, Jamorski C, Casida KC, Salahub DR: Molecular excitation energies to high-lying bound states from time-dependent density-functional response theory: Characterization and correction of the time-dependent local density approximation ionization threshold. J Chem Phys. 1998, 108: 4439-4449. 10.1063/1.475855.
Gorelsky SI: Swizard program. [http://www.sg-chem.net/]
ChemCraft 1.6. [http://www.chemcraftprog.com]
Gorelsky SI: AOMix program. [http://www.sg-chem.net/]
Mendoza-Wilson AM, Glossman-Mitnik D: CHIH-DFT deternmination of the molecular structure, infrared and ultraviolet spectra of the flavonoid quercetin. J Mol Struct - THEOCHEM. 2004, 681: 71-76. 10.1016/j.theochem.2004.04.054.
Mendoza-Wilson AM, Glossman-Mitnik D: CHIH-DFT study of the electronic properties and chemical reactiviy of quercetin. J Mol Struct - THEOCHEM. 2005, 716: 67-72. 10.1016/j.theochem.2004.10.083.
Fossen T, Andersen OM: Spectroscopic techniques applied to flavonoids. Flavonoids - Chemistry, Biochemistry and Applications. Edited by: Andersen OM, Markham KR, Boca Ratón. 2006, FL: Taylor and Francis, LLC, 1: 37-142. 1
Zevallos J, Toro-Labbé A: A theoretical analysis of the Kohn-Sham and Hartree-Fock orbitals and their use in the determination of electronic properties. J Chil Chem Soc. 2003, 48: 39-47. 10.4067/S0717-97072003000400007.
This work has been partially supported by Consejo Nacional de Ciencia y Tecnología (CONACYT, Mexico) and by Fondo Mixto del Estado de Baja California (FOMIX-BC) through Project 69363. SAPG gratefully acknowledges a fellowship from CONACYT. NFH, APH and DGM are researchers of CONACYT and CIMAV.
The authors declare that they have no competing interests.
SAPG carried out the DFT calculations and wrote some portions of the paper. NFH analysed the IR and UV-Vis spectra. APH and DGM conceived of the study, analysed all the data, discussed the results and wrote the manuscript. MPM discussed the results and revised the manuscript. All authors read and approved the final manuscript.