DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

Non-covalent interactions – QTAIM and NBO analysis

MP2(full)/6-311++G(3df,3pd) calculations were carried out on complexes linked through various non-covalent Lewis acid – Lewis base interactions. These are: hydrogen bond, dihydrogen bond, hydride bond and halogen bond. The quantum theory of ´atoms in molecules´ (QTAIM) as well as the natural bond orbitals (NBO) method were applied to analyze properties of these interactions. It was found that for the A-H…B hydrogen bond as well as for the A-X…B halogen bond (X designates halogen) the complex formation leads to the increase of s-character in the A-atom hybrid orbital aimed toward the H or X atom. In opposite, for the A…H-B hydride bond, where the H-atom possesses negative charge, the decrease of s-character in the B-atom orbital is observed. All these changes connected with the redistribution of the electron charge being the effect of the complex formation are in line with Bent´s rule. The numerous correlations between energetic, geometrical, NBO and QTAIM parameters were also found.

Ab Initio Quantum Mechanics Analysis of Imidazole C-H…O Water Hydrogen Bonding and a Molecular Mechanics Forcefield Correction

Abstract

While it is well established that classical hydrogen bonds play an important role in enzyme structure, function and dynamics, the role of weaker, but ‘activated’ C-H donor hydrogen bonds is poorly understood. The most important such case involves histidine which often plays a direct role in enzyme catalysis and possesses the most acidic C-H donor group of the standard amino acids. In the present study, we obtained optimized geometries and hydrogen bond interaction energies for C-H…O hydrogen bonded complexes between methane, ethylene, benzene, acetylene, and imidazole with water at the MP2-FC/6-31++G(2d,2p) and MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) levels of theory. A strong linear relationship is obtained between the stability of the various hydrogen bonded complexes and both separation distances for H…0 and C—O. In general, these calculations indicate that C-H…0 interactions can be classified as hydrogen bonding interactions, albeit significantly weaker than the classical hydrogen bonds, but significantly stronger than just van der Waals interactions. For instance, while the electronic energy of stabilization at the MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) level of theory of a water C-H…O water hydrogen bond is 4.36 kcal/mol more stable than the methane C-H…O water interaction, the water-water hydrogen bond is only 2.06 kcal/mol more stable than the imidazole C_e?H…O water hydrogen bond. Neglecting this latter hydrogen bonding interaction is obviously unacceptable. We next compare the potential energy surfaces for the imidazole C_e?H…O water and imidazole N_d?H…O hydrogen bonded complexes computed at the MP2/6-31++G(2d,2p) level of theory with the potential energy surface computed using the AMBER molecular mechanics program and forcefields. While the Weiner et al and Cornell et al AMBER forcefields reasonably account for the imidazole N-H…O water interaction, these forcefields do not adequately account for the imidazole C_e?H…O water hydrogen bond. A forcefield modification is offered that results in excellent agreement between the ab initio and molecular mechanics geometry and energy for this C-H…O hydrogen bonded complex.     

Theoretical investigations into the intermolecular hydrogen-bonding interactions of <Emphasis Type="Italic">N</Emphasis>-(hydroxymethyl)acetamide dimers

The structures of the N-(hydroxymethyl)acetamide (model molecule of ceramide) dimers have been fully optimized at B3LYP/6–311++G** level. The intermolecular hydrogen bonding interaction energies have been calculated using the B3LYP/6–311++G**, B3LYP/6–311++G(2df,2p), MP2(full)/6–311++G** and MP2(full)/6–311++G(2df,2p) methods, respectively. The results show that the O–H···O, N–H···O, O–H···N, and C–H···O hydrogen bonding interactions could exist in N-(hydroxymethyl)acetamide dimers, and the O–H···O, N–H···O, and O–H···N hydrogen bonding interactions could be stronger than C–H···O. The three-dimensional network structure formed by ceramide molecules through intermolecular hydrogen bonding interactions may be the main reason why the stratum corneum of skin could prevent foreign substances from entering our body, as is in accordance with the experimental results. The stability of hydrogen-bonding interactions follow the order of (a)?>?(b)?≈?(c)?>?(d)?>?(e)?≈?(f)?>?(g)?>?(h). The analyses of the energy decomposition, frequency, atoms in molecules (AIM), natural bond orbital (NBO), and electron density shift are used to further reveal the nature of the complex formation. In the range of 263.0–328.0 K, the complex is formed via an exothermic reaction, and the solvent with lower temperature and dielectric constant is favorable to this process.

Open image in new window

Graphical abstract The structures and the O–H···O=C, N–H···O=C and C–H···O=C H-bonding interactions in the N-(hydroxymethyl)acetamide (model molecule of ceramide) dimers were investigated using the B3LYP and MP2(full) methods.

     

A B3LYP and MP2(full) theoretical investigation into the strength of the C-NO2 bond upon the formation of the molecule-cation interaction between Na+ and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the molecule-cation interaction between Na⁺ and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced in comparison with that in the isolated nitrotriazole molecule upon the formation of molecule-cation interaction. The increment of the C–NO₂ bond dissociation energy (ΔBDE) correlated well with the molecule-cation interaction energy. Electron density shifts analysis showed that the electron density shifted toward the C-NO₂ bond upon complex formation, leading to the strengthened C-NO₂ bond and the possibly reduced explosive sensitivity.

Low-energy conformers of pamidronate and their intramolecular hydrogen bonds: a DFT and QTAIM study

Extensive DFT and ab initio calculations were performed to characterize the conformational space of pamidronate, a typical pharmaceutical for bone diseases. Mono-, di- and tri-protic states of molecule, relevant for physiological pH range, were investigated for both canonical and zwitterionic tautomers. Semiempirical PM6 method were used for prescreening of the single bond rotamers followed by geometry optimizations at the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels. For numerous identified low energy conformers the final electronic energies were determined at the MP2/6-311++G(2df,2p) level and corrected for thermal effects at B3LYP level. Solvation effects were also considered via the COSMO and C-PCM implicit models. Reasonable agreement was found between bond lengths and angle values in comparison with X-ray crystal structures. Relative equilibrium populations of different conformers were determined from molecular partition functions and the role of electronic, vibrational and rotational degrees of freedom on the stability of conformers were analyzed. For no level of theory is a zwitterionic structure stable in the gas-phase while solvation makes them available depending on the protonation state. Geometrically identified intramolecular hydrogen bonds were analyzed by QTAIM approach. All conformers exhibit strong inter-phosphonate hydrogen bonds and in most of them the alkyl-amine side chain is folded on the P-C-P backbone for further hydrogen bond formation.

A B3LYP and MP2(full) theoretical investigation into the strength of the C–NO2 bond upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced and the charge of nitro group turned more negative in complex in comparison with those in isolated nitrotriazole molecule. The increment of the C–NO₂ bond dissociation energies correlated well with the intermolecular H-bonding interaction energies. Electron density shifts analyses showed that the electron density shifted toward the C–NO₂ bond upon complex formation, leading to the strengthened C–NO₂ bond and the possibly reduced explosive sensitivity.

A B3LYP and MP2(full) theoretical investigation into the cooperativity effect between dihydrogen-bonding and H–M∙∙∙π (M = Li, Na, K) interactions among HF, MH with the π-electron donor C2H2, C2H4 or C6H6

The DFT-B3LYP/6-311++G(3df,2p) and MP2(full)/6-311++G(3df,2p) calculations were carried out on the binary complex formed by HM (M?=?Li, Na, K) and HF or the π-electron donor (C₂H₂, C₂H₄, C₆H₆), as well as the ternary system FH???HM???C₂H₂/C₂H₄/C₆H₆. The cooperativity effect between the dihydrogen-bonding and H–M???π interactions was investigated. The result shows that the equilibrium distances R _H???H and R _M???π in the ternary complex decrease and both the H???H and H–M???π interactions are strengthened when compared to the corresponding binary complex. The cooperativity effect of the dihydrogen bond on the H–M???π interaction is more pronounced than that of the M???π bond on the H???H interaction. Furthermore, the values of cooperativity effect follow the order of FH???HNa???π?>?FH???HLi???π?>?FH???HK???π and FH???HM???C₆H₆?>?FH???HM???C₂H₄?>?FH???HM???C₂H₂. The nature of the cooperativity effect was revealed by the analyses of the charge of the hydrogen atoms in H???H moiety, atom in molecule (AIM) and electron density shifts methods.

Theoretical studies of the interaction between enflurane and water

Increase of the atmospheric concentration of halogenated organic compounds is partially responsible for a change of the global climate. In this work we have investigated the interaction between halogenated ether and water, which is one of the most important constituent of the atmosphere. The structures of the complexes formed by the two most stable conformers of enflurane (a volatile anaesthetic) with one and two water molecules were calculated by means of the counterpoise CP-corrected gradient optimization at the MP2/6–311++G(d,p) level. In these complexes the CH…O_w hydrogen bonds are formed, with the H…O_w distances varying between 2.23 and 2.32 Å. A small contraction of the CH bonds and the blue shifts of the ν(CH) stretching vibrations are predicted. There is also a weak interaction between one of the F atoms and the H atom of water, with the H_w…F distances between 2.41 and 2.87 Å. The CCSD(T)/CBS calculated stabilization energies in these complexes are between ?5.89 and ?4.66 kcal?mol^?1, while the enthalpies of formation are between ?4.35 and ?3.22 kcal?mol^?1. The Cl halogen bonding between enflurane and water has been found in two complexes. The intermolecular (Cl···O) distance is smaller than the sum of the corresponding van der Waals radii. The CCSD(T)/CBS stabilization energies for these complexes are about ?2 kcal?mol^?1.

Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6–311++G(2df,2pd) basis set and CCSD/6–31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

Structural characterization of the (MeSH)4 potential energy surface

A random walk on the PES for (MeSH)₄ clusters produced 50 structural isomers held together by hydrogen-bonding networks according to calculations performed at the B3LYP/6–311++G** and MP2/6–311++G** levels. The geometric motifs observed are somewhat similar to those encountered for the methanol tetramer, but the interactions responsible for cluster stabilization are quite different in origin. Cluster stabilization is not related to the number of hydrogen bonds. Two distinct, well-defined types of hydrogen bonds scattered over a wide range of distances are predicted.

Spectroscopic investigations and hydrogen bond interactions of 8-aza analogues of xanthine, theophylline and caffeine: a theoretical study

The structure, spectral properties and the hydrogen bond interactions of 8-aza analogues of xanthine, theophylline and caffeine have been studied by using quantum chemical methods. The time-dependent density functional theory (TD-DFT) and the singly excited configuration interaction (CIS) methods are employed to optimize the excited state geometries of isolated 8-azaxanthine, 8-azatheophylline tautomers and 8-azacaffeine in both the gas and solvent phases. The solvent phase calculations are performed using the polarizable continuum model (PCM). The absorption and emission spectra are calculated using the time-dependent density functional theory (TD-DFT) method. The results from the TD-DFT calculations reveal that the excitation spectra are red shifted relative to absorption in aqueous medium. These changes in the transition energies are qualitatively comparable to the experimental data. The examination of molecular orbital reveals that the molecules with a small H→L energy gap possess maximum absorption and emission wavelength. The relative stability and hydrogen bonded interactions of mono and heptahydrated 8-azaxanthine, 8-azatheophylline tautomers and 8-azacaffeine have been studied using the density functional theory (DFT) and Møller Plesset perturbation theory (MP2) implementing the 6-311++G(d,p) basis set. The formation of strong N-H…O bond has resulted in the highest interaction energy among the monohydrates. Hydration does not show any significant impact on the stability of heptahydrated complexes. The atoms in molecule (AIM) and natural bonding orbital (NBO) analyses have been performed to elucidate the nature of the hydrogen bond interactions in these complexes.

A comparative theoretical investigation into the strength of the trigger-bond in the Na+, Mg2+ and HF complexes involving the nitro group of R–NO2 (R = –CH3, –NH2 and –OCH3) or the C = C bond of (E)-O2N–CH = CH–NO2

A comparative theoretical investigation into the change in strength of the trigger-bond upon formation of the Na⁺, Mg²⁺ and HF complexes involving the nitro group of RNO₂ (R?=? –CH₃, –NH₂, –OCH₃) or the C?=?C bond of (E)-O₂N–CH?=?CH–NO₂ was carried out using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. Except for the Mg²⁺?π system with (E)-O₂N–CH?=?CH–NO₂ (i.e., C₂H₂N₂O₄?Mg²⁺), the strength of the trigger-bond X–NO₂ (X?=?C, N or O) was enhanced upon complex formation. Furthermore, the increment of bond dissociation energy of the X–NO₂ bond in the Na⁺ complex was far greater than that in the corresponding HF system. Thus, the explosive sensitivity in the former might be lower than that in the latter. For C₂H₂N₂O₄?Mg²⁺, the explosive sensitivity might also be reduced. Therefore, it is possible that introducing cations into the structure of explosives might be more efficacious at reducing explosive sensitivity than the formation of an intermolecular hydrogen-bonded complex. AIM, NBO and electron density shifts analyses showed that the electron density shifted toward the X–NO₂ bond upon complex formation, leading to a strengthened X–NO₂ bond and possibly reduced explosive sensitivity.

A theoretical investigation into the cooperativity effect between the H∙∙∙O and H∙∙∙F– interactions and electrostatic potential upon 1:2 (F–:N-(Hydroxymethyl)acetamide) ternary-system formation

The cooperativity effects between the O/N–H???F^– anionic hydrogen-bonding and O/N–H???O hydrogen-bonding interactions and electrostatic potentials in the 1:2 (F^–:N-(Hydroxymethyl)acetamide (signed as “ha”)) ternary systems are investigated at the B3LYP/6-311++G** and MP2/6-311++G** levels. A comparison of the cooperativity effect in the “F^–???ha???ha” and “FH???ha^–???ha” systems is also carried out. The result shows that the increase of the H???O interaction energy in the O–H???O–H, N–H???O–H or N–H???O?=?C link is more notable than that in the O–H???O?=?C contact upon ternary-system formation. The cooperativity effect is found in the complex formed by the O/N–H???F^– and O/N–H???O interactions, while the anti-cooperativity effect is present in the system with only the O/N–H???F^– H-bond or the “FH???ha^–???ha” complex by the N^–???H–F contact. Atoms in molecules (AIM) analysis and shift of electron density confirm the existence of cooperativity. The most negative surface electrostatic potential (V _S,min ) correlates well with the interaction energy E′ _{int.(ha???F–)} and synergetic energy E _syn., respectively. The relationship between the change of V _S,min (i.e., ΔV _S,min ) and E _syn. is also found.

Prototropic tautomerism and basic molecular principles of hypoxanthine mutagenicity: an exhaustive quantum-chemical analysis

The molecular structures, relative stability order, and dipole moments of a complete family of 21 planar hypoxanthine (Hyp) prototropic molecular–zwitterionic tautomers including ylidic forms were computationally investigated at the MP2/6–311++G(2df,pd)//B3LYP/6–311++G(d,p) level of theory in vacuum and in three different surrounding environments: continuum with a low dielectric constant (??=?4) corresponding to a hydrophobic interface of protein–nucleic acid interactions, dimethylsulfoxide (DMSO), and water. The keto-N1HN7H tautomer was established to be the global minimum in vacuum and in continuum with ??=?4, while Hyp molecule exists as a mixture of the keto-N1HN9H and keto-N1HN7H tautomers in approximately equal amounts in DMSO and in water at T?=?298.15?K. We found out that neither intramolecular tautomerization by single proton transfer in the Hyp base, nor intermolecular tautomerization by double proton transfer in the most energetically favorable Hyp·Hyp homodimer (symmetry C _2h ), stabilized by two equivalent N1H…O6 H-bonds, induces the formation of the enol tautomer (marked with an asterisk) of Hyp with cis-oriented O6H hydroxyl group relative to neighboring N1C6 bond. We first discovered a new scenario of the keto–enol tautomerization of Hyp?·?Hyp homodimer (C _2h ) via zwitterionic near-orthogonal transition state (TS), stabilized by N1⁺H…N1^? and O6⁺H…N1^? H-bonds, to heterodimer Hyp^??·?Hyp (C _s ), stabilized by O6H…O6 and N1H…N1 H-bonds. We first showed that Hyp^??·?Thy mispair (C _s ), stabilized by O6H…O4, N3H…N1, and C2H…O2 H-bonds, mimicking Watson–Crick base pairing, converts to the wobble Hyp?·?Thy base pair (C _s ), stabilized by N3H…O6 and N1H…O2 H-bonds, via high- and low-energy TSs and intermediate Hyp?·?Thy^?, stabilized by O4H…O6, N1H…N3, and C2H…O2 H-bonds. The most energetically favorable TS is the zwitterionic pair Hyp⁺?·?Thy^? (C _s ), stabilized by O6⁺H…O4^?, O6⁺H…N3^?, N1⁺H…N3^?, and N1⁺H…O2^? H-bonds. The authors expressed and substantiated the hypothesis, that the keto tautomer of Hyp is a mutagenic compound, while enol tautomer Hyp^? does not possess mutagenic properties. The lifetime of the nonmutagenic tautomer Hyp^? exceeds by many orders the time needed to complete a round of DNA replication in the cell. For the first time purine–purine planar H-bonded mispairs containing Hyp in the anti-orientation with respect to the sugar moiety – Hyp?·?Ade _syn , Hyp?·?Gua^? _syn , and Hyp?·?Gua _syn , that closely resembles the geometry of the Watson–Crick base pairs, have been suggested as the source of transversions. An influence of the surrounding environment (??=?4) on the stability of studied complexes and corresponding TSs was estimated by means of the conductor-like polarizable continuum model. Electron-topological, structural, vibrational, and energetic characterictics of all conventional and nonconventional H-bonds in the investigated structures are presented. Presented data are key to understanding elementary molecular mechanisms of mutagenic action of Hyp as a product of the adenine deamination in DNA.     

DFT studies of the adsorption and dissociation of H2O on the Al13 cluster: origins of this reactivity and the mechanism for H2 release

A theoretical study of the chemisorption and dissociation pathways of water on the Al₁₃ cluster was performed using the hybrid density functional B3LYP method with the 6-311+G(d, p) basis set. The activation energies, reaction enthalpies, and Gibbs free energy of activation for the reaction were determined. Calculations revealed that the H₂O molecule is easily adsorbed onto the Al₁₃ surface, forming adlayers. The dissociation of the first H₂O molecule from the bimolecular H₂O structure via the Grotthuss mechanism is the most kinetically favorable among the five potential pathways for O–H bond breaking. The elimination of H₂ in the reaction of an H₂O molecule with a hydrogen atom on the Al cluster via the Eley–Rideal mechanism has a lower activation barrier than the elimination of H₂ in the reaction of two adsorbed H atoms or the reaction of OH and H. Following the adsorption and dissociation of H₂O, the structure of Al₁₃ is distorted to varying degrees.

Quantum-chemical investigation of tautomerization ways of Watson-Crick DNA base pair guanine-cytosine

A novel physico-chemical mechanism of the Watson-Crick DNA base pair Gua.Cyt tautomerization Gua.Cyt*<---->Gua.Cyt<---->Gua*.Cyt (mutagenic tautomers of bases are marked by asterisks) have been revealed and realized in a pathway of single proton transfer through two mutual isoenergetic transition states with Gibbs free energy of activation 30.4 and 30.6 kcal/mol and they are ion pairs stabilized by three (N2H...N3, N1H...N4- and O6+H...N4-) and five (N2H...O2, N1H...O2, N1H...N3, O6+H...N4- and 06+H...N4-) H-bonds accordingly. Stable base pairs Gua-Cyt* and Gua*.Cyt which dissociate comparably easy into monomers have acceptable relative Gibbs energies--12.9 and 14.3 kcal/mol--for the explanation of the nature of the spontaneous transitions of DNA replication. Results are obtained at the MP2/6-311++G(2df,pd)//B3LYP/6-31 1++G(d,p) level of theory in vacuum approach.     

Quantum-chemical investigation of the elementary molecular mechanisms of pyrimidine-purine transversions

Purine-purine mispairs of DNA (thus involving template base in anti-conformation along the glycosidic bond and base of the incoming nucleotide - in syn-conformation) leading to pyrimidine-purine "transversions"-type point mutations were revealed and characterized at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory in vacuum approach adequately modeling hydrophobic environment of the active centre of high-fidelity replicative DNA-polymerases.     

Cooperativity effect involving drug–DNA/RNA intermolecular interaction: A B3LYP-D3 and MP2 theoretical investigation on ketoprofen⋯cytosine⋯H2O system

In order to examine the origin of the drug action and design new DNA/RNA-targeted drugs, the cooperativity effect involving drug–DNA/RNA intermolecular interaction in ketoprofen?cytosine?H₂O ternary system were investigated by the B3LYP, B3LYP-D3, and MP2 methods with the 6-311++G(2d,p) basis set. The thermodynamic cooperativity was also evaluated at 310.15 K. The N–H?O, O–H?O, O–H?N, C–H?N, and C–H?O H bonds coexist in ternary complexes. The intermolecular interactions obtained by B3LYP-D3 are close to those calculated by MP2. The steric effects and van der Waals interactions have little influence on the cooperativity effects. The anti-cooperativity effect in ket?cyt?H₂O is far more notable than the cooperativity effect, and the stability of the cyclic structure with anti-cooperativity effect is higher than that of the linear structure with cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that, in the presence of H₂O, the anti-cooperativity effect plays a dominant role in the drug–DNA/RNA interaction, and the nature of the hydration in the binding of drugs to DNA/RNA bases is the H-bonding anti-cooperativity effect. Furthermore, the drug always links simultaneously with DNA/RNA base and H₂O, and only in this way can the biological activity of drugs play a role. In most cases, the enthalpy change is the major factor driving the cooperativity, as is different from most of biomacromolecule complexes.     

Density functional theory study of N-H...O, O-H...O and C-H...O hydrogen-bonding effects on the 14N and 2H nuclear quadrupole coupling tensors of N-acetyl-valine

Hydrogen-bonding effects in the crystalline structure of N-acetyl-valine, NAV, were studied using the (14)N and (2)H quadrupole coupling tensors via density functional theory. The calculations were carried out at the B3LYP level with the 6-311++G(d,p) and 6-311+G(d) basis sets. The theoretical quadrupole coupling components and their relative orientation in the molecular frame axes at the nitrogen site are compared to experimental values. This nucleus is involved in a rather strong intermolecular O=CNH...O=CNH hydrogen bond, r(N-H...O(1))=2.04 A and angleN-H...O(1)=171.53 degrees. A reasonably good agreement between the experimentally obtained (2)H quadrupole coupling tensors and the B3LYP/6-311++G(d,p) calculations is achievable only in molecular model where a complete hydrogen-bonding network is considered.