Structural Properties, Order-Disorder Phenomena, and Phase Stability of Orotic Acid Crystal Forms

Bifonazole caffeate: The first molecular salt of bifonazole with enhanced biopharmaceutical property based on experiments and quantum chemistry research

In order to make novel breakthroughs in molecular salt studies of BCS class-IV antifungal medication bifonazole (BIF), a salification-driven strategy towards ameliorating attributes and aiding augment efficiency is raised. This strategy fully harnesses structural characters together attributes and benefits of caffeic acid (CAF) to concurrently enhance dissolvability and permeability of BIF by introducing the two ingredients into the identical molecular salt lattice through the salification reaction, which, coupled with the aroused potential activity of CAF significantly amplifies the antifungal efficacy of BIF. Guided by this route, the first BIF-organic molecular salt, BIF-CAF, is directionally designed and synthesized with satisfactorily structural characterizations and integrated theoretical and experimental explorations on the pharmaceutical properties. Single-crystal X-ray diffraction resolving confirms that there is a lipid-water amphiphilic sandwich structure constructed by robust charge-assistant hydrogen bonds in the salt crystal, endowing the molecular salt with the potential to enhance both dissolvability and permeability relative to the parent drug, which is validated by experimental evaluations. Remarkably, the comprehensive DFT-based theoretical investigations covering frontier molecular orbital, molecular electrostatic potential, Hirshfeld surface analysis, reduced density gradient, topology, sphericity and planarity analysis strongly support these observations, thereby allowing some positive relationships between macroscopic properties and microstructures of the molecular salt can be made. Intriguingly, the optimal properties, together with the stimulated activity of CAF markedly augment in vitro antifungal ability of the molecular salt, with magnifying inhibition zones and reducing minimum inhibitory concentrations. These findings fill in the gaps on researches of BIF-organic molecular salt, and adequately exemplify the feasibility and validity by integrating theoretical and experimental approaches to resolve BIF's problems via the salification-driven tactic.

Cocrystallization-driven self-assembly with vanillic acid offers a new opportunity for surmounting fast and excessive absorption issues of antifungal drug 5-fluorocytosine: a combined theoretical and experimental research

The cocrystal of 5-fluorocytosine (FCY) with vanillic acid (VAA) was assembled via a cocrystallization technique, giving a novel understanding for conquering the dose-limited hepatotoxicity caused by the rapid and almost complete absorption of FCY.

Single-Crystal X-ray and Solid-State NMR Characterisation of AND-1184 and Its Hydrochloride Form

In this study, we report on a structural investigation of AND-1184, with the chemical name N-[3-[4-(6-fluoro-1,2-benzoxazol-3-yl)piperidin-1-yl]propyl]-3-methylbenzenesulfonamide (MBS), and its hydrochloride form (MBS^HCl); AND-1184 is a potential API for the treatment of dementia. The single-crystal X-ray investigation of both forms results in monoclinic crystal systems with P2₁/c and C2/c symmetry for MBS and MBS^HCl, respectively. This solid-state NMR study, combined with quantum-chemical calculations, allowed us to assign all ¹³C and most ¹H signals. The MBS structure was defined as a completely rigid system without significant dynamic behaviours, whereas MBS^HCl exhibited limited dynamic motion of the aromatic part of the molecule.

Pharmaceutical Hydrates Analysis-Overview of Methods and Recent Advances

This review discusses a set of instrumental and computational methods that are used to characterize hydrated forms of APIs (active pharmaceutical ingredients). The focus has been put on highlighting advantages as well as on presenting some limitations of the selected analytical approaches. This has been performed in order to facilitate the choice of an appropriate method depending on the type of the structural feature that is to be analyzed, that is, degree of hydration, crystal structure and dynamics, and (de)hydration kinetics. The presented techniques include X-ray diffraction (single crystal X-ray diffraction (SCXRD), powder X-ray diffraction (PXRD)), spectroscopic (solid state nuclear magnetic resonance spectroscopy (ssNMR), Fourier-transformed infrared spectroscopy (FT-IR), Raman spectroscopy), thermal (differential scanning calorimetry (DSC), thermogravimetric analysis (TGA)), gravimetric (dynamic vapour sorption (DVS)), and computational (molecular mechanics (MM), Quantum Mechanics (QM), molecular dynamics (MD)) methods. Further, the successful applications of the presented methods in the studies of hydrated APIs as well as studies on the excipients' influence on these processes have been described in many examples.

NMR crystallography of molecular organics

Developments of NMR methodology to characterise the structures of molecular organic structures are reviewed, concentrating on the previous decade of research in which density functional theory-based calculations of NMR parameters in periodic solids have become widespread. With a focus on demonstrating the new structural insights provided, it is shown how "NMR crystallography" has been used in a spectrum of applications from resolving ambiguities in diffraction-derived structures (such as hydrogen atom positioning) to deriving complete structures in the absence of diffraction data. As well as comprehensively reviewing applications, the different aspects of the experimental and computational techniques used in NMR crystallography are surveyed. NMR crystallography is seen to be a rapidly maturing subject area that is increasingly appreciated by the wider crystallographic community.

Crystal structures of two furazidin polymorphs revealed by a joint effort of crystal structure prediction and NMR crystallography

This work presents the crystal structure determination of two elusive polymorphs of furazidin, an antibacterial agent, employing a combination of crystal structure prediction (CSP) calculations and an NMR crystallography approach. Two previously uncharacterized neat crystal forms, one of which has two symmetry-independent molecules (form I), whereas the other one is a Z' = 1 polymorph (form II), crystallize in P2₁/c and P1 space groups, respectively, and both are built by different conformers, displaying different intermolecular interactions. It is demonstrated that the usage of either CSP or NMR crystallography alone is insufficient to successfully elucidate the above-mentioned crystal structures, especially in the case of the Z' = 2 polymorph. In addition, cases of serendipitous agreement in terms of ¹H or ¹³C NMR data obtained for the CSP-generated crystal structures different from the ones observed in the laboratory (false-positive matches) are analyzed and described. While for the majority of analyzed crystal structures the obtained agreement with the NMR experiment is indicative of some structural features in common with the experimental structure, the mentioned serendipity observed in exceptional cases points to the necessity of caution when using an NMR crystallography approach in crystal structure determination.

Periodic DFT Calculations-Review of Applications in the Pharmaceutical Sciences

In the introduction to this review the complex chemistry of solid-state pharmaceutical compounds is summarized. It is also explained why the density functional theory (DFT) periodic calculations became recently so popular in studying the solid APIs (active pharmaceutical ingredients). Further, the most popular programs enabling DFT periodic calculations are presented and compared. Subsequently, on the large number of examples, the applications of such calculations in pharmaceutical sciences are discussed. The mentioned topics include, among others, validation of the experimentally obtained crystal structures and crystal structure prediction, insight into crystallization and solvation processes, development of new polymorph synthesis ways, and formulation techniques as well as application of the periodic DFT calculations in the drug analysis.

Enhancement of the properties of biosourced poly(3-hydroxybutyrate-co-4-hydroxybutyrate) by the incorporation of natural orotic acid

The aim of this study is to use natural orotic acid (OA) as a sustainable, environmentally friendly additive to improve the crystallization, rheological, thermal, mechanical, and biodegradation properties of bacterially synthesized poly(3-hydroxybutyrate-co-4-hydroxybutyrate) (P34HB). OA was found to be an efficient nucleating agent for P34HB, and dramatically enhanced both non-isothermal and isothermal crystallization rates. The incorporation of OA increased nucleation density and decreased spherulite size, but had little effect on the crystalline structure. The rheological properties of the P34HB were greatly improved by the solid filler OA, particularly when a percolation network structure was formed in the blends. The thermal stability of P34HB was strongly enhanced, as exemplified by the ~23 °C increase in the onset thermal decomposition temperature (T_o) for the blend loaded with 5 wt% OA compared to that of pure P34HB. Moreover, the yield strength and elongation at break of P34HB containing 0.5 wt% OA increased by 25% and 119%, respectively. The most intriguing result was the clear enhancement in the enzymatic hydrolysis rates of the P34HB/OA blends compared to that of neat P34HB. The synergetic improvement in these properties may be of significant importance for the wider practical application of biosourced P34HB.

Solvent driven phase transitions of acyclovir – the role of water and solvent polarity

The pathways of transformations of acyclovir forms I and V induced by organic solvents and water have been identified. Significant differences in the thermal dehydration process of forms V and VI were observed.

Highly Water-Soluble Orotic Acid Nanocrystals Produced by High-Energy Milling

Orotic acid (OA), a heterocyclic compound also known as vitamin B13, has shown potent antimalarial and cardiac protection activities; however, its limited water solubility has posed a barrier to its use in therapeutic approaches. Aiming to overcome this drawback, OA freeze-dried nanocrystal formulations (FA, FB, and FC) were developed by using the high-energy milling method. Polysorbate 80 (FA) and povacoat^® (FC) were used alone and combined (FB) as stabilizers. Nanocrystals were fully characterized by dynamic light scattering, laser diffraction, transmission electron microscopy, thermal analysis (thermogravimetry and derivative thermogravimetry, and differential scanning calorimetry), and X-ray powder diffraction revealing an acceptable polydispersity index, changes in the crystalline state with hydrate formation and z-average of 100-200 nm, a remarkable 200-time reduction compared to the OA raw material (44.3 μm). Furthermore, saturation solubility study showed an improvement of 13 times higher than the micronized powder. In addition, cytotoxicity assay revealed mild toxicity for the FB and FC formulations prepared with povacoat^®. OA nanocrystal platform can deliver innovative products allowing untapped the versatile potential of this drug substance candidate.

GAtor: A First-Principles Genetic Algorithm for Molecular Crystal Structure Prediction

We present the implementation of GAtor, a massively parallel, first-principles genetic algorithm (GA) for molecular crystal structure prediction. GAtor is written in Python and currently interfaces with the FHI-aims code to perform local optimizations and energy evaluations using dispersion-inclusive density functional theory (DFT). GAtor offers a variety of fitness evaluation, selection, crossover, and mutation schemes. Breeding operators designed specifically for molecular crystals provide a balance between exploration and exploitation. Evolutionary niching is implemented in GAtor by using machine learning to cluster the dynamically updated population by structural similarity and then employing a cluster-based fitness function. Evolutionary niching promotes uniform sampling of the potential energy surface by evolving several subpopulations, which helps overcome initial pool biases and selection biases (genetic drift). The various settings offered by GAtor increase the likelihood of locating numerous low-energy minima, including those located in disconnected, hard to reach regions of the potential energy landscape. The best structures generated are re-relaxed and re-ranked using a hierarchy of increasingly accurate DFT functionals and dispersion methods. GAtor is applied to a chemically diverse set of four past blind test targets, characterized by different types of intermolecular interactions. The experimentally observed structures and other low-energy structures are found for all four targets. In particular, for Target II, 5-cyano-3-hydroxythiophene, the top ranked putative crystal structure is a Z' = 2 structure with P1̅ symmetry and a scaffold packing motif, which has not been reported previously.

Insight into the role of bound water of a nucleating agent in polymer nucleation: a comparative study of anhydrous and monohydrated orotic acid on crystallization of poly(l-lactic acid)

The bound water of orotic acid and its dehydration transition play a negative role in nucleation effects on PLLA crystallization.

Accurate and efficient representation of intramolecular energy in ab initio generation of crystal structures. I. Adaptive local approximate models

The global search stage of crystal structure prediction (CSP) methods requires a fine balance between accuracy and computational cost, particularly for the study of large flexible molecules. A major improvement in the accuracy and cost of the intramolecular energy function used in the CrystalPredictor II [Habgood et al. (2015). J. Chem. Theory Comput. 11, 1957-1969] program is presented, where the most efficient use of computational effort is ensured via the use of adaptive local approximate model (LAM) placement. The entire search space of the relevant molecule's conformations is initially evaluated using a coarse, low accuracy grid. Additional LAM points are then placed at appropriate points determined via an automated process, aiming to minimize the computational effort expended in high-energy regions whilst maximizing the accuracy in low-energy regions. As the size, complexity and flexibility of molecules increase, the reduction in computational cost becomes marked. This improvement is illustrated with energy calculations for benzoic acid and the ROY molecule, and a CSP study of molecule (XXVI) from the sixth blind test [Reilly et al. (2016). Acta Cryst. B72, 439-459], which is challenging due to its size and flexibility. Its known experimental form is successfully predicted as the global minimum. The computational cost of the study is tractable without the need to make unphysical simplifying assumptions.

Why do Hydrates (Solvates) Form in Small Neutral Organic Molecules? Exploring the Crystal Form Landscapes of the Alkaloids Brucine and Strychnine

Computational methods were used to generate and explore the crystal structure landscapes of the two alkaloids strychnine and brucine. The computed structures were analyzed and rationalized by correlating the modelling results to a rich pool of available experimental data. Despite their structural similarity, the two compounds show marked differences in the formation of solid forms. For strychnine only one anhydrous form is reported in the literature and two new solvates from 1,4-dioxane were detected in the course of this work. In contrast, 22 solid forms are so far known to exist for brucine, comprising two anhydrates, four hydrates (HyA - HyC and a 5.25-hydrate), twelve solvates (alcohols and acetone) and four heterosolvates (mixed solvates with water and alcohols). For strychnine it is hard to produce any solid form other than the stable anhydrate while the formation of specific solid state forms of brucine is governed by a complex interplay between temperature and relative humidity/water activity and it is rather a challenging to avoid hydrate formation. Differences in crystal packing and the high tendency for brucine to form hydrates are not intuitive from the molecular structure alone, as both molecules have hydrogen bond acceptor groups but lack hydrogen bond donor groups. Only the evaluation of the crystal energy landscapes, in particular the close-packed crystal structures and high-energy open frameworks containing voids of molecular (water) dimensions, allowed us to unravel the diverse solid state behavior of the two alkaloids at a molecular level. In this study we demonstrate that expanding the analysis of anhydrate crystal energy landscapes to higher energy structures and calculating the solvent-accessible volume can be used to estimate non-stoichiometric or channel hydrate (solvate) formation, without explicitly computing the hydrate/solvate crystal energy landscapes.

Can computed crystal energy landscapes help understand pharmaceutical solids?

Computational crystal structure prediction (CSP) methods can now be applied to the smaller pharmaceutical molecules currently in drug development. We review the recent uses of computed crystal energy landscapes for pharmaceuticals, concentrating on examples where they have been used in collaboration with industrial-style experimental solid form screening. There is a strong complementarity in aiding experiment to find and characterise practically important solid forms and understanding the nature of the solid form landscape.

Computational and Experimental Characterization of Five Crystal Forms of Thymine: Packing Polymorphism, Polytypism/Disorder and Stoichiometric 0.8-Hydrate

New polymorphs of thymine emerged in an experimental search for solid forms, which was guided by the computationally generated crystal energy landscape. Three of the four anhydrates (AH) are homeoenergetic (A° - C) and their packing modes differ only in the location of oxygen and hydrogen atoms. AHs A° and B are ordered phases, whereas AH C shows disorder (X-ray diffuse scattering). Anhydrates AHs A° and B are ordered phases, whereas AH C shows disorder (X-ray diffuse scattering). Analysis of the crystal energy landscape for alternative AH C hydrogen bonded ribbon motifs identified a number of different packing modes, whose 3D structures were calculated to deviate by less than 0.24 kJ mol-1 in lattice energy. These structures provide models for stacking faults. The three anhydrates A° - C show strong similarity in their powder X-ray diffraction, thermoanalytical and spectroscopic (IR and Raman) characteristics. The already known anhydrate AH A° was identified as the thermodynamically most stable form at ambient conditions; AH B and AH C are metastable but show high kinetic stability. The hydrate of thymine is stable only at water activities (aw) > 0.95 at temperatures ≤ 25 °C. It was found to be a stoichiometric hydrate despite being a channel hydrate with an unusual water:thymine ratio of 0.8:1. Depending on the dehydration conditions, either AH C or AH D is obtained. The hydrate is the only known precursor to AH D. This study highlights the value and complementarity of simultaneous explorations of computationally and experimentally generated solid form landscapes of a small molecule anhydrate ↔ hydrate system.