Variable Temperature Neutron Diffraction Analysis of a Very Short O−H···O Hydrogen Bond in 2,3,5,6-Pyrazinetetracarboxylic Acid Dihydrate: Synthon-Assisted Short Oacid−H···Owater Hydrogen Bonds in a Multicenter Array

Combined X-ray Crystallographic, IR/Raman Spectroscopic, and Periodic DFT Investigations of New Multicomponent Crystalline Forms of Anthelmintic Drugs: A Case Study of Carbendazim Maleate

Synthesis of multicomponent solid forms is an important method of modifying and fine-tuning the most critical physicochemical properties of drug compounds. The design of new multicomponent pharmaceutical materials requires reliable information about the supramolecular arrangement of molecules and detailed description of the intermolecular interactions in the crystal structure. It implies the use of a combination of different experimental and theoretical investigation methods. Organic salts present new challenges for those who develop theoretical approaches describing the structure, spectral properties, and lattice energy Elatt. These crystals consist of closed-shell organic ions interacting through relatively strong hydrogen bonds, which leads to Elatt > 200 kJ/mol. Some technical problems that a user of periodic (solid-state) density functional theory (DFT) programs encounters when calculating the properties of these crystals still remain unsolved, for example, the influence of cell parameter optimization on the Elatt value, wave numbers, relative intensity of Raman-active vibrations in the low-frequency region, etc. In this work, various properties of a new two-component carbendazim maleate crystal were experimentally investigated, and the applicability of different DFT functionals and empirical Grimme corrections to the description of the obtained structural and spectroscopic properties was tested. Based on this, practical recommendations were developed for further theoretical studies of multicomponent organic pharmaceutical crystals.

Water adsorption on carbon - A review

Water adsorption on carbonaceous materials has been studied increasingly in the recent years, not only because of its impact on many industrial processes, but also motivated by a desire to understand, at a fundamental level, the distinctive character of directional interactions between water molecules, and between water molecules and other polar groups, such as the functional groups (FGs) at the surfaces of graphene layers. This paper presents an extensive review of recent experimental and theoretical work on water adsorption on various carbonaceous materials, with the aim of gaining a better understanding of how water adsorption in carbonaceous materials relates to the concentration of FGs, their topology (arrangement of the groups) and the structure of the confined space in porous carbons. Arising from this review we are able to propose mechanisms for water adsorption in carbonaceous materials as the adsorbate density increases. The intricate interplay between the roles of FGs and confinement makes adsorption of water on carbon materials very different from that of other simple molecules.

Supramolecular networks in molecular complexes of pyridine boronic acids and polycarboxylic acids: synthesis, structural characterization and fluorescence properties

3- and 4-pyridineboronic acids have been combined with trimesic and pyromellitic acids to give three molecular complexes.

Proton transfer aiding phase transitions in oxalic acid dihydrate under pressure

Oxalic acid dihydrate, an important molecular solid in crystal chemistry, ecology and physiology, has been studied for nearly 100 years now. The most debated issues regarding its proton dynamics have arisen due to an unusually short hydrogen bond between the acid and water molecules. Using combined in situ spectroscopic studies and first-principles simulations at high pressures, we show that the structural modification associated with this hydrogen bond is much more significant than ever assumed. Initially, under pressure, proton migration takes place along this strong hydrogen bond at a very low pressure of 2 GPa. This results in the protonation of water with systematic formation of dianionic oxalate and hydronium ion motifs, thus reversing the hydrogen bond hierarchy in the high pressure phase II. The resulting hydrogen bond between a hydronium ion and a carboxylic group shows remarkable strengthening under pressure, even in the pure ionic phase III. The loss of cooperativity of hydrogen bonds leads to another phase transition at ∼ 9 GPa through reorientation of other hydrogen bonds. The high pressure phase IV is stabilized by a strong hydrogen bond between the dominant CO2 and H2O groups of oxalate and hydronium ions, respectively. These findings suggest that oxalate systems may provide useful insights into proton transfer reactions and assembly of simple molecules under extreme conditions.

Born-Oppenheimer Molecular Dynamics Study on Proton Dynamics of Strong Hydrogen Bonds in Aspirin Crystals, with Emphasis on Differences between Two Crystal Forms

In this study, the proton dynamics of hydrogen bonds for two forms of crystalline aspirin was investigated by the Born-Oppenheimer molecular dynamics (BOMD) method. Analysis of the geometrical parameters of hydrogen bonds using BOMD reveals significant differences in hydrogen bonding between the two crystalline forms of aspirin, Form I and Form II. Analysis of the trajectory for Form I shows spontaneous proton transfer in cyclic dimers, which is absent in Form II. Quantization of the O-H stretching modes allows a detailed discussion on the strength of hydrogen-bonding interactions. The focal point of our study is examination of the hydrogen bond characteristics in the crystal structure and clarification of the influence of hydrogen bonding on the presence of the two crystalline forms of aspirin. In the BOMD method, thermal motions were taken into account. Solving the Schrödinger equation for the snapshots of 2D proton potentials, extracted from MD, gives the best agreement with IR spectra. The character of medium-strong hydrogen bonds in Form I of aspirin was compared with that of weaker hydrogen bonds in aspirin Form II. Two proton minima are present in the potential function for the hydrogen bonds in Form I. The band contours, calculated by using one- and two-dimensional O-H quantization, reflect the differences in the hydrogen bond strengths between the two crystalline forms of aspirin, as well as the strong hydrogen bonding in the cyclic dimers of Form I and the medium-strong hydrogen bonding in Form II.

A comprehensive classification and nomenclature of carboxyl-carboxyl(ate) supramolecular motifs and related catemers: implications for biomolecular systems

Carboxyl and carboxylate groups form important supramolecular motifs (synthons). Besides carboxyl cyclic dimers, carboxyl and carboxylate groups can associate through a single hydrogen bond. Carboxylic groups can further form polymeric-like catemer chains within crystals. To date, no exhaustive classification of these motifs has been established. In this work, 17 association types were identified (13 carboxyl-carboxyl and 4 carboxyl-carboxylate motifs) by taking into account the syn and anti carboxyl conformers, as well as the syn and anti lone pairs of the O atoms. From these data, a simple rule was derived stating that only eight distinct catemer motifs involving repetitive combinations of syn and anti carboxyl groups can be formed. Examples extracted from the Cambridge Structural Database (CSD) for all identified dimers and catemers are presented, as well as statistical data related to their occurrence and conformational preferences. The inter-carboxyl(ate) and carboxyl(ate)-water hydrogen-bond properties are described, stressing the occurrence of very short (strong) hydrogen bonds. The precise characterization and classification of these supramolecular motifs should be of interest in crystal engineering, pharmaceutical and also biomolecular sciences, where similar motifs occur in the form of pairs of Asp/Glu amino acids or motifs involving ligands bearing carboxyl(ate) groups. Hence, we present data emphasizing how the analysis of hydrogen-containing small molecules of high resolution can help understand structural aspects of larger and more complex biomolecular systems of lower resolution.

Complex three-dimensional lanthanide metal–organic frameworks with variable coordination spheres based on pyrazine-2,3,5,6-tetracarboxylate

An open 3-D MOF with complex connectivity of multi-topic pyrazine-2,3,5,6-tetracarboxylate linker and Ln(iii) ions.

Towards full-color-tunable emission of two component Eu(iii)-doped Gd(iii) coordination frameworks by the variation of excitation light

Eu(iii)-doped Gd(iii) pyrazine-2,3,5,6-tetracarboxylate alkali–lanthanide metal–organic frameworks (LnMOF) with a unique (4¹¹·6⁸·8²)(4³·6²·8)(4³) topology present full color-tunable luminescence and white emission by the variation of excitation wavelengths.

Cl···Cl interactions in molecular crystals: insights from the theoretical charge density analysis

The structure, IR harmonic frequencies and intensities of normal vibrations of 20 molecular crystals with the X-Cl···Cl-X contacts of different types, where X = C, Cl, and F and the Cl···Cl distance varying from ~3.0 to ~4.0 Å, are computed using the solid-state DFT method. The obtained crystalline wave functions have been further used to define and describe quantitatively the Cl···Cl interactions via the electron-density features at the Cl···Cl bond critical points. We found that the electron-density at the bond critical point is almost independent of the particular type of the contact or hybridization of the ipso carbon atom. The energy of Cl···Cl interactions, E(int), is evaluated from the linking E(int) and local electronic kinetic energy density at the Cl···Cl bond critical points. E(int) varies from 2 to 12 kJ/mol. The applicability of the geometrical criterion for the detection of the Cl···Cl interactions in crystals with two or more intermolecular Cl···Cl contacts for the unique chlorine atom is not straightforward. The detection of these interactions in such crystals may be done by the quantum-topological analysis of the periodic electron density.

Intermolecular hydrogen bond energies in crystals evaluated using electron density properties: DFT computations with periodic boundary conditions

The hydrogen bond (H-bond) energies are evaluated for 18 molecular crystals with 28 moderate and strong O-H···O bonds using the approaches based on the electron density properties, which are derived from the B3LYP/6-311G** calculations with periodic boundary conditions. The approaches considered explore linear relationships between the local electronic kinetic G(b) and potential V(b) densities at the H···O bond critical point and the H-bond energy E(HB). Comparison of the computed E(HB) values with the experimental data and enthalpies evaluated using the empirical correlation of spectral and thermodynamic parameters (Iogansen, Spectrochim. Acta Part A 1999, 55, 1585) enables to estimate the accuracy and applicability limits of the approaches used. The V(b)-E(HB) approach overestimates the energy of moderate H-bonds (E(HB) < 60 kJ/mol) by ~20% and gives unreliably high energies for crystals with strong H-bonds. On the other hand, the G(b)-E(HB) approach affords reliable results for the crystals under consideration. The linear relationship between G(b) and E(HB) is basis set superposition error (BSSE) free and allows to estimate the H-bond energy without computing it by means of the supramolecular approach. Therefore, for the evaluation of H-bond energies in molecular crystals, the G(b) value can be recommended to be obtained from both density functional theory (DFT) computations with periodic boundary conditions and precise X-ray diffraction experiments.

5,6-Dimethyl-pyrazine-2,3-dicarb-oxy-lic acid

The asymmetric unit of the title compound, C₈H₈N₂O₄, consists of one complete mol­ecule and a second mol­ecule generated by the application of twofold axis. The mean planes of the two carboxyl groups attached to the pyrazine ring at neighboring positions are twisted by 10.8 (1) and 87.9 (1)° in the complete molecule and 43.0 (1)° in the symmetry-generated molecule. The crystal packing features O—H⋯N hydrogen bonds, which link the mol­ecules into layers along [101].

Diisopropyl pyrazine-2,5-dicarboxyl-ate

The mol-ecule of the title compound, C(12)H(16)N(2)O(4), is located on an inversion center. The carboxyl-ate groups are twisted slightly with respect to the pyrazine ring, making a dihedral angle of 6.4 (3)°.

Bond lengths in organic and metal-organic compounds revisited: X-H bond lengths from neutron diffraction data

The number of structures in the Cambridge Structural Database (CSD) has increased by an order of magnitude since the preparation of two major compilations of standard bond lengths in mid-1985. It is now of interest to examine whether this huge increase in data availability has implications for the mean bond-length values published in the late 1980s. Those compilations reported mean X-H bond lengths derived from rather sparse information and for rather few chemical environments. During the intervening years, the number of neutron studies has also increased, although only by a factor of around 2.25, permitting a new analysis of X-H bond-length distributions for (a) organic X = C, N, O, B, and (b) a variety of terminal and homometallic bridging transition metal hydrides. New mean values are reported here and are compared with earlier results. These new overall means are also complemented by an analysis of X-H distances at lower temperatures (T < or = 140 K), which indicates the general level of librational effects in X-H systems. The study also extends the range of chemical environments for which statistically acceptable mean X-H bond lengths can be obtained, although values from individual structures are also collated to further extend the chemical range of this compilation. Updated default 'neutron-normalization' distances for use in hydrogen-bond and deformation-density studies are also proposed for C-H, N-H and O-H, and the low-temperature analysis provides specific values for certain chemical environments and hybridization states of X.

Comparison of cationic, anionic and neutral hydrogen bonded dimers

Short Strong Hydrogen Bonds (SSHBs) play an important role in many fields of physics, chemistry and biology. Since it is known that SSHBs exist in many biological systems, the role of hydrogen bonding motifs has been particularly interesting in enzyme catalysis, bio-metabolism, protein folding and proton transport phenomena. To explore the characteristic features of neutral, anionic and cationic hydrogen bonds, we have carried out theoretical studies of diverse homogeneous and heterogeneous hydrogen bonded dimers including water, peroxides, alcohols, ethers, aldehydes, ketones, carboxylic acids, anhydrides, and nitriles. Geometry optimization and harmonic frequency calculations are performed at the levels of Density Functional Theory (DFT) and Møller-Plesset second order perturbation (MP2) theory. First principles Car-Parrinello molecular dynamics (CPMD) simulations are performed to obtain IR spectra derived from velocity- and dipole-autocorrelation functions. We find that the hydrogen bond energy is roughly inversely proportional to the fourth power of the r(O/N-H) distance. Namely, the polarization of the proton accepting O/N atom by the proton-donating H atom reflects most of the binding energy in these diverse cation/anion/neutral hydrogen bonds. The present study gives deeper insight into the nature of hydrogen-bonded dimers including SSHBs.

Tautomerism and isotopic multiplets in the 13C NMR spectra of partially deuterated 3-arylpyrimido[4,5-c]pyridazine-5,7(6H,8H)-diones and their sulfur analogs--evidence for elucidation of the structure backbone and tautomeric forms

The tautomerism of the synthesized 3-arylpyrimido[4,5-c]pyridazine-5,7(6H,8H)-diones (1a-d) and 3-aryl-7-thioxo-7,8-dihydro-6H-pyrimido[4,5-c]pyridazine-5-ones (2a-d) was studied in dimethyl sulfoxide (DMSO)-d(6). (1)H NMR spectra of 1a-d showed a clustered water molecule in the structure backbone that is attached by strong intermolecular H bonding. The relation between the temperature and H bonding of the clustered water molecule with 1a was also studied as representative. The relation between the electronegativity (chi) of the substituent on phenyl ring and the chemical shifts of clustered water protons in 1a-d was also studied. All of 1a-d and also 2d compounds existed in lactam (I) form, whereas 2a-c compounds have two distinguished tautomers in DMSO-d(6) [lactam (I) and lactim (II) forms]. The solvent-substrate proton exchange was examined in compounds 1a-d and 2a-d by adding one drop of D(2)O. All compounds (except 1d) showed proton/deuterium exchange of the clustered water protons in DMSO by adding one drop of D(2)O. Some compounds (but not all of them) that are easily soluble in DMSO-d(6) containing D(2)O showed isotopic splitting (beta-isotope effect) in their (13)C NMR spectra. Among them, compound 1a was the best evidence to help the spectral assignments and structure determination of predominant tautomer by carbon-13 splitting (beta-isotope effect).

Crystal engineering: a holistic view

Crystal engineering, the design of molecular solids, is the synthesis of functional solid-state structures from neutral or ionic building blocks, using intermolecular interactions in the design strategy. Hydrogen bonds, coordination bonds, and other less directed interactions define substructural patterns, referred to in the literature as supramolecular synthons and secondary building units. Crystal engineering has considerable overlap with supramolecular chemistry, X-ray crystallography, materials science, and solid-state chemistry and yet it is a distinct discipline in itself. The subject goes beyond the traditional divisions of organic, inorganic, and physical chemistry, and this makes for a very eclectic blend of ideas and techniques. The purpose of this Review is to highlight some current challenges in this rapidly evolving subject. Among the topics discussed are the nature of intermolecular interactions and their role in crystal design, the sometimes diverging perceptions of the geometrical and chemical models for a molecular crystal, the relationship of these models to polymorphism, knowledge-based computational prediction of crystal structures, and efforts at mapping the pathway of the crystallization reaction.

QTAIM study of strong H-bonds with the O-H...a fragment (A=O, N) in three-dimensional periodical crystals

The relationship between the d(H...A) distance (A=O, N) and the topological properties at the H...A bond critical point of 37 strong (short) hydrogen bonds occurring in 26 molecular crystals are analyzed using the quantum theory of atoms in molecules (QTAIM). Ground-state wave functions of the three-dimensional periodical structures representing the accurate experimental geometries calculated at the B3LYP/6-31G** level of approximation were used to obtain the QTAIM electron density characteristics. The use of an electron-correlated method allowed us to reach the quantitatively correct values of electron density rhob at the H...A bond critical point. However, quite significant differences can appear for small absolute values of the Laplacian (<0.5 au). The difference between the H...O and H...N interactions is described using the rhob versus d(H...A) dependence. It is demonstrated that the values of parameters in this dependence are defined by the nature of the heavy atom forming the H...A bond. An intermediate (or transit) region separating the shared and closed-shell interactions is observed for the H-bonded crystals in which the bridging proton can move from one heavy atom to another. The crystalline environment changes the location of the bridging proton in strong H-bonded systems; however, the d(O-H)/d(H...O) ratio is approximately the same for both the gas-phase complexes and molecular crystals with a linear or near-linear O-H...O bond.

Interactions of apo cytochrome C with alternating copolymers of maleic acid and alkene

Apo cytochrome c (apo cyt c) tends to aggregate at alkali pH. Poly(isobutylene-alt-maleic acid) (PIMA) is soluble molecularly, whereas poly(1-tetradecene-alt-maleic acid) (PTMA) forms particles that tend to dissociate by increasing pH and decreasing concentration. Dynamic light scattering and surface plasmon resonance are used to investigate the interactions of PIMA and PTMA with apo cyt c at different pH values to understand the mechanism of the interactions. When the positive or negative charges are in excess, the copolymer-protein complex particles can be stabilized by the charges on the surface. When the ratio of the positive to negative charges is close to the stoichiometric value, precipitation occurs. At pH 11.8, both PTMA and apo cyt c carry negative charges, but the hydrophobic interaction makes them form complexes. A competition exists between the interaction of the copolymer with apo cyt c and the self-aggregation of PTMA or apo cyt c alone. The interaction of PIMA or PTMA with apo cyt c at neutral and alkali pH destroys the aggregation of PTMA or apo cyt c and forms new complex particles.

Structural Transformation of Apocytochrome c Induced by Alternating Copolymers of Maleic Acid and Alkene

Apocytochrome c interacts with two copolymers: poly(isobutylene-alt-maleic acid) (PIMA) and poly(1-tetradecene-alt-maleic acid) (PTMA). The interaction leads to apocytochrome c, a conformational change from random coil to alpha-helical structure. The alpha-helix content is influenced by the copolymer concentration, the length of alkyl chain of the copolymers, and pH of the medium. The electrostatic attraction between the copolymer and protein is an indispensable factor for the folding of the protein at acid pH. The hydrophobic interaction is an important factor over the entire pH range, especially when both the copolymer and protein carry negative charges at alkaline pH. The electrostatic and hydrophobic attractions between the copolymer and protein exclude water molecules, promoting the formation of hydrogen bonds within the helical structure. On the other hand, the hydrogen bonds formed between the ionized carboxyl of the copolymer and the amide of the protein partly restrain the formation of hydrogen bonds within the helical structure when the copolymer concentration is higher at pH 6.5 and 10.5.

Short hydrogen bonds in proteins

Short hydrogen bonds are present in many chemical and biological systems. It is well known that these short hydrogen bonds are found in the active site of enzymes and aid enzyme catalysis. This study aims to systematically characterize all short hydrogen bonds from a nonredundant dataset of protein structures. The study has revealed that short hydrogen bonds are commonly found in proteins and are widely present in different regions of the protein chain, such as the backbone or side chain, and in different secondary structural regions such as helices, strands and turns. The frequency of occurrence of donors and acceptors from the charged side chains as well as from the neutral backbone atoms is equally high. This suggests that short hydrogen bonds in proteins occur either due to increased strength or due to geometrical constraints and this has been illustrated from several examples.