Next Article in Journal
Autophagy Modulation as a Treatment of Amyloid Diseases
Previous Article in Journal
Investigation of Metal-Organic Framework-5 (MOF-5) as an Antitumor Drug Oridonin Sustained Release Carrier
Previous Article in Special Issue
Tetrel Interactions from an Interacting Quantum Atoms Perspective
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 24
Issue 18
10.3390/molecules24183370
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intermolecular Non-Covalent Carbon-Bonding Interactions with Methyl Groups: A CSD, PDB and DFT Study

by
Tiddo J. Mooibroek
van ’t Hoff Institute for Molecular Sciences, Universiteit van Amsterdam, Science Park 904, 1098 XH, Amsterdam, The Netherlands
Molecules 2019, 24(18), 3370; https://doi.org/10.3390/molecules24183370
Submission received: 28 July 2019 / Revised: 5 September 2019 / Accepted: 12 September 2019 / Published: 16 September 2019
(This article belongs to the Special Issue Tetrel Bonds)
Browse Figures
Review Reports Versions Notes

Abstract

:
A systematic evaluation of the CSD and the PDB in conjunction with DFT calculations reveal that non-covalent Carbon-bonding interactions with X–CH3 can be weakly directional in the solid state (P ≤ 1.5) when X = N or O. This is comparable to very weak CH hydrogen bonding interactions and is in line with the weak interaction energies calculated (≤ –1.5 kcal·mol−1) of typical charge neutral adducts such as [Me3N-CH3···OH2] (2a). The interaction energy is enhanced to ≤–5 kcal·mol−1 when X is more electron withdrawing such as in [O2N-CH3··O=Cdme] (20b) and to ≤18 kcal·mol−1 in cationic species like [Me3O+-CH3···OH2]+ (8a).

Graphical Abstract

1. Introduction

The manner in which molecules interact with one another is largely determined by non-covalent interactions.c [1] So-called ‘σ-Hole interactions’c [2,3,4,5] like hydrogen bonding are prominent identifiable interactions that bear biological significance [6]. Such σ-hole interactions have also been identified with other non-metals [7,8,9,10,11] like halogen atoms to generate halogen bonding interaction [12,13]. The impact of halogen bonding interactions on molecular biology has come into focus since about 2004 [14]. Indeed, evaluations of the protein data bank (PDB) [15] revealed that halogen bonding is structurally very similar to hydrogen bonding [12,14,16,17,18] and can be functionally relevant [19,20,21,22]. Relatively weak π-hole interactions [4,23,24,25,26,27,28,29,30,31] involving organic carbonyls, [26,32,33,34,35,36] π-acidic aromatics, [37,38] metal carbonyls [33,34,36,39] and nitro-compounds [40,41,42,43,44,45] are increasingly acknowledged as relevant drivers of molecular aggregation such as in ligand-protein complexes.
The impact of a novel type of weak interaction on molecular recognition phenomena naturally leads one to speculate that other non-canonical interactions may play a similar role. One interesting candidate are σ-hole interactions involving sp3-hybridized C-atoms. Such interactions have been studied since about 2013 [7,46] and are particularly interesting because sp3-C is abundant in living systems. More specifically, the methyl group (X–CH3, where X = any atom or group) is frequently encountered in natural and synthetic compounds and ‘non-covalent Carbon bonding’ involving methyl groups has thus been studied by various researchers [47,48,49,50,51,52,53,54,55,56,57,58,59]. Most of these contributions are computational inquiries, while a small amount of these articles also deals with an analysis of non-covalent Carbon bonding interactions in protein structures present in the Protein Data Bank (PDB) [47,50,56]. Interestingly, none of the studies so far have systematically evaluated the crystal structure data present in the Cambridge Structure Database (CSD) [60,61]. What is more, evaluations of the PDB were largely anecdotal or only considered structures that comply to the (rather strict) geometric criteria of a Carbon bonding geometry. Some also included intramolecular contact distances (which are notoriously difficult to evaluate).
In this contribution a combined CSD and PDB evaluation is presented aimed at elucidating whether electron rich entities have a preferential orientation around a methyl group within a rather large envelope, i.e., whether intermolecular non-covalent Carbon bonding interactions with methyl groups are directional. For evaluative purposes, several Density Functional Theory (DFT) computations were conducted as well. This combined database/DFT study reaffirms that non-covalent Carbon-bonding interactions with X-CH3 can be significant, although the interaction is hardly directional, in particular when the methyl group is poorly polarized such as most C–CH3 structures.

2. Materials and Methods

2.1. General Information on Database Analyses

The CSD [60,61] version 5.40 including two updates (until May 2019) was inspected using ConQuest [62] version 2.0.2 (build 246353, 2019). X-ray powder structures were omitted from the searches, which were further limited to structures containing 3D coordinates and those with an R-factor ≤ 0.1. The PDB was queried using Relibase [63] 3.2.3 and restricted to protein and DNA crystal structures where the packing environment was also searched. No other restrictions were imposed on the PDB search. Datasets were obtained using the general query shown in Figure 1a. The methyl groups were split in those connected to a C, N, O, P, or S atom (X in the figure, in the PDB search specified as part of a ligand). The interacting ‘electron rich’ partners (ElR in the figure) considered were a water, amide or carboxy-O atom, a sulphur atom or the centroid of an aryl ring (in the PDB search always specified as part of the protein). The geometric constraints imposed on the searches were that the intermolecular distance d between the methyl C-atom and ElR was ≤5 Å and that the X–CH3∙∙∙ElR angle (α) was 90°–180°. All the data were thus confined within a hemisphere with a basal radius of 5 Å, centered on the methyl C-atoms as is shown in Figure 1b.

2.2. Methodology to Generate P(α) Plots

The datasets obtained as described above (2.1) were analysed to assess whether the distribution of ElR within the methyl-centered hemisphere reflects any directionality. This method has been successfully applied to assess the directional behaviors (in the solid state) of various other weak non-covalent interactions such as anion/lone-pair-π, [29,64] CH-π, [11,65] halogen-π [66,67] and nitro π-hole interactions [42,68]. The method works by first computing the freely accessible volumes at each α-value (αfree) by subtracting the volume of a model methyl group from a spherical cone with 5 Å height and a cone angle of 180-α. This can be achieved by using the 3D-drawing program Autodesk® Inventor® Pro [29]. This is illustrated in Figure 1b, where the spherical cones are shown at 10° intervals. The model methyl group was generated by using standard aliphatic C–H bond distances (1.06 Å) [69] and the van der Waals radius of C (1.70 Å) and H (1.09 Å) [70]. The interfering volume between each spherical cone and the model methyl group can be obtained using the ‘inspect interference’ option in Autodesk® Inventor® Pro; the red part in Figure 1b is the interfering volume involving a spherical cone with a cone angle of 60° (i.e., at α = 120°). The volume differences between such ‘free’ volumes with increasing values of α thus give the absolute volume distribution of freely accessible volume around a methyl model within the hemisphere, as a function of α: Δαfree(α). Dividing each volume (Δαfree) in this distribution by the total freely accessible volume (i.e., the volume of a hemisphere minus the interfering volume of the model methyl group in that hemisphere) thus gives the relative volume distribution as a function of α: Δrelαfree(α). This distribution is the random (or volume) distribution. The data retrieved form the CSD and the PDB can be binned as a function of α. Relating this binned data to all the data in a dataset thus gives the observed relative distribution as a function of α: Δrelαdata(α). The quotient of this relative data distribution over the random distribution is a measure for the actual probability (P) of finding data at a certain value of α. That is, P(α) is unity for a random distribution of data, while P-values larger than unity reveal a relative concentration of data, which is indicative of attractive interactions.

2.3. Methodology to Generate N(d’) Plots

A second analysis involved plotting the hit fraction (N) for a subset with α = 160°–180° as a function of the van der Waals corrected H3C∙∙∙ElR distance d’ (i.e., d – vdW(C) – vdW(ElR)): [70] N(d’). Such distributions show how much of the data is involved in van der Waals overlap with the methyl C-atom along the vector of the X–CH3 bond and how such data is distributed. For attractive interacting pairs this distribution is expected to exhibit a peak-like feature, or an S-like curvature when the cumulative hit fraction is used.

2.4. Computational Methods

DFT geometry optimization calculations were performed with Spartan 2016 at the B3LYP [71,72]-D3 [73]/def2-TZVP [74,75] level of theory, which is known to give accurate results at reasonable computational cost and a very low basis set superposition error (BSSE) [74,75]. The typical starting geometry for possible Carbon bonding adducts was set to d’ = −0.1 Å and α = 180°, and in the case of dimethylacetamide the C···O=C angle was also set to 180°. The geometry optimizations were performed without any constraints. For other geometries (e.g., a H-bonded geometry), the molecular fragments were manually oriented in a suitable constellation before starting an unconstrained geometry optimization. The Amsterdam Density Functional (ADF) [76] modelling suite at the B3LYP [71,72]-D3 [73]/TZ2P [74,75] level of theory (no frozen cores) was used for energy decomposition and ‘atoms in molecules’ [77] analyses. Details of the Morokuma-Ziegler inspired energy decomposition scheme used in the ADF-suite have been reported elsewhere [76,78] and the scheme has proven useful to evaluate hydrogen bonding interactions [79].

3. Results and Discussion

3.1. P(α) Plots

A numerical overview of the amount of crystallographic information files (CIFs) and protein data bank files (PDBs) for each search query is given in Table S1, together with the amount of hits found in each dataset (a .cif or .pdb file can contain multiple hits). Shown in Figure 2 are the P(α) plots for the CSD (left) and PDB (right) data plotted at 5° intervals involving X = C, N, O and ElR = water-O (top) or amide-O (bottom). These datasets were chosen because they all contained a large number of hits (>7,500) and thus allow for the most reliable comparison. A complete set of P(α) plots is provided in Figure S1.
The data plotted in Figure 1 largely trace the line at P = 1 (highlighted in green), which is indicative of a random distribution of data. For X = O and N, these values are somewhat above unity around α = 160°–180° for water O-atoms in both databases and for amide O-atoms in the CSD. The maximum P-values are very small at about 1.5, which indicates a very small amount of directionality. Indeed, maximum P-values for several weak inter-molecular interactions are: ~2.5 for CH-π; [11,65] ~3 for π interactions with nitro compounds; [42,68] ~2.5–5 for anion-π and lone-pair-π; [29,64] ~2.5–10 for halogen-π [67] and also about 2.5–10 for halogen bonding with aryl-halogens. [66] Interestingly, the P-values did not peak near α = 90°–120°, an angle congruent with hydrogen bonding. These data thus suggest that the Carbon binding geometry is more directional than a hydrogen bonding geometry, although this directionality is very weak. In all cases for X = C, the P-values around α = 160°–180° are below unity, suggesting that the Carbon bonding geometry is least favored in these instances. The data for ElR = RCO2 and RyCS are very similarly distributed and the data for aryl rings is skewed towards α = 160°–180° only for the CSD data (for X = C, N, O, P and S, see Figure S1). For all other cases where X = P, very few hits were obtained (most numerous was RCO2 in the CSD with N = 1,454). While some datasets with X = S were of a reasonable size, too many were well below N = 7,500 and these data will thus not be discussed in the main text (there is a small discussion in the caption of Figure S1).

3.2. N(d’) plots

The data characterized by α = 160°–180° was inspected further by means of N(d’) plots as described in the methods and materials section. A numerical overview of the amount of data in each dataset as well as the (relative) amount of van der Waals overlap found in each dataset is given in Table S2. Shown in Figure 3 are plots of the hit fractions (in %) as a function of d’ (the van der Waals corrected d) for X = C, N, O and ElR = water or amide. The same data plotted as cumulative hit fractions is shown in Figure S2 and the N(d’) plots for all datasets containing >500 hits are shown in Figure S3.
In the CSD (left), The data involving N/O–CH3∙∙∙Oamides are very similarly distributed and grouped near d’ ≈ 0 Å with about 30% of all the data involved in van der Waals overlap. Likewise, data involving N/O–CH3∙∙∙Owaters are also very similarly distributed but group near the larger d’ ≈ 0.25 Å with 15% van der Waals overlap. The C–CH3∙∙∙Oamides/waters data hardly displays van der Waals overlap (3–4%) and is broadly grouped around d’ = 0.4 Å for waters and not grouped at all for amides. Similar trends are present in the PDB (right), albeit the features are much less pronounced and the N/O–CH3∙∙∙O datasets with amides and waters are very similar.
The data presented in Figure 3 thus imply a somewhat directional nature of X–CH3∙∙∙Oamides/waters interactions for X = N and O, but not at all for X = C. These findings are in line with the lack of directionality observed in the P(α) plots for X = C and the somewhat directional behavior for X = N or O (see Figure 2). A likely explanation for this is the larger (Pauling) electronegativity of N (3.04) and O (3.44) compared to C (2.55), resulting in a larger degree of polarization of the X–CH3 bond for X = N or O. Another conceivable manner to make a methyl group pore electropositive is to bind it to a cationic fragment such as in protonated or quaternary R3N+–CH3 fragments. Thus, an additional dataset was retrieved from the CSD involving R3N+–CH3···ElR pairs that fulfilled the d ≤ 5 Å and α = 160°–180° criteria (see bottom entries in Tables S1 and S2). Shown in Figure 4 are the N(d’) plots of the most numerous datasets involving R3N+–CH3 (hexagonals, N+) together with similar datasets involving all possible N–CH3 fragments (diamonds, N). In all three cases (ElR = water, carboxy or aryl), the distribution is shifted to lower d’ distances for cationic R3N+–CH3, which means that relatively more van der Waals overlap is present in these datasets. As can be seen especially from the cumulative N(d’) plots (left), the grouping is tightest with carboxy O-atoms (green), followed by water O-atoms (red) and aryl rings centroids (grey) are not grouped at all (nearly linear).

3.3. Computations

In order to gain insight into the nature and energetics of possible non-covalent Carbon bonding interactions involving various methyl groups and water or amide O-atoms, DFT calculations were performed of X–CH3 adducts with water and with dimethylacetamide (dma, see methods section for details). An overview of these adducts is given in Table 1, together with α of the optimized structures, the total interaction energy of the adducts in kcal·mol−1 and the percentages of electrostatic (E), orbital (O), and dispersion (D) interactions that contribute to this total energy [79]. Perspective views and atoms-in-molecules analyses of all converged structures are shown in Figure S4 and several representative examples are shown in Figure 5.
For comparison purposes, adducts with ethane were computed as shown in entries 1 of Table 1. Both structures converged to a hydrogen bonding interaction. ΔE = −1.2 kcal·mol−1 in 1a and a methyl acts as electron donating site; i.e., an O–H···C hydrogen bonding interaction. This can be understood due to the polarization in ethane, where both C’s are most electronegative and are polarized by the H-atoms. Adduct 1b is about twice as stable with ΔE = −2.8 kcal·mol−1 and also features a hydrogen bonding interaction, but now between a methyl CH and a π-bond in dma. Both 1a and 1b are stabilized mainly by dispersion (46–58%), then electrostatics (32–26%) and least by orbital interactions (22–16%). The neutral water adducts where X = permethylated N, O, P, or S are energetically nearly identical to the ethane adduct (2a–5a). Like with ethane, 4a converged in a O–H···C hydrogen bonding interaction, which can be rationalized by the lower (Pauling) electronegativity of P (2.19) compared to C (2.55). 2a, 3a and 5a converged at a geometry consistent with a Carbon bonding interaction. This is illustrated for structure 3a in Figure 5, where a single bond critical point (bcp) is located between C and O with a bond density of 0.60 · 102 a.u.. This can be rationalized by the higher (Pauling) electronegativity of N (3.05), O (3.44) and S (2.58) compared to C (2.55). Electrostatics or dispersion are the main energetic stabilizing factor in adducts 2a–5a which is typical for weak and non-directional interactions like in adduct 1a. A similar series with dma was computed as adducts 2b–5b. Only 2b converged in a geometry consistent with a Carbon bonding interaction (with dispersion as main driver) while the others are C–H···O hydrogen bonding interactions. Carbon bonding interactions with regular permethylated main group elements are thus comparable to very weak C–H hydrogen bonding interactions and less than about −1.5 kcal·mol−1 in strength. These energies are in line with earlier computations with the adducts [H2N-CH3···OCH2] [47] and [HO-CH3···OH2] [48] of −0.7 and −1.0 kcal·mol-1 respectively.
The cationic adducts 6–11a were computed as well and the most stable of these involved the Me3C+ carbocation in 6 (adducts with pentamethylated Carbon are unstable). The bonding interaction in 6b is largely covalent, as evidenced by the interaction energy of −82.4 kcal·mol−1, the large orbital contribution (55%), a dense bond critical point (18.4 · 102 a.u.) and a clear pyramidalization of the central C-atom (see Figure S4). Of the other adducts, all except 6a (Me3C+···O interaction) and 9 (with the least electronegative P) converged into an X–CH3···O Carbon bonding geometry. This is illustrated for 8a and 11a in Figure 4, where a clear bcp can be seen in between methylC and Owater with a bond density of 1.15 · 102 and 0.98 · 102 a.u. for 8a and 11a respectively. The bonding energies in 7–11a are mainly electrostatic in origin (~70%) and about −8 kcal·mol−1 for water and −15 kcal·mol−1 for dma. The most stable adducts in both series involved the most electronegative O (3.44) in Me3O+ (8). Two alternative configurations with N-methylpyridinium were also computed (12 and 13). In 12, the O points in between two CH hydrogens as is illustrated for 12a in Figure 4. In adducts 13 the O atom is located directly above the cationic N+. Both 12 and 13 are more stable than the Carbon bonding geometry found in 11, suggesting that hydrogen bonding interactions are most preferred. The interaction energies of Carbon bonding interaction with cationic species is similar to previous data of the adducts: [H3N+-CH3···OCH2] (−9.7 kcal·mol−1); [47] [Me3N+-CH3···OC(H)NH2] (−13 kcal·mol−1); [49] [Me2S+-CH3···OH2/NH3/OCH2] (about −8–9 kcal·mol−1); [47,49] and [R2S+-CH3···various lone-pairs] (about −9.0 kcal·mol−1) [50].
As the calculations with cationic species imply that electron withdrawing substituents amplify the Carbon bonding interaction, it was decided to compute adducts with small molecules that have an electron withdrawing group: Iodomethane (14), 1,1,1-trifluoroethane (16), acetonitrile (18) and nitromethane (20). All these adducts converged as a Carbon bonding geometry and are energetically favorable by 1.6–2.8 kcal·mol−1 for water and 3.0–4.9 kcal·mol−1 for dma. The adducts involving nitromethane (20) were most stable and are shown in Figure 5, together with an aim analysis revealing a single C···O bcp (ρ = 0.79 · 102 and 0.83 · 102 a.u. for 20a and 20a respectively). For the water adducts, H-bonding geometries were also optimized: 15a, 17a, 19a and 21a. All these adducts were about twice as stable and the Carbon bonding geometry. This can be ascribed to the fact that besides a C–H···O hydrogen bonding interaction, another weak hydrogen bonding interaction is present as well (i.e., C–H···I in 15a, C–H···F in 17a and C–H···π in 19a, see also see also Figure S4). For example, a dimer of HCF3 is estimated at −2.6 kcal·mol−1 and exhibits two weak C–H···F hydrogen bonding interactions (not shown). The interaction energies with neutral yet polarized methyl groups (14–20) liken those reported by others for: [F-CH3···OH2/NH3/PH3] (about −2–3 kcal·mol−1); [48,51,52,53] [F-CH3···C2H2] (−1.2 kcal·mol−1); [54] [Hlg-CH3···C2H4/NH3/PH3] (−1–4 kcal·mol−1); [53,55] [NC-CH3···C2H4/dma] (−3–5 kcal·mol−1); [55,56] and [O2N-CH3···dma] (−4.9 kcal·mol−1) [56].

3.4. General Discussion

From all the calculations collected in Table 1 it is evident that the interaction energies of Carbon bonding geometries with the sp2-O in dma (adducts ‘b’) is consistently about twice as strong as the interaction with sp3-O in water (adducts ‘a’). This is in line with the larger amount of van der Waals overlap observed in the N(d’) plots (Figure 2, ~30% for amides vs ~15% for water). The interaction energies of adducts with a Carbon bonding geometry range from very weak (below −1.5 kcal·mol−1 in 2, 3, 5), to moderately weak (between −1.5 and −5 kcal·mol−1 in 14, 16, 18 and 20) to fairly strong in the cationic adducts (between −7 and −18 kcal·mol−1 in 6–11). ΔE becomes smaller (more stable) in the order 2 < 14 < 16 < 18 < 20 < 7 < 8. Within this series, the orbital contribution remains constant at about 15–20%, while the electrostatic component increases from 30–35% in 2 to about 65% in 8. This implies that stronger Carbon bonding interactions are mainly driven by electrostatic interactions and that weaker such adducts are driven by dispersion. These computational results are consistent with recent literature reports [47,48,49,50,51,52,53,54,55,56,57,58,59] and the database analyses presented here; neutral adducts are very weak and thus hardly (or non) directional but can be made stronger (and thus presumably more directional) when X in X–CH3 is strongly polarized (see especially Figure 4). The relevance of Carbon bonding interactions with methyl groups is thus likely limited to highly polarized and/or cationic species. While this limits the scope considerably, it is worth pointing out that ligands with methyl groups related to those in adducts 2–10 and 14–21 are abundant within proteins structures and that cationic methyl groups also occur. For example, methylated methionine residues in methyl transferases [80] and nicotinamide derivatives such as nicotinamide adenine dinucleotide [81,82].

4. Summary and Conclusions

The CSD and the PDB were systematically evaluated for potential directional behavior of intermolecular non-covalent Carbon bonding interactions involving X–CH3 and electron rich entities such as O/S atoms or an aryl ring (ElR) within a hemisphere of 5 Å basal radius (centered on C). It was found that X–CH3···ElR interactions can be as directional as very weak hydrogen bonding interaction involving C–H (Pmax ≤ 1.50) but not directional at all when X = C. Grouping of data with significant amounts of van der Waals overlap (up to ~30%) was observed in various sub-datasets in the region where the X–CH3···ElR angle α is 160°–180°. These distributions were significantly shifted to shorter distances (i.e., more van der Waals overlap) in the case of cationic R3N+–CH3···Owater/amide compared to charge-neutral R2N–CH3···Owater/amide interactions.
Model DFT calculations revealed that charge neutral X–CH3···O adducts with water and dimethylacetamide are very weak (≤ –1.5 kcal·mol−1 in 2, 3a, 5a) and are often not the energy minima of the adducts (1, 3b, 4, 5b). The interaction energies can be increased by deploying a more electron withdrawing X (–1.5 to –5 kcal·mol−1 in 14, 16, 18 and 20). Rendering X cationic leads to even more stable adducts (–7.0 to –18 kcal·mol−1) in 7, 8, 10 and 11). Carbon-bonding adducts with dimethylacetamide are consistently twice as stable as those with water. Energy decomposition analyses showed that increased stability is driven by electrostatics and atom-in-molecule analyses regularly gave a clear bond critical point involving the methyl C-atom.
It is thus concluded that this combined database / DFT study reaffirms that intermolecular non-covalent Carbon interactions with X–CH3 is electrostatically driven and can be significant. The interaction can even by mildly directional in the solid state (comparable to weak CH hydrogen bonding interactions), provided X is sufficiently electron withdrawing.

Supplementary Materials

The following are available online, Table S1. Numerical overview of datasets retrieved from the CSD and the PDB. See Figure S1 and Figure 2 for P (α) directionality plots. For data of refined datasets see Table S2 (α ≥ 160° and van der Waals overlap). ‘n.a.’ stands for ‘not assessed’. The data for quaternary (cationic) N-atoms is not used in the paper for P(α) plot but they were used in N(d’) plots. Table S2. Numerical overview of the amount of data in each dataset from Table S1 characterized by α ≥ 160° (Nα ≥ 160°), together with the amount of data within that dataset where van der Waals radii overlap (NΣvdW). ‘n.a.’ stands for ‘not assessed’. Datasets with Nα ≥ 160° larger than 100 were inspected further by means of the N(d’) plots shown in Figures S2–S3. As a guide to the eye, datasets with less than 100 hits are in grey, those between 100 and 500 hits are in blue and those above 500 are in regular black. Figure S1. P (α) directionality plots for the data retrieved from the CSD (left) and the PDB (right) using the general query shown in the top-left inset figure for X–CH3∙∙∙ElR pairs. X can be C, N, O, P, or S and ‘ElR’ can be a water, amide or carboxyl O-atom, an RyCS S-atom (y = 2 or 3, R = any non-metal) or the centroid of an aryl ring, as is indicated in the right-hand side of the figure. The insert figure in the top right is intended as a guide to the eye to interpret the spatial location of data with a certain value of α. Due to the amount of data per dataset (see Table S1 for numerical overview), the top four P (α) plots are given at a 5° resolution for α and the bottom six at a 10° resolution. NB: Interestingly, in the P(α) plots for X = S and ElR = water, P is above unity at α = 160°–180° for the CSD data, while P < 1 for the PDB data in this same region. In both databases, the P-values are ≥ 1 around α = 105°. This indicates that the H-bonding geometry (α ≈ 105°) is somewhat directional on both databases but that the carbon bonding geometry is more directional only in the CSD. However, the dataset retrieved from the CSD is much smaller (1,546 hits) than the dataset from the PDB (30,725 hits) and the observed feature α = 160°–180° in the CSD might well be an artefact. Similarly, the P(α) plots with X = S and an amide O-atom reveal that P ≥ 1 at α = 90°–105°, again congruent with a hydrogen bonding geometry. P ≥ 1 also at α = 170°–180, but only for the PDB data. As this dataset is more voluminous (N = 11,215 vs 2,175 in the CSD), this implies that a carbon bonding geometry is more directional (at least in protein structures). A possible reason for this discrepancy might be that many cysteine and methionine residues are involved in metal coordination or directly methylated thus polarizing the S–C bond. For example, iron-sulphur clusters are held in place by cysteine–Fe coordination bonds [1] and methylated methionine residues are a (crystallographically known) intermediate in methyl transferases. [2]. Figure S2. Cumulative hit fraction (in %) as a function of the van der Waals corrected XCH3···ElR distance d’ (in Å) for several datasets from Table S2 (α ≥ 160°, as illustrated by the inset figure). The interacting pairs involve water-O (red, empty) or amide-O (blue, half-filled) with X = C (squares), N (diamonds), or O (circles). See Figure 2 for the same data plotted as a regular hit fraction. Figure S3. Cumulative (top) and regular (bottom) hit fraction (in %) as a function of the van der Waals corrected XCH3···ElR distance d’ (in Å) for all the datasets from Table S2 (α ≥ 160°, illustrated with the inset figure) that contain ≥ 500 data points. The inset legends show the nature of X (horizontally) and ElR (vertically). NB: In addition to the trends observed and discussed in the main text according to Figure 2, It is interesting to note that in all cases, for X = C the plot is shifted most to longer d’ and has the least van der Waals overlap. Moreover, carboxy O-atoms are distributed about the same as amide O-atoms and thio S-atoms about the same as water. Aryls are always randomly distributed with the least amount of van der Waals overlap. Figure S4. Ball and stick representations of perspective views of all molecular adducts listed in Table 1 that were optimized by DFT (B3LYP-D3/def2-TZVP). The thin lines are bond paths (bp’s) and the small red spheres are bond critical points (bcp’s) obtained from an ‘atoms-in-molecules’ analysis. The bond density (ρ) is in arbitrary units 102 and bcp’s indicative of non-covalent carbon bonding have been highlighted in yellow. Figure S5. Ball and stick representations of molecular adducts selected from Table 1 that were optimized by DFT (B3LYP-D3/def2-TZVP). The thin lines are bond paths (bp’s) and the small red spheres are bond critical points (bcp’s) obtained from an ‘atoms-in-molecules’ analysis. The bond density (ρ) is in arbitrary units ·102 and bcp’s indicative of non-covalent carbon bonding have been highlighted in yellow.

Funding

This research was funded by the Netherlands Organization for Scientific Research (NWO) grant number 723.015.006.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Lehn, J.M. Supramolecular Chemistry: Concepts and Perspectives, 1st ed.; Wiley VCH: Weinheim, Germany, 1995. [Google Scholar]
  2. Clark, T. σ-Holes. WIREs Comput. Mol. Sci. 2013, 3, 13–20. [Google Scholar] [CrossRef]
  3. Politzer, P.; Murray, J.S.; Clark, T. Halogen bonding and other σ-hole interactions: A perspective. Phys. Chem. Chem. Phys. 2013, 15, 11178–11189. [Google Scholar] [CrossRef] [PubMed]
  4. Bauza, A.; Mooibroek, T.J.; Frontera, A. The bright future of unconventional σ/π-hole interactions. Chem. Phys. Chem. 2015, 16, 2496–2517. [Google Scholar] [CrossRef] [PubMed]
  5. Crabtree, R.H. Hypervalency, secondary bonding and hydrogen bonding: Siblings under the skin. Chem. Soc. Rev. 2017, 46, 1720–1729. [Google Scholar] [CrossRef] [PubMed]
  6. Grabowski, S.J. Hydrogen Bonding: New Insights; Springer: Heidelberg, Germany, 2006. [Google Scholar]
  7. Bauza, A.; Mooibroek, T.J.; Frontera, A. Tetrel-bonding interaction: Rediscovered supramolecular force? Angew. Chem. Int. Ed. 2013, 52, 12317–12321. [Google Scholar] [CrossRef] [PubMed]
  8. Bauza, A.; Mooibroek, T.J.; Frontera, A. Tetrel bonding interactions. Chem. Rec. 2016, 16, 473–487. [Google Scholar] [CrossRef] [PubMed]
  9. Bauza, A.; Frontera, A.; Mooibroek, T.J. 1,1,2,2-Tetracyanocyclopropane (TCCP) as supramolecular synthon. Phys. Chem. Chem. Phys. 2016, 18, 1693–1698. [Google Scholar] [CrossRef]
  10. Bauza, A.; Mooibroek, T.J.; Frontera, A. Influence of ring size on the strength of carbon bonding complexes between anions and perfluorocycloalkanes. Phys. Chem. Chem. Phys. 2014, 16, 19192–19197. [Google Scholar] [CrossRef]
  11. Bauza, A.; Mooibroek, T.J.; Frontera, A. Non-covalent sp3 carbon bonding with ArCF3 is analogous to CH-π interactions. Chem. Commun. 2014, 50, 12626–12629. [Google Scholar] [CrossRef]
  12. Scholfield, M.R.; van der Zanden, C.M.; Carter, M.; Ho, P.S. Halogen bonding (X-bonding): A biological perspective. Protein Sci. 2013, 22, 139–152. [Google Scholar] [CrossRef]
  13. Cavallo, G.; Metrangolo, P.; Milani, R.; Pilati, T.; Priimagi, A.; Resnati, G.; Terraneo, G. The halogen bond. Chem. Rev. 2016, 116, 2478–2601. [Google Scholar] [CrossRef] [PubMed]
  14. Auffinger, P.; Hays, F.A.; Westhof, E.; Ho, P.S. Halogen bonds in biological molecules. Proc. Natl. Acad. Sci. USA 2004, 101, 16789–16794. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Voth, A.R.; Khuu, P.; Oishi, K.; Ho, P.S. Halogen bonds as orthogonal molecular interactions to hydrogen bonds. Nat. Chem. 2009, 1, 74–79. [Google Scholar] [CrossRef] [PubMed]
  17. Hardegger, L.A.; Kuhn, B.; Spinnler, B.; Anselm, L.; Ecabert, R.; Stihle, M.; Gsell, B.; Thoma, R.; Diez, J.; Benz, J.; et al. Systematic investigation of halogen bonding in protein-ligand interactions. Angew. Chem. Int. Ed. 2011, 50, 314–318. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, Q.; Xu, Z.J.; Shi, J.Y.; Zhu, W.L. Underestimated halogen bonds forming with protein backbone in protein data bank. J. Chem Inf. Model. 2017, 57, 1529–1534. [Google Scholar] [CrossRef] [PubMed]
  19. Ho, P.S. Halogen bonding in medicinal chemistry: From observation to prediction. Future Med. Chem. 2017, 9, 637–640. [Google Scholar] [CrossRef] [PubMed]
  20. Mendez, L.; Henriquez, G.; Sirimulla, S.; Narayan, M. Looking back, looking forward at halogen bonding in drug discovery. Molecules 2017, 22, 1397. [Google Scholar] [CrossRef] [PubMed]
  21. Voth, A.R.; Hays, F.A.; Ho, P.S. Directing macromolecular conformation through halogen bonds. Proc. Natl. Acad. Sci. USA 2007, 104, 6188–6193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Danelius, E.; Andersson, H.; Jarvoll, P.; Lood, K.; Grafenstein, J.; Erdelyi, M. Halogen bonding: A powerful tool for modulation of peptide conformation. Biochemistry 2017, 56, 3265–3272. [Google Scholar] [CrossRef] [PubMed]
  23. Murray, J.S.; Lane, P.; Clark, T.; Riley, K.E.; Politzer, P. σ-holes, π-holes and electrostatically-driven interactions. J. Mol. Model. 2012, 18, 541–548. [Google Scholar] [CrossRef] [PubMed]
  24. Murrayrust, P.; Burgi, H.B.; Dunitz, J.D. Chemical reaction paths. 5. SN1 reaction of tetrahedral molecules. J. Am. Chem. Soc. 1975, 97, 921–922. [Google Scholar] [CrossRef]
  25. Choudhary, A.; Gandla, D.; Krow, G.R.; Raines, R.T. Nature of Amide Carbonyl-Carbonyl Interactions in Proteins. J. Am. Chem. Soc. 2009, 131, 7244–7246. [Google Scholar] [CrossRef] [PubMed]
  26. Bartlett, G.J.; Choudhary, A.; Raines, R.T.; Woolfson, D.N. n → π* interactions in proteins. Nat. Chem. Biol. 2010, 6, 615–620. [Google Scholar] [CrossRef] [PubMed]
  27. Newberry, R.W.; Bartlett, G.J.; VanVeller, B.; Woolfson, D.N.; Raines, R.T. Signatures of n → π* interactions in proteins. Protein Sci. 2014, 23, 284–288. [Google Scholar] [CrossRef]
  28. Schottel, B.L.; Chifotides, H.T.; Dunbar, K.R. Anion-π interactions. Chem. Soc. Rev. 2008, 37, 68–83. [Google Scholar] [CrossRef]
  29. Mooibroek, T.J.; Gamez, P. Anion-arene and lone pair-arene interactions are directional. Cryst. Eng. Comm. 2012, 14, 1027–1030. [Google Scholar] [CrossRef]
  30. Frontera, A.; Gamez, P.; Mascal, M.; Mooibroek, T.J.; Reedijk, J. Putting anion-π interactions into perspective. Angew. Chem. Int. Ed. 2011, 50, 9564–9583. [Google Scholar] [CrossRef]
  31. Mooibroek, T.J.; Gamez, P.; Reedijk, J. Lone pair-π interactions: A new supramolecular bond? Cryst. Eng. Comm. 2008, 10, 1501–1515. [Google Scholar] [CrossRef]
  32. Burgi, H.B. Stereochemistry of reaction paths as determined from crystal-structure data - relationship between structure and energy. Angew. Chem. Int. Ed. 1975, 14, 460–473. [Google Scholar] [CrossRef]
  33. Doppert, M.T.; van Overeem, H.; Mooibroek, T.J. Intermolecular π-Hole / n→π* interactions with carbon monoxide ligands in crystal structures. Chem. Commun. 2018, 54, 12049–12052. [Google Scholar] [CrossRef] [PubMed]
  34. Ruigrok van der Werve, A.; van Dijk, Y.R.; Mooibroek, T.J. π-Hole / n→π* interactions with acetonitrile in crystal structures. Chem. Commun. 2018, 54, 10742–10745. [Google Scholar] [CrossRef] [PubMed]
  35. Lewinski, J.; Bury, W.; Justyniak, W. Significance of intermolecular S center dot center dot center dot C(π) interaction involving M-S and -C=O centers in crystal structures of metal thiolate complexes. Eur. J. Inorg. Chem. 2005. [Google Scholar] [CrossRef]
  36. Echeverria, J. The n → π* interaction in metal complexes. Chem. Commun. 2018, 54, 3061–3064. [Google Scholar] [CrossRef]
  37. Gamez, P.; Mooibroek, T.J.; Teat, S.J.; Reedijk, J. Anion binding involving π-acidic heteroaromatic rings. Acc. Chem. Res. 2007, 40, 435–444. [Google Scholar] [CrossRef] [PubMed]
  38. Gamez, P. The anion-π interaction: Naissance and establishment of a peculiar supramolecular bond. Inorg. Chem. Front. 2014, 1, 35–43. [Google Scholar] [CrossRef]
  39. Echeverria, J. Intermolecular Carbonyl ··· Carbonyl Interactions in Transition-Metal Complexes. Inorg. Chem. 2018, 57, 5429–5437. [Google Scholar] [CrossRef]
  40. Mooibroek, T.J. Coordinated nitrate anions can be directional π-hole donors in the solid state: A CSD study. Cryst. Eng. Comm. 2017, 19, 4485–4488. [Google Scholar] [CrossRef]
  41. Bauza, A.; Frontera, A.; Mooibroek, T.J. NO3– anions can act as Lewis acid in the solid state. Nat. Commun. 2017, 8, 14522. [Google Scholar] [CrossRef]
  42. Bauza, A.; Frontera, A.; Mooibroek, T.J. π-Hole interactions involving nitro compounds: Directionality of nitrate esters. Cryst. Growth Des. 2016, 16, 5520–5524. [Google Scholar] [CrossRef]
  43. Franconetti, A.; Frontera, A.; Mooibroek, T.J. Intramolecular π-hole interactions with nitro aromatics. Cryst. Eng. Comm. 2019. [Google Scholar] [CrossRef]
  44. Bauza, A.; Frontera, A.; Mooibroek, T.J. π-Hole Interactions Involving Nitro Aromatic Ligands Within Protein Structures. Chem. Eur. J. 2019. [Google Scholar] [CrossRef]
  45. Li, W.; Spada, L.; Tasinato, N.; Rampino, S.; Evangelisti, L.; Gualandi, A.; Cozzi, P.G.; Melandri, S.; Barone, V.; Puzzarini, C. Theory Meets Experiment for Noncovalent Complexes: The Puzzling Case of Pnicogen Interactions. Angew. Chem. Int. Ed. 2018, 57, 13853–13857. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Grabowski, S.J. Tetrel bond-sigma-hole bond as a preliminary stage of the S(N)2 reaction. Phys. Chem. Chem. Phys. 2014, 16, 1824–1834. [Google Scholar] [CrossRef] [PubMed]
  47. Bauza, A.; Frontera, A. RCH3 center dot center dot center dot O Interactions in Biological Systems: Are They Trifurcated H-Bonds or Noncovalent Carbon Bonds? Crystals 2016, 6, 26. [Google Scholar] [CrossRef]
  48. Mani, D.; Arunan, E. The X-C ··· Y (X = O/F, Y = O/S/F/Cl/Br/N/P) ’carbon bond’ and hydrophobic interactions. Phys. Chem. Chem. Phys. 2013, 15, 14377–14383. [Google Scholar] [CrossRef] [PubMed]
  49. Scheiner, S. Ability of IR and NMR Spectral Data to Distinguish between a Tetrel Bond and a Hydrogen Bond. J. Phys. Chem. A 2018, 122, 7852–7862. [Google Scholar] [CrossRef] [Green Version]
  50. Trievel, R.C.; Scheiner, S. Crystallographic and Computational Characterization of Methyl Tetrel Bonding in S-Adenosylmethionine-Dependent Methyltransferases. Molecules 2018, 23, 2965. [Google Scholar] [CrossRef]
  51. Sethio, D.; Oliveira, V.; Kraka, E. Quantitative Assessment of Tetrel Bonding Utilizing Vibrational Spectroscopy. Molecules 2018, 23, 2763. [Google Scholar] [CrossRef]
  52. Scheiner, S. Systematic Elucidation of Factors That Influence the Strength of Tetrel Bonds. J. Phys. Chem. A 2017, 121, 5561–5568. [Google Scholar] [CrossRef] [Green Version]
  53. Laconsay, C.J.; Galbraith, J.M. A valence bond theory treatment of tetrel bonding interactions. Comput. Theor. Chem. 2017, 1116, 202–206. [Google Scholar] [CrossRef]
  54. Wei, Y.X.; Li, H.B.; Cheng, J.B.; Li, W.Z.; Li, Q.Z. Prominent enhancing effects of substituents on the strength of π···σ-hole tetrel bond. Int. J. Quantum Chem. 2017, 117, e25448. [Google Scholar] [CrossRef]
  55. Mani, D.; Arunan, E. The X-C ··· π (X = F, Cl, Br, CN) Carbon Bond. J. Phys. Chem. A 2014, 118, 10081–10089. [Google Scholar] [CrossRef] [PubMed]
  56. Mundlapati, V.R.; Sahoo, D.K.; Bhaumik, S.; Jena, S.; Chandrakar, A.; Biswal, H.S. Noncovalent Carbon-Bonding Interactions in Proteins. Angew. Chem. Int. Ed. 2018, 57, 16496–16500. [Google Scholar] [CrossRef] [PubMed]
  57. Varadwaj, P.R.; Varadwaj, A.; Jin, B.Y. Significant evidence of C ··· O and C center dot center dot center dot C long-range contacts in several heterodimeric complexes of CO with CH3-X, should one refer to them as carbon and dicarbon bonds! Phys. Chem. Chem. Phys. 2014, 16, 17238–17252. [Google Scholar] [CrossRef] [PubMed]
  58. Garcia-Llinas, X.; Bauza, A.; Seth, S.K.; Frontera, A. Importance of R-CF3 center dot center dot center dot O Tetrel Bonding Interactions in Biological Systems. J. Phys. Chem. A 2017, 121, 5371–5376. [Google Scholar] [CrossRef] [PubMed]
  59. Thomas, S.P.; Pavan, M.S.; Row, T.N.G. Experimental evidence for ’carbon bonding’ in the solid state from charge density analysis. Chem. Commun. 2014, 50, 49–51. [Google Scholar] [CrossRef] [PubMed]
  60. Allen, F.H. The Cambridge Structural Database: A quarter of a million crystal structures and rising. Acta Crystallogr. Sect. B-Struct. Sci. 2002, 58, 380–388. [Google Scholar] [CrossRef] [PubMed]
  61. Groom, C.R.; Bruno, I.J.; Lightfoot, M.P.; Ward, S.C. The Cambridge Structural Database. Acta Crystallogr. Sect. B-Struct. Sci. 2016, 72, 171–179. [Google Scholar] [CrossRef] [PubMed]
  62. Bruno, I.J.; Cole, J.C.; Edgington, P.R.; Kessler, M.; Macrae, C.F.; McCabe, P.; Pearson, J.; Taylor, R. New software for searching the Cambridge Structural Database and visualizing crystal structures. Acta Crystallogr. Sect. B-Struct. Sci. 2002, 58, 389–397. [Google Scholar] [CrossRef]
  63. Hendlich, M.; Bergner, A.; Gunther, J.; Klebe, G. Relibase: Design and development of a database for comprehensive analysis of protein-ligand interactions. J. Mol. Biol. 2003, 326, 607–620. [Google Scholar] [CrossRef]
  64. Mooibroek, T.J.; Gamez, P. Directional character of solvent- and anion-pentafluorophenyl supramolecular interactions. Cryst. Eng. Comm. 2012, 14, 3902–3906. [Google Scholar] [CrossRef]
  65. Mooibroek, T.J.; Gamez, P. How directional are D-H center dot center dot center dot phenyl interactions in the solid state (D = C, N, O)? CrystEngComm 2012, 14, 8462–8467. [Google Scholar] [CrossRef]
  66. Mooibroek, T.J.; Gamez, P. Halogen ··· phenyl supramolecular interactions in the solid state: Hydrogen versus halogen bonding and directionality. CrystEngComm 2013, 15, 1802–1805. [Google Scholar] [CrossRef]
  67. Mooibroek, T.J.; Gamez, P. Halogen bonding versus hydrogen bonding: What does the Cambridge Database reveal? Cryst. Eng. Comm. 2013, 15, 4565–4570. [Google Scholar] [CrossRef]
  68. Bauza, A.; Mooibroek, T.J.; Frontera, A. Directionality of π-holes in nitro compounds. Chem. Commun. 2015, 51, 1491–1493. [Google Scholar] [CrossRef] [PubMed]
  69. Allen, F.H.; Kennard, O.; Watson, D.G.; Brammer, L.; Orpen, A.G.; Taylor, R. Tables of bond lengths determined by X-ray and neutron-diffraction. 1. Bond lengths in organic compounds. J. Chem. Soc. Perkin Trans. 2 1987, S1–S19. [Google Scholar] [CrossRef]
  70. Bondi, A. van der Waals Volumes and Radii. J. Phys. Chem. 1964, 68, 441–452. [Google Scholar] [CrossRef]
  71. Lee, C.T.; Yang, W.T.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef] [Green Version]
  72. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic-behaviour. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  73. 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]
  74. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef]
  75. Weigend, F. Accurate Coulomb-fitting basis sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057–1065. [Google Scholar] [CrossRef]
  76. te Velde, G.; Bickelhaupt, F.M.; Baerends, E.J.; Guerra, C.F.; van Gisbergen, S.J.A.; Snijders, J.G.; Ziegler, T. Chemistry with ADF. J. Comput. Chem. 2001, 22, 931–967. [Google Scholar] [CrossRef]
  77. Bader, R.F.W. Atoms in Molecules. Acc. Chem. Res. 1985, 18, 9–15. [Google Scholar] [CrossRef]
  78. Bickelhaupt, F.M.; Baerends, E.J. Kohn-Sham density functional theory: Predicting and understanding chemistry. In Reviews in Computational Chemistry, Vol 15; Lipkowitz, K.B., Boyd, D.B., Eds.; Wiley-Vch, Inc.: New York, NY, USA, 2000; Volume 15, pp. 1–86. [Google Scholar]
  79. van der Lubbe, S.C.C.; Guerra, C.F. The Nature of Hydrogen Bonds: A Delineation of the Role of Different Energy Components on Hydrogen Bond Strengths and Lengths. Chem. Asian J. 2019, 14, 2760–2769. [Google Scholar] [CrossRef] [Green Version]
  80. Schubert, H.L.; Blumenthal, R.M.; Cheng, X.D. Many paths to methyltransfer: A chronicle of convergence. Trends Biochem.Sci. 2003, 28, 329–335. [Google Scholar] [CrossRef]
  81. Foster, J.W.; Moat, A.G. Nicotinamide Adenosine-dinucleotide Biosynthesis and Pyridine-nucleotide Cycle Metabilism in Microbial Systems. Microbiol. Rev. 1980, 44, 83–105. [Google Scholar]
  82. Pollak, N.; Dolle, C.; Ziegler, M. The power to reduce: Pyridine nucleotides - small molecules with a multitude of functions. Biochem. J. 2007, 402, 205–218. [Google Scholar] [CrossRef]
Molecules 24 03370 g001 550
Figure 1. Representation of the method used to retrieve and analyse data from the CSD and the PDB. (a) general query to obtain data with d ≤ 5Å, α = 90°–180°, X = C, N, O, P or S and ElR (electron rich entity) is as indicated. (b) Illustration of the method used to assess directionality (see text for details).
Figure 1. Representation of the method used to retrieve and analyse data from the CSD and the PDB. (a) general query to obtain data with d ≤ 5Å, α = 90°–180°, X = C, N, O, P or S and ElR (electron rich entity) is as indicated. (b) Illustration of the method used to assess directionality (see text for details).
Molecules 24 03370 g001
Molecules 24 03370 g002 550
Figure 2. P (α) directionality plots for the data retrieved from the CSD (left) and the PDB (right) using the general query shown in the top-right inset figure for X–CH3∙∙∙ElR pairs. X can be C, N or O and ‘ElR’ can be a water or an amide O-atom. The insert figure in the top left is intended as a guide to the eye to interpret the spatial location of data with a certain value of α. Due to the amount of data per dataset (see Table S1 for numerical overview), the plots are given at a 5° resolution for α. A full set of P(α) plots (i.e., for all the X vs ElR pairs in Figure 1) is given in Figure S1. The P value of 1 is highlighted in green and indicates a random distribution of data. N (CSD/PDB) = 46,000/29,508 (C, water); 18,170/17,101 (N, water); 7190/11,392 (O, water); 53,473/22,538 (C, amide); 7663/10,855 (N, amide); 9,158/11,064 (O, amide).
Figure 2. P (α) directionality plots for the data retrieved from the CSD (left) and the PDB (right) using the general query shown in the top-right inset figure for X–CH3∙∙∙ElR pairs. X can be C, N or O and ‘ElR’ can be a water or an amide O-atom. The insert figure in the top left is intended as a guide to the eye to interpret the spatial location of data with a certain value of α. Due to the amount of data per dataset (see Table S1 for numerical overview), the plots are given at a 5° resolution for α. A full set of P(α) plots (i.e., for all the X vs ElR pairs in Figure 1) is given in Figure S1. The P value of 1 is highlighted in green and indicates a random distribution of data. N (CSD/PDB) = 46,000/29,508 (C, water); 18,170/17,101 (N, water); 7190/11,392 (O, water); 53,473/22,538 (C, amide); 7663/10,855 (N, amide); 9,158/11,064 (O, amide).
Molecules 24 03370 g002
Molecules 24 03370 g003 550
Figure 3. Hit fraction (in %) as a function of the van der Waals corrected XCH3···ElR distance d’ (in Å) for several datasets from Table S2 (α ≥ 160°, as illustrated by the inset figure). The interacting pairs involve water-O (red, empty) or amide-O (blue, half-filled) with X = C (squares), N (diamonds), or O (circles). See Figure S2 for the same data plotted as a cumulative hit fraction. N (CSD/PDB) = 2,376/1,622 (C, water); 1,089/1,186 (N, water); 452/757 (O, water); 2,599/1,179 (C, amide); 504/650 (N, amide); 640/483 (O, amide).
Figure 3. Hit fraction (in %) as a function of the van der Waals corrected XCH3···ElR distance d’ (in Å) for several datasets from Table S2 (α ≥ 160°, as illustrated by the inset figure). The interacting pairs involve water-O (red, empty) or amide-O (blue, half-filled) with X = C (squares), N (diamonds), or O (circles). See Figure S2 for the same data plotted as a cumulative hit fraction. N (CSD/PDB) = 2,376/1,622 (C, water); 1,089/1,186 (N, water); 452/757 (O, water); 2,599/1,179 (C, amide); 504/650 (N, amide); 640/483 (O, amide).
Molecules 24 03370 g003
Molecules 24 03370 g004 550
Figure 4. Cumulative (left) and regular (right) hit fraction (in %) as a function of the van der Waals corrected XCH3···ElR distance d’ (in Å) for the datasets from Table S2 (α = 160°–180° as illustrated in the inset figure). The interacting pairs involve water-O (red, empty), carboxy-O (green, right-filled) or the centroid of an aryl ring (grey, left-filled). X can be any N (diamond) or a cationic (tetravalent) N+ (hexagonal). N (N/N+) = 1,809/290 (water, red); 2,062/274 (carboxy, green); 10,585/527 (aryl, grey).
Figure 4. Cumulative (left) and regular (right) hit fraction (in %) as a function of the van der Waals corrected XCH3···ElR distance d’ (in Å) for the datasets from Table S2 (α = 160°–180° as illustrated in the inset figure). The interacting pairs involve water-O (red, empty), carboxy-O (green, right-filled) or the centroid of an aryl ring (grey, left-filled). X can be any N (diamond) or a cationic (tetravalent) N+ (hexagonal). N (N/N+) = 1,809/290 (water, red); 2,062/274 (carboxy, green); 10,585/527 (aryl, grey).
Molecules 24 03370 g004
Molecules 24 03370 g005 550
Figure 5. Ball and stick representations of molecular adducts selected from Table 1 that were optimized by DFT (B3LYP-D3/def2-TZVP). The thin lines are bond paths (bp’s) and the small red spheres are bond critical points (bcp’s) obtained from an ‘atoms-in-molecules’ analysis. The bond density (ρ) is in arbitrary units ·102 and bcp’s indicative of non-covalent Carbon bonding have been highlighted in yellow.
Figure 5. Ball and stick representations of molecular adducts selected from Table 1 that were optimized by DFT (B3LYP-D3/def2-TZVP). The thin lines are bond paths (bp’s) and the small red spheres are bond critical points (bcp’s) obtained from an ‘atoms-in-molecules’ analysis. The bond density (ρ) is in arbitrary units ·102 and bcp’s indicative of non-covalent Carbon bonding have been highlighted in yellow.
Molecules 24 03370 g005
Table
Table 1. Numerical overview of adducts computed with DFT between an indicated X-CH3 methyls and water (adducts ‘a’) or dimethylacetamide (dma, adducts ‘b’). Using ADF at the B3LYP-D3/def2-TZVP level of theory, interaction energies (in kcal·mol−1) were computed and an energy decomposition analyses is shown as a percentage of the total amount of interaction energies split up as electrostatic (E), orbital (O) and dispersion (D) interactions. Entries in grey are not consistent with a Carbon bonding geometry. See Figure S4 for perspective views and atoms in molecules analyses.
Table 1. Numerical overview of adducts computed with DFT between an indicated X-CH3 methyls and water (adducts ‘a’) or dimethylacetamide (dma, adducts ‘b’). Using ADF at the B3LYP-D3/def2-TZVP level of theory, interaction energies (in kcal·mol−1) were computed and an energy decomposition analyses is shown as a percentage of the total amount of interaction energies split up as electrostatic (E), orbital (O) and dispersion (D) interactions. Entries in grey are not consistent with a Carbon bonding geometry. See Figure S4 for perspective views and atoms in molecules analyses.
Interacting X–CH3:Water–ODimethylacetamide (dma) – C=O
adductα (°)ΔEE/O/D in %adductα (°)dEE/O/D in %
Me–CH31a170a−1.232/22/46%1b122f−2.826/16/58%
Me2N–CH32a174b−0.836/18/45%2b173b−1.633/22/44%
MeO–CH33a175 b−1.248/17/35%3b164f−3.227/16/56%
Me2P–CH34a167a−1.228/32/41%4b91d−4.245/21/33%
MeS–CH35a177 b−1.044/17/38%5b100d−3.746/21/32%
Me2C+–CH36a75c−12.666/22/12%6b109c−82.443/55/2%
Me3N+–CH37a173b−7.873/15/12%7b170b−14.770/19/10%
Me2O+–CH38a176b−9.172/18/10%8b175b−17.669/22/8%
Me3P+–CH39a55d−8.269/14/18%9b71d−21.764/23/13%
Me2S+–CH310a168b−7.673/15/12%10a170b−15.070/20/10%
PyN+–CH311a170b−7.373/15/12%11b175b−14.271/19/10%
12a93d−10.670/22/9%12b91d−20.365/27/8%
13a68e−9.070/14/16%13b74e−18.264/21/15%
I–CH3g14a175b−1.653/18/28%14b175b−3.052/20/28%
15a84d−3.855/23/21%
F3C–CH3 g16a179b−1.859/13/27%16b175b−3.459/15/26%
17a76d−3.860/15/24%
N≡C–CH3 g18a158b−2.365/13/23%18b172b−4.362/15/23%
19a68d−5.163/19/18%
O2N–CH3 g20a165b−2.866/13/21%20b179b−4.964/15/21%
21a86d−6.463/25/12%
a One of the water H-atoms and one of the RCH3 C-atoms are closest to each other; b Carbon bonding interaction geometry; cInteraction with the cationic C; d Hydrogen bonding interaction(s); e Interaction with cationic N; f CH-π interaction; g The interaction energy with benzene was also computed, starting from a geometry with X–CH3···benzene centroid = 180°. All adducts converged at a geometry where this angle was about 90° degrees. Interaction energies are about 4–5 kcal·mol−1 and dominated by electrostatics (35–40%) and dispersion (40–50%, see Figure S5 for details).

Share and Cite

MDPI and ACS Style

Mooibroek, T.J. Intermolecular Non-Covalent Carbon-Bonding Interactions with Methyl Groups: A CSD, PDB and DFT Study. Molecules 2019, 24, 3370. https://doi.org/10.3390/molecules24183370

AMA Style

Mooibroek TJ. Intermolecular Non-Covalent Carbon-Bonding Interactions with Methyl Groups: A CSD, PDB and DFT Study. Molecules. 2019; 24(18):3370. https://doi.org/10.3390/molecules24183370

Chicago/Turabian Style

Mooibroek, Tiddo J. 2019. "Intermolecular Non-Covalent Carbon-Bonding Interactions with Methyl Groups: A CSD, PDB and DFT Study" Molecules 24, no. 18: 3370. https://doi.org/10.3390/molecules24183370

APA Style

Mooibroek, T. J. (2019). Intermolecular Non-Covalent Carbon-Bonding Interactions with Methyl Groups: A CSD, PDB and DFT Study. Molecules, 24(18), 3370. https://doi.org/10.3390/molecules24183370

Article Metrics

Back to TopTop