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.     

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.     

A density functional study of (17)O, (14)N and (2)H electric field gradient tensors in the real crystalline structure of alpha-glycine

A density functional theory (DFT) study was carried out to calculate (17)O, (14)N and (2)H electric field gradient (EFG) tensors in accurate neutron diffraction structures of alpha-glycine at 288 and 427 K. B3LYP is the used method and 6-311+G(*) and 6-311++G(**) are the basis sets in the calculations of EFG tensors at the sites of (17)O, (14)N and (2)H nuclei in the monomer and the octameric cluster of alpha-glycine at two temperatures. Quadrupole coupling constants and asymmetry parameters are the converted parameters of calculated EFG tensors to experimentally measurable ones. The calculated results of monomer and the target molecule in octameric cluster reveal that hydrogen-bonding interactions play an important role in the crystalline structure of alpha-glycine where the results of the target molecule in octameric cluster are in good agreement with the experiments.     

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.     

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.     

<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.     

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.     

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.     

Insight into the copper coordination environment in the prion protein through density functional theory calculations of EPR parameters

Density functional theory (DFT) calculations of Cu(II) electron paramagnetic resonance (EPR) parameters for the octarepeat unit of the prion protein were conducted. Model complexes were constructed and optimized using the crystal structure of the octarepeat unit of the prion protein. Copper g and A tensors and nitrogen hyperfine and quadrupole coupling constants were calculated using DFT. Solvent effects were incorporated using the conductor-like screening model as well as through the inclusion of explicit water molecules. Calculations using the model with an additional axial water molecule added to the coordination sphere of the Cu(II) metal center give the best qualitative agreement for the copper g and A tensors. The S-band experimental EPR spectra were interpreted in light of the DFT calculations of the directly coordinated nitrogen hyperfine coupling constants which indicate that the three directly coordinated nitrogen atoms in the octarepeat unit are not equivalent. These results demonstrate that DFT calculations of EPR parameters can provide important insight with respect to the structural interpretation of experimental EPR data. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

N-H...O, O-H...O, and C-H...O hydrogen bonds in protein-ligand complexes: strong and weak interactions in molecular recognition

The characteristics of N-H...O, O-H...O, and C-H...O hydrogen bonds are examined in a group of 28 high-resolution crystal structures of protein-ligand complexes from the Protein Data Bank and compared with interactions found in small-molecule crystal structures from the Cambridge Structural Database. It is found that both strong and weak hydrogen bonds are involved in ligand binding. Because of the prevalence of multifurcation, the restrictive geometrical criteria set up for hydrogen bonds in small-molecule crystal structures may need to be relaxed in macromolecular structures. For example, there are definite deviations from linearity for the hydrogen bonds in protein-ligand complexes. The formation of C-H...O hydrogen bonds is influenced by the activation of the C(alpha)-H atoms and by the flexibility of the side-chain atoms. In contrast to small-molecule structures, anticooperative geometries are common in the macromolecular structures studied here, and there is a gradual lengthening as the extent of furcation increases. C-H...O bonds formed by Gly, Phe, and Tyr residues are noteworthy. The numbers of hydrogen bond donors and acceptors agree with Lipinski's "rule of five" that predicts drug-like properties. Hydrogen bonds formed by water are also seen to be relevant in ligand binding. Ligand C-H...O(w) interactions are abundant when compared to N-H...O(w) and O-H...O(w). This suggests that ligands prefer to use their stronger hydrogen bond capabilities for use with the protein residues, leaving the weaker interactions to bind with water. In summary, the interplay between strong and weak interactions in ligand binding possibly leads to a satisfactory enthalpy-entropy balance. The implications of these results to crystallographic refinement and molecular dynamics software are discussed.     

Perturbation of hydrogen bonds in the adenine ... thymine base pair by Na+, Mg2+, Ca2+ and NH4+ cations

Perturbation of the hydrogen bonds in the adenine ... thymine base pair by Na+, Mg2+, Ca2+ and NH4+ cations has been investigated by means of ab initio SCF calculations with the STO-3G basis set. The geometry of adenine...thymine, as well as those of the perturbed pairs were optimized. Approach of any cation to thymine at O6 leads to destabilization of the adenine...thymine pair; divalent cations (Mg2+, Ca2+) have a profound effect on the structure of the base pair. The approach of a cation to other available sites (thymine: O2, adenine N1 and N3) leads, on the other hand, to stabilization of the base pair. If a water molecule is placed between the cation and the base pair, the structure and stability of the base pair are changed only negligibly.     

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.     

Solution and solid-state structure of N-acetamido-3,4,6-tri-O-acetyl-2-azido-2-deoxy-alpha-D-galactopyranosylamine

High resolution 1H NMR and 13C NMR spectroscopic and single crystal X-ray structural analyses of N-acetamido-3,4,6-tri-O-acetyl-2-azido-2-deoxy-alpha-D-galactopyranosylamine (1), a minor product of azidonitration reaction of 3,4,6-tri-O-acetyl galactal, are reported. The solution phase studies of 1 reflect that the compound exists in 4C1 conformation with cis-orientation of the substituents at C-1 and C-2. The solid-state structure of 1 reveals that a molecule of water is entrapped in the solid state of 1 and this water molecule serves to mediate N-H...O and C-H...O interactions.     

Simulation of interactions between nucleic acid bases by refined atom-atom potential functions

Energy of interaction between nitrogen bases of nucleic acid has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H . . . O and (or) N-H . . . N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H . . . O or N-H . . . N and one C-H . . . O or C-H . . . N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum-mechanical approach.     

Simulation of Interactions between Nucleic Acid Bases by Refined Atom-Atom Potential Functions

Abstract

Energy of interaction between nitrogen bases of nucleic acids has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H … O and (or) N-H … N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component, in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H … O or N-H … N and one C-H … O or C-H … N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum- mechanical approach.     

C-H ... O hydrogen bonding in solutions of methylated nucleic acid base analogs as revealed by NMR

Formation and thermodynamic characteristics of C-H ... O hydrogen bonding of methylated uracils and caffeine have been studied by nmr along two lines. 1. The concentration and temperature dependencies of the PMR spectra of 1,3-dimethyluracil (m2 1,3Ura), 1,3-dimethylthymine (m2 1,3Thy), and 1,3,6-trimethyluracil (m3 1,3,6Ura) in chloroform at high concentrations of base analogs indicated the self-association of m2 1,3Ura and m2 1,3Thy via C(6)H ... O hydrogen bonding and the competitive formation of C-H ... O bonds between carbonyl oxygens and chloroform. The intermolecular interaction energy and the arrangement of molecules in the local minima of various m2 1,3Ura dimers were calculated by the method of atom-atom potentials. The deepest minimum for the m2 1,3Ura coplanar dimer corresponds to a C(6)-H ... O hydrogen-bond formation. 2. At low concentration of m2 1,3Ura and caffeine in CCl4, C(6)-H ... O bonding for m2 1,3Ura and C(8)-H ... O bonding for caffeine with oxygens of dimethyl sulfoxide (DMSO) and acetone were observed. The association constants of these complexes were obtained at different temperatures. The enthalpies delta H, of the m2 1,3Ura-DMSO, m2 1,3Ura-accetone, caffeine-DMSO, and caffeine-acetone complexes were -2 +/- 0.1 kcal/mol. The calculations showed that the deepest minimum of the caffeine-acetone coplanar complex corresponds to C(8)-H ... O bonding with energy of -3.5 kcal/mol and that of the m2 1,3Ura-acetone complexes corresponds to C(6)-H ... O bonding with energy of -3.4 kcal/mol. The approximate correction for the solvent effect provides good agreement of the experimental data with the calculations.     

C9 conformation of N-(N alpha-[(tert.-butyloxy)-carbonyl]-L-alanyl)-N,N'-dicyclohexylu rea in solid and solution.

An X-ray diffraction study was carried out on a single crystal of N-(N alpha-[(tert.-butyloxy)-carbonyl]-L-alanyl)-N,N'-dicyclohexylur ea belonging to the tetragonal space group P4(1)2(1)2, having cell dimensions a = b = 10.102(3) A, c = 46.067(7) A, V = 4701.2 A3, Z = 8. The crystal structure was solved by direct methods and refined to an R value of 0.056 for 1602 unique reflections with I greater than 2.5 sigma(I). Crystal structure analysis shows the presence of an intramolecular N-H ... O=C H-bond stabilizing the molecule in a folded form similar to that of a beta turn, forming a nine-membered ring. IR and 1H-NMR studies in CDCl3 solution confirm the stable folded conformation found in the crystalline state, as well as the existence of N-H ... O=C H-bonds in the title compound, as in peptides.     

Conformational properties of N-acetyl-L-alanine N',N'-dimethylamide

Ab initio/DFT analysis of the conformational properties of free Ac-Ala-NMe(2) (N-acetyl-L-alanine-N',N'-dimethylamide) in terms of the N-H.O, N-H.N, C-H.O hydrogen bonds and C(delta+) = O(delta-) dipole attractions was performed. The Ala residue combined with the C-terminal tertiary amide prefers an extended conformation and that characteristic of the (i + 1)th position of the betaVIb turn. These can be easily remodelled into a structure compatible with the (i + 1)th position of the betaII/betaVIa turn. The residue has also the potential to adopt the conformation accommodated at both central positions of the betaIII/betaIII' turn or the (i + 1)th position of the betaI/beta'I turn.     

Observation of a unique pattern of bifurcated hydrogen bonds in the crystal structures of the N-glycoprotein linkage region models

Elucidation of the intra- and intermolecular carbohydrate-protein interactions would greatly contribute toward obtaining a better understanding of the structure-function correlations of the protein-linked glycans. The weak interactions involving C-H...O have recently been attracting immense attention in the domain of biomolecular recognition. However, there has been no report so far on the occurrence of C-H...O hydrogen bonds in the crystal structures of models and analogs of N-glycoproteins. We present herein an analysis of C-H...O interactions in the crystal structures of all N-glycoprotein linkage region models and analogs. The study reveals a cooperative network of bifurcated hydrogen bonds consisting of N-H...O and C-H...O interactions seen uniquely for the models. The cooperative network consists of two antiparallel chains of bifurcated hydrogen bonds, one involving N1-H, C2'-H and O1' of the aglycon moiety and the other involving N2-H, C1-H and O1' of the sugar. Such bifurcated hydrogen bonds between the core glycan and protein are likely to play an important role in the folding and stabilization of proteins.