Next Article in Journal
Deacetylated Konjac Glucomannan with a Slower Hydration Rate Delays Rice Digestion and Weakens Appetite Response
Previous Article in Journal
Effect of Metal Elements on Microstructure and Mechanical Properties of Ultrafine Cemented Carbide Prepared by SPS
Previous Article in Special Issue
Theoretical Insights on ORR Activity of Sn-N-C Single-Atom Catalysts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

An Update in Computational Methods for Environmental Monitoring: Theoretical Evaluation of the Molecular and Electronic Structures of Natural Pigment–Metal Complexes

by
Gabriella Josephine Maranata
1,
Sandra Megantara
1,2 and
Aliya Nur Hasanah
1,2,*
1
Department of Pharmaceutical Analysis and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Jl. Raya Bandung Sumedang KM 21, 5, Jatinangor, Sumedang 45363, Indonesia
2
Drug Development Study Centre, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang 45363, Indonesia
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(7), 1680; https://doi.org/10.3390/molecules29071680
Submission received: 2 March 2024 / Revised: 3 April 2024 / Accepted: 5 April 2024 / Published: 8 April 2024
(This article belongs to the Special Issue Computational and Theoretical Chemistry for Material Research)

Abstract

:
Metals are beneficial to life, but the presence of these elements in excessive amounts can harm both organisms and the environment; therefore, detecting the presence of metals is essential. Currently, metal detection methods employ powerful instrumental techniques that require a lot of time and money. Hence, the development of efficient and effective metal indicators is essential. Several synthetic metal detectors have been made, but due to their risk of harm, the use of natural pigments is considered a potential alternative. Experiments are needed for their development, but they are expensive and time-consuming. This review explores various computational methods and approaches that can be used to investigate metal–pigment interactions because choosing the right methods and approaches will affect the reliability of the results. The results show that quantum mechanical methods (ab initio, density functional theory, and semiempirical approaches) and molecular dynamics simulations have been used. Among the available methods, the density functional theory approach with the B3LYP functional and the LANL2DZ ECP and basis set is the most promising combination due to its good accuracy and cost-effectiveness. Various experimental studies were also in good agreement with the results of computational methods. However, deeper analysis still needs to be carried out to find the best combination of functions and basis sets.

Graphical Abstract

1. Introduction

Metals are relatively abundant and can be found in almost all aspects of life. They provide many benefits, including serving as enzyme cofactors, participation in regulatory systems, and the ability to support protein structures [1,2]. However, some metals can be very dangerous and cause serious problems, especially when present at high concentrations for a long period of time [3,4,5,6,7,8,9]. Based on the danger they pose, metals can be classified into two distinct groups, namely heavy metals and non-heavy metals.
Heavy metals have high densities and atomic weights and are toxic even in very small amounts [4]. These elements are naturally present within Earth’s crust [10,11] and released into the environment through various naturally occurring processes, anthropogenic activities, and bio-geo processes [12,13]. The term heavy metal usually refers to metals and metalloids that are related to toxicity and pollution [14,15]. Almost all heavy metals are environmental pollutants and have a markedly negative impact on biota. Today, the toxicity of heavy metals is becoming an increasingly rampant problem from ecological, evolutionary, health, and, especially, environmental perspectives [15,16,17,18,19]. Heavy metal contamination can disturb the balance of both soil and aquatic ecosystems, which can lead to a reduction in biota populations and even extinction [20,21]. Apart from the environment, heavy metals are dangerous for the body. For example, mercury (Hg) can disrupt the membrane potential and also disturb intracellular calcium (Ca) homeostasis [22], while cadmium (Cd) causes iron (Fe) deficiency by binding to aspartate, cysteine, glutamate, and histidine ligands [23]. Another example is arsenic (As), which affects the sulfhydryl groups in cells and thus disrupts cell respiration, mitosis, and enzyme function [9].
Non-heavy metals are generally not dangerous in small amounts, but if the amounts are excessive, they can also have negative impacts on the body. Fe, cobalt (Co), manganese (Mn), zinc (Zn), and copper (Cu) can actually be useful and are needed by our body, but toxic effects may occur when they are present in excess [11,24]. For example, Fe is a cofactor for various important proteins and enzymes but can also cause liver fibrosis and cirrhosis [25]; increase the potential for heart failure and arrhythmias [26,27]; and increase the possibility of osteoporosis, hypothyroidism, and hypoparathyroidism [28]. In addition, Cu and chromium (Cr), among other non-heavy metals, have a very small beneficial concentration range [4].
To protect the health of people, regulatory organisations have proposed maximum amounts of certain metals (Table 1) [7,24,29,30]. In reality, the dangers of metals are very difficult to control because often the amount of exposure is not clearly known. At present, environmental pollution by metals, especially heavy metals, is an unavoidable threat and challenge, especially with increasing industrial activity without proper waste processing [31,32,33,34,35] and an increase in traffic volume [36]. Pollution spreads through the air, water, animal products, agricultural products, and various other components into the surroundings [37,38,39,40].
Metals can enter plants, animals, and humans in various ways, such as through the respiratory tract, food and drink, and even skin contact. Due to their abundance and the fact that they cannot be seen with the naked eye, metals can enter tissues unnoticed over long periods of time and accumulate in various organs, including the brain, kidneys, liver, and heart. After accumulation, they can disrupt normal biological systems and functions [11,14,16,17,41].
Detecting the presence of metals is crucial to ensuring that organisms and the environment are protected from their negative effects. Many researchers detect metal levels by using atomic absorption spectroscopy [42,43], flame spectroscopy [44], X-ray fluorescence spectrometry [45], inductively coupled plasma mass spectroscopy [36], neutron activation analysis, and inductively coupled plasma-optical emission spectrometry [46,47], among others. These methods show good selectivity and sensitivity; however, they lack efficiency when viewed in terms of time, energy, and costs. Moreover, the availability of instruments is still limited [20,39]. Due to existing limitations, the development of portable detectors that are efficient in terms of time, cost, and energy could represent a solution. While many portable detectors have been developed, some of them are prepared using reagents and materials that may cause danger and produce contamination. Therefore, using natural pigments in the development of portable metal detectors is a potential solution to create detectors that are safe for the environment, simple, do not require chemical reagents, and are efficient in terms of energy and time [48].
Natural pigments—coloured compounds from various sources in nature—can be a solution to this need [49]. These pigments provide colour to fruits and vegetables—for example, betalains produce a reddish-yellow colour, anthocyanins give a reddish-blue colour, carotenoids give a yellow-to-orange colour, and chlorophylls give a green colour [50]. Pigments have chromophores in their structure, so they absorb light in the ultraviolet and visible regions [51]. The chromophore group captures energy by attracting electrons from lower to higher energy orbitals, and then the remaining energy that is not absorbed is reflected and appears as colour to the human eye [50]. Examples of natural pigments that are often found and used in the development of metal detectors from natural pigments can be seen in Figure 1.
Natural pigments can form complexes with metals, and these complexes result in colour changes. Hence, portable metal detectors could incorporate natural pigments to detect metals based on colour changes [52,53,54]. Several studies have shown that pigments and metals are capable of forming complexes. Wybraniec et al. (2013) [55] reported that betalains can interact with nickel (Ni) and Cu, decreasing the absorbance of the pigment. Saithongdee et al. (2014) [56] reported that curcumin forms a complex with Fe and provides a shift from yellow to brown. Khaodee et al. (2014) [57] reported that cyanidin forms colour complexes with Cu, lead (Pb), and aluminium (Al), showing a shift from purple to blue, as well as with Fe, which shows a shift from purplish pink to blue. Several other studies have also been conducted related to the development of natural pigments as metal detectors. In Table 2, data from some studies related to this topic are shown, along with the colour changes and detection limits that can be visually observed.
Computational methods can be used to efficiently study the interaction of metals with natural pigments. Direct experiments in the laboratory may take more time, cost more, and produce more waste. Computational studies are more efficient and able to determine the interactions that occur and other parameters [64,65]. Computational calculations will be very useful in predicting the possibility of colour shifts because changes in the ground-state geometry will affect the absorption peak shift [66]. A small change in the basic geometry of a molecule can trigger larger shape and chemical changes, including optical properties, distribution of electrons, and electronic structure, that shift the wavelengths of light it absorbs, which is directly related to the colour change that can be observed by the naked eye [39,67,68]. Until now, there has been no review that discusses this field, even though this preliminary research is very important and useful. Therefore, we provide an overview of the computational methods that can be utilised to study metal–natural pigment interactions. We also discuss the validation of these computational methods with experimental results.

2. Computational Methods for Studying Metal–Pigment Interactions in General

Computational methods have been used to study metal–natural pigment interactions for several years. During this time, they have improved and become increasingly well known. Studying the bonds and interactions between metals and pigments using computational methods is an important step before laboratory research; it can lead to more targeted studies that require fewer resources, costs, and energy. In addition, computational methods can overcome several limitations, and consistent, detailed experimental information at the structural and energetic levels can be obtained [69].
Due to its various advantages, computational study is now widely considered for use in various fields. Various calculations and predictions can be carried out using computational methods, for example in studying molecular system energy [70,71], gibbs free energy [64], binding energy [72], bond length and bond angle [73,74], natural bond orbitals [65], HOMO and LUMO [71,75,76], stability of the complexes [64,77], optical properties [75], and molecular dynamics studies [78]. Apart from that, this computational method can also be used to predict the electronic structure of a system and, thus, the distribution of electrons in the system [79,80].
In the development of natural pigments as metal detectors, one of the parameters that is directly related to metal detector development is excitation energy. Excitation energy is related to colour changes that may occur before and after the pigment interacts with the metal. The colour that an object possesses is determined by the wavelengths of light that are absorbed and reflected by that object, where the wavelength will depend on the structure of the molecules. When bonds are formed between natural pigments and metals, the electron configuration will undergo changes that also affect molecular orbits and energy levels. This is what will affect the occurrence of wavelength shifts and observed colour changes. In developing a metal detector from natural pigments, the excitation energy can be used to predict whether there is a change in colour and how the colour changes [71,81]. However, before studying the excitation energy of a molecule, studying the interactions that occur through analysing bond energy, total energy, affinity, and ground-state energy calculations is also important to do as a preliminary study to find out whether the bond will be formed between natural pigments and metals because, to produce colour changes, natural pigments must be able to bind to metal [17,82]. Apart from that, various other calculations such as the stability of the complexes, optical properties, and the geometry of the complexes formed can also be very helpful in the development of metal detectors from natural pigments [64,75,83].
Various computational methods and approaches are available to study the parameters mentioned above, and each method used and each approach applied will impact the results of the study carried out. Therefore, it is important to know what method is most appropriate and what kind of approach can be taken. It should be noted that the selection of an appropriate method will also be influenced by various factors and must be based on the specific needs and goals of the study, including what parameters are to be studied. In addition, it is also important to find an optimal combination of basis sets, methods, and approaches through multi-level approaches so that the results obtained have a balance between efficiency, accuracy, and robustness [84,85].
As mentioned earlier, there are many factors that can influence the outcome of computational calculations, which must also be taken into consideration. The main thing that must be considered is, of course, the molecule to be investigated [86], but apart from that, there are several other factors, such as the presence of other molecules in the system like co-pigments, pH [78,86], temperature [87,88], and also the solvent effect [89,90]. Solvents are capable of causing the occurrence of spectral shifts and, as explained earlier, the wavelength of the complex formed becomes a direct parameter that can be studied in the development of metal detectors from natural pigments. Solvents can cause changes in the ground and the excited state energy of molecules, so they can cause shifts in the absorption or emission spectrum of molecules. Thus, it can be concluded that the optimal wavelength in one solvent will be different than in another, as well as in the gas and solution phases [91]. The interactions between the molecules and the solvent can be studied using various methods. Explicit evaluation can be carried out using various methods, including free energy methods such as thermodynamic integration or metadynamics, where molecular dynamic-based sampling techniques and Monte Carlo methods can be carried out. This method can provide more detailed data compared to the implicit method, but the computational cost is relatively high [92]. On the other hand, implicit evaluation can also be carried out, where, in this model, the solvent is treated as a structureless continuum with certain dielectric and interface properties. In this approach, the contributions to the solvation-free energy of several models, such as conductor-like screening (COSMO) and polarizable continuum models (PCMs), can be applied [92,93]. Usually, method selection is based on a balance between calculation speed, cost, and expected accuracy. Data from various methods is also need to be taken into consideration in selecting the appropriate method [91]. These two approaches can also be combined to obtain more optimal results with better efficiency [94,95,96].

3. Quantum Mechanical Methods to Study Metal–Pigment Interactions

Quantum mechanical methods describe the electronic structure of a system to estimate the distribution of electrons in the system. These methods also allow one to evaluate the ability of pigments and metals to form complexes, including the possibility of bond formation, the structures formed, bond distances, bond energies, and the stability of the complex formed [80,97]. The major quantum mechanical approaches are the ab initio method, DFT, and semiempirical methods. Several studies have used these methods to study the bonds between metals and pigments. Below, we discuss the quantum mechanical methods in greater detail. The study flow of computational chemistry using quantum mechanical methods to study the interaction between natural pigments and metals can be seen in Figure 2.

3.1. Ab Initio Methods

Ab initio methods involve calculations that are derived directly from theoretical principles without requiring experimental data. There are several systematic approaches, including solving differential equations or Born–Oppenheimer approximation. One of the most widely used approaches to studying metal–pigment interactions is the HF calculation with the central field approximation method [98]. Table 3 presents studies that have used ab initio methods to study metal–pigment interactions.
Linnanto and Korppi (2004) [101] used the HF/6-31G* approach to calculate fully optimised structures and to calculate atomic charges of methyl bacteriochlorophyllides a, b, g, and h, which is a pigment responsible for photosynthesis in bacteria [105], and the magnesium (Mg)–bacteriochlorin complex. The authors also carried out optimisation using (DFT/B3LYP/6-31G*). In addition, the CIS/6-31G* and CIS/6-311G** configurations, as well as the HF/6-31G* and HF/6-311G** methods, have also been utilised to calculate the associated spectroscopic transition energy of the chromophores. The HF approach provides the structure in a vacuum; the four coordinated Mg atoms of bacteriochlorin are positioned centrally, nearly in the plane of the ring. The outcomes were consistent with calculations by other methods (e.g., DFT/B3LYP/6-31G*) and some of the earlier X-ray structure study results [106,107,108]. These results are consistent with semiempirical [101] and time-dependent DFT [109] calculations. The ab initio and DFT calculations produced slight differences in energy transitions compared with the experimental results. The bond length and coordinates have a strong influence on transition energy calculations, and because there are slight differences in the results from the methods used, the transition energy results are slightly different. The electronic transition calculations using ab initio methods are systematically higher than DFT calculations. Moreover, DFT with B3LYP and ZINDO/S transition energies produced the finest correlation coefficient between the calculated and experimentally determined transition energies [101].
Cao et al. (2014) [103] used the ab initio method with the HF/6-31G* level of theory and the Gaussian 09 program for full geometry optimisation and to determine whether there is an influence of additional functional groups in the induction of atomic charge distribution and the influence of methylation position on luteolin and quercetin. Quercetin and luteolin are flavonoids, which are pigments with a phenolic structure that can be found in vegetables, fruits, and medicinal plants [80,110]. Before carrying out calculations, the authors compared luteolin and quercetin regarding functional groups, the type of bond that occurs, and the possibility of electron movement. Quercetin has an additional 3-OH group in its structure, while luteolin does not; the three -OH groups will donate electrons in quercetin with a delocalised π bond. With these considerations, the authors carried out calculations on a single molecule to determine whether the additional -OH group in quercetin induces a different atomic charge distribution on the two oxygen atoms in the luteolin and quercetin catechol groups. It is important to perform this initial assessment before the calculations to choose the most appropriate approach and get more information from the calculations. Following the preliminary analysis, the authors calculated the partial charges of the meta-O and para-O structures in luteolin and quercetin and noted that the meta-O structure is more nucleophilic because both luteolin and quercetin have more negative charges in their meta-O structure than in their para-O structure. Furthermore, the equivalent hydrogen in meta-OH is more attracted to the oxygen atom than it is in para-OH, meaning that meta-OH would have a larger reaction capacity for methylation than the para-OH structure [103].
Apart from being used to study the bonds that occur between pigments and metals, ab initio methods can also be used to see the effect of hydration on metal–pigment interactions. Hydration can impact the physicochemical properties of metals, pigments, and metal–pigment complexes; the interactions that occur; and the reactivity and stability of various systems. In addition, water molecules can impact the compounds and the complexes formed in diverse ways [39,111].
Overall, ab initio methods can be used to study metal–pigment interactions, including to optimise complex structures, analyse the dynamic properties, solvation capabilities, and hydration behaviour of a compound, band structures, polarisabilities, and atomic charges, and calculate energies. The HF method is widely used along with the 6-31G, 6-31G*, 6-31G**, and 6-31G(d) basis sets, both single basis sets and combinations of several basis sets. Deeper development of the method with various combinations of other approaches could be carried out. However, ab initio methods are not the main choice to study metal–natural pigment interactions due to their relatively poor level of accuracy. Although ab initio methods have not been used to study metal–pigment interactions in depth, they have been used to study the interactions between heavy metals and other compounds [112,113,114]. Thus, the use of this method to study heavy metal–natural pigment interactions can still be developed to determine its uses, advantages, and disadvantages for studying metal–pigment complexes.

3.2. Density Functional Theory (DFT)

DFT has been used since 1990 to calculate the electronic structure of various molecules; it is considered a powerful method for simulating chemical systems [115,116]. In addition, it has a good level of efficiency and accuracy [115]. DFT can be used to calculate atomisation energy or binding energy [17,103,116], bond lengths and angles [73,115], electronic configuration [116], electron affinity [115], ionisation potential [65,115], heat of formation [115], several types of nonbonding interactions [115], and geometric optimisation of a complex, especially for a system involving transition metals, which needs to consider electron correlation. In this method, the system of electrons that undergo interactions is mapped specifically into an effective non-interacting system with a similar total density [116,117]. Table 4 presents studies that have used DFT to study metal–pigment interactions.
DFT calculations are the most widely used method, including to study metal–pigment interactions [65,124,126,127,128,129,130,131]. This approach continues to be developed functionally, and to date, there are many methods with various applications that can be used. Sousa et al. (2007) [115] collected data related to the various DFT functionals available. They found that B3LYP was the most widely used DFT functional, with a frequency of 80% from 1990 to 2006. BLYP, B3PW9, BP86, and M05-2X were also relatively widely used functionals (3% each). B3LYP is widely used because it has several advantages: it is considered an accurate method with wide coverage and relatively affordable calculation costs. However, it has several shortcomings, such as the possibility of a decrease in the level of accuracy as the size of the system being investigated increases, as well as an underestimation of weak interactions [132]. This method can still be developed further to achieve higher accuracy and wider use because each available functional is very flexible, making it possible to combine one or more exchange functionals and correlation functionals [115,129,130,131,133].
DFT has been widely applied to study interactions between natural pigments and non-heavy metals (such as Al, Co, Cr, Cu, Fe, Mg, Mn, Ni, palladium [Pd], and Zn) and heavy metals (such as As, Cd, Hg, Pb, and tin [Sn]). Sun et al. (2008) [72] studied the chelation process between quercetin and Cr using DFT with the B3LYP functional combined with the LANL2DZ pseudopotential and the 6-31+G (d) basis set to determine the appropriate chelation site and to produce a stable metal–pigment complex. They calculated the bond energy and total energy and performed natural bond orbital analysis to analyse the charge transfer between pigment and metal. They found that Cr(III) tends to form a complex with quercetin at the 5-hydroxy-4-keto site. In addition, deprotonated quercetin attracts Cr better than its natural form. The authors also compared the computational results with direct synthesis and analysis of the bond characteristics using infrared spectroscopy and ultraviolet–visible (UV–Vis) spectroscopy. The experimental and computational experiments led to the same conclusions. Of the three accessible chelating sites, Cr(III) ions tended to chelate at the 5-hydroxy-4-keto site, where deprotonated quercetin has a higher chelating capacity than natural quercetin. The UV–Vis spectrum showed that deprotonated 5-OH coordinates the Cr(III) ion to form a complex. This study showed that a computational method provides accurate initial testing that is less expensive and relatively more efficient than direct experiments [72].
Guo et al. (2008) [117] studied the interaction between the butein anion and Mg, Cr, Fe, and Cu cations. Butein is a chalcone compound that can be found in various plants, such as dahlia, coreopsis, and lacquer tree [134]. Gou et al. (2008) used the B3LYP functional combined with the LANL2DZ pseudopotential and the 6-31+G (d) basis set, similarly to the approach employed by Sun et al. (2008) [72]. They performed metal–butein structure optimisation without any symmetry constraints, followed by a single-point calculation utilising a different basis set (6-311++G(d,p)). This method produced binding energies and total energy at the B3LYP/6-311++G(d,p) level to establish the most appropriate coordination and binding sites for a stable complex. Apart from that, the authors calculated natural bond orbitals with the same functional (B3LYP) but with a different basis set (6-311+G(d,p)) to study the charge transfer between butein and the metal cation. They found that the oxygen atom in positions 2′ and 9 was the most preferred. The calculations revealed that the strength of metal against metal is Mg2+ < Fe2+ < Cr2+ < Cu2+ in both the gas and liquid phases [117].
Leopoldini et al. (2006) [76] use the B3LYP/6-31G*/LANL2DZ level of DFT for geometric and electronic optimisation, followed by natural bond orbital analysis for better characterisation of Fe(II)–quercetin complexes. They performed single-point calculations with the 6-311++G** basis set for the nonmetal atoms to refine the electronic energies. Both neutral and deprotonated quercetin can form stable complexes with Fe, with an optimal binding energy at a metal-to-ligand ratio of 1:2. The most favoured bond positions are the oxygen atoms on carbons 3 and 4, and 5 and 4 with high binding energy values, which indicate that quercetin and Fe(II) are capable of forming strong complexes. In another study, the authors compared the calculations with experimental results based on a UV–Vis spectrum obtained computationally (B3LYP/6-31G*) with experimental results in methanol [135] and dimethylsulfoxide [136]. The authors identified absorption bands in the 370–390 nm region.
Eno et al. (2023) [65] studied the interaction between several metal ions (sodium [Na+], potassium [K+], Mg2+, Ca2+, and Al3+) with quercetin at the B3LYP/6-31+G(d) level of DFT. They examined reactivity, stability, origin of interaction, and the use of metal–quercetin complexes as potential antioxidants. The metal–pigment interactions were purely electrostatic. The K+–quercetin complex was the most stable, and the Ca2+–quercetin complex had the lowest values in the neutral form and was the most reactive. Based on its reactivity as an electron donor, in neutral conditions, the K+–quercetin complex is the most reactive. It has the lowest ionisation potential (382.78 kcal/mol). The Ca2+–quercetin and Mg2+–quercetin complexes have the lowest reactivity as electron donors, with an ionisation potential of 395.33 kcal/mol. In the deprotonated condition, the Mg2+–quercetin and Ca2+–quercetin complexes have a lower ionisation potential (363.95 kcal/mol), which indicates that in the deprotonated condition, the Ca2+–quercetin and Mg2+–quercetin complexes have a highest electron transfer reactivity [65].
Several studies have also used the M05-2X functional, which can be combined with several different basis sets to study the interaction between natural pigments and non-heavy metals. Malacaria et al. (2022) [64] analysed the complex between luteolin and Al, Fe, and Cu. They evaluated the properties and optimised geometry with the M05-2X hybrid function and the 6-31+G(d) basis set for carbon (C), hydrogen (H), oxygen (O), and Al atoms, and the relativistic compact Stuttgart/Dresden effective core potential with its related split valence for Cu and Fe atoms. They performed optimisation in a solvent by applying the solvent-based density model with a dielectric constant of 78.0. The authors calculated the free energy from the complexation reaction of the substitution of water molecules in the hexaaqua complex to compare the bonding abilities of various metals with pigments. Then, they refined the energy by single-point calculations using the same functional but with a larger basis set (6-311++G(d,p)). The computational approach shows that the Cu ion does not prefer specific sites in the pigment, while the Al and Fe atoms appear to prefer to bind to the 4-carbonyl-5-hydroxyl site on the A and C rings of the ligand [64]. Corrente et al. (2021) [17] also conducted the same research while comparing the results obtained with experimental experiments with the aid of 1HNMR (Proton Nuclear Magnetic Resonance) and 13CNMR (Carbon-13 Nuclear Magnetic Resonance) spectroscopy obtained during the titration experiment. The computational and experimental analysis provided consistent results. Quercetin tends to bind three cations at different coordination sites. Fe and Al are known to have no tendency to bind to certain sites from the three complexation sites on quercetin, whereas for Cu, sites 4–5 can be excluded [17].
Tanui et al. (2022) [123] used DFT with another approach to study the complex formed between morin, which is a yellowish pigment obtained from Chlorophora tinctoria wood [137], or quercetin, and Fe and Cu ions. They used the M06-2X/def2-SVP level of DFT to evaluate the stability of the complex formed. They found that morin has better stability than quercetin. It is also known that six-membered structures for morin and quercetin complexes with metals have slightly better stability compared with their respective five-membered counterparts. The two pigment complexes with Fe showed better stability compared with the complexes with Cu (based on a much lower energy) [123]. The M06-2X/def2-SVP functional can also be used to investigate the effect of the solvent by examining possible interactions between the solvent and the compound. Mollaamin et al. (2020) [138] investigated van der Waals density functions to examine the influence of solvents on delphinidin complexes with several non-heavy metals (Al, gallium (Ga), Cr, Fe, and Mg). They carried out accurate calculations by using m062x, m06-L, and m06-HF. These three methods have appropriate correspondence for non-bonded calculations between the compounds and the solvents used [138].
DFT can also be used to study heavy metal–natural pigment interactions (Table 5). Cornard et al. (2005) [135] studied the interaction between Pb(II) ions and quercetin in methanol, with the B3LYP/6-31G(d,p) level and LANL2DZ for Pb to optimise the geometry of the complex formed. Based on this research, three potential locations are associated with the possibility of bond formation. The highest complexation strength for the Pb(II) ion is in the catechol functional group, where, after complexation, the Pb atom is coordinated to the catechol group, which is on the same plane as quercetin, and the ligand stays completely planar. These results are consistent with the experimental electronic spectrum and complex vibrations. Based on the theoretical calculation results, there is a bathochromic shift of the wavelength band in the UV–Vis spectrum when complexation occurs. This is thought to occur due to the transfer of the ligand charge to the metal. By comparing the experimental and theoretical results, it is also possible to compare the chelating strength of Pb ions from three potential sites in quercetin, and the results show agreement that the catechol group has the best complex formation strength [135].
Pratiwi et al. (2020) [139] also used the B3LYP level and the 6-31G and LANL2DZ basis sets to study interactions and complex formation and to calculate the binding constants of three cationic porphyrin compounds, namely meso-tetrakis (N-methylpyridinium-4-yl) porphyrin (TMPyP), meso-tetrakis (1,3-dimethylimidazolium-2-yl) porphyrin (TDMImP), and meso-tetrakis (1,2-dimethylpyrazolium-4-yl) porphyrin (TDMPzP) with five heavy metals (Pb, Cd, Hg, Sn, and As). The authors discovered that the complexes between porphyrin cations and the five heavy metals occur spontaneously, as indicated by the lower energy compared with the free base form. The metal–TDMPzP complexes are less electrophilic and have greater chemical hardness compared with the metal–TMPyP or metal–TDMImP complexes. Hence, the metal–TDMPzP complexes are most stable, especially the complex with Cd, which has a lower energy level and electronic chemical potential, lower electrophilicity, and a higher binding constant as well as chemical hardness. On the other hand, TMPyP has a higher energy compared with the other two porphyrin cations, which shows that this compound has the lowest stability. The stability of the heavy metal–porphyrin complexes is also influenced by the ionic radii of the heavy metals, where it is known that Cd with a radius of 95 and Sn with a radius of 69 can enter the core of the porphyrin and form a complex, while Pb and Hg, which have larger radii (118 and 102, respectively), are unable to enter the porphyrin core because they are too large, thus forming the metalloporphyrin located above. On the other hand, As has a smaller radius than all other metals but tends to move out of the porphyrin plane, a phenomenon that is estimated to occur because the valence electrons of As ions are not suitable for interacting with the nitrogen in the porphyrin core. This result is also in line with experimental results, where it is known that TDMPzP will form a stabler metalloporphyrin complex when reacting with Cd2+ ions, with a binding constant of 5.1 × 107 M−1 [139].
Overall, DFT methods can be used to study metal–natural pigment interactions with a fairly good level of accuracy. Among the many approaches that can be used, the B3LYP level has been the main choice to study natural pigment interactions with non-heavy metals and heavy metals. This is thought to be related to its good accuracy and wide coverage at an inexpensive cost. The LANL2DZ ECP and basis set are also often chosen and combined with other basis sets. A combination of various basis sets is often used to obtain better results and to complement the shortcomings of each basis set. Future development and research related to combinations in computational systems using DFT methods are very promising.

3.3. Semiempirical Methods to Study Metal–Pigment Interactions

Quantum mechanical calculations using semiempirical methods are generally similar to ab initio methods. These two methods can complement each other and overcome each other’s shortcomings. Compared with the ab initio approach, semiempirical calculations are relatively faster, but the results will depend greatly on the availability of parameters that correspond to the molecules analysed in the database. If the molecules are significantly different from the molecules in the database, then the calculations may deviate from the actual results. When carrying out inorganic analysis, semiempirical methods are used relatively rarely because the potential for errors that occur in the analysis tends to be higher, especially when transition elements are involved. Therefore, it is a good idea to compare semiempirical calculations with calculations using other methods or with experimental results [98,141]. Table 6 presents research that has used semiempirical methods to assess metal–pigment interactions.
Linnanto and Korppi (2004) [101] used a semiempirical method—at the ZINDO/S CIS (40,40) or (45,45) level in the ArgusLab (version 2.0.0) software—to calculate the transition energies and oscillation strengths of the Mg–bacteriochlorin complex and methyl bacteriochlorophyllides a, b, g, and h. They also combined the ZINDO/S CIS (40,40) level with the simple self-consistent reaction field (SCRF) method. The authors used a semiempirical method to assist in the estimation of the corresponding spectroscopic transition energy values of the chromophore and also studied the existence of dark electronic states in the system. Based on this research, pigment complexation will induce a dark electronic state below the main Soret transition, where this condition may have an important role in the energy transfer process [101].
Linnanto and Korppi (2004) [141] also used the semiempirical PM5 method to study orbitals on chlorophylls and bacteriochlorophylls. They optimised the structure of Mg–bacteriochlorin; Mg–chlorin; Mg–porphin; mesochlorophyll a; chlorophylls a, b, c1, c2, c3, and d; and bacteriochlorophylls a, b, c, d, e, f, g, and h, all with homologous structures. To ensure that the results were truly accurate, the authors compared the semiempirical PM5 method with several other quantum mechanical methods. The experimental, ab initio, and density functional results were suitable. Based on calculations and an experimental study with X-rays, all pigments are predicted to have a planar structure of porphyrin with four Mg atoms coordinated in a position almost in the centre of the porphyrin ring plane. More clearly, Mg is in the centre for Mg–bacteriochlorin, Mg–chlorin, and Mg–porphin and is slightly shifted out of the porphyrin plane for the chlorophylls and bacteriochlorophylls. The authors also used the PM5 method to investigate minimum energy and atomic charges. They found only a slight difference between the minimum energy obtained with PM5 and other methods in the form of B3LYP/6-31G* or HF/6-31G*. These results indicate that the efficiency of PM5 is quite good and computationally efficient, considering that the HF and DFT methods are more intensive [141].
Cornard and Merlin (2002) [74] conducted a semiempirical study using the AM1 Hamiltonian in the Hyperchem (version 5.0) program to evaluate structural modifications caused by the ligand via chelation of one and then two Al ions. Previously, quercetin’s shape and the complex it forms with Al had also been optimised by using this similar methodology. In an isolated state, the molecule adopts a staggered conformation with a O (O1–C2–C1′–C6′) angle of 26.7° [74,148]. Cornard and Merlin (2002) used the same semiempirical method to determine the geometric parameters of quercetin and two complexes that are generated after energy minimisation with the AM1 method and several stoichiometric ratios in different solvents (pure methanol, methanol + AcO, and methanol + water + HCl). The authors also generated experimental electronic and vibrational spectra and compared these results with the calculations. They found that the 3-hydroxychromone group is the initial site implicated in the complex formation process in acidic media and pure methanol, whereas the catechol group has maximum chelating power in alkaline media. It is known that for large amounts of Al(III), with an Al(III)-to-quercetin ratio greater than 0.5, the same complex [Al2QR4]3+ is formed in methanol and alkaline media. However, in acidic media, the ortho-dihydroxyl group is never involved in complexation with Al(III). The authors also carried out a spectroscopic study to validate the calculations. The UV–Vis spectrophotometry results showed wavelength conformity with the calculation results. They also validated the calculations with Raman spectroscopy. The complex exhibits important spectral changes, especially in the range of 1500–1700 cm−1, and the mass formed by the two broad bands at 1612 and 1576 cm−1 in the quercetin spectrum is profoundly affected, indicating that the structure of both rings is modified by chelation [74].
These studies have shown that semiempirical methods can be used to study metal–natural pigment complexes, including how the complex forms, modifications that occur due to the chelation process, intermolecular behaviour, oscillation strengths, and corresponding spectroscopic transition energy values. In addition, semiempirical methods can be used to optimise structures and solvation. The calculations are consistent with several experimental results, but in terms of accuracy, the ab initio and DFT methods are still superior. Semiempirical methods have not yet been used to study heavy metal–natural pigment interactions in greater depth, even though they have been used to study interactions between heavy metals and other components [149,150].
For quantum mechanics, DFT methods are selected most often to study metal–pigment interactions. They have several advantages over ab initio and semiempirical methods, including a good level of efficiency and accuracy with a high level of implementation flexibility [115,129]. However, the suitability of DFT can also differ depending on the expected objectives. Linnanto and Korppi (2004) compared ab initio, DFT, and semiempirical methods to predict the structure of Mg–bacteriochlorin complexes. They found that the calculations using DFT with B3LYP/6-31G* had the smallest difference in transition energy compared with the experimental results. However, this study also showed that the size of the model system chosen and the environmental conditions for calculations greatly influence the reliability of the results obtained. Indeed, the pigment environment itself can alter the geometry of the chromophore and the electron density properties, which will change the energy level position. Therefore, it can be concluded that, apart from choosing the method that best suits the expected goals, it is crucial to choose the appropriate calculation environment conditions and system size [101].

4. Molecular Dynamics Methods to Study Metal–Pigment Interactions

Molecular dynamics is a computational method that can be used in simulation to predict the position and movement of each atom in the system at each point in time based on a general model of physics [98,151]. In this simulation, the future speed and position of the atom are estimated based on its current speed and position by considering thermal motion [78]. As an illustration, the steps that can be carried out in molecular dynamics simulations to study interactions between natural pigments and metals and the influence of different parameters can be seen in Figure 3.
Recently, molecular dynamics has become increasingly popular because it is powerful and easy to access [152]. Its wide use is also related to its capabilities; it includes tests of various kinds of molecular behaviour as well as chemical systems as a function of time [98]. Molecular dynamics simulations can model important biomolecular processes, including conformational changes and conformation theory, binding of ligand molecules, and thermodynamic parameters, revealing the positions of all atoms at femtosecond temporal resolution. These simulations also allow researchers to describe the structure, thermodynamic properties, kinetic energy, surface potential energy, and dynamic properties of a system [78,98,152]. Table 7 presents the details of studies that have used molecular dynamics simulations to study metal–pigment interactions.
Moradi et al. (2021) [158] used molecular dynamics simulations to study the adsorption behaviour of the quercetin from Satureja hortensis L. extract (SHE) complexed with Zn on metal surfaces. SHE is derived from S. hortensis (synonym Satureja postii Arzn.), an aromatic plant from the Lamiaceae family [159]. The authors performed these simulations to understand the adsorption behaviour of Zn cation complexes with pigments on the surface of an Fe substrate. The authors performed modelling by using the Forcite module of the Material Studio software in a simulation box (2.73 × 2.73 × 3.31 nm). They attached a bottom Fe plate and an upper solvent layer, containing one complex and 800 water molecules, at a speed of 1500 ps in the NVT ensemble using a universal force field with a timeline of 1 fs and a temperature of 298 K. They observed interactions between the Zn(II)–quercetin complex and the substrate (Fe) in the water phase. The authors considered non-bonding interactions in the system to be an ‘atom-based simulation’ and electrostatic interactions adapted as an ‘ewald simulation’ approach. By balancing temperature and energy, the entire system under investigation reaches a state of equilibrium. The authors reported an absorption energy (Eads) of −42.4 kJ mol−1 for Zn(II)–quercetin, indicating that there is good adsorption strength between the adsorbate and the metal substrate. Eads will greatly determine the strength of adsorption on the surface of the substrate: an increasingly negative value indicates a stronger interaction as well as stable and spontaneous absorption. The authors compared the theoretical results with the experimental results. The calculations could satisfactorily explain the experimental findings [160,161]. Based on the field emission scanning electron microscopy results, it is known that the addition of an inhibitor can coat Fe to prevent corrosion. Apart from that, the Fourier-transform infrared spectrum showed additional peaks that indicate the formation of Fe-O bonds from the complex and Zn-O bonds from the Zn hydroxide and its adsorption on metal surfaces. This is also in line with the UV–Vis spectroscopic results: a comparison of adsorption rates shows that the complex will be more adsorbed on the steel surface and form an inhibitory layer than dispersed in the electrolyte with and without the addition of inhibitors [158].
Moin and Hofer (2014) [155] conducted ab initio quantum mechanical charge field–molecular dynamics simulations using the HF method along with 6-31G** basis sets to study structural properties and dynamic properties and to observe the hydration behaviour of porphyrin and Mg–porphyrin complexes. Porphyrin is one of the most prevalent types of biochromes; this macrocyclic compound is the result of haem breakdown [139,162]. This simulation was done for Mg–porphyrin complexes with porphyrins inserted into a cubic box with a side length of 39.28 Å with B2000 explicit water molecules utilising periodic boundary conditions. The system was run in the NPT (isothermal–isobaric) ensemble using the Velocity–Verlet algorithm and at 298.15 K. The charge of the water molecules was adopted from the SPC/E water model. Based on the binding free energy and surface area exposed to solvents, these two confirmations have different hydration behaviours. Complexation of Mg ions with porphyrin will change the hydration pattern significantly, where the Mg–porphyrin complex has much better stability, as indicated by the stability of the penta-coordinated axial water molecules throughout the simulation. Compared with uncomplexed porphyrin, the Mg–porphyrin complex shows differences in molecular coordination that produce different H-bonding patterns. This is also supported by the vibrational power spectrum evaluated for both solutions obtained via the Fourier-transform velocity autocorrelation function observed in the quantum mechanical charge field–molecular dynamics simulation [155].
Singh et al., 2019 [156] also carried out molecular dynamics simulations with the AMBER14SB force field for standard groups and the GAFF force field for the remaining atoms in studying metallophilic interactions between closed-shell metal ions. In this research, a molecular dynamics study was carried out to evaluate the possibility of increasing or decreasing internuclear distance upon low-energy conformational changes in Ag in the middle of the porphyrin ring. To evaluate changes in the internuclear distance between Ag molecules, an analysis was carried out on the flexibility of the complex (Porphyrin-Ag) in vacuum for 100 ns. The root-mean-square deviation of the system was calculated during the simulation by using the minimised initial structure as a reference. The results obtained show that the trajectory attains relatively stable RMSD after approximately the first 20 ns, with an average value of 0.71 (0.36). In addition, the results in this study also show estimates of internuclear distances Ag···Ag, i.e., in the range 3.0–4.2 in the first cluster (average: 3.87(0.59)) and 3.9–4.9 in the second (average: 4.37(0.51)), where the results of molecular dynamics simulation are perfectly compatible with those predicted by EPR spectroscopy [156].
Huamán et al. (2021) [157] carried out molecular dynamics simulations of cyanidin-3-glucoside and bixin analysed in an isolated titanium dioxide (TiO2) dye/cluster system. Cyanidin is an anthocyanin that is responsible for the red and purple colours in berries, red sweet potatoes, and purple corn [163]. During calculations, the authors applied charge balance throughout the molecular system. They used the NVT ensemble at 300 K to maintain the temperature during the simulation under environmental conditions with mass conservation. The authors determined the molecular system’s energy by solving for potential and kinetic energy and for the binding energy of the dye–TiO2 complex. They subtracted the total energy of each component from the total energy of the complex. In addition, to study the dye–TiO2 electrode interaction, the authors employed complex molecular dynamics simulations using Reax FF. They found that chemisorption is carried out through a successive deprotonation process, and the formation of a bond between Ti-O and cyanidin-3-glucoside occurs by an OH retaining group and by a COOH retaining group for the bond between Ti-O and bixin, where bixin also shows greater affinity with the overly folded TiO2. Bixin also shows a tendency for dimer formation due to intermolecular π–π interactions that facilitate aggregation at the surface. The binding between Ti-O and cyanidin-3-glucoside is related to the knowledge that cyanidin-3-glucoside sensitisation is a monodentate anchor via the C4′ hydroxyl group of the benzene–diol molecule, where this interaction causes hydrogen bonding and passivation of TiO2 with protons from the glycosidic part of the molecule. The authors validated the calculations with UV–Vis spectroscopy and Fourier-transform infrared spectroscopy. They noted a broadening of the UV–Vis absorbance spectrum, evidence that electronic transfer occurs in the context of the formation of Ti-O bonds. Furthermore, Fourier-transform infrared spectroscopy showed that after sensitisation to TiO2, the peak associated with free OH bonds (3600 cm−1) disappear, forming a single band associated with hydrogen bonds, which includes inter- and intramolecular interactions, then C-H stretching from benzene becomes weaker and shifts to the red (2910 cm−1). Apart from that, one of the strongest pieces of evidence of an interaction is the presence of the Ti-O absorption band at 559 cm−1 [157,164].
Asadi et al. (2020) [120] carried out Monte Carlo/molecular dynamics simulations to study the interaction between the Zn–luteolin complex and Fe. Using the Adsorption Locator module in the Materials Studio software, the authors initially performed Monte Carlo simulation, using the Monte Carlo cells of the complexes as the initial configurations for the aqueous-phase molecular dynamics simulation. Next, the authors performed the molecular dynamics simulation in the NVT ensemble with a time step of 1 fs, a controlled temperature of 298 K, and in the COMPASS force field. The results showed that the Zn(luteolin)n (n = 2–5) complexes adsorb on the Fe surface. However, stability near the surface can only be achieved by the Zn(luteolin)2 and Zn(luteolin)3 complexes because of the tetrahedral-like geometry of large complexes. Subsequent investigation revealed that chemicals adsorb with parallel alignment to the surface. The authors calculated the adsorption energies: −691.78 kcal mol−1 for Zn(luteolin)2, −862.68 kcal mol−1 for Zn(luteolin)3, −1000.19 kcal mol−1 for Zn(luteolin) 4, and −1049.55 kcal mol−1 for Zn(luteolin)5. This negative value indicates that the adsorption that occurs is stable and spontaneous. The adsorption energy increases as the number of luteolin ligands bound to the Zn atom increases. The field emission scanning electron microscopy and atomic force microscopy results show that after mixing, there is a uniform film of a combination of Zn2+ cations and organic compounds of lemon balm extract in the specimen formed. Apart from that, grazing incidence X-ray diffraction showed that the formation of lower corrosion products (Fe oxide/hydroxide) occurs on the surface, a phenomenon that can be connected to the protective inhibitory film generation on the steel surface with the addition of an inhibitor, in this case a compound complexed with Zn [120].
Molecular dynamics simulations can be widely used to study metal–pigment interactions, including the simultaneous interactions between several metals and pigments. Moreover, various conditions can be simulated to evaluate their impact on the experimental results. Because molecular dynamics considers Newton’s laws of motion—where the spatial position of each atom as a function of time is predicted and then used to calculate the force of each atom—it can be used to predict how a molecular system will experience movement over time based on general models of physics [151,152]. The structure and properties of the system can also be described dynamically and thermodynamically. In simulations carried out in a solvent system, the solute and the simulated environment in the system can be controlled, and their position and movement can be known at each point in time; this is very difficult, or perhaps even impossible, to perform experimentally [165,166]. In addition, the components involved in the simulation are known in detail, such as the initial conformation of the compound used, the ligands that are bound, and what other molecules are present in the simulation systems. In this way, various impacts or influences from variations in conditions can also be identified in the simulation. Another advantage of molecular dynamics simulations is their ability to reveal the dynamic behaviour of water molecules and various ions in the system, and the accuracy of the modelled structure can also be checked and even refined. Based on Table 6, molecular dynamics simulations are consistent with experimental results. Therefore, the use of molecular dynamics simulations to study metal–pigment interactions is very useful and promising [152,167]. Although the application of molecular dynamics simulations is increasingly widespread, they have not been used to study heavy metal–natural pigment interactions. However, they have been used to analyse interactions between heavy metals and other components [168,169,170].

5. Conclusions

Quantum mechanics, including ab initio, DFT, and semiempirical methods, as well as molecular dynamics, have specific advantages and can be combined to complement each other. Among the available quantum mechanical methods, DFT with the B3LYP functional and the LANL2DZ ECP and basis set have been the most widely used due to their good accuracy and efficiency and relatively low cost. Molecular dynamics simulations allow one to study interactions in greater depth in various system conditions and to determine changes or predict dynamic movements over time. These simulations have been widely validated with various experimental studies. Considering their time, cost, and reliability, computational methods for preliminary studies and even for further studies of a system are very promising. Additional development and analysis of each method is still required, especially to determine the optimal combination of methods, approaches, functions, and basis sets to complement each other and increase accuracy and efficiency. Additional experiments are needed to find the best combination to study the bonds and interactions between heavy and non-heavy metals and natural pigments. Matters related to spin states and dispersion effects of metals on the development of metal indicators could also be a good part to evaluate in future research as the amount is very limited.

Author Contributions

Conceptualisation, A.N.H. and S.M.; methodology, A.N.H.; data curation, G.J.M.; writing—original draft preparation, G.J.M.; writing—review and editing, A.N.H. and S.M.; visualisation, G.J.M.; supervision, A.N.H. and S.M.; project administration, A.N.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mendel, R.R.; Smith, A.G.; Marquet, A.; Warren, M.J. Metal and Cofactor Insertion. Nat. Prod. Rep. 2007, 24, 963. [Google Scholar] [CrossRef] [PubMed]
  2. Permyakov, E.A. Metal Binding Proteins. Encyclopedia 2021, 1, 261–292. [Google Scholar] [CrossRef]
  3. Fu, F.; Wang, Q. Removal of Heavy Metal Ions from Wastewaters: A Review. J. Environ. Manag. 2011, 92, 407–418. [Google Scholar] [CrossRef] [PubMed]
  4. Tchounwou, P.B.; Yedjou, C.G.; Patlolla, A.K.; Sutton, D.J. Heavy Metal Toxicity and the Environment. In Molecular, Clinical and Environmental Toxicology; Springer: Berlin/Heidelberg, Germany, 2012; pp. 133–164. [Google Scholar]
  5. Wieczorek-Dąbrowska, M.; Tomza-Marciniak, A.; Pilarczyk, B.; Balicka-Ramisz, A. Roe and Red Deer as Bioindicators of Heavy Metals Contamination in North-Western Poland. Chem. Ecol. 2013, 29, 100–110. [Google Scholar] [CrossRef]
  6. Agwu, K.K.; Okoye, C.M.I.; Okeji, M.C.; Clifford, E.O. Potential Health Impacts of Heavy Metal Concentrations in Fresh and Marine Water Fishes Consumed in Southeast, Nigeria. Pak. J. Nutr. 2018, 17, 647–653. [Google Scholar] [CrossRef]
  7. Fu, Z.; Xi, S. The Effects of Heavy Metals on Human Metabolism. Toxicol. Mech. Methods 2020, 30, 167–176. [Google Scholar] [CrossRef] [PubMed]
  8. Manoj, S.R.; Karthik, C.; Kadirvelu, K.; Arulselvi, P.I.; Shanmugasundaram, T.; Bruno, B.; Rajkumar, M. Understanding the Molecular Mechanisms for the Enhanced Phytoremediation of Heavy Metals through Plant Growth Promoting Rhizobacteria: A Review. J. Environ. Manag. 2020, 254, 109779. [Google Scholar] [CrossRef]
  9. Jaishankar, M.; Tseten, T.; Anbalagan, N.; Mathew, B.B.; Beeregowda, K.N. Toxicity, Mechanism and Health Effects of Some Heavy Metals. Interdiscip. Toxicol. 2014, 7, 60–72. [Google Scholar] [CrossRef] [PubMed]
  10. Al-Samman, T. Effect of Heavy Metal Impurities in Secondary Mg Alloys on the Microstructure and Mechanical Properties during Deformation. Mater. Des. 2015, 65, 983–988. [Google Scholar] [CrossRef]
  11. Rehman, K.; Fatima, F.; Waheed, I.; Akash, M.S.H. Prevalence of Exposure of Heavy Metals and Their Impact on Health Consequences. J. Cell Biochem. 2018, 119, 157–184. [Google Scholar] [CrossRef]
  12. He, T.; Feng, X.; Guo, Y.; Qiu, G.; Li, Z.; Liang, L.; Lu, J. The Impact of Eutrophication on the Biogeochemical Cycling of Mercury Species in a Reservoir: A Case Study from Hongfeng Reservoir, Guizhou, China. Environ. Pollut. 2008, 154, 56–67. [Google Scholar] [CrossRef] [PubMed]
  13. Chowdhury, S.; Mazumder, M.A.J.; Al-Attas, O.; Husain, T. Heavy Metals in Drinking Water: Occurrences, Implications, and Future Needs in Developing Countries. Sci. Total Environ. 2016, 569–570, 476–488. [Google Scholar] [CrossRef]
  14. Asati, A.; Pichhode, M.; Nikhil, K. Effect of Heavy Metals on Plants: An Overview. Int. J. Appl. Innov. Eng. Manag. 2016, 5, 56–66. [Google Scholar]
  15. Ali, H.; Khan, E. What Are Heavy Metals? Long-Standing Controversy over the Scientific Use of the Term ‘Heavy Metals’—Proposal of a Comprehensive Definition. Toxicol. Environ. Chem. 2018, 100, 6–19. [Google Scholar] [CrossRef]
  16. Nurchi, V.; Crespo-Alonso, M.; Toso, L.; Lachowicz, J.; Crisponi, G. Chelation Therapy for Metal Intoxication: Comments from a Thermodynamic Viewpoint. Mini-Rev. Med. Chem. 2013, 13, 1541–1549. [Google Scholar] [CrossRef]
  17. Corrente, G.A.; Malacaria, L.; Beneduci, A.; Furia, E.; Marino, T.; Mazzone, G. Experimental and Theoretical Study on the Coordination Properties of Quercetin towards Aluminum(III), Iron(III) and Copper(II) in Aqueous Solution. J. Mol. Liq. 2021, 325, 115171. [Google Scholar] [CrossRef]
  18. Abd Elnabi, M.K.; Elkaliny, N.E.; Elyazied, M.M.; Azab, S.H.; Elkhalifa, S.A.; Elmasry, S.; Mouhamed, M.S.; Shalamesh, E.M.; Alhorieny, N.A.; Abd Elaty, A.E.; et al. Toxicity of Heavy Metals and Recent Advances in Their Removal: A Review. Toxics 2023, 11, 580. [Google Scholar] [CrossRef] [PubMed]
  19. Singh, V.; Singh, N.; Rai, S.N.; Kumar, A.; Singh, A.K.; Singh, M.P.; Sahoo, A.; Shekhar, S.; Vamanu, E.; Mishra, V. Heavy Metal Contamination in the Aquatic Ecosystem: Toxicity and Its Remediation Using Eco-Friendly Approaches. Toxics 2023, 11, 147. [Google Scholar] [CrossRef]
  20. Luo, A.; Wang, H.; Wang, Y.; Huang, Q.; Zhang, Q. A Novel Colorimetric and Turn-on Fluorescent Chemosensor for Iron (III) Ion Detection and Its Application to Cellular Imaging. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2016, 168, 37–44. [Google Scholar] [CrossRef]
  21. Barbieri, M. The Importance of Enrichment Factor (EF) and Geoaccumulation Index (Igeo) to Evaluate the Soil Contamination. J. Geol. Geophys. 2016, 5, 1–4. [Google Scholar] [CrossRef]
  22. Patrick, L. Mercury Toxicity and Antioxidants: Part I: Role of Glutathione and Alpha-Lipoic Acid in the Treatment of Mercury Toxicity (Mercury Toxicity). Altern. Med. Rev. 2002, 7, 456–472. [Google Scholar] [PubMed]
  23. Castagnetto, J.M.; Hennessy, S.W.; Roberts, V.A.; Getzoff, E.D.; Tainer, J.A.; Pique, M.E. MDB: The Metalloprotein Database and Browser at the Scripps Research Institute. Nucleic Acids Res. 2002, 30, 379–382. [Google Scholar] [CrossRef] [PubMed]
  24. Lane, T.W.; Morel, F.M.M. A Biological Function for Cadmium in Marine Diatoms. Proc. Natl. Acad. Sci. USA 2000, 97, 4627–4631. [Google Scholar] [CrossRef] [PubMed]
  25. Lambing, A.; Kachalsky, E.; Mueller, M.L. The Dangers of Iron Overload: Bring in the Iron Police. J. Am. Acad. Nurse Pract. 2012, 24, 175–183. [Google Scholar] [CrossRef] [PubMed]
  26. Kwiatkowski, J.L.; Cohen, A.R. Iron Chelation Therapy in Sickle-Cell Disease and Other Transfusion-Dependent Anemias. Hematol. Oncol. Clin. 2004, 18, 1355–1377. [Google Scholar] [CrossRef] [PubMed]
  27. Cappellini, M.D.; Taher, A. Deferasirox (Exjade®) for the Treatment of Iron Overload. Acta Haematol. 2009, 122, 165–173. [Google Scholar] [CrossRef]
  28. Abetz, L.; Baladi, J.-F.; Jones, P.; Rofail, D. The Impact of Iron Overload and Its Treatment on Quality of Life: Results from a Literature Review. Health Qual. Life Outcomes 2006, 4, 1–6. [Google Scholar] [CrossRef] [PubMed]
  29. World Health Organization Guidelines for Drinking-Water Quality, 4th ed.; Incorporating the 1st Addendum; WHO: Geneva, Switzerland, 2017.
  30. USEPA (U.S. Environmental Protection Agency). National Primary Drinking Water Regulations. Available online: https://www.epa.gov/ground-water-and-drinking-water/national-primary-drinking-water-regulations (accessed on 25 November 2023).
  31. Khan, F.U.; Rahman, A.U.; Jan, A.; Riaz, M. Toxic and Trace Metals (Pb, Cd, Zn, Cu, Mn, Ni, Co and Cr) in Dust, Dustfall/Soil. J. Chem. Soc. Pak. 2004, 26, 453–456. [Google Scholar]
  32. Ali, H.; Khan, E.; Sajad, M.A. Phytoremediation of Heavy Metals—Concepts and Applications. Chemosphere 2013, 91, 869–881. [Google Scholar] [CrossRef] [PubMed]
  33. Ali, H.; Khan, E. Environmental Chemistry in the Twenty-First Century. Environ. Chem. Lett. 2017, 15, 329–346. [Google Scholar] [CrossRef]
  34. Hashem, M.A.; Nur-A-Tomal, M.S.; Mondal, N.R.; Rahman, M.A. Hair Burning and Liming in Tanneries Is a Source of Pollution by Arsenic, Lead, Zinc, Manganese and Iron. Environ. Chem. Lett. 2017, 15, 501–506. [Google Scholar] [CrossRef]
  35. Ali, H.; Khan, E.; Ilahi, I. Environmental Chemistry and Ecotoxicology of Hazardous Heavy Metals: Environmental Persistence, Toxicity, and Bioaccumulation. J. Chem. 2019, 2019, 6730305. [Google Scholar] [CrossRef]
  36. Wang, G.; Yinglan, A.; Jiang, H.; Fu, Q.; Zheng, B. Modeling the Source Contribution of Heavy Metals in Surficial Sediment and Analysis of Their Historical Changes in the Vertical Sediments of a Drinking Water Reservoir. J. Hydrol. 2015, 520, 37–51. [Google Scholar] [CrossRef]
  37. Miedico, O.; Pompa, C.; Moscatelli, S.; Chiappinelli, A.; Carosielli, L.; Chiaravalle, A.E. Lead, Cadmium and Mercury in Canned and Unprocessed Tuna: Six-Years Monitoring Survey, Comparison with Previous Studies and Recommended Tolerable Limits. J. Food Compos. Anal. 2020, 94, 103638. [Google Scholar] [CrossRef]
  38. Abedin, M.J.; Khan, R.; Siddique, M.A.B.; Khan, A.H.A.N.; Islam, M.T.; Rashid, M.B. Metal (Loid) s in Tap-Water from Schools in Central Bangladesh (Mirpur): Source Apportionment, Water Quality, and Health Risks Appraisals. Heliyon 2023, 9, e15747. [Google Scholar] [CrossRef] [PubMed]
  39. Mulyaningsih, R.D.; Pratiwi, R.; Hasanah, A.N. An Update on the Use of Natural Pigments and Pigment Nanoparticle Adducts for Metal Detection Based on Colour Response. Biosensors 2023, 13, 554. [Google Scholar] [CrossRef]
  40. Zhou, H.; Chen, Y.; Yue, X.; Ren, D.; Liu, Y.; Yang, K. Identification and Hazard Analysis of Heavy Metal Sources in Agricultural Soils in Ancient Mining Areas: A Quantitative Method Based on the Receptor Model and Risk Assessment. J. Hazard. Mater. 2023, 445, 130528. [Google Scholar] [CrossRef]
  41. Singh, M.R. Impurities-Heavy Metals: IR Perspective. Int. J. Phys. Sci. 2007, 5, 1045–1058. [Google Scholar]
  42. Aragay, G.; Merkoci, A. Nanomaterials Application in Electrochemical Detection of Heavy Metals. Electrochim. Acta 2012, 84, 49–61. [Google Scholar] [CrossRef]
  43. Gong, T.; Liu, J.; Liu, X.; Liu, J.; Xiang, J.; Wu, Y. A Sensitive and Selective Sensing Platform Based on CdTe QDs in the Presence of L-Cysteine for Detection of Silver, Mercury and Copper Ions in Water and Various Drinks. Food Chem. 2016, 213, 306–312. [Google Scholar] [CrossRef]
  44. Chen, J.; Teo, K.C. Determination of Cadmium, Copper, Lead and Zinc in Water Samples by Flame Atomic Absorption Spectrometry after Cloud Point Extraction. Anal. Chim. Acta 2001, 450, 215–222. [Google Scholar] [CrossRef]
  45. Sitko, R.; Janik, P.; Zawisza, B.; Talik, E.; Margui, E.; Queralt, I. Green Approach for Ultratrace Determination of Divalent Metal Ions and Arsenic Species Using Total-Reflection X-ray Fluorescence Spectrometry and Mercapto-Modified Graphene Oxide Nanosheets as a Novel Adsorbent. Anal. Chem. 2015, 87, 3535–3542. [Google Scholar] [CrossRef] [PubMed]
  46. Pöykiö, R.; Perämäki, P. Acid Dissolution Methods for Heavy Metals Determination in Pine Needles. Environ. Chem. Lett. 2003, 1, 191–195. [Google Scholar] [CrossRef]
  47. Losev, V.N.; Buyko, O.V.; Trofimchuk, A.K.; Zuy, O.N. Silica Sequentially Modified with Polyhexamethylene Guanidine and Arsenazo I for Preconcentration and ICP–OES Determination of Metals in Natural Waters. Microchem. J. 2015, 123, 84–89. [Google Scholar] [CrossRef]
  48. Gałuszka, A.; Migaszewski, Z.M.; Namieśnik, J. Moving Your Laboratories to the Field—Advantages and Limitations of the Use of Field Portable Instruments in Environmental Sample Analysis. Environ. Res. 2015, 140, 593–603. [Google Scholar] [CrossRef] [PubMed]
  49. Sigurdson, G.T.; Tang, P.; Giusti, M.M. Natural Colorants: Food Colorants from Natural Sources. Annu. Rev. Food Sci. Technol. 2017, 8, 261–280. [Google Scholar] [CrossRef] [PubMed]
  50. Rodríguez-Mena, A.; Ochoa-Martínez, L.A.; González-Herrera, S.M.; Rutiaga-Quiñones, O.M.; González-Laredo, R.F.; Olmedilla-Alonso, B. Natural Pigments of Plant Origin: Classification, Extraction and Application in Foods. Food Chem. 2023, 398, 133908. [Google Scholar] [CrossRef] [PubMed]
  51. Giuliani, A.A.; Cichelli, A.; Lorenzo, C. Colors: Properties and Determination of Natural Pigments. In Encyclopedia of Food and Health; Elsevier: Amsterdam, The Netherlands, 2016; pp. 273–284. ISBN 0123849535. [Google Scholar]
  52. Oyama, K.; Yamada, T.; Ito, D.; Kondo, T.; Yoshida, K. Metal Complex Pigment Involved in the Blue Sepal Color Development of Hydrangea. J. Agric. Food Chem. 2015, 63, 7630–7635. [Google Scholar] [CrossRef]
  53. Fedenko, V.S.; Shemet, S.A.; Landi, M. UV–Vis Spectroscopy and Colorimetric Models for Detecting Anthocyanin-Metal Complexes in Plants: An Overview of in Vitro and in Vivo Techniques. J. Plant Physiol. 2017, 212, 13–28. [Google Scholar] [CrossRef]
  54. Christie, R.; Abel, A. Metal Complex Pigments. Phys. Sci. Rev. 2021, 6, 375–389. [Google Scholar] [CrossRef]
  55. Wybraniec, S.; Starzak, K.; Skopińska, A.; Szaleniec, M.; Słupski, J.; Mitka, K.; Kowalski, P.; Michałowski, T. Effects of Metal Cations on Betanin Stability in Aqueous-Organic Solutions. Food Sci. Biotechnol. 2013, 22, 353–363. [Google Scholar] [CrossRef]
  56. Saithongdee, A.; Praphairaksit, N.; Imyim, A. Electrospun Curcumin-Loaded Zein Membrane for Iron (III) Ions Sensing. Sens. Actuators B Chem. 2014, 202, 935–940. [Google Scholar] [CrossRef]
  57. Khaodee, W.; Aeungmaitrepirom, W.; Tuntulani, T. Effectively Simultaneous Naked-Eye Detection of Cu(II), Pb(II), Al(III) and Fe(III) Using Cyanidin Extracted from Red Cabbage as Chelating Agent. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2014, 126, 98–104. [Google Scholar] [CrossRef]
  58. Yılmaz, D.D.; Demirezen, D.A.; Mıhçıokur, H. Colorimetric Detection of Mercury Ion Using Chlorophyll Functionalized Green Silver Nanoparticles in Aqueous Medium. Surf. Interfaces 2021, 22, 100840. [Google Scholar] [CrossRef]
  59. Mohan, A.; Prakash, J. Fabrication of Eco-Friendly Hydrogel Strips for the Simultaneous Quantification of Heavy Metal Ions in Aqueous Environment. Dye. Pigment. 2022, 199, 110045. [Google Scholar] [CrossRef]
  60. Raj, S.; Shankaran, D.R. Curcumin Based Biocompatible Nanofibers for Lead Ion Detection. Sens. Actuators B Chem. 2016, 226, 318–325. [Google Scholar] [CrossRef]
  61. Sheikhzadeh, E.; Naji-Tabasi, S.; Verdian, A.; Kolahi-Ahari, S. Equipment-Free and Visual Detection of Pb2+ Ion Based on Curcumin-Modified Bacterial Cellulose Nanofiber. J. Iran. Chem. Soc. 2022, 19, 283–290. [Google Scholar] [CrossRef]
  62. Kumar, P.; Paul, W.; Sharma, C.P. Curcumin Stabilized Gold Nanoparticle-Based Colorimetric Sensing of Mercury (II). Trends Biomater. Artif. Organs 2014, 28, 121–123. [Google Scholar]
  63. Zin, M.M.; Márki, E.; Bánvölgyi, S. Conventional Extraction of Betalain Compounds from Beetroot Peels with Aqueous Ethanol Solvent. Acta Aliment. 2020, 49, 163–169. [Google Scholar] [CrossRef]
  64. Malacaria, L.; La Torre, C.; Furia, E.; Fazio, A.; Caroleo, M.C.; Cione, E.; Gallelli, L.; Marino, T.; Plastina, P. Aluminum (III), Iron (III) and Copper (II) Complexes of Luteolin: Stability, Antioxidant, and Anti-Inflammatory Properties. J. Mol. Liq. 2022, 345, 117895. [Google Scholar] [CrossRef]
  65. Eno, E.A.; Shagal, M.H.; Godfrey, O.C.; Ngana, O.C.; Ekong, J.E.; Benjamin, I.; Ugi, B.U.; Louis, H. Computational Study of the Interaction of Metal Ions (Na+, K+, Mg2+, Ca2+, and Al3+) with Quercetin and Its Antioxidant Properties. J. Indian. Chem. Soc. 2023, 100, 101059. [Google Scholar] [CrossRef]
  66. Swart, M.; Gruden, M. Spinning around in Transition-Metal Chemistry. Acc. Chem. Res. 2016, 49, 2690–2697. [Google Scholar] [CrossRef] [PubMed]
  67. Alzahrani, E. Colorimetric Detection Based on Localized Surface Plasmon Resonance Optical Characteristics for Sensing of Mercury Using Green-Synthesized Silver Nanoparticles. J. Anal. Methods Chem. 2020, 2020, 6026312. [Google Scholar] [CrossRef] [PubMed]
  68. Mahimwalla, Z.; Yager, K.G.; Mamiya, J.I.; Shishido, A.; Priimagi, A.; Barrett, C.J. Azobenzene Photomechanics: Prospects and Potential Applications. Polym. Bull. 2012, 69, 967–1006. [Google Scholar] [CrossRef]
  69. Sousa, S.F.; Ribeiro, A.J.M.; Neves, R.P.P.; Brás, N.F.; Cerqueira, N.M.F.S.A.; Fernandes, P.A.; Ramos, M.J. Application of Quantum Mechanics/Molecular Mechanics Methods in the Study of Enzymatic Reaction Mechanisms. WIREs Comput. Mol. Sci. 2017, 7, 1281. [Google Scholar] [CrossRef]
  70. Starikov, A.G.; Starikova, A.A.; Minkin, V.I. Quantum-Chemical Study of (Z)-6,8-Di-Tert-Butyl-N-(4-Methoxyphenyl)-3-((4-Methoxyphenyl)Imino)-3H-Phenoxazine-2-Amine Complexation with Cobalt Bis(Chelate)s. Doklady Chem. 2017, 476, 215–218. [Google Scholar] [CrossRef]
  71. Ajaykamal, T.; Palaniandavar, M. Mononuclear Nickel(Ii)-Flavonolate Complexes of Tetradentate Tripodal 4N Ligands as Structural and Functional Models for Quercetin 2,4-Dioxygenase: Structures, Spectra, Redox and Dioxygenase Activity. RSC Adv. 2023, 13, 24674–24690. [Google Scholar] [CrossRef] [PubMed]
  72. Sun, S.; Chen, W.; Cao, W.; Zhang, F.; Song, J.; Tian, C. Research on the Chelation between Quercetin and Cr (III) Ion by Density Functional Theory (DFT) Method. J. Mol. Struct. THEOCHEM 2008, 860, 40–44. [Google Scholar] [CrossRef]
  73. Furia, E.; Marino, T.; Russo, N. Insights into the Coordination Mode of Quercetin with the Al(III) Ion from a Combined Experimental and Theoretical Study. Dalton Trans. 2014, 43, 7269–7274. [Google Scholar] [CrossRef]
  74. Cornard, J.P.; Merlin, J.C. Spectroscopic and Structural Study of Complexes of Quercetin with Al(III). J. Inorg. Biochem. 2002, 92, 19–27. [Google Scholar] [CrossRef]
  75. Ren, J.; Meng, S.; Lekka, C.E.; Kaxiras, E. Complexation of Flavonoids with Iron: Structure and Optical Signatures. J. Phys. Chem. B 2008, 112, 1845–1850. [Google Scholar] [CrossRef] [PubMed]
  76. Leopoldini, M.; Russo, N.; Chiodo, S.; Toscano, M. Iron Chelation by the Powerful Antioxidant Flavonoid Quercetin. J. Agric. Food Chem. 2006, 54, 6343–6351. [Google Scholar] [CrossRef] [PubMed]
  77. Krishnan, R.; Yang, K.; Hashem, K.; Jiang, J. Metallated Porphyrinic Metal− Organic Frameworks for CO2 Conversion to HCOOH: A Computational Screening and Mechanistic Study. Mol. Catal. 2022, 527, 112407. [Google Scholar] [CrossRef]
  78. Mollaamin, F.; Mohammadian, N.T.; Najaflou, N.; Monajjemi, M. Iranian Qara Qat Fruit (Redcurrant) in Arasbaran Forests as the Resource of Anthocyanin Pigments in Formation of [ACN-Mg2+/Al3+/Ga3+/Sn2+/Cr3+/Fe3+] Chelation Clusters. SN Appl. Sci. 2021, 3, 404. [Google Scholar] [CrossRef]
  79. Hong, B.; Zhao, D.; Li, W.; Li, S. Generalized Energy-Based Fragmentation Approach for Structures and Properties of Periodic Condensed Phase Systems. In Comprehensive Computational Chemistry; Elsevier: Amsterdam, The Netherlands, 2024; pp. 129–138. [Google Scholar]
  80. Spiegel, M.; Sroka, Z. Quantum-Mechanical Characteristics of Apigenin: Antiradical, Metal Chelation and Inhibitory Properties in Physiologically Relevant Media. Fitoterapia 2023, 164, 105352. [Google Scholar] [CrossRef]
  81. Renger, T. Theory of Excitation Energy Transfer: From Structure to Function. Photosynth. Res. 2009, 102, 471–485. [Google Scholar] [CrossRef] [PubMed]
  82. Estévez, L.; Sánchez-Lozano, M.; Mosquera, R.A. Complexation of Common Metal Cations by Cyanins: Binding Affinity and Molecular Structure. Int. J. Quantum Chem. 2019, 119, 25834. [Google Scholar] [CrossRef]
  83. Bechaieb, R.; Fredj, A.B.; Akacha, A.B.; Gérard, H. Interactions of Copper (II) and Zinc (II) with Chlorophyll: Insights from Density Functional Theory Studies. New J. Chem. 2016, 40, 4543–4549. [Google Scholar] [CrossRef]
  84. Frison, G.; Ohanessian, G. A Comparative Study of Semiempirical, Ab Initio, and DFT Methods in Evaluating Metal-Ligand Bond Strength, Proton Affinity, and Interactions between First and Second Shell Ligands in Zn-Biomimetic Complexes. J. Comput. Chem. 2008, 29, 416–433. [Google Scholar] [CrossRef]
  85. Bursch, M.; Mewes, J.; Hansen, A.; Grimme, S. Best-Practice DFT Protocols for Basic Molecular Computational Chemistry. Angew. Chem. Int. Ed. 2022, 61, 5735. [Google Scholar] [CrossRef]
  86. Enaru, B.; Drețcanu, G.; Pop, T.D.; Stǎnilǎ, A.; Diaconeasa, Z. Anthocyanins: Factors Affecting Their Stability and Degradation. Antioxidants 2021, 10, 1967. [Google Scholar] [CrossRef]
  87. Ghiasi, M.; Molaei, M.; Zahedi, M. QM Study of Complexation between Natural Bilirubin and Poly-Terthiophene Carboxylic Acid-Mn(II) as a Biosensor: Temperature and Interferences Effect. J. Theor. Comput. Chem. 2015, 14, 1550048. [Google Scholar] [CrossRef]
  88. Fazary, A.E.; Alshihri, A.S.; Alfaifi, M.Y.; Saleh, K.A.; Elbehairi, S.E.I.; Fawy, K.F.; Abd-Rabboh, H.S.M. Gibbs Energies of Protonation and Complexation of Platinum and Vanadate Metal Ions with Naringenin and Phenolic Acids: Theoretical Calculations Associated with Experimental Values. J. Chem. Thermodyn. 2016, 100, 7–21. [Google Scholar] [CrossRef]
  89. Decherchi, S.; Masetti, M.; Vyalov, I.; Rocchia, W. Implicit Solvent Methods for Free Energy Estimation. Eur. J. Med. Chem. 2015, 91, 27–42. [Google Scholar] [CrossRef]
  90. Orzeł, Ł.; Waś, J.; Kania, A.; Susz, A.; Rutkowska-Zbik, D.; Staroń, J.; Witko, M.; Stochel, G.; Fiedor, L. Factors Controlling the Reactivity of Divalent Metal Ions towards Pheophytin a. J. Biol. Inorg. Chem. 2017, 22, 941–952. [Google Scholar] [CrossRef]
  91. Eilmes, A. Solvatochromic Probe in Molecular Solvents: Implicit versus Explicit Solvent Model. Theor. Chem. Acc. 2014, 133, 1–13. [Google Scholar] [CrossRef]
  92. Ehlert, S.; Stahn, M.; Spicher, S.; Grimme, S. Robust and Efficient Implicit Solvation Model for Fast Semiempirical Methods. J. Chem. Theory Comput. 2021, 17, 4250–4261. [Google Scholar] [CrossRef]
  93. Cossi, M.; Rega, N.; Scalmani, G.; Barone, V. Energies, Structures, and Electronic Properties of Molecules in Solution with the C-PCM Solvation Model. J. Comput. Chem. 2003, 24, 669–681. [Google Scholar] [CrossRef]
  94. Okur, A.; Wickstrom, L.; Layten, M.; Geney, R.; Song, K.; Hornak, V.; Simmerling, C. Improved Efficiency of Replica Exchange Simulations through Use of a Hybrid Explicit/Implicit Solvation Model. J. Chem. Theory Comput. 2006, 2, 420–433. [Google Scholar] [CrossRef]
  95. Brancato, G.; Rega, N.; Barone, V. A Hybrid Explicit/Implicit Solvation Method for First-Principle Molecular Dynamics Simulations. J. Chem. Phys. 2008, 128, 144501. [Google Scholar] [CrossRef]
  96. Lee, M.S.; Salsbury, F.R., Jr.; Olson, M.A. An Efficient Hybrid Explicit/Implicit Solvent Method for Biomolecular Simulations. J. Comput. Chem. 2004, 25, 1967–1978. [Google Scholar] [CrossRef]
  97. Wang, Z.; Zhao, J.; Cai, Q. CO 2 Electroreduction Performance of a Single Transition Metal Atom Supported on Porphyrin-like Graphene: A Computational Study. Phys. Chem. Chem. Phys. 2017, 19, 23113–23121. [Google Scholar] [CrossRef]
  98. Pranowo, H.D. Kimia Komputasi; Universitas Gadjah Mada: Yogyakarta, Indonesia, 2000. [Google Scholar]
  99. Meyer, M. Ab Initio Study of Flavonoids. Int. J. Quantum Chem. 2000, 76, 724–732. [Google Scholar] [CrossRef]
  100. Yan, L.K.; Pomogaeva, A.; Gu, F.L.; Aoki, Y. Theoretical Study on Nonlinear Optical Properties of Metalloporphyrin Using Elongation Method. Theor. Chem. Acc. 2010, 125, 511–520. [Google Scholar] [CrossRef]
  101. Linnanto, J.; Korppi-Tommola, J. Structural and Spectroscopic Properties of Mg-Bacteriochlorin and Methyl Bacteriochlorophyllides a, b, g, and h Studied by Semiempirical, Ab Initio, and Density Functional Molecular Orbital Methods. J. Phys. Chem. A 2004, 108, 5872–5882. [Google Scholar] [CrossRef]
  102. Freitas, A.A.; Shimizu, K.; Dias, L.G.; Quina, F.H. A Computational Study of Substituted Flavylium Salts and Their Quinonoidal Conjugate-Bases: S0-> S1 Electronic Transition, Absolute PKa and Reduction Potential Calculations by DFT and Semiempirical Methods. J. Braz. Chem. Soc. 2007, 18, 1537–1546. [Google Scholar] [CrossRef]
  103. Cao, Y.; Chen, Z.-J.; Jiang, H.-D.; Chen, J.-Z. Computational Studies of the Regioselectivities of COMT-Catalyzed Meta-/Para.-O Methylations of Luteolin and Quercetin. J. Phys. Chem. B 2014, 118, 470–481. [Google Scholar] [CrossRef]
  104. Cavalheiro, E.P.; Cardoso, M.V.P.; De Souza, S.P.; Sousa, M.E.S.; Allan Leite Dos Reis, M. Correlation Between Clar Sextet Theory and Electronic Properties of Chromophores: An Ab Initio Study. J. Comput. Theor. Nanosci. 2017, 14, 1524–1527. [Google Scholar] [CrossRef]
  105. Yang, C.-H.; Huang, K.-S.; Wang, Y.-T.; Shaw, J.-F. A Review of Bacteriochlorophyllides: Chemical Structures and Applications. Molecules 2021, 26, 1293. [Google Scholar] [CrossRef]
  106. Tronrud, D.E.; Schmid, M.F.; Matthews, B.W. Structure and X-ray Amino Acid Sequence of a Bacteriochlorophyll a Protein from Prosthecochloris Aestuarii Refined at 1.9 Å Resolution. J. Mol. Biol. 1986, 188, 443–454. [Google Scholar] [CrossRef]
  107. Freer, A.; Prince, S.; Sauer, K.; Papiz, M.; Lawless, A.H.; McDermott, G.; Cogdell, R.; Isaacs, N.W. Pigment–Pigment Interactions and Energy Transfer in the Antenna Complex of the Photosynthetic Bacterium Rhodopseudomonas acidophila. Structure 1996, 4, 449–462. [Google Scholar] [CrossRef] [PubMed]
  108. Hofmann, E.; Wrench, P.M.; Sharples, F.P.; Hiller, R.G.; Welte, W.; Diederichs, K. Structural Basis of Light Harvesting by Carotenoids: Peridinin-Chlorophyll-Protein from Amphidinium carterae. Science 1996, 272, 1788–1791. [Google Scholar] [CrossRef] [PubMed]
  109. Lanzo, I.; Russo, N.; Sicilia, E. First-Principle Time-Dependent Study of Magnesium-Containing Porphyrin-like Compounds Potentially Useful for Their Application in Photodynamic Therapy. J. Phys. Chem. B 2008, 112, 4123–4130. [Google Scholar] [CrossRef] [PubMed]
  110. Tian, C.; Liu, X.; Chang, Y.; Wang, R.; Lv, T.; Cui, C.; Liu, M. Investigation of the Anti-Inflammatory and Antioxidant Activities of Luteolin, Kaempferol, Apigenin and Quercetin. S. Afr. J. Bot. 2021, 137, 257–264. [Google Scholar] [CrossRef]
  111. Tanaka, M.; Hayashi, T.; Morita, S. The Roles of Water Molecules at the Biointerface of Medical Polymers. Polym. J. 2013, 45, 701–710. [Google Scholar] [CrossRef]
  112. Guo, X.; Zhang, X.; Zhao, S.; Huang, Q.; Xue, J. High Adsorption Capacity of Heavy Metals on Two-Dimensional MXenes: An Ab Initio Study with Molecular Dynamics Simulation. Phys. Chem. Chem. Phys. 2016, 18, 228–233. [Google Scholar] [CrossRef] [PubMed]
  113. Narin, P.; Abbas, J.M.A.; Atmaca, G.; Kutlu, E.; Lisesivdin, S.B.; Ozbay, E. Ab Initio Study of Electronic Properties of Armchair Graphene Nanoribbons Passivated with Heavy Metal Elements. Solid. State Commun. 2019, 296, 8–11. [Google Scholar] [CrossRef]
  114. Qi, C.; Xu, X.; Chen, Q.; Liu, H.; Min, X.; Fourie, A.; Chai, L. Ab Initio Calculation of the Adsorption of As, Cd, Cr, and Hg Heavy Metal Atoms onto the Illite (001) Surface: Implications for Soil Pollution and Reclamation. Environ. Pollut. 2022, 312, 120072. [Google Scholar] [CrossRef] [PubMed]
  115. Sousa, S.F.; Fernandes, P.A.; Ramos, M.J. General Performance of Density Functionals. J. Phys. Chem. A 2007, 111, 10439–10452. [Google Scholar] [CrossRef]
  116. Liu, S.-B. Conceptual Density Functional Theory and Some Recent Developments. Acta Phys. Chim. Sin. 2009, 25, 590–600. [Google Scholar]
  117. Guo, P.; Chen, W.; Song, J.; Cao, W.; Tian, C. A DFT Study of the Interaction between Butein Anion and Metal Cations (M= Mg2+, Cr2+, Fe2+, and Cu2+): Taking an Insight into Its Chelating Property. J. Mol. Struct. THEOCHEM 2008, 849, 33–36. [Google Scholar] [CrossRef]
  118. Muñoz, V.A.; Ferrari, G.V.; Sancho, M.I.; Montaña, M.P. Spectroscopic and Thermodynamic Study of Chrysin and Quercetin Complexes with Cu(II). J. Chem. Eng. Data 2016, 61, 987–995. [Google Scholar] [CrossRef]
  119. Amat, A.; Clementi, C.; Miliani, C.; Romani, A.; Sgamellotti, A.; Fantacci, S. Complexation of Apigenin and Luteolin in Weld Lake: A DFT/TDDFT Investigation. Phys. Chem. Chem. Phys. 2010, 12, 6672–6684. [Google Scholar] [CrossRef] [PubMed]
  120. Asadi, N.; Ramezanzadeh, M.; Bahlakeh, G.; Ramezanzadeh, B. Theoretical MD/DFT Computer Explorations and Surface-Electrochemical Investigations of the Zinc/Iron Metal Cations Interactions with Highly Active Molecules from Lemon Balm Extract toward the Steel Corrosion Retardation in Saline Solution. J. Mol. Liq. 2020, 310, 113220. [Google Scholar] [CrossRef]
  121. Ekennia, A.C.; Osowole, A.A.; Onwudiwe, D.C.; Babahan, I.; Ibeji, C.U.; Okafor, S.N.; Ujam, O.T. Synthesis, Characterization, Molecular Docking, Biological Activity and Density Functional Theory Studies of Novel 1,4-naphthoquinone Derivatives and Pd(II), Ni(II) and Co(II) Complexes. Appl. Organomet. Chem. 2018, 32, e4310. [Google Scholar] [CrossRef]
  122. Fiorucci, S.; Golebiowski, J.; Cabrol-Bass, D.; Antonczak, S. DFT Study of Quercetin Activated Forms Involved in Antiradical, Antioxidant, and Prooxidant Biological Processes. J. Agric. Food Chem. 2007, 55, 903–911. [Google Scholar] [CrossRef] [PubMed]
  123. Tanui, H.K.; Nkabyo, H.A.; Pearce, B.H.; Hussein, A.A.; Lopis, A.S.; Luckay, R.C. Iron (III) and Copper (II) Complexes Derived from the Flavonoids Morin and Quercetin: Chelation, Crystal Structure and DFT Studies. J. Mol. Struct. 2022, 1257, 132591. [Google Scholar] [CrossRef]
  124. Janeeshma, E.; Rajan, V.K.; Puthur, J.T. Spectral Variations Associated with Anthocyanin Accumulation; an Apt Tool to Evaluate Zinc Stress in Zea mays L. Chem. Ecol. 2021, 37, 32–49. [Google Scholar] [CrossRef]
  125. Beneduci, A.; Corrente, G.A.; Marino, T.; Aiello, D.; Bartella, L.; Di Donna, L.; Napoli, A.; Russo, N.; Romeo, I.; Furia, E. Insight on the Chelation of Aluminum (III) and Iron (III) by Curcumin in Aqueous Solution. J. Mol. Liq. 2019, 296, 111805. [Google Scholar] [CrossRef]
  126. Argaman, N.; Makov, G. Density Functional Theory: An Introduction. Am. J. Phys. 2000, 68, 69–79. [Google Scholar] [CrossRef]
  127. Zhang, G.; Musgrave, C.B. Comparison of DFT Methods for Molecular Orbital Eigenvalue Calculations. J. Phys. Chem. A 2007, 111, 1554–1561. [Google Scholar] [CrossRef] [PubMed]
  128. Hernández-Marin, E.; Barbosa, A.; Martínez, A. The Metal Cation Chelating Capacity of Astaxanthin. Does This Have Any Influence on Antiradical Activity? Molecules 2012, 17, 1039–1054. [Google Scholar] [CrossRef] [PubMed]
  129. Yang, G. Density Functional Calculations: Recent Progresses of Theory and Application; BoD–Books on Demand; Intech Open: Rijeka, Croatia, 2018; ISBN 1789231329. [Google Scholar]
  130. Sengupta, P.; Pal, U.; Roy, P.; Samanta, T.; Chattopadhyay, N.; Sen, K.; Bose, A. Effect of a Metal Ion in Modulating the Binding Interaction of a Dietary Flavonoid with Bovine Serum Albumin and DNA: A Spectroscopic and Theoretical Approach. ACS Food Sci. Technol. 2022, 2, 114–124. [Google Scholar] [CrossRef]
  131. Attarki, J.; Khnifira, M.; Boumya, W.; Mahsoune, A.; Lemdek, E.M.; Sadiq, M.; Abdennouri, M.; Kaya, S.; Barka, N. Adsorption Mechanism Investigation of Methylthioninium Chloride Dye onto Some Metal Phosphates Using Monte Carlo Dynamic Simulations and DFT Calculations. Inorg. Chem. Commun. 2023, 149, 110436. [Google Scholar] [CrossRef]
  132. Tirado-Rives, J.; Jorgensen, W.L. Performance of B3LYP Density Functional Methods for a Large Set of Organic Molecules. J. Chem. Theory Comput. 2008, 4, 297–306. [Google Scholar] [CrossRef] [PubMed]
  133. Šimunková, M.; Valko, M.; Bučinský, L.; Malček, M. Structure Functionality Relationship of Flavonoids (Myricetin, Morin, Taxifolin and 3′, 4′-Dihydroxyflavone). A Computational Study via the Cupric Ion Probe. J. Mol. Struct. 2020, 1222, 128923. [Google Scholar] [CrossRef]
  134. Semwal, R.B.; Semwal, D.K.; Combrinck, S.; Viljoen, A. Butein: From Ancient Traditional Remedy to Modern Nutraceutical. Phytochem. Lett. 2015, 11, 188–201. [Google Scholar] [CrossRef]
  135. Cornard, J.P.; Dangleterre, L.; Lapouge, C. Computational and Spectroscopic Characterization of the Molecular and Electronic Structure of the Pb(II)−Quercetin Complex. J. Phys. Chem. A 2005, 109, 10044–10051. [Google Scholar] [CrossRef] [PubMed]
  136. Bodini, M.E.; Copia, G.; Tapia, R.; Leighton, F.; Herrera, L. Iron Complexes of Quercetin in Aprotic Medium. Redox Chemistry and Interaction with Superoxide Anion Radical. Polyhedron 1999, 18, 2233–2239. [Google Scholar] [CrossRef]
  137. Wu, T.-W.; Zeng, L.-H.; Wu, J.; Fung, K.-P. Morin: A Wood Pigment That Protects Three Types of Human Cells in the Cardiovascular System against Oxyradical Damage. Biochem. Pharmacol. 1994, 47, 1099–1103. [Google Scholar] [CrossRef]
  138. Mollaamin, F.; Monajjemi, M. Thermodynamic and IR Spectral Study of Metal Cations–Anthocyanin Chelation: Mechanism of Formation of Pigments. Russ. J. Phys. Chem. A 2020, 94, 1887–1901. [Google Scholar] [CrossRef]
  139. Pratiwi, R.; Ibrahim, S.; Tjahjono, D.H. Reactivity and Stability of Metalloporphyrin Complex Formation: DFT and Experimental Study. Molecules 2020, 25, 4221. [Google Scholar] [CrossRef] [PubMed]
  140. Martínez, A.; Vargas, R. Electron Donor–Acceptor Properties of Metal Atoms Interacting with Pterins. New J. Chem. 2010, 34, 2988. [Google Scholar] [CrossRef]
  141. Linnanto, J.; Korppi-Tommola, J. Semiempirical PM5 Molecular Orbital Study on Chlorophylls and Bacteriochlorophylls: Comparison of Semiempirical, Ab Initio, and Density Functional Results. J. Comput. Chem. 2004, 25, 123–138. [Google Scholar] [CrossRef] [PubMed]
  142. Dangles, O.; Elhabiri, M.; Brouillard, R. Kinetic and Thermodynamic Investigation of the Aluminium–Anthocyanin Complexation in Aqueous Solution. J. Chem. Soc. Perkin Trans. 2 1994, 2587–2596. [Google Scholar] [CrossRef]
  143. Elhabiri, M.; Figueiredo, P.; Toki, K.; Saito, N.; Brouillard, R. Anthocyanin–Aluminium and –Gallium Complexes in Aqueous Solution. J. Chem. Soc. Perkin Trans. 2 1997, 355–362. [Google Scholar] [CrossRef]
  144. Guaratini, T.; Vessecchi, R.; Pinto, E.; Colepicolo, P.; Lopes, N.P. Balance of Xanthophylls Molecular and Protonated Molecular Ions in Electrospray Ionization. J. Mass. Spectrom. 2005, 40, 963–968. [Google Scholar] [CrossRef] [PubMed]
  145. Chen, M.; Cai, Z.-L. Theoretical Study on the Thermodynamic Properties of Chlorophyll D-Peptides Coordinating Ligand. Biochim. Biophys. Acta (BBA)—Bioenerg. 2007, 1767, 603–609. [Google Scholar] [CrossRef]
  146. Chassaing, S.; Lefeuvre, D.; Jacquet, R.; Jourdes, M.; Ducasse, L.; Galland, S.; Grelard, A.; Saucier, C.; Teissedre, P.; Dangles, O.; et al. Physicochemical Studies of New Anthocyano-Ellagitannin Hybrid Pigments: About the Origin of the Influence of Oak C -Glycosidic Ellagitannins on Wine Color. Eur. J. Org. Chem. 2010, 2010, 55–63. [Google Scholar] [CrossRef]
  147. Phan, K.; Van Den Broeck, E.; Van Speybroeck, V.; De Clerck, K.; Raes, K.; De Meester, S. The Potential of Anthocyanins from Blueberries as a Natural Dye for Cotton: A Combined Experimental and Theoretical Study. Dye. Pigment. 2020, 176, 108180. [Google Scholar] [CrossRef]
  148. Cornard, J.P.; Merlin, J.C.; Boudet, A.C.; Vrielynck, L. Structural Study of Quercetin by Vibrational and Electronic Spectroscopies Combined with Semiempirical Calculations. Biospectroscopy 1997, 3, 183–193. [Google Scholar] [CrossRef]
  149. Daoud, I.; Mesmoudi, M.; Ghalem, S. MM/QM Study: Interactions of Copper(II) and Mercury(II) with Food Dyes in Aqueous Solutions. Int. J. Chem. Anal. Sci. 2013, 4, 49–56. [Google Scholar] [CrossRef]
  150. Visa, A.; Maranescu, B.; Lupa, L.; Crisan, L.; Borota, A. New Efficient Adsorbent Materials for the Removal of Cd(II) from Aqueous Solutions. Nanomaterials 2020, 10, 899. [Google Scholar] [CrossRef] [PubMed]
  151. Karplus, M.; McCammon, J.A. Molecular Dynamics Simulations of Biomolecules. Nat. Struct. Biol. 2002, 9, 646–652. [Google Scholar] [CrossRef]
  152. Hollingsworth, S.A.; Dror, R.O. Molecular Dynamics Simulation for All. Neuron 2018, 99, 1129–1143. [Google Scholar] [CrossRef] [PubMed]
  153. March, R.; Brodbelt, J. Analysis of Flavonoids: Tandem Mass Spectrometry, Computational Methods, and NMR. J. Mass. Spectrom. 2008, 43, 1581–1617. [Google Scholar] [CrossRef] [PubMed]
  154. Teixeira, N.; Cruz, L.; Brás, N.F.; Mateus, N.; Ramos, M.J.; de Freitas, V. Structural Features of Copigmentation of Oenin with Different Polyphenol Copigments. J. Agric. Food Chem. 2013, 61, 6942–6948. [Google Scholar] [CrossRef] [PubMed]
  155. Moin, S.T.; Hofer, T.S. Hydration of Porphyrin and Mg-Porphyrin: Ab Initio Quantum Mechanical Charge Field Molecular Dynamics Simulations. Mol. Biosyst. 2014, 10, 117–127. [Google Scholar] [CrossRef]
  156. Singh, A.K.; Usman, M.; Sciortino, G.; Garribba, E.; Rath, S.P. Through-Space Spin Coupling in a Silver(II) Porphyrin Dimer upon Stepwise Oxidations: AgII···AgII, AgII···AgIII, and AgIII···AgIII Metallophilic Interactions. Chem. Eur. J. 2019, 25, 10098–10110. [Google Scholar] [CrossRef]
  157. Huamán Aguirre, A.A.; Salazar Salinas, K.; Quintana Cáceda, M. Molecular Interaction of Natural Dye Based on Zea mays and Bixa Orellana to Nanocrystalline TiO2 into Dye Sensitized Solar Cells. J. Electrochem. Sci. Eng. 2021, 11, 179–195. [Google Scholar] [CrossRef]
  158. Moradi, A.; Ranjbar, Z.; Guo, L.; Javadpour, S.; Ramezanzadeh, B. Molecular Dynamic (MD) Simulation and Electrochemical Assessments of the Satureja Hortensis Extract for the Construction of Effective Zinc-Based Protective Film on Carbon Steel. J. Mol. Liq. 2021, 338, 116606. [Google Scholar] [CrossRef]
  159. Hassanzadeh, M.K.; Tayarani Najaran, Z.; Nasery, M.; Emami, S.A. Summer Savory (Satureja hortensis L.) Oils. In Essential Oils in Food Preservation, Flavor and Safety; Elsevier: Amsterdam, The Netherlands, 2016; pp. 757–764. [Google Scholar]
  160. Dehghani, A.; Bahlakeh, G.; Ramezanzadeh, B.; Mostafatabar, A.H. Construction of a Zinc-Centered Metal–Organic Film with High Anti-Corrosion Potency through Covalent-Bonding between the Natural Flavonoid-Based Molecules (Quercetin)/Divalent-Zinc: Computer Modeling (Integrated-DFT&MC/MD)/Electrochemical-Surface Assessments. J. Ind. Eng. Chem. 2020, 88, 382–395. [Google Scholar]
  161. Liu, H.; Fan, B.; Fan, G.; Zhao, X.; Liu, Z.; Hao, H.; Yang, B. Long-Term Protective Mechanism of Poly (N-Methylaniline)/Phosphate One-Step Electropolymerized Coatings for Copper in 3.5% NaCl Solution. J. Alloys Compd. 2021, 872, 159752. [Google Scholar] [CrossRef]
  162. D’Ambrosio, M.; Santos, A.C.; Alejo-Armijo, A.; Parola, A.J.; Costa, P.M. Light-Mediated Toxicity of Porphyrin-like Pigments from a Marine Polychaeta. Mar. Drugs 2020, 18, 302. [Google Scholar] [CrossRef]
  163. Khoo, H.E.; Azlan, A.; Tang, S.T.; Lim, S.M. Anthocyanidins and Anthocyanins: Colored Pigments as Food, Pharmaceutical Ingredients, and the Potential Health Benefits. Food Nutr. Res. 2017, 61, 1361779. [Google Scholar] [CrossRef] [PubMed]
  164. Ahliha, A.H.; Nurosyid, F.; Supriyanto, A.; Kusumaningsih, T. Optical Properties of Anthocyanin Dyes on TiO 2 as Photosensitizers for Application of Dye-Sensitized Solar Cell (DSSC). IOP Conf. Ser. Mater. Sci. Eng. 2018, 333, 012018. [Google Scholar] [CrossRef]
  165. Hsissou, R.; Azogagh, M.; Benhiba, F.; Echihi, S.; Galai, M.; Shaim, A.; Bahaj, H.; Briche, S.; Kaya, S.; Serdaroğlu, G.; et al. Insight of Development of Two Cured Epoxy Polymer Composite Coatings as Highly Protective Efficiency for Carbon Steel in Sodium Chloride Solution: DFT, RDF, FFV and MD Approaches. J. Mol. Liq. 2022, 360, 119406. [Google Scholar] [CrossRef]
  166. Hsissou, R.; Benhiba, F.; El Aboubi, M.; Abbout, S.; Benzekri, Z.; Safi, Z.; Rafik, M.; Bahaj, H.; Kaba, M.; Galai, M.; et al. Synthesis and Performance of Two Ecofriendly Epoxy Resins as a Highly Efficient Corrosion Inhibition for Carbon Steel in 1 M HCl Solution: DFT, RDF, FFV and MD Approaches. Chem. Phys. Lett. 2022, 806, 139995. [Google Scholar] [CrossRef]
  167. Burg, J.S.; Ingram, J.R.; Venkatakrishnan, A.J.; Jude, K.M.; Dukkipati, A.; Feinberg, E.N.; Angelini, A.; Waghray, D.; Dror, R.O.; Ploegh, H.L.; et al. Structural Basis for Chemokine Recognition and Activation of a Viral G Protein–Coupled Receptor. Science 2015, 347, 1113–1117. [Google Scholar] [CrossRef]
  168. Dash, B.; Dash, B.; Rath, S.S. A Thorough Understanding of the Adsorption of Ni (II), Cd (II) and Zn (II) on Goethite Using Experiments and Molecular Dynamics Simulation. Sep. Purif. Technol. 2020, 240, 116649. [Google Scholar] [CrossRef]
  169. Giri, A.K.; Cordeiro, M.N.D.S. Heavy Metal Ion Separation from Industrial Wastewater Using Stacked Graphene Membranes: A Molecular Dynamics Simulation Study. J. Mol. Liq. 2021, 338, 116688. [Google Scholar] [CrossRef]
  170. Liu, L.; Zhang, C.; Jiang, W.; Li, X.; Dai, Y.; Jia, H. Understanding the Sorption Behaviors of Heavy Metal Ions in the Interlayer and Nanopore of Montmorillonite: A Molecular Dynamics Study. J. Hazard. Mater. 2021, 416, 125976. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Chemical structure of natural pigments: (a) betanin, a derivative of the betalain pigment; (b) cyanidin; (c) luteolin; (d) quercetin; (e) chlorophyll; (f) curcumin.
Figure 1. Chemical structure of natural pigments: (a) betanin, a derivative of the betalain pigment; (b) cyanidin; (c) luteolin; (d) quercetin; (e) chlorophyll; (f) curcumin.
Molecules 29 01680 g001
Figure 2. Illustration of the computational study flow to study the interaction between natural pigments and metals.
Figure 2. Illustration of the computational study flow to study the interaction between natural pigments and metals.
Molecules 29 01680 g002
Figure 3. Illustration of Molecular Dynamics Simulation Flow to study the interaction between natural pigments and metals.
Figure 3. Illustration of Molecular Dynamics Simulation Flow to study the interaction between natural pigments and metals.
Molecules 29 01680 g003
Table 1. The maximum permissible limit of metals based on several regulatory agencies.
Table 1. The maximum permissible limit of metals based on several regulatory agencies.
ElementsLimits (mg L−1)
FDAWHOEU
Heavy metals
Antimony0.0060.020.005
Arsenic0.010.010.01
Cadmium0.0050.0030.005
Lead0.0050.010.01
Mercury-0.006-
Non-heavy metals
Chromium0.100.050.05
Cobalt-0.05-
Copper1.002.002.00
Iron0.302.000.20
Manganese0.050.400.30
Nickel0.100.020.02
Zinc5.003.005.00
Note: EU, European Union; USFDA, United States Food and Drug Administration; WHO, World Health Organisation.
Table 2. Research related to the use of natural pigments as metal indicators.
Table 2. Research related to the use of natural pigments as metal indicators.
PigmentMetalColor ChangeLimit of Detection (Visually Visible)Reference
Chlorophyll-based
silver nanoparticle
HgBrown to light brown or colourless60 µM[58]
Curcumin-anthocyanin (hydrogel strips)CdWhite to bluish green0.2 µM[59]
HgWhite to blue0.2 µM
Curcumin
(cellulose acetate sensor strip)
PbYellow to orange20 µM[60]
Curcumin
(cellulose nanofiber)
PbOrange to red9 µM[61]
Curcumin-gold nanoparticleHgReddish wine to light blue2–10 µM[62]
Curcumin
(zein membrane)
FeYellow to brown7.16 µM[56]
Cyanidin
(solution)
AlPurple to violet to blue50 µM[57]
CuViolet to blue50 µM
FePink to violet to blue200 µM
PbPurple to violet to blue80 µM
Cyanidin
(dipstick sensor)
FeWhite to pink179–7162 µM[63]
Table 3. Research on metal–natural pigment interactions have used ab initio methods.
Table 3. Research on metal–natural pigment interactions have used ab initio methods.
MethodApproachPurposePigmentMetalCompatibility with
Experimental
Results
Reference
Ab initioHF method with the
6-31G(d) basis set
Optimise structures in
various conformations
Flavone and
flavylium
-Compatible in terms of internal rotation barrier[99]
Ab initioHF method with the
6-31G basis set
Calculate molecular
(hyper) polarisabilities and band structures
PorphyrinMg, Ni, Zn-[100]
Ab initioHF method with the
6-31G* basis set
Calculate fully optimised structures and atomic chargesBacteriochlorinMgCompatible in terms of structures, slightly differences in terms of
transition energy
[101]
IEFPCM-ab initioIEFPCM at the HF/6-31G(d)//mPW1PW91/6-31G(d) levelCalculate solvation free
energies
Flavylium--[102]
Ab initioHF method with the
6-31G* basis set
Full geometry
optimisation, study the
influence of methylation position, and study the
influence of additional functional groups in the
induction of atomic charge distribution
QuercetinMgCompatible in terms of binding free
energy
[103]
Luteolin Mg
Ab initioHF method with the
6-31G basis set
Electronic characterisationAnthracene,
naphthalene,
naphthacene, and pentacene
--[104]
Note: HF, Hartree–Fock; QMCF-MD, quantum mechanical charge field–molecular dynamics; IEFPCM, Integral Equation Formalism for the Polarizable Continuum Model.
Table 4. Research on metal–natural pigment interactions that has used density functional theory methods.
Table 4. Research on metal–natural pigment interactions that has used density functional theory methods.
ApplicationApproachConditionPurposePigmentMetalComplexBE
(kcal mol−1)
ΔBE
(kcal mol−1)
ΔE
(kcal mol−1)
ΔG
(kcal mol−1)
ΔG (binding)
(kcal mol−1)
Molecular
Energies
(Hartree)
Ground State
Energies
(kcal mol−1)
Compatibility with
Experimental
Results
Reference
Gaussian 09DFT/M052X/6-31+G(d) basis set (C, O, Al and H) and atoms
- relativistic compact Stuttgart/Dres den effective core potential in conjunction with its split valence basis set for Cu and Fe atoms
Aqueous phaseDetermine coordination propertiesQuercetinAl, Fe, and Cu(4–5) Al(OH)H4Que+ 13.70 Compatible[17]
(3′-4′) AlH3Que+ 0.04.5
(3′-4′s-cis) Fe(OH)2H3Que 0.6
(4–5) CuH4Que+ 28.20
(3′-4′H) CuH4Que+ 53.44.7
(3–4s-cis) CuH4Que+ 0.02
Gaussian 03DFT B3LYP/6-31+G(d)/LANL2DZ followed by single-point calculations using
the 6-311+G(2d,2p) basis set
Gas phaseLocate the exact chelation siteQuercetin (3,30,40,5,7-pentahydroxylflavone)CrNatural quercetin-bare Cr(III) ion1238.703333 Compatible[72]
Ethanol phaseNatural quercetin-bare Cr(III) ion668.98
Gas phaseDeprotonated quercetin-bare Cr(III) ion1514.885
Ethanol phaseDeprotonated quercetin-bare Cr(III) ion721.735
Gas phaseDeprotonated quercetin-hydrated Cr(III)469.3525
Ethanol phaseDeprotonated quercetin-hydrated Cr(III)78.4225
Gaussian 16DFT/M05–2X/6-311+G(d,p) Study chelating propertiesApigeninCu and Feapigenin/H3A − Cu(II) −2.2 -[80]
apigenin/H2A − Cu(II) −4.6
apigenin/HA2− − Cu(II) −6.9
apigenin/H3A − Fe(III) −7.0
apigenin/H2A − Fe(III) −12.9
apigenin/HA2− − Fe(III) −18.2
Gaussian 09DFT/UB3LYP/6-31+G(d,p)/LANL2DZGas phaseGeometry
optimisation
QuercetinCui-quercetin − Cu(II) −18.00 (kJ mol−1) Compatible[118]
Solution phase −41.03 (kJ mol−1)
Gas phaseii-quercetin − Cu(II) −9.31 (kJ mol−1)
Solution phase −43.29 (kJ mol−1)
Gas phaseiii-quercetin − Cu(II) 334 (kJ mol−1)
Solution phase −22.44 (kJ mol−1)
Gas phaseChrysinchrysin − Cu(II) −8.31 (kJ mol−1)
Solution phase −33.22 (kJ mol−1)
Gaussian 03DFT/B3LYP/6-31g**/CPCMMethanol and water phaseInvestigate
interactions and geometry optimisation
ApigeninAlAl1B4,5Ap2+ −4.95 Compatible[119]
Al1B4,5Ap12+ 14.10
LuteolinAl1B4,5Lu12+ 14.51
DMol3 codeMD/DFT/GGA/PBE functional and DNP basis setSolution phase Investigate metal
cation–pigment interactions
LuteolinZnZn(luteolin)2 (hydroxyl complexation) 2.585 Compatible[120]
Zn(luteolin)3 (hydroxyl complexation) 2.550
Zn(luteolin)4 (hydroxyl complexation) 2.526
Zn(luteolin)5 (hydroxyl complexation) 2.438
Zn(luteolin)2 (carbonyl complexation) 2.035
Zn(luteolin)3 (carbonyl complexation) 2.093
Zn(luteolin)4 (carbonyl complexation) 2.170
Zn(luteolin)5 (carbonyl complexation) 2.131
Gaussian 16DFT/B3LYP/6-31+G(d,p)/LANL2DZ Structure
optimisation
NaphthoquinonePd, Ni, and Co[Co(L1)2(H2O)2].2H2O 0.29 (eV)−581.21 kJ mol−1 -[121]
[NiL1(H2O)2(CH3COO)] 0.18 (eV)−631.88 kJ mol−1
[PdL1 (H2O)Cl] 0.39 (eV)−501.43 kJ mol−1
[CoL2(H2O)2(CH3COO)] 0.32 (eV)−434.27 kJ mol−1
[NiL2(H2O)2
(CH3COO)].H2O
0.71 (eV)−414.27 kJ mol−1
[PdL2 (H2O)Cl] 0.34 (eV)−301.43 kJ mol−1
Gaussian 09DFT/B3LYP/6-31++G(d,p)/PCMAqueous phaseStructure
optimisation
CyaninNa+ −12.21 1.16 Compatible[82]
K+ −6.36 10.58
Mg(II) −25.26 −13.50
Ca(II) −17.62 −2.71
Cr(II) −42.46 −28.22
Mn(II) −30.69 −17.85
Fe(II) −32.40 −19.75
Co(II) −31.26 −19.56
Ni(II) −32.40 −20.02
Cu(II) −52.04 −37.59
Zn(II) -31.74 −19.65
Al(III) −63.42 −53.97
Cr(III) −72.62 −63.54
Fe(III) −80.64 −73.44
Co(III) −108.66 −100.94
Gaussian 98DFT/B3LYP/6-31G* basis set for C and H atoms and 6-31+G* basis set for O atoms/LanL2DZ/PCM continuum modelSolution with ε = 78.4 (corresponding to bulk water)Structure optimisationQuercetinCuOQ1Cu −302.7 Compatible[122]
OQ2Cu −336.1
OQ3Cu −308.2
SQ1Cu −322.8 (−312.5)
SQ2Cu −323.3 (−328.8)
SQ3Cu −302.7 (−301.5)
D1Cu −500.2
D2Cu −499.8
DD3Cu −676.7
Gaussian 03DFT/M052 × 29/6-31+G(d) b/PCMEthanol
(ε = 24.85)
For full optimisationQuercetinAla3-4eq H4QueAl(H2O)2(OH)2 2.9 Compatible[73]
a3-4ax H4QueAl(H2O)2(OH)2 3.6
a4-5eq H4QueAl(H2O)2(OH)2 1.2
a4-5ax H4QueAl(H2O)2(OH)2 0
b3-4eq H3QueAl(H2O)3(OH) 5.1
b3-4ax H3QueAl(H2O)3(OH) 7.4
b4-5eq H3QueAl(H2O)3(OH) 5.9
b4-5ax H3QueAl(H2O)3(OH) 7.4
Gaussian 03DFT/B3LYP/6-31+G (d)/LANL2DZ followed by a single-point calculation using a different basis set (6-311++G(d,p))Water phaseGeometry optimisationChalcone (butein)Mg, Cr, Fe, and CuCu2+–O2′O9127 -[117]
Cu2+–O4O3125
Fe2+–O2′O932
Fe2+–O4O321
Mg2+–O2′O930
Mg2+–O4O322
Cr2+–O2′O946
Cr2+–O4O326
Gas phaseChalcone (butein)Mg, Cr, Fe, and CuCu2+–O2′O9517 0.03187 (Hertee)
Cu2+–O4O3493
Fe2+–O2′O9461 0.0716 (Hertee)
Fe2+–O4O3406
Mg2+–O2′O9423 0.06548 (Hertee)
Mg2+–O4O3370
Cr2+–O2′O9462 0.09428 (Hertee)
Cr2+–O4O3397
Gaussian09W and GaussView 6.0.16DFT/6-31+G(d)/B3LYP Investigate stability, reactivity, nature of interaction, and application of the complexes Quercetin (5-hydroxy-4-keto group)Al 2.8297 (eV) Compatible[65]
Mg 1.2679 (eV)
Na 3.1149 (eV)
K 3.2172 (eV)
Ca 0.6134 (eV)
Al 2.603
(eV)
Mg 1.3749 (eV)
Na 3.3362 (eV)
K 3.0782 (eV)
Ca 1.3983 (eV)
Quercetin (O3′/O4′
ortho-dihydroxyl (catechol) group)
Al 2.8828 (eV)
Mg 2.3431 (eV)
Na 2.8488 (eV)
K 3.1038 (eV)
Ca 1.156 (eV)
Al 1.051 (eV)
Mg 1.076 (eV)
Na 1.739 (eV)
K 1.697 (eV)
Ca 0.854 (eV)
Gaussian 16DFT/M06-2X/def2-SVPGas phaseInvestigate electronic and structural properties of morinQuercetinFe and CuCu(II)M2 6 m −2,411,782Compatible[123]
Fe(III)M2 6 m −2,175,328
Cu(II)Q2 6 m −2,411,784
Fe(III)Q2 6 m −2,175,331
Cu(II)M2 5 m −2,411,780
Fe(III)M2 5 m −2,175,325
Cu(II)Q2 5 m −2,411,769
Fe(III)Q2 5 m −2,175,326
Gaussian 09DFT/B3LYP/6-31+G-(d, p) Structural analysisCyanidinZn 1.47 (eV) -[124]
Gaussian 03DFT/B3LYP/6-31G*/LANL2DZ followed by single-point calculations using the extended 6-311++G** basis setGas phaseGeometry optimisationQuercetinFeI-Q-Fe2+ 12.3 Compatible[76]
II-Q-Fe2+ 10.1
III-Q-Fe2+ 0
I-Q-Fe2+ 27.4
II-Q-Fe2+ 23.6
III-Q-Fe2+ 5.6
IV-Q-Fe2+ 0
V-Q-Fe2+ 34.2
VI-Q-Fe2+ 35.3
VII-Q-Fe2+ 48.3
VIII-Q-Fe2+ 41.7
I-Q-Fe2+(H2O)4 13.8
II-Q-Fe2+(H2O)4 12.9
III-Q-Fe2+(H2O)4 0
IV-Q-Fe2+(H2O)4 1.8
V-Q-Fe2+(H2O)4 41
VI-Q-Fe2+(H2O)4 50.5
VII-Q-Fe2+(H2O)4 50.5
VIII-Q-Fe2+(H2O)4 42.3
I-2Q-Fe2+ 0
II-2Q-Fe2+ 5.9
III-2Q-Fe2+ 1.6
IV-2Q-Fe2+ 2.5
I-2Q-Fe2+(H2O)2 12.9
II-2Q-Fe2+(H2O)2 5.2
III-2Q-Fe2+(H2O)2 0
IV-2Q-Fe2+(H2O)2 9.3
Gaussian 16DFT/M05-2X/6-31+G(d) (C, H, O, and Al) and relativistic
compact Stuttgart/
Dresden
effective core potential with its related split valence (Cu and Fe)
Aqueous phaseProperties and geometry optimisationLuteolinAl, Fe, and CuAl(OH)2(H2O)2(H3Lu) (3′-4′) 4.9−89.9 -[64]
Al(OH)2(H2O)2(H3Lu) (4–5) 0−94.8
Fe(OH)2(H2O)2(H3Lu) (3′-4′) 1.4−82.4
Fe(OH)2(H2O)2(H3Lu) (4–5) 0−83.9
[Fe(H2O)3(OH)(H3Lu)](3′-4′) 2.8−61.7
[Fe(H2O)3(OH)(H3Lu)](4–5) 0−64.4
Cu(OH)2(H2O)2(H3Lu)(3′-4′) 0.2−61.1
Cu(OH)2(H2O)2(H3Lu)(4–5) 0−61.2
Gaussian 09DFT/M052X/6-31+G(d)Water phaseStudied the complexationCurcuminAl and Fe[Al(H2O)3(OH)(LA)]+ −135.1 Compatible[125]
[Al(H2O)3(OH)(LB)] −124.9
[Fe(H2O)(OH)3(LA)] −57.1
[Fe(H2O)(OH)3(LB)] −55.5
Gaussian 09DFT/U-B3LYP/6-31G* and relativistic effective core potential with a valence basis set/LANL2DZ Coordination and geometry optimisationQuercetinNi[Ni(L1)(fla)]ClO4 3.056 (eV) Compatible[71]
[Ni(L1)(fla)]ClO4 3.045 (eV)
[Ni(L3)(fla)]ClO4 3.029 (eV)
[Ni(ntb)(fla)]+ 3.033 (eV)
Gaussian 09DFT/UB3LYP*/6-311++G(d,p) Full geometry optimisationPhenoxazinesCoI(R = R1 = CH3)LSCoIII-SQ 4.5 -[70]
I(R = R1 = CH3)HSCoIII-Q
I(R = CH3, R1 = CF3)LSCoIII-SQ −1.7
I(R = CH3, R1 = CF3)HSCoIII-Q
I(R = R1 = CF3)LSCoIII-SQ −8.7
I(R = R1 = CF3)HSCoIII-Q
II(R2 = H)LSCoIII-SQ 14.4
II(R2 = H)HSCoIII-Q
II(R2 = CH3)LSCoIII-SQ 12.1
II(R2 = CH3)HSCoIII-Q
II(R2 = Ph)LSCoIII-SQ 3.4
II(R2 = Ph)HSCoIII-Q
III(R2 = H)LSCoIII-SQ 14.4
III(R2 = H)HSCoIII-Q
III(R2 = CH3)LSCoIII-SQ 11.9
III(R2 = CH3)HSCoIII-Q
III(R2 = Ph)LSCoIII-SQ 3.3
III(R2 = Ph)HSCoIII-Q
Note: BE, energy associated with the formation of bonds between two or more molecules in a system; ΔBE, change in energy associated with the formation of bonds between two or more molecules in a system; ΔE, change in molecular system energy; ΔG, change in Gibbs free energy of the molecular system; ΔG (binding), change in Gibbs free energy associated with the formation of a bond between two or more molecules.
Table 5. Research on heavy metal–natural pigment interactions that has used density functional theory methods.
Table 5. Research on heavy metal–natural pigment interactions that has used density functional theory methods.
ApplicationApproachConditionPurposePigmentMetalComplexEE
(eV)
ΔE
(eV)
ΔEDIS
(kJ mol−1)
Molecular
Energy
(Hartree)
Compatibility with
Experimental Results
Reference
Gaussian 09DFT/B3LYP/6-
31G/LANL2DZ
Determine the
reactivity and stability of the complex
PorphyrinPb, Cd, Hg, Sn, and AsPb-TMPyP −2138.3Compatible[139]
Cd-TMPyP −2182.9
Hg-TMPyP −2177.5
Sn-TMPyP −2137.1
As3+-TMPyP −2140.4
As5+-TMPyP −2138.9
Pb-TDMImP −2207.3
Cd-TDMImP −2251.9
Hg-TDMImP −246.5
Sn-TDMImP −2205.9
As3+-TDMImP −2209.3
As5+-TDMImP −2207.8
Pb-TDMPzP −2207.2
Cd-TDMPzP −2251.8
Hg-TDMPzP −2246.4
Sn-TDMPzP −2206.0
As3+-TDMPzP −2209.3
As5+-TDMPzP −2207.9
Gaussian 03DFT/B3LYP/6-31G(d,p)/LANL2DZMethanol phaseGeometry optimisationQuercetinPbPb(II)-quercetin 0.124 Compatible[135]
Gaussian 09DFT/B3LYP/LANL2DZ followed by TD-DFT/LC-wPBEGas phaseGeometry optimisation and energy
calculations
AstaxanthinPb, Cd, and Hg[ASTA-Pb]+22.05 Compatible[128]
[ASTA-Pb2]+41.84
[ASTA-Cd(H2O)2]+21.85
[ASTA-Cd2(H2O)4]+42.08
[ASTA-Hg(H2O)2]+21.93
[ASTA-Hg2(H2O)4]+42.06
Ethanol phase[ASTA-Pb]+21.96
[ASTA-Pb2]+41.82
[ASTA-Cd(H2O)2]+21.68
[ASTA-Cd2(H2O)4]+42.41
[ASTA-Hg(H2O)2]+21.82
[ASTA-Hg2(H2O)4]+42.40
Gaussian 03DFT/B3LYP/6-311G(d,p)/LANL2DZ Calculation of
complete optimisation
Pterins,
isoxanthopterin,
sepiapterin
Cd, HgPterins-Cd 7.5 Compatible[140]
Pterins-Cd+1 270.3
Pterins-Cd+2 762.7
Pterins Hg 4.6
Pterins-Hg+1 252.7
Pterins-Hg+2 797.1
Isoxanthopterin-Cd 7.5
Isoxanthopterin-Cd+1 276.6
Isoxanthopterin-Cd+2 778.2
Isoxanthopterin-Hg 5.9
Isoxanthopterin-Hg+1 258.6
Isoxanthopterin-Hg+2 807.5
Sepiapterin-Cd 12.9
Sepiapterin-Cd+1 366.9
Sepiapterin-Cd+2 1018.4
Sepiapterin-Hg 9.2
Sepiapterin-Hg+1 346.4
Sepiapterin-Hg+2 1032.2
EE, energy related to the excitation process or energy given to a system related to changes in state from the ground state to the excited state; ΔE, change in molecular system energy; ΔEDIS, change in energy during the dissociation process.
Table 6. Research on metal–natural pigment interactions that has used semiempirical methods.
Table 6. Research on metal–natural pigment interactions that has used semiempirical methods.
MethodApproachPurposePigmentMetalCompatibility with
Experimental
Results
Reference
Semiempirical methodAM1 using HyperChem programStudy complexation processesAnthocyaninAlCompatible[142]
Semiempirical methodMM+ and AM1Molecular calculationAnthocyaninAl and GaCompatible[143]
Semiempirical methodAM1 Hamiltonian Calculate the structural modifications caused by the chelation processQuercetin AlCompatible[74]
Semiempirical methodZINDO/S CIS (40,40) or (45, 45) levelsCalculate transition energies and oscillation strengths, estimate corresponding spectroscopic transition energy values, and study the existence of dark electronic states in the systemBacteriochlorinMgCompatible[106]
Semiempirical methodPM5 Study orbitals Bacteriochlorin, chlorin, and porphinMgCompatible[141]
Semiempirical methodAM1 using MOPAC17 version 6Geometry optimisation and electronic structure analysisXanthophylls, antheraxanthin, lutein, neoxanthin,
violaxanthin, and zeaxanthin
-Compatible[144]
Semiempirical methodPM5 using the MOPAC 2002 packageStructure optimisationChlorophyll and chlorophyll d peptides MgCompatible[145]
Semiempirical methodAM1 (available in the AMPAC package)Structure optimisationMalvidin-Compatible[146]
Semiempirical molecular dynamicsPM6 with Grimme D3 dispersion correctionInvestigate intermolecular behaviourAnthocyanins (malvidin-3-glucoside)Al and SnCompatible[147]
ONIOM semiempirical methodPM6 and Pm3MMSolvation optimisationAnthocyaninsAl, Ga, Cr, Fe, and Mg-[138]
QM/MM semiempirical methodPm3MMSolvation optimisationAnthocyaninsMg, Al, Ga, Sn, Cr, and FeCompatible[78]
Notes: ONIOM, Our Own N-layered Integrated molecular Orbital and molecular Mechanics; QM/MM, Quantum Mechanics/Molecular Mechanics.
Table 7. Research on metal–natural pigment interactions that has used molecular dynamics simulation.
Table 7. Research on metal–natural pigment interactions that has used molecular dynamics simulation.
MethodApproachPurposePigmentMetalCompatibility with
Experimental
Results
Reference
Molecular
dynamics
simulations
-Calculate energy-minimised structuresHesperidin, rutin,
neodiosmin, diosmin, and neohesperidin
CoCompatible[153]
Molecular
dynamics
simulations
AMBER force field (GAFF) and TIP3P model using the Sander module in the Amber 10.0 simulation packageIdentify several
conformations of the
co-pigmentation complex
Oenin-Compatible[154]
Molecular
dynamics simulations
restrained electrostatic
potential (RESP) protocol
implemented in the ANTECHAMBER module of AMBER 11
Partial atomic charges
determination
Quercetin and luteolinMgCompatible[108]
AMBER force field (GAFF)Describe the force field
parameters of the
substrates
AMBER ff99SB force fieldPotential determination
Ab initio QMCF-MD simulationsHF method with the
6-31G** basis set; Amber force field (GAFF) and restrained electrostatic potential (RESP)
Study structural and
dynamic properties, and observe hydration
behaviour
PorphyrinMgCompatible[155]
Molecular
dynamics
simulations
AMBER14SB force field and GAFF force fieldEvaluate the possibility of increasing or decreasing internuclear distance upon low-energy conformational changes and assessing the flexibility of the complexPorphyrinAgCompatible[156]
Monte Carlo/
molecular
dynamics
simulations
Materials Studio in the NVT ensemble with the COMPASS force fieldAnalyse the interaction of complexesLuteolin–Zn complexFeCompatible[120]
Monte Carlo/
Molecular
dynamics
simulations
AMBERStudy conformation
theories, thermodynamic parameters and movement rules of the molecular
machine and kinetic
energy to the potential
energy surface
Cyanidin, delphinidin, petunidinMg, Al, Ga, Sn, Cr, and FeCompatible[78]
Molecular
dynamics
simulations
Reax FF force field using NVT ensembleAnalyse interactions
between components
Cyanidin-3-glucosideTiCompatible[157]
Molecular
dynamics
simulations
Forcite module of a Material Studio software in NVT
ensemble employing universal forcefield
Study adsorption
behaviour
Satureja hortensis extract (isoferulic acid, caffeic acid, kuersetin,
rosmarinic acid,
apigenin glucoside, and chlorogenic acid)
ZnCompatible[158]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Maranata, G.J.; Megantara, S.; Hasanah, A.N. An Update in Computational Methods for Environmental Monitoring: Theoretical Evaluation of the Molecular and Electronic Structures of Natural Pigment–Metal Complexes. Molecules 2024, 29, 1680. https://doi.org/10.3390/molecules29071680

AMA Style

Maranata GJ, Megantara S, Hasanah AN. An Update in Computational Methods for Environmental Monitoring: Theoretical Evaluation of the Molecular and Electronic Structures of Natural Pigment–Metal Complexes. Molecules. 2024; 29(7):1680. https://doi.org/10.3390/molecules29071680

Chicago/Turabian Style

Maranata, Gabriella Josephine, Sandra Megantara, and Aliya Nur Hasanah. 2024. "An Update in Computational Methods for Environmental Monitoring: Theoretical Evaluation of the Molecular and Electronic Structures of Natural Pigment–Metal Complexes" Molecules 29, no. 7: 1680. https://doi.org/10.3390/molecules29071680

APA Style

Maranata, G. J., Megantara, S., & Hasanah, A. N. (2024). An Update in Computational Methods for Environmental Monitoring: Theoretical Evaluation of the Molecular and Electronic Structures of Natural Pigment–Metal Complexes. Molecules, 29(7), 1680. https://doi.org/10.3390/molecules29071680

Article Metrics

Back to TopTop