Intra- and intermolecular interactions in small bioactive molecules: cooperative features from experimental and theoretical charge-density analysis

2006 ◽  
Vol 62 (4) ◽  
pp. 612-626 ◽  
Author(s):  
Parthapratim Munshi ◽  
Tayur N. Guru Row

The topological features of the charge densities, ρ(r), of three bioactive molecules, 2-thiouracil [2,3-dihydro-2-thioxopyrimidin-4(1H)-one], cytosine [4-aminopyrimidin-2(1H)-one] monohydrate and salicylic acid (2-hydroxybenzoic acid), have been determined from high-resolution X-ray diffraction data at 90 K. The corresponding results are compared with the periodic theoretical calculations, based on theoretical structure factors, performed using DFT (density-functional theory) at the B3LYP/6-31G** level. The molecules pack in the crystal lattices via weak intermolecular interactions as well as strong hydrogen bonds. All the chemical bonds, including the intra- and intermolecular interactions in all three compounds, have been quantitatively described by topological analysis based on Bader's quantum theory of `Atoms In Molecules'. The roles of interactions such as C—H...O, C—H...S, C—H...π and π...π have been investigated quantitatively in the presence of strong hydrogen bonds such as O—H...O, N—H...O and N—H...S, based on the criteria proposed by Koch and Popelier to characterize hydrogen bonds and van der Waals interactions. The features of weak intermolecular interactions, such as S...S in 2-thiouracil, the hydrogen bonds generated from the water molecule in cytosine monohydrate and the formation of the dimer via strong hydrogen bonds in salicylic acid, are highlighted on a quantum basis. Three-dimensional electrostatic potentials over the molecular surfaces emphasize the preferable binding sites in the structure and the interaction features of the atoms in the molecules, which are crucial for drug–receptor recognition.

2006 ◽  
Vol 84 (5) ◽  
pp. 804-811 ◽  
Author(s):  
David Wolstenholme ◽  
Manuel AS Aquino ◽  
T Stanley Cameron ◽  
Joseph D Ferrara ◽  
Katherine N Robertson

The tetraphenylphosphonium squarate salt crystallizes with a number of diverse interactions, which all have the potential to be classified as hydrogen bonds. The squarate anions are found as dimers linked by O-H···O interactions. The multipole refinement of the tetraphenylphosphonium squarate was performed using the Hansen–Coppens model followed by topological analysis of its intermolecular interactions. A total of 28 interactions were found among the symmetry related molecules, which include a number of C-H···Cπ, C-H···O, and C-H···H-C interactions, along with the O-H···O interaction. With the criteria for hydrogen bonding proposed by Popelier and Koch, it is possible to determine which of these interactions are hydrogen bonds and which are van der Waals interactions. Both linear and exponentially dependent correlations can be seen for the properties of the bond critical points involving the intermolecular interactions that fulfill these criteria. All this leads to a better understanding of the role that hydrogen bonds play in the formation of small organic compounds.Key words: electron density, multiple refinement, hydrogen bonds.


2009 ◽  
Vol 65 (6) ◽  
pp. 757-769 ◽  
Author(s):  
Louis J. Farrugia ◽  
Pavel Kočovský ◽  
Hans Martin Senn ◽  
Štěpán Vyskočil

The charge density in (±)-8′-benzhydrylideneamino-1,1′-binaphthyl-2-ol (1) has been studied experimentally using Mo Kα X-ray diffraction at 100 K, and by theory using density-functional thoery (DFT) calculations at the B3LYP/6-311++G** level. The nature of the weak intramolecular peri-C...N, CH...π, H...H and C(π)...C(π) interactions has been examined by topological analysis using the Quantum Theory of Atoms in Molecules (QTAIM) approach. An analysis of the density ρ(r), the Laplacian of the density ∇2ρ(r b) and other topological properties at the bond-critical points were used to classify these interactions. The study confirms the presence of the intramolecular CH...π interaction in (1), which was previously suspected on geometrical grounds. An analysis of the ellipticity profiles along the bond paths unambiguously shows the π-delocalization between the imine unit and one N-phenyl group. The weak intermolecular interactions in the crystal of (1) were examined experimentally and theoretically through the pairwise interactions of the seven independent dimeric pairs of (1) responsible for the set of unique intermolecular interactions, and also through examination of the Hirshfeld surface d norm property. The theoretical dimeric-pair calculations used the BLYP-D functional which supplements the exchange-correlational functional with an empirical dispersion term to provide a more accurate determination of the energies for the weak intermolecular interactions.


2021 ◽  
Vol 36 (1) ◽  
pp. 35-42
Author(s):  
Shivang Bhaskar ◽  
Joseph T. Golab ◽  
James A. Kaduk ◽  
Amy M. Gindhart ◽  
Thomas N. Blanton

The crystal structure of pimecrolimus Form B has been solved and refined using synchrotron X-ray powder diffraction data and optimized using density functional techniques. Pimecrolimus crystallizes in the space group P21 (#4) with a = 15.28864(7), b = 13.31111(4), c = 10.95529(5) Å, β = 96.1542(3)°, V = 2216.649(9) Å3, and Z = 2. Although there are an intramolecular six-ring hydrogen bond and some larger chain and ring patterns, the crystal structure is dominated by van der Waals interactions. There is a significant difference between the conformation of the Rietveld-refined and the DFT-optimized structures in one portion of the macrocyclic ring. Although weak, intermolecular interactions are apparently important in determining the solid-state conformation. The powder pattern is included in the Powder Diffraction File™ (PDF®) as entry 00-066-1619. This study provides the atomic coordinates to be added to the PDF entry.


Author(s):  
Sehrish Akram ◽  
Arshad Mehmood ◽  
Sajida Noureen ◽  
Maqsood Ahmed

Thermal-induced transformation of glutamic acid to pyroglutamic acid is well known. However, confusion remains over the exact temperature at which this happens. Moreover, no diffraction data are available to support the transition. In this article, we make a systematic investigation involving thermal analysis, hot-stage microscopy and single-crystal X-ray diffraction to study a one-pot thermal transition of glutamic acid to pyroglutamic acid and subsequent self-cocrystallization between the product (hydrated pyroglutamic acid) and the unreacted precursor (glutamic acid). The melt upon cooling gave a robust cocrystal, namely, glutamic acid–pyroglutamic acid–water (1/1/1), C5H7NO3·C5H9NO4·H2O, whose structure has been elucidated from single-crystal X-ray diffraction data collected at room temperature. A three-dimensional network of strong hydrogen bonds has been found. A Hirshfeld surface analysis was carried out to make a quantitative estimation of the intermolecular interactions. In order to gain insight into the strength and stability of the cocrystal, the transferability principle was utilized to make a topological analysis and to study the electron-density-derived properties. The transferred model has been found to be superior to the classical independent atom model (IAM). The experimental results have been compared with results from a multipolar refinement carried out using theoretical structure factors generated from density functional theory (DFT) calculations. Very strong classical hydrogen bonds drive the cocrystallization and lend stability to the resulting cocrystal. Important conclusions have been drawn about this transition.


2021 ◽  
Vol 77 (2) ◽  
pp. 116-122
Author(s):  
Maria Marczak ◽  
Kinga Biereg ◽  
Beata Zadykowicz ◽  
Artur Sikorski

The synthesis and structural characterization of the monohydrated 1:2 cocrystal salt of acriflavine with 3,5-dinitrobenzoic acid [systematic name: 3,6-diamino-10-methylacridin-10-ium 3,5-dinitrobenzoate–3,5-dinitrobenzoic acid–water (1/1/1), C14H14N3 +·C7H3N2O6 −·C7H4N2O6·H2O] are reported. Single-crystal X-ray diffraction measurements show that the title solvated monohydrate salt crystalizes in the monoclinic space group P21 with one acriflavine cation, a 3,5-dinitrobenzoate anion, a 3,5-dinitrobenzoic acid molecule and a water molecule in the asymmetric unit. The neutral and anionic forms of 3,5-dinitrobenzoic acid are linked via O—H...O hydrogen bonds to form a monoanionic dimer. Neighbouring monoanionic dimers of 3,5-dinitrobenzoic acid are linked by nitro–nitro N—O...N and nitro–acid N—O...π intermolecular interactions to produce a porous organic framework. The acriflavine cations are linked with carboxylic acid molecules directly via amine–carboxy N—H...O, amine–nitro N—H...O and acriflavine–carboxy C—H...O hydrogen bonds, and carboxy–acriflavine C—O...π, nitro–acriflavine N—O...π and acriflavine–nitro π–π interactions, or through the water molecule by amino–water N—H...O and water–carboxy O—H...O hydrogen bonds, and are located in the voids of the porous organic framework. The intermolecular interactions were studied using the CrystalExplorer program to provide information about the interaction energies and the dispersion, electrostatic, polarization and repulsion contributions to the lattice energy.


2016 ◽  
Vol 18 (42) ◽  
pp. 29249-29257 ◽  
Author(s):  
Chengqian Yuan ◽  
Haiming Wu ◽  
Meiye Jia ◽  
Peifeng Su ◽  
Zhixun Luo ◽  
...  

Utilizing dispersion-corrected density functional theory (DFT) calculations, we demonstrate the weak intermolecular interactions of phenylenediamine dimer (pdd) clusters, emphasizing the local lowest energy structures and decomposition of interaction energies by natural bond orbital (NBO) and atoms in molecule (AIM) analyses.


Author(s):  
Piyush Panini ◽  
K. N. Venugopala ◽  
Bharti Odhav ◽  
Deepak Chopra

A new polymorph belonging to the tetrahydropyrimidinium class of compounds, namely 6-(4-chlorophenyl)-5-(methoxycarbonyl)-4-methyl-2-(3-(trifluoromethylthio)phenylamino)-3,6-dihydropyrimidin-1-ium chloride, and a hydrate of 2-(3-bromophenylamino)-6-(4-chlorophenyl)-5-(methoxycarbonyl)-4-methyl-3,6-dihydropyrimidin-1-ium chloride, have been isolated and characterized using single-crystal X-ray diffraction (XRD). A detailed comprehensive analysis of the crystal packing in terms of the associated intermolecular interactions and a quantification of their interaction energies have been performed for both forms of the two different organic salts (AandB) using X-ray crystallography and computational methods such as density functional theory (DFT) quantum mechanical calculations, PIXEL lattice-energy calculations (with decomposition of total lattice energy into the Coulombic, polarization, dispersion and repulsion contribution), the calculation of the Madelung constant (the EUGEN method), Hirshfeld and two-dimensional fingerprint plots. The presence of ionic [N—H]+...Cl−and [C—H]+...Cl−hydrogen bonds mainly stabilizes the crystal packing in both formsAandB, while in the case ofB·H2O [N—H]+...Owaterand Owater—H...Cl−hydrogen bonds along with [N—H]+...Cl−and [C—H]+...Cl−provide stability to the crystal packing. The lattice-energy calculations from both PIXEL and EUGEN methods revealed that in the case ofA, form (I) (monoclinic) is more stable whereas forBit is the anhydrous form that is more stable. The analysis of the `Madelung mode' of crystal packing of two forms ofAandBand its hydrates suggest that differences exist in the position of the charged ions/atoms in the organic solid state. TheR/E(distance–energy) plots for all the crystal structures show that the molecular pairs in their crystal packing are connected with either highly stabilizing (due to the presence of organicR+and Cl−) or highly destabilizing Coulombic contacts. The difference in crystal packing and associated intermolecular interactions between polymorphs (in the case ofA) or the hydrates (in the case ofB) have been clearly elucidated by the analysis of Hirshfeld surfaces and two-dimensional fingerprint plots. The relative contributions of the various interactions to the Hirshfeld surface for the cationic (dihydropyrimidinium) part and anionic (chloride ion) part for the two forms ofAandBand its hydrate were observed to be different.


Author(s):  
Chinnasamy Kalaiarasi ◽  
Mysore S Pavan ◽  
Poomani Kumaradhas

An experimental charge density distribution of 2-nitroimidazole was determined from high-resolution X-ray diffraction and the Hansen–Coppens multipole model. The 2-nitroimidazole compound was crystallized and a high-angle X-ray diffraction intensity data set has been collected at low temperature (110 K). The structure was solved and further, an aspherical multipole model refinement was performed up to octapole level; the results were used to determine the structure, bond topological and electrostatic properties of the molecule. In the crystal, the molecule exhibits a planar structure and forms weak and strong intermolecular hydrogen-bonding interactions with the neighbouring molecules. The Hirshfeld surface of the molecule was plotted, which explores different types of intermolecular interactions and their strength. The topological analysis of electron density at the bond critical points (b.c.p.) of the molecule was performed, from that the electron density ρbcp(r) and the Laplacian of electron density ∇2ρbcp(r) at the b.c.p.s of the molecule have been determined; these parameters show the charge concentration/depletion of the nitroimidazole bonds in the crystal. The electrostatic parameters like atomic charges and the dipole moment of the molecule were calculated. The electrostatic potential surface of the molecule has been plotted, and it displays a large electronegative region around the nitro group. All the experimental results were compared with the corresponding theoretical calculations performed usingCRYSTAL09.


2009 ◽  
Vol 65 (5) ◽  
pp. 639-646 ◽  
Author(s):  
Edward E. Ávila ◽  
Asiloé J. Mora ◽  
Gerzon E. Delgado ◽  
Ricardo R. Contreras ◽  
Luis Rincón ◽  
...  

The molecular and crystalline structure of ethyl 1′,2′,3′,4′,4a′,5′,6′,7′-octahydrodispiro[cyclohexane-1,2′-quinazoline-4′,1′′-cyclohexane]-8′-carbodithioate (I) was solved and refined from powder synchrotron X-ray diffraction data. The initial model for the structural solution in direct space using the simulated annealing algorithm implemented in DASH [David et al. (2006). J. Appl. Cryst. 39, 910–915] was obtained performing a conformational study on the fused six-membered rings of the octahydroquinazoline system and the two spiran cyclohexane rings of (I). The best model was chosen using experimental evidence from 1H and 13C NMR [Contreras et al. (2001). J. Heterocycl. Chem. 38, 1223–1225] in combination with semi-empirical AM1 calculations. In the refined structure the two spiran rings have the chair conformation, while both of the fused rings in the octahydroquinazoline system have half-chair conformations compared with in-vacuum density-functional theory (DFT) B3LYP/6-311G*, DFTB (density-functional tight-binding) theoretical calculations in the solid state and other related structures from X-ray diffraction data. Compound (I) presents weak intramolecular hydrogen bonds of the type N—H...S and C—H...S, which produce delocalization of the electron density in the generated rings described by graph symbols S(6) and S(5). Packing of the molecules is dominated by van der Waals interactions.


Sign in / Sign up

Export Citation Format

Share Document