Molecular Dynamics of the Supercooled Pharmaceutical Agent Posaconazole Studied via Differential Scanning Calorimetry and Dielectric and Mechanical Spectroscopies

Water-mediated phase transformations of posaconazole: An intricate jungle of crystal forms

Posaconazole is a broad-spectrum antifungal agent exhibiting rich polymorphism. Up to now, a total of fourteen different crystal forms have been reported, sometimes with an ambiguous nomenclature, but less is known about their properties and stability relationships. Investigating the solid-state of a drug compound is essential to identify the most stable form under working conditions and to prevent the risk of undesired solid-phase transformations under processing and storage. In this paper, we study posaconazole polymorphism by providing a description of its polymorphs, hydrates, and solvates. Powder X-ray diffraction (PXRD), dynamic vapor sorption (DVS), spectroscopic and thermal techniques were employed to characterize the different forms. In addition, the solid-phase transformations of posaconazole in aqueous suspensions were studied by means of Raman microscopy. Surprisingly, we found that Form S, the crystal form contained in the marketed oral suspension, is not the most stable form in water. Form S readily converts to a more stable hydrate, i.e. Form A, after storage in water for two weeks. In the commercial oral formulation the conversion between the two forms is prevented by the presence of polysorbate 80. Such insights into the stabilizing excipient effects beyond particle dispersion are critical to formulators.

Exploring the Cocrystal Landscape of Posaconazole by Combining High-Throughput Screening Experimentation with Computational Chemistry

The development of multicomponent crystal forms, such as cocrystals, represents a means to enhance the dissolution and absorption properties of poorly water-soluble drug compounds. However, the successful discovery of new pharmaceutical cocrystals remains a time- and resource-consuming process. This study proposes the use of a combined computational-experimental high-throughput approach as a tool to accelerate and improve the efficiency of cocrystal screening exemplified by posaconazole. First, we employed the COSMOquick software to preselect and rank cocrystal candidates (coformers). Second, high-throughput crystallization experiments (HTCS) were conducted on the selected coformers. The HTCS results were successfully reproduced by liquid-assisted grinding and reaction crystallization, ultimately leading to the synthesis of thirteen new posaconazole cocrystals (7 anhydrous, 5 hydrates, and 1 solvate). The posaconazole cocrystals were characterized by PXRD, ¹H NMR, Fourier transform-Raman, thermogravimetry-Fourier transform infrared spectroscopy, and differential scanning calorimetry. In addition, the prediction performance of COSMOquick was compared to that of two alternative knowledge-based methods: molecular complementarity (MC) and hydrogen bond propensity (HBP). Although HBP does not perform better than random guessing for this case study, both MC and COSMOquick show good discriminatory ability, suggesting their use as a potential virtual tool to improve cocrystal screening.

Solubilizing temperature of crystalline drug in polymer carrier: A rheological investigation on a posaconazole-copovidone system with low drug load

Measuring the solubility of a crystalline active pharmaceutical ingredient (API) in a polymer-rich system is challenging due to the high viscosity of the polymer which kinetically impedes reaching the solubility equilibrium. In this study, a rheological approach of determining the solubilizing temperature of a crystalline API in a polymeric carrier has been developed. To validate the method, a model physical mixture of crystalline posaconazole and copovidone with a relatively low API load (25 wt%) was utilized. First, a comparison between conventional differential scanning calorimetry (DSC) and a rheological temperature ramp was conducted to illustrate that the rheological method could capture the melting point depression behavior similarly to the more well-known DSC technique. Second, to further understand the dissolution process of the crystalline posaconazole into the copovidone carrier and precisely measure the solubilizing temperature, a series of isothermal rheological time sweeps were carried out at various temperatures selected based on the rheological temperature sweep. Because the dissolved API molecule imparted a plasticizing effect to the polymeric carrier, the complex viscosity of the API-polymer system decreased gradually over time and correlated well to an exponential decay function. Moreover, dependent on the applied temperature, the API-polymer system eventually accomplished distinct equilibrium states (complex viscosities) within different time frames. The obtained time constants at different temperatures were fitted to the Arrhenius equation, allowing the determination of the activation energy of the mixing process. The results indicated that once the processing temperature exceeded a critical point below the melting point of the crystalline API, the API-polymer solubilization process switched from a surface dominated dispersive mechanism to a molecular-level solubilization mode, manifested by the significantly increased activation energy. To the best of our knowledge, the currently developed rheological approach was the first successful measurement of the solubilizing temperature of a crystalline drug in a polymer-rich system.

The impact of chemical structure on the formation of the medium-range order and dynamical properties of selected antifungal APIs

In this paper, we have analyzed structural, thermal, and dynamical properties of four azole antifungals: itraconazole (ITZ), posaconazole (POS), terconazole (TER) and ketoconazole (KET), differing mainly in the length of the rod-like backbone and slightly in side groups. Our investigations clearly demonstrated that the changes in the chemical structure result in a different ability to form the medium-range order (MRO) and variation in thermal and dynamical properties of these pharmaceuticals. Direct comparison of the diffractograms collected for glassy and crystalline materials indicated that the MRO observed in the former phases is related to maintaining the local molecular arrangement of the crystal structure. Moreover, it was shown that once the MRO-related diffraction peaks appear, additional mobility (δ- or α' relaxation), slower than the structural (α)-process, is also detected in dielectric spectra. This new mode is connected to the motions within supramolecular nanoaggregates. Detailed analysis of dielectric and calorimetric data also revealed that the variation in the internal structure and MRO of the examined pharmaceuticals have an impact on the glass transition temperature (Tg) shape of the α-process, isobaric fragility, molecular dynamics in the glassy state and number of dynamically correlated molecules. These findings could be helpful in an understanding the influence of different types of intermolecular MRO on the properties of substances having a similar chemical backbone.

Surface diffusion in glasses of rod-like molecules posaconazole and itraconazole: effect of interfacial molecular alignment and bulk penetration

The method of surface grating decay has been used to measure surface diffusion in the glasses of two rod-like molecules posaconazole (POS) and itraconazole (ITZ). Although structurally similar antifungal medicines, ITZ forms liquid-crystalline phases while POS does not. Surface diffusion in these systems is significantly slower than in the glasses of quasi-spherical molecules of similar volume when compared at the glass transition temperature T_g. Between the two systems, ITZ has slower surface diffusion. These results are explained on the basis of the near-vertical orientation of the rod-like molecules at the surface and their deep penetration into the bulk where mobility is low. For molecular glasses without extensive hydrogen bonds, we find that the surface diffusion coefficient at T_g decreases smoothly with the penetration depth of surface molecules and the trend has the double-exponential form for the surface mobility gradient observed in simulations. This supports the view that these molecular glasses have a similar mobility vs. depth profile and their different surface diffusion rates arise simply from the different depths at which molecules are anchored. Our results also provide support for a previously observed correlation between the rate of surface diffusion and the fragility of the bulk liquid.

Probing the Molecular-Level Interactions in an Active Pharmaceutical Ingredient (API) - Polymer Dispersion and the Resulting Impact on Drug Product Formulation

PURPOSE

An investigation of underlying mechanisms of API-polymer interaction patterns has the potential to provide valuable insights for selecting appropriate formulations with superior physical stability and processability.

MATERIALS AND METHODS

In this study, copovidone was used as a polymeric carrier for several model compounds including clotrimazole, nifedipine, and posaconazole. The varied chemical structures conferred the ability for the model compounds to form distinct interactions with copovidone. Rheology and nuclear magnetic resonance (NMR) were combined to investigate the molecular pattern and relative strength of active pharmaceutical ingredient (API)-polymer interactions. In addition, the impact of the interactions on formulation processability via hot melt extrusion (HME) and physical stability were evaluated.

RESULTS

The rheological response of an API-polymer system was found to be highly sensitive to API-polymer interaction, depending both on API chemistry and API-polymer miscibility. In the systems studied, dispersed API induced a stronger plasticizer effect on the polymer matrix compared to crystalline/aggregated API. Correspondingly, the processing torque via HME showed a proportional relationship with the maximum complex viscosity of the API-polymer system. In order to quantitatively evaluate the relative strength of the API-polymer interaction, homogeneously dispersed API-polymer amorphous samples were prepared by HME at an elevated temperature. DSC, XRD, and rheology were employed to confirm the amorphous integrity and homogeneity of the resultant extrudates. Subsequently, the homogeneously dispersed API-polymer amorphous dispersions were interrogated by rheology and NMR to provide a qualitative and quantitative assessment of the nature of the API-polymer interaction, both macroscopically and microscopically. Rheological master curves of frequency sweeps of the extrudates exhibited a strong dependence on the API chemistry and revealed a rank ordering of the relative strength of API-copovidone interactions, in the order of posaconazole > nifedipine > clotrimazole. NMR data provided the means to precisely map the API-polymer interaction pattern and identify the specific sites of interaction from a molecular perspective. Finally, the impact of API-polymer interactions on the physical stability of the resultant extrudates was studied.

CONCLUSION

Qualitative and quantitative evaluation of the relative strength of the API-polymer interaction was successfully accomplished by utilizing combined rheology and NMR. Graphical Abstract.

Vapor deposition of a nonmesogen prepares highly structured organic glasses

We show that glasses with aligned smectic liquid crystal-like order can be produced by physical vapor deposition of a molecule without any equilibrium liquid crystal phases. Smectic-like order in vapor-deposited films was characterized by wide-angle X-ray scattering. A surface equilibration mechanism predicts the highly smectic-like vapor-deposited structure to be a result of significant vertical anchoring at the surface of the equilibrium liquid, and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy orientation analysis confirms this prediction. Understanding of the mechanism enables informed engineering of different levels of smectic order in vapor-deposited glasses to suit various applications. The preparation of a glass with orientational and translational order from a nonliquid crystal opens up an exciting paradigm for accessing extreme anisotropy in glassy solids.

Evaluation and optimized selection of supersaturating drug delivery systems of posaconazole (BCS class 2b) in the gastrointestinal simulator (GIS): An in vitro-in silico-in vivo approach

Supersaturating drug delivery systems (SDDS) have been put forward in the recent decades in order to circumvent the issue of low aqueous solubility. Prior to the start of clinical trials, these enabling formulations should be adequately explored in in vitro/in silico studies in order to understand their in vivo performance and to select the most appropriate and effective formulation in terms of oral bioavailability and therapeutic outcome. The purpose of this work was to evaluate the in vivo performance of four different oral formulations of posaconazole (categorized as a biopharmaceutics classification system (BCS) class 2b compound) based on the in vitro concentrations in the gastrointestinal simulator (GIS), coupled with an in silico pharmacokinetic model to predict their systemic profiles. Recently published intraluminal and systemic concentrations of posaconazole for these formulations served as a reference to validate the in vitro and in silico results. Additionally, the morphology of the formed precipitate of posaconazole was visualized and characterized by optical microscopy studies and thermal analysis. This multidisciplinary work demonstrates an in vitro-in silico-in vivo approach that provides a scientific basis for screening SDDS by a user-friendly formulation predictive dissolution (fPD) device in order to rank these formulations towards their in vivo performance.

Molecular dynamics, physical and thermal stability of neat amorphous amlodipine besylate and in binary mixture

In this paper, a stable amorphous solid dispersion of an antihypertensive drug, amlodipine besylate (AMB) was prepared by entrapping it in a polymer matrix, polyvinyl pyrrollidone, in different weight ratios (AMB/PVP 05:95, 10:90, 20:80, 30:70). The glass forming ability of all binary dispersions were studied by means of differential scanning calorimetry and found good correlation between experimental T_g and Fox Flory's prediction. By considering the daily dosage limit of 5 mg, a weight ratio of 05:95 was further considered for the study. The structures of neat and binary of AMB were characterized by density functional theory, Fourier transform infrared spectroscopy, Fourier transform Raman spectroscopy and UV-visible spectroscopy. Further, detailed molecular dynamics of both pure and binary were investigated using broadband dielectric spectroscopy to judge the physical stability of the amorphous dispersions. Translation-rotation coupling of AMB possibly explains the dual conductivity and dipolar nature of the secondary relaxation in neat AMB. Thus, the binary dispersion of AMB with commercially acceptable weight ratio with strong glass forming behaviour and better shelf life was prepared and characterized for practical applications.

Interplay between the static ordering and dynamical heterogeneities determining the dynamics of rotation and ordinary liquid phases in 1,6-anhydro-β-D-glucose

In this letter, we reported thorough the structural and molecular dynamics studies on 1,6-anhydro-β-D-glucose, the second compound reported so far that is capable to form rotator and supercooled liquid phases. In contrast to the data presented for ethanol, temperature dependences of structural dynamics in both phases are very comparable. On the other hand, X ray measurements revealed unusually long range ordering/correlations between molecules in the ODIC (d ≈ 95 Å) and supercooled phases (d ≈ 30-40 Å) of this carbohydrate. Our consideration clearly demonstrated that the interplay between length scales of static range ordering and dynamical heterogeneities as well as internal molecular arrangement seem to be the key to understanding the molecular dynamics of different materials characterized by varying degree of disorder in the vicinity of the glass transition temperature.

Molecular Relaxations in Supercooled Liquid and Glassy States of Amorphous Quinidine: Dielectric Spectroscopy and Density Functional Theory Approaches

In this article, we conduct a comprehensive molecular relaxation study of amorphous Quinidine above and below the glass-transition temperature (Tg) through broadband dielectric relaxation spectroscopy (BDS) experiments and theoretical density functional theory (DFT) calculations, as one major issue with the amorphous state of pharmaceuticals is life expectancy. These techniques enabled us to determine what kind of molecular motions are responsible, or not, for the devitrification of Quinidine. Parameters describing the complex molecular dynamics of amorphous Quinidine, such as Tg, the width of the α relaxation (βKWW), the temperature dependence of α-relaxation times (τα), the fragility index (m), and the apparent activation energy of secondary γ relaxation (Ea-γ), were characterized. Above Tg (> 60 °C), a medium degree of nonexponentiality (βKWW = 0.5) was evidenced. An intermediate value of the fragility index (m = 86) enabled us to consider Quinidine as a glass former of medium fragility. Below Tg (< 60 °C), one well-defined secondary γ relaxation, with an apparent activation energy of Ea-γ = 53.8 kJ/mol, was reported. From theoretical DFT calculations, we identified the most reactive part of Quinidine moieties through exploration of the potential energy surface. We evidenced that the clearly visible γ process has an intramolecular origin coming from the rotation of the CH(OH)C9H14N end group. An excess wing observed in amorphous Quinidine was found to be an unresolved Johari-Goldstein relaxation. These studies were supplemented by sub-Tg experimental evaluations of the life expectancy of amorphous Quinidine by X-ray powder diffraction and differential scanning calorimetry. We show that the difference between Tg and the onset temperature for crystallization, Tc, which is 30 K, is sufficiently large to avoid recrystallization of amorphous Quinidine during 16 months of storage under ambient conditions.

Universal Behavior of Dielectric Responses of Glass Formers: Role of Dipole-Dipole Interactions

From an exhaustive examination of the molecular dynamics in practically all van der Waals molecular glass formers ever probed by dielectric spectroscopy, we found that the width of the α-loss peak at or near the glass transition temperature T_{g} is strongly anticorrelated with the polarity of the molecule. The larger the dielectric relaxation strength Δε(T_{g}) of the system, the narrower is the α-loss peak. This remarkable property is explained by the contribution from the dipole-dipole interaction potential V_{dd}(r)=-Dr^{-6} to the attractive part of the intermolecular potential, making the resultant potential more harmonic, and the effect increases rapidly with the dipole moment μ and Δε(T_{g}) in view of the relation, D∝(μ^{4}/kT_{g})∝kT_{g}[Δε(T_{g})]^{2}. Since the novel correlation discovered encompasses practically all van der Waals molecular glass formers studied by dielectric spectroscopy, it impacts the large dielectric research community as well as those engaged in solving the glass transition problem.