Characterizing Hydrogen Bonding and Proton Transfer in 2:1 FH:NH3 and FH:Collidine Complexes through One- and Two-Bond Spin−Spin Coupling Constants across Hydrogen Bonds

A combined molecular dynamics simulation and quantum mechanics study on mercaptopurine interaction with the cucurbit [6,7] urils: Analysis of electronic structure

In the current study, the probability of complex formation between mercaptopurine drug with cucurbit[6]urils and cucurbit[7]urils has been investigated. The calculations for geometry optimization of complexes have been carried out by means of DFT (B3LYP), DFT-D (B3LYP-D) and M06-2X methods. The Atoms In Molecules (AIM), Natural Bond Orbital (NBO), NMR, the density of states (DOSs) and frontier molecular orbital (MO) analyses have been done on the inclusion complexes. In addition, the UV-Vis spectra of the first eight states have been obtained by CAM-B3LYP/TD-DFT calculation. The obtained results of the complexation process reveal that CB[7]-DRG complexes are more favorable than that of CB[6]-DRG interactions. Furthermore, our theoretical results show that configurations III and I are the most stable configurations related to the CB[6]/DRG and CB[7]/DRG interactions, respectively. The positive ∇²_ρ(r) and HC values at the bond critical points indicate that exist the weak H-bonds between CB[6] and CB[7] with H atoms of the drug molecule. The obtained negative binding energy values of CB[7]-DRG interaction in solution phase show the stability of these complexes in the aqueous medium. Also, all of the observed parameters of molecular dynamics simulation such as the number of contacts, hydrogen bonding, center-of-mass distance and van der Waals energy values confirm the encapsulation of mercaptopurine molecule inside the cucurbit[7]urils cavity at about 3.2ns.

Dynamic and static behavior of the H...π and E...π interactions in EH₂ adducts of benzene π-system (E = O, S, Se and Te), elucidated by QTAIM dual functional analysis

Dynamic and static behavior of the interactions in the EH2 adducts of a benzene π-system (E = O, S, Se and Te) is elucidated by applying QTAIM-DFA (QTAIM dual functional analysis). Two types of H-*-π and E-*-π interactions are detected in the adducts, where the asterisk (*) emphasizes the existence of the bond critical point (BCP) on the interaction in question. Total electron energy densities Hb(rc) are plotted versus Hb(rc) -Vb(rc)/2 [=(ℏ(2)/8m)∇(2)ρb(rc)] at BCPs in QTAIM-DFA, where Vb(rc) are the potential energy densities at BCPs. Data from the fully optimized structures are analyzed by polar (R, θ) coordinate representation. Each plot for an interaction, containing data from the perturbed structures with those of the fully optimized one, shows a specific curve, which provides important information. The plot is expressed by (θp, κp): θp corresponds to the tangent line for the plot and κp is the curvature. θ and θp are measured from the y-axis and y-direction, respectively. Moreover, (R, θ) corresponds to the static nature, (θp, κp) represents the dynamic nature of interactions. While θ classifies the interaction in question, θp characterizes it. Both values are less than 90° for all H-*-π and E-*-π interactions examined in this study; therefore, they are all classified by the pure closed-shell interactions and predicted to have the character of vdW nature. However, it is suggested that E-*-π has the nature of the stronger interaction than the case of H-*-π for dynamic behavior in the same species evaluated at the MP2 and M06-2X levels. The nature of the interactions is well analyzed and specified by applying QTAIM-DFA.

Dynamic behavior of hydrogen bonds from pure closed shell to shared shell interaction regions elucidated by AIM dual functional analysis

The dynamic behavior of hydrogen bonds (HBs) was clarified for the wide range of interactions applying AIM dual functional analysis. Plots of H(b)(r(c)) versus H(b)(r(c)) - V(b)(r(c))/2 are analyzed in the polar (R, θ) representation, where H(b)(r(c)) and V(b)(r(c)) are total electron and potential energy densities at bond critical points, respectively, for the fully optimized structures. Data of the fully optimized structure and four perturbed ones around it are plotted for each interaction, which give a specific curve. The curve is analyzed by (θ(p), κ(p)): θ(p) corresponds to the tangent line from the y-direction and κ(p) is the curvature. Whereas (R, θ) correspond to the static nature, (θ(p), κ(p)) represent the dynamic nature of interactions. Indeed, HBs can be classified only by one parameter of θ, but θ(p) supplies more information necessary for better understanding of HBs. Although H(2)Se-*-HSeH and H(3)N-*-HNH(2) show the nature of pure CS (closed shell) of the vdW-type, H(2)S-*-HSH and H(2)O-*-HOH contain the nature of pure CS other than the vdW-type (HB-typical). The regular CS nature is observed for B-*-HF (B = HF, H(2)Se, H(2)S, H(2)O, and H(2)C═O). The HF-*-HF interaction is described as HB-typical, whereas others are by CT(MC)-type. The nature of H(3)N-*-HX (X = F, Cl, Br) is regular CS of the CT(TBP)-type. HBs in charged species, such as [HOH-*-OH](-) and [H(2)O-*-H-*-OH(2)], show the weak covalent nature of SS (shard shell). The dynamic behavior of HBs helps us to understand HBs in more detail, in addition to the static behavior.

Theoretical study of bifurcated hydrogen bonding effects on the 1J(N,H), 1hJ(N,H), 2hJ(N,N) couplings and 1H, 15N shieldings in model pyrroles

According to the density functional theory calculations, the X...H...N (X=N, O) intramolecular bifurcated (three-centered) hydrogen bond with one hydrogen donor and two hydrogen acceptors causes a significant decrease of the (1h)J(N,H) and (2h)J(N,N) coupling constants across the N-H...N hydrogen bond and an increase of the (1)J(N,H) coupling constant across the N-H covalent bond in the 2,5-disubstituted pyrroles. This occurs due to a weakening of the N-H...N hydrogen bridge resulting in a lengthening of the N...H distance and a decrease of the hydrogen bond angle at the bifurcated hydrogen bond formation. The gauge-independent atomic orbital calculations of the shielding constants suggest that a weakening of the N-H...N hydrogen bridge in case of the three-centered hydrogen bond yields a shielding of the bridge proton and deshielding of the acceptor nitrogen atom. The atoms-in-molecules analysis shows that an attenuation of the (1h)J(N,H) and (2h)J(N,N) couplings in the compounds with bifurcated hydrogen bond is connected with a decrease of the electron density rho(H...N) at the hydrogen bond critical point and Laplacian of this electron density nabla(2)rho(H...N). The natural bond orbital analysis suggests that the additional N-H...X interaction partly inhibits the charge transfer from the nitrogen lone pair to the sigma*(N-H) antibonding orbital across hydrogen bond weakening of the (1h)J(N,H) and (2h)J(N,N) trans-hydrogen bond couplings through Fermi-contact mechanism. An increase of the nitrogen s-character percentage of the N-H bond in consequence of the bifurcated hydrogen bonding leads to an increase of the (1)J(N,H) coupling constant across the N-H covalent bond and deshielding of the hydrogen donor nitrogen atom.

Comparative analysis of hydrogen bonding with participation of the nitrogen, oxygen and sulfur atoms in the 2(2'-heteroaryl)pyrroles and their trifluoroacetyl derivatives based on the 1H, 13C, 15N spectroscopy and DFT calculations

The N-H...X (X = N,O,S) intramolecular hydrogen bond in the series of 2(2'-heteroaryl)pyrroles and their trifluoroacetyl derivatives is examined by the (1)H, (13)C, (15)N spectroscopy and density functional theory (DFT) calculations. The influence of the hydrogen bond on coupling and shielding constants is considered. It is shown that the N-H...N intramolecular hydrogen bond causes a larger increase in the absolute size of the (1)J(N,H) coupling constant and a larger deshielding of the bridge proton than the N-H...O hydrogen bond. The effect of the N-H...S interaction on the (1)J(N,H) coupling constant and the shielding of the bridge proton is small. The NMR parameter changes in the series of the 2(2'-heteroaryl)pyrroles due to N-H...X hydrogen bond and the series of the 1-vinyl-2-(2'-heteroaryl)-pyrroles due to C-H...X hydrogen bond have the same order. The proximity of the nitrogen, oxygen or sulfur lone pair to the F...H hydrogen bridge quenches the trans-hydrogen bond spin-spin couplings (1h)J(F,H-1) and (2h)J(F,N).

Systematic ab initio study of 15N-15N and 15N-1H spin-spin coupling constants across N-H+-N hydrogen bonds: predicting N-N and N-H coupling constants and relating them to hydrogen bond type

A systematic ab initio EOM-CCSD study of 15N-15N and 15N-1H spin-spin coupling constants has been carried out for a series of complexes formed from 11 nitrogen bases with experimentally measured proton affinities. When these complexes are arranged in order of increasing proton affinity of the proton-acceptor base and, for each proton acceptor, increasing order of proton affinity of the protonated N-H donor, trends in distances and signs of coupling constants are evident that are indicative of the nature of the hydrogen bond. All two-bond spin-spin coupling constants (2hJ(N-N)) are positive and decrease as the N-N distance increases. All one-bond N-H coupling constants (1J(N-H)) are negative (1K(N-H) are positive). 1J(N-H) is related to the N-H distance and the hybridization of the donor N atom. One-bond H...N coupling constants (1hJ(H-N)) are positive (1hK(H-N) are negative) for traditional hydrogen bonds, but 1hJ(H-N) becomes negative when the hydrogen bond acquires sufficient proton-shared character. The N-N and H...N distances at which 1hJ(H-N) changes sign are approximately 2.71 and 1.62 A, respectively. Predictions are made of the values of 2hJ(N-N) and 1J(N-H), and the signs of 1hJ(H-N), for those complexes that are too large for EOM-CCSD calculations.

Complexes with N-H(+)-P hydrogen bonds: structures, binding energies, and spin-spin coupling constants

Ab-initio MP2/aug'-cc-pVTZ calculations have been performed to determine the structures and binding energies of proton-bound complexes stabilized by N-H+-P hydrogen bonds and to investigate the nature of the proton-transfer coordinate in these systems. Double minima are found only when the difference between the protonation energies of the N and P bases is less than about 4 kcal/mol. The isomer in which the protonated nitrogen base is the donor lies lower on the potential surface and also has a greater binding energy relative to the corresponding isolated monomers. Equation-of-motion coupled cluster singles and doubles (EOM-CCSD) calculations have been employed to obtain one- and two-bond spin-spin coupling constants across these hydrogen bonds. Two-bond coupling constants (2h)J(N-P) correlate with N-P distances, irrespective of whether the donor ion is N-H+ or P-H+. One-bond coupling constants (1)J(N-H) and (1h)J(H-P) for complexes stabilized by N-H+...P hydrogen bonds correlate with corresponding distances, but similar correlations are not found for (1)J(P-H) and (1h)J(H-N) for complexes with P-H+...N hydrogen bonds. Negative values of (1h)K(H-N) and (1h)K(H-P) indicate that the hydrogen bonds in these complexes are traditional. Comparisons are made with complexes stabilized by N-H+-N and P-H+-P hydrogen bonds.

Probing P−H+−P Hydrogen Bonds: Structures, Binding Energies, and Spin−Spin Coupling Constants

Ab initio MP2/aug'-cc-pVTZ calculations have been performed to determine the structures and binding energies of 22 open and 3 cyclic complexes formed from the sp2 [H(2)C=PH and HP=PH (cis and trans)] and sp3 [PH2(CH3) and PH3] hybridized phosphorus bases and their corresponding protonated ions. EOM-CCSD calculations have been carried out to obtain (31)P-(31)P and (31)P-(1)H coupling constants across P-H+-P hydrogen bonds. Two equilibrium structures with essentially linear hydrogen bonds have been found along the proton-transfer coordinate, except for complexes with P(CH3)H3+ as the proton donor to the sp2 bases. Although the isomer having the conjugate acid of the stronger base as the proton donor lies lower on the potential energy surface, it has a smaller binding energy relative to the corresponding isolated monomers than the isomer with the conjugate acid of the weaker base as the donor. The hydrogen bond of the latter has increased proton-shared character. All of the complexes are stabilized by traditional hydrogen bonds, as indicated by positive values of the reduced coupling constants (2h)K(P-P) and (1)K(P-H), and negative values of (1h)K(H-P). (2h)J(P-P) correlates with the P-P distance, a correlation determined primarily by the nature of the proton donor. For open complexes, (1)J(P-H) always increases relative to the isolated monomer, while (1h)J(H-P) is relatively small and negative. (2h)J(P-P) values are quite large in open complexes, but are much smaller in cyclic complexes in which the P-H+-P hydrogen bonds are nonlinear. Thus, experimental measurements of (2h)J(P-P) should be able to differentiate between open and cyclic complexes.

Computed coupling constants in X(CH3)nH(4-n) moieties where X = 13C and 15N+, and n = 0-4: comparisons with experimental data

Seventy-three unique spin-spin coupling constants have been analyzed for the ten species in the two series X(CH3)nH(4-n), where the central atom X is 13C or 15N+. Thirty-seven experimental values have been obtained from the literature, and several new coupling constants have been measured for the methyl-substituted ammonium ions. Both DFT with the B3LYP functional and ab initio EOM-CCSD calculations have been carried out on these same systems. Coupling constants computed by these two methods are in agreement with experimental values. Some problems related to coupling constants for the cationic ammonium systems have been resolved when these were recomputed at EOM-CCSD for complexes in which NH4+ is hydrogen-bonded to H2O molecules.

The C60 Formation Puzzle “Solved”: QM/MD Simulations Reveal the Shrinking Hot Giant Road of the Dynamic Fullerene Self-Assembly Mechanism

The dynamic self-assembly mechanism of fullerene molecules is an irreversible process emerging naturally under the nonequilibrium conditions of hot carbon vapor and is a consequence of the interplay between the dynamics and chemistry of polyyne chains, pi-conjugation and corresponding stabilization, and the dynamics of hot giant fullerene cages. In this feature article we briefly present an overview of experimental findings and past attempts to explain fullerene formation and show in detail how our recent quantum chemical molecular dynamics simulations of the dynamics of carbon vapor far from thermodynamic equilibrium have assisted in the discovery of the combined size-up/size-down "shrinking hot giant" road that leads to the formation of buckminsterfullerene C60, C70, and larger fullerenes. This formation mechanism is the first reported case of order created out of chaos where a distinct covalent bond network of an entire molecule is spontaneously self-assembled to a highly symmetric structure and fully explains the fullerene formation process consistently with all available experimental observations a priori. Experimental evidence suggests that it applies universally to all fullerene formation processes irrespective of the carbon source.

Ab Initio Study of Hydrogen Bonding and Proton Transfer in 3:1 FH:NH3 and FH:Collidine Complexes: Structures and One- and Two-Bond Coupling Constants across Hydrogen Bonds

Ab initio EOM-CCSD calculations have been performed on 3:1 FH:NH3 complexes at their own optimized MP2/6-31+G(d,p) geometries and at the optimized geometries in the hydrogen-bonding regions of corresponding 3:1 FH:collidine complexes. The isolated gas-phase equilibrium 3:1 FH:NH3 complex has an open structure with a proton-shared Fa-Ha-N hydrogen bond, while the isolated equilibrium 3:1 FH:collidine complex has a perpendicular structure with an Fa-Ha-N hydrogen bond that is on the ion-pair side of proton-shared. The Fa-N coupling constant ((2h)J(Fa-N)) for the equilibrium 3:1 FH:NH3 complex is large and negative, consistent with a proton-shared Fa-Ha-N hydrogen bond; (2h)JFb-Fa is positive, reflecting a short Fb-Fa distance and partial proton transfer from Fb to Fa across the Fb-Hb-Fa hydrogen bond. In contrast, (2h)JFa-N has a smaller absolute value and (2h)JFb-Fa is greater for the 3:1 FH:NH3 complex at the equilibrium 3:1 FH:collidine geometry, consistent with the structural characteristics of the Fa-Ha-N and Fb-Hb-Fa hydrogen bonds. Coupling constants computed at proton-transferred 3:1 FH:collidine perpendicular geometries are consistent with experimental coupling constants for the 3:1 FH:collidine complex in solution and indicate that the role of the solvent is to promote further proton transfer from Fa to N across the Fa-Ha-N hydrogen bond, and from Fb to Fa across the two equivalent Fb-Hb-Fa hydrogen bonds. The best correlations between experimental and computed coupling constants are found for complexes with perpendicular proton-transferred structures, one having the optimized geometry of a 3:1 FH:collidine complex at an Fa-Ha distance of 1.80 A, and the other at the optimized 3:1 FH:collidine geometry with distances derived from the experimental coupling constants. These calculations provide support for the proposed perpendicular structure of the 3:1 FH:collidine complex as the structure which exists in solution.