Dioxygen-halogen bonding exemplified by crystalline peroxosolvates of N,N'-bis(haloacetyl) bispidines

The halogen bonding in molecular crystals and supramolecular assemblies has been widely investigated. Special attention is given to the molecular structures capable of simultaneously exhibiting different types of non-covalent interactions, including conventional hydrogen bonds and halogen bonds. This paper systematically analyzes crystalline peroxosolvates of bispidine-based bis-amide derivatives, containing haloacetic acid residues, namely previously reported 1,1'-(1,5-dimethyl-3,7-diazabicyclo[3.3.1]nonane-3,7-diyl)bis(2-iodooethanone) peroxosolvate C13H20I2N2O2·H2O2 (1) and four new crystalline compounds, 1,1'-(1,5-dimethyl-3,7-diazabicyclo[3.3.1]nonane-3,7-diyl)bis(2-bromoethanone) peroxosolvate C13H20Br2N2O2·H2O2 (2), 1,1'-(9-hydroperoxy-9-hydroxy-1,5-dimethyl-3,7-diazabicyclo[3.3.1]nonane-3,7-diyl)bis(2-iodoethanone) peroxosolvate C13H20I2N2O5·0.5H2O2 (3), 1,1'-(9-hydroperoxy-9-hydroxy-1,5-dimethyl-3,7-diazabicyclo[3.3.1]nonane-3,7-diyl)bis(2-bromoethanone) peroxosolvate C13H20Br2N2O5·H2O2 (4), and 1,1'-(9-hydroperoxy-9-hydroxy-1,5-dimethyl-3,7-diazabicyclo[3.3.1]nonane-3,7-diyl)bis(2-chloroethanone) peroxosolvate C13H20Cl2N2O5·H2O2 (5). Compounds 2-5 were synthesized for the first time and their crystal structures were determined by single-crystal X-ray diffractometry (SCXRD). To the best of our knowledge, 3-5 are unprecedented crystalline hydrogen peroxide adducts of organic hydroperoxides (R-OOH). Short intermolecular contacts between halogen and hydroperoxo oxygen atoms were found in 1-3. The halogen bonding of C-I(Br) fragments with dioxygen species in compounds 1-3 as well as in the previously reported cocrystal of diacetone diperoxide with triodotrinitrobenzene (6) was identified through reduced density gradient analysis, Hirshfeld surface analysis, and Bader analysis of crystalline electron density. The interactions were quantified using the electron density topological properties acquired from the periodic DFT calculations and evaluated to lie in the range of 9-19 kJ mol-1. A distinctive spectral feature was revealed for this type of interaction, involving a red shift of the characteristic O-O stretching vibration by about 6 cm-1, which appeared in IR spectra as a narrow low-intensity band in the region 837-872 cm-1.

Characterizing hydrogen bonds in crystalline form of guanidinium salicylate in the terahertz range

For pharmaceutical compounds with poor solubility, there is an effective method to address this dilemma without tampering their intrinsic chemical properties by forming weak hydrogen bonds. Guanidinium salicylate, which is a typical pharmaceutical salt with a complex crystal structure, was systematically investigated by terahertz time-domain spectroscopy combined with density functional theory in order to obtain the complete information of weak hydrogen bonds. As a result of the influence of weak hydrogen bonds, there are substantial differences between guanidinium salicylate and its parent molecule (salicylic acid) in the experimental fingerprint spectra in the range of 0.2-2.5 THz, such as the number, amplitude and frequency positions of absorption peaks. With the help of isolated molecule density functional theory calculations, the possible sites of weak hydrogen bonds were determined by natural bond orbital analysis. It can be concluded that there is an intricate hydrogen bond network due to the polar distribution of molecular electrostatic potential. Furthermore, all THz absorption peaks were assigned to their corresponding vibrational modes and the complete information of the related hydrogen bonds (including type, role, angle, and bond length) was determined by using dispersion-corrected density functional theory. The results laid a good foundation for further study on the enhancement of solubility of pharmaceutical salts by forming weak hydrogen bonds.

Substituent effect on the packing architecture of adamantane and memantine derivatives of sulfonamide molecular crystals

A specific number of structures described in this paper with adamantane and memantine fragments have been synthesised and characterised.

Local-structure effects on 31P NMR chemical shift tensors in solid state

The effect of the local structure on the ³¹P NMR chemical shift tensor (CST) has been studied experimentally and simulated theoretically using the density functional theory gauge-independent-atomic-orbital approach. It has been shown that the dominating impact comes from a small number of noncovalent interactions between the phosphorus-containing group under question and the atoms of adjacent molecules. These interactions can be unambiguously identified using the Bader analysis of the electronic density. A robust and computationally effective approach designed to attribute a given experimental ³¹P CST to a certain local morphology has been elaborated. This approach can be useful in studies of surfaces, complex molecular systems, and amorphous materials.

Cyclic dipeptide peroxosolvates: first direct evidence for hydrogen bonding between hydrogen peroxide and a peptide backbone

The crystal structures of cyclic dipeptide peroxosolvates provide valuable insight into the non-redox interaction of hydrogen peroxide with the peptide backbone.

Specific features of supramolecular organisation and hydrogen bonding in proline cocrystals: a case study of fenamates and diclofenac

New zwitterionic cocrystals of fenamate drugs and diclofenac with the naturally occurring amino acid l-proline have been obtained and thoroughly characterised by a variety of experimental and theoretical techniques.

Structural and energetic aspects of adamantane and memantine derivatives of sulfonamide molecular crystals: experimental and theoretical characterisation

A number of new sulfonamide compounds with adamantane and memantine fragments were synthesised and characterised.

Quantification of Conventional and Nonconventional Charge-Assisted Hydrogen Bonds in the Condensed and Gas Phases

Charge-assisted hydrogen bonds (CAHBs) play critical roles in many systems from biology through to materials. In none of these areas has the role and function of CAHBs been explored satisfactorily because of the lack of data on the energy of CAHBs in the condensed phases. We have, for the first time, quantified three types of CAHBs in both the condensed and gas phases for 1-(2'-hydroxylethyl)-3-methylimidazolium acetate ([C2OHmim][OAc]). The energy of conventional OH···[OAc](-) CAHBs is ∼10 kcal·mol(-1), whereas nonconventional C(sp2)H···[OAc](-) and C(sp3)H···[OAc](-) CAHBs are weaker by ∼5-7 kcal·mol(-1). In the gas phase, the strength of the nonconventional CAHBs is doubled, whereas the conventional CAHBs are strengthened by <20%. The influence of cooperativity effects on the ability of the [OAc](-) anion to deprotonate the imidazolium cation is evaluated. The ability to quantify CAHBs in the condensed phase on the basis of easier accessible gas-phase estimates is highlighted.

Influence of Secondary Interactions on the Structure, Sublimation Thermodynamics, and Solubility of Salicylate:4-Hydroxybenzamide Cocrystals. Combined Experimental and Theoretical Study

Cocrystal screening of 4-hydroxybenzamide with a number of salicylates (salicylic acid, SA; 4-aminosalicylic acid, PASA; acetylsalicylic acid, ASA; and salicylsalicylic acid, SSA) was conducted to confirm the formation of two cocrystals, [SA+4-OHBZA] (1:1) and [PASA+4-OHBZA] (1:1). Their structures were determined using single-crystal X-ray diffraction, and the hydrogen-bond network topology was studied. Thermodynamic characteristics of salicylic acid cocrystal sublimation were obtained experimentally. It was proved that PASA cocrystallization with 4-OHBZA makes the drug more stable and prevents the irreversible process of decarboxylation of PASA resulting in formation of toxic 3-aminophenol. The pattern of non-covalent interactions in the cocrystals is described quantitatively using solid-state density functional theory followed by Bader analysis of the periodic electron density. It has been found that the total energy of secondary interactions between synthon atoms and the side hydroxyl group of the acid molecule in [SA+4-OHBZA] (1:1) and [PASA+4-OHBZA] (1:1) cocrystals is comparable to the energy of the primary acid-amide heterosynthon. The theoretical value of the sublimation enthalpy of [SA+4-OHBZA], 231 kJ/mol, agrees fairly well with the experimental one, 272 kJ/mol. The dissolution experiments with [SA+4-OHBZA] have proved that the relatively large cocrystal stability in relation to the stability of its components has a negative effect on the dissolution rate and equilibrium solubility. The [PASA+4-OHBZA] (1:1) cocrystal showed an enhancement of apparent solubility compared to that of the corresponding pure active pharmaceutical ingredient, while their intrinsic dissolution rates are comparable.

Saccharin salts of biologically active hydrazone derivatives

Crystal structures, solubility and formation thermodynamics of saccharin salts with biologically active hydrazone derivatives were investigated.

Salicylamide cocrystals: screening, crystal structure, sublimation thermodynamics, dissolution, and solid-state DFT calculations

A new cocrystal of 2-hydroxybenzamide (A) with 4-acetamidobenzoic acid (B) has been obtained by the DSC screening method. Thermophysical analysis of the aggregate [A:B] has been conducted and a fusion diagram has been plotted. Cocrystal formation from melts was studied by using thermomicroscopy. A cocrystal single-crystal was grown and its crystal structure was determined. The pattern of noncovalent interactions has been quantified using the solid-state DFT computations coupled with the Bader analysis of the periodic electron density. The sublimation processes of A-B cocrystal have been studied and its thermodynamic functions have been calculated. The classical method of substance transfer by inert gas-carrier was chosen to investigate sublimation processes experimentally. The lattice energy is found to be 143 ± 4 kJ/mol. It is lower than the sum of the corresponding values of the cocrystal pure components. The theoretical value of the lattice energy, 156 kJ/mol, is in reasonable agreement with the experimental one. A ternary phase diagram of solubility (A-B-ethanol) has been plotted and the areas with solutions for growing thermodynamically stable cocrystals have been determined.

Capping amyloid β-sheets of the tau-amyloid structure VQIVYK with hexapeptides designed to arrest growth. An ONIOM and density functional theory study

We present ONIOM calculations using density functional theory (DFT) as the high and AM1 as the medium level that explore the abilities of different hexapeptide sequences to terminate the growth of a model for the tau-amyloid implicated in Alzheimer’s disease. We delineate and explore several design principles (H-bonding in the side chains, using antiparallel interactions on the growing edge of a parallel sheet, using all-d residues to form rippled interactions at the edge of the sheet, and replacing the H-bond donor N–H’s that inhibit further growth) that can be used individually and in combination to design such peptides that will have a greater affinity for binding to the parallel β-sheet of acetyl-VQIVYK-NHCH₃ than the natural sequence and will prevent another strand from binding to the sheet, thus providing a cap to the growing sheet that arrests further growth. We found peptides in which the Q is replaced by an acetyllysine (aK) residue to be particularly promising candidates, particularly if the reverse sequence (KYVIaKV) is used to form an antiparallel interaction with the sheet.

First Example of Hydrogen Bonding to Platinum Hydride

The interaction between various proton donors (indole, CF3CH2OH, (CF3)2CHOH, (CF3)3COH) and (NNC)PtH hydrido complex (NNC-H=6-(1,1'-dimethylbenzyl)-2,2'-bipyridine) was investigated through low-temperature IR and NMR spectroscopy in combination with density functional theory calculations at the M06 level of theory. The experiment shows formation of very weak hydrogen bonded complexes (Δ HºHB ca. -1.0 kcal mol-1), which undergo subsequent proton transfer surprisingly easy. Computational analysis of the hydrogen bonded complexes geometry, electronic parameters (obtained by NBO and AIM analysis), and orbital interaction energies shows that all the complexes are better described as bonded to the metal atom. At that in case of weak alcohols (CH3OH, TFE) there is also the additional interaction with the hydride ligand. These computational results allow explaining the observed experimental trends and give the first example of hydrogen bonding to a metal atom in the presence of hydride ligand.

Aqueous solvation of polyalanine α-helices with specific water molecules and with the CPCM and SM5.2 aqueous continuum models using density functional theory

We present density functional theory (DFT) calculations at the X3LYP/D95(d,p) level on the solvation of polyalanine α-helices in water. The study includes the effects of discrete water molecules and the CPCM and AMSOL SM5.2 solvent continuum model both separately and in combination. We find that individual water molecules cooperatively hydrogen-bond to both the C- and N-termini of the helix, which results in increases in the dipole moment of the helix/water complex to more than the vector sum of their individual dipole moments. These waters are found to be more stable than in bulk solvent. On the other hand, individual water molecules that interact with the backbone lower the dipole moment of the helix/water complex to below that of the helix itself. Small clusters of waters at the termini increase the dipole moments of the helix/water aggregates, but the effect diminishes as more waters are added. We discuss the somewhat complex behavior of the helix with the discrete waters in the continuum models.

Comparative Study of Vibrational Spectra of Two Bioactive Natural Products Lupeol and Lupenone Using MM/QM Method

This work deals with the theoretical study on the molecular structure and vibrational spectra of two well-known natural products: lupeol and lupenone. The spectra were interpreted with the aid of normal mode analysis following full-structure optimization carried out with the hybrid two-level ONIOM (B3LYP/6-31G: PM3) method. A detailed interpretation of the infrared spectra of Lupeol and Lupenone is also reported in the present work. The similarities and differences between the vibrational spectra of the two molecules studied have been highlighted. The scaled theoretical wave numbers are in perfect agreement with the experimental values. The thermodynamic calculations related to the title compounds were also performed at B3LYP/6-31G: PM3 level of theory. Quantum chemical calculations have been carried out to understand the dynamical behavior of the bioactive molecules Lupeol and Lupenone.

H-bond network in amino acid cocrystals with H2O or H2O2. The DFT study of serine-H2O and serine-H2O2

The structure, IR spectrum, and H-bond network in the serine-H(2)O and serine-H(2)O(2) crystals were studied using DFT computations with periodic boundary conditions. Two different basis sets were used: the all-electron Gaussian-type orbital basis set and the plane wave basis set. Computed frequencies of the IR-active vibrations of the titled crystals are quite different in the range of 10-100 cm(-1). Harmonic approximation fails to reproduce IR active bands in the 2500-2800 frequency region of serine-H(2)O and serine-H(2)O(2). The bands around 2500 and 2700 cm(-1) do exist in the anharmonic IR spectra and are caused by the first overtone of the OH bending vibrations of H(2)O and a combination vibration of the symmetric and asymmetric bendings of H(2)O(2). The quantum-topological analysis of the crystalline electron density enables us to describe quantitatively the H-bond network. It is much more complex in the title crystals than in a serine crystal. Appearance of water leads to an increase of the energy of the amino acid-amino acid interactions, up to ~50 kJ/mol. The energy of the amino acid-water H-bonds is ~30 kJ/mol. The H(2)O/H(2)O(2) substitution does not change the H-bond network; however, the energy of the amino acid-H(2)O(2) contacts increases up to 60 kJ/mol. This is caused by the fact that H(2)O(2) is a much better proton donor than H(2)O in the title crystals.

Comparison of β-sheets of capped polyalanine with those of the tau-amyloid structures VQIVYK and VQIINK. A density functional theory study

We present ONIOM calculations using B3LYP/d95(d,p) as the high and AM1 as the low level on parallel β-sheets containing from two to ten strands of Ac-VQIVYK-NHMe and Ac-VQIINK-NHMe, as well as both parallel and antiparallel Ac-AAAAAA-NHMe. We find that the first two sequences form more stable sheets due to the additional H-bonding between the Q's in the side chains of both and the N's in the side chain of Ac-VQIINK-NHMe. However, the H-bonds in the amyloid chains are significantly weakened by attractive strain, which prevents all the interstrand H-bonds from achieving their optimal geometries simultaneously and requires high distortion energies for the individual strands in the sheets. The antiparallel Ac-AAAAAA-NHMe's are generally more stable and more cooperative than the parallel sheets, principally due to the higher distortion energies of the latter.

A model proton-transfer system in the condensed phase: NH4(+)OOH(-), a crystal with short intermolecular H-bonds

The crystal structure of NH(4)(+)OOH(-) is determined from single-crystal x-ray data obtained at 150 K. The crystal belongs to the space group P2(1)/c and has four molecules in a unit cell. The structure consists of discrete NH(4)(+) and OOH(-) ions. The OOH(-) ions are linked by short hydrogen bonds (2.533 Å) to form parallel infinite chains. The ammonium ions form links between these chains (the N⋯O distances vary from 2.714 to 2.855 Å) giving a three-dimensional network. The harmonic IR spectrum and H-bond energies are computed at the Perdew-Burke-Ernzerhof (PBE)/6-31G(∗∗) level with periodic boundary conditions. A detailed analysis of the shared (bridging) protons' dynamics is obtained from the CPMD simulations at different temperatures. PBE functional with plane-wave basis set (110 Ry) is used. At 10 K the shared proton sits near the oxygen atom, only a few proton jumps along the chain are detected at 70 K while at 270 K numerous proton jumps exist in the trajectory. The local-minimum structure of the space group Cc is localized. It appears as a result of proton transfer along a chain. This process is endothermic (∼2 kJ/mol) and is described as P2(1)/c↔2Cc. The computed IR spectrum at 10 K is close to the harmonic one, the numerous bands appear at 70 K while at 270 K it shows a very broad absorption band that covers frequencies from about 1000 to 3000 cm(-1). The advantages of the NH(4)(+)OOH(-) crystal as a promising model for the experimental and DFT based molecular dynamics simulation studies of proton transfer along the chain are discussed.

Reactivity of a phospholipid monolayer model under periodic boundary conditions: a density functional theory study of the Schiff base formation between phosphatidylethanolamine and acetaldehyde

A mechanism for the formation of the Schiff base between an acetaldehyde and an amine-phospholipid monolayer model based on Dmol3/density functional theory calculations under periodic boundary conditions was constructed. This is the first time such a system has been modeled to examine its chemical reactivity at this computation level. Each unit cell contains two phospholipid molecules, one acetaldehyde molecule, and nine water molecules. One of the amine-phospholipid molecules in the cell possesses a neutral amino group that is used to model the nucleophilic attack on the carboxyl group of acetaldehyde, whereas the other has a charged amino group acting as a proton donor. The nine water molecules form a hydrogen bond network along the polar heads of the phospholipids that facilitates very fast proton conduction at the interface. Using periodic boundary conditions afforded proton transfer between different cells. The reaction takes place in two steps, namely, (1) formation of a carbinolamine and (2) its dehydration to the Schiff base. The carbinolamine is the primary reaction intermediate, and dehydration is the rate-determining step of the process, consistent with available experimental evidence for similar reactions. On the basis of the results, the cell membrane surface environment may boost phospholipid glycation via a neighboring catalyst effect.

Investigation of protein conformation and interactions with salts via X-ray absorption spectroscopy

Nitrogen K-edge spectra of aqueous triglycine were measured using liquid microjets, and the effects of Hofmeister-active salts on the spectra were observed. Spectra simulated using density functional theory, sampled from room temperature classical molecular dynamics trajectories, capture all major features in the measured spectra. The spectrum of triglycine in water is quite similar to that in the presence of chaotropic sodium bromide (and other halides), which raises the solubility of proteins. However, a new feature is found when kosmotropic Na(2)SO(3), which lowers solubility, is present; this feature results from excitations of the nitrogen atom in the terminal amino group of triglycine. Both direct interactions between this salt and the protonated amino terminus, as well as corresponding changes in the conformational dynamics of the system, contribute to this new feature. These molecular measurements support a different mechanism for the Hofmeister effect than has previously been suggested based on thermodynamic measurements. It is also shown that near edge X-ray absorption fine structure (NEXAFS) is sensitive to strong direct interaction between certain salts and charged peptides. However, by investigating the sensitivity of NEXAFS to the extreme structural differences between model beta-sheets and alpha-helices, we conclude that this technique is relatively insensitive to secondary structure of peptides and proteins.

Density functional study of the proton transfer effect on vibrations of strong (short) intermolecular O-H...N/O-...H-N+ hydrogen bonds in aprotic solvents

The structure and spectroscopic properties of the 1:1 complexes of substituted pyridines with benzoic acid and phenol derivatives in aprotic solvents are studied using B3LYP functional combined with the polarizable continuum model approximation. Two extreme structures are investigated: the state without (HB) and with proton transfer (PT). In the presence of an external electric field the O...N distance is contracted and the PT state does appear. The PT state of both the pyridine-benzoic and the pyridine-phenol complexes displays the only IR-active band in the 2800-1800 frequency region, which is located around 2000 cm(-1). However, the nature of the band is different for these two complexes. In the pyridine-benzoic acid complex it is practically a pure stretching vibration of the HN(+) group, while in the pyridine-phenol complex it is the mixed vibration of the bridging proton. A specific feature of the PT state in the pyridine-phenol complex is an IR-intensive band near 600 cm(-1), associated with the asymmetric stretching vibrations of the O(-)...HN(+) fragment. Its intensity is reciprocally proportional to the O...N distance. The appearance of this band provides an efficient criterion to differentiate between the HB and PT states of the 1:1 complexes of phenols with pyridines in aprotic solvents.

Computational study on the conformation and vibration frequencies of β-sheet of ε-polylysine in vacuum

Two oligomers, each containing 3 l-lysine residues, were used as model molecules for the simulation of the β-sheet conformation of ɛ-polylysine (ɛ-PLL) chains. Their C terminals were capped with ethylamine and N terminals were capped with α-l-aminobutanoic acid, respectively. The calculations were carried out with the hybrid two-level ONOIM (B3LYP/6-31G:PM3) computational chemistry method. The optimized conformation was obtained and IR frequencies were compared with experimental data. The result indicated that the two chains were winded around each other to form a distinct cyclohepta structure through bifurcated hydrogen bonds. The groups of amide and α-amidocyanogen coming from one chain and the carbonyl group from the other chain were involved in the cyclohepta structure. The bond angle of the bifurcated hydrogen bonds was 66.6°. The frequency analysis at ONIOM [B3LYP/6-31G (d):PM3] level showed the IR absorbances of the main groups, such as the amide and amidocyanogen groups, were in accordance with the experimental data.