Influence of Solvent Selection on the Crystallizability and Polymorphic Selectivity Associated with the Formation of the "Disappeared" Form I Polymorph of Ritonavir.

The comparative crystallizability and polymorphic selectivity of ritonavir, a novel protease inhibitor for the treatment of acquired immune-deficiency syndrome, as a function of solvent selection are examined through an integrated and self-consistent experimental and computational molecular modeling study. Recrystallization at high supersaturation by rapid cooling at 283.15 K is found to produce the metastable "disappeared" polymorphic form I from acetone, ethyl acetate, acetonitrile, and toluene solutions in contrast to ethanol which produces the stable form II. Concomitant crystallization of the other known solid forms is not found under these conditions. Isothermal crystallization studies using turbidometric detection based upon classical nucleation theory reveal that, for an equal induction time, the required driving force needed to initiate solution nucleation decreases with solubility in the order of ethanol, acetone, acetonitrile, ethyl acetate, and toluene consistent with the expected desolvation behavior predicted from the calculated solute solvation free energies. Molecular dynamics simulations of the molecular and intermolecular chemistry reveal the presence of conformational interplay between intramolecular and intermolecular interactions within the solution phase. These encompass the solvent-dependent formation of intramolecular O-H...O hydrogen bonding between the hydroxyl and carbamate groups coupled with differing conformations of the hydroxyl's shielding phenyl groups. These conformational preferences and their relative interaction propensities, as a function of solvent selection, may play a rate-limiting role in the crystallization behavior by not only inhibiting to different degrees the nucleation process but also restricting the assembly of the optimal intermolecular hydrogen bonding network needed for the formation of the stable form II polymorph.

The structural pathway from its solvated molecular state to the solution crystallisation of the α- and ß-polymorphic forms of para amino benzoic acid.

Para amino benzoic acid (PABA) has two well-characterised α- and ß-polymorphic forms and, whilst both crystallise in the monoclinic space group P21/n, they have quite different crystal chemistry and crystallisability behaviour. Previous work has shown that the molecular conformation deformation energy in the crystalline state is higher for the ß-form than for the α-form and that the lattice energy for the former converges more slowly than for the latter overall. This suggests that not only is there a higher barrier to crystallisation for the ß-form but also that low solution supersaturations might be needed for it to preferentially nucleate. Additionally, solute cluster propensity and solute solvation energetic analysis highlight the importance of an aqueous solvation environment in inhibiting the α-form's strong OH⋯O carboxylic acid hydrogen bond (H-bond) dimer. Despite this, the detailed molecular-scale pathway from solvated molecules to 3D crystallographic structure still remains unclear, most notably regarding how the nucleation process is activated and how, as a result, this mediates the preferential formation of either of the two polymorphic forms. Molecular dynamics (MD) simulations coupled with FTIR studies and intermolecular synthon analysis address this issue through characterisation of the propensity of the incipient bulk synthons that are important in the crystallisation of the two polymorphic forms within the solution state. MD molecular trajectory analysis within crystallisation solutions reveals a greater propensity for OH⋯O synthons (both single H-bonds and homodimers) typical of the α-form and NH⋯O synthons found in both the α- and ß-forms when compared to aqueous solution but much lower propensities for the ß-form's "fingerprinting" OH⋯N and π-π stacking synthons. In contrast, data from the aqueous solution environment reveals a much greater propensity for the ß-form's π-π interaction synthons. IR dilution studies in acetonitrile in the carbonyl region reveal the presence of two CO vibrational stretching bands, whose relative intensities vary as a function of solution dilution. These were assigned to the solvated PABA monomer and a COOH dimer of PABA. Similar data in ethanol shows a main CO stretching band with a shoulder peak suggesting a similar monomer vs. dimer speciation may exist in this solvent. The IR data is consistent with the organic solvent MD data, albeit the corresponding analysis for the aqueous solution was precluded due to the latter's strong OH vibrational mode which restricted validation in aqueous solutions.

Molecular, Solid-State and Surface Structures of the Conformational Polymorphic Forms of Ritonavir in Relation to their Physicochemical Properties.

PURPOSE: Application of multi-scale modelling workflows to characterise polymorphism in ritonavir with regard to its stability, bioavailability and processing. METHODS: Molecular conformation, polarizability and stability are examined using quantum mechanics (QM). Intermolecular synthons, hydrogen bonding, crystal morphology and surface chemistry are modelled using empirical force fields. RESULTS: The form I conformation is more stable and polarized with more efficient intermolecular packing, lower void space and higher density, however its shielded hydroxyl is only a hydrogen bond donor. In contrast, the hydroxyl in the more open but less stable and polarized form II conformation is both a donor and acceptor resulting in stronger hydrogen bonding and a more stable crystal structure but one that is less dense. Both forms have strong 1D networks of hydrogen bonds and the differences in packing energies are partially offset in form II by its conformational deformation energy difference with respect to form I. The lattice energies converge at shorter distances for form I, consistent with its preferential crystallization at high supersaturation. Both forms exhibit a needle/lath-like crystal habit with slower growing hydrophobic side and faster growing hydrophilic capping habit faces with aspect ratios increasing from polar-protic, polar-aprotic and non-polar solvents, respectively. Surface energies are higher for form II than form I and increase with solvent polarity. The higher deformation, lattice and surface energies of form II are consistent with its lower solubility and hence bioavailability. CONCLUSION: Inter-relationship between molecular, solid-state and surface structures of the polymorphic forms of ritonavir are quantified in relation to their physical-chemical properties.

Understanding the Effects of a Polymer on the Surface Dissolution of Pharmaceutical Cocrystals Using Combined Experimental and Molecular Dynamics Simulation Approaches.

The molecular interactions between the surfaces of cocrystals [i.e., flufenamic acid and theophylline (FFA-TP), flufenamic acid and nicotinamide (FFA-NIC), and carbamazepine and nicotinamide (CBZ-NIC)] and the polymers [i.e., polyethylene glycol (PEG), polyvinylpyrrolidone (PVP), and copolymer of vinylpyrrolidone (60%)/vinyl acetate (40%) (PVP-VA)] were investigated through combined experimental and molecular dynamics simulation approaches to resolve the mechanisms of cocrystal dissolution and precipitation. It was found that adsorption of the polymers on the surfaces of cocrystals might prevent the precipitation of the parent drug and alter the dissolution rate. The effect of polymers on precipitation could be determined by the cocrystal dissolution rate, the interactions of polymers with the surfaces of cocrystals, the characters of the noncovalent bonds formed between the polymers and the cocrystal surfaces, and the mobility and conformation of the polymers. The etching experiments of single cocrystals revealed that FFA-NIC and CBZ-NIC appeared as surface precipitation cocrystals while FFA-TP could lead to bulk precipitation. Both PVP and PVP-VA were good precipitation inhibitors for FFA-NIC, and they could completely inhibit the recrystallization of FFA III on the surfaces of dissolving cocrystals. In addition, as the adsorption of the polymer was slower than dissolution rate of the cocrystals, PVP and PVP-VA could only partially inhibit the recrystallization of CBZ dihydrate on the surface of CBZ-NIC. While PEG had no inhibitory effect on the surface crystallization of FFA-NIC and CBZ-NIC, due to its weak interactions with the surfaces of the cocrystals, it enhanced the dissolution performance of FFA-TP. In contrast, PVP and PVP-VA reduced the dissolution rate of FFA-TP and subsequently undermined the performance of cocrystals. Taken together, the approach of combining experimental and molecular dynamics simulation provided insights into the mechanisms of cocrystal dissolution as well as the polymers acting as inhibitory excipients for precipitation/recrystallization, making contribution to the development of novel formulations.

Off-the-shelf DFT-DISPersion methods: Are they now "on-trend" for organic molecular crystals?

Organic molecular crystals contain long-range dispersion interactions that can be challenging for solid-state methods such as density functional theory (DFT) to capture, and in some industrial sectors are overlooked in favor of classical methods to calculate atomistic properties. Hence, this publication addresses the critical question of whether dispersion corrected DFT calculations for organic crystals can reproduce the structural and energetic trends seen from experiment, i.e., whether the calculations can now be said to be truly "on-trend." In this work, we assess the performance of three of the latest dispersion-corrected DFT methods, in calculating the long-range, dispersion energy: the pairwise methods of D3(0) and D3(BJ) and the many-body dispersion method, MBD@rsSCS. We calculate the energetics and optimized structures of two homologous series of organic molecular crystals, namely, carboxylic acids and amino acids. We also use a classical force field method (using COMPASS II) and compare all results to experimental data where possible. The mean absolute error in lattice energies is 9.59 and 343.85 kJ/mol (COMPASS II), 10.17 and 16.23 kJ/mol (MBD@rsSCS), 10.57 and 18.76 kJ/mol [D3(0)], and 8.52 and 14.66 kJ/mol [D3(BJ)] for the carboxylic acids and amino acids, respectively. MBD@rsSCS produces structural and energetic trends that most closely match experimental trends, performing the most consistently across the two series and competing favorably with COMPASS II.

Crystal regeneration - a unique growth phenomenon observed in organic crystals post breakage.

Crystal regeneration has been observed in macroscopic paracetamol crystals post breakage along their cleavage plane. High resolution imaging confirmed regeneration rates to be 3-fold faster than growth prior to breakage. Further analysis of the solute-solvent interactions is required to elucidate the process which currently lacks linearity with traditional growth theories.

Navigating the Complex Solid Form Landscape of the Quercetin Flavonoid Molecule.

Quercetin, a naturally occurring bioflavonoid substance widely used in the nutraceutical and food industries, exists in various solid forms that can have different physicochemical properties, thus impacting this compound's performance in various applications. In this work, we will clarify the complex solid-form landscape of this molecule. Two elusive isostructural solvates of quercetin were obtained from ethanol and methanol. The obtained crystals were characterized experimentally, but the crystallographic structure could not be solved due to their high instability. Nevertheless, the desolvated structure resulting from a high-temperature treatment (or prolonged storage at ambient conditions) of both these two labile crystals was characterized and solved via powder X-ray diffraction and solid-state nuclear magnetic resonance (SSNMR). This anhydrous crystal structure was compared with another anhydrous quercetin form obtained in our previous work, indicating that, at least, two different anhydrous polymorphs of quercetin exist. Navigating the solid-form landscape of quercetin is essential to ensure accurate control of the functional properties of food, nutraceutical, or pharmaceutical products containing crystal forms of this substance.

Relating Crystal Structure to Surface Properties: A Study on Quercetin Solid Forms.

The surface energy and surface chemistry of a crystal are of great importance when designing particles for a specific application, as these will impact both downstream manufacturing processes as well as final product quality. In this work, the surface properties of two different quercetin solvates (quercetin dihydrate and quercetin DMSO solvate) were studied using molecular (synthonic) modeling and experimental techniques, including inverse gas chromatography (IGC) and contact angle measurements, to establish a relationship between crystal structure and surface properties. The attachment energy model was used to predict morphologies and calculate surface properties through the study of their growth synthons. The modeling results confirmed the surface chemistry anisotropy for the two forms. For quercetin dihydrate, the {010} facets were found to grow mainly by nonpolar offset quercetin-quercetin stacking interactions, thus being hydrophobic, while the {100} facets were expected to be hydrophilic, growing by a polar quercetin-water hydrogen bond. For QDMSO, the dominant facet {002} grows by a strong polar quercetin-quercetin hydrogen bonding interaction, while the second most dominant facet {011} grows by nonpolar π-π stacking interactions. Water contact angle measurements and IGC confirmed a greater overall surface hydrophilicity for QDMSO compared to QDH and demonstrated surface energy heterogeneity for both structures. This work shows how synthonic modeling can help in the prediction of the surface nature of crystalline particles and guide the choice of parameters that will determine the optimal crystal form and final morphology for targeted surface properties, for example, the choice of crystallization conditions, choice of solvent, or presence of additives or impurities, which can direct the crystallization of a specific crystal form or crystal shape.

Triglycine (GGG) Adopts a Polyproline II (pPII) Conformation in Its Hydrated Crystal Form: Revealing the Role of Water in Peptide Crystallization.

Polyproline II (pPII) is a left-handed 31-helix conformation, which has been observed to be the most abundant secondary structure in unfolded peptides and proteins compared to α-helix and ß-sheet. Although pPII has been reported as the most stable conformation for several unfolded short chain peptides in aqueous solution, it is rarely observed in their solid state. Here, we show for the first time a glycine homopeptide (gly-gly-gly) adopting the pPII conformation in its crystalline dihydrate structure. The single crystal X-ray structure with molecular dynamic simulation suggests that a network of water and the charged carboxylate group is critical in stabilizing the pPII conformation in solid state, offering an insight into the structures of unfolded regions of proteins and the role of water in peptide crystallization.

Crystallographic Structure, Intermolecular Packing Energetics, Crystal Morphology and Surface Chemistry of Salmeterol Xinafoate (Form I).

Single crystals of salmeterol xinafoate (form I), prepared from slow cooled supersaturated propan-2-ol solutions, crystallize in a triclinic P1¯ symmetry with 2 closely related independent salt pairs within the asymmetric unit, with an approximately double-unit cell volume compared with the previously published crystal structure. Synthonic analysis of the bulk intermolecular packing confirms the similarity in packing energetics between the 2 salt pairs. The strongest synthons, as expected, are dominated by coulombic interactions. Morphologic prediction reveals a plate-like morphology, dominated by the {001}, {010}, and {100} surfaces, consistent with experimentally grown crystals. Although surface chemistry of the slow-growing {001} face comprises large sterically hindering phenyl groups, although weaker coulombic interactions still prevail from the alcohol group present on the phenyl and hydroxymethyl groups. The surface chemistry of the faster growing {010} and {100} faces are dominated by the significantly stronger cation/anion interactions occurring between the carboxylate and protonated secondary ammonium ion groups. The importance of understanding the cohesive and adhesive nature of the crystal surfaces of an active pharmaceutical ingredient, with respect to their interaction with other active pharmaceutical ingredient crystals and how that may affect formulation design, is highlighted.