Surface Analysis-From Crystal Structures to Particle Properties

Understanding the surface properties of particles is crucial for optimizing the performance of formulated products in various industries. However, acquiring this understanding often requires expensive trial-and-error studies. Here, we present advanced surface analysis tools that enable the visualization and quantification of chemical and topological information derived from crystallographic data. By employing functional group analysis, roughness calculations, and statistical interaction data, we facilitate direct comparisons of surfaces. We further demonstrate the practicality of our approach by correlating the sticking propensity of distinct ibuprofen morphologies with surface and particle descriptors calculated from a single crystal structure. Our findings support and expand upon previous work, demonstrating that the presence of a carboxylic acid group on the {011} facet leads to significant differences in particle properties and explains the higher electrostatic potential observed in the block-like morphology. While our surface analysis tools are not intended to replace the importance of chemical intuition and expertise, they provide valuable insights for formulators and particle engineers, facilitating informed, data-driven decisions to mitigate formulation risks. This research represents a significant step toward a comprehensive understanding of particle surfaces and their impact on products.

Predicting the Strength of Cohesive and Adhesive Interparticle Interactions for Dry Powder Inhalation Blends of Terbutaline Sulfate with α-Lactose Monohydrate

Grid-based systematic search methods are used to investigate molecule-molecule, molecule-surface, and surface-surface contributions to interparticle interactions in order to identify the crystal faces that most strongly affect particle behavior during powder blend formulation and delivery processes. The model system comprises terbutaline sulfate (TBS) as an active pharmaceutical ingredient (API) and α-form lactose monohydrate (LMH). A combination of systematic molecular modeling and X-ray computed tomography (XCT) is used to determine not only the adhesive and cohesive interparticle energies but, also the agglomeration behavior during manufacturing and de-agglomeration behavior during delivery after inhalation. This is achieved through a detailed examination of the balance between the adhesive and cohesive energies with the XCT results confirming the blend segregation tendencies, through the particle-particle de-agglomeration process. The results reveal that the cohesive interaction energies of TBS-TBS are higher than the adhesive energies between TBS and LMH, but that the cohesive energies of LMH-LMH are the smallest between molecule and molecule, molecule and surface, and surface and surface. This shows how systematic grid-search molecular modeling along with XCT can guide the digital formulation design of inhalation powders in order to achieve optimum aerosolization and efficacy for inhaled medicines. This will lead to faster pharmaceutical design with less variability, higher quality, and enhanced performance.

Exploring the CSD Drug Subset: An Analysis of Lattice Energies and Constituent Intermolecular Interactions for the Crystal Structures of Pharmaceuticals

Intermolecular (synthonic) modelling is used for a statistical analysis of crystal lattice energies, together with their contributing intermolecular interactions for the crystallographic structures selected from the CCDC's Drug Subset (https://doi.org/10.1016/j.xphs.2018.12.011). Analysis of this selected subset reveal similarities in packing compared to other organic crystals in the CSD with linear relationships between molecular weight and unit cell volume, void space, and packing coefficient. Crystal lattice energy calculations converge within a 30 Å intermolecular radius characterised by a mean lattice energy of ca. -36 kcal mol^-1 with ca. 85% and 15% due to dispersive and electrostatic interactions, respectively. The distribution of the strongest synthons within the individual structures reveals an average strength of -5.79 kcal mol^-1. The diversity of chemical space within the drug molecules is in agreement with the analysis of atom types across the selected subset with phenyl groups being found to contribute the highest mean energy of -11.28 kcal mol^-1, highlighting the importance of aromatic interactions within pharmaceutical compounds. Despite an initial focus on Z' = 1 structures, this automated approach enables rapid and consistent quantitative analysis of lattice energy, synthon strength and functional group contributions, providing solid-form informatics for pharmaceutical R&D and a helpful basis for further investigations.

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 has two well-characterised polymorphic forms and, whilst both crystallise in the monoclinic space group, they have quite different crystal chemistry and crystallisability behaviour. The molecular conformation...

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.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-021-03048-2.

Crystallographic tomography and molecular modelling of structured organic polycrystalline powders

Novel combination of crystallographic tomography and molecular modelling is used to examine the powder packing behaviour and crystal interactions for an organic polycrystalline powder bed.