Abstract

Five density functionals, CAM-B3LYP, LC-ωHPBE, MN12SX, N12SX, and ωB97XD, in connection with the Def2TZVP basis set were assessed together with the SMD solvation model for the calculation of the molecular properties, chemical reactivities, and solubilities of some pigments derived from astaxanthin, β-cryptoxanthin, fucoxanthin, myxol, siphonaxanthin, siphonein, and zeaxanthin marine carotenoids in the presence of different solvents (hexane, methanol, ethanol, and water). All the chemical reactivity descriptors for the systems were calculated via conceptual density functional theory (CDFT). Finally, the potential bioavailability and druggability as well as the bioactivity scores for the marine carotenoid pigments were predicted through different methodologies already reported in the literature, which have been previously validated during the study of other natural products obtained from marine sources.

1. Introduction

Carotenoids are pigments widely found in nature, either in plants, animals, or microorganisms, forming a group consisting of over 700 compounds. Carotenoids stand as the major group of compounds used as color additives. These natural pigments are responsible for many of the colors seen in edible fruits, vegetables, mushrooms, flowers, and even lobster and trout hues from the animal kingdom [1].

Due to their inherent coloration, these pigments have been proven to be useful as food colorants [2] and as sensitizers for applications in organic photovoltaic solar cells and artificial photosynthesis [3]. Our group has earlier performed computational chemistry and molecular modeling studies on some carotenoid pigments that could be potentially considered to be the energy absorbing material of dye-sensitized solar cells [46]. Additionally, these compounds function as sources of provitamin A, are capable of absorbing sunlight, and can act as oxygen transporters, powerful quenchers of singlet oxygen, and free radical scavengers.

Marine-derived bioactive compounds offer an abundant source of pharmacologically active agents with great chemical diversity and complexity and the potential to produce valuable therapeutic entities [7]. Due to this, a lot of research has been performed during the last two decades dedicated to the search for new natural products that can be acquired from the knowledge of species of marine origin through the field of marine pharmacognosy [8].

Among the chemical species that can be obtained from natural products of marine origin, carotenoid pigments are molecules of relatively small size. The therapeutic application of these pigments is currently an active field of research because they could have a strong potential for the treatment of several ailments [9, 10].

In order to consider these carotenoid pigments as the starting point for the development of medical drugs, it is mandatory to acquire knowledge about their chemical reactivity properties as well as the bioactivity associated with them. From the basics of medicinal chemistry, it is known that drugs exert their effect by interacting with the active site of a receptor, which is generally a protein [11]. These interactions rely on the different kinds of bindings between the pharmacophore and the chemical groups present in the active site and thus are intimately related to their chemical reactivity from a molecular perspective [12, 13]. One of the most powerful tools to understand the chemical reactivity of interacting molecular systems within computational chemistry is probably conceptual DFT (density functional theory) [14, 15], also called chemical reactivity theory, which allows us to accomplish this task by resorting to several global and local descriptors which are in turn related to variations in the electronic densities of the studied systems.

On the basis of previous considerations, the objective of this work is to study the chemical reactivity of a group of carotenoids of marine origin, astaxanthin, β-cryptoxanthin, fucoxanthin, myxol, siphonaxanthin, siphonein, and zeaxanthin, which are representative of the most abundant xanthophylls present in ocean plants and animals, using the techniques of conceptual DFT, determining their global reactivity properties, that is, of the molecule as a whole. Moreover, during the process of development of new drugs, there is a need to learn about the drug-like properties of the involved molecular systems [16]. Thus, the descriptors of bioavailability and bioactivity (bioactivity scores) will be calculated through different procedures described in the literature [17, 18], trying to relate them with the calculated conceptual DFT descriptors.

2. Computational Methodology

In the same way as we have proceeded in our recent studies [1928], the computational tasks in this work have been performed by considering the popular Gaussian 09 software [29]. Following the conclusions obtained from those studies, five density functionals have been chosen, CAM-B3LYP, LC-ωHPBE, MN12SX [30], N12SX, and ωB97XD, because they can be considered to be well suited for our purposes according to our proposed KID (for Koopmans in DFT) criteria [1928]. For the calculation of the electronic properties, several model chemistries have been considered, based on the mentioned density functionals in connection with the Def2TZVP basis set, while a smaller Def2SVP was considered for the prediction of the most stable structures [31, 32]. In order to obtain accurate results, all calculations were performed using four different solvents, namely hexane, methanol, ethanol, and water, simulated with the SMD model [33]. While water is the universal biological solvent, the others were chosen because a good solvent selection is necessary for easy and effective synthesis of compounds used in drugs.

3. Results and Discussion

The molecular structures of the optimized conformers of the seven carotenoids of marine origin as presented in Figure 1 were optimized in the gas phase by means of the DFTBA model available in the software and then reoptimized with the five density functionals described previously, the Def2SVP basis set, and hexane, methanol, ethanol, or water as the solvent. The calculation of the electronic properties was performed by using the same model chemistries but changing the basis set with the Def2TZVP one.

In order to verify the fulfillment of our proposed KID procedure, it is necessary to perform a comparison of the orbital energies with the results obtained by means of vertical I and A through the SCF criterion. To this end, the three main descriptors are linked as with and with , and their behavior in describing the HOMO-LUMO gap is defined as , , and . Another descriptor, , the difference between the SOMO and LUMO, was also designed to guide the verification of the accuracy of approximation [1925, 27, 28]. The results of this analysis are presented in Tables 14.

The overall conclusion that can be extracted from the inspection of the results presented in Tables 14 is that in agreement with our previous studies on melanoidins and peptides, the model chemistries involving the MN12SX and N12SX density functionals are the best for verifying our proposed criteria of good behavior, that is, the values of , , , and are close to zero. Moreover, it can be seen that the KID procedure is better verified for the case of the polar solvents rather than for hexane.

It is possible to obtain additional information about the solubility of these carotenoid pigments by means of the calculation of the free energy of solvation in every solvent according to each of the model chemistries. The results are presented in Table 5, where it can be seen that while all the carotenoid pigments are soluble in all solvents, their relative solubility is greater in methanol, ethanol, and hexane than in water. As it can be observed from this table, β-cryptoxanthin has the lowest solvation energy in water. This particular molecule has only one hydroxyl and its polarity is lower than that from the other analyzed systems.

4. Calculation of Global Reactivity Descriptors

By taking into account the KID procedure presented in our previous works together with the finite difference approximation, the global reactivity descriptors can be expressed as follows:Electronegativity: [14, 15]Global Hardness: [14, 15]Electrophilicity: ω =  [34]Electrodonating power:  =  [35]Electroaccepting power:  =  [35]Net electrophilicity: [36]

Where and are the energies of the HOMO and LUMO, respectively.

The results for the global reactivity descriptors for the carotenoids pigments considered in this study based on the values of the HOMO and LUMO energies calculated with the five density functionals are presented in Tables 612:

As expected from the molecular structure of these carotenoids, their electrodonating ability is more important that their electroaccepting character. However, on comparing the results for global hardness and electrophilicity of all the systems in the four solvents, it can be seen that MN12SX and N12SX density functionals (which verify the KID criteria) give results different from those obtained from calculation with the other three density functionals. It is also well known that the MN12SX density functional has the best across-the-board performance of all the considered functionals for chemistry and physics and for energies and structures.

The electrophilicity (ω) index encompasses the balance between the tendency of an electrophile to acquire an extra magnitude of electron density and the resistance of a molecule to exchange electron density with the environment [37]. By studying the electrophilicity of a series of reagents involved in Diels–Alder reactions [38], through quantitative characterization of the global electrophilicity pattern of some reagents involved in 1,3-dipolar cycloaddition reactions [39], and by means of understanding of the mechanism of polar Diels–Alder reactions [40], Domingo et al. established an electrophilicity (ω) scale for the classification of organic molecules as strong electrophiles with ω > 1.5 eV, moderate electrophiles with 0.8 < ω < 1.5 eV, and marginal electrophiles with ω < 0.8 eV [38].

By inspection of Tables 712, it can be seen that all the carotenoids considered in this study can be regarded as strong electrophiles, irrespective of the solvents in which the calculations have been performed. The values for the global electrophilicity obtained by using the N12SX density functional are a little smaller than those achieved with the MN12SX density functional, but the conclusions are the same. Besides the electrophilicity classification, an electrophilicity scale for the marine carotenoids can be displayed as follows: myxol > astaxanthin > siphonaxanthin > siphonein > fucoxanthin > zeaxanthin > β-cryptoxanthin, which is the same for the case of both density functionals and the four solvents considered here.

The nucleophilicity N is another important chemical reactivity descriptor. 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]. However, these indices fail for more complex molecules which display concurrently both electrophilic and nucleophilic behaviors, and in these cases, the value of electrophilicity (ω) does not correlate well with their expected nucleophilicity [37, 42]. This has compelled Domingo and his collaborators [37, 4043] to propose a new nucleophilicity index, N, simply based on the highest occupied molecular orbital (HOMO) energy obtained within the Kohn–Sham scheme with an arbitrary shifting of the origin with tetracyanoethylene (TCE) taken as a reference. The corresponding definition for the nucleophilicity (N) index is [37, 4043], and the results for the calculation of this index for the marine carotenoids using the MN12SX and N12SX density functionals are presented in Table 13.

On the basis of the previous definition and the scale established in the mentioned study [43], it can be concluded that all the marine carotenoids considered within this research can be regarded as strong nucleophiles because their values are greater than 3 eV, independent of the methodology employed for the calculations and the solvent considered in each case. However, the nucleophilic character for all the marine carotenoids is slightly greater in hexane than in the polar solvents (methanol, ethanol, and water).

5. Bioactivity Scores

The molecular properties that are related to the concept of drug-likeness and in particular, those associated with the criteria proposed by Lipinski et al. [44, 45] for the prediction of oral bioavailability have been calculated by feeding the \ponding SMILES notations into the MolInspiration software readily available online (Slovensky Grob, Slovak Republic (http://www.molinspiration.com)). The results are presented in Table 14.

However, what the Lipinski's rule of five really measures is the oral bioavailability of a potential drug because this is the desired property for a molecule having drug-like character. Then, a different approach was followed by considering similarity searches in the chemical space of compounds with structures that can be compared to those that are being studied and with known pharmacological properties. The same software was used for the calculation of the bioactivity scores which are a measure of the ability of the potential drug to interact with different receptors, that is, to act as GPCR ligands or kinase inhibitors, to perform as ion channel modulators, or to interact with enzymes and nuclear receptors. These predicted bioactivity scores were compared with those resulting from the use of alternative software like MolSoft from Molsoft LLC. (San Diego, CA, USA) (http://molsoft.com/mprop/) as well as ChemDoodle Version 9.02 from iChemLabs LLC, Richmond, VA, USA (http://www.chemdoodle.com). The values of the bioactivity scores for the seven carotenoid pigments are presented in Table 15.

These bioactivity scores for organic molecules can be interpreted as active (when the bioactivity score 0), moderately active (when the bioactivity score lies between −5.0 and 0.0), and inactive (when the bioactivity score −5.0). Thus, astaxanthin, β-cryptoxanthin, myxol, siphonaxanthin, and zeaxanthin will be prone to act as nuclear receptor ligands, while fucoxanthin and siphonein will be moderately bioactive, regarding all the interactions.

6. Conclusions

In this paper, we have presented a new study performed on the chemical reactivity of seven carotenoid pigments of marine origin based on conceptual DFT as a tool to explain molecular interactions.

The knowledge of the values of the global and local descriptors of the molecular reactivity of the carotenoid pigments studied could be useful in the development of new drugs based on these compounds or some analogs. On the basis of the analysis of the electrophilicity index ω and Domingo’s nucleophilicity index N, it could be concluded that the studied carotenoids display concurrently both strong electrophilic and nucleophilic behaviors.

Finally, the molecular properties related to bioavailability have been predicted using different methodologies already described in the literature, and the descriptors used for the quantification of the bioactivity allowed us to characterize some of the studied molecules as being bioactive by acting as ligands for a nuclear receptor.

Data Availability

All the data related to this submission are available upon request.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.

Authors’ Contributions

Daniel Glossman-Mitnik conceived and designed the research and headed, wrote, and revised the manuscript, while Juan Frau and Norma Flores-Holguín contributed to the writing and revision of the article.

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. Daniel Glossman-Mitnik conducted this work while being a visiting lecturer at the University of the Balearic Islands, from which support is gratefully acknowledged. Norma Flores-Holguín and Daniel Glossman-Mitnik are researchers of CIMAV and CONACYT. This work was cofunded by the Ministerio de Economía y Competitividad (MINECO) and the European Fund for Regional Development (FEDER) (CTQ2014-55835-R).