Probing phase transitions in simvastatin with terahertz time-domain spectroscopy

Progress in application of terahertz time-domain spectroscopy for pharmaceutical analyses

Terahertz time-domain spectroscopy is an analytical method using terahertz time-domain pulses to study the physical and chemical properties of substances. It has strong potential for application in pharmaceutical analyses as an original non-destructive, efficient and convenient technology for spectral detection. This review briefly introduces the working principle of terahertz time-domain spectroscopy technology, focuses on the research achievements of this technology in analyses of chemical drugs, traditional Chinese medicine and biological drugs in the past decade. We also reveal the scientific feasibility of practical application of terahertz time-domain spectroscopy for pharmaceutical detection. Finally, we discuss the problems in practical application of terahertz time-domain spectroscopy technology, and the prospect of further development of this technology in pharmaceutical analyses. We hope that this review can provide a reference for application of terahertz time-domain spectroscopy technology in pharmaceutical analyses in the future.

Prediction and Preparation of Coamorphous Phases of a Bislactam

The effectiveness of a partial least squares-discriminant analysis coamorphous prediction model was tested using coamorphous screening data for a promising coamorphous former, the dimer of N-vinyl(caprolactam) (bisVCap) with a range of active pharmaceutical ingredients. The prediction model predicted 71% of the systems correctly. An experimental coamorphous screen was performed with this coformer with 13 different active pharmaceutical ingredients, and the results were compared to the predictions from the model. A total of 85% of the systems were correctly predicted. Stability assessments of three coamorphous systems showed that the prediction model score did not strongly correlate with the stability of the coamorphous material. The model performed well with small-molecule coformers, such as bisVCap, despite the difference in structure and properties compared to the amino-acid-based model training set.

Selected Applications of Terahertz Pulses in Medicine and Industry

This article contains a brief summary of areas where terahertz technology is making an impact in research and industrial applications. We cover some of its uses in the pharmaceutical setting, where both imaging and spectroscopy play important roles. Medical applications are also being pursued in many research laboratories, primarily for imaging purposes and following on from the first results just over 20 years ago. The three-dimensional imaging capability of pulsed terahertz allows for the observation of tumours below the surface of tissue, such as basal cell carcinoma of skin. The recent use of the technology in studies of cultural heritage has shown to increase our understanding of the past. The power of terahertz is exemplified by the discussion on its importance in different industries, such as semiconductor circuit manufacturing and automotive assembly.

New Insights into the Crystallographic Disorder in the Polymorphic Forms of Aspirin from Low-Frequency Vibrational Analysis

Terahertz time-domain spectroscopy (THz-TDS) is applied to two polymorphs of acetylsalicylic acid (aspirin), and the experimental spectra are compared to lattice dynamical calculations using high accuracy density functional theory. The calculations confirm that forms I and II have very close energetic and thermodynamic properties and also that they show similar spectral features in the far-infrared region, reflecting the high degree of similarity in their crystal structures. Unique vibrational modes are identified for each polymorph which allow them to be distinguished using THz-TDS measurements. The observation of spectral features attributable to both polymorphic forms in a single sample, however, provides further evidence to support the hypothesis that crystalline aspirin typically comprises intergrown domains of forms I and II. Differences observed in the baseline of the measured THz-TDS spectra indicate a greater degree of structural disorder in the samples of form II. Calculated Gibbs free-energy curves show a turning point at 75 K, inferring that form II is expected to be more stable than form I above this temperature as a result of its greater vibrational entropy. The calculations do not account for any differences in configurational entropy that may arise from expected structural defects. Further computational work on these structures, such as ab initio molecular dynamics, would be very useful to further explore this perspective. Here, aspirin is a model system to show how the additional insight from the low-frequency vibrational information complements the structural data and allows for quantitative thermodynamic information of pharmaceutical polymorphs to be extracted. The methodology is directly applicable to other polymorphic systems.

A terahertz time-domain super-resolution imaging method using a local-pixel graph neural network for biological products

The low image acquisition speed of terahertz (THz) time-domain imaging systems limits their application in biological products analysis. In the current study, a local pixel graph neural network was built for THz time-domain imaging super-resolution. The method could be applied to the analysis of any heterogeneous biological products as it only required a small number of sample images for training and particularly it focused on THz feature frequencies. The graph network applied the Fourier transform to graphs extracted from low-resolution (LR) images bringing an invariance of rotation and flip for local pixels, and the network then learnt the relationship between the state of graphs and the corresponding pixels to be reconstructed. With wood cores and seeds as examples, the images of these samples were captured by a THz time-domain imaging system for training and analysed by the method, achieving the root mean square error (RMSE) of pixels of 0.0957 and 0.1061 for the wood core and seed images, respectively. In addition, the reconstructed high-resolution (HR) images, LR images and true HR images at several feature frequencies were also compared in the current study. Results indicated that the method could not only reconstruct the spatial details and the useful signals from high noise signals at high feature frequencies but could also operate super-resolution in both spatial and spectral aspects.

The necessity of periodic boundary conditions for the accurate calculation of crystalline terahertz spectra

Terahertz vibrational spectroscopy has emerged as a powerful spectroscopic technique, providing valuable information regarding long-range interactions - and associated collective dynamics - occurring in solids. However, the terahertz sciences are relatively nascent, and there have been significant advances over the last several decades that have profoundly influenced the interpretation and assignment of experimental terahertz spectra. Specifically, because there do not exist any functional group or material-specific terahertz transitions, it is not possible to interpret experimental spectra without additional analysis, specifically, computational simulations. Over the years simulations utilizing periodic boundary conditions have proven to be most successful for reproducing experimental terahertz dynamics, due to the ability of the calculations to accurately take long-range forces into account. On the other hand, there are numerous reports in the literature that utilize gas phase cluster geometries, to varying levels of apparent success. This perspective will provide a concise introduction into the terahertz sciences, specifically terahertz spectroscopy, followed by an evaluation of gas phase and periodic simulations for the assignment of crystalline terahertz spectra, highlighting potential pitfalls and good practice for future endeavors.

Probing lattice vibration of alkali halide crystals by broadband terahertz spectroscopy

Terahertz spectral features of alkali halide crystals were studied with the combination of broadband terahertz time-domain spectroscopy and the solid-state-based density functional theory calculations. To understand the particular modes of the observed terahertz features of the alkali halide crystals, the resonant modes of KCl and CsCl were analyzed using face-centered cubic and body-centered cubic lattice models, respectively. The results show that the characteristic terahertz absorption peaks could be assigned to the lattice vibration of the ionic crystals. Furthermore, the terahertz responses of a series of alkali halides were recorded, and obvious absorption peaks were observed in each salt in the frequency region below 8.5 THz. What is more interestingly is that the frequencies of these observed peaks are red-shifted with the increases of the mass and radius of the ions. This correlation between the resonant frequency of the lattice vibration, the reduced atomic mass, and the equilibrium distance between the ions agrees well with the harmonic oscillator model.

Advanced therapeutic inhalation aerosols of a Nrf2 activator and RhoA/Rho kinase (ROCK) inhibitor for targeted pulmonary drug delivery in pulmonary hypertension: design, characterization, aerosolization, in vitro 2D/3D human lung cell cultures, and in vivo efficacy

Inhalable nanostructured microparticles of simvastatin, a Nrf2 activator and RhoA/Rho kinase (ROCK) inhibitor, were rationally designed for targeted pulmonary delivery as dry powder inhalers (DPIs) for the treatment of pulmonary hypertension (PH). Advanced particle engineering design technology was employed to develop inhalable dry powders using different dilute feed concentrations and spray drying pump rates. Several analytical techniques were used comprehensively to characterize the physicochemical properties of the resulting powders. Scanning electron microscopy (SEM) was used to visualize particle morphology (shape), surface structure, size, and size distribution. Karl Fischer titration (KFT) was employed to quantify the residual water content in the powders. X-ray powder diffraction (XRPD) was used to determine crystallinity. Hot-stage microscopy (HSM) under cross-polarizing lens was used to observe the presence or absence of birefringence characteristic of crystallinity. Differential scanning calorimetry (DSC) was employed to quantify thermotropic phase behavior. Attenuated total reflectance (ATR)-Fourier-transform infrared (FTIR) spectroscopy and Raman spectroscopy were used to determine the molecular fingerprint of simvastatin powders before and after particle engineering design. In vitro aerosol dispersion performance was performed with three different Food and Drug Administration (FDA)-approved human DPI devices. Cell viability and transepithelial electrical resistance (TEER) were demonstrated using different in vitro human pulmonary cell two and three-dimensional models at the air-liquid interface, and in vivo safety in healthy rats by inhalation. Efficacy was demonstrated in the in vivo lamb model of PH. Four different inhalable powders of simvastatin were successfully produced. They possessed nanostructured surfaces and were in the inhalable size range. Simvastatin retained its crystallinity following particle engineering design. The more dilute feed concentration spray dried at the lower pump rate produced the smallest particles. All powders successfully aerosolized with all three DPI human devices. Inhaled simvastatin as an aerosol restored the endothelial function in the shunt lamb model of PH, as demonstrated by the reduction of pulmonary vascular resistance (PVR) in response to the endothelium-dependent vasodilator acetylcholine.The reviews of this paper are available via the supplemental material section.

Polymorphism in Simvastatin: Twinning, Disorder, and Enantiotropic Phase Transitions

Simvastatin is one of the most widely used active pharmaceutical ingredients for the treatment of hyperlipidemias. Because the compound is employed as a solid in drug formulations, particular attention should be given to the characterization of different polymorphs, their stability domains, and the nature of the phase transitions that relate them. In this work, the phase transitions delimiting the stability domains of three previously reported simvastatin forms were investigated from structural, energetics, and dynamical points of view based on single crystal X-ray diffraction (SCXRD), hot stage microscopy (HSM), and differential scanning calorimetry (DSC) experiments (conventional scans and heat capacity measurements), complemented with molecular dynamics (MD) simulations. Previous assignments of the crystal forms were confirmed by SCXRD: forms I and II were found to be orthorhombic ( P2₁2₁2₁, Z'/ Z = 1/4) and form III was monoclinic ( P2₁, Z'/ Z = 2/4). The obtained results further indicated that (i) the transitions between different forms are observed at 235.9 ± 0.1 K (form III → form II) and at 275.2 ± 0.2 K (form II → form I) in DSC runs carried out at 10 K min^-1 and close to these values when other types of techniques are used (e.g., HSM). (ii) They are enantiotropic (i.e., there is a transition temperature relating the two phases before fusion at which the stability order is reversed), fast, reversible, with very little hysteresis between heating and cooling modes, and occur under single crystal to single crystal conditions. (iii) A nucleation and growth mechanism seems to be followed since HSM experiments on single crystals evidenced the propagation of an interface, accompanied by a change of birefringence and crystal contraction or expansion (more subtle in the case of form III → form II), when the phase transitions are triggered. (iv) Consistent with the reversible and small hysteresis nature of the phase transitions, the SCXRD results indicated that the molecular packing is very similar in all forms and the main structural differences are associated with conformational changes of the "ester tail". (v) The MD simulations further suggested that the tail is essentially "frozen" in two conformations below the III → II transition temperature, becomes progressively less hindered throughout the stability domain of form II, and acquires a large conformational freedom above the II → I transition. Finally, the fact that these transitions were found to be fast and reversible suggests that polymorphism is unlikely to be a problem for pharmaceutical formulations employing crystalline simvastatin because, if present, the III and II forms will readily convert to form I at ambient temperature.

Pharmaceutical cocrystals, salts and polymorphs: Advanced characterization techniques

The main goal of a novel drug development is to obtain it with optimal physiochemical, pharmaceutical and biological properties. Pharmaceutical companies and scientists modify active pharmaceutical ingredients (APIs), which often are cocrystals, salts or carefully selected polymorphs, to improve the properties of a parent drug. To find the best form of a drug, various advanced characterization methods should be used. In this review, we have described such analytical methods, dedicated to solid drug forms. Thus, diffraction, spectroscopic, thermal and also pharmaceutical characterization methods are discussed. They all are necessary to study a solid API in its intrinsic complexity from bulk down to the molecular level, gain information on its structure, properties, purity and possible transformations, and make the characterization efficient, comprehensive and complete. Furthermore, these methods can be used to monitor and investigate physical processes, involved in the drug development, in situ and in real time. The main aim of this paper is to gather information on the current advancements in the analytical methods and highlight their pharmaceutical relevance.

Direct measurement of molecular mobility and crystallisation of amorphous pharmaceuticals using terahertz spectroscopy

Despite much effort in the area, no comprehensive understanding of the formation and behaviour of amorphous solids has yet been achieved. This severely limits the industrial application of such materials, including drug delivery where, in principle, amorphous solids have demonstrated their great usefulness in increasing the bioavailability of poorly aqueous soluble active pharmaceutical ingredients. Terahertz time-domain spectroscopy is a relatively novel analytical technique that can be used to measure the fast molecular dynamics of molecules with high accuracy in a non-contact and non-destructive fashion. Over the past decade a number of applications for the characterisation of amorphous drug molecules and formulations have been developed and it has been demonstrated how this technique can be used to determine the onset and strength in molecular mobility that underpins the crystallisation of amorphous drugs. In this review we provide an overview of the history, fundamentals and future perspective of pharmaceutical applications related to the terahertz dynamics of amorphous systems.

Classification of THz pulse signals using two-dimensional cross-correlation feature extraction and non-linear classifiers

This work provides a performance comparison of four different machine learning classifiers: multinomial logistic regression with ridge estimators (MLR) classifier, k-nearest neighbours (KNN), support vector machine (SVM) and naïve Bayes (NB) as applied to terahertz (THz) transient time domain sequences associated with pixelated images of different powder samples. The six substances considered, although have similar optical properties, their complex insertion loss at the THz part of the spectrum is significantly different because of differences in both their frequency dependent THz extinction coefficient as well as differences in their refractive index and scattering properties. As scattering can be unquantifiable in many spectroscopic experiments, classification solely on differences in complex insertion loss can be inconclusive. The problem is addressed using two-dimensional (2-D) cross-correlations between background and sample interferograms, these ensure good noise suppression of the datasets and provide a range of statistical features that are subsequently used as inputs to the above classifiers. A cross-validation procedure is adopted to assess the performance of the classifiers. Firstly the measurements related to samples that had thicknesses of 2mm were classified, then samples at thicknesses of 4mm, and after that 3mm were classified and the success rate and consistency of each classifier was recorded. In addition, mixtures having thicknesses of 2 and 4mm as well as mixtures of 2, 3 and 4mm were presented simultaneously to all classifiers. This approach provided further cross-validation of the classification consistency of each algorithm. The results confirm the superiority in classification accuracy and robustness of the MLR (least accuracy 88.24%) and KNN (least accuracy 90.19%) algorithms which consistently outperformed the SVM (least accuracy 74.51%) and NB (least accuracy 56.86%) classifiers for the same number of feature vectors across all studies. The work establishes a general methodology for assessing the performance of other hyperspectral dataset classifiers on the basis of 2-D cross-correlations in far-infrared spectroscopy or other parts of the electromagnetic spectrum. It also advances the wider proliferation of automated THz imaging systems across new application areas e.g., biomedical imaging, industrial processing and quality control where interpretation of hyperspectral images is still under development.

Characteristics of Low-Frequency Molecular Phonon Modes Studied by THz Spectroscopy and Solid-State ab Initio Theory: Polymorphs I and III of Diflunisal

THz absorption spectra of two polymorphs of diflunisal, form I and form III, exhibit distinct features due to the influence of packing conformations on the frequency distributions and IR activities of gamma point phonon modes within the 100 cm(-1) region. In order to understand the origins of these THz modes, we perform a detailed mode analysis. The result shows that although the spectral features are different, these low-frequency phonon modes of the two molecular polymorphs have similar vibrational characteristics in terms of harmonic couplings of intermolecular and intramolecular vibrations.

Investigation of the terahertz vibrational modes of ZIF-8 and ZIF-90 with terahertz time-domain spectroscopy

We present experimental and computational evidence that gate-opening modes for zeolitic imidazole frameworks can be observed at terahertz frequencies.