Next Article in Journal
Phytochemical Investigation and Biological Activities of Lantana rhodesiensis
Next Article in Special Issue
Research Progress on Synthesis and Application of Cyclodextrin Polymers
Previous Article in Journal
Lignins Isolated via Catalyst-Free Organosolv Pulping from Miscanthus x giganteus, M. sinensis, M. robustus and M. nagara: A Comparative Study
Previous Article in Special Issue
Inclusion Compound of Efavirenz and γ-Cyclodextrin: Solid State Studies and Effect on Solubility
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Water Microsolvation on the Excited-State Proton Transfer of 3-Hydroxyflavone Enclosed in γ-Cyclodextrin

by
Khanittha Kerdpol
1,
Rathawat Daengngern
2,3,
Chanchai Sattayanon
4,
Supawadee Namuangruk
4,
Thanyada Rungrotmongkol
5,6,7,
Peter Wolschann
8,9,
Nawee Kungwan
10,11,* and
Supot Hannongbua
1,7,*
1
Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
2
Department of Chemistry, Faculty of Science, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand
3
Integrated Applied Chemistry Research Unit, Faculty of Science, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand
4
National Nanotechnology Center (NANOTEC), NSTDA, 111 Thailand Science Park, Pahonyothin Road, Klong Luang, Pathum Thani 12120, Thailand
5
Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
6
Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok 10330, Thailand
7
Molecular Sensory Science Center, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
8
Department of Pharmaceutical Chemistry, University of Vienna, 1090 Vienna, Austria
9
Institute of Theoretical Chemistry, University of Vienna, 1090 Vienna, Austria
10
Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
11
Center of Excellence in Materials Science and Technology, Chiang Mai University, Chiang Mai 50200, Thailand
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(4), 843; https://doi.org/10.3390/molecules26040843
Submission received: 10 January 2021 / Revised: 31 January 2021 / Accepted: 2 February 2021 / Published: 5 February 2021
(This article belongs to the Special Issue Cyclodextrin Chemistry and Toxicology)

Abstract

:
The effect of microsolvation on excited-state proton transfer (ESPT) reaction of 3-hydroxyflavone (3HF) and its inclusion complex with γ-cyclodextrin (γ-CD) was studied using computational approaches. From molecular dynamics simulations, two possible inclusion complexes formed by the chromone ring (C-ring, Form I) and the phenyl ring (P-ring, Form II) of 3HF insertion to γ-CD were observed. Form II is likely more stable because of lower fluctuation of 3HF inside the hydrophobic cavity and lower water accessibility to the encapsulated 3HF. Next, the conformation analysis of these models in the ground (S0) and the first excited (S1) states was carried out by density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations, respectively, to reveal the photophysical properties of 3HF influenced by the γ-CD. The results show that the intermolecular hydrogen bonding (interHB) between 3HF and γ-CD, and intramolecular hydrogen bonding (intraHB) within 3HF are strengthened in the S1 state confirmed by the shorter interHB and intraHB distances and the red-shift of O–H vibrational modes involving in the ESPT process. The simulated absorption and emission spectra are in good agreement with the experimental data. Significantly, in the S1 state, the keto form of 3HF is stabilized by γ-CD, explaining the increased quantum yield of keto emission of 3HF when complexing with γ-CD in the experiment. In the other word, ESPT of 3HF is more favorable in the γ-CD hydrophobic cavity than in aqueous solution.

Graphical Abstract

1. Introduction

Fluorescent organic molecules with strong intramolecular hydrogen bonds (intraHBs) connected by proton donating and accepting groups have gained considerable attention in recent years owing to their unique fluorescent emission in term of large Stokes shift without self-absorption [1]. Their unique properties harnessing from the excited state intramolecular proton transfer (ESIntraPT) could be typically described by the characteristic four-level photocycle. Initially, the molecule in an enol form (E) in the ground (S0) state absorbing light in the shorter wavelength region results in the photoexcitation process from the S0 state into the excited (S1) state. The intraHB of E is strengthened in the S1 state because of the charge redistribution upon photoexcitation, leading to the transfer of proton from the donor (D: −NH2, −OH) to the acceptor (A: C=O, −N=), which changes the enol form (E*) to the keto form (K*) in the S1 state. After that, the K* emits the fluorescence in the remarkably longer wavelength than the absorption and relaxes to the S0 state, resulting in the notably large Stokes shift (the difference between positions of absorption and emission peaks). Then the K changes to the E through the back proton transfer (BPT) process spontaneously in the S0 state due to the PT barrierless and high exothermic reaction. Generally, their photophysical properties can be easily modulated using many strategies [2,3] such as introducing electron-donating and withdrawing substituents, the heteroatom substitution, and the π-conjugation into the main core structure, to give the desirable absorption and emission spectra as well as large Stokes shift. The ESIntraPT molecules with such tunable photophysical properties including derivatives of salicylates [4,5,6], salicylideneanilines [7,8,9], flavones [10,11,12,13], benzazoles [14,15,16,17,18], and chalcones [19,20,21] have been reported and widely used in various applications ranging from chemical sensing to light-emitting diodes [22,23,24,25,26,27].
Among various ESIntraPT molecules, 3-hydroxyflavone (3HF), which consists of a chromone ring (C-ring) and a phenyl ring (P-ring), is one of the best-known molecular systems. 3HF exhibits ESIntraPT and gives dual fluorescence corresponding to its E* and K* forms with a large Stokes shift and photostability [28,29,30]. Thus, 3HF has been used as a prototype for the ESIntraPT processes and as sensitive fluorescence probes for discovering binding sites in various bio-relevant targets such as DNA, protein, and bio-membranes [31,32,33]. The photophysical properties and ESPT processes of 3HF in organic solvents have been extensively studied [34,35,36,37,38,39,40,41,42]. In non-polar solvent, only the K* emission peak of 3HF in toluene was observed at 530 nm with large Stokes shift [33] because the ESIntraPT process effectively occurs, giving only K* form. However, in a protic solvent, the dual emission peaks from E* and K* of 3HF in methanol were observed at 409 and 528 nm, respectively [33], because the IntraHB of 3HF is disrupted and the intermolecular hydrogen bonds (interHBs) between 3HF and protic solvents is formed depending on the nature of solvents and the arrangement of protic solvent around 3HF. This favorable formation of interHBs could reduce the formation of K*, resulting in the low quantum yield in protic solvents or aqueous solution [40,41,42,43,44].
One of the strategies to enhance the K* emission intensity of 3HF in an aqueous solution is disrupting the interHBs of 3HF-water cage-like network using cyclodextrin (CD) [45,46,47]. The reduction of polarity and restricted environment inside CD’s cavity is essential for many aspects of photophysical phenomena by inclusion complexes between 3HF and CD [45,46,47]. CDs are the cyclic oligosaccharides consisting of the crucial α-d-glucose unit, which exhibit the conical shape and the α-d-glucose units of 6, 7, and 8 are represented as α-, β-, and γ-cyclodextrins (α-, β-, and γ-CDs), respectively [48]. S. Das and N. Chattopadhyay experimentally studied the fluorescence anisotropy of the inclusion complexes of 3HF in α-, β-, and γ-CDs compared to 3HF in the aqueous medium. The fluorescence anisotropy of these probes decreased following the order α-, β-, and γ-CD, which is attributed to the disruption of the 3HF-water network in an aqueous medium [46]. It can be stated that the micro-environment of 3HF derivatives was able to alter and prevent the self-aggregation effectively by using CDs, especially γ-CD. From the investigation on ESPT processes of 3HF in β-, and γ-CDs [45,47], the intensity of K* of 3HF/γ-CD inclusion complex is higher than that of 3HF/β-CD inclusion complex. Consequently, the encapsulation in CDs should be a possible method to tune the fluorescence emission of 3HF derivatives and other hydrophobic compounds [45,46,47,49,50,51,52].
From the studies of the multi-spectroscopic approaches and molecular docking of the encapsulation of 3HF in different small-ring CDs, the 3HF/γ-CD complex showed the strongest interaction, and it provided a higher fluorescent yield than that in water medium [45,46,47]. Understanding the role of CD in enhancement of the fluorescent yield at atomic level might help us to be able to adjust the fluorescent wavelength to fit for fluorescence probes for bio-labeling in aqueous medium. However, to the best of our knowledge, the detailed information of ESPT reaction in the S0 and S1 states of 3HF in γ-CD at atomic level has not been reported. In this work, we aimed to systematically investigate the effect of a water molecule on the photophysical properties and ESPT reactions of an isolated 3HF and its encapsulation. All-atom molecular dynamics (MD) simulations for 300 ns were performed to study the structure and dynamics properties of the two possible 3HF/γ-CD complexes. The detailed information of each complex both in the S0 and S1 states was then investigated using density functional theory (DFT) and time-dependent DFT (TD-DFT) methods. The important distances and simulated infrared (IR) vibrational spectra from the optimized structures as well as the topology analysis were used to describe the hydrogen-bonded strength. The frontier molecular orbitals (MOs) were analyzed to provide the charge distribution of the complex. The simulated absorption and emission spectra were calculated and compared with the experimental data. Moreover, the energies of E and K forms of 3HF in each system at the S0 and S1 states were discussed to explain why the fluorescent yield of K* in aqueous medium increases when encapsulating 3HF into γ-CD.

2. Results and Discussion

2.1. Possible Inclusion Complexes

From the docking study, 3HF could form two possible inclusion complexes with γ-CD through its chromone ring (C-ring, Form I) or phenyl ring (P-ring, Form II) insertion into the hydrophobic cavity as depicted in Figure 1B. Although a higher occurrence (65%) was found for Form I, the interaction energies of both forms were likely comparable (Form I: −21.42 kcal/mol and Form II: −20.78 kcal/mol). Additionally, the previous docking studies [46,47] suggested that Form II was more stable. Thus, in the present work the 3HF/γ-CD complexes in both forms were further studied by MD simulations and DFT calculations.

2.2. 3HF Mobility in γ-CD Cavity and Water Accessibility

To study the inclusion complexes in solution, the four different initial structures of Form I and Form II obtained from molecular docking and QM calculation were simulated by 300-ns MD simulations (MD1-MD4). All trajectories were analyzed and discussed as follows. The plots of RMSD and Rgyr of complex in Supplemental Figure S1 suggested that in Form I 3HF spontaneously released from the γ-CD cavity at ~68 ns, ~54 ns, and ~198 ns for MD1, MD2 and MD3. Interestingly, it feasibly moved back to form a complexation with γ-CD, resulting in Form I (MD3 ~201 ns) and Form II (MD1 ~70 ns and MD2 ~57 ns) as considerably seen by the plot of distance between the center of mass (Cm) of each 3HF ring and the Cm of the primary rim of the γ-CD in Figure 2, and the plot of distance between the Cm of each 3HF ring and the Cm of the secondary rim of the γ-CD in Supplementary Materials Figure S2. This is in contrast for Form II, in which the 3HF was well encapsulated inside the hydrophobic cavity in all MD1-MD4 systems throughout the simulation time (RMSD of 1.2–4.0 Å, and Rgyr of 6.0−7.0 Å). Moreover, it can be noticeable that the both rings of 3HF in Form I (MD3 and MD4) were fluctuated higher than those in Form II, suggesting that the complex in Form II was found to be more stable in aqueous solution; in other words, the P-ring insertion is the suitable binding mode of 3HF for encapsulation with γ-CD in consistent with the lower water accessibility to the encapsulated 3HF (n(r) of 2.0 ± 0.7 and 2.5 ± 1.2 at O1 and O2, respectively, in Figure 3). In Form I, the hydroxyl oxygen O1 of 3HF positioned closer to the wider rim of γ-CD had a significantly higher interaction with waters (4.4 ± 0.8), and in vice versa less waters (2.0 ± 0.2) can access to the carbonyl oxygen O2. No peak detected within ~2.8 Å of the O3 of 3HF, suggesting that this atom had a very weak hydration interaction as found in some flavonoids/CD complexes [53,54]. It is worth noting that such accessible water molecules at O1 and O2 sites may involve into proton transfer processes either blocking ESIntraPT or assisting ESInterPT in 3HF/γ-CD inclusion complex. To study the water assisted PT in 3HF/γ-CD either ground state or excited state, the model of 3HF/γ-CD with a water molecule placed between these two sites of 3HF was further investigated by DFT and TD-DFT calculations and discussed in the next sections.

2.3. Structural Optimizations

All optimized E forms of 3HF, 3HFW, and 3HF encapsulated in the γ-CD cavity (Form I, Form II, Form I-W and Form II-W) with important labeled atoms and distances (O−H covalent bonds, intraHB, and interHBs) in the S0 state are shown in Figure 4, where the measured distances in the S0 and S1 states are summarized in Table 1.
From Table 1, for the compounds without a water molecule (3HF, Form I, and Form II), the O1–H1 covalent bond at the 3HF hydroxyl group slightly increases from the S0 to S1 states around 0.015–0.034 Å. Whereas, the length of O2⋯H1 intraHB between the hydroxyl group and carbonyl group of 3HF significantly decreases in the S1 state at 0.214, 0.174, and 0.097 Å for 3HF, Form I, and Form II, respectively. The length of O2⋯H1 intraHB of Form I and Form II is longer than that of the isolated 3HF because the O2 of 3HF forms a stronger interHB with a hydroxy group at the primary rim of γ-CD (O6–H) for Form I and at the secondary rim of γ-CD (O2–H) for Form II in the S1 state. Overall, these results indicate that the strength of an intraHB in the S1 state of 3HF, Form I, and Form II is stronger than that in the S0 state. Consequently, the ESIntraPT process might easily occur in the S1 state.
For the compounds with a water molecule (3HFW, Form I-W, and Form II-W), a water molecule forms interHBs with 3HF, therefore an O2⋯H1 intraHB of these complexes is dramatically longer than that of the compounds without a water molecule. It can be stated that a water molecule might block ESIntraPT process and support ESInterPT process. The O1–H1 and Ow–Hw covalent bonds of 3HFW/Form II-W are slightly increased from the S0 to S1 states at 0.069/0.074 Å for O1–H1 and 0.024/0.023 Å for Ow–Hw covalent bonds. While the lengths of Ow⋯H1 and O2⋯Hw interHBs of 3HFW/Form II-W are significantly decreased from the S0 to S1 states around 0.236/0.237 Å for Ow⋯H1 and 0.151/0.152 Å for O2⋯Hw interHBs, so ESInterPT may easily take place more than ESIntraPT for 3HFW and Form II-W due to the strong interHBs in the S1 state. For Form I-W, the length of Ow⋯H1 interHB decreases (0.082 Å), while the length of O2⋯Hw interHB increases (0.095 Å) in the S1 state because Hw of a water molecule forms a stronger interHB with O6 at a primary rim of γ-CD instead of O2 acceptor of 3HF in the S1 state. Consequently, ESInterPT of Form I-W is quite difficult to occur because the O2⋯Hw interHB between 3HF and a water molecule is weaker (elongated length) in the S1 state. Furthermore, the strength of intraHBs and interHBs in the S1 state will be further discussed in the next section.

2.4. O–H Stretching and Topology Analysis

In general, the strength of the intraHBs and interHBs in the S1 state could be revealed based on monitoring the red-shift of vibrational modes involving the hydrogen-bonded formation and the topology analysis at bond critical points (BCPs) using the QTAIM method. The calculated IR spectra of all studied compounds in the conjunct vibrational regions of the O–H stretching modes related with PT process both in the S0 and S1 states are listed in Table 2. These O–H stretching modes can be classified into the O1–H1 stretching mode of 3HF and Ow–Hw stretching mode of a water molecule.
For 3HF, the O1–H1 stretching vibrational mode of 3HF is located at 3544 cm−1 and 3003 cm−1 for the S0 and S1 states, respectively, giving a large red-shift of 541 cm−1. Moreover, the large red-shift of the O1–H1 stretching vibrational mode is also observed for Form I (367 cm−1) and Form II (244 cm−1). Therefore, the strength of the O2⋯H1 intraHB for these compounds is increased in the S1 state providing ESIntraPT process.
For the compounds with a water molecule, the O1–H1 stretching vibrational modes of 3HFW, Form I-W, and Form II-W are located around 2967–3129 cm−1 in the S0 state. Note that these vibrational modes changed to be around 2575–2599 cm−1 in the S1 state, which evidently demonstrates that the red-shift is induced by the strengthening of the O1–H1⋯Ow interHB after photoexcitation. Similarly, the Ow–Hw stretching modes of these compounds are also red-shifted. In addition, the red-shift value of 3HFW and Form II-W is larger than that of Form I-W, indicating that the O1–H1⋯Ow and Ow–Hw⋯O2 intermolecular hydrogen-bonded strength of 3HFW and Form II-W in the S1 state is stronger than that of Form I-W. Overall, these results show that O–H stretching vibrational frequencies shift to lower frequencies in the S1 state compared with the S0 state, which confirms that the hydrogen bonding interaction is stronger in the S1 state.
A topology analysis of the electron density was used to further determine the strength of the intraHB and interHBs in the excited-state structures (E* form) of all compounds. The following parameters at BCPs were analyzed: the electron density ρ(r), the potential energy density V(r), the Laplacian of the electron density ∇2ρ(r), the Lagrangian kinetic energy G(r), the Hamiltonian kinetic energy density H(r), and the electron delocalization index (DI) between the proton acceptor and transferred proton, which are an O2⋯H1 intraHB, and the Ow⋯H1 and O2⋯Hw interHBs for all studied compounds. Additionally, the hydrogen-bonded energy (EHB) can be calculated by using the Espinosa’s equation: E H B = 1 2 | V ( r B C P ) | . These results are summarized in Table S1 in the Supplementary Materials. From Table S1, EHB of O2⋯H1 intraHB of 3HF (0.0224 a.u.) is the highest among the compounds without a water molecule. Then, EHB is slightly decreased for Form I (0.0175 a.u.) and Form II (0.0111 a.u.). These results indicate that the intraHB of isolated 3HF is stronger than that the inclusion complexes. Nevertheless, this intraHB of all compounds can facilitate ESIntraPT processes. For the compounds with a water molecule, EHB of Ow⋯H1 and O2⋯Hw of 3HFW/Form II-W are 0.0696/0.0649 and 0.0306/0.0317 a.u., respectively. The Ow⋯H1 interHB of 3HFW/Form II-W is stronger than the O2⋯Hw interHB, indicating that the proton might transfer via a water molecule from O1–H1 bond-breaking before Ow–Hw bond-breaking. However, only EHB of Ow⋯H1 of Form I-W is obviously dropped (0.0253 a.u.). The result implies that the ESInterPT process might occur in 3HFW/Form II-W better than Form I-W. Overall, it can be observed that the intraHB and interHBs are strengthened in the S1 state confirmed by the shorter distances of important bonds involving in the ESPT process, the red-shift observed by IR vibrational spectral calculations, and a high value of EHB from topology analysis of the electron density.

2.5. Frontier MOs Analysis and Simulated Spectra

To further investigate behaviors of charge distribution and charge transfer in the S1 state, the frontier MOs of the highest occupied molecular orbitals (HOMO) and the lowest unoccupied molecular orbitals (LUMO) of all studied compounds were analyzed because the main electronic transition is only related with these orbitals in range of 98% (HOMO → LUMO), which is assigned as π to π* characters, and illustrated in Figure 5. It can be noted that electron density of both HOMO and LUMO is fully localized on the 3HF moieties and no electron density is located on water or γ-CD, indicating that no intramolecular charge transfer within 3HF and no intermolecular charge transfer between 3HF and water or γ-CD. Moreover, the HOMO and LUMO are localized on different parts of 3HF. For the HOMO orbitals, the electron density is distributed more on P-ring and partially on C-ring of 3HF. Whereas, that of LUMO is distributed completely on the whole molecule of 3HF.
The UV-Vis absorption and emission spectra of all studied compounds were simulated based on their optimized S0 and S1 structures, respectively and then plotted in Figure S3. Furthermore, the absorption band maxima of E form (λabs of E), emission band maxima of E* (λemis of E*) and K* forms (λemis of K*), the excitation energy (Eex), and the oscillator strength (f) as well as the major MOs contribution (%) of the absorption band for all compounds are reported in Table 3. From the detailed information in Table 3, the simulated absorption peaks of the complexes without and with a water molecule are around 341–345 nm and 350–353 nm, respectively, which are in good agreement with the experimental value of 341 nm for 3HF in water and 340 nm for 3HF in γ-CD [46,47]. Moreover, the predicted maximum wavelength for dual emission spectra of all studied compounds is also consistent with the experimental data [46,47], in which the λemis of E* in water and in γ-CD are reported at 410 and 404 nm, respectively, while the λemis of K* in water and in γ-CD are 511 and 538 nm. The deviations from the experimental data around 59–69 nm (0.52–0.61 eV) indicate that the chosen method at TD-PBE0/def2-SVP level of theory is adequate to describe the electronic spectra and provide the insight understanding of the ESPT process.
The relative energy of E and K forms of all studied complexes at the S0 and S1 states was investigated to explain the ESPT phenomena as illustrated in Table 4. The results show that the E form is more stable than the K form in the S0 state for all complexes with the energy differences at 6.63–16.85 kcal/mol. Moreover, most of K forms were stabilized in the γ-CD cavity except only Form I. However, in the S1 state, the K* form is more stable than the E* form for all complexes. It is predicted that both ESIntraPT and ESInterPT processes are favorable in the S1 state but not in the S0 state. In case of the complexes without a water molecule (Form I and Form II), K* of Form II is more stable than Form I. So, the ESIntraPT process of Form II may be more effective than that of Form I, related to the MD results that Form II is favorably more stable. For the complexes with a water molecule, the ESInterPT via interHBs network is feasible to occur in both Form I-W and Form II-W especially Form II-W because of the slightly higher oscillator strength and the lower energy of K*. In addition, from our previous work, the ESInterPT of 3HFW is hard to occur due to the higher ESInterPT barrier and the rearrangement of a water molecule surrounding 3HF [40]. Therefore, it can be predicted that the encapsulating 3HF into γ-CD assists the disruption of the 3HF-water network in aqueous solution leading to an increment of the fluorescent yield of K* in aqueous solution from the ESIntraPT process.

3. Methodology

3.1. Model Preparation

The 3D structure of 3HF (Figure 1A) was generated by GaussView version 6.0 [55]. To investigate the effect of water on PT process, a water molecule near the PT site of 3HF was placed between the H1 and O2 atoms to model 3HF with a water molecule (3HFW). The 3HF and 3HFW structures were fully optimized with DFT and TD-DFT methods using Gaussian 16, Revision C.01 [56]. The γ-CD structure extracted from the co-crystal structure of cyclo/maltodextrin-binding protein in complex with γ-CD (PDB code: 2ZYK) was used in this study. To obtain the inclusion complex, encapsulation of the optimized 3HF into γ-CD’s cavity was performed by molecular docking with 100 runs using the CDOCKER module in the Accelrys Discovery Studio 2.5 (Accelrys Software Inc., San Diego, CA, USA). The three docked inclusion complexes with the lowest CDOCKER interaction energy for the chromone ring (C-ring, Form I) and phenyl ring (P-ring, Form II) insertion into the hydrophobic cavity of the γ-CD were selected for MD simulations (MD1-MD3 for Form I and Form II). The most stable 3HF/γ-CD inclusion complexes in Form I and Form II from a docking study were fully optimized at PBE0/def2-SVP level of theory for investigating PT reaction. Note that, these Form I and Form II were used as the additional initial structures for MD simulations, MD4. To study the effect of water associated in PT for 3HF encapsulated in γ-CD, a water molecule was added to form interHBs with 3HF both in Form I and Form II namely Form I-W and Form II-W, respectively. Then, Form I-W and Form II-W were fully optimized with DFT and TD-DFT methods.

3.2. Molecular Dynamics Simulations

MD simulations of the 3HF/γ-CD complexes in Form I and Form II using the four different initial structures (eight complexes in total) were carried out using the AMBER16 program package [57]. The partial atomic charges and parameters of 3HF were generated in accordance with the previous standard procedures [58,59,60]. The general AMBER force field [61] and the Glycam06j carbohydrate force field [62] were applied for 3HF and γ-CD, respectively. The models were solvated by a truncated octahedral box of TIP3P water molecules with a spacing distance of 15 Å from complex. Afterward, all added water molecules were minimized using 1000 steps of steepest descent (SD) and continued by 3000 steps of conjugated gradient (CG). Next, the minimizations with the same steps were performed on the whole system. All studied models were heated up from 10 K to 298 K with a constant volume ensemble (NVT) for 100 ps and followed by all-atom MD simulations with a constant pressure ensemble (NPT) at 1 atm and 298 K for 300 ns. All chemical bonds involving hydrogen were constrained using the SHAKE algorithm [63]. The particle mesh Ewald’s method [64] was performed for the treatment of the long-range interactions. The cpptraj module of AMBER16 program was used to calculate the root-mean-square displacement (RMSD), the radius of gyration (Rgyr), and the distance between guest and host molecules as well as the radial distribution function (RDF).

3.3. Quantum Chemical Calculations

All E and K forms of 3HF, 3HFW, and 3HF/γ-CD inclusion complexes in Form I and Form II in the S0 and S1 states were studied by PBE0 and TD-PBE0 methods, respectively with def2-SVP basis set by using Gaussian 16, Revision C.01 [56]. The solvation effect was taken into account by means of the non-equilibrium implementation of the conductor polarized continuum model (C-PCM) framework [65,66], so-called PCM-LR [67]. To confirm that the optimized structures (E, K, E*, and K*) are located at the local minimum, no imaginary frequency from vibrational calculations was found for all optimized structures both in S0 and S1 states. A hydrogen-bonded strength was determined by the important distance parameters involving ESIntraPT and ESInterPT processes namely the covalent O–H bond of 3HF and a water molecule, the intraHB between proton donor and proton acceptor of 3HF, and the interHBs between 3HF and a water molecule for the case study of water assisted effect. The strength of all intraHB and interHBs was further confirmed by the red-shift values of the O–H stretching vibrational frequencies between the S0 and S1 states from the simulated IR spectra, together with the topology analysis at bond critical point (BCP) from quantum theory of atoms in molecules (QTAIM) performed by Multiwfn [68], which was employed in previous studies [69,70,71]. The electronic spectra and frontier MOs were also calculated. Additionally, the absorption and emission spectra of all complexes were simulated to investigate the photophysical properties.

4. Conclusions

The effect of water microsolvation on the strength of hydrogen bonding, ESIntraPT and ESInterPT reactions as well as photophysical behaviors of 3HF and its inclusion complexes in γ-CD has been systematically studied using MD simulations and DFT/TD-DFT at PBE0/TD-PBE0 with def2-SVP basis set. Two possible 3HF/γ-CD inclusion complexes; C-ring and P-ring insertions (Form I and Form II), were observed from molecular docking. From MD results, a lower 3HF mobility in the hydrophobic cavity and a lower water accessibility to the encapsulated 3HF in Form II suggest that this form is favorably more stable. From the static calculation results, the strength of hydrogen bonding of all studied compounds increases upon photoexcitation into the S1 state leading to being easier deprotonation, confirmed by the change of important bond lengths (the increasing of the covalent O–H bond length of proton donor, together with decreasing of the O⋯H intraHBs and interHBs), the red shift of the O–H stretching modes, and the bond energy from the topology analysis. In addition, frontier MOs confirm that the main contribution for vertical S0 to S1 transition is π to π* attributed from HOMO (π) to LUMO (π*). For simulated spectra, the λabs of E, the λemis of E* and the λemis of K* are in good agreement with the experimental data (in the range of 0.52–0.61 eV relative differences) indicating that the present method is adequate to provide the information on their spectra and the possibility of ESPT processes. Besides, the ESIntraPT processes of the inclusion complexes (Form I and Form II) can occur easily like in the case of 3HF in aprotic solvents. Furthermore, the ESInterPT processes via interHBs of Form I-W and Form II-W inclusion complexes are feasible to take place. In addition, K* of Form II/Form II-W is more stable than that of Form I/Form I-W due to the lower energy and the higher oscillator strength. Consequently, the ESIntraPT and ESInterPT might be likely to occur in P-ring insertion in accordance with the MD results, in which P-ring insertion is the majority of the 3HF/γ-CD inclusion complexes with lower water accessibility. However, it is already known that the ESInterPT of 3HF in aqueous solution is hard to occur due to the higher ESInterPT barrier and the higher fluctuation of the water-rearrangement surrounding 3HF, which leads to a decrease of the fluorescence intensity. Thus, from the present work, we found that 3HF is stable inside the γ-CD hydrophobic cavity and promotes ESIntraPT by suppressing the 3HF-water network. This leads to the increment of fluorescent intensity. In the other word, the fluorescence intensity of K* could be efficiently tuned via host-guest complexation.

Supplementary Materials

The following are available online. Figure S1: RMSD plots for all atoms and Rgyr of the two orientations of C-ring insertion (Form I), and P-ring insertion (Form II) inclusion complexes for four different MD runs of each system, Figure S2: The plots of distance measured from the Cm of each 3HF ring to the Cm of the secondary rim of γ-CD (all 7 O2 atoms) for the four MD simulations MD1-MD4 with different initial structures of complexes in Form I and Form II, Figure S3: The simulated absorption spectra of E, and the simulated emission spectra of E* and K* for all studied compounds computed at TD-PBE0/def2-SVP level of theory, Table S1: Electron density ρ(r), the Lagrangian kinetic energy G(r), potential energy density V(r), the Hamiltonian kinetic energy density H(r), the Laplacian of the electron density ∇2ρ(r), the electron delocalization index (DI), and hydrogen-bonded energy (EHB) at selected BCPs in the S1 state (a.u.) for all compounds.

Author Contributions

Conceptualization, N.K., T.R., and S.H.; methodology, K.K.; validation, N.K., T.R., and S.N.; formal analysis, K.K., and C.S.; data curation, K.K.; writing—original draft preparation, K.K.; writing—review and editing, R.D., N.K., T.R., and S.N.; visualization, K.K.; supervision, P.W., N.K., T.R., S.N., and S.H.; project administration, N.K., and S.H.; funding acquisition, N.K., T.R., and S.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by Ratchadapisek Somphot Fund for Postdoctoral Fellowship, Chulalongkorn University (CU).

Acknowledgments

The authors gratefully acknowledge the financial support from the Center of Excellence in Computational Chemistry (CECC) from Chulalongkorn University, the Center of Excellence in Materials Science and Technology, Chiang Mai University and the Office of National Higher Education Science Research and Innovation Policy Council (NXPO) in Global Partnership Project. N.K. thanks the Thailand Research Fund (Grant No. RSA6180044). Computational resources are provided by the Center of Excellence in Computational Chemistry (CECC), Computational Chemistry Laboratory Chiang Mai University (CCL-CMU), and NSTDA Supercomputer center (ThaiSC).

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Zhou, P.; Han, K. Unraveling the detailed mechanism of excited-state proton transfer. Acc. Chem. Res. 2018, 51, 1681–1690. [Google Scholar] [CrossRef]
  2. Uzhinov, B.M.; Khimich, M.N. Conformational effects in excited state intramolecular proton transfer of organic compounds. Russ. Chem. Rev. 2011, 80, 553–577. [Google Scholar] [CrossRef]
  3. Padalkar, V.S.; Seki, S. Excited-state intramolecular proton-transfer (ESIPT)-inspired solid state emitters. Chem. Soc. Rev. 2016, 45, 169–202. [Google Scholar] [CrossRef]
  4. Weller, A. Innermolekularer protonenübergang im angeregten zustand. Z. Elektrochem. Ber. Bunsenges. Phys. Chem. 1956, 60, 1144–1147. [Google Scholar]
  5. Goodman, J.; Brus, L.E. Proton transfer and tautomerism in an excited state of methyl salicylate. J. Am. Chem. Soc. 1978, 100, 7472–7474. [Google Scholar] [CrossRef]
  6. Yoon, M.; Kim, M.; Kim, M.H.; Kang, J.-G.; Sohn, Y.; Kim, I.T. Synthesis and photophysical properties of S, N and Se-modified methyl salicylate derivatives. Inorg. Chim. Acta 2019, 495, 119008. [Google Scholar] [CrossRef]
  7. Lee, J.; Kwon, J.E.; You, Y.; Park, S.Y. Wholly π-Conjugated low-molecular-weight organogelator that displays triple-channel responses to fluoride ions. Langmuir 2014, 30, 2842–2851. [Google Scholar] [CrossRef]
  8. Mitra, S.; Tamai, N. Dynamics of photochromism in salicylideneaniline: A femtosecond spectroscopic study. Phys. Chem. Chem. Phys 2003, 5, 4647–4652. [Google Scholar] [CrossRef] [Green Version]
  9. Borbone, F.; Tuzi, A.; Panunzi, B.; Piotto, S.; Concilio, S.; Shikler, R.; Nabha, S.; Centore, R. On–Off Mechano-responsive Switching of ESIPT Luminescence in Polymorphic N-Salicylidene-4-amino-2-methylbenzotriazole. Cryst. Growth Des. 2017, 17, 5517–5523. [Google Scholar] [CrossRef]
  10. Liu, B.; Wang, J.; Zhang, G.; Bai, R.; Pang, Y. Flavone-based ESIPT ratiometric chemodosimeter for detection of cysteine in living cells. ACS Appl. Mater. Interfaces 2014, 6, 4402–4407. [Google Scholar] [CrossRef] [PubMed]
  11. Chou, P.-T.; Chen, Y.-C.; Yu, W.-S.; Cheng, Y.-M. Spectroscopy and dynamics of excited-state intramolecular proton-transfer reaction in 5-hydroxyflavone. Chem. Phys. Lett. 2001, 340, 89–97. [Google Scholar] [CrossRef]
  12. Yang, Y.; Zhao, J.; Li, Y. Theoretical study of the ESIPT process for a new natural product quercetin. Sci. Rep. 2016, 6, 32152. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, Y.; Deng, Y.; Ji, N.; Zhang, J.; Fan, C.; Ding, T.; Cao, Z.; Li, Y.; Fang, Y. A rationally designed flavone-based ESIPT fluorescent chemodosimeter for highly selective recognition towards fluoride and its application in live-cell imaging. Dyes Pigm. 2019, 166, 473–479. [Google Scholar] [CrossRef]
  14. Henary, M.M.; Fahrni, C.J. Excited state intramolecular proton transfer and metal ion complexation of 2-(2‘-Hydroxyphenyl)benzazoles in aqueous solution. J. Phys. Chem. A 2002, 106, 5210–5220. [Google Scholar] [CrossRef]
  15. Das, S.; Chattopadhyay, N. Heteroatom controlled probe-water cluster formation of a series of ESIPT probes: An exploration with fluorescence anisotropy. Chem. Phys. Lett. 2018, 708, 37–41. [Google Scholar] [CrossRef]
  16. Liu, C.; Wang, F.; Xiao, T.; Chi, B.; Wu, Y.; Zhu, D.; Chen, X. The ESIPT fluorescent probes for N2H4 based on benzothiazol and their applications for gas sensing and bioimaging. Sens. Actuators B 2018, 256, 55–62. [Google Scholar] [CrossRef]
  17. Mai, V.T.N.; Shukla, A.; Mamada, M.; Maedera, S.; Shaw, P.E.; Sobus, J.; Allison, I.; Adachi, C.; Namdas, E.B.; Lo, S.-C. Low amplified spontaneous emission threshold and efficient electroluminescence from a carbazole derivatized excited-state intramolecular proton transfer dye. ACS Photonics 2018, 5, 4447–4455. [Google Scholar] [CrossRef]
  18. Munch, M.; Curtil, M.; Vérité, P.M.; Jacquemin, D.; Massue, J.; Ulrich, G. Ethynyl-Tolyl Extended 2-(2′-Hydroxyphenyl)benzoxazole Dyes: Solution and Solid-state Excited-State Intramolecular Proton Transfer (ESIPT) Emitters. Eur. J. Org. Chem. 2019, 2019, 1134–1144. [Google Scholar] [CrossRef]
  19. Song, Z.; Kwok, R.T.K.; Zhao, E.; He, Z.; Hong, Y.; Lam, J.W.Y.; Liu, B.; Tang, B.Z. A ratiometric fluorescent probe based on ESIPT and AIE processes for alkaline phosphatase activity assay and visualization in living cells. ACS Appl. Mater. Interfaces 2014, 6, 17245–17254. [Google Scholar] [CrossRef]
  20. Dommett, M.; Crespo-Otero, R. Excited state proton transfer in 2′-hydroxychalcone derivatives. Phys. Chem. Chem. Phys. 2017, 19, 2409–2416. [Google Scholar] [CrossRef]
  21. Luo, Z.; Liu, B.; Qin, T.; Zhu, K.; Zhao, C.; Pan, C.; Wang, L. Cyclization of chalcone enables ratiometric fluorescence determination of hydrazine with a high selectivity. Sens. Actuators B 2018, 263, 229–236. [Google Scholar] [CrossRef]
  22. Kwon, J.E.; Park, S.Y. Advanced organic optoelectronic materials: Harnessing excited-state intramolecular proton transfer (ESIPT) process. Adv. Mater. 2011, 23, 3615–3642. [Google Scholar] [CrossRef]
  23. Ashton, T.D.; Jolliffe, K.A.; Pfeffer, F.M. Luminescent probes for the bioimaging of small anionic species in vitro and in vivo. Chem. Soc. Rev. 2015, 44, 4547–4595. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lee, M.H.; Kim, J.S.; Sessler, J.L. Small molecule-based ratiometric fluorescence probes for cations, anions, and biomolecules. Chem. Soc. Rev. 2015, 44, 4185–4191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Azarias, C.; Budzák, Š.; Laurent, A.D.; Ulrich, G.; Jacquemin, D. Tuning ESIPT fluorophores into dual emitters. Chem. Sci. 2016, 7, 3763–3774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sedgwick, A.C.; Wu, L.; Han, H.-H.; Bull, S.D.; He, X.-P.; James, T.D.; Sessler, J.L.; Tang, B.Z.; Tian, H.; Yoon, J. Excited-state intramolecular proton-transfer (ESIPT) based fluorescence sensors and imaging agents. Chem. Soc. Rev. 2018, 47, 8842–8880. [Google Scholar] [CrossRef] [Green Version]
  27. Massue, J.; Pariat, T.M.; Vérité, P.; Jacquemin, D.; Durko, M.; Chtouki, T.; Sznitko, L.; Mysliwiec, J.; Ulrich, G. Natural born laser dyes: Excited-State Intramolecular Proton Transfer (ESIPT) Emitters and their use in random lasing studies. Nanomaterials 2019, 9, 1093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Gunduz, S.; Goren, A.C.; Ozturk, T. Facile Syntheses of 3-Hydroxyflavones. Org. Lett. 2012, 14, 1576–1579. [Google Scholar] [CrossRef] [PubMed]
  29. Sengupta, P.K.; Kasha, M. Excited state proton-transfer spectroscopy of 3-hydroxyflavone and quercetin. Chem. Phys. Lett. 1979, 68, 382–385. [Google Scholar] [CrossRef]
  30. Kasha, M. Proton-transfer spectroscopy. Perturbation of the tautomerization potential. J. Chem. Soc. Faraday Trans. 2 1986, 82, 2379–2392. [Google Scholar] [CrossRef]
  31. Ash, S.; De, S.P.; Pyne, S.; Misra, A. Excited state intramolecular proton transfer in 3-hydroxy flavone and 5-hydroxy flavone: A DFT based comparative study. J. Mol. Model. 2010, 16, 831–839. [Google Scholar] [CrossRef] [PubMed]
  32. Das, S.; Chakrabarty, S.; Chattopadhyay, N. Origin of unusually high fluorescence anisotropy of 3-hydroxyflavone in water: Formation of probe–solvent cage-like cluster. J. Phys. Chem. B 2020, 124, 173–180. [Google Scholar] [CrossRef]
  33. Lazzaroni, S.; Dondi, D.; Mezzetti, A.; Protti, S. Role of solute-solvent hydrogen bonds on the ground state and the excited state proton transfer in 3-hydroxyflavone. A systematic spectrophotometric study. Photochem. Photobiol. Sci. 2018, 17, 923–933. [Google Scholar] [CrossRef] [PubMed]
  34. Klymchenko, A.S.; Demchenko, A.P. Multiparametric probing of intermolecular interactions with fluorescent dye exhibiting excited state intramolecular proton transfer. Phys. Chem. Chem. Phys 2003, 5, 461–468. [Google Scholar] [CrossRef]
  35. Klymchenko, A.S.; Pivovarenko, V.G.; Ozturk, T.; Demchenko, A.P. Modulation of the solvent-dependent dual emission in 3-hydroxychromones by substituents. New J. Chem. 2003, 27, 1336–1343. [Google Scholar] [CrossRef]
  36. Klymchenko, A.S.; Kenfack, C.; Duportail, G.; Mély, Y. Effects of polar protic solvents on dual emissions of 3-hydroxychromones. J. Chem. Sci. 2007, 119, 83–89. [Google Scholar] [CrossRef]
  37. Klymchenko, A.S.; Demchenko, A.P. Chapter 3 multiparametric probing of microenvironment with solvatochromic fluorescent dyes. Meth. Enzymol. 2008, 450, 37–58. [Google Scholar]
  38. Ghosh, D.; Pradhan, A.K.; Mondal, S.; Begum, N.A.; Mandal, D. Proton transfer reactions of 4′-chloro substituted 3-hydroxyflavone in solvents and aqueous micelle solutions. Phys. Chem. Chem. Phys 2014, 16, 8594–8607. [Google Scholar] [CrossRef] [PubMed]
  39. Ghosh, D.; Batuta, S.; Das, S.; Begum, N.A.; Mandal, D. Proton transfer dynamics of 4′-N,N-dimethylamino-3-hydroxyflavone observed in hydrogen-bonding solvents and aqueous micelles. J. Phys. Chem. B 2015, 119, 5650–5661. [Google Scholar] [CrossRef]
  40. Salaeh, R.; Prommin, C.; Chansen, W.; Kerdpol, K.; Daengngern, R.; Kungwan, N. The effect of protic solvents on the excited state proton transfer of 3-hydroxyflavone: A TD-DFT static and molecular dynamics study. J. Mol. Liq. 2018, 252, 428–438. [Google Scholar] [CrossRef]
  41. Zhao, J.; Ji, S.; Chen, Y.; Guo, H.; Yang, P. Excited state intramolecular proton transfer (ESIPT): From principal photophysics to the development of new chromophores and applications in fluorescent molecular probes and luminescent materials. Phys. Chem. Chem. Phys. 2012, 14, 8803–8817. [Google Scholar] [CrossRef]
  42. Santos, F.S.; Ramasamy, E.; Ramamurthy, V.; Rodembusch, F.S. Excited state chemistry of flavone derivatives in a confined medium: ESIPT emission in aqueous media. Photochem. Photobiol. Sci. 2014, 13, 992–996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Bartl, K.; Funk, A.; Gerhards, M. IR/UV spectroscopy on jet cooled 3-hydroxyflavone (H2O)n (n = 1,2) clusters along proton transfer coordinates in the electronic ground and excited states. J. Chem. Phys. 2008, 129, 234306. [Google Scholar] [CrossRef] [PubMed]
  44. Stamm, A.; Maué, D.; Gerhards, M. Structural rearrangement by isomer-specific infrared excitation in the neutral isolated dihydrated cluster of 3-hydroxyflavone. J. Phys. Chem. Lett. 2018, 9, 4360–4366. [Google Scholar] [CrossRef]
  45. Banerjee, A.; Sengupta, P.K. Encapsulation of 3-hydroxyflavone and fisetin in β-cyclodextrins: Excited state proton transfer fluorescence and molecular mechanics studies. Chem. Phys. Lett. 2006, 424, 379–386. [Google Scholar] [CrossRef]
  46. Das, S.; Chattopadhyay, N. Supramolecular inclusion-assisted disruption of probe-solvent network. ChemistrySelect 2017, 2, 6078–6081. [Google Scholar] [CrossRef]
  47. Pahari, B.; Chakraborty, S.; Sengupta, P.K. Encapsulation of 3-hydroxyflavone in γ-cyclodextrin nanocavities: Excited state proton transfer fluorescence and molecular docking studies. J. Mol. Struct. 2011, 1006, 483–488. [Google Scholar] [CrossRef]
  48. Szejtli, J. Introduction and general overview of cyclodextrin chemistry. Chem. Rev. 1998, 98, 1743–1754. [Google Scholar] [CrossRef]
  49. Arunkumar, E.; Forbes, C.C.; Smith, B.D. Improving the properties of organic dyes by molecular encapsulation. Eur. J. Org. Chem. 2005, 2005, 4031. [Google Scholar] [CrossRef]
  50. Hou, X.; Ke, C.; Bruns, C.J.; McGonigal, P.R.; Pettman, R.B.; Stoddart, J.F. Tunable solid-state fluorescent materials for supramolecular encryption. Nat. Commun. 2015, 6, 6884. [Google Scholar] [CrossRef] [Green Version]
  51. Pahari, B.; Sengupta, B.; Chakraborty, S.; Thomas, B.; McGowan, D.; Sengupta, P.K. Contrasting binding of fisetin and daidzein in γ-cyclodextrin nanocavity. J. Photochem. Photobiol. B 2013, 118, 33–41. [Google Scholar] [CrossRef] [Green Version]
  52. Douhal, A. Ultrafast guest dynamics in cyclodextrin nanocavities. Chem. Rev. 2004, 104, 1955–1976. [Google Scholar] [CrossRef]
  53. Kicuntod, J.; Khuntawee, W.; Wolschann, P.; Pongsawasdi, P.; Chavasiri, W.; Kungwan, N.; Rungrotmongkol, T. Inclusion complexation of pinostrobin with various cyclodextrin derivatives. J. Mol. Graph. Model. 2016, 63, 91–98. [Google Scholar] [CrossRef] [PubMed]
  54. Nutho, B.; Khuntawee, W.; Rungnim, C.; Pongsawasdi, P.; Wolschann, P.; Karpfen, A.; Kungwan, N.; Rungrotmongkol, T. Binding mode and free energy prediction of fisetin/β-cyclodextrin inclusion complexes. Beilstein J. Org. Chem. 2014, 10, 2789–2799. [Google Scholar] [CrossRef] [Green Version]
  55. Dennington, R.; Keith, T.A.; Millam, J.M. (Eds.) Gauss View; Version 6; Semichem Inc.: Shawnee Mission, KS, USA, 2016. [Google Scholar]
  56. Gaussian 16, Revision, C.01; Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A. (Eds.) Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  57. Case, D.A.; Betz, R.M.; Cerutti, D.S.; Cheatham, T.E., III; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; Homeyer, N.; et al. AMBER; University of California: San Francisco, CA, USA, 2016. [Google Scholar]
  58. Kaiyawet, N.; Rungrotmongkol, T.; Hannongbua, S. Effect of Halogen substitutions on dUMP to stability of thymidylate synthase/dUMP/mTHF ternary complex using molecular dynamics simulation. J. Chem. Inf. Model. 2013, 53, 1315–1323. [Google Scholar] [CrossRef]
  59. Mahalapbutr, P.; Chusuth, P.; Kungwan, N.; Chavasiri, W.; Wolschann, P.; Rungrotmongkol, T. Molecular recognition of naphthoquinone-containing compounds against human DNA topoisomerase IIα ATPase domain: A molecular modeling study. J. Mol. Liq. 2017, 247, 374–385. [Google Scholar] [CrossRef]
  60. Mahalapbutr, P.; Wonganan, P.; Charoenwongpaiboon, T.; Prousoontorn, M.; Chavasiri, W.; Rungrotmongkol, T. Enhanced solubility and anticancer potential of mansonone G By β-cyclodextrin-based host-guest complexation: A computational and experimental study. Biomolecules 2019, 9, 545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a general amber force field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef] [PubMed]
  62. Kirschner, K.N.; Yongye, A.B.; Tschampel, S.M.; González-Outeiriño, J.; Daniels, C.R.; Foley, B.L.; Woods, R.J. GLYCAM06: A generalizable biomolecular force field. Carbohydrates. J. Comput. Chem. 2008, 29, 622–655. [Google Scholar] [CrossRef] [Green Version]
  63. Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H.J.C. Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 1977, 23, 327–341. [Google Scholar] [CrossRef] [Green Version]
  64. Luty, B.A.; van Gunsteren, W.F. Calculating electrostatic interactions using the particle−particle particle−mesh method with nonperiodic long-range interactions. J. Phys. Chem. 1996, 100, 2581–2587. [Google Scholar] [CrossRef]
  65. Jacquemin, D.; Perpète, E.A.; Scalmani, G.; Frisch, M.J.; Assfeld, X.; Ciofini, I.; Adamo, C. Time-dependent density functional theory investigation of the absorption, fluorescence, and phosphorescence spectra of solvated coumarins. J. Chem. Phys. 2006, 125, 164324. [Google Scholar] [CrossRef] [PubMed]
  66. Corni, S.; Cammi, R.; Mennucci, B.; Tomasi, J. Electronic excitation energies of molecules in solution within continuum solvation models: Investigating the discrepancy between state-specific and linear-response methods. J. Chem. Phys. 2005, 123, 134512. [Google Scholar] [CrossRef]
  67. Jacquemin, D.; Mennucci, B.; Adamo, C. Excited-state calculations with TD-DFT: From benchmarks to simulations in complex environments. Phys. Chem. Chem. Phys 2011, 13, 16987–16998. [Google Scholar] [CrossRef]
  68. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef] [PubMed]
  69. Li, Y.; Wen, K.; Feng, S.; Yuan, H.; An, B.; Zhu, Q.; Guo, X.; Zhang, J. Tunable excited-state intramolecular proton transfer reactions with NH or OH as a proton donor: A theoretical investigation. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2017, 187, 9–14. [Google Scholar] [CrossRef] [PubMed]
  70. Prommin, C.; Kerdpol, K.; Saelee, T.; Kungwan, N. Effects of π-expansion, an additional hydroxyl group, and substitution on the excited state single and double proton transfer of 2-hydroxybenzaldehyde and its relative compounds: TD-DFT static and dynamic study. New J. Chem. 2019, 43, 19107–19119. [Google Scholar] [CrossRef]
  71. Sun, C.; Su, X.; Zhou, Q.; Shi, Y. Regular tuning of the ESIPT reaction of 3-hydroxychromone-based derivatives by substitution of functional groups. Org. Chem. Front. 2019, 6, 3093–3100. [Google Scholar] [CrossRef]
Figure 1. (A) Chemical structure of 3-Hydroxyflavone, 3HF. (B) Docked structures of the two possible 3HF/γ-CD complexes, Form I and Form II, where their percentage of occurrence and the lowest interaction energy retrieved from 100 independent docking runs are also given.
Figure 1. (A) Chemical structure of 3-Hydroxyflavone, 3HF. (B) Docked structures of the two possible 3HF/γ-CD complexes, Form I and Form II, where their percentage of occurrence and the lowest interaction energy retrieved from 100 independent docking runs are also given.
Molecules 26 00843 g001
Figure 2. The plots of distance measured from the Cm of each 3HF ring to the Cm of the primary rim of γ-CD (all 7 O6 atoms) for the four MD simulations MD1-MD4 with different initial structures of complexes in Form I and Form II.
Figure 2. The plots of distance measured from the Cm of each 3HF ring to the Cm of the primary rim of γ-CD (all 7 O6 atoms) for the four MD simulations MD1-MD4 with different initial structures of complexes in Form I and Form II.
Molecules 26 00843 g002
Figure 3. (A) Radial distribution function (RDF or g(r)) of water oxygen atoms around the O1-O3 atoms of 3HF encapsulated in the γ-CD cavity over the last 50-ns MD simulations for Form I (MD3-MD4) and Form II (MD1-MD4 with initial structures in this form and additional MD1-MD2 from re-encapsulation process found in Form I). (B) Average integration number, n(r), up to the first minimum derived from RDF plots, corresponding to the number of water molecules pointing toward the focused oxygens of 3HF.
Figure 3. (A) Radial distribution function (RDF or g(r)) of water oxygen atoms around the O1-O3 atoms of 3HF encapsulated in the γ-CD cavity over the last 50-ns MD simulations for Form I (MD3-MD4) and Form II (MD1-MD4 with initial structures in this form and additional MD1-MD2 from re-encapsulation process found in Form I). (B) Average integration number, n(r), up to the first minimum derived from RDF plots, corresponding to the number of water molecules pointing toward the focused oxygens of 3HF.
Molecules 26 00843 g003
Figure 4. All S0 optimized structures of 3HF, 3HFW, and the different conformations of the 3HF/γ-CD inclusion complexes (Form I and Form II) as well as its inclusion complex with a water molecule (Form I-W and Form II-W) computed by PBE0/def2-SVP level of theory. The blue and green dot lines represent intraHB in 3HF and interHBs between 3HF and a water molecule, respectively.
Figure 4. All S0 optimized structures of 3HF, 3HFW, and the different conformations of the 3HF/γ-CD inclusion complexes (Form I and Form II) as well as its inclusion complex with a water molecule (Form I-W and Form II-W) computed by PBE0/def2-SVP level of theory. The blue and green dot lines represent intraHB in 3HF and interHBs between 3HF and a water molecule, respectively.
Molecules 26 00843 g004
Figure 5. Frontier MOs of all studied compounds.
Figure 5. Frontier MOs of all studied compounds.
Molecules 26 00843 g005
Table 1. A summary of the important bonds and distances involving PT process of enol form for all complexes in the S0 and S1 states.
Table 1. A summary of the important bonds and distances involving PT process of enol form for all complexes in the S0 and S1 states.
CompoundImportant Bond Distance (Å)
S0 StateS1 State
O1–H1O2⋯H1Ow⋯H1Ow–HwO2⋯HwO1–H1O2⋯H1Ow⋯H1Ow–HwO2⋯Hw
3HF0.9831.920 1.0171.706
3HFW1.0002.3871.6300.9841.7231.0692.4391.3941.0081.572
Form I0.9791.978 1.0021.804
Form II0.9812.177 0.9962.080
Form I-W1.0092.3901.5850.9771.8081.0302.4151.5030.9781.903
Form II-W1.0042.4291.6120.9831.7331.0782.4521.3751.0061.581
Table 2. A summary of the values of the O1–H1 and the Ow–Hw stretching vibrational modes of enol form for all compounds in both S0 and S1 states and their spectral shifts ( Δ υ in cm−1).
Table 2. A summary of the values of the O1–H1 and the Ow–Hw stretching vibrational modes of enol form for all compounds in both S0 and S1 states and their spectral shifts ( Δ υ in cm−1).
CompoundWavenumber (cm−1)
O1–H1Ow–Hw
S0S1 Δ υ S0S1 Δ υ
3HF35443002541
3HFW3129259953035103074436
Form I36173250367
Form II35743330244
Form I-W2967257539236623499163
Form II-W3063257748635263111415
Table 3. UV/Vis absorption band maxima of enol form (λabs of E), emission band maxima of enol (λemis of E) and emission band maxima of keto forms (λemis of K), the excitation energy (Eex), and the oscillator strength (f) as well as their major contribution (%) calculated by TD-PBE0/def2-SVP level of theory.
Table 3. UV/Vis absorption band maxima of enol form (λabs of E), emission band maxima of enol (λemis of E) and emission band maxima of keto forms (λemis of K), the excitation energy (Eex), and the oscillator strength (f) as well as their major contribution (%) calculated by TD-PBE0/def2-SVP level of theory.
CompoundsUV/Vis AbsorptionEmission
λabs of E
(nm)
Eex
(eV)
fMOs (%contribution)λemis of E*
(nm)
Eex
(eV)
fλemis of K*
(nm)
Eex
(eV)
f
3HF3413.630.4907HOMO→LUMO (98%)3903.180.54545142.410.4632
3HFW3513.530.4857HOMO→LUMO (98%)408 3.040.53355042.460.4619
Form I3453.590.4618HOMO→LUMO (98%)3923.160.51235242.360.3692
Form II3453.590.4081 HOMO→LUMO (98%)3933.150.46895282.350.3792
Form I-W3533.510.4288HOMO→LUMO (98%)3993.110.52135292.340.3708
Form II-W3503.540.4174HOMO→LUMO (98%)4123.010.45755082.440.3899
Table 4. The relative energy and computed energy differences between the E and K forms (ΔE = Eenol − Eketo) in the S0 and S1 states for all complexes.
Table 4. The relative energy and computed energy differences between the E and K forms (ΔE = Eenol − Eketo) in the S0 and S1 states for all complexes.
ComplexRelative Energy (kcal/mol)
S0S1ΔE at S0ΔE at S1
EKE*K*
3HF010.7475.1965.8110.74−9.38
3HFW010.2774.9263.4110.27−11.51
Form I016.8574.7971.1616.85−3.63
Form II08.4772.3964.688.47−7.71
Form I-W06.6371.2463.436.63−7.81
Form II-W08.3573.1264.728.35−8.40
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kerdpol, K.; Daengngern, R.; Sattayanon, C.; Namuangruk, S.; Rungrotmongkol, T.; Wolschann, P.; Kungwan, N.; Hannongbua, S. Effect of Water Microsolvation on the Excited-State Proton Transfer of 3-Hydroxyflavone Enclosed in γ-Cyclodextrin. Molecules 2021, 26, 843. https://doi.org/10.3390/molecules26040843

AMA Style

Kerdpol K, Daengngern R, Sattayanon C, Namuangruk S, Rungrotmongkol T, Wolschann P, Kungwan N, Hannongbua S. Effect of Water Microsolvation on the Excited-State Proton Transfer of 3-Hydroxyflavone Enclosed in γ-Cyclodextrin. Molecules. 2021; 26(4):843. https://doi.org/10.3390/molecules26040843

Chicago/Turabian Style

Kerdpol, Khanittha, Rathawat Daengngern, Chanchai Sattayanon, Supawadee Namuangruk, Thanyada Rungrotmongkol, Peter Wolschann, Nawee Kungwan, and Supot Hannongbua. 2021. "Effect of Water Microsolvation on the Excited-State Proton Transfer of 3-Hydroxyflavone Enclosed in γ-Cyclodextrin" Molecules 26, no. 4: 843. https://doi.org/10.3390/molecules26040843

APA Style

Kerdpol, K., Daengngern, R., Sattayanon, C., Namuangruk, S., Rungrotmongkol, T., Wolschann, P., Kungwan, N., & Hannongbua, S. (2021). Effect of Water Microsolvation on the Excited-State Proton Transfer of 3-Hydroxyflavone Enclosed in γ-Cyclodextrin. Molecules, 26(4), 843. https://doi.org/10.3390/molecules26040843

Article Metrics

Back to TopTop