Next Article in Journal
Soil-Borne Neosartorya spp.: A Heat-Resistant Fungal Threat to Horticulture and Food Production—An Important Component of the Root-Associated Microbial Community
Next Article in Special Issue
Excess Absorbance as a Novel Approach for Studying the Self-Aggregation of Vital Dyes in Liquid Solution
Previous Article in Journal
The Circadian Clock, Nutritional Signals and Reproduction: A Close Relationship
Previous Article in Special Issue
Predicting the Electronic Absorption Band Shape of Azobenzene Photoswitches
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 2
10.3390/ijms24021542
ijms-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

Unraveling the Nature of Hydrogen Bonds of “Proton Sponges” Based on Car-Parrinello and Metadynamics Approaches

by
Beata Kizior
1,
Mariusz Michalczyk
1,
Jarosław J. Panek
2,
Wiktor Zierkiewicz
1 and
Aneta Jezierska
2,*
1
Faculty of Chemistry, Wrocław University of Science and Technology, ul. Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
2
Faculty of Chemistry, University of Wrocław, ul. F. Joliot-Curie 14, 50-383 Wrocław, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(2), 1542; https://doi.org/10.3390/ijms24021542
Submission received: 13 November 2022 / Revised: 26 December 2022 / Accepted: 5 January 2023 / Published: 12 January 2023
(This article belongs to the Special Issue Feature Papers in Physical Chemistry and Chemical Physics 2022)
Browse Figures
Versions Notes

Abstract

:
The nature of intra- and intermolecular non-covalent interactions was studied in four naphthalene derivatives commonly referred to as “proton sponges”. Special attention was paid to an intramolecular hydrogen bond present in the protonated form of the compounds. The unsubstituted “proton sponge” served as a reference structure to study the substituent influence on the hydrogen bond (HB) properties. We selected three compounds substituted by methoxy, amino, and nitro groups. The presence of the substituents either retained the parent symmetry or rendered the compounds asymmetric. In order to reveal the non-covalent interaction properties, the Hirshfeld surface (HS) was computed for the crystal structures of the studied compounds. Next, quantum-chemical simulations were performed in vacuo and in the crystalline phase. Car–Parrinello molecular dynamics (CPMD), Path Integral Molecular Dynamics (PIMD), and metadynamics were employed to investigate the time-evolution changes of metric parameters and free energy profile in both phases. Additionally, for selected snapshots obtained from the CPMD trajectories, non-covalent interactions and electronic structure were studied. Quantum theory of atoms in molecules (QTAIM) and the Density Overlap Regions Indicator (DORI) were applied for this purpose. It was found based on Hirshfeld surfaces that, besides intramolecular hydrogen bonds, other non-covalent interactions are present and have a strong impact on the crystal structure organization. The CPMD results obtained in both phases showed frequent proton transfer phenomena. The proton was strongly delocalized in the applied time-scale and temperature, especially in the PIMD framework. The use of metadynamics allowed for tracing the free energy profiles and confirming that the hydrogen bonds present in “proton sponges” are Low-Barrier Hydrogen Bonds (LBHBs). The electronic and topological analysis quantitatively described the temperature dependence and time-evolution changes of the electronic structure. The covalency of the hydrogen bonds was estimated based on QTAIM analysis. It was found that strong hydrogen bonds show greater covalency, which is additionally determined by the proton position in the hydrogen bridge.

1. Introduction

Molecular dynamics (MD) methods play a key role in studying changes in system properties as a function of time at the molecular level [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. The classical MD simulations are realized by solving the corresponding Newton’s equations of motion for systems of interacting particles. The inter-particle forces and related potential energies are often calculated based on interatomic potentials or classical molecular mechanics force fields. These methods are used in chemical physics, biophysics, chemistry and materials science [15,16,17]. One of the most important methods belonging to ab initio time-evolution schemes is Car–Parrinello Molecular Dynamics (CPMD) [2]. The main idea of the CPMD method is to separate the fast motion of electrons from the slow motion of atomic nuclei on a time scale, and then replace the dynamics of real electrons with the dynamics of fictitious particles—Kohn-Sham orbital coefficients.The slow motion of atomic nuclei is then characterized based on the classical (Newtonian) mechanics. In turn, the fast motion of electrons (or rather changes in the orbital composition) is determined by the Density Functional Theory (DFT) with reference to the Kohn–Sham equations [18]. This approach allows a detailed discussion of effects related to temperature or pressure as well as entropic changes in the investigated chemical species. This method allows a satisfactory description to be obtained of the nature of interactions in the studied systems. Thus, it is one of the methods applied very often to obtain a detailed characteristic of intra- and intermolecular interactions; see, e.g., Refs. [19,20,21,22,23,24,25,26,27,28,29,30,31]. These interactions play a fundamental role in chemistry, biochemistry, and biology in determining the structure and physico-chemical properties of systems [32,33,34]. In order to describe the quantum nature of nuclei, Path Integral Molecular Dynamics (PIMD) was developed [6,35]. The method was successfully applied in our numerous studies of intra- and intermolecular interactions [23,36,37].
The presence of intramolecular hydrogen bond results very often in the quasi-ring formation, which stabilizes the spatial structure of compounds [31,38]. Other examples could be the catalysis of reactions involving enzymes and proteins [39]; the formation of secondary, tertiary, and quaternary structures in proteins [40,41]; and the modulation of the electronic structure and aromaticity of compounds [42,43]. It is well known that hydrogen bonds play an important role in proton transfer (PT) processes [44,45], as well as in interactions during charge transfer [46]. Moreover, they can change and affect the resonance capabilities of the substituents (electron-donating or -accepting) in the molecules [42]. Hence, it can be said that the role of intramolecular hydrogen bonding is very significant.
In the current study, quantum-chemical results for the protonated form of 1,8-bis(dimethyl- amino)naphthalene (DMANH + ) and its three derivatives (see Figure 1 and Figure S1 of the Supplementary Materials) are reported [47,48,49,50]. These compounds are known as “proton sponges” in the literature. The first report characterizing 1,8-bis(dimethylamino)naphthalene (DMAN) was published in 1968 by Alder et al. [51]. The DMAN molecule contains two NMe 2 groups in close proximity, so the attachment of the proton leads to the formation of a short intramolecular N-H⋯N hydrogen bond. Recent studies indicate that N⋯N distance can change in the range of 2.522–2.621 Å [52,53,54,55,56,57,58]. It should be noted that the shortest N⋯N distances are found in systems with a symmetric HB [53,54,55,57,59,60]. The intramolecular HB in the protonated form of 1,8-bis(dimethylamino)naphthalene has been described as a strong Low-Barrier Hydrogen Bond (LBHB) [27,61,62]. In the literature it was reported that LBHB can be of a single or double potential minimum type [63]. The IR spectra recorded for “proton sponges” exhibit an extremely low frequency absorption band in the 500–600 cm − 1 region, indicating the presence of N-H⋯N moiety [64]. In turn, the NMR spectroscopy 15 N– 15 N chemical coupling constant ( 2 J N N ) of the intramolecular symmetric HB of the “proton sponge” has a value of 11 Hz [64,65,66]. In addition to their spectroscopic features, these systems exhibit excellent base properties [67,68,69]. The first reported value of pKa for 1,8-bis(dimethylamino)naphthalene (DMAN) was 12.34 [51]. Various proton sponges have been synthesized with much better base properties than DMAN until now [62,68,70]. Therefore, it can be concluded that proton sponges are an interesting group of compounds to study. These features often receive credit in materials chemistry [71], as well as in catalysis, such as in the activation of CO 2 [72]. A relatively small number of available experimental and theoretical studies, as well as short interatomic separation in the N-H⋯N moiety [64], large values of the 15 N– 15 N coupling constants [64,65,66], and other interesting physico-chemical features of DMANH + and its derivatives, prompted us to investigate the dynamic nature of the intramolecular hydrogen present in the investigated “proton sponges” (see Figure 1).
These types of studies are a valuable source of information used in crystal engineering for the rational design of new systems with the desired physico-chemical properties and structural features. Following this line of thinking, our work focused on the theoretical investigation of intramolecular (N-H⋯N) hydrogen bonds in the gas and crystalline phases. The Hirshfeld surface (HS) [73,74] and fingerprint [75,76] approaches were used to reflect intra- and intermolecular interactions in the crystal structures. Quantum-chemical simulations were carried out using Car–Parrinello Molecular Dynamics (CPMD) [2]. Finally, the Quantum Theory of Atoms in Molecules (QTAIM) and Density Overlap Regions Indicator (DORI) approaches were used to obtain detailed characteristics of intramolecular hydrogen bonding and non-covalent interactions in the prototypic symmetric DMAN compound, its methoxy derivative, and asymmetric derivatives—substituted by amino and nitro group, respectively. This study aims to deepen the understanding of the nature of intramolecular interactions in the particular “proton sponges” and to answer the following questions:
(i)
What is the mechanism of proton motion in an intramolecular N-H⋯N bond in symmetric/asymmetric “proton sponges”?
(ii)
Does the proton transfer play a determinant role in changing the physico-chemical features in selected “proton sponges”?
(iii)
What is the environmental impact on the intramolecular hydrogen bond features?
(iv)
What is the temperature impact on the dynamics of the intramolecular hydrogen bond?
To the best of our knowledge, this is one of the few attempts to use Car–Parrinello Molecular Dynamics (CPMD) [2] and Path Integrals Molecular Dynamics (PIMD) [6,35] to study HB dynamics in symmetric and asymmetric “proton sponges”. One of our earlier studies [29] was devoted to the CPMD investigation of proton dynamics in symmetric “proton sponges” including DMANH + , but in the gas phase only and with emphasis on the vibrational signatures. The particular feature of the current work is a comparison of dynamical nature of the intramolecular hydrogen bond and its properties as a function of time in two different environments—gas phase vs. crystalline phase—and of temperatures.

2. Results and Discussion

2.1. Hirshfeld Surface Analysis

The Hirshfeld surface (HS) analysis of the investigated symmetric and asymmetric “proton sponges” was carried out and are illustrated in Figure 2 and Figure 3 (see Figure S2 for relevant unit cells). The size and shape of the Hirshfeld surface allow us to study and show qualitatively intra- and intermolecular interactions in molecular crystals [73,77]. This analysis describes the d e (distance from the Hirshfeld surface to the nearest atom outside the surface), d i (distance from the Hirshfeld surface to the nearest atom inside the surface) as well as d n o r m (normalized contact distance) in crystal structures [73,74]. In Hirshfeld maps, the interactions with a distance equal to the sum of the van der Waals (VDW) radii are illustrated with the white surface. In turn, red areas characterize shorter interactions, and blue indicates interactions with distances longer than the VDW radii. The red surface around the two nitrogen atoms represents the negative potential (the acceptor), and the blue area around the hydrogen atoms illustrates the positive potential (the donor). The formation of a quasi-ring due to the presence of intramolecular HB has a significant impact on the structure and molecular properties of the studied compounds. It is worth noting that blue areas were revealed for compound (2), indicating the presence of interactions with distances greater than VDW radii. From the HS analysis results, we can conclude that the intermolecular interactions in the crystal structures play a key role in the crystal stability, as presented in Figure 3 and Figure S3 in the Supplementary Information.
In the analyzed symmetric and asymmetric “proton sponges”, the molecules are organized in the crystal through Cl⋯H and O⋯H, H⋯C and O⋯H, O⋯H, H⋯C and H⋯Br, O⋯H short contacts for compounds (1)–(4), respectively, (see Figure 3 and Table 1).
It was found that the shortest O⋯H contacts were present in compounds (3) and (4), as indicated by the red areas (see Figure 3). The different trend was noted for symmetrical compounds, where the areas characterizing the O⋯H interaction are much lighter. Thus, the molecular packing of the studied crystals is strongly influenced by the presence of intra- and intermolecular interactions. In the next step, the fingerprint plot approach was used to investigate short and long contacts in the studied crystals, as presented in Figure S3 in the Supplementary Information. The fingerprint maps show a picture similar to HS. The plot in Figure S3 shows the usual features confirming the presence of an intramolecular hydrogen bond. The obtained plots for the “proton sponges” contain expected intramolecular HB spikes at the bottom left of the maps (the lower spike, d i > d e ) (see Figure S3). All crystal structures exhibit diverse O⋯H contacts. Interestingly, the crystal structures of compounds (2) and (3) show “wings” in the upper left corner which proves the presence of C-H ⋯ π interactions in the studied structures and is consistent with the results obtained for the naphthalene ring [75]. Compounds (1) and (4), on the other hand, exhibit this property to a far lesser extent. In summary, the presented results allowed us to conclude that intramolecular hydrogen bonding is present in all symmetric/asymmetric “proton sponges” taken into consideration in the study. In addition, intermolecular interactions are fundamental for maintaining stability in the crystal structures of the compounds.

2.2. Car–Parrinello Molecular Dynamics (CPMD) in the Gas and Crystalline Phases

The electronic structure framework behind the Car–Parrinello Molecular Dynamics (CPMD) scheme, namely Density Functional Theory (DFT), has profound advantage over the classical force fields: it allows for a description of the bond breaking and forming processes without any special provisions. Thus, the DFT-based CPMD is ideally suited for the studies of proton dynamics in hydrogen bonds, where short-lived or permanent proton transfer (PT) events can occur, and where a priori knowledge of proton positions is not always possible.
The data on the evolution of the relevant N-H⋯N bond lengths as functions of time are shown in Figure 4, Figure 5, Figure 6, Figure 7 and Figures S4–S7 (for the OLYP functional) in the Supplementary Information. In order to obtain detailed differences in the dynamics of the bridge protons for asymmetric and symmetric “proton sponges”, the CPMD simulations were performed in the gas and crystal phases at 100 K and 300 K using two functionals (PBE and OLYP). Figure 4 and Figure 5 show that for compounds (1)–(3), the bridge proton is located on the donor (N) side in the gas phase simulations. Interestingly, in the investigated compounds (1)–(3), the bridged proton does exhibit very often short contacts with the acceptor nitrogen (N) atom; therefore, it can be indeed described as very labile. The situation is different for an asymmetric “proton sponge” with the nitro group as the substituent (compound (4)), in which the proton present in the intramolecular hydrogen bridge is more often observed on the acceptor side (see Figure 4 and Figure 5). At this point, it is worthwhile to notice that the notions of the “donor” and “acceptor” are not well defined, because the neutral “proton sponge” catches the proton between the two nitrogen atoms, so we have chosen to denote as “donor” the nitrogen atom at position 1 of the naphthalene moiety. In the case of the asymmetric proton sponges, one could use the diffraction data to assign the “donor” and “acceptor” labels, but we prefer to use a consistent notation throughout this work. One last observation from the gas phase results: the events of proton entering the acceptor site are more frequent for the symmetric “proton sponges” (1) and (2).
The crystal environment seems to promote instantaneous proton sharing or proton transfer events, which are visibly more frequent in Figure 6 and Figure 7 than in Figure 4 and Figure 5. Compound (4) consistently exhibits dominant proton position differently from the compounds (1)–(3). However, the semi-quantitative nature of the distance timeline graphs prevents further detailed analysis solely on their basis. Before proceeding with statistical analysis of the trajectories, it is interesting to observe that the time evolution data for the hydrogen bridge structural parameters are very similar for the two proposed functionals in both environments and temperatures. Thus, the applied computational level satisfactorily reproduced the evolution of the proton in the hydrogen bridge as a function of time. In summary, the influence of environment and temperature plays a key role in our CPMD study in the gas and crystal phases, as they were able to alter the dynamics of bridged protons, introducing events not observed for the isolated molecule. A quantitative description of these phenomena was carried out on the basis of proton possession statistics.
Proton possession, defined as the percentage of the simulation time spent by the proton within either the donor or the acceptor site region, is an important characteristics of the symmetry of the proton potential energy well. In the current study, the simplest and most intuitive definition is used: the proton is in (temporary) possession by that of the two nitrogen atoms to which it is closest. As already noted above, the donor and acceptor concepts are not well established for the strong hydrogen bonding in the “proton sponges”; therefore, the data in Table 2 and Table S1 was calculated with reference to the nitrogen atom at the position 1 of the aromatic moiety. At the beginning of this part of data analysis, it should be emphasized that the results for the PBE functional reported in the article (Table 2) are in very close agreement with the OLYP functional calculations (reported in the Supporting Information, Table S1).
Comparison of the statistics in Table 2 for different temperatures revealed that the barrier between two potential energy minima is sufficient to severely hinder the proton mobility at 100 K. There are two possible sources of the perceived asymmetry of the proton possession data: either the minima are equivalent but the barrier cannot be easily crossed at 100 K (so that the proton is trapped in one of the minima), or the minima are in fact not equivalent. The former option happens in the gas phase, while the latter can result from the asymmetry of the crystal environment. The gas phase results for the symmetric compounds (1) and (2) reflect the effect of the barrier height by deviating from the ideal 50–50% proton possession. This deviation is larger at 100 K, up to 6% and up to 2.1% at 300 K. The influence of the crystalline environment is very significant at 100 K—the compounds which are symmetric by themselves suddenly lock the proton in one of the potential wells for ca. two thirds of the simulation time. This effect is much smaller at 300 K, where again the proton potential profile is approaching the symmetric double well shape. The situation is very different when the asymmetric compounds (3) and (4) are analyzed. The gas phase results are indicative of the substituent effect only, and the impact of two selected substituents at position 4, namely the -NH 2 group in (3) and the -NO 2 group in (4), is clearly significant indeed. While substitution by another amino group keeps the proton at the same ring where the additional -NH 2 functional group is present, the -NO 2 substitution strongly promotes the proton transition to the other half of the naphthalene moiety. An increase in the temperature from 100 K to 300 K makes the proton distribution more balanced by ca. 20%, e.g., from 91.5% to 74.2% for (3). The asymmetry of the preferred proton location agrees well with the substituent properties. The electron-withdrawing nitro group in position 4 makes position 1 more positively charged and thus less preferable to the proton. On the other hand, the electron-donating -NH 2 group at position 4 increases electron density at position 1, making the proton more willing to reside at this site. Interestingly, the crystalline phase results for (3) and (4) at 300 K are so similar to their gas phase counterparts that the effect of neighboring moieties seems to be non-existent. However, upon lowering the simulation temperature to 100 K, we found out that this effect is very significant for (3), counteracting the intramolecular modulation of the proton position. These findings underline the delicate nature of the interplay between intra- and inter-molecular factors determining the proton potential function in the studied compounds.

2.3. Path Integral Molecular Dynamics—Inclusion of Nuclear Quantum Effects

Car–Parrinello Molecular Dynamics is inherently classical (Newtonian) when the nuclear degrees of freedom are considered. However, light nuclei such as protons can exhibit strong nuclear quantum effects, including tunneling and probability distribution broadening. These effects can be efficiently sampled within the Path Integral Molecular Dynamics (PIMD) [6] framework. A serious drawback of PIMD is that it does not provide access to real-time quantum dynamics of the system; the ensemble averages are the properties that reflect the inclusion of quantum effects. Therefore, the results of PIMD cannot be presented as time evolution data such as in Figure 4, Figure 5, Figure 6 and Figure 7. We have converted both CPMD and PIMD trajectory data from the gas phase simulations at 300 K into probability distributions ρ (r), and then—via the formula F(r) = – k B T ln ρ (r)—we have arrived at the free energy landscapes presented in Figure 8. The landscapes are prepared as functions of the N-H distances in the intramolecular hydrogen bonds of the studied systems. The horizontal axis in each subfigure of Figure 8 corresponds to the “donor”-proton N–H distance, while the vertical one corresponds to the proton–“acceptor” distance.
Figure 8 shows that the asymmetry of the proton positions is visible for the molecules of (3) and (4) with classical (CPMD) description of nuclei. There is also a small difference between the CPMD results for (1) and (2)—the prototypic DMANH + (1) seems to have a small barrier for the proton transfer. When the nuclear quantum effects are included via PIMD, a significant broadening of the proton distribution is observed, and careful insight allows us to state that the N⋯N distance was shortened by quantum effects because the outer regions of the probability distributions for PIMD end at lower N—H distances than for CPMD. The quantized protons are located closer to the center of the bridge than in the CPMD, especially for the compounds (1)–(2). In these cases, the free energy lansdcape seems to be barrierless and effectively corresponds to a single-well profile. The asymmetry of the proton distribution is still visible for (3)–(4) and is in agreement with the proton possession data presented in the previous section.
The reported free-energy landscapes exhibit small coverage of the phase space, but, given the sterical restraints of the tight hydrogen bond of proton sponges, this is not surprising. The well depths were estimated at ca. 3 kcal/mol, which definitely does not correspond to the energy of strong, sterically and electronically assisted hydrogen bonds. This is, however, the result of relying on thermal energy only. Reproducing the hydrogen bond energy in detail requires the use of “accelerated sampling” techniques. One of such schemes is metadynamics, and its application to the proton sponges is presented below.

2.4. Metadynamics

The CPMD and PIMD results presented in the previous sections showed a significant range of modulation of the proton position in the hydrogen bond depending on substitution, the electrostatic field of the crystal, and the temperature. The metadynamics [78] study is aimed at reproducing the free energy surface for proton motion with a history-dependent potential so that the barriers can be easily crossed and their heights estimated more accurately than for standard CPMD. The latter case can utilize only the thermal energy to cross the barriers, while metadynamics is an extended theory aimed at overcoming the barriers via added Gaussian hills gradually filling up the current potential well.
The results for the gas phase metadynamics for compounds (1)–(4) are shown in Figure 9. The metadynamics was carried out in the phase space of three variables, e.g., N-H, H⋯N, and N⋯N distances, while the figure shows slices across the N⋯N axis. The chosen N⋯N value of 2.6 Å corresponds well to the average N-H distance in the studied compounds. The depths of the well vary significantly, from only −17.7 kcal/mol for (1) to −34 kcal/mol for (4). However, we have chosen in the Figure to describe the range of free energy from −20 to 0 kcal/mol, that is the fragment of the potential hypersurfaces corresponding to the most abundant N⋯N region. This has some bearing on the shape of the surfaces, especially for (4) where the true depth is larger, but this convention allows for easy comparison of the energy profiles. Starting from the prototypic DMANH + , compound (1), one can note that the free energy wells are different between 100 K and 300 K. The minima at 100 K are deeper, which is not surprising—three times smaller kinetic energy values than for 300 K allow for much less frequent visits of the proton to the acceptor side. Crossing the minima (lying at −17.7 kcal/mol) requires increase in the free energy to −10 kcal/mol, which turns out not to be too problematic even at 100 K, as seen in Figure 4 and in the proton possession data described above. In case of the compound (2), the minima are deeper, −22 kcal/mol, but still symmetrical. Raising the well bottoms at 300 K reveals the barrier at the middle of the bridge, i.e., at equal N-H and H⋯N distances (1.4 Å). The most curious are the cases of (3) and (4). Compound (3) has a strongly asymmetric well, biased towards the “donor” atom. On the contrary, the molecule of (4) tends to place the proton at the “acceptor” side, e.g., the ring devoid of substituents, with the depth reaching −34 kcal/mol. The reported depths of the free energy surface for all four compounds, from −17.7 to −34 kcal/mol, show indeed that the proton is strongly held by the NMe 2 groups forming characteristic “clamps”, and the shallow barriers indicate that the involved hydrogen bonds are of the LBHB class. It should be mentioned at the end of this section that even after introducing a significant amount of energy into the system, the metadynamics did not lead to the event of the proton being pushed out of the hydrogen bond area. The potential energy pumped into the system via the metadynamics quite often flows into unexpected channels, and the fact that the proton stays in the N-H⋯N bridge is one more manifestation of the specific mechanism of basicity of the “proton sponges”.

2.5. Static Density Functional Theory (DFT)

For the discussed “proton sponges”, the reaction path describing the proton transfer in the hydrogen bond was investigated. The simulations were performed at the PBE0-D3/def2-TZVP level of theory. The obtained results are shown in Figure S8. The relative energy barrier is less than 1 kcal/mol in all studied cases. Two energy minima were detected for each of the studied compounds. In the asymmetric molecules, the secondary minimum (located at the “acceptor” site for (3) and the “donor” site for (4), respectively) is ca. 0.5 kcal/mol above the primary minimum. Such a low energy barrier indicates that the bridged proton can easily move between the donor and acceptor atoms. The static DFT results correspond well with the CPMD, PIMD, and metadynamics findings.
The QTAIM parameters for the studied molecules are gathered in Table S2. Two bond critical points (BCPs) for N⋯H and H⋯N interactions, including protons fluctuating between dimethylamino (-NMe 2 ) groups, were considered for these systems. Differences between these parameters demonstrate which NMe 2 moiety the proton belongs to at a given time. At first glance, one may notice that all values of ρ at BCPs are in approximate range of 0.1–0.3 a.u., which signifies the certain amount of covalency in all the investigated interactions. The behavior of protons, which balances between two nitrogen atoms, is expressed by the difference in values presented in single-table cells for a given compound at 100 K or 300 K temperature. It must be recalled that in compounds (1) and (2), the dimethylamino groups attached to the ring are indistinguishable due to the high symmetry of the molecules. In turn, in compounds (3) and (4), the data in the upper row of a single table cell refer to the N⋯H interaction with the N atom of the -NMe 2 group attached to the substituted ring (in the sections above denoted as the “donor”), while the data in lower row concerns H⋯N interaction where N atom is derived from the moiety incorporated into the unsubstituted ring (denoted as the “acceptor”). For clarity, the QTAIM molecular diagram for compound (4) in four snapshot time projections is captured in Figure 10. For this compound at 100 K, one can see that in the first three snapshots (0, 8, and 16 ps), the proton is shifted into the acceptor nitrogen closer to the ring without the nitro functional group. It is represented by the values of ρ of about 0.2 a.u. in a lower row towards the ones of about 0.1 a.u. appearing in upper row. This is followed by other QTAIM parameters: for the stronger H⋯N interactions, the value of ∇ 2 ρ is negative (as opposed to a positive or minimally negative value of this parameter for a weaker N⋯H interaction), the ϵ (ellipticity) is lower (from 0.006 to 0.009), which means a lower amount of π component in this bond, and its higher dynamic stability, H (electron energy density), is more negative and finally the |V|/G ratio (potential energy density and kinetic energy density, respectively) is greater (significantly exceeds 1 while for the weaker contact, it is only slightly higher than 1). The combination of H, V, and G parameters is in line with the remaining QTAIM descriptors and accords with Bader’s QTAIM protocol indicating that the stronger H⋯N interaction is covalent in nature, while the weaker one admittedly bears some hallmarks of covalency but is much weaker. The situation is completely reversed for the fourth snapshot, where the QTAIM parameters point to the supremacy of the H⋯N interactions with hydrogen directed to the nitrogen inserted in amino group of the substituted ring. When switching temperature to 300 K, the proton transfer between -NMe 2 groups can be easily identified. At 0 and 16 ps, the one type of hydrogen bond is favored, while at 8 and 24 ps, the second one dominates. Similar patterns can be deduced for remaining compounds in different temperatures. As a general conclusion regarding QTAIM studies at 300 K, the proton motion is faster in this temperature as changes in the affiliation of hydrogen to interact with a particular nitrogen atom become more rapid than at 100 K. Only in a few cases can one consider the proton in a state of “suspension” between both competing groups. The compound (3) at 100 K captured at 8 ps can serve as such an example. All relevant quantities are identical or very similar for both BCPs in this situation.
These results can be put into the context of earlier studies of intermolecular [NHN] + bridges. Literature reports classified such structural motif as covalent bond between the one of the N atoms and hydrogen atom (N-H) along with the second, significantly weaker contact, which was treated as a relevant hydrogen bond (H⋯N). Recently, in the work of El-Emam et al. [79] the N-H⋯N hydrogen bond was found in the crystal structure of adamantane-1,3,4-thiadiazole hybrid derivatives. Despite being formally classified as a non-covalent interaction, the ρ at BCP connecting the analyzed atoms was in the range of 0.173–0.220 a.u. (near to our findings for the majority of the current cases). Therefore, it exhibits distinguishing characteristics of a shared bond. Other factors, including the positive Laplacian of electron density, H(r) < 0, and -V/G ratio > 1, show that the intermolecular N-H-N contact in this system has a bonding type that is halfway between shared and closed-shell interaction. Two additional investigations that are comparable to ours in terms of the type of N-H contact examined solely theoretical models. Deepa et al. [80] found that the N-H⋯N hydrogen bonds between the AT and GC nucleic acid base pairs and selected amino acid side chains were relatively weak. The ρ was between 0.032 and 0.049 a.u., and other parameters also took on values typical of weaker non-covalent interactions, namely ∇ 2 ρ >0 and an ε oscillating around 0.050–0.060, which are very far from the results obtained in the current work. Finally, Bavafa et al. [81] studied the interactions between nitrosamine and formamide or formic acid. The QTAIM analysis unveiled the presence of an NHN bridge in which the proton was probably shifted into one of the nitrogen atoms. The partition of this fragment into a strong covalent bond and a weak non-covalent contact was clearly indicated by the ratio of the ρ values at the proper BCPs, which was 0.317 to 0.0213 a.u. The Laplacian of ρ had a negative sign for a stronger one, while a positive sign for a weaker one. The H value had the reverse relationship; it was positive for the unshared interaction.
Unfortunately, the non-covalent index (NCI) analysis was unable to identify the studied interactions because of their high strength. Therefore, alternatively, the DORI approach [82] was applied for compound (4). The DORI molecular diagrams in Figure 11 illustrate the covalent (blue) and non-covalent (green) interacting regions for compound (4) at 100 K temperature at 0 and 24 ps. As can be seen from this figure, the whole N-H-N bridge is treated as stabilized by covalent interactions of comparable magnitude as those between carbon atoms in the ring or N-O bonds in the nitro group. The green regions signifying non-covalent contacts are shown in the case of weak O⋯H hydrogen interactions.

3. Materials and Methods

Hirshfeld surface (HS) [73,74] and fingerprint plot [75,76] methods were used to visualize and verify intra- and intermolecular interactions in the crystal structures of the investigated compounds—symmetric/asymmetric “proton sponges” [47,48,49,50]. These analyses were performed in the MultiWFN program [83,84]. Car–Parrinello molecular dynamics [2] was carried out using the CPMD 4.3 program [85]. Simulations were performed in the gas and crystalline phases. The models for the CPMD crystalline phase studies were prepared on the basis of crystal structures deposited in the Cambridge Crystallographic Data Centre (CCDC) [86] denoted as TAPCES with the deposition number 1266338 (1-dimethylammonio-8-dimethylaminonaphthlene pentachlorophenolate bis(pentachlorophenol) [47], RISBEA with the deposition number 122006 (1,8-bis(dimethylamino)-4,5-dimethoxynaphthalene hydrobromide) [48], XUCKAH with the deposition number 170012 (4-amino-1,8-bis(dimethylamino)naphthalene hydrobromide monohydrate) [49], and ZOSKEX with the deposition number 1315231 (4-nitro-1-dimethylamino-8-dimethylammonionaphthalene perchlorate) [50] (see Figure S2 in Supplementary Information and Table 3 for details). The gas-phase models (single molecules) were prepared using X-ray coordinates as well and were placed in cubic boxes with a = 12 Å (TAPCES), a = 15 Å (RISBEA), a = 14 Å (XUCKAH), and a = 14 Å (ZOSKEX), respectively. The DMANH + derivatives include either electron-donating (-NH 2 and -OCH 3 ) or -accepting (-NO 2 ) substituents attached in the para position to the naphthalene ring relative to the NMe 2 group (see Figure 1).
In the first step, all prepared models were optimized with the Hessian matrix initialized according to the Schlegel scheme [87]. Then, the CPMD runs were performed. They were divided into equilibation and production runs. Electronic exchange and correlation effects were treated using the PBE [88] and OLYP [89,90] functionals in both phases. The simulations were performed using pseudopotentials proposed by N. Troullier and J. L. Martins [91]. The dispersion corrections by Grimme (DFT-D2 method) were included to reproduce weak interactions in the studied structures [92]. The kinetic energy cutoff of 100 Ry was used to determine the upper limit of the plane-wave basis set. The fictitious orbital mass parameter was equal to 400 a.u. The time step was set to 3 a.u. Molecular dynamics simulations were carried out at 300 K temperature controlled by the Nosé–Hoover chains thermostat [93,94]. The initial equilibration period of the CPMD runs was accrued for 30,000 steps for the gas and crystalline phases and was excluded from the further trajectories analyses. The trajectories were collected for 24 ps in both phases.
In the next step of the investigatiob of molecular dynamics, nuclear quantum effects were modeled within the Path Integrals Molecular Dynamics (PIMD) approach [5,6] for the studied cationic forms of the proton sponges in the gas phase. The electronic structure was propagated according to the Car–Parrinello scheme with the same computational setup as given above for the CPMD runs with the PBE functional. Eight Trotter replicas and staging transformation were used as parts of the PIMD setup. A separate Nosé–Hoover chains thermostat was coupled with each degree of freedom (“massive” thermostatting), and the simulations were carried out at 100 K and 300 K. As in the CPMD run, the first 30,000 steps were taken as equilibration phase—a time much larger than needed for the convergence of the primitive and virial kinetic energy estimators—and the PIMD data collection lasted for a further 12 ps.
In addition, metadynamics [78] based on the CPMD computational setup was employed to study the proton potential profile in the N-H⋯N hydrogen-bonded moiety. The gas phase models of the studied “proton sponges” were investigated via metadynamics at 100 K and 300 K lasting for 24 ps. The N–H, H⋯N, and N⋯N interatomic distances were taken as collective variables (CVs). The time-dependent metadynamics potential consisted of Gaussian hills with an initial height of 0.001 a.u. and transversal width of 0.05 a.u., placed every 150–500 CPMD steps (depending on the actual displacement of the CVs from the previous position). The 3D volumetric free energy maps were then prepared using the Vreco utility written by Nisanth Nair and distributed along with the CPMD package [85].
The CPMD study was followed by analysis of snapshots from the dynamics to characterize the intramolecular hydrogen bonding. The Quantum Theory of Atoms in Molecules (QTAIM) [95] and Density Overlap Regions Indicator (DORI) [82] approaches were used to obtain an in-depth characterization of the intramolecular hydrogen bridge in the studied compounds. For this purpose, four snapshots of the CPMD trajectory in the gas phase for each compound and each temperature were taken (at 0, 8, 16, and 24 ps of simulation time) and analyzed, yielding 32 snapshots in total. This allowed us to track changes in the electron structure during the simulation time. Kohn–Sham equations were solved for the snapshot structures at the PBE0-D3/def2-TZVP level [88,96,97], and the corresponding wavefunctions were saved for subsequent QTAIM processing. These electronic structure calculations were carried out with the Gaussian16 Rev. C.01 program [98]. The QTAIM analysis (determination of electron density critical points and their local properties) was carried out with assistance of the AimAll [99] software, while the DORI surfaces were generated with the MultiWFN program [83,84]. The QTAIM studies were carried out at non-equilibrium structures, but the usual nomenclature, associated with standard equilibrium-structure QTAIM calculations, was retained for convenience, following the practice common in the literature [100,101,102,103]. In addition, the proton reaction path within the intramolecular hydrogen bridges was studied by means of the scan with geometry optimization approach (with 0.05 Å increment of the N-H distance, freezing of the NHN valence angle describing the HB and optimization of the remaining structural parameters of the investigated compounds) using the PBE0-D3/def2-TZVP level of theory. The data were analyzed and visualized using the VMD [104], Vesta [105], Chimera [106], Gnuplot [107], and GaussView6 [108] programs.

4. Conclusions

In the present work, we focused on characterizing intra- and intermolecular interactions (N-H⋯N) in symmetric and asymmetric “proton sponges” based on theoretical methods. The Hirshfeld surface and fingerprint analyses were applied to give insight into intra- and intermolecular interactions in the crystal structures of the studied compounds. The CPMD results obtained for the gas phase as well as crystalline phase revealed the time-evolution in hydrogen bridge metric parameters. Metadynamics simulations were employed to reconstruct free energy surfaces in the investigated “proton sponges”. In addition, the QTAIM and DORI approaches made it possible to determine the strength and nature of interactions based on electron density. Therefore, based on the research methods used and the results obtained, we came to the following conclusions:
  • The Hirshfeld Surface (HS) and fingerprint data show intra- and intermolecular interactions in the crystals of the studied “proton sponges”.
  • Car–Parrinello Molecular Dynamics results enabled the time-evolution analysis of the bridged proton in the intramolecular hydrogen bond. The proton is very labile, and its movement is governed by a double-well potential, which can be asymmetric. The broken symmetry in the symmetric “proton sponges” can have internal origin (not enough energy to cross the barriers-vanishing at 300 K) or external cause (asymmetry of the crystal environment). For the asymmetric compounds (3) and (4), the intramolecular factor is more important than the environment. Quantum effects reproduced via the PIMD approach lead to significant broadening of the bridge proton probability distribution.
  • Metadynamics has allowed us to estimate the depths of free energy wells, ranging from −17.7 to −34 kcal/mol. This biased molecular dynamics scheme has proven successful in investigating the double-well free energy surface, but the proton remained within the space between the two nitrogen atoms.
  • QTAIM studies showed that the stronger H⋯N interactions are covalent in the investigated compounds (1)–(3), while the weaker contacts exhibit only subtle covalency. A different situation is observed for compound (4) with the supremacy of the H⋯N interactions with the proton directed to the “donor” nitrogen atom of the amine group of the substituted ring.
  • DORI molecular plots revealed regions indicating the presence of covalent and non-covalent interactions for compound (4) at 100 K with simulation times of 0 and 24 ps.

Supplementary Materials

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

Author Contributions

Conceptualization, A.J.; methodology, A.J., W.Z.; validation, B.K., M.M.; formal analysis, B.K., M.M., J.J.P., W.Z., and A.J.; investigation, B.K., M.M., and J.J.P.; writing—original draft preparation, B.K., M.M., J.J.P., W.Z., and A.J.; writing—review and editing, B.K., M.M., J.J.P., W.Z., and A.J.; visualization, B.K. and M.M.; supervision, A.J. and W.Z.; project administration, A.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financed in part by a statutory activity subsidy from the Polish Ministry of Science and Higher Education for the Faculty of Chemistry of Wroclaw University of Science and Technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data gathered in this study are presented in the manuscript and the Supplementary Material.

Acknowledgments

The authors gratefully acknowledge the Academic Computing Centre Cyfronet-Kraków (Prometheus supercomputer, part of the PL-Grid infrastructure), as well as the Poznań Supercomputing and Networking Center (PSNC) for generous grants of computer time and facilities. In addition, the ARCHER2 UK National Supercomputing Service (https://www.archer2.ac.uk) is acknowledged for generous CPU time and facilities in the framework of the DECI–17 access program.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BCPBond Critical Point
CPMDCar–Parrinello Molecular Dynamics
CVsCollective Variables
DFTDensity Functional Theory
DORIDensity Overlap Regions Indicator
HBHydrogen Bond
IRInfrared Spectroscopy
LBHBLow-Barrier Hydrogen Bond
MDMolecular Dynamics
NCINon-covalent Interactions Index
NMRNuclear magnetic resonance
PIMDPath Integrals Molecular Dynamics
PTProton transfer
QTAIMQuantum Theory of Atoms in Molecules
RCPRing Critical Point
RDGReduced Density Gradient
VDWvan der Waals

References

  1. Rapaport, D. The Art of Molecular Dynamics Simulation; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar] [CrossRef] [Green Version]
  2. Car, R.; Parrinello, M. Unified Approach for Molecular Dynamics and Density-Functional Theory. Phys. Rev. Lett. 1985, 55, 2471–2474. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Svishchev, I.M.; Kusalik, P.G. Crystallization of Liquid Water in a Molecular Dynamics Simulation. Phys. Rev. Lett. 1994, 73, 975–978. [Google Scholar] [CrossRef] [Green Version]
  4. Head-Gordon, M.; Tully, J.C. Molecular dynamics with electronic frictions. J. Chem. Phys. 1995, 103, 10137–10145. [Google Scholar] [CrossRef]
  5. Marx, D.; Parrinello, M. The Effect of Quantum and Thermal Fluctuations on the Structure of the Floppy Molecule C2H3+. Science 1996, 271, 179–181. [Google Scholar] [CrossRef]
  6. Tuckerman, M.E.; Marx, D.; Klein, M.L.; Parrinello, M. Efficient and general algorithms for path integral Car–Parrinello molecular dynamics. J. Chem. Phys. 1996, 104, 5579–5588. [Google Scholar] [CrossRef]
  7. Tuckerman, M.E.; Marx, D.; Klein, M.L.; Parrinello, M. On the Quantum Nature of the Shared Proton in Hydrogen Bonds. Science 1997, 275, 817–820. [Google Scholar] [CrossRef]
  8. Frank, I.; Hutter, J.; Marx, D.; Parrinello, M. Molecular dynamics in low-spin excited states. J. Chem. Phys. 1998, 108, 4060–4069. [Google Scholar] [CrossRef]
  9. Tuckerman, M.E.; Marx, D. Heavy-Atom Skeleton Quantization and Proton Tunneling in “Intermediate-Barrier” Hydrogen Bonds. Phys. Rev. Lett. 2001, 86, 4946–4949. [Google Scholar] [CrossRef]
  10. Carloni, P.; Rothlisberger, U.; Parrinello, M. The Role and Perspective of Ab Initio Molecular Dynamics in the Study of Biological Systems. Acc. Chem. Res. 2002, 35, 455–464. [Google Scholar] [CrossRef]
  11. Doltsinis, N.L.; Marx, D. Nonadiabatic Car–Parrinello Molecular Dynamics. Phys. Rev. Lett. 2002, 88, 166402. [Google Scholar] [CrossRef] [Green Version]
  12. Scherlis, D.A.; Fattebert, J.L.; Gygi, F.; Cococcioni, M.; Marzari, N. A unified electrostatic and cavitation model for first-principles molecular dynamics in solution. J. Chem. Phys. 2006, 124, 074103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Berkelbach, T.C.; Lee, H.S.; Tuckerman, M.E. Concerted Hydrogen-Bond Dynamics in the Transport Mechanism of the Hydrated Proton: A First-Principles Molecular Dynamics Study. Phys. Rev. Lett. 2009, 103, 238302. [Google Scholar] [CrossRef] [PubMed]
  14. Curchod, B.F.E.; Martínez, T.J. Ab Initio Nonadiabatic Quantum Molecular Dynamics. Chem. Rev. 2018, 118, 3305–3336. [Google Scholar] [CrossRef] [PubMed]
  15. Maginn, E.J.; Elliott, J.R. Historical Perspective and Current Outlook for Molecular Dynamics As a Chemical Engineering Tool. Ind. Eng. Chem. Res. 2010, 49, 3059–3078. [Google Scholar] [CrossRef]
  16. Matthews, B.L.C. Molecular Dynamics: With Deterministic and Stochastic Numerical Methods; Springer International Publishing: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  17. Alabugin, I.V.; Gilmore, K.M.; Peterson, P.W. Hyperconjugation. WIREs Comput. Mol. Sci. 2011, 1, 109–141. [Google Scholar] [CrossRef]
  18. Marx, D.; Hütter, J. Ab Initio Molecular Dynamics: Theory and Implementation in Modern Methods and Algorithms of Quantum Chemistry; Grotendorst, J., Ed.; NIC Series; John von Neumann Institute for Computing: Jülich, Germany, 2000. [Google Scholar]
  19. Fois, E.; Gamba, A. Dynamical Host/Guest Interactions in Zeolites: Framework Isotope Effects on Proton Transfer Studied by Car–Parrinello Molecular Dynamics. J. Phys. Chem. B 1999, 103, 1794–1799. [Google Scholar] [CrossRef]
  20. Boero, M.; Terakura, K.; Ikeshoji, T.; Liew, C.C.; Parrinello, M. Hydrogen Bonding and Dipole Moment of Water at Supercritical Conditions: A First-Principles Molecular Dynamics Study. Phys. Rev. Lett. 2000, 85, 3245–3248. [Google Scholar] [CrossRef] [PubMed]
  21. Stirling, A.; Bernasconi, M.; Parrinello, M. Ab initio simulation of water interaction with the (100) surface of pyrite. J. Chem. Phys. 2003, 118, 8917–8926. [Google Scholar] [CrossRef]
  22. Costanzo, F.; Valle, R.G.D. Car–Parrinello MD Simulations for the Na+-Phenylalanine Complex in Aqueous Solution. J. Phys. Chem. B 2008, 112, 12783–12789. [Google Scholar] [CrossRef]
  23. Jezierska, A.; Panek, J.J. Investigations of an O-H⋯S hydrogen bond via Car–Parrinello and path integral molecular dynamics. J. Comput. Chem. 2008, 30, 1241–1250. [Google Scholar] [CrossRef]
  24. Durlak, P.; Latajka, Z. Car–Parrinello and path integral molecular dynamics study of the intramolecular hydrogen bond in the novel class of anionic H-chelates: 6-Nitro-2, 3-dipyrrol-2-ylquinoxaline anion. Chem. Phys. Lett. 2009, 480, 173–177. [Google Scholar] [CrossRef]
  25. Shishkin, O.V.; Dopieralski, P.; Omelchenko, I.V.; Gorb, L.; Latajka, Z.; Leszczynski, J. Dynamical Nonplanarity of Benzene. Evidences from the Car–Parrinello Molecular Dynamics Study. J. Phys. Chem. Lett. 2011, 2, 2881–2884. [Google Scholar] [CrossRef]
  26. Dopieralski, P.; Perrin, C.L.; Latajka, Z. On the Intramolecular Hydrogen Bond in Solution: Car–Parrinello and Path Integral Molecular Dynamics Perspective. J. Chem. Theory Comput. 2011, 7, 3505–3513. [Google Scholar] [CrossRef] [PubMed]
  27. Jezierska, A.; Panek, J.J. “Zwitterionic Proton Sponge” Hydrogen Bonding Investigations on the Basis of Car–Parrinello Molecular Dynamics. J. Chem. Inf. Model. 2015, 55, 1148–1157. [Google Scholar] [CrossRef] [PubMed]
  28. Durlak, P.; Latajka, Z. Car–Parrinello and Path Integral Molecular Dynamics Study of the Proton Transfer in the Intramolecular Hydrogen Bonds in the Ketohydrazone–Azoenol System. J. Phys. Chem. B 2018, 122, 7862–7873. [Google Scholar] [CrossRef]
  29. Jezierska, A.; Panek, J.J. Theoretical study of intramolecular hydrogen bond in selected symmetric “proton sponges” on the basis of DFT and CPMD methods. J. Mol. Model. 2020, 26, 2–10. [Google Scholar] [CrossRef]
  30. Panek, J.J.; Zasada, J.; Szyja, B.M.; Kizior, B.; Jezierska, A. Sensitivity of Intra- and Intermolecular Interactions of Benzo[h]quinoline from Car–Parrinello Molecular Dynamics and Electronic Structure Inspection. Int. J. Mol. Sci. 2021, 22, 5220. [Google Scholar] [CrossRef]
  31. Kizior, B.; Panek, J.J.; Szyja, B.M.; Jezierska, A. Structure-Property Relationship in Selected Naphtho- and Anthra-Quinone Derivatives on the Basis of Density Functional Theory and Car–Parrinello Molecular Dynamics. Symmetry 2021, 13, 564. [Google Scholar] [CrossRef]
  32. Müller-Dethlefs, K.; Hobza, P. Noncovalent Interactions: A Challenge for Experiment and Theory. Chem. Rev. 2000, 100, 143–168. [Google Scholar] [CrossRef]
  33. Scheiner, S. Hydrogen Bonding; Oxford University Press: Oxford, UK, 1997. [Google Scholar]
  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] [Green Version]
  35. Feynman, R.P. Space-Time Approach to Non-Relativistic Quantum Mechanics. Rev. Mod. Phys. 1948, 20, 367–387. [Google Scholar] [CrossRef] [Green Version]
  36. Jezierska-Mazzarello, A.; Vuilleumier, R.; Panek, J.J.; Ciccotti, G. Molecular Property Investigations of an ortho-Hydroxy Schiff Base Type Compound with the First-Principle Molecular Dynamics Approach. J. Phys. Chem. B 2009, 114, 242–253. [Google Scholar] [CrossRef] [PubMed]
  37. Wojtkowiak, K.; Michalczyk, M.; Zierkiewicz, W.; Jezierska, A.; Panek, J.J. Chalcogen Bond as a Factor Stabilizing Ligand Conformation in the Binding Pocket of Carbonic Anhydrase IX Receptor Mimic. Int. J. Mol. Sci. 2022, 23, 13701. [Google Scholar] [CrossRef] [PubMed]
  38. Jezierska, A.; Kizior, B.; Szyja, B.M.; Panek, J.J. On the nature of inter- and intramolecular interactions involving benzo[h]quinoline and 10-hydroxybenzo[h]quinoline: Electronic ground state vs excited state study. J. Mol. Struct. 2021, 1234, 130126. [Google Scholar] [CrossRef]
  39. Simón, L.; Goodman, J.M. Enzyme Catalysis by Hydrogen Bonds: The Balance between Transition State Binding and Substrate Binding in Oxyanion Holes. J. Org. Chem. 2010, 75, 1831–1840. [Google Scholar] [CrossRef]
  40. Abeln, S.; Vendruscolo, M.; Dobson, C.M.; Frenkel, D. A Simple Lattice Model That Captures Protein Folding, Aggregation and Amyloid Formation. PLoS ONE 2014, 9, e85185. [Google Scholar] [CrossRef]
  41. Wieczorek, R.; Dannenberg, J.J. H-Bonding Cooperativity and Energetics of α-Helix Formation of Five 17-Amino Acid Peptides. J. Am. Chem. Soc. 2003, 125, 8124–8129. [Google Scholar] [CrossRef]
  42. Jezierska-Mazzarello, A.; Panek, J.J.; Szatyłowicz, H.; Krygowski, T.M. Hydrogen Bonding as a Modulator of Aromaticity and Electronic Structure of Selected ortho-Hydroxybenzaldehyde Derivatives. J. Phys. Chem. A 2012, 116, 460–475. [Google Scholar] [CrossRef]
  43. Pocheć, M.; Kułacz, K.; Panek, J.J.; Jezierska, A. How Substitution Combines with Non-Covalent Interactions to Modulate 1, 4-Naphthoquinone and Its Derivatives Molecular Features—Multifactor Studies. Int. J. Mol. Sci. 2021, 22, 10357. [Google Scholar] [CrossRef]
  44. Zeng, Y.; Li, A.; Yan, T. Hydrogen Bond Dynamics in the Solvation Shell on Proton Transfer in Aqueous Solution. J. Phys. Chem. B 2020, 124, 1817–1823. [Google Scholar] [CrossRef]
  45. Durlak, P.; Latajka, Z. Proton transfer dynamics in the propionic acid dimer from path integral molecular dynamics calculations. J. Mol. Model. 2011, 17, 2159–2168. [Google Scholar] [CrossRef] [PubMed]
  46. Pullanchery, S.; Kulik, S.; Rehl, B.; Hassanali, A.; Roke, S. Charge transfer across C–H⋯O hydrogen bonds stabilizes oil droplets in water. Science 2021, 374, 1366–1370. [Google Scholar] [CrossRef] [PubMed]
  47. Kanters, J.A.; ter Horst, E.H.; Kroon, J.; Grech, E. Complexes of the ’proton sponge’ 1, 8-bis(dimethylamino)naphthalene (DMAN). III. Structure of [DMANH]+.[pentachlorophenolate]−.[pentachlorophenol]2 at 100 K. Acta Crystallogr. C Struct. Chem. 1992, C48, 328–332. [Google Scholar] [CrossRef]
  48. Staab, H.A.; Kriege, C.; Hieber, G.; Oberdorf, K. 1,8-Bis(dimethylamino)4, 5-dihydroxynaphthalene, a Natural, Intramolecularly Protonated“Proton Sponge” with Zwitterionic Structure. Angew. Chem. Int. Ed. EngI. 1997, 36, 1884–1886. [Google Scholar] [CrossRef]
  49. Ozeryanskii, V.; Pozharskii, A.; Głowiak, T.; Majerz, I.; Sobczyk, L.; Grech, E.; Nowicka-Scheibe, J. X-ray diffraction and IR-spectroscopic studies on protonated 4-amino-1, 8-bis(dimethylamino)naphthalene. J. Mol. Struct. 2002, 607, 1–8. [Google Scholar] [CrossRef]
  50. Pozharskii, A.F.; Kuźmenko, V.V.; Aleksandrov, G.G.; Dmitrienko, D.V. 1,8-Bis(dimethylamino)naphthalene. XIII. Solvatochromism and molecular structure of 4-nitro-1,8-bis(dimethylamino)naphthalene and its salt with perchloric acid. Zhurnal Organicheskoi Khimii 1995, 31, 570. [Google Scholar]
  51. Alder, R.W.; Bowman, P.S.; Steele, W.R.S.; Winterman, D.R. The remarkable basicity of 1, 8-bis(dimethylamino)naphthalene. Chem. Commun. 1968, 13, 723–724. [Google Scholar] [CrossRef]
  52. Bakshi, P.K.; Cameron, T.S.; Knop, O. Crystal chemistry of tetraradial species. Part 8. Mix and match: Cation geometry, ion packing, hydrogen bonding, and π–π interactions in cis-2, 2′-bipyridinium(1+) and 1, 10-phenanthrolinium(1+) tetraphenylborates—and what about proton sponges? Can. J. Chem. 1996, 74, 201–220. [Google Scholar] [CrossRef]
  53. Pyżalska, D.; Pyżalski, R.; Borowiak, T. Structure of 1, 8-bis(dimethylamino)naphthalene hydrobromide dihydrate. J. Crystallogr. Spectrosc. Res. 1983, 13, 211–220. [Google Scholar] [CrossRef]
  54. Woźniak, K.; Krygowski, T.; Kariuki, B.; Jones, W.; Grech, E. Crystallographic studies on sterically affected chemical species Part II. Molecular and crystal structure of 1, 8-bis(dimethylamino)- naphthalene tetrafluoroborate. Analysis of distortion of geometry in the aromatic part due to intramolecular hydrogen bonding. J. Mol. Struct. 1990, 240, 111–118. [Google Scholar] [CrossRef]
  55. Kellett, P.J.; Anderson, O.P.; Strauss, S.H.; Abney, K.D. Solid-state motion of OTeF5− compounds: Detection by 19F NMR and IR spectroscopy and correlation with the X-ray structure of an orthorhombic crystalline modification of [C14H19N2+][OTeF5−] (C14H19N2+ = protonated 1, 8-bis(dimethylamino)naphthalene). Can. J. Chem. 1989, 67, 2023–2029. [Google Scholar] [CrossRef] [Green Version]
  56. Bartoszak, E.; Dega-Szafran, Z.; Grundwald-Wyspiańska, M.; Jaskólski, M.; Szafran, M. X-ray, Fourier-transform infrared, 1H and 13C nuclear magnetic resonance, and PM3 studies of (N—H⋯N)+ and (O—H⋯O)− intramolecular hydrogen bonds in a complex of 1, 8-bis(dimethylamino)naphthalene with maleic acid. J. Chem. Soc., Faraday Trans. 1993, 89, 2085–2094. [Google Scholar] [CrossRef]
  57. Bartoszak, E.; Jaskólski, M.; Grech, E.; Gustafsson, T.; Olovsson, I. Structure of thiocyanate salt of 1, 8-bis(dimethylamino)naphthalene (dmanH+.SCN-) at 188 and 290 K. Acta. Crystallogr. B. 1994, B50, 358–363. [Google Scholar] [CrossRef]
  58. Grech, E.; Malarski, Z.; Sawka-Dobrowolska, W.; Sobczyk, L. The structure and IR spectra of the 1:1 and 1:2 adducts of 1, 8-bis(dimethylamino)naphthalene (DMAN) with 4, 5-dicyanoimidazole (DCI). J. Mol. Struct. 1997, 406, 107–117. [Google Scholar] [CrossRef]
  59. Degtyarev, A.V.; Ryabtsova, O.V.; Pozharskii, A.F.; Ozeryanskii, V.A.; Starikova, Z.A.; Sobczyk, L.; Filarowski, A. 2, 7-Disubstituted proton sponges as borderline systems for investigating barrier-free intramolecular hydrogen bonds. Protonated 2, 7-bis(trimethylsilyl)- and 2, 7-di(hydroxymethyl)-1, 8-bis(dimethylamino)naphthalenes. Tetrahedron 2008, 64, 6209–6214. [Google Scholar] [CrossRef]
  60. Ozeryanskii, V.A.; Marchenko, A.V.; Pozharskii, A.F.; Filarowski, A.; Spiridonova, D.V. Combination of “Buttressing” and “Clothespin” Effects for Reaching the Shortest NHN Hydrogen Bond in Proton Sponge Cations. J. Org. Chem. 2021, 86, 3637–3647. [Google Scholar] [CrossRef]
  61. Frey, P.A.; Whitt, S.A.; Tobin, J.B. A Low-Barrier Hydrogen Bond in the Catalytic Triad of Serine Proteases. Science 1994, 264, 1927–1930. [Google Scholar] [CrossRef]
  62. Raab, V.; Kipke, J.; Gschwind, R.M.; Sundermeyer, J. 1,8-Bis(tetramethylguanidino)naphthalene (TMGN): A New, Superbasic and Kinetically Active “Proton Sponge”. Chem. Eur. J. 2002, 8, 1682–1693. [Google Scholar] [CrossRef]
  63. Hibbert, F.; Emsley, J. Hydrogen Bonding and Chemical Reactivity. Adv. Phys. Org. Chem. 1990, 26, 255–379. [Google Scholar] [CrossRef]
  64. Latajka, Z.; Sobczyk, L. The potential energy shape for the proton motion in protonated naphthalene proton sponges (DMANs) and its manifestations. In Practical Aspects of Computational Chemistry; Leszczyński, J., Shukla, M., Eds.; Springer: Berlin, Germany, 2009; Chapter 19. [Google Scholar]
  65. Pietrzak, M.; Wehling, J.; Limbach, H.H.; Golubev, N.S.; López, C.; Claramunt, R.M.; Elguero, J. 13C Detected Scalar Nitrogen-Nitrogen Couplings Across the Intramolecular Symmetric NHN Hydrogen Bond of Proton Sponge. J. Am. Chem. Soc. 2001, 123, 4338–4339. [Google Scholar] [CrossRef]
  66. Lloyd-Jones, G.C.; Harvey, J.N.; Hodgson, P.; Murray, M.; Woodward, R.L. Scalar Coupling between the 15N Centres in Methylated 1, 8-Diaminonaphthalenes and 1, 6-Diazacyclodecane: To What Extent is 2HJNN a Reliable Indicator of N—N Distance? Chem. Eur. J. 2003, 9, 4523–4535. [Google Scholar] [CrossRef] [PubMed]
  67. Raab, V.; Gauchenova, E.; Merkoulov, A.; Harms, K.; Sundermeyer, J.; Kovačević, B.; Maksić, Z.B. 1,8-Bis(hexamethyl-triaminophosphazenyl)naphthalene, HMPN: A Superbasic Bisphosphazene “Proton Sponge”. J. Am. Chem. Soc. 2005, 127, 15738–15743. [Google Scholar] [CrossRef] [PubMed]
  68. Alder, R.W. Strain effects on amine basicities. Chem. Rev. 1989, 89, 1215–1223. [Google Scholar] [CrossRef]
  69. Boiko, L.; Sorokin, V.; Filatova, E.; Starikova, Z.; Ozeryanskii, V.; Pozharskii, A. Three examples of naphthalene proton sponges with extreme or unusual structural parameters. General view on factors influencing proton sponge geometry. J. Mol. Struct. 2011, 1005, 12–16. [Google Scholar] [CrossRef]
  70. Raab, V.; Harms, K.; Sundermeyer, J.; Kovačević, B.; Maksić, Z.B. 1, 8-Bis(dimethylethyleneguanidino)naphthalene: Tailoring the Basicity of Bisguanidine “Proton Sponges” by Experiment and Theory. J. Org. Chem. 2003, 68, 8790–8797. [Google Scholar] [CrossRef]
  71. Martínez-Franco, R.; Sun, J.; Sastre, G.; Yun, Y.; Zou, X.; Moliner, M.; Corma, A. Supra-molecular assembly of aromatic proton sponges to direct the crystallization of extra-large-pore zeotypes. Proc. R. Soc. A 2014, 470, 20140107. [Google Scholar] [CrossRef] [Green Version]
  72. Sabet-Sarvestani, H.; Izadyar, M.; Eshghi, H.; Noroozi-Shad, N.; Bakavoli, M. Proton sponge as a new efficient catalyst for carbon dioxide transformation to methanol: Theoretical approach. Fuel 2018, 221, 491–500. [Google Scholar] [CrossRef]
  73. McKinnon, J.J.; Jayatilaka, D.; Spackman, M.A. Towards quantitative analysis of intermolecular interactions with Hirshfeld surfaces. Chem. Commun. 2007, 33, 3814–3816. [Google Scholar] [CrossRef]
  74. Spackman, M.A.; Jayatilaka, D. Hirshfeld surface analysis. CrystEngComm 2009, 11, 19–32. [Google Scholar] [CrossRef]
  75. Spackman, M.A.; McKinnon, J.J. Fingerprinting intermolecular interactions in molecular crystals. CrystEngComm 2002, 4, 378–392. [Google Scholar] [CrossRef]
  76. McKinnon, J.J.; Spackman, M.A.; Mitchell, A.S. Novel tools for visualizing and exploring intermolecular interactions in molecular crystals. Acta Crystallogr. B 2004, B60, 627–668. [Google Scholar] [CrossRef]
  77. Hachani, A.; Dridi, I.; Elleuch, S.; Roisnel, T.; Kefi, R. Crystal structure, spectroscopic and biological study of a new inorganic-organic hybrid compound [Cd4Cl12(H2O)2]n (C10N4H28)n. Inorg. Chem. Commun. 2019, 100, 134–143. [Google Scholar] [CrossRef]
  78. Laio, A.; Parrinello, M. Escaping free-energy minima. Proc. Natl. Acad. Sci. USA 2002, 99, 12562–12566. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. El-Emam, A.A.; Kumar, E.S.; Janani, K.; Al-Wahaibi, L.H.; Blacque, O.; El-Awady, M.I.; Al-Shaalan, N.H.; Percino, M.J.; Thamotharan, S. Quantitative assessment of the nature of noncovalent interactions in N-substituted-5-(adamantan-1-yl)-1, 3, 4-thiadiazole-2-amines: Insights from crystallographic and QTAIM analysis. RSC Adv. 2020, 10, 9840–9853. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Deepa, P.; Kolandaivel, P.; Senthilkumar, K. Hydrogen-bonding studies of amino acid side-chains with DNA base pairs. Mol. Phys. 2011, 109, 1995–2008. [Google Scholar] [CrossRef]
  81. Bavafa, S.; Behjatmanesh-Ardakani, R.; Mashhadi, F.F. Density functional study of the nitrosamine–formic acid and nitrosamine–formamide interactions. J. Mol. Struct. Theochem 2010, 960, 15–21. [Google Scholar] [CrossRef]
  82. de Silva, P.; Corminboeuf, C. Simultaneous Visualization of Covalent and Noncovalent Interactions Using Regions of Density Overlap. J. Chem. Theory Comput. 2014, 10, 3745–3756. [Google Scholar] [CrossRef]
  83. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  84. Lu, T.; Chen, F. Quantitative analysis of molecular surface based on improved Marching Tetrahedra algorithm. J. Molec. Graph. Model. 2012, 38, 314–323. [Google Scholar] [CrossRef]
  85. (1990–2004) Copyright MPI für Festkoerperforschung Stuttgart (1997–2001), CPMD Version 4.3-4610; Copyright IBM Corp.
  86. CCDC Structural Database. Available online: https://www.ccdc.cam.ac.uk/ (accessed on 7 May 2021).
  87. Schlegel, H.B. Estimating the hessian for gradient-type geometry optimizations. Theor. Chem. Acc. 1984, 66, 333–340. [Google Scholar] [CrossRef]
  88. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Handy, N.C.; Cohen, A.J. Left-right correlation energy. Mol. Phys. 2001, 99, 403–412. [Google Scholar] [CrossRef]
  90. Lee, C.; Yang, W.; 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] [PubMed] [Green Version]
  91. Troullier, N.; Martins, J.L. Efficient pseudopotentials for plane-wave calculations. Phys. Rev. B 1991, 43, 1993–2006. [Google Scholar] [CrossRef]
  92. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef]
  93. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef] [Green Version]
  94. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [Green Version]
  95. Bader, R. Atoms in Molecules: A Quantum Theory; International Ser. of Monogr. on Chem.; Clarendon Press: Oxford, UK, 1994. [Google Scholar]
  96. Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158–6170. [Google Scholar] [CrossRef]
  97. 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]
  98. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16 Revision C.01; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  99. Keith, T.A.; Gristmill, T. AIMAll (Version 19.10.12); Software: Overland Park, KS, USA, 2019. [Google Scholar]
  100. Gálvez, O.; Gómez, P.C.; Pacios, L.F. Variation with the intermolecular distance of properties dependent on the electron density in hydrogen bond dimers. J. Chem. Phys. 2001, 115, 11166–11184. [Google Scholar] [CrossRef]
  101. Pacios, L.F. Change with the Intermolecular Distance of Electron Properties of Hydrogen Bond Dimers at Equilibrium and Non-equilibrium Geometries. Struct. Chem. 2005, 16, 223–241. [Google Scholar] [CrossRef]
  102. Brovarets’, O.O.; Voiteshenko, I.S.; Hovorun, D.M. Physico-chemical profiles of the wobble–Watson–Crick G*·2AP(w)–G·2AP(WC) and A·2AP(w)–A*·2AP(WC) tautomerisations: A QM/QTAIM comprehensive survey. Phys. Chem. Chem. Phys. 2018, 20, 623–636. [Google Scholar] [CrossRef] [PubMed]
  103. Pichierri, F. Theoretical insights into the nature of the bonding between carbon monoxide and iron(II) phthalocyanine: How do QTAIM descriptors change as a function of the Fe–CO distance? Chem. Phys. Lett. 2022, 804, 139901. [Google Scholar] [CrossRef]
  104. Humphrey, W.; Dalke, A.; Schulten, K. VMD–Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  105. Momma, K.; Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 2011, 44, 1272–1276. [Google Scholar] [CrossRef]
  106. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera-A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef]
  107. Williams, T.; Kelley, C.; Lang, R.; Kotz, D.; Campbell, J.; Elber, G.; Woo, A. Gnuplot 5.8.2: An Interactive Plotting Program. Available online: http://www.gnuplot.info (accessed on 13 November 2022).
  108. Dennington, R.; Keith, T.A.; Millam, J.M. GaussView Version 6; Semichem Inc.: Shawnee, KS, USA, 2019. [Google Scholar]
Ijms 24 01542 g001 550
Figure 1. Structure with CCDC code and relevant atoms numbering of the protonated forms of (1) 1,8-bis(dimethylamino)naphthalene (TAPCES); (2) 1,8-bis(dimethylamino)-4,5-dimethoxynaphthalene (RISBEA); (3) 4-amino-1,8-bis(dimethylamino)naphthalene (XUCKAH); (4) 4-nitro-1,8-bis(dimethylamino)naphthalene (ZOSKEX). The dotted lines indicate the presence of intramolecular hydrogen bonds.
Figure 1. Structure with CCDC code and relevant atoms numbering of the protonated forms of (1) 1,8-bis(dimethylamino)naphthalene (TAPCES); (2) 1,8-bis(dimethylamino)-4,5-dimethoxynaphthalene (RISBEA); (3) 4-amino-1,8-bis(dimethylamino)naphthalene (XUCKAH); (4) 4-nitro-1,8-bis(dimethylamino)naphthalene (ZOSKEX). The dotted lines indicate the presence of intramolecular hydrogen bonds.
Ijms 24 01542 g001
Ijms 24 01542 g002 550
Figure 2. Hirshfeld surfaces obtained for the studied “proton sponge” crystals with emphasis on intramolecular hydrogen bonding. Color coding: red—oxygen; blue–nitrogen; gray–carbon; white—hydrogen; green—chlorine; orange—bromine. (1)–(4) denote the studied compounds.
Figure 2. Hirshfeld surfaces obtained for the studied “proton sponge” crystals with emphasis on intramolecular hydrogen bonding. Color coding: red—oxygen; blue–nitrogen; gray–carbon; white—hydrogen; green—chlorine; orange—bromine. (1)–(4) denote the studied compounds.
Ijms 24 01542 g002
Ijms 24 01542 g003 550
Figure 3. Hirshfeld surfaces obtained for the studied “proton sponge” crystals with emphasis on intermolecular interactions. Color coding: red—oxygen; blue—nitrogen; gray—carbon; white—hydrogen; green—chlorine; orange—bromine. (1)–(4) indicate the studied compounds.
Figure 3. Hirshfeld surfaces obtained for the studied “proton sponge” crystals with emphasis on intermolecular interactions. Color coding: red—oxygen; blue—nitrogen; gray—carbon; white—hydrogen; green—chlorine; orange—bromine. (1)–(4) indicate the studied compounds.
Ijms 24 01542 g003
Ijms 24 01542 g004 550
Figure 4. Time evolution of the intramolecular hydrogen bond metric parameters. Gas phase CPMD results at 100 K for DMANH + and its derivatives.
Figure 4. Time evolution of the intramolecular hydrogen bond metric parameters. Gas phase CPMD results at 100 K for DMANH + and its derivatives.
Ijms 24 01542 g004
Ijms 24 01542 g005 550
Figure 5. Time evolution of the intramolecular hydrogen bond metric parameters. Gas phase CPMD results at 300 K for DMANH + and its derivatives.
Figure 5. Time evolution of the intramolecular hydrogen bond metric parameters. Gas phase CPMD results at 300 K for DMANH + and its derivatives.
Ijms 24 01542 g005
Ijms 24 01542 g006 550
Figure 6. Time evolution of the intramolecular hydrogen bond metric parameters. Crystalline phase CPMD results at 100 K for DMANH + and its derivatives.
Figure 6. Time evolution of the intramolecular hydrogen bond metric parameters. Crystalline phase CPMD results at 100 K for DMANH + and its derivatives.
Ijms 24 01542 g006
Ijms 24 01542 g007 550
Figure 7. Time evolution of the hydrogen bond metric parameters. Crystalline-phase CPMD results at 300 K for DMANH + and its derivatives.
Figure 7. Time evolution of the hydrogen bond metric parameters. Crystalline-phase CPMD results at 300 K for DMANH + and its derivatives.
Ijms 24 01542 g007
Ijms 24 01542 g008 550
Figure 8. Free energy surface for proton motion in the compounds (1)–(4)-comparison of gas phase simulations at 300 K with classical (CPMD) and quantum (PIMD) nuclear regime. The axes are the N⋯H distances—the horizontal axis for the donor-H and the vertical axis for the acceptor-H distances in Å.
Figure 8. Free energy surface for proton motion in the compounds (1)–(4)-comparison of gas phase simulations at 300 K with classical (CPMD) and quantum (PIMD) nuclear regime. The axes are the N⋯H distances—the horizontal axis for the donor-H and the vertical axis for the acceptor-H distances in Å.
Ijms 24 01542 g008
Ijms 24 01542 g009 550
Figure 9. Free energy surface for proton motion in the compounds (1)–(4). Two-dimensional slices across the three-dimensional volumetric data resulting from the CPMD-based gas phase metadynamics. The N⋯N value for the slice is 2.6 Å.
Figure 9. Free energy surface for proton motion in the compounds (1)–(4). Two-dimensional slices across the three-dimensional volumetric data resulting from the CPMD-based gas phase metadynamics. The N⋯N value for the slice is 2.6 Å.
Ijms 24 01542 g009
Ijms 24 01542 g010 550
Figure 10. QTAIM molecular diagrams for compound (4) at 100 K and 300 K. Green dots represent bond critical points (BCPs) while red – ring critical points (RCPs). Values of electron densities at the left and right N-H BCPs are given in a.u. Color coding: red—oxygen; blue—nitrogen; gray—carbon; white—hydrogen.
Figure 10. QTAIM molecular diagrams for compound (4) at 100 K and 300 K. Green dots represent bond critical points (BCPs) while red – ring critical points (RCPs). Values of electron densities at the left and right N-H BCPs are given in a.u. Color coding: red—oxygen; blue—nitrogen; gray—carbon; white—hydrogen.
Ijms 24 01542 g010
Ijms 24 01542 g011 550
Figure 11. DORI isosurfaces (green spheres represent non-covalent interaction regions, blue s spheres represent covalent ones, and red spheres represent steric repulsions) in the RDG 0.5 a.u. isovalue for compound (4) at 100 K at 0 and 24 ps.
Figure 11. DORI isosurfaces (green spheres represent non-covalent interaction regions, blue s spheres represent covalent ones, and red spheres represent steric repulsions) in the RDG 0.5 a.u. isovalue for compound (4) at 100 K at 0 and 24 ps.
Ijms 24 01542 g011
Table
Table 1. The interatomic contacts present in the crystal structures of the investigated “proton sponges”. Symbol “+” means the presence of the contact, while “−” means the absence of the contact.
Table 1. The interatomic contacts present in the crystal structures of the investigated “proton sponges”. Symbol “+” means the presence of the contact, while “−” means the absence of the contact.
CompoundContact
O⋯HCl⋯HC⋯HH⋯BrC-H⋯π
(1)++−−−
(2)+−+−+
(3)+−+++
(4)+−−−−
Table
Table 2. Proton possession percentage during the CPMD simulation of the DMANH + and its derivatives. The reference nitrogen atoms are those at the position 1 of the aromatic system, i.e., the left amine functional group in Figure 1. The simulations were carried out using the PBE functional.
Table 2. Proton possession percentage during the CPMD simulation of the DMANH + and its derivatives. The reference nitrogen atoms are those at the position 1 of the aromatic system, i.e., the left amine functional group in Figure 1. The simulations were carried out using the PBE functional.
CompoundProton Possession (%)
Gas PhaseCrystalline Phase
100 K300 K100 K300 K
(1)50.248.984.460.8
(2)51.850.829.748.4
(3)91.574.255.472.5
(4)18.337.810.936.8
Table
Table 3. Cambridge Crystallographic Data Centre (CCDC) code, number, and experimental unit cell data for the investigated “proton sponges”.
Table 3. Cambridge Crystallographic Data Centre (CCDC) code, number, and experimental unit cell data for the investigated “proton sponges”.
NumberCCDC Code and NumberUnit Cell DataRef.
1TAPCES (1266338)orthorhombic, a = 11.363 Å, b = 16.676 Å,[47]
c = 20.307 Å, Z = 4
2RISBEA (122006)tetragonal, a = 11.410 Å, b = 11.410 Å,[48]
c = 13.126 Å, Z = 4
3XUCKAH (170012)monoclinic, a = 11.182 Å, b = 14.236 Å,[49]
c = 19.935 Å, Z = 8
4ZOSKEX (1315231)triclinic, a= 7.986 Å, b = 12.463 Å,[50]
c = 8.663 Å, Z = 2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kizior, B.; Michalczyk, M.; Panek, J.J.; Zierkiewicz, W.; Jezierska, A. Unraveling the Nature of Hydrogen Bonds of “Proton Sponges” Based on Car-Parrinello and Metadynamics Approaches. Int. J. Mol. Sci. 2023, 24, 1542. https://doi.org/10.3390/ijms24021542

AMA Style

Kizior B, Michalczyk M, Panek JJ, Zierkiewicz W, Jezierska A. Unraveling the Nature of Hydrogen Bonds of “Proton Sponges” Based on Car-Parrinello and Metadynamics Approaches. International Journal of Molecular Sciences. 2023; 24(2):1542. https://doi.org/10.3390/ijms24021542

Chicago/Turabian Style

Kizior, Beata, Mariusz Michalczyk, Jarosław J. Panek, Wiktor Zierkiewicz, and Aneta Jezierska. 2023. "Unraveling the Nature of Hydrogen Bonds of “Proton Sponges” Based on Car-Parrinello and Metadynamics Approaches" International Journal of Molecular Sciences 24, no. 2: 1542. https://doi.org/10.3390/ijms24021542

APA Style

Kizior, B., Michalczyk, M., Panek, J. J., Zierkiewicz, W., & Jezierska, A. (2023). Unraveling the Nature of Hydrogen Bonds of “Proton Sponges” Based on Car-Parrinello and Metadynamics Approaches. International Journal of Molecular Sciences, 24(2), 1542. https://doi.org/10.3390/ijms24021542

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

Article Metrics

Back to TopTop