A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

An integrated approach combining experimental, informatics and energetic methods for solid form derisking of PF-06282999

The landscapes of observed and predicted three-dimensional crystal packing arrangements of small-molecule drug candidates can be complex. The possible appearance of a more thermodynamically stable solid form during drug development has led to the digital workflow of informatics-based risk assessments, named a Solid Form Health Check. Herein, we describe the use of a combined approach consisting of experiments, informatics together with energetic calculations in analysis of four competing polymorphs of PF-06282999, a myeloperoxidase (MPO) inhibitor with conformational flexibility and multiple plausible hydrogen bond networks. This combined approach offered a comprehensive understanding of the solid form structure, properties, and performance, ensuring robust solid form derisking and selection.

The Persistence of Hydrogen Bonds in Pyrimidinones: From Solution to Crystal

Pyrimidinone scaffolds are present in a wide array of molecules with synthetic and pharmacological utility. The inherent properties of these compounds may be attributed to intermolecular interactions analogous to the interactions that molecules tend to establish with active sites. Pyrimidinones and their fused derivatives have garnered significant interest due to their structural features, which resemble nitrogenous bases, the foundational building blocks of DNA and RNA. Similarly, pyrimidinones are predisposed to forming N-H···O hydrogen bonds akin to nitrogenous bases. Given this context, this study explored the supramolecular features and the predisposition to form hydrogen bonds in a series of 18 substituted 4-(trihalomethyl)-2(1H)-pyrimidinones. The formation of hydrogen bonds was observed in solution via nuclear magnetic resonance (NMR) spectroscopy experiments, and subsequently confirmed in the crystalline solid state. Hence, the 18 compounds were crystallized through crystallization assays by slow solvent evaporation, followed by single-crystal X-ray diffraction (SC-XRD). The supramolecular cluster demarcation was employed to evaluate all intermolecular interactions, and all crystalline structures exhibited robust hydrogen bonds, with an average energy of approximately -21.64 kcal mol-1 (∼19% of the total stabilization energy of the supramolecular clusters), irrespective of the substituents at positions 4, 5, or 6 of the pyrimidinone core. To elucidate the nature of these hydrogen bonds, an analysis based on the quantum theory of atoms in molecules (QTAIM) revealed that the predominant intermolecular interactions are N-H···O (average of -16.55 kcal mol-1) and C-H···O (average of -6.48 kcal mol-1). Through proposing crystallization mechanisms based on molecular stabilization energy data and contact areas between molecules and employing the supramolecular cluster and retrocrystallization concepts, it was determined that altering the halogen (F/Cl) at position 4 of the pyrimidinone nucleus modifies the crystallization mechanism pathway. Notably, the hydrogen bonds present in the initial proposed steps were confirmed by 1H NMR experiments using concentration-dependent techniques.

Expedient Decagram-Scale Synthesis of Robust Organic Cages That Bind Sulfate Strongly and Selectively in Water

Selective anion recognition remains a key challenge in supramolecular chemistry: only a very small number of systems that can function in water are known, and these nearly always preferentially bind hydrophobic anions. In this work, we report three robust hexa-cationic cages that can be prepared on scales up to 14 g in two simple and high-yielding steps from commercially available materials. One of these cages displays unusually strong sulfate binding in water (Ka = 12,000 M-1), and demonstrates high selectivity for this anion over H2PO4-/HPO42- in DMSO/buffer mixtures. These results demonstrate that relatively large, three-dimensional supramolecular hosts can be prepared in high yields and on large scales, and can be highly potent receptors.

Harnessing conformational drivers in drug design

This review article explores the pivotal role of conformational drivers in the discovery of drug-like molecules and illustrates their significance through real-life examples. Understanding molecular conformation is paramount to drug hunting as it can impact on- and off-target potency, metabolism, permeability, and solubility. Each conformational driver or effector is described and exemplified in a separate section. The final section is dedicated to NMR spectroscopy and illustrates its utility as an essential tool for conformational design.

Going beyond the Ordered Bulk: A Perspective on the Use of the Cambridge Structural Database for Predictive Materials Design

When Olga Kennard founded the Cambridge Crystallographic Data Centre in 1965, the Cambridge Structural Database was a pioneering attempt to collect scientific data in a standard format. Since then, it has evolved into an indispensable resource in contemporary molecular materials science, with over 1.25 million structures and comprehensive software tools for searching, visualizing and analyzing the data. In this perspective, we discuss the use of the CSD and CCDC tools to address the multiscale challenge of predictive materials design. We provide an overview of the core capabilities of the CSD and CCDC software and demonstrate their application to a range of materials design problems with recent case studies drawn from topical research areas, focusing in particular on the use of data mining and machine learning techniques. We also identify several challenges that can be addressed with existing capabilities or through new capabilities with varying levels of development effort.

Sublimation of pyridine derivatives: fundamental aspects and application for two-component crystal screening

The saturated vapour pressures of five heterocyclic compounds containing the pyridine fragment, namely, three isomers of aminopyridine (2-aminopyridine (2AmPy), 3-aminopyridine (3AmPy), and 4-aminopyridine (4AmPy)); 3-hydroxypyridine (3OHPy) and 2-(1H-imidazol-2-yl)pyridine (ImPy), were measured at appropriate temperature intervals using a transpiration (inert gas flow) method. The standard molar enthalpies, entropies, and Gibbs energies of sublimation for all the studied substances were determined. Among the compounds studied, the largest value of ΔH298sub was observed for ImPy. The influence of substitution and the effects of hydrogen bonds in the crystal lattices on sublimation parameters are discussed herein. The reliable dependences relating ΔG298sub to Tfus and ΔH298sub to ΔG298sub were plotted. A comparative analysis of several calculation schemes for the estimation of sublimation enthalpy and Gibbs free energy was carried out. Thermodynamic parameters obtained in this study were applied for the evaluation of cocrystallisation thermodynamic functions for two-component crystals (virtual screening) on the basis of the studied substituted pyridines.

Exfoliatable Layered 2D Honeycomb Crystals of Host-guest Complexes Networked by CH-π Hydrogen Bonds

Studies of graphene show that robust chemical bonds such as covalent bonds with trigonal-planar atoms afford layered atomic 2D crystals possessing unique properties. Although layered molecular crystals are of interest to diversify elements and structures of 2D materials, the structural diversity of molecules as well as weak intermolecular interactions inevitably makes the design to be one-off and individual. We herein report a versatile method to assemble layered molecular crystals. By developing a D3-symmetry host at vertices to form a honeycomb layer, a diverse range of layered 2D host-guest crystals were obtained. Substituents on the host, elements/structures of the guest, the stereochemistry of the host and types of intercalants were diversified, which should allow for 6×32×3×2 combinations for structural diversification.

Halogen Bonding in Perfluoroalkyl Adsorption

Adsorption is a promising technology to remove perfluoroalkyl substances (PFAS), including perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA), from contaminated water. Although a large number of materials have been evaluated for PFAS adsorption, guidelines that can facilitate the rational design and selection of adsorbents have not been established due to the lack of a mechanistic understanding on the molecular level. Using a novel interpretation of the Freundlich isotherm, this study identifies halogen bonding as the main mechanism controlling perfluoroalkyl adsorption by using a materiomic approach that compares the electrostatic polarities of a variety of carbon, polymer, and mineral-based materials reported in the literature. Comparisons show that both PFOS and PFOA are favorably adsorbed by materials containing high densities of π electrons, lone electron pairs, and negative charges, consistent with the formation of halogen bonding between the positive σ-hole of fluorine as a Lewis acid and a nucleophilic solid as a Lewis base. The identification of this previously unappreciated noncovalent bonding mechanism offers fresh insight into the search of suitable materials for perfluoroalkyl adsorption.

The Peculiar H-Bonding Network of 4-Methylcatechol: A Coupled Diffraction and In Silico Study

The crystal structure of 4-methylcatechol (4MEC) has, to date, never been solved, despite its very simple chemical formula C7O2H8 and the many possible applications envisaged for this molecule. In this work, this gap is filled and the structure of 4MEC is obtained by combining X-ray powder diffraction and first principle calculations to carefully locate hydrogen atoms. Two molecules are present in the asymmetric unit. Hirshfeld analysis confirmed the reliability of the solved structure, since the two molecules show rather different environments and H-bond interactions of different directionality and strength. The packing is characterised by a peculiar hydrogen bond network with hydroxyl nests formed by two adjacent octagonal frameworks. It is noteworthy that the observed short contacts suggest strong inter-molecular interactions, further confirmed by strong inter-crystalline aggregation observed by microscopic images, indicating the growth, in many crystallization attempts, of single aggregates taller than half a centimetre and, often, with spherical shapes. These peculiarities are induced by the presence of methyl group in 4MEC, since the parent compound catechol, despite its chemical similarity, shows a standard layered packing alternating hydrophobic and polar layers. Finally, the complexity and peculiarity of the packing and crystal growth features explain why a single crystal could not be obtained for a standard structural analysis.

Linking solid-state phenomena via energy differences in `archetype crystal structures'

Categorization underlies understanding. Conceptualizing solid-state structures of organic molecules with `archetype crystal structures' bridges established categories of disorder, polymorphism and solid solutions and is herein extended to special position and high-Z' structures. The concept was developed in the context of disorder modelling [Dittrich, B. (2021). IUCrJ, 8, 305-318] and relies on adding quantum chemical energy differences between disorder components to other criteria as an explanation as to why disorder - and disappearing disorder - occurs in an average structure. Part of the concept is that disorder, as probed by diffraction, affects entire molecules, rather than just the parts of a molecule with differing conformations, and the finding that an R·T energy difference between disorder archetypes is usually not exceeded. An illustrative example combining disorder and special positions is the crystal structure of oestradiol hemihydrate analysed here, where its space-group/subgroup relationship is required to explain its disorder of hydrogen-bonded hydrogen atoms. In addition, we show how high-Z' structures can also be analysed energetically and understood via archetypes: high-Z' structures occur when an energy gain from combining different rather than overall alike conformations in a crystal significantly exceeds R·T, and this finding is discussed in the context of earlier explanations in the literature. Twinning is not related to archetype structures since it involves macroscopic domains of the same crystal structure. Archetype crystal structures are distinguished from crystal structure prediction trial structures in that an experimental reference structure is required for them. Categorization into archetype structures also has practical relevance, leading to a new practice of disorder modelling in experimental least-squares refinement alluded to in the above-mentioned publication.

Synthesis and crystal structure of an iron triazole complex resulting from the unexpected ligand cleavage of a triazolium carbene precursor

Typically reactions of N-heterocyclic carbenes with transition metals are straightforward and require a carbene salt, a base strong enough to deprotonate such a salt and a metal. Yet when carbene precursors are in the form of triazolium salts, reaction may not proceed as easily as expected. In our work, we intended to obtain a triazolylidene complex of iron(II) chloride, but due to the presence of small amounts of water in the tetrahydrofuran solvent used, bis(acetonitrile)tetrakis(1-benzyl-1H-1,2,4-triazole-κN4)iron(II) μ-oxido-bis[trichloridoferrate(III)] acetonitrile disolvate, [Fe(C9H9N3)4(CH3CN)2][Fe2Cl6O]·2CH3CN - an interesting anion with a linear geometry of the O atom - was formed instead of the iron carbene complex. Reaction proceeded via cleavage of the alkyl N-substituent of the triazolium salt. The formation of the product was confirmed by X-ray crystallography. The crystal structure and possible reaction pathways are discussed.

Benzoyl Valine Quasiracemates: Pairing CF3 Quasienantiomers with H to t-Butyl

Understanding the interplay of structural features responsible for molecular assembly is essential for molecular crystal engineering. When assembling molecules with encoded motifs, first choice supramolecular strategies almost always include robust directional nonbonded contacts. Quasiracemic materials, considered near racemates since cocrystallization occurs with chemically unique components, lack a molecular framework or functional group restrictions, highlighting the importance of molecular shape to molecular assembly. Recently, our group reported quasiracemates derived from benzoyl leucine/phenylalanine derivatives with two points of chemical difference. In this study, we modified the chemical framework with valine and increased the scope of the work by imposing a larger variance in the side chain substituents. Pairing a CF3 component with quasienantiomers that differ iteratively from hydrogen to t-butyl offers an important view into the supramolecular landscape of these materials. Single-crystal X-ray crystallography and lattice energy assessments, coupled with conformational and crystal structure similarity searches, show an elevated degree of isomorphism for many of the targeted 17 racemates and quasiracemates. These benzoyl amino acid molecular architectures create extended hydrogen-bond patterns in the crystal that provide enhanced opportunities to study the shape space and molecular recognition profiles for a diverse family of quasienantiomeric components.

Pyramidalization of the Carboxamide sp2-Center in Peptide Structures

A CSD search in the Cambridge Crystallographic Database for the substructure N-CαH-C'(═O)-N gave 24,180 peptide structures for analysis of the pyramidalization of the sp2-hybridized carboxamide group C'(═O)NCα, which had not been investigated before. The dependence of the pyramidalization θ = O-N-C'-Cα on the rotation angle ψ = O═C'-Cα-N about bond C'-Cα resulted in a curve with three maxima, three minima, and six zero-crossings. Surprisingly, the ψ/θ analysis of the individual amino acid building blocks showed that all of them exhibited similar curves, irrespective of their different R substituents. This unusual behavior is explained by a 3-fold short-range potential set up by the three covalent bonds, emanating from Cα. The tie-up of the rotation angle ψ and the pyramidalization θ in a rigid coupling is remarkable. In the 24,180 peptide structures, subjected to X-ray crystallography, there is no dynamics. For peptides in solution, the rotation/pyramidalization curve ψ/θav determines the degree of pyramidalization θ, when the rotation angle ψ runs through a full 360° circle. Density functional theory (DFT) calculations of alaninamide supported the analysis.

Organic crystal structure prediction via coupled generative adversarial networks and graph convolutional networks

Organic crystal structures exert a profound impact on the physicochemical properties and biological effects of organic compounds. Quantum mechanics (QM)-based crystal structure predictions (CSPs) have somewhat alleviated the dilemma that experimental crystal structure investigations struggle to conduct complete polymorphism studies, but the high computing cost poses a challenge to its widespread application. The present study aims to construct DeepCSP, a feasible pure machine learning framework for minute-scale rapid organic CSP. Initially, based on 177,746 data entries from the Cambridge Crystal Structure Database, a generative adversarial network was built to conditionally generate trial crystal structures under selected feature constraints for the given molecule. Simultaneously, a graph convolutional attention network was used to predict the density of stable crystal structures for the input molecule. Subsequently, the distances between the predicted density and the definition-based calculated density would be considered to be the crystal structure screening and ranking basis, and finally, the density-based crystal structure ranking would be output. Two such distinct algorithms, performing the generation and ranking functionalities, respectively, collectively constitute the DeepCSP, which has demonstrated compelling performance in marketed drug validations, achieving an accuracy rate exceeding 80% and a hit rate surpassing 85%. Inspiringly, the computing speed of the pure machine learning methodology demonstrates the potential of artificial intelligence in advancing CSP research.

Stereoselective Single Step Cyclization to Give Belt-Functionalized Pillar[6]arenes

Two macrocycles were synthesized through cyclization reactions of secondary benzylic alcohols, giving pillar[6]arenes with a methyl substituent at each belt position. These macrocycles form stereoselectively with only the rtctct isomer with alternating up and down orientations of the belt methyl groups definitively identified. Isolated yields were modest (7 and 9%), but the macrocycles are prepared in a single step from either a commercially available alcohol or a very readily prepared precursor. X-ray crystal structures of the macrocycles indicate they have a capsule-like structure, which is far from the conventional pillar shape. Density functional theory calculations reveal that the energy barrier required to obtain the pillar conformation is significantly higher for these belt-functionalized macrocycles than for conventional belt-unfunctionalized pillar[6]arenes.

A Cationic Catechol Derivative Binds Anions in Competitive Aqueous Media

A simple dihydroxy isoquinolinium molecule (3+ ) was prepared by a modification of a literature procedure. Interestingly, during optimisation of the synthesis a small amount of the natural product pseudopalmatine was isolated, and characterised for the first time by X-ray crystallography. Compound 3+ contains a catechol motif and positive charge on the same scaffold and was found to be a potent anion receptor, binding sulfate strongly in 8 : 2 d6 -acetone:D2 O and 7 : 3 d6 -acetone:D2 O (Ka >104 and 2,100 M-1 , respectively). Unsurprisingly, chloride binding was much weaker, even in the less polar solvent mixture 9 : 1 d6 -acetone:D2 O. The sulfate binding is remarkably strong for such a simple molecule, however anion binding studies were complicated by the tendency of the molecule to react with BPh4 - or BF4 - species during anion metathesis reactions. This gave two unusual zwitterions containing tetrahedral boronate centres, which were both characterised by X-ray crystallography.

Crystal Clear: Decoding Isocyanide Intermolecular Interactions through Crystallography

The isocyanide group is the chameleon among the functional groups in organic chemistry. Unlike other multiatom functional groups, where the electrophilic and nucleophilic moieties are typically separated, isocyanides combine both functionalities in the terminal carbon. This unique feature can be rationalized using the frontier orbital concept and has significant implications for its intermolecular interactions and the reactivity of the functional group. In this study, we perform a Cambridge Crystallographic Database-supported analysis of isocyanide intramolecular interactions to investigate the intramolecular interactions of isocyanides in the solid state, excluding isocyanide-metal complexes. We discuss examples of different interaction classes, including the isocyanide as a hydrogen bond acceptor (RNC···HX), halogen bonding (RNC···X), and interactions involving the isocyanide and carbon atoms (RNC···C). The latter interaction serves as an intriguing illustration of a Bürgi-Dunitz trajectory and represents a crucial experimental detail in the well-known multicomponent reactions such as the Ugi- and Passerini-type mechanisms. Understanding the spectrum of intramolecular interactions that isocyanides can undergo holds significant implications in fields such as medicinal chemistry, materials science, and asymmetric catalysis.

Graph Neural Networks with Multi-features for Predicting Cocrystals using APIs and Coformers Interactions

INTRODUCTION

Active pharmaceutical ingredients (APIs) have gained direct pharmaceutical interest, along with their in vitro properties, and thus utilized as auxiliary solid dosage forms upon FDA guidance and approval on pharmaceutical cocrystals when reacting with coformers, as a potential and attractive route for drug substance development.

METHODS

However, screening and selecting suitable and appropriate coformers that may potentially react with APIs to successfully form cocrystals is a time-consuming, inefficient, economically expensive, and labour-intensive task. In this study, we implemented GNNs to predict the formation of cocrystals using our introduced API-coformers relational graph data. We further compared our work with previous studies that implemented descriptor-based models (e.g., random forest, support vector machine, extreme gradient boosting, and artificial neural networks).

RESULTS

All built graph-based models show compelling performance accuracies (i.e., 91.36, 94.60 and 95. 95% for GCN, GraphSAGE, and RGCN respectively). RGCN demonstrated effectiveness and prevailed among the built graph-based models due to its capability to capture intricate and learn nuanced relationships between entities such as non-ionic and non-covalent interactions or link information between APIs and coformers which are crucial for accurate predictions and representations.

CONCLUSION

These capabilities allows the model to adeptly learn the topological structure inherent in the graph data.

ReaLigands: A Ligand Library Cultivated from Experiment and Intended for Molecular Computational Catalyst Design

Computational catalyst design requires identification of a metal and ligand that together result in the desired reaction reactivity and/or selectivity. A major impediment to translating computational designs to experiments is evaluating ligands that are likely to be synthesized. Here, we provide a solution to this impediment with our ReaLigands library that contains >30,000 monodentate, bidentate (didentate), tridentate, and larger ligands cultivated by dismantling experimentally reported crystal structures. Individual ligands from mononuclear crystal structures were identified using a modified depth-first search algorithm and charge was assigned using a machine learning model based on quantum-chemical calculated features. In the library, ligands are sorted based on direct ligand-to-metal atomic connections and on denticity. Representative principal component analysis (PCA) and uniform manifold approximation and projection (UMAP) analyses were used to analyze several tridentate ligand categories, which revealed both the diversity of ligands and connections between ligand categories. We also demonstrated the utility of this library by implementing it with our building and optimization tools, which resulted in the very rapid generation of barriers for 750 bidentate ligands for Rh-hydride ethylene migratory insertion.

Systematic Comparison of Experimental Crystallographic Geometries and Gas-Phase Computed Conformers for Torsion Preferences

We performed exhaustive torsion sampling on more than 3 million compounds using the GFN2-xTB method and performed a comparison of experimental crystallographic and gas-phase conformers. Many conformer sampling methods derive torsional angle distributions from experimental crystallographic data, limiting the torsion preferences to molecules that must be stable, synthetically accessible, and able to be crystallized. In this work, we evaluate the differences in torsional preferences of experimental crystallographic geometries and gas-phase computed conformers from a broad selection of compounds to determine whether torsional angle distributions obtained from semiempirical methods are suitable priors for conformer sampling. We find that differences in torsion preferences can be mostly attributed to a lack of available experimental crystallographic data with small deviations derived from gas-phase geometry differences. GFN2 demonstrates the ability to provide accurate and reliable torsional preferences that can provide a basis for new methods free from the limitations of experimental data collection. We provide Gaussian-based fits and sampling distributions suitable for torsion sampling and propose an alternative to the widely used "experimental torsion and knowledge distance geometry" (ETKDG) method using quantum torsion-derived distance geometry (QTDG) methods.

Metronidazole Cocrystal Polymorphs with Gallic and Gentisic Acid Accessed through Slurry, Atomization Techniques, and Thermal Methods

In this study, key features of metronidazole (MNZ) cocrystal polymorphs with gallic acid (GAL) and gentisic acid (GNT) were elucidated. Solvent-mediated phase transformation experiments in 30 solvents with varying properties were employed to control the polymorphic behavior of the MNZ cocrystal with GAL. Solvents with relative polarity (RP) values above 0.35 led to cocrystal I°, the thermodynamically stable form. Conversely, solvents with RP values below 0.35 produced cocrystal II, which was found to be only 0.3 kJ mol-1 less stable in enthalpy. The feasibility of electrospraying, including solvent properties and process conditions required, and spray drying techniques to control cocrystal polymorphism was also investigated, and these techniques were found to facilitate exclusive formation of the metastable MNZ-GAL cocrystal II. Additionally, the screening approach resulted in a new, high-temperature polymorph I of the MNZ-GNT cocrystal system, which is enantiotropically related to the already known form II°. The intermolecular energy calculations, as well as the 2D similarity between the MNZ-GAL polymorphs and the 3D similarity between MNZ-GNT polymorphs, rationalized the observed transition behaviors. Furthermore, the evaluation of virtual cocrystal screening techniques identified molecular electrostatic potential calculations as a supportive tool for coformer selection.

Carbon-Bonding Metal Catalysis (CBMC): A Supramolecular Complex Directs Structural-Isomer Selection in Gold-Catalyzed Reactions

Carbon is a primary element to constitute organic molecules, while metal catalysis is a basic tool in organic synthesis. The establishment of a link between the ubiquitous carbon bonding and metal catalysis is thus a fundamentally important problem. However, there is yet no experimental example to introduce the role of carbon bonding in a metal catalysis process. Herein, we merged the topics of carbon bonding and metal catalysis together and demonstrated that a supramolecular carbon-bonding metal complex can not only give rise to catalytic activity but, more remarkably, direct structural-isomer selection events in gold-catalyzed reactions. The experimental results unveil the fact that the imposing of weak carbon-bonding interactions on a gold complex can alter the carbene as well as the Lewis acid property of these catalysts. These results illustrate a non-negligible role of weak carbon-bonding interactions in the modulation of metal catalysis. As such, carbon-bonding metal catalysis is suggested to be used as a routine tool not only in the development of reactions but more frequently in analyzing reaction processes in metal catalysis.

A structural comparison of salt forms of dopamine with the structures of other phenylethylamines

The structures of four salt forms of dopamine are reported. These are dopamine [2-(3,4-dihydroxyphenyl)ethan-1-aminium] benzoate, C8H12NO2+·C7H5O2-, I, dopamine 4-nitrobenzoate, C8H12NO2+·C7H4NO4-, II, dopamine ethanedisulfonate, 2C8H12NO2+·C2H4O6S22-, III, and dopamine 4-hydroxybenzenesulfonate monohydrate, C8H12NO2+·C6H5O4S-·H2O, IV. In all four structures, the dopamine cation adopts an extended conformation. Intermolecular interaction motifs that are common in the salt forms of tyramine can be found in related dopamine structures, but hydrogen bonding in the dopamine structures appear to be more variable and less predictable than for tyramine. Packing analysis discovered three dopamine-containing groups of structures that can be described as isostructural with regards to the cation positions. Two of these groups contain both dopamine and tyramine species, and one of these is also highly variable in other ways too, containing anhydrous and hydrated forms, different anion types and ionized and neutral phenylethylamine species. As such, the group illustrates that packing behaviour can be robust and similar even where intermolecular interactions such as hydrogen bonds are very different.

Hydrogen atoms in supramolecular chemistry: a structural perspective. Where are they, and why does it matter?

Hydrogen bonding interactions are ubiquitous across the biochemical and chemical sciences, and are of particular interest to supramolecular chemists. They have been used to assemble hydrogen bonded polymers, cages and frameworks, and are the functional motif in many host-guest systems. Single crystal X-ray diffraction studies are often used as a key support for proposed structures, although this presents challenges as hydrogen atoms interact only weakly with X-rays. In this Tutorial Review, we discuss the information that can be gleaned about hydrogen bonding interactions through crystallographic experiments, key limitations of the data, and emerging techniques to overcome these limitations.

Thiolated α-cyclodextrin: The likely smallest drug carrier providing enhanced cellular uptake and endosomal escape

This study aimed to evaluate the effect of thiolated α-cyclodextrin (α-CD-SH) on the cellular uptake of its payload. For this purpose, α-CD was thiolated using phosphorous pentasulfide. Thiolated α-CD was characterized by FT-IR and ¹H NMR spectroscopy, differential scanning calorimetry (DSC), and powder X-ray diffractometry (PXRD). Cytotoxicity of α-CD-SH was evaluated on Caco-2, HEK 293, and MC3T3 cells. Dilauryl fluorescein (DLF) and coumarin-6 (Cou) serving as surrogates for a pharmaceutical payload were incorporated in α-CD-SH, and cellular uptake was analyzed by flow cytometry and confocal microscopy. Endosomal escape was investigated by confocal microscopy and hemolysis assay. Results showed no cytotoxic effect within 3 h, while dose-dependent cytotoxicity was observed within 24 h. The cellular uptake of DLF and Cou was up to 20- and 11-fold enhanced by α-CD-SH compared to native α-CD, respectively. Furthermore, α-CD-SH provided an endosomal escape. According to these results, α-CD-SH is a promising carrier to shuttle drugs into the cytoplasm of target cells.

Efficient, Formal, Material, and Final Causes in Biology and Technology

This paper considers how a classification of causal effects as comprising efficient, formal, material, and final causation can provide a useful understanding of how emergence takes place in biology and technology, with formal, material, and final causation all including cases of downward causation; they each occur in both synchronic and diachronic forms. Taken together, they underlie why all emergent levels in the hierarchy of emergence have causal powers (which is Noble's principle of biological relativity) and so why causal closure only occurs when the upwards and downwards interactions between all emergent levels are taken into account, contra to claims that some underlying physics level is by itself causality complete. A key feature is that stochasticity at the molecular level plays an important role in enabling agency to emerge, underlying the possibility of final causation occurring in these contexts.

Pharmaceutical Cocrystal Discovery via 3D-SMINBR: A New Network Recommendation Tool Augmented by 3D Molecular Conformations

Cocrystals have significant potential in various fields such as chemistry, material, and medicine. For instance, pharmaceutical cocrystals have the ability to address issues associated with physicochemical and biopharmaceutical properties. However, it can be challenging to find proper coformers to form cocrystals with drugs of interest. Herein, a new in silico tool called 3D substructure-molecular-interaction network-based recommendation (3D-SMINBR) has been developed to address this problem. This tool first integrated 3D molecular conformations with a weighted network-based recommendation model to prioritize potential coformers for target drugs. In cross-validation, the performance of 3D-SMINBR surpassed the 2D substructure-based predictive model SMINBR in our previous study. Additionally, the generalization capability of 3D-SMINBR was confirmed by testing on unseen cocrystal data. The practicality of this tool was further demonstrated by case studies on cocrystal screening of armillarisin A (Arm) and isoimperatorin (iIM). The obtained Arm-piperazine and iIM-salicylamide cocrystals present improved solubility and dissolution rate compared to their parent drugs. Overall, 3D-SMINBR augmented by 3D molecular conformations would be a useful network-based tool for cocrystal discovery. A free web server for 3D-SMINBR can be freely accessed at http://lmmd.ecust.edu.cn/netcorecsys/.

Polymorphism of Carbamazepine Pharmaceutical Cocrystal: Structural Analysis and Solubility Performance

Polymorphism is a common phenomenon among single- and multicomponent molecular crystals that has a significant impact on the contemporary drug development process. A new polymorphic form of the drug carbamazepine (CBZ) cocrystal with methylparaben (MePRB) in a 1:1 molar ratio as well as the drug's channel-like cocrystal containing highly disordered coformer molecules have been obtained and characterized in this work using various analytical methods, including thermal analysis, Raman spectroscopy, and single-crystal and high-resolution synchrotron powder X-ray diffraction. Structural analysis of the solid forms revealed a close resemblance between novel form II and previously reported form I of the [CBZ + MePRB] (1:1) cocrystal in terms of hydrogen bond networks and overall packing arrangements. The channel-like cocrystal was found to belong to a distinct family of isostructural CBZ cocrystals with coformers of similar size and shape. Form I and form II of the 1:1 cocrystal appeared to be related by a monotropic relationship, with form II being proven to be the thermodynamically more stable phase. The dissolution performance of both polymorphs in aqueous media was significantly enhanced when compared with parent CBZ. However, considering the superior thermodynamic stability and consistent dissolution profile, the discovered form II of the [CBZ + MePRB] (1:1) cocrystal seems a more promising and reliable solid form for further pharmaceutical development.

Halogen bonds, chalcogen bonds, pnictogen bonds, tetrel bonds and other σ-hole interactions: a snapshot of current progress

We report here on the status of research on halogen bonds and other σ-hole interactions involving p-block elements in Lewis acidic roles, such as chalcogen bonds, pnictogen bonds and tetrel bonds. A brief overview of the available literature in this area is provided via a survey of the many review articles that address this field. Our focus has been to collect together most review articles published since 2013 to provide an easy entry into the extensive literature in this area. A snapshot of current research in the area is provided by an introduction to the virtual special issue compiled in this journal, comprising 11 articles and entitled `Halogen, chalcogen, pnictogen and tetrel bonds: structural chemistry and beyond.'

Modulating Thermal Properties of Polymers through Crystal Engineering

Crystal engineering has exclusively focused on the development of advanced materials based on small organic molecules. We now demonstrate how the cocrystallization of a polymer yields a material with significantly enhanced thermal stability but equivalent mechanical flexibility. Isomorphous replacement of one of the cocrystal components enables the formation of solid solutions with melting points that can be readily fine-tuned over a usefully wide temperature range. The results of this study credibly extend the scope of crystal engineering and cocrystallization from small molecules to polymers.

Supramolecular Organization in Salts of Riluzole with Dihydroxybenzoic Acids—The Key Role of the Mutual Arrangement of OH Groups

Intermolecular interactions, in particular hydrogen bonds, play a key role in crystal engineering. The ability to form hydrogen bonds of various types and strengths causes competition between supramolecular synthons in pharmaceutical multicomponent crystals. In this work, we investigate the influence of positional isomerism on the packing arrangements and the network of hydrogen bonds in multicomponent crystals of the drug riluzole with hydroxyl derivatives of salicylic acid. The supramolecular organization of the riluzole salt containing 2,6-dihydroxybenzoic acid differs from that of the solid forms with 2,4- and 2,5-dihydroxybenzoic acids. Because the second OH group is not at position 6 in the latter crystals, intermolecular charge-assisted hydrogen bonds are formed. According to periodic DFT calculations, the enthalpy of these H-bonds exceeds 30 kJ·mol−1. The positional isomerism appears to have little effect on the enthalpy of the primary supramolecular synthon (65–70 kJ·mol−1), but it does result in the formation of a two-dimensional network of hydrogen bonds and an increase in the overall lattice energy. According to the results of the present study, 2,6-dihydroxybenzoic acid can be treated as a promising counterion for the design of pharmaceutical multicomponent crystals.

Flexible Aliphatic Diammonioacetates as Zwitterionic Ligands in UO22+ Complexes: Diverse Topologies and Interpenetrated Structures

N,N,N',N'-Tetramethylethane-1,2-diammonioacetate (L1) and N,N,N',N'-tetramethylpropane-1,3-diammonioacetate (L2) are two flexible zwitterionic dicarboxylates which have been used as ligands for the uranyl ion, 12 complexes having been obtained from their coupling to diverse anions, mostly anionic polycarboxylates, or oxo, hydroxo and chlorido donors. The protonated zwitterion is a simple counterion in [H₂L1][UO₂(2,6-pydc)₂] (1), where 2,6-pydc^2- is 2,6-pyridinedicarboxylate, but it is deprotonated and coordinated in all the other complexes. [(UO₂)₂(L2)(2,4-pydcH)₄] (2), where 2,4-pydc^2- is 2,4-pyridinedicarboxylate, is a discrete, binuclear complex due to the terminal nature of the partially deprotonated anionic ligands. [(UO₂)₂(L1)(ipht)₂]·4H₂O (3) and [(UO₂)₂(L1)(pda)₂] (4), where ipht^2- and pda^2- are isophthalate and 1,4-phenylenediacetate, are monoperiodic coordination polymers in which central L1 bridges connect two lateral strands. Oxalate anions (ox^2-) generated in situ give [(UO₂)₂(L1)(ox)₂] (5) a diperiodic network with the hcb topology. [(UO₂)₂(L2)(ipht)₂]·H₂O (6) differs from 3 in being a diperiodic network with the V₂O₅ topological type. [(UO₂)₂(L1)(2,5-pydc)₂]·4H₂O (7), where 2,5-pydc^2- is 2,5-pyridinedicarboxylate, is a hcb network with a square-wave profile, while [(UO₂)₂(L1)(dnhpa)₂] (8), where dnhpa^2- is 3,5-dinitro-2-hydroxyphenoxyacetate, formed in situ from 1,2-phenylenedioxydiacetic acid, has the same topology but a strongly corrugated shape leading to interdigitation of layers. (2R,3R,4S,5S)-Tetrahydrofurantetracarboxylic acid (thftcH₄) is only partially deprotonated in [(UO₂)₃(L1)(thftcH)₂(H₂O)] (9), which crystallizes as a diperiodic polymer with the fes topology. [(UO₂)₂Cl₂(L1)₃][(UO₂Cl₃)₂(L1)] (10) is an ionic compound in which discrete, binuclear anions cross the cells of the cationic hcb network. 2,5-Thiophenediacetate (tdc^2-) is peculiar in promoting self-sorting of the ligands in the ionic complex [(UO₂)₅(L1)₇(tdc)(H₂O)][(UO₂)₂(tdc)₃]₄·CH₃CN·12H₂O (11), which is the first example of heterointerpenetration in uranyl chemistry, involving a triperiodic, cationic framework and diperiodic, anionic hcb networks. Finally, [(UO₂)₇(O)₃(OH)_4.3Cl_2.7(L2)₂]Cl·7H₂O (12) crystallizes as a 2-fold interpenetrated, triperiodic framework in which chlorouranate undulating monoperiodic subunits are bridged by the L2 ligands. Complexes 1, 2, 3, and 7 are emissive with photoluminescence quantum yields in the range of 8-24%, and their solid-state emission spectra show the usual dependence on number and nature of donor atoms.

Virtual Screening, Structural Analysis, and Formation Thermodynamics of Carbamazepine Cocrystals

In this study, the existing set of carbamazepine (CBZ) cocrystals was extended through the successful combination of the drug with the positional isomers of acetamidobenzoic acid. The structural and energetic features of the CBZ cocrystals with 3- and 4-acetamidobenzoic acids were elucidated via single-crystal X-ray diffraction followed by QTAIMC analysis. The ability of three fundamentally different virtual screening methods to predict the correct cocrystallization outcome for CBZ was assessed based on the new experimental results obtained in this study and data available in the literature. It was found that the hydrogen bond propensity model performed the worst in distinguishing positive and negative results of CBZ cocrystallization experiments with 87 coformers, attaining an accuracy value lower than random guessing. The method that utilizes molecular electrostatic potential maps and the machine learning approach named CCGNet exhibited comparable results in terms of prediction metrics, albeit the latter resulted in superior specificity and overall accuracy while requiring no time-consuming DFT computations. In addition, formation thermodynamic parameters for the newly obtained CBZ cocrystals with 3- and 4-acetamidobenzoic acids were evaluated using temperature dependences of the cocrystallization Gibbs energy. The cocrystallization reactions between CBZ and the selected coformers were found to be enthalpy-driven, with entropy terms being statistically different from zero. The observed difference in dissolution behavior of the cocrystals in aqueous media was thought to be caused by variations in their thermodynamic stability.

Supramolecular chemistry of two new bis(1,2,3-triazolyl)pyridine macrocycles: metal complexation, self-assembly and anion binding

Two new macrocycles containing the bis(1,2,3-triazolyl)pyridine (btp) motif were prepared in high yields from a btp diazide precursor (1). Solution ¹H NMR studies show that this diazide undergoes self-assembly with divalent transition metal ions to form ML₂ complexes with pendant azide groups, apparently suitable for conversion into metal-templated catenanes; however attempts to form these catenanes were unsuccessful. Instead a new macrocycle containing two btp motifs was prepared, which forms a nanotube structure in the solid state. Reduction of the azide groups to amines followed by amide bond formation was used to convert 1 into macrocycle 8 containing btp and isophthalamide functionalities. This macrocycle binds halide and oxalate anions in acetonitrile solely through the isophthalamide motif, and binds aromatic dicarboxylates very strongly through both the isophthalamide amide donors and the btp triazole donors. The macrocycle was complexed with Pd(II) and the resulting complexes were shown to bind strongly to halide anions. The solid state structures of [Pd·8·X]BF₄ (X = Cl^-, Br^-, I^-) were investigated by X-ray crystallography, which showed that [Pd·8·Br] forms an unusual "chain of dimers" structure assembled by metal complexation, N-H⋯Br^- hydrogen bonding and short Pd⋯Pd contacts.

Ab Initio Study of H-Bond Dynamics in Three-Component Crystals Comprising (DABCOH+ )n Polycationic Chains

In this work, the dynamic character of hydrogen-bond (H-bond) networks in two three-component crystals comprising polycationic chains was described. The first studied system was 1,4-diazabicyclo[2.2.2]octan-1-ium (DABCOH⁺ ) sulfamate monohydrate, known for its large negative linear compressibility. The second analyzed material was the newly obtained polar salt co-crystal: 1,4-diazabicyclo[2.2.2]octan-1-ium sulfamate urea. X-ray diffraction measurements enabled us to study the H-bond systems in both crystals using the graph set analysis. Obtained structures served as the initial models for Born-Oppenheimer molecular dynamics computations. A detailed study of intermolecular interactions and power spectra was conducted. The analysis of time and space correlations between the changes in H-bonds enabled the detection of proton transfer occurring in both systems at 300 K. Further study of those dynamic phenomena was done using the Energy Decomposition Analysis for selected trajectory fragments. Our work should improve the understanding of dielectric and ferroelectric properties of hybrid organic-inorganic materials.

The HEALED SBU Library of Chemically Realistic Building Blocks for Construction of Hypothetical Metal-Organic Frameworks

Advancements in hypothetical metal-organic framework (hMOF) databases and construction tools have resulted in a rapidly expanding chemical design space for nanoporous materials. The bulk of these hypothetical structures are constructed using structural building units (SBUs) derived from experimental MOF structures, often collected from the CoRE-MOF database. Recent investigations into the state of these deposited experimental structures' chemical accuracy identified an array of common structural errors, including omitted protons, missing counterions, and disordered structures. These structural errors propagate into the SBUs mined from experimental MOFs, culminating in inaccurate hMOF structures possessing net charges or missing atoms which were not accounted for previously. This work demonstrates how manual investigation was applied to diagnose structural errors in SBUs obtained from several popular hMOF construction tools and databases. An analysis of the prevailing errors discovered during the examination process is provided along with representative cases to aid with error detection in future studies involving SBU extraction and hMOF construction. A novel repair protocol was established and employed to generate a library of SBUs that are hand-examined and labeled with enhanced detail (HEALED). This repaired library of SBUs contains 952 inorganic SBUs and 568 organic SBUs ideally suited for the generation of hypothetical frameworks that are chemically accurate and properly charge labeled. Additionally, case studies following the effects of SBU errors on electrostatic potential-fitted charges and GCMC-simulated gas adsorption predictions are presented to highlight the significance of using chemically accurate hMOF structures exclusively in all screening efforts going forward.

Post-Synthetic Modification of a Porous Hydrocarbon Cage to Give a Discrete Co24 Organometallic Complex

A new alkyne-based hydrocarbon cage was synthesized in high overall yield using alkyne-alkyne coupling in the cage forming step. The cage is porous and displays a moderately high BET surface area (546 m2  g-1 ). The cage loses crystallinity on activation and thus is porous in its amorphous form, while very similar cages have been either non-porous, or retained crystallinity on activation. Reaction of the cage with Co2 (CO)8 results in exhaustive metalation of its 12 alkyne groups to give the Co24 (CO)72 adduct of the cage in good yield.

Zwitterionic and Anionic Polycarboxylates as Coligands in Uranyl Ion Complexes, and Their Influence on Periodicity and Topology

The three zwitterionic di- and tricarboxylate ligands 1,1'-[(2,3,5,6-tetramethylbenzene-1,4-diyl)bis(methylene)]bis(pyridin-1-ium-4-carboxylate) (pL1), 1,1'-[(2,3,5,6-tetramethylbenzene-1,4-diyl)bis(methylene)]bis(pyridin-1-ium-3-carboxylate) (mL1), and 1,1',1″-[(2,4,6-trimethylbenzene-1,3,5-triyl)tris(methylene)]tris(pyridin-1-ium-4-carboxylate) (L2) have been used as ligands to synthesize a series of 15 uranyl ion complexes involving various anionic coligands, in most cases polycarboxylates. [(UO₂)₂(pL1)₂(cbtc)(H₂O)₂]·10H₂O (1, cbtc^4- = cis,trans,cis-1,2,3,4-cyclobutanetetracarboxylate) is a discrete, dinuclear ring-shaped complex with a central cbtc^4- pillar. While [UO₂(pL1)(NO₃)₂] (2), [UO₂(pL1)(OAc)₂] (3), and [UO₂(pL1)(HCOO)₂] (4) are simple chains, [(UO₂)₂(mL1)(1,3-pda)₂] (5, 1,3-pda^2- = 1,3-phenylenediacetate) is a daisy chain and [UO₂(pL1)(pdda)]₃·10H₂O (6, pdda^2- = 1,2-phenylenedioxydiacetate) is a double-stranded, ribbon-like chain. Both [UO₂(pL1)(pht)]·5H₂O (7, pht^2- = phthalate) and [(UO₂)₃(mL1)(pht)₂(OH)₂] (8) crystallize as diperiodic networks with the sql topology, the latter involving hydroxo-bridged trinuclear nodes. [(UO₂)₂(pL1)(c/t-1,3-chdc)₂] (9, c/t-1,3-chdc^2- = cis/trans-1,3-cyclohexanedicarboxylate) and [UO₂(pL1)(t-1,4-chdc)]·1.5H₂O (10, t-1,4-chdc^2- = trans-1,4-cyclohexanedicarboxylate) are also diperiodic, with the V₂O₅ and sql topologies, respectively. Both [(UO₂)₂(mL1)(c/t-1,4-chdc)₂] (11) and [(UO₂)₂(pL1)(1,2-pda)₂] (12, 1,2-pda^2- = 1,2-phenylenediacetate) crystallize as diperiodic networks with hcb topology, and they display threefold parallel interpenetration. [HL2][(UO₂)₃(L2)(adc)₃]Br (13, adc^2- = 1,3-adamantanedicarboxylate) contains a very corrugated hcb network with two different kinds of cells, and the uncoordinated HL2⁺ molecule associates with the coordinated L2 to form a capsule containing the bromide anion. [(UO₂)₂(pL1)(kpim)₂] (14, kpim^2- = 4-ketopimelate) is a three-periodic framework with pL1 molecules pillaring fes diperiodic subunits, whereas [(UO₂)₂(L2)₂(t-1,4-chdc)](NO₃)_1.7Br_0.3·6H₂O (15), the only cationic complex in the series, is a triperiodic framework with dmc topology and t-1,4-chdc^2- anions pillaring fes diperiodic subunits. Solid-state emission spectra and photoluminescence quantum yields are reported for all complexes.

Data format standards in analytical chemistry

Research data is an essential part of research and almost every publication in chemistry. The data itself can be valuable for reuse if sustainably deposited, annotated and archived. Thus, it is important to publish data following the FAIR principles, to make it findable, accessible, interoperable and reusable not only for humans but also in machine-readable form. This also improves transparency and reproducibility of research findings and fosters analytical work with scientific data to generate new insights, being only accessible with manifold and diverse datasets. Research data requires complete and informative metadata and use of open data formats to obtain interoperable data. Generic data formats like AnIML and JCAMP-DX have been used for many applications. Special formats for some analytical methods are already accepted, like mzML for mass spectrometry or nmrML and NMReDATA for NMR spectroscopy data. Other methods still lack common standards for data. Only a joint effort of chemists, instrument and software vendors, publishers and infrastructure maintainers can make sure that the analytical data will be of value in the future. In this review, we describe existing data formats in analytical chemistry and introduce guidelines for the development and use of standardized and open data formats.

Applications of Artificial Intelligence and Machine Learning Algorithms to Crystallization

Artificial intelligence and specifically machine learning applications are nowadays used in a variety of scientific applications and cutting-edge technologies, where they have a transformative impact. Such an assembly of statistical and linear algebra methods making use of large data sets is becoming more and more integrated into chemistry and crystallization research workflows. This review aims to present, for the first time, a holistic overview of machine learning and cheminformatics applications as a novel, powerful means to accelerate the discovery of new crystal structures, predict key properties of organic crystalline materials, simulate, understand, and control the dynamics of complex crystallization process systems, as well as contribute to high throughput automation of chemical process development involving crystalline materials. We critically review the advances in these new, rapidly emerging research areas, raising awareness in issues such as the bridging of machine learning models with first-principles mechanistic models, data set size, structure, and quality, as well as the selection of appropriate descriptors. At the same time, we propose future research at the interface of applied mathematics, chemistry, and crystallography. Overall, this review aims to increase the adoption of such methods and tools by chemists and scientists across industry and academia.

Crystal Engineering of Pharmaceutical Cocrystals in the Discovery and Development of Improved Drugs

The subject of crystal engineering started in the 1970s with the study of topochemical reactions in the solid state. A broad chemical definition of crystal engineering was published in 1989, and the supramolecular synthon concept was proposed in 1995 followed by heterosynthons and their potential applications for the design of pharmaceutical cocrystals in 2004. This review traces the development of supramolecular synthons as robust and recurring hydrogen bond patterns for the design and construction of supramolecular architectures, notably, pharmaceutical cocrystals beginning in the early 2000s to the present time. The ability of a cocrystal between an active pharmaceutical ingredient (API) and a pharmaceutically acceptable coformer to systematically tune the physicochemical properties of a drug (i.e., solubility, permeability, hydration, color, compaction, tableting, bioavailability) without changing its molecular structure is the hallmark of the pharmaceutical cocrystals platform, as a bridge between drug discovery and pharmaceutical development. With the design of cocrystals via heterosynthons and prototype case studies to improve drug solubility in place (2000-2015), the period between 2015 to the present time has witnessed the launch of several salt-cocrystal drugs with improved efficacy and high bioavailability. This review on the design, synthesis, and applications of pharmaceutical cocrystals to afford improved drug products and drug substances will interest researchers in crystal engineering, supramolecular chemistry, medicinal chemistry, process development, and pharmaceutical and materials sciences. The scale-up of drug cocrystals and salts using continuous manufacturing technologies provides high-value pharmaceuticals with economic and environmental benefits.

Mining Knowledge from Crystal Structures: Oxidation States of Oxygen-Coordinated Metal Atoms in Ionic and Coordination Compounds

We propose a universal scheme for predicting the oxidation states of metal atoms in ionic and coordination compounds with a small set of structural descriptors, which include the parameters of atomic Voronoi polyhedra. The scheme has been trained and checked with more than 35,000 crystal structures containing more than 90,000 metal atoms in the oxygen environment. The accuracy of the prediction exceeded 95%; we have detected a number of wrong oxidation states and incorrect chemical compositions in the crystallographic databases using this scheme. The scheme is easily extendable to any kind of atomic environment and can be used to search for correlations between geometrical and physical properties of crystal structures.

Evaluation of Piperine as Natural Coformer for Eutectics Preparation of Drugs Used in the Treatment of Cardiovascular Diseases

Piperine (PIP) was evaluated as a natural coformer in the preparation of multicomponent organic materials for enhancing solubility and dissolution rate of the poorly water-soluble drugs: curcumin (CUR), lovastatin (LOV), and irbesartan (IBS). A screening based on liquid assisted grinding technique was performed using 1:1 drug-PIP molar ratio mixtures, followed by differential scanning calorimetry (DSC) and powder X-ray diffraction (PXRD) analyses. Three eutectic mixtures (EMs) composed of CUR-PIP, LOV-PIP, and IBS-PIP were obtained. Therefore, binary phase and Tamman's diagrams were constructed for each system to obtain the exact eutectic composition, which was 0.41:0.59, 0.29:0.71, and 0.31:0.69 for CUR-PIP, LOV-PIP, and IBS-PIP, respectively. Further, bulk materials of each system were prepared to characterize them through DSC, PXRD fully, Fourier transform infrared spectroscopy (FT-IR), and solution-state nuclear magnetic resonance (NMR) spectroscopy. In addition, the contact angle, solubility, and dissolution rate of each system were evaluated. The preserved characteristic in the PXRD patterns and FT-IR spectra of the bulk material of each system confirmed the formation of EM mixture without molecular interaction in solid-state. The formation of EM resulted in improved aqueous solubility and dissolution rate associated with the increased wettability observed by the decrease in contact angle. In addition, solution NMR analyses of CUR-PIP, LOV-PIP, and IBS-PIP suggested no significant intermolecular interactions in solution between the components of the EM. Hence, this study concludes that PIP could be an effective coformer to improve the solubility and dissolution rate of CUR, LOV, and IBS.