Next Article in Journal
Co-Formulation of Amphiphilic Cationic and Anionic Cyclodextrins Forming Nanoparticles for siRNA Delivery in the Treatment of Acute Myeloid Leukaemia
Next Article in Special Issue
Ultrafast Continuum IR Generation and Its Application in IR Spectroscopy
Previous Article in Journal
Bioenergetic and Metabolic Impairments in Induced Pluripotent Stem Cell-Derived Cardiomyocytes Generated from Duchenne Muscular Dystrophy Patients
Previous Article in Special Issue
Mechanism of Ir(ppy)3 Guest Exciton Formation with the Exciplex-Forming TCTA:TPBI Cohost within a Phosphorescent Organic Light-Emitting Diode Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermodynamics of π–π Interactions of Benzene and Phenol in Water

1
Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea
2
Department of Chemistry, Pusan National University, Busan 46241, Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(17), 9811; https://doi.org/10.3390/ijms23179811
Submission received: 25 July 2022 / Revised: 20 August 2022 / Accepted: 22 August 2022 / Published: 29 August 2022
(This article belongs to the Special Issue State-of-the-Art Physical Chemistry and Chemical Physics in Korea)

Abstract

:
The π–π interaction is a major driving force that stabilizes protein assemblies during protein folding. Recent studies have additionally demonstrated its involvement in the liquid–liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs). As the participating residues in IDPs are exposed to water, π–π interactions for LLPS must be modeled in water, as opposed to the interactions that are often established at the hydrophobic domains of folded proteins. Thus, we investigated the association of free energies of benzene and phenol dimers in water by integrating van der Waals (vdW)-corrected density functional theory (DFT) and DFT in classical explicit solvents (DFT-CES). By comparing the vdW-corrected DFT and DFT-CES results with high-level wavefunction calculations and experimental solvation free energies, respectively, we established the quantitative credibility of these approaches, enabling a reliable prediction of the benzene and phenol dimer association free energies in water. We discovered that solvation influences dimer association free energies, but not significantly when no direct hydrogen-bond-type interaction exists between two monomeric units, which can be explained by the enthalpy–entropy compensation. Our comprehensive computational study of the solvation effect on π–π interactions in water could help us understand the molecular-level driving mechanism underlying the IDP phase behaviors.

1. Introduction

The π–π interaction, a typical noncovalent bond involving π-electronic aromatic systems, is one of the prevalent molecular interactions in nature [1]. The association of large biomolecules such as proteins [2,3] and nucleic acids [4] typically involves π–π interaction. Additionally, it plays a significant role in biological processes such as charge transfer [5,6,7], molecular recognition [8,9], and drug delivery [10].
Since Burley identified an interaction between aromatic rings that stabilizes protein structure in 1985 [2], there have been extensive attempts in the realm of quantum chemistry to elucidate the nature of the π–π interaction. The Schlag group obtained benzene dimer interaction energies using Møller–Plesset perturbation theories (MP2 and MP4) and coupled-cluster single-double and perturbative triple (CCSD(T)) theory [11,12,13,14]. They identified three stable geometries of benzene dimer, namely sandwich (S), parallel displaced (PD), and T-shape (T) structures, and discovered that PD and T are the most stable structures with similar binding energies. The Sherrill group further investigated the physical nature and substituent effect of the π–π interaction and provided benchmark-quality potential energy curves for a benzene dimer calculated at the CCSD (T) level with a complete basis set (CBS) limit [15,16,17,18]. These studies provide a fundamental understanding of the π–π interaction and highlight the significance of its role in stabilizing the structure of folded proteins.
Intrinsically disordered proteins (IDPs) are a subclass of proteins lacking cooperatively folded structures in their native state [19], and their biological significance has been extensively studied in recent years. IDPs may undergo liquid–liquid phase separation (LLPS), which is sequence-dependent, and the significance of the π–π interaction in determining their LLPS behavior has been highlighted [20,21]. For example, Gabryelczyk et al. discovered that interactions between tyrosine residues initiate the LLPS of histidine-rich proteins, forming microdroplets [22]. The Pappu group modeled the inter-residue interactions using the stickers-and-spacers framework to predict the phase behavior of IDPs and discovered that the uniform patterning of aromatic residues promotes LLPS and inhibits aggregation [23,24,25]. At present, however, many molecular features of inter-residue interactions, governing the LLPS of IDPs, remain obscure.
While the residues of folded proteins develop the π–π interaction primarily in the hydrophobic pocket, which could be modeled in vacuo, as was likely done in numerous previous quantum chemical studies, the residues of IDPs frequently protrude from the protein backbone toward solvents to develop an interchain π–π interaction. Assuming that the solvation effect’s presence is a major distinction between the inter-residue interactions of folded proteins and IDPs, we explored aromatic π–π interactions in aqueous environments. We examined how water solvation modifies the π–π interaction of benzene and phenol dimers by combining two modern computational chemistry methods: (1) the van der Waals (vdW)-corrected density functional theory (DFT) for an accurate description of direct monomer–monomer interaction energy; (2) DFT in classical explicit solvents (DFT-CES), which is a mean-field quantum mechanics/molecular mechanics (QM/MM) method [26,27,28] enabling an explicit treatment of solvent molecules for a reliable description of solvation free energy. We found that the solvation effect renormalizes the interaction energies between benzene and phenol molecules, the conformational dependence of which can be understood in terms of the changes in the number of hydrogen bonds (HB) before and after dimer formation: water–benzene HB for benzene dimers and water–hydroxyl HB for phenol dimers. The lack of a significant energy difference between in vacuo and in aqua is ascribed to the enthalpy–entropy compensation, which preserves the most stable configuration of a benzene dimer as a T-shape and a phenol dimer as a hydrogen-bonded pair between hydroxyl groups. We expect that our study, which quantifies the interaction strength between aromatic π systems in water, will provide a molecular understanding of the π–π interaction mediated phase behavior of IDPs.

2. Methods

2.1. Brief Review of DFT-CES/2PT

DFT-CES is a grid-based mean-field QM/MM method that was recently developed by the Kim group [26,27,28]. The mean-field QM/MM method iteratively solves QM optimizations and molecular dynamics (MD) simulations until a self-consistent solution is obtained [29,30,31]. Each QM optimization step relaxes the electron density ( ρ QM ) and structure ( r QM ) of the QM particles in the presence of the ensemble-averaged electrostatic potential of the MM component ( V MM ), whereas each MD simulation propagates the positions ( r MM ) and momenta ( p MM ) of MM particles in the presence of frozen QM particles and ρ QM . Throughout this iterative procedure, known as the self-consistent ensemble-averaged reaction field (SCERF), the free energy functional of the total QM/MM system, A tot , is minimized for the QM nuclear and electronic degrees of freedom [32]:
A tot = E KS [ ρ QM ; r QM ] k B T e β H MM d r MM d p MM
Here, the first term indicates the QM subsystem’s internal energy ( E QM int ), where E KS is the Kohn–Sham DFT total-energy functional. The second term represents the Helmholtz free energy of the MM subsystem ( A MM ), where H MM is the Hamiltonian regarding the dynamics of MM particles subject to the external potential applied by the frozen QM subsystem (and, thereby, it includes the interaction between the QM and MM regions). E QM int can be directly evaluated by DFT calculations, whereas free energy A MM cannot be directly computed by conventional MD simulations, since such a calculation is often accomplished by turning on and off the constant solute potential [27].
To determine A MM , we utilized the two-phase thermodynamics (2PT) model established by Lin et al. [33], which permits a direct computation of free energy quantities from MD trajectories. The 2PT model separates the total degrees of freedom in a liquid system into gaseous and solid components. The total free energy of the system can then be calculated from the linear combination of the gaseous and solid components, each of which has been studied theoretically using the hard-sphere theory and quasi-harmonic oscillator model, respectively. Notably, the 2PT model has been effective in estimating the thermodynamics and phase behavior of Lennard–Jones particles at various densities [33], the absolute entropy of water molecules [34] and organic solvents [35], and the surface tension of typical liquids [36]. Furthermore, the 2PT model permits the direct decomposition of free energy into enthalpic and entropic components, hence, enhancing the comprehension of the investigated thermodynamic process.
By treating the solutes with QM and the solvents with MM, the DFT-CES/2PT method can be used to evaluate the solvation free energies Δ G sol . Considering that the solvation process incurs negligible pressure–volume work, Δ G sol can be evaluated by computing the A tot . Our previous DFT-CES experiments paired with the 2PT method demonstrated high accuracy in predicting the Δ G sol of polar [26] and nonpolar molecules [27] and the thermodynamics and structure of solid–liquid interfaces [37,38,39,40].

2.2. Computational Details

We implemented DFT-CES by coupling two open-source programs: the plane-wave (PW) DFT code Quantum Espresso [41] and the large-scale atomic/molecular massively parallel simulator (LAMMPS) [42]. For the DFT calculation, we employed the Perdew–Burke–Ernzerhof (PBE) exchange-correlation functional [43] with Grimme’s D3-type vdW correction [44] and Becke–Johnson damping [45]. The electron–nucleus interaction was described using the projector augmented-wave (PAW) method [46,47] with a 50 Ry kinetic energy cutoff.
For the MD simulation at each SCERF iteration, we performed a canonical ensemble (NVT) MD simulation for 2.5 ns at 300 K using the Nosé–Hoover thermostat [48,49] with a 100-fs damping constant. The final 2 ns of the MD trajectory were used to evaluate the V MM . The intermolecular interactions of water molecules were represented using a modified TIP3P water potential [50], with bond lengths and angles constrained to equilibrium values using the RATTLE algorithm [51]. Long-range Coulomb interactions were treated using the Ewald summation method [52] with a real-space cut-off of 15 Å.
The simulation box, consisting of one or two QM solute molecules (either benzene or phenol) solvated by 1000 water molecules, has a cubic shape with a box length of 31.2 Å. The vdW pair interaction between QM and MM particles was described using the OPLS-AA force field [53], and SCERF iterations were conducted until the difference in E QM int was less than 0.1 kcal/mol. DFT structural optimization may place the QM system too close to the site of averaged MM particles because the DFT-CES involves Pauli repulsion between the QM and MM regions only in the MD component utilizing the pairwise potential, but not in the DFT part. To prevent such a technical issue with phenol systems, the H atom from the hydroxyl group in phenol was maintained at its original position in the structure optimized in vacuo, while no constraints were imposed on the benzene system.
We then calculated the free energy of association in water, Δ G aq , using the thermodynamic cycle shown in Figure 1. The free energy of association in vacuo, Δ G vac , was calculated using vdW-corrected DFT energetics, wherein the vibrational free energies, including the zero-point energy (ZPE) contribution, were corrected using harmonic oscillator partition functions. We estimated the vibrational frequencies for computational efficiency using the Gaussian orbital code, Jaguar version 10.8 [54], at the level of PBE-D3/AUG-cc-pVTZ-PP(-F). Then, using separate DFT-CES/2PT calculations, the solvation free energies of the monomers and their dimer complexes, i.e., Δ G mon sol and Δ G dim sol , were calculated, leading to Δ G aq = Δ G vac +   Δ Δ G sol , where Δ Δ G sol = Δ G dim sol 2 Δ G mon sol .

2.3. Exploration of Different Conformations

Herein, we focused on the π–π interactions in benzene and phenol dimers. As in earlier quantum chemical investigations [11,12,13,14,15,16,17,18], stacked (ST), parallel-displaced (PD), and T-shape (T) conformations were explored for benzene dimers (Figure 2a). It is to be noted that two T conformations with different symmetries (labeled as T1 and T2) were investigated for completeness.
Phenol dimers can assume an extensive variety of different conformations depending on the relative locations of the hydroxyl (OH) groups attached to the phenyl rings. To systematically investigate the impact of the relative locations of OH groups, we first constructed two distinct ST conformations of the phenol dimer: one by positioning OH groups with the closest distance, and the other by positioning OH groups with the farthest distance, which are labeled ST (1,1′) and ST (1,4′), respectively (see Figure 2b). By horizontally shifting one of the phenyl rings of ST (1,1′) and ST (1,4′), we constructed three distinct PD conformations of phenol dimers, labeled PD (1,1′), PD (4,1′), and PD (1,4′), as shown in Figure 2b.
For the T conformations, the OH group location in the horizontally laid phenyl ring (whose carbons are labeled using primed numbers in Figure 2a) is thought to be less important than the OH group location in the vertically standing phenyl ring (whose carbons are labeled using unprimed numbers in Figure 2a). Indeed, we discovered virtually no difference in Δ G vac with different OH group locations in the horizontally laid phenyl ring (Figure S1). Therefore, we chose two representative cases, T1 (1,2′) and T2 (1,1′), as shown in Figure 2b. To explore the effect of the OH group location on the vertically standing phenyl ring, we further chose the T2 (2,1′), T2 (3,1′), and T2 (4,1′) conformations.
Finally, as a hydrogen-bonded pair (HBP) can develop between the hydroxyl groups of two phenol molecules, the corresponding geometry was also added for further investigation, although it cannot be strictly classified as a π–π interaction.

3. Results and Discussion

3.1. Benchmark: Dimer Binding Energy in Vacuo & Monomer Solvation Energy

Table 1 shows the binding energies of the benzene and phenol dimers in vacuo ( Δ E vac ), which were calculated using the difference in the DFT self-consistent field (SCF) energies between the dimer and monomer, Δ E vac = E dimer SCF 2 E monomer SCF . Incorporating the entropic corrections and the ZPE, the calculation of the corresponding free energy, Δ G vac , is allowed as listed in Table 1. To evaluate the accuracy of PBE-D3, we compared the Δ E vac values to the results of previous studies using high-level wavefunction theory, such as coupled cluster theory. For the S, PD, and T1 conformations of the benzene dimer, the results from the PBE-D3/PW level and the estimated CCSD(T)/aug-cc-pVQZ* level agreed within approximately 0.2 kcal/mol [18]. When available, we further confirmed that the Δ E vac values of phenol dimers correlate well with the MP2 and coupled-cluster level data from previous studies [55,56].
To reliably predict Δ G aq , the solvation free energy, Δ G sol , must be accurately calculated using the thermodynamic cycle depicted in Figure 1. In this regard, we compared the estimated and observed values of Δ G sol for the benzene and phenol monomers. Several implicit solvation methods implemented in the Jaguar suite [54], including the Poisson–Boltzmann finite element method (PBF) [57]; polarizable continuum model (PCM) [58]; and Minnesota solvation models (SM6 and SM8) [59,60] were also evaluated. The results are summarized in Table 2, and we found that the DFT-CES values agree with the experimental values [61], in stark contrast with the implicit solvation methods, which predict Δ G sol with gross errors. It is clear that the existence of the OH group hugely determines the solvation behavior in water, accounting for significant differences of Δ G sol between benzene and phenol, and this behavior can be extrapolated to polyphenols such as quercetin [62].

3.2. Free Energy of Association in Aqueous Environment

Based on the successful descriptions of Δ E vac and Δ G sol enabled by PBE-D3 and DFT-CES, respectively, we investigated the Δ G aq of a benzene dimer for the ST, PD, T1, and T2 conformations (Figure 3a). We observed that the solvation effect slightly weakens the benzene dimer interaction for ST and PD conformations by 0.5–1 kcal/mol, while it remains almost unchanged for the T1 and T2 conformations. This alters the most favorable conformation of the benzene dimer; the PD and T conformations are almost equally stable in vacuum, but the T conformations preferentially attain stability in water.
Intriguingly, the solvent-accessible surface areas (SASAs) of the two T conformations were greater than those of the other two conformations, resulting in more water molecules surrounding the solute molecules (Figure 3b). In contrast, hydrophobic interactions occur between nonpolar and hydrophobic molecules, such as benzene, which tend to aggregate in water while minimizing their solvent-exposed surface. This is because the aromatic ring of benzene is not simply nonpolar and can develop two benzene–water hydrogen bonds on both sides of the phenyl ring (Figure 3c) [63]. As the association into T conformations only breaks one benzene–water hydrogen bond, while that into ST or PD conformations breaks two hydrogen bonds, the benzene dimer association in water favors the T conformations. Indeed, our earlier research indicated that breaking the benzene-water hydrogen bond requires approximately 1 kcal/mol [64], which nearly corresponds to the current difference in Δ G sol between the PD and T conformations.
Figure 4a shows the Δ G aq values of phenol dimers in various conformations. Note that the T2 (2,1′) and T2 (4,1′) conformations were excluded from further discussion because they were discovered to be extremely unstable in water, resulting in either a structure largely deviated from its original conformation or a positive value of Δ G aq (Figure S2). Except for the T2 (3,1′) and HBP conformations, we found smaller differences among different conformations and hydroxyl group positions. Moreover, Δ G aq values were comparable to Δ G vac values within a 1 kcal/mol margin. However, for the T2 (3,1′) and HBP conformations, the solvation effect significantly weakened the dimer association.
The major reason for the strongest binding in vacuum of the T2 (3,1′) and HBP conformations can be found in the direct interaction between the two OH groups. The HBP conformation creates an intermolecular HB between the two phenols. The T2 (3,1′) conformation also causes the OH group at the vertically upright phenyl ring to point toward the O atom of the OH group at the horizontally placed phenyl ring (see Figure 2), forming a weak HB. Because OH groups are solvated by water molecules via strong hydrogen bonds in the aqueous phase, when the association process develops a direct interaction between two OH groups (e.g., for the case of T2 (3,1′) and HBP), the loss of water–hydroxyl hydrogen bonds is inevitable. Figure 4b depicts the variation in the number of water–hydroxyl HB during the phenol dimer association. HB connectivity is shown in Figure S3. As the OH groups in the ST and PD conformations protrude toward the water phase, they can form ~2 HBs with water molecules for each OH group, as depicted in Figure 4c for the PD (4,1′) case. On the other hand, the OH groups in the T2 (3,1′) and HBP conformations mutually interact in the buried space, and only one solute–solvent hydrogen bond per OH group is permitted, as illustrated in Figure 4d. This is ascribed to the weakened interaction between T2 (3,1′) and HBP phenol dimers in water, whereas Δ G aq and Δ G vac are comparable for the other conformations. Nevertheless, despite the loss of water–hydroxyl hydrogen bonds, the HBP phenol dimer remains the most stable phenol dimer conformation in water.
We compared the values of Δ G vac and Δ G aq with the ring-to-ring distances between the dimers of benzene and phenol (Figure S5). Both Δ G vac and Δ G aq have negative relationships with the distances, indicating that the most stable configurations (i.e., T-shape for benzene dimer and HBP for phenol dimer) have the longest R-R distances. We also investigated the electrostatic potential surfaces of their configurations (Figure S6) and conclude that the ring components of benzene or phenol tend to become apart from each other to avoid the overlapping of charge densities, hence, developing strong π–π interactions.

3.3. Enthalpy–Entropy Compensation in Solvation Effect

Despite being non-negligible, the association free energy difference in vacuo and in aqua is mostly limited within 1.4 kcal/mol for all benzene and phenol dimer cases, except for the T2 (3,1′) and HBP phenol dimers, which have solute–solute hydrogen bonds (as also supported by the natural bond orbital analyses shown in Figure S7). As the association free energy difference is determined by the Δ Δ G sol , as illustrated in the thermodynamic cycle (Figure 1), we decomposed Δ Δ G sol into the enthalpic contribution, Δ Δ H sol , and the entropic contribution, T Δ Δ S sol , which are plotted in Figure 5.
Except for the two outliers, T2 (3,1′) and HBP phenol dimers, the enthalpy–entropy compensation behavior in Δ Δ G sol was well-defined. Moreover, as the corresponding Δ Δ H sol and T Δ Δ S sol values are distributed around the y = −x line, we discovered that the solvation effects of enthalpy and entropy cancel each other out at 300 K, yielding small differences between Δ G aq and Δ G vac within 1.4 kcal/mol. Interestingly, both Δ Δ H sol and Δ Δ S sol were negative for phenol dimers (except for T2 (3,1′) and HBP conformations), while Δ Δ H sol and Δ Δ S sol were positive for benzene dimers. Thus, it is inferred that enthalpy drives the coalescence of two hydration shells of phenol monomers into one, whereas entropy drives the fusion of benzene hydration shells. Consequently, the former can be classified as a nonclassical hydrophobic association, whereas the latter can be considered as an example of classical association.
During the dimerization of phenol into the T2 (3,1′) or HBP conformation, the enthalpy disfavors the coalescence of two hydration shells ( Δ Δ H sol > 0 ) owing to the loss of the hydroxyl–water hydrogen bond, but the entropy is insufficient to compensate ( Δ Δ S sol ~0). This causes the solvation thermodynamics of the T2 (3,1′) and HBP conformations to deviate from the trend established by the enthalpy–entropy compensation, resulting in a significantly weakened association free energy in water. However, despite this energy cost in water, the HBP conformation of phenol dimers has the strongest Δ G aq among all considered conformations; it is also stronger than the Δ G aq of the benzene dimer.

3.4. Relationships with Phase Behaviors of IDPs

Benzene and phenol are analogs of phenylalanine (Phe) and tyrosine (Tyr) side chains, respectively. Even with a small number of charged residues, a subset of IDPs containing a significant number of Phe/Tyr residues can drive liquid–liquid phase separation. One well-known example is the N-terminal low-complexity domain of FUS (fused in sarcoma), which contains 24 Tyr residues and few charged residues (2 Asp) [65].
The strongest benzene dimer interaction is −3.44 kcal/mol for the T1 conformation, and the strongest phenol dimer interaction is −4.93 kcal/mol for the HBP conformation; both are negative and sufficiently large, so they can promote the association, but the phenol dimer has a stronger binding. Other theoretical investigations [66,67] have demonstrated that the interaction strength of the associating units is positively correlated with the LLPS propensity. Indeed, previous experiments [23] have demonstrated that Phe and Tyr can replace each other without significantly altering the phase behaviors, but Tyr induces a stronger tendency toward LLPS than Phe [68,69].

4. Conclusions

In this study, the thermodynamics of benzene and phenol dimer association were investigated by determining the free energy of π–π interactions in water. Using vdW-corrected DFT and DFT-CES methods, we quantitatively estimated the free energy of dimer association in water, Δ G aq . Each component term (free energy of dimer association in vacuum and solvation free energy) was sufficiently accurate to be compared with available high-level coupled cluster-level calculations and experimental values. Interestingly, we discovered that the solvation effect only marginally modifies the π–π interaction free energies, unless a direct hydrogen-bond-like interaction develops between the two monomers. This trend can be explained in terms of enthalpy–entropy compensation, which clarifies that entropy and enthalpy drive the association of benzene and phenol dimers, respectively. We expect that this study will yield quantitatively reliable free-energetic values for π–π interactions in water and provide a fundamental understanding of the molecular driving force for the liquid–liquid phase separation of IDPs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23179811/s1.

Author Contributions

Conceptualization, H.K. and J.-M.C.; investigation, D.P. and H.L.; writing—original draft preparation, D.P.; writing—review and editing, H.L., H.K. and J.-M.C.; supervision, H.K. and J.-M.C.; funding acquisition, H.K. and J.-M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Research Foundation of Korea (NRF) grants funded by the Korean government (Grant No. 2017R1A5A1015365, 2021R1A2C2009643 and 2021R1C1C1010943).

Data Availability Statement

Not applicable.

Acknowledgments

We appreciate constructive feedback from two anonymous reviewers and valuable discussion with the members of the Kim group.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hunter, C.A.; Sanders, J.K. The nature of π-π interactions. J. Am. Chem. Soc. 1990, 112, 5525–5534. [Google Scholar] [CrossRef]
  2. Burley, S.; Petsko, G.A. Aromatic-aromatic interaction: A mechanism of protein structure stabilization. Science 1985, 229, 23–28. [Google Scholar] [CrossRef] [PubMed]
  3. McGaughey, G.B.; Gagné, M.; Rappé, A.K. π-Stacking Interactions: Alive and well in proteins. J. Biol. Chem. 1998, 273, 15458–15463. [Google Scholar] [CrossRef] [PubMed]
  4. Wilson, K.A.; Kung, R.W.; D’souza, S.; Wetmore, S.D. Anatomy of noncovalent interactions between the nucleobases or ribose and π-containing amino acids in RNA–protein complexes. Nucleic Acids Res. 2021, 49, 2213–2225. [Google Scholar] [CrossRef]
  5. Holmlin, R.E.; Dandliker, P.J.; Barton, J.K. Charge Transfer through the DNA Base Stack. Angew. Chem. Int. Ed. 1997, 36, 2714–2730. [Google Scholar] [CrossRef]
  6. Lewis, F.D.; Letsinger, R.L.; Wasielewski, M.R. Dynamics of Photoinduced Charge Transfer and Hole Transport in Synthetic DNA Hairpins. Acc. Chem. Res. 2001, 34, 159–170. [Google Scholar] [CrossRef]
  7. O’Neil, M.A.; Barton, J.K. DNA Charge Transport:  Conformationally Gated Hopping through Stacked Domains. J. Am. Chem. Soc. 2004, 126, 11471–11483. [Google Scholar] [CrossRef]
  8. Rutledge, L.R.; Campbell-Verduyn, L.S.; Hunter, K.C.; Wetmore, S.D. Characterization of Nucleobase−Amino Acid Stacking Interactions Utilized by a DNA Repair Enzyme. J. Phys. Chem. B 2006, 110, 19652–19663. [Google Scholar] [CrossRef]
  9. Rutledge, L.R.; Campbell-Verduyn, L.S.; Wetmore, S.D. Characterization of the stacking interactions between DNA or RNA nucleobases and the aromatic amino acids. Chem. Phys. Lett. 2007, 444, 167–175. [Google Scholar] [CrossRef]
  10. Zhuang, W.R.; Wang, Y.; Cui, P.F.; Xing, L.; Lee, J.; Kim, D.; Jiang, H.L.; Oh, Y.K. Applications of π-π stacking interactions in the design of drug-delivery systems. J. Control Release 2019, 294, 311–326. [Google Scholar] [CrossRef]
  11. Hobza, P.; Selzle, H.; Schlag, E. Floppy structure of the benzene dimer: Ab initio calculation on the structure and dipole moment. J. Chem. Phys. 1990, 93, 5893–5897. [Google Scholar] [CrossRef]
  12. Hobza, P.; Selzle, H.L.; Schlag, E.W. New structure for the most stable isomer of the benzene dimer: A quantum chemical study. J. Phys. Chem. 1993, 97, 3937–3938. [Google Scholar] [CrossRef]
  13. Hobza, P.; Selzle, H.L.; Schlag, E.W. Potential energy surface of the benzene dimer: Ab initio theoretical study. J. Am. Chem. Soc. 1994, 116, 3500–3506. [Google Scholar] [CrossRef]
  14. Hobza, P.; Selzle, H.L.; Schlag, E.W. Potential energy surface for the benzene dimer. Results of ab initio CCSD (T) calculations show two nearly isoenergetic structures: T-shaped and parallel-displaced. J. Phys. Chem. 1996, 100, 18790–18794. [Google Scholar] [CrossRef]
  15. Sinnokrot, M.O.; Valeev, E.F.; Sherrill, C.D. Estimates of the ab initio limit for π–π interactions: The benzene dimer. J. Am. Chem. Soc. 2002, 124, 10887–10893. [Google Scholar] [CrossRef] [PubMed]
  16. Sinnokrot, M.O.; Sherrill, C.D. Unexpected substituent effects in face-to-face π-stacking interactions. J. Phys. Chem. A 2003, 107, 8377–8379. [Google Scholar] [CrossRef]
  17. Sinnokrot, M.O.; Sherrill, C.D. Highly accurate coupled cluster potential energy curves for the benzene dimer: Sandwich, T-shaped, and parallel-displaced configurations. J. Phys. Chem. A 2004, 108, 10200–10207. [Google Scholar] [CrossRef]
  18. Sinnokrot, M.O.; Sherrill, C.D. High-accuracy quantum mechanical studies of π–π interactions in benzene dimers. J. Phys. Chem. A 2006, 110, 10656–10668. [Google Scholar] [CrossRef]
  19. Burger, V.M.; Gurry, T.; Stultz, C.M. Intrinsically Disordered Proteins: Where Computation Meets Experiment. Polymers 2014, 6, 2684–2719. [Google Scholar] [CrossRef]
  20. Kato, M.; Han, T.W.; Xie, S.; Shi, K.; Du, X.; Wu, L.C.; McKnight, S.L. Cell-free formation of RNA granules: Low complexity sequence domains form dynamic fibers within hydrogels. Cell 2012, 149, 753–767. [Google Scholar] [CrossRef] [Green Version]
  21. Jiang, H.; Wang, S.; Huang, Y.; He, X.; Cui, H.; Zhu, X.; Zheng, Y. Phase transition of spindle-associated protein regulate spindle apparatus assembly. Cell 2015, 163, 108–122. [Google Scholar] [CrossRef] [PubMed]
  22. Gabryelczyk, B.; Cai, H.; Shi, X.; Sun, Y.; Swinkels, P.J.M.; Salentinig, S.; Pervushin, K.; Miserez, A. Hydrogen bond guidance and aromatic stacking drive liquid-liquid phase separation of intrinsically disordered histidine-rich peptides. Nat. Commun. 2019, 10, 5465. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, J.; Choi, J.M.; Holehouse, A.S.; Lee, H.O.; Zhang, X.; Jahnel, M.; Maharana, S.; Lemaitre, R.; Pozniakovsky, A.; Drechsel, D.; et al. A Molecular Grammar Governing the Driving Forces for Phase Separation of Prion-like RNA Binding Proteins. Cell 2018, 174, 688–699. [Google Scholar] [CrossRef]
  24. Martin, E.W.; Holehouse, A.S.; Peran, I.; Farag, M.; Incicco, J.J.; Bremer, A.; Grace, C.R.; Soranno, A.; Pappu, R.V.; Mittag, T. Valence and patterning of aromatic residues determine the phase behavior of prion-like domains. Science 2020, 367, 694–699. [Google Scholar] [CrossRef] [PubMed]
  25. Bremer, A.; Farag, M.; Borcherds, W.M.; Peran, I.; Martin, E.W.; Pappu, R.V.; Mittag, T. Deciphering how naturally occurring sequence features impact the phase behaviours of disordered prion-like domains. Nat. Chem. 2022, 14, 196–207. [Google Scholar] [CrossRef]
  26. Lim, H.K.; Lee, H.; Kim, H. A Seamless Grid-Based Interface for Mean-Field QM/MM Coupled with Efficient Solvation Free Energy Calculations. J. Chem. Theory Comput. 2016, 12, 5088–5099. [Google Scholar] [CrossRef]
  27. Lee, H.; Lim, H.K.; Kim, H. Hydration Thermodynamics of Non-Polar Aromatic Hydrocarbons: Comparison of Implicit and Explicit Solvation Models. Molecules 2018, 23, 2927. [Google Scholar] [CrossRef]
  28. Jang, T.; Paik, D.; Shin, S.J.; Kim, H. Density functional theory in classical explicit solvents: Mean-field QM/MM method for simulating solid–liquid interfaces. Bull. Korean Chem. Soc. 2022, 43, 476–483. [Google Scholar] [CrossRef]
  29. Nakano, H.; Yamamoto, T. Accurate and efficient treatment of continuous solute charge density in the mean-field QM/MM free energy calculation. J. Chem. Theory Comput. 2013, 9, 188–203. [Google Scholar] [CrossRef] [PubMed]
  30. Galván, I.F.; Sánchez, M.; Martın, M.; del Valle, F.O.; Aguilar, M. ASEP/MD: A program for the calculation of solvent effects combining QM/MM methods and the mean field approximation. Comput. Phys. Commun. 2003, 155, 244–259. [Google Scholar] [CrossRef]
  31. Sanchez, M.L.; Aguilar, M.; Olivares del Valle, F. Study of solvent effects by means of averaged solvent electrostatic potentials obtained from molecular dynamics data. J. Comput. Chem. 1997, 18, 313–322. [Google Scholar] [CrossRef]
  32. Yamamoto, T. Variational and perturbative formulations of quantum mechanical/molecular mechanical free energy with mean-field embedding and its analytical gradients. J. Chem. Phys. 2008, 129, 244104. [Google Scholar] [CrossRef] [PubMed]
  33. Lin, S.-T.; Blanco, M.; Goddard, W.A., III. The two-phase model for calculating thermodynamic properties of liquids from molecular dynamics: Validation for the phase diagram of Lennard-Jones fluids. J. Chem. Phys. 2003, 119, 11792–11805. [Google Scholar] [CrossRef]
  34. Lin, S.-T.; Maiti, P.K.; Goddard, W.A., III. Two-phase thermodynamic model for efficient and accurate absolute entropy of water from molecular dynamics simulations. J. Phys. Chem. B 2010, 114, 8191–8198. [Google Scholar] [CrossRef]
  35. Pascal, T.A.; Lin, S.-T.; Goddard, W.A., III. Thermodynamics of liquids: Standard molar entropies and heat capacities of common solvents from 2PT molecular dynamics. Phys. Chem. Chem. Phys. 2011, 13, 169–181. [Google Scholar] [CrossRef] [PubMed]
  36. Pascal, T.A.; Goddard, W.A., III. Interfacial thermodynamics of water and six other liquid solvents. J. Phys. Chem. B 2014, 118, 5943–5956. [Google Scholar] [CrossRef]
  37. Lim, H.K.; Kwon, Y.; Kim, H.S.; Jeon, J.; Kim, Y.H.; Lim, J.A.; Kim, B.S.; Choi, J.; Kim, H. Insight into the Microenvironments of the Metal-Ionic Liquid Interface during Electrochemical CO2 Reduction. ACS Catal. 2018, 8, 2420–2427. [Google Scholar] [CrossRef]
  38. Gim, S.; Lim, H.K.; Kim, H. Multiscale Simulation Method for Quantitative Prediction of Surface Wettability at the Atomistic Level. J. Phys. Chem. Lett. 2018, 9, 1750–1758. [Google Scholar] [CrossRef]
  39. Gim, S.; Cho, K.J.; Lim, H.K.; Kim, H. Structure, Dynamics, and Wettability of Water at Metal Interfaces. Sci. Rep. 2019, 9, 14805. [Google Scholar] [CrossRef]
  40. Shin, S.J.; Kim, D.H.; Bae, G.; Ringe, S.; Choi, H.; Lim, H.K.; Choi, C.H.; Kim, H. On the importance of the electric double layer structure in aqueous electrocatalysis. Nat. Commun. 2022, 13, 174. [Google Scholar] [CrossRef]
  41. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.L.; Cococcioni, M. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 2009, 21, 39. [Google Scholar] [CrossRef] [PubMed]
  42. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef]
  43. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  44. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [PubMed]
  45. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef] [PubMed]
  46. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B Condens. Matter Mater. Phys. 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  47. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B Condens. Matter Mater. Phys. 1994, 50, 17953–17979. [Google Scholar] [CrossRef]
  48. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef]
  49. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A At. Mol. Opt. Phys. 1985, 31, 1695–1697. [Google Scholar] [CrossRef]
  50. Price, D.J.; Brooks, C.L., III. A modified TIP3P water potential for simulation with Ewald summation. J. Chem. Phys. 2004, 121, 10096–10103. [Google Scholar] [CrossRef]
  51. Andersen, H.C. Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations. J. Comput. Phys. 1983, 52, 24–34. [Google Scholar] [CrossRef]
  52. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N· log (N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef]
  53. Kaminski, G.A.; Friesner, R.A.; Tirado-Rives, J.; Jorgensen, W.L. Evaluation and reparametrization of the OPLS-AA force field for proteins via comparison with accurate quantum chemical calculations on peptides. J. Phys. Chem. B 2001, 105, 6474–6487. [Google Scholar] [CrossRef]
  54. Bochevarov, A.D.; Harder, E.; Hughes, T.F.; Greenwood, J.R.; Braden, D.A.; Philipp, D.M.; Rinaldo, D.; Halls, M.D.; Zhang, J.; Friesner, R.A. Jaguar: A High-Performance Quantum Chemistry Software Program with Strengths in Life and Materials Sciences. Int. J. Quantum Chem. 2013, 113, 2110–2142. [Google Scholar] [CrossRef]
  55. Ali, I.; Sharma, S.; Bezbaruah, B. Quantum Mechanical Study on the π-π Stacking Interaction and Change in Conformation of Phenolic Systems with Different Intermolecular Rotations. Comput. Chem. 2018, 6, 71–86. [Google Scholar] [CrossRef]
  56. Kołaski, M.; Kumar, A.; Singh, N.J.; Kim, K.S. Differences in structure, energy, and spectrum between neutral, protonated, and deprotonated phenol dimers: Comparison of various density functionals with ab initio theory. Phys. Chem. Chem. Phys. 2011, 13, 991–1001. [Google Scholar] [CrossRef]
  57. Chen, L.; Holst, M.J.; Xu, J. The Finite Element Approximation of the Nonlinear Poisson–Boltzmann Equation. SIAM J. Numer. Anal. 2007, 45, 2298–2320. [Google Scholar] [CrossRef]
  58. Mennucci, B. Polarizable continuum model. WIREs Comput. Mol. Sci. 2012, 2, 386–404. [Google Scholar] [CrossRef]
  59. Kelly, C.P.; Cramer, C.J.; Truhlar, D.G. SM6: A Density Functional Theory Continuum Solvation Model for Calculating Aqueous Solvation Free Energies of Neutrals, Ions, and Solute-Water Clusters. J. Chem. Theory Comput. 2005, 1, 1133–1152. [Google Scholar] [CrossRef]
  60. Marenich, A.V.; Olson, R.M.; Kelly, C.P.; Cramer, C.J.; Truhlar, D.G. Self-Consistent Reaction Field Model for Aqueous and Nonaqueous Solutions Based on Accurate Polarized Partial Charges. J. Chem. Theory Comput. 2007, 3, 2011–2033. [Google Scholar] [CrossRef]
  61. Rankin, K.N.; Sulea, T.; Purisima, E.O. On the transferability of hydration-parametrized continuum electrostatics models to solvated binding calculations. J. Comput. Chem. 2003, 24, 954–962. [Google Scholar] [CrossRef] [PubMed]
  62. Yu, B.; Ren, H.; Piao, X. Towards Adsorptive Enrichment of Flavonoids from Honey Using h-BN Monolayer. ChemPhysChem 2022, 23, e202100828. [Google Scholar] [CrossRef] [PubMed]
  63. Suzuki, S.; Green, P.G.; Bumgarner, R.E.; Blake, G.A.; Dasgupta, S.; Goddard, W.A., III. Benzene Forms Hydrogen Bonds with Water. Science 1992, 257, 942–945. [Google Scholar] [CrossRef] [PubMed]
  64. Lee, H.; Dehez, F.; Chipot, C.; Lim, H.K.; Kim, H. Enthalpy–Entropy Interplay in π-Stacking Interaction of Benzene Dimer in Water. J. Chem. Theory Comput. 2019, 15, 1538–1545. [Google Scholar] [CrossRef]
  65. Burke, K.A.; Janke, A.M.; Rhine, C.L.; Fawzi, N.L. Residue-by-Residue View of In Vitro FUS Granules that Bind the C-Terminal Domain of RNA Polymerase II. Molecular Cell 2015, 60, 231–241. [Google Scholar] [CrossRef] [PubMed]
  66. Choi, J.M.; Holehouse, A.S.; Pappu, R.V. Physical Principles Underlying the Complex Biology of Intracellular Phase Transitions. Annu. Rev. Biophys. 2020, 49, 107–133. [Google Scholar] [CrossRef]
  67. Choi, J.M.; Hyman, A.A.; Pappu, R.V. Generalized models for bond percolation transitions of associative polymers. Phys. Rev. E 2020, 102, 042403. [Google Scholar] [CrossRef] [PubMed]
  68. Frey, S.; Richter, R.P.; Görlich, D. FG-rich repeats of nuclear pore proteins form a three-dimensional meshwork with hydrogel-like properties. Science 2006, 314, 815–817. [Google Scholar] [CrossRef]
  69. Lin, Y.; Currie, S.L.; Rosen, M.K. Intrinsically disordered sequences enable modulation of protein phase separation through distributed tyrosine motifs. J. Biol. Chem. 2017, 292, 19110–19120. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Thermodynamic cycle describing aromatic π–π interaction in aqueous solution. The free energy of association in water ( Δ G aq ) can be calculated using the free energy of association in vacuum ( Δ G vac ) and the solvation free energy difference ( Δ Δ G sol ) between a dimer ( Δ G dim sol ) and two monomers ( Δ G mon sol ). Different colors indicate different atom types: white—hydrogen; brown—carbon; and red—oxygen (same for the following figures).
Figure 1. Thermodynamic cycle describing aromatic π–π interaction in aqueous solution. The free energy of association in water ( Δ G aq ) can be calculated using the free energy of association in vacuum ( Δ G vac ) and the solvation free energy difference ( Δ Δ G sol ) between a dimer ( Δ G dim sol ) and two monomers ( Δ G mon sol ). Different colors indicate different atom types: white—hydrogen; brown—carbon; and red—oxygen (same for the following figures).
Ijms 23 09811 g001
Figure 2. (a) Geometric conformations of benzene dimer. Each carbon from the upper (1~6) and lower (1′~6′) benzene is assigned a number. (b) Geometric conformations of phenol dimer. Each structure is designated by the position of the ring component (ST, PD, T1, T2, or HBP) and the number of substituted carbons with a hydroxyl group from the upper (1~6) and lower (1′~6′) phenol.
Figure 2. (a) Geometric conformations of benzene dimer. Each carbon from the upper (1~6) and lower (1′~6′) benzene is assigned a number. (b) Geometric conformations of phenol dimer. Each structure is designated by the position of the ring component (ST, PD, T1, T2, or HBP) and the number of substituted carbons with a hydroxyl group from the upper (1~6) and lower (1′~6′) phenol.
Ijms 23 09811 g002
Figure 3. (a) Association free energy of benzene dimers in vacuum ( Δ G vac ) and in solution ( Δ G aq ). (b) Average number of water molecules in the first solvent shell ( N wat ) and solvent accessible surface area (SASA) corresponding to the benzene dimer association in water. (c) Schematic illustrating the benzene–water hydrogen bonds.
Figure 3. (a) Association free energy of benzene dimers in vacuum ( Δ G vac ) and in solution ( Δ G aq ). (b) Average number of water molecules in the first solvent shell ( N wat ) and solvent accessible surface area (SASA) corresponding to the benzene dimer association in water. (c) Schematic illustrating the benzene–water hydrogen bonds.
Ijms 23 09811 g003
Figure 4. (a) Association free energy of phenol dimers in vacuum ( Δ G vac ) and in water ( Δ G aq ). (b) Average number of water–hydroxyl hydrogen bonds present in phenol dimers. The dashed line represents twice the average number of hydrogen bonds of phenol monomer. (c,d) Schematic illustrating the water–hydroxyl hydrogen bonds of (c) PD (4,1′) and (d) HBP conformation. The numbers are H-bond distances in angstroms, which are taken from the first peak locations of radial distribution functions shown in Figure S4.
Figure 4. (a) Association free energy of phenol dimers in vacuum ( Δ G vac ) and in water ( Δ G aq ). (b) Average number of water–hydroxyl hydrogen bonds present in phenol dimers. The dashed line represents twice the average number of hydrogen bonds of phenol monomer. (c,d) Schematic illustrating the water–hydroxyl hydrogen bonds of (c) PD (4,1′) and (d) HBP conformation. The numbers are H-bond distances in angstroms, which are taken from the first peak locations of radial distribution functions shown in Figure S4.
Ijms 23 09811 g004
Figure 5. Enthalpic ( Δ Δ H sol ) and entropic contributions ( T Δ Δ S sol ) of the solvation free energy difference ( Δ Δ G sol ) between a dimer and two monomers. The dashed line represents y = x .
Figure 5. Enthalpic ( Δ Δ H sol ) and entropic contributions ( T Δ Δ S sol ) of the solvation free energy difference ( Δ Δ G sol ) between a dimer and two monomers. The dashed line represents y = x .
Ijms 23 09811 g005
Table 1. Estimated binding energies ( Δ E vac ) and binding free energies ( Δ G vac ) of benzene and phenol dimer in vacuum. The listed values are in kcal/mol. Δ E vac from other studies were calculated using (a) estimated CCSD(T)/aug-cc-pVQZ*, which is a modified aug-cc-pVQZ basis set lacking g functions on carbon and f functions on hydrogen [18]. (b) MP2/M062X/6-311++G(d,p) [55]. (c) CCSD(T)/CBS [56].
Table 1. Estimated binding energies ( Δ E vac ) and binding free energies ( Δ G vac ) of benzene and phenol dimer in vacuum. The listed values are in kcal/mol. Δ E vac from other studies were calculated using (a) estimated CCSD(T)/aug-cc-pVQZ*, which is a modified aug-cc-pVQZ basis set lacking g functions on carbon and f functions on hydrogen [18]. (b) MP2/M062X/6-311++G(d,p) [55]. (c) CCSD(T)/CBS [56].
Present StudyOthers
Δ G vac Δ E vac Δ E vac
Benzene
Dimer
ST−2.34−1.70−1.70 (a)
PD−3.11−2.43−2.63 (a)
T1−3.22−2.66−2.61 (a)
T2−3.25−2.66
Phenol
Dimer
ST (1,1′)−2.46−1.69−1.61 (b)
ST (1,4′)−3.11−2.42−2.76 (b)
PD (1,1′)−3.19−2.45−3.62 (b)
PD (1,4′)−3.71−3.00
PD (4,1′)−3.14−2.47−4.26 (b)
T1 (1,1′)−3.12−2.57
T1 (1,2′)−3.18−2.64
T1 (1,3′)−3.34−2.83
T2 (1,1′)−3.31−2.79
T2 (1,2′)−3.10−2.56
T2 (1,3′)−3.26−2.75
T2 (1,4′)−3.28−2.77
T2 (2,1′)−3.50−2.99
T2 (3,1′)−5.90−5.66
T2 (4,1′)−3.65−3.12
HBP−7.29−6.87−6.81 (c)
Table 2. Calculated solvation free energies ( Δ G sol ) of benzene and phenol monomer using PBF, PCM, SM6, SM8, and DFT-CES method. The listed values are in kcal/mol. The experimental values of Δ G sol are from the hydration dataset of a previous study about solvated binding free energies [61].
Table 2. Calculated solvation free energies ( Δ G sol ) of benzene and phenol monomer using PBF, PCM, SM6, SM8, and DFT-CES method. The listed values are in kcal/mol. The experimental values of Δ G sol are from the hydration dataset of a previous study about solvated binding free energies [61].
BenzenePhenol
PBF−0.09−6.72
PCM−2.52−5.45
SM6−3.27−5.60
SM8−3.61−5.63
DFT-CES−0.85−6.72
Experiment−0.87−6.62
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Paik, D.; Lee, H.; Kim, H.; Choi, J.-M. Thermodynamics of π–π Interactions of Benzene and Phenol in Water. Int. J. Mol. Sci. 2022, 23, 9811. https://doi.org/10.3390/ijms23179811

AMA Style

Paik D, Lee H, Kim H, Choi J-M. Thermodynamics of π–π Interactions of Benzene and Phenol in Water. International Journal of Molecular Sciences. 2022; 23(17):9811. https://doi.org/10.3390/ijms23179811

Chicago/Turabian Style

Paik, Dooam, Hankyul Lee, Hyungjun Kim, and Jeong-Mo Choi. 2022. "Thermodynamics of π–π Interactions of Benzene and Phenol in Water" International Journal of Molecular Sciences 23, no. 17: 9811. https://doi.org/10.3390/ijms23179811

APA Style

Paik, D., Lee, H., Kim, H., & Choi, J. -M. (2022). Thermodynamics of π–π Interactions of Benzene and Phenol in Water. International Journal of Molecular Sciences, 23(17), 9811. https://doi.org/10.3390/ijms23179811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop