In Silico Prediction of Growth and Dissolution Rates for Organic Molecular Crystals: A Multiscale Approach

Spiral growth of multicomponent crystals: theoretical aspects

This paper presents recent advances in the theory of multicomponent crystal growth from gas or solution, focusing on the most common step-flow mechanisms: Burton-Cabrera-Frank, Chernov, and Gilmer-Ghez-Cabrera. Analytical expressions for the spiral crystal growth rate are presented, taking into account the properties of all species involved in the growth process. The paper also outlines theoretical approaches to consider these mechanisms in multicomponent systems, providing a foundation for future developments and exploration of previously unexplored effects. Some special cases are discussed, including the formation of nanoislands of pure components on the surface and their self-organization, the impact of applied mechanical stress on the growth rate, and the mechanisms of its influence on growth kinetics. The growth due to chemical reactions on the surface is also considered. Possible future directions for developing the theory are outlined. A brief overview of numerical approaches and software codes that are useful in theoretical studies of crystal growth is also given.

A Digital Mechanistic Workflow for Predicting Solvent-Mediated Crystal Morphology: The α and β Forms of l-Glutamic Acid

The solvent-mediated crystal morphologies of the α and β polymorphic forms of l-glutamic acid are presented. This work applies a digital mechanistically based workflow that encompasses calculation of the crystal lattice energy and its constituent intermolecular synthons, their interaction energies, and their key role in understanding and predicting crystal morphology as well as assessing the surface chemistry, topology, and solvent binding on crystal habit growth surfaces. Through a comparison between the contrasting morphologies of the conformational polymorphs of l-glutamic acid, this approach highlights how the interfacial chemistry of organic crystalline materials and their inherent anisotropic interactions with their solvation environments direct their crystal habit with potential impact on their further downstream processing behavior.

Interferometric Probing of Physical and Chemical Properties of Solutions: Noncontact Investigation of Liquids

Interferometry is a highly versatile tool for probing physical and chemical phenomena. In addition to the benefit of noncontact investigations, even spatially resolved information can be obtained by choosing a suitable setup. This review presents the evolution of the various setups that have evolved since the first interferometers were developed in the mid-nineteenth century and highlights the benefits, limitations, and typical areas of application. This review focuses on interferometry based on electromagnetic waves in the near-infrared and visible range applied to liquid samples, categorizes the chemical/physical properties (e.g., pressure, temperature, composition) and phenomena (e.g., evaporation, crystal growth, diffusion, thermophoresis) that can be assessed, and presents a comprehensive literature review of specific existing applications. Finally, it discusses some fundamental open questions with respect to geometric considerations and overlapping effects. Expected final online publication date for the Annual Review of Chemical and Biomolecular Engineering, Volume 13 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Intermolecular Forces Driving Hexamethylenetetramine Co-Crystal Formation, a DFT and XRD Analysis

Interest in co-crystals formation has been constantly growing since their discovery, almost a century ago. Such success is due to the ability to tune the physical-chemical properties of the components in solid state by avoiding a change in their molecular structure. The properties influenced by the co-crystals formation range from an improvement of mechanical features and chemical stability to different solubility. In the scientific research area, the pharmacological field is undoubtedly one of those in which an expansion of the co-crystal knowledge can offer wide benefits. In this work, we described the crystalline structure of hexamethylenetetramine co-crystallized with the isophthalic acid, and we compared it with another co-crystal, showing the same components but different stoichiometry. To give a wider overview on the nature of the interactions behind the observed crystal packing and to rationalize the reasons of its formation, a computational analysis on such structures was carried out.

Identifying Conformational Isomers of Organic Molecules in Solution via Unsupervised Clustering

We present a systematic approach for the identification of statistically relevant conformational macrostates of organic molecules from molecular dynamics trajectories. The approach applies to molecules characterized by an arbitrary number of torsional degrees of freedom and enables the transferability of the macrostates definition across different environments. We formulate a dissimilarity measure between molecular configurations that incorporates information on the characteristic energetic cost associated with transitions along all relevant torsional degrees of freedom. Such metric is employed to perform unsupervised clustering of molecular configurations based on the Fast Search and Find of Density Peaks algorithm. We apply this method to investigate the equilibrium conformational ensemble of Sildenafil, a conformationally complex pharmaceutical compound, in different environments including the crystal bulk, the gas phase, and three different solvents (acetonitrile, 1-butanol, and toluene). We demonstrate that while Sildenafil can adopt more than 100 metastable conformational configurations, only 12 are significantly populated across all of the environments investigated. Despite the complexity of the conformational space, we find that the most abundant conformers in solution are the closest to the conformers found in the most common Sildenafil crystal phase.

Capturing Crystal Shape Evolution from Molecular Simulations

A simple and efficient algorithm for tracking shape evolution of small-molecule organic crystals during molecular simulations is described. It is based on the reconstruction of a crystal surface from molecular coordinates using an alpha-shape triangulation algorithm followed by the DBSCAN clustering of neighboring triangles with similar normal vectors to crystal faces. No information except the unit cell parameters is needed beforehand, enabling the user to automatically detect not only existing but also new forming crystal faces and edges, which is valuable for prediction of growth and dissolution kinetics. The results are demonstrated for aspirin and paracetamol crystals.

First-Principles and Empirical Approaches to Predicting In Vitro Dissolution for Pharmaceutical Formulation and Process Development and for Product Release Testing

This manuscript represents the perspective of the Dissolution Working Group of the International Consortium for Innovation and Quality in Pharmaceutical Development (IQ) and of two focus groups of the American Association of Pharmaceutical Scientists (AAPS): Process Analytical Technology (PAT) and In Vitro Release and Dissolution Testing (IVRDT). The intent of this manuscript is to show recent progress in the field of in vitro predictive dissolution modeling and to provide recommended general approaches to developing in vitro predictive dissolution models for both early- and late-stage formulation/process development and batch release. Different modeling approaches should be used at different stages of drug development based on product and process understanding available at those stages. Two industry case studies of current approaches used for modeling tablet dissolution are presented. These include examples of predictive model use for product development within the space explored during formulation and process optimization, as well as of dissolution models as surrogate tests in a regulatory filing. A review of an industry example of developing a dissolution model for real-time release testing (RTRt) and of academic case studies of enabling dissolution RTRt by near-infrared spectroscopy (NIRS) is also provided. These demonstrate multiple approaches for developing data-rich empirical models in the context of science- and risk-based process development to predict in vitro dissolution. Recommendations of modeling best practices are made, focused primarily on immediate-release (IR) oral delivery products for new drug applications. A general roadmap is presented for implementation of dissolution modeling for enhanced product understanding, robust control strategy, batch release testing, and flexibility toward post-approval changes.

Polarized fluorescence of a crystal having uniaxially oriented molecules by a carbazole-diyl-bridged macrocage

An axially oriented π-electron system is achieved in a single crystal of a macrocage molecule, and polarized fluorescence of the single crystal was observed.

Dynamics and Thermodynamics of Ibuprofen Conformational Isomerism at the Crystal/Solution Interface

Conformational flexibility of molecules involved in crystal growth and dissolution is rarely investigated in detail and usually considered to be negligible in the formulation of mesoscopic models of crystal growth. In this work, we set out to investigate the conformational isomerism of ibuprofen as it approaches and is incorporated in the morphologically dominant {100} crystal face, in a range of different solvents: water, 1-butanol, toluene, cyclohexanone, cyclohexane, acetonitrile, and trichloromethane. To this end, we combine extensive molecular dynamics and well-tempered metadynamics simulations to estimate the equilibrium distribution of conformers, compute conformer-conformer transition rates, and extract the characteristic relaxation time of the conformer population in solution, adsorbed at the solid/liquid interface, incorporated in the crystal in contact with the mother solution, and in the crystal bulk. We find that, while the conformational equilibrium distribution is weakly dependent on the solvent, relaxation times are instead significantly affected by it. Furthermore, differences in the relaxation dynamics are enhanced on the crystal surface, where conformational transitions become slower and specific conformational transition pathways are hindered. This leads us to observe that the dominant mechanisms of conformational transition can also change significantly moving from the bulk solution to the crystal interface, even for a small molecule with limited conformational flexibility such as ibuprofen. Our findings suggest that understanding conformational flexibility is key to provide an accurate description of the solid/liquid interface during crystal dissolution and growth, and therefore, its relevance should be systematically assessed in the formulation of mesoscopic growth models.

Computer Simulations: Essential Tools for Crystal Growth Studies

This special issue discusses recent advances in computer simulation studies of crystal growth. Crystal growth is a key to innovation in science and technology. Owing to recent progress in computer performance, computer simulation studies of crystal growth have become increasingly important. This special issue covers a variety of simulation methods, including the Monte Carlo, molecular dynamics, first-principles, multiscale, and continuum simulation methods, which are used for studies on the fundamentals and applications of crystal growth and related phenomena for different materials, such as hard-sphere systems, ice, organic crystals, semiconductors, and graphene.