Accurate Lattice Energies for Molecular Crystals from Experimental Crystal Structures

Cellulose nanocrystal-stabilized Pickering emulsion gels as vehicles for follicular delivery of minoxidil

Minoxidil (MXD) is the only topical over-the-counter medicine approved by the United States Food and Drug Administration for the treatment of androgenetic alopecia. For the purpose of targeting the delivery of MXD to dermal papilla in the hair follicle, MXD Pickering emulsion gels were fabricated based on the designability of deep eutectic solvent (DES) and the versatility of cellulose nanocrystal (CNC) and sodium carboxymethyl cellulose (CMC-Na). Structural studies and theoretical calculations results suggest that CNC can stabilize the interface between the MXD-DES and water, leading to the formation of Pickering emulsions. The rheological properties and stabilities of MXD Pickering emulsions were enhanced through gelation using CMC-Na, which highlights the good compatibility and effectiveness of natural polysaccharides in emulsion gels. Due to the particle size of emulsion droplets (679 nm) and the rheological properties of emulsion gel, the fabricated MXD formulations show in vivo hair regrowth promotion and hair follicle targeting capabilities. Interestingly, the MXD Pickering emulsion-based formulations exert therapeutic effects by upregulating the expression of hair growth factors. The proposed nanodrug strategy based on supramolecular strategies of CNC and CMC-Na provides an interesting avenue for androgenetic alopecia treatment.

Designing Organic π-Conjugated Molecules for Crystalline Solid Solutions: Adamantane-Substituted Naphthalenes

We showcase herein organic crystalline solid solutions (CSSs) based on the simplest polycyclic aromatic hydrocarbon (PAH) scaffold, naphthalene, stabilized by dispersion forces induced by adamantane substitution. High thermal stability of the host and guest molecules synthesized by cross-coupling of dibromonaphthalene derivatives and 4-(1-adamantyl)phenyl boronic ester enabled formation of crystals by sublimation. We could generate binary monocrystalline solid solution systems proven by X-ray crystallography, the first system of designed CSSs stabilized exclusively via dispersion forces with structural evidence. These observations are additionally supported by lattice energy calculations and spectroscopic examinations. For the generation of CSSs, it is of utmost importance that the host and guest molecules have similar lattice energies and spatial compatibility. We anticipate that the thermostable organic CSS design demonstrated herein would be beneficial for functional materials and further investigation towards materials with unique properties.

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.

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.

(S)-2-Carb-oxy-ethyl l-cysteinyl sulfone

The title compound {systematic name: (2S)-2-aza-niumyl-3-[(2-carb-oxy-ethane)-sulfon-yl]propano-ate}, C6H11NO6S, forms enanti-opure crystals in the monoclinic space group P21 and exists as a zwitterion, with a protonated α-amino group and a deprotonated α-carboxyl group. Both the carboxyl groups and the amino group are involved in an extensive multicentered inter-molecular hydrogen-bonding scheme. In the crystal, the diperiodic network of hydrogen bonds propagates parallel to (101) and involves inter-connected heterodromic R 4 3(10) rings. Electrostatic forces are major contributors to the structure energy, which was estimated by DFT calculations as E total = -333.5 kJ mol-1.

First-Principles Models of Polymorphism of Pharmaceuticals: Maximizing the Accuracy-to-Cost Ratio

Accuracy and sophistication of in silico models of structure, internal dynamics, and cohesion of molecular materials at finite temperatures increase over time. Applicability limits of ab initio polymorph ranking that would be feasible at reasonable costs currently represent crystals of moderately sized molecules (less than 20 nonhydrogen atoms) and simple unit cells (containing rather only one symmetry-irreducible molecule). Extending the applicability range of the underlying first-principles methods to larger systems with a real-life significance, and enabling to perform such computations in a high-throughput regime represent additional challenges to be tackled in computational chemistry. This work presents a novel composite method that combines the computational efficiency of density-functional tight-binding (DFTB) methods with the accuracy of density-functional theory (DFT). Being rooted in the quasi-harmonic approximation, it uses a cheap method to perform all of the costly scans of how static and dynamic characteristics of the crystal vary with respect to its volume. Such data are subsequently corrected to agree with a higher-level model, which must be evaluated only at a single volume of the crystal. It thus enables predictions of structural, cohesive, and thermodynamic properties of complex molecular materials, such as pharmaceuticals or organic semiconductors, at a fraction of the original computational cost. As the composite model retains the solid physical background, it suffers from a minimum accuracy deterioration compared to the full treatment with the costly approach. The novel methodology is demonstrated to provide consistent results for the structural and thermodynamic properties of real-life molecular crystals and their polymorph ranking.

A multi-chromic boron trifluoride-pyridyl Lewis adduct

A boron trifluoride-pyridyl Lewis adduct is reported, which exhibits various types of chromism and high solid-state photoluminescence quantum yields, as well as excitation-dependent emission in the mechanically ground form. The facile synthetic approach offers a simple and potentially versatile strategy for inducing chromism in pyridyl ligands with donor moieties. We envisage this approach as having a dual benefit: simplicity and extensive applicability.

Releasing a bound molecular spring with light: a visible light-triggered photosalient effect tied to polymorphism

Here we present a study on the solid state properties of trans tetra-ortho-bromo azobenzene (4Br-Azo). Two distinct crystal polymorphs were identified: the α-phase and β-phase. Notably, only the β-phase exhibited an extraordinary photosalient effect (jumping/breaking) upon exposure to a wide range of visible light. Powder X-ray diffraction and Raman spectroscopy revealed that the β-phase is metastable and can transition to the α-phase when subjected to specific stimuli like heat and light. Furthermore, single crystal X-ray diffraction and density functional theory calculations highlighted the significance of a highly strained conformer in the β-phase, showing that the metastability of the phase potentially arises from relieving this strain. This metastability leads to a light induced phase transition, which appears to be the cause of the photosalient effect in these crystals. Interestingly the polymorphism at the core of 4Br-Azo's dynamic behavior is based on different arrangements of halogen based intermolecular interactions. It is possible that continued study on combining visible light capturing chromophores with halogen interaction-based polymorphism will lead to the discovery of even more visible light controlled dynamic crystal materials.

Structural Landscape of α-Acetamidocinnamic Acid Cocrystals with Bipyridine-Based Coformers: Influence of Crystal Packing on Their Thermal and Photophysical Properties

Controlling the supramolecular synthon outcome in systems with different functionalities has been a key factor for the design of supramolecular materials, which also affected their physicochemical properties. In this contribution, we have analyzed the structural landscape of α-acetamidocinnamic acid (HACA) aiming to find its synthon outcome from the competitivity between its acidic and amidic groups. We prepared four multicomponent forms including one dihydrate (HACA·2H2O) and three cocrystals bearing different bipyridine coformers with formulas (HACA)2(1,2-bpe) (1), (HACA)2(4,4'-azpy) (2), and (HACA)2(4,4'-bipy)3 (3) (1,2-bpe = 1,2-bis(4-pyridyl)ethylene; 4,4'-azpy = 4,4'-azopyridine; 4,4'-bipy = 4,4'-bipyridine). First, we applied a virtual screening approach to assess the feasibility of cocrystal formation. Then, we synthesized the cocrystals, via liquid-assisted grinding (LAG) (1 and 2) or solvothermal (3) techniques, and single crystals of HACA, and their four multicomponent forms were obtained showing different synthons and crystal packings. Besides, a Cambridge Structural Database (CSD) search of the cocrystals presenting bipyridine-type coformers and molecules with acid and amide functionalities was performed, and the observed synthon occurrences as well as the possibility of synthon modification by tuning the H-donor/H-acceptor propensity of the acidic and amidic groups were shown. Finally, we measured their thermal and photophysical properties, which were correlated with their structural features.

Crystallization from solution versus mechanochemistry to obtain double-drug multicomponent crystals of ethacridine with salicylic/acetylsalicylic acids

Salicylic and acetylsalicylic acids and ethacridine have complementary bioactive properties. They can be combined to obtain double-drug multicomponent crystals. Their reactivity in different environments was explored to obtain the possible compounds, stable at different hydration degrees. Solution, liquid-assisted grinding, and dry preparation approaches were applied to the couples of reactants in different stoichiometric ratios. Four compounds were obtained, and three out of them were stable and reproducible enough to determine their structures using SCXRD or PXRD methods. When coupled to ethacridine, salicylic acid gave two stable structures (1 and 3, both showing 1:1 ratio but different hydration degree) and a metastable one (5), while acetylsalicylic acid only one structure from solution (2 in 1:1 ratio), while LAG caused hydrolysis and formation of the same compound obtained by LAG of ethacridine with salicylic acid. While solution precipitation gave dihydrated (1) or monohydrated (2) structures with low yields, LAG of salicylic acid and ethacridine allowed obtaining an anhydrous salt complex (3) with a yield close to 1. The structures obtained by solution crystallizations maximize π(acridine)-π(acridine) contacts with a less compact packing, while the LAG structure is more compact with a packing driven by hydrogen bonds. For all compounds, NMR, ATR-FTIR, and Hirshfeld surface analysis and energy framework calculations were performed.

Crystal-structure simulation of molecular semiconductors: brickwork-related crystal structures of methylthiolated peri-condensed polycyclic aromatic hydrocarbons

Despite the critical importance to carrier transport properties, studies on the control and prediction of crystal structures of molecular semiconductors have not been well-matured. To tackle this issue, we have developed "in silico crystallization" (ISC) protocols for simulating the brickwork (BW) crystal structures of methylchalcogenolated polycyclic aromatic hydrocarbons (PAHs). In this study, by carefully analyzing a BW-related polymorph of experimental crystal structures, an inclined brickwork (iBW) structure, we further extend the ISC protocol to simulate various BW-related crystal structures including iBW structures. Rational conditional branching in the simulation not only makes it possible to simulate eight polymorph candidates of methylchalcogenolated PAHs but also helps understand the relationship between the polymorphs. Furthermore, the relative favorability of each polymorphic candidate, i.e., the likelihood of the appearance among the polymorph candidates, can also be evaluated.

Can Spin-Component Scaled MP2 Achieve kJ/mol Accuracy for Cohesive Energies of Molecular Crystals?

Attaining kJ/mol accuracy in cohesive energy for molecular crystals is a persistent challenge in computational materials science. In this study, we evaluate second-order Møller-Plesset perturbation theory (MP2) and its spin-component scaled models for calculating cohesive energies for 23 molecular crystals (X23 data set). Using periodic boundary conditions and Brillouin zone sampling, we converge results to the thermodynamic and complete basis set limits, achieving an accuracy of about 2 kJ/mol (0.5 kcal/mol), which is rarely achieved in previous MP2 calculations for molecular crystals. When compared to experimental data, our results have a mean absolute error of 12.9 kJ/mol, comparable to Density Functional Theory with the PBE functional and TS dispersion correction. By separately scaling the opposite-spin and same-spin correlation energy components, using predetermined parameters, we reduce the mean absolute error to 9.5 kJ/mol. Further fine-tuning of these scaling parameters specifically for the X23 data set brings the mean absolute error down to 7.5 kJ/mol.

A transferable quantum mechanical energy model for intermolecular interactions using a single empirical parameter

The calculation of intermolecular interactions in molecular crystals using model energies provides a unified route to understanding the complex interplay of driving forces in crystallization, elastic properties and more. Presented here is a new single-parameter interaction energy model (CE-1p), extending the previous CrystalExplorer energy model and calibrated using density functional theory (DFT) calculations at the ωB97M-V/def2-QZVP level over 1157 intermolecular interactions from 147 crystal structures. The new model incorporates an improved treatment of dispersion interactions and polarizabilities using the exchange-hole dipole model (XDM), along with the use of effective core potentials (ECPs), facilitating application to molecules containing elements across the periodic table (from H to Rn). This new model is validated against high-level reference data with outstanding performance, comparable to state-of-the-art DFT methods for molecular crystal lattice energies over the X23 set (mean absolute deviation 3.6 kJ mol-1) and for intermolecular interactions in the S66x8 benchmark set (root mean-square deviation 3.3 kJ mol-1). The performance of this model is further examined compared to the GFN2-xTB tight-binding model, providing recommendations for the evaluation of intermolecular interactions in molecular crystal systems.

CrystalClear: an open, modular protocol for predicting molecular crystal growth from solution

We present a new protocol for the prediction of free energies that determine the growth of sites in molecular crystals for subsequent use in Monte Carlo simulations using tools such as CrystalGrower [Hill et al., Chemical Science, 2021, 12, 1126-1146]. Key features of the proposed approach are that it requires minimal input, namely the crystal structure and solvent only, and provides automated, rapid generation of the interaction energies. The constituent components of this protocol, namely interactions between molecules (growth units) in the crystal, solvation contributions and treatment of long-range interactions are described in detail. The power of this method is shown via prediction of crystal shapes for ibuprofen grown from ethanol, ethyl acetate, toluene and acetonitrile, adipic acid grown from water, and five polymorphs (ON, OP, Y, YT04 and R) of ROY (5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile), with promising results. The predicted energies may be used directly or subsequently refined against experimental data, facilitating insight into the interactions governing crystal growth, while also providing a prediction of the solubility of the material. The protocol has been implemented in standalone, open-source software made available alongside this publication.

Novel Adamantane-Linked Isothiourea Derivatives as Potential Chemotherapeutic Agents: Synthesis, Structural Insights, and Antimicrobial/Anti-Proliferative Profiles

In this study, two adamantane-linked isothiourea derivatives containing a common 4-chlorophenyl substituent coupled with 4-nitrobenzyl or 4-bromobenzyl moieties were synthesized. Both derivatives were characterized, in the solid state and in solution, through a synergistic combination of experimental and in silico techniques, and the results are of great value for the chemical and structural characterization of related compounds. The crystal structures of both derivatives were analyzed in depth, including Hirshfeld surface analysis and lattice energy calculations, revealing a predominant dispersive component of the total energy that stabilizes crystal packing. Both compounds showed potent broad-spectrum antibacterial activity and moderate activity against the pathogenic fungus Candida albicans. In addition, in vitro anti-proliferative activity assays showed that the 4-bromobenzyl analogue displays higher activity than the 4-nitrobenzyl one, with IC₅₀ values under 30 μM against five human cancer cell lines. Our results give evidence of the potential of the adamantane/isothiourea combination to render auspicious scaffolds for new potential chemotherapeutic agents.

Polymorphism of Pradofloxacin: Crystal Structure Analysis, Stability Study, and Phase Transformation Behavior

PURPOSE

Pradofloxacin is an important antibiotic with poor physical stability. At present, there is no systematic study on its polymorphic form. The purpose of this study is to develop new crystal forms to improve the stability of Pradofloxacin and systematically study the crystal transformation relationships to guide industrial production.

METHOD

In this work, three solvent-free forms (Form A, Form B and Form C), a new dimethyl sulfoxide solvate (Form PL-DMSO) and a new hydrate (Form PL-H) were successfully obtained and the single crystal data of Form A, Form B and Form PL-DMSO were solved for the first time. Various solid state analysis techniques and slurry experiments have been used to evaluate the stability and determine phase transformation relationships of five crystal forms, the analysis of crystal structure provided theoretical support for the results.

RESULT

The water vapor adsorption and desorption experiences of Forms A, B, C and Form PL-H were studied, and the results show that the new hydrate has good hygroscopic stability and certain development potential. The thermal stability of different forms was determined by thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC) and the crystal structure shows that there are more hydrogen bonds and C - H···π interactions in form B, which is the reason why Form B is more stable than form A. Finally, the phase transformation relationships of the five crystal forms were systematically studied and discussed.

CONCLUSION

These results are helpful to provide guiding methods in the production and storage of pradofloxacin.

XDM-corrected hybrid DFT with numerical atomic orbitals predicts molecular crystal lattice energies with unprecedented accuracy

Molecular crystals are important for many applications, including energetic materials, organic semiconductors, and the development and commercialization of pharmaceuticals. The exchange-hole dipole moment (XDM) dispersion model has shown good performance in the calculation of relative and absolute lattice energies of molecular crystals, although it has traditionally been applied in combination with plane-wave/pseudopotential approaches. This has limited XDM to use with semilocal functional approximations, which suffer from delocalization error and poor quality conformational energies, and to systems with a few hundreds of atoms at most due to unfavorable scaling. In this work, we combine XDM with numerical atomic orbitals, which enable the efficient use of XDM-corrected hybrid functionals for molecular crystals. We test the new XDM-corrected functionals for their ability to predict the lattice energies of molecular crystals for the X23 set and 13 ice phases, the latter being a particularly stringent test. A composite approach using a XDM-corrected, 25% hybrid functional based on B86bPBE achieves a mean absolute error of 0.48 kcal mol^-1 per molecule for the X23 set and 0.19 kcal mol^-1 for the total lattice energies of the ice phases, compared to recent diffusion Monte-Carlo data. These results make the new XDM-corrected hybrids not only far more computationally efficient than previous XDM implementations, but also the most accurate density-functional methods for molecular crystal lattice energies to date.

Crystal structure of hydrazinium methanesulfonate, CH8N2O3S

CH8N2O3S, monoclinic, P21/c (no. 14), a = 9.7583 (16) Å, b = 5.4033 (9) Å, c = 10.4729 (17) Å, β = 110.483 (4)°, V = 517.29 (15) Å3, Z = 1, R gt(F) = 0.0225, wR ref(F 2) = 0.0649, T = 150 K.

Low-Molecular-Weight Organogelators Based on N-dodecanoyl-L-amino Acids-Energy Frameworks and Supramolecular Synthons

Lauric acid was used to synthesize the low-molecular-weight organogelators (LMOGs), derivatives of two endogenous (L)-alanine, (L)-leucine, and three exogenous (L)-valine, (L)-phenylalanine, and (L)-proline amino acids. The nature of processes responsible for the gel formation both in polar and in apolar solvents of such compounds is still under investigation. Knowing that the organization of surfactant molecules affects the properties of nano scale materials and gels, we decided to elucidate this problem using crystallographic diffraction and energy frameworks analysis. The single crystals of the mentioned compounds were produced successfully from heptane/^tBuOMe mixture. The compounds form lamellar self-assemblies in crystals. The energetic landscapes of single crystals of a series of studied amphiphilic gelators have been analyzed to explore the gelling properties. The presented results may be used as model systems to understand which supramolecular interactions observed in the solid state and what energy contributions are desired in the designing of new low-molecular-weight organic gelators.

Co-Crystallization Approach to Enhance the Stability of Moisture-Sensitive Drugs

Stability is an essential quality attribute of any pharmaceutical formulation. Poor stability can change the color and physical appearance of a drug, directly impacting the patient's perception. Unstable drug products may also face loss of active pharmaceutical ingredients (APIs) and degradation, making the medicine ineffective and toxic. Moisture content is known to be the leading cause of the degradation of nearly 50% of medicinal products, leading to impurities in solid dose formulations. The polarity of the atoms in an API and the surface chemistry of API particles majorly influence the affinity towards water molecules. Moisture induces chemical reactions, including free water that has also been identified as an important factor in determining drug product stability. Among the various approaches, crystal engineering and specifically co-crystals, have a proven ability to increase the stability of moisture-sensitive APIs. Other approaches, such as changing the salt form, can lead to solubility issues, thus making the co-crystal approach more suited to enhancing hygroscopic stability. There are many reported studies where co-crystals have exhibited reduced hygroscopicity compared to pure API, thereby improving the product's stability. In this review, the authors focus on recent updates and trends in these studies related to improving the hygroscopic stability of compounds, discuss the reasons behind the enhanced stability, and briefly discuss the screening of co-formers for moisture-sensitive drugs.

Unusual shape-preserved pathway of a core-shell phase transition triggered by orientational disorder

The ubiquitous presence of crystal defects provides great potential and opportunities to construct the desired structure (hence with the desired properties) and tailor the synthetic process of crystalline materials. However, little is known about their regulation role in phase transition and crystallization pathways. It was generally thought that a phase transition in solution proceeds predominantly via the solvent-mediated phase-transformation pathway due to energetically high-cost solid-state phase transitions (if any). Herein, we report an unprecedented finding that an orientational disorder defect present in the crystal structure triggers an unusual pathway of a core-shell phase transition with apparent shape-preserved evolution. In the pathway, the solid-state dehydration phase transition occurs inside the crystal prior to its competitive transformation approach mediated by solvent, forming an unconventional core-shell structure. Through a series of combined experimental and computational techniques, we revealed that the presence of crystal defects, introduced by urate tautomerism over the course of crystallization, elevates the metastability of uric acid dihydrate (UAD) crystals and triggers UAD dehydration to the uric acid anhydrate (UAA) phase in the crystal core which precedes with surface dissolution of the shell UAD crystal and recrystallization of the core phase. This unique phase transition could also be related to defect density, which appears to be influenced by the thickness of UAD crystals and crystallization driving force. The discovery of an unusual pathway of the core-shell phase transition suggests that the solid-state phase transition is not necessarily slower than the solvent-mediated phase transformation in solution and provides an alternative approach to constructing the core-shell structure. Moreover, the fundamental role of orientational disorder defects on the phase transition identified in this study demonstrates the feasibility to tailor phase transition and crystallization pathways by strategically importing crystal defects, which has broad applications in crystal engineering.

Polar and non-polar stacking of perfectly aligned parallel beloamphiphile monolayers (PBAMs) of (PhO, F)-azine. The interplay of non-covalent interlayer interactions and unit cell polarity

The differences are discussed of the antiferroelectric and ferroelectric stacking of the PBAMs of the polymorphs of (PhO, F)-azine. We will show how non-covalent interlayer H⋯F and F⋯F interactions between the PBAM surfaces affect their stacking.

Quantitative Analysis the Weak Non-Covalent Interactions of the Polymorphs of Donepezil

Donepezil has polymorphism. Different crystalline forms can exhibit different physicochemical properties and biological activities. Exploration of intermolecular interactions is essential to reveal the formation mechanism and differences in properties of polymorphs. This study explores the weak non-covalent intermolecular interactions of donepezil polymorphs through fully ab initio quantum mechanical methods, semi-empirical methods, and Hirshfeld surface analysis. The results show that the Hirshfeld surface analysis method can clearly and intuitively reveal the intermolecular interactions. Theoretical calculations using the atom-atom Coulomb-London-Pauli (AA-CLP) method were also performed to understand the interaction energies toward the total lattice energy. The value of the lattice energy was in accordance with the melting points of the donepezil polymorphs and brought to light the nature of thermal stability. In the specific energy distribution, the contribution of the dispersion force is the most prominent. Further interaction energy analysis found that within a distance of 3.8 Å from the center of the donepezil molecule, different crystalline forms of donepezil molecules have different interaction energies with surrounding molecules. The different interaction energies between polymorphs may lead to polymorphs with different physical-chemical properties.

An experimental and theoretical charge density study of theophylline and malonic acid cocrystallization

The pharmaceutical agent theophylline (THEO) is primarily used as a bronchodilator and is commercially available in both tablet and liquid dosage forms. THEO is highly hygroscopic, reducing its stability, overall shelf-life, and therefore usage as a drug. THEO and dicarboxylic acid cocrystals were designed by Trask et al. in an attempt to decrease the hygroscopic behaviour of THEO; cocrystallisation of THEO with malonic acid (MA) did not improve the hygroscopic stability of THEO in simulated atmospheric humidity testing. The current study employed high-resolution X-ray crystallography, and Density Functional Theory (DFT) calculations to examine the electron density distribution (EDD) changes between the cocrystal and its individual components. The EED changes identified the reasons why the THEO:MA cocrystal did not alter the hygroscopic profile of THEO. The cocrystal was equally porous, with atomic packing factors (APF) similar to those of THEO 0.73 vs. 0.71, respectively. The THEO:MA (1) cocrystal structure is held together by an array of interactions; a heterogeneous synthon between the imidazole and a carboxylic fragment stabilising the asymmetric unit, a pyrimidine-imidazole homosynthon, and an aromatic cycle stack between two THEO moieties have been identified, providing 9.7-12.9 kJ mol-1 of stability. These factors did not change the overall relative stability of the cocrystal relative to its individual THEO and MA components, as shown by cocrystal (1) and THEO being equally stable, with calculated lattice energies within 2.5 kJ mol-1 of one other. The hydrogen bond analysis and fragmented atomic charge analysis highlighted that the formation of (1) combined both the EDD of THEO and MA with no net chemical change, suggesting that the reverse reaction - (1) back to THEO and MA - is of equal potential, ultimately producing THEO hydrate formation, in agreement with the work of Trask et al. These results highlight that a review of the EDD change associated with a chemical reaction can aid in understanding cocrystal design. In addition, they indicate that cocrystal design requires further investigation before becoming a reliable process, with particular emphasis on identifying the appropriate balance of synthon engineering, weak interactions, and packing dynamics.

Investigation of crystal structures, energetics and isostructurality in halogen-substituted phosphoramidates

A total of 14 compounds, one unsubstituted and 13 halogen-substituted phosphoramidates, have been synthesized from unsubstituted and halogenated (fluoro-, difluoro-, chloro-, bromo-, iodo-substituted) aniline and diphenyl phosphoryl chloride to investigate their molecular assembly in solid-state structures. Amongst them, six groups were formed based on similarities in unit-cell dimensions, space group and molecular assembly of the crystal. The analysis reveals that all the crystal structures contain robust N-H...O hydrogen bonds which are the primary building blocks with ancillary interactions such as C-H...O, C-H...π, C-H...F/Cl/Br/I, F...F, F...π, I...π, Br...π, I...O and Br...O. The role of short and directional C-H...O and C-H...π interactions providing significant stabilization to the densely packed crystalline arrangement is discussed. The contribution of these interactions in stabilizing the crystalline assembly was deduced via computing total interaction energy between dimers and the overall lattice energies using the computer programs Crystal Explorer 17.5 and PIXELC, respectively. Additionally, the occurrence of 3D isostructurality in phosphoradimates and their halogenated analogs was investigated using the XPac program. A comparison of the magnitudes of the torsion angles in the compounds substantiates the role of conformational flexibility in the solid state.

Structural studies of various olmesartan solvates

Seven solvates of the angiotensin II receptor blocker agent olmesartan (C24H26N6O3), namely, the methanol (C24H26N6O3·CH4O), ethanol (C24H26N6O3·C2H6O), isopropanol (C24H26N6O3·C3H8O), isobutanol (C24H26N6O3·C4H10O), 2-ethoxyethanol (C24H26N6O3·C4H10O2), chloroform (C24H26N6O3·CHCl3) and acetonitrile (C24H26N6O3·C2H3N) solvates, were successfully obtained. The crystal structures were determined using the single-crystal X-ray diffraction technique and the structural features are described, each solvate containing one molecule of olmesartan and one of solvent in the asymmetric unit. The samples were also analyzed by powder X-ray diffraction. Total lattice energies and binding energies between the olmesartan and solvent molecules were evaluated, which can be partitioned into electrostatic, polarization, dispersion and repulsion components. Hirshfeld and fingerprint plot analysis was performed to highlight the intermolecular contacts. Hydrogen bonding and supramolecular arrangements were comparatively studied for the seven solvates.

Long-Range Supramolecular Synthon Isomerism: Insight from a Case Study of Vinylic Tellurium Trihalides Cl(Ph)C=C(Ph)TeX3 (X = Cl, I)

A slight modification of the synthetic procedure resulted in a new (Cc) polymorph of vinylic tellurium-trichloride Z-Cl(Ph)C=C(Ph)TeCl3 (1, β-form) which is stabilized by Te⋯Cl chalcogen bonds, assembling its molecules into the zigzag chains. Such a packing motive is in contrast to the known (Pca21) polymorph of Z-Cl(Ph)C=C(Ph)TeCl3 (1, α-form, CCDC refcode: BESHOW), which is built upon Te⋯π(Ph) chalcogen bonded chains. We noted a similar case of [Te⋯halogen] vs. [Te⋯π(Ph)] supramolecular synthon polymorphism in its triiodide congener Z-Cl(Ph)C=CPh(TeI3) (2, α and β-polymorphic forms). Quantum chemical calculations of the intermolecular interaction and lattice energies for 1α–β and 2α–β supported the assumption that α is thermodynamic while β is a kinetic form. Kinetic forms 1β and 2β are isostructural (Cc), while the thermodynamic forms 1α (Pca21) and 2α (P21/c) are not and feature an unusual example of long-range supramolecular synthon module isomerism. In other words, 1α–2α pairs demonstrate very similarly to isostructural Te⋯πPh ChB stabilized chains, which are further packed differently relative to each other, following different angular geometry of type-I Cl⋯Cl and type-II I⋯I halogen bonding. These structural considerations are backed by quantum chemical calculations that support the proposed hierarchy of primary and secondary supramolecular synthons and the assignment of α and β as thermodynamic and kinetic forms, respectively.

Macrocyclic "tet a"-Derived Cobalt(III) Complex with a N,N'-Disubstituted Hexadentate Ligand: Crystal Structure, Photonuclease Activity, and as a Photosensitizer

A cobalt(III) complex, [Co(L)]Cl (complex 1, where L = 1,8-[N,N-bis{(3-formyl-2-hydroxy-5-methyl)benzyl}]-1,4,8,11-tetraaza-5,5,7,12,12,14-hexamethylcyclotetradecane) with distorted octahedral geometry has been synthesized and characterized using various spectroscopic techniques. The structure of the ligand has remarkably rich hydrogen intermolecular interactions such as H···H, H···C/C···H, and H···O/O···H that vary with the presence of the metal ion, and the structure of complex 1 has Cl···H interactions; this result has been proved by Hirshfeld surface and two-dimensional (2D) fingerprint maps analyses. The complex exhibits a quasi-reversible Co(III)/Co(II) redox couple with E _1/2 = -0.76 V. Calf thymus DNA (CT DNA) binding abilities of the ligand and complex 1 were confirmed by spectroscopic and electrochemical analyses. According to absorption studies, the ligand and complex 1 bind to CT DNA via intercalative binding mode, with intrinsic binding strengths of 1.41 × 10³ and 8.64 × 10³ M^-1, respectively. A gel electrophoresis assay shows that complex 1 promotes the pUC19 DNA cleavage under dark and light irradiation conditions. Complex 1 has superior antimicrobial activity than the ligand. The cytotoxicity of complex 1 was tested against MDA-MB-231 breast cancer cells with values of IC₅₀ of 1.369 μg mL^-1 in the dark and 0.9034 μg mL^-1 after light irradiation. Besides, cell morphological studies confirmed the morphological changes with AO/EB dual staining, reactive oxygen species (ROS) staining, mitochondria staining, and Hoechst staining on MDA-MB-231 cancer cells by fluorescence microscopy. Complex 1 was found to be a potent antiproliferative agent against MDA-MB-231 cells, and it can induce mitochondrial-mediated and caspase-dependent apoptosis with activation of downregulated caspases. The biotoxicity assay of complex 1 on the development of Artemia nauplii was evaluated at an IC₅₀ value of 200 μg mL^-1 and with excellent biocompatibility.

Crystal structure of (E)-N′-(1-(2-hydroxy-4-methoxyphenyl)ethylidene) isonicotinohydrazide, C15H15N3O3

C15H15N3O3, triclinic, P 1 ‾ $\overline{1}$ (no. 2), a = 9.0752(2) Å, b = 9.1969(2) Å, c = 9.9519(2) Å, α = 112.7059(8)°, β = 92.4309(8)°, γ = 115.2776(7)°, V = 671.16(3) Å3, Z = 2, R gt (F) = 0.0347, wR ref (F 2) = 0.1113, T = 100 K.

Structural systematics in the isomorphous binary co-crystal solvates comprising 2,2'-dithiodibenzoic acid, 4-halobenzoic acid and dimethylformamide (1:1:1), for halide = chloride, bromide and iodide

The 1:1:1 binary co-crystal solvates formulated as 2,2'-dithiodibenzoic acid (DTBA), 4-halobenzoic acid (4-XBA) and dimethylformamide (DMF) for X = Cl (1), Br (2) and I (3) are isomorphous and the...

Effect of substituents in novel bioactive tavaborole derivatives on the intermolecular interaction hierarchy

Tavaborole, a molecule based on the benzoxaborole scaffold, is an effective antifungal drug marketed under the Kerydin® trademark.

Assembling photoactive materials from polycyclic aromatic hydrocarbons (PAHs): room temperature phosphorescence and excimer-emission in co-crystals with 1,4-diiodotetrafluorobenzene

Co-crystallization of PAHs with a polyhalogenated co-former afforded three novel co-crystals, which display remarkable features such as mechanochemical interconversion, photoreactivity, excimer fluorescence, and RTP phosphorescence in the solid state.

Conformational preferences in a series of α-hydroxy ketone derivatives: interplay of conformational energies and lattice cohesive energies

A detailed structural analysis of a series of α-hydroxy ketone derivatives in terms of lattice energies and conformational energies reveals the competing factors of inter- and intramolecular interactions that direct conformational preferences.

Exploring the molecular landscape of multicomponent crystals formed by naproxen drug and acridines

Three cocrystals were obtained by naproxen and acridines, optimizing the yield to more than 99% with LAG. The two structures by solution show a host-guest structure, while that by LAG a layered one, with no interconversion between parent structures.

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.

Metallization-Prompted Robust Porphyrin-Based Hydrogen-Bonded Organic Frameworks for Photocatalytic CO2 Reduction

Under topological guidance, the self-assembly process based on a tetratopic porphyrin synthon results in a hydrogen-bonded organic framework (HOF) with the predicted square layers topology (sql) but unsatisfied stability. Strikingly, simply introducing a transition metal in the porphyrin center does not change the network topology but drastically causes noticeable change on noncovalent interaction, orbital overlap, and molecular geometry, therefore ultimately giving rise to a series of metalloporphyrinic HOFs with high surface area, and excellent stability (intact after being soaked in boiling water, concentrated HCl, and heated to 270 °C). On integrating both photosensitizers and catalytic sites into robust backbones, this series of HOFs can effectively catalyze the photoreduction of CO₂ to CO, and their catalytic performances greatly depend on the chelated metal species in the porphyrin centers. This work enriches the library of stable functional HOFs and expands their applications in photocatalytic CO₂ reduction.

Understanding Hygroscopicity of Theophylline via a Novel Cocrystal Polymorph: A Charge Density Study

The charge density distribution in a novel cocrystal (1) complex of 1,3-dimethylxanthine (theophylline) and propanedioic acid (malonic acid) has been determined. The molecules crystallize in the triclinic, centrosymmetric space group P1̅, with four independent molecules (Z = 4) in the asymmetric unit (two molecules each of theophylline and malonic acid). Theophylline has a notably high hygroscopic nature, and numerous cocrystals have shown a significant improvement in stability to humidity. A charge density study of the novel polymorph has identified interesting theoretical results correlating the stability enhancement of theophylline via cocrystallization. Topological analysis of the electron density highlighted key differences (up to 17.8) in Laplacian (∇²ρ) between the experimental (EXP) and single-point (SP) models, mainly around intermolecular-bonded carbonyls. Further investigation via molecular electrostatic potential maps reaffirmed that the charge redistribution enhanced intramolecular hydrogen bonding, predominantly for N(2') and N(2) (61.2 and 61.8 kJ mol^-1, respectively). An overall weaker lattice energy of the triclinic form (-126.1 kJ mol^-1) compared to that of the monoclinic form (-133.8 kJ mol^-1) suggests a lower energy threshold to overcome to initiate dissociation. Future work via physical testing of the novel cocrystal in both dissolution and solubility will further solidify the correlation between theoretical and experimental results.

Enantioseparation of P-Stereogenic Secondary Phosphine Oxides and Their Stereospecific Transformation to Various Tertiary Phosphine Oxides and a Thiophosphinate

Secondary phosphine oxides incorporating various aryl and alkyl groups were synthesized in racemic form, and these products formed the library reported in this study. TADDOL derivatives were used to obtain the optical resolution of these P-stereogenic secondary phosphine oxides. The developed resolution method showed a good scope under the optimized reaction conditions, as 9 out of 14 derivatives could be prepared with an enantiomeric excess (ee) ≥ 79% and 5 of these derivatives were practically enantiopure >P(O)H compounds (ee ≥ 98%). The scalability of this resolution method was also demonstrated. Noncovalent interactions responsible for the formation of diasteromeric complexes were elucidated by single-crystal XRD measurements. (S)-(2-Methylphenyl)phenylphosphine oxide was transformed to a variety of P-stereogenic tertiary phosphine oxides and a thiophosphinate in stereospecific Michaelis–Becker, Hirao, or Pudovik reactions.

A Universal Force Field for Materials, Periodic GFN-FF: Implementation and Examination

In this study, the adaption of the recently published molecular GFN-FF for periodic boundary conditions (pGFN-FF) is described through the use of neighbor lists combined with appropriate charge sums to handle any dimensionality from 1D polymers to 2D surfaces and 3D solids. Numerical integration over the Brillouin zone for the calculation of π bond orders of periodic fragments is also included. Aside from adapting the GFN-FF method to handle periodicity, improvements to the method are proposed in regard to the calculation of topological charges through the inclusion of a screened Coulomb term that leads to more physical charges and avoids a number of pathological cases. Short-range damping of three-body dispersion is also included to avoid collapse of some structures. Analytic second derivatives are also formulated with respect to both Cartesian and strain variables, including prescreening of terms to accelerate the dispersion/coordination number contribution to the Hessian. The modified pGFN-FF scheme is then applied to a wide range of different materials in order to examine how well this universal model performs.

Dynamics of proton, ion, molecule, and crystal lattice in functional molecular assemblies

Dynamic molecular processes, such as short- or long-range proton (H⁺) and ion (M⁺) motions, and molecular rotations in electrical conducting and magnetic molecular assemblies enable the fabrication of electron-H⁺ (or M⁺) coupling systems, while crystal lattice dynamics and molecular conformation changes in hydrogen-bonded molecular crystals have been utilised in external stimuli responsive reversible gas-induced gate opening and molecular adsorption/desorption behavior. These dynamics of the polar structural units are responsible for the dielectric measurements. The H⁺ dynamics are formed from ferroelectrics and H⁺ conductors, while the dynamic M⁺ motions of Li⁺ and Na⁺ involve ionic conductors and coupling to the conduction electrons. In n-type organic semiconductors, the crystal lattices are modulated by replacing M⁺ cations, with cations such as Li⁺, Na⁺, K⁺, Rb⁺, and Cs⁺. The use of polar rotator or inversion structures such as alkyl amides, m-fluoroanilinium cations, and bowl-shaped trithiasumanene π-cores enables the formation of ferroelectric molecular assemblies. The host-guest molecular systems of ESIPT fluorescent chromic molecules showed interesting molecular sensing properties using various bases, where the dynamic transformation of the crystal lattice and the molecular conformational change were coupled to each other.