Next Article in Journal / Special Issue
Cell-Adhesive Bioinspired and Catechol-Based Multilayer Freestanding Membranes for Bone Tissue Engineering
Previous Article in Journal / Special Issue
Size Control and Fluorescence Labeling of Polydopamine Melanin-Mimetic Nanoparticles for Intracellular Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Noncovalent Interactions in the Catechol Dimer

1
Scuola Normale Superiore di Pisa, Piazza dei Cavalieri, I-56126 Pisa, Italy
2
Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
3
Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy
*
Author to whom correspondence should be addressed.
Biomimetics 2017, 2(3), 18; https://doi.org/10.3390/biomimetics2030018
Submission received: 21 June 2017 / Revised: 4 September 2017 / Accepted: 5 September 2017 / Published: 13 September 2017
(This article belongs to the Special Issue Bioinspired Catechol-Based Systems: Chemistry and Applications)

Abstract

:
Noncovalent interactions play a significant role in a wide variety of biological processes and bio-inspired species. It is, therefore, important to have at hand suitable computational methods for their investigation. In this paper, we report on the contribution of dispersion and hydrogen bonds in both stacked and T-shaped catechol dimers, with the aim of delineating the respective role of these classes of interactions in determining the most stable structure. By using second-order Møller–Plesset (MP2) calculations with a small basis set, specifically optimized for these species, we have explored a number of significant sections of the interaction potential energy surface and found the most stable structures for the dimer, in good agreement with the highly accurate, but computationally more expensive coupled cluster single and double excitation and the perturbative triples (CCSD(T))/CBS) method.

Graphical Abstract

1. Introduction

Nowadays, there is a general consensus about the primary role played by noncovalent interactions, in particular those involving aromatic rings, in molecular, life, and materials sciences. In addition to being responsible for key biological processes that range from base stacking in deoxyribonucleic acid (DNA) [1], to the color of red wine [2] and, more generally, food quality [3], it is of the foremost importance to understand, rationalize and, hence, exploit their features in cutting-edge applications as advanced catalysis [4,5], biomedical materials [6,7] and novel drugs design [8], advanced organic photovoltaics [9,10,11,12,13], complex self-assembled structures [14], or bio-nano-materials [15,16]. Such ubiquity of the aromatic interactions has often inspired multidisciplinary research [17], aimed to exploit their peculiar features in the design and construction of biomimetic materials. From a physical point of view, noncovalent interactions among molecules bearing aromatic moieties originate from a variety of different forces, including π-stacking, XH–π or charge-transfer (CT), besides the ubiquitous dispersion. Furthermore, the presence of additional functional groups can introduce other kinds of interactions (like e.g., hydrogen (HB) or halogen bonds), leading to nontrivial interference effects, which tune both the structure and the properties of the resulting material. In this framework, computational methods can play a crucial role for rational design and interpretation, provided that they are able to couple reliability, feasibility, and ability to unravel the different contributions [18,19]. It should be also mentioned that, although the embedding environment is often neglected, or only roughly approximated in most computational studies, its effect can be significant or even decisive in biomimetic processes. However, comprehensive studies of pairs of interacting species in the gas phase are a mandatory starting point for unraveling the weight of the different effects.
In the past few years, catechol has attracted increasing attention as a precursor of bio-inspired materials [20,21,22,23,24,25,26]. From a theoretical point of view, catechol is an ideal candidate to test the capability of new computational approaches to accurately represent the delicate balance among the different kinds of noncovalent interactions, occurring in the presence of catechol units. In fact, apart from the π-stacking and XH–π interactions due to the aromatic core, interactions between these species are also characterized by the insurgence of strong (OH–H) and weak (OH–π) HB patterns, which may play an important role in the supramolecular assembling. The main problem is that aromatic interactions are dominated by dispersion forces that standard electronic calculations have difficulty to reproduce. Indeed, in the past ten years, much effort has been devoted to the development of approaches that overcome the problem [27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43]. Within the framework of density functional theory (DFT), attempts have been made to set appropriate functionals which incorporate the effects of dispersion, such as that of Truhlar et al. [43] or to introduce semi-empirical atomistic corrections, as suggested by Grimme and coworkers [30,32,33]. Among wave function (WF)-based approaches, the most accurate but also computationally most expensive method is the coupled cluster approach including a full account of single and double excitations together with perturbative inclusion of connected triple excitations, and extrapolation to the complete basis set limit (CCSD(T)/CBS) [4,19,34,37,38,41,44,45,46,47,48,49,50]. Still within a WF framework, perturbative second-order Møller–Plesset (MP2) calculations could be carried out at a much lower computational cost, yet it is well known [40] that they tend to overestimate aromatic binding energies, especially when employed with large basis sets. These inaccuracies can be overcome by resorting to an idea proposed almost forty years ago by Kroon-Batenburg and Van Duijneveldt [51] and successively refined by Hobza and Zahradnik [52], based on the use MP2 calculations with the small 6-31G* basis set, modified by reducing to 0.25 the exponent of the d polarization function placed on each carbon atom of the benzene dimer. Such an approach, often referred to as MP2/6-31G*(0.25), was then fully validated with reference to interaction energies of benzene and a few other aromatic dimers computed at the CCSD(T)/CBS level [53,54,55,56,57,58,59,60,61]. More recently, the method has been generalized to different basis sets, and applied to several molecular prototypes, including liquid crystals [62,63], pyridine [64], quinhydrone [27], dihydroxyindole derivatives relevant in eumelanin formation [65], and, very recently, to small aromatic heterocycles [66], where the procedure to find the suitable modified basis sets, labeled MP2mod, has been automated and extended to the optimization of the orbital exponents of d functions on heteroatoms and p functions on hydrogen, within the 6-31G** basis set.
Here, the MP2mod method is applied to the catechol dimer in the gas phase. First, MP2mod accuracy is validated against high-quality CCSD(T)/CBS predictions, purposely carried out for a number of selected geometries of catechol dimers. Next, MP2mod is employed in the exploration of the catechol’s interaction potential energy surface (IPES), with the aim of finding the optimal structure of the dimer by a comparison of different possible arrangements. This allows us to investigate the different roles played by HB and π-stacking interactions in the dimer formation. Incidentally, it might also be of interest, following Wheeler group’s suggestions [4,44,45], to verify if noncovalent interactions in catechol can be correlated to the simple direct interaction between the (hydroxyl) substituents, or if, on the contrary, a rationalization of the resulting interaction patterns requires a more complex analysis, taking into account the specific role of each contribution.
The catechol dimer has also been studied at the DFT level by Estévez et al. [67], who considered structures determined either by X-ray measurements or by geometry optimizations at the MPW1B95/6-311++G(2d,2p) level. In the following these results will also be discussed in comparison with our findings.

2. Computational Details

The full geometry optimization of the catechol monomer has been performed by DFT, at the B3LYP/aug-cc-pvDZ level, by minimizing the energy with respect to all internal coordinates. Unless otherwise stated, the internal monomer’s geometry was kept unaltered in all subsequent calculations.
As far as the intermolecular energy is concerned, reference CCSD(T)/CBS calculations have been carried out on catechol dimers following the protocol adopted in previous works [27,36,66], which can be summarized as follows:
  • The difference ΔCC-MP2 between CCSD(T) and MP2 interaction energy is evaluated using for both calculations the Dunning’s correlated aug-cc-pvDz basis sets:
    Δ C C M P 2 = | Δ E C C S D ( T ) | a u g c c p v D z | Δ E M P 2 | a u g c c p v D z
  • The MP2 energy in the CBS limit, Δ E C B S M P 2 , is computed through the extrapolation scheme proposed by Halkier et al. [68], making use of the aug-cc-pvDz and aug-cc-pvTz basis sets. Despite the state-of-the-art extrapolation procedure [37,41,50] is often carried out with the larger aug-cc-pvTz and aug-cc-pvQz basis sets, it has been recently shown that, for similar aromatic dimers, the use of the smaller aug-cc-pvDz and aug-cc-pvTz affects the computed interaction energies by few hundredths of kcal/mol [66]. In consideration of the fairly large number of dimers investigated and the computational cost of a CCSD(T) calculation at the aug-cc-pvQz level, the smaller sets (Dz and Tz) were chosen as the best compromise between accuracy and feasibility.
  • Finally, the CCSD(T)/CBS interaction energy, Δ E C B S C C S D ( T ) , is recovered as:
    Δ E C B S C C S D ( T ) =   Δ E C B S M P 2 +   Δ C C M P 2
  • All energies were corrected for the basis set superposition error (BSSE) with the standard counterpoise (CP) correction [69].
The MP2mod exponent optimization was performed by means of the Exopt code [27,36,66], by minimizing the objective function I:
I ( P ¯ ) =   1 N g e o m k = 1 N g e o m [ Δ E C B S C C S D ( T ) Δ E M P 2 m o d ( P ¯ ) ] 2
where Ngeom is the number of considered dimer geometries and P ¯ the vector containing the basis sets exponents to be optimized. All the MP2mod calculations were carried out with the 6-31G** basis set, and the exponents of the d functions on heavy atoms and the p functions on H were optimized. Further details on the optimization protocol can be found in [66] and are also briefly commented in the next section. In all MP2mod calculations, the CP correction was applied to take care of the basis set superposition error.
Finally, to better compare with the results reported by Estévez et al. [67], the interaction energy of selected dimer arrangements was also computed at the DFT level, using the same procedure employed in [67]: the MPW1B95 functional was employed, together with the 6-311++G(2d,2p), while no correction was applied to take care of the BSSE.
All CCSD(T), MP2, MP2mod and DFT calculations were carried out with the Gaussian09 software package [70].

3. Results and Discussion

3.1. MP2mod Tuning and Validation

After geometry optimization, the catechol monomer is planar with the two hydroxyl hydrogens pointing in the same direction (see Figure 1a).
Based on the results recently achieved for several heteroaromatic dimers, where stacked and T-shaped (TS) conformers where found to be the most stable, four starting arrangements have been set up by placing the two monomers at different distances and relative orientation. Namely, the face-to-face (FF, Figure 1b), the antiparallel face-to-face (AFF, Figure 1c), and two TS conformations, one with both hydroxyls (TS1, Figure 1d) and one with only one hydroxyl (TS2, Figure 1e) pointing towards the other ring. Following the protocol recently developed in our group [66], the MP2mod best exponents were determined as follows: starting from each of the four selected conformations, a set of dimer arrangements was created by displacing one monomer along a selected coordinate R, defined as the line connecting the centers of the two rings, as shown in the insets of Figure 2. Next, an estimate (data not shown) of the interaction energy (ΔE) of the resulting dimer geometries was obtained at the MP2mod level, employing the basis set recently optimized by us for quinhydrone [27], thus obtaining preliminary interaction energy profiles. Three points (displayed as blue squares in Figure 2) were selected for each profile (namely one in the minimum, one in the short distance range and one in the attractive branch of the curve) and the corresponding CCSD(T)/CBS interaction energies were computed and used to build a reference database containing 12 elements. This database was then used for the optimization of the exponents of the polarization functions of the 6-31G** basis sets suitable for MP2mod calculations. The starting exponents of the standard 6-31G** basis set are 0.80 for d functions on carbon and oxygen and 1.1 for p functions on hydrogen. After optimization, the best exponents were found to be 0.27 and 0.34 for the d functions on carbon and oxygen, respectively, and 0.36 for p functions on hydrogen. The final standard deviation, I , see Equation (3), resulted to be less than 0.3 kcal/mol with respect to the CCSD(T)/CBS energies.
The resulting MP2mod curves are displayed in Figure 2, together with the reference values. The excellent agreement between the two methods, in line with the results previously obtained for similar molecules, allows us to apply rather confidently the MP2mod method to the study of the catechol dimer. According to both CCSD(T)/CBS and MP2mod results, the most stable structure is the TS2 one (around −5.0 kcal/mol), with the minimum at a slightly smaller value of R (5.4 Å), with respect to the similar TS1 conformer (5.6 Å), which is in turn almost as stable (≈−4.0 kcal/mol) as the antiparallel stacked conformer (AFF, −3.8 kcal/mol). Among the two stacked conformations, FF and AFF, the second one is more stable, in agreement with the repulsive interaction between the OH dipoles in the FF form.

3.2. Stacked Cathecol Dimers

Due to its importance, the stacked arrangement has been studied with some care as a function of the ring–ring distance R and of the angle β, which expresses, as shown in Figure 3b, the relative rotation of the two rings with respect to the line connecting their centers. The relevant results are reported in Figure 3. In the left panel, the interaction energy, reported vs. R for assigned rotation angles, shows minima at similar R values for all angles, and a marked dependence on β at low vales (from 0 to 60°), whereas for β > 90° the curves are close to each other: at the minimum the interaction energy changes by only ≈0.25 kcal/mol in the range 90–180°. Although this behavior seems roughly consistent with a dipole–dipole interaction, the resemblance of the 90, 120, 150, and 180° curves is an indication that higher multipoles, or, equivalently, local dipoles, should play a role in an electrostatic rationalization of the observed energy curves. This is in agreement with the idea of Wheeler and coworkers [45,46] that stacking interaction in substituted aromatic species is strongly influenced by the local interaction of the substituents, rather than to changes induced in the π electronic density upon substitution, as suggested by older models.
Figure 3b shows the energy variation as a function of β and connects the FF (β = 0°) to the AFF (β = 180°) arrangement at a fixed ring-ring distance (R = 3.5 Å). The curve shows a not monotonic behavior, probably due to the presence of two functional groups, with an absolute minimum near 110°, rather than at 180°, as could be expected for single substituted benzene rings. However, despite the perturbations triggered by the specific interaction among the two strong local dipoles of the monomers, the transition from FF to AFF arrangements along β is rather marked and clearly indicates a preference for antiparallel stacked arrangements, as already put in evidence in Figure 2.
In order to gain a deeper insight into the orientation dependence of the stacking forces in the catechol dimer, taking advantage from the low computational cost of the MP2mod method, we can explore different sections of the catechol IPES. For instance, in Figure 4 a two-dimensional contour plot of the interaction energy (ΔE) is reported as a function of the horizontal displacement of the two rings (R) and of the rotation angle (β) of one of the two rings around the perpendicular axis, at the inter-ring distance of 3.5 Å (i.e., the position of the minimum for the stacked energy curves reported in Figure 3).
Figure 4 clearly shows that the dimer is much more stable when displaced and rotated with respect to the FF arrangement, with a minimum at R ≈ 1.2 Å and β ≈ 130°. It is noteworthy that the effects of horizontal displacement (i.e., varying R) and β rotation can be ascribed to different origins, closely related to the catechol molecular structure. In fact, the increase of the binding energy upon displacement closely resembles the well-known behavior of the benzene dimer [47,49,50] originated from a “pure” aromatic interaction: shifting one monomer along the R coordinate diminishes the quadrupolar repulsion between the two rings [49], whereas the attractive dispersion interaction decreases to a lesser extent, hence resulting in a global increase of the binding energy [47,49]. As discussed above, the energy profile vs. β rotation is strictly connected with the presence of OH substituents, as suggested by the net increase of the interaction energy in going from a parallel (β = 0°) to an antiparallel (β = 180°) arrangement.
This simple picture is consistent with the minimum of −5.2 kcal/mol (R = 1.2 Å, β = 120°) in a displaced near antiparallel configuration, not coincident with the perfect antiparallel arrangement (β = 180°) where the MP2mod interaction energy is −4.7 kcal/mol. This subtle difference can find a rationale at a closer look of the molecular structure, embracing Wheeler’s idea that unexpected substituents effects can be explained by considering their direct interaction with the neighboring cloud of the other ring [44,45,46]. The β = 120° and β = 180° conformers are displayed in Figure 5. In Figure 5b,d, where a top view of both dimers is shown, the positions of the oxygen atoms are marked with colored circles, to put in evidence the differences between the two arrangements. It appears as in the β = 180° geometry all oxygen atoms lie approximately above a C=C bond of the other ring, resulting in an unfavorable electrostatic interaction with the carbon π orbitals, while at β = 120° only three oxygen atoms contribute to such repulsive term. Consistently, the Hartree–Fock contribution to the total MP2mod energy, which is repulsive in both cases, increases by 1 kcal/mol, in going from β = 120° (3.3 kcal/mol) to β = 180° (4.3 kcal/mol). Finally, another possible source of attractive interaction comes from the HB interaction between the hydrogen atom of one hydroxyl group and the closest oxygen of the other ring, as evidenced in Figure 5a,c, where it appears as in the β = 120° conformer the hydrogen atoms lie at much closer distances (3.7 Å).

3.3. T-Shaped Cathecol Dimers

As shown in Figure 2, another kind of arrangement which can compete with the stacked geometries discussed above is the TS configuration. In this case, most of the interaction energy is expected to come from XH–π forces, in particular when two or one hydroxyl groups point towards the other ring’s plane, as in the TS1 and TS2 geometries.
In order to verify this assumption, the MP2mod computational feasibility has been exploited once again to explore an additional IPES section, related to the TS conformers and shown in Figure 6. At small inter-ring distances, the dependence on β-rotation is striking and the most favorite conformer at R = 4.9 Å is found at β = 0° (i.e., the TS arrangement shown in Figure 6b), with the interaction energy (−2.4 kcal/mol) very similar to the value reported for the benzene dimer in the same configuration [49,50,53,57,71]. Conversely, due to the small distance between the H hydroxyl atom and the other catechol ring (see for instance TS1 in Figure 2), the interaction energy in the 180–300° range is repulsive, with a maximum of almost 25 kcal/mol at β = 270°. The situation changes dramatically by increasing R, as in the 180–300° range the interaction energy shows a much steeper gradient. In fact, the IPES section minimum is found in a TS conformation at β = 270° and R = 5.5 Å, where the hydroxyl group points towards the other ring plane similarly to the TS2 arrangement shown in the right bottom panel of Figure 2, resulting in a total interaction energy of −5.1 kcal/mol.

3.4. Effect of the Hydrogen Bond

The above described competition between stacked and TS geometries misses although another player, which could significantly alter the delicate balance between them. In fact, apart for a small contribution to the stability of the β = 120° conformer in the stacked conformations, the HB contribution was never decisive to the total interaction, due to too large distances between the involved hydrogen and oxygen atoms, which could be reduced by allowing internal rotation around C–O bonds. In order to find even more stable structures, we have released such constraint and performed a full optimization at MP2mod level, starting from four different conformations (see Figure 7, top panels). The first starting geometry is a displaced AFF (AFFD). Next, two TS structures were prepared, with one or both hydroxyl groups pointing down towards the other ring (TSd and TSu, respectively). Notice that the latter is very similar to that taken from crystallographic data and investigated by Estévez et al. [67]. Finally, a fourth arrangement was built from scratch, where the two rings are placed in side-by-side (SS) conformation, with both hydroxyl groups resulting at close distance, thus maximizing the effect of HBs. All optimizations ended up successfully in four different local minima, as confirmed by a frequency calculation purposely carried out for each of the resulting structures. The corresponding optimized structures are shown in the bottom panel of Figure 7 as I, II, III, and IV, respectively.
Dimer formation does not result in large changes in the internal geometry of each catechol monomer. Bond lengths within each monomer change by less than 0.03 Å and the backbone remains planar. For each ring, only one hydroxyl hydrogen moves out of plane, establishing OH–O or OH–π interactions, while the other O–H bond remains nearly coplanar with the ring, due to the formation of an intramolecular OH–O HB with the closest oxygen atom (in Figure 7, geometry II it is out of plane by only 13°). The dihedral angle which drives the position of the out of plane hydrogen is 66° for I, 71.5° for II, and 68.7° for III, whereas the other ring hydrogen does not rotate. The final conformations reveal that the internal rotation has a significant effect on the interplay of the different interaction terms. In fact, as evident from Figure 7, both TS conformations are not stable upon a full optimization, and eventually end up in a stacked arrangement, whereas the AFFD conformer undergoes to the expected rotation from β = 180° to β ≈ 120°, but maintains the stacking arrangements. The OH–O interaction plays the major role in TSd, which becomes II, while OH–π weak HBs guide the hydroxyl rotation and are prevalent in AFFD, which becomes I. Although less stable, the last optimized conformer III, is characterized by a single hydroxyl rotation, which allows the insurgence of a HB (green dashed line in Figure 7), while the other hydrogen remains coplanar to the ring, yet interacting with the other monomer establishing a OH–π noncovalent bond. For the conformation IV, geometry optimization results in a structure which is again less stable than I and II and somewhat more stable than III (see Table 1). In this case, the hydroxyl hydrogens undergo a small rotation with respect to the initial conformation and the two rings are slightly displaced out of the plane that initially (see SS) contained both monomers.
The interaction energies for the four final structures are reported in Table 1, along with the value computed at the same geometries with the MPW1B95/6-311++G(2d,2p) in [67], as well as with the “gold standard” CCSD(T)/CBS. From these data, it is clear that the most stable structure is II, which differs by only 1.6 kcal/mol from I, whereas III and IV are far higher in energy. The agreement between the MP2mod values and their CCSD(T)/CBS counterparts is very good, especially considering that these latter geometries are outside the MP2mod training set, while the computational advantage of using MP2mod with small basis sets is apparent from the last three columns. Surprisingly, the MPW1B95 functional severely underestimates the reference CCSD(T)/CBS interaction energies, yielding, in the present case, only a qualitative correct description, at least according to the protocol provided in [67].
Finally, it is interesting to investigate the different HB contributions in the two most stable conformations I and II. This can be done by performing a rigid scan of the rotation angle δ of the two hydrogen atoms with respect to the C–O bond in both conformations (δ). The results are shown in Figure 8. For δ = 0° (i.e., when each hydrogen is coplanar to the aromatic ring), dispersion interactions are the main source of attraction, although perturbed by the electrostatic interaction between the dipoles, which favors dimer I (in an antiparallel alignment) by ≈1 kcal/mol with respect to dimer II. As δ increases, both the hydrogen atoms involved in the rotation come to closer distances from the other monomer, and may establish HBs. These noncovalent interactions remarkably stabilize both complexes, by almost 7 kcal/mol in I and more than 10 kcal/mol in II. These differences can find a rationale by looking at the insets of Figure 8. In dimer I, since each hydrogen points approximately towards the center of the neighboring ring, two weak HBs of the OH–π type are settled whereas, in dimer II, both hydrogen atoms are involved in a stronger OH–O HB. As a consequence, the minimum of the latter conformer is stabilized by ≈2 kcal/mol with respect to I.

4. Conclusions

In this paper, we have reported our study of the intermolecular landscape of a catechol dimer with a two-fold interest. On the one hand, noncovalent interactions, and especially those involving aromatic rings, govern many biological processes and it is, therefore, of basic importance to reach a good comprehension of the different role that the various forces play in specific systems. On the other hand, noncovalent interactions are still a challenging benchmark for standard computational methods, hence, it can be significant to exploit dedicated approaches.
Catechol is well known to be a precursor of many bioinspired materials and it is, therefore, a good candidate to investigate on the interplay between dispersion interactions, essentially due to aromaticity, and strong (OH–O) or weak (OH–π) HBs, settled by the hydroxyl substituents. The employed MP2mod computational route consists in MP2 calculations with a small 6-31G** basis set, in which the exponents of the polarization functions are suitably modified. This has been done through a validation procedure based on the comparison with the highly accurate CCSD(T)/CBS calculations, resulting in new exponents for polarization functions on carbon (0.27), hydrogen (0.36), and oxygen (0.34).
Within the IPES sections explored, two minima were identified, held together by a network of stacking, OH–O, and OH–π interactions, whose relative weight has been analyzed in some detail. The two catechol units tend to aggregate in stacked conformation, which eventually result more stable than the TS ones, thanks to their ability to establish strong and weak HBs.
A final remark should be made concerning the effects that solvation can have in these systems. Despite most computational approaches designed for noncovalent interactions only focus on two isolated molecules, we are aware that water might affect the results and change the picture that we report here (see, for instance, [27,68]). It is, however, important to have a preliminary reference to guide the more complex study in solution, which is a natural continuation of the one presented here.

Author Contributions

The present article has been conceived by the four authors during a nice dinner. Computations have been shared to be more efficient as it was with the writing. All the authors have therefore equally contributed to the final work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mak, C.H. Unraveling base stacking driving forces in DNA. J. Phys. Chem. B 2016, 120, 6010–6020. [Google Scholar] [CrossRef] [PubMed]
  2. Trouillas, P.; Sancho-García, J.C.; De Freitas, V.; Gierschner, J.; Otyepka, M.; Dangles, O. Stabilizing and modulating color by copigmentation: Insights from theory and experiment. Chem. Rev. 2016, 116, 4937–4982. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Lv, C.; Zhao, G.; Ning, Y. The interactions between plant proteins/enzymes and other food components and their effects on food quality. Crit. Rev. Food Sci. Nutr. 2015, 57, 1718–1728. [Google Scholar] [CrossRef] [PubMed]
  4. Wheeler, S.E.; Seguin, T.J.; Guan, Y.; Doney, A.C. Noncovalent interactions in organocatalysis and the prospect of computational catalyst design. Acc. Chem. Res. 2016, 49, 1061–1069. [Google Scholar] [CrossRef] [PubMed]
  5. Neel, A.J.; Hilton, M.J.; Sigman, M.S.; Toste, F.D. Exploiting non-covalent π interactions for catalyst design. Nature 2017, 543, 637–646. [Google Scholar] [CrossRef] [PubMed]
  6. Ding, X.; Wang, Y. Weak bond-based injectable and stimuli responsive hydrogels for biomedical applications. J. Mater. Chem. B 2017, 5, 887–906. [Google Scholar] [CrossRef]
  7. Yilmazer, N.; Korth, M. Prospects of applying enhanced semi-empirical QM methods for 2101 virtual drug design. Curr. Med. Chem. 2016, 23, 2101–2111. [Google Scholar] [CrossRef] [PubMed]
  8. Meanwell, N.A. A synopsis of the properties and applications of heteroaromatic rings in medicinal chemistry. Adv. Heterocycl. Chem. 2017, 123, 245–361. [Google Scholar]
  9. Ghosh, T.; Panicker, J.; Nair, V. Self-assembled organic materials for photovoltaic application. Polymers 2017, 9, 112. [Google Scholar] [CrossRef]
  10. Würthner, F. Dipole–dipole interaction driven self-assembly of merocyanine dyes: From dimers to nanoscale objects and supramolecular materials. Acc. Chem. Res. 2016, 49, 868–876. [Google Scholar] [CrossRef] [PubMed]
  11. Chandra, B.K.C.; D’Souza, F. Design and photochemical study of supramolecular donor–acceptor systems assembled via metal–ligand axial coordination. Coord. Chem. Rev. 2016, 322, 104–141. [Google Scholar]
  12. Bredas, J.-L.; Marder, S.R.; Reichmanis, E. Preface to the chemistry of materials special issue on π-functional materials. Chem. Mater. 2011, 23, 309. [Google Scholar] [CrossRef]
  13. Shewmon, N.T.; Watkins, D.L.; Galindo, J.F.; Zerdan, R.B.; Chen, J.; Keum, J.; Roitberg, A.E.; Xue, J.; Castellano, R.K. Enhancement in organic photovoltaic efficiency through the synergistic interplay of molecular donor hydrogen bonding and π-stacking. Adv. Funct. Mater. 2015, 25, 5166–5177. [Google Scholar] [CrossRef]
  14. Rest, C.; Kandanelli, R.; Fernández, G. Strategies to create hierarchical self-assembled structures via cooperative non-covalent interactions. Chem. Soc. Rev. 2015, 44, 2543–2572. [Google Scholar] [CrossRef] [PubMed]
  15. Hammann, F.; Schmid, M. Determination and quantification of molecular interactions in protein films: A review. Materials 2014, 7, 7975–7996. [Google Scholar] [CrossRef] [PubMed]
  16. Makwana, K.M.; Mahalakshmi, R. Implications of aromatic-aromatic interactions: From protein structures to peptide models. Protein Sci. 2015, 24, 1920–1933. [Google Scholar] [CrossRef] [PubMed]
  17. Cragg, P.J. Supramolecular Chemistry—From Biological Inspiration to Biomedical Applications; Springer: Amsterdam, The Netherlands, 2010. [Google Scholar]
  18. Wun Hwang, J.; Li, P.; Shimizu, K.D. Synergy between experimental and computational studies of aromatic stacking interactions. Org. Biomol. Chem. 2017, 15, 1554–1564. [Google Scholar] [CrossRef] [PubMed]
  19. Sherrill, C.D. Energy component analysis of π interactions. Acc. Chem. Res. 2013, 46, 1020–1028. [Google Scholar] [CrossRef] [PubMed]
  20. Heinzmann, C.; Weder, C.; de Espinosa, L.M. Supramolecular polymer adhesives: Advanced materials inspired by nature. Chem. Soc. Rev. 2015, 342, 342–358. [Google Scholar] [CrossRef] [PubMed]
  21. Jenkins, C.L.; Siebert, H.M.; Wilker, J.J. Integrating mussel chemistry into a bio-based polymer to create degradable adhesives. Macromolecules 2017, 50, 561–568. [Google Scholar] [CrossRef]
  22. Kord Forooshani, P.; Lee, B.P. Recent approaches in designing bioadhesive materials inspired by mussel adhesive protein. J. Polym. Sci. Part A: Polym. Chem. 2017, 55, 9–33. [Google Scholar] [CrossRef] [PubMed]
  23. North, M.A.; Del Grosso, C.A.; Wilker, J.J. High strength underwater bonding with polymer mimics of mussel adhesive proteins. ACS Appl. Mater. Interfaces 2017, 9, 7866–7872. [Google Scholar] [CrossRef] [PubMed]
  24. Saiz-Poseu, J.; Sedó, J.; García, B.; Benaiges, C.; Parella, T.; Alibés, R.; Hernando, J.; Busqué, F.; Ruiz-Molina, D. Versatile nanostructured materials via direct reaction of functionalized catechols. Adv. Mater. 2013, 25, 2066–2070. [Google Scholar] [CrossRef] [PubMed]
  25. Sedó, J.; Saiz-Poseu, J.; Busqué, F.; Ruiz-Molina, D. Catechol-based biomimetic functional materials. Adv. Mater. 2013, 25, 653–701. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, X.; Jing, S.; Liu, Y.; Liu, S.; Tan, Y. Diblock copolymer containing bioinspired borneol and dopamine moieties: Synthesis and antibacterial coating applications. Polymer 2017, 116, 314–323. [Google Scholar] [CrossRef]
  27. Barone, V.; Cacelli, I.; Crescenzi, O.; D’Ischia, M.; Ferretti, A.; Prampolini, G.; Villani, G. Unraveling the interplay of different contributions to the stability of the quinhydrone dimer. RSC Adv. 2014, 4, 876–885. [Google Scholar] [CrossRef]
  28. Burns, L.A.; Vázquez-Mayagoitia, A.; Sumpter, B.G.; Sherrill, C.D. Density-functional approaches to noncovalent interactions: A comparison of dispersion corrections (DFT-D), exchange-hole dipole moment (XDM) theory, and specialized functionals. J. Chem. Phys. 2011, 134, 084107. [Google Scholar] [CrossRef] [PubMed]
  29. Corminboeuf, C. Minimizing density functional failures for non-covalent interactions beyond van der Waals complexes. Acc. Chem. Res. 2014, 47, 3217–3224. [Google Scholar] [CrossRef] [PubMed]
  30. Goerigk, L.; Grimme, S. Efficient and accurate double-hybrid-meta-GGA density functionals-evaluation with the extended GMTKN30 database for general main group thermochemistry, kinetics, and noncovalent interactions. J. Chem. Theory Comput. 2011, 7, 291–309. [Google Scholar] [CrossRef] [PubMed]
  31. Goldey, M.B.; Belzunces, B.; Head-Gordon, M. Attenuated MP2 with a long-range dispersion correction for treating nonbonded interactions. J. Chem. Theory Comput. 2015, 11, 4159–4168. [Google Scholar] [CrossRef] [PubMed]
  32. Grimme, S. Density functional theory with London dispersion corrections. WIREsComput. Mol. Sci. 2011, 1, 211–228. [Google Scholar] [CrossRef]
  33. Grimme, S.; Hansen, A.; Brandenburg, J.G.; Bannwarth, C. Dispersion-corrected mean-field electronic structure methods. Chem. Rev. 2016, 116, 5105–5154. [Google Scholar] [CrossRef] [PubMed]
  34. Hobza, P.; Zahradník, R.; Müller-Dethlefs, K. The world of non-covalent interactions: 2006. Collect. Czechoslov. Chem. Commun. 2006, 71, 443–531. [Google Scholar] [CrossRef]
  35. Piton, M. Accurate intermolecular interaction energies from a combination of MP2 and TDDFT response theory. J. Chem. Theory Comput. 2010, 6, 168–178. [Google Scholar] [CrossRef] [PubMed]
  36. Prampolini, G.; Cacelli, I.; Ferretti, A. Intermolecular interactions in eumelanins: A computational bottom-up approach. I. Small building blocks. RSC Adv. 2015, 5, 38513–38526. [Google Scholar] [CrossRef]
  37. Řezáč, J.; Hobza, P. Describing noncovalent interactions beyond the common approximations: How accurate is the “gold standard”, CCSD(T) at the complete basis set limit? J. Chem. Theory Comput. 2013, 9, 2151–2155. [Google Scholar] [CrossRef] [PubMed]
  38. Řezáč, J.; Hobza, P. Benchmark calculations of interaction energies in noncovalent complexes and their applications. Chem. Rev. 2016, 116, 5038–5071. [Google Scholar] [CrossRef] [PubMed]
  39. Richard, R.M.; Lao, K.U.; Herbert, J.M. Achieving the CCSD(T) basis-set limit in sizable molecular clusters: Counterpoise corrections for the many-body expansion. J. Phys. Chem. Lett. 2013, 4, 2674–2680. [Google Scholar] [CrossRef] [PubMed]
  40. Riley, K.E.; Platts, J.A.; Řezáč, J.; Hobza, P.; Hill, J.G. Assessment of the performance of MP2 and MP2 variants for the treatment of noncovalent interactions. J. Phys. Chem. A 2012, 116, 4159–4169. [Google Scholar] [CrossRef] [PubMed]
  41. Sherrill, C.D.; Takatani, T.; Hohenstein, E.G. An assessment of theoretical methods for nonbonded interactions: Comparison to complete basis set limit coupled-cluster potential energy curves for the benzene dimer, the methane dimer, benzene−methane, and benzene−H2S. J. Phys. Chem. A 2009, 113, 10146–10159. [Google Scholar] [CrossRef] [PubMed]
  42. Tkatchenko, A.; DiStasio, R.A.; Head-Gordon, M.; Scheffler, M. Dispersion-corrected Møller–Plesset second-order perturbation theory. J. Chem. Phys. 2009, 131, 94106. [Google Scholar] [CrossRef] [PubMed]
  43. Zhao, Y.; Truhlar, D.G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four M06-class functionals and 12 other function. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar] [CrossRef]
  44. An, Y.; Doney, A.C.; Andrade, R.B.; Wheeler, S.E. Stacking interactions between 9-methyladenine and heterocycles commonly found in pharmaceuticals. J. Chem. Inf. Model. 2016, 56, 906–914. [Google Scholar] [CrossRef] [PubMed]
  45. Wheeler, S.E.; Bloom, J.W.G. Toward a more complete understanding of noncovalent interactions involving aromatic rings. J. Phys. Chem. A 2014, 118, 6133–6147. [Google Scholar] [CrossRef] [PubMed]
  46. Wheeler, S.E.; Houk, K.N. Substituent effects in the benzene dimer are due to direct interactions of the substituents with the unsubstituted benzene. J. Am. Chem. Soc. 2008, 130, 10854–10855. [Google Scholar] [CrossRef] [PubMed]
  47. Podeszwa, R.; Bukowski, R.; Szalewicz, K. Potential energy surface for the benzene dimer and perturbational analysis of π−π interactions. J. Phys. Chem. A 2006, 110, 10345–10354. [Google Scholar] [CrossRef] [PubMed]
  48. 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]
  49. Tsuzuki, S.; Honda, K.; Uchimaru, T.; Mikami, M.; Tanabe, K. Origin of attraction and directionality of the π/π interaction: Model chemistry calculations of benzene dimer interaction. J. Am. Chem. Soc. 2002, 124, 104–112. [Google Scholar] [CrossRef] [PubMed]
  50. 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]
  51. Kroon-Batenburg, L.; Van Duijneveldt, F. The use of a moment-optimized DZP basis set for describing the interaction in the water dimer. J. Mol. Struct.: THEOCHEM 1985, 22, 185–199. [Google Scholar] [CrossRef]
  52. Hobza, P.; Zahradnik, R. Intermolecular interactions between medium-sized systems. Nonempirical and empirical calculations of interaction energies. Successes and failures. Chem. Rev. 1988, 88, 871–897. [Google Scholar] [CrossRef]
  53. 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]
  54. Sponer, J.; Leszczynski, J.; Hobza, P. Base stacking in cytosine dimer. A comparison of correlated ab initio calculations with three empirical potential models and density functional theory calculations. J. Comput. Chem. 1996, 17, 841–850. [Google Scholar] [CrossRef]
  55. Šponer, J.; Leszczyński, J.; Hobza, P. Nature of nucleic acid−base stacking: Nonempirical ab initio and empirical potential characterization of 10 stacked base dimers. Comparison of stacked and H-bonded base pairs. J. Phys. Chem. 1996, 100, 5590–5596. [Google Scholar] [CrossRef]
  56. Hobza, P.; Šponer, J. Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations. J. Am. Chem. Soc. 2002, 124, 11802–11808. [Google Scholar] [CrossRef] [PubMed]
  57. Cacelli, I.; Cinacchi, G.; Prampolini, G.; Tani, A. Computer simulation of solid and liquid benzene with an atomistic interaction potential derived from ab initio calculations. J. Am. Chem. Soc. 2004, 126, 14278–14286. [Google Scholar] [CrossRef] [PubMed]
  58. Mignon, P.; Loverix, S.; De Proft, F.; Geerlings, P. Influence of stacking on hydrogen bonding: Quantum chemical study on pyridine–benzene model complexes. J. Phys.Chem. A 2004, 108, 6038–6044. [Google Scholar] [CrossRef]
  59. Řeha, D.; Kabeláč, M.; Ryjáček, F.; Šponer, J.; Šponer, J.E.; Elstner, M.; Suhai, S.; Hobza, P. Intercalators. 1. Nature of stacking interactions between intercalators (ethidium, daunomycin, ellipticine, and 4′,6-diaminide-2-phenylindole) and DNA base pairs. Ab initio quantum chemical, density functional theory, and empirical potential study. J. Am. Chem. Soc. 2002, 124, 3366–3376. [Google Scholar] [CrossRef] [PubMed]
  60. Sponer, J.; Gabb, H.A.; Leszczynski, J.; Hobza, P. Base–base and deoxyribose–base stacking interactions in B-DNA and Z-DNA: A quantum-chemical study. Biophys. J. 1997, 73, 76–87. [Google Scholar] [CrossRef]
  61. Hobza, P.; Kabeláč, M.; Šponer, J.; Mejzlík, P.; Vondrášek, J. Performance of empirical potentials (AMBER, CFF95, CVFF, CHARMM, OPLS, POLTEV), semiempirical quantum chemical methods (AM1, MNDO/M, PM3), andab initio Hartree–Fock method for interaction of DNA bases: Comparison with nonempirical beyond Hartree–Fock results. J. Comput. Chem. 1997, 18, 1136–1150. [Google Scholar]
  62. Cacelli, I.; Prampolini, G.; Tani, A. Atomistic simulation of a nematogen using a force field derived from quantum chemical calculations. J. Phys. Chem. B 2005, 109, 3531–3538. [Google Scholar] [CrossRef] [PubMed]
  63. Cacelli, I.; Lami, C.F.; Prampolini, G. Force-field modeling through quantum mechanical calculations: Molecular dynamics simulations of a nematogenic molecule in its condensed phases. J. Comput. Chem. 2009, 30, 366–378. [Google Scholar] [CrossRef] [PubMed]
  64. Cacelli, I.; Cimoli, A.; Livotto, P.R.; Prampolini, G. An automated approach for the parameterization of accurate intermolecular force-fields: Pyridine as a case study. J. Comput. Chem. 2012, 33, 1055–1067. [Google Scholar] [CrossRef] [PubMed]
  65. Micillo, R.; Panzella, L.; Iacomino, M.; Prampolini, G.; Cacelli, I.; Ferretti, A.; Crescenzi, O.; Koike, K.; Napolitano, A.; D’Ischia, M. Eumelanin broadband absorption develops from aggregation-modulated chromophore interactions under structural and redox control. Sci. Rep. 2017, 7, 41532. [Google Scholar] [CrossRef] [PubMed]
  66. Prampolini, G.; Greff da Silveira, L.; Jacobs, M.; Livotto, P.R.; Cacelli, I. Interaction energy landscapes of aromatic heterocycles through a reliable yet affordable computational approach. J. Chem. Theory Comput. 2017. submitted. [Google Scholar]
  67. Estévez, L.; Otero, N.S.; Mosquera, R.A. Computational study on the stacking interaction in catechol complexes. J. Phys. Chem. A 2009, 113, 11051–11058. [Google Scholar] [CrossRef] [PubMed]
  68. Barone, V.; Cacelli, I.; Ferretti, A.; Prampolini, G.; Villani, G. Proton and electron transfer mechanisms in the formation of neutral and charged quinhydrone-like complexes: A multilayered computational study. J. Chem.Theory Comput. 2014, 10, 4883–4895. [Google Scholar] [CrossRef] [PubMed]
  69. Boys, S.F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  70. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  71. Prampolini, G.; Livotto, P.R.; Cacelli, I. Accuracy of quantum mechanically derived force-fields parameterized from dispersion-corrected DFT data: The benzene dimer as a prototype for aromatic interactions. J. Chem. Theory Comput. 2015, 11, 5182–5196. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) Catechol structural formula (top) and graphical representation (bottom). Stacked dimers: (b) face-to-face (FF) and (c) antiparallel face-to-face (AFF); T-shaped (TS) dimers: (d) TS1 and (e) TS2. C: Cyan; H: White; O: Red.
Figure 1. (a) Catechol structural formula (top) and graphical representation (bottom). Stacked dimers: (b) face-to-face (FF) and (c) antiparallel face-to-face (AFF); T-shaped (TS) dimers: (d) TS1 and (e) TS2. C: Cyan; H: White; O: Red.
Biomimetics 02 00018 g001
Figure 2. Comparison between the ‘best exponent’ and CCSD(T)/CBS for the interaction energy profiles obtained by displacement of the four structures shown in Figure 1.
Figure 2. Comparison between the ‘best exponent’ and CCSD(T)/CBS for the interaction energy profiles obtained by displacement of the four structures shown in Figure 1.
Biomimetics 02 00018 g002
Figure 3. MP2mod results for the stacked configurations. (a) Interaction energy as a function of the inter-ring separation R for different β angles. (b) Interaction energy as a function of the β angle at the ring–ring separation (R = 3.5 Å) corresponding to the minimum energy.
Figure 3. MP2mod results for the stacked configurations. (a) Interaction energy as a function of the inter-ring separation R for different β angles. (b) Interaction energy as a function of the β angle at the ring–ring separation (R = 3.5 Å) corresponding to the minimum energy.
Biomimetics 02 00018 g003
Figure 4. (a) Two-dimensional scan of the catechol interaction potential energy surface (IPES) in stacked conformations, performed at the MP2mod level. (b) The IPES section was sampled by varying the angle β; the displacement R is also shown.
Figure 4. (a) Two-dimensional scan of the catechol interaction potential energy surface (IPES) in stacked conformations, performed at the MP2mod level. (b) The IPES section was sampled by varying the angle β; the displacement R is also shown.
Biomimetics 02 00018 g004
Figure 5. Stacked displaced (R = 1.2 Å) geometries at (a,b) β = 120° and (c,d) β = 180°. (a,c): Side view, H–O distances of (a) 3.7 and (c) 5.0 Å are indicated with a green arrow ; (b,d): Top view, the position of oxygen atoms is shown with colored circles, distinguishing more (yellow) or less (green) interacting ones.
Figure 5. Stacked displaced (R = 1.2 Å) geometries at (a,b) β = 120° and (c,d) β = 180°. (a,c): Side view, H–O distances of (a) 3.7 and (c) 5.0 Å are indicated with a green arrow ; (b,d): Top view, the position of oxygen atoms is shown with colored circles, distinguishing more (yellow) or less (green) interacting ones.
Biomimetics 02 00018 g005
Figure 6. (a) Two-dimensional scan of the catechol IPES in TS conformation, performed at MP2mod level. The IPES section was sampled by varying the angle β and the displacement R shown in (b), where the TS arrangement at β = 0° is displayed. The β rotation is performed as indicated by the black arrow (e.g., for β = 240° the dimer is found in the TS1 geometry shown in the right top panel of Figure 2).
Figure 6. (a) Two-dimensional scan of the catechol IPES in TS conformation, performed at MP2mod level. The IPES section was sampled by varying the angle β and the displacement R shown in (b), where the TS arrangement at β = 0° is displayed. The β rotation is performed as indicated by the black arrow (e.g., for β = 240° the dimer is found in the TS1 geometry shown in the right top panel of Figure 2).
Biomimetics 02 00018 g006
Figure 7. MP2mod geometry optimization starting from the displaced AFF (AFFD), TSd, and TSu conformations (top panels). The corresponding optimized structures, I, II, III, and IV, are displayed in the bottom row. The rotated hydroxyl groups are evidenced in the top panel with a blue arrow, while the atoms involved in OH–O and OH–π interactions are connected in the bottom panels by green and orange dashed lines, respectively.
Figure 7. MP2mod geometry optimization starting from the displaced AFF (AFFD), TSd, and TSu conformations (top panels). The corresponding optimized structures, I, II, III, and IV, are displayed in the bottom row. The rotated hydroxyl groups are evidenced in the top panel with a blue arrow, while the atoms involved in OH–O and OH–π interactions are connected in the bottom panels by green and orange dashed lines, respectively.
Biomimetics 02 00018 g007
Figure 8. MP2mod scans of the HOCC dihedral (δ) for conformers I (blue) and II (red), highlighting the role of OH–O and OH–π interactions, indicated with orange and green dashed lines, respectively.
Figure 8. MP2mod scans of the HOCC dihedral (δ) for conformers I (blue) and II (red), highlighting the role of OH–O and OH–π interactions, indicated with orange and green dashed lines, respectively.
Biomimetics 02 00018 g008
Table 1. Interaction energies, in kcal/mol, for the four optimized conformations shown in Figure 7, computed with MP2mod, CCSD(T)/CBS and MPW1B95/6-311++G(2d,2p). Central processing unit (CPU) times on a single 2.60 GHz Intel® Xeon CPU are also given for an evaluation of the computational cost of the different methods.
Table 1. Interaction energies, in kcal/mol, for the four optimized conformations shown in Figure 7, computed with MP2mod, CCSD(T)/CBS and MPW1B95/6-311++G(2d,2p). Central processing unit (CPU) times on a single 2.60 GHz Intel® Xeon CPU are also given for an evaluation of the computational cost of the different methods.
GeometryEnergies (kcal/mol)CPU Time (min)
MP2modCCSD(T)/CBSMPW1B95MP2modCCSD(T)/CBSMPW1B95
I−10.7−11.1−8.1 2750,640145
II−12.4−12.6−8.32549,740180
III−5.3−5.7−2.22750,82079
IV−6.1−7.3−5.71851,720142

Share and Cite

MDPI and ACS Style

Barone, V.; Cacelli, I.; Ferretti, A.; Prampolini, G. Noncovalent Interactions in the Catechol Dimer. Biomimetics 2017, 2, 18. https://doi.org/10.3390/biomimetics2030018

AMA Style

Barone V, Cacelli I, Ferretti A, Prampolini G. Noncovalent Interactions in the Catechol Dimer. Biomimetics. 2017; 2(3):18. https://doi.org/10.3390/biomimetics2030018

Chicago/Turabian Style

Barone, Vincenzo, Ivo Cacelli, Alessandro Ferretti, and Giacomo Prampolini. 2017. "Noncovalent Interactions in the Catechol Dimer" Biomimetics 2, no. 3: 18. https://doi.org/10.3390/biomimetics2030018

APA Style

Barone, V., Cacelli, I., Ferretti, A., & Prampolini, G. (2017). Noncovalent Interactions in the Catechol Dimer. Biomimetics, 2(3), 18. https://doi.org/10.3390/biomimetics2030018

Article Metrics

Back to TopTop