Control of Molecular Conformation and Crystal Packing of Biphenyl Derivatives

This Review presents a discussion of the conformation of biphenyl derivatives in different chemical environments. The interplay between aromatic stabilization and steric repulsion, normally considered to explain the conformation of the molecule, is contrasted with the interpretation provided by models not based on molecular orbitals. The electronic control of conformation by means of appropriate hydrogen substitution is discussed by examples taken from chemistry and molecular electronics. Supramolecular synthons involving biphenyl are critically analyzed in terms of the molecular conformation, crystal packing and intermolecular forces. Some directions for future research on the control of the conformation of biphenyls are also presented.

Two Faces of Water in the Formation and Stabilization of Multicomponent Crystals of Zwitterionic Drug-Like Compounds

Two new hydrated multicomponent crystals of zwitterionic 2-aminonicotinic acid with maleic and fumaric acids have been obtained and thoroughly characterized by a variety of experimental (X-ray analysis and terahertz Raman spectroscopy) and theoretical periodic density functional theory calculations, followed by Bader analysis of the crystalline electron density) techniques. It has been found that the Raman-active band in the region of 300 cm−1 is due to the vibrations of the intramolecular O-H...O bond in the maleate anion. The energy/enthalpy of the intermolecular hydrogen bonds was estimated by several empirical approaches. An analysis of the interaction networks reflects the structure-directing role of the water molecule in the examined multicomponent crystals. A general scheme has been proposed to explain the proton transfer between the components during the formation of multicomponent crystals in water. Water molecules were found to play the key role in this process, forming a “water wire” between the COOH group of the dicarboxylic acid and the COO– group of the zwitterion and the rendering crystal lattice of the considered multicomponent crystals.

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.

London dispersion-driven hetero-aryl-aryl interactions in 1,2-diaryldisilanes

Non-covalent aryl-aryl interactions for the molecular structures of three 1,2-diaryltetramethyldisilanes X5C6-(SiMe2)2-C6Y5 (X ≠ Y; X, Y = H, F, Cl) were studied by single crystal X-ray and gas electron diffraction as well as by SAPT calculations. Aryl-aryl interactions cause all to adopt exclusively rather untypical eclipsed syn-conformers in the gas phase, and C6F5-(SiMe2)2-C6Cl5 also in the solid state.

Interplay of π-stacking and inter-stacking interactions in two-component crystals of neutral closed-shell aromatic compounds: periodic DFT study

This paper bridges the gap between high-level ab initio computations of gas-phase models of 1 : 1 arene–arene complexes and calculations of the two-component (binary) organic crystals using atom–atom potentials. The studied crystals consist of electron-rich and electron-deficient compounds, which form infinite stacks (columns) of heterodimers. The sublimation enthalpy of crystals has been evaluated by DFT periodic calculations, while intermolecular interactions have been characterized by Bader analysis of the periodic electronic density. The consideration of aromatic compounds without a dipole moment makes it possible to reveal the contribution of quadrupole–quadrupole interactions to the π-stacking energy. These interactions are significant for heterodimers formed by arenes with more than 2 rings, with absolute values of the traceless quadrupole moment (Q_zz) larger than 10 D Å. The further aggregation of neighboring stacks is due to the C–H⋯F interactions in arene/perfluoroarene crystals. In crystals consisting of arene and an electron-deficient compound such as pyromellitic dianhydride, aggregation occurs due to the C–H⋯O interactions. The C–H⋯F and C–H⋯O inter-stacking interactions make the main contribution to the sublimation enthalpy, which exceeds 150 kJ mol⁻¹ for the two-component crystals formed by arenes with more than 2 rings.

The interplay of π-stacking and inter-stacking interactions in two-component organic crystals without conventional hydrogen bonds.

Solving the enigma of weak fluorine contacts in the solid state: a periodic DFT study of fluorinated organic crystals

The nature and strength of weak interactions with organic fluorine in the solid state are revealed by periodic density functional theory (periodic DFT) calculations coupled with experimental data on the structure and sublimation thermodynamics of crystalline organofluorine compounds. To minimize other intermolecular interactions, several sets of crystals of perfluorinated and partially fluorinated organic molecules are considered. This allows us to establish the theoretical levels providing an adequate description of the metric and electron-density parameters of the C–F⋯F–C interactions and the sublimation enthalpy of crystalline perfluorinated compounds. A detailed comparison of the C–F⋯F–C and C–H⋯F–C interactions is performed using the relaxed molecular geometry in the studied crystals. The change in the crystalline packing of aromatic compounds during their partial fluorination points to the structure-directing role of C–H⋯F–C interactions due to the dominant electrostatic contribution to these contacts. C–H⋯F–C and C–H⋯O interactions are found to be identical in nature and comparable in energy. The factors that determine the contribution of these interactions to the crystal packing are revealed. The reliability of the results is confirmed by considering the superposition of the electrostatic potential and electron density gradient fields in the area of the investigated intermolecular interactions.

The nature and strength of weak C–H⋯F–C and C–F⋯F–C interactions and their role in organofluorine molecular crystals were studied using periodic DFT coupled with CSD data mining and experimental sublimation enthalpies.

CrystalExplorer model energies and energy frameworks: extension to metal coordination compounds, organic salts, solvates and open-shell systems

The application domain of accurate and efficient CE-B3LYP and CE-HF model energies for intermolecular interactions in molecular crystals is extended by calibration against density functional results for 1794 molecule/ion pairs extracted from 171 crystal structures. The mean absolute deviation of CE-B3LYP model energies from DFT values is a modest 2.4 kJ mol-1 for pairwise energies that span a range of 3.75 MJ mol-1. The new sets of scale factors determined by fitting to counterpoise-corrected DFT calculations result in minimal changes from previous energy values. Coupled with the use of separate polarizabilities for interactions involving monatomic ions, these model energies can now be applied with confidence to a vast number of molecular crystals. Energy frameworks have been enhanced to represent the destabilizing interactions that are important for molecules with large dipole moments and organic salts. Applications to a variety of molecular crystals are presented in detail to highlight the utility and promise of these tools.

The Role of the DNA Backbone in Minor-Groove Ligand Binding

Molecular recognition between ligands and nucleic acids plays a key role in therapeutic activity. Some molecules interact with DNA in a nonbonded manner through intercalation or through the DNA grooves. The recognition of minor-groove binders is attributed to a set of hydrogen-bonding interactions between the binders and the hydrogen-bond-acceptor groups on the groove floor and walls. It is commonly considered that interactions with the sugar groups of the DNA backbone are insignificant and do not contribute to the binding affinity or the specificity. However, our group has found that the deoxyribose rings have a central function in the recognition and the intercalation of metal complexes into DNA. Herein, we determined the specific interactions between the binder CGP 40215A and the minor-groove atoms, based on the local properties of electron density. We found that specific interactions between the deoxyribose moiety within the backbone of DNA and the binder are essential for molecular recognition, and they are responsible for one third of the interaction energy.