"Higher density does not mean higher stability" mystery of paracetamol finally unraveled

Advancing ecotoxicity assessment: Leveraging pre-trained model for bee toxicity and compound degradability prediction

The prediction of ecological toxicity plays an increasingly important role in modern society. However, the existing models often suffer from poor performance and limited predictive capabilities. In this study, we propose a novel approach for ecological toxicity assessment based on pre-trained models. By leveraging pre-training techniques and graph neural network models, we establish a highperformance predictive model. Furthermore, we incorporate a variational autoencoder to optimize the model, enabling simultaneous discrimination of toxicity to bees and molecular degradability. Additionally, despite the low similarity between the endogenous hormones in bees and the compounds in our dataset, our model confidently predicts that these hormones are non-toxic to bees, which further strengthens the credibility and accuracy of our model. We also discovered the negative correlation between the degradation and bee toxicity of compounds. In summary, this study presents an ecological toxicity assessment model with outstanding performance. The proposed model accurately predicts the toxicity of chemicals to bees and their degradability capabilities, offering valuable technical support to relevant fields.

An unexpected high-pressure stability domain for a lower density polymorph of benzophenone

For over a century, it was thought that the crystalline polymorph II of benzophenone does not possess a stable domain in the pressure-temperature phase diagram. With a combination of new experimental results and literature data, this case of crystalline dimorphism has finally been solved and it is shown that form II possesses a stable domain at high pressure and high temperature, even though its density is lower than that of form I, the stable form under ordinary pressure and temperature conditions. The phase diagram of benzophenone is a clear demonstration of the fact that to understand the phase behaviour of a chemical substance both the exchange of heat (due to the change in intermolecular interactions) and work (due to the change of volume at a given pressure) need to be taken into account.

Advances in the development of amorphous solid dispersions: The role of polymeric carriers

Amorphous solid dispersion (ASD) is one of the most effective approaches for delivering poorly soluble drugs. In ASDs, polymeric materials serve as the carriers in which the drugs are dispersed at the molecular level. To prepare the solid dispersions, there are many polymers with various physicochemical and thermochemical characteristics available for use in ASD formulations. Polymer selection is of great importance because it influences the stability, solubility and dissolution rates, manufacturing process, and bioavailability of the ASD. This review article provides a comprehensive overview of ASDs from the perspectives of physicochemical characteristics of polymers, formulation designs and preparation methods. Furthermore, considerations of safety and regulatory requirements along with the studies recommended for characterizing and evaluating polymeric carriers are briefly discussed.

Weak, Broken, but Working-Intramolecular Hydrogen Bond in 2,2'-bipyridine

From an academic and practical point of view, it is desirable to be able to assess the possibility of the proton exchange of a given molecular system just by knowing the positions of the proton acceptor and the proton donor. This study addresses the difference between intramolecular hydrogen bonds in 2,2'-bipyridinium and 1,10-phenanthrolinium. Solid-state ¹⁵N NMR and model calculations show that these hydrogen bonds are weak; their energies are 25 kJ/mol and 15 kJ/mol, respectively. Neither these hydrogen bonds nor N-H stretches can be responsible for the fast reversible proton transfer observed for 2,2'-bipyridinium in a polar solvent down to 115 K. This process must have been caused by an external force, which was a fluctuating electric field present in the solution. However, these hydrogen bonds are the grain that tips the scales precisely because they are an integral part of a large system of interactions, including both intramolecular interactions and environmental influence.

Thermal and Structural Characterization of Two Crystalline Polymorphs of Tafamidis Free Acid

Tafamidis, chemical formula C14H7Cl2NO3, is a drug used to delay disease progression in adults suffering from transthyretin amyloidosis, and is marketed worldwide under different tradenames as a free acid or in the form of its meglumine salt. The free acid (CAS no. 594839-88-0) is reported to crystallize as distinct (polymorphic) crystal forms, the thermal stability and structural features of which remained thus far undisclosed. In this paper, we present-by selectively isolating highly pure batches of Tafamidis Form 1 and Tafamidis Form 4-the full characterization of these solids, in terms of crystal structures (determined using state-of-the-art structural powder diffraction methods) and spectroscopic and thermal properties. Beyond conventional thermogravimetric and calorimetric analyses, variable-temperature X-ray diffraction was employed to measure the highly anisotropic response of these (poly)crystalline materials to thermal stimuli and enabled the determination of the linear and volumetric thermal expansion coefficients and of the related indicatrix. Both crystal phases are monoclinic and contain substantially flat and π-π stacked Tafamidis molecules, arranged as centrosymmetric dimers by strong O-H···O bonds; weaker C-H···N contacts give rise, in both polymorphs, to infinite ribbons, which guarantee the substantial stiffness of the crystals in the direction of their elongation. Complete knowledge of the structural models will foster the usage of full-pattern quantitative phase analyses of Tafamidis in drug and polymorphic mixtures, an important aspect in both the forensic and the industrial sectors.

The Relevance of Experimental Charge Density Analysis in Unraveling Noncovalent Interactions in Molecular Crystals

The work carried out by our research group over the last couple of decades in the context of quantitative crystal engineering involves the analysis of intermolecular interactions such as carbon (tetrel) bonding, pnicogen bonding, chalcogen bonding, and halogen bonding using experimental charge density methodology is reviewed. The focus is to extract electron density distribution in the intermolecular space and to obtain guidelines to evaluate the strength and directionality of such interactions towards the design of molecular crystals with desired properties. Following the early studies on halogen bonding interactions, several "sigma-hole" interaction types with similar electrostatic origins have been explored in recent times for their strength, origin, and structural consequences. These include interactions such as carbon (tetrel) bonding, pnicogen bonding, chalcogen bonding, and halogen bonding. Experimental X-ray charge density analysis has proved to be a powerful tool in unraveling the strength and electronic origin of such interactions, providing insights beyond the theoretical estimates from gas-phase molecular dimer calculations. In this mini-review, we outline some selected contributions from the X-ray charge density studies to the field of non-covalent interactions (NCIs) involving elements of the groups 14-17 of the periodic table. Quantitative insights into the nature of these interactions obtained from the experimental electron density distribution and subsequent topological analysis by the quantum theory of atoms in molecules (QTAIM) have been discussed. A few notable examples of weak interactions have been presented in terms of their experimental charge density features. These examples reveal not only the strength and beauty of X-ray charge density multipole modeling as an advanced structural chemistry tool but also its utility in providing experimental benchmarks for the theoretical studies of weak interactions in crystals.

Precipitation dominated thin films of acetaminophen fabricated by meniscus guided coating

Crystallization above the solvent boiling point facilitates the identification of a new precipitation dominant morphology during meniscus guided coating.

Long-range supramolecular synthon polymorphism: a case study of two new polymorphic cocrystals of Ph2Te2–1,4-C6F4I2

Two new polymorphic forms of Ph2Te2–1,4-C6F4I2 cocrystals feature an unusual packing of Ph2Te2 molecules, which is typical for native Ph2Se2 but not Ph2Te2. This suggests the existence the yet unknown, Ph2Se2-like polymorph of Ph2Te2.

Polymorphic forms of antiandrogenic drug nilutamide: structural and thermodynamic aspects

Attempts to obtain new cocrystals of nonsteroidal antiandrogenic drug nilutamide produced alternative polymorphic forms of the compound (Form II and Form III) and their crystal structures were elucidated by single-crystal X-ray diffraction. Apart from the cocrystallization technique, lyophilization was found to be an effective strategy for achieving polymorph control of nilutamide, which was difficult to obtain by other methods. The physicochemical properties and relative stability of the commercial Form I and newly obtained Form II were comprehensively investigated by a variety of analytical methods (thermal analysis, solution calorimetry, solubility, and sublimation), whereas for Form III, only a handful of experimental parameters were obtained due to the elusive nature of the polymorph. Form I and Form II were found to be monotropically related, with Form I being confirmed as the thermodynamically most stable solid phase. In addition, the performance of different DFT-D and semi-empirical schemes for lattice energy calculation and polymorph energy ranking was compared and analysed. Lattice energy calculations using periodic DFT at B3LYP-D3/6-31(F+)G(d,p) and PBEh-3c/def2-mSVP levels of theory were found to provide the most accurate lattice energy values for Form I against experimental data, while PIXEL and PBEh-3c/def2-mSVP were the only methods that predicted the correct order of stability of Forms I and II.

Three differently coloured polymorphs of 3,6-bis(4-chlorophenyl)-2,5-dipropyl-2,5-dihydropyrrolo[3,4-c]pyrrole-1,4-dione

In this paper, the conformational polymorphism of a chlorinated diketopyrrolopyrrole (DPP) dye having flexible substituents in a non-hydrogen-bonding system is reported. The propyl-substituted DPP derivative (PR3C) has three polymorphic forms, each showing a different colour (red, orange and yellow). All polymorphs could be obtained concomitantly under various crystallization conditions. The results of the crystal structure analysis indicate that PR3C adopts different conformations in each polymorph. The packing effect caused by the difference in the arrangement of neighbouring molecules was found to play an important role in the occurrence of the observed polymorphism. The thermodynamic stability relationship between the three polymorphs was identified by thermal analysis and indicates that the yellow polymorph is the thermally stable form. The results indicate that the yellow form and orange form are enantiotropically related, and the other polymorph is monotropically related to the others.

Bulk and Surface-Stabilized Structures of Paracetamol Revisited by Raman Confocal Microscopy

We revisit the polymorphism of paracetamol by means of a micro-Raman technique, which has proved to be a powerful tool for structure recognition. Distinct lattice phonon spectra clearly identified the pure phases. Confocality enabled us to detect phase mixing between form II and either I or III on a micrometric scale in the same crystallite. Following the most recent findings on surface-mediated structures, we also investigated spin-coated films grown on glass, gold, and polystyrene substrates, confirming the selectivity of these surfaces for the metastable form III, which shows an unprecedented stability over a time span of several months. A mechanism of its transformation to phase II, via a partially ordered intermediate state, is suggested by polarized Raman measurements.

Printing of small molecular medicines from the vapor phase

There is growing need to develop efficient methods for early-stage drug discovery, continuous manufacturing of drug delivery vehicles, and ultra-precise dosing of high potency drugs. Here we demonstrate the use of solvent-free organic vapor jet printing to deposit nanostructured films of small molecular pharmaceutical ingredients, including caffeine, paracetamol, ibuprofen, tamoxifen, BAY 11-7082 and fluorescein, with accuracy on the scale of micrograms per square centimeter, onto glass, Tegaderm, Listerine tabs, and stainless steel microneedles. The printed films exhibit similar crystallographic order and chemistry as the original powders; controlled, order-of-magnitude enhancements of dissolution rate are observed relative to powder-form particles. In vitro treatment of breast and ovarian cancer cell cultures in aqueous media by tamoxifen and BAY 11-7082 films shows similar behavior to drugs pre-dissolved in dimethyl sulfoxide. The demonstrated precise printing of medicines as films, without the use of solvents, can accelerate drug screening and enable continuous manufacturing, while enhancing dosage accuracy.

Traditional approaches used in the pharmaceutical industry are not precise or versatile enough for customized medicine formulation and manufacture. Here the authors produce a method to form coatings, with accurate dosages, as well as a means of closely controlling dissolution kinetics.

Conformational and vibrational reassessment of solid paracetamol

This work provides an answer to the urge for a more detailed and accurate knowledge of the vibrational spectrum of the widely used analgesic/antipyretic drug commonly known as paracetamol. A comprehensive spectroscopic analysis - including infrared, Raman, and inelastic neutron scattering (INS) - is combined with a computational approach which takes account for the effects of intermolecular interactions in the solid state. This allows a full reassessment of the vibrational assignments for Paracetamol, thus preventing the propagation of incorrect data analysis and misassignments already found in the literature. In particular, the vibrational modes involving the hydrogen-bonded NH and OH groups are correctly reallocated to bands shifted by up to 300cm^-1 relatively to previous assignments.

Recent progress of structural study of polymorphic pharmaceutical drugs

This review considers advances in the understanding of active pharmaceutical ingredient polymorphism since around 2010 mainly from a structural view point, with a focus on twelve model drugs. New polymorphs of most of these drugs have been identified despite that the polymorphism of these old drugs has been extensively studied so far. In addition to the conventional modifications of preparative solvents, temperatures, and pressure, more strategic structure-based methods have successfully yielded new polymorphs. The development of analytical techniques, including X-ray analyses, spectroscopy, and microscopy has facilitated the identification of unknown crystal structures and also the discovery of new polymorphs. Computational simulations have played an important role in explaining and predicting the stability order of polymorphs. Furthermore, these make significant contributions to the design of new polymorphs by considering structure and energy. The new technologies and insights discussed in this review will contribute to the control of polymorphic forms, both during manufacture and in the drug formulation.

High-pressure high-temperature phase diagram of organic crystal paracetamol

High-pressure high-temperature (HPHT) Raman spectroscopy studies have been performed on the organic crystal paracetamol in a diamond anvil cell utilizing boron-doped heating diamond anvil. Isobaric measurements were conducted at pressures up to 8.5 GPa and temperature up to 520 K in five different experiments. Solid state phase transitions from monoclinic Form I  →  orthorhombic Form II were observed at various pressures and temperatures as well as transitions from Form II  →  unknown Form IV. The melting temperature for paracetamol was observed to increase with increasing pressures to 8.5 GPa. This new data is combined with previous ambient temperature high-pressure Raman and x-ray diffraction data to create the first HPHT phase diagram of paracetamol.

Triclinic and monoclinic polymorphs of meso-(E,E)-1,1'-[1,2-bis(4-chlorophenyl)ethane-1,2-diyl]bis(phenyldiazene): the high-yield synthesis of an unexpected product, concomitant polymorphism and configurational disorder

Pyrazolidine-3,5-diones and their derivatives exhibit a wide range of biological activities. Seeking to explore the effect of combining a hydrocarbyl ring substituent, as present in sulfinpyrazone (used to treat gout), with a chlorinated aryl ring, as present in muzolimine (a diuretic), we explored the reaction between 1-phenylpyrazolidine-3,5-dione and 4-chlorobenzaldehyde under mildly basic conditions in the expectation of producing the simple condensation product 4-(4-chlorobenzylidene)-1-phenylpyrazolidine-3,5-dione. However, the reaction product proved to be meso-(E,E)-1,1'-[1,2-bis(4-chlorophenyl)ethane-1,2-diyl]bis(phenyldiazene), C26H20Cl2N4, and a tentative mechanism is proposed. Crystallization from ethanol produces two concomitant polymorphs, i.e. a triclinic form, (I), in the space group P-1, and a monoclinic form, (II), in the space group C2/c. In both polymorphs, the molecules lie across centres of inversion, but in (II), the molecules are subject to whole-molecule disorder equivalent to configurational disorder with occupancies of 0.6021 (19) and 0.3979 (19). There are no hydrogen bonds in the crystal structure of polymorph (I), but the molecules of polymorph (II) are linked by C-H...π(arene) hydrogen bonds into complex chains, which are further linked into sheets by C-H...N interactions.

Vibrational analysis and chemical activity of paracetamol–oxalic acid cocrystal based on monomer and dimer calculations: DFT and AIM approach

Hydrogen bonding network present in monomer and dimer + 2OXA models of cocrystal.

Invariom approach to electron density studies of open-shell compounds: the case of an organic nitroxide radical

The study reports the successful modeling of electron density with invarioms for an open-shell compound, an organic nitroxide radical.

Stabilisation of metastable polymorphs: the case of paracetamol form III

A unique approach is applied in preparing the first ‘air-stable’ formulation of form III of paracetamol, long regarded as a highly unstable form.

Probing Ionic Crystals by the Invariom Approach: An Electron Density Study of Guanidinium Chloride and Carbonate

A comparative study of two guanidinium salts, chloride and carbonate, is carried out to test the performance of the invariom approach for ionic crystals. Although treating them as formed by isolated ions with no charge transfer between them, the invariom approach provides features of interionic contacts that are amazingly similar to those obtained from conventional charge density analysis of high-resolution X-ray diffraction data, thus emerging as an easy way towards reliable description of chemical bonding peculiarities in ionic crystals.

Invariom approach as a new tool in electron density studies of ionic liquids: a model case of 1-butyl-2,3-dimethylimidazolium chloride BDMIM[Cl]

In this comparative study of an ionic liquid (IL), BDMIM[Cl], the invariom approximation emerges as a useful tool in search for ‘structure–property’ relations in ILs, as entangling complex interplay of interionic interactions that operate in them.

Adamantyl analogues of paracetamol as potent analgesic drugs via inhibition of TRPA1

Paracetamol also known as acetaminophen, is a widely used analgesic and antipyretic agent. We report the synthesis and biological evaluation of adamantyl analogues of paracetamol with important analgesic properties. The mechanism of nociception of compound 6a/b, an analog of paracetamol, is not exerted through direct interaction with cannabinoid receptors, nor by inhibiting COX. It behaves as an interesting selective TRPA1 channel antagonist, which may be responsible for its analgesic properties, whereas it has no effect on the TRPM8 nor TRPV1 channels. The possibility of replacing a phenyl ring by an adamantyl ring opens new avenues in other fields of medicinal chemistry.

Contributions of charge-density research to medicinal chemistry

This article reviews efforts in accurate experimental charge-density studies with relevance to medicinal chemistry. Initially, classical charge-density studies that measure electron density distribution via least-squares refinement of aspherical-atom population parameters are summarized. Next, interaction density is discussed as an idealized situation resembling drug-receptor interactions. Scattering-factor databases play an increasing role in charge-density research, and they can be applied both to small-molecule and macromolecular structures in refinement and analysis; software development facilitates their use. Therefore combining both of these complementary branches of X-ray crystallography is recommended, and examples are given where such a combination already proved useful. On the side of the experiment, new pixel detectors are allowing rapid measurements, thereby enabling both high-throughput small-molecule studies and macromolecular structure determination to higher resolutions. Currently, the most ambitious studies compute intermolecular interaction energies of drug-receptor complexes, and it is recommended that future studies benefit from recent method developments. Selected new developments in theoretical charge-density studies are discussed with emphasis on its symbiotic relation to crystallography.

Insight to the thermodynamic stability of molecular crystals through crystallographic studies of a multipolymorph system

Five solvent-free polymorphs of a pharmaceutical compound were discovered during polymorph screening. Out of the five polymorphs, only one has strong intermolecular N-H···N hydrogen bonding, whereas the others exhibit only weak C-H···N and π-π stacking interactions in addition to all the other weak C-H···X and van der Waals interactions. The relative thermodynamic stability relationships among the polymorphs are not intuitive and quite complex due to enantiotropic phase behavior. For instance, the polymorph with the most efficient packing (i.e., highest density) is not always the most thermodynamically stable form, and the polymorph with strong intermolecular interactions is not thermodynamically more stable than the polymorph with weak intermolecular interactions at all temperatures. Nevertheless, systematic examination and comparison of the molecular packing and intermolecular interactions of these polymorphs provide insight into the importance of H-bonding and packing efficiency to the thermodynamic stability of a crystalline form, and how these effects are dependent on temperature. This study seeks to correlate single-crystal structure features with experimentally established thermodynamic stability, and provides an example where a polymorph with only van der Waals forces and weak intermolecular interactions can be more stable than a polymorph that displays strong H-bonding in its structural make-up.

Polymorphism in paracetamol: evidence of additional forms IV and V at high pressure

The structural phase stability of N-(4-hydroxyphenyl) acetamide (paracetamol) has been studied at ambient temperature up to 23 GPa using Raman spectroscopy. Spectral changes have provided further evidence for a highly kinetically driven Form I → II transition that occurs as a mixed phase from 4.8 to 6.5 GPa, and might complete as early as 7 GPa. Upon further compression to 8.1 GPa, a drastic shift in spectral signature was observed providing the first evidence for a previously undiscovered Form IV of paracetamol. Additional shifts in mode intensities were observed near 11 GPa indicating a potential restructuring of the hydrogen bonding network and/or structural modification to a potentially new Form V. Phase boundaries at 7 and 8 GPa were confirmed under hydrostatic conditions using Raman spectroscopy. Spectral changes indicate that the transition Form IV → V occurs near 11 GPa. Multiple ab initio harmonic frequency calculations at different levels of theory were performed with a B3LYP/6-31G** being used to provide a more robust mode assignment to our experimentally obtained Raman modes. High pressure X-ray diffraction (XRD) was performed up to 21 GPa, which provided further evidence for a highly kinetically driven Form I → II transition in agreement with our Raman measurements. In addition, the XRD provided further evidence for the existence of Form IV near 8 GPa and Form V near 11 GPa with Form V persisting up to 21 GPa.

Extremely long Cu···O contact as a possible pathway for magnetic interactions in Na2Cu(CO3)2

Chemical binding in a mixed copper sodium carbonate Na2Cu(CO3)2, a layered material showing ferromagnetic intralayer exchange and weak antiferromagnetic interlayer coupling, was examined within the topological analysis of experimental (from high-resolution X-ray diffraction) and theoretical (from periodic quantum chemical calculations) electron density functions in its crystal. Together with modeling of a superexchange pathway within the LSDA and DFT+U approach, the results obtained reveal a very weak Cu···O interaction (0.5 kcal/mol worth) between the copper-carbonate layers that is nevertheless stabilizing (bonding) and may serve as a possible pathway for antiferromagnetic interactions.

Energy Analysis of Competing Non-Covalent Interaction in 1:1 and 1:2 Adducts of Collidine with Benzoic Acids by Means of X-Ray Diffraction

The hydrogen bond pattern and the types of non-covalent interactions in the crystals of the 1:1 and 1:2 adducts of 2,4,6-trimethylpyridine and benzoic acids are studied using high-resolution X-ray diffraction. The geometries of the hydrogen bonds are estimated using a combined XRD/DFT approach that provides the geometrical parameters within the margin of error of neutron diffraction studies. The energies of the non-covalent interactions are estimated on the base of the experimental electron density distribution function. It is shown that the structures of the adducts are governed by the NOH and OHO hydrogen bonds. In turn, C-H...O contacts and stacking interactions define the packing of the adducts in the crystal. On the other hand, it is important to note that the latter interactions affect the competition of the former hydrogen bonds in some 1:2 adducts.

Probing stereoelectronic interactions in an O-N-O unit by the atomic energies: experimental and theoretical electron density study

Stereoelectronic interaction lp(O1) → σ*(N1-O2) in a O-N-O unit was analyzed by means of R. Bader's Atoms in Molecule theory on the basis of X-ray diffraction data for dimethyl-(2R,4aR,5S,7R)-2,5,7-triphenylhexahydro-4H-[1,2]oxazino[2,3-b][1,2]oxazine-4,4-dicarboxylate. Atomic energies obtained by applying this approach to both the experimental and theoretical electron densities were used to probe the energy of this strong stereoelectronic interaction, giving consistent results with the NBO analysis, although showing its destabilizing character.