Investigation of intermolecular hydrogen bond interactions in crystalline l-Cysteine by DFT calculations of the oxygen-17, nitrogen-14, and hydrogen-2 EFG tensors and AIM analysis

A systematic computational study is carried out to investigate hydrogen bond (HB) interactions in the real crystalline structures of l-Cysteine at 30 and 298 K by density functional theory (DFT) calculations of electric field gradient (EFG) tensors at the sites of O-17, N-14, and H-2 nuclei. One-molecule (monomer) and nine-molecule (cluster) models of l-Cysteine are created by available crystal coordinates at both temperatures and the EFG tensors are calculated for both models to indicate the effect of HB interactions on the tensors. The calculated EFG tensors at the level of B3LYP and B3PW91 DFT methods and 6-311++G?? and cc-pVTZ basis sets are converted to those experimentally measurable nuclear quadrupole resonance (NQR) parameters i.e. quadrupole coupling constants (qcc) and asymmetry parameters (η_Q). The evaluated NQR parameters reveal that the EFG tensors of ¹⁷O, ¹⁴N, and ²H are influenced and show particular trends from monomer to the target molecule in the cluster due to the contribution of target molecule to classic N–H…O, and non-classic S–H…O and S–H…S types of HB interactions. On the other hand, atoms in molecules (AIM) analyses confirm the presence of HB interactions and rationalize the observed EFG trends. The results indicate different contribution of various nuclei to HB interactions in the cluster where O2 and N1 have major contributions. The EFG tensors as well as AIM analysis at the H6 site show that the N1-H6…O2 HB undergoes a significant change from 30 to 298 K where changes in other N–H…O interactions are almost negligible. There is a good agreement between the calculated ¹⁴N NQR parameters and reported experimental data.     

Investigation of C-H...O=C and N-H...OC hydrogen-bonding interactions in crystalline thymine by DFT calculations of O-17, N-14 and H-2 NQR parameters

A computational study at the level of density functional theory (DFT) was carried out to investigate C-H...O=C and N-H...O=C hydrogen-bonding interactions (HBs) in the real crystalline cluster of thymine by O-17, N-14 and H-2 calculated nuclear quadrupole resonance (NQR) parameters. To perform the calculations, a hydrogen-bonded pentameric cluster of thymine was created using X-ray coordinates where the hydrogen atoms positions are optimized and the electric field gradient (EFG) tensors were calculated for the target molecule. Additional EFG calculations were also performed for crystalline monomer and an optimized isolated gas-phase thymine. The calculated EFG tensors at the level of B3LYP and B3PW91 DFT methods and 6-311++G**and CC-pVTZ basis sets were converted to those experimentally measurable NQR parameters, quadrupole coupling constants and asymmetry parameters. The results reveal that because of strong contribution to N-H...O=C HBs, NQR parameters of O2, N1 and N3 undergo significant changes from monomer to the target molecule in cluster. Furthermore, the NQR parameters of O2 also undergo some changes because of non-classical C-H...O=C HBs.     

Theoretical investigation of hydrogen bonding effects on oxygen, nitrogen, and hydrogen chemical shielding and electric field gradient tensors of chitosan/HI salt

A density functional theory study has been carried out to calculate the (17)O, (15)N, (13)C, and (1)H chemical shielding as well as (17)O, (14)N, and (2)H electric field gradient tensors of chitosan/HI type I salt. These calculations were performed using the B3LYP functional and 6-311++G (d,p) and 6-31++G (d,p) basis sets. Calculated EFG and chemical shielding tensors were used to evaluate the (17)O, (14)N, and (2)H nuclear quadruple resonance, NQR, and (17)O, (15)N, (13)C, and (1)H nuclear magnetic resonance, NMR, parameters in the cluster model, which are in good agreement with the available experimental data. The difference in the isotropic shielding (sigma(iso)) and quadrupole coupling constant (C(Q)) between monomer and target molecule in the cluster was analyzed in detail. It was shown that both EFG and CS tensors are sensitive to hydrogen-bonding interactions, and calculating both tensors is an advantage. A different influence of various hydrogen bond types, N-Hcdots, three dots, centeredI, O-Hcdots, three dots, centeredI, and N-Hcdots, three dots, centeredO was observed on the calculated CS and EFG tensors. On the basis of this study, nitrogen and O-6 are the most important nuclei to confirm crystalline structure of chitosan/HI. These nuclei have large change in their CS and EFG tensors because of forming intermolecular hydrogen bonds. Moreover, the quantum chemical calculations indicated that the intermolecular hydrogen-bonding interactions play an essential role in determining the relative orientation of CS and EFG tensors of O-6 and nitrogen atoms in the molecular frame axes.     

Influence of N-H...O and O-H...O hydrogen bonds on the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen: a density functional theory study

A computational investigation was carried out to characterize the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen. We found that N-H...O and O-H...O hydrogen bonds around the acetaminophen molecule in the crystal lattice have different influences on the calculated (17)O, (15)N and (13)C chemical shielding eigenvalues and their orientations in the molecular frame of axes. The calculations were performed with the B3LYP method and 6-311++G(d, p) and 6-311+G(d) standard basis sets using the Gaussian 98 suite of programs. Calculated chemical shielding tensors were used to evaluate the (17)O, (15)N, and (13)C NMR chemical shift tensors in crystalline acetaminophen, which are in reasonable agreement with available experimental data. The difference between the calculated NMR parameters of the monomer and molecular clusters shows how much hydrogen-bonding interactions affect the chemical shielding tensors of each nucleus. The computed (17)O chemical shielding tensor on O(1), which is involved in two intermolecular hydrogen bonds, shows remarkable sensitivity toward the choice of the cluster model, whereas the (17)O chemical shielding tensor on O(2) involved in one N-H...O hydrogen bond, shows smaller improvement toward the hydrogen-bonding interactions. Also, a reasonably good agreement between the experimentally obtained solid-state (15)N and (13)C NMR chemical shifts and B3LYP/6-311++G(d, p) calculations is achievable only in molecular cluster model where a complete hydrogen-bonding network is considered. Moreover, at the B3LYP/6-311++G(d, p) level of theory, the calculated (17)O, (15)N and (13)C chemical shielding tensor orientations are able to reproduce the experimental values to a reasonably good degree of accuracy.     

A theoretical study of repeating sequence in HRP II: a combination of molecular dynamics simulations and (17)O quadrupole coupling tensors

Histidine rich protein II derived peptide (HRP II 169-182) was investigated by molecular dynamics, MD, simulation and (17)O electric field gradient, EFG, tensor calculations. MD simulation was performed in water at 300 K with alpha-helix initial structure. It was found that peptide loses its initial alpha-helix structure rapidly and is converted to random coil and bent secondary structures. To understand how peptide structure affects EFG tensors, initial structure and final conformations resulting from MD simulations were used to calculate (17)O EFG tensors of backbone carbonyl oxygens. Calculations were performed using B3LYP method and 6-31+G basis set. Calculated (17)O EFG tensors were used to evaluate quadrupole coupling constants, QCC, and asymmetry parameters, eta(Q). Difference between the calculated QCC and eta(Q) values revealed how hydrogen-bonding interactions affect EFG tensors at the sites of each oxygen nucleus.     

<Emphasis Type="Italic">Ab initio</Emphasis> DFT study of bisphosphonate derivatives as a drug for inhibition of cancer: NMR and NQR parameters

DFT computations were carried out to characterize the ¹⁷Oand ²H electric field gradient, EFG, in various bisphosphonate derivatives. The computations were performed at the B3LYP level with 6-311++G (d,P) standard basis set. Calculated EFG tensors were used to determine the ¹⁷O and ²H nuclear quadrupole coupling constant, χ and asymmetry parameter, η. For better understanding of the bonding and electronic structure of bisphosphonates, isotropic and anisotropic NMR chemical shieldings were calculated for the ¹³C, ¹⁷O and ³¹P nuclei using GIAO method for the optimized structure of intermediate bisphosphonates at B3LYP level of theory using 6-311++G (d, p) basis set. The results showed that various substituents have a strong effect on the nuclear quadrupole resonance (NQR) parameters (χ, η) of ¹⁷O in contrast with ²H NQR parameters. The NMR and NQR parameters were studied in order to find the correlation between electronic structure and the activity of the desired bisphosphonates. In addition, the effect of substitutions on the bisphosphonates polarity was investigated. Molecular polarity was determined via the DFT calculated dipole moment vectors and the results showed that substitution of bromine atom on the ring would increase the activity of bisphosphonates.     

Role of spin state on the geometry and nuclear quadrupole resonance parameters in hemin complex

Theoretical calculations of structural parameters, 57Fe, 14N and 17 O electric field gradient (EFG) tensors for full size-hemin group have been carried out using density functional theory. These calculations are intended to shed light on the difference between the geometry parameters, nuclear quadrupole coupling constants (QCC), and asymmetry parameters (eta Q) found in three spin states of hemin; doublet, quartet and sextet. The optimization results reveal a significant change for propionic groups and porphyrin plane in different spin states. It is found that all principal components of EFG tensor at the iron site are sensitive to electronic and geometry structures. A relationship between the EFG tensor at the 14N and 17 O sites and the spin state of hemin complex is also detected.     

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.     

A theoretical study of cooperative and anticooperative effects on hydrogen-bonded clusters of water and the cyanuric acid

Ab initio and density functional calculations are used to analyse the interaction between a molecule of the cyanuric acid and one, two and three molecules of water at B3LYP/6-311++ G(d,p) and MP2/6-311++ G(d,p) computational levels. Also, the cooperative effect (CE) in terms of the stabilisation energy of clusters is calculated and discussed. Depending on the geometry of clusters under study, the cooperative, non- or anti-CE was found with an increasing cluster size. Red shifts of N–H and C = O stretching frequencies illustrate a good dependence on the CE. The atoms in molecules theory is used to analyse the CE on topological parameters.     

Correlation between NQR parameters and residue size of aliphatic amino acids and their dimers

Nuclear quadrupole coupling constants (NQCC), chi, and asymmetry parameters, eta, of 2D, 14N and 17O nuclei have been calculated for aliphatic amino acids and their dimers using MP2/6-311++G** method to shed some light on the differences between the structural parameters in the aliphatic amino acids and their dimers. For this purpose, electric field gradient (EFG) at the sites of quadrupolar nuclei have been calculated and evaluated for each compound. A correlation is observed between the calculated NQCC parameters and the conformational structures of the compounds, showing that extraction of structural data from the NQR spectra might be promising. Our results showed that 17O NQCC of terminal carboxylic acid and 14N NQCC of the terminal amino groups are, respectively, the least and the most sensitive parameters to the variation of the size of the residue. It is found also that conformation of R (i.e. values of the dihedral angles) plays a very effective role in the determination of the values of the calculated NQCC parameters. Sensitivity of the NQR parameters to the changes in the conformational structure is significantly greater (nearly 20-fold) than that to the changes in the other structural parameters such as bond lengths.     

NMR studies of exocyclic 1,N2-propanodeoxyguanosine adducts (X) opposite purines in DNA duplexes: protonated X(syn).A(anti) pairing (acidic pH) and X(syn).G(anti) pairing (neutral pH) at the lesion site

Proton and phosphorus two-dimensional NMR studies are reported for the complementary d(C1-A2-T3-G4-X5-G6-T7-A8-C9).d(G10-T11-A12-C13-A14-C15-A 16-T17-G18) nonanucleotide duplex (designated X.A 9-mer) that contains a 1,N2-propanodeoxyguanosine exocyclic adduct, X5, opposite deoxyadenosine A14 in the center of the helix. The NMR studies detect a pH-dependent conformational transition; this paper focuses on the structure present at pH 5.8. The two-dimensional NOESY studies of the X.A 9-mer duplex in H2O and D2O solution establish that X5 adopts a syn orientation while A14 adopts an anti orientation about the glycosidic bond at the lesion site. The large downfield shift of the amino protons of A14 demonstrates protonation of the deoxyadenosine base at pH 5.8 such that the protonated X5(syn).A14(anti) pair is stabilized by two hydrogen bonds at low pH. At pH 5.8, the observed NOE between the H8 proton of X5 and the H2 proton of A14 in the X.A 9-mer duplex demonstrates unequivocally the formation of the protonated X5(syn).A14(anti) pair. The 1,N2-propano bridge of X5(syn) is located in the major groove. Selective NOEs from the exocyclic methylene protons of X5 to the major groove H8 proton of flanking G4 but not G6 of the G4-X5-G6 segment provide additional structural constraints on the local conformation at the lesion site. A perturbation in the phosphodiester backbone is detected at the C13-A14 phosphorus located at the lesion site by 31P NMR spectroscopy. The two-dimensional NMR studies have been extended to the related complementary X.G 9-mer duplex that contains a central X5.G14 lesion in a sequence that is otherwise identical with the X.A 9-mer duplex. The NMR experimental parameters are consistent with formation of a pH-independent X5(syn).G14(anti) pair stabilized by two hydrogen bonds with the 1,N2-propano exocyclic adduct of X5(syn) located in the major groove.     

Theoretical studies on the structures,intra- and inter-molecular hydrogen bonding interactions in HNF and HNF–H2O clusters in the gaseous,aqueous and solid phases

Hydrazimium nitroformate ([N₂H₅]⁺[C(NO₂)₃]^? , HNF) is an ionic oxidiser used in solid propellants. Its properties are easily affected by H₂O because of its hygroscopicity. In this article, density functional theory (DFT) and molecular dynamics (MD) were employed to study the isolated HNF molecule and the HNF–H₂O cluster in gas phase and in the aqueous solution. Three stable conformations were obtained for HNF in the gas phase and in the aqueous solution, respectively, and each conformation can form several different HNF–H₂O clusters. Irrespective of whether it is in gas phase or in solution, intramolecular hydrogen bond interactions and other interactions (e.g. the binding energy, the dispersion energy, the second-order perturbation energy and the energy gap between frontier orbitals) of HNF are weaker in the clusters than in the isolated state. The initial decomposition energy of the cluster is lower than that of the isolated HNF molecule in both gaseous and aqueous phases, while the dissociation processes are the same. Molecular dynamic simulations showed that the clustered H₂O elongates and weakens the C–NO₂ bond in the solid HNF–H₂O cluster compared with that in the solid HNF. H₂O reduces and weakens intramolecular N–HΛO bonds too, and O–HΛN is the dominant intermolecular hydrogen bond between HNF and H₂O.     

Complete conformational family of 2',3'-didehydro-2',3'-dideoxyguanosine: quantum chemical and electron density topological study

Comprehensive conformational analysis of the biologically active nucleoside 2',3'-didehydro-2',3'-dideoxyaguanosine (d4G) has been performed at the MP2/6-311++G(d,p)//DFT B3LYP/6-31G(d,p) level of theory. The energetic, geometrical and polar characteristics of twenty d4G conformers as well as their conformational equilibrium were investigated. The electron density topological analysis allowed us to establish that the d4G molecule is stabilized by nine types of intramolecular interactions: O5'H...N3, O5'H...C8, C8H...O5', C2'H...N3, C5'H1...N3, C5'H2...N3, C8H...H1C5', C8H...H2'C5' and N2H1...O5'. The obtained results of conformational analysis permit us to think that d4G may be a terminator of the DNA chain synthesis in the 5'-3' direction. Thus it can be inferred that d4G competes with canonical 2'-deoxyaguanosine in binding an active site of the corresponding enzyme.     

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.

     

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.     

Dissection of the structural and functional role of a conserved hydration site in RNase T1

The reoccurrence of water molecules in crystal structures of RNase T1 was investigated. Five waters were found to be invariant in RNase T1 as well as in six other related fungal RNases. The structural, dynamical, and functional characteristics of one of these conserved hydration sites (WAT1) were analyzed by protein engineering, X-ray crystallography, and (17)O and 2H nuclear magnetic relaxation dispersion (NMRD). The position of WAT1 and its surrounding hydrogen bond network are unaffected by deletions of two neighboring side chains. In the mutant Thr93Gln, the Gln93N epsilon2 nitrogen replaces WAT1 and participates in a similar hydrogen bond network involving Cys6, Asn9, Asp76, and Thr91. The ability of WAT1 to form four hydrogen bonds may explain why evolution has preserved a water molecule, rather than a side-chain atom, at the center of this intricate hydrogen bond network. Comparison of the (17)O NMRD profiles from wild-type and Thr93Gln RNase T1 yield a mean residence time of 7 ns at 27 degrees C and an orientational order parameter of 0.45. The effects of mutations around WAT1 on the kinetic parameters of RNase T1 are small but significant and probably relate to the dynamics of the active site.     

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.

Protonation and deprotonation enthalpies of alloxan and implications for the structure and energy of its complexes with water: a computational study

The optimized geometries, harmonic vibrational frequencies, and energies of the structures of monohydrated alloxan were computed at the DFT/ωB97X-D and B3LYP/6–311++G** level of theory. Results confirm that the monohydrate exists as a dipolar alloxan–water complex which represents a global minimum on the potential energy surface (PES). Trajectory dynamics simulations show that attempt to reorient this monohydrate, to a more favorable orientation for H-bonding, is opposed by an energy barrier of 25.07?kJ/mol. Alloxan seems to prefer acting as proton donor than proton acceptor. A marked stabilization due to the formation of N–H–OH₂ bond is observed. The concerted proton donor–acceptor interaction of alloxan with one H₂O molecule does not increase the stability of the alloxan–water complex. The proton affinity of the O and N atoms and the deprotonation enthalpy of the NH bond of alloxan are computed at the same level of theory. Results are compared with recent data on uracil, thymine, and cytosine. The intrinsic acidities and basicities of the four pyrimidines were discussed. Results of the present study reveal that alloxan is capable of forming stronger H-bonds and more stable cyclic complex with water; yet it is of much lower basicity than other pyrimidines.     

Nonorthogonal tight-binding model with H–C–N–O parameterisation

A parametric nonorthogonal tight-binding model (NTBM1) with the set of parameters for H–C–N–O systems is presented. This model compares well with widely used semi-empirical AM1 and PM3/PM7 models but contains less fitting parameters per atom. All NTBM1 parameters are derived based on a criterion of the best agreement between the calculated and experimental values of bond lengths, valence angles and binding energies for various H–C–N–O molecules. Results for more than 200 chemical compounds are reported. Parameters are currently available for hydrogen, carbon, nitrogen, oxygen atoms and corresponding interatomic interactions. The model has a good transferability and can be used for both relaxation of large molecular systems (e.g., high-molecular compounds or covalent cluster complexes) and long-timescale molecular dynamics simulation (e.g., modelling of thermal decomposition processes). The program package based on this model is available for download at no cost from http://ntbm.info.     

Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.